A common origin of magnetism from planets to white dwarfs

Abstract: Isolated magnetic white dwarfs have field strengths ranging from kilogauss to gigagauss. However, the origin of the magnetic field has not been hitherto elucidated. Whether these fields are fossil, hence the remnants of original weak magnetic fields amplified during the course of the evolution of their progenitor stars, or are the result of binary interactions or, finally, they are produced by other internal physical mechanisms during the cooling of the white dwarf itself, remains a mystery. At sufficiently low temperatures white dwarfs crystallize. Upon solidification, phase separation of its main constituents, 12C and 16O, and of the impurities left by previous evolution occurs. This process leads to the formation of a Rayleigh-Taylor unstable liquid mantle on top of a solid core. This convective region, as it occurs in Solar System planets like the Earth and Jupiter, can produce a dynamo able to yield magnetic fields of strengths of up to 0.1 MG, thus providing a mechanism that could explain magnetism in single white dwarfs.

\section{Introduction}

White  dwarfs  are the  most  common  end-point of  stellar  evolution
\citep{alth10}.   Some of  them have  measurable magnetic  fields with
strengths  ranging  from  $10^3$   to  $10^9$~G  \citep{ferr15}.   The
incidence  of magnetism  is  a matter  of  debate.  In  volume-limited
surveys -- those that select stars  within a maximum distance from the
Sun  --  the incidence  of  magnetism  is $\sim  20\%$  \citep{kawk07,
giam12,  sion14},   while  for  magnitude-limited  surveys   --  those
selecting  stars brighter  than a  given apparent  magnitude --  it is
$\sim 8\%$  \citep{lieb03, kepl13,  kepl16}.  Nevertheless,  because a
population  of white  dwarfs with  magnetic fields  weaker than  $\sim
1$~kG may exist, the fraction of magnetic white dwarfs could be larger
\citep{kawk12,   koes11,    ferr15}.    Observations    also   suggest
\citep{valy99,  lieb03,  sion14,  kawk14,  valy14,  holl15}  that  the
fraction of white dwarfs with strong  magnetic fields is larger at low
effective temperatures.  This could  indicate that magnetic fields are
amplified during the evolution  of white dwarfs.  This interpretation,
however, has  been questioned \citep{ferr15} because  apparently there
is no clear correlation between the strength of the magnetic field and
the luminosity,  except for low  magnetic field ($B \lesssim  0.1$ MG)
white  dwarfs. Interestingly  enough,  the observed  paucity of  white
dwarfs with magnetic fields between 0.1 and 1~MG \citep{kawk12,koes11}
suggests a bimodal distribution \citep{ferr15}.

The origin  of magnetic white  dwarfs remains  unknown, and up  to now
three scenarios have  been proposed.  Within the first  of them, white
dwarf  magnetic fields  are  the  fossil remnants  of  those of  their
progenitors. Thus, the  progenitors of magnetic white  dwarfs could be
main  sequence Ap/Bp  stars \citep{ange81,  wick05}, or  could be  the
result  of   field  amplification  during  the   helium-burning  phase
\citep{levy74}.  Within the second scenario, magnetic white dwarfs are
the  result of  the evolution  of binary  systems.  In  this case  the
magnetic field  is amplified by  a dynamo that operates  either during
the common envelope phase \citep{tout08,  nord11} or in the hot corona
produced  during the  merger of  two white  dwarfs \citep{garc12}.   A
difficulty  of  these  two   scenarios  is  that  detailed  population
synthesis studies cannot reproduce the number of stars observationally
found \citep{ferr15}.  Finally, within the third scenario the field is
generated in the  outer convective envelope that is  formed during the
evolution of  single white dwarfs.  However,  theoretical calculations
show that  the strength of  the resulting field is  $\lesssim 0.01$~MG
\citep{font73,trem15}, well  below the observed magnetic  fields, so a
more efficient mechanism should be invoked.

The  evolution  of   white  dwarfs  can  be  described   as  a  simple
gravothermal process \citep{alth10}. The core  of the vast majority of
white dwarfs is  made of a completely ionized mixture  of $^{12}$C and
$^{16}$O,  and   some  minor  chemical  species   like  $^{22}$Ne  and
$^{56}$Fe. During  their evolution, two physical  processes modify the
internal  chemical  profiles  of  white  dwarfs.   The  first  one  is
gravitational  settling of  neutron-rich species  in the  liquid phase
\citep{brav92,bild01,garc10,cami16}, which  occurs at  moderately high
luminosities.  The second one is phase separation upon crystallization
\citep{iser97,iser00,garc10},    and    takes     place    at    lower
luminosities.  In both  cases,  the energy  involved  is large,  $\sim
2\times 10^{46}$~erg.

For the  sake of simplicity, here  we will only consider  white dwarfs
made of $^{12}$C and $^{16}$O.  The phase diagram of the carbon-oxygen
mixture is of the spindle  form \citep{segr93,horo10}. Since the solid
phase is oxygen-rich,  hence denser than the liquid  phase, when white
dwarfs crystallize it settles down,  while the carbon-rich liquid left
behind    is    redistributed   by    Rayleigh-Taylor    instabilities
\citep{iser97,iser00}.  This configuration, a solid core surrounded by
a convective  mantle driven  by compositional  buoyancy is  similar to
that found in  the core of the Earth, where  the light element release
associated  with the  inner core  growth is  a primary  driver of  the
dynamo \citep{list95}.  Given this analogy,  in this Letter we analyze
if the  convective mantle  of crystallizing  white dwarfs  can produce
magnetic fields of the observed strengths.

\section{The stellar dynamo}

For  a   typical  white  dwarf   of  mass  $\simeq   0.61\,  M_{\sun}$
solidification starts at $\log(L/L_{\sun}) \simeq -3.5$, whereas for a
massive one  of $1.0 \,  M_{\sun}$ crystallization sets in  at $\simeq
-2.6$. At this point of the evolution the temperatures of their nearly
isothermal cores are, respectively, $\log T_{\rm c}\, ({\rm K}) = 6.3$
and $6.6$,  whereas their central  densities are $\log  \rho_{\rm c}\,
({\rm g\, cm^{-3}})=  6.6$ and $7.5$ \citep{sala10}.  For  the sake of
conciseness, and unless specified otherwise, here we will only discuss
the case of the heavier white dwarf, since it is observationally found
that magnetic white dwarfs are more massive than usual \citep{ferr15}.

\begin{figure}
\plotone{fig1.eps}
\caption{Evolution  of the  inner and  outer radii  of the  convective
  mantle  of  a  carbon-oxygen  white  dwarf  as  a  function  of  the
  luminosity. The  upper and lower  panels correspond to  white dwarfs
  with  masses  $1.0$  and   $0.6\,  M_{\sun}$,  respectively.   Their
  respective total  radii are $R  = 4.7  \times 10^8$ and  $7.5 \times
  10^8$~cm.
\label{f1}}
\end{figure}

As can be  seen in Fig.~\ref{f1}, for a $1.0  \,M_{\sun}$ white dwarf,
when  solidification   starts  the  Rayleigh-Taylor   unstable  region
encompasses a  large fraction  of the  radius of  the star.  This also
holds for  the $0.6 \,M_{\sun}$  model star.  As the  $1.0 \,M_{\sun}$
white  dwarf cools,  the outer  edge of  the convective  mantle barely
moves,  while  the  inner   edge  moves  progressively  outwards.   At
$\log(L/L_{\sun})\simeq -2.9$, the outer edge of the mantle is located
at $r/R  \simeq 0.65$, and  the inner edge  at $\simeq 0.4$.   At this
moment,   the   outer   edge   moves  abruptly   outwards   and   when
$\log(L/L_{\sun})\simeq -3.1$ it almost reaches the surface, while the
inner  edge continues  its progression  at  a slower  pace.  The  main
difference with the $0.6 \,M_{\sun}$ white  dwarf is that in this case
the  outer radius  stabilizes at  $r/R\simeq 0.9$.   Finally, in  both
cases,  the convective  region disappears  at about  $\log(L/L_{\sun})
\simeq -4.5$.

The  density contrast  between the  carbon enriched  material and  the
ambient liquid mixture is $\delta \rho /\rho \sim 10^{-3}$, leading to
effective  accelerations $a_{\rm  eff} =  g(\delta \rho/\rho)\simeq  2
\times  10^6$~cm~s$^{-2}$. However,  since  the  viscosity of  Coulomb
plasmas  is very  small, drag  cannot be  neglected and  results in  a
limiting    velocity     of    the    turbulent     eddies,    $v_{\rm
s}=\sqrt{(3/8)C_{\rm b} a_{\rm eff} D_{\rm  b}}$, where $D_{\rm b}$ is
the radius of the bubble, and $C_{\rm b} \approx 0.2$ in the spherical
case \citep{garc95}.  The  value of $D_{\rm b}$ is not  yet known, but
it  cannot be  larger than  the curvature  radius at  the edge  of the
crystallized core.  When the same value found for the Earth is adopted
$D_{\rm b} \sim 0.1 R_{\rm core}$ \citep{moff94}, we obtain velocities
of $\sim 35$~km~s$^{-1}$.

For  the thermodynamical  conditions found  in white  dwarf interiors,
viscosities   are   very   small  and   conductivities   rather   high
\citep{nand84}.   In  particular, for  the  heavier  white dwarf,  the
electrical conductivity is $\sigma  = 1.3\times 10^{21}$~s$^{-1}$, the
magnetic diffusivity is  $\eta=5.6 \times 10^{-2}$~cm$^2$~s$^{-1}$ and
the kinematic viscosity is $\nu =3.13 \times 10^{-2}$~cm$^2$~s$^{-1}$.
Adopting $\sim  2\times 10^8$~cm  for the  characteristic size  of the
convective mantle  (see Fig.~\ref{f1}, upper panel),  the Reynolds and
magnetic Reynolds numbers\footnote{The  magnetic Reynolds number, Rm$=
ul/\eta$, where  $u$ and  $l$ are typical  velocities and  lengths and
$\eta$ is  the magnetic diffusivity, measures  the relative importance
of induction versus dissipation, while  the magnetic Prandtl number is
the  ratio between  the magnetic  and the  ordinary Reynolds  numbers,
Pm=Rm/Re.}  range from  $\sim 10^{14}$ to $10^{15}$,  and the magnetic
Prandtl  number is  $\approx 0.58$.   Furthermore, the  characteristic
Ohmic decay time $\tau_\Omega  \sim R^2_{\rm WD}/\eta\simeq 10^{18}$~s
is much  longer than the age  of the white dwarf.   Consequently, once
the magnetic  field is generated,  it will remain almost  constant for
long times \citep{cumm02}.

