⚛️ lattice

The Neutron Electric Dipole Moment from Lattice QCD using a Background Electric Field

Using 2+1 flavor domain wall fermion lattice QCD ensembles and a novel method involving a background electric field with single time-slice topological charge sampling to mitigate statistical noise and excited-state contamination, the authors calculate the neutron electric dipole moment to be $d_n=-0.0050(4)^\text{stat}(8)^\text{sys}\bar{\theta}$ $e$ fm after extrapolation to the physical point.

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

Published 2026-07-14

📖 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

Imagine the neutron not as a tiny, solid marble, but as a squishy, invisible cloud of quarks and gluons, buzzing with energy. Now, imagine that this cloud has a secret: it might have a tiny, permanent "lopsidedness" in its electric charge. If it does, it would have an Electric Dipole Moment (nEDM). Think of this like a tiny internal compass needle that points in a specific direction, even when the neutron is just sitting still.

Why do we care? Because finding this lopsidedness is like finding a fingerprint of a cosmic crime scene. It would tell us that the universe has a hidden preference for "left" over "right" in how it breaks the rules of physics (a violation called CP violation). The Standard Model of physics predicts this lopsidedness should be so incredibly tiny—about a million times smaller than our best microscopes can see—that it's basically zero. But if we find it, it means there's a whole new layer of reality we haven't discovered yet.

The Problem: The Ghost in the Machine

To measure this, the scientists in this paper decided to build a virtual universe inside a supercomputer. They used a method called Lattice QCD, which is like drawing a 3D grid (a lattice) and simulating the rules of the strong force that holds quarks together.

But here's the tricky part: The neutron is a quantum object. It doesn't just sit there; it vibrates, wiggles, and constantly pops in and out of existence with other particles. When you try to measure its "lopsidedness," you get a signal that is drowning in noise. It's like trying to hear a whisper in a hurricane.

Previous attempts to find this whisper used a method that was like trying to measure the wind by looking at the entire ocean at once. They summed up the "topological charge" (a fancy way of counting how twisted the quantum fields are) over the whole computer simulation. The problem? This global sum was so noisy that the signal was lost. It was like trying to find a specific person in a stadium by counting every single person in the building at the same time; the crowd noise was too loud.

The New Trick: The Single-Page Snapshot

The authors of this paper came up with a clever new way to listen to the whisper. Instead of looking at the whole stadium, they decided to take a snapshot of just one row of seats (a single time-slice).

They realized that if you apply a gentle, uniform electric field to their virtual neutron, the energy of the neutron shifts slightly. This shift is directly related to the electric dipole moment. But to get the right answer, they needed to measure the "topological charge" right where the neutron is, not everywhere else.

They used a mathematical tool called the Feynman-Hellmann theorem, which is like a rule that says: "If you know how the energy of a system changes when you tweak a knob, you know exactly what the system is doing." By tweaking the "theta-term" (a parameter that controls the CP violation) and looking at the energy shift, they could calculate the nEDM.

The "Ghost" Problem: Excited States

There was another hurdle. The neutron isn't just one thing; it has "excited states." Imagine the neutron as a guitar string. It can vibrate in its lowest, calmest note (the ground state), but it can also vibrate in higher, noisier notes (excited states).

In their simulations, the "noisy notes" were contaminating the signal. It was like trying to hear the guitar's main note, but the higher, shriller notes were screaming over it. The paper highlights a critical flaw in the traditional approach: using the standard "positive parity" operator to create the neutron was subject to massive excited-state contamination. This meant that if you used this old method, the results looked messy and inconsistent, changing wildly depending on which mathematical "guitar string" (interpolating operator) you used to describe the neutron.

To fix this, the team used a sophisticated math trick called the Generalized Eigenvalue Problem (GEVP). Think of this as a super-powered noise-canceling headphone. Instead of just listening to the sound, the math analyzes the pattern of the noise and the signal together, allowing them to isolate the pure, calm "ground state" of the neutron and ignore the screaming excited states. Crucially, once they used this new method to filter out the contamination, they found that different types of neutron operators (whether they were "covariant" or "non-covariant") all agreed on the same answer.

