Search

Context. HD140283, or the Methuselah star, is a well-known reference object in stellar evolution. Its peculiar chemical composition, proximity and absence of reddening makes it an interesting case-study of Pop II stars. Thanks to recent observational efforts, we now have precise interferometric and spectroscopic constraints, as well as revised astrometric parallaxes from the Gaia mission. Aims. We aim at determining the age of HD140283 with these lastest constraints, as well as quantifying the impact of systematics from physical inaccuracies in the stellar evolution models. Methods. Using recent spectroscopic abundances from the literature, including 3D non-LTE values for C, O, and Fe, we compute opacity tables specific to HD140283. We then use them in grids of stellar evolution models coupled to a Markov Chain Monte Carlo tool to determine the age of HD140283. Results. With our tailored models we find an age of 12.3Gy. Using a solar-scaled mixture instead results in an age value of 14Gy, in tension with the age of the universe ($13.77\pm0.06$Gy). We also find that reducing the mixing length parameter from its solar calibrated value will lead to an even lower age, in agreement with other recent studies. However, we find no direct evidence to favour a lower mixing length parameter value from our modelling. Conclusions. Taking into account the specific elemental abundances is crucial for the modelling of HD140283, as it leads to significant differences in the inferred age. However, this effect is degenerate with a lowering of the mixing length parameter. In this respect, asteroseismic constraints might play a key role in accurately deriving the mass of HD140283, therefore strongly constraining its age., Comment: Accepted for publication in A&A

The stars of the Milky Way carry the chemical history of our Galaxy in their atmospheres as they journey through its vast expanse. Like barcodes, we can extract the chemical fingerprints of stars from high-resolution spectroscopy. The fourth data release (DR4) of the Galactic Archaeology with HERMES (GALAH) Survey, based on a decade of observations, provides the chemical abundances of up to 32 elements for 917 588 stars that also have exquisite astrometric data from the $Gaia$ satellite. For the first time, these elements include life-essential nitrogen to complement carbon, and oxygen as well as more measurements of rare-earth elements critical to modern-life electronics, offering unparalleled insights into the chemical composition of the Milky Way. For this release, we use neural networks to simultaneously fit stellar parameters and abundances across the full spectrum, leveraging synthetic grids computed with Spectroscopy Made Easy. These grids account for atomic line formation in non-local thermodynamic equilibrium for 14 elements. In a two-iteration process, we first fit stellar labels for all 1 085 520 spectra, then co-add repeated observations and refine these labels using astrometric data from $Gaia$ and 2MASS photometry, improving the accuracy and precision of stellar parameters and abundances. Our validation thoroughly assesses the reliability of spectroscopic measurements and highlights key caveats for catalogue users. GALAH DR4 represents yet another milestone in Galactic archaeology, combining detailed chemical compositions from multiple nucleosynthetic channels with kinematic information and age estimates. The resulting dataset, covering nearly a million stars, opens new avenues for understanding not only the chemical and dynamical history of the Milky Way, but also the broader questions of the origin of elements and the evolution of planets, stars, and galaxies., Comment: 43 pages, 38 figures to be submitted to PASA. Accompanying the GALAH Data Release 4, see https://www.galah-survey.org and https://cloud.datacentral.org.au/teamdata/GALAH/public/GALAH_DR4/. All code available on http://github.com/svenbuder/GALAH_DR4/ and https://github.com/svenbuder/galah_dr4_paper. Comments welcome

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 II 313.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 $\pm$ 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 $\pm$ 0.06 dex, or 22 $\pm$ 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., Comment: 14 pages, 9 figures; accepted for publication in Astronomy & Astrophysics

The negative hydrogen ion H$^-$ is, almost without exception, treated in local thermodynamic equilibrium (LTE) in the modelling of F, G, and K stars, where it is the dominant opacity source in the visual spectral region. This assumption rests in practice on a study from the 1960s. Since that work, knowledge of relevant atomic processes and theoretical calculations of stellar atmospheres and their spectra have advanced significantly, but this question has not been reexamined. We present calculations based on a slightly modified analytical model that includes H, H$_2$, and H$^-$, together with modern atomic data and a grid of 1D LTE theoretical stellar atmosphere models with stellar parameters ranging from T$_\mathrm{eff} = 4000$ to 7000~K, $\log{g} = 1$ to 5 cm/s$^2$, and [Fe/H]$=-3$ to 0. We find direct non-LTE effects on populations in spectrum-forming regions, continua, and spectral lines of about 1-2% in stars with higher T$_\mathrm{eff}$ and/or lower $\log g$. Effects in models for solar parameters are smaller by a factor of 10, about 0.1-0.2%, and are practically absent in models with lower T$_\mathrm{eff}$ and/or higher $\log g$. These departures from LTE found in our calculations originate from the radiative recombination of electrons with hydrogen to form H$^-$ exceeding photodetachment, that is, overrecombination. Modern atomic data are not a source of significant differences compared to the previous work, although detailed data for processes on H$_2$ resolved with vibrational and rotational states provide a more complete and complex picture of the role of H$_2$ in the equilibrium of H$^-$. In the context of modern studies of stellar spectra at the percent level, our results suggest that this question requires further attention, including a more extensive reaction network, and indirect effects due to non-LTE electron populations., Comment: Accepted by Astronomy and Astrophysics. Replaced with new version following language editing, as well as correct font, and correction of some minor typos

Building on the success of previous missions, asteroseismic modelling will play a key role in future space-based missions, such as PLATO, CubeSpec, and Roman. Despite remarkable achievements, asteroseismology has revealed significant discrepancies in the physics of theoretical stellar models, which have the potential to bias stellar characterisation at the precision level demanded by PLATO. The current modelling strategies largely overlook magnetic activity, assuming that its effects are masked by filtering the so-called surface effects. Given the presence of activity cycles in multiple solar-like oscillators, and activity variations in a significant fraction of Kepler observations of main-sequence stars (Santos et al. 2019b, 2021, 2023), we measured the impact of magnetic activity on the asteroseismic characterisation of the Sun based on 26.5 years of GOLF and BiSON observations. While magnetic activity is partially absorbed in the treatment of surface effects, we found a discernible imprint of the activity cycle in the determination of the solar age. Notably, this imprint persists across both BiSON and GOLF datasets, with significant variations of up to 6.5% observed between solar minima and maxima. Considering that the Sun exhibits low levels of activity, our study underscores the looming challenge posed by magnetic activity for future photometry missions, and prompts a potential reevaluation of the asteroseismic characterisation of Kepler's most active targets., Comment: Accepted for publication as a Letter in Astronomy and Astrophysics

Magnesium is one of the important elements in stellar physics as an electron donor and in Galactic Archaeology as a discriminator of different stellar populations. However, previous studies of Mg I and Mg II lines in metal-poor benchmark stars have flagged problems with magnesium abundances inferred from one-dimensional (1D), hydrostatic models of stellar atmospheres, both with or without the local thermodynamic equilibrium (LTE) approximation. We here present 3D non-LTE calculations for magnesium in FG-type dwarfs, and provide corrections for 1D LTE abundances. The 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 then applied our abundance corrections to 1D LTE literature results for stars in the thin disc, thick disc, $\alpha$-rich halo, and $\alpha$-poor halo. We find that the 3D non-LTE results show a richer substructure in [Mg/Fe]-[Fe/H] in the $\alpha$-poor halo, revealing two subpopulations at the metal-rich end. These two subpopulations are also separated in kinematics, supporting the astrophysical origin of the separation. While the more magnesium-poor subpopulation is likely to be debris from a massive accreted galaxy, Gaia-Enceladus, the other subpopulation may be related to a previous identified group of stars, called Eos. The presence of additional separation in [Mg/Fe] suggests that previous Mg abundance measurements may have been limited in the precision by the 1D and LTE approximations, highlighting the importance of 3D non-LTE modelling., Comment: accepted to A\&A, 8 pages

The A-type metallic-line (Am) stars are typically considered to be non-magnetic or possessing very weak sub-G magnetic fields. This view has been repeatedly challenged in the literature, most commonly for the bright hot Am star o Peg. Several studies claimed to detect 1-2 kG field of unknown topology in this object, possibly indicating a new process of magnetic field generation in intermediate-mass stars. In this study, we revisit the evidence of a strong magnetic field in o Peg using new high-resolution spectropolarimetric observations and advanced spectral fitting techniques. The mean magnetic field strength in o Peg is estimated from the high-precision CRIRES+ measurement of near-infrared sulphur lines. This observation is modelled with a polarised radiative transfer code, including treatment of the departures from local thermodynamic equilibrium. In addition, the least-squares deconvolution multi-line technique is employed to derive longitudinal field measurements from archival optical spectropolarimetric observations of this star. Our analysis of the near-infrared S I lines reveals no evidence of Zeeman broadening, ruling out magnetic field with a strength exceeding 260 G. This null result is compatible with the relative intensification of Fe II lines in the optical spectrum taking into account blending and uncertain atomic parameters of the relevant diagnostic transitions. Longitudinal field measurements at three different nights also yield null results with a precision of 2 G. This study refutes the claims of kG-strength dipolar or tangled magnetic field in o Peg. This star is effectively non-magnetic, with the surface magnetic field characteristics no different from those of other Am stars., Comment: 11 pages, 8 figures; accepted for publication in A&A

