Äîêóìåíò âçÿò èç êýøà ïîèñêîâîé ìàøèíû. Àäðåñ îðèãèíàëüíîãî äîêóìåíòà : http://star.arm.ac.uk/preprints/2013/647.pdf Äàòà èçìåíåíèÿ: Wed Dec 4 17:47:51 2013 Äàòà èíäåêñèðîâàíèÿ: Fri Feb 28 00:55:04 2014 Êîäèðîâêà: Ïîèñêîâûå ñëîâà: http astrokuban.info astrokuban

Astronomy & Astrophysics manuscript no. Bailey-Landstreet-Bagnulo De c e m b e r 4 , 2 0 1 3

c ESO 2013

Discovery of secular variations in the atmospheric abundances of magnetic Ap stars
J. D. Bailey1 , J. D. Landstreet2,1 , and S. Bagnulo2
1 2

Department of Physics & Astronomy, The University of Western Ontario, London, Ontario, N6A 3K7, Canada Armagh Observatory, College Hill, Armagh, Northern Ireland BT61 9DG

Accepted 26 November 2013
ABSTRACT

Context. The stars of the middle main sequence have relatively quiescent outer layers, and unusual chemical abundance patterns may develop in their atmospheres. The presence of chemical peculiarities reveal the action of such subsurface phenomena as gravitational settling and radiatively driven levitation of trace elements, and their competition with mixing processes such as turbulent diffusion. At present, little is known about the time evolution of these anomalous abundances, nor about the role that diffusion may play in maintaining them, during the main sequence lifetime of such a star. Aims. We want to establish whether abundance peculiarities change as stars evolve on the main sequence, and provide observational constraints to diffusion theory. Methods. We have performed spectral analysis of 15 magnetic Bp stars that are members of open clusters (and thus have well-known ages), with masses between about 3 and 4 M . For each star, we measured the abundances of He, O, Mg, Si, Ti, Cr, Fe, Pr and Nd. Results. We have discovered the systematic time evolution of trace elements through the main-sequence lifetime of magnetic chemically peculiar stars as their atmospheres cool and evolve toward lower gravity. During the main sequence lifetime, we observe clear and systematic variations in the atmospheric abundances of He, Ti, Cr, Fe, Pr and Nd. For all these elements, except He, the atmospheric abundances decrease with age. The abundances of Fe-peak elements converge toward solar values, while the rare-earth elements converge toward values at least 100 times more abundant than in the Sun. Helium is always underabundant compared to the Sun, evolving from about 1% up to 10% of the solar He abundance. We have attempted to interpret the observed abundance variations in the context of radiatively driven diffusion theory, which appears to provide a framework to understand some, but not all, of the anomalous abundance levels and variations that we observe.
Key words. stars:magnetic field - stars: chemically peculiar - stars: abundances - stars: atmospheres - stars: evolution

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 outward radiation flux drives a strong stellar wind, stripping maBased in part on observations made with the European Southern Observatory (ESO) telescopes under the ESO programmes 072.D0410(A) and 086.D-0449(A). It is also based in part on observations carried out at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council of Canada, the ´ Institut National des Science de lUniverse of the Centre National de la Recherche Scientifique of France and the University of Hawaii.

terial 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 downward 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


2

J. D. Bailey et al.: Abundance variations in magnetic Ap stars Table 1. Stars analysed in this study. Listed are the star designations, instrument used (with the number of spectra indicated in parentheses), spectral resolution, and spectral range. Star HD HD HD HD HD Instrument FEROS (2) PaDOnS (2) FEROS (2) FEROS (2) PaDOnS (1) FEROS (1) PaDOnS (2) PaDOnS (2) PaDOnS (2) PaDOnS (2) FEROS (2) UVES (1) PaDOnS (1) PaDOnS (2) PaDOnS (2) PaDOnS (2) R 48000 65000 48000 48000 65000 48000 65000 65000 65000 65000 48000 110000 65000 65000 65000 65000 ( å) 3528-9217 3690-10481 3528-9217 3528-9217 3690-10481 3528-9217 3690-10481 3690-10481 3690-10481 3690-10481 3528-9217 3070-10398 3690-10481 3690-10481 3690-10481 3690-10481

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 NX /NH 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 upward and downward diffusion (Landstreet 2004). It is therefore of obvious importance to observe and characterise the nature and time variations of these p h en o m en a. 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 108 - 109 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. Sect. 3 describes the modelling technique. Sect. 4 details the observational results. Sect. 5 explores possible mechanisms to explain the observed trends and Sect. 6 summarises the main conclusions.

