A High Resolution Search for Dark-Matter Axions

A High Resolution Search for Dark-Matter Axions

L. D. Duffy, P. SikivieB and D. B. Tanner
Department of Physics, University of Florida, Gainesville, Florida 32611, USA

S. J. Asztalos, C. Hagmann, D. Kinion, L. J Rosenberg, K. van Bibber, and D. B. YvQ

Lawrence Livermore National Laboratory, Livermore, California 94550, USA

R. F. Bradley
National Radio Astronomy Observatory, Charlottesville,

(Dated: March 3, 2006)

Vtrgima 22903, USA

O
O

(N

>
00

O

m
o

o

:^

o

X

We have performed a high resolution search for galactic halo axions in cold flows using a microwave
cavity detector. The analysis procedure and other details of this search are described. No axion
signal was found in the mass range 1.98-2.17 /xeV. We place upper limits on the density of axions
in local discrete flows based on this result.

PACS numbers: 14.80.Mz, 95.35.+d, 98.35.Gi

I. INTRODUCTION

In the current concordance cosmology, 23% of the uni-
verse's total energy density is contributed by exotic dark
matter 0. The axion, arising from the Peccei-Quinn
solution to the strong CP problem satisfies the two
criteria necessary for cold dark matter Q]: (1) a very cold
population of axions could be present in our universe in
sufficient quantities to provide the required dark matter
energy density and (2) axions are effectively coUisionless;
i.e., their only significant long-range interaction is gravi-
tational. The Axion Dark Matter eXperiment (ADMX)
uses a Sikivie microwave cavity detector to search
for axions in our galactic halo.

The power in an axion signal observed by a microwave
cavity detector is proportional to the local axion density.
The signal width is caused by the velocity dispersion of
dark-matter axions. Therefore, in searching for axions
it is necessary to make some assumptions about their
velocity distribution in our galactic halo. A variety of
galactic halo models have been put forward: the isother-
mal model, results from N-body simulations and the
caustic ring model The predictions of these models

are used to guide ADMX's search.

In the isothermal model it is expected that a signifi-
cant fraction of the dark matter halo will have an isother-
mal velocity distribution resulting from a period of "vio-
lent relaxation" of the early galaxy This component
of the halo will have velocities described by a MaxwcU-
Boltzmann distribution. The "medium resolution" (MR)
channel 10] searches for these axions, assuming that

*On leave during Fall 2005 at School of Natural Sciences, Institute
for Advanced Study, Princeton, New Jersey 08540, USA, and dur-
ing Spring 2006 at Theoretical Physics Division, CERN, CH1211
Geneve 23, Switzerland

tPresent address: Piper Jaffray, LLC, 245 Park Ave., 33rd floor.
New York, New York 10167, USA

the velocity dispersion is 0(10 ^c) or less. (The es-
cape velocity from our galaxy for axions is approximately
2 X lO-^c.)

N-body simulations and the caustic ring model both
predict substructure within halos. Numerical simulations
indicate that hundreds of smaller clumps, or subhalos,
exist within the larger halo 6] . Tidal disruption of these
subhalos leads to flows in the form of "tidal tails" or
"streams" . The Earth may currently be in a stream of
dark matter from the Sagittarius A dwarf galaxy .

