Showing 1–50 of 57 results for author: Amarsi, A

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

A distinct halo population revealed from 3D non-LTE magnesium abundances

Authors: T. Matsuno, A. M. Amarsi, M. Carlos, P. E. Nissen

Abstract: 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 loca… ▽ More 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, $α$-rich halo, and $α$-poor halo. We find that the 3D non-LTE results show a richer substructure in [Mg/Fe]-[Fe/H] in the $α$-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. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: accepted to A\&A, 8 pages

A conclusive non-detection of magnetic field in the Am star o Peg with high-precision near-infrared spectroscopy

Authors: O. Kochukhov, A. M. Amarsi, A. Lavail, H. L. Ruh, A. Hahlin, A. Hatzes, E. Nagel, N. Piskunov, K. Pouilly, A. Reiners, M. Rengel, U. Seemann, D. Shulyak

Abstract: 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-mas… ▽ More 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. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: 11 pages, 8 figures; accepted for publication in A&A

An extended and refined grid of 3D STAGGER model atmospheres. Processed snapshots for stellar spectroscopy

Authors: Luisa F. Rodríguez Díaz, Cis Lagae, Anish M. Amarsi, Lionel Bigot, Yixiao Zhou, Víctor Aguirre Børsen-Koch, Karin Lind, Regner Trampedach, Remo Collet

Abstract: 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 work… ▽ More 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. △ Less

Submitted 13 May, 2024; originally announced May 2024.

On the Pair-Instability Supernova origin of J1010+2358

Authors: Ása Skúladóttir, Ioanna Koutsouridou, Irene Vanni, Anish M. Amarsi, Romain Lucchesi, Stefania Salvadori, David Aguado

Abstract: 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… ▽ More 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 $(χ^2/χ^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. △ Less

Submitted 3 May, 2024; v1 submitted 29 April, 2024; originally announced April 2024.

Comments: Accepted to ApJ Letters (3rd of May 2024). 1 table, 4 figures, and Appendix

Spectroscopic characterisation of gravitationally lensed stars at high redshifts

Authors: Emma Lundqvist, Erik Zackrisson, Calum Hawcroft, Anish M. Amarsi, Brian Welch

Abstract: 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… ▽ More 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 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 $λ$ 1550 Å line and the Si IV $λλ$1393, 1403 Å-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. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: 11 pages, 6 figures, submitted to A&A. Comments are welcome

In-depth analysis of solar models with high-metallicity abundances and updated opacity tables

Authors: G. Buldgen, A. Noels, R. Scuflaire, A. M. Amarsi, N. Grevesse, P. Eggenberger, J. Colgan, C. J. Fontes, V. A. Baturin, A. V. Oreshina, S. V. Ayukov, P. Hakel, D. P. Kilcrease

Abstract: 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… ▽ More 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. △ Less

Submitted 16 April, 2024; originally announced April 2024.

Comments: Accepted for publication in Astronomy and Astrophysics, version before language editing

Titanium abundances in late-type stars, II. Grid of departure coefficients and application to a sample of $70\,000$ stars

Authors: J. W. E. Mallinson, K. Lind, A. M. Amarsi, K. Youakim

Abstract: 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 wit… ▽ More 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, $β_ν$, 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 $Δ$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. △ Less

Submitted 28 March, 2024; originally announced March 2024.

Targeted optimization in small-scale atomic structure calculations: application to Au I

Authors: Sema Caliskan, Jon Grumer, Anish M. Amarsi

Abstract: 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 comp… ▽ More 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. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 14 pages, 3 figures, 6 tables; published in Journal of Physics B

Journal ref: 2024, J. Phys. B: At. Mol. Opt. Phys., 57, 055003

3D non-LTE abundance analyses of late-type stars

Authors: Karin Lind, Anish Mayur Amarsi

Abstract: 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 abundances 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-i… ▽ More 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 abundances 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 compared to traditional one-dimensional (1D) hydrostatic models. As evident from a multitude of observational tests, the coupling of 3D RHD models with line-formation in non-local 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: 1) 3D and non-LTE effects are intricately coupled and consistent modelling thereof is necessary for high-precision abundances, which is currently feasible for individual elements in large surveys. Mean 3D (<3D>) models are not adequate as substitutes. 2) The solar abundance debate is presently dominated by choices and systematic uncertainties that are not specific to 3D non-LTE modelling. 3) 3D non-LTE abundance corrections have a profound impact on our understanding of FGK-type stars, exoplanets, and the nucleosynthetic origins of the elements. △ Less

