Source author record

K. Jansen

K. Jansen 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

81works
15topics
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

81 published item(s)

preprint2024arXiv

Pion Transition Form Factor from Twisted-Mass Lattice QCD and the Hadronic Light-by-Light $π^0$-pole Contribution to the Muon $g-2$

The neutral pion generates the leading pole contribution to the hadronic light-by-light tensor, which is given in terms of the nonperturbative transition form factor $\mathcal{F}_{π^0γγ}(q_1^2,q_2^2)$. Here we present an ab-initio lattice calculation of this quantity in the continuum and at the physical point using twisted-mass lattice QCD. We report our results for the transition form factor parameterized using a model-independent conformal expansion valid for arbitrary space-like kinematics and compare it with experimental measurements of the single-virtual form factor, the two-photon decay width, and the slope parameter. We then use the transition form factors to compute the pion-pole contribution to the hadronic light-by-light scattering in the muon $g-2$, finding $a_μ^{π^0\text{-pole}} = 56.7(3.2) \times 10^{-11}$.

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.

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

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$.

preprint2019arXiv

Simulating Lattice Gauge Theories within Quantum Technologies

Lattice gauge theories, which originated from particle physics in the context of Quantum Chromodynamics (QCD), provide an important intellectual stimulus to further develop quantum information technologies. While one long-term goal is the reliable quantum simulation of currently intractable aspects of QCD itself, lattice gauge theories also play an important role in condensed matter physics and in quantum information science. In this way, lattice gauge theories provide both motivation and a framework for interdisciplinary research towards the development of special purpose digital and analog quantum simulators, and ultimately of scalable universal quantum computers. In this manuscript, recent results and new tools from a quantum science approach to study lattice gauge theories are reviewed. Two new complementary approaches are discussed: first, tensor network methods are presented - a classical simulation approach - applied to the study of lattice gauge theories together with some results on Abelian and non-Abelian lattice gauge theories. Then, recent proposals for the implementation of lattice gauge theory quantum simulators in different quantum hardware are reported, e.g., trapped ions, Rydberg atoms, and superconducting circuits. Finally, the first proof-of-principle trapped ions experimental quantum simulations of the Schwinger model are reviewed.

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.

preprint2016arXiv

Overcoming the sign problem in 1-dimensional QCD by new integration rules with polynomial exactness

In this paper we describe a new integration method for the groups $U(N)$ and $SU(N)$, for which we verified numerically that it is polynomially exact for $N\le 3$. The method is applied to the example of 1-dimensional QCD with a chemical potential. We explore, in particular, regions of the parameter space in which the sign problem appears due the presence of the chemical potential. While Markov Chain Monte Carlo fails in this region, our new integration method still provides results for the chiral condensate on arbitrary precision, demonstrating clearly that it overcomes the sign problem. Furthermore, we demonstrate that our new method leads to orders of magnitude reduced errors also in other regions of parameter space.

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)$.

preprint2015arXiv

On the efficient numerical solution of lattice systems with low-order couplings

We apply the Quasi Monte Carlo (QMC) and recursive numerical integration methods to evaluate the Euclidean, discretized time path-integral for the quantum mechanical anharmonic oscillator and a topological quantum mechanical rotor model. For the anharmonic oscillator both methods outperform standard Markov Chain Monte Carlo methods and show a significantly improved error scaling. For the quantum mechanical rotor we could, however, not find a successful way employing QMC. On the other hand, the recursive numerical integration method works extremely well for this model and shows an at least exponentially fast error scaling.

preprint2015arXiv

Thermal evolution of the one-flavour Schwinger model using Matrix Product States

The Schwinger model, or 1+1 dimensional QED, offers an interesting object of study, both at zero and non-zero temperature, because of its similarities to QCD. In this proceeding, we present the a full calculation of the temperature dependent chiral condensate of this model in the continuum limit using Matrix Product States (MPS). MPS methods, in general tensor networks, constitute a very promising technique for the non-perturbative study of Hamiltonian quantum systems. In the last few years, they have shown their suitability as ansatzes for ground states and low-lying excita- tions of lattice gauge theories. We show the feasibility of the approach also for finite temperature, both in the massless and in the massive case.

preprint2015arXiv

Thermal evolution of the Schwinger model with Matrix Product Operators

We demonstrate the suitability of tensor network techniques for describing the thermal evolution of lattice gauge theories. As a benchmark case, we have studied the temperature dependence of the chiral condensate in the Schwinger model, using matrix product operators to approximate the thermal equilibrium states for finite system sizes with non-zero lattice spacings. We show how these techniques allow for reliable extrapolations in bond dimension, step width, system size and lattice spacing, and for a systematic estimation and control of all error sources involved in the calculation. The reached values of the lattice spacing are small enough to capture the most challenging region of high temperatures and the final results are consistent with the analytical prediction by Sachs and Wipf over a broad temperature range.

