abstract: The NMR technique is quite essential as it permits atomic scale measurements in materials. The aim of this article is to present the main parameters accessible to NMR experiments and to relate them to the electronic properties of simple solids in which electronic correlations are not dominant. It will be shown that in metals the Pauli paramagnetism, that is the electronic density of states is accessible through NMR shift and spin lattice relaxation data. Those permit as well to determine electronic gaps in some electronic structures, or in the SC state. The capability to reveal by NMR the incidence of impurities and disorder is also underlined. A complementary article will be dedicated to the specific properties of strongly correlated materials which have been revealed and studied by NMR experiments.
author: Henri Alloul [^1]
title: **NMR studies of electronic properties of solids**

# Introduction

The electronic properties of solids were on the first half of the twentieth century considered mostly in the frame of an independent electron approximation with spin degeneracy. Their resulting electronic band structure is such that each electronic level could be doubly occupied. In such an approach one expects metals or insulators with no significant magnetic properties.

The traditional experimental studies of the electronic states in such solids usually require a determination of their electronic band structure and how it affects the physical properties. Those are usually determined by measurements taken at the macroscopic scale such as optical, transport and magnetic data. To go beyond these approaches the Nuclear Magnetic Resonance is thoroughly used. This technique, discovered in the mid of the 20th century is quite essential as it permits to perform atomic scale measurements in the materials, differentiating the properties which can be attributed to specific phases or sites in the structure.

Let us recall that the physical properties of solids, and in particular their magnetic properties, are determined by the electronic states. On the other hand, the nuclear spin moments, which do not affect these properties, provide an extremely useful probe for the electronic properties. Atomic nuclei are made up of neutrons and protons, which are spin 1/2 particles. They are assembled into quantum states in which the nuclear ground state has a total spin $\vec{I}$ that may be integer or half-integer. The associated magnetic moment ${\bf \mu }_{\mathrm{n}}$ is proportional to the magnetic moment of the proton $\mu _{\mathrm{p}}$, where the multiplicative factor is analogous to the Landé factor for an atomic electronic moment. Each atomic nucleus thus has a specific magnetic moment ${\bf \mu }_{\mathrm{n}}=\hbar \gamma _{\mathrm{n}}\vec{I}$. The gyromagnetic ratio $\gamma _{\mathrm{n}}$ is known to great accuracy for each of the stable isotopes in the periodic table.

Since $\mu _{\mathrm{p}}\simeq 10^{-3}\mu _{\mathrm{B}}$, the nuclear moments are extremely small, as are their mutual interactions. As a consequence, they are almost always in a paramagnetic state with a Curie magnetization ${\bf \mu _{z}=N(\hbar {\gamma _\mathrm{n} })^2 I(I+1)}B_{0}/3k_{B}T$. We see immediately that, in a given applied field, the nuclear magnetisation is about $10^{6}$ times smaller than the electronic magnetisation. The nuclear magnetism is practically impossible to detect using just macroscopic magnetization measurements taken on bulk samples. But, although the nuclear magnetic susceptibilities are weak, they can be detected by magnetic resonance which is a spectroscopy that permits a selective detection of the nuclear spin response. The weakness of the interaction between electronic and nuclear spins permits to consider the nuclear spins as somewhat ideal probes of the electronic properties of materials.

As the NMR technique is described at length in various articles in Wikipedia and Scholarpedia for its applications in chemistry, biology and medical sciences, we shall only recall below very briefly in section the technical principles of NMR, but will describe in more detail in section the couplings between nuclear spins and electron spins. Those allow one to probe the electronic properties of solids and their magnetic response through the determination of the shift of the NMR signal (section). In section , we recall the spin lattice relaxation processes which permit to use the dynamics of the nuclear magnetisation to probe the excited electronic states in condensed matter. The pairing of electrons in the superconducting state of simple metals also yield remarkable effects on the NMR (section ). We shall recall then in section that while spin $1/2$ nuclei only probe electronic spin excitations, nuclear spins probes with $I>1/2$ sense as well an electrostatic quadrupole interaction which affects the nuclear spin levels and gives access to charge ordering in materials. Finally, in section , we evidence the local effects induced by impurities or disorder on the NMR spectra.

# The magnetic resonance phenomenon and some of its basic applications

Discovered at the outset of the second world War by F. Bloch and E.M. Purcell, the nuclear magnetic resonance technique has become quite immediately a unique method to investigate the chemical and physical properties of condensed matter (Abragam ,1961, Slichter ,1963). Its success results from the fact that it resolves spectroscopicaly the properties of the nuclear spins of the distinct atomic species present in materials. A homogeneous applied external magnetic field $B_{0}$ induces a splitting of the nuclear spin energy level $h\nu _{L}=\hbar \gamma _{\mathrm{n}}B_{0}$ which usually falls in the radio-frequency range of the order of 10 MHz per Tesla. The absorption of a radio-frequency field at the adequate frequency $\nu _{L}$ permits to detect the presence of the corresponding nuclei.

One highlight of NMR is the acquired possibility to provide images of the spatial distribution of $^{1}$H nuclei in-vivo in biological matter, which is the basis of medical Magnetic Resonance Imaging. Though this is the popular application of NMR known from a large audience, NMR is an even more powerful technique when one uses the interactions of given nuclear spins in a material with their neighboring atomic states. This results in rich spectroscopic splittings of the NMR lines which permit to locate the atoms in molecular states and therefore to determine the molecular structures in chemistry or in the solid state ("Nuclear Magnetic Resonance" in Wikipedia). Such spectroscopic techniques have been revealed since the 1980s but their impact has been tremendously highlighted by the improvement of the SC magnet industry which has allowed to produce extremely homogeneous magnetic fields $B_{0}$ as large as 21 Tesla with negligible drift in time. In the corresponding range of frequencies $\nu _{L}$, exceptionally stable coherent sources are available, with narrower spectral widths than the transitions to be observed. These are obtained by electronic oscillators with frequency stabilized on the resonant mode of a piezoelectric quartz crystal. Moreover, for such frequencies, very powerful amplifiers are also available. With radio-frequency pulses, these can be used to significantly modify the populations of the spin quantum states. Such setups are useful for studying nuclear relaxation. Therefore on the electronics side, NMR has highly benefited of all the developments of fast semi-conductor and computing capabilities associated with the expansion of information technologies.

#  Electronic hyperfine couplings

It is clear that changes in the magnetic induction in a material can be detected directly by a change in the nuclear Larmor frequency. In weakly magnetic materials, for which the magnetization is negligible, $B_{0}=\mu
_{0}H_{a}$ and the nuclear Larmor frequency $\nu _{L}$ should be determined solely by the applied external field $H_{a}$. It would be difficult to obtain information about the physical properties of materials in such a limit. But we have to recall that the nucleus is a kind of atomic scale microscopic probe, coupled to the electrons. Interactions like the dipole interactions between nuclear and electronic spins are such that the nuclear spin feels a magnetic field associated with the polarization of the electronic magnetic moments. This means that the magnetic field felt by the nuclear spins is modified with respect to the applied field. It is the spectroscopy of these fields that provides atomic scale information about the immediate vicinity of the nuclei in the material. Let us examine the different interactions between the nuclear spins and the magnetic moments of electronic origins, known collectively as hyperfine interactions (Abragam 1961, Slichter 1963).

The dipole interaction between the moments associated with a nuclear spin $\mathbf{I}$ and an electron spin $\mathbf{S}$ separated by a displacement $\mathbf{r}$ is

$$\label{eq:Hdd}
H_{dd}=-\frac{\mu_{0}}{4\pi }
\frac{\gamma _{e}\gamma _{n}\ \hbar ^{2}}{r^{3}}
\left[ \mathbf{I}.\mathbf{S}-\frac{3(\mathbf{I}.\mathbf{r})(\mathbf{S}.\mathbf{r})}{r^{2}}\right]$$

where $\gamma _{n}$ and $\gamma _{e}$ are the nuclear and electronic gyromagnetic moments, respectively, and $\mathbf{S}$ and$\mathbf{\ I}$ are here dimensionless quantities. This dipole interaction diverges when $r$ tends to zero, and is only therefore valid for electrons with zero probability of being at the site of the nucleus. This is the case for electrons in the $p,d$, or $f$ shells. On the other hand, the$\ s$ electrons have nonzero probability of being at the site of the nucleus. The Dirac Hamiltonian can be used to show that the corresponding interaction, called the contact interaction $H_{c}$ is scalar, and is given in this case by $$\label{eq:H Contact}
H_{c}=\frac{%
%TCIMACRO{\U{3bc} }%
%BeginExpansion
\mu
%EndExpansion
_{0}}{4\pi}8\pi ^{3}\ \gamma _{e}\gamma _{n}\ \hbar ^{2\ }\ \mathbf{%
I}.\mathbf{S\ }\delta (\mathbf{r}).$$

