Source author record

Constantia Alexandrou

Constantia Alexandrou 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

66works
7topics
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

66 published item(s)

preprint2024arXiv

Generalized form factors of the pion and kaon using twisted mass fermions

We present an update on our lattice calculations of the Mellin moments of PDFs and GPDs for the pion and kaon, using momentum-boosted meson states. In particular, we focus on the calculation of the scalar and tensor local operators, and the vector operator with up to three-covariant derivatives. The corresponding matrix elements allow us to extract the scalar and tensor charges, as well as $\langle x^n \rangle$ with $n=1,2,3$. In addition, we introduce momentum transfer between the initial and final meson state that leads to the scalar, vector, and tensor form factors, and the generalized form factors up to three covariant derivatives. The above results are obtained using two ensembles of maximally twisted mass fermions with clover improvement with two degenerate light, a strange, and a charm quark $(N_f=2+1+1)$ at lattice spacings of 0.093 fm and 0.081 fm. The pion mass of the ensembles is about 260 MeV. We study excited-states effects by analyzing four values of the source-sink time separation (1.12 - 1.67 fm). We also examine discretization and volume effects. The lattice data are renormalized non-perturbatively, and the results are presented in the MS scheme at a scale of 2 GeV.

preprint2024arXiv

Investigation of two-particle contributions to nucleon matrix elements

We investigate contributions of excited states to nucleon matrix elements by studying the two- and three-point functions using nucleon and pion-nucleon interpolating fields. This study is carried out using twisted mass fermion ensembles with pion masses 346 MeV and 131 MeV. We compare the results obtained using these two ensembles and show preliminary results for nucleon charges.

preprint2024arXiv

Low-lying baryon masses using twisted mass fermions ensembles at the physical pion mass

We investigate the low-lying baryon spectrum using three $N_f$=2+1+1 ensembles simulated with physical values of the quark masses and lattice spacings of 0.080, 0.069, and 0.057 fm. The ensembles are generated using twisted mass clover-improved fermions and the Iwasaki gauge action. The spatial length is kept approximately the same at about 5.1 fm to 5.5 fm fulfilling the condition $m_πL$> 3.6. We investigate isospin splitting within isospin multiples and verify that for most cases the isospin splitting for these lattice spacing is consistent with zero. In the couple of cases, for which there is a non-zero value, in the continuum limit, the mass splitting goes to zero. The baryon masses are extrapolated to the continuum limit using the three $N_f$=2+1+1 ensembles and are compared to other recent lattice QCD results. For the strange and charm quark masses, we find, respectively, $m_s$(2 GeV)=99.2(2.7) MeV and $m_c$(3 GeV)=1.015(39) GeV. The values predicted for the masses of the doubly charmed $Ξ_{cc}^\star$, $Ω_{cc}$ and $Ω_{cc}^\star$ baryons are 3.676(55) GeV, 3.703(51) GeV and 3.803(50) GeV, respectively, and for the triply charmed $Ω_{ccc}$ baryon is 4.785(71) GeV.

preprint2022arXiv

Lattice QCD Study of Transverse-Momentum Dependent Soft Function

In this work, we perform a lattice QCD study of the intrinsic, rapidity-independent soft function within the framework of large momentum effective theory. The computation is carried out using a gauge ensemble of $N_f=2+1+1$ clover-improved twisted mass fermion. After applying an appropriate renormalization procedure and the removal of significant higher-twist contamination, we obtain the intrinsic soft function that is comparable to the one-loop perturbative result at large external momentum. The determination of the nonperturbative soft function from first principles is crucial to sharpen our understanding of the processes with small transverse momentum such as the Drell-Yan production and the semi-inclusive deep inelastic scattering. Additionally, we calculate the Collins-Soper evolution kernel using the quasi-transverse-momentum-dependent wave function as input.

preprint2022arXiv

Pion-pole contribution to HLbL from twisted mass lattice QCD at the physical point

We report on our computation of the pion transition form factor ${\cal F}_{P\rightarrow γ^*γ^*}$ from twisted mass lattice QCD in order to determine the numerically dominant light pseudoscalar pole contribution in the hadronic light-by-light scattering contribution to the anomalous magnetic moment of the muon $a_μ=(g-2)_μ$. The pion transition form factor is computed directly at the physical point. We present first results for our estimate of the pion-pole contribution with kinematic setup for the pion at rest.

preprint2022arXiv

Proton generalized parton distributions from lattice QCD

Momentum and spatial distributions of quarks and gluons inside hadrons are typically encoded in the so-called generalized parton distributions (GPDs). GPDs are multi-dimensional quantities that are very challenging to extract, both experimentally and within lattice QCD. We present the first lattice results on the $x$-dependence of isovector unpolarized, helicity and transversity GPDs of the proton, obtained from lattice QCD using an ensemble of $N_f=2+1+1$ maximally twisted mass fermions, with pion mass $M_π=260$ MeV and lattice spacing $a\simeq 0.093$ fm. Our calculations use the quasi-distribution formalism and the final distributions are presented in the MS-bar scheme at a renormalization scale of 2 GeV.