The Results: A Precise Simulation, Not a Final Measurement

After all this hard work, filtering out the noise, and isolating the ground state, they got a result.

They found that the neutron's electric dipole moment is:
dn = −0.0050(4)stat(8)sys ¯θ e fm

Let's break that down:

−0.0050: This is the value derived from their specific simulation setup. The negative sign indicates the direction of the lopsidedness in this model.
¯θ: This is a variable representing the strength of the CP violation. The result is proportional to this.
e fm: This is the unit (electron charge times femtometers).
(4)stat and (8)sys: These are the "error bars." The first number (4) is the statistical uncertainty (how much the data wiggles). The second (8) is the systematic uncertainty (how much the method might be slightly off).

Crucially, the authors emphasize that this is a simulation result, not a final measurement of the real world. The value they obtained is consistent with zero within the current uncertainties, and it represents a step toward the physical point rather than the final answer. They used "2+1 flavor domain wall fermion ensembles" (a specific type of virtual matter) with a lattice spacing of 0.11 fm and pion masses of 340, 420, and 576 MeV.

What They Ruled Out (and What They Didn't)

The paper explicitly argues against the idea that you can just use the "global" topological charge (the whole stadium count) to get a clean signal. They showed that this old method is too noisy and gets contaminated by excited states. They also showed that using different types of "neutron operators" (different ways of mathematically describing the neutron) gives the same answer only if you properly filter out the excited states. If you don't filter them, the results look messy and inconsistent.

They also tested two different ways of defining "topological charge": one based on the gluon fields (the "glue" holding quarks together) and one based on the quark density itself. They found that while the global versions of these two definitions matched, the local versions behaved differently depending on the mass of the quarks. However, once they used their new "noise-canceling" method, the results from both definitions converged to the same answer.

How Sure Are They?

The authors are confident in their method and their simulation results, but they are very clear that this is not the final word on the real world.

They have simulated the nEDM with high precision using their new technique.
They have measured the value within their virtual universe, finding a result that is currently consistent with zero.
They have not yet accounted for all the real-world errors, like the specific size of their virtual grid (discretization) or the fact that their pion masses were still a bit heavier than in the real universe. They explicitly state that these "conventional systematic errors" will be addressed in future work.

So, while they haven't "solved" the mystery of the neutron's lopsidedness for the real universe yet, they have built a much better, quieter microphone to listen for it. They've proven that with the right math and a single-time-slice snapshot, the whisper of the neutron's electric dipole moment can finally be heard above the cosmic hurricane.

Problem Statement

The neutron electric dipole moment (nEDM) is a sensitive probe for CP violation beyond the Standard Model. While the Standard Model predicts an nEDM far below current experimental limits, the QCD $\theta$ -term (the "strong CP problem") could generate a measurable signal. Calculating the $\theta$ -induced nEDM from first principles using Lattice QCD is challenging due to two primary issues:

Statistical Noise: Traditional methods using the global topological charge ( $Q_{top} = \int d^4x \, q_{top}$ ) suffer from large statistical fluctuations in nucleon correlation functions.
Excited State Contamination: Extracting the nEDM requires isolating the ground state matrix element. In the presence of a background electric field, the nucleon ground state mixes with negative-parity excited states and its own parity partner. Traditional positive-parity projected operators are subject to significant excited-state contamination, leading to unreliable results.
Operator Dependence: Previous studies suggested that the choice of nucleon interpolating operators (covariant vs. non-covariant under chiral transformations) might significantly affect the extracted nEDM values.

Methodology

The authors employ a background field method combined with a novel approach to topological charge and excited state isolation:

Background Field Method: Instead of extrapolating the electric dipole form factor $F_3(Q^2)$ to $Q^2=0$ , the nEDM is extracted from the energy shift of the neutron in a static, uniform external electric field ( $\vec{E}$ ). The field is introduced via a $U(1)$ gauge link modification on the lattice, analytically continued to imaginary time to preserve periodic boundary conditions.
Local Topological Charge & Feynman-Hellmann Theorem: Motivated by the Feynman-Hellmann theorem, the authors calculate the energy shift induced by the $\theta$ -term. They propose sampling the local topological charge density ( $q_{top}$ ) on a single time-slice rather than the global topological charge. This reduces the stochastic noise associated with "disconnected" diagrams. The nEDM is related to the matrix element of this local charge density between nucleon states deformed by the electric field:
$d_n/\bar{\theta} = \frac{1}{E_z} \langle N_0^L | q_{top}(0) | N_0^R \rangle_{E_z}$
Generalized Eigenvalue Problem (GEVP): To address excited state contamination and the mixing of positive- and negative-parity states induced by the electric field, the authors construct a correlation matrix using both positive-parity ( $N$ ) and negative-parity ( $N^*$ ) interpolating operators. They solve a non-Hermitian generalized eigenvalue problem to isolate the true ground state eigenvectors ( $v_{L,0}, v_{R,0}$ ) in the presence of the anti-Hermitian perturbation caused by the electric field.
Computational Setup: The calculation uses $2+1$ flavor domain wall fermion ensembles with fixed lattice spacing $a \approx 0.11$ fm and pion masses of 340, 420, and 576 MeV. To improve statistical precision, the authors utilize All Mode Averaging (AMA) and Low Mode Averaging (LMA), leveraging the low-lying eigenmodes of the Dirac operator which are crucial for topological charge sensitivity.
Topological Charge Definitions: The study compares two definitions of topological charge:
1. Gluonic: Constructed from gradient-flowed gauge fields.
2. Fermionic: Constructed from the pseudoscalar quark density ( $\bar{\psi}\gamma_5\psi$ ) via the ABJ anomaly, using low-mode approximations.

Key Contributions and Results

Validation of Local Topological Charge: The authors demonstrate that using the local topological charge on a single time-slice, combined with the GEVP, dramatically improves the statistical precision of the $\theta$ -induced nEDM compared to using the global topological charge.
Excited State Removal: By solving the non-Hermitian GEVP, the authors successfully isolate the ground state, removing the large excited-state contamination that plagued traditional positive-parity projections.
Operator Independence: The study compares four different nucleon interpolating operators (two non-covariant and two chiral-covariant). The results show that consistent nEDM values are obtained regardless of whether the operators are covariant or non-covariant under chiral transformations, provided the operators have good overlap with the relevant nucleon eigenstates in the background field.
Topological Charge Consistency:
- Results based on the gluonic topological charge (with varying gradient flow times) are consistent with each other.
- Results based on the pseudoscalar quark density converge to the same value as the probing quark mass increases. For smaller quark masses, the pseudoscalar definition shows stronger $N\pi$ contamination.
Final Result: After extrapolating to the physical point (though full systematic errors like discretization and finite-volume effects are reserved for future work), the authors obtain:
$d_n = -0.0050(4)_{\text{stat}}(8)_{\text{sys}} \, \bar{\theta} \, e \cdot \text{fm}$
The systematic uncertainty includes excited-state effects (estimated via Euclidean-time fit variations) and dependence on the electric field strength.

Significance and Claims

The paper claims that this methodology provides a more statistically precise and systematically reliable determination of the nEDM compared to previous lattice approaches. By shifting from global to local topological charge and rigorously treating the non-Hermitian nature of the background field problem via GEVP, the authors overcome the "statistical noise" and "excited state contamination" barriers that have historically hindered direct lattice calculations of the nEDM at physical quark masses.

The authors explicitly state that while their result is consistent with other phenomenological estimates, the primary contribution is the methodological framework that enables the extraction of the nEDM with controlled uncertainties. They note that conventional systematic errors (discretization, finite-volume, and chiral extrapolation) are not fully addressed in this work and will be the subject of future publications. The paper does not claim to resolve the strong CP problem or predict a specific experimental outcome but rather establishes a robust computational tool for future precision studies.