⚛️ lattice

Perturbative quantum electrodynamics with generalized domain wall fermions

This paper derives the expansion of the generalized domain-wall fermion Dirac operator including electromagnetic corrections up to $\mathcal{O}(e^2)$ , identifying necessary contact terms and new local operator insertions (seagull vertices) to enable gauge-covariant lattice computations of radiative corrections for hadronic processes involving chiral fermions.

Original authors: Matteo Di Carlo, Antonin Portelli

Published 2026-07-24

📖 8 min read🧠 Deep dive

CC BY 4.0

Original authors: Matteo Di Carlo, Antonin Portelli

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.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 universe as a giant, cosmic LEGO set. At the very bottom, the tiniest bricks are particles called quarks, which snap together to build protons and neutrons—the stuff that makes up your body, the stars, and everything you see. But these bricks don't just sit there; they are constantly buzzing with energy and interacting with invisible forces. One of the most important forces is electromagnetism, the same thing that makes magnets stick to your fridge or gives you a static shock when you touch a doorknob.

For decades, scientists have been trying to build a perfect digital model of this LEGO universe using supercomputers. They call this "lattice QCD." It's like trying to simulate a complex video game, but instead of pixels, the game world is made of a grid of points in space and time. The goal is to predict exactly how these particles behave. However, there's a catch: the rules of the game are incredibly tricky. The particles have a property called "chirality," which is a fancy way of saying they have a specific "handedness" (like a left hand vs. a right hand) that affects how they interact. When scientists tried to add the rules of electromagnetism to their digital models, they found that the "left-handed" particles were behaving strangely, breaking the symmetry of the universe in their simulations. To fix this, they invented a clever trick called "Domain Wall Fermions," which essentially builds a tiny, extra dimension into their computer model to keep the particles' handedness intact. But even with this clever trick, there was a missing piece of the puzzle when it came to calculating how light interacts with these particles.

This paper is about finding that missing piece. The authors, M. Di Carlo and A. Portelli, have figured out exactly how to calculate the effects of electromagnetism on these special "domain wall" particles, specifically looking at how the particles react when you tweak the electric charge. They discovered that when you try to add these electromagnetic effects to the simulation, you can't just use the standard rules; you have to add some special "glue" terms to the equations. Without these extra terms, the simulation would break the fundamental laws of physics, specifically the rule that says electric charge must be conserved. They didn't just guess these terms; they mathematically derived them from scratch, showing exactly how the particles' behavior changes when you turn on the electromagnetic force. They also proved that their new formulas work by running computer tests that showed the math holds up perfectly, even when they pushed the electric charge to different levels. This work is a crucial update to the rulebook for anyone trying to simulate the universe with high precision, ensuring that future predictions about particle behavior are as accurate as possible.

The Story of the Missing Glue

Think of the universe's particles as dancers in a massive, choreographed ballet. For a long time, scientists have been trying to film this dance on a computer grid. But there's a problem: the dancers have a "handedness." Some are left-handed, some are right-handed, and the laws of physics say they must dance in a very specific way to keep the show running smoothly. If you try to simulate this dance on a computer, the grid itself can mess up the rhythm, making the dancers trip.

To fix this, physicists invented a special stage setup called "Generalized Domain Wall Fermions." Imagine building a second floor onto the dance stage. The dancers (the quarks) live on this 5th dimension, but we only watch the performance on the ground floor (our 4D world). This extra floor acts like a safety net, keeping the dancers' handedness perfect and preventing them from tripping over the grid.

Now, the scientists wanted to add a new element to the show: the electromagnetic force. This is like adding a spotlight that follows the dancers, changing how they move. In the real world, this is easy; in the computer simulation, it's a nightmare. When you try to turn on the spotlight in the simulation, the dancers start to stumble. Why? Because the way the dancers are connected to the extra floor (the 5th dimension) depends on the spotlight itself.

The authors of this paper realized that previous attempts to simulate this were missing a crucial step. They were trying to describe the dancers' movements using a simple script, but the script didn't account for the fact that the spotlight changes the dancers' connection to the safety net. When you expand the math to include the electric charge (the spotlight), you get extra terms that pop up out of nowhere. The authors call these "seagull vertices."

What are "Seagull Vertices"?
Imagine a seagull landing on a beach. It doesn't just walk; it lands, flaps its wings, and takes off all in one quick motion. In physics, a "seagull vertex" is a term in the math that represents a particle interacting with two photons (particles of light) at the exact same spot, instantly. It's a "contact term"—a sudden, direct interaction that happens without the particle traveling anywhere in between.

