Stellar abundances of beryllium and CUBES

Stellar abundances of beryllium are useful in different areas of astrophysics, including studies of the Galactic chemical evolution, of stellar evolution, and of the formation of globular clusters. Determining Be abundances in stars is, however, a challenging endeavor. The two Be II resonance lines useful for abundance analyses are in the near UV, a region strongly affected by atmospheric extinction. CUBES is a new spectrograph planned for the VLT that will be more sensitive than current instruments in the near UV spectral region. It will allow the observation of fainter stars, expanding the number of targets where Be abundances can be determined. Here, a brief review of stellar abundances of Be is presented together with a discussion of science cases for CUBES. In particular, preliminary simulations of CUBES spectra are presented, highlighting its possible impact in investigations of Be abundances of extremely metal-poor stars and of stars in globular clusters.


Å (wavelengths from
. However, it has been shown that the identification of the Be I lines was erroneous (Grevesse 1968;Hauge and Engvold 1968). Only the Be II 2 S-2 P 0 resonance lines at ∼3130 Å are now used to determine Be abundances in latetype stars. Other Be lines can be found in the UV below 3000 Å, but are only observable from space and seem to be of limited use (Garcia Lopez 1996). Isotopic shifted lines of unstable 7 Be and 10 Be have been searched for, but were never detected (Warner 1968;Boesgaard and Heacox 1973;Chmielewski et al. 1975).
High signal-to-noise (S/N ) spectra in the near UV is hard to obtain because of atmospheric extinction. In addition, the spectral region of the Be lines is crowded with other atomic and molecular lines (see the solar spectrum in Fig. 1). Some of these blending lines are still unidentified, as for example the one contaminating the blue wing of the 3131.067 Å line (see e.g. Ross and Aller 1974;Garcia Lopez et al. 1995;Primas et al. 1997;Smiljanic et al. 2011).
Early analyses of Be abundances were reviewed by Wallerstein and Conti (1969) and Boesgaard (1976). The solar Be abundance seems to have been first determined by Russell (1929). The approximation of local thermodynamic equilibrium (LTE) has been shown to result in correct abundances for the Sun (Chmielewski et al. 1975;Shipman and Auer 1979;Asplund 2005). The solar Be abundance also seems to be largely insensitive to 3D hydrodynamical effects (Asplund 2004;Asplund et al. 2009). The current value of the solar meteoritic abundance of Be is A(Be) 1 = 1.32 (Makishima and Nakamura 2006;Lodders 2010)  while the photospheric abundance is A(Be) = 1.38 (Asplund et al. 2009). Nevertheless, most 1D LTE model atmosphere analyses of Be in the Sun tend to find values between A(Be) = 1.10 and 1.15. The difference is ascribed to near-UV continuum opacity missing in the computations (Balachandran and Bell 1998;Bell et al. 2001).
Beryllium is not produced in significant amounts by the primordial nucleosynthesis, because there are no stable elements with mass number 5 or 8 to act as an intermediate step in synthesising 9 Be (Wagoner et al. 1967;Thomas et al. 1993;Coc et al. 2012). In addition, Be is rapidly destroyed by proton capture reactions when in regions inside a star with a temperature above ∼3.5 × 10 6 K (Greenstein and Tandberg Hanssen 1954;Salpeter 1955). Therefore, stars do not produce Be.
In stars like the Sun, for example, Be is only present in the external regions of lower temperature (see Fig. 1 of Barnes et al. 1999). In evolved stars, where the convective envelope has increased in size and mixed the surface material with the interior, Be, unlike Li, is usually not detected (Boesgaard et al. 1977;García Pérez and Primas 2006;Smiljanic et al. 2010). Beryllium has never been detected in Li-rich giants (de Medeiros et al. 1997;Castilho et al. 1999;Melo et al. 2005).
Long ago, it was understood that Be can only be produced by the spallation of heavier nuclei, mostly from carbon, nitrogen, and oxygen (Fowler et al. 1955;Bernas et al. 1967). The only known way to produce significant amounts of Be is by cosmic-ray induced spallation in the interstellar medium (ISM), as first shown by Reeves (1970) and Meneguzzi et al. (1971).

Fig. 2
The Be abundances as a function of metallicity in the metal-poor stars analyzed by Smiljanic et al. (2009). Also shown are the slopes expected for primary and secondary elements. Clearly, the inverse process dominates the Be production and it behaves as a primary element Two channels of cosmic-ray spallation might work to produce Be. In the so-called direct process, Be is produced by accelerated protons and α-particles that collide with CNO nuclei of the ISM (Meneguzzi and Reeves 1975;Vangioni-Flam et al. 1990;Prantzos et al. 1993). In the inverse process, accelerated CNO nuclei collide with protons and α-particles of the ISM (Ramaty et al. 1997;Lemoine et al. 1998;Vangioni-Flam et al. 1998;Valle et al. 2002;Prantzos 2012).
If the first channel dominates in the early Galaxy, Be should behave as a secondary element, as its production rate would be proportional to the metallicity of the ISM. If the second channel dominates, Be should instead behave as a primary element. The two behaviors can be distinguished by the analysis of Be abundances in metal-poor stars. For primary elements, one should observe a linear correlation between log(Be/H) and the metallicity [Fe/H] 2 with a slope close to one. For secondary elements, the slope should be around two.
Linear relations with the two slopes are compared in Fig. 2. Also shown are the Be abundances of metal-poor stars determined by Smiljanic et al. (2009). It is clear then that Be behaves as a primary element, meaning that the inverse process dominates (as first suggested by Duncan et al. 1992). The determination of abundances of Be in metal-poor stars was first attempted by Molaro and Beckman (1984) and Rebolo et al. (1988). But it was with the works of Ryan et al. (1992) and Gilmore et al. (1992) that the linear correlation with slope of one became well established. Beryllium abundances in metal-poor stars have been further investigated in several works since then (e.g. Boesgaard and King 1993;Molaro et al. 1997;Boesgaard et al. 1999Boesgaard et al. , 2011Smiljanic et al. 2009;Tan et al. 2009;Rich and Boesgaard 2009;Tan and Zhao 2011).
With the inverse process dominating, and considering that cosmic-rays are globally transported across the Galaxy, then the Be production in the early Galaxy should be a widespread process. It follows that, at a given time, the abundance of Be across the Galaxy should have a smaller scatter than the products of stellar nucleosynthesis (such as Fe and O), as suggested by Beers et al. (2000) and Suzuki and Yoshii (2001). In this case, Be abundances could be used as a time scale for the early Galaxy (Pasquini et al. 2005;Smiljanic et al. 2009). This is an interesting application that still needs to be further constrained. Barbuy et al. this volume, for more details). The properties of CUBES are still not finalized, but it will have a resolution of at least R = 20 000 and should provide access to the wavelength range between 3000-4000 Å.

