We present phosphorus abundance measurements for a total of 102 giant stars, including 82 stars in 24 open clusters and 20 Cepheids, based on high-resolution near-infrared spectra obtained with GIANO-B. Evolution of phosphorus abundance, despite its astrophysical and biological significance, remains poorly understood due to a scarcity of observational data. By combining precise stellar parameters from the optical, a robust line selection and measurement method, we measure phosphorus abundances using available P I lines. Our analysis confirms a declining trend in [P/Fe] with increasing [Fe/H] around solar metallicity for clusters and Cepheids, consistent with previous studies. We also report a [P/Fe]-age relation among open clusters older than 1 Gyr, indicating a time-dependent enrichment pattern. Such a pattern can be explained by the different stellar formation history of their parental gas, with more efficient stellar formation in the gas of older clusters (thus with higher phosphorus abundances). [P/Fe] shows a flat trend among Cepheids and clusters younger than 1 Gyr (along with three Cepheids inside open clusters), possibly hinting at the phosphorus contribution from the previous-generation low-mass stars. Such a trend suggests that the young clusters share a nearly common chemical history, with a mild increase in phosphorus production by low-mass stars

Recent studies show that the Milky Way stellar halo is composed of populations of different origins, shaped by multiple accretion events. To better understand the formation of the Milky Way and other spiral galaxies, we characterize the chemical and kinematic properties of halo substructures using GALAH DR4 and Gaia data. We apply wavelet transforms in the space of sqrt(J_r) and azimuthal action (L_z) to identify kinematic overdensities. Stars in the detected structures are analyzed in elemental abundance space to determine their origin. We further assess contamination using the unsupervised machine-learning algorithm t-distributed stochastic neighbor embedding (t-SNE), performing chemical tagging with 15 elemental abundances. We recover five structures: the Galactic disk, the Splash, Gaia-Sausage-Enceladus (GSE), Thamnos1, and Thamnos2. GSE shows two peaks; one at sqrt(J_r) ~ 25 kpc km s^-1 is due to disk contamination, while the other above sqrt(J_r) ~ 40 kpc km s^-1 represents the cleanest GSE population. Thamnos exhibits three peaks linked to Thamnos1 and Thamnos2. Thamnos2 shows higher [alpha/Fe], iron-peak elements are enhanced in the Splash, and halo groups retain a stronger r-process signature. The multiply peaked structures suggest that the splashed disk extends beyond prograde orbits. The distinct chemo-dynamical properties of the halo groups support their extragalactic origin.

This White Paper presents the scientific rationale and instrument concept for HRMOS (High-Resolution Multi-Object Spectrograph), a next-generation instrument proposed for the ESO Very Large Telescope within the VLT 2030 roadmap. Current and planned facilities offer either multi-object spectroscopy or ultra-high spectral resolution, but not both. HRMOS fills this gap by combining very high spectral resolution, multi-object capability, and radial-velocity stability, enabling transformative studies in Galactic and extragalactic astrophysics. The baseline design provides a resolving power of R = 80000, radial-velocity precision of 10 m s-1 (goal: 5 m s-1), simultaneous observations of 50-60 targets, and broad optical coverage down to 385 nm. These capabilities enable precise measurements of elemental abundances, isotopic ratios, line profiles, and radial velocities for large stellar samples, including crowded fields, star clusters, the Galactic bulge, and nearby dwarf galaxies. HRMOS will address key questions on the age of the oldest stellar populations through nucleocosmochronology, the formation and survival of planetary systems, the assembly history of the Milky Way and satellites, the origin of the heaviest elements, stellar evolution, and the chemical and dynamical properties of the interstellar and circumgalactic medium. It will bridge large spectroscopic surveys and the next generation of extremely large telescopes, with strong synergies with 4MOST, Gaia, TESS, PLATO, the proposed Haydn mission, and future ELT instruments. Building on VLT/FLAMES heritage, HRMOS represents a strategic investment for European astronomy in the 2030s.

