Source author record

Terrence Draper

Terrence Draper appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2023arXiv

Nucleon Electric Dipole Moment from the $θ$ Term with Lattice Chiral Fermions

We calculate the nucleon electric dipole moment (EDM) from the $θ$ term with overlap fermions on three domain wall lattices with different sea pion masses at lattice spacing 0.11 fm. Due to the chiral symmetry conserved by the overlap fermions, we have well defined topological charge and chiral limit for the EDM. Thus, the chiral extrapolation can be carried out reliably at nonzero lattice spacings. We use three to four different partially quenched valence pion masses for each sea pion mass and find that the EDM dependence on the valence and sea pion masses behaves oppositely, which can be described by partially quenched chiral perturbation theory. With the help of the cluster decomposition error reduction (CDER) technique, we determine the neutron and proton EDM at the physical pion mass to be $d_{n}=-0.00148\left(14\right)\left(31\right)\barθ$ e$\cdot$fm and $d_{p}=0.0038\left(11\right)\left(8\right)\barθ$ e$\cdot$fm. This work is a clear demonstration of the advantages of using chiral fermions in the nucleon EDM calculation and paves the road to future precise studies of the strong $CP$ violation effects.

preprint2022arXiv

Proton momentum and angular momentum decompositions with overlap fermions

We present a calculation of the proton momentum and angular momentum decompositions using overlap fermions on a $2+1$-flavor RBC/UKQCD domain-wall lattice at 0.143 fm with a pion mass of 171 MeV which is close to the physical one. A complete determination of the momentum and angular momentum fractions carried by up, down, strange and glue inside the proton has been done with valence pion masses varying from 171 to 391 MeV. We have utilized fast Fourier transform on the stochastic-sandwich method for connected-insertion parts and the cluster-decomposition error reduction technique for disconnected-insertion parts has been used to reduce statistical errors. The full nonperturbative renormalization and mixing between the quark and glue operators are carried out. The final results are normalized with the momentum and angular momentum sum rules and reported at the physical valence pion mass at ${\overline{\rm {MS}}}\, (μ= 2\ {\rm{GeV}})$. The renormalized momentum fractions for the quarks and glue are $\langle x \rangle^q = 0.491(20)(23)$ and $\langle x \rangle^g = 0.509(20)(23)$, respectively, and the renormalized total angular momentum fractions for quarks and glue are $2 J^q = 0.539(22)(44)$ and $2 J^g = 0.461(22)(44)$, respectively. The quark spin fraction is $Σ= 0.405(25)(37)$ from our previous work and the quark orbital angular momentum fraction is deduced from $2 L^q = 2 J^q - Σ$ to be $0.134(22)(44)$.

preprint2022arXiv

RI/MOM and RI/SMOM renormalization of quark bilinear operators using overlap fermions

We present the vector, scalar and tensor renormalization constants (RCs) using overlap fermions with either regularization independent momentum subtraction (RI/MOM) or symmetric momentum subtraction (RI/SMOM) as the intermediate scheme on the lattice with lattice spacings $a$ from 0.04 fm to 0.12 fm. Our gauge field configurations from the MILC and RBC/UKQCD collaborations include sea quarks using either the domain wall or the HISQ action, respectively. The results show that RI/MOM and RI/SMOM can provide consistent renormalization constants to the $\overline{\textrm{MS}}$ scheme, after proper $a^2p^2$ extrapolations. But at $p\sim 2$\,GeV, both RI/MOM and RI/SMOM suffer from nonperturbative effects which cannot be removed by the perturbative matching. The comparison between the results with different sea actions also suggests that the renormalization constant is discernibly sensitive to the lattice spacing but not to the bare gauge coupling in the gauge action.

preprint2020arXiv

Constraints on charm-anticharm asymmetry in the nucleon from lattice QCD

We present the first lattice QCD calculation of the charm quark contribution to the nucleon electromagnetic form factors $G^c_{E,M}(Q^2)$ in the momentum transfer range $0\leq Q^2 \leq 1.4$ $\rm GeV^2$. The quark mass dependence, finite lattice spacing and volume corrections are taken into account simultaneously based on the calculation on three gauge ensembles including one at the physical pion mass. The nonzero value of the charm magnetic moment $μ^c_M=-0.00127(38)_{\rm stat}(5)_{\rm sys}$, as well as the Pauli form factor, reflects a nontrivial role of the charm sea in the nucleon spin structure. The nonzero $G^c_{E}(Q^2)$ indicates the existence of a nonvanishing asymmetric charm-anticharm sea in the nucleon. Performing a nonperturbative analysis based on holographic QCD and the generalized Veneziano model, we study the constraints on the $[c(x)-\bar{c}(x)]$ distribution from the lattice QCD results presented here. Our results provide complementary information and motivation for more detailed studies of physical observables that are sensitive to intrinsic charm and for future global analyses of parton distributions including asymmetric charm-anticharm distribution.

