Physics:Neutron-antineutron oscillations

Short description: Hypothetical conversion between particles

Neutron-antineutron oscillations are a hypothetical process in which a neutron can spontaneously convert into an antineutron and vice versa. This process would violate baryon-number (B) conservation by two units ($\mathrm{\Delta }B=2$). No experimental evidence for neutron–antineutron oscillations has been found to date.

The possibility of baryon-number violation was identified by Andrei Sakharov in 1967 as one of the conditions required to explain the matter–antimatter asymmetry of the Universe. In 1970, Vadim Kuzmin proposed neutron–antineutron oscillations as a concrete realization of such violation. The conceptual foundation dates back to Ettore Majorana's 1937 hypothesis that neutrons could be identical to their own antiparticles. The theoretical framework was later developed by Sheldon Glashow in 1979, and formalized by Rabindra Mohapatra and Robert Marshak in 1980. These developments motivated experimental searches with both free neutrons and neutrons bound in nuclei.

Neutron–antineutron oscillations complement proton decay as a probe of baryon-number violation. While proton decay would violate baryon number by one unit, neutron–antineutron oscillations would violate it by two. The two processes explore very different mass scales (see Theoretical motivation). Baryon number conservation is predicted by the Standard Model of particle physics, although it is an accidental global symmetry, not associated with any fundamental gauge symmetry. Grand Unified Theories, which aim to unify the electromagnetic, weak and strong interactions into a single framework, naturally predict such violations.

Phenomenology

Neutron–antineutron ($n-\overline{n}$) oscillations can be described using a $2\times 2$ Hamiltonian governing the time evolution of a neutron–antineutron system:

(1)     $i\hslash \frac{\partial }{\partial t}\left(\begin{array}{c}n\\ \overline{n}\end{array}\right)=\left(\begin{array}{cc}{E}_n& \delta m\\ \delta m& E_{\overline{n}}\end{array}\right)\left(\begin{array}{c}n\\ \overline{n}\end{array}\right)$.

Here, $\delta m$ denotes the $n-\overline{n}$ mixing parameter, induced by the underlying baryon number violating processes, while $E_n$ and $E_{\overline{n}}$ are the energies of the neutron and antineutron, respectively, which may be affected by external fields.

From equation (1), the probability of observing an antineutron at time t, given an initial neutron state at t=0, is:

(2)     $P_{\overline{n}}(t)=\frac{4\delta m^2}{\mathrm{\Delta }{E}^2+4\delta m^2}{\sin }^2\left(\sqrt{\mathrm{\Delta }{E}^2+4\delta m^2}\right)t;$

where $\mathrm{\Delta }E=E_n-E_{\overline{n}}$. Neutron decay, which would introduce a factor $e^{-t/{\tau }_n}$ (with ${\tau }_n\simeq 880$ s), is neglected, as the relevant evolution times are much shorter than the neutron lifetime ${\tau }_n$.

Two limiting cases can be considered:

$n-\overline{n}$ oscillations of free neutrons

Even for free neutrons, the amplitude of the $P_{\overline{n}}(t)$ evolution is strongly suppressed. For oscillation times ${\tau }_{n\overline{n}}$ of order 10⁸ s, corresponding to the range probed experimentally, the mixing parameter $\delta m$ (with $\delta m=1/{\tau }_{n\overline{n}}$ in natural units) is of order 10^-29 MeV.

Neutron and antineutrons have opposite magnetic moments (${\mu }_n\simeq 9.7\cdot 10^{-27}$J/T), and in presence of the Earth's magnetic field ($B\simeq 50$ μT), the energy splitting is $\mathrm{\Delta }E=2{\mu }_{n}B\simeq 3\cdot 10^{-20}$MeV, leading to a suppression factor of order 10⁹.

A practical regime is the quasi-free limit , $|\mathrm{\Delta }E|t\ll 1$, in which the Taylor expansion of the oscillatory term in equation (2) compensates for the suppression arising from the denominator. In this limit, it results

(3)     $P_{\overline{n}}(t)\simeq [(\delta m)t{]}^2={\left(\frac{t}{{\tau }_{n\overline{n}}}\right)}^2$.

The quasi-free limit can be interpreted in terms of the energy–time uncertainty principle applied to the transition.

$n-\overline{n}$ oscillations in matter

Inside nuclei, neutron and antineutron potentials differ by approximately 100 MeV, which suppresses the oscillation probability by roughly 31 orders of magnitude. Despite this suppression, experimental searches remain feasible, necessitating a rigorous treatment of the underlying equation (2).

The total energies can be defined as $E_n=m_n+V_n$ and $E_{\overline{n}}=m_n+V_{\overline{n}}$. While the neutron nuclear potential $V_n$ is predominantly real, the antineutron potential $V_{\overline{n}}$ incorporates a substantial imaginary component to account for antineutron annihilation processes: $V_{\overline{n}}=V_{\overline{n}R}-i{V}_{\overline{nI}}$ where the real part satisfies $V_{\overline{n}R}\simeq V_n$ and the imaginary part scales as $V_{\overline{nI}}\sim O(100)$ MeV.

Diagonalization of the mass mixing matrix yields the following energy eigenvalues:

(4)     $E_{1,2}=1/2[E_n+E_{\overline{n}}\pm \sqrt{(E_n-E_{\overline{n}}{)}^2+4(\delta m{)}^2}]$.

A perturbative expansion of the E₁ state, which consists predominantly of the neutron component, leads to:

(5)     $E_1\simeq m_n+V_n-i\frac{(\delta m{)}^2{V}_{\overline{n}I}}{{(V_{nR}-V_{\overline{n}R})}^2+V_{\overline{n}I}^2}$

The imaginary part of E₁ describes matter instability via antineutron annihilation, characterized by the decay rate:

(6)     ${\mathrm{\Gamma }}_m=\frac{1}{T_{n\overline{n}}}=\frac{2(\delta m{)}^2{V}_{\overline{n}I}}{{(V_{nR}-V_{\overline{n}R})}^2+V_{\overline{n}I}^2}$

where the subscript m denotes "matter". Consequently, the characteristic lifetime scales as $T_{n\overline{n}}=1/{\mathrm{\Gamma }}_m\propto (\delta m{)}^{-2}$, that can be expressed as

(7)     $T_{n\overline{n}}=R{\tau }_{n\overline{n}}^2$.

Equation (7) relates the characteristic nuclear lifetime $T_{n\overline{n}}$, representing the time scale for an antineutron to be produced and subsequently annihilate within the nucleus, to the free neutron-antineutron oscillation time ${\tau }_{n\overline{n}}$. The conversion factor R, which has dimensions of s^-1, is a nucleus-dependent quantity that cannot be derived from first principles; instead, it must be evaluated using suitable nuclear models.

The overall theoretical uncertainty for these one-nucleon processes is approximately 10%–15%. This precision represents a significant reduction compared to the 50%–100% uncertainty margins characteristic of early calculations from the 1980s and 1990s. While these improvements cover single-nucleon mechanisms, an additional 15%–30% systematic uncertainty related to two-nucleon processes within the nuclear medium should be must be accounted for. The primary driver behind this reduction in systematic error is the integration of extensive and precise data from antiprotonic atoms, which became available after the earlier calculations were published.

Furthermore, several studies in the literature have raised questions regarding whether free neutron oscillations and oscillations in matter are mediated by the same underlying effective operators.

Experimental searches

$n-\overline{n}$ oscillations with free neutron beams

In an experiment with free neutron beams, the oscillation time ${\tau }_{n\overline{n}}$ can be expressed as:

(8)     ${\tau }_{n\overline{n}}=\sqrt{\frac{ITϵ}{\overline{N}}}t$;

where $\overline{N}$ is the number of detected antineutrons, $I$ the neutron intensity, T the running time, $ϵ$ the antineutron detection efficiency and t the neutron propagation time under quasi-free conditions. Based on Poisson statistics, an experiment detecting no antineutron events sets $\overline{N}=2.3$ at the 90% confidence level for limits on ${\tau }_{n\overline{n}}$. The most intense sources for cold neutrons are research nuclear reactors and spallation neutron sources (see also neutron sources).

In the early 1980s several experiments were proposed at neutron facilities such as Oak Ridge National Laboratory, the Omega West Reactor in Los Alamos, the Los Alamos Meson Physics Facility,, the Moscow Meson Factory, the Triga Mark II reactor at Pavia University and the nuclear reactor at the Institute Laue-Langevin (ILL) in Grenoble (see also the reviews in ). Only the last two experiments were ultimately carried out.

The ILL experiment reported the first experimental limit on ${\tau }_{n\overline{n}}$ with free neutrons in 1985, ${\tau }_{n\overline{n}}\ge 10^6$ s at 90% confidence level, while the NADIR collaboration in Pavia reported ${\tau }_{n\overline{n}}\ge 0.5\cdot 10^6$s in 1990.

Subsequently, groups from both collaborations joined to propose a new experiment at ILL ($N\overline{N}$) with a projected sensitivity of ${\tau }_{n\overline{n}}\ge 10^8$.

The experimental layout is detailed in the figure below.

File:Nnbar-ILL-layout.jpg

Layout of the $N\overline{N}$ neutron-antineutron experiment at the ILL reactor.

The experiment took data for a time $T=2.4\cdot 10^7$ s, less than the designed 1 year because of ILL reactor breakdown. With an antineutron detection efficiency $ϵ=0.52$ and no candidate events, it established a lower limit ${\tau }_{n\overline{n}}\ge 0.86\cdot 10^8$ s with a 90% confidence level.

$n-\overline{n}$ oscillations in matter

As noted in phenomenology, equation (7), experiments looking for neutron–antineutron oscillations in matter measure the intranuclear oscillation time $T_{n\overline{n}}$ and infer the free oscillation time ${\tau }_{n\overline{n}}$ by applying a nuclear factor R. Searches for $n-\overline{n}$ oscillations in matter are generally carried out by the same experiments used in proton decay searches (they are better described in proton decay). The difficulty of these searches arises from several factors: antineutron annihilation in nuclei leads to a wide variety of final states with different branching ratios, so there is no single distinctive experimental signature; in addition, the annihilation process typically yields four to five low-momentum pions that current detectors struggle to reconstruct efficiently; in addition, these pions may rescatter within the nucleus before detection, further degrading the signal.

Consequently, detection efficiencies for annihilation events are low, making it difficult to distinguish signal from atmospheric neutrino backgrounds. As a result, experimental lower limits on the $n-\overline{n}$ lifetime $T_{n\overline{n}}$ are substantially weaker than proton lifetime limits.

The experimental limits published so far are reported in the table below. In each case, $T_{n\overline{n}}$ is the one published by the experiment, whereas ${\tau }_{n\overline{n}}$ is evaluated by applying the most recent computation for the nuclear factor R (and can differ by factors 2-3 from the original published values). The earliest results were published in 1983 by the water Cherenkov Homestake experiment, and by the tracking calorimeter Nusex, installed in the Mont Blanc Tunnel, Italy.

The most stringent limit has been published by the Super-Kamiokande experiment in 2021. The experiment analyzed an exposure of 0.37 Mton-years, corresponding to approximately 16.5 years of data taking with a fiducial volume of 22.5 kton. The total signal efficiency was 4.1% (with a 33% systematic uncertainty), and the expected background was 9.3 events over the full data set (with a 28% systematic uncertainty). The experiment observed 11 candidate events establishing a limit $T_{n\overline{n}}\ge 3.6\cdot 10^{32}$ yr at 90% confidence level, corresponding to ${\tau }_{n\overline{n}}\ge 4.7\cdot 10^8$ s.

Results of $n-\overline{n}$ oscillation searches from bound neutrons.

Year	Nucleus	Experiment	$T_{n\bar{n}}$ (1032 yr)	R (1023/s)	${\tau }_{n\bar{n}}$ (108 s)
1983	16O	Homestake	0.014	0.52	0.07
1983	56Fe	Nusex	0.6	1.4	1.0
1984	16O	0.24	0.52	1.2
1986	16O	0.4	0.52	1.6
1986	56Fe	0.3	1.4	0.5
1990	56Fe	0.6	1.4	1.2
2002	56Fe	0.7	1.4	1.3
2017	2H	0.1	0.25	1.4
2021	16O	Super-K	3.6	0.52	4.7

Future initiatives

Proposals for experiments with free-neutron beams were published after the conclusion of the $N\overline{N}$ experiment at facilities as the HFIR nuclear reactor at Oak Ridge National Laboratory, the WWR-M nuclear reactor at Saint Petersburg, or the polarised cold neutron beam at the Institute Laue-Langevin at Grenoble. Unfortunately, none of them has been realized.

The HIBEAM/NNBAR collaboration has proposed a two-stage experiment at the European Spallation Source (ESS). The first stage, HIBEAM, is designed as a pilot program during the early commissioning phase of the ESS and will search for $n-\overline{n}$ transitions without using the facility's full planned beam power. Its sensitivity is not expected to surpass the current limit on ${\tau }_{n\overline{n}}$. The second stage, NNBAR, would employ the full ESS beam power together with high-reflectivity supermirror reflectors, either ellipsoidal or differential, that could collect a larger fraction of the neutron flux and focus it onto the target. Although the final configuration remains under development, the projected sensitivity is ${\tau }_{n\overline{n}}\gtrsim 2.6\cdot 10^9$s.

Experiments searching for $n-\overline{n}$ oscillations in matter have largely reached the limits of their potential for significant improvement, as low efficiencies and substantial background subtraction limit further sensitivity gains. The Deep Underground Neutrino Experiment (DUNE) experiment in the United States scheduled to begin data collection in 2031, utilizes liquid argon time projection chamber (LArTPC) technology. This approach is expected to improve the efficiency and purity of the collected sample due to its advanced tracking capabilities. DUNE's expected sensitivity is estimated at ${\tau }_{n\overline{n}}\ge 5.5\cdot 10^8$s after ten years of data taking in its full configuration. The successor of Super-Kamiokande, Hyper-Kamiokande, has not yet released a prediction for its sensitivity on ${\tau }_{n\overline{n}}$; however, a baseline extrapolation based on ten years of exposure yields an estimated sensitivity of ${\tau }_{n\overline{n}}\gtrsim 10^9$s.

Theoretical motivation

File:Nnosc-so10.png

Feynman diagram for neutron-antineutron oscillation in SO(10)

File:Nnosc-susy.png

Feynman diagram for neutron-antineutron oscillation in supersymmetry

At the quark level, the $n-\overline{n}$ transition converts three quarks into three antiquarks (udd → u^cd^cd^c ). This process violates baryon number conservation by 2 units ($\mathrm{\Delta }B=2$) while conserving lepton number ($\mathrm{\Delta }L=0$). It requires six-quark operators; the corresponding amplitude has mass dimension 9 and scales as ${\lambda }_{B-L}^{-5}$, where ${\lambda }_{B-L}$denotes the energy scale of (B−L) violation. The diquark scalars required to mediate this process are not present in the Standard Model but arise naturally in some grand unified theories (GUT). In contrast to proton decay, GUTs do not provide robust predictions for the neutron–antineutron oscillation time ${\tau }_{n\overline{n}}$ (see,e.g., for a review).

The SU(5) group, introduced in 1974 by Georgi and Glashow, cannot accommodate $\mathrm{\Delta }B=2$ processes. In the minimal SU(5) model, the difference between baryon number (B) and lepton number (L),denoted as B−L, remains an exact global symmetry. Extensions of the model involving higher-dimensional Higgs multiplets can break this symmetry, allowing neutrinos to acquire Majorana masses while simultaneously providing the operators needed to mediate neutron oscillations.

SO(10) GUTs are a natural framework for $n-\overline{n}$ oscillations because they allow B−L to be a gauged symmetry. Spontaneous breaking of this symmetry by two units (Δ(B−L)=2) creates a deep theoretical link between Majorana neutrino masses (via the seesaw mechanism) and $n-\overline{n}$ transitions. While standard GUT scales are near 10¹⁶ GeV, a restricted class of SO(10) models can support intermediate scales ( ${\lambda }_{B-L}\sim 10^2-10^3$ TeV) where oscillations become experimentally observable. In SO(10) $n-\overline{n}$ oscillations can be mediated by color-sextet scalar diquark fields, as illustrated in figure. Specifically, the post-sphaleron baryogenesis scenario within these models predicts an upper limit for the oscillation time of 5×10¹⁰ seconds.

Supersymmetry (SUSY) significantly modifies the operators mediating $n-\overline{n}$ oscillations by introducing superpartners like squarks and gluinos. These additional states can alleviate the strong suppression present in the Standard Model by allowing lower-dimensional operators, for example, of dimension 4 or 5. As a consequence, the oscillation time scales as ${\tau }_{n\overline{n}}\propto {\lambda }_{B-L}^2{\lambda }_{SM}^3$ rather than ${\lambda }_{B-L}^5$ (where ${\lambda }_{SM}$ denotes the Standard Model energy scale). This behavior can lead to detectable oscillation times (∼10¹⁰ s) even at very high scales of ${\lambda }_{B-L}\sim 10^8-10^{11}$ GeV. A representative Feynman diagram for such a transition is shown in Figure.

Models that propagate Standard Model fields into extra dimensions predict ${\tau }_{n\overline{n}}$ of the order of 10⁹ s.

Further reading

• Mohapatra, R N (2009-09-16). "Neutron–anti-neutron oscillation: theory and phenomenology". Journal of Physics G: Nuclear and Particle Physics 36 (10). doi:10.1088/0954-3899/36/10/104006. ISSN 0954-3899. Bibcode: 2009JPhG...36j4006M. 
• Phillips, D. G. et al. (2016-02-11). "Neutron-antineutron oscillations: Theoretical status and experimental prospects". Physics Reports 612: 1–45. doi:10.1016/j.physrep.2015.11.001. ISSN 0370-1573. Bibcode: 2016PhR...612....1P. https://www.sciencedirect.com/science/article/pii/S0370157315004457. 
• Dev, P S B et al. (2024-01-23). "Searches for baryon number violation in neutrino experiments: a white paper". Journal of Physics G: Nuclear and Particle Physics 51 (3): 033001. doi:10.1088/1361-6471/ad1658. ISSN 0954-3899. Bibcode: 2024JPhG...51c3001D. 
• Broussard, Leah J. (2025). "Baryon number violation: from nuclear matrix elements to BSM physics". Journal of Physics G: Nuclear and Particle Physics 52 (8): 083001. doi:10.1088/1361-6471/adf081. Bibcode: 2025JPhG...52h3001B.