New method of galactic axion search

An appealing candidate of the galactic dark matter is the axion, which was postulated to solve the strong CP (Charge-conjugation Parity) violation problem in the standard particle theory. A new experimental method is proposed to determine the axion mass. The method uses collectively and coherently excited atoms or molecules, the trigger laser inducing galactic axion absorption along with signal photon emission to be detected.


Introduction
The galactic dark matter problem may be solved by the QCD (Quantum Chromo-Dynamics) axion [1][2][3][4][5][6][7], which provides simultaneously a solution to the conundrum of the strong CP problem [8][9][10][11] in the standard particle theory. The necessary galactic dark matter mass density is ρ G (0.3-0.45) GeV/c 2 cm −3 which implies by far the largest ambient number density n a ∼ 10 13 cm −3 (its precise value depending on the Peceei-Quinn (PQ) symmetry breaking scale f a ). The current parameter region explored in major experiments is in the range, f a = 10 8 − 10 13 GeV, the corresponding axion mass in the range, m a = 10 2 − 10 −3 meV (inversely proportional to f a ).
The cosmological origin of the galactic axion is traced back to the QCD epoch of cosmic temperature 100 MeV/k B ∼ 10 12 K, since the axion mass is generated by the QCD chiral symmetry breaking. Its couplings to ordinary quarks, leptons, and gauge bosons are all suppressed by the PQ symmetry breaking scale ∝ 1/ f a , hence axions are born in cold, namely their initial velocities at creation are much less than the velocity of light. Cold axions decouple from radiation and matter throughout the later epoch of cosmological evolution [5][6][7]. After the radiation energy density drops below the axion matter density, gravitational self-interaction begins a e-mail: sasao@okayama-u.ac.jp b e-mail: yoshim@okayama-u.ac.jp to dominate and axions become virialized (process of violent relaxation [12]), leaving a velocity dispersion of order 10 −3 × the light velocity.
The PQ-symmetry breaking scale may be raised if it is higher than the reheat temperature due to inflation [13]. This widens the parameter space to be explored up to f a ∼ 10 16 GeV.
Ongoing and proposed experiments [14][15][16][17], of galactic axion search sensitive to low mass axions use magnetic axion conversion into microwave in cavity [18,19]. In this paper we propose an entirely new experimental method of axion detection using atoms or molecules, which is sensitive to the axion mass. It provides a powerful tool to directly link physics beyond the standard theory and the dark matter cosmology/astrophysics.
We use the natural unit ofh = c = 1 unless otherwise stated.

Microscopic process of galactic axion detection
Our method proposed here uses atomic or molecular process, a G + γ t + |i → | f + γ s : a galactic axion (a G ) collides against an atom or a molecule in an excited state |i , which then de-excites to a lower energy state | f , emitting a signal photon γ s under the trigger photon γ t . The process is called TRiggered Absorption of Cosmic Axion (TRACA for short). Relevant diagrams are depicted in Fig. 1.
The virtual photon of four momentum k may be off or on the mass shell, depending on its invariant mass squared k 2 = k 2 0 − k 2 = 0 or 0. The two cases have different experimental merits. In the on-shell case one may view TRACA occurring in two steps of real processes: the triggered decay axion produces a real photon, which is inelastically scattered off excited atom/molecule to generate a signal photon γ s . The on-shell TRACA may occur in the trigger photon frequency ω t at m a /2, hence in the rf (radio frequency) or microwave region, while the trigger frequency of the off-shell TRACA may be in the optical to the infrared (IR) range depending on the level spacing i f .
We also considered TRACA arising from diagrams containing direct electron coupling with the axion (both trigger laser and axion are connected to atomic line directly). In most axion parameter region we found that rates of axion emission are smaller than those of Fig. 1 even for DFSZ model [3,4] of large axion-electron coupling.
The method can be applied to axion-like particles [20] as well.
In laboratory experiments on Earth, axions collide with its momentum around m a β a = − m a β e , where β a (β e ) is the velocity of axion (Earth). The axion momentum distribution function F a (q) is assumed to have a momentum dispersion ), allowing for a possible BEC component [21], this component being separately detectable. The uniform − β e wind is superimposed by diurnal and sidereal modulation, which can be detected, too. One of our ultimate objectives of the proposed experiment is to determine these important parameters of cosmology.
To maximize TRACA rates, we fully utilize effects of large axion occupation and coherence (denoted by ρ i f ) between the initial and final states: the total TRACA rate for entire atoms is proportional to the product of three factors, (1) axion number density n a , giving the enhancement ), (2) field power |E t | 2 = ω t n t (its photon number density denoted by n t ) driving the axion decay, (3) macro-coherently prepared atomic or molecular target number squared N 2 T ρ 2 i f /V with N T the total number of excited atoms within a volume V . The principle of amplification by macro-coherence [22,23] related to the second and the third issues has been experimentally demonstrated in [24][25][26] for a weak QED (Quantum ElectroDynamics) process of two-photon emission in para-hydrogen (p-H 2 ) vibrational transition (the quantum number change of v = 1 → 0). The process was amplified by a factor ∼ 10 18 in rate, with a realization of atomic coherence ρ i f ∼ 0.03 [26]. The degree of achievable coherence was limited by various de-coherence effects including Doppler and pressure broadenings as well as laser line width. The process was termed PSR (Paired Super-Radiance) in an obvious analogy to Dicke's super-radiance [27][28][29][30].

