text:
# Spectrum and Regge-trajectories in QCD

## 1 Introduction

Starting in sixties an active group of physicists under the guidance of Prof. K.A.Ter-Martirosyan was creating the theory of high-energy processes in QCD. From the beginning the key element of this theory is the notion of Regge trajectories and in particular of the pomeron trajectory, which have been introduced phenomenologically. In this talk I review the problem of spectrum and Regge trajectories as it can be derived from the nonperturbative QCD dynamics.

The problem of the QCD spectrum is the central issue in the nonperturbative QCD and is intimately connected to problems of confinement and mass generation in QCD. These latter issue make QCD so much different from QED, in addition explicit form of Nonperturbative Interaction (NPI) was not known for a long time. Recently with the introduction of the Field Correlator Method (FCM)  the situation has changed favorably, since the NPI is defined there in a gauge-invariant way and the simple field correlators, which are sufficient for all dynamical calculations are known from lattice data  and analytic results .

In view of all this one can put the problem of the QCD spectrum in the most general framework, as the many-channel problem of bound states of quarks, antiquarks and valence gluons with all possible mixings, the states being stable in the limit $`N_c\mathrm{}`$, and acquiring decay widths when $`N_c=3`$. It is the purpose of this talk to describe the QCD spectrum in the Hamiltonian language, using FCM and the $`N_c=\mathrm{}`$ limit. In addition FCM provides another convenient limit – the gluon correlation length $`T_g`$ tending to zero, while the string tension $`\sigma `$ is kept constant. As a result the universal Hamiltonian becomes local and has a transparent structure for any number of constituents. As will be seen, it contains only two fixed input parameters:$`\sigma `$ and $`\alpha _s`$ (or $`\mathrm{\Lambda }_{QCD})`$ in addition to current quark masses and is able to predict any meson, baryon, hybrid and glueball states and their mixings (for a review see ).

In doing so one solves the following puzzles: 1) constituent mass of quarks and gluons 2) meson Regge trajectories with a correct slope and intercept 3) the explicit notion of the valence gluon 4) hybrid spectra and Regge trajectories.

Remarkably all calculated spectra are in very good agreement with lattice data and experiment which gives an additional justification for FCM. However immediate consequences of the whole approach are far more reaching. They include a new formulation of the perturbation theory, the so-called Background Perturbation theory (BPTh)  with $`\alpha _s`$ saturating at small Euclidean momenta instead of diverging near the Landau ghost pole. The whole structure of QCD becomes interconnected with the spectrum via the quark-hadron duality and hybrids and glueballs play a very important role in scattering and decay.

## 2 Hamiltonian

There are two possible approaches to incorporating nonperturbative field correlators in the quark-antiquark (or $`3q`$) dynamics. The first has to deal with the effective nonlocal quark Lagrangian containing field correlators . From this one obtains first-order Dirac-type integro-differential equations for heavy-light mesons , light mesons  and baryons . these equations contain the effect of chiral symmetry breaking  which is directly connected to confinement.

The second approach is based on the effective Hamiltonian for any gauge-invariant quark-gluon system. In the limit $`T_g0`$ this Hamiltoniam is simple and local, and in most cases when spin interaction can be considered as a perturbation one obtains results for the spectra in an analytic form, which is transparent.

For this reason we choose below the second, Hamiltonian approach . We start with the exact Fock-Feynman-Schwinger Representation for the $`q\overline{q}`$ Green’s function (for a review see ), taking for simplicity nonzero flavor case

$$G_{q\overline{q}}^{(x,y)}=_0^{\mathrm{}}𝑑s_1_0^{\mathrm{}}𝑑s_2(z)_{xy}(D\overline{z})_{xy}r^{K_1K_2}tr\mathrm{\Gamma }_{in}(m_1\widehat{D}_1)W_\sigma (C)\mathrm{\Gamma }_{out}(m_2\widehat{D}_2)_A$$
(1)
where $`K_i=_0^{s_1}𝑑\tau _i(m_i+\frac{1}{4}(\dot{z}_\mu ^{(i)})^2),`$ $`\mathrm{\Gamma }_{in,out}=1,\gamma _5,\mathrm{}`$ are meson vertices, and $`W_\sigma (C)`$ is the Wilson loop with spin insertions, taken along the contour $`C`$ formed by paths $`(Dz)_{xy}`$ and $`(D\overline{z})_{xy}`$,

$$W_\sigma (C)=P_FP_A\mathrm{exp}(ig_CA_\mu 𝑑z_\mu )\mathrm{exp}(g_0^{s_1}\sigma _{\mu \nu }^{(1)}F_{\mu \nu }𝑑\tau _1g_0^{s_2}\sigma _{\mu \nu }^{(2)}F_{\mu \nu }𝑑\tau _2).$$
(2)
The last factor in (2) defines the spin interaction of quark and antiquark. The average $`W_\sigma _A`$ in (1) can be computed exactly through field correlators $`F(1)\mathrm{}F(n)_A`$, and keeping only the lowest one,$`F(1)F(2)`$, which yields according to lattice calculation  accuracy around 1%,  one obtains

$$W_\sigma (C)_A\mathrm{exp}(\frac{1}{2}[_{S_{min}}ds_{\mu \nu }(1)_{S_{min}}ds_{\lambda \sigma }(2)+$$

$$+\underset{i,j=1}{\overset{2}{}}_0^{s_i}\sigma _{\mu \nu }^{(i)}d\tau _i_0^{s_j}\sigma _{\lambda \sigma }^{(j)}d\tau _j]F_{\mu \nu }(1)F_{\lambda \sigma }(2)).$$
(3)
The Gaussian correlator $`F_{\mu \nu }(1)F_{\lambda \sigma }(2)D_{\mu \nu ,\lambda \sigma }(1,2)`$ can be rewritten identically in terms of two scalar functions $`D(x)`$ and $`D_1(x)`$ , which have been computed on the lattice  to have the exponential form $`D(x)D_1(x)\mathrm{exp}(|x|/T_g`$ with the gluon correlation length $`T_g0.2`$ fm.

Thus the first term in the exponent (3) yields the area law $`W_\sigma `$ $`\mathrm{exp}(\sigma S_{\mathrm{min}})`$, with the string tension $`\sigma `$

$$\sigma =\frac{1}{2}D(x)d^2x.$$
(4)
We concentrate now on this confining term in (3) and quote the result for the spin-dependent term at the end of this section.

The important point to be discussed now is the character of dynamics one gets from (3), (4). To this end one should compare two characteristic lengths (times) $`T_g`$ and $`T_q`$ , where $`T_q`$ is the typical period of quark motion, e.g. the classical period of motion along the Coulomb orbit (for heavy quarks) or along the orbit in the linear potential (for light quarks). In all cases one gets $`T_q\begin{array}{c}>\hfill \\ \hfill \end{array}1`$ fm and hence $`T_qT_g`$. Thus one has a local dynamics however relativistic for light quarks, but in any case the excitation of gluonic vacuum degrees of freedom can be neglected in the first approximation, so that the dynamics is of the potential type.

As the next step one introduces the einbein variables $`\mu _i`$ and $`\nu `$; the first one to transform the proper times $`s_i,\tau _i`$ into the actual (Euclidean) times $`t_iz_4^{(i)}`$. One has

$$2\mu _i(t_i)=\frac{dt_i}{d\tau _i},_0^{\mathrm{}}𝑑s_i(D^4z^{(i)})_{xy}=constD\mu _i(t_i)(D^3z^{(i)})_{xy}.$$
(5)
The variable $`\nu `$ enters in the Gaussian representation of the Nambu-Goto form for $`S_{min}`$ and its stationary value $`\nu _0`$ has the physical meaning of the energy density along the string. In case of several strings, as in the baryon case or the hybrid case, each piece of string has its own parameter $`\nu ^{(i)}.`$

To get rid of the path integration in (1) one can go over to the effective Hamiltonian using the identity

$$G_{q\overline{q}}(x,y)=x|\mathrm{exp}(HT)|y$$
(6)
where $`T`$ is the evolution parameter corresponding to the hypersurface chosen for the Hamiltonian: it is the hyperplane $`z_4=const`$ in the c.m. case  and $`z_+=const`$ in the light-cone case .