preprint2020arXiv

Ratio of strange to $u/d$ momentum fraction in disconnected insertions

The ratio of the strange quark momentum fraction $\langle x\rangle_{s+\bar{s}}$ to that of light quark $u$ or $d$ in disconnected insertions (DI) is calculated on the lattice with overlap fermions on four domain wall fermion ensembles. These ensembles cover three lattice spacings, three volumes and several pion masses including the physical one, from which a global fitting is carried out. A complete nonperturbative renormalization and the mixing between the quark and glue operators are taken into account. We find the ratio to be $\langle x\rangle_{s+\bar{s}}/\langle x\rangle_{u+\bar{u}} ({\rm DI})=0.795(79)(77)$ at $μ= 2$ GeV in the $\overline{\rm MS}$ scheme. This ratio can be used as a constraint to better determine the strange parton distribution especially in the small $x$ region in the global fittings of PDFs when the connected and disconnected sea are fitted and evolved separately, demonstrating a new way that connects lattice calculations with global analyses.

preprint2020arXiv

Roper State from Overlap Fermions

The Roper state is extracted with valence overlap fermions on a $2+1$-flavor domain-wall fermion lattice (spacing $a = 0.114$ fm and $m_π = 330$ MeV) using both the Sequential Empirical Bayes (SEB) method and the variational method. The results are consistent, provided that a large smearing-size interpolation operator is included in the variational calculation to have better overlap with the lowest radial excitation. Similar calculations carried out for an anisotropic clover lattice with similar parameters find the Roper $\approx 280$ MeV higher than that of the overlap fermion. The fact that the prediction of the Roper state by overlap fermions is consistently lower than those of clover fermions, chirally improved fermions, and twisted-mass fermions over a wide range of pion masses has been dubbed a "Roper puzzle." To understand the origin of this difference, we study the hairpin $Z$-diagram in the isovector scalar meson ($a_0$) correlator in the quenched approximation. Comparing the $a_0$ correlators for clover and overlap fermions, at a pion mass of 290 MeV, we find that the spectral weight of the ghost state with clover fermions is smaller than that of the overlap at $a = 0.12$ fm and $0.09$ fm, whereas the whole $a_0$ correlators of clover and overlap at $a = 0.06$ fm coincide within errors. This suggests that chiral symmetry is restored for clover at $a \le 0.06$ fm and that the Roper should come down at and below this $a$. We conclude that this work supports a resolution of the "Roper puzzle" due to $Z$-graph type chiral dynamics. This entails coupling to higher components in the Fock space (e.g. $Nπ$, $Nππ$ states) to induce the effective flavor-spin interaction between quarks as prescribed in the chiral quark model, resulting in the parity-reversal pattern as observed in the experimental excited states of $N, Δ$ and $Λ$.

preprint2020arXiv

The nucleon isovector tensor charge from lattice QCD using chiral fermions

In this work we present the isovector flavor combination for the nucleon tensor charge extracted from lattice QCD simulations using overlap fermions on $N_f=2+1$ domain-wall configurations. The pion mass dependence is studied using six valence quark masses, each reproducing a value for the pion mass in the valence sector between 147 and 330 MeV. We investigate and eliminate systematic uncertainties due to contamination by excited states, by employing several values for the source-sink separation that span from 1 fm to 1.6 fm. We apply a chiral extrapolation in the valence sector using a quadratic and a logarithmic term to fit the pion mass dependence, which describes well the lattice data. The lattice matrix element is renormalized non-perturbatively, and the final result is $g_T=1.096(30)$ in the $\overline{\rm MS}$ scheme at a renormalization scale of 2 GeV.

preprint2019arXiv

Towards the nucleon hadronic tensor from lattice QCD

