© ESO, 2014

1. Introduction

Our empirical knowledge of the chemical history of the galaxy is based on the assumption that the atmospheres of the main sequence stars that trace its evolution have chemical compositions that reflect that of the interstellar gas at the time when they formed. This is not always the case, and it is essential to understand situations in which stellar surface chemistry does not reflect the bulk or initial chemical composition of the star.

There are physical processes, such as diffusional gravitational settling of trace elements, and levitation of specific ions by radiative acceleration, that can lead to substantial evolution of atmospheric chemistry during the main sequence lifetime of a star. In cool stars (with atmospheric effective temperatures below about T_eff ≤ 7000 K), these separation processes are overwhelmed by deep convective mixing of the outer envelope of the star, which maintains chemical homogeneity, and nearly the initial chemistry. In hot stars (with T_eff ≥ 20 000 K), the intense outwards radiation flux drives a strong stellar wind, stripping material off the stellar surface so rapidly that separation processes act too slowly to be able to compete.

Between these two effective temperature ranges lie the A- and late B-type main sequence stars. Such stars do not have any single powerful process acting to force the surface layers to retain essentially their initial chemical composition. In such stars, the atmospheric chemistry may evolve under the influence of loss of heavy trace ions downwards by diffusive gravitational settling. Chemical evolution may also occur if the atmosphere acquires ions that are driven upwards from the invisible subsurface layers by radiative acceleration due to the outflowing radiative energy flux. Atoms that are driven into the atmosphere by this effect may also be driven out into space and thus lost (Michaud et al. 1976).

Other processes interact with or modify the effects of diffusion in the A and B stars. For example, a magnetic field can greatly impede the flow of material through the upper atmosphere of the star into interstellar space. Large-scale circulation currents in the stellar interior can lead to shear-induced turbulence that mixes subsurface layers. This mixing is strongly dependent on the stellar rotation rate, and is probably ineffective in stars with rotation periods of more than a few days. Accretion from dense interstellar clouds, through which the star passes from time to time, may also alter the surface chemistry (e.g. Havnes & Conti 1971; Havnes 1974; Kamp & Paunzen 2002). Finally, the star may be a member of a double system, and may acquire fresh surface material that is ejected by the evolution of its companion.

As a result, many A and B stars show chemical abundance ratios log N_X/N_H that are remarkably different from those measured in the Sun. Among A and B stars there are several families of distinctive compositional patterns that reflect different sets of physical conditions. The compositional peculiarities vary rather strongly with effective temperature (and thus with the momentum carried by the outflowing radiation that levitates some trace elements). The peculiarities are very different, depending on whether the star has a strong magnetic field or not. There is also a family of peculiarities due to recent accretion from an interstellar cloud. These chemical peculiarities provide powerful probes of invisible processes occurring beneath the visible layers of all kinds of star, particularly of upwards and downwards diffusion (Landstreet 2004). It is therefore of obvious importance to observe and characterise the nature and time variations of these phenomena.

Until now, the study of atmospheric abundances of magnetic A and B (Ap and Bp) stars has been restricted mainly to individual studies of specific stars. The detail of these studies vary from coarse models that roughly describe the magnetic field geometry and abundance variations over the stellar surface (e.g. Landstreet 1988; Bailey et al. 2011, 2012) to more detailed maps of both the magnetic field structure and abundance distributions using spectroscopic observations in all four Stokes parameters (e.g. Kochukhov et al. 2004; Kochukhov & Wade 2010). In this paper, we present the results of a large project to study the time evolution of average surface characteristics in the Ap and Bp stars during the 10⁸ − 10⁹ yrs of their main sequence phase.

A major difficulty of determining how the observed chemical signatures vary with time has been that the ages of isolated (field) Ap and Bp stars can be estimated only very roughly. To solve this problem we have studied surface abundance evolution in a sample of magnetic peculiar A and B stars that are open cluster members. The ages of cluster members can be determined with much better accuracy than those of isolated stars, and the cluster age applies to all its members, which are presumed to have formed essentially contemporaneously (Meynet et al. 1993; Bagnulo et al. 2006; Landstreet et al. 2007).

The following section discusses the observations. Section 3 describes the modelling technique. Section 4 details the observational results. Section 5 explores possible mechanisms to explain the observed trends and Sect. 6 summarises the main conclusions.

2. Observations

Our study used high dispersion spectra of 15 stars, listed in Table 1. The majority of spectra were acquired using the ESPaDOnS spectropolarimeter and the FEROS spectrograph located at the Canada-France-Hawaii-Telescope (CFHT) and the European Southern Observatory’s (ESO) La Silla Observatory, respectively. One spectrum was acquired using the UVES spectrograph at ESO’s Paranal Observatory.

We are studying a sample of stars with a limited range of masses (between 3 and 4 M_⊙) because the physical effects driving atmospheric abundance changes, and the main sequence time-scale over which these effects can act, are expected to depend on mass. Our sample is selected to allow us to reconstruct the atmospheric abundance evolution of a magnetic Bp star of about 3.5 M_⊙.