The final form of the c.m. Hamiltonian (apart from the spin and perturbative terms to be discussed later) for the $`q\overline{q}`$ case is

$$H_0=\underset{i=1}{\overset{2}{}}\left(\frac{m_i^2+𝐩_i^2}{2\mu _i}+\frac{\mu _i}{2}\right)+\frac{\widehat{L}^2/r^2}{2[\mu _1(1\zeta )^2+\mu _2\zeta ^2+_0^1𝑑\beta (\beta \zeta )^2\nu (\beta )]}+$$

$$+\frac{\sigma ^2r^2}{2}_0^1\frac{d\beta }{\nu (\beta )}+_0^1\frac{\nu (\beta )}{2}𝑑\beta .$$
(7)
Here $`\zeta =(\mu _1+_0\beta \nu 𝑑\beta )/(\mu _1+\mu _2+_0^1\beta \nu 𝑑\beta )`$ and $`\mu _i`$ and $`\nu (\beta )`$ are to be found from the stationary point of the Hamiltonian

$$\frac{H_0}{\mu _i}|_{\mu _i=\mu _i^{(0)}}=0,\frac{H_0}{\nu }|_{\nu =\nu ^{(0)}}=0.$$
(8)
Note that $`H_0`$ contains as input only $`m_1,m_2`$ and $`\sigma `$, where $`m_i`$ are current masses defined at the scale 1 GeV. The further analysis is simplified by the observation that for $`L=0`$ one finds $`\nu ^{(0)}=\sigma r`$ from (8) and $`\mu _i=\sqrt{m^2+𝐩^2}`$, hence $`H_0`$ becomes the usual Relativistic Quark Model (RQM) Hamiltonian

$$H_0(t=0)=\underset{i=1}{}\sqrt{m_i^2+𝐩^2}+\sigma r.$$
(9)
For large $`L`$ however one can neglect $`\mu _i`$ as compared to $`\nu `$ and one gets

$$H_0^22\pi \sigma \sqrt{L(l+1)},\nu ^{(0)}(\beta )=\sqrt{\frac{8\sigma l}{\pi }}\frac{1}{\sqrt{14(\beta \frac{1}{2})^2}}.$$
(10)
From (9), (10) one can see that $`\nu ^{(0)}(\beta )`$ is indeed the energy density along the string and $`\mu ^{(0)}`$ is the c.m. energy of the quark, which plays the role of constituent quark mass as will be seen below.

To proceed one can use two approximations. First, replace $`H_0`$ in (8) by its eigenvalue $`M_0`$, which is accurate within 5% . Second, approximate the $`L`$-dependent term in (7), introducing the correction $`\delta H_{str}`$ namely

$$H_0H_R+\mathrm{\Delta }H_{str},H_R=\underset{i=1}{\overset{2}{}}\left(\frac{𝐩^2+m_i^2}{2\mu _i}+\frac{\mu _i}{2}\right)+\sigma r$$
(11)
and the mass correction due to $`\mathrm{\Delta }H_{str}`$ for equal mass case is

$$\mathrm{\Delta }_{str}(nL)=\mathrm{\Delta }H_{str}=\frac{16}{3}\frac{\sigma ^2L(L+1)}{M_0^3}$$
(12)
where $`M_0`$ is the eigenvalue of $`H_R`$; a more accurate approximation is given by

$$\mathrm{\Delta }_{str}(nL)=\frac{2\sigma L(L+1)1/r}{M_0^2}.$$
(13)
But $`H_0`$ is not the whole story, one should take into account 3 additional terms: spin terms in (3) which produce two types of contributions: self-energy correction

$$H_{self}=\underset{i=1}{\overset{2}{}}\frac{\mathrm{\Delta }m_q^2(i)}{2\mu _i},\mathrm{\Delta }m_q^2=\frac{4\sigma }{\pi }\eta (m_i),\eta (0)1÷0.9$$
(14)
and spin-dependent interaction between quark and antiquark $`H_{spin}`$  which is entirely described by the field correlators $`D(x),D_1(x)`$, including also the one-gluon exchange part present in $`D_1(x)`$.

Finally one should take into account gluon exchange contributions , which can be divided into the Coulomb part $`H_{Coul}=\frac{4}{3}\frac{\alpha _s(r)}{r},`$ and $`H_{rad}`$ including space-like gluon exchanges and perturbative self-energy corrections (we shall systematically omit these corrections since they are small for light quarks to be discussed below). In addition there are gluon contributions which are nondiagonal in number of gluons $`n_g`$ and quarks (till now only the sector $`n_g=0`$ was considered) and therefore mixing meson states with hybrids and glueballs . we call these terms $`H_{mix}`$ and refer the reader to  and the cited there references for more discussion. Assembling all terms together one has the following total Hamiltonian in the limit of large $`N_c`$ and small $`T_g`$:

$$H=H_0+H_{self}+H_{spin}+H_{Coul}+H_{rad}+H_{mix}.$$
(15)
We start with $`H_0=H_R+H_{string}`$. The eigenvalues $`M_0`$ of $`H_R`$ can be given with 1% accuracy by

$$M_0^28\sigma L+4\pi \sigma (n+\frac{3}{4})$$
(16)
where $`n`$ is the radial quantum number, $`n=0,1,2,\mathrm{}`$ Remarkably $`M_04\mu _0`$, and for $`L=n=0`$ one has $`\mu _0(0,0)=0.35`$ GeV for $`\sigma =0.18`$ GeV2, and $`\mu _0`$ is fast increasing with growing $`n`$ and $`L`$. This fact explains that spin interactions become unimportant beyond $`L=0,1,2`$ since they are proportional to $`d\tau _1d\tau _2\frac{1}{4\mu _1\mu _2}dt_1dt_2`$ (see (3) and ). Thus constituent mass (which is actually ”constituent energy”) $`\mu _0`$ is ”running”. The validity of $`\mu _0`$ as a socially accepted ”constituent mass” is confirmed by its numerical value given above, the spin splittings of light  and heavy-light mesons  and by baryon magnetic moments expressed directly through $`\mu _0`$, and being in agreement with experimental values .

The next topic is Regge trajectories in QCD. As it is clear from (10), one has the correct asymptotic Regge slope coinciding with the string picture, while for small $`L,L4`$, the approximation (11), (12) holds  with almost the same slope. The Regge intercept depends strongly on the term $`H_{self}`$ , since (16) yields too large value for $`M_0`$, e.g. $`M_0(0,0)1.4`$ GeV is almost twice the $`\rho `$ mass. However the self-energy term (14) defined unambigiously through $`\sigma `$  has a negative sign and a magnitude which brings the mass back near the experimental value. Thus one can understand the origin of the large negative phenomenological constant, which is usually introduced in RQM, but it is also rewarding that it is not actually constant, but depends on $`n,L`$ via $`\mu _0(nL)`$ so that the linear Regge behaviour is preserved.

To compare with the experiment and disentangle the contribution of spin interaction we shall consider the center-of-gravity (cog) masses for each meson multiplet as in . Then the masses of all orbital excitations $`(n=0)`$ can be nicely described by the linear Regge trajectory, which we call the Regge $`L`$-trajectory with experimental parameters

$$\overline{M}^2(L)=(1.23\pm 0.02)L+0.37\pm 0.02(\mathrm{GeV}^2)$$
(17)
or

$$L=0.81\overline{M}^2(L)0.30.$$
(18)
This is different from the leading $`\rho `$-trajectory

$$J=\alpha _J^{}M^2(J)+0.48,\alpha _J^{}=0.88\mathrm{GeV}^2,\alpha _J(0)=0.48$$
(19)
since its parameters depend on spin interactions. Now using (11-16) one obtains

$$L=0.80\overline{M}^2(L)0.34$$
(20)
which agrees with experiment (18) within 10%, the accuracy being in accordance with the estimates of the neglected terms in (15), namely $`H_{rad}`$ and $`H_{mix}`$. In this way one solves the problem of Regge trajectories for orbital excitations in QCD, thus supporting the foundations o of the wide and fruitful activity undertaken by Professor K.A.Ter-Martirosyan and his group to describe the high-energy scattering and production processes in the framework of the Regge theory. This important contribution was reviewed in .

