Skip to main content

Full text of "A High Resolution Search for Dark-Matter Axions"

See other formats


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.