520 Astronomie und zugeordnete Wissenschaften
Refine
Year of publication
Document Type
- Article (352)
- Doctoral Thesis (88)
- Conference Proceeding (66)
- Postprint (23)
- Other (7)
- Review (3)
- Habilitation Thesis (2)
- Master's Thesis (1)
Keywords
- stars: massive (20)
- cosmic rays (17)
- methods: numerical (17)
- subdwarfs (17)
- stars: atmospheres (16)
- stars: evolution (16)
- stars: mass-loss (14)
- stars: winds, outflows (13)
- stars: early-type (12)
- MHD (11)
Magnetic fields can play an important role in stellar evolution. Among white dwarfs, the most common stellar remnant, the fraction of magnetic systems is more than 20 per cent. The origin of magnetic fields in white dwarfs, which show strengths ranging from 40 kG to hundreds of MG, is still a topic of debate. In contrast, only one magnetic hot subdwarf star has been identified out of thousands of known systems. Hot subdwarfs are formed from binary interaction, a process often associated with the generation of magnetic fields, and will evolve to become white dwarfs, which makes the lack of detected magnetic hot subdwarfs a puzzling phenomenon. Here we report the discovery of three new magnetic hot subdwarfs with field strengths in the range 300-500 kG. Like the only previously known system, they are all helium-rich O-type stars (He-sdOs). We analysed multiple archival spectra of the three systems and derived their stellar properties. We find that they all lack radial velocity variability, suggesting formation via a merger channel. However, we derive higher than typical hydrogen abundances for their spectral type, which are in disagreement with current model predictions. Our findings suggest a lower limit to the magnetic fraction of hot subdwarfs of 0.147(+0.143)(-0.047) per cent, and provide evidence for merger-induced magnetic fields which could explain white dwarfs with field strengths of 50-150 MG, assuming magnetic flux conservation.
Context.
Thermonuclear supernovae (SNe), a subset of which are the highly important SNe of Type Ia and Iax, are relatively poorly understood phenomena. One of the more promising scenarios leading up to the creation of a thermonuclear SN involves accretion of helium-rich material from a binary companion. Following the SN, the binary companion is then ejected from the location of the progenitor binary at velocities possibly large enough to unbind it from the gravitational potential of the Galaxy. Ejected companion stars should form a detectable population, if their production mechanism is not exceedingly rare.
Aims.
This study builds on previous works, producing the most extensive prediction of the properties of such a hypothetical population to date, taking both Chandrasekhar and non-Chandrasekhar mass events into account. These results are then used to define criteria for membership of this population and characterise putative subpopulations.
Methods.
This study contains 6 x 10(6) individual ejection trajectories out of the Galactic plane calculated with the stellar kinematics framework SHyRT, which are analysed with regard to their bulk observational properties. These are then put into context with the only previously identified population member US 708 and applied to a number of other possible candidate objects.
Results.
We find that two additional previously observed objects possess properties to warrant a designation as candidate objects. Characterisation of these object with respect to the predicted population finds all of them to be extreme in at least one astrometric observable. Higher mass ( >0 :7 M-circle dot) objects should be over-represented in the observationally accessible volume, with the ratio of bound to unbound objects being an accessible observable for the determination of the dominant terminal accretor mass. We find that current observations of runaway candidates within 10 kpc support a Galactic SN rate of the order of similar to 3 x 10(-7) yr(-1) to similar to 2 x 10(-6) yr(-1), three orders of magnitude below the inferred Galactic SN Ia rate and two orders of magnitude below the formation rate of predicted He-donor progenitors.
Conclusions.
The number of currently observed population members suggests that the He-donor scenario, as suspected before, is not a dominant contributor to the number of observed SNe Ia. However, even at the low event rate suggested, we find that the majority of possibly detectable population members is still undetected. The extreme nature of current population members suggests that a still larger number of objects has simply evaded detection up to this point, hinting at a higher contribution than is currently supported by observation.
Losses of radiation belt energetic particles by encounters with four of the inner Moons of Jupiter
(2022)
Based on an improved model of the moon absorption of Jovian radiation belt particles, we investigate quantitatively and comprehensively the absorption probabilities and particle lifetimes due to encounters with four of the inner moons of Jupiter (Amalthea, Thebe, Io, and Europa) inside L < 10. Our results demonstrate that the resultant average lifetimes of energetic protons and electrons vary dramatically between similar to 0.1 days and well above 1,000 days, showing a strong dependence on the particle equatorial pitch angle, kinetic energy and moon orbit. The average lifetimes of energetic protons and electrons against moon absorption are shortest for Io (i.e., similar to 0.1-10 days) and longest for Thebe (i.e., up to thousands of days), with the lifetimes in between for Europa and Amalthea. Due to the diploe tilt angle absorption effect, the average lifetimes of energetic protons and electrons vary markedly below and above alpha eq ${\alpha }_{\mathrm{e}\mathrm{q}}$ = 67 degrees. Overall, the average electron lifetimes exhibit weak pitch angle dependence, but the average proton lifetimes are strongly dependent on equatorial pitch angle. The average lifetimes of energetic protons decrease monotonically and substantially with the kinetic energy, but the average lifetimes of energetic electrons are roughly constant at energies <similar to 10 MeV, increase substantially around the Kepler velocities of the moons (similar to 10-50 MeV), and decrease quickly at even higher energies. Compared with the averaged electron lifetimes, the average proton lifetimes are longer at energies below a few MeV and shorter at energies above tens of MeV.
Context. The Large Magellanic Cloud (LMC) is the most luminous satellite galaxy of the Milky Way and, owing to its companion, the Small Magellanic Cloud (SMC), represents an excellent laboratory to study the interaction of dwarf galaxies. Aims. The aim of this study is to investigate the kinematics of the outer regions of the LMC by using stellar proper motions to understand the impact of interactions, for example with the SMC about 250 Myr ago. Methods. We calculate proper motions using multi-epoch K s -band images from the VISTA survey of the Magellanic Cloud system (VMC). Observations span a time baseline of 2-5 yr. We combine the VMC data with data from the Gaia Early Data Release 3 and introduce a new method to distinguish between Magellanic and Milky Way stars based on a machine learning algorithm. This new technique enables a larger and cleaner sample selection of fainter sources as it reaches below the red clump of the LMC. Results. We investigate the impact of the SMC on the rotational field of the LMC and find hints of stripped SMC debris. The southeastern region of the LMC shows a slow rotational speed compared to the overall rotation. N-body simulations suggest that this could be caused by a fraction of stripped SMC stars located in that particular region that move opposite to the expected rotation.
With the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure that uncertainties inherent in the gravitational-wave models will not lead to systematic biases when information from multiple detections is combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyze them with a set of four different waveform models. For our analysis with 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to 1250 m (2% at 90% credible interval) but that systematic differences between currently employed waveform models can be twice as large. Hence, it will be essential to ensure that systematic biases will not become dominant in inferences of the neutron-star equation of state when capitalizing on future developments.
Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudospectral code ELLIPTICA, an infrastructure for construction of initial data for various binary and single gravitational systems of all kinds. The elliptic equations under consideration are solved on a single spatial hypersurface of the spacetime manifold. Using coordinate maps, the hypersurface is covered by patches whose boundaries can adapt to the surface of the compact objects. To solve elliptic equations with arbitrary boundary condition, ELLIPTICA deploys a Schur complement domain decomposition method with a direct solver. In this version, we use cubed sphere coordinate maps and the fields are expanded using Chebyshev polynomials of the first kind. Here, we explain the building blocks of ELLIPTICA and the initial data construction algorithm for a black hole-neutron star binary system. We perform convergence tests and evolve the data to validate our results. Within our framework, the neutron star can reach spin values close to breakup with arbitrary direction, while the black hole can have arbitrary spin with dimensionless spin magnitude ∼0.8.
Recent observations have revealed remarkable insights into the gas reservoir in the circumgalactic medium (CGM) of galaxy haloes. In this paper, we characterise the gas in the vicinity of Milky Way and Andromeda analogues in the hestia (High resolution Environmental Simulations of The Immediate Area) suite of constrained Local Group (LG) simulations. The hestia suite comprise of a set of three high-resolution arepo-based simulations of the LG, run using the Auriga galaxy formation model.
For this paper, we focus only on the 𝑧 = 0 simulation datasets and generate mock skymaps along with a power spectrum analysis to show that the distributions of ions tracing low-temperature gas (H i and Si iii) are more clumpy in comparison to warmer gas tracers (O vi, O vii and O viii). We compare to the spectroscopic CGM observations of M31 and low-redshift galaxies.
hestia under-produces the column densities of the M31 observations, but the simulations are consistent with the observations of low-redshift galaxies. A possible explanation for these findings is that the spectroscopic observations of M31 are contaminated by gas residing in the CGM of the Milky Way.
Effect of temperature on the densification of silicate melts to lower earth's mantle conditions
(2022)
Physical properties of silicate melts play a key role for global planetary dynamics, controlling for example volcanic eruption styles, mantle convection and elemental cycling in the deep Earth. They are significantly modified by structural changes at the atomic scale due to external parameters such as pressure and temperature or due to chemistry. Structural rearrangements such as 4- to 6-fold coordination change of Si with increasing depth may profoundly influence melt properties, but have so far mostly been studied at ambient temperature due to experimental difficulties. In order to investigate the structural properties of silicate melts and their densification mechanisms at conditions relevant to the deep Earth's interior, we studied haplo basaltic glasses and melts (albite-diopside composition) at high pressure and temperature conditions in resistively and laser-heated diamond anvil cells using X-ray absorption near edge structure spectroscopy. Samples were doped with 10 wt% of Ge, which is accessible with this experimental technique and which commonly serves as a structural analogue for the network forming cation Si. We acquired spectra on the Ge K edge up to 48 GPa and 5000 K and derived the average Ge-O coordination number NGe-O, and bond distance RGe-O as functions of pressure. Our results demonstrate a continuous transformation from tetrahedral to octahedral coordination between ca. 5 and 30 GPa at ambient temperature. Above 1600 K the data reveal a reduction of the pressure needed to complete conversion to octahedral coordination by ca. 30 %. The results allow us to determine the influence of temperature on the Si coordination number changes in natural melts in the Earth's interior. We propose that the complete transition to octahedral coordination in basaltic melts is reached at about 40 GPa, corresponding to a depth of ca. 1200 km in the uppermost lower mantle. At the core-mantle boundary (2900 km, 130 GPa, 3000 K) the existence of non-buoyant melts has been proposed to explain observed low seismic wave velocity features. Our results highlight that the melt composition can affect the melt density at such extreme conditions and may strongly influence the structural response.
The H alpha spectral line is a well-studied absorption line revealing properties of the highly structured and dynamic solar chromosphere. Typical features with distinct spectral signatures in H alpha include filaments and prominences, bright active-region plages, superpenumbrae around sunspots, surges, flares, Ellerman bombs, filigree, and mottles and rosettes, among others. This study is based on high-spectral resolution H alpha spectra obtained with the Echelle spectrograph of the Vacuum Tower Telescope (VTT) located at Observatorio del Teide, Tenerife, Spain. The t-distributed stochastic neighbor embedding (t-SNE) is a machine-learning algorithm, which is used for nonlinear dimensionality reduction. In this application, it projects H alpha spectra onto a two-dimensional map, where it becomes possible to classify the spectra according to results of cloud model (CM) inversions. The CM parameters optical depth, Doppler width, line-of-sight velocity, and source function describe properties of the cloud material. Initial results of t-SNE indicate its strong discriminatory power to separate quiet-Sun and plage profiles from those that are suitable for CM inversions. In addition, a detailed study of various t-SNE parameters is conducted, the impact of seeing conditions on the classification is assessed, results for various types of input data are compared, and the identified clusters are linked to chromospheric features. Although t-SNE proves to be efficient in clustering high-dimensional data, human inference is required at each step to interpret the results. This exploratory study provides a framework and ideas on how to tailor a classification scheme toward specific spectral data and science questions.
Context
Thermonuclear supernovae (SNe), a subset of which are the highly important SNe Type Ia, remain one of the more poorly understood phenomena known to modern astrophysics.
In recent years, the single degenerate helium (He) donor channel, where a white dwarf star accretes He-rich matter from a hydrogen-depleted companion, has emerged as a promising candidate progenitor scenario for these events.
An unresolved question in this scenario is the fate of the companion star, which would be evident as a runaway hot subdwarf O/B stars (He sdO/B) in the aftermath of the SN event.
Aims
Previous studies have shown that the kinematic properties of an ejected companion provide an opportunity to closer examine the properties of an SN progenitor system. However, with the number of observed objects not matching predictions by theory, the viability of this mechanism is called into question.
In this study, we first synthesize a population of companion stars ejected by the aforementioned mechanism, taking into account predicted ejection velocities, the inferred population density in the Galactic mass distribution, and subsequent kinematics in the Galactic potential. We then discuss the astrometric properties of this population.
Methods
We present 10(6) individual ejection trajectories, which were numerically computed with a newly developed, lightweight simulation framework. Initial conditions were randomly generated, but weighted according to the Galactic mass density and ejection velocity data. We then discuss the bulk properties (Galactic distribution and observational parameters) of our sample.
Results
Our synthetic population reflects the Galactic mass distribution.
A peak in the density distribution for close objects is expected in the direction of the Galactic centre. Higher mass runaways should outnumber lower mass ones. If the entire considered mass range is realised, the radial velocity distribution should show a peak at 500 km s(-1).
If only close US 708 analogues are considered, there should be a peak at (similar to 750-850) km s(-1). In either case, US 708 should be a member of the high-velocity tail of the distribution.
Conclusions
We show that the puzzling lack of confirmed surviving companion stars of thermonuclear SNe, though possibly an observation-related selection effect, may indicate a selection against high mass donors in the SD He donor channel.
Solar H alpha excess during Solar Cycle 24 from full-disk filtergrams of the Chromospheric Telescope
(2022)
Context
The chromospheric H alpha spectral line is a strong line in the spectrum of the Sun and other stars. In the stellar regime, this spectral line is already used as a powerful tracer of stellar activity. For the Sun, other tracers, such as Ca II K, are typically used to monitor solar activity. Nonetheless, the Sun is observed constantly in H alpha with globally distributed ground-based full-disk imagers.
Aims
The aim of this study is to introduce the imaging H alpha excess and deficit as tracers of solar activity and compare them to other established indicators. Furthermore, we investigate whether the active region coverage fraction or the changing H alpha excess in the active regions dominates temporal variability in solar H alpha observations.
Methods
We used observations of full-disk H alpha filtergrams of the Chromospheric Telescope and morphological image processing techniques to extract the imaging H alpha excess and deficit, which were derived from the intensities above or below 10% of the median intensity in the filtergrams, respectively. These thresholds allowed us to filter for bright features (plage regions) and dark absorption features (filaments and sunspots). In addition, the thresholds were used to calculate the mean intensity I-mean(E/D) for H alpha excess and deficit regions. We describe the evolution of the H alpha excess and deficit during Solar Cycle 24 and compare it to the mean intensity and other well established tracers: the relative sunspot number, the F10.7 cm radio flux, and the Mg II index. In particular, we tried to determine how constant the H alpha excess and number density of H alpha excess regions are between solar maximum and minimum. The number of pixels above or below the intensity thresholds were used to calculate the area coverage fraction of H alpha excess and deficit regions on the Sun, which was compared to the imaging H alpha excess and deficit and the respective mean intensities averaged for the length of one Carrington rotation. In addition, we present the H alpha excess and mean intensity variation of selected active regions during their disk passage in comparison to the number of pixels of H alpha excess regions.
Results.
The H alpha excess and deficit follow the behavior of the solar activity over the course of the cycle. They both peak around solar maximum, whereby the peak of the H alpha deficit is shortly after the solar maximum. Nonetheless, the correlation of the monthly averages of the H alpha excess and deficit is high with a Spearman correlation of rho = 0.91. The H alpha excess is closely correlated to the chromospheric Mg II index with a correlation of 0.95. The highest correlation of the H alpha deficit is found with the F10.7 cm radio flux, with a correlation of 0.89, due to their peaks after the solar activity maximum. Furthermore, the H alpha deficit reflects the cyclic behavior of polar crown filaments and their disappearance shortly before the solar maximum. We investigated the mean intensity distribution for H alpha excess regions for solar minimum and maximum. The shape of the distributions for solar minimum and maximum is very similar, but with different amplitudes. Furthermore, we found that the area coverage fraction of H alpha excess regions and the H alpha excess are strongly correlated with an overall Spearman correlation of 0.92. The correlation between the H alpha excess and the mean intensity of H alpha excess regions is 0.75. The correlation of the area coverage fraction and the mean intensity of H alpha excess regions is in general relatively low (rho = 0.45) and only for few active regions is this correlation above 0.7. The weak correlation between the area coverage fraction and mean intensity leaves us pessimistic that the degeneracy between these two quantities can be broken for the modeling of unresolved stellar surfaces.
We study the first-arrival (first-hitting) dynamics and efficiency of a one-dimensional random search model performing asymmetric Levy flights by leveraging the Fokker-Planck equation with a delta-sink and an asymmetric space-fractional derivative operator with stable index alpha and asymmetry (skewness) parameter beta.
We find exact analytical results for the probability density of first-arrival times and the search efficiency, and we analyse their behaviour within the limits of short and long times.
We find that when the starting point of the searcher is to the right of the target, random search by Brownian motion is more efficient than Levy flights with beta <= 0 (with a rightward bias) for short initial distances, while for beta>0 (with a leftward bias) Levy flights with alpha -> 1 are more efficient.
When increasing the initial distance of the searcher to the target, Levy flight search (except for alpha=1 with beta=0) is more efficient than the Brownian search. Moreover, the asymmetry in jumps leads to essentially higher efficiency of the Levy search compared to symmetric Levy flights at both short and long distances, and the effect is more pronounced for stable indices alpha close to unity.
EC 22536-5304
(2021)
Helium-burning hot subdwarf stars of spectral types O and B (sdO/B) are thought to be produced through various types of binary interactions. The helium-rich hot subdwarf star EC 22536-5304 was recently found to be extremely enriched in lead. Here, we show that EC 22536-5304 is a binary star with a metal-poor subdwarf F-type (sdF) companion. We performed a detailed analysis of high-resolution SALT/HRS and VLT/UVES spectra, deriving metal abundances for the hot subdwarf, as well as atmospheric parameters for both components. Because we consider the contribution of the sdF star, the derived lead abundance for the sdOB, + 6.3 +/- 0.3 dex relative to solar, is even higher than previously thought. We derive T-eff = 6210 +/- 70 K, log g = 4.64 +/- 0.10, [FE/H] = - 1.95 +/- 0.04, and [alpha/Fe] = + 0.40 +/- 0.04 for the sdF component. Radial velocity variations, although poorly sampled at present, indicate that the binary system has a long orbital period of about 457 days. This suggests that the system was likely formed through stable Roche lobe overflow (RLOF). A kinematic analysis shows that EC 22536-5304 is on an eccentric orbit around the Galactic centre. This, as well as the low metallicity and strong alpha enhancement of the sdF-type companion, indicate that EC 22536-5304 is part of the Galactic halo or metal-weak thick disc. As the first long-period hot subdwarf binary at [FE/H] less than or similar to- 1, EC 22536-5304 may help to constrain the RLOF mechanism for mass transfer from low-mass, low-metallicity red giant branch (RGB) stars to main-sequence companions.
Thermal electrons have gyroradii many orders of magnitude smaller than the finite width of a shock, thus need to be pre-accelerated before they can cross it and be accelerated by diffusive shock acceleration. One region where pre-acceleration may occur is the inner foreshock, which upstream electrons must pass through before any potential downstream crossing. In this paper, we perform a large-scale particle-in-cell simulation that generates a single shock with parameters motivated from supernova remnants. Within the foreshock, reflected electrons excite the oblique whistler instability and produce electromagnetic whistler waves, which comove with the upstream flow and as nonlinear structures eventually reach radii of up to 5 ion-gyroradii. We show that the inner electromagnetic configuration of the whistlers evolves into complex nonlinear structures bound by a strong magnetic field around four times the upstream value. Although these nonlinear structures do not in general interact with cospatial upstream electrons, they resonate with electrons that have been reflected at the shock. We show that they can scatter, or even trap, reflected electrons, confining around 0.8% of the total upstream electron population to the region close to the shock where they can undergo substantial pre-acceleration. This acceleration process is similar to, yet approximately three times more efficient than, stochastic shock drift acceleration.
Assessing the impact of hydrogen absorption on the characteristics of the Galactic center excess
(2022)
We present a new reconstruction of the distribution of atomic hydrogen in the inner Galaxy that is based on explicit radiation transport modeling of line and continuum emission and a gas-flow model in the barred Galaxy that provides distance resolution for lines of sight toward the Galactic center.
The main benefits of the new gas model are (a) the ability to reproduce the negative line signals seen with the HI4PI survey and (b) the accounting for gas that primarily manifests itself through absorption.
We apply the new model of Galactic atomic hydrogen to an analysis of the diffuse gamma-ray emission from the inner Galaxy, for which an excess at a few GeV was reported that may be related to dark matter.
We find with high significance an improved fit to the diffuse gamma-ray emission observed with the Fermi-LAT, if our new H i model is used to estimate the cosmic-ray induced diffuse gamma-ray emission.
The fit still requires a nuclear bulge at high significance. Once this is included there is no evidence of a dark-matter signal, be it cuspy or cored. But an additional so-called boxy bulge is still favored by the data.
This finding is robust under the variation of various parameters, for example, the excitation temperature of atomic hydrogen, and a number of tests for systematic issues.
Gravitational waves from the collision of binary neutron stars provide a unique opportunity to study the behaviour of supranuclear matter, the fundamental properties of gravity and the cosmic history of our Universe. However, given the complexity of Einstein's field equations, theoretical models that enable source-property inference suffer from systematic uncertainties due to simplifying assumptions. We develop a hypermodel approach to compare and measure the uncertainty of gravitational-wave approximants. Using state-of-the-art models, we apply this new technique to the binary neutron star observations GW170817 and GW190425 and to the sub-threshold candidate GW200311_103121. Our analysis reveals subtle systematic differences (with Bayesian odds of similar to 2) between waveform models. A frequency-dependence study suggests that this may be due to the treatment of the tidal sector. This new technique provides a proving ground for model development and a means to identify waveform systematics in future observing runs where detector improvements will increase the number and clarity of binary neutron star collisions we observe.
The current paradigm of cosmic-ray (CR) origin states that the greater part of galactic CRs is produced by supernova remnants. The interaction of supernova ejecta with the interstellar medium after a supernova's explosions results in shocks responsible for CR acceleration via diffusive shock acceleration (DSA). We use particle-in-cell (PIC) simulations and a combined PIC-magnetohydrodynamic (PIC-MHD) technique to investigate whether DSA can occur in oblique high Mach number shocks. Using the PIC method, we follow the formation of the shock and determine the fraction of the particles that gets involved in DSA. With this result, we use PIC-MHD simulations to model the large-scale structure of the plasma and the magnetic field surrounding the shock and find out whether or not the reflected particles can generate upstream turbulence and trigger DSA. We find that the feasibility of this process in oblique shocks depends strongly on the Alfvenic Mach number, and the DSA process is more likely to be triggered at high Mach number shocks.
Most image restoration methods in astronomy rely upon probabilistic tools that infer the best solution for a deconvolution problem. They achieve good performances when the point spread function (PSF) is spatially invariant in the image plane.
However, this condition is not always satisfied in real optical systems. We propose a new method for the restoration of images affected by static and anisotropic aberrations using Deep Neural Networks that can be directly applied to sky images.
The network is trained using simulated sky images corresponding to the T80-S Telescope optical model, a 80-cm survey imager at Cerro Tololo (Chile), which are synthesized using a Zernike polynomial representation of the optical system.
Once trained, the network can be used directly on sky images, outputting a corrected version of the image that has a constant and known PSF across its field of view. The method is to be tested on the T80-S Telescope.
We present the method and results on synthetic data.
3A 1954+319 has been classified for a long time as a symbiotic X-ray binary, hosting a slowly rotating neutron star and an aged M red giant. Recently, this classification has been revised thanks to the discovery that the donor star is an M supergiant. This makes 3A 1954+319 a rare type of high-mass X-ray binary consisting of a neutron star and a red supergiant donor. In this paper, we analyse two archival and still unpublished XMM-Newton and NuSTAR observations of the source. We perform a detailed hardness ratio-resolved spectral analysis to search for spectral variability that could help investigating the structures of the inhomogeneous M supergiant wind from which the neutron star is accreting. We discuss our results in the context of wind-fed supergiant X-ray binaries and show that the newest findings on 3A 1954+319 reinforce the hypothesis that the neutron star in this system is endowed with a magnetar-like magnetic field strength (greater than or similar to 10(14) G).
Context. The spectroscopic class of subdwarf A-type (sdA) stars has come into focus in recent years because of their possible link to extremely low-mass white dwarfs, a rare class of objects resulting from binary evolution. Although most sdA stars are consistent with metal-poor halo main-sequence stars, the formation and evolution of a fraction of these stars are still matters of debate. Aims. The identification of photometric variability can help to put further constraints on the evolutionary status of sdA stars, in particular through the analysis of pulsations. Moreover, the binary ratio, which can be deduced from eclipsing binaries and ellipsoidal variables, is important as input for stellar models. In order to search for variability due to either binarity or pulsations in objects of the spectroscopic sdA class, we have extracted all available high precision light curves from the Kepler K2 mission.
Methods. We have performed a thorough time series analysis on all available light curves, employing three different methods. Frequencies with a signal-to-noise ratio higher than four have been used for further analysis.
Results. From the 25 targets, 13 turned out to be variables of different kinds (i.e., classical pulsating stars, ellipsoidal and cataclysmic variables, eclipsing binaries, and rotationally induced variables). For the remaining 12 objects, a variability threshold was determined.
The mean free path of ionizing photons, lambda(mfp), is a key factor in the photoionization of the intergalactic medium (IGM). At z greater than or similar to 5, however, lambda(mfp) may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of lambda(mfp) that address this bias and extend up to z similar to 6 for the first time. Our measurements at z similar to 5 are based on data from the Giant Gemini GMOS survey and new Keck LRIS observations of low-luminosity QSOs. At z similar to 6 we use QSO spectra from Keck ESI and VLT X-Shooter. We measure lambda(mfp) = 9.09(-1.28)(+1.62) proper Mpc and 0.75(-0.45)(+0.65) proper Mpc (68 percent confidence) at z = 5.1 and 6.0, respectively. The results at z = 5.1 are consistent with existing measurements, suggesting that bias from the proximity effect is minor at this redshift. At z = 6.0, however, we find that neglecting the proximity effect biases the result high by a factor of two or more. Our measurement at z = 6.0 falls well below extrapolations from lower redshifts, indicating rapid evolution in lambda(mfp) over 5 < z < 6. This evolution disfavours models in which reionization ended early enough that the IGM had time to fully relax hydrodynamically by z = 6, but is qualitatively consistent with models wherein reionization completed at z = 6 or even significantly later. Our mean free path results are most consistent with late reionization models wherein the IGM is still 20 percent neutral at z = 6, although our measurement at z = 6.0 is even lower than these models prefer.
Symbiotic X-ray binaries are systems hosting a neutron star accreting form the wind of a late-type companion. These are rare objects and so far only a handful of them are known. One of the most puzzling aspects of the symbiotic X-ray binaries is the possibility that they contain strongly magnetized neutron stars. These are expected to be evolutionary much younger compared to their evolved companions and could thus be formed through the (yet poorly known) accretion induced collapse of a white dwarf. In this paper, we perform a broad-band X-ray and soft gamma-ray spectroscopy of two known symbiotic binaries, Sct X-1 and 4U 1700+24, looking for the presence of cyclotron scattering features that could confirm the presence of strongly magnetized NSs. We exploited available Chandra, Swift, and NuSTAR data. We find no evidence of cyclotron resonant scattering features (CRSFs) in the case of Sct X-1 but in the case of 4U 1700+24 we suggest the presence of a possible CRSF at similar to 16 keV and its first harmonic at similar to 31 keV, although we could not exclude alternative spectral models for the broad-band fit. If confirmed by future observations, 4U 1700+24 could be the second symbiotic X-ray binary with a highly magnetized accretor. We also report about our long-term monitoring of the last discovered symbiotic X-ray binary IGR J17329-2731 performed with Swift/XRT. The monitoring revealed that, as predicted, in 2017 this object became a persistent and variable source, showing X-ray flares lasting for a few days and intriguing obscuration events that are interpreted in the context of clumpy wind accretion.
The investigation of metabolic fluxes and metabolite distributions within cells by means of tracer molecules is a valuable tool to unravel the complexity of biological systems. Technological advances in mass spectrometry (MS) technology such as atmospheric pressure chemical ionization (APCI) coupled with high resolution (HR), not only allows for highly sensitive analyses but also broadens the usefulness of tracer-based experiments, as interesting signals can be annotated de novo when not yet present in a compound library. However, several effects in the APCI ion source, i.e., fragmentation and rearrangement, lead to superimposed mass isotopologue distributions (MID) within the mass spectra, which need to be corrected during data evaluation as they will impair enrichment calculation otherwise. Here, we present and evaluate a novel software tool to automatically perform such corrections. We discuss the different effects, explain the implemented algorithm, and show its application on several experimental datasets. This adjustable tool is available as an R package from CRAN.
The nature of the sources powering nebular He II emission in star-forming galaxies remains debated, and various types of objects have been considered, including Wolf-Rayet stars, X-ray binaries, and Population III stars.
Modern X-ray observations show the ubiquitous presence of hot gas filling star-forming galaxies. We use a collisional ionization plasma code to compute the specific He II ionizing flux produced by hot gas and show that if its temperature is not too high (less than or similar to 2.5 MK), then the observed levels of soft diffuse X-ray radiation could explain He II ionization in galaxies.
To gain a physical understanding of this result, we propose a model that combines the hydrodynamics of cluster winds and hot superbubbles with observed populations of young massive clusters in galaxies. We find that in low-metallicity galaxies, the temperature of hot gas is lower and the production rate of He II ionizing photons is higher compared to high-metallicity galaxies. The reason is that the slower stellar winds of massive stars in lower-metallicity galaxies input less mechanical energy in the ambient medium.
Furthermore, we show that ensembles of star clusters up to similar to 10-20 Myr old in galaxies can produce enough soft X-rays to induce nebular He II emission. We discuss observations of the template low-metallicity galaxy I Zw 18 and suggest that the He II nebula in this galaxy is powered by a hot superbubble.
Finally, appreciating the complex nature of stellar feedback, we suggest that soft X-rays from hot superbubbles are among the dominant sources of He II ionizing flux in low-metallicity star-forming galaxies.
Starting from the observation that the reduced state of a system strongly coupled to a bath is, in general, an athermal state, we introduce and study a cyclic battery-charger quantum device that is in thermal equilibrium, or in a ground state, during the charge storing stage. The cycle has four stages: the equilibrium storage stage is interrupted by disconnecting the battery from the charger, then work is extracted from the battery, and then the battery is reconnected with the charger; finally, the system is brought back to equilibrium. At no point during the cycle are the battery-charger correlations artificially erased. We study the case where the battery and charger together comprise a spin-1/2 Ising chain, and show that the main characteristics-the extracted energy and the thermodynamic efficiency-can be enhanced by operating the cycle close to the quantum phase transition point. When the battery is just a single spin, we find that the output work and efficiency show a scaling behavior at criticality and derive the corresponding critical exponents. Due to always present correlations between the battery and the charger, operations that are equivalent from the perspective of the battery can entail different energetic costs for switching the battery-charger coupling. This happens only when the coupling term does not commute with the battery's bare Hamiltonian, and we use this purely quantum leverage to further optimize the performance of the device.
The influence of the process gas, laser scan speed, and sample thickness on the build-up of residual stresses and porosity in Ti-6Al-4V produced by laser powder bed fusion was studied. Pure argon and helium, as well as a mixture of those (30% helium), were employed to establish process atmospheres with a low residual oxygen content of 100 ppm O-2. The results highlight that the subsurface residual stresses measured by X-ray diffraction were significantly lower in the thin samples (220 MPa) than in the cuboid samples (645 MPa). This difference was attributed to the shorter laser vector length, resulting in heat accumulation and thus in-situ stress relief. The addition of helium to the process gas did not introduce additional subsurface residual stresses in the simple geometries, even for the increased scanning speed. Finally, larger deflection was found in the cantilever built under helium (after removal from the baseplate), than in those produced under argon and an argon-helium mixture. This result demonstrates that complex designs involving large scanned areas could be subjected to higher residual stress when manufactured under helium due to the gas's high thermal conductivity, heat capacity, and thermal diffusivity.
In light of substantial new discoveries of hot subdwarfs by ongoing spectroscopic surveys and the availability of the Gaia mission Early Data Release 3 (EDR3), we compiled new releases of two catalogues of hot subluminous stars: the data release 3 (DR3) catalogue of the known hot subdwarf stars contains 6616 unique sources and provides multi-band photometry, and astrometry from Gaia EDR3 as well as classifications based on spectroscopy and colours.
This is an increase of 742 objects over the DR2 catalogue.
This new catalogue provides atmospheric parameters for 3087 stars and radial velocities for 2791 stars from the literature. In addition, we have updated the Gaia Data Release 2 (DR2) catalogue of hot subluminous stars using the improved accuracy of the Gaia EDR3 data set together with updated quality and selection criteria to produce the Gaia EDR3 catalogue of 61 585 hot subluminous stars, representing an increase of 21 785 objects.
The improvements in Gaia EDR3 astrometry and photometry compared to Gaia DR2 have enabled us to define more sophisticated selection functions.
In particular, we improved hot subluminous star detection in the crowded regions of the Galactic plane as well as in the direction of the Magellanic Clouds by including sources with close apparent neighbours but with flux levels that dominate the neighbourhood.
As part of a survey to find close binary systems among central stars of planetary nebula, we present two newly discovered binary systems.
GALEX J015054.4+310745 is identified as the central star of the possible planetary nebula Fr 2-22. We find it to be a single-lined spectroscopic binary with an orbital period of 0.2554435(10) d.
We support the previous identification of GALEX J015054.4+310745 as an sdB star and provide physical parameters for the star from spectral modelling. We identify its undetected companion as a likely He white dwarf. Based on this information, we find it unlikely that Fr 2-22 is a true planetary nebula. In addition, the central star of the true planetary nebula Hen 2-84 is found to be a photometric variable, likely due to the irradiation of a cool companion. The system has an orbital period of 0.485645(30) d.
We discuss limits on binary parameters based on the available light-curve data.
Hen 2-84 is a strongly shaped bipolar planetary nebula, which we now add to the growing list of axially or point-symmetric planetary nebulae with a close binary central star.
In this paper, we present a study comparing the depth to diameter (d/D) ratio of small simple craters (200-1000 m) of an area between -88.5 degrees to -90 degrees latitude at the lunar south pole containing Permanent Shadowed Regions (PSRs) versus craters without PSRs. As PSRs can reach temperatures of 110 K and are capable of harboring volatiles, especially water ice, we analyzed the relationship of depth versus diameter ratios and its possible implications for harboring water ice. Variations in the d/D ratios can also be caused by other processes such as degradation, isostatic adjustment, or differences in surface properties. The conducted d/D ratio analysis suggests that a differentiation between craters containing PSRs versus craters without PSRs occurs. Thus, a possible direct relation between d/D ratio, PSRs, and water ice harboring might exist. Our results suggest that differences in the target's surface properties may explain the obtained results. The resulting d/D ratios of craters with PSRs can help to select target areas for future In-Situ Resource Utilization (ISRU) missions.
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.