Search

We present a comparison of different quantum state preparation algorithms and their overall efficiency for the Schwinger model with a theta term. While adiabatic state preparation (ASP) is proved to be effective, in practice it leads to large CNOT gate counts to prepare the ground state. The quantum approximate optimization algorithm (QAOA) provides excellent results while keeping the CNOT counts small by design, at the cost of an expensive classical minimization process. We introduce a ``blocked'' modification of the Schwinger Hamiltonian to be used in the QAOA that further decreases the length of the algorithms as the size of the problem is increased. The rodeo algorithm (RA) provides a powerful tool to efficiently prepare any eigenstate of the Hamiltonian, as long as its overlap with the initial guess is large enough. We obtain the best results when combining the blocked QAOA ansatz and the RA, as this provides an excellent initial state with a relatively short algorithm without the need to perform any classical steps for large problem sizes., Comment: 11 pages, 9 figures, 4 tables

We present a new method to calculate moments of parton distribution functions of any order with lattice QCD computations. This method leverages the gradient flow for fermion and gauge fields. The flowed matrix elements of twist-2 operators renormalize multiplicatively, and the matching with physical matrix elements is achieved through the use of continuum symmetries. We derive the matching coefficients at one-loop in perturbation theory for moments of any order in the flavor non-singlet case and provide specific examples of operators suitable for lattice QCD computations. The multiplicative renormalization and matching are independent of the choice of Lorentz indices, allowing the use of temporal indices for twist-2 operators of any dimension. This approach should then also significantly enhance the signal-to-noise ratio in the computation of moments., Comment: 6 pages - Proceedings Contribution to the 31st International Workshop on Deep Inelastic Scattering - DIS2024

Fermionic gradient flow in combination with the short-flow-time expansion provides a computational method where the renormalisation of hadronic matrix elements on the lattice can be simplified to address e.g. the issue that operators with different mass dimension can mix. We demonstrate our gradient flow renormalisation procedure by determining matrix elements of four-quark operators describing neutral meson mixing or meson lifetimes. While meson mixing calculations are well-established on the lattice and serve to validate our procedure, a lattice calculation of matrix elements for heavy meson lifetimes is still outstanding. Preliminary results for mesons formed of a charm and strange quark are presented., Comment: 11 pages, 7 figures, 1 table, Contribution to the 41st International Symposium on Lattice Field Theory, July 28th - August 3rd 2024, Liverpool, UK

We discuss a machine learning (ML) regression model to reduce the computational cost of disconnected diagrams in lattice QCD calculations. This method creates a mapping between the results of fermionic loops computed at different quark masses and flow times. The ML mapping, trained with just a small fraction of the complete data set, makes use of translational invariance and provides consistent result with comparable uncertainties over the calculation done over the whole ensemble, resulting in a significant computational gain., Comment: 17 pages, 11 figures. Version accepted for publication

The continued generation of $N_f=2+1$ quark flavor gauge configurations using stabilized Wilson fermions by the open lattice initiative (OpenLat) is reported. We present the status of our ongoing production and show updates on increasing statistics at the four lattice spacings $a=0.12, 0.094, 0.077$ and $0.064$ fm. Aside from the $SU(3)$ flavor symmetric point we discuss advancements in going towards physical pion masses. We show preliminary results of the pion decay constants, extending previous results, and discuss further validation observables on the available ensembles., Comment: 9 pages, 4 figures, proceedings contribution to "The 40th International Symposium on Lattice Field Theory (Lattice 2023)"

We describe a procedure to determine moments of parton distribution functions of any order in lattice quantum chromodynamics (QCD). The procedure is based on the gradient flow for fermion and gauge fields. The flowed matrix elements of twist-2 operators renormalize multiplicatively, and the matching with the physical matrix elements can be obtained using continuum symmetries and the irreducible representations of Euclidean 4-dimensional rotations. We calculate the matching coefficients at one-loop in perturbation theory for moments of any order in the flavor nonsinglet case. We also give specific examples of operators that could be used in lattice QCD computations. It turns out that it is possible to choose operators with identical Lorentz indices and still have a multiplicative matching. One can thus use twist-2 operators exclusively with temporal indices, thus substantially improving the signal-to-noise ratio in the computation of the hadronic matrix elements., Comment: 10 pages, 2 figures - Version matches published version

