Lithium abundances in globular clusters

Abstract: Lithium is created during the Big Bang nucleosynthesis and it is destroyed in stellar interiors at relatively low temperatures. However, it should be preserved in the stellar envelopes of unevolved stars and progressively diluted during mixing processes. In particular, after the first dredge-up along the RGB, lithium should be completely destroyed, but this is not what we observe today in globular clusters. This element allows to test stellar evolutionary models, as well as different types of polluters for second population stars in the multiple population scenarios. Due to the difficulty in the measurement of the small available lithium line, few GCs have been studied in details so far. Literature results are not homogeneous for what concerns type of stars, sample sizes, and chemical analysis methods. The Gaia-ESO survey allows us to study the largest sample of GCs stars (about 2000, both dwarfs and giants) for which the lithium has been analysed homogeneously.

\section{Introduction}

Lithium (Li) is one of the elements created during the Big Bang nucleosynthesis
and it is destroyed in the stellar interiors at relatively low temperatures
($\sim2.5$x$10^6K$). However, it should be preserved in the stellar envelopes
of unevolved stars and diluted by diffusion and mixing processes.

\begin{table*}[htbp!]
\caption{\footnotesize GCs included in GES, with the main
evolutionary phases covered by the members sample.}
\label{table:gc}
\begin{center}
\begin{tabular}{|ll|ll|ll|}
\hline \hline
Cluster&Phase&Cluster&Phase&Cluster&Phase\\
\hline
M12&RGB, AGB&M15&RGB, AGB&M2&RGB, AGB\\
NGC104&MS, RGB, AGB&NGC362&RGB, AGB&NGC1261&RGB, AGB\\
NGC1851&RGB, AGB&NGC1904&RGB, AGB&NGC2808&RGB, AGB\\
NGC4372&RGB, AGB&NGC4590&RGB, AGB&NGC4833&RGB, AGB\\
NGC5927&RGB, AGB&NGC6553&RGB, AGB&NGC6752&MS, RGB, AGB\\
\hline
\end{tabular}
\end{center}
\end{table*}

The study of Li abundance in Globular Clusters (hereafter
GCs), among the oldest objects in the Galaxy, allows to investigate the cosmic
lithium problem, as well as to  give important information on
stellar evolution and mixing processes and on the formation of
multiple populations. Few GCs have been investigated so far in detail, and
different evolutionary phases have  been considered, but the emeriging scenario
is complex. In general, the dwarf stars share the same lithium abundance,
$A(Li)\sim2.2$ dex as found by \citet{spitea, spiteb}, forming
the so called Spite plateau.  During the sub-giant branch phase
diffusion processes dilute Li, as well as
mixing during the first dredge-up process in the red giant
branch (hereafter RGB). However, \citet{mucciarelli12} found that stars above
the first dredge-up and below the RGB bump form another plateau
at about $A(Li)\simeq1.0$ dex, while after the bump Li is
diluted again. The above framework is further complicated by
the phenomenon of multiple populations in globular clusters. In fact, in some
clusters (e. g., NGC~362) first and second population stars
have the same Li abundance, while in others (e. g., NGC~6752)
they do not, and in the case of $\omega$ Cen, both behaviors are observed at the
same time \citep[see][and reference therein, for details]{mucciarelli18}.

Moreover, 14 Li-rich stars have been discovered so far in 12
different GCs \citep{mucciarelli19}. The majority are RGB
stars, with only two main sequence (hereafter MS) stars and two
that belong to the asymptotic giant branch (hereafter AGB). Different
enrichment scenarios have been proposed to explain their
existance, such as engulfment of sub-stellar systems,
self-enrichment based on the Cameron-Fowler mechanism \citep{cameronfowler} and
mass-transfer, but a consensus has not been reached yet.

The existence of stars enhanced in lithium has also important
implications in terms of the formation of multiple populations in GCs. Being so
fragile, Li is a powerful tracer of the physical conditions of the polluting
source, and poses strict limits to the temperatures involved in the process.

\section{Data sets and analysis}

This work is based on spectra gathered by the Gaia-ESO survey
\citep[hereafter GES]{gilmore12,randich13} and on its recommended parameters and
abundances. The spectra have been acquired using FLAMES-GIRAFFE
\citep{pasquini} and FLAMES-UVES \citep{dekker00}, combining new observations
and ESO archival data obtained with the GES instrumental
setups. The reduction of the GIRAFFE and UVES spectra was performed as
described in \citet{jackson15} and \citet{smiljanic}, respectively. The
atmospheric parameters and  abundances were determined with the GES
multiple-pipeline strategy
\citep{recioblanco14,smiljanic,pancino17}. The lithium
abundance has been derived using the Li 6707 \AA~doublet, taking into account
the blend with the neighboring Fe line. Within GES, 15 GCs were observed as part
of the survey calibration strategy \citep{pancino17,pancino17a}, covering a
range of metallicities and global properties. We list them in
Table~\ref{table:gc} together with the evolutionary stages sampled by the member
stars. In particular, RGB stars have been observed
in all GCs, while for two GCs MS stars from the archive sample
have been analyzed as well: this represents the largest sample homogeneously
analyzed ever used to spectroscopically investigate Lithium in GCs.
Of the about 2500 stars analyzed in GES, we selected almost
2000 bona-fide members using the GES [Fe/H] and radial velocity determinations,
following \citet{pancino17a}.

\begin{figure*}[htbp!]
\resizebox{\hsize}{!}{\includegraphics[clip=true]{li_plot.eps}}
\caption{\footnotesize Li abundance as a function of surface
gravity for the member stars of the 15 {GCs included in GES},
color-coded by metallicity.  The Spite plateau is well defined, at
{A(Li)~$\simeq 2.2$ dex}, as well as the plateau
defined by stars above the first dredge-up and below the bump
along the RGB, at A(Li)~$\simeq 1.1$ dex. The evolutionary
phases of the stars included in the sample are clearly identified.}
\label{li}
\end{figure*}

\section{Discussion}

Our knowledge about Li in GCs is based on inhomogeneous
studies, focused on a few clusters, investigated with varying sample sizes, and
analyzed with different methods. For example, in the case of NGC~6752 only 9 MS
stars have been investigated \citep{pasquini05}, while $\simeq 350$ stars (from
below the MS turnoff up to the bright end of the RGB) have been studied in the
case of NGC~6397 \citep{lind09}.

GES provides an unprecedented data set, that allows us to
effectively attack the lithium problem in GCs, thanks to the very high quality
of the spectra and the strategy used to homogenize the abundances.
Typical uncertainties in the abundances of metal-poor stars
amount to about 0.1~dex or less \citep{pancino17a}. In this respect, we can also
take advantage of the relatively high number of GCs included in the survey and
on their different properties, covering a range of global GC parameters like
metallicity, total mass, concentration, and so on.

Figure~\ref{li} shows the Li abundances of the $\simeq$2000
selected member stars in the 15 GCs, determined by GES, as a function of surface
gravity (log$g$). The Spite plateau defined by dwarf stars
appears well defined at A(Li)~$\simeq 2.2$ dex, as well as
the plateau defined by stars above the first dedge-up and below
the RGB bump, at A(Li)~$\simeq 1.1$ dex. The RGB stars brighter than the
RGB bump show the expected behaviour, whit Lithium
progressively diluted due to mixing processes. The stars
in Figure~\ref{li} are color-coded by metallicity
showing several effects that are worth further investigation.

\subsection{Li-rich stars}

As already mentioned above, 14 Li-rich stars in GCs are known
in the literature, and their origin mechanism is still not well
understood. Two of them are included in the GES sample
presented here: the RGB bump star discovered in  NGC~362 by
\citet{dorazi} and the bright RGB or AGB star discovered in
NGC~4590 by \citet{ruchti} and \citet{kirby}.

Here we take the occasion and use these stars to demonstrate
the exquisite quality of the GES spectra.  Figure~\ref{li_rich} shows the
comparison between the {spectra of the two Li-rich stars in
NGC~362 and NGC~4590 observed by GES, compared with Li-normal member stars of
the two GCs, having similar atmospheric parameters. The typical Li enhancement
with respect to similar cluster stars is of at least one order of magnitude. The
two Li-rich stars are well visible also in Figure~\ref{li}, where they display
A(Li)$\simeq$2.7~dex in the case of NGC~362 and A(Li)$\simeq$3.3~dex in the case
of NGC~4590. These preliminary abundance estimates are not yet corrected for
NLTE effects.}