The stars of our sample are all definite or probable members of open clusters or associations. The age of each cluster was determined by fitting a theoretical isochrone to the cluster members in a Hertzsprung-Russell (HR) diagram (log (L/L_⊙) versus T_eff). For each cluster, the brightest members (those closest to turn-off from the main sequence) are used as the main indicators of cluster age. It has been found that fitting isochrones to the observed stars in the theoretical HR diagram provides a more precise cluster age than doing this fitting in an observational HR diagram, such as V versus B − V. This is because in the theoretical HR diagram the isochrones have a strong hook near the cluster turn-off, while in the observational HR diagram the isochrones tend to a vertical line that does not discriminate clearly between various ages (Landstreet et al. 2007).

For each star, accurate values of effective temperature (T_eff), gravity (log g), evolutionary mass, magnetic field strength and age were adopted from the literature (Landstreet et al. 2007; Bailey & Landstreet 2013b). When values of T_eff and log g were not available, we derived them from available Geneva and uvbyβ photometry. For the Geneva photometry the fortran program described by Kunzli et al. (1997) was used. For the Strömgren uvbyβ photometry we used a version of the fortran program “UVBYBETANEW” (see Napiwotzki et al. 1993) that corrects the T_eff of the magnetic Bp stars to the appropriate temperature scale (see Landstreet et al. 2007). As per the discussions of Landstreet et al. (2007) and Bailey & Landstreet (2013b), the uncertainties in T_eff and log g were taken to be about ±500 K and 0.2 dex, respectively. Table 2 lists the properties of each star, including its designation, associated cluster, age, mass, T_eff, log g, and vsini. The final column shows the root-mean-square magnetic field strength (B_rms) as computed from whatever modern measurements of mean longitudinal magnetic field ⟨B_z⟩ are available (Landstreet et al. 2007, 2008); this quantity is the most useful value we have to characterise the magnetic field of each star. Where applicable, the appropriate reference is given for each parameter. Unreferenced parameters were derived for this study. We note that the photometrically determined values of T_eff and log g vary systematically with increasing age from about 13 500 K and 4.4 to 10 000 K and 3.5, as expected for the evolution of a single star of about 3.5 M_⊙ from ZAMS to TAMS.

3. Modelling technique

To determine the atmospheric abundances of the magnetic Bp stars, the fortran program zeeman was used (Landstreet 1988; Wade et al. 2001a). zeeman is a spectrum synthesis program for stars having magnetic fields. It assumes a magnetic field geometry which is modelled as a simple co-linear multipole expansion, with the strength of the multipole components, the inclination i of the rotation axis, and the obliquity β of the magnetic field axis to the rotation axis specified. zeeman interpolates an appropriate stellar atmospheric structure from a pre-tabulated grid of ATLAS 9 atmospheric models, based on the T_eff and log g assumed. The atomic data for individual spectral lines are taken from the Vienna Atomic Line Database (VALD; Kupka et al. 1999, 2000; Ryabchikova et al. 1997; Piskunov et al. 1995). For all stars, a uniform atmospheric abundance distribution vertically and over the stellar surface was assumed. The microturbulence parameter has been set to zero for all abundance determinations for two reasons. Firstly, it has been found that even non-magnetic stars show no evidence of microturbulence between about 11 000 and 14 000 K (Landstreet et al. 2009); and secondly, it is very probable that the magnetic field is able to suppress convective motions in the atmosphere because of the large energy density in the field.

Fig. 1 Spectrum synthesis of the region 5000−5080 Å for two spectra of HIP 109911. The observed spectra are in black and the model fits are in red.

zeeman searches for an optimal fit between the synthetic and observed spectra by means of a reduced χ² fit of the computed spectrum to the observed one. Multiple spectral windows can be synthesised simultaneously and zeeman automatically provides as output the best values for the radial velocity v_R and vsini. The abundance (log N_X/N_H) of one element at a time is optimised by identifying unblended lines of that element and then fitting, as well as possible, the observed spectral lines of that element. The stars modelled vary in T_eff from about 10 000 to 14 000 K and therefore share many spectral lines in common. For consistency, we endeavoured to deduce elemental abundances from the same sets of lines for all stars. Table 3 lists the lines used for modelling all the stars. As far as possible, lines with a range of strengths were used to deduce the final abundances. For each element, we adopted an uncertainty consistent with the change in abundance from the best-fit model that was necessary to produce an unsatisfactory fit in the spectral window (determined by visual inspection). zeeman cannot model simultaneously spectral lines that are widely separated. In such instances, at least two lines (in different spectral windows) were fit separately and the average abundance deduced from the lines was adopted. The uncertainty in this case was estimated from the observed scatter between the computed values of the different spectral lines. All of our stars have measurements of the line-of-sight magnetic field ⟨B_z⟩, either previously published (Bagnulo et al. 2006; Kochukhov & Bagnulo 2006; Landstreet et al. 2008) or from our own unpublished results. For two of the stars in our sample, magnetic field models are available from detailed studies: HD 133880 (Bailey et al. 2012) and HD 147010 (Bailey & Landstreet 2013a). In these cases, we adopted the published geometries. When no detailed magnetic field model was available, we simply assumed a dipolar magnetic field that was approximately three times B_rms and computed the spectrum with the line of sight parallel to the magnetic axis. For the majority of the stars in our sample we have multiple observations. In the cases where more than two spectra were available, we chose the two that exhibited the greatest differences. In all cases, the derived abundances are the average of results from the two modelled spectra, with the associated uncertainties propagated accordingly.