Context: Traditional one-dimensional (1D) hydrostatic model atmospheres introduce systematic modelling errors into spectroscopic analyses of FGK-type stars. Aims: We present an updated version of the STAGGER-grid of 3D model atmospheres, and explore the accuracy of post-processing methods in preparation for spectral synthesis. Methods: New and old models were (re)computed following an updated workflow, including an updated opacity binning technique. Spectroscopic tests were performed in 3D LTE for a grid of 216 fictitious Fe I lines, spanning a wide range in oscillator strength, excitation potential and central wavelength, and eight model atmospheres that cover the stellar atmospheric parameter range (Teff, log g, [Fe/H]) of FGK-type stars. Using this grid, the impact of vertical and horizontal resolution, and temporal sampling of model atmospheres on spectroscopic diagnostics was tested. Results: We find that downsampling the horizontal mesh from its original size of 240 x 240 grid cells to 80 x 80 cells, i.e. sampling every third grid cell, introduces minimal errors on the equivalent width and normalized line flux across the line and stellar parameter space. Regarding temporal sampling, we find that sampling ten statistically independent snapshots is sufficient to accurately model the shape of spectral line profiles. For equivalent widths, a subsample consisting of only two snapshots is sufficient, introducing an abundance error of less than 0.015 dex. Conclusions: We have computed 32 new model atmospheres and recomputed 116 old model atmospheres present in the original grid. The public release of the STAGGER-grid contains 243 models, excluding models with [Fe/H] = -4.00, and the processed snapshots can be used to improve the accuracy of spectroscopic analyses.

The first (Pop III) stars formed only out of H and He and were likely more massive than present-day stars. Massive Pop III stars in the range 140-260 M$_\odot$ 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$_{\odot}$) was discovered by the LAMOST survey, the star J1010+2358. However, the 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 VLT/UVES spectrum, correcting for 3D and/or non-LTE effects. Our measurements of both C and Al give much higher values (~1 dex) than expected from a 260 M$_{\odot}$ PISN. Furthermore, we find significant discrepancies with the previous analysis, and therefore a much less pronounced odd-even pattern. Thus, we show that J1010+2358 cannot be a pure descendant of a 260 M$_{\odot}$ PISN. Instead, we find that the best fit model consists of a 13 M$_{\odot}$ Pop II core-collapse supernova combined with a Pop III supernova. Alternative, less favoured solutions $(\chi^2/\chi^2_{\rm best}\approx2.3)$ include a 50% contribution from a 260 M$_{\odot}$ PISN, or a 40% contribution from a Pop 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., Comment: Accepted to ApJ Letters (3rd of May 2024). 1 table, 4 figures, and Appendix

Deep imaging of galaxy cluster fields have in recent years revealed tens of candidates for gravitationally lensed stars at redshifts $z\approx$ 1-6, and future searches are expected to reveal highly magnified stars from even earlier epochs. Multi-band photometric observations may be used to constrain the redshift, effective temperature $T_\mathrm{eff}$ and dust attenuation along the line of sight to such objects. When combined with an estimate of the likely magnification, these quantities may be converted into a constraint on the stellar luminosity and, for an adopted set of stellar evolutionary tracks, the initial stellar mass. Further characterization is, however, difficult without spectroscopic observations, which at the typical brightness levels of high-redshift lensed stars becomes extremely challenging for even the largest existing telescopes. Here, we explore what spectral features one can realistically hope to detect in lensed stars with peak brightness in the range 26-28 AB mag, $T_\mathrm{eff}=$ 4000-50 000 K and redshifts $z=$1-10, using spectroscopy with the James Webb Space Telescope (JWST) and the forthcoming Extremely Large Telescope. We find that a majority of detectable lines appear in the rest UV-range for stars with $T_\mathrm{eff}\geq$15 000 K. The strongest detectable spectral lines are the C IV $\lambda$ 1550 \r{A} line and the Si IV $\lambda\lambda$1393, 1403 \r{A}-doublet at $T_\mathrm{eff}=$30 000 K. For lower temperatures, the calcium H- and K-lines at $T_\mathrm{eff}=$6000 K are among the most readily detectable. In limited wavelength ranges, ELT is expected to provide more sensitive spectroscopic observations, and with higher resolution than JWST. We find that variations of both mass loss rate and metallicity lead to noticeable effects in the detectability of certain spectral lines with both JWST and ELT., Comment: 12 pages, 6 figures, accepted in A&A