Neutral meson mixing and meson lifetimes are theory-side parametrised in terms four-quark operators which can be determined by calculating weak decay matrix elements using lattice Quantum Chromodynamics. While calculations of meson mixing matrix elements are standard, determinations of lifetimes typically suffer from complications in renormalisation procedures because dimension-6 four-quark operators can mix with operators of lower mass dimension and, moreover, quark-line disconnected diagrams contribute. We present work detailing the idea to use fermionic gradient flow to non-perturbatively renormalise matrix elements describing meson mixing or lifetimes, and combining it with a perturbative calculation to match to the $\overline{\rm MS}$ scheme using the shoft-flow-time expansion., Comment: 11 pages, 3 figures, 1 table, Contribution to the 40th International Symposium on Lattice Field Theory, 31st July - 4th August, 2023, Fermilab, USA

The calculation of the neutron electric dipole moment within effective field theories for physics beyond the Standard Model requires non-perturbative hadronic matrix elements of effective operators composed of quark and gluon fields. In order to use input from lattice computations, these matrix elements must be translated from a scheme suitable for lattice QCD to the minimal-subtraction scheme used in the effective-field-theory framework. The accuracy goal in the context of the neutron electric dipole moment necessitates at least a one-loop matching calculation. Here, we provide the one-loop matching coefficients for the $CP$-odd three-gluon operator between two different minimally subtracted 't Hooft-Veltman schemes and the gradient flow. This completes our program to obtain the one-loop gradient-flow matching coefficients for all $CP$-violating and flavor-conserving operators in the low-energy effective field theory up to dimension six., Comment: 6 pages, 4 figures; version published in PLB

Molecules containing short-lived, radioactive nuclei are uniquely positioned to enable a wide range of scientific discoveries in the areas of fundamental symmetries, astrophysics, nuclear structure, and chemistry. Recent advances in the ability to create, cool, and control complex molecules down to the quantum level, along with recent and upcoming advances in radioactive species production at several facilities around the world, create a compelling opportunity to coordinate and combine these efforts to bring precision measurement and control to molecules containing extreme nuclei. In this manuscript, we review the scientific case for studying radioactive molecules, discuss recent atomic, molecular, nuclear, astrophysical, and chemical advances which provide the foundation for their study, describe the facilities where these species are and will be produced, and provide an outlook for the future of this nascent field.

I analyze cutoff effects of the gradient flow for Wilson-type fermions. I show that with a proper choice of the higher dimensional fields in the Symanzik effective theory, O($a$) improvement of the action is achieved changing the initial conditions of the gradient flow equation., Comment: 8 pages. Proceedings contribution to the "39th International Symposium on Lattice Field Theory (Lattice 2022)"

We review the gradient flow for gauge and fermion fields and its applications to lattice gauge theory computations. Using specific examples, we discuss the interplay between perturbative and non-perturbative calculations in the context of renormalization with the gradient flow., Comment: 12 pages, 3 figures, review presented at the XVth "Quark confinement and the Hadron spectrum" conference. This version fixes few typos of the published version

The OpenLat initiative presents its results of lattice QCD simulations using Stabilized Wilson Fermions (SWF) using 2+1 quark flavors. Focusing on the $\mathrm{SU}(3)$ flavor symmetric point $m_\pi=m_K=412$ MeV, four different lattice spacings ($a=0.064,0.077,0.094,0.12$ fm) are used to perform the continuum limit to study cutoff effects. We present results on light hadron masses; for the determination we use a Bayesian analysis framework with constraints and model averaging to minimize the bias in the analysis., Comment: 8 pages, 5 figure, to be published in Proceedings of Science, contribution to "The 39th International Symposium on Lattice Field Theory (LATTICE 2022). Corrected and reviewed version

In this contribution we report the status and plans of the open lattice initiative to generate and share new gauge ensembles using the stabilised Wilson fermion framework. The production strategy is presented in terms of a three stage plan alongside summaries of the data management as well as access policies. Current progress in completing the first stage of generating ensembles at four lattice spacings at the flavor symmetric point is given., Comment: 8 pages, 1 figure, to be published in Proceedings of Science, contribution to "The 39th International Symposium on Lattice Field Theory (LATTICE 2022)", this article supersedes arXiv:2201.03874, corrected and reviewed version

