A New Limit on Time-Reversal Violation in Beta Decay

Abstract

We report the results of an improved determination of the triple correlation $D𝐏\cdot ({𝐩}_{𝐞}\times {𝐩}_{𝐯})$ that can be used to limit possible time-reversal invariance in the beta decay of polarized neutrons and constrain extensions to the Standard Model. Our result is $D=(-0.96\pm 1.89(stat)\pm 1.01(sys))\times {10}^{-4}$. The corresponding phase between $g_A$ and $g_V$ is ${\varphi }_{AV}={180.013}^{\circ }\pm {0.028}^{\circ }$ (68% confidence level). This result represents the most sensitive measurement of $D$ in nuclear beta decay.

The existence of charge-parity (CP) symmetry violation in nature is particularly important in that it is necessary to explain the preponderance of matter over antimatter in the universe (1). Thus far, CP violation has been observed only in the K and B meson systems (2; 3; 4) and can be entirely accounted for by a phase in the Cabbibo-Kobayashi-Maskawa matrix in the electroweak Lagrangian. This phase is insufficient to account for the known baryon asymmetry in the context of Big Bang cosmology (5), so there is good reason to search for CP violation in other systems. As CP and time-reversal (T) violation can be related to each other through the CPT theorem, experimental limits on electric dipole moments and T-odd observables in nuclear beta decay place strict constraints on some, but not all, possible sources of new CP violation.

Figure 1: A schematic of the emiT detector illustrating the alternating electron and proton detector segments. The darker shaded proton detectors indicate the the paired-ring at $z=\pm 10$ cm. The cross section view illustrates, in an exaggerated manner, the effect of the magnetic field on the particle trajectories and average opening angle. A P${}_2$E${}_3$ coincidence event is shown.

The decay probability distribution for neutron beta decay, $dW$, can be written in terms of the beam polarization $𝐏$ and the momenta (energies) of the electron ${𝐩}_e$ ($E_e$) and antineutrino ${𝐩}_{\nu }$ ($E_{\nu }$) as (6)

	$dW\propto 1+a{\displaystyle \frac{{𝐩}_e\cdot {𝐩}_{\nu }}{E_e{E}_{\nu }}}+b{\displaystyle \frac{m_e}{E_e}}+$						(1)
			$\mathrm{ }𝐏\cdot \left(A{\displaystyle \frac{{𝐩}_e}{E_e}}+B{\displaystyle \frac{{𝐩}_{\nu }}{E_{\nu }}}+D{\displaystyle \frac{{𝐩}_e\times {𝐩}_{\nu }}{E_e{E}_{\nu }}}\right).$

A contribution of the parity-even triple correlation $D𝐏\cdot ({𝐩}_e\times {𝐩}_{\nu })$ above the level of calculable final-state interactions (FSI) implies T-violation. The PDG average of recent measurements is $D=(-4\pm 6)\times {10}^{-4}$ (7; 8; 9), while the FSI for the neutron are $\sim {10}^{-5}$ (10; 11). Complementary limits can be set on other T-violating correlations, and recently a limit on $R$ has been published (12). Various theoretical models that extend the SM, such as left-right symmetric theories, leptoquarks, and certain exotic fermions could give rise to observable effects that are as large as the present experimental limits (13). Calculations performed within the Minimal Supersymmetric Model, however, predict $D\lesssim {10}^{-7}$ (14).