We now come to the gluon-containing systems, hybrids and glueballs. Referring the reader to the original papers - one can recapitulate the main results for the spectrum. In both cases the total Hamiltonian has the same form as in (15), however the contribution of corrections differs.

For glueballs it was argued in  that $`H_0`$ (11) has the same form, but with $`m_i=0`$ and $`\sigma \sigma _{adj}=\frac{9}{4}\sigma `$ while $`H_{self}=0`$ due to gauge invariance, and $`H_{coul}`$ is small due to strong cancellation between tree graps and loop corrections . Thus glueball masses are expressed through only $`\sigma =0.18`$ GeV2 (fixed by meson Regge trajectories) for the center-of-gravity and in addition through $`\alpha _s`$ for the spin splittings. One can see in Table 1 of  the theoretical c.o.g. mass values computed in  in comparison with lattice data. The agreement is striking, especially if one takes into account that in theoretical calculations there are no parameters at all –$`\sigma `$ was fixed beforehand at the same value, as in lattice, and $`\alpha _s`$ was neglected altogether for the reason stated above.

We now come to the delicate and very important topic of glueball trajectories and especially of the pomeron trajectory. Since the glueball Hamiltonian is basically the same as for mesons, one expects that the asymptotic slope of all Regge trajectories would be

$$\alpha _G^{}(M^2\mathrm{})\frac{1}{2\pi \sigma _{adj}}=\frac{4}{9}\alpha _M^{}.$$
(21)
One expects that the pomeron trajectory is passing through the states $`2^{++}`$ (2.29 GeV for $`\sigma =0.18`$ GeV2) and $`4^{++}`$ (around 3.2 GeV) in which case the pomeron intercept appears too low $`\alpha _P(0)<1`$. A possible way out was suggested in , where the intersection of two lowest meson trajectories with the pomeron trajectory was introduced yielding the correct intercept $`\alpha _p(0)=1.1.÷1.2`$ for reasonable parameters of trajectory interactions. However this topic is far from clarity, and the BFKL perturbative results for $`\alpha _P(0)`$ seem to be not stable with respect to inclusion to higher orders and/or nonperturbative effects. It seems that the pomeron trajectory depends on both perturbative and nonperturbative contributions and their possible interference, and the problem was never considered in that fullness.

We now coming to the final topic of this talk -hybrids and their role in hadron dynamics. We start with the hybrid Hamiltonian and spectrum. This topic in the framework of FCM was considered in  The Hamiltonian $`H_0`$ for hybrid looks like

$$H_0^{(hyb)}=\frac{m_1^2}{2\mu _1}+\frac{m_2^2}{2\mu _2}+\frac{\mu _1+\mu _2+\mu _g}{2}+\frac{𝐩_\xi ^2+𝐩_\eta ^2}{2\mu }+\sigma \underset{i=1}{\overset{2}{}}|𝐫_g𝐫_i|+H_{str}.$$
(22)
Here $`𝐩_\xi ,𝐩_\eta `$ are Jacobi momenta of the 3-body system, $`H_{self}`$ is the same as for meson, while $`H_{spin}`$ and $`H_{coul}`$ have different structure .

The main feature of the present approach based on the BPTh, is that valence gluon in the hybrid is situated at some point on the string connecting quark and antiquark, and the gluon creates a kink on the string so that two pieces of the string move independently (however connected at the point of gluon). This differs strongly from the flux-tube model where hybrid is associated with the string excitation as a whole.

The difference between two approaches is especially pronounced in the case of the hybrid with static quarks  where the flux-tube model predicts for large $`R`$ the terms $`E_n(R)\frac{\pi n}{R}`$, while in FCM there is another branch corresponding to the longitudinal d.o.f. of gluon $`E(R)\frac{const}{R^{1/3}}`$ , and the lattice results certainly prefer the latter and contradict the flux-tube asymptotics. Also at intermediate interquark distances $`R`$ the spectrum of FCM approach is in much better agreement than only another model, see  for details and discussion.

Results for light and heavy exotic $`1^+`$ hybrids are also given in  and are in agreement with lattice calculations. Typically an additional gluon in the exotic $`(L=1)`$ state ”weights” 1.2$`÷`$1.5 GeV for light to heavy quarks, while nonexotic gluon $`(L=0)`$ brings about 1 GeV to the mass of the total $`q\overline{q}g`$ system.

It is understandable now that hybrids play a very special role in QCD, namely they describe the excitation of the film – the string world sheet – which is between $`q`$ and $`\overline{q}`$ and in particular covers the gluon and quark loops appearing in the $`\alpha _s`$ renormalization. The ”additional gluon” mass $`M_g1`$ GeV2 cited above enters as a screening mass in the one-loop running $`\alpha _s`$ (see   for details and derivations

$$\alpha _s(Q)=\frac{4\pi }{b_0\mathrm{ln}\left(\frac{Q^2+M_g^2}{\mathrm{\Lambda }_{QCD}^2}\right)}$$
(23)
and this form (in $`r`$-space) is in a perfect agreement with a recent calculation of $`\alpha _s`$ on the lattice .

Thus one can say that the perturbative QCD in the IR region is defined by the hybrid physics. A similar conclusion can be drawn with respect to the DIS diagrams at small to moderate $`Q^2,Q^2\begin{array}{c}<\hfill \\ \hfill \end{array}1`$ GeV2, where familiar ladder diagrams with gluon exchanges are replaced by (multi)-hybrid Green’s functions, and to the hadron-hadron scattering, where e.g. a sequence of BFKL ladders is replaced by multihybrid diagrams. The progress in this direction is highly desirable.

All discussion above refers to the leading large $`N_c`$ limit, which has a remarkable accuracy, as it was recently demonstrated on the lattice  and by comparison of calculated and experimental masses in , , , . We now turn to the $`1/N_c`$ effects which as shown in  strongly modify the masses of radially excited mesons. Typically those mesons have radii larger than 1.5 fm and as it was argued in  at the interquark distances $`R1.4`$ fm there appear sea-quark holes in the confining film connecting $`q`$ and $`\overline{q}`$ trajectories. It was assumed in  that the film with the virtual holes establishes the quasistationary state which can be called the ”predecay state”. It is clear that the effective string tension decreases due to holes and the confining potential is partly screened at $`R1.4`$ fm. The masses of radially excited mesons for $`L=0,1,2,3`$ and $`n_r=0,1,2,3,4`$ have been calculated in  with the help of this potential and are shown in Fig.1. One can see a remarkable agreement of theoretical masses (black dots) with experimental candidates. Another interesting feature is the almost exact linearity of trajectories $`M^2(n_r,L)`$ vs $`n_r`$. The effect of the mass decrease due to the sea-quark holes is significant (for high $`n_r`$ it is around $`0.5`$ GeV).

Concluding the talk I would like to stress the simplicity of the method (FCM) which solves at least 4 important problems in the QCD spectrum: 1) Regge slope 2) Regge intercept 3) the problem of constituent masses 4) problem of radial Regge trajectories.

The author is indebted for useful discussions to A.B.Kaidalov, Yu. S. Kalashnikova. During all years of the development of the present method the author felt a constant interest, a helpful attention and support from K.A.Ter-Martirosyan, which are very gratefully acknowledged. The present work was partially supported by RFBR grants 00-02-17836 and 00-15-96786 and by INTAS grants 00-110 and 00-00366.


format: text

text:
# Automated Spectroscopic abundances of A and F-type stars using echelle spectrographs

## 1 Introduction

A-F type stars present a large number of chemical peculiarities: Ap, Am, $`\lambda `$ Bootis, $`\rho `$ Puppis (see Gray & Garrison 1989b , for a discussion of the $`\delta `$ Del stars and a definition of the $`\rho `$ Pup stars). Radiative diffusion is generally invoked most of these peculiarities, and the theory of this mechanism, developed mostly by Michaud and collaborators (Michaud michaud (1991)) seems indeed largely successful, even though some adjustments seem to remain necessary. For instance, stellar wind was added to radiative diffusion in order to make more realistic predictions of abundances of Am stars (Michaud et al. michaud83 (1983)).