\begin{figure}[]
\resizebox{\hsize}{!}{\includegraphics[clip=true]{lirich.eps}}
\caption{\footnotesize Top panel: comparison between the
spectrum of the Li-rich star in NGC~362 (red line) and a reference member star
with similar atmospheric parameters. Bottom panel: a similar comparison in the
case of the Li-rich star in NGC~4590.}
\label{li_rich}
\end{figure}

\section{Conclusions}

Using prelimnary abundances from the Gaia-ESO Survey, we have
demonstrated the power of the GES dataset for a global understanding of lithium
in GCs. Almost 2000 member stars at different evolutionary stages were analyzed
homogeneously, belonging to 15 GCs with different metallicity, mass,
concentration, orbit, and selected to cover well the space of GC global
properties as calibrators for GES. This sample offers the possibility  to
investigate stars with different properties, such as temperature, gravity and
metallicity, in the largest sample homogeneously analyzed so far. The exquisite
quality of the spectra allows us to well determine the Li abundance and identify
rare Li-rich stars in this kind of systems.
The GES dataset will enable a full statistical investigation of
the Li behaviour in GCs, and it will undoubtedly provide an invaluable reference for further
theoretical studies, allowing for a deeper understanding of stellar evolution
and providing hopefully new insights on the multiple stellar population problem
in GCs.

\bibliographystyle{aa}

Solar and Stellar Astrophysics (astro-ph.SR); Astrophysics of Galaxies (astro-ph.GA)

Lithium abundances in globular clusters


Initial Chemical Enrichment in Galaxies

Abstract: We present evidence that damped Lyman-alpha galaxies detected in spectra of quasars may not have started forming stars until the redshift z~3. If damped Lyman-alpha absorbers are the progenitors of disk galaxies, then the above result may indicate that star formation in galactic disks first began at z~3.

\section{Introduction}

  One definition  of the epoch of galaxy formation is when galaxies first began
to form stars. Different types of galaxies (e.g., ellipticals and spirals)
or different parts of the same type of galaxies (e.g., bulge and
disk of spirals) may have formed at different epochs. A important goal of
observational cosmology is to identify these different formation
epochs.

  Conventional wisdom suggests that ellipticals and the spheroidal
component of  spirals
formed very early, followed by disks. The population of galaxies
at $z>3$ identified using the Lyman limit drop-out
technique may very well be the progenitors of the spheroidal component
of massive galaxies [1,2].
Here we discuss evidence for a sharp rise in the metallicity
distribution of damped Ly$\alpha$ absorbers at $z\leq 3$,
which may signify the onset of star formation in galactic disks.

\section{Results}

   Damped Ly$\alpha$ (DLA) absorption systems seen in spectra of background
quasars are widely accepted to be the progenitors of present-day galaxies
[3], although their exact nature (dwarfs, spheroids, or disks?) is still
unclear.  A program is carried out using the Keck telescopes to study the
chemical compositions of the absorbing gas in DLA systems.
One of the goals is to (hopefully) identify the epoch of the first episode
of star formation in these galaxies,
hence constraining theories  of galaxy formation.

   Figure 1 shows the distribution of [Fe/H] in DLA systems as a function
of redshift. Detailed descriptions of the data and analyses are given in
refs [4,5]. The low
metallicities of DLAs testify the youth of these galaxies: they have yet to
make the bulk of their stars. Remarkably, all 6 of the
highest redshift absorbers
have [Fe/H]$\leq -2$; while many absorbers have reached ten times higher
metallicity at just slightly lower redshifts. This indicates an epoch of rapid
star formation at $z\sim 3$. The effect is likely to be real:
if DLA systems at $z>3$ have [Fe/H] that is uniformly distributed between
$-1$ and $-2.5$ (i.e.,, similar to the distribution at $2<z<3$),
then the {\it posterior} probability for all six of the highest
redshift systems to have [Fe/H]$\leq -2$ by chance is $1.4\times 10^{-3}$.

\begin{figure}
\psfig{figure=figure1.ps,height=2.5in,width=3.5in}
\caption{Metallicity distribution of damped Ly$\alpha$ absorbers.
   \label{fig:figure 1}}
\end{figure}

   Coincidentally, the metallicities of DLA systems at $z>3$, [Fe/H]$=-2$ to
$-2.5$, is identical (within the uncertainties) to those found for the
IGM clouds at similar redshifts, as inferred from the C IV absorption
associated with Ly$\alpha$ forest clouds [6,7,8].  This coincidence
suggests that the metals in DLA galaxies at $z>3$ may simply reflect
those in the IGM, however they were made (e.g., Pop III stars, ejected from
protogalaxies); DLA galaxies did not start making
their own stars (hence metals) until $z\sim 3$.

\section{Discussion}

The implications of the above result for the general question of
galaxy formation and evolution depend on the nature of the DLA galaxies.

   It was suggested [9] that DLA systems may represent the progenitors
of disk galaxies. This is supported by the very recent finding [10] that
the kinematics of DLA absorbers as inferred from the metal absorption line
profiles appears to be dominated by rotation with
large circular velocities ($>200$ km s$^{-1}$).
However, the mean metallicities of DLAs at $z>1.6$ are significantly
below that of the Milky Way disk at the corresponding epoch [4,11].
The problem with the metallicity
distribution may be lessened if DLAs represent a thick disk phase of galaxies
[3] or if low surface brightness disk galaxies (which have substantially
sub-solar metallicities) make up a significant
fraction of DLA absorbers [12]. {\it If} the
disk hypothesis for DLA absorbers is correct, we may have identified the
epoch of initial star formation in disk galaxies.

   Alternatively, DLAs may represent dwarf galaxies or the spheroidal
component of massive
galaxies; this conjecture stems from  the similarity
between the metallicity distribution of DLAs and those in
halo globular clusters and local gas-rich dwarf galaxies [4].
In this case, however, one has to explain the kinematics of DLAs
[10] by other means.

\section*{Acknowledgments}
LL appreciates support from a Hubble Fellowship (HF1062-01-94A). WWS
was supported by NSF grant AST95-29073.

\section*{References}

Initial Chemical Enrichment in Galaxies


Binarity as the Origin of Long Secondary Periods in Red Giant Stars

Abstract: Long secondary periods (LSPs), observed in a third of pulsating red giant stars, are the only unexplained type of large-amplitude stellar variability known at this time. Here we show that this phenomenon is a manifestation of a substellar or stellar companion orbiting the red giant star. Our investigation is based on a sample of about 16,000 well-defined LSP variables detected in the long-term OGLE photometric database of the Milky Way and Magellanic Clouds, combined with the mid-infrared data extracted from the NEOWISE-R archive. From this collection, we selected about 700 objects with stable, large-amplitude, well-sampled infrared light curves and found that about half of them exhibit secondary eclipses, thus presenting an important piece of evidence that the physical mechanism responsible for LSPs is binarity. Namely, the LSP light changes are due to the presence of a dusty cloud orbiting the red giant together with the companion and obscuring the star once per orbit. The secondary eclipses, visible only in the infrared wavelength, occur when the cloud is hidden behind the giant. In this scenario, the low-mass companion is a former planet that has accreted a significant amount of mass from the envelope of its host star and grown into a brown dwarf.

\section{Introduction}

As a red giant star evolves along the asymptotic giant branch (AGB) and it exceeds a certain luminosity, it becomes a variable star. The photometric variability of red giants manifests itself in several ways. Single or (more frequently) multiperiodic brightness changes are caused by stellar pulsations, irregular variations signal the presence of circumstellar dust, while stochastic noise in the light curves may result from the supergranular convection in the envelopes of giant stars. However, there is one type of red giant variability that has not yet found a satisfactory explanation -- the so-called long secondary periods (LSPs).

The LSPs range from several months to several years -- an order of magnitude longer than the pulsation periods. This phenomenon can be detected in at least a third of luminous AGB stars and supergiants. It is likely that low-amplitude LSP variations occur also in the brightest first-ascent red giant branch (RGB) stars\footnote{In this paper we use the term ``red giant'' for both, AGB and RGB, stars} \citep{pawlak2021}. The visual amplitudes of the LSP variations reach 1~mag, while the peak-to-peak radial velocity amplitudes range from 2 to 7~km~s$^{-1}$, with a clustering around 3.5~km~s$^{-1}$ \citep{hinkle2002,wood2004,nicholls2009}. In the period--luminosity plane, LSP variables obey a relation called sequence~D \citep{wood1999}, which partly overlaps with sequence~E populated by red giants in close binary systems showing ellipsoidal variability \citep{soszynski2004,nicholls2010}.