3.1. Abundance analysis

Atmospheric abundances of He, O, Mg, Si, Ti, Cr, Fe, Pr and Nd were determined for the 15 cluster magnetic Bp stars in our 3.5 M_⊙ sample. In Table 4, we tabulate the mean abundances log (N_X/N_H) for each star. The first three columns recall the star designation, age (log t) and T_eff from Table 2. The subsequent columns list the average abundance values, with their associated uncertainties, for He, O, Mg, Si, Ti, Cr, Fe, Pr and Nd. For reference, the solar abundance ratios are also shown (Asplund et al. 2009). Note that the abundance scale we use can be converted to the common scale, having the logarithmic abundance of H as +12, by adding 12 to all our values. Figure 1 provides an example of the quality of fits we achieved for each star.

We know that serious discrepancies (of the order of 1 dex) are found between abundances derived from lines of Si ii compared to those derived using lines of Si iii, and that this situation exists in most, or all, magnetic stars within the T_eff range of this study (Bailey & Landstreet 2013b). In general, for magnetic Bp stars the abundance of Si deduced from lines of Si iii is significantly larger than the value found using lines of Si ii. The study by Bailey & Landstreet (2013b) suggests that this discrepancy is mostly due to strong vertical stratification of Si, similar to the stratification profiles computed by (Leblanc et al. 2009) for stellar atmospheres with T_eff = 8000 and 12 000 K, and shown in their Fig. 9. We have determined Si abundances using only lines of Si ii.

Stratification has been directly deduced from observed spectra in the vertical distribution of Fe and other elements (e.g., Bagnulo et al. 2001; Wade et al. 2001b; Ryabchikova et al. 2005), mainly in cool magnetic Ap stars. We expect this phenomenon to occur in the T_eff range discussed in this paper as well. Because both observation (Bagnulo et al. 2001; Ryabchikova et al. 2004a) and theory Leblanc et al. (2009) suggest that the vertical abundance variations may be of the order of 1 dex or more, this phenomenon raises quite serious questions about the precise meaning of atmospheric chemical abundances derived assuming that the elements are uniformly distributed in the vertical direction. This is similar to the ambiguous meaning of abundances derived for (horizontally) patchy stars using models that assume uniform horizontal abundance.

Although these complications make the meaning of abundance determinations rather inexact, such simplified abundance measurements nevertheless form a useful “instrumental” system, which reveals differences between different stars satisfactorily. The typical symptom of vertical abundance stratification is that when we find a homogeneous model spectrum that fits lines of an element of intermediate strength, that same model predicts lines that are stronger than the strongest observed lines, and weaker than the weakest observed lines (see Figs. 2 of Bagnulo et al. 2001). Thus, by simply fitting on average the spectral lines of the element under study, we hypothesise that we are deriving a fairly stable mean abundance which has a similar meaning in stars that do not differ greatly in T_eff. However, the ambiguity due to probable vertical stratification and horizontal inhomogeneities should be kept in mind in evaluating our results.

Fig. 2 Atmospheric abundances of He, O, Mg, Si, Ti, Cr, Fe, Pr, and Nd plotted against log t. Each plot contains the solar abundance ratio for that element (solid red line) and the best-fit linear regression (dashed blue line).

4. Results

4.1. Abundance variations with time

Figure 2 shows nine plots, one for each of the nine elements studied in this paper. In every panel we show the mean abundance found for each of the 15 stars of our sample as a function of the age of that star, log t. Uncertainties in age are taken from Table 2, and in abundance from Table 4. Although the error bars are large, in six of the nine panels a trend is clearly present. We can immediately see that on average several of the elements studied decrease or increase in abundance as our average 3.5 M_⊙ star ages. It is essentially because we have good age resolution that we are able to detect these variations.

For each panel we show the best-fit linear regression to log N_X/N_H versus log t (dashed blue line). Table 5 tabulates the slope of the best-fit line with its uncertainty, as well as the significance of the slope, σ_D. This table confirms the visual impression that the slopes are significantly (or nearly significantly) different from zero for He, the Fe-peak elements, and the rare earths, but not for the light elements: these elements evolve substantially in abundance through the main sequence life of magnetic Bp stars of about 3.5 M_⊙.

