Source author record

Krzysztof Cichy

Krzysztof Cichy 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

49works
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

49 published item(s)

preprint2023arXiv

Generalized Parton Distributions from Lattice QCD with Asymmetric Momentum Transfer: Unpolarized Quarks

Traditionally, lattice QCD computations of generalized parton distributions (GPDs) have been carried out in a symmetric frame, where the transferred momentum is symmetrically distributed between the incoming and outgoing hadrons. However, such frames are inconvenient since they require a separate calculation for each value of the momentum transfer, increasing significantly the computational cost. In this work, by focusing on the quasi-distribution approach, we lay the foundation for faster and more effective lattice QCD calculations of GPDs exploiting asymmetric frames, with freedom in the transferred momentum distribution. An important ingredient of our approach is the Lorentz covariant parameterization of the matrix elements in terms of Lorentz-invariant amplitudes, which allows one to relate matrix elements in different frames. We also use this amplitude approach to propose a new definition of quasi-GPDs that is frame-independent and, more importantly, may lead to smaller power corrections in the matching relations to the light-cone GPDs. We demonstrate the efficacy of the formalism through numerical calculations using one ensemble of $N_f$=2+1+1 twisted mass fermions with a clover improvement. The value of the light-quark masses lead to a pion mass of about 260 MeV. Concentrating on the proton, and limiting ourselves to a vanishing longitudinal momentum transfer to the target, we extract the invariant amplitudes from matrix element calculations in both the symmetric and asymmetric frame, and obtain results for the twist-2 light-cone GPDs for unpolarized quarks, that is, $H$ and $E$.

preprint2023arXiv

GPDs in asymmetric frames

It is often taken for granted that Generalized Parton Distributions (GPDs) are defined in the "symmetric" frame, where the transferred momentum is symmetrically distributed between the incoming/outgoing hadrons. However, such frames pose computational challenges for the lattice QCD practitioners. In these proceedings, we lay the foundation for lattice QCD calculations of GPDs in "asymmetric" frames, where the transferred momentum is not symmetrically distributed between the incoming/outgoing hadrons. The novelty of our work relies on the parameterization of the matrix elements in terms of Lorentz-invariant amplitudes, which not only helps in establishing relations between the said frames but also helps in isolating higher-twist contaminations. As an example, we focus on the unpolarized GPDs for spin-1/2 particles.

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

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

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.

preprint2021arXiv

Flavor non-singlet parton distribution functions from lattice QCD at physical quark masses via the pseudo-distribution approach

One of the great challenges of QCD is to determine the partonic structure of the nucleon from first principles. In this work, we provide such a determination of the flavor non-singlet ($u-d$) unpolarized parton distribution function (PDF), utilizing the non-perturbative formulation of QCD on the lattice. We apply Radyushkin's pseudo-distribution approach to lattice results obtained using simulations with the light quark mass fixed to its physical value; this is the first ever attempt for this approach directly at the physical point. The extracted coordinate-space matrix elements are used to find the relevant physical Ioffe time distributions from a matching procedure. The full Bjorken-$x$ dependence of PDFs is resolved using several reconstruction methods to tackle the ill-conditioned inverse problem encountered when using discrete lattice data. We consider both the valence distribution $q_v$ and the combination with antiquarks $q_v+2\bar{q}$, related to, respectively, the real and imaginary part of extracted matrix elements. Good agreement is found with PDFs from global fits already within statistical uncertainties and it is further improved by quantifying several systematic effects. The results presented here are the first ever \emph{ab initio} determinations of PDFs fully consistent with global fits in the whole $x$-range. Thus, they pave the way to investigating a wider class of partonic distributions, such as e.g.\ singlet PDFs and generalized parton distributions. Therefore, essential and yet missing first-principle insights can be achieved, complementing the rich experimental programs dedicated to the structure of the nucleon.

preprint2021arXiv

Overview of lattice calculations of the $x$-dependence of PDFs, GPDs and TMDs

For a long time, lattice QCD was unable to address the $x$-dependence of partonic distributions, direct access to which is impossible in Euclidean spacetime. Recent years have brought a breakthrough for such calculations when it was realized that partonic light-cone correlations can be accessed through spatial correlations computable on the lattice. Appropriately devised observables can be factorized into physical PDFs via a perturbative procedure called matching, analogous to the standard factorization of experimental cross sections. In this short review, aimed at a broader high-energy and nuclear physics community, we discuss the recent highlights of this research program. Key concepts are outlined, followed by a case study illustrating the typical stage of current lattice extractions and by a brief review of the most recent explorations. We finalize with a number of messages for the prospects of lattice determinations of partonic structure.

