non-LTE Model Research Articles

The solar beryllium abundance revisited with 3D non-LTE models

The present-day abundance of beryllium in the solar atmosphere provides clues about mixing mechanisms within stellar interiors. However, abundance determinations based on the Be II313.107 nm line are prone to systematic errors due to imperfect model spectra. These errors arise from missing continuous opacity in the UV, a significant unidentified blend at 313.102 nm, departures from local thermodynamic equilibrium (LTE), and microturbulence and macroturbulence fudge parameters associated with one-dimensional (1D) hydrostatic model atmospheres. Although these factors have been discussed in the literature, no study has yet accounted for all of them simultaneously. To address this, we present 3D non-LTE calculations for neutral and ionised beryllium in the Sun. We used these models to derive the present-day solar beryllium abundance, calibrating the missing opacity on high resolution solar irradiance data and the unidentified blend on the centre-to-limb variation. We find a surface abundance of 1.21 ± 0.05 dex, which is significantly lower than the value of 1.38 dex that has been commonly adopted since 2004. Taking the initial abundance via CI chondrites, our result implies that beryllium has been depleted from the surface by an extra 0.11 ± 0.06 dex, or 22 ± 11%, on top of any effects of atomic diffusion. This is in tension with standard solar models, which predict negligible depletion, as well as with contemporary solar models that have extra mixing calibrated on the abundances of helium and lithium, which predict excessive depletion. These discrepancies highlight the need for further improvements to the physics in solar and stellar models.

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.

Imaging the Jet of MWC 349A with Resolved Radio Recombination Line Emission from ALMA

Jets and disk winds arise from materials with excess angular momentum ejected from the accretion disks in forming stars. How these structures are launched and how they impact the gas within the innermost regions of these objects remains poorly understood. MWC349A is a massive star that has a circumstellar disk that rotates in accord with Kepler’s Law, with an ionized wind and a high-velocity jet launched from the disk surface. The strongly maser-amplified emission of hydrogen radio recombination lines (RRLs) observed toward this system provides a comprehensive picture of its ionized environment with exquisite detail. In this Letter, we present Atacama Large Millimeter/submillimeter Array observations of the H26α RRL and continuum emission obtained with the highest angular resolution ever used toward this source (beam of ∼0.″02). The maser RRL emission is resolved for the first time and clearly delineates the ionized disk, wind, and jet. We analyzed the RRL data cubes with the 3D non-LTE radiative transfer model MORELI, confirming that the jet is poorly collimated. We found that the jet orientation is closer to the rotation axis of the system than derived from spatially unresolved data. This study confirms that hydrogen RRL masers are powerful probes of the physical structure and kinematics of the innermost ionized material around massive stars.

3D non-local thermodynamic equilibrium magnesium abundances reveal a distinct halo population

Magnesium is one of the most important elements in stellar physics as an electron donor; in Galactic archaeology, magnesium serves to distinguish different stellar populations. However, previous studies of Mg I and Mg II lines in metal-poor benchmark stars indicate that magnesium abundances inferred from one-dimensional (1D), hydrostatic models of stellar atmospheres, both with and without the local thermodynamic equilibrium (LTE) approximation, can be problematic. Here, we present three-dimensional (3D) non-LTE calculations for magnesium in FG-type dwarfs and provide corrections for 1D-LTE abundances. 3D non-LTE corrections reduce the ionisation imbalances in the benchmark metal-poor stars HD84937 and HD140283 from −0.16 dex and −0.27 dex in 1D LTE to just −0.02 dex and −0.09 dex, respectively. We applied our abundance corrections to 1D LTE literature results for stars in the thin disc, thick disc, α-rich halo, and α-poor halo. We observed that 3D non-LTE results had a richer substructure in [Mg/Fe] − [Fe/H] in the α-poor halo, revealing two sub-populations at the metal-rich end. These two sub-populations also differ in kinematics, supporting the astro-physical origin of the separation. While the more magnesium-poor sub-population is likely to be debris from a massive accreted galaxy, Gaia-Enceladus, the other sub-population may be related to a previously identified group of stars, called Eos. The additional separation in [Mg/Fe] suggests that previous Mg abundance measurements may have been imprecise due to the 1D and LTE approximations, highlighting the importance of 3D non-LTE modelling.

Chemical inventory of the envelope of the Class I protostar L1551 IRS 5