preprint2022arXiv

The $πγ\to ππ$ transition and the $ρ$ radiative decay width from lattice QCD

We report a lattice QCD determination of the $πγ\to ππ$ transition amplitude for the $P$-wave, $I=1$ two-pion final state, as a function of the photon virtuality and $ππ$ invariant mass. The calculation was performed with $2+1$ flavors of clover fermions at a pion mass of approximately $320$ MeV, on a $32^3 \times 96$ lattice with $L\approx 3.6$ fm. We construct the necessary correlation functions using a combination of smeared forward, sequential and stochastic propagators, and determine the finite-volume matrix elements for all $ππ$ momenta up to $|\vec{P}|= \sqrt{3} \frac{2π}{L}$ and all associated irreducible representations. In the mapping of the finite-volume to infinite-volume matrix elements using the Lellouch-Lüscher factor, we consider two different parametrizations of the $ππ$ scattering phase shift. We fit the $q^2$ and $s$ dependence of the infinite-volume transition amplitude in a model-independent way using series expansions, and compare multiple different truncations of this series. Through analytic continuation to the $ρ$ resonance pole, we also determine the $πγ\to ρ$ resonant transition form factor and the $ρ$ meson photocoupling, and obtain $|G_{ρπγ}| = 0.0802(32)(20)$.

preprint2022arXiv

Transversity GPDs of the proton from lattice QCD

We present the first calculation of the $x$-dependence of the isovector transversity generalized parton distributions (GPDs) for the proton within lattice QCD. We compute the matrix elements with non-local operators containing a Wilson line. The calculation implements the Breit symmetric frame. The proton momenta are chosen as $0.83,\,1.25,\,1.67$ GeV, and the values of the momentum transfer squared are $0.69,\,1.02$ GeV$^2$. These combinations include cases with zero and nonzero skewness. The calculation is performed using one ensemble of two degenerate-mass light, a strange and a charm quark of maximally twisted mass fermions with a clover term. The lattice results are renormalized non-perturbatively and finally matched to the light-cone GPDs using one-loop perturbation theory within the framework of large momentum effective theory. The final GPDs are given in the $\overline{\rm MS}$ scheme at a scale of 2 GeV. In addition to the individual GPDs, we form the combination of the transversity GPDs that is related to the transverse spin structure of the proton. Finally, we extract the lowest two moments of GPDs and draw a number of important qualitative conclusions.

preprint2022arXiv

Twisted mass gauge ensembles at physical values of the light, strange and charm quark masses

Lattice QCD simulations directly at physical masses of dynamical light, strange and charm quarks are highly desirable especially to remove systematic errors due to chiral extrapolations. However such simulations are still challenging. We discuss the adaption of efficient algorithms, like multi-grid methods or higher order integrators, within the molecular dynamic steps of the Hybrid Monte Carlo algorithm, that are enabling simulations of a new set of gauge ensembles by the Extended Twisted Mass collaboration (ETMC). We present the status of the on-going ETMC simulation effort that aims to enabling studies of finite size and discretization effects. We work within the twisted mass discretization which is free of odd-discretization effects at maximal twist and present our tuning procedure.

preprint2022arXiv

Using classical bit-flip correction for error mitigation including 2-qubit correlations

We present an error mitigation scheme which corrects readout errors on Noisy Intermediate-Scale Quantum (NISQ) computers [1,2]. After a short review of applying the method to one qubit, we proceed to discuss the case when correlations between different qubits occur. We demonstrate how the readout error can be mitigated in this case. By performing experiments on IBMQ hardware, we show that such correlations do not have a strong effect on the results, justifying to neglect them.

preprint2021arXiv

Ruling Out the Massless Up-Quark Solution to the Strong CP Problem by Computing the Topological Mass Contribution with Lattice QCD

The infamous strong CP problem in particle physics can in principle be solved by a massless up quark. In particular, it was hypothesized that topological effects could substantially contribute to the observed nonzero up-quark mass without reintroducing CP violation. Alternatively to previous work using fits to chiral perturbation theory, in this Letter, we bound the strength of the topological mass contribution with direct lattice QCD simulations, by computing the dependence of the pion mass on the dynamical strange-quark mass. We find that the size of the topological mass contribution is inconsistent with the massless up-quark solution to the strong CP problem.

preprint2021arXiv

The Mellin moments $\langle x \rangle$ and $\langle x^2 \rangle$ for the pion and kaon from lattice QCD