\begin{figure}
\plotone{fig2.eps}
\caption{Convective  energy  density available  for  the  dynamo as  a
  function of the  luminosity, for white dwarfs of  masses 1.00, 0.87,
  0.77, 0.68, 0.61,  and $0.55\, M_{\sun}$. The  top panel corresponds
  to  white dwarfs  with H-deficient  atmospheres (non-DA),  while the
  bottom one corresponds to stars with H-rich atmospheres (DA).
  \label{f2}}
\end{figure}

\begin{table*}[t]
\centering
\caption{Properties of magnetic white dwarfs within 20~pc of the Sun.}
\label{tb1}
\begin{threeparttable}
\begin{small}
\begin{tabular}{lccccccccc}
\hline
\hline
Name          &  Type  & $B$~(MG)           & Period~(h)           & $M\,(M_{\sun})$ & $\log(L/L_{\sun})$ & $M_{\rm s}/M_{\rm WD}$  & $\rho_{\rm b}$   & $r_{o}-r_{i}$      & $q_{\rm top}$       \\
\hline
WD 0009+501   &     DA &         0.2        &  2 -- 20             & $0.73 \pm 0.04$ & $-3.72$            & 0.43                    & $3.0\times 10^6$ & $2.3\times 10^8$   & $7.1\times 10^{10}$ \\
WD 0011$-$134 &     DA & $16.7\pm 0.6$      &    ---               & $0.72 \pm 0.07$ & $-3.85$            & 0.52                    & $2.3\times 10^6$ & $2.6\times 10^8$   & $5.7\times 10^{10}$ \\
WD 0322$-$019 &     DA & 0.120              &    ---               & $0.63 \pm 0.05$ & $-4.02$            & 0.52                    & $1.6\times 10^6$ & $5.2\times 10^8$   & $3.0\times 10^{10}$ \\
WD 0413$-$077 &     DA & $0.0023\pm 0.0007$ &    ---               & $0.59 \pm 0.03$ & $-1.85$            & 0.00                    & ---              & ---                & ---                 \\
WD 0548$-$001 & non-DA & $\sim 10$          &   4.117              & $0.69 \pm 0.03$ & $-3.80$            & 0.82                    & $7.0\times 10^5$ & $1.2\times 10^8$   & $6.5\times 10^{10}$ \\
WD 0553+053   &     DA & $20\pm 3$          &   0.97               & $0.72 \pm 0.03$ & $-3.91$            & 0.35                    & $1.9\times 10^6$ & $2.5\times 10^8$   & $4.6\times 10^{10}$ \\
WD 0728+642   &     DA & $0.0396\pm 0.0116$ &    ---               & $0.58 \pm 0.00$ & $-4.00$            & 0.36                    & $1.5\times 10^6$ & $2.0\times 10^8$   & $3.7\times 10^{10}$ \\
WD 0912+536   & non-DA & $\sim 100$         &  31.9224             & $0.75 \pm 0.02$ & $-3.57$            & 0.71                    & $1.5\times 10^6$ & $1.6\times 10^8$   & $1.1\times 10^{11}$ \\
WD 1008+290   & non-DA & $\sim 100$         &    ---               & $0.68 \pm 0.01$ & $-4.31$            & 0.93                    & $9.8\times 10^4$ & $3.0\times 10^6$   & $2.1\times 10^{10}$ \\
WD 1036$-$204 & non-DA & 50                 & $636 \pm 36$         & $0.60 \pm 0.01$ & $-4.19$            & 0.76                    & $1.9\times 10^5$ & $1.1\times 10^7$   & $2.5\times 10^{10}$ \\
WD 1309+853   &     DA & $4.9\pm0.5$        &    ---               & $0.71 \pm 0.02$ & $-4.01$            & 0.60                    & $1.6\times 10^6$ & $2.3\times 10^8$   & $3.4\times 10^{10}$ \\
WD 1748+708   & non-DA & $\ga 100$          & $\ga 1.7\times 10^4$ & $0.79 \pm 0.01$ & $-4.07$            & 0.79                    & $2.0\times 10^5$ & $4.5\times 10^6$   & $4.6\times 10^{10}$ \\
WD 1829+547   & non-DA & $170-180$          &    ---               & $0.90 \pm 0.07$ & $-3.94$            & 0.85                    & $2.7\times 10^5$ & $7.2\times 10^6$   & $8.3\times 10^{10}$ \\
WD 1900+705   & non-DA & $320\pm20$  & $\ga 8.5 \times 10^5$       & $0.93 \pm 0.02$ & $-2.88$            & 0.18                    & $1.4\times 10^7$ & $2.4\times 10^8$   & $1.2\times 10^{12}$ \\
WD 1953$-$011 &     DA & $0.1-0.5$          & 34.60224             & $0.73 \pm 0.03$ & $-3.38$            & 0.13                    & $5.6\times 10^6$ & $2.4\times 10^8$   & $3.3\times 10^{11}$ \\
WD 2105$-$820 &     DA & $0.010\pm 0.001$   &    ---               & $0.74 \pm 0.13$ & $-2.93$            & 0.00                    & ---              & ---                & ---                 \\
WD 2359$-$434 &     DA & $0.0031\pm0.0005$  &    ---               & $0.78 \pm 0.03$ & $-3.26$            & 0.12                    & $7.2\times 10^6$ & $1.8\times 10^8$   & $7.3\times 10^{11}$ \\
\hline
\end{tabular}
\end{small}
\begin{tablenotes}[para,flushleft]
Data obtained from \citet{kawk07}, \citet{giam12} and references therein.
\end{tablenotes}
\end{threeparttable}
\end{table*}

The aim  of dynamo  scaling theories  is to  explain the  intensity of
magnetic fields in  terms of the properties of the  region hosting the
dynamo.  The most  comprehensive of these theories  takes into account
the  balance  between  the  ohmic  dissipation  and  the  energy  flux
available to  the dynamo  \citep{chri10}. The  magnetic fields  of the
Earth and Jupiter  are generated by convective dynamos  powered by the
cooling and  chemical segregation of  their interiors and it  has been
shown that the scaling law based  on the energy fluxes can be extended
to fully convective  stars like T~Tauri and  rapidly rotating M~dwarfs
\citep{chri09}.  Therefore,  it is  natural to  use these  theories to
compute  the  magnetic field  produced  by  the convective  mantle  of
crystallizing white dwarfs.

If  the dynamo  is saturated,  the magnetic  Reynolds number  is large
enough, and  convection is described  by the mixing  length formalism,
this scaling law can be written as \citep{chri10}:
\begin{equation}
\frac{{B^2 }}{{2\mu _0 }} = cf_\Omega
\frac{1}{V}\int\limits_{r_i }^{r_o }
{\left[ {\frac{{q_c \left( r \right)\lambda\left( r \right)}}
{{H\left( r \right)}}} \right]}^{{2\mathord{\left/{\vphantom {2 3}} \right. \kern- \nulldelimiterspace} 3}}
\,\rho \left( r \right)^{{1 \mathord{\left/
 {\vphantom {1 3}} \right.
 \kern-\nulldelimiterspace} 3}} 4\pi r^2 dr
\label{chris_rel}
\end{equation}
where the integral is the energy of the convective mantle, $E$, $c$ is
an adjustable parameter, $\mu_0$ is the vacuum permeability, $f_\Omega
\le 1$ is the ratio of the Ohmic dissipation to the total dissipation,
$q_{\rm c}$  is the  convected energy  flux, $H$  is the  scale height
(temperature and compositional),  $V$ is the volume  of the convective
region, $r_o$ and  $r_i$ are its outer and inner  radii, and $\lambda$
is the mixing length (the minimum between the density scale height and
the  size  of the  convective  zone,  $r_o  -  r_i$).  Here  we  adopt
$f_\Omega  = 1$  and  $\lambda=H$.   Also, we  use  the BasTI  cooling
sequences \citep{sala10} which,  in addition to the  release of latent
heat,  consider  the  energy   released  by  chemical  differentiation
\citep{moch83,segr94}.

\section{Results and discussion}

Figure~\ref{f2} displays the density of convective energy available to
the dynamo, namely the  r.h.s.  of Eq.~(\ref{chris_rel}) excluding the
factor $cf_\Omega$, for several white dwarf masses.  The dynamo starts
when  solidification sets  in, then  reaches a  maximum and,  finally,
slowly declines as the bottom  of the convective mantle moves outwards
at fainter luminosities.   This is the result of  two effects.  First,
the mass  of the  star, that  determines when  crystallization starts.
White  dwarfs with  larger masses  crystallize at  larger luminosities
because their  central densities are  also larger.  The  luminosity of
the star  also plays  a role,  since it determines  the rate  at which
crystallization takes place.  This luminosity, in turn, depends on the
transparency of the  atmosphere and the temperature of  the core.  The
maximum energy available  to the dynamo ranges from $\log  E = 8.2$ to
9.6,  and from  8.3 to  9.7 (c.g.s.   units) for  DA and  non-DA white
dwarfs,   respectively,    when   masses   ranging   from    0.54   to
$1.00\  M_{\sun}$  are  considered.   The  scaling  law  relating  the
magnetic fields of  the Earth, Jupiter, T~Tauri and  M~dwarf stars can
be approximated by $\log B = 0.5  \log E - 5.42$ -- see Fig.~\ref{f3},
solid  line  -- and  predicts  that  the  maximum  field that  can  be
generated at the top of the dynamo range from $\sim 0.05$ to 0.25~MG.

