# Quasi-PDFs from Twisted mass fermions at the physical point

###### Abstract

We present results for the flavor non-singlet parton distribution functions within the nucleon using the quasi-PDF approach. The lattice calculation is performed by employing the twisted mass formulation and two gauge ensembles, having and dynamical fermions with masses tuned to their physical value. For the physical point ensemble, the unpolarized, helicity and transversity distributions are computed for three values of the nucleon momentum, namely corresponding to [0.83,1.11,1.38] GeV. Upon renormalization, we find that, as the nucleon momentum increases, the lattice results approach the phenomenological distributions resulting from analyses of deep inelastic scattering data, opening a promising path for a direct evaluation of parton distributions from the QCD Lagrangian. For the physical point ensemble, we present preliminary results for the unpolarized distribution extracted from a nucleon boosted by or 0.97 GeV.

Quasi-PDFs from Twisted mass fermions at the physical point

Constantia Alexandrou,
Krzysztof Cichy,
Martha Constantinou,
Kyriakos Hadjiyiannakou,
Karl Jansen,
Aurora Scapellato^{†}^{†}thanks: Speaker. ,
Fernanda Steffens

Department of Physics, University of Cyprus, P.O. Box 20537, 1678 Nicosia, Cyprus

Computation-based Science and Technology Research Center, Cyprus Institute, 20 Kavafi Str., Nicosia 2121, Cyprus

Faculty of Physics, Adam Mickiewicz University, Umultowska 85, 61-614 Poznań, Poland

Temple University, 1925 N. 12th Street, Philadelphia, PA 19122, USA

John von Neumann Institute for Computing (NIC), DESY, Platanenallee 6, D-15738 Zeuthen, Germany

University of Wuppertal, Gaußstr. 20, 42119 Wuppertal, Germany

Institut für Strahlen- und Kernphysik, Rheinische Friedrich-Wilhelms-Universität Bonn, Nussallee 14-16, 53115 Bonn

E-mail: scapellato.aurora@ucy.ac.cy

\abstract@cs

## 1 Introduction

Hadrons are complex systems whose internal structure is determined by the strong interactions among quarks and gluons, known as partons. An open question in particle physics is how the decomposition of the hadron momentum among its constituents emerges from Quantum Chromodynamics (QCD). In this framework, one of the most essential quantities are the parton distribution functions (PDFs), which describe the probability density that a parton carries a collinear momentum fraction of the total momentum of the hadron. PDFs are intrinsically non-perturbative objects and require light-cone dynamics, which however cannot be accomodated on a Euclidean lattice, where a finite lattice spacing is employed. For this reason, global QCD analyses of deep inelastic scattering data, aided by phenomenological models, have been the only possible source of information up to date [Lin:2017snn]. In this work, we report our effort in extracting the -dependence of the isovector () distribution functions for unpolarized and polarized nucleons by using the quasi-PDF approach [Ji:2014gla]. Other proposals to access PDFs from first principles have also been recently suggested, like pseudo-PDFs [Radyushkin:2017cyf] implemented within lattice QCD in Ref. [Orginos:2017kos] and the construction of good Lattice Cross Sections (LCSs) [Ma:2017pxb]. Here we show our results for unpolarized, helicity and transversity PDFs extracted from physical point simulations, also presented in [Alexandrou:2018pbm, Alexandrou:2018eet], and report preliminary results for unpolarized PDFs obtained analyzing a ensemble, where all quark masses are tuned to their physical value. The simulation parameters of the ensembles are listed in Table 1.

[fm] | ||||||||
---|---|---|---|---|---|---|---|---|

2 | 96, 48 | 2.10 | 1.57751 | 0.137290 | 0.0009 | 0.06208(2) | 0.0938(2) | 2.98 |

2+1+1 | 128, 64 | 1.778 | 1.69 | 0.1394265 | 0.00072 | 0.05658(6) | 0.08029(41) | 3.62 |

## 2 Results from physical point simulations

For the ensemble we compute the matrix elements relevant for quasi-PDFs at three nucleon momenta, , and or 0.83, 1.11 and 1.38 GeV, using sequential inversions through the sink and the momentum smearing technique [Bali:2016lva]. More details about the setup and statistics can be found in Refs. [Alexandrou:2018pbm, Alexandrou:2018eet] and [Krzysztof_proc:18], where the effects of excited states on the bare matrix elements are also discussed. The results presented in the following are obtained from our largest source-sink time separation, namely fm, where ground state convergence is observed within our statistical uncertainties, of around . Having computed the lattice matrix elements (see [Alexandrou:2018pbm, Alexandrou:2018eet, Krzysztof_proc:18]), we renormalize them by adopting a RI-type scheme. The renormalization functions (Z-factors) are extracted by imposing the renormalization conditions described in Ref. [Alexandrou:2017huk]. We study the RI-scale dependence on the vertex functions exploring 17 choices for the RI scales, having the form with in order to suppress discretization effects [Alexandrou:2015sea]. We take as final choices the ones for which to avoid contamination from non-perturbative effects. The Z-factors are then converted from RI to the scheme and evolved to GeV using the 1-loop formulae of Ref. [Constantinou:2017sej]. To remove the residual dependence on , we perform a linear fit in for each length of the Wilson line, of the kind , and we take as final Z-factor. In this work, we also study how the vertex functions depend on the smearing applied to the gauge field, as the smearing mitigates the power divergences emerging for non-local operators. In our lattice computation, we apply 3-dimensional stout smearing only to the Wilson line connecting the fermion fields and the result for the Z-factors is shown in Fig.1 for the helicity operator.