The paper shows that for these special 5D dancers, when you turn on the electromagnetic force, you must include these seagull terms. If you leave them out, the simulation loses its "gauge invariance." That's a fancy physics way of saying the simulation breaks the law of conservation of charge. It's like if your video game suddenly allowed you to create electricity out of thin air; the whole game would crash.

The Discovery and the Proof

The authors didn't just say, "Hey, we think we need these terms." They went into the math lab and derived them step-by-step.

The First Order (The Single Flash): They started by looking at what happens when the electric charge is turned on just a tiny bit (the first order). They re-derived a known result, showing that the dancers interact with the light in a way that creates a "conserved current." This is like the dancers moving in a perfect line, respecting the rules.
The Second Order (The Double Flash): Then, they looked at what happens when the charge is squared (the second order). This is where the magic happened. They found that because the dancers' connection to the 5th dimension depends on the light, new terms appear. These are the seagull vertices. They are local, meaning they happen right where the particle is, and they are essential to keep the math consistent.

The paper explicitly argues against the idea that you can just ignore these terms or replace them with a simple "renormalization factor" (a fudge factor to fix the numbers). The authors show that while some previous studies did this, it technically breaks the U(1) gauge invariance (the rule that charge is conserved) in the lattice theory. They insist that for a truly accurate simulation, you must include these specific, derived seagull terms.

How Sure Are They?

The authors are very confident in their results, but they back it up with more than just math. They wrote a computer program (using a library called Grid) to test their formulas.

The Test: They took a random snapshot of the universe (a gauge configuration) and a random starting point for a particle. They then calculated the particle's behavior using their new formulas.
The Comparison: To check if they were right, they also calculated the behavior by slightly changing the electric charge and measuring the difference (a numerical derivative). This is like checking your math by doing the problem a different way.
The Result: They compared the two methods. The difference between their formula and the numerical test was tiny and scaled perfectly as they made the electric charge smaller. The errors vanished as the charge went to zero.

This simulation proves that their derivation is correct. They didn't just find a theoretical possibility; they built a working tool that passes the test.

Why Does This Matter?

You might wonder, "So what? Who cares about these seagull vertices?"

The answer is precision. Scientists are trying to measure things like the magnetic moment of the muon (a heavy cousin of the electron) or the mass of protons and neutrons with incredible accuracy. They want to know if the Standard Model of physics is perfect or if there are cracks in it that point to new, undiscovered physics.

To find those cracks, the simulations need to be perfect. If you leave out the seagull vertices, your simulation has a hidden error. It's like trying to measure the weight of a feather on a scale that has a piece of gum stuck to it. You might get a number, but it's not the true number.

By providing the correct formulas for these seagull vertices, the authors have given the scientific community a cleaner, more accurate scale. This allows future experiments to calculate radiative corrections (the tiny tweaks caused by light) with full confidence. It paves the way for testing the "CKM unitarity," a fundamental rule about how particles mix and change, with even stricter precision. If the rules hold up under this new, more accurate scrutiny, we know our understanding of the universe is solid. If they don't, we might just be on the verge of discovering something entirely new.

In short, this paper is about fixing the glue in the universe's LEGO set so that when we build our digital models, they don't fall apart when we turn on the lights. It's a small piece of math with a huge impact on how we understand the fundamental building blocks of reality.

Technical Summary: Perturbative Quantum Electrodynamics with Generalized Domain Wall Fermions

Problem Statement
The inclusion of isospin-breaking effects and electromagnetic (QED) corrections in lattice QCD is essential for achieving sub-percent precision in hadronic observables. While perturbative approaches like the RM123 method allow for the reuse of existing isosymmetric QCD gauge ensembles by expanding correlation functions in powers of the electric charge $e$ , they require the identification of specific operator insertions at each order to preserve gauge invariance. For Wilson-like fermions and the standard Shamir implementation of Domain-Wall Fermions (DWF), these insertions (conserved currents at $O(e)$ and "seagull" vertices at $O(e^2)$ ) are well-established. However, for generalized DWF formulations (where the coefficients $b_s$ and $c_s$ in the 5D Dirac operator are arbitrary complex numbers), a rigorous derivation of the $O(e^2)$ local operator insertions has been missing. Consequently, previous studies using generalized DWF often neglected these terms or substituted the conserved current with a local current, resulting in a lattice theory that is not strictly $U(1)$ gauge invariant.