\begin{figure}
\plotone{fig3.eps}
\caption{Magnetic field intensity  as a function of  the dynamo energy
  density. Earth and Jupiter are  represented using black symbols with
  their corresponding error bars. T~Tauri and M~dwarfs are shown using
  cyan and  magenta symbols, respectively \citep{chri09}.   The DA and
  non-DA white dwarfs listed in Table~\ref{tb1} are plotted as red and
  blue  symbols, respectively.   The  solid line  is the  relationship
  relating  the magnetic  fields of  the Earth,  Jupiter, T~Tauri  and
  M~dwarf stars, while dotted lines  allow for an additional deviation
  of a factor of  3 from it. The dashed lines  help to represent where
  non-DA  stars  cluster.  The  top  and  the bottom  panels  display,
  respectively, the intensity  of the magnetic field as  a function of
  the present and of the maximum energy densities of the dynamo.
\label{f3}}
\end{figure}

The X-ray  emission of stars  with convective regions  correlates with
the rotation period,  and with magnetic activity.  The  existence of a
plateau  in  the  relationship  linking magnetic  activity  and  X-ray
luminosity  suggests  that dynamos  saturate  when  the Rossby  number
Ro$=P_{\rm rot}/\tau<0.1$,  being $\tau$ the convective  turnover time
\citep{wrig11,wrig16}.   For a  star of  mass $1.0  \, M_{\sun}$,  and
adopting the  typical values  found in Sect.~2,  the turnover  time is
$\sim 90$~s.   Since the critical  stability rotation period  is $\sim
6$~s,  we  conclude  that  rapidly  rotating  white  dwarfs  can  have
saturated  dynamos.  Because  white  dwarf spectral  lines are  broad,
measuring the rotation periods is a difficult task, and have only been
measured  in  pulsating  and   magnetic  white  dwarfs.   Thus,  their
distribution   of  angular   velocities   is   not  well   determined.
Observations  indicate  that  white  dwarfs  are  slow  rotators,  but
strictly  speaking, the  existence  of rapidly  rotating white  dwarfs
cannot be discarded \citep{kawa15}.   Within our scenario white dwarfs
with  high magnetic  fields would  rotate rapidly  and have  saturated
dynamos, while  stars with  weak fields would  rotate slowly  and have
non-saturated dynamos.  Therefore, our  mechanism could easily explain
white dwarfs with magnetic fields $B \lesssim 0.1$~MG.

Table~\ref{tb1}   shows   the   properties  of   white   dwarfs   with
carbon-oxygen cores within 20~pc  of the Sun \citep{kawk07,giam12}. We
list the name of the star,  the spectral type, the magnetic field, the
rotation period in hours, the mass  and luminosity in solar units, the
relative  mass of  the crystallized  core  and the  properties of  the
convective mantle  (the density at the  inner edge, the size,  and the
flux at  the top, all  in c.g.s. units).   It is important  to realize
that all  of them,  except WD~0413$-$077 (40~Eri~B)  and WD~2105$-$820
(L24-52) have crystallized cores.  For each of these stars we computed
the   energy   density   available   to  the   dynamo   according   to
Eq.~(\ref{chris_rel})  employing the  value of  $c$ derived  using the
scaling law for the Earth and Jupiter, T~Tauri stars and M~dwarfs.  In
the top  panel of  Fig.~\ref{f3} we  show the  set of  calculations in
which the energy density of the  dynamo has been computed adopting the
actual size  of the  convective mantle  corresponding to  the observed
luminosity  of the  star,  whereas in  the bottom  panel  we show  the
magnetic field strength when the  maximum energy density of the dynamo
shown in  Fig.~\ref{f2} is  adopted. This  is a  reasonable assumption
given the very  long ohmic decay timescale.  In both  cases, to relate
the  outer and  inner  magnetic  fields we  used  the prescription  of
\citet{chri09}.   The  distribution  of magnetic  field  strengths  is
clearly bimodal. Most  DA white dwarfs cluster around  the scaling law
(the  solid line  in Fig.~\ref{f3}),  while non-DA  white dwarfs  have
magnetic fields  larger than DA  stars with similar energy  density of
the dynamos.  This points towards to a different origin of the dynamo,
and/or to  other scaling  laws like  the Coriolis-Inertial-Archimedian
balance or  the Magnetic-Archimedian-Coriolis  balance \citep{chri10},
because the  mixing length scaling  law considered here  demands $\sim
10^{15}$~erg~s$^{-1}$ to produce fields of $\sim100$~MG.

It is  important to realize that  the time necessary for  the magnetic
field to diffuse through the radiative  outer layers is very long, and
only the  convective mantle of  massive white  dwarfs is close  to the
surface, as  clearly shown in  Fig~\ref{f1}. The final  magnetic field
strength will depend on how  this mantle interacts with the convective
envelope of  cool white  dwarfs. Furthermore,  only white  dwarfs with
progenitors  more  massive than  $2.3\,  M_{\sun}$  can have  rotation
periods of minutes \citep{kawa15}.  This,  together with the fact that
the  energy of  the dynamo  induced by  crystallization is  larger for
massive white dwarfs could naturally explain why magnetic white dwarfs
are more massive  than average.  This scenario could  also explain why
magnetic white dwarfs are preferentially cool.  Moreover, it is likely
that low-mass white dwarfs would  have weaker magnetic fields confined
in  their   interiors.   Therefore,  this  scenario   puts  forward  a
tantalizing possibility, namely, that  the magnetic fields observed in
planets, non-evolved  stars and white  dwarfs share a  common physical
origin.  In summary, our calculations indicate that the magnetic field
observed in a sizable fraction of all white dwarfs could be the result
of a dynamo generated by phase separation upon crystallization.

Furthermore, our  scenario does not preclude  other possibilities, but
instead alleviates one of the main drawbacks of the current hypotheses
to explain  magnetic white dwarfs,  since both the fossil  field model
and  the binary  scenario predict  an insuficient  number of  magnetic
white dwarfs.  For instance, it has  been long suspected that at least
a fraction  of high-field magnetic  white dwarfs could  originate from
the merger of  two white dwarfs \citep{wick00}.  During  the merger, a
hot,   differentially-rotating,   convective    corona   forms.    The
temperatures reached during the coalescence  are so high that hydrogen
is burned during the early phases of the merger.  This corona is prone
to  magnetorotational instability,  and  it has  been  shown that  can
produce  magnetic  fields with  the  energy  required to  explain  the
magnetic fields of non-DA stars shown in Fig.~\ref{f3} \citep{garc12}.
This  could explain  why many  white  dwarfs with  very high  magnetic
fields  are  H-deficient.  In  this  case  a  dynamo of  a  completely
different nature would be operating.

\acknowledgements
This work has been supported by MINECO grants ESP2013- 47637-P, ESP2015-66134-R (J.I.), and AYA2014-59084-P (E.G.-B.), by the European Union FEDER funds, by grants 2014SGR1458 (J.I.), 2014SGR0038 (E.G.-B.) of the AGAUR, and by the CERCS program of the Generalitat de Catalunya.

Solar and Stellar Astrophysics (astro-ph.SR); Earth and Planetary Astrophysics (astro-ph.EP)


Magnetism and its microscopic origin in iron-based high-temperature superconductors

Abstract: High-temperature superconductivity in the iron-based materials emerges from, or sometimes coexists with, their metallic or insulating parent compound states. This is surprising since these undoped states display dramatically different antiferromagnetic (AF) spin arrangements and N$\rm \acute{e}$el temperatures. Although there is general consensus that magnetic interactions are important for superconductivity, much is still unknown concerning the microscopic origin of the magnetic states. In this review, progress in this area is summarized, focusing on recent experimental and theoretical results and discussing their microscopic implications. It is concluded that the parent compounds are in a state that is more complex than implied by a simple Fermi surface nesting scenario, and a dual description including both itinerant and localized degrees of freedom is needed to properly describe these fascinating materials.

\section{Introduction}

Soon after the discovery of high critical
temperature (high-$T_c$) superconductivity in copper oxides  \cite{bednorz},
neutron scattering studies revealed that
the parent compounds of these superconductors have an antiferromagnetic
(AF) ground state with a simple collinear spin structure (Fig. 1a) \cite{vaknin,tranquada88}.
Because the associated AF spin fluctuations
may be responsible for electron pairing and superconductivity \cite{scalapino,dagotto,palee},
over the past 25 years a tremendous effort has focused on
characterizing the interplay between magnetism and superconductivity in these materials \cite{fujita}.
In the undoped state, the parent compounds of copper oxide superconductors
are Mott insulators and have exactly one valence fermion with spin 1/2 for each copper atom,
leading to robust electronic correlations and
localized magnetic moments \cite{dagotto,palee}.  Superconductivity emerges after introducing charge
carriers that suppress the static AF order. Although the strong Coulomb repulsion in the parent
compounds is screened by the doped charge carriers,
the electronic correlations are certainly important for the physics of the doped cuprates,
particularly in the underdoped regime \cite{palee}.

Consider now the iron-based superconductors \cite{johnston,stewart,paglione10}. Several
parent compounds of these materials, such as LaFeAsO, BaFe$_2$As$_2$, NaFeAs, and FeTe, are not
insulators but semimetals \cite{kamihara,rotter,cwchu,mawkuen2}. In these cases,
electronic band structure calculations have revealed that their Fermi Surfaces (FS) are
composed of nearly cylindrical
hole and electron pockets at the $\Gamma(0,0)$ and $M(1,0)/M(0,1)$ points,
respectively \cite{mazin2011n,hirschfeld}.
The high density of states resulting from the
extended momentum space with nearly parallel FS
between the hole and electron pockets leads to an enhancement of the particle-hole susceptibility. This  suggests that
FS nesting among those pockets
could induce spin-density-wave (SDW) order at the in-plane AF wavevector ${\bf Q}_{AF} = (1,0)$
with a collinear spin structure (Fig. 1b) \cite{dong}, much like the FS-nesting induced SDW in
  pure chromium \cite{fawcett}.