Static electric dipole moments of nondegenerate systems probe mass scales for physics beyond the Standard Model well beyond those reached directly at high energy colliders. Discrimination between different physics models, however, requires complementary searches in atomic-molecular-and-optical, nuclear and particle physics. In this report, we discuss the current status and prospects in the near future for a compelling suite of such experiments, along with developments needed in the encompassing theoretical framework., Comment: Contribution to Snowmass 2021; updated with community edits and endorsements

The gradient flow exponentially suppresses ultraviolet field fluctuations and removes ultraviolet divergences (up to a multiplicative fermionic wavefunction renormalization). It can be used to describe real-space Wilsonian renormalization group transformations and determine the corresponding beta function. We propose a new nonperturbative renormalization scheme for local composite fermionic operators that uses the gradient flow and is amenable to lattice QCD calculations. We present preliminary nonperturbative results for the running of quark bilinear operators in this scheme and outline the calculation of perturbative matching to the MS-bar scheme., Comment: Contribution to The 28th International Symposium on Lattice Field Theory. 12 pages, 4 figures. Minor clarifications for v2

In this joint contribution we announce the formation of the "OPEN LATtice initiative", https://openlat1.gitlab.io, to study Stabilised Wilson Fermions (SWF). They are a new avenue for QCD calculations with Wilson-type fermions and we report results on our continued study of this framework: Tuning the clover improvement coefficient, and extending the reach of lattice spacings to $a=0.12$ fm. We fix the flavor symmetric points $m_\pi=m_K=412$ MeV at $a=0.055,0.064, 0.077, 0.094, 0.12$ fm and define the trajectories to the physical point by fixing the trace of the quark mass matrix. Currently our pion mass range extends down to $m_\pi\sim200$ MeV. We outline our tuning goals and strategy as well as our future planned ensembles. First scaling studies are performed on $f_\pi$ and $m_\pi$. Additionally results of a preliminary continuum extrapolation of $m_N$ at the flavor symmetric point are presented. Going further a first determination of the light and strange hadron spectrum chiral dependence is shown, which serves to check the quality of the action for precision measurements. We also investigate other quantities such as flowed gauge observables to study how the continuum limit is approached. Taken together we observe the SWF enable us to perform stable lattice simulations across a large range of parameters in mass, volume and lattice spacing. Pooling resources our new initiative has made our reported progress possible and through it we will share generated gauge ensembles under an open science philosophy., Comment: 18 pages, 8 figures, to be published in Proceedings of Science, contribution to "The 38th International Symposium on Lattice Field Theory (LATTICE 2021)"

The quark chromoelectric dipole (qCEDM) operator is a CP-violating operator describing, at hadronic energies, beyond-the-standard-model contributions to the electric dipole moment of particles with nonzero spin. In this paper we define renormalized dipole operators in a regularization-independent scheme using the gradient flow, and we perform the matching at one loop in perturbation theory to renormalized operators of the same and lower dimension in the more familiar MS scheme. We also determine the matching coefficients for the quark chromomagnetic dipole operator (qCMDM), which contributes, for example, to matrix elements relevant to CP-violating and CP-conserving kaon decays. The calculation provides a basis for future lattice QCD computations of hadronic matrix elements of the qCEDM and qCMDM operators., Comment: 24 pages, 4 figures Published version: added references, added discussion of nonperturbative effects in section 5.2. all results unchanged

It is not possible, using standard lattice techniques in Euclidean space, to calculate the complete fermionic spectrum of a quantum field theory. Algorithms running on quantum computers have the potential to access the theory with real-time evolution, enabling a direct computation. As a testing ground we consider the 1 + 1-dimensional Schwinger model with the presence of a {\theta} term using a staggered fermions discretization. We study the convergence properties of two different algorithms - adiabatic evolution and the Quantum Approximate Optimization Algorithm - with an emphasis on their cost in terms of CNOT gates. This is crucial to understand the feasibility of these algorithms, because calculations on near-term quantum devices depend on their rapid convergence. We also propose a blocked algorithm that has the first indications of a better scaling behavior with the dimensionality of the problem., Comment: 9 pages, 2 figures, 38th International Symposium on Lattice Field Theory, LATTICE2021 26th-30th July, 2021 Zoom/Gather@Massachusetts Institute of Technology

