-
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.
-
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.
-
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.
-
Shock-heated radiation-driven outflows as a solution to the weak-wind problem of late O-type stars
Authors:
C. Lagae,
F. A. Driessen,
L. Hennicker,
N. D. Kee,
J. O. Sundqvist
Abstract:
Context. Radiation-driven mass loss is key to our understanding of massive-star evolution. However, for low-luminosity O-type stars there are big discrepancies between theoretically predicted and empirically derived mass-loss rates (called the weak-wind problem). Aims. We compute radiation-line-driven wind models of a typical weak-wind star to determine its temperature structure and the correspond…
▽ More
Context. Radiation-driven mass loss is key to our understanding of massive-star evolution. However, for low-luminosity O-type stars there are big discrepancies between theoretically predicted and empirically derived mass-loss rates (called the weak-wind problem). Aims. We compute radiation-line-driven wind models of a typical weak-wind star to determine its temperature structure and the corresponding impact on ultra-violet (UV) line formation. Methods. We carried out hydrodynamic simulations of the line-deshadowing instability (LDI) for a weak-wind star in the Galaxy. Subsequently, we used this LDI model as input in a short-characteristics radiative transfer code to compute synthetic UV line profiles. Results. We find that the line-driven weak wind is significantly shock heated to high temperatures and is unable to cool down effciently. This results in a complex temperature structure where more than half of the wind volume has temperatures significantly higher than the stellar effective temperature. Therefore, a substantial portion of the weak wind will be more ionised, resulting in a reduction of the UV line opacity and therefore in weaker line profiles for a given mass-loss rate. Quantifying this, we find that weak-wind mass-loss rates derived from unsaturated UV lines could be underestimated by a factor of between 10 and 100 if the high-temperature gas is not properly taken into account in the spectroscopic analysis. This offers a tentative basic explanation for the weak-wind problem: line-driven weak winds are not really weaker than theoretically expected, but rather a large portion of their wind volume is much hotter than the stellar effective temperature.
△ Less
Submitted 29 March, 2021;
originally announced March 2021.