preprint2021arXiv

Parton distribution functions beyond leading twist from lattice QCD: The $h_L(x)$ case

We report the first-ever calculation of the isovector flavor combination of the chiral-odd twist-3 parton distribution $h_L(x)$ for the proton from lattice QCD. We employ gauge configurations with two degenerate light, a strange and a charm quark ($N_f=2+1+1$) of maximally twisted mass fermions with a clover improvement. The lattice has a spatial extent of 3 fm and lattice spacing of 0.093 fm. The values of the quark masses lead to a pion mass of $260$ MeV. We use a source-sink time separation of 1.12 fm to control contamination from excited states. Our calculation is based on the quasi-distribution approach, with three values for the proton momentum: 0.83 GeV, 1.25 GeV, and 1.67 GeV. The lattice data are renormalized non-perturbatively using the RI$'$ scheme, and the final result for $h_L(x)$ is presented in the $\overline{\rm MS}$ scheme at the scale of 2 GeV. Furthermore, we compute in the same setup the transversity distribution, $h_1(x)$, which allows us, in particular, to compare $h_L(x)$ to its Wandzura-Wilczek approximation. We also combine results for the isovector and isoscalar flavor combinations to disentangle the individual quark contributions for $h_1(x)$ and $h_L(x)$, and address the Wandzura-Wilczek approximation in that case as well.

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

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

Running coupling constant from position-space current-current correlation functions in three-flavor lattice QCD

In this Letter, we provide a determination of the coupling constant in three-flavor quantum chromodynamics (QCD), $α^{\overline{\mathrm{MS}}}_s(μ)$, for $\overline{\mathrm{MS}}$ renormalization scales $μ\in (1,\,2)$ GeV. The computation uses gauge field configuration ensembles with $\mathcal{O}(a)$-improved Wilson-clover fermions generated by the Coordinated Lattice Simulations (CLS) consortium. Our approach is based on current-current correlation functions and has never been applied before in this context. We convert the results perturbatively to the QCD $Λ$-parameter and obtain $Λ_{\overline{\mathrm{MS}}}^{N_f=3} = 342 \pm 17$ MeV, which agrees with the world average published by the Particle Data Group and has competing precision. The latter was made possible by a unique combination of state-of-the-art CLS ensembles with very fine lattice spacings, further reduction of discretization effects from a dedicated numerical stochastic perturbation theory simulation, combining data from vector and axial-vector channels and matching to high-order perturbation theory.

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.

preprint2019arXiv

Parton distributions from lattice data: the nonsinglet case

We revise the relation between Parton Distribution Functions (PDFs) and matrix elements computable from lattice QCD, focusing on the quasi-Parton Distribution Functions (qPDFs) approach. We exploit the relation between PDFs and qPDFs in the case of the unpolarized isovector parton distribution to obtain a factorization formula relating the real and imaginary part of qPDFs matrix elements to specific nonsinglet distributions, and we propose a general framework to extract PDFs from the available lattice data, treating them on the same footing as experimental data. We implement the proposed approach within the NNPDF framework, and we study the potentiality of such lattice data in constraining PDFs, assuming some plausible scenarios to assess the unknown systematic uncertainties. We finally extract the two nonsinglet distributions involved in our analysis from a selection of the available lattice data.

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

Chiral condensate in the Schwinger model with Matrix Product Operators

Tensor network (TN) methods, in particular the Matrix Product States (MPS) ansatz, have proven to be a useful tool in analyzing the properties of lattice gauge theories. They allow for a very good precision, much better than standard Monte Carlo (MC) techniques for the models that have been studied so far, due to the possibility of reaching much smaller lattice spacings. The real reason for the interest in the TN approach, however, is its ability, shown so far in several condensed matter models, to deal with theories which exhibit the notorious sign problem in MC simulations. This makes it prospective for dealing with the non-zero chemical potential in QCD and other lattice gauge theories, as well as with real-time simulations. In this paper, using matrix product operators, we extend our analysis of the Schwinger model at zero temperature to show the feasibility of this approach also at finite temperature. This is an important step on the way to deal with the sign problem of QCD. We analyze in detail the chiral symmetry breaking in the massless and massive cases and show that the method works very well and gives good control over a broad range of temperatures, essentially from zero to infinite temperature.