Finally, the interaction with the magnetic field associated with the orbital angular momentum of the electron is $$\label{eq:Horb}
H_{orb}=-\frac{%
%TCIMACRO{\U{3bc} }%
%BeginExpansion
\mu
%EndExpansion
_{0}}{4\pi }\gamma _{e}\gamma _{n}\  \hbar ^{2\ }\ \frac{\mathbf{I}.%
\mathbf{\ell }}{r^{3}}\mathbf{\ }.$$

These Hamiltonians can all be written in the form $$\label{eq:Heff}
H_{eff}=-\gamma _{n}\ \hbar \mathbf{\ \ I}.\mathbf{B}_{eff},$$

and we may consider that each electron induces a magnetic field $\mathbf{B}_{eff}$ at the nuclear site. As the temporal fluctuations of the electronic moments are very fast compared with the nuclear Larmor frequency, the static component of $\mathbf{B}_{eff}$ is its time average. The position of the NMR for a given nucleus is thus determined by the time average $\mathbf{B}_{eff}$ of the resultant of the fields due to the different electrons in the material. It is easy to see that the hyperfine interaction will vanish for filled electronic shells, because they have zero total spin and total orbital angular momentum. When there is no applied field, $\mathbf{B}_{eff}$ can only be nonzero for materials in which there is a static spin or orbital magnetic moment. This will be the case for magnetically ordered materials.

# NMR shifts

In a material, each of the distinct hyperfine couplings listed above induces a specific contribution to the NMR shift. Those are of course quite dependent on the magnitude of the corresponding hyperfine coupling and of the electronic state of the considered material.

## Chemical shifts

In substances where the electrons are paired in atomic or molecular orbitals, the static part of the hyperfine coupling is only nonzero in the presence of an applied field $\mathbf{B}_{0}$, and is proportional to $\mathbf{B}_{0}$, like the magnetization. The resonance is shifted with respect to that of the free atom in a gas. The relative shift $\mathbf{B}_{eff}/\mathbf{B}_{0}$ may be due to the orbital part of the hyperfine coupling. This is the case, for example, for the displacement due to the orbital currents induced by the external magnetic field in electronic or molecular shells close to the nucleus. Since this shift depends on the electronic charge distributions, it is highly sensitive to the chemical environment of the given atom, hence the name chemical shift. These effects are generally small, and expressed in in parts per million (ppm) but can be used to distinguish the nuclear spin resonances of the different atoms depending on their environment. This has become a very powerful tool, used universally in chemistry and biology. Routine chemical analyses are carried out by NMR. It also helps one to determine the 3D structures of biological molecules, using multidimensional methods, which have reached an exceedingly high level of refinement.

## Knight shifts in Metals

When the electron states are not paired in molecular states or in bonds, the spin degeneracy of the electronic states might be lifted by the applied field, as for electron states at the Fermi level in a metallic band. In that case a $\mathbf{B}_{eff}$ component due to the electronic atomic moment may arise via the contact hyperfine term. In a metal the corresponding frequency shift $K$ of the Larmor frequency is called Knight shift (Knight 1967) and is directly proportional to the Pauli spin susceptibility $\chi _{P0}$ of the metallic band. In usual metallic systems such as alkali metals the early day studies by NMR permitted to demonstrate that this technique gives indeed the best evaluation of the electronic spin susceptibility. Assuming that the main hyperfine coupling is the direct on-site contact interaction of Equ(2), which can be written

$$H_{c} =A_{0}\mathbf{I.s}\delta(\mathbf{r})
\label{eq:Hc}$$

this yields an NMR Knight shift

$$\label{eq:KSki}
K=A_{0}\chi _{P0}/\ (g\mu _{B}\hbar \gamma _{\mathrm{n}})$$

In such simple metals the spin susceptibility $\chi _{P0}=(g\mu _{B})^{2}\rho (E_{F})/2$ measures the actual density of states $\rho (E_{F})$ taken per spin direction at the Fermi level, which is typically $T$ independent as the conduction electron bandwidth is usually quite large, and the Fermi level much higher than $k_{B}T$. The Knight shift is usually a large quantity which is measured in %. This comes about because the contact coupling $A_{0}$ is usually much larger than the corresponding dipole or orbital couplings, which permits to sense very effectively the Pauli susceptibility.

## NMR in Magnetic Materials

In magnetic materials, the electronic moments are static at low temperatures, as compared with their behavior at the Curie or Néel temperatures. It follows that the static effective fields are nonzero even in the absence of any applied field. For atomic nuclei carrying an electronic moment, this field will be very large (several Tesla in general), and will give rise to a resonance at the Larmor frequency $h\nu _{L}=\hbar \gamma _{\mathrm{n}}B_{eff}$, which can be detected in the absence of an applied external field. One speaks then of Zero Field Nuclear Magnetic Resonance (ZFNMR). The fields induced on the nuclear spins of non magnetic atomic sites are generally weaker but can still give valuable information on the properties of the magnetic state as will be seen later on in various magnetic materials.

## 1D metals

The incidence of lattice dimensionality can give rather original effects, the one dimensional case (1D) being really an important one especially near half filling, that is with one electron per atom in atomic chains. One does usually expect a dimerisation of the chains at low $T$ which induces the opening of a gap at the Fermi level at low $T$ and a corresponding so called Peierls transition from a metal to an insulator (see "Peierls transition" in Wikipedia). This effect has been observed in many materials with linear metallic chains. One of the first evidence for such a Peierls transition has been found from Pt NMR data taken (Niedoba *et. al.* 1973) in the Low-Temperature Insulating State in the One-Dimensional conductor K$_{2}$Pt(CN)$_{4}$Br$_{0.3}$, 3H$_{2}$O. A similar effect is seen in the case of the organic salt TTF-TCNQ, although that case is slightly more complicated as two successive ordering transitions occur for the two type of chains involved in that structure (Jerome and Schulz 2002). Let us point out that such a metal insulator transition often occurs solely because the elastic energy loss produced by the Peierls dimerisation is encompassed by the energy gain associated with the opening of the electronic gap in the band structure. This can be unrelated to the existence of electronic correlations. More important phenomena are expected to happen in metallic 1D systems in presence of strong electronic correlations.

# Spin lattice relaxation

The local susceptibility measurements are giving pertinent information on the electronic properties of the material in its ground state. But NMR also permits to probe the excited electronic states through the fluctuations of the local field $\mathbf{B}_{eff}$. This occurs through the nuclear spin lattice relaxation (NSLR) processes which drive back the nuclear spins towards their thermodynamic equilibrium once the latter has been disturbed intentionally.

Indeed the nuclear spin magnetization is not established immediately if an external magnetic field is applied instantaneously to the material. The very interactions between the nuclear spins and the electronic degrees of freedom govern the spin lattice relaxation time $T_{1}$ which is required to establish thermodynamic equilibrium. One typically needs transverse local field components at the Larmor frequency $h\nu _{L}$ to induce the difference of population of the nuclear spin levels. Therefore $T_{1}$ is directly linked with the transverse field Fourier component at $\nu _{L}$ of $\mathbf{B}_{eff}(t)$ . One can see that, for the hyperfine couplings considered above, this resumes in a measurement of the electronic dynamic susceptibility of the electron spin system.

So, in systems with unpaired spins the dominant $T_{1}$ process is due to local field fluctuations induced by the dynamics of the local electronic magnetization. Theoretically, the spin contributions to $(T_{1}T)^{-1}$ may be written using the imaginary part of the dynamical electron spin-susceptibility ${\chi ^{\prime \prime }(\mathbf{q},\nu _{n})}$ as

$$(T_{1}T)^{-1}=\frac{2\gamma _{n}^{2}k_{B}}{g^{2}\mu _{B}^{2}}\sum_{\mathbf{q}%
}|A_0|^{2}\frac{\chi ^{\prime \prime }(\mathbf{q},\nu _{n})}{\nu
_{n}}.  \label{eq:1/T1T}$$

