Source author record

Martin Asplund

Martin Asplund appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

75works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

75 published item(s)

preprint2022arXiv

The SAPP pipeline for the determination of stellar abundances and atmospheric parameters of stars in the core program of the PLATO mission

We introduce the SAPP (Stellar Abundances and atmospheric Parameters Pipeline), the prototype of the code that will be used to determine parameters of stars observed within the core program of the PLATO space mission. The pipeline is based on the Bayesian inference and provides effective temperature, surface gravity, metallicity, chemical abundances, and luminosity. The code in its more general version can have a much wider range of applications. It can also provide masses, ages, and radii of stars and can be used for stars of stellar types not targeted by the PLATO core program, such as red giants. We validate the code on a set of 27 benchmark stars that includes 19 FGK-type dwarfs, 6 GK-type sub-giants, and 2 red giants. Our results suggest that combining various observables is the optimal approach, as it allows to break degeneracies between different parameters and yields more accurate values of stellar parameters and more realistic uncertainties. For the PLATO core sample, we obtain a typical uncertainty of 27 ($\rm{syst.}$) $\pm$ 37 ($\rm{stat.}$) K for T$_{\rm{eff}}$, 0.00 $\pm$ 0.01 dex for log$g$, 0.02 $\pm$ 0.02 dex for metallicity [Fe/H], -0.01 $\pm$ 0.03 R$_\odot$ for radii, -0.01 $\pm$ 0.05 M$_\odot$ for stellar masses, and -0.14 $\pm$ 0.63 Gyrs for ages. We also show that the best results are obtained by combining the $ν_{max}$ scaling relation and stellar spectra. This resolves the notorious problem of degeneracies, which is particularly important for F-type stars.

preprint2021arXiv

Explosive nucleosynthesis of a metal-deficient star as the source of a distinct odd-even effect in the solar twin HIP 11915

The abundance patterns observed in the Sun and in metal-poor stars show a clear odd-even effect. An important question is whether the odd-even effect in solar-metallicity stars is similar to the Sun, or if there are variations that can tell us about different chemical enrichment histories. In this work, we report for the first time observational evidence of a differential odd-even effect in the solar twin HIP 11915, relative to the solar odd-even abundance pattern. The spectra of this star were obtained with high resolving power (140 000) and signal-to-noise ratio ($\sim$420) using the ESPRESSO spectrograph and the VLT telescope. Thanks to the high spectral quality, we obtained extremely precise stellar parameters ($σ(T_{\rm eff})$ = 2 K, $σ(\rm{[Fe/H]})$ = 0.003 dex, and $σ(\log g)$ = 0.008 dex). We determine the chemical abundance of 20 elements ($Z\leq39$) with high precision ($\sim$0.01 dex), which shows a strong pattern of the odd-even effect even after performing Galactic Chemical Evolution corrections. The odd-even effect is reasonably well-reproduced by a core-collapse supernova of 13 $\rm{M_{\odot}}$ and metallicity Z = 0.001 diluted into a metal-poor gas of 1 $\rm{M_{\odot}}$. Our results indicate that HIP 11915 has an odd-even effect slightly different than the Sun, thus confirming a different supernova enrichment history.

preprint2021arXiv

The COMBS Survey -- II. Distinguishing the Metal-Poor Bulge from the Halo Interlopers

The metal-poor stars in the bulge are important relics of the Milky Way's formation history, as simulations predict that they are some of the oldest stars in the Galaxy. In order to determine if they are truly ancient stars, we must understand the origins of this population. Currently, it is unclear if the metal-poor stars in the bulge ([Fe/H] < -1 dex) are merely halo interlopers, a unique accreted population, part of the boxy/peanut-shaped (B/P) bulge or a classical bulge population. In this work, we use spectra from the VLT/FLAMES spectrograph to obtain metallicity estimates using the Ca-II triplet (CaT) of 473 bulge stars (187 of which have [Fe/H]<-1 dex), targeted using SkyMapper photometry. We also use Gaia DR2 parallaxes and proper motions to infer the Galactic positions and velocities along with orbital properties for 523 bulge stars. We employ a probabilistic orbit analysis and find that about half of our sample has a > 50\% probability of being bound to the bulge, and half are halo interlopers. We also see that the occurrence rate of halo interlopers increases steadily with decreasing metallicity across the full range of our sample (-3 < [Fe/H] < 0.5). Our examination of the kinematics of the confined compared to the unbound stars indicates the metal-poor bulge ([Fe/H] < -1 dex) comprises at least two populations; those confined to the boxy/peanut bulge and halo stars passing through the inner galaxy. We conclude that an orbital analysis approach, as we have employed, is important to uncover and understand the composite nature of the metal-poor stars in the inner region.

preprint2021arXiv

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

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.

preprint2020arXiv

Atomic data for the Gaia-ESO Survey

We describe the atomic and molecular data that were used for the abundance analyses of FGK-type stars carried out within the Gaia-ESO Survey. We present an unprecedented effort to create a homogeneous line list, which was used by several abundance analysis groups to calculate synthetic spectra and equivalent widths. The atomic data are accompanied by quality indicators and detailed references to the sources. The atomic and molecular data are made publicly available in electronic form. In general experimental transition probabilities were preferred but theoretical values were also used. Astrophysical gf-values were avoided due to the model-dependence of such a procedure. For elements whose lines are significantly affected by hyperfine structure or isotopic splitting a concerted effort has been made to collate the necessary data for the individual line components. We also performed a detailed investigation of available data for line broadening due to collisions with neutral hydrogen atoms. Synthetic spectra calculated for the Sun and Arcturus were used to assess the blending properties of the lines. Among a subset of over 1300 lines of 35 elements in the wavelength ranges from 475 nm to 685 nm and from 850 nm to 895 nm we identified about 200 lines of 24 species which have accurate gf-values and are free of blends in the spectra of the Sun and Arcturus. For the broadening due to collisions with neutral hydrogen we recommend data based on Anstee-Barklem-O'Mara theory, where available, and to avoid lines of neutral species otherwise. Theoretical broadening data by R.L. Kurucz should be used for Sc II, Ti II, and Y II lines. For ionised rare-earth species the Unsöld approximation with an enhancement factor of 1.5 for the line width can be used. Desirable improvements in atomic data were identified for a number of species, including Al I, S I, Cr II, Na I, Si I, Ca II, and Ni I.

preprint2020arXiv

Convective Excitation and Damping of Solar-like Oscillations

The last decade has seen a rapid development in asteroseismology thanks to the CoRoT and Kepler missions. With more detailed asteroseismic observations available, it is becoming possible to infer exactly how oscillations are driven and dissipated in solar-type stars. We have carried out three-dimensional (3D) stellar atmosphere simulations together with one-dimensional (1D) stellar structural models of key benchmark turn-off and subgiant stars to study this problem from a theoretical perspective. Mode excitation and damping rates are extracted from 3D and 1D stellar models based on analytical expressions. Mode velocity amplitudes are determined by the balance between stochastic excitation and linear damping, which then allows the estimation of the frequency of maximum oscillation power, $ν_{\max}$, for the first time based on ab initio and parameter-free modelling. We have made detailed comparisons between our numerical results and observational data and achieved very encouraging agreement for all of our target stars. This opens the exciting prospect of using such realistic 3D hydrodynamical stellar models to predict solar-like oscillations across the HR-diagram, thereby enabling accurate estimates of stellar properties such as mass, radius and age.

preprint2020arXiv

K2-HERMES II. Planet-candidate properties from K2 Campaigns 1-13

Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present complete results for planet-candidate hosts from the K2-HERMES survey, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra for more than 30,000 K2 stars. We present complete host-star parameters and planet-candidate radii for 224 K2 candidate planets from C1-C13. Our results cast severe doubt on 30 K2 candidates, as we derive unphysically large radii, larger than 2 $R_{Jup}$. This work highlights the importance of obtaining accurate, precise, and self-consistent stellar parameters for ongoing large planet search programs - something that will only become more important in the coming years, as TESS begins to deliver its own harvest of exoplanets.

preprint2020arXiv

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

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.

preprint2019arXiv

Discovery of a 21 Myr old stellar population in the Orion complex

The Orion complex is arguably the most studied star-forming region in the Galaxy. While stars are still being born in the Orion nebula, the oldest part was believed to be no more than 13 Myr old. In order to study the full hierarchy of star formation across the Orion complex, we perform a clustering analysis of the Ori OB1a region using new stellar surveys and derive robust ages for each identified stellar aggregate. We use Gaia DR2 parameters supplemented with radial velocities from the GALAH and APOGEE surveys to perform clustering of the Ori OB1a association. Five overdensities are resolved in a six-dimensional parameter space (positions, distance, proper motions, and radial velocity). Most correspond to previously known structures (ASCC 16, 25 Orionis, ASCC 20, ASCC 21). We use Gaia DR2, Pan-STARRS1 and 2MASS photometry to fit isochrones to the colour-magnitude diagrams of the identified clusters. The ages of the clusters can thus be measured with ~10% precision. While four of the clusters have ages between 11 and 13 Myr, the ASCC 20 cluster stands out at an age of 21 $\pm$ 3 Myr. This is significantly greater than the age of any previously known component of the Orion complex. To some degree, all clusters overlap in at least one of the six phase-space dimensions. We argue that the formation history of the Orion complex, and its relation to the Gould belt, must be reconsidered. A significant challenge in reconstructing the history of the Ori OB1a association is to understand the impact of the newly discovered 21 Myr old population on the younger parts of the complex, including their formation.

preprint2019arXiv

The GALAH Survey: Chemodynamics of the Solar Neighbourhood

We present the chemodynamic structure of the solar neighbourhood using 62 814 stars within a 500 pc sphere of the Sun observed by GALAH and with astrometric parameters from Gaia DR2. We measure the velocity dispersion for all three components (vertical, radial, and tangential) and find that it varies smoothly with [Fe/H] and [$α$/Fe] for each component. The vertical component is especially clean, with $σ_{v_z}$ increasing from a low of $8$ km s$^{-1}$ at solar-[$α$/Fe] and [Fe/H] to a high of more than 50 km s$^{-1}$ for more metal-poor and [$α$/Fe] enhanced populations. We find no evidence of a large decrease in the velocity dispersion of the higher-[$α$/Fe] populations as claimed in analysis prior to Gaia DR2, although the trend of increasing velocity dispersion with [$α$/Fe] for the same metallicity does significantly flatten at high-[$α$/Fe]. The eccentricity distribution for local stars varies most strongly as a function of [$α$/Fe], where stars with [$α$/Fe]$<0.1$ dex having generally circular orbits ($e<0.15$), while the median eccentricity increases rapidly for more [$α$/Fe]-enhanced stellar populations up to $e\sim0.35$. These [$α$/Fe]-enhanced populations have guiding radii consistent with origins in the inner Galaxy. Of the stars with metallicities much higher than the local ISM ([Fe/H]>0.1 dex), we find that more than 70\% have $e<0.2$ and are likely observed in the solar neighbourhood through churning/migration rather than blurring effects, as the epicyclic motion for these stars is not large enough to reach the radii at which they were likely born based on their metallicity.

preprint2019arXiv

The GALAH Survey: Temporal Chemical Enrichment of the Galactic Disk

