Hints on Galaxy Formation from the Metal Content of Galaxy Cluster and Other Fossil Evidence

Abstract: Clusters of galaxies allow a direct estimate of the metallicity and metal production yield on the largest scale so far. It is argued that cluster metallicity ($\sim 1/3$ solar) should be taken as representative of the low-$z$ universe as a whole. There is now compelling evidence that the bulk of stars not only in cluster ellipticals but also in field ellipticals and bulges formed at high redshifts ($z\gsim 3$). Since such stars account for at least $\sim 30%$ of the baryons now locked into stars, it is argued that at least 30% of stars and metals formed before $z\simeq 3$, and correspondingly the metallicity of the universe at $z=3$ is predicted to be $\sim 1/10$ solar. This requires the cosmic star formation rate to run at least flat from $z\sim 1$ to $\sim 5$, which appears to agree with the most recent derect determinations of the star formation rate in Lyman-break galaxies.

\section{Introduction}

The observation of forming galaxies at high redshift is certainly the most
direct way to look at the {\it Birth of Galaxies}, and a great
observational effort is currently being made in this direction. Yet,
high redshift galaxies are very faint, and at least for now only few
of their global properties can be measured. Nearby and moderate
redshift galaxies can instead be studied in far greater detail, and
their {\it fossil evidence} can provide a view of galaxy
formation and evolution that is fully complementary to that given
by high redshift observations. By fossil evidence I mean those
observables that do not refer to ongoing, {\it active} star formation,
and which are instead the result of the integrated past star formation
history.  This evidence includes the global metallicity of the
universe, and the distribution of stellar ages in $z\simeq 0$ galaxies
as well as in high redshift, but passively evolving galaxies.

The paper is organized as follows. Section 2 presents the current
evidence for the chemical composition of local clusters of galaxies
at low redshift, separately for the intracluster medium (ICM) and for
cluster galaxies.
In Section 3 the current evidence  is presented for
the bulk of stars in galactic spheroids (i.e. ellipticals and bulges
alike) being very old, formed at high redshift, no matter whether they
reside in rich clusters or in the low density environment we usually
refer to as the {\it field}. In Section 4 these evidences are used to
set constraints on the main epoch of star formation, metal production,
and on the metallicity of the high-$z$ universe.

\section{Clusters as Archives of the Past Star and Metal Production}
Individual galaxies are not good examples of the {\it closed box}
model of chemical evolution: enriched matter certainly flows out of
them,
and quasi-pristine material may occasionally infall and be accreted by them.
Theoretical simulations predict instead that the baryon fraction of rich
clusters
cannot change appreciably in the course of their evolution [43].
We can then expect within a cluster
to find confined in the same place
all the dark matter, all the baryons, all the galaxies, and all the
metals that have participated in the play. Clusters are then good {\it
archives} of their past star formation and metal production history,
an ideal laboratory for the study of the fossil evidence.

Metals in clusters are partly spread through their ICM,
partly locked into galaxies and stars. The mass
of metals in the
ISM of galaxies is instead negligible compared to that in the two
other components. ICM abundances can be obtained from X-ray
observations, while optical observations combined to population
synthesis models provide estimates for the metallicity of the stellar
component of galaxies.

\subsection{Iron and $\alpha$-Elements in the Intracluster Medium}

Among abundant elements, iron has the most accurately  measured ICM abundance.
The
so-called iron-K emission complex at $\sim 7$ keV is prominent in the
X-ray spectrum of clusters, and allows fairly precise measurements.
Fig. 1 shows the
iron abundance of clusters and groups as a function of ICM
temperature, while
Fig. 2 shows the iron-mass-to-ligh-ratio (Fe$M/L$) of the ICM [30]. This is
defined as the ratio $\hbox{$M_{\rm Fe}^{\rm ICM}$}/\hbox{$L_{\rm B}$}$ of the total iron mass in
the ICM over the total $B$-band luminosity of the galaxies in the cluster.

\begin{figure}
\vskip-5.5truecm
\centerline{\psfig{file=renzini-fig1.eps,width=15cm,angle=0,height=14cm}}

\vskip-1truecm
\caption{\piedi A compilation of the iron abundance in the ICM as a
function of ICM temperature
for a sample of clusters and groups [30], including several clusters at
moderately high redshift with $<z>\simeq 0.35$, represented by small
filled circles.}
\end{figure}

