The Quantum Chromodynamics (QCD) is based on the SU(3)-color interaction, and is in general considered as a theory, which is both mathematically complicated and numerically hard-to-handle.

We present here a new calculation method (on-lattice minimization of action), which is numerically simpler than the other methods, and which uses a new ansatz for wavefunctions.

We also present calculation results for energy-mass of first-generation hadrons, which agree well with the observed values, and new results for internal component distribution, which give interesting insights into the symmetry and internal structure of these hadrons.

In Chapter 2 we compare the different calculation methods.

In Chapter 3 we describe the ansatz for the component quark and gluon wavefunctions.

In Chapter 4 the numerical algorithm is described.

In Chapter 5 the calculation results for energy-mass and component distribution are presented for the three families of first-generation hadrons: nucleons, pseudo-scalar mesons, and vector mesons.

[1] [2] [3] [4]

Basically, there are four solution methods in lattice-QCD (LQCD):

Perturbative analytic Feynman solution

[5] [6] [7]

Here one calculates the reaction cross-sections from Feynman diagrams evaluating the corresponding Feynman-integrals in analogy to the QED. As the QCD is renormalizable, all Feynman integrals can be made finite.

However, this works only for convergent Feynman series, i.e. if the interaction constant g_c < 1. This is the case for large energies E > E_Λ = 220 MeV.

Non-perturbative on-lattice Wilson-loop method

[8] [9] [10]

Here the expectation value of an operator (e.g. energy = Hamilton operator) is calculated using path integrals

〈 O 〉 = T r ∫ Π μ d U μ ( x ) O det ( M [ U ] ) exp ( − S [ U ( x ) ] ) (1)

with interaction matrix M[U], under U μ ( x ) = exp ( i   g a   A μ ( x ) ) the local gauge transformation, with coupling constant g, lattice step size a, gluon field A μ ( x ) , action S [ U ( x ) ] , on closed loops on the lattice.

We get the Wilson action on equidistant lattice L ( x k , t k ) = t k × r k × θ k × ϕ k with lattice constant a

S ( A i , ψ i ) = S f + S g the fermionic action becomes (2)

S f 0 = m q ∑ x ψ ¯ ( x ) ψ ( x ) + 1 2 a ∑ x ψ ¯ ( x ) γ μ ( U μ ( x ) ψ ( x + μ ) − U μ † ( x − μ ) ψ ( x − μ ) )

with lattice translations y = x + μ   a U μ ( x ) ≡ U ( x , x + μ ) = e i g a A μ ( x + μ 2 ) .

And the gluon action is

S g = 2 g 2 ∑ x ∑ μ < ν R e   T r ( 1 − W μ ν 1 × 1 ( x ) ) (3)

with Wilson plaquette action

W μ ν 1 × 1 ( x ) = U μ ( x ) U ν ( x + μ ) U μ † ( x + μ ) U ν † ( x ) = e i g a ( A μ ( x + μ 2 ) + A ν ( x + μ + ν 2 ) − A μ ( x + ν + μ 2 ) − A ν ( x + ν 2 ) )

The Wilson action S ( A i , ψ i ) is minimized in A i , μ ( x , E i ) ψ i μ ( x , E i ) on lattice L ( x k , t k ) .

The solution A i , μ ( x , E i ) ψ i μ ( x , E i ) yields the corresponding masses m q , i for fermions and energies E i for gluons.

Non-perturbative on-lattice eom solution

[11]

The QCD equations-of-motion (eom) are derived from the minimal-action-principle as the Euler-Lagrange-equations corresponding to the QCD Lagrangian. They are ∂ μ F a μ ν + g f a b c A b μ F c μ ν = 0 the Yang-Mills-equations for the gluon wavefunction A a μ ( x ) and the color-field-tensor F a μ ν ( x )

F a μ ν = ∂ μ A a ν − ∂ ν A a μ + g f a b c A b μ A c ν (4)