preprint2016arXiv

Computing the topological susceptibility from fixed topology QCD simulations

The topological susceptibility is an important quantity in QCD, which can be computed using lattice methods. However, at a fine lattice spacing, or when using high quality chirally symmetric quarks, algorithms which proceed in small update steps --- in particular the HMC algorithm --- tend to get stuck in a single topological sector. In such cases, the computation of the topological susceptibility is not straightforward. Here, we explore two methods to extract the topological susceptibility from lattice QCD simulations restricted to a single topological sector. The first method is based on the correlation function of the topological charge density, while the second method relies on measuring the topological charge within spacetime subvolumes. Numerical results for two-flavor QCD obtained by using both methods are presented.

preprint2016arXiv

Continuum limit of the $D$ meson, $D_s$ meson and charmonium spectrum from $N_f=2+1+1$ twisted mass lattice QCD

We compute masses of $D$ meson, $D_s$ meson and charmonium states using $N_f=2+1+1$ Wilson twisted mass lattice QCD. All results are extrapolated to physical light quark masses, physical strange and charm quark masses and to the continuum. Our analysis includes states with spin $J = 0,1,2$, parity $\mathcal{P} = -,+$ and in case of charmonium also charge conjugation $\mathcal{C} = -,+$. Computations are based on a large set of quark-antiquark meson creation operators. We investigate and quantify all sources of systematic errors, including fitting range uncertainties, finite volume effects, isospin breaking effects and the choice of the fitting ansatz for the combined chiral and continuum extrapolation such that the resulting meson masses can be compared directly and in a meaningful way to experimental results. Within combined statistical and systematic errors, which are between below two per mille and three percent, our results agree with available experimental results for most of the states. In the few cases where we observe discrepancies, we discuss possible reasons.

preprint2016arXiv

Investigation of $B\bar B$ four-quark systems using lattice QCD

We investigate $B \bar B$ systems by computing potentials of two static quarks in the presence of two quarks of finite mass using lattice QCD. By solving the Schrödinger equation we check whether these potentials are sufficiently attractive to host bound states. Particular focus is put on the experimentally most promising bottomonium-like tetraquark candidate $Z_b^\pm$ with quantum numbers $I(J^P)=1(1^+)$.

preprint2016arXiv

Non-perturbative running of renormalization constants from correlators in coordinate space using step scaling

Working in a quenched setup with Wilson twisted mass valence fermions, we explore the possibility to compute non-perturbatively the step scaling function using the coordinate (X-space) renormalization scheme. This scheme has the advantage of being on-shell and gauge invariant. The step scaling method allows us to calculate the running of the renormalization constants of quark bilinear operators. We describe here the details of this calculation. The aim of this exploratory study is to identify the feasibility of the X-space scheme when used in small volume simulations required by the step scaling technique. Eventually, we translate our final results to the continuum MSbar scheme and compare against four-loop analytic formulae finding satisfactory agreement.

preprint2016arXiv

Step scaling in coordinate space: running of the quark mass

We perform a benchmark study of the step scaling procedure for the ratios of renormalization constants extracted from position space correlation functions. We work in the quenched approximation and consider the pseudoscalar, scalar, vector and axial vector bilinears. The pseudoscalar/scalar cases allow us to obtain the non-perturbative running of the quark mass over a wide range of energy scales - from around 17 GeV to below 1.5 GeV - which agrees well with the 4-loop prediction of continuum perturbation theory. We find that step scaling is feasible in X-space and we discuss its advantages and potential problems.

preprint2016arXiv

The multi-flavor Schwinger model with chemical potential - Overcoming the sign problem with Matrix Product States

During recent years there has been an increasing interest in the application of matrix product states, and more generally tensor networks, to lattice gauge theories. This non-perturbative method is sign problem free and has already been successfully used to compute mass spectra, thermal states and phase diagrams, as well as real-time dynamics for Abelian and non-Abelian gauge models. In previous work we showed the suitability of the method to explore the zero-temperature phase structure of the multi-flavor Schwinger model at non-zero chemical potential, a regime where the conventional Monte Carlo approach suffers from the sign problem. Here we extend our numerical study by looking at the spatially resolved chiral condensate in the massless case. We recover spatial oscillations, similar to the theoretical predictions for the single-flavor case, with a chemical potential dependent frequency and an amplitude approximately given by the homogeneous zero density condensate value.