Submitted 1 January, 2024; originally announced January 2024.

Comments: To appear in Annual Reviews of Astronomy and Astrophysics (65 pages, 13 figures)

Abundances of iron-peak elements in accreted and in situ born Galactic halo stars

Authors: P. E. Nissen, A. M. Amarsi, Á. Skúladóttir, W. J. Schuster

Abstract: 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 e… ▽ More 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. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Comments: Accepted for publication in A&A. 12 pages, 4 tables, and 7 figures

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

Authors: G. Canocchi, K. Lind, C. Lagae, A. G. M. Pietrow, A. M. Amarsi, D. Kiselman, O. Andriienko, H. J. Hoeijmakers

Abstract: 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 tha… ▽ More 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 $μ=0.1$ and $μ=1.0$. Results. Our 3D non-LTE modeling of the Na I D$_1$ 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. 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. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Comments: 23 pages, accepted for publication by A&A

Helioseismic determination of the solar metal mass fraction

Authors: G. Buldgen, A. Noels, V. A. Baturin, A. V. Oreshina, S. V. Ayukov, R. Scuflaire, A. M. Amarsi, N. Grevesse

Abstract: 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… ▽ More 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. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: Accepted for publication in A&A, online materials available on demand or on publisher website

3D Stagger model atmospheres with FreeEOS I. Exploring the impact of microphysics on the Sun

Authors: Yixiao Zhou, Anish M. Amarsi, Victor Aguirre Børsen-Koch, Klara G. Karlsmose, Remo Collet, Thomas Nordlander

Abstract: 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 t… ▽ More 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 $α$-element abundance, carbon-enhanced metal-poor stars and population II stars with peculiar chemical compositions using 3D Stagger model atmospheres. △ Less

Submitted 11 July, 2023; originally announced July 2023.

Comments: 24 pages, 20 figures, accepted for publication in A&A

Journal ref: A&A 677, A98 (2023)

Extended atomic data for oxygen abundance analyses

Authors: W. Li, P. Jönsson, A. M. Amarsi, M. C. Li, J. Grumer

Abstract: 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 r… ▽ More 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ε_{\mathrm{O}}=8.70\pm0.04$. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Comments: 13 pages, 5 figures; Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 674, A54 (2023)

Surface structure of 45 Hercules: An otherwise unremarkable Ap star with a surprisingly weak magnetic field

Authors: O. Kochukhov, H. Gürsoytrak Mutlay, A. M. Amarsi, P. Petit, I. Mutlay, B. Gürol

Abstract: 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 a… ▽ More 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. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 15 pages, 19 figures; Accepted for publication in MNRAS

Raising the observed metallicity floor with a 3D non-LTE analysis of SDSS J102915.14+172927.9

Authors: C. Lagae, A. M. Amarsi, L. F. Rodríguez Díaz, K. Lind, T. Nordlander, T. T. Hansen, A. Heger

Abstract: 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… ▽ More 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. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: 14 pages, 9 figures, Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 672, A90 (2023)

Extended MCDHF calculations of energy levels and transition data for N I

Authors: Meichun Li, Wenxian Li, P. Jönsson, Anish M. Amarsi, Jon Grumer

Abstract: 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 method… ▽ More 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%). △ Less

Submitted 31 January, 2023; originally announced February 2023.

Comments: 17 pages, 7 figures; Accepted for publication in The Astrophysical Journal Supplement Series

PISN-explorer: hunting the descendants of very massive first stars

Authors: D. S. Aguado, S. Salvadori, A. Skúladóttir, E. Caffau, P. Bonifacio, I. Vanni, V. Gelli, I. Koutsouridou, A. M. Amarsi

Abstract: 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… ▽ More 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. △ Less

Submitted 9 January, 2023; originally announced January 2023.

Comments: 13 pages, 7 figures Accepted for publication in MNRAS

Higher metal abundances do not solve the solar problem

