Source author record

P. E. L. Rakow

P. E. L. Rakow 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

92works
5topics
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

92 published item(s)

preprint2022arXiv

Exploring interpolating momentum schemes

We compute the renormalisation factors of the quark mass and wave function using IMOM (Interpolating MOMenta) schemes. The framework is the Rome-Southampton non-renormalisation method, but the momentum transfer in the quark bilinears is not restricted to zero or to the symmetric point. We study the scale dependence, infrared contamination and lattice artefacts for different values of this momentum transfer and for two different kinds of projectors. For the numerical simulations, we use data generated by the RBC-UKQCD collaborations, with $N_f = 2+1$ flavours of Domain-Wall fermions, and inverse lattice spacing of $1.79 $ and $2.38$ GeV.

preprint2022arXiv

Measurements of $SU(3)_f$ symmetry breaking in $B$ meson decay constants

We present updates from QCDSF/UKQCD/CSSM on the $SU(3)_f$ breaking in $B$ meson decay constants. The $b$-quarks are generated with an anisotropic clover-improved action, and are tuned to match properties of the physical $B$ and $B^*$ mesons. Configurations are generated with $\overline{m}=(1/3)(2m_l+m_s)$ kept constant to control symmetry breaking effects. Various sources of systematic uncertainty will be discussed, including those from continuum extrapolations and extrapolations to the physical point. We also present new efforts to calculate $f_B$ and $f_{B_s}$ using weighted averages across multiple time fitting regions. The use of an automated weighted averaging technique over multiple fitting ranges allows for timely tuning of the $b$-quark and reduces the impact of systematic errors from fitting range biases in calculations of $f_B$ and $f_{B_s}$

preprint2022arXiv

Nucleon Form Factors from the Feynman-Hellmann Method in Lattice QCD

Lattice QCD calculations of the nucleon electromagnetic form factors are of interest at both the high and low momentum transfer regions. For high momentum transfers especially there are open questions which require more intense study, such as the potential zero crossing in the proton's electric form factor. We will present recent progress from the QCDSF/UKQCD/CSSM collaboration on the calculation of these form factors using the Feynman-Hellmann method in lattice QCD. The Feynman-Hellmann method allows for greater control over excited states which we take advantage of by going to high values of the momentum transfer. In this proceeding we present results of the form factors up to $6 \textrm{GeV}^{2}$, using $N_{f}=2+1$ flavour fermions for three different pion masses in the range 310-470 $\textrm{MeV}$. The results are extrapolated to the physical pion mass through the use of a flavour breaking expansion.

preprint2022arXiv

The Compton Amplitude, lattice QCD and the Feynman-Hellmann approach

A major objective of lattice QCD is the computation of hadronic matrix elements. The standard method is to use three-point and four-point correlation functions. An alternative approach, requiring only the computation of two-point correlation functions is to use the Feynman-Hellmann theorem. In this talk we develop this method up to second order in perturbation theory, in a context appropriate for lattice QCD. This encompasses the Compton Amplitude (which forms the basis for deep inelastic scattering) and hadron scattering. Some numerical results are presented showing results indicating what this approach might achieve.

preprint2021arXiv

State mixing and masses of the $π^0$, $η$ and $η^\prime$ mesons from $n_f=1+1+1$ lattice QCD+QED

We present a lattice analysis of the light pseudoscalar mesons with consideration for the mixing between the flavour-neutral states $π^0$, $η$ and $η^\prime$. We extract the masses and flavour compositions of the pseudoscalar meson nonet in $n_f=1+1+1$ lattice QCD+QED around an SU(3)-flavour symmetric point, and observe flavour-symmetry features of the extracted data, along with preliminary extrapolation results for the flavour compositions at the physical point. A key result of this work is the observed mass splitting between the $π^0$ and $η$ on our ensembles, which is found to exhibit behaviour that is simply related to the corresponding flavour compositions.

preprint2020arXiv

Determining the glue component of the nucleon

Computing the gluon component of momentum in the nucleon is a difficult and computationally expensive problem, as the matrix element involves a quark-line-disconnected gluon operator which suffers from ultra-violet fluctuations. But also necessary for a successful determination is the non-perturbative renormalisation of this operator. As a first step we investigate here this renormalisation in the RI-MOM scheme. Using quenched QCD as an example, a statistical signal is obtained in a direct calculation using an adaption of the Feynman-Hellmann technique.

preprint2020arXiv

Electromagnetic contribution to $Σ$-$Λ$ mixing using lattice QCD+QED

Mixing in the $Σ^0$-$Λ^0$ system is a direct consequence of broken isospin symmetry and is a measure of both isospin-symmetry breaking as well as general SU(3)-flavour symmetry breaking. In this work we present a new scheme for calculating the extent of $Σ^0$-$Λ^0$ mixing using simulations in lattice QCD+QED and perform several extrapolations that compare well with various past determinations. Our scheme allows us to easily contrast the QCD-only mixing case with the full QCD+QED mixing.

preprint2020arXiv

Scaling and higher twist in the nucleon Compton amplitude

The partonic structure of hadrons plays an important role in a vast array of high-energy and nuclear physics experiments. It also underpins the theoretical understanding of hadron structure. Recent developments in lattice QCD offer new opportunities for reliably studying partonic structure from first principles. Here we report on the use of the Feynman-Hellmann theorem to study the forward Compton amplitude in the unphysical region. We demonstrate how this amplitude provides direct constraint on hadronic inelastic structure functions. The use of external momentum transfer allows us to study the $Q^2$ evolution to explore the onset of asymptotic scaling and reveal higher-twist effects in partonic structure.

preprint2019arXiv

Anomalous magnetic moment of the muon with dynamical QCD+QED

The current $3.5σ$ discrepancy between experimental and Standard Model determinations of the anomalous magnetic moment of the muon $a_μ=(g-2)/2$ can only be extended to the discovery $5σ$ regime through a reduction of both experimental and theoretical uncertainties. On the theory side, this means a determination of the hadronic vacuum polarisation (HVP) contribution to better than 0.5%, a level of precision that demands the inclusion of QCD + QED effects to properly understand how the behaviour of quarks are modified when their electric charges are turned on. The QCDSF collaboration has generated an ensemble of configurations with dynamical QCD and QED fields with the specific aim of studying flavour breaking effects arising from differences in the quark masses and charges in physical quantities. Here we study these effects in a calculation of HVP around the SU(3) symmetric point. Furthermore, by performing partially-quenched simulations we are able to cover a larger range of quark masses and charges on these configurations and then fit the results to an SU(3) flavour breaking expansion. Subsequently, this allows for an extrapolation to the physical point.

