Chiral Anomaly of Kogut-Susskind Fermion in the (3+1)-dimensional Hamiltonian formalism

This paper investigates the chiral anomaly of Kogut-Susskind fermions in a (3+1)-dimensional Hamiltonian formalism by defining a non-onsite axial charge that, while generally non-conserved, satisfies the continuum anomalous conservation law for specific $U(1)$ link configurations, as confirmed by numerical simulations.

The Gist

Imagine you are trying to build a perfect digital simulation of the universe, specifically the world of subatomic particles called fermions (like electrons). In the real world, these particles have a mysterious property called "chirality" or "handedness." Think of it like your hands: your left hand is a mirror image of your right, but you can't rotate your left hand to make it look exactly like your right. In physics, this "handedness" is crucial because the laws of nature treat left-handed and right-handed particles differently.

The problem is that when physicists try to simulate this on a computer grid (a lattice), a famous rule called the Nielsen-Ninomiya theorem says you can't just put a single "handed" particle on the grid without accidentally creating a bunch of unwanted "ghost" particles that mess up the simulation.

This paper is about a clever workaround using a specific type of grid particle called Kogut-Susskind (or Staggered) fermions. Here is the story of what the authors discovered, explained simply:

1. The "Checkerboard" Trick

Imagine a 3D checkerboard. Instead of putting a particle on every single square, the Kogut-Susskind method puts particles in a staggered pattern. It's like a dance where partners are always one step apart.

The authors realized that on this specific grid, there is a special move called a "Diagonal Shift."

• The Analogy: Imagine you are walking through a city grid. Usually, you walk North, then East, then South. But this special move tells you to take a giant diagonal step that jumps over the corners of the blocks.
• The Discovery: The authors found that if you perform this specific diagonal jump, it acts exactly like flipping the "handedness" (chirality) of the particle. It's a discrete, digital version of the continuous "chiral symmetry" found in the real world.

2. The "Ghost Charge" (The Anomaly)

In the real world, there is a phenomenon called the Chiral Anomaly. It's like a conservation law that usually says "Handedness is preserved," but under certain conditions (like strong magnetic and electric fields), it breaks down. The universe "leaks" handedness.

The authors asked: Does this digital "Diagonal Shift" trick also leak handedness when we turn on digital magnetic fields?

• The Challenge: In their digital world, the "Diagonal Shift" usually gets messy when you add magnetic fields. It's like trying to do a perfect dance step while the floor is shaking.
• The Solution: They found a very specific, special arrangement of magnetic and electric fields (like a perfectly tuned magnetic "wind" blowing diagonally) where the dance step works perfectly again.
• The Result: When they ran the numbers, they found that the "handedness charge" in their simulation did leak exactly the way it should in the real world. The digital simulation reproduced the famous "anomaly" equation of the real universe.

3. The "Onsager Algebra" (The Rulebook)

The paper also dives into the mathematical "rulebook" (algebra) that governs these particles.

• The Analogy: Think of the particles as players in a game with a strict set of rules (the Onsager algebra). Some players have "quantized" scores (whole numbers only), while others have "non-quantized" scores (fractions allowed).
• The Discovery: The authors defined a new "handedness score" (called $Q_A$) based on their diagonal shift. They found that this score is non-quantized (it can be a fraction) and it plays nicely with the other rules of the game. It acts like a "central charge"—a special, constant value that sits at the heart of the system's symmetry, much like the center of a spinning top.

4. Why This Matters

Why should a non-physicist care?

• Better Simulations: Currently, simulating the strong nuclear force (which holds atoms together) is incredibly hard because of these "handedness" problems. If we can define a "chiral symmetry" that works perfectly on a computer grid, we can build much more accurate simulations of the universe.
• The "Theta" Term: The authors mention that this could help simulate "theta terms," which are mysterious parameters in physics that might explain why the universe has more matter than antimatter.
• The Bridge: This paper builds a bridge between the messy, discrete world of computer grids and the smooth, continuous world of real physics. They showed that even on a pixelated grid, you can capture the deep, magical laws of the universe.

Summary

In short, the authors took a specific type of digital particle, found a special "diagonal dance move" that flips its handedness, and proved that when you add magnetic fields, this digital system behaves exactly like the real universe does—leaking handedness in a very specific, predictable way. This is a major step toward creating perfect computer models of the fundamental forces of nature.

Technical Summary

1. Problem Statement

The paper addresses the challenge of defining and understanding chiral symmetry and the associated chiral anomaly for Kogut-Susskind (KS) fermions (also known as staggered fermions) within the Hamiltonian formalism in (3+1) dimensions.

While significant progress has been made in (1+1) dimensions regarding discrete shift symmetries acting as chiral transformations, the extension to (3+1) dimensions presents complexities:

• In (1+1)D, specific shift operators generate quantized axial charges that form an Onsager algebra.
• In (3+1)D, previous works identified quantized charges related to shift symmetries, but the relationship between these charges, the chiral anomaly, and the continuum limit (specifically for a two-flavor theory) remained to be fully clarified.
• The authors aim to define a non-on-site axial charge in (3+1)D that reproduces the correct chiral anomaly equation under background gauge fields, analogous to the (1+1)D case but distinct from previously defined quantized charges.

2. Methodology

The authors employ a combination of analytical lattice field theory and numerical simulations:

• Hamiltonian Formalism: They work with continuous time and discrete space. The KS fermion Hamiltonian is defined on a 3D lattice with shift operators $S_k$ ($k=1,2,3$) in each spatial direction.
• Definition of Chiral Operator: They define a diagonal shift operator $\Gamma = \sqrt{-1} S_1 S_2 S_3$. Unlike the shift operators in (1+1)D, this operator is unitary (not Hermitian).
• Charge Construction:
  • They decompose $\Gamma$ into Hermitian ($\Gamma + \Gamma^\dagger$) and anti-Hermitian ($\Gamma - \Gamma^\dagger$) parts.
  • The Hermitian part defines the axial charge $Q_A = \frac{1}{2} \sum_x \chi^\dagger (\Gamma + \Gamma^\dagger) \chi$. This is a non-on-site, non-quantized charge that commutes with the vector charge $Q_V$.
  • They analyze the algebraic structure of these charges within the generalized Onsager algebra generated by shift operators acting on Majorana components of the KS fermion.
• Gauge Field Coupling:
  • They introduce background $U(1)$ link variables (gauge fields). Generally, the chiral operator $\Gamma$ does not commute with the Hamiltonian in the presence of gauge fields.
  • They identify a specific class of link configurations (uniform magnetic and electric fields aligned diagonally or along specific axes) where the commutativity $[H, \Gamma] = 0$ is restored.
• Numerical Simulation:
  • They perform numerical diagonalization of the one-particle Hamiltonian on a lattice ($N=8$).
  • They simulate an adiabatic evolution of the electric field to track the time evolution of the vacuum expectation value (VEV) of the axial charge $\langle Q_A(t) \rangle$.
  • They compare the lattice results against the continuum prediction for the chiral anomaly equation: $\partial_\mu j^\mu_A = -\frac{E \cdot B}{2\pi^2}$.

3. Key Contributions

A. Definition of a New Axial Charge ($Q_A$)

The authors define a new axial charge $Q_A$ based on the Hermitian part of the diagonal shift operator $\Gamma$.

• Properties: $Q_A$ is non-on-site, non-quantized, and commutes with the vector charge $Q_V$.
• Distinction: Unlike the quantized charges $Q_{ST}$ found in previous literature (which generate the Onsager algebra), $Q_A$ acts as a central charge of the generalized Onsager algebra in (3+1)D. It is not related to the quantized charges in the same way the (1+1)D analog $\tilde{Q}_A$ is related to $Q_1 + Q_{-1}$.
• Continuum Limit: In the continuum limit, the diagonal shift operator $\Gamma$ converges to $\gamma_5 \otimes \mathbb{1}$ on the taste basis, corresponding to the chiral operator for a two-flavor Dirac fermion.

B. Algebraic Structure (Onsager Algebra)

The paper constructs the generalized Onsager algebra for the (3+1)D KS system.

• They demonstrate that the regularized chiral charge $Q_{reg}$ belongs to this algebra.
• Crucially, they show that $Q_A$ commutes with all elements of the Onsager algebra, identifying it as a central charge. This supports the argument that the previously identified quantized charge $Q_{ST}$ is not the chiral charge responsible for the anomaly.

C. Restoration of Symmetry with Specific Gauge Fields

The authors prove that while the chiral symmetry is generally broken by arbitrary gauge fields, it is preserved for a specific class of $U(1)$ link configurations carrying non-trivial uniform magnetic and electric fields. This allows for a rigorous definition of the vacuum state and the calculation of the anomaly.

D. Numerical Verification of the Anomaly

Through numerical simulations, the authors verify that the time derivative of the vacuum expectation value of $Q_A$ satisfies the anomalous conservation law:
$$\frac{d}{dt} \langle Q_A(t) \rangle = - \sum_x \frac{E_i B_i}{2\pi^2}$$

• The results match the continuum prediction for a two-flavor Dirac fermion system (a factor of 2 difference from the single-flavor case).
• They also analyze the anti-Hermitian part ($\tilde{Q}$), showing its expectation value corresponds to the imaginary part of the continuum current, further confirming the consistency of the lattice regularization.

4. Results

• Anomaly Equation: The numerical data for $\langle Q_A(t) \rangle$ under adiabatic electric field evolution exhibits discontinuities at times corresponding to zero-mode crossings. The slope of the linear regions between discontinuities perfectly matches the theoretical prediction for the chiral anomaly in a two-flavor theory.
• Continuum Convergence: The lattice definitions of the chiral charge density and current converge to the continuum expressions involving $\gamma_5$ and the electromagnetic field tensor.
• Algebraic Insight: The identification of $Q_A$ as a central charge clarifies the structure of chiral symmetries in KS fermions, distinguishing the "true" chiral charge from other shift-generated charges.

5. Significance

• Hamiltonian Lattice QCD: This work provides a robust framework for studying chiral symmetry and anomalies in the Hamiltonian formalism, which is crucial for real-time dynamics and quantum simulation (e.g., on quantum computers) where Euclidean path integrals are not directly applicable.
• Two-Flavor Consistency: It confirms that KS fermions in (3+1)D correctly reproduce the physics of two-flavor Dirac fermions, resolving ambiguities about the "taste" structure in the chiral limit.
• Quantum Simulation: By constructing a discrete chiral operator that explicitly possesses chiral symmetry in specific gauge configurations, the authors suggest a pathway to constructing Hamiltonians for Quantum Electrodynamics (QED) that explicitly preserve chiral symmetry. This could significantly improve the convergence of numerical calculations to the continuum limit, aiding in the simulation of physical phenomena involving $\theta$-terms and topological effects.
• Theoretical Clarification: The paper resolves the relationship between shift symmetries, Onsager algebras, and chiral anomalies in higher dimensions, correcting previous interpretations that conflated quantized shift charges with the chiral charge.