Neutron scattering experiments on
LaFeAsO \cite{cruz}, BaFe$_2$As$_2$ \cite{qhuang},
and NaFeAs \cite{slli09} have reported results compatible with the theoretically predicted
 AF spin structure, albeit with an ordered magnetic moment smaller than expected from
first-principles calculations \cite{mazin08}.
In addition, quasiparticle excitations between the hole and electron FS
can induce $s^{\pm}$-wave superconductivity \cite{mazin2011n,hirschfeld,kuroki08,chubukov,fwang09}.
One of the consequences of this superconducting state is that the imaginary
part of the dynamic susceptibility, $\chi^{\prime\prime}(Q,\omega)$ should have a sharp peak,
termed spin resonance in copper oxide superconductors \cite{eschrig},
at ${\bf Q}_{AF} = (1,0)$ below $T_c$ \cite{maier,korshunov}.  This prediction is also confirmed by inelastic
neutron scattering (INS) experiments in iron-based superconductors such as hole-doped Ba$_{1-x}$K$_x$Fe$_2$As$_2$ \cite{christianson,chenglinzhang,castellan}, electron-doped BaFe$_{2-x}T_x$As$_2$ ($T=$Co, Ni) \cite{lumsden,schi09,dsinosov09,jtpark,clester,hfli,hqluo12},
and FeTe$_{1-x}$Se$_x$ \cite{hamook,qiu09,lumsden2}.
Finally, angle resolved photoemission spectroscopy (ARPES) experiments find
that the general characterization of the FS and
the superconducting order parameter are consistent with the band structure
calculations and with isotropic $s$-wave superconducting gaps \cite{prichard}.
Therefore, at first sight it may appear that antiferromagnetism in the iron-based materials
originates from FS nesting of itinerant electrons, superconductivity must have
a $s^{\pm}$-wave symmetry for related reasons, and electron correlations or local moments
do not play an important role for magnetism and superconductivity \cite{mazin2011n}.

However, although the parent compounds of iron pnictide superconductors
have metallic ground states consistent with
band structure calculations, there are reasons to believe
that electron correlations could be sufficiently strong
to produce an ``incipient''  Mott physics \cite{si2008,qmsi09},
where local moments are as important as itinerant electrons for magnetic, transport,
and superconducting properties in these materials \cite{cfang,cxu}.
In fact, the $s^{\pm}$ pairing symmetry is also naturally derived  in
multi-orbital $t-J$-type models~\cite{seo2008,Fang2011} and recent diagonalization calculations~\cite{nicholson11}
have shown that the AF state, as well as the $A_{1g}$ $s$-wave pairing state,
evolve smoothly from weak to strong coupling, suggesting
that the physics of the pnictides could also be
rationalized based on short length scale concepts not rooted in weak-coupling nesting.
After all, in the context of the copper oxide superconductors, weak coupling studies of
the one-orbital Hubbard model also led to the correct checkerboard AF state and $d$-wave pairing,
showing that these problems can be attacked from a variety of view points.
In addition,  the newly discovered  $A_y$Fe$_{2-x}$Se$_2$ ($A=$ K, Rb, Cs, Tl)
iron-chalcogenide superconductors~\cite{jgguo,mhfang} do {\it not}
exhibit hole pockets~\cite{Wang_122Se, Zhang_122Se,Mou_122Se},
but have strong AF ordered insulating phases with extremely
high N\'{e}el transition temperatures~\cite{wbao1,fye}.
Such a strong magnetism and high superconducting transition temperature ($T_c\approx 33$ K) cannot be explained
by FS nesting since this is based on the enhancement of the particle-hole susceptibility due to an extended
momentum space with nearly parallel Fermi surfaces, i.e. it applies only to particle and hole FS's and not
to purely electronic Fermi pockets.

Since iron-based superconductors have
six electrons occupying the nearly degenerate $3d$ Fe orbitals, the system is intrinsically
multi-orbital and therefore it is technically difficult to
define and study a simple microscopic Hamiltonian to describe the electronic
properties of these materials and characterize
the strength of the electronic correlations. From optical conductivity
measurements \cite{qazilbash09}, it has been argued that electronic
correlations in Fe pnictides are weaker than in underdoped copper-oxides
but are stronger than those of Fermi liquid metals, contrary to the
conclusion based on local density approximation calculations \cite{mazin2011n}.  Therefore,
it is important to determine whether magnetism in Fe-based materials arises from weakly correlated
itinerant  electrons~\cite{mazin2011n}, as in the case of the SDW in chromium \cite{fawcett}, or whether it requires some degree of electron
  correlations~\cite{HG}, or if magnetism is dominated by the contributions
of quasi-localized moments induced by incoherent electronic excitations~\cite{qmsi09} such as in the AF insulating state
of Cu oxides \cite{palee}.

In this review, recent experimental and theoretical progress in the study of
iron-based superconductors is summarized, with focus on the undoped parent compounds.
In section II, the magnetically ordered states in nonsuperconducting
iron pnictides, iron chalcogenides, and iron selenides are discussed.  Section III describes the effect of
electron and hole doping on static AF order and their associated spin excitations.  In section IV, we provide several examples
where deviations from the simple SDW FS nesting picture are prominent.
Finally in section V, we present our perspective on
the importance of electron correlations in these materials.

\section{Magnetic order arrangements and spin waves in the parent compounds}

Although the overall crystal structures and chemical formulas of the copper-oxide superconductors can be quite different,
their parent compounds are all AF
Mott insulators characterized by
the Cu spin structure shown in Fig. 1a, where the tetragonal or pseudo-tetragonal unit cells have
a nearest-neighbor Cu-Cu spacing with $a\approx b\approx 3.8$ \AA.  In the notation
of reciprocal lattice units (rlu) $(2\pi/a,2\pi/b,2\pi/c)$,
the AF Bragg peaks occur at the in-plane ordering wave vectors ${\bf Q}_{AF}=(\pm 1/2+m,\pm 1/2+n)$, where
$m,n=0,\pm 1, \pm 2, \cdots$ rlu, shown as red circles in Fig.~1e \cite{vaknin,tranquada88}.
Time-of-flight INS experiments \cite{coldea,headings} have mapped out spin waves of
the insulating La$_2$CuO$_4$ throughout the Brillouin zone and found no evidence for spin-wave broadening at high energies.
The dispersions of spin waves are well described by a Heisenberg Hamiltonian
with nearest-neighbor (NN) exchange coupling
$J_1=111.8\pm 4$ meV and next-nearest-neighbor (NNN) exchange $J_2=-11.4\pm 3$ meV \cite{coldea}.
Therefore, the dominant magnetic exchange coupling
in La$_2$CuO$_4$ is the NN magnetic interaction and the higher-order interactions amount
to only  $\sim$10\% of the total magnetic energy with a bandwidth of $\sim$320 meV (Fig. 1i).

\begin{figure}[t]
\includegraphics[scale=.7]{Fig1R2}
\caption{Antiferromagnetic structure and spin-wave dispersions for the insulating copper oxide
La$_2$CuO$_4$ and the parent compounds of iron-based superconductors BaFe$_2$As$_2$, FeTe,
and $A_y$Fe$_{1.6+x}$Se$_2$.  The chemical unit cells are marked light green.  The dark and light brown As/Te/Se atoms indicate their vertical positions above and below the Fe-layer, respectively.
(a) The AF structure of La$_2$CuO$_4$, where the chemical unit cell is marked light green.
(b) The collinear AF structure of nonsuperconducting iron pnictides in the FeAs-layer, where spins are aligned anti-parallel along the orthorhombic $a_o$-axis \cite{cruz,qhuang,slli09}. (c) The bi-collinear AF structure of FeTe \cite{webao08,slli08}. (d) The block AF order of the insulating $A_y$Fe$_{1.6+x}$Se$_2$, where
the $\sqrt{5}\times\sqrt{5}$ superlattice structure is marked by solid line with lattice parameter $a_s=8.663$ \AA\ and the orthorhombic lattice cell is shaded green \cite{wbao1,fye}.  The iron vacancies are marked as yellow squares. (e) The wave vector dependence of the AF order
in the $(H,K)$ plane of the reciprocal space for La$_2$CuO$_4$ \cite{coldea};
(f) BaFe$_2$As$_2$ \cite{qhuang};
(g) Fe$_{1.05}$Te \cite{webao08,slli08}; and
(h)  the insulating $A_y$Fe$_{2-x}$Se$_2$ \cite{wbao1,fye}.
(i) Spin-wave dispersions along two high symmetry directions for La$_2$CuO$_4$ \cite{coldea}.  The overall energy
scale of spin waves for copper oxides is about $320$ meV and spin waves are instrumental
resolution limited. (j) Spin-wave dispersions for BaFe$_2$As$_2$ and they broad considerably for energies
above $\sim$100 meV \cite{lharriger}.  (k) Spin-wave dispersions for Fe$_{1.05}$Te, and spin waves are very broad for energies above 30 meV \cite{lipscombe}.
(l) Spin waves for the insulating Rb$_{0.89}$Fe$_{1.58}$Se$_2$ \cite{mywang11}.  In spite of dramatically different dispersions for
various iron-based materials, their spin wave overall energy scales are similar and about 220 meV, less than that
of the insulating copper oxides. Twinning is considered.}
\end{figure}