preprint2014arXiv

Baryon spectrum using Nf=2+1+1 ensembles of twisted mass fermions

We present results on the masses of the low-lying baryons using ten ensembles of gauge configurations with $N_f =2+1+1$ dynamical twisted mass fermions, at three values of the lattice spacing, spanning a pion mass range from about 210 MeV to about 430 MeV. The strange and charm quark masses are tuned to approximately their physical values. We examine isospin symmetry breaking effects on the baryon mass and the dependence on the lattice spacing. After taking the continuum limit we use chiral perturbation theory to extrapolate to the physical vlaue of the pion mass for all forty baryons. We provide predictions for the masses of doubly and triply charmed baryons that have not yet been measured experimentally.

preprint2014arXiv

Baryon spectrum with $N_f=2+1+1$ twisted mass fermions

The masses of the low lying baryons are evaluated using a total of ten ensembles of dynamical twisted mass fermion gauge configurations. The simulations are performed using two degenerate flavors of light quarks, and a strange and a charm quark fixed to approximately their physical values. The light sea quarks correspond to pseudo scalar masses in the range of about 210~MeV to 430~MeV. We use the Iwasaki improved gluonic action at three values of the coupling constant corresponding to lattice spacing $a=0.094$~fm, 0.082~fm and 0.065~fm determined from the nucleon mass. We check for both finite volume and cut-off effects on the baryon masses. We examine the issue of isospin symmetry breaking for the octet and decuplet baryons and its dependence on the lattice spacing. We show that in the continuum limit isospin breaking is consistent with zero, as expected. We performed a chiral extrapolation of the forty baryon masses using SU(2) $χ$PT. After taking the continuum limit and extrapolating to the physical pion mass our results are in good agreement with experiment. We provide predictions for the mass of the doubly charmed $Ξ_{cc}^*$, as well as of the doubly and triply charmed $Ω$s that have not yet been determined experimentally.

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

Progress in Simulations with Twisted Mass Fermions at the Physical Point

In this contribution, results from $N_f=2$ lattice QCD simulations at one lattice spacing using twisted mass fermions with a clover term at the physical pion mass are presented. The mass splitting between charged and neutral pions (including the disconnected contribution) is shown to be around $20(20)~\mathrm{MeV}$. Further, a first measurement using the clover twisted mass action of the average momentum fraction of the pion is given. Finally, an analysis of pseudoscalar meson masses and decay constants is presented involving linear interpolations in strange and charm quark masses. Matching to meson mass ratios allows the calculation of quark mass ratios: $μ_s/μ_l=27.63(13)$, $μ_c/μ_l=339.6(2.2)$ and $μ_c/μ_s=12.29(10)$. From this mass matching the quantities $f_K=153.9(7.5)~\mathrm{MeV}$, $f_D=219(11)~\mathrm{MeV}$, $f_{D_s}=255(12)~\mathrm{MeV}$ and $M_{D_s}=1894(93)~\mathrm{MeV}$ are determined without the application of finite volume or discretization artefact corrections and with errors dominated by a preliminary estimate of the lattice spacing.

preprint2014arXiv

Sigma-terms and axial charges for hyperons and charmed baryons

We present results for the $σ$-terms and axial charges for various hyperons and charmed baryons using $N_f=2+1+1$ twisted mass fermions. For the computation of the three-point function we use the fixed current method. For one of the $N_f=2+1+1$ ensembles with pion mass of 373 MeV we compare the results of the fixed current method with those obtained with a stochastic method for computing the all-to-all propagator involved in the evaluation of the three point functions.

preprint2013arXiv

A first look at maximally twisted mass lattice QCD calculations at the physical point

In this contribution, a first look at simulations using maximally twisted mass Wilson fermions at the physical point is presented. A lattice action including clover and twisted mass terms is presented and the Monte Carlo histories of one run with two mass-degenerate flavours at a single lattice spacing are shown. Measurements from the light and heavy-light pseudoscalar sectors are compared to previous $N_f = 2$ results and their phenomenological values. Finally, the strategy for extending simulations to $N_f = 2 + 1 + 1$ is outlined.

preprint2013arXiv

Experiences with OpenMP in tmLQCD

An overview is given of the lessons learned from the introduction of multi-threading using OpenMP in tmLQCD. In particular, programming style, performance measurements, cache misses, scaling, thread distribution for hybrid codes, race conditions, the overlapping of communication and computation and the measurement and reduction of certain overheads are discussed. Performance measurements and sampling profiles are given for different implementations of the hopping matrix computational kernel.

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

