pith. machine review for the scientific record. sign in

arxiv: 2604.10038 · v1 · submitted 2026-04-11 · ✦ hep-lat · hep-ph· nucl-th

Recognition: unknown

Understanding the structure of nucleon excitations from their wavefunctions

Derek B. Leinweber, Finn M. Stokes, Jackson A. Mickley, Waseem Kamleh

Authors on Pith no claims yet

Pith reviewed 2026-05-10 16:08 UTC · model grok-4.3

classification ✦ hep-lat hep-phnucl-th
keywords nucleon excitationslattice QCDwavefunction nodesinterpolating fieldsDirac componentsparity spectrumradial wavefunctionsquark model
0
0 comments X

The pith

Lattice wavefunctions of nucleon excitations display two node types, one built into single interpolating fields.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper computes position-space wavefunctions for the ground state and lowest excited nucleons using lattice QCD. It employs two local interpolating fields that carry the same quantum numbers but differ in spin-flavour structure, one vanishing in the nonrelativistic limit. By examining the d-quark distribution around two u quarks fixed at the origin, the calculation isolates nodes that arise only when the fields are combined from those already present inside the s-wave Dirac part of one field. A sympathetic reader would care because this separation shows how the nucleon spectrum is shaped by the choice of basic lattice operators rather than emerging solely from state mixing. The work is performed at zero momentum on a heavy-pion ensemble and covers both positive- and negative-parity channels.

Core claim

Relativistic wavefunctions of nucleon excitations are scrutinised to understand their node structure and the underlying role of local interpolating fields in generating the nucleon spectrum. Approximately 4000 propagators are employed on the heaviest PACS-CS ensemble at m_π ≃ 702 MeV. We examine the ground and four lowest-lying excited states at zero momentum for both positive- and negative-parity spectra, where the proton's d-quark wavefunction is calculated about the two u quarks at the origin. This is achieved using two local interpolating fields that each carry the quantum numbers of the nucleon but with differing spin-flavour structures, one of which vanishes in the nonrelativistic. We

What carries the argument

Local interpolating fields with differing spin-flavour structures that generate nucleon wavefunctions on the lattice, separating superposition nodes from built-in nodes inside s-wave Dirac components.

If this is right

  • The single-particle nucleon spectrum can be decomposed into contributions traceable to individual lattice operators.
  • Nodes observed in excited states arise either from linear combinations of interpolators or from intrinsic Dirac structure.
  • Radial wavefunction profiles quantify the location and character of both node types across parity channels.
  • Visual volume and surface renderings confirm the qualitative distinction between the two node classes.
  • The approach supplies a concrete link between quark-model node expectations and the output of lattice operators.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same distinction between built-in and superposition nodes could be tested in other baryons by applying analogous pairs of interpolators.
  • If the built-in nodes survive at lighter pion masses, they may point to a feature of QCD bound states independent of the heavy-pion approximation.
  • Improved interpolator design might exploit the built-in nodes to enhance overlap with specific radial excitations.
  • The method offers a template for examining whether similar intrinsic nodes appear in meson wavefunctions under comparable lattice setups.

Load-bearing premise

The two chosen local interpolating fields have enough overlap with the physical eigenstates and the heavy-pion ensemble yields node structures representative of the physical world.

What would settle it

Recompute the same wavefunctions on an ensemble with physical or near-physical pion mass and check whether the built-in nodes in the s-wave components persist unchanged.

Figures

Figures reproduced from arXiv: 2604.10038 by Derek B. Leinweber, Finn M. Stokes, Jackson A. Mickley, Waseem Kamleh.