preprint2019arXiv

Patterns of flavour symmetry breaking in hadron matrix elements involving u, d and s quarks

By considering a flavour expansion about the SU(3)-flavour symmetric point, we investigate how flavour-blindness constrains octet baryon matrix elements after SU(3) is broken by the mass difference between quarks. Similarly to hadron masses we find the expansions to be constrained along a mass trajectory where the singlet quark mass is held constant, which provides invaluable insight into the mechanism of flavour symmetry breaking and proves beneficial for extrapolations to the physical point. Expansions are given up to third order in the expansion parameters. Considering higher orders would give no further constraints on the expansion parameters. The relation of the expansion coefficients to the quark-line-connected and quark-line disconnected terms in the 3-point correlation functions is also given. As we consider Wilson clover-like fermions, the addition of improvement coefficients is also discussed and shown to be included in the formalism developed here. As an example of the method we investigate this numerically via a lattice calculation of the flavour-conserving matrix elements of the vector first class form factors.

preprint2016arXiv

Isospin splittings of meson and baryon masses from three-flavor lattice QCD + QED

Lattice QCD simulations are now reaching a precision where isospin breaking effects become important. Previously, we have developed a program to systematically investigate the pattern of flavor symmetry beaking within QCD and successfully applied it to meson and baryon masses involving up, down and strange quarks. In this Letter we extend the calculations to QCD + QED and present our first results on isospin splittings in the pseudoscalar meson and baryon octets. In particular, we obtain the nucleon mass difference of $M_n-M_p=1.35(18)(8)\,\mbox{MeV}$ and the electromagnetic contribution to the pion splitting $M_{π^+}-M_{π^0}=4.60(20)\,\mbox{MeV}$. Further we report first determination of the separation between strong and electromagnetic contributions in the $\bar{MS}$ scheme.

preprint2016arXiv

Partially conserved axial vector current and applications

We investigate implications of the use of the point-split axial vector current derived from a Wilson like fermionic action. We compute the corresponding renormalization factor nonperturbatively for one beta value. The axial charge gA calculated from this nonlocal current is found to be nearer to the physical value than computed with the local axial vector current -- computed both on the same lattice with the same action.

preprint2016arXiv

QED effects in the pseudoscalar meson sector

We present results on the pseudoscalar meson masses from a fully dynamical simulation of QCD+QED. We concentrate particularly on violations of isospin symmetry. We calculate the $π^+$-$π^0$ splitting and also look at other isospin violating mass differences. We have presented results for these isospin splittings in arXiv:1508.06401 [hep-lat]. In this paper we give more details of the techniques employed, discussing in particular the question of how much of the symmetry violation is due to QCD, arising from the different masses of the $u$ and $d$ quarks, and how much is due to QED, arising from the different charges of the quarks. This decomposition is not unique, it depends on the renormalisation scheme and scale. We suggest a renormalisation scheme in which Dashen's theorem for neutral mesons holds, so that the electromagnetic self-energies of the neutral mesons are zero, and discuss how the self-energies change when we transform to a scheme such as $\bar{MS}$, in which Dashen's theorem for neutral mesons is violated.

preprint2015arXiv

A lattice determination of Sigma - Lambda mixing

Isospin breaking effects in baryon octet (and decuplet) masses are due to a combination of up and down quark mass differences and electromagnetic effects and lead to small mass splittings. Between the Sigma and Lambda this mass splitting is much larger, this being mostly due to their different wavefunctions. However when isospin is broken, there is a mixing between between these states. We describe the formalism necessary to determine the QCD mixing matrix and hence find the mixing angle and mass splitting between the Sigma and Lambda particles due to QCD effects.

preprint2015arXiv

Applications of the Feynman-Hellmann theorem in hadron structure

The Feynman-Hellmann (FH) relation offers an alternative way of accessing hadronic matrix elements through artificial modifications to the QCD Lagrangian. In particular, a FH-motivated method provides a new approach to calculations of disconnected contributions to matrix elements and high-momentum nucleon and pion form factors. Here we present results for the total nucleon axial charge, including a statistically significant non-negative total disconnected quark contribution of around $-5\%$ at an unphysically heavy pion mass. Extending the FH relation to finite-momentum transfers, we also present calculations of the pion and nucleon electromagnetic form factors up to momentum transfers of around 7-8 GeV$^2$. Results for the nucleon are not able to confirm the existence of a sign change for the ratio $\frac{G_E}{G_M}$, but suggest that future calculations at lighter pion masses will provide fascinating insight into this behaviour at large momentum transfers.

preprint2015arXiv

Charge Symmetry Violation in the Electromagnetic Form Factors of the Proton

Experimental tests of QCD through its predictions for the strange-quark content of the proton have been drastically restricted by our lack of knowledge of the violation of charge symmetry (CSV). We find unexpectedly tiny CSV in the proton's electromagnetic form factors by performing the first extraction of these quantities based on an analysis of lattice QCD data. The resulting values are an order of magnitude smaller than current bounds on proton strangeness from parity violating electron-proton scattering experiments. This result paves the way for a new generation of experimental measurements of the proton's strange form factors to challenge the predictions of QCD.

preprint2015arXiv

Determination of the strange nucleon form factors

The strange contribution to the electric and magnetic form factors of the nucleon is determined at a range of discrete values of $Q^2$ up to $1.4$ GeV$^2$. This is done by combining recent lattice QCD results for the electromagnetic form factors of the octet baryons with experimental determinations of those quantities. The most precise result is a small negative value for the strange magnetic moment: $G_M^s(Q^2=0) = -0.07\pm0.03\,μ_N$. At larger values of $Q^2$ both the electric and magnetic form factors are consistent with zero to within $2$-sigma.

preprint2015arXiv

Disconnected contributions to the spin of the nucleon

The spin decomposition of the proton is a long-standing topic of much interest in hadronic physics. Lattice QCD has had much success in calculating the connected contributions to the quark spin. However, complete calculations, which necessarily involve gluonic and strange-quark contributions, still present some challenges. These "disconnected" contributions typically involve small signals hidden against large statistical backgrounds and rely on computationally intensive stochastic techniques. In this work we demonstrate how a Feynman-Hellmann approach may be used to calculate such quantities, by measuring shifts in the proton energy arising from artificial modifications to the QCD action. We find a statistically significant non-zero result for the disconnected quark spin contribution to the proton of about -5% at a pion mass of 470 MeV.

preprint2015arXiv

Improving the lattice axial vector current