We present isochrone ages and initial bulk metallicities ($\rm [Fe/H]_{bulk}$, by accounting for diffusion) of 163,722 stars from the GALAH Data Release 2, mainly composed of main sequence turn-off stars and subgiants ($\rm 7000 K>T_{eff}>4000 K$ and $\rm log g>3$ dex). The local age-metallicity relationship (AMR) is nearly flat but with significant scatter at all ages; the scatter is even higher when considering the observed surface abundances. After correcting for selection effects, the AMR appear to have intrinsic structures indicative of two star formation events, which we speculate are connected to the thin and thick disks in the solar neighborhood. We also present abundance ratio trends for 16 elements as a function of age, across different $\rm [Fe/H]_{bulk}$ bins. In general, we find the trends in terms of [X/Fe] vs age from our far larger sample to be compatible with studies based on small ($\sim$ 100 stars) samples of solar twins but we now extend it to both sub- and super-solar metallicities. The $α$-elements show differing behaviour: the hydrostatic $α$-elements O and Mg show a steady decline with time for all metallicities while the explosive $α$-elements Si, Ca and Ti are nearly constant during the thin disk epoch (ages $\lessapprox $ 12 Gyr). The s-process elements Y and Ba show increasing [X/Fe] with time while the r-process element Eu have the opposite trend, thus favouring a primary production from sources with a short time-delay such as core-collapse supernovae over long-delay events such as neutron star mergers.

preprint2016arXiv

3D NLTE analysis of the most iron-deficient star, SMSS0313-6708

Models of star formation in the early universe depend on the details of accretion, fragmentation and radiative feedback. Different simulations predict different initial mass functions of the first stars, ranging from predominantly low mass (0.1-10 Msol), to massive (10-100 Msol), or even supermassive (100-1000 Msol). The mass distribution of the first stars should lead to unique chemical imprints on the low-mass second and later generation metal-poor stars still in existence. The chemical composition of SMSS0313-6708, which has the lowest abundances of Ca and Fe of any star known, indicates it was enriched by a single massive supernova. However, even weak spectral lines may be affected by strong 3D and NLTE effects in metal-poor stars. If these effects are ignored or treated incorrectly, errors in the inferred abundances may significantly bias the inferred properties of the polluting supernovae. We redetermine the chemical composition of SMSS0313-6708 using 3D NLTE radiative transfer to obtain accurate abundances for Li, Na, Mg, Al, Ca and Fe. The model atoms employ realistic collisional rates, with no calibrated free parameters. We find significantly higher abundances in 3D NLTE than 1D LTE by 0.8 dex for Fe, and 0.5 dex for Mg, Al and Ca, while Li and Na are unaffected to within 0.03 dex. In particular, our upper limit for [Fe/H] is now a factor ten larger, at [Fe/H] < -6.53 (3 sigma), than previous estimates based on <3D> NLTE (i.e. using averaged 3D models). This higher estimate is due to a conservative upper limit estimation, updated NLTE data, and 3D-<3D> NLTE differences, all of which lead to a higher abundance determination. We find that the revised chemical composition of SMSS0313-6708 matches supernova yields for massive progenitors of 20-60 Msol exploding with low energies (1-2 x 10^51 erg), as well as progenitors of 10 Msol with very low explosion energies (< 10^51 erg).

preprint2016arXiv

Atomic diffusion and mixing in old stars VI: The lithium content of M30

The prediction of the PLANCK-constrained primordial lithium abundance in the Universe is in discordance with the observed Li abundances in warm Population II dwarf and subgiant stars. Among the physically best motivated ideas, it has been suggested that this discrepancy can be alleviated if the stars observed today had undergone photospheric depletion of lithium. The cause of this depletion is investigated by accurately tracing the behaviour of the lithium abundances as a function of effective temperature. Globular clusters are ideal laboratories for such an abundance analysis as the relative stellar parameters of their stars can be precisely determined. We performed a homogeneous chemical abundance analysis of 144 stars in the metal-poor globular cluster M30, ranging from the cluster turnoff point to the tip of the red giant branch. NLTE abundances for Li, Ca, and Fe were derived where possible. Stellar parameters were derived by matching isochrones to the observed V vs V-I colour-magnitude diagram. Independent effective temperatures were obtained from automated profile fitting of the Balmer lines and by applying colour-T_eff calibrations to the broadband photometry. Li abundances of the turnoff and early subgiant stars form a thin plateau that is broken off abruptly in the middle of the SGB as a result of the onset of Li dilution caused by the first dredge-up. Abundance trends with effective temperature for Fe and Ca are observed and compared to predictions from stellar structure models including atomic diffusion and ad hoc additional mixing below the surface convection zone. The comparison shows that the stars in M30 are affected by atomic diffusion and additional mixing. After applying a conservative correction for atomic diffusion, we find an initial Li abundance of A(Li) = $2.48\pm0.10$ for the globular cluster M30.

preprint2016arXiv

GRACES observations of young [alpha/Fe]-rich stars

We measure chemical abundance ratios and radial velocities in four massive (i.e., young) [alpha/Fe]-rich red giant stars using high-resolution high-S/N spectra from ESPaDOnS fed by Gemini-GRACES. Our differential analysis ensures that our chemical abundances are on the same scale as the Alves-Brito et al. (2010) study of bulge, thin and thick disk red giants. We confirm that the program stars have enhanced [alpha/Fe] ratios and are slightly metal poor. Aside from lithium enrichment in one object, the program stars exhibit no chemical abundance anomalies when compared to giant stars of similar metallicity throughout the Galaxy. This includes the elements Li, O, Si, Ca, Ti, Cr, Ni, Cu, Ba, La, and Eu. Therefore, there are no obvious chemical signatures that can help to reveal the origin of these unusual stars. While our new observations show that only one star (not the Li-rich object) exhibits a radial velocity variation, simulations indicate that we cannot exclude the possibility that all four could be binaries. In addition, we find that two (possibly three) stars show evidence for an infrared excess, indicative of a debris disk. This is consistent with these young [alpha/Fe]-rich stars being evolved blue stragglers, suggesting their apparent young age is a consequence of a merger or mass transfer. We would expect a binary fraction of ~50% or greater for the entire sample of these stars, but the signs of the circumbinary disk may have been lost since these features can have short timescales. Radial velocity monitoring is needed to confirm the blue straggler origin.

preprint2016arXiv

High precision analysis of the solar twin HIP 100963

Methods. We used the HIRES spectrograph on the Keck I telescope to acquire high-resolution (R $\approx$ 70000) spectra with a high signal-to-noise ratio (S/N $\approx$ 400 - 650 per pixel) of HIP 100963 and the Sun for a differential abundance analysis. We measured the equivalent widths (EWs) of iron lines to determine the stellar parameters by employing the differential spectroscopic equilibrium. We determined the composition of volatile, refractory, and neutron-capture elements through a differential abundance analysis with respect to the Sun. Results. The stellar parameters we found are $T_{\rm{eff}}=5818 \pm 4$ K, log $g = 4.49 \pm 0.01$ dex, $v_{t} = 1.03 \pm 0.01 $ $\rm{km\ {s}}^{-1}$ , and [Fe/H] $ = -\ 0.003 \pm 0.004$ dex. These low errors allow us to compute a precise mass ($1.03^{+0.02}_{-0.01}$ M$_{\odot}$) and age (2.0 $\pm$ 0.4 Gyr), obtained using Yonsei-Yale isochrones. Using our [Y/Mg] ratio, we have determined an age of $2.1 \pm 0.4$ Gyr, in agreement with the age computed using isochrones. Our isochronal age also agrees with the age determined from stellar activity (2.4 $\pm$ 0.3 Gyr). We study the abundance pattern with condensation temperature ($\rm{T_{cond}}$) taking corrections by the GCE into account. We show that the enhancements of neutron-capture elements are explained by contributions from both the $s$- and $r$-process. The lithium abundance follows the tight Li-age correlation seen in other solar twins. Conclusions. We confirm that HIP 100963 is a solar twin and demonstrate that its abundance pattern is about solar after corrections for GCE. The star also shows enrichment in $s-$ and $r$-process elements, as well as depletion in lithium that is caused by stellar evolution.

preprint2016arXiv

Implications of solar wind measurements for solar models and composition

We critically examine recent claims of a high solar metallicity by von Steiger \& Zurbuchen (2016) based on \textit{in situ} measurements of the solar wind, rather than the standard spectroscopically-inferred abundances (Asplund et al. 2009). We test the claim by Vagnozzi et al. (2016) that a composition based on the solar wind enables one to construct a standard solar model in agreement with helioseismological observations and thus solve the decades-old solar modelling problem. We show that, although some helioseismological observables are improved compared to models computed with spectroscopic abundances, most are in fact worse. The high abundance of refractory elements leads to an overproduction of neutrinos, with a predicted $^8$B flux that is nearly twice its observed value, and $^7$Be and CNO fluxes that are experimentally ruled out at high confidence. A combined likelihood analysis shows that models using the vSZ16 abundances fare worse than AGSS09 despite a higher metallicity. We also present astrophysical and spectroscopic arguments showing the vSZ16 composition to be an implausible representation of the solar interior, identifying the first ionisation potential effect in the outer solar atmosphere and wind as the likely culprit.

preprint2016arXiv

Nominal values for selected solar and planetary quantities: IAU 2015 Resolution B3

In this brief communication we provide the rationale for, and the outcome of the International Astronomical Union (IAU) resolution vote at the XXIX-th General Assembly in Honolulu, Hawaii, in 2015, on recommended nominal conversion constants for selected solar and planetary properties. The problem addressed by the resolution is a lack of established conversion constants between solar and planetary values and SI units: a missing standard has caused a proliferation of solar values (e.g., solar radius, solar irradiance, solar luminosity, solar effective temperature and solar mass parameter) in the literature, with cited solar values typically based on best estimates at the time of paper writing. As precision of observations increases, a set of consistent values becomes increasingly important. To address this, an IAU Working Group on Nominal Units for Stellar and Planetary Astronomy formed in 2011, uniting experts from the solar, stellar, planetary, exoplanetary and fundamental astronomy, as well as from general standards fields to converge on optimal values for nominal conversion constants. The effort resulted in the IAU 2015 Resolution B3, passed at the IAU General Assembly by a large majority. The resolution recommends the use of nominal solar and planetary values, which are by definition exact and are expressed in SI units. These nominal values should be understood as conversion factors only, not as the true solar/planetary properties or current best estimates. Authors and journal editors are urged to join in using the standard values set forth by this resolution in future work and publications to help minimize further confusion.

preprint2016arXiv

The Detailed Science Case for the Maunakea Spectroscopic Explorer: the Composition and Dynamics of the Faint Universe

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

preprint2016arXiv

The EMBLA Survey -- Metal-poor stars in the Galactic bulge

Cosmological models predict the oldest stars in the Galaxy should be found closest to the centre of the potential well, in the bulge. The EMBLA Survey successfully searched for these old, metal-poor stars by making use of the distinctive SkyMapper photometric filters to discover candidate metal-poor stars in the bulge. Their metal-poor nature was then confirmed using the AAOmega spectrograph on the AAT. Here we present an abundance analysis of 10 bulge stars with -2.8<[Fe/H]<-1.7 from MIKE/Magellan observations, in total determining the abundances of 22 elements. Combining these results with our previous high-resolution data taken as part of the Gaia-ESO Survey, we have started to put together a picture of the chemical and kinematic nature of the most metal-poor stars in the bulge. The currently available kinematic data is consistent with the stars belonging to the bulge, although more accurate measurements are needed to constrain the stars' orbits. The chemistry of these bulge stars deviates from that found in halo stars of the same metallicity. Two notable differences are the absence of carbon-enhanced metal-poor bulge stars, and the alpha-element abundances exhibit a large intrinsic scatter and include stars which are underabundant in these typically enhanced elements.

preprint2016arXiv

