Quasar Structure Emerges from the Three Forms of Radiation Pressure

Abstract: All quasar spectra show the same atomic features in the optical, UV, near-IR and soft X-rays over all of cosmic time, luminosity black hole mass and accretion rate. This is a puzzle. Here I show that it is possible that all of these atomic features can be accounted for by gas from an accretion disk driven by the three forms of radiation pressure: electron scattering, line driving and dust driving. The locations where they successfully drive an escaping wind, and those where they produce only a failed wind are both needed to produce the observed features.

\section{The Puzzle of Quasar Spectra}
It is remarkable that, to a good approximation, the spectra of all quasars and
active galactic nuclei (called quasars hereinafter) are the same over 13~Gyr of
cosmic time, 6 orders of magnitude in luminosity, and 3 orders of magnitude in
both black hole mass and Eddington ratio. This constancy poses a
problem. Whatever structures produce these atomic features they must be formed
by some robust physics that does not depend strongly on any of these parameters.
Accretion disk winds hint at a solution, but whether this can be realized has
not been clear.

The features are these: All show (1) permitted broad (FWHM$\sim$ 5000
km/s) emission lines (BELs); (2) narrow (FWHM$\sim$ 1000 km/s) emission lines
(NELs) that often extend into resolved bi-cones; (3) an Fe-K X-ray emission line
at 6.4 keV. Half or more show: high ionization (4) UV narrow (FWHM $<$200 km/s)
absorption lines (NALs) and (5) X-ray narrow (FWHM $<$300km/s) absorption lines
(known as 'Warm Absorbers, or "WAs"'), while about 15\% show (6) broad
($\Delta$v$\sim$ 10,000 km/s) absorbers, known as Broad Absorption Lines, or
"BALs").

There are two exceptions: in both "type 2" quasars and "blazars" some
or all of these features are covered up, either by reddening of the central
region (type 2s), or by overwhelming them with a beamed continuum from a
jet (blazars). But in type 2s the missing BELs can be seen in polarized light,
while in blazars a normal quasar spectrum can be glimpsed when the jet is
temporarily in a low state. So we believe that, intrinsically, all quasars have
the same features in their inner regions.

Here I show that it is possible that these atomic features can be completely
accounted for by the three forms of radiation pressure: electron scattering,
line driving and dust driving, both in the locations where they are effective in
driving an escaping wind, and where they produce only a failed wind.

\section{Funnel Wind Model for Quasars}

In 2000 I proposed a phenomenological model for quasar structure that could
account for all the atomic features in quasars (\cite{2000ApJ...545...63E},
'E00'). The model determined the geometry and kinematics of the quasar.
The model posits gas thrown off the accretion disk vertically in a narrow band
of radii. As this wind is hit by the continuum radiation it bends radially
outward, forming a funnel-shaped flow. BALs are seen looking along the radial
flow, which reaches $\tau_{es}\sim$1 so giving the Fe-K line. NALs and WAs are
seen looking almost perpendicularly through the flow, at low angles to the
accretion disk; viewed from above the wind no absorption features are seen. The
relative numbers of NALs, BALs and no absorption quasars set the angle of the
flow to the disk to be $\sim$30$^{\circ}$ [ignoring any possible equatorial
obscuration (see Elvis 2012, in prep.)]. The BELs are formed from a cool phase
of the wind.

\subsection{Funnel Wind Predictions and Tests}

A good model makes specific, testable, predictions. Table~1 lists predictions of
the E00 model, and the papers which tested them. In each case the model
survived the test, though some are still debated.

\begin{table}
\caption{Predictions and Tests of the Elvis (2000) Funnel Wind Model}
\begin{small}
\begin{tabular}{|r|l|}
\hline
Ref \# & Prediction  \\
\hline
1 & WA is in outflow \\
2 & WA has narrow absorption lines \\
3 & WAs \& NALs have same outflow velocities  \\
4 & WAs \& NALs have consistent ionization states \& column densities  \\
5 & WAs \& NALs occur in the same objects  \\
6 & NAL bi-cones will be hollow and matter bounded outflows \\
7 & WAs \& NALs are common in high luminosity objects  \\
8 & WAs \& NALs are common in edge-on AGNs  \\
9 & WA has a few distinct phases  \\
10& WA phases are in pressure balance  \\
11& WA is radially thin  \\
12& NALs arise at accretion disk radii  \\
13& WAs arise at accretion disk radii  \\
14& BELs have a large scale height  \\
15& BELs are dominated by rotation  \\
16& BAL regions exist in all quasars \\
17& BAL region rotates  \\
18& BAL scattering region creates narrow Fe-K emission  \\
\hline
\end{tabular}

{\footnotesize References for each test: (1) \cite{2000ApJ...535L..17K};
(2)\cite{2001ApJ...557....2C};
(3) \cite{2005ApJ...631..741G},
(4) \cite{1994ApJ...434..493M}, \cite{2000ApJ...536..101H},
    \cite{2003ApJ...597..832K};
(5) \cite{2001ASPC..224...45K};
(6) \cite{2000ApJ...532L.101C}, \cite{2011ApJ...727...71F};
(7) \cite{2005A&A...432...15P}, \cite{2003ApJ...599..116V};
(8) \cite{1997ApJ...489L..25L};
(9) \cite{2010ApJ...711..888A};
(10) \cite{2005ApJ...620..165K};
(11) \cite{2007ApJ...659.1022K};
(12) \cite{2005ApJ...623...85G}, \cite{2005ApJ...620..665A};
(13) \cite{2007ApJ...659.1022K}, \cite{2010ApJ...709..611D};
(14) \cite{2000ApJ...536..284K};
(15) \cite{2005MNRAS.359..846S};
(16) \cite{1999MNRAS.303..227Y};
(17) \cite{2007Natur.450...74Y};
(18) \cite{2008MNRAS.388..611S}.
}
\end{small}
\end{table}

The first 3 predictions were quickly confirmed by {\em Chandra} and {\em
XMM-Newton} grating spectra.  The 4th required the realization that partial
covering applied to the UV lines, while the 5th was demonstrated by explicitly
testing if NALs predicted WAs and vice versa. The 6th was demonstrated by {\em
Hubble} STIS long slit spectroscopy of narrow line region bi-cones. Number 7 was
shown with larger UV and X-ray samples, while \#~8 had already been shown, but
was not widely known.

Predictions~9 and 10 have been controversial. It is now clear that every good
X-ray spectrum {\em can} be fitted with a 2-3 phase WA in pressure equilibrium,
but more complex solutions are allowed.  Higher spectral resolution is needed to
be definitive.

Tests 11-13 are closely connected. Initial results showed small radii, but a
number of newer observations clearly put certain NAL systems at large,
kiloparsec, distances. This is the greatest challenge so far to the funnel wind
model. How general this conclusion is remains unknown.

BELs (\#14-15) have been less tested. Polarimetry indicates rotation in some
bright objects. BALs (\#17-18) are also poorly tested. Again, polarimetery shows
that BAL velocity gas is present in non-BAL objects and that the BAL region
rotates.

These results gave us the confidence to search for a physical understanding
based on this structure, one that can explain the more peculiar features in E00:
the thin wind region, and a flow neither polar nor equatorial.

\section{A Non-hydrodynamic Radiation Driving model}

Radiation can transfer momentum to matter in only three ways\footnotemark, via:
(1) electron scattering, (2) atomic absorption, (3) molecular/surface physics
absorption.  Normally we talk of these as Compton scattering, line driving and
dust driving. Below the Eddington limit gravity exceeds the Compton scattering
force, preventing a Compton scattered wind escaping; but few quasars reach
Eddington (\cite{2006ApJ...648..128K},\cite{2010MNRAS.402.2637S}).  Line driving
is hundreds of times more effective in O-stars, driving escaping winds well
below Eddington. A dust grain effectively absorbs all the incident light, and so
a medium with sufficient charged dust can even more effectively accelerate a
wind. Which of these mechanisms work in quasars?
\footnotetext{Neglecting pair production.}

\subsection{The QWIND Code}
We built a {\em non-hydrodynamic} code, QWIND (\cite{2010A&A...516A..89R},
'RE10').  This is justified as the gas is always moving supersonically. NAL/WA
blueshifts are $\sim$1000~km/s and the gas, at T$\le$10$^6$K, has a thermal
velocity of $\le$100~km/s, so this assumption is reasonable. Avoiding hydro let
us quickly explore parameter space.

QWIND assumes that packets of gas are ejected at {\em all} radii from the disk,
a generalization of E00, with some velocity and density, and follows their
equation of motion as they are irradiated (as in \cite{1980AJ.....85..329I},
\cite{1999PASJ...51..725W}).
CLOUDY is used to determine the ionization of the gas packet, and performs the
radiative transfer to filter the quasar continuum through the inner gas
packets. We used the \cite{1975ApJ...195..157C} 'CAK') formalism to determine
$\alpha(CAK)$, the multiplier, or gain, of line driving over pure Compton
scattering.  Importantly, the UV continuum has a cos$\theta$ dependence on angle
above the disk, while the X-ray continuum is considered isotropic. The UV/X-ray
ratio is an input parameter. Only the supermassive black hole gravity was
included.

For typical conditions QWIND produced the flow lines in Figure~1. A number of
distinct zones form naturally.

\begin{figure}
\centering
\includegraphics[width=2.5in, angle=-90]{ELVIS_Figure1.ps}
\caption{\small Flowlines of gas from a quasar accretion disk generated by QWIND
  (annotated and expanded from \cite{2005ApJ...630L.129R}. Note the axes are
  linear in units of Schwartzchild radii.}
\end{figure}

\subsection{Compton Scattering Zone}
The innermost zone is overionized, so only Compton scattering is effective, but
below L(Edd) the wind cannot escape. Once $\tau_{es}>1$ this gas preferentially
scatters X-rays out of the disk plane, reducing the ionization parameter in the
next zone. This is the "hitchhiking gas" of \cite{1998ApJ...494..125M}). The gas
in this region still reaches a few 1000 km/s) before falling back. Such high
ionization, turbulent, variable velocity gas has been seen in NGC 1365
(\cite{2005ApJ...630L.129R}), with $\tau_{es}\sim$1.  Some L(Edd) objects show
high velocity winds (\cite{2003MNRAS.345..705P}, \cite{2002ApJ...579..169C}).

