Source author record

M. Constantinou

M. Constantinou 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

44works
8topics
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

44 published item(s)

preprint2025arXiv

Three-dimensional imaging of hadrons with hard exclusive reactions: advances in experiment, theory, phenomenology, and lattice QCD

Generalized Parton Distributions (GPDs) have emerged as a powerful framework for exploring the internal structure of hadrons in terms of their partonic constituents. Over the past three decades, the field has witnessed significant theoretical and experimental advancements. The interpretation of GPDs in impact parameter space offers a vivid three-dimensional visualization of hadron structure, correlating longitudinal momentum and transverse spatial distributions, thereby enabling tomographic imaging of hadrons. Furthermore, the link between GPDs and the matrix elements of the QCD energy-momentum tensor provides access to fundamental properties of hadrons, including spin decomposition and internal pressure distributions. Notably, recent analyses of Deeply Virtual Compton Scattering (DVCS) data have enabled the empirical extraction of the quark pressure profile inside the proton. This white paper presents an overview of recent developments in GPD theory and phenomenology, as well as progress in lattice QCD studies. It outlines the prospects for advancing our understanding of hadron structure through the next generation of dedicated experiments, including the extension of the Jefferson Lab 12~GeV program (and its potential 22~GeV upgrade), J-PARC, COMPASS/AMBER, LHC ultra-peripheral collisions, and the future electron-ion colliders EIC and EicC.

preprint2021arXiv

$\langle x\rangle$ and $\langle x^2\rangle$ of the pion PDF from Lattice QCD with $N_f=2+1+1$ dynamical quark flavours

Using Nf=2+1+1 lattice QCD, we determine the fermionic connected contributions to the second and third moment of the pion PDF. Based on gauge configurations from the European Twisted Mass Collaboration, chiral and continuum extrapolations are performed using pion masses in the range of 230 to 500 MeV and three values of the lattice spacing. Finite volume effects are investigated using different volumes. In order to avoid mixing under renormalisation for the third moment, we use an operator with two non-zero spatial components of momentum. Momenta are injected using twisted boundary conditions. Our final values read $\langle x\rangle=0.2075(106)$ and $\langle x^2\rangle=0.163(33)$, determined at 2 GeV in the $\overline{MS}$-scheme and with systematic and statistical uncertainties summend in quadrature.

preprint2021arXiv

Charge radii of the nucleon from its flavor dependent Dirac form factors

We have determined the proton and the neutron charge radii from a global analysis of the proton and the neutron elastic form factors, after first performing a flavor decomposition of these form factors under charge symmetry in the light cone frame formulation. We then extracted the transverse mean-square radii of the flavor dependent quark distributions. In turn, these are related in a model-independent way to the proton and neutron charge radii but allow us to take into account motion effects of the recoiling nucleon for data at finite but high momentum transfer. In the proton case we find $\langle r_p \rangle = 0.852 \pm0.002_{\rm (stat.)} \pm0.009_{\rm (syst.)}~({\rm fm})$, consistent with the proton charge radius obtained from muonic hydrogen spectroscopy \cite{pohl:2010,antog2013}. The current method improves on the precision of the $\langle r_p \rangle$ extraction based on the form factor measurements. Furthermore, we find no discrepancy in the $\langle r_p \rangle$ determination among the different electron scattering measurements, all of which, utilizing the current method of extraction, result in a value that is consistent with the smallest $\langle r_p \rangle$ extraction from the electron scattering measurements \cite{Xiong:2019umf}. Concerning the neutron case, past results relied solely on the neutron-electron scattering length measurements, which suffer from an underestimation of underlying systematic uncertainties inherent to the extraction technique. Utilizing the present method we have performed the first extraction of the neutron charge radius based on nucleon form factor data, and we find $\langle r_n^2 \rangle = -0.122 \pm0.004_{\rm (stat.)} \pm0.010_{\rm (syst.)}~({\rm fm}^2)$.

preprint2020arXiv

Complete flavor decomposition of the spin and momentum fraction of the proton using lattice QCD simulations at physical pion mass

We evaluate the gluon and quark contributions to the spin of the proton using an ensemble of gauge configuration generated at physical pion mass. We compute all valence and sea quark contributions to high accuracy. We perform a non-perturbative renormalization for both quark and gluon matrix elements. We find that the contribution of the up, down, strange and charm quarks to the proton intrinsic spin is $\frac{1}{2}\sum_{q=u,d,s,c}ΔΣ^{q^+}=0.191(15)$ and to the total spin $\sum_{q=u,d,s,c}J^{q^+}=0.285(45)$. The gluon contribution to the spin is $J^g=0.187(46)$ yielding $J=J^q+J^g=0.473(71)$ confirming the spin sum. The momentum fraction carried by quarks in the proton is found to be $0.618(60)$ and by gluons $0.427(92)$, the sum of which gives $1.045(118)$ confirming the momentum sum rule. All scale and scheme dependent quantities are given in the $\mathrm{ \overline{MS}}$ scheme at 2 GeV.