We present the first calculation of the hadronic tensor on the lattice for the nucleon. The hadronic tensor can be used to extract the structure functions in deep inelastic scatterings and also provide information for the neutrino-nucleon scattering which is crucial to the neutrino-nucleus scattering experiments at low energies. The most challenging part in the calculation is to solve an inverse problem. We have implemented and tested three algorithms using mock data, showing that the Bayesian Reconstruction method has the best resolution in extracting peak structures while the Backus-Gilbert and Maximum Entropy methods are somewhat more stable for the flat spectral function. Numerical results are presented for both the elastic case (clover fermions on domain wall configuration with $m_π\sim$ 370 MeV and $a\sim$ 0.06 fm) and a case (anisotropic clover lattice with $m_π\sim$ 380 MeV and $a_t\sim$ 0.035 fm) with large momentum transfer. For the former case, the reconstructed Minkowski hadronic tensor gives precisely the vector charge which proves the feasibility of the approach. While for the latter case, the nucleon resonances and possibly shallow inelastic scattering contributions around $ν=1$ GeV are clearly observed but no information is obtained for higher excited states with $ν>2$ GeV. A check of the effective masses of $ρ$ meson with different lattice setups indicates that, in order to reach higher energy transfers, using lattices with smaller lattice spacings is essential.

preprint2016arXiv

$π$N and strangeness sigma terms at the physical point with chiral fermions

Lattice QCD calculations with chiral fermions of the $π$N sigma term $σ_{πN}$ and strangeness sigma term $σ_{sN}$ including chiral interpolation with continuum and volume corrections are provided in this work, with the excited-state contaminations subtracted properly. We calculate the scalar matrix element for the light/strange quark directly and find $σ_{πN}=45.9(7.4)(2.8)$ MeV, with the disconnected insertion part contributing 20(12)(4)\%, and $σ_{sN}=40.2(11.7)(3.5)$ MeV, which is somewhat smaller than $σ_{πN}$. The ratio of the strange/light scalar matrix elements is $y$ = 0.09(3)(1).

preprint2016arXiv

Glue Spin of the Proton

We report the progress on the lattice QCD calculation of the glue spin contribution to proton spin. This calculation is carried out with valence overlap fermion on 2+1 flavor DWF gauge configurations at two lattice spacings with the momentum of the frame in the kinematic range $0\leq p^2 \leq 2$ GeV$^2$. A mild frame dependence is observed. The matching and mixing with large-momentum effective field theory are in progress. The unrenormalized result at $p^2=4$ GeV$^2$ with $O(a^2)$ correction gives $S_G$ = 0.13(3).

preprint2016arXiv

Stochastic method with low mode substitution for nucleon isovector matrix elements

We introduce a stochastic sandwich method with low-mode substitution to evaluate the connected three-point functions. The isovector matrix elements of the nucleon for the axial-vector coupling $g_A^3$, scalar couplings $g_S^3$ and the quark momentum fraction $\langle x\rangle_{u -d}$ are calculated with overlap fermion on 2+1 flavor domain-wall configurations on a $24^3 \times 64$ lattice at $m_π = 330$ MeV with lattice spacing $a = 0.114$ fm.

preprint2015arXiv

Charm and strange quark masses and $f_{D_s}$ from overlap fermions

We use overlap fermions as valence quarks to calculate meson masses in a wide quark mass range on the $2+1$-flavor domain-wall fermion gauge configurations generated by the RBC and UKQCD Collaborations. The well-defined quark masses in the overlap fermion formalism and the clear valence quark mass dependence of meson masses observed from the calculation facilitate a direct derivation of physical current quark masses through a global fit to the lattice data, which incorporates $O(a^2)$ and $O(m_c^4a^4)$ corrections, chiral extrapolation, and quark mass interpolation. Using the physical masses of $D_s$, $D_s^*$ and $J/ψ$ as inputs, Sommer's scale parameter $r_0$ and the masses of charm quark and strange quark in the $\overline{\rm MS}$ scheme are determined to be $r_0=0.465(4)(9)$ fm, $m_c^{\overline{\rm MS}}(2\,{\rm GeV})=1.118(6)(24)$ GeV (or $m_c^{\overline{\rm MS}}(m_c)=1.304(5)(20)$ GeV), and $m_s^{\overline{\rm MS}}(2\,{\rm GeV})=0.101(3)(6)\,{\rm GeV}$, respectively. Furthermore, we observe that the mass difference of the vector meson and the pseudoscalar meson with the same valence quark content is proportional to the reciprocal of the square root of the valence quark masses. The hyperfine splitting of charmonium, $M_{J/ψ}-M_{η_c}$, is determined to be 119(2)(7) MeV, which is in good agreement with the experimental value. We also predict the decay constant of $D_s$ to be $f_{D_s}=254(2)(4)$ MeV. The masses of charmonium $P$-wave states $χ_{c0}, χ_{c1}$ and $h_c$ are also in good agreement with experiments.

preprint2015arXiv

Meson Mass Decomposition from Lattice QCD