In the neutron rest frame, the triple correlation can be expressed as $D𝐏\cdot ({𝐩}_p\times {𝐩}_e)$, where ${𝐩}_p$ is the proton momentum. Thus one can extract $D$ from the spin dependence of proton-electron coincidences in the decay of cold polarized neutrons. Our measurement was carried out at the National Institute of Standards and Technology Center for Neutron Research (15). The detector, shown schematically in Fig. 1, consisted of an octagonal array of four electron-detection planes and four proton-detection planes concentric with a longitudinally polarized beam. The beam, with a neutron capture fluence rate at the detector of $1.7\times {10}^8$ cm${}^{-2}$ s${}^{-1}$, was defined using a series of ${}^6$LiF apertures and polarized to $>91\%$ (95% C.L.) by a double-sided bender-type supermirror (15). A 560 $\mu$T guide field maintained the polarization direction throughout the fiducial volume and a current-sheet spin-flipper was used to reverse the neutron spin direction every 10 s. The symmetric octagonal geometry was chosen to maximize sensitivity to $D$ while approximately canceling systematic effects stemming from detector efficiency variations or coupling to the spin correlations $A$ and $B$ (16; 8). Each of the four proton segments consisted of a $2\times 8$ array of silicon surface-barrier detectors (SBDs) with an active layer 300 mm${}^2$$\times$ 300 $\mu$m. Each SBD was contained within an acceleration and focusing cell consisting of a 94% transmitting grounded wire-mesh box through which the recoil protons entered. Each SBD, situated within a field-shaping cylindrical tube, was held at a fixed voltage in the range $-25$ kV to $-32$ kV. The sensitive regions of the beta detectors were plastic scintillator measuring 50 cm by 8.4 cm by 0.64 cm thick, with photomultiplier tube (PMT) readout at both ends. This thickness is sufficient to stop electrons at the decay endpoint energy of 782 keV. The proton and beta detectors were periodically calibrated in situ with gamma and beta sources respectively. Details of the apparatus are presented elsewhere (17; 15; 8).

Data were acquired in a series of runs from October 2002 through December 2003. Typical count rates were 3 s${}^{-1}$ and 100 s${}^{-1}$ for single proton and beta detectors, respectively, while the coincidence rate for the entire array was typically 25 s${}^{-1}$. Of the raw events, 12% were eliminated by filtering on various operational parameters (e.g. coil currents) and by requiring equal counting time in each spin-flip state. A beta-energy software threshold of 90 keV eliminated detection efficiency drifts due to changes in PMT gain coupled with the hardware threshold. This was the largest single cut, eliminating 14% of the raw events. A requirement that a single beta be detected in coincidence with each proton eliminated 7% of events. All cuts were varied to test for systematic effects.

Figure 2: Intensity log plot of SBD-scintillator coincidence data showing proton energy vs delay time. Events near $\mathrm{\Delta }t=0$ are prompt coincidences due primarily to beam-related backgrounds.

The remaining coincidence events were divided into two timing windows: a preprompt window from -12.3 $\mu$s to -0.75 $\mu$s that was used to determine the background from random coincidences, and the decay window from -0.5 $\mu$s to 6.0 $\mu$s as shown in Fig. 2. The recoil proton has an endpoint of 750 eV. On average it is delayed by $\sim$ 0.5 $\mu$s. The average signal-to-background ratio was $\sim$ 30/1. The energy-loss spectrum produced by minimum ionizing particles in 300 $\mu$m of silicon is peaked at approximately 100 keV and, being well separated from the proton energy spectrum, yielded an estimated contamination below 0.1%. The final data set consisted of approximately 300 million accepted coincidence events.

A detailed Monte Carlo simulation was used to estimate a number of systematic effects. The program Penelope (18), which has been tested against data in a variety of circumstances of relevance to neutron decay (19), was embedded within a custom tracking code. All surfaces visible to decay particles were included. The Monte Carlo was based on the measured beam distribution upstream and downstream of the fiducial volume (15) and incorporated the magnetic field and electron energy threshold. A separate Monte Carlo based on the package SIMION (20), incorporating the detailed geometry of the proton cells, was used to model the proton detection response function.

Achieving the desired sensitivity to $D$ in the presence of the much larger spin-asymmetries due to $A$ and $B$ depends critically on the measurement symmetry. To the extent that this symmetry is broken, corrections must be applied to the measured result. These corrections are listed in Table 1 and are discussed below. To extract $D$, coincident events were first combined into approximately efficiency-independent asymmetries