Non-thermalized flows from late infall of dark matter
onto the halo are also expected . Insufficient time has
elapsed for dark matter that has fallen into the gravita-
tional potential of the galaxy after violent relaxation to
thermalize with the rest of the halo. Matter which has
fallen onto the galaxy only recently will be present in the
halo in the form of discrete fiows. There will be one fiow
of particles falling into the gravitational potential for the
first time, one flow of particles falling out for the first
time, one due to particles falling in for the second time,
etc. Furthermore, where the gradient of the particle ve-
locity diverges, particles "pile up" and form caustics. In
the limit of zero flow velocity dispersion, caustics have
inflnite particle density. The velocity dispersion of cold
axions at a time, t, prior to galaxy formation is approxi-
mately (Jwa ~ 3 X 10-1^(10-5 eV/ma)[to/tf/^ 1^, where
to is the present age of the universe and rria is the axion
mass, constrained to lie between IQ-^ and 10^^ eV by
cosmology and astrophysical processes Thus, a flow
of dark matter axions will have a small velocity disper-
sion, leading to large, but finite density at the location
of a caustic.

The caustic ring model predicts that the Earth is lo-
cated near a caustic feature 0|. Fitting the model
to bumps in the Milky Way rotation curve and a tri-
angular feature seen in the IRAS maps predicts that
the flows falling in and out of the halo for the fifth
time contain a significant fraction of the local halo den-
sity. The predicted densities are 1.7 x 10"^** g/crn^ and
1.5 X IQ-^^ g/cm^ IT^, comparable to the local dark

2

matter density of 9.2 x 10~^^ g/cw? predicted in [l5l |.
The flow of the greatest density is referred to as the
"Big Flow" . The possible existence of discrete flows, or
streams, provides an opportunity to increase the ADMX
discovery potential. A discrete axion flow produces a nar-
row peak in the spectrum of microwave photons in the ex-
periment and such a peak can be searched for with higher
signal-to-noise than a signal from axions in an isother-
mal model halo. The "high resolution" (HR) channel
was built to take advantage of this opportunity. Further-
more, if a signal is found, the HR channel will provide us
with detailed information on the structure of the Milky
Way halo.

The HR channel is the most recent addition to ADMX,
implemented as a simple addition to the receiver chain,
running in parallel with the MR channel. This channel
and the possible existence of discrete flows can improve
ADMX's sensitivity by a factor of 3 T^, signiflcantly
enhancing its discovery potential. The full ADMX de-
tector is described in Section ^1 Each discrete flow of
cold axions with small velocity dispersion will be seen as
a narrow peak in the detector's output spectrum. Our ex-
pectations for a signal are discussed in Section [llll which
also contains the details of the HR analysis, the primary
topic of this paper. After a full search of the frequency
range 478-525 MHz, no axion signals were found and
we place limits on the density of cold axions in discrete
flows in Section llVI This limit is compared to our pre-
vious results for the MR channel and halo substructure
predictions in Section Ivl

II. AXION DARK MATTER EXPERIMENT

ADMX uses a microwave cavity detector to search for
axions in our galactic halo. We outline the principle of
the detector and briefly describe ADMX. Further details
of the experiment can be found in ^ .19.] .

The microwave cavity detector utilizes the axion-
electromagnetic coupling to induce resonant conversion
of axions to photons. The relevant interaction is

*^ay^ — y7 « XL/ ■ J3 ,

TT Ja

(1)

where E and B are the electric and magnetic fields,
a is the fine structure constant, fa is the axion de-
cay constant, a{x) is the axion field and g-y is a model-
dependent coupling, of order one. In the Kim-Shifman-
Vainshtein-Zakharov (KSVZ) model T?!, = -0.97,
and in the Dine-Fischler-Srednicki-Zhitnitsky (DFSZ)
model g-y = 0.36. Axions in the galactic halo are
non-relativistic, i.e. the energy of a single axion with
mass rua and velocity v is

Ea = maC^ + -niaV^

(2)

where c is the speed of light. The axion-to-photon con-
version process conserves energy, i.e. an axion of energy

Ea converts to a photon of frequency v = Ea/h. When
h' falls within the bandwidth of a cavity mode, the con-
version process is resonantly enhanced. The signal is a
peak in the spectrum output by the detector. The power,
P, developed in the cavity due to resonant axion-photon
conversion is

P =

\7rfa

VBlpaC

min((3, Qa)

(3)

where V is the cavity volume, Bq is the magnetic field
strength, pa is the local density of axions with energy
corresponding to the cavity frequency, Q is the loaded
quality factor of the cavity, Qa is the ratio of the energy
of the halo axions to their energy spread, equivalent to a
"quality factor" for the halo axion signal, and C is a mode
dependent form factor which is largest for the fundamen-
tal transverse magnetic mode, TMoiq. The quantity C
is given by

C =

\jy (fix'E^ ■ Bo I
BlV jy<Pxe\E^

(4)

in which Eij(x)e*'^* is the time dependent electric field of
the mode under consideration, Bo(x) is the static mag-
netic field in the cavity and e is the dielectric constant of
the medium inside the cavity. The frequency-dependent
form factor is evaluated numerically. Eq. can be re-
cast in the convenient form

P= 0.5x10"

Wi

500

O.SxlO-^* g.cm-

(Bo_Y r ( g-' ^

inin[Q,Qn

V IGHz;

(5)

where Va is the axion mass frequency. As the experi-
ment operates with the cavity at critical coupling, half
the power developed in the cavity is lost to its walls and
half is passed to the receiver chain.

The HR channel became fully operational in August,
2002. A schematic of ADMX, showing both the MR and
HR channels, is given in Fig.^ A more detailed illustra-
tion of the magnet, cavity and cryogenic components is
shown in Fig. |21

The microwave cavity has an inner volume, V, of 189 L.
The frequency of the TMqiq mode can be tuned by mov-
ing a pair of metal or dielectric rods inside the cavity.
The cavity is located in the bore of a superconducting
solenoid, which generates a static magnetic field, Bq, of
7.8 T. The voltage developed across a probe coupled to
the electromagnetic field inside the cavity is passed to the
receiver chain. During typical operation, Q is approxi-
mately 7 X lO'* and the total noise temperature for the
experiment, T„, is conservatively estimated to be 3.7 K,
including contributions from both the cavity and the re-
ceiver chain.

The first segment of the receiver chain is common to
both the MR and HR channels. It consists of a cryo-
genic GaAs HFET amplifier built by NRAO, a crystal
bandpass filter and mixers. At the end of this segment.

3

@@rH>

IMAGE

REJECT 10.7MHz MIXER
MIXER #1 IF #2

1.3K

A

HFET
AMPLIFIER

CAVITY &
TUNING RODS

MAGNET

35kHz AF

MIXER

#3

A

5kHz AF

125Hz
FR

BIN

0.02Hz
FR

BIN

DISK

FIG. 1: Schematic diagram of the receiver chain.

Vacuum Pump

Stepping motors

Cryostat vessel
Magnet support
Cavity LHe reservoir

Magnet LHe reservoir

1.3K J-T refrigerator

Cavity vacuum ctiamber
AmpIiHers

Tuning mectianism
Microwave Cavity

Dielectric tuning rod
Metal tuning rod

Superconducting magnet

FIG. 2: Sketch of the ADMX detector.

the signal is centered at 35 kHz, with a 50 kHz span.
The MR signal is sampled directly after this part of the
receiver chain. The HR channel contains an additional
bandpass filter and mixer, resulting in a spectrum cen-
tered at 5 kHz with a 6 kHz span.

Time traces of the voltage output from the HR channel,
consisting of 2^" data points, are taken with a sampling
frequency of 20 kHz. This results in a data stream of
52.4 s in length, corresponding to 0.019 Hz resolution in
the frequency spectrum. The data were primarily taken
in parallel with the operations of the MR channel over
a period beginning in November, 2002 and ending May,
2004. Continuous HR coverage has been obtained and

candidate peak elimination performed for the frequency
range 478-525 MHz. Data with Q less than 40 000 and/or
cavity temperature above 5 K were discarded. In these
cases, additional data were taken to ensure coverage of
the entire range.

III. HIGH RESOLUTION ANALYSIS

We use the HR channel to search for narrow peaks
caused by flows of cold axions through the detector. The
background is thermal and electronic noise plus narrow
lines from the rf environment of the experiment. These
environmental peaks are signals from nearby emitters
(e.g. computer clocks) which leak into the cavity by a
variety of means.

When placing limits on cold flows of axions, we as-
sume that the flows are steady, i.e. the rates of change
of velocity, velocity dispersion and density of the flows
are slow compared to the time scale of the experiment.
The assumption of a steady flow implies that the signal
we are searching for is always present. Even so, the ki-
netic energy term in Eq. Q and the corresponding signal
frequency change over time due to the Earth's rotational
and orbital motions. In addition to a signal frequency
shift in data taken at different times, apparent broaden-
ing of the signal occurs because its frequency shifts while
the data are being taken.

In this section, we first describe the signal frequency
shifts due to the Earth's rotation and orbital motion, and
the associated signal broadening (subsection A). Next,
we describe the properties of the noise in the HR channel
(subsection B). In subsection C, we describe how we cor-
rect each spectrum for systematic effects introduced by
the receiver chain. Finally, in subsection D, we describe
how candidate peaks are selected and shown not to be
caused by cold flows of dark matter axions.

4

A. Axion signal properties

An axion signal will undergo diurnal and annual mod-
ulations due to the Earth's rotation and orbital motion,
respectively [20|. Thus, the frequency at which axions are
resonantly converted to photons will shift. We show here
that this ft'equency modulation cannot move a signal by
an amount which is larger than the detector bandwidth.
As the energy of an axion in the ADMX detector is non-
relativistic (Eq. Q), the shift in frequency of the signal,
A/, due to a change in velocity of the axion flow relative
to the detector, Au, is

A/

fvAv

(6)

We have investigated the magnitude of both annual and
diurnal signal modulation at / = 500 MHz. The velocity
of a dark matter flow relative to the Earth will be in the
range 100-1000 km/s. We chose 600 km/s as a represen-
tative value for the purpose of estimation. For the daily
modulation, we have assumed that the detector is located
at the Earth's equator and that, in the frame in which the
axis of rotation of the Earth is stationary, the flow veloc-
ity is first aligned and then anti-aligned (or vice versa)
with the detector's motion due to the Earth's rotation
over the course of a day. These assumptions result in the
largest possible change in relative velocity between the
detector and the flow due to the Earth's rotation. The
Earth's rotational velocity is 0.4 km/s at the equator.
The resulting daily signal modulation is of order 1 Hz.
For the annual modulation, to again maximize the change
in relative velocity, we have also considered the case of
extreme flow velocity alignment with the Earth's orbital
motion. The Earth's orbital velocity is 30 km/s. The an-
nual modulation produces the larger frequency shift, of
order 100 Hz within a year. The HR channel has a 6 kHz
bandwidth, ensuring that an axion peak will appear in
spectra taken with center frequency equal to a previously
observed axion signal frequency.

The signal broadening, Sf, due to a change 6v in flow
velocity while data are being taken is

fvSv

(7)

The most significant signal broadening is due to the
Earth's rotation. Using the same assumptions as before,
we find that the broadening is at most 4x 10^^ Hz during
the 52 s taken to acquire a single time trace. This is less
than the spectral resolution of 0.019 Hz. The broadening
due to the Earth's orbital motion is only of order 10^** Hz
in this same time interval.

As the signal broadening due to the Earth's rotation
and orbital motion is negligible, we can use Eq. iQ to
relate the width, Sf, of a signal peak to the velocity dis-
persion, 6v, of the axion flow that causes it. In general,
we do not know the velocity dispersion of the cold ax-
ion flows which we search for, although we note that [l^

claims an upper limit of 53 m/s on the velocity disper-
sion of the Big Flow. Subsequently, we do not know the
signal width. To compensate, we perform our search at
multiple resolutions by combining 0.019 Hz wide bins.
These searches are referred to as n-bin searches, where
n = 1, 2, 4, 8, 64, 512 and 4096. For / = 500 MHz and
V = 300 km/s, the corresponding flow velocity disper-
sions are

Sv„ = 12n m/s

300km/s

(8)

Further details on the n-bin searhces are given in Sec-
tion |mD|

B. Noise in the HR Channel

The power output from the HR channel is expressed
in units of ct, the rms noise power. This noise power is
related to the noise temperature, T„, via

cr = kBT„ \l -

(9)

where fcs is Boltzmann's constant, b is the frequency res-
olution and t is the acquisition time. The total noise tem-
perature Tn = Tc + Tel , where Tc is the physical cavity
temperature and T^i is the electronic noise contribution
from the receiver chain. As no averaging is performed in
HR sampling, b = 1/t. Thus, the rms noise power is

cr = ksbTn

(10)

Output power is normalized to a and r„ is used to deter-
mine this power. We verified Eq. (|10|) experimentally by
allowing the cavity to warm and observing that a is pro-
portional to Tc- As this is our calibration of the power
output from the cavity, it is important that we under-
stand the noise in the HR channel.

The noise in the HR channel is observed to have an
exponential distribution. We now explain why this is ex-
pected. The noise in a single bin is the sum of indepen-
dent sine and cosine components, as no averaging occurs.
We expect that the noise amplitude, a, for a single com-
ponent (i.e. sine or cosine) has a Gaussian probability
distribution.

da

1

2TTaa

exp

a

2ai

(11)

where CTq is the standard deviation. Indeed, the energy
distribution should be proportional to a Boltzmann fac-
tor, cxp(— E / kT) , and non-relativistic and classical ener-
gies, such as E — mv'^/2 oi E — kx^ /2 are proportional
to squares of the amplitude.

As there are two components per bin, the addition of
n bins is that of 2n independent contributions. The sum
of 2n independent normal-distributed components is de-
scribed by a chi-square distribution with 2n degrees of

5

freedom. Thus, the probabihty distribution for an n-bin
is a x^(2n) distribution. We demonstrate this explicitly
in the following.

The probability distribution, dP/dpn, of observing
noise power p„ in an n-bin is

dP

dpn

' 2n

n

dai

I 1 v~>^" 2

exp(-2^Ej=i«j

Evaluating the above expression.

Pn

dP _

d^ ~ (n-l)!c72"

, P7l

(12)
(13)

For n — 1,

dP
dpi

1 / Pi

which is indeed a simple exponential, as expected.

Using this noise distribution, we can easily see that
the average (rms) noise power in the one bin search, a, is
a — (Ja^. Substituting this in Eq. 1)14(1 . the noise power
distribution function becomes

dP
dpi

1

■ exp

(15)

For each individual spectrum, the baseline noise level,
cr, is determined by plotting the number of frequency
bins, Np^ with power between p and p + Ap against p.
According to Eq. 1(15(1 .

NAp

■ exp

where N is the total number of data points. As

a

In

(16)

(17)

a is the inverse of the slope of the IniVp versus p plot.
Fig. O demonstrates that the data is in good agreement
with this relation for p less than 20a. The deviation of
the data from Eq. l(T7| for p greater than 20a is due
to the fact that our background is not pure noise, but
also contains environmental signals of a non-statistical
nature.

As we combine an increasing number of bins, the noise
power probability distribution approaches a Gaussian, in
accordance with the central limit theorem. The right-
hand side of Eq. I(13() approaches a Gaussian in the limit
of large n. We have examined a large sample of noise
in each n-bin search and verified that it is distributed
according to Eq. ((T^ . Figure ^ illustrates statistics for
the 8-bin search and Fig. [S] shows statistics for the 4096-
bin search, a near Gaussian curve. Figures 13 through [S]
illustrates the progression from exponential to Gaussian
noise power probability distribution.

108

106

o 104

102

10°

10 20 30

Power (a)

40

(14) FIG. 3: Power distribution for a large sample of 1-bin data.

8x10

6x10

CO

.E 4x1 0'-
O
D.

2x10'-

••

<

•

>

•

•

<

»

•

•

•

•

•

•

<

5 10 15 20 25 30 35 40

Power (a)

FIG. 4: Power distribution for a large sample of 8-bin data.

6000

4000

Q.

1000

t

s

t

I

\

f

t

%
\

/

t

i
1

r

3800 3900 4000 4100 4200 4300 4400

Power (tj)

FIG. 5: Power distribution for a large sample of 4096-bin
data.

6

BW=1.21 Hz
^ 8- t = 52sec
g

6 ' '

3

1 1 1 1 1 '

479.9966 479.9970 479.9974

Frequency (MHz)

FIG. 6: An environmental peak as it appears in the MR search
(top) and the 64-bin HR search. The unit for the vertical axis
is the rms power fluctuation in each case.

In addition to examining the behavior of the noise
statistics, we have performed a cross-calibration between
the HR and MR channels. The signal power of an en-
vironmental peak, observed at 480 MHz and shown in
Fig. El was examined in both the HR and MR channels.
The observed HR signal power was (1.8±0.1) x 10~^^ W,
where the error quoted is the statistical uncertainty. The
MR channel observed signal power 1.7 x 10~^^ W, in
agreement with the HR channel. Note that the MR sig-
nal was acquired with a much longer integration time
than that of the HR signal (2000 s for MR versus 52 s
for HR).

The combination of the calibration of the noise power
with cavity temperature, the consistency between ex-
pected and observed noise statistics, and the agreement
of signal power observed in both the HR and MR chan-
nels, makes us confident that the signal power is accu-
rately determined in the HR channel.

C. Removal of systematic effects

There are two systematic effects introduced in the re-
ceiver chain shown in Fig. ^ Two passband filters are
present on the HR receiver chain: one with bandwidth
35 kHz on the shared MR-HR section and a passive LC

1.0-
0.8-

0.6-
0.4-

0.2-
0.0-

2000 4000 6000 8000 10000

FFT frequency (Hz)

FIG. 7: HR fllter response calibration data (512 bin aver-
age). The power has been normalized to the maximum power
output.

filter of bandwidth 6 kHz, seen by the HR channel only.
The combined response of both these filters has been ana-
lyzed and removed from the data. The second systematic
effect is due to the frequency-dependent response of the
coupling between the cavity and the first cryogenic am-
plifier. This effect is removed using the equivalent circuit
model described later.

The combined passband filter response was determined
by taking data with a white noise source at the rf input
of the receiver chain. A total of 872 time traces were
recorded over a two day period. In order to achieve a
reasonably smooth calibration curve, 512 bins in the fre-
quency spectrum for each time trace were averaged giving
9.77 Hz resolution. The combined average of all data is
shown in Fig. |7| This measured response was removed
from all data used in the HR search, as follows. The raw
power spectra have frequency 0-10 kHz, where the center
frequency of 5 kHz has been mixed down from the cav-
ity frequency. Each raw power spectrum is cropped to
the region 2-8 kHz to remove the frequencies not within
the LC filter bandwidth. Each remaining frequency bin
is then weighted by a factor equal to the receiver chain
response at the given frequency divided by the maxi-
mum receiver chain response. Interpolation for frequency
points not specifically included in the calibration curve
is performed by assuming that each point on the cali-
bration curve was representative of 512 bins centered on
that frequency, so all power corresponding to frequencies
within that range is normalized by the same factor. As
the calibration curve varies slowly with frequency within
the window to which each spectra is cropped, this is an
adequate treatment of the normalization.

In the MR channel, the effect of the cavity-amplifier
coupling is described using an equivalent-circuit model
PH . This model has been adapted for use in the HR
channel. The frequency dependent response of the cavity

7

4600

4500

4400-

-S- 4300
(U

I 4200

a.

4100-

4000

3900

494.8100 494.8120 494.8140

Frequency (Hz)

494.8160

FIG. 8: Sample 4096-bin spectrum before correction for the
cavity-amplifier coupling. The line is the fit obtained using
the equivalent circuit model.

04 as it is an overall factor which results from normalizing
the power to the single bin noise baseline. In practice, the
parameters ai through 05 are established by fitting. The
line in Fig. |S1 shows the fit obtained using the equivalent
circuit model.

Large peaks in the data, e.g. an axion signal or environ-
mental peak, are removed before fitting to prevent bias.
The 4096-bin spectrum is used to perform the fit and then
the original 1-bin spectrum is corrected to remove the
systematic effect. The weighting factors are calculated
using Eq. I|18|) and the fitted parameters, oi through 05,
at the center of each bin of width &4096- These factors
are the ratio of the fit at a given point to the maximum
value of the fit. Each 1-bin is multiplied by the factor
calculated for the bin of width 64096 within which it falls.

The removal of the cavity-amplifier coupling and the
passband filter response using the techniques described
above has been demonstrated to result in flat HR spectra.

D. Axion signal search procedure

amplifier coupling is most evident in the 4096-bin search,
thus this is the data used to apply the equivalent circuit
model. A sample spectrum before correction is shown in
Fig. El

In the equivalent-circuit model, each frequency is given
by A, the number of bins it is offset from the bin of
the center frequency, measured in units of the 4096-bin
resolution, i.e. 64095 = 78.1 Hz. The equivalent-circuit
model predicts that the power (in units of the rms noise)
at the NRAO amplifier output (the point labelled "RF"
in Fig. ^ in the 4096-bin search at the frequency offset
A is

«i+8a3(^)'+4a4(^)
P(A) = "^-^ rr^^^, (18)

where the parameters ai through are

ai = (64096/6) (^c + Tj + Tv)/Tn , (19)

a2 = /o/(&4096 Q) , (20)

03 = (64096/6)(T/ +Ty + (T,-Ty )cos(2fcL))/r^21)

04 = (64096/6)((T/-ry)sin(2fcL))/T„and (22)

05 — ifo — /cen)/64096 • (23)

In the above expressions, Tc is the physical temperature
of the microwave cavity, Tj and Ty are the current and
voltage noise, respectively, contributed by the amplifier,
Tn is the noise temperature contributed from all compo-
nents, 6 is the frequency resolution of the HR channel,
i.e. 0.019 Hz, L is the electrical (cable) length from the
cavity to the HFET amplifer, /g is the cavity resonant fre-
quency, fcen is the center frequency of the spectrum and k
is the wavenumber corresponding to frequency fcen + 6A.
The factor 64096/6 appears in the parameters ai, 03 and

We now describe the search for an axion signal and
summarize the analysis performed on each time trace.

The width of an axion signal is determined by the sig-
nal frequency, axion velocity and flow velocity dispersion
(Eq. {T))), the latter being the most uncertain variable.
n-hin searches, where n is the number of adjacent 1-bins
added together {n = 1, 2, 4, 8, 64, 512 and 4096), are
conducted to allow for various velocity dispersions. For
searches with n > 1, there is an overlap between succes-
sive n-bins such that each n-bin overlaps with the last
half of the previous and first half of the following n-bin.
This scheme is illustrated for the 2, 4 and 8-bin searches
in Fig. El

The search for an axion signal is performed by scanning
each spectrum for peaks above a certain threshold. All
such peaks are considered candidate axion signals. The
thresholds are set at a level where there is only a small
probability that a pure noise peak will occur and such
that the number of frequencies considered as candidate
axion peaks is manageable. The candidate thresholds
used were 20, 25, 30, 40, 120, 650 and 4500 a, in increas-
ing order of n.

All time traces are analyzed in the same manner. A
fast Fourier transform is performed and an initial esti-
mate of a is obtained by fitting the 1-bin noise distribu-
tion to Eq. (|17l) . Systematic effects are then removed,
i.e. the corrections described in Section Fill ("I for the fil-
ter passband response and cavity-amplifier coupling are
performed. "Large" peaks not included in the equivalent
circuit model fit for the cavity-amplifier response are de-
fined to be those greater than 120% of the search thresh-
old for each n-bin search. After the removal of system-
atic effects, the 1-bin noise distribution is again fitted to
Eq. ((T7|) to obtain the true value of a and the search for
peaks above the thresholds takes place.

The axion mass is not known, requiring that a range

8

1-bin search: 1 2

3

4

5

6

7

8

9

10

11

12

13

14

15 16

2-bin search:

1 2

3 4

5 6

7 8

9 10

11 12

13 14

15 16

2 3

4 5

6 7

8 9

10 11 ]

L2 13

14 15

4-bin search:

12 3 4

5

6 7 8

9 10 11 12

13 14 15 16

3 4 5

6

7 8 9 10

11 12 13 14

8-bin search:

12 3 4 5 6 7 8

9 10 11 12 13 14 15 16

5 6 7 8 9 10 11 12

FIG. 9: Illustration of the addition scheme for the 2, 4 and 8-bin searches. The numbers correspond to the data points of the
1-bin search. Numbers within the same box are bins added together to form a single datum in the n-bin searches with n > 1.

of frequencies must be examined. Full HR coverage has
been obtained for the region 478-525 MHz, correspond-
ing to axion masses between 1.98 and 2.17 /^eV. The se-
lected frequency range is examined in three stages for
axion peaks, as follows:

Stage 1: Data for the entire selected frequency range
is taken. The frequency step between successive spectra
is approximately 1 kHz, i.e. the center frequency of each
spectrum differs from the previous spectrum by 1 kHz.
Frequencies at which candidate axion peaks occur are
recorded for further examination during stage 2.

Stage 2: Multiple time traces are taken at each candi-
date frequency from stage 1. The steady flow assumption
described in Section [ill Al means that a peak will appear
in spectra taken with center frequency equal to the can-
didate frequency from stage 1 if such a peak is an axion
signal. The frequencies of persistent peaks, i.e. peaks
that appear during both stage 1 and 2 are examined fur-
ther in stage 3.

Stage 3: Frequencies of persistent peaks undergo a
three-part examination. The first step is to repeat stage
2, to ensure the peaks still persist. Secondly, the warm
port attenuator is removed from the cavity and multiple
time traces taken. If the peak is due to external radio
signals entering the cavity (an environmental peak), the
signal power will increase dramatically. If the signal orig-
inates in the cavity due to axion-photon conversion, the
power developed in the cavity will remain the same as
that for the normal configuration. The third step is to
use an external antenna probe as a further confirmation
that the signal is environmental. Some difficulties were
encountered with the antenna probe, due to polarization
of environmental signals. However, the second step is
adequate to confirm that peaks are environmental. If a
persistent peak is determined to not be environmental,
a final test will confirm that it is an axion signal. The
power in such a signal must grow proportionally with the
square of the magnetic field {Bo in Eq. jSJ) and disappear
when the magnetic field is switched off.

No axion peaks were found in the range 478-525 MHz
using this approach. The exclusion limit calculated from
this data is discussed in the following section.

IV. RESULTS

Over the frequency range 478-525 MHz, we derive an
upper limit on the density of individual fiows of axion
dark matter as a function of the velocity dispersion of
the flow. The corresponding axion mass range is 1.97-
2.17 /xeV. Each n-bin search places an upper limit on the
density of a flow with maximum velocity dispersion,
as given by Eq. ||SJ|.

Several factors reduce the power developed in an ax-
ion peak from that given in Eq. JSJ. The experiment is
operated near critical coupling of the cavity to the pream-
plifier, so that half this power is observed when the cavity
resonance frequency, /o, is precisely tuned to the axion
energy. If /o is not at the center of a 1-bin, the power is
spread into adjacent bins, as discussed below. When the
axion energy is off-resonance, but still within the cav-
ity bandwidth at a frequency /, the Lorentzian cavity
response reduces the power developed by an additional
factor of

Kf)

1

1

1

(24)

To be conservative, we calculate the limits at points
where successive spectra overlap, i.e. at the frequency
offset from /q that minimizes h{f).

If a narrow axion peak falls at the center of a 1-bin,
all power is deposited in that 1-bin. However, if such
a peak does not fall at the center of a 1-bin, the power
will be spread over several 1-bins. We now calculate the
minimum power in a single n-bin caused by a randomly
situated, infinitely narrow axion line. The data recorded
is the voltage output from the cavity as a function of
time. The voltage as a function of frequency is obtained
by Fourier transformation and then squared to obtain a
raw "power" spectrum. The actual power is obtained
by comparison to the rms noise power. The data are
sampled for a finite amount of time and thus, the Fourier
transformation of the output, ^{f), will be of the voltage
multiplied by a windowing function, i.e.

v{t)w{t) exp{i2TTft)dt ,

(25)

where v(t) is the measured output voltage and w{t) is

9

the windowing function for a sampling period T,

^^^■^ I otherwise . ^"^^^
Eq. is equivalent to

/•oo

(27)

where V{f) and M^(/) are the Fourier transforms of the
output voltage, v{t), and the windowing function, w{t),
i.e. !F{f) is the convolution of V{f) and W{f), given by

sin(^/r)

(28)

Discretizing Eq. H27|l and inserting Eq. H28|) . we have

AT

^0 2 ^(i-("^+2))

where b is the frequency resolution of the HR channel,
2N points are taken in the original time trace, and the
center frequency of the jth 1-bin is {j + 1/2)6. Thus,
for an axion signal of frequency / falling in 1-bin j, a
fraction of the power

aim)

sin(m7r + 5)
rmr + S

(30)

is lost to the rnth 1-bin from 1-bin j, where 6 — Trim +
1/2 — f/b). If (5 = 0, i.e the axion signal frequency is
exactly equal to a 1-bin center frequency, all the power
is deposited in a single 1-bin. However, if this is not the
case, power is lost to other 1-bins. In setting limits, we
assume that the power loss is maximal.

The maximum power loss occurs when a signal in the
1-bin search falls exactly between the center frequency of
two adjacent 1-bins. In this case, when 6 = tt/2, Eq. H30|)
shows that 40.5% of the power will be deposited in each
of two 1-bins. In n-bin seaches with n > 2, not as much
power is lost to other n-bins, due to the overlap between
successive n-bins. The minimum power deposited in an
n-bin is 81% for n = 2, 87% for n = 4 and 93% for
n — 8. For n = 64, 512 and 4096, the amount of power
not deposited in a single n-bin is negligible.

For the n-bin searches with n — 64, 512 and 4096, a
background noise subtraction was performed which will
lead to exclusion limits at the 97.7% confidence level.
These limits are derived using the power at which the
sum of the signal power and background noise power
have a 97.7% probability to exceed the candidate thresh-
olds. We call this power the "effective" threshold for each
search. The effective thresholds are obtained by integrat-
ing the noise probability distribution, Eq. 113() . numeri-
cally solving for the background noise power correspond-
ing to the 97.7% confidence level for each n and subtract-
ing these values from the original candidate thresholds.

TABLE I: Effective power thresholds for all n-bin searches,
with the frequency resolutions, 6„ and corresponding max-
imum flow velocity dispersions, Sv„, for a flow velocity of
600 km/s.

n

Effective

threshold (a)

(Hz)

(m/s)

I

20

0.019

6

2

25

0.038

10

4

30

0.076

20

8

40

0.15

50

64

71

1.2

400

512

182

9.8

3000

4096

531

78

20000

TABLE II:

Numerically calculated values of the form factor,

C, and amplifler noise temperatures, T^i, from NRAO speci-

flcations.

Frequency (MHz)

C

450

0.43

1.9

475

0.42

1.9

500

0.41

1.9

520

0.38

1.9

550

0.36

2.0

For n = 64, 512 and 4096, the effective thresholds are 71,
182 and 531 cr, respectively. For smaller values of n, back-
ground noise subtraction does not significantly improve
the limits and the effective threshold was taken to be the
candidate threshold. Table summarizes this informa-
tion and shows the frequency resolution of each search
with the corresponding maximum flow velocity disper-
sion from Eq. for w = 600 km/s.

Our exclusion limits were calculated for an axion signal
with power above the effective threshold reduced by the
appropriate factors. These factors arise from the critical
coupling, the Lorentzian cavity response and the maxi-
mum power loss due to the peak not falling in the center
of an n-bin, as outlined above. Equations (O and H10|)
were used, for both KSVZ and DFSZ axion couplings.
The cavity volume, V, is 189 L. Measured values of the
quality factor, Q, the magnetic field, Bq, and the cavity
temperature, Tc, are recorded in each data file. Nu-
merically determined values of the form factor, C are
given in Table UTI The electronic noise temperature, Tej,
was conservatively taken from the specifications of the
NRAO amplifier, the dominant source of noise in the re-
ceiver chain, although our measurements indicate that
Tel is less than specified. These values are also given in
Table Ull Linear interpolation between values at the fre-
quencies specified was used to obtain values of C and T^i
at all frequencies.

The 2-bin search density exclusion limit obtained using
these values is shown in Fig. ^] For values of n other
than n = 2, the exclusion limits differ by only constant
factors. The constant factors are 1.60, 1.00, 1.12, 1.39,
2.53, 5.90 and 17.2 for n = 1, 2, 4, 8, 64, 512 and 4096,

10

2.00 2.05 2.10 2.15

axion mass (^eV)

FIG. 10: 97.7% confidence level limits for the HR 2-bin search
on the density of any local axion dark matter flow as a func-
tion of axion mass, for the DFSZ and KSVZ 077 coupling
strengths. Also shown is the previous ADMX limit using the
MR channel. The HR limits assume that the flow velocity
dispersion is less than Sv2 given by Eq. (|HJ.

respectively.

V. DISCUSSION

to a signal occurring away from the center of a bin than
in the n = 2 case. For n > 2, the limits are more general,
but the larger power threshold of the searches make them
less stringent.

The largest flow predicted by the caustic ring model
has density 1.7 x 10^24 g/^^^s 95 GeV/cm^), velocity
of approximately 300 km/s relative to the detector, and
velocity dispersion less than 53 m/s Using Eq. (0)
with Table ^ and the information displayed in Fig. [TUI
multiplied by the appropriate factors of 1.12 to obtain
the 4-bin limit, it can be seen that the 4- bin search, cor-
responding to maximum velocity 50 m/s for v — 300
km/s, would detect this flow if it consisted of KSVZ ax-
ions. For DFSZ axions, this flow would be detected for
approximately half the search range.

Figure EH demonstrates that the high resolution anal-
ysis improves the detection capabilities of ADMX when
a significant fraction of the local dark matter density
is due to flows from the incomplete thermalization
of matter that has only recently fallen onto the halo.
The addition of this channel to ADMX provides an
improvement of a factor of 3 over our previous medium
resolution analysis.

We have obtained exclusion limits on the density in
local flows of cold axions over a wide range of velocity
dispersions. The most stringent limit, shown in Fig. 1101
is from the 2-bin search. For a flow velocity of 600 km/s
relative to the detector, the 2-bin search corresponds to
a maximum flow velocity dispersion of 10 m/s. The 1-
bin search limit is less general, in that the corresponding
flow velocity dispersion is half that of the 2-bin limit. It
is also less stringent; much more power may be lost due

Acknowledgements

This research is supported in part by the U.S. De-
partment of Energy under Contract W-7405-ENG-48 at
Lawrence Livermore National Laboratory, under grant
DE-FG02-97ER41029 at the University of Florida, and
by an IBM Einstein Endowed Fellowship at the Institute
for Advanced Study.

[1] C. L. Bennett et ai, Astrophys. J. Suppl. 148, 1 (2003).
[2] R. D. Peccei and H. R. Quinn, Phys. Rev. D 16, 1791

(1977) ; R. D. Peccei and H. R. Quinn, Phys. Rev. Lett.
38, 1440 (1977); S. Weinberg, Phys. Rev. Lett. 40, 223

(1978) ; F. Wilczek, Phys. Rev. Lett. 40, 279 (1978).

[3] L. F. Abbott and P. Sikivie, Phys. Lett. B 120, 133
(1983). J. Preskill, M. B. Wise and F. Wilczek, Phys.
Lett. B 120, 127 (1983). M. Dine and W. Fischler, Phys.
Lett. B 120, 137 (1983).

[4] S. Asztalos et al., Phys. Rev. D64 092003 (2001).

[5] P. Sikivie, Phys. Rev. Lett. 51, 1415 (1983) [Erratum-
ibid. 52, 695 (1984)].

[6] J. F. Navarro, C. S. Frenk and S. D. M. White, As-
trophys. J. 462, 563 (1996); B. Moore, F. Governato,
T. Quinn, J. Stadel and G. Lake, Astrophys. J. 499, L5
(1998).

[7] P. Sikivie, Phys. Lett. B432 139 (1998).
[8] P. Sikivie, Phys. Rev. D 60, 063501 (1999).
[9] D. Lynden-BeU, Mon. Not. Roy. Astron. Soc. 136, 101
(1967).

[10] S. J. Asztalos et al, Phys. Rev. D 69, OlllOl(R) (2004).

[11] K. Freese, P. Gondolo, H. J. Newberg and M. Lewis,
Phys. Rev. Lett. 92, 111301 (2004); K. Freese, P. Gon-
dolo and H. J. Newberg, Phys. Rev. D 71, 043516 (2005).

[12] P. Sikivie and J. Ipser, Phys. Lett. B291 288 (1992);
A. Natarajan and P. Sikivie, Phys. Rev. D 72, 083513
(2005).

[13] J.E. Kim, Phys. Rep. 150 1 (1987); H.-Y. Cheng, Phys.

Rep. 158 1 (1988); M.S. Turner, Phys. Rep. 197 67

(1990); G.G. Raffelt, Phys. Rep. 198 1 (1990).
[14] P. Sikivie, Phys. Lett. B 567, 1 (2003).
[15] E. I. Gates, G. Gyuk and M. S. Turner, Astrophys. J.

449, L123 (1995).
[16] L. Duflty et al, Phys. Rev. Lett. 95, 091304 (2005).
[17] J. E. Kim, Phys. Rev. Lett. 43, 103 (1979); M. A. Shif-

man, A. I. Vainshtein and V. I. Zakharov, Nucl. Phys. B

166, 493 (1980).
[18] M. Dine, W. Fischler and M. Srednicki, Phys. Lett. B

104, 199 (1981). A. R. Zhitnitsky, Sov. J. Nucl. Phys.

31, 260 (1980) [Yad. Fiz. 31, 497 (1980)].

11

[19] S. De Panfilis et al, Phys. Rev. Lett. 59, 839 (1987);
C. Hagmann, P. Sikivio, N. S. Sullivan and D. B. Tanner,
Phys. Rev. D 42, R1297 (1990); C. Hagmann et al., Phys.
Rev. Lett. 80, 2043 (1998); H. Peng et al., Nucl. Instrum.
Meth. A 444, 569 (2000); R. Bradley et al.. Rev. Mod.
Phys. 75, 777 (2003).

[20] M. S. Turner, Phys. Rev. D 42, 3572 (1990); F. S. Ling,
P. Sikivic and S. Wick, Phys. Rev. D 70, 123503 (2004).

[21] E. J. Daw, PhD thesis, Massachusetts Institute of Tech-
nology, 1998.