Studies of the oldest and most metal-poor stars in the Milky Way have confirmed a common abundance signature of high carbon enrichment, coupled with a subsolar abundance pattern of neutron-capture elements. The so-called CEMP-no stars have been speculated to be bona fide population II stars, potentially tracing the nucleosynthesis of the very first stars formed in the universe. However, constraining the binary nature of CEMP-no stars is crucial for understanding their abundance patterns. In previous radial-velocity monitoring of CEMP-no stars, the binary fraction has tentatively been found to vary with carbon enhancement. Here we present the results of radial-velocity monitoring of 30 CEMP-no stars over five years. Combined with literature data, this yields a total sample of 90 CEMP-no stars with constrained binary statuses, providing a larger statistical sample to investigate the CEMP-no binary fraction as a function of carbon enrichment. We find an overall binary frequency of 50^+13_-13% among high-carbon (A(C) ≥7.3) CEMP-no stars, compared to 18^+5_-4% for low-carbon (A(C) <7.3) stars, establishing for the first time a statistically significant increase in the CEMP-no binary frequency as a function of carbon at the 2σ confidence level. Of the confirmed binary systems, we derive orbital parameters for four new ones, which, combined with literature data, amount to a total of 12 CEMP-no binaries with constrained orbits. We discuss these results in the context of the progenitors of CEMP-no stars; in particular, we explore the nature of the companion in these binary systems and the possibility of mass transfer.

Context. Recent studies have revealed that the Milky Way's stellar halo is a composite of stellar populations of different origins, including multiple accretion events. To better understand how the Milky Way and other spiral galaxies were formed, it is necessary to thoroughly characterize the chemical and kinematic properties of these structures. Aims. We search for kinematic structures of the stellar halo to find any substructures within them (when indeed present) and characterize the chemodynamical properties of the identified groups with the GALAH DR4 and Gaia surveys. Methods. We applied wavelet transforms in the space defined by a square root of radial action (J(r)) and azimuthal action (L-z) to search for kinematic overdensities. Then, we selected stars in the detected structures and investigated their elemental abundances to determine their origin. Additionally, we checked for any contamination from other stellar populations within the detected groups with the unsupervised machine learning algorithm t-Distributed Stochastic Neighbor Embedding (t-SNE), for which we performed chemical tagging in a high-dimensional parameter space using 15 elemental abundances as input. Results. We recovered five kinematic structures in the action space with the wavelet transform. These groups are the Galactic disk, Splash, Gaia-Sausage-Enceladus (GSE), Thamnos1, and Thamnos2. We found that GSE has two peaks with the wavelet transform. One of these peaks is located at root J(r )similar or equal to 25 kpc km s(-1 )and is a result of contamination from disk stars. The other peak corresponds to the 'cleanest' GSE population and is located above root J(r) similar or equal to 40 kpc km s(-1) . We also detected three peaks in Thamnos. We linked two of them to Thamnos1, while the peak with the stars on the most retrograde orbits was linked to Thamnos2. The t-SNE algorithm confirmed these findings. We also analyzed individual elemental abundances of each group and found that Thamnos2 has a higher [alpha/Fe] ratio than the other groups and that iron-peak elements are more abundant in the Splash than in the halo groups, while the halo structures retain a higher r-process signature than the splashed disk. Conclusions. A multiply peaked substructure we observe in action space in GSE and Thamnos suggests that the splashed disk extends beyond the borders of prograde orbits. Each of the four halo groups studied in this paper have unique chemodynamical properties that confirm their extra-galactic origin.

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