45583 61045 63401 74535 133652

ES ES E E E E E E E E S S S S S S S S

HD 133880 HD 147010 HD 162576 HD 162725 HD 304842 NGC 2169 12 BD+00 1659 BD-19 5044L BD+49 3789 HIP 109911

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 ab o u t 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 an 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 Stromgren ¨ 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 v sin i. The final column shows the root-mean-square magnetic field strength ( Brms ) as computed from whatever modern measurements of mean longitudinal magnetic field Bz 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 13500 K and 4.4 to 10000 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


J. D. Bailey et al.: Abundance variations in magnetic Ap stars Table 2. Physical properties of the stars studied. Star Cluster HD 147010 Upper Sco NGC 2169 12 NGC 2169 HD 133652 Upper Cen Lup HD 133880 HD 45583 NGC 2232 HD 63401 NGC 2451 HD 74535 IC 2391 BD-19 5044L IC 4725 BD+49 3789 NGC 7243 HIP 109911 HD 61045 NGC 2422 HD 304842 NGC 3114 BD+00 1659 NGC 2301 HD 162576 NGC 6475 HD 162725 references ­ (1) Landstreet et al. ( log t 6.70 ± 0.101 6.97 ± 0.101 7.20 ± 0.101 7. 7. 7. 8. 8. 8. 8. 8. 8. 55 70 70 02 06 08 13 22 41 ± ± ± ± ± ± ± ± ± 0. 0. 0. 0. 0. 0. 0. 0. 0. 101 101 151 081 10 111 151 10 131 M/M 3.15 ± 0.201 3.65 ± 0.151 3.35 ± 0.151 3.20 ± 0.151 3.30 ± 0.151 3.70 ± 0.201 3.85 ± 0.151 3.55 ± 0.151 3.55 ± 0.15 3.65 ± 0.15 3.85 ± 0.201 3.55 ± 0.151 3.65 ± 0.15 3.10 ± 0.15 3.30 ± 0.201 ey & Landstreet ( T eff (K) 13000 ± 5002 13800 ± 5001 13000 ± 500 13000 ± 6003 12700 ± 5002 13500 ± 5001 13600 ± 5001 12800 ± 5002 12900 ± 5002 13000 ± 500 13000 ± 5002 12500 ± 5002 12500 ± 5002 10300 ± 500 10000 ± 500 2013b); (3) Bail log g v sin i (km s-1 ) 4.40 ± 0.202 15 ± 22 4.30 ± 0.20 56 ± 5 4.30 ± 0.20 48 ± 2 4.34 ± 0.163 103 ± 103 2 4.20 ± 0.20 70 ± 62 4.20 ± 0.20 52 ± 4 4.30 ± 0.20 45 ± 4 4.50 ± 0.202 15 ± 32 2 4.20 ± 0.20 85 ± 52 4.30 ± 0.20 60 ± 2 4.10 ± 0.202 64 ± 32 3.90 ± 0.202 65 ± 52 4.00 ± 0.202 7.0 ± 12 3.70 ± 0.20 28 ± 3 3.50 ± 0.20 31 ± 3 ey et al. (2012); (*) This work B
r ms

3

( G) 48251 34101 11201 23001 27301 3651 951 2351 561 348 4301 201 394 15 69

2007); (2) Bail

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 pretabulated 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. 2000; Ryabchikova et al. 1997; Piskunov et al. 1995; Kupka et al. 1999). 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 11000 and 14000 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. zeeman searches for an optimal fit between the synthetic and observed spectra by means of a reduced 2 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 vR and v sin i. The abundance (log NX /NH ) 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 10000 to 14000 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