preprint2016arXiv

The Slab Method to Measure the Topological Susceptibility

In simulations of a model with topological sectors, algorithms which proceed in small update steps tend to get stuck in one sector, especially on fine lattices. This distorts the numerical results; in particular it is not straightforward to measure the topological susceptibility chi_t. We test a method to measure chi_t even if configurations from only one sector are available. It is based on the topological charges in sub-volumes, which we denote as "slab". This enables the evaluation of chi_t, as we demonstrate with numerical results for non-linear sigma-models and for 2-flavour QCD. In the latter case, the gradient flow is applied for the smoothing of the gauge configurations, and the slab method results for chi_t are stable over a broad range of flow times.

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

BB interactions with static bottom quarks from Lattice QCD

The isospin, spin and parity dependent potential of a pair of $B$ mesons is computed using Wilson twisted mass lattice QCD with two flavours of degenerate dynamical quarks. The $B$ meson is addressed in the static-light approximation, i.e.\ the $b$ quarks are infinitely heavy. From the results of the $B\,B$ meson-meson potentials, a simple rule can be deduced stating which isospin, spin and parity combinations correspond to attractive and which to repulsive forces. We provide fits to the ground state potentials in the attractive channels and discuss the potentials in the repulsive and excited channels. The attractive channels are most important since they can possibly lead to a bound four-quark state, i.e.\ a $\bar{b}\bar{b}ud$ tetraquark. Using these attractive potentials in the Schrödinger equation, we find indication for such a tetraquark state of two static bottom antiquarks and two light $u/d$ quarks with mass extrapolated down to the physical value.

preprint2015arXiv

Evidence for the existence of $u d \bar{b} \bar{b}$ and the non-existence of $s s \bar{b} \bar{b}$ and $c c \bar{b} \bar{b}$ tetraquarks from lattice QCD

We combine lattice QCD results for the potential of two static antiquarks in the presence of two quarks $q q$ of finite mass and quark model techniques to study possibly existing $q q \bar{b} \bar{b}$ tetraquarks. While there is strong indication for a bound four-quark state for $q q = (ud-du) / \sqrt{2}$, i.e. isospin $I=0$, we find clear evidence against the existence of corresponding tetraquarks with $q q \in \{ uu , (ud+du) / \sqrt{2} , dd \}$, i.e. isospin $I=1$, $q q = s s$ and $q q = c c$.

preprint2015arXiv

Exploring possibly existing $q q \bar b \bar b$ tetraquark states with $q q = ud, ss, cc$

We compute potentials of two static antiquarks in the presence of two quarks $qq$ of finite mass using lattice QCD. In a second step we solve the Schrödinger equation, to determine, whether the resulting potentials are sufficiently attractive to host a bound state, which would indicate the existence of a stable $q q \bar b \bar b$ tetraquark. We find a bound state for $qq=(ud-du)/\sqrt{2}$ with corresponding quantum numbers $I(J^ P)=0(1^+)$ and evidence against the existence of bound states with isospin $I=1$ or $qq \in \{cc,ss \}$

preprint2015arXiv

Mass spectra of mesons containing charm quarks - continuum limit results from twisted mass fermions

We present results from an ongoing computation of masses of $D$ mesons, $D_s$ mesons and charmonium, including both ground states and several parity and angular momentum excitations. We employ 2+1+1 flavours of dynamical maximally twisted mass fermions at three lattice spacings and three $u/d$ quark masses at each lattice spacing. We consider different combinations of valence quark discretizations, with either identical or opposite signs in front of the twisted mass terms. In the end, our setup allows for a good control of different kinds of systematic effects, in particular the quark mass dependence of the resulting meson masses and cut-off effects. We obtain good agreement with experiment for the majority of states and we discuss improvements that will be made to finalize the analysis.

preprint2015arXiv

Non-perturbative Test of the Witten-Veneziano Formula from Lattice QCD

We compute both sides of the Witten-Veneziano formula using lattice techniques. For the one side we perform dedicated quenched simulations and use the spectral projector method to determine the topological susceptibility in the pure Yang-Mills theory. The other side we determine in lattice QCD with $N_f=2+1+1$ dynamical Wilson twisted mass fermions including for the first time also the flavour singlet decay constant. The Witten-Veneziano formula represents a leading order expression in the framework of chiral perturbation theory and we also employ leading order chiral perturbation theory to relate the flavor singlet decay constant to the relevant decay constant parameters in the quark flavor basis and flavor non-singlet decay constants. After taking the continuum and the SU$(2)$ chiral limits we compare both sides and find good agreement within uncertainties.

