Non-local Thermodynamic Equilibrium Research Articles

An improved and extended parameterization of the CO2 15 µm cooling in the middle and upper atmosphere (CO2_cool_fort-1.0)

Abstract. The radiative infrared cooling of CO2 in the middle atmosphere, where it emits under non-local thermodynamic equilibrium (non-LTE) conditions, is a crucial contribution to the energy balance of this region and hence to establishing its thermal structure. The non-LTE computation is too CPU time-consuming to be fully incorporated into climate models, and hence it is parameterized. The most used parameterization of the CO2 15 µm cooling for Earth's middle and upper atmosphere was developed by Fomichev et al. (1998). The valid range of this parameterization with respect to CO2 volume mixing ratios (VMRs) is, however, exceeded by the CO2 of several scenarios considered in the Coupled Climate Model Intercomparison Projects, in particular the abrupt-4×CO2 experiment. Therefore, an extension, as well as an update, of that parameterization is both needed and timely. In this work, we present an update of that parameterization that now covers CO2 volume mixing ratios in the lower atmosphere from ∼0.5 to over 10 times the CO2 pre-industrial value of 284 ppmv (i.e. 150 to 3000 ppmv). Furthermore, it is improved by using a more contemporary CO2 line list and the collisional rates that affect the CO2 cooling rates. Overall, its accuracy is improved when tested for the reference temperature profiles as well as for measured temperature fields covering all expected conditions (latitude and season) of the middle atmosphere. The errors obtained for the reference temperature profiles are below 0.5 K d−1 for the present-day and lower CO2 VMRs. Those errors increase to ∼1–2K d−1 at altitudes between 110 and 120 km for CO2 concentrations of 2 to 3 times the pre-industrial values. For very high CO2 concentrations (4 to 10 times the pre-industrial abundances), those errors are below ∼1 K d−1 for most regions and conditions, except at 107–135 km, where the parameterization overestimates them by ∼1.2 %. These errors are comparable to the deviation of the non-LTE cooling rates with respect to LTE at about 70 km and below, but they are negligible (several times smaller) above that altitude. When applied to a large dataset of global (pole to pole and four seasons) temperature profiles measured by MIPAS (Michelson Interferometer for Passive Atmospheric Spectroscopy) (middle- and upper-atmosphere mode), the errors of the parameterization for the mean cooling rate (bias) are generally below 0.5 K d−1, except between 5×10-3 and 3×10-4 hPa (∼85–98 km), where they can reach biases of 1–2 K d−1. For single-temperature profiles, the cooling rate error (estimated by the root mean square – rms – of a statistically significant sample) is about 1–2 K d−1 below 5×10-3 hPa (∼85 km) and above 2×10-4 hPa (∼102 km). In the intermediate region, however, it is between 2 and 7 K d−1. For elevated stratopause events, the parameterization underestimates the mean cooling rates by 3–7 K d−1 (∼10 %) at altitudes of 85–95 km and the individual cooling rates show a significant rms (5–15 K d−1). Further, we have also tested the parameterization for the temperature obtained by a high-resolution version of the Whole Atmosphere Community Climate Model (WACCM-X), which shows a large temperature variability and wave structure in the middle atmosphere. In this case, the mean (bias) error of the parameterization is very small, smaller than 0.5 K d−1 for most atmospheric layers, reaching only maximum values of 2 K d−1 near 5×10-4 hPa (∼ 96 km). The rms has values of 1–2 K d−1 (∼20 %) below ∼2×10-2 hPa (∼80 km) and values smaller than 4 K d−1 (∼2 %) above 10−4 hPa (∼105 km). In the intermediate region between ∼5×10-3 and ∼2×10-4 hPa (85–102 km), the rms is in the range of 5–12 K d−1. While these values are significant in percentage at ∼5×10-3–5×10-4 hPa, they are very small above ∼5×10-4 hPa (96 km). The routine is very fast, taking (1.5–7.5) ×10-5 s, depending on the extension of the atmospheric profile, the processor and the Fortran compiler.

Peculiarities of the chemical enrichment of metal-poor stars in the Milky Way Galaxy