Hadron masses can be decomposed as a sum of quark and glue components which are defined through hadronic matrix elements of QCD operators. The components consist of the quark mass term, the quark energy term, the glue energy term, and the trace anomaly term. We calculate these components for mesons with lattice QCD for the first time. The calculation is carried out with overlap fermion on $2+1$ flavor domain-wall fermion gauge configurations. We confirm that $\sim 50\%$ of the light pion mass comes from the quark mass term and $\sim 10\%$ comes from the quark energy; whereas, while for the $ρ$ meson, the quark energy contributes roughly half of its mass but the quark mass term contributes little. The combined glue components contribute $\sim 40 - 50\%$ for both mesons. It is interesting to observe that the quark mass contribution to the mass of the vector meson is almost linear in quark mass over a large quark mass region below the charm quark mass. For heavy mesons, the quark mass term dominates the masses, while the contribution from the glue components is about $200$ MeV (a bare value around 2GeV) for the heavy pseudoscalar and vector mesons. The charmonium hyperfine splitting is found to be dominated by the quark energy term which is consistent with the picture of the quark potential model.

preprint2014arXiv

Meson Mass Decomposition

Hadron masses can be decomposed as a sum of components which are defined through hadronic matrix elements of QCD operators. The components consist of the quark mass term, the quark energy term, the glue energy term and the trace anomaly term. We calculate these components of mesons with lattice QCD for the first time. The calculation is carried out with overlap fermion on $2+1$ flavor domain-wall fermion gauge configurations. We confirm that $\sim 50\%$ of the light pion mass comes from the quark mass and $\sim 10\%$ comes from the quark energy, whereas, the contributions are found to be the other way around for the $ρ$ mass. The combined glue components contribute $\sim 40 - 50\%$ for both mesons. It is interesting to observe that the quark mass contribution to the mass of the vector meson is almost linear in quark mass over a large quark mass region below the charm quark mass. For heavy mesons, the quark mass term dominates the masses, while the contribution from the glue components is about $400\sim500$ MeV for the heavy pseudoscalar and vector mesons. The charmonium hyperfine splitting is found to be dominated by the quark energy term which is consistent with the quark potential model.

preprint2011arXiv

The Analysis of Space-Time Structure in QCD Vacuum II: Dynamics of Polarization and Absolute X-Distribution

We propose a framework for quantitative evaluation of dynamical tendency for polarization in arbitrary random variable that can be decomposed into a pair of orthogonal subspaces. The method uses measures based on comparisons of given dynamics to its counterpart with statistically independent components. The formalism of previously considered X-distributions is used to express the aforementioned comparisons, in effect putting the former approach on solid footing. Our analysis leads to definition of a suitable correlation coefficient with clear statistical meaning. We apply the method to the dynamics induced by pure-glue lattice QCD in local left-right components of overlap Dirac eigenmodes. It is found that, in finite physical volume, there exists a non-zero physical scale in the spectrum of eigenvalues such that eigenmodes at smaller (fixed) eigenvalues exhibit convex X-distribution (positive correlation), while at larger eigenvalues the distribution is concave (negative correlation). This chiral polarization scale thus separates a regime where dynamics enhances chirality relative to statistical independence from a regime where it suppresses it, and gives an objective definition to the notion of "low" and "high" Dirac eigenmode. We propose to investigate whether the polarization scale remains non-zero in the infinite volume limit, in which case it would represent a new kind of low energy scale in QCD.

preprint2010arXiv

Absolute Measure of Local Chirality and the Chiral Polarization Scale of the QCD Vacuum

The use of the absolute measure of local chirality is championed since it has a uniform distribution for randomly reshuffled chiral components so that any deviations from uniformity in the associated "X-distribution" are directly attributable to QCD-induced dynamics. We observe a transition in the qualitative behavior of this absolute X-distribution of low-lying eigenmodes which, we propose, defines a chiral polarization scale of the QCD vacuum.

preprint2010arXiv

Nucleon strangeness form factors and moments of PDF

The calculation of the nucleon strangeness form factors from N_f=2+1 clover fermion lattice QCD is presented. Disconnected insertions are evaluated using the Z(4) stochastic method, along with unbiased subtractions from the hopping parameter expansion. We find that increasing the number of nucleon sources for each configuration improves the signal significantly. We obtain G_M^s(0) = -0.017(25)(07), which is consistent with experimental values, and has an order of magnitude smaller error. Preliminary results for the strangeness contribution to the second moment of the parton distribution function are also presented.

preprint2009arXiv

Nucleon strangeness form factors from N_f=2+1 clover fermion lattice QCD