A little explored category of chemically peculiar stars in this range of spectral types is that of metallic giants, discovered by Hauck (hauck (1986)) on the basis of textscGeneva photometry and further studied by Berthet (berthet90 (1990, 1991)) using spectroscopy. They are A-F stars leaving the main sequence and having the Geneva photometric parameter $`\mathrm{\Delta }\mathrm{m}_20.013`$ (see Hauck, hauck73 (1973) for the definition of $`\mathrm{\Delta }\mathrm{m}_2`$). Berthet (berthet90 (1990, 1991)) confirmed that these stars have metal abundances higher than solar ones, as suggested by their large $`\mathrm{\Delta }\mathrm{m}_2`$ values. However, he analyzed only metallic giants and relied on published abundances of normal giants for comparison. Berthet found that the abundance pattern of metallic giants is similar to that of Am stars, except for calcium and scandium, which are notoriously underabundant in the latter, while they are roughly normal or slightly overabundant in the former. On this basis, Berthet (berthet92 (1992)) proposed that metallic giants might be evolved Am stars, an idea theoretically justified by the following scenario: in the beginning of the life of an Am star, Ca and Sc sink and accumulate below the line forming region because of their too small radiative acceleration. During stellar evolution on the main sequence, however, the effective temperature decreases and the outer convection zone develops and thickens, until it reaches the depth where Ca and Sc have accumulated and dredges them up, so that their surface abundances becomes normal again. Künzli & North (kunzli98 (1998)) had challenged this scenario, arguing, in particular, that the $`V\mathrm{sin}i`$ distribution of metallic giants is not compatible with that of Am stars, especially evolved ones. The fact that metallic giants are more evolved and more massive than normal ones (North et al. north98 (1998)) does not fit Berthet’s scenario either.

The initial purpose of this work was to focus our attention on metallic giants by comparing their abundances with those of normal giants, in order to better define their characteristics and determine how far they form a well defined spectroscopic group. However, since Berthet’s scenario remains unsatisfactory, the question of their origin remains open, so we decided to examine these stars in a more general context and try to understand what are their progenitors. So we judged important to study not only giant A-F stars (normal and metallic) but also their possible progenitors, i.e. normal A-F type stars of similar masses on the main sequence. To this purpose, we chose to observe a sample of A-F stars that cover the HR diagram as homogeneously as possible, in order to map the abundance pattern in the whole region of interest. In addition to the homogeneous coverage of the HR diagram, it is important to ensure also the homogeneity of the method of spectral analysis. In the latter sense, homogeneity is crucial, since it is very delicate to compare abundance analyses done by different authors. There is always at least one assumption that differs. The most common problem is the choice of atmospheric parameters, $`T_{\mathrm{eff}}`$ and $`\mathrm{log}g`$. Indeed, different techniques, using different photometric systems or spectroscopy, are used to derive atmospheric parameters. Abundance differences arising from these techniques are generally not very important since they remain within error boxes, but they significantly increase the scatter in the results. Another common difference is the choice of atmosphere model theory. In the field of A-F stars, ATLAS9 models from Kurucz (kurucz1 (1993)) are the rule. However, these models exist with and without overshooting, and both are used depending on the authors. Finally, methods to derive the abundances differ. Some people use equivalent widths while others prefer synthetic spectra adjustments. The oscillator strengths and other line parameters may come from different sources and these sources may vary from one element to the next. The adopted solar iron abundance may also differ from one author to the other. Therefore, even for iron, a large scatter is found. For example, even for a very bright star such as Procyon, \[Fe/H\] values vary from -0.18 to 0.05 (Cayrel de Strobel et al. cayrel (2001)). All this shows how well justified is the extensive and homogeneous analysis of a large sample of A-F type stars that we present in this paper.

In Sect. 2, we describe the observations and the reduction procedure, while Sect. 3 gives more details on the definition of the sample. The determination of atmospheric parameters is described in Sect. 4, and sources of uncertainties on the derived abundances are discussed in Sect. 5. The results are presented in Sect. 6 and discussed in Sect. 7.

## 2 Observations and data reduction

The spectra used in this work were obtained with two almost identical spectrographs, ELODIE and CORALIE (see Baranne et al., baranne (1996), for technical details about ELODIE and Queloz et al., queloz (2001), for CORALIE). Both spectrographs have been built as a joint project between the Geneva Observatory and the Observatoire de Haute-Provence. ELODIE is attached to the 1.93 m telescope of the Observatoire de Haute-Provence whereas CORALIE is attached to the Swiss 1.2 m Leonard Euler telescope at the ESO la Silla Observatory. The optical design of both spectrographs is the same. The major difference is the better resolving power of CORALIE (50 000 compared to 42 000 for ELODIE) that comes from a slightly different optical combination at the entrance of the spectrograph and the use of a 2k$`\times `$2k CCD camera (against a 1k$`\times `$1k for ELODIE) with smaller pixels (15 $`\mu `$m vs 24 $`\mu `$m) (Queloz et al. queloz99 (1999)).

Two observing runs were held in June and August 1999 in Haute-Provence, and one in July 2000 in Chile for a total of 23 nights (10 and 13 respectively). Target stars were bright ones, because we wanted a high signal to noise ratio (S/N $`>`$ 200). The total number of observed stars is 140, after elimination of all spectroscopic binaries or stars presenting strange cross-correlation functions due e.g. to pulsations.

Data reduction was done as presented in Paper I. As the optical design of CORALIE is the same as the one of ELODIE, only a few changes in the parameters used to characterize the CCD were needed to adapt the IRAF reduction procedure which was initially designed for ELODIE spectra. However, after order extraction, the different orders merge perfectly without scaling, contrary to the case of ELODIE spectra. The whole spectrum is obtained by a simple concatenation of consecutive orders with a detection of the overlapping region. As the S/N ratio is better in the red part of the order than in the blue one, the flux of the first order is retained for 3/4 of the overlapping region, and the flux of the following for the remaining 1/4. The subsequent operations are done using the same procedure as described in Paper I.

## 3 Sample

The mass range is 1.3 to about 2.8 M, encompassing that of most metallic giants (1.7 to 3 M). As we wanted to observe as many stars as possible during the observing runs at our disposal, our choice naturally went to the brighter stars. Two additional requirements were used to define the sample: availability of Geneva photometry (Golay G80 (1980); Rufener & Nicolet RN88 (1988); Cramer C99 (1999)) and of Hipparcos parallaxes (ESA esa (1997)), in order to be able to determine the atmospheric parameters with the method presented in Sect. 4. Known binaries were exluded too. Finally, we chose to keep only stars with $`V\mathrm{sin}i`$ $`<150`$ $`\mathrm{km}\mathrm{s}^1`$ or that have no determination in the bright star catalogue (Hoffleit & Jaschek bsc (1991)). These were the only restrictions before observing.

During the observing runs, as our main targets were metallic giants, we focused a little more on them and on their progenitors, i.e. on stars falling in the mass range 1.7 to 2.5 M. Low mass stars had been included in the sample essentially because the early Montréal models (see Turcotte et al. (turcotte98 (1998)), Richer et al. (richer00 (2000)) and Richard et al. (richard01 (2001)) and reference therein) did not consider masses larger than 1.5 M.

As observations were done with spectrographs designed to get radial velocities, the cross-correlation peak was routinely available. we used it to exclude stars presenting a “dubious” peak, i.e. an asymmetric one that may come either from binarity or from pulsation, or a very broad one resulting from $`V\mathrm{sin}i`$ larger than 200 $`\mathrm{km}\mathrm{s}^1`$.

The final sample is composed of 140 stars, out of which 35 are giants (18 metallic and 17 normal), 10 are Am and 6 are Ap. Spectral types are taken from Hauck (hauck (1986)) for stars in common, and from Gray & Garrison (gray87 (1987); 1989a ; 1989b ) otherwise. The chemically peculiar Am and Ap stars are included only because their positions in the HR diagram correspond to the studied region. It is interesting that without taking care of this detail, the proportion of Am and Ap stars, relatively to stars of the main sequence, is representative.

## 4 Atmospheric parameters determination