The oldest stars in the Milky Way are metal-poor with Fe/H $<$ -- 1.0, displaying peculiar elemental abundances compared to solar values. The relative variations in the chemical compositions among stars is also increasing with decreasing stellar metallicity, allowing for the pure signature of unique nucleosynthesis processes to be revealed. The study of the $r$-process is, for instance, one of the main goals of stellar archaeology and metal-poor stars exhibit an unexpected complexity in the stellar production of the $r$-process elements in the early Galaxy. In this work, we report the atmospheric parameters, main dynamic properties, and the abundances of four metal-poor stars: HE 1523-–0901, HD 6268, HD 121135, and HD 195636 (--1.5 $>$ Fe/H The abundances were derived from spectra obtained with the HRS echelle spectrograph at the Southern African Large Telescope, using both local and non-local thermodynamic equilibrium (LTE and NLTE) approaches, with the average error between 0.10 and 0.20 dex. Based on their kinematical properties, we show that HE 1523--0901 and HD 195636 are halo stars with typical high velocities. In particular, HD 121135 displays a peculiar kinematical behaviour, making it unclear whether it is a halo or an accreted star. Furthermore, HD 6268 is possibly a rare prototype of very metal-poor thick disk stars. The abundances derived for our stars are compared with theoretical stellar models and with other stars with similar metallicity values from the literature. HD 121135 is Al-poor and Sc-poor, compared to stars observed in the same metallicity range (--1.62 $>$ Fe/H $>$--1.12). The most metal-poor stars in our sample, HE 1523 -- 0901, HD 6268, and HD 195636, exhibit anomalies that are better explained by supernova models from fast-rotating stellar progenitors for elements up to the Fe group. Compared to other stars in the same metallicity range, their common biggest anomaly is represented by the low Sc abundances. If we consider the elements beyond Zn, HE 1523--0901 can be classified as an r-II star, HD 6268 as an r-I candidate, and HD 195636 and HD 121135 exhibiting a borderline $r$-process enrichment between limited-r and r-I star. Significant relative differences are observed between the $r$-process signatures in these stars.

M3DIS - A grid of 3D radiation-hydrodynamics stellar atmosphere models for stellar surveys. I. Procedure, validation, and the Sun

Large-scale stellar surveys, such as SDSS-V, 4MOST, WEAVE, and PLATO, require accurate atmospheric models and synthetic spectra of stars for accurate analyses of fundamental stellar parameters and chemical abundances. The primary goal of our work is to develop a new approach to solve radiation-hydrodynamics (RHD) and generate model stellar spectra in a self-consistent and highly efficient framework. We build upon the Copenhagen legacy RHD code, the MULTI3D non-local thermodynamic equilibrium (NLTE) code, and the DISPATCH high-performance framework. The new approach allows us to calculate 3D RHD models of stellar atmospheres on timescales of a few thousand CPU hours and to perform subsequent spectrum synthesis in local thermodynamic equilibrium (LTE) or NLTE for the desired physical conditions within the parameter space of FGK-type stars. We compare the 3D RHD solar model with other available models and validate its performance against solar observations, including the centre-to-limb variation of intensities and key solar diagnostic lines of H and Fe. We show that the performance of the new code allows to overcome the main bottleneck in 3D NLTE spectroscopy and enables calculations of multi-dimensional grids of synthetic stellar observables for comparison with modern astronomical observations.

EMISSA: Exploring millimetre indicators of solar-stellar activity III. Comparison of Ca II indices and millimetre continua in a 3D model atmosphere

Amongst several spectral lines, some of the strongest chromospheric diagnostics are offered by the Ca II H K lines. These lines can be used to gauge the temperature stratification of the atmosphere since the line core and wings are formed in different regions of the solar atmosphere. Furthermore, the Ca II lines act as tracers for the magnetic structure of the solar atmosphere, as the line cores are formed in the upper chromosphere even though they are formed in non-local thermodynamic equilibrium (NLTE). In contrast, the formation of millimetre (mm) continuum radiation occurs under local thermodynamic equilibrium (LTE) conditions. As a result, the brightness temperatures obtained from observations with the Atacama Large Millimetre/Submillimetre Array (ALMA) offer a complementary perspective on the activity and thermal structure of stellar atmospheres. The overall aim is to establish more robust solar/stellar activity indicators using ALMA observations in comparison with classical diagnostics, such as the s index and infrared triplet (IRT) index. We employed the 1.5D radiative transfer codes RH1.5D and advanced radiative transfer (ART) to compute the synthetic spectra for the Ca II lines and the millimetre (mm) continua, respectively. These calculations were performed using an enhanced network atmosphere model, which incorporates non-equilibrium hydrogen ionisation generated by the state-of-the-art 3D radiation magnetohydrodynamics (rMHD) Bifrost code. To account for the limited spatial resolution of ALMA, we simulated the effect using a Gaussian point spread function (PSF). Additionally, we analysed the correlations and slopes of scatter plots between the Ca II indices and mm continuum for the original and degraded resolutions, focusing on the entire simulation box, quiet Sun regions, and enhanced network patches separately. The activity indices generated from these lines could further be used to compare the spectra of Sun-like stars with the solar spectrum. We present a comparative study between synthetic continuum brightness temperature maps at mm wavelengths (0.3 mm to 8.5 mm) and the Ca II activity indices; namely, the s index and infrared triplet (IRT) index. The Ca II activity indices and mm brightness temperatures are weakly correlated at the high resolution, with the highest correlation observed at a wavelength of 0.3 mm, corresponding to ALMA band 10. As the resolution decreases, the correlation consistently increases. Conversely, the slopes exhibit a decreasing trend with increasing wavelength, while the degradation of resolution does not noticeably affect the calculated slopes. As the spatial resolution decreases, the standard deviations of the Ca II activity indices and brightness temperatures decrease, while the correlations between them increase. However, the slopes do not exhibit significant changes. Consequently, these relationships could be valuable for calibrating the mm continuum maps obtained through ALMA observations.

OCCASO V. Chemical-abundance trends with Galactocentric distance and age

Open clusters provide valuable information on stellar nucleosynthesis and the chemical evolution of the Galactic disk, as their age and distances can be measured more precisely with photometry than in the case of field stars. Our aim is to study the chemical distribution of the Galactic disk using open clusters by analyzing the existence of gradients with Galactocentric distance, azimuth, or height from the plane and dependency with age. We used the high-resolution spectra ($R$>60\,000) of 194 stars belonging to 36 open clusters to determine the atmospheric parameters and chemical abundances with two independent methods: equivalent widths and spectral synthesis. The sample was complemented with 63 clusters with high-resolution spectroscopy from literature. We measured LTE abundances for 21 elements: alpha (Mg, Si, Ca, and Ti), odd-Z (Na and Al), Fe-peak (Fe, Sc, V, Cr, Mn, Co, Ni, Cu, and Zn), and neutron-capture (Sr, Y, Zr, Ba, Ce, and Nd). We also provide non-local thermodynamic equilibrium abundances for elements when corrections are available. We find inner disk young clusters enhanced in Mg/Fe and Si/Fe compared to other clusters of their age. For Ba/Fe we report an age trend flattening for older clusters (age$<$2.5\,Ga). The studied elements follow the expected radial gradients as a function of their nucleosynthesis groups, which are significantly steeper for the oldest systems. For the first time, we investigate the existence of an azimuthal gradient, finding some hints of its existence among the old clusters (age$>$2\,Ga).

Testing analytical methods to derive the cosmic-ray ionisation rate in cold regions via synthetic observations

Context. Cosmic rays (CRs) heavily impact the chemistry and physics of cold and dense star-forming regions. However, the characterisation of their ionisation rate continues to pose a challenge from the observational point of view. Aims. In the past, a few analytical formulas have been proposed to infer the cosmic-ray ionisation rate, ζ2, from molecular line observations. These have been derived from the chemical kinetics of the involved species, but they have not yet been validated using synthetic data processed with a standard observative pipeline. In this work, we aim to bridge this gap. Methods. We performed a radiative transfer on a set of three-dimensional magneto-hydrodynamical simulations of prestellar cores, exploring different initial ζ2, evolutionary stages, types of radiative transfer (for instance assuming local-thermodynamic-equilibrium conditions), and telescope responses. We then computed the column densities of the involved tracers to determine ζ2, employing a recently proposed method based on the detection of H2D+. We compared this approach with a previous method, based on more common tracers. Both approaches are commonly used. Results. Our results confirm that the equation based on the detection of H2D+ accurately retrieves the actual ζ2 within a factor of two to three in the physical conditions explored in our tests. Since we have also explored a non-local thermodynamic equilibrium (non-LTE) radiative transfer, this work indirectly offers insights into the excitation temperatures of common transitions at moderate volume densities (n ≈ 105 cm−3). We also performed a few tests using a previous methodology that is independent of H2D+, which overestimates the actual ζ2 by at least two orders of magnitude. We considered a new derivation of this method, however, we found that it still leads to high over-estimations. Conclusions. The method based on H2D+ is further validated in this work and demonstrates a reliable method for estimating ζ2 in cold and dense gas. On the contrary, the former analytical equation, as already pointed out by its authors, has no global domain of application. Thus, we find that it ought to be employed with caution.

One-dimensional, geometrically stratified semi-empirical models of the quiet-Sun photosphere and lower chromosphere

One-dimensional, semi-empirical models of the solar atmosphere are widely employed in numerous contexts within solar physics, ranging from the determination of element abundances and atomic parameters to studies of the solar irradiance and from Stokes inversions to coronal extrapolations. These models provide the physical parameters (i.e. temperature, gas pressure, etc.) in the solar atmosphere as a function of the continuum optical depth $ c $. The transformation to the geometrical $z$ scale (i.e. vertical coordinate) is provided via vertical hydrostatic equilibrium. Our aim is to provide updated, one-dimensional, semi-empirical models of the solar atmosphere as a function of $z,$ but employing the more general case of three-dimensional magneto-hydrostatic equilibrium (MHS) instead of vertical hydrostatic equilibrium (HE). We employed a recently developed Stokes inversion code that, along with non-local thermodynamic equilibrium effects, considers MHS instead of HE. This code is applied to spatially and temporally resolved spectropolarimetric observations of the quiet Sun obtained with the CRISP instrument attached to the Swedish Solar Telescope. We provide average models for granules, intergranules, dark magnetic elements, and overall quiet-Sun as a function of both $ c $ and $z$ from the photosphere to the lower chromosphere. We demonstrate that, in these quiet-Sun models, the effect of considering MHS instead of HE is negligible. However, employing MHS increases the consistency of the inversion results before averaging. We surmise that in regions with stronger magnetic fields (i.e. pores, sunspots, network) the benefits of employing the magneto-hydrostatic approximation will be much more palpable.

NLTE modelling of water-rich exoplanet atmospheres. Cooling and heating rates

The hydrogen and water molecules respond very differently to the collisional-radiative processes taking place in planetary atmospheres. Naturally, the question arises whether H2O-rich atmospheres are more (or less) resilient to long-term mass loss than H2-dominated ones if they radiate away the incident stellar energy more (or less) efficiently. If confirmed, the finding would have implications on our understanding of the evolution of sub-Neptune exoplanets. As a key step towards answering this question, we present a non-local thermodynamic equilibrium (NLTE) model of H2O for the atmospheric region where the gas accelerates to escape the planet and conditions relevant to close-in sub-Neptunes. Our exploratory calculations for isothermal gas composed of H2, H2O and e− reveal that: (1) In the pressure region ∼10−2–10−4 dyn cm−2 where the stellar extreme-ultraviolet (XUV) photons are typically deposited in the atmosphere, H2O is in rotational LTE but vibrational NLTE. Vibrational LTE is facilitated by high H2O abundances and fractional ionizations, and we report critical densities for the LTE-NLTE transition; (2) Vibrational cooling may locally dominate over rotational cooling, partly because of the comparatively small opacities of ro-vibrational lines; (3) Even low H2O abundances notably enhance the cooling, foreseeably offsetting some of the stellar heating; (4) Heating due to the deposition of stellar infrared (IR) photons is significant at pressures ≳0.1 dyn cm−2. We estimate the contribution of H2O excitation to the internal energy of the gas and speculate on the photodissociation from the excited vibrational states. Ultimately, our findings motivate the consideration of NLTE in the mass loss rate calculations of H2O-rich atmospheres.

A possibly solar metallicity atmosphere escaping from HAT-P-32b revealed by Halpha and He absorption

This paper presents a hydrodynamic simulation that couples detailed non-local thermodynamic equilibrium (NLTE) calculations of the helium and hydrogen level populations to model the Halpha and He 10830 transmission spectra of the hot Jupiter HAT-P-32b. A Monte Carlo simulation was applied to calculate the number of Lyalpha resonance scatterings, which is the main process for populating H(2). In the examined parameter space, only models with H/He geq 99.5/0.5, $(0.5 3.0)$ times the fiducial value of XUV $, and spectral index $ 0.3)$, can explain the Halpha and He 10830 lines simultaneously. We found a mass-loss rate of $ $ g s$^ $, consistent with previous studies. Moreover, we found that the stellar Lyalpha flux should be as high as $4 $ erg cm$^ $ s$^ $, indicating high stellar activity during the observation epoch of the two absorption lines. Despite the fact that the metallicity in the lower atmosphere of HAT-P-32b may be super-solar, our simulations tentatively suggest it is close to solar in the upper atmosphere. Understanding the difference in metallicity between the lower and upper atmospheres is essential for future atmospheric characterisations.