Table 3. List of spectral lines modelled. Element He i 4437. 4713. 5015. 5047. 5875. 5875. 5875. 5875. 5875. 5875. 6155. 6155. 6155. 6156. 6156. 6156. 6158. 6158. 7771. 7774. 7775. 4481. 4481. 4481. 4621. 4621. 5041. ( å) 551 139 678 738 599 614 615 625 640 966 961 966 986 736 755 776 146 172 941 161 388 126 150 325 418 722 024 Element Si ii Ti ii Cr ii 5055. 5056. 4533. 4563. 4571. 4558. 4565. 4588. 4592. 4541. 4555. 4556. 4583. 4583. 5029. 5030. 5032. 5035. 6160. 6161. 7781. 4911. 4912. 4914. 5050. 6145. ( å) 984 317 960 757 968 650 739 199 049 524 893 392 837 999 097 630 712 708 233 194 983 653 944 094 695 068

Fe ii

Oi

Pr iii Nd iii

Mg ii Si ii

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 Bz , 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,


4

J. D. Bailey et al.: Abundance variations in magnetic Ap stars

1

Intensity

Fe II He I Fe II

Ti II

Fe II Fe II

Fe II

0.8

Si II

0.6 5000

5020

5040 Wavelength (å)

Si II

5060

Fe II Fe II Ti II Fe II 5080

He I

Fe II

Fe II

Fe II Fe II

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.

Fe II

we simply assumed a dipolar magnetic field that was approximately three times Brms 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(NX /NH ) 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 12000 K, and shown in their Figure 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) sug-

gest 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.

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 NX /NH 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

Fe II


J. D. Bailey et al.: Abundance variations in magnetic Ap stars

Table 4. Average derived abundances for the stars studied. Star HD 147010 NGC 2169 12 HD 133652 HD 133880 HD 45583 HD 63401 HD 74535 BD-19 5044L BD+49 3789 HIP 109911 HD 61045 HD 304842 BD+00 1659 HD 162576 HD 162725 Sun references ­ (1) log t 6.70 ± 0.101 6.97 ± 0.101 7.20 ± 0.101 7. 7. 7. 8. 8. 8. 8. 8. 8. 55 70 70 02 06 08 13 22 41 ± ± ± ± ± ± ± ± ± 0. 0. 0. 0. 0. 0. 0. 0. 0. 101 101 151 081 10 111 151 10 131 log(He/H) -3.11 ± 0.28 -3.09 ± 0.30 -3.16 ± 0.28 -2.00 -2.40 -2.64 ± 0.28 -2.52 ± 0.28 -2.16 ± 0.21 -2.52 ± 0.35 -2.34 ± 0.28 -2.06 ± 0.28 -1.54 ± 0.28 -2.34 ± 0.15 -1.51 ± 0.28 -1.74 ± 0.36 -1.07 al. (2007); (2) Bailey & Landstreet 5002 5001 500 6003 5002 5001 5001 5002 5002 600 5002 5002 5002 500 500 T eff (K 13000 13800 13000 13000 12700 13500 13600 12800 12900 13000 13000 12500 12500 10300 10000 ) ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± log(O/H -4.49 ± -3.88 ± -3.87 ± -2.90 ± -3.61 ± -3.73 ± -4.26 ± -3.99 ± -3.60 ± -3.69 ± -3.72 ± -3.78 ± -3.58 ± -3.77 ± -3.77 ± -3.31 (2013b); ) 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. log(Mg/H) -5.67 ± 0.14 -5.56 ± 0.15 -5.18 ± 0.14 -4.12 ± 0.28 -4.19 ± 0.28 -5.86 ± 0.14 -5.44 ± 0.28 -5.44 ± 0.21 -5.25 ± 0.25 -5.38 ± 0.14 -5.10 ± 0.14 -5.76 ± 0.21 -5.62 ± 0.10 -5.34 ± 0.14 -5.23 ± 0.14 -4.40 (3) Bailey et al. (2012) 28 15 28 42 42 21 28 28 22 21 28 25 10 28 28 log(Si ii/H) -3.72 ± 0.14 -3.87 ± 0.15 -3.36 ± 0.28 -2.94 ± 0.28 -3.34 ± 0.28 -3.96 ± 0.14 -4.28 ± 0.14 -3.94 ± 0.28 -3.54 ± 0.14 -3.49 ± 0.14 -3.94 ± 0.28 -3.64 ± 0.28 -3.53 ± 0.20 -4.44 ± 0.14 -3.54 ± 0.28 -4.49 ; (*) This work log(Ti -5.68 -5.90 -5.01 -5.44 -5.68 -6.44 -6.29 -6.82 -6.33 -5.90 -6.24 -6.83 -6.69 -7.72 -7.30 -7.05 /H ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ) 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. 0. log(Cr/H) -4.04 ± 0. -4.80 ± 0. -4.17 ± 0. -4.65 ± 0. -4.74 ± 0. -5.32 ± 0. -5.65 ± 0. -5.88 ± 0. -5.44 ± 0. -5.35 ± 0. -5.47 ± 0. -6.05 ± 0. -5.02 ± 0. -5.12 ± 0. -4.64 ± 0. -6.36 log(Fe/H) -3.33 ± 0. -3.76 ± 0. -2.90 ± 0. -3.56 ± 0. -3.43 ± 0. -3.84 ± 0. -4.02 ± 0. -4.24 ± 0. -4.06 ± 0. -3.59 ± 0. -3.98 ± 0. -4.35 ± 0. -3.79 ± 0. -3.95 ± 0. -3.66 ± 0. -4.50 log(Pr/H) -6.15 ± 0.28 -7.31 ± 0.20 -6.98 ± 0.28 -6.96 ± 0.28 -7.49 ± 0.42 -7.14 ± 0.28 -8.15 ± 0.21 -8.12 ± 0.28 -7.72 ± 0.28 -7.70 ± 0.28 -7.74 ± 0.28 -7.65 ± 0.36 -8.11 ± 0.20 -10.15 ± 0.42 -8.32 ± 0.42 -11.28 log(Nd/H) -6.50 ± 0. -7.21 ± 0. -6.42 ± 0. -7.04 ± 0. -7.04 ± 0. -7.44 ± 0. -7.52 ± 0. -7.88 ± 0. -7.72 ± 0. -7.37 ± 0. -7.77 ± 0. -7.04 ± 0. -7.34 ± 0. -9.10 ± 0. -7.44 ± 0. -10.58