Quasi-Monte Carlo methods for lattice systems: a first look

We investigate the applicability of Quasi-Monte Carlo methods to Euclidean lattice systems for quantum mechanics in order to improve the asymptotic error behavior of observables for such theories. In most cases the error of an observable calculated by averaging over random observations generated from an ordinary Markov chain Monte Carlo simulation behaves like N^{-1/2}, where N is the number of observations. By means of Quasi-Monte Carlo methods it is possible to improve this behavior for certain problems to N^{-1}, or even further if the problems are regular enough. We adapted and applied this approach to simple systems like the quantum harmonic and anharmonic oscillator and verified an improved error scaling.

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}.

preprint2012arXiv

A first look at quasi-Monte Carlo for lattice field theory problems

In this project we initiate an investigation of the applicability of Quasi-Monte Carlo methods to lattice field theories in order to improve the asymptotic error behavior of observables for such theories. In most cases the error of an observable calculated by averaging over random observations generated from an ordinary Monte Carlo simulation behaves like 1/sqrt(N), where N is the number of observations. By means of Quasi-Monte Carlo methods it is possible to improve this behavior for certain problems to up to 1/N. We adapted and applied this approach to simple systems like the quantum harmonic and anharmonic oscillator and verified an improved error scaling.

preprint2012arXiv

A quenched study of the Schroedinger functional with chirally rotated boundary conditions: applications

In a previous paper [1], we have discussed the non-perturbative tuning of the chirally rotated Schroedinger functional (XSF). This tuning is required to eliminate bulk O(a) cutoff effects in physical correlation functions. Using our tuning results obtained in [1] we perform scaling and universality tests analyzing the residual O(a) cutoff effects of several step-scaling functions and we compute renormalization factors at the matching scale. As an example of possible application of the XSF we compute the renormalized strange quark mass using large volume data obtained from Wilson twisted mass fermions at maximal twist.

preprint2012arXiv

A quenched study of the Schroedinger functional with chirally rotated boundary conditions: non-perturbative tuning

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 up to O(a) boundary contributions. The elimination of bulk O(a) effects requires the non-perturbative tuning of the critical mass and one additional boundary counterterm. We present the results of such a tuning in a quenched setup for several values of the renormalized gauge coupling, from perturbative to non-perturbative regimes, and for a range of lattice spacings. We also check that the correct boundary conditions and symmetries are restored in the continuum limit.

preprint2012arXiv

Evaluation of disconnected contributions using GPUs

We calculate on GPUs the disconnected diagrams associated with the nucleon form factors and moments of generalized parton distributions using Nf=2+1+1 twisted mass fermions. We employ the truncated solver method (TSM) for estimating the all-to-all propagators. Due to the fact that the TSM involves many low precision stochastic estimators, the usage of GPUs is essential to perform efficiently the contractions and the inversions.

preprint2012arXiv

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

We investigate scalar matrix elements of the nucleon using $N_f=2+1+1$ flavors of maximally twisted mass fermions at a fixed value of the lattice spacing of $a\approx 0.078{\rm fm}$. We compute disconnected contributions to the relevant three-point functions using an efficient noise reduction technique. Using these methods together with an only multiplicative renormalization applicable for twisted mass fermions, allows us to obtain accurate results in the light and strange sector.

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.

preprint2012arXiv

Topological susceptibility and chiral condensate with $N_f=2+1+1$ dynamical flavors of maximally twisted mass fermions

We study the 'spectral projector' method for the computation of the chiral condensate and the topological susceptibility, using $N_f=2+1+1$ dynamical flavors of maximally twisted mass Wilson fermions. In particular, we perform a study of the quark mass dependence of the chiral condensate $Σ$ and topological susceptibility $χ_{top}$ in the range $270 MeV < m_π < 500 MeV$ and compare our data with analytical predictions. In addition, we compute $χ_{top}$ in the quenched approximation where we match the lattice spacing to the $N_f=2+1+1$ dynamical simulations. Using the Kaon, $η$ and $η^{\prime}$ meson masses computed on the $N_f=2+1+1$ ensembles, we then perform a preliminary test of the Witten-Veneziano relation.

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

Hadronic Vacuum Polarization Contribution to g-2 from the Lattice

We give a short description of the present situation of lattice QCD simulations. We then focus on the computation of the anomalous magnetic moment of the muon using lattice techniques. We demonstrate that by employing improved observables for the muon anomalous magnetic moment, a significant reduction of the lattice error can be obtained. This provides a promising scenario that the accuracy of lattice calculations can match the experimental errors.