Precise photoionisation treatment and hydrodynamic effects in atmospheric modelling of warm and hot Neptunes

Context. Observational breakthroughs in the field of exoplanets in the last decade have motivated the development of numerous theoretical models, such as those describing atmospheres and mass loss, which is believed to be one of the main drivers of planetary evolution. Aims. We outline for which types of close-in planets in the Neptune-mass range the accurate treatment of photoionisation effects is most relevant, particularly with respect to atmospheric mass loss and the parameters relevant for interpreting observations, such as temperature and ion fraction. Methods. We developed the Cloudy e Hydro Ancora INsieme (CHAIN) model, which combines our 1D hydrodynamic upper atmosphere model with the non-local thermodynamical equilibrium (NLTE) photoionisation and radiative transfer code Cloudy accounting for ionisation, dissociation, detailed atomic level populations, and chemical reactions for all chemical elements up to zinc. The hydro-dynamic code is responsible for describing the outflow, while Cloudy solves the photoionisation and heating of planetary atmospheres. We applied CHAIN to model the upper atmospheres of a range of Neptune-like planets with masses between 1 and 50 M⊕, also varying the orbital parameters. Results. For the majority of warm and hot Neptunes, we find slower and denser outflows, with lower ion fractions, compared to the predictions of the hydrodynamic model alone. Furthermore, we find significantly different temperature profiles between CHAIN and the hydrodynamic model alone, though the peak values are similar for similar atmospheric compositions. The mass-loss rates predicted by CHAIN are higher for hot strongly irradiated planets and lower for more moderate planets. All differences between the two models are strongly correlated with the amount of high-energy irradiation. Finally, we find that the hydrodynamic effects significantly impact ionisation and heating. Conclusions. The impact of the precise photoionisation treatment provided by Cloudy strongly depends on the system parameters. This suggests that some of the simplifications typically employed in hydrodynamic modelling might lead to systematic errors when studying planetary atmospheres, even at a population-wide level.