Due to the high quality constraints available for the Sun, we can carry out combined analyses using neutrino, spectroscopic and helioseismic observations. Such studies lay the ground for future improvements of key physical components of solar and stellar models, such as the equation of state, radiative opacities or prescriptions for macroscopic transport processes of chemicals which are then used to study other stars in the Universe. We study the existing degeneracies in solar models using the recent high-metallicity spectroscopic abundances by comparing them to helioseismic and neutrino data and discuss how their properties are impacted by changes in various physical ingredients. We carry out a detailed study of solar models computed with a high-metallicity composition from the literature based on averaged-3D models that was claimed to solve the solar problem. The properties of the solar models are significantly affected by using the recent OPLIB opacities and the inclusion of macroscopic transport. The properties of the standard solar models computed using the OPAL opacities are similar to those using the OP opacities. We show that a modifying the temperature gradient just below the base of the convective zone is required to erase the discrepancies in solar models, particularly in the presence of macroscopic mixing. This can be simulated by a local increase of opacity of a few percent. We conclude that the existing degeneracies and issues in solar modelling are not erased by an increase in the solar metallicity in contradiction to was suggested in recent papers. Therefore, standard solar models cannot be used as an argument for a high metallicity composition. While further work is required to improve solar models, we note that direct helioseismic inversions indicate a low metallicity in the convective envelope, in agreement with spectroscopic analyses based on full 3D models., Comment: Accepted for publication in Astronomy and Astrophysics, version before language editing

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. 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. A large grid of departure coefficients, $\beta_\nu$, were computed on standard MARCS model atmospheres. The grid extends from 3000K to 8000K in T$_{\mathrm{eff}}$, -0.5 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 and additional metal-poor dwarfs. The non-LTE effects grow towards lower [Fe/H], lower log(g), and higher T$_{\mathrm{eff}}$, with a minimum and maximum $\Delta$A(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 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. 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.

The lack of reliable atomic data can be a severe limitation in astrophysical modelling, in particular of events such as kilonovae that require information on all neutron-capture elements across a wide range of ionization stages. Notably, the presence of non-orthonormalities between electron orbitals representing configurations that are close in energy can introduce significant inaccuracies in computed energies and transition probabilities. Here, we propose an explicit targeted optimization method that can effectively circumvent this concern while retaining an orthonormal orbital basis set. We illustrate this method within the framework of small-scale atomic structure models of Au I, using the GRASP2018 multiconfigurational Dirac-Hartree-Fock atomic structure code. By comparing to conventional optimization schemes we show how a targeted optimization approach improves the energy level positioning and ordering. Targeted optimization also leads to better agreement with experimental data for the strongest E1 transitions. This illustrates how small-scale models can be significantly improved with minor computational costs if orbital non-orthonormalities are considered carefully. These results should prove useful to multi-element atomic structure calculations in, for example, astrophysical opacity applications involving neutron-capture elements., Comment: 14 pages, 3 figures, 6 tables; published in Journal of Physics B

Previous work on the abundances of C, O, Na, Mg, Si, Ca, Ti, Cr, Mn, Fe, Ni, Cu, and Zn in low-alpha (accreted) and high-alpha (in situ born) halo stars is extended to include the abundances of Sc, V, and Co, enabling us to study the nucleosynthesis of all iron-peak elements along with the lighter elements. The Sc, V, and Co abundances were determined from a 1D MARCS model-atmosphere analysis of equivalent widths of atomic lines in high signal-to-noise, high resolution spectra assuming local thermodynamic equilibrium (LTE). In addition, new 3D and/or non-LTE calculations were used to correct the 1D LTE abundances for several elements including consistent 3D non-LTE calculations for Mg. The two populations of accreted and in situ born stars are well separated in diagrams showing [Sc/Fe], [V/Fe], and [Co/Fe] as a function of [Fe/H]. The [X/Mg] versus [Mg/H] trends for high-alpha and low-alpha stars were used to determine the yields of core-collapse and Type Ia supernovae. The largest Type Ia contribution occurs for Cr, Mn, and Fe, whereas Cu is a pure core-collapse element. Sc, Ti, V, Co, Ni, and Zn represent intermediate cases. A comparison with yields calculated for supernova models shows poor agreement for the core-collapse yields. The Ia yields suggest that sub-Chandrasekhar-mass Type Ia supernovae provide a dominant contribution to the chemical evolution of the host galaxies of the low-alpha stars. A substructure in the abundances and kinematics of the low-alpha stars suggests that they arise from at least two different satellite accretion events, Gaia-Sausage-Enceladus and Thamnos., Comment: Accepted for publication in A&A. 12 pages, 4 tables, and 7 figures

Context. Transmission spectroscopy is one of the most powerful techniques to characterize transiting exoplanets since it allows to measure the abundance of the atomic and molecular species in the planetary atmosphere. However, the stellar lines can 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 non-local thermodynamic equilibrium (non-LTE) line formation are required for an accurate modeling of the stellar CLV of the Na I D$_1$ and K I resonance lines on transmission spectra. Methods. We model the CLV of the Na I D$_1$ and K I resonance lines in the Sun with 3D non-LTE radiative transfer. The synthetic spectra are 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 $\mu=0.1$ and $\mu=1.0$. Results. Our 3D non-LTE modeling of the Na I D$_1$ resonance line at 5896 {\AA} and the K I 7699 {\AA} resonance line in the Sun is in good agreement with the observed CLV in the solar spectrum. The simulated CLV curve for a Jupiter-Sun system inferred with a 3D non-LTE analysis shows significant differences from that obtained from a 1D atmosphere. The latter tends to overestimate the amplitude of the transmission curve by a factor that is of the same order of magnitude as a planetary absorption depth (up to 0.2 %). Conclusions. In order to correctly characterize exoplanetary atmospheres, 3D non-LTE synthetic spectra should be used to estimate the stellar CLV effect in transmission spectra of solar-like planet hosts. The work will be extended to other lines and FGK-type stars, allowing synthetic high-resolution spectra to mitigate the stellar contamination of low-resolution planetary spectra, e.g. those from JWST., Comment: 23 pages, accepted for publication by A&A

Context. The metal mass fraction of the Sun Z is a key constraint in solar modelling, but its value is still under debate. The standard solar chemical composition of the late 2000s have the ratio of metals to hydrogen Z/X = 0.0181, with a small increase to 0.0187 in 2021, as inferred from 3D non-LTE spectroscopy. However, more recent work on a horizontally and temporally averaged <3D> model claim Z/X = 0.0225, consistent with the high values of twenty-five years ago based on 1D LTE spectroscopy. Aims. We aim to determine a precise and robust value of the solar metal mass fraction from helioseismic inversions, thus providing independent constraints from spectroscopic methods. Methods. We devise a detailed seismic reconstruction technique of the solar envelope, combining multiple inversions and equations of state to accurately and precisely determine the metal mass fraction value. Results. We show that a low value of the solar metal mass fraction corresponding to Z/X = 0.0187 is favoured by helioseismic constraints and that a higher metal mass fraction corresponding to Z/X = 0.0225 are strongly rejected by helioseismic data. Conclusions. We conclude that direct measurement of the metal mass fraction in the solar envelope favours a low metallicity, in line with the 3D non-LTE spectroscopic determination of 2021. A high metal mass fraction as measured using a <3D> model in 2022 is disfavoured by helioseismology for all modern equations of state used to model the solar convective envelope., Comment: Accepted for publication in A&A, online materials available on demand or on publisher website

Three-dimensional radiation-hydrodynamics (3D RHD) simulations of stellar surface convection provide valuable insights into many problems in solar and stellar physics. However, almost all 3D near-surface convection simulations to date are based on solar-scaled chemical compositions, which limit their application on stars with peculiar abundance patterns. To overcome this difficulty, we implement the robust and widely-used FreeEOS equation of state and our Blue opacity package into the Stagger 3D radiation-magnetohydrodynamics code. We present a new 3D RHD model of the solar atmosphere, and demonstrate that the mean stratification as well as the distributions of key physical quantities are in good agreement with those of the latest Stagger solar model atmosphere. The new model is further validated by comparing against solar observations. The new model atmospheres reproduce the observed flux spectrum, continuum centre-to-limb variation, and hydrogen line profiles at a satisfactory level, thereby confirming the realism of the model and the underlying input physics. These implementations open the prospect for studying other stars with different $\alpha$-element abundance, carbon-enhanced metal-poor stars and population II stars with peculiar chemical compositions using 3D Stagger model atmospheres., Comment: 24 pages, 20 figures, accepted for publication in A&A

As the most abundant element in the universe after hydrogen and helium, oxygen plays a key role in planetary, stellar, and galactic astrophysics. Its abundance is especially influential on stellar structure and evolution, and as the dominant opacity contributor at the base of the Sun's convection zone it is central to the discussion around the solar modelling problem. However, abundance analyses require complete and reliable sets of atomic data. We present extensive atomic data for O I, by using the multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods. Lifetimes and transition probabilities for radiative electric dipole transitions are given and compared with results from previous calculations and available measurements. The accuracy of the computed transition rates is evaluated by the differences between the transition rates in Babushkin and Coulomb gauges, as well as by a cancellation factor analysis. Out of the 989 computed transitions in this work, 205 are assigned to the accuracy classes AA-B, that is, with uncertainties less than 10%, following the criteria defined by the National Institute of Standards and Technology Atomic Spectra Database. We discuss the influence of the new log(gf) values on the solar oxygen abundance and ultimately advocate $\log\epsilon_{\mathrm{O}}=8.70\pm0.04$., Comment: 13 pages, 5 figures; Accepted for publication in Astronomy & Astrophysics

The origin of magnetic fields and their role in chemical spot formation on magnetic Ap stars is currently not understood. Here we contribute to solving this problem with a detailed observational characterisation of the surface structure of 45 Her, a weak-field Ap star. We find this object to be a long-period, single-lined spectroscopic binary and determine the binary orbit as well as fundamental and atmospheric parameters of the primary. We study magnetic field topology and chemical spot distribution of 45 Her with the help of the Zeeman Doppler imaging technique. Magnetic mapping reveals the stellar surface field to have a distorted dipolar topology with a surface-averaged field strength of 77 G and a dipolar component strength of 119 G - confirming it as one of the weakest well-characterised Ap-star fields known. Despite its feeble magnetic field, 45 Her shows surface chemical inhomogeneities with abundance contrasts of up to 6 dex. Of the four chemical elements studied, O concentrates at the magnetic equator whereas Ti, Cr and Fe avoid this region. Apart from this trend, the positions of Fe-peak element spots show no apparent correlation with the magnetic field geometry. No signs of surface differential rotation or temporal evolution of chemical spots on the time scale of several years were detected. Our findings demonstrate that chemical spot formation does not require strong magnetic fields to proceed and that both the stellar structure and the global field itself remain stable for sub-100 G field strengths contrary to theoretical predictions., Comment: 15 pages, 19 figures; Accepted for publication in MNRAS

Context: The first stars produced the first heavy elements and set the stage for the formation of the first galaxies. Accurate chemical abundances of ultra metal-poor stars ([Fe/H]<-4) can be used to infer properties of the first stars, and thus the formation mechanism for low-mass second generation stars in the early universe. Spectroscopic studies have shown that most second generation stars are carbon-enhanced with one notable exception SDSS J102915.14+172927.9. Aims: We reanalyse the composition of SDSS J102915.14+172927.9. Methods: We developed a tailored 3D model atmosphere for SDSS J102915.14+172927.9 with the Stagger-code, making use of an improved surface gravity estimate based on the Gaia DR3 parallax. This model was used as input in the radiative transfer code Balder to compute 3D non-LTE synthetic spectra. These spectra were then used to infer abundances for Mg, Si, Ca, Fe and Ni, and upper limits on Li, Na and Al. 3D LTE synthetic spectra were computed with Scate to infer the abundance of Ti and upper limits on C and N. Results: In contrast to earlier works based on 1D non-LTE corrections to 3D LTE results, we are able to achieve ionisation balance for Ca I and Ca II when employing our consistent 3D non-LTE treatment. Moreover, the elemental abundances are systematically higher than those found in earlier works. In particular, [Fe/H] increases by 0.57 dex, and the upper limits of C and N increase by 0.90 dex and 1.82 dex, respectively. Conclusions: We find that Population III progenitors with masses 10-20 M_sun exploding with energy E<=3*10^{51} erg can reproduce our 3D non-LTE abundance pattern. Contrary to previous work, we obtain higher upper limits on the carbon abundance that are ``marginally consistent'' with star formation through atomic line cooling, and as such, prevent strong conclusions about the formation mechanism of this low mass star., Comment: 14 pages, 9 figures, Accepted for publication in Astronomy & Astrophysics

Accurate and extensive atomic data are essential for spectroscopic analyses of stellar atmospheres and other astronomical objects. We present energy levels, lifetimes, and transition probabilities for neutral nitrogen, the sixth most abundant element in the cosmos. The calculations employ the fully relativistic multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods, and span the 103 lowest states up to and including 2s$^2$2p$^2$5s. Our theoretical energies are in excellent agreement with the experimental data, with an average relative difference of 0.07%. In addition, our transition probabilities are in good agreement with available experimental and theoretical data. We further verify the agreement of our data with experimental results via a re-analysis of the solar nitrogen abundance, with the results from the Babushkin and Coulomb gauges consistent to 2% or 0.01 dex. We estimated the uncertainties of the computed transition data based on a statistical analysis of the differences between the transition rates in Babushkin and Coulomb gauges. Out of the 1701 computed electric dipole transitions in this work, 83 (536) are associated with uncertainties less than 5% (10%)., Comment: 17 pages, 7 figures; Accepted for publication in The Astrophysical Journal Supplement Series

The very massive first stars ($m>100\rm M_{\odot}$) were fundamental to the early phases of reionization, metal enrichment, and super-massive black hole formation. Among them, those with $140\leq\rm m/\rm M_{\odot}\leq260$ are predicted to evolve as Pair Instability Supernovae (PISN) leaving a unique chemical signature in their chemical yields. Still, despite long searches, the stellar descendants of PISN remain elusive. Here we propose a new methodology, the PISN-explorer, to identify candidates for stars with a dominant PISN enrichment. The PISN-explorer is based on a combination of physically driven models, and the FERRE code; and applied to data from large spectroscopic surveys (APOGEE, GALAH, GES, MINCE, and the JINA database). We looked into more than 1.4 million objects and built a catalogue with 166 candidates of PISN descendants. One of which, 2M13593064+3241036, was observed with UVES at VLT and full chemical signature was derived, including the killing elements, Cu and Zn. We find that our proposed methodology is efficient in selecting PISN candidates from both the Milky Way and dwarf satellite galaxies such as Sextans or Draco. Further high-resolution observations are highly required to confirm our best selected candidates, therefore allowing us to probe the existence and properties of the very massive First Stars., Comment: 13 pages, 7 figures Accepted for publication in MNRAS

Context. The Sun acts as a cornerstone of stellar physics. Thanks to spectroscopic, helioseismic and neutrino flux observations, we can use the Sun as a laboratory of fundamental physics in extreme conditions. The conclusions we draw are then used to inform and calibrate evolutionary models of all other stars in the Universe. However, solar models are in tension with helioseismic constraints. The debate on the ``solar problem'' has hitherto led to numerous publications discussing potential issues with solar models and abundances. Aims. Using the recently suggested high-metallicity abundances for the Sun, we investigate whether standard solar models, as well as models with macroscopic transport reproducing the solar surface lithium abundances and analyze their properties in terms of helioseismic and neutrino flux observations. Methods. We compute solar evolutionary models and combine spectroscopic and helioseismic constraints as well as neutrino fluxes to investigate the impact of macroscopic transport on these measurements. Results. When high-metallicity solar models are calibrated to reproduce the measured solar lithium depletion, tensions arise with respect to helioseismology and neutrino fluxes. This is yet another demonstration that the solar problem is also linked to the physical prescriptions of solar evolutionary models and not to chemical composition alone. Conclusions. A revision of the physical ingredients of solar models is needed in order to improve our understanding of stellar structure and evolution. The solar problem is not limited to the photospheric abundances if the depletion of light elements is considered. In addition, tighter constraints on the solar beryllium abundance will play a key role in the improvement of solar models., Comment: Accepted for publication in Astronomy and Astrophysics

Pseudo two-colour diagrams or Chromosome maps (ChM) indicate that NGC 2808 host five different stellar populations. The existing ChMs have been derived by the Hubble Space Telescope photometry, and comprise of stars in a small field of view around the cluster centre. To overcome these limitations, we built a ChM with U,B,I photometry from ground-based facilities that disentangle the multiple stellar populations of NGC 2808 over a wider field of view. We used spectra collected by GIRAFFE@VLT in a sample of 70 red giant branch (RGB) and seven asymptotic giant branch (AGB) stars to infer the abundances of C, N, O, Al, Fe, and Ni, which combined with literature data for other elements (Li, Na, Mg, Si, Ca, Sc, Ti, Cr and Mn), and together with both the classical and the new ground-based ChMs, provide the most complete chemical characterisation of the stellar populations in NGC 2808 available to date. As typical of the multiple population phenomenon in globular clusters, the light elements vary from one stellar population to another; whereas the iron peak elements show negligible variation between the different populations (at a level of $\lesssim0.10$~dex). Our AGB stars are also characterised by the chemical variations associated with the presence of multiple populations, confirming that this phase of stellar evolution is affected by the phenomenon as well. Intriguingly, we detected one extreme O-poor AGB star (consistent with a high He abundance), challenging stellar evolution models which suggest that highly He-enriched stars should avoid the AGB phase and evolve as AGB-manqu\'e star., Comment: Accepted for publication in MNRAS. 19 pages and 13 figures

Star formation models predict that the metal-poor initial mass function (IMF) can be substantially different from that observed in the metal-rich Milky Way. This changeover occurs because metal-poor gas clouds cool inefficiently due to their lower abundance of metals and dust. However, predictions for the metal-poor IMF to date rely on assuming Solar-scaled abundances, that is, [X/O] = 0 at all [O/H]. There is now growing evidence that elements such as C and O that dominate metal line cooling in the ISM do not follow Solar scaling at low metallicities. In this work, we extend models that predict the variation in the characteristic (or, the peak) IMF mass as a function of metallicity using [C/O] ratios derived from observations of metal-poor Galactic stars and of H II regions in dwarf galaxies. These data show [C/O] < 0 at sub-Solar [O/H], which leads to a substantially different metal-poor IMF in the metallicity range where C I and C II cooling dominate ISM thermodynamics, resulting in an increase in the characteristic mass by a factor as large as 7. An important consequence of this difference is a shift in the location of the transition from a top- to a bottom-heavy IMF upwards by 0.5 $-$ 1 dex in metallicity. Our findings indicate that the IMF is very sensitive to the assumptions around Solar-scaled ISM compositions in metal-poor systems (e.g., dwarf galaxies, the Galactic halo and metal-poor stars) that are a key focus of JWST., Comment: 14 pages, 8 figures. Accepted by MNRAS

The titanium abundances of late-type stars are important tracers of Galactic formation history. However, abundances inferred from Ti I and Ti II lines can be in stark disagreement in very metal-poor giants. Departures from local thermodynamic equilibrium (LTE) have a large impact on the minority neutral species and thus influences the ionisation imbalance, but satisfactory non-LTE modelling for both dwarfs and giants has not been achieved in previous literature. The reliability of titanium abundances is reassessed in benchmark dwarfs and giants using a new non-LTE model and one-dimensional (1D) model atmospheres. A comprehensive model atom was compiled with a more extended level structure and newly published data for inelastic collisions between Ti I and neutral hydrogen. In 1D LTE, the Ti I and Ti II lines agree to within $0.06$ dex for the Sun, Arcturus, and the very metal-poor stars HD84937 and HD140283. For the very metal-poor giant HD122563, the Ti I lines give an abundance that is $0.47$ dex lower than that from Ti II. The 1D non-LTE corrections can reach $+0.4$ dex for individual Ti I lines and $+0.1$ dex for individual Ti II lines, and reduce the overall ionisation imbalance to $-0.17$ dex for HD122563. However, it also increases the imbalance for the very metal-poor dwarf and sub-giant to around $0.2$ dex. Using 1D non-LTE reduces the ionisation imbalance in very metal-poor giants but breaks the balance of other very metal-poor stars, consistent with the conclusions in earlier literature. To make further progress, consistent 3D non-LTE models are needed., Comment: 9 pages plus appendix, 6 figures; accepted for publication in Astronomy & Astrophysics

Spectroscopic measurements of iron abundances are prone to systematic modelling errors. We present 3D non-LTE calculations across 32 STAGGER-grid models with effective temperatures from 5000 K to 6500 K, surface gravities of 4.0 dex and 4.5 dex, and metallicities from $-$3 dex to 0 dex, and study the effects on 171 Fe I and 12 Fe II optical lines. In warm metal-poor stars, the 3D non-LTE abundances are up to 0.5 dex larger than 1D LTE abundances inferred from Fe I lines of intermediate excitation potential. In contrast, the 3D non-LTE abundances can be 0.2 dex smaller in cool metal-poor stars when using Fe I lines of low excitation potential. The corresponding abundance differences between 3D non-LTE and 1D non-LTE are generally less severe but can still reach $\pm$0.2 dex. For Fe II lines the 3D abundances range from up to 0.15 dex larger, to 0.10 dex smaller, than 1D abundances, with negligible departures from 3D LTE except for the warmest stars at the lowest metallicities. The results were used to correct 1D LTE abundances of the Sun and Procyon (HD 61421), and of the metal-poor stars HD 84937 and HD 140283, using an interpolation routine based on neural networks. The 3D non-LTE models achieve an improved ionisation balance in all four stars. In the two metal-poor stars, they remove excitation imbalances that amount to 250 K to 300 K errors in effective temperature. For Procyon, the 3D non-LTE models suggest [Fe/H] = 0.11 $\pm$ 0.03, which is significantly larger than literature values based on simpler models. We make the 3D non-LTE interpolation routine for FG-type dwarfs publicly available, in addition to 1D non-LTE departure coefficients for standard MARCS models of FGKM-type dwarfs and giants. These tools, together with an extended 3D LTE grid for Fe II from 2019, can help improve the accuracy of stellar parameter and iron abundance determinations for late-type stars., Comment: 17 pages, 11 figures, 5 tables; arXiv abstract abridged; accepted for publication in Astronomy & Astrophysics

Carbon, nitrogen, and oxygen are the fourth, sixth, and third most abundant elements in the Sun. Their abundances remain hotly debated due to the so-called solar modelling problem that has persisted for almost $20$ years. We revisit this issue by presenting a homogeneous analysis of $408$ molecular lines across $12$ diagnostic groups, observed in the solar intensity spectrum. Using a realistic 3D radiative-hydrodynamic model solar photosphere and LTE (local thermodynamic equilibrium) line formation, we find $\log\epsilon_{C} = 8.47\pm0.02$, $\log\epsilon_{N} = 7.89\pm0.04$, and $\log\epsilon_{O} = 8.70\pm0.04$. The stipulated uncertainties mainly reflect the sensitivity of the results to the model atmosphere; this sensitivity is correlated between the different diagnostic groups, which all agree with the mean result to within $0.03$ dex. For carbon and oxygen, the molecular results are in excellent agreement with our 3D non-LTE analyses of atomic lines. For nitrogen, however, the molecular indicators give a $0.12$ dex larger abundance than the atomic indicators, and our best estimate of the solar nitrogen abundance is given by the mean: $7.83$ dex. The solar oxygen abundance advocated here is close to our earlier determination of $8.69$ dex, and so the present results do not significantly alleviate the solar modelling problem., Comment: 12 pages, 6 figures, 3 tables; accepted for publication in Astronomy & Astrophysics

Although true metal-free "Population III" stars have so-far escaped discovery, their nature, and that of their supernovae, is revealed in the chemical products left behind in the next generations of stars. Here we report the detection of an ultra-metal poor star in the Sculptor dwarf spheroidal galaxy, AS0039. With [Fe/H]$_{\rm LTE}=-4.11$, it is the most metal-poor star so far discovered in any external galaxy. Contrary to the majority of Milky Way stars at this metallicity, AS0039 is clearly not enhanced in carbon, with [C/Fe]$_{\rm LTE}=-0.75$ and A(C)=+3.60, making it the lowest detected carbon abundance in any star to date. It furthermore lacks $\alpha$-element uniformity, having extremely low [Mg/Ca]$_{\rm NLTE}=-0.60$ and [Mg/Ti]$_{\rm NLTE}=-0.86$, in stark contrast with the near solar ratios observed in C-normal stars within the Milky Way halo. The unique abundance pattern indicates that AS0039 formed out of material that was predominantly enriched by a $\sim$20$ M_\odot$ progenitor star with an unusually high explosion energy $E=10\times10^{51}$ erg. The star AS0039 is thus one of the first observational evidence for zero-metallicity hypernovae and provides a unique opportunity to investigate the diverse nature of Population III stars., Comment: Accepted in ApJL, 7 pages, 3 figures, 1 table, and an appendix (1 figure and 3 tables)

The chemical composition of the Sun is a fundamental yardstick in astronomy, relative to which essentially all cosmic objects are referenced. We reassess the solar abundances of all 83 long-lived elements, using highly realistic solar modelling and state-of-the-art spectroscopic analysis techniques coupled with the best available atomic data and observations. Our new improved analysis confirms the relatively low solar abundances of C, N, and O obtained in our previous 3D-based studies: $\log\epsilon_{\text{C}}=8.46\pm0.04$, $\log\epsilon_{\text{N}}=7.83\pm0.07$, and $\log\epsilon_{\text{O}}=8.69\pm0.04$. The revised solar abundances for the other elements also typically agree well with our previously recommended values with just Li, F, Ne, Mg, Cl, Kr, Rb, Rh, Ba, W, Ir, and Pb differing by more than $0.05$ dex. The here advocated present-day photospheric metal mass fraction is only slightly higher than our previous value, mainly due to the revised Ne abundance from Genesis solar wind measurements: $X_{\rm surface}=0.7438\pm0.0054$, $Y_{\rm surface}=0.2423\pm 0.0054$, $Z_{\rm surface}=0.0139\pm 0.0006$, and $Z_{\rm surface}/X_{\rm surface}=0.0187\pm 0.0009$. Overall the solar abundances agree well with those of CI chondritic meteorites but we identify a correlation with condensation temperature such that moderately volatile elements are enhanced by $\approx 0.04$ dex in the CI chondrites and refractory elements possibly depleted by $\approx 0.02$ dex, conflicting with conventional wisdom of the past half-century. Instead the solar chemical composition resembles more closely that of the fine-grained matrix of CM chondrites. The so-called solar modelling problem remains intact with our revised solar abundances, suggesting shortcomings with the computed opacities and/or treatment of mixing below the convection zone in existing standard solar models., Comment: 31 pages, arXiv abstract abridged; accepted for publication in Astronomy & Astrophysics

We establish a quantitative relationship between photometric and spectroscopic detections of solar-like oscillations using ab initio, three-dimensional (3D), hydrodynamical numerical simulations of stellar atmospheres. We present a theoretical derivation as proof of concept for our method. We perform realistic spectral line formation calculations to quantify the ratio between luminosity and radial velocity amplitude for two case studies: the Sun and the red giant $\epsilon$ Tau. Luminosity amplitudes are computed based on the bolometric flux predicted by 3D simulations with granulation background modelled the same way as asteroseismic observations. Radial velocity amplitudes are determined from the wavelength shift of synthesized spectral lines with methods closely resembling those used in BiSON and SONG observations. Consequently, the theoretical luminosity to radial velocity amplitude ratios are directly comparable with corresponding observations. For the Sun, we predict theoretical ratios of 21.0 and 23.7 ppm/[m/s] from BiSON and SONG respectively, in good agreement with observations 19.1 and 21.6 ppm/[m/s]. For $\epsilon$ Tau, we predict K2 and SONG ratios of 48.4 ppm/[m/s], again in good agreement with observations 42.2 ppm/[m/s], and much improved over the result from conventional empirical scaling relations which gives 23.2 ppm/[m/s]. This study thus opens the path towards a quantitative understanding of solar-like oscillations, via detailed modelling of 3D stellar atmospheres., Comment: 16 pages, 7 figures, accepted for publication in MNRAS

Accurate atomic data are essential for opacity calculations and for abundance analyses of the Sun and other stars. The aim of this work is to provide accurate and extensive results of energy levels and transition data for C I - IV. The Multiconfiguration Dirac-Hartree-Fock and relativistic configuration interaction methods were used in the present work. To improve the quality of the wave functions and reduce the relative differences between length and velocity forms for transition data involving high Rydberg states, alternative computational strategies were employed by imposing restrictions on the electron substitutions when constructing the orbital basis for each atom and ion. Transition data, e.g., weighted oscillator strengths and transition probabilities, are given for radiative electric dipole (E1) transitions involving levels up to 1s$^2$2s$^2$2p6s for C I, up to 1s$^2$2s$^2$7f for C II, up to 1s$^2$2s7f for C III, and up to 1s$^2$8g for CIV. Using the difference between the transition rates in length and velocity gauges as an internal validation, the average uncertainties of all presented E1 transitions are estimated to be 8.05%, 7.20%, 1.77%, and 0.28%, respectively, for C I - IV. Extensive comparisons with available experimental and theoretical results are performed and good agreement is observed for most of the transitions. In addition, the C I data were employed in a reanalysis of the solar carbon abundance. The new transition data give a line-by-line dispersion similar to the one obtained when using transition data that are typically used in stellar spectroscopic applications today., Comment: 21 pages, 4 figures

Advances in merged-beams instruments have allowed experimental studies of the mutual neutralisation (MN) processes in collisions of both Li$^+$ and Na$^+$ ions with D$^-$ at energies below 1 eV. These experimental results place constraints on theoretical predictions of MN processes of Li$^+$ and Na$^+$ with H$^-$, important for non-LTE modelling of Li and Na spectra in late-type stars. We compare experimental results with calculations for methods typically used to calculate MN processes, namely the full quantum (FQ) approach, and asymptotic model approaches based on the linear combination of atomic orbitals (LCAO) and semi-empirical (SE) methods for deriving couplings. It is found that FQ calculations compare best overall with the experiments, followed by the LCAO, and the SE approaches. The experimental results together with the theoretical calculations, allow us to investigate the effects on modelled spectra and derived abundances and their uncertainties arising from uncertainties in the MN rates. Numerical experiments in a large grid of 1D model atmospheres, and a smaller set of 3D models, indicate that neglect of MN can lead to abundance errors of up to 0.1 dex (26\%) for Li at low metallicity, and 0.2 dex (58\%) for Na at high metallicity, while the uncertainties in the relevant MN rates as constrained by experiments correspond to uncertainties in abundances of much less than 0.01~dex (2\%). This agreement for simple atoms gives confidence in the FQ, LCAO and SE model approaches to be able to predict MN with the accuracy required for non-LTE modelling in stellar atmospheres., Comment: Accepted by ApJ

The ensemble of chemical element abundance measurements for stars, along with precision distances and orbit properties, provides high-dimensional data to study the evolution of the Milky Way. With this third data release of the Galactic Archaeology with HERMES (GALAH) survey, we publish 678 423 spectra for 588 571 mostly nearby stars (81.2% of stars are within <2 kpc), observed with the HERMES spectrograph at the Anglo-Australian Telescope. This release (hereafter GALAH+ DR3) includes all observations from GALAH Phase 1 (bright, main, and faint survey, 70%), K2-HERMES (17%), TESS-HERMES (5%), and a subset of ancillary observations (8%) including the bulge and >75 stellar clusters. We derive stellar parameters $T_\text{eff}$, $\log g$, [Fe/H], $v_\text{mic}$, $v_\text{broad}$ & $v_\text{rad}$ using our modified version of the spectrum synthesis code Spectroscopy Made Easy (SME) and 1D MARCS model atmospheres. We break spectroscopic degeneracies in our spectrum analysis with astrometry from $Gaia$ DR2 and photometry from 2MASS. We report abundance ratios [X/Fe] for 30 different elements (11 of which are based on non-LTE computations) covering five nucleosynthetic pathways. We describe validations for accuracy and precision, flagging of peculiar stars/measurements and recommendations for using our results. Our catalogue comprises 65% dwarfs, 34% giants, and 1% other/unclassified stars. Based on unflagged chemical composition and age, we find 62% young low-$\alpha$, 9% young high-$\alpha$, 27% old high-$\alpha$, and 2% stars with $\mathrm{[Fe/H]} \leq -1$. Based on kinematics, 4% are halo stars. Several Value-Added-Catalogues, including stellar ages and dynamics, updated after $Gaia$ eDR3, accompany this release and allow chrono-chemodynamic analyses, as we showcase., Comment: 48 pages, 27 figures, 5 pages appendix, accepted for publication in MNRAS, comments welcome, GALAH+ DR3 can be found at https://docs.datacentral.org.au/galah/dr3/overview/, Catalogs can be downloaded from https://cloud.datacentral.org.au/teamdata/GALAH/public/GALAH_DR3/

In this work we combine spectroscopic information from the \textit{SkyMapper survey for Extremely Metal-Poor stars} and astrometry from Gaia DR2 to investigate the kinematics of a sample of 475 stars with a metallicity range of $ -6.5 \leq \rm [Fe/H] \leq -2.05$ dex. Exploiting the action map, we identify 16 and 40 stars dynamically consistent with the \textit{Gaia Sausage} and \textit{Gaia Sequoia} accretion events, respectively. The most metal-poor of these candidates have metallicities of $\rm [Fe/H]=-3.31$ and $\rm [Fe/H]=-3.74$, respectively, helping to define the low-metallicity tail of the progenitors involved in the accretion events. We also find, consistent with other studies, that $\sim$21\% of the sample have orbits that remain confined to within 3~kpc of the Galactic plane, i.e., |Z$_{max}$| $\leq$ 3~kpc. Of particular interest is a sub-sample ($\sim$11\% of the total) of low |Z$_{max}$| stars with low eccentricities and prograde motions. The lowest metallicity of these stars has [Fe/H] = --4.30 and the sub-sample is best interpreted as the very low-metallicity tail of the metal-weak thick disk population. The low |Z$_{max}$|, low eccentricity stars with retrograde orbits are likely accreted, while the low |Z$_{max}$|, high eccentricity pro- and retrograde stars are plausibly associated with the \textit{Gaia Sausage} system. We find that a small fraction of our sample ($\sim$4\% of the total) is likely escaping from the Galaxy, and postulate that these stars have gained energy from gravitational interactions that occur when infalling dwarf galaxies are tidally disrupted., Comment: 23 pages, 15 figures, accepted for publication in MNRAS

Accurately known stellar lithium abundances may be used to shed light on a variety of astrophysical phenomena such as Big Bang nucleosynthesis, radial migration, ages of stars and stellar clusters, and planet engulfment events. We present a grid of synthetic lithium spectra that are computed in non-local thermodynamic equilibrium (NLTE) across the STAGGER grid of three-dimensional (3D) hydrodynamic stellar atmosphere models. This grid covers three Li lines at 610.4 nm, 670.8 nm, and 812.6 nm for stellar parameters representative of FGK-type dwarfs and giants, spanning $T_{\rm{eff}}=4000$-7000 K, $\log g=1.5$-5.0, $[\rm{Fe}/\rm{H}] = -4.0$-0.5, and $\textrm{A(Li)} = -0.5$-4.0. We find that our abundance corrections are up to 0.15 dex more negative than in previous work, due to a previously overlooked NLTE effect of blocking of UV lithium lines by background opacities, which has important implications for a wide range of science cases. We derive a new 3D NLTE solar abundance of $\textrm{A(Li)} = 0.96 \pm 0.05$, which is 0.09 dex lower than the commonly used value. We make our grids of synthetic spectra and abundance corrections publicly available through the Breidablik package. This package includes methods for accurately interpolating our grid to arbitrary stellar parameters through methods based on Kriging (Gaussian process regression) for line profiles, and MLP (Multi-Layer Perceptrons, a class of fully connected feedforward neural networks) for NLTE corrections and 3D NLTE abundances from equivalent widths, achieving interpolation errors of the order 0.01 dex., Comment: 20 pages, 12 figures, accepted for publication in MNRAS

Massive sets of stellar spectroscopic observations are rapidly becoming available and these can be used to determine the chemical composition and evolution of the Galaxy with unprecedented precision. One of the major challenges in this endeavour involves constructing realistic models of stellar spectra with which to reliably determine stellar abundances. At present, large stellar surveys commonly use simplified models that assume that the stellar atmospheres are approximately in local thermodynamic equilibrium (LTE). To test and ultimately relax this assumption, we have performed non-LTE calculations for $13$ different elements (H, Li, C, N, O, Na, Mg, Al, Si, K, Ca, Mn, and Ba), using recent model atoms that have physically-motivated descriptions for the inelastic collisions with neutral hydrogen, across a grid of $3756$ 1D MARCS model atmospheres that spans $3000\leq T_{\mathrm{eff}}/\mathrm{K}\leq8000$, $-0.5\leq\log{g/\mathrm{cm\,s^{-2}}}\leq5.5$, and $-5\leq\mathrm{[Fe/H]}\leq1$. We present the grids of departure coefficients that have been implemented into the GALAH DR3 analysis pipeline in order to complement the extant non-LTE grid for iron. We also present a detailed line-by-line re-analysis of $50126$ stars from GALAH DR3. We found that relaxing LTE can change the abundances by between $-0.7\,\mathrm{dex}$ and $+0.2\,\mathrm{dex}$ for different lines and stars. Taking departures from LTE into account can reduce the dispersion in the $\mathrm{[A/Fe]}$ versus $\mathrm{[Fe/H]}$ plane by up to $0.1\,\mathrm{dex}$, and it can remove spurious differences between the dwarfs and giants by up to $0.2\,\mathrm{dex}$. The resulting abundance slopes can thus be qualitatively different in non-LTE, possibly with important implications for the chemical evolution of our Galaxy., Comment: 19 pages, 25 figures, 2 tables, arXiv abstract abridged; accepted for publication in A&A

An unprecedented number of exoplanets are being discovered by the Transiting Exoplanet Survey Satellite (TESS). Determining the orbital parameters of these exoplanets, and especially their mass and radius, will depend heavily upon the measured physical characteristics of their host stars. We have cross-matched spectroscopic, photometric, and astrometric data from GALAH Data Release 2, the TESS Input Catalog and Gaia Data Release 2, to create a curated, self-consistent catalog of physical and chemical properties for 47,285 stars. Using these data we have derived isochrone masses and radii that are precise to within 5\%. We have revised the parameters of three confirmed, and twelve candidate, TESS planetary systems. These results cast doubt on whether CTOI-20125677 is indeed a planetary system since the revised planetary radii are now comparable to stellar sizes. Our GALAH-TESS catalog contains abundances for up to 23 elements. We have specifically analysed the molar ratios for C/O, Mg/Si, Fe/Si and Fe/Mg, to assist in determining the composition and structure of planets with $R_p < 4R_\oplus$. From these ratios, 36% fall within 2 sigma of the Sun/Earth values, suggesting that these stars may host rocky exoplanets with geological compositions similar to planets found within our own Solar system., Comment: Submitted to MNRAS, comments from the community welcomed. GALAH-TESS catalog is available upon request from the corresponding authour

Lithium depletion and enrichment in the cosmos is not yet well understood. To help tighten constraints on stellar and Galactic evolution models, we present the largest high-resolution analysis of Li abundances A(Li) to date, with results for over 100 000 GALAH field stars spanning effective temperatures $5900\,\mathrm{K} \lesssim \rm{T_{eff}} \lesssim7000\,\mathrm{K}$ and metallicities $-3 \lesssim \rm[Fe/H] \lesssim +0.5$. We separated these stars into two groups, on the warm and cool side of the so-called Li-dip, a localised region of the Kiel diagram wherein lithium is severely depleted. We discovered that stars in these two groups show similar trends in the A(Li)-[Fe/H] plane, but with a roughly constant offset in A(Li) of 0.4 dex, the warm group having higher Li abundances. At $\rm[Fe/H]\gtrsim-0.5$, a significant increasing in Li abundance with increasing metallicity is evident in both groups, signalling the onset of significant Galactic production. At lower metallicity, stars in the cool group sit on the Spite plateau, showing a reduced lithium of around 0.4 dex relative to the primordial value predicted from Big Bang nucleosynthesis (BBN). However, stars in the warm group between [Fe/H] = -1.0 and -0.5, form an elevated plateau that is largely consistent with the BBN prediction. This may indicate that these stars in fact preserve the primordial Li produced in the early Universe., Comment: 5 pages, 2 figures, Accepted for publication by MNRAS

Nitrogen is an important element in various fields of stellar and Galactic astronomy, and the solar nitrogen abundance is crucial as a yardstick for comparing different objects in the cosmos. In order to obtain a precise and accurate value for this abundance, we carried out N i line formation calculations in a 3D radiative-hydrodynamic STAGGER model solar atmosphere, in full 3D non-local thermodynamic equilibrium (non-LTE), using a model atom that includes physically-motivated descriptions for the inelastic collisions of N i with free electrons and with neutral hydrogen. We selected five N i lines of high excitation energy to study in detail, based on their strengths and on their being relatively free of blends. We found that these lines are slightly strengthened from non-LTE photon losses and from 3D granulation effects, resulting in negative abundance corrections of around $-0.01$ dex and $-0.04$ dex respectively. Our advocated solar nitrogen abundance is $\log\epsilon_{\mathrm{N}} = 7.77$, with the systematic $1\sigma$ uncertainty estimated to be $0.05$ dex. This result is consistent with earlier studies after correcting for differences in line selections and equivalent widths., Comment: 11 pages, 7 figures, 3 tables; accepted for publication in Astronomy & Astrophysics

The abundances of carbon, oxygen, and iron in late-type stars are important parameters in exoplanetary and stellar physics, as well as key tracers of stellar populations and Galactic chemical evolution. We carried out three-dimensional (3D) non-LTE radiative transfer calculations for CI and OI, and 3D LTE radiative transfer calculations for FeII, across the STAGGER-grid of 3D hydrodynamic model atmospheres. The absolute 3D non-LTE versus 1D LTE abundance corrections can be as severe as $-0.3$ dex for CI lines in low-metallicity F dwarfs, and $-0.6$ dex for OI lines in high-metallicity F dwarfs. The 3D LTE versus 1D LTE abundance corrections for FeII lines are less severe, typically less than $+0.15$ dex. We used the corrections in a re-analysis of carbon, oxygen, and iron in $187$ F and G dwarfs in the Galactic disk and halo. Applying the differential 3D non-LTE corrections to 1D LTE abundances visibly reduces the scatter in the abundance plots. The thick disk and high-$\alpha$ halo population rise in carbon and oxygen with decreasing metallicity, and reach a maximum of [C/Fe]$\approx0.2$ and a plateau of [O/Fe]$\approx0.6$ at [Fe/H]$\approx-1.0$. The low-$\alpha$ halo population is qualitatively similar, albeit offset towards lower metallicities and with larger scatter. Nevertheless, these populations overlap in the [C/O] versus [O/H] plane, decreasing to a plateau of [C/O]$\approx-0.6$ below [O/H]$\approx-1.0$. In the thin-disk, stars having confirmed planet detections tend to have higher values of C/O at given [O/H]; this potential signature of planet formation is only apparent after applying the abundance corrections to the 1D LTE results. Our grids of line-by-line abundance corrections are publicly available and can be readily used to improve the accuracy of spectroscopic analyses of late-type stars., Comment: 20 pages, 14 figures, 7 online tables; arXiv abstract abridged; published in Astronomy & Astrophysics

Older GCE models predict [K/Fe] ratios as much as 1 dex lower than those inferred from stellar observations. Abundances of potassium are mainly based on analyses of the 7698 $\AA$ resonance line, and the discrepancy between models and observations is in part caused by the LTE assumption. We study the statistical equilibrium of KI, focusing on the non-LTE effects on the $7698 \ \AA$ line. We aim to determine how non-LTE abundances of K can improve the analysis of its chemical evolution, and help to constrain the yields of models. We construct a model atom that employs the most up-to-date data. In particular, we calculate and present inelastic e+K collisional excitation cross-sections from the convergent close-coupling and the $B$-Spline $R$-matrix methods, and H+K collisions from the two-electron model. We constructed a fine grid of non-LTE abundance corrections that span $4000<\teff / \rm{K}<8000$, $0.50<\lgg<5.00$, $-5.00<\feh<+0.50$, and applied the corrections to abundances from the literature. In concordance with previous studies, we find severe non-LTE effects in the $7698 \ \AA$ line, which is stronger in non-LTE with abundance corrections that can reach $\sim-0.7\,\dex$. We explore the effects of atmospheric inhomogeneity by computing a full 3D non-LTE stellar spectrum of KI for a test star. We find that 3D is necessary to predict a correct shape of the resonance 7698 $\AA$ line, but the line strength is similar to that found in 1D non-LTE. Our non-LTE abundance corrections reduce the scatter and change the cosmic trends of literature K abundances. In the regime [Fe/H]$\lesssim-1.0$ the non-LTE abundances show a good agreement with the GCE model with yields from rotating massive stars. The reduced scatter of the non-LTE corrected abundances of a sample of solar twins shows that line-by-line differential analysis techniques cannot fully compensate for systematic modelling errors., Comment: 18 pages, 15 Figures

Low-energy inelastic collisions with neutral hydrogen atoms are important processes in stellar atmospheres, and a persistent source of uncertainty in non-LTE modelling of stellar spectra. We have calculated and studied excitation and charge transfer of C i and of N i due to such collisions. We used a previously presented method that is based on an asymptotic two-electron linear combination of atomic orbitals (LCAO) model of ionic-covalent interactions for the adiabatic potential energies, combined with the multichannel Landau-Zener model for the collision dynamics. We find that charge transfer processes typically lead to much larger rate coefficients than excitation processes do, consistent with studies of other atomic species. Two-electron processes were considered and lead to non-zero rate coefficients that can potentially impact statistical equilibrium calculations. However, they were included in the model in an approximate way, via an estimate for the two-electron coupling that was presented earlier in the literature: the validity of these data should be checked in a future work., Comment: 11 pages, 5 figures, 4 tables; accepted for publication in A&A

Carbon abundances in late-type stars are important in a variety of astrophysical contexts. However C i lines, one of the main abundance diagnostics, are sensitive to departures from local thermodynamic equilibrium (LTE). We present a model atom for non-LTE analyses of C i lines, that uses a new, physically-motivated recipe for the rates of neutral hydrogen impact excitation. We analyse C i lines in the solar spectrum, employing a three-dimensional (3D) hydrodynamic model solar atmosphere and 3D non-LTE radiative transfer. We find negative non-LTE abundance corrections for C i lines in the solar photosphere, in accordance with previous studies, reaching up to around 0.1 dex in the disk-integrated flux. We also present the first fully consistent 3D non-LTE solar carbon abundance determination: we infer log $\epsilon_{\text{C}}$ = $8.44\pm0.02$, in good agreement with the current standard value. Our models reproduce the observed solar centre-to-limb variations of various C i lines, without any adjustments to the rates of neutral hydrogen impact excitation, suggesting that the proposed recipe may be a solution to the long-standing problem of how to reliably model inelastic collisions with neutral hydrogen in late-type stellar atmospheres., Comment: 14 pages, 8 figures, 2 tables; accepted for publication in A&A

Oxygen abundance measurements are important for understanding stellar structure and evolution. Measured in Cepheids, they further provide clues on the metallicity gradient and chemo-dynamical evolution in the Galaxy. However, most of the abundance analyses of Cepheids to date have been based on one-dimensional (1D) hydrostatic model atmospheres. Here, we test the validity of this approach for the key oxygen abundance diagnostic, the \ion{O}{I} $777\,\mathrm{nm}$~triplet lines. We carry out 2D non-LTE radiative transfer clculations across two different 2D radiation hydrodynamics simulations of Cepheid atmospheres, having stellar parameters of $T_\mathrm{eff}= 5600$ K, solar chemical compositions, and $\log\,g= 1.5$ and $2.0$, corresponding to pulsation periods of 9 and 3 days, respectively. We find that the 2D non-LTE versus 1D LTE abundance differences range from $-1.0$~dex to $-0.25$~dex depending on pulsational phase. The 2D non-LTE versus 1D non-LTE abundance differences range from $-0.2$~dex to $0.8$~dex. The abundance differences are smallest when the Cepheid atmospheres are closest to hydrostatic equilibrium, corresponding to phases of around $0.3$ to $0.8$, and we recommend these phases for observers deriving the oxygen abundance from \ion{O}{I} $777\,\mathrm{nm}$ triplet with 1D hydrostatic models., Comment: 9 pages, 10 figures; Published in Astronomy and Astrophysics

Carbon and oxygen are key tracers of the Galactic chemical evolution; in particular, a reported upturn in [C/O] towards decreasing [O/H] in metal-poor halo stars could be a signature of nucleosynthesis by massive Population III stars. We reanalyse carbon, oxygen, and iron abundances in thirty-nine metal-poor turn-off stars. For the first time, we take into account three-dimensional (3D) hydrodynamic effects together with departures from local thermodynamic equilibrium (LTE) when determining both the stellar parameters and the elemental abundances, by deriving effective temperatures from 3D non-LTE H$\beta$ profiles, surface gravities from Gaia parallaxes, iron abundances from 3D LTE Feii equivalent widths, and carbon and oxygen abundances from 3D non-LTE Ci and Oi equivalent widths. We find that [C/Fe] stays flat with [Fe/H], whereas [O/Fe] increases linearly up to $0.75$ dex with decreasing [Fe/H] down to $-3.0$ dex. As such [C/O] monotonically decreases towards decreasing [O/H], in contrast to previous findings, mainly by virtue of less severe non-LTE effects for Oi at low [Fe/H] with our improved calculations., Comment: 5 pages, 2 figures; published in A&A Letters

Several recent spectroscopic investigations have presented conflicting results on the existence of Na-rich asymptotic giant branch (AGB) stars in the Galactic globular cluster M4 (NGC6121). The studies disagree on whether or not Na-rich red giant branch (RGB) stars evolve to the AGB. For a sample of previously published HER- MES/AAT AGB and RGB stellar spectra we present a re-analysis of O, Na, and Fe abundances, and a new analysis of Mg and Al abundances; we also present CN band strengths for this sample, derived from low-resolution AAOmega spectra. Following a detailed literature comparison, we find that the AGB samples of all studies consistently show lower abundances of Na and Al, and are weaker in CN, than RGB stars in the cluster. This is similar to recent observations of AGB stars in NGC 6752 and M 62. In an attempt to explain this result, we present new theoretical stellar evolutionary models for M 4; however, these predict that all stars, including Na-rich RGB stars, evolve onto the AGB. We test the robustness of our abundance results using a variety of atmospheric models and spectroscopic methods; however, we do not find evidence that systematic modelling uncertainties can explain the apparent lack of Na- rich AGB stars in M4. We conclude that an unexplained, but robust, discordance between observations and theory remains for the AGB stars in M 4., Comment: 25 pages, 21 figures, 15 tables, accepted for publication in MNRAS

Open cluster members are coeval and share the same initial bulk chemical compositions. Consequently, differences in surface abundances between members of a cluster that are at different evolutionary stages can be used to study the effects of mixing and internal chemical processing. We carry out an abundance analysis of seven elements (Li, O, Na, Mg, Al, Si, Fe) in 66 stars belonging to the open cluster M67, based on high resolution GALAH spectra, 1D MARCS model atmospheres, and, for the first time for a large sample of M67 stars, non-local thermodynamic equilibrium (non-LTE) radiative transfer. From the non-LTE analysis, we find a typical star-to-star scatter in the abundance ratios of around 0.05 dex; this scatter is slightly but systematically larger when LTE is assumed instead. We find trends in the abundance ratios with effective temperature, indicating systematic differences in the surface abundances between turn-off and giant stars; these trends are more pronounced when LTE is assumed. However, in the non-LTE analysis, most of the element trends have been flattened. Two species are exceptions to this behaviour, namely Al and Si, which both clearly display remaining trends in the non-LTE analysis. We comment on the possible origin of these trends, by comparing them with recent stellar models that include atomic diffusion., Comment: 14 pages, 10 figures; resubmitted to MNRAS

The Galactic Archaeology with HERMES (GALAH) survey is a large-scale stellar spectroscopic survey of the Milky Way and designed to deliver chemical information complementary to a large number of stars covered by the $Gaia$ mission. We present the GALAH second public data release (GALAH DR2) containing 342,682 stars. For these stars, the GALAH collaboration provides stellar parameters and abundances for up to 23 elements to the community. Here we present the target selection, observation, data reduction and detailed explanation of how the spectra were analysed to estimate stellar parameters and element abundances. For the stellar analysis, we have used a multi-step approach. We use the physics-driven spectrum synthesis of Spectroscopy Made Easy (SME) to derive stellar labels ($T_\mathrm{eff}$, $\log g$, $\mathrm{[Fe/H]}$, $\mathrm{[X/Fe]}$, $v_\mathrm{mic}$, $v \sin i$, $A_{K_S}$) for a representative training set of stars. This information is then propagated to the whole survey with the data-driven method of $The~Cannon$. Special care has been exercised in the spectral synthesis to only consider spectral lines that have reliable atomic input data and are little affected by blending lines. Departures from local thermodynamic equilibrium (LTE) are considered for several key elements, including Li, O, Na, Mg, Al, Si, and Fe, using 1D MARCS stellar atmosphere models. Validation tests including repeat observations, Gaia benchmark stars, open and globular clusters, and K2 asteroseismic targets lend confidence in our methods and results. Combining the GALAH DR2 catalogue with the kinematic information from $Gaia$ will enable a wide range of Galactic Archaeology studies, with unprecedented detail, dimensionality, and scope., Comment: 34+7 pages, 31 Figures, 5 tables, 1 catalog, submitted to MNRAS

Hydrogen Balmer lines are commonly used as spectroscopic effective temperature diagnostics of late-type stars. However, the absolute accuracy of classical methods that are based on one-dimensional (1D) hydrostatic model atmospheres and local thermodynamic equilibrium (LTE) is still unclear. To investigate this, we carry out 3D non-LTE calculations for the Balmer lines, performed, for the first time, over an extensive grid of 3D hydrodynamic STAGGER model atmospheres. For H$\alpha$, H$\beta$, and H$\gamma$, we find significant 1D non-LTE versus 3D non-LTE differences (3D effects): the outer wings tend to be stronger in 3D models, particularly for H$\gamma$, while the inner wings can be weaker in 3D models, particularly for H$\alpha$. For H$\alpha$, we also find significant 3D LTE versus 3D non-LTE differences (non-LTE effects): in warmer stars ($T_{\text{eff}}\approx6500$K) the inner wings tend to be weaker in non-LTE models, while at lower effective temperatures ($T_{\text{eff}}\approx4500$K) the inner wings can be stronger in non-LTE models; the non-LTE effects are more severe at lower metallicities. We test our 3D non-LTE models against observations of well-studied benchmark stars. For the Sun, we infer concordant effective temperatures from H$\alpha$, H$\beta$, and H$\gamma$; however the value is too low by around 50K which could signal residual modelling shortcomings. For other benchmark stars, our 3D non-LTE models generally reproduce the effective temperatures to within $1\sigma$ uncertainties. For H$\alpha$, the absolute 3D effects and non-LTE effects can separately reach around 100K, in terms of inferred effective temperatures. For metal-poor turn-off stars, 1D LTE models of H$\alpha$ can underestimate effective temperatures by around 150K. Our 3D non-LTE model spectra are publicly available, and can be used for more reliable spectroscopic effective temperature determinations., Comment: 19 pages, 10 figures, abstract abridged; accepted for publication in Astronomy & Astrophysics