preprint2020arXiv

Confronting lattice parton distributions with global QCD analysis

We present the first Monte Carlo based global QCD analysis of spin-averaged and spin-dependent parton distribution functions (PDFs) that includes nucleon isovector matrix elements in coordinate space from lattice QCD. We investigate the degree of universality of the extracted PDFs when the lattice and experimental data are treated under the same conditions within the Bayesian likelihood analysis. For the unpolarized sector, we find rather weak constraints from the current lattice data on the phenomenological PDFs, and difficulties in describing the lattice matrix elements at large spatial distances. In contrast, for the polarized PDFs we find good agreement between experiment and lattice data, with the latter providing significant constraints on the spin-dependent isovector quark and antiquark distributions.

preprint2020arXiv

Moments of nucleon generalized parton distributions from lattice QCD simulations at physical pion mass

We present results for the moments of nucleon isovector vector and axial generalised parton distribution functions computed within lattice QCD. Three ensembles of maximally twisted mass clover-improved fermions simulated with a physical value of the pion mass are analyzed. Two of these ensembles are generated using two degenerate light quarks. A third ensemble is used having, in addition to the light quarks, strange and charm quarks in the sea. A careful analysis of the convergence to the ground state is carried out that is shown to be essential for extracting the correct nucleon matrix elements. This allows a controlled determination of the unpolarised, helicity and tensor second Mellin moments. The vector and axial-vector generalised form factors are also computed as a function of the momentum transfer square up to about 1 GeV$^2$. The three ensembles allow us to check for unquenching effects and to assess lattice finite volume effects.

preprint2020arXiv

Nucleon axial and pseudoscalar form factors from lattice QCD at the physical point

We compute the nucleon axial and induced pseudoscalar form factors using three ensembles of gauge configurations, generated with dynamical light quarks with mass tuned to approximately their physical value. One of the ensembles also includes the strange and charm quarks with their mass close to physical. The latter ensemble has large statistics and finer lattice spacing and it is used to obtain final results, while the other two are used for assessing volume effects. The pseudoscalar form factor is also computed using these ensembles. We examine the momentum dependence of these form factors as well as relations based on pion pole dominance and the partially conserved axial-vector current hypothesis.

preprint2019arXiv

Nucleon strange electromagnetic form factors

The role of the strange quarks on the low-energy interactions of the proton can be probed through the strange electromagnetic form factors. Knowledge of these form factors provides essential input for parity-violating processes and contributes to the understanding of the sea quark dynamics. We determine the strange electromagnetic form factors of the nucleon within the lattice formulation of Quantum Chromodynamics using simulations that include light, strange and charm quarks in the sea all tuned to their physical mass values. We employ state-of-the-art techniques to accurately extract the form factors for values of the momentum transfer square up to 0.8~GeV$^2$. We find that both the electric and magnetic form factors are statistically non-zero. We obtain for the strange magnetic moment $μ^s=-0.017(4)$, the strange magnetic radius $\langle r^2_M \rangle^s=-0.015(9)$~fm$^2$, and the strange charge radius $\langle r^2_E \rangle^s=-0.0048(6)$~fm$^2$.

preprint2016arXiv

Direct Evaluation of the Quark Content of Nucleons from Lattice QCD at the Physical Point

We evaluate the light, strange and charm scalar content of the nucleon using one lattice QCD ensemble generated with two degenerate light quarks with mass fixed to their physical value. We use improved techniques to evaluate the disconnected quark loops to sufficient accuracy to determine the strange and charm nucleon $σ$-terms. in addition to the light quark content $σ_{πN}$. We find $σ_{πN}$=$37.2(2.6)(^{4.7}_{2.9})$ MeV $σ_{s}$=$41.1(8.2)(^{7.8}_{5.8})$ MeV and $σ_c$=$79(21)(^{12}_{\phantom{1}8})$ MeV, where the first error is statistical and the second is the systematic error due to the determination of the lattice spacing, the assessment of finite volume and residual excited state effects.

preprint2016arXiv

Neutron electric dipole moment using $N_f{=}2{+}1{+}1$ twisted mass fermions