The CP-violating quark chromoelectric dipole moment (qCEDM) operator, contributing to the electric dipole moment (EDM), mixes under renormalization and particularly on the lattice with the pseudoscalar density. The mixing coefficient is power-divergent with the inverse lattice spacing squared, $1/a^2$, regardless of the lattice action used. We use the gradient flow to define a multiplicatively renormalized qCEDM operator and study its behavior at small flow time. We determine nonperturbatively the linearly divergent coefficient with the flow time, $1/t$, up to subleading logarithmic corrections, and compare it with the 1-loop perturbative expansion in the bare and renormalized strong coupling. We also discuss the O($a$) improvement of the qCEDM defined at positive flow time., Comment: 43 pages, 18 figures Added few sentences for clarification in the main text. Changed notation for the power divergent contribution in Eq.(51), Eq.(C23) and Tab.3 to avoid confusion. Version accepted for publication on PRD

Measurements of a permanent neutron electric dipole moment (EDM) potentially probe Beyond-the-Standard Model (BSM) sources of CP-violation. At low energy the CP-violating BSM interactions are parametrized by flavor-conserving CP-violating operators of dimension higher than four. QCD calculations of the nucleon matrix elements of these operators are required to fully reconstruct the sources and magnitudes of the different CP-violating contributions to the nucleon EDM. Herein we study the quark-chromo electric dipole moment (qCEDM) operator and the three-gluon Weinberg operator. The non-perturbative determination, using lattice QCD, of the nucleon matrix elements of these CP-violating operators is hampered by their short-distance behavior. Under renormalization these operators mix with lower dimensional operators, which induces power divergences in the lattice spacing, as the continuum limit is approached. We study the short-distance behavior of the qCEDM and the Weinberg operators using the gradient flow. We perform a short flow time expansion and determine, in perturbation theory, the expansion coefficients of the linearly-divergent terms stemming from the mixing with the pseudoscalar density and the topological charge, confirming the expectations of the operator product expansion. We introduce a new method to perform calculations at non-zero flow-time for arbitrary values of the external momenta. This method allows us to work in four dimensions for most of the calculations described in this paper, avoiding the complications associated with defining $\gamma_5$ in generic d dimensions. We show that leading contributions in the external momenta can be reproduced by defining $\gamma_5$ using the 't Hooft-Veltman-Breitenlohner-Maison scheme., Comment: 23 pages, 3 figures

We calculate the electric dipole moment of the nucleon induced by the QCD theta term. We use the gradient flow to define the topological charge and use $N_f = 2+1$ flavors of dynamical quarks corresponding to pion masses of $700$, $570$, and $410$ MeV, and perform an extrapolation to the physical point based on chiral perturbation theory. We perform calculations at $3$ different lattice spacings in the range of $0.07~{\rm fm} < a < 0.11$ fm at a single value of the pion mass, to enable control on discretization effects. We also investigate finite size effects using $2$ different volumes. A novel technique is applied to improve the signal-to-noise ratio in the form factor calculations. The very mild discretization effects observed suggest a continuum-like behavior of the nucleon EDM towards the chiral limit. Under this assumption our results read $d_{n}=-0.00152(71)\ \bar\theta\ e~\text{fm}$ and $d_{p}=0.0011(10)\ \bar\theta\ e~\text{fm}$. Assuming the theta term is the only source of CP violation, the experimental bound on the neutron electric dipole moment limits $\left|\bar\theta\right| < 1.98\times 10^{-10}$ ($90\%$ CL). A first attempt at calculating the nucleon Schiff moment in the continuum resulted in $S_{p} = 0.50(59)\times 10^{-4}\ \bar\theta\ e~\text{fm}^3$ and $S_{n} = -0.10(43)\times 10^{-4}\ \bar\theta\ e~\text{fm}^3$., Comment: Updated version submitted to journal. Added extraction of CP-odd pion-nucleon coupling and details about chiral extrapolation. Fixed bug in analysis code. Results mostly unchanged. Minor changes in the wording