Here, as for the Knight shift, we assumed that the dominant hyperfine coupling is the contact term. We shall see later many examples which permit to evidence that $T_{1}$ measurements permit one to monitor the occurrence of phase transitions, to give relevant informations on energy gaps between the ground state and excited states in many electronic systems. In cases where some ionic species are mobile in a material, as for instance in ionic conductors, the atomic diffusion processes can govern the local field fluctuations sensed on some nuclear spin sites, and the $T_{1}$ measurements may permit to monitor these ionic diffusion motions.

## Spin lattice relaxation in standard 3D metals

For a simple metallic band, the dynamic electronic susceptibility response is simple enough and one writes

$$[\label{eq:kiomega}{\chi }_{0}(\omega )=\sum\nolimits_{\mathbf{q}}{\chi }_{0}(\mathbf{q}),\omega )=\chi _{P0}\left[ 1+i\pi \hbar {\omega \ }\rho (E_{F})\right] .$$

One can immediately see that this yields a simple expression for the spin lattice relaxation from Equ (7)

$$\label{eq:T1T}(T_{1}T)^{-1}=\pi k_{B}\ A_{0}^{2}\ \rho ^{2}(E_{F})/\hbar$$

so that a universal relation holds between the Knight shift and $T_{1}$.

$$\label{eq:Korringa}K^{2}T_{1}T=\mathcal{S=}\left( \hbar /4\pi k_{B}\right) \left( g\mu
_{B}/\hbar\gamma _{\mathrm{n}}\right) ^{2}.$$

As $K$ is $T$ independent, this so called "Korringa" relation applies rather well in the absence of electronic correlations. As an example one can see in Fig.  that $T_{1}T$ of $^{27}$Al is constant in pure aluminium metal on a $T$ range which extends over more than three orders of magnitude. The $T_{1}$ value in a metal is quite often used to define an empirical temperature scale especially in the very low $T$ regime below 1 K.

## Incidence of weak electron correlations

So far we did not consider any influence of electronic correlations though even in simple alkali metals electron-electron interactions play a role in the electronic scattering processes. We also do know that some electronic systems are on the verge to become magnetic. Those quasi AF or quasi ferromagnetic metals can be identified by the very fact that the Korringa relation does not apply then straightforwardly as the spin susceptibility does not behave as described above for free electron Fermi liquid systems. In such cases the dynamic spin susceptibility is not uniform in $\mathbf{q}$ space as was shown initially by (Moriya 1956), and exhibits enhanced values either for $\mathbf{q}=0$ for nearly ferromagnetic metals or for an AF wave vector $%
\mathbf{q}_{\mathbf{AF}}\ $for nearly AF materials.

In the former case the $\mathbf{q}=0$ spin susceptibility is enhanced by a factor ${S}=1/(1-{I}\chi_{P0})$ usually called the Stoner factor, and $\chi_{P}=S \chi_{P0}$ exhibits a large increase with decreasing $T$ and only saturates at very low $T$, as has been illustrated in the nearly ferromagnetic metals like TiBe$_{2}$ or in elemental Pd metal. In those cases, the enhancement of the dynamic spin susceptibility $\chi_{0}(\mathbf{q}\omega)$ is not uniform in $\mathbf{q\ }$ space and is weaker for $\mathbf{q}\neq 0$, therefore the Knight shift is more enhanced than $(T_{1}T)^{-1}$. The Korringa relation does only apply when $\chi _{P}\ $saturates (Alloul and Mihaly 1982), with a $T=0$ Korringa constant $K^{2}T_{1}T=\Delta \ \mathcal{S}$ increased by an ${S}$ dependent factor $\Delta({S})$. On the contrary in nearly AF metals $\chi_{0}(\omega)$ is peaked for $\mathbf{q}= \mathbf{q}_\mathbf{AF}$ which means that the static spin susceptibility and the Knight shift $K$ are less unhanced than $\chi_{0}(\mathbf{q}_\mathbf{AF},\omega)$. Correspondingly $K^{2}T_{1}T=\Delta \ \mathcal{S}$ corresponds in that case to$\ \Delta <1$, that is a decreased Korringa constant, as has been seen for instance in the compound MnSi (Corti et al 2007).

#  NMR in superconductors

Obviously, the establishment of a SC state yields profound transformations of the electronic properties which will be seen in the NMR response. NMR experiments do not only evidence the occurrence of SC. They also permit to characterize the properties of the SC electronic state (Mac Laughlin 1976).

## Knight shift, relaxation and gap in the SC state

One of the major effects which occur for phonon mediated SC in usual metals is the pairing of electrons in a singlet state. Such a pairing suppresses totally the normal state spin susceptibility at $T=0$. This is seen quite simply as a full suppression of the spin contribution $K_{s}$ to the Knight shift in NMR. In type I superconductors, the magnetic induction vanishes in the Meissner state which by itself forbids observation of the NMR signal in the SC state. But in type II superconductors the field penetrates as an array of vortices which becomes so dense near the upper critical field $H_{c2}$ that it becomes possible to detect the NMR signal in that regime, and to see the suppression of the Knight shift. Taking into account the variation of the SC gap and the thermal population at temperatures near $T_{c}$, yields a specific $T$ dependence of the spin susceptibility, that is of $K_{s}(T)$, which has been computed by (Yosida 1956), and which is given by

$$K_{s}(T)/K_{n}=\int_\Delta^\infty[N(0)|E|/(E^2-\Delta^2)^{1/2}](df/dE)dE
\label{eq:Yosida}$$ Here $f$ is the Fermi function. The actual variation of the knight shift with decreasing temperature can be measured and displays an agreement with this Yosida function as can be seen in Fig. .

As for spin lattice relaxation data, it can be taken in type I SC using ingenious tricks such as experiments in which the external field is cycled from a field exceeding the critical field $H_{c}$ down to a field $H<H_{c}$ in which the nuclear spin magnetization is let free to evolve under the influence of the electronic system.

The opening of the SC gap yields an activated exponential increase of the NMR spin lattice relaxation rate as $T$ approaches $0$ , which permits a determination of the SC gap magnitude from the corresponding low $T$ variation of $T_{1}^{-1}$ (Fig. ). However the great advance of BCS theory has been its ability to describe the excited states in the SC state up to $T_{c}$. Indeed in such a BCS SC state subtle effects are revealed by $T_{1}T$ data taken near $T_{c}$. An increase of the spin lattice relaxation rate above the normal state Korringa value takes place below $T_{c}$. This so called coherence peak (Hebel and Slichter, 1957) results partly from the thermal population of the increased density of electronic states which piles up above the SC gap. The $T_{1}$ only lengthens at somewhat lower temperatures than $T_{c}$ (see Fig. ).

Both these discoveries of the decrease of the spin susceptibility and of the occurrence of a Hebel-Slichter coherence peak have given the early evidences for the applicability of BCS theory of superconductivity in usual metallic systems.

## Field distribution in the mixed state of type II superconductors

In type II superconductors the magnetic induction varies significantly in space in the mixed state. This leads to a distribution of Larmor frequencies for the nuclear spins in the material. The shape of the NMR spectrum reconstitutes the histogram of the magnetic fields. Close to $H_{c2}$, singularities appear in the spectrum for values of the magnetic field corresponding to the extrema of the field distribution. The shape and width of the observed resonance can be used to deduce $\lambda$, the magnetic field penetration depth. However, for experimental reasons, NMR is not the best method for studying the superconducting state. A related technique uses elementary particles called muons. These behave like heavy electrons (or light protons), and have the property of decaying by emission of positrons in the direction of their spin. A muon whose spin is initially polarised perpendicularly to the field $B_{0}$ is implanted in the sample at time zero. One then observes the direction of the emitted positrons when it decays. By repeating this experiment for a large number of events, the free precession signal of the muon spin can be reconstructed statistically. This experiment is equivalent to an NMR experiment, and can be used to determine $\lambda$. Since muons can be implanted in almost any sample, it has been possible to make comparative measurements of $\lambda$ in a wide range of superconducting materials.

#  Quadrupolar electrostatic interactions and NQR

Generally in solid state NMR an atomic nucleus with spin $\mathbf{I}$ has its spin energy levels determined by the Zeeman interaction with the external magnetic field $\mathbf{B}_{0}$. But nuclei with $I>1/2$ , in which the nuclear charge has a non spherical shape are also characterized by the quadrupole moment $Q$ of charges . This charge distribution interacts with the electrostatic electric fields (in fact the electric field gradient EFG) induced on the nucleus by the distribution of electronic charges in the material. Therefore the nuclear spins are influenced both by the Zeeman and the quadrupolar Hamiltonian $\mathcal{H}_{Q}$, which can be written (Abragam 1961, Slichter 1963)