We evaluate the neutron electric dipole moment $\vert \vec{d}_N\vert$ using lattice QCD techniques. The gauge configurations analyzed are produced by the European Twisted Mass Collaboration using $N_f{=}2{+}1{+}1$ twisted mass fermions at one value of the lattice spacing of $a \simeq 0.082 \ {\rm fm}$ and a light quark mass corresponding to $m_π \simeq 373 \ {\rm MeV}$. Our approach to extract the neutron electric dipole moment is based on the calculation of the $CP$-odd electromagnetic form factor $F_3(Q^2)$ for small values of the vacuum angle $θ$ in the limit of zero Euclidean momentum transfer $Q^2$. The limit $Q^2 \to 0$ is realized either by adopting a parameterization of the momentum dependence of $F_3(Q^2)$ and performing a fit, or by employing new position space methods, which involve the elimination of the kinematical momentum factor in front of $F_3(Q^2)$. The computation in the presence of a $CP$-violating term requires the evaluation of the topological charge ${\cal Q}$. This is computed by applying the cooling technique and the gradient flow with three different actions, namely the Wilson, the Symanzik tree-level improved and the Iwasaki action. We demonstrate that cooling and gradient flow give equivalent results for the neutron electric dipole moment. Our analysis yields a value of $\vert \vec{d}_N\vert=0.045(6)(1)\ \barθ \ e \cdot {\rm fm}$ for the ensemble with $m_π=373$ MeV considered.

preprint2016arXiv

Nucleon and pion structure with lattice QCD simulations at physical value of the pion mass

We present results on the nucleon scalar, axial and tensor charges as well as on the momentum fraction, and the helicity and transversity moments. The pion momentum fraction is also presented. The computation of these key observables is carried out using lattice QCD simulations at a physical value of the pion mass. The evaluation is based on gauge configurations generated with two degenerate sea quarks of twisted mass fermions with a clover term. We investigate excited states contributions with the nucleon quantum numbers by analyzing three sink-source time separations. We find that, for the scalar charge, excited states contribute significantly and to a less degree to the nucleon momentum fraction and helicity moment. Our analysis yields a value for the nucleon axial charge agrees with the experimental value and we predict a value of 1.027(62) in the $\overline{\text{MS}}$ scheme at 2 GeV for the isovector nucleon tensor charge directly at the physical point. The pion momentum fraction is found to be $\langle x\rangle_{u-d}^{π^\pm}=0.214(15)(^{+12}_{-9})$ in the $\overline{\rm MS}$ at 2 GeV.

preprint2016arXiv

Nucleon spin and quark content at the physical point

We present results on the spin and quark content of the nucleon using $N_f=2$ twisted mass clover-improved fermion simulations with a pion mass close to its physical value. We use recently developed methods to obtain accurate results for both connected and disconnected contributions. We provide results for the axial charge, quark and gluon momentum fraction as well as the light, strange and charm $σ$-terms.

preprint2015arXiv

Electrical analysis of hysteresis in solution processed silicon nanowire field effect transistors

Silicon nanowires (Si NW) are ideal candidates for solution processable field effect transistors (FETs). The interface between the nanowire channel and the gate dielectric plays a crucial role in the FET performance, and it can be responsible for unwanted effects such as hysteresis of the I-V characteristics due to threshold voltage shift when the gate voltage is applied. Using gate-voltage bias stress measurements we show that a large hysteresis of up to 40V in Si NW FETs with SiO2 dielectric is mainly due to the holes traps at the nanowire/SiO2 interface. An approach for reducing this hysteresis to just 2.5V using solution processable fluoropolymer dielectric Cytop in the top-gate configuration is demonstrated. Experimental results suggest that the density of surface traps in Si nanowire transistors is dictated mainly by the nature of the dielectric layer. The influence of the gate dielectric was studied by assessing the field effect transport behaviour of a representative double gate FETs based on SiO2 bottom dielectric and top Cytop dielectric layer. Such devices were characterised, revealing an order of magnitude higher hole traps density at the nanowire/SiO2 interface (1x10^13cm-2) compared to that of nanowire/fluoropolymer interface (7.5x10^11cm-2).

preprint2015arXiv

First moment of the flavour octet nucleon parton distribution function using lattice QCD

We perform a lattice computation of the flavour octet contribution to the average quark momentum in a nucleon, $\la x\ra^{(8)}_{μ^2 = 4 \gev^2}$. In particular, we fully take the disconnected contributions into account in our analysis for which we use a generalization of the technique developed in \cite{Dinter:2012tt}. We investigate systematic effects with a particular emphasis on the excited states contamination. We find that in the renormalization free ratio $\frac{\la x \ra^{(3)}}{\la x \ra^{(8)}}$ (with $\la x \ra^{(3)}$ the non-singlet moment) the excited state contributions cancel to a large extend making this ratio a promising candidate for a comparison to phenomenological analyses. Our final result for this ratio is in agreement with the phenomenological value and we find, including systematic errors, $\frac{\la x \ra^{(3)}}{\la x \ra^{(8)}} = 0.39(1)(4)$.

preprint2014arXiv

Disconnected quark loop contributions to nucleon observables in lattice QCD

