Non-local Thermodynamic Equilibrium Research Articles

Three-Dimensional Nonlocal Thermodynamic Equilibrium Abundance Analyses of Late-Type Stars

The chemical compositions of stars encode the history of the universe and are thus fundamental for advancing our knowledge of astrophysics and cosmology. However, measurements of elemental abundance ratios, and our interpretations of them, strongly depend on the physical assumptions that dictate the generation of synthetic stellar spectra. Three-dimensional radiation-hydrodynamic (3D RHD) box-in-a-star simulations of stellar atmospheres offer a more realistic representation of surface convection occurring in late-type stars than do traditional one-dimensional (1D) hydrostatic models. As evident from a multitude of observational tests, the coupling of 3D RHD models with line formation in nonlocal thermodynamic equilibrium (non-LTE) today provides a solid foundation for abundance analysis for many elements. This review describes the ongoing and transformational work to advance the state of the art and replace 1D LTE spectrum synthesis with its 3D non-LTE counterpart. In summary: ▪3D and non-LTE effects are intricately coupled, and consistent modeling thereof is necessary for high-precision abundances; such modeling is currently feasible for individual elements in large surveys. Mean 3D (〈3D〉) models are not adequate as substitutes.▪The solar abundance debate is presently dominated by choices and systematic uncertainties that are not specific to 3D non-LTE modeling.▪3D non-LTE abundance corrections have a profound impact on our understanding of FGK-type stars, exoplanets, and the nucleosynthetic origins of the elements.

H_3^+ absorption and emission in local (U)LIRGs with JWST NIRSpec: Evidence for high H_2 ionization rates

We study the $3.4-4.4$ fundamental rovibrational band of H$_3^+$, a key tracer of the ionization of the molecular interstellar medium (ISM), in a sample of 12 local ($d$$<$400\,Mpc) (ultra)luminous infrared galaxies ((U)LIRGs) observed with JWST NIRSpec. The P, Q, and R branches of the band are detected in 13 out of 20 analyzed regions within these (U)LIRGs, which increases the number of extragalactic H$_3^+$ detections by a factor of 6. For the first time in the ISM, the H$_3^+$ band is observed in emission; we detect this emission in three regions. In the remaining ten regions, the band is seen in absorption. The absorptions are produced toward the $3.4-4.4$ hot dust continuum rather than toward the stellar continuum, indicating that they likely originate in clouds associated with the dust continuum source. The H$_3^+$ band is undetected in Seyfert-like (U)LIRGs where the mildly obscured X-ray radiation from the active galactic nuclei might limit the abundance of this molecule. For the detections, the H$_3^+$ abundances, $N($H$_3^+$) H (0.5--5.5)times 10$^ imply relatively high ionization rates, $ H_2 $, of between 3times 10$^ $ and $>$4times 10$^ $, which is lower than the molecular outflow velocities measured using other tracers such as OH\,119 or rotational CO lines. This suggests that H$_3^+$ traces gas close to the outflow-launching sites before it has been fully accelerated. We used nonlocal thermodynamic equilibrium models to investigate the physical conditions of these clouds. In seven out of ten objects, the H$_3^+$ excitation is consistent with inelastic collisions with H$_2$ in warm translucent molecular clouds kin $sim 250--500\,K and H_2 )$sim 10$^ $). In three objects, dominant infrared pumping excitation is required to explain the absorptions from the (3,0) and (2,1) levels of H$_3^+$ detected for the first time in the ISM.

Enabling the Assimilation of CrIS Shortwave Infrared Observations in Global NWP at NOAA. Part I: Background and Methods

Abstract Radiance observations from Earth-observing satellites have a significant positive impact on numerical weather prediction (NWP) forecasts, but some spectral regions are not fully exploited. Observations from hyperspectral infrared (IR) sounders in the longwave region (650-1100 cm−1), for instance, are routinely assimilated in many NWP models, but observations in the shortwave region (2155-2550 cm−1) are not. Each of these regions provides information on the temperature structure of the atmosphere, but the shortwave IR (SWIR) region is considered challenging to assimilate due to noise equivalent delta temperature (NEDT) that is highly variable depending on scene brightness temperature and to phenomena that are difficult to model, like non-Local Thermodynamic Equilibrium (NLTE) and solar reflectance. With recent advances in small-satellite technology, SWIR temperature sounders may provide an agile and cost-effective complement to the current constellation of IR sounders. Therefore, a better understanding of the use and impact of SWIR observations in data assimilation for NWP is warranted. In part one of this study, as presented here, the amount of unique information (as determined by Empirical Orthogonal Decomposition (EOD)) made available to a data assimilation system by Cross-track Infrared Sounder (CrIS) SWIR observations is reviewed, recent advancements to the Community Radiative Transfer Model (CRTM) for the simulation of CrIS shortwave radiances are tested, and enhancements to NOAA’s Global Data Assimilation System (GDAS) for the assimilation of CrIS SWIR observations are implemented and evaluated. Part two of this study, which seeks to assess the value of assimilating shortwave IR observations in global NWP, is also introduced.