$$\label{eq:HQ}
\mathcal{H}_Q = \frac{eQ}{2l(I-1)}\sum_{\alpha }V_{\alpha \alpha }I_{\alpha}^2$$

where $V_{\alpha \alpha }$ are the second derivatives of the electric potential, that is the EFG tensor which has $\alpha =X,Y,Z$ as principal axes. The three quadrupolar frequencies

$$\label{eq:nuQ}
\nu _{\alpha }=\frac{3eQV_{\alpha \alpha }}{ 2I(2I-1)h}$$

are linked by Laplace equation $\sum_{\alpha }V_{\alpha \alpha}=\sum_{\alpha }\nu _{\alpha }=0$. Therefore it is more common to use two parameters - the quadrupolar frequency $\nu _{Q}=\nu _{Z}$ corresponding to the largest principal axis component $V_{ZZ}$ of the EFG tensor and the asymmetry parameter $\eta =(V_{XX}-V_{YY})/V_{ZZ}$ (here the principal axes of the EFG tensor are chosen following $\left\vert V_{ZZ}\right\vert \geq \left\vert V_{YY}\right\vert \geq \left\vert V_{XX}\right\vert$).

##  Quadrupole split NMR spectrum

The quadrupolar term modifies the NMR spectrum and can be treated as a perturbation if the quadrupole term is small as compared to the Zeeman splitting. Let us consider as an example the $^{23}$Na NMR in the layered cobaltates Na$_{x}$CoO$_{2}$. The $^{23}$Na nuclear spin is $I=3/2$ and therefore for a given direction of applied magnetic field $H_{0}$ relative to the crystal the NMR spectrum for a single Na site consists of 3 lines - a central line which corresponds to the $-1/2\leftrightarrow 1/2$ transition and two satellites corresponding to the $-3/2\leftrightarrow -1/2$ and $1/2\leftrightarrow 3/2\ $transitions. The position of the central line is determined by the applied field and the distance between the two satellite lines depends on the orientation of the external field with respect to the principal axis of EFG tensor described by the spherical angular coordinates $\theta$ and $\varphi$ which can be expressed as

$$\Delta \nu =\nu _{Q}(3\cos ^{2}\theta -1+\eta \sin ^{2}\theta \cos 2\varphi).
\label{eq:SatelDist}$$

In Fig. (a), we show such a spectrum obtained with the field $B$ perpendicular to the CoO$_{2}$ layers in Na$_{x}$CoO$_{2}$ for $x=2/3$. One can see that the spectrum is somewhat complicated and can be fitted by three spectra with different $\nu _{Q}$ that is three different Na sites. Such a spectrum permits to establish that the Na sites are in an ordered pattern with three sites per unit cell for this Na composition (Alloul et al 2009).

## Pure quadrupole resonance : NQR

One can immediately see that in the absence of an applied external field the quadrupole Hamiltonian lifts partially the degeneracy of the spin states. For instance for a spin $I=3/2$ one gets two levels which are separated in frequency by (Abragam 1961)

$$\nu =\nu _{Q}\sqrt{1+\eta ^{2}/3} \label{eq:nuNQR}$$

so that the quadrupole frequency can be determined directly in a pure nuclear quadrupole resonance experiment taken in zero magnetic field. It is for instance illustrative to see that the $^{23}$Na NQR of Na$_{2/3}$CoO$_{2}$ shown in Fig. (b) indeed displays the three frequencies expected from the analysis of the NMR spectrum, which perfectly confirms the existence of the three Na lattice sites in the unit cell (Platova et al 2009).

This example permits us to evidence that the nuclear spin resonance spectra of nuclei with $I>1/2$ contain information relevant to the atomic structure or electronic charge distributions in materials. Thus the quadrupole parameters may be regarded as complementary tools with respect to diffraction techniques for investigating structural changes or charge ordering in solids. Let us point out that the NQR, being a zero magnetic field experiment, is quite useful in the studies of electronic properties in the SC state. This permits for instance to restrict the technical limitations in measurements of the spin lattice relaxation in type I superconductors.

# Impurities and disorder

All real crystalline materials contain structural defects. Those are often impurity atoms substituted to some atoms in the ideal structure, or disorder induced by vacancies on some sites of the atomic structure or by deviations to the ideal structural arrangement of atoms. The incidence of specific substituted impurities on the physical properties of the material is sometimes well understood. The essence of the observed phenomena is that an impurity is a local screened Coulomb potential, which ideally is a uniform perturbation in **q** space, inducing a response which is inhomogeneous in real space but which reflects the response to all **q** values. In a classical metallic system, since the response is homogeneous up to $\left\vert \mathbf{q}\right\vert =k_{F}$, the main detected feature comes from this truncation of the response at $k_{F}$, which yields the well known Friedel oscillations in the local density of states for a charge defect, and the RKKY oscillation for a spin defect. Such effects differentiate the electronic response on sites in the vicinity of the defect and can be detected by NMR experiments. We shall first illustrate this in the case of magnetic defects in sp metals and show that NMR spectra permit to probe directly the RKKY oscillations.

##  Local RKKY spin density oscillations induced by magnetic impurities in metals

If one substitutes a magnetic impurity on a lattice site of a metal such as Cu, the neighbouring sites are differentiated and the magnetic response depends on the distance to the impurity. For local moment impurities, Yosida (1956) calculated explicitly the spin density oscillations assuming that the free electron spin $s$ and local moment $S$ interact by an exchange interaction

$$H=-J\ S.\ s\ \delta (r).
\label{eq:Exchange}$$

The resulting local spin density, calculated in perturbation theory at a position $R_{n}$ with respect to the impurity is given by

$$n(R_{n})=-\frac{1}{4\pi }J\rho (E_{F})\ \frac{cos(2k_{F}R_{n})}{R_{n}^{3}}%
<S_{z}>, \label{eq:RKKY}$$

for a field applied in the $z$ direction. For local moment impurities such as Mn in Cu the NMR Knight shift of a Cu nuclear spin at position $%
R_{n}$ with respect to the impurity acquires an extra shift $\Delta K$ given in an applied field $H$ by

$$H\Delta K(R_n)=A_{hf}\ n(R_{n})=A(R_{n})<S_{z}>,
\label{eq:Shift(R)}$$ where $A_{hf}$ is the on site Cu hyperfine coupling. In most dilute alloys of transition elements only a few near-neighbor shells of the impurity could be resolved. However in Cu-Mn, the impurity magnetization $<S_{z}>$ becomes so large at low $T$ and high fields, that $n(R_{n})$ becomes sizable on many neighboring sites of the impurity. Up to 17 distinct shells of neighbors could be detected in that case, as can be seen in the spectra of Fig. , which gives a straightforward illustration of the occurrence of spin density oscillations. One can see that there are about as many extra lines (we call them satellite lines) on the right and on the left of the central line. The technical details about the assignment of the different lines to specific shells of neighbors, and the analysis of the spatial dependence of the spin polarization which can be deduced from those data are summarized in (Alloul 2012). One could push the analysis at a stage permitting to confirm the overall $R_{n}^{-3}$ dependence, but also to evidence deviations from the asymptotic limit at short distances. Such deviations with respect to the spatial dependence of (17) could be analyzed by using a more reliable model than an exchange coupling between the impurity and electron spins.

## Transferred hyperfine couplings

So far we have only considered hyperfine couplings between nuclear spins and the electrons involved in the atomic orbitals of the corresponding site. However the illustration of the RKKY interaction given in the section above did allow us to demonstrate that electrons on a given atomic site interact as well with the neighboring sites. Indeed if we consider Equ.(18) we see that it permits to define a transferred hyperfine coupling $A(R_{n})$ between the local moment at the origin and the nuclear spin at $R_{n}$.

Therefore in systems where different atomic sites are involved in a unit cell such transferred hyperfine couplings do play an important role. In particular, if one considers atomic sites which are not displaying the most important magnetic response they will still sense the response of the magnetic sites through such transferred hyperfine couplings. In the case of Cu-Mn the transferred hyperfine coupling with the Mn magnetism extends on a large number of Cu sites. But, as the transferred hyperfine coupling decreases strongly with distance, it is often sufficient to consider solely transferred hyperfine couplings with the first nearest neighbors. For instance in a square lattice such as that of the CuO$_{2}$ planes of cuprate superconductors for which the electronic magnetic response is located on the Cu sites, the $^{63}$Cu nuclei will be coupled to the on site magnetic response and with that of the near neighbors. This can be cast in a combined wave-vector$\ \mathbf{q}$ dependence of the hyperfine coupling $A(\mathbf{q})=A_{0}+\sum A_{i}e^{i%
\mathbf{q}\mathbf{r}_{i}}$ in which $A_{0}$ is the on site local hyperfine interaction between the observed nuclear and electron spins and $A_{i}$ is the hyperfine interaction with electron spins at neighboring sites at $\mathbf{r}_{i}$. Similarly one has to substitute $A_{0}$ by $A(\mathbf{q})$ in the Equ (7) for the spin lattice relaxation (Moriya 1956).