$$w^{p_i{e}_j}=\frac{N_{+}^{p_i{e}_j}-N_{-}^{p_i{e}_j}}{N_{+}^{p_i{e}_j}+N_{-}^{p_i{e}_j}},$$

(2)

where $N_{+}^{p_i{e}_j}$ is the integrated number of coincident events in proton detector $i=\mathrm{1...64}$, beta detector $j=\mathrm{1...4}$, with neutron spin + ($-$) aligned (anti-aligned) with the guide field. For uniform polarization, $𝐏$, the asymmetries, $w^{p_i{e}_j}$, can be written in terms of decay correlations as

$$w^{p_i{e}_j}\approx 𝐏\cdot (A{\widetilde{𝐊}}_A^{p_i{e}_j}+B{\widetilde{𝐊}}_B^{p_i{e}_j}+D{\widetilde{𝐊}}_D^{p_i{e}_j}),$$

(3)

where the $𝐊$’s are obtained from Eqn. 1 by integrating the normalized kinematic terms over the phase space of the decay, the neutron beam volume, and the acceptance of the indicated detectors (8). ${\widetilde{𝐊}}_A\propto ⟨{𝐩}_e/E_e⟩$ and ${\widetilde{𝐊}}_B\propto ⟨{𝐩}_{\nu }/E_{\nu }⟩$ are primarily transverse to the detector axis but have roughly equal longitudinal components for coincidence events involving the two beta detectors opposite from the indicated proton detector (E${}_2$ and E${}_3$ for P${}_2$ as shown in Fig. 1). The ${\widetilde{𝐊}}_D$’s, however, are primarily along the detector axis and are opposite in sign for the two beta detectors. Thus for each proton detector we can choose an appropriate combination of detector pairs that is sensitive to the $D$-correlation but that largely cancels the parity-violating $A$ and $B$ correlations. One such combination is

$$v^{p_i}=\frac{1}{2}(w^{p_i{e}_R}-w^{p_i{e}_L}),$$

(4)

where $e_R$ and $e_L$ label the electron-detector at approximately 135${}^{\circ }$ giving a positive and negative cross-product ${𝐩}_p\times {𝐩}_e$ respectively (P${}_2$E${}_3$ vs P${}_2$E${}_2$ as shown in Fig. 1). Proton cells at the detector ends accept decays with larger longitudinal components of ${\widetilde{𝐊}}_A$ and are more sensitive to a range of effects that break the detector symmetry. We therefore define $\overline{v}$ as the average of the values of $v$ from the sixteen proton-cells at the same $|z|$, i.e. $z=\pm$2, $\pm$6, $\pm$10, and $\pm$14 cm. Each set of detectors corresponds to a paired-ring with the same symmetry as the full detector, e.g. the shaded detectors in Fig. 1. We then define

$$\tilde{D}=\frac{\overline{v}}{P{\overline{K}}_D},$$

(5)

where ${\overline{K}}_D=0.378$ is the average of $\widehat{z}\cdot ({\widetilde{𝐊}}_D^{p_i{e}_R}-{\widetilde{𝐊}}_D^{p_i{e}_L})$ determined by Monte Carlo. The experiment provides four independent measurements corresponding to each of the four paired-rings. The systematic corrections to $\tilde{D}$ presented in Table I yield our final value for $D$. Eqn. 5 is based on the following: 1) accurate background corrections, 2) uniform proton and electron detection efficiencies, 3) cylindrical symmetry of the neutron beam and polarization, and 4) accurate determination of ${\overline{K}}_D$, $P$, and spin state.

Table 1: Systematic corrections and combined standard uncertainties (68% confidence level). Values should be multiplied by ${10}^{-4}$.