or taking into account self-interaction with the self-current J μ   a = g ψ ¯ a γ μ ψ a the Yang-Mills field equations become

∂ μ F a μ ν + g f a b c A μ   b F c μ ν = g ψ ¯ a γ μ ψ a (5)

and the Dirac equation

( i   ℏ   D μ γ μ − m c ) ψ a = 0 (6)

or ( i   D μ γ μ − m ) ψ a = 0 with the usual convention ℏ = c = 1 with the color-covariant-derivative D μ = I   ∂ μ − i g T a A a μ and the quark-wavefunction ψ a μ ( x ) .

These are n_q + n_g partial differential equations (pdeq) first order in x μ , for the n_q = 2 or n_q = 3 quarks and n_g = 8 gluons, adding a gauge condition and a boundary condition for A a μ ( x ) .

They must be solved numerically on a lattice as an eigenvalue problem of the Dirac equation, which is very difficult and time-consuming for a one-dimensional lattice of, say, n₁ = 100 points (total number of points n = n₁⁴ = 10⁸).

Non-perturbative on-lattice minimization of action (new)

The starting point is the minimum-action-principle for QCD:

S = ∫ L Q C D ( x μ , q i , A g i ) d x = min (7)

with a gauge condition and a boundary condition for A g i ( x ) .

It can be extended to QCD + QED

S = ∫ ( L Q C D ( x μ , q i , A g i ) + L Q E D ( x μ , q i , A e i ) ) d x = min

for the quarks q i , QCD-gluons A g i , QED-photons A e i .

In order to carry out the minimization numerically, we introduce an equidistant 4-dimensional lattice

L ( t k , r k , θ k , ϕ k ) = t k × r k × θ k × φ k ,

extract a small random sub-lattice L_sub.

And approximate the integral by a sum over L_sub:

S ˜ = ∑ x ∈ L s u b L Q C D ( x , q i , A g i ) Δ V (8)

where Δ V = Δ t   Δ r   Δ θ   Δ φ is the elementary integration volume in spherical coordinates, and model the quark wavefunctions as parameterized Gauss functions q = q ( x ,   p a r ( q ) ) .

And the gluon-wavefunctions as Ritz-Galerkin series on a function system f k ( x ) with coefficients α k :

A g = A g ( ∑ α k f k ( x ) ,   p a r ( A g ) ) ,

accordingly the photon-wavefunctions A e = A e ( x , { α k } , p a r ( A e ) ) .

We impose the gauge condition for A g i : ∂ μ A g i μ = 0 and a boundary condition: A g i ( r = r 0 ) = 0 , the quark-wavefunctions are normalized ∫ q i ( x ) d 3 x = 1 .

The minimization is carried out in dependence on

p a r = ( p a r ( q ) , { α k } , p a r ( A g ) )

p a r 0 = min ( S ˜ , p a r ) ,

where par₀ yields information about the energy (= mass), the sizes and the inner structure of the considered hadron.

Gluon wavefunction

For the gluon wavefunction we apply here the full Ritz-Galerkin series on the function system

f k ( r , θ ) = { b f u n c ( r , r 0 , d r 0 ) r k 1 , k 1 = 0 , ⋯ , n r }     × { ( cos k 2 θ , cos k 2 θ sin θ ) , k 2 = 0 , ⋯ , n θ } with coefficients α k ,

where b f u n c ( r , r 0 , d r 0 ) = 1 1 + exp ( r − r 0 d r 0 ) is a Fermi-step-function which limits the region r ≤ r 0 of the hadron with “smearing width” dr₀.

A g ( t , r , θ ) = { ( A g i ( t , r , θ ) cos a A i A g i ( t , r , θ ) sin a A i ) , i = 1 , ⋯ , 8 } ,

where a A i is the phase angle between the particle and the anti-particle part of the gluon, and with the Ritz-Galerkin-expansion

A g k ( t , r , θ ) = ∑ j α [ k , j ] f j ( r , θ ) exp ( − i   t   E A k ) with energies E A k .