Episodic accretion in protostars leads to luminosity outbursts that end up heating their surroundings. This rise in temperature pushes the snow lines back, enabling the desorption of chemical species from dust grain surfaces, which may significantly alter the chemical history of the accreting envelope. However, a limited number of extensive chemical surveys of eruptive young stars have been performed thus far. In the present study, we carry out a large spectral survey of the binary Class I protostar L1551 IRS 5, known to be a FUor-like object, in the 3 mm and 2mm bands with the IRAM-30m telescope. As a result, we detected more than 400 molecular lines. The source displays a great chemical richness with the detection of 75 species, including isotopologues. Among these species, there are 13 hydrocarbons, 25 N-bearing species, 30 O-bearing species, 15 S-bearing species, 12 deuterated molecules, and a total of 10 complex organic molecules (l-C4H2, CH3CCH, CH2DCCH, CH3CHO, CH3CN, CH3OCH3, CH3OCHO, CH3OH, CH2DOH, and HC5N). With the help of local thermodynamic equilibrium (LTE) and non-LTE models, we determined the column densities of most molecules as well as excitation and kinetic temperatures. While most of those molecules trace the cold envelope (≲20 K), the OCS and CH3OH emission arise from the warm (>100 K) innermost (<2″) regions. We compared the chemical inventory of L1551 IRS 5 and its column density ratios, including isotopic ratios, with other protostellar sources. A broad chemical diversity is seen among Class I objects. More observations with both single-dish telescopes and interferometers are needed to characterize the diversity in a larger sample of protostars, while more astrochemical models would help explain this diversity, in addition to the impact of luminosity outbursts on the chemistry of protostellar envelopes.

Titanium abundances in late-type stars

Context. Rapidly growing datasets from stellar spectroscopic surveys are providing unprecedented opportunities to analyse the chemical evolution history of our Galaxy. However, spectral analysis requires accurate modelling of synthetic stellar spectra for late-type stars, for which the assumption of local thermodynamic equilibrium (LTE) has been shown to be insufficient in many cases. Errors associated with LTE can be particularly large for Ti I, which is susceptible to over-ionisation, particularly in metal-poor stars. Aims. The aims of this work are to study and quantify the 1D non-LTE effects on titanium abundances across the Hertzsprung-Russell diagram for a large sample of stars. Methods. A large grid of departure coefficients, βν, were computed on standard MARCS model atmospheres. The grid extends from 3000 K to 8000 K in Teff, −0.5 dex to +5.5 dex in log g, and −5.0 to +1.0 in [Fe/H], with non-LTE effects in this grid reaching up to 0.4 dex. This was used to compute abundance corrections that were subsequently applied to the LTE abundances of over 70 000 stars selected from the GALAH survey in addition to a smaller sample of literature Keck data for metal-poor dwarfs. Results. The non-LTE effects grow towards lower [Fe/H], lower log g, and higher Teff, with a minimum and maximum ΔA(Ti)Ti I of 0.02 and 0.19 in the GALAH sample. For metal-poor giants, the non-LTE modelling reduces the average ionisation imbalance (ΔI−II) from −0.11 dex to −0.01 dex at [Fe/H] = −1.7, and the enhancement in titanium abundances from Ti I lines results in a [Ti/Fe] versus [Fe/H] trend that more closely resembles the behaviour of Ti II at low metallicities. At higher metallicities, the results are limited by the precision of the GALAH DR3 LTE abundances and the effects are within the errors. For the most metal-poor dwarfs from the Keck sample, the average ionisation imbalance increases from −0.1 dex to +0.2 dex, a shortcoming that is consistent with previous 1D non-LTE studies and which we speculate could be related to 3D effects. Conclusions. Non-LTE effects on titanium abundances are significant. Neglecting them may alter our understanding of Galactic chemical evolution. We have made our grid of departure coefficients publicly available, with the caveat that the Ti abundances of metal-poor dwarfs need further study in 3D non-LTE.

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.

Maser Activity of Organic Molecules toward Sgr B2(N)