We present a calculation of the pion quark momentum fraction, $\langle x \rangle$, and its third Mellin moment $\langle x^2 \rangle$. We also obtain directly, for the first time, $\langle x \rangle$ and $\langle x^2 \rangle$ for the kaon using local operators. We use an ensemble of two degenerate light, a strange and a charm quark ($N_f=2+1+1$) of maximally twisted mass fermions with clover improvement. The quark masses are chosen so that they reproduce a pion mass of about 260 MeV, and a kaon mass of 530 MeV. The lattice spacing of the ensemble is 0.093 fm and the lattice has a spatial extent of 3 fm. We analyze several values of the source-sink time separation within the range of $1.12-2.23$ fm to study and eliminate excited-states contributions. The necessary renormalization functions are calculated non-perturbatively in the RI$'$ scheme, and are converted to the $\overline{\rm MS}$ scheme at a scale of 2 GeV. The final values for the momentum fraction are $\langle x \rangle^π_{u^+}=0.261(3)_{\rm stat}(6)_{\rm syst}$, $\langle x \rangle^K_{u^+}=0.246(2)_{\rm stat}(2)_{\rm syst}$, and $\langle x \rangle^K_{s^+}=0.317(2)_{\rm stat}(1)_{\rm syst}$. For the third Mellin moments we find $\langle x^2 \rangle^π_{u^+}=0.082(21)_{\rm stat}(17)_{\rm syst}$, $\langle x^2 \rangle^K_{u^+}=0.093(5)_{\rm stat}(3)_{\rm syst}$, and $\langle x^2 \rangle^K_{s^+}=0.134(5)_{\rm stat}(2)_{\rm syst}$. The reported systematic uncertainties are due to excited-state contamination. We also give the ratio $\langle x^2 \rangle/\langle x \rangle$ which is an indication of how quickly the PDFs lose support at large $x$.

preprint2021arXiv

The scalar, vector and tensor form factors for the pion and kaon from lattice QCD

We present a calculation of the scalar, vector, and tensor form factors for the pion and kaon in lattice QCD. We use an ensemble of two degenerate light, a strange and a charm quark ($N_f=2+1+1$) of maximally twisted mass fermions with clover improvement. The corresponding pion and kaon masses are about 265 MeV and 530 MeV, respectively. The calculation is done in both rest and boosted frames obtaining data for four-vector momentum transfer squared up to $-q^2=2.5$ GeV$^2$ for the pion and 3 GeV$^2$ for the kaon. The excited-states effects are studied by analyzing six values of the source-sink time separation for the rest frame ($1.12-2.23$ fm) and for four values for the boosted frame ($1.12-1.67$ fm). The lattice data are renormalized non-perturbatively and the results for the scheme- and scale-dependent scalar and tensor form factors are presented in the $\overline{\rm MS}$ scheme at a scale of 2 GeV. We apply different parametrizations to describe $q^2$-dependence of the form factors to extract the scalar, vector, and tensor radii, as well as the tensor anomalous magnetic moment. We compare the pion and kaon form factors to study SU(3) flavor symmetry breaking effects. By combining the data for the vector and tensor form factors we also obtain the lowest moment of the densities of transversely polarized quarks in the impact parameter space. Finally, we give an estimate for the average transverse shift in the $y$ direction for polarized quarks in the $x$ direction.

preprint2020arXiv

Comparison of topological charge definitions in Lattice QCD

In this paper, we show a comparison of different definitions of the topological charge on the lattice. We concentrate on one small-volume ensemble with 2 flavours of dynamical, maximally twisted mass fermions and use three more ensembles to analyze the approach to the continuum limit. We investigate several fermionic and gluonic definitions. The former include the index of the overlap Dirac operator, the spectral flow of the Wilson--Dirac operator and the spectral projectors. For the latter, we take into account different discretizations of the topological charge operator and various smoothing schemes to filter out ultraviolet fluctuations: the gradient flow, stout smearing, APE smearing, HYP smearing and cooling. We show that it is possible to perturbatively match different smoothing schemes and provide a well-defined smoothing scale. We relate the smoothing parameters for cooling, stout and APE smearing to the gradient flow time $τ$. In the case of hypercubic smearing the matching is performed numerically. We investigate which conditions have to be met to obtain a valid definition of the topological charge and susceptibility and we argue that all valid definitions are highly correlated and allow good control over topology on the lattice.

preprint2020arXiv

Model-independent determination of the nucleon charge radius from lattice QCD

Lattice QCD calculations of nucleon form factors are restricted to discrete values of the Euclidean four-momentum transfer. Therefore, the extraction of radii typically relies on parametrizing and fitting the lattice QCD data to obtain its slope close to zero momentum transfer. We investigate a new method, which allows to compute the nucleon radius directly from existing lattice QCD data, without assuming a functional form for the momentum dependence of the underlying form factor. The method is illustrated for the case of the isovector mean square charge radius of the nucleon $\langle r^2_\mathrm{isov} \rangle$ and the quark-connected contributions to $\langle r^2_p\rangle$ and $\langle r^2_n \rangle$ for the proton and neutron, respectively. Computations are performed using a single gauge ensemble with $N_f=2+1+1$ maximally twisted mass clover-improved fermions at physical quark mass and a lattice spacing of $a=0.08\mathrm{fm}$.

preprint2020arXiv

Parton distribution functions of $Δ^+$ on the lattice