preprint2015arXiv

Short distance singularities and automatic O($a$) improvement: the cases of the chiral condensate and the topological susceptibility

Short-distance singularities in lattice correlators can modify their Symanzik expansion by leading to additional O($a$) lattice artifacts. At the example of the chiral condensate and the topological susceptibility, we show how to account for these lattice artifacts for Wilson twisted mass fermions and show that the property of automatic O($a$) improvement is preserved at maximal twist.

preprint2015arXiv

Testing the Witten-Veneziano Formula on the Lattice

In this proceeding contribution we report on a test of the famous Witten-Veneziano formula using lattice techniques. We perform dedicated quenched simulations and apply the spectral projector method to determine the topological susceptibility in pure Yang-Mills theory. In order to compute the relevant meson masses and the flavor singlet decay constant we employ lattice QCD with $N_f=2+1+1$ dynamical Wilson twisted mass fermions. Taking the continuum and the $SU(2)$ chiral limits we find good agreement within uncertainties for both sides of the formula.

preprint2015arXiv

The microscopic Twisted Mass Dirac spectrum and the spectral determination of the LECs of Wilson $χ$-PT

We present the comparison of the analytical microscopic spectral density for lattice QCD with $N_{\rm f}=2$ twisted mass fermions with the one obtained on the lattice utilizing configurations produced by the ETM collaboration. We extract estimates for the chiral condensate as well as the low-energy constant $W_8$ of Wilson $χ$-PT by employing spectral information of the Wilson Dirac operator with fixed index at finite volume.

preprint2014arXiv

Comparison of different lattice definitions of the topological charge

We present a comparison of different definitions of the topological charge on the lattice, using a small-volume ensemble with 2 flavours of dynamical twisted mass fermions. The investigated definitions are: index of the overlap Dirac operator, spectral projectors, spectral flow of the Hermitian Wilson-Dirac operator and field theoretic with different kinds of smoothing of gauge fields (HYP and APE smearings, gradient flow, cooling). We also show some results on the topological susceptibility.

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

Lattice Hamiltonian approach to the Schwinger model: further results from the strong coupling expansion

We employ exact diagonalization with strong coupling expansion to the massless and massive Schwinger model. New results are presented for the ground state energy and scalar mass gap in the massless model, which improve the precision to nearly $10^{-9} %$. We also investigate the chiral condensate and compare our calculations to previous results available in the literature. Oscillations of the chiral condensate which are present while increasing the expansion order are also studied and are shown to be directly linked to the presence of flux loops in the system.

preprint2014arXiv

Quark mass anomalous dimension and $Λ_{\overline{\textrm{MS}}}$ from the twisted mass Dirac operator spectrum

We investigate whether it is possible to extract the quark mass anomalous dimension and its scale dependence from the spectrum of the twisted mass Dirac operator in Lattice QCD. The answer to this question appears to be positive, provided that one goes to large enough eigenvalues, sufficiently above the non-perturbative regime. The obtained results are compared to continuum perturbation theory. By analyzing possible sources of systematic effects, we find the domain of applicability of the approach, extending from an energy scale of around 1.5 to 4 GeV. The lower limit is dictated by physics (non-perturbative effects at low energies), while the upper bound is set by the ultraviolet cut-off of present-day lattice simulations. The information about the scale dependence of the anomalous dimension allows also to extract the value of the $Λ_{\overline{\textrm{MS}}}$-parameter of 2-flavour QCD, yielding the value $303(13)(25)$ MeV, where the first error is statistical and the second one systematic. We use gauge field configuration ensembles generated by the European Twisted Mass Collaboration (ETMC) with 2 flavours of dynamical twisted mass quarks, at 4 lattice spacings in the range between around 0.04 and 0.08 fm.

preprint2014arXiv

The temperature dependence of the chiral condensate in the Schwinger model with Matrix Product States

We present our recent results for the tensor network (TN) approach to lattice gauge theories. TN methods provide an efficient approximation for quantum many-body states. We employ TN for one dimensional systems, Matrix Product States, to investigate the 1-flavour Schwinger model. In this study, we compute the chiral condensate at finite temperature. From the continuum extrapolation, we obtain the chiral condensate in the high temperature region consistent with the analytical calculation by Sachs and Wipf.