At centimeter wavelengths, single-dish observations have suggested that the Sagittarius (Sgr) B2 molecular cloud at the Galactic Center hosts weak maser emission from several organic molecules, including CH2NH, HNCNH, and HCOOCH3. However, the lack of spatial distribution information on these new maser species has prevented us from assessing the excitation conditions of the maser emission as well as their pumping mechanisms. Here, we present a mapping study toward Sgr B2 north (N) to locate the region where the complex maser emission originates. We report the first detection of the Class I methanol (CH3OH) maser at 84 GHz and the first interferometric map of the methanimine (CH2NH) maser at 5.29 GHz toward this region. In addition, we present a tool for modeling and fitting the unsaturated molecular maser signals with non-LTE radiative transfer models and Bayesian analysis using the Markov Chain Monte Carlo approach. These enable us to quantitatively assess the observed spectral profiles. The results suggest a two-chain-clump model for explaining the intense CH3OH Class I maser emission toward a region with low continuum background radiation. By comparing the spatial origin and extent of maser emission from several molecular species, we find that the 5.29 GHz CH2NH maser has a close spatial relationship with the 84 GHz CH3OH Class I masers. This relationship serves as observational evidence to suggest a similar collisional pumping mechanism for these maser transitions.

Molecular isotopologue measurements toward super star clusters and the relation to their ages in NGC 253 with ALCHEMI

Context. Determining the evolution of the CNO isotopes in the interstellar medium (ISM) of starburst galaxies can yield important constraints on the ages of super star clusters (SSCs), or on other aspects and factors contributing to their evolution, such as the initial mass function (IMF). Due to the time-dependent nature of the abundances of isotopes within the ISM – as they are supplied from processes such as nucleosynthesis or chemical fractionation –, this provides the opportunity to test whether or not isotope ratios trace the ages of highly star-forming regions, such as SSCs. Aims. The goal of this study is to investigate whether the isotopic variations in SSC regions within NGC 253 are correlated with their different ages as derived from stellar population modelling. Methods. We measured abundance ratios of CO, HCN, and HCO+ isotopologues in six regions containing SSCs within NGC 253 using high-spatial-resolution (1.6″, ∼28 pc) data from the ALCHEMI (ALma Comprehensive High-resolution Extragalactic Molecular Inventory) ALMA Large program. We then analysed these ratios using RADEX radiative transfer modelling, with the parameter space sampled using the nested sampling Monte Carlo algorithm MLFriends. These abundance ratios were then compared to ages predicted in each region via the fitting of observed star-formation tracers (such as Brγ) to Starburst99 starburst stellar population evolution models. Results. We determined the isotopic column density ratios across multiple regions of SSC activity in NGC 253 using non-LTE radiative transfer modelling. We do not find any significant trend with age for the CO and HCN isotopologue ratios on timescales of the ages of the SSC* regions observed. However, HCO+ may show a correlation with age over these timescales in 12C/13C. Conclusions. The driving factors of these ratios within SSCs could be the IMF or fractionation effects. To further probe these effects in SSCs over time, a larger sample of SSCs must be observed spanning a larger age range.

The ALMA-QUARKS Survey. II. The ACA 1.3 mm Continuum Source Catalog and the Assembly of Dense Gas in Massive Star-Forming Clumps

Leveraging the high resolution, sensitivity, and wide frequency coverage of the Atacama Large Millimeter/submillimeter Array (ALMA), the QUARKS survey, standing for “Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures”, is observing 139 massive star-forming clumps at ALMA Band 6 (λ ∼ 1.3 mm). This paper introduces the Atacama Compact Array (ACA) 7 m data of the QUARKS survey, describing the ACA observations and data reduction. Combining multi-wavelength data, we provide the first edition of QUARKS atlas, offering insights into the multiscale and multiphase interstellar medium in high-mass star formation. The ACA 1.3 mm catalog includes 207 continuum sources that are called ACA sources. Their gas kinetic temperatures are estimated using three formaldehyde transitions with a non-LTE radiation transfer model, and the mass and density are derived from a dust emission model. The ACA sources are massive (16–84 percentile values of 6–160 M ⊙), gravity-dominated (M ∝ R 1.1) fragments within massive clumps, with supersonic turbulence ( M>1 ) and embedded star-forming protoclusters. We find a linear correlation between the masses of the fragments and the massive clumps, with a ratio of 6% between the two. When considering fragments as representative of dense gas, the ratio indicates a dense gas fraction (DGF) of 6%, although with a wide scatter ranging from 1% to 10%. If we consider the QUARKS massive clumps to be what is observed at various scales, then the size-independent DGF indicates a self-similar fragmentation or collapsing mode in protocluster formation. With the ACA data over four orders of magnitude of luminosity-to-mass ratio (L/M), we find that the DGF increases significantly with L/M, which indicates clump evolutionary stage. We observed a limited fragmentation at the subclump scale, which can be explained by a dynamic global collapse process.