Probing atmospheric escape through metastable He I triplet lines in 15 exoplanets observed with SPIRou

For several years, the metastable helium triplet line has been successfully used as a tracer to probe atmospheric escape in transiting exoplanets. This absorption in the near-infrared (1083.3 nm) can be observed from the ground using high-resolution spectroscopy, providing new constraints on the mass-loss rate and the temperature characterizing the upper atmosphere of close-in exoplanets. The aim of this work is to search for the He triplet signature in 15 transiting exoplanets – ranging from super-Earths to ultrahot Jupiters – observed with SPIRou, a high-resolution (R ~ 70 000) near-infrared spectropolarimeter at the CFHT, in order to bring new constraints or to improve existing ones regarding atmospheric escape through a homogeneous study. We developed a full data processing and analysis pipeline to correct for the residual telluric and stellar contributions. We then used two different 1D models based on the Parkerwind equations and nonlocal thermodynamic equilibrium (NLTE) radiative transfer to interpret the observational results. We confirm published He triplet detections for HAT-P-11 b, HD 189733 b, and WASP-69 b. We tentatively detect the signature of escaping He in HD 209458 b, GJ 3470 b, and WASP-76 b. We report new constraints on the mass-loss rate and temperature for our three detections and set upper limits for the tentative and nondetections. We notably report improved constraints on the mass-loss rate and temperature of the escaping gas for TOI-1807 b, and report a nondetection for the debated atmospheric escape in GJ 1214 b. We also conducted the first search for the He signature in GJ 486 b since its discovery and report a nondetection of the He triplet. Finally, we studied the impact of important model assumptions on our retrieved parameters, notably the limitations of 1D models and the influence of the H/He ratio on the derived constraints.

Simultaneous Spectropolarimetric Observations in the Hα and Ca ii 8662 Å Lines of an Active Region

We present spectropolarimetric observations of an active region recorded simultaneously in the Hα and the Ca ii 8662 Å lines. The sunspot exhibits multiple structures, including a light bridge and a region where the Ca ii 8662 Å line core is in emission. Correspondingly, the Hα line core image displays brightening in the emission region, with the spectral profiles showing elevated line cores. The stratification of the line-of-sight magnetic field is inferred through nonlocal thermodynamic equilibrium multiline inversions of the Ca ii 8662 Å line and the weak field approximation (WFA) over the Hα line. The field strength inferred from the Hα line core is consistently smaller than that inferred from inversions at logτ500 = −4.5. However, the study finds no correlation between the WFA over the core of the Hα line and that inferred from inversions at logτ500 = −4.5. In regions exhibiting emission features, the morphology of the magnetic field at logτ500 = −4.5 resembles that at logτ500 = −1, with slightly higher or comparable field strengths. The magnetic field morphology inferred from the core of the Hα line is also similar to that inferred from the full spectral range of the Hα line in the emission region. The field strength inferred in the light bridge at logτ500 = −1 is smaller than the surrounding umbral regions and comparable at logτ500 = −4.5. Similarly, the field strength inferred in the light bridge from the WFA over the Hα line appears lower compared to the surrounding umbral regions.

Final Moments. II. Observational Properties and Physical Modeling of Circumstellar-material-interacting Type II Supernovae

We present ultraviolet/optical/near-infrared observations and modeling of Type II supernovae (SNe II) whose early time (δ t < 2 days) spectra show transient, narrow emission lines from shock ionization of confined (r < 1015 cm) circumstellar material (CSM). The observed electron-scattering broadened line profiles (i.e., IIn-like) of H i, He i/ii, C iv, and N iii/iv/v from the CSM persist on a characteristic timescale (t IIn) that marks a transition to a lower-density CSM and the emergence of Doppler-broadened features from the fast-moving SN ejecta. Our sample, the largest to date, consists of 39 SNe with early time IIn-like features in addition to 35 “comparison” SNe with no evidence of early time IIn-like features, all with ultraviolet observations. The total sample includes 50 unpublished objects with a total of 474 previously unpublished spectra and 50 multiband light curves, collected primarily through the Young Supernova Experiment and Global Supernova Project collaborations. For all sample objects, we find a significant correlation between peak ultraviolet brightness and both t IIn and the rise time, as well as evidence for enhanced peak luminosities in SNe II with IIn-like features. We quantify mass-loss rates and CSM density for the sample through the matching of peak multiband absolute magnitudes, rise times, t IIn, and optical SN spectra with a grid of radiation hydrodynamics and non-local thermodynamic equilibrium radiative-transfer simulations. For our grid of models, all with the same underlying explosion, there is a trend between the duration of the electron-scattering broadened line profiles and inferred mass-loss rate: tIIn≈3.8[Ṁ/ (0.01 M ⊙ yr−1)] days.