The Gaia-ESO survey: Hydrogen lines in red giants directly trace stellar mass

Red giant stars are perhaps the most important type of stars for Galactic and extra-galactic archaeology: they are luminous, occur in all stellar populations, and their surface temperatures allow precise abundance determinations for many different chemical elements. Yet, the full star formation and enrichment history of a galaxy can be traced directly only if two key observables can be determined for large stellar samples - age and chemical composition. While spectroscopy is a powerful method to analyse the detailed abundances of stars, stellar ages are the "missing link in the chain", since they are not a direct observable. However, spectroscopy should be able to estimate stellar masses, which for red giants directly infer ages provided their chemical composition is known. Here we establish a new empirical relation between the shape of the hydrogen line in the observed spectra of red giants and stellar mass determined from asteroseismology. The relation allows to determine stellar masses and ages with the accuracy of 10-15%. The method can be used with confidence for stars in the following range of stellar parameters: 4000 < Teff < 5000 K, 0.5 < log g < 3.5, -2.0 < [Fe/H] < 0.3, and luminosities log L/LSun < 2.5. Our analysis provides observational evidence that the Halpha spectral characteristics of red giant stars are tightly correlated with their mass and therefore their age. We also show that the method samples well all stellar populations with ages above 1 Gyr. Targeting bright giants, the method allows to obtain simultaneous age and chemical abundance information far deeper than would be possible with asteroseismology, extending the possible survey volume to remote regions of the Milky Way and even to neighbouring galaxies like Andromeda or the Magellanic Clouds already with present instrumentation, like VLT and Keck facilities.

preprint2016arXiv

The GALAH survey: The data reduction pipeline

We present the data reduction procedures being used by the GALAH survey, carried out with the HERMES fibre-fed, multi-object spectrograph on the 3.9~m Anglo-Australian Telescope. GALAH is a unique survey, targeting 1 million stars brighter than magnitude V=14 at a resolution of 28,000 with a goal to measure the abundances of 29 elements. Such a large number of high resolution spectra necessitates the development of a reduction pipeline optimized for speed, accuracy, and consistency. We outline the design and structure of the Iraf-based reduction pipeline that we developed, specifically for GALAH, to produce fully calibrated spectra aimed for subsequent stellar atmospheric parameter estimation. The pipeline takes advantage of existing Iraf routines and other readily available software so as to be simple to maintain, testable and reliable. A radial velocity and stellar atmospheric parameter estimator code is also presented, which is used for further data analysis and yields a useful verification of the reduction quality. We have used this estimator to quantify the data quality of GALAH for fibre cross-talk level ($\lesssim0.5$%) and scattered light ($\sim5$ counts in a typical 20 minutes exposure), resolution across the field, sky spectrum properties, wavelength solution reliability (better than $1$ $\mathrm{km\ s^{-1}}$ accuracy) and radial velocity precision.

preprint2016arXiv

The Solar Twin Planet Search: IV. The Sun as a typical rotator and evidence for a new rotational braking law for Sun-like stars

It is still unclear how common the Sun is when compared to other similar stars in regards to some of its physical properties, such as rotation. Considering that gyrochronology relations are widely used today to estimate ages of stars in the main sequence, and that the Sun is used to calibrate it, it is crucial to assess if these procedures are acceptable. We analyze the rotational velocities -- limited by the unknown rotation axis inclination angle -- of an unprecedented large sample of solar twins in order to study the rotational evolution of Sun-like stars, and assess if the Sun is a typical rotator. We use high-resolution ($R = 115000$) spectra obtained with the HARPS spectrograph and ESO's 3.6 m telescope at La Silla Observatory. The projected rotational velocities for 82 solar twins are estimated by line profile fitting with synthetic spectra. Macroturbulence velocities are inferred from a prescription that accurately reflects their dependence with effective temperature and luminosity of the stars. Our sample of solar twins include some spectroscopic binaries with enhanced rotational velocities, and we do not find any non-spectroscopic binaries with unusually high rotation velocities. We verified that the Sun does not have a peculiar rotation, but the solar twins exhibit rotational velocities that depart from the Skumanich relation. The Sun is a regular rotator when compared to solar twins with a similar age. Additionally, we obtain a rotational braking law that better describes the stars in our sample ($v \propto t^{-0.6}$) in contrast to previous, often-used scalings.

preprint2016arXiv

The Solar Twin Planet Search. V. Close-in, low-mass planet candidates and evidence of planet accretion in the solar twin HIP 68468