Four years after the initial discovery of superconductivity in LaFeAsO$_{1-x}$F$_x$ \cite{kamihara},
there are now three major families of
iron-based superconductors: the iron pnictides \cite{johnston,stewart},
iron chalcogenides \cite{mawkuen2,mhfang08},
and alkaline iron selenides \cite{jgguo,mhfang}.  The parent compounds of the pnictides,
such as $A$FeAsO ($A=$ La, Ce, Sm, Pr, etc.),
$A$Fe$_2$As$_2$ ($A=$ Ba, Sr, Ca), and NaFeAs, all have the same collinear AF structure as shown in Fig.~1b,
with a small ordered moment ($<1\ \mu_B$/Fe) and N$\rm \acute{e}$el temperature $T_N\le 200$ K
\cite{cruz,qhuang,cwchu}. The AF spin moments  are aligned along the
weak orthorhombic unit cell $a$-axis direction ($a\approx5.62$, $b\approx 5.57$ \AA).
In reciprocal space, the AF Bragg peaks occur at in-plane ordering
wave vectors ${\bf Q}_{AF}=(\pm 1+m, n)$ and at
${\bf Q}_{AF}\approx (m, \pm 1+n)$ due to twinning (red circles in Fig. 1f),
consistent with the
$\Gamma(0,0)\leftrightarrow M(1,0)/M(0,1)$
FS nesting picture \cite{mazin2011n}.
However, although the calculated FS of the chalcogenides Fe$_{1+y}$Te$_{1-x}$Se$_x$
is similar to that of iron pnictides \cite{subedi}, surprisingly its parent compound
Fe$_{1+y}$Te actually has a bi-collinear spin structure
(Fig. 1c) \cite{webao08,slli08}.  Here, the AF Bragg peaks occur
at ${\bf Q}_{AF}=(\pm 1/2+m,\pm 1/2+n)$ (Fig. 1g) in
the pseudo-tetragonal notation ($a\approx b\approx 5.41$ \AA),  suggesting that
FS nesting cannot induce such AF order.  Finally, the parent compounds of the
alkaline iron selenide $A$Fe$_{1.6+x}$Se$_2$ superconductors
are insulators \cite{jgguo,mhfang}
and form a $\sqrt{5}\times\sqrt{5}$ block AF structure as shown in Fig. 1d
with a large ordered moment ($\sim$$3\ \mu_B$/Fe) along the $c$-axis and
$T_N\approx 500$ K \cite{wbao1,fye}.
 In reciprocal space, defined using the pseudo-tetragonal unit cell
of iron pnictides ($a\approx b\approx 5.41$ \AA),
the block AF Bragg peaks appear at ${\bf Q}_{AF}=(\pm 0.2+m,\pm 0.6+n)$ and
$(\pm 0.6+m,\pm 0.2+n)$ combining left and right chiralities (red circles in Fig. 1h).

Since the parent compounds of iron-based superconductors can have
different AF spin structures and either
metallic or insulating ground states \cite{johnston,stewart,jgguo,mhfang},
the microscopic origin of the
AF order cannot be induced by a simple FS nesting.
If magnetism is relevant for high-$T_c$ superconductivity,
then it would be important to determine magnetic exchange couplings for different
 classes of Fe-based superconductors and compare the results with those
of the copper-oxides \cite{coldea}. For pnictides, INS
experiments have mapped out spin waves
on single crystals of CaFe$_2$As$_2$ \cite{diallo09,jzhao}, SrFe$_2$As$_2$ \cite{raewings},
and BaFe$_2$As$_2$ \cite{lharriger} throughout the Brillouin zone.  Although there are still
debates concerning whether spin waves in
these materials can be described by a pure itinerant picture \cite{diallo09,raewings} or require
local moments \cite{jzhao,lharriger}, the overall spin-wave energy scales
are around 220 meV.
Therefore, magnetic exchange couplings in iron pnictides are clearly smaller than those of copper
oxides (Figs. 1i and 1j). Although spin waves are broadened at high energies, the
spin-wave dispersion curves (Fig. 1j)  can still be described by a Heisenberg Hamiltonian with
strong anisotropic NN exchange couplings ($J_{1a}\gg J_{1b}$) and fairly large NNN exchange
coupling ($J_2$) \cite{jzhao,lharriger}. This large in-plane magnetic exchange coupling
anisotropy has been interpreted as due to possible electronic nematic
phase and/or orbital ordering \cite{jzhao,lharriger}. Table I compares the effective magnetic exchange
couplings of the Fe-based systems studied thus far against those of the insulating copper-oxide La$_2$CuO$_4$.

For the chalcogenides Fe$_{1+y}$Te, the commensurate bi-collinear AF spin structure in Fig. 1c
becomes incommensurate for concentration $y>0.12$ \cite{rodriguez}.
The overall spin-waves energy scale (Fig. 1k) is similar to those of the iron pnictides.
Although the large static ordered moment of $\sim$$2\ \mu_B$/Fe in Fe$_{1+y}$Te \cite{webao08,slli08}
suggests that local moments may be important, spin waves are rather broad in energy
and difficult to fit using a Heisenberg Hamiltonian with only NN and NNN exchange couplings \cite{lipscombe}.
By including third-neighbor (NNNN) exchange couplings, a Heisenberg Hamiltonian can fit the
spin-wave dispersion with an anisotropic ferromagnetic NN exchange
couplings and strong AF NNN exchange coupling (Table I). In a separate INS
 experiment on  Fe$_{1.1}$Te, the total integrated Fe magnetic moment was
 found to increase with increasing temperature from 10 K to 80 K \cite{zaliznyak}.
These results suggest that in the temperature range relevant for superconductivity,
there is a remarkable redistribution of the magnetism arising from both
itinerant and localized electrons.

\begin{figure}[t]
\includegraphics[scale=.55]{Fig2R2}
\caption{The electronic phase diagrams and the evolution of
FS's, static AF order, and spin excitations upon electron or hole doping to BaFe$_2$As$_2$.
(a) The AF and superconducting phase diagram for hole-doped Ba$_{1-x}$K$_x$Fe$_2$As$_2$.  In the underdoped
regime, there is a region of coexisting AF order and superconductivity \cite{ewiesenmayer11}.
Incommensurate spin excitations appear for $x\ge 0.4$ \cite{castellan} and persist till $x=1$
at KFe$_2$As$_2$ \cite{chlee11}. (b) Phase diagram for electron-doped
BaFe$_{2-x}$Ni$_x$As$_2$ \cite{hluo12}.  The long commensurate AF order changes into short-range incommensurate
AF order near $x=0.092$. The inset shows the transverse incommensurate AF order. Superconductivity
in the electron-doped materials only extends to $x\approx 0.25$.
(c) Schematics of FS's correspond to 35\%
hole-doped BaFe$_2$As$_2$ \cite{chenglinzhang} with possible nesting vectors marked with arrows.
The $d_{xz}$, $d_{yz}$, and $d_{xy}$ orbitals for different Fermi
surfaces are colored as red, green, and blue, respectively. (d) FS's of BaFe$_2$As$_2$ with
orbital characters \cite{graser10}. (e) Fermi surfaces of 8\% electron-doped
BaFe$_2$As$_2$ \cite{chenglinzhang}.  For all three cases,
FS's are plotted with zero wave vector transfers along the $c$-axis.
(f) Longitudinally elongated spin excitations at $E=20$ meV seen in the
optimally hole-doped Ba$_{0.67}$K$_{0.33}$Fe$_2$As$_2$ \cite{chenglinzhang}.
(g) Transversely elongated spin waves at $E=20$ meV for BaFe$_2$As$_2$  \cite{lharriger}.
(h) Transversely elongated spin excitations at $E=20$ meV for BaFe$_{1.9}$Ni$_{0.1}$As$_2$ \cite{jtpark,msliu12}. (i) Energy dependence of $\chi^{\prime\prime}(\omega)$ for BaFe$_2$As$_2$ (blue solid line)
and BaFe$_{1.9}$Ni$_{0.1}$As$_2$ below (red dashed line) and above (red solid circles) $T_c$ in absolute units of
$\mu_B^2$eV$^{-1}$f.u.$^{-1}$. The sharp peak near $E\approx 8$ meV below $T_c$ is the neutron spin resonance coupled directly to superconductivity \cite{lumsden,schi09,dsinosov09,jtpark,clester,hfli,hqluo12}.
}
\end{figure}

In the case of insulating $A$Fe$_{1.6+x}$Se$_2$, spin waves have an acoustic mode
and two optical modes separated by spin gaps (Fig. 1l) \cite{mywang11}.
In contrast to iron pnictide $A$Fe$_2$As$_2$ \cite{diallo09,jzhao,raewings,lharriger} and
iron chalcogenide Fe$_{1+y}$Te \cite{lipscombe,zaliznyak}, spin waves in insulating $A$Fe$_{1.6+x}$Se$_2$
can be well-described by a Heisenberg Hamiltonian with NN, NNN, and NNNN exchange couplings \cite{mywang11}. Comparing effective exchange couplings for different iron-based materials (Table I),
it is clear that the NN exchange couplings are quite different, but the NNN exchange couplings are AF and rather similar. In addition, spin waves for iron-based materials are much broader at high energies.
This is different from the insulating copper oxides, where the NN exchange coupling dominates the magnetic interactions and spin waves are instrumental resolution limited throughout the Brillouin zone \cite{coldea,headings}.  These results suggest that itinerant electrons play a role in spin waves of metallic iron-based materials.

\begin{table}[h]
\begin{center}
     \begin{tabular}{ |l | l | l | l | l | l | l |}
     \hline
           Parent compounds  & $T_N$ (K)    & $SJ_{1a}$ (meV)       & $SJ_{1b}$ (meV)       & $SJ_{2a}$ (meV)         & $SJ_{2b}$ (meV) &  $SJ_{3}$ (meV) \\ \hline
