PhysRevD 108 034512
PhysRevD 108 034512
PhysRevD 108 034512
(RQCD Collaboration)
1
Institut für Theoretische Physik, Universität Regensburg, 93040 Regensburg, Germany
2
European Spallation Source ERIC, Data Management and Software Centre (DMSC),
Box 176, 22100 Lund, Sweden
We determine the axial, scalar and tensor isovector charges of the nucleon, sigma and cascade baryons as
well as the difference between the up and down quark masses, mu − md . We employ gauge ensembles with
N f ¼ 2 þ 1 nonperturbatively improved Wilson fermions at six values of the lattice spacing in the range
a ≈ ð0.039–0.098Þ fm, generated by the coordinated lattice simulations (CLS) effort. The pion mass Mπ
ranges from around 430 MeV down to a near physical value of 130 MeV and the linear spatial lattice extent
L varies from 6.5M −1 −1
π to 3.0M π , where LM π ≥ 4 for the majority of the ensembles. This allows us to
perform a controlled interpolation/extrapolation of the charges to the physical mass point in the infinite
volume and continuum limit. Investigating SU(3) flavor symmetry, we find moderate symmetry breaking
effects for the axial charges at the physical quark mass point, while no significant effects are found for the
other charges within current uncertainties.
DOI: 10.1103/PhysRevD.108.034512
034512-2
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
Note that we do not include the Λ baryon here since in utilized to ensure leading Oða2 Þ discretization effects
this case the isovector combination trivially gives zero. are then given. Finally, the strategy for interpolation/
We consider vector (V), axial-vector (A), scalar (S) and extrapolation to the physical point in the infinite volume
tensor (T) operators which are defined through the Dirac and continuum limit is outlined.
matrices ΓJ ¼ γ 4 ; γ i γ 5 ; 1; σ ij for J ∈ fV; A; S; Tg, with
σ μν ¼ 12 ½γ μ ; γ ν , where i; j ∈ f1; 2; 3g and i < j. A. Gauge ensembles
The axial charges in the ms ¼ ml ¼ 0 chiral limit are We employ ensembles generated with N f ¼ 2 þ 1 fla-
important parameters in SU(3) ChPT and enter the expan- vors of nonperturbatively OðaÞ improved Wilson fermions
sion of every baryonic quantity. These couplings can be and a tree-level Symanzik improved gauge action, which
decomposed into two LECs F and D which appear in the were mostly produced within the coordinated lattice sim-
first order meson-baryon Lagrangian for three light quark ulations (CLS) [32] effort. Either periodic or open boun-
flavors (see, e.g., Ref. [27]): dary conditions in time [33] are imposed, where the latter
choice is necessary for ensembles with a < 0.06 fm in
gNA ¼ F þ D; gΣA ¼ 2F; gΞA ¼ F − D: ð7Þ order to avoid freezing of the topological charge and thus to
ensure ergodicity [34]. The hybrid Monte Carlo (HMC)
Due to group theoretical constraints, see, e.g., simulations are stabilized by introducing a twisted mass
Refs. [28,29], such a decomposition also holds for term for the light quarks [35], whereas the strange quark is
ms ¼ ml > 0, for the axial as well as for the other charges. included via the rational hybrid Monte Carlo (RHMC)
We define for m ¼ ms ¼ ml algorithm [36]. The modifications of the target action are
corrected for by applying the appropriate reweighting, see
gNJ ðmÞ ¼ FJ ðmÞ þ DJ ðmÞ; ð8Þ Refs. [17,32,37] for further details.
In total 47 ensembles were analyzed spanning six lattice
gΣJ ðmÞ ¼ 2FJ ðmÞ; ð9Þ spacings a in the range 0.039 fm ≲ a ≲ 0.098 fm, with
pion masses between 430 MeV and 130 MeV (below the
gΞJ ðmÞ ¼ FJ ðmÞ − DJ ðmÞ; ð10Þ physical pion mass), as shown in Fig. 2. The lattice spatial
extent L is kept sufficiently large, where LMπ ≥ 4 for the
where F ¼ FA ð0Þ, D ¼ DA ð0Þ. The vector Ward identity majority of the ensembles. A limited number of smaller
(conserved vector current, CVC relation) implies that gNV ¼ volumes are employed to enable finite volume effects to be
gΞV ¼ FV ¼ 1 and gΣV ¼ 2FV , i.e., in this case the above investigated, with the spatial extent varying across all the
relations also hold for ms ≠ ml, with FV ðmÞ ¼ 1 ensembles in the range 3.0 ≤ LMπ ≤ 6.5. Further details
and DV ðmÞ ¼ 0. are given in Table I. The ensembles lie along three
In this article we determine the charges at many different trajectories in the quark mass plane, as displayed in Fig. 3:
positions in the quark mass plane and investigate SU(3) (i) The symmetric line: the light and strange quark
flavor symmetry breaking, i.e., the extent of violation of masses are degenerate (ml ¼ ms ) and SU(3) flavor
Eqs. (8)–(10). Due to this, other than for J ≠ V where symmetry is exact.
DV =FV ¼ 0, the functions DJ ðmÞ and FJ ðmÞ are not (ii) The tr M ¼ const line: starting at the ml ¼ ms flavor
uniquely determined at the physical point, where symmetric point, the trajectory approaches the
ms ≫ ml . At this quark mass point we will find the physical point holding the trace of the quark mass
approximate ratios DA =FA ≈ ð1.55 − 1.95Þ, DS =FS ≈ matrix (2ml þ ms , i.e., the flavor averaged quark
−ð0.2 − 0.5Þ and DT =FT ≈ 1.5. The first ratio can be mass) constant such that 2M2K þ M2π is close to its
compared to the SU(6) quark model expectation physical value.
DA ðmÞ=FA ðmÞ ¼ 3=2 (see, e.g., Ref. [30]), which is (iii) The ms ¼ const line: the renormalized strange quark
consistent with the large-N c limit [31]. mass is kept near to its physical value [38].
The latter two trajectories intersect close to the physical
point, whereas the symmetric line approaches the SU(3)
III. LATTICE SETUP
chiral limit. In figures where data from different lines are
In this section we present the details of our lattice setup. shown, we will distinguish these, employing the symbol
First, we describe the gauge ensembles employed and the shapes of Fig. 2 (triangle, circle, diamond). The excellent
construction of the correlation functions. The computation coverage of the quark mass plane enables the interpolation/
of the three-point correlation functions via a stochastic extrapolation of the results for the charges to the physical
approach is briefly discussed. Following this, we present point to be tightly constrained. In addition, considering the
the fitting analysis and treatment of excited state contri- wide range of lattice spacings and spatial volumes and the
butions employed to extract the ground state baryon matrix high statistics available for most ensembles, all sources of
elements. The renormalization factors used to match to systematic uncertainty associated with simulating at
the continuum MS scheme and the improvement factors unphysical quark mass, finite lattice spacing and finite
034512-3
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
FIG. 2. Parameter landscape of the ensembles listed in Table I. The ensembles are grouped according to the three quark mass
trajectories (see the text and Fig. 3): (left) the symmetric line (ml ¼ ms ), (middle) the tr M ¼ const line and (right) the ms ¼ const line.
TABLE I. List of the gauge ensembles analyzed in this work. The rqcd xyz ensembles were generated by the RQCD group using the
BQCD code [39], whereas all other ensembles were created within the CLS effort. Note that for H102 there are two replicas. These have
the same parameters but were generated with slightly different algorithmic setups and, therefore, have to be analyzed separately. N401
and N451 differ in terms of the boundary conditions (bc) imposed in the time direction (open (o) and antiperiodic (p), respectively). The
lattice spacings a are determined in Ref. [17]. In the second to last column t denotes the source-sink separation of the connected three-
point functions. The subscript (superscript) given for each separation indicates the number of measurements performed using the
sequential source (stochastic) method on each configuration (see Secs. III B and III C). The last column gives N cnfg , the number of
configurations analyzed.
034512-4
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
TABLE I. (Continued)
volume can be investigated. Our strategy for performing a N α ¼ ϵijk uiα ðujT Cγ 5 dk Þ; ð13Þ
simultaneous continuum, quark mass and infinite volume
extrapolation is given in Sec. III F.
Σα ¼ ϵijk uiα ðsjT Cγ 5 uk Þ; ð14Þ
B. Correlation functions
The baryon octet charges are extracted from two- and Ξα ¼ ϵijk siα ðsjT Cγ 5 uk Þ; ð15Þ
three-point correlation functions of the form
with spin index α, color indices i, j, k and C being the
X
CB2pt ðtÞ ¼ P αβ hBα ðx0 ; tÞB̄β ð0; 0Þi; ð11Þ charge conjugation matrix. The Λ baryon is not considered
þ
x0 here since three-point functions with the currents OJ ¼
X ūΓJ u − d̄ΓJ d vanish in this case. Without loss of generality,
CB3pt ðt; τÞ ¼ P αβ hBα ðx0 ; tÞOJ ðy; τÞB̄β ð0; 0Þi: ð12Þ we place the source space-time position at the origin (0, 0)
x0 ;y and the sink at ðx0 ; tÞ such that the source-sink separation in
time equals t. The current is inserted at ðy; τÞ with
Spin-1=2 baryon states are created (annihilated) using 0 ≤ τ ≤ t.1 The annihilation interpolators are projected
suitable interpolators B̄ (B): for the nucleon, Σ and Ξ, onto zero-momentum via the sums over the spatial sink
we employ interpolators corresponding to the proton, Σþ
and Ξ0 , respectively, Note that in practice we only analyze data with 2a ≤ τ ≤ t−2a.
1
034512-5
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
FIG. 3. Position of the ensembles in the quark mass plane. The tr M ¼ const line (indicated as a blue line) and ms ¼ const line (green
line) intersect close to the physical point (black cross). The symmetric line (orange line), which approaches the SU(3) chiral limit,
crosses the tr M ¼ const line around M π ¼ 411 MeV.
position, while momentum conservation (and the sum over combination of the current. The two-point functions are
y for the current) means the source is also at rest. constructed in the standard way using point-to-all propa-
We ensure positive parity via the projection operator gators. For the three-point functions either a stochastic
P þ ¼ 12 ð1 þ γ 4 Þ. For the three-point functions, P ¼ P þ for approach (described in the next subsection) or the sequen-
J ¼ V, S and P ¼ iγ i γ 5 P þ for J ¼ A, T. The latter tial source method [41] (on some ensembles in combination
corresponds to taking the difference of the polarizations with the coherent sink technique [42]) is employed. The
(in the i direction). The interpolators are constructed from stochastic approach provides a computationally cost effi-
spatially extended quark fields in order to increase the cient way of evaluating the three-point functions for the
overlap with the ground state of interest and minimize whole of the baryon octet, however, additional noise is
contributions to the correlation functions from excited introduced. The relevant measurements for the nucleon
states. Wuppertal smearing is employed [20], together with (which has the worst signal-to-noise ratio of the octet) have
APE-smeared [40] gauge transporters. The number of already been performed with the sequential source method
smearing iterations is varied with the aim of ensuring that as part of other projects by our group, see, e.g., Ref. [43].
ground state dominance sets in for moderate time separa- We use these data in our analysis and the stochastic
tions. The root mean squared light quark smearing approach for the correlation functions of the Σ and the Ξ
radii range from about 0.6 fm (for M π ≈ 430 MeV) up baryons. Note that along the symmetric line (ml ¼ ms ) the
to about 0.8 fm (for Mπ ≈ 130 MeV), whereas the hyperon three-point functions can be obtained as linear
strange quark radii range from about 0.5 fm (for the combinations of the contractions carried out for the currents
physical strange quark mass) up to about 0.7 fm (for ūΓJ u and d̄ΓJ d within the proton. Therefore, no stochastic
MK ¼ M π ≈ 270 MeV). See Sec. E.1 (and in particular three-point functions are generated in these cases.
Table 15) of Ref. [17] for further details. Figure 4 dem- We typically realize four source-sink separations with
onstrates that, when keeping for similar pion and kaon t=fm ≈ f0.7; 0.8; 1.0; 1.2g in order to investigate excited
masses the smearing radii fixed in physical units, the state contamination and reliably extract the ground state
ground state dominates the two-point functions at similar baryon octet charges. Details of our fitting analysis are
physical times across different lattice spacings. presented in Sec. III D. Multiple measurements are per-
Performing the Wick contractions for the two- and three- formed per configuration, in particular for the larger
point correlation functions leads to the connected quark- source-sink separations to improve the signal, see
line diagrams displayed in Fig. 5. Note that there are no Table I. The source positions are chosen randomly on each
disconnected quark-line diagrams for the three-point func- configuration in order to reduce autocorrelations. On
tions as these cancel when forming the isovector flavor ensembles with open boundary conditions in time only
034512-6
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
the spatial positions are varied and the source and sink time
slices are restricted to the bulk of the lattice (sufficiently
away from the boundaries), where translational symmetry
is effectively restored.
034512-7
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
Furthermore, no additional inversions are needed for f. Note that su;i ¼ sd;i since our light quarks are mass-
different sink momenta. degenerate.
As depicted in Fig. 7, we simultaneously compute the Using γ 5 -Hermiticity (γ 5 Df γ 5 ¼ D†f ) and the properties
three-point functions for a baryon propagating (forwards) given in Eqs. (17) and (18), the time slice-to-all propagator
from source time slice x4 to sink time slice x0;fwd 4 and connecting all points of the (forward and backward) sink
0;bwd
propagating (backwards) from x4 to x4 . We start with the time slices x04 to all points of any insertion time slice y4 can
definition of the stochastic source and solution vectors be estimated as
which can be used to construct the time slice-to-all
1 X
propagator (shown as a green wiggly line in Fig. 6). In sto N
0 0
the following i ∈ f1; …; N sto g is the “stochastic index,” we Gf3 ðx 0
; yÞγc0μd ≈ ½γ η ðx0 Þγc0 ½γ 5 sf3 ðyÞμd : ð20Þ
N sto i¼1 5 i
denote spin indices with Greek letters, color indices with
Latin letters (other than f or i) and we use flavor indices
f n ∈ fu; d; sg. We introduce (time partitioned) complex Z2 Combining this time slice-to-all propagator with point-
noise vectors [51,52] to-all propagators for the source position x4 , the baryonic
three-point correlation functions Eq. (12) can be factorized
pffiffiffi 0;fwdjbwd
ðZ2 ⊗ iZ2 Þ= 2 if x4 ¼ x4 ; as visualized in Fig. 6 into a spectator part (S) and an
α
ηi ðxÞa ¼ ð16Þ insertion part (I), leaving all flavor and spin indices open:
0 otherwise;
0 0 0
where the noise vector has support on time slices x0;fwd
4 and C3pt ðp0 ; qjx04 ; y4 ; x4 Þfα1αβ βγ μνγ
f2 f3 f4
0;bwd
x4 . The noise vectors have the properties 1 X sto X
3
N
0 0 0
≈ ðS ðp0 ; x04 ; x4 Þαicαβ βγ
N sto i¼1 c¼1 f1 f2
1 X
N sto
α 1
η ðxÞ ¼ O pffiffiffiffiffiffiffiffi ; ð17Þ × I f3 f4 ðq; y4 ; x4 Þμνγ
N sto i¼1 i a N sto ic Þ: ð21Þ
where we sum over repeated indices (other than f) and I f3 f4 ðq; y4 ; x4 Þμνγ
ic
X
Df ðx; yÞαβ
ab is the Wilson-Dirac operator for the quark flavor ≔ ½γ 5 sf3 ;i ðyÞμd Gf4 ðy; xÞνγ
þiq·ðy−xÞ
ð23Þ
dc e :
y
034512-8
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
FIG. 8. Left: polarized nucleon three-point correlation function propagating in the forward and backward directions for a ūγ y γ z u
current insertion obtained from one stochastic measurement on ensemble N200 (a ¼ 0.064 fm). Note that two different source positions
were needed to obtain the same source-sink separation t ¼ 14a in the two directions, see Fig. 7. Right: forward-backward average of the
stochastic three-point function shown on the left compared to that obtained using the sequential source method (with one source
position, propagating in the forward direction).
034512-9
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
three-point function (or the ratio RBJ ) to extract the charge J ∈ fA; S; Tg. As the axial, scalar and tensor chan-
gBJ .2 However, the three-quark baryon interpolators we use nels are the main focus of this study, we only
by design have only a small overlap with the multi-particle consider excluding the vector channel data.
states containing five or more quarks and antiquarks and it (ii) The parametrization: either one (“ES ¼ 1”) or two
is difficult to extract the lower lying excited state spectrum (“ES ¼ 2”) excited states are included in the fits. In
from the two-point function. Nonetheless, multiparticle the latter case, in order to stabilize the fit, we use a
states can significantly contribute to the three-point func- prior for ΔE1 corresponding to the energy gap for
tion if the transition matrix elements hnjOJ j0i are large. the lowest lying multi-particle state. As a cross-
Furthermore, depending on the current, different matrix check we repeat these fits using the average re-
elements, and hence excited state contributions, will sult obtained for ΔE2 in fits 5–8 as a prior and
dominate. In particular, one would expect the axial and leaving ΔE1 as a free parameter (fits 13–20). The
scalar currents to couple to the Bπ P-wave and Bππ S-wave widths of the priors are set to E1 =100 and to E2 =100,
states, respectively, while the tensor and vector currents respectively. In general, the contributions from
may enhance transitions between B and Bππ states when excited state to excited state transitions could not
ππ is in a P-wave. be resolved and the parameters bJ2;4 are set to zero.
The summation method [41] is an alternative approach, We also found that the tensor and vector currents
which involvesP summing the ratio over the operator insertion couple more strongly to the second excited state,
time SBJ ðtÞ ¼ t−τ τ¼τ0 RJ ðt; τÞ, where one can show that the
0 B consistent with the expectations mentioned above,
leading excited state contributions to SBJ ðtÞ only depend on t and the first excited state contributions are omitted
[rather than also on t − τ and τ as for RBJ ðt; τÞ]. However, one for these channels in the ES ¼ 2 fits. Furthermore,
needs a large number of source-sink separations (more than due to the large statistical error of the stochastic
the four values of t that are realized in this study) in order to three-point functions for the Σ and Ξ baryons in the
extract reliable results from this approach. vector channel (see Fig. 9 and the discussion in
These considerations motivate us to extract the charges Appendix A), we are not able to resolve bV1 (and
by fitting to the ratio of correlation functions using a fit analogously bV3 ). For these baryons we also set
form which takes into account contributions from up to two bV1;3 ¼ 0 in all the fits.
excited states, (iii) The fit range: two fit intervals τ ∈ ½δtj ; t − δtj are
used with δt1 ¼ n1 a ≈ 0.15 fm and δt2 ¼ n2 a ≈
RBJ ðt; τÞ ¼ bJ0 þ bJ1 ðe−ΔE1 ðt−τÞ þ e−ΔE1 τ Þ þ bJ2 e−ΔE1 t 0.25 fm.3
A typical fit to the ratios for the cascade baryon is
þ bJ3 ðe−ΔE2 ðt−τÞ þ e−ΔE2 τ Þ þ bJ4 e−ΔE2 t ; ð27Þ
shown in Fig. 9 for ensemble N302 (M π ¼ 348 MeV
and a ¼ 0.049 fm). The variation in the ground state matrix
where ΔEn ¼ EBn − EB0 denotes the energy gap between the
elements extracted from the 20 different fits is shown in
ground state and the nth excited state of baryon B and we Fig. 10, also for the nucleon on the same ensemble. See
have not included transitions between the first and the Appendix C for the analogous plot for the Σ baryon.
second excited state. The amplitude bJ0 ¼ gB;latt J gives the Overall, the results are consistent within errors, however,
charge, while bJ1;3 and bJ2;4 are related to the ground state to some trends in the results can be seen across the different
excited state and excited state to excited state transition ensembles. In the axial channel, in particular the results for
matrix elements, respectively. In practice, even when the fits involving a single excited state (fits 1–4), tend to be
simultaneously fitting to all available source-sink separa- lower than those involving two excited states (fits 5–20).
tions, it is difficult to determine the energy gaps (and The former are, in general, statistically more precise than
amplitudes) for a particular channel J. Similar to the the latter due to the smaller number of parameters in the fit.
strategy pursued in Ref. [53], we simultaneously fit to In order to study the systematics arising from any residual
all four channels J ∈ fV; A; S; Tg for a given baryon. As the excited state contamination in the final results at the physical
same energy gaps are present, the overall number of fit point (in the continuum limit at infinite volume), the
parameters is reduced and the fits are further constrained. extrapolations, detailed in Sec. III F, are performed for
To ensure that the excited state contributions are suffi- the results obtained from fits 1–4 (ES ¼ 1) and fits 5–8
ciently under control, we carry out a variety of different fits, (ES ¼ 2), separately. For each set of fits, 500 samples are
summarized in Table II. We vary drawn from the combined bootstrap distributions of the four
(i) The datasets included in the fit: simultaneous fits
are performed to the data for J ∈ fA; S; T; Vg and 3
Due to nj ≥ 2 and its quantization, δt1 and δt2 depend slightly
on the lattice spacing: δtj ≈ 0.20 fm; 0.29 fm (β ¼ 3.34), δtj ≈
2
Given the precision of the two-point function relative to that 0.17 fm; 0.26 fm (β ¼ 3.40), δtj ≈ 0.15 fm; 0.23 fm (β ¼ 3.46),
of the three-point function, this strategy is very similar to fitting δtj ≈ 0.13 fm; 0.26 fm (β ¼ 3.55), δtj ≈ 0.15 fm; 0.25 fm
CB2pt and CB3pt simultaneously. (β ¼ 3.70), δtj ≈ 0.16 fm; 0.27 fm (β ¼ 3.85).
034512-10
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
FIG. 9. Unrenormalized ratios RΞJ ðt; τÞ, J ∈ fA; S; T; Vg [defined in Eq. (26)] for the cascade baryon on ensemble N302
(M π ¼ 348 MeV and a ¼ 0.049 fm), where t ≈ f0.7; 0.8; 1.0; 1.2g fm. The gray horizontal lines and bands show the results for
the ground state matrix elements h0jūΓJ u − d̄ΓJ dj0i ¼ gΞ;latt
J , obtained from a simultaneous fit to the ratios for all channels and source-
sink separations using parametrization 7 [see Eq. (27) and Table II]. The data points with τ ∈ ½δt; t − δt, where δt ¼ 2a, are included in
the fit (the faded data points are omitted), which is the maximum fit range possible for our action. The colored curves show the
expectation from the fit for each source-sink separation.
FIG. 10. Results for the four unrenormalized charges of the nucleon (top) and cascade baryon (bottom) obtained from the fits listed in
Table II for ensemble N302 (M π ¼ 348 MeV and a ¼ 0.049 fm). The green (blue) horizontal lines and bands indicate the final results
and errors obtained from the median and 68% confidence level interval of the combined bootstrap distributions determined from the fits
indicated by the green (blue) data points which include one (two) excited state(s). On the right the energy gaps determined in the fits and
those corresponding to the lowest lying multiparticle states are displayed using the same color coding as in Fig. 11.
034512-11
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
9 A, S, T δt1 2 ΔE1 bJ2 , bJ4 FIG. 11. Results for the first and second excited state energy
10 A, S, T δt2 2 ΔE1 bJ2 , bJ4 gaps of the cascade baryon, ΔE1 (brown data points) and ΔE2
11 A, S, T, V δt1 2 ΔE1 bJ2 , bJ4 (orange data points), respectively, determined on ensembles lying
12 A, S, T, V δt2 2 ΔE1 bJ2 , bJ4 on the tr M ¼ const trajectory with a ¼ 0.064 fm. The pion mass
δt1 ΔE2 decreases from left to right with Mπ ¼ 414 MeV for ensemble
13 A, S, T 2 bT1 , bJ2 , bJ4
N202 and M π ¼ 202 MeV for ensemble D200, see Table I. For
14 A, S, T δt2 2 ΔE2 bT1 , bJ2 , bJ4
each ensemble, the ΔE1 obtained using fits 1–4 of Table II are
15 A, S, T, V δt1 2 ΔE2 bT;V J
1 , b2 , b4
J
shown on the left and the ΔE1 (fixed with a prior to the lowest
16 A, S, T, V δt2 2 ΔE2 T;V
b1 , b2 , bJ4
J
multi-particle energy gap) and ΔE2 resulting from fits 5–8 are
17 A, S, T δt1 2 ΔE2 bJ2 , bJ4 displayed on the right. For comparison, the energy gaps of the
18 A, S, T δt2 2 ΔE2 bJ2 , bJ4 lower lying noninteracting multiparticle states with the quantum
19 A, S, T, V δt1 2 ΔE2 bJ2 , bJ4 numbers of the cascade baryon are shown as horizontal lines,
20 A, S, T, V δt2 2 ΔE2 bJ2 , bJ4 where the momenta utilized for each hadron are indicated in
lattice units.
fit variations. The final result and error, shown as the green two excited states. The ΔE2 energy gaps from the two
and blue bands in Fig. 10, correspond to the median and the excited state fits (with the first excited state fixed with a
68% confidence interval, respectively. Note that we take the prior to the lowest multi-particle level) are consistent with
same 500 bootstrap samples for all the baryons to preserve the next level that is significantly above the first excited
correlations. The final results for all the ensembles are listed state, although for ensemble D200 the errors are too large to
in Tables XVI–XVIII of Appendix B for the nucleon, sigma draw a conclusion. Given that more than one excited state is
and cascade baryons, respectively. contributing significantly, we expect that the latter fits
In terms of the energy gaps extracted, Fig. 10 shows that isolate the ground state contribution more reliably. We
we find consistency across variations in the fit range and remark that within present statistics, two-exponential fits to
whether the vector channel data is included or not. the two-point functions alone give energy gaps aΔE ¼
However, the first excited energy gap ΔE1 obtained 0.390ð37Þ; 0.371ð34Þ; 0.430ð37Þ and 0.312(46) for N202,
from the single excited state fits tends to be higher than N203, N200 and D200, respectively, that are all larger than
the lowest multiparticle level, in particular, as the pion mass
is decreased, suggesting that contributions from higher
excited states are significant. This can be seen in Fig. 11,
where we compare the results for the energy gaps for the TABLE III. Improvement coefficients bJ for J ∈ fA; S; T; Vg
cascade baryon with the lower lying noninteracting Ξπ and from Refs. [54,55].
Ξππ states for four ensembles with a ¼ 0.064 fm and pion
β bA bS bT bV
masses ranging from 414 MeV down to 202 MeV. Note that
the multi-particle levels are modified in a finite volume, 3.34 1.249(16) 1.622(47) 1.471(11) 1.456(11)
although the corresponding energy shifts may be small for 3.4 1.244(16) 1.583(62) 1.4155(48) 1.428(11)
the large volumes realized here. There are a number of 3.46 1.239(15) 1.567(74) 1.367(12) 1.410(13)
3.55 1.232(15) 1.606(98) 1.283(14) 1.388(17)
levels within roughly 500 MeV of the first excited state.
3.7 1.221(13) 1.49(11) 1.125(15) 1.309(22)
Some levels lie close to each other and one would not 3.85 1.211(12) 1.33(16) 0.977(38) 1.247(26)
expect that the difference can be resolved by fits with one or
034512-12
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
aΔE2 , with the exception of D200, where the two gaps For the renormalization factors, we employ the values
agree within errors. obtained in Ref. [56]. The factors are determined non-
perturbatively in the RI0 -SMOM scheme [57,58] and then
E. Nonperturbative renormalization (for ZS and ZT ) converted to the MS scheme using three-
and improvement loop matching [59–61]. We remark that the techniques for
implementing the Rome-Southampton method were
The isovector lattice charges, gB;latt
J , extracted in the extended in Ref. [56] to ensembles with open boundary
previous subsection need to be matched to the continuum conditions in time. This development enables us to utilize
MS scheme. The renormalized matrix elements suffer ensembles with a < 0.06 fm, where only open boundary
from discretization effects, however, the leading order conditions in time are available due to the need to maintain
effects are reduced to Oða2 Þ when implementing full ergodicity. A number of different methods are employed in
OðaÞ improvement. In the forward limit, in addition to Ref. [56] to determine the renormalization factors. In order
using a nonperturbatively OðaÞ improved fermion action, to assess the systematic uncertainty arising from the
this involves taking mass dependent terms into account. matching in the final results for the charges at the physical
The following multiplicative factors are applied, point in the continuum limit, we make use of two sets of
results, collected in Tables V and VI and referred to as Z1J
gBJ ¼ ZkJ ð1 þ aml bJ þ 3am̄b̃J ÞgB;latt
J þ Oða2 Þ; ð28Þ and Z2J , respectively, in the following. The first set of results
are extracted using the fixed-scale method, where the
for J ∈ fV; A; S; Tg, where ZJ are the renormalization RI0 -SMOM factors are determined at a fixed scale (ignoring
factors and bJ and b̃J are the OðaÞ improvement coef- discretization effects), while the second set are obtained by
ficients. Note that the renormalization factors for the scalar fitting the factors as a function of the scale and the lattice
and tensor currents depend on the scale, ZS;T ¼ ZS;T ðμÞ, spacing, the “fit method.” See Ref. [56] for further details.
where we take μ ¼ 2 GeV. The vector Ward identity lattice In both cases, lattice artefacts are reduced by subtracting
quark mass amq is obtained from the hopping parameter κq the perturbative one-loop expectation. For the axial and
(q ¼ l; s) and the critical hopping parameter κcrit via
amq ¼ ð1=κq − 1=κ crit Þ=2. m̄ ¼ ð2ml þ ms Þ=3 denotes TABLE V. Set of renormalization factors taken from Ref. [56],
the flavor averaged quark mass. The hopping parameters denoted as Z1J in the text. The factors are determined using
for all ensembles used within this work are tabulated in the RI0 -SMOM scheme and the “fixed-scale method” with the
Table XV of Appendix B. For κcrit we utilize the inter- perturbative subtraction of lattice artefacts. For ZA and ZV , the
polation formula [17] values correspond to those listed under Z0A and Z0V , respectively,
which are obtained using renormalization conditions consistent
with the respective Ward identities. The statistical and systematic
1 1 þ 0.28955g2 − 0.1660g6
¼ 8 − 0.402454g2 : ð29Þ errors have been added in quadrature.
κcrit 1 þ 0.22770g2 − 0.2540g4
β ZA S ð2 GeVÞ ZT ð2 GeVÞ
ZMS MS
ZV
The improvement coefficients bJ and b̃J are determined 3.34 0.77610(58) 0.6072(26) 0.8443(35) 0.72690(71)
nonperturbatively in Ref. [54]. We make use of updated 3.4 0.77940(36) 0.6027(25) 0.8560(35) 0.73290(67)
preliminary values, which will appear in a future publica- 3.46 0.78240(32) 0.5985(25) 0.8665(36) 0.73870(71)
tion [55]. These are listed in Tables III and IV, respectively. 3.55 0.78740(22) 0.5930(25) 0.8820(37) 0.74740(82)
Note that no estimates of b̃J are available for β ¼ 3.85. 3.7 0.79560(98) 0.5846(24) 0.9055(42) 0.76150(94)
Considering the size of the statistical errors, the general 3.85 0.8040(13) 0.5764(25) 0.9276(42) 0.77430(76)
reduction of the jb̃J j values with increasing β (and the
decreasing a), at this lattice spacing we set b̃J ¼ 0 for all J.
TABLE VI. Set of renormalization factors denoted as Z2J in the
TABLE IV. Improvement coefficients b̃J for J ∈ fA; S; T; Vg text. These are determined as in Table V but now using the “fit
from Refs. [54,55]. Note that no results are available for method.”
β ¼ 3.85.
β ZA S ð2 GeVÞ
ZMS T ð2 GeVÞ
ZMS ZV
β b̃A b̃S b̃T b̃V
3.34 0.7579(42) 0.6115(93) 0.8321(95) 0.7072(60)
3.34 −0.06ð28Þ −0.24ð55Þ 1.02(16) 1.05(13) 3.4 0.7641(35) 0.6068(86) 0.8462(88) 0.7168(49)
3.4 −0.11ð13Þ −0.36ð23Þ 0.49(17) 0.41(11) 3.46 0.7695(36) 0.6025(79) 0.8585(84) 0.7250(43)
3.46 0.08(11) −0.421ð83Þ 0.115(19) 0.158(28) 3.55 0.7774(36) 0.5968(66) 0.8756(76) 0.7367(37)
3.55 −0.03ð13Þ −0.25ð12Þ 0.000(37) 0.069(42) 3.7 0.7895(32) 0.5880(45) 0.9010(63) 0.7544(30)
3.7 −0.047ð75Þ −0.274ð65Þ −0.0382ð60Þ −0.031ð18Þ 3.85 0.8006(25) 0.5793(35) 0.9243(55) 0.7699(38)
034512-13
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
TABLE VII. Renormalization factors ZA and ZV obtained from works, see, e.g., Refs. [68–70]. Finite volume effects
the interpolation formulas in Eqs. (C.7) and (C.8) in Ref. [63], appear at Oðp3 Þ with no additional LECs appearing in
denoted as Z3J in the text. the coefficients. Again the signs of the corrections are the
β ZA ZV
opposite to the trend seen in the data and, when included, it
is difficult to resolve the effects of the decuplet baryons. As
3.34 0.7510(11) 0.7154(11) is shown in Sec. IV, the data for all the charges are well
3.4 0.75629(65) 0.72221(65) described when the fit form is restricted to the dominant
3.46 0.76172(39) 0.72898(39) terms, with free coefficients c0 , cπ , cK and cV .
3.55 0.76994(34) 0.73905(35)
We remark that the same set of LECs appear in the Oðp2 Þ
3.7 0.78356(32) 0.75538(33)
3.85 0.79675(45) 0.77089(47) SU(3) ChPT expressions for the three different octet
baryons (for a particular charge). Ideally, one would carry
out a simultaneous fit to the whole baryon octet (taking the
vector currents, we also consider a third set of renormal- correlations between the gBJ determined on the same
ization factors, Z3J , listed in Table VII, that are obtained ensemble into account). However, we obtain very similar
with the chirally rotated Schrödinger functional approach results when fitting the gBJ individually compared to fitting
[62], see Ref. [63]. We emphasize that employing the the results for all the octet baryons simultaneously. For
different sets of renormalization factors should lead to simplicity, we choose to do the former, such that the
consistent results for the charges in the continuum limit. coefficients cX for the different baryons are independent of
one another.
F. Extrapolation strategy Lattice spacing effects also need to be taken into account
and we add both mass independent and mass dependent
In the final step of the analysis the renormalized charges terms to the continuum fit ansatz to give
gBJ determined at unphysical quark masses and finite lattice
spacing and spatial volume are extrapolated to the physical
point in the continuum and infinite volume limits. We gBJ ðM π ; MK ; L; aÞ ¼ gBJ ðMπ ; MK ; L; 0Þ
employ a similar strategy to the one outlined in Ref. [64] þ ca a2 þ c̄a M̄2 a2 þ δca δM2 a2
and choose continuum fit functions of the form
þ ca;3 a3 ; ð31Þ
gBJ ðM π ; MK ; L; a ¼ 0Þ
e−LMπ where M̄2 ¼ ð2M2K þ M2π Þ=3 and δM2 ¼ M2K − M2π .
¼ c0 þ cπ M 2π þ cK M2K þ cV M 2π pffiffiffiffiffiffiffiffiffiffi ; ð30Þ
LMπ The meson masses areprescaledffiffiffiffiffiffi with the Wilson flow
scale t0 [71], Mπ;K ¼ 8t0 M π;K to form dimensionless
to parameterize the quark mass and finite volume depend- combinations. This rescaling is required to implement full
ence, where L is the spatial lattice extent and the coef- OðaÞ improvement (along with employing a fermion action
ficients cX , X ∈ f0; π; K; Vg are understood to depend on and isovector currents that are nonperturbatively OðaÞ
the baryon B and the current J. The leading order improved) when simulating at fixed bare lattice coupling
coefficients c0 give the charges in the SU(3) chiral limit, instead of at fixed lattice spacing, see Sec. 4.1 of Ref. [17]
which can be expressed in terms of two LECs, e.g., F and for a detailed discussion of this issue. The values of t0 =a2
D, for the axial charges, see Eq. (7). and the pion and kaon masses in lattice units for our
Equation (30) is a phenomenological fit form based on set of ensembles are given in Table XV of Appendix B.
the SU(3) ChPT expressions for the axial charge. It contains We translate between different lattice spacings using t⋆0 ,
the expected Oðp2 Þ terms for the quark mass dependence the value of t0 along the symmetric line where 12t0 M2π ¼
and the dominant finite volume corrections. The Oðp3 Þ pffiffiffiffiffiffiffi
1.110 [72], i.e., a ¼ a= 8t⋆0 . The values, determined in
expressions for gBA [65–67] contain log terms with coef-
Ref. [17], are listed in Table VIII. Note that we include a
ficients completely determined by the LECs F and D. In an
earlier study of the axial charges on the ms ¼ ml subset of term that is cubic in the lattice spacing in the fit form,
the ensembles used here [64], we found that including these however, this term is only utilized in the analysis of the
terms did not provide a satisfactory description of the data. vector charge, for which we have the most precise data.
When terms arising from loop corrections that contain
decuplet baryons are taken into account, additional LECs TABLE VIII. Values for t⋆0 =a2 at each β-value as determined in
enter that are difficult to resolve. If the coefficient of the log Ref. [17].
term is left as a free parameter, one finds that the coefficient
has the opposite sign to the ChPT expectation without β 3.34 3.4 3.46 3.55 3.7 3.85
decuplet loops. We made similar observations in this study t⋆0 2.204(6) 2.888(8) 3.686(11) 5.157(15) 8.617(22) 13.988(34)
a2
and this is also consistent with the findings of previous
034512-14
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
To obtain results at the physical quark mass point, we IV. EXTRAPOLATIONS TO THE CONTINUUM,
make use of the scale setting parameter INFINITE VOLUME, PHYSICAL QUARK
MASS LIMIT
qffiffiffiffiffiffiffiffiffiffiffiffiffiffi
ð20Þ
8t0;phys ¼ 0.4098ð25Þ fm; ð32Þ We present the extrapolations to the physical point in the
continuum and infinite volume limits of the isovector
vector (V), axial (A), scalar (S) and tensor (T) charges
determined in Ref. [17] and take the isospin corrected pion
for the nucleon (N), sigma (Σ) and cascade (Ξ) octet
and kaon masses quoted in the FLAG 16 review [73] to
baryons.
define the physical point in the quark mass plane,
034512-15
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
pffiffiffiffiffiffi
parameter t0, while the lattice spacing a ¼ a= 8t0 . See the
previous section for further details of the extrapolation
procedure. We implement full OðaÞ improvement and
leading discretization effects are quadratic in the lattice
spacing. However, the data for the nucleon vector charge
are statistically very precise and higher order effects can be
resolved. This motivates the addition of the cubic term in
Eq. (36). The data for gΣV and gΞV are less precise as they are
determined employing the stochastic approach outlined in
Sec. III C which introduces additional noise, see Appendix A
for further discussion.
The data are well described by Eq. (36), as demonstrated
by the fit, shown in Fig. 13, for gNV which has a goodness of FIG. 14. Results for the nucleon vector charge gNV in the
fit of χ 2 =N dof ¼ 0.92. The data are extracted using two continuum limit at the physical point obtained using Z3V (see
excited states in the fitting analysis (see Sec. III D) and we Table VII) and five different parametrizations applied to dataset
employ the most precise determination of the renormaliza- DSðM <400
π
MeV
Þ (fits to Eq. (31) with different coefficients set to
tion factors (Z3V , see Table VII). A cut of M π < 400 MeV is zero labeled 1, …, 5, see the text) and DSðM <400 π
MeV ; a<0.1 fm Þ
imposed on the ensembles entering the fit, however, fits (6, …, 10). See Sec. III F for the definitions of the data sets. The
including all data points are also performed, as detailed data were extracted including two excited states in the fitting
below. When the data are corrected for the discretization analysis, see Sec. III D. The model average is shown as the red
effects according to the fit, we see consistency with gNV ¼ 1, data point and the green horizontal line and band. On the right the
model averaged distribution is displayed as a histogram where
for all pion and kaon masses. Using the fit to shift the data
also the median and the 68% confidence level interval, which
points to the physical point, we observe that the lattice form the final result, are indicated (green lines). The top panel
shows the weights (gray points) assigned to the individual fits,
with the corresponding χ 2 =N dof values given above.
034512-16
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
TABLE IX. Results for gBV , B ∈ fN; Σ; Ξg, obtained with three within a combined statistical and systematic uncertainty of
different sets of renormalization factors. The errors include the about 1‰.
statistical and all the systematic uncertainties. The above analysis is also performed utilizing the sets
of renormalization factors Z1V and Z2V , determined via the
Renormalization gNV gΣV gΞV
RI0 -SMOM scheme [56]. The results for the nucleon vector
Z1V (Table V) 0.9975ð20Þ
ð22Þ
2.012ð16Þ
ð26Þ
1.012ð11Þ
ð13Þ
charge are compared in Fig. 15. The uncertainties on these
Z2V (Table VI) 0.9945ð41Þ
ð66Þ
2.014ð16Þ
ð36Þ
1.008ð8Þ
ð19Þ factors are larger, in particular for Z2V , than those of set Z3V ,
which is derived using the chirally rotated Schrödinger
Z3V (Table VII) 1.0012ð11Þ
ð12Þ
2.021ð27Þ
ð21Þ
1.015ð11Þ
ð10Þ
functional approach [63]. This translates into larger errors
for gNV for those fits. The lattice spacing dependence is
model averaging procedure) give values that are incon- somewhat different: the first quadratic mass dependent term
sistent with the continuum expectation. The results are in Eq. (36) and the cubic term can no longer be fully
stable under the removal of the coarsest ensembles. resolved and also parametrization ð2; fc0 ; ca ; δca gÞ gives a
Performing the model averaging procedure, the final χ 2 =N dof ¼ 1.00 (0.95) when employing Z1V (Z2V ).
result for the nucleon, given in the last row of the first The systematic uncertainty of the results due to residual
column of Table IX, agrees with the expectation gNV ¼ 1 excited state contamination and the range of pion masses
FIG. 15. Results for the nucleon vector charge gNV as in Fig. 14 but now also including those obtained employing Z1V and Z2V .
FIG. 16. Overview of the results for the vector charges gBV , B ∈ fN; Σ; Ξg, obtained from different datasets. These are labeled by the
number of excited states used in the fitting analysis (ES ¼ 1, 2), the pion mass cut imposed (denoted A or B) and the set of
renormalization factors employed (ZkV , k ¼ 1, 2, 3). Each data point represents a model averaged result. The label A indicates that 15 fits
[5 fit variations applied to three datasets, DS(0), DSðM π<400 MeV Þ and DSða<0.1 fm Þ] are averaged, while the results labeled with B are
based on the set of 10 fits utilized in Fig. 14 (5 fit variations applied to two datasets, DSðMπ<400 MeV Þ and DS[M <400
π
MeV , a<0.1 fm )]. See
k
Sec. III F for the definitions of the datasets. The final results for each ZV (filled squares) are listed in Table IX.
034512-17
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
employed in the extrapolations is also considered. employing the parametrization in Eq. (31) [with the
Figure 16 shows the model averaged results discussed so continuum form in Eq. (30)]. The discretization effects
far, displayed as filled squares, and also those obtained are found to be fairly mild and we are not able to resolve the
using several other sets of fits. These are labeled in terms of quadratic mass dependent terms or a cubic term. These
the number of excited states (one or two) included in the terms are omitted throughout. As already mentioned in
fitting analysis, the cuts imposed on the pion mass (A or B) Sec. III F we are also not able to resolve any higher order
and the renormalization factors utilized. For the results ChPT terms in the continuum parametrization.
from pion mass cut A, 15 fits enter the model average, the A five parameter fit, with free coefficients
five different parametrizations are applied to three datasets fc0 ; cπ ; cK ; cV ; ca g, describes the data well, as demon-
DS(0), DSðM<400 π
MeV
Þ and DSða<0.1 fm Þ. Note that the first strated in Fig. 17 for the nucleon (with χ 2 =N dof ¼ 0.86) and
and third dataset include ensembles with pion masses up to
430 MeV. For mass cut B, datasets DSðM <400 π
MeV
Þ and DS
<400 MeV <0.1 fm
[M π ,a )] are used, giving 10 fits in total. The
results only depend on the choice of renormalization
factors, suggesting that the systematic uncertainties due
to excited state contamination and the cut made on the pion
mass are very small.
Repeating the whole procedure for the sigma and the
cascade baryons gives vector charges which are also con-
sistent with the expected values to within 1.5σ, as shown in
Fig. 16 (see Figs. 41 and 42 in Appendix C for the individual
fits for mass cut B). The statistical noise introduced by the
stochastic approach dominates, leading to much less precise
values and very little variation between the results for the
different hyperon data sets. We take the values obtained from
the datasets (2, B, ZkV ), listed in Table IX, as our estimates of
the vector charges as these datasets give the most reliable
determinations of the charges across the different channels
(as discussed in the following subsections).
Overall, the results demonstrate that the systematics
arising from excited state contamination, renormalization
and finite lattice spacing are under control in our analysis in
this channel (to within an error of 1‰ for the nucleon).
B. Axial charges
In the following we present our results for the
nucleon, sigma and cascade isovector axial charges gBA , FIG. 17. Simultaneous quark mass, continuum and finite vol-
B ∈ fN; Σ; Ξg. The nucleon axial charge is very precisely ume extrapolation of the nucleon isovector axial charge gNA
extracted on ensembles with Mπ < 400 MeV using two excited
measured in experiment, λ ¼ gNA =gNV ¼ 1.2754ð13Þ [75],
states in the fitting analysis (see Sec. III D) and renormalization
and serves as another benchmark quantity when assessing factors Z3A (see Table VII). A five parameter fit form is employed,
the size of the systematics of the final results. Note, see the text. (Top) Pion mass dependence of gNA , where the data
however, that possible differences of up to 2%, due to points are corrected, using the fit, for finite volume and discre-
radiative corrections, between λ computed in QCD and an tization effects and shifted (depending on the ensemble) to kaon
effective λ measured in experiment have been discussed masses corresponding to the tr M ¼ const and ms ¼ const tra-
recently [76,77]. Lattice determinations of gNA are known to jectories. The fit is shown as a gray band with the three trajectories
be sensitive to excited state contributions, finite volume distinguished by blue (tr M ¼ const, circles), green (ms ¼ const,
effects and other systematics. Whereas there is a long diamonds) and orange (ms ¼ ml , triangles) lines, respectively.
history of lattice QCD calculations of gNA , see, e.g., the The vertical dashed line indicates the physical point. (Middle)
Lattice spacing dependence at the physical point in the infinite
FLAG 21 review [78], there are very few lattice compu- volume limit. (Bottom) Finite volume dependence at the physical
tations of hyperon axial charges [22–26] and only few point in the continuum limit. The dashed blue line (band) indicates
phenomenological estimates exist from measurements of the infinite volume result. For better visibility, the data points for
semileptonic hyperon decay rates. ensembles D150, E250, and D452, which have relatively large
We carry out simultaneous continuum, quark mass and errors (see Table XVI), are not displayed. The black cross at the
finite volume fits to the individual baryon charges physical point indicates the experimental value [75].
034512-18
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
FIG. 18. The same as Fig. 17 for the isovector axial charges gBA of the sigma baryon (left) and the cascade baryon (right). For better
visibility, the data points for ensembles D452 and D451, which have relatively large errors (see Table XVIII), are not displayed for the
cascade baryon. Compared to the nucleon, for the analysis of the hyperon charges a reduced set of ensembles is employed, see
Tables XVII and XVIII for a complete list of ensembles.
Fig. 18 for the sigma and cascade baryons (with χ 2 =N dof ¼ The results of the eight fits and their model averages for
0.85 and 1.25, respectively). The data are extracted using the three different determinations of the renormalization
two excited states (ES ¼ 2) in the fitting analysis (see factors are shown in Fig. 19 for the nucleon and in Fig. 43
Sec. III D) and renormalized with factors Z3A (that are the of Appendix C for the hyperon axial charges. In all cases,
most precise of the three determinations considered, see we find consistent results across the different fits and choice
Table VII). For the cascade baryon, with two strange of renormalization factor suggesting that the statistical
quarks, the data on the three quark mass trajectories errors dominate. The additional lattice spacing term is
(tr M ¼ const, ms ¼ const and ml ¼ ms ) are clearly delin- not really resolved with the goodness of fit only changing
eated, however, note the different scale on the right of slightly, while the errors on the coefficients increase. For
Fig. 18. The availability of ensembles on two trajectories the nucleon and sigma baryon, all fits have a χ 2 =N dof < 1
which intersect at the physical point helps to constrain the and are given a similar weight in the model average, while
physical value of the axial charge. In terms of the finite for the cascade baryon, the cut M<300 π
MeV
is needed to
volume effects, only the nucleon shows a significant achieve a goodness of fit around 1 and these fits have the
dependence on the spatial extent. The quark mass depend- highest weight factors.
ence is also pronounced in this case. In order to further explore the systematics, additional
As in the vector case, we quantify the systematics datasets are considered. We assess the sensitivity of the
associated with the extraction of the charges at the physical results to excited state contributions by performing extrap-
point (in the continuum and infinite volume limits) by olations of the data extracted using only one excited state
varying the parametrization and the set of ensembles that (“ES=1”) in the fitting analysis. In addition, as only the
are included in the fit. We consider two fit forms Oðp2 Þ ChPT terms are included in the continuum para-
ð1; fc0 ; cπ ; cK ; cV ; ca gÞ and ð2; fc0 ; cπ ; cK ; cV ; ca ; δca g) metrization, we test the description of the quark mass
and four datasets, DS(M<400 π
MeV
), DSðM <300 π
MeV
Þ, dependence by performing 10 fits, involving the two
DSðM<400
π
MeV <0.1 fm
; a Þ and DSðM <400 MeV
π ; LM >4
π Þ, see parametrization variations above, applied to five datasets,
Sec. III F for their definitions. DS(0), DSðM<400π
MeV
Þ, DSðM <300
π
MeV
Þ, DSða<0.1 fm Þ and
034512-19
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
FIG. 19. The same as Fig. 15 for the nucleon axial charge gNA . The eight fits correspond to two fit variations, see the text, applied to four
datasets, DSðM <400
π
MeV
Þ, DSðM π<300 MeV Þ, DSðM<400
π
MeV <0.1 fm
;a Þ and DSðM <400π
MeV
π Þ. The data are extracted using two
; LM >4
excited states in the fitting analysis, see Sec. III D.
DSðLM>4 π Þ. The first, fourth and fifth datasets include As detailed in Sec. III D, more than one excited state is
ensembles with pion masses up to 430 MeV. contributing significantly to the ratio of three-point over
The results for the axial charges from model averaging two-point correlation functions and including two excited
the 10 fits (denoted A) employing the 5 datasets and also states in the fitting analysis enables the ground state matrix
from the 8 fits (denoted B) using the 4 datasets given above, element to be isolated more reliably. Considering the pion
for the ES ¼ 1 and ES ¼ 2 data and the different renorm- mass cuts, to be conservative we take the results of the
alization factors are displayed in Fig. 20. Very little model averages of the B datasets (where all the ensembles
variation is seen in the results in terms of the range of have M π < 400 MeV) as only the dominant mass depen-
pion masses included and, as before, the renormalization dent terms are included in the continuum parametrization.
factors employed, suggesting the associated systematics are Our estimates, corresponding to the (ES ¼ 2, B, ZkA ) results
accounted for within the combined statistical and system- in Fig. 20, are listed in Table X.
atic error (which includes the uncertainty due to lattice
spacing and finite volume effects). However, the results are
C. Scalar charges
sensitive to the number of excited states included in the
fitting analysis. This is only a significant effect for the As there is no isovector scalar current interaction at tree-
nucleon, for which the ES ¼ 1 results lie around 2.5σ level in the Standard Model, the scalar charges cannot be
below experiment. Similar underestimates of gNA have been measured directly in experiment. However, the conserved
observed in many earlier lattice studies [78]. vector current (CVC) relation can be used to estimate the
FIG. 20. The same as Fig. 16 for the nucleon, sigma and cascade axial charges. The label A indicates that 10 fits enter the model
average corresponding to two fit variations, see the text, applied to 5 datasets, DS(0), DSðM <400
π
MeV
Þ, DSðM π<300 MeV Þ, DSða<0.1 fm Þ and
DSðLM π Þ. For the data points labeled with B, the results of the 8 fits employed in Fig. 19 are averaged. For the nucleon, the FLAG 21
>4
average for N f ¼ 2 þ 1 [53,79] and the experimental value [75] are indicated (black diamonds).
034512-20
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
034512-21
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
FIG. 22. The same as Fig. 16 for the nucleon, sigma and cascade scalar charges. The label A indicates that 8 fits enter the model
average corresponding to two fit variations, see the text, applied to 4 datasets, DS(0), DSðM π<400 MeV Þ, DSða<0.1 fm Þ and DSðLM >4π Þ. For
the data points labeled with B, the two fit variations are performed on 3 datasets, DSðM <400 π
MeV Þ, DSðM <400 MeV ; a<0.1 fm Þ and
π
DSðM <400 MeV ; LM >4 Þ, giving a total of 6 fits for the average. For the nucleon, the FLAG 21 result for N ¼ 2 þ 1 [53] is also shown
π π f
(black diamond).
D. Tensor charges fits have very large uncertainties, unless some assumptions
In the isosymmetric limit, the nucleon tensor charge is are made. In fact, in some analyses, the fit is constrained to
equal to the first moment of the nucleon isovector trans- reproduce the lattice results for the isovector charge,
versity parton distribution function. Due to the lack of see Refs. [14,15]. The FLAG 21 review [78] gives as
experimental data, estimates of gNT from phenomenological the N f ¼ 2 þ 1 value for the nucleon tensor charge the
result of Ref. [53], gNT ¼ 0.965ð61Þ, whereas, as far as we
know, there is only one previous study of the hyperon
tensor charges [26].
The extraction of the octet baryon tensor charges at the
physical point follows the analysis of the axial charges in
Sec. IV B. In particular, the parametrizations employed
and the datasets considered are the same. Figure 23
displays a typical example of an extrapolation for the
nucleon tensor charge for a five parameter fit with a
χ 2 =N dof ¼ 0.63. See Fig. 47 in Appendix C for the
analogous figures for the sigma and cascade baryons.
The variation of the fits with the parametrization and the
datasets utilized and the corresponding model averages,
for the datasets with pion mass cut B (see Sec. IV B), are
shown in Fig. 45.
An overview of the model averaged results for all
variations of the input data is given in Fig. 24. The
agreement between the different determinations suggests
the systematics associated with the extrapolation are under
control. Although the results utilizing data extracted with
two excited states (ES ¼ 2) in the fitting analysis are
consistently above or below those extracted from the
ES ¼ 1 data, considering the size of the errors of the
model averages (which combine the statistical and system-
atic uncertainties), the differences are not significant. Our
estimates for the tensor charges, corresponding to the
FIG. 23. The same as Fig. 17 for the nucleon tensor charge gNT . (ES ¼ 2, B, ZkT ) values, are listed in Table XII.
The factors Z1T are used for the matching (see Table V). For
orientation, the FLAG 21 result for N f ¼ 2 þ 1 [53] is indicated V. DISCUSSION OF THE RESULTS
(black diamond) at the physical point. For better visibility, the
data points for ensembles D150 and D452, which have relatively Our values for the vector, axial, scalar and tensor charges
large errors (see Table XVI), are not displayed. of the nucleon, sigma and cascade baryons are given in
034512-22
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
FIG. 24. The same as Fig. 16 for the nucleon, sigma and cascade tensor charges. For the nucleon, the FLAG 21 result for N f ¼ 2 þ 1
[53] is also shown (black diamond).
Tables IX–XII, respectively. In each case, we take the most dynamical fermions employed, including charm quarks in
precise value as our final result, i.e., the one obtained using the sea is not expected to lead to a discernible effect.
Z3V and Z3A for the vector and axial channels, respectively, Regarding the hyperon axial charges, far fewer works
and Z1S and Z1T for the scalar and the tensor. In the following exist. Lin et al. [22,92] performed the first study, utiliz-
we compare with previous determinations of the charges ing N f ¼ 2 þ 1 ensembles with pion masses ranging
taken from the literature and discuss the SU(3) flavor between 350 MeV and 750 MeV and a single lattice
symmetry breaking effects in the different channels. We use spacing of 0.12 fm. After an extrapolation to the physical
the conserved vector current relation and our result for the pion mass they obtain gΣA ¼ 0.900ð42Þstat ð54Þsys and gΞA ¼
scalar charge of the sigma baryon to determine the up −0.277ð15Þstat ð19Þsys , where estimates of finite volume and
and down quark mass difference. We compute the QCD discretization effects are included in the systematic uncer-
contributions to baryon isospin mass splittings and evaluate tainty. Note that we have multiplied their result for gΣA by
the isospin breaking effects on the pion baryon σ terms. a factor of two to match our normalization convention.
In Refs. [24,93] ETMC determined all octet and decuplet
A. Individual charges
We first consider the axial charges. Our final values read
TABLE XII. Results for gBT, B ∈ fN; Σ; Ξg, obtained with two
different sets of renormalization factors. The errors include the
statistical and all the systematic uncertainties.
FIG. 25. Compilation of recent lattice determinations of the
Renormalization gNT gΣT gΞT nucleon axial charge gNA with N f ¼ 2 þ 1 [26,53,79,81–86] and
N f ¼ 2 þ 1 þ 1 [68,69,87–91] dynamical fermions. Values with
Z1T (Table V) 0.984ð29Þ
ð19Þ
0.798ð21Þ
ð15Þ
−0.1872ð41Þ
ð59Þ
filled symbols were obtained via a chiral, continuum and finite
Z2T (Table VI) 0.979ð27Þ
ð19Þ
0.793ð21Þ
ð17Þ
−0.1872ð42Þ
ð59Þ
volume extrapolation. The vertical black line gives the exper-
imental result [75].
034512-23
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
034512-24
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
mB B
gBP ¼ g ; ð38Þ
ml A
034512-25
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
034512-26
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
renormalization factors and improvement coefficients and, carry out the same comparison for the tensor charges.
possibly, other systematics cancel from Eq. (43). For the In this case, within errors of a few percent, no flavor
ratio of the Ξ over the Σ axial charge we see no significant symmetry violation is seen. Moreover, DT ðmÞ=FT ðmÞ ¼
difference between the physical point value and that DT ð0Þ=FT ð0Þ within errors.
obtained for the same average quark mass at the flavor In order to quantify the symmetry breaking effect
symmetric point. The symmetry breaking effect of the between matrix elements involving the current J as a
combination involving gNA =gΣA can be attributed to the pion function of the quark mass splitting ms − ml , we define
mass dependence of gNA , see Fig. 17. The red symbols at the
physical point (dashed vertical line) correspond to our gΞJ þ gNJ − gΣJ
δJSUð3Þ ¼ ; ð44Þ
continuum, infinite volume limit extrapolated results, listed gΞJ þ gNJ þ gΣJ
in Table XIII for the combinations Eq. (43).
In Fig. 31 the combinations Eq. (43) are shown for the where for ms ¼ ml, δJSUð3Þ ¼ ð2FJ −2FJ Þ=ð2FJ þ2FJ Þ ¼ 0,
isovector scalar charges. These are compared to our SU(3)
see Eqs. (8)–(10). Also from these ratios some of the
chiral limit extrapolated results (yellow bands) and the
systematics as well as the renormalization factors and
continuum, infinite volume limit results at the physical
improvement terms will cancel. We define a dimensionless
point (red diamonds). We find no statistically significant
SU(3) breaking parameter x ¼ ðM2K −M 2π Þ=ð2M2K þM 2π Þ ∼
symmetry breaking in this case. However, the statistical
ms −ml and assume a polynomial dependence:
errors are larger than for the axial case and also FS > FA .
Therefore, we cannot exclude symmetry breaking of a X
similar size as for the axial charges, in particular, in the ratio δJSUð3Þ ¼ aJn xn : ð45Þ
of the Ξ over the Σ baryon charge. Finally, in Fig. 32 we n>0
The data for δASUð3Þ ðxÞ depicted in Fig. 33 become more and
more positive as the physical point (vertical dashed line) is
approached. This observation agrees with findings from
earlier studies [22–26]. We fit to data for which the average
quark mass is kept constant (blue circles). However, there is
no significant difference between these and the ms ≈ const
points (black squares). Both linear and quadratic fits in x
(aAn ¼ 0 for n ≠ 1 and aAn ¼ 0 for n ≠ 2, respectively) give
FIG. 31. The same as Fig. 30 for the scalar channel. The yellow
bands depict the ms ¼ ml ¼ 0 predictions DS ð0Þ=FS ð0Þ ob-
tained from the extrapolations of the individual charges.
034512-27
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
adequate descriptions of the data and agree with our C. The up and down quark mass difference
continuum, infinite volume limit extrapolated physical Our results on the scalar charges, in particular, gΣS , enable
point result (red diamond) us to determine the quark mass splitting δm ¼ mu − md .
While we simulate the isosymmetric theory, in Nature
ð23Þ this symmetry is broken. The extent of isospin symmetry
δASUð3Þ ¼ 0.075ð27Þ ; ð46Þ breaking is determined by two small parameters, δm =ΛQCD
and the fine structure constant αQED , which are similar in
derived from the values for the individual charges. Effects size. The vector Ward identity relates δm to the QCD
of this sign and magnitude were also reported previously. contributions to baryon mass splittings within an isomul-
ETMC [24] find gNA þ gΞA − gΣA ¼ 0.147ð24Þ, whereas tiplet. In particular, to leading order in δm =ΛQCD and αQED ,
Savanur and Lin [25] quote ðgNA þgΞA −gΣA Þ=gNA ¼ 0.087ð15Þ. the difference between the Σþ and Σ− baryon masses is a
For J ≠ A no statistically significant effects were pure QCD effect from which, with our knowledge of gΣS , we
observed. Nevertheless, for completeness we carry out can extract δm without additional assumptions.
the same analysis for J ¼ S and J ¼ T, see Fig. 34. Our We consider isospin multiplets of baryons BQ ∈
continuum, infinite volume limit extrapolated physical fN Q ; ΣQ ; ΞQ g with electric charges Q ¼ I 3 þ 12 ð1 þ SÞ ∈
point results f0; 1g (N þ ¼ p, N 0 ¼ n) and define the mass differences
ΔmBQþ1 ¼ mBQþ1 − mBQ . Note that for the Σ there are two
differences [75],
ð37Þ ð16Þ
δSSUð3Þ ¼ −0.040ð41Þ ; δTSUð3Þ ¼ −0.001ð23Þ ð47Þ
ΔmΣþ ¼ mΣþ − mΣ0 ¼ −3.27ð7Þ MeV; ð48Þ
provide upper limits on the relative size of SU(3) flavor
violation at the physical point. ΔmΣ0 ¼ mΣ0 − mΣ− ¼ −4.81ð4Þ MeV: ð49Þ
The mass differences can be split into QCD (∼δm ) and QED
(∼αQED ΛQCD ) contributions:
ΔmB ¼ ΔmQCD
B þ ΔmQED
B : ð51Þ
1
ΔmQCD
Σ ¼ ðmΣþ − mΣ− Þ
2
¼ −4.04ð4Þ MeV; ð52Þ
1
ΔmQED
Σþ
¼ ðmΣþ þ mΣ− Þ − mΣ0
2
FIG. 34. The same as Fig. 33 for the scalar (top) and tensor
(bottom) charges. ¼ 0.77ð5Þ MeV ¼ −ΔmQED
Σ0
: ð53Þ
034512-28
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
From the Ademollo-Gatto theorem [74] we know that the TABLE XIV. Comparison of the light quark mass difference.
leading isospin breaking effects on the vector charges
gNV ¼ gΞV ¼ 1 and gΣV ¼ 2 are quadratic functions of Nf δm /MeV
δm =ΛQCD and αQED , whereas the scalar charges gBS are RM123 [121] 2þ1þ1 −2.38ð18Þ
subject to linear corrections in αQED and δm =ΛQCD . FNAL-MILCa [122] 2þ1þ1 ð9Þ
−2.55ð7Þ
The Lorentz decomposition of the on-shell QCD matrix MILCa [123] 2þ1þ1 ð11Þ
−2.57ð6Þ
element for the isovector vector current between baryons BMWc [108] 2þ1 −2.41ð12Þ
B0 ¼ BQþ1 and B ¼ BQ (that differ by ΔI 3 ¼ 1 in their This work 2þ1 −2.03ð12Þ
isospin) gives [see Eq. (1)] a
FNAL-MILC [122] and MILC [123] only quote the ratios
ð131Þ ð157Þ
B0 B
i∂μ hB0 ðp0 Þjd̄γ μ ujBðpÞi ¼ gV i∂μ ūB0 ðp0 Þγ μ uB ðpÞ mu =md ¼ 0.4556ð93Þ and mu =md ¼ 0.4529ð82Þ , respectively
0
(all errors added in quadrature). Using the FLAG 21 [78]
¼ gVB B ΔmQCD
B ½1 þ Oðδm =ΛQCD Þ; average ml ¼ 3.410ð43Þ, we combine these results to form
δm ¼ 2ml ðmu =md − 1Þ=ðmu =md þ 1Þ and compute the error
ð54Þ by error propagation.
034512-29
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
ΔmQED
N ¼ ΔmN − gNS ðmu − md Þ by BMWc [108] from simulations of QCD with quenched
QED, differs by 2.2 standard deviations from our result in
2gNS
¼ ΔmN − ΔmQCD
Σ ; ð58Þ Eq. (57). Also other lattice results on the QCD contribution
gΣS to the mass-splittings (summarized in Fig. 35), obtained at a
single lattice spacing from Endres et al. [125] using QEDTL
2gΞS and QEDM , Brantley et al. [126] and CSSM-QCDSF-
ΔmQED
Ξ ¼ ΔmΞ − ΔmQCD
Σ : ð59Þ
gΣS UKQCD [127] agree within errors.
We mention the possibility of an enhancement of the
This gives (higher order) δ2m =ΛQCD correction to the Σ0 mass due to the
possibility of mixing with the Λ0 , which, however, appears
ð31Þ
ΔmQED
N ¼ 0.97ð36Þ MeV; to be a very small effect [128]. A positive contribution to the
ð31Þ Σ0 mass would increase ΔmQED Σ but leave ΔmQCD Σ [and
ΔmQCD
N ¼ −2.26ð36Þ MeV; ð60Þ
therefore the quark mass difference Eq. (57)] invariant.
The electromagnetic contributions to the p, Σ and Ξ−
ΔmQED
Σ ¼ 0.77ð05Þ MeV; masses are all similar to 1 MeV, with an enhancement for
ΔmQCD ¼ −4.04ð04Þ MeV; ð61Þ the heavier, more compact cascade baryon. Recently,
Σ
combining the Cottingham formula [118] with experimen-
ð37Þ tal input from elastic scattering and parton distribution
ΔmQED
Ξ ¼ −1.65ð39Þ MeV;
functions, the value ΔmQED N ¼ 0.58ð16Þ MeV was deter-
ð42Þ
ΔmQCD
Ξ ¼ −5.20ð44Þ MeV: ð62Þ mined in Ref. [124]. While within errors our result Eq. (60)
agrees with this value, the number obtained in Ref. [124] is
For completeness we included the values for the Σ baryons more inline with the suggested ordering ΔmQED N <
that we determined from the experimental masses alone, ΔmQED
Σ < −Δm QED
Ξ . Combining their value with our deter-
without lattice input. The above mass splittings agree with ð28Þ
mination of gNS gives δm ¼ −1.69ð26Þ MeV, somewhat
the BMWc [109] continuum limit results from simulations
smaller in modulus than our result Eq. (57) and certainly
of QCD plus QED, see Fig. 35 (errors added in quadrature).
in tension with, e.g., δm ¼ −2.41ð12Þ MeV [108].
Nevertheless, as mentioned above, the value of δm , reported
We find that the effect of md > mu on the Ξ and Σ mass
splittings is much bigger than for the nucleon since this is
proportional to gBS =gBV and gNS < gΣS =2 < gΞS . This hierarchy
is due to gNS ≈ FS þ DS , gΣS =2 ≈ FS and gΞS ≈ FS − DS with
FS > 0 and DS < 0. Interestingly, the pion baryon σ terms
σ πB ¼ σ uB þ σ dB that encode the up plus down quark mass
contribution to the baryon masses exhibit the opposite
ordering [17], σ πN > σ πΣ > σ πΞ .
δm
mu ūu þ md d̄d ¼ ml ðūu þ d̄dÞ þ ðūu − d̄dÞ; ð63Þ
2
or one can start from the Feynman-Hellmann theorem
FIG. 35. Comparison of the QCD contributions to isospin mass hBjq̄qjBi ∂m
splittings [75,109,124–127]. Note that in our normalization σ qB ¼ mq ¼ mq B : ð64Þ
mΣþ − mΣ− ≈ 2ΔmQCD . Values with filled symbols were obtained
hBjBi ∂mq
Σ
via a quark mass, continuum and finite volume extrapolation.
Endres et al. [125] only quote values for ΔmQED from which we Writing mp ¼ mN þ ΔmQCD N =2 þ ΔmN
QED
and mn ¼
N
compute ΔmN QCD
employing the experimental proton-neutron mN − ΔmQCD N =2, where Δm QCD
N ¼ δ g
m S
N N
=g V , and realizing
mass splitting. that the dependence of the QED contributions on the quark
034512-30
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
034512-31
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
research and innovation programme under the Marie positions of the point-to-all propagators. To form the three-
Skłodowska-Curie Grant Agreement No. 813942 (ITN point functions for the nucleon an additional N sto ¼ 100
EuroPLEx) and Grant Agreement No. 824093 light stochastic solves (for the time slice-to-all propagators
(STRONG 2020) is gratefully acknowledged, as well as connecting the sink and current time slices, the wiggly line
initial stage funding through the German Research in Fig. 6) are performed. This setup provides 8 measure-
Foundation (DFG) collaborative research centre SFB/ ments of the nucleon three-point function (as shown in
TRR-55. The authors gratefully acknowledge the Gauss Fig. 7, with the source-sink separations t=a ¼ 11, 14, 16,
Centre for Supercomputing (GCS) for providing computing 19) and includes all polarizations (and the unpolarized
time through the John von Neumann Institute for case) as well as a range of sink momenta (almost) for free.
Computing (NIC) on the supercomputer JUWELS [132] In principle, decuplet baryon three-point functions can also
and in particular on the Booster partition of the supercom- be constructed at the analysis stage. This setup is evaluated
puter JURECA [133] at Jülich Supercomputing Centre twice on each configuration leading to a total of 296
(JSC). GCS is the alliance of the three national super- inversions. Similarly, an additional ð4 × 12 þ 100Þ × 2
computing centres HLRS (Universität Stuttgart), JSC strange solves are performed in order to form the three-
(Forschungszentrum Jülich), and LRZ (Bayerische point functions for all the (octet and decuplet) hyperons,
Akademie der Wissenschaften), funded by the BMBF including strangeness changing currents that we did not
and the German State Ministries for Research of Baden- consider here.
Württemberg (MWK), Bayern (StMWFK) and Nordrhein- In the sequential source setup, we compute the three-
Westfalen (MIWF). Additional simulations were carried point function for the nucleon at rest, again for source-sink
out on the QPACE 3 Xeon Phi cluster of SFB/TRR-55 and separations t=a ¼ 11, 14, 16 and 19. Ten measurements are
the Regensburg Athene 2 Cluster. The authors also thank carried out per configuration (corresponding to 1, 2, 3 and 4
the JSC for their support and for providing services and measurements for each t, respectively), where in each case
computing time on the HDF Cloud cluster [134] at JSC, the two light quark flavors (u and d) of the current and the
funded via the Helmholtz Data Federation (HDF) pro- four possible polarizations of the nucleon require 2 × 4
gramme. Most of the ensembles were generated using sequential sources to be constructed. This amounts to
openQCD [35] within the CLS effort. A few additional performing ð4 þ 10 × 2 × 4Þ × 12 ¼ 1008 light solves.
ensembles were generated employing the BQCD-code [39] The additional 4 × 12 inversions refer to the point-to-all
on the QPACE supercomputer of SFB/TRR-55. For the propagators for 4 different source positions that connect the
computation of hadronic two- and three-point functions we source to the sink (and the current). This is three-times the
used a modified version of the Chroma [135] software cost of the stochastic approach (for the nucleon three-point
package along with the LibHadronAnalysis library and the functions), which realizes a range of sink momenta.
multigrid solver implementation of Refs. [136,137] (see The ratios of the three-point over two-point functions for
also Ref. [138]) as well as the IDFLS solver [139] of the nucleon obtained from the two different approaches are
openQCD. We used Matplotlib [140] to create the figures. compared in Fig. 36. A significant part of the gauge noise
cancels in the ratio, while the (additional) stochastic noise
remains. For our setup, this leads to larger statistical errors
APPENDIX A: FURTHER DETAILS OF THE
for the stochastic data compared to the sequential source
THREE-POINT FUNCTION MEASUREMENTS
results. This difference can clearly be seen for the ratio in
1. Comparison of the stochastic the vector channel, for which the gauge noise is minimal,
and sequential source methods however, the difference is less pronounced for the other
We computed the connected three-point functions for all charges. For the sigma and cascade baryons we generally
the octet baryons utilizing the computationally efficient find a good statistical signal in the ratios employing the
stochastic approach outlined in Sec. III C. This approach stochastic approach, see Fig. 37, although, also in this case
introduces additional stochastic noise on top of the gauge the ratios for the hyperon vector charges suffer from large
noise. In the analysis presented, for the nucleon, we make errors.
use of statistically more precise three-point correlation In the case of a large-scale analysis effort including high
function measurements determined via the sequential statistics, as presented in this article, the disk space required
source method as part of other projects. In the following, to store the stochastic three-point function data is signifi-
we compare the computational costs of the stochastic and cant. The individual spectator (S) and insertion (I) parts as
the sequential source methods and the results for the ratios defined in Eqs. (22) and (23), respectively, are stored
of the three-point over two-point functions for the nucleon. with all indices open. In general this amounts to N ¼
As a typical example, we consider the measurements N½S þ N½I complex double precision floating point num-
performed on ensemble N200 (M π ¼ 286 MeV and bers for each gauge field configuration where
a ¼ 0.064 fm). For our setup, illustrated in Figs. 6 and
7, a total of 4 × 12 solves are needed for the 4 source N½S ¼ N SF · N p0 · N snk · N src · N sto · N c · N 5s ; ðA1Þ
034512-32
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
FIG. 36. Unrenormalized ratio of three-point over two-point functions for the nucleon vector, axial, scalar and tensor charge (from top
to bottom) on ensemble N200. The left hand side shows the data from the sequential source method from 1, 2, 3 and 4 measurements (for
increasing values of t) compared to the same measurements obtained from the stochastic approach with four measurements for each
value of t on the right-hand side.
034512-33
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
FIG. 37. Unrenormalized ratio of three-point over two-point functions for the Σ (left) and Ξ (right) scalar charge on ensemble N200.
Here N S=I
F denotes the number of flavor combinations for
the spectator and insertion parts (typically 4 and 2,
respectively), N p gives the number of momentum combi-
nations for a maximum momentum jpj (with p either being
the sink momentum p0 or the momentum transfer q),
N src=snk corresponds to the number of source (4) and sink
(2) positions, N c ¼ 3 and N s ¼ 4 are the dimensions of
color and spin space and N τ is the number of current
insertion time slices, usually the distance between the two
sink time slices, see Figs. 6 and 7. N ∂μ refers to the number
of derivatives included in the current insertion. We consider
all currents including up to one derivative (N ∂μ ¼ 1 þ 4),
FIG. 38. Nucleon three-point function for the scalar channel
although only the currents without derivatives are presented
with a source-sink separation t=a ¼ 16 for a single source position
in this work. This adds up to a file size of the order of GBs on ensemble D200. The average and standard error is computed
for a single gauge field configuration and disk space usage including and excluding one particular gauge configuration.
of the order of TBs for a typical CLS gauge ensemble.
Storing the data with all indices open allows for a very
flexible analysis. Octet or decuplet baryon three-point
functions can be constructed from the spectator and
insertion parts for different polarizations, current insertions
as well as for a large number of momentum combinations.
2. Treatment of outliers
When analysing the three-point functions on some of the
ensembles, we observe a small number of three-point
function results that are by many orders of magnitude
larger than the rest. These outliers, whose origin currently
remains unexplained, would have a significant impact on
the analysis, as illustrated in Fig. 38 for ensemble D200
(M π ¼ 202 MeV and a ¼ 0.064 fm). The three-point func-
tion for the scalar charge with source-sink separation t=a ¼
16 for a single source position is displayed. In this case one FIG. 39. The numbers of configurations separated by a distance
outlier is identified (according to the criterion given below) jOwi − Ôj=ΔO from the central value Ô for the three-point
and one sees a substantial change in the configuration function data of Fig. 38 on the insertion time slice τ=a ¼ 15.
average and a reduction in the standard deviation if this The blue dashed line indicates the cutoff K ¼ 100, which is
measurement (on a particular configuration) is excluded. exceeded by a single measurement. Note the cut and different
The Mainz group [141] reported similar outliers when scales of the ordinate. Due to this representation a few histogram
determining the nucleon σ terms on a subset of the CLS entries are not shown.
034512-34
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
TABLE XV. The gauge ensembles utilized in this work: the light and strange hopping parameters κ l and κ s , the gradient flow scale
parameter t0 =a2 and the pion (M π ) and kaon masses (M K ) [17].
Ensemble κl κs t0 =a2 aM π aM K
A653 0.1365716 0.1365716 2.1729(50) 0.21235(94) 0.21235(94)
A650 0.136600 0.136600 2.2878(72) 0.1833(13) 0.1833(13)
A654 0.136750 0.136216193 2.1950(77) 0.1669(11) 0.22714(91)
H101 0.13675962 0.13675962 2.8545(81) 0.18283(57) 0.18283(57)
U103 0.13675962 0.13675962 2.8815(57) 0.18133(61) 0.18133(61)
H107 0.13694566590798 0.136203165143476 2.7193(76) 0.15913(73) 0.23745(53)
H102r002 0.136865 0.136549339 2.8792(90) 0.15490(92) 0.19193(77)
U102 0.136865 0.136549339 2.8932(63) 0.15444(84) 0.19235(61)
H102r001 0.136865 0.136549339 2.8840(89) 0.15311(98) 0.19089(78)
rqcd021 0.136813 0.136813 3.032(15) 0.14694(88) 0.14694(88)
H105 0.136970 0.13634079 2.8917(65) 0.1213(14) 0.20233(64)
N101 0.136970 0.13634079 2.8948(39) 0.12132(58) 0.20156(30)
H106 0.137015570024 0.136148704478 2.8227(68) 0.1180(21) 0.22471(67)
C102 0.13705084580022 0.13612906255557 2.8682(47) 0.09644(77) 0.21783(36)
C101 0.137030 0.136222041 2.9176(38) 0.09586(64) 0.20561(33)
D101 0.137030 0.136222041 2.910(10) 0.0958(11) 0.20572(45)
S100 0.137030 0.136222041 2.9212(91) 0.0924(31) 0.20551(57)
D150 0.137088 0.13610755 2.9476(30) 0.05497(79) 0.20834(17)
B450 0.136890 0.136890 3.663(11) 0.16095(49) 0.16095(49)
S400 0.136984 0.136702387 3.6919(74) 0.13535(42) 0.17031(38)
B452 0.1370455 0.136378044 3.5286(66) 0.13471(47) 0.20972(34)
rqcd030 0.1369587 0.1369587 3.914(15) 0.12202(68) 0.12202(68)
N451 0.1370616 0.1365480771 3.6822(46) 0.11067(32) 0.17828(20)
N401 0.1370616 0.1365480771 3.6844(52) 0.10984(57) 0.17759(37)
N450 0.1370986 0.136352601 3.5920(42) 0.10965(31) 0.20176(18)
X450 0.136994 0.136994 3.9935(92) 0.10142(62) 0.10142(62)
D451 0.137140 0.136337761 3.6684(36) 0.08370(31) 0.19385(15)
D450 0.137126 0.136420428639937 3.7076(75) 0.08255(41) 0.18354(12)
D452 0.137163675 0.136345904546 3.7251(37) 0.05961(50) 0.18647(13)
(Table continued)
034512-35
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
Ensemble κl κs t0 =a2 aM π aM K
N202 0.137000 0.137000 5.165(14) 0.13388(35) 0.13388(35)
N204 0.137112 0.136575049 4.9473(79) 0.11423(33) 0.17734(29)
X250 0.137050 0.137050 5.283(28) 0.11319(39) 0.11319(39)
N203 0.137080 0.136840284 5.1465(63) 0.11245(30) 0.14399(24)
S201 0.137140 0.13672086 5.1638(91) 0.09379(47) 0.15220(37)
N201 0.13715968 0.136561319 5.0427(75) 0.09268(31) 0.17040(22)
N200 0.137140 0.13672086 5.1600(71) 0.09236(29) 0.15061(24)
X251 0.137100 0.137100 5.483(26) 0.08678(40) 0.08678(40)
D200 0.137200 0.136601748 5.1793(39) 0.06540(33) 0.15652(15)
D201 0.1372067 0.136546844 5.1378(66) 0.06472(42) 0.16302(18)
E250 0.137232867 0.136536633 5.2027(41) 0.04227(23) 0.159370(61)
N300 0.137000 0.137000 8.576(21) 0.10642(38) 0.10642(38)
N304 0.137079325093654 0.136665430105663 8.322(20) 0.08840(33) 0.13960(31)
N302 0.137064 0.1368721791358 8.539(19) 0.08701(41) 0.11370(36)
J304 0.137130 0.1366569203 8.497(12) 0.06538(18) 0.13181(14)
J303 0.137123 0.1367546608 8.615(14) 0.06481(19) 0.11975(16)
E300 0.137163 0.1366751636177327 8.6241(74) 0.04402(20) 0.12397(15)
J500 0.136852 0.136852 14.013(34) 0.08116(34) 0.08116(34)
J501 0.1369032 0.136749715 13.928(39) 0.06589(26) 0.08798(23)
1 2
#ES
Ensemble gN;latt
A gN;latt
S gN;latt
T gN;latt
V gN;latt
A gN;latt
S gN;latt
T gN;latt
V
034512-36
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
1 2
#ES
Ensemble gN;latt
A gN;latt
S gN;latt
T gN;latt
V gN;latt
A gN;latt
S gN;latt
T gN;latt
V
TABLE XVII. Results for the unrenormalized hyperon charges gJΣ;latt for J ∈ fA; S; T; Vg. #ES labels the number of excited states used
to determine the ground state matrix element, see the discussion in Sec. III D. The three-point functions are computed employing the
“stochastic” approach.
1 2
#ES
Ensemble gAΣ;latt gΣ;latt
S gTΣ;latt gVΣ;latt gAΣ;latt gΣ;latt
S gΣ;latt
T gVΣ;latt
A653 1.180(10) 4.40(14) 1.018(14) 2.8354(32) 1.189(16) 4.45(14) 1.023(12) 2.8347(27)
A650 1.1764(90) 4.46(12) 0.9968(95) 2.8381(39) 1.190(16) 4.49(15) 1.000(10) 2.8380(37)
A654 1.173(12) 4.01(16) 1.028(13) 2.852(12) 1.204(37) 4.11(30) 1.028(16) 2.850(14)
H101 1.1877(81) 5.03(12) 0.982(10) 2.78711(66) 1.198(13) 5.12(14) 0.987(10) 2.78686(65)
U103 1.137(13) 4.69(18) 0.983(16) 2.78535(89) 1.147(21) 4.75(18) 0.989(16) 2.78507(80)
H107 1.207(18) 4.25(19) 1.030(12) 2.823(15) 1.232(28) 4.35(27) 1.036(13) 2.820(11)
H102r002 1.177(13) 4.67(17) 0.988(13) 2.807(14) 1.199(25) 4.77(22) 0.992(14) 2.805(12)
H102r001 1.157(14) 4.73(19) 0.973(16) 2.788(11) 1.146(27) 4.79(28) 0.978(17) 2.788(11)
rqcd021 1.149(15) 5.27(24) 0.958(19) 2.790(14) 1.145(27) 5.35(33) 0.962(19) 2.790(13)
H105 1.147(22) 4.83(44) 0.958(23) 2.812(20) 1.162(40) 5.25(55) 0.969(21) 2.809(16)
H106 1.139(15) 4.24(20) 0.968(16) 2.796(15) 1.149(38) 4.27(48) 0.970(19) 2.796(15)
C102 1.140(19) 4.32(29) 0.950(27) 2.791(23) 1.133(45) 3.94(76) 0.952(22) 2.792(18)
C101 1.174(21) 4.58(43) 0.957(22) 2.837(21) 1.240(51) 4.80(94) 0.963(21) 2.834(17)
B450 1.1743(76) 4.80(14) 0.992(11) 2.7464(30) 1.195(15) 4.90(18) 0.997(11) 2.7460(28)
S400 1.159(16) 5.25(29) 0.945(17) 2.751(13) 1.180(23) 5.43(27) 0.962(16) 2.752(11)
B452 1.147(13) 4.33(18) 0.990(16) 2.764(13) 1.163(28) 4.45(33) 0.994(16) 2.763(12)
rqcd030 1.140(12) 5.20(22) 0.944(16) 2.7576(35) 1.157(22) 5.41(33) 0.950(15) 2.7573(33)
N451 1.168(15) 4.90(34) 0.976(14) 2.795(18) 1.193(30) 5.43(38) 0.988(11) 2.789(12)
N401 1.130(52) 5.22(68) 0.933(50) 2.781(42) 1.100(67) 5.23(82) 0.949(33) 2.776(31)
(Table continued)
034512-37
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
1 2
#ES
Ensemble gAΣ;latt gΣ;latt
S gTΣ;latt gVΣ;latt gAΣ;latt gΣ;latt
S gΣ;latt
T gVΣ;latt
N450 1.148(18) 4.11(17) 0.994(20) 2.738(24) 1.154(43) 3.54(65) 0.992(17) 2.737(21)
X450 1.168(23) 6.05(38) 0.924(25) 2.7700(84) 1.195(49) 6.38(58) 0.933(25) 2.7690(79)
D451 1.143(53) 4.79(95) 0.955(38) 2.67(10) 1.18(15) 3.8(3.0) 0.948(45) 2.624(56)
D450 1.129(63) 4.9(1.1) 0.938(52) 2.721(77) 1.09(11) 3.9(2.1) 0.944(42) 2.730(48)
D452 1.115(52) 4.8(1.1) 0.923(52) 2.789(28) 1.07(14) 5.3(2.6) 0.935(44) 2.787(29)
N202 1.166(10) 5.68(20) 0.939(12) 2.69777(45) 1.182(16) 5.85(15) 0.955(13) 2.69721(41)
N204 1.1465(91) 4.30(12) 0.980(11) 2.686(11) 1.160(26) 4.22(30) 0.979(12) 2.686(12)
X250 1.1549(88) 5.85(20) 0.928(12) 2.7051(32) 1.167(16) 6.12(21) 0.940(12) 2.7046(27)
N203 1.156(12) 5.60(26) 0.941(10) 2.7246(96) 1.182(19) 6.05(22) 0.9623(87) 2.7184(69)
N201 1.150(17) 5.22(43) 0.953(18) 2.737(17) 1.163(42) 5.70(53) 0.965(17) 2.733(14)
N200 1.119(13) 5.03(21) 0.947(13) 2.721(11) 1.116(26) 5.37(36) 0.954(12) 2.720(10)
X251 1.145(16) 6.59(33) 0.903(13) 2.7165(69) 1.195(24) 7.06(37) 0.916(14) 2.7145(55)
D200 1.082(48) 6.4(1.4) 0.841(59) 2.774(26) 1.058(82) 7.7(1.8) 0.867(38) 2.777(29)
D201 1.125(34) 5.7(1.3) 0.924(37) 2.730(31) 1.133(62) 6.8(1.2) 0.944(24) 2.726(19)
N300 1.121(12) 5.52(19) 0.918(17) 2.63208(25) 1.126(19) 5.59(18) 0.925(16) 2.63197(23)
N304 1.106(18) 4.84(31) 0.925(21) 2.642(17) 1.130(30) 4.87(45) 0.928(20) 2.641(15)
N302 1.086(26) 5.41(35) 0.864(34) 2.668(25) 1.118(42) 5.38(49) 0.868(26) 2.666(21)
J304 1.091(29) 5.95(65) 0.864(30) 2.646(25) 1.149(54) 6.61(74) 0.879(23) 2.645(18)
J303 1.100(19) 5.23(59) 0.914(22) 2.671(31) 1.198(61) 5.88(96) 0.927(17) 2.669(26)
E300 1.047(93) 7.9(3.8) 0.853(57) 2.728(87) 1.04(12) 9.1(3.1) 0.884(41) 2.705(55)
J500 1.100(10) 5.78(20) 0.887(13) 2.58250(15) 1.106(15) 5.82(16) 0.894(11) 2.58241(12)
J501 1.053(35) 6.01(86) 0.857(46) 2.629(22) 1.038(54) 6.26(66) 0.878(22) 2.628(19)
TABLE XVIII. Results for the unrenormalized hyperon charges gΞ;latt J for J ∈ fA; S; T; Vg. #ES labels the number of excited states
used to determine the ground state matrix element, see the discussion in Sec. III D. The three-point functions are computed employing
the “stochastic” approach.
1 2
#ES
Ensemble gΞ;latt
A gΞ;latt
S gΞ;latt
T gΞ;latt
V gΞ;latt
A gΞ;latt
S gΞ;latt
T gΞ;latt
V
034512-38
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
1 2
#ES
Ensemble gΞ;latt
A gΞ;latt
S gΞ;latt
T gΞ;latt
V gΞ;latt
A gΞ;latt
S gΞ;latt
T gΞ;latt
V
034512-39
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
FIG. 41. The same as Fig. 13 for the isovector vector charges gBV of the sigma baryon (left) and the cascade baryon (right). For better
visibility, the data points for ensemble D451, which have large errors (see Tables XVII and XVIII), are not displayed. See Tables XVII
and XVIII for the set of ensembles used.
FIG. 42. The same as Fig. 15 for the vector charge gBV of the sigma (right) and cascade (left) baryon.
FIG. 43. The same as Fig. 19 for the axial charge gBA of the sigma (right) and cascade (left) baryon.
034512-40
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
FIG. 44. The same as Fig. 19 for the scalar charge gBS of the nucleon (top, left), sigma (top, right) and cascade (bottom) baryons. The
six fits correspond to two fit variations, see the text, applied to three datasets, DSðM π<400 MeV Þ, DSðM π<400 MeV ; a<0.1 fm Þ and
DSðM <400 π Þ. The data are extracted using two excited states in the fitting analysis, see Sec. III D. Only two different sets of
MeV
π ; LM >4
renormalization factors are used.
FIG. 45. The same as Fig. 19 for the tensor charge gBT for the nucleon (top, left), sigma (top, right) and cascade (bottom) baryons. The
data are extracted using two excited states in the fitting analysis, see Sec. III D. Only two different sets of renormalization
factors are used.
034512-41
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
FIG. 46. The same as Fig. 21 for the isovector scalar charges gBS of the sigma baryon (left) and the cascade baryon (right). For better
visibility, the data point for ensemble E300, which has a relatively large error (see Table XVII), is not displayed for the sigma baryon.
See Tables XVII and XVIII for the set of ensembles used.
FIG. 47. The same as Fig. 23 for the isovector tensor charges gBT of the sigma baryon (left) and the cascade baryon (right). For better
visibility, the data point for ensemble E300, which has a relatively large error (see Table XVIII), is not displayed for the cascade baryon.
See Tables XVII and XVIII for the set of ensembles used.
034512-42
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
[1] X.-D. Ji, Gauge-Invariant Decomposition of Nucleon Spin, (ultra)cold neutrons to the LHC, Phys. Rev. D 85, 054512
Phys. Rev. Lett. 78, 610 (1997). (2012).
[2] R. L. Jaffe and A. Manohar, The g1 problem: Fact and [17] G. S. Bali, S. Collins, P. Georg, D. Jenkins, P. Korcyl, A.
fantasy on the spin of the proton, Nucl. Phys. B337, 509 Schäfer, E. E. Scholz, J. Simeth, W. Söldner, and S.
(1990). Weishäupl (RQCD Collaboration), Scale setting and the
[3] G. S. Bali et al. (QCDSF Collaboration), Strangeness light baryon spectrum in N f ¼ 2 þ 1 QCD with Wilson
Contribution to the Proton Spin from Lattice QCD, Phys. fermions, J. High Energy Phys. 05 (2023) 035.
Rev. Lett. 108, 222001 (2012). [18] H. Leutwyler, Theoretical aspects of chiral dynamics, Proc.
[4] M. Engelhardt, Strange quark contributions to nucleon Sci. CD15 (2015) 022 [arXiv:1510.07511].
mass and spin from lattice QCD, Phys. Rev. D 86, 114510 [19] F. Fucito, G. Parisi, and S. Petrarca, First evaluation of
(2012). gA =gV in lattice QCD in the quenched approximation,
[5] M. Gong, Y.-B. Yang, J. Liang, A. Alexandru, T. Draper, Phys. Lett. 115B, 148 (1982).
and K.-F. Liu (χQCD Collaboration), Strange and charm [20] S. Güsken, U. Löw, K.-H. Mütter, R. Sommer, A. Patel,
quark spins from the anomalous Ward identity, Phys. Rev. and K. Schilling, Nonsinglet axial vector couplings of the
D 95, 114509 (2017). baryon octet in lattice QCD, Phys. Lett. B 227, 266 (1989).
[6] C. Alexandrou, M. Constantinou, K. Hadjiyiannakou, K. [21] R. M. Woloshyn and K.-F. Liu, A study of the nucleon
Jansen, C. Kallidonis, G. Koutsou, A. Vaquero Avilés- axial current coupling in quenched lattice QCD, Nucl.
Casco, and C. Wiese (ETM Collaboration), Nucleon Spin Phys. B311, 527 (1989).
and Momentum Decomposition Using Lattice QCD Sim- [22] H.-W. Lin and K. Orginos, First calculation of hyperon
ulations, Phys. Rev. Lett. 119, 142002 (2017). axial couplings from lattice QCD, Phys. Rev. D 79, 034507
[7] J. Green, N. Hasan, S. Meinel, M. Engelhardt, S. Krieg, J. (2009).
Laeuchli, J. Negele, K. Orginos, A. Pochinsky, and S. [23] G. Erkol, M. Oka, and T. T. Takahashi, Axial charges of
Syritsyn, Up, down, and strange nucleon axial form factors octet baryons in two-flavor lattice QCD, Phys. Lett. B 686,
from lattice QCD, Phys. Rev. D 95, 114502 (2017). 36 (2010).
[8] S. Bürger, T. Wurm, M. Löffler, M. Göckeler, G. Bali, S. [24] C. Alexandrou, K. Hadjiyiannakou, and C. Kallidonis
Collins, A. Schäfer, and A. Sternbeck (RQCD Collabora- (ETM Collaboration), Axial charges of hyperons and
tion), Lattice results for the longitudinal spin structure and charmed baryons using N f ¼ 2 þ 1 þ 1 twisted mass
color forces on quarks in a nucleon, Phys. Rev. D 105, fermions, Phys. Rev. D 94, 034502 (2016).
054504 (2022). [25] A. Savanur and H.-W. Lin, Lattice-QCD determination of
[9] D. de Florian, R. Sassot, M. Stratmann, and W. Vogelsang, the hyperon axial couplings in the continuum limit, Phys.
Global Analysis of Helicity Parton Densities and Their Rev. D 102, 014501 (2020).
Uncertainties, Phys. Rev. Lett. 101, 072001 (2008). [26] R. E. Smail et al. (QCDSF-UKQCD-CSSM Collabora-
[10] J. Blümlein and H. Böttcher, QCD analysis of polarized tion), Constraining beyond the Standard Model nucleon
deep inelastic scattering data, Nucl. Phys. B841, 205 isovector charges, arXiv:2304.02866.
(2010). [27] B. Borasoy and U.-G. Meissner, Chiral expansion of baryon
[11] E. R. Nocera, R. D. Ball, S. Forte, G. Ridolfi, and J. Rojo masses and σ-terms, Ann. Phys. (N.Y.) 254, 192 (1997).
(NNPDF Collaboration), A first unbiased global determi- [28] A. N. Cooke, R. Horsley, Y. Nakamura, D. Pleiter, P. E. L.
nation of polarized PDFs and their uncertainties, Nucl. Rakow, G. Schierholz, and J. M. Zanotti, The effects of
Phys. B887, 276 (2014). flavor symmetry breaking on hadron matrix elements,
[12] N. Sato, W. Melnitchouk, S. E. Kuhn, J. J. Ethier, and A. Proc. Sci. LATTICE2012 (2012) 116 [arXiv:1212.2564].
Accardi (JAM Collaboration), Iterative Monte Carlo analy- [29] J. M. Bickerton, R. Horsley, Y. Nakamura, H. Perlt, D.
sis of spin-dependent parton distributions, Phys. Rev. D 93, Pleiter, P. E. L. Rakow, G. Schierholz, H. Stüben, R. D.
074005 (2016). Young, and J. M. Zanotti (QCDSF-UKQCD-CSSM Col-
[13] J. J. Ethier, N. Sato, and W. Melnitchouk (JAM Collabo- laboration), Patterns of flavor symmetry breaking in
ration), First Simultaneous Extraction of Spin-Dependent hadron matrix elements involving u, d, and s quarks,
Parton Distributions and Fragmentation Functions from a Phys. Rev. D 100, 114516 (2019).
Global QCD Analysis, Phys. Rev. Lett. 119, 132001 [30] E. E. Jenkins and A. V. Manohar, Chiral corrections to the
(2017). baryon axial currents, Phys. Lett. B 259, 353 (1991).
[14] H.-W. Lin, W. Melnitchouk, A. Prokudin, N. Sato, and H. [31] R. F. Dashen and A. V. Manohar, 1=N c corrections to the
Shows (JAM Collaboration), First Monte Carlo Global baryon axial currents in QCD, Phys. Lett. B 315, 438
Analysis of Nucleon Transversity with Lattice QCD (1993).
Constraints, Phys. Rev. Lett. 120, 152502 (2018). [32] M. Bruno et al. (CLS Effort Collaboration), Simulation of
[15] L. Gamberg, M. Malda, J. A. Miller, D. Pitonyak, A. QCD with N f ¼ 2 þ 1 flavors of nonperturbatively im-
Prokudin, and N. Sato (JAM Collaboration), Updated proved Wilson fermions, J. High Energy Phys. 02 (2015)
QCD global analysis of single transverse-spin asymme- 043.
tries: Extracting H̃, and the role of the Soffer bound and [33] M. Lüscher and S. Schaefer, Lattice QCD without top-
lattice QCD, Phys. Rev. D 106, 034014 (2022). ology barriers, J. High Energy Phys. 07 (2011) 036.
[16] T. Bhattacharya, V. Cirigliano, S. D. Cohen, A. Filipuzzi, [34] S. Schaefer, R. Sommer, and F. Virotta (ALPHA Collabo-
M. Gonzalez-Alonso, M. L. Graesser, R. Gupta, and H.-W. ration), Critical slowing down and error analysis in lattice
Lin, Probing novel scalar and tensor interactions from QCD simulations, Nucl. Phys. B845, 93 (2011).
034512-43
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
[35] M. Lüscher and S. Schaefer, Lattice QCD with open [51] S.-J. Dong and K.-F. Liu, Stochastic estimation with Z2
boundary conditions and twisted-mass reweighting, noise, Phys. Lett. B 328, 130 (1994).
Comput. Phys. Commun. 184, 519 (2013). [52] S. Bernardson, P. McCarty, and C. Thron, Monte Carlo
[36] M. A. Clark and A. D. Kennedy, Accelerating Dynamical methods for estimating linear combinations of inverse
Fermion Computations Using the Rational Hybrid matrix entries in lattice QCD, Comput. Phys. Commun.
Monte Carlo (RHMC) Algorithm with Multiple Pseudo- 78, 256 (1993).
fermion Fields, Phys. Rev. Lett. 98, 051601 (2007). [53] T. Harris, G. von Hippel, P. Junnarkar, H. B. Meyer, K.
[37] D. Mohler and S. Schaefer, Remarks on strange-quark Ottnad, J. Wilhelm, H. Wittig, and L. Wrang, Nucleon
simulations with Wilson fermions, Phys. Rev. D 102, isovector charges and twist-2 matrix elements with
074506 (2020). N f ¼ 2 þ 1 dynamical Wilson quarks, Phys. Rev. D
[38] G. S. Bali, E. E. Scholz, J. Simeth, and W. Söldner (RQCD 100, 034513 (2019).
Collaboration), Lattice simulations with N f ¼ 2 þ 1 im- [54] P. Korcyl and G. S. Bali, Non-perturbative determination
proved Wilson fermions at a fixed strange quark mass, of improvement coefficients using coordinate space corre-
Phys. Rev. D 94, 074501 (2016). lators in N f ¼ 2 þ 1 lattice QCD, Proc. Sci. LAT-
[39] Y. Nakamura and H. Stüben, BQCD—Berlin quantum TICE2016 (2016) 190 [arXiv:1609.09477].
chromodynamics program, Proc. Sci. LATTICE2010 [55] G. S. Bali, K. G. Chetyrkin, P. Korcyl, and J. Simeth, Non-
(2010) 040 [arXiv:1011.0199]. perturbative determination of quark-mass independent
[40] M. Falcioni, M. L. Paciello, G. Parisi, and B. Taglienti improvement coefficients in nf ¼ 2 þ 1 lattice QCD (to
(APE Collaboration), Again on SU(3) glueball mass, Nucl. be published).
Phys. B251, 624 (1985). [56] G. S. Bali et al. (RQCD Collaboration), Nonperturbative
[41] L. Maiani, G. Martinelli, M. L. Paciello, and B. Taglienti, renormalization in lattice QCD with three flavors of clover
Scalar densities and baryon mass differences in lattice fermions: Using periodic and open boundary conditions,
QCD with Wilson fermions, Nucl. Phys. B293, 420 Phys. Rev. D 103, 094511 (2021); 107, 039901(E) (2023).
(1987). [57] G. Martinelli, C. Pittori, C. T. Sachrajda, M. Testa, and
[42] J. D. Bratt et al. (LHP Collaboration), Nucleon structure A. Vladikas, A general method for nonperturbative re-
from mixed action calculations using 2 þ 1 flavors of normalization of lattice operators, Nucl. Phys. B445, 81
asqtad sea and domain wall valence fermions, Phys. Rev. D (1995).
82, 094502 (2010). [58] C. Sturm, Y. Aoki, N. H. Christ, T. Izubuchi, C. T. C.
[43] G. S. Bali, L. Barca, S. Collins, M. Gruber, M. Löffler, A. Sachrajda, and A. Soni, Renormalization of quark bilinear
Schäfer, W. Söldner, P. Wein, S. Weishäupl, and T. Wurm operators in a momentum-subtraction scheme with a
(RQCD Collaboration), Nucleon axial structure from nonexceptional subtraction point, Phys. Rev. D 80, 014501
lattice QCD, J. High Energy Phys. 05 (2020) 126. (2009).
[44] R. Evans, G. Bali, and S. Collins, Improved semileptonic [59] B. A. Kniehl and O. L. Veretin, Bilinear quark operators in
form factor calculations in lattice QCD, Phys. Rev. D 82, the RI/SMOM scheme at three loops, Phys. Lett. B 804,
094501 (2010). 135398 (2020).
[45] G. S. Bali, S. Collins, B. Gläßle, M. Göckeler, J. Najjar, R. [60] A. Bednyakov and A. Pikelner, Quark masses: N3LO
Rödl, A. Schäfer, A. Sternbeck, and W. Söldner, Nucleon bridge from RI=SMOM to MS scheme, Phys. Rev. D 101,
structure from stochastic estimators, Proc. Sci. LAT- 091501 (2020).
TICE2013 (2014) 271 [arXiv:1311.1718]. [61] B. A. Kniehl and O. L. Veretin, Moments n ¼ 2 and n ¼ 3
[46] G. S. Bali, S. Collins, P. Korcyl, R. Rödl, S. Weishäupl, of the Wilson twist-two operators at three loops in the
and T. Wurm, Hyperon couplings from N f ¼ 2 þ 1 lattice RI0 =SMOM scheme, Nucl. Phys. B961, 115229 (2020).
QCD, Proc. Sci. LATTICE2019 (2019) 099 [arXiv:1907 [62] S. Sint, The chirally rotated Schrödinger functional with
.13454]. Wilson fermions and automatic O(a) improvement, Nucl.
[47] G. S. Bali, S. Collins, B. Gläßle, S. Heybrock, P. Korcyl, Phys. B847, 491 (2011).
M. Löffler, R. Rödl, and A. Schäfer, Baryonic and mesonic [63] M. Dalla Brida, T. Korzec, S. Sint, and P. Vilaseca, High
3-point functions with open spin indices, EPJ Web Conf. precision renormalization of the flavour nonsinglet
175, 06014 (2018). Noether currents in lattice QCD with Wilson quarks,
[48] M. Löffler, P. Wein, T. Wurm, S. Weishäupl, D. Jenkins, R. Eur. Phys. J. C 79, 23 (2019).
Rödl, A. Schäfer, and L. Walter (RQCD Collaboration), [64] G. S. Bali, S. Collins, W. Söldner, and S. Weishäupl
Mellin moments of spin dependent and independent PDFs (RQCD Collaboration), Leading order mesonic and bar-
of the pion and rho meson, Phys. Rev. D 105, 014505 yonic SU(3) low energy constants from N f ¼ 3 lattice
(2022). QCD, Phys. Rev. D 105, 054516 (2022).
[49] C. Alexandrou, S. Dinter, V. Drach, K. Jansen, K. [65] E. E. Jenkins and A. V. Manohar, Baryon chiral perturba-
Hadjiyiannakou, and D. B. Renner (ETM Collaboration), tion theory using a heavy fermion Lagrangian, Phys.
A stochastic method for computing hadronic matrix Lett. B 255, 558 (1991).
elements, Eur. Phys. J. C 74, 2692 (2014). [66] J. Bijnens, H. Sonoda, and M. B. Wise, On the validity of
[50] Y.-B. Yang, A. Alexandru, T. Draper, M. Gong, and K.-F. chiral perturbation theory for weak hyperon decays, Nucl.
Liu (χQCD Collaboration), Stochastic method with low Phys. B261, 185 (1985).
mode substitution for nucleon isovector matrix elements, [67] T. Ledwig, J. Martin Camalich, L. S. Geng, and M. J.
Phys. Rev. D 93, 034503 (2016). Vicente Vacas, Octet-baryon axial-vector charges and
034512-44
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
SU(3)-breaking effects in the semileptonic hyperon de- N f ¼ 2 þ 1 ensembles, Proc. Sci. LATTICE2021 (2022)
cays, Phys. Rev. D 90, 054502 (2014). 343 [arXiv:2110.10500].
[68] C. C. Chang et al., A per-cent-level determination of the [86] R. Tsuji, N. Tsukamoto, Y. Aoki, K.-I. Ishikawa, Y.
nucleon axial coupling from quantum chromodynamics, Kuramashi, S. Sasaki, E. Shintani, and T. Yamazaki (PACS
Nature (London) 558, 91 (2018). Collaboration), Nucleon isovector couplings in N f ¼ 2þ1
[69] R. Gupta, Y.-C. Jang, B. Yoon, H.-W. Lin, V. Cirigliano, lattice QCD at the physical point, Phys. Rev. D 106,
and T. Bhattacharya (PNDME Collaboration), Isovector 094505 (2022).
charges of the nucleon from 2 þ 1 þ 1-flavor lattice QCD, [87] T. Bhattacharya, V. Cirigliano, S. Cohen, R. Gupta, H.-W.
Phys. Rev. D 98, 034503 (2018). Lin, and B. Yoon (PNDME Collaboration), Axial, scalar
[70] M. F. M. Lutz, U. Sauerwein, and R. G. E. Timmermans, and tensor charges of the nucleon from 2 þ 1 þ 1-flavor
On the axial-vector form factor of the nucleon and chiral lattice QCD, Phys. Rev. D 94, 054508 (2016).
symmetry, Eur. Phys. J. C 80, 844 (2020). [88] E. Berkowitz et al., An accurate calculation of the nucleon
[71] M. Lüscher, Properties and uses of the Wilson flow in axial charge with lattice QCD, arXiv:1704.01114.
lattice QCD, J. High Energy Phys. 08 (2010) 071; 03 [89] C. Alexandrou, S. Bacchio, M. Constantinou, J.
(2014) 92(E). Finkenrath, K. Hadjiyiannakou, K. Jansen, G. Koutsou,
[72] M. Bruno, T. Korzec, and S. Schaefer, Setting the scale for and A. Vaquero Aviles-Casco (ETM Collaboration), Nu-
the CLS 2 þ 1 flavor ensembles, Phys. Rev. D 95, 074504 cleon axial, tensor, and scalar charges and σ-terms in lattice
(2017). QCD, Phys. Rev. D 102, 054517 (2020).
[73] S. Aoki et al., Review of lattice results concerning low- [90] A. Walker-Loud et al., Lattice QCD determination of gA ,
energy particle physics, Eur. Phys. J. C 77, 112 (2017). Proc. Sci. CD2018 (2020) 020 [arXiv:1912.08321].
[74] M. Ademollo and R. Gatto, Nonrenormalization Theorem [91] Y.-C. Jang, R. Gupta, T. Bhattacharya, B. Yoon, and H.-W.
for the Strangeness Violating Vector Currents, Phys. Rev. Lin (PNDME Collaboration), Nucleon isovector axial
Lett. 13, 264 (1964). form factors, arXiv:2305.11330.
[75] R. L. Workman et al. (Particle Data Group), Review of [92] H.-W. Lin, Nucleon structure and hyperon form factor
particle physics, Prog. Theor. Exp. Phys. 2022, 083C01 from lattice QCD, arXiv:0707.3844.
(2022). [93] C. Alexandrou, K. Hadjiyiannakou, K. Jansen, and C.
[76] L. Hayen, Standard model OðαÞ renormalization of gA and Kallidonis, Sigma-terms and axial charges for hyperons
its impact on new physics searches, Phys. Rev. D 103, and charmed baryons, Proc. Sci. Lattice2013 (2014) 279
113001 (2021). [arXiv:1404.6355].
[77] V. Cirigliano, J. de Vries, L. Hayen, E. Mereghetti, and [94] M. Göckeler, P. Hägler, R. Horsley, Y. Nakamura, D. Pleiter,
A. Walker-Loud, Pion-Induced Radiative Corrections to P. E. L. Rakow, A. Schäfer, G. Schierholz, H. Stüben, and
Neutron β Decay, Phys. Rev. Lett. 129, 121801 (2022). J. M. Zanotti (QCDSF-UKQCD Collaboration), Baryon
[78] Y. Aoki et al. (Flavour Lattice Averaging Group (FLAG)), axial charges and momentum fractions with N f ¼ 2 þ 1
FLAG review 2021, Eur. Phys. J. C 82, 869 (2022). dynamical fermions, Proc. Sci. LATTICE2010 (2010) 163
[79] J. Liang, Y.-B. Yang, T. Draper, M. Gong, and K.-F. Liu [arXiv:1102.3407].
(χQCD Collaboration), Quark spins and anomalous Ward [95] K.-S. Choi, W. Plessas, and R. F. Wagenbrunn, Axial
identity, Phys. Rev. D 98, 074505 (2018). charges of octet and decuplet baryons, Phys. Rev. D 82,
[80] M. González-Alonso and J. Martin Camalich, Isospin Break- 014007 (2010).
ing in the Nucleon Mass and the Sensitivity of β Decays to [96] X. Y. Liu, D. Samart, K. Khosonthongkee, A. Limphirat,
New Physics, Phys. Rev. Lett. 112, 042501 (2014). K. Xu, and Y. Yan, Axial charges of octet and decuplet
[81] N. Yamanaka, S. Hashimoto, T. Kaneko, and H. Ohki baryons in a perturbative chiral quark model, Phys. Rev. C
(JLQCD Collaboration), Nucleon charges with dynamical 97, 055206 (2018).
overlap fermions, Phys. Rev. D 98, 054516 (2018). [97] R. Qi, J.-B. Wang, G. Li, C.-S. An, C.-R. Deng, and J.-J.
[82] E. Shintani, K.-I. Ishikawa, Y. Kuramashi, S. Sasaki, and T. Xie, Investigations on the flavor-dependent axial charges
Yamazaki (PACS Collaboration), Nucleon form factors of the octet baryons, Phys. Rev. C 105, 065204 (2022).
and root-mean-square radii on a ð10.8 fmÞ4 lattice at the [98] G.-S. Yang and H.-C. Kim, Hyperon semileptonic decay
physical point, Phys. Rev. D 99, 014510 (2019); 102, constants with flavor SU(3) symmetry breaking, Phys.
019902(E) (2020). Rev. C 92, 035206 (2015).
[83] N. Hasan, J. Green, S. Meinel, M. Engelhardt, S. Krieg, J. [99] F. Herren and M. Steinhauser, Version 3 of RunDec and
Negele, A. Pochinsky, and S. Syritsyn (LHPC Collabora- CRunDec, Comput. Phys. Commun. 224, 333 (2018).
tion), Nucleon axial, scalar, and tensor charges using lattice [100] M. Bruno, M. Dalla Brida, P. Fritzsch, T. Korzec, A.
QCD at the physical pion mass, Phys. Rev. D 99, 114505 Ramos, S. Schaefer, H. Simma, S. Sint, and R. Sommer
(2019). (ALPHA Collaboration), QCD Coupling from a Non-
[84] S. Park, R. Gupta, B. Yoon, S. Mondal, T. Bhattacharya, perturbative Determination of the Three-Flavor Λ Param-
Y.-C. Jang, B. Jo’o, and F. Winter (NME Collaboration), eter, Phys. Rev. Lett. 119, 102001 (2017).
Precision nucleon charges and form factors using (2 þ 1)- [101] P. A. Baikov, K. G. Chetyrkin, and J. H. Kühn, Five-Loop
flavor lattice QCD, Phys. Rev. D 105, 054505 (2022). Running of the QCD Coupling Constant, Phys. Rev. Lett.
[85] K. Ottnad, D. Djukanovic, T. Harris, H. B. Meyer, G. 118, 082002 (2017).
von Hippel, and H. Wittig, Improved analysis of nucleon [102] P. A. Baikov, K. G. Chetyrkin, and J. H. Kühn, Five-loop
isovector charges and twist-2 matrix elements on CLS fermion anomalous dimension for a general gauge group
034512-45
GUNNAR S. BALI et al. PHYS. REV. D 108, 034512 (2023)
from four-loop massless propagators, J. High Energy Phys. [118] W. N. Cottingham, The neutron proton mass difference and
04 (2017) 119. electron scattering experiments, Ann. Phys. (N.Y.) 25, 424
[103] Y. Schröder and M. Steinhauser, Four-loop decoupling (1963).
relations for the strong coupling, J. High Energy Phys. 01 [119] R. F. Dashen, Chiral SUð3Þ × SUð3Þ as a symmetry of the
(2006) 051. strong interactions, Phys. Rev. 183, 1245 (1969).
[104] T. Liu and M. Steinhauser, Decoupling of heavy quarks at [120] J. Gasser, A. Rusetsky, and I. Scimemi, Electromagnetic
four loops and effective Higgs-fermion coupling, Phys. corrections in hadronic processes, Eur. Phys. J. C 32, 97
Lett. B 746, 330 (2015). (2003).
[105] P. Marquard, A. V. Smirnov, V. A. Smirnov, M. [121] D. Giusti, V. Lubicz, C. Tarantino, G. Martinelli, F.
Steinhauser, and D. Wellmann, MS-on-shell quark mass Sanfilippo, S. Simula, and N. Tantalo (RM123 Collabo-
relation up to four loops in QCD and a general SUðNÞ ration), Leading isospin-breaking corrections to pion, kaon
gauge group, Phys. Rev. D 94, 074025 (2016). and charmed-meson masses with twisted-mass fermions,
[106] M. Abramczyk, T. Blum, T. Izubuchi, C. Jung, M. Lin, A. Phys. Rev. D 95, 114504 (2017).
Lytle, S. Ohta, and E. Shintani (RBC and UKQCD [122] A. Bazavov et al. (Fermilab Lattice and MILC Collabo-
Collaborations), Nucleon mass and isovector couplings rations), B- and D-meson leptonic decay constants from
in 2 þ 1-flavor dynamical domain-wall lattice QCD near four-flavor lattice QCD, Phys. Rev. D 98, 074512 (2018).
physical mass, Phys. Rev. D 101, 034510 (2020). [123] S. Basak et al. (MILC Collaboration), Lattice computation
[107] L. Liu, T. Chen, T. Draper, J. Liang, K.-F. Liu, G. Wang, of the electromagnetic contributions to kaon and pion
and Y.-B. Yang (χQCD Collaboration), Nucleon isovector masses, Phys. Rev. D 99, 034503 (2019).
scalar charge from overlap fermions, Phys. Rev. D 104, [124] J. Gasser, H. Leutwyler, and A. Rusetsky, On the mass
094503 (2021). difference between proton and neutron, Phys. Lett. B 814,
[108] Z. Fodor, C. Hoelbling, S. Krieg, L. Lellouch, T. Lippert, 136087 (2021).
A. Portelli, A. Sastre, K. K. Szabo, and L. Varnhorst, Up [125] M. G. Endres, A. Shindler, B. C. Tiburzi, and A. Walker-
and Down Quark Masses and Corrections to Dashen’s Loud, Massive Photons: An Infrared Regularization
Theorem from Lattice QCD and Quenched QED, Phys. Scheme for Lattice QCD þ QED, Phys. Rev. Lett. 117,
Rev. Lett. 117, 082001 (2016). 072002 (2016).
[109] S. Borsanyi et al., Ab initio calculation of the neutron- [126] D. A. Brantley, B. Joó, E. V. Mastropas, E. Mereghetti, H.
proton mass difference, Science 347, 1452 (2015). Monge-Camacho, B. C. Tiburzi, and A. Walker-Loud,
[110] D. Horkel, Y. Bi, M. Constantinou, T. Draper, J. Liang, Strong isospin violation and chiral logarithms in the
K.-F. Liu, Z. Liu, and Y.-B. Yang (χQCD Collaboration), baryon spectrum, arXiv:1612.07733.
Nucleon isovector tensor charge from lattice QCD using [127] R. Horsley et al. (CSSM-QCDSF-UKQCD Collaboration),
chiral fermions, Phys. Rev. D 101, 094501 (2020). Isospin splittings in the decuplet baryon spectrum from
[111] T. Bhattacharya, V. Cirigliano, S. Cohen, R. Gupta, A. dynamical QCD þ QED, J. Phys. G 46, 115004 (2019).
Joseph, H.-W. Lin, and B. Yoon (PNDME Collaboration), [128] R. Horsley, J. Najjar, Y. Nakamura, H. Perlt, D. Pleiter,
Iso-vector and iso-scalar tensor charges of the nucleon P. E. L. Rakow, G. Schierholz, A. Schiller, H. Stüben, and
from lattice QCD, Phys. Rev. D 92, 094511 (2015). J. M. Zanotti (QCDSF-UKQCD Collaboration), Lattice
[112] C. Alexandrou et al. (ETM Collaboration), Moments of the determination of Sigma-Lambda mixing, Phys. Rev. D
nucleon transverse quark spin densities using lattice QCD, 91, 074512 (2015).
Phys. Rev. D 107, 054504 (2023). [129] S. Dürr et al., Lattice Computation of the Nucleon Scalar
[113] M. Anselmino, M. Boglione, U. D’Alesio, S. Melis, F. Quark Contents at the Physical Point, Phys. Rev. Lett. 116,
Murgia, and A. Prokudin, Simultaneous extraction of 172001 (2016).
transversity and Collins functions from new SIDIS and [130] N. Salone, P. Adlarson, V. Batozskaya, A. Kupsc, S.
eþ e− data, Phys. Rev. D 87, 094019 (2013). Leupold, and J. Tandean, Study of CP violation in hyperon
[114] Z.-B. Kang, A. Prokudin, P. Sun, and F. Yuan, Extraction decays at super-charm-tau factories with a polarized
of quark transversity distribution and Collins fragmenta- electron beam, Phys. Rev. D 105, 116022 (2022).
tion functions with QCD evolution, Phys. Rev. D 93, [131] P. L. J. Petrak, G. Bali, S. Collins, J. Heitger, D. Jenkins,
014009 (2016). and S. Weishäupl, Sigma terms of the baryon octet in N f ¼
[115] M. Radici and A. Bacchetta, First Extraction of Trans- 2 þ 1 QCD with Wilson quarks, Proc. Sci. LATTICE2022
versity from a Global Analysis of Electron-Proton and (2023) 112 [arXiv:2301.03871].
Proton-Proton Data, Phys. Rev. Lett. 120, 192001 (2018). [132] Jülich Supercomputing Centre, JUWELS: Modular tier-
[116] J. Benel, A. Courtoy, and R. Ferro-Hernandez, A con- 0=1 supercomputer at the Jülich Supercomputing Centre,
strained fit of the valence transversity distributions from J. Large-Scale Res. Facil. 5, A135 (2019).
dihadron production, Eur. Phys. J. C 80, 465 (2020). [133] Jülich Supercomputing Centre, JURECA: Modular super-
[117] J. Cammarota, L. Gamberg, Z.-B. Kang, J. A. Miller, D. computer at Jülich Supercomputing Centre, J. Large-Scale
Pitonyak, A. Prokudin, T. C. Rogers, and N. Sato (JAM Res. Facil. 4, A132 (2018).
Collaboration), Origin of single transverse-spin asymme- [134] Jülich Supercomputing Centre, HDF cloud—Helmholtz
tries in high-energy collisions, Phys. Rev. D 102, 054002 data federation cloud resources at the Jülich Supercomput-
(2020). ing Centre, J. Large-Scale Res. Facil. 5, A137 (2019).
034512-46
OCTET BARYON ISOVECTOR CHARGES FROM N f ¼ 2 þ 1 … PHYS. REV. D 108, 034512 (2023)
[135] R. G. Edwards and B. Joó (SciDAC, LHP, and UKQCD multigrid for the lattice Wilson Dirac operator, SIAM J.
Collaborations), The Chroma software system for lattice Sci. Comput. 36, A1581 (2014).
QCD, Nucl. Phys. B, Proc. Suppl. 140, 832 (2005). [139] M. Lüscher, Deflation acceleration of lattice QCD simu-
[136] S. Heybrock, M. Rottmann, P. Georg, and T. Wettig, lations, J. High Energy Phys. 12 (2007) 011.
Adaptive algebraic multigrid on SIMD architectures, Proc. [140] J. D. Hunter, Matplotlib: A 2d graphics environment, Comput. Sci.
Sci. LATTICE2015 (2016) 036 [arXiv:1512.04506]. Eng. 9, 90 (2007).
[137] P. Georg, D. Richtmann, and T. Wettig, DD-αAMG on [141] A. Agadjanov, D. Djukanovic, G. von Hippel, H. B.
QPACE 3, EPJ Web Conf. 175, 02007 (2018). Meyer, K. Ottnad, and H. Wittig, The nucleon sigma
[138] A. Frommer, K. Kahl, S. Krieg, B. Leder, and M. Rottmann, terms with N f ¼ 2 þ 1 OðaÞ-improved Wilson fermions,
Adaptive aggregation based domain decomposition arXiv:2303.08741.
034512-47