Search

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 ≤ Teff∕K ≤ 8000, − 0.5 ≤log g∕cm s−2 ≤ 5.5, and − 5 ≤ [Fe/H] ≤ 1. 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 50 126 stars from GALAH DR3. We found that relaxing LTE can change the abundances by between − 0.7 dex and + 0.2 dex for different lines and stars. Taking departures from LTE into account can reduce the dispersion in the [A/Fe] versus [Fe/H] plane by up to 0.1 dex, and it can remove spurious differences between the dwarfs and giants by up to 0.2 dex. The resulting abundance slopes can thus be qualitatively different in non-LTE, possibly with important implications for the chemical evolution of our Galaxy. The grids of departure coefficients are publicly available and can be implemented into LTE pipelines to make the most of observational data sets from large spectroscopic surveys. [ABSTRACT FROM AUTHOR]

We present a line-by-line differential analysis of a sample of 16 planet-hosting stars and 68 comparison stars using high-resolution, high signal-to-noise ratio spectra gathered using Keck. We obtained accurate stellar parameters and high-precision relative chemical abundances with average uncertainties in T eff, log  g , [Fe/H], and [X/H] of 15 K, 0.034 cm s−2 , 0.012 dex, and 0.025 dex, respectively. For each planet host, we identify a set of comparison stars and examine the abundance differences (corrected for Galactic chemical evolution effect) as a function of the dust condensation temperature, T cond, of the individual elements. While we confirm that the Sun exhibits a negative trend between abundance and T cond, we also confirm that the remaining planet hosts exhibit a variety of abundance– T cond trends with no clear dependence upon age, metallicity, or T eff. The diversity in the chemical compositions of planet-hosting stars relative to their comparison stars could reflect the range of possible planet-induced effects present in these planet hosts, from the sequestration of rocky material (refractory poor) to the possible ingestion of planets (refractory rich). Other possible explanations include differences in the time-scale, efficiency and degree of planet formation, or inhomogeneous chemical evolution. Although we do not find an unambiguous chemical signature of planet formation among our sample, the high-precision chemical abundances of the host stars are essential for constraining the composition and structure of their exoplanets. [ABSTRACT FROM AUTHOR]

Context. The benchmark open cluster M 67 is known to have solar metallicity and an age similar to that of the Sun. It thus provides us with a great opportunity to study the properties of solar twins, as well as the evolution of Sun-like stars. Aims. Previous spectroscopic studies of M 67 reported possible subtle changes in stellar surface abundances throughout the stellar evolutionary phase, namely the effect of atomic diffusion. In this study we attempt to confirm and quantify more precisely the effect of atomic diffusion, and to explore the level of chemical (in)homogeneity in M 67. Methods. We presented a strictly line-by-line differential chemical abundance analysis of two groups of stars in M 67: three turn-off stars and three subgiants. Stellar atmospheric parameters and elemental abundances were obtained with very high precision using the Keck/HIRES spectra. Results. The subgiants in our sample show negligible abundance variations (≤0.02 dex), which implies that M 67 was born chemically homogeneous. We note that there is a significant abundance difference (~0.1–0.2 dex) between subgiants and turn-off stars, which can be interpreted as the signature of atomic diffusion. Qualitatively stellar models with diffusion agree with the observed abundance results. Some turn-off stars do not follow the general pattern, which suggests that in some cases diffusion can be inhibited, or they might have undergone some sort of mixing event related to planets. Conclusions. Our results pose additional challenges for chemical tagging when using turn-off stars. In particular, the effects of atomic diffusion, which could be as large as 0.1–0.2 dex, must be taken into account in order for chemical tagging to be successfully applied. [ABSTRACT FROM AUTHOR]