Atmospheric parameters were determined from Geneva photometry for $`T_{\mathrm{eff}}`$ and from Hipparcos data (ESA, esa (1997)) for $`\mathrm{log}g`$. The calib code (Künzli et al., kunzli (1997)) was used to derive effective temperature from textscGeneva photometry. Although this code also provides $`\mathrm{log}g`$, we prefer to determine that value from Hipparcos data. All stars of the sample are bright and have an accurate trigonometric parallax measurement by Hipparcos. They all have a parallax larger than 7.4 mas (i.e. are closer than 135 pc except two, HD 175510 and HD 187764 that are at 166 and 171 pc respectively). Therefore reddening was neglected 111Although this is completely justified for the 121 objects closer than 70 pc on the basis of $`uvby\beta `$ photometry, a closer examination requested by the referee shows that the remaining objects are slightly more reddened on average. In particular, 4 stars have $`E(by)>0.020`$: HD 115604 ($`E(by)=0.025`$), HD 181333 (0.035), HD 200723 (0.023) and HD 214441 (0.027), implying that $`T_{\mathrm{eff}}`$ might be underestimated by $`300`$ K in the worst case. Among the 121 nearby stars, 8 have also $`E(by)>0.020`$, but they constitute the queue of a fairly gaussian distribution, 6 of them being closer than 50 pc. The colour excesses were estimated using the usual calibrations by Crawford (C75 (1975), C79 (1979)), incorporated in a code kindly made available by Dr. C. Jordi. Bolometric luminosity was calculated using parallax, bolometric correction interpolated in Flower’s table (flower (1996)) and $`M_{\text{bol}}=4.75`$ for the Sun. The bolometric luminosity is then given by

$$\mathrm{log}\frac{L}{L_{}}=0.4(M_\text{v}4.75+B.C.)$$

The $`\mathrm{log}g`$ value was then interpolated in the grids of models of Schaller et al (schaller (1992)), Schaerer et al. (1993a ), Charbonnel et al. (charbonnel (1993)) and Schaerer et al. (1993b ) according to the method described by North et al. (north (1997)). Finally, an ATLAS9 atmosphere model was computed without overshooting (Castelli et al., castelli (1996)) with temperature rounded to 25K and $`\mathrm{log}g`$ to 0.05 dex using the version of ATLAS9 adapted by M. Lemke222http://www.sternwarte.uni-erlangen.de/ftp/michael/ atlas-lemke.tgz for UNIX system.

Derived parameters of our sample stars are presented in Table 5 available at CDS. It contains HD, $`T_{\mathrm{eff}}`$, $`\mathrm{log}g`$, $`V_{\mathrm{rad}}`$, $`V\mathrm{sin}i`$ and $`\xi _{\mathrm{micro}}`$. The last three values are those derived with our method. Fig. 1 shows their position in the HR diagram. Stars cover the whole main sequence and extend slightly beyond.

## 5 Uncertainty estimation

### 5.1 Errors resulting from wrong atmospheric parameters

In order to check the influence of erroneous atmospheric parameters ($`T_{\mathrm{eff}}`$ $`\pm 150K`$, $`\mathrm{log}g`$ $`\pm 0.2`$ dex) on derived abundances, simulations were done on a synthetic spectrum. These values do not represent the internal error of the methods used to derive these parameters, but rather the scatter in the estimations done with different methods (photometric or spectroscopic). The synthetic spectrum was computed with a model atmosphere of $`T_{\mathrm{eff}}`$ = 6950 K, $`\mathrm{log}g`$ = 4.0 and solar chemical composition. This model corresponds to a typical star of our sample lying in the middle of both temperature and gravity ranges. The other parameters, namely radial, microturbulent and especially rotational velocities, are $`V_{\mathrm{rad}}`$ = 0 $`\mathrm{km}\mathrm{s}^1`$, $`\xi _{\mathrm{micro}}`$ = 1.5 $`\mathrm{km}\mathrm{s}^1`$ and $`V\mathrm{sin}i`$ = 10 $`\mathrm{km}\mathrm{s}^1`$. The latter was chosen so that the code returns the input abundances within $`\pm 0.01`$ dex when abundances are estimated with the input atmosphere model, the S/N ratio being infinite except for rounding errors. Thus all variations in derived abundances can be attributed to errors on $`T_{\mathrm{eff}}`$ or $`\mathrm{log}g`$, even though $`V_{\mathrm{rad}}`$, $`\xi _{\mathrm{micro}}`$ and $`V\mathrm{sin}i`$ were adjusted, as well as the abundances. The input abundances are solar. The results are presented in Fig. 2 and listed in Table 1

An erroneous $`T_{\mathrm{eff}}`$ leads to an error in abundances $``$ 0.1 dex with almost all elements affected in the same way except C, O, and S which react in the opposite direction (see Fig. 2 b and c).

It is interesting to note that an erroneous $`\mathrm{log}g`$ has almost no effect on abundances of iron peak elements (Cr to Zn), and leads to an error $``$ 0.1 dex for the other elements.

When combining both errors, the total error is always $``$ 0.15 dex, as elements that are very sensitive to temperature are not very sensitive to gravity and conversely, at least for the adopted temperature and gravity.

### 5.2 Errors resulting from wrong continuum adjustment

Another source of error comes from bad continuum adjustment. When rotational velocity increases, it becomes harder to find continuum points. The use of echelle spectra is helpful because the whole spectrum from 3900Å to 6820Å is at our disposal. The normalization method is presented in Paper I but can be summarized as follows : the spectrum is divided into 6 parts, with one common order between parts, and each part is normalized separately. Each part has a length of more that 400Å and was chosen so that there are some continuum points near the borders. The continuum is fitted with the continuum function of IRAF whose parameters were defined using slow rotators, mainly the Sun and Vega.

Even when working with high S/N ratios, the main problem for continuum estimation is rotational velocity. It is relatively easy to normalize a spectrum for a star with a $`V\mathrm{sin}i`$ $``$40 $`\mathrm{km}\mathrm{s}^1`$, because lines are not too broad. When rotational velocity increases, lines become broader and it is hard to find lines that are not blended. Blends become so important that continuum points vanish. Selecting wide parts allows to preserve continuum points even for large rotational velocities.

Simulations were done on synthetic spectra having the same parameters as in previous section, but with $`V\mathrm{sin}i`$ = 10, 50, 100, 150 and 200 $`\mathrm{km}\mathrm{s}^1`$. Moreover, in order to simulate as closly as possible real spectra, gaussian noise was added to obtain a reference spectrum with S/N = 200, which corresponds to the minimum observed value. Finally, the spectra were scaled by a factor 0.99 and 1.01, so that the continuum level is respectively too high (by a factor 1.01) and too low (by a factor 0.99). For larger differences, it becomes obvious that the continuum is false. Results of these test are presented in Fig. 3 and in Table 2.

It is interesting to note that for the correct continuum definition, almost all adjusted abundances are well within 0.1 dex of their input values, even for high $`V\mathrm{sin}i`$ (see central column of Fig. 3). Moreover the adjustment with the correct value allows to select “sensitive” elements, such as Al, Si, S, V, Cu, Y, and Nd. These elements have very few and very weak lines, generally blended with much stronger lines.

Another interesting point is the shape of the abundance pattern obtained with a scaled continuum. This shape remains very constant when the rotational velocity increases, but its amplitude increases with this velocity.

Finally, this simulation can be interpreted as yielding the maximum possible error, because all adjusted parts were scaled by one and the same factor. In reality, even for $`V\mathrm{sin}i`$ = 150 $`\mathrm{km}\mathrm{s}^1`$, only the first two parts ($`\lambda <4900`$Å) can be affected by the lack of continuum points, the 5 remaining parts being much easier to normalize. Therefore it is possible to conclude that the influence of the normalization on real spectra will be well below the values listed in Table 2.

## 6 Results

All results are listed in Tables 5 and 6 available at the CDS. The first table lists, for each star, the radial velocity, the microturbulence and the projected rotational velocity. The second one lists the abundances. In 12 cases, there are two spectra per star – hence two lines per star in Tables 5 and 6 – either because the same object was observed once with ELODIE and once with CORALIE (6 cases), or because the star was observed twice with a given instrument (6 more cases).

### 6.1 Abundances

As the sample covers a large temperature domain, it is difficult to summarize the results. For early A-type stars, only a few elements show lines in the spectrum while for F-type stars, it is possible to derive abundances for 28 elements. Table 6, available at CDS, contains the detailed abundances for each star and the number of lines used for the determination.