Because of color-symmetry, the active (non-zero) gluons are (9)

A g = { A g 1 , ⋯ , A g 8 } all gluons for nucleons.

A g = { A g 1 , A g 2 , A g 4 , A g 5 , A g 6 , A g 7 } 6 non-diagonal gluons for vector-mesons.

A g = { A g 2 , A g 5 , A g 7 } 3 quark-antiquark gluons for for pseudo-scalar mesons.

Quark wavefunction

The first-generation (u, d)-hadrons consist of three quarks (nucleons) or three color-symmetric quark-antiquark-combinations (vector-mesons) or two quark-antiquark-combinations (pseudo-scalar mesons) (10).

For nucleons

q = { ( q 1 0 )   ,   ( q 2 0 ) , ( q 3 0 ) }

For vector-mesons

q = { ( ( q 1 q ¯ 1 ) ± ( q 2 q ¯ 2 ) 2 )   ,   ( ( q 1 q ¯ 1 ) ± ( q 2 q ¯ 2 ) 2 )   ,   ( ( q 1 q ¯ 1 ) ± ( q 2 q ¯ 2 ) 2 ) }

or

q = { ( q 1 q ¯ 2 )   ,   ( q 1 q ¯ 2 )   ,   ( q 1 q ¯ 2 ) }

For pseudo-scalar mesons (pi+, pi0)

q = { ( q 1 0 )   ,   ( 0 q ¯ 2 )   ,   0 } or q = { ( ( q 1 0 ) + ( 0 q ¯ 2 ) 2 )   ,   ( ( q 1 0 ) − ( 0 q ¯ 2 ) 2 )   ,   0 }

A Ritz-Galerkin series for quarks would blow up the complexity of calculation, therefore we use here a simpler model, based on the asymptotic-freedom property of quarks: gaussian “blobs”

q k ( t , r , θ ) = exp ( − i   t   E u k ) exp ( − ( r → − r → u , k ) 2 2   d r u , k ) cos a k ,

where E u k is the energy, r → u , k = ( r u k , θ u k ) and d r u , k is the position (r, θ) and its width, a k is the quark-antiquark phase and the antiquark is

q ¯ k ( t , r ) = exp ( − i   t   E u k ) exp ( − ( r → − r → u , k ) 2 2   d r u , k ) sin a k

The ansatz and the color symmetry

The form of the quark color-wavefunction and the corresponding set of active gluons are enforced by the color-symmetry and the number of particles equal to the number of combinations.

The 8 gluons of the SU(5) form 3 families:

• The diagonal { A g 3 , A g 8 } , which map color indices into itself,

• The non-diagonal { A g 1 , A g 4 , A g 6 } , which exchange color-index with a different color index, and

• The non-diagonal { A g 2 , A g 5 , A g 7 } , which exchange color-index with a different anti-color index.

The nucleons consist of three quarks with color (r, g, b), and the color wavefunction q is mapped into itself under color-permutations, therefore the full set of 8 gluons Ag_i is required, and there are only two possibilities for first-generation hadrons: p = uud and n = ddu.

The vector mesons consist of quark-antiquark pairs, where the color wavefunction q has three identical components.

q is mapped into itself under the corresponding set of 6 non-diagonal gluons A g = { A g 1 , A g 2 , A g 4 , A g 5 , A g 6 , A g 7 } (each flips two color indices).

It is seen immediately that the three combinations listed above are the only possible ones, which is confirmed by the existence of the three v-mesons omega0, rho0, rho+.

The pseudo-scalar mesons consist of quark-antiquark pairs, where the color wavefunction q has two non-zero components. The corresponding gluon set are the 3 non-diagonal color-anticolor gluons A g = { A g 2 , A g 5 , A g 7 } , which exchange a color-quark with a different anti-color-quark. For example, A g 2 flips color-indices (3, 1) and transforms q 12 = { ( q 1 c 0 )   ,   ( 0 q ¯ 2 c ¯ )   ,   0 } into q 23 = { 0 ,   ( q 1 c ¯ 0 )   ,   ( 0 q ¯ 2 c ) } . So in reality, the wavefunction is a superposition of the three ( q 12 ,   q 23 ,   q 31 ) and is mapped by the gluon set into itself.