This modification is important as it shows that the NMR experiments permit to sense the $\ \mathbf{q}$ dependence of the dynamic susceptibility in complex materials involving diverse atomic sites or a multiband electronic structure. Combining NMR shift and macroscopic susceptibility data gives access to an experimental determination of $A_{i}$ which contains information on the electronic structure.

##  Disorder as viewed from quadrupolar effects

We have seen in sec. that the nuclear spin resonance spectra of nuclei with $I>1/2$ contain information relevant to the atomic structure or electronic charge distributions in materials. Thus the quadrupole parameters may be regarded as complementary tools with respect to diffraction techniques for investigating structural changes or charge ordering in solids. The NMR and NQR techniques are also important for instance to characterize the local disorder in solids whereas the interpretation of the data on disordered materials from usual scattering experiments, such as X-ray or neutron scattering, is complicated by the absence of the translation symmetry of long-range order. One will have the opportunity to discuss these aspects in relation with the physical study of correlated electron systems, which, being complicated materials, do involve many such defects.

# References

1.  A. Abragam, "The Principles of nuclear magnetism" (Oxford: Clarendon Press, London, 1961).

2.  Alloul H and Mihaly L . "Ti NMR study of the nearly ferromagnetic system TiBe2". Phys.Rev.Letters, 48, 1420 (1982).

3.  H. Alloul, I. R. Mukhamedshin, T. A. Platova and A. V. Dooglav , "Na ordering imprints a metallic kagomé lattice onto the Co planes of Na2/3CoO2", Europhysics Letters 85 , 47006 (2009).

4.  H. Alloul, "From Friedel oscillations and Kondo effect to the pseudogap in cuprates", in J Supercond. Nov. Mag. 25, 385 (2012) ; (DOI) 10.1007/s10948-012-1472-x, arXiv:1204.3804.

5.  M. Corti, F. Carbone, M. Filibian, Th. Jarlborg, A. A. Nugroho, and P. Carretta , "Spin dynamics in a weakly itinerant magnet from Si29 NMR in MnSi", Phys. Rev. B 75, 115111 (2007)

6.  Hebel L.C. and Slichter C.P, "Nuclear Spin Relaxation in Normal and Superconducting Aluminum", Phys Rev 113, 1504 (1957)

7.  Jerome D. and Schulz H. J.,"Organic conductors and superconductors", Adv. Phys., 51, 293 (2002).DOI: 10.1080/00018730110116362

8.  Knight W.D. "Manetic Resonance and relaxation", Ed R. Blinc, North Holland 311 (1967)

9.  D.E. Mac Laughlin "Magnetic Resonance in the Superconducting state" in Solid State Physics, Academic Press , New York 31 ,1 (1976)

10. Moriya T., "Nuclear Magnetic Relaxation in Antiferromagnetics" Prog. Theor. Phys., 16, 23 (1956).

11. H. Niedoba, H. Launois, D. Brinkmann, R. Brugger and H. R. Zeller, "NMR Study of the Low-Temperature Insulating State in the One-Dimensional "Conductor" $K_2Pt (CN) 4Br0. 3.3H_2O$", Phys. Status Solidi, 58b 309(1973) DOI: 10.1002/pssb.2220580130

12. T. A. Platova, I. R. Mukhamedshin, H. Alloul, A. V. Dooglav and G. Collin, "NQR and X-ray investigation of the structure of Na2/3CoO2 compound",cond-mat/0908.0834, Phys. Rev. B 80, 224106 (2009). DOI:10.1103/PhysRevB.80.224106

13. C. P. Slichter, "Principles of Magnetic Resonance", Harper and Row (1963)(Springer-Verlag, New York, (1989), 3rd ed.

14. P. Wzietek, T. Mito, H. Alloul, D. Pontiroli, M. Aramini and M. Riccò, "NMR study of the Superconducting gap variation near the Mott transition in Cs3C60" arXiv:1310.5529, Phys. Rev. Lett. 112, 066401 (2014).

15. Yosida K , "Paramagnetic Susceptibility in Superconductors" Phys. Rev. 110 , 769 (1956)

# Further reading

1.  Abragam A. Bleaney B. "Electronic Paramagnetic Resonance of Transition ions" (Dover publications)(1986)

2.  Alloul H. "Introduction to the physics of Electrons in Solids ", Graduate texts in Physics, Springer (Heidelberg) (2011), ISBN 978-3-642-13564-4 DOI:10.1007/978-3-642-13565-1

3.  Ashcroft Neil W., Mermin N. David, Saunders College (1976), ISBN 0030493463, 9780030493461

4.  Kittel C., "Introduction to Solid State Physics", 8th Edition, Wiley (2005), ISBN 047141526X, 9780471415268

5.  Tinkham M. , "Introduction to Superconductivity", Gordon and Breach (New York 1965), Dover Publications (1996), ISBN 0486134725, 9780486134727

# See also

1.  "Bardeen-Cooper-Schrieffer theory", Leon Cooper and Dimitri Feldman (2009), Scholarpedia, 4(1):6439.

2.  "Magnetic resonance imaging", Joan Dawson and Paul C. Lauterbur Scholarpedia, 3(7):3381(2008)).

3.  "NMR in strongly correlated materials", H. Alloul , Scholarpedia, 10(1):30632 (2015).

4.  <a href="http://en.wikipedia.org/wiki/Nuclear_magnetic_resonance Nuclear magnetic resonance" class="uri">http://en.wikipedia.org/wiki/Nuclear_magnetic_resonance Nuclear magnetic resonance</a> in Wikipedia

5.  <a href="http://en.wikipedia.org/wiki/Peierls_transition Peierls transition" class="uri">http://en.wikipedia.org/wiki/Peierls_transition Peierls transition</a> in Wikipedia

[^1]: Henri Alloul (2014), *NMR studies of electronic properties of solids*, Scholarpedia, 9(9):32069.


author: J. D. Bailey; J. D. Landstreet; S. Bagnulo
bibliography: abund-evol.bib
date: Accepted 26 November 2013
title: Discovery of secular variations in the atmospheric abundances of magnetic Ap stars[^1]

We want to establish whether abundance peculiarities change as stars evolve on the main sequence, and provide observational constraints to diffusion theory. 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$_{\odot}$. For each star, we measured the abundances of He, O, Mg, Si, Ti, Cr, Fe, Pr and Nd. 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.