As can be seen, smearing clearly reduces the value of Z-factors for and the change from to is more significant than from 5 to 10. In order to renormalize the matrix elements properly, the bare matrix elements have to be produced by applying the same number of iterations of stout smearing. The smearing dependence is considerable at the level of the bare matrix elements [Alexandrou:2017dzj], but the renormalized matrix elements agree as expected and demonstrated in Fig.2.

In the quasi-PDF approach, the first step towards the reconstruction of the physical PDFs is the Fourier transform of the renormalized matrix elements. The renormalized Fourier transform, known as quasi-PDF, is shown in Fig.3 (left plot) by the green curve, for nucleon momentum . The quasi-PDF does not represent yet a physical result, until a matching procedure and target mass corrections (TMCs) are applied to account for a finite value of the momentum employed on the lattice. In this work, we use the perturbative matching formulae of Refs. [Alexandrou:2018pbm, Alexandrou:2018eet] that we developed for the unpolarized, helicity and transversity PDFs, and for the TMCs we use Ref. [Chen:2016utp]. The matching procedure drastically alters the slope of the quasi-PDF in the whole x-range, as can be seen from the red band in Fig.3 (left plot). It shifts and enhances the peak towards the small x-region, reproducing the antisymmetry between quark-antiquark distributions expected from global QCD analyses. TMCs are found to have a small effect, mostly evident in the positive -region.

After applying the matching and TMCs to the quasi-PDFs computed at different momenta, one is ready to check how the lattice PDFs depend on the value of the boost and how they compare with phenomenology. As concerns the momentum dependence we refer to Refs. [Alexandrou:2018pbm, Alexandrou:2018eet], where we find that, as the momentum increases, lattice PDFs move towards the phenomenological PDFs, a feature that enhances our confidence in the quasi-PDF approach. Here we show our results for the largest boost employed in this work. For the unpolarized PDF (right plot in Fig.3) we observe a similar slope as compared with global QCD analyses, while for the helicity we obtain agreement with phenomenology for and for small negative , as can be seen in the left plot of Fig.4. Thus, a fully reliable reconstruction of PDFs from the lattice can become possible for even higher momenta, maintaining control over a number of systematic effects [Krzysztof_proc:18], such as contamination by excited states and cut-off effects. For the transversity PDF shown in Fig.4, we find that the statistical errors of the lattice PDFs are strikingly smaller than the phenomenological fits when use only SIDIS data and this also holds for the SIDIS data constrained with the tensor charge computed from lattice.

## 3 Quasi-PDFs from physical point simulations

As a next step towards a better control of systematic effects, we compute PDFs on a physical point ensemble simulated with mass-degenerate up and down quarks and a strange and charm quark in the sea, since eventually our aim is to take the continuum limit of hadron observables using as complete description of QCD as possible. We employ gauge configurations recently generated by the ETM collaboration [Alexandrou:2018egz] with parameters listed in Table 1. This ensemble, as compared to the previously analyzed by us, has a smaller lattice spacing, bigger volume and a larger value for . Thus, we expect discretization and finite volume effects to be less significant for this setup. The lattice computation of PDFs proceeds along the lines of the one performed on the ensemble and uses the same optimized setup of Ref. [Alexandrou:2018egz], including momentum smearing. In Fig.5 we show the bare matrix elements for the unpolarized PDF obtained from the analysis of configurations with total statistics of 4500 measurements, at fm and nucleon boost GeV.

For the renormalization we use, at this preliminary stage, one RI scale which is set to , with and . A detailed study on the RI scale dependence of the renormalization factors will be carried out in the near future. Having then converted the Z-factors to at 2 GeV and applied the matching kernel of Ref. [Alexandrou:2018pbm] and the TMCs of Ref. [Chen:2016utp], we obtain the result depicted in Fig.6. A nice feature is that even at the value of the nucleon boost employed here, the lattice results exhibit a similar slope as PDFs from global analyses for and approaches zero at . The disagreement with phenomenological data is however expected for a number of reasons, most notably a small value of the nucleon boost, for which higher twist effects may play an important role in the matching of the quasi PDFs. The calculation at higher boosts is one of our essential directions in pursuing the extraction of quark distributions from this physical point ensemble.

## 4 Conclusions and outlook

In this work we report our results on the extraction of PDFs from lattice QCD simulations performed with and ensembles at the physical pion mass. Our results have several nice features, the most important of which is the tendency of the lattice PDFs to approach the phenomenological ones as the nucleon boost increases. This is motivating further studies that may improve the lattice estimates on several points, such as the subtraction of discretization effects in the renormalization functions and the calculation of the conversion factor RI to up to 2-loops. The computation at higher momentum is currently ongoing. The investigation of finite volume effects and cut-off effects at the physical value of quark masses is also planned within our long-term program. Availability of large computational resources to study these systematic effects would be needed for a reliable check.

## 5 Acknowledgment

This work has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Skodowska-Curie grant agreement No 642069 (HPC-LEAP). K.C. is supported by the National Science Centre grant SONATA BIS no. 2016/22/E/ST2/00013. F.S. is funded by the Deutsche Forschungsgemeinschaft (DFG) project number 392578569. M.C. acknowledges financial support by the U.S. Department of Energy, Office of Nuclear Physics, within the framework of the TMD Topical Collaboration, as well as, by the National Science Foundation under Grant No. PHY-1714407. This research used resources of the Oak Ridge Leadership Computing Facility, which is a DOE Office of Science User Facility supported under Contract DE-AC05-00OR22725, Jülich Supercomputing Centre, Prometheus supercomputer at the Academic Computing Centre Cyfronet AGH in Cracow (grant ID quasipdfs), Okeanos supercomputer at the Interdisciplinary Centre for Mathematical and Computational Modelling in Warsaw (grant IDs gb70-17, ga71-22), Eagle supercomputer at the Poznan Supercomputing and Networking Center (grant no. 346).