We present a study on ab-initio calculations of three-neutron correlators from Lattice QCD. We extend the method of baryon blocks to systems of three spacially displaced baryons. This allows the measurement of three-neutron $p$-wave correlators with total spin $S=1/2$ and $3/2$. In addition, we use automatic code generation that has high flexibility and allows for easy inclusion of additional channels in the future while optimizing the evaluation of contractions. Our measurements were performed on a newly generated $96\times 48^3$ Clover-Wilson gauge field ensemble with $m_{\pi} \approx 370\,\text{MeV}$. We present preliminary results of our calculations of one pion and nucleon as well as two nucleon ($2N$) and three neutron ($3n$) correlators., Comment: Proceedings from Lattice2018, 7 pages, 4 figures

In this proceedings, we lay the foundation for computing the contribution of quark chromo-electric dipole moment (qCEDM) operator to the nucleon electric dipole moment. By applying the gradient flow technique, we can parameterize the renormalization and operator mixing issues associated with the qCEDM operator on the lattice. As the nucleon mixing angle $\alpha_N$ is a key component for determining the neutron and proton electric dipole moments induced by the qCEDM operator, we present the formalism and preliminary results for $\alpha_N$ with respect to the gradient flow time $t_f$. The results are computed on $N_f=2+1$ Wilson-clover lattices provided by PACS-CS. The 3 ensembles have lattice spacing values of $a=\lbrace 0.1095,\,0.0936,\,0.0684\rbrace $fm, whilst keeping a similar $m_{\pi}\approx701$MeV, and a fixed box size of $L\approx1.9$fm., Comment: 7 pages, 9 figures, presented at the 36th International Symposium on Lattice Field Theory (Lattice 2018)

We use the Yang-Mills gradient flow to study the mixing of CP-violating pure gauge operators in continuum QCD with special attention to Weinberg's d=6 purely gluonic operator. The gradient flow allows for a relatively clear derivation of the Wilson coefficients of the CP-violating effective Hamiltonian. This calculation is the first step towards a high-energy matching of matrix elements involving the CP-violating operators between the perturbative and lattice regimes., Comment: 7 pages, 1 figure, Proceedings from The 36th Annual International Symposium on Lattice Field Theory

Using the gradient flow, we calculated the nucleon mixing angle $\alpha_{N}$ and the nucleon electric dipole moment (EDM) induced by the QCD $\theta$-term. To do so, we computed the topological charge, and the nucleon two-point and three-point correlation functions. The purpose of these proceedings is to explore how the topological charge density interacts with the nucleon interpolation operators. By understanding this relation, we can try to suppress noise contributions to the $\alpha_{N}$ and EDM signals by selecting specific regions where the signal dominates. Using gauge fields provided by PACS-CS at $N_{f}=2+1$, a first collection of ensembles were selected at a fixed lattice spacing $a=0.0907$ fm ($\beta=1.90$), fixed dimensions $32^{3}\times 64$ and varying $m_{\pi}\approx\lbrace 411,\,570,\,701\rbrace$ MeV. A second collection was selected at fixed $m_{\pi}\approx701$ MeV, fixed box size of $L\approx1.9$ fm and varying $a=\lbrace 0.1215,\,0.0980,\,0.0685\rbrace $ fm., Comment: 7 pages, 12 figures, presented at the 36th International Symposium on Lattice Field Theory (Lattice 2018)

We utilize the gradient flow to define and calculate electric dipole moments induced by the strong QCD $\theta$-term and the dimension-6 Weinberg operator. The gradient flow is a promising tool to simplify the renormalization pattern of local operators. The results of the nucleon electric dipole moments are calculated on PACS-CS gauge fields (available from the ILDG) using $N_{f}=2+1$, of discrete size $32^{3}\times 64$ and spacing $a \simeq 0.09$ fm. These gauge fields use a renormalization-group improved gauge action and a non-perturbatively $O(a)$ improved clover quark action at $\beta = 1.90$, with $c_{SW} = 1.715$. The calculation is performed at pion masses of $m_{\pi} \simeq 411,701$ MeV., Comment: 8 pages, 13 figures, presented at the 35th International Symposium on Lattice Field Theory (Lattice 2017)