21 20 21 28 28 28 28 28 35 28 28 14 10 14 28

14 20 21 21 21 28 14 28 28 28 28 32 10 14 22

14 10 14 14 21 14 28 28 28 28 28 32 10 14 28

28 20 28 42 42 28 28 42 25 28 42 28 20 42 28

Landstreet et

5


6

J. D. Bailey et al.: Abundance variations in magnetic Ap stars

Table 5. The linear fit parameters for the abundance of each element versus age and magnetic field strength. Shown are the slopes with their respective uncertainties and the significance of the slope, D y-axis log(He/H) log(O/H) log(Mg/H) log(Si/H) log(Ti/H) log(Cr/H) log(Fe/H) log(Pr/H) log(Nd/H) x-axis log t log Brm log t log Brm log t log Brm log t log Brm log t log Brm log t log Brm log t log Brm log t log Brm log t log Brm Slope ± 0.18 ± 0.11 ± 0.17 ± 0.13 ± 0.26 ± 0.18 ± 0.20 ± 0.13 ± 0.27 ± 0.12 ± 0.24 ± 0.16 ± 0.16 ± 0.08 ± 0.29 ± 0.20 ± 0.24 ± 0.17 D 4.2 4.4 0.76 0.31 0.58 1.1 0.65 1.5 4.0 5.8 4.0 2.5 2.6 3.0 4.4 4.4 3.4 2.7

out the main sequence evolution. Moderate Si overabundance is common for magnetic Bp stars between effective temperatures o f a b o u t 1 0 0 0 0 a n d 1 5 0 0 0 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

s s

- -

s

-
s

-
s

-
s

-
s

-
s

-
s

0. 0. 0. 0. 0. 0. 0. 0. 1. 0. 0. 0. 0. 0. 1. 0. 0. 0.