For Wilson and clover fermions traditional formulations of the axial vector current do not respect the continuum Ward identity which relates the divergence of that current to the pseudoscalar density. Here we propose to use a point-split or one-link axial vector current whose divergence exactly satisfies a lattice Ward identity, involving the pseudoscalar density and a number of irrelevant operators. We check in one-loop lattice perturbation theory with SLiNC fermion and gauge plaquette action that this is indeed the case including order $O(a)$ effects. Including these operators the axial Ward identity remains renormalisation invariant. First preliminary results of a nonperturbative check of the Ward identity are also presented.

preprint2015arXiv

Reply to "Comment on `Lattice determination of Sigma - Lambda mixing' "

In this Reply, we respond to the above Comment. Our computation [Phys. Rev. D 91 (2015) 074512] only took into account pure QCD effects, arising from quark mass differences, so it is not surprising that there are discrepancies in isospin splittings and in the Sigma - Lambda mixing angle. We expect that these discrepancies will be smaller in a full calculation incorporating QED effects.

preprint2015arXiv

SU(3) breaking in hyperon transition vector form factors

We present a calculation of the SU(3)-breaking corrections to the hyperon transition vector form factors to $\mathcal{O}(p^4)$ in heavy baryon chiral perturbation theory with finite-range regularisation. Both octet and decuplet degrees of freedom are included. We formulate a chiral expansion at the kinematic point $Q^2=-(M_{B_1}-M_{B_2})^2$, which can be conveniently accessed in lattice QCD. The two unknown low-energy constants at this point are constrained by lattice QCD simulation results for the $Σ^-\rightarrow n$ and $Ξ^0\rightarrow Σ^+$ transition form factors. Hence we determine lattice-informed values of $f_1$ at the physical point. This work constitutes progress towards the precise determination of $|V_{us}|$ from hyperon semileptonic decays.

preprint2015arXiv

The electric dipole moment of the neutron from 2+1 flavor lattice QCD

We compute the electric dipole moment d_n of the neutron from a fully dynamical simulation of lattice QCD with 2+1 flavors of clover fermions and nonvanishing theta term. The latter is rotated into the pseudoscalar density in the fermionic action using the axial anomaly. To make the action real, the vacuum angle theta is taken to be purely imaginary. The physical value of d_n is obtained by analytic continuation. We find d_n = -3.8(2)(9) x 10^{-16} [theta e cm], which, when combined with the experimental limit on d_n, leads to the upper bound theta < 7.6 x 10^{-11}.

preprint2015arXiv

Wilson flow and scale setting from lattice QCD

We give a determination of the phenomenological value of the Wilson (or gradient) flow scales t0 and w0 for 2+1 flavours of dynamical quarks. The simulations are performed keeping the average quark mass constant, which allows the approach to the physical point to be made in a controlled manner. O(a) improved clover fermions are used and together with four lattice spacings this allows the continuum extrapolation to be taken.

preprint2014arXiv

A Feynman-Hellmann approach to the spin structure of hadrons

We perform a Nf = 2 + 1 lattice QCD simulation to determine the quark spin fractions of hadrons using the Feynman-Hellmann theorem. By introducing an external spin operator to the fermion action, the matrix elements relevant for quark spin fractions are extracted from the linear response of the hadron energies. Simulations indicate that the Feynman-Hellmann method offers statistical precision that is comparable to the standard three-point function approach, with the added benefit that it is less susceptible to excited state contamination. This suggests that the Feynman-Hellmann technique offers a promising alternative for calculations of quark line disconnected contributions to hadronic matrix elements. At the SU(3)-flavour symmetry point, we find that the connected quark spin fractions are universally in the range 55-70% for vector mesons and octet and decuplet baryons. There is an indication that the amount of spin suppression is quite sensitive to the strength of SU(3) breaking.

preprint2014arXiv

A novel approach to nonperturbative renormalization of singlet and nonsinglet lattice operators

A novel method for nonperturbative renormalization of lattice operators is introduced, which lends itself to the calculation of renormalization factors for nonsinglet as well as singlet operators. The method is based on the Feynman-Hellmann relation, and involves computing two-point correlators in the presence of generalized background fields arising from introducing additional operators into the action. As a first application, and test of the method, we compute the renormalization factors of the axial vector current $A_μ$ and the scalar density $S$ for both nonsinglet and singlet operators for $N_f=3$ flavors of SLiNC fermions. For nonsinglet operators, where a meaningful comparison is possible, perfect agreement with recent calculations using standard three-point function techniques is found.

preprint2014arXiv

Connected and disconnected quark contributions to hadron spin

By introducing an external spin operator to the fermion action, the quark spin fractions of hadrons are determined from the linear response of the hadron energies using the Feynman-Hellmann (FH) theorem. At our SU(3)-flavour symmetric point, we find that the connected quark spin fractions are universally in the range 55-70\% for vector mesons and octet and decuplet baryons. There is an indication that the amount of spin suppression is quite sensitive to the strength of SU(3) breaking. We also present first preliminary results applying the FH technique to calculations of quark-line disconnected contributions to hadronic matrix elements of axial and tensor operators. At the SU(3)-flavour symmetric point we find a small negative contribution to the nucleon spin from disconnected quark diagrams, while the corresponding tensor matrix elements are consistent with zero.

preprint2014arXiv

Determining Sigma - Lambda mixing

SU2 isospin breaking effects in baryon octet (and decuplet) masses are due to a combination of up and down quark mass differences and electromagnetic effects. These mass differences are small. Between the Sigma and Lambda the mass splitting is much larger, but this is mostly due to their different wavefunctions. However there is now also mixing between these states. We determine the QCD mixing matrix and hence find the mixing angle and mass splitting.

preprint2014arXiv

Nucleon axial charge and pion decay constant from two-flavor lattice QCD

The axial charge of the nucleon $g_A$ and the pion decay constant $f_π$ are computed in two-flavor lattice QCD. The simulations are carried out on lattices of various volumes and lattice spacings. Results are reported for pion masses as low as $m_π=130\,\mbox{MeV}$. Both quantities, $g_A$ and $f_π$, suffer from large finite size effects, which to leading order ChEFT and ChPT turn out to be identical. By considering the naturally renormalized ratio $g_A/f_π$, we observe a universal behavior as a function of decreasing quark mass. From extrapolating the ratio to the physical point, we find $g_A^R=1.29(5)(3)$, using the physical value of $f_π$ as input and $r_0=0.50(1)$ to set the scale. In a subsequent calculation we attempt to extrapolate $g_A$ and $f_π$ separately to the infinite volume. Both volume and quark mass dependencies of $g_A$ and $f_π$ are found to be well decribed by ChEFT and ChPT. We find at the physical point $g_A^R=1.24(4)$ and $f_π^R=89.6(1.1)(1.8)\,\mbox{MeV}$. Both sets of results are in good agreement with experiment. As a by-product we obtain the low-energy constant $\bar{l}_4=4.2(1)$.