CUBES (Cassegrain U-Band Brazil-ESO Spectrograph) is a new medium-resolution ground based near UV spectrograph for use at the VLT (Very Large Telescope). CUBES is being constructed by Brazilian institutions together with ESO (see
CUBES will be more efficient than UVES (Dekker et al. 2000) in the near UV. UVES (Ultraviolet and Visual Echelle Spectrograph) is the current spectrograph at the VLT that is able to obtain high-resolution spectra in this region. There is an expected gain of about three magnitudes at 3200 Å. This gain in sensitivity will expand the number of targets that can be observed. Among the main science cases for CUBES are the study of abundances of Be, C, N, O, and of heavy neutron-capture elements in metal-poor stars (see also Siqueira-Mello, this volume and Bonifacio, this volume).
Most of the future new instruments and facilities (such as the E-ELT, European Extremely Large Telescope) are now being optimized for the red and near-infrared spectral regions. Thus, CUBES has the potential of being an unique instrument in terms of sensitivity, spectral range, and resolution. With respect to current ESO instruments, CUBES will outperform both UVES and X-Shooter (Vernet et al. 2011). CUBES should also be more efficient than other groundbased near-UV capable spectrographs, such as HIRES (High Resolution Echelle Spectrometer) at the Keck Observatory and HDS (High Dispersion Spectrograph) at the Subaru telescope.
Regarding space-based telescopes, the HST (Hubble Space Telescope) has two UV spectrographs, the Space Telescope Imaging Spectrograph (STIS) and the Cosmic Origins Spectrograph (COS). Nevertheless, the HST will likely not be operational when CUBES comes online (∼2017/2018). The WSO-UV (World Space Observatory-Ultraviolet), a 170 cm space based telescope, that should be launched in 2016 will provide access to the UV region (see Shustov, this volume). The WSO-UV is however optimized for the region between 1150-3200 Å. It is thus not a competitor but complementary to CUBES.

Extremely metal-poor stars
The linear relation with slope close to one between log (Be/H) and [Fe/H] is well established down to [Fe/H] = −2.5/ − 3.0 (but see also Rich and Boesgaard 2009, for a slightly different view). Going to the extremely metalpoor regime ([Fe/H] ≤ −3.00), however, the situation is less clear.
The detection of Be in two stars (LP 815-43 and G64-12) with [Fe/H] ∼ −3.00 by Primas et al. (2000a,b) seems to suggest some deviation from the linear trend, with a possible flattening of the relation between log(Be/H) and [Fe/H]. The Be abundance of a third star (G64-37), analyzed by Boesgaard and Novicki (2006), seems to be consistent with this flattening, although these authors argue that there is only Fig. 4 Comparison of three spectra with R = 20 000 of a metal-poor star with [Fe/H] = −3.00. The black spectrum was calculated without Be, the blue with log(Be/H) = −13.10 (simulating an abundance plateau), and the red one with an abundance decreased by 0.60 dex. The lines at the red spectra can never be detected evidence for a dispersion of the Be abundances. The Be abundances of these three stars, together with the stars from Smiljanic et al. (2009), and of stars BD+03 740 and BD-13 3442 ([Fe/H] = −2.85 and −2.80, respectively) are shown in Fig. 3. The Be abundances and metallicities of these extremely metal-poor stars were redetermined by Smiljanic et al. (2012).
A Be upper limit of log(Be/H) < −13.80 was determined for the carbon-enhanced metal-poor star BD+44°4 93 with [Fe/H] = −3.80 by Ito et al. (2009Ito et al. ( , 2013. This limit is in principle consistent with an extension of the linear trend between Be and metallicity to lower values of Fe (see Fig. 3). Care in the interpretation is needed because the origin of the carbon enhancement in this type of stars is also not well established. If any kind of transfer of material from evolved stars (which are depleted in Be) is involved, the surface abundance of Be would be diluted.
Whether the flattening exists or not is thus still not clear, because Be abundances have been determined in only a few stars with metallicity around or below [Fe/H] = −3.00 (see Boesgaard et al. 2011;Smiljanic et al. 2012). A number of different scenarios could cause such flattening, each one with its own astrophysical implication. Orito et al. (1997), for example, suggest that a Be plateau, although an order of magnitude below of the current detections, could be the result of an inhomogeneous primordial nucleosynthesis. The identification of such a primordial plateau of Be, similar to the one observed for Li (Spite and Spite 1982), might be an important test of the standard big bang model. Other possibilities include: (1) the accretion of metal-enriched interstellar gas onto metal-poor halo stars, while crossing the Galactic plane (see Yoshii et al. 1995 for the case of beryllium, but see also Hattori et al. 2014 for a more recent assessment of this effect on stellar metallicities); (2) pre-Galactic production by cosmic-rays in the intergalactic medium (Rollinde et al. 2008;Kusakabe 2008); and (3) although a plateau is not predicted, Be enriched stars with [Fe/H] ≥ −3.2 might be associated to population III quark-novae (Ouyed 2013).
To address the existence of the flattening, it is important to expand the number of stars with [Fe/H] < −3.00 where Be abundances have been determined. All additional known stars with this or lower metallicity are too faint to be observed with current instrumentation. CUBES will then offer an unique opportunity to obtain the spectra needed to test this possibility. Comparison of three spectra with R = 20 000 of a metal-poor star with [Fe/H] = −3.50. The black spectrum was calculated without Be, the blue with log(Be/H) = −13.10 (simulating an abundance plateau), and the red one with an abundance decreased by 1.00 dex. The lines at the red spectra can never be detected For this science case, it is important to establish whether a sample of stars has Be abundances at the level of a given plateau or whether the Be abundances follow an extension of the linear relation between log(Be/H) and [Fe/H]. Therefore, the key requirement is to detect the Be lines at lowmetallicities and derive proper abundances, not only upperlimits.
To test if this will be possible with CUBES, spectra of two different stars were simulated, one with [Fe/H] = −3.00 and the other with [Fe/H] = −3.50. Synthetic spectra were computed with the same codes and line lists used in Smiljanic et al. (2009). The spectra of each star was computed with two different Be abundances. The higher abundance for the two stars is the same, chosen to simulate the presence of a plateau. This abundance is log(Be/H) = −13.10, i.e. the Be abundance found for star G64-12 by Primas et al. (2000b). The smaller Be abundance, was decreased by 0.60 dex for the star with The Be abundances reflect a scale where A(Be) = 1.10, as derived by Smiljanic et al. (2009). In addition, all the spectra were computed with the abundance of the alphaelements increased, [α/Fe] = +0.40.
The spectra were computed with three different spectral resolutions, R = 15 000, 20 000, and 25 000, although only the case of R = 20 000 is shown here, as this currently seems to be the favored CUBES value. To simulate the observations, noise was added to the spectra with four different levels in 10 different realizations (kindly computed by H. Kuntschner). The signal-to-noise levels were S/N = 20, 50, 100, and 150 for the case with R = 25 000. For the other cases the S/N was scaled to try to simulate the gain in S/N given by the decrease in resolution. Thus, the values used were S/N = 22, 56, 112, 168 for the R = 20 000 case and S/N = 26, 65, 129, 194 for the R = 15 000 case.
Example are shown in Figs. 4 and 5. It is possible to see that detecting low Be abundances (log(Be/H) ≤ −13.70) is not possible, in the simulated spectra. This is also valid for the cases with different resolutions not shown here. Nev- Fig. 6 Comparison between turn-off stars of NGC6752 with R = 20 000. A difference in abundances of Be = 0.60 dex can be detected in all cases except the lowest S/N case, where a detection is marginal at best ertheless, it is also shown that if there is a plateau around log(Be/H) = −13.10, the Be lines can be detected in both metal-poor stars, if S/N 100. Therefore, with CUBES data it would be possible to detect the flattening, if it has a level of Be abundances similar to the one currently found for stars with [Fe/H] ∼ −3.00. The exact limit of Be detection for CUBES remains to be determined. To calculate that, more information is needed also to understand which level of S/N can be reached.

Stars in globular clusters
Globular clusters host multiple generations of stars (see e.g. Gratton et al. 2012, and references therein). Second generation stars are formed out of material contaminated by proton-capture nucleosynthetic products. The polluters are usually thought to be either first-generation massive asymptotic-giant-branch stars (AGBs) or first-generation fast-rotating massive stars (Ventura and D'Antona 2009;Decressin et al. 2007).
That these peculiarities have a pristine origin, and are not the result of deep mixing inside the stars themselves, became clear when the same chemical pattern was found in turn-off stars (Gratton et al. 2001). The most obvious observational result of mixing pristine and processed material is the appearance of the correlations and anti-correlations between light-element abundances. Proton-capture reactions can decrease the abundances of Li, C, O, and Mg and enhance the abundances of N, Na, and Al. All the globular clusters analyzed so far with sufficient depth show signs of these chemical inhomogeneities (Carretta et al. 2009).
Beryllium, another light element, should also be affected by the polluting material. In particular, and as discussed before, Be is not produced inside stars but only destroyed. Thus the polluting material from the first generation of evolved stars, regardless of the exact nature of the polluting star, is completely devoid of Be. A star formed with only pristine material should have the original Be abundance of that material. Stars with different amounts of polluted material should have diluted the surface Be abundance to different levels. Stars with bigger fraction of polluted material should have stronger correlations and anti-correlations between the light elements and should be strongly depleted in Be.
Lithium abundances should in principle behave in a similar way. However, AGB stars can produce Li via the Fig. 7 Comparison between stars at the metal-poor turn-off of Cen with R = 20 000. A difference in abundances of Be = 0.60 dex can not be easily detected in the lowest S/N case. It seems possible in the three other cases Cameron-Fowler mechanism (Cameron and Fowler 1971). Even though some extreme fine tuning might be needed to explain the full range of oxygen abundances in stars with seemly normal Li, this option can not be completely excluded as of yet (see e.g. Shen et al. 2010). In this situation, Be offers an important complementary test. As Be is never produced in stars, it is an unique tracer of the amount of pollution suffered by the globular cluster stars.
In this science case for CUBES, the objective is to detect differences in the Be abundance among different turn-off stars of the same cluster. These different stars should have been polluted to different levels. If different Be abundances are detected, then one can use them to quantify the fraction of polluting material. On the other hand, if stars with different levels of pollution are found to have the same Be abundance, this would definitely show that all the simple dilution-pollution scenarios currently proposed to explain the chemical properties of globular clusters are too simplistic and not valid. It would show that the true formation scenario was more complex than that.
Measurements of Be abundances have been attempted in two globular clusters so far: (i) In NGC6397, Be was detected in two turn-off stars by Pasquini et al. (2004), and an upper limit was derived for the extremely Li-rich main-sequence star NGC6397 1657 (Pasquini et al. 2014) and (ii) two turn-off stars in NGC6752 were analyzed by Pasquini et al. (2007), resulting in one detection and one upper limit. The results seem to indicate that the stars have the same Be abundance, even though they have different abundances of oxygen (and thus were polluted to different levels). Nevertheless, the results are uncertain due to the low S/N of the spectra (S/N ∼ 10-15), even though they were obtained with long total exposure time (∼900 min). CUBES will allow better spectra to be acquired with shorter exposure times, increasing also the number of globular clusters were Be abundances can be investigated.
To test this science case, synthetic spectra of three different stars were computed. One star represents turn-off stars in the cluster NGC6752, the other stars represent the metalpoor and metal-rich turn-off members of Cen. For each star, two different Be abundances were adopted. The higher abundance in NGC6752 is the one derived by Pasquini et al. (2007), log(Be/H) = −12.05. The smaller abundance was decreased by 0.60 dex. For the stars in Cen, the Be abundances were scaled according to the metallicity Fig. 8 Comparison between stars at the metal-rich turn-off of Cen with R = 20 000. A difference in abundances of Be = 0.60 dex can be detected in all cases -log(Be/H) = −12.6 and −13.2 for turn-off star 1 and log(Be/H) = −11.7 and −12.3 for turn-off star 2: 1. Turn-off star in NGC 6752: T eff = 6400 K; log g = 4.30; [Fe/H] = −1.50; ξ = 1.20 km s −1 2. Turn-off star 1 in Cen: T eff = 6700 K; log g = 4.20; [Fe/H] = −2.00; ξ = 1.20 km s −1 3. Turn-off star 2 in Cen: T eff = 6550 K; log g = 4.20; [Fe/H] = −1.10; ξ = 1.20 km s −1 The same combinations of resolving power and S/N used before were adopted. The results are shown in Figs. 6, 7, and 8. It can be seen that basically in all cases a difference in Be abundances can be detected, the only exception being the metal-poor turn-off of Cen with low-S/N . Therefore, with CUBES it will be possible to use Be as a tracer of the dilution-pollution processes in globular clusters. More sound quantitative predictions have to wait until the CUBES properties are better defined. Then, it will be possible to define the list of clusters where the measurements will be possible and better determine the limit of detection of the Be abundances in these stars.

Summary
Abundances of the fragile element Be can be used in different areas of astrophysics, including studies of the Galactic chemical evolution, of stellar evolution, and of the formation of globular clusters. Only the Be II resonance lines at λ 3130.422 and 3131.067 Å are used for abundance studies. The lines are in the near UV spectral region, a region strongly affected by atmospheric absorption. It is very time demanding to obtain the high-resolution, high-S/N spectra needed to study Be.
Because of that, Be abundances have not been used to its full potential yet. To increase the number of stars where Be abundances can be determined, and to decrease the time cost of obtaining high-quality spectra, new UV sensitive instruments are needed. CUBES is one such instrument. It has an expected sensitivity gain of about three magnitudes at 3200 Å over UVES.
In the near future, CUBES will likely be the only instrument offering the opportunity to measure Be abundances in samples of extremely metal-poor stars and in turn-off stars of globular clusters. Here, preliminary simulations of CUBES-like spectra for these types of stars were presented.
For the first case, the simulations indicate that with CUBES spectra it will be possible to investigate the suggested flattening of the relation between Be and Fe at low metallicities ([Fe/H] < −3.00). This flattening might have implications, for example, to primordial nucleosynthesis models.
With CUBES, it will also be possible to investigate Be abundances in turn-off stars of globular clusters. Beryllium can be used as a tracer of the fraction of polluting material accreted by the second generation stars. If stars with different levels of pollution are found to have the same Be abundances, this would argue against the currently favored dilution-pollution scenarios invoked to explain the chemical properties of globular clusters.