[Methods]. We obtained high-precision radial velocities with HARPS on the ESO 3.6 m telescope and determined precise stellar elemental abundances (~0.01 dex) using MIKE spectra on the Magellan 6.5m telescope. [Results]. Our data indicate the presence of a planet with a minimum mass of 26 Earth masses around the solar twin HIP 68468. The planet is a super-Neptune, but unlike the distant Neptune in our solar system (30 AU), HIP 68468c is close-in, with a semi-major axis of 0.66 AU, similar to that of Venus. The data also suggest the presence of a super-Earth with a minimum mass of 2.9 Earth masses at 0.03 AU; if the planet is confirmed, it will be the fifth least massive radial velocity planet discovery to date and the first super-Earth around a solar twin. Both isochrones (5.9 Gyr) and the abundance ratio [Y/Mg] (6.4 Gyr) indicate an age of about 6 billion years. The star is enhanced in refractory elements when compared to the Sun, and the refractory enrichment is even stronger after corrections for Galactic chemical evolution. We determined a NLTE Li abundance of 1.52 dex, which is four times higher than what would be expected for the age of HIP 68468. The older age is also supported by the low log(R'HK) (-5.05) and low jitter. Engulfment of a rocky planet of 6 Earth masses can explain the enhancement in both lithium and the refractory elements. [Conclusions]. The super-Neptune planet candidate is too massive for in situ formation, and therefore its current location is most likely the result of planet migration that could also have driven other planets towards its host star, enhancing thus the abundance of lithium and refractory elements in HIP 68468. The intriguing evidence of planet accretion warrants further observations to verify the existence of the planets that are indicated by our data and to better constrain the nature of the planetary system around this unique star.

preprint2016arXiv

Using the Multi-Object Adaptive Optics demonstrator RAVEN to observe metal-poor stars in and towards the Galactic Centre

The chemical abundances for five metal-poor stars in and towards the Galactic bulge have been determined from H-band infrared spectroscopy taken with the RAVEN multi-object adaptive optics science demonstrator and the IRCS spectrograph at the Subaru 8.2-m telescope. Three of these stars are in the Galactic bulge and have metallicities between -2.1 < [Fe/H] < -1.5, and high [alpha/Fe] ~+0.3, typical of Galactic disk and bulge stars in this metallicity range; [Al/Fe] and [N/Fe] are also high, whereas [C/Fe] < +0.3. An examination of their orbits suggests that two of these stars may be confined to the Galactic bulge and one is a halo trespasser, though proper motion values used to calculate orbits are quite uncertain. An additional two stars in the globular cluster M22 show [Fe/H] values consistent to within 1 sigma, although one of these two stars has [Fe/H] = -2.01 +/- 0.09, which is on the low end for this cluster. The [alpha/Fe] and [Ni/Fe] values differ by 2 sigma, with the most metal-poor star showing significantly higher values for these elements. M22 is known to show element abundance variations, consistent with a multi-population scenario (i.e. Marino et al. 2009, 2011; Alves-Brito et al. 2012) though our results cannot discriminate this clearly given our abundance uncertainties. This is the first science demonstration of multi-object adaptive optics with high resolution infrared spectroscopy, and we also discuss the feasibility of this technique for use in the upcoming era of 30-m class telescope facilities.

preprint2015arXiv

First Light Results from the Hermes Spectrograph at the AAT

The High Efficiency and Resolution Multi Element Spectrograph, HERMES, is a facility-class optical spectrograph for the Anglo-Australian Telescope (AAT). It is designed primarily for Galactic Archaeology, the first major attempt to create a detailed understanding of galaxy formation and evolution by studying the history of our own galaxy, the Milky Way. The goal of the GALAH survey is to reconstruct the mass assembly history of the Milky Way through a detailed chemical abundance study of one million stars. The spectrograph is based at the AAT and is fed by the existing 2dF robotic fiber positioning system. The spectrograph uses volume phase holographic gratings to achieve a spectral resolving power of 28,000 in standard mode and also provides a high-resolution mode ranging between 40,000 and 50,000 using a slit mask. The GALAH survey requires an SNR greater than 100 for a star brightness of V ?= 14 in an exposure time of one hour. The total spectral coverage of the four channels is about 100 nm between 370 and 1000 nm for up to 392 simultaneous targets within the 2-degree field of view. HERMES has been commissioned over three runs, during bright time in October, November, and December 2013, in parallel with the beginning of the GALAH pilot survey, which started in November 2013. We present the first-light results from the commissioning run and the beginning of the GALAH survey, including performance results such as throughput and resolution, as well as instrument reliability.

preprint2015arXiv

High-Resolution Spectroscopic Study of Extremely Metal-Poor Star Candidates from the SkyMapper Survey

The SkyMapper Southern Sky Survey is carrying out a search for the most metal-poor stars in the Galaxy. It identifies candidates by way of its unique filter set that allows for estimation of stellar atmospheric parameters. The set includes a narrow filter centered on the Ca II K 3933A line, enabling a robust estimate of stellar metallicity. Promising candidates are then confirmed with spectroscopy. We present the analysis of Magellan-MIKE high-resolution spectroscopy of 122 metal-poor stars found by SkyMapper in the first two years of commissioning observations. 41 stars have [Fe/H] <= -3.0. Nine have [Fe/H] <= -3.5, with three at [Fe/H] ~ -4. A 1D LTE abundance analysis of the elements Li, C, Na, Mg, Al, Si, Ca, Sc, Ti, Cr, Mn, Co, Ni, Zn, Sr, Ba and Eu shows these stars have [X/Fe] ratios typical of other halo stars. One star with low [X/Fe] values appears to be "Fe-enhanced," while another star has an extremely large [Sr/Ba] ratio: >2. Only one other star is known to have a comparable value. Seven stars are "CEMP-no" stars ([C/Fe] > 0.7, [Ba/Fe] < 0). 21 stars exhibit mild r-process element enhancements (0.3 <=[Eu/Fe] < 1.0), while four stars have [Eu/Fe] >= 1.0. These results demonstrate the ability to identify extremely metal-poor stars from SkyMapper photometry, pointing to increased sample sizes and a better characterization of the metal-poor tail of the halo metallicity distribution function in the future.

preprint2015arXiv

HIP 10725: The First Solar Twin/Analogue Field Blue Straggler

[Context]. Blue stragglers are easy to identify in globular clusters, but are much harder to identify in the field. Here we present the serendipitous discovery of one field blue straggler, HIP 10725, that closely matches the Sun in mass and age, but with a metallicity slightly lower than the Sun's. [Methods]. We employ high resolution (R $\sim 10^5$) high S/N (330) VLT/UVES spectra to perform a differential abundance analysis of the solar analogue HIP 10725. Radial velocities obtained by other instruments were also used to check for binarity. We also study its chromospheric activity, age and rotational velocity. [Results]. We find that HIP 10725 is severely depleted in beryllium ([Be/H] <= -1.2 dex) for its stellar parameters and age. The abundances relative to solar of the elements with Z <= 30 show a correlation with condensation temperature and the neutron capture elements produced by the s-process are greatly enhanced, while the r-process elements seem normal. We found its projected rotational velocity (v sin i = 3.3+/-0.1 km/s) to be significantly larger than solar, and incompatible with its isochrone-derived age. Radial velocity monitoring shows that the star has a binary companion. [Conclusions]. Based on the high s-process element enhancements and low beryllium abundance, we suggest that HIP 10725 has been polluted by mass-transfer from an AGB star, probably with initial mass of about 2 M_Sun. The radial velocity variations suggest the presence of an unseen binary companion, probably the remnant of a former AGB star. Isochrones predict a solar-age star, but this is in disagreement with the high projected rotational velocity and high chromospheric activity. We conclude that HIP 10725 is a field blue straggler, rejuvenated by the mass transfer process of its former AGB companion.

preprint2015arXiv

Nucleosynthesis in a Primordial Supernova: Carbon and Oxygen Abundances in SMSS J031300.36-670839.31

SMSS J031300.36-670839.3 (hereafter SM0313-6708) is a sub-giant halo star, with no detectable Fe lines and large overabundances of C and Mg relative to Ca. We obtained VLT-UVES spectra extending to 3060 Angstroms showing strong OH A-X band lines enabling an oxygen abundance to be derived. The OH A-X band lines in SM0313-6708 are much stronger than the CH C-X band lines. Spectrum synthesis fits indicate an [O/C] ratio of 0.02 +- 0.175. Our high S/N UVES data also enabled us to lower the Fe abundance limit to [Fe/H]{3D},NLTE < -7.52 (3 sigma). These data support our previous suggestion that the star formed from the iron-poor ejecta of a single massive star Population III supernova.

preprint2015arXiv

The first high precision differential abundance analysis of Extremely Metal Poor stars

Context: Studies of extremely metal-poor stars indicate that chemical abundance ratios [X/Fe] have an rms scatter as low as 0.05 dex (12 \%). It remains unclear whether this reflects observational uncertainties or intrinsic astrophysical scatter arising from physical conditions in the ISM at early times. Aims: Measure differential chemical abundance ratios in extremely metal-poor stars to investigate the limits of precision and to understand whether cosmic scatter or observational errors are dominant. Methods: We used high resolution (R $\sim 95,000$) and high S/N (S/N $= 700$ at 5000$Å$) HIRES/Keck spectra, to determine high precision differential abundances between two extremely metal-poor stars through a line-by-line differential approach. We determined stellar parameters for the star G64-37 with respect to the standard star G64-12. We performed EW measurements for the two stars for the lines recognized in both stars and performed spectral synthesis to study the carbon abundances. Results: The differential approach allowed us to obtain errors of $σ$(T$_{eff}$ ) $=$ 27 K, $σ$(log $g$) $=$ 0.06 dex, $σ$([Fe/H]) $=$ 0.02 dex and $σ$(v$_{t}$ ) $=$ 0.06 kms$^{-1}$. We estimated relative chemical abundances with a precision as low as $σ$([X/Fe]) $\approx$ 0.01 dex. The small uncertainties demonstrate that there are genuine abundance differences larger than the measurement errors. The observed Li difference can not be explained by the difference in mass, because the less massive star has more Li. Conclusions: It is possible to achieve an abundance precision around $\approx$ 0.01-0.05 dex for extremely metal-poor stars, opening new windows on the study of the early chemical evolution of the Galaxy.

preprint2014arXiv

Does SEGUE/SDSS indicate a dual Galactic halo?

We re-examine recent claims of observational evidence for a dual Galactic halo in SEGUE/SDSS data, and trace them back to improper error treatment and neglect of selection effects. In particular, the detection of a vertical abundance gradient in the halo can be explained as a metallicity bias in distance. A similar bias, and the impact of disk contamination, affect the sample of blue horizontal branch stars. These examples highlight why non-volume complete samples require forward-modelling from theoretical models or extensive bias-corrections. We also show how observational uncertainties produce the specific non-Gaussianity in the observed azimuthal velocity distribution of halo stars, which can be erroneously identified as two Gaussian components. A single kinematic component yields an excellent fit to the observed data, when we model the measurement process including distance uncertainties. Furthermore, we show that sample differences in proper motion space are the direct consequence of kinematic cuts, and are enhanced when distance estimates are less accurate. Thus, their presence is neither a proof for a separate population, nor a measure of reliability for the applied distances. We conclude that currently there is no evidence from SEGUE/SDSS that would favour a dual Galactic halo over a single halo full of substructure.

preprint2014arXiv

HIP 114328: a new refractory-poor and Li-poor solar twin

[Context]. The standard solar model fails to predict the very low lithium abundance in the Sun, which is much lower than the proto-solar nebula. This Li problem has been debated for decades, and it has been ascribed either to planet formation or to secular stellar depletion. In order to test the evolution of Li, it is important to find solar twins in a range of ages. Also, the study of stars similar to the Sun is relevant in relation to the signature of terrestrial planet formation around the Sun. [Methods]. We acquired high-resolution (R=110,000), high S/N (~300) ESO/VLT UVES spectra of several solar twin candidates and the Sun (as reflected from the asteroid Juno). Among the solar twin candidates we identify HIP 114328 as a solar twin and perform a differential line-by-line abundance analysis of this star relative to the Sun. [Results]. HIP 114328 has stellar parameters Teff = 5785+/-10 K, log g = 4.38+/-0.03, [Fe/H] = -0.022+/-0.009, and a microturbulent velocity 0.05+/-0.03 km/s higher than solar. The differential analysis shows that this star is chemically very similar to the Sun. The refractory elements seem even slightly more depleted than in the Sun, meaning that HIP 114328 may be as likely to form terrestrial planets as the Sun. HIP 114328 is about 2 Gyr older than the Sun, and is thus the second oldest solar twin analyzed at high precision. It has a Li abundance of A(Li)NLTE <= 0.46, which is about 4 times lower than in the Sun (A(Li)NLTE = 1.07 dex), but close to the oldest solar twin known, HIP 102152. [Conclusions]. Based on the lower abundances of refractory elements when compared to other solar twins, HIP 114328 seems an excellent candidate to host rocky planets. The low Li abundance of this star is consistent with its old age and fits very well the emerging Li-age relation among solar twins of different ages.

preprint2014arXiv

Improvements to Stellar Structure Models, Based on a Grid of 3D Convection Simulations. I. $T(τ)$-Relations

Relations between temperature, T, and optical depth, tau, are often used for describing the photospheric transition from optically thick to optically thin in stellar structure models. We show that this is well justified, but also that currently used T(tau) relations are often inconsistent with their implementation. As an outer boundary condition on the system of stellar structure equations, T(tau) relations have an undue effect on the overall structure of stars. In this age of precision asteroseismology, we need to re-assess both the method for computing and for implementing T(tau) relations, and the assumptions they rest on. We develop a formulation for proper and consistent evaluation of T(tau) relations from arbitrary 1D or 3D stellar atmospheres, and for their implementation in stellar structure and evolution models. We extract radiative T(tau) relations, as described by our new formulation, from 3D simulations of convection in deep stellar atmospheres of late-type stars from dwarfs to giants. These simulations employ realistic opacities and equation of state, and account for line-blanketing. For comparison, we also extract T(tau) relations from 1D MARCS model atmospheres using the same formulation. T(tau)-relations from our grid of 3D convection simulations display a larger range of behaviours with surface gravity, compared with those of conventional theoretical 1D hydrostatic atmosphere models. Based on this, we recommend no longer to use scaled solar T(tau) relations. Files with T(tau) relations for our grid of simulations are made available to the community, together with routines for interpolating in this irregular grid. We also provide matching tables of atmospheric opacity, for consistent implementation in stellar structure models.

preprint2014arXiv

Improvements to stellar structure models, based on a grid of 3D convection simulations. II. Calibrating the mixing-length formulation

We perform a calibration of the mixing length of convection in stellar structure models against realistic 3D radiation-coupled hydrodynamics (RHD) simulations of convection in stellar surface layers, determining the adiabat deep in convective stellar envelopes. The mixing-length parameter $α$ is calibrated by matching averages of the 3D simulations to 1D stellar envelope models, ensuring identical atomic physics in the two cases. This is done for a previously published grid of solar-metallicity convection simulations, covering from 4200 K to 6900 K on the main sequence, and 4300-5000 K for giants with logg=2.2. Our calibration results in an $α$ varying from 1.6 for the warmest dwarf, which is just cool enough to admit a convective envelope, and up to 2.05 for the coolest dwarfs in our grid. In between these is a triangular plateau of $α$ ~ 1.76. The Sun is located on this plateau and has seen little change during its evolution so far. When stars ascend the giant branch, they largely do so along tracks of constant $α$, with $α$ decreasing with increasing mass.

preprint2014arXiv

NGC 6522: A typical globular cluster in the Galactic bulge without signatures of rapidly rotating Population III stars

We present an abundance analysis of eight potential member stars of the old Galactic bulge globular cluster NGC6522. The same stars have previously been studied by Chiappini et al. (2011), who found very high abundances of the slow neutron capture elements compared with other clusters and field stars of similar metallicity, which they interpreted as reflecting nucleosynthesis in rapidly rotating, massive Population III stars. In contrast to their analysis, we do not find any unusual enhancements of the neutron capture elements Sr, Y, Ba and Eu and conclude that previous claims result mainly from not properly accounting for blending lines. Instead we find NGC6522 to be an unremarkable globular cluster with comparable abundance trends to other Galactic globular clusters at the same metallicity ([Fe/H] = -1.15 +/- 0.16). The stars are also chemically similar to halo and bulge field stars at the same metallicity, spanning a small range in [Y/Ba] and with normal α-element abundances. We thus find no observational evidence for any chemical signatures of rapidly rotating Population III stars in NGC 6522.

preprint2014arXiv

Planet transit and stellar granulation detection with interferometry

Aims. We used realistic three-dimensional (3D) radiative hydrodynamical (RHD) simulations from the Stagger-grid and synthetic images computed with the radiative transfer code Optim3D to provide interferometric observables to extract the signature of stellar granulation and transiting planets. Methods. We computed intensity maps from RHD simulations for twelve interferometric instruments covering wavelengths ranging from optical to infrared. The stellar surface asymmetries in the brightness distribution mostly affect closure phases. We compared the closure phases of the system star with a transiting planet and the star alone and considered the impact of magnetic spots constructing a hypothetical starspots image. Results. All the simulations show departure from the axisymmetric case at all wavelengths. We presented two possible targets (Beta Com and Procyon) and found that departures up to 16 deg can be detected on the 3rd lobe and higher. In particular, MIRC is the most appropriate instrument because it combines good UV coverage and long baselines. Moreover, we explored the impact of convection on interferometric planet signature for three prototypes of planets. It is possible to disentangle the signature of the planet at particular wavelengths (either in the infrared or in the optical) by comparing the closure phases of the star at difference phases of the planetary transit. Conclusions. The detection and characterisation of planets must be based on a comprehensive knowledge of the host star; this includes the detailed study of the stellar surface convection with interferometric techniques. In this context, RHD simulations are crucial to reach this aim. We emphasize that interferometric observations should be pushed at high spatial frequencies by accumulating observations on closure phases at short and long baselines.

preprint2014arXiv

The elemental composition of the Sun I. The intermediate mass elements Na to Ca

The composition of the Sun is an essential piece of reference data for astronomy, cosmology, astroparticle, space and geo-physics. This article, dealing with the intermediate-mass elements Na to Ca, is the first in a series describing the comprehensive re-determination of the solar composition. In this series we severely scrutinise all ingredients of the analysis across all elements, to obtain the most accurate, homogeneous and reliable results possible. We employ a highly realistic 3D hydrodynamic solar photospheric model, which has successfully passed an arsenal of observational diagnostics. To quantify systematic errors, we repeat the analysis with three 1D hydrostatic model atmospheres (MARCS, MISS and Holweger & Müller 1974) and a horizontally and temporally-averaged version of the 3D model ($\langle$3D$\rangle$). We account for departures from LTE wherever possible. We have scoured the literature for the best transition probabilities, partition functions, hyperfine and other data, and stringently checked all observed profiles for blends. Our final 3D+NLTE abundances are: $\logε_{\mathrm{Na}}=6.21\pm0.04$, $\logε_{\mathrm{Mg}}=7.59\pm0.04$, $\logε_{\mathrm{Al}}=6.43\pm0.04$, $\logε_{\mathrm{Si}}=7.51\pm0.03$, $\logε_{\mathrm{P}}=5.41\pm0.03$, $\log ε_{\mathrm{S}}=7.13\pm0.03$, $\logε_{\mathrm{K}}=5.04\pm0.05$ and $\logε_{\mathrm{Ca}}=6.32\pm0.03$. The uncertainties include both statistical and systematic errors. Our results are systematically smaller than most previous ones with the 1D semi-empirical Holweger & Müller model. The $\langle$3D$\rangle$ model returns abundances very similar to the full 3D calculations. This analysis provides a complete description and a slight update of the Na to Ca results presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), with full details of all lines and input data.