We also show for comparison the solar abundance of each element studied (solid horizontal red line). It is clear that overabundance or underabundance relative to solar abundance is systematically present for all the elements studied. We discuss individual elements below.

4.1.1. Helium

One of the defining properties of magnetic Bp stars is the fact that they are visibly underabundant in He compared to normal stars of the same T_eff values. This is true of all of the magnetic Bp stars in our sample. Unexpectedly, we observe a clear evolutionary increase in the abundance of He with stellar age. Substantial underabundances (of the order of 2 dex compared to the Sun) are found in stars near the ZAMS and those with ages up to about log t ~ 8. Older stars in our sample (nearer the TAMS) show only modest underabundances of He, about 1−0.5 dex smaller than the solar ratio.

4.1.2. The light elements: oxygen, magnesium and silicon

Oxygen, magnesium and silicon are systematically different from solar abundance ratios by values up to about 1 dex. No apparent trend in abundance with age is seen in any of these three elements. Typically, both oxygen and magnesium are underabundant compared to the solar abundance ratios with a few exceptions. For oxygen, only HD 133880 is overabundant compared to in the Sun; this anomalous behaviour is well-documented by Bailey et al. (2012). HD 133880 and HD 45583 have nearly solar abundances of magnesium, nearly 0.5 dex larger than the other stars of the sample. It is not clear why these two (very similar) Bp stars depart so far from the general behaviour of other stars in our sample. In general, silicon appears to be overabundant compared to the solar abundance ratio throughout the main sequence evolution. Moderate Si overabundance is common for magnetic Bp stars between effective temperatures of about 10 000 and 15 000 K.

Although there are no clear trends in abundance for any of these three elements, all show a substantial star-to-star variations around the mean.

4.1.3. The Fe-peak elements: Ti, Cr and Fe

All three Fe-peak elements studied here, Ti, Cr and Fe, are quite overabundant compared to their solar ratios early in the main sequence phase. However, the level of excess abundance is larger for the lower abundance elements Ti and Cr than for Fe. In fact, Cr reaches almost the same level of total abundance as Fe close to the ZAMS. For these elements, definite decreases in atmospheric abundance are observed with time and by the end of the main sequence all three elements are not much more abundant than the solar values.

4.1.4. The rare earth elements: Pr and Nd

Close to the ZAMS, both these elements appear very overabundant (more than 4 dex) relative to the Sun. There is a clear decrease in the abundances of each element with age, although even at the TAMS both elements are still typically 2–3 dex overabundant. The apparently smaller scatter around the regression lines compared to other elements is produced mainly by the very large range seen in the y-axis of these plots.

The single outlier with low abundance near the TAMS is HD 162576. This star also has the smallest Si abundances.

Fig. 3 Same as Fig. 2 but for Brms.

Fig. 4 Shown are Cr (left) and Fe (right) abundances verses Teff. The black triangles are data taken from Ryabchikova et al. (2004b) and Ryabchikova (2005), the green squares are from Bailey & Landstreet (2013b) and the blue circles are data from this paper.

4.2. Abundance variations with magnetic field strength

Landstreet et al. (2007, 2008) discovered that magnetic field strength decreases during the main sequence lifetime of magnetic Bp stars. Since we see trends in abundance with time for He, Ti, Cr, Fe, Pr and Nd, we expect that there will be correlations between abundances and stellar magnetic field strength. Figure 3 plots abundances versus log B_rms in the same manner as for log t. In all the elements for which a trend was observed versus time (He, Ti, Cr, Fe, Pr and Nd) a correlation is also found for abundance versus log B_rms, but in the opposite sense: helium abundance decreases with increasing magnetic field strength and the Fe-peak and rare-earth elements increase in abundance with increasing magnetic field strength. It is clear that the correlations seen in Fig. 3 are those expected from the known decline of ⟨B_z⟩ with age together with the newly discovered variations of abundance with age. However, the observed correlations may also contain some information about how the time evolution of atmospheric abundance is modified by magnetic fields of various strengths. This may become clear with more detailed modelling.

We note that the results shown in Fig. 3 seem to be the explanation of the correlation between the Geneva photometric measurements of the 5200 Å depression (using the Z index) and field strength, as discussed by Cramer & Maeder (1980). This effect has also been exploited recently by the Special Astrophysical Observatory (SAO) group (e.g. Kudryavtsev et al. 2006, 2011; Kudryavtsev & Romanyuk 2008) to search for particularly large magnetic fields in Ap/Bp stars.

4.3. Cr and Fe abundances versus effective temperature