Lithium Cepheid V708 Car with an unusual chemical composition

Aims. The purpose of this work is to spectroscopically analyse the classical Cepheid V708 Car. A preliminary check of the spectrum of V708 Car showed that this is a lithium-rich supergiant. We also found that V708 Car has an unusual chemical composition in that the abundances of various elements correlate with their condensation temperatures. We tried to find an explanation of this feature, which is unusual for classical Cepheids. Methods. For the spectroscopic analysis, we used methods based on the assumption of local and non-local thermodynamic equilibrium. Results. We determined the fundamental parameters of our program star V708 Car. This long-period Cepheid has a mass of about 12 M⊙. We derived the abundances of 27 chemical elements in this star. They are clearly correlated with their condensation temperature: the higher the condensation temperature, the lower the abundance (there are exceptions for sodium and barium, however). We explain this peculiar chemical composition of the V708 Car atmosphere by the gas–dust separation in the envelope of this star. A similar mechanism leads to the observed peculiarities of the chemical composition of λ Boo, W Vir, and asymptotic giant branch stars.

Simulations of radiatively cooled magnetic reconnection driven by pulsed power

Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and $P_{1/3}$ multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number $M_A \approx 1.5$ ), supersonic (Sonic Mach number $M_S \approx 4-5$ ) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio $L/\delta \approx 100$ , Lundquist number $S_L \approx 400$ ). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate = $\tau _{\text {cool}}^{-1}/\tau _{H}^{-1}\approx 100$ ); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Finally, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