It is now beyond doubt that the presence of the LSPs in red giants is associated with unusual amounts of circumstellar dust. \citet{wood2009} found a mid-infrared (mid-IR) excess in the LSP variables compared to stars without LSPs, which indicates the presence of circumstellar matter absorbing the energy in the optical light and reemitting it at longer wavelengths. Moreover, \citet{wood2009} noticed that dust around LSP variables is not spherically symmetric distributed, but it is rather in a clumpy or a disk-like configuration. \citet{mcdonald2019} linked the onset of the LSP phenomenon in asymptotic giant branch (AGB) stars with the onset of strong mass-loss due to dust-driven winds. \citet{pawlak2021} analyzed spectral energy distributions of a large number of LSP and non-LSP variables and confirmed these conclusions.

Numerous authors have explored various scenarios for the origin of LSPs \citep[e.g.][]{hinkle2002,olivier2003,wood2004,derekas2006,nicholls2009,nie2010,stothers2010,takayama2015,percy2016}, but were unable to give a final solution to this problem. In the recent years, two the most favored explanations have been discussed in the literature: oscillatory convective non-radial modes \citep[e.g.][]{saio2015,takayama2020} and eclipses by a dusty cloud surrounding and following a close-orbiting brown-dwarf or low-mass stellar companion. The latter hypothesis was proposed by \citet{wood1999} and developed by \citet{soszynski2007} and \citet{soszynski2014}. In this scenario, the presence of the brown-dwarf companion is necessary to produce the observed radial velocity amplitudes, but the large fraction of the LSP variables among red giant stars can only be explained by the fact that the companion was initially an exoplanet that accreted a significant amount of mass from the envelope of its host.

\begin{figure}
\centering
\includegraphics[width=8.3cm]{fig1.eps}
\caption{Example OGLE {\it I}-band light curves of LSP variables from the Galactic bulge and LMC. Tick marks on the vertical axis are spaced at 0.1~mag intervals.\label{fig1}}
\end{figure}

There have been surprisingly few papers that addressed the characteristic light curve shapes of LSP variables. The LSP light curves have significantly different morphology than the light curves produced by, both, radial and non-radial pulsation modes in red giants. The sample of {\it I}-band light curves of fifteen LSP variables from the Galactic bulge and Large Magellanic Cloud (LMC) are displayed in Figure~\ref{fig1}. The light curves are folded with the LSPs and sorted by their amplitudes. Here, we selected objects with relatively small amplitudes of the red giant pulsations (visible as short-period oscillations in the light curves) to highlight the characteristic features of the LSP modulation.

A typical LSP cycle can be divided into two parts, each lasting approximately half of the period. During the first part, the brightness does not vary at all (except pulsations) or it increases slowly with time. The second part of the light curve has a triangular shape -- the luminosity decreases, reaches a minimum, and then increases. In the largest-amplitude LSP variables, the duration of this triangle minimum may extend and exceed half the cycle, while the maxima often take a rounded shape. In some stars, the minima of the light curves significantly vary in depth from cycle to cycle, sometimes the LSP modulation completely disappears for a period of time.

\citet{soszynski2014} demonstrated that the LSP light curves can be produced by a dusty cloud with a comet-like tail trailing behind a low-mass companion. It is worth noting that although the radial velocity amplitudes in the vast majority of LSP variables are consistent with the presence of a brown-dwarf companion, the characteristic LSP light curves can also be observed in stellar binary systems, for example in some symbiotic stars \citep[e.g.][]{gromadzki2013}. Additionally, \citet{soszynski2014} showed that the non-sinusoidal shapes of the radial velocity curves \citep[e.g.][]{hinkle2009,nicholls2009} can be explained by the Rossiter-McLaughlin effect, even if a companion is in a circular orbit. On the other hand, the stellar pulsation models are capable of obtaining periods \citep[e.g.][]{saio2015} or amplitudes \citep[e.g.][]{takayama2020} of the LSP modulation, but are unable to reproduce the characteristic light curve shapes shown in Figure~\ref{fig1}.

In this paper, we study the mid-infrared (mid-IR) light curves of LSP giants detected in the Optical Gravitational Lensing Experiment (OGLE) photometric databases. The binary model of the LSP phenomenon predicts an additional effect that should be visible in the mid-IR regime -- the secondary eclipses that occur when the dusty cloud hides behind the red giant. The discovery of the secondary eclipses in the LSP light curves would be a strong argument supporting the binary explanation of the LSPs in red giants stars.

\section{Sample Selection}

The OGLE project \citep{udalski2015} is an ongoing sky survey conducting photometric observations of the Magellanic Clouds and the Galactic bulge and disk, about 3500 square degrees in total. The observations are carried out with the 1.3-meter Warsaw Telescope at Las Campanas Observatory, Chile. The OGLE database contains long-term light curves for about 2~billion stars obtained in the {\it I} and {\it V}-band passbands from the Johnson-Cousins photometric system. This dataset is ideally suited to search for variable red giants and study their photometric properties.

Based on observations obtained between 1997 and 2009, during the second and third phases of the OGLE project (OGLE-II and OGLE-III), the OGLE team found over 340 thousand long-period variables (pulsating red giants) in the Magellanic System and Galactic bulge \citep{soszynski2009,soszynski2011,soszynski2013}. At least 100 thousand of these objects should be potential LSP variables, however for the vast majority of them the photometric amplitudes associated with the LSP changes are smaller than the amplitudes caused by the stellar pulsation. To avoid confusing the LSP modulation with other types of red giant variability (e.g. ellipsoidal variations or irregular dimming events due to dust obscuration), we retrieved from the OGLE catalogs only relatively large-amplitude, characteristic LSP light curves, similar to those shown in Figure~\ref{fig1}. Additionally, we carried out a limited search for LSP variables in the region of the sky not observed during the OGLE-III survey, but monitored during the fourth phase of the OGLE project (OGLE-IV) from 2010 to 2020.

\begin{figure*}
\centering
\includegraphics[width=17.7cm]{fig2.eps}
\caption{Period--Wesenheit index diagrams for long period variables in the LMC (left panel) and SMC (right panel). The NIR Wesenheit index is an extinction-free quantity defined as $W_{JK}=K_\mathrm{s}-0.686(J-K_\mathrm{s})$, where {\it J}- and $K_\mathrm{s}$-band magnitudes originate from the IRSF Catalog \citep{kato2007} or from the 2MASS Catalog \citep{cutri2003}. The period--luminosity sequences are labeled with letters according to the \citet{wood1999} nomenclature. Different colors of the points refer to different types of variable red giants: orange -- OGLE small-amplitude red giants (OSARGs; sequences A$'$, A, and B), red -- semiregular variables (SRVs; sequences C$'$ and C), brown -- Miras (sequence C), blue -- LSP variables (sequence D), and green -- eclipsing and ellipsoidal close binary systems containing red giants as one of the components (sequence E). Periods and classification of the stars were taken from the OGLE Collection of Variable Stars \citep{soszynski2009,soszynski2011}. \label{fig2}}
\end{figure*}

As a result, we found about 16,000 well defined LSP variables: $\sim$9000 in the Milky Way and $\sim$7000 in the Magellanic Clouds. The OGLE-IV {\it I}-band light curves of these stars contain from a few hundred to more than 10,000 data points, depending on the field. Then, we fine-tuned the LSPs of our sample using the {\sc Tatry} code which implements the multiharmonic analysis of variance algorithm \citep{schwarzenberg1996}. The classification of the LMC and SMC variables was additionally verified by the inspection of the near-infrared (NIR) period--luminosity diagrams for both galaxies, namely $\log{P}$--$W_{JK}$ diagram (Figure~\ref{fig2}), where $W_{JK}=K_\mathrm{s}-0.686(J-K_\mathrm{s})$ is the extinction-free Wesenheit index, while {\it J}- and $K_\mathrm{s}$-band magnitudes were taken from the IRSF Magellanic Clouds Point Source Catalog \citep{kato2007} or (if not available) from the 2MASS All-Sky Catalog of Point Sources \citep{cutri2003}. Objects located outside sequence~D \citep{wood1999} in the period--luminosity diagram were removed from our list of LSP stars. The vast majority of the selected LSP variables are undoubtedly AGB stars, since no significant surplus of objects can be detected at luminosities below the RGB tip. In the next step, we extracted the mid-IR light curves of our sample.

\section{Mid-IR Photometry}