Our results can be compared to previous studies of magnetic Bp stars. Ryabchikova et al. (2004b) presented a comprehensive comparison of Cr and Fe abundances in roAp (rapidly oscillating Ap) stars versus effective temperature. However, their study was mostly restricted to stars at or below around 10 000 K. Ryabchikova (2005) added stars between about 10 000 and 15 000 K and our current study, as well as a previous study by Bailey & Landstreet (2013a), adds 23 stars in that same temperature range. The results are shown in Fig. 4. Our observations confirm the tendency, already visible in the data of Ryabchikova, for the star with T_eff above 8000 K to exhibit substantially larger abundance dispersion than is shown by cooler stars. Our new results suggest that this extra scatter might be due to mixing stars of quite different ages at a given value of T_eff, although we do not yet have enough abundance measurements of magnetic stars with both known mass and known age to confirm or reject this idea. It appears that for Cr the largest observed abundances peak at around 10 500 K and decline with increasing effective temperature. A similar trend is seen with Fe, but to a lesser extent. There is only a broad flat peak in the abundance of Fe with T_eff. These results support the conclusions of Ryabchikova et al. (2004b) that the maximum abundance for Cr and Fe both approach the same value of between − 3 and − 4 dex.

Fig. 5 Position of the stars of this study in the Hertzsprung-Russell diagram, compared to the main sequence evolution tracks for stars with 3.0, 3.5 and 4.0 M⊙. Also shown are the ZAMS, the TAMS, and the points where half the main sequence lifetime is completed.

5. Discussion

It is well known that during its evolution on the main sequence, a star undergoes an extensive process of climate change. The atmospheric effective temperature decreases by about 30%, the radius increases by a factor of order three, and the gravity decreases by a factor of ten. In effect, by studying the changes in average atmospheric chemical composition of stars in our sample, with its limited mass range, as a function of cluster age, we are observing the time evolution of the atmospheric chemistry of a 3.5 M_⊙ star through its main sequence life. Figure 5 shows an HR diagram for the stars in this study, compared to theoretical stellar evolutionary models (Girardi et al. 2000).

What we have observed, for the first time, is a clear evidence of the secular evolution of the atmospheric chemistry of the magnetic peculiar B-type (Bp) stars during the main sequence phase. That is, we have observed how the consequences of diffusion processes are modified by time and by stellar climate change. Figure 2 shows that the abundances of He, Ti, Cr, Fe, Pr, and Nd change monotonically with time. In contrast, the light metals O, Mg and Si also have peculiar abundances, but show no significant trends with age.

It is generally believed that the peculiar chemical abundances found in the atmospheres of magnetic Ap and Bp stars are the result of microscopic diffusion, in competition with other processes such as turbulent diffusion, convection, meridional circulation, well-mixed mass loss, accretion, and the effects of a magnetic field.

The basic ideas of how diffusion can lead to anomalous atmospheric abundances have been discussed for some years (Michaud et al. 1976; Vauclair 1983). Essentially, in a quiescent stellar plasma, the gravitational field tends to cause trace atoms to diffuse downwards relative to the dominant hydrogen gas, but the outwards flow of radiation from centre to surface is responsible for a force which tends to levitate some trace atoms up through the atmosphere, if they absorb at many wavelengths and are not too abundant. In slowly rotating middle main sequence stars, which have relatively weak mixing near the surface, this phenomenon can lead to anomalous atmospheric chemical abundance.

Detailed and quantitative modelling of diffusion is actually quite complicated. The average upwards acceleration per ion of a particular species depends on the intensity of the radiation (which increases with effective temperature T_eff), on the specific atomic transitions that the ion can undergo (especially those arising from low-lying energy levels), and on the number density of the trace ion. In general, as the fractional number density of the trace ion increases, the acceleration per ion decreases. For many trace ions of low abundance (having, say, number density less than about 10^-6 that of H) in stars of T_eff above about 7000 K, the upwards acceleration g_R due to the radiation field is larger than the local acceleration of gravity g, and the ion diffuses upwards rather than downwards. In slowly rotating middle main sequence stars, which have relatively weak mixing near the surface, this phenomenon can lead to anomalous chemical abundances in the stellar atmospheres.

The actual vertical variation of radiative acceleration for various ions in specific stellar models has been calculated by a number of groups. The calculations fall into two general types. Because the acceleration per ion decreases as the fractional abundance of the ion increases (due to saturation of the spectral lines that the ion absorbs), an ion which is levitated at low abundance can become locally abundant enough that g_R = g, a situation in which the tendency of that trace ion to diffuse upwards or downwards vanishes. It is possible to determine the fractional abundance for a specific atom as a function of position in a region (e.g. with height in the stellar atmosphere) for which this condition is satisfied everywhere in the region, and diffusion of this atom ceases. The run of abundance through the region satisfying this condition is known as the “equilibrium abundance distribution”. Some recent computations of equilibrium abundance distributions in the atmospheres of the late B stars are re reported by LeBlanc (2003), Leblanc et al. (2009) and Stift & Alecian (2012).

Recent work in this field has focussed on re-computing the structure of atmospheres in which equilibrium abundance distribution computations suggest the presence of strong vertical variations in abundance, in order to have the atmospheric structure be consistent with the abundance stratification (LeBlanc 2003; Leblanc & Monin 2004; Shulyak et al. 2009; Leblanc et al. 2009; Stift & Alecian 2012). It is found that when large variations in abundances of abundant elements with altitude are present, the atmospheric structure is substantially perturbed.