The evolution of continuum polarization in type II supernovae as a diagnostic of ejecta morphology

The linear polarization of the optical continuum of type II supernovae (SNe), together with its temporal evolution is a promising source of information about the large-scale geometry of their ejecta. To help access this information, we undertook 2D polarized radiative transfer calculations to map the possible landscape of type II SN continuum polarization (Pcont) from 20 to 300 days after explosion. Our simulations were based on crafted 2D axisymmetric ejecta constructed from 1D nonlocal thermodynamic equilibrium time-dependent radiative transfer calculations for the explosion of a red supergiant star. Following the approach used in our previous work on SN 2012aw, we considered a variety of bipolar explosions in which spherical symmetry is broken by material within ~30° of the poles that has a higher kinetic energy (up to a factor of two) and higher 56Ni abundance (up to a factor of about five, allowing for 56Ni at high velocity). Our set of eight 2D ejecta configurations produced considerable diversity in Pcont (λ ~ 7000 Å), although its maximum of 1–4% systematically occurs around the transition to the nebular phase. Before and after this transition, Pcont may be null, constant, rising, or decreasing, which is caused by the complex geometry of the depth-dependent density and ionization and also by optical depth effects. Our modest angle-dependent explosion energy can yield a Pcont of 0.5–1% at early times. Residual optical-depth effects can yield an angle-dependent SN brightness and constant polarization at nebular times. The observed values of Pcont tend to be lower than obtained here. This suggests that more complicated geometries with competing large-scale structures cancel the polarization. Extreme asymmetries seem to be excluded.