The drop of the derived Fe$M/L$ \ in poor clusters and groups (i.e. for
 $kT\ifmmode{\mathrel{\mathpalette\@versim<} 2$ keV) can be traced back to a drop in both factors
 entering in its definition, i.e., in the iron abundance (cf. Fig. 1)
{\it and} in
 the ICM mass to light ratio [30].  We don't know whether this drop is a
 real effect.  Groups may not behave as closed boxes, and may be
 subject to baryon and metal losses due to the feedback effect of star
 formation in individual galaxies driving galactic winds [34] [9].
In addition, there may be a
 diagnostic problem, since in groups iron is derived
 from iron-L transitions, which involve atomic physics which is
 far more complicated than that of the iron-K complex [2][30].
For this reason I will not further
 discuss poor clusters, i.e. those cooler than $kT\ifmmode{\mathrel{\mathpalette\@versim<} 2$ keV.

Fig. 1 and 2 show that both the iron abundance and the
Fe$M/L$ \ in rich clusters ($kT\ifmmode{\mathrel{\mathpalette\@versim>} 2$ keV) are independent of
cluster temperature, hence of cluster richness and
optical luminosity. For these clusters one has
$\hbox{$Z^{\rm Fe}_{\rm ICM}$} =0.3\pm 0.1$ solar, and $\hbox{$M_{\rm Fe}^{\rm ICM}$}/\hbox{$L_{\rm B}$} = (0.02\pm 0.01)$ for
$\hbox{$H_\circ$}=50$.
The most straightforward interpretation is that clusters did not lose
iron (hence baryons), nor selectively acquired pristine baryonic material, and
that the conversion of baryonic gas into stars and galaxies has
proceeded with the same efficiency and stellar IMF in all clusters [30].
Otherwise, we should observe cluster to cluster variations of the iron
abundance and of the Fe$M/L$.
The theoretical prediction of the constancy of the baryon fraction in clusters
[43] is nicely supported by these evidences.

\begin{figure}
\vskip-5truecm
\centerline{\psfig{file=renzini-fig2.eps,width=15cm,angle=0,height=14cm}}
\vskip-1truecm
\caption{\piedi The iron mass to light ratio  of  the  ICM of clusters
and groups (for $\hbox{$H_\circ$}=50$) as a function of the ICM temperature [30].}
\end{figure}

X-ray observations have also allowed to measure the abundance of
other elements in the ICM, especially that of the $\alpha$-elements such as
O, Ne, Mg, and Si, with the {\it ASCA} X-Ray telescope having
superseded all
previous attempt in this respect. A fairly high $\alpha$-element
enhancement, with
$<\![\alpha$/Fe]$>\simeq +0.4$, was initially reported [23]. Later  this
estimate has been revised down to
$<\![\alpha$/Fe]$>\simeq +0.2$  [24],
using a global average among observed clusters, and averaging
O, Ne, Mg, and Si for the $\alpha$-elements [30].
In this way a more robust estimate of the global [$\alpha$/Fe] ratio
is obtained, given the errors affecting the abundances
of individual elements in individual galaxies. This may still suggest a modest
$\alpha$-element enhancement, with the ICM enrichment being dominated
by SNII products [23].

However, it has been pointed out that this
small apparent $\alpha$-element enhancement in the ICM comes from
having assumed the reference solar iron abundance from the
``photospheric''
model atmosphere analysis [15].
The ``meteoritic'' iron abundance is instead
now generally adopted, and since this is $\sim 0.16$
dex  lower than the photospheric value, one can conclude that
there is virtually no $\alpha$-element enhancement at all in
the ICM (formally $<\![\alpha$/Fe]$>\simeq +0.04\;\pm\sim 0.2$, [15], [30]).
Clusters of galaxies
are therefore nearly {\it solar} as far as the
elemental ratios are concerned, which argues for stellar
nucleosynthesis having proceeded in quite the same way in the solar
neighborhood as well as at the galaxy cluster scale.
This implies a similar ratio of the number of Type Ia to Type II SNs,
as well as
a similar IMF [30]. This result speaks for the universality of the
star formation process (IMF, binary fraction, etc.), and
may help limiting the number of free
parameters to play with.

\subsection{The Iron Content of Galaxies and the Iron Share}

The metal abundance of the stellar component of cluster galaxies can
only be inferred from integrated spectra coupled to synthetic stellar
populations. Much of the stellar mass in clusters is confined to
passively evolving spheroids (ellipticals and bulges) for which the
iron abundance may range from $\sim 1/3$ solar to a few times solar.
Hence, the average iron abundance cannot be much different from solar,
even when taking into account the presence of radial gradients
[2], $\alpha$-element enhancements [8], and the luminosity bias [14].

The  global iron
abundance of a whole cluster is  therefore given by:
$$Z_{\rm CL}^{\rm Fe}={\hbox{$Z^{\rm Fe}_{\rm ICM}$}\hbox{$M_{\rm ICM}$} + \hbox{$Z^{\rm Fe}_*$} M_* \over \hbox{$M_{\rm ICM}$} + M_*}=
      {5.5\hbox{$Z^{\rm Fe}_{\rm ICM}$} h^{-5/2} + \hbox{$Z^{\rm Fe}_*$} h^{-1}\over 5.5h^{-5/2} +
      h^{-1}},\eqno(1)$$
where $\hbox{$Z^{\rm Fe}_*$}$ is the average abundance of stars in galaxies and $ M_*$
is the mass in stars. For the second equality
I have assumed as prototypical the Coma cluster values
adopted by White et al. [43]: $\hbox{$M_{\rm ICM}$}\simeq 5.5\times
10^{13}h^{-5/2}\hbox{$M_\odot$}$ and $M_*\simeq 10^{13}h^{-1}\hbox{$M_\odot$}$.
With $\hbox{$Z^{\rm Fe}_{\rm ICM}$}=0.3$ solar and $\hbox{$Z^{\rm Fe}_*$}=1$ solar,
equation (1) gives a global cluster abundance of 0.34,
0.37, and 0.41 times solar, respectively for $h=0.5$, 0.75, and 1.
Under the same assumptions, the ratio of the iron mass in the ICM to
      the iron mass locked into stars is:
$${\hbox{$Z^{\rm Fe}_{\rm ICM}$}\hbox{$M_{\rm ICM}$}\over\hbox{$Z^{\rm Fe}_*$} M_*}\simeq 1.65 h^{-3/2},\eqno(2)$$ or 4.6,
2.5, and 1.65, respectively for $h=0.5$, 0.75, and 1. Note that with
the adopted values for the quantities in equation (2) most of the
cluster iron resides
in the ICM, rather than being locked into stars, especially for low
values of $\hbox{$H_\circ$}$. These estimates could be somewhat decreased if
clusters contain a sizable population of stars not bound to censed
individual galaxies, if the average iron abundance in stars is
supersolar (luminosity-weighted determinations underestimate true
abundances [14], or if the galaxy $M_*/L$ ratio is higher
than adopted here, i.e., $<\!M_*/\hbox{$L_{\rm B}$}\!>=6.4h$ [43].
However, the bottom line is that there is at least as much metal mass inside
cluster galaxies, as there is out of them in the ICM. This must be
taken as
a strong constraint by models of the chemical evolution of galaxies:
clearly galaxies do not evolve as a closed box, and outflows must play a
leading role.

With the adopted masses and iron abundances for the two baryonic
components  one can also evaluate the total
cluster Fe$M/L$:
$${\hbox{$M_{\rm Fe}^{\rm ICM}$} +\hbox{$M_{\rm Fe}^*$}\over\hbox{$L_{\rm B}$}}\simeq 1.3\times 10^{-2}(1.65\, h^{-1/2}+h)
\; (\hbox{$M_\odot$}/\hbox{$L_\odot$}),\eqno(3)$$
or Fe$M/L$=0.037 or 0.034 $\hbox{$M_\odot$}/\hbox{$L_\odot$}$, respectively for $h=0.5$ and
1. The total
Fe$M/L$ \ is therefore fairly insensitive to the adopted distance scale.
Simple calculations
show that to reproduce this value one needs either
a fairly flat IMF ($x\simeq 0.9$) if all iron is attributed to SNII's,
or a major contribution from SNIa's, if one adopts a Salpeter IMF
($x=1.35$) [34]. The former option
dictates a substantial $\alpha$-element enhancement, similar to the
values observed in the Galactic halo ([$\alpha$/Fe]$\simeq +0.5$).
The latter option instead predicts near solar proportions for the
cluster as a whole, and requires the SNIa rate to have been much
higher in the past [34]. The evidence presented in Section 2.1  favors
the second option.

From the near solar proportions of cluster
abundances one obtains the total metal mass to light ratio of a typical
cluster as $M_{\rm Z}/\hbox{$L_{\rm B}$}\simeq 10\times M_{\rm Fe}/\hbox{$L_{\rm B}$}\simeq  0.3\pm 0.1\;
(\hbox{$M_\odot$}/\hbox{$L_\odot$})$.
It is worth noting that this is an interesting, fully empirical
estimate of the metal
yield of stellar populations. Following
Tinsley [41], the metal yield is usually defined per unit mass of
stars, a quantity which theoretical counterpart depends on the poorly known
low mass end of the
IMF. The estimate above gives instead the yield per unit luminosity of
present day cluster galaxies, a quantity that depends on the IMF only
for $M\ifmmode{\mathrel{\mathpalette\@versim>}\hbox{$M_\odot$}$.  Theoretical mass-related yields have been recently
estimated by Thomas, Greggio, \& Bender
[40] based on massive star models and supernova explosion [39], [44].
These yields can be purged
from their mass dependence, and transformed into luminosity-related
yields. For this purpose I assume an age of 15 Gyr for the bulk of
stars in clusters (cf. Section 4), and use the proper luminosity-IMF
normalization [29] [31]: i.e. $\psi(M)=AM^{-(1+x)}$ for the IMF with
$A\simeq 3.0\hbox{$L_{\rm B}$}$. Thus, theoretical yields turn out to be $M_{\rm
Z}/\hbox{$L_{\rm B}$}=0.08$, 0.24, and $0.33\; \hbox{$M_\odot$}/\hbox{$L_\odot$}$, respectively for
$x=1.7$, 1.35, and 1.00, which compares to $M_{\rm Z}/\hbox{$L_{\rm B}$}\simeq 0.3\pm
0.1 \hbox{$M_\odot$}/\hbox{$L_\odot$}$ for the empirical cluster value. One can conclude
that current stellar yields do not require a very flat IMF to account
for the cluster metals.

\subsection{Clusters as Fair Samples of the Local Universe}

To what extent the cluster global metallicity, and
the ICM to galaxies iron share are representative of the low-$z$
universe as a whole? For example, Madau et al. [21] adopt $\hbox{$H_\circ$}=50$, a stellar
mass density parameter $\Omega_*=0.0036$, and a baryon mass density
parameter $\Omega_{\rm b}=0.05$. With these values the fraction of
baryons that have been locked into stars is $\sim 7\%$.
This compares to  $\sim
1/(1+5.5h^{-3/2})$ in clusters, or $\sim 6\%$ and $\sim 10\%$,
respectively for $h=0.5$ and 0.75.
Therefore, it appears that the efficiency of baryon conversion into
galaxies and stars adopted  by Madau et
al. [21] is nearly the same as that observed in clusters, which
supports the notion of clusters being representative of the low-$z$
universe $(\Omega_*/\Omega_{\rm b})_\circ$.

For lack of direct evidence on the metallicity of the local
 intergalactic medium (IGM),  Madau et
al. [21] assumed all the metals in the low redshift universe to be
 locked into stars (galaxies),  with
a negligible metal content for the IGM, which however
comprises the vast majority of all the baryons. With these
assumptions, and adopting the average metallicity of all stars to be
 solar,  the metallicity of the present day universe is $\sim
1\times 0.0036/0.05= 0.07$ solar, or $\sim 5$ times lower than the measured
value in clusters of galaxies.

Is this difference real, or does it just reflect the assumption of a
zero metallicity IGM? There is no obvious reason why the metal yield
of stellar populations should be $\sim 5$ times lower in field
galaxies compared to cluster galaxies.
On the contrary, assuming the
difference to be real would force us to accept
a drastic difference between the behavior of
 galaxies and stellar populations in  clusters and in the field, which
 would require quite contrived explanations [30]. Much simpler, hence
more attractive, appears to be the option according to which no major
difference in metal enrichment exists between
 field and clusters, and  the global
metallicity of the present day universe is nearly the same as
that of the only place where we can thoroughly measure it, i.e. in galaxy
clusters where it is $\sim 1/3$ times solar.
If so, there should be a comparable share of metals in the IGM,
as there is in the cluster ICM, i.e., just like most of the baryons,
most of the metals should reside
in the IGM rather than within field galaxies [30].
\vfil\newpage

\section{The Age of Spheroids and the Main Epoch of Metal Production}

The cluster abundances as illustrated in the previous section don't say
much about the cosmic epoch when the bulk of the cluster metals were
produced and dispersed through the ICM. The only constraint comes from the
iron abundance in moderate redshift clusters ($z\simeq 0.5$) being
virtually identical to that of local clusters (see Fig. 1),
hence the bulk of iron had to be manufacted at $z\ifmmode{\mathrel{\mathpalette\@versim>} 0.5$.
A much more stringent constraint comes from current age estimates of
the dominant stellar populations in cluster ellipticals, the likely
producers of the bulk of the metals, and from other
fossil evidences which are schematically reported in this section.

A first tight constraint on the formation epoch of stars in cluster
ellipticals came from Bower, Lucey, \& Ellis [6] who noted that the very tight
color-$\sigma$ relation followed by galaxies in the Virgo and Coma
clusters demostrates that at least
{\it cluster} ellipticals are made of very old stars, with the bulk
of them having formed at $z\ifmmode{\mathrel{\mathpalette\@versim>} 2$.
This result had the merit to cut short inconclusive discussions on the
age of ellipticals based on matching synthetic spectra to individual
galaxies, and showed instead that the homogeneity of elliptical
populations sets tight, almost model independent age constraints.
Along these lines, evidence supporting an early formation of ellipticals
 has greatly expanded over the last few years. This
came from the tightness of the fundamental plane
relation for ellipticals in local clusters [33],
from the tightness of the color-magnitude relation for ellipticals in
clusters up to $z\sim 1$ (e.g., [1], [37]), and from the
modest shift with increasing redshift in the zero-point of the fundamental
plane, Mg$_2-\sigma$, and color-magnitude relations of cluster
ellipticals (e.g., [4], [10], [12], [18], [26], [37], [42]).
All these studies
agree in concluding that most stars in cluster ellipticals formed at $z\ifmmode{\mathrel{\mathpalette\@versim>}
3$, though the precise value depends on the adopted cosmological
parameters (e.g. in a $\Lambda$-dominated universe this constraint
would become $z\ifmmode{\mathrel{\mathpalette\@versim>} 2$).

\begin{figure}
\vskip-1truecm
\centerline{\psfig{file=renzini-fig3.eps,width=15cm,angle=0,height=20cm}}

\caption{The Mg$_2-\sigma$ relation for a sample of early-type
galaxies (upper panel), as well as for the field, group and cluster
subsamples (lower panels). The corresponding number of objects,
the slope, and the zero-point (z.p.) are shown in the upper
left corner of each panel. he least squares fits to the Mg$_2-\sigma$ relation
are also shown. Note that for the three subsamples the slope as
derived for the total sample was retained, and only the zero-point was
determined [5].}
\end{figure}

This cogent  result has been established for cluster ellipticals, but
additional fossil evidence argues for its validity for field
ellipticals as well as for most bulges of spirals, i.e. in general
for the vast majority of galactic {\it spheroids}.  Indeed, as shown
in Fig. 3, field early-type galaxies (ellipticals and S0's) are found
to follow virtually the same Mg$_2-\sigma$ relation of their cluster
counterparts, with the age difference being less than $\sim 1$ Gyr
[5], while most bulges follow the same
Mg$_2-\sigma$ and color-magnitude relations of ellipticals [17].
Even more directly, the bulge of our own
Galaxy is found to be dominated by stars which age is
indistinguishable from that of the Galactic halo, or $\ifmmode{\mathrel{\mathpalette\@versim>} 12$ Gyr
[25], and similar conclusions have been reached for
the bulge of M31 [31], [16].
 It is also worth noting that the close
resemblance of cluster and field early-type galaxies illustrated in
Fig. 3 indicates that metal enrichment is not detectably
different in cluster and in field galaxies.

 With spheroids containing at least 30\% of all
stars in the local universe [27], [35], or even more [13], one can
conclude that at least
 30\% of all stars and metals have formed at $z\ifmmode{\mathrel{\mathpalette\@versim>} 3$ (Renzini [32];
see also Dressler \& Gunn [11]).
This is several times more than suggested by a conservative
 interpretation of the early attempt at tracing the cosmic history of
 star formation, either empirically [21] or from
theoretical simulations, e.g., [3].
Yet, it is is more in line with
recent direct estimates from the spectroscopy of Lyman-break galaxies
[38], where the cosmic SFR runs flat for $z\ifmmode{\mathrel{\mathpalette\@versim>} 1$,
as in one of the options offered by the models of Madau et al. [22].

At first sight these inferences from the local, fossil evidence appear
to conflict with some model predictions, and with some of the current
 interpretations of the direct observations at $z\ifmmode{\mathrel{\mathpalette\@versim>}
3$.  For example, the standard CDM model of the Durham group predicts
 that only $\ifmmode{\mathrel{\mathpalette\@versim<} 5\%$ of stars have
formed by $z=3$ [3], [7]. The supernova-feedback free parameter of this
specific CDM model was tuned to reproduce a $z=0$ galaxy luminosity
function that runs flat at the faint end. Tuning instead to the Zucca
et al. [45] luminosity function, which is steep at the faint end,
allows to make much more stars at early times (Frenk,
this conference), in somewhat better agreement
with the constraint sets by the old age of galactic spheroids.

With $\sim 30\%$ of all stars having formed at $z\ifmmode{\mathrel{\mathpalette\@versim>} 3$, and the
metallicity of the $z=0$ universe being $\sim 1/3$ solar, it is
straightforward
to conclude that the global metallicity of the $z=3$ universe had to
be $\sim 1/3\times 1/3\sim 1/10$ solar,
or more [32]. Damped Ly$_\alpha$
systems (DLA) may offer an opportunity to check this prediction, though they
may provide a vision of the early
universe that is biased in favor of cold, metal-poor gas that has been
only marginally affected by star formation and metal pollution.
Metal rich objects such as giant starbursts that would be dust obscured,
the metal rich
passively evolving spheroids, and the hot ICM/IGM obviously do not
enlist among DLAs.
Yet, in spite of these limitations the average metallicity of the DLAs
at $z=3$ appears to be
$\sim 1/20$ solar (Pettini et al. [28], see their Fig. 4), just a factor of
2  below the expected value from the {\it fossil evidence}.
This is still much higher than the extreme lower limit
$Z\simeq 10^{-3}Z_\odot$ at $z=3$ as
inferred from Ly$_\alpha$ forest observations [36]. Ly$_\alpha$ forest material
is believed to contain a major fraction of cosmic baryons at high $z$, hence
(perhaps) of metals. There is therefore a potential conflict with
the estimated global metallicity at $z\simeq 3$, and the notion of
Ly$_\alpha$ forest metallicity being representative of the the
universe metallicity at this redshift.
Scaling
down from the cluster yield, such low metallicity was achieved when
only $\sim 0.3\%$ of stars had formed, which may be largely
insufficient to
ionize the universe and keep it ionized up to this redshift [20].
This suggests that Ly$_\alpha$ forest may not trace the
mass-averaged metallicity
of high redshift universe, and that the universe was very
inhomogeneous at that epoch. The bulk of metals would be partly locked
into stars in the young spheroidals, partly would reside
in a yet undetected hot
IGM, a phase hotter than the Ly$_\alpha$ forest phase.
\par\bigskip
 \noindent
I would like to thank Len Cowie and Sandra Savaglio for useful
 discussions on the chemestry of DLAs and the Ly$_\alpha$ forest.

\begin{bloisbib}
\bibitem{} Aragon-Salamanca, A., Ellis, R.S., Couch, W.J. \&
     Carter, D. 1993, MNRAS, 262, 764
\bibitem{} Arimoto, N., Matsushita, K., Ishimaru, Y., Ohashi, T., \&
    Renzini, 1997, ApJ, 477, 128
\bibitem{} Baugh, C.M., Cole, S.,  Frenk, C.S., \& Lacey, C.G. 1997,
    Astro-ph 9703111
\bibitem{} Bender, R., Saglia, R.P., Ziegler, B., Belloni, P., Greggio, L.,
     Hopp, U. \& Bruzual, G.A. 1997, ApJ, 493, 529
\bibitem{} Bernardi, M., Renzini, A., da Costa, L.N., Wegner, G.,
           Alonso, M.V., Pellegreini, P.S., Rit\'e, C., \& Willmer,
           C.N.A. 1998, ApJL, in press (astro-ph/9810066)
\bibitem{} Bower, R.G., Lucey, J.R. \& Ellis, R.S. 1992, MNRAS, 254, 613
\bibitem{} Cole, S., Aragon-Salamanca, A., Frenk, C.S., Navarro, J.S., \&
    Zepf, S.E. 1994, MNRAS, 271, 781
\bibitem{} Davies, R.L., Sadler, E.M., \& Peletier, R.F. 1993, MNRAS,
    262, 650
\bibitem{} Davis, D.S., Mulchaey, J.S., \& Mushotzky, R.F. 1998,
           astro-ph/9808085
\bibitem{} Dickinson, M. 1995, in Fresh Views of Elliptical Galaxies,
         ed. A. Buzzoni, A. Renzini, \& A. Serrano, ASP Conf. Ser. 86,
         283
\bibitem{} Dressler, A., \& Gunn, J. E. 1990, ASP. Conf. Ser. 10, 204
\bibitem{} Ellis, R.S., Smail, I., Dressler, A., Couch, W.J., Oemler,
     A. Jr., Butcher, H., \& Sharples, R.M. 1997, ApJ, 483, 582
\bibitem{} Fukujita, M., Hogan, C.J., \& Peebles, P.J.E. 1998,
          astro-ph/9712020v2
\bibitem{} Greggio, L. 1997, MNRAS, 285, 151
\bibitem{} Ishimaru, Y., \& Arimoto, N. 1997, PASJ, 49, 1
\bibitem{} Jablonka, P., Bridges, T.J., Sarajedini, A., Meylan, G.,
         Maeder, A., \& Meynet, G. 1998
\bibitem{} Jablonka, P., Martin, P., \& Arimoto, N. 1996, AJ, 112, 1415
\bibitem{} Kodama, T., Arimoto, N., Barger, A.J., \&
      Aragon-Salamanca, A. 1998, astro-ph/9802245
\bibitem{} Kauffmann, G. 1996, MNRAS, 281, 487
\bibitem{} Madau, P. 1998, astro-ph/9807200
\bibitem{} Madau, P., Ferguson, H.C., Dickinson, M.E., Giavalisco, M.,
     Steidel, C.C., \& Fruchter, A. 1996, MNRAS, 283, 1388
\bibitem{} Madau, P., Pozzetti, L., \& Dickinson, M. 1997, Astro-ph 9708220
\bibitem{} Mushotzky, R.F. 1994, in Clusters of Galaxies, ed. F. Durret,
     A. Mazure, \& J. Tran Thanh Van (Gyf-sur-Yvette: Editions
     Fronti\`eres), p. 167
\bibitem{} Mushotzky, R.F., et al. 1996, ApJ, 466, 686
\bibitem{} Ortolani, S., Renzini, A., Gilmozzi, R., Marconi, G., Barbuy, B.,
     Bica, E., \& Rich, R.M., 1995, Nat, 377, 701
\bibitem{} Pahre, M.A., Djorgovski, S.G., \& de Carvalho, R.R. 1997, in
     Galaxy Scaling Relations: Origins, Evolution and Applications,
     ed. L. da Costa \& A. Renzini (Berlin: Springer), p. 197
\bibitem{} Persic, M., \& Salucci, P. 1992, MNRAS, 258, 14p
\bibitem{} Pettini, M., Smith, L.J., King, D.L., \& Hunstead,
     R.W. 1997, ApJ, 486, 665
\bibitem{} Renzini, A. 1994, Galaxy Formation, ed. J. Silk \&
     N. Vittorio (Amsterdam: North Holland), p. 303
\bibitem{} Renzini, A. 1997, ApJ, 488, 35
\bibitem{} Renzini, A. 1998a, AJ, 115, 2459
\bibitem{} Renzini, A. 1998b, astro-ph/9801209, in The Young Universe,
     ed. S. D'Odorico, A. Fontana, E. Giallongo, ASP Conf. Ser. 146, 298
\bibitem{} Renzini, A., \& Ciotti, L. 1993, ApJ, 416, L49
\bibitem{} Renzini, A., Ciotti, L., D'Ercole, A., \& Pellegrini,
     S. 1993, ApJ, 419, 52
\bibitem{} Schechter, P.L., \& Dressler, A. 1987, AJ, 94, 563
\bibitem{} Songaila, A. 1997, ApJ, 490, L1
\bibitem{} Stanford, S.A., Eisenhardt, P.R., \& Dickinson, M. 1998,
           ApJ, 492, 461
\bibitem{} Steidel, C.C., Adelberger, K.L., Giavalisco, M., Dickinson,
           M., \& Pettini, M. 1998, preprint
\bibitem{}  Thielemann, F.-K., Nomoto, K., \& Hashimoto, M. 1996, ApJ,
     460, 408
\bibitem{}  Thomas, D., Greggio, L., \& Bender, R. 1997, Astro-ph 9710004
\bibitem{}  Tinsley, B.M. 1980, Fund. Cosmic Phys. 5, 287
\bibitem{} van Dokkum, P. G., Franx, M., Kelson, D. D.,
     \& Illingworth, G. D. 1998, ApJ, 504, L17
\bibitem{} White, S.D.M., Navarro, J.F., Evrard, A.E., \& Frenk,
     C.S. 1993,   Nat, 366, 429
\bibitem{} Woosley, S.E., \& Weaver, T.A. 1995, ApJS, 101, 181
\bibitem{} Zucca, E., et al. 1997, A\&A, 326, 477
\end{bloisbib}

\end{document}

\bibitem{GG} Genius G., 1995, {\it preprint}
\bibitem{jk} Kirk J., 1995, \aa {100} {44}
\bibitem{wd} Mouse M., Duck D., 1935, \apj {1} {2}
\bibitem{conf} Smith G., 1995, in {\it Proceedings of the nth Conference}
     p. 129, eds Durand et al., World Editions
\bibitem{} Lennon, J., McCartney, P., Harrison, G. \& Starr, R. 1962, \mnras
{100} {1}
\end{bloisbib}
\vfill
\end{document}

Now, what is the main epoch of metal production in the
universe? or, equivalently, when the cosmic SFR has reached its peak
value?  The most straightforward approach to answer these questions is
certainly offered by the direct determination of the global SFR as a
function of redshift (the
Madau diagram). However, as well known, direct determinations of the global
SFR at high redshift encounters two major difficulties: 1) dust
obscuration, which on average may be higher than in low redshift
galaxies since star formation may predominantly take  place in major dust
enshrouded, {\it bulge-forming} starbursts, rather than more
quiescently as later in disks and
irregulars; and 2) alternatively, a sizable fraction of the global SFR
may take place in small entities below detection threshold, as
predicted by current hierarchical models (e.g. Kauffmann 1996). Therefore, the
complementary approach based on the fossil record at low redshift
will provide a vision of the early universe from a different point of
view, hence subject to different biases and uncertainties. Ultimately, the two
approaches will have to lead to converging results.

As far as direct observations are concerned, the fossil evidence
supports the interpretation as {\it lower limits} of the SFR estimates
at $z\ifmmode{\mathrel{\mathpalette\@versim>} 3$ (e.g. Madau et al. 1996).  Of the two models presented
by Madau et al. (1997) the one in which the SFR slowly increases all
the way to $z=5$ is favored. It remains to be ascertained
whether the undetected star formation activity is obscured by dust, or
dispersed in small fragments below the  threshold of the Hubble
Deep Field.


Extended ionized-gas structures in Seyfert 2 galaxy Mrk 78

Abstract: Search for and study of extended emission-line regions (EELRs) related to AGN in early-type galaxies is interesting to probe the history of nuclear ionization activity and also to understand the process of external gas accretion. In this work, we present observations of the EELR in Mrk 78 obtained at the 6-m Russian telescope using the long-slit and 3D spectroscopy methods. We show that ionized-gas clouds at the 12-16 kpc projected distances from the nucleus are ionized by the AGN radiation. Also we have checked if the galaxy appearing in the optical images in the immediate neighbourhood of Mrk 78 near the external clouds is a dwarf companion or a part of a tidal structure. However, the spectrum of this galaxy, SDSS J074240.37+651021.4, obtained at the 6-m telescope corresponds to the distant background galaxy.

\section{Introduction}
\label{intr}

According to the unified model of active galactic nuclei (AGN), the radiation of the central accretion engine comes out in  broad axisymmetrical cones collimated by the circumnuclear obscuring torus. The impact of AGN ionizing radiation together with radio jet  mechanical energy  on the surrounding interstellar medium appears in the shape of ionizing cones from hundreds pc to a few kpc in typical size.  In some fraction of Seyfert galaxies,  extended emission-line regions (EELRs) are also detected at radial distances of tens of kpc. EELRs can be used to
measure the AGN energetic output at both time and spatial scales and also to study the gas distribution in the galactic environment. Indeed, a significant faction of EELR traces the off-plane gas related with tidal debris \citep{Keel2012} or even the gas in a companion galaxy \citep{Keel2019}. One of the  intriguing cases is the well-known AGN  Mrk~6: deep imaging and spectroscopic observations revealed an extended system of gaseous filaments up to 40 kpc from the nucleus accreted by the Seyfert galaxy \citep{Smirnova2018}. Therefore, the detailed study of EELRs is topical even around well-studied nearby active galaxies.

The Seyfert 2 galaxy Mrk~78 is a good illustration of the picture described above. The spectrophotometric properties of the inner ($r<4-5''$, or 3--4\,kpc) elongated emission-line structure studied at the HST \citep{2004AJ....127..606W,2006MNRAS.372..961K} can be explained in terms of jet-cloud interaction together with the mass outflow and a cone of UV-radiation emerging from the Sy nucleus \citep{Fischer2011}, see the corresponding [OIII] HST image in Fig.~\ref{MPFS}. Ground-based  narrow-band images obtained by \citet{DeRobertis1987}  and \citet{Pedlar1989} also reveal a weak extended ionized-gas emission outside the central region with an  extensive asymmetric structure at $r>15''$ in the south-west. The presence of the gas emitting  in  [OIII] in the western direction up to 12-14\,$''$ was  confirmed in spectroscopic observations by \citet{Unger1987}. Also,   \citet{AfanasievSilchenko1991} detected the [OIII] line emission even at $r=18-20''$. \citet{Pedlar1989} have shown that the weak emitting narrow-line region extending out to 10 kpc to the west is consistent with AGN ionizing cones, however, the actual ionization  state of the most distant part of the EELR is still unknown. Mrk~78 was included in the sample of the confirmed EELR selected from the Galaxy Zoo survey, but  the follow-up spectra were obtained along  $PA=90\,\degr$\citep{Keel2012}, therefore, the properties of the southwestern external emission knots were not measured.

In this paper, we analysed the  large-scale ionized gas distribution and spectral properties of the most distant parts of the EELR in Mrk~78 derived from new long-slit and 3D spectroscopic observations at the  6-m telescope of Special Astrophysical Observatory of Russian Academy of Sciences (SAO RAS). The gas and stellar kinematics will be considered in the forthcoming paper. In the present paper, we accepted the distance to the galaxy as $165$\,Mpc with a corresponding scale of 0.80 \,kpc/$''$ (according to  NED\footnote{http://ned.ipac.caltech.edu/}).

\begin{figure}[h]
\centerline{
\includegraphics[width=0.55\textwidth]{HST.eps}
\includegraphics[width=0.35\textwidth]{MPFS.eps}
}
\caption{Left: the HST/FOC image with the F520M filter, shows the  [OIII]  emission \citep[see][]{Fischer2011}; right: the MPFS image in the [OIII]${\lambda}5007$ line with the superimposed continuum isophotes.}
\label{MPFS}
\end{figure}

\section{Observations and data reduction}
\label{obs}

Observations were carried out at the prime focus of the 6-m SAO RAS telescope with the multi-mode focal reducers: SCORPIO \citep{Afanasiev2005} and SCORPIO-2 \citep{Afanasiev2011}. We used the 3D spectroscopy with the Fabry-Perot interferometer (FPI), long-slit spectroscopy, and direct-image modes. Also, we used the data obtained with the integral-field Multi-Pupil Field Spectrograph \citep{2001sdcm.conf..103A}. The log of observations and parameters of the used instruments are listed in the Table\,\ref{t1}.

SCORPIO and SCORPIO-2 both have the $6\farcm1$ fields of view with the 0.36-arcsec pixel size. The slit has the $6\farcm1$ length and the $1''$ width. The scanning FPI mapped the spectral region around the redshifted [OIII]${\lambda}4959$ emission line. We made the binning of these data to reach a higher signal-to-noise ratio  in the resulting data cube with 0.7\,$''$/px. The MPFS data cube had the 16x16 arcsec$^2$ field-of-view centered at the nucleus with the $1''/$spaxel spatial scale and covered a wide spectral range.

Calibration of the SCORPIO-2 and MPFS spectra into the absolute energy flux units was carried out using spectrophotometric standard stars observed in the same night as Mrk~78. The data reduction was performed in a  standard way (see \citet{Smirnova2018} for references on the software and algorithms).

Figure~\ref{fields} shows the maps derived from the Voight fitting of the [OIII] spectra in the FPI data cube: the monochromatic image, the velocity field, and the velocity dispersion map free from instrumental effects. The emission-line fluxes in the MPFS and long-slit data were derived from the single-component Gaussian fitting. The map of the [OIII] emission in the circumnuclear region according to the MPFS data is shown in Fig.~\ref{MPFS}. For the internal region, where well-known multi-component  lines were observed \citep[see][and references therein]{AfanasievSilchenko1991}, only a brighter peak was fitted in the MPFS and FPI data.

\begin{table}[t]
\small
\begin{center}
\caption{Log of observations at the SAO RAS telescope}
\label{t1}
\begin{tabular}{llllll}
\hline\hline
Date  & Instrument/mode & Exp. time, s & Sp. range & FWHM & Seeing\\
\hline
14/15 Mar 2007 & SCORPIO/FPI  & 32$\times$250 & [OIII]${\lambda}4959$ & 1.2\,\AA& 1\,\farcs2 \\
22/23 Oct 2008 & SCORPIO/DI  & 5$\times$120 & R$_C$-filter & & 1\,\farcs2 \\
03/04 Feb 2008 & MPFS  & $6\times1200$ & 3650--5900\,\AA & 5\AA& 1\,\farcs3 \\
03/04 Feb 2008 & MPFS  & $6\times1200$ & 4300--7380\,\AA & 5\AA& 1\,\farcs3 \\
12/13 Feb 2019 & SCORPIO-2/LS & 3x900 & 3650--8530\,\AA  & 7\AA & 2\,\farcs7  \\
\hline\hline
\end{tabular}
\end{center}
\end{table}

\section{External gas clouds}
\label{prop}

The ionization state of the EELR inside  the galaxy disc ($r<10-15''$) was analysed in the papers cited above in Sec.~\ref{intr}. However, the FPI map also revealed  separate ionized-gas clouds (labeled as C1 and C2 in Fig.~\ref{fields}) at a distance of ${\sim}$15-20${\farcs}$ (${\sim}$12-16 kpc) southwest of the galaxy. The line-of-sight  velocity field shows that the cloud velocities lie in the velocities range of the disc gas (10900--11500~ km\,s$^{-1}$). Whereas the velocities of C1 and C2 by 50--200~ km\,s$^{-1}$ exceed the nearest S-W side of the galactic disc. This indicates that both clouds are moving in the Mrk~78 gravitational potential, but rotate out of the galaxy plane on retrograde orbits possibly.

To determine the C1 and C2 ionization state, the SCORPIO-2 long-slit spectra were obtained along PA=127$\degr$. The slit position has been selected according to the centers of both clouds and a neighbouring galaxy with an unknown spectroscopic redshift (see Sec.~\ref{back}).

\begin{figure}[h]
\centerline{
\includegraphics[height=0.43\textwidth]{R_slit_ccg.eps}
\includegraphics[height=0.43\textwidth]{mono_color.eps}
}
\centerline{\includegraphics[height=0.53\textwidth]{vr_color.eps}
\includegraphics[height=0.53\textwidth]{disp_color.eps}
}
\caption{Mrk~78 observations. Top left: the SCORPIO $R$-band image with the~long-slit position, a candidate  companion galaxy marked as `G'; top right: the map of the $[OIII]{\lambda}4959$ emission derived from the FPI data; bottom: the line-of-sight velocity field  in the [OIII] emission line (left) and the  velocity dispersion map with the marked ionization cone position (right). The colour bars are in the km\,s$^{-1}$, the external gaseous clouds are marked as `C1' and `C2'.}
\label{fields}
\end{figure}

\begin{figure}[h]
\centerline{
\includegraphics[width=\textwidth]{ident_lines.eps}
}
\centerline{
\includegraphics[width=\textwidth]{integ_spec.eps}
}
\caption{Top: identification of the emission lines in the ionized-gas clouds (C1 and C2) and the nearby galaxy (G) long-slit spectra; bottom: the integrated spectrum of both clouds.}
\label{spec}
\end{figure}

To determine the ionization mechanism of the gas in C1 and C2, the emission-line ratio diagnostic diagrams \citep[BPT-diagrams after][]{1981PASP...93....5B} were used with the branches from \citet{2006MNRAS.372..961K} to separate different mechanisms of  gas excitation: the HII-regions and AGN (Seyfert or LINER). These diagrams involve   the following  line flux  ratios independent from the interstellar extinction: [N II]${\lambda}6583/H{\alpha}$,  [OIII]${\lambda}5007/H{\beta}$, and $[S II]{\lambda}6717+6731/H{\alpha}$.

The integrated spectrum of the C1+C2 clouds is shown in Fig.~\ref{spec}, the BPT diagrams are presented in Fig. \ref{BPT}. We designated the EELR ionization mechanism according to the long-slit spectroscopy. To compare their ionization state to that of the circumnuclear region, the MPFS intergral-field spectroscopic data were used. It is clearly seen that the C1 and C2 ionization state is the same as that in the internal part of the EELR: the  Seyfert-type AGN radiation predominates. This result is in the agreement with the diagrams based on the  long-slit spectroscopy along the major galactic axis (PA=90$\degr$) presented in \cite{Keel2012}.

\begin{figure}[h]
\centerline{
\includegraphics[width=0.471\textwidth]{SII_color.eps}
\includegraphics[width=0.53\textwidth]{NII_color.eps}
}
\caption{BPT diagrams. The separate lines are taken from \citet{2006MNRAS.372..961K}. The circles denote the MPFS data for the circumnuclear region, the colour bar shows the projected distance of pixels on the MPFS maps from the nucleus. The diamond shows the line ratio of C1+C2 derived from the long-slit spectrum.}
\label{BPT}
\end{figure}

\cite{Fischer2011} estimated the parameters of the  ionized cones according to the HST  data for the inner ($r<4''$) brightest part of the emission-line structure. We overlaid the western side of this bicone (having an opening angle of ${\sim}$75${\degr}$) on the FPI dispersion map (Fig.~ \ref{fields}): C1 and C2 both fall into the projected cone borders. Therefore, the ionization state of the distant gas in the EELR is consistent with the model based on the analysis of the circumnuclear gas kinematics and ionization state. The ionized-gas  velocity dispersion in C1 and C2  (${40<\sigma<100}$~km\,s$^{-1}$, see Fig.~ \ref{fields})  is similar to the gas in the outer part of the galaxy disc and significantly smaller than that observed in the inner region, where jet-cloud interaction takes place. It means that the external part of the EELR is dynamically cold, i.e., the AGN highlight off-plane gas clouds are related  with the external accretion or tidal structures, but are not the result of the AGN outflow. The similar picture is observed in  most EELRs in Seyfert galaxies \citep[see] [and references therein]{Keel2015}.

\section{Possible dwarf companion}
\label{back}

Assumptions that Mrk~78 is a post-merging system were considered in a number of papers since the eighties of the last century. Morphological peculiarities on different spatial scales were considered as the galaxy interaction imprints: the nuclear dust lane \citep{2004AJ....127..606W},  the disturbed external isophotes in the S-W part of the galactic disc \citep{DeRobertis1987, Pedlar1989,  Keel2015}, the asymmetric distribution of the external [OIII] emission \citep{Unger1987,AfanasievSilchenko1991}.

\citet{DeRobertis1987} suggested that the disturbed external isophotes are caused by passing through a companion galaxy and the secondary emission of this region due to the companion galaxy nucleus, but the forthcoming HST data disproved this idea \citep{Fischer2011}.

 The SDSS DR15 and SCORPIO $R$-band  images reveal the galaxy SDSS J074240.37+651021.4 with an unknown redshift  in the nearest neighbourhood of the C1 and C2 clouds at a distance of ${\sim}$18$''$ from the Mrk~78 nucleus (marked as `G' in  Fig.~\ref{fields}).  We have checked the possibility that this galaxy being a dwarf companion  of  Mrk~78   or remnant of a tidal structure. The obtained long-slit spectra (see the slit position in Fig.~\ref{fields}, PA=127$\degr$) show the continuum of  `G' with weak emission lines ($\lambda$=5230\AA~   and $\lambda$=7016\AA) marked with arrows in Fig.~\ref{spec}. We interpreted them as $[OII]{\lambda}3727$ and $[OIII]{\lambda}5007$ at the redshift $z=0.382$. In this case, the  $H{\beta}$ should be blended with  [S II]${\lambda}6731$ in the C1 and C2 clouds. This spectroscopic redshift is comparable to the photometric SDSS DR15 \citep{SDSSDR15} redshift: $z=0.367\pm 0.130$. Therefore, it is a distant background galaxy.

\section{Summary}
\label{sum}

The 6-m telescope observations   reveal a pair of the ionized-gas  clouds   southwest from the stellar disc of  Mrk~78 at the 12--16\,kpc projected distance from the AGN located beyond the galactic disc plane. The BPT diagrams  have demonstrated  that the AGN radiation makes the main contribution to the ionization state of these clouds. Moreover, the location of both clouds together with other external [OIII]-emitting structures in the west side of the Mrk~78 disc is in a good agreement  with the  ionization cone borders proposed by \cite{Fischer2011} according only to the circumnuclear region structure and kinematics.

The AGN outflow as a source of the described ionized-gas clouds   is the contrary to the dynamically cold gas in this part of the EELR. The most possible scenario is an external accretion, a minor merging, or a tidal interaction with the companion. The galaxy  SDSS~J074240.37+651021.4 at a projected distance of ${\sim}$\,18\,$''$ from the Mrk~78 nucleus in the southwest direction was   considered  a candidate. However, the  long-slit spectroscopy data gave a redshift of 0.382;  it is a distant background galaxy. Thereby, the origin of the gas in the outer  part of the Mrk~78 EELRs is still unclear.

\acknowledgements
This study was supported by the Russian Science Foundation, project no. 17-12-01335 `Ionized gas in galaxy discs and beyond the optical radius' and based on the observations conducted with the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences carried out with the financial support of the Ministry of Science and Higher Education of the Russian Federation.  Authors are  grateful
to Victor Afanasiev for his interest to this work.
This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.


The solar photospheric abundance of europium. Results from CO5BOLD 3-D hydrodynamical model atmospheres

Abstract: Context. Europium is an almost pure r-process element, which may be useful as a reference in nucleocosmochronology. Aims. To determine the photospheric solar abundance using CO5BOLD 3-D hydrodynamical model atmospheres. Methods. Disc-centre and integrated-flux observed solar spectra are used. The europium abundance is derived from the equivalent width measurements. As a reference 1D model atmospheres have been used, in addition. Results. The europium photospheric solar abundance is 0.52 +- 0.02 in agreement with previous determinations. We also determine the photospheric isotopic fraction of Eu(151) to be 49 % +- 2.3 % from the intensity spectra and 50% +-2.3 from the flux spectra. This compares well to the the meteoritic isotopic fraction 47.8%. We explore the 3D corrections also for dwarfs and sub-giants in the temperature range ~5000 K to ~6500 K and solar and 1/10--solar metallicities and find them to be negligible for all the models investigated. Conclusions. Our photospheric Eu abundance is in good agreement with previous determinations based on 1D models. This is in line with our conclusion that 3D effects for this element are negligible in the case of the Sun.

\section{Introduction}
\label{intro}

Europium (Z=63) is formed through neutron captures
on {\sl seed}-nuclei; in the seminal paper of
\citet{burb} it was assigned both to the {\sl s}-process
(slow neutron capture) and to the
{\sl r}-process (rapid neutron capture).
According to the current understanding of
{\sl r}-process nucleosynthesis, europium
is  an
almost {\sl pure} {\sl r}-process element.
About 95-97\% of the Eu
in the solar system is
contributed by the {\sl r-}process \citep{arl,burris}.
See the introduction in the paper by \citet{MG2000}
on the significance of the Eu/Ba ratio for assessing
the relative contributions of the {\sl r-}process
and {\sl s}-process.

The trend of the [Eu/Fe] abundance ratio as a function of [Fe/H] in
the Milky Way stars (both halo and disk) seems to mimic the same
behaviour of the [$\alpha$/Fe] ratio, with an enhanced value in the
metal-poor stars and a decrease for [Fe/H]$>$-1 dex down to solar
values. Measurements of this abundance ratio in extra-galactic
stars highlight a decoupling between [Eu/Fe] and [$\alpha$/Fe] ratios,
with over-solar [Eu/Fe] at high metallicity and solar or sub-solar
[$\alpha$/Fe] values, as observed in Sagittarius \citep{boni00}, in
the Small Magellanic Cloud \citep{hill97} and Large Magellanic Cloud
\citep{hill95}.

The ratio of the two {\sl r-}process elements Eu/Th
\footnote{The isotope $^{232}$Th is radioactive with a half-life of 14.05 Gyr.}
is potentially an interesting chronometer, provided
the production ratio of the two nuclei can be reliably
predicted theoretically
\citep[see][for an extensive review of the {\sl r-}process]{cowan}.

In this paper we reconsider the Eu solar abundance in the light of the
recent progress of the 3-dimensional model atmosphere computations,
measuring the solar abundance both with 1D and 3D models
in order to assess the impact of this new generation of models on
the solar Eu abundance.

The Eu solar abundance reported in the compilation \citet{gs98}
is  A(Eu)=0.51$\pm$0.08
\footnote{We adopt the spectroscopic notation A(X)=log(N(X)/N(H))+12.}
and A(Eu)=0.55$\pm$0.02 dex, from the solar photosphere and
from meteorites, respectively. The analysis
of \citet{lawler} yields A(Eu)=0.52$\pm$0.01 dex, adopting new atomic data for the
Eu lines and taking into account the hyperfine structure of the lines employed.

\section{Solar europium chemical analysis}

\subsection{Theoretical tools}

In this paper we derived the solar Eu photospheric abundance by
using 3D model atmospheres computed with the \ensuremath{\mathrm{CO}^5\mathrm{BOLD}}\ code
(COnservative COde for the COmputation of COmpressible COnvection in a
BOx of L Dimensions with L=2,3) \citep{freytag02,wed}
and we compared it to the 1D model results.
\ensuremath{\mathrm{CO}^5\mathrm{BOLD}}\ solves the coupled non-linear equations of compressible
hydrodynamics including non-local frequency-dependent radiation transport
for a small volume located at the stellar surface
(for technical issues, see the on-line manual available on
http://www.astro.uu.se/$\sim$bf/co5bold/index.html).
The atmospheric flow field is sampled in
equal temporal intervals each of which we call
a ``snapshot''. In total,
25 snapshots were selected from a \ensuremath{\mathrm{CO}^5\mathrm{BOLD}}\ simulation to represent the
solar photosphere.\\
As reference  we also adopted several 1D
solar models:

\begin{itemize}
\item  The 1D solar model computed by
F. Castelli\footnote{http://www.user.oats.inaf.it/castelli/sun/ap00t5777g44377k1asp.dat}
with the ATLAS9 code
adopting the solar abundances of \citet{asplund},
\ensuremath{T_\mathrm{eff}} =5777 K, $\log{g}$=4.4377
and micro-turbulence velocity of 1 km/s.
\item  The semi-empirical solar atmosphere model derived by \citet{hm}.
\item  A solar 1D atmospheric model
obtained by the temporal and horizontal average of the 3D structure over surfaces
of equal Rosseland optical depth. The comparison between the 3D \ensuremath{\mathrm{CO}^5\mathrm{BOLD}}\
model and this kind of {\sl average} 3D model allows to estimate the
influence of the fluctuations around the mean stratification
on the line formation process.
\item  A solar 1D atmospheric model computed with a Lagrangian hydrodynamical
code LHD (see \citealt{caffau07}), including the same opacities and equation-of-state
adopted by the \ensuremath{\mathrm{CO}^5\mathrm{BOLD}}\ 3D code. LHD code treats the convection
with the standard mixing length theory (MLT) in the
formulation given by \citet{mih}.
The use of this kind of models
allows to compare directly these 1D models with the 3D \ensuremath{\mathrm{CO}^5\mathrm{BOLD}}\ models,
erasing the systematics due to different physical assumptions.
\end{itemize}

The spectral synthesis from ATLAS and HM models
are performed by using the SYNTHE code \citep{K93,kurucz} in
its Linux version \citep{sbordone04,sbordone05}.
For \ensuremath{\mathrm{CO}^5\mathrm{BOLD}}\ and LHD models the
Linfor3D\footnote{http://www.aip.de/~mst/Linfor3D/linfor\_3D\_manual.pdf}
code is used.

\subsection{Observational material}
The present study is based on two sets of high-resolution, high signal-to-noise ratio
spectra of solar flux and disc-centre intensity:

\begin{itemize}
\item {\sl Solar Flux --} We adopted the solar flux spectra
of \citet{neckel84}
and of Kurucz 2005\footnote{See http://kurucz.harvard.edu/sun.html}.
\item {\sl Solar Intensity --} As centre disc solar intensity spectra,
we used the intensity of \citet{neckel84} and that of
\citet{deb}\footnote{http://bass2000.obspm.fr/solar\_spect.php}.
\end{itemize}

We selected 5 \ion{Eu}{ii} optical spectral lines from the list of \citet{lawler}.
It is worth to note that not all these solar spectra are useful to measure
the selected features.
In the following, we describe briefly the single \ion{Eu}{ii} features used in this
analysis and the corresponding adopted solar spectrum:

\begin{itemize}
\item 412.972 nm -- Strong \ion{Eu}{ii} line, with a weak blending on
the blue side; this feature exhibits the same shape in all the adopted solar spectra,
without telluric contamination.
\item 604.951 nm -- Clean feature without particular difficulties,
it is not blended or contaminated by nearby lines to cause problems in the
placement of the continuum.
\item 664.519 nm -- This is one of the strongest optical transitions and
is commonly used to infer
the \ion{Eu}{ii} abundance. In all the four solar observed spectra we consider
this line results not blended with telluric features.
This feature exhibits in the solar spectrum a
blending in its red wing, due to the presence of the weak features of
Cr I and Si I.
\item 707.709 nm -- This weak line shows an H$_{2}$O telluric line
contamination on its red side and only in the Delbouille spectrum this blending
seems to be less severe.
\item 721.756 nm -- This line is measurable only
in the Delbouille spectrum, instead the other adopted solar spectra
exhibit a strong blending of this feature with a H$_{2}$O telluric line,
completely absent in the Delbouille spectrum.
\end{itemize}

\subsection{Chemical analysis}

\subsubsection{Analysis}

The chemical analysis of the selected Eu II features has
been performed adopting the atomic parameters
for the \ion{Eu}{ii} lines by \citet{lawler} and summarised
in Table 1. The \ion{Eu}{ii} spectral lines display
significant hyperfine structures. We included in the
line list hyperfine structure and isotopic splitting,
adopting the meteoritic isotopic ratio
\footnote{This element has two isotopes,
with Z=151 and Z=153 with meteoritic
abundance of 47.8\% and 52.2\% respectively \citep{ag89}}
and the hyperfine constants A and B measured by \citet{lawler}.
The calculation of the hyperfine structure was done using the
code LINESTRUC, described by \citet{wal}.
All the hyperfine components
for each \ion{Eu}{ii} feature, computed without
the assumption of a specific isotopic ratio, are
available in the on-line version.
We did not take into account possible
NLTE effects which are different between intensity and flux spectra
and could explain the small positive difference between intensity
and flux abundances: \citet{mash} analysed the NLTE effects for the
resonance \ion{Eu}{ii} line at 421.9 nm in solar-like stars, finding
a NLTE correction of $\sim$0.04 dex.

The solar Eu abundance was derived from the curve of growth of each line
calculated with Linfor3D, adopting a meteoritic isotopic ratio.
The equivalent width (EW) of the \ion{Eu}{ii} lines was measured with a Gaussian fit by
using the IRAF task SPLOT, adopting the deblending option.
The 3D models include only the \ion{Eu}{ii} lines, without the contribution
of possible blending features. The choice to infer the abundance by using the
EW measurement comes from the inefficiency of the line profile fitting
with a 3D grid due to the lack of the weak blending components in the 3D
synthetic spectra.
This is due to the inability of the current version of Linfor3D\ to handle a large
number of lines.
In Table 2  we  provide our results for both 1D and
3D models, the 3D correction defined by \citet{caffau07} as
$A(X)_\mathrm{3D}$-$A(X)_\mathrm{\ensuremath{\mathrm{1D}_{\mathrm{LHD}}}}$, and the difference between
3D and \ensuremath{\left\langle\mathrm{3D}\right\rangle}\ models. We reported also the error ($\sigma_\mathrm{EW}$)
in the Eu abundance due to the uncertainty in the EW measurement (in order to
estimate this latter issue we performed EW measurements with different
continuum placements and deblending assumptions for each line),
typically of $\sim$0.02-0.03 dex (only the
\ion{Eu}{ii} line at 412.972 nm shows an error in the abundance
of $\sim$0.05 dex, probably due to the blending on the blue side).

To place solar 3D correction results in a wider context,
we computed 3D corrections of the 664.5 nm \ion{Eu}{ii} line
for flux spectra in F and G-type atmospheric stellar models.
We explored a parameter grid including \ensuremath{T_\mathrm{eff}}\ between 4980
and 6460 K, log~g=3.5,4.0,4.5 and [M/H]=0.0,--1.0.
The Eu abundance is scaled with respect to the metallicity
of the model, according to the solar ratio.
The reference solar Eu abundance is 0.52. The
results are listed in Table \ref{3dcor}. The majority of the
3D corrections (3D-\ensuremath{\mathrm{1D}_{\mathrm{LHD}}}) are negligible, and the largest is
just 0.011\,dex. The 3D correction related to the average
temperature profile (3D-\ensuremath{\left\langle\mathrm{3D}\right\rangle}) is in the range 0.01-0.02\,dex
for all models and it is larger
than the complete 3D correction.

As additional check to test the consistency of our results, we performed
a {\sl classical} 1D analysis on these 5 features.
This step is necessary to compare the results obtained by the LHD models
and the 1D models usually used in the chemical analysis.
To compute the abundance we used line profile fitting and employed
the line list from the Kurucz database, updated including the atomic parameters
for the \ion{Eu}{ii} lines.
This was done by using a code \citep{caffau} that
performs a linear interpolation in a synthetic spectra grid with
the Eu abundance as a free fitting parameter:
the final best-fit is obtained by
the numerical $\chi^2$ minimisation,  using MINUIT \citep{james}.
Even the line shift and the continuum placement can be a free
parameter to be adjust to optimise the fit.
Only for the two strong features (namely 412.9 and 664.5 nm) we adopted
a different version of this code, including as free fitting
parameters both the Eu abundance and
the fraction of the Eu isotope $^{151}$Eu with respect
to the total abundance, $\log(N(^{151}{\rm Eu})/N({\rm Eu}_\mathrm{tot}))$.

\section{Results and discussion}

The main results of this analysis are:
\begin{enumerate}

\item The 3D analysis, based on different high-resolution high
signal-to-noise solar spectra and by using the \ensuremath{\mathrm{CO}^5\mathrm{BOLD}}\ model, provided
a mean Eu photospheric abundance of A(Eu)=0.506 dex with a standard
deviation $\sigma$=0.008 for the flux spectra (by using the first three
spectral features) and A(Eu)=0.527 with $\sigma$=0.024 for the intensity
spectra at the disk-centre (by using all the five lines). As a final Eu
solar photospheric abundance we recommend A(Eu)=0.518 dex ($\sigma$=0.024).
This value comes from the average of all the measurements (both flux
and intensity, since they are very close), and
it is consistent with the previous 1D determinations.\\

\item The difference 3D-\ensuremath{\left\langle\mathrm{3D}\right\rangle}\  allows to estimate the
3D corrections due to the horizontal temperature
fluctuations (a component not taken into account in the classical
1D models). This correction is negligible for all  of the lines
considered, with an average value of
--0.009 dex ($\sigma$=0.016) and 0.010 dex ($\sigma$=0.003) for
flux and intensity respectively. This difference between the two
solar data-sets has been already observed in a previous 3D analysis
for sulphur \citep{caffau} and phosphorus \citep{caffau07} and
can be ascribed to the different atmospheric layers where intensity
and flux spectra originate (the centre disc intensity spectra
arise from deeper layers, where the temperature
fluctuations are more pronounced).\\

\item The difference 3D-\ensuremath{\mathrm{1D}_{\mathrm{LHD}}}\ allows to compare
3D and 1D models which employ the same physical assumptions,
like equation of state and opacities, and provides a
{\sl 3D correction}.
These
values are near to zero both   for flux and
intensity; with an average difference
of 0.004 dex ($\sigma$=0.013) and
0.021 dex ($\sigma$=0.004) respectively.
The 3D-\ensuremath{\mathrm{1D}_{\mathrm{LHD}}}\ corrections appear
to be systematically higher than  3D-\ensuremath{\left\langle\mathrm{3D}\right\rangle}\ corrections  with
a difference of about 0.010 dex.\\

\item Finally, as consistency check, we performed a classical
1D analysis by using the \citet{hm} and ATLAS models. We derived a
mean photospheric abundance for europium of A(Eu)=0.515 dex
($\sigma$=0.022) and 0.523 dex ($\sigma$=0.014) for disc-centre
intensity and flux, respectively, in good agreement with the
previous ones by \citet{ag89} of 0.51 dex and by \citet{lawler}
of 0.52 dex. Moreover, also the isotopic ratio
$(N(^{151}{\rm Eu})/N({\rm Eu}_\mathrm {tot}))$  computed from
the two strongest Eu lines results in good agreement with the
meteoritic ratio: 0.49 ($\sigma$=0.023)
and 0.50 ($\sigma$=0.023) for intensity and flux, respectively.
\end{enumerate}

There is no way to know {\em a priori} if the 3D effects
are important for any given line. A detailed calculation
has to be done in each case.
Solar abundances are widely used as a reference and their
implication goes beyond the pure chemical composition, but
touches field such as helioseismology and solar neutrino production.
The low solar abundances of \citet{asplund} have put some strain
on our understanding of both.
As suggested by \citet{bahcall05}, different measurements of
solar abundances, obtained using different observed
spectra and different solar models, allow a better
estimation of the systematic uncertainties.
In the case of Eu
we conclude that the 3D effects are
negligible in the Sun and  solar-like stars.
The scenario is very coherent, when we consider 3D
of the non-solar models. 3D corrections are negligible for both
solar and slightly metal-poor models.
Moreover, this is in line with the findings discussed
by \citet{steffen} that investigated granulation corrections in the
Sun for several elements. Despite Eu is not included in this
study, we can compare our results for Eu with their findings
for Sr. These two elements show very similar line formation properties
for spectral lines with similar excitation potential and oscillator
strength. Also for Sr, the corrections are negligible,
typically between --0.02 and +0.02.
Therefore also for Eu, like for S \citep{caffauS}
and P \citep{caffauP} we conclude that the use of
3D models does not imply a substantial downward
revision of the solar abundances with respect
to what was deduced from the use of 1D models.

\acknowledgements

We warmly thank the referee, A. J. Sauval, for his useful
suggestions.
The authors E.C., H.-G.L., P.B. acknowledge financial
   support from EU contract MEXT-CT-2004-014265 (CIFIST).
   We acknowledge use of the supercomputing centre CINECA,
   which has granted us time to compute part of the hydrodynamical
   models used in this investigation, through the INAF-CINECA
   agreement 2006,2007.