We perform a first calculation for the unpolarized parton distribution function of the $Δ^+$ baryon using lattice QCD simulations within the framework of Large Momentum Effective Theory. Two ensembles of $N_f=2+1+1$ twisted mass fermions are utilized with a pion mass of 270 MeV and 360 MeV, respectively. The baryon, which is treated as a stable single-particle state, is boosted with momentum $P_3$ with values $\{0.42,0.83,1.25\}$ GeV, and we utilize momentum smearing to improve the signal. The unpolarized parton distribution function of $Δ^+$ is obtained using a non-perturbative renormalization and a one-loop formula for the matching, with encouraging precision. In particular, we compute the $\overline{d}(x)-\overline{u}(x)$ asymmetry and compare it with the same quantity in the nucleon, in a first attempt towards resolving the physical mechanism responsible for generating such asymmetry.

preprint2020arXiv

Tetraquark interpolating fields in a lattice QCD investigation of the $D_{s0}^\ast(2317)$ meson

We investigate the $D_{s0}^\ast(2317)$ meson using lattice QCD and considering correlation functions of several $\bar{c} s$ two-quark and $\bar{c} s (\bar{u} u + \bar{d} d)$ four-quark interpolating fields. These interpolating fields generate different structures in color, spin and position space including quark-antiquark pairs, tetraquarks and two-meson scattering states. For our computation we use an ensemble simulated with pion mass $m_π\approx 0.296 \, \textrm{GeV}$ and spatial volume of extent $2.90 \, \textrm{fm}$. We find in addition to the expected spectrum of two-meson scattering states another state around $60 \, \textrm{MeV}$ below the $D K$ threshold, which we interpret as the $D_{s0}^\ast(2317)$ meson. This state couples predominantly to a quark-antiquark interpolating field and only weakly to a $D K$ two-meson interpolating field. The coupling to the tetraquark interpolating fields is essentially zero, rendering a tetraquark interpretation of the $D_{s0}^\ast(2317)$ meson rather unlikely. Moreover, we perform a scattering analysis using Lüscher's method and the effective range approximation to determine the $D_{s0}^\ast(2317)$ mass for infinite spatial volume. We find this mass $51 \, \textrm{MeV}$ below the $D K$ threshold, rather close to both our finite volume result and the experimentally observed value.

preprint2020arXiv

Unpolarized and helicity generalized parton distributions of the proton within lattice QCD

We present the first calculation of the $x$-dependence of the proton generalized parton distributions (GPDs) within lattice QCD. Results are obtained for the isovector unpolarized and helicity GPDs. We compute the appropriate matrix elements of fast-moving protons coupled to non-local operators containing a Wilson line. We present results for proton momenta $0.83,\,1.25,\,1.67$ GeV, and momentum transfer squared $0.69,\,1.38$ GeV$^2$. These combinations include cases with zero and nonzero skewness. The calculation is performed using one ensemble of two degenerate mass light, a strange and a charm quark of maximally twisted mass fermions with a clover term. The lattice results are matched to the light-cone GPDs using one-loop perturbation theory within the framework of large momentum effective theory. The final GPDs are given in the $\overline{\rm MS}$ scheme at a scale of 2 GeV.

preprint2016arXiv

A Lattice Calculation of Parton Distributions

We present results for the $x$ dependence of the unpolarized, helicity, and transversity isovector quark distributions in the proton using lattice QCD, employing the method of quasi-distributions proposed by Ji in 2013. Compared to a previous calculation by us, the errors are reduced by a factor of about 2.5. Moreover, we present our first results for the polarized sector of the proton, which indicate an asymmetry in the proton sea in favor of the $u$ antiquarks for the case of helicity distributions, and an asymmetry in favor of the $d$ antiquarks for the case of transversity distributions.

preprint2016arXiv

A study of the radiative transition $ππ\to πγ^{*}$ with lattice QCD

Lattice QCD calculations of radiative transitions between hadrons have in the past been limited to processes of hadrons stable under the strong interaction. Recently developed methods for $1\to2$ transition matrix elements in a finite volume now enable the determination of radiative decay rates of strongly unstable particles. Our lattice QCD study focuses on the process $ππ\to πγ^{*}$, where the $ρ$ meson is present as an enhancement in the cross-section. We use $2+1$ flavors of clover fermions at a pion mass of approximately $320$ MeV and a lattice size of approximately $3.6$ fm. The required $2$-point and $3$-point correlation functions are constructed from a set of forward, sequential and stochastic light quark propagators. In addition to determining the $ρ$ meson resonance parameters via the Lüscher method, the scattering phase shift is used in conjunction with the $1\to2$ transition matrix element formalism of Briceño, Hansen and Walker-Loud to compute the $ππ\toπγ^{*}$ amplitude at several values of the momentum transfer and $ππ$ invariant mass.

preprint2016arXiv

Adaptive Aggregation-based Domain Decomposition Multigrid for Twisted Mass Fermions

The Adaptive Aggregation-based Domain Decomposition Multigrid method (arXiv:1303.1377) is extended for two degenerate flavors of twisted mass fermions. By fine-tuning the parameters we achieve a speed-up of the order of hundred times compared to the conjugate gradient algorithm for the physical value of the pion mass. A thorough analysis of the aggregation parameters is presented, which provides a novel insight into multigrid methods for lattice QCD independently of the fermion discretization.

preprint2016arXiv