preprint2011arXiv

Mass loss out of close binaries. II

Liberal evolution of interacting binaries has been proposed previously by several authors in order to meet various observed binary characteristics better than conservative evolution does. Since Algols are eclipsing binaries the distribution of their orbital periods is precisely known. The distribution of their mass ratios contains however more uncertainties. We try to reproduce these two distributions theoretically using a liberal scenario in which the gainer star can lose mass into interstellar space as a consequence of its rapid rotation and the energy of a hot spot. In a recent paper (Van Rensbergen et al. 2010, A&A) we calculated the liberal evolution of binaries with a B-type primary at birth where mass transfer starts during core hydrogen burning of the donor. In this paper we include the cases where mass transfer starts during hydrogen shell burning and it is our aim to reproduce the observed distributions of the system parameters of Algol-type semi-detached systems. Our calculations reveal the amount of time that an Algol binary lives with a well defined value of mass ratio and orbital period. We use these data to simulate the distribution of mass ratios and orbital periods of Algols. Binaries with a late B-type initial primary hardly lose any mass whereas those with an early B primary evolve in a non-conservative way. Conservative binary evolution predicts only ~ 12 % of Algols with a mass ratio q above 0.4. This value is raised up to ~ 17 % using our scenario of liberal evolution, which is still far below the ~ 45 % that is observed. Observed orbital periods of Algol binaries larger than one day are faithfully reproduced by our liberal scenario. Mass ratios are reproduced better than with conservative evolution, but the resemblance is still poor.

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

Renormalisation constants of quark bilinears in lattice QCD with four dynamical Wilson quarks

We present preliminary results of the non-perturbative computation of the RI-MOM renormalisation constants in a mass-independent scheme for the action with Iwasaki glue and four dynamical Wilson quarks employed by ETMC. Our project requires dedicated gauge ensembles with four degenerate sea quark flavours at three lattice spacings and at several values of the standard and twisted quark mass parameters. The RI-MOM renormalisation constants are obtained from appropriate O(a) improved estimators extrapolated to the chiral limit.

preprint2011arXiv

Renormalization constants for Wilson fermion lattice QCD with four dynamical flavours

We report on an ongoing non-perturbative computation of RI-MOM scheme renormalization constants for the lattice action with four dynamical flavours currently in use by ETMC. For this goal dedicated simulations with four degenerate sea quark flavours are performed at several values of the standard and twisted quark mass parameters. We discuss a method for removing possible O(a) artifacts at all momenta and extrapolating renormalization constant estimators to the chiral limit. We give preliminary results at one lattice spacing.

preprint2010arXiv

A proposal for B-physics on current lattices

A method to extract B-physics parameters (b-quark mass and fB, fBs decay constants) from currently available lattice data is presented and tested. The approach is based on the idea of constructing appropriate ratios of heavy-light meson masses and decay constants, respectively, possessing a precisely known static limit, and evaluating them at various pairs of heavy quark masses around the charm. Via a smooth interpolation in the heavy quark mass from the easily accessible charm region to the asymptotic point, B-physics parameters are computed with a few percent (statistical + systematic) error using recently produced Nf=2 maximally twisted Wilson fermions data.

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

Effects of a potential fourth fermion generation on the Higgs boson mass bounds

We study the effect of a potential fourth fermion generation on the upper and lower Higgs boson mass bounds. This investigation is based on the numerical evaluation of a chirally invariant lattice Higgs-Yukawa model emulating the same Higgs-fermion coupling structure as in the Higgs sector of the electroweak Standard Model. In particular, the considered model obeys a Ginsparg-Wilson version of the underlying $SU(2)_L\times U(1)_Y$ symmetry, being a global symmetry here due to the neglection of gauge fields in this model. We present our results on the modification of the upper and lower Higgs boson mass bounds induced by the presence of a hypothetical very heavy fourth quark doublet. Finally, we compare these findings to the standard scenario of three fermion generations.

preprint2010arXiv

Higgs boson mass bounds in the presence of a very heavy fourth quark generation

We study the effect of a potential fourth quark generation on the upper and lower Higgs boson mass bounds. This investigation is based on the numerical evaluation of a chirally invariant lattice Higgs-Yukawa model emulating the same Higgs-fermion coupling structure as in the Higgs sector of the electroweak Standard Model. In particular, the considered model obeys a Ginsparg-Wilson version of the underlying ${SU}(2)_L\times {U}(1)_Y$ symmetry, being a global symmetry here due to the neglection of gauge fields in this model. We present our results on the modification of the upper and lower Higgs boson mass bounds induced by the presence of a hypothetical very heavy fourth quark doublet. Finally, we compare these findings to the standard scenario of three fermion generations.

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