The Wide-field Infrared Survey Explorer \citep[WISE;][]{wright2010} is a 40-cm telescope onboard a satellite in a Sun-synchronous polar orbit around Earth. In 2010, the WISE telescope completed its primary mission by mapping the entire sky in four mid-IR bands -- W1, W2, W3, and W4 -- centered at 3.4, 4.6, 12, and 22~$\mu$m, respectively. After the depletion of the solid hydrogen cryostat, the spacecraft was placed into hibernation in early 2011 and reactivated in 2013. Since then, WISE continues to survey the sky in the W1 and W2 channels. Every sky location is visited by the WISE telescope every six months and twelve or more independent exposures are made during each visit. The exception to this rule are regions around the ecliptic poles which are observed much more frequently due to the polar trajectory of the WISE satellite. The South Ecliptic Pole is located in the outskirts of the Large Magellanic Cloud (LMC) -- the area also covered by the OGLE survey, so a number of our LSP variables are located in this region.

We crossmatched our sample of LSP variables to the Near-Earth Object WISE Reactivation Mission \citep[NEOWISE-R;][]{mainzer2014} archive using a matching radius of 6~arcsec and we found counterparts for 89\% objects. The NEOWISE-R light curves cover the time span of six years (2013--2019) and were obtained in two bands: W1 and W2. We decided not to include the photometry from the WISE primary mission (2010--2011), because we noticed that both datasets are not always compatible for individual stars, in particular there are magnitude shifts which would be difficult to account for.

The NEOWISE-R data points for a typical star are distributed over about 12 epochs at intervals of six months. Each epoch consists of 10--24 single-exposure observations obtained within one to a few days. We converted these individual measurements into fluxes, averaged (with outlier rejection) within each epoch, and converted back to the magnitude scale. As a result, we obtained W1 and W2 light curves containing typically 12 points (median value), however about 2000 of LSP variables, mostly in the LMC, have much better-sampled light curves, consisting of several dozen or even over 100 points.

\section{Analysis of the Light Curves}

We inspected about 8500 phase-folded mid-IR light curves of the LSP stars with at least 12 epochs in each of the W1 and W2 bands. For a significant fraction of the sample, the phase coverage turned out to be too sparse to draw definite conclusions about the characteristics of the mid-IR time series. In particular, most stars with LSPs around 1, 1.5, or 2~yr had wide gaps in the folded light curves due to the six-month intervals between the WISE observations. More than half of the examined LSP variables had peak-to-peak amplitudes of the mid-IR light curves below 0.1~mag, which in many cases was comparable to the noise caused by the (``short-period'') red giant pulsation and the errors of the WISE magnitudes. We also rejected some stars exhibiting irregular mean-magnitude variations and significant amplitude changes of the LSP modulation. As a result, we left with about 2500 objects for which the LSP variations could be easily recognized in the mid-IR data, however a detailed analysis of the light curve morphology was possible only for several hundred ($\sim$700) of them.

The general characteristics of the WISE light curves turned out to be very similar to the characteristics of the optical time series, i.e. a broad triangular minimum lasting about half an LSP cycle is visible in the mid-IR data for the majority of large-amplitude LSP stars. The amplitudes of the W1 and W2 light curves are usually smaller than the amplitudes in the {\it I}-band, however the $A_\mathrm{W1}/A_I$ and $A_\mathrm{W2}/A_I$ amplitude ratios range from about 0.3 to over 1 in various objects.

\begin{figure*}
\centering
\includegraphics[width=17.2cm]{fig3.eps}
\caption{Example light curves of LSP variables showing secondary minima in the mid-IR wavelengths. In each panel, blue, orange, and red points show OGLE-IV {\it I}-band, WISE W1, and WISE W2 photometry, respectively. Solid lines show the Fourier series fits to the mid-IR data. The LSPs are given in parentheses above the panels.\label{fig3}}
\end{figure*}

The primary minima in the optical and mid-IR wavelengths occur at the same time in the vast majority of LSP variables. We did not find any systematic phase lags between {\it I}-band and mid-IR light curves for the general sample of LSP stars, although positive or negative phase shifts can be detected in some individual cases. This is in contrast to the findings of \citet{takayama2020} who noticed systematic phase lags between the LSP light curves in the optical and NIR bands.

The most striking feature that distinguishes mid-IR from optical LSP light curves is the presence of secondary minima. In Figure~\ref{fig3}, we present {\it I}, W1, and W2-band light curves of nine LSP variables from the LMC -- all these stars exhibit well pronounced secondary minima in the W1 and W2 channels but this feature is not visible in the optical wavelengths. The LMC sources have the best-sampled WISE light curves (especially in the region of the sky around the South Ecliptic Pole), but we found such a behavior also in about 120 LSP variables in other environments: the bulge and disk of the Milky Way and the Small Magellanic Cloud.

As can be seen in Figure~\ref{fig3}, the secondary minima in the mid-IR light curves have different widths and depths in different LSP variables. Some stars show secondary minima as wide as the primary ones (half the LSP cycle), others exhibit secondary minima lasting only 0.2 of the LSP cycle or less. The secondary minima are generally much shallower than the primary ones, with the exception of very few LSP variables with secondary and primary minima of comparable depths. It is interesting that some of these extreme cases seem to exhibit shallow secondary minima also in the optical range. OGLE and WISE light curves of two such LSP variables in the LMC are shown in Figure~\ref{fig4}.

The efficiency of detecting secondary minima in the mid-IR light curves is strongly correlated with the phase coverage by the observations, as well as amplitudes, periods, and stability of the LSP light curves. Taking into account only relatively well-sampled W1 and W2 time series, we found about 350 candidates for LSP variables exhibiting the secondary minima. On the other hand, we also indicated a similar number of LSP variables with no distinct secondary minima in their mid-IR light curves. The remaining objects from our sample had too sparse photometry to draw strict conclusions about their IR light curve morphology.

\section{Discussion}
As mentioned, the secondary minima in the mid-IR light curves are one of the predictions of the binary model of the LSPs in red giant stars. In this scenario, the red giant has a close-orbiting substellar or stellar companion surrounded by a comet-like dusty cloud. Hydrodynamical simulations of AGB binary systems \citep[e.g.][]{saladino2019,chen2020} show such spiral dusty structures formed by the matter lost by the giant due to stellar wind. The cloud absorbs the stellar radiation over a wide range of the electromagnetic spectrum and re-emits this energy at infrared wavelengths. When the red giant star is obscured by the cloud, we can observe the primary minima in all spectral ranges. When the dusty cloud is hidden behind the red giant, we can detect the secondary minima only in the mid-IR range. The shapes of the WISE light curves confirm that binarity is responsible for the LSP modulation observed in at least 30\% of pulsating red giants.

The amplitudes of the radial velocity curves measured by \citet{nicholls2009} for several dozen LSP variables suggest the presence of a brown dwarf or a very low-mass main-sequence star companion. Thus, the high prevalence of LSPs among long-period variables imply a large fraction of red giant stars with brown-dwarf companions. This is at odds with the so-called ``brown dwarf desert'' -- a paucity of brown-dwarf companions relative to planets around main-sequence stars \citep{marcy2000,grether2006}. To explain this inconsistency, it should be assumed that the low-mass companions were initially planets that increased their mass accreting material from the stellar wind or via direct mass transfer. This means that the majority of LSP giants are orbited by objects formed from the planets and the total population of these variable stars can be used to study the populations of exoplanets in different regions of the Milky Way and in other galaxies.

\begin{figure*}
\centering
\includegraphics[width=12.0cm]{fig4.eps}
\caption{Light curves of two LSP variables showing deep secondary minima in the mid-IR wavelengths (middle and lower panels) and shallow secondary minima in the {\it I}-band (upper panels).\label{fig4}}
\end{figure*}

The idea that planets orbiting red giants can be transformed into brown dwarfs or even low-mass stellar companions is not new \citep[e.g.][]{livio1984,retter2005}. It is also known that planets around red giant stars are on average more massive than those around main-sequence and subgiant stars \citep[e.g.][]{jones2014,niedzielski2015} which suggests that planets grow during the latest stages of their stellar hosts' evolution. According to \citet{wittenmyer2020}, about 7\% of main-sequence stars have a Jupiter-mass planet on wide orbits, but \citet{jones2021} estimated that $33^{+9}_{-7}$\% of giant stars host at least one Jovian planet within 5~AU, which is in agreement with the observed occurence rate of LSP stars among long-period variables. Note also that the majority of planetary nebulae do not have spherical symmetry, which is associated with the presence of stellar or substellar companions that survived the AGB phase of their stars \citep{decin2020}. Thus, the binary explanation of the LSP phenomenon is compatible with theoretical and observational arguments.