Derivation of rate formula
We shall first treat the off mass shell case and later discuss briefly the on mass shell case. Suppose that the intermediate atomic state |n related by electric dipole transitions to |i and | f is far away in energy from these states; both of ni , n f i f . This holds for the example of p-H 2 . We fur- Fig. 1 TRACA diagram of axion detection. Virtual or real photon emitted by the axion a G decay due to trigger γ t induces inelastic photon or Raman scattering producing a signal photon γ s . Diagram of exchanged atomic/molecular vertexes (s-channel |n here being replaced by uchannel |n not depicted here) adds coherently to this contribution. This diagram is further doubled by γ t ↔ γ s exchange. The quantities in (··) denote the 4-vector symbols used in the text ther assume the frequency relation, both of ω t , ω s (trigger and signal frequencies) m a for the off-shell TRACA. The axion-photon-photon vertex in Fig. 1 has a form g aγ γ aE · B in terms of electric E and magnetic B field operators, since the axion field a is a pseudo-scalar. Its coupling constant is suppressed by g aγ γ α/(π f a ) , α = 1/137. The dominant atomic couping along the electron line is of the usual electric dipole form d · E , giving rise to an effective, external field coupling of the form, B t · 0|T (E E)|0 · d with propagation taken into account (B t : the trigger photon magnetic field). The central part of the probability amplitude is thus Here d n f , d ni are dipole operators and their product divided by the energy denominator is related to i f off-diagonal polarizabiity when summed over states |n . The propagator factor is approximately (k t − q) 2 0 /(k t − q) 2 ∼ −ω t /(2m a ) with q the 4-vector of axion. Model dependent coefficient c aγ γ [31][32][33] is typically −0.97 for KSVZ-model [1,2], and 0.36 for DFSZ-model [3,4].
For simplicity we take the atomic transition of no spin parity change, 0 + → 0 + . Excluding the squared probability amplitude given above, the rate contains a product of coherence and position-dependent phase factor of the target system integrated over a target volume V prior to taking the square, with Q = k t + k s − p i f and p i f the phase imprinted at the excitation to |i [34]. The coherence ρ i f and the excited target number density n T = N T /V may in principle be a function of target atom positions, but we take its averaged constant values within the entire target. We perform the axion momentum integration first, and use the fact that q 0 = m a β 0 = 10 neV( m a 10μeV ) much smaller than the inverse of linear size of target volume. The infinite target volume limit then gives Eq. (1) ∼ ρ 2 i f n T N T (2π) 3 δ(k t + k s − p i f − q e ) , q e = −m a β e , which cannot be used however due to the following reason.
Due to kinematic reasons one needs to calculate TRACA rates in a finite target volume, for which we take a cylinder of radius R and length L along excitation lasers. This gives the geometric factor defined in Eq. (1), with J 1 (x) the Bessel function of the first order. In our parameter range the dimensionless A is of order 10 −7 when we integrate over all solid angles. The differential off-shell TRACA rate is thus calculated as μ i f the atomic or molecular polarizability. The relation E 2 t = ω t n t of the trigger field power to the photon number density n t and various level spacings, ab = a − b were used. We assumed that only one dominant upper level |n contributes, but extension to include many contributing levels is straightforward.
A striking feature of the off-shell TRACA formula (3) is the dependence sin 2 θ pol with the relative angle θ pol between polarizations of two photons, γ t and γ s . The rate is maximal at the perpendicular θ pol = π/2 and vanishes at the parallel configuration.
We shall make an order of magnitude estimate of on-shell rate. The easiest way of calculating the on-shell amplitude is to replace the real propagator of virtual photon by the imaginary part in the zero-width approximation, to give for In practice, the squared delta-function in the rate of on-shell amplitude should be convoluted with the trigger frequency distribution of a finite width . Moreover, (2π) 2 × the squared delta-function is replaced by 2π T × the delta-function where T is the time of trigger irradiation. This time T is equal to the target length/2.