\subsection{Line Driving Zone}
The next zone out is where line driving becomes effective and an escaping wind
is formed.  Line driving has become widely accepted for NALs/WAs, following
Murray et al. (1995, \cite{1995ApJ...451..498M}). We find that in this zone
$\alpha(CAK)$ is only $\sim$10. The escaping wind still reaches BAL-like speeds
exceeding 10,000~km/s, but the low $\alpha(CAK)$ allows the wind to make a large
angle to the disk, contrary to Murray \& Chiang (1995
\cite{1995ApJ...454L.105M}), and consistent with E00.  Once the column
density reaches N$_H\sim$10$^{22}$cm$^{-2}$, the line driving UV photons are
absorbed, so this zone is physically thin, as in E00.  This N$_H$ is
characteristic of WAs.  High ionization BELs (CIV, OVI, NV, HeII, 'HiBELs') are
often considered part of a line-driven wind, and can produce the observed line
profiles (\cite{1997ApJ...474...91M}). Some multi-phase WA solutions include
a $\sim$10$^4$K branch, which may produce the HiBELs.  This should be explored
in detail.

\subsection{Failed Wind Zone}
Beyond the line driving zone the gas packets no longer reach escape velocity,
but do achieve a significant scale height.  Gas in this region could be
responsible for the low ionization BELs (FeII, MgII, 'LoBELs'), a modification
of E00. The LoBELs have long been located in the accretion disk
(\cite{1988MNRAS.232..539C}). This needs investigating using the LOC
formalism (\cite{1995ApJ...455L.119B}).  X-ray N$_H$ changes in days are
increasingly common (\cite{2004ApJ...615L..25E},
\cite{2007MNRAS.377..607P}), and have radii and densities suggestive of
LoBEL clouds.

\section{Dust Driven Zone}
RE10 do not include dust driving, as dust should be absent within the dust
sublimation radius, $r_{sub}$. However, \cite{2011A&A...525L...8C} point out
that the pressure and temperature conditions at which AGB star winds produce
dust are also found on quasar disk surfaces, from a few 1000 Schwartzchild radii
outward. Interior to $r_{sub}$ these clouds will initially be accelerated, but
their dust will soon evaporate, and they will cease to accelerate, probably
forming a failed wind. There are changes in X-ray N$_H$ on timescales of years
(e.g. Marinucci et al. in prep.), which might be caused by dusty clouds at these
radii crossing the line of sight. Beyond $r_{sub}$ these clouds will form a
wind, which may be one form of the obscuring "torus" needed to explain type~1
and type~2 quasars (Elvis 2012, in prep.).

\section{Line-Driven Wind Parameter Space}
Varying the parameters in QWIND lets us map out the conditions under which an
escaping wind forms. Clouds with BEL densities and temperatuers make a winds
hard to drive at high black hole mass; while an X-ray quiet continuum produces
winds over a much wider range of conditions, as expected
(\cite{1995ApJ...454L.105M}).

Interestingly high L/L(Edd) does not, as often assumed, imply faster winds.  At
higher L/L(Edd) the wind begins further out, and the 1/r$^2$ weakening of the
continuum lowers the final speed of the wind.

\section{Is Matter Launched at All Disk Radii?}
RE10 explicitly separate the problem of {\em accelerating} the clouds, which is
treated by QWIND, from the initial {\em launching} of clouds, by assuming that
launching occurs at all radii. There is no necessary connection between the two
processes. In fact, it seems unlikely that radiation driving alone can launch
material at all disk radii.