preprint2014arXiv

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

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

preprint2013arXiv

Electromagnetic splitting of quark and pseudoscalar meson masses from dynamical QCD + QED

Lattice QCD simulations are now reaching a precision where electromagnetic corrections from QED become important. In investigating the effects of SU(3) breaking due to quark mass differences within QCD, a group-theoretical analysis of the mass dependence greatly helped us organize our results. We now do the same with electromagnetic charge effects by extending the calculations to dynamical 1+1+1 flavor QCD + QED.

preprint2013arXiv

Perturbatively improving renormalization constants

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

preprint2013arXiv

Perturbatively improving RI-MOM renormalization constants

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

preprint2013arXiv

SU(3) flavour breaking and baryon structure

We present results from the QCDSF/UKQCD collaboration for hyperon electromagnetic form factors and axial charges obtained from simulations using Nf=2+1 flavours of O(a)-improved Wilson fermions. We also consider matrix elements relevant for hyperon semileptonic decays. We find flavour-breaking effects in hyperon magnetic moments which are consistent with experiment, while our results for the connected quark spin content indicates that quarks contribute more to the spin of the Xi baryon than they do to the proton.

preprint2013arXiv

SU(3) flavour symmetry breaking and charmed states

By extending the SU(3) flavour symmetry breaking expansion from up, down and strange sea quark masses to partially quenched valence quark masses we propose a method to determine charmed quark hadron masses including possible QCD isospin breaking effects. Initial results for some open charmed pseudoscalar meson states and singly and doubly charmed baryon states are encouraging and demonstrate the potential of the procedure. Essential for the method is the determination of the scale using singlet quantities, and to this end we also give here a preliminary estimation of the recently introduced Wilson flow scales.

preprint2012arXiv

Charge Symmetry Breaking in Spin Dependent Parton Distributions and the Bjorken Sum Rule

We present the first determination of charge symmetry violation (CSV) in the spin-dependent parton distribution functions of the nucleon. This is done by determining the first two Mellin moments of the spin-dependent parton distribution functions of the octet baryons from N_f = 2 + 1 lattice simulations. The results are compared with predictions from quark models of nucleon structure. We discuss the contribution of partonic spin CSV to the Bjorken sum rule, which is important because the CSV contributions represent the only partonic corrections to the Bjorken sum rule.

preprint2012arXiv

First moments of the nucleon generalized parton distributions from lattice QCD

We report on our lattice calculations of the nucleon's generalized parton distributions (GPDs), concentrating on their first moments for the case of N_f=2. Due to recent progress on the numerical side we are able to present results for the generalized form factors at pion masses as low as 260 MeV. We perform a fit to one-loop covariant baryon chiral perturbation theory with encouraging results.

preprint2012arXiv

Hyperon sigma terms for 2+1 quark flavours

QCD lattice simulations determine hadron masses as functions of the quark masses. From the gradients of these masses and using the Feynman-Hellmann theorem the hadron sigma terms can then be determined. We use here a novel approach of keeping the singlet quark mass constant in our simulations which upon using an SU(3) flavour symmetry breaking expansion gives highly constrained (i.e. few parameter) fits for hadron masses in a multiplet. This is a highly advantageous procedure for determining the hadron mass gradient as it avoids the use of delicate chiral perturbation theory. We illustrate the procedure here by estimating the light and strange sigma terms for the baryon octet.

preprint2012arXiv

Isospin breaking in octet baryon mass splittings

Using an SU(3) flavour symmetry breaking expansion in the quark mass, we determine the QCD component of the nucleon, Sigma and Xi mass splittings of the baryon octet due to up-down (and strange) quark mass differences in terms of the kaon mass splitting. Provided the average quark mass is kept constant, the expansion coefficients in our procedure can be determined from computationally cheaper simulations with mass degenerate sea quarks and partially quenched valence quarks. Both the linear and quadratic terms in the SU(3) flavour symmetry breaking expansion are considered; it is found that the quadratic terms only change the result by a few percent, indicating that the expansion is highly convergent.

preprint2012arXiv

Nucleon mass and sigma term from lattice QCD with two light fermion flavors

We analyze Nf=2 nucleon mass data with respect to their dependence on the pion mass down to mpi = 157 MeV and compare it with predictions from covariant baryon chiral perturbation theory (BChPT). A novel feature of our approach is that we fit the nucleon mass data simultaneously with the directly obtained pion-nucleon sigma-term. Our lattice data below mpi = 435 MeV is well described by O(p^4) BChPT and we find sigma=37(8)(6) MeV for the sigma-term at the physical point. Using the nucleon mass to set the scale we obtain a Sommer parameter of r_0=0.501(10)(11) fm.

preprint2012arXiv

Octet baryon mass splittings from up-down quark mass differences

Using an SU(3) flavour symmetry breaking expansion in the quark mass, we determine the QCD component of the neutron-proton, Sigma and Xi mass splittings of the baryon octet due to up-down (and strange) quark mass differences. Provided the average quark mass is kept constant, the expansion coefficients in our procedure can be determined from computationally cheaper simulations with mass degenerate sea quarks and partially quenched valence quarks. Full details and numerical results are given in ref 1.

preprint2012arXiv

Perturbative and Nonperturbative Renormalization in Lattice QCD

We investigate the perturbative and nonperturbative renormalization of composite operators in lattice QCD restricting ourselves to operators that are bilinear in the quark fields (quark-antiquark operators). These include operators which are relevant to the calculation of moments of hadronic structure functions. The nonperturbative computations are based on Monte Carlo simulations with two flavors of clover fermions and utilize the Rome-Southampton method also known as the RI-MOM scheme. We compare the results of this approach with various estimates from lattice perturbation theory, in particular with recent two-loop calculations.

preprint2012arXiv

Perturbative subtraction of lattice artifacts in the computation of renormalization constants

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

preprint2012arXiv

Scattering phases for meson and baryon resonances on general moving-frame lattices

A proposal by Lüscher enables one to compute the scattering phases of elastic two-body systems from the energy levels of the lattice Hamiltonian in a finite volume. In this work we generalize the formalism to S--, P-- and D--wave meson and baryon resonances, and general total momenta. Employing nonvanishing momenta has several advantages, among them making a wider range of energy levels accessible on a single lattice volume and shifting the level crossing to smaller values of $m_πL$.