Some properties of the dusty cloud orbiting the red giant star can be learned from the shapes of the mid-IR light curves. The widths and depths of the secondary eclipses in relation to the primary eclipses depend on several factors: the structure and size of the obscuring cloud, its maximum temperature and temperature distribution inside the cloud, its transparency in various passbands, and the inclination of its orbit. As a result, the secondary eclipses may be too shallow to be detectable in some LSP stars, which we actually observe. Additionally, the pulsations of the host star, its irregular light variations, and sparsely sampled WISE light curves certainly reduce the effectiveness of our secondary eclipses detection. However, the detailed calculations of the secondary-eclipses detection efficiency must be very complex and are beyond the scope of this paper.

Several of the LSP light curves shown in Figure~\ref{fig3} have secondary eclipses much narrower than the primary ones. Such a behavior can be interpreted as a result of specific density and temperature profiles of the orbiting dusty cloud. Simply, the densest and hottest part of the cloud that emits the bulk of mid-IR light may be significantly smaller than the whole cloud with its comet-like tail. The secondary eclipses can be observed only when this compact bright region of the cloud is hidden behind the star which takes much less time than the primary eclipse, when the star is obscured by the entire cloud. In the case of the W2 light curve of OGLE-LMC-LPV-03550 (lower left panel of Figure~\ref{fig3}), the secondary eclipse seems to begin before the primary eclipse is ended. This may happen when the dusty tail is exceptionally long.

The secondary minima are most often visible in both, W1 and W2, light curves, and the W2 eclipses are usually deeper. Relative depths of the secondary eclipses in both mid-IR channels can give us a hint about the temperature of the dusty cloud. For this purpose, we performed a simple exercise. Using the PHOEBE2 package \citep{prsa2005}, we derived a series of synthetic light curves of eclipsing binary systems containing a $T=3500$~K primary (red giant) and a secondary with temperatures ranging from 1000 to 3500~K (equivalent of the dusty cloud). ATLAS9 \citep{castelli2003} stellar atmospheres were used for the primary. Since the effective temperature values of the secondary fall outside the ATLAS9 range, we used the blackbody stellar atmospheres for this component. Obviously, hotter dusty clouds produce deeper secondary eclipses in both WISE filters, but we also found that the temperature of the cloud is strongly correlated with the W2-to-W1 ratios of the eclipse depths. Namely, for the lowest checked temperature (1000~K), the W2 amplitudes are larger by a factor of 1.5 than the W1 amplitudes, while for the temperatures around 3500~K, the depths of the secondary minima in the W1 and W2 bands are nearly equal. The W2-to-W1 ratios of the secondary eclipse depths derived for the LSP light curves shown in Figure~\ref{fig3} revealed that the temperature of the dusty cloud in most of these objects is in the range 1000--1500~K.

In most LSP light curves, the secondary minima are located symmetrically between the primary minima (at the phase around 0.5) which indicates the circular orbits of the low-mass companions. The tidal circularization theory \citep{zahn1977,soker2000} predicts that the majority of close binary systems containing red giant components should have low eccentricities, because the typical circularization timescales are only several thousand years \citep{nicholls2010}. However, we also found examples of the LSP variables exhibiting clearly non-symmetrically positioned secondary minima which suggests eccentric orbits. We estimate that such asymmetric mid-IR light curves appear in about 10--15\% of the LSP giants with detectable secondary minima. It is striking that similar fraction of eccentric systems was found by \citet{soszynski2004} in the sample of ellipsoidal red giants -- close binary system populating sequence~E in the PL diagram.

\section{Conclusions}
\citet{wood2004} stated that LSPs ``represent the only form of unexplained, large-amplitude stellar variability known at this time''. This is no longer true. We demonstrated that a large fraction of LSP variables exhibit secondary minima in their mid-IR light curves which proves that the physical mechanism responsible for the LSP modulation is the obscuration by a dusty cloud orbiting the red giant with a substellar or stellar companion. In this scenario, this object is a former planet that accreted a significant amount of mass from the envelope of its host star and appear as a brown dwarf or a very low-mass main-sequence star. Thus, the phenomenon known for many decades \citep{oconnell1933,payne1954,houk1963} seems to be a tracer of extrasolar planets.

Careful examination of the mid-IR and optical light curves can unveil properties of the LSP variables. We showed that the majority of LSP systems have circular orbits, as is expected from the tidal circularization theory. Note that the observed non-sinusoidal shapes of the radial velocity curves \citep{hinkle2009,nicholls2009} are not in conflict with this finding \citep{soszynski2014}. Based on the secondary eclipse widths and depths in the two mid-IR passbands, we demonstrated that the orbiting dusty cloud has a non-homogeneous structure and its temperature may reach $\sim$1500~K in the densest part. For lower orbital inclinations or cooler dusty clouds the secondary eclipses may not be visible in the WISE photometry which also happens in some LSP stars.

\acknowledgments

We thank the anonymous referee for the constructive comments that helped us to improve the manuscript. We are grateful to Andrzej Niedzielski for helpful suggestions which improved the paper. This work has been supported by the National Science Centre, Poland, grant MAESTRO no. 2016/22/A/ST9/00009. The OGLE project has received funding from the Polish National Science Centre grant MAESTRO no. 2014/14/A/ST9/00121. This publication makes use of data products from the Near-Earth Object Wide-field Infrared Survey Explorer (NEOWISE), which is a joint project of the Jet Propulsion Laboratory/California Institute of Technology and the University of Arizona. NEOWISE is funded by the National Aeronautics and Space Administration.

Solar and Stellar Astrophysics (astro-ph.SR)


A cosmic abundance standard: chemical homogeneity of the solar neighbourhood and the ISM dust-phase composition

Abstract: A representative sample of unevolved early B-type stars in nearby OB associations and the field is analysed to unprecedented precision using NLTE techniques. The resulting chemical composition is found to be more metal-rich and much more homogeneous than indicated by previous work. A rms scatter of ~10% in abundances is found for the six stars (and confirmed by six evolved stars), the same as reported for ISM gas-phase abundances. A cosmic abundance standard for the present-day solar neighbourhood is proposed, implying mass fractions for hydrogen, helium and metals of X=0.715, Y=0.271 and Z=0.014. Good agreement with solar photospheric abundances as reported from recent 3D radiative-hydrodynamical simulations of the solar atmosphere is obtained. As a first application we use the cosmic abundance standard as a proxy for the determination of the local ISM dust-phase composition, putting tight observational constraints on dust models.

\section{Introduction}\label{intro}
The Sun is unique among the stars because independent~in\-dicators allow its
chemical composition to be constrained with a precision unmatched for any
other star. This can
be done by spectroscopic analysis of its photosphere and by
measurement of solar wind and solar energetic particles. Solar nebula
abundances can be determined from CI
chon\-drites, which are unaltered since the formation of the system.
The wealth of information established the Sun~as~the principal
standard for the chemical composition of cosmic matter
(e.g. Grevesse \& Sauval 1998, GS98; Holweger 2001; Asplund et al 2005,
AGS05).
However, is a 4.6\,Gyr old star indeed representative of the cosmic matter in its
neighbourhood\footnote{
We consider the region
at distances shorter than $\sim$1\,kpc (and $\pm$500 pc in
Galactocentric direction)
as solar neighbourhood in order to minimize bias
due to Galactic abundance gradients, see Fig.~\ref{evol} for a schematic
overview.}
at present?

Ideal indicators for pristine abundances
are unevolved early B-stars of spectral types B0--B2.
Slowly rotating stars are preferred as
their photospheres should be essentially unaffected by mixing of
CN-processed material
\citep[][]{maeder00}. The
atmospheres of early B-stars are also unaffected by atomic diffusion
that gives rise to peculiarities of metal abundances in many
later-type stars \citep[e.g.][]{smith96}.
A major practical advantage is also their relatively simple photospheric
physics, which is represented well by
classical model atmospheres,
unaffected by complications such as stellar winds
or convection.

As a consequence, early B-stars in the solar neighbourhood were subject of
several NLTE studies in the past
\citep[e.g.][]{gies92,kilian92,kilian94,cunha94,daflon99,daflon01a,daflon01b,daflon03,lyubimkov04,lyubimkov05}.
Overall, they found a wide range of abundances, by about a factor $\sim$10,
and an average metallicity of only $\sim$2/3 solar (GS98).
Hence, the impression arose that the solar neighbourhood is chemically highly
heterogeneous, and the Sun anomalously metal-rich compared to young stars.