There is evidence that the disk supplies gas at many radii. In NGC~5548, the
H$\beta$ BEL radius 'breathes' in and out by a factor $\sim$3 as the $<$3500\AA\
luminosity changes by a factor $>$6 (\cite{1999ApJ...510..659P}). The few year
timescale is much shorter than the viscous timescale, but comparable to the
dynamical time. Plausibly there is pre-existing gas at all these radii, and the
continuum luminosity picks out the radius where the H$\beta$ emissivity is a
maximum, as in the LOC model, but now applied to a disk geometry. BELs form a
layered sequence from high ionization at small radii out to low ionization at
factor $\sim$20 larger radii (\cite{1999ApJ...521L..95P}, again as though the
gas is pre-existing and emitting whatever BEL is optimal.

The X-ray 'eclipses' by (probable) BEL clouds seen in NGC~1365, have a
comet-like form: small, dense heads with larger, lower density and higher
ionization, tails (\cite{2010A&A...517A..47M}).  The head lifetimes are
$\sim$2 months, substantially shorter than the orbital timescale.  The clouds
must then be constantly replenished. The disk is the only plausible source.

Magnetic reconnection could launch gas at all radii. Solar Coronal mass
ejections reach $\sim$1000~km/s, more than adequate for quasars. The
magneto-rotational instability (MRI), likely responsible for the viscosity in
quasar disks, can produce reconnection events. However, there may be a 'dead
zone' where the gas is too neutral to support MRI, as in proto-planetary disks
(\cite{2011ApJ...740L...6M}). Would there be a corresponding gap in the gas
supply for accelaration? MHD winds in general are another means of accelerating
quasar winds (\cite{2008ApJ...685..160N}). Different acceleration processes
may act in different quasars.  Tests are needed to discriminate between
them.

\section{Conclusions}
There is an elegance to this view of quasar structure. To form so many observed
phenomena from such simple underlying physics is economical. That initially
strange features of the funnel wind model are natural consequences of the onset
and failure of the three forms of radiation driving is particularly appealing.

This physics is robust to changes in parameters, including L/L(Edd) and black
hole mass, as boundary conditions are not important.  Anything beyond the
immediate environs of the black hole does not change the wind physics, so cosmic
time is not important. Hence the puzzle of quasar spectra, with which I began,
can be understood: all quasars are much the same because radiative acceleration
in a disk setting is always the same.

It will now be important to build on this success by seeing if the model can
produce second-order effects, such as eigenvector 1
(\cite{1992ApJS...80..109B}), and whether a predictive model for the energy and
momentum carried by the wind can be produced, in order to have a physical model
of quasar feedback.

\acknowledgements The author would especially like to the thank Guido Risaliti,
and also Bozena Czerny for discussions of her work. He is grateful for the
opportunity to present this work in Charleston.

Cosmology and Nongalactic Astrophysics (astro-ph.CO)


Tile and crystal calorimeter for KLOE-2 experiment

Abstract: The upgrade of the DA$\Phi$NE machine layout requires a modification of the size and position of the inner focusing quadrupoles of KLOE$^2$ thus asking for the realization of two new calorimeters covering the quadrupoles area. To improve the reconstruction of $K_L\to 2\pi^0$ events with photons hitting the quadrupoles, a tile calorimeter, QCALT, with high efficiency to low energy photons (20-300 MeV), time resolution of less than 1 ns and space resolution of few cm, is needed. We propose a tile calorimeter with a high granularity readout corresponding to about 2500 silicon photomultipliers (SiPM) of $1\times 1$ mm$^2$ area. Moreover, the low polar angle regions need the realization of a dense crystal calorimeter with very high time resolution performances to extend the acceptance for multiphotons events. Best candidates for this calorimeter are LYSO crystals with APD readout or PbWO$_4$ crystals with large area SIPM readout.

\section{The KLOE$^2$ proposal}

In the last decade a wide experimental program has
been carried out at Da$\Phi$ne\cite{dafne}, the $e^+e^-$
collider of the Frascati National Laboratories,
running at a center of mass energy of 1020 MeV, the
$\Phi$ resonance mass. During KLOE run, Da$\Phi$ne
delivered a peak luminosity of 1.5$\times$10$^{32}$
cm$^{-2}$s$^{-1}$
which granted about 1 fb$^{-1}$ per year in the last
data taking period.

A new machine scheme has been recently proposed
by the Frascati accelerator group aiming at increasing
the luminosity of the machine up to a factor 5.
This scheme has been succesfully tested at Da$\Phi$ne,
and these encouraging results push for
a new data taking compaign
for the KLOE experiment to complete its physics program
and to perform a new interesting set of measurements.

The new experiment, named KLOE$^2$,
expects to collect 5 fb$^{-1}$/year.
We are now working to
improve the performances
of our detector\cite{kloe_all} adding: an inner tracker, a tagger
system to study $\gamma\gamma$ physics, a new small angle calorimeter
and a new quadrupole
calorimeter.
In this paper we explain the project and the
R\&D for these last two items.

\section{Quadrupole tile calorimeter, QCALT}

In Fig.\ref{figkloe}, we show a section of the KLOE detector in which
is visible the old position of the focalizing quadrupoles and
the surrounding
calorimeters QCAL \cite{oldqcal}
which have a polar angle
coverage of 0.94$\,<|\cos\theta|<\break$0.99.
Each calorimeter consists of  16 azimuthal sectors composed
by alternating layers of 2 mm lead and 1 mm BC408
scintillator tiles, for a total thickness of $\sim$5X$_0$.
The readout is done by wavelength shifter fibers
 (Kuraray Y11-200) and photomultipliers. The fiber arrangement allows
the measurement of the longitudinal coordinate by time differences.
These calorimeters are characterized by a low light response
($\sim$3 pe/mip/tile at zero
distance from the photomultiplier) due
to the  coupling in air, to the fiber lenght
($\sim$2 m for each tile) and
to the quantum efficiency of the used photomultipliers (standard bialkali
with $\sim$20\% QE).

\begin{figure}[htb]
\psfig{file=kloe.eps,width=8cm}
\vspace{-1.5cm}
\caption{Section of the KLOE detector.}
\label{figkloe}
\end{figure}

The project of the new QCAL consists in a dodecagonal structure,
1 m long, covering the region of
the quadrupoles. The structure
consists in a sampling of 5 layers of
5 mm thick scintillator plates alternated
with 3.5 mm thick tungsten plates,
for a total depth of 4.75 cm (5.5 X$_0$).
The active part of each plane is divided into
twenty tiles of 5$\times$5 cm$^2$ area
with 1 mm diameter WLS fibers embedded
in circular grooves. Each fiber is then
optically connected to a silicon
photomultiplier of 1 mm$^2$ area, SiPM,
for a total of 2400 channels.

We have done some R\&D studies on SiPM, fibers and tiles
to choose the combination which optimizes
the response of our system.

\subsection{Test on MPPC}

We have compared the characteristics of two different SiPM
produced by Hamamatsu (multi pixel photon counter, MPPC):
100 (S10362-11-100U) and 400 pixels
(S10362-11-050U), both with 1$\times 1$ mm$^2$ active area.
To manage the signals, the electronic
service of the Frascati Laboratory (SELF) has
developed a custom electronics composed by
a 1$\times$2 cm$^2$ card, containing the pre-amplifier, and
a multifunction NIM board.
For these tests, we have set the
pre-amplifier gain to 50. The NIM board
supplies the voltage to the photodetector (Vbias) with
a precision of 2 mV and a stability at the level of
0.03 permill. A low threshold
discriminator and a fanout are also present in the board.

To determine the gain, we have prepared a setup
based on a blue light LED and a polaroid filter to change
light intensity.
We have measured the gain and the dark rate
variation as a function of the applied Vbias
and the temperature of the photodetector.

The readout electronics was based on CAMAC,
with a charge sensitivity of 0.25 pC/count
and a time of 125 ps/count.

Our tests confirm the performances declared by Hamamatsu and
show a significative variation of the detector gain
as a function of the temperature.
The 400 pixels shows a temperature dependence of the gain which is
a factor four smaller than the 100 pixels (3\% versus
12\%), with a gain reduction of
a factor five.

\subsection{Tests on fibers}

We have studied the light response of three different,
1 mm$^2$, fibers optically connected to MPPC:
\begin{itemize}
\item Kuraray SCSF81 (blue scintillating)
\item Saint Gobain BCF92 single cladding
(WLS from blue to green)
\item Saint Gobain BCF92 multi cladding
(WLS from blue to green)
\end{itemize}

The test is done firing the fiber with a Sr$^{90}$ source.
The trigger is provided by a NE110 scintillator finger
(1$\times$5$\times$0.5
cm$^3$) connected to a bialkali photomultiplier positioned below
the fiber.

\begin{figure}[htb]
\psfig{file=sgmc.eps,width=7cm}
\vspace{-1cm}
\caption{Charge response to Sr$^{90}$ for a BCF92
multicladding fiber.The different
spectra are obtained setting threshold at 0.5, 1.5 and 2.5
photoelectrons.}
\label{figsgmc}
\end{figure}

As expected, a large light yield is shown for SCSF81 while the WLS fibers
have a reduced response. However, the BCF92 multi cladding has a
reasonable light yield as shown in Fig.\ref{figsgmc}.
For this fiber we have: maximum
light yield, fast response (5 ns/pe) and high
attenuation length (3.5 m).

\subsection{Tests on tiles}

Light response and time resolution of tiles have been
measured using cosmic rays. The system was prepared connecting
the fiber to the MPPC and using two NE110 fingers to trigger
the signal. We have prepared two different tiles:
\begin{description}
\item[A]  3 mm thick tile with 400 pixels MPPC,
\item[B]  5 mm thick tile with 100 pixels MPPC.
\end{description}
Using ADC distribution we find: 14 pe/mip for tile ``A'' and
26 pe/mip for tile ``B''
(See Fig.\ref{figtile}). These results
are comparable taking into account the thickness ratio between tiles
and the photon detection efficiency of the two
detectors (40\% for 400 pixels and 45\% for 100 pixels).

\begin{figure}[htb]
\vspace{9pt}
\psfig{file=tile.eps,width=7cm}
\vspace{-1cm}
\caption{ADC distribution, in number of photoelectrons,
for two different type of tiles.}
\label{figtile}
\end{figure}

Correcting for the time dependence
on pulse height, we find a preliminary
time resolution of 1000 ps (750 ps) for tile ``A'' (``B'').

\subsection{Next plans}

We are now assembling two small
dimension multi-tiles
prototypes of the QCAL, to study
signal transportation and to
measure the effective radiation length.
In 2009,
we plan also to construct a ``module 0''
consisting of a complete slice of the dodecagon
(1/12 of one calorimeter) with final material and electronics.

\section{Crystal calorimeter with timing, CCALT}

In the new design of Da$\Phi$ne
interaction region, the
position of the quadrupoles increases the acceptance
of the central calorimeter from 21$^{\circ}$ to
18$^{\circ}$. Below this limit we can safely insert
few crystals to improve the acceptance for photons
coming from $\eta$ and $K_S$ decays. This detector could
work as veto detector for photons
down to 8$^{\circ}$. The particular region is visible in
Fig.\ref{figkloe} and is delimited between the
focalizing quadrupoles and the spherical interaction region
of the KLOE detector.

The proposed solution is to insert an
homogeneus calorimeter based on LYSO
($Lu_{18}Y_{.2}SiO_5:Ce$) crystals.
The most important characteristics
of these crystals are a very high light yield,
a time emission, $\tau$, of 40 ns,
high density and X$_0$ without beeing
hygroscopic (See Tab.\ref{tablfs}).
These crystals well match the request
of high efficiency to low energy
photons and excellent time resolution for
the CCALT, which will help in
fighting the high level of machine
background events present in the low energy
region.

The preliminary project consists in a dodecagonal
barrel for each side of the interaction region,
composed by LYSO crystals (2$\times$2$\times$13 cm$^3$)
readout using
avalanche photodiodes (APD).

\subsection{Preliminary tests on LYSO}

We have tested one LYSO crystals from Saint Gobain
(2$\times$2$\times$15 cm$^3$)
using both
cosmic rays and electrons from Frascati beam
test facility (BTF).
For cosmic rays test, the crystal was
readout using standard bialkali photomultiplier,
while in the test beam
we used a 5$\times$5 mm$^2$ avalanche
photodiode provided
by Hamamatsu (S8664-55).

In both cases, the crystal was simply wrapped with mylar and
the photodevice coupled with optical grease and
simple mechanical arrangement.

After correcting the time dependence on
pulse height,
we obtain, from cosmic rays, a preliminary time resolution
$\sigma_T=360$ ps,
which corresponds to 12000 pe/mip
($N_{pe}=(\tau/\sigma_T)^2$).
Assuming a MIP to deposit 10 MeV/cm
this correponds to 600 pe/MeV.

At the BTF, we have used
electrons from 50 to 400 MeV to
measure the time resolution
of the crystals which is well
parametrized by:
$$
\sigma_T=\frac{82\,ps}{\sqrt{E(GeV)}}\oplus 293\,ps
$$

The constant term is probably related to the different
arrival times of the showering photons along
the crystal axis.
From the statistical term we derive a light yield
of 240 pe/MeV which is in rough agreement with
the results from cosmic rays taking into account
quantum and collection efficiency.

\subsection{Next plans}

A dedicated test beam
with electrons will be held in Frascati beam test
facility in the first
months of 2009.
The aim of the test is
to characterize energy and time
resolution of a 3$\times$3 matrix of
high quality crystals surrounded by
an outer leakage detector
done with PbWO.
This test will also allow to compare the response
of different kind of LYSO crystals (Saint Gobain
and Scionix) and a similar crystal (LFS by
Zecotek) which could be an alternative candidate
(See Tab.\ref{tablfs}).

\begin{table}
\begin{tabular}{|c|c|c|}
\hline
 & LYSO & LFS \\
\hline
Density & 7.1 & 7.2-7.3 \\
Attenuation lenght (cm) & 1.2 & 1.12 \\
Decay constant (ns) & 41 & 35-36 \\
Max emission (nm) & 420 & 435-438 \\
Light yield (relative NaI) & 75 & 80-85 \\
Energy resolution & 8 & 9-12 \\
Hygroscopic & no & no \\
Refractive index & 1.81 & 1.78 \\
\hline
\end{tabular}
\caption{Comparison between LYSO and LFS characteristics}
\label{tablfs}
\end{table}

\section{Conclusions}

The new scheme proposed for the Da$\Phi$ne machine allows
a factor 5 increase in the delivered luminosity.
Some R\&D are in progress to add new
components to the KLOE apparatus. We have presented the
proposal for two low angle calorimeters. QCALT is a tile
calorimeter surronding the focalizing
quadrupoles to increase the coverage of the
electromagnetic calorimeter. CCALT is a LYSO calorimeter
aimed to be a good veto for low angle photons.
We have presented the preliminary measurement on tiles and
crystals and the characterization of the
photodevices used with tiles.

\section{Acknowledgement}

We want to thank G.Corradi, D.Tagnani and C.Paglia to
have developed the readout electronics of the SiPM; M.Cordelli
for his help during the preparation of test setup;
M.Arpaia, G.Bisogni, A.Cassar\`a, A.Di Virgilio, U.Martini
and A.Olivieri for their help in the mechanical preparation
of the setup and in the
preparation of the tiles.

Instrumentation and Detectors (physics.ins-det); High Energy Physics - Experiment (hep-ex)


Does self-replication imply evolvability?

Abstract: The most prominent property of life on Earth is its ability to evolve. It is often taken for granted that self-replication--the characteristic that makes life possible--implies evolvability, but many examples such as the lack of evolvability in computer viruses seem to challenge this view. Is evolvability itself a property that needs to evolve, or is it automatically present within any chemistry that supports sequences that can evolve in principle? Here, we study evolvability in the digital life system Avida, where self-replicating sequences written by hand are used to seed evolutionary experiments. We use 170 self-replicators that we found in a search through 3 billion randomly generated sequences (at three different sequence lengths) to study the evolvability of generic rather than hand-designed self-replicators. We find that most can evolve but some are evolutionarily sterile. From this limited data set we are led to conclude that evolvability is a likely--but not a guaranteed-- property of random replicators in a digital chemistry.

\section{Introduction}
For life of the type as we experience it on Earth to emerge from an initially abiotic state requires two seemingly independent things to happen. First, a self-replicator has to emerge (or else, to arrive on Earth from extraterrestrial sources~\citep{Arrhenius1908,HoyleWickramasinghe1981,Wickramasinghe2011}). Secondly, this self-replicator must be able to evolve and thus diversify into the complexity we see today. In general we think of self-replicators as extremely rare~\citep{vonNeumann1966}. But in order to jump-start the process of evolution, these rare self-replicators now also must be endowed with another property--evolvability~\citep{WagnerAltenberg1996,Kirschner1998,Wagner2005}. By evolvability here we mean the ability to produce variants or mutants (produced by a faulty self replication process or else by external noise) that can also self-replicate. Without evolvability the self-replicator would just multiply but not adapt. While it is safe to assume that one of the first self-replicators on Earth was evolvable ($n=1$), it is not at all clear whether evolvability is a general property of self-replicators, or else is an additional constraint that renders the emergence of life even more improbable. Here we use the computational evolution system Avida~\citep{AdamiBrown1994,Adami1998,Ofriaetal2009} to investigate whether random self-replicators, that is, randomly generated sequences of code written in the avidian instruction set that happen to be able to self-replicate, also are automatically endowed with the capacity to evolve. Even though evolvability appears to be inherent to avidians, we cannot rule out a priori whether the observed evolvability of avidians is a consequence of the evolvabilty of the hand-written ancestor or is instead a germane property of all self-replicators that can exist within this digital chemistry.  To resolve this question, we searched for self-replicating sequences that we generated randomly, within the three size classes $L=8$, $L=15$ (the size of the standard avidian self-replicator), and $L=30$. We found 170 such sequences after generating 3 billion random sequences, and report their evolvabilty below.

\section{Self-replicators in Avida}
In Avida, small self-replicating computer programs (``avidians") compete for limited memory space and limited CPU resources needed to successfully self-replicate (see \citealt{Ofriaetal2009} for a complete description of the Avida system). This ability to self-replicate is contained within an individual avidian's genome of instructions (see Fig.~\ref{avida}). Because this genome is then passed on to an avidian's descendants, the ability to self-replicate is heritable. Selection in Avida is implicit, since faster replicators have a higher chance of making offspring.  During the process of self-replication, mutations may be introduced, resulting in error-prone replication and variation within the population. Thus, Avida is an {\em instance} of evolution by natural selection~\citep{Pennock2007}, because it contains inheritance, variation, and differential fitness among individuals. Avida has been used to explore a diverse set of topics in evolutionary biology such as the evolution of genomic complexity~\citep{lenski1999genome}, the ``survival of the flattest" effect at high mutation rates~\citep{wilke2001evolution}, the evolution of complex features~\citep{lenski2003evolutionary}, the evolution of reproductive division of labor~\citep{goldsby2014evolutionary}, and the role of coevolution in the origin of complexity~\citep{zaman2014coevolution}.
\begin{figure}[!htbp]
   \centering
   \includegraphics[width=3.5in]{Fig1.pdf}
   \caption{The avidian CPU in the process of executing a segment of code. The CPU uses three registers (AX,BX,CX) as well as an instruction pointer (IP) that reads the program into the CPU.  A read-head, a write-head, and a flow-head are used to specify positions in the CPU's memory. The `copy' command reads from the read-head and writes to the write-head, while `jump'-type statements move the instruction pointer to the flow-head. The CPU uses two stacks to simulate an ``infinite Turing tape", while input/output buffers serve to communicate between the CPU and its environment (reproduced from~\cite{Ofriaetal2009}, with permission). }
   \label{avida}
\end{figure}

Avidian genomes are composed of usually 26 different instructions, typically rendered as a string of lowercase letters, where each letter corresponds to one command (see Table 1). These instructions can be considered analogous to the 20 amino acids common to all biological organisms.
\begin{table*}[htbp]
   \centering

   \begin{tabular}{@{} lll @{}}
      \toprule
      Instruction    & Description & Symbol\\
      \midrule
nop-A    & no operation (type A) & a \\
nop-B   & no operation (type B) & b \\
nop-C   & no operation (type C) & c \\
if-n-equ & Execute next instruction only-if ?BX? does not equal complement & d\\
if-less & Execute next instruction only if ?BX? is less than its complement & e\\
if-label & Execute next instruction only if template complement was just copied & f\\
mov-head & Move instruction pointer to same position as flow-head & g\\
jmp-head & Move instruction pointer by fixed amount found in register CX & h\\
get-head & Write position of instruction pointer into register CX & i\\
set-flow & Move the flow-head to the memory position specified by ?CX? & j\\
shift-r & Shift all the bits in ?BX? one to the right & k\\
shift-l & Shift all the bits in ?BX? one to the left & l\\
inc & Increment ?BX? & m\\
dec & Decrement ?BX? & n\\
push &       Copy value of ?BX? onto top of  current stack & o\\
pop & Remove number from current stack and place in ?BX? & p\\
swap-stk & Toggle the active stack & q\\
swap & Swap the contents of ?BX? with its complement & r\\
add & Calculate  sum of BX and CX; put  result in ?BX? & s\\
sub & Calculate  BX minus CX; put result in ?BX? & t\\
nand & Perform bitwise NAND on BX and CX; put  result in ?BX? & u\\
h-copy &  Copy instruction from read-head to write-head and advance both & v\\
h-alloc & Allocate memory for offspring & w\\
h-divide & Divide off an offspring located between read-head and write-head & x \\
IO &  Output value ?BX? and replace with new input & y\\
h-search & Find complement template and place flow-head after it & z\\
      \bottomrule
   \end{tabular}
   \caption{Instruction set of the avidian programming language used in this study. The notation ?BX? implies that the command operates on a register specified by the subsequent nop instruction (for example, nop-A specifies the AX register, and so forth). If no nop instruction follows, use the register BX as a default. More details about this instruction set can be found in~\cite{Ofriaetal2009}.}
   \label{avidainst}
\end{table*}

A self-replicator in Avida must have a set of instructions in the right order to: first allocate memory, then copy itself into this newly allocated memory, and finally it has to split the new memory from the old one in order to crate a new organism. These instructions and their interactions are complex and in general it is not possible to predict if an organism will replicate just by examining the sequence. As a consequence, the ability to self-replicate has to be tested directly by allowing the instructions to execute (we surmise that this problem of predicting the ability to self-replicate is similar to the Halting Problem~\citep{turing1936computable}, where it is necessary to run the algorithm to find out if it halts or not).

In the majority of all Avida experiments a default hand-designed ``start" organism is used to seed a population. This choice is to some extent historical: the designers assumed that self-replicators were too rare to be found by a random process~\citep{Adami1998}. Indeed, any particular sequence of length $L$=15 for example (the length of one of the standard hand-written ancestors) has a likelihood of $p=26^{-15}\approx 6\times 10^{-22}$. If a million of such sequences could be checked per second, on one thousand CPUs running in parallel, it would take about 50,000 years to find it. As a remedy, the designers wrote one instead (most of the common ancestors in Avida were written by Charles Ofria). However, it is clear that the density of self-replicators in the space of all sequences depends on the chemistry (here, the instruction set) used. Since the inception of Avida in 1993, the standard instruction set has changed, and it appears that using the current set (that is, the current ``chemistry") self-replicators can be found among randomly generated sequences (as reported here) because the information content of the sequences is significantly smaller than 15.