M giants with IGRINS

Context. To be able to investigate the chemical history of the entire Milky Way, it is imperative to also study its dust-obscured regions in detail, as this is where most of the mass lies. The Galactic Center is an example of such a region. Due to the intervening dust along the line of sight, near-infrared spectroscopic investigations are necessary to study this region of interest. Aims. The aim of this work is to demonstrate that M giants observed at high spectral resolution in the H- and K-bands (1.5–2.4 μm) can yield useful abundance ratio trends versus metallicity for 21 elements. These elements can then also be studied for heavily dust-obscured regions of the Galaxy, such as the Galactic Center. The abundance ratio trends will be important for further investigation of the Galactic chemical evolution in these regions. Methods. We observed near-infrared spectra of 50 M giants in the solar neighborhood at high signal-to-noise and at a high spectral resolution with the IGRINS spectrometer on the Gemini South telescope. The full H- and K-bands were recorded simultaneously at R = 45 000. Using a manual spectral synthesis method, we determined the fundamental stellar parameters for these stars and derived the stellar abundances for 21 atomic elements, namely, F, Mg, Si, S, Ca, Na, Al, K, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, Ce, Nd, and Yb. We systematically studied useful spectral lines of all these elements in the H- and K-bands. Results. We demonstrate that elements can be analyzed from H- and K-band high-resolution spectra, and we show which spectral lines can be used for an abundance analysis, identifying them line by line. We discuss the 21 abundance ratio trends and compare them with those determined from APOGEE and from the optical Giants in the Local Disk (GILD) sample. From high-resolution H- and K-band spectra, the trends of the heavy elements Cu, Zn, Y, Ce, Nd, and Yb can be retrieved. This opens up the nucleosynthetic channels, including the s-process and the r-process in dust-obscured populations. The [Mn/Fe] versus [Fe/H] trend is shown to be more or less flat at low metallicities, implying that existing non-local thermodynamic equilibrium correction is relevant. Conclusions. With high-resolution near-infrared spectra, it is possible to determine reliable abundance ratio trends versus metallicity for 21 elements, including elements formed in several different nucleosynthetic channels. It is also possible to determine the important neutron-capture elements, both s- and r-dominated elements. This opens up the possibility to study the chemical evolution in detail of dust-obscured regions of the Milky Way, such as the Galactic Center. The M giants are useful bright probes for these regions and for future studies of extra-galactic stellar populations. A careful analysis of high-quality spectra is needed to retrieve all of these elements, which are often from weak and blended lines. A spectral resolution of R ≳ 40 000 is a further quality that helps in deriving precise abundances for this range of elements. In comparison to APOGEE, we can readily obtain the abundances for Cu, Ce, Nd, and Yb from the H-band, demonstrating an advantage of analyzing high-resolution spectra.

A sequence of Type Ib, IIb, II-L, and II-P supernovae from binary-star progenitors with varying initial separations