Disconnected diagrams with twisted-mass fermions

The latest results from the Twisted-Mass collaboration on disconnected diagrams at the physical value of the pion mass are presented. In particular, we focus on the sigma terms, the axial charges and the momentum fraction, all of them for the nucleon. A detailed error analysis for each observable follows, showing the strengths and weaknesses of the one-end trick. Alternatives are discussed.

preprint2016arXiv

Importance of closed quark loops for lattice QCD studies of tetraquarks

To investigate the light scalar tetraquark candidate $a_0(980)$ (quantum numbers $J^P =0^+$), a correlation matrix including a variety of two- and four-quark interpolating operators has to be computed. We discuss efficient techniques to compute the elements of this correlation matrix, in particular diagrams with closed quark loops. Furthermore, we present evidence that such diagrams are not negligible given our precision, and their contribution is essential to obtain physically meaningful results. In particular, we find indications of the existence of an "additional" state around the two-particle thresholds of $K+ \bar K$ and $η+π$, which could correspond to the $a_0(980)$ meson.

preprint2016arXiv

Position space method for the nucleon magnetic moment in lattice QCD

The extraction of the magnetic form factor of the nucleon at zero momentum transfer is usually performed by adopting a parametrization for its momentum dependence and fitting the results obtained at finite momenta. We present a position space method that allows us to remove the momentum prefactor in the form factor decomposition and hence compute the magnetic form factor directly at zero momentum without the need to assume a functional form for its momentum dependence. The method is explored on one ensemble using $N_f=2+1+1$ Wilson twisted mass fermions with a light quark mass corresponding to $M_π\approx373\mathrm{GeV}$ and a lattice spacing of $a\approx0.082\mathrm{fm}$. We obtain results for the isovector magnetic moment and for the proton and neutron magnetic moments. The value we find for the isovector magnetic moment is larger as compared to fitting the form factor at the discrete values of the lattice momentum transfers using a dipole Ansatz, bringing it closer to the experimental value.

preprint2016arXiv

Recent results for the proton spin decomposition from lattice QCD

The exact decomposition of the proton spin has been a much debated topic, on the experimental as well as the theoretical side. In this talk we would like to report on recent non-perturbative results and ongoing efforts to explore the proton spin from lattice QCD. We present results for the relevant generalized form factors from gauge field ensembles that feature a physical value of the pion mass. These generalized form factors can be used to determine the total spin and angular momentum carried by the quarks. In addition we present first results for our ongoing effort to compute the angular momentum of the gluons in the proton.

preprint2016arXiv

Study of decuplet baryon resonances from lattice QCD

A lattice QCD study of the strong decay width and coupling constant of decuplet baryons to an octet baryon - pion state is presented. The transfer matrix method is used to obtain the overlap of lattice states with decuplet baryon quantum numbers on the one hand and octet baryon-pion quantum numbers on the other as an approximation to the matrix element of the corresponding transition. By making use of leading order effective field theory, the coupling constants, as well as the widths for the various decay channels are determined. The transitions studied are $ Δ\to π\,N$, $Σ^* \to Λ\,π$, $Σ^* \to Σ\,π$ and $Ξ^* \to Ξ\,π$. We obtain results for two ensembles of $N_f=2+1$ dynamical fermion configurations, one using domain wall valence quarks on a staggered sea at a pion mass of $350\,\mathrm{MeV}$ and a box size of $3.4\,\mathrm{fm}$ and a second one using domain wall sea and valence quarks at pion mass $180\,\mathrm{MeV}$ and box size $4.5\,\mathrm{fm}$.

preprint2015arXiv

A Lattice Calculation of Parton Distributions

We report on our exploratory study for the direct evaluation of the parton distribution functions from lattice QCD, based on a recently proposed new approach. We present encouraging results using Nf = 2 + 1 + 1 twisted mass fermions with a pion mass of about 370 MeV. The focus of this work is a detailed description of the computation, including the lattice calculation, the matching to an infinite momentum and the nucleon mass correction. In addition, we test the effect of gauge link smearing in the operator to estimate the influence of the Wilson line renormalization, which is yet to be done.

preprint2015arXiv

Calculation of the decay width of decuplet baryons

We calculate the coupling constant and decay width of the decuplet to octet baryon transitions in lattice QCD using the transfer matrix method. The transition amplitude is related to the coupling constant and via the Fermi's Golden Rule to the decay width. The method is applicable for near-degeneracy of the energy levels of initial and final states and, when this condition is fulfilled, yields a good estimate of the decay width. We present results using a hybrid action with domain wall valence quarks on a staggered sea with $350$ MeV pion mass as well as for a domain wall fermion action with $180$ MeV pion mass. We find $Γ\left( Δ\to π\,N \right) = 119\,( 8)\,( 8)$ MeV for the transition of Delta to pion-nucleon within the unitary domain wall setup. We also report values for the decay widths of the $Σ^*$ and $Ξ*$ baryons.

preprint2015arXiv

Computation of correlation matrices for tetraquark candidates with $J^P = 0^+$ and flavor structure $q_1 \bar{q_2} q_3 \bar{q}_3$