preprint2014arXiv

Topological susceptibility from the twisted mass Dirac operator spectrum

We present results of our computation of the topological susceptibility with $N_f=2$ and $N_f=2+1+1$ flavours of maximally twisted mass fermions, using the method of spectral projectors. We perform a detailed study of the quark mass dependence and discretization effects. We make an attempt to confront our data with chiral perturbation theory and extract the chiral condensate from the quark mass dependence of the topological susceptibility. We compare the value with the results of our direct computation from the slope of the mode number. We emphasize the role of autocorrelations and the necessity of long Monte Carlo runs to obtain results with good precision. We also show our results for the spectral projector computation of the ratio of renormalization constants $Z_P/Z_S$.

preprint2013arXiv

Chiral condensate from the twisted mass Dirac operator spectrum

We present the results of our computation of the chiral condensate with $N_f=2$ and $N_f=2+1+1$ flavours of maximally twisted mass fermions. The condensate is determined from the Dirac operator spectrum, applying the spectral projector method proposed by Giusti and Luscher. We use 3 lattice spacings and several quark masses at each lattice spacing to perform the chiral and continuum extrapolations. We study the effect of the dynamical strange and charm quarks by comparing our results for $N_f=2$ and $N_f=2+1+1$ dynamical flavours.

preprint2013arXiv

Computation of the chiral condensate using $N_f=2$ and $N_f=2+1+1$ dynamical flavors of twisted mass fermions

We apply the spectral projector method, recently introduced by Giusti and Lüscher, to compute the chiral condensate using $N_f=2$ and $N_f=2+1+1$ dynamical flavors of maximally twisted mass fermions. We present our results for several quark masses at three different lattice spacings which allows us to perform the chiral and continuum extrapolations. In addition we report our analysis on the $O(a)$ improvement of the chiral condensate for twisted mass fermions. We also study the effect of the dynamical strange and charm quarks by comparing our results for $N_f=2$ and $N_f=2+1+1$ dynamical flavors.

preprint2013arXiv

Lattice Hamiltonian approach to the massless Schwinger model: precise extraction of the mass gap

We present results of applying the Hamiltonian approach to the massless Schwinger model. A finite basis is constructed using the strong coupling expansion to a very high order. Using exact diagonalization, the continuum limit can be reliably approached. This allows to reproduce the analytical results for the ground state energy, as well as the vector and scalar mass gaps to an outstanding precision better than 10^{-6} %.

preprint2013arXiv

Matrix Product States for Lattice Field Theories

The term Tensor Network States (TNS) refers to a number of families of states that represent different ansätze for the efficient description of the state of a quantum many-body system. Matrix Product States (MPS) are one particular case of TNS, and have become the most precise tool for the numerical study of one dimensional quantum many-body systems, as the basis of the Density Matrix Renormalization Group method. Lattice Gauge Theories (LGT), in their Hamiltonian version, offer a challenging scenario for these techniques. While the dimensions and sizes of the systems amenable to TNS studies are still far from those achievable by 4-dimensional LGT tools, Tensor Networks can be readily used for problems which more standard techniques, such as Markov chain Monte Carlo simulations, cannot easily tackle. Examples of such problems are the presence of a chemical potential or out-of-equilibrium dynamics. We have explored the performance of Matrix Product States in the case of the Schwinger model, as a widely used testbench for lattice techniques. Using finite-size, open boundary MPS, we are able to determine the low energy states of the model in a fully non-perturbative manner and to accurately extract the mass spectrum. In this work we extend the analysis to the determination of the chiral condensate, both for massless and massive fermions. The method allows for accurate finite size and continuum limit extrapolations and produces remarkably precise results, thus showing the feasibility of these techniques for gauge theory problems.

preprint2013arXiv

Properties of pseudoscalar flavour-singlet mesons from 2+1+1 twisted mass lattice QCD

We study properties of pseudoscalar flavour-singlet mesons from Wilson twisted mass lattice QCD with Nf=2+1+1 dynamical quark flavors. Results for masses are presented at three values of the lattice spacing and light quark masses corresponding to values of the pion mass from 230 MeV to 500 MeV. We briefly discuss scaling effects and the light and strange quark mass dependence of the eta mass. In addition we present an exploratory study using Osterwalder-Seiler type strange and charm valence quarks. This approach avoids some of the complications of the twisted mass heavy doublet. We present first results for matching valence and unitary actions and a comparison of statistical uncertainties.

