# Cosmic Neutrinos and their Detection

## I Introduction

According to the standard Big Bang model, a cosmic neutrino background, similar to the familiar microwave background, should exist in our universe. Because of their weak interactions, cosmic neutrinos fell out of thermal equilibrium at $`t1\mathrm{sec}`$, and have been redshifting since then. Neutrinos of mass $`10^3`$ eV would still be relativistic today, have a Fermi-Dirac spectrum with $`T1.9`$ K, and be spatially uniform with $`n_{\nu _L}=n_{\overline{\nu }_R}50/\mathrm{cm}^3`$ per flavor. In contrast, neutrinos of mass $`10^3`$ eV would be nonrelativistic today, and be clustered around galaxies with a typical velocity $`v_\nu 10^3`$. They contribute to the cosmological energy density an amount $`\mathrm{\Omega }_\nu _im_{\nu _i}/(92h^2\mathrm{eV})`$, where $`h0.65`$ is the Hubble expansion rate in units of $`100\mathrm{km}\mathrm{s}^1\mathrm{Mpc}^1`$, and where the sum includes all neutrino flavors that are nonrelativistic today. Massive neutrinos are a natural candidate for the hot component in currently favored “Mixed Hot+Cold Dark Matter (HCDM)”  models of galaxy formation. In this scenario, neutrinos would contribute $``$ 20 % ($`_im_{\nu _i}8\mathrm{eV}`$), and CDM (e.g. Wimps and axions) the remainder of the dark matter. Another popular model involving a cosmological constant is $`\mathrm{\Lambda }\mathrm{HCDM}`$  with $`\mathrm{\Omega }_\nu 0.050.1`$. The available phase space for massive neutrinos  restricts the local neutrino number density to

$$n_\nu <\mathrm{\hspace{0.17em}2}\times 10^6\mathrm{cm}^3\left(\frac{v_m}{10^3}\right)^3\underset{i}{}\left(\frac{m_{\nu _i}}{10\mathrm{eV}}\right)^3.$$
(1)   
where $`v_m`$ is the maximum neutrino velocity in the halo. Hence, neutrinos of $`m_\nu `$ eV, cannot account for all of the local halo matter density $`300\mathrm{MeV}/\mathrm{cm}^3`$.

The detection of cosmic neutrinos by interactions with single electrons or nucleons is unlikely because of the extremely small rates and energy deposits. Past proposals have therefore focused on detecting the coherent mechanical effects on macroscopic targets due to the “neutrino wind”. In 1974, Stodolsky  suggested to use the energy splitting of an electron,

$$\mathrm{\Delta }E=\sqrt{2}G_Fv_e\left((n_{\nu _e}n_{\overline{\nu }_e})(n_{\nu _\mu }n_{\overline{\nu }_\mu })(n_{\nu _\tau }n_{\overline{\nu }_\tau })\right)$$
(2)   
moving through the neutrino background with velocity $`v_e`$. However even for large neutrino number asymmetries, e.g. $`n_{\overline{\nu }_e}=0`$, $`n_{\nu _{\mu ,\tau }}=n_{\overline{\nu }_{\mu ,\tau }}`$, $`\mathrm{\Delta }E1.3\times 10^{33}\mathrm{eV}(v_e/10^3)(n_{\nu _e}/10^7\mathrm{cm}^3)`$ is still tiny. In principle, this effect could be measured by observing a torque on a permanent magnet, which is shielded against magnetic noise with superconductors. A benefit of this detection method is that it works equally well for Dirac and Majorana neutrinos. Other forces $`G_F`$ arising from the reflection or refraction of neutrinos by macroscopic objects have been proposed by several authors . However, all of these ideas were later found to be flawed .

Another approach is to consider forces $`G_F^2`$ due to random elastic neutrino scattering . Here, spatial coherence dramatically increases the cross section of targets smaller than the neutrino wavelength $`\lambda _\nu `$. In the nonrelativistic limit, one must distinguish between Majorana and Dirac neutrinos. For Dirac $`\mu `$ or $`\tau `$ neutrinos, the cross section is dominated by the vector neutral current contribution