24 stars of the sample are common with the catalogue of $`\left[\frac{Fe}{H}\right]`$ done by Cayrel de Strobel et al. (cayrel (2001)). Comparisons for $`T_{\mathrm{eff}}`$, $`\mathrm{log}g`$ and $`\left[\frac{Fe}{H}\right]`$ are illustrated in Fig. 4. The agreement is very good. It is intersting to notice that 12 stars come from a single, homogeneous source, Balachandran (balach (1990)). This sample alone shows a slight shift in $`T_{\mathrm{eff}}`$ of less than 100 K and $`\mathrm{log}g`$ values which are in good agreement with ours. $`\left[\frac{Fe}{H}\right]`$ values are in good agreement, except for two cases where the difference amounts to about 0.3 dex.

### 6.2 Microturbulence

Microturbulence is an ad hoc parameter which relates to turbulent motions on scales smaller than the mean free path of a photon. As commonly done in practice, we assumed it is constant with optical depth even though it should be allowed to vary in principle. Our estimates are based on the the four bluer spectral segments only, since there are too few appropriate lines in the three redder segments to constrain the adjustment; the microturbulence value is an average of the four individual fitted values. Because of the large size and homogeneity of our sample, the relation between microturbulence, effective temperature and luminosity is worth examining.

The values are reported in Table 5 available at CDS. Fig. 5 shows the relation between $`T_{\mathrm{eff}}`$ and $`\xi _{\mathrm{micro}}`$ for reliable $`\xi _{\mathrm{micro}}`$ estimates, defined as those for which the standard deviation of the 4 individual fitted values is smaller than 0.5 $`\mathrm{km}\mathrm{s}^1`$. This relation is compared with the average, hand-drawn one presented in Coupry & Burkhart (coupry92 (1992)), which is the continuous curve in Fig. 5. The agreement is good for cool stars but, beyond $`T_{\mathrm{eff}}`$ $`=7400`$ K, the relation of Coupry & Burkhart becomes an upper limit to our values.

The relation of Nissen (N81 (1981)), which is linear in both $`T_{\mathrm{eff}}`$ and $`\mathrm{log}g`$, is represented by the three straight lines in its domain of validity. The agreement seems not too bad for low luminosity stars, but our data seem to indicate that all stars share the same $`\xi _{\mathrm{micro}}`$ at low temperature, whatever their $`\mathrm{log}g`$, while evolved stars separate from the unevolved ones very rapidly as $`T_{\mathrm{eff}}`$ increases. The $`T_{\mathrm{eff}}`$ dependence of $`\xi _{\mathrm{micro}}`$ is much steeper in our case than in Nissen’s, and the slope seems to increase with decreasing $`\mathrm{log}g`$. Qualitatively, however, we confirm that evolved stars have higher $`\xi _{\mathrm{micro}}`$ in this range of $`T_{\mathrm{eff}}`$. This is also in qualitative agreement with the results of Gray et al. (gray (2001)) for stars cooler than A8.

Strangely enough, the trend reverses for $`T_{\mathrm{eff}}`$ between 8000 and 8800 K, evolved stars having significantly lower $`\xi _{\mathrm{micro}}`$ on average than unevolved ones. The latter trend is largely due to the group of stars with $`T_{\mathrm{eff}}`$ larger than 8000 K and $`\xi _{\mathrm{micro}}`$ larger than 2.4, which correspond to stars with mass of 1.8 to 2.0 M that present iron and nickel overabundances (see Sect. 7.3 for details). These stars are very close to the Zero Age Main Sequence (ZAMS) in the HR diagram. This behaviour was already noted by Gray et al. (gray (2001)). Whether $`\xi _{\mathrm{micro}}`$ depends on $`[M/H]`$ remains an open question, and more high quality data would be needed to better understand what is happening in this diagram. Let us notice that the observed overabundance of iron cannot be the result of a bias in the estimate of $`\xi _{\mathrm{micro}}`$, since only an underestimate of it would lead to an overestimate of iron abundance; what we see is, on the contrary, rather high $`\xi _{\mathrm{micro}}`$ values for these metal rich stars.

It is difficult to draw any conclusion about a possible relation between $`\xi _{\mathrm{micro}}`$ and $`\mathrm{log}g`$ for the stars hotter than 8800 K because their number is too small. In addition, the $`\xi _{\mathrm{micro}}`$ estimates are less reliable than for cool stars because of the much smaller number of lines.

## 7 Discussion

### 7.1 Rotational velocities

#### 7.1.1 Comparison with other works

Fig. 6 illustrates the comparison of our determinations of $`V\mathrm{sin}i`$ with works of Uesugi & Fukuda (uesugi (1982)), Abt & Morell (abt (1995)), and Royer et al. (royer (2002)). Our values are slightly larger than those of Abt & Morell, especially for $`V\mathrm{sin}i`$ larger than 100 $`\mathrm{km}\mathrm{s}^1`$. However, the trend is in the other direction when comparing with Uesugi & Fukuda. Finally, comparison with Royer et al. is excellent if we neglect one point ($``$150 $`\mathrm{km}\mathrm{s}^1`$ for our estimate). This point corresponds to HD 12311. This star is still classified as a dwarf in Simbad even though it has the atmospheric parameters of a giant. This was already noted by Hauck (hauck (1986)) using Geneva photometry, and Gray & Garrison (1989b ) give a spectral type of F0 III-IVn. Royer et al. (royer (2002)) use Fourier transforms of several line profiles to estimate $`V\mathrm{sin}i`$. This technique is very efficient for stars with well separated lines. In cooler stars, it becomes harder to find isolated lines to analyze. HD 12311 is precisely one of the coolest stars of Royer et al. (royer (2002)). In this case, all lines of their list except one were unusable. Therefore the difference can be explained by a poor precision of their determination for this star. In this comparison, there are less points because Royer et al. focus on A type stars.

We conclude that our determination of $`V\mathrm{sin}i`$ are probably even better than those of Royer et al., because we make use of all lines present in the spectra, whether they are blended or not. On the other hand, we neglect macroturbulence, so our $`V\mathrm{sin}i`$ values are slightly biased towards large values in very slow rotators ($`V\mathrm{sin}i`$ $`<810`$ $`\mathrm{km}\mathrm{s}^1`$).

#### 7.1.2 Systematic abundance error due to rotation

When dealing with large rotational velocities, one can easily underestimate the continuum level. That leads to underestimated abundances because lines seem shallower.

The apparent abundances of three elements are strongly correlated with $`V\mathrm{sin}i`$. These elements are Sr, Sc, and Na. This is illustrated in Fig. 7, where iron was also included to show the behavior of an uncorrelated element. The correlations involving Sr, Sc and Na are spurious, however, for the following reasons :

##### Strontium:

This element shows the most significant correlation (correlation coefficient $`r=0.61`$). The problem arises when only one line, Sr ii $`\lambda `$4215.519, is visible. This happens mainly in the hotter stars of the sample. The fact is that the continuum is slightly underestimated in this region. When rotational velocity increases, this small error leads to large underestimates of the abundances. As there is no other line of Sr for those stars, this error cannot be corrected. In Fig. 7, almost all stars with $`V\mathrm{sin}i`$ $`>150`$ $`\mathrm{km}\mathrm{s}^1`$ are of spectral type A, which reflects the well-known drop of rotational velocity for spectral types F and later. It is for these relatively hot stars that only one Sr ii line is used, and in this case, the result is much poorer than in the simulation illustrated by Fig. 3, which involves an F star where two Sr lines are used.

This is the only negative correlation for all adjusted elements. It is very interesting to illustrate the accuracy of the continuum adjustment, and the stability of this method. Indeed, an underestimate of the continuum leads to a negative correlation. The fact that only one element shows such a behaviour is a good justification for our choice of continuum adjustment.

##### Scandium:

This element yields the largest positive correlation ($`r=0.5`$). It should be remembered that during the analysis, spectra are cut into 7 parts, each part is adjusted, and finally the different estimates are put together (see Paper I for details). In the case of Sc, there is a large scatter between abundances of different spectral ranges. Moreover, since Sc is ionized, the adjustment is very sensitive to gravity (see Sect. 5.1). However, there is no correlation between $`\mathrm{log}g`$ and $`V\mathrm{sin}i`$. Unfortunately, the reason of this behaviour was not found, but it is probable that there are several causes adding up.