We perform a high statistics calculation of disconnected fermion loops on Graphics Processing Units for a range of nucleon matrix elements extracted using lattice QCD. The isoscalar electromagnetic and axial vector form factors, the sigma-terms and the momentum fraction and helicity are among the quantities we evaluate. We compare the disconnected contributions to the connected ones and give the physical implications on nucleon observables that probe its structure.

preprint2014arXiv

Evaluation of disconnected quark loops for hadron structure using GPUs

A number of stochastic methods developed for the calculation of fermion loops are investigated and compared, in particular with respect to their efficiency when implemented on Graphics Processing Units (GPUs). We assess the performance of the various methods by studying the convergence and statistical accuracy obtained for observables that require a large number of stochastic noise vectors, such as the isoscalar nucleon axial charge. The various methods are also examined for the evaluation of sigma-terms where noise reduction techniques specific to the twisted mass formulation can be utilized thus reducing the required number of stochastic noise vectors.

preprint2014arXiv

K -> pi matrix elements of the chromagnetic operator on the lattice

We present preliminary results of the first lattice QCD calculation of the K -> pi matrix elements of the chromomagnetic operator O_{CM}=g sbar sigma_{munu} G_{munu} d, which appears in the effective Hamiltonian describing Delta S=1 transitions in and beyond the Standard Model. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with operators of equal and lower dimensionality. The multiplicative renormalization factor as well as the mixing coefficients with the operators of equal dimension have been computed at one-loop in perturbation theory. The power divergent coefficients controlling the mixing with operators of lower dimension have been computed non-perturbatively, by imposing suitable subtraction conditions. The numerical simulations have been carried out using the gauge field configurations produced by the European Twisted Mass Collaboration with N_f=2+1+1 dynamical quarks at three values of the lattice spacing. Our preliminary result for the B-parameter of the chromomagnetic operator is B_{CMO}=0.29(11), which can be compared with the estimate B_{CMO}~1-4 currently used in phenomenological analyses.

preprint2014arXiv

Perturbative and non-perturbative renormalization results of the Chromomagnetic Operator on the Lattice

The Chromomagnetic operator (CMO) mixes with a large number of operators under renormalization. We identify which operators can mix with the CMO, at the quantum level. Even in dimensional regularization (DR), which has the simplest mixing pattern, the CMO mixes with a total of 9 other operators, forming a basis of dimension-five, Lorentz scalar operators with the same flavor content as the CMO. Among them, there are also gauge noninvariant operators; these are BRST invariant and vanish by the equations of motion, as required by renormalization theory. On the other hand using a lattice regularization further operators with $d \leq 5$ will mix; choosing the lattice action in a manner as to preserve certain discrete symmetries, a minimul set of 3 additional operators (all with $d<5$) will appear. In order to compute all relevant mixing coefficients, we calculate the quark-antiquark (2-pt) and the quark-antiquark-gluon (3-pt) Green's functions of the CMO at nonzero quark masses. These calculations were performed in the continuum (dimensional regularization) and on the lattice using the maximally twisted mass fermion action and the Symanzik improved gluon action. In parallel, non-perturbative measurements of the $K-π$ matrix element are being performed in simulations with 4 dynamical ($N_f = 2+1+1$) twisted mass fermions and the Iwasaki improved gluon action.

preprint2014arXiv

Renormalization of Flavor Singlet and Nonsinglet Fermion Bilinear Operators

We compute the difference in the renormalization of flavor singlet and nonsinglet fermion bilinear operators, to two loops in perturbation theory. Our results are applicable to a rather wide class of lattice actions with Symanzik improved gluons, stout links and clover fermions, including the Twisted Mass and SLiNC actions. A more detailed presentation of our results, along with relevant references, will appear in a forthcoming publication.

preprint2014arXiv

Renormalization of local quark-bilinear operators for Nf=3 flavors of SLiNC fermions

The renormalization factors of local quark-bilinear operators are computed non-perturbatively for $N_f=3$ flavors of SLiNC fermions, with emphasis on the various procedures for the chiral and continuum extrapolations. The simulations are performed at a lattice spacing $a=0.074$ fm, and for five values of the pion mass in the range of 290-465 MeV, allowing a safe and stable chiral extrapolation. Emphasis is given in the subtraction of the well-known pion pole which affects the renormalization factor of the pseudoscalar current. We also compute the inverse propagator and the Green's functions of the local bilinears to one loop in perturbation theory. We investigate lattice artifacts by computing them perturbatively to second order as well as to all orders in the lattice spacing. The renormalization conditions are defined in the RI$'$-MOM scheme, for both the perturbative and non-perturbative results. The renormalization factors, obtained at different values of the renormalization scale, are translated to the ${\bar{\rm MS}}$ scheme and are evolved perturbatively to 2 GeV. Any residual dependence on the initial renormalization scale is eliminated by an extrapolation to the continuum limit. We also study the various sources of systematic errors. Particular care is taken in correcting the non-perturbative estimates by subtracting lattice artifacts computed to one loop perturbation theory using the same action. We test two different methods, by subtracting either the ${\cal O}(g^2\,a^2)$ contributions, or the complete (all orders in $a$) one-loop lattice artifacts.

