520 Astronomie und zugeordnete Wissenschaften
Refine
Year of publication
Document Type
- Article (315)
- Doctoral Thesis (87)
- Conference Proceeding (66)
- Postprint (23)
- Other (7)
- Review (3)
- Habilitation Thesis (2)
- Master's Thesis (1)
Keywords
- stars: massive (18)
- cosmic rays (17)
- methods: numerical (17)
- stars: atmospheres (15)
- stars: evolution (15)
- stars: mass-loss (14)
- stars: winds, outflows (13)
- stars: early-type (12)
- MHD (11)
- hydrodynamics (11)
Institute
Relativistic pair beams produced in the cosmic voids by TeV gamma rays from blazars are expected to produce a detectable GeV-scale cascade emission missing in the observations. The suppression of this secondary cascade implies either the deflection of the pair beam by intergalactic magnetic fields (IGMFs) or an energy loss of the beam due to the electrostatic beam-plasma instability. IGMF of femto-Gauss strength is sufficient to significantly deflect the pair beams reducing the flux of secondary cascade below the observational limits. A similar flux reduction may result in the absence of the IGMF from the beam energy loss by the instability before the inverse Compton cooling. This dissertation consists of two studies about the instability role in the evolution of blazar-induced beams.
Firstly, we investigated the effect of sub-fG level IGMF on the beam energy loss by the instability. Considering IGMF with correlation lengths smaller than a few kpc, we found that such fields increase the transverse momentum of the pair beam particles, dramatically reducing the linear growth rate of the electrostatic instability and hence the energy-loss rate of the pair beam. Our results show that the IGMF eliminates beam plasma instability as an effective energy-loss agent at a field strength three orders of magnitude below that needed to suppress the secondary cascade emission by magnetic deflection. For intermediate-strength IGMF, we do not know a viable process to explain the observed absence of GeV-scale cascade emission and hence can be excluded.
Secondly, we probed how the beam-plasma instability feeds back on the beam, using a realistic two-dimensional beam distribution. We found that the instability broadens the beam opening angles significantly without any significant energy loss, thus confirming a recent feedback study on a simplified one-dimensional beam distribution. However, narrowing diffusion feedback of the beam particles with Lorentz factors less than 1e6 might become relevant even though initially it is negligible. Finally, when considering the continuous creation of TeV pairs, we found that the beam distribution and the wave spectrum reach a new quasi-steady state, in which the scattering of beam particles persists and the beam opening angle may increase by a factor of hundreds. This new intrinsic scattering of the cascade can result in time delays of around ten years, thus potentially mimicking the IGMF deflection. Understanding the implications on the GeV cascade emission requires accounting for inverse Compton cooling and simulating the beam-plasma system at different points in the IGM.
Ammonium salts are a reservoir of nitrogen on a cometary nucleus and possibly on some asteroids
(2020)
The measured nitrogen-to-carbon ratio in comets is lower than for the Sun, a discrepancy which could be alleviated if there is an unknown reservoir of nitrogen in comets. The nucleus of comet 67P/Churyumov-Gerasimenko exhibits an unidentified broad spectral reflectance feature around 3.2 micrometers, which is ubiquitous across its surface. On the basis of laboratory experiments, we attribute this absorption band to ammonium salts mixed with dust on the surface. The depth of the band indicates that semivolatile ammonium salts are a substantial reservoir of nitrogen in the comet, potentially dominating over refractory organic matter and more volatile species. Similar absorption features appear in the spectra of some asteroids, implying a compositional link between asteroids, comets, and the parent interstellar cloud.
Ammonium salts are a reservoir of nitrogen on a cometary nucleus and possibly on some asteroids
(2020)
The measured nitrogen-to-carbon ratio in comets is lower than for the Sun, a discrepancy which could be alleviated if there is an unknown reservoir of nitrogen in comets. The nucleus of comet 67P/Churyumov-Gerasimenko exhibits an unidentified broad spectral reflectance feature around 3.2 micrometers, which is ubiquitous across its surface. On the basis of laboratory experiments, we attribute this absorption band to ammonium salts mixed with dust on the surface. The depth of the band indicates that semivolatile ammonium salts are a substantial reservoir of nitrogen in the comet, potentially dominating over refractory organic matter and more volatile species. Similar absorption features appear in the spectra of some asteroids, implying a compositional link between asteroids, comets, and the parent interstellar cloud.
Macro-economic assessments of climate impacts lack an analysis of the distribution of daily rainfall, which can resolve both complex societal impact channels and anthropogenically forced changes(1-6). Here, using a global panel of subnational economic output for 1,554 regions worldwide over the past 40 years, we show that economic growth rates are reduced by increases in the number of wet days and in extreme daily rainfall, in addition to responding nonlinearly to the total annual and to the standardized monthly deviations of rainfall. Furthermore, high-income nations and the services and manufacturing sectors are most strongly hindered by both measures of daily rainfall, complementing previous work that emphasized the beneficial effects of additional total annual rainfall in low-income, agriculturally dependent economies(4,7). By assessing the distribution of rainfall at multiple timescales and the effects on different sectors, we uncover channels through which climatic conditions can affect the economy. These results suggest that anthropogenic intensification of daily rainfall extremes(8-10) will have negative global economic consequences that require further assessment by those who wish to evaluate the costs of anthropogenic climate change.
Transparent conductive materials based on indium oxide remain yet irreplaceable in various optoelectronic applications. Amorphous oxides appear especially attractive for technology as they are isotropic, demonstrate relatively high electron mobility and can be processed at low temperatures. Among them is indium zinc oxide (IZO) with a large zinc content that is crucial for keeping the amorphous state but redundant for the doping. In this work we investigated water-free and water containing IZO films obtained by radio frequency sputtering. The correlation between temperature driven changes of the chemical state, the optical and electrical properties as well as the progression of crystallization was in focus. Such characterization methods as: scanning electron microscopy, X-ray diffraction, X-ray photoelectron spectroscopy, Raman spectroscopy, temperature dependent Hall-effect measurements and others were applied. Temperature dependent electrical properties of amorphous IZO and IZO:H2O films were found to evolve similarly. Based on our experience in In2O3:H2O (In2O3:H or IOH) we proposed an explanation for the changes observed. Water admixture was found to decrease crystallization temperature of IZO significantly from similar to 550 degrees C to similar to 280 degrees C. Herewith, the presence and concentration of water and/or hydroxyls was found to determine Zn distribution in the film. In particular, Zn enrichment was detected at the film's surface respective to the high water and/or hydroxyl amount. Raman spectra revealed a two-dimensional crystallization of w-ZnO which precedes regardless water presence an extensive In2O3 crystallization. An abrupt loss of electron mobility as a result of crystallization was attributed to the formation of ZnO interlayer on grain boundaries.
We investigate the initiation and early evolution of 12 solar eruptions, including six active-region hot channel and six quiescent filament eruptions, which were well observed by the Solar Dynamics Observatory, as well as by the Solar Terrestrial Relations Observatory for the latter. The sample includes one failed eruption and 11 coronal mass ejections, with velocities ranging from 493 to 2140 km s(-1). A detailed analysis of the eruption kinematics yields the following main results. (1) The early evolution of all events consists of a slow-rise phase followed by a main-acceleration phase, the height-time profiles of which differ markedly and can be best fit, respectively, by a linear and an exponential function. This indicates that different physical processes dominate in these phases, which is at variance with models that involve a single process. (2) The kinematic evolution of the eruptions tends to be synchronized with the flare light curve in both phases. The synchronization is often but not always close. A delayed onset of the impulsive flare phase is found in the majority of the filament eruptions (five out of six). This delay and its trend to be larger for slower eruptions favor ideal MHD instability models. (3) The average decay index at the onset heights of the main acceleration is close to the threshold of the torus instability for both groups of events (although, it is based on a tentative coronal field model for the hot channels), suggesting that this instability initiates and possibly drives the main acceleration.
We investigate the initiation and early evolution of 12 solar eruptions, including six active-region hot channel and six quiescent filament eruptions, which were well observed by the Solar Dynamics Observatory, as well as by the Solar Terrestrial Relations Observatory for the latter. The sample includes one failed eruption and 11 coronal mass ejections, with velocities ranging from 493 to 2140 km s(-1). A detailed analysis of the eruption kinematics yields the following main results. (1) The early evolution of all events consists of a slow-rise phase followed by a main-acceleration phase, the height-time profiles of which differ markedly and can be best fit, respectively, by a linear and an exponential function. This indicates that different physical processes dominate in these phases, which is at variance with models that involve a single process. (2) The kinematic evolution of the eruptions tends to be synchronized with the flare light curve in both phases. The synchronization is often but not always close. A delayed onset of the impulsive flare phase is found in the majority of the filament eruptions (five out of six). This delay and its trend to be larger for slower eruptions favor ideal MHD instability models. (3) The average decay index at the onset heights of the main acceleration is close to the threshold of the torus instability for both groups of events (although, it is based on a tentative coronal field model for the hot channels), suggesting that this instability initiates and possibly drives the main acceleration.
Diffraction enhanced imaging (DEI) is an advanced digital radiographic imaging technique employing the refraction of X-rays to contrast internal interfaces. This study aims to qualitatively and quantitatively evaluate images acquired using this technique and to assess how different fitting functions to the typical rocking curves (RCs) influence the quality of the images. RCs are obtained for every image pixel. This allows the separate determination of the absorption and the refraction properties of the material in a position-sensitive manner. Comparison of various types of fitting functions reveals that the Pseudo-Voigt (PsdV) function is best suited to fit typical RCs. A robust algorithm was developed in the Python programming language, which reliably extracts the physically meaningful information from each pixel of the image. We demonstrate the potential of the algorithm with two specimens: a silicone gel specimen that has well-defined interfaces, and an additively manufactured polycarbonate specimen.
Over the past 3 yr, the fading non-thermal emission from the GW170817 remained generally consistent with the afterglow powered by synchrotron radiation produced by the interaction of the structured jet with the ambient medium. Recent observations by Hajela et al. indicate the change in temporal and spectral behaviour in the X-ray band. We show that the new observations are compatible with the emergence of a new component due to non-thermal emission from the fast tail of the dynamical ejecta of ab-initio binary neutron star merger simulations. This provides a new avenue to constrain binary parameters. Specifically, we find that equal mass models with soft equations of state (EOSs) and high-mass ratio models with stiff EOSs are disfavoured as they typically predict afterglows that peak too early to explain the recent observations. Moderate stiffness and mass ratio models, instead, tend to be in good overall agreement with the data.
We present the analysis of the optical variability of the early, nitrogen-rich Wolf-Rayet (WR) star WR 7. The analysis of multisector Transiting Exoplanet Survey Satellite (TESS) light curves and high-resolution spectroscopic observations confirm multiperiodic variability that is modulated on time-scales of years. We detect a dominant period of 2.6433 +/- 0.0005 d in the TESS sectors 33 and 34 light curves in addition to the previously reported high-frequency features from sector 7. We discuss the plausible mechanisms that may be responsible for such variability in WR 7, including pulsations, binarity, co-rotating interaction regions (CIRs), and clumpy winds. Given the lack of strong evidence for the presence of a stellar or compact companion, we suggest that WR 7 may pulsate in quasi-coherent modes in addition to wind variability likely caused by CIRs on top of stochastic low-frequency variability. WR 7 is certainly a worthy target for future monitoring in both spectroscopy and photometry to sample both the short (less than or similar to 1 d) and long (greater than or similar to 1000 d) variability time-scales.
We here present the results from a detailed analysis of nebular abundances of commonly observed ions in the collisional ring galaxy Cartwheel using the Very Large Telescope (VLT) Multi-Unit Spectroscopic Explorer (MUSE) data set. The analysis includes 221 H II regions in the star-forming ring, in addition to 40 relatively fainter H a-emitting regions in the spokes, disc, and the inner ring. The ionic abundances of He, N, O, and Fe are obtained using the direct method (DM) for 9, 20, 20, and 17 ring H II regions, respectively, where the S++ temperature-sensitive line is detected. For the rest of the regions, including all the nebulae between the inner and the outer ring, we obtained O abundances using the strong-line method (SLM). The ring regions have a median 12 + log O/H = 8.19 +/- 0.15, log N/O = -1.57 +/- 0.09 and log Fe/O = -2.24 +/- 0.09 using the DM. Within the range of O abundances seen in the Cartwheel, the N/O and Fe/O values decrease proportionately with increasing O, suggesting local enrichment of O without corresponding enrichment of primary N and Fe. The O abundances of the disc H II regions obtained using the SLM show a well-defined radial gradient. The mean O abundance of the ring H II regions is lower by similar to 0.1 dex as compared to the extrapolation of the radial gradient. The observed trends suggest the preservation of the pre-collisional abundance gradient, displacement of most of the processed elements to the ring, as predicted by the recent simulation by Renaud et al., and post-collisional infall of metal-poor gas in the ring.
The increasing number of known exoplanets raises questions about their demographics and the mechanisms that shape planets into how we observe them today. Young planets in close-in orbits are exposed to harsh environments due to the host star being magnetically highly active, which results in high X-ray and extreme UV fluxes impinging on the planet. Prolonged exposure to this intense photoionizing radiation can cause planetary atmospheres to heat up, expand and escape into space via a hydrodynamic escape process known as photoevaporation. For super-Earth and sub-Neptune-type planets, this can even lead to the complete erosion of their primordial gaseous atmospheres. A factor of interest for this particular mass-loss process is the activity evolution of the host star. Stellar rotation, which drives the dynamo and with it the magnetic activity of a star, changes significantly over the stellar lifetime. This strongly affects the amount of high-energy radiation received by a planet as stars age. At a young age, planets still host warm and extended envelopes, making them particularly susceptible to atmospheric evaporation. Especially in the first gigayear, when X-ray and UV levels can be 100 - 10,000 times higher than for the present-day sun, the characteristics of the host star and the detailed evolution of its high-energy emission are of importance.
In this thesis, I study the impact of stellar activity evolution on the high-energy-induced atmospheric mass loss of young exoplanets. The PLATYPOS code was developed as part of this thesis to calculate photoevaporative mass-loss rates over time. The code, which couples parameterized planetary mass-radius relations with an analytical hydrodynamic escape model, was used, together with Chandra and eROSITA X-ray observations, to investigate the future mass loss of the two young multiplanet systems V1298 Tau and K2-198. Further, in a numerical ensemble study, the effect of a realistic spread of activity tracks on the small-planet radius gap was investigated for the first time. The works in this thesis show that for individual systems, in particular if planetary masses are unconstrained, the difference between a young host star following a low-activity track vs. a high-activity one can have major implications: the exact shape of the activity evolution can determine whether a planet can hold on to some of its atmosphere, or completely loses its envelope, leaving only the bare rocky core behind. For an ensemble of simulated planets, an observationally-motivated distribution of activity tracks does not substantially change the final radius distribution at ages of several gigayears. My simulations indicate that the overall shape and slope of the resulting small-planet radius gap is not significantly affected by the spread in stellar activity tracks. However, it can account for a certain scattering or fuzziness observed in and around the radius gap of the observed exoplanet population.
Synchronization regimes in an ensemble of phase oscillators coupled through a diffusion field
(2022)
We consider an ensemble of identical phase oscillators coupled through a common diffusion field. Using the Ott-Antonsen reduction, we develop dynamical equations for the complex local order parameter and the mean field. The regions of the existence and stability are determined for the totally synchronous, partially synchronous, and asynchronous spatially homogeneous states. A procedure of searching for inhomogeneous states as periodic trajectories of an auxiliary system of the ordinary differential equations is demonstrated. A scenario of emergence of chimera structures from homogeneous synchronous solutions is described.
We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh AREPO code. Therein, the anisotropic transport of cosmic rays (CRs) along magnetic field lines is discretized using a path-conservative finite volume method on the unstructured time-dependent Voronoi mesh of AREPO. The interaction of CRs and gyroresonant Alfven waves is described by short time-scale source terms in the CRHD equations. We employ a custom-made semi-implicit adaptive time stepping source term integrator to accurately integrate this interaction on the small light-crossing time of the anisotropic transport step. Both the transport and the source term integration step are separated from the evolution of the magnetohydrodynamical equations using an operator split approach. The new algorithm is tested with a variety of test problems, including shock tubes, a perpendicular magnetized discontinuity, the hydrodynamic response to a CR overpressure, CR acceleration of a warm cloud, and a CR blast wave, which demonstrate that the coupling between CR and magnetohydrodynamics is robust and accurate. We demonstrate the numerical convergence of the presented scheme using new linear and non-linear analytic solutions.
To undergo diffusive shock acceleration, electrons need to be preaccelerated to increase their energies by several orders of magnitude, else their gyroradii will be smaller than the finite width of the shock. In oblique shocks, where the upstream magnetic field orientation is neither parallel nor perpendicular to the shock normal, electrons can escape to the shock upstream, modifying the shock foot to a region called the electron foreshock. To determine the preacceleration in this region, we undertake particle-in-cell simulations of oblique shocks while varying the obliquity and in-plane angles. We show that while the proportion of reflected electrons is negligible for theta (Bn) = 74.degrees 3, it increases to R similar to 5% for theta (Bn) = 30 degrees, and that, via the electron acoustic instability, these electrons power electrostatic waves upstream with energy density proportional to R (0.6) and a wavelength approximate to 2 lambda (se), where lambda (se) is the electron skin length. While the initial reflection mechanism is typically a combination of shock-surfing acceleration and magnetic mirroring, we show that once the electrostatic waves have been generated upstream, they themselves can increase the momenta of upstream electrons parallel to the magnetic field. In less than or similar to 1% of cases, upstream electrons are prematurely turned away from the shock and never injected downstream. In contrast, a similar fraction is rescattered back toward the shock after reflection, reinteracts with the shock with energies much greater than thermal, and crosses into the downstream.
The strong chromospheric absorption lines Ca ii H & K are tightly connected to stellar surface magnetic fields. Only for the Sun, spectral activity indices can be related to evolving magnetic features on the solar disk. The Solar Disk-Integrated (SDI) telescope feeds the Potsdam Echelle Polarimetric and Spectroscopic Instrument (PEPSI) of the Large Binocular Telescope (LBT) at Mt. Graham International Observatory, Arizona, U.S.A. We present high-resolution, high-fidelity spectra that were recorded on 184 & 82 days in 2018 & 2019 and derive the Ca ii H & K emission ratio, that is, the S-index. In addition, we compile excess brightness and area indices based on full-disk Ca ii K-line-core filtergrams of the Chromospheric Telescope (ChroTel) at Observatorio del Teide, Tenerife, Spain and full-disk ultraviolet (UV) 1600 angstrom images of the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO). Thus, Sun-as-a-star spectral indices are related to their counterparts derived from resolved images of the solar chromosphere. All indices display signatures of rotational modulation, even during the very low magnetic activity in the minimum of Solar Cycle 24. Bringing together different types of activity indices has the potential to join disparate chromospheric datasets yielding a comprehensive description of chromospheric activity across many solar cycles.
Levy walks are continuous-time random-walk processes with a spatiotemporal coupling of jump lengths and waiting times. We here apply the Hermite polynomial method to study the behavior of LWs with power-law walking time density for four different cases. First we show that the known result for the infinite density of an unconfined, unbiased LW is consistently recovered. We then derive the asymptotic behavior of the probability density function (PDF) for LWs in a constant force field, and we obtain the corresponding qth-order moments. In a harmonic external potential we derive the relaxation dynamic of the LW. For the case of a Poissonian walking time an exponential relaxation behavior is shown to emerge. Conversely, a power-law decay is obtained when the mean walking time diverges. Finally, we consider the case of an unconfined, unbiased LW with decaying speed v(r ) = v0/./r. When the mean walking time is finite, a universal Gaussian law for the position-PDF of the walker is obtained explicitly.
Gas flows in galaxy mergers
(2022)
In major galaxy mergers, the orbits of stars are violently perturbed, and gas is torqued to the centre, diluting the gas metallicity and igniting a starburst. In this paper, we study the gas dynamics in and around merging galaxies using a series of cosmological magnetohydrodynamical zoom-in simulations. We find that the gas bridge connecting the merging galaxies pre-coalescence is dominated by turbulent pressure, with turbulent Mach numbers peaking at values of 1.6-3.3. This implies that bridges are dominated by supersonic turbulence, and are thus ideal candidates for studying the impact of extreme environments on star formation. We also find that gas accreted from the circumgalactic medium (CGM) during the merger significantly contributes (27-51 percent) to the star formation rate (SFR) at the time of coalescence and drives the subsequent reignition of star formation in the merger remnant. Indeed, 19-53 percent of the SFR at z = 0 originates from gas belonging to the CGM prior the merger. Finally, we investigate the origin of the metallicity-diluted gas at the centre of merging galaxies. We show that this gas is rapidly accreted on to the Galactic Centre with a time-scale much shorter than that of normal star-forming galaxies. This explains why coalescing galaxies are not well-captured by the fundamental metallicity relation.
Constraints on the time-scales of quasar activity are key to understanding the formation and growth of supermassive black holes (SMBHs), quasar triggering mechanisms, and possible feedback effects on their host galaxies. However, observational estimates of this so-called quasar lifetime are highly uncertain (t(Q) similar to 10(4)-10(9) yr), because most methods are indirect and involve many model-dependent assumptions. Direct evidence of earlier activity is gained from the higher ionization state of the intergalactic medium (IGM) in the quasar environs, observable as enhanced Ly alpha transmission in the so-called proximity zone. Due to the similar to 30 Myr equilibration time-scale of He II in the z similar to 3 IGM, the size of the He II proximity zone depends on the time the quasar had been active before our observation t(on) <= t(Q), enabling up to +/- 0.2 dex precise measurements of individual quasar on-times that are comparable to the e-folding time-scale t(S) <= 44 Myr of SMBH growth. Here we present the first statistical sample of 13 quasars whose accurate and precise systemic redshifts allow for measurements of sufficiently precise He II quasar proximity zone sizes between similar or equal to 2 and similar or equal to 15 proper Mpc from science-grade Hubble Space Telescope (HST) spectra. Comparing these sizes to predictions from cosmological hydrodynamical simulations post-processed with 1D radiative transfer, we infer a broad range of quasar on-times from t(on) less than or similar to 1Myr to t(on) > 30 Myr that does not depend on quasar luminosity, black hole mass, or Eddington ratio. These results point to episodic quasar activity over a long duty cycle, but do not rule out substantial SMBH growth during phases of radiative inefficiency or obscuration.
Owing to global warming and particularly high regional ocean warming, both Thwaites and Pine Island Glaciers in the Amundsen region of the Antarctic Ice Sheet could lose their buttressing ice shelves over time. We analyse the possible consequences using the parallel ice sheet model (PISM), applying a simple cliff-calving parameterization and an ice melange-buttressing model. We find that the instantaneous loss of ice-shelf buttressing, due to enforced ice-shelf melting, initiates grounding-line retreat and triggers marine ice sheet instability (MISI). As a consequence, the grounding line progresses into the interior of the West Antarctic Ice Sheet and leads to a sea level contribution of 0.6 m within 100 a. By subjecting the exposed ice cliffs to cliff calving using our simplified parameterization, we also analyse marine ice cliff instability (MICI). In our simulations it can double or even triple the sea level contribution depending on the only loosely constrained parameter that determines the maximum cliff-calving rate. The speed of MICI depends on this upper bound of the calving rate, which is given by the ice melange buttressing the glacier. However, stabilization of MICI may occur for geometric reasons. Because the embayment geometry changes as MICI advances into the interior of the ice sheet, the upper bound on calving rates is reduced and the progress of MICI is slowed down. Although we cannot claim that our simulations bear relevant quantitative estimates of the effect of ice-melange buttressing on MICI, the mechanism has the potential to stop the instability. Further research is needed to evaluate its role for the past and future evolution of the Antarctic Ice Sheet.