No self-replicator in nature (or within Avida) replicates without error, because noise in the system is inevitable and will always affect replication by increasing variation. Variation in Avida can have two different causes. First, it  is possible that the sequence of commands performs in such ways that the resulting copy is modified. This typically leads to repeated commands, insertions, or deletions. Such changes are inherent to the self-replicator (they are genetically encoded and thus deterministic) and referred to as ``implicit" mutations. A second mechanism is less deterministic, and occurs during the processes of division and reproduction, processes that are inherently made to be error-prone to a degree that can be specified by the user. These copy-errors, as well as the probabilistic insertion and deletion of instructions or code snippets reflect the noisiness of the system and are not under control of the organism itself. Note that even the most sophisticated polymerases that also have proof reading abilities cannot create perfect replicates every time, while this is in principle possible for avidians without implicit mutations by turning the mutation rate to zero. While biochemical-based life does not have a mutational mechanism similar to the implicit mutations seen in avidians, it is not unreasonable to assume that early replicators also underwent large mutations due to a lack of error-correction during replication.

\section{Evolvability in Avida}
Evolvability describes the ability of an organism to undergo evolutionary adaptation~\citep{WagnerAltenberg1996,Kirschner1998,Wagner2005}. There are two main pathways to adaptation in Avida. The first is {\em optimization}, where organisms with a minimal replication time outcompete others in the population (akin to $r$-selection~\citep{Pianka1970}). The second is {\em innovation}, where avidians can evolve new phenotypic traits that enable a fitness increase (the $K$-selection mode). These phenotypic traits are the ability to perform certain Boolean logic operations; these logic operations allow an individual to execute a greater proportion of its genome than other avidians that do not perform such logic operations. This indirectly allows them to self-replicate at a greater rate (see, for example,~\citealt{Adami1998,Adami2006a}). These different modes of survival can be thought of as different niches that avidians can inhabit~\citep{WhiteAdami2004}.

It is possible in general that an organism (digital or otherwise) carries such a specific sequence of nucleotides or (as in Avida) commands that every possible mutation prevents self-replication. From a fitness landscape point of view, such replicators would represent an isolated fitness peak where self replication is only possible at the top, and where each mutation is lethal. Such self-replicators, born at the top of the fitness peak, so to speak, would be unevolvable. A self-replicator that can evolve, on the other hand,  would never emerge at peak in the landscape, but rather on a lower level and subsequently find positive and or neutral mutations that ultimately lead to higher fitness levels in the landscape. We know that the hand-written replicators in Avida are evolvable, but it is not clear how likely evolvability is among randomly generated replicators.

\section{Methods}
All experiments are performed using Avida version 2.14, which can be downloaded from https://github.com/devosoft/avida. We first randomly generated Avida sequences to discover individuals that could self-replicate, using a uniform random distribution for each command at each site, which ensures that each sequence has the same likelihood to be generated, given by
\begin{eqnarray}{}
p=\frac{1}{D^{L}}\;,
\end{eqnarray}{}
where $D$ is the size of the alphabet ($D$=26 here) and $L$ is the length of the sequence generated~\citep{AdamiLaBar2015}. To decide whether a sequence could successfully self-replicate, it must pass two tests. First, we tested whether the organism could successfully divide within its lifespan. Here, we used standard  Avida parameters for an organism's lifespan: it must divide before it executes $20\times L$ instructions. This test indicates that an avidian can successfully reproduce, it does not imply that its descendants also can reproduce. In our search we discovered many viable avidians that were able to successfully divide into two non-viable organisms. Therefore, we only counted sequences that could replicate and produce offspring that could also replicate as true self-replicators (in other words, they are ``colony-forming"). This does not imply that every self-replicator produces a perfect copy of itself in the absence of mutation. Indeed, most of these replicators undergo implicit mutations solely due to their genome sequence, and their offspring differ in length from the parent. In analyzing a genome's ability to self-replicate, we used the default Avida settings, described for example in~\citep{Ofriaetal2009}.
For our study of evolvability, we also included the default length 15 hand-written ancestor in our set of self-replicators.