preprint2014arXiv

The elemental composition of the Sun II. The iron group elements Sc to Ni

We redetermine the abundances of all iron group nuclei in the Sun, based on neutral and singly-ionised lines of Sc, Ti, V, Mn, Fe, Co and Ni in the solar spectrum. We employ a realistic 3D hydrodynamic model solar atmosphere, corrections for departures from local thermodynamic equilibrium (NLTE), stringent line selection procedures and high quality observational data. We have scoured the literature for the best quality oscillator strengths, hyperfine constants and isotopic separations available for our chosen lines. We find $\log ε_\mathrm{Sc}=3.16\pm0.04$, $\log ε_\mathrm{Ti}=4.93\pm0.04$, $\log ε_\mathrm{V}=3.89\pm0.08$, $\log ε_\mathrm{Cr}=5.62\pm0.04$, $\log ε_\mathrm{Mn}=5.42\pm0.04$, $\log ε_\mathrm{Fe}=7.47\pm0.04$, $\log ε_\mathrm{Co}=4.93\pm0.05$ and $\log ε_\mathrm{Ni}=6.20\pm0.04$. Our uncertainties factor in both statistical and systematic errors (the latter estimated for possible errors in the model atmospheres and NLTE line formation). The new abundances are generally in good agreement with the CI meteoritic abundances but with some notable exceptions. This analysis constitutes both a full exposition and a slight update of the preliminary results we presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), including full line lists and details of all input data we employed.

preprint2014arXiv

The elemental composition of the Sun III. The heavy elements Cu to Th

We re-evaluate the abundances of the elements in the Sun from copper ($Z=29$) to thorium ($Z=90$). Our results are mostly based on neutral and singly-ionised lines in the solar spectrum. We use the latest 3D hydrodynamic solar model atmosphere, and in a few cases also correct for departures from local thermodynamic equilibrium (LTE) using non-LTE (NLTE) calculations performed in 1D. In order to minimise statistical and systematic uncertainties, we make stringent line selections, employ the highest-quality observational data and carefully assess oscillator strengths, hyperfine constants and isotopic separations available in the literature, for every line included in our analysis. Our results are typically in good agreement with the abundances in the most pristine meteorites, but there are some interesting exceptions. This analysis constitutes both a full exposition and a slight update of the relevant parts of the preliminary results we presented in Asplund, Grevesse, Sauval & Scott (arXiv:0909.0948), including full line lists and details of all input data that we have employed.

preprint2014arXiv

The Stagger-grid: A grid of 3D stellar atmosphere models - III. The relation to mixing-length convection theory

We investigate the relation between 1D atmosphere models that rely on the mixing length theory and models based on full 3D radiative hydrodynamic (RHD) calculations to describe convection in the envelopes of late-type stars. The adiabatic entropy value of the deep convection zone, s_bot, and the entropy jump, Δs, determined from the 3D RHD models, are matched with the mixing length parameter, α_MLT, from 1D hydrostatic atmosphere models with identical microphysics (opacities and equation-of-state). We also derive the mass mixing length, α_m, and the vertical correlation length of the vertical velocity, C[v_z,v_z], directly from the 3D hydrodynamical simulations of stellar subsurface convection. The calibrated mixing length parameter for the Sun is α_MLT (s_bot) = 1.98. For different stellar parameters, α_MLT varies systematically in the range of 1.7 - 2.4. In particular, α_MLT decreases towards higher effective temperature, lower surface gravity and higher metallicity. We find equivalent results for α_MLT (Δs). Also, we find a tight correlation between the mixing length parameter and the inverse entropy jump. We derive an analytical expression from the hydrodynamic mean field equations that motivates the relation to the mass mixing length, α_m, and find that it exhibits qualitatively a similar variation with stellar parameter (between 1.6 and 2.4) with a solar value of α_m = 1.83. The vertical correlation length scaled with the pressure scale height yields for the Sun 1.71, but displays only a small systematic variation with stellar parameters, the correlation length slightly increasing with Teff. We derive mixing length parameters for various stellar parameters that can be used to replace a constant value. Within any convective envelope, α_m and related quantities vary a lot.

preprint2014arXiv

The Stagger-grid: A grid of 3D stellar atmosphere models - IV. Limb darkening coefficients

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. We have computed full 3D LTE synthetic spectra based on 3D radiative hydrodynamic atmosphere models from the Stagger-grid. From the resulting intensities at different wavelength, 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. The 3D models are often found to predict steeper limb darkening 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. 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.

preprint2014arXiv

The Stagger-grid: A grid of 3D stellar atmosphere models - V. Fe line shapes, shifts and asymmetries

We present a theoretical study of the effects and signatures of realistic velocity field and atmospheric inhomogeneities associated with convective motions at the surface of cool late-type stars on the emergent profiles of iron spectral lines for a large range in stellar parameters. We compute 3D spectral line flux profiles under the assumption of local thermodynamic equilibrium (LTE) by employing state-of-the-art, time-dependent, 3D, radiative-hydrodynamical atmosphere models from the Stagger-grid. A set of 35 real unblended, optical FeI and FeII lines of varying excitation potential are considered. Additionally, fictitious Fe i and Fe ii lines (5000A and 0, 2, 4 eV) are used to construct general curves of growth and enable comparison of line profiles with the same line strength to illustrate systematical trends stemming from the intrinsic structural differences among 3D model atmospheres with different stellar parameters. Theoretical line shifts and bisectors are derived to analyze the shapes, shifts, and asymmetries imprinted in the full 3D line profiles emerging self-consistently from the convective simulations with velocity fields and atmospheric inhomogeneities. We find systematic variations in line strength, shift, width, and bisectors, that can be related to the respective physical conditions at the height of the line formation in the stellar atmospheric environment, in particular the amplitude of the vertical velocity field. Line shifts and asymmetries arise due to the presence of convective velocities and the granulation pattern that are ubiquitously found in observed stellar spectra of cool stars.

preprint2014arXiv

The Stagger-grid: A grid of 3D stellar atmosphere models - VI. Surface appearance of stellar granulation

In the surface layers of late-type stars, stellar convection is manifested with its typical granulation pattern due to the presence of convective motions. The resulting photospheric up- and downflows leave imprints in the observed spectral line profiles. We perform a careful statistical analysis of stellar granulation and its properties for different stellar parameters. We employ realistic 3D radiative hydrodynamic (RHD) simulations of surface convection from the Stagger-grid, a comprehensive grid of atmosphere models that covers a large parameter space in terms of Teff, logg, and [Fe/H]. Individual granules are detected from the (bolometric) intensity maps at disk center with an efficient granulation pattern recognition algorithm. From these we derive their respective properties: diameter, fractal dimension (area-perimeter relation), geometry, topology, variation of intensity, temperature, density and velocity with granule size. Also, the correlation of the physical properties at the optical surface are studied. We find in all of our 3D RHD simulations stellar granulation patterns imprinted, which are qualitatively similar to the solar case, despite the large differences in stellar parameters. The granules exhibit a large range in size, which can be divided into two groups - smaller and larger granules - by the mean granule size. These are distinct in their properties: smaller granules are regular shaped and dimmer, while the larger ones are increasingly irregular and more complex in their shapes and distribution in intensity contrast. This is reflected in their fractal dimensions, which is close to unity for the smaller granules, and close to two for larger granules, which is due to the fragmentation of granules. Stellar surface convection seems to operate scale-invariant over a large range in stellar parameters, which translates into a self-similar stellar granulation pattern.

preprint2013arXiv

A Grid of 3D Stellar Atmosphere Models of Solar Metallicity: I. General Properties, Granulation and Atmospheric Expansion

Present grids of stellar atmosphere models are the workhorses in interpreting stellar observations, and determining their fundamental parameters. These models rely on greatly simplified models of convection, however, lending less predictive power to such models of late type stars. We present a grid of improved and more reliable stellar atmosphere models of late type stars, based on deep, 3D, convective, stellar atmosphere simulations. This grid is to be used in general for interpreting observations, and improve stellar and asteroseismic modeling. We solve the Navier Stokes equations in 3D and concurrent with the radiative transfer equation, for a range of atmospheric parameters, covering most of stellar evolution with convection at the surface. We emphasize use of the best available atomic physics for quantitative predictions and comparisons with observations. We present granulation size, convective expansion of the acoustic cavity, asymptotic adiabat, as function of atmospheric parameters. These and other results are also available in electronic form.

preprint2013arXiv

High Precision Abundances of the Old Solar Twin HIP 102152: Insights on Li Depletion from the Oldest Sun

We present the first detailed chemical abundance analysis of the old 8.2 Gyr solar twin, HIP 102152. We derive differential abundances of 21 elements relative to the Sun with precisions as high as 0.004 dex ($\lesssim$1%), using ultra high-resolution (R = 110,000), high S/N UVES spectra obtained on the 8.2-m Very Large Telescope. Our determined metallicity of HIP 102152 is [Fe/H] = -0.013 $\pm$ 0.004. The atmospheric parameters of the star were determined to be 54 K cooler than the Sun, 0.09 dex lower in surface gravity, and a microturbulence identical to our derived solar value. Elemental abundance ratios examined vs. dust condensation temperature reveal a solar abundance pattern for this star, in contrast to most solar twins. The abundance pattern of HIP 02152 appears to be the most similar to solar of any known solar twin. Abundances of the younger, 2.9 Gyr solar twin, 18 Sco, were also determined from UVES spectra to serve as a comparison for HIP 102152. The solar chemical pattern of HIP 102152 makes it a potential candidate to host terrestrial planets, which is reinforced by the lack of giant planets in its terrestrial planet region. The following non-local thermodynamic equilibrium Li abundances were obtained for HIP 102152, 18 Sco, and the Sun: log $ε$ (Li) = 0.48 $\pm$ 0.07, 1.62 $\pm$ 0.02, and 1.07 $\pm$ 0.02, respectively. The Li abundance of HIP 102152 is the lowest reported to date for a solar twin, and allows us to consider an emerging, tightly constrained Li-age trend for solar twin stars.

preprint2013arXiv

High Precision Differential Abundance Measurements in Globular Clusters: Chemical Inhomogeneities in NGC 6752