preprint2012arXiv

Strangeness Contribution to the Proton Spin from Lattice QCD

We compute the strangeness and light-quark contributions Delta s, Delta u and Delta d to the proton spin in nf=2 lattice QCD at a pion mass of about 285 MeV and at a lattice spacing a approx 0.073 fm, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson action. We carry out the renormalization of these matrix elements which involves mixing between contributions from different quark flavours. Our main result is the small negative value MSbar scheme value Delta s = -0.020(10)(4), at a scale mu = 7.4 GeV, of the strangeness contribution to the nucleon spin. The second error is an estimate of the uncertainty, due to the missing extrapolation to the physical point.

preprint2012arXiv

The effects of flavour symmetry breaking on hadron matrix elements

By considering a flavour expansion about the SU(3)-flavour symmetric point, we investigate how flavour-blindness constrains octet baryon matrix elements after SU(3) is broken by the mass difference between the strange and light quarks. We find the expansions to be highly constrained along a mass trajectory where the singlet quark mass is held constant, which proves beneficial for extrapolations of 2+1 flavour lattice data to the physical point. We investigate these effects numerically via a lattice calculation of the flavour-conserving and flavour-changing matrix elements of the vector and axial operators between octet baryon states.

preprint2012arXiv

The strange and light quark contributions to the nucleon mass from Lattice QCD

We determine the strangeness and light quark fractions of the nucleon mass by computing the quark line connected and disconnected contributions to the matrix elements m_q <N|qbar q|N> in lattice QCD, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson Fermionic action. We simulate n_F=2 mass degenerate sea quarks with a pion mass of about 285 MeV and a lattice spacing a approx 0.073 fm. The renormalization of the matrix elements involves mixing between contributions from different quark flavours. The pion-nucleon sigma-term is extrapolated to physical quark masses exploiting the sea quark mass dependence of the nucleon mass. We obtain the renormalized values σ_{πN} = 38(12) MeV at the physical point and f_{T_s}=σ_s/m_N= 0.012(14)^{+10}_{-3} for the strangeness contribution at our larger than physical sea quark mass.

preprint2012arXiv

Wilson loops to 20th order numerical stochastic perturbation theory

We calculate Wilson loops of various sizes up to 20 loops in SU(3) pure lattice gauge theory at different lattice sizes for Wilson gauge action using the technique of numerical stochastic perturbation theory. This allows us to investigate the perturbative series for various Wilson loops at high loop orders. We observe differences in the behavior of those series as function of the loop order. Up to $n=20$ we do not find evidence for the factorial growth of the expansion coefficients often assumed to characterize an asymptotic series. Based on the actually observed behavior we sum the series in a model parametrized by hypergeometric functions. Alternatively we estimate the total series in boosted perturbation theory using information from the first 14 loops. We introduce generalized ratios of Wilson loops of different sizes. Together with the corresponding Wilson loops from standard Monte Carlo measurements they enable us to assess their non-perturbative parts.

preprint2011arXiv

A lattice study of the strangeness content of the nucleon

We determine the quark contributions to the nucleon spin Delta s, Delta u and Delta d as well as their contributions to the nucleon mass, the sigma-terms. This is done by computing both, the quark line connected and disconnected contributions to the respective matrix elements, using the non-perturbatively improved Sheikholeslami-Wohlert Wilson Fermionic action. We simulate n_F=2 mass degenerate sea quarks with a pion mass of about 285 MeV and a lattice spacing a = 0.073 fm. The renormalization of the matrix elements involves mixing between contributions from different quark flavours. The pion-nucleon sigma-term is extrapolated to physical quark masses exploiting the sea quark mass dependence of the nucleon mass. We obtain the renormalized value sigma_{piN}=38(12) MeV at the physical point and the strangeness fraction f_{Ts}=sigma_s/m_N=0.012(14)(+10-3) at our larger than physical sea quark mass. For the strangeness contribution to the nucleon spin we obtain in the MSbar scheme at the renormalization scale of 2.71 GeV Delta s = -0.020(10)(2).

preprint2011arXiv

An Update on Distribution Amplitudes of the Nucleon and its Parity Partner

The calculation of baryon wave functions at small inter-quark separations is an ongoing effort within the QCDSF collaboration. In this update on normalization constants and distribution amplitudes of the nucleon and its negative parity partner, N* (1535), we present new lattice data which helps us controlling finite size effects. We use new chiral perturbation theory results to perform the extrapolation to the physical point.

preprint2011arXiv

Baryon Axial Charges and Momentum Fractions with N_f=2+1 Dynamical Fermions

We report on recent results of the QCDSF/UKQCD Collaboration on investigations of baryon structure using configurations generated with N_f=2+1 dynamical flavours of O(a) improved Wilson fermions. With the strange quark mass as an additional dynamical degree of freedom in our simulations we avoid the need for a partially quenched approximation when investigating the properties of particles containing a strange quark, e.g. the hyperons. In particular, we will focus on the nucleon and hyperon axial charges and quark momentum fractions.

preprint2011arXiv

Charge Symmetry Breaking in Parton Distribution Functions from Lattice QCD

By determining the quark momentum fractions of the octet baryons from N_f=2+1 lattice simulations, we are able to predict the degree of charge symmetry violation in the parton distribution functions of the nucleon. This is of importance, not only as a probe of our understanding of the non-perturbative structure of the proton but also because such a violation constrains the accuracy of global fits to parton distribution functions and hence the accuracy with which, for example, cross sections at the LHC can be predicted. A violation of charge symmetry may also be critical in cases where symmetries are used to guide the search for physics beyond the Standard Model.

preprint2011arXiv

Dirac and Pauli form factors from lattice QCD

We present a comprehensive analysis of the electromagnetic form factors of the nucleon from a lattice simulation with two flavors of dynamical O(a)-improved Wilson fermions. A key feature of our calculation is that we make use of an extensive ensemble of lattice gauge field configurations with four different lattice spacings, multiple volumes, and pion masses down to m_π~ 180 MeV. We find that by employing Kelly-inspired parametrizations for the Q^2-dependence of the form factors, we are able to obtain stable fits over our complete ensemble. Dirac and Pauli radii and the anomalous magnetic moments of the nucleon are extracted and results at light quark masses provide evidence for chiral non-analytic behavior in these fundamental observables.

preprint2011arXiv

Disconnected contributions to D-meson semi-leptonic decay form factors

We calculate the disconnected contribution to the form factor for the semileptonic decay of a D-meson into a final state, containing a flavor singlet eta meson. We use QCDSF n_f=2+1 configurations at the flavor symmetric point m_u=m_d=m_s and the partially quenched approximation for the relativistic charm quark. Several acceleration and noise reduction techniques for the stochastic estimation of the disconnected loop are tested.

