We describe the first implementation of a Josephson Traveling Wave Parametric Amplifier (JTWPA) in an axion dark matter search. The operation of the JTWPA for a period of about two weeks achieved sensitivity to axion-like particle dark matter with axion–photon couplings above 10−13 Ge V−1 over a narrow range of axion masses centered around 19.84 µeV by tuning the resonant frequency of the cavity over the frequency range of 4796.7–4799.5 MHz. The JTWPA was operated in the insert of the axion dark matter experiment as part of an independent receiver chain that was attached to a 0.56-l cavity. The ability of the JTWPA to deliver high gain over a wide (3 GHz) bandwidth has engendered interest from those aiming to perform broadband axion searches, a longstanding goal in this field.

We aim to compose 85%2 of the matter content of the universe, dark matter appears necessary to explain a variety of observations, from galactic rotation curves and gravitational lensing to the anisotropies of the cosmic microwave background and the properties of galaxy cluster collisions.2–5 Despite the scope of these observations, the nature of dark matter remains a mystery. One particularly promising candidate is the axion, which has the unique ability not only to account for all the dark matter6–8 but also to solve the so-called strong-CP problem.9–12 

One method to render the axion visible is called the resonant cavity haloscope, first proposed by Pierre Sikivie.13,14 The Axion Dark Matter eXperiment (ADMX) is one such experiment that is unique in its ability to detect Dine–Fischler–Srednicki–Zhitnitsky (DFSZ) axions.1,15,16 Other haloscopes have achieved sensitivity to the band around the Kim–Schifman–Vainshtein–Zakharov (KSVZ) model of the axion.17–27 Of the two, the DFSZ axion is the more difficult to detect, with signals of an order of magnitude smaller than the KSVZ axion. Sensitivity to the DFSZ axion is therefore a noteworthy achievement. Such progress is attributable to the advent of ultra-low-noise quantum amplifiers, such as the Microstrip SQUID Amplifier (MSA)28 and the Josephson Parametric Amplifier (JPA).29,30 A limitation of these amplifiers is their narrowband nature. Recent developments have enabled the fabrication of modified JPAs that can provide bandwidth over as much as 640 MHz.31 Nevertheless, the subsequently developed Josephson Traveling Wave Parametric Amplifier (JTWPA)32 provides power gain over a much wider bandwidth of several GHz, which could further enable broadband axion searches. We note that one other experiment implemented a JTWPA in a receiver chain intended for axion searches,33 though no axion search was performed as it was not placed in a magnetic field.

The JTWPA consists of a lumped element transmission line, where Josephson junctions are the non-linear inductive element. When a microwave signal travels down the line, the non-linear inductance causes four-wave mixing. This feature enables the JTWPA to provide broadband power gain.32 Because of this feature, it is recognized that the JTWPA is generally well-suited for axion searches that aim to cover a wider range of the axion parameter space.34 Broadband experiments are ideal for axion searches involving transient signals appearing off-resonance and allow for the possibility of multimode searches.35,36 Another advantage of the JTWPA is that it does not require a circulator to separate incoming and outgoing modes, enabling a compact receiver design. This helpful feature is ideal for axion searches aiming to take advantage of any extra space in a volume that is typically constrained by the diameter of a magnet bore. This characteristic may be especially advantageous in the case of multi-cavity searches that require a scaling up of microwave components contained within a field-cancellation coil.

In an effort to integrate a JTWPA into an axion search, a prototype experiment was operated within the ADMX magnet bore. The “sidecar” cavity is mounted on top of the ADMX main cavity,1 residing in a magnetic field of about 3.83 T just above the main experiment’s superconducting magnet.37 Although it uses an independent receiver chain, its operation is secondary to the main experiment. Data-taking is therefore likely to cease when technical issues pertaining to the main experiment are encountered. The 0.56- sidecar cavity can be disassembled into two halves, thus allowing a single tuning rod to be mounted inside with ease. Typical measurements of the sidecar loaded quality factor, QL, were ∼700. The quality factor was lower than expected due to larger than expected material losses in the alumina axles as well as leakage through the RF feedthrough holding the readout antenna. The form factor, representing the overlap of the static magnetic and axion electric field, was computed to be 0.41 ± 0.06 using detailed simulations performed in both CST (Computer Simulation Technology) Microwave Studio Suite38 and COMSOL.39 The cavity is coupled by an antenna to a receiver chain that exists independently of the main experiment.

Axion conversion in the magnetic field is expected to produce power in the cavity given by14 
(1)
where V is the volume of the cavity; B is the magnitude of the magnetic field; C010 is the form factor; gγ is the model-dependent numerical constant −0.97 (0.36) for the KSVZ (DFSZ) model that determines, along with the axion decay constant fa, the axion coupling to two photons gaγγ = αgγ/πfa; ρa is the expected dark matter density in the cavity; and f is the frequency of the photon induced by the axion field. The mass of the axion is related to the frequency of the conversion photon by fmc2/h.

At present, the current iteration of the sidecar uses a JTWPA as its first-stage amplifier and a heterostructure field-effect transistor (HFET) amplifier as its second-stage amplifier.40 The diagram of the receiver chain for the sidecar cavity is shown in Fig. 1. The receiver chain allows for transmission and reflection measurements to determine the cavity frequency, quality factor, and antenna coupling. The JTWPA is located on the output line, acting as the first-stage amplifier and conduit for power coming from the cavity. The JTWPA is operated by providing a pump tone coming from a local oscillator in the warm electronics space. The particular JTWPA used in this search was fabricated at and acquired from MIT Lincoln Laboratory.41 It is made from niobium, with a critical temperature of 9.3 K.