We report on a strictly differential line-by-line analysis of high quality UVES spectra of bright giants in the metal-poor globular cluster NGC 6752. We achieved high precision differential chemical abundance measurements for Fe, Na, Si, Ca, Ti, Cr, Ni, Zn, Y, Zr, Ba, La, Ce, Pr, Nd, Sm, Eu and Dy with uncertainties as low as ~0.01 dex (~2%). We obtained the following main results. (1) The observed abundance dispersions are a factor of ~2 larger than the average measurement uncertainty. (2) There are positive correlations, of high statistical significance, between all elements and Na. (3) For any pair of elements, there are positive correlations of high statistical significance, although the amplitudes of the abundance variations are small. Removing abundance trends with effective temperature and/or using a different pair of reference stars does not alter these results. These abundance variations and correlations may reflect a combination of (a) He abundance variations and (b) inhomogeneous chemical evolution in the pre- or proto-cluster environment. Regarding the former, the current constraints on Delta Y from photometry likely preclude He as being the sole explanation. Regarding the latter, the nucleosynthetic source(s) must have synthesised Na, alpha, Fe-peak and neutron-capture elements and in constant amounts for species heavier than Si; no individual object can achieve such nucleosynthesis. We speculate that other, if not all, globular clusters may exhibit comparable abundance variations and correlations to NGC 6752 if subjected to a similarly precise analysis.

preprint2013arXiv

How realistic are solar model atmospheres?

Recently, new solar model atmospheres have been developed to replace classical 1D LTE hydrostatic models and used to for example derive the solar chemical composition. We aim to test various models against key observational constraints. In particular, a 3D model used to derive the solar abundances, a 3D MHD model (with an imposed 10 mT vertical magnetic field), 1D models from the PHOENIX project, the 1D MARCS model, and the 1D semi-empirical model of Holweger & Müller. We confront the models with observational diagnostics of the temperature profile: continuum centre-to-limb variations (CLV), absolute continuum fluxes, and the wings of hydrogen lines. We also test the 3D models for the intensity distribution of the granulation and spectral line shapes. The predictions from the 3D model are in excellent agreement with the continuum CLV observations, performing even better than the Holweger & Müller model (constructed largely to fulfil such observations). The predictions of the 1D theoretical models are worse, given their steeper temperature gradients. For the continuum fluxes, predictions for most models agree well with the observations. No model fits all hydrogen lines perfectly, but again the 3D model comes ahead. The 3D model also reproduces the observed continuum intensity fluctuations and spectral line shapes very well. The excellent agreement of the 3D model with the observables reinforces the view that its temperature structure is realistic. It outperforms the MHD simulation in all diagnostics, implying that recent claims for revised abundances based on MHD modelling are premature. Several weaknesses in the 1D models are exposed. The differences between the PHOENIX LTE and NLTE models are small. We conclude that the 3D hydrodynamical model is superior to any of the tested 1D models, which gives further confidence in the solar abundance analyses based on it.

preprint2013arXiv

The lithium isotopic ratio in very metal-poor stars

Un-evolved, very metal-poor stars are the most important tracers of the cosmic abundance of lithium in the early universe. Combining the standard Big Bang nucleosynthesis model with Galactic production through cosmic ray spallation, these stars at [Fe/H]<-2 are expected to show an undetectably small 6Li/7Li isotopic signature. Evidence to the contrary may necessitate an additional pre-galactic production source or a revision of the standard model of Big Bang nucleosynthesis. We revisit the isotopic analysis of four halo stars, two with claimed 6Li-detections in the literature, to investigate the influence of improved model atmospheres and line formation treatment. For the first time, a combined 3D, NLTE (non-local thermodynamic equilibrium) modelling technique for Li, Na, and Ca lines is utilised to constrain the intrinsic line-broadening and to determine the Li isotopic ratio. We discuss the influence of 3D NLTE effects on line profile shapes and assess the realism of our modelling using the Ca excitation and ionisation balance. By accounting for NLTE line formation in realistic 3D hydrodynamical model atmospheres, we can model the Li resonance line and other neutral lines with a consistency that is superior to LTE, with no need for additional line asymmetry caused by the presence of 6Li. Contrary to the results from 1D and 3D LTE modelling, no star in our sample has a significant 2-sigma detection of the lighter isotope in NLTE. Over a large parameter space, NLTE modelling systematically reduces the best-fit Li isotopic ratios by up to five percentage points. As a bi-product, we also present the first ever 3D NLTE Ca and Na abundances of halo stars, which reveal significant departures from LTE. The observational support for a significant and non-standard 6Li production source in the early universe is substantially weakened by our findings.

preprint2013arXiv

The Stagger-grid - A Grid of 3D Stellar Atmosphere Models

Theoretical atmosphere models provide the basis for a variety of applications in astronomy. In simplified one-dimensional (1D) atmosphere models, convection is usually treated with the mixing length theory despite its well-known insufficiency, therefore, the superadiabatic regime is poorly rendered. Due to the increasing computational power over the last decades, we are now capable to compute large grids of realistic three-dimensional (3D) hydrodynamical model atmospheres with the realistic treatment of the radiative transfer. We have computed the Stagger-grid, a comprehensive grid of 3D atmosphere models for late-type stars. In the presented contribution, we discuss initial results of the grid by exploring global properties and mean stratifications of the 3D models. Furthermore, we also depict the differences to classic 1D atmosphere models.

preprint2013arXiv

The Stagger-grid: A Grid of 3D Stellar Atmosphere Models - II. Horizontal and Temporal Averaging and Spectral Line Formation

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. We explore methods for computing mean <3D> stratifications from the Stagger-grid time-dependent 3D radiative hydro- dynamical atmosphere models by considering four different reference depth scales (geometrical depth, column-mass density, and two optical depth scales). Furthermore, we investigate the influence of alternative averages (logarithmic or enforced hydrostatic equilibrium, flux-weighted temperatures). For the line formation we compute curves of growth for Fe i and Fe ii lines in LTE . The resulting <3D> stratifications for the four reference depth scales can be considerably different. We find typically 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 logg, 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 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 the increasing interference of the reversed granulation towards higher effective temperature, in particular at low metallicity.

preprint2012arXiv

NLTE line formation of Fe for late-type stars. I. Standard stars with 1D and <3D> model atmospheres

We investigate departures from LTE in the line formation of Fe for a number of well-studied late-type stars in different evolutionary stages. A new model of Fe atom was constructed from the most up-to-date theoretical and experimental atomic data available so far. Non-local thermodynamic equilibrium (NLTE) line formation calculations for Fe were performed using 1D hydrostatic MARCS and MAFAGS-OS model atmospheres, as well as the spatial and temporal average stratifications from full 3D hydrodynamical simulations of stellar convection computed using the Stagger code. It is shown that the Fe I/Fe II ionization balance can be well established with the 1D and mean 3D models under NLTE including calibrated inelastic collisions with H I calculated from the Drawin's (1969) formulae. Strong low-excitation Fe I lines are very sensitive to the atmospheric structure; classical 1D models fail to provide consistent excitation balance, particularly so for cool metal-poor stars. A better agreement between Fe I lines spanning a range of excitation potentials is obtained with the mean 3D models. Mean NLTE metallicities determined for the standard stars using the 1D and mean 3D models are fully consistent. Also, the NLTE spectroscopic effective temperatures and gravities from ionization balance agree with that determined by other methods, e.g., infrared flux method and parallaxes, if one of the stellar parameters is constrained independently.

preprint2012arXiv

NLTE line formation of Fe in late-type stars II: 1D spectroscopic stellar parameters

We investigate departures from local thermodynamic equilibrium (NLTE) in the line formation of neutral and singly ionised iron lines and their impact on spectroscopic stellar parameters. The calculations are performed for an extensive grid of 1D MARCS models of metal-rich and metal-poor late-type dwarfs and giants. We find that iron abundances derived from FeI lines are increasingly underestimated in hotter, lower surface-gravity, and more metal-poor stars, in a simple and well-defined pattern, while LTE is usually a realistic approximation for FeII lines. For the vast majority of dwarfs and giants, the perturbed ionisation balance of FeI and FeII is the main relevant NLTE effect to consider in the determination of spectroscopic parameters, while for extremely metal-poor stars and hot giant stars significant impact is seen also on the excitation balance and on the microturbulence determination from FeI lines.

preprint2012arXiv

The detection and treatment of distance errors in kinematic analyses of stars

We present a new method for detecting and correcting systematic errors in the distances to stars when both proper motions and line-of-sight velocities are available. The method, which is applicable for samples of 200 or more stars that have a significant extension on the sky, exploits correlations between the measured U, V and W velocity components that are introduced by distance errors. We deliver a formalism to describe and interpret the specific imprints of distance errors including spurious velocity correlations and shifts of mean motion in a sample. We take into account correlations introduced by measurement errors, Galactic rotation and changes in the orientation of the velocity ellipsoid with position in the Galaxy. Tests on pseudodata show that the method is more robust and sensitive than traditional approaches to this problem. We investigate approaches to characterising the probability distribution of distance errors, in addition to the mean distance error, which is the main theme of the paper. Stars with the most overestimated distances bias our estimate of the overall distance scale, leading to the corrected distances being slightly too small. We give a formula that can be used to correct for this effect. We apply the method to samples of stars from the SEGUE survey, exploring optimal gravity cuts, sample contamination, and correcting the used distance relations.

preprint2012arXiv

The remarkable solar twin HIP 56948: a prime target in the quest for other Earths

We study HIP 56948, the best solar twin known to date, to determine with an unparalleled precision how similar is to the Sun in its physical properties, chemical composition and planet architecture. We explore whether the abundances anomalies may be due to pollution from stellar ejecta or to terrestrial planet formation. We perform a differential abundance analysis (both in LTE and NLTE) using high resolution (R = 100,000) high S/N (600) Keck HIRES spectra of the Sun and HIP 56948. We use precise radial velocity data from the McDonald and Keck observatories to search for planets around this star. We achieve a precision of sigma = 0.003 dex for several elements. Including errors in stellar parameters the total uncertainty is as low as sigma = 0.005 dex (1 %), which is unprecedented in elemental abundance studies. The similarities between HIP 56948 and the Sun are astonishing. HIP 56948 is only 17+/-7 K hotter than the Sun, and log g, [Fe/H] and microturbulence are only +0.02+/-0.02 dex, +0.02+/-0.01 dex and +0.01+/-0.01 km/s higher than solar, respectively. HIP 56948 has a mass of 1.02+/-0.02M_Sun and is 1 Gyr younger than the Sun. Both stars show a chemical abundance pattern that differs from most solar twins. The trend with T_cond in differential abundances (twins - HIP56948) can be reproduced very well by adding 3 M_Earth of a mix of Earth and meteoritic material, to the convection zone of HIP 56948. From our radial velocity monitoring we find no indications of giant planets interior to or within the habitable zone of HIP 56948. We conclude that HIP 56948 is an excellent candidate to host a planetary system like our own, including the possible presence of inner terrestrial planets. Its striking similarity to the Sun and its mature age makes HIP 56948 a prime target in the quest for other Earths and SETI endeavors.

preprint2011arXiv

Is the solar spectrum latitude dependent? An investigation with SST/TRIPPEL