preprint2013arXiv

Nucleon form factors and moments of generalized parton distributions using $N_f=2+1+1$ twisted mass fermions

We present results on the axial and the electromagnetic form factors of the nucleon, as well as, on the first moments of the nucleon generalized parton distributions using maximally twisted mass fermions. We analyze two N_f=2+1+1 ensembles having pion masses of 210 MeV and 354 MeV at two values of the lattice spacing. The lattice scale is determined using the nucleon mass computed on a total of 18 N_f=2+1+1 ensembles generated at three values of the lattice spacing, $a$. The renormalization constants are evaluated non-perturbatively with a perturbative subtraction of ${\cal O}(a^2)$-terms. The moments of the generalized parton distributions are given in the $\bar{\rm MS}$ scheme at a scale of $ μ=2$ GeV. We compare with recent results obtained using different discretization schemes. The implications on the spin content of the nucleon are also discussed.

preprint2013arXiv

Nucleon generalized form factors with twisted mass fermions

We present results on the nucleon form factors, momentum fraction and helicity moment for $N_f=2$ and $N_f=2+1+1$ twisted mass fermions for a number of lattice volumes and lattice spacings. First results for a new $N_f=2$ ensemble at the physical pion mass are also included. The implications of these results on the spin content of the nucleon are discussed taking into account the disconnected contributions at one pion mass.

preprint2013arXiv

Perturbative renormalization functions of local operators for staggered fermions with stout improvement

In this paper we present the perturbative computation of the renormalization functions for the quark field and for a complete set of ultra-local fermion bilinears. The computation of the relevant Green's functions was carried out at 1-loop level for the staggered action using massive fermions. The gluon links which appear both in the fermion action and in the definition of the bilinears have been improved by applying a stout smearing procedure up to 2 times, iteratively. In the gluon sector we employed the Symanzik improved gauge action for different sets of values of the Symanzik coefficients. The renormalization functions are presented in (two variants of) the RI' and in the MSbar renormalization scheme; the dependence on all stout parameters, as well as on the fermion mass, the gauge fixing parameter and the renormalization scale, is shown explicitly. This work is related to our recent paper [Phys. Rev. D86 (2012) 094512, arXiv:1209.6015]. To make our results easily accessible to the reader, we include them in the distribution package of this paper, as a Mathematica input file, Staggered.m.

preprint2013arXiv

Perturbative renormalization of staggered fermion operators with stout improvement: Application to the magnetic susceptibility of QCD

We calculate the fermion propagator and the quark-antiquark Green's functions for a complete set of ultralocal fermion bilinears, ${{\cal O}_Γ}$ [$Γ$: scalar (S), pseudoscalar (P), vector (V), axial (A) and tensor (T)], using perturbation theory up to one-loop and to lowest order in the lattice spacing. We employ the staggered action for fermions and the Symanzik Improved action for gluons. From our calculations we determine the renormalization functions for the quark field and for all ultralocal taste-singlet bilinear operators. The novel aspect of our calculations is that the gluon links which appear both in the fermion action and in the definition of the bilinears have been improved by applying a stout smearing procedure up to two times, iteratively. Compared to most other improved formulations of staggered fermions, the above action, as well as the HISQ action, lead to smaller taste violating effects. The renormalization functions are presented in the RI$'$ scheme; the dependence on all stout parameters, as well as on the coupling constant, the number of colors, the lattice spacing, the gauge fixing parameter and the renormalization scale, is shown explicitly. We apply our results to a nonperturbative study of the magnetic susceptibility of QCD at zero and finite temperature. In particular, we evaluate the "tensor coefficient", $τ$, which is relevant to the anomalous magnetic moment of the muon.

preprint2013arXiv

Perturbatively improving renormalization constants

Renormalization factors relate the observables obtained on the lattice to their measured counterparts in the continuum in a suitable renormalization scheme. They have to be computed very precisely which requires a careful treatment of lattice artifacts. In this work we present a method to suppress these artifacts by subtracting one-loop contributions proportional to the square of the lattice spacing calculated in lattice perturbation theory.

preprint2013arXiv

Perturbatively improving RI-MOM renormalization constants

The determination of renormalization factors is of crucial importance in lattice QCD. They relate the observables obtained on the lattice to their measured counterparts in the continuum in a suitable renormalization scheme. Therefore, they have to be computed as precisely as possible. A widely used approach is the nonperturbative Rome-Southampton method. It requires, however, a careful treatment of lattice artifacts. In this paper we investigate a method to suppress these artifacts by subtracting one-loop contributions to renormalization factors calculated in lattice perturbation theory. We compare results obtained from a complete one-loop subtraction with those calculated for a subtraction of contributions proportional to the square of the lattice spacing.