preprint2013arXiv

Topological susceptibility from twisted mass fermions using spectral projectors

We discuss the computation of the topological susceptibility using the method of spectral projectors and dynamical twisted mass fermions. We present our analysis concerning the O(a)-improvement of the topological susceptibility and we show numerical results for Nf=2 and Nf=2+1+1 flavours, performing a study of the quark mass dependence in terms of leading order chiral perturbation theory.

preprint2012arXiv

Non-perturbative renormalization in coordinate space for $N_f=2$ maximally twisted mass fermions with tree-level Symanzik improved gauge action

We present results of a lattice QCD application of a coordinate space renormalization scheme for the extraction of renormalization constants for flavour non-singlet bilinear quark operators. The method consists in the analysis of the small-distance behaviour of correlation functions in Euclidean space and has several theoretical and practical advantages, in particular: it is gauge invariant, easy to implement and has relatively low computational cost. The values of renormalization constants in the X-space scheme can be converted to the MSbar scheme via 4-loop continuum perturbative formulae. Our results for $N_f=2$ maximally twisted mass fermions with tree-level Symanzik improved gauge action are compared to the ones from the RI-MOM scheme and show full agreement with this method.

preprint2012arXiv

Overlap valence quarks on a twisted mass sea: a case study for mixed action Lattice QCD

We discuss a Lattice QCD mixed action investigation employing Wilson maximally twisted mass sea and overlap valence fermions. Using four values of the lattice spacing, we demonstrate that the overlap Dirac operator assumes a point-like locality in the continuum limit. We also show that by adopting suitable matching conditions for the sea and valence theories a consistent continuum limit for the pion decay constant and light baryon masses can be obtained. Finally, we confront results for sea-valence mixed meson masses and the valence scalar correlator with corresponding expressions of chiral perturbation theory. This allows us to extract low energy constants of mixed action chiral perturbation which characterize the strength of unitarity violations in our mixed action setup.

preprint2011arXiv

$N_f = 2+1+1$ flavours of twisted mass quarks: cut-off effects at tree-level of perturbation theory

We present a calculation of cut-off effects at tree-level of perturbation theory for the K and D mesons using the twisted mass formulation of lattice QCD. The analytical calculations are performed in the time-momentum frame. The relative sizes of cut-off effects are compared for the pion, the kaon and the D meson masses. In addition, different realizations of maximal twist condition are considered and the corresponding cut-off effects are analyzed.

preprint2010arXiv

Continuum Limit of Overlap Valence Quarks on a Twisted Mass Sea

We study a lattice QCD mixed action with overlap valence quarks on two flavours of Wilson maximally twisted mass sea quarks. Employing three different matching conditions to relate both actions to each other, we investigate the continuum limit by using three values of the lattice spacing ranging from 0.05 fm to 0.08 fm. A particular emphasis is put on the effect on physical observables of the topological zero modes appearing in the valence overlap operator. We estimate the region of parameter space where the contribution from these zero modes is sufficiently small such that their effects can be safely controlled and a restoration of unitarity of the mixed action in the continuum limit is reached.

preprint2010arXiv

Overlap Valence Quarks on a Twisted Mass Sea

We present the results of an investigation of a mixed action approach of overlap valence and maximally twisted mass sea quarks. Employing a particular matching condition on the pion mass, we analyze the continuum limit scaling of the pion decay constant and the role of chiral zero modes of the overlap operator in this process. We employ gauge field configurations generated by the European Twisted Mass Collaboration with linear lattice size $L$ ranging from 1.3 to 1.9 fm. The continuum limit is taken at a fixed value of $L=1.3$ fm, employing three values of the lattice spacing and two values of the pion mass constructed from sea quarks only.

preprint2009arXiv

Continuum-limit scaling of overlap fermions as valence quarks

We present the results of a mixed action approach, employing dynamical twisted mass fermions in the sea sector and overlap valence fermions, with the aim of testing the continuum limit scaling behaviour of physical quantities, taking the pion decay constant as an example. To render the computations practical, we impose for this purpose a fixed finite volume with lattice size $L\approx1.3$ fm. We also briefly review the techniques we have used to deal with overlap fermions.