FIG. 1.

RF schematic of the sidecar cavity for this run. The JTWPA resides between two circulators on the output line. These two circulators are used to prevent power from reflecting into the cavity. All these components are mounted within a so-called quantum amplifier package that is thermally sunk and bolted to the main ADMX cavity via a cold finger. The package itself resides in a low-field region that is maintained using a field cancellation coil, and the amplifiers are protected with an additional mu metal shielding.

FIG. 1.

RF schematic of the sidecar cavity for this run. The JTWPA resides between two circulators on the output line. These two circulators are used to prevent power from reflecting into the cavity. All these components are mounted within a so-called quantum amplifier package that is thermally sunk and bolted to the main ADMX cavity via a cold finger. The package itself resides in a low-field region that is maintained using a field cancellation coil, and the amplifiers are protected with an additional mu metal shielding.

Close modal

Data were acquired for about two weeks with the sidecar cavity until the main experiment required a magnet ramp down. Our first implementation investigated only a narrow frequency range that was accessible with the sidecar cavity due to a sticking tuning rod. Nevertheless, these results represent the first demonstration of an axion search with a JTWPA.

Throughout the course of data-taking, the JTWPA gain and signal-to-noise ratio improvement (SNRI) were monitored and optimized. The SNRI is computed using42,
(2)
where Gon/Goff is the relative gain of the receiver chain when the JTWPA pump is turned on rather than off, and Poff/Pon is the relative power coming from the cavity with the JTWPA pump is turned on as compared to off.42 Occasional rebiases of the JTWPA were necessary because of changing temperature conditions and mechanical vibrations of the cavity and receiver chain. Therefore, the JTWPA pump frequency and the local oscillator frequency were adjusted as the frequency shifted. A JTWPA SNRI1 greater than 9 dB was maintained over the course of 2 weeks, as well as a measured power gain between 12 and 17 dB over the same time frame, at the frequency of the TM010 mode. This power gain was evident over the range of several GHz. Figure 2 shows the power gain measured over a 3-GHz range.
FIG. 2.

Wideband power gain of the JTWPA from 4 to 7 GHz, measured during data-taking operations. The magnetic field is on, and the JTWPA temperature is 100 mK. The JTWPA provided gain with pump tone frequencies roughly between 7 and 8.5 GHz.

FIG. 2.

Wideband power gain of the JTWPA from 4 to 7 GHz, measured during data-taking operations. The magnetic field is on, and the JTWPA temperature is 100 mK. The JTWPA provided gain with pump tone frequencies roughly between 7 and 8.5 GHz.

Close modal
FIG. 3.

Transmission measurements through the JTWPA, acquired by routing power from the cavity bypass line through the output line, with two different pump frequencies applied to the JTWPA, shown in solid (blue) and dashed (orange) curves. The gain peaks are spaced roughly 38 MHz apart. The JTWPA was biased so that the 2-MHz power spectrum was centered on a region with maximum gain.

FIG. 3.

Transmission measurements through the JTWPA, acquired by routing power from the cavity bypass line through the output line, with two different pump frequencies applied to the JTWPA, shown in solid (blue) and dashed (orange) curves. The gain peaks are spaced roughly 38 MHz apart. The JTWPA was biased so that the 2-MHz power spectrum was centered on a region with maximum gain.

Close modal
To calibrate the sensitivity of the sidecar experiment, it is necessary to understand the system noise temperature, which can be written as43 
(3)
Here, THFET is the noise temperature of the HFET and downstream warm receiver components, hereafter referred to as “downstream receiver noise temperature,” and ɛ < 1 is the transmission efficiency between the JTWPA and the cavity. The value of THFET was measured during a period in which the cavity was warmed up from about 0.2 to 0.6 K. During this measurement, the JTWPA was deactivated by removing the pump tone. The hot load, as shown in Fig. 1, was not used because of an unreliable RF switch. In this scenario, the power measured off cavity resonance is defined by
(4)
where GHFET is the gain of the HFET, Tattn is the physical temperature from the final stage attenuator, labeled A1 in Fig. 1, on the sidecar bypass RF line, ɛc is the transmission efficiency from the final stage attenuator, A1, through the JTWPA, and b is the bandwidth. The temperature Tcav is the physical temperature of the lines leading from the attenuator to the cavity and through the quantum amplifier package. All of the components in this portion of the receiver chain maintained very close temperatures throughout data-taking, as verified with multiple temperature sensors. Our measurement of the downstream receiver noise temperature was acquired during the period over which the experiment was being warmed. Because the final stage attenuator is thermally sunk to the main cavity along with the sidecar cavity, so that Tattn = Tcav, Eq. (4) reduces to
(5)
During the warm-up, power from the receiver was sampled and the temperature of the attenuator was measured. We fit the data using Eq. (5) to extract THFET. Two examples of this fit at frequencies on either side of the resonant frequency can be seen in Fig. 4. The parameter THFET was determined over a range of frequencies by fitting data acquired every 300 MHz. A plot of THFET as a function of frequency can be seen in Fig. 5. We interpolated between nearby points in frequency space to determine the value for THFET at the resonant frequency of 4.798 GHz. We used this value of 3.7 ± 0.2 K as the downstream receiver noise temperature in the analysis. This noise temperature is typical of HFET amplifiers in the presence of a magnetic field.44,45 The system noise was computed using Eq. (3), with this value for THFET.
FIG. 4.