Context: In studies of the solar spectrum relative to spectra of solar twin stars, it has been found that the chemical composition of the Sun seems to depart systematically from those of the twins. One possible explanation is that the effect is due to the special aspect angle of the Sun when observed from Earth, as compared with the aspect angles of the twins. Thus, a latitude dependence of the solar spectrum, even with the heliocentric angle constant, could lead to effects of the type observed. Aim: We explore a possible variation in the strength of certain spectral lines, used in the comparisons between the composition of the Sun and the twins, at loci on the solar disk with different latitudes but at constant heliocentric angle. Methods: We use the TRIPPEL spectrograph at the Swedish 1-m Solar Telescope on La Palma to record spectra in five spectral regions in order to compare different locations on the solar disk at a heliocentric angle of 45 deg. Equivalent widths and other parameters are measured for fifteen different lines representing nine atomic species. Results: The relative variations in equivalent widths at the equator and at solar latitude 45 deg are found to be less than 1.5 % for all spectral lines studied. Translated to elemental abundances as they would be measured from a terrestrial and a hypothetical pole-on observer, the difference is estimated to be within 0.005 dex in all cases. Conclusion: It is very unlikely that latitude effects could cause the reported abundance difference between the Sun and the solar twins. The accuracy obtainable in measurements of small differences in spectral line strengths between different solar disk positions is very high.

preprint2011arXiv

New constraints on the chemical evolution of the solar neighbourhood and Galactic disc(s). Improved astrophysical parameters for the Geneva-Copenhagen Survey

We present a re-analysis of the Geneva-Copenhagen survey, which benefits from the infrared flux method to improve the accuracy of the derived stellar effective temperatures and uses the latter to build a consistent and improved metallicity scale. Metallicities are calibrated on high-resolution spectroscopy and checked against four open clusters and a moving group, showing excellent consistency. The new temperature and metallicity scales provide a better match to theoretical isochrones, which are used for a Bayesian analysis of stellar ages. With respect to previous analyses, our stars are on average 100 K hotter and 0.1 dex more metal rich, which shift the peak of the metallicity distribution function around the solar value. From Stromgren photometry we are able to derive for the first time a proxy for alpha elements, which enables us to perform a tentative dissection of the chemical thin and thick disc. We find evidence for the latter being composed of an old, mildly but systematically alpha-enhanced population that extends to super solar metallicities, in agreement with spectroscopic studies. Our revision offers the largest existing kinematically unbiased sample of the solar neighbourhood that contains full information on kinematics, metallicities, and ages and thus provides better constraints on the physical processes relevant in the build-up of the Milky Way disc, enabling a better understanding of the Sun in a Galactic context.

preprint2011arXiv

NLTE effects on Fe I/II in the atmospheres of FGK stars and application to abundance analysis of their spectra

We describe the first results from our project aimed at large-scale calculations of NLTE abundance corrections for important astrophysical atoms and ions. In this paper, the focus is on Fe which is a proxy of stellar metallicity and is commonly used to derive effective temperature and gravity. We present a small grid of NLTE abundance corrections for Fe I lines and discuss how NLTE effects influence determination of effective temperature, surface gravity, and metallicity for late-type stars.

preprint2011arXiv

On the alleged duality of the Galactic halo

We examine the kinematics of the Galactic halo based on SDSS/SEGUE data by Carollo et al. (2007, 2010). We find that their claims of a counter-rotating halo are the result of substantial biases in distance estimates (of order 50%): the claimed retrograde component, which makes up only a tiny fraction of the entire sample, prone to contaminations, is identified as the tail of distance overestimates. The strong overestimates also result in a lift in the vertical velocity component, which explains the large altitudes those objects were claimed to reach. Errors are worst for the lowest metallicity stars, which explains the metal-poor nature of the artificial component. We also argue that measurement errors were not properly accounted for and that the use of Gaussian fitting on intrinsically non-Gaussian Galactic components invokes the identification of components that are distorted or even artificial. Our evaluation of the data leads to a revision of the estimated velocity ellipsoids and does not yield any reliable evidence for a counterrotating halo component. If a distinct counterrotating halo component exists it must be far weaker than claimed by Carollo et al. Finally we note that their revised analysis presented in Beers et al. (2011) does not alleviate our main concerns.

preprint2011arXiv

The chemical evolution of the Galactic Bulge seen through micro-lensing events

Galactic bulges are central to understanding galaxy formation and evolution. Here we report on recent studies using micro-lensing events to obtain spectra of high resolution and moderately high signal-to-noise ratios of dwarf stars in the Galactic bulge. Normally this is not feasible for the faint turn-off stars in the Galactic bulge, but micro-lensing offers this possibility. Elemental abundance trends in the Galactic bulge as traced by dwarf stars are very similar to those seen for dwarf stars in the solar neighbourhood. We discuss the implications of the ages and metallicity distribution function derived for the micro-lensed dwarf stars in the Galactic bulge.

preprint2011arXiv

The StaggerGrid Project: a Grid of 3D Model Atmospheres for High-Precision Spectroscopy

In this contribution, we present the StaggerGrid, a collaborative project for the construction of a comprehensive grid of time-dependent, three-dimensional (3D), hydrodynamic model atmospheres of solar- and late-type stars with different effective temperatures, surface gravities, and chemical compositions. We illustrate the main characteristics of these 3D models and their effects on the predicted strengths, wavelength-shifts, and shapes of spectral lines, highlighting the differences with respect to calculations based on classical, one-dimensional, hydrostatic models, and discuss some of their possible applications to elemental abundance analysis of stellar spectra in the context of large observational surveys.

preprint2010arXiv

Lithium depletion in solar-like stars: no planet connection

We have determined precise stellar parameters and lithium abundances in a sample of 117 stars with basic properties very similar to the Sun. This sample selection reduces biasing effects and systematic errors in the analysis. We estimate the ages of our sample stars mainly from isochrone fitting but also from measurements of rotation period and X-ray luminosity and test the connection between lithium abundance, age, and stellar parameters. We find strong evidence for increasing lithium depletion with age. Our sample includes 14 stars that are known to host planets and it does not support recent claims that planet-host stars have experienced more lithium depletion than stars without planets. We find the solar lithium abundance normal for a star of its age, mass, and metallicity. Furthermore, we analyze published data for 82 stars that were reported to support an enhanced lithium depletion in planet hosts. We show that those stars in fact follow an age trend very similar to that found with our sample and that the presence of giant planets is not related to low lithium abundances. Finally, we discuss the systematic biases that led to the incorrect conclusion of an enhanced lithium depletion in planet-host stars.

preprint2010arXiv

Observational evidence for a broken Li Spite plateau and mass-dependent Li depletion

We present NLTE Li abundances for 88 stars in the metallicity range -3.5 < [Fe/H] < -1.0. The effective temperatures are based on the infrared flux method with improved E(B-V) values obtained mostly from interstellar NaI D lines. The Li abundances were derived through MARCS models and high-quality UVES+VLT, HIRES+Keck and FIES+NOT spectra, and complemented with reliable equivalent widths from the literature. The less-depleted stars with [Fe/H] < -2.5 and [Fe/H] > -2.5 fall into two well-defined plateaus of A_{Li} = 2.18 (sigma = 0.04) and A_{Li} = 2.27 (sigma = 0.05), respectively. We show that the two plateaus are flat, unlike previous claims for a steep monotonic decrease in Li abundances with decreasing metallicities. At all metallicities we uncover a fine-structure in the Li abundances of Spite plateau stars, which we trace to Li depletion that depends on both metallicity and mass. Models including atomic diffusion and turbulent mixing seem to reproduce the observed Li depletion assuming a primordial Li abundance A_{Li} = 2.64, which agrees well with current predictions (A_{Li} = 2.72) from standard Big Bang nucleosynthesis. Adopting the Kurucz overshooting model atmospheres increases the Li abundance by +0.08 dex to A_{Li} = 2.72, which perfectly agrees with BBN+WMAP.

preprint2010arXiv

The light elements in the light of 3D and non-LTE effects

In this review we discuss possible systematic errors inherent in classical 1D LTE abundance analyses of late-type stars for the light elements (here: H, He, Li, Be and B). The advent of realistic 3D hydrodynamical model atmospheres and the availability of non-LTE line formation codes place the stellar analyses on a much firmer footing and indeed drastically modify the astrophysical interpretations in many cases, especially at low metallicities. For the Teff-sensitive hydrogen lines both stellar granulation and non-LTE are likely important but the combination of the two has not yet been fully explored. A fortuitous near-cancellation of significant but opposite 3D and non-LTE effects leaves the derived 7Li abundances largely unaffected but new atomic collisional data should be taken into account. We also discuss the impact on 3D non-LTE line formation on the estimated lithium isotopic abundances in halo stars in light of recent claims that convective line asymmetries can mimic the presence of 6Li. While Be only have relatively minor non-LTE abundance corrections, B is sensitive even if the latest calculations imply smaller non-LTE effects than previously thought.

preprint2009arXiv

A Puzzle Involving Galactic Bulge Microlensing Events

We study a sample of 16 microlensed Galactic bulge main sequence turnoff region stars for which high dispersion spectra have been obtained with detailed abundance analyses. We demonstrate that there is a very strong and highly statistically significant correlation between the maximum magnification of the microlensed bulge star and the value of the [Fe/H] deduced from the high resolution spectrum of each object. Physics demands that this correlation, assuming it to be real, be the result of some sample bias. We suggest several possible explanations, but are forced to reject them all,and are left puzzled. To obtain a reliable metallicity distribution in the Galactic bulge based on microlensed dwarf stars it will be necessary to resolve this issue through the course of additional observations.

preprint2009arXiv

New Solar Composition: The Problem With Solar Models Revisited

We construct updated solar models with different sets of solar abundances, including the most recent determinations by Asplund et al. (2009). The latter work predicts a larger ($\sim 10%$) solar metallicity compared to previous measurements by the same authors but significantly lower ($\sim 25%$) than the recommended value from a decade ago by Grevesse & Sauval (1998). We compare the results of our models with determinations of the solar structure inferred through helioseismology measurements. The model that uses the most recent solar abundance determinations predicts the base of the solar convective envelope to be located at $R_{\rm CZ}= 0.724{\rm R_\odot}$ and a surface helium mass fraction of $Y_{\rm surf}=0.231$. These results are in conflict with helioseismology data ($R_{\rm CZ}= 0.713\pm0.001{\rm R_\odot}$ and $Y_{\rm surf}=0.2485\pm0.0035$) at 5$-σ$ and 11$-σ$ levels respectively. Using the new solar abundances, we calculate the magnitude by which radiative opacities should be modified in order to restore agreement with helioseismology. We find that a maximum change of $\sim 15%$ at the base of the convective zone is required with a smooth decrease towards the core, where the change needed is $\sim 5%$. The required change at the base of the convective envelope is about half the value estimated previously. We also present the solar neutrino fluxes predicted by the new models. The most important changes brought about by the new solar abundances are the increase by $\sim 10%$ in the predicted $^{13}$N and $^{15}$O fluxes that arise mostly due to the increase in the C and N abundances in the newly determined solar composition.

preprint2009arXiv

Oxygen lines in solar granulation. I. Testing 3D models against new observations with high spatial and spectral resolution

Aims: we seek to provide additional tests of the line formation of theoretical 3D solar photosphere models. In particular, we set out to test the spatially-resolved line formation at several viewing angles, from the solar disk-centre to the limb and focusing on atomic oxygen lines. The purpose of these tests is to provide additional information on whether the 3D model is suitable to derive the solar oxygen abundance. We also aim to empirically constrain the NLTE recipes for neutral hydrogen collisions, using the spatially-resolved observations of the OI 777 nm lines. Methods: using the Swedish 1-m Solar Telescope we obtained high-spatial-resolution observations of five atomic oxygen lines (along with lines for other species) for five positions on the solar disk. These observations have a high spatial and spectral resolution, and a continuum intensity contrast up to 9% at 615 nm. The theoretical line profiles were computed using the 3D model, with a full 3D NLTE treatment for oxygen and LTE for the other lines. Results: at disk-centre we find an excellent agreement between predicted and observed line shifts, strengths, FWHM and asymmetries. At other viewing angles the agreement is also good, but the smaller continuum intensity contrast makes a quantitative comparison harder. We use the disk-centre observations we constrain S_H, the scaling factor for the efficiency of collisions with neutral hydrogen. We find that S_H=1 provides the best match to the observations. Conclusions: overall there is a very good agreement between predicted and observed line properties over the solar granulation. This further reinforces the view that the 3D model is realistic and a reliable tool to derive the solar oxygen abundance.