76 48 13 04 15 20 13 19 07 69 97 40 41 24 29 88 81 46

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 that the solar values.
4.1.4. The rare ear th elements: Pr and Nd

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 o f ab o u t 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 . He l i u m

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.
4.2. Abundance variations with magnetic field strength

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 welldocumented 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 through-

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. Fig. 3 plots abundances versus log Brms 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 Brms , 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 Figure 3 are those expected from the known decline of Bz 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 Figure 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; Kudryavtsev & Romanyuk 2008; Kudryavtsev et al. 2011) to search for particularly large magnetic fields in Ap/Bp stars.


J. D. Bailey et al.: Abundance variations in magnetic Ap stars
-1 -1.5 -2.5 -4 -3

7

log (Mg/H) log (Cr/H) log (Nd/H)
7

log (He/H)

-2 -2.5 -3 -3.5

log (O/H)

-4.5 -5 -5.5

-3.5 -4 -4.5

-6

-4 -3 -5 -5.5 -4.5 -5 -5.5 -6

log (Ti/H)

log (Si/H)

-3.5

-6 -6.5 -7 -7.5

-4

-4.5

-3

-3.5

-4

-4.5 6.5 7

-6 -6.5 -7 -7.5 -8 -8.5 -9 -9.5 -10 -10.5 -11 6.5

-6.5 -7 -7.5 -8 -8.5 -9 -9.5 -10 7.5 8 8.5 -10.5 6.5 7 7.5 8 8.5

log (Fe/H)

log t (yrs)

7.5

8

8.5

log (Pr/H)

log t (yrs)

log t (yrs)

Fig. 2. Shown are the 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.3. Cr and Fe abundances versus effective temperature