preprint2011arXiv

Flavour blindness and patterns of flavour symmetry breaking in lattice simulations of up, down and strange quarks

QCD lattice simulations with 2+1 flavours (when two quark flavours are mass degenerate) typically start at rather large up-down and strange quark masses and extrapolate first the strange quark mass and then the up-down quark mass to its respective physical value. Here we discuss an alternative method of tuning the quark masses, in which the singlet quark mass is kept fixed. Using group theory the possible quark mass polynomials for a Taylor expansion about the flavour symmetric line are found, first for the general 1+1+1 flavour case and then for the 2+1 flavour case. This ensures that the kaon always has mass less than the physical kaon mass. This method of tuning quark masses then enables highly constrained polynomial fits to be used in the extrapolation of hadron masses to their physical values. Numerical results for the 2+1 flavour case confirm the usefulness of this expansion and an extrapolation to the physical pion mass gives hadron mass values to within a few percent of their experimental values. Singlet quantities remain constant which allows the lattice spacing to be determined from hadron masses (without necessarily being at the physical point). Furthermore an extension of this programme to include partially quenched results is given.

preprint2011arXiv

Hyperon Form Factors from N_f=2+1 QCD

We present results from the QCDSF/UKQCD collaboration for the electromagnetic and semi-leptonic form factors for the hyperons. The simulations are performed on our new ensembles generated with 2+1 flavours of dynamical O(a)-improved Wilson fermions. A unique feature of these configurations is that the quark masses are tuned so that the singlet quark mass is held fixed at its physical value. We use 5 such choices of the individual quark masses on 24^3x48 lattices with a lattice spacing of about 0.078 fm.

preprint2011arXiv

Nucleon form factors and structure functions

We give an update on our ongoing efforts to compute the nucleon's form factors and moments of structure functions using Nf=2 flavours of non-perturbatively improved Clover fermions. We focus on new results obtained on gauge configurations where the pseudo-scalar meson mass is in the range of 170-270 MeV. We will compare our results with various estimates obtained from chiral effective theories since we have some overlap with the quark mass region where results from such theories are believed to be applicable.

preprint2011arXiv

Nucleon sigma terms for 2+1 quark flavours

QCD lattice simulations yield hadron masses as functions of the quark masses. From the gradients of the hadron masses the sigma terms can then be determined. We consider here dynamical 2+1 flavour simulations, in which we start from a point of the flavour symmetric line and then keep the singlet or average quark mass fixed as we approach the physical point. This leads to highly constrained fits for hadron masses in a multiplet. The gradient of this path for a hadron mass then gives a relation between the light and strange sigma terms. A further relation can be found from the change in the singlet quark mass along the flavour symmetric line. This enables light and strange sigma terms to be estimated for the baryon octet.

preprint2011arXiv

QCD in the delta-Regime

The delta-regime of QCD is characterised by light quarks in a small spatial box, but a large extent in (Euclidean) time. In this setting a specific variant of chiral perturbation theory - the delta-expansion - applies, based on a quantum mechanical treatment of the quasi one-dimensional system. In particular, for vanishing quark masses one obtains a residual pion mass M_pi^R, which has been computed to the third order in the delta-expansion. A comparison with numerical measurements of this residual mass allows for a new determination of some Low Energy Constants, which appear in the chiral Lagrangian. We first review the attempts to simulate 2-flavour QCD directly in the delta-regime. This is very tedious, but results compatible with the predictions for M_pi^R have been obtained. Then we show that an extrapolation of pion masses measured in a larger volume towards the delta-regime leads to good agreement with the theoretical predictions. From those results, we also extract a value for the (controversial) sub-leading Low Energy Constant \bar l_3.

preprint2010arXiv

Exploring the Nucleon Structure from First Principles of QCD

Quantum Chromodynamics (QCD) is generally assumed to be the fundamental theory underlying nuclear physics. In recent years there is progress towards investigating the nucleon structure from first principles of QCD. Although this structure is best revealed in Deep Inelastic Scattering, a consistent analysis has to be performed in a fully non-perturbative scheme. The only known method for this purpose are lattice simulations. We first sketch the ideas of Monte Carlo simulations in lattice gauge theory. Then we comment in particular on the issues of chiral symmetry and operator mixing. Finally we present our results for the Bjorken variable of a single quark, and for the second Nachtmann moment of the nucleon structure functions.

preprint2010arXiv

Flavour symmetry breaking and tuning the strange quark mass for 2+1 quark flavours

QCD lattice simulations with 2+1 flavours typically start at rather large up-down and strange quark masses and extrapolate first the strange quark mass to its physical value and then the up-down quark mass. An alternative method of tuning the quark masses is discussed here in which the singlet quark mass is kept fixed, which ensures that the kaon always has mass less than the physical kaon mass. Using group theory the possible quark mass polynomials for a Taylor expansion about the flavour symmetric line are found, which enables highly constrained fits to be used in the extrapolation of hadrons to the physical pion mass. Numerical results confirm the usefulness of this expansion and an extrapolation to the physical pion mass gives hadron mass values to within a few percent of their experimental values.

preprint2010arXiv

Nucleon and $N^* (1535)$ Distribution Amplitudes

The QCDSF collaboration has investigated the distribution amplitudes and wavefunction normalization constants of the nucleon and its parity partner, the $N^* (1535)$. We report on recent progress in the calculation of these quantities on configurations with two dynamical flavors of $\mathcal{O}(a)$-improved Wilson fermions. New data at pion masses of approximately 270 MeV helps in significantly reducing errors in the extrapolation to the physical point.

preprint2010arXiv

Probing the finite temperature phase transition with Nf=2 nonperturbatively improved Wilson fermions

The critical temperature and the nature of the QCD finite temperature phase transition are determined for Nf=2 dynamical flavors of nonperturbatively improved Wilson fermions. The calculations are performed on large lattices with temporal extents Nt=12, 10 and 8, and lattice spacings down to a = 0.075 fm. We find the deconfinement and chiral phasetransitions to take place at the same temperature. Our results are in broad agreement with a second order phase transition in the chiral limit. The critical temperature at the physical quark mass is found to be Tc = 174(3)(6) MeV.

preprint2010arXiv

Renormalisation of composite operators in lattice QCD: perturbative versus nonperturbative

The perturbative and nonperturbative renormalisation of quark-antiquark operators in lattice QCD with two flavours of clover fermions is investigated within the research programme of the QCDSF collaboration. Operators with up to three derivatives are considered. The nonperturbative results based on the RI-MOM scheme are compared with estimates from one- and two-loop lattice perturbation theory.