Upper Higgs boson mass bounds from a chirally invariant lattice Higgs-Yukawa model

We establish the cutoff-dependent upper Higgs boson mass bound by means of direct lattice computations in the framework of a chirally invariant lattice Higgs-Yukawa model emulating the same chiral Yukawa coupling structure as in the Higgs-fermion sector of the Standard Model. As expected from the triviality picture of the Higgs sector, we observe the upper mass bound to decrease with rising cutoff parameter $Λ$. Moreover, the strength of the fermionic contribution to the upper mass bound is explored by comparing to the corresponding analysis in the pure $Φ^4$-theory.

preprint2009arXiv

f_B and f_Bs with maximally twisted Wilson fermions

We present a lattice QCD calculation of the heavy-light decay constants f_B and f_Bs performed with Nf=2 maximally twisted Wilson fermions, at four values of the lattice spacing. The decay constants have been also computed in the static limit and the results are used to interpolate the observables between the charm and the infinite-mass sectors, thus obtaining the value of the decay constants at the physical b quark mass. Our preliminary results are f_B=191(14) MeV, f_Bs=243(14) MeV, f_Bs/f_B=1.27(5). They are in good agreement with those obtained with a novel approach, recently proposed by our Collaboration (ETMC), based on the use of suitable ratios having an exactly known static limit.

preprint2009arXiv

Light Meson Physics from Maximally Twisted Mass Lattice QCD

We present a comprehensive investigation of light meson physics using maximally twisted mass fermions for two mass-degenerate quark flavours. By employing four values of the lattice spacing, spatial lattice extents ranging from 2.0 fm to 2.5 fm and pseudo scalar masses in the range 280 MeV to 650 MeV we control the major systematic effects of our calculation. This enables us to confront our data with chiral perturbation theory and extract low energy constants of the effective chiral Lagrangian and derived quantities, such as the light quark mass, with high precision.

preprint2009arXiv

Mass loss out of close binaries

In a liberal evolutionary scenario, mass can escape from a binary during eras of fast mass transfer. We calculate the mass lost by binaries with a B-type primary at birth where mass transfer starts during hydrogen core burning of the donor. We simulate the distribution of mass-ratios and orbital periods for those interacting binaries. The amount of time the binary shows Algol characteristics within different values of mass-ratio and orbital period has been fixed from conservative and liberal evolutionary calculations. We use these data to simulate the distribution of mass-ratios and orbital periods of Algols with the conservative as well as the liberal model. We compare mass-ratios and orbital periods of Algols obtained by conservative evolution with those obtained by our liberal model. Since binaries with a late B-type primary evolve almost conservatively, the overall distribution of mass-ratios will only yield a few Algols more with high mass-ratios than conservative calculations do. Whereas the simulated distribution of orbital periods of Algols fits the observations well, the simulated distribution of mass-ratios produces always too few systems with large values.

preprint2009arXiv

Monte Carlo approach to turbulence

The behavior of the one-dimensional random-force-driven Burgers equation is investigated in the path integral formalism on a discrete space-time lattice. We show that by means of Monte Carlo methods one may evaluate observables, such as structure functions, as ensemble averages over different field realizations. The regularization of shock solutions to the zero-viscosity limit (Hopf-equation) eventually leads to constraints on lattice parameters required for the stability of the simulations. Insight into the formation of localized structures (shocks) and their dynamics is obtained.

preprint2009arXiv

Nucleon form factors with Nf=2 dynamical twisted mass fermions

We present results on the electromagnetic and axial nucleon form factors using two degenerate flavors of twisted mass fermions on lattices of spatial size 2.1 fm and 2.7 fm and a lattice spacing of about 0.09 fm. We consider pion masses in the range of 260-470 MeV. We chirally extrapolate results on the nucleon axial ch arge, the isovector Dirac and Pauli root mean squared radii and magnetic moment to the physical point and co mpare to experiment.

preprint2009arXiv

Partially quenched study of strange baryon with Nf = 2 twisted mass fermions

We present results on the mass of the baryon octet and decuplet using two flavors of light dynamical twisted mass fermions. The strange quark mass is fixed to its physical value from the kaon sector in a partially quenched set up. Calculations are performed for light quark masses corresponding to a pion mass in the range 270-500 MeV and lattice sizes of 2.1 fm and 2.7 fm. We check for cut-off effects and isospin breaking by evaluating the baryon masses at two different lattice spacings. We carry out a chiral extrapolation for the octet baryons and discuss results for the Omega.

preprint2009arXiv

Phase structure of thermal lattice QCD with N_f=2 twisted mass Wilson fermions