Mean digitized power of sidecar as a function of the temperature of the final stage attenuator for two frequencies off resonance: 4.6 GHz (left) and 4.9 GHz (right). Data were acquired on either side of the resonance because THFET can easily be determined using Eq. (4) off-cavity resonance. The value of THFET on resonance can then be computed by interpolating between the two off-resonance values for THFET. The actual temperature sensor for the x-axis was mounted to the quantum amplifier package, which is thermally connected to the final stage attenuator on the bypass line. The gain of the second-stage and downstream electronics and their associated noise temperature were extracted from the fit. Each data point consisted of an integration time of 10 s. Data are not evenly spaced in temperature space because of the non-linear warming of the ADMX insert. This value of THFET, along with the JTWPA SNRI, was used to compute the system noise.

FIG. 4.

Mean digitized power of sidecar as a function of the temperature of the final stage attenuator for two frequencies off resonance: 4.6 GHz (left) and 4.9 GHz (right). Data were acquired on either side of the resonance because THFET can easily be determined using Eq. (4) off-cavity resonance. The value of THFET on resonance can then be computed by interpolating between the two off-resonance values for THFET. The actual temperature sensor for the x-axis was mounted to the quantum amplifier package, which is thermally connected to the final stage attenuator on the bypass line. The gain of the second-stage and downstream electronics and their associated noise temperature were extracted from the fit. Each data point consisted of an integration time of 10 s. Data are not evenly spaced in temperature space because of the non-linear warming of the ADMX insert. This value of THFET, along with the JTWPA SNRI, was used to compute the system noise.

Close modal
FIG. 5.

Downstream receiver noise temperature (THFET) as a function of frequency. To obtain the downstream receiver noise temperature at any given frequency, power was sampled and integrated for 10 s while the final stage attenuator on the bypass line was heated. The resulting data were fit using Eq. (4). The value of THFET at 4.798 GHz, shown in red, was computed by interpolating between the two nearest values for THFET. This value of THFET, along with the JTWPA SNRI, were used to compute the system noise. The uncertainties on THFET are a combination of the statistical uncertainty and the temperature sensor uncertainty.

FIG. 5.

Downstream receiver noise temperature (THFET) as a function of frequency. To obtain the downstream receiver noise temperature at any given frequency, power was sampled and integrated for 10 s while the final stage attenuator on the bypass line was heated. The resulting data were fit using Eq. (4). The value of THFET at 4.798 GHz, shown in red, was computed by interpolating between the two nearest values for THFET. This value of THFET, along with the JTWPA SNRI, were used to compute the system noise. The uncertainties on THFET are a combination of the statistical uncertainty and the temperature sensor uncertainty.

Close modal

The SNRI is continuously optimized throughout data-taking to minimize the system noise. The wide bandwidth of the JTWPA simplifies the amplifier optimization process. A JTWPA requires only the adjustment of the pump frequency and power, an advantage over a JPA that requires the adjustment of the former together with a bias current to tune its resonant frequency. A feature of the JTWPA is that its gain is not flat across a wide range. Instead, there are oscillations in the gain of about 10 dB, with peaks spaced roughly 38 MHz apart. This feature arises from imperfect impedance matching to the Josephson junctions in the transmission line that constitutes the JTWPA. The exact spacing and height of the peaks depend on the individual JTWPA. Adjusting the frequency of the JTWPA pump tone adjusts the frequencies of these peaks and, therefore, optimizes the gain at a particular frequency. A demonstration of this effect is shown in the plot of the receiver gain as a function of frequency, as shown in Fig. 3. In-situ measurements of the JTWPA gain were acquired by reversing the warm RF switch in Fig. 1 toward the VNA. Likewise, power measurements were made by reversing it toward the warm receiver. The SNRI associated with a given spectrum was computed by smoothing and interpolating between measurements of the JTWPA SNRI before and after sampling power from the cavity to acquire a reasonable value over the time duration of data acquisition.

We measured the attenuation of all cables and components between the cavity and the JTWPA in liquid nitrogen. We then scaled the measured attenuation to 100 mK, the temperature stage of the quantum amplifier package, using the temperature scaling ratios described in the cable specifications sheet.46 The value of ɛ, which we measured to be 0.40 ± 0.04 in linear units, was then used to compute the system noise. This measurement includes attenuation from the JTWPA, which has an associated insertion loss of −3.0 ± 0.3 dB47 at 4.798 GHz. Assuming an average JTWPA SNRI of 8.25 dB, the average system noise was measured to be 1.38 ± 0.08 K, with the dominant uncertainty coming from the uncertainty in the attenuation from the cavity to the JTWPA. This value is within the expected range of the system noise using a JTWPA.48 At a frequency of about 4.798 GHz, the standard quantum limit is 226 mK. The measured total system noise is therefore about 6 times the standard quantum limit. The JTWPA contributes about 1.5 ± 1.7 noise photons in this setup.33 