$$\sigma _D\frac{G_F^2m_\nu ^2}{8\pi }N_n^2=2\times 10^{55}\left(\frac{m_\nu }{10\mathrm{eV}}\right)^2N_n^2\mathrm{cm}^2$$
(3)   
where $`N_n`$ is the number of neutrons in the target of size $`<\lambda _\nu /2\pi =20\mu \mathrm{m}(10\mathrm{eV}/m_\nu )(10^3/v_\nu )`$. For Majorana neutrinos the vector contribution to the cross section is suppressed by a factor $`(v_\nu /c)^2`$ and potentially the largest contribution arises from the axial current. The cross section of a spin-polarized target is $`\sigma _MN_s^2`$ , where $`N_s`$ is the number of aligned spins in the grain.

Due to the Sun’s peculiar motion $`(v_s220\mathrm{km}/\mathrm{s})`$ through the galaxy, a test body will experiences a neutrino wind force through random neutrino scattering events. The wind direction is modulated with the sidereal period of 23hrs+56min due to the combined effects of the Earth’s diurnal rotation and annual revolution around the Sun. For Dirac neutrinos, the acceleration of a target of density $`\rho `$ and optimal radius $`\lambda _\nu /2\pi `$ has an amplitude 

$$a=8\times 10^{24}\left(\frac{AZ}{A}\right)^2\left(\frac{v_\mathrm{s}}{10^3}\right)^2\left(\frac{n_\nu }{10^7\mathrm{cm}^3}\right)\left(\frac{\rho }{20\mathrm{gcm}^3}\right)\frac{\mathrm{cm}}{\mathrm{s}^2}$$
(4)   
and is independent of $`m_\nu `$. For clustered Majorana neutrinos, the acceleration is suppressed by a factor $`10^6(10^3)`$ for an unpolarized (polarized) target. For unclustered relativistic neutrinos, the force should be aligned with the direction of the microwave background dipole. Neutrinos of Dirac or Majorana type would have the same cross section, giving rise to an acceleration $`10^{34}\mathrm{cm}/\mathrm{s}^2`$. A target size much larger than $`\lambda _\nu `$ can be assembled, while avoiding destructive interference, by using foam-like  or laminated materials . Alternatively, grains of size $`\lambda _\nu `$ could be randomly embedded (with spacing $`\lambda _\nu `$) in a low density host material.

## II Detector Concept

The measurement of the neutrino-induced forces and torques will require major improvements in sensor technology. At present, the most sensitive detector of small forces is the “Cavendish-type” torsion balance (see Figure 1), which has been widely used for measurements of the gravitational constant and searches for new forces. A typical arrangement consists of a dumbbell shaped test mass suspended by a tungsten fiber inside a vacuum chamber at room temperature. The angular deflection is usually read out with an optical lever. The torsional oscillation frequency is constrained by the yield strength of the fiber to $`>`$ 1mHz. Internal friction in the fiber limits the $`Q`$ to $`<10^6`$. Nonlinearities in the suspension make the balance quite sensitive to seismic noise. Additional background forces arise from thermal noise and time varying gravity gradients due to tides, weather, traffic, people etc. The angular resolution of optical levers is $`10^{10}\mathrm{rad}`$ and the smallest measurable acceleration is $`10^{13}\mathrm{cm}/\mathrm{s}^2`$ .

Several improvements seem possible (see Figure 2): Thermal noise can be decreased by lowering the temperature and by employing a low dissipation (high $`Q`$) suspension, as seen from the expression for the effective thermal noise acceleration

