The total infrared luminosity may significantly overestimate the star formation rate of recently quenched galaxies

Abstract: The total infrared (IR) luminosity is very useful for estimating the star formation rate (SFR) of galaxies, but converting the IR luminosity into an SFR relies on assumptions that do not hold for all galaxies. We test the effectiveness of the IR luminosity as an SFR indicator by applying it to synthetic spectral energy distributions generated from three-dimensional hydrodynamical simulations of isolated disc galaxies and galaxy mergers. In general, the SFR inferred from the IR luminosity agrees well with the true instantaneous SFR of the simulated galaxies. However, for the major mergers in which a strong starburst is induced, the SFR inferred from the IR luminosity can overestimate the instantaneous SFR during the post-starburst phase by greater than two orders of magnitude. Even though the instantaneous SFR decreases rapidly after the starburst, the stars that were formed in the starburst remain dust-obscured and thus produce significant IR luminosity. Consequently, use of the IR luminosity as an SFR indicator may cause one to conclude that post-starburst galaxies are still star-forming, whereas in reality, star formation was recently quenched.

\section*{Acknowledgments}\end{small}}

\voffset=-0.6in

\title[IR luminosity as an SFR indicator]{The total infrared luminosity may significantly overestimate
the star formation rate of quenching and recently quenched galaxies}
\author[C.~C. Hayward et al.]{
\parbox[t]{\textwidth}{
Christopher C. Hayward$^{1,2}$\thanks{E-mail: christopher.hayward@h-its.org}\thanks{Moore Prize Postdoctoral Scholar in Theoretical Astrophysics},
Lauranne Lanz$^3$, Matthew L. N. Ashby$^4$, Giovanni Fazio$^4$, Lars Hernquist$^4$, Juan Rafael Mart\'{i}nez-Galarza$^4$,
Kai Noeske$^5$, Howard A. Smith$^4$, \\ Stijn Wuyts$^6$ and Andreas Zezas$^{4,7,8}$}
\vspace*{6pt} \\
$^1$TAPIR, Mailcode 350-17, California Institute of Technology, 1200 E. California Boulevard, Pasadena, CA 91125, USA \\
$^2$Heidelberger Institut f\"ur Theoretische Studien, Schloss--Wolfsbrunnenweg 35, 69118 Heidelberg, Germany \\
$^3$Infrared Processing and Analysis Center, California Institute of Technology, 1200 E. California Boulevard, Pasadena, CA 91125, USA \\
$^4$Harvard--Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA \\
$^5$Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA\\
$^6$Max-Planck-Institut f\"ur extraterrestrische Physik, Postfach 1312, Giessenbachstrasse 1, 85741 Garching, Germany \\
$^7$University of Crete, Physics Department \& Institute of Theoretical \&
Computational Physics, 71003 Heraklion, Crete, Greece \\
$^8$Foundation for Research and Technology-Hellas, 71110 Heraklion, Crete, Greece}

\usepackage{hyperref}

\begin{document}

\date{Submitted to MNRAS}

\pagerange{\pageref{firstpage}--\pageref{lastpage}} \pubyear{2014}

\maketitle

\label{firstpage}

\begin{abstract}
The total infrared (IR) luminosity is very useful for estimating the star formation rate (SFR) of galaxies,
but converting the IR luminosity into an SFR relies on assumptions that do not hold for all galaxies. We
test the effectiveness of the IR luminosity as an SFR indicator by applying it to synthetic spectral energy
distributions generated from three-dimensional hydrodynamical simulations of isolated disc galaxies and
galaxy mergers. In general, the SFR inferred from the IR luminosity agrees well with the true instantaneous
SFR of the simulated galaxies. However, for the major mergers in which a strong starburst is induced,
the SFR inferred from the IR luminosity can overestimate the instantaneous SFR during the post-starburst
phase by greater than two orders of magnitude. Even though the instantaneous SFR decreases rapidly
after the starburst, the stars that were formed in the starburst can remain dust-obscured and thus produce
significant IR luminosity. Consequently, use of the IR luminosity as an SFR indicator may cause one to
conclude that post-starburst galaxies are still star-forming, whereas in reality, star formation was recently
quenched.
\end{abstract}

\begin{keywords}
dust, extinction -- galaxies: interactions -- galaxies: starburst -- infrared: galaxies -- radiative transfer --
stars: formation.
\end{keywords}

\section{Introduction} \label{S:intro}

The total infrared (IR) luminosity ($L_{\rm IR}$) of galaxies has been used to infer the star formation rate (SFR)
of galaxies for decades (see \citealt{Kennicutt:1998review} and \citealt{Kennicutt:2012} for reviews). Since
the launch of the
\textit{Herschel Space Observatory}, the IR luminosity has been often used to infer the SFRs of diverse
types of local and high-redshift galaxies \citep[e.g.][]{Rodighiero:2010,Elbaz:2011,Wuyts:2011a,
Magnelli:2012,Rosario:2012,Rosario:2013,Lee:2013}.

Using the IR luminosity as an SFR indicator has several advantages. For example, other SFR
indicators, such as recombination lines (e.g. H$\alpha$\xspace) and ultraviolet (UV) emission, must be corrected
for dust attenuation to recover the intrinsic SFR. Correcting for dust is inherently uncertain because the
observed emission is
dominated by the less-obscured lines of sight; consequently, the dust attenuation can be underestimated.
In contrast, the IR is a promising SFR indicator precisely because of dust: in the limit of high obscuration,
essentially all of the UV-optical light from young stars is absorbed and re-emitted into the IR. Thus, the IR
luminosity should be the best SFR indicator in those situations, such as for starburst galaxies. Furthermore,
the IR luminosity can be combined
with tracers of unobscured star formation (uncorrected for dust attenuation); such combinations may
yield more accurate SFR estimates (e.g. \citealt{Kennicutt:2007,Kennicutt:2009}; \citealt*{Relano:2009};
\citealt{Wuyts:2011a,Wuyts:2011b,Reddy:2012, Lanz:2013}).

Naturally, there are also caveats to using the IR luminosity as an SFR indicator because multiple
assumptions must be made to convert the IR luminosity into an SFR (see \citealt{Kennicutt:1998review}
for a more thorough discussion). As with all SFR indicators,
the IR luminosity is sensitive to the initial mass function (IMF; e.g. \citealt{Narayanan:2012IMF}). The
sensitivity to the high end of the IMF is less severe for the
IR than for recombination lines because the IR light originates from dust that is heated by both lower-mass
stars and the massive stars that photoionise nebulae and thereby generate recombination line emission.
However, this effect can also be a disadvantage of using the IR luminosity as an SFR indicator: it is
commonly stated \citep[e.g.][]{Kennicutt:1998review} that, similar to the UV, the IR luminosity is sensitive
to the SFR averaged over the past $\sim 100$ Myr (in contrast, recombination lines are sensitive to the
SFR averaged over the past $\sim 10$ Myr). Thus, the IR
luminosity may not trace the instantaneous SFR accurately if it varies over timescales of $\la 100$ Myr.