Authors: G. Buldgen, P. Eggenberger, A. Noels, R. Scuflaire, A. M. Amarsi, N. Grevesse, S. Salmon

Abstract: 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… ▽ More 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. △ Less

Submitted 13 December, 2022; originally announced December 2022.

Comments: Accepted for publication in Astronomy and Astrophysics

Journal ref: A&A 669, L9 (2023)

The chemical compositions of multiple stellar populations in the globular cluster NGC 2808

Authors: M. Carlos, A. F. Marino, A. P. Milone, E. Dondoglio, S. Jang, M. V. Legnardi, A. Mohandasan, G. Cordoni, E. P. Lagioia, A. M. Amarsi, H. Jerjen

Abstract: 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 stell… ▽ More 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é star. △ Less

Submitted 2 December, 2022; originally announced December 2022.

Comments: Accepted for publication in MNRAS. 19 pages and 13 figures

The impact of carbon and oxygen abundances on the metal-poor initial mass function

Authors: Piyush Sharda, Anish M. Amarsi, Kathryn Grasha, Mark R. Krumholz, David Yong, Gen Chiaki, Arpita Roy, Thomas Nordlander

Abstract: 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… ▽ More 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. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 14 pages, 8 figures. Accepted by MNRAS

Titanium abundances in late-type stars I. 1D non-local thermodynamic equilibrium modelling in benchmark dwarfs and giants

Authors: J. W. E. Mallinson, K. Lind, A. M. Amarsi, P. S. Barklem, J. Grumer, A. K. Belyaev, K. Youakim

Abstract: 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 bo… ▽ More 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. △ Less

Submitted 7 February, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 9 pages plus appendix, 6 figures; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 668, A103 (2022)

3D non-LTE iron abundances in FG-type dwarfs

Authors: A. M. Amarsi, S. Liljegren, P. E. Nissen

Abstract: 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 abundance… ▽ More 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. △ Less

Submitted 9 December, 2022; v1 submitted 27 September, 2022; originally announced September 2022.

Comments: 17 pages, 11 figures, 5 tables; arXiv abstract abridged; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 668, A68 (2022)

The solar carbon, nitrogen, and oxygen abundances from a 3D LTE analysis of molecular lines

Authors: A. M. Amarsi, N. Grevesse, M. Asplund, R. Collet