Source	Correction	Uncertainty
Background asymmetry	$0$ ${}^{\mathrm{a}}$	$0.30$
Background subtraction	$0.03$	$0.003$
Electron backscattering	$0.11$	$0.03$
Proton backscattering	$0$ ${}^{\mathrm{a}}$	$0.03$
Beta threshold	$0.04$	$0.10$
Proton threshold	$-0.29$	$0.41$
Beam expansion, magnetic field	$-1.50$	$0.40$
Polarization non-uniformity	$0$ ${}^{\mathrm{a}}$	$0.10$
ATP - misalignment	$-0.07$	$0.72$
ATP - Twist	$0$ ${}^{\mathrm{a}}$	$0.24$
Spin-correlated flux ${}^{\mathrm{b}}$	$0$ ${}^{\mathrm{a}}$	$3\times {10}^{-6}$
Spin-correlated pol.	$0$ ${}^{\mathrm{a}}$	$5\times {10}^{-4}$
Polarization ${}^{\mathrm{c}}$		$0.04$ ${}^{\mathrm{d}}$
${\overline{K}}_D$ ${}^{\mathrm{c}}$		$0.05$
Total systematic corrections	$-1.68$	$1.01$

${}^{\mathrm{a}}$ Zero indicates no correction applied. ${}^{\mathrm{b}}$ Includes spin-flip time, cycle asymmetry, and flux variation. ${}^{\mathrm{c}}$ Included in the definition of $\tilde{D}$. ${}^{\mathrm{d}}$ Assumed polarization uncertainty of 0.05.

Backgrounds not properly accounted for contribute two systematic errors: 1) multiplicative errors due to dilution of the asymmetries, and 2) spin-dependent backgrounds that can lead to a false $D$. Errors in background subtraction, as well as possible spin-dependent asymmetries in this background, have a small effect. The multiplicative correction to the value of $w^{p_i{e}_j}$ due to backscattered electrons was determined using the Monte Carlo. The uncertainty given in Table 1 reflects the 20% uncertainty assigned to the backscattering fractions due to limitations of the detector and beam model and due to limited knowledge of backscattering at energies below a few hundred keV. Proton backscattering, though observable, produces a negligible effect on $\tilde{D}$.

In principle, the values of $w^{p_i{e}_j}$ are independent of the absolute efficiencies of the proton and electron detectors; however, they do depend on any energy dependence of the efficiencies through the factors $⟨{𝐩}_e/E_e⟩$. Spatial variation of the efficiencies breaks the symmetry assumed in combining proton-cell data into paired-rings. Beta energy thresholds were observed to vary less than 20 keV across the detector, implying the almost negligible correction given in Table 1. Proton detector efficiency variations, however, were more significant. Lower energy proton thresholds varied across the detector and over the course of the experiment. These thresholds combined with the spin-dependence of the accelerated proton energy spectra can result in significant deviations in the values of $w^{p_i{e}_j}$, though the effect on the the value of $v^{p_i}$ is largely mitigated because the low-energy portion of the proton energy spectrum is roughly the same for the $e_R$ and $e_L$ coincidence pairs. To estimate the proton-threshold-nonuniformity effect on $\tilde{D}$, spin-dependent proton energy spectra were generated by Monte Carlo for all proton-detector-electron-detector pairings and convoluted with model detector response functions based on fits to the average proton-SBD spectra. The average fit parameters were varied over a range characteristic of the observed variations during the run. Representative thresholds were then applied to determine the effect on $\tilde{D}$. An alternative and consistent estimate was derived by correcting the values of $w^{p_i{e}_j}$ on a day-by-day basis using the spectrum centroid shift and empirically determined functional form of the spectrum at the threshold.

Beam expansion from a radius of 2.5 cm to 2.75 cm combined with the magnetic field breaks the symmetry of the detector because the average proton-electron opening angle for each proton-electron detector pair is modified. Monte Carlo calculations using measured upstream and downstream density profile maps were used to calculate the correction given in Table 1. Possible inaccuracies in the determination of the beam density were estimated and their implications explored with the Monte Carlo code. The effect on the value of $v^{p_i}$ as a function of position is illustrated in Fig. 3.

