⚛️ lattice

Calculation of neutron electric dipole moment from Lattice QCD

This paper presents a novel lattice QCD calculation of the neutron electric dipole moment induced by the QCD theta term using the Feynman-Hellmann theorem and background electric fields, yielding a result that confirms the strong-CP problem by constraining the theta angle to $|\bar{\theta}|\lesssim 10^{-11}$ .

Original authors: Thomas Blum, Fangcheng He, Taku Izubuchi, Luchang Jin, Hiroshi Ohki, Sergey Syritsyn

Published 2026-08-25

📖 7 min read🧠 Deep dive

CC0 1.0

Original authors: Thomas Blum, Fangcheng He, Taku Izubuchi, Luchang Jin, Hiroshi Ohki, Sergey Syritsyn

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The universe we inhabit is built on a delicate balance of forces, but one of the most profound mysteries in physics concerns a subtle imbalance between matter and antimatter. In the early moments after the Big Bang, the universe should have produced equal amounts of both, which would have annihilated each other instantly, leaving nothing but light. Yet, we exist, made entirely of matter. To explain this, physicists know that the laws of nature must have treated matter and antimatter slightly differently, a phenomenon known as symmetry violation. While scientists have observed this difference in the weak nuclear force, the amount they have found is far too small to account for the vast abundance of matter we see today. This leaves a gap in our understanding, suggesting that another, hidden source of this difference must exist within the strong nuclear force, the glue that holds atomic nuclei together.

The search for this missing piece focuses on a property called the electric dipole moment of the neutron. A neutron is a neutral particle found in the heart of every atom, but if it possesses an electric dipole moment, it means its internal charge distribution is slightly shifted, creating a tiny separation between positive and negative charge. Such a shift would be a direct signature of the symmetry violation needed to explain our existence. However, experiments have never found this shift; they have only set increasingly strict limits on how small it could possibly be. The fact that it remains undetected, despite the theoretical need for it, points to a deep puzzle in our understanding of the strong force, often called the strong-CP problem.

A team of researchers has now taken a significant step toward solving this puzzle by calculating how the neutron should behave under the influence of this theoretical symmetry violation, using the most powerful tool available to study the subatomic world: a supercomputer simulation of the strong force itself. By recreating the conditions of the early universe in a digital grid, the scientists were able to measure the neutron's response to a specific type of symmetry breaking that had previously been too difficult to calculate with precision. Their work provides a clear, reliable prediction for the size of the neutron's electric dipole moment, confirming that the theoretical problem is real and that the universe's current state requires a very specific, tiny adjustment to the laws of physics to exist.

To understand what the researchers did, one must first visualize the strong force not as a smooth field, but as a turbulent, churning sea of interactions. In the standard model of particle physics, the strong force is described by a theory called quantum chromodynamics, which governs how quarks and gluons interact. At the low energies found inside a neutron, these interactions are so intense that traditional mathematical tools fail; the equations become too complex to solve with pen and paper. Instead, physicists use lattice quantum chromodynamics, a method that breaks space and time into a grid of tiny points, allowing a computer to calculate the behavior of particles by summing up the interactions at each point. This approach is the only way to study the strong force from first principles without relying on approximations that might miss crucial details.

The specific challenge the team tackled was how to measure the neutron's electric dipole moment in this digital environment. In previous attempts, researchers tried to calculate this value by looking at the total "topological charge" of the entire simulated universe. Topological charge is a measure of how the gluon fields twist and turn within the vacuum. However, calculating this total charge across the entire grid introduced massive amounts of random noise, making it nearly impossible to distinguish the tiny signal of the dipole moment from the background static. It was like trying to hear a whisper in a hurricane; the signal was there, but it was drowned out.

The new study introduced a clever workaround by changing how the neutron was observed. Instead of looking at the entire universe at once, the researchers placed the neutron in a uniform electric field, much like placing a compass needle in a magnetic field. This field does not change the fundamental laws of physics but forces the neutron to align in a specific way, mixing its internal states. In this mixed state, the neutron becomes sensitive to the local, or single-point, topological charge rather than the total charge of the whole system. By focusing on this local density, the team effectively turned down the volume on the noise, allowing the signal to emerge clearly. They used a mathematical technique called the Feynman-Hellmann theorem, which relates the energy shift of a particle in a field to the force acting on it, to extract the dipole moment from these local measurements.