Again, one can see immediately that there are only two possible combinations, which correspond to the two known ps-mesons pi+ and pi0.

The corrected coupling constant

In the original Callan-Symanzik relation the QCD coupling constant has the (asymptotic) energy dependence

α s ( μ ) = g 2 ( μ ) 4 π = 1 8 π β 0 log ( μ Λ ) = 12 π ( 33 − 2 n f ) log ( μ 2 Λ 2 )

g ( μ ) = 4 π 3 ( 11 N − 2 n f ) 2 log ( μ Λ )

where

Λ ≈ 220 MeV critical QCD energy;

n_f = 3 generations, N = 3 QCD charges.

For energies μ ≈ Λ it must be modified to avoid the singularity

g c ( μ ) = 4 π 3 54 ( log ( μ Λ ) ) 2 + c G E 0 2 (11)

for the numerical calculation we set

c G E 0 = 1 log ( m ( p ) Λ Q C D ) = 0.683 ≈ log 2 ,

which is consistent with the Callan-Symanzik relation for μ > 2 Λ , as shown below (Figure 1).

The energy, length, and time are made dimensionsless by using the units: E ( E 0 = ℏ c 1   fm = 0.196   GeV ), r(fm), t(fm/c) fm = 10⁻¹⁵ m. The hadrons have axial symmetry, so we can set φ = 0 and use the spherical coordinates (t, r, θ).

We choose the equidistant lattice for the intervals ( t , r , θ ) ∈ [ 0 , 1 ] × [ 0 , 1 ] × [ 0 , π ] with 21 × 21 × 11 points and, for the minimization 8× in parallel, 8 random sublattices:

l [ i x , j ] = { { ( t i 1 , r i 2 , t i 3 ) | ( i 1 , i 2 , i 3 ) = r a n d o m ( l a t t i c e ,   j = 1 , ⋯ , 100 ) } | i x = 1 , ⋯ , 8 } .

For the Ritz-Galerkin expansion we use the 12 functions

f k ( r , θ ) = { b f u n c ( r , r 0 , d r 0 ) r k 1 , k 1 = 0 , ⋯ , n r } × { ( cos k 2 θ , cos k 2 θ sin θ ) , k 2 = 0 , ⋯ , n θ }

The action

S = ∫ L Q C D ( x μ , q i , A g i )   r 2 sin θ   d t   d r   d θ   d φ (12)

becomes a mean-value on the sublattice l [ i x ]

S ˜ [ i x ] = 1 N ( l [ i x ] ) ∑ x ∈ l [ i x ] s u b L Q C D ( x ,   q i , A g i )   2 π V t r θ (12a)

where V t r θ = π the ( t , r , θ ) -volume and N ( l [ i x ] ) is number of points of the sublattice, we set N ( l [ i x ] ) = 100 .

We impose the gauge condition and the boundary condition for A g i via penalty-function (imposing exact conditions is possible, but slows down the minimization process enormously).

S ˜ is minimized 8× in parallel with the Mathematica-minimization method “simulated annealing”, the execution time on a 2.7 GHz Xeon E5 is 9100s for the proton p = uud, the complexity K ( S ˜ [ i x ] ) = 8.4 million terms.

The minimization is performed in the parameters p a r = ( p a r ( q ) , { α k } , p a r ( A g ) ) , for the proton is the number of parameters N ( { α k } ) = 16 × 12 = 164 , N ( p a r ( q ) ) = 3 × 5 = 15 , N ( p a r ( A g ) ) = 8 × 2 = 16 .

The proper parameters of the quarks and the gluons are:

p a r ( q i ) = { E u i ,   a i ,   r u i ,   θ u i ,   d r u i } , p a r ( A g i ) = { E A i ,   a A i }

Criteria for correctness of the ansatz

1) Convergence of minimization

As we found out during the computation, a wrong ansatz, e.g. lacking color symmetry, leads to a non-convergent minimization. We chose a high goal precision of prec = 10⁻⁴, so there was a high probability that a convergent minimizations hits a real (global) minimum.

2) High relative deviation between solutions

Strongly differing solutions indicate a non-correct ansatz, as we found out e.g. for the nucleons with too many degrees-of-freedom for the gluons: the relative deviation for crucial variables, like energy, should be no more than 2% for the nucleons and 6% for the ps-mesons.

3) Vanishing parameter-derivatives

A true minimum must satisfy the derivative-condition ∂ S ∂ p i = 0 , where p i is one of the minimization parameters, Normally, the parameter-derivatives are close to zero, otherwise the minimum is not genuine, or the ansatz is wrong.

4) Boundary condition and gauge condition

The boundary and gauge condition must have values close to zero, otherwise the weight for the penalty function is too low.

5) Minimum value

The minimum value should be −30, …, 30 for the considered parameter range. Very large positive values result in the case of too high penalty weights. Very large negative values may come out, if the Ritz-Galerkin parameters α i are not bounded appropriately.

6) Correct energy scale and number of particles

The three types of first-generation hadrons have energy scales: E (nucleon) ≈ 0.98 GeV, E (v-meson) ≈ 0.78 GeV, E (ps-meson) ≈ 0.14 GeV, and these values emerge automatically with 8, 6 and 3 gluons respectively.

Furthermore, with the above ansatz, the number of possible particles is 2, 3, 2 respectively.

[13] [14]

nucleons n, p quarks (3), gluons (8), spin = 1/2.

Masses (Table 1).

Energy for quark-number (n = 1, 2, 3), gluon-number (n = 4, …, 11), both sorted with increasing energy (Figure 2).

Distribution quarks (r[fm], θ) sorted with increasing energy (Figure 3).

The quark distribution differs largely between the nucleons: the proton is ring-symmetric (no θ-component), the neutron has two orbitals with an angle of α = π/2. The small mass difference is due to the electromagnetic contribution, which is about 1% of the total mass.

The mass of the nucleons, as is the case for all first-generation-hadrons is generated almost exclusively by the energy of the gluons and the quarks, the rest masses of u and d (m_u = 2.3 MeV, m_d = 4.8 MeV) contribute very little to the total mass.

The gluon distribution is practically the same for both nucleons, which is to be expected, since the two particles are identical for the color interaction.

The radius of the nucleons can be assessed from the above diagram: r(p) ≈ 0.8 fm, r(n) ≈ 1 fm.

Proton p = uud.

m = 0.938 GeV, r₀ = 0.84 fm.

E_tot = 0.945 GeV, ΔE_tot = 0.032, dE_em = −0.013 (Table 2).

Mean calculation error ΔEu_i ΔEA_i Δa_i ΔaA_i Δdru_i Δru_i Δθu_i [12] .

Gluons Agi (Figure 4).

The proton p has one rotation plane (orbital), the two quarks (u, d) are close at r = 0.15 low energy E ≤ 0.03, the second u-quark further outside r = 0.4, and high energy E = 0.2. The “smearing” width is comparable, δr ≈ 0.3.

The electromagnetic correction is negative and much larger than with the neutron, dE_em = −0.013 GeV, which is apparently the reason for the proton’s smaller mass.

Neutron n = ddu

m = 0.939 GeV, r₀ = 0.84 fm.

E_tot = 0.945 GeV, ΔE_tot = 0.018, dE_em = +0.0017 (Table 3).

Mean calculation error ΔEu_i ΔEA_i Δa_i ΔaA_i Δdru_i Δru_i Δθu_i [12] .

Gluons Agi (Figure 5)