Abundances of neutron-capture elements in selected solar-type stars

The primary objective of this study is to accurately determine the abundances of Cu, Sr, Y, Zr, Ba, La, and Ce in selected solar-type stars. This will allow us to establish observational abundance--metallicity and abundance--age relations and to explore the reasons for the excess of Ba compared to other $s$-elements in younger solar-type stars. The chosen $s$-process elements are critical diagnostics for understanding the chemical evolution of our Galaxy. We analysed HARPS spectra with a high resolution ($R$ = 115\,000) and high signal-to-noise ratio (close to 100) of main-sequence solar-type FGK stars with metallicities from $-0.15$ to +0.35 dex and ages from 2 to 14 Gyr using one-dimensional (1D) local thermodynamic equilibrium (LTE) synthesis and MARCS atmospheric models. In the procedure of fitting synthetic to observed line profiles, the free parameters included abundance and microturbulent and macroturbulent velocity. The macroturbulent velocity can substantially compensate for non-local thermodynamic equilibrium (NLTE) effects in the line core. The resulting elemental abundance X/H increases with metallicity and age for solar-type stars. The ratio of the abundances of $s$-process elements s/Fe increases with decreasing metallicity and age, while the Cu/Fe ratio increases with both metallicity and age. These observed trends agree well with published observational data and with predictions from Galactic chemical evolution (GCE) models. A small Ba/Fe enhancement of 0.08 ± 0.08 dex has been detected in seven younger stars with an average age of $2.8 0.6$ Gyr. Compared to the abundances of other $s$-process elements Ba/Fe is 0.07 and 0.08 dex higher than La and Ce on average, respectively. Furthermore, we find that the Ba/Fe ratio increases with increasing chromospheric activity. The average Ba/Fe for the three most active stars is $0.15 0.10$ dex higher than that of the other stars. Chromospheric activity, characterised by stronger magnetic fields found in active regions such as pores, spots, plages, and networks, can significantly alter the physical conditions in the formation layers of the Ba lines. Our primary conclusion is that to account for the observed excess of Ba/Fe abundance in younger stars, it is essential to use more complex atmospheric models that incorporate magnetic structures.

Discovery of an Extremely r-process-enhanced Thin-disk Star with [Eu/H] = +0.78

Highly r-process-enhanced (RPE) stars are rare and usually metal poor ([Fe/H] < −1.0), and they mainly populate the Milky Way halo and dwarf galaxies. This study presents the discovery of a relatively bright (V = 12.72), highly RPE (r-II) star ([Eu/Fe] = +1.32, [Ba/Eu] = −0.95), LAMOST J020623.21+494127.9. This star was selected from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope medium-resolution (R ∼ 7500) spectroscopic survey; follow-up high-resolution (R ∼ 25,000) observations were conducted with the High Optical Resolution Spectrograph installed on the Gran Telescopio Canarias. The stellar parameters (T eff = 4130 K, logg = 1.52, [Fe/H] = −0.54, ξ = 1.80 km s−1) have been inferred taking into account nonlocal thermodynamic equilibrium effects. The abundances of [Ce/Fe], [Pr/Fe], and [Nd/Fe] are +0.19, +0.65, and +0.64, respectively, relatively low compared to the Solar r-process pattern normalized to Eu. This star has a high metallicity ([Fe/H] = −0.54) compared to most other highly RPE stars and has the highest measured abundance ratio of Eu to H ([Eu/H] = +0.78). It is classified as a thin-disk star based on its kinematics and does not appear to belong to any known stream or dwarf galaxy.

Understanding the deficiency in inertial confinement fusion hohlraum x-ray flux predictions using experiments at the National Ignition Facility.