Hydrogen Balmer lines are commonly used as spectroscopic effective temperature diagnostics of late-type stars. However, reliable inferences require accurate model spectra, and 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 (Teff ≈ 6500 K) the inner wings tend to be weaker in non-LTE models, while at lower effective temperatures (Teff ≈ 4500 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 50 K 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 100 K, in terms of inferred effective temperatures. For metal-poor turn-off stars, 1D LTE models of Hα can underestimate effective temperatures by around 150 K. Our 3D non-LTE model spectra are publicly available, and can be used for more reliable spectroscopic effective temperature determinations. [ABSTRACT FROM AUTHOR]

We present an abundance analysis of the low-metallicity benchmark red giant star HD122563 based on realistic, state-of-the-art, high-resolution, three-dimensional (3D) model stellar atmospheres including non-grey radiative transfer through opacity binning with 4, 12, and 48 bins. The 48-bin 3D simulation reaches temperatures lower by ~300-500K than the corresponding 1D model in the upper atmosphere. Small variations in the opacity binning, adopted line opacities, or chemical mixture can cool the photospheric layers by a further ~100-300K and alter the effective temperature by ~100 K. A 3D local thermodynamic equilibrium (LTE) spectroscopic analysis of Fe I and Fe II lines gives discrepant results in terms of derived Fe abundance, which we ascribe to non-LTE effects and systematic errors on the stellar parameters. We also determine C, N, and O abundances by simultaneously fitting CH, OH, NH, and CN molecular bands and lines in the ultraviolet, visible, and infrared. We find a small positive 3D-1D abundance correction for carbon (+0.03 dex) and negative ones for nitrogen (-0.07 dex) and oxygen (-0.34 dex). From the analysis of the [OI] line at 6300.3 Å, we derive a significantly higher oxygen abundance than from molecular lines (+0.46 dex in 3D and +0.15 dex in 1D). We rule out important OH photodissociation effects as possible explanation for the discrepancy and note that lowering the surface gravity would reduce the oxygen abundance difference between molecular and atomic indicators. [ABSTRACT FROM AUTHOR]

Context. The surface structures and dynamics of cool stars are characterised by the presence of convective motions and turbulent flows which shape the emergent spectrum. Aims. We used realistic three-dimensional (3D) radiative hydrodynamical simulations from the STAGGER-grid to calculate synthetic spectra with the radiative transfer code OPTIM3D for stars with different stellar parameters to predict photometric colours and convective velocity shifts. Methods. We calculated spectra from 1000 to 200 000Å with a constant resolving power of λ/Ξλ = 20 000 and from 8470 and 8710Å (Gaia Radial Velocity Spectrometer - RVS - spectral range) with a constant resolving power of λ/Ξλ = 300 000. Results. We used synthetic spectra to compute theoretical colours in the Johnson-Cousins UBV(RI)C, SDSS, 2MASS, Gaia, SkyMapper, Strömgren systems, and HST-WFC3. Our synthetic magnitudes are compared with those obtained using 1D hydrostatic models. We showed that 1D versus 3D differences are limited to a small percent except for the narrow filters that span the optical and UV region of the spectrum. In addition, we derived the effect of the convective velocity fields on selected Fe I lines. We found the overall convective shift for 3D simulations with respect to the reference 1D hydrostatic models, revealing line shifts of between -0.235 and +0.361 km s-1.We showed a net correlation of the convective shifts with the effective temperature: lower effective temperatures denote redshifts and higher effective temperatures denote blueshifts. We conclude that the extraction of accurate radial velocities from RVS spectra need an appropriate wavelength correction from convection shifts. Conclusions. The use of realistic 3D hydrodynamical stellar atmosphere simulations has a small but significant impact on the predicted photometry compared with classical 1D hydrostatic models for late-type stars.We make all the spectra publicly available for the community through the POLLUX database. [ABSTRACT FROM AUTHOR]

Context. Red giant branch (RGB) stars are very bright objects in galaxies and are often used as standard candles. Interferometry is the ideal tool to characterize the dynamics and morphology of their atmospheres. Aims. We aim at precisely characterising the surface dynamics of a sample of RGB stars. Methods. We obtained interferometric observations for three RGB stars with the MIRC instrument mounted at the CHARA interferometer. We looked for asymmetries on the stellar surfaces using limb-darkening models. Results. We measured the apparent diameters of HD197989 (ε Cyg) = 4.61 ± 0.02 mas, HD189276 (HR 7633) = 2.95 ± 0.01 mas, and HD161096 (β Oph) = 4.43 ± 0.01 mas. We detected departures from the centrosymmetric case for all three stars with the tendency of a greater effect for lower logg of the sample. We explored the causes of this signal and conclude that a possible explanation to the interferometric signal is the convection-related and/or the magnetic-related surface activity. However, it is necessary to monitor these stars with new observations, possibly coupled with spectroscopy, in order to firmly establish the cause. [ABSTRACT FROM AUTHOR]

As one of the most important elements in astronomy, iron abundance determinations need to be as accurate as possible. We investigate the accuracy of spectroscopic iron abundance analyses using archetypal metal-poor stars. We perform detailed 3D non-LTE radiative transfer calculations based on 3D hydrodynamic STAGGER model atmospheres, and employ a new model atom that includes new quantum-mechanical neutral hydrogen collisional rate coefficients. With the exception of the red giant HD122563, we find that the 3D non-LTE models achieve Fe I/Fe II excitation and ionization balance as well as not having any trends with equivalent width to within modelling uncertainties of 0.05 dex, all without having to invoke any microturbulent broadening; for HD122563 we predict that the current best parallax-based surface gravity is overestimated by 0.5 dex. Using a 3D non-LTE analysis, we infer iron abundances from the 3D model atmospheres that are roughly 0.1 dex higher than corresponding abundances from 1D MARCS model atmospheres; these differences go in the same direction as the non-LTE effects themselves. We make available grids of departure coefficients, equivalent widths and abundance corrections, calculated on 1D MARCS model atmospheres and horizontally and temporally averaged 3D STAGGER model atmospheres. [ABSTRACT FROM AUTHOR]

Stars in open clusters are expected to share an identical abundance pattern. Establishing the level of chemical homogeneity in a given open cluster deserves further study as it is the basis of the concept of chemical tagging to unravel the history of the Milky Way. M67 is particularly interesting given its solar metallicity and age as well as being a dense cluster environment. We conducted a strictly line-by-line differential chemical abundance analysis of two solar twins in M67: M67-1194 and M67-1315. Stellar atmospheric parameters and elemental abundances were obtained with high precision using Keck/High Resolution Echelle Spectrometer spectra. M67-1194 is essentially identical to the Sun in terms of its stellar parameters. M67-1315 is warmer than M67-1194 by ≈150 K as well as slightly more metal-poor than M67-1194 by ≈0.05 dex. M67-1194 is also found to have identical chemical composition to the Sun, confirming its solar-twin nature. The abundance ratios [X/Fe] of M67-1315 are similar to the solar abundances for elements with atomic number Z ≤ 30, while most neutron-capture elements are enriched by ≈0.05 dex, which might be attributed to enrichment from a mixture of asymptotic giant branch ejecta and r-process material. The distinct chemical abundances for the neutron-capture elements in M67-1315 and the lower metallicity of this star compared to M67-1194, indicate that the stars in M67 are likely not chemically homogeneous. This poses a challenge for the concept of chemical tagging since it is based on the assumption of stars forming in the same star-forming aggregate. [ABSTRACT FROM AUTHOR]

We present a high-precision differential abundance analysis of 16 solar-type stars in the Hyades open cluster based on high-resolution, high signal-to-noise ratio (S/N ≈ 350-400) spectra obtained from the McDonald 2.7-m telescope. We derived stellar parameters and differential chemical abundances for 19 elements (C, O, Na, Mg, Al, Si, S, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn and Ba) with uncertainties as low as ~0.01-0.02 dex. Our main results include: (1) there is no clear chemical signature of planet formation detected among the sample stars, i.e. no correlations in abundances versus condensation temperature; (2) the observed abundance dispersions are a factor of ≈1.5-2 larger than the average measurement errors for most elements; (3) there are positive correlations, of high statistical significance, between the abundances of at least 90 per cent of pairs of elements. We demonstrate that none of these findings can be explained by errors due to the stellar parameters. Our results reveal that the Hyades is chemically inhomogeneous at the 0.02 dex level. Possible explanations for the abundance variations include (1) inhomogeneous chemical evolution in the proto-cluster environment, (2) supernova ejection in the proto-cluster cloud and (3) pollution of metal-poor gas before complete mixing of the proto-cluster cloud. Our results provide significant new constraints on the chemical composition of open clusters and a challenge to the current view of Galactic archaeology. [ABSTRACT FROM AUTHOR]

The O I 777 nm lines are among the most commonly used diagnostics for the oxygen abundances in the atmospheres of FGK-type stars. However, they form in conditions that are far from local thermodynamic equilibrium (LTE).We explore the departures from LTE of atomic oxygen, and their impact on O I lines, across the STAGGER-grid of three-dimensional hydrodynamic model atmospheres. For the O I 777 nm triplet, we find significant departures from LTE. These departures are larger in stars with larger effective temperatures, smaller surface gravities, and larger oxygen abundances. We present grids of predicted 3D non-LTE based equivalent widths for the O I 616 nm, [O I] 630 nm, [O I] 636 nm, and O I 777 nm lines, as well as abundance corrections to 1D LTE based results. [ABSTRACT FROM AUTHOR]

We revisit the Galactic chemical evolution of oxygen, addressing the systematic errors inherent in classical determinations of the oxygen abundance that arise from the use of one-dimensional (1D) hydrostatic model atmospheres and from the assumption of local thermodynamic equilibrium (LTE). We perform detailed 3D non-LTE radiative-transfer calculations for atomic oxygen lines across a grid of 3D hydrodynamic STAGGER model atmospheres for dwarfs and subgiants. We apply our grid of predicted line strengths of the [O I] 630 nm and O I 777 nm lines using accurate stellar parameters from the literature. We infer a steep decay in [O/Fe] for [Fe/H] ≳ -1.0, a plateau [O/Fe] ≈ 0.5 down to [Fe/H] ≈ -2.5, and an increasing trend for [Fe/H] ≲ -2.5. Our 3D non-LTE calculations yield overall concordant results from the two oxygen abundance diagnostics. [ABSTRACT FROM AUTHOR]

Aims. We compute the emergent stellar spectra from the UV to far infrared for different viewing angles using realistic 3D model atmospheres for a large range in stellar parameters to predict the stellar limb darkening. Methods. We have computed full 3D LTE synthetic spectra based on 3D radiative hydrodynamic atmosphere models from the Stagger-grid in the ranges: Teff from 4000 to 7000 K, log g from 1.5 to 5.0, and [Fe/H], from -4.0 to +0.5. From the resulting intensities, we derived coefficients for the standard limb darkening laws considering a number of often-used photometric filters. Furthermore, we calculated theoretical transit light curves, in order to quantify the differences between predictions by the widely used 1D model atmosphere and our 3D models. Results. The 3D models are often found to predict steeper darkening towards the limb compared to the 1D models, mainly due to the temperature stratifications and temperature gradients being different in the 3D models compared to those predicted with 1D models based on the mixing length theory description of convective energy transport. The resulting differences in the transit light curves are rather small; however, these can be significant for high-precision observations of extrasolar transits, and are able to lower the residuals from the fits with 1D limb darkening profiles. Conclusions. We advocate the use of the new limb darkening coefficients provided for the standard four-parameter non-linear power law, which can fit the limb darkening more accurately than other choices. [ABSTRACT FROM AUTHOR]

Aims. We study the implications of averaging methods with different reference depth scales for 3D hydrodynamical model atmospheres computed with the Stagger-code. The temporally and spatially averaged (hereafter denoted as <3D>) models are explored in the light of local thermodynamic equilibrium (LTE) spectral line formation by comparing spectrum calculations using full 3D atmosphere structures with those from <3D> averages. Methods. We explored methods for computing mean <3D> stratifications from the Stagger-grid time-dependent 3D radiative hydrodynamical atmosphere models by considering four different reference depth scales (geometrical depth, column-mass density, and two optical depth scales). Furthermore, we investigated the influence of alternative averages (logarithmic, enforced hydrostatic equilibrium, flux-weighted temperatures). For the line formation we computed curves of growth for Fe i and Fe ii lines in LTE. Results. The resulting {3D} stratifications for the four reference depth scales can be very different. We typically find that in the upper atmosphere and in the superadiabatic region just below the optical surface, where the temperature and density fluctuations are highest, the differences become considerable and increase for higher Teff, lower log g, and lower [Fe/H]. The differential comparison of spectral line formation shows distinctive differences depending on which <3D> model is applied. The averages over layers of constant column-mass density yield the best mean <3D> representation of the full 3D models for LTE line formation, while the averages on layers at constant geometrical height are the least appropriate. Unexpectedly, the usually preferred averages over layers of constant optical depth are prone to increasing interference by reversed granulation towards higher effective temperature, in particular at low metallicity. [ABSTRACT FROM AUTHOR]

We compare limb darkening laws derived from 3D hydrodynamical model atmospheres and 1D hydrostatic MARCS models for the host stars of two well-studied transiting exoplanet systems, the late-type dwarfs HD 209458 and HD 189733. The surface brightness distribution of the stellar disks is calculated for a wide spectral range using 3D LTE spectrum formation and opacity sampling". We test our theoretical predictions using least-squares fits of model light curves to wavelength-integrated primary eclipses that were observed with the Hubble Space Telescope (HST). The limb darkening law derived from the 3D model of HD 209458 in the spectral region between 2900 Å and 5700 Å produces significantly better fits to the HST data, removing systematic residuals that were previously observed for model light curves based on 1D limb darkening predictions. This difference arises mainly from the shallower mean temperature structure of the 3D model, which is a consequence of the explicit simulation of stellar surface granulation where 1D models need to rely on simplified recipes. In the case of HD 189733, the model atmospheres produce practically equivalent limb darkening curves between 2900 Å and 5700 Å, partly due to obstruction by spectral lines, and the data are not sufficient to distinguish between the light curves. We also analyze HST observations between 5350 Å and 10 500 Å for this star; the 3D model leads to a better fit compared to 1D limb darkening predictions. The significant improvement of fit quality for the HD 209458 system demonstrates the higher degree of realism of 3D hydrodynamical models and the importance of surface granulation for the formation of the atmospheric radiation field of late-type stars. This result agrees well with recent investigations of limb darkening in the solar continuum and other observational tests of the 3D models. The case of HD 189733 is no contradiction as the model light curves are less sensitive to the temperature stratification of the stellar atmosphere and the observed data in the 2900-5700 Å region are not sufficient to distinguish more clearly between the 3D and 1D limb darkening predictions. [ABSTRACT FROM AUTHOR]

We present a high-precision, differential elemental abundance analysis of the HAT-P-1 stellar binary based on high-resolution, high signal-to-noise ratio Keck/HIRES (High Resolution Echelle Spectrometer) spectra. The secondary star in this double system is known to host a transiting giant planet while no planets have yet been detected around the primary star. The derived metallicities ([Fe/H]) of the primary and secondary stars are identical within the errors: 0.146 ± 0.014 dex (σ = 0.033 dex) and 0.155 ± 0.007 dex (σ = 0.023 dex), respectively. Extremely precise differential abundance ratios of 23 elements have been measured (mean error of σ([X/Fe]) = 0.013 dex) and are found to be indistinguishable between the two stars: Δ[X/Fe] (secondary − primary) = +0.001 ± 0.006 dex (σ = 0.008 dex). The striking similarity in the chemical composition of the two stellar components in HAT-P-1 is contrary to the possible 0.04 dex level difference seen in 16 Cyg A+B, which also hosts a giant planet, at least three times more massive than the one around HAT-P-1 secondary star. We conclude that the presence of giant planets does not necessarily imply differences in the chemical compositions of the host stars. The elemental abundances of each star in HAT-P-1 relative to the Sun show an identical, positive correlation with the condensation temperature of the elements; their abundance patterns are thus very similar to those observed in the majority of solar twins. In view of the Meléndez et al. interpretation of the peculiar solar abundance pattern, we conclude that HAT-P-1 experienced less efficient formation of terrestrial planets than the Sun. This is in line with the expectation that the presence of close-in giant planets preventing the formation or survival of terrestrial planets. [ABSTRACT FROM AUTHOR]