The neutron n has two orbitals with an angle of α = π/4, the u-quark is at the center with low energy E = 0.05, the two d-quarks sit in the orbitals with higher energies E = 0.09, 0.013. The “smearing” width is comparable, δr ≈ 0.4 and higher than with the proton.

The gluon distribution is practically the same as for the proton, which is to be expected, since the two particles are identical for the color interaction.

The electromagnetic correction is positive and much smaller than with the proton, dE_em = +0.0017 GeV, which is apparently the reason for the proton’s smaller mass.

Magnetic moment of nucleons

The magnetic moment is μ = q 2 m L = q 2 m m ω   r 2 , for a rotating charge distribution:

μ = ω 2 ∑ i q i r i 2 = ω 2 I q ,

where I q = ∑ i q i r i 2 → ∫ r 2 d q is the momentum of charge, in analogy to the momentum of inertia I m = ∫ r 2 d m .

For a rotating solid sphere with radius r₀ with constant charge density I q = 2 5 q r 0 2 .

The magnetic moment of the nucleons is measured in nuclear magnetons μ N = e ℏ 2 m , which is the magnetic moment of a rotating solid sphere with constant charge density μ N = ω 2 I q ( s p h e r e ) = ω 2 2 5 e r 0 2 .

The actual momentum of charge is therefore:

I q = ∑ i q ( q i ) r ( q i ) 2

We have to take into account the “smearing” Δr_i of radius r_i

〈 r i 2 〉 = ∫ r 2 exp ( − ( r − r i ) 2 2 Δ r ) d r ∫ exp ( − ( r − r i ) 2 2 Δ r ) d r ,

so it becomes

I q = ∑ i q ( q i ) 〈 r ( q i ) 2 〉

We get for the neutron.

I_qn = −0.1766e, I_qNn = 0.106e, so I q n I q N n = −1.766, measured μ μ N = −1.91.

And for the proton

I_qp = + 0.2226e, I_qNp = 0.0909e, so I q p I q N p = + 2.448, measured μ μ N = +2.793.

The calculation does not take into account the orbitals, and there is also the statistical uncertainty of the order 7%, so the results are satisfactory.

Pseudo-scalar mesons pi+, pi0 quarks (2), gluons (3), spin = 0-

[6]

Masses (Table 4)

Energy for quark-number (n = 1, 2), gluon-number (n = 3, 4, 5), both sorted with increasing energy (Figure 6)

Distribution quarks (r[fm], θ): independent (θ) = spherical (Figure 7)

The pseudo-scalar mesons are spherically-symmetric, there is no θ-dependence: θ ≈ 0 in the quark-distribution, the gluon-wavefunctions show little θ-dependence, and the gluon amplitudes are much smaller (factor 30) for pi0 than for pi+.

For the pi0, u u ¯ and d d ¯ sit at r = 0.4 E ≈ 0, and at r = 0.75 E ≈ 0.1.

For the pi+, the u and d ¯ have practically equal radii, but different energies: r = 0.6 E ≈ 0.001, and r = 0.6 E ≈ 0.01.

The measured masses of the ps-mesons (0.135, 0.139) are reproduced by the calculation (0.155 ± 0.025, 0.129 ± 0.026), but only roughly within the error bounds.

Ps-meson pi0 = ( u u ¯ − d d ¯ ) / 2 .

m = 0.135 GeV, r₀ = 0.66 fm.

E_tot = 0.155 GeV, ΔE_tot = 0.025, dE_em = +0.007 (Table 5).

Mean calculation error ΔEu_i ΔEA_i Δa_i ΔaA_i Δdru_i Δru_i Δθu_i [12] .

Gluons Agi (Figure 8).

Ps-meson pi+ = u d ¯ .

m = 0.139 GeV, r₀ = 0.66 fm.

E_tot = 0.129 GeV, ΔE_tot = 0.026, dE_em = +0.0014 (Table 6).