The predicted implosion performance of deuterium-tritium fuel capsules in indirect-drive inertial confinement fusion experiments relies on precise calculations of the x-ray drive in laser-heated cavities (hohlraums). This requires accurate, spectrally dependent simulations of laser to x-ray conversion efficiencies and x-ray absorption losses to the hohlraum wall. A set of National Ignition Facility experiments have identified a cause for the long-standing hohlraum "drive deficit" as the overprediction of gold emission at ∼2.5 keV in nonlocal thermodynamic equilibrium coronal plasma regions within the hohlraum. Reducing the emission and absorption opacity in this spectral region by ∼20% brings simulations into agreement with measured x-ray fluxes and spectra.

New rotational rate coefficients computation of the linear NaC3N(X1Σ+) by collision with He(1S) and astrophysical implication

ABSTRACT For a better understanding of the physico-chemistry in the interstellar medium, collisional data are needed. In this work, we provide rate coefficients for the recently detected by Cabezas and collaborators in 2023 sodium cyanoacetylide molecule NaC$_{3}$N($X^{1}\Sigma ^{+}$) induced by collisions with He. A new two-dimensional potential energy surface (2D-PES) is derived by adopting the high-level theory of the explicitly correlated coupled cluster with single, double, and perturbative triple excitations (CCSD(T)-F12) in conjunction with the aug-cc-pVTZ basis sets. Our 2D-PES presents two minima: the global one with a well depth of $-86.58$ cm$^{-1}$ and the second well with a $-47.01$ cm$^{-1}$ depth. This was expected by the Natural Bond Orbital analysis showing two electrophilic zones located on Na and N linked Carbon with respective atomic charges : $+0.96575$ and $+0.28910$. This PES was used to calculate inelastic cross-sections within the framework of the close-coupling theory for total energies $\le 170$ cm$^{-1}$ and $J\le 22$. Using the Maxwell–Boltzmann distribution of kinetic energies, these cross-sections were integrated to generate the collisional (de)-excitation rate coefficients for temperatures below 30 K. In order to estimate the impact of the computed collisional rates, a non-local thermodynamical equilibrium radiative transfer calculation was performed. The new collisional data can allow accurate NaC$_{3}$N abundance determination in the interstellar medium.

OCCASO