Both findings are problematic in
terms of Galactic chemical evolution. Dispersal of stellar nucleosynthesis
products increases the metallicity
over time \citep[e.g.][]{chiappini03} and hydrodynamic mixing tends to
homogenize the interstellar medium (ISM)
locally \citep{edmunds75}. Characteristic timescales
for homogenization are short, ranging from 10$^6$--10$^8$\,yrs on scales of
100-1000\,pc \citep{roy95}.

In contrast to the young stars the interstellar gas shows a high
degree of chemical homogeneity in the solar neighbourhood \citep{sofia04},
with the
rms scatter of mean abundances being $\sim$10\%.
However, the ISM gas phase is not suitable as a tracer for cosmic
abundances because of selective depletion
of elements onto dust grains.
Here we reinvestigate the conundrum of inhomogeneous stellar vs. homogeneous
ISM gas-phase abundances in the solar neighbourhood, motivated by the
finding of homogeneous B-star abundances for carbon
\citep[see detailed analysis by][NP08]{nieva08}.

\section{Sample Analysis}
Six bright and apparently slow-rotating early B stars in the solar
neighbourhood -- randomly distributed in
OB associations and in the field, and covering a wide range of stellar parameters --
were observed in early 2005 at ESO/La Silla,
using FEROS on the 2.2\,m telescope.
Spectra with broad wavelength coverage and resolving power $\lambda/\Delta
\lambda$\,$\approx$\,48\,000 were
obtained, at very high-S/N (up to $\sim$800 in the $B$-band).

\begin{table*}
\footnotesize
\caption{Stellar parameters \& elemental abundances\label{params}}
\setlength{\tabcolsep}{1.5mm}
\begin{center}
\vspace{-5mm}
\begin{tabular}{lrrrrrr}
\tableline\\[-2.5mm]
\tableline
                            & \object[HR 6165]{HR\,6165} & \object[HR 3055]{HR\,3055} & \object[HR 1861]{HR\,1861} & \object[HR 2928]{HR\,2928} & \object[HR 3468]{HR\,3468} & \object[HR 5285]{HR\,5285}\\
\tableline
Sp. Type                    &              B0.2\,V &              B0\,III &               B1\,IV &               B1\,IV &            B1.5\,III &                B2\,V\\
Association                 &              Sco Cen &                Field &             Ori OB1b &                Field &                Field &              Sco Cen\\[.5mm]
$d$\,(pc)                   &           152$\pm$20 &           438$\pm$57 &           450$\pm$59 &            481$\pm$63 &           319$\pm$41 &            155$\pm$20 \\
$T_{\rm eff}$\,(K)          &        32000$\pm$300 &        31200$\pm$300 &        27000$\pm$300 &        26300$\pm$300 &        22900$\pm$300 &  20800$\pm$300\\
$\log g$\,(cgs)             &        4.30$\pm$0.05 &        3.95$\pm$0.05 &        4.12$\pm$0.05 &        4.15$\pm$0.05 &        3.60$\pm$0.05 &        4.22$\pm$0.05\\
$\xi$\,(km\,s$^{-1}$)       &              5$\pm$1 &              8$\pm$1 &              3$\pm$1 &              3$\pm$1 &              5$\pm$1 &              3$\pm$1\\
$v\sin i$\,(km\,s$^{-1}$)   &              4$\pm$4 &             29$\pm$4 &             12$\pm$1 &             14$\pm$1 &             11$\pm$2 &             18$\pm$1\\
$\zeta$\,(km\,s$^{-1}$)     &              4$\pm$4 &             37$\pm$8 &              \nodata &             20$\pm$2 &             20$\pm$1 &              \nodata\\
\tableline
$\varepsilon$(He)$^{\rm a}$           &  10.99$\pm$0.05~(20) &  10.94$\pm$0.05~(16) &  10.99$\pm$0.05~(14) &  10.99$\pm$0.05~(14) &  10.99$\pm$0.05~(14) &  10.99$\pm$0.05~(13)\\
$\varepsilon$(C\,{\sc ii})$^{\rm b}$  &   8.27$\pm$0.14~(13) &   8.35$\pm$0.08~(10) &   8.32$\pm$0.10~(19) &   8.28$\pm$0.08~(18) &   8.36$\pm$0.10~(17) &   8.32$\pm$0.08~(20)\\
$\varepsilon$(C\,{\sc iii})$^{\rm b}$ &   8.31$\pm$0.11~(17) &  8.30$\pm$0.05~~~(7) &   8.36$\pm$0.03~(11) &   8.27$\pm$0.02~~~(5) &   8.47$\pm$0.04~~~(2) &   8.42$\pm$0.06~~~(2) \\
$\varepsilon$(C\,{\sc iv})$^{\rm b}$  &           8.34~~~(2) &           8.45~~~(2) &   \nodata & \nodata & \nodata & \nodata\\
$\varepsilon$(N\,{\sc ii})$^{\rm c}$  &   8.16$\pm$0.12~(73) &   7.77$\pm$0.08~(23) &   7.75$\pm$0.09~(61) &   8.00$\pm$0.12~(61) &   7.92$\pm$0.10~(56) &   7.76$\pm$0.08~(47) \\
$\varepsilon$(O\,{\sc i})$^{\rm d}$   &              \nodata &              \nodata &  8.82$\pm$0.03~~~(3) & 8.83$\pm$0.05~~~(5) & 8.82$\pm$0.03~~~(7) & 8.79$\pm$0.05~~~(7) \\
$\varepsilon$(O\,{\sc ii})$^{\rm e}$  &   8.77$\pm$0.08~(51) &   8.79$\pm$0.10~(41) &   8.74$\pm$0.11~(52) &   8.74$\pm$0.09~(46) &   8.80$\pm$0.09~(40) &   8.71$\pm$0.05~(45) \\
$\varepsilon$(Ne\,{\sc i})$^{\rm f}$  &  8.12$\pm$0.05~~~(2) &              \nodata &  8.12$\pm$0.08~~~(9) & 8.11$\pm$0.09~~~(9) & 8.05$\pm$0.09~(10) & 8.07$\pm$0.07~(14)\\
$\varepsilon$(Ne\,{\sc ii})$^{\rm f}$ &   8.14$\pm$0.07~(16) &  8.07$\pm$0.07~~~(8) &   8.08$\pm$0.09~(14) & 8.03$\pm$0.12~~~(8) & 8.06$\pm$0.03~~~(2) & \nodata\\
$\varepsilon$(Mg\,{\sc ii})$^{\rm g}$ &  7.62$\pm$0.03~~~(3) & 7.60$\pm$0.01~~~(2)  & 7.58$\pm$0.10~~~(6)  & 7.56$\pm$0.03~~~(3)  & 7.51$\pm$0.10~~~(6)  & 7.50$\pm$0.05~~~(4) \\
$\varepsilon$(Si\,{\sc ii})$^{\rm h}$ &              \nodata &              \nodata &  7.47$\pm$0.17~~~(2) & 7.56$\pm$0.08~~~(2) & 7.51$\pm$0.10~~~(5) & 7.22$\pm$0.13~~~(6) \\
$\varepsilon$(Si\,{\sc iii})$^{\rm h}$&  7.50$\pm$0.08~~~(8) &  7.48$\pm$0.08~~~(6) &   7.46$\pm$0.11~~~(9) & 7.52$\pm$0.11~~~(8) & 7.53$\pm$0.17~~~(7) &7.29$\pm$0.05~~~(9) \\
$\varepsilon$(Si\,{\sc iv})$^{\rm h}$ &   7.50$\pm$0.04~(10) & 7.51$\pm$0.18~~\,(5) & 7.50$\pm$0.08~~\,(3) & 7.48$\pm$0.14~~\,(2) & 7.50$\pm$0.04~~\,(2)  & \nodata \\
$\varepsilon$(Fe\,{\sc ii})$^{\rm i}$ &              \nodata &              \nodata &              \nodata &              \nodata &          7.38~~~(1) &          7.38~~~(1) \\
$\varepsilon$(Fe\,{\sc iii})$^{\rm j}$&   7.38$\pm$0.12~(17) &  7.49$\pm$0.12~~~(5) &   7.44$\pm$0.09~(33) &   7.48$\pm$0.10~(30) &   7.42$\pm$0.12~(36) & 7.40$\pm$0.09~(32) \\
\tableline
\end{tabular}
\end{center}
\vspace{-3mm}
$\varepsilon({\rm El})$\,$=$\,$\log({\rm El}/{\rm H})$\,$+$\,12, with rms uncertainties and number of analysed lines in parentheses.
NLTE model atoms: H: \citet{przybilla04};
$^{\rm a}$\,\citet{przybilla05};
$^{\rm b}$\,\citet{nieva06,nieva08}; $^{\rm c}$\,\citet{przybilla01}; $^{\rm
d}$\,\citet{przybillaetal00}; $^{\rm e}$\,\citet{becker88},updated; $^{\rm
f}$\,\citet{morel08}, $gf$-values of \citet{fft04} for \ion{Ne}{1}; $^{\rm g}$\,\citet{przybillaetal01}; $^{\rm
h}$\,\citet{becker90}, extended \& updated; $^{\rm i}$\,\citet{becker98};
$^{\rm j}$\,\citet{morel07}\\
\end{table*}}