We present numerical results for the phase diagram of lattice QCD at finite temperature in the formulation with twisted mass Wilson fermions and a tree-level Symanzik-improved gauge action. Our simulations are performed on lattices with temporal extent N_τ=8, and lattice coupling βranging from strong coupling to the scaling domain. Covering a wide range in the space spanned by the lattice coupling βand the hopping and twisted mass parameters κand μ, respectively, we obtain a comprehensive picture of the rich phase structure of the lattice theory. In particular, we verify the existence of an Aoki phase in the strong coupling region and the realisation of the Sharpe-Singleton scenario at intermediate couplings. In the weak coupling region we identify the phase boundary for the physical finite temperature phase transition/crossover. Its shape in the three-dimensional parameter space is consistent with Creutz's conjecture of a cone-shaped thermal transition surface.

preprint2009arXiv

Pseudoscalar decay constants of kaon and D-mesons from Nf=2 twisted mass Lattice QCD

We present the results of a lattice QCD calculation of the pseudoscalar meson decay constants fpi, fK, fD and fDs, performed with Nf=2 dynamical fermions. The simulation is carried out with the tree-level improved Symanzik gauge action and with the twisted mass fermionic action at maximal twist. We have considered for the final analysis three values of the lattice spacing, a~0.10 fm, 0.09 fm and 0.07 fm, with pion masses down to mpi~270 MeV. Our results for the light meson decay constants are fK=158.1(2.4) MeV and fK/fpi=1.210(18). From the latter ratio, by using the experimental determination of Gamma(K-->mu nu_mu (gamma))/ Gamma(pi--> mu nu_mu (gamma)) and the average value of |Vud| from nuclear beta decays, we obtain |Vus|=0.2222(34), in good agreement with the determination from semileptonic Kl3 decays and the unitarity constraint. For the D and Ds meson decay constants we obtain fD=197(9) MeV, fDs=244(8) MeV and fDs/fD=1.24(3). Our result for fD is in good agreement with the CLEO experimental measurement. For fDs our determination is smaller than the PDG 2008 experimental average but in agreement with a recent improved measurement by CLEO at the 1.4 sigma level.

preprint2009arXiv

The epsilon regime of chiral perturbation theory with Wilson-type fermions

In this proceeding contribution we report on the ongoing effort to simulate Wilson-type fermions in the so called epsilon regime of chiral perturbation theory. We present results for the chiral condensate and the pseudoscalar decay constant obtained with Wilson twisted mass fermions employing two lattice spacings, two different physical volumes and several quark masses. With this set of simulations we make a first attempt to estimate the systematic uncertainties.

preprint2009arXiv

tmLQCD: a program suite to simulate Wilson Twisted mass Lattice QCD

We discuss a program suite for simulating Quantum Chromodynamics on a 4-dimensional space-time lattice. The basic Hybrid Monte Carlo algorithm is introduced and a number of algorithmic improvements are explained. We then discuss the implementations of these concepts as well as our parallelisation strategy in the actual simulation code. Finally, we provide a user guide to compile and run the program.

preprint2008arXiv

Dynamical Twisted Mass Fermions with Light Quarks: Simulation and Analysis Details

In a recent paper [hep-lat/0701012] we presented precise lattice QCD results of our European Twisted Mass Collaboration (ETMC). They were obtained by employing two mass-degenerate flavours of twisted mass fermions at maximal twist. In the present paper we give details on our simulations and the computation of physical observables. In particular, we discuss the problem of tuning to maximal twist, the techniques we have used to compute correlators and error estimates. In addition, we provide more information on the algorithm used, the autocorrelation times and scale determination, the evaluation of disconnected contributions and the description of our data by means of chiral perturbation theory formulae.

preprint2008arXiv

The eta' meson from lattice QCD

We study the flavour singlet pseudoscalar mesons from first principles using lattice QCD. With N_f=2 flavours of light quark, this is the so-called eta_2 meson and we discuss the phenomenological status of this. Using maximally twisted-mass lattice QCD, we extract the mass of the eta_2 meson at two values of the lattice spacing for lighter quarks than previously discussed in the literature. We are able to estimate the mass value in the limit of light quarks with their physical masses.

preprint2008arXiv

The finite-temperature phase structure of lattice QCD with twisted-mass Wilson fermions

We report progress in our exploration of the finite-temperature phase structure of two-flavour lattice QCD with twisted-mass Wilson fermions and a tree-level Symanzik-improved gauge action for a temporal lattice size N_τ=8. Extending our investigations to a wider region of parameter space we gain a global view of the rich phase structure. We identify the finite temperature transition/crossover for a non-vanishing twisted-mass parameter in the neighbourhood of the zero-temperature critical line at sufficiently high β. Our findings are consistent with Creutz's conjecture of a conical shape of the finite temperature transition surface. Comparing with NLO lattice chiral perturbation theory we achieve an improved understanding of this shape.