##### Sodium:

The last significant correlation ($`r=0.43`$) is found for Na. For this element, the reason was found to be the telluric lines. The Na I doublet ($`\lambda `$ 5889 and $`\lambda `$ 5895) lies within telluric lines. Even if the adjustment is not very sensitive to these lines, a small error can lead to a large abundance overestimate, especially for fast rotators, because these lines are broad and saturated.

There is no significant correlation between other elements and $`V\mathrm{sin}i`$.

#### 7.1.3 Correlation between rotation and true abundances

The decrease of $`\left[\frac{Fe}{H}\right]`$ with increasing $`V\mathrm{sin}i`$ found by Takeda & Sadakane (takeda97 (1997)) and confirmed by Varenne & Monier (varenne99 (1999)) is firmly denied by our results. Fig. 8 shows clearly that there is no dependence of iron or nickel abundance with $`V\mathrm{sin}i`$. In our case, there are 60 stars with $`V\mathrm{sin}i`$ $``$ 100 $`\mathrm{km}\mathrm{s}^1`$ while there was only 10 in Varenne & Monier and 5 in Takeda and Sadakane. It is important to notice that the anticorrelation was found for stars of the Hyades, but, a priori, there is no reason why only these stars should be affected. The result of Takeda et al. (takeda97 (1997)) is discussed in more details in Sect. 7.3

### 7.2 Metallic giants

Fig. 9 shows the position of all giant stars of our sample in the HR diagram. There are 35 giants, out of which 18 are metallic. Their characteristics are reported in Tables 3 and 4.

Abundances of observed giants are illustrated in Fig. 10. There are no very clear-cut peculiarities such as the underabundance of Ca and/or Sc in Am stars. However, abundances of some elements (Al, Ca, Ti, Fe, Ba for example) are clearly different.

To ensure that the abundances are really different, the Kolmogorov-Smirnov test was used. This test gives the probability that two distributions are drawn from the same parent population.

It was already shown in North et al. (north98 (1998)) that the $`\mathrm{log}g`$ distribution of metallic and normal giants are different with a probability of 99.99%. However, the $`T_{\mathrm{eff}}`$ distributions are comparable. The distributions of $`V\mathrm{sin}i`$ and $`\xi _{\mathrm{micro}}`$ estimated by our method are also comparable.

Fig. 11 shows the cumulative distribution of abundances for each of the 6 elements that have the largest differences. They are Ca, Ti, Mn, Fe, Ni and Ba. Al and Cr also have different distributions. Sc is not added, even though the corresponding distributions are clearly different, because of its correlation with $`V\mathrm{sin}i`$.

Thus, we conclude that metallic giants have abundances of Al, Ca, Ti, Cr, Mn, Fe, Ni and Ba that are significantly larger than those of normal giants.

### 7.3 Abundance variation with evolution

Recently, the Montréal group has published a new study for normal stars. Models including all effects of atomic diffusion and radiative accelerations have been presented for stars of different masses in Turcotte et al. (turcotte98 (1998)), Richer et al. (richer00 (2000)) and Richard et al. (richard01 (2001)).

The models of Turcotte et al. predict the evolution of the abundances for an F star of a given mass consistently with the internal structure. They are computed in the diffusion framework with only the mixing length as free-parameter. They include the effect of 28 chemical elements (Z $``$ 28) but no macroscopic mixing. Atomic diffusion plays an important role for stars more massive than 1.3 M. We agree with Varenne & Monier (varenne99 (1999)) that the theoretical prediction of Turcotte et al. (turcotte98 (1998)) for F stars are not confirmed by observations. In our case, abundance of C, O, Si, Ca, Fe, and Ni do not show any variations with effective temperature (see Fig. 12).

Richard et al. (richard01 (2001)) present stellar models that are evolved from the pre-main sequence to the giant branch for stars of 1.3 to 4.0 M and metallicities of Z=0.01 to 0.03. In order to compare our observations with these models, our sample was split into 7 mass groups with masses below 1.4 M, 1.4 to 1.6 M, 1.6 to 1.8 M, …, 2.4 to 2.6 M and above 2.6 M. Fig. 13 shows the behaviour of Ca, Fe and Ni for stars of 1.6 to 1.8 M and 1.8 to 2.0 M with respect to $`\mathrm{log}g`$. It can be directly compared with Fig. 10 of Richard et al. (richard01 (2001)). They predict a large abundance variation of at least 1 dex with decreasing $`\mathrm{log}g`$ (or equivalently, with time ). In our case, there is no variation with $`\mathrm{log}g`$ except for stars of 1.8 to 2.0 M, and even in this case, the amplitude is much smaller than the theoretical predictions. In order to confirm the statistical significance of the correlation, $`2\times 2`$ contingency tables have been computed, segregating stars with gravities larger or lower than $`\mathrm{log}g`$ $`=4.0`$ and with \[Fe/H\] or \[Ni/H\] larger or lower than $`0.0`$. The corresponding $`\chi ^2`$ were computed taking Yates’ correction into account (Crow et al. CDM60 (1960)) and the results were $`8.98`$ and $`10.94`$ for \[Fe/H\] and \[Ni/H\] respectively, which comfortably exceeds the $`7.88`$ value corresponding to a 0.5% probability that abundances and $`\mathrm{log}g`$ are unrelated. Thus, taken at face value, the correlation is statistically significant. If this variation is physically true, it suggests that additional hydrodynamical processes are at play in these stars, which lower the effects of microscopic diffusion. However, before concluding to the physical reality of this correlation, it is necessary to test how far this behaviour can be explained by a correlation of abundances with $`V\mathrm{sin}i`$ or by NLTE effect.

Takeda & Sadakane (takeda97 (1997)) also found a variation of iron abundance with $`\mathrm{log}g`$, which was qualitatively the same, for Hyades’ members. They explain this variation as a superficial phenomenon due to the existence of a mutual relationship between $`\mathrm{log}g`$ and $`V\mathrm{sin}i`$. They found an anticorrelation between $`V\mathrm{sin}i`$ and $`\mathrm{log}g`$, in the sense that high-gravity stars tend to be slow rotators. They consider that the cause of this relation is rotation, by which the apparent and effective values of $`\mathrm{log}g`$ are affected. In their case, $`\mathrm{log}g`$ is estimated from Strömgren photometry, and rotation affects this determination: it appears that $`\mathrm{log}g`$ values progressively decrease with increasing $`V\mathrm{sin}i`$.

However, in our case, $`\mathrm{log}g`$ values were estimated from Hipparcos data, and they are not correlated with $`V\mathrm{sin}i`$ (see Fig. 14). The effect of rotation is to increase the $`B2V1`$ index. Maeder (maeder (1971)) has shown that, to a good approximation, this index is essentially sensitive to the product $`V\mathrm{sin}i`$. For a star of 1.4 M (2.0 M respectively), a $`V\mathrm{sin}i`$ of 200 $`\mathrm{km}\mathrm{s}^1`$ leads to a decrease in log($`T_{\mathrm{eff}}`$) of 0.016 (0.010 respectively). The effect on $`M_V`$ is of a few hundredths of a magnitude ($`0.1`$ mag respectively) (Maeder & Peytremann, maeder70 (1970)). This affects slightly the mass estimate when interpolating through evolutionary tracks of non-rotating models, but the error cannot exceed 0.1 M (5% at 2 M), propagating into no more than $`0.02`$ dex in $`\mathrm{log}g`$. In fact our estimates of $`\mathrm{log}g`$ are close to the effective gravity, i.e. the sum of gravitational and centrifugal accelerations average over the visible disk; we verified that they agree quite well with photometric estimates computed from the $`uvby\beta `$ colours and the calibration by Moon & Dworetsky (MD85 (1985)). The interesting thing with effective $`\mathrm{log}g`$ is that for $`V\mathrm{sin}i`$ $`200`$ $`\mathrm{km}\mathrm{s}^1`$, it differs by no more, on the equator, than 0.1 dex from a non-rotating 2 M star near the ZAMS, while the difference is at least twice as large at the turn-off. Thus the observed $`\mathrm{log}g`$ axis (e.g. in Fig. 13) is slightly stretched for fast rotators, relative to the $`\mathrm{log}g`$ non-rotating stars would have in the same mass range. But since slow rotators are not affected, the net effect is just to blur very slightly the relation metallicity vs. $`\mathrm{log}g`$ shown in Fig. 13. Therefore, this cannot explain the behaviour of $`\left[\frac{Fe}{H}\right]`$ with $`\mathrm{log}g`$. Moreover, in Takeda & Sadakane (takeda97 (1997)), $`\mathrm{log}g`$ varies on a range that is no more than 0.2 dex; under the assumption that this variation is due entirely to a bias in the $`\mathrm{log}g`$ estimate, it can hardly explain a 0.4 dex variation of $`\left[\frac{Fe}{H}\right]`$ for the single line of neutral iron they used.