To acquire axion search data, the output line in Fig. 1 is connected to the warm receiver and digitizer, while all input lines to the cavity are terminated. Power is sampled for 10 ms on the output line 104 times. Each dataset is Fourier transformed and combined to form a single spectrum representative of a 100 s integration time. Each 2-MHz-wide spectrum consisted of 20 000 bins, each 100 Hz wide. The first and last 5000 bins were removed from all spectra because the cavity Lorentzian shape was distorted by intrinsic distortions in the gain of the JTWPA. The bins at the edge of the Lorentzian are least valuable to the data as they exist far from resonance. In a truly broadband experiment that was not limited by the bandwidth of the cavity and receiver electronics, these distortions could be removed, enabling a wide frequency range to be probed. Individual power spectra from the sidecar were processed by first dividing by the measured warm electronics transfer function and then applying a Padé approximant filter to remove the transfer function of the cold receiver chain.49 To acquire the warm electronics transfer function, we terminated the input of the warm receiver, then sampled the power from the warm receiver. An example of the Padé approximant fit (accounting for the warm receiver shape) to the raw spectrum is shown in Fig. 6. Removing the transfer functions of the warm and cold electronics resulted in power spectra having the expected Gaussian noise distribution, as shown in Fig. 7. Each bin in the spectrum was weighted according to the cavity Lorentzian given by a fit to the network analyzer transmission data, acquired by injecting power on the weak port and examining the output on the signal line (see Fig. 1). The transmission power as a function of frequency acquires the Lorentzian shape given by
(6)
where f0 is the resonant frequency, Q is the quality factor, f is the frequency, C is a constant offset, and δy is the depth of the Lorentzian. Individual spectra were then combined to form a grand spectrum in the process described in Ref. 1. A filter based on the Maxwell–Boltzmann line shape was applied to search for axion candidates. Such a candidate is expected to have a linewidth of about 15 kHz in this frequency range.
FIG. 6.

Padé approximant fit an individual raw spectrum for the sidecar. A strong peak in the spectrum was determined to be interference coming from IF electronics, as it appeared always in the same bins, even in power spectra acquired far off cavity resonance. It was removed by masking bins 2800–4000, explaining the gap in the data, above.

FIG. 6.

Padé approximant fit an individual raw spectrum for the sidecar. A strong peak in the spectrum was determined to be interference coming from IF electronics, as it appeared always in the same bins, even in power spectra acquired far off cavity resonance. It was removed by masking bins 2800–4000, explaining the gap in the data, above.

Close modal
FIG. 7.

Noise distribution (blue) and Gaussian fit (red) of the data after removing the receiver shape for the sidecar. After background subtraction, the distribution of powers measured is well represented by a Gaussian, as expected for white noise. The mean of the fit was computed to be −0.007, with a sigma of 0.990.

FIG. 7.

Noise distribution (blue) and Gaussian fit (red) of the data after removing the receiver shape for the sidecar. After background subtraction, the distribution of powers measured is well represented by a Gaussian, as expected for white noise. The mean of the fit was computed to be −0.007, with a sigma of 0.990.

Close modal

A strong peak between 4797.5 and 4797.75 MHz was determined to be due to radio-frequency interference by examining scans far from cavity resonance and observing that the peak occupied the same IF bins of the power spectrum regardless of the cavity frequency. A mask was created to exclude these frequencies from the final limit plot. The total bandwidth that was excluded from any individual raw spectrum due to the mask was 120 kHz.

Axion candidates were defined as any power excesses greater than a 3σ fluctuation above the average power level from the flattened spectrum. Two candidates were flagged but were not identified as axion-like due to a lack of persistence between scans. The systematic uncertainties are shown in Table I. The dominant systematic uncertainties come from the quality factor. The sensitivity achieved by this experiment can be seen in Fig. 8, which excludes axion dark matter with an assumption that axions constitute all of the dark matter, taken to have a density of 0.45 Ge V/cm3.50 The sensitivity can be seen in the larger context of axion searches in Fig. 9. Optimal filtering was used with the assumption of a Maxwell–Boltzmann line shape.51 The shape of the excluded region shown in Fig. 8 arises from the fact that the frequency of the TM010 mode shifted during the course of data-taking. These frequency shifts coincided with a filling of the liquid helium reservoir, which is known to induce mechanical vibrations in the insert. The mechanical vibrations cause small shifts in the positions of the antenna and tuning rod. These vibrations resulted in slight shifts in the resonant frequency of the cavity, despite the fact that the cavity was not intentionally tuned during the period of this data acquisition. The resonant frequency of the cavity drifted from about 4.7975 to 4.7990 GHz throughout the course of data-taking.

TABLE I.

Dominant sources of systematic uncertainty. The uncertainties were added in quadrature to attain the uncertainty on the total axion power from the cavity, shown in the bottom row. For the first entry, B is the magnetic field, V is the volume, and C010 is the form factor. The last row shows the total uncertainty of the axion power from the cavity.

SourceFractional
uncertainty
B2VC010 0.15 
Q 0.2 
Antenna coupling 0.01 
THFET/ɛ 0.11 
SNRI measurement 0.11 
Total on power 0.26 
SourceFractional
uncertainty
B2VC010 0.15 
Q 0.2 
Antenna coupling 0.01 
THFET/ɛ 0.11 
SNRI measurement 0.11 
Total on power 0.26 
FIG. 8.

Sensitivity plot from the sidecar data with a JTWPA, assuming an axion density of 0.45 GeV/cm3, equal to the total dark matter density at the Earth. The two lines represent DFSZ and KSVZ models. The prior limit from the CAST experiment is also shown.52 

FIG. 8.

Sensitivity plot from the sidecar data with a JTWPA, assuming an axion density of 0.45 GeV/cm3, equal to the total dark matter density at the Earth. The two lines represent DFSZ and KSVZ models. The prior limit from the CAST experiment is also shown.52 