preprint2009arXiv

Oxygen lines in solar granulation. II. Centre-to-limb variation, NLTE line formation, blends and the solar oxygen abundance

Context: There is a lively debate about the solar oxygen abundance and the role of 3D models in its recent downward revision. The models have been tested using high resolution solar atlases. Further testing can be done using centre-to-limb variations. Aims: Using high quality observations of oxygen lines across the solar surface we seek to test if the 3D and 1D models reproduce their observed centre-to-limb variation (CLV). In particular we seek to assess whether the 3D model is appropriate to derive the solar oxygen abundance. Methods: We use our recent observations of OI 777 nm, OI 615.81 nm, [OI] 630.03 nm and nine lines of other elements for five viewing angles 0.2<mu<1 of the quiet solar disk. We compare them with the predicted line profiles from the 3D and 1D models computed using the most up-to-date line formation codes, line data and allowing for departures of LTE. The CLV of the OI 777 nm lines is also used to obtain an empirical correction for the poorly known efficiency of the inelastic collisions with H. Results: The 3D model generally reproduces the CLV observations of the lines very well, particularly the oxygen lines. From the OI 777 nm lines we find that the classical Drawin recipe slightly overestimates H collisions. The limb observations of the OI 615.82 nm line allow us to identify a previously unknown contribution of molecules for this line, prevalent at the solar limb. A detailed treatment of the [OI] 630.03 nm line shows that the 3D modeling provides an excellent agreement with the observations. The derived oxygen abundances with the 3D model are 8.68 (OI 777 nm), 8.66 ([OI] 630.03 nm) and 8.62 (OI 615.82 nm). Conclusions: These additional tests have reinforced the trustworthiness of the 3D model and line formation for abundance analyses.

preprint2009arXiv

The chemical composition of the Sun

The solar chemical composition is an important ingredient in our understanding of the formation, structure and evolution of both the Sun and our solar system. Furthermore, it is an essential reference standard against which the elemental contents of other astronomical objects are compared. In this review we evaluate the current understanding of the solar photospheric composition. In particular, we present a re-determination of the abundances of nearly all available elements, using a realistic new 3-dimensional (3D), time-dependent hydrodynamical model of the solar atmosphere. We have carefully considered the atomic input data and selection of spectral lines, and accounted for departures from LTE whenever possible. The end result is a comprehensive and homogeneous compilation of the solar elemental abundances. Particularly noteworthy findings are significantly lower abundances of carbon, nitrogen, oxygen and neon compared with the widely-used values of a decade ago. The new solar chemical composition is supported by a high degree of internal consistency between available abundance indicators, and by agreement with values obtained in the solar neighborhood and from the most pristine meteorites. There is, however, a stark conflict with standard models of the solar interior according to helioseismology, a discrepancy that has yet to find a satisfactory resolution.

preprint2009arXiv

The peculiar solar composition and its possible relation to planet formation

We have conducted a differential elemental abundance analysis of unprecedented accuracy (0.01 dex) of the Sun relative to 11 solar twins from the Hipparcos catalogue and 10 solar analogs from planet searches. We find that the Sun shows a characteristic signature with a ~20% depletion of refractory elements relative to the volatile elements in comparison with the solar twins. The abundance differences correlate strongly with the condensation temperatures of the elements. This peculiarity also holds in comparisons with solar analogs known to have close-in giant planets while the majority of solar analogs found not to have such giant planets in radial velocity monitoring show the solar abundance pattern. We discuss various explanations for this peculiarity, including the possibility that the differences in abundance patterns are related to the formation of planetary systems like our own, in particular to the existence of terrestrial planets.

preprint2007arXiv

Very Large Excesses of 18O in Hydrogen-Deficient Carbon and R Coronae Borealis Stars: Evidence for White Dwarf Mergers

We have found that at least seven hydrogen-deficient carbon (HdC) and R Coronae Borealis (RCB) stars, have 16O/18O ratios close to and in some cases less than unity, values that are orders of magnitude lower than measured in other stars (the Solar value is 500). Greatly enhanced 18O is evident in every HdC and RCB we have measured that is cool enough to have detectable CO bands. The three HdC stars measured have 16O/18O < 1, lower values than any of the RCB stars. These discoveries are important clues in determining the evolutionary pathways of HdC and RCB stars, for which two models have been proposed: the double degenerate (white dwarf (WD) merger), and the final helium-shell flash (FF). No overproduction of 18O is expected in the FF scenario. We have quantitatively explored the idea that HdC and RCB stars originate in the mergers of CO- and He-WDs. The merger process is estimated to take only a few days, with accretion rates of 150 Msun/ yr producing temperatures at the base of the accreted envelope of 1.2 - 1.9 x 10^8 K. Analysis of a simplified one-zone calculation shows that nucleosynthesis in the dynamically accreting material may provide a suitable environment for a significant production of 18O, leading to very low values of 16O/18O, similar to those observed. We also find qualitative agreement with observed values of 12C/13C and with the CNO elemental ratios. H-admixture during the accretion process from the small H-rich C/O WD envelope may play an important role in producing the observed abundances. Overall our analysis shows that WD mergers may very well be the progenitors of O18-rich RCB and HdC stars, and that more detailed simulations and modeling are justified.

preprint2005arXiv

An Extremely Large Excess of 18O in the Hydrogen-Deficient Carbon Star, HD 137613

We report the discovery of a uniquely large excess of 18O in the hydrogen-deficient carbon (HdC) star, HD 137613, based on a spectrum of the first overtone bands of CO at 2.3-2.4 micron in which three strong absorption bands of 12C18O are clearly present. Bands of 12C16O also are present but no bands of 13C16O or 12C17O are seen. We estimate an isotopic ratio 16O/18O ~1. The Solar value of this ratio is ~500. Neither He-core burning nor He-shell flash burning can produce the isotopic ratios of oxygen and carbon observed in HD 137613. However, a remarkable similarity exists between the observed abundances and those found in the outer layers of the broad He-shell of early-AGB stars, soon after the end of He-core burning. It is not known how the outer envelope down to the He-shell could be lost but some mechanism of enhanced mass loss must be involved. HD 137613 may be a post-early-AGB star with the outer layers of the former He-burning shell as its photosphere. The unusual elemental abundances of the HdC stars resemble those of the R Coronae Borealis (RCB) stars, but HdC stars do not produce clouds of dust that produce declines in brightness. None of the other RCB or HdC stars observed shows significant 18O.

preprint2005arXiv

Effects of line-blocking on the non-LTE Fe I spectral line formation

The effects of background line opacity (line-blocking) in statistical equilibrium calculations for Fe in late-type stellar atmospheres have been investigated using an extensive and up-to-date model atom with radiative data primarily from the IRON Project. The background metal line opacities have been computed using data from the MARCS stellar model atmospheres. While accounting for this line opacity is important at solar metallicity, the differences between calculations including and excluding line-blocking at low metallicity are insignificant for the non-local thermodynamic equilibrium (non-LTE) abundance corrections for Fe I lines. The line-blocking has no impact on the non-LTE effects of Fe II lines. The dominant uncertainty in Fe non-LTE calculations for metal-poor stars is still the treatment of the inelastic H I collisions, which have here been included using scaling factors to the classical Drawin formalism, and whether or not thermalisation of the high Fe I levels to Fe II ground state should be enforced. Without such thermalisation, the Fe I non-LTE abundance corrections are substantial in metal-poor stars: about 0.3 dex with efficient (i.e. Drawin-like) H I collisions and about 0.5 dex without. Without both thermalisation and H I collisions, even Fe II lines show significant non-LTE effects in such stars.

preprint2005arXiv

Non-LTE neutral carbon spectral line formation in late-type stars

We present non-Local Thermodynamic Equilibrium (non-LTE) calculations for neutral carbon spectral line formation, carried out for a grid of model atmospheres covering the range of late-type stars. The results of our detailed calculations suggest that the carbon non-LTE corrections in these stars are higher than usually adopted, remaining substantial even at low metallicity. For the most metal-poor stars in the sample of Akerman et al. (2004), the non-LTE abundance corrections are of the order of -0.35...-0.45 dex (when neglecting H collisions). Applying our results to those observations, the apparent [C/O] upturn seen in their LTE analysis is no longer present, thus revealing no need to invoke contributions from Pop. III stars to the carbon nucleosynthesis.

preprint2004arXiv

The solar chemical composition

We review our current knowledge of the solar chemical composition as determined from photospheric absorption lines. In particular we describe the recent significant revisions of the solar abundances as a result of the application of a time-dependent, 3D hydrodynamical model of the solar atmosphere instead of 1D hydrostatic models. This has decreased the metal content in the solar convection zone by almost a factor of two compared with the widely used compilation by Anders & Grevesse (1989). While resolving a number of long-standings problems, the new 3D-based element abundances also pose serious challenges, most notably for helioseismology.

preprint2001arXiv

Spectroscopic Signatures of Convection in the Spectrum of Procyon. Fundamental Parameters and Iron Abundance

We have observed the spectrum of Procyon A (F5IV) from 4559 to 5780 A with a S/N of ~ 1e3 and a resolving power of 2e5. We have measured the line bisectors and relative line shifts of a large number of Fe I and Fe II lines, comparing them to those found in the Sun. A three-dimensional(3D) hydrodynamical model atmosphere has been computed and is tested against observations. The model reproduces in detail most of the features observed, although we identify some room for improvement. At all levels, the comparison of the 3D time-dependent calculations with the observed spectral lines shows a much better agreement than for classical homogeneous models, making it possible to refine previous estimates of the iron abundance, the projected rotational velocity, the limb-darkening, and the systemic velocity of the Procyon binary system. The difference between the iron abundance determined with the 3D model and its 1D counterpart is <~ 0.05 dex. We find consistency between the iron abundance derived from Fe I and Fe II lines, suggesting that departures from LTE in the formation of the studied lines are relatively small. The scatter in the iron abundance determined from different lines still exceeds the expectations from the uncertainties in the atomic data, pointing out that one or more components in the modeling can be refined further.

preprint2001arXiv

The 'Forbidden' Abundance of Oxygen in the Sun

We reexamine closely the solar photospheric line at 6300 A, which is attributed to a forbidden line of neutral oxygen, and is widely used in analyses of other late-type stars. We use a three-dimensional time-dependent hydrodynamical model solar atmosphere which has been tested successfully against observed granulation patterns and an array of absorption lines. We show that the solar line is a blend with a Ni I line, as previously suggested but oftentimes neglected. Thanks to accurate atomic data on the [O I] and Ni I lines we are able to derive an accurate oxygen abundance for the Sun: log epsilon (O) = 8.69 +/- 0.05 dex, a value at the lower end of the distribution of previously published abundances, but in good agreement with estimates for the local interstellar medium and hot stars in the solar neighborhood. We conclude by discussing the implication of the Ni I blend on oxygen abundances derived from the [O I] 6300 A line in disk and halo stars.