La$_2$CuO$_4$, Ref. \cite{coldea} & $317\pm 3$ & $111.8\pm4$ & $111.8\pm4$ & $-11.4\pm 3$ & $-11.4\pm 3$  & 0 \\ \hline
CaFe$_2$As$_2$, Ref. \cite{jzhao} & $\sim$170 & $49.9\pm 9.9$ & $-5.7\pm4.5$ & $18.9\pm 3.4$ & $18.9\pm 3.4$ & 0 \\
\hline
BaFe$_2$As$_2$, Ref. \cite{lharriger} & $\sim$138 & $59.2\pm2.0$ & $-9.2\pm 1.2$ & $13.6\pm 1$ & $13.6\pm 1$ & 0 \\
\hline
Fe$_{1.05}$Te, Ref. \cite{lipscombe} & $\sim$68 & $-17.5\pm 5.7$ & $-51.0\pm 3.4$ & $21.7\pm 3.5$ & $21.7\pm 3.5$ & $6.8\pm 2.8$ \\
\hline
Rb$_{0.89}$Fe$_{1.58}$Se$_2$, Ref. \cite{mywang11} & $\sim$475 & $-36\pm 2$ & $15\pm 8$ & $12\pm 2$ & $16\pm 5$ & $9\pm 5$ \\
\hline
     \end{tabular}
     \caption{\label{tab:5/tc}
Comparison of effective magnetic exchange couplings for parent compounds of copper-based
and iron-based superconductors obtained by fitting spin waves with a Heisenberg Hamiltonian with
NN ($J_{1a}, J_{1b}$), NNN ($J_{2a},J_{2b}$), and NNNN ($J_3$).
The N$\rm \acute{e}$el temperatures for different materials are also listed.
     }

 \end{center}
\end{table}

\section{The effects of hole and electron doping on the magnetic correlations and excitations}

As discussed before  \cite{johnston,stewart,paglione10}, superconductivity
in Fe-based materials can be induced via electron/hole doping, pressure, and isoelectronic substitution.
Figures 2a and 2b show the electronic phase diagrams of
hole and electron  doping on BaFe$_2$As$_2$, respectively.
In the undoped state, BaFe$_2$As$_2$ exhibits simultaneous structural and magnetic phase transitions below $\sim$138 K, changing
from the high-temperature paramagnetic tetragonal phase to the low-temperature orthorhombic phase with the collinear AF structure (Fig. 1b) \cite{qhuang}.
Upon electron-doping BaFe$_2$As$_2$ by partially replacing Fe by Co or Ni to form BaFe$_{2-x}T_x$As$_2$, the static AF order is suppressed and
superconductivity emerges \cite{johnston,stewart,paglione10}.  From systematic transport and magnetic measurements of
single crystals \cite{nni08,jhchu09}, the phase diagram for BaFe$_{2-x}$Co$_x$As$_2$ was established, where the single
structural/magnetic phase transition in BaFe$_2$As$_2$ splits with increasing Co-doping.  Neutron diffraction experiments
on BaFe$_{2-x}$Co$_x$As$_2$ \cite{lester09} confirm
that the commensurate AF order appears below the structural transition
temperature and superconductivity coexists with AF order for $0.06\leq x\leq 0.102$.
Neutron scattering measurements on BaFe$_{2-x}$Co$_x$As$_2$
with coexisting AF order and superconductivity reveal that the intensity of AF Bragg peaks
actually decreases below $T_c$ without changing the spin-spin correlation lengths \cite{pratt09,adchristianson09}.  While these results indicate
that the static AF order competes with superconductivity, it remains unclear whether the long-range AF order truly coexists microscopically
with superconducting regions \cite{mywang,mywang11b}.  Recently,
for electron-doped samples near optimal superconductivity it has been shown that
 the commensurate
static AF order changes into transversely incommensurate short-range AF order that coexists and competes with superconductivity
(see inset in Fig. 2b) \cite{dkpratt11,hluo12}.  Taking the temperature dependence of the orthorhombic
lattice distortion of BaFe$_{2-x}$Co$_x$As$_2$ into account \cite{snandi10}, the AF order, structure, and
superconductivity phase diagrams for BaFe$_{2-x}T_x$As$_2$ are shown in Fig.~2b.

Although the superconducting transition temperature for hole-doped
Ba$_{1-x}$K$_x$Fe$_2$As$_2$ can reach up to $T_c=38$ K \cite{rotter} as compared to the
$T_c\approx 25$ K for electron-doped BaFe$_{2-x}T_x$As$_2$ \cite{johnston,stewart}, those materials
are much less studied because of the difficulty in growing
high-quality single crystals. The initial transport and neutron scattering
experiments on powder samples indicated a gradual suppression of the concurrent structural and
magnetic phase transitions with increasing K-doping.  For the underdoped regime
$0.2\leq x\leq 0.4$, commensurate AF order appears to microscopically coexist with
superconductivity \cite{hchen09}.  Subsequent neutron scattering
and muon spin rotation ($\mu$SR) measurements on
single crystals grown in Sn-flux suggested mesoscopic separation of the AF and superconducting phases \cite{jtpark09}.  However, recent
neutron \cite{savci}, X-ray scattering, and $\mu$SR work \cite{ewiesenmayer11} on high-quality powder samples confirm the microscopic coexistence of the commensurate AF order with superconductivity in the underdoped region between $0.2\leq x\leq 0.3$ and the suppression of the
orthorhombic phase below $T_c$ (Fig. 2a).
Since currently there is no neutron diffraction work on high-quality single crystals of Ba$_{1-x}$K$_x$Fe$_2$As$_2$ grown using FeAs-flux, it is unclear if there are also short-range incommensurate AF order in Ba$_{1-x}$K$_x$Fe$_2$As$_2$ near optimal superconductivity.

The appearance of static incommensurate AF order along the transverse direction of
the collinear AF ordering wave vector ${\bf Q}_{AF}=(\pm 1,0)$
in BaFe$_{2-x}T_x$As$_2$ suggests that such order arises from the electron doping effect of
FS nesting \cite{dkpratt11,hluo12}. Based on
a five-orbital tight-binding model, fitted to
the density functional theory (DFT) band structure for BaFe$_2$As$_2$ \cite{graser10},
there should be five FS pockets with different orbital contributions in the two-dimensional reciprocal space at ${\bf Q}_z=0$ (Fig. 2d).
The intraorbital, but interband, scattering process between $\Gamma(0,0)\leftrightarrow M(1,0)$ shown in Fig. 2d favors
the transversely lengthened vertices \cite{jhzhang10}. This momentum anisotropy is compatible with the experimentally observed
elliptically shaped low-energy spin excitations in superconducting BaFe$_{2-x}T_x$As$_2$ \cite{jtpark,clester,hfli,hqluo12} and spin waves
in BaFe$_2$As$_2$ (Fig. 2g) \cite{lharriger}.  Upon electron-doping to enlarge the electron pockets near
$M(1,0)/(0,1)$ and shrink the hole pockets near $\Gamma(0,0)$,  the mismatch
between the electron and hole Fermi pockets
becomes larger (Figs. 2e), resulting in a more transversely elongated ellipse in the low-energy
 magnetic response (Fig. 2h).  Indeed, this is qualitatively
 consistent with the doping evolution of the low-energy spin excitations \cite{jtpark,hqluo12,msliu12}.

For hole-doped Ba$_{1-x}$K$_x$Fe$_2$As$_2$, one should expect enlarged hole Fermi pockets near $\Gamma(0,0)$ and reduced electron pockets near
$M(1,0)/(0,1)$, as shown in Fig. 2c.
Based on first principles calculations, spin excitations for optimally hole-doped
Ba$_{1-x}$K$_x$Fe$_2$As$_2$ at $x=0.4$ should have longitudinally elongated ellipses \cite{jtpark},
and gradually evolve into incommensurate magnetic scattering (elastic and/or inelastic) with increasing $x$
due to poor nesting between the hole and electron Fermi pockets \cite{castellan}.
INS experiments on single crystal
Ba$_{0.67}$K$_{0.33}$Fe$_2$As$_2$ \cite{chenglinzhang} indeed confirm that the low-energy spin excitations
are longitudinally elongated ellipses that are rotated 90$^\circ$ from
that of the electron-doped BaFe$_{2-x}T_x$As$_2$ (Fig. 2f) \cite{jtpark,clester,hfli,hqluo12}. Furthermore, INS
measurements on powder samples of Ba$_{1-x}$K$_x$Fe$_2$As$_2$ reveal that spin excitations change from commensurate
to incommensurate for $x\ge 0.4$, although their exact line shape and incommensurability in reciprocal space are unknown \cite{castellan}.
Finally, INS experiments on hole-overdoped KFe$_2$As$_2$ found incommensurate spin fluctuations along
the longitudinal direction (inset in Fig.~2a), again consistent with the FS nesting picture \cite{chlee11}.
Figure 2a shows the electronic phase diagram of hole-doped Ba$_{1-x}$K$_x$Fe$_2$As$_2$ based on the present understanding of these materials.

Although FS nesting is compatible with a number of experimental observations
for the evolution of spin excitations in electron/hole-doped iron-based
superconductors, there are several
aspects of the problem where such a scenario cannot be reconciled with experiments.
In a recent INS experiment on optimally electron-doped
BaFe$_{1.9}$Ni$_{0.1}$As$_2$, magnetic excitations throughout the Brillouin zone have been measured in absolute units and compared with spin waves for
AF BaFe$_2$As$_2$  \cite{msliu12}. In the fully localized (insulating) case, the formal
Fe$^{2+}$ oxidation state in BaFe$_2$As$_2$ would give a $3d^6$ electronic configuration and
Hund's rules would yield $S=2$.  The
total fluctuating moments should be $\left\langle m^2\right\rangle=(g\mu_B)^2
S(S+1)=24\ \mu_B^2$ per Fe assuming $g=2$ \cite{clester,msliu12}.  For spin waves in the insulating Rb$_{0.89}$Fe$_{1.58}$Se$_2$, the total moment sum rule
appears to be satisfied \cite{mywang11}.  The fluctuating moments for BaFe$_2$As$_2$ and BaFe$_{1.9}$Ni$_{0.1}$As$_{2}$ are  $\left\langle
m^2\right\rangle=3.17\pm 0.16$ and $3.2\pm 0.16\ \mu_B^2$ per Fe(Ni), respectively \cite{msliu12}.  While these values
are considerably smaller than those of the fully localized case,
they are much larger than expected from the fully itinerant SDW using the random phase approximation \cite{hpark12}. A calculation combining DFT
and dynamical mean field theory (DMFT)
suggests that both the band structure and the local
moment aspects (e.g. Hunds coupling) of the iron electrons are needed for a good description of the magnetic responses \cite{msliu12}.  Figure 2i
shows the energy dependence of $\chi^{\prime\prime}(\omega)$ for BaFe$_2$As$_2$ and BaFe$_{1.9}$Ni$_{0.1}$As$_{2}$, and it is clear that
the impact of electron doping and superconductivity are limited to spin excitation energies below 100 meV.  These results suggest that high-energy spin excitations are likely to arise from the
local moments instead of FS nesting effects.