Discovery of a new blue large-amplitude pulsator in the SkyMapper DR2: SMSS J184506.82−300804.7

ABSTRACT We report the discovery of a new Blue Large-Amplitude Pulsator (BLAP), SMSS J184506−300804 (SMSS-BLAP-1) in Data Release 2 of the SkyMapper Southern Sky Survey. We conduct high-cadence photometric observations in the u band to confirm a periodic modulation of the light curve. SMSS-BLAP-1 has a ∼19-min pulsation period with an amplitude of 0.2 mag in u band, and is similar to the classical BLAPs found by OGLE. From spectroscopic observations with the Wide-Field Spectrograph on the Australian National University 2.3m telescope, we confirm it as a low-gravity BLAP: best-fitting parameters from the non-LTE tlusty model are estimated as Teff = 29 020$^{+193}_{-34}$ K, logg = 4.661$^{+0.008}_{-0.143}$ (cm s−2), and logn(He)/n(H) = −2.722$^{+0.057}_{-0.074}$ dex. However, our BLAP exhibits a very He-deficient atmosphere compared to both low- and high-gravity BLAPs, which have logn(He)/n(H) in the range −0.41∼−2.4.

3D non-LTE modeling of the stellar center-to-limb variation for transmission spectroscopy studies

Context. Transmission spectroscopy is one of the most powerful techniques used to characterize transiting exoplanets, since it allows for the abundance of the atomic and molecular species in the planetary atmosphere to be measured. However, stellar lines may bias the determination of such abundances if their center-to-limb variations (CLVs) are not properly accounted for. Aims. This paper aims to show that three-dimensional (3D) radiation hydrodynamic models and the assumption of non-local ther-modynamic equilibrium (non-LTE) line formation are required for an accurate modeling of the stellar CLV of the Na I D1 and K I resonance lines on transmission spectra. Methods. We modeled the CLV of the Na I D1 and K I resonance lines in the Sun with 3D non-LTE radiative transfer. The synthetic spectra were compared to solar observations with high spatial and spectral resolution, including new data collected with the CRISP instrument at the Swedish 1-m Solar Telescope between µ = 0.1 and µ = 1.0. Results. Our 3D non-LTE modeling of the Na I D1 resonance line at 5896 Å and the K I 7699 Å resonance line in the Sun is in good agreement with the observed CLV in the solar spectrum. Moreover, the simulated CLV curve for a Jupiter-Sun system inferred with a 3D non-LTE analysis shows significant differences from the one obtained from a 1D atmosphere. The latter does indeed tend to overestimate the amplitude of the transmission curve by a factor that is on the same order of magnitude as a planetary absorption depth (i.e., up to 0.2%). Conclusions. This work highlights the fact that to correctly characterize exoplanetary atmospheres, 3D non-LTE synthetic spectra ought to be used to estimate the stellar CLV effect in transmission spectra of solar-like planet hosts. Moreover, since different spectral lines show different CLV curves for the same geometry of the planet-star system, it is fundamental to model the CLV individually for each line of interest. The work will be extended to other lines and FGK-type stars, allowing for synthetic high-resolution spectra to mitigate the stellar contamination of low-resolution planetary spectra, for example, those drawn from JWST.

ZPEKTR: A code for spectral synthesis of fast-rotating stars