TRACA rates for p-H 2 vibrational transition
A typical experimental scheme we consider uses two excitation lasers and a trigger field, all of which may be irradiated near the same axis direction. The target is of a cylindrical shape whose length and volume are denoted by L , V . Consider para-hydrogen molecule as a high density target, its gas density at the liquid nitrogen temperature 77 K being ∼ 5 × 10 19 cm −3 . Solid p-H 2 target is also promising from a number of reasons [25]; (1) a larger molecule density like 2.6 × 10 22 cm −3 , (2) expected small relaxation caused by phonon emission. Relevant phonon emission involves the axion-proton coupling predominantly of the spin operator. The spatial extension of axion wave function is of order 1 cm for 10 μeV mass, hence the axion field feels a collective body of para-hydrogen molecule. Since the molecule has a total nuclear spin 0, this coupling, hence the phonon relaxation, is highly suppressed. The vibrational transition of electronically ground state, v = 1 → 0 of level spacing 0.52 eV, may be used. The electronically excited level (1s)(2 p) is far away from the ground vibrational levels, separated by ∼ 10 eV.
To calculate the signal rate for molecules, we need to take into account vibrational modes associated with electronically excited states such as (1s)(2 p). The effect of Franck-Condon factor related to the molecular potential curve of this first electronically excited state may be automatically avoided by using p-H 2 off-diagonal polarizability μ i f ∼ 1.43 × 10 −24 cm 3 [35].
The off-shell TRACA rate for solid p-H 2 is numerically The on-shell TRACA rate is obtained by multiplying ∼ 2 × We illustrate total rates of the off-shell and the on-shell TRACA contributions in Fig. 2. There are higher sensitivities for smaller axion masses less than 100 μeV.

Raman excitation and trigger
We consider Raman-type of excitation from the ground v = 0 to v = 1 state, with frequencies, ω 1 > ω 2 , ω 1 − ω 2 = 0.52eV. In the on-shell process the trigger frequency must coincide with the axion mass in search (in the rf or microwave frequency range with photon density of n t = 10 18 cm −3 ), and the momentum conservation dictates the signal photon µ Fig. 2 TRACA rate for p-H 2 . Assumed parameters are the target number density n T = 2.6 × 10 22 cm −3 , the target volume V = 1 cm 3 , the target length L = 10 cm, the coherence ρ i f = 0.1, ω t = 0.124 eV trigger for the off-shell case and ω t = m a /2, = 300 MHz for the on-shell case and of a common photon number density n t = 10 18 cm −3 and a common A = 1. Experimental parameters such as detection efficiency or duty factor (if a pulsed laser is used) are not included. We used c aγ γ values for the KSVZ and DFSZ models, and the perpendicular polarization of the trigger and the signal was taken, using the formula (5) emission very near the excitation axis. On the other hand, in the off-shell case it can be in a near infra-red range. For example, using CO 2 laser of ω t = 0.124eV for the trigger, one can relate the angles of trigger irradiation and signal emission as functions of the axion mass: θ t = 11 mrad m a 10μeV and θ s = −3.5 mrad m a 10μeV (formulas valid for small axion mass less than of order 1 meV) away from both the trigger and the excitation lasers. The off-axis configuration and its kinematic constraint is illustrated in Fig. 3. This much of the off-axis configuration may be crucial in rejecting various QED backgrounds (see below on more of this).

Backgrounds
By far the largest background may arise from (macrocoherently amplified) PSR process. The ratio of the signal (off-shell part) to this PSR rate is found to be ∼  Fig(2). The last factor comes from the fact that polarization vectors of γ t and γ s must be parallel in the case of the PSR process. The geometric factor A would suppress the background (since A(m a ) A(m a = 0)), but not entirely due to diffuseness of the finite-size volume effects. Reduction of Fig. 3 Kinematic configuration of initial Raman excitation in the left, the final decay processes in the right. Both energy and momentum conservations are satisfied for PSR only when the trigger (ω t ) is injected colinearly with the excitation lasers, while they are satisfied for TRACA when the trigger (ω t ) is injected with angle due to the axion absorption. These different configurations help separate two processes the background is also possible by choosing θ pol π/2, but not enough with current technologies. It is crucial to devise a clever method of PSR background rejection for a realistic experimental proposal.

Summary
In summary, we proposed a novel experimental method of axion detection using macro-coherently excited atoms or molecules which can be implemented in small scale laboratories. The method may clarify how the galactic axion dark matter is distributed and serve to determine the fundamental symmetry breaking scale beyond the standard particle theory.