preprint2010arXiv

Tuning the strange quark mass in lattice simulations

QCD lattice simulations with 2+1 flavours typically start at rather large up-down and strange quark masses and extrapolate first the strange quark mass to its physical value and then the up-down quark mass. An alternative method of tuning the quark masses is discussed here in which the singlet quark mass is kept fixed, which ensures that the kaon always has mass less than the physical kaon mass. It can also take into account the different renormalisations (for singlet and non-singlet quark masses) occurring for non-chirally invariant lattice fermions and so allows a smooth extrapolation to the physical quark masses. This procedure enables a wide range of quark masses to be probed, including the case with a heavy up-down quark mass and light strange quark mass. Results show the correct order for the baryon octet and decuplet spectrum and an extrapolation to the physical pion mass gives mass values to within a few percent of their experimental values.

preprint2010arXiv

Very high order lattice perturbation theory for Wilson loops

We calculate perturbative Wilson loops of various sizes up to loop order $n=20$ at different lattice sizes for pure plaquette and tree-level improved Symanzik gauge theories using the technique of Numerical Stochastic Perturbation Theory. This allows us to investigate the behavior of the perturbative series at high orders. We observe differences in the behavior of perturbative coefficients as a function of the loop order. Up to $n=20$ we do not see evidence for the often assumed factorial growth of the coefficients. Based on the observed behavior we sum this series in a model with hypergeometric functions. Alternatively we estimate the series in boosted perturbation theory. Subtracting the estimated perturbative series for the average plaquette from the non-perturbative Monte Carlo result we estimate the gluon condensate.

preprint2009arXiv

A Non-Perturbative Operator Product Expansion

Nucleon structure functions can be observed in Deep Inelastic Scattering experiments, but it is an outstanding challenge to confront them with fully non-perturbative QCD results. For this purpose we investigate the product of electromagnetic currents (with large photon momenta) between quark states (of low momenta). By means of an Operator Product Expansion the structure function can be decomposed into matrix elements of local operators, and Wilson coefficients. For consistency both have to be computed non-perturbatively. Here we present precision results for a set of Wilson coefficients. They are evaluated from propagators for numerous quark momenta on the lattice, where the use of chiral fermions suppresses undesired operator mixing. This over-determines the Wilson coefficients, but reliable results can be extracted by means of a Singular Value Decomposition.

preprint2009arXiv

Lattice Investigations of Nucleon Structure at Light Quark Masses

Lattice simulations of hadronic structure are now reaching a level where they are able to not only complement, but also provide guidance to current and forthcoming experimental programmes at, e.g. Jefferson Lab, COMPASS/CERN and FAIR/GSI. By considering new simulations at low quark masses and on large volumes, we review the recent progress that has been made in this exciting area by the QCDSF/UKQCD collaboration. In particular, results obtained close to the physical point for several quantities, including electromagnetic form factors and moments of ordinary parton distribution functions, show some indication of approaching their phenomenological values.

preprint2009arXiv

Results from 2+1 flavours of SLiNC fermions

QCD results are presented for a 2+1 flavour fermion clover action (which we call the SLiNC action). A method of tuning the quark masses to their physical values is discussed. In this method the singlet quark mass is kept fixed, which solves the problem of different renormalisations (for singlet and non-singlet quark masses) occuring for non-chirally invariant lattice fermions. This procedure enables a wide range of quark masses to be probed, including the case with a heavy up-down quark mass and light strange quark mass. Preliminary results show the correct splittings for the baryon (octet and) decuplet spectrum.

preprint2009arXiv

Wilson loops in very high order lattice perturbation theory

We calculate Wilson loops of various sizes up to loop order $n=20$ for lattice sizes of $L^4 (L=4, 6, 8, 12)$ using the technique of Numerical Stochastic Perturbation Theory in quenched QCD. This allows to investigate the behaviour of the perturbative series at high orders. We discuss three models to estimate the perturbative series: a renormalon inspired fit, a heuristic fit based on an assumed power-law singularity and boosted perturbation theory. We have found differences in the behavior of the perturbative series for smaller and larger Wilson loops at moderate $n$. A factorial growth of the coefficients could not be confirmed up to $n=20$. From Monte Carlo measured plaquette data and our perturbative result we estimate a value of the gluon condensate $<\fracαπGG>$.

preprint2008arXiv

Clover improvement for stout-smeared 2+1 flavour SLiNC fermions: non-perturbative results

We discuss an action in which the fermion matrix has single level stout smearing for the hopping terms together with unsmeared links for the clover term. With the (tree level) Symanzik improved gluon action this constitutes the Stout Link Non-perturbative Clover or SLiNC action. To cancel O(a) terms the clover coefficient, csw, has to be tuned. We present here preliminary results of a non-perturbative determination of csw using the Schrodinger functional and as a by-product also a determination of the critical hopping parameter. A determination of the renormalisation constant for the local vector current is also given. Comparisons of the results are made with lowest order perturbation theory results.

preprint2008arXiv

Clover improvement for stout-smeared 2+1 flavour SLiNC fermions: perturbative results

For the Stout Link Non-perturbative Clover (SLiNC) action we determine in one-loop lattice perturbation theory the critical hopping parameter $κ_c$ and the clover parameter $c_{SW}$ which is needed for $\mathcal{O}(a)$ improvement. Performing this calculation off-shell we are also able to compute the non gauge invariant quark field improvement coefficient $c_{NGI}$. Additionally, we present first results for the renormalization factors of the scalar, pseudoscalar, vector and axial vector currents. We discuss mean field improvement for the SLiNC action.

preprint2008arXiv

Nucleon distribution amplitudes and proton decay matrix elements on the lattice

Baryon distribution amplitudes (DAs) are crucial for the theory of hard exclusive reactions. We present a calculation of the first few moments of the leading-twist nucleon DA within lattice QCD. In addition we deal with the normalization of the next-to-leading (twist-four) DAs. The matrix elements determining the latter quantities are also responsible for proton decay in Grand Unified Theories. Our lattice evaluation makes use of gauge field configurations generated with two flavors of clover fermions. The relevant operators are renormalized nonperturbatively with the final results given in the MSbar scheme. We find that the deviation of the leading-twist nucleon DA from its asymptotic form is less pronounced than sometimes claimed in the literature.

preprint2008arXiv

Renormalisation of Composite Operators in Lattice Perturbation Theory with Clover Fermions: Non-forward Matrix Elements

