# 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.