However, an equilibrium abundance distribution may not be achievable in a region if not enough atoms of the species is available from below. In this second case, atoms may diffuse slowly up into the region from below and at the same time be removed from the top. This situation could settle into a stationary state in which the flux of atoms through the region is constant with height, and the abundance distribution is unchanging but the actual abundance is smaller than the value leading to equilibrium. In principle, if a large enough volume of the star is considered, the evolution of abundance of an atom with radius, including through the atmosphere, could be computed as a function of time. Because the time-scale for evolution varies strongly with density, this is a stiff problem. Solutions on the scale of a stellar envelope, to try to explain the chemical abundances in the atmospheres of metallic-line (Am) stars, including either a (mixed) stellar wind or deep, but weak, turbulent mixing have been reported by Michaud et al. (2011) and Vick et al. (2010, earlier work on abundance evolution is cited there).

The time variation with height of abundance of a very unabundant element diffusing up through a stellar atmosphere has recently been studied by Alecian et al. (2011), who find that a stationary state of constant particle flux with height is generally achieved after a few hundred years, provided that the atoms continue to diffuse out of the top of the atmosphere. The actual value of the flux, and the run of abundance with height, is set by the abundance of the atom supplied at the bottom of the atmosphere from the reservoir below. We refer to this situation as “stationary flow-through”.

We now apply these general ideas to the elements studied here, in stellar atmospheres of log g ~ 4 and T_eff ~ 10−13 000 K. In the panels of Fig. 2, it appears that we can identify several different cases.

5.1. Abundant light elements: He, O

The general behaviour of these elements was already predicted by Michaud et al. (1976). For the very abundant light elements He and O, the radiative acceleration for T_eff ~ 12 000 K is much too weak to support a solar abundance of these elements (cf. g_rad calculations for T_eff = 12 000 K by Hui-Bon-Hoa et al. 2002). These two elements are expected to diffuse downwards into the stellar envelope below the atmosphere until the relative abundance has fallen enough that radiative acceleration can support them in the atmosphere. Consequently, He and O are expected to have abundances well below the solar values, and this is observed.

The diffusion of He in stellar atmospheres with well-mixed winds (including H) has been studied in more detail by Vauclair et al. (1991), Krticka & Kubát (2004) and Théado et al. (2005); all groups confirm that without a mixed wind from the star, that He should diffuse downwards until quite low relative abundance is reached. However, we know of no published predictions of how low the He abundance should drop before an equilibrium distribution is achieved. It may well be that some degree of turbulent mixing with sub-atmospheric layers is required to keep the abundance of He as large as is observed (Fig. 2).

The increase in He abundance with stellar age that we observe is particularly unexpected. Radiative levitation is so weak (due to the shadowing of resonance lines of He i by the very strong H continuum bound-free absorption below the Lyman limit) that it is not expected to play a significant role in the observed He abundance, and especially not in its increase with decreasing T_eff. The observed increase of He abundance towards the solar mixing ratio suggests either that there is significant and increasing mixing upwards of envelope He by some process unrelated to diffusion, or that there is accretion of (He-rich) interstellar gas.

Landstreet et al. (1998) have studied the diffusion of O in stars in our mass range, using a simple approximation to estimate the radiative acceleration. Extrapolating their results, it appears that radiation may be able to support O at an abundance one to two dex below the solar abundance. This conclusion is supported by the very small g_rad value for O at 12000 K found by Hui-Bon-Hoa et al. (2002). The rather mild underabundance of O that we observe (about 0.5 dex) is probably higher than the level that would be found by an equilibrium abundance calculation. Furthermore, as the stars in our sample age, T_eff and log g both decrease. The change in T_eff typically reduces g_R. However, the decrease in log g means that a smaller value of g_R is required to support ions of O, which in turn means that a larger abundance can be supported. It is not clear which of these two effects dominates, and in fact the observed lack of significant variation in the O abundance with age suggests that if the observed abundance is the result of radiative levitation, the two effects roughly cancel. We do not have any explanation at present for the slight overabundance of O observed for a single star, HD 133880.

It would clearly be of interest to have available some published results of computed equilibrium abundances of He and O in our temperature range, and especially to have such calculations follow the evolution of a star of 3.5M_⊙ from T_eff = 13 000 K and log g = 4.5 to T_eff = 10 000 K and log g = 3.5, in order to determine the importance of radiative levitation of He and O in the observed abundance evolution of the stars of our study.

5.2. Light metals: Mg and Si