We consider the renormalisation of lattice QCD operators with one and two covariant derivatives related to the first and second moments of generalised parton distributions and meson distribution amplitudes. Employing the clover fermion action we calculate their non-forward quark matrix elements in one-loop lattice perturbation theory. For some representations of the hypercubic group commonly used in simulations we determine the sets of all possible mixing operators and compute the matrices of renormalisation factors in one-loop approximation. We describe how tadpole improvement is applied to the results.

preprint2007arXiv

Moments of generalized parton distributions and quark angular momentum of the nucleon

The internal structure of hadrons is important for a variety of topics, including the hadron form factors, proton spin and spin asymmetry in polarized proton scattering. For a systematic study generalized parton distributions (GPDs) encode important information on hadron structure in the entire impact parameter space. We report on a computation of nucleon GPDs based on simulations with two dynamical non-perturbatively improved Wilson quarks with pion masses down to 350MeV. We present results for the total angular momentum of quarks with chiral extrapolation based on covariant baryon chiral perturbation theory.

preprint2005arXiv

A determination of the strange quark mass for unquenched clover fermions using the AWI

Using the O(a) Symanzik improved action an estimate is given for the strange quark mass for unquenched (nf=2) QCD. The determination is via the axial Ward identity (AWI) and includes a non-perturbative evaluation of the renormalisation constant. Numerical results have been obtained at several lattice spacings, enabling the continuum limit to be taken. Results indicate a value for the strange quark mass (in the MSbar-scheme at a scale of 2GeV) in the range 100 - 130MeV.

preprint2005arXiv

Meson decay constants from Nf=2 clover fermions

We present recent results for meson decay constants calculated on configurations with two flavours of O(a)-improved Wilson fermions. Non-perturbative renormalisation is applied and quark mass dependencies as well as finite volume and discretisation effects are investigated. In this work we also present a computation of the coupling of the light vector mesons to the tensor current using dynamical fermions.

preprint2005arXiv

Non-perturbative renormalisation for overlap fermions

Using non-perturbative techniques we have found the renormalisation factor, Z, in the RI-MOM scheme for quark bilinear operators in quenched QCD. We worked with overlap fermions using the Luescher-Weisz gauge action. Our calculation was performed at beta=8.45 at a lattice spacing of 1/a=2.1 GeV using a value of rho=1.4. Our results show good agreement between the vector and the axial vector in the zero mass limit. This shows that overlap fermions have good chiral properties. To attempt to improve the discretisation errors in our results we subtracted the O(a^2) terms in one-loop lattice perturbation theory from the Monte Carlo Green functions. In particular we paid attention to the operators for the observable <x>. We found a value for the renormalisation constants Z^msbar_(v_2b) and Z^msbar_(v_2a) just less than 1.9 at mu=1/a=2.1 GeV.

preprint2005arXiv

Perturbative renormalisation of quark bilinear operators for overlap fermions with and without stout links and improved gauge action

We calculate lattice renormalisation constants of local and one-link quark operators for overlap fermions and improved gauge actions in one-loop perturbation theory. For the local operators we stout smear the SU(3) links in the fermionic action. Using the popular tadpole improved Lüscher-Weisz actions at $β=8.45$ and $β=8.0$ we present numerical values for the Z factors in the $\bar{MS}$ scheme (partly as function of the stout smearing strength). We compare various levels of mean field (tadpole) improvement which have been applied to our results.

preprint2004arXiv

The nucleon mass in N_f=2 lattice QCD: finite size effects from chiral perturbation theory

In the framework of relativistic SU(2)_f baryon chiral perturbation theory we calculate the volume dependence of the nucleon mass up to and including O(p^4). Since the parameters in the resulting finite size formulae are fixed from the pion mass dependence of the large volume nucleon masses and from phenomenology, we obtain a parameter-free prediction of the finite size effects. We present mass data from the recent N_f=2 simulations of the UKQCD and QCDSF collaborations and compare these data as well as published mass values from the dynamical simulations of the CP-PACS and JLQCD collaborations with the theoretical expectations. Remarkable agreement between the lattice data and the predictions of chiral perturbation theory in a finite volume is found.

preprint2002arXiv

Excited nucleon spectrum using a non-perturbatively improved clover fermion action

We discuss the extraction of negative-parity baryon masses from lattice QCD calculations. The mass of the lowest-lying negative-parity $J = 1/2^{-}$ state is computed in quenched lattice QCD using an ${\cal O}(a)$-improved clover fermion action, and a splitting found with the nucleon mass. The calculation is performed on two lattice volumes, and three lattice spacings enabling a study of both finite-volume and finite-lattice-spacing uncertainties. A measurement of the first excited radial excitation of the nucleon finds a mass considerably larger than that of the negative-parity ground state, in accord with other lattice determinations but in disagreement with experiment. Results are also presented for the lightest negative-parity $I=3/2$ state.

preprint2002arXiv

Lattice Calculation of Baryon Masses using the Clover Fermion Action

We present a calculation of the lowest-lying baryon masses in the quenched approximation to QCD. The calculations are performed using a non-perturbatively improved clover fermion action, and a splitting is found between the masses of the nucleon and its parity partner. An analysis of the mass of the first radial excitation of the nucleon finds a value considerably larger than that of the parity partner of the nucleon, and thus little evidence for the Roper resonance as a simple three-quark state

preprint2002arXiv

Negative-parity Baryon Masses using an O(a)-improved Fermion Action

We present a calculation of the mass of the lowest-lying negative-parity J=1/2- state in quenched QCD. Results are obtained using a non-perturbatively O(a)-improved clover fermion action, and a splitting is found between the mass of the nucleon and its parity partner. The calculation is performed on two lattice volumes and at three lattice spacings, enabling a study of both finite-volume and finite lattice-spacing uncertainties. A comparison is made with results obtained using the unimproved Wilson fermion action.

preprint1999arXiv

Beyond the Thouless energy

The distribution and the correlations of the small eigenvalues of the Dirac operator are described by random matrix theory (RMT) up to the Thouless energy $E_c\propto 1/\sqrt{V}$, where $V$ is the physical volume. For somewhat larger energies, the same quantities can be described by chiral perturbation theory (chPT). For most quantities there is an intermediate energy regime, roughly $1/V<E<1/\sqrt{V}$, where the results of RMT and chPT agree with each other. We test these predictions by constructing the connected and disconnected scalar susceptibilities from Dirac spectra obtained in quenched SU(2) and SU(3) simulations with staggered fermions for a variety of lattice sizes and coupling constants. In deriving the predictions of chPT, it is important to take into account only those symmetries which are exactly realized on the lattice.