Understanding the structure of nucleon excitations from their wavefunctions

This paper utilizes approximately 4000 lattice QCD propagators on a heavy pion ensemble to analyze relativistic nucleon wavefunctions, revealing that the observed node structure in the spectrum arises from two distinct mechanisms: "superposition nodes" created by combining interpolating fields and novel "built-in nodes" inherent to the s-wave Dirac components of individual fields.

The Gist

Imagine the proton (the particle inside an atom's nucleus) not as a solid marble, but as a tiny, vibrating cloud of three quarks dancing together. For decades, physicists have known what these particles weigh and how they spin, but they've struggled to see how they move and arrange themselves inside. It's like knowing a song's pitch and tempo, but not being able to see the sheet music.

This paper is like finally turning on the lights to see the sheet music of the proton's excited states. The authors used a supercomputer to simulate the universe's fundamental forces and took a high-resolution "snapshot" of these dancing quarks. Here is what they found, explained simply:

1. The Setup: A Digital Sandbox

To study these particles, the researchers built a digital universe on a grid (a lattice). They didn't just look at one type of particle; they looked at the proton and its "excited" versions (like a guitar string vibrating in different patterns).

They used two different "flashlights" (mathematical tools called interpolating fields) to shine on the quarks:

• Flashlight A (The Standard One): This sees the quarks moving slowly, like a classic cartoon character.
• Flashlight B (The Relativistic One): This sees the quarks moving at near-light speeds. It's a more complex tool that "disappears" if you try to look at slow-moving particles, but it reveals hidden details when things get fast.

2. The Big Discovery: The "Node" Mismatch

In quantum mechanics, a node is a spot where the probability of finding a particle drops to zero. Think of a vibrating guitar string: the ends are fixed, but there's a point in the middle that stays still while the rest vibrates. That still point is a node.

Usually, physicists expect all parts of a particle's wave (the "upper" and "lower" parts of the mathematical description) to have the same number of nodes. If the string has one bump, the whole string has one bump.

But this paper found something weird:
For certain excited protons, the "upper" part of the wave had a node (a bump), but the "lower" part did not. It was as if the top half of the guitar string was vibrating in a complex pattern, while the bottom half was just a simple curve.

They found two types of these "missing" bumps:

• The "Superposition" Node: This happens when two different wave patterns cancel each other out in the middle, creating a gap. This affects the whole particle evenly.
• The "Built-in" Node: This is the surprise. It's a gap that is hard-wired into the math of one specific flashlight tool. It only appears in certain parts of the wave, creating that mismatch.

3. The "Role Reversal" Analogy

Here is the most fascinating part. The authors found that the two flashlights (A and B) play opposite roles depending on whether the proton is spinning "up" or "down" (positive or negative parity).

• For "Positive" Protons: The "Standard" flashlight (A) shows simple waves, but the "Relativistic" flashlight (B) has the built-in node. If the proton is mostly made of the "Relativistic" type, you see that weird mismatch.
• For "Negative" Protons: The roles flip! Now the "Standard" flashlight (A) has the built-in node, and the "Relativistic" one (B) is simple.

It's like a pair of twins who wear different colored shirts. In the morning, Twin A wears red and Twin B wears blue. But in the afternoon, they swap shirts. The paper figured out exactly why they swap and how it changes the shape of the dance.

4. Why Does This Matter?

Imagine trying to understand a complex machine by only looking at its shadow. For a long time, we only saw the "shadow" (the mass and spin) of these particles.

By mapping out the actual "shape" of the wavefunctions (the dance moves), the authors found that:

1. The "Built-in" nodes explain the mass: The extra bump in the wave pushes the quarks further out, which costs more energy. This explains why some excited protons are heavier than others.
2. It confirms our models: They compared their computer results to the "MIT Bag Model" (a theory where quarks are trapped in a tiny balloon). The math showed that the "balloon" naturally creates these extra bumps in specific ways, confirming that their computer simulation is accurate.

The Takeaway

This paper is like a detective story where the clues were hidden in the "noise" of the data. By using a clever mix of two different mathematical tools and looking at the data from every angle, they discovered that the internal structure of the proton is more complex than we thought.

They found that nature uses a "switch" (the $\gamma_5$ operator in the math) to decide which part of the particle gets the extra "bump" (node). This helps us understand not just what the proton is, but how the fundamental forces of the universe stitch it together. It's a step toward understanding the very fabric of matter.

Technical Summary

1. Problem Statement

The paper addresses the fundamental challenge of understanding the internal structure of nucleon excitations (resonances) within Quantum Chromodynamics (QCD). While experimental data provides masses and quantum numbers, the underlying wavefunction structure—specifically the node structure of quark distributions—remains theoretically complex.

• The Gap: Previous lattice QCD studies focused primarily on the upper (large) Dirac components of positive-parity nucleon wavefunctions using standard scalar-diquark interpolating fields. These studies confirmed a "standard" node structure (zero nodes for the ground state, one node per excitation).
• The Question: How do relativistic effects, specifically the inclusion of pseudoscalar-diquark interpolating fields (which vanish in the non-relativistic limit) and the lower (small) Dirac components, alter the node structure of both positive- and negative-parity spectra? Do the upper and lower components always share the same node count, or can they differ?

2. Methodology

The authors employed Lattice QCD calculations using a variational analysis approach to extract wavefunctions.

• Ensemble: They utilized the heaviest available (2+1)-flavour PACS-CS dynamical ensemble with a pion mass of $m_\pi \simeq 702$ MeV and lattice spacing $a \simeq 0.0907$ fm. This heavy mass regime is chosen to minimize finite-volume effects and isolate quark-model-like states from multi-particle scattering states.
• Interpolating Fields: A 10$\times$10 correlation matrix was constructed using two distinct types of local interpolating fields for the proton, each at five levels of gauge-invariant Gaussian smearing (12, 24, 48, 96, 192 sweeps):
  1. $\chi_1$ (Scalar-diquark): Standard field with a non-relativistic reduction overlapping strongly with the SU(6) quark model.
  2. $\chi_2$ (Pseudoscalar-diquark): Contains a $\gamma_5$ matrix acting on the uncontracted quark field. It vanishes in the non-relativistic limit and couples strongly to higher-energy states.
• Wavefunction Extraction:
  • The calculation was performed in Landau gauge.
  • To extract the full relativistic wavefunction, the authors shifted the spatial position of the $d$-quark relative to two fixed $u$-quarks at the origin.
  • They exploited parity and spin mismatches between the source and sink Dirac indices. By selecting specific combinations of source ($\alpha$) and sink ($\beta$) indices, they isolated specific relativistic components (upper/lower, s-wave/p-wave) and angular momentum dependencies.
• Analysis:
  • Variational Analysis: Solved the Generalized Eigenvalue Problem (GEVP) to extract energy eigenstates and eigenvectors (coefficients of interpolating fields).
  • Visualization: Generated volume and surface renderings of wavefunction amplitudes.
  • Quantitative Analysis: Calculated radial wavefunctions by averaging over lattice sites at fixed radii after dividing out the predicted spherical harmonics ($Y_{\ell}^{m}$).

3. Key Contributions

1. First Full Relativistic Wavefunction Study: This is the first study to visualize and quantify both the upper and lower Dirac components for negative-parity nucleon excitations, as well as the lower components for positive-parity states.
2. Discovery of "Node Mismatch": The paper identifies a novel phenomenon where the number of nodes in the s-wave components differs from the p-wave components within the same energy eigenstate.
3. Identification of Two Node Types:
   • Superposition Nodes: Formed by the linear combination (interference) of different interpolating fields (e.g., $\chi_1$ and $\chi_2$). These appear in all spinor components.
   • Built-in Nodes: Intrinsic nodes present in the wavefunction of a single interpolating field due to its specific spin-flavour structure (specifically the presence or absence of $\gamma_5$). These appear only in specific components (s-wave or p-wave).
4. Role Reversal of Interpolators: The study demonstrates that $\chi_1$ and $\chi_2$ play inverse roles for positive- and negative-parity spectra regarding which field carries the "built-in" node.

4. Key Results

A. Positive-Parity Spectrum

• Ground State ($1^+$): Zero nodes in all components.
• First Excitation ($2^+$): One node in both upper (s-wave) and lower (p-wave) components. This node arises from the superposition of $\chi_1$ and $\chi_2$ fields with opposite signs.
• Second Excitation ($3^+$): A node mismatch is observed. There is one node in the upper (s-wave) components, but zero nodes in the lower (p-wave) components.
  • Cause: This state is dominated by the $\chi_2$ interpolator. The $\chi_2$ field possesses a "built-in" node in its upper (s-wave) Dirac components but none in its lower components.
• Higher States: The pattern repeats. States dominated by $\chi_1$ show matching nodes; states dominated by $\chi_2$ show a mismatch (extra node in s-wave).

B. Negative-Parity Spectrum

• Ground State ($1^-$): Zero nodes in all components.
• First Excitation ($2^-$): One node in both upper (p-wave) and lower (s-wave) components.
• Second Excitation ($3^-$): A node mismatch is observed. There is one node in the lower (s-wave) components, but zero nodes in the upper (p-wave) components.
  • Cause: This state is dominated by the $\chi_1$ interpolator. The $\chi_1$ field possesses a "built-in" node in its lower (s-wave) components.
• Pattern: The roles are reversed compared to positive parity. $\chi_1$-dominated states exhibit the mismatch (extra node in s-wave), while $\chi_2$-dominated states have matching nodes.

C. Theoretical Consistency

• The authors compared their findings with the MIT Bag Model and the Dirac equation with a linear scalar potential ($S(r) = \sigma r$).
• Both models predict that for negative-parity solutions, the lower (s-wave) components naturally possess an extra node compared to the upper (p-wave) components. This confirms that the "built-in" nodes observed in the lattice data are a fundamental feature of relativistic confinement, not an artifact of the lattice method.

5. Significance

• Refining the Quark Model: The results challenge the simple non-relativistic quark model view where nodes are strictly determined by the principal quantum number $n$. Instead, the node structure is a complex interplay of relativistic spin-flavour dynamics and the specific interpolating fields used to construct the state.
• Understanding Mass Splittings: The node mismatch (where the s-wave component is pushed to larger radii due to the extra node) provides a mechanism for the mass splitting between closely spaced resonances (e.g., the $N(1535)$ and $N(1650)$). The extra node in the s-wave component increases the energy of the state.
• Methodological Advancement: The paper establishes a robust framework for extracting and visualizing full relativistic wavefunctions, including the often-neglected lower components and negative-parity states. This opens the door for future studies on non-zero momentum and the separation of quark flavors.
• Fundamental QCD Insight: The identification of "built-in" nodes linked to the $\gamma_5$ structure of interpolators offers a deeper understanding of how chiral symmetry breaking and confinement manifest in the spatial distribution of quarks.

In summary, the paper provides a comprehensive picture of the nucleon spectrum, revealing that the wavefunction node structure is not uniform across all Dirac components and is fundamentally dictated by the relativistic properties of the interpolating fields used to generate the states.