$$a_{\mathrm{th}}2\times 10^{23}\left(\frac{T}{1\mathrm{K}}\right)^{1/2}\left(\frac{1\mathrm{k}\mathrm{g}}{m}\right)^{1/2}\left(\frac{1\mathrm{d}\mathrm{a}\mathrm{y}}{\tau _0}\right)^{1/2}\left(\frac{10^6\mathrm{s}}{\tau }\right)^{1/2}\left(\frac{10^{16}}{Q}\right)^{1/2}\frac{\mathrm{cm}}{\mathrm{s}^2}$$
(5)   
where $`m`$ is the target mass, $`\tau `$ is the measurement time, $`\tau _0`$ is the oscillator period, and $`T`$ is the operating temperature. A promising high-$`Q`$ suspension method is a Meissner suspension consisting of a superconducting body floating above a superconducting coil. Niobium or NbTi in bulk or film form have been used in the past in such diverse applications as gravimeters, gyros, and gravitational wave antennas. Generally the magnetic field applied to the superconductor is limited to $`B0.2\mathrm{T}`$ to avoid flux penetration or loss of superconductivity. However, even for small fields the magnetic flux exclusion is usually incomplete in polycrystalline Nb superconductors and flux creep will cause noise . The lifting pressure $`B^2`$ is limited to about $`100\mathrm{g}/\mathrm{cm}^2`$. This could be dramatically increased by replacing the bulk superconductor with a persistent-mode superconducting magnet, as shown in Figure 2. Commonly used NbTi wire has a critical field of several T. Moreover, the flux lines are strongly pinned by artificial lattice defects in the magnet wire, leading to low levels of flux noise and dissipation. Very low resistance $`10^{14}\mathrm{\Omega }`$  wire joints can be made by cold welding, and the magnetic field decay rate can be as low as $`\dot{B}/B=R/L10^9/\mathrm{day}`$, where $`R`$ is the joint resistance and $`L`$ is the inductance. The current decay due to losses in the joint will cause the suspended magnet to sink gradually. The concomitant change of magnetization will produce a slowly varying torque on the float due to the gyromagnetic (“Einstein-deHaas”) effect. Due to unavoidable deviations from cylindrical symmetry, the oscillator will have a very long, but finite rotational oscillation period. For small amplitudes, the oscillating supercurrents in the wire will cause the fluxoids to move elastically about their pinning centers, and very little dissipation is expected. In order to match the rotation period to the signal period $`(1\mathrm{day})`$, an additional restoring potential may be added. One possibility is a superconducting or low-loss dielectric (e.g. single-crystal sapphire) ellipsoid in the constant electric field of a parallel plate capacitor. Both negative and positive torsion coefficients can be produced by aligning the long axis of the ellipsoid perpendicular or parallel to the field.

The effects of time-varying gravity gradients can be reduced with a highly cylindrically symmetric target (see Figure 1). With the c.m. of the target centered below the suspension support, the leading order gravitational torque arise from the dipole moment of the torsion balance. For a source mass $`M`$ at a distance $`R`$ from the balance c.m., the torque is

$$\tau _z\frac{G_NM}{R^3}(I_xI_y)$$
(6)   
where $`I_x,I_y`$ are the moments of inertia about the two horizontal axes. Diurnally varying mass distributions, $`e.g.`$ changes of the atmospheric air mass due to solar heating, are especially worrisome. These effects can be minimized with a careful mechanical design and subsequent balancing steps using a movable laboratory source mass. In addition, locating the experiment in a deep mine would be beneficial. The gravitational interaction of the balance with the Sun and the Moon are less dangerous, because they produce torques with a 12 hr period to first order.

Low frequency seismic noise can also be a performance limiting factor in torsion balance experiments. At a period of $``$ 24 hrs, the noise is dominated by tidal effects, with typical variations in the strain of $`\mathrm{\Delta }L/L10^7`$ , in the tilt about the local (g-referenced) vertical axis of $`\mathrm{\Delta }\theta 10^6`$ rad, and in the local gravitational acceleration of $`\mathrm{\Delta }g/g10^7`$. Small deviations from ideal symmetry will couple these motions into the rotational mode of the balance. For example, the tilt-rotation coupling of a tungsten fiber suspension is typically $`0.02`$. Two possible remedies are (1) improving the symmetry of the balance, and (2) employing active anti-seismic isolation systems. A remaining worry is low frequency rotational seismic noise for which no reliable data is available. It will directly mask the signal and needs to be compensated.