Mean calculation error ΔEu_i ΔEA_i Δa_i ΔaA_i Δdru_i Δru_i Δθu_i [12] .

Gluons Agi (Figure 9).

Vector mesons rho0, rho+, omega0 quarks (2), gluons (6 non-diagonal), spin = 1.

[15]

Masses (Table 7).

Energy for quark-number (n = 1, 2), gluon-number (n = 3, …, 8), both sorted with increasing energy (Figure 10).

Distribution quarks (r[fm], θ) (Figure 11).

The vector mesons are spin-1 bosons but only rho+ shows an explicit θ-dependence of quark-distribution: it is ellipsoidal. The gluons show explicit θ-dependence and are, as for the nucleons, practically equal for all three particles.

For rho0: the quarks u u ¯ and d d ¯ have identical parameters r = 0.5, δr = 0.3, E = 0.1.

For omega0: the quarks u u ¯ and d d ¯ again have identical parameters, are at center, δr = 0.25, E = 0.1.

For rho+: the heavier quark d ¯ has r = 0.5, δr = 0.05, E = 0.05, the light quark u has r = 0.9, δr = 0.5, E = 0.07, rho+ has two orthogonal orbitals. Its two quarks have completely different width; the d ¯ quark closer to the center has a small bandwidth, the light u quark is strongly “smeared” like all the other quarks in the 3 particles.

The measured masses of the v-mesons (0.775, 0.775, 0.782) are reproduced correctly by the calculation (0.771 ± 0.0052, 0.779 ± 0.012, 0.782 ± 0.007).

V-meson rho0 = ( u u ¯ − d d ¯ ) / 2 .

m = 0.775 GeV, r₀ = 0.75 fm.

E_tot = 0.771 GeV, ΔE_tot = 0.0052, dE_em = +0.002 (Table 8).

Mean calculation error ΔEu_i ΔEA_i Δa_i ΔaA_i Δdru_i Δru_i Δθu_i [12] .

Gluons Agi (Figure 12).

V-meson rho+ = u d ¯ .

m = 0.775 GeV, r₀ = 0.75 fm.

E_tot = 0.779 GeV, ΔE_tot = 0.012, dE_em = +0.002 (Table 9).

Mean calculation error ΔEu_i ΔEA_i Δa_i ΔaA_i Δdru_i Δru_i Δθu_i [12] .

Gluons Agi (Figure 13).

V-meson omega0 = ( u u ¯ + d d ¯ ) / 2 .

m = 0.782 GeV, r₀ = 0.75 fm.

E_tot = 0.782 GeV, ΔE_tot = 0.007, dE_em = +0.002 (Table 10).

Mean calculation error ΔEu_i ΔEA_i Δa_i ΔaA_i Δdru_i Δru_i Δθu_i [12] .

Gluons Agi (Figure 14).

We present here (Chapter 2) a new calculation method for on-lattice QCD, namely on-lattice minimization of action.

It works by direct minimization on parameters of action

S = ∫ L Q C D ( x μ , q i , A g j ) d x = min ,

where the action integral is approximated as a summation on a random sublattice of an equidistant lattice.

This has the following advantages:

• The minimal-action principle is a fundamental principle, from which the equation-of-motion (eom), i.e. the Dirac equation for QCD is derived. The parameters q i , A g j of quarks q i and gluons A g j in the Lagrangian L Q C D ( x μ , q i , A g i ) (e.g. energy-mass of a quark) can be calculated from it in principle exactly without solving the eom. In comparison, the Wilson loop method uses an approximation in order to make the path integral tractable numerically.

• The calculation is a simple summation, which is very fast numerically, as opposed to the analytical integral calculation of the perturbative analytic Feynman solution.

• The calculation is scalable, i.e. the precision can be increased arbitrarily, simply by making the step size of the lattice smaller, or the size of the sublattice larger.

• The calculation can be carried-out in parallel by N_p processes on N_p different sublattices with the same number of points (in this implementation we have N_p = 8). The mean of the N_p resulting values is then the calculation result for a parameter, whereas the standard-deviation assesses the calculation error.