Figure 1
Figure 1. Figure 1: FIG. 1. The extracted mass spectrum for the ground state [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The Euclidean-time evolution of the (1 [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. The RMS radius [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. As in Fig [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Volume and surface renderings of the (1 [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Volume and surface renderings of the (1 [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Volume and surface renderings of the (1 [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Scatter plots of the ( [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. The upper ( [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. As in Fig [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Volume and surface renderings of the (1 [PITH_FULL_IMAGE:figures/full_fig_p017_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. Volume and surface renderings of the (1 [PITH_FULL_IMAGE:figures/full_fig_p018_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13. The upper ( [PITH_FULL_IMAGE:figures/full_fig_p019_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: FIG. 14. As in Fig [PITH_FULL_IMAGE:figures/full_fig_p020_14.png] view at source ↗
Figure 16
Figure 16. Figure 16: FIG. 16. An illustration of how two normalised Gaussians, [PITH_FULL_IMAGE:figures/full_fig_p021_16.png] view at source ↗
Figure 15
Figure 15. Figure 15: FIG. 15. The eigenvector components [PITH_FULL_IMAGE:figures/full_fig_p021_15.png] view at source ↗
Figure 17
Figure 17. Figure 17: FIG. 17. The upper ( [PITH_FULL_IMAGE:figures/full_fig_p022_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: FIG. 18. As with Fig [PITH_FULL_IMAGE:figures/full_fig_p023_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: FIG. 19. The upper and lower radial wavefunctions from the [PITH_FULL_IMAGE:figures/full_fig_p024_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: FIG. 20. Same as in Fig [PITH_FULL_IMAGE:figures/full_fig_p025_20.png] view at source ↗
read the original abstract

Relativistic wavefunctions of nucleon excitations are scrutinised to understand their node structure and the underlying role of local interpolating fields in generating the nucleon spectrum. In addressing quark model perspectives, approximately 4000 propagators are employed on the heaviest PACS-CS ensemble at $m_\pi \simeq$ 702 MeV. We examine the ground and four lowest-lying excited states at zero momentum for both positive- and negative-parity spectra, where the proton's d-quark wavefunction is calculated about the two u quarks at the origin. This is achieved using two local interpolating fields that each carry the quantum numbers of the nucleon but with differing spin-flavour structures, one of which vanishes in the nonrelativistic limit. We find that two distinct types of wavefunction nodes are manifest: "superposition nodes" formed through a linear combination of interpolating fields, and novel "built-in nodes" that are fundamentally built in to the s-wave Dirac components of an individual interpolating field. These are investigated qualitatively through visualisations in the form of both volume and surface renderings, and quantitatively by the calculation of radial wavefunctions. Combined, these findings build a comprehensive picture of the single-particle nucleon spectrum and how its properties derive from fundamental lattice operators.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a lattice QCD study on the PACS-CS ensemble (m_π ≃ 702 MeV) that computes relativistic wavefunctions of the nucleon ground state and its four lowest-lying excitations (zero momentum, both parities) using a 2×2 variational basis of two local interpolating fields. It identifies two classes of nodes in the d-quark wavefunction about the two u quarks at the origin—'superposition nodes' arising from linear combinations of the operators and 'built-in nodes' intrinsic to the s-wave Dirac components of a single operator—supported by volume/surface visualizations and radial profiles derived from approximately 4000 propagators.

Significance. If the extracted states are sufficiently close to the true eigenstates, the distinction between the two node types offers a concrete link between the choice of lattice interpolators and the nodal structure of the nucleon spectrum, complementing quark-model intuition with first-principles QCD information. The scale of the computation and the combination of qualitative renderings with quantitative radial wavefunctions constitute a clear technical strength.

major comments (2)
  1. [Abstract/Methods] Abstract and Methods: The central claim that the observed nodes reflect intrinsic properties of the nucleon excitations presupposes that the variational procedure with only two local interpolators isolates the ground and four excited states with acceptable purity. A 2×2 correlation matrix is minimal; without additional operators, explicit checks for excited-state contamination, or overlap factors with known physical states, the apparent node locations in the excited states could be shifted by mixing rather than representing true eigenstate structure.
  2. [Results] Results: The radial wavefunction profiles and node classifications are presented without error bars or any quantification of statistical or systematic uncertainties. In addition, the analysis is performed exclusively at m_π ≃ 702 MeV; no assessment is given of how the heavy-pion spectrum ordering or wavefunction overlaps might differ from the physical point, which directly affects the interpretation of the built-in versus superposition nodes.
minor comments (2)
  1. [Abstract] The abstract refers to 'approximately 4000 propagators' but provides no further details on the number of configurations, source-sink separations, or smearing parameters used to generate the wavefunctions.
  2. Figure captions and text should explicitly state the normalization convention for the radial profiles and the precise definition of the radial coordinate used in the surface renderings.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their thorough review and constructive feedback on our manuscript. We address each major comment below, indicating planned revisions where appropriate while providing our honest assessment of the work's scope and limitations.

read point-by-point responses
  1. Referee: [Abstract/Methods] The central claim that the observed nodes reflect intrinsic properties of the nucleon excitations presupposes that the variational procedure with only two local interpolators isolates the ground and four excited states with acceptable purity. A 2×2 correlation matrix is minimal; without additional operators, explicit checks for excited-state contamination, or overlap factors with known physical states, the apparent node locations in the excited states could be shifted by mixing rather than representing true eigenstate structure.

    Authors: We agree that a 2×2 variational basis is minimal and that the purity of the extracted states is central to interpreting the nodes as intrinsic. The two local interpolators were selected for their complementary spin-flavour structures (one vanishing in the non-relativistic limit) to span the low-lying spectrum effectively. In the revised manuscript we will add effective-mass plots of the principal correlators, the eigenvalues of the 2×2 matrix, and a quantitative discussion of the overlap factors to demonstrate state isolation and address possible mixing. We maintain that the observed distinction between superposition and built-in nodes remains robust within this controlled setup, but we will qualify the claims accordingly. revision: partial

  2. Referee: [Results] The radial wavefunction profiles and node classifications are presented without error bars or any quantification of statistical or systematic uncertainties. In addition, the analysis is performed exclusively at m_π ≃ 702 MeV; no assessment is given of how the heavy-pion spectrum ordering or wavefunction overlaps might differ from the physical point, which directly affects the interpretation of the built-in versus superposition nodes.

    Authors: We accept that the radial profiles lack error bars and that the single heavy-pion ensemble limits the generality of the conclusions. In the revision we will include statistical uncertainties on the radial wavefunctions computed from the ~4000 propagators. We will also add an explicit discussion of the heavy-mass limitation, noting that the PACS-CS ensemble was chosen to maximise signal quality for these computationally demanding wavefunction measurements. A quantitative assessment at the physical point is not possible with the present data set. revision: yes

standing simulated objections not resolved
  • Quantitative assessment of spectrum ordering and node structure at the physical pion mass, which would require additional ensembles beyond the single heavy-mass ensemble used here.

Circularity Check

0 steps flagged

No circularity: node classifications extracted directly from lattice wavefunction computations

full rationale

The paper performs a direct lattice QCD calculation of nucleon wavefunctions on a fixed ensemble using a 2x2 variational basis of local interpolators. The distinction between superposition nodes (arising from linear combinations in the eigenstate extraction) and built-in nodes (intrinsic to individual interpolator Dirac structures) is obtained by explicit computation of overlaps, radial profiles, and visualizations of the resulting states. No parameters are fitted to define or predict the node types, no self-citations load-bear the classification, and the results do not reduce to the inputs by construction. The derivation chain remains self-contained as empirical output from the correlator matrix diagonalization and wavefunction sampling.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The analysis rests on standard lattice QCD assumptions that local interpolating fields overlap with physical states and that the chosen ensemble provides insight into the nucleon spectrum; no new free parameters, axioms, or invented entities are introduced beyond conventional lattice techniques.

axioms (2)
  • domain assumption Local interpolating fields with the correct quantum numbers have non-zero overlap with the physical nucleon eigenstates.
    Invoked when the authors attribute observed nodes to the structure of the nucleon spectrum generated by these fields.
  • domain assumption The heavy-pion ensemble at m_π ≃ 702 MeV yields node structures representative of the physical spectrum.
    Implicit in the choice of the PACS-CS ensemble for the study.

pith-pipeline@v0.9.0 · 5534 in / 1242 out tokens · 31365 ms · 2026-05-10T16:08:32.269482+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

51 extracted references · 26 canonical work pages

  1. [1]

    interpolator wavefunctions

    As such, these components can be appropriately combined to ob- tain an improved estimate without the need to explicitly calculate the wavefunctions forU ∗ configurations. In fu- ture work that considers nonzero momentum, computing wavefunctions for the{U ∗}ensemble will provide a gen- uine improvement. With this explained, we now present our plots of the ...

  2. [2]

    The lowest-lyingχ 1- andχ 2-dominated states both have zero superposition nodes

  3. [3]

    Separately add one node for eachχ1-dominated and χ2-dominated radial excitation

  4. [4]

    Add an extra node to thes-wave components ofχ 2- dominated positive-parity states andχ 1-dominated negative-parity states (the upper and lower compo- nents, respectively). C. MIT bag model Although which interpolating field carries the built-in node can be easily understood, we have yet to account for its presence in thes-wave components of the wavefunc- ...

  5. [5]

    E. E. Salpeter and H. A. Bethe, A Relativistic Equation for Bound-State Problems, Phys. Rev.84, 1232 (1951)

  6. [6]

    D. S. Roberts, W. Kamleh, and D. B. Leinweber, Wave function of the Roper from lattice QCD, Phys. Lett. B 725, 164 (2013), arXiv:1304.0325 [hep-lat]

  7. [7]

    D. S. Roberts, W. Kamleh, and D. B. Leinweber, Nucleon excited state wave functions from lattice QCD, Phys. Rev. D89, 074501 (2014), arXiv:1311.6626 [hep-lat]

  8. [8]

    Aokiet al.(PACS-CS), 2 + 1 flavor lattice QCD to- ward the physical point, Phys

    S. Aokiet al.(PACS-CS), 2 + 1 flavor lattice QCD to- ward the physical point, Phys. Rev. D79, 034503 (2009), arXiv:0807.1661 [hep-lat]

  9. [9]

    M. G. Beckett, B. Joo, C. M. Maynard, D. Pleiter, O. Tatebe, and T. Yoshie, Building the International Lattice Data Grid, Comput. Phys. Commun.182, 1208 (2011), arXiv:0910.1692 [hep-lat]

  10. [10]

    Renormalization Group Analysis of Lattice Theories and Improved Lattice Action. II -- four-dimensional non-abelian SU(N) gauge model

    Y. Iwasaki, Renormalization Group Analysis of Lat- tice Theories and Improved Lattice Action. II. Four- dimensional non-Abelian SU(N) gauge model, (1983), arXiv:1111.7054 [hep-lat]

  11. [11]

    Iwasaki, K

    Y. Iwasaki, K. Kanaya, T. Kaneko, and T. Yoshie, Scal- ing in SU(3) pure gauge theory with a renormalization- group-improved action, Phys. Rev. D56, 151 (1997), arXiv:hep-lat/9610023

  12. [12]

    Sheikholeslami and R

    B. Sheikholeslami and R. Wohlert, Improved continuum limit lattice action for QCD with wilson fermions, Nucl. Phys. B259, 572 (1985)

  13. [13]

    S. Aokiet al.(CP-PACS, JLQCD), Nonperturbative O(a) improvement of the Wilson quark action with the renormalization-group-improved gauge action using the Schr¨ odinger functional method, Phys. Rev. D73, 034501 27 (2006), arXiv:hep-lat/0508031

  14. [14]

    Chung, H

    Y. Chung, H. G. Dosch, M. Kremer, and D. Schall, Baryon sum rules and chiral symmetry breaking, Nucl. Phys. B197, 55 (1982)

  15. [15]

    Brommel, P

    D. Brommel, P. Crompton, C. Gattringer, L. Y. Gloz- man, C. B. Lang, S. Schaefer, and A. Schafer (Bern-Graz- Regensburg), Excited nucleons with chirally improved fermions, Phys. Rev. D69, 094513 (2004), arXiv:hep- ph/0307073

  16. [16]

    D. B. Leinweber, Nucleon properties from unconven- tional interpolating fields, Phys. Rev. D51, 6383 (1995), arXiv:nucl-th/9406001

  17. [17]

    M. S. Mahbub, A. O. Cais, W. Kamleh, D. B. Leinwe- ber, and A. G. Williams, Positive-parity excited states of the nucleon in quenched lattice QCD, Phys. Rev. D82, 094504 (2010), arXiv:1004.5455 [hep-lat]

  18. [18]

    M. S. Mahbub, W. Kamleh, D. B. Leinweber, P. J. Moran, and A. G. Williams, Structure and flow of the nucleon eigenstates in lattice QCD, Phys. Rev. D87, 094506 (2013), arXiv:1302.2987 [hep-lat]

  19. [19]

    J. M. Zanotti, D. B. Leinweber, A. G. Williams, J. B. Zhang, W. Melnitchouk, and S. Choe (CSSM Lattice), Spin-3/2 nucleon and ∆ baryons in lattice QCD, Phys. Rev. D68, 054506 (2003), arXiv:hep-lat/0304001

  20. [20]

    B. G. Lasscock, J. N. Hedditch, W. Kamleh, D. B. Lein- weber, W. Melnitchouk, A. G. Williams, and J. M. Zan- otti, Even parity excitations of the nucleon in lattice QCD, Phys. Rev. D76, 054510 (2007), arXiv:0705.0861 [hep-lat]

  21. [21]

    M. S. Mahbub, A. O. Cais, W. Kamleh, B. G. Lasscock, D. B. Leinweber, and A. G. Williams, Isolating excited states of the nucleon in lattice QCD, Phys. Rev. D80, 054507 (2009), arXiv:0905.3616 [hep-lat]

  22. [22]

    M. S. Mahbub, A. O. Cais, W. Kamleh, B. G. Lasscock, D. B. Leinweber, and A. G. Williams, Isolating the Roper resonance in lattice QCD, Phys. Lett. B679, 418 (2009), arXiv:0906.5433 [hep-lat]

  23. [23]

    Gusken, A study of smearing techniques for hadron correlation functions, Nucl

    S. Gusken, A study of smearing techniques for hadron correlation functions, Nucl. Phys. B Proc. Suppl.17, 361 (1990)

  24. [24]

    Michael, Adjoint sources in lattice gauge theory, Nucl

    C. Michael, Adjoint sources in lattice gauge theory, Nucl. Phys. B259, 58 (1985)

  25. [25]

    Luscher and U

    M. Luscher and U. Wolff, How to calculate the elastic scattering matrix in two-dimensional quantum field the- ories by numerical simulation, Nucl. Phys. B339, 222 (1990)

  26. [26]

    M. S. Mahbub, W. Kamleh, D. B. Leinweber, P. J. Moran, and A. G. Williams (CSSM Lattice), Roper reso- nance in 2+1 flavor QCD, Phys. Lett. B707, 389 (2012), arXiv:1011.5724 [hep-lat]

  27. [27]

    Draper, R

    T. Draper, R. M. Woloshyn, W. Wilcox, and K.-F. Liu, Electromagnetic form factors of hadrons, Nucl. Phys. B Proc. Suppl.9, 175 (1989)

  28. [28]

    C. D. Abell, D. B. Leinweber, Z.-W. Liu, A. W. Thomas, and J.-J. Wu, Low-lying odd-parity nucleon resonances as quark-model-like states, Phys. Rev. D108, 094519 (2023), arXiv:2306.00337 [hep-lat]

  29. [29]

    Kamleh, D

    W. Kamleh, D. Leinweber, Z.-w. Liu, F. Stokes, A. Thomas, S. Thomas, and J.-j. Wu, Structure of the Nucleon and its Excitations, EPJ Web Conf.175, 06019 (2018), arXiv:1711.05413 [hep-lat]

  30. [30]

    Melnitchouk, S

    W. Melnitchouk, S. O. Bilson-Thompson, F. D. R. Bon- net, J. N. Hedditch, F. X. Lee, D. B. Leinweber, A. G. Williams, J. M. Zanotti, and J. B. Zhang, Excited baryons in lattice QCD, Phys. Rev. D67, 114506 (2003), arXiv:hep-lat/0202022

  31. [31]

    Chodos, R

    A. Chodos, R. L. Jaffe, K. Johnson, C. B. Thorn, and V. F. Weisskopf, New extended model of hadrons, Phys. Rev. D9, 3471 (1974)

  32. [32]

    Chodos, R

    A. Chodos, R. L. Jaffe, K. Johnson, and C. B. Thorn, Baryon structure in the bag theory, Phys. Rev. D10, 2599 (1974)

  33. [33]

    C. L. Critchfield, Scalar binding of quarks, Phys. Rev. D 12, 923 (1975)

  34. [34]

    C. L. Critchfield, Scalar potentials in the Dirac equation, J. Math. Phys.17, 261 (1976)

  35. [35]

    D. W. Rein, Low-energy properties of a relativistic quark model with linear scalar potential, Nuovo Cim. A38, 19 (1977)

  36. [36]

    Ram and R

    B. Ram and R. Halasa, Dirac equation with a quadratic scalar potential—are quarks confined?, Lett. Nuovo Cim. 26, 551 (1979)

  37. [37]

    Ram, Comments on the Dirac equation with a quadratic scalar potential, Lett

    B. Ram, Comments on the Dirac equation with a quadratic scalar potential, Lett. Nuovo Cim.28, 476 (1980)

  38. [38]

    Ram and M

    B. Ram and M. Arafah, Bound-state dirac eigenvalues for scalar potentials, Lett. Nuovo Cim.30, 5 (1981)

  39. [39]

    V. P. Iyer and L. K. Sharma, Scalar potentials in the Dirac equation and quark confinement, Phys. Lett. B 102, 154 (1981)

  40. [40]

    Ravndal, A harmonic quark bag model, Phys

    F. Ravndal, A harmonic quark bag model, Phys. Lett. B 113, 57 (1982)

  41. [41]

    Su and Z.-Q

    R.-k. Su and Z.-Q. Ma, The Confinement Properties for Dirac Equation With the Scalar Like and Vector Like Potentials, J. Phys. A19, 1739 (1986)

  42. [42]

    F. M. Stokes, W. Kamleh, D. B. Leinweber, M. S. Mah- bub, B. J. Menadue, and B. J. Owen, Parity-expanded variational analysis for nonzero momentum, Phys. Rev. D92, 114506 (2015), arXiv:1302.4152 [hep-lat]

  43. [43]

    ’t Hooft, On the phase transition towards permanent quark confinement, Nucl

    G. ’t Hooft, On the phase transition towards permanent quark confinement, Nucl. Phys. B138, 1 (1978)

  44. [44]

    ’t Hooft, A property of electric and magnetic flux in non-Abelian gauge theories, Nucl

    G. ’t Hooft, A property of electric and magnetic flux in non-Abelian gauge theories, Nucl. Phys. B153, 141 (1979)

  45. [45]

    H. B. Nielsen and P. Olesen, A quantum liquid model for the QCD vacuum: Gauge and rotational invariance of domained and quantized homogeneous color fields, Nucl. Phys. B160, 380 (1979)

  46. [46]

    Greensite, The confinement problem in lattice gauge theory, Prog

    J. Greensite, The confinement problem in lattice gauge theory, Prog. Part. Nucl. Phys.51, 1 (2003), arXiv:hep- lat/0301023

  47. [47]

    Greensite, Confinement from Center Vortices: A re- view of old and new results, EPJ Web Conf.137, 01009 (2017), arXiv:1610.06221 [hep-lat]

    J. Greensite, Confinement from Center Vortices: A re- view of old and new results, EPJ Web Conf.137, 01009 (2017), arXiv:1610.06221 [hep-lat]

  48. [48]

    O’Malley, W

    E.-A. O’Malley, W. Kamleh, D. Leinweber, and P. Moran,SU(3) centre vortices underpin confinement and dynamical chiral symmetry breaking, Phys. Rev. D 86, 054503 (2012), arXiv:1112.2490 [hep-lat]

  49. [49]

    Trewartha, W

    A. Trewartha, W. Kamleh, and D. Leinweber, Evidence that centre vortices underpin dynamical chiral symmetry breaking in SU(3) gauge theory, Phys. Lett. B747, 373 (2015), arXiv:1502.06753 [hep-lat]

  50. [50]

    Trewartha, W

    A. Trewartha, W. Kamleh, and D. Leinweber, Centre vortex removal restores chiral symmetry, J. Phys. G44, 125002 (2017), arXiv:1708.06789 [hep-lat]

  51. [51]

    Kamleh, Evolving the COLA software library, PoS LATTICE2022, 339 (2023), arXiv:2302.00850 [hep-lat]

    W. Kamleh, Evolving the COLA software library, PoS LATTICE2022, 339 (2023), arXiv:2302.00850 [hep-lat]