preprint2013arXiv

The chromomagnetic operator on the lattice

We study matrix elements of the "chromomagnetic" operator on the lattice. This operator is contained in the strangeness-changing effective Hamiltonian which describes electroweak effects in the Standard Model and beyond. Having dimension 5, the chromomagnetic operator is characterized by a rich pattern of mixing with other operators of equal and lower dimensionality, including also non gauge invariant quantities; it is thus quite a challenge to extract from lattice simulations a clear signal for the hadronic matrix elements of this operator. We compute all relevant mixing coefficients to one loop in lattice perturbation theory; this necessitates calculating both 2-point (quark-antiquark) and 3-point (gluon-quark-antiquark) Green's functions at nonzero quark masses. We use the twisted mass lattice formulation, with Symanzik improved gluon action. For a comprehensive presentation of our results, along with detailed explanations and a more complete list of references, we refer to our forthcoming publication [1].

preprint2013arXiv

The quark contents of the nucleon and their implication for dark matter search

We present results concerning the light and strange quark contents of the nucleon using $N_f=2+1+1$ flavours of maximally twisted mass fermions. The corresponding $σ$-terms are casting light on the origin of the nucleon mass and their values are important to interpret experimental data from direct dark matter searches. We discuss our strategy to estimate systematic uncertainties arising in our computations. Our preliminary results for the $σ-$terms read $σ_{πN} = 37(2.6)(24.7) \mev$ and $σ_s=28(8)(10) \mev$. We present our recent final analysis of the $y_N$ parameter and found $y_N=0.135(46)$ including systematics\cite{Alexandrou:2013nda}.

preprint2013arXiv

Thermodynamic properties of QCD in external magnetic fields

We consider the effect of strong external electromagnetic fields on thermodynamic observables in QCD, through lattice simulations with 1+1+1 flavors of staggered quarks at physical quark masses. Continuum extrapolated results are presented for the light quark condensates and for their tensor polarizations, as functions of the temperature and the magnetic field. We find the light condensates to undergo inverse magnetic catalysis in the transition region, in a manner that the transition temperature decreases with growing magnetic field. We also compare the results to other approaches and lattice simulations. Furthermore, we relate the tensor polarization to the spin part of the magnetic susceptibility of the QCD vacuum, and show that this contribution is diamagnetic.

preprint2012arXiv

Magnetic susceptibility of QCD at zero and at finite temperature from the lattice

The response of the QCD vacuum to a constant external (electro)magnetic field is studied through the tensor polarization of the chiral condensate and the magnetic susceptibility at zero and at finite temperature. We determine these quantities using lattice configurations generated with the tree-level Symanzik improved gauge action and N_f=1+1+1 flavors of stout smeared staggered quarks with physical masses. We carry out the renormalization of the observables under study and perform the continuum limit both at T>0 and at T=0, using different lattice spacings. Finite size effects are studied by using various spatial lattice volumes. The magnetic susceptibilities χ_f reveal a diamagnetic behavior; we obtain at zero temperature χ_u=-(2.08 +/- 0.08) 1/GeV^2, χ_d=-(2.02 +/- 0.09) 1/GeV^2 and χ_s=-(3.4 +/- 1.4) 1/GeV^2 for the up, down and strange quarks, respectively, in the MSBar scheme at a renormalization scale of 2 GeV. We also find the polarization to change smoothly with the temperature in the confinement phase and then to drastically reduce around the transition region.

preprint2012arXiv

Perturbative subtraction of lattice artifacts in the computation of renormalization constants

The determination of renormalization factors is of crucial importance. They relate the observables obtained on finite, discrete lattices to their measured counterparts in the continuum in a suitable renormalization scheme. Therefore, they have to be computed as precisely as possible. A widely used approach is the nonperturbative Rome-Southampton method. It requires, however, a careful treatment of lattice artifacts. They are always present because simulations are done at lattice spacings $a$ and momenta $p$ with $ap$ not necessarily small. In this paper we try to suppress these artifacts by subtraction of one-loop contributions in lattice perturbation theory. We compare results obtained from a complete one-loop subtraction with those calculated for a subtraction of $O(a^2)$.

preprint2012arXiv

Renormalization constants of local operators for Wilson type improved fermions