Figure 3: Solid (open) squares show the values of $v$ averaged over the four planes for proton cells on the even (odd) side of the proton detection plane (the side with P${}_2$ vs P${}_1$ in Fig. 1). Monte Carlo results are indicated by lines. The broken symmetry, due to the combination of the magnetic field and beam expansion, is evident in the shift in the crossing point from the detector center.

For a symmetric beam, contributions to $v^{p_i}$ due to transverse polarization cancel for opposing proton planes; however small azimuthal beam asymmetries can affect this cancellation. This asymmetric-beam transverse-polarization (ATP) effect is proportional to both $\sin {\theta }_P$, the angle of the average neutron-spin orientation with respect to the detector axis, and $\sin ({\varphi }_P-{\varphi }^{p_i})$, where ${\varphi }^{p_i}$ is the effective azimuthal position of the proton cell, and ${\varphi }_P$ is the azimuthal direction of the neutron polarization. To study this effect transverse-polarization calibration runs with ${\theta }_P={90}^{\circ }$ and several values of ${\varphi }_P$ were taken over the course of the experiment. In these runs the ATP effect was amplified by $\approx 200$. The values of $\sin {\theta }_P$ and ${\varphi }_P$ for the experiment were determined using the calibration runs and Monte Carlo corrections for the beam density variations. To estimate the effect, the extreme value of sin${\theta }_P=12.8\times {10}^{-3}$ and the range of were used. The uncertainty is due to uncertainties in the angles ${\theta }_P$ and ${\varphi }_P$. The effect of nonuniform beam polarization is also given in Table 1. Time-dependent variations in flux, polarization, and the spin-flipper, as well as the uncertainty in the instrumental constant ${\overline{K}}_D$, can be shown to produce asymmetries proportional to $\tilde{D}$. These effects are listed in Table 1.

Correlations of $\tilde{D}$ with a variety of experimental parameters were studied by varying the cuts and by breaking the data up into subsets taken under different conditions of proton acceleration voltage and number of live SBDs as shown in Fig. 4. A linear correlation of $\tilde{D}$ with high-voltage, revealed by the cuts study, yields ${\chi }^2=5.6$ with 11 DOF compared to 10.4 for 12 DOF for no correlation. The acceleration-voltage dependence of the focusing properties was extensively studied by Monte Carlo with no expected effect, and we intrepret the 2.1 sigma slope as a statistical fluctuation.

To improve the symmetry of the the detector, we combine counts for the entire run to determine the values of $v^{p_i}$ and to extract $\tilde{D}$ for each paired-ring. A blind analysis was performed by adding a constant hidden factor to Eqn. 2. The blind was only removed once all analyses of systematic effects were complete and had been combined into a final result. The weighted average of $\tilde{D}$ for the four paired-rings is $0.72\pm 1.89$ with ${\chi }^2=0.73$ for 3 DOF.

Figure 4: Results for $\tilde{D}$ by run subset. Uncertainties are statistical. The fully functioning paired-rings used for each subset are indicated: 1-4 indicates all paired-rings were used.

The result including all corrections from Table 1 is

$$D=(-0.96\pm 1.89(stat)\pm 1.01(sys))\times {10}^{-4}.$$

Our result represents the most sensitive measurement of the $D$ coefficient in nuclear beta decay. Assuming purely vector and axial-vector currents, ${\varphi }_{AV}={180.013}^{\circ }\pm {0.028}^{\circ }$ which is the best direct determination of a possible CP-violating phase between the axial and vector currents in nuclear beta decay. Previously the most sensitive measurement was in ${}^{19}$Ne, with $D=(1\pm 6)\times {10}^{-4}$ (21).

The authors acknowledge the support of the National Institute of Standards and Technology, U.S. Department of Commerce, in providing the neutron facilities used in this work. This research was made possible in part by grants from the U.S. Department of Energy Office of Nuclear Physics (DE-FG02-97ER41020, DE-AC02-05CH11231, and DE-FG02-97ER41041) and the National Science Foundation (PHY-0555432, PHY-0855694, PHY-0555474, and PHY-0855310).