520 Astronomie und zugeordnete Wissenschaften
Refine
Year of publication
Document Type
- Article (324)
- Doctoral Thesis (88)
- Conference Proceeding (66)
- Postprint (23)
- Other (7)
- Review (3)
- Habilitation Thesis (2)
- Master's Thesis (1)
Keywords
- stars: massive (19)
- cosmic rays (17)
- methods: numerical (17)
- stars: atmospheres (15)
- stars: evolution (15)
- stars: mass-loss (14)
- stars: winds, outflows (13)
- stars: early-type (12)
- subdwarfs (12)
- MHD (11)
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.
Context:
Gyrochronology allows the derivation of ages for cool main sequence stars based on their observed rotation periods and masses, or a suitable proxy thereof. It is increasingly well-explored for FGK stars, but requires further measurements for older ages and K - M-type stars.
Aims:
We study the 2.7 Gyr-old open cluster Ruprecht 147 to compare it with the previously-studied, but far more distant, NGC 6819 cluster, and especially to measure cooler stars than was previously possible there.
Methods:
We constructed an inclusive list of 102 cluster members from prior work, including Gaia DR2, and for which light curves were also obtained during Campaign 7 of the Kepler/K2 space mission. We placed them in the cluster color-magnitude diagram and checked the related information against appropriate isochrones. The light curves were then corrected for data systematics using Principal Component Analysis on all observed K2 C07 stars and subsequently subjected to periodicity analysis.
Results:
Periodic signals are found for 32 stars, 21 of which are considered to be both highly reliable and to represent single, or effectively single, Ru 147 stars. These stars cover the spectral types from late-F to mid-M stars, and they have periods ranging from 6 d - 33 d, allowing for a comparison of Ruprecht 147 to both other open clusters and to models of rotational spindown. The derived rotation periods connect reasonably to, overlap with, and extend to lower masses the known rotation period distribution of the 2.5 Gyr-old cluster NGC 6819.
Conclusions:
The data confirm that cool stars lie on a single surface in rotation period-mass-age space, and they simultaneously challenge its commonly assumed shape. The shape at the low mass region of the color-period diagram at the age of Ru 147 favors a recently-proposed model which requires a third mass-dependent timescale in addition to the two timescales required by a former model, suggesting that a third physical process is required to model rotating stars effectively.
Stellar atmosphere modeling and chemical abundance determinations require the knowledge of spectral line shapes. Spectral lines of chromium in various ionization stages are common in stellar spectra but detailed data on Stark broadening for them is scarce. Recently we reported on the first calculations of Stark widths for several 4s-4p transitions of double-ionized chromium, employing the Modified Semi-Empirical approach (MSE). In this work we present applications of the data to spectrum synthesis of Cr III lines in the ultraviolet region of hot stars. The Atlas9 model atmosphere code and the line-formation code Surface were used with the assumption of local thermodynamic equilibrium. The improvements of adopting the MSE broadening tables instead of approximate Stark broadening coefficients are investigated for a total of 56 Cr III lines visible in HST/STIS spectra of the B3 subgiant star Iota Herculis and the subdwarf B-star Feige 66.
In recent decades, astronomy has seen a boom in large-scale stellar surveys of the Galaxy. The detailed information obtained about millions of individual stars in the Milky Way is bringing us a step closer to answering one of the most outstanding questions in astrophysics: how do galaxies form and evolve? The Milky Way is the only galaxy where we can dissect many stars into their high-dimensional chemical composition and complete phase space, which analogously as fossil records can unveil the past history of the genesis of the Galaxy. The processes that lead to large structure formation, such as the Milky Way, are critical for constraining cosmological models; we call this line of study Galactic archaeology or near-field cosmology.
At the core of this work, we present a collection of efforts to chemically and dynamically characterise the disks and bulge of our Galaxy. The results we present in this thesis have only been possible thanks to the advent of the Gaia astrometric satellite, which has revolutionised the field of Galactic archaeology by precisely measuring the positions, parallax distances and motions of more than a billion stars. Another, though not less important, breakthrough is the APOGEE survey, which has observed spectra in the near-infrared peering into the dusty regions of the Galaxy, allowing us to determine detailed chemical abundance patterns in hundreds of thousands of stars. To accurately depict the Milky Way structure, we use and develop the Bayesian isochrone fitting tool/code called StarHorse; this software can predict stellar distances, extinctions and ages by combining astrometry, photometry and spectroscopy based on stellar evolutionary models. The StarHorse code is pivotal to calculating distances where Gaia parallaxes alone cannot allow accurate estimates.
We show that by combining Gaia, APOGEE, photometric surveys and using StarHorse, we can produce a chemical cartography of the Milky way disks from their outermost to innermost parts. Such a map is unprecedented in the inner Galaxy. It reveals a continuity of the bimodal chemical pattern previously detected in the solar neighbourhood, indicating two populations with distinct formation histories. Furthermore, the data reveals a chemical gradient within the thin disk where the content of 𝛼-process elements and metals is higher towards the centre. Focusing on a sample in the inner MW we confirm the extension of the chemical duality to the innermost regions of the Galaxy. We find stars with bar shape orbits to show both high- and low-𝛼 abundances, suggesting the bar formed by secular evolution trapping stars that already existed. By analysing the chemical orbital space of the inner Galactic regions, we disentangle the multiple populations that inhabit this complex region. We reveal the presence of the thin disk, thick disk, bar, and a counter-rotating population, which resembles the outcome of a perturbed proto-Galactic disk. Our study also finds that the inner Galaxy holds a high quantity of super metal-rich stars up to three times solar suggesting it is a possible repository of old super-metal-rich stars found in the solar neighbourhood.
We also enter into the complicated task of deriving individual stellar ages. With StarHorse, we calculate the ages of main-sequence turn-off and sub-giant stars for several public spectroscopic surveys. We validate our results by investigating linear relations between chemical abundances and time since the 𝛼 and neutron capture elements are sensitive to age as a reflection of the different enrichment timescales of these elements. For further study of the disks in the solar neighbourhood, we use an unsupervised machine learning algorithm to delineate a multidimensional separation of chrono-chemical stellar groups revealing the chemical thick disk, the thin disk, and young 𝛼-rich stars. The thick disk is shown to have a small age dispersion indicating its fast formation contrary to the thin disk that spans a wide range of ages.
With groundbreaking data, this thesis encloses a detailed chemo-dynamical view of the disk and bulge of our Galaxy. Our findings on the Milky Way can be linked to the evolution of high redshift disk galaxies, helping to solve the conundrum of galaxy formation.
The evolution of a galaxy is pivotally governed by its pattern of star formation over a given period of time. The star formation rate at any given time is strongly dependent on the amount of cold gas available in the galaxy. Accretion of pristine gas from the Intergalactic medium (IGM) is thought to be one of the primary sources for star-forming gas. This gas first passes through the virial regions of the galaxy before reaching the Interstellar medium (ISM), the hub of star formation. On the other hand, owing to the evolutionary course of young and massive stars, energetic winds are ejected from the ISM to the virial regions of the galaxy. A bunch of interlinked, complex astrophysical processes, arising from the concurrent presence of both infalling as well as outbound gas, play out over a range of timescales in the halo region or the Circumgalactic medium (CGM) of a galaxy. It would not be incorrect to say that the CGM has a stronghold over the gas reserves of a galaxy and thus, plays a backhand, yet, rather pivotal role in shaping many galactic properties, some of which are also readily observable. Observing the multi-phase CGM (via spectral-line ion measurements), however, remains a non-trivial effort even today. Low particle densities as well as the CGM’s vast spatial extent, coupled with likely deviations from a spherical distribution, marr the possibility of obtaining complete, unbiased, high-quality spectral information tracing the full extent of the gaseous halo. This often incomplete information leads to multiple inferences about the CGM properties that give rise to multiple contradicting models. In this regard, computer simulations offer a neat solution towards testing and, subsequently, falsifying many of these existing CGM models. Thanks to their controlled environments, simulations are able to not only effortlessly transcend several orders of magnitude in time and space, but also get around many of the observational limitations and provide some unique views on many CGM properties. In this thesis, I focus on effectively using different computer simulations to understand the role of CGM in various astrophysical contexts, namely, the effect of Local Group (LG) environment, major merger events and satellite galaxies. In Chapter 2, I discuss the approach used for modeling various phases of the simulated z = 0 LG CGM in Hestia constrained simulations. Each of the three realizations contain a Milky Way (MW)–Andromeda (M31) galaxy pair, along with their corresponding sets of satellite galaxies, all embedded within the larger cosmological context. For characterizing the different temperature–density phases within the CGM, I model five tracer ions with cloudy ionization modeling. The cold and cool–ionized CGM (H i and Si iii respectively) in Hestia is very clumpy and distributed close to the galactic centers, while the warm-hot and hot CGM (O vi, O vii and O viii) is tenuous and volume-filling. On comparing the H i and Si iii column densities for the simulated M31 with observational measurements from Project AMIGA survey and other low-z galaxies, I found that Hestia galaxies produced less gas in the outer CGM, unlike observations. My carefully designed observational bias model subsequently revealed the possibility that some MW gas clouds might be incorrectly associated with the M31 CGM in observations, and hence, may be partly responsible for giving rise to the detected mismatch between simulated data and observations. In Chapter 3, I present results from four zoom–in, major merger, gas–rich simulations and the subsequent role of the gas, originally situated in the CGM, in influencing some of the galactic observables. The progenitor parameters are selected such that the post–merger remnants are MW–mass galaxies. We generally see a very clear gas bridge joining the merging galaxies in case of multiple passage mergers while such a bridge is mostly absent when a direct collision occurs. On the basis of particle–to–galaxy distance computations and tracer particle analysis, I found that about 33–48 percent of the cold gas contributing to the merger–induced star formation in the bridge originated from the CGM regions. In Chapter 4, I used a sample of 234 MW-mass, L* galaxies from the TNG50 cosmological simulations, with an aim of characterizing the impact of their global satellite populations on the extended cold CGM properties of their host L* halos. On the basis of halo mass and number of satellite galaxies (N_sats ), I categorized the sample into low and high mass bins, and subsequently into bottom, inter and top quartiles respectively. After confirming that satellites indeed influence the extended cold halo gas density profiles of the host galaxies, I investigated the effects of different satellite population parameters on the host halo cold CGMs. My analysis showed that there is hardly any cold gas associated with the satellite population of the lowest mass halos. The stellar mass of the most massive satellite (M_*mms ) impacted the cold gas in low mass bin halos the most, while N_sats (followed by M_*mms ) was the most influential factor for the high mass halos. In any case, how easily cold gas was stripped off the most massive satellite did not play much role. The number of massive (Stellar mass, M* > 10^8 M_solar) satellites as well as the M_*mms associated with a galaxy are two of the most crucial parameters determining how much cold gas ultimately finds its way from the satellites to the host halo. Low mass galaxies are found rather lacking on both these fronts unlike their high mass counterparts. This work highlights some aspects of the complex gas physics that constitute the basic essence of a low-z CGM. My analysis proved the importance of a cosmological environment, local surroundings and merger history in defining some key observable properties of a galactic CGM. Furthermore, I found that different satellite properties were responsible for affecting the cold–dense CGM of the low and high-mass parent galaxies. Finally, the LG emerged as an exciting prospect for testing and pinning down several intricate details about the CGM.
Classical Wolf-Rayet (cWR) stars are at a crucial evolutionary stage for constraining the fates of massive stars. The feedback of these hot, hydrogen-depleted stars dominates their surrounding by tremendous injections of ionizing radiation and kinetic energy. The strength of a Wolf-Rayet (WR) wind decides the eventual mass of its remnant, likely a massive black hole. However, despite their major influence and importance for gravitational wave detection statistics, WR winds are particularly poorly understood. In this paper, we introduce the first set of hydrodynamically consistent stellar atmosphere models for cWR stars of both the carbon (C) and the nitrogen (N) sequence, i.e. WC and WN stars, as a function of stellar luminosity-to-mass ratio (or Eddington Gamma) and metallicity. We demonstrate the inapplicability of the CAK wind theory for cWR stars and confirm earlier findings that their winds are launched at the (hot) iron (Fe) opacity peak. For log Z/Z(circle dot) > -2, Fe is also the main accelerator throughout the wind. Contrasting previous claims of a sharp lower mass-loss limit forWR stars, we obtain a smooth transition to optically thin winds. Furthermore, we find a strong dependence of the mass-loss rates on Eddington Gamma, both at solar and subsolar metallicity. Increases inWCcarbon and oxygen abundances turn out to slightly reduce the predicted mass-loss rates. Calculations at subsolar metallicities indicate that below the metallicity of the Small Magellanic Cloud, WR mass-loss rates decrease much faster than previously assumed, potentially allowing for high black hole masses even in the local Universe.
X-ray observations of star-planet systems are important to grow our understanding of exoplanets; these observations allow for studies of photoevaporation of the exoplanetary atmosphere, and in some cases even estimations of the size of the outer planetary atmosphere. The German-Russian eROSITA instrument onboard the SRG (Spectrum Roentgen Gamma) mission is performing the first all-sky X-ray survey since the 1990s, and provides X-ray fluxes and spectra of exoplanet host stars over a much larger volume than was accessible before. Using new eROSITA data as well as archival data from XMM-Newton, Chandra, and ROSAT, we estimate mass-loss rates of exoplanets under an energy-limited escape scenario and identify several exoplanets with strong X-ray irradiation and expected mass loss that are amenable to follow-up observations at other wavelengths. We model sample spectra using a toy model of an exoplanetary atmosphere to predict what exoplanet transit observations with future X-ray missions such as Athena will look like and estimate the observable X-ray transmission spectrum for a typical hot Jupiter-type exoplanet.
Alone but not lonely
(2020)
Context. Hot subdwarfs are core-helium burning stars that show lower masses and higher temperatures than canonical horizontal branch stars. They are believed to be formed when a red giant suffers an extreme mass-loss episode. Binary interaction is suggested to be the main formation channel, but the high fraction of apparently single hot subdwarfs (up to 30%) has prompted single star formation scenarios to be proposed.Aims. We investigate the possibility that hot subdwarfs could form without interaction by studying wide binary systems. If single formation scenarios were possible, there should be hot subdwarfs in wide binaries that have undergone no interaction.Methods. Angular momentum accretion during interaction is predicted to cause the hot subdwarf companion to spin up to the critical velocity. The effect of this should still be observable given the timescales of the hot subdwarf phase. To study the rotation rates of companions, we have analysed light curves from the Transiting Exoplanet Survey Satellite for all known hot subdwarfs showing composite spectral energy distributions indicating the presence of a main sequence wide binary companion. If formation without interaction were possible, that would also imply the existence of hot subdwarfs in very wide binaries that are not predicted to interact. To identify such systems, we have searched for common proper motion companions with projected orbital distances of up to 0.1 pc to all known spectroscopically confirmed hot subdwarfs using Gaia DR2 astrometry.Results. We find that the companions in composite hot subdwarfs show short rotation periods when compared to field main sequence stars. They display a triangular-shaped distribution with a peak around 2.5 days, similar to what is observed for young open clusters. We also report a shortage of hot subdwarfs with candidate common proper motion companions. We identify only 16 candidates after probing 2938 hot subdwarfs with good astrometry. Out of those, at least six seem to be hierarchical triple systems, in which the hot subdwarf is part of an inner binary.Conclusions. The observed distribution of rotation rates for the companions in known wide hot subdwarf binaries provides evidence of previous interaction causing spin-up. Additionally, there is a shortage of hot subdwarfs in common proper motion pairs, considering the frequency of such systems among progenitors. These results suggest that binary interaction is always required for the formation of hot subdwarfs.
High-resolution spectroscopy of an erupting minifilament and its impact on the nearby chromosphere
(2020)
We study the evolution of a minifilament eruption in a quiet region at the center of the solar disk and its impact on the ambient atmosphere. We used high spectral resolution imaging spectroscopy in H alpha acquired by the echelle spectrograph of the Vacuum Tower Telescope, Tenerife, Spain; photospheric magnetic field observations from the Helioseismic Magnetic Imager; and UV/EUV imaging from the Atmospheric Imaging Assembly of the Solar Dynamics Observatory. The H alpha line profiles were noise-stripped using principal component analysis and then inverted to produce physical and cloud model parameter maps. The minifilament formed between small-scale, opposite-polarity magnetic features through a series of small reconnection events, and it erupted within an hour after its appearance in H alpha. Its development and eruption exhibited similarities to large-scale erupting filaments, indicating the action of common mechanisms. Its eruption took place in two phases, namely, a slow rise and a fast expansion, and it produced a coronal dimming, before the minifilament disappeared. During its eruption, we detected a complicated velocity pattern, indicative of a twisted, thread-like structure. Part of its material returned to the chromosphere, producing observable effects on nearby low-lying magnetic structures. Cloud model analysis showed that the minifilament was initially similar to other chromospheric fine structures, in terms of optical depth, source function, and Doppler width, but it resembled a large-scale filament on its course to eruption. High spectral resolution observations of the chromosphere can provide a wealth of information regarding the dynamics and properties of minifilaments and their interactions with the surrounding atmosphere.
A conducting Taylor-Couette flow with quasi-Keplerian rotation law containing a toroidal magnetic field serves as a mean-field dynamo model of the Tayler-Spruit type. The flows are unstable against non-axisymmetric perturbations which form electromotive forces defining a effect and eddy diffusivity. If both degenerated modes with m = +/- 1 are excited with the same power then the global a effect vanishes and a dynamo cannot work. It is shown, however, that the Tayler instability produces finite alpha effects if only an isolated mode is considered but this intrinsic helicity of the single-mode is too low for an alpha(2) dynamo. Moreover, an alpha Omega dynamo model with quasi-Keplerian rotation requires a minimum magnetic Reynolds number of rotation of Rm similar or equal to 2000 to work. Whether it really works depends on assumptions about the turbulence energy. For a steeper-than-quadratic dependence of the turbulence intensity on the magnetic field, however, dynamos are only excited if the resulting magnetic eddy diffusivity approximates its microscopic value, eta(T) similar or equal to eta. By basically lower or larger eddy diffusivities the dynamo instability is suppressed.
Stars are uniform spheres, but only to first order. The way in which stellar rotation and magnetism break this symmetry places important observational constraints on stellar magnetic fields, and factors in the assessment of the impact of stellar activity on exoplanet atmospheres. The spatial distribution of flares on the solar surface is well known to be nonuniform, but elusive on other stars. We briefly review the techniques available to recover the loci of stellar flares, and highlight a new method that enables systematic flare localization directly from optical light curves. We provide an estimate of the number of flares we may be able to localize with the Transiting Exoplanet Survey Satellite, and show that it is consistent with the results obtained from the first full sky scan of the mission. We suggest that nonuniform flare latitude distributions need to be taken into account in accurate assessments of exoplanet habitability.
We explore the possibility to systematically study the extended, hot gaseous halos of low-redshift galaxies with coronal broad Lya absorbers (CBLAs). These are weak, thermally broadenend H I absorption lines arising from the tiny fraction of neutral hydrogen that resides in the collisionally ionized, million-degree halo gas in these galaxies. Using a semi-analytic approach, we model the spatial density and temperature distribution of hot coronal gas to predict strength, spectral shape, and cross section of CBLAs as a function of galaxy-halo mass and line-of-sight impact parameter. For virial halo masses in the range log M M = 10.6 12.6, the characteristic logarithmic CBLA H I column densities and Doppler parameters are log N(H I) = 12.4- 13.4 and b(H I).=.70-200 km s-1, indicating that CBLAs represent weak, shallow spectral features that are difficult to detect. Yet, the expected number density of CBLAs per unit redshift in the above given mass range is d. dz(CBLA). 3, implying that CBLAs have a substantial absorption cross section. We compare the model predictions with a combined set of UV absorption-line spectra from the Hubble Space Telescope (HST)/Cosmic Origins Spectrograph and HST/Space Telescope Imaging Spectrograph that trace the halos of four low-redshift galaxies. We demonstrate that CBLAs might already have been detected in these spectra, but the complex multi-component structure and the limited signal-to-noise ratio complicate the interpretation of these CBLA candidate systems. Our study suggests that CBLAs represent a very interesting absorber class that potentially will allow us to further explore the hot coronae of galaxies with UV spectral data.
Context
The recently claimed discovery of a massive (M-BH = 68(-13)(+11) M-circle dot) black hole in the Galactic solar neighborhood has led to controversial discussions because it severely challenges our current view of stellar evolution.
Aims
A crucial aspect for the determination of the mass of the unseen black hole is the precise nature of its visible companion, the B-type star LSV +22 25. Because stars of different mass can exhibit B-type spectra during the course of their evolution, it is essential to obtain a comprehensive picture of the star to unravel its nature and, thus, its mass.
Methods
To this end, we study the spectral energy distribution of LSV +22 25 and perform a quantitative spectroscopic analysis that includes the determination of chemical abundances for He, C, N, O, Ne, Mg, Al, Si, S, Ar, and Fe.
Results
Our analysis clearly shows that LSV +22 25 is not an ordinary main sequence B-type star. The derived abundance pattern exhibits heavy imprints of the CNO bi-cycle of hydrogen burning, that is, He and N are strongly enriched at the expense of C and O. Moreover, the elements Mg, Al, Si, S, Ar, and Fe are systematically underabundant when compared to normal main-sequence B-type stars. We suggest that LSV +22 25 is a stripped helium star and discuss two possible formation scenarios. Combining our photometric and spectroscopic results with the Gaia parallax, we infer a stellar mass of 1.1 +/- 0.5 M-circle dot. Based on the binary system's mass function, this yields a minimum mass of 2-3 M-circle dot for the compact companion, which implies that it may not necessarily be a black hole but a massive neutron- or main sequence star.
Conclusions
The star LSV +22 25 has become famous for possibly having a very massive black hole companion. However, a closer look reveals that the star itself is a very intriguing object. Further investigations are necessary for complete characterization of this object.
Multiphase galaxy winds, the accretion of cold gas through galaxy haloes, and gas stripping from jellyfish galaxies are examples of interactions between cold and hot gaseous phases. There are two important regimes in such systems. A sufficiently small cold cloud is destroyed by the hot wind as a result of Kelvin-Helmholtz instabilities, which shatter the cloud into small pieces that eventually mix and dissolve in the hot wind. In contrast, stripped cold gas from a large cloud mixes with the hot wind to intermediate temperatures, and then becomes thermally unstable and cools, causing a net accretion of hot gas to the cold tail. Using the magneto-hydrodynamical code AREPO, we perform cloud crushing simulations and test analytical criteria for the transition between the growth and destruction regimes to clarify a current debate in the literature. We find that the hot-wind cooling time sets the transition radius and not the cooling time of the mixed phase. Magnetic fields modify the wind-cloud interaction. Draping of wind magnetic field enhances the field upstream of the cloud, and fluid instabilities are suppressed by a turbulently magnetized wind beyond what is seen for a wind with a uniform magnetic field. We furthermore predict jellyfish galaxies to have ordered magnetic fields aligned with their tails. We finally discuss how the results of idealized simulations can be used to provide input to subgrid models in cosmological (magneto-)hydrodynamical simulations, which cannot resolve the detailed small-scale structure of cold gas clouds in the circumgalactic medium.
The central gas in half of all galaxy clusters shows short cooling times. Assuming unimpeded cooling, this should lead to high star formation and mass cooling rates, which are not observed. Instead, it is believed that condensing gas is accreted by the central black hole that powers an active galactic nuclei jet, which heats the cluster. The detailed heating mechanism remains uncertain. A promising mechanism invokes cosmic ray protons that scatter on self-generated magnetic fluctuations, i.e. Alfvén waves. Continuous damping of Alfvén waves provides heat to the intracluster medium. Previous work has found steady state solutions for a large sample of clusters where cooling is balanced by Alfvénic wave heating. To verify modeling assumptions, we set out to study cosmic ray injection in three-dimensional magnetohydrodynamical simulations of jet feedback in an idealized cluster with the moving-mesh code arepo. We analyze the interaction of jet-inflated bubbles with the turbulent magnetized intracluster medium.
Furthermore, jet dynamics and heating are closely linked to the largely unconstrained jet composition. Interactions of electrons with photons of the cosmic microwave background result in observational signatures that depend on the bubble content. Those recent observations provided evidence for underdense bubbles with a relativistic filling while adopting simplifying modeling assumptions for the bubbles. By reproducing the observations with our simulations, we confirm the validity of their modeling assumptions and as such, confirm the important finding of low-(momentum) density jets.
In addition, the velocity and magnetic field structure of the intracluster medium have profound consequences for bubble evolution and heating processes. As velocity and magnetic fields are physically coupled, we demonstrate that numerical simulations can help link and thereby constrain their respective observables. Finally, we implement the currently preferred accretion model, cold accretion, into the moving-mesh code arepo and study feedback by light jets in a radiatively cooling magnetized cluster. While self-regulation is attained independently of accretion model, jet density and feedback efficiencies, we find that in order to reproduce observed cold gas morphology light jets are preferred.
The Magellanic Bridge, stretching between the Small and the Large Magellanic Cloud (SMC and LMC), is the nearest tidally stripped intergalactic environment. The Bridge has a significantly low average metallicity of Z less than or similar to 0.1 Z(circle dot). Here we report the first discovery of O-type stars in the Magellanic Bridge. Three massive O stars were identified thanks to the archival spectra obtained by the ESO's Very Large Telescope FLAMES instrument. We analyze the spectra of each star using the Potsdam Wolf-Rayet (PoWR) non-local thermodynamic equilibrium model atmosphere code, which provides the physical parameters, ionizing photon fluxes, and surface abundances. The ages of the newly discovered O stars suggest that star formation in the Bridge is ongoing. Furthermore, the discovery of O stars in the Bridge implies that tidally stripped galactic tails containing low-density but highly dynamical gas are capable of producing massive O stars. The multi-epoch spectra indicate that all three O stars are binaries. Despite their spatial proximity to one another, these O stars are chemically distinct. One of them is a fast-rotating giant with nearly LMC-like abundances. The other two are main-sequence stars that rotate extremely slowly and are strongly metal depleted. We discover the most nitrogen-poor O star known to date. Taking into account the previous analyses of B stars in the Bridge, we interpret the various metal abundances as the signature of a chemically inhomogeneous interstellar medium (ISM), suggesting that the Bridge gas might have accreted during multiple episodes of tidal interaction between the Clouds. Attributing the lowest derived metal content to the primordial gas, the time of the initial formation of the Bridge may date back several billion years. Using the Gaia and Galex color-magnitude diagrams, we roughly estimate the total number of O stars in the Bridge and their total ionizing radiation. Comparing this with the energetics of the diffuse ISM, we find that the contribution of the hot stars to the ionizing radiation field in the Bridge is less than 10% and conclude that the main sources of ionizing photons are leaks from the LMC and SMC. We estimate a lower limit for the fraction of ionizing radiation that escapes from these two dwarf galaxies.
We present an efficient technique for control of synchrony in a globally coupled ensemble by pulsatile action. We assume that we can observe the collective oscillation and can stimulate all elements of the ensemble simultaneously. We pay special attention to the minimization of intervention into the system. The key idea is to stimulate only at the most sensitive phase. To find this phase, we implement an adaptive feedback control. Estimating the instantaneous phase of the collective mode on the fly, we achieve efficient suppression using a few pulses per oscillatory cycle. We discuss the possible relevance of the results for neuroscience, namely, for the development of advanced algorithms for deep brain stimulation, a medical technique used to treat Parkinson's disease.
In light of substantial new discoveries of hot subdwarfs by ongoing spectroscopic surveys and the availability of new all-sky data from ground-based photometric surveys and the Gaia mission Data Release 2, we compiled an updated catalogue of the known hot subdwarf stars. The catalogue contains 5874 unique sources including 528 previously unknown hot subdwarfs and provides multi-band photometry, astrometry from Gaia, and classifications based on spectroscopy and colours. This new catalogue provides atmospheric parameters of 2187 stars and radial velocities of 2790 stars from the literature. Using colour, absolute magnitude, and reduced proper motion criteria, we identified 268 previously misclassified objects, most of which are less luminous white dwarfs or more luminous blue horizontal branch and main-sequence stars.
In the last century, several astronomical measurements have supported that a significant percentage (about 22%) of the total mass of the Universe, on galactic and extragalactic scales, is composed of a mysterious ”dark” matter (DM). DM does not interact with the electromagnetic force; in other words it does not reflect, absorb or emit light. It is possible that DM particles are weakly interacting massive particles (WIMPs) that can annihilate (or decay) into Standard Model (SM) particles, and modern very- high-energy (VHE; > 100 GeV) instruments such as imaging atmospheric Cherenkov telescopes (IACTs) can play an important role in constraining the main properties of such DM particles, by detecting these products. One of the most privileged targets where to look for DM signal are dwarf spheroidal galaxies (dSphs), as they are expected to be high DM-dominated objects with a clean, gas-free environment. Some dSphs could be considered as extended sources, considering the angular resolution of IACTs; their angu- lar resolution is adequate to detect extended emission from dSphs. For this reason, we performed an extended-source analysis, by taking into account in the unbinned maximum likelihood estimation both the energy and the angular extension dependency of observed events. The goal was to set more constrained upper limits on the velocity-averaged cross-section annihilation of WIMPs with VERITAS data. VERITAS is an array of four IACTs, able to detect γ-ray photons ranging between 100 GeV and 30 TeV. The results of this extended analysis were compared against the traditional spectral analysis. We found that a 2D analysis may lead to more constrained results, depending on the DM mass, channel, and source. Moreover, in this thesis, the results of a multi-instrument project are presented too. Its goal was to combine already published 20 dSphs data from five different experiments, such as Fermi-LAT, MAGIC, H.E.S.S., VERITAS and HAWC, in order to set upper limits on the WIMP annihilation cross-section in the widest mass range ever reported.
Cosmic rays (CRs) constitute an important component of the interstellar medium (ISM) of galaxies and are thought to play an essential role in governing their evolution. In particular, they are able to impact the dynamics of a galaxy by driving galactic outflows or heating the ISM and thereby affecting the efficiency of star-formation. Hence, in order to understand galaxy formation and evolution, we need to accurately model this non-thermal constituent of the ISM. But except in our local environment within the Milky Way, we do not have the ability to measure CRs directly in other galaxies. However, there are many ways to indirectly observe CRs via the radiation they emit due to their interaction with magnetic and interstellar radiation fields as well as with the ISM.
In this work, I develop a numerical framework to calculate the spectral distribution of CRs in simulations of isolated galaxies where a steady-state between injection and cooling is assumed. Furthermore, I calculate the non-thermal emission processes arising from the modelled CR proton and electron spectra ranging from radio wavelengths up to the very high-energy gamma-ray regime.
I apply this code to a number of high-resolution magneto-hydrodynamical (MHD) simulations of isolated galaxies, where CRs are included. This allows me to study their CR spectra and compare them to observations of the CR proton and electron spectra by the Voyager-1 satellite and the AMS-02 instrument in order to reveal the origin of the measured spectral features.
Furthermore, I provide detailed emission maps, luminosities and spectra of the non-thermal emission from our simulated galaxies that range from dwarfs to Milk-Way analogues to starburst galaxies at different evolutionary stages. I successfully reproduce the observed relations between the radio and gamma-ray luminosities with the far-infrared (FIR) emission of star-forming (SF) galaxies, respectively, where the latter is a good tracer of the star-formation rate. I find that highly SF galaxies are close to the limit where their CR population would lose all of their energy due to the emission of radiation, whereas CRs tend to escape low SF galaxies more quickly. On top of that, I investigate the properties of CR transport that are needed in order to match the observed gamma-ray spectra.
Furthermore, I uncover the underlying processes that enable the FIR-radio correlation (FRC) to be maintained even in starburst galaxies and find that thermal free-free-emission naturally explains the observed radio spectra in SF galaxies like M82 and NGC 253 thus solving the riddle of flat radio spectra that have been proposed to contradict the observed tight FRC.
Lastly, I scrutinise the steady-state modelling of the CR proton component by investigating for the first time the influence of spectrally resolved CR transport in MHD simulations on the hadronic gamma-ray emission of SF galaxies revealing new insights into the observational signatures of CR transport both spectrally and spatially.