The quantitative analysis of the sample stars was carried
out following the hybrid NLTE approach discussed by
\citet[][NP07]{nieva07} and NP08.
In brief, line-blanketed LTE model atmospheres were computed with ATLAS9
\citep{kurucz93} and NLTE line-formation calculations were performed using updated
versions of our codes {\sc Detail} and {\sc Surface}.
State-of-the-art model atoms were adopted (see Table~\ref{params}), which
allow atmospheric parameters and
elemental abundances to be obtained with high accuracy.

\begin{figure}[b]
\plotone{f1.eps}
\caption{Comparison with stellar evolution tracks \citep{meynet03}. Observed
and predicted N/C ratios (by mass) are indicated. The inset shows the
location of the sample stars as projected on the Galactic plane, with the
solar neighbourhood (as considered here) schematically outlined.
Crosses denote positions of BA-SGs from the control sample (see text).
\label{evol}}
\end{figure}

Multiple independent spectroscopic indicators were considered simultaneously
for the determination of the atmospheric parameters, effective temperature
$T_{\rm eff}$ and
surface gravity $\log g$: all Stark-broadened Balmer lines
and 4--6 ionization equilibria, of \ion{He}{1/II},
\ion{C}{2/III/IV}, \ion{O}{1/II}, \ion{Ne}{1/II}, \ion{Si}{2/III/IV}
and \ion{Fe}{2/III}. Also, the observed spectral energy distributions
were reproduced \citep{nieva06}. The
redundancy
helps to avoid systematic errors. The microturbulent velocity $\xi$ was determined
by demanding that abundances be independent of line
equivalent widths. Elemental abundances $\varepsilon$(El),
rotational velocity $v \sin i$ and macroturbulence $\zeta$ were determined from
fits to individual line profiles. The results are summarized in
Table~\ref{params}. Stellar parameters and He and C abundances
are identical with those derived by NP07/NP08,
except for HR\,5285, where consideration of additional ionization equilibria
indicated small revisions, though agreement is obtained
within the mutual uncertainties. Spectral types
and spectroscopic distances are also given in Table~\ref{params},
agreeing well with HIPPARCOS parallaxes (HR\,6165, HR\,5285)
or with the association distance (HR\,1861).
The positions of the stars in the $T_{\rm
eff}$--$\log g$-plane are indicated in Fig.~\ref{evol}, where a comparison with
evolution tracks
is made. An overview of the location of the stars
in the solar vicinity is also given there.

The uncertainties in the atmospheric parameters were determined from the
quality of the simultaneous fits to all diagnostic indicators. Statistical
uncertainties for abundances were obtained from the
individual line data (rms values). Systematic errors in the abundances due
to uncertainties in atmospheric parameters, atomic data and the quality of the
spectrum are $\sim$0.1\,dex \citep[NP08,][]{przybillaetal06}, i.e. about as large
as the statistical errors.

\begin{figure*}
\centering
\includegraphics[width=.97\linewidth]{f2.eps}
\caption{Comparison of chemical abundance studies (NLTE) of B-stars in the
solar neighbourhood.
Red bars: present work; full line and thick full line (for the same stars
as in our work): \citet{kilian92,kilian94}; dotted: \citet{cunha94},
\citet{cunha06} for Ne;
short-dashed: \citet{gies92}, excluding bright giants and SGs;
long-dashed: \citet{daflon99,daflon01a,daflon01b,daflon03}; dot-dashed:
\citet{morel08};
triple-dot-dashed: \citet{lyubimkov04,lyubimkov05}. Bin width is $\sigma$/2
of the individual studies. Iron abundances in LTE in all
previous work. HR\,5285 is Si-peculiar and therefore excluded from the
silicon mean abundance. Solar abundances ($\odot$): GS98 \& AGS05 (lower values).
The panel for C is reproduced from NP08 for completeness.
See the text for details.\\ \label{histograms}}
\end{figure*}