Context. 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. Aims. 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. Methods. We used the high-resolution spectra (R &gt; 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. Results. We measured LTE abundances for 21 elements: α (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 &lt; 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 &gt; Ga).

Diagnostics of 3D explosion asymmetries of stripped-envelope supernovae by nebular line profiles

ABSTRACT Understanding the explosion mechanism and hydrodynamic evolution of core-collapse supernovae (SNe) is a long-standing quest in astronomy. The asymmetries caused by the explosion are encoded into the line profiles which appear in the nebular phase of the SN evolution – with particularly clean imprints in He star explosions. Here, we carry out nine different supernova simulations of He-core progenitors, exploding them in 3D with parametrically varied neutrino luminosities using the prometheus-hotb code, hydrodynamically evolving the models to the homologous phase. We then compute nebular phase spectra with the 3D Non-Local Thermodynamic Equilibrium spectral synthesis code extrass (EXplosive TRAnsient Spectral Simulator). We study how line widths and shifts depend on progenitor mass, explosion energy, and viewing angle. We compare the predicted line profile properties against a large set of Type Ib observations, and discuss the degree to which current neutrino-driven explosions can match observationally inferred asymmetries. With self-consistent 3D modelling – circumventing the difficulties of representing $^{56}$Ni mixing and clumping accurately in 1D models – we find that neither low-mass He cores exploding with high energies nor high-mass cores exploding with low energies contribute to the Type Ib SN population. Models which have line profile widths in agreement with this population give sufficiently large centroid shifts for calcium emission lines. Calcium is more strongly affected by explosion asymmetries connected to the neutron star kicks than oxygen and magnesium. Lastly, we turn to the near-infrared spectra from our models to investigate the potential of using this regime to look for the presence of He in the nebular phase.

Carbon monoxide cooling in radiative transfer modelling of supernovae

ABSTRACT Carbon monoxide (CO) emission has been observed in a number of core-collapse supernovae (SNe) and is known to be an important coolant at late times. We have implemented a chemical reaction network in the radiative-transfer code cmfgen to investigate the formation of CO and its impact on SN ejecta. We calculate two 1D SN models with and without CO: a BSG explosion model at one nebular epoch and a full time-sequence (50–300 d) for a red supergiant explosion. In both models, CO forms at nebular times in the dense, inner regions at velocities &amp;lt;2000 km s$^{-1}$ where line emission from CO can dominate the cooling and reduce the local temperature by as much as a factor of 2, weakening emission lines and causing the optical light curve to fade faster. That energy is instead emitted in CO bands, primarily the fundamental band at $\sim 4.5 \mathrm{\mu m}$, which accounts for up to 20 per cent of the total luminosity at late times. However, the non-monotonic nature of the CO cooling function can cause numerical difficulties and introduce multiple temperature solutions. This issue is compounded by the sensitivity of the CO abundance to a few reaction rates, many of which have large uncertainties or disparate values across literature sources. Our results also suggest that, in many SNe, CO-level populations are far from their local thermodynamic equilibrium values. Unfortunately, accurate collisional data, necessary to compute non-local thermodynamic equilibrium populations, are limited to a few transitions.

On the Pair-instability Supernova Origin of J1010+2358* * Based on observations made with the Very Large Telescope (VLT) of the ESO at the La Silla Paranal Observatory under the program ID 112.26WJ.001

The first (Population III) stars formed only out of H and He and were likely more massive than present-day stars. Massive Population III stars in the range 140–260 M ⊙ are predicted to end their lives as pair-instability supernovae (PISNe), enriching the environment with a unique abundance pattern, with high ratios of odd to even elements. Recently, the most promising candidate for a pure descendant of a zero-metallicity massive PISN (260 M ⊙) was discovered by the LAMOST survey, the star J1010+2358. However, key elements to verify the high PISN contribution, C and Al, were missing from the analysis. To rectify this, we obtained and analyzed a high-resolution Very Large Telescope/UVES spectrum, correcting for 3D and/or non-local thermodynamic equilibrium effects. Our measurements of both C and Al give much higher values (∼1 dex) than expected from a 260 M ⊙ PISN. Furthermore, we find significant discrepancies with the previous analysis and therefore a much less pronounced odd–even pattern. Our results show that J1010+2358 cannot be a pure descendant of a 260 M ⊙ PISN. Instead, we find that the best-fit model consists of a 13 M ⊙ Population II core-collapse supernova combined with a Population III supernova. Alternative, less favored solutions ( χ2/χbest2≈2.3 ) include a 50% contribution from a 260 M ⊙ PISN or a 40% contribution from a Population III Type Ia supernova. Ultimately, J1010+2358 is certainly a unique star giving insights into the earliest chemical enrichment; however, this star is not a pure PISN descendant.

Bias versus variance when fitting multi-species molecular lines with a non-LTE radiative transfer model

Context. Robust radiative transfer techniques are requisite for efficiently extracting the physical and chemical information from molecular rotational lines. Aims. We study several hypotheses that enable robust estimations of the column densities and physical conditions when fitting one or two transitions per molecular species. We study the extent to which simplifying assumptions aimed at reducing the complexity of the problem introduce estimation biases and how to detect them. Methods. We focus on the CO and HCO+ isotopologues and analyze maps of a 50 square arcminutes field. We used the RADEX escape probability model to solve the statistical equilibrium equations and compute the emerging line profiles, assuming that all species coexist. Depending on the considered set of species, we also fixed the abundance ratio between some species and explored different values. We proposed a maximum likelihood estimator to infer the physical conditions and considered the effect of both the thermal noise and calibration uncertainty. We analyzed any potential biases induced by model misspecifications by comparing the results on the actual data for several sets of species and confirmed with Monte Carlo simulations. The variance of the estimations and the efficiency of the estimator were studied based on the Cramér-Rao lower bound. Results. Column densities can be estimated with 30% accuracy, while the best estimations of the volume density are found to be within a factor of two. Under the chosen model framework, the peak 12CO (1 – 0) is useful for constraining the kinetic temperature. The thermal pressure is better and more robustly estimated than the volume density and kinetic temperature separately. Analyzing CO and HCO+ isotopologues and fitting the full line profile are recommended practices with respect to detecting possible biases. Conclusions. Combining a non-local thermodynamic equilibrium model with a rigorous analysis of the accuracy allows us to obtain an efficient estimator and identify where the model is misspecified. We note that other combinations of molecular lines could be studied in the future.

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 $&lt;$ -- 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 $&gt;$ 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 $&gt;$ Fe/H $&gt;$--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.

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

Context. 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. Aims. We aim to determine the 14N/15N and 12C/13C 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 12C/13C ratio. Methods. 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 14N/15N 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 12C/13C ratio, which we compared to the local interstellar value. Results. The derived 14N/15N 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 ~250, mainly resulting from a lower 12C/13C ratio in the range ~20–40, as compared to the local interstellar medium value of 68. Conclusions. 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.