The Mg abundance, based generally only on the 4481 Å line of Mg ii, appears to be about 1 dex below the solar abundance, or log (N_Mg/N_H) ≈ −5.2. This may be compared with the predicted Mg abundance profile for a star of T_eff = 12 000 K of Leblanc et al. (2009) and of Alecian & Stift (2010). The abundance predicted by the equilibrium calculation (g_R = g throughout the atmosphere) is mildly non-uniform through the line-forming region, but is of about this same value. Thus it appears that Mg may be an element in which the supply of atoms from below has been large enough to allow the development of an equilibrium stratification.

No strong trend of abundance with age is observed. We have no explanation at present for the two stars (HD 133880 and HD 45583) that deviate strongly from the mean behaviour, with Mg abundances about 1 dex larger than other, younger and older, stars.

Since Mg may well be described in the stars of our sample by the equilibrium abundance distribution, it would be of great interest to have calculations of the equilibrium atmospheric abundance of Mg following the evolution of T_eff and log g for our 3.5 M_⊙ stars.

Silicon equilibrium atmospheric abundance has been studied by a number of authors, including Alecian & Vauclair (1981), Leblanc et al. (2009), and Alecian & Stift (2010). While detailed results differ, the overall conclusion is that Si is expected to be of order 1 dex underabundant in the atmosphere at T_eff ~ 12 000 K. Instead, as has been found in the past, we observe an overabundance of about 1 dex at all ages. This situation has been a long-standing puzzle. Because it is not clear how to obtain an atmospheric abundance that is nearly 2 dex larger that the maximum value that can be supported by radiation pressure, various explanations have been probed, such as support by a horizontal magnetic field, non-LTE effects, etc., but none have been found to offer convincing explanations of the observed overabundance at all ages.

5.3. Iron peak elements: Ti, Cr and Fe

Expected abundances for these elements based on equilibrium have recently been computed for stars in our mass range by Leblanc et al. (2009), Alecian & Stift (2010), and (for Fe) by Stift & Alecian (2012). The results are most extensive for iron. The computed equilibrium abundances are only qualitatively in agreement with one another, and depend on still uncertain physics, particularly on how to treat redistribution of momentum between ionisation stages, but also on the effects of the magnetic field, and the chemical composition assumed for the computation. However, all the computations of equilibrium abundances agree that Fe in the line forming region of stars in our mass range should reach equilibrium at an abundance of around −3 to −3.5 dex relative to H. This is reasonably consistent with the values that we observe (Fig. 2), so this may well be an element for which the equilibrium assumption of zero diffusion is approximately correct. This implies that an adequate supply of Fe ions is available from below the atmosphere to replenish atmospheric atoms lost to space during the initial period of equilibration, and to keep the abundance high enough to satisfy the equilibrium condition.

It is clear from the few values of T_eff for which equilibrium Fe abundances have been computed that the radiative acceleration and the equilibrium abundance decrease with decreasing effective temperature. However, there is no computational information on how the equilibrium abundance varies as log g decreases from 4.5 to 3.5, except for the qualitative result that the effect of the decrease in log g should be to make it possible for a given radiative force to support more atoms. Since the observed evolution of Fe abundance is that the abundance decreases with time, it appears that the decrease in radiative acceleration with decreasing T_eff may dominate. This is certainly a question that could be studied by computing an appropriate series of equilibrium abundance models following the evolution of a 3.5 M_⊙ star.

The available equilibrium calculations for Ti and Cr in the mass range of our observations by Leblanc et al. (2009) have been cut off in the line-forming region because of artificial limits of 1000 times the solar abundances imposed on the calculations. Thus we may suppose that without these limits, the equilibrium calculations would imply equilibrium abundances of Ti and Cr of log (N_el/N_H) similar to those of Fe, in the range of −3 to −4. This is confirmed by the calculations of Alecian & Stift (2010). Although we observe overabundances of these two elements of up to 2 dex when the stars are young, the Cr abundance seems to be generally less than the equilibrium value, and the Ti abundances almost certainly are below equilibrium even at young ages. These may well be elements for which the abundance would be closer to that predicted by the assumption of a stationary flow-through state, with atoms fed into the atmosphere from below and lost from the top (cf. Alecian et al. 2011). In this case the limiting factor determining the atmospheric abundances is the available number density of Cr or Ti brought up from the envelope by diffusion.

As for iron, it is not known whether the probable decrease in radiative acceleration with declining T_eff, or the increased abundance that can be supported as log g decreases, is the dominant effect during main sequence evolution. It appears that the effect of decreasing radiative acceleration, which may well reduce the supply of Cr and Ti at the bottom of the atmosphere and thus may reduce the supply if these elements are in a stationary flow-through state but not in equilibrium, may dominate, as the abundances of both these elements are observed to decrease strongly with increasing stellar age (cf. Alecian & Stift 2010). It would be of great interest to have computations of equilibrium abundances of Ti and Cr following the evolution of a 3.5 M_⊙ star for more direct comparison with our results.

5.4. Trace heavy elements: Nd and Pr