In order to test whether these self-replicators could optimize their fitness, we evolved them in an environment where {\em only} decreased self-replication speed was under positive selection (the $r$-selection niche). We evolved these replicators for $10^3$ generations at a population size of $3,600$ individuals (10 replicates). Instruction mutations occurred at a genomic rate of 0.1 at division (meaning the likelihood to incur an error is length-independent), and both insertions and deletions occurred at a genomic rate of 0.005 per division. We measured an organism's evolvability as its gain in relative fitness at the end of the experiment.

Finally, we tested the self-replicators' ability to evolve greater complexity by evolving new phenotypic traits (that is, in the $K$-selection mode). The traits that are rewarded are logical functions that an avidian can solve by stringing together code involving the {\tt nand} function (given by the letter {\tt u}, see Table 1). In this setting, nine different logic tasks may be rewarded, denoted as the bit-wise {\tt NOT, NAND, AND, OR-NOT,OR, AND-NOT, NOR, XOR}, and {\tt EQU}. We evolved the self-replicators of length 15 in an environment where nine Boolean logic operations, and hence nine phenotypic traits, are under positive selection; this is often referred to as the logic-9 environment~\citep{lenski2003evolutionary}. For this experiment, we evolved these replicators for $10^4$ generations at a population size of $10^4$ individual (ten replicates). Increased population size and experiment time was used to allow time for trait evolution. The mutation rates were the same as in the optimization experiments. We quantified an organism's evolvability in these experiments by measuring the number of evolved phenotypic traits.

\section{Results}
Of the $10^9$ randomly-generated Avidian genomes in each length class, we found 6 self-replicators of length 8, 58 self-replicators of length 15, and 106 self-replicators of length 30 (see Table~\ref{tab:seq} for their sequences). However, it is unlikely that each self-replicator needs every instruction in its genome in order to self-replicate~\citep{Adami2015}. Many of these self-replicators could be unstable and would generate implicit mutations upon replication. Therefore, we tested the replication fidelity of each self-replicator (Fig.~\ref{fidelityAll}). All replicators of length 8 were able to undergo error-free replication without external noise (mutation rate set to zero). This is likely due to the fact that avidian genomes cannot be shorter than length 8; however, that does not prevent the occurrence of genomes that replicate to a length larger than 8, which we did not see among this set. Self-replicators of length 15 and 30 could not replicate error-free, even without external noise, and most of their genome sequences decreased upon successful replication.

\begin{figure}[htb]
\begin{center}
\includegraphics[width=8cm]{fidelityAll.pdf}
\caption{Change in length after the first round of replication for all self replicators we found of length 8, 15, and 30. The self-replicators are ordered by length increase within their individual groups. The dashed lines indicate the start organism's length.}
\label{fidelityAll}
\end{center}
\end{figure}

Next, we tested the evolvability of these replicators in the sense of how well they could  optimize their genomes as a way to increase replication speed, and thus their fitness. We found that the majority of these spontaneous self-replicators of length 8, 15, and 30 possess the ability to optimize their replication algorithm (see Figure~\ref{relativeFitnessAll}). While the increase in fitness across the different self-replicators varied due to differing ancestral fitness, most increased in fitness by more than two-fold over the course of $10^3$ generations of evolution. However, we found a few replicators (~3\%) that were evolutionary sterile (defined as a relative fitness $<2$), demonstrating that not all self-replicators can easily undergo significant adaptation.

Because the ability to optimize replication speed does not automatically imply the ability to evolve greater complexity, we also tested at the ability of the length 15 replicators to evolve new phenotypic traits (see Methods for a description of those traits). All 58 self-replicators were able to evolve at least one phenotypic trait in one replicate, although the likelihood of the emergence of phenotypic traits varied greatly across the different replicators and within the replicates for each specific replicator (see Fig.~\ref{traitsLength15}). Moreover, most self-replicators were able to evolve multiple traits. The hand-written Avida ancestor of length 15 displayed the {\em least} evolvability in regards to trait evolution compared to the 58 randomly-generated self-replicators. Only one out of ten replicates evolved any phenotypic traits in the allotted time, although that replicate did manage the evolution of two traits.

\begin{figure}[htb]
\begin{center}
\includegraphics{relativeFitnessAll.pdf}
\caption{Relative fitness after 1,000 generatios of evolution for all replicators we found of length 8, 15, and 30. The replicators are ordered by relative fitness increase within their groups of similar length. The fitness of all start organisms is normalized to $1$, indicated by the dashed line. The error bars indicate the standard deviation over $10$ replicate evolution trials with the same ancestor.}\label{relativeFitnessAll}
\end{center}
\end{figure}

\begin{figure}[htb]
\begin{center}
\includegraphics[width=3.5in]{phenotypic_traits_revised.pdf}
\caption{The number of logic tasks that spontaneous replicators of length 15 were able to perform after $10,000$ generations of evolution in the Avida logic-9 environment. Replicators are ordered by number of tasks achieved. Error bars indicate the standard error of the mean for 10 replicates per ancestral organism.}
\label{traitsLength15}
\end{center}
\end{figure}

\section{Discussion}
Here, we asked if spontaneous self-replicators in Avida would all possess the additional ability to change in such a way that they could initiate evolution. We showed that the majority of self-replicators are robust enough to tolerate mutations as well as changes to their genome size. This suggests that self-replicators in this digital chemistry are robust and thus will most likely be able to jump-start evolution. Further, this result holds both when we looked at the ability of self-replicators to increase fitness through optimization or through the evolution of new phenotypes. In addition, many spontaneous self-replicators made faulty copies of themselves in the absence of external noise, which further supports the idea that variability is an inherent property of spontaneous self-replication.

These results confirm that not every spontaneous self-replicator necessarily is evolvable. On the other hand, we also find that the majority of spontaneous self-replicators deterministically do not make exact copies even in the absence of mutation, questioning the term ``self-replicator" for these genotypes. Indeed, while these spontaneously generated sequences are ``colony-forming" (in the sense that they produce sequences that can also replicate themselves), these are not self-replicators in the strictest sense of the definition. Instead, these replicators mutate themselves to become a self-replicator: we call these individuals ``proto-self-replicators" or ``proto-replicators". This means that the number of potential sequences that can start evolution is increased by existence of these proto-replicators. In natural chemistries, this phenomenon has been discussed widely~\citep{bernal1949physical,miller1953production,eigen1971,miller1974origins,eschenmoser1992chemistry}, suggesting that the first self-replicator might have needed either self organisation or some other forms of catalysis or autocatalysis in order to create the right chemistry to start polymerization in the first place. Figure \ref{overview} gives an illustration of this concept.

Leo Tolstoy famously begun his ``Anna Karenina" by remarking that {\em``Happy families are all alike; every unhappy family is unhappy in its own way."}. In the same manner we asked here whether, with respect to evolvability, all self-replicators are alike (while of course all non-replicators are non-replicating in their own way). Our results show that not every self-replicator is suitable as the ancestor for evolution. But we also find that, with respect to the capacity to evolve functional complexity, random replicators are better than even those designed by thinking humans.

Our results further suggest that life, if possibly found elsewhere, does not necessarily experience evolution. We might at some point in the future find self-replicators that maintained their ancestral form due to their inability to tolerate enough variation to evolve. Similarly, it seems possible that, if such an mutationally-fragile self-replicator exists, it might be outcompeted by other self replicators that can take advantage of mutations, and thus attain higher fitness and ultimately outcompete those  mutationally-fragile self-replicators. Thus, as long as evolvability is evolvable, complexity should ensue even when sparked by the most humble and awkward replicators.