Context . The Milky Way halo has been built-up over cosmic time through the accretion and dissolution of star clusters and dwarf galaxies as well as through their complex interactions with the Galactic disc. Traces of these accreted structures persist to the present day in the chemical and kinematic properties of stars and their orbits and allow for the disentangling of the accretion history of the Galaxy through observations of Milky Way stars. Aims . We utilised 6D phase-space information in combination with [Fe/H] measurements to facilitate a clustering analysis of stars using their kinematics and chemistry simultaneously, a technique known as chemo-kinematic tagging. We aim to associate large halo substructure groups with Milky Way halo globular clusters, stellar streams, and satellite galaxies in order to investigate the common origins of these groups of structures. Methods . We implemented t-distributed stochastic neighbour embedding (t-SNE) to perform dimensionality reduction and identify stars from clusters and streams that are co-localised in the kinematic and chemical parameter space. We used the orbital parameters E, J r , J z , L z , r apo , r peri , and eccentricity as well as [Fe/H] as input into the algorithm, and we performed a clustering analysis for a sample of 5347 stars from 229 individual Milky Way substructures. Results . Most notably, we recovered several large-scale structures that have been reported in the literature, including GSE, Thamnos, Sequoia, I’itoi, LMS-1/Wukong, Sagittarius, Kraken/Koala, the splashed disc, and a candidate structure recently found in the literature. We assigned globular cluster populations to each of these accreted structures and find that 44% of Milky Way globular clusters are consistent with having an accreted origin. In addition, we find that the chemo-dynamic properties of Omega Cen are consistent with a common accretion with the Thamnos structure. Finally, we identified many small-scale structures, including several stream-progenitor associations, and a connection between the Orphan-Chenab stream and the Grus II ultra-faint dwarf galaxy, which supports previous findings that these two objects were brought into the Galaxy in the same accretion event.

Modern spectroscopic surveys output large data volumes. Theoretical models provide a means to transform the information encoded in these data to measurements of physical stellar properties. However, in detail, the models are incomplete and simplified, and prohibit interpretation of the fine details in spectra. Instead, the available data provide an opportunity to use data-driven, differential analysis techniques, as a means towards understanding spectral signatures. We deploy such an analysis to examine core helium-fusing red clump (RC) and shell hydrogen-fusing red giant branch (RGB) stars, to uncover signatures of evolutionary state imprinted in optical stellar spectra. We exploit 786 pairs of RC and RGB stars from the GALAH survey, chosen to minimize spectral differences, with evolutionary state classifications from TESS and K2 asteroseismology. We report sub-per cent residual, systematic spectral differences between the two classes of stars, and show that these residuals are significant compared to a reference sample of RC-RC and RGB-RGB pairs selected using the same criteria. First, we report systematic differences in the Swan (C (2) ) band and CN bands caused by stellar evolution and a difference in mass, where RGB stars at similar stellar parameters have higher masses than RC stars. Secondly, we observe systematic differences in the line-width of the H (alpha) and H (beta) lines caused by a difference in microturbulence, as measured by GALAH, where we measure higher microturbulence in RC stars than RGB stars. This work demonstrates the ability of large surveys to uncover the subtle spectroscopic signatures of stellar evolution using model-free, data-driven methods.

The Ca II resonance doublet (HK) and the near-infrared triplet (CaT) are among the strongest features in stellar spectra of FGK-type stars. These spectral lines remain prominent down to extremely low metallicities and are thus useful for providing stellar parameters via ionisation balance and as radial velocity diagnostics. However, the majority of studies that model these lines in late-type stars still rely on one dimensional (1D) hydrostatic model atmospheres and the assumption of local thermodynamic equilibrium (LTE). We present 3D non-LTE radiative transfer calculations of the CaT and HK lines in an extended grid of 3D model atmospheres of metal-poor FGK-type. We investigate the impact of 3D non-LTE effects on abundances, line bisectors and radial velocities. We used a subset of 3D model atmospheres from the recently published STAGGER-grid to synthesize spectra in 3D (non-)LTE. For comparison, similar calculations were performed in 1D (non-)LTE using models from the MARCS grid. Abundance corrections for the CaT lines relative to 1D LTE range from +0.1 to -1.0 dex, with more severe corrections for strong lines in giants. With fixed line strength, the abundance corrections become more negative with increasing effective temperature and decreasing surface gravity. Radial velocity corrections relative to 1D LTE based on cross-correlation of the whole line profile range from -0.2 km/s to +1.5 km/s, with more severe corrections where the CaT lines are strongest. The corrections are even more severe if the line core alone is used to infer the radial velocity. The line strengths and shapes, and consequently the abundance and radial velocity corrections, are strongly affected by the chosen radiative transfer assumption, 1/3D (non)-LTE. We release grids of theoretical spectra that can be used to improve the accuracy of stellar spectroscopic analyses based on the Ca II triplet lines.