The conjecture that several recently observed mesons have a structure, which is not dominated by an ordinary quark-antiquark pair, but by a four-quark structure, is being actively investigated both theoretical and experimentally. Such a state may be characterized as a mesonic molecule or as a diquark-antidiquark pair. Lattice QCD provides a theoretically sound framework to study such states. To quantitatively investigate the internal structure of such mesons, one needs to precisely compute correlation matrices containing several interpolating operators including two and four quarks. Here we discuss certain technical aspects of such correlation matrices suited to study tetraquark candidates with $J^P = 0^+$ and flavor structure $q_1 \bar{q_2} q_3 \bar{q}_3$, e.g.\ the $a_0(980)$ meson, the $D_{s0}^\ast$ meson and some of the charged $c \bar{c}$ $X$ states. Some numerical results for the $a_0(980)$ meson are presented.

preprint2015arXiv

Disconnected quark loop contributions to nucleon observables using $N_f=2$ twisted clover fermions at the physical value of the light quark mass

We compute the disconnected quark loops contributions entering the determination of nucleon observables, by using a $N_f = 2$ ensemble of twisted mass fermions with a clover term at a pion mass $m_π= 133$ MeV. We employ exact deflation and implement all calculations in GPUs, enabling us to achieve large statistics and a good signal.

preprint2015arXiv

Extraction of the isovector magnetic form factor of the nucleon at zero momentum

The extraction of the magnetic form factor of the nucleon at zero momentum transfer is usually performed by adopting a parametrization for its momentum dependence and fitting the results obtained at finite momenta. We present position space methods that rely on taking the derivative of relevant correlators to extract directly the magnetic form factor at zero momentum without the need to assume a functional form for its momentum dependence. These methods are explored on one ensemble using $N_f=2+1+1$ Wilson twisted mass fermions.

preprint2015arXiv

Nucleon electromagnetic form factors from twisted mass lattice QCD

Results on the electromagnetic form factors of the nucleon using twisted mass fermion configurations are presented. These include a gauge field ensemble simulated with two degenerate light quarks yielding a pion mass of around 130 MeV, as well as two ensembles that include strange and charm quarks in the sea yielding pion masses of 210 MeV and 373 MeV. Details of the methods used and systematic errors are discussed, such as noise reduction techniques and the effect of excited states contamination.

preprint2015arXiv

Strangeness of the nucleon from Lattice Quantum Chromodynamics

We present a non-perturbative calculation of the strangeness of the nucleon $y_N$ within the framework of lattice QCD. This observable is known to be an important cornerstone to interpret results from direct dark matter detection experiments. We perform a lattice computation for $y_N$ with an analysis of systematic effects originating from discretization, finite size, chiral extrapolation and excited state effects leading to the value of $y_N= 0.173(50)$. The rather large uncertainty of this value of $y_N$ is dominated by systematic uncertainties which we were able to quantify in this work.

preprint2015arXiv

The electric dipole moment of the neutron from $N_f=2+1+1$ twisted mass fermions

We extract the neutron electric dipole moment (nEDM) $ \vert \vec{d}_n \vert $ on configurations produced with $N_f=2+1+1$ twisted mass fermions with lattice spacing of $a \simeq 0.082$ fm and a light quark mass that corresponds to $M_π \simeq 373$ MeV. We do so by evaluating the $CP$-odd form factor $F_3$ for small values of the $CP$-violation parameter $θ$ in the limit of zero momentum transfer. This limit is extracted using the usual parametrization but in addition position space methods. The topological charge is computed via cooling and gradient flow using the Wilson, Symanzik tree-level improved and Iwasaki actions for smoothing. We obtain consistent results for all choices of smoothing procedures and methods to extract $F_3$ at zero momentum transfer. For the ensemble analyzed we find a value of nEDM of $\vert \vec{d}_n \vert / θ= -0.045(6)(1) {\rm e} \cdot {\rm fm}$.

preprint2014arXiv

A Stochastic Method for Computing Hadronic Matrix Elements

We present a stochastic method for the calculation of baryon three-point functions that is more versatile compared to the typically used sequential method. We analyze the scaling of the error of the stochastically evaluated three-point function with the lattice volume and find a favorable signal-to-noise ratio suggesting that our stochastic method can be used efficiently at large volumes to compute hadronic matrix elements.

preprint2014arXiv

Computation of disconnected contributions to nucleon observables

We compare several methods for computing disconnected fermion loops contributing to nucleon three-point functions. The comparison is carried out using one ensemble of $N_f=2+1+1$ twisted mass fermions with pion mass of 373 MeV. The complete set of operators up to one-derivative are examined by developing optimized code for mutli-GPUs. Simple guidelines are given as to the preferable method for each class of operators.

preprint2014arXiv

First results with twisted mass fermions towards the computation of parton distribution functions on the lattice

We report on our exploratory study for the evaluation of the parton distribution functions from lattice QCD, based on a new method proposed in Ref.~arXiv:1305.1539. Using the example of the nucleon, we compare two different methods to compute the matrix elements needed, and investigate the application of gauge link smearing. We also present first results from a large production ensemble and discuss the future challenges related to this method.