# 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_{\rm eff} \le 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_{\rm eff} \ge 20\,000$ K), the intense outward 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 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 .

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 . 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_{\rm X}/N_{\rm 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 upward and downward diffusion . 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 to more detailed maps of both the magnetic field structure and abundance distributions using spectroscopic observations in all four Stokes parameters . 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^8 - 10^9$ 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 .

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.

# Observations

| Star        |   Instrument |      R | $\lambda$ (Å) |
|:------------|-------------:|-------:|--------------:|
| HD 45583    |    FEROS (2) |  48000 |     3528-9217 |
| HD 61045    | ESPaDOnS (2) |  65000 |    3690-10481 |
| HD 63401    |    FEROS (2) |  48000 |     3528-9217 |
| HD 74535    |    FEROS (2) |  48000 |     3528-9217 |
| HD 133652   | ESPaDOnS (1) |  65000 |    3690-10481 |
|             |    FEROS (1) |  48000 |     3528-9217 |
| HD 133880   | ESPaDOnS (2) |  65000 |    3690-10481 |
| HD 147010   | ESPaDOnS (2) |  65000 |    3690-10481 |
| HD 162576   | ESPaDOnS (2) |  65000 |    3690-10481 |
| HD 162725   | ESPaDOnS (2) |  65000 |    3690-10481 |
| HD 304842   |    FEROS (2) |  48000 |     3528-9217 |
| NGC 2169 12 |     UVES (1) | 110000 |    3070-10398 |
| BD+00 1659  | ESPaDOnS (1) |  65000 |    3690-10481 |
| BD-19 5044L | ESPaDOnS (2) |  65000 |    3690-10481 |
| BD+49 3789  | ESPaDOnS (2) |  65000 |    3690-10481 |
| HIP 109911  | ESPaDOnS (2) |  65000 |    3690-10481 |
|             |              |        |               |

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.

Our study used high dispersion spectra of 15 stars, listed in Table . 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$_{\odot}$) 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_\odot$.

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_{\odot})$ versus $T_{\rm 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 .

For each star, accurate values of effective temperature ($T_{\rm eff}$), gravity ($\log g$), evolutionary mass, magnetic field strength and age were adopted from the literature . When values of $T_{\rm eff}$ and $\log g$ were not available, we derived them from available Geneva and $uvby\beta$ photometry. For the Geneva photometry the fortran program described by was used. For the Strömgren $uvby\beta$ photometry we used a version of the fortran program "UVBYBETANEW" that corrects the $T_{\mathrm{eff}}$ of the magnetic Bp stars to the appropriate temperature scale . As per the discussions of and , the uncertainties in $T_{\mathrm{eff}}$ and $\log g$ were taken to be about $\pm$`<!-- -->`{=html}500 K and 0.2 dex, respectively. Table  lists the properties of each star, including its designation, associated cluster, age, mass, $T_{\rm eff}$, $\log g$, and $v \sin i$. The final column shows the root-mean-square magnetic field strength ($B_{\rm rms}$) as computed from whatever modern measurements of mean longitudinal magnetic field $B_{\rm z}$ are available ; 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_{\rm 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$_{\odot}$ from ZAMS to TAMS.

# Modelling technique

```latex
\begin{table*}
\caption{Physical properties of the stars studied.}
\begin{tabular}{lrlrrrrr}
\hline\hline
Star & Cluster & $\log t$ & M/M$_{\odot}$ & \ensuremath{T_{\mathrm{eff}}}\ (K) & $\log g$\ & \ensuremath{v \sin i}\ (km s$^{-1}$) & B$_{\rm rms}$ (G)\\
\hline
HD~147010 & Upper Sco & 6.70 $\pm$ 0.10$^{1}$ & 3.15 $\pm$ 0.20$^{1}$ & 13000 $\pm$ 500$^{2}$ & 4.40 $\pm$ 0.20$^{2}$ & 15 $\pm$ 2$^{2}$ & 4825$^{1}$ \\
NGC~2169~12 & NGC~2169 & 6.97 $\pm$ 0.10$^{1}$ & 3.65 $\pm$ 0.15$^{1}$ & 13800 $\pm$ 500$^{1}$ & 4.30 $\pm$ 0.20 & 56 $\pm$ 5 & 3410$^{1}$\\
HD~133652 & Upper Cen Lup & 7.20 $\pm$ 0.10$^{1}$ & 3.35 $\pm$ 0.15$^{1}$ & 13000 $\pm$ 500 & 4.30 $\pm$ 0.20 & 48 $\pm$ 2 & 1120$^{1}$\\
HD~133880 & & & 3.20 $\pm$ 0.15$^{1}$ & 13000 $\pm$ 600$^{3}$ & 4.34 $\pm$ 0.16$^{3}$ & 103 $\pm$ 10$^{3}$ & 2300$^{1}$ \\
HD~45583 & NGC~2232 & 7.55 $\pm$ 0.10$^{1}$ & 3.30 $\pm$ 0.15$^{1}$ & 12700 $\pm$ 500$^{2}$ & 4.20 $\pm$ 0.20$^{2}$ & 70 $\pm$ 6$^{2}$ & 2730$^{1}$\\
HD~63401 & NGC~2451 & 7.70 $\pm$ 0.10$^{1}$ & 3.70 $\pm$ 0.20$^{1}$ & 13500 $\pm$ 500$^{1}$ & 4.20 $\pm$ 0.20 & 52 $\pm$ 4 & 365$^{1}$ \\
HD~74535 & IC~2391 & 7.70 $\pm$ 0.15$^{1}$ & 3.85 $\pm$ 0.15$^{1}$ & 13600 $\pm$ 500$^{1}$ & 4.30 $\pm$ 0.20 & 45 $\pm$ 4 & 95$^{1}$\\
BD-19~5044L & IC~4725 & 8.02 $\pm$ 0.08$^{1}$ & 3.55 $\pm$ 0.15$^{1}$ & 12800 $\pm$ 500$^{2}$ & 4.50 $\pm$ 0.20$^{2}$ & 15 $\pm$ 3$^{2}$ & 235$^{1}$\\
BD+49~3789 & NGC~7243 & 8.06 $\pm$ 0.10$^{*}$ & 3.55 $\pm$ 0.15$^{*}$ & 12900 $\pm$ 500$^{2}$ & 4.20 $\pm$ 0.20$^{2}$ & 85 $\pm$ 5$^{2}$ & 561$^{*}$\\
HIP~109911 & & & 3.65 $\pm$ 0.15$^{*}$ & 13000 $\pm$ 500 & 4.30 $\pm$ 0.20 & 60 $\pm$ 2 & 348$^{*}$\\
HD~61045 & NGC~2422 & 8.08 $\pm$ 0.11$^{1}$ & 3.85 $\pm$ 0.20$^{1}$ & 13000 $\pm$ 500$^{2}$ & 4.10 $\pm$ 0.20$^{2}$ & 64 $\pm$ 3$^{2}$ & 430$^{1}$\\
HD~304842 & NGC~3114 & 8.13 $\pm$ 0.15$^{1}$ & 3.55 $\pm$ 0.15$^{1}$ & 12500 $\pm$ 500$^{2}$ & 3.90 $\pm$ 0.20$^{2}$ & 65 $\pm$ 5$^{2}$ & 20$^{1}$\\
BD+00~1659 & NGC~2301 & 8.22 $\pm$ 0.10$^{*}$ & 3.65 $\pm$ 0.15$^{*}$ & 12500 $\pm$ 500$^{2}$ & 4.00 $\pm$ 0.20$^{2}$ & 7.0 $\pm$ 1$^{2}$ & 394$^{*}$\\
HD~162576 & NGC~6475 & 8.41 $\pm$ 0.13$^{1}$ & 3.10 $\pm$ 0.15$^{*}$ & 10300 $\pm$ 500 & 3.70 $\pm$ 0.20 & 28 $\pm$ 3 & 15$^{*}$\\
HD~162725 & & & 3.30 $\pm$ 0.20$^{1}$ & 10000 $\pm$ 500 & 3.50 $\pm$ 0.20 & 31 $\pm$ 3 & 69$^{*}$\\
\hline
\multicolumn{8}{p{0.85\textwidth}}{{\sc references} -- (1) \citet{paper2}; (2) \citet{BL2013}; (3) \citet{Bailey2012}; (*) This work} \\
\label{properties}
\end{tabular}
\end{table*}
```

To determine the atmospheric abundances of the magnetic Bp stars, the fortran program zeeman was used . 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 $\beta$ 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_{\rm eff}$ and $\log g$ assumed. The atomic data for individual spectral lines are taken from the Vienna Atomic Line Database . 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 ; 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 $\chi^{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 $v_{\rm R}$ and $v \sin i$. The abundance ($\log N_{\rm X}/N_{\rm 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_{\rm 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  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.

| Element | $\lambda$ (Å) | Element | $\lambda$ (Å) |
|:--------|--------------:|:--------|--------------:|
| He i    |      4437.551 | Si ii   |      5055.984 |
|         |      4713.139 |         |      5056.317 |
|         |      5015.678 | Ti ii   |      4533.960 |
|         |      5047.738 |         |      4563.757 |
|         |      5875.599 |         |      4571.968 |
|         |      5875.614 | Cr ii   |      4558.650 |
|         |      5875.615 |         |      4565.739 |
|         |      5875.625 |         |      4588.199 |
|         |      5875.640 |         |      4592.049 |
|         |      5875.966 | Fe ii   |      4541.524 |
| O i     |      6155.961 |         |      4555.893 |
|         |      6155.966 |         |      4556.392 |
|         |      6155.986 |         |      4583.837 |
|         |      6156.736 |         |      4583.999 |
|         |      6156.755 |         |      5029.097 |
|         |      6156.776 |         |      5030.630 |
|         |      6158.146 |         |      5032.712 |
|         |      6158.172 |         |      5035.708 |
|         |      7771.941 | Pr iii  |      6160.233 |
|         |      7774.161 |         |      6161.194 |
|         |      7775.388 |         |      7781.983 |
| Mg ii   |      4481.126 | Nd iii  |      4911.653 |
|         |      4481.150 |         |      4912.944 |
|         |      4481.325 |         |      4914.094 |
| Si ii   |      4621.418 |         |      5050.695 |
|         |      4621.722 |         |      6145.068 |
|         |      5041.024 |         |               |
|         |               |         |               |

List of spectral lines modelled.

All of our stars have measurements of the line-of-sight magnetic field $\langle B_z \rangle$, either previously published or from our own unpublished results. For two of the stars in our sample, magnetic field models are available from detailed studies: HD 133880 and HD 147010 . 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_{\rm 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.

## 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$_{\odot}$ sample. In Table , we tabulate the mean abundances **$\log (N_{\rm X}/N_{\rm H})$** for each star. The first three columns recall the star designation, age ($\log t$) and $T_{\rm eff}$ from Table . 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 . 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  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_{\rm eff}$ range of this study . 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 suggests that this discrepancy is mostly due to strong vertical stratification of Si, similar to the stratification profiles computed by for stellar atmospheres with $T_{\rm 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 , mainly in cool magnetic Ap stars. We expect this phenomenon to occur in the $T_{\rm eff}$ range discussed in this paper as well. Because both observation and theory 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 ). 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_{\mathrm{eff}}$. However, the ambiguity due to probable vertical stratification and horizontal inhomogeneities should be kept in mind in evaluating our results.

```latex
\begin{sidewaystable*}
\vspace{20cm}
\caption{Average derived abundances for the stars studied.}
\centering
\begin{tabular}{llllllllllll}
\hline\hline
Star & $\log t$ & \ensuremath{T_{\mathrm{eff}}}\ (K) & $\log$(He/H) & $\log$(O/H) & $\log$(Mg/H) & $\log$(Si~{\sc ii}/H) & $\log$(Ti/H) & $\log$(Cr/H) & $\log$(Fe/H) & $\log$(Pr/H) & $\log$(Nd/H) \\
\hline
HD~147010 & 6.70 $\pm$ 0.10$^{1}$ & 13000 $\pm$ 500$^{2}$ & $-3.11 \pm 0.28$ & $-4.49 \pm 0.28$ & $-5.67 \pm 0.14$ & $-3.72 \pm 0.14$ & $-5.68 \pm 0.21$ & $-4.04 \pm 0.14$ & $-3.33 \pm 0.14$ & $-6.15 \pm 0.28$ & $-6.50 \pm 0.28$ \\
NGC~2169~12 & 6.97 $\pm$ 0.10$^{1}$ & 13800 $\pm$ 500$^{1}$ & $-3.09 \pm 0.30$ & $-3.88 \pm 0.15$ & $-5.56 \pm 0.15$ & $-3.87 \pm 0.15$ & $-5.90 \pm 0.20$ & $-4.80 \pm 0.20$ & $-3.76 \pm 0.10$ & $-7.31 \pm 0.20$ & $-7.21 \pm 0.20$ \\
HD~133652 & 7.20 $\pm$ 0.10$^{1}$ & 13000 $\pm$ 500 & $-3.16 \pm 0.28$ & $-3.87 \pm 0.28$ & $-5.18 \pm 0.14$ & $-3.36 \pm 0.28$ & $-5.01 \pm 0.21$ & $-4.17 \pm 0.21$ & $-2.90 \pm 0.14$ & $-6.98 \pm 0.28$ & $-6.42 \pm 0.28$ \\
HD~133880 & & 13000 $\pm$ 600$^{3}$ & $\leq-2.00$ & $-2.90 \pm 0.42$ & $-4.12 \pm 0.28$ & $-2.94 \pm 0.28$ & $-5.44 \pm 0.28$ & $-4.65 \pm 0.21$ & $-3.56 \pm 0.14$ & $-6.96 \pm 0.28$ & $-7.04 \pm 0.42$ \\
HD~45583 & 7.55 $\pm$ 0.10$^{1}$ & 12700 $\pm$ 500$^{2}$ & $\leq-2.40$ & $-3.61 \pm 0.42$ & $-4.19 \pm 0.28$ & $-3.34 \pm 0.28$ & $-5.68 \pm 0.28$ & $-4.74 \pm 0.21$  & $-3.43 \pm 0.21$ & $-7.49 \pm 0.42$ & $-7.04 \pm 0.42$ \\
HD~63401 & 7.70 $\pm$ 0.10$^{1}$ & 13500 $\pm$ 500$^{1}$ & $-2.64 \pm 0.28$ & $-3.73 \pm 0.21$ & $-5.86 \pm 0.14$ & $-3.96 \pm 0.14$ & $-6.44 \pm 0.28$ & $-5.32 \pm 0.28$ & $-3.84 \pm 0.14$ & $-7.14 \pm 0.28$ & $-7.44 \pm 0.28$ \\
HD~74535 & 7.70 $\pm$ 0.15$^{1}$ & 13600 $\pm$ 500$^{1}$ & $-2.52 \pm 0.28$ & $-4.26 \pm 0.28$ & $-5.44 \pm 0.28$ & $-4.28 \pm 0.14$ & $-6.29 \pm 0.28$ & $-5.65 \pm 0.14$ & $-4.02 \pm 0.28$ & $-8.15 \pm 0.21$ & $-7.52 \pm 0.28$ \\
%7.80 $\pm$ 0.15$^{1}$ & HD~318107 & 11800 $\pm$ 500$^{1}$ & $\leq-2.50$ & $-3.65 \pm 0.20$ & $-5.70 \pm 0.20$ & $-3.65 \pm 0.15$ & $-5.05 \pm 0.10$ & $-4.50 \pm 0.15$ & $-3.00 \pm 0.20$ & $-6.35 \pm 0.20$ & $-6.20 \pm 0.15$ \\
BD-19~5044L & 8.02 $\pm$ 0.08$^{1}$ & 12800 $\pm$ 500$^{2}$ & $-2.16 \pm 0.21$ & $-3.99 \pm 0.28$ & $-5.44 \pm 0.21$ & $-3.94 \pm 0.28$ & $-6.82 \pm 0.28$ & $-5.88 \pm 0.28$ & $-4.24 \pm 0.28$ & $-8.12 \pm 0.28$ & $-7.88 \pm 0.42$ \\
BD+49~3789 & 8.06 $\pm$ 0.10$^{*}$ & 12900 $\pm$ 500$^{2}$ & $-2.52 \pm 0.35$ & $-3.60 \pm 0.22$ & $-5.25 \pm 0.25$ & $-3.54 \pm 0.14$ & $-6.33 \pm 0.35$ & $-5.44 \pm 0.28$ & $-4.06 \pm 0.28$ & $-7.72 \pm 0.28$ & $-7.72 \pm 0.25$ \\
HIP~109911 & & 13000 $\pm$ 600 & $-2.34 \pm 0.28$ & $-3.69 \pm 0.21$ & $-5.38 \pm 0.14$ & $-3.49 \pm 0.14$ & $-5.90 \pm 0.28$ & $-5.35 \pm 0.28$ & $-3.59 \pm 0.28$ & $-7.70 \pm 0.28$ & $-7.37 \pm 0.28$ \\
HD~61045 & 8.08 $\pm$ 0.11$^{1}$ & 13000 $\pm$ 500$^{2}$ & $-2.06 \pm 0.28$ & $-3.72 \pm 0.28$ & $-5.10 \pm 0.14$ & $-3.94 \pm 0.28$ & $-6.24 \pm 0.28$ & $-5.47 \pm 0.28$ & $-3.98 \pm 0.28$ & $-7.74 \pm 0.28$ & $-7.77 \pm 0.42$ \\
HD~304842 & 8.13 $\pm$ 0.15$^{1}$ & 12500 $\pm$ 500$^{2}$ & $-1.54 \pm 0.28$ & $-3.78 \pm 0.25$ & $-5.76 \pm 0.21$ & $-3.64 \pm 0.28$ & $-6.83 \pm 0.14$ & $-6.05 \pm 0.32$ & $-4.35 \pm 0.32$ & $-7.65 \pm 0.36$ & $-7.04 \pm 0.28$ \\
BD+00~1659 & 8.22 $\pm$ 0.10$^{*}$ & 12500 $\pm$ 500$^{2}$ & $-2.34 \pm 0.15$ & $-3.58 \pm 0.10$ & $-5.62 \pm 0.10$ & $-3.53 \pm 0.20$ & $-6.69 \pm 0.10$ & $-5.02 \pm 0.10$ & $-3.79 \pm 0.10$ & $-8.11 \pm 0.20$ & $-7.34 \pm 0.20$ \\
HD~162576 & 8.41 $\pm$ 0.13$^{1}$ & 10300 $\pm$ 500 & $-1.51 \pm 0.28$ & $-3.77 \pm 0.28$ & $-5.34 \pm 0.14$ & $-4.44 \pm 0.14$ & $-7.72 \pm 0.14$ & $-5.12 \pm 0.14$ & $-3.95 \pm 0.14$  & $-10.15 \pm 0.42$ & $-9.10 \pm 0.42$ \\
HD~162725 & & 10000 $\pm$ 500 & $-1.74 \pm 0.36$ & $-3.77 \pm 0.28$ & $-5.23 \pm 0.14$ & $-3.54 \pm 0.28$  & $-7.30 \pm 0.28$ & $-4.64 \pm 0.22$ & $-3.66 \pm 0.28$ & $-8.32 \pm 0.42$ & $-7.44 \pm 0.28$ \\
Sun & & & $-1.07$ & $-3.31$ & $-4.40$ & $-4.49$ & $-7.05$ & $-6.36$ & $-4.50$ & $-11.28$ & $-10.58$\\
\hline
\multicolumn{12}{p{0.95\textwidth}}{{\sc references} -- (1) \citet{paper2}; (2) \citet{BL2013}; (3) \citet{Bailey2012}; (*) This work}\\
\label{abundances}
\end{tabular}
\end{sidewaystable*}
```

# Results

| y-axis           | x-axis             |              Slope | $\sigma_{D}$ |
|:-----------------|:-------------------|-------------------:|-------------:|
| $\log(\rm He/H)$ | $\log t$           |    0.76 $\pm$ 0.18 |          4.2 |
|                  | $\log B_{\rm rms}$ | $-0.48$ $\pm$ 0.11 |          4.4 |
| $\log(\rm O/H)$  | $\log t$           |    0.13 $\pm$ 0.17 |         0.76 |
|                  | $\log B_{\rm rms}$ |    0.04 $\pm$ 0.13 |         0.31 |
| $\log(\rm Mg/H)$ | $\log t$           | $-0.15$ $\pm$ 0.26 |         0.58 |
|                  | $\log B_{\rm rms}$ |    0.20 $\pm$ 0.18 |          1.1 |
| $\log(\rm Si/H)$ | $\log t$           | $-0.13$ $\pm$ 0.20 |         0.65 |
|                  | $\log B_{\rm rms}$ |    0.19 $\pm$ 0.13 |          1.5 |
| $\log(\rm Ti/H)$ | $\log t$           | $-1.07$ $\pm$ 0.27 |          4.0 |
|                  | $\log B_{\rm rms}$ |    0.69 $\pm$ 0.12 |          5.8 |
| $\log(\rm Cr/H)$ | $\log t$           | $-0.97$ $\pm$ 0.24 |          4.0 |
|                  | $\log B_{\rm rms}$ |    0.40 $\pm$ 0.16 |          2.5 |
| $\log(\rm Fe/H)$ | $\log t$           | $-0.41$ $\pm$ 0.16 |          2.6 |
|                  | $\log B_{\rm rms}$ |    0.24 $\pm$ 0.08 |          3.0 |
| $\log(\rm Pr/H)$ | $\log t$           | $-1.29$ $\pm$ 0.29 |          4.4 |
|                  | $\log B_{\rm rms}$ |    0.88 $\pm$ 0.20 |          4.4 |
| $\log(\rm Nd/H)$ | $\log t$           | $-0.81$ $\pm$ 0.24 |          3.4 |
|                  | $\log B_{\rm rms}$ |    0.46 $\pm$ 0.17 |          2.7 |
|                  |                    |                    |              |

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, $\sigma_{D}$

## Abundance variations with time

Figure  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 , and in abundance from Table . 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_\odot$ 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_{\rm X}/ N_{\rm H}$ versus $\log t$ (dashed blue line). Table  tabulates the slope of the best-fit line with its uncertainty, as well as the significance of the slope, $\sigma_{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_\odot$.

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.

### 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_{\mathrm{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 \sim 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.

### 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 . 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 10000 and 15000 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.

### 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 that the solar values.

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

## Abundance variations with magnetic field strength

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.  plots abundances versus $\log B_{\rm 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_{\rm 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 Figure  are those expected from the known decline of $\langle B_z \rangle$ 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  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 . This effect has also been exploited recently by the Special Astrophysical Observatory (SAO) group to search for particularly large magnetic fields in Ap/Bp stars.

## Cr and Fe abundances versus effective temperature

Our results can be compared to previous studies of magnetic Bp stars. 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. added stars between about 10000 and 15000 K and our current study, as well as a previous study by , adds 23 stars in that same temperature range. The results are shown in Fig. . Our observations confirm the tendency, already visible in the data of Ryabchikova, for the star with $T_{\mathrm{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_{\mathrm{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_{\mathrm{eff}}$. These results support the conclusions of that the maximum abundance for Cr and Fe both approach the same value of between $-3$ and $-4$ dex.

# 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_\odot$ star through its main sequence life. Figure  shows a Hertzsprung-Russell (HR) diagram for the stars in this study, compared to theoretical stellar evolutionary models .

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

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_{\mathrm{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_{\mathrm{eff}}$ above about 7000 K, the upward acceleration $g_{\rm R}$ 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 $g_{\rm 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 , and .

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

The time variation with height of abundance of a very unabundant element diffusing up through a stellar atmosphere has recently been studied by , 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 \sim 4$ and $\ensuremath{T_{\mathrm{eff}}}\sim 10-13\,000$ K. In the panels of Figure , it appears that we can identify several different cases.

## Abundant light elements: He, O

The general behaviour of these elements was already predicted by . For the very abundant light elements He and O, the radiative acceleration for $\ensuremath{T_{\mathrm{eff}}}\sim 12000$ K is much too weak to support a solar abundance of these elements (cf. $g_{\rm rad}$ calculations for $\ensuremath{T_{\mathrm{eff}}}= 12000$ K by ). 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 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 , and ; 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 ).

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_{\mathrm{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.

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_{\rm rad}$ value for O at $12000$ K found by . 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_{\mathrm{eff}}$ and $\log g$ both decrease. The change in $T_{\mathrm{eff}}$ typically reduces $g_{\rm R}$. However, the decrease in $\log g$ means that a smaller value of $g_{\rm 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.5 M_\odot$ from $\ensuremath{T_{\mathrm{eff}}}= 13000$ K and $\log g = 4.5$ to $\ensuremath{T_{\mathrm{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.

## 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_{\rm Mg}/n_{\rm H}) \approx -5.2$. This may be compared with the predicted Mg abundance profile for a star of $\ensuremath{T_{\mathrm{eff}}}= 12000$ K of and of . The abundance predicted by the equilibrium calculation ($g_{\rm R} = 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_{\mathrm{eff}}$ and $\log g$ for our $3.5 M_\odot$ stars.

Silicon equilibrium atmospheric abundance has been studied by a number of authors, including , , and . While detailed results differ, the overall conclusion is that Si is expected to be of order 1 dex underabundant in the atmosphere at $\ensuremath{T_{\mathrm{eff}}}\sim 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 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.

## 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 , , and (for Fe) by . 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_{\mathrm{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_{\mathrm{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_\odot$ star

The available equilibrium calculations for Ti and Cr in the mass range of our observations by 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_{\rm el}/n_{\rm H})$ similar to those of Fe, in the range of $-3$ to $-4$. This is confirmed by the calculations of . 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 . 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_{\mathrm{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 . It would be of great interest to have computations of equilibrium abundances of Ti and Cr following the evolution of a $3.5 M_\odot$ star for more direct comparison with our results.

## 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_{\mathrm{eff}}$ temperature range of interest to us here, except for the equilibrium stratification calculated for a simplified artificial low-abundance element with $\ensuremath{T_{\mathrm{eff}}}= 12000$ K by . 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 . 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 upward is **probably** declining with decreasing $T_{\mathrm{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.

# 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$^{4}$ 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$_{\odot}$ 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_{\mathrm{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$_{\odot}$) and less massive stars (2-3 M$_{\odot}$) to see if similar trends are observed.

[^1]: Based in part on observations made with the European Southern Observatory (ESO) telescopes under the ESO programmes 072.D-0410(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 lÚniverse of the Centre National de la Recherche Scientifique of France and the University of Hawaii.