Close modal
FIG. 9.

Sensitivity plot in the context of nearby searches. Green shows prior exclusions from the ADMX Collaboration. Red shows this work.

FIG. 9.

Sensitivity plot in the context of nearby searches. Green shows prior exclusions from the ADMX Collaboration. Red shows this work.

Close modal

In conclusion, we have demonstrated the sustained operation of a JTWPA for axion searches. While the data analyzed in this run were acquired in about 2 weeks, the JTWPA provided reasonable gain for a period of several months during engineering runs prior to data acquisition. Furthermore, an incident occurred in which an accidental increase in the magnet current led to a magnetic field of about 0.41 T in the vicinity of the JTWPA. As a result, the JTWPA was temporarily disabled, but its operation was fully restored after the insert temperature was raised above the critical temperature of niobium. We note that a similar recovery is expected for both the MSA and JPA. Such resilience against inadvertently applied fields bodes well for future axion experiments. This work sets the stage for broadband axion searches, in which we take full advantage of the ability of the JTWPA to operate with high gain over a wider bandwidth. The wide bandwidth of the JTWPA may be useful in frequency-multiplexed axion searches.53 It is clear that improvements to the noise performance can be made, and studies of JTWPA performance in small magnetic fields are warranted. Future searches that aim to use a JTWPA should continue to improve the means by which we measure the system noise temperature, and consider including an RF line to bypass the JTWPA. This bypass would provide some ability to isolate the JTWPA and, therefore, study its impact on the receiver chain. Finally, increasing the gain of the JTWPA would allow for a lower system noise temperature. It remains true that a big advantage of JPA over JTWPA is its near quantum-limited noise performance, which is valuable in extremely sensitive measurements such as the axion search. Improvements to the gain may be possible if magnetic field flux penetration into the amplifier package could be further reduced, and the gain optimization script improved. Although the JTWPA noise performance is less optimal than that of other quantum amplifiers, it remains true that it achieves wideband sensitivity that could be used in future axion search applications.

This work was supported by the U.S. Department of Energy through Grants Nos. DE-SC0009800, DESC0009723, DE-SC0010296, DE-SC0010280, DE-SC0011665, DEFG02-97ER41029, DEFG02-96ER40956, DEAC52-07NA27344, DEC03-76SF00098, and DE-SC0017987. Fermilab is a U.S. Department of Energy, Office of Science, HEP User Facility. Fermilab is managed by Fermi Research Alliance, LLC (FRA), acting under Contract No. DE-AC02-07CH11359. Additional support was provided by the Heising-Simons Foundation and the Lawrence Livermore National Laboratory and Pacific Northwest National Laboratory LDRD offices (LLNL Release No. LLNL-JRNL-825283). UWA was funded by the ARC Center of Excellence for Engineered Quantum Systems (Grant No. CE170100009) and Dark Matter Particle Physics (Grant No. CE200100008). Ben McAllister was funded by the Forrest Research Foundation. Tatsumi Nitta was supported by JSPS Overseas Research Fellowships No. 202060305. MIT Lincoln Laboratory acknowledges support from the Office of the Director of National Intelligence (ODNI), Intelligence Advanced Research Projects Activity (IARPA), under Air Force Contract No. FA8721-05-C-0002. The views and conclusions contained herein are those of the authors and should not be interpreted as necessarily representing the official policies or endorsements, either expressed or implied, of ODNI, IARPA, or the U.S. Government.

The authors have no conflicts to disclose.