preprint2007arXiv

Dynamical Twisted Mass Fermions with Light Quarks

We present results of dynamical simulations with 2 flavours of degenerate Wilson twisted mass quarks at maximal twist in the range of pseudo scalar masses from 300 to 550 MeV. The simulations are performed at one value of the lattice spacing a \lesssim 0.1 fm. In order to have O(a) improvement and aiming at small residual cutoff effects, the theory is tuned to maximal twist by requiring the vanishing of the untwisted quark mass. Precise results for the pseudo scalar decay constant and the pseudo scalar mass are confronted with chiral perturbation theory predictions and the low energy constants F, \bar{l}_3 and \bar{l}_4 are evaluated with small statistical errors.

preprint2007arXiv

Light quark masses and pseudoscalar decay constants from Nf=2 Lattice QCD with twisted mass fermions

We present the results of a lattice QCD calculation of the average up-down and strange quark masses and of the light meson pseudoscalar decay constants with Nf=2 dynamical fermions. The simulation is carried out at a single value of the lattice spacing with the twisted mass fermionic action at maximal twist, which guarantees automatic O(a)-improvement of the physical quantities. Quark masses are renormalized by implementing the non-perturbative RI-MOM renormalization procedure. Our results for the light quark masses are m_ud^{msbar}(2 GeV)= 3.85 +- 0.12 +- 0.40 MeV, m_s^{msbar}(2 GeV) = 105 +- 3 +- 9 MeV and m_s/m_ud = 27.3 +- 0.3 +- 1.2. We also obtain fK = 161.7 +- 1.2 +- 3.1 MeV and the ratio fK/fpi=1.227 +- 0.009 +- 0.024. From this ratio, by using the experimental determination of Gamma(K-> mu nu (gamma))/Gamma(pi -> mu nu (gamma)) and the average value of |Vud| from nuclear beta decays, we obtain |Vus|=0.2192(5)(45), in agreement with the determination from Kl3 decays and the unitarity constraint.

preprint2007arXiv

The phase structure of a chirally invariant lattice Higgs-Yukawa model for small and for large values of the Yukawa coupling constant

We consider a chirally invariant lattice Higgs-Yukawa model based on the Neuberger overlap operator. As a first step towards the eventual determination of Higgs mass bounds we study the phase diagram of the model analytically in the large Nf-limit. We present an expression for the effective potential at tree-level in the regime of small Yukawa and quartic coupling constants and determine the order of the phase transitions. In the case of strong Yukawa couplings the model effectively becomes an O(4)-symmetric non-linear sigma-model for all values of the quartic coupling constant. This leads to the existence of a symmetric phase also in the regime of large values of the Yukawa coupling constant. On finite and small lattices, however, strong finite volume effects prevent the expectation value of the Higgs field from vanishing thus obscuring the existence of the symmetric phase at strong Yukawa couplings.

preprint2006arXiv

Iterative methods for overlap and twisted mass fermions

We present a comparison of a number of iterative solvers of linear systems of equations for obtaining the fermion propagator in lattice QCD. In particular, we consider chirally invariant overlap and chirally improved Wilson (maximally) twisted mass fermions. The comparison of both formulations of lattice QCD is performed at four fixed values of the pion mass between 230MeV and 720MeV. For overlap fermions we address adaptive precision and low mode preconditioning while for twisted mass fermions we discuss even/odd preconditioning. Taking the best available algorithms in each case we find that calculations with the overlap operator are by a factor of 30-120 more expensive than with the twisted mass operator.

preprint2005arXiv

Light quarks with twisted mass fermions

We investigate Wilson twisted mass fermions in the quenched approximation using different definitions of the critical bare quark mass m_c to realize maximal twist and, correspondingly, automatic O(a) improvement for physical observables. A particular definition of m_c is given by extrapolating the value of m_c obtained from the PCAC relation at non-vanishing bare twisted quark mass mu to mu=0. Employing this improved definition of the critical mass the Wilson twisted mass formulation provides the possibility to perform reliable simulations down to very small quark masses with correspondingly small pion masses of m_pi \simeq 250 MeV, while keeping the cutoff effects of O(a^2) under control.

preprint2005arXiv

Nf=2 Lattice QCD and Chiral Perturbation Theory