\section{Deviations from the simple SDW Fermi surface nesting picture}

After the early research efforts on Fe-based
superconductors~\cite{johnston,stewart,paglione10}, recent experimental
and theoretical investigations are providing a more refined perspective of
these materials. Below,
several selected examples will be discussed, supplementing
those presented in the neutrons sections.

{\it Strength of electronic correlations}. The strength of electronic correlations
are often characterized via the ratio
between the on-site Hubbard repulsion coupling $U$ and the bandwidth $W$ of the
hole or electron carriers.
Early on, it was assumed that pnictides were in the weak-interaction
limit $U/W\ll 1$. However, recent investigations revealed that the electronic
correlations induce large enhancements between the
effective and bare electronic masses, signaling that
correlation effects cannot be neglected. For instance,
Haas-van Alphen experiments for
KFe$_2$As$_2$ unveiled discrepancies between the band-theory
calculated and observed FS's, including a large electronic mass
enhancement 3-7 caused by band narrowing \cite{tterashima10}. Similar ratios
for the overdoped Tl$_2$Ba$_2$CuO$_{6+\delta}$ copper-oxides have been
reported~\cite{pmcrourke}, suggesting that the undoped parent compounds of the pnictides
resemble the overdoped copper oxides.

Additional insight is provided by optical conductivity experiments,
since the ratio $R$ between the experimentally measured kinetic energy
and that of band-theory calculations can be measured and
contrasted against other compounds \cite{qazilbash09}. $R\approx 1$ signals
a good metal such as Ag. LaFePO presents a ratio $R\approx 0.5$ which is borderline
between weak and moderate coupling. However, pnictides such as BaFe$_{2-x}T_x$As$_2$
are characterized by an even stronger correlation that induces a ratio $R\approx 0.3$
which is similar to results for overdoped La$_{2-x}$Sr$_x$CuO$_4$, widely considered to be a ``correlated metal''.
Other studies have arrived to similar conclusions with regards
to the correlation strength \cite{nakamura,dsinosov11}. In agreement
with experiments, DFT+DMFT predicts a mass enhancement
$m^\ast/m_{band}\sim$2-3 for BaFe$_{2-x}T_x$As$_2$ and $\sim$7
for FeTe \cite{zpyin10}. Moreover, ARPES studies
of NaFeAs revealed band reconstructions in the magnetic state involving
bands well below the FS~\cite{che10}, contrary to a weak coupling picture.

\begin{figure}[t]
\includegraphics[scale=.48]{Fig3R2}
\caption{Summary of the phase diagram of multiorbital Hubbard models
and the electronic state of Fe near the FS.
(a) Sketch of the phase diagram of a typical multiorbital Hubbard model
in the undoped limit, varying the on-site same-orbital repulsion $U$ and
the ratio between the Hund coupling $J_H$ and $U$.  Highlighted is a region
dubbed ``physical region'' where the properties of the model are in good
agreement with experiments. Note the location of this region in the intermediate
magnetic-metallic phase, with magnetic order at ${\bf Q}_{AF}=(1,0)$, at similar
distance from the paramagnetic state and from the insulator orbitally-ordered state \cite{qluo10}.
(b) Sketch of the DOS illustrating the phenomenon of FS orbital order, which
is a weight redistribution at the FS of the states associated with the $xz$
and $yz$ $d$-orbitals. Even though the integral over energy gives similar values
for both orbitals, at the FS there are drastic differences that influence
on several properties such as transport \cite{daghofer10,daghofer12}.
(c) Sketch of the anisotropy found in transport experiments for detwinned
Ba(Fe$_{1-x}$Co$_x$)$_2$As$_2$. Note that this anisotropy is present
at temperatures substantially larger than $T_N$ \cite{irfisher11}.
(d) Orbitals of relevance for the discussion of the Fe-based superconductors
and their splitting at the FS.
(e) Sketch of the ARPES results of Ref. \cite{shimojima11}, illustrating
the absence of a nesting partner for one of the hole pockets. The material
still displays a nearly uniform superconducting gap at this and all the other
hole and electron pockets.
(f) Sketch of the magnetic moment at wave vector ${\bf Q}_{AF}=(1,0)$ for
two models. On the left is the result for a traditional model of pnictides,
with the $xz$, $yz$, and $xy$ $d$-orbitals active at the FS.
This model displays magnetic order in a broad range of couplings from very
weak to strong. On the right, results of a model with the same FS but totally
different orbital composition. While at small
$U$ there is no order, at larger couplings this model
converges to the same ${\bf Q}_{AF}=(1,0)$ order \cite{nicholson}.
}
\end{figure}

Hubbard model investigations provide additional insight
on this subject. When compared with similar efforts
for the cuprates, the study of Hubbard models for the pnictides
is far more challenging because several Fe orbitals are needed. For this reason, many efforts are restricted to
mean-field Hartree-Fock approximations. For the undoped three-orbital
Hubbard model, employing the $d_{xz}$, $d_{yz}$, and $d_{xy}$ orbitals
of relevance at the FS, a sketch of a typical mean-field
phase diagram varying $U$ and the Hund coupling $J_H$ \cite{qluo10} is in Fig.~3a.
Three regimes are identified: a small $U$ phase where the state is paramagnetic,
followed with increasing $U$ by an intermediate regime
simultaneously metallic and magnetic~\cite{qluo10}, and finally a large-$U$ phase
where a gap in the
density-of-states (DOS) is induced leading to an insulator (with concomitant orbital order). Comparing
the theoretical predictions for the magnetic moment in the ${\bf Q}_{AF}=(1,0)$
wave-vector channel against neutrons, and the one-particle spectral function
$A({\bf k},\omega)$ against ARPES, the intermediate-coupling region dubbed
``physical region'' in yellow in Fig.~3a represents qualitatively
the undoped BaFe$_{2-x}T_x$As$_2$ compounds \cite{qluo10}. In this regime,
$U/W\sim$0.3-0.4, and similar results were reported for the two- and five-orbital models~\cite{mDaghofer}.
Note that Hartree-Fock usually produces critical couplings
smaller than they truly are because of the neglect of quantum fluctuations.
In fact, recent investigations beyond Hartree-Fock ~\cite{kkubo11}
suggest that the relevant $U$ may be larger than those found in Hartree-Fock \cite{qluo10} by
approximately a factor two. The study of
effective low-energy Hamiltonians starting from first-principles calculations
also led to the conclusion that $U/W$ is between 0.5 and 1.0 for the pnictides depending on the particular
compound~\cite{arita-imada}. Thus,
the regime of relevance is neither very weak coupling nor strong
coupling but the more subtle, and far less explored, intermediate  region.
Previous efforts converged to similar conclusions~\cite{Johannes09}.
This is also compatible with the notion
that the parent compound is close to a Mott insulator~\cite{si2008,qmsi09}.
In the ``physical region'' the ratio $J_H/U$ is approximately 1/4 ~\cite{qluo10},
as in other estimations \cite{zpyin10}, highlighting the
importance of $J_H$ in these materials that are sometimes referred to as Hund metals~\cite{hund-metal}.
Finally, it is very important to remark that
the above described analysis of $U/W$ holds for pnictides but the recent discovery of the alkaline
iron selenides \cite{jgguo,mhfang} has opened a new chapter in this field and it is
conceivable that for these materials $U/W$ will be larger than in pnictides explaining, for
example, the large values of the iron moments.

{\it Role of the orbital degree of freedom}.
The ``physical region'' in Fig.~3a is not only close to the paramagnetic regime,
but also similarly close to the insulator, which in the mean-field approximation
is also orbitally ordered \cite{mDaghofer}. The potential relevance of the orbital
degree of freedom in pnictides has been discussed \cite{lv10,wgyin10}.
The orbital can be of relevance not only in its long-range-ordered
form, but also via its coupling to the spin and its influence near the FS.
In fact, polarized ARPES experiments on BaFe$_2$As$_2$ \cite{shimojima} reported that
at the FS there was an asymmetry between the populations of the $d_{xz}$
and $d_{yz}$ orbitals. Theoretical studies showed that this effect indeed
occurs in the ${\bf Q}_{AF}=(1,0)$ magnetic state, and it is linked to an
orbital-dependent reduction in the DOS at
the FS \cite{daghofer10}, sketched in Fig. 3b, phenomenon
dubbed ``Fermi surface orbital order''.

This effect, while not sufficiently strong to induce long-range order
as in manganites, can still severely influence the
properties of the material. Consider for example the transport anisotropy
observed in detwinned BaFe$_{2-x}T_x$As$_2$ single-crystals \cite{irfisher11,matanatar},
sketched in Fig. 3c. At low temperatures the difference between the
$a$-axis (spins antiparallel, Fig. 1b) and $b$-axis (spin parallel, Fig. 1b)
directions can be rationalized based on the magnetic state, since
the different spin arrangements along the $a$ and $b$ break rotational
invariance~\cite{xtzhang}. However, both in the undoped case and particularly
in the lightly-doped regime, the asymmetry persists well above the N$\rm \acute{e}$el
temperature, $T_N$, into a new temperature scale $T^\ast$ that may be associated
with the onset of nematic order \cite{cfang,cxu}, similarly as
in some ruthenates and copper oxides \cite{efradkin}.  ARPES experiments
on the same materials \cite{myi} reported a $d_{xz}$ and $d_{yz}$ band
splitting (Fig.~3d) that occurs above $T_N$ in the same region where transport
anisotropies were found. Although the splitting is too small to be a canonical
long-range orbital order, it reveals the importance of fluctuations above the
critical temperatures. Optical spectra studies also unveiled anisotropies in
the spectra persisting up to 2 eV, incompatible with SDW
scenarios \cite{Nakajima}. Note that the discussion on this subject is still
fluid. While neutron diffraction investigations showed that $T_N$ actually
substantially increases as the pressure needed to detwin the crystals increases,
potentially explaining the observed resistivity anisotropies \cite{dhital},
magnetic torque measurements without external pressure revealed clear evidence for
 electronic nematicity \cite{kasahara}.