C. Bartram: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Software (equal); Supervision (equal); Writing – original draft (equal); Writing – review & editing (equal). T. Braine: Data curation (equal); Investigation (equal); Writing – review & editing (equal). R. Cervantes: Conceptualization (equal); Data curation (equal); Investigation (equal); Software (equal); Writing – review & editing (equal). N. Crisosto: Data curation (equal); Investigation (equal); Writing – review & editing (equal). N. Du: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Software (equal); Writing – review & editing (equal). G. Leum: Data curation (equal); Investigation (equal). P. Mohapatra: Data curation (equal); Investigation (equal). T. Nitta: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Investigation (equal); Methodology (equal); Software (equal); Writing – review & editing (equal). L. J. Rosenberg: Conceptualization (equal); Data curation (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Supervision (equal). G. Rybka: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Writing – review & editing (equal). J. Yang: Data curation (equal); Investigation (equal); Methodology (equal); Writing – review & editing (equal). John Clarke: Writing – review & editing (equal). I. Siddiqi: Resources (equal). A. Agrawal: Writing – review & editing (equal). A. V. Dixit: Writing – review & editing (equal). M. H. Awida: Writing – review & editing (equal). A. S. Chou: Formal analysis (equal); Writing – review & editing (equal). M. Hollister: Writing – review & editing (equal). S. Knirck: Writing – review & editing (equal). A. Sonnenschein: Project administration (equal); Resources (equal); Writing – review & editing (equal). W. Wester: Writing – review & editing (equal). J. R. Gleason: Conceptualization (equal); Writing – review & editing (equal). A. T. Hipp: Writing – review & editing (equal). S. Jois: Writing – review & editing (equal). P. Sikivie: Conceptualization (equal); Methodology (equal); Writing – review & editing (equal). N. S. Sullivan: Methodology (equal); Writing – review & editing (equal). D. B. Tanner: Conceptualization (equal); Investigation (equal); Methodology (equal); Writing – review & editing (equal). E. Lentz: Formal analysis (equal); Writing – review & editing (equal). R. Khatiwada: Investigation (equal); Writing – review & editing (equal). C. Cisneros: Writing – review & editing (equal). N. Robertson: Writing – review & editing (equal). N. Woollett: Conceptualization (equal); Data curation (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Resources (equal); Software (equal); Writing – review & editing (equal). L. D. Duffy: Writing – review & editing (equal). C. Boutan: Conceptualization (equal); Investigation (equal); Methodology (equal); Resources (equal); Writing – review & editing (equal). M. Jones: Writing – review & editing (equal). B. H. LaRoque: Conceptualization (equal); Data curation (equal); Methodology (equal); Software (equal); Writing – review & editing (equal). N. S. Oblath: Software (equal); Writing – review & editing (equal). M. S. Taubman: Writing – review & editing (equal). E. J. Daw: Writing – review & editing (equal). M. G. Perry: Writing – review & editing (equal). J. H. Buckley: Writing – review & editing (equal). C. Gaikwad: Writing – review & editing (equal). J. Hoffman: Conceptualization (equal); Investigation (equal); Resources (equal); Validation (equal); Writing – review & editing (equal). K. Murch: Writing – review & editing (equal). M. Goryachev: Writing – review & editing (equal). B. T. McAllister: Writing – review & editing (equal). A. Quiskamp: Writing – review & editing (equal). C. Thomson: Investigation (equal); Writing – review & editing (equal). M. E. Tobar(ADMX Collaboration): Writing – review & editing (equal). V. Bolkhovsky: Resources (equal); Writing – review & editing (equal). G. Calusine: Resources (equal); Writing – review & editing (equal). W. Oliver: Resources (equal); Writing – review & editing (equal). K. Serniak: Resources (equal); Writing – review & editing (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

1.
C.
Bartram
,
T.
Braine
,
R.
Cervantes
,
N.
Crisosto
,
N.
Du
,
G.
Leum
,
L.
Rosenberg
,
G.
Rybka
,
J.
Yang
,
D.
Bowring
et al, “
Axion dark matter experiment: Run 1b analysis details
,”
Phys. Rev. D
103
,
032002
(
2021
).
2.
P. A. R.
Ade
et al, (
Planck), “Planck 2013 results. XVI. Cosmological parameters
,”
Astron. Astrophys.
571
,
A16
(
2014
).
3.
V. C.
Rubin
,
W. K.
Ford
, Jr.
, and
N.
Thonnard
, “
Rotational properties of 21 SC galaxies with a large range of luminosities and radii, from NGC 4605 (R=4kpc) to UGC 2885 (R=122kpc)
,”
Astrophys. J.
238
,
471
487
(
1980
).
4.
N.
Kaiser
and
G.
Squires
, “
Mapping the dark matter with weak gravitational lensing
,”
Astrophys. J.
404
,
441
450
(
1993
).
5.
A.
Robertson
,
R.
Massey
, and
V.
Eke
, “
What does the bullet cluster tell us about self-interacting dark matter?
,”
Monthly Notices of the Royal Astronomical Society
,
stw2670
(
2016
).
6.
J.
Preskill
,
M. B.
Wise
, and
F.
Wilczek
, “
Cosmology of the invisible axion
,”
Phys. Lett. B
120
,
127
132
(
1983
).
7.
L. F.
Abbott
and
P.
Sikivie
, “
A cosmological bound on the invisible axion
,”
Phys. Lett. B
120
,
133
136
(
1983
).
8.
M.
Dine
and
W.
Fischler
, “
The not-so-harmless axion
,”
Phys. Lett. B
120
,
137
141
(
1983
).
9.
R. D.
Peccei
and
H. R.
Quinn
, “
CP conservation in the presence of pseudoparticles
,”
Phys. Rev. Lett.
38
,
1440
1443
(
1977
).
10.
R. D.
Peccei
and
H. R.
Quinn
, “
Constraints imposed by CP conservation in the presence of pseudoparticles
,”
Phys. Rev. D
16
,
1791
1797
(
1977
).
11.
S.
Weinberg
, “
A new light boson?
Phys. Rev. Lett.
40
,
223
226
(
1978
).
12.
F.
Wilczek
, “
Problem of strong p and t invariance in the presence of instantons
,”
Phys. Rev. Lett.
40
,
279
282
(
1978
).
13.
P.
Sikivie
, “
Detection rates for ‘invisible’-axion searches
,”
Phys. Rev. D
32
,
2988
2991
(
1985
).
14.
P.
Sikivie
, “
Experimental tests of the invisible axion
,”
Phys. Rev. Lett.
51
,
1415
1417
(
1983
).
15.
M.
Dine
,
W.
Fischler
, and
M.
Srednicki
, “
A simple solution to the strong CP problem with a harmless axion
,”
Phys. Lett. B
104
,
199
(
1981
).
16.
A. P.
Zhitnitskii
, “
Possible suppression of axion-hadron interactions
,”
Sov. J. Nucl. Phys.
31
,
2
(
1980
) (English Translation).
17.
J. E.
Kim
, “
Weak interaction singlet and strong CP invariance
,”
Phys. Rev. Lett.
43
,
103
(
1979
).
18.
M. A.
Shifman
,
A. I.
Vainshtein
, and
V. I.
Zakharov
, “
Can confinement ensure natural CP invariance of strong interactions?
,”
Nucl. Phys. B
166
,
493
(
1980
).
19.
L.
Zhong
,
S.
Al Kenany
,
K. M.
Backes
,
B. M.
Brubaker
,
S. B.
Cahn
,
G.
Carosi
,
Y. V.
Gurevich
,
W. F.
Kindel
,
S. K.
Lamoreaux
,
K. W.
Lehnert
,
S. M.
Lewis
,
M.
Malnou
,
R. H.
Maruyama
,
D. A.
Palken
,
N. M.
Rapidis
,
J. R.
Root
,
M.
Simanovskaia
,
T. M.
Shokair
,
D. H.
Speller
,
I.
Urdinaran
, and
K. A.
van Bibber
, “
Results from phase 1 of the haystac microwave cavity axion experiment
,”
Phys. Rev. D
97
,
092001
(
2018
).
20.
B. M.
Brubaker
,
L.
Zhong
,
S. K.
Lamoreaux
,
K. W.
Lehnert
, and
K. A.
van Bibber
, “
Haystac axion search analysis procedure
,”
Phys. Rev. D
96
,
123008
(
2017
).
21.
B. M.
Brubaker
,
L.
Zhong
,
Y. V.
Gurevich
,
S. B.
Cahn
,
S. K.
Lamoreaux
,
M.
Simanovskaia
,
J. R.
Root
,
S. M.
Lewis
,
S.
Al Kenany
,
K. M.
Backes
,
I.
Urdinaran
,
N. M.
Rapidis
,
T. M.
Shokair
,
K. A.
van Bibber
,
D. A.
Palken
,
M.
Malnou
,
W. F.
Kindel
,
M. A.
Anil
,
K. W.
Lehnert
, and
G.
Carosi
, “
First results from a microwave cavity axion search at 24 μeV
,”
Phys. Rev. Lett.
118
,
061302
(
2017
).
22.
C.
Hagmann
,
D.
Kinion
,
W.
Stoeffl
,
K.
van Bibber
,
E.
Daw
,
H.
Peng
,
L. J.
Rosenberg
,
J.
LaVeigne
,
P.
Sikivie
,
N. S.
Sullivan
,
D. B.
Tanner
,
F.
Nezrick
,
M. S.
Turner
,
D. M.
Moltz
,
J.
Powell
, and
N. A.
Golubev
, “
Results from a high-sensitivity search for cosmic axions
,”
Phys. Rev. Lett.
80
,
2043
2046
(
1998
).
23.
S. J.
Asztalos
,
R. F.
Bradley
,
L.
Duffy
,
C.
Hagmann
,
D.
Kinion
,
D. M.
Moltz
,
L. J.
Rosenberg
,
P.
Sikivie
,
W.
Stoeffl
,
N. S.
Sullivan
,
D. B.
Tanner
,
K.
van Bibber
, and
D. B.
Yu
, “
Improved rf cavity search for halo axions
,”
Phys. Rev. D
69
,
011101(R)
(
2004
).
24.
S. J.
Asztalos
,
E.
Daw
,
H.
Peng
,
L. J.
Rosenberg
,
D. B.
Yu
,
C.
Hagmann
,
D.
Kinion
,
W.
Stoeffl
,
K.
van Bibber
,
J.
LaVeigne
,
P.
Sikivie
,
N. S.
Sullivan
,
D. B.
Tanner
,
F.
Nezrick
, and
D. M.
Moltz
, “
Experimental constraints on the axion dark matter halo density
,”
Astrophys. J. Lett.
571
,
L27
(
2002
).
25.
O.
Kwon
,
D.
Lee
,
W.
Chung
,
D.
Ahn
,
H.
Byun
,
F.
Caspers
,
H.
Choi
,
J.
Choi
,
Y.
Chong
,
H.
Jeong
et al, “
First results from an axion haloscope at capp around 10.7 μ ev
,”
Phys. Rev. Lett.
126
,
191802
(
2021
).
26.
J.
Choi
,
S.
Ahn
,
B. R.
Ko
,
S.
Lee
, and
Y. K.
Semertzidis
, “
Capp-8tb: Axion dark matter search experiment around 6.7 μev
,”
Nucl. Instrum. Methods Phys. Res., Sect. A
1013
,
165667
(
2021
).
27.
D.
Alesini
,
C.
Braggio
,
G.
Carugno
,
N.
Crescini
,
D.
D’Agostino
,
D.
Di Gioacchino
,
R.
Di Vora
,
P.
Falferi
,
U.
Gambardella
,
C.
Gatti
,
G.
Iannone
,
C.
Ligi
,
A.
Lombardi
,
G.
Maccarrone
,
A.
Ortolan
,
R.
Pengo
,
A.
Rettaroli
,
G.
Ruoso
,
L.
Taffarello
, and
S.
Tocci
, “
Search for invisible axion dark matter of mass ma = 43 μeV with the quax– experiment
,”
Phys. Rev. D
103
,
102004
(
2021
).
28.
M.
Mück
,
M.-O.
André
,
J.
Clarke
,
J.
Gail
, and
C.
Heiden
, “
Radio-frequency amplifier based on a niobium dc superconducting quantum interference device with microstrip input coupling
,”
Appl. Phys. Lett.
72
,
2885
2887
(
1998
).
29.
M.
Hatridge
,
R.
Vijay
,
D. H.
Slichter
,
J.
Clarke
, and
I.
Siddiqi
, “
Dispersive magnetometry with a quantum limited squid parametric amplifier
,”
Phys. Rev. B
83
,
134501
(
2011
).
30.
S.
Boutin
,
D. M.
Toyli
,
A. V.
Venkatramani
,
A. W.
Eddins
,
I.
Siddiqi
, and
A.
Blais
, “
Effect of higher-order nonlinearities on amplification and squeezing in josephson parametric amplifiers
,”
Phys. Rev. Applied
8
,
054030
(
2017
).
31.
T.
Roy
,
S.
Kundu
,
M.
Chand
,
A. M.
Vadiraj
,
A.
Ranadive
,
N.
Nehra
,
M. P.
Patankar
,
J.
Aumentado
,
A. A.
Clerk
, and
R.
Vijay
, “
Broadband parametric amplification with impedance engineering: Beyond the gain-bandwidth product
,”
Appl. Phys. Lett.
107
,
262601
(
2015
).
32.
J.
Aumentado
, “
Superconducting parametric amplifiers: The state of the art in josephson parametric amplifiers
,”
IEEE Microwave Magazine
21
,
45
59
(
2020
).
33.
C.
Braggio
,
G.
Cappelli
,
G.
Carugno
,
N.
Crescini
,
R.
Di Vora
,
M.
Esposito
,
A.
Ortolan
,
L.
Planat
,
A.
Ranadive
,
N.
Roch
et al, “
A haloscope amplification chain based on a traveling wave parametric amplifier
,”
Rev. Sci. Instrum.
93
,
094701
(
2022
).
34.
H.
Johnston
, “
Us invests $249 m in quantum tech
,”
Phys. World
31
,
9
(
2018
).
35.
A.
Arvanitaki
,
S.
Dimopoulos
,
M.
Galanis
,
L.
Lehner
,
J. O.
Thompson
, and
K.
Van Tilburg
, “
Large-misalignment mechanism for the formation of compact axion structures: Signatures from the qcd axion to fuzzy dark matter
,”
Phys. Rev. D
101
,
083014
(
2020
).
36.
L. J.
Rosenberg
, “
Final report for the admx phase 2a project at the university of Washington
,” Tech. Rep. (
Univ. of Washington
,
Seattle, WA
,
2015
).
37.
C.
Boutan
,
A piezoelectrically tuned RF-cavity search for dark matter axions
, Ph.D. thesis,
University of Washington
(
2017
).
39.
COMSOL AB, Stockholm, Sweden, “Comsol multiphysics®.”
40.
“Low noise factory hfet amplifier,” https://www.lownoisefactory.com/files/7015/7825/6000/LNF-LNC4_8C.pdf, accessed: 2021-09-23.
41.
K.
O’Brien
,
C.
Macklin
,
I.
Siddiqi
, and
X.
Zhang
, “
Resonant phase matching of josephson junction traveling wave parametric amplifiers
,”
Phys. Rev. Lett.
113
,
157001
(
2014
).
42.
H. A.
Haus
and
J. A.
Mullen
, “
Quantum noise in linear amplifiers
,”
Phys. Rev.
128
,
2407
2413
(
1962
).
43.
T.
Braine
et al (
ADMX Collaboration
), “
Extended search for the invisible axion with the axion dark matter experiment
,”
Phys. Rev. Lett.
124
,
101303
(
2020
).
44.
R.
Khatiwada
,
D.
Bowring
,
A. S.
Chou
,
A.
Sonnenschein
,
W.
Wester
,
D. V.
Mitchell
,
T.
Braine
,
C.
Bartram
,
R.
Cervantes
,
N.
Crisosto
et al, “
Axion dark matter experiment: Detailed design and operations
,”
Rev. Sci. Instrum.
92
,
124502
(
2021
).
45.
E.
Daw
and
R. F.
Bradley
, “
Effect of high magnetic fields on the noise temperature of a heterostructure field-effect transistor low-noise amplifier
,”
J. Appl. Phys.
82
,
1925
1929
(
1997
).
46.
“Keycom,” https://www.keycom.co.jp/eproducts/upj/upj2/page.htm, accessed: 2021-09-26.
47.
C.
Macklin
,
K.
O’Brien
,
D.
Hover
,
M. E.
Schwartz
,
V.
Bolkhovsky
,
X.
Zhang
,
W. D.
Oliver
, and
I.
Siddiqi
, “
A near–quantum-limited josephson traveling-wave parametric amplifier
,”
Science
350
,
307
310
(
2015
).
48.
S.
Simbierowicz
,
V.
Vesterinen
,
J.
Milem
,
A.
Lintunen
,
M.
Oksanen
,
L.
Roschier
,
L.
Grönberg
,
J.
Hassel
,
D.
Gunnarsson
, and
R. E.
Lake
, “
Characterizing cryogenic amplifiers with a matched temperature-variable noise source
,”
Rev. Sci. Instrum.
92
,
034708
(
2021
).
49.
E. J.
Daw
, “
A search for halo axions
,” Ph.D. thesis (
MIT
,
2022
).
50.
M. S.
Turner
, “
Cosmic and local mass density of ‘invisible’ axions
,”
Phys. Rev. D
33
,
889
896
(
1986
).
51.
M. S.
Turner
, “
Periodic signatures for the detection of cosmic axions
,”
Phys. Rev. D
42
,
3572
(
1990
).
52.
C.
Collaboration
et al, “
New cast limit on the axion–photon interaction
,”
Nature Physics
13
,
584
590
(
2017
).
53.
See https://www.snowmass21.org/ for “Snowmass 2021-letter of interest frequency multiplexed dark matter axionsearches,” by
G.
Carosi
.