preprint2014arXiv

Investigation of the tetraquark candidate $a_0(980)$: technical aspects and preliminary results

We discuss technical aspects and first results of a lattice QCD study of the $a_0(980)$ state. We employ various interpolating operators of quark-antiquark, mesonic molecule, diquark-antidiquark and two-meson type. Both connected and disconnected contributions including diagrams with closed fermion loops are computed. To keep statistical errors small, it is essential to optimize the computation of these diagrams by choosing that combination of techniques most appropriate for each type of diagram from the correlation matrix of interpolating operators. We illustrate, how this can be done, by discussing certain diagrams in detail. We also present preliminary results corresponding to a $4\times 4$ submatrix computed with 2+1 flavors of clover fermions.

preprint2014arXiv

Nucleon observables and axial charges of other baryons using twisted mass fermions

We present results on the nucleon scalar, axial and tensor charges, as well as, on the first moments of the unpolarized, polarized and transversity parton distributions using $N_f=2$ and $N_f=2+1+1$ twisted mass fermions. These include an ensemble that yields the physical value of the ratio of the nucleon to the pion mass. Results on the axial charges of hyperons and charmed baryons are also presented for a range of pion masses including the physical one.

preprint2014arXiv

Nucleon structure from lattice QCD - recent achievements and perspectives

We present recent developments in lattice QCD simulations as applied in the study of hadron structure. We discuss the challenges and perspectives in the evaluation of benchmark quantities such as the nucleon axial charge and the isovector parton momentum fraction, as well as, in the computation of the nucleon $σ$-terms, which involve the calculation of disconnected quark loop contributions.

preprint2014arXiv

Results on the disconnected contributions for hadron structure

We present results on the disconnected contributions to three point functions entering in studies of hadron structure. We use $N_F = 2+1+1$ twisted mass fermions and give a detailed description on the results of the nucleon σ-terms, isoscalar axial charge and first moments of bare parton distributions for a range of pions masses. In addition we give the σ-terms and the computations are performed using QUDA code implemented on GPUs.

preprint2013arXiv

Implementation of the twisted mass fermion operator in the QUDA library

We discuss an extension of the QUDA library for the Wilson twisted mass operator. A performance analysis is presented for both degenerate and non-degenerate flavor doublets. The degenerate twisted mass fermion operator runs at up to 190, 487 and 856 Gflops, for double, single and half precisions respectively on recent NVIDIA Kepler GPUs, while our implementation for the non-degenerate flavor doublet allows to reach 163, 516 and 879 GFlops, respectively. The code is currently in production for the hadron structure study.

preprint2013arXiv

Looking at the gluon moment of the nucleon with dynamical twisted mass fermions

To understand the structure of hadrons it is important to know the PDF of their constituents, the quarks and gluons. In our work we aim to compute the first moment of the gluon PDF $\langle x \rangle_g$ for the nucleon. We follow two possible approaches in order to extract the gluon moment: the Feynman-Hellmann theorem and a direct method with smearing of the gluon operator. We present preliminary results computed on $24^3 \times 48$ lattices for the case where the Feynman-Hellman theorem is used and $32^3 \times 64$ lattices for the direct method, employing $N_f=2+1+1$ maximally twisted mass fermions.

preprint2013arXiv

Nucleon transversity generalized form factors with twisted mass fermions

We present results on the nucleon tensor form factors and first moment of the transversity distribution using maximally twisted mass fermions. We analyze two Nf=2+1+1 ensembles having pion masses of 213 MeV and 373 MeV with lattice spacing a=0.064 fm and a=0.082 fm, respectively. First results using an Nf=2 ensemble of twisted mass fermions with a clover term at a physical pion mass and lattice spacing a=0.094 fm are also presented. The renormalization function for the local tensor form factors is evaluated non-perturbatively with a perturbative subtraction of O(a^2)-terms, while for the first moment of the transversity we use a perturbative estimate. Results are given in the MS-bar scheme at a scale of mu=2 GeV, and are compared with recent results obtained using different discretization schemes.

preprint2013arXiv

Scalar mesons and tetraquarks from twisted mass lattice QCD

We study light scalar mesons with particular focus on the a_0(980) using lattice QCD with 2+1+1 dynamical quark flavors. To investigate the structure of these scalar mesons and to identify, whether a sizeable tetraquark component is present, we use a large set of operators, including diquark-antidiquark, mesonic molecule and two-meson operators. We find that the low-lying states overlap essentially exclusively with two-meson states. This indicates that in the channels investigated no tightly bound four-quark states of either molecular or diquark-antidiquark type exist.

preprint2012arXiv

Lattice investigation of the tetraquark candidates a0(980) and kappa

It is a long discussed issue whether light scalar mesons have sizeable four-quark components. We present an exploratory study of this question using N_f=2+1+1 twisted mass lattice QCD. A mixed action approach ignoring disconnected contributions is used to calculate correlator matrices consisting of mesonic molecule, diquark-antidiquark and two-meson interpolating operators with quantum numbers of the scalar mesons a0(980) (1(0^++)) and kappa (1/2(0+)). The correlation matrices are analyzed by solving the generalized eigenvalue problem. The theoretically expected free two-particle scattering states are identified, while no additional low lying states are observed. We do not observe indications for bound four-quark states in the channels investigated.