Inclusion of QED in lattice QCD calculations can lead to power-law volume artifacts as a consequence of the long-range nature of the interaction. Such artifacts must be removed by extrapolation in order to attain reliable infinite volume estimates of observables and quantities derived from them. As an alternative to this methodology, we consider the use of a photon mass term as an infrared regulator for QCD+QED, and explore the viability of its use in determining hadron mass shifts and splittings., Comment: 7 pages, 4 figures; presented at The 33rd International Symposium on Lattice Field Theory, Kobe International Conference Center, Kobe, Japan (July 14 -18, 2015)

Standard methods for including electromagnetic interactions in lattice quantum chromodynamics calculations result in power-law finite-volume corrections to physical quantities. Removing these by extrapolation requires costly computations at multiple volumes. We introduce a photon mass to alternatively regulate the infrared, and rely on effective field theory to remove its unphysical effects. Electromagnetic modifications to the hadron spectrum are reliably estimated with a precision and cost comparable to conventional approaches that utilize multiple larger volumes. A significant overall cost advantage emerges when accounting for ensemble generation. The proposed method may benefit lattice calculations involving multiple charged hadrons, as well as quantum many-body computations with long-range Coulomb interactions., Comment: 6 pages, 4 figures, 2 tables; significant revisions to abstract and main text; revised presentation of results for clarity (results unchanged); acknowledgements updated; matches published version

We propose a new method to calculate electric dipole moments induced by the strong QCD $\theta$-term. The method is based on the gradient flow for gauge fields and is free from renormalization ambiguities. We test our method by computing the nucleon electric dipole moments in pure Yang-Mills theory at several lattice spacings, enabling a first-of-its-kind continuum extrapolation. The method is rather general and can be applied for any quantity computed in a $\theta$ vacuum. This first application of the gradient flow has been successful and demonstrates proof-of-principle, thereby providing a novel method to obtain precise results for nucleon and light nuclear electric dipole moments., Comment: 32 pages, 14 figures, 2 tables. v2: added 1 plot, 1 table and 1 reference. Typos corrected. Published version

Non-singlet chiral Ward identities for fermionic operators at positive flow-time are derived using standard techniques based on local chiral variations of the action and of local operators. The gradient flow formalism is applied to twisted mass fermions and it is shown that automatic O(a) improvement for Wilson twisted mass fermions at maximal twist is a property valid also at positive flow-time. A definition of the chiral condensate that is multiplicatively renormalizable and automatically O(a) improved is then derived., Comment: 25 pages; v2: few typos corrected, one reference added, it matches the published version

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., Comment: 7 pages, 3 figures; presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We apply the spectral projector method, recently introduced by Giusti and L\"uscher, 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., Comment: 7 pages, 3 figures; presented at the 31st International Symposium on Lattice Field Theory (Lattice 2013), 29 July - 3 August 2013, Mainz, Germany

We investigate the continuum limit of the low lying static-light meson spectrum using Wilson twisted mass lattice QCD with N_f = 2 dynamical quark flavours. We consider three values of the lattice spacing a ~ 0.051 fm, 0.064 fm, 0.080 fm and various values of the pion mass in the range 280 MeV < m_PS < 640 MeV. We present results in the continuum limit for light cloud angular momentum j = 1/2, 3/2, 5/2 and for parity P = +, -. We extrapolate our results to physical quark masses, make predictions regarding the spectrum of B and B_s mesons and compare with available experimental results., Comment: 18 pages, 3 figures

The use of chirally rotated boundary conditions provides a formulation of the Schroedinger functional that is compatible with automatic O(a)-improvement of Wilson fermions in the bulk. The elimination of bulk O(a) terms requires the non-perturbative tuning of the critical mass and one additional boundary counterterm. We present the results of such a tuning in the quenched approximation at three values of the renormalised gauge coupling and for a range of lattice spacings., Comment: 7 pages, 4 figures, Proceeding-Lattice 2009