We present the N_f=2+1 clover fermion lattice QCD calculation of the nucleon strangeness form factors. We evaluate disconnected insertions using the Z(4) stochastic method, along with unbiased subtractions from the hopping parameter expansion. We find that increasing the number of nucleon sources for each configuration improves the signal significantly. We obtain G_M^s(0) = -0.017(25)(07), where the first error is statistical, and the second is the uncertainties in Q^2 and chiral extrapolations. This is consistent with experimental values, and has an order of magnitude smaller error.

preprint2009arXiv

The calculation of nucleon strangeness form factors from N_f=2+1 clover fermion lattice QCD

We study the strangeness electromagnetic form factors of the nucleon from the N_f=2+1 clover fermion lattice QCD calculation. The disconnected insertions are evaluated using the Z(4) stochastic method, along with unbiased subtractions from the hopping parameter expansion. In addition to increasing the number of Z(4) noises, we find that increasing the number of nucleon sources for each configuration improves the signal significantly. We obtain G_M^s(0) = -0.017(25)(07), where the first error is statistical, and the second is the uncertainties in Q^2 and chiral extrapolations. This is consistent with experimental values, and has an order of magnitude smaller error. We also study the strangeness second moment of the partion distribution function of the nucleon, <x^2>_{s-\bar{s}}.

preprint2008arXiv

2+1 flavor QCD calculation of <x> and <x^2>

We calculate the connected insertions of the nucleon three-point function to study the first few moments of the unpolarized structure functions of the nucleon. (The disconnected insertions are discussed elsewhere in these proceedings). The calculation employs the CP-PACS/JLQCD 2+1 dynamical clover fermions on a 16^3x32 lattice with lattice spacing a=0.1219 fm. The sequential source technique, using non-zero and zero momentum point nucleon field as the secondary source, is applied enabling a study of different currents at various momentum transfer.

preprint2008arXiv

Light scalar mesons in 2+1 flavor full QCD

We study the a_0 and K_0^* light scalar mesons in 2+1 flavor full QCD. Particular attention is paid to fitting excited states, with an eye toward determining whether scattering states are revealed. An ultimate goal will be to see how dynamical quarks affect the picture outlined with an earlier quenched study using overlap fermions, namely, that it is the a_0(1450), not the a_0(980), which is the lowest $\bar{q}q$ isovector scalar state.

preprint2008arXiv

Strangeness and glue in the nucleon from lattice QCD

We study the strangeness contribution to nucleon matrix elements using Nf=2+1 dynamical clover fermion configurations generated by the CP-PACS/JLQCD collaboration. In order to evaluate the disconnected insertion (DI), we use the Z(4) stochastic method, along with unbiased subtraction from the hopping parameter expansion which reduces the off-diagonal noises in the stochastic method. Furthermore, we find that using many nucleon sources for each configuration is effective in improving the signal. Our results for the quark contribution to the first moment <x>_q in the DI, and the strangeness magnetic moment show that the statistical errors are under control with these techniques. We also study the gluonic contribution to the nucleon using the overlap operator to construct the gauge field tensor, F_{mu,nu}. The application to the calculation of first moment, <x>_G, gives a good signal in quenched lattice QCD.

preprint2005arXiv

Locality and Scaling of Quenched Overlap Fermions

The overlap fermion offers the tremendous advantage of exact chiral symmetry on the lattice, but is numerically intensive. This can be made affordable while still providing large lattice volumes, by using coarse lattice spacing, given that good scaling and localization properties are established. Here, using overlap fermions on quenched Iwasaki gauge configurations, we demonstrate directly that the overlap Dirac operator's range is comfortably small in lattice units for each of the lattice spacings 0.20 fm, 0.17 fm, and 0.13 fm (and scales to zero in physical units in the continuum limit). In particular, our direct results contradict recent speculation that an inverse lattice spacing of $1 {\rm GeV}$ is too low to have satisfactory localization. Furthermore, hadronic masses (available on the two coarser lattices) scale very well.

preprint2004arXiv

Improved Measure of Local Chirality

It is popular to probe the structure of the QCD vacuum indirectly by studying individual fermion eigenmodes, because this provides a natural way to filter out UV fluctuations. The double-peaking in the distribution of the local chiral orientation parameter (X) has been offered as evidence, by some, in support of a particular model of the vacuum. Here we caution that the X-distribution peaking varies significantly with various versions of the definition of X. Furthermore, each distribution varies little from that resulting from a random reshuffling of the left-handed (and independently the right-handed) fields, which destroys any QCD-induced left-right correlation; that is, the double-peaking is mostly a phase-space effect. We propose a new universal definition of the X parameter whose distribution is uniform for randomly reshuffled fields. Any deviations from uniformity for actual data can then be directly attributable to QCD-induced dynamics. We find that the familiar double peak disappears.