Methodology
The authors derive the perturbative expansion of the generalized DWF Dirac operator and the physical quark propagator up to $O(e^2)$ using a background-field approach. The methodology proceeds as follows:

Formalism: The generalized DWF action is defined via a 5-dimensional Dirac operator $D_{dw}$ acting on fields $\psi$ and $\bar{\psi}$ . The physical 4D quark fields $q$ and $\bar{q}$ are related to these 5D fields via projection operators $P$ and $R_5$ . Crucially, in the generalized formulation, the relation for the antiquark field $\bar{q}$ depends explicitly on the gauge links through the operator $D_-$ .
Expansion: The gauge links $U_\mu$ are promoted to $U(3)$ links via $U_\mu = \exp(ieA_\mu)U_\mu$ . The authors expand the physical propagator $S_q$ with respect to $e$ .
First-Order ( $O(e)$ ): Using the background-field method, the first derivative of the propagator is derived. This recovers the known conserved vector current insertion but reveals an additional contact term arising from the $e$ -dependence of the $D_-$ operator in the definition of $\bar{q}$ .
Second-Order ( $O(e^2)$ ): The second derivative is computed. This yields two types of contributions:
- Double insertions of the vector current (and single insertions of the second-order current).
- New local operator insertions (seagull vertices) arising from the expansion of the $D_-$ operator and the second derivative of the gauge links.
Ward Identities: The authors verify that these new terms satisfy the lattice Ward-Takahashi identities, ensuring that the contact terms cancel unphysical contributions when current positions coincide.
Numerical Validation: The derived expressions for the first and second-order currents were implemented in the Grid library. Validation tests compared the perturbative expansions against numerical finite-difference derivatives of the quark propagator for various DWF parameterizations (Shamir, scaled Shamir, and zMöbius), confirming the linear scaling of errors with $e^2$ .

Key Contributions

Explicit Derivation of Seagull Vertices: The paper provides the first explicit derivation of the $O(e^2)$ local operator insertions (seagull vertices) for generalized domain-wall fermions. These terms are shown to originate from the gauge-link dependence of the physical antiquark field definition in generalized formulations.
Re-derivation of $O(e)$ Terms: The first-order correction is re-derived, clarifying the structure of the conserved current and the additional contact term for generalized DWF.
Gauge Invariance Restoration: The work demonstrates that neglecting these $O(e^2)$ terms breaks $U(1)$ gauge invariance on the lattice. The derived vertices are necessary to restore this invariance and ensure the renormalizability of the lattice theory order-by-order.
Extension to Thermodynamics: The formalism is extended to the calculation of quark number susceptibilities at finite temperature by analytically continuing the electromagnetic potential to an imaginary chemical potential ( $eA_\nu \to i\mu_f \delta_{\nu,4}$ ).
Open-Source Implementation: The authors provide a public implementation of these currents in the Grid C++ library, including validation routines.

Results

The expansion of the physical quark propagator $S_q$ up to $O(e^2)$ is given by Eq. (90), which includes terms involving the conserved current $V^\mu$ , the first-order contact term $\bar{Q}^\mu_1$ , and the second-order seagull current $T^\mu$ and contact term $\bar{Q}^\mu_2$ .
The Ward-Takahashi identity for the domain-wall fermion currents (Eq. 97) explicitly links the divergence of the two-point current correlator to the expectation value of the seagull current $T^\mu$ , proving that the seagull term is required to cancel contact singularities.
Numerical tests on an $8^4$ lattice for Shamir, scaled Shamir ( $\alpha=2$ ), and zMöbius DWF actions show that the relative errors between the perturbative expansion and numerical derivatives scale linearly with $e^2$ and vanish as $e \to 0$ , validating the analytical derivation.

Significance
The paper claims that this work provides the necessary theoretical foundation for performing fully consistent, gauge-invariant lattice calculations of radiative corrections ( $O(\alpha_{em})$ ) to hadronic processes involving chiral fermions using generalized domain-wall fermions. By providing the missing $O(e^2)$ vertices, the authors enable future lattice studies to eliminate systematic errors associated with gauge non-invariance. This is particularly relevant for precision tests of CKM unitarity and the determination of hadronic contributions to the muon anomalous magnetic moment. The work does not propose new experimental setups but rather supplies the rigorous lattice machinery required to interpret high-precision experimental data.

The Story of the Missing Glue

The Discovery and the Proof

How Sure Are They?

Why Does This Matter?

More like this