that the maximum abundance for Cr and Fe both approach the same value of between -3 and -4 dex.

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 10000 K. Ryabchikova (2005) added stars between about 10000 and 15000 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 10500 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)

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 a Hertzsprung-Russell (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


8
-1 -1.5

J. D. Bailey et al.: Abundance variations in magnetic Ap stars
-2.5 -4 -3

log (Mg/H) log (Cr/H) log (Nd/H)
2 3 4

log (He/H)

-2 -2.5 -3 -3.5

log (O/H)

-4.5 -5 -5.5

-3.5 -4 -4.5

-6

-4 -3 -5 -5.5 -4.5 -5 -5.5 -6

log (Ti/H)

log (Si/H)

-3.5

-6 -6.5 -7 -7.5

-4

-4.5

-3

-3.5

-4

-4.5 1 1.5 2

-6 -6.5 -7 -7.5 -8 -8.5 -9 -9.5 -10 -10.5 -11 1 1.5

-6.5 -7 -7.5 -8 -8.5 -9 -9.5 -10 2.5 3.5 -10.5 1 1.5 2 2.5 3 3.5 4

log (Fe/H)

log Brms (G)

2.5

3

3.5

4

log (Pr/H)

log Brms (G)

log Brms (G)

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

r ms

.

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 downward relative to the dominant hydrogen gas, but the outward 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 abund an ce. Detailed and quantitative modelling of diffusion is actually quite complicated. The average upward 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 upward acceleration gR due to the radiation field is larger than the local acceleration of gravity g, and the ion diffuses upwards rather than downward. 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 gR = 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


J. D. Bailey et al.: Abundance variations in magnetic Ap stars

9

-2.5 -3 -3.5 -4 log (Cr/H) -4.5 -5 -5.5 -6 -6.5 -7 log (Fe/H) 6000 8000 10000 12000 14000 Effective Temperature (K)

-2.5 -3 -3.5 -4 -4.5 -5 -5.5 -6 -6.5 -7

6000 8000 10000 12000 14000 Effective Temperature (K)

Fig. 4. Shown are Cr (left) and Fe (right) abundances verses T eff . 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.

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 (2 0 1 2 ). 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 atmo-

spheres 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 will 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 Figure 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 12000 K is much too weak to support a solar abundance of these elements (cf. grad calculations for T eff = 12000 K by Hui-Bon-Hoa et al. 2002). These two elements are expected to diffuse downward 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


10

J. D. Bailey et al.: Abundance variations in magnetic Ap stars

3

4.0

ZAMS TAMS 1/2 MS

log (L/L )

3.5

3.0

2

4.2

4.1

4

3.9

log (Teff)
Fig. 5. The 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.0M . Also shown are the ZAMS, the TAMS, and the points where half the main sequence lifetime is completed.

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 & Kubat (2004) and Theado et al. (2005); ´ ´ all groups confirm that without a mixed wind from the star, that He should diffuse downward 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 (Figure 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 upward 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 ap-

pears 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 grad 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 gR . However, the decrease in log g means that a smaller value of gR 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.5 M from T eff = 13000 K and log g = 4.5 to T eff = 10000 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.


J. D. Bailey et al.: Abundance variations in magnetic Ap stars

11

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(nMg /nH ) -5.2. This may be compared with the predicted Mg abundance profile for a star of T eff = 12000 K of Leblanc et al. (2009) and of Alecian & Stift (2010). The abundance predicted by the equilibrium calculation (gR = g throughout the atmosphere) is mildly non-uniform through the line-forming region, but is of about this same values. 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 12000 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 longstanding 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 (Figure 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 infor-

mation 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(nel /nH ) 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 flowthrough 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 = 12000 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.


12

J. D. Bailey et al.: Abundance variations in magnetic Ap stars

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 upward 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.

References
Alecian, G. & Stift, M. J. 2010, A&A, 516, A53 Alecian, G., Stift, M. J., & Dorfi, E. A. 2011, MNRAS, 418, 986 Alecian, G. & Vauclair, S. 1981, A&A, 101, 16 Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 Bagnulo, S., Landstreet, J. D., Mason, E., et al. 2006, A&A, 450, 777 Bagnulo, S., Wade, G. A., Donati, J.-F., et al. 2001, A&A, 369, 889 Bailey, J. D., Grunhut, J., Shultz, M., et al. 2012, MNRAS, 423, 328 Bailey, J. D. & Landstreet, J., D. 2013a, MNRAS Bailey, J. D. & Landstreet, J. D. 2013b, A&A, 551, A30 Bailey, J. D., Landstreet, J. D., Bagnulo, S., et al. 2011, A&A, 535, A25 Cramer, N. & Maeder, A. 1980, A&A, 88, 135 Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&A, 141, 371 Havnes, O. 1974, A&A, 32, 161 Havnes, O. & Conti, P. S. 1971, A&A, 14, 1 Hui-Bon-Hoa, A., LeBlanc, F., Hauschildt, P. H., & Baron, E. 2002, A&A, 381, 197 Kamp, I. & Paunzen, E. 2002, MNRAS, 335, L45 Kochukhov, O. & Bagnulo, S. 2006, A&A, 450, 763 Kochukhov, O., Bagnulo, S., Wade, G. A., et al. 2004, A&A, 414, 613 Kochukhov, O. & Wade, G. A. 2010, A&A, 513, A13 Krticka, J. & Kubat, J. 2004, in IAU Symposium, Vol. 224, The A-Star Puzzle, ´ ed. J. Zverko, J. Ziznovsky, S. J. Adelman, & W. W. Weiss, 201­207 Kudryavtsev, D. O. & Romanyuk, I. I. 2008, Contributions of the Astronomical Observatory Skalnate Pleso, 38, 229 Kudryavtsev, D. O., Romanyuk, I. I., Elkin, V. G., & Paunzen, E. 2006, MNRAS, 372, 1804 Kudryavtsev, D. O., Romanyuk, I. I., Semenko, E. A., & Yakunin, I. A. 2011, Astronomische Nachrichten, 332, 961 Kunzli, M., North, P., Kurucz, R. L., & Nicolet, B. 1997, A&AS, 122, 51 Kupka, F., Piskunov, N. E., Ryabchikova, T. A., Stempels, N. C., & Weiss, W. W. 1999, A&A, 138, 119 Kupka, F. G., Ryabchikova, T. A., Piskunov, N. E., Stempels, H. C., & Weiss, W. W. 2000, Baltic Astronomy, 9, 590 Landstreet, J. D. 1988, ApJ, 326, 967 Landstreet, J. D. 2004, in IAU Symposium, Vol. 224, The A-Star Puzzle, ed. J. Zverko, J. Ziznovsky, S. J. Adelman, & W. W. Weiss, 423­432 Landstreet, J. D., Bagnulo, S., Andretta, V., et al. 2007, A&A, 470, 685 Landstreet, J. D., Dolez, N., & Vauclair, S. 1998, A&A, 333, 977 Landstreet, J. D., Kupka, F., Ford, H. A., et al. 2009, A&A, 503, 973 Landstreet, J. D., Silaj, J., Andretta, V., et al. 2008, A&A, 481, 465 LeBlanc, F. 2003, in Astronomical Society of the Pacific Conference Series, Vol. 305, Magnetic Fields in O, B and A Stars: Origin and Connection to Pulsation, Rotation and Mass Loss, ed. L. A. Balona, H. F. Henrichs, & R. Medupe, 206 Leblanc, F. & Monin, D. 2004, in IAU Symposium, Vol. 224, The A-Star Puzzle, ed. J. Zverko, J. Ziznovsky, S. J. Adelman, & W. W. Weiss, 193­200 Leblanc, F., Monin, D., Hui-Bon-Hoa, A., & Hauschildt, P. H. 2009, A&A, 495, 937 Mashonkina, L. I., Ryabchikova, T. A., & Ryabtsev, A. N. 2004, in IAU Symposium, Vol. 224, The A-Star Puzzle, ed. J. Zverko, J. Ziznovsky, S. J. Adelman, & W. W. Weiss, 315­322 Meynet, G., Mermilliod, J.-C., & Maeder, A. 1993, A&AS, 98, 477 Michaud, G., Charland, Y., Vauclair, S., & Vauclair, G. 1976, ApJ, 210, 447 Michaud, G., Richer, J., & Vick, M. 2011, A&A, 534, A18 Napiwotzki, R., Schoenberner, D., & Wenske, V. 1993, A&A, 268, 653 Nesvacil, N., Shulyak, D., Ryabchikova, T. A., et al. 2013, A&A, 552, A28 Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&A, 112, 525 Ryabchikova, T., Leone, F., & Kochukhov, O. 2005, A&A, 438, 973 Ryabchikova, T., Leone, F., Kochukhov, O., & Bagnulo, S. 2004a, in IAU Symposium, Vol. 224, The A-Star Puzzle, ed. J. Zverko, J. Ziznovsky, S. J. Adelman, & W. W. Weiss, 580­586 Ryabchikova, T., Nesvacil, N., Weiss, W. W., Kochukhov, O., & Stutz, C. 2004b, ¨ A&A, 423, 705 Ryabchikova, T. A. 2005, Astronomy Letters, 31, 388 Ryabchikova, T. A., Piskunov, N. E., Kupka, F., & Weiss, W. W. 1997, Baltic Astronomy, 6, 244 Shulyak, D., Ryabchikova, T., Mashonkina, L., & Kochukhov, O. 2009, A&A, 499, 879 Stift, M. J. & Alecian, G. 2012, MNRAS, 425, 2715 Theado, S., Vauclair, S., & Cunha, M. S. 2005, A&A, 443, 627 ´ Vauclair, S. 1983, in Saas-Fee Advanced Course 13: Astrophysical Processes in Upper Main Sequence Stars, ed. A. N. Cox, S. Vauclair, & J. P. Zahn, 167 Vauclair, S., Dolez, N., & Gough, D. O. 1991, A&A, 252, 618 Vick, M., Michaud, G., Richer, J., & Richard, O. 2010, A&A, 521, A62 Wade, G. A., Bagnulo, S., Kochukhov, O., et al. 2001a, A&A, 374, 265 Wade, G. A., Ryabchikova, T. A., Bagnulo, S., & Piskunov, N. 2001b, in

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 104 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 (45 M ) and less massive stars (2-3 M ) to see if similar trends are observed.
Acknowledgements. JDB and JDL are grateful for support by the Natural Sciences and Engineering Research Council of Canada. JDL is pleased to acknowledge a Leverhulme Visiting Professorship at the Armagh Observatory. JDB thanks Dr. Tanya Ryabchikova for generously sharing her improved data to help create Figure 4.


J. D. Bailey et al.: Abundance variations in magnetic Ap stars
Astronomical Society of the Pacific Conference Series, Vol. 248, Magnetic Fields Across the Hertzsprung-Russell Diagram, ed. G. Mathys, S. K. Solanki, & D. T. Wickramasinghe, 373

13