Our analysis of each {\em individual} star differs from standard studies
in two main respects: {\sc i}) practically all (un\-blen\-ded) lines of
the ion spectra are analyzed instead of a few selected `good' lines
thus avoiding selection effects, and
{\sc ii}) all parameter indicators (in particular the
ionization equilibria) are closely matched simultaneously,
which has never been achieved before. As a result, practically the {\em entire} observed
stellar spectrum is reproduced closely by the spectral synthesis,
see NP07, NP08, \citet{przybillaetal08} for examples.
This is facilitated by the use of critically evaluated data in the model
atom construction and a (time-consuming) iterative approach for a precise determination of
the stellar parameters
(NP08).
Less accurate photometric $T_{\rm eff}$-estimates as adopted in most previous work
are avoided, as
$T_{\rm eff}$-uncertainties are often the most important sources of systematic
error in the abundance derivation, next to ill-chosen atomic data (NP08) and
$\log g$-uncertainties (NP07).

\section{Chemical homogeneity of the solar vicinity}
The status of previous NLTE abundance studies of early B-stars in the solar
neighbourhood (covering clusters, associations and field stars)
is illustrated in Fig.~\ref{histograms}. A wide range of abundance values
is found for most elements, typically spanning $\sim$1\,dex
\citep[for comparison, such a range is bridged by
the cumulative effect of $\sim$13\,Gyrs of Galactochemical
evolution, see e.g. Fig.~2 of][]{chiappini03}. Moreover,
the abundance distributions peak in most cases at sub-solar values, in particular
when referring to
the solar composition of GS98.
Exceptions are He, where most previous studies find values on average larger than
solar, and Ne (about solar, GS98, from two very recent studies; see also
Lanz et al. 2008 for the case of Ar).
Several of these older B-star studies were combined by \citet{snow96} and
Sofia \& Meyer (2001, SM01, see Table~\ref{abus}) to derive a reference composition,
inevitably resulting in sub-solar average values and a large rms scatter.
The former
discrepancy has since been largely removed from a re-evaluation of solar abundances
(AGS05). However, the status quo in
terms of Galactochemical evolution can only be understood by invoking and fine-tuning extra
processes such as infall/outflow of material and local retention of supernova
products by large~amounts.

On the other hand, our sample of early B-stars implies a high degree of homogeneity
for elemental abundances in the solar neighbourhood, with a scatter of $\sim$10\%,
and absolute values of about solar (GS98 and/or AGS05, see Fig.~\ref{histograms} and
Table~\ref{abus}). The only
exception is N, which is most sensitive to mixing of the atmospheric
layers with CN-processed material
\citep[e.g.][]{maeder00}. In this case the
pristine N abundance may be indicated by the 3 objects with the lowest
value, implying a pristine N/C ratio of 0.31$\pm$0.05 (by mass; error bar
accounts for additional uncertainties).

Although our sample is small, we regard it as representative for the early
B-star population in the solar neighbourhood.
The stars sample the relevant portion of the H-burning phase of the objects in
the HRD in terms of $T_{\rm eff}$ and $\log g$ (see Fig.~\ref{evol}). They also sample
one hemisphere of the solar neighbourhood (inset of Fig.~\ref{evol}), half of them
located in OB associations and the other half in the field.
All six stars were analyzed by \citet{kilian92,kilian94} before,
typically spanning the entire abundance range in her sample of 21 stars (see
Fig.~\ref{histograms}). We therefore also find a chance
selection of stars with similar chemical composition for our sample unlikely.
The wide abundance ranges found in previous work reflect the lower
accuracy of the analyses, while shifts of the abundance distributions
relative to each other reflect systematics, with different temperature
scales being the most important among these.
Our results are supported further by a control sample of six BA-type
supergiants (BA-SGs, Fig.~\ref{evol}), for which mean
values of $\varepsilon$(O)\,=\,8.80$\pm$0.02 and
$\varepsilon$(Mg)\,=\,7.55$\pm$0.07
were derived using the same analysis methodology as applied
here \citep{przybillaetal06,firnstein06}.

The finding of chemical homogeneity for our sample is in
excellent accordance with results from the analysis of the ISM gas-phase in
the solar neighbourhood \citep[][and references therein]{sofia04} and with theory regarding the
efficiency of hydrodynamic mixing in the ISM \citep[][]{edmunds75,roy95}.
Excellent agreement is also obtained with elemental abundances in the Orion
nebula \citep[][E04, see Table~\ref{abus}]{esteban04}, with the exception of
C, which may be a consequence of the atomic data used in the Orion analysis
(see NP08 for the stellar case) plus overestimated dust corrections.

In the following we briefly investigate the impact of this {\em cosmic abundance standard}
on important topics of contemporary astrophysics.

\begin{table*}
\footnotesize
\caption{Chemical composition of different object classes in the solar neighbourhood and of the Sun\label{abus}}
\begin{center}
\vspace{-5mm}
\begin{tabular}{lr@{/}r@{~~~~~~~}rrrrrr@{/}r}
\tableline\\[-2.5mm]
\tableline
        & \multicolumn{2}{c}{cosmic standard}     & Orion              &                   & Young                &  ISM                     & ISM\\
Elem.   & \multicolumn{2}{c}{B stars -- this work}& gas+dust$^{\rm b}$  & B stars$^{\rm c}$ & F\&G stars$^{\rm c}$     &  gas                     & dust$^{\rm d}$ & \multicolumn{2}{c}{Sun$^{\rm e/f}$}\\
\tableline
He      & 10.98$\pm$0.02 & \nodata$^{\rm a}$       &    10.988$\pm$0.003 &           \nodata &              \nodata &  \nodata                 & \nodata        & \multicolumn{2}{c}{10.99$\pm$0.02}\\
C       &  8.32$\pm$0.03 & 209$\pm$15             &      8.52$\pm$0.02  &     8.28$\pm$0.17 &        8.55$\pm$0.10 &  8.15$\pm$0.06$^{\rm g}$ & 68$\pm$26      & 8.52$\pm$0.06 & 8.39$\pm$0.05\\
N       &  7.76$\pm$0.05 &  58$\pm$\phantom{4}7   &      7.73$\pm$0.09  &     7.81$\pm$0.21 &              \nodata &  7.79$\pm$0.03$^{\rm h}$ & \nodata        & 7.92$\pm$0.06 & 7.78$\pm$0.06\\
O       &  8.76$\pm$0.03 & 575$\pm$41             &      8.73$\pm$0.03  &     8.54$\pm$0.16 &        8.65$\pm$0.15 &  8.59$\pm$0.01$^{\rm i}$ & 186$\pm$42     & 8.83$\pm$0.06 & 8.66$\pm$0.05\\
Ne      &  8.08$\pm$0.03 & 120$\pm$\phantom{4}9   &      8.05$\pm$0.07  &           \nodata &              \nodata &  \nodata                 & \nodata        & 8.08$\pm$0.06 & 7.84$\pm$0.06\\
Mg      &  7.56$\pm$0.05 &  36$\pm$\phantom{4}4   &            \nodata  &     7.36$\pm$0.13 &        7.63$\pm$0.17 &  6.17$\pm$0.02$^{\rm j}$ & 34.8$\pm$4.4   & 7.58$\pm$0.05 & 7.53$\pm$0.09\\
Si      &  7.50$\pm$0.02 &  32$\pm$\phantom{4}1   &            \nodata  &     7.27$\pm$0.20 &        7.60$\pm$0.14 &  6.35$\pm$0.05$^{\rm j}$ & 29.6$\pm$2.2   & 7.55$\pm$0.05 & 7.51$\pm$0.04\\
Fe      &  7.44$\pm$0.04 &  28$\pm$\phantom{4}3   &            \nodata  &     7.45$\pm$0.26 &        7.45$\pm$0.12 &  5.41$\pm$0.04$^{\rm j}$ & 27.3$\pm$2.7   & 7.50$\pm$0.05 & 7.45$\pm$0.05\\
\tableline
\end{tabular}
\end{center}
\vspace{-3mm}
$^{\rm a}$\,in units of $\log ({\rm El}/{\rm H})+12$\,/\,atoms per 10$^6$ H
nuclei -- computed from average star abundances (mean values over all individual
lines {\em per element}, equal weight per line);
$^{\rm b}$\,E04;
$^{\rm c}$\,SM01;
$^{\rm d}$\,difference between the cosmic standard
and ISM gas-phase abundances, in units of atoms per 10$^6$ H nuclei;
$^{\rm e/f}$\,GS98/AGS05, photospheric values;
$^{\rm g}$\,\citet{sofia04};
$^{\rm h}$\,\citet{meyer97}, corrected accordingly to \citet{jensen07};
$^{\rm i}$\,\citet{cartledge04};
$^{\rm j}$\,\citet{cartledge06}\\
\end{table*}}

\section{The cosmic abundance standard: first applications}
In general, excellent agreement of our B-star abundances with solar values from recent
3D radiative-hydrodynamical simulations of the solar atmosphere
(AGS05) is obtained. The oxygen value falls between GS98 and AGS05 values
and neon is compatible with GS98.
This opens up new perspectives in the ongoing discussion on helioseismic constraints,
chemical abundances and the solar interior model as reviewed by \citet{basu08}.

Our cosmic abundance standard also facilitates a precise determination of
dust depletion in the local ISM for the primary constituents.
The amount of material incorporated into
dust grains is determined by the difference between our B-star abundances and
the ISM gas-phase abundances,
see Table~\ref{abus}. Accordingly, a composition poor in carbon but rich in oxygen
and refractory elements is indicated.

Such studies were undertaken previously,
using e.g.\ abundances of the Sun, of B stars and of young F \& G stars
\citep[e.g.][SM01, see Table~\ref{abus}]{snow96} as proxies for the determination
of the dust-phase composition, however with mixed success.
In particular, B stars were rejected as reliable indicators as the derived
abundances of material in dust at that time were too low to produce the
observed interstellar extinction. Our study revives B stars as proxies of
the ISM dust-phase composition, and even more so because of the extremely
low abundance scatter compared to all other standards considered so far,
except for the Sun.

The present results imply tight observational constraints on dust models in terms of
carbon abundance. The observed properties of dust grains, as
inferred from the interstellar extinction law, have to be produced by a
rather small amount of carbon, posing a challenge to most dust models
\citep[see e.g.][]{snow95}.
We can carry out an important consistency check, following
\citet{cartledge06}:
the O predicted to be incorporated in grains
from the observed Mg, Si and Fe dust abundances and a
rudimentary dust model agrees with the derived O dust abundance within
the mutual (small) uncertainties. For the rudimentary dust model we assume
silicates to be predominantly MgSiO$_3$, with only a small fraction of Fe
bound in silicates and only a small fraction being of olivine-like
composition. The remaining Mg and Fe fraction~is considered to be
in oxide form (MgO, FeO, Fe$_2$O$_3$, Fe$_3$O$_4$), see e.g. \citet{draine03} for a
discussion of observational evidence.

Finally, we combine our B-star abundances with data for S, Cl and Ar
from the analysis of the Orion nebula (E04) and solar meteoritic
values for other abundant refractory elements (with
$\varepsilon$(El)\,$\gtrsim$\,5, AGS05)
to derive mass fractions for hydrogen, helium and the metals.
Values of $X$\,$=$\,0.715, $Y$\,$=$\,0.271, $Z$\,$=$\,0.014 and
$Z/X$\,$=$\,0.020 characterize the present-day cosmic matter in the
solar neighbourhood \citep[to be compared to protosolar values
$X_0$\,=\,0.7133, $Y_0$\,=\,0.2735 and $Z_0$=0.0132,][]{GAS07}.
These combined abundances are our recommended values
for a wide range of applications requiring an accurate knowledge of the
chemical composition at present (e.g.\ for opacity calculations), examples
being models of star/planet formation or stellar evolution (in particular of
short-lived massive stars), or for the empirical calibration of
Galactochemical evolution models.

\acknowledgments
We express our deep gratitude to U. Heber for support and useful
comments on the manuscript, and thank
M. Asplund and A. Serenelli
for stimulating discussion.
M.F.N. acknowledges support by DFG (grant HE\,1356/45-1).