The proposed angular rotation readout has very high immunity against vertical, horizontal, and tilt seismic noise. It consists of a parametric transducer, which converts the angle to an optical frequency. As shown in Figure 3, the transducer consists of a high-$`Q`$ optical cavity of length $`l`$, tuned by a Brewster-angled low loss dielectric plate of thickness $`d`$. Using high-reflectivity mirrors and a cavity length $`10\mathrm{cm}`$ should give linewidths of order 100 kHz for the Gaussian $`\mathrm{TEM}_{00p}`$ modes, with a frequency tuning sensitivity of $`df/d\theta f(d/l)`$ $`10^{14}\mathrm{Hz}/\mathrm{rad}`$. The angular measurement precision depends on the number of photons $`N`$ via $`\mathrm{\Delta }\theta \lambda /(dF\sqrt{N})`$, where $`F`$ is the cavity finesse, and $`\lambda `$ is the laser wavelength. This is a factor $`F`$ better than the optical lever readout for equal laser power. Cryogenic optical resonators have excellent long-term stability and have been proposed as secondary frequency standards. The measured frequency drifts range from $`1\mathrm{H}\mathrm{z}`$ over minutes to $`100\mathrm{Hz}`$ over days. For the measurement of the rotation angle, a stable reference frequency will be required. This can be easily implemented using a laser locked to a second cavity. The described angle readout has little sensitivity to lateral, tilt, and vertical motion, but couples to rotational noise. A possible solution would be to suspend the target as well as the optical cavity in order to suppress the common rotation mode. The cavity suspension should have a much longer natural period than the target suspension to avoid its excitation at the signal frequency. A different approach to suppress rotational noise would employ two identical torsion balances, but rotated by 180 degrees with respect to each other.

Additional background forces will arise from gas collisions, cosmic ray hits, radioactivity, solar neutrino and WIMP interactions resulting in a Brownian motion of the target. The equivalent acceleration is $`a(\overline{p}/m)\sqrt{n/\tau }`$, where $`\overline{p}`$ is the average momentum transfer, and $`n`$ is the collision or decay rate. The requirements on the residual gas pressures are very severe with only a few collisions with the target allowed per second. Cryopumping and the use of getters will be essential. The cosmic muon flux at sea level is about $`0.01\mathrm{cm}^2\mathrm{s}^1`$ with a mean energy of $`1\mathrm{GeV}`$. Thus, a target of area $`100\mathrm{cm}^2`$ will experience a collision rate $`𝒪(1)\mathrm{s}^1`$, causing an acceleration comparable to the signal in Equation 4. Cosmic rays can also lead to a net charge buildup on the torsion balance and spurious electrostatic forces. This background can be reduced by several orders of magnitude by going underground. Further disturbing forces may be caused by time-varying electric and magnetic background fields, which can be shielded with superconducters. Blackbody radiation and radiometric effects would be greatly reduced in a temperature controlled cryogenic environment.

Finally, there is a fundamental limit imposed by the uncertainty principle. The accuracy of a force measurement obtainable by continously recording the test mass position (standard quantum limit) is

$$a_{\mathrm{SQL}}=5\times 10^{24}\left(\frac{10\mathrm{k}\mathrm{g}}{m}\right)^{1/2}\left(\frac{1\mathrm{d}\mathrm{a}\mathrm{y}}{\tau _0}\right)^{1/2}\left(\frac{10^6\mathrm{s}}{\tau }\right)\frac{\mathrm{cm}}{\mathrm{s}^2}$$
(7)   
where $`m`$ is the test mass, $`\tau _0`$ is the oscillator period, and $`\tau `$ is the measurement time. In our proposed position readout, the disturbing back action force arises from spatial fluctuations in the photon flux passing through the central tuning plate. A laser beam pulse containing $`N`$ photons, and centered on the plate rotation axis, will produce a random change in the angular momentum of the plate of $`\mathrm{\Delta }LF\sqrt{N}\mathrm{}\omega d/c`$. Here, $`d`$ is the average transverse photon displacement with respect to the rotation axis, which is also of order the plate thickness. A slightly higher sensitivity can be obtained with a stroboscopic “quantum nondemolition” (QND) measurement , where the test mass position is recorded every half period. Here,

$$a_{\mathrm{QND}}=a_{\mathrm{SQL}}\sqrt{\omega _0\mathrm{\Delta }t}$$
(8)   
and $`\mathrm{\Delta }t`$ is the strobe duration.

There is hope that cosmic neutrinos will be detected in the laboratory early in the next century, especially if they have masses in the eV range and are of Dirac type. The most viable way seems to be a greatly improved torsion balance operated in a deep mine.

Naturally, the proposed torsion balance could also be applied to tests of the weak equivalence principle with the Sun as the source mass, and to searches for new short range $`(1\mathrm{m})`$ forces with a movable laboratory mass.

## ACKNOWLEDGMENTS

This research was performed by LLNL under the auspices of the U.S. Department of Energy under contract no. W-7405-ENG-48. I thank E. Adelberger, S. Baessler, I. Ferreras, A. Melissinos, R. Newman, P. Smith, and W. Stoeffl for useful discussions.