Perturbative and non-perturbative results are presented on the renormalization constants of the quark field and the vector, axial-vector, pseudoscalar, scalar and tensor currents. The perturbative computation, carried out at one-loop level and up to second order in the lattice spacing, is performed for a fermion action, which includes the clover term and the twisted mass parameter yielding results that are applicable for unimproved Wilson fermions, as well as for improved clover and twisted mass fermions. We consider ten variants of the Symanzik improved gauge action corresponding to ten different values of the plaquette coefficients. Non-perturbative results are obtained using the twisted mass Wilson fermion formulation employing two degenerate dynamical quarks and the tree-level Symanzik improved gluon action. The simulations are performed for pion masses in the range of 480 MeV to 260 MeV and at three values of the lattice spacing, a, corresponding to beta=3.9, 4.05, 4.20. For each renormalization factor computed non-perturbatively we subtract its perturbative O(a^2) terms so that we eliminate part of the cut-off artifacts. The renormalization constants are converted to MS-bar at a scale of mu=2 GeV. The perturbative results depend on a large number of parameters and are made easily accessible to the reader by including them in the distribution package of this paper, as a Mathematica input file.

preprint2012arXiv

Sigma terms and strangeness content of the nucleon with $Nf=2+1+1$ twisted mass fermions

We investigate excited state contaminations in a direct computation of the nucleon $σ$-terms. This is an important source of systematic effects that needs to be controlled besides the light quark mass dependence and lattice artefacts. We use maximally twisted mass fermions with dynamical light ($u$,$d$), strange and charm degrees of freedom. Employing an efficient stochastic evaluation of the disconnected contribution available for twisted mass fermions, we show that the effect of excited states is large in particular for the strange $σ$-terms, where it can be as big as $O(\gtrsim 40%$). This leads to the unfortunate conclusion that even with a source-sink separation of $\sim 1.5 \fm$ and a good statistical accuracy it is not clear, whether excited state effects are under control for this quantity.

preprint2011arXiv

BK-parameter from Nf = 2 twisted mass lattice QCD

We present an unquenched $N_f=2$ lattice computation of the $B_{K}$ parameter which controls $K^0-\bar K^0$ oscillations. A partially quenched setup is employed with two maximally twisted dynamical (sea) light Wilson quarks, and valence quarks of both the maximally twisted and the Osterwalder--Seiler variety. Suitable combinations of these two kinds of valence quarks lead to a lattice definition of the $B_{K}$ parameter which is both multiplicatively renormalizable and O($a$) improved. Employing the non-perturbative RI-MOM scheme, in the continuum limit and at the physical value of the pion mass we get $B^{\rm RGI}_K=0.729\pm 0.030$, a number well in line with the existing quenched and unquenched determinations.

preprint2011arXiv

Moments of nucleon generalized parton distributions from lattice QCD

We present results on the lower moments of the nucleon generalized parton distri butions within lattice QCD using two dynamical flavors of degenerate twisted mass fermions. Our simulations are performed on lattices with three different values of the lattice spacings, namely $a=0.089$ fm, $a=0.070$ fm and $a=0.056$ fm, allowing the investigation of cut-off effects. The volume dependence is examined using simulations on two lattices of spatial length $L=2.1$ fm and $L=2.8$ fm. The simulations span pion masses in the range of 260-470 MeV. Our results are renormalized non-perturbatively and the values are given in the $\bar{\rm MS}$ scheme at a scale $ μ=2$ GeV. They are chirally extrapolated to the physical point in order to compare with experiment. The consequences of these results on the spin carried by the quarks in the nucleon are investigated.

preprint2011arXiv

Nucleon electromagnetic form factors in twisted mass lattice QCD

We present results on the nucleon electromagnetic form factors within lattice QCD using two flavors of degenerate twisted mass fermions. Volume effects are examined using simulations at two volumes of spatial length L=2.1 fm and L=2.8 fm. Cut-off effects are investigated using three different values of the lattice spacings, namely a=0.089 fm, a=0.070 and a=0.056 fm. The nucleon magnetic moment, Dirac and Pauli radii are obtained in the continuum limit and chirally extrapolated to the physical pion mass allowing for a comparison with experiment.

preprint2011arXiv

Nucleon form factors and moments of parton distributions in twisted mass lattice QCD

We present results on the electroweak form factors and on the lower moments of parton distributions of the nucleon, within lattice QCD using two dynamical flavors of degenerate twisted mass fermions. Results are obtained on lattices with three different values of the lattice spacings, namely a=0.089 fm, a=0.070 fm and a=0.056 fm, allowing the investigation of cut-off effects. The volume dependence is examined by comparing results on two lattices of spatial length L=2.1 fm and L=2.8 fm. The simulations span pion masses in the range of 260-470 MeV. Our results are renormalized non-perturbatively and the values are given in the MS-scheme at a scale mu=2 GeV.

preprint2011arXiv

Renormalization constants for 2-twist operators in twisted mass QCD