Estimating the physical states of the surfaces of fast-rotating stars is challenging due to several intrinsic processes, which include radiative flux inhomogeneities on the photosphere induced by rotation and circumstellar signatures in their spectra. The analysis of their spectra ultimately requires the use of synthetic grids of spectra accounting for all these physical processes. In this paper, we present the `von ZeiPEl's code for gravity darKening specTRal synthesis' (ZPEKTR) code, which is designed to perform the spectral synthesis of fast-rotating stars, accounting for gravity darkening, limb-darkening effects in the continuum and geometrical deformation induced by fast rotation. We consider colatitudinal temperature and surface-gravity variations, assuming both the classical prescription developed by von Zeipel and the new formulation by Espinosa-Lara. The code runs either with a rectangular or a triangular mesh on the stellar surface. We compare the temperature and gravitational distribution as a function of the stellar latitude arising from both models. The line profiles of 4388, 4471, 4922, and 6678 produced with both formalisms are compared at three different rotation rates and illustrate differences in shape and central intensity. We also illustrate the fittings of 31 line spectra of classical Be stars averaged from the Be Stars Observation Survey (BeSOS) database and make a comparison among their apparent physical parameters and ages determined from plane-parallel non-local thermodynamical equilibrium (non-LTE) models and parameters determined from classical von Zeipel models, finding a displacement of more evolved objects towards the zero-age main sequence. We also compare the distributions of projected rotation velocities of these objects obtained with and without the inclusion of gravity-darkening effects with ZPEKTR. We observe a shift of the histogram of rotation velocities calculated accounting for effects of gravity darkening concerning rotation velocities obtained through the fittings with classical plane-parallel non-LTE models. We show that models that do not account for gravity darkening can underestimate the rotation velocity, because the stellar latitudes that contribute the higher velocities are those in the equator with the least radiative flux. We envisage near-future improvements to the code, such as the inclusion of differential rotation and treatment of tidal forces in binary stellar systems.

A survey of molecular line emission towards Herbig Be star V645 Cyg

ABSTRACT We present a survey of molecular line emission towards the molecular cloud surrounding Herbig Be star V645 Cyg. The survey was performed with the 20-m Onsala space telescope at 3 and 4 mm. We detected emission lines of 33 molecules and their isotopologues from diatomic molecules to four COMs up to seven atoms. Using detected lines, we estimated molecular column densities and abundances relative to molecular hydrogen in local thermodynamic equilibrium (LTE) approximation for all molecules except for methanol, for which we obtained physical parameters using a non-LTE model. Moreover, in the basement of the non-LTE model of methanol line emission, we consider that there is a weak maser effect in the additional spectral component of 51–40 E methanol line at 84.521 GHz. We compared the molecular abundances with values found in several astrochemical templates: molecular clouds, hot cores, and photodissociation regions, and found that signatures of these different types can be found towards V645 Cyg. We also obtained maps of the cloud in several molecular emission lines. The peaks of CO and CH3OH emission are shifted from the direction of the star, but the CS, HCO+, HNC, HCN, and N2H+ emission peaks are observed directly towards the star. Exploring the gas kinematics around V645 Cyg, we found that velocity structure in the ambient molecular cloud on the scale ≈1.6–2.0 pc is not the same as within ≈0.5 pc found previously by other authors.

Version 8 IMK–IAA MIPAS temperatures from 12–15 µm spectra: Middle and Upper Atmosphere modes

Abstract. Motivated by an improved European Space Agency (ESA) version of calibrated Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) spectra (version 8.03), we have released version 8 of MIPAS temperatures and pointing information retrieved from 2005–2012 MIPAS measurements at 12–15 µm in the Middle Atmosphere (MA), Upper Atmosphere (UA) and Noctilucent Cloud (NLC) measurement modes. The Institute of Meteorology and Climate Research–Instituto de Astrofisica de Andalucia (IMK–IAA) retrieval processor in use considers non-local thermodynamic equilibrium (non-LTE) emission explicitly for each limb scan. This non-LTE treatment is essential to obtain accurate temperatures above the mid-mesosphere because at the altitudes covered, up to 115 km, the simplified climatology-based non-LTE treatment employed for the Nominal (NOM) measurements is insufficient. Other updates in MA/UA/NLC version 8 non-LTE temperature retrievals from previous data releases include more realistic atomic oxygen and carbon dioxide abundances, an updated set of spectroscopic data, an improved spectral shift retrieval, a continuum retrieval extended to altitudes up to 58 km, consideration of an altitude-dependent radiance offset retrieval, the use of wider microwindows above 85 km to capture the offset, an improved accuracy in forward model calculations, new a priori temperature information, improved temperature horizontal gradient retrievals and the use of MIPAS version 5 interfering species where available. The resulting MIPAS MA/UA/NLC IMK–IAA temperature dataset is reliable for scientific analysis in the full measurement vertical range for the MA (18–102 km) and the NLC (39–102 km) observations and from 42 to 115 km for the UA observations. The random temperature errors, dominated by the instrumental noise, are typically less than 1 K below 60 km, 1–3 K at 60–70 km, 3–5 K at 70–90 km, 6–8 K at 90–100 km, 8–12 K at 100–105 km and 12–20 K at 105–115 km. Random pointing correction errors, also mainly arising from instrumental noise, are on average 50 m for tangent altitudes up to 60 km and decrease linearly to values smaller than 20 m for altitudes above 95 km. The vertical resolution is 3 km at altitudes below 50 km, 3–5 km at 50–70 km, 4–6 km at 70–90 km, 6–10 km at 90–100 km and 8–11 km at 100–115 km. The systematic errors in retrieved temperatures below 75 km are driven by uncertainties in the CO2 spectroscopic data and, above 80 km, by uncertainties in the non-LTE model parameters (including collisional rates and atomic oxygen abundance) and the CO2 abundance. These lead to systematic temperature errors of less than 0.7 K below 55 km, 1 K at 60–80 km, 1–2 K at 80–90 km, 3 K at 95 km, 6–8 K at 100 km, 10–20 K at 105 km and 20–30 K at 115 km. Systematic errors in the tangent altitude correction, mainly arising from CO2 spectroscopic uncertainties, are 250 m at 20 km, 200 m at 40–60 km, 100 m at 80 km and smaller than 50 m above 90 km. The consistency between the MA/UA/NLC and the NOM IMK–IAA datasets is excellent below 70 km (typical 0.5–1 K differences). The comparison of this temperature dataset with co-located Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) temperature measurements shows excellent agreement, with differences typically within 1.5 K below 90 km, 1–3 K at 90–95 km, 1–5 K at 95–100 km, 1–8 K at 100–105 km and 1–10 K above. The agreement with SABER improves with respect to previous MIPAS IMK–IAA data versions.

Correction to: An accurate set of H3O+ − H2 collisional rate coefficients for non-LTE modelling of warm interstellar clouds

Correction to: An accurate set of H3O+ − H2 collisional rate coefficients for non-LTE modelling of warm interstellar clouds

TITANS metal-poor reference stars

Context. Representative samples of F-, G-, K-type stars located outside of the solar neighbourhood have started to become available in spectroscopic surveys. The fraction of metal-poor ([Fe/H] ≲ −0.8 dex) giants becomes increasingly relevant towards greater distances. In metal-poor stars, effective temperatures (Teff) based on local thermodynamic equilibrium (LTE) spectroscopy and on former colour–Teff relations – which are still widely used – have been reported to be inaccurate. It is necessary to recalibrate chemical abundances based on these Teff scales in the multiple available surveys in order to bring them to the same standard scale for their simultaneous use. To this end, a complete sample of standards is required, which until now has been restricted to a few stars with quasi-direct Teff measurements. Aims. We aim to provide a legacy sample of metal-poor standards with proven accurate atmospheric parameters. We add 47 giants to the TITANS metal-poor reference stars. Methods. We derived Teff using 3D non-LTE Hα modelling, the accuracy of which was tested against interferometry and with the Infra Red Flux Method (IRFM). We derived surface gravity (log ɡ) by fitting Mg I b triplet lines, and tested their accuracy against asteroseismology. Metallicity was derived using Fe II lines, and we find our results to be identical to the [Fe/H] derived from non-LTE spectral synthesis. Results. The Teff that we find using 3D non-LTE Hα is equivalent to interferometric and IRFM temperatures within a ±46 K uncertainty. We achieve precision of ~50 K for 34 stars with spectra with the highest signal-to-noise ratio (S/N). For log ɡ, we achieve a total uncertainty of ±0.15 dex. For [Fe/H], we obtain a total uncertainty of ±0.09 dex. We find that the ionisation equilibrium of Fe lines under LTE is not valid in metal-poor giants. LTE leads to a small but significant metallicity underestimation of ~0.1 dex when derived from weak Fe I lines, and only provided accurate Teff and log ɡ. This bias totally disappears under non-LTE.

Collisional excitation of PO+ by para-H2: potential energy surface, scattering calculations, and astrophysical applications

ABSTRACT We report the derivation of rate coefficients for the rotational (de-)excitation of PO+ induced by collisions with H2. The calculations were performed on a 4D potential energy surface, obtained on top of highly accurate ab initio energy points. Preliminary tests pointed out the low influence of the coupling between j = 0 and the higher rotational levels of H2 on the cross-sections values, thus allowing to neglect the rotational structure of H2. On this basis, state-to-state collisional rate coefficients were derived for temperatures ranging from 5 to 200 K. Radiative transfer calculations have been used to model the recent observation of PO+ in the G+0.693–0.027 molecular cloud, in order to evaluate the possible impact of non-LTE models on the determination of its physical conditions. The derived column density was found to be approximately ∼ 3.7 × 1011 cm−2, which is 60% (a factor of ∼ 1.7) smaller than the previously LTE-derived value. Extensive simulations show that PO+ low-j rotational lines exhibit maser behaviour at densities between 104 and 106 cm−3, thus highlighting the importance of a proper treatment of the molecular collisions to accurately model PO+ emissions in the interstellar medium.

Modeling of the High-velocity Jet Powered by the Massive Star MWC 349A

MWC 349A is a massive star with a well-known circumstellar disk rotating following a Keplerian law, and an ionized wind launched from the disk surface. Recent observations with the Atacama Large Millimeter/submillimeter Array (ALMA) carried out toward this system, however, have revealed an additional high-velocity component in the strong, maser emission of hydrogen radio recombination lines (RRLs), suggesting the presence of a high-velocity ionized jet. In this work, we present 3D non-LTE radiative transfer modeling of the emission of the H30α and H26α maser lines, and of their associated radio continuum emission, toward MWC 349A. By using the MORELI code, we reproduce the spatial distribution and kinematics of the high-velocity emission of the H30α and H26α maser lines with a high-velocity ionized jet expanding at a velocity of ∼250 km s−1, surrounded by MWC 349A’s wide-angle ionized wind. The bipolar jet, which is launched from MWC 349A’s disk, is poorly collimated and slightly misaligned with respect to the disk rotation axis. Thanks to the unprecedented sensitivity and spatial accuracy provided by ALMA, we also find that the already known, wide-angle ionized wind decelerates as it expands radially from the ionized disk. We briefly discuss the implications of our findings for understanding the formation and evolution of massive stars. Our results show the huge potential of RRL masers as powerful probes of the innermost ionized regions around massive stars and of their high-velocity jets.

A high HDO/H2O ratio in the Class I protostar L1551 IRS5

Context. Water is a very abundant molecule in star-forming regions. Its deuterium fractionation provides an important tool for understanding its formation and evolution during the star and planet formation processes. While the HDO/H2O abundance ratio has been determined toward several young Class 0 protostars and comets, the number of studies toward Class I protostars is limited. Aims. Our aim is to study the water deuteration toward the Class I protostellar binary L1551 IRS5 and to investigate the effect of evolutionary stage and environment on variations in the water D/H ratio. Methods. Observations were carried out toward L1551 IRS5 using the NOrthern Extended Millimeter Array (NOEMA) interferometer. The HDO 31, 2–22, 1 transition at 225.9 GHz and the H218O 31, 3–22, 0 transition at 203.4 GHz were covered with a spatial resolution of 0.5″× 0.8″, while the HDO 42, 2–42, 3 transition at 143.7 GHz was observed with a resolution of 2.0″ × 2.5″. We constrained the water D/H ratio using both local thermodynamic equilibrium (LTE) and non-LTE models. Results. The three transitions are detected. The line profiles display two peaks, one at ∼6 km s−1 and one at ∼9 km s−1. We derive an HDO/H2O ratio of (2.1 ± 0.8) × 10−3 for the redshifted component and a lower limit of > 0.3 × 10−3 for the blueshifted component. This lower limit is due to the blending of the blueshifted H218O component with redshifted CH3OCH3 emission. Conclusions. The HDO/H2O ratio in L1551 IRS5 is similar to the values in Class 0 isolated sources and in the disk of the Class I protostar V883 Ori, while it is significantly higher than in the previously studied clustered Class 0 sources and the comets. This result suggests that the chemistry of protostars that belong to molecular clouds with relatively low source densities, such as L1551, share more similarities with the isolated sources than the protostars of very dense clusters. If the HDO/H2O ratios in Class 0 protostars with few sources around are comparable to those found to date in isolated Class 0 objects, it would mean that there is little water reprocessing from the Class 0 to Class I protostellar stage.