• On-lattice minimization of action uses parameter minimization algorithms instead of solving partial differential equations (as in on-lattice eom solution) or instead of calculation of parameterized integrals (as in Wilson-loop method or in analytic Feynman method). Nowadays, there exists a large selection of powerful algorithms for parameter minimization which can be used for this purpose.

• On-lattice minimization of action, as opposed to the other solution methods, yields information about radial and axial distribution within hadrons.

The calculation ansatz in Chapter 3 for the quark wavefunctions:

For nucleons

q = { ( q 1 0 )   ,   ( q 2 0 ) , ( q 3 0 ) }

For vector-mesons

q = { ( ( q 1 q ¯ 1 ) ± ( q 2 q ¯ 2 ) 2 )   ,   ( ( q 1 q ¯ 1 ) ± ( q 2 q ¯ 2 ) 2 )   ,   ( ( q 1 q ¯ 1 ) ± ( q 2 q ¯ 2 ) 2 ) }

or

q = { ( q 1 q ¯ 2 )   ,   ( q 1 q ¯ 2 )   ,   ( q 1 q ¯ 2 ) }

For pseudo-scalar mesons (pi+, pi0)

q = { ( q 1 0 )   ,   ( 0 q ¯ 2 )   ,   0 } or q = { ( ( q 1 0 ) + ( 0 q ¯ 2 ) 2 )   ,   ( ( q 1 0 ) − ( 0 q ¯ 2 ) 2 )   ,   0 }

And gluon wavefunctions

A g = { A g 1 , ⋯ , A g 8 } all 8 gluons for nucleons.

A g = { A g 1 , A g 2 , A g 4 , A g 5 , A g 6 , A g 7 } 6 non-diagonal gluons for vector-mesons.

A g = { A g 2 , A g 5 , A g 7 } 3 quark-antiquark gluons for pseudo-scalar mesons.

Explains effectively the mass scale for the three types of first-generation hadrons.

Nucleons M n u c ≈ 940   MeV .

Vector-mesons M v m ≈ 780   MeV .

Pseudo-scalar mesons M p s m ≈ 150   MeV .

Simply by using gluon configurations compatible with the SU(3) symmetry, i.e. corresponding to subgroups of the SU(3).

In Chapter 5 we present the results of calculations for first-generation hadrons.

The calculated masses agree well with the observed values, and the calculation shows that individual gluons on average contribute as much as individual quarks to the total energy-mass of hadrons.

On-lattice minimization of action, as opposed to the other solution methods, yields information about radial and axial distribution within hadrons, which gives interesting insights into their symmetry, “smearing-out” of components, and internal energy-mass distribution.

• Structure of nucleons

The proton is spherically symmetric, the neutron is axial with two orbitals.

The gluon distribution is practically the same for both.

The small mass difference is due to the electromagnetic contribution, which is about 1% of the total mass.

• Structure of pseudo-scalar mesons

The pseudo-scalar mesons are spherically-symmetric.

The gluon amplitudes are much smaller (factor 30) for pi0 than for pi+.

• Structure of vector mesons

The neutral v-mesons are spherically-symmetric.

For neutral rho0: the quarks u u ¯ and d d ¯ have identical parameters with distribution peak at r = 0.5.

For neutral omega0: the quarks u u ¯ and d d ¯ again have identical parameters, with distribution peak at the center.

The charged rho+ is axial, has two orthogonal orbitals, its two quarks have completely different width; the d ¯ quark closer to the center has a small bandwidth, the light u quark is strongly “smeared-out”.

The author declares no conflicts of interest regarding the publication of this paper.

Helm, J. (2024) Quantum Chromodynamics on Lattice: Direct Minimization of QCD-QED-Action with New Results. Journal of High Energy Physics, Gravitation and Cosmology, 10, 228-256. https://doi.org/10.4236/jhepgc.2024.101019