Perturbative and non-perturbative results on the renormalization constants of the fermion field and the twist-2 fermion bilinears are presented with emphasis on the non-perturbative evaluation of the one-derivative twist-2 vector and axial vector operators. Non-perturbative results are obtained using the twisted mass Wilson fermion formulation employing two degenerate dynamical quarks and the tree-level Symanzik improved gluon action. The simulations have been performed for pion masses in the range of about 450-260 MeV and at three values of the lattice spacing $a$ corresponding to $β=3.9, 4.05, 4.20$. Subtraction of ${\cal O}(a^2)$ terms is carried out by performing the perturbative evaluation of these operators at 1-loop and up to ${\cal O}(a^2)$. The renormalization conditions are defined in the RI$'$-MOM scheme, for both perturbative and non-perturbative results. The renormalization factors, obtained for different values of the renormalization scale, are evolved perturbatively to a reference scale set by the inverse of the lattice spacing. In addition, they are translated to ${\bar{\rm MS}}$ at 2 GeV using 3-loop perturbative results for the conversion factors.

preprint2010arXiv

Axial Nucleon form factors from lattice QCD

We present results on the nucleon axial form factors within lattice QCD using two flavors of degenerate twisted mass fermions. Volume effects are examined using simulations at two volumes of spatial length $L=2.1$ fm and $L=2.8$ fm. Cut-off effects are investigated using three different values of the lattice spacings, namely $a=0.089$ fm, $a=0.070$ fm and $a=0.056$ fm. The nucleon axial charge is obtained in the continuum limit and chirally extrapolated to the physical pion mass enabling comparison with experiment.

preprint2010arXiv

Non-perturbative renormalization of quark bilinear operators with Nf=2 (tmQCD) Wilson fermions and the tree-level improved gauge action

We present results for the renormalization constants of bilinear quark operators obtained by using the tree-level Symanzik improved gauge action and the Nf=2 twisted mass fermion action at maximal twist, which guarantees automatic O(a)-improvement. Our results are also relevant for the corresponding standard (un-twisted) Wilson fermionic action since the two actions only differ, in the massless limit, by a chiral rotation of the quark fields. The scale-independent renormalization constants ZV, ZA and the ratio ZP/ZS have been computed using the RI-MOM approach, as well as other alternative methods. For ZA and ZP/ZS, the latter are based on both standard twisted mass and Osterwalder-Seiler fermions, while for ZV a Ward Identity has been used. The quark field renormalization constant Zq and the scale dependent renormalization constants ZS, ZP and ZT are determined in the RI-MOM scheme. Leading discretization effects of O(g^2 a^2), evaluated in one-loop perturbation theory, are explicitly subtracted from the RI-MOM estimates.

preprint2010arXiv

Renormalization constants for one-derivative fermion operators in twisted mass QCD

We present perturbative and non-perturbative results on the renormalization constants of the local and one-derivative vector and axial vector operators. Non-perturbative results are obtained using the twistedmassWilson fermion formulation employing two degenerate dynamical quarks and the tree-level Symanzik improved gluon action for pion masses in the range of about 450-260 MeV and at there values of the lattice spacing, namely 0.055 fm, 0.070 fm and 0.089 fm. Subtraction of O(a^2) terms is carried out by performing the perturbative evaluation of these operators at 1- loop and up to O(a^2). The renormalization conditions are defined in the RI'-MOM scheme, for both perturbative and non-perturbative results. The Z-factors, obtained for different values of the renormalization scale, are evolved perturbatively to a reference scale set by the inverse of the lattice spacing. In addition, they are translated to MS-bar at 2 GeV using 3-loop perturbative results for the conversion factors.

preprint2009arXiv

O(a^2) corrections to the one-loop propagator and bilinears of clover fermions with Symanzik improved gluons

We calculate corrections to the fermion propagator and to the Green's functions of all fermion bilinear operators of the form $\barΨΓΨ$, to one-loop in perturbation theory. We employ the Wilson/clover action for fermions and a family of Symanzik improved actions for gluons. The novel aspect of our calculations is that they are carried out to second order in the lattice spacing, $O(a^2)$. Consequently, they have addressed a number of new issues, most notably the appearance of loop integrands with strong IR divergences (convergent only beyond 6 dimensions). Such integrands are not present in $O(a^1)$ improvement calculations; there, IR divergent terms are seen to have the same structure as in the $O(a^0)$ case, by virtue of parity under integration, and they can thus be handled by well-known techniques. We explain how to correctly extract the full $O(a^2)$ dependence; in fact, our method is generalizable to any order in $a$. The $O(a^2)$ corrections to the quark propagator and Green's functions computed in this paper are useful to improve the nonperturbative RI-MOM determination of renormalization constants for quark bilinear operators. Our results depend on a large number of parameters: coupling constant, number of colors, lattice spacing, external momentum, clover parameter, Symanzik coefficients, gauge parameter. To make these results most easily accessible to the reader, we have included them in the distribution package of this paper, as an ASCII file named: Oa2results.m ; the file is best perused as Mathematica input.