By employing a twisted mass term, we compare recent results from lattice calculations of Nf=2 dynamical Wilson fermions with Wilson Chiral Perturbation Theory (WChPT). The final goal is to determine some combinations of Gasser-Leutwyler Low Energy Constants (LECs). A wide set of data with different lattice spacings (a ~ 0.2-0.12 fm), different gauge actions (Wilson plaquette, DBW2) and different quark masses (down to the lowest pion mass allowed by lattice artifacts and including negative quark masses) provide a strong check of the applicability of WChPT in this regime and the scaling behaviours in the continuum limit.

preprint2005arXiv

Non-perturbative Pion Matrix Element of a twist-2 operator from the Lattice

We give a continuum limit value of the lowest moment of a twist-2 operator in pion states from non-perturbative lattice calculations. We find that the non-perturbatively obtained renormalization group invariant matrix element is <x>_{RGI} = 0.179(11), which corresponds to <x>^{MSbar}(2 GeV) = 0.246(15). In obtaining the renormalization group invariant matrix element, we have controlled important systematic errors that appear in typical lattice simulations, such as non-perturbative renormalization, finite size effects and effects of a non-vanishing lattice spacing. The crucial limitation of our calculation is the use of the quenched approximation. Another question that remains not fully clarified is the chiral extrapolation of the numerical data.

preprint2005arXiv

Scaling test of fermion actions in the Schwinger model

We discuss the scaling behaviour of different fermion actions in dynamical simulations of the 2-dimensional massive Schwinger model. We have chosen Wilson, hypercube, twisted mass and overlap fermion actions. As physical observables, the pion mass and the scalar condensate are computed for the above mentioned actions at a number of coupling values and fermion masses. We also discuss possibilities to simulate overlap fermions dynamically avoiding problems with low-lying eigenvalues of the overlap kernel.

preprint2005arXiv

Testing Topology Conserving Gauge Actions for Lattice QCD

We explore gauge actions for lattice QCD, which are constructed such that the occurrence of small plaquette values is strongly suppressed. Such actions originate from the admissibility condition in order to conserve the topological charge. The suppression of small plaquette values is expected to be advantageous for numerical studies in the $ε$-regime and also for simulations with dynamical quarks. Performing simulations at a lattice spacing of about 0.1 fm, we present numerical results for the static potential, the physical scale $r_0$, the stability of the topological charge history, the condition number of the kernel of the overlap operator and the acceptance rate against the step size in the local HMC algorithm.

preprint2005arXiv

The phase structure of lattice QCD with two flavours of Wilson quarks and renormalization group improved gluons

The effect of changing the lattice action for the gluon field on the recently observed [1] first order phase transition near zero quark mass is investigated by replacing the Wilson plaquette action by the DBW2 action. The lattice action for quarks is unchanged: it is in both cases the original Wilson action. It turns out that Wilson fermions with the DBW2 gauge action have a phase structure where the minimal pion mass and the jump of the average plaquette are decreased, when compared to Wilson fermions with Wilson plaquette action at similar values of the lattice spacing. Taking the DBW2 gauge action is advantageous also from the point of view of the computational costs of numerical simulations.

preprint2004arXiv

Going chiral: overlap versus twisted mass fermions

We compare the behavior of overlap fermions, which are chirally invariant, and of Wilson twisted mass fermions at full twist in the approach to the chiral limit. Our quenched simulations reveal that with both formulations of lattice fermions pion masses of O(250 MeV) can be reached in practical applications. Our comparison is done at a fixed value of the lattice spacing a=0.123 fm. A number of quantities are measured such as hadron masses, pseudoscalar decay constants and quark masses obtained from Ward identities. We also determine the axial vector renormalization constants in the case of overlap fermions.

preprint2003arXiv

Finite-Size Scaling of Vector and Axial Current Correlators

Using quenched chiral perturbation theory, we compute the long-distance behaviour of two-point functions of flavour non-singlet axial and vector currents in a finite volume, for small quark masses, and at a fixed gauge-field topology. We also present the corresponding predictions for the unquenched theory at fixed topology. These results can in principle be used to measure the low-energy constants of the chiral Lagrangian, from lattice simulations in volumes much smaller than one pion Compton wavelength. We show that quenching has a dramatic effect on the vector correlator, which is argued to vanish to all orders, while the axial correlator appears to be a robust observable only moderately sensitive to quenching.

preprint2001arXiv

Speeding up the Hybrid-Monte-Carlo algorithm for dynamical fermions

We propose a modification of the Hybrid-Monte-Carlo algorithm that allows for a larger step-size of the integration scheme at constant acceptance rate. The key ingredient is the splitting of the pseudo-fermion action into two parts. We test our proposal at the example of the two-dimensional lattice Schwinger model and four-dimensional lattice QCD with two degenerate flavours of Wilson- fermions.