Large spectroscopic surveys aim to consistently compute stellar parameters of very diverse stars while minimizing systematic errors. We explore the use of stellar clusters as benchmarks to verify the precision of spectroscopic parameters in the 4. data release (DR4) of the GALAH survey. We examine 58 open and globular clusters and associations to validate measurements of temperature, gravity, chemical abundances, and stellar ages. We focus on identifying systematic errors and understanding trends between stellar parameters, particularly temperature and chemical abundances. We identify trends by stacking measurements of chemical abundances against effective temperature and modelling them with splines. We also refit spectra in three clusters with the Spectroscopy Made Easy and Korg packages to reproduce the trends in DR4 and to search for their origin by varying temperature and gravity priors, linelists, and spectral continuum. Trends are consistent between clusters of different ages and metallicities, can reach amplitudes of 0.5 dex and differ for dwarfs and giants. We use the derived trends to correct the DR4 abundances of 24 and 31 chemical elements for dwarfs and giants, and publish a detrended catalogue. While the origin of the trends could not be pinpointed, we found that: i) photometric priors affect derived abundances, ii) temperature, metallicity, and continuum levels are degenerate in spectral fitting, and it is hard to break the degeneracy even by using independent measurements, iii) the completeness of the linelist used in spectral synthesis is essential for cool stars, and iv) different spectral fitting codes produce significantly different iron abundances for stars of all temperatures. We conclude that clusters can be used to characterise the systematic errors of parameters produced in large surveys, but further research is needed to explain the origin of the trends.

Neutral sodium was the first atom detected in an exoplanetary atmosphere via transmission spectroscopy and remains the most frequently detected species due to its strong doublet in the optical. However, the center-to-limb variation (CLV) of these lines in the host star can bias the Na detection.When combined with the Rossiter-McLaughlin (RM) effect, the CLV can mimic or obscure planetary absorption features. This work investigates the impact of 3D radiation hydrodynamic stellar atmospheres and non-local thermodynamic equilibrium (NLTE) radiative transfer on the modeling of the CLV+RM effect in single-line transmission spectroscopy, to improve the characterization of exoplanet atmospheres. We produced a grid of 3D NLTE synthetic spectra for Na I for FGK stars within the following parameter space: Teff=4500-6500 K, log g=4.0-5.0 and [Fe/H]=-0.5, 0, 0.5.This grid was then interpolated to match the stellar parameters of four stars hosting giant exoplanets, to correct for the CLV+RM effect in their transmission spectra. We used ESPRESSO archival observations. Our work confirms the Na detections in three systems, namely WASP-52b, WASP-76b, and WASP-127b, improving the accuracy of the measured absorption depth. Furthermore, we find that 3D NLTE stellar models can explain the spectral features in HD 209458b without the need for any planetary absorption. In the grid of synthetic spectra, we observe that the CLV effect is stronger for stars with low Teff and high log g. However, the combined effect of CLV and RM is highly dependent on the orbital geometry of the planet-star system. With the continuous improvement of instrumentation, it is crucial to use the most accurate stellar models to correct for the CLV+RM effect in high-resolution transmission spectra, to achieve the best possible characterization of exoplanet atmospheres. We make our grid of 3D NLTE spectra for Na publicly available.

Context. Traditional one-dimensional 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 three-dimensional model atmospheres, and explore the accuracy of postprocessing 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 three-dimensional local thermodynamic equilibrium for a grid of 216 fictitious Fe I lines, spanning a wide range of oscillator strengths, excitation potentials, and central wavelengths, and eight model atmospheres that cover the stellar atmospheric parameter range (T-eff, log g, and [Fe/H]) of FGK-type stars. Using this grid, the impact of vertical and horizontal resolutions, 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, in other words, sampling every third grid cell, introduces minimal errors on the equivalent width and normalised 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 ones present in the original grid. The public release of the STAGGER-grid contains 243 models and the processed snapshots can be used to improve the accuracy of spectroscopic analyses.