preprint2012arXiv

Scalar mesons and tetraquarks by means of lattice QCD

We study the light scalar mesons a_0(980) and kappa using N_f = 2+1+1 flavor lattice QCD. In order to probe the internal structure of these scalar mesons, and in particular to identify, whether a sizeable tetraquark component is present, we use a large set of operators, including diquark-antidiquark, mesonic molecule and two-meson operators. The inclusion of disconnected diagrams, which are technically rather challenging, but which would allow us to extend our work to e.g. the f_0(980) meson, is introduced and discussed.

preprint2011arXiv

Excited State Effects in Nucleon Matrix Element Calculations

We perform a high-statistics precision calculation of nucleon matrix elements using an open sink method allowing us to explore a wide range of sink-source time separations. In this way the influence of excited states of nucleon matrix elements can be studied. As particular examples we present results for the nucleon axial charge $g_A$ and for the first moment of the isovector unpolarized parton distribution $<x>_{u-d}$. In addition, we report on preliminary results using the generalized eigenvalue method for nucleon matrix elements. All calculations are performed using $N_f=2+1+1$ maximally twisted mass Wilson fermions.

preprint2011arXiv

Nucleon matrix elements with $N_f=2+1+1$ maximally twisted fermions

We present the first lattice calculation of nucleon matrix elements using four dynamical flavors. We use the $N_f=2+1+1$ maximally twisted mass formulation. The renormalization is performed non-perturbatively in the RI$^\prime$-MOM scheme and results are given for the vector and axial vector operators with up to one-derivative. Our calculation of the average momentum of the unpolarized non-singlet parton distribution is presented and compared to our previous results obtained from the $N_f=2$ case.

preprint2011arXiv

Precision Study of Excited State Effects in Nucleon Matrix Elements

We present a dedicated analysis of the influence of excited states on the calculation of nucleon matrix elements. This calculation is performed at a fixed value of the lattice spacing, volume and pion mass that are typical of contemporary lattice computations. We focus on the nucleon axial charge, g_A, for which we use about 7,500 measurements, and on the average momentum of the unpolarized isovector parton distribution, <x>_{u-d}, for which we use about 23,000 measurements. All computations are done employing N_f=2+1+1 maximally-twisted-mass Wilson fermions and using non-perturbatively calculated renormalization factors. Excited state effects are shown to be negligible for g_A whereas they lead to an O(10%) downward shift for <x>_{u-d}.

preprint2010arXiv

Low lying baryon spectrum with $N_f = 2+1+1$ dynamical twisted quarks

We present first results on the octet and decuplet strange baryon spectrum with $N_f=2+1+1$ twisted mass quarks. We use an Osterwalder Seiler valence strange quark with a mass tuned to the kaon and compare the results with those obtained in the unitary setup. This comparison allows to perform a first study of the lattice artefacts introduced by the mixed action approach. We investigate the effect of the strange and charm quarks in the sea by using two lattice spacings and comparing with preceding $N_f = 2$ twisted mass fermion calculations.

preprint2010arXiv

Strange and charmed baryons using N_f=2 twisted mass QCD

We compute the mass spectrum for strange/charmed baryons in the partially quenched approach using N_f=2 twisted mass QCD configurations. We investigate two main issues: the size of lattice artefacts using three values of the lattice spacing (the smallest of which is approximately 0.05 fm) and the dependence of baryon masses on meson (or quark) masses. We thus perform a global fit in order to extrapolate simultaneously to the continuum limit and to the physical point. We estimate the masses of Omega_{sss}, Xi_{dss}, Lambda_{uds}, Omega_{ccc}, Xi_{dcc}, Lambda_{udc}.

preprint2009arXiv

Quark transverse charge densities in the $Δ(1232)$ from lattice QCD

We extend the formalism relating electromagnetic form factors to transverse quark charge densities in the light-front frame to the case of a spin-3/2 baryon and calculate these transverse densities for the $Δ(1232)$ isobar using lattice QCD. The transverse charge densities for a transversely polarized spin-3/2 particle are characterized by monopole, dipole, quadrupole, and octupole patterns representing the structure beyond that of a pure point-like spin-3/2 particle. We present lattice QCD results for the $Δ$-isobar electromagnetic form factors for pion masses down to approximatively 350 MeV for three cases: quenched QCD, two-degenerate flavors of dynamical Wilson quarks, and three flavors of quarks using a mixed action that combines domain wall valence quarks and dynamical staggered sea quarks. We extract transverse quark charge densities from these lattice results and find that the $Δ$ is prolately deformed, as indicated by the fact that the quadrupole moment $G_{E2}(0$) is larger than the value -3 characterizing a point particle and the fact that the transverse charge density in a $Δ^+$ of maximal transverse spin projection is elongated along the axis of the spin.