Abstract: 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… ▽ More 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ε_{C} = 8.47\pm0.02$, $\logε_{N} = 7.89\pm0.04$, and $\logε_{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. △ Less

Submitted 10 September, 2021; originally announced September 2021.

Comments: 12 pages, 6 figures, 3 tables; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 656, A113 (2021)

Zero-metallicity hypernova uncovered by an ultra metal-poor star in the Sculptor dwarf spheroidal galaxy

Authors: Ása Skúladóttir, Stefania Salvadori, Anish M. Amarsi, Eline Tolstoy, Michael J. Irwin, Vanessa Hill, Pascale Jablonka, Giuseppina Battaglia, Else Starkenburg, Davide Massari, Amina Helmi, Lorenzo Posti

Abstract: 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 e… ▽ More 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 $α$-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. △ Less

Submitted 24 June, 2021; v1 submitted 22 June, 2021; originally announced June 2021.

Comments: Accepted in ApJL, 7 pages, 3 figures, 1 table, and an appendix (1 figure and 3 tables)

The chemical make-up of the Sun: A 2020 vision

Authors: M. Asplund, A. M. Amarsi, N. Grevesse

Abstract: 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… ▽ More 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ε_{\text{C}}=8.46\pm0.04$, $\logε_{\text{N}}=7.83\pm0.07$, and $\logε_{\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. △ Less

Submitted 14 July, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 31 pages, arXiv abstract abridged; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 653, A141 (2021)

The relationship between photometric and spectroscopic oscillation amplitudes from 3D stellar atmosphere simulations

Authors: Yixiao Zhou, Thomas Nordlander, Luca Casagrande, Meridith Joyce, Yaguang Li, Anish M. Amarsi, Henrique Reggiani, Martin Asplund

Abstract: 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… ▽ More 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 $ε$ 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 $ε$ 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. △ Less

Submitted 3 February, 2021; originally announced February 2021.

Comments: 16 pages, 7 figures, accepted for publication in MNRAS

Extended theoretical transition data in C I - IV

Authors: W. Li, A. M. Amarsi, A. Papoulia, J. Ekman, P. Jönsson

Abstract: 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… ▽ More 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. △ Less

Submitted 22 January, 2021; originally announced January 2021.

Comments: 21 pages, 4 figures

Mutual neutralisation in Li$^+$+H$^-$/D$^-$ and Na$^+$+H$^-$/D$^-$ collisions: Implications of experimental results for non-LTE modelling of stellar spectra

Authors: Paul S. Barklem, Anish M. Amarsi, Jon Grumer, Gustav Eklund, Stefan Rosén, MingChao Ji, Henrik Cederquist, Henning Zettergren, Henning T. Schmidt

Abstract: 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… ▽ More 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. △ Less

Submitted 22 December, 2020; originally announced December 2020.

Comments: Accepted by ApJ

The GALAH+ Survey: Third Data Release

Authors: Sven Buder, Sanjib Sharma, Janez Kos, Anish M. Amarsi, Thomas Nordlander, Karin Lind, Sarah L. Martell, Martin Asplund, Joss Bland-Hawthorn, Andrew R. Casey, Gayandhi M. De Silva, Valentina D'Orazi, Ken C. Freeman, Michael R. Hayden, Geraint F. Lewis, Jane Lin, Katharine J. Schlesinger, Jeffrey D. Simpson, Dennis Stello, Daniel B. Zucker, Tomaz Zwitter, Kevin L. Beeson, Tobias Buck, Luca Casagrande, Jake T. Clark , et al. (22 additional authors not shown)

Abstract: 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 spe… ▽ More 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-$α$, 9% young high-$α$, 27% old high-$α$, 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. △ Less

Submitted 28 April, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 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/

Exploring the Galaxy's halo and very metal-weak thick disk with SkyMapper and Gaia DR2

Authors: G. Cordoni, G. S. Da Costa, D. Yong, A. D. Mackey, A. F. Marino, S. Monty, T. Nordlander, J. E. Norris, M. Asplund, M. S. Bessell, A. R. Casey, A. Frebel, K. Lind, S. J. Murphy, B. P. Schmidt, X. D. Gao, T. Xylakis-Dornbusch, A. M. Amarsi, A. P. Milone

Abstract: 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 Sequoi… ▽ More 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. △ Less

Submitted 2 November, 2020; originally announced November 2020.

Comments: 23 pages, 15 figures, accepted for publication in MNRAS

3D NLTE spectral line formation of lithium in late-type stars

Authors: E. Wang, T. Nordlander, M. Asplund, A. M. Amarsi, K. Lind, Y. Zhou

Abstract: 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) hydrodynami… ▽ More 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. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: 20 pages, 12 figures, accepted for publication in MNRAS

The GALAH Survey: Non-LTE departure coefficients for large spectroscopic surveys

Authors: A. M. Amarsi, K. Lind, Y. Osorio, T. Nordlander, M. Bergemann, H. Reggiani, E. X. Wang, S. Buder, M. Asplund, P. S. Barklem, A. Wehrhahn, Á. Skúladóttir, C. Kobayashi, A. I. Karakas, X. D. Gao, J. Bland-Hawthorn, G. M. De Silva, J. Kos, G. F. Lewis, S. L. Martell, S. Sharma, J. D. Simpson, D. B. Zucker, K. Čotar, J. Horner , et al. (1 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 21 August, 2020; originally announced August 2020.

Comments: 19 pages, 25 figures, 2 tables, arXiv abstract abridged; accepted for publication in A&A

Journal ref: A&A 642, A62 (2020)

The GALAH Survey: Using Galactic Archaeology to Refine our Knowledge of TESS Target Stars

Authors: Jake T. Clark, Mathieu Clerte, Natalie R. Hinkel, Cayman T. Unterborn, Robert A. Wittenmyer, Jonathan Horner, Duncan J. Wright, Brad Carter, Timothy D. Morton, Lorenzo Spina, Martin Asplund, Sven Buder, Joss Bland-Hawthorn, Andy Casey, Gayandhi De Silva, Valentina D'Orazi, Ly Duong, Michael Hayden, Ken Freeman, Janez Kos, Geraint Lewis, Jane Lin, Karin Lind, Sarah Martell, Sanjib Sharma , et al. (7 additional authors not shown)

Abstract: 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 Inp… ▽ More 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. △ Less

Submitted 12 August, 2020; originally announced August 2020.

Comments: Submitted to MNRAS, comments from the community welcomed. GALAH-TESS catalog is available upon request from the corresponding authour

The GALAH Survey: A new constraint on cosmological lithium and Galactic lithium evolution from warm dwarf stars

Authors: Xudong Gao, Karin Lind, Anish M. Amarsi, Sven Buder, Joss Bland-Hawthorn, Simon W. Campbell, Martin Asplund, Andrew R. Casey, Gayandhi M. De Silva, Ken C. Freeman, Michael R. Hayden, Geraint F. Lewis, Sarah L. Martell, Jeffrey D. Simpson, Sanjib Sharma, Daniel B. Zucker, Tomaž Zwitter, Jonathan Horner, Ulisse Munari, Thomas Nordlander, Dennis Stello, Yuan-Sen Ting, Gregor Traven, Robert A. Wittenmyer, the GALAH collaboration

Abstract: 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… ▽ More 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. △ Less

Submitted 9 June, 2020; originally announced June 2020.

Comments: 5 pages, 2 figures, Accepted for publication by MNRAS

The 3D non-LTE solar nitrogen abundance from atomic lines

Authors: A. M. Amarsi, N. Grevesse, J. Grumer, M. Asplund, P. S. Barklem, R. Collet

Abstract: 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 thermodyn… ▽ More 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ε_{\mathrm{N}} = 7.77$, with the systematic $1σ$ uncertainty estimated to be $0.05$ dex. This result is consistent with earlier studies after correcting for differences in line selections and equivalent widths. △ Less

Submitted 31 March, 2020; v1 submitted 27 March, 2020; originally announced March 2020.

Comments: 11 pages, 7 figures, 3 tables; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 636, A120 (2020)

Carbon, oxygen, and iron abundances in disk and halo stars. Implications of 3D non-LTE spectral line formation

Authors: A. M. Amarsi, P. E. Nissen, Á. Skúladóttir

Abstract: 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 at… ▽ More 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-$α$ 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-$α$ 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. △ Less

Submitted 8 October, 2019; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: 20 pages, 14 figures, 7 online tables; arXiv abstract abridged; published in Astronomy & Astrophysics

Journal ref: A&A 630, A104 (2019)

Non-LTE analysis of K I in late-type stars

Authors: Henrique Reggiani, Anish M. Amarsi, Karin Lind, Paul S. Barklem, Oleg Zatsarinny, Klaus Bartschat, Dmitry V. Fursa, Igor Bray, Lorenzo Spina, Jorge Meléndez

Abstract: 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 $Å$ 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 \ Å$ line. We aim to det… ▽ More 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 $Å$ 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 \ Å$ 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 \ Å$ 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 $Å$ 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. △ Less

Submitted 19 June, 2019; originally announced June 2019.

Comments: 18 pages, 15 Figures

Journal ref: A&A 627, A177 (2019)

Excitation and charge transfer in low-energy hydrogen atom collisions with neutral carbon and nitrogen

Authors: A. M. Amarsi, P. S. Barklem

Abstract: 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 atom… ▽ More 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. △ Less

Submitted 23 May, 2019; v1 submitted 30 March, 2019; originally announced April 2019.

Comments: 11 pages, 5 figures, 4 tables; accepted for publication in A&A

Journal ref: A&A 625, A78 (2019)

3D non-LTE line formation of neutral carbon in the Sun

Authors: A. M. Amarsi, P. S. Barklem, R. Collet, N. Grevesse, M. Asplund

Abstract: 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 i… ▽ More 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 $ε_{\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. △ Less

Submitted 21 March, 2019; originally announced March 2019.

Comments: 14 pages, 8 figures, 2 tables; accepted for publication in A&A

Journal ref: A&A 624, A111 (2019)

4MOST: Project overview and information for the First Call for Proposals

Authors: R. S. de Jong, O. Agertz, A. Agudo Berbel, J. Aird, D. A. Alexander, A. Amarsi, F. Anders, R. Andrae, B. Ansarinejad, W. Ansorge, P. Antilogus, H. Anwand-Heerwart, A. Arentsen, A. Arnadottir, M. Asplund, M. Auger, N. Azais, D. Baade, G. Baker, S. Baker, E. Balbinot, I. K. Baldry, M. Banerji, S. Barden, P. Barklem , et al. (313 additional authors not shown)

Abstract: We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolut… ▽ More We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = λ/Δλ\sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations. △ Less

Submitted 1 April, 2019; v1 submitted 6 March, 2019; originally announced March 2019.

Comments: Part of the 4MOST issue of The Messenger, published in preparation of 4MOST Community Workshop, see http://www.eso.org/sci/meetings/2019/4MOST.html

Journal ref: The Messenger, 2019, 175, 3

Two-dimensional non-LTE \ion{O}{I} 777\,nm line formation in radiation hydrodynamics simulations of Cepheid atmospheres

Authors: V. Vasilyev, A. M. Amarsi, H. -G. Ludwig, B. Lemasle

Abstract: 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… ▽ More 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. △ Less

Submitted 17 April, 2019; v1 submitted 5 March, 2019; originally announced March 2019.

Comments: 9 pages, 10 figures; Published in Astronomy and Astrophysics

Journal ref: A&A 624, A85 (2019)

Carbon and oxygen in metal-poor halo stars

Authors: A. M. Amarsi, P. E. Nissen, M. Asplund, K. Lind, P. S. Barklem

Abstract: 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) hydrodynam… ▽ More 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$β$ 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. △ Less

Submitted 1 February, 2019; v1 submitted 11 January, 2019; originally announced January 2019.

Comments: 5 pages, 2 figures; published in A&A Letters

Journal ref: A&A 622, L4 (2019)

On the AGB stars of M 4: A robust disagreement between spectroscopic observations and theory

Authors: B. T. MacLean, S. W. Campbell, A. M. Amarsi, T. Nordlander, P. L. Cottrell, G. M. De Silva, J. Lattanzio, T. Constantino, V. DOrazi, L. Casagrande

Abstract: 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, an… ▽ More 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. △ Less

Submitted 20 August, 2018; originally announced August 2018.

Comments: 25 pages, 21 figures, 15 tables, accepted for publication in MNRAS

The GALAH Survey: Verifying abundance trends in the open cluster M67 using non-LTE spectroscopy

Authors: Xudong Gao, Karin Lind, Anish M. Amarsi, Sven Buder, Aaron Dotter, Thomas Nordlander, Martin Asplund, Joss Bland-Hawthorn, Gayandhi M. De Silva, Valentina D'Orazi, Ly Duong, Ken C. Freeman, Janez Kos, Geraint F. Lewis, Jane Lin, Sarah L. Martell, Katharine. J. Schlesinger, Sanjib Sharma, Jeffrey D. Simpson, Daniel B. Zucker, Tomaz Zwitter, Gary Da Costa, Borja Anguiano, Jonathan Horner, Elaina A. Hyde , et al. (6 additional authors not shown)

Abstract: 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 cl… ▽ More 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. △ Less

Submitted 23 August, 2018; v1 submitted 17 April, 2018; originally announced April 2018.

Comments: 14 pages, 10 figures; resubmitted to MNRAS

The GALAH Survey: Second Data Release

Authors: S. Buder, M. Asplund, L. Duong, J. Kos, K. Lind, M. K. Ness, S. Sharma, J. Bland-Hawthorn, A. R. Casey, G. M. De Silva, V. D'Orazi, K. C. Freeman, G. F. Lewis, J. Lin, S. L. Martell, K. J. Schlesinger, J. D. Simpson, D. B. Zucker, T. Zwitter, A. M. Amarsi, B. Anguiano, D. Carollo, K. Cotar, P. L. Cottrell, G. Da Costa , et al. (17 additional authors not shown)

Abstract: 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 abundance… ▽ More 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. △ Less

Submitted 17 April, 2018; originally announced April 2018.

Comments: 34+7 pages, 31 Figures, 5 tables, 1 catalog, submitted to MNRAS

Effective temperature determinations of late-type stars based on 3D non-LTE Balmer line formation

Authors: A. M. Amarsi, T. Nordlander, P. S. Barklem, M. Asplund, R. Collet, K. Lind

Abstract: 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 tim… ▽ More 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$α$, H$β$, and H$γ$, 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$γ$, while the inner wings can be weaker in 3D models, particularly for H$α$. For H$α$, 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$α$, H$β$, and H$γ$; 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σ$ uncertainties. For H$α$, 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$α$ 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. △ Less

Submitted 2 August, 2018; v1 submitted 6 April, 2018; originally announced April 2018.

Comments: 19 pages, 10 figures, abstract abridged; accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 615, A139 (2018)

Inelastic O+H collisions and the OI 777nm solar centre-to-limb variation

Authors: A. M. Amarsi, P. S. Barklem, M. Asplund, R. Collet, O. Zatsarinny

Abstract: The OI 777 nm triplet is a key diagnostic of oxygen abundances in the atmospheres of FGK-type stars; however it is sensitive to departures from local thermodynamic equilibrium (LTE). The accuracy of non-LTE line formation calculations has hitherto been limited by errors in the inelastic O+H collisional rate coefficients: several recent studies have used the so-called Drawin recipe, albeit with a c… ▽ More The OI 777 nm triplet is a key diagnostic of oxygen abundances in the atmospheres of FGK-type stars; however it is sensitive to departures from local thermodynamic equilibrium (LTE). The accuracy of non-LTE line formation calculations has hitherto been limited by errors in the inelastic O+H collisional rate coefficients: several recent studies have used the so-called Drawin recipe, albeit with a correction factor $\mathrm{S_{H}}$ that is calibrated to the solar centre-to-limb variation of the triplet. We present a new model oxygen atom that incorporates inelastic O+H collisional rate coefficients using an asymptotic two-electron model based on linear combinations of atomic orbitals, combined with a free electron model, based on the impulse approximation. Using a 3D hydrodynamic stagger model solar atmosphere and 3D non-LTE line formation calculations, we demonstrate that this physically-motivated approach is able to reproduce the solar centre-to-limb variation of the triplet to 0.02 dex, without any calibration of the inelastic collisional rate coefficients or other free parameters. We infer $\logε_{\mathrm{O}}=8.69\pm0.03$ from the triplet alone, strengthening the case for a low solar oxygen abundance. △ Less

Submitted 28 October, 2018; v1 submitted 28 March, 2018; originally announced March 2018.

Comments: 13 pages, 8 figures; published in Astronomy & Astrophysics

Journal ref: A&A 616, A89 (2018)

Non-LTE line formation of Fe in late-type stars IV: Modelling of the solar centre-to-limb variation in 3D

Authors: K. Lind, A. M. Amarsi, M. Asplund, P. S. Barklem, M. Bautista, M. Bergemann, R. Collet, D. Kiselman, J. Leenaarts, T. M. D. Pereira

Abstract: Our ability to model the shapes and strengths of iron lines in the solar spectrum is a critical test of the accuracy of the solar iron abundance, which sets the absolute zero-point of all stellar metallicities. We use an extensive 463-level Fe atom with new photoionisation cross-sections for FeI as well as quantum mechanical calculations of collisional excitation and charge transfer with neutral h… ▽ More Our ability to model the shapes and strengths of iron lines in the solar spectrum is a critical test of the accuracy of the solar iron abundance, which sets the absolute zero-point of all stellar metallicities. We use an extensive 463-level Fe atom with new photoionisation cross-sections for FeI as well as quantum mechanical calculations of collisional excitation and charge transfer with neutral hydrogen; the latter effectively remove a free parameter that has hampered all previous line formation studies of Fe in non-local thermodynamic equilibrium (NLTE). For the first time, we use realistic 3D NLTE calculations of Fe for a quantitative comparison to solar observations. We confront our theoretical line profiles with observations taken at different viewing angles across the solar disk with the Swedish 1-m Solar Telescope. We find that 3D modelling well reproduces the observed centre-to-limb behaviour of spectral lines overall, but highlight aspects that may require further work, especially cross-sections for inelastic collisions with electrons. Our inferred solar iron abundance is log(eps(Fe))=7.48+-0.04. △ Less

Submitted 11 March, 2017; originally announced March 2017.

Comments: 13 pages, 8 figures, accepted for publication in MNRAS