Stellar mergers and accretion events have been crucial in shaping the evolution of the Milky Way (MW). These events have been dynamically identified and chemically characterised using red giants and main-sequence stars. RR Lyrae (RRL) variables can play a crucial role in tracing the early formation of the MW since they are ubiquitous, old (t$\ge$10 Gyr) low-mass stars and accurate distance indicators. We exploited Data Release 3 of the GALAH survey to identify 78 field RRLs suitable for chemical analysis. Using synthetic spectra calculations, we determined atmospheric parameters and abundances of Fe, Mg, Ca, Y, and Ba. Most of our stars exhibit halo-like chemical compositions, with an iron peak around [Fe/H]$\approx -$1.40, and enhanced Ca and Mg content. Notably, we discovered a metal-rich tail, with [Fe/H] values ranging from $-$1 to approximately solar metallicity. This sub-group includes almost ~1/4 of the sample, it is characterised by thin disc kinematics and displays sub-solar $\alpha$-element abundances, marginally consistent with the majority of the MW stars. Surprisingly, they differ distinctly from typical MW disc stars in terms of the s-process elements Y and Ba. We took advantage of similar data available in the literature and built a total sample of 535 field RRLs for which we estimated kinematical and dynamical properties. We found that metal-rich RRLs (1/3 of the sample) likely represent an old component of the MW thin disc. We also detected RRLs with retrograde orbits and provided preliminary associations with the Gaia-Sausage-Enceladus, Helmi, Sequoia, Sagittarius, and Thamnos stellar streams.

Lithium's susceptibility to burning in stellar interiors makes it an invaluable tracer for delineating the evolutionary pathways of stars, offering insights into the processes governing their development. Observationally, the complex Li production and depletion mechanisms in stars manifest themselves as Li plateaus, and as Li-enhanced and Li-depleted regions of the HR diagram. The Li-dip represents a narrow range in effective temperature close to the main-sequence turn-off, where stars have slightly super-solar masses and strongly depleted Li. To study the modification of Li through stellar evolution, we measure 3D non-local thermodynamic equilibrium (NLTE) Li abundance for 581 149 stars released in GALAH DR3. We describe a novel method that fits the observed spectra using a combination of 3D NLTE Li line profiles with blending metal line strength that are optimized on a star-by-star basis. Furthermore, realistic errors are determined by a Monte Carlo nested sampling algorithm which samples the posterior distribution of the fitted spectral parameters. The method is validated by recovering parameters from a synthetic spectrum and comparing to 26 stars in the Hypatia catalogue. We find 228 613 Li detections, and 352 536 Li upper limits. Our abundance measurements are generally lower than GALAH DR3, with a mean difference of 0.23 dex. For the first time, we trace the evolution of Li-dip stars beyond the main sequence turn-off and up the subgiant branch. This is the first 3D NLTE analysis of Li applied to a large spectroscopic survey, and opens up a new era of precision analysis of abundances for large surveys.

The Gaia -ESO Survey is an European Southern Observatory (ESO) public spectroscopic survey that targeted 10 5 stars in the Milky Way covering the major populations of the disk, bulge and halo. The observations were made using FLAMES on the VLT obtaining both UVES high ( R ~ 47 000) and GIRAFFE medium ( R ~ 20 000) resolution spectra. The analysis of the Gaia -ESO spectra was the work of multiple analysis teams (nodes) within five working groups (WG). The homogenisation of the stellar parameters within WG11 (high resolution observations of FGK stars) and the homogenisation of the stellar parameters within WG10 (medium resolution observations of FGK stars) is described here. In both cases, the homogenisation was carried out using a Bayesian Inference method developed specifically for the Gaia -ESO Survey by WG11. The method was also used for the chemical abundance homogenisation within WG11, however, the WG10 chemical abundance data set was too sparsely populated so basic corrections for each node analysis were employed for the homogenisation instead. The WG10 homogenisation primarily used the cross-match of stars with WG11 as the reference set in both the stellar parameter and chemical abundance homogenisation. In this way the WG10 homogenised results have been placed directly onto the WG11 stellar parameter and chemical abundance scales. The reference set for the metal-poor end was sparse which limited the effectiveness of the homogenisation in that regime. For WG11, the total number of stars for which stellar parameters were derived was 6 231 with typical uncertainties for T eff , log g and [Fe/H] of 32 K, 0.05 and 0.05 respectively. One or more chemical abundances out of a possible 39 elements were derived for 6 188 of the stars. For WG10, the total number of stars for which stellar parameters were derived was 76 675 with typical uncertainties for T eff , log g and [Fe/H] of 64 K, 0.15 and 0.07 respectively. One or more chemical abundances out of a possible 30 elements were derived for 64177 of the stars.

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