We make several observations concerning the low quark mass region with Wilson fermions and how this is connected with the epsilon regime in the continuum. A transition from tiny cutoff effects to rather large discretization errors would take place in general with Wilson fermions if we lower the quark mass at finite lattice spacing. We argue that these two regions exhibit rather different behaviours concerning the coupling between cutoff effects and zero-modes. We interpolate between these two regimes adding to the continuum epsilon regime formulae, in the spirit of the Symanzik expansion, the relevant operators parametrising the leading cutoff effects. We compute the partition function, the chiral condensate, scalar and pseudo-scalar correlation functions. The final formulae can be used to fit lattice data to extract physical low energy constants, and to estimate systematic uncertainties coming from discretization errors. Moreover they suggest ways on how to remove these cutoff effects, the core of which are captured by the continuum zero modes., Comment: 17 pages, 1 figure. v2: few typos corrected. Published on PLB

We compute the static-light meson spectrum with N_f = 2 flavours of sea quarks using Wilson twisted mass lattice QCD. We consider five different values for the light quark mass corresponding to 300 MeV < m_PS < 600 MeV and we present results for angular momentum j = 1/2, j = 3/2 and j = 5/2 and for parity P = + and P = -. We extrapolate our results to physical quark masses and make predictions regarding the spectrum of B and B_s mesons., Comment: 20 pages, 2 figures

We compute the static-light meson spectrum using two-flavor Wilson twisted mass lattice QCD. We have considered five different values for the light quark mass corresponding to 300 MeV < m_PS < 600 MeV. We have extrapolated our results, to make predictions regarding the spectrum of B and B_s mesons., Comment: 7 pages, 2 figures, talk given at the XXVI International Symposium on Lattice Field Theory, July 14 - 19 2008, Williamsburg, Virginia, USA

I report on first results of the ETMC obtained simulating lattice QCD with two degenerate flavours of Wilson twisted mass fermions, for which physical observables are automatically O(a) improved. Recent improvements of the HMC algorithm allow simulations of pseudoscalar masses below 300 MeV on volumes L^3 x 2L with L>2 fm and at values of the lattice spacing a<0.1 fm., Comment: 4 pages, 2 figures, talk given at ICHEP 2006

I review recent theoretical developments and numerical results of twisted mass QCD. I argue that, combined with an efficient algorithm, twisted mass QCD can be an attractive QCD lattice action, to perform large scale simulations at small pion masses, where a matching with chiral perturbation theory can be performed. Open issues like flavour breaking effects are also addressed., Comment: Plenary talk given at the 23rd International Symposium on Lattice Field Theory, Dublin 25-30 July 2005, 22 pages

We introduce new discretizations of the action for static quarks. They achieve an exponential improvement (compared to the Eichten-Hill regularization) on the signal to noise ratio in static-light correlation functions. This is explicitly checked in a quenched simulation and it is understood quantitatively in terms of the self energy of a static quark and the lattice heavy quark potential at zero distance. We perform a set of scaling tests in the Schroedinger functional and find scaling violations in the O(a) improved theory to be rather small -- for one observable significantly smaller than with the Eichten-Hill regularization. In addition we compute the improvement coefficients of the static light axial current up to O(g_0^4) corrections and the corresponding renormalization constants non-perturbatively. The regularization dependent part of the renormalization of the b-quark mass in static approximation is also determined., Comment: 30 pages, 9 figures

We present a first scaling test of twisted mass QCD with pure Wilson quarks for a twisting angle of pi/2. We have computed the vector meson mass and the pseudoscalar decay constant for different values of beta at fixed value of r_0 m_PS. The results obtained in the quenched approximation are compared with data for pure Wilson and non-perturbatively O(a) improved Wilson computations. We show that our results from Wilson twisted mass QCD show clearly reduced lattice spacing errors, consistent with O(a) improvement and without the need of any improvement terms added. These results thus provide numerical evidence of the prediction in ref. [1]., Comment: 15 pages, 4 figures; v2: two typos corrected, accepted for publication in Phys. Lett. B

We define, within the Schr\"odinger functional (SF) scheme, the matrix elements of the twist-2 operators corresponding to the first two moments of non-singlet parton density, and the first moment of singlet parton densities. We perform a lattice one-loop calculation that fixes the relation between the SF scheme and other common schemes and shows the main source of lattice artefacts. Few remarks on the improvement case are added., Comment: Presented at LATTICE99, 3 pages