Furthermore, dust can also be heated by older stars; if such dust heating is non-negligible, as has
been suggested by some authors (e.g. \citealt*{Sauvage:1992}; \citealt{Smith:1994};
\citealt*{Smith:1996,Walterbos:1996};
\citealt{Kennicutt:2009,Salim:2009}; \citealt{Calzetti:2010};
\citealt{Kelson:2010,Murphy:2011,Totani:2011}; \citealt{Groves:2012,
Leroy:2012}; \citealt{Fumagalli:2013,Utomo2014})
converting the IR luminosity into an SFR using a standard calibration will overestimate the true SFR.
This effect becomes increasingly important with decreasing SFR \citep{Calzetti:2010,Smith:2012,
Rowlands2014,Utomo2014}.
When the contribution of older stellar populations to dust heating is non-negligible, the implicit
averaging timescale is $>> 100$ Myr \citep{Salim:2009,Kelson:2010}.
Emission from active galactic nuclei (AGN) can also
{heat dust and thus cause the IR-inferred SFR to overestimate the true SFR.
The IR-inferred SFR can also under-predict the true SFR
if a significant fraction of the photons from young stars is not absorbed by dust or if stellar
populations younger than 100 Myr dominate the bolometric luminosity \citep[e.g.][]{Schaerer:2013,
Sklias2014}.

In this paper, we quantify the effectiveness of the IR luminosity as an SFR indicator by applying it to mock
spectral energy distributions (SEDs) calculated from three-dimensional hydrodynamical simulations of
isolated disc galaxies and galaxy mergers. Because the simulated galaxies have star formation histories
that are significantly more complex than those assumed to calibrate the $L_{\rm IR}$-SFR relation (i.e., a single
burst or continuous star formation history), realistically treat dust attenuation from both
stellar birth clouds and the diffuse interstellar medium (ISM), and include AGN emission, they can be used
to characterise the effects of some of the aforementioned assumptions that may cause the IR luminosity to
poorly trace the instantaneous SFR. The remainder of this work is organised as follows: in Section
\ref{S:methods}, we briefly summarise the details of the simulations. In Section \ref{S:results}, we compare
the $L_{\rm IR}$-inferred SFR and true instantaneous SFR and the resulting inferred and true SFR--stellar mass
($M_{\star}$) relations. In Section \ref{S:discussion}, we discuss some implications of our results and ways
forward.

\section{Methods} \label{S:methods}

The simulations used in this work are a subset of those used in
\citet[][hereafter L14]{Lanz2014}. We briefly describe them here, but we refer
the reader to that work for full details. We performed three-dimensional
hydrodynamical simulations of four isolated disc galaxies and
binary mergers of all possible combinations of those galaxies for a single
(non-special) orbit (thereby yielding ten merger simulations)
using the $N$-body/smoothed-particle hydrodynamics (SPH) code \textsc{gadget-3}\xspace
\citep{Springel:2005gadget}.\footnote{Recently,
various authors have highlighted issues with the standard formulation of SPH
that may cause the results of simulations performed using the
technique to be inaccurate. However, for the type of idealised simulations
used in this work, the standard form of SPH yields results that are very
similar to those of the more-accurate moving-mesh hydrodynamics technique
\citep{Hayward2014arepo}.} The progenitor disc galaxies are referred to
as M0, M1, M2, and M3, where a higher number indicates a greater mass;
the masses of the stellar discs of (M0, M1, M2, M3) are $(0.061, 0.38, 1.18, 4.22) \times 10^{10} ~\mathrm{M_{\odot}}$,
and other galaxy properties, such as the gas fraction and disc scale length, are varied
such that the galaxies represent typical local galaxies \citep{Jonsson:2006,Cox:2008}.
The major (equal-mass) merger simulations used in this work are referred to as
M0M0e, M1M1e, M2M2e, and M3M3e, which correspond to mergers of two
M0, M1, M2, and M3 disc galaxies, respectively, on the \textit{e} orbit of \citet{Cox:2006}.

The simulations include gas cooling and simple models for star formation, the multiphase ISM
\citep{Springel:2003}, and black hole accretion and AGN feedback
\citep{Springel:2005feedback}. The star formation rate of the simulations is calculated
on a particle-by-particle basis: each gas particle with density greater than
a threshold density ($n_{\rm H} \sim 0.1$ cm$^{-3}$) is assigned an SFR according
to a volume-density-dependent Kennicutt--Schmidt
prescription \citep{Schmidt:1959,Kennicutt:1998}, $\rho_{\rm SFR} \propto \rho^N$;
we assume $N = 1.5$. Summing the SFRs of the individual gas particles yields
the total instantaneous SFR of the galaxy, SFR$_{\mathrm{inst}}$. In practice, star formation is
implemented by stochastically converting gas particles into star particles,
where the probability of a gas particle being converted into a star particle
depends on its SFR.

\begin{figure}
\centering
\includegraphics[width=\columnwidth]{./figures/sfr_comp_iso.eps}
\caption{SFR inferred from the IR luminosity vs. the instantaneous
SFR for simulated isolated disc galaxies. The points are coloured according to the
simulation from which they originate. The black line indicates equality. The inferred and actual SFR
typically agree well. For the lower-mass galaxies, a significant fraction of the photons emitted
by young stars is not absorbed; consequently, the SFR is slightly underestimated.}
\label{fig:sfr_comp_iso}
\end{figure}

To calculate mock SEDs of the simulated
galaxies, we perform three-dimensional dust radiative transfer in
post-processing at multiple times throughout each \textsc{gadget-3}\xspace simulation
using the Monte Carlo dust radiative transfer code \textsc{sunrise}\xspace
\citep*{Jonsson:2006sunrise,Jonsson:2010sunrise}.
The star and black hole particles in the \textsc{gadget-3}\xspace simulations are
the sources of radiation, and the metal distribution determines the
dust distribution; we assume that 40 per cent of the metals are in
dust \citep[e.g.][]{James:2002,Sparre2014b}.
We use the default ISM model of L14, which is a conservative
choice because it exhibits less dust attenuation than the alternative model,
and we have checked that our conclusions are insensitive to whether we use the default
or alternative ISM treatments.
Given these inputs, radiative transfer is performed to calculate the
absorption, scattering, and re-emission of stellar and AGN light by dust.

These or similar simulations have been demonstrated to provide good
matches to the SEDs of local disc galaxies \citep{Jonsson:2010sunrise},
local interacting galaxies (L14), and high-redshift obscured
starbursts and AGN (\citealt{Narayanan:2009,Narayanan:2010dog,Narayanan:2010smg,Snyder:2013};
E. Roebuck et al., in preparation). Furthermore, at
different times in their evolution, they reproduce various properties of
other populations, including submillimetre galaxies
\citep[e.g.][]{Hayward:2011smg_selection,
Hayward:2012smg_bimodality,Hayward:2013number_counts},
post-starburst galaxies \citep{Snyder:2011}, and compact quiescent galaxies
\citep{Wuyts:2010}, among others. The success of the models supports our use
of them to test how well $L_{\rm IR}$ traces the instantaneous SFR.

\section{Results} \label{S:results}

To infer the SFR from the IR luminosity of our simulated galaxies, we first
integrate the (spatially integrated) galaxy SEDs over the wavelength range $8-1000$
\micron~to calculate $L_{\rm IR}$. We adopt this wavelength range following
\citet{Kennicutt:1998review}, but the precise definition (e.g. if we use only the far-IR, FIR)
is unimportant for our conclusions because the FIR luminosity dominates $L_{\rm IR}$.
Then, we convert $L_{\rm IR}$ to an SFR using the following relation:
\begin{equation} \label{eq:sfr_cal}
{\rm SFR}_{\rm IR} = 3.0 \times 10^{-37} \left(\frac{L_{\rm IR}}{\mathrm{W}}\right) ~\mathrm{M_{\odot} {\rm ~yr}^{-1}},
\end{equation}
which is the \citet{Kennicutt:1998review} conversion factor converted to the \citet{Kroupa:2001}
IMF by dividing by 1.5 \citep[e.g.][]{Schiminovich:2007} .
We have made this conversion because the \citeauthor{Kroupa:2001} IMF is assumed
when calculating the {\sc starburst99} \citep{Leitherer:1999} single-age stellar population SEDs
that are used as input for the radiative transfer calculations; thus, we must use the same IMF for
the calibration for consistency.
Also, we note that this factor is calculated by assuming continuous bursts of age 10-100 Myr, and
\citet{Kennicutt:1998review} explicitly states that strictly speaking, it is valid only for
starbursts that are less than 100 Myr old. Nevertheless, this conversion is routinely
applied to diverse galaxy types, including those for which detailed information about their SFHs
is unavailable (e.g. for high-redshift galaxies).

\begin{figure}
\centering
\includegraphics[width=\columnwidth]{./figures/sfr_comp_mergers.eps} \\
\includegraphics[width=\columnwidth]{./figures/ssfr_comp_mergers.eps}
\caption{SFR (top) and SSFR (bottom) inferred from the IR luminosity vs. the instantaneous
SFR (SSFR) for the major merger simulations.
The data points are connected to illustrate the time evolution. Near coalescence
of the black holes, which is indicated by a star, the data points are separated by 10 Myr.
At other times, the separation is 100 Myr.
The large filled symbols indicate the start of each simulation. In the simulations in which a strong starburst
is induced, the IR-inferred SFR (SSFR) can overestimate the instantaneous SFR (SSFR)
in the post-starburst phase by greater than two orders of magnitude.
Note that for fixed SSFR (of $\ga 10^{-11}$ yr$^{-1}$), the overestimation occurs only in the
post-starburst phase of the simulations. This indicates that the time evolution of the SSFR, not
just the value of the instantaneous SSFR, determines whether the IR luminosity overestimates the SSFR.}
\label{fig:sfr_comp_mergers}
\end{figure}

\subsection{Isolated disc galaxies}

Fig. \ref{fig:sfr_comp_iso} compares the SFR inferred from the IR luminosity
using Equation (\ref{eq:sfr_cal}) and the instantaneous
SFR\footnote{The `instantaneous SFR' of the simulated galaxies is
the sum of the SFRs of the individual gas particles, which are calculated based on their gas
densities and the assumed sub-resolution star formation prescription. Consequently, the instantaneous
SFR value for the simulations corresponds to an average over a timescale that is less than the maximum
time step, 5 Myr. The exact timescale is not well-defined because adaptive time steps are used.
However, the salient feature of the SFR value of the simulations is that the averaging timescale is
less than that of any of the observable SFR tracers.}
for simulated isolated disc galaxies. The data are coloured according to the
simulation from which they originate. For the major mergers, the data points are connected to
illustrate the time evolution. The start of each simulation is marked with a larger
filled symbol, and the time of coalescence of the central black holes is marked with a star.
In this and the following figures, we show the results for a single viewing angle because
the SFR is independent of viewing angle by definition, and
the variation in $L_{\rm IR}$ with viewing angle is $\la 10$ percent
(we note that the variation with viewing angle can be more significant in high-redshift, gas-rich mergers;
\citealt{Hayward:2012smg_bimodality}).

The inferred and true SFRs agree well for the isolated discs, although the
IR luminosity tends to underestimate the SFR of the least-massive disc galaxy
(M0). The reason for the underestimation is that this galaxy is not sufficiently
dust-obscured. Thus, a significant fraction of the light from young stars
escapes the galaxy without being absorbed by dust.

\subsection{Galaxy mergers}

Fig. \ref{fig:sfr_comp_mergers} compares the IR-inferred and true instantaneous SFR (top)
and specific SFR (SSFR $\equiv \mathrm{SFR}/M_{\star}$; bottom) values
for the major merger simulations. The inferred and actual (S)SFR values agree very well
for much of the mergers' evolution. However, there is an important exception: in the
simulations in which a strong starburst is induced near coalescence (all but M0M0e),
the IR-inferred (S)SFR significantly overestimates the
instantaneous (S)SFR in the post-starburst phase. This overestimation first occurs
immediately after the maximum SFR of the starburst is attained, and it can persist
for hundreds of Myr (until the end of the simulations).
At the lowest (S)SFRs, the instantaneous (S)SFR can be overestimated by greater than
two orders of magnitude.

One might expect that the SSFR determines whether older stellar populations contribute
significantly to the IR luminosity: for actively star-forming galaxies (i.e. those with high SSFR values), young
stars should dominate the dust heating, whereas for galaxies with low SSFR values, dust heating from older
stellar populations may be significant. Indeed, the bottom panel of Fig. \ref{fig:sfr_comp_mergers} indicates
that the overestimation is more significant for lower SSFR values. However, this panel also indicates that
the instantaneous value of the SSFR is not the sole determinant of how well the IR luminosity can be used
to recover the (S)SFR: for fixed SSFR (of $\ga 10^{-11}$ yr$^{-1}$, which is the minimum SSFR
value during the pre-starburst phase of the simulations), the overestimation occurs only during
the post-starburst phase of the simulations. Thus, the time evolution of the SSFR, not just the
instantaneous value, is important for determining how well the IR luminosity traces the (S)SFR.

\begin{figure}
\centering
\includegraphics[width=\columnwidth]{./figures/M3e_sfr_vs_t.eps}
\caption{SFR (solid blue line) and $L_{\rm IR}$ (dashed red line) vs. time
relative to coalescence for the most-massive major merger simulation (M3M3e).
The values on the left (SFR) and right ($L_{\rm IR}$) axes are related
through Eq. (\ref{eq:sfr_cal}); thus, the IR-inferred SFR can be read
off from the $L_{\rm IR}$ curve by referring to the left axis.
Until the peak of the starburst, $L_{\rm IR}$ has a similar
evolution to the SFR, and thus the IR-inferred SFR agrees with the true SFR.
However, after the starburst, $L_{\rm IR}$ decreases much
more gradually than the SFR because there is significant dust heating from
stars that formed during the burst; thus, $L_{\rm IR}$ is a poor tracer of the SFR in this regime.}
\label{fig:sfr_vs_t}
\end{figure}

To clarify the reason for the overestimation, we show the time evolution of the SFR and $L_{\rm IR}$ for
one of the major merger simulations, M3M3e (the most-massive simulation in our suite), in Fig.
\ref{fig:sfr_vs_t}. Near the time of final coalescence of the two galaxies, a strong starburst is induced
because tidal forces drive gas into the nuclear region
and consequently, the gas density and thus SFR increase sharply. In the starburst,
the instantaneous SFR increases by a factor of $\sim 30$ in $<200$ Myr.
Subsequently, because most of the
gas has been consumed by star formation or heated by shocks and AGN feedback
(i.e., the star formation is `quenched'), the SFR plummets by greater than four orders of magnitude.
Afterward, the SFR remains very low, and the galaxy is a passively evolving spheroid.

Until the peak of the starburst, the evolution of the SFR and $L_{\rm IR}$
is essentially identical. However, after the peak of the starburst, when the SFR
decreases very rapidly, this similarity no longer holds. Rather, $L_{\rm IR}$
declines more gradually; this difference is the origin of the aforementioned
discrepancy between the SFR inferred from $L_{\rm IR}$ and the instantaneous SFR.
Part of the reason for the difference is that during the starburst, the SFR varies
significantly on timescales of less than 100 Myr. Consequently, the instantaneous SFR
and SFR averaged over the past 100 Myr differ significantly at this time. Thus, one of the fundamental
assumptions of the $L_{\rm IR}$-SFR conversion does not hold. However, the overestimation persists and
actually becomes more severe well after the starburst, when the instantaneous SFR is
relatively constant on a 100-Myr timescale.

The fundamental reason for the overestimation is that although the SFR is extremely low
($\sim 10^{-3} ~\mathrm{M_{\odot} {\rm ~yr}^{-1}}$) after the starburst, a significant number
of young stars (with ages of 100s of Myr) remain. These stars are still
very luminous and obscured and hence can effectively heat the dust and thus
yield significant IR emission. Consequently, the effective averaging timescale for the IR-inferred
SFR is $>>100$ Myr during the post-starburst phase of the simulations.

In principle, the AGN can also
heat the dust and result in IR emission. By re-running the radiative
transfer calculations with the AGN disabled, we have confirmed that
dust heated by AGN is not the dominant cause of the overestimation (in all but a few snapshots of
these simulations, $L_{\rm IR}$ decreases by $\la 50$ per cent when the AGN emission is disabled). Nevertheless,
FIR emission from AGN-heated dust can be important in some regimes and
will be a subject of future work (L. Rosenthal et al., in preparation).

\subsection{Implications for the SFR-stellar mass relation}

\begin{figure}
\centering
\includegraphics[width=\columnwidth]{./figures/sfr-mstar.eps}
\caption{SFR vs. $M_{\star}$ for the major merger simulations.
The blue circles indicate that the true instantaneous
SFR is used, whereas the red triangles indicate that the SFR value is inferred
from $L_{\rm IR}$. In both cases, the true $M_{\star}$ is used (but for clarity, the
$M_{\star}$ values for the red triangles have been increased by 0.03 dex).
The solid black line is the best fit to the $z \sim 0.2-0.4$ star-forming galaxies of Karim et al. (2011).
For the recently quenched galaxies, the inferred and actual locations in the
SFR--$M_{\star}$ plane differ significantly (i.e., there are many blue
circles that lie significantly below all of the red triangles).}
\label{fig:sfr-mstar}
\end{figure}

The overestimation demonstrated above has significant consequences for interpreting
IR observations (from e.g. \textit{Herschel} and the Atacama Large Millimeter/submillimeter Array,
ALMA) of galaxies that have
recently undergone strong starbursts. One topic for which this overestimation is potentially
significant is the SFR--$M_{\star}$ relation (aka `main sequence') of star-forming
galaxies \citep{Brinchmann:2004,Salim:2005,Noeske:2007a,Daddi:2007}.
Fig.~\ref{fig:sfr-mstar} shows SFR vs. $M_{\star}$ for the major merger simulations
for both the true instantaneous SFR and the SFR inferred from $L_{\rm IR}$.
Because the galaxies have low initial gas fractions (see L14 for details), the stellar mass
increases by a relatively small amount over the course of a simulation; consequently,
a given simulation primarily moves up and then down in the SFR--$M_{\star}$ plot as
the SFR is increased during the starburst and then quenched.

In most cases, the inferred and actual positions of
the simulated galaxies in the SFR--$M_{\star}$ plot agree because the IR-inferred SFR
and instantaneous SFR values are similar and the true
$M_{\star}$ values are used for both.\footnote{Uncertainties in
how the stellar mass is inferred from observations can also have important
consequences for the resulting SFR--$M_{\star}$ relation \citep[e.g.][]{Michalowski:2012}.}
However, for the post-starburst phase of the
strong-starburst mergers, this is not the case; instead, the inferred and
true locations in the SFR--$M_{\star}$ plot differ significantly (i.e., there are many blue
circles that lie significantly below all of the red triangles). Because the
IR-inferred SFR overestimates the true SFR in this phase, use of the IR-inferred
SFR leads to the conclusion that the recently quenched galaxies are much
closer to the SFR--$M_{\star}$ relation of actively star-forming galaxies (the black line in Fig. \ref{fig:sfr-mstar}) than they
actually are. This effect may also cause the apparent scatter in the SFR--$M_{\star}$ relation to be less
then the true scatter (see also \citealt{Schaerer:2013}),
but cosmological simulations are necessary to quantify the
effect on the scatter, which depends on the abundance of sufficiently strong
starbursts and the duration of the overestimate. Unfortunately, even state-of-the-art
cosmological simulations are not yet suitable for such an analysis
because of their relatively limited resolution \citep{Sparre2014}.

\section{Discussion} \label{S:discussion}

\subsection{Implications and scope of the overestimation}

Because of the SFR overestimation demonstrated above, when one uses $L_{\rm IR}$ to infer the SFR of galaxies
that have been recently quenched, inaccurate conclusions are likely.
For example, as we have demonstrated, the IR-inferred SFR may suggest
that a galaxy lies closer to the SFR--$M_{\star}$ relation than it actually does.
Similarly, FIR emission from AGN host galaxies
is often considered evidence that AGN do not quench star formation, which would
contradict many theoretical models \citep[e.g.][]{Springel:2005c,Hopkins:2006unified_model}.
However, the results presented above suggest that if
quenching happens rapidly (on timescales of order 100 Myr), as is suggested by
simulations, detection of FIR emission from AGN hosts does not (necessarily)
imply that star formation is ongoing in those galaxies. Other topics for
which this overestimation can have significant implications include the FIR--radio
correlation (see \citealt{Condon:1992} for a review) and the Kennicutt--Schmidt relation.

Because of the typical depths of FIR imaging, the potential bias during the post-starburst
phase that we have demonstrated is most serious for observations of low-redshift galaxies and stacking analyses
of high-redshift galaxies; i.e., for many post-starburst galaxies, SFR$_\mathrm{IR}$ may overestimate
the instantaneous SFR, but $L_{\rm IR}$ may be less than the detection threshold and thus
render the overestimate irrelevant.
We stress that for SFR (SSFR) values of as high as $\sim 10 ~\mathrm{M_{\odot} {\rm ~yr}^{-1}}$
($\sim 10^{-10}$ yr$^{-1}$), the instantaneous SFR can be overestimated by factors
of a few or more. This is consistent with the results of \citet{Utomo2014}, who presented observational evidence
that indicates that the SFRs inferred from a combination of UV and IR luminosities overestimate
the true values for the same (S)SFR range as in the simulations. Thus, this overestimation is
clearly relevant even for existing observational surveys and will become even more important as detection
limits are pushed lower using e.g. ALMA.

Furthermore, it is important to note that the overestimation problem is
not limited to the specific situation studied here (low-redshift galaxy mergers); rather,
whenever the SFR declines sufficiently rapidly, regardless of the absolute SFR values
and the quenching mechanism, this overestimation should occur. Indeed, this overestimation
is also present in the post-starburst phase of simulations of $z \sim 2-3$ galaxy mergers
(see fig. 1 of \citealt{Hayward:2011smg_selection}).

\subsection{Possible solutions}

One method to correct IR-inferred SFRs for contamination from older stellar populations is to
subtract the diffuse component of the dust emission, which is likely heated predominantly by the diffuse
interstellar radiation field and thus connected with older stellar populations rather than recently formed stars
\citep[e.g.][]{Smith:1994,Smith:1996,Groves:2012,Leroy:2012,LoFaro:2013}.
However, such a correction requires spatially resolved, multi-wavelength IR data, which are typically
available only for low-redshift galaxies. It is perhaps possible to calibrate corrections using local
galaxies and apply them to high-redshift galaxies, but it is unclear that such a technique would be
effective because of the significant differences between the properties of low- and high-redshift galaxies.

Another possible approach for avoiding this overestimation is to use SFR tracers other than the IR
or a `ladder of SFR indicators' \citep{Wuyts:2011a,Wuyts:2011b}.
Radio emission is particularly promising because it is insensitive to dust obscuration and contributions
from older stellar populations \citep{Murphy:2011}. However, one must be careful to account for other
sources of radio emission.
Recombination lines are another type of SFR tracer that is insensitive to contamination from older
stellar populations. However, the most commonly used recombination line, H$\alpha$, must be corrected
for dust attenuation, and our results suggest that in the post-starburst regime, the effects of dust
are particularly nefarious: suppose that the SFR inferred from the IR luminosity is greater
than that inferred from the UV emission or recombination lines. The usual interpretation
is that dust attenuation causes the UV-optical tracers to underestimate the true SFR;
thus, the IR estimate is more accurate. However, the effect that we have demonstrated
works in the opposite direction: dust attenuation of older stars can cause the IR
to overestimate the SFR. Consequently, when UV-optical and IR SFR indicators disagree, it
is unclear a priori which is closer to the truth.

Because of this degeneracy, it is best to use broadband SED modelling, a multi-wavelength
approach that accounts for both unobscured and obscured star formation and that can place
meaningful constraints on the star formation history rather than just a single
value of the SFR that corresponds to a somewhat ill-defined time-averaged SFR.
Multiple potentially suitable tools, including {\sc magphys} \citep{daCunha:2008},
{\sc cigale} \citep{Noll:2009,Serra:2011}, and {\sc chiburst}
(\citealt{Dopita:2005,Dopita:2006c,Dopita:2006b,Groves:2008,MG2011};
R. Mart\'inez-Galarza et al., in preparation), have been developed for this task.
Application of one of these tools ({\sc magphys}) to simulated SEDs suggests that the
inferred SFRs are indeed more accurate than the IR-inferred SFRs (C. Hayward \& D. Smith,
submitted), which has also been suggested based on observational evidence \citep{Smith:2012,
Rowlands2014,Utomo2014}. The aforementioned tools and similar methods may be able to overcome the
shortcomings of simple SFR indicators such as the IR luminosity and
thus fully leverage the wealth of multiwavelength data that has become available
in recent years.

\section{Conclusions}

We have used a set of mock SEDs generated by performing dust radiative transfer
on hydrodynamical simulations of isolated disc galaxies and galaxy mergers to
investigate the effectiveness of the IR luminosity at recovering the true
instantaneous SFR of the simulated galaxies. Our principal conclusions are the
following:
\begin{enumerate}
\item For most galaxies (i.e., quiescently star-forming disc galaxies and minor mergers), the SFR inferred
from $L_{\rm IR}$ agrees very well with the true SFR.
\item In the simulations in which a strong starburst occurs, $L_{\rm IR}$
can significantly overestimate the SFR after the peak of
the starburst. The reason for this overestimation is that although the SFR
decreases rapidly from the peak, there is still significant IR emission
from dust heated by the stars formed during, and even before, the starburst.
\item The magnitude of the overestimate is greater for lower SSFR values. However, the instantaneous
SSFR value is not the sole determinant of the contribution of older stellar populations to the dust heating
because for fixed SSFR, the overestimation is more severe in the post-starburst phase of the simulations.
Thus, the time evolution of the SSFR, not just the current value, is an important determinant of how well
the IR-inferred SFR traces the true SFR.
\item This overestimation may have significant implications for e.g. the SFR-stellar mass
relation. In particular, it may cause the number of quenched galaxies and degree
to which galaxies are quenched to be underestimated.
\end{enumerate}
The results presented in this work highlight the need for caution when applying simple SFR tracers
to galaxies for which the underlying assumptions may not hold. More-sophisticated techniques, such
as SED modelling, should be used when possible, but one must also mind the caveats and potential
inherent degeneracies in these techniques.

\begin{small}\section*{Acknowledgments}\end{small}

We thank Caitlin Casey, Mattia Fumagalli, Xavier Koenig, Barry Rothberg,
Daniel Schaerer, Beverly Smith, Dan Smith, and Tomo Totani for providing useful comments
on the manuscript, and we especially thank Samir Salim and the anonymous referee for their very detailed
comments, which led to significant improvements to the manuscript.
CCH is grateful to the Klaus Tschira Foundation for financial support and
acknowledges the hospitality of the Aspen Center for Physics, which is supported
by the National Science Foundation Grant No. PHY-1066293.
HAS and LL acknowledge partial support from NASA grants NNX12AI55G and
NNX10AD68G, and JPL RSA contract 1369566. The simulations in
this paper were performed on the Odyssey cluster supported by the FAS Research
Computing Group at Harvard University. This research has made use of NASA's Astrophysics Data System.
\\

\footnotesize{
\bibliographystyle{mn2e}

Astrophysics of Galaxies (astro-ph.GA); Cosmology and Nongalactic Astrophysics (astro-ph.CO)


Metals in Star-Forming Galaxies at High Redshift

Abstract: The chemical composition of high-redshift galaxies is an important property that gives clues to their past history and future evolution. Measuring abundances in distant galaxies with current techniques is often a challenge, and the canonical metallicity indicators can often not be applied. I discuss currently available metallicity indicators based on stellar and interstellar absorption and emission lines, and assess their limitations and systematic uncertainties. Recent studies suggest that star-forming galaxies at redshift around 3 have heavy-element abundances already close to solar, in agreement with predictions from cosmological models.

\section{Introduction}

Until recently, abundance determinations at cosmological redshifts ($z \gtrsim 2$) were limited to bright quasars and the material along their sightlines. Over the past 5~--~10 years, the number and nature of accessible objects has expanded dramatically, and ``normal'' star-forming galaxies at high redshift have become targets of abundance studies (Pettini 2004). This progress was made possible by the use of efficient search techniques for finding large numbers of candidate galaxies (Steidel et al. 1999) and by the availability of 8~--~10-m class telescopes for spectroscopic studies. I will discuss the general properties of these galaxies, introduce the techniques used for determining abundances, and provide an overview of their chemical composition. I will conclude with a cosmological perspective.

\section{Sample}

The galaxies discussed in this paper are commonly referred to as {\em Lyman-break galaxies} (LBG; Giavalisco 2002) because of the search technique employed to detect them. They are color-selected via their strong Lyman discontinuity and/or their blue restframe ultraviolet (UV) colors. These galaxies are sufficiently bright for obtaining spectroscopic data ($25^m \gtrsim R \gtrsim 21^m$) --- a necessary prerequisite for abundance analyses. Owing to the selection criterion, they are relatively dust-poor with typical visual attenuations of $A_{\rm V} \simeq 1^m$ and are fairly luminous, but not ultraluminous ($L \simeq 10^{11}$~L$_\odot$). LBGs are actively star-forming at rates of order $10^2$~M$_\odot$~yr$^{-1}$ (Pettini et al. 2001) and have correspondingly short recycling times. Therefore the chemical composition of the newly formed stars and the interstellar gas are identical.

Most LBGs are in the redshift range $2 \lesssim z \lesssim 3.5$. The lower and upper limits are imposed by atmospheric cut-offs in ground-based observations and decreasing efficiencies of the search techniques. The corresponding lookback times in a standard universe are 11~Gyr~$\lesssim t \lesssim$~12.5~Gyr, or about 15\% of the age of the universe after its formation.

\begin{figure}
\begin{center}
 \includegraphics[angle=90,width=0.75\textwidth]{Leitherer_Paris05_fig1.eps}
\end{center}
  \caption{Comparison between the observed spectrum of MS1512--cB58 (solid) and two synthetic models
with 1/4~Z$_\odot$ (lower; dashed) and Z$_\odot$ (upper; dotted).
The models have continuous star formation, age 100 Myr, and Salpeter IMF between 1 and 100 M$_\odot$. The stellar lines
are weaker in the metal-poor model (from Leitherer et al. 2001).}
\end{figure}

\section{Techniques --- Restframe Optical versus UV}

Abundance determinations typically fall into two categories, either relying on indicators in the restframe optical, or on those in the restframe UV. The restframe optical wavelength region has traditionally been used to determine galaxy abundances from nebular emission lines. At a redshift of $z = 3$, the restframe optical is observed in the near-infrared (IR) H and K bands. Spectroscopic observations of LBGs in the near-IR have become technically feasible (e.g., Pettini et al. 2001) but abundance analyses are still challenging. Only the strongest lines such as, e.g., H$\alpha$, H$\beta$, [N~II] $\lambda$6584, or [O~III] $\lambda$5007 are detectable at sufficient S/N. Even when good-quality spectra are available, the atmospheric windows usually restrict the wavelengths to a narrow range, which precludes commonly used techniques such as the classical R23 strong-line method (McGaugh 1991).

The need for alternative variants of the classical strong-line method led Pettini \& Pagel (2004) to readdress the usefulness of the N2 and O3N2 ratios. The former is defined as the ratio [N~II] $\lambda$6584 over H$\alpha$ and was recently discussed by Denicol\'o et al. (2002); the latter includes the oxygen line for the ratio ([O~III]~$\lambda$5007/H$\beta$)/([N~II]~$\lambda$6584/H$\alpha$) and was originally introduced by Alloin et al. (1979). After calibrating the two abundance indicators with a local H~II region sample, Pettini \& Pagel find that O3N2 and N2 predict O/H to within 0.25~dex and 0.4~dex at the 2~$\sigma$ confidence level, respectively.

The observed frame optical wavelength region corresponds to the restframe UV of LBGs. The UV contains few nebular emission lines in star-forming galaxies (Leitherer 1997) and has rarely been used for chemical composition studies in {\em local} galaxies of this type. Fig.~1 compares the UV spectrum of the LBG MS1512--cB58 with theoretical spectra (Leitherer et al. 2001). Three groups of lines can be distinguished: (i) Interstellar absorption lines, most of which are strong and heavily saturated. Only in very few cases can unsaturated absorption lines in LBGs be used for an abundance analysis. (ii) Broad stellar-wind lines with emission and blueshifted absorption. These lines are the telltales of massive OB stars whose stellar winds are metallicity dependent. (iii) Weak photospheric absorption lines which can only be seen in high-quality spectra.
Abundance studies from stellar lines in restframe UV spectra must rely either on suitable template stars or on extensive non-LTE radiation-hydrodynamic models which are only beginning to become available (Rix et al. 2004).

\section{The Chemical Composition of LBGs}

An initial, rough estimate of the heavy-element abundances can be obtained from the equivalent widths of the strong UV absorption lines. Heckman et al. (1998) pointed out the close correlation of the Si~IV $\lambda$1400 and C~IV $\lambda$1550 equivalent widths with O/H in a sample of {\em local} star-forming galaxies. This correlation seems surprising, as these stellar-wind lines are deeply saturated. The reason for the metallicity dependence is the behavior of stellar winds in different chemical environments. At lower abundance, the winds are weaker and have lower velocity, and the lines become weaker and narrower. As a result, the equivalent widths are smaller at lower O/H. If the same correlation holds at high redshift, the observed equivalent widths in LBGs suggest [O/H]~$\simeq$~--0.5 (Leitherer 1999). A similar, somewhat weaker correlation exists between O/H and the equivalent widths of the strongest {\em interstellar} lines. This is even more unexpected because the equivalent widths of saturated lines have essentially no dependence on the column density $N_{\rm ion}$: $W \propto b[\ln(N_{\rm ion}/b)]^{0.5}$. Therefore the correlation must be caused by the $b$ factor, and therefore by velocity. More metal-rich galaxies are thought to host more powerful starbursts with correspondingly larger mechanical energy release by stellar winds and supernovae. The energy input leads to increased macroscopic turbulence and higher gas velocities at higher O/H (Heckman et al. 1998). If the same applies to star-forming galaxies in the high-redshift universe, their measured equivalent widths again indicate an oxygen abundance of about 1/3 the solar value.

Pettini et al. (2001) determined oxygen abundances in five LBGs from emission lines in restframe optical spectra. The redshift range of the sample dictated the use of the R23 method. The galaxies turned out to be rather metal-rich, with O/H somewhat below the solar value. This is roughly in agreement with restframe UV results, and an order of magnitude above the metallicities found in damped Lyman-$\alpha$ absorbers (DLA) which are found at the same redshift. Because of the double-valued nature of the R23 method, the possibility exists but is deemed less likely that the sample has oxygen abundances of only 1/10 the solar value.

The lensed LBG MS1512--cB58 and its bright restframe UV spectrum can be studied at sufficiently high S/N and resolution to detect and resolve faint, unsaturated interstellar absorption lines. Pettini et al. (2002) measured numerous transitions from H to Zn covering several ionization stages. Abundances of several key elements could be derived. The $\alpha$-elements O, Mg, Si, P, and S all have abundances of about 40\% solar, indicating that the interstellar medium  is highly enriched in the chemical elements produced by type II supernovae. In contrast, N and the Fe-peak elements Mn, Fe, and Ni are all less abundant than expected by factors of several. In standard chemical evolution models, most of the nitrogen is produced by intermediate-mass stars, whereas type Ia supernovae contribute most of the Fe-peak elements. Since the evolutionary time scales of intermediate- and low-mass stars are significantly longer than those of massive stars producing the $\alpha$-elements, the release of N and the Fe-group elements into the interstellar medium is delayed by $\sim$10$^9$~yr. MS1512--cB58 may be an example of a star-forming galaxy in its early stage of chemical enrichment, consistent with its cosmological age of only about 15\% of the age of the universe. Mehlert et al. (2002) provided similar arguments to explain variation of the  C~IV $\lambda$1550 line relative to Si~IV $\lambda$1400 in a small sample of LBGs. C~IV appears to decrease in strength relative to Si~IV from lower to higher redshift, which may reflect the time delay of the carbon release by intermediate-mass stars.

\begin{figure}
\includegraphics[width=0.5\textwidth]{Leitherer_Paris05_fig2a.eps}\includegraphics[width=0.5\textwidth]{Leitherer_Paris05_fig2b.eps}
  \caption{Left pair of panels: comparison of the observed spectrum of MS1512--cB58 (thick) with fully synthetic spectra (thin) for five different metallicities, from twice solar to 1/20 solar. First panel: region around 1425~\AA; second: region of the Fe III blend near 1978~\AA. Each pair of panels is labeled with the metallicity of the synthetic spectrum shown. Right pair of panels: same as left pair, but for Q1307--BM1163 (from Rix et al. 2004).}
\end{figure}

The interstellar lines in LBGs have blueshifts with velocities of up to several hundred km~s$^{-1}$ indicating large-scale outflows. The associated galactic mass-loss rates of $\sim$10$^2$~M$_\odot$~yr$^{-1}$ are comparable to the rates of star formation. The newly formed heavy elements are removed from their birth sites by stellar winds and supernovae and are
transported into the halo and possibly into the intergalactic medium (Pettini et al. 2002).

Detailed studies of weak interstellar lines such as that done for MS1512--cB58 remain a technical challenge, even for high-throughput spectrographs at the largest telescopes. Furthermore, the results for Fe-peak elements carry some uncertainty because of the a priori unknown depletion corrections. Abundance analyses using stellar lines are not affected by depletion uncertainties. However, the existence of non-standard element ratios precludes the use of locally observed template spectra for spectral synthesis. Therefore our group (F. Bresolin, R. Kudritzki, C. Leitherer, M. Pettini, S. Rix) has embarked on a project to model the spectra of hot stars and link them with a spectral synthesis code to predict the emergent UV spectrum of a composite stellar population as a function of metallicity. We generated a grid of hydrodynamic non-LTE atmospheres with the WM-basic code (Pauldrach et al. 2001) and calculated the corresponding UV line spectra. The resulting library was incorporated into the Starburst99 code (Leitherer et al. 1999) which then allowed us to compute a suite of model spectra for appropriate stellar population parameters. As a first application, we used several faint stellar blends around 1425~\AA\ and 1978~\AA\ as a metallicity indicator (Rix et al. 2004). The 1425~\AA\ feature is a blend of Si~III, C~III, and Fe~V, and the 1978~\AA\ absorption is mainly Fe~III. The synthesized spectra for five metallicities are compared to the observed restframe UV spectra of MS1512--cB58 and Q1307--BM1163 in Fig.~2. The model having 40\% solar metal abundance provides the best fit to the data, in agreement with the results from other methods.

A variety of independent techniques lead to consistent results for the chemical composition of LBGs. While each method by itself is subject to non-negligible uncertainties, the overall agreement of the results gives confidence in the derived abundances. LBGs at $z \simeq 3$ have heavy-element abundances of about 1/3 the solar value.

\section{Cosmological Perspective}

\begin{figure}
\begin{center}
 \includegraphics[width=0.65\textwidth]{Leitherer_Paris05_fig3.eps}
\end{center}
  \caption{Metallicity-luminosity relationship. Data for {\em local} spiral and irregular galaxies are from Garnett (2002). The $z = 2$ objects are overluminous for their (O/H) abundances, derived using the N2 calibration of Pettini \& Pagel (2004) but lie closer to the relationship for the local galaxies than $z=3$ LBGs (from Shapley et al. 2004).}
\end{figure}

Star-forming galaxies at $z \simeq 3$, at an epoch when the universe's age was only 15\% the present value, display a high level of chemical enrichment. What does their chemical composition tell us about their relation to other galaxies at lower redshift and to other structures found at $z = 3$?

Galaxies at somewhat lower redshift have only recently become accessible for detailed study due to the combined challenges of instrumentation and the galactic spectral properties. Shapley et al. (2004) obtained K-band spectroscopy of seven UV-selected star-forming galaxies at redshifts between 2 and 2.5. The N2 method calibrated by Pettini \& Pagel (2004) was used as an abundance diagnostic. When compared to the original higher-$z$ LBGs, the $z \simeq 2$ sample is more metal-rich. This can be seen in Fig.~3, where O/H of the $z = 2$ galaxies is compared with that of LBGs at $z \simeq 3$ and of local star-forming galaxies over a range of blue luminosities. The latter were analyzed with the R23 method. The $z = 2$ sample has almost solar chemical composition but is still less metal-rich than local late-type galaxies {\em with comparable luminosities}. As a caveat, the comparison rests on the assumption that the N2 and R23 calibrations have no significant offset. The difference between the average redshift of the LBG sample and of the $z = 2$ galaxies translates into a mean age difference of about 1~Gyr. Both the chemical properties and the masses of the $z = 2$ galaxies and LBGs are consistent with standard passive evolution models.

Kewley \& Kobulnicky (2005) followed the metallicity evolution of star-forming galaxies with comparable luminosities from $z = 0$ to 3.5. O/H was determined from restframe optical emission lines using the strong-line method in four homogeneous galaxy samples. The samples were taken from the CfA2 survey, from the GOODS field, from Shapley et al. (2004), and from the LBG sample, covering $z \approx 0$, 0.7, $2.1 - 2.5$, and $2.5 - 3.5$, respectively. The average oxygen abundance in the local universe, as defined by the CfA2 sample is about solar. O/H decreases with redshift to approximately 1/3 solar at $z=3$.

It is instructive to compare the heavy-element abundances of LBGs to those of DLAs and to the Lyman-forest at the same redshift (Pettini 2004). DLA systems have metallicities of about 1/15~Z$_\odot$ and are thought to be the cross sections of the outer regions and halos of (proto)-galaxies seen along the sightlines of quasars. Although the properties of LBGs and DLAs do not immediately support a close relation between the two classes of objects, at least some link seems likely. If so, the observed outflows in LBGs may provide the metal enrichment of the halos.
The Lyman-forest is predicted by cold dark matter models to result from structure formation in the
presence of an ionizing background. The Lyman-forest had long been
thought to be truly primordial, but metal enrichment of 1/100~--1/1000~Z$_\odot$ has recently been detected (Aguirre et al. 2004).
This relatively high metal abundance early in the evolution of the universe could have been produced by a first generation of Population III stars. Such stars can account for
the amount of metals, and at the same time could have provided copious ionizing photons, as metal and photon
production are closely correlated. Alternatively, star-forming galaxies at high redshift could be the
production sites of the metals seen in the intergalactic medium {\em if} superwinds are capable of removing the newly formed metals from galactic disks.

\begin{acknowledgments}
I would like to thank Max Pettini for a careful reading of the manuscript and for providing invaluable comments.
\end{acknowledgments}


Detecting the Most Distant (z>7) Objects with ALMA

Abstract: Detecting and studying objects at the highest redshifts, out to the end of Cosmic Reionization at z>7, is clearly a key science goal of ALMA. ALMA will in principle be able to detect objects in this redshift range both from high-J (J>7) CO transitions and emission from ionized carbon, [CII], which is one of the main cooling lines of the ISM. ALMA will even be able to resolve this emission for individual targets, which will be one of the few ways to determine dynamical masses for systems in the Epoch of Reionization. We discuss some of the current problems regarding the detection and characterization of objects at high redshifts and how ALMA will eliminate most (but not all) of them.

\section{Introduction: The highest redshift galaxies}
\vspace*{-2.5mm}

In recent years, deep narrow band surveys have revealed a major
population of Lyman Alpha Emitters (LAE) out to very high redshifts
(e.g.  Hu et al.\ 2002, Kurk et al.\ 2004, Stern et al.\ 2005,
Murayama et al.\ 2007). In particular, Taniguchi et al.\ (2005) report
the detection of 9 spectroscopically confirmed LAE at redshifts of
z$\sim$6.6 in the Subaru Deep Field (currently, the published LAE
redshift record holder is at z=6.98, Iye et al.\ 2006).  The mere
presence of Lyman alpha emission in these sources provides strong
evidence that they are undergoing bursts of star formation: the star
formation rates of individual objects are
$\sim$10\,M$_{\odot}$\,yr$^{-1}$ (based on their FUV luminosities) and
their redshifts place them well within the end of cosmic reionization.
They also appear to be very numerous: Tanaguchi et al.\ 2005 find
$\sim$30 LAEs in only a quarter degree field (e.g., compared to
$\sim$10 QSOs at z$>$6 which are distributed over a quarter of the
sky! Fan et al.\ 2004). This implies that LAEs may play an important
role in reionizing the universe at z$>$6 (for a review see Fan,
Carilli \& Keating, 2006).  Investigating the physical properties of
these sources are thus of great interest and ALMA will play a critical
role in studying these objects, as discussed in the following.

\begin{figure}
\begin{center}
  \includegraphics[scale=0.45]{walter_fig1.eps}
\end{center}
\caption{ALMA CO `discocery space': The horizontal lines indicate
  which CO transition (plotted on the y--axis) can be observed with
  which ALMA band as a function of redshift (plotted on the x--axis).
  For objects with z$>$7, only the higher--J CO transitions can be
  observed with ALMA. The [CII] `discovery space' is also indicated.}
\label{fig:1}
\end{figure}

\section{Interstellar Medium: CO vs. [CII] emission}
\vspace*{-2.5mm}

\subsection{Carbon Monoxide (CO)}
\vspace*{-2.5mm}

Constraining the properties of the molecular gas in objects at the end
of cosmic reionization is clearly of key importance as such
observations 1) will measure the available 'fuel' for star formation,
2) will help to constrain the dynamical mass of the system and will
thus 3) allow to put these objects in an evolutionary context for
early galaxy formation.  Typically, at low and high redshifts, CO
emission is used as a tracer for the molecular gas phase (e.g. review
by Solomon \& Vanden Bout 2004). It is important to keep in mind
though, that, at the highest z, only the very high rotational lines of
CO will be observable with ALMA.  E.g. even in the lowest (currently
funded) frequency band of ALMA (band 3, 84--119\,GHz), only CO
transitions with J$>$7 (i.e., CO(7--6), CO(8--7), etc.) will be
observable at z$>$7. This is graphically illustrated in Fig.~1 where
we plot ALMA's `CO discovery space' (i.e., which line can be observed
at which redshift using which ALMA band). The high--J transitions
correspond to highly excited gas (either due to high kinetic
temperatures, high densities, or both) which may not be excited in
normal starforming environments.  This is shown in Fig.~2 (taken from
Weiss et al.\ 2005) where measured CO line strengths (as a function of
J, this is sometimes referred to as CO line ladders/SEDs) are plotted
for a number of key sources. What is immediately obvious from this
plot is that most objects have sharply decreasing CO line strengths
beyond J$>8$, in particular starforming systems such as NGC\,253, or
the sub--millimeter galaxy plotted in this diagram (the quasars appear
to be more excited, but their CO line SED still turns over at
J$\sim$7, for an exceptional object see APM\,07279, Weiss et al.\
2007). This comparison immediately implies that emission from the CO
molecule will typically be very difficult to observe with ALMA at
z$>$7 as the observable lines will simply not be excited.

\begin{figure}
\begin{center}
  \includegraphics[scale=0.45]{walter_fig2.ps}
\end{center}
\caption{Comparison of various normalized (by their CO(1--0) flux
  density) CO line SEDs at low and high redshift (figure taken from
  Weiss et al.\ 2005). The CO line SEDs decline rapidly beyond
  J=6--8.}
\label{fig:2}
\end{figure}

\subsection{Ionized Carbon ([CII]) to the Rescue!}
\vspace*{-2.5mm}

An alternative tracer of the interstellar medium is one of the main
cooling line of the ISM, the $^2$P$_{3/2}$ $\to$ $^2$P$_{1/2}$
fine--structure line of C$^+$ (or [CII]).  In brief, the [CII] line is
expected to be much stronger than any of the CO lines.  Given its high
frequency (157.74\,$\mu$m, corresponding to 1900.54\,GHz) [CII]
studies in the local universe are limited to airborne or satellite
missions (e.g.  Stacey et al.\ 1991, Malhotra et al.\ 1997, Madden et
al.\ 1997).  These studies have demonstrated that this single line can
indeed carry a good fraction of the total infrared luminosity (L$_{\rm
  FIR}$) of an entire galaxy.  In the local universe, the ratio
L$_{\rm CII}$/L$_{\rm FIR}$ has been found to be 2--5$\times$10$^{-4}$
in the case of ULIRGS (e.g. Gerin \& Phillips 2000), but is more like
5--10$\times$10$^{-3}$ in more typical starforming galaxies (for a
discussion on possible reasons for the supressed ratio in ULIRGs see,
e.g., Luhman et al.\ 1998). Notably, the ratio has been found to be
1\% or even higher in low metallicity environments. E.g., in the low-
metallicity galaxy IC10, L$_{\rm CII}$/L$_{\rm FIR}$ reaches values as
high as 4\%, with an average value of 2\% (Madden et al.\ 1997, see
Israel et al.\ 1996 for a similar result for the LMC). This is the
reason why it has long been argued (e.g., Stark 1997) that observation
of the [CII] line of prestine systems at the highest redshifts will
likely be the key to study molecular gas in the earliest starforming
systems, in particular in the era of ALMA. The ALMA [CII] `discovery
space' is also indicated in Fig.~1.

\begin{figure}
\begin{center}
  \includegraphics[scale=0.4]{walter_fig3.ps}
\end{center}
\caption{{\em Top:} First detection of [CII] at high redshift in the
  z=6.42 QSO J1148+5251 (Maiolino et al.\ 2005). {\em Bottom:}
  brightest CO transition (J=6) in the same source (Bertoldi et al.\
  2003, Walter et al.\ 2003). Note that the [CII] line is brighter by a
  factor of $\sim$5.}
\label{fig:3}
\end{figure}

\section{Expected [CII] line strengths}
\vspace*{-2.5mm}

At the redshifts of the LAEs, the [CII] line is shifted to the 1\,mm
band of ALMA (band 6, 211--275\,GHz). [CII] emission has recently been
successfully detected using the IRAM 30\,m in the highest redshift
quasar J1148+5251 at z=6.42 (Maiolino et al.\ 2005, see Fig.~3).  The
noteable difference between J1148+5251 and the z$>$6 LAEs is that the
ratio L$_{\rm CII}$/L$_{\rm FIR}$ has been found to be very low
($\sim$5$\times$10$^{-4}$) in J1148+5251, i.e. in perfect agreement
with studies of low redshift ULIRGs that show a central AGN.  On the
contrary, the LAE are presumably pure starbursts (no evidence for an
AGN is found, Taniguchi et al.\ 2005) and they likely have lower
metallicities compared to the highly overdense regions in which the
luminous quasars are supposedly present. All these arguments point
towards a L$_{\rm CII}$/L$_{\rm FIR}$ ratio in LAE that is close to
what is found for nearby normal galaxies, or perhaps even for the
metal--poor dwarf galaxies (i.e., around 1\% or even higher). In other
words, the [CII] luminosity of the LAE may well be an order of
magnitude stronger (for a given IR luminosity) than what has been
found in the z=6.4 QSO. In the following we present a quick
back--of--the envelope calculation based on the detected [CII] line
strength in J\,1148+5251 ($\sim$10~mJy) which has a SFR of a few 1000
M$_\odot$\,yr$^{-1}$. This SFR is more than two order of magnitudes
higher than the SFR found in a typical LAE, but as the L$_{\rm
CII}$/L$_{\rm FIR}$ may be higher by an order of magnitude in the
LAEs, the expected [CII] line strength of the LAE may be as high as
1\,mJy. Such a line should be easily detectable with ALMA at high
significance in a few hours.

\section{Resolving the ISM}
\vspace*{-2.5mm}

{\em Detecting} the [CII] (or CO) emission is critical to estimate the
reservoir of the (molecular) gas in these early systems. A second step
is then to  spatially resolve the molecular gas distribution. In
particular, given the typical diameters of galaxies of many kpc, a
linear resolution of $\sim 1$\,kpc is needed to resolve the structure
of the underlying galaxy. Such measurements are needed 1) to get an
estimate for the size of the galaxy (and thus a better estimate
for the dynamical mass), 2) to resolve potentially merging systems and
3) to better constrain the physical properties of the gas (e.g.,
by measuring the brightness temperature of the hosts).  A linear
resolution of 1\,kpc corresponds to a resolution of 0.15$''$ at the
redshifts under consideration (1$''\sim5.8$\,kpc at $z=6$).  Such
observations can then in turn be used to constrain the predictions by
CDM simulations of early galaxy formation, and, if a large sample was
available, put limits on the frequency of mergers at high redshift. In
addition, such studies can be used to constrain the possible
redshift-evolution of the M$_{\rm BH}$--$\sigma_{\rm v}$ relation in
high--z quasars. Such observations will clearly be feasible with ALMA
in the extended arrays. High--resolution CO imaging is already
possible with the current generation of telescopes: we have used the
VLA to resolve the molecular gas in the host galaxy of the z=6.42 QSO
J\,1148+5251 (see Fig.~4, Walter et al.\ 2004).

\section{The case for ALMA band 5}
\vspace*{-2.5mm}

As a technical note: the ALMA redshift coverage for the [CII] line is
not ideal as its frequency lies between the CO(17--16) and CO(16--15)
transition (see Fig.~1). One concern is that the critical redshift
range (8$<$z$<$10.5) is currently not fully covered: This frequency
range corresponds to the ALMA band 5 which is only partly funded by
the European Union as part of the Sixth Framework Programme (FP6) for
up to 8 antennas. Clearly, it would be highly desireable to equip as
many ALMA antennas with band 5 receivers as possible.

\begin{figure}
\begin{center}
  \includegraphics[scale=0.5]{walter_fig4.eps}
\end{center}
\caption{High--resolution CO image of the z=6.42 QSO J\,1148+5251
  obtained at the VLA (Walter et al.\ 2004). The resolution achieved
  in these observations (0.15$''$, corresponding to $\sim$1\,kpc) will
  be routinely reached with ALMA.}
\label{fig:4}
\end{figure}

\section{Concluding remarks}

ALMA observations of the [CII] line will play a fundamental role in
studying the youngest galaxies in the Epoch of Cosmic Reionization at
z$>$7. Given the expected line strengths it should be possible to
resolve these galaxies in the [CII] line emission on kpc scales. Such
measurements would not only constrain the sizes but would also help to
derive the dynamical masses in these early starforming systems. Given
the typical CO excitation in starforming galaxies (i.e. drop in
excitation around the J$\sim$6 transitions), ALMA will likely act as a
[CII]-- rather than a CO--machine for objects at these extreme
redshifts.


New clues on non standard mixing on the RGB

Abstract: We discuss the observational evidences for the non-standard mixing that occurs in low mass stars while they are ascending the Red Giant Branch.

\section{Chemical anomalies in low mass red giants}
It has been known for a long time that low mass red giants
(i.e., with initial mass lower than $\sim$ 2 M$_{\odot}$)
present chemical anomalies which are not predicted by the standard
evolution theory.
In particular, the models which neglect the transport of chemicals
in the stellar radiative regions can not explain the observed behavior
of the carbon isotopic ratio both in Pop I and II stars on the red giant
branch (RGB), of the lithium and carbon abundances (respectively in halo and
globular cluster giants) after the completion of the first dredge-up,
neither the evidences of O versus N anticorrelation and Na and Al versus N
correlation seen in a large number of globular cluster giants
(see Charbonnel, Brown \& Wallerstein 1998, hereafter CBW98, for references).
These observations suggest that low mass stars undergo a non-standard
mixing process which adds to the first dredge-up and modifies the surface
abundances.

\section{When does the extra-mixing occur?}

To throw a new light on the problem of the onset of the extra-mixing, and of
its possible metallicity dependance, we assembled accurate abundance
observations for the $\rm ^{12}C/^{13}C$ ratio as well as for Li, C, and N
for five field giants with [Fe/H] $\simeq - 0.6$ and two stars in the globular
cluster 47 Tuc.
Using the HIPPARCOS parallaxes, we could constrain the evolutionary
status of the sample stars.

Our data
can be viewed as an evolutionary sequence, as shown in Fig.1 where
the $\rm ^{12}C/^{13}C$ ratio is plotted against M$_{\rm bol}$ and where
the position of the luminosity function bump of 47 Tuc
(at M$_{\rm bol}=+1.05 \pm 0.2$) is indicated.
Along this sequence, the behavior of the $\rm ^{12}C/^{13}C$ ratio is
very similar to the one observed at solar metallicity in the open
cluster M67 (Gilroy \& Brown 1991) :
Below the bump, the observations are in agreement with standard predictions
for dilution ($\sim 20$). Then between $\rm M_{bol}$ = +1 and +0.5, the
$\rm ^{12}C/^{13}C$ ratios drop from $\sim 20$ to near 7, i.e., well below
the standard predicted ratio (at the same time, Li disappears and the
$\rm ^{12}C/^{14}N$ ratio diminishes by 0.2 to 0.4 dex
{\footnote{see Tables 1 and 2 in CBW98}}).
Finaly, there is no further change in the $\rm ^{12}C/^{13}C$ ratio from
$\rm M_{bol}$ = +0.4 to $-$2.

\begin{figure}[ht]
\plotfiddle{fig1STScI98.eps}{70truemm}{0}{40}{40}{-160}{-75}
\caption{The ratio of $\rm ^{12}C/^{13}C$ is plotted against the
bolometric magnitude of our stars with moderate metal deficiency.
The absolute magnitudes are derived from the HIPPARCOS parallaxes.
The horizontal line indicates the position of the bump in the
luminosity function of 47 Tuc (King et al. 1985).}
\end{figure}

These results show that the extra-mixing process which leads to very
low $\rm ^{12}C/^{13}C$ ratios in metal-deficient low-mass evolved stars
becomes efficient exactly after the luminosity function bump, i.e., the
evolutionary point where the hydrogen-burning shell crosses the chemical
discontinuity created by the outward moving convective envelope.
This confirms the inhibiting effect of molecular weight (or $\mu$)
barriers against the development of the extra-mixing (see Sweigart \&
Mengel 1979; Charbonnel 1994, 1995).
Since stars more massive than $\sim$ 1.7 - 2.2 M$_{\odot}$ do not reach
the bump, they are not expected to experience this extra-mixing on the RGB
(this is confirmed by the observations in open clusters; Gilroy 1989).

\section{How many low mass stars do undergo the extra-mixing on the RGB?}
To answer this question, we collected in the literature all the
observations of the carbon isotopic ratio in field and cluster giant stars
(see Charbonnel \& Dias 1998).
For the whole sample stars, M$_{\rm V}$ values were determined thanks to
the HIPPARCOS parallaxes.
In the [Fe/H] range we consider, M$_{\rm V}^{\rm bump}$ is higher than
$\sim$ -0.2.

\begin{figure}[ht]
\plotfiddle{fig2STScI98.eps}{70truemm}{0}{40}{40}{-160}{-75}
\caption{Carbon isotopic ratio collected in the literature as a function
of [Fe/H], for field and cluster giants (black and white symbols
respectively). The circles correspond to stars with M$_{\rm V}$ higher
than 2, the triangles to with M$_{\rm V}$ between 0 and 2, and the
squares to the brightest stars with M$_{\rm V}$ lower than 0}
\end{figure}

Above the bump (squares in Fig.2), i.e., the evolutionary point
where the extra-mixing can occur, the disagreement between the post
dredge-up standard predictions ($^{12}$C/$^{13}$C around 20-30) and the
observations appears in 96$\%$ of the low mass giants.
This high value for the fraction of low mass stars that experience an
extra-mixing process on the RGB and are thus expected to destroy their
$^3$He at this evolutionary phase (Hogan 1995, Charbonnel 1995)
satisfies the galactic requirements for the evolution of the $^3$He
abundance (see Galli et al. 1997 and Charbonnel \& Dias 1998).

\section{Critical $\mu$-gradient}
As can be seen in Fig.2, the final $\rm ^{12}C/^{13}C$ ratios are lower
in the more metal-poor stars. The observations at various metallicities
provide precise clues on the extension of the extra-mixed region,
and on the $\mu$-barriers that may shield the central regions of a star
from extra-mixing.

From evolutionary models of stars with various masses and metallicities,
we showed that, in order to account for the $^{12}$C/$^{13}$C ratios
observed in stars of different metallicities, the extra-mixing must extend
from the base of the external convection zone down to a region where the
gradient of molecular weight is equal to $ \sim 1.5 \times 10^{-13}$.
This value for the ``observational" critical $\mu$-gradient, ($\nabla
\ln \mu)_{c,obs}$, which appears to shield the central regions of the star
from extra-mixing agrees with the one which is expected to stabilize
meridional circulation.
This result indicates that the extra-mixing on the RGB could be related to
rotation-induced processes.
If we assume that such a $\mu$-barrier can not be penetrated, then
the energy production in the hydrogen-burning shell should not be
affected.

Let us note that, on the main sequence, $\nabla \ln \mu$ becomes higher
than ($\nabla \ln \mu)_{c,obs}$ in the very external part of the low
mass stars, and no extra-mixing is expected to occur in the stellar region
of energy production.
This explains the perfect agreement between standard theoretical
dilution and observations of $\rm ^{12}C/^{13}C$ for stars that have not
yet reached the bump.
This is also in agreement with the best solar models
(for what concerns helioseismological comparison and agreement with Li and
Be observations; Richard et al. 1996)
that include both element segregation and rotation-induced mixing which
is cut-off when the $\mu$-gradient becomes $\geq$ to 1.5 - 4 $\times 10^{-13}$.

\section{Conclusions}

The physics of the extra-mixing process in low mass stars has to be better
understood. Detailed simulations, with a consistent treatment of the
transport of matter and angular momentum, have to be carried out for
different stellar masses and metallicities, and various mass loss and
rotation histories.
The impact of this process on the behavior of various elements in RGB stars
(C $\searrow$, Na $\nearrow$, O $\searrow$, Al $\nearrow$, Mg O $\searrow$)
and on the precise yields of $^3$He has to be investigated
in details. Consequences for the energy production in the HBS and for
the HB morphology may not be neglected.


NMR study of the Mott transitions to superconductivity in the two Cs_3C_60 phases

Abstract: We report an NMR and magnetometry study on the expanded intercalated fulleride Cs_3C_60 in both its A15 and face centered cubic structures. NMR allowed us to evidence that both exhibit a first-order Mott transition to a superconducting (SC) state, occuring at distinct critical pressures p_c and temperatures T_c. Though the ground state magnetism of the Mott phases differs, their high $T$ paramagnetic and SC properties are found similar, and the phase diagrams versus unit volume per C_60 are superimposed. Thus, as expected for a strongly correlated system, the inter-ball distance is the relvevant parameter driving the electronic behavior and quantum transitions of these systems.

\section{}

High transition temperature ($T_{c}$) superconductivity (SC) in the vicinity
of magnetic phases is rather commonly found nowadays, especially in
materials involving transition metal ions. There, the incidence of the spin
fluctuations on the SC state and the symmetry of its order parameter is
being highly debated and various possibilities have been considered
depending whether the electronic states involved at the Fermi level reside
on a single orbital (as in the cuprates) or multiorbital occupancy occurs
(as for Fe pnictides).

Alkali doped fullerides $A_{3}$C$_{60}$ ($A=$ alkali ion) represent a
distinct family of HTSC in a multiorbital case, as the lowest unoccupied C
_{60}$ molecular orbitals (named $t_{1u}$) are sixfold degenerate \cit
{gunnarsson-RMP69}. The importance of electron correlations in $A_{n}$C$_{60}
$ compounds has been suggested first from the detection of insulating states
for even $n$, contrary to the expected metallicity for any $1\leq n\leq 5$
in an independent electron picture. On-ball localization of pairs of
electrons is favored by the energy gain due to Jahn-Teller (JT) distorsions
of the charged C$_{60}$ molecules which adds to the coulomb energy $U$. The
occurrence of this Mott Jahn-Teller insulator \cite{capone-PRB62}, in which
Hund's rule is violated, has been confirmed by the observation by NMR of a
spin gap due to singlet-triplet excitations, smaller than the optical gap,
for both crystal structures of $A_{4}$C$_{60}$ \cite{kerkoud-JPCS57} and Na
_{2}$C$_{60}$ \cite{brouet-PRB66}. Furthermore for $n=1$, in fcc-CsC$_{60}$,
the observation of charge segregation of singlet electron pairs on a sizable
fraction of the C$_{60}$ balls at low $T$ has been an even more direct
evidence that electron pairs are favored by JT\ distortions \cit
{brouet-PRL82}.

For $n=3,$ JT distortions were expected to be less effective. So the $s
-wave SC evidenced in the fcc-$A_{3}$C$_{60}$, and the scaling of $T_{c}$
with the distance between C$_{60}$ balls has been interpreted by most
researchers as purely BCS, driven by on ball phonons, with weak incidence of
electronic correlations \cite{gunnarsson-RMP69}. The effort to increase
T_{c}$ by expanding the C$_{60}$ lattice has, however, led to the discovery
of various magnetic compounds, such as $($NH$_{3})_{x}$K$_{3}$C$_{60}$,
which displays a Mott transition to a SC state under pressure \cit
{prassides-JACS121}. But one was still led to suspect that this behavior
could be attributed to a lifting of the degeneracy of the $t_{1u}$ levels by
the peculiar lattice structure required to expand the C$_{60}$ lattice \cit
{rosseinsky-Nature364}. The most expanded fulleride Cs$_{3}$C$_{60}$ had
been found to become SC under pressure ($p$), with $T_{c}\simeq 40$ K \cit
{palstra-SSC93}, but it has only recently been shown that this occurs in a
cubic A15 structure \cite{ganin-NM7}. A renewed interest arises then as this
phase is antiferromagnetic (AF) at ambient $p$, with $T_{N}=47$ K \cit
{takabayashi-Science323}, and undergoes a Mott transition to a metallic
state for $p\sim 4$ kbar. It is then of great interest to find out whether
the electronic properties of A15-Cs$_{3}$C$_{60}$ exhibit any difference
with those of other fcc-$A_{3}$C$_{60}$ phases.

Ganin \textit{et al.} indicated that the low-temperature reaction process
\cite{ganin-NM7} used to synthesize A15-Cs$_{3}$C$_{60}$ produces mixed
phases including fcc-Cs$_{3}$C$_{60}$ and body centered orthorhombic bco-Cs
_{4}$C$_{60}$. Using the spectroscopic capabilities of $^{133}$Cs NMR
experiments, we sorted out the signals from the two Cs$_{3}$C$_{60}$ isomers
in such mixed-phase samples.\ Taking that in advantage, we report in this
letter the first direct comparison, and demonstrate that a Mott transition
to SC occurs as well in fcc-Cs$_{3}$C$_{60}$. At ambient pressure, we
evidence the decrease of spin freezing temperature for the fcc as compared
to that of the A15 phase and associate it with the geometrical frustration
of the former lattice. The occurrence of a Mott transiton is shown to be
independent of the crystal structure and of the specific C$_{60}$ ball
orientational (merohedral) disorder present in fcc-Cs$_{3}$C$_{60}$. This
confers then a very important place to these phases in helping to reach an
understanding of SC in the vicinity of magnetic phases.

\begin{figure}[tbp]
\begin{center}
\includegraphics[width=7.5cm]{FIG1.eps}
\end{center}
\caption{a) The A15 $^{133}$Cs NMR spectrum taken at ambient pressure
displays a seven-peak powder NMR spectrum (arrows), with a quadrupole
splitting of $\sim 39$ kHz at $100$ K, and a slight asymmetric spectral
shape originating from a Knight shift anisotropy. b) $^{133}$Cs NMR spectra
for the fcc phase displaying the octahedral O and tetrahedral sites peaks, T
and T'. c) and d) The spectra are found nearly unmodified at $7.6$ kbar when
both phases are metallic (see text). }
\label{fig1}
\end{figure}

$^{133}$Cs \textit{NMR spectra of the two phases.}--- We could synthesize
mixed-phase samples and selected three of them with significantly differing
phase contents, labelled A1 and A2 for A15 rich and F1 for the fcc rich.
Their compositions are A1\ (58.4, 12, 29.5), A2\ (41.7, 12, 46.5) , F1( 34,
55, 11), where the \% contents in formula units are given respectively for
the A15, fcc and bco phases.
As $^{133}$Cs has a nuclear spin $I=7/2$, its NMR spectrum is
sensitive to the local site symmetry through the coupling of the nuclear
quadrupole moment with the electric field gradient (EFG) induced by the
local charge distribution. So in the A15 phase the single NMR site displays
a quadrupole-split seven-line spectrum (Fig.~\ref{fig1}a), as the unique Cs
site displays a non-cubic local symmetry \cite{jeglic-PRB80}. In the fcc
phase, the unit cell contains two alkali sites, with occupancy ratio 1:2 for
the octahedral (O) and tetrahedral (T) sites. Their local symmetry being
cubic, the EFG vanishes and each site has a narrow non-split signal. A
peculiarity evidenced by NMR \cite{walstedt-Nature362} in all formerly known
fcc-$A_{3}$C$_{60}$ is that the tetrahedral site splits into two sites (T
and T') which have been assigned to the merohedral disorder of C$_{60}$
balls \cite{matus-PRB74}. The detection in sample F1 of these three lines
(Fig.~\ref{fig1}b) establishes then the identical structure of fcc-Cs$_{3}$C
_{60}$. This difference in NMR spectra allowed us, as done indeed in Fig.
\ref{fig1}, to detect selectively the $^{133}$Cs NMR of a given phase.
Let us point out, as will be discussed
later, that the data of Fig.~\ref{fig1}c demonstrates that these structures
are not modified under pressure.
\begin{figure}[tbp]
\begin{center}
\includegraphics[width=7cm]{FIG2.eps}
\end{center}
\caption{(color on line) Ambient pressure data taken in samples A1 and F1
for the A15 and fcc phases. a) Similar variations above 100K\ of the
anisotropic shift contribution $K^{\mathrm{ax}}$ to the $^{13}$C spectra
(definition recalled in the inset), The scale chosen on the right for
\protect\chi $ (black dots, SQUID\ data on A1) emphasizes the linear
relation with $K^{\mathrm{ax}}$. b) $^{133}$Cs $(T_{1}T)^{-1}$ data for the
fcc phase T site and the A15 phase have similar high $T$ variations but
differ markedly near and below their ordered magnetic states. The spin
freezing is monitored in the inset by the variation of the $^{13}$C and
^{133}$Cs signal intensities taken on sample F1(see text). }
\label{fig2}
\end{figure}

\textit{Paramagnetism at ambient pressure.}--- SQUID data on all samples do
not display any superconductivity for $p=1$ bar, and only exhibit a
paramagnetic susceptibility
\begin{equation*}
\chi (T)=\chi _{orb}+\chi _{s}(T),
\end{equation*
which includes an orbital term and a $T$ dependence due to the spin
magnetism of unpaired electrons. The data were similar to those attributed
to the A15 phase \cite{takabayashi-Science323}. Comparisons between the two
phases are possible from analyses of the $^{13}$C NMR spectra, as the NMR
shift involves also orbital and spin components
\begin{equation*}
K^{\alpha }=K_{orb}^{\alpha }+K_{s}^{\alpha }(T)=K_{orb}^{\alpha }+\
A^{\alpha }\chi _{s}(T).
\end{equation*}

Here, index $\alpha $ refers to the direction of the appplied field $B$ with
respect to local axes on the $^{13}$C site.\ Indeed, both $K_{orb}^{\alpha }$
due to the orbital magnetism of the $sp^{2}$ bonding electrons \cit
{tycko-PRL67}, and $K_{s}^{\alpha }(T)$ associated with the spin
magnetization of electrons in the $t_{1u}\ $orbitals are anisotropic \cit
{pennington-RMP68}. The random orientation of the balls with respect to $B$
gives a typical $^{13}$C powder NMR\ spectrum, with two singularities for
K^{\perp }$ and $K^{||}\ $which correspond to $B$\ directions $\perp $ and
||$ to the tangential plane to the C$_{60}$\ ball at the $^{13}$C\ site, as
shown in the inset of Fig.~\ref{fig2}a.\ One can notice there that, at
T=160 $ K, the spectra are identical for the A15\ and fcc rich samples.

Furthermore, in Fig.~\ref{fig2}a the $T$ variations of the anisotropy $K^
\mathrm{ax}}=2(K^{\perp }-K^{||})/3$ , which is obtained from fits of the
spectra, \textit{cannot be differenciated for the two phases above $100$ K
and track those of SQUID data for $\chi $ taken on sample A1. The deduced
data for $\chi _{s}(T)$ are then intrinsic and similar for both phases and
can be fitted above $100$ K with a Curie-Weiss law with an effective moment
p_{\mathrm{eff}}=1.52(5)\ \mu _{B}$ per C$_{60}$ and a Weiss temperature
\Theta _{W}=-70\pm 5$ K. Such a value for $p_{\mathrm{eff}}$ let us suggest
that the C$_{60}^{3-}$ ion is \textit{in a low spin state} in the fcc phase
as established before for the A15 phase \cite{takabayashi-Science323,
jeglic-PRB80}.

\textit{Magnetic ordering at ambient pressure.}--- For all samples the SQUID
data exhibits a sharp increase of  magnetization at $T=47$ K due to the AF
state of the fraction of A15\  phase as shown in Fig.~\ref{fig2}a. We also
detected the increased linewidth of the A15 phase $^{133}$Cs NMR signal
found in Ref.~\cite{takabayashi-Science323, jeglic-PRB80}. We could as well
find in the A15 rich sample that the magnetization on the C$_{60}$ ball
induces such a large broadening of the $^{13}$C NMR that its intensity,
detected within a small frequency window of $100$ kHz (that is $0.01$ T),
vanishes below $T_{N}=47$ K. This allowed us then in sample F1, after the
loss of the A15 fraction ($\sim 40$ \%),\  to isolate below 47K the  fcc
^{13}$C NMR.\ The paramagnetism of this phase,  given by  $K^{\mathrm{ax}}$,
decreases slightly with respect to its high $T$ value, as seen in Fig.~\re
{fig2}a.

\begin{figure}[tbp]
\begin{center}
\includegraphics[width=6.5cm]{FIG3.eps}
\end{center}
\caption{(color on line) SC diamagnetism is detected at a lower $p_{c}$ for
the sample F1 than for A2 (nearly pure A15). The arrows point the values of
T_{c}(p)$ deduced by extrapolation of the sharp $T$ dependent diamagnetism
(slightly smaller than the onset values). The diamagnetic magnetization at
10K (extrapolated for some $p$ values) is displayed in the inset for both
samples. The magnetic volume fraction is also
shown for the A15 phase. }
\label{fig3}
\end{figure}

Spin freezing in the fcc is only detected below $T_{f}\simeq 10$ K from the
sharp intensity drops of both the$^{13}$C and  the fcc- specific $^{133}$Cs
NMR signal, shown in the insert of Fig.~\ref{fig2}b. This latter observation
establishes that in the fcc-frozen spin state the internal field on the
^{133}$Cs site is much larger than that in the A15\ AF\ phase for which the
integrated intensity over a $100$ kHz window only slightly declines below
T_{N}$. This confirms that the transferred internal field on $^{133}$Cs is
partly compensated in the A15 AF phase due to the bipartite\ body centered
lattice symmetry \cite{jeglic-PRB80}.\ On the contrary this large internal
fields on $^{133}$Cs and the small $T_{f}\ $\ value give evidence that the
fcc phase magnetic state is influenced by the inherent frustration of the
fcc lattice.

To compare the dynamical magnetic properties, $^{133}$Cs spin lattice
relaxation $T_{1}$ data have been taken on the two phases. In the fcc, the
less shifted O site could only be resolved above $50$ K, and its $T_{1}$
scales by a factor three with that of the T site, as seen in Fig.~\ref{fig2
b. So both sites sense the same magnetic fluctuations through distinct
hyperfine couplings. For the T site, $(T_{1}T)^{-1}$ follows in the
paramagnetic regime a similar variation as that seen in the A15 phase.
There, below $T_{N}=47$ K \ a sharp gap in the spin excitations is detected
and the data can be fitted with $(T_{1}T)^{-1}\propto \exp (-\Delta /k_{B}T)
, with $\Delta \simeq 50$ K (full line on Fig.~\ref{fig2}b). On the contrary
in the fcc phase $(T_{1}T)^{-1}\ $ remains nearly constant down to $10$ K
and even displays a fast increase for the $^{133}$Cs nuclear spins which are
not submitted to a large static field. This persistence of spin fluctuations
at low $T$ is also to be linked to frustration effects.
\begin{figure}[tbp]
\includegraphics[width=7.5cm]{FIG4.eps}
\caption{(color on line) Phase diagram reporting the critical temperatures
(a) versus $p$, (b) versus $V_{C_{60}}$. The pressures $p_{c}$ of the Mott
transitions shown as dotted lines in (a) merge into the gray area in (b).
There we reported for comparison the former data on A15 phase \protect\cit
{takabayashi-Science323} as a dotted line, and on the other fcc-$A_{3}$C
_{60}$ as empty downward triangles \protect\cite{maniwa-JPSJ63}. }
\label{fig4}
\end{figure}

\textit{SC phases under pressure and phase diagram.}--- SC diamagnetism
could be probed \textit{in situ} by monitoring the shift of the NMR coil
tuning frequency ($f\simeq 28$ MHz). This allowed us to estimate the SC
volume fractions, as we kept the geometry unchanged during the pressure
sweeps. We find here that the fcc phase displays a transition to a SC state
at $3.9$ kbar (sample F1), while this occurs only above $5.5$ kbar in the
A15 rich (A2), which points out that the critical pressure $p_{c}$ for the
magnetic-SC transition is lower for the fcc phase than for the A15. This can
be seen in the insert of Fig.~\ref{fig3}, where the magnitude of the
diamagnetic signal measured at $10$ K is plotted versus $p$. It is as well
compared for the A15 phase with the magnetic volume fraction estimated from
the $^{133}$Cs NMR. The abrupt loss of magnetism above
7 $ kbar and onset of SC above $6$ kbar points for the \textit{first order
AF/SC phase boundaries} at $p_{c}=6.5(5)$ kbar in the A15 and $p_{c}=3.5(5)$
kbar in the fcc phase.

We can then report the $(p,T)$ phase diagrams in Fig.~\ref{fig4}a, and do
find in Fig.~\ref{fig4}b that $p_{c}$ and $T_{c}$ merge together for the two
Cs$_{3}$C$_{60}$ phases if plotted versus $V_{C_{60}}$, the unit volume per
$_{60}$ ball \cite{footnote}. There, we can as well compare the present
T_{c} $ data with former results on the other fcc-$A_{3}$C$_{60}$, and we
evidence that a similar maximum of $T_{c}$ versus $V_{C_{60}}$ applies for
the two structures.

Let us now consider the NMR data taken at high enough $p,$ for which both
phases become SC at low $T$.\ It is first clear, as seen in Fig.~\ref{fig1},
that the spectra above $T_{c}$ do not differ from those taken at one bar.
This absence of structural modification for both phases establishes then
that the evolution with $p$ only implies electronic degrees of freedom.
\begin{figure}[tbp]
\begin{center}
\includegraphics[width=7.5cm]{FIG5.eps}
\end{center}
\caption{$T$ variations of $(T_{1}T)^{-1}$ for $^{133}$Cs and $^{13}$C (left
scale) and of the $^{133}$Cs shift $K$ (blue empty diamonds right scale) in
the A15 phase above $p_{c}$. The magnitude of the T site $^{133}$Cs
relaxation rate in the fcc compound is found similar to that obtained in the
A15 compound (see text).}
\label{fig5}
\end{figure}

In the metallic state, a Korringa-like $T$ independent $T_{1}T$ is seen
above $T_{c}$ for the two phases as shown in Fig.~\ref{fig5}. The constant
values are similar for $^{133}$Cs in the A15\ as for the T site of the fcc,
which displayed similar $T_{1}$ in the paramagnetic phases as well (Fig.~\re
{fig2}b). This points that the spin dynamics has comparable evolution with
pressure in the two phases. In the A15, the opening of the SC gap at $T_{c}
, as detected from the onset of decrease in the $^{133}$Cs Knight shift
^{133}K,$ occurs slightly above the observed decrease in $(T_{1}T)^{-1}$for
both $^{13}$C and $^{133}$Cs nuclei. This result perfectly mimics the
observations done in high applied fields in the other fcc-$A_{3}$C$_{60}$
compounds \cite{stenger-PRL74}. There, such a persistence of spin
excitations slightly below $T_{c}$ are remnants of the $s$-wave BCS like
Hebel-Slichter coherence peak \cite{coherencepeak}, which is damped in high
fields and could only be fully revealed from low field data \cit
{gunnarsson-RMP69}.

\textit{Summary and discussion.}--- In conclusion we have shown here that
the magnetic ground states of the Cs$_{3}$C$_{60}$\ phases are quite
distinct at ambient pressure. The reduction of ordering temperature and the
persistence of low-energy spin fluctuations at low $T$ in the fcc phase can
be assigned to the frustration effects inherent to this structure and its
merohedral disorder. We further evidenced that no major crystal structure
modification occurs under pressure as shown as well from x-ray spectra in
the A15 phase \cite{takabayashi-Science323}. So the low $T$ transition from
a magnetic to a SC state, which appears to be of first order, is fully
determined by electronic parameters in both cases. Comparison of the phase
diagrams demonstrates that the crititical pressure $p_{c}$ for the
transition occurs for a similar value of the volume per C$_{60}$ ball
V_{C_{60}}$, which highlights the Mottness of the transition to be opposed
to a CDW/SDW transition, which should sensitively depend on lattice and
Fermi surface symmetry.

Our data evidences that fcc-$A_{3}$C$_{60}$ compounds exhibit a dome
behavior of $T_{c}(V_{C_{60}})$, identical to that found in the A15\ phase
under pressure.\ This unexpected feature in a purely BCS $s$-wave scenario
dominated by a density of state variation is then quite generic.\ Our result
gives then weight to the theoretical attempts to take correlation and JT\
effects into account in these systems \cite{capone-Science296}, which did
suggest such a behavior beforehand.\ We are presently investigating the
evolution of the spin dynamics accross the Mott transition, which should
permit more thorough comparisons with such theoretical approaches, and hope
that the present work will trigger diverse other experimental studies of the
electronic properties of the Cs$_{3}$C$_{60}$ phases across the Mott
transition.

We would like to thank V.~Brouet, M.~Fabrizio and E.~Tosatti for their
interest and for stimulating exchanges. Y.~I ackowledges financial support
from JSPS for his post-doctoral stay in Orsay.

Superconductivity (cond-mat.supr-con); Strongly Correlated Electrons (cond-mat.str-el)