There are no calculations of the equilibrium abundances of heavy elements such as Nd and Pr near the T_eff temperature range of interest to us here, except for the equilibrium stratification calculated for a simplified artificial low-abundance element with T_eff = 12 000 K by Alecian et al. (2011). This artificial element is expected to behave somewhat like Hg, and the equilibrium abundance computed for it is about − 8 dex without a magnetic field, and possibly even smaller high in the atmosphere in the presence of a strong magnetic field. It is not clear to what extent this result is applicable to Nd or Pr. More realistic computations would clearly be of great interest.

In considerably cooler magnetic Ap stars, it is known from modelling of observed spectra that these elements are strongly stratified, with the abundance of Nd as much as 4 dex more abundant high in the atmosphere than near optical depth unity (Mashonkina et al. 2004; Nesvacil et al. 2013). However, we know of no similar stratification models of rare earths in Ap stars in our temperature range. The stars in our sample are hot enough that we have not been able to identify the Nd ii lines that would make stratification analysis possible.

We thus have no real evidence as to whether Nd and Pr at the temperatures of our stellar sample have approximately equilibrium abundance distributions, or are in a stationary flow-through state in which these elements are entering the bottom of the atmosphere from below and are being lost to space from the top.

Like the Fe-peak elements, these rare earths are observed to have mean abundances that clearly decrease with stellar age. This could be because the radiative acceleration upwards is probably declining with decreasing T_eff, and thus either decreasing the equilibrium abundance in the atmosphere, or decreasing the supply at the bottom of the atmosphere if the flow-through case applies.

6. Summary and conclusions

In this paper, we report the discovery of the time variations of atmospheric abundances during the main sequence lifetime of a magnetic Bp star. Large overabundances are observed for Ti, Cr, and Fe near the ZAMS. As the star evolves on the main sequence, the abundances of these Fe-peak elements clearly decrease, approaching nearly solar values closer to the TAMS. The rare-earth elements Pr and Nd show drastic overabundances in young Bp stars, with values nearly 10⁴ times greater than in the Sun. Near the TAMS, magnetic Bp stars remain overabundant in rare-earth elements, still exhibiting abundances of Pr and Nd that are at least a factor of 100 larger than the solar ratios. The light elements O, Mg and Si show no evidence of time variations. However, O and Mg are, in general, underabundant and Si is always overabundant compared to the Sun. Remarkably, we found that the abundance of He increases during the main sequence lifetime from about 1% to 10% of the solar He abundance. We conclude that the observed increase is either the result of significant mixing in the stellar envelope of He that is not related to diffusion, or that there is accretion of He-rich gas from the interstellar medium.

As the climate in magnetic Bp stars changes with stellar age, important and systematic changes in surface chemistry result. These changes are due to both the evolving current climate of the outer layers of the star, and to the previous history of diffusion and competing effects. The systematic changes we have discovered present a major challenge to theoretical modelling. Efforts to reproduce theoretically the observed evolutionary changes should lead to greatly improved understanding of the physics at work in the envelopes of magnetic Bp stars. Specifically, more detailed calculations on the equilibrium abundances for elements are necessary to fully interpret the results we present. Further, efforts to follow the evolution of a 3.5 M_⊙ star with time from a series of equilibrium abundance models will help differentiate between competing effects that influence the amount of radiative support for atoms in the stellar atmosphere. These competing effects include, for example, the decrease in gravity and the probable decrease in radiative acceleration with decreasing T_eff.

Future work will increase the number of stars in the current mass bin to enhance the current sample as well as increase the number of mass bins studied to include more massive stars (4−5 M_⊙) and less massive stars (2−3 M_⊙) to see if similar trends are observed.

Acknowledgments

J.D.B. and J.D.L. are grateful for support by the Natural Sciences and Engineering Research Council of Canada. J.D.L. is pleased to acknowledge a Leverhulme Visiting Professorship at the Armagh Observatory. J.D.B. thanks Dr. Tanya Ryabchikova for generously sharing her improved data to help create Fig. 4.

References

All Tables

All Figures

	Fig. 1 Spectrum synthesis of the region 5000−5080 Å for two spectra of HIP 109911. The observed spectra are in black and the model fits are in red.
In the text

	Fig. 2 Atmospheric abundances of He, O, Mg, Si, Ti, Cr, Fe, Pr, and Nd plotted against log t. Each plot contains the solar abundance ratio for that element (solid red line) and the best-fit linear regression (dashed blue line).
In the text

	Fig. 3 Same as Fig. 2 but for Brms.
In the text

	Fig. 4 Shown are Cr (left) and Fe (right) abundances verses Teff. The black triangles are data taken from Ryabchikova et al. (2004b) and Ryabchikova (2005), the green squares are from Bailey & Landstreet (2013b) and the blue circles are data from this paper.
In the text

	Fig. 5 Position of the stars of this study in the Hertzsprung-Russell diagram, compared to the main sequence evolution tracks for stars with 3.0, 3.5 and 4.0 M⊙. Also shown are the ZAMS, the TAMS, and the points where half the main sequence lifetime is completed.
In the text