ABSTRACT WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, saw first light in late 2022. WEAVE comprises a new 2-deg field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable ‘mini’ integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366–959 nm at R ∼ 5000, or two shorter ranges at $R\sim 20\, 000$. After summarizing the design and implementation of WEAVE and its data systems, we present the organization, science drivers, and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy’s origins by completing Gaia’s phase-space information, providing metallicities to its limiting magnitude for ∼3 million stars and detailed abundances for ∼1.5 million brighter field and open-cluster stars; (ii) survey ∼0.4 million Galactic-plane OBA stars, young stellar objects, and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey ∼400 neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionized gas in z < 0.5 cluster galaxies; (vi) survey stellar populations and kinematics in ${\sim} 25\, 000$ field galaxies at 0.3 ≲ z ≲ 0.7; (vii) study the cosmic evolution of accretion and star formation using >1 million spectra of LOFAR-selected radio sources; and (viii) trace structures using intergalactic/circumgalactic gas at z > 2. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.

Context. Neutral sodium was the first atom that was detected in an exoplanetary atmosphere using the transmission spectroscopy technique. To date, it remains the most successfully detected species due to its strong doublet in the optical at 5890 & Aring; and 5896 & Aring;. However, the center-to-limb variation (CLV) of these lines in the host star can bias the Na I detection. When combined with the Rossiter-McLaughlin (RM) effect, the CLV can mimic or obscure a planetary absorption feature if it is not properly accounted for. Aims. This work aims to investigate the impact of three-dimensional (3D) radiation hydrodynamic stellar atmospheres and non-local thermodynamic equilibrium (NLTE) radiative transfer on the modeling of the CLV+RM effect in single-line transmission spectroscopy to improve the detection and characterization of exoplanet atmospheres. Methods. We produced a grid of 3D NLTE synthetic spectra for Na I for FGK-type dwarfs within the following parameter space: Teff = 4500-6500 K, log g = 4.0-5.0, and [Fe/H] = [-0.5, 0, 0.5]. This grid was then interpolated to match the stellar parameters of four stars hosting well-known giant exoplanets, generating stellar spectra to correct for the CLV+RM effect in their transmission spectra. We used archival observations taken with the high-resolution ESPRESSO spectrograph. Results. Our work confirms the Na I detections in three systems, namely WASP-52b, WASP-76b, and WASP-127b, also improving the accuracy of the measured absorption depth. Furthermore, we find that 3D NLTE stellar models can explain the spectral features in the transmission spectra of HD 209458b without the need for any planetary absorption. In the grid of stellar synthetic spectra, we observe that the CLV effect is stronger for stars with low Teff and high log g. However, the combined effect of CLV and RM is highly dependent on the orbital geometry of the planet-star system. Conclusions. With the continuous improvement of instrumentation, it is crucial to use the most accurate stellar models available to correct for the CLV+RM effect in high-resolution transmission spectra to achieve the best possible characterization of exoplanet atmospheres. This will be fundamental in preparation for instruments such as ANDES at the Extremely Large Telescope to fully exploit its capabilities in the near future. We make our grid of 3D NLTE synthetic spectra for Na I publicly available.

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

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