Recent calculations addressing transport indeed find an important role of the orbital
states above $T_N$ \cite{rmfernandes}. The
orbital degree of freedom, closely entangled to the spin and the lattice, may lead
to a more complex ``normal'' state than anticipated from weak coupling
particularly because of the FS orbital order~\cite{daghofer10}.
In fact, neutron scattering shows that although the low-energy magnetic
dispersion changes substantially when crossing critical temperatures, the higher
energy features remain the same over a large doping and temperature range \cite{msliu12},
suggesting that spin, orbital, and lattice are closely entangled. Establishing who
is the ``driver'' and who is the ``passenger'' may define an important area of
focus of future research.

{\it Local moments at room temperature}. Another deviation from a simple weak coupling picture
is the observation of local magnetic moments at room temperature.
Within the SDW scenario, magnetic moments are formed upon
cooling simultaneously with the development of long-range magnetic order.
But recent Fe X-ray emission
spectroscopy experiments unveiled the existence of local moments in the
room-temperature paramagnetic state~\cite{hgretarsson}. In fact, with the only
exception of FeCrAs, for all the pnictides
and chalcogenides investigated a sizable room temperature magnetic moment was found.
This includes LiFeAs, that actually does not order magnetically at any
temperature~\cite{cwchu}, and  $A$Fe$_{1.6+x}$Se$_2$ with
a regular arrangement of Fe vacancies (Fig. 1d). These observed local moments
are similar in magnitude to those reported in the low-temperature neutron
scattering experiments reviewed in previous sections. Similar conclusions
to those of~\cite{hgretarsson} were reached in a study of $3s$ core level
emission for CeFeAsO$_{0.89}$F$_{0.11}$~\cite{bondino} and also in LDA+DMFT
investigations~\cite{hansmann}.

{\it Polarized ARPES results and orbital composition}.
While research using ARPES techniques applied to pnictides has
already been reviewed~\cite{prichard},
 some intriguing recent results addressing the influence of nesting
are included in our discussion. Using bulk-sensitive laser ARPES on
BaFe$_2$(As$_{0.65}$P$_{0.35}$)$_2$ and Ba$_{0.6}$K$_{0.4}$Fe$_2$As$_2$,
an orbital-independent superconducting gap magnitude was found for the
hole-pockets FS's \cite{shimojima11}. These results are
incompatible with nesting where the FS nested portions
must have a robust component of the same orbital to be effective.
Actually, the red hole pocket shown in the sketch in Fig. 3e, that experimentally
displays a robust and nearly wave-vector-independent superconducting gap similar to those
found in the other hole pockets, does not have a matching electron pocket with the
same orbital composition and, thus, it cannot develop its superconductivity via a
nesting pairing mechanism \cite{mazin2011n}. Perhaps inter-orbital
pairing~\cite{amoreo} or orbital fluctuations could be relevant to explain
this paradox. Recent theoretical work~\cite{nicholson} addressed the importance
of orbital composition via two models: one with nested electron- and
hole-pocket Fermi surfaces with the standard orbital composition of pnictide models,
and another with
the same FS shape but with electron and hole pockets having totally different
orbital compositions. As sketched in Fig. 3f, the former develops magnetic order
at smaller values of $U$ than the latter. However,
 with sufficiently large $U$ both have magnetic ground states with the same
wavevector ${\bf Q}_{AF}=(1,0)$ (Fig. 3f). At large $U$ it is clear that
the ${\bf Q}_{AF}=(1,0)$
order can be understood within a local picture,
based on the similar magnitude of the super-exchange
interactions between NN and NNN spins using a simple
Heisenberg model.

{\it Additional experimental results}. De Haas-van Alphen
studies \cite{bjarnold} in non-superconducting BaFe$_2$P$_2$, the end member of the
superconducting series BaFe$_2$(As$_{1-x}$P$_x$)$_2$, indicate that the differences
in the pairing susceptibility varying $x$ are caused by increases in
$U$ and $J_H$ rather than improved geometric
nesting. Moreover, ARPES studies of LiFeAs, without
long-range magnetic order at low temperatures, report a strong renormalization
of the band structure by a factor $\sim$3 and the absence of
nesting \cite{Borisenko}. Yet, at $T_c=18$ K \cite{cwchu} LiFeAs
still becomes superconducting suggesting that nesting is not necessary for
superconductivity to develop. Similarly, ARPES experiments on
superconducting $A$Fe$_{1.6+x}$Se$_2$ \cite{yzhang,tqian,Dxmou} revealed
the absence of the hole-like FS's necessary
for the $\Gamma(0,0)\leftrightarrow M(1,0)$
$s^{\pm}$-wave superconductivity.
Also note that related materials such as LaFePO with a well-nested
FS also do not order magnetically. Why weak coupling arguments would work in some
cases and not others? Finally, scanning tunneling microscopy (STM)
experiments~\cite{tmchuang} on Ca(Fe$_{1-x}$Co$_x$)$_2$As$_2$
shows an exotic ``nematic'' electronic structure,
similar to those found for strongly coupled copper-oxides.

{\it Additional theoretical results}. In
fluctuation-exchange approximation studies it was concluded
that the nesting results are not robust against the addition of self-energy
corrections \cite{arita09}. Other calculations have suggested
that magnetic order in pnictides is neither fully localized nor fully itinerant:
the $J_H$ coupling forms the local moments, while the particular ground state
is selected by itinerant one-electron interactions~\cite{Johannes09}. Moreover,
studies of a spin-fermion model for the pnictides~\cite{lv10,wgyin10,shliang11}
revealed the crucial role played by the Hund's rule coupling and suggested that
the Fe superconductors are closer kin to manganites, where similar
spin-fermion models were extensively studied~\cite{dagotto-cmr},
than to copper-oxides with regards to their diverse magnetism and incoherent
normal-state electron transport.

\begin{figure}[t]
\includegraphics[scale=.48]{Fig4R2}
\caption{
Sketch of the expected phase diagram of the Hubbard model varying temperature
and $U/W$ in the undoped limit. Highlighted are the regimes of weak coupling
where nesting dominates and strong coupling where localized spins approaches
are suitable. At temperatures above $T_N$, there are regions with and without preformed
local moments. The dashed line tentatively locates the pnictides and chalcogenides
in the ``middle'', with a physics involving itinerant electrons coexisting with localized
moments. Note that this phase diagram is guided by results
known for the one-orbital case, while the true multiorbital Hubbard model phase
diagram may display an even richer structure. In particular, a second critical $U/W$
at low temperatures separating the metallic AF state from the insulating AF state
is not shown for simplicity.}
\end{figure}

\section{Conclusions}

Recent studies of Fe-based superconductors
are unveiling a perspective of these exciting materials
that is far richer than previously anticipated. While in the early days weak coupling
approaches seemed sufficient to understand
these compounds, several recent efforts, reviewed in part here, suggest
that understanding the physics of these materials may require more
refined concepts, better many-body theoretical calculations, and additional
sophisticated experiments for a more in-depth rationalization
of their properties. In fact, evidence
is building that pnictides and chalcogenides inhabit the mostly unexplored
``intermediate'' region of Hubbard $U/W$ couplings, which is neither very weak
coupling, where FS nesting concepts apply, nor strong coupling, where localized
spins provide a good starting point as it occurs in the undoped copper oxides.
The situation is qualitatively summarized
in Fig. 4 where a crude sketch of a plausible phase
diagram for a generic undoped Hubbard model is provided varying temperature $T$ and
$U/W$ at, e.g., a fixed $J_H/U$ such as 1/4. In weak coupling, first a critical
value of $U$ must be crossed before magnetic order develops at low temperatures.
In this region, nesting works properly. As $U$ increases, $T_N$ first
increases, reaches
a broad maximum, and then eventually in the regime of localized spins $T_N$ starts to
decrease since it becomes regulated by the Heisenberg superexchange that scales
as $1/U$. Above $T_N$ a ``crossover'' temperature that roughly grows like $U$ is
shown separating regions with and without ``preformed'' local moments. Since the
pnictides have local moments at room temperature,
then a tentative location
for these materials is provided by the dashed line.
However, whether this line coincides with the
maximum $T_N$ or is shifted to the left or the right is too early to say, but it
cannot be too far from optimal otherwise local moments would be absent, if far left,
or an insulator should be found at low temperatures, if far right.
Theoretical mean-field
estimates reviewed here using the multiorbital Hubbard model find that $U/W$$\sim$$0.3-0.5$
could work for pnictides. However, for chalcogenides and alkaline iron selenides,
and also after including quantum fluctuations, the ratio $U/W$ may increase further,
and it may reach the $U/W\sim1$ threshold widely consider to mark the starting point for a strong coupling
description. Note also that the
sketch in Fig.~4 is based on our knowledge on the one-orbital
Hubbard model and a proper multiorbital analysis will lead to an even richer phase diagram.
In fact, a critical $U$ for the transition between the magnetic metallic state
and the magnetic insulating state at low temperatures
should also be present, but it is
not shown in the sketch for simplicity: this transition should occur at a
$U$ larger than the pnictides dashed line since these materials are metallic at low temperatures.

In summary, the Fe-based superconductors continue surprising us with their
exotic properties that do not fit into the simple limits of weak or strong coupling $U$.
Additional experimental and theoretical efforts are
needed to unveil the secrets of this intriguing family of materials.

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