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.
Thanks to the combined effort s of scientist s in several research fields, the preceding decade has witnessed considerable progress in the use of conjugated polymers as emerging thermoelectric materials leading to significant improvements in performance and demonstration of a number of diverse applications. Despite these recent advances, systematic assessments of the impact of molecular design on thermoelectric properties are scarce. Although several reviews marginally highlight the role of chemical structure, the understanding of structure-performance relationships is still fragmented. An in-depth understanding of the relationship between molecular structure and thermoelectric properties will enable the rational design of next-generation thermoelectric polymers. To this end, this review showcases the state-of-the-art thermoelectric polymers, discusses structure-performance relationships, suggests strategies for improving thermoelectric performance that go beyond molecular design, and highlights some of the most impressive applications of thermoelectric polymers.
Planets that closely orbit magnetically active stars are thought to be able to interact with their magnetic fields in a way that modulates stellar activity. This modulation in phase with the planetary orbit, such as enhanced X-ray activity, chromospheric spots, radio emission, or flares, is considered the clearest sign of magnetic star-planet interaction (SPI). However, the magnitude of this interaction is poorly constrained, and the intermittent nature of the interaction is a challenge for observers. AU Mic is an early M dwarf, and the most actively flaring planet host detected to date. Its innermost companion, AU Mic b, is a promising target for magnetic SPI observations. We used optical light curves of AU Mic obtained by the Transiting Exoplanet Survey Satellite to search for signs of flaring SPI with AU Mic b using a customized Anderson-Darling test. In the about 50 d of observations, the flare distributions with orbital, rotational, and synodic periods were generally consistent with intrinsic stellar flaring. We found the strongest deviation (p = 0.07, n = 71) from intrinsic flaring with the orbital period of AU Mic b, in the high-energy half of our sample (ED > 1 s). If it reflects the true SPI signal from AU Mic b, extending the observing time by a factor of 2-3 will yield a >3 sigma detection. Continued monitoring of AU Mic may therefore reveal flaring SPI with orbital phase, while rotational modulation will smear out due to the star's strong differential rotation.
Hot subdwarf stars represent a late and peculiar stage in the evolution of low-mass stars, since they are likely formed by close binary interactions. In this work, we perform a radial velocity (RV) variability study of a sample of 646 hot subdwarfs with multi-epoch radial velocities based on spectra from Sloan Digital Sky Survey (SDSS) and Large Sky Area Multi-Object Fibre Spectroscopic Telescope (LAMOST). The atmospheric parameters and RVs were taken from the literature. For stars with archival spectra but without literature values, we determined the parameters by fitting model atmospheres. In addition, we redetermined the atmospheric parameters and RVs for all the He-enriched sdO/Bs. This broad sample allowed us to study RV-variability as a function of the location in the T-eff - log g- and T-eff - log n(He)/n(H) diagrams in a statistically significant way. We used the fraction of RV-variable stars and the distribution of the maximum RV variations Delta RVmax as diagnostics. Both indicators turned out to be quite inhomogeneous across the studied parameter ranges. A striking feature is the completely dissimilar behaviour of He-poor and He-rich hot subdwarfs. While the former have a high fraction of close binaries, almost no significant RV variations could be detected for the latter. This has led us to the conclusion that there is likely no evolutionary connection between these subtypes. On the other hand, intermediate He-rich- and extreme He-rich sdOB/Os are more likely to be related. Furthermore, we conclude that the vast majority of this population is formed via one or several binary merger channels. Hot subdwarfs with temperatures cooler than similar to 24 000 K tend to show fewer and smaller RV-variations. These objects might constitute a new subpopulation of binaries with longer periods and late-type or compact companions. The RV-variability properties of the extreme horizontal branch (EHB) and corresponding post-EHB populations of the He-poor hot subdwarfs match and confirm the predicted evolutionary connection between them. Stars found below the canonical EHB at somewhat higher surface gravities show large RV variations and a high RV variability fraction. These properties are consistent with most of them being low-mass EHB stars or progenitors of low-mass helium white dwarfs in close binaries.
In this study, we model a sequence of a confined and a full eruption, employing the relaxed end state of the confined eruption of a kink-unstable flux rope as the initial condition for the ejective one. The full eruption, a model of a coronal mass ejection, develops as a result of converging motions imposed at the photospheric boundary, which drive flux cancellation. In this process, parts of the positive and negative external flux converge toward the polarity inversion line, reconnect, and cancel each other. Flux of the same amount as the canceled flux transfers to a flux rope, increasing the free magnetic energy of the coronal field. With sustained flux cancellation and the associated progressive weakening of the magnetic tension of the overlying flux, we find that a flux reduction of approximate to 11% initiates the torus instability of the flux rope, which leads to a full eruption. These results demonstrate that a homologous full eruption, following a confined one, can be driven by flux cancellation.
Context.
About 10% of all stars exhibit absorption lines of ultra-highly excited (UHE) metals (e.g., O VIII) in their optical spectra when entering the white dwarf cooling sequence. This is something that has never been observed in any other astrophysical object, and poses a decades-long mystery in our understanding of the late stages of stellar evolution. The recent discovery of a UHE white dwarf that is both spectroscopically and photometrically variable led to the speculation that the UHE lines might be created in a shock-heated circumstellar magnetosphere.
Aims.
We aim to gain a better understanding of these mysterious objects by studying the photometric variability of the whole population of UHE white dwarfs, and white dwarfs showing only the He II line problem, as both phenomena are believed to be connected.
Methods.
We investigate (multi-band) light curves from several ground- and space-based surveys of all 16 currently known UHE white dwarfs (including one newly discovered) and eight white dwarfs that show only the He II line problem.
Results.
We find that 75(-13)(+8) % of the UHE white dwarfs, and 75(-19)(+9)% of the He II line problem white dwarfs are significantly photometrically variable, with periods ranging from 0.22 d to 2.93 d and amplitudes from a few tenths to a few hundredths of a magnitude. The high variability rate is in stark contrast to the variability rate amongst normal hot white dwarfs (we find 9(2)(+4)%), marking UHE and He II line problem white dwarfs as a new class of variable stars. The period distribution of our sample agrees with both the orbital period distribution of post-common-envelope binaries and the rotational period distribution of magnetic white dwarfs if we assume that the objects in our sample will spin-up as a consequence of further contraction.
Conclusions.
We find further evidence that UHE and He II line problem white dwarfs are indeed related, as concluded from their overlap in the Gaia HRD, similar photometric variability rates, light-curve shapes and amplitudes, and period distributions. The lack of increasing photometric amplitudes towards longer wavelengths, as well as the nondetection of optical emission lines arising from the highly irradiated face of a hypothetical secondary in the optical spectra of our stars, makes it seem unlikely that an irradiated late-type companion is the origin of the photometric variability. Instead, we believe that spots on the surfaces of these stars and/or geometrical effects of circumstellar material might be responsible.
Recently, glasses, a subset of amorphous solids, have gained attention in various fields, such as polymer chemistry, optical fibers, and pharmaceuticals. One of their characteristic features, the glass transition temperature (T-g) which is absent in 100% crystalline materials, influences several material properties, such as free volume, enthalpy, viscosity, thermodynamic transitions, molecular motions, physical stability, mechanical properties, etc. In addition to T-g, there may be several other temperaturedependent transitions known as sub-T-g transitions (or beta-, gamma-, and delta-relaxations) which are identified by specific analytical techniques. The study of T-g and sub-T-g transitions occurring in amorphous solids has gained much attention because of its importance in understanding molecular kinetics, and it requires the combination of conventional and novel characterization techniques. In the present study, three different analytical techniques [modulated differential scanning calorimetry (mDSC), dynamic mechanical analysis (DMA), and dielectric relaxation spectroscopy (DRS)] were used to perform comprehensive qualitative/quantitative characterization of molecular relaxations, miscibility, and molecular interactions present in an amorphous polymer (PVPVA), a model drug (indomethacin, IND), and IND/PVPVA-based amorphous solid dispersions (ASDs). This is the first ever reported DMA study on PVPVA in its powder form, which avoids the contribution of solvent to the mechanical properties when a selfstanding polymer film is used. A good correlation between the techniques in determining the T-g value of PVPVA, IND, and IND/ PVPVA-based ASDs is established, and the negligible difference (within 10 degrees C) is attributed to the different material properties assessed in each technique. However, the overall T-g behavior, the decrease in T-g with increase in drug loading in ASDs, is universally observed in all the above-mentioned techniques, which reveals their complementarity. DMA and DRS techniques are used to study the different sub-T-g transitions present in PVPVA, amorphous IND, and IND/PVPVA-based ASDs because these transitions are normally too weak or too broad for mDSC to detect. For IND/PVPVA-based ASDs, both techniques show a shift of sub-T-g transitions (or secondary relaxation peaks) toward the high-temperature region from -140 to -45 degrees C. Thus, this paper outlines the usage of different solid-state characterization techniques in understanding the different molecular dynamics present in the polymer, drug, and their interactions in ASDs with the integrated information obtained from individual techniques.
Theory predicts that the temperature of the X-ray-emitting gas (similar to 10(6) K) detected from planetary nebulae (PNe) is a consequence of mixing or thermal conduction when in contact with the ionized outer rim (similar to 10(4) K). Gas at intermediate temperatures (similar to 10(5) K) can be used to study the physics of the production of X-ray-emitting gas, via C iv, N v, and O vi ions. Here, we model the stellar atmosphere of the CSPN of NGC 1501 to demonstrate that even this hot H-deficient [WO4]-type star cannot produce these emission lines by photoionization. We use the detection of the C iv lines to assess the physical properties of the mixing region in this PNe in comparison with its X-ray-emitting gas, rendering NGC 1501 only the second PNe with such characterization. We extend our predictions to the hottest [WO1] and cooler [WC5] spectral types and demonstrate that most energetic photons are absorbed in the dense winds of [WR] CSPN and highly ionized species can be used to study the physics behind the production of hot bubbles in PNe. We found that the UV observations of NGC 2452, NGC 6751, and NGC 6905 are consistent with the presence mixing layers and hot bubbles, providing excellent candidates for future X-ray observations.
Context.
The supergiant ionized shell SMC-SGS 1 (DEM 167), which is located in the outer Wing of the Small Magellanic Cloud (SMC), resembles structures that originate from an energetic star-formation event and later stimulate star formation as they expand into the ambient medium. However, stellar populations within and surrounding SMC-SGS 1 tell a different story. Aims. We present a photometric study of the stellar population encompassed by SMC-SGS 1 in order to trace the history of such a large structure and its potential influence on star formation within the low-density, low-metallicity environment of the SMC.
Methods.
For a stellar population that is physically associated with SMC-SGS 1, we combined near-ultraviolet (NUV) photometry from the Galaxy Evolution Explorer with archival optical (V-band) photometry from the ESO Danish 1.54 m Telescope. Given their colors and luminosities, we estimated stellar ages and masses by matching observed photometry to theoretical stellar isochrone models. Results. We find that the investigated region supports an active, extended star-formation event spanning similar to 25-40 Myr ago, as well as continued star formation into the present. Using a standard initial mass function, we infer a lower bound on the stellar mass from this period of similar to 3 x 10(4) M-circle dot, corresponding to a star-formation intensity of similar to 6 x 10(-3) M-circle dot kpc(-2) yr(-1).
Conclusions.
The spatial and temporal distributions of young stars encompassed by SMC-SGS 1 imply a slow, consistent progression of star formation over millions of years. Ongoing star formation, both along the edge and interior to SMC-SGS 1, suggests a combined stimulated and stochastic mode of star formation within the SMC Wing. We note that a slow expansion of the shell within this low-density environment may preserve molecular clouds within the volume of the shell, leaving them to form stars even after nearby stellar feedback expels local gas and dust.
A comprehensive photometric and spectroscopic analysis of the variable TYC 5532-1333-1 (TYC) along with an investigation of its orbital period variation is presented for the first time. The B- and V-band photometric study indicates that TYC is an intermediate contact binary with degree of contact and mass ratio of 34 per cent and similar to 0.24, respectively. The derived equivalent widths from the spectroscopic study of H alpha and Na-I lines reveal phase-dependent variation and mutual correlation. Using the available times of minimum light, an investigation of orbital period variation shows a long-term decrease at a rate of 3.98 x 10(-6) d yr(-1). Expected causes for such decline in the orbital period could be angular momentum loss and a quasi-sinusoidal variation due to light-time effect probably caused by a third-body companion. The minimum mass of the third body (M-3) was derived to be 0.65 M-circle dot. Our presented study is an attempt to evaluate and understand the evolutionary state of above-mentioned neglected contact binary.
An extinction-free estimator of the star formation rate (SFR) of galaxies is critical for understanding the high-redshift universe. To this end, the nearly linear, tight correlation of far-infrared (FIR), and radio luminosity of star-forming galaxies is widely used. While the FIR is linked to massive star formation, which also generates shock-accelerated cosmic-ray (CR) electrons and radio synchrotron emission, a detailed understanding of the underlying physics is still lacking. Hence, we perform three-dimensional magnetohydrodynamical (MHD) simulations of isolated galaxies over a broad range of halo masses and SFRs using the moving-mesh code AREPO, and evolve the CR proton energy density self-consistently. In post-processing, we calculate the steady-state spectra of primary, shock-accelerated and secondary CR electrons, which result from hadronic CR proton interactions with the interstellar medium. The resulting total radio luminosities correlate with the FIR luminosities as observed and are dominated by primary CR electrons if we account for anisotropic CR diffusion. The increasing contribution of secondary emission up to 30 per cent in starbursts is compensated by the larger bremsstrahlung and Coulomb losses. CR electrons are in the calorimetric limit and lose most of their energy through inverse Compton interactions with star light and cosmic microwave background (CMB) photons while less energy is converted into synchrotron emission. This implies steep steady-state synchrotron spectra in starbursts. Interestingly, we find that thermal free-free emission flattens the total radio spectra at high radio frequencies and reconciles calorimetric theory with observations while free-free absorption explains the observed low-frequency flattening towards the central regions of starbursts.
SwSt 1 (PN G001.5-06.7) is a bright and compact planetary nebula containing a late [WC]-type central star. Previous studies suggested that the nebular and stellar lines are slowly changing with time. We studied new and archival optical and ultraviolet spectra of the object. The [O III] 4959 and 5007 angstrom to H beta line flux ratios decreased between about 1976 and 1997/2015. The stellar spectrum also shows changes between these epochs. We modelled the stellar and nebular spectra observed at different epochs. The analyses indicate a drop of the stellar temperature from about 42 kK to 40.5 kK between 1976 and 1993. We do not detect significant changes between 1993 and 2015. The observations show that the star performed a loop in the H-R diagram. This is possible when a shell source is activated during its post-AGB evolution. We infer that a late thermal pulse (LTP) experienced by a massive post-AGB star can explain the evolution of the central star. Such a star does not expand significantly as the result of the LTP and does not became a born-again red giant. However, the released energy can remove the tiny H envelope of the star.
Cosmic-ray acceleration at non-relativistic shocks relies on scattering by turbulence that the cosmic rays drive upstream of the shock. We explore the rate of energy transfer from cosmic rays to non-resonant Bell modes and the spectral softening it implies. Accounting for the finite time available for turbulence driving at supernova-remnant shocks yields a smaller spectral impact than found earlier with steady-state considerations. Generally, for diffusion scaling with the Bohm rate by a factor eta, the change in spectral index is at most eta divided by the Alfvenic Mach number of the thermal sub-shock. For M (A) less than or similar to 50 it is well below this limit. Only for very fast shocks and very efficient cosmic-ray acceleration can the change in spectral index reach 0.1. For standard SNR parameters, it is negligible. Independent confirmation is derived by considering the synchrotron energy losses of electrons: if intense nonthermal multi-keV emission is produced, the energy loss, and hence the spectral steepening, is very small for hadronic cosmic rays that produce TeV-band gamma-ray emission.
On the toroidal-velocity antidynamo theorem under the presence of nonuniform electric conductivity
(2022)
Laminar electrically conducting Couette flows with the hydrodynamically stable quasi-Keplerian rotation profile and nonuniform conductivity are probed for dynamo instability. In spherical geometry, the equations for the poloidal and the toroidal field components completely decouple, resulting in free decay, regardless of the spatial distribution of the electric conductivity. In cylindrical geometry the poloidal and toroidal components do not decouple, but here also we do not find dynamo excitations for the cases that the electric conductivity only depends on the radius or - much more complex- that it only depends on the azimuthal or the axial coordinate. The transformation of the plane-flow dynamo model of Busse and Wicht (1992) to cylindrical or spherical geometry therefore fails. It is also shown that even the inclusion of axial flows of both directions does not support the dynamo mechanism. The Elsasser toroidal-velocity antidynamo theorem, according to which dynamos without any radial velocity component cannot work, is thus not softened by nonuniform conductivity distributions.
The investigation of heavy ions dynamics and properties in the Earth's magnetosphere is still an important field of research as they play an important role in several space weather aspects. We present a statistical survey of the average ion mass in the dayside magnetosphere made comparing plasma mass density with electron number density measurements and focusing on both spatial and geomagnetic activity dependence. Field line resonance frequency observations across the European quasi-Meridional Magnetometer Array, are used to infer the equatorial plasma mass density in the range of magnetic L-shells 1.6-6.2. The electron number density is derived from local electric field measurements made on Van Allen Probes using the Neural-network-based Upper-hybrid Resonance Determination algorithm. The analysis is conducted separately for the plasmasphere and the plasmatrough during favorable periods for which both the plasma parameters are observed simultaneously. We found that throughout the plasmasphere the average ion mass is similar or equal to 1 amu for a wide range of geomagnetic activity conditions, suggesting that the plasma mainly consist of hydrogen ions, without regard to the level of geomagnetic activity. Conversely, the plasmatrough is characterized by a variable composition, highlighting a heavy ion mass loading that increases with increasing levels of geomagnetic disturbance. During the most disturbed conditions, the average radial structure shows a broad maximum around 3-4 Earth radii, probably correlated with the accumulation of oxygen ions near the plasmapause. Those ions are mostly observed in the post-dawn and pre-dusk longitudinal sectors.
In recent years, gravitational-wave astronomy has motivated increasingly accurate perturbative studies of gravitational dynamics in compact binaries. This in turn has enabled more detailed analyses of the dynamical black holes in these systems. For example, Pound et al. [Phys. Rev. Lett. 124, 021101 (2020)] recently computed the surface area of a Schwarzschild black hole's apparent horizon, perturbed by an orbiting body, to second order in the binary's mass ratio. In this paper, we take that as the starting point for a comprehensive study of a perturbed Schwarzschild black hole's apparent and event horizon at second perturbative order, deriving generic formulas for the first- and second-order corrections to the horizons' radial profiles, surface areas, Hawking masses, and intrinsic curvatures. We find that the two horizons are remarkably similar, and that any teleological behavior of the event horizon is suppressed in several ways. Critically, we establish that at all orders, the perturbed event horizon in a small-mass-ratio binary is effectively localized in time. Even more pointedly, the event horizon is identical to the apparent horizon at linear order regardless of the source of perturbation, implying that the seemingly teleological "tidal lead," previously observed in linearly perturbed event horizons, is not genuinely teleological in origin. The two horizons do generically differ at second order, but their Hawking masses remain identical, implying that the event horizon obeys the same energy-flux balance law as the apparent horizon. At least in the case of a binary system, the difference between their surface areas remains extremely small even in the late stages of inspiral. In the course of our analysis, we also numerically illustrate puzzling behavior in the black hole's motion around the binary's center of mass.
Uppsala and Berkeley
(2022)
The development of modern photoelectron spectroscopy is reviewed with a special focus on the importance of research at Uppsala University and at Berkeley. The influence of two pioneers, Kai Siegbahn and Dave Shirley, is underlined. Early interaction between the two centers helped to kick-start the field. Both laboratories have continued to play an important role in the field, both in terms of creating new experimental capabilities and developing the theoretical understanding of the spectroscopic processes.
Context.
Even after the Rosetta mission, some of the mechanical parameters of comet 67P/Churyumov-Gerasimenko's surface material are still not well constrained. They are needed to improve our understanding of cometary activity or for planning sample return procedures.
Aims.
We discuss the physical process dominating the formation of aeolian-like surface features in the form of moats and wind taillike bedforms around obstacles and investigate the mechanical and geometrical parameters involved.
Methods.
By applying the discrete element method (DEM) in a low-gravity environment, we numerically simulated the dynamics of the surface layer particles and the particle stream involved in the formation of aeolian-like morphological features. The material is composed of polydisperse spherical particles that consist of a mixture of dust and water ice, with interparticle forces given by the Hertz contact model, cohesion, friction, and rolling friction. We determined a working set of parameters that enables simulations to be reasonably realistic and investigated morphological changes when modifying these parameters.
Results.
The aeolian-like surface features are reasonably well reproduced using model materials with a tensile strength on the order of 0.1-1 Pa. Stronger materials and obstacles with round shapes impede the formation of a moat and a wind tail. The integrated dust flux required for the formation of moats and wind tails is on the order of 100 kg m(-2), which, based on the timescale of morphological changes inferred from Rosetta images, translates to a near-surface particle density on the order of 10(-6)-10(-4) kg m(-3).
Conclusions.
DEM modeling of the aeolian-like surface features reveals complex formation mechanisms that involve both deposition of ejected material and surface erosion. More numerical work and additional in situ measurements or sample return missions are needed to better investigate mechanical parameters of cometary surface material and to understand the mechanics of cometary activity.
Current-day cosmic ray (CR) propagation studies use static Milky Way models and fit parametrized source distributions to data. Instead, we use three-dimensional magnetohydrodynamic (MHD) simulations of isolated galaxies with the moving-mesh code arepo that self-consistently accounts for hydrodynamic effects of CR protons. In post-processing, we calculate their steady-state spectra, taking into account all relevant loss processes. We show that this steady-state assumption is well justified in the disc and generally for regions that emit non-thermal radio and gamma rays. Additionally, we model the spectra of primary electrons, accelerated by supernova remnants, and secondary electrons and positrons produced in hadronic CR proton interactions with the gas. We find that proton spectra above 10 GeV only weakly depend on galactic radius, while they acquire a radial dependence at lower energies due to Coulomb interactions. Radiative losses steepen the spectra of primary CR electrons in the central galactic regions, while diffusive losses dominate in the outskirts. Secondary electrons exhibit a steeper spectrum than primaries because they originate from the transported steeper CR proton spectra. Consistent with Voyager-1 and AMS-02 data, our models (i) show a turnover of proton spectra below GeV energies due to Coulomb interactions so that electrons start to dominate the total particle spectra and (ii) match the shape of the positron fraction up to 10 GeV. We conclude that our steady-state CR modelling in MHD CR galaxy simulations is sufficiently realistic to capture the dominant transport effects shaping their spectra, arguing for a full MHD treatment to accurately model CR transport in the future.
Nanostructured silicon and silicon-aluminum compounds are synthesized by a novel synthesis strategy based on spark plasma sintering (SPS) of silicon nanopowder, mesoporous silicon (pSi), and aluminum nanopowder. The interplay of metal-assisted crystallization and inherent porosity is exploited to largely suppress thermal conductivity. Morphology and temperature-dependent thermal conductivity studies allow us to elucidate the impact of porosity and nanostructure on the macroscopic heat transport. Analytic electron microscopy along with quantitative image analysis is applied to characterize the sample morphology in terms of domain size and interpore distance distributions. We demonstrate that nanostructured domains and high porosity can be maintained in densified mesoporous silicon samples. In contrast, strong grain growth is observed for sintered nanopowders under similar sintering conditions. We observe that aluminum agglomerations induce local grain growth, while aluminum diffusion is observed in porous silicon and dispersed nanoparticles. A detailed analysis of the measured thermal conductivity between 300 and 773 K allows us to distinguish the effect of reduced thermal conductivity caused by porosity from the reduction induced by phonon scattering at nanosized domains. With a modified Landauer/Lundstrom approach the relative thermal conductivity and the scattering length are extracted. The relative thermal conductivity confirms the applicability of Kirkpatrick's effective medium theory. The extracted scattering lengths are in excellent agreement with the harmonic mean of log-normal distributed domain sizes and the interpore distances combined by Matthiessen's rule.
Even after decades of research, the origin of multiple stellar populations in globular clusters (GCs) remains enigmatic. The question as to whether the galaxy environment plays a role in their formation remains unanswered. To that extent, we analysed two classical (>10 Gyr old) Large Magellanic Cloud GCs, NGC 1786 and NGC 1898, using imaging data from Hubble Space Telescope to compare and contrast them with ancient Galactic GCs to assess systematic differences that might exist between their abundance variations. We calculated their red giant branch width, subtracted the effect of metallicity, and compared it with the available data on Galactic GCs by plotting them against initial and current cluster mass. We see that the two clusters follow the same general trend as that of the Galactic GCs, and Galactic GCs from different progenitors follow the same general trend as one another, indicating that the galaxy environment may only play a minor role in the formation of multiple stellar populations within GCs.
We revisited 10 known exoplanetary systems using publicly available data provided by the transiting exoplanet survey satellite (TESS). The sample presented in this work consists of short period transiting exoplanets, with inflated radii and large reported uncertainty on their planetary radii. The precise determination of these values is crucial in order to develop accurate evolutionary models and understand the inflation mechanisms of these systems. Aiming to evaluate the planetary radius measurement, we made use of the planet-to-star radii ratio, a quantity that can be measured during a transit event. We fit the obtained transit light curves of each target with a detrending model and a transit model. Furthermore, we used emcee, which is based on a Markov chain Monte Carlo approach, to assess the best fit posterior distributions of each system parameter of interest. We refined the planetary radius of WASP-140 b by approximately 12%, and we derived a better precision on its reported asymmetric radius uncertainty by approximately 86 and 67%. We also refined the orbital parameters of WASP-120 b by 2 sigma. Moreover, using the high-cadence TESS datasets, we were able to solve a discrepancy in the literature, regarding the planetary radius of the exoplanet WASP-93 b. For all the other exoplanets in our sample, even though there is a tentative trend that planetary radii of (near-) grazing systems have been slightly overestimated in the literature, the planetary radius estimation and the orbital parameters were confirmed with independent observations from space, showing that TESS and ground-based observations are overall in good agreement.
Cloud model inversions of strong chromospheric absorption lines using principal component analysis
(2020)
High-resolution spectroscopy of strong chromospheric absorption lines delivers nowadays several millions of spectra per observing day, when using fast scanning devices to cover large regions on the solar surface. Therefore, fast and robust inversion schemes are needed to explore the large data volume. Cloud model (CM) inversions of the chromospheric H alpha line are commonly employed to investigate various solar features including filaments, prominences, surges, jets, mottles, and (macro-) spicules. The choice of the CM was governed by its intuitive description of complex chromospheric structures as clouds suspended above the solar surface by magnetic fields. This study is based on observations of active region NOAA 11126 in H alpha, which were obtained November 18-23, 2010 with the echelle spectrograph of the vacuum tower telescope at the Observatorio del Teide, Spain. Principal component analysis reduces the dimensionality of spectra and conditions noise-stripped spectra for CM inversions. Modeled H alpha intensity and contrast profiles as well as CM parameters are collected in a database, which facilitates efficient processing of the observed spectra. Physical maps are computed representing the line-core and continuum intensity, absolute contrast, equivalent width, and Doppler velocities, among others. Noise-free spectra expedite the analysis of bisectors. The data processing is evaluated in the context of "big data," in particular with respect to automatic classification of spectra.
Context
Transmission spectroscopy is a promising tool for the atmospheric characterization of transiting exoplanets. Because the planetary signal is faint, discrepancies have been reported regarding individual targets.
Aims
We investigate the dependence of the estimated transmission spectrum on deviations of the orbital parameters of the star-planet system that are due to the limb-darkening effects of the host star. We describe how the uncertainty on the orbital parameters translates into an uncertainty on the planetary spectral slope.
Methods
We created synthetic transit light curves in seven different wavelength bands, from the near-ultraviolet to the near-infrared, and fit them with transit models parameterized by fixed deviating values of the impact parameter b. First, we performed a qualitative study to illustrate the effect by presenting the changes in the transmission spectrum slope with different deviations of b. Then, we quantified these variations by creating an error envelope (for centrally transiting, off-center, and grazing systems) based on a derived typical uncertainty on b from the literature. Finally, we compared the variations in the transmission spectra for different spectral types of host stars.
Results
Our simulations show a wavelength-dependent offset that is more pronounced at the blue wavelengths where the limb-darkening effect is stronger. This offset introduces a slope in the planetary transmission spectrum that becomes steeper with increasing b values. Variations of b by positive or negative values within its uncertainty interval introduce positive or negative slopes, thus the formation of an error envelope. The amplitude from blue optical to near-infrared wavelength for a typical uncertainty on b corresponds to one atmospheric pressure scale height and more. This impact parameter degeneracy is confirmed for different host types; K stars present prominently steeper slopes, while M stars indicate features at the blue wavelengths.
Conclusions
We demonstrate that transmission spectra can be hard to interpret, basically because of the limitations in defining a precise impact parameter value for a transiting exoplanet. This consequently limits a characterization of its atmosphere.
Stars under influence: evidence of tidal interactions between stars and substellar companions
(2023)
Tidal interactions occur between gravitationally bound astrophysical bodies. If their spatial separation is sufficiently small, the bodies can induce tides on each other, leading to angular momentum transfer and altering of evolutionary path the bodies would have followed if they were single objects. The tidal processes are well established in the Solar planet-moon systems and close stellar binary systems. However, how do stars behave if they are orbited by a substellar companion (e.g. a planet or a brown dwarf) on a tight orbit?
Typically, a substellar companion inside the corotation radius of a star will migrate toward the star as it loses orbital angular momentum. On the other hand, the star will gain angular momentum which has the potential to increase its rotation rate. The effect should be more pronounced if the substellar companion is more massive. As the stellar rotation rate and the magnetic activity level are coupled, the star should appear more magnetically active under the tidal influence of the orbiting substellar companion. However, the difficulty in proving that a star has a higher magnetic activity level due to tidal interactions lies in the fact that (I) substellar companions around active stars are easier to detect if they are more massive, leading to a bias toward massive companions around active stars and mimicking the tidal interaction effect, and that (II) the age of a main-sequence star cannot be easily determined, leaving the possibility that a star is more active due to its young age.
In our work, we overcome these issues by employing wide stellar binary systems where one star hosts a substellar companion, and where the other star provides the magnetic activity baseline for the host star, assuming they have coevolved, and thereby provides the host's activity level if tidal interactions have no effect on it. Firstly, we find that extrasolar planets can noticeably increase the host star's X-ray luminosity and that the effect is more pronounced if the exoplanet is at least Jupiter-like in mass and close to the star. Further, we find that a brown dwarf will have an even stronger effect, as expected, and that the X-ray surface flux difference between the host star and the wide stellar companion is a significant outlier when compared to a large sample of similar wide binary systems without any known substellar companions. This result proves that substellar hosting wide binary systems can be good tools to reveal the tidal effect on host stars, and also show that the typical stellar age indicators as activity or rotation cannot be used for these stars. Finally, knowing that the activity difference is a good tracer of the substellar companion's tidal impact, we develop an analytical method to calculate the modified tidal quality factor Q' of individual host stars, which defines the tidal dissipation efficiency in the convective envelope of a given main-sequence star.
The low-mass star GJ 1151 has been reported to display variable low-frequency radio emission, which has been interpreted as a signpost of coronal star-planet interactions with an unseen exoplanet. Here we report the first X-ray detection of GJ 1151's corona based on the XMM-Newton data. We find that the star displays a small flare during the X-ray observation. Averaged over the observation, we detect the star with a low coronal temperature of 1.6 MK and an X-ray luminosity of L-X = 5.5 x 10(26) erg s(-1). During the quiescent time periods excluding the flare, the star remains undetected with an upper limit of L-X,L- qui <= 3.7 x 10(26) erg s(-1). This is compatible with the coronal assumptions used in a recently published model for a star-planet interaction origin of the observed radio signals from this star.
Tidal interaction between an exoplanet and its host star is a possible pathway to transfer angular momentum between the planetary orbit and the stellar spin. In cases where the planetary orbital period is shorter than the stellar rotation period, this may lead to angular momentum being transferred into the star's rotation, possibly counteracting the intrinsic stellar spin-down induced by magnetic braking. Observationally, detecting altered rotational states of single, cool field stars is challenging, as precise ages for such stars are rarely available. Here we present an empirical investigation of the rotation and magnetic activity of a sample of planet-hosting stars that are accompanied by wide stellar companions. Without needing knowledge about the absolute ages of the stars, we test for relative differences in activity and rotation of the planet hosts and their co-eval companions, using X-ray observations to measure the stellar activity levels. Employing three different tidal interaction models, we find that host stars with planets that are expected to tidally interact display elevated activity levels compared to their companion stars. We also find that those activity levels agree with the observed rotational periods for the host stars along the usual rotation-activity relationships, implying that the effect is indeed caused by a tidal interaction and not a purely magnetic interaction that would be expected to affect the stellar activity, but not necessarily the rotation. We conclude that massive, close-in planets have an impact on the stellar rotational evolution, while the smaller, more distant planets do not have a significant influence.
Transit observations in the helium triplet around 10 830 Angstrom are a successful tool to study exoplanetary atmospheres and their mass loss. Forming those lines requires ionization and recombination of helium in the exoplanetary atmosphere. This ionization is caused by stellar photons at extreme ultraviolet (EUV) wavelengths; however, no currently active telescopes can observe this part of the stellar spectrum. The relevant part of the stellar EUV spectrum consists of individual emission lines, many of them being formed by iron at coronal temperatures. The stellar iron abundance in the corona is often observed to be depleted for high-activity low-mass stars due to the first ionization potential (FIP) effect. I show that stars with high versus low coronal iron abundances follow different scaling laws that tie together their X-ray emission and the narrow-band EUV flux that causes helium ionization. I also show that the stellar iron to oxygen abundance ratio in the corona can be measured reasonably well from X-ray CCD spectra, yielding similar results to high-resolution X-ray observations. Taking coronal iron abundance into account, the currently observed large scatter in the relationship of EUV irradiation with exoplanetary helium transit depths can be reduced, improving the target selection criteria for exoplanet transmission spectroscopy. In particular, previously puzzling non-detections of helium for Neptunic exoplanets are now in line with expectations from the revised scaling laws.
Context
Line-driven wind instability is expected to cause small-scale wind inhomogeneities, X-ray emission, and wind line profile variability. The instability can already develop around the sonic point if it is initiated close to the photosphere due to stochastic turbulent motions. In such cases, it may leave its imprint on the light curve as a result of wind blanketing.
Aims
We study the photometric signatures of the line-driven wind instability.
Methods
We used line-driven wind instability simulations to determine the wind variability close to the star. We applied two types of boundary perturbations: a sinusoidal one that enables us to study in detail the development of the instability and a stochastic one given by a Langevin process that provides a more realistic boundary perturbation. We estimated the photometric variability from the resulting mass-flux variations. The variability was simulated assuming that the wind consists of a large number of independent conical wind sectors. We compared the simulated light curves with TESS light curves of OB stars that show stochastic variability.
Results
We find two typical signatures of line-driven wind instability in photometric data: a knee in the power spectrum of magnitude fluctuations, which appears due to engulfment of small-scale structure by larger structures, and a negative skewness of the distribution of fluctuations, which is the result of spatial dominance of rarefied regions. These features endure even when combining the light curves from independent wind sectors.
Conclusions
The stochastic photometric variability of OB stars bears certain signatures of the line-driven wind instability. The distribution function of observed photometric data shows negative skewness and the power spectra of a fraction of light curves exhibit a knee. This can be explained as a result of the line-driven wind instability triggered by stochastic base perturbations.
Meter-sized ground-based telescopes are frequently used today for the follow-up of extrasolar planet candidates. While the transit signal of a Jupiter-sized object can typically be detected to a high level of confidence with small telescope apertures as well, the shallow transit dips of planets with the size of Neptune and smaller are more challenging to reveal. We employ new observational data to illustrate the photometric follow-up capabilities of meter-sized telescopes for shallow exoplanet transits. We describe in detail the capability of distinguishing the photometric signal of an exoplanet transit from an underlying trend in the light curve. The transit depths of the six targets we observed, Kepler-94b, Kepler-63b, K2-100b, K2-138b, K2-138c, and K2-138e, range from 3.9 ppt down to 0.3 ppt. For five targets of this sample, we provide the first ground-based photometric follow-up. The timing of three targets is precisely known from previous observations, and the timing of the other three targets is uncertain and we aim to constrain it. We detect or rule out the transit features significantly in single observations for the targets that show transits of 1.3 ppt or deeper. The shallower transit depths of two targets of 0.6 and 0.8 ppt were detected tentatively in single light curves, and were detected significantly by repeated observations. Only for the target of the shallowest transit depth of 0.3 ppt were we unable to draw a significant conclusion despite combining five individual light curves. An injection-recovery test on our real data shows that we detect transits of 1.3 ppt depth significantly in single light curves if the transit is fully covered, including out-of-transit data toward both sides, in some cases down to 0.7 ppt depth. For Kepler-94b, Kepler-63b, and K2-100b, we were able to verify the ephemeris. In the case of K2-138c with a 0.6 ppt deep transit, we were able to refine it, and in the case of K2-138e, we ruled out the transit in the time interval of more than ±1.5 σ of its current literature ephemeris.
We present observations of three-dimensional magnetic power spectra in wavevector space to investigate the anisotropy and scalings of sub-Alfvenic solar wind turbulence at magnetohydrodynamic (MHD) scale using the Magnetospheric Multiscale spacecraft. The magnetic power distributions are organized in a new coordinate determined by wavevectors ((kappa) over cap) and background magnetic field ((b) over cap (0)) in Fourier space. This study utilizes two approaches to determine wavevectors: the singular value decomposition method and multispacecraft timing analysis. The combination of the two methods allows an examination of the properties of magnetic field fluctuations in terms of mode compositions without any spatiotemporal hypothesis. Observations show that fluctuations (delta B-perpendicular to 1) in the direction perpendicular to (kappa) over cap and (b) over cap (0) prominently cascade perpendicular to (b) over cap (0), and such anisotropy increases with wavenumbers. The reduced power spectra of 6.8 11 follow Goldreich-Sridhar scalings: (P) over cap (k(perpendicular to)) proportional to k(perpendicular to)(-5/3) and (P) over cap (k(parallel to)) proportional to k(parallel to)(-2). In contrast, fluctuations within the (k) over cap(b) over cap (0) plane show isotropic behaviors: perpendicular power distributions are approximately the same as parallel distributions. The reduced power spectra of fluctuations within the (k) over cap(b) over cap (0) plane follow the scalings (P) over cap (k(perpendicular to)) proportional to k(perpendicular to)(-3/2) and (P) over cap (k(parallel to)) proportional to k(parallel to)(-3/2). Comparing frequency-wavevector spectra with theoretical dispersion relations of MHD modes, we find that delta B-perpendicular to 1 are probably associated with Alfven modes. On the other hand, magnetic field fluctuations within the (k) over cap(b) over cap (0) plane more likely originate from fast modes based on their isotropic behaviors. The observations of anisotropy and scalings of different magnetic field components are consistent with the predictions of current compressible MHD theory. Moreover, for the Alfvenic component, the ratio of cascading time to the wave period is found to be a factor of a few, consistent with critical balance in the strong turbulence regime. These results are valuable for further studies of energy compositions of plasma turbulence and their effects on energetic particle transport.
Late-type stars are by far the most frequent stars in the universe and of fundamental interest to various fields of astronomy – most notably to Galactic archaeology and exoplanet research. However, such stars barely change during their main sequence lifetime; their temperature, luminosity, or chemical composition evolve only very slowly over the course of billions of years. As such, it is difficult to obtain the age of such a star, especially when it is isolated and no other indications (like cluster association) can be used. Gyrochronology offers a way to overcome this problem.
Stars, just like all other objects in the universe, rotate and the rate at which stars rotate impacts many aspects of their appearance and evolution. Gyrochronology leverages the observed rotation rate of a late-type main sequence star and its systematic evolution to estimate their ages. Unlike the above-mentioned parameters, the rotation rate of a main sequence star changes drastically throughout its main sequence lifetime; stars spin down. The youngest stars rotate every few hours, whereas much older stars rotate only about once a month, or – in the case of some late M-stars – once in a hundred days. Given that this spindown is systematic (with an additional mass dependence), it gave rise to the idea of using the observed rotation rate of a star (and its mass or a suitable proxy thereof) to estimate a star’s age. This has been explored widely in young stellar open clusters but remains essentially unconstrained for stars older than the sun, and K and M stars older than 1 Gyr.
This thesis focuses on the continued exploration of the spindown behavior to assess, whether gyrochronology remains applicable for stars of old ages, whether it is universal for late-type main sequence stars (including field stars), and to provide calibration mileposts for spindown models. To accomplish this, I have analyzed data from Kepler space telescope for the open clusters Ruprecht 147 (2.7 Gyr old) and M 67 (4 Gyr). Time series photometry data (light curves)
were obtained for both clusters during Kepler’s K2 mission. However, due to technical limitations and telescope malfunctions, extracting usable data from the K2 mission to identify (especially long) rotation periods requires extensive data preparation.
For Ruprecht 147, I have compiled a list of about 300 cluster members from the literature and adopted preprocessed light curves from the Kepler archive where available. They have been cleaned of the gravest of data artifacts but still contained systematics. After correcting them for said artifacts, I was able to identify rotation periods in 31 of them.
For M 67 more effort was taken. My work on Ruprecht 147 has shown the limitations imposed by the preselection of Kepler targets. Therefore, I adopted the time series full frame image directly and performed photometry on a much higher spatial resolution to be able to obtain data for as many stars as possible. This also means that I had to deal with the ubiquitous artifacts in Kepler data. For that, I devised a method that correlates the artificial flux variations with the ongoing drift of the telescope pointing in order to remove it. This process was a large success and I was able to create light curves whose quality match and even exceede those that were created by the Kepler mission – all while operating on higher spatial resolution and processing fainter stars. Ultimately, I was able to identify signs of periodic variability in the (created) light curves for 31 and 47 stars in Ruprecht 147 and M 67, respectively. My data connect well to bluer stars of cluster of the same age and extend for the first time to stars redder than early-K and older than 1 Gyr. The cluster data show a clear flattening in the distribution of Ruprecht 147 and even a downturn for M 67, resulting in a somewhat sinusoidal shape. With that, I have shown that the systematic spindown of stars continues at least until 4 Gyr and stars continue to live on a single surface in age-rotation periods-mass space which allows gyrochronology to be used at least up to that age. However, the shape of the spindown – as exemplified by the newly discovered sinusoidal shape of the cluster sequence – deviates strongly from the expectations.
I then compiled an extensive sample of rotation data in open clusters – very much including my own work – and used the resulting cluster skeleton (with each cluster forming a rip in color-rotation period-mass space) to investigate if field stars follow the same spindown as cluster stars. For the field stars, I used wide binaries, which – with their shared origin and coevality – are in a sense the smallest possible open clusters. I devised an empirical method to evaluate the consistency between the rotation rates of the wide binary components and found that the vast majority of them are in fact consistent with what is observed in open clusters. This leads me to conclude that gyrochronology – calibrated on open clusters – can be applied to determine the ages of field stars.
We quantitatively address the conjecture that magnetic helicity must be shed from the Sun by eruptions launching coronal mass ejections in order to limit its accumulation in each hemisphere. By varying the ratio of guide and strapping field and the flux rope twist in a parametric simulation study of flux rope ejection from approximately marginally stable force-free equilibria, different ratios of self- and mutual helicity are set and the onset of the torus or helical kink instability is obtained. The helicity shed is found to vary over a broad range from a minor to a major part of the initial helicity, with self helicity being largely or completely shed and mutual helicity, which makes up the larger part of the initial helicity, being shed only partly. Torus-unstable configurations with subcritical twist and without a guide field shed up to about two-thirds of the initial helicity, while a highly twisted, kink-unstable configuration sheds only about one-quarter. The parametric study also yields stable force-free flux rope equilibria up to a total flux-normalized helicity of 0.25, with a ratio of self- to total helicity of 0.32 and a ratio of flux rope to external poloidal flux of 0.94. These results numerically demonstrate the conjecture of helicity shedding by coronal mass ejections and provide a first account of its parametric dependence. Both self- and mutual helicity are shed significantly; this reduces the total initial helicity by a fraction of ∼0.4--0.65 for typical source region parameters.
High-energy irradiation is a driver for atmospheric evaporation and mass loss in exoplanets. This work is based on data from eROSITA, the soft X-ray instrument on board the Spectrum Roentgen Gamma mission, as well as on archival data from other missions. We aim to characterise the high-energy environment of known exoplanets and estimate their mass-loss rates. We use X-ray source catalogues from eROSITA, XMM-Newton, Chandra, and ROSAT to derive X-ray luminosities of exoplanet host stars in the 0.2–2 keV energy band with an underlying coronal, that is, optically thin thermal spectrum. We present a catalogue of stellar X-ray and EUV luminosities, exoplanetary X-ray and EUV irradiation fluxes, and estimated mass-loss rates for a total of 287 exoplanets, 96 of which are characterised for the first time based on new eROSITA detections. We identify 14 first-time X-ray detections of transiting exoplanets that are subject to irradiation levels known to cause observable evaporation signatures in other exoplanets. This makes them suitable targets for follow-up observations.
Mergers play an important role in galaxy evolution. In particular, major mergers are able to have a transformative effect on galaxy morphology. In this paper, we investigate the role of magnetic fields in gas-rich major mergers. To this end, we run a series of high-resolution magnetohydrodynamic (MHD) zoom-in simulations with the moving-mesh code arepo and compare the outcome with hydrodynamic simulations run from the same initial conditions. This is the first time that the effect of magnetic fields in major mergers has been investigated in a cosmologically consistent manner. In contrast to previous non-cosmological simulations, we find that the inclusion of magnetic fields has a substantial impact on the production of the merger remnant. Whilst magnetic fields do not strongly affect global properties, such as the star formation history, they are able to significantly influence structural properties. Indeed, MHD simulations consistently form remnants with extended discs and well-developed spiral structure, whilst hydrodynamic simulations form more compact remnants that display distinctive ring morphology. We support this work with a resolution study and show that whilst global properties are broadly converged across resolution and physics models, morphological differences only develop given sufficient resolution. We argue that this is due to the more efficient excitement of a small-scale dynamo in higher resolution simulations, resulting in a more strongly amplified field that is better able to influence gas dynamics.
A very small fraction of (runaway) massive stars have masses exceeding 60-70 M-circle dot and are predicted to evolve as luminous blue variable and Wolf-Rayet stars before ending their lives as core-collapse supernovae. Our 2D axisymmetric hydrodynamical simulations explore how a fast wind (2000 km s(-1)) and high mass-loss rate (10(-5)M(circle dot) yr(-1)) can impact the morphology of the circumstellar medium. It is shaped as 100 pc-scale wind nebula that can be pierced by the driving star when it supersonically moves with velocity 20-40 km s(-1) through the interstellar medium (ISM) in the Galactic plane. The motion of such runaway stars displaces the position of the supernova explosion out of their bow shock nebula, imposing asymmetries to the eventual shock wave expansion and engendering Cygnus-loop-like supernova remnants. We conclude that the size (up to more than 200 pc) of the filamentary wind cavity in which the chemically enriched supernova ejecta expand, mixing efficiently the wind and ISM materials by at least 10 per cent in number density, can be used as a tracer of the runaway nature of the very massive progenitors of such 0.1Myr old remnants. Our results motivate further observational campaigns devoted to the bow shock of the very massive stars BD+43 degrees 3654 and to the close surroundings of the synchrotron-emitting Wolf-Rayet shell G2.4+1.4.
Context.
Tycho's supernova remnant (SNR) is associated with the historical supernova (SN) event SN 1572 of Type Ia. The explosion occurred in a relatively clean environment, and was visually observed, providing an age estimate. This SNR therefore represents an ideal astrophysical test-bed for the study of cosmic-ray acceleration and related phenomena. A number of studies suggest that shock acceleration with particle feedback and very efficient magnetic-field amplification combined with Alfvenic drift are needed to explain the rather soft radio spectrum and the narrow rims observed in X-rays.
Aims.
We show that the broadband spectrum of Tycho's SNR can alternatively be well explained when accounting for stochastic acceleration as a secondary process. The re-acceleration of particles in the turbulent region immediately downstream of the shock should be efficient enough to impact particle spectra over several decades in energy. The so-called Alfvenic drift and particle feedback on the shock structure are not required in this scenario. Additionally, we investigate whether synchrotron losses or magnetic-field damping play a more profound role in the formation of the non-thermal filaments.
Methods.
We solved the full particle transport equation in test-particle mode using hydrodynamic simulations of the SNR plasma flow. The background magnetic field was either computed from the induction equation or follows analytic profiles, depending on the model considered. Fast-mode waves in the downstream region provide the diffusion of particles in momentum space.
Results.
We show that the broadband spectrum of Tycho can be well explained if magnetic-field damping and stochastic re-acceleration of particles are taken into account. Although not as efficient as standard diffusive shock acceleration, stochastic acceleration leaves its imprint on the particle spectra, which is especially notable in the emission at radio wavelengths. We find a lower limit for the post-shock magnetic-field strength similar to 330 mu G, implying efficient amplification even for the magnetic-field damping scenario. Magnetic-field damping is necessary for the formation of the filaments in the radio range, while the X-ray filaments are shaped by both the synchrotron losses and magnetic-field damping.
We present the results from Chandra X-ray observations, and near- and mid-infrared analysis, using VISTA/VVV and Spitzer/GLIMPSE catalogs, of the high-mass star-forming region IRAS 16562-3959, which contains a candidate for a high-mass protostar. We detected 249 X-ray sources within the ACIS-I field of view. The majority of the X-ray sources have low count rates (<0.638 cts/ks) and hard X-ray spectra. The search for YSOs in the region using VISTA/VVV and Spitzer/GLIMPSE catalogs resulted in a total of 636 YSOs, with 74 Class I and 562 Class II YSOs. The search for near- and mid-infrared counterparts of the X-ray sources led to a total of 165 VISTA/VVV counterparts, and a total of 151 Spitzer/GLIMPSE counterparts. The infrared analysis of the X-ray counterparts allowed us to identify an extra 91 Class III YSOs associated with the region. We conclude that a total of 727 YSOs are associated with the region, with 74 Class I, 562 Class II, and 91 Class III YSOs. We also found that the region is composed of 16 subclusters. In the vicinity of the high-mass protostar, the stellar distribution has a core-halo structure. The subcluster containing the high-mass protostar is the densest and the youngest in the region, and the high-mass protostar is located at its center. The YSOs in this cluster appear to be substantially older than the high-mass protostar.
Planets around young stars are thought to undergo atmospheric evaporation due to the high magnetic activity of the host stars. Here we report on X-ray observations of V1298 Tau, a young star with four transiting exoplanets. We use X-ray observations of the host star with Chandra and ROSAT to measure the current high-energy irradiation level of the planets and employ a model for the stellar activity evolution together with exoplanetary mass-loss to estimate the possible evolution of the planets. We find that V1298 Tau is X-ray bright with log L-X [erg s(-1)] = 30.1 and has a mean coronal temperature of approximate to 9 MK. This places the star amongst the more X-ray luminous ones at this stellar age. We estimate the radiation-driven mass-loss of the exoplanets and find that it depends sensitively on the possible evolutionary spin-down tracks of the star as well as on the current planetary densities. Assuming the planets are of low density due to their youth, we find that the innermost two planets can lose significant parts of their gaseous envelopes and could be evaporated down to their rocky cores depending on the stellar spin evolution. However, if the planets are heavier and follow the mass-radius relation of older planets, then even in the highest XUV irradiation scenario none of the planets is expected to cross the radius gap into the rocky regime until the system reaches an age of 5 Gyr.
A second peak in the extreme ultraviolet sometimes appears during the gradual phase of solar flares, which is known as the EUV late phase (ELP). Stereotypically ELP is associated with two separated sets of flaring loops with distinct sizes, and it has been debated whether ELP is caused by additional heating or extended plasma cooling in the longer loop system. Here we carry out a survey of 55 M-and-above GOES-class flares with ELP during 2010-2014. Based on the flare-ribbon morphology, these flares are categorized as circular-ribbon (19 events), two-ribbon (23 events), and complex-ribbon (13 events) flares. Among them, 22 events (40%) are associated with coronal mass ejections, while the rest are confined. An extreme ELP, with the late-phase peak exceeding the main-phase peak, is found in 48% of two-ribbon flares, 37% of circular-ribbon flares, and 31% of complex-ribbon flares, suggesting that additional heating is more likely present during ELP in two-ribbon than in circular-ribbon flares. Overall, cooling may be the dominant factor causing the delay of the ELP peak relative to the main-phase peak, because the loop system responsible for the ELP emission is generally larger than, and well separated from, that responsible for the main-phase emission. All but one of the circular-ribbon flares can be well explained by a composite "dome-plate" quasi-separatrix layer (QSL). Only half of these show a magnetic null point, with its fan and spine embedded in the dome and plate, respectively. The dome-plate QSL, therefore, is a general and robust structure characterizing circular-ribbon flares.
We study populations of globally coupled noisy rotators (oscillators with inertia) allowing a nonequilibrium transition from a desynchronized state to a synchronous one (with the nonvanishing order parameter). The newly developed analytical approaches resulted in solutions describing the synchronous state with constant order parameter for weakly inertial rotators, including the case of zero inertia, when the model is reduced to the Kuramoto model of coupled noise oscillators. These approaches provide also analytical criteria distinguishing supercritical and subcritical transitions to the desynchronized state and indicate the universality of such transitions in rotator ensembles. All the obtained analytical results are confirmed by the numerical ones, both by direct simulations of the large ensembles and by solution of the associated Fokker-Planck equation. We also propose generalizations of the developed approaches for setups where different rotators parameters (natural frequencies, masses, noise intensities, strengths and phase shifts in coupling) are dispersed.
Context.
Wide hot subdwarf B (sdB) binaries with main-sequence companions are outcomes of stable mass transfer from evolved red giants. The orbits of these binaries show a strong correlation between their orbital periods and mass ratios. The origins of this correlation have, so far, been lacking a conclusive explanation.
Aims.
We aim to find a binary evolution model which can explain the observed correlation.
Methods.
Radii of evolved red giants, and hence the resulting orbital periods, strongly depend on their metallicity. We performed a small but statistically significant binary population synthesis study with the binary stellar evolution code MESA. We used a standard model for binary mass loss and a standard metallicity history of the Galaxy. The resulting sdB systems were selected based on the same criteria as was used in observations and then compared with the observed population.
Results.
We have achieved an excellent match to the observed period-mass ratio correlation without explicitly fine-tuning any parameters. Furthermore, our models produce a very good match to the observed period-metallicity correlation. We predict several new correlations, which link the observed sdB binaries to their progenitors, and a correlation between the orbital period, metallicity, and core mass for subdwarfs and young low-mass helium white dwarfs. We also predict that sdB binaries have distinct orbital properties depending on whether they formed in the Galactic bulge, thin or thick disc, or the halo.
Conclusions
We demonstrate, for the first time, how the metallicity history of the Milky Way is imprinted in the properties of the observed post-mass transfer binaries. We show that Galactic chemical evolution is an important factor in binary population studies of interacting systems containing at least one evolved low-mass (M-init< 1.6 M-circle dot) component. Finally, we provide an observationally supported model of mass transfer from low-mass red giants onto main-sequence stars.
Undifferentiated meteorites, like primitive chondrites, can contain presolar and solar nebula materials which would provide information about the origin and initial conditions of the solar system, whereas differentiated meteorites like iron meteorites, can show early phases of planetary accretion. They also provide the possibility to receive information about core properties and planetary bodies. In addition to the gain in such fundamental scientific knowledge both types are of interest for the exploration of critical raw materials (CRMs) and precious elements.
The Santa Rosa de Viterbo meteorite shower, discovered 1810 in the Boyaca province of Colombia, represents a typical iron-nickel meteorite. The present study presents new structural, textural and geochemical results of one fragment of this meteorite, using reflecting microscopy, electron probe micro analyses (EPMA) and electron microscopy (SEM) with energy dispersive X-ray spectroscopy (EDX). The present study presents trace element concentrations of the meteorite's minerals for the first time.
The sample is dominated by kamacite (alpha-FeNi). Schreibersite (FeNi3P), taenite (gamma-FeNi) and plessite (mixture of kamacite and taenite) are minor constituents. The occurrence of cohenite ((Fe,Ni,Co)(3)C) and troilite (FeS) are likely. The meteorite sample contains classical Neuman bands passing through kamacite and frequent Widmanstadtten pattern. The bandwidth of kamacite defines the meteorite as finest octahedrite. Geochemically, it is characterized as a "Type IC meteorite".
While improving the characterization and classification of the Santa Rosa de Viterbo Iron Meteorite, notable concentrations of Au (>400 ppm) and Ge (>230 ppm) alongside major elements such as Fe, Ni and Co in the bulk composition of that meteorite, were proven. Major and rock-forming minerals such as kamacite and taenite incorporate hundreds of ppm of Ge whereas schreibersite, itself a minor component in that particular meteorite, is the major source for Au (>1400 ppm). In kamacite and taenite also Ir, Pd and Ga were found in minor amounts. Nano-scale inclusions or atomic clusters called nano-nuggets may have been responsible for the high concentrations of Au, Ir, Pd and Ga. Raman and Laser-induced plasma spectroscopes installed in in space probes seems suitable exploration methods for Fe-Ni meteorites, containing Ni-concentrations > 5.8 wt% defining the meteorite as octaedrites.
Deep observation of the High Energy Stereoscopic System (HESS) on the most extended pulsar wind nebula HESS J1825-137 reveals an enhanced energy-dependent morphology, providing useful information on the particle transport mechanism in the nebula. We find that the energy-dependent morphology is consistent with a diffusion-dominated transport of electrons/positrons. It provides an alternative possible interpretation for the unusually large spatial extent (i.e. greater than or similar to 100 pc) of the nebula, which could then be attributed to the diffusion of escaping electrons/positrons from a compact plerion. The influence of various model parameters on the energy-dependent extent of the nebula is studied in the diffusion-dominated scenario. We also show that the energy-dependent morphology of the nebula may also be used to study the spin-down history of the pulsar.
Numerical-relativity simulations are essential for studying the last stages of the binary neutron star coalescence. Unfortunately, for stable simulations there is the need to add an artificial low-density atmosphere. Here we discuss a new framework in which we can effectively set the density surrounding the neutron stars to zero to ensure a more accurate simulation. We test our method with a number of single star test cases and for an equal-mass binary neutron star simulation. While the bulk motion of the system is not influenced, and hence there is no improvement with respect to the emitted gravitational-wave signal, we find that the new approach is superior with respect to mass conservation and it allows a much better tracking of outward moving material. This will allow a more accurate simulation of the ejected material and supports the interpretation of present and future multimessenger observations with more accurate numerical-relativity simulations.
Context.
Blue horizontal-branch stars evolve from low-mass stars that have completed their main-sequence lifetimes and undergone a helium flash at the end of their red-giant phase. As such, blue horizontal-branch stars are very old objects that can be used as markers in studies of the Galactic structure and formation history. To create a clean sky catalogue of blue horizontal-branch stars, we cross-matched the Gaia data release 2 (DR2) dataset with existing reference catalogues to define selection criteria based on Gaia DR2 parameters. Following the publication of Gaia early data release 3 (EDR3), these methods were verified and subsequently applied to this latest release.
Aims.
Previous catalogues of blue horizontal-branch stars were developed using spectral analyses or were restricted to individual globular clusters. The purpose of this catalogue is to identify a set of blue horizontal-branch star candidates that have been selected using photometric and astrometric observations and exhibits a low contamination rate. This has been deemed important as the success of the Gaia mission has changed the way that targets are selected for large-scale spectroscopic surveys, meaning that far fewer spectra will be acquired for blue horizontal-branch stars in the future unless they are specifically targeted. <br /> Methods. We cross-matched reference blue horizontal-branch datasets with the Gaia DR2 database and defined two sets of selection criteria. Firstly, in Gaia DR2 - colour and absolute G magnitude space, and secondly, in Gaia DR2 - colour and reduced proper motion space. The main-sequence contamination in both subsets of the catalogue was reduced, at the expense of completeness, by concentrating on the Milky Way's Galactic halo, where relatively young main-sequence stars were not expected. The entire catalogue is limited to those stars with no apparent neighbours within 5 arcsec. These methods were verified and subsequently applied to the Gaia EDR3.
Results.
We present a catalogue, based on Gaia EDR3, of 57 377 blue horizontal-branch stars. The Gaia EDR3 parallax was used in selecting 16 794 candidates and the proper motions were used to identify a further 40 583 candidates.
Planets outside our solar system, so-called "exoplanets", can be detected with different methods, and currently more than 5000 exoplanets have been confirmed, according to NASA Exoplanet Archive. One major highlight of the studies on exoplanets in the past twenty years is the characterization of their atmospheres usingtransmission spectroscopy as the exoplanet transits. However, this characterization is a challenging process and sometimes there are reported discrepancies in the literature regarding the atmosphere of the same exoplanet. One potential reason for the observed atmospheric inconsistencies is called impact parameter degeneracy, and it is highly driven by the limb darkening effect of the host star. A brief introductionto those topics in presented in chapter 1, while the motivation and objectives of thiswork are described in chapter 2.The first goal is to clarify the origin of the transmission spectrum, which is anindicator of an exoplanet’s atmosphere; whether it is real or influenced by the impactparameter degeneracy. A second goal is to determine whether photometry from space using the Transiting Exoplanet Survey Satellite (TESS), could improve on the major parameters, which are responsible for the aforementioned degeneracy, of known exoplanetary systems. Three individual projects were conducted in order toaddress those goals. The three manuscripts are presented, in short, in the manuscriptoverview in chapter 3.More specifically, in chapter 4, the first manuscript is presented, which is an ex-tended investigation on the impact parameter degeneracy and its application onsynthetic transmission spectra. Evidently, the limb darkening of the host star isan important driver for this effect. It keeps the degeneracy persisting through different groups of exoplanets, based on the uncertainty of their impact parameter and on the type of their host star. The second goal, was addressed in the second and third manuscripts (chapter 5 and chapter 6 respectively). Using observationsfrom the TESS mission, two samples of exoplanets were studied; 10 transiting inflated hot-Jupiters and 43 transiting grazing systems. Potentially, the refinement or confirmation of their major system parameters’ measurements can assist in solving current or future discrepancies regarding their atmospheric characterization.In chapter 7 the conclusions of this work are discussed, while in chapter 8 itis proposed how TESS’s measurements can be able to discern between erroneousinterpretations of transmission spectra, especially on systems where the impact parameter degeneracy is likely not applicable.
Unveiling the Local Universe
(2023)
In order to explain the variance of the solar rotation law during the activity minima and maxima, the angular momentum transport by rotating magnetoconvection is simulated in a convective box penetrated by an inclined azimuthal magnetic field. Turbulence-induced kinetic and magnetic stresses and the Lorentz force of the large-scale magnetic background field are the basic transporters of angular momentum. Without rotation, the sign of the magnetic stresses naturally depends on the signs of the field components as positive (negative) B theta B phi transport the angular momentum poleward (equatorward). For fast enough rotation, however, the turbulence-originated Reynolds stresses start to dominate the transport of the angular momentum flux. The simulations show that positive ratios of the two meridional magnetic field components to the azimuthal field reduce the inward radial as well as the equatorward latitudinal transport, which result from hydrodynamic calculations. Only for B theta B phi>0 (generated by solar-type rotation laws with an accelerated equator) does the magnetic-influenced rotation at the solar surface prove to be flatter than the nonmagnetic profile together with the observed slight spin-down of the equator. The latter phenomenon does not appear for antisolar rotation with polar vortex as well as for rotation laws with prevailing radial shear.
Discrete element modeling of boulder and cliff morphologies on comet 67P/Churyumov-Gerasimenko
(2020)
Context:
Even after the Rosetta mission, some of the mechanical parameters of comet 67P/Churyumov-Gerasimenko's surface material are not yet well constrained. These parameters are needed to improve our understanding of cometary activity or for planning sample return missions.
Aims:
We study some of the physical processes involved in the formation of selected surface features and investigate the mechanical and geometrical parameters involved.
Methods:
Applying the discrete element method (DEM) in a low-gravity environment, we numerically simulated the surface layer particle dynamics involved in the formation of selected morphological features. The material considered is a mixture of polydisperse ice and dust spheres with inter-particle forces given by the Hertz contact model, translational friction, rolling friction, cohesion from unsintered contacts, and optionally due to bonds from ice sintering. We determined a working set of parameters that enables the simulations to be reasonably realistic and investigated morphological changes due to modifications thereof.
Results:
The selected morphological features are reasonably well reproduced using model materials with a tensile strength on the order of 1-10 Pa. Increasing the diameters of the spherical particles decreases the material strength, and increasing the friction leads to a more brittle but somewhat stronger material. High friction is required to make the material sufficiently brittle to match observations, which points to the presence of very rough, even angular particles. Reasonable seismic activity does not suffice to trigger the collapses of cliffs without material heterogeneities or structural defects.
Conclusions:
DEM modeling can be a powerful tool to investigate mechanical parameters of cometary surface material. However, many uncertainties arise from our limited understanding of particle shapes, spatial configurations, and size distributions, all on multiple length scales. Further numerical work, in situ measurements, and sample return missions are needed to better understand the mechanics of cometary material and cometary activity.