\begin{figure}[htb]
\begin{center}
\includegraphics[width=3.5in]{overview.pdf}
\caption{Illustration how an Avida proto-replicator (left) ``copies" itself into a self-replicator (middle). This replicator can now experience external mutations, which leads to a new self replicator. The repetition of such process is what jump-starts evolution (right). The text string is a program example from Avida.}
\label{overview}
\end{center}
\end{figure}

\footnotesize

Populations and Evolution (q-bio.PE); Adaptation and Self-Organizing Systems (nlin.AO); Biomolecules (q-bio.BM)


Non-WIMP Candidates

Abstract: Non-WIMP dark matter candidates include particles motivated by minimality, candidates motivated by experimental anomalies, and exotic possibilities motivated primarily by the desire of clever iconoclasts to highlight how truly ignorant we are about the nature of dark matter. In this review, I discuss candidates that are not WIMPs, but nevertheless share the same theoretical motivations as WIMPs and also naturally have the correct relic density. There are two classes: superWIMP dark matter, where the desired relic density is inherited through decays, and WIMPless dark matter, where the dark matter's mass and couplings scale together to maintain the desired thermal relic density.

\part{Candidates}

\author[J.~L.~Feng]{Jonathan L.~Feng \\ Department of Physics and
  Astronomy, University of California, Irvine, CA 92697, USA}

\setcounter{chapter}{9}
\chapter[Non-WIMP Candidates]{Non-WIMP Candidates$^\dagger$}

\section{Motivations}

\footnote[0]{\hspace*{-.07in}$^\dagger$Published as Chapter 10,
pp.~190-204, in {\em Particle Dark Matter: Observations, Models and
Searches}, edited by Gianfranco Bertone (Cambridge University Press,
2010), available at
http://cambridge.org/us/catalogue/catalogue.asp?isbn=9780521763684.}There
are many non-WIMP dark matter candidates.  Two prominent and highly
motivated examples are axions and sterile neutrinos, which are
reviewed in Chapters~11\rem{\ref{Chap:Sikivie}} and
12\rem{\ref{Chap:Shaposhnikov}}, respectively.  In addition, there are
candidates motivated by minimality, particles motivated by
experimental anomalies, and exotic possibilities motivated primarily
by the desire of their inventors to highlight how truly ignorant we
are about the nature of dark matter.

In this brief Chapter, we focus on dark matter candidates that are not
WIMPs, but which nevertheless share the most important virtues of
WIMPs.  As discussed in Chapters~7\rem{\ref{Chap:Gelmini}},
8\rem{\ref{Chap:Servant}}, and 9\rem{\ref{Chap:Ellis}}, WIMPs have
several nice properties:
\begin{itemize}
\item They exist in well-motivated particle theories.
\item They are naturally produced with the correct thermal relic
  density (the ``WIMP miracle'').
\item They predict signals that may be seen in current and near future
  experiments.
\end{itemize}
The candidates we discuss also have all three of these properties.
They fall naturally into two classes: superWIMP candidates, which
inherit the correct relic density through decays, and WIMPless
candidates, which have neither weak-scale masses nor weak
interactions, but which nevertheless have the correct thermal relic
density.  These possibilities appear in the same particle physics
frameworks as WIMPs, but they imply very different cosmological
histories for our Universe, as well as qualitatively new dark matter
signals for both astrophysical observatories and particle physics
experiments.

\section{SuperWIMP Dark Matter}

\subsection{Candidates and Relic Densities}

In the superWIMP framework for dark matter, WIMPs freeze out as usual
in the early Universe, but later decay to superWIMPs,
superweakly-interacting massive particles that form the dark matter
that exists today.  Because superWIMPs are very weakly-interacting,
they have no impact on WIMP freeze out in the early universe, and the
WIMPs decouple, as usual, with a thermal relic density
$\Omega_{\text{WIMP}}$ that is naturally near the required density
$\Omega_{\text{DM}} \approx 0.23$.  Assuming that each WIMP decay
produces one superWIMP, the relic density of superWIMPs is
\begin{equation}
\Omega_{\text{SWIMP}} = \frac{m_{\text{SWIMP}}}{m_{\text{WIMP}}}
\Omega_{\text{WIMP}} \ .
\end{equation}
SuperWIMPs therefore inherit their relic density from WIMPs, and for
$m_{\text{SWIMP}} \sim m_{\text{WIMP}}$, they are also naturally
produced in the desired amount to be much or all of dark matter.  The
evolution of number densities is shown in \figref{freezeout_swimp}.

\begin{figure*}[tb]
\centering
\includegraphics[height=10cm,angle=-90]{freezeout_swimp_rot.eps}
\caption{In superWIMP scenarios, WIMPs freeze out as usual, but then
decay to superWIMPs, superweakly-interacting massive particles that
form dark matter.}
\label{fig:freezeout_swimp}
\end{figure*}

The superWIMP scenario is realized in many particle physics models.
The prototypical example is gravitinos, which exist in all
supersymmetric theories~\cite{Feng:2003xh,Feng:2003uy,%
Feng:2003nr,Ellis:2003dn,Buchmuller:2004rq,Wang:2004ib,Feng:2004zu,%
Feng:2004mt,Feng:2004gn,Ellis:2004bx,Roszkowski:2004jd}.  In the
simplest supersymmetric models, supersymmetry is transmitted to
standard model superpartners through gravitational interactions, and
supersymmetry is broken at a high scale.  The mass of the gravitino
$\tilde{G}$ is
\begin{equation}
m_{\tilde{G}} = \frac{F}{\sqrt{3} M_{*}} \ ,
\label{gravitinomass}
\end{equation}
and the masses of standard model superpartners are
\begin{equation}
\tilde{m} \sim \frac{F}{M_{*}} \ ,
\label{tildem}
\end{equation}
where $M_{*} = (8 \pi G_N)^{-1/2} \simeq 2.4 \times 10^{18}~\,\mbox{GeV}$ is
the reduced Planck scale and $F \sim (10^{11}~\,\mbox{GeV})^2$ is the
supersymmetry-breaking scale squared.  The precise ordering of masses
depends on unknown, presumably ${\cal O}(1)$, constants in
\eqref{tildem}. It is, then, perfectly possible that the gravitino is
the lightest supersymmetric particle (LSP) and a candidate for
superWIMP dark matter. The role of the decaying WIMP is played by the
next-to-lightest supersymmetric particle (NLSP), typically a slepton,
sneutrino, or neutralino.  As required, the gravitino couples very
weakly, with interactions suppressed by $M_{*}$, and so it is
irrelevant during the WIMP's thermal freeze out.

The gravitino superWIMP scenario differs markedly from other gravitino
dark matter scenarios.  In previous
frameworks~\cite{Pagels:1981ke,Weinberg:1982zq,Krauss:1983ik,%
Nanopoulos:1983up,Khlopov:1984pf,Ellis:1984eq,Ellis:1984er,%
Juszkiewicz:1985gg,Ellis:1990nb,Moroi:1993mb,Bolz:2000fu}, gravitinos
were expected to be produced either thermally, with
$\Omega_{\tilde{G}} \sim 0.1$ obtained by requiring $m_{\tilde{G}} \sim
\,\mbox{keV}$, or through reheating, with $\Omega_{\tilde{G}} \sim 0.1$
obtained by tuning the reheat temperature to $T_{\text{RH}} \sim
10^{10}~\,\mbox{GeV}$.  In the superWIMP scenario, the desired amount of dark
matter is obtained without relying on the introduction of new,
fine-tuned energy scales.

Other examples of superWIMPs include Kaluza-Klein gravitons in
scenarios with universal extra
dimensions~\cite{Feng:2003xh,Feng:2003uy,Feng:2003nr},
axinos~\cite{Rajagopal:1990yx,Covi:1999ty,Covi:2001nw} and
quintessinos~\cite{Bi:2003qa,Bi:2004ys} in supersymmetric theories,
and many other scenarios in which a metastable particle decays to the
true dark matter particle through highly suppressed interactions, with
lifetimes ranging from fractions of a second to beyond the age of the
Universe.

\subsection{Astrophysical Signals}

Because superWIMPs are very weakly interacting, they are impossible to
detect in conventional direct and indirect dark matter search
experiments.  At the same time, the extraordinarily weak couplings of
superWIMPs imply that the decays of WIMPs to superWIMPs may be very
late and have an observable impact on, for example, Big Bang
nucleosynthesis (BBN), the Planckian spectrum of the cosmic microwave
background (CMB), small scale structure, the diffuse photon flux, and
cosmic ray experiments.

In the prototypical case of a slepton decaying to a gravitino
superWIMP, the decay width is
\begin{equation}
 \Gamma(\tilde{l} \to l \tilde{G}) =\frac{1}{48\pi M_{*}^2}
 \frac{m_{\tilde{l}}^5}{m_{\tilde{G}}^2}
 \left[1 -\frac{m_{\tilde{G}}^2}{m_{\tilde{l}}^2} \right]^4 \ ,
\label{sfermionwidth}
\end{equation}
assuming the lepton mass is negligible. This decay width depends on
only the slepton mass, the gravitino mass, and the Planck mass.  For
$m_{\tilde{G}} / m_{\tilde{l}} \approx 1$, the slepton decay lifetime is
\begin{eqnarray}
 \tau(\tilde{l} \to l \tilde{G})
\simeq 3.6\times 10^8~\text{s}
\left[\frac{100~\,\mbox{GeV}}{m_{\tilde{l}} - m_{\tilde{G}}}\right]^4
\left[\frac{m_{\tilde{G}}}{\,\mbox{TeV}}\right]\ .
\label{decaylifetime}
\end{eqnarray}
This expression is valid only when the gravitino and slepton are
nearly degenerate, but usefully illustrates that decay lifetimes of
the order of days or months are perfectly natural. Similar expressions
hold for the decay of a neutralino NLSP to a gravitino.

\subsubsection{BBN and CMB}

Signals in BBN and the CMB are determined primarily by the WIMP
lifetime and the energy released in visible decay products when the
WIMP decays.  This energy release destroys and creates light elements,
distorting the predictions of standard BBN.  In addition, the
injection of electromagnetic energy may also distort the frequency
dependence of the CMB away from its ideal black body spectrum.  For
the decay times of interest with redshifts $z \sim 10^5$ to $10^7$,
the resulting photons interact efficiently through $\gamma e^- \to
\gamma e^-$ and $e X \to e X \gamma$, where $X$ is an ion, but photon
number is conserved, since double Compton scattering $\gamma e^- \to
\gamma \gamma e^-$ is inefficient.  The spectrum therefore relaxes to
statistical but not thermodynamic equilibrium, resulting in a
Bose-Einstein distribution function
\begin{equation}
f_{\gamma}(E) = \frac{1}{e^{E/(kT) + \mu} - 1} \ ,
\end{equation}
with chemical potential $\mu \ne 0$.

The energy release is conveniently expressed in terms of
\begin{eqnarray}
\xi_{\text{EM}} \equiv \epsilon_{\text{EM}} B_{\text{EM}}
Y_{\text{NLSP}}
\label{eq:xi_EM}
\end{eqnarray}
for electromagnetic energy, with a similar expression for hadronic
energy.  Here $\epsilon_{\text{EM}}$ is the initial EM energy released
in NLSP decay, and $B_{\text{EM}}$ is the branching fraction of NLSP
decay into EM components.  $Y_{\text{NLSP}} \equiv
n_{\text{NLSP}}/n_{\gamma}$ is the NLSP number density just before
NLSP decay, normalized to the background photon number density
$n_{\gamma} = 2 \zeta(3) T^3 / \pi^2$.  It can be expressed in terms
of the superWIMP abundance:
\begin{equation}
Y_{\text{NLSP}}\simeq 3.0 \times 10^{-12}
\left[\frac{\,\mbox{TeV}}{m_{\tilde{G}}}\right]
\left[\frac{\Omega_{\tilde{G}}}{0.23}\right] \ .
\label{eq:def_Y}
\end{equation}
Once an NLSP candidate is specified, and assuming superWIMPs make up
all of the dark matter, with $\Omega_{\tilde{G}} = \Omega_{\text{DM}}
= 0.23$, both the lifetime and energy release are determined by only
two parameters: $m_{\tilde{G}}$ and $m_{\text{NLSP}}$.  The results for
slepton and neutralino NLSPs are given in \figref{prediction}.

\begin{figure}[tb]
\includegraphics*[width=10cm]{prediction.eps}
\caption{Predicted values of WIMP lifetime $\tau$ and electromagnetic
energy release $\zeta_{\text{EM}} \equiv \varepsilon_{\text{EM}} Y_{\WIMP}$ in the $\tilde{B}$ (left)
and $\tilde{\tau}$ (right) NLSP scenarios for $m_{\SWIMP} = 1~\,\mbox{GeV}$, $10~\,\mbox{GeV}$,
\ldots, $100~\,\mbox{TeV}$ (top to bottom) and $\Delta m \equiv m_{\WIMP} -
m_{\SWIMP} = 1~\,\mbox{TeV}$, $100~\,\mbox{GeV}$, \ldots, $100~\,\mbox{MeV}$ (left to right).
For the $\tilde{\tau}$ NLSP scenario, we assume $\varepsilon_{\text{EM}} = \frac{1}{2}
E_{\tau}$.  {}From Ref.~\cite{Feng:2003uy}.
\label{fig:prediction} }
\end{figure}

In \figref{mu}, these predictions are compared with BBN and CMB
constraints.  The shaded regions are excluded by an analysis of BBN
constraints on EM energy release~\protect\cite{Cyburt:2002uv}.  This
analysis has been strengthened by including hadronic constraints and
updated and refined in many ways in recent years, as described in
Chapter~28\rem{\ref{Chap:Jedamzik}}.  Although the excluded region has
shifted around, the basic features remain: some of the gravitino
superWIMP parameter space is excluded, and some remains.  In addition,
late decays to superWIMPs may in fact improve the current disagreement
of standard BBN predictions with the observed $^7$Li and $^6$Li
abundances~\cite{Cumberbatch:2007me,Bailly:2008yy}.

\begin{figure}[tbp]
\includegraphics*[width=10cm]{mu.eps}
\caption{The grid gives predicted values of WIMP lifetime $\tau$ and
electromagnetic energy release $\zeta_{\text{EM}} \equiv \varepsilon_{\text{EM}} Y_{\WIMP}$ in the
$\tilde{B}$ (left) and $\tilde{\tau}$ (right) WIMP scenarios for $m_{\SWIMP} =
100~\,\mbox{GeV}$, $300~\,\mbox{GeV}$, $500~\,\mbox{GeV}$, $1~\,\mbox{TeV}$, and $3~\,\mbox{TeV}$ (top to
bottom) and $\Delta m \equiv m_{\WIMP} - m_{\SWIMP} = 600~\,\mbox{GeV}$, $400~\,\mbox{GeV}$,
$200~\,\mbox{GeV}$, and $100~\,\mbox{GeV}$ (left to right).  For the $\tilde{\tau}$ WIMP
scenario, we assume $\varepsilon_{\text{EM}} = \frac{1}{2} E_{\tau}$. The shaded
regions are excluded in one analysis of BBN
constraints~\protect\cite{Cyburt:2002uv}; the circle gives a region in
which $^7$Li is reduced to observed levels.  The contours are for
$\mu$, which parameterizes the distortion of the CMB from a Planckian
spectrum. {}From Ref.~\cite{Feng:2003uy}.
\label{fig:mu} }
\end{figure}

Figure \ref{fig:mu} also includes contours of the chemical potential
$\mu$, as determined by updating the analysis of Ref.~\cite{Hu:1993gc}.
The current bound is $\mu < 9\times
10^{-5}$~\cite{Fixsen:1996nj,Eidelman:2004wy}. Although
there are at present no indications of deviations from black body,
current limits are already sensitive to the superWIMP scenario, and
future improvements will further probe superWIMP parameter space.

\subsubsection{Small Scale Structure}

In contrast to WIMPs, superWIMPs are produced with large velocities at
late times.  This has two effects.  First, the velocity dispersion
reduces the phase space density, smoothing out cusps in DM halos.
Second, such particles damp the linear power spectrum, reducing power
on small
scales~\cite{Cembranos:2005us,Kaplinghat:2005sy,Borzumati:2008zz}.

Depending on the particular decay time and decay kinematics,
superWIMPs may be cold or warm.  As seen in \figref{small_scale},
superWIMPs may suppress small scale structure as effectively as a 1
keV sterile neutrino.  Some superWIMP scenarios may therefore been
differentiated from standard cold DM scenarios by studies of halo
profiles, and may even be favored by indications that cold DM predicts
halos that are too cuspy, as discussed in
Chapter~3\rem{\ref{Chap:Kaplinghat}}.

\begin{figure}[tb]
\includegraphics[width=8cm]{small_scale.eps}
\vspace*{-.2in}
\caption{The power spectrum for scenarios in which dark matter is
completely composed of WIMPs (solid), half WIMPs and half superWIMPs
(dashed), and completely composed of superWIMPs (dotted).  For
comparison, the this solid curve is the transfer function for a 1 keV
warm DM model.  {}From Ref.~\cite{Kaplinghat:2005sy}.
\label{fig:small_scale}
}
\end{figure}

\subsection{Astroparticle and Collider Signals}

The possibility of long-lived charged particles in superWIMP scenarios
also has many implications for astroparticle and particle physics
experiments.

\subsubsection{Cosmic Rays}

In superWIMP (and other similar) scenarios, long-lived charged
particles may be produced by cosmic rays, resulting in exotic signals
in cosmic ray and cosmic neutrino
experiments~\cite{Albuquerque:2003mi,Bi:2004ys,%
Albuquerque:2006am,Ahlers:2006pf,Ando:2007ds,Albuquerque:2008zs,%
Canadas:2008ey}.  As an example, ultra-high energy neutrinos may
produce events with two long-lived sleptons through $\nu q \to
\tilde{l} \tilde{q}'$ followed by the decay $\tilde{q}' \to
\tilde{l}$.  The sleptons are metastable and propagate to neutrino
telescopes~\cite{Huang:2006ie}, where they have a typical transverse
separation of hundreds of meters.  They may therefore be detected
above background as events with two upward-going, extremely high
energy charged tracks in experiments such as IceCube.

\subsubsection{Colliders}
\label{sec:colliders}

As evident in \eqref{decaylifetime}, in supersymmetric superWIMP
scenarios, the NLSP decays to the gravitino with lifetimes that may be
of the order of seconds to months.  Such particles are effectively
stable in collider experiments, and this scenario therefore implies
that each supersymmetric event will be characterized not by missing
energy, but by two charged, heavy metastable particles.  This is a
spectacular signal that will be cannot escape notice at the
LHC~\cite{Drees:1990yw,Goity:1993ih,%
Nisati:1997gb,Feng:1997zr,Feng:2005gj}.  In addition, given the
possibility of long lifetimes, it suggests that decays to gravitinos
may be observed by capturing slepton NLSPs and detecting their decays.

The sleptons may be captured in water tanks placed outside collider
detectors~\cite{Feng:2004yi}, in the detectors
themselves~\cite{Hamaguchi:2004df}, or by mining LHC cavern walls for
sleptons~\cite{DeRoeck:2005bw}.  In the first case, shown in
\figref{trap_diagram}, the water tanks may be drained periodically to
underground reservoirs where slepton decays may be observed in quiet
environments.  As many as $10^4$ sleptons per year may be stopped in 1
meter thick water tanks, opening up the possibility of a precise
measurement of slepton lifetime and the first study of a gravitational
process at high energy colliders, along with many other
implications~\cite{Feng:2004gn}.

\begin{figure*}[tb]
\centering
\includegraphics[width=7cm]{trap_diagram.eps}
\caption{Configuration for slepton trapping in gravitino superWIMP
scenarios.  {}From Ref.~\protect\cite{Feng:2004yi}.}
\label{fig:trap_diagram}
\end{figure*}

\section{WIMPless Dark Matter}

\subsection{Candidates and Relic Densities}

Under general conditions, the thermal relic density of a particle $X$
is~\cite{Zeldovich:1965,Chiu:1966kg,Steigman:1979kw,Scherrer:1985zt}
\begin{equation}
\Omega_X \propto {1\over \langle \sigma v \rangle}
\sim \frac{m_X^2}{g_X^4} \ ,
\label{omega}
\end{equation}
where $\langle \sigma v \rangle$ is its thermally-averaged
annihilation cross section, and $m_X$ and $g_X$ are the characteristic
mass scale and coupling entering this cross section.  The last step
follows from dimensional analysis.  The WIMP miracle is the statement
that, for $m_X \sim \mweak \sim 100~\,\mbox{GeV} - 1~\,\mbox{TeV}$ and $g_X \sim
\gweak \simeq 0.65$, $\Omega_X$ is roughly $\Omega_{\text{DM}} \approx
0.23$.

Equation (\ref{omega}) makes clear, however, that the thermal relic
density fixes only one combination of the dark matter's mass and
coupling, and other combinations of $(m_X, g_X)$ can also give the
correct $\Omega_X$.  WIMPless models~\cite{Feng:2008ya} are those in
which the correct thermal relic density is achieved with parameters
$(m_X, g_X) \ne (\mweak, \gweak)$.

Because WIMPless dark matter does not have weak interactions, and
existing constraints effectively exclude electromagnetic and strong
interactions, WIMPless dark matter is necessarily hidden dark matter,
that is, dark matter that has no standard model gauge interactions.
Hidden sectors have a long history, and hidden sector dark matter has
been discussed for decades, beginning with work on mirror matter and
related ideas. For a general discussion and references, see
Chapter~9\rem{\ref{Chap:Servant}}.  Here we note only that, counter to
conventional wisdom, existing constraints place only weak bounds on
hidden sectors.  For example, light degrees of freedom change the
expansion rate of the Universe and thereby impact BBN.  The constraint
from BBN is highly sensitive to the temperature of the hidden sector,
however.  Current bounds from BBN on the number of light and heavy
degrees of freedom are given in \figref{gstarlog}.  For hidden sector
temperatures within a factor of 2 of the observable sector, hundreds
of degrees of freedom, equivalent to several copies of the standard
model or the minimal supersymmetric standard model (MSSM), may be
accommodated.

\begin{figure}[tb]
\begin{center}
\includegraphics*[width=10cm]{gstarlog.eps}
\end{center}
\vspace*{-.29in}
\caption{Bounds from BBN in the $(g^{h\, \text{BBN}}_{\text{light}}, g^{h\, \text{BBN}}_{\text{heavy}})$ plane,
where $g^{h\, \text{BBN}}_{\text{light}}$ and $g^{h\, \text{BBN}}_{\text{heavy}}$ are the hidden degrees of
freedom with masses $m < T^h_{\text{BBN}}$ and $T^h_{\text{BBN}} < m < T^h_{\text{RH}}$,
respectively, for hidden to observable sector reheat temperature ratios
$T^h_{\text{RH}}/T_{\text{RH}} = 0.5$, 0.7, 0.8, 1.0 (from top to bottom).  The regions
above the contours are excluded.  The values of $(g^{h\, \text{BBN}}_{\text{light}},
g^{h\, \text{BBN}}_{\text{heavy}})$ are marked for four example hidden sectors: (A)
1-generation and (B) 3-generation flavor-free versions of the MSSM
with $T^h_{\text{BBN}} < m_X < T^h_{\text{RH}}$, and (C) 1-generation and (D) 3-generation
flavor-free versions of the MSSM with $m_X < T^h_{\text{BBN}}/2$.  {}From
Ref.~\cite{Feng:2008mu}.  }
\label{fig:gstarlog}
\end{figure}

Of course, WIMPless dark matter requires hidden sectors with
additional structure to guarantee that the hidden sector's dark matter
has the desired thermal relic density.  Remarkably, this structure may
be found in well-motivated models that have been explored previously
for many other reasons~\cite{Feng:2008ya}.  As an example, consider
supersymmetric models with gauge-mediated supersymmetry breaking
(GMSB)~\cite{Dine:1981za,Dimopoulos:1981au,Nappi:1982hm,%
AlvarezGaume:1981wy,Dine:1994vc,Dine:1995ag}.  These models
necessarily have several sectors, as shown in \figref{sectors}.  The
supersymmetry-breaking sector includes the fields that break
supersymmetry dynamically and mediate this breaking to other sectors
through gauge interactions.  The MSSM sector includes the fields of
the minimal supersymmetric standard model.  In addition, supersymmetry
breaking may also be mediated to one or more hidden sectors.  The
hidden sectors are not strictly necessary, but given the discussion
above, there is no reason to prevent them, and hidden sectors are
ubiquitous in such models originating in string theory.

\begin{figure}[tb]
\begin{center}
\includegraphics*[width=9cm]{sectors.eps}
\caption{Sectors of supersymmetric models.  Supersymmetry breaking is
mediated by gauge interactions to the MSSM and the hidden sector,
which contains the dark matter particle $X$.  An optional connector
sector contains fields $Y$, charged under both MSSM and hidden sector
gauge groups, which induce signals in direct and indirect searches and
at colliders.  There may also be other hidden sectors, leading to
multi-component dark matter.  {}From Ref.~\protect\cite{Feng:2008ya}.
}
\label{fig:sectors}
\end{center}
\end{figure}

As is well-known, GMSB models generate superpartner masses
proportional to gauge couplings squared.  Slightly more precisely, the
MSSM superpartner masses are
\begin{equation}
\label{mmass}
m \sim \frac{g^2}{16 \pi^2} \frac{F}{M} \ ,
\end{equation}
where $g$ is the largest relevant standard model gauge coupling, and
$F$ and $M$ are the vacuum expectation values of the
supersymmetry-breaking sector's chiral field $S$, with $\langle S
\rangle = M + \theta^2 F$.  With analogous couplings of the hidden
sector fields to hidden messengers, the hidden sector superpartner
masses are
\begin{equation}
\label{mxmass}
m_X \sim \frac{g_X^2}{16 \pi^2} \frac{F}{M} \ ,
\end{equation}
where $g_X$ is the relevant hidden sector gauge coupling.  As a
result,
\begin{equation}
\frac{m_X}{g_X^2} \sim \frac{m}{g^2} \sim
\frac{F}{16 \pi^2 M} \ ;
\label{mxgx}
\end{equation}
that is, $m_X/g_X^2$ is determined solely by the
supersymmetry-breaking sector.  As this is exactly the combination of
parameters that determines the thermal relic density of \eqref{omega},
the hidden sector automatically includes a dark matter candidate that
has the desired thermal relic density, irrespective of its mass.

The freeze out of hidden sector dark matter in such GMSB models has
been studied numerically.  As an example, in Ref.~\cite{Feng:2008mu}
the hidden sector was assumed to be a copy of the MSSM, but with a
free superpartner mass scale $m_X$ and all Yukawa couplings $\sim
{\cal O}(1)$, so that the only light hidden sector particles are the
hidden gluon, photon, and neutrinos.  The results are given in
\figref{gxmx}.  The criterion that the standard model weak scale be
between 100 GeV and 1 TeV requires values of $(m_X, g_X)$ within the
band.  The solid curves, where the thermal relic density of hidden
dark matter is consistent with dark matter, are seen to lie within this
band, confirming the scaling arguments and rough estimates described
above.

\begin{figure}[tb]
\begin{center}
\includegraphics*[width=10cm]{gxmx.eps}
\end{center}
\vspace*{-.29in}
\caption{Contours of $\Omega_X h^2 = 0.11$ in the $(m_X,g_X)$ plane
for hidden to observable reheat temperature ratios $T^h_{\text{RH}}/T_{\text{RH}}=0.8$
(upper solid) and 0.3 (lower solid), where the hidden sector is a
1-generation flavor-free version of the MSSM.  Also plotted are lines
of $\mweak \equiv (m_X/g^2_X)g'^2 = 100~\,\mbox{GeV}$ (upper dashed) and 1 TeV
(lower dashed).  {}From Ref.~\cite{Feng:2008mu}.  }
\label{fig:gxmx}
\end{figure}

In summary, well-known frameworks for hidden sectors include models in
which the hidden sector contains a particle whose thermal relic
density is automatically in the desired range to be dark matter, even
when the particle's mass is not at the weak scale.  This property
relies on the relation $m_X \propto g_X^2$, which is common to other
frameworks for new physics that avoid flavor-changing problems, such
as anomaly-mediated supersymmetry breaking.  The ``coincidence''
required for WIMPless dark matter may also be found in other settings;
see, for example, Ref.~\cite{Hooper:2008im}.

\subsection{Direct and Indirect Detection Signals}

The decoupling of the WIMP miracle from WIMPs has many possible
implications and observable consequences.  In the case that the dark
matter is truly hidden, it implies that there are no prospects for
direct or indirect detection.  Signals must be found in astrophysical
observations, as in the case of superWIMPs.  Alternatively, there may
be connector sectors containing particles that mediate interactions
between the standard model and the hidden sector through non-gauge
(Yukawa) interactions (see \figref{sectors}).  Such connectors may
generate many signals with energies and rates typically unavailable to
WIMPs.

As an example, first consider direct detection.  The DAMA signal,
interpreted as spin-independent, elastic scattering, has
conventionally favored a region in the mass-cross section plane with
$(m_X, \sigma_{\text{SI}}) \sim (20 - 200~\,\mbox{GeV},
10^{-5}~\,\mbox{pb})$~\cite{Belli:1999nz}.  This is now excluded, most
stringently by XENON10~\cite{Angle:2007uj} and CDMS
(Ge)~\cite{Ahmed:2008eu}, which require $\sigma_{\text{SI}} < 10^{-7}~\,\mbox{pb}$
throughout this range of $m_X$.  Gondolo and Gelmini have noted,
however, that an alternative region with $(m_X, \sigma_{\text{SI}}) \sim
(1-10~\,\mbox{GeV}, 10^{-3}~\,\mbox{pb})$ may explain the DAMA results without
violating other known bounds~\cite{Gondolo:2005hh}. DAMA's relative
sensitivity to this region follows from its low energy threshold and
the lightness of Na nuclei.  This region may be extended to lower
masses and cross sections by the effects of
channeling~\cite{Bernabei:2007hw,Petriello:2008jj,Bottino:2008mf,%
Savage:2008er,Fairbairn:2008gz} and may also be broadened by dark
matter streams in the solar neighborhood~\cite{Gondolo:2005hh}.

The acceptable DAMA-favored region with $m_X \sim 5~\,\mbox{GeV}$ has masses
that are low for WIMPs, given that their masses are expected to be
around $\mweak \sim 100~\,\mbox{GeV} - 1~\,\mbox{TeV}$.  However, in WIMPless models,
where the thermal relic density is achieved for a variety of dark
matter masses, such masses are perfectly natural.  A WIMPless particle
$X$ may couple to the standard model through Yukawa interactions
\begin{equation}
{\cal L} = \lambda_f X \bar{Y}_L f_L
+ \lambda_f X \bar{Y}_R f_R \ ,
\label{connector}
\end{equation}
where $Y$ is a vector-like connector fermion, and $f$ is a standard
model fermion.  Taking $f$ to be the $b$ quark, and the $Y$ mass to be
400 GeV, consistent with current bounds, these couplings generate
spin-independent scattering cross sections given in
\figref{superkdirect2_mod}.  We see that WIMPless dark matter may
explain the DAMA results without difficulty.

\begin{figure}[tb]
\includegraphics*[width=10cm]{superkdirect2_mod.eps}
\caption{Direct detection cross sections for spin-independent
$X$-nucleon scattering as a function of dark matter mass $m_X$.  The
magenta shaded region is DAMA-favored given channeling and no
streams~\cite{Petriello:2008jj}, and the medium green shaded region is
DAMA-favored at 3$\sigma$ given streams and no
channeling~\cite{Gondolo:2005hh}.  The light yellow shaded region is
excluded by the direct detection experiments indicated.  The blue
cross-hatched region is the parameter space of WIMPless models with
connector mass $m_Y = 400~\,\mbox{GeV}$ and $0.3 <\lambda_b < 1.0$.  The black
solid line is the published Super-K exclusion
limit~\cite{Desai:2004pq}, and the black dashed line is a projection
of future Super-K sensitivity.  {}From Ref.~\cite{Feng:2008qn}.
\label{fig:superkdirect2_mod}
}
\end{figure}

WIMPless dark matter also provides new target signals for indirect
detection.  For WIMPs, annihilation cross sections determine both the
thermal relic density and indirect detection signals.  The thermal
relic density therefore constrains the rates of indirect detection
signals.  In the WIMPless case, however, this connection is weakened,
since the thermal relic density is governed by hidden sector
annihilation and gauge interactions, while the indirect detection
signals are governed by the interactions of \eqref{connector}.

This provides a wealth of new opportunities for indirect detection.
As an example, WIMPless dark matter may be detected through its
annihilation to neutrinos in the Sun by experiments such as
Super-Kamiokande.  Although such rates depend on the competing cross
sections for capture and annihilation, the Sun has almost certainly
reached its equilibrium state, and the annihilation rate is determined
by the scattering cross section~\cite{Desai:2004pq}.  The prospects
for Super-Kamiokande may therefore be compared to direct detection
rates~\cite{Desai:2004pq,Hooper:2008cf,Feng:2008qn}.  The results are
given in \figref{superkdirect2_mod}.  In the near future,
Super-Kamiokande may be able to probe the low mass regions
corresponding to the DAMA signal.

WIMPless dark matter also provides additional targets for indirect
detection experiments looking for photons, positrons, and other
annihilation products.  The connectors may also play an interesting
role in collider experiments.  Further details may be found in
Refs.~\cite{Feng:2008dz,Hooper:2008cf,Feng:2008qn}.

High Energy Physics - Phenomenology (hep-ph); Cosmology and Nongalactic Astrophysics (astro-ph.CO)

Non-WIMP Candidates