Given the variation of $`T_{\mathrm{eff}}`$ and $`\mathrm{log}g`$ for these stars, NLTE effects might play an important role. Rentzsch-Holm (rentzsch (1996)) has estimated NLTE effects on iron abundances for A-type stars. NLTE correction increase with $`T_{\mathrm{eff}}`$ going from 7 000 to 10 000 K (see Fig. 4 of Rentzsch-Holm (rentzsch (1996))). On the other hand, at given effective temperature, they increase with decreasing gravity. Since both $`T_{\mathrm{eff}}`$ and $`\mathrm{log}g`$ decrease with time when a star is evolving on the main sequence, this tends to flatten the NLTE correction. For our stars of 1.8 to 2.0 M, there is a small NLTE correction of 0.1 to 0.15 dex. However, this correction is almost constant (variation of 0.05 dex) for the $`T_{\mathrm{eff}}`$ and $`\mathrm{log}g`$ considered333Note that this is valid only if $`T_{\mathrm{eff}}`$ and $`\mathrm{log}g`$ are correlated. The Pearson product-moment correlation coefficient is 0.972 (0.968) for stars of 1.8 to 2.0 M (1.6 to 1.8 M respectively), and in the same order for the other mass interval.. Moreover, with our method, we do not have separate abundance values for Fe i and Fe ii. Let us recall that Fe ii lines are almost not affected by NLTE effects. When $`T_{\mathrm{eff}}`$ increases, the ratio of line numbers of Fe i and Fe ii decreases. Therefore our determination becomes less sensitive to NLTE effects. This also tends to flatten the variation of NLTE correction with $`T_{\mathrm{eff}}`$. Therefore, we conclude that, for stars of 1.8 to 2.0 M, only a small NLTE correction of about 0.1 dex, that is roughly constant with $`T_{\mathrm{eff}}`$, has to be taken into account when iron abundance is estimated with our method. To summarize, all points in the $`\left[\frac{Fe}{H}\right]`$ vs $`\mathrm{log}g`$ diagram for stars of 1.8 - 2.0 M in Fig. 13 should be shifted upward by about 0.1 to obtain a more realistic picture.

Thévenin & Idiart (thevenin99 (1999)) discuss the problem of $`\mathrm{log}g`$ determination when dealing with NLTE correction. Their work focuses on cooler stars, but the problem of $`\mathrm{log}g`$ determination is comparable. When the NLTE correction is applied to LTE abundances, abundances of Fe i and Fe ii are no longer the same (assuming that $`\mathrm{log}g`$ had been adjusted to satisfy the LTE ionization equilibrium). Therefore, $`\mathrm{log}g`$ has to be corrected until ionization equilibrium is satisfied after NLTE correction. Nissen et al. (nissen97 (1997)) have shown that there exists a discrepancy between spectroscopic $`\mathrm{log}g`$ (adjusted assuming LTE) and those deduced from Hipparcos parallaxes. Thévenin & Idiart (thevenin99 (1999)) have solved this discrepancy by adjusting $`\mathrm{log}g`$ after NLTE correction; the latter are in close agreement with those deduced from Hipparcos (our $`\mathrm{log}g`$ are estimated in this way). However, in order to apply NLTE corrections precisely, it is necessary to have separate abundance estimates of Fe i and Fe ii, which is not our case. Roughly, our determination can be seen as an average of both abundances weighted by the number of lines of each ionization stage.

To summarize, NLTE effects on iron is no more than 0.1 dex and remains constant to better than 0.1 dex, for stars of 1.8 to 2.0 M when abundance are estimated with our method.

In view of the above discussion, it seems that the correlation between $`\left[\frac{Fe}{H}\right]`$ and $`\mathrm{log}g`$ for 1.8-2.0 M is indeed physical, and due neither to systematic errors nor to selection effects. It agrees qualitatively with the prediction of Richard et al. (richard01 (2001)), but it is strange that only the stars of 1.8 to 2.0 M follow their models. It is true that more massive stars are less numerous, so their statistics is less good. For less massive stars, models predict a variation that occurs on a shorter $`\mathrm{log}g`$ range. It is possible that this variation is too quick to be observed, but then we would expect to see at least an increase of the scatter for large $`\mathrm{log}g`$ values, while this is not the case.

The comparison with the Montréal models shows that it is necessary to add hydrodynamical processes in order to reproduce observations, like mass loss and meridional circulation, as already suggested by the Montréal group. Maybe turbulence will also have to be modelled in a slightly different way yet to be defined.

We also note that the metal-rich 1.8 to 2.0 M stars lie just around the blue border of the instability strip or outside of it: the coolest of these (HD 144197) has $`T_{\mathrm{eff}}`$ $`=7691`$ K and the next cooler one (HD 148367) has $`T_{\mathrm{eff}}`$ $`=7962`$ K, while the blue limit of the instability strip is defined by $`B2V1=0.01`$ on the ZAMS (Jasniewicz jas84 (1984)) which correspond to $`T_{\mathrm{eff}}`$ $`=8016`$ K according to the calibration of Hauck & Künzli (HK96 (1996)). The hottest metal-rich star (HD 195943) has $`T_{\mathrm{eff}}`$ $`=9057`$ K. Almost all stars with $`[Fe/H]<0.0`$ are cooler and lie, therefore, either inside the instability strip or to the red side of it. This might be more than a coincidence, and is reminiscent of the well known exclusion between $`\delta `$ Scuti type pulsations and the Am phenomenon. The difference is that while most Am stars do lie inside the instability strip even though they do not pulsate, the metal-rich 1.8 to 2.0 M stars are found either outside or on the verge of it.

## 8 Conclusion

The method of analysis presented in Paper I was used to derive detailed abundances of 140 A and F-type stars with $`V\mathrm{sin}i`$ up to about 200 $`\mathrm{km}\mathrm{s}^1`$. $`V\mathrm{sin}i`$ derived with this method are consistent with those in the literature. Only Cr, Sc, and Na are very sensitive to rotational velocity in the sense that their abundance determination is altered when $`V\mathrm{sin}i`$ $``$ 100 $`\mathrm{km}\mathrm{s}^1`$, for various reasons.

This study allows us to specify more exactly the nature of metallic A-F giants. They are more evolved than normal giants and they show overabundances of at least 8 elements with respect to normal giants. These elements are Al, Ca, Ti, Cr, Mn, Fe, Ni and Ba.

Our results do not confirm, in general, the predictions of the Montréal models. However, stars of 1.8 to 2.0 M show a clear correlation of iron abundances with $`\mathrm{log}g`$ which is reminiscent of theoretical predictions by Richard et al. (richard01 (2001)). Neither rotational velocity nor NLTE effects can explain this behaviour, which appears, therefore, to be real. The observed $`\left[\frac{Fe}{H}\right]`$ and $`\left[\frac{Ni}{H}\right]`$ variations qualitatively agree with those predicted by the Montréal models, but are much weaker. The smaller amplitude indicates that their models need more hydrodynamical processes (meridional circulation, mass loss, …) and/or a better tuning of convection. Our large homogeneous set of abundances hopefully can help this work.

###### Acknowledgements.

We express our warm thanks to P. Bartholdi and C. Briner for their help in using the computer network of the Geneva observatory and for the instructive discussions. Constructive criticism by the referee, Dr. Richard Gray, is gratefully acknowledged. This work has been partly supported by the Swiss National Science Foundation.