Over the last decade, evidence has accumulated that massive stars do not typically evolve in isolation but instead follow a tumultuous journey with a companion star on their way to core collapse. While Roche-lobe overflow appears instrumental for the production of a large fraction of Type Ib and Ic supernovae (SNe), variations in the initial orbital period, $P_ init $, of massive interacting binaries may also produce a wide diversity of case B, BC, or C systems, with pre-SN stars endowed from minute to massive H-rich envelopes. Focusing here on the explosion of the primary donor star, originally 12.6\ we used radiation hydrodynamics and nonlocal thermodynamic equilibrium time-dependent radiative transfer to document the gas and radiation properties of such SNe, covering Types Ib, IIb, II-L, and II-P. Variations in init $ are the root cause of the wide diversity of our SN light curves, which present single-peak, double-peak, fast-declining, or plateau-like morphologies in the $V$ band. The different ejecta structures, expansion rates, and relative abundances (e.g., H, He, and can lead to a great deal of diversity in terms of spectral line shapes (absorption versus emission strength and width) and evolution. We emphasize that Halpha is a key tracer of these modulations, and that He is an enduring optical diagnostic for the presence of He. Our grid of simulations fares well against representative Type Ib, IIb, and II-P SNe, but interaction with circumstellar material, which is ignored in this work, is likely at the origin of the tension between our Type II-L SN models and observations (e.g., of SN\,2006Y). Remaining discrepancies in the rise time to bolometric maximum of our models call for a proper account of both small-scale and large-scale structures in core-collapse SN ejecta. Discrepant Type II-P SN models, with a high plateau brightness but small spectral line widths, can be fixed by adopting more compact red-supergiant star progenitors.

Searching for NLTE effects in the high-resolution transmission spectrum of WASP-121 b with cloudy for exoplanets

ABSTRACT Ultrahot Jupiters (UHJs) undergo intense irradiation by their host stars and are expected to experience non-local thermodynamic equilibrium (NLTE) effects in their atmospheres. Such effects are computationally intensive to model but, at the low pressures probed by high-resolution cross-correlation spectroscopy (HRCCS), can significantly impact the formation of spectral lines. The UHJ WASP-121 b exhibits a highly inflated atmosphere, making it ideal for investigating the impact of NLTE effects on its transmission spectrum. Here, we formally introduce cloudy for exoplanets, a cloudy-based modelling code, and use it to generate 1D LTE and NLTE atmospheric models and spectra to analyse archival HARPS WASP-121 b transmission spectra. We assessed the models using two HRCCS methods: (i) Pearson cross-correlation, and (ii) a method that aims to match the average observed line depth for given atmospheric species. All models result in strong detections of Fe i (7.5 < S/N < 10.5). However, the highest S/N model (LTE) does not agree with the best-matching model of the average line depth (NLTE). We also find degeneracy, such that increasing the isothermal temperature and metallicity of the LTE models can produce average line depths similar to cooler, less metal rich NLTE models. Thus, we are unable to conclusively remark on the presence of NLTE effects in the atmosphere of WASP-121 b. We instead highlight the need for standardized metrics in HRCCS that enable robust statistical assessment of complex physical models, e.g. NLTE or 3D effects, that are currently too computationally intensive to include in HRCCS atmospheric retrievals.

An improved methodology for modeling short pulse buried layer x-ray emission spectra

Radiation-hydrodynamic and spectroscopic modeling are important aspects of high energy density experimental design. In this paper, we improve the performance and capabilities over those obtainable with a previous methodology used for simulating x-ray emission spectra from buried layer targets heated by short pulse lasers. The improvement incorporates post-processing HYDRA radiation-hydrodynamic output with a non-local thermodynamic equilibrium atomic-kinetics radiation transport code, Cretin. Each code uses an independent radiation field which allows decoupling HYDRA's radiation group structure from Cretin's spectral output to improve the speed and flexibility of the design methodology. The execution time decreases from 2–3 days to a few hours while the flexibility of the improved methodology allows for performing sensitivity studies including a comparison of steady-state and time-dependent atomic kinetics and differences in the radiation group structure.

Fractionation in young cores: Direct determinations of nitrogen and carbon fractionation in HCN

Nitrogen fractionation is a powerful tracer of the chemical evolution during star and planet formation. It requires robust determinations of the nitrogen fractionation across different evolutionary stages. We aim to determine the $ $N/$ $N and $ $C/$ $C ratios for HCN in six starless and prestellar cores and to compare the results between the direct method using radiative transfer modeling and the indirect double isotope method, assuming a fixed $ $C/$ $C ratio. We present IRAM observations of the HCN $1$-$0$, HCN 3-2, HC15N 1-0 and H13CN $1$-$0$ transitions toward six embedded cores. The $ $N/$ $N ratio was derived using both the indirect double isotope method and directly through non-local thermodynamic equilibrium (NLTE) 1D radiative transfer modeling of the HCN emission. The latter also provides the $ $C/$ $C ratio, which we compared to the local interstellar value. The derived $ $N/$ $N ratios using the indirect method are generally in the range of 300-550. This result could suggest an evolutionary trend in the nitrogen fractionation of HCN between starless cores and later stages of the star formation process. However, the direct method reveals lower fractionation ratios of around sim 250, mainly resulting from a lower $ $C/$ $C ratio in the range sim 20--40, as compared to the local interstellar medium value of 68. This study reveals a significant difference between the nitrogen fractionation ratio in HCN derived using direct and indirect methods. This can influence the interpretation of the chemical evolution and reveal the pitfalls of the indirect double isotope method for fractionation studies. However, the direct method is challenging, as it requires well-constrained source models to produce accurate results. No trend in the nitrogen fractionation of HCN between earlier and later stages of the star formation process is evident when the results of the direct method are considered.

Detailed chemical composition of the globular cluster Sextans A GC-1 on the outskirts of the Local Group

The chemical composition of globular clusters (GCs) across the Local Group provides information on chemical abundance trends. Studying GCs in isolated systems in particular provides us with important initial conditions plausibly unperturbed by mergers and tidal forces from the large Local Group spirals. We present a detailed chemical abundance analysis of Sextans A GC-1. The host galaxy, Sextans A, is a low-surface-brightness dwarf irregular galaxy located on the edge of the Local Group. We derive the dynamical mass of the GC together with the mass-to-light ratio and the abundances of the alpha , Fe-peak, and heavy elements. Abundance ratios were determined from the analysis of an optical integrated-light spectrum of Sextans A GC-1, obtained with UVES on the VLT. We apply non-local thermodynamic equilibrium (NLTE) corrections to Mg, Ca, Ti, Fe, and Ni. The GC appears to be younger and more metal-poor than the majority of the GCs of the Milky Way, with an age of 8.6pm 2.7 Gyr and $ Fe/H dex. The calculated dynamical mass is dyn odot $, which results in an atypically high value of the mass-to-light ratio, 4.35pm 1.40 M$_ odot $/L$_ V odot $. Sextans A GC-1 has varying alpha elements -- the Mg abundance is extremely low, Ca and Ti are solar-scaled or mildly enhanced, and Si is enhanced. The measured values are $ Mg/Fe Ca/Fe Ti/Fe and $ Si/Fe which makes the mean alpha abundance (excluding Mg) to be enhanced $ <Si,Ca,Ti>/Fe NLTE The Fe-peak elements are consistent with scaled-solar or slightly enhanced abundances: $ Cr/Fe Mn/Fe Sc/Fe and $ Ni/Fe The heavy elements measured are Ba, Cu, Zn, and Eu. Ba and Cu have sub-solar abundance ratios ($ Ba/Fe and $ Cu/Fe <-0.343$), while Zn and Eu are consistent with their upper limits being solar-scaled and enhanced, $ Zn/Fe <+0.171$ and $ Eu/Fe The composition of Sextans A GC-1 resembles the overall pattern and behaviour of GCs in the Local Group. The anomalous values are the mass-to-light ratio and the depleted abundance of Mg. There is no definite explanation for such an extreme abundance value. Variations in the initial mass function or the presence of an intermediate-mass black hole might explain the high mass-to-light ratio value.

Observational constraints on the origin of the elements

Context. The chemical evolution history of slow neutron-capture elements in the Milky Way is still a matter of debate, especially in the metal-poor regime ([Fe/H] < −1). Aims. Based on Gaia-ESO spectroscopic data, a recent study investigated the chemical evolution of neutron-capture elements in the regime [Fe/H] > −1. Here, we aim to complement this study down to [Fe/H] = −3, and focus on Ba, Y, and Sr, along with the abundance ratios of [Ba/Y] and [Sr/Y], which give comprehensive views on s-process nucleosynthesis channels. Methods. We measured the local thermodynamic equilibrium (LTE) and non-local thermodynamic equilibrium (NLTE) abundances of Ba, Y, and Sr in 323 Galactic metal-poor stars using high-resolution optical spectra with high signal-to-noise ratios. We used the spectral fitting code TSFitPy together with 1D model atmospheres, using previously determined LTE and NLTE atmospheric parameters. Results. We find that the NLTE effects are on the order of ∼ − 0.1 to ∼0.2 dex, depending on the element. We find that stars enhanced (deficient) in [Ba/Fe] and [Y/Fe] are also enhanced (deficient) in [Sr/Fe], suggesting a common evolution channel for these three elements. We find that the ratio between heavy and light s-process elements [Ba/Y] varies weakly with [Fe/H] even in the metal-poor regime, which is consistent with the behaviour in the metal-rich regime. The [Ba/Y] scatter at a given metallicity is larger than the abundance measurement uncertainties. Homogeneous chemical evolution models with different yield prescriptions are not able to accurately reproduce the [Ba/Y] scatter in the low-[Fe/H] regime. Adopting the stochastic chemical evolution model by Cescutti & Chiappini allows us to reproduce the observed scatter in the abundance pattern of [Ba/Y] and [Ba/Sr]. Based on our observations, we have ruled out the need for an arbitrary scaling of the r-process contribution, as previously suggested by the authors behind the construction of the model. Conclusions. We show how important it is to properly include NLTE effects when measuring chemical abundances, especially in the metal-poor regime. This work demonstrates that the choice of the Galactic chemical evolution model (stochastic versus one-zone) is key when comparing models to observations. Upcoming large-scale spectroscopic surveys such as 4MOST and WEAVE are poised to deliver high-quality data for many thousands of metal-poor stars and this work gives a typical case study of what could be achieved with such surveys in the future.