The team ran these simulations on three different sets of digital grids, each with a different size and level of detail, corresponding to different masses for the particles involved. They used a specific type of particle interaction model known as domain wall fermions, which is known for preserving certain symmetries that other models struggle with. By carefully analyzing the data from these simulations, they were able to isolate the contribution of the symmetry-breaking term to the neutron's electric dipole moment. The result was a statistically significant signal, a clear measurement that had been elusive in previous attempts. They found that the neutron's electric dipole moment is proportional to a specific parameter in the theory, with a value of approximately minus 0.0050 times that parameter, measured in units of electron times femtometer.

This finding is crucial because it allows physicists to connect the theoretical prediction directly to experimental reality. The current best experimental limit on the neutron's electric dipole moment is incredibly small, far smaller than what the strong force alone would produce if the symmetry-breaking parameter were large. By comparing their calculated value with the experimental limit, the researchers determined that the symmetry-breaking parameter must be smaller than one part in one hundred billion. This tiny, unnatural number is the essence of the strong-CP problem: the laws of physics seem to conspire to make this parameter vanishingly small, yet there is no obvious reason why it should be so. The study confirms that the problem is not a flaw in the theory but a genuine feature of our universe that requires an explanation.

The reliability of this result stems from the method's ability to control errors that have plagued previous studies. The researchers tested their approach by varying the types of particles they used to represent the neutron and by checking the results across different simulation parameters. They found that as long as they properly accounted for the excited states of the neutron—temporary, higher-energy configurations that can contaminate the measurement—their results remained consistent. This robustness gives them confidence that the signal they found is real and not an artifact of the simulation. Unlike earlier methods that required extrapolating data to a limit that was hard to reach, this approach measures the quantity directly in a stable state, reducing the uncertainty associated with the calculation.

The implications of this work extend beyond just the neutron. The method developed by the team is versatile and can be applied to other sources of symmetry violation that might exist in nature, such as interactions involving multiple gluons or four-quark operators. These are other potential candidates for the physics that lies beyond our current understanding, and having a reliable way to calculate their effects is a major step forward. The ability to simulate these interactions with high precision means that future experiments can be guided by more accurate theoretical predictions, narrowing the search for new physics.

Ultimately, this study does not solve the mystery of why the universe is made of matter, but it provides a precise map of the terrain where the answer lies. It confirms that the strong-CP problem is a real, quantifiable constraint on our theories, and it offers a new, reliable tool for exploring the hidden symmetries of the subatomic world. The researchers have shown that by looking at the local details of the quantum vacuum rather than the global picture, we can hear the whisper of the neutron's dipole moment clearly. This clarity brings us closer to understanding the fundamental rules that allowed our existence, even as it deepens the mystery of why those rules are so finely tuned.

Technical Summary: Calculation of Neutron Electric Dipole Moment from Lattice QCD

Problem Statement
The neutron electric dipole moment (nEDM) is a critical observable for probing Charge–Parity (CP) violation, which is a necessary condition for the Baryonic Asymmetry of the Universe. While the Standard Model's weak interactions provide insufficient CP violation to explain observed baryon asymmetry, the QCD $\Theta$ -term offers a potential source. However, determining the contribution of strong interactions to the nEDM is theoretically challenging due to the non-perturbative nature of low-energy QCD.

Previous Lattice QCD calculations of the nEDM induced by the $\Theta$ -term have faced significant hurdles. Traditional methods extract the nEDM from the forward limit ( $Q^2 \to 0$ ) of the electric dipole form factor $F_3(Q^2)$ , computed via three-point correlation functions involving CP-odd operators. These approaches suffer from severe statistical noise due to the fluctuation of the global topological charge operator. Consequently, obtaining a statistically significant signal at the physical point has remained elusive, and even at heavier quark masses, results often require truncated-volume summations that introduce systematic uncertainties. Additionally, the extrapolation of $F_3(Q^2)$ to $Q^2=0$ introduces further systematic errors.

Methodology
This work presents a novel determination of the nEDM using the background field method combined with the Feynman-Hellmann theorem. Instead of computing global topological charge correlations, the authors calculate the nEDM from the matrix element of the local (single time-slice) topological charge density between nucleon ground states that are spatially polarized by a uniform background electric field.

Key methodological features include:

Background Field Implementation: A uniform electric field $\vec{E}$ is imposed on the lattice by modifying SU(3) gauge links with U(1) phases, satisfying 't Hooft quantization conditions. The field is treated as an infinitesimal perturbation.
Feynman-Hellmann Relation: The nEDM is related to the derivative of the nucleon energy with respect to the physical theta angle ( $\bar{\theta}$ ). This derivative is expressed as the matrix element of the local topological charge operator $q_{top}$ between the left and right eigenstates of the nucleon ground state in the presence of the electric field.
Variational Analysis (GEVP): The electric field breaks parity, causing the physical nucleon ground state to become a mixture of positive- and negative-parity states. To isolate this ground state and suppress excited state contamination, the authors employ a Generalized Eigenvalue Problem (GEVP) analysis. They construct a non-Hermitian correlation matrix using both positive- and negative-parity projected nucleon interpolating operators. Solving the GEVP yields optimal left and right eigenvectors that project onto the mixed-parity ground state.
Local vs. Global: By utilizing the local topological charge density rather than a global volume sum, the method avoids the severe statistical noise associated with global topological charge fluctuations.

Numerical Setup
The calculations were performed using 2+1 flavor dynamical domain wall fermions (DWF) with the Iwasaki gauge action. The study utilized three gauge ensembles with pion masses of 340, 420, and 576 MeV and a lattice spacing of $a \approx 0.11$ fm. The electro-quenched approximation was employed, where only valence quarks interact with the background field. The topological charge was constructed from gauge fields evolved with gradient flow.

Key Results

Statistical Significance: The authors obtained a statistically significant signal for the $\bar{\theta}$ -induced nEDM, overcoming the noise issues that plagued previous global-sum approaches.
Control of Systematics: The GEVP analysis successfully suppressed excited state contamination. The results were shown to be consistent across different nucleon interpolating operators and electric field strengths ( $|n_z|=1$ and $|n_z|=2$ ), indicating that higher-order corrections in the electric field are negligible.
Physical Point Extrapolation: Using a chiral extrapolation ansatz involving $m_\pi^2$ and $m_\pi^2 \log(m_\pi^2)$ , the results were extrapolated to the physical pion mass.
Final Value: The calculated nEDM sensitivity is:
$d_n = -0.0050(4)(8) \bar{\theta} \, e \cdot \text{fm}$
where the first uncertainty is statistical and the second is systematic (arising from fit range choices and the difference between electric field strengths).

Significance and Claims
The paper claims that this pioneering work demonstrates that the neutron EDM can be reliably determined from the local density of topological charge with robust control of systematic effects. The primary significance lies in:

Resolution of Statistical Noise: The method successfully bypasses the statistical noise inherent in global topological charge summations, enabling a clear signal extraction.
Systematic Control: By avoiding the $Q^2 \to 0$ extrapolation and using variational analysis to handle parity mixing, the method reduces specific systematic uncertainties associated with previous form-factor-based approaches.
Strong-CP Constraint: Comparing the calculated value with current experimental bounds ( $d_n \sim 10^{-26} e \cdot \text{cm}$ ), the authors derive a constraint on the QCD theta angle of $|\bar{\theta}| \lesssim 10^{-11}$ . This confirms the existence of the strong-CP problem in QCD.
Extensibility: The authors state that this framework can be directly extended to study other CP-violating interactions, such as the Weinberg three-gluon operator and isosinglet four-quark operators, where global summation methods typically fail due to noise.

The authors remain modest regarding remaining uncertainties, acknowledging that finite volume effects, discretization effects, and the extrapolation from the lowest pion mass ($340$ MeV) to the physical point require further investigation with lighter masses and finer lattice spacings.

More like this