Heavy Flavor Probes of Quark-Gluon Plasma

Abstract: Some aspects of heavy flavor probes of quark-gluon plasma (QGP) including quarkonia and open heavy flavor D- and B-mesons, aimed to study the properties of QGP, are discussed in this lecture note.

\section{Introduction}

 Lattice QCD (LQCD) calculations predict that at a critical temperature
  $T_c\simeq 170\,\mbox{MeV}$, corresponding to an
energy density  $\varepsilon_c \simeq 1\,\mathrm{GeV/fm^3}$,
nuclear matter undergoes a phase transition to a deconfined state
of quarks and gluons, called Quark-Gluon Plasma (QGP). At the
modern collider facilities such as CERN supersynchrotron (SPS)
(the nucleon-nucleon (NN) centre-of-mass energy for collisions of
the heaviest ions  is $\sqrt{s} =17.3\,\mbox{GeV}$), Relativistic
Heavy Ion Collider (RHIC) at Brookhaven ($\sqrt{s}
=200\,\mbox{GeV}$), and Large Hadron Collider (LHC) at CERN
($\sqrt{s} =5.5\,\mbox{TeV}$), whose heavy-ion program  will start
soon, heavy-ion collisions are used to attain the energy density,
exceeding $\varepsilon_c$. This makes the QCD phase transition
potentially realizable within the reach of the laboratory
experiments. The objective is then to identify and to assess
suitable QGP signatures, allowing to study the properties of QGP.
To that end, a variety of observables (probes) can be
used~\cite{Isayev:JPG1,Isayev:JPG2}. Further we will be mainly
interested in heavy-flavour probes of QGP, i.e., utilizing the
particles having c- and b-quarks.

A special role of heavy $Q = (c, b)$ quarks as probes of the
medium created in heavy-ion collision (HIC) resides on the fact
that their masses ($m_c\approx 1.3\,\mbox{GeV},\, m_b\approx
4.2\,\mbox{GeV}$) are significantly larger than the typically
attained ambient temperatures or other nonperturbative scales,
$m_Q\gg T_c, \Lambda_{QCD}=0.2\,\mbox{GeV}$~\cite{Isayev:RH}. This
has several implications: (i) The production of heavy quarks is
essentially constrained to the early, primordial stages of  HIC.
Hence, heavy quarks can probe the properties of the dense matter
produced early in the collision. (ii)~Thermalization of heavy
quarks is "delayed" relative to light quarks. One could expect
that heavy quarks could "thermalize" to a certain extent, but not
fully on a timescale of the lifetime of the QGP. Therefore, their
spectra could be significantly modified, but still would retain
memory about their interaction history, and, hence, represent a
valuable probe. (iii) RHIC, and especially LHC experiments allow
to reach very low parton momentum fractions $x$, where gluon
saturation effects become important. Heavy quarks are useful tools
to study gluon saturation, since, due to their large masses, charm
and bottom cross sections are calculable via perturbative QCD  and
their yield is sensitive to the initial gluon density.

The heavy-flavor hadrons we will be interested in include: 1) open
charm $D=(c \bar q)$ and open beauty $B=(b \bar q)$ mesons
composed of a heavy quark $Q=(c, b)$ and a light antiquark $\bar
q=(\bar u, \bar d)$. These mesons could be sensitive to the energy
density of the medium through the mechanism of in-medium energy
loss; 2)  hidden charm [charmonia=($c \bar c$)] and hidden beauty
[bottomonia=($b \bar b$)] mesons (called collectively heavy
quarkonia) being the bound states of the charm quark-antiquark, or
bottom quark-antiquark pairs, respectively. Heavy quarkonia could
be sensitive to the initial temperature of the system through the
dissociation due to color screening of the color charge that will
be discussed later.

  For detecting heavy flavor hadrons, different decay channels
are used. At LHC, the ALICE experiment for detecting open charm
will use  the reconstructed hadronic decays like $D^0\rightarrow
K^-\pi^+$, $D^+\rightarrow K^-\pi^+\pi^+$, etc., for detecting
open beauty ALICE will use inclusive single lepton decays
$B\rightarrow  e+X$ (midrapidity), $B\rightarrow \mu+X $ (forward
rapidity) and inclusive displaced charmonia decays $B\rightarrow
J/\psi (\rightarrow e^+e^-)+X$. At RHIC, in PHENIX and STAR
experiments  the measurement of the spectra of open heavy flavors
is based on the measurement of the spectra of heavy flavor (HF)
electrons and positrons [($e^++e^-)/2$] from the semileptonic
decays like $D^0\rightarrow K^-e^+\nu_e$, $D^+\rightarrow \bar K^0
e^+ \nu_e$, etc. These measurements are based on the fact that the
decay kinematics of HF electrons/positrons largely conserves the
spectral properties of the parent particles. Besides, the STAR
experiment has the capability to directly reconstruct open heavy
flavor mesons through the hadronic decay channels ($D^0\rightarrow
K\pi$, $D^\pm\rightarrow K\pi\pi$, etc.). At low $p_T$, the STAR
experiment also uses heavy flavor decay muons to provide open
heavy flavor measurements. In all of the above experiments,
quarkonia are detected through their dilepton decays $Q\bar
Q\rightarrow e^+e^-$ (midrapidity), $Q\bar Q\rightarrow
\mu^+\mu^-$ (forward rapidity).

\section{Heavy Flavor Probes of QGP: Quarkonia}

We begin the discussion of heavy flavor probes of QGP with heavy
quarkonia. The question we would like to address is: What could
happen with quarkonium yields in HIC  if QGP is really formed? In
the above mensioned experiments, quarkonia are detected in the
dielectron channel at midrapidity and in the dimuon channel at
forward rapidity. For example, in the dimuon channel, signals of
about $7\cdot 10^5$ $J/\psi$s and $10^4$ $\Upsilon$s are expected
in the most central Pb-Pb collisions at LHC/ALICE for one year of
data taking at nominal luminosity (see Fig.~\ref{Isayev:fig:1}
from~\cite{Isayev:JPG2}). If quarkonia will be completely
suppressed by QGP, these signals will disappear, otherwise, in the
case of enhancement of quarkonium yields, the peaks will be more
pronounced as compared to the signals obtained in $pp$ collisions
after scaling them with the corresponding number of independent
binary NN  collisions.

\begin{figure}[tb]\centering
\includegraphics[height=5.5cm,keepaspectratio]{fig1.eps}
\caption{Left: Dimuon invariant mass spectra (after subtraction of
the uncorrelated background) in central Pb-Pb collisions at
LHC/ALICE for one year of data taking ($10^6$ s) and luminosity
$5\cdot 10^{26}\,\mathrm{cm}^{-2}\mathrm{s}^{-1}$ in the  $J/\psi$
mass region (PYTHIA calculations). Right:~Same as  on the left,
but in the $\Upsilon$ mass region.
 } \label{Isayev:fig:1}
\end{figure}

Let us note that if the deconfinement phase transition really
takes place in a dense medium created in  HIC, then a color charge
in the QGP will be screened analogously to the Debye screening of
an electric charge in the electromagnetic plasma. As a result of
color Debye screening of the heavy quark interaction in QGP, the
binding energy of a bound state decreases and one can expect that
this would lead to the suppression of quarkonium yields in HIC.
This idea was first suggested by H. Matsui and
H.~Satz~\cite{Isayev:MS} who predicted that color Debye screening
will result in the suppression of $J/\psi$ meson ($c\bar c$ in the
$^3S_1$ state, $M=3.097\, \mathrm{GeV}$) yields. After that, the
$J/\psi$ suppression  was considered as one of the key probes for
the QGP formation in heavy ion collisions. $J/\psi$ is especially
promising because of the large production cross-section and
dilepton decay channels which make it easily detectable.

However, soon it was realized that besides melting $J/\psi$ mesons
in QGP due to the screening of the color charge, there are also a
few competing mechanisms which could explain the suppression of
$J/\psi$ production in heavy-ion collisions. These mechanisms are
referred to as cold nuclear matter (CNM) effects. The first CNM
 effect is the absorption of $J/\psi$ by
nuclear fragments from colliding nuclei. Let us consider, e.g.,
the proton-nucleus collision. Once produced in the hard primary
parton processes, $J/\psi$ has to cross the length $L$ of nuclear
matter, before exiting the nucleus, and, when traversing nuclear
matter, it can be absorbed by forthcoming nucleons of a nucleus.
The production cross section of $J/\psi$ is $p$-$A$ collision can
be parameterized as \begin{equation}
\sigma_{pA}^{J/\psi}=A\sigma_{pp}^{J/\psi}e^{-\sigma_{abs}^{J/\psi}\varrho
L},
\end{equation}
 where $\sigma_{pp}^{J/\psi}$ is the production cross section of
$J/\psi$ in $pp$ collisions and $\sigma_{abs}^{J/\psi}$ is the
nuclear absorption cross section. From  the global fit to the data
on charmonium  production in $p$-$A$ collisions the value of
$\sigma_{abs}^{J/\psi}$ can be extracted, in particular, at SPS
(NA50 experiment) it was obtained that $\sigma_{abs}^{J/\psi}=
4.2\pm 0.5$ mb~\cite{Isayev:Al}.

The second CNM effect is related to shadowing of low momentum
partons. This means the depletion of low momentum partons in
nucleons bound in nuclei as compared to free nucleons. This effect
can be accounted for in terms of the modification of the parton
distribution functions in nucleon within the nucleus with respect
to the parton distribution functions in a free nucleon:
\begin{equation} R_i^A(x,Q^2)=\frac{f_i^A(x,Q^2)}{
f_i^N(x,Q^2)}<1,\quad i=q_v,q_{sea},g.\end{equation}  Here "i"
denotes valence quarks, sea quarks, and gluons, $x$ is the parton
momentum fraction, $Q^2$ is the momentum transfer squared. At high
energies, $J/\psi$s are dominantly produced through the gluon
fusion,  and the $J/\psi$ yield is therefore sensitive to gluon
shadowing. The underlying idea explaining the occurrence  of gluon
shadowing  is that the gluon density strongly rises at small $x$
to the point where gluon fusion, $gg\rightarrow g$, becomes
significant. In the case of proton-nucleus and nucleus-nucleus
collisions, where nuclei with large mass number $A$ are involved,
the nonlinear effects are enhanced by the larger density of gluons
per unit transverse area of the colliding nuclei. A direct
consequence of nuclear shadowing is the reduction of
hard-scattering cross sections in the phase-space region
characterized by small-$x$ incoming partons. For gluons, e.g.,
shadowing becomes important at $x\lesssim5\times10^{-2}$, and,
hence, is relevant for the conditions of RHIC and LHC. Note
however that the strength of the reduction is constrained by the
current experimental data only for $x\gtrsim10^{-3}$.

After discussing CNM effects on the $J/\psi$ production, let us
consider the $J/\psi$ production at CERN SPS experiments.
Fig.~\ref{Isayev:fig:3} shows the $J/\psi$ production cross
section, where $J/\psi$ are registered through their decay in the
dimuon channel, normalized to the cross section of the Drell-Yan
(DY) process~\cite{Isayev:S}. The nuclear thickness parameter $L$
parameterizes the number of participating nucleons and increases
with centrality.  At SPS energies,  DY process is not modified by
the medium and its production cross section scales linearly with
the number of binary NN collisions. Hence, this process can be
used as unsuppressed reference process for the $J/\psi$
production. The current interpretation of the obtained results is
the following. The effects of cold nuclear matter are quantified
in proton-nucleus collisions,  where the energy density is not
enough to reach the critical value for the formation of QGP. There
is no anomalous suppression, beyond CNM effects, of the $J/\psi$
yield in the collisions of light ions of sulfur ($A=32$) with the
heavy ions of uranium ($A=238$). However, there is the anomalous
suppression beyond the CNM effects in the central collisions of
lead nuclei when the number of participating nucleons is the
largest, and, to a certain extent, in the central collisions of
indium nuclei ($A=115$)~\cite{Isayev:Arn}. This could be
considered as the signature of melting the charmonium states in
quark-gluon plasma. At this point, it is also worth to note that
some level of the $J/\psi$ suppression could originate with the
reduced feed-down to $J/\psi$ from excited charmonium states
($\psi',\chi_c$), which melt just above the QGP transition
temperature. One needs further to investigate  what is really
suppressed, the directly produced $J/\psi$s or the $J/\psi$s
coming from the heavier  charmonium states.

\begin{figure}[tb]\begin{center}
\includegraphics[width=6.8cm,keepaspectratio]{fig2_sps.eps}
\end{center} \caption{$J/\psi$ yields normalized by Drell-Yan
process, as a function of the nuclear thickness $L$, as measured
at the SPS.
 } \label{Isayev:fig:3}
\end{figure}

Let us consider the $J/\psi$ production at RHIC experiments, when
the nucleon-nucleon center-of-mass energy for collisions of the
heaviest ions increases by a factor 10 compared to the CERN SPS
experiments. The $J/\psi$ suppression can be characterized by a
ratio called the nuclear modification factor
\begin{equation}R_{AB}(p_{\,T},y)=
\frac{d^2N_{J/\psi}^{AB}/dp_{\,T}dy}{N_{coll}d^2N_{J/\psi}^{pp}/dp_{\,T}dy},\end{equation}
 obtained by normalizing the $J/\psi$ yield in nucleus-nucleus
collision  by the $J/\psi$ yield in $p+p$ collision at the same
energy per nucleon pair times the average number  of binary
inelastic NN collisions.  This ratio characterizes the impact of
the medium on the particle spectrum. If heavy ion collision is a
superposition of independent $N_{coll}$ inelastic NN collisions,
then $R_{AB}=1$, whereas $R_{AB}<1$ ($R_{AB}>1$) corresponds to
the case of the $J/\psi$ suppression (enhancement). As we
discussed already, at first, it is necessary to clarify the role
of  CNM effects on $J/\psi$ production. At RHIC, CNM effects are
studied in collisions of light deuteron and heavy gold nuclei,
when the energy density reached in the collision is not enough for
the formation of QGP.   At RHIC energies, shadowing of partons is
important and in the model calculations is implemented in two
shadowing schemes for the nuclear parton distribution functions,
the EKS model~\cite{Isayev:EKS} and NDSG model~\cite{Isayev:NDSG}.
The $J/\psi$ break-up cross sections obtained for two shadowing
schemes from the best fit to data are $\sigma_{\rm
breakup}=2.8^{+2.3}_{-2.1}$ mb (EKS) and $\sigma_{\rm
breakup}=2.6^{+2.2}_{-2.6}$ mb (NDSG)~\cite{Isayev:Adare} (in
Ref.~\cite{Isayev:Adare}, the term "break-up cross section" is
used  instead of the term "absorption cross section"). Although
these values are consistent, within large uncertainties, with the
corresponding value obtained at CERN SPS, a recent
analysis~\cite{Isayev:L}  shows that, in fact, the level of
$J/\psi$ CNM break-up significantly decreases with the collision
energy.

Let us now consider charmonium production in heavy ion collisions
at RHIC. Fig.~\ref{Isayev:fig:4}  shows the $p_T$-integrated
$J/\psi$ nuclear modification factor obtained in Au-Au collisions
at RHIC/PHENIX experiment as a function of centrality,
parametrized by the number of participating nucleons at mid- (left
panel) and forward (right panel) rapidity~\cite{Isayev:A}. The
PHENIX data are shown by blue symbols. The $R_{AA}$ approaches
unity for the peripheral collisions (small $N_{part}$) and goes
down to approximately 0.2 at most central collisions (large
$N_{part}$). To see the level of the anomalous suppression beyond
the cold nuclear matter effects it is necessary to extrapolate the
CNM effects obtained in d-Au collisions to the Au-Au collisions
within the given shadowing scheme and the $J/\psi$ break-up cross
section. The results are shown by black and red curves with the
corresponding error bands. It is seen that $J/\psi$ production is
significantly suppressed beyond CNM effects at forward rapidity
(right) and suppression is less pronounced at midrapidity (left)
in most central Au-Au collisions.

\begin{figure}[tb]
\begin{center}
  \begin{tabular}{cc}
  \includegraphics[bb= 0 0 567 479,width=6.8cm]{fig3a.eps} &
  \includegraphics[bb= 0 0 567 479,width=6.8cm]{fig3b.eps}
  \end{tabular}
  \end{center}
\caption{Left: $J/\psi$'s $R_{AA}$ for $Au$-$Au$ collisions at
midrapidity compared to a band of theoretical curves for the
breakup values found to be consistent with the $d$-$Au$ data. Both
EKS and NDSG shadowing schemes are included. Right: Same as on the
left but at forward rapidity.} \label{Isayev:fig:4}
\end{figure}

However, not all is still clear. Measurements of the $J/\psi$
suppression by PHENIX collaboration at RHIC lead to some
surprising features. Fig.~\ref{Isayev:fig6} shows compiled data
for the nuclear modification factor obtained in CERN SPS and RHIC
PHENIX experiments. There are two surprising results in these
measurements. First, the mid rapidity suppression in PHENIX (the
red boxes) is lower than the forward rapidity suppression (blue
boxes) despite the experimental evidence that the energy density
is higher at midrapidity than at forward rapidity, and, hence, one
could expect that at midrapidity the $J/\psi$s will be more
suppressed due to higher density of color charges. Secondly, the
nuclear modification factor $R_{AA}$ at midrapidity in PHENIX (red
boxes) and SPS (black crosses) are in agreement within error bars,
a surprising result considering that the energy density reached at
RHIC is larger than the one reached at SPS. This indicates that at
RHIC energies additional mechanisms countering the suppression,
could be operative.

\begin{figure}[tb]\begin{center}
\includegraphics[width=6.8cm]{fig4_rhic.eps}
\end{center}
\caption{$J/\psi$ nuclear modification factor for the most
energetic SPS (Pb-Pb) and RHIC (Au-Au) collisions, as a function
of the number of participants $N_{part}$. }
\label{Isayev:fig6
\end{figure}

Let us consider possible explanations of the above features. 1.
Regeneration of $J/\psi$s in the hot partonic phase from initially
uncorrelated $c$ and $\bar c$ quarks (quark coalescence model). If
to compare the $J/\psi$ suppression pattern at RHIC and SPS,
$J/\psi$s could be indeed more suppressed at RHIC than at SPS, but
then regenerated during (or at the boundary of) the hot partonic
phase from initially uncorrelated $c$ and $\bar c$ quarks. If to
compare the results at RHIC (midrapidity vs. forward rapidity), at
midrapidity, due to higher energy density, there are more $c$ and
$\bar c$ to regenerate than at forward rapidity that could explain
the stronger suppression at forward rapidity. Note that the total
number of initial $c\bar c$ pairs is larger than 10 in the most
central Au-Au collisions. Certainly, if regeneration is important
at the RHIC conditions, it will be even more important at the LHC
conditions where more than 100 $c\bar c$ pairs is expected to be
produced in the central Pb-Pb collisions. 2. $J/\psi$ production
could be more suppressed at forward rapidity due to the nuclear
shadowing effects. Standard gluon shadowing parametrizations do
not tend to produce such an effect but they are poorly constrained
by the data and further saturation effects are not
excluded~\cite{Isayev:dC}.

To show the complexity of the problem, let us consider some
theoretical models for the charmonium production, which quite
satisfactory describe the RHIC data but whose predictions for LHC
are drastically different. First, in the {\it statistical
hadronization model} (SHM)~\cite{Isayev:ABRS}, it is assumed that:
1. All heavy quarks (charm and bottom) are produced in primary
hard collisions and their total number stays constant until
hadronization. 2. Heavy quarks reach thermal equilibrium in the
QGP before the chemical freeze-out (hadronization). 3.~All
quarkonia are produced (nonperturbatively) through the statistical
coalescence of heavy quarks at hadronization. Multiplicities of
various hadrons are calculated with the grand canonical ensemble.
The generation of $J/\psi$ proceeds effectively if $c, \bar c$
quarks are free to travel over large distances implying
deconfinement.

Fig.~\ref{Isayev:fig:7} shows the rapidity dependence of the
nuclear modification factor, obtained in this model and the
comparison with the rapidity dependence at PHENIX for two
centrality bins. Two theoretical curves correspond to two fitting
procedures, with one and two Gaussians of the $J/\psi$ data in
$pp$ collisions. In both cases, calculations reproduce rather well
(considering the systematic errors) the $R_{AA}$ data. The model
describes the larger suppression away from midrapidity. The
maximum of $R_{AA}$ at midrapidity in this model is due to the
enhanced generation of charmonium around midrapidity, determined
by the rapidity dependence of the charm production cross section.
The centrality dependence of $R_{AA}$ at $y=0$ is shown in
Fig.~\ref{Isayev:fig:8}. The model reproduces quite well the
decreasing trend with centrality seen in the RHIC data.
Fig.~\ref{Isayev:fig:8} also shows the prediction of the model for
the LHC. At much higher LHC energies, the charm production cross
section is expected to be larger by about an order of magnitude.
As a result, a totally opposite trend as a function of centrality
is predicted, with $R_{AA}$ exceeding unity for central
collisions.

\begin{figure}[tb]
\begin{center}
  \includegraphics[bb= 0 0 527 307,width=10cm]{fig5.eps}
  \end{center}
  \caption{Rapidity dependence of  $R_{AA}^{J/\psi}$
for two centrality classes in the statistical hadronization model.
The data from the PHENIX experiment (symbols with errors) are
compared to calculations (lines, see text).}\label{Isayev:fig:7}
\end{figure}

\begin{figure}[tb]
\begin{center}

 \includegraphics[bb= 0 0  507 507,width=0.42\textwidth]{raa_rhic_lhc.eps}
   \end{center}
  \caption{Centrality dependence of the   $R_{AA}^{J/\psi}$
 at midrapidity, according to the statistical hadronization model.}\label{Isayev:fig:8}
\end{figure}

Let us consider {\it the comovers interaction model}
(CIM)~\cite{Isayev:CBFK}. This model does not assume the
deconfinement phase transition. Anomalous suppression of $J/\psi$
(beyond CNM effects) is the result of the final state interaction
of the $c \bar c$ pair with the dense medium produced in the
collision (comovers interaction). The model consistently treats
the initial and final state effects. The initial state effects
include: 1) nuclear absorption of the pre-resonant $c \bar c$
pairs by nucleons of the colliding nuclei, 2) consistent treatment
of nuclear shadowing for hard production of charmonium. The final
state effects include absorption of the $c \bar c$ pairs by the
dense medium created in the collision (interaction with comoving
hadrons produced in the collision). The model does not assume
thermodynamic equilibrium and, thus, does not use thermodynamic
concepts. The density of charmonium is governed by the
differential rate equation
\begin{equation}
\tau\frac{dn_{J/\psi}}{d\tau}=-\sigma_{co}[n_{co}(b,s,y)n_{J/\psi}(b,s,y)-n_c(b,s,y)
n_{\bar c}(b,s,y)],\label{Isayev:2}
\end{equation}
supposing a pure longitudinal expansion of the system and boost
invariance. In Eq.~(\ref{Isayev:2}), $n_{co}$ is the density of
comovers, which is found in the dual parton
model~\cite{Isayev:CSTT} together with the proper shadowing
correction, $\sigma_{co}$ is the cross section of $J/\psi$
dissociation due to interactions with comovers, taken such as to
reproduce the low energy SPS experimental data (with $\sigma_{co}
= 0.65$ mb). The first term on the right describes dissociation of
charmonium due to interaction with comovers. The second term
describes the recombination of charmonium and is proportional to
the product of densities of charm quarks and antiquarks. The
important feature of the CIM is that recombination of $c$-$\bar c$
quarks proceeds only locally, when the densities of quarks and
antiquarks are taken at the same transverse coordinate $s$. This
is different from the recombination in the SHM, where recombining
quarks can be separated by large distance that implies
deconfinement. The effective recombination cross section in the
CIM is equal to the dissociation cross section due to the detailed
balance. The results for the centrality dependence of the $R_{AA}$
are shown in Fig.~\ref{Isayev:fig:9}. At midrapidity, the
experimental data are well reproduced by full theoretical
calculations (solid curve) taking into account nuclear shadowing,
dissociation by comovers and recombination from charm quark and
antiquark pairs. At forward rapidity, the results also well agree
with data, in particular, the $J/\psi$ suppression at forward
rapidity is somewhat larger that the suppression at midrapidity.
\begin{figure}[tb]
\begin{center}
  \begin{tabular}{cc}
  \includegraphics[bb= 0 0 567 542,width=6.8cm]{fig7_left.eps} &
  \includegraphics[bb= 0 0 567 542,width=6.8cm]{fig7_right.eps}
  \end{tabular}
  \end{center}
\caption{Results for $J/\psi$ suppression in Au-Au collisions at
RHIC  at mid- (left panel), and at forward rapidities (right
panel) in the comovers interaction model. The solid curves are the
final results. The dashed-dotted ones are the results without
recombination ($C = 0$). The dashed line is the total
initial-state effect. The dotted line in the right panel is the
result of shadowing. In the left panel the last two lines
coincide. }
\label{Isayev:fig:9
\end{figure}
 Fig.~\ref{Isayev:fig:10} shows the
predictions of the model for LHC. The parameter $C$ encodes the
recombination from $c$-$\bar c$ pairs and is vanishing in the
absence of recombination. Although the density of charm grows
substantially from RHIC to LHC, the combined effect of
initial-state shadowing, absorption and comovers dissociation
overcomes the effect of parton recombination. This is in sharp
contrast with the predictions of the statistical hadronization
model where a strong enhancement of the $J/\psi$ yield with
increasing centrality was predicted.

\begin{figure}[tb]
\begin{center}
  \includegraphics[bb= 0 0 567 539,width=.4\textwidth]{fig8.eps}
  \end{center}
  \caption{Results for $J/\psi$ suppression in Pb+Pb at LHC ($\sqrt{s} = 5.5$ TeV) at
  midrapidity
  for different
values of the parameter $C$ in the comovers interaction model. The
upper line is the suppression due to initial-state effects
(shadowing).}\label{Isayev:fig:10}
\end{figure}

Thus, the $J/\psi$ suppression is an important characteristic to
search for QGP. But if $J/\psi$ anomalous suppression (beyond CNM
effects) was observed in HIC, there are a few competing mechanisms
to explain that: 1. Charmonium is dissociated due to the genuine
color screening in the deconfined medium. 2. Charmonium is
dissociated through interactions with comoving hadrons in the
medium formed in HIC. How to differ these mechanisms? In comovers
interaction model, the anomalous suppression sets in smoothly from
peripheral to central collisions rather than in a sudden way when
the corresponding dissociation temperature in the deconfined
medium is reached. At SPS, current experimental errors still do
not allow to disentangle these two mechanisms. However, even if
the color screening mechanism is dominating, it is unclear what is
really melted, directly produced $J/\psi$s or originating from the
feed-down of less bound charmonium states, $\chi_c$ ($\rightarrow
J/\psi+X$), $\psi'$ ($\rightarrow J/\psi+X$), which have lower
dissociation temperatures.  At RHIC and LHC conditions, the
feed-down from B-meson decays becomes also important.
Recombination enhances $J/\psi$ production and much complicates
the picture but its effect may be different depending on whether
the deconfinement phase transition happened or not. It is even
possible that after all $J/\psi$s were melted in QGP they can be
statistically regenerated at hadronization. True operative
mechanisms of $J/\psi$ production can be established only after
studying all important dependences (from centrality, rapidity,
collision energy $\sqrt{s}$,...) of all relevant observables with
sufficient accuracy.

So far we considered the suppression (or enhancement) patterns for
the $J/\psi$ production. Now let us briefly consider the
perspective for bottomonia. One could expect that bottomonium
production might be easier to understand than charmonium
production due to the following reasons. 1. Since less than one
$b\bar b$ pair is produced in one central Au-Au collision, the
regeneration is negligible at the conditions of RHIC. Besides,
only about 5 $b\bar b$ pairs are expected to be produced in a
single central Pb-Pb collision at LHC. Hence, regeneration should
play much less role in the beauty sector than in the charm sector.
2. Having higher masses, bottomonia originate from higher momentum
partons and will less suffer from shadowing effects. 3. The
absorption cross section for $\Upsilon$ is by $40-50\,\%$
 smaller than the corresponding cross
section for $J/\psi$ and $\psi'$.

These features should ease the separation of the anomalous
suppression in the $\Upsilon$'s family.

\section{Heavy Flavor Probes of QGP: Open Heavy Flavor Mesons}

What qualitative effects could one expect to obtain when probing
the dense matter by heavy quarks (charm or bottom)? As  well known
from electrodynamics, the bremsstrahlung off an accelerated heavy
quark $Q$ is suppressed by the large power of its mass $\sim
(m_q/m_Q)^4$
 as compared to light quarks. Therefore,
gluon radiation off heavy quarks (i.e., radiative energy loss) is
much suppressed relative to light quarks. Consequently, one could
expect a decrease of high $p_T$ suppression and of the elliptic
flow coefficient $v_2$~\cite{Isayev:I} from light to charm to
bottom quarks. Or, that the energy loss and coupling to matter of
heavy quarks is smaller than for light quarks as well as that the
thermalization time for heavy quarks is longer than for light
quarks. Due to the above features, one should observe a pattern of
gradually increasing $R_{AA}$ when going from the mostly
gluon-originated light-flavor hadrons ($h^\pm$ and $\pi^0$) to $D$
  to $B$ mesons: $R_{AA}^h\lesssim R_{AA}^D\lesssim
R_{AA}^B$~\cite{Isayev:Ar} (gluons lose more energy than quarks
since gluons have a higher color charge). The enhancement above
the unity of the heavy-to-light ratio
$R_{AA}^{D/h}=R_{AA}^D/R_{AA}^h$ probes the color charge
dependence of the parton energy loss while the ratio
$R_{AA}^{B/D}=R_{AA}^B/R_{AA}^D$ probes the mass dependence of the
parton energy loss.

Let us now consider what the experiment tells us about the open
heavy flavor $p_T$ suppression and elliptic flow. As mentioned
earlier, at RHIC, open heavy flavors are studied through
measurements with the electrons and positrons originating from the
semileptonic decays of $D$- and $B$-mesons.
Fig.~\ref{Isayev:fig:11} shows the nuclear modification factor and
elliptic flow coefficient for HF electrons as functions of $p_T$,
obtained in central collisions of gold nuclei at RHIC (closed
circles)~\cite{Isayev:A2}-\cite{Isayev:R}. In contrast to the
above expectations, the results for $R_{AA}^{HF}$ show a strong
suppression of HF decay electrons at $p_T>2$ GeV/c, approaching at
high $p_T$ the level of suppression for $\pi^0$. This evidences
that produced medium is quite dense for heavy quarks to lose
energy as efficiently as light quarks do. The measurement of
elliptic flow gives rather large value for $v_2^{HF}$. This means
that HF electrons are involved in a collective motion being
indicative of the collective flow of their parent particles as
well.

\begin{figure}[tb]
\begin{center}
  \includegraphics[bb= 10 15 550 545,width=.4\textwidth]{fig9.eps}
  \end{center}
  \caption{Nuclear modification factor (upper panel, central Au-Au) and elliptic flow
(lower panel, minimum-bias Au-Au) of non-photonic electrons at
RHIC, compared to theory. The band corresponds to the Langevin
simulations based on an expanding fireball with effective heavy
quark resonance interactions~[22].} \label{Isayev:fig:11}
\end{figure}

In Fig.~\ref{Isayev:fig:11}, the results of some model
considerations are also shown. The best description is provided by
the model assuming the Brownian motion of heavy quarks within the
framework of Langevin dynamics~\cite{Isayev:HGR}. Let us consider
the basic assumptions of this theory. Firstly, the thermal heavy
quark momentum $p^2\sim mT$ (in nonrelativistic approximation) is
much larger than the typical momentum transfer $Q^2\sim T^2$ from
a thermal medium to a heavy quark. Hence, the motion of a heavy
quark in the QGP can be represented as the Brownian motion, which
can be described using the  Langevin equation. Secondly, heavy
quark loses its energy in elastic scattering processes with light
partons. Besides, as evidenced from calculations of heavy and
light meson correlators within LQCD, in QGP the $D$- and $B$-meson
like resonant states exist up to $T<2T_c$. Rescattering on these
resonant states plays an important role in thermalizing heavy
quarks. Thirdly, in order to get the spectrum of HF electrons,
$c$- and $b$-quarks are to be hadronized into $D$- and $B$-mesons
via quark coalescence (at low $p_T$) and fragmentation (at high
$p_T$).

The analysis shows  that resonance scattering decreases nuclear
modification factor $R_{AA}^{HF}$ and increases azimuthal
asymmetry $v^{HF}_2$. Heavy-light quark coalescence in subsequent
hadronization significantly amplifies $v^{HF}_2$ and increases
$R_{AA}^{HF}$, especially in the $p_T\simeq2$ GeV/c region. The
contribution from B-mesons to $R_{AA}^{HF}$ and $v^{HF}_2$ is
estimated providing full calculations with $c$+$b$ quarks and with
only $c$ quarks. The result is that the B-meson contribution
increases $R_{AA}^{HF}$ and decreases $v^{HF}_2$, and becomes
important above $p_T\simeq3$ GeV/c. Thus, one can conclude, that
the combined effects of coalescence of heavy quarks $Q$ with light
quarks $q$, and of the resonant heavy-quark interactions are
essential in generating strong elliptic flow $v^{HF}_2$ of up to
$10\%$, together with strong suppression of heavy flavor electrons
with $R^{HF}_{AA}$ about 0.5.

\section*{References}

High Energy Physics - Phenomenology (hep-ph)


The evolution of C and O abundances in stellar populations

Abstract: Carbon and oxygen abundances in F and G main-sequence stars ranging in metallicity from [Fe/H] = -1.6 to +0.5 are determined from a non-LTE analysis of CI and OI atomic lines in high-resolution spectra. Both C and O are good tracers of stellar populations; distinct trends of [C/Fe] and [O/Fe] as a function of [Fe/H] are found for high- and low-alpha halo stars and for thick- and thin-disk stars. These trends and that of [C/O] provide new information on the nucleosynthesis sites of carbon and the time-scale for the chemical enrichment of the various Galactic components.

\section{Introduction}

Recent $\Lambda$CDM, hydrodynamical simulations of the formation of the Galaxy
(e.g. Zolotov et al. \cite{zolotov10}; Font et al. \cite{font11}; McCarthy et al.
\cite{mccarthy12})
predict the existence of two populations of halo stars. The first one is formed
{\em in situ} in merging and dissipating gas clumps, whereas the other one is
{\em accreted} from
satellite galaxies. For the upper end of the halo metallicity distribution,
Zolotov et al. (\cite{zolotov10}) find that \alphafe\
in the accreted population decreases with increasing [Fe/H] relative to the
near-constant \alphafe\ in the in-situ population.
This is due to a difference in star formation rate (SFR);
the chemical enrichment proceeds at a slower rate in satellite galaxies, so that
Type Ia SNe start contributing with iron at a lower metallicity.

Evidence for two halo populations\footnote{Selected as stars having a
total space velocity $V_{\rm total} >  180$\,\kmprs\ relative to the
local standard of rest (LSR).} in the solar neighborhood with different
\alphafe\ trends has been found by Nissen \& Schuster (\cite{nissen10}, hereafter NS10).
The ``high-alpha" stars have $\alphafe \simeq 0.3$, similar to
thick-disk stars, and extend in metallicity up to $\feh \simeq -0.4$.
The ``low-alpha" stars show a declining \alphafe -trend from $\sim 0.3$\,dex at
$\feh \simeq -1.6$ to $\sim 0.1$\,dex at $\feh \simeq -0.8$, which is the maximum
\feh\ reached by this population. Furthermore, the Toomre diagram shows that
the high-alpha stars tend to move on prograde orbits, whereas the low-alpha
stars have higher space velocities with respect to the LSR
and an excess of retrograde orbits. These data are consistent with a
scenario where the high-alpha stars have been formed in dissipational
mergers of gas clouds, and the low-alpha stars have been accreted
from satellite galaxies.

The high- and low-alpha halo stars also separate in \nafe , \nife , \cufe ,
and \znfe , but not in \mnfe\ (Nissen \& Schuster \cite{nissen11}).
Furthermore, the high-alpha stars seem to be on average 2 - 3\,Gyr older
than the low-alpha stars (Schuster et al. \cite{schuster12}).

In this paper we report on carbon and oxygen abundances in the two halo populations
and compare with abundances in stars having thin- and thick-disk kinematics.

\section{Carbon abundances}

The abundance of carbon was determined from equivalent widths (EWs) of the
high-excitation
\CI\ lines at 5052.2 and 5380.3\,\AA\ as measured in high-resolution,
high-S/N VLT/UVES and NOT/FIES spectra. For each star, a 1D model atmosphere
was obtained from the MARCS grid (Gustafsson et al. \cite{gustafsson08})
by interpolating to the stellar values of \teff , \logg , \feh , and \alphafe .
The model was used to to derive an LTE carbon abundance. Non-LTE corrections from
Takeda \& Honda (\cite{takeda05}) were afterwards applied, but they are small,
i.e. less than 0.02\,dex.

The atmospheric parameters of the stars were determined spectroscopically
by analysing \FeI\ and \FeII\ lines relative to the same lines in the spectra of
two nearby thick-disk stars, HD\,22879 and HD\,76932, for which \teff\ and \logg\
are well known from colour indices and Hipparcos parallaxes. This enable us to
determine very precise differential parameters and abundances for stars
that belong to the same region of the HR-diagram as the standard stars.
The absolute values may be more uncertain. Thus we have increased the \teff\ scale
by +100\,K relative to the scale used in NS10
to take into account the new IRFM \teff -- colour calibration by
Casagrande et al. (\cite{casagrande10}). This has only a small effect on
\logg , \feh\ (derived from \FeII\ lines), and \alphafe\ (derived from
neutral lines), but it has a significant effect on C and O abundances
derived from high-excitation atomic lines.

\begin{figure}[t]
\begin{center}
 \includegraphics[width=9cm]{Nissen_PE_fig1.eps}
 \caption{[C/Fe] versus [Fe/H] with the same classification in high- and
low-alpha halo and thick-disk stars as in NS10
and with thin-disk stars from Nissen (\cite{nissen13}) included.}
   \label{Nissen:fig1}
\end{center}
\end{figure}

The derived \cfe\ values are shown in Fig. \ref{Nissen:fig1}.
In addition to stars from NS10, we have included  stars with thin-disk kinematics
from Nissen (\cite{nissen13}), who also used
the $\lambda \lambda 5052, 5380$ \CI\ lines (measured in high-quality
HARPS spectra) to derive C abundances. The error bars shown were
estimated by adding, in quadrature, the errors arising from the EWs
and those corresponding to the internal
uncertainties of \teff\ and \logg . At low metallicity and low \alphafe ,
the \CI\ lines are very weak, i.e. $EW \sim$ 2 - 4\,m\AA\ only, so
the error of \cfe\ due to the uncertainty of the EW measurements
becomes large.

\section{Oxygen abundances}

Oxygen abundances were determined from equivalent widths of the $\lambda 7774$
\OI\ triplet lines. The NOT/FIES spectra do not cover this wavelength region
and only a few of the UVES spectra have the \OI\ triplet
included. We have therefore adopted EW measurements from
Ram\'{\i}rez et al. (\cite{ramirez12}), who derived oxygen abundances
for a subset of stars from NS10 using high-resolution spectra obtained
with the 2.7\,m telescope at the McDonald observatory, Keck/HIRES, and
the MIKE spectrograph at the Magellan Telescope. For the thin-disk stars
we used FEROS spectra obtained with the ESO 2.2\,m telescope. Altogether,
\OI\ triplet data are available for 101 stars out of the 117 stars for which
\CI\ abundances were determined.

In deriving O abundances from the 7774\,\AA\ triplet, non-LTE corrections
from Fabbian et al. (\cite{fabbian09a}) were applied as described
by Nissen  (\cite{nissen13}). These corrections are important, even
for differential determinations in
the \teff -range of our sample of stars; the correction of [O/H]
ranges from about +0.1\,dex at 5300\,K to $-0.1$\,dex at 6300\,K.

The derived \ofe\ values are shown in Fig. \ref{Nissen:fig2}.
The distribution of stars looks much the same as that obtained by
Ram\'{\i}rez et al. (\cite{ramirez12}, Fig. 1, lower panel), except
that they find a flatter distribution of [O/Fe] for the
thick-disk and high-alpha halo stars. This is due to a somewhat stronger dependence
of their non-LTE corrections on [Fe/H] than those of  Fabbian et al. (\cite{fabbian09a}).

\begin{figure}[t]
\begin{center}
 \includegraphics[width=9cm]{Nissen_PE_fig2.eps}
 \caption{[O/Fe] versus [Fe/H].}
   \label{Nissen:fig2}
\end{center}
\end{figure}

\section{Discussion}

As seen from Figs. \ref{Nissen:fig1} and \ref{Nissen:fig2},
the low-alpha halo stars are separated in
\cfe\ and \ofe\ from thick-disk and high-alpha halo stars.
For the latter two populations, the dispersion in \cfe\ and \ofe\
at a given metallicity can be explained in terms of the errors of
the derived abundances, whereas there seems to be an additional cosmic
dispersion in \cfe\ and \ofe\ for the low-alpha halo population.
In support of a real cosmic scatter,
Fig. \ref{Nissen:fig3} shows \cfe\ versus \nafe\ (from NS10) for
the metallicity  range $-1.2 < \feh < -0.7$, where the largest
dispersion is present. As seen, there is a striking correlation
between \cfe\ and \nafe\ except for two strongly deviating stars.

As suggested by the $\Lambda$CDM simulations of Zolotov et al. (\cite{zolotov10}),
a possible explanation of the abundance differences seen
in these figures is that the high-alpha stars were born in
the innermost part of the Galaxy
in a deep gravitational potential with such a high SFR
that Type Ia SNe did not contribute significantly with iron
until a metallicity of $\feh \simeq -0.4$.
Later these stars were dispersed to the halo by merging satellite galaxies.
The low-alpha stars, on the other hand, were formed in
dwarf galaxies with a relatively shallow potential, and hence low SFR,
so that Type Ia SNe started contributing with iron at a metallicity
around or below $\feh = -1.6$. The reason for the scattter
in \cfe , \ofe , \nafe , and \alphafe\ at a given metallicity could then be that
the various satellite galaxies from which the low-alpha stars were accreted
had different masses and therefore different star formation rates.

\begin{figure}[t]
\begin{center}
 \includegraphics[width=9cm]{Nissen_PE_fig3.eps}
 \caption{[C/Fe] versus [Na/Fe] for stars with $-1.2 < \feh < -0.7$.}
   \label{Nissen:fig3}
\end{center}
\end{figure}

The two deviating low-alpha stars in Fig. \ref{Nissen:fig3},
G\,53-41 and G\,150-40, require a special explanation. They
have high \nafe\ but low \cfe\ and also very low \ofe . Hence,
they share the so-called Na-O  anti-correlation in second generation
stars in globular clusters, which are thought to be made of gas
polluted by first-generation AGB stars in which the Ne-Na cycle
has occured. This suggests that some of the halo stars originate
from disrupted globular clusters, as also noted by Ram\'{\i}rez et al. (\cite{ramirez12}).

Figures \ref{Nissen:fig1} and \ref{Nissen:fig2} also give a hint
for a difference in \cfe\ and \ofe\ between thin- and
thick-disk stars in the overlapping metallicity range
$-0.7 < \feh < -0.2$ although we have only two thin-disk stars
in this range. A systematic difference is well established
in the case of \ofe\ (e.g. Bensby \& Feltzing \cite{bensby06}, Fig. 11b)
and have also been found for the alpha-capture elements,
Mg, Si, Ca, and Ti (e.g. Adibekyan \cite{adibekyan12}, Fig. 8),
but it was not seen in the case of \cfe\ by
Bensby \& Feltzing (\cite{bensby06}, Fig. 11a), Clearly, we need
more thin-disk stars in the overlapping metallicity range to verify
the possible difference in \cfe\ between the thin- and thick-disk
populations.

Fig. \ref{Nissen:fig4} shows \co\ versus \oh\ for
the four populations studied in this paper with thin- and
thick-disk stars from Bensby \& Feltzing (\cite{bensby06}) included.
Their C and O abundances are derived from the weak forbidden
$\lambda 8727$ [\CI ] and $\lambda 6300$ [\OI ] lines but agree
very well with our abundances derived from high-excitation
\CI\ and \OI\ lines.

\begin{figure}[t]
\begin{center}
 \includegraphics[width=9cm]{Nissen_PE_fig4.eps}
 \caption{[C/O] versus [O/H]. Thin- and thick-disk stars from
  Bensby \& Feltzing (\cite{bensby06}) are included without error bars.}
   \label{Nissen:fig4}
\end{center}
\end{figure}

As seen from Fig. \ref{Nissen:fig4}, there is no systematic
shift in \co\ between high- and low-alpha halo stars.
This is perhaps surprising, because the time-scale for
the chemical enrichment of the low-alpha population is
long enough to allow Type Ia SNe to contribute with iron
(our explanation for the low $\alpha$/Fe ratios),
and one could therefore have expected that low- and
intermediate-mass AGB stars had enough time to contribute with
carbon and raise \co\ in the low-alpha stars
to higher values than in high-alpha stars.
The explanation may be that intermediate-mass (4 - 8 $M_{\rm Sun}$)
AGB stars contribute very little to $^{12}$C
(Kobayashi et al.  \cite{kobayashi11}), and
that the evolution time-scale of low-mass (1 - 3 $M_{\rm Sun}$) AGB stars
(which do have a high
$^{12}$C yield according to Kobayashi et al.) is longer than the
chemical enrichment time-scale of the low-alpha
population. We conclude that carbon in both high- and low-alpha
halo stars was made primarely in high-mass stars ($> 10 \,  M_{\rm Sun}$),
but with a metallicity-dependent yield to account
for the slight increase of \co\ with increasing \oh .

\begin{figure}[b]
\begin{center}
 \includegraphics[width=12cm]{Nissen_PE_fig5.eps}
 \caption{Part of the CH-band in a high- and a low-alpha star with similar
atmospheric parameters. All unmarked lines are due to CH.}
   \label{Nissen:fig5}
\end{center}
\end{figure}

As also seen from Fig. \ref{Nissen:fig4}, thick-disk stars
have on average the same \co\ as high-alpha halo stars,
but thin-disk stars have systematically higher \co\
in the overlapping metallicity range $-0.5 < \oh < +0.1$.
This offset in  \co\ between thin- and thick-disk stars
was already detected by Bensby \& Feltzing (\cite{bensby06}, Fig. 12),
but is even more clearly seen in our Fig. \ref{Nissen:fig4}.
Probably, the higher \co\ ratios in the Galactic thin disk
are due to low-mass AGB stars, i.e. the
chemical enrichment time-scale of the thin disk
is long enough to allow low-mass stars to contribute.

As a final remark, we note that the $\lambda \lambda 5052, 5380$ \CI\ lines
used to derive C abundances in this paper are very weak in metal-poor stars.
Thus, very high S/N, high-resolution spectra are required to derive precise
C abundances. The somewhat stronger \CI\ lines in the 9000 - 9500\,\AA\ region may
be used instead (Fabbian et al. \cite{fabbian09b}) but they are blended
by telluric lines. As an interesting alternative one may use the CH band.
Fig. \ref{Nissen:fig5} shows part of this band
in NOT/FIES spectra of a high-alpha star, G\,31-55 (\teff , \logg , \feh , \alphafe)
= (5738\,K, 4.33, $-$1.10, 0.29) and a low-alpha star HD\,193901 (5756\,K, 4.39,
$-$1.09, 0.16). As seen, the three \FeI\ lines have similar strengths in the
two stars, whereas all CH lines are significantly weaker in the low-alpha star.
Hence, it may be possible to use the CH band
in large forthcoming surveys such as HERMES/GALAH and
ESO/4MOST for high precision differential studies of carbon
abundances in stellar populations.

\begin{acknowledgements}
Funding for the Stellar Astrophysics Centre is provided by the
Danish National Research Foundation (Grant agreement no.: DNRF106).
This paper is based on observations made with the Nordic Optical Telescope
on La Palma, and on data from the European Southern Observatory
ESO/ST-ECF Science Archive Facility.
\end{acknowledgements}

Astrophysics of Galaxies (astro-ph.GA); Solar and Stellar Astrophysics (astro-ph.SR)


Inverted orbital polarization in strained correlated oxide films

Abstract: Manipulating the orbital occupation of valence electrons via epitaxial strain in an effort to induce new functional properties requires considerations of how changes in the local bonding environment affect the band structure at the Fermi level. Using synchrotron radiation to measure the x-ray linear dichroism of epitaxially strained films of the correlated oxide CaFeO3, we demonstrate that the orbital polarization of the Fe valence electrons is opposite from conventional understanding. Although the energetic ordering of the Fe 3d orbitals is confirmed by multiplet ligand field theory analysis to be consistent with previously reported strain-induced behavior, we find that the nominally higher energy orbital is more populated than the lower. We ascribe this inverted orbital polarization to an anisotropic bandwidth response to strain in a compound with nearly filled bands. These findings provide an important counterexample to the traditional understanding of strain-induced orbital polarization and reveal a new method to engineer otherwise unachievable orbital occupations in correlated oxides.

\section{Introduction}

The use of epitaxial strain to induce occupation of specific electron orbitals by removing orbital degeneracies has been pursued in transition metal oxides in an effort to engineer new electronic and magnetic properties \cite{Tokura_Manganites_XLD, Aruta_LSMO_XLD, Strain_OP_Csiszar, Hansmann_Held_theory_nickelate_SL_OP, Nickelate_holes_reduced_OP, Chak_asymmetric_XLD, Freeland_LNO_XLD, Wu_nickelate_SL_XLD, Wu_PrNiO3_orbital_polarization, Bruno_nickelate_XLD, Pesquera_LSMO_XLD}. Such strain-induced orbital polarization has been very successfully described by ligand field theory, which considers the overlap of electron orbitals between a central cation and its surrounding anions \cite{Cox_oxide_book, Khomskii_book}. For transition metal perovskite oxides, the metal cation is octahedrally coordinated by six oxygen anions, or ligands. This $O_h$ symmetry splits the five degenerate $d$-levels into two groups: a lower, triply degenerate group ($t_{2g}$) and a doubly degenerate group ($e_g$) higher in energy by an amount $10Dq$. Whereas the lobes of the O $p$ orbitals point in between the $t_{2g}$ lobes, they directly overlap with the $e_g$ lobes, which comes at a coulombic energy cost that raises the $e_g$ orbitals in energy. Epitaxial strain alters the local crystal field and lifts the $t_{2g}$ and $e_g$ degeneracies. For example, tensile strain reduces the overlap between the $e_g$ orbital of $d_{x^2-y^2}$ symmetry and its ligands, thus lowering its energy relative to the other $e_g$ orbital, $d_{3z^2-r^2}$, by an amount $\Delta e_g$ [see Fig. 1(a) inset]. Unless the $e_g$ orbitals are fully filled, one subsequently expects $d_{x^2-y^2}$ to be preferentially occupied; the converse applies for compressive strain. This simple picture has been used to explain strain-induced orbital polarization in many systems, particularly \textit{AB}O\textsubscript{3} perovskite oxides \cite{Tokura_Manganites_XLD, Aruta_LSMO_XLD, Hansmann_Held_theory_nickelate_SL_OP, Nickelate_holes_reduced_OP, Chak_asymmetric_XLD, Freeland_LNO_XLD, Wu_nickelate_SL_XLD, Wu_PrNiO3_orbital_polarization, Bruno_nickelate_XLD, Pesquera_LSMO_XLD}. In this Letter, we find that this model fails to explain orbital polarization in strained films of CaFeO\textsubscript{3}, which exhibit orbital polarization opposite to that described above.

To quantify the electron occupation of specific $e_g$ orbitals, we measure x-ray absorption across the Fe $L$- and O $K$-edge resonance energies using linearly polarized photons, which allows us to differentiate between $d_{x^2-y^2}$ and $d_{3z^2-r^2}$ occupations. Analyzing the x-ray linear dichroism using multiplet ligand field theory reveals that the effect of epitaxial strain on the energetic ordering of the $e_g$ orbitals is consistent with the aforementioned considerations--stretched bonds are lower in energy than unstretched. Given this energetic landscape, however, the expected orbital \textit{occupations} do not follow: The out-of-plane ($d_{3z^2-r^2}$) orbitals are more populated under tensile strain (and \textit{vice versa} for compressive strain). We propose that this inverted orbital polarization arises from strain-induced anisotropic changes in the Fe-O-Fe bond angles and the resulting anisotropic bandwidths in bands that are more-than-half-filled. Such conditions are not limited to ferrates but could arise in other strongly hybridized systems, such as the rare-earth nickelates \cite{Sawatzky_bond_disproportionation}.

 \begin{figure}
 \includegraphics{Fig1.eps}
 \caption{Polarization-dependent x-ray absorption measured by total electron yield (TEY) across the Fe $L$-edge for CaFeO\textsubscript{3} under (a) tensile strain and (b) compressive strain. Inset:  Octahedral crystal field splitting of transition metal $d$ levels for a (001)-oriented film under no strain ($c=a$) and under biaxial tensile strain ($c<a$). \label{Fig1_Fe_XLD}}
 \end{figure}

CaFeO\textsubscript{3} films of 40 pseudocubic unit cells (${\sim}$15 nm thick) were deposited by oxygen plasma-assisted molecular beam epitaxy. Epitaxial strain was achieved by deposition on single crystal, (001)-oriented substrates: YAlO$_3$ (YAO, -2.0\% strain), SrLaAlO$_4$ (SLAO, -0.7\%), LaAlO$_3$ (LAO, 0.2\%), (La$_{0.18}$Sr$_{0.82}$)(Al$_{0.59}$Ta$_{0.41}$)O$_3$ (LSAT, 2.3\%), and SrTiO$_3$ (STO, 3.3\%). As previously reported, the films are coherently strained and exhibit bulk-like electrical transport, indicating high-quality, stoichiometric films \cite{Rogge_PRM}. Prior to all measurements, the films were reoxidized by heating to ${\sim}$600 $^\circ$C in oxygen plasma (200 Watts, 1x$10^{-5}$ Torr chamber pressure) and then slowly cooled to room temperature in oxygen plasma. X-ray absorption spectroscopy was performed at the Advanced Light Source, Beamline 4.0.2 and at the National Synchrotron Light Source-II, Beamline 23-ID-1. The spectra were recorded at 290 K, where CaFeO\textsubscript{3} is paramagnetic with metallic conductivity \cite{Woodward_CFO}. The x-ray incident angle was $20^{\circ}$ from the film plane, and a geometric correction was applied to the absorption measured with photons polarized out of the film plane \cite{SI_XLD}.

Although CaFeO\textsubscript{3} has an unusually high formal oxidation state of Fe\textsuperscript{4+}, its ground state exhibits a significant self-doped ligand hole density due to its negative charge transfer energy, $\Delta$ \cite{Kawasaki_CFO_first_transport, Bocquet_SFO_ligand_holes, Woodward_CFO, Matsuno_CFO_dispro, Takeda_CFO, Rogge_PRM}. In this regime, the transition metal cation does not adopt its formal oxidation state but instead keeps an extra electron that results in a hole ($\underline{L}^1$) on the oxygen ligand \cite{ZSA, Sawatzky_neg_charge_trans_1, Matsuno_CFO_dispro}. So while CaFeO\textsubscript{3} has a nominal Fe configuration of $d^4$ ($e_g^1$), its ground state has a strong $d^5\underline{L}^1$ contribution. Because of the half-filled $d$-shell, this $d^5\underline{L}^1$ ($e_g^2$) state has no significant orbital polarization and is expected to decrease the degree of orbital polarization achievable in the Fe states.

X-ray absorption across the Fe $L$-edge for a CaFeO\textsubscript{3} film under tensile strain is shown in Fig. \ref{Fig1_Fe_XLD}(a). The $L_3$ peak exhibits primarily a single, broad peak (with a small shoulder) that is consistent with nominal Fe\textsuperscript{4+} \cite{Abbate_SFO_XAS, Reduced_SFO_XAS} and significantly contrasts with the well-separated double peak structure seen in Fe\textsuperscript{3+} perovskites, such as LaFeO\textsubscript{3} and EuFeO\textsubscript{3} \cite{Reduced_SFO_XAS, Amber_EuFeO3}. This spectral signature as well as the bulk-like electrical transport indicate that oxygen vacancies have been sufficiently suppressed. As seen in Fig. \ref{Fig1_Fe_XLD}(a), the x-ray absorption is polarization-dependent. The difference in absorption measured with photons polarized parallel to the film plane, $I_x$, and photons polarized out of the film plane, $I_z$, is termed x-ray linear dichroism (XLD = $I_x-I_z$). The XLD shows areas of both positive and negative intensity, and this lineshape is similar but nearly opposite in sign for the compressively strained film, CaFeO\textsubscript{3}/SLAO, shown in Fig. \ref{Fig1_Fe_XLD}(b). Because $I_x$ preferentially probes empty states in $d_{x^2-y^2}$ and $I_z$ probes $d_{3z^2-r^2}$, their difference in total integrated intensity is a measure of the orbital polarization \cite{Thole_dichroism_prediction, van_der_Laan_XLD}, and indeed the XLD integrals are non-zero.

Evaluating the sign of the integrated XLD, however, uncovers a surprising result: the $e_g$ electron occupation does not follow the conventional ligand field model. For tensile strain the positive XLD integral implies more empty $d_{x^2-y^2}$ states. Thus under tensile strain CaFeO\textsubscript{3} has more electrons in $d_{3z^2-r^2}$, which is opposite of that predicted by ligand field theory. Under compressive strain, the integrated XLD sign implies that $d_{x^2-y^2}$ has more electrons. This behavior is consistent among the other films: The integrated XLD for tensile CaFeO\textsubscript{3} on STO (+3.3\% strain) is positive, compressed CaFeO\textsubscript{3} on YAO (-2.0\%) is negative, and the relatively unstrained CaFeO\textsubscript{3} film on LAO (+0.2\%) is approximately zero \cite{SI_XLD}.

In order to verify these relative $e_g$ occupations, we repeated the XLD measurements at the O $K$-edge. This transition probes unoccupied states with O $2p$ character, which are strongly hybridized with Fe $3d$ states due to the negative charge transfer energy \cite{Abbate_SFO_XAS}. Because these ligand states have the same symmetry as the Fe $3d$ states that they hybridize with \cite{Sawatzky_bond_disproportionation}, they are expected to mimic the Fe $e_g$ occupation. We particularly focus on the O $K$-edge prepeak feature between 526 and 529 eV because it directly probes the oxygen ligand hole states \cite{Abbate_SFO_XAS, Chen_cuprate_O_prepeak, Suntivich_O_Kedge_holes, Pellegrin_holes_prepeak}. We note that oxygen in the substrates contributes only at energies above the prepeak. As seen in Fig. \ref{O_XAS}, the oxygen prepeak exhibits linear dichroism, where the tensile strained film, CaFeO\textsubscript{3}/LSAT, has positive dichroism and the compressively strained film, CaFeO\textsubscript{3}/SLAO, has negative dichroism. A positive integrated XLD indicates more empty states in the $p_x$ and $p_y$ orbitals compared to $p_z$--that is, under tensile strain, more electrons have $p_z$ character than $p_x$ and $p_y$; the opposite situation exists for the film under compressive strain. This precisely mirrors the $e_g$ occupation measured for the Fe $3d$ states. The other strained films (CaFeO\textsubscript{3}/STO, CaFeO\textsubscript{3}/LAO, CaFeO\textsubscript{3}/YAO) exhibit an O prepeak XLD consistent with the two films highlighted here \cite{SI_XLD}.

 \begin{figure}
 \includegraphics{Fig2.eps}
 \caption{Polarization-dependent x-ray absorption of the O $K$-edge prepeak (arrow) for CaFeO\textsubscript{3} under (a) tensile and (b) compressive strain measured by total fluorescence yield. \label{O_XAS}}
 \end{figure}

With the qualitative occupation of the Fe $e_g$ orbitals confirmed, we now quantitatively estimate the orbital polarization by computing the hole ratio, $\int I_x dE$/$\int I_z dE$ \cite{Thole_hole_ratio, Haverkort_dissertation}. The hole ratio depends on the Fe valence filling, and for high-spin CaFeO\textsubscript{3} there are three $t_{2g}$ electrons and, as will be shown below, we find that the total $e_g$ occupation is 1.85 electrons. For CaFeO\textsubscript{3}/LSAT (tensile), the hole ratio is 1.018, which gives 0.90 electrons in $d_{x^2-y^2}$ and 0.95 electrons in $d_{3z^2-r^2}$ \cite{SI_XLD}, or an orbital polarization of ${\sim}$6\%. Repeating for compressively strained CaFeO\textsubscript{3}/SLAO, we find 0.93 electrons in $d_{x^2-y^2}$ and 0.91 electrons in $d_{3z^2-r^2}$, or ${\sim}$2\% polarized. This smaller orbital polarization is consistent with its lower strain state (-0.7\%) compared to CaFeO\textsubscript{3}/LSAT (+2.3\%).

To help interpret these findings, we analyze the x-ray absorption for CaFeO\textsubscript{3} using multiplet ligand field theory of a FeO\textsubscript{6} cluster \cite{Groot_XAS_multiplet_review}. We begin with the formal Fe$^{4+}$ ($3d^4$) occupation and full ligand orbitals while including a negative charge transfer energy \cite{Rogge_PRM} such that the configuration interaction ground state has primarily a $d^5\underline{L}^1$ character but still exhibits the same $S=2$ high spin symmetry of the $3d^4$ case \cite{Sawatzky_bond_disproportionation}. This $S=2$ configuration has a two-fold degeneracy due to the hybridized $e_g$ orbitals that is lifted by the imposed strain (corresponding to preferential occupation of $d_{x^2-y^2}$ and preferential occupation of $d_{3z^2-r^2}$, respectively). Further, each of these $S=2$ states has a five-fold spin degeneracy, which is lifted via the atomic spin-orbit interaction and non-tetragonal local crystal field distortions. We neglect the latter and hence label the spin-orbit split states by $J_z = 0,\pm 1,\pm 2$. Thus our XAS and XLD spectra are expected to be linear combinations of two sets of 5 spectra, one set corresponding to preferential $d_{x^2-y^2}$ occupation and one for preferential $d_{3z^2-r^2}$ occupation \cite{Wu_nickelate_SL_XLD}. The model parameters were optimized by comparing the calculated XLD to the experimental XLD \cite{SI_XLD}.

The XLD from these two sets of five calculated spectra are shown in Fig. \ref{Fig_Jz}(a) for moderate tensile strain ($\Delta e_g = +40$ meV). At finite temperature, the experimental spectrum is expected to be a combination of these $J_z$ spectra \cite{Haverkort_LaTiO3_Jz_boltzmann}, depending on their relative energies due to the spin-orbit splitting and low symmetry crystal field distortions. Therefore, a least-squares fitting procedure was used to determine a coefficient value for each of the $J_z$ XLD spectra such that the resulting combination produces the best fit with experiment. This calculated XLD spectrum has a corresponding x-ray absorption spectrum, and the XLD fitting was constrained such that the resulting calculated x-ray absorption spectral weight ($I_x + I_z$) is within $\pm1$\% of the experimental spectral weight.

 \begin{figure*}
 \includegraphics{Fig3.eps}
 \caption{(a) Linear dichroism spectra calculated using multiplet ligand field theory for a FeO$_6$ cluster with preferential $d_{x^2-y^2}$ occupation and preferential $d_{3z^2-r^2}$ occupation under moderate tensile strain ($\Delta e_g = +40$ meV; $\Delta = -2.0$ eV). Because the $J_z$ doublets exhibit nearly identical spectra [\textit{i.e.}, $(J_z = +1) \equiv (J_z = -1$)], their averaged spectrum is shown, reducing the number of XLD spectra from 10 to six. A combination of the calculated XLD spectra were fit to the experimental XLD for CaFeO\textsubscript{3} under (b) tensile and (c) compressive strain. (d) From the best XLD fit for CaFeO\textsubscript{3}/LSAT, the resulting x-ray absorption spectrum is compared to experiment. (e) The $\chi^2$ value (filled squares) for the $J_z$ fit to CaFeO\textsubscript{3}/LSAT is minimized for negative values of $\Delta$. The total number of $e_g$ electrons (open circles) increases as $\Delta$ decreases. (f) The best XLD fit occurs for $\Delta e_g > 0$ for tensile strained CaFeO\textsubscript{3}/LSAT and for $\Delta e_g < 0$ for compressively strained CaFeO\textsubscript{3}/SLAO. (g) Simplified schematic of the proposed effect of anisotropic bandwidths on the resulting orbital polarization under tensile strain ($\Delta e_g>0$) for a band with greater than half-filling (Fermi level indicated by the vertical, dashed line; filled states are shaded). \label{Fig_Jz}}
 \end{figure*}

As seen in Figs. \ref{Fig_Jz}(b) and \ref{Fig_Jz}(c), the experimental XLD is well-captured by the $J_z$ fit for both CaFeO\textsubscript{3}/LSAT and CaFeO\textsubscript{3}/SLAO. All major features of the $L_3$ and $L_2$ XLD peaks are replicated. The corresponding x-ray absorption spectrum of the optimized XLD fit for CaFeO\textsubscript{3}/LSAT, shown in Fig. \ref{Fig_Jz}(d), also has excellent agreement with experiment. The coefficients for each $J_z$ spectrum are listed in the Supplemental Material \cite{SI_XLD}. Fig. \ref{Fig_Jz}(e) highlights that the goodness of fit, $\chi^2$, is a strong function of $\Delta$, and the lowest $\chi^2$ values are obtained for $\Delta < 0$, further confirming that CaFeO\textsubscript{3} is a negative charge transfer material. We find that $\Delta = -2.0$ eV provides the best fit to experiment, which is in good agreement with previously reported values for formal Fe\textsuperscript{4+} SrFeO\textsubscript{3} \cite{Bocquet_SFO_ligand_holes}, and is more negative than the rare-earth nickelates \cite{Wu_nickelate_SL_XLD, Robert} but not so negative that the $t_{2g}$ and $e_g$ levels are inverted, as in some compounds \cite{Ushakov_crystal_field}. This value sets the number of self-doped ligand holes, and as seen in Fig. \ref{Fig_Jz}(e), for $\Delta = -2.0$ eV the Fe $e_g$ occupation is 1.85 electrons. This large Fe $e_g$ occupation is consistent with the small measured Fe $3d$ orbital polarization.

The XLD fits also reproduce the measured Fe orbital polarization. Converting the preferential $x^2-y^2$ and the preferential $3z^2-r^2$ fit contributions to a $d_{x^2-y^2}$ and $d_{3z^2-r^2}$ occupation, we find that the orbital occupation for tensile CaFeO\textsubscript{3}/LSAT exhibits a small preference for $d_{3z^2-r^2}$, where $d_{x^2-y^2}$ has 0.91 electrons and $d_{3z^2-r^2}$ has 0.94 electrons \cite{SI_XLD}. This difference of 0.03 electrons agrees well with the difference of 0.05 electrons determined by the sum rule analysis of the XLD integrals. For the compressively strained film, CaFeO\textsubscript{3}/SLAO, the best $J_z$ fit is with equal occupation of 0.93 electrons in both $d_{x^2-y^2}$ and $d_{3z^2-r^2}$.

Of particular note is the sign of the strain-induced crystal field $e_g$ splitting, $\Delta e_g$, that produces the best agreement with experiment. As seen in Fig. \ref{Fig_Jz}(f), for tensile strained CaFeO\textsubscript{3}/LSAT, the lowest $\chi^2$ value occurs for +40 meV; for compressively strained CaFeO\textsubscript{3}/SLAO, -30 meV produces the best fit. These magnitudes are of the same order as other similarly strained perovskite oxides \cite{Aruta_LSMO_XLD, Wu_nickelate_SL_XLD, Gilberto_Nickelate_XLD}. Importantly, the respective signs indicate an $e_g$ splitting consistent with the traditional ligand field model: $\Delta e_g > 0$ implies that $d_{x^2-y^2}$ is lower in energy than $d_{3z^2-r^2}$, which would be expected for tensile strain, and \textit{vice versa} for compressive strain. This provides a critical insight: The energetic landscape of the Fe $3d$ orbitals follows the typical ligand field understanding, where, for example, tensile strain lowers $d_{x^2-y^2}$ in energy relative to $d_{3z^2-r^2}$. Despite this, the $d_{x^2-y^2}$ orbital has fewer electrons than $d_{3z^2-r^2}$ in the film under tensile strain, indicating an inversion in orbital polarization.

What, then, overrides the $\Delta e_g$ splitting and produces the inverted $e_g$ orbital occupation? Although oxygen vacancies can be equatorially or apically ordered under epitaxial strain \cite{Spaldin_O_vacancy_ordering}, the resulting preferential orbital occupation would be opposite of the results here. Moreover, because our experimental findings are not replicated by previous density functional theory calculations \cite{Antonio_CFO_strain}, we propose a new mechanism. It is well known that perovskites can accommodate epitaxial strain by changes in both bond lengths and rotations of the octahedral complexes surrounding the transition metal (TM) cation \cite{Miniotas_strain_rotations, Xie_strain_rotations, Steve_octahedral_rotations_strain, Rondinelli_Spaldin_Adv_Mater}. Rotations alter the TM-O-TM bond angle, and angles less than $180^{\circ}$ have reduced orbital overlap and thus narrower bands. For a perovskite that exhibits rotations in its bulk form, such as CaFeO\textsubscript{3}, biaxial tensile strain increases the in-plane TM-O-TM bond angle towards $180^{\circ}$, whereas the out-of-plane angle decreases further and is typically more strongly affected than the in-plane angles \cite{Rondinelli_Spaldin_Adv_Mater, Antonio_CFO_strain}. Thus in the simplest approximation where strain is accommodated predominantly by octahedral rotations, under tensile strain one would expect the in-plane (x, y) bandwidth to increase and the out-of-plane (z) bandwidth to decrease.

Such anisotropic bandwidth effects can lead to an inverted orbital polarization in compounds with greater-than-half-filled bands. As illustrated in Fig. \ref{Fig_Jz}(g) for the case under tensile strain, $\Delta e_g > 0$ shifts the band center of masses, but a broadening of the $x^2-y^2$ band and a narrowing of the $3z^2-r^2$ band can result in $3z^2-r^2$ being more occupied than $x^2-y^2$. The precise orbital polarization is expected to depend on the specific bandstructure and Fermi level position. For bands with half-filling or less, the same anisotropic bandwidths result in the conventional orbital polarization and thus do not replicate our findings \cite{SI_XLD}. We further note that this effect does not require metallicity and indeed when repeating the Fe $L$-edge XLD measurements at lower temperatures (180 K) where CaFeO\textsubscript{3} is insulating, the inverted orbital polarization is maintained \cite{SI_XLD}.

In summary, we have shown that epitaxially strained films of CaFeO\textsubscript{3} exhibit orbital polarization that responds to the strain state in a way that requires considerations beyond the commonly assumed ligand field model. By analyzing the x-ray linear dichroism with multiplet ligand field simulations, we find that under tensile strain the $e_g$ electronic population is weighted towards $d_{3z^2-r^2}$ orbitals, despite being ${\sim}$40 meV higher in energy than $d_{x^2-y^2}$. The opposite is observed under compressive strain. We propose an explanation for this behavior by considering anisotropic modifications of the bandwidth of the $e_g$ states, in which under tensile strain a broadened $d_{x^2-y^2}$ band and a narrowed $d_{3z^2-r^2}$ lead to this inverted orbital polarization configuration. This scenario is consistent with the orbital energetic landscape as determined by ligand field theory, as well as the measured film strain, under the assumption that strain is accommodated primarily by octahedral bond rotations. More generally, our results demonstrate that effects typically not considered in the conventional understanding of strain-induced orbital polarization can mitigate or even invert the orbital polarization. This highlights that the interpretation of orbital polarization in ultrathin films and short-period superlattices \cite{Chak_asymmetric_XLD, Freeland_LNO_XLD, Keimer_Nickelate_reflectometry, Cao_XLD, Wu_PrNiO3_orbital_polarization, Wu_nickelate_SL_XLD, Disa_Ahn_tricolor_SL}, where non-bulk octahedral rotations can be induced, should include such considerations. Additionally, these results demonstrate that bandwidth control is a potentially new way to engineer orbital polarization in correlated oxides.

\begin{acknowledgments}
We thank G. Sawatzky and A. Fujimori for helpful discussions. PCR and SJM were supported by the Army Research Office, grant number W911NF-15-1-0133, and film synthesis at Drexel utilized deposition instrumentation acquired through an Army Research Office DURIP grant (W911NF-14-1-0493). RJG was supported by the Natural Sciences and Engineering Research Council of Canada. Work at Brookhaven National Laboratory was supported by the U.S. Department of Energy (DOE), Office of Basic Energy Sciences under Contract No. DE-SC0012704 and Early Career Award Program under Award No. 1047478. This work used resources at the Advanced Light Source, which is a DOE Office of Science User Facility under contract No. DE-AC02-05CH11231, and at Beamline 23-ID-1 of the National Synchrotron Light Source II, a DOE Office of Science User Facility operated for the DOE Office of Science by Brookhaven National Laboratory under Contract No. DE-SC0012704.
\end{acknowledgments}

Strongly Correlated Electrons (cond-mat.str-el); Materials Science (cond-mat.mtrl-sci)