Researcher profile

A. D. Kennedy

A. D. Kennedy contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2012arXiv

Algorithms for Dynamical Fermions

This is the write-up of three lectures on algorithms for dynamical fermions that were given at the ILFTN workshop 'Perspectives in Lattice QCD' in Nara during November 2005. The first lecture is on the fundamentals of Markov Chain Monte Carlo methods and introduces the Hybrid Monte Carlo (HMC) algorithm and symplectic integrators; the second lecture covers topics in approximation theory and thereby introduces the Rational Hybrid Monte Carlo (RHMC) algorithm and ways of evading integrator instabilities by means of multiple pseudofermion fields; the third lecture introduces on-shell chiral (Ginsparg-Wilson) lattice fermions and discusses five-dimensional formulations for computing fermion propagators for such fermions.

preprint2011arXiv

Bounds on the convergence of Ritz values from Krylov subspaces to interior eigenvalues of Hermitean matrices

We consider bounds on the convergence of Ritz values from a sequence of Krylov subspaces to interior eigenvalues of Hermitean matrices. These bounds are useful in regions of low spectral density, for example near voids in the spectrum, as is required in many applications. Our bounds are obtained by considering the usual Kaniel-Paige-Saad formalism applied to the shifted and squared matrix.

preprint2010arXiv

Time symmetric Bloch-Nordsieck infrared cancellation

It is usually assumed that the infrared divergences that occur for virtual photons in Feynman diagrams cancel against the corresponding divergences in phase space integrals for the emission of real soft photons. If this were the case then why do the corresponding phase space integrals for the absorption of real soft photons not diverge too, as one might expect on grounds of time-reversal symmetry, and if so what cancels these divergences? We show that this apparent time asymmetry is an artifact of the way that the diagrams are calculated, and we introduce a procedure for such calculations that is manifestly time symmetric. Moreover, we show that it is not mathematically meaningful to distinguish arbitrarily soft massless particles of positive and negative energy.

preprint2009arXiv

Force Gradient Integrators

We present initial results of the use of Force Gradient integrators for lattice field theories. These promise to give significant performance improvements, especially for light fermions and large lattices. Our results show that this is indeed the case, indicating a speed-up of more than a factor of two, which is expected to increase as the integration step size becomes smaller for larger lattices and smaller fermion masses.

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.

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

Tuning HMC using Poisson brackets

We discuss how the integrators used for the Hybrid Monte Carlo (HMC) algorithm not only approximately conserve some Hamiltonian $H$ but exactly conserve a nearby shadow Hamiltonian (\tilde H), and how the difference $ΔH \equiv \tilde H - H $ may be expressed as an expansion in Poisson brackets. By measuring average values of these Poisson brackets over the equilibrium distribution $\propto e^{-H}$ generated by HMC we can find the optimal integrator parameters from a single simulation. We show that a good way of doing this in practice is to minimize the variance of $ΔH$ rather than its magnitude, as has been previously suggested. Some details of how to compute Poisson brackets for gauge and fermion fields, and for nested and force gradient integrators are also presented.

preprint2007arXiv

2+1 flavor domain wall QCD on a (2 fm)^3 lattice: light meson spectroscopy with Ls = 16

We present results for light meson masses and pseudoscalar decay constants from the first of a series of lattice calculations with 2+1 dynamical flavors of domain wall fermions and the Iwasaki gauge action. The work reported here was done at a fixed lattice spacing of about 0.12 fm on a 16^3\times32 lattice, which amounts to a spatial volume of (2 fm)^3 in physical units. The number of sites in the fifth dimension is 16, which gives m_{res} = 0.00308(4) in these simulations. Three values of input light sea quark masses, m_l^{sea} \approx 0.85 m_s, 0.59 m_s and 0.33 m_s were used to allow for extrapolations to the physical light quark limit, whilst the heavier sea quark mass was fixed to approximately the physical strange quark mass m_s. The exact rational hybrid Monte Carlo algorithm was used to evaluate the fractional powers of the fermion determinants in the ensemble generation. We have found that f_π= 127(4) MeV, f_K = 157(5) MeV and f_K/f_π= 1.24(2), where the errors are statistical only, which are in good agreement with the experimental values.

preprint1994arXiv

The heavy quark potential in QCD with 2 flavors of dynamical quarks

We compute the heavy quark potential on configurations generated by the HEMCGC collaboration with dynamical staggered fermions at $6/g^2 = 5.6$ and with dynamical Wilson fermions at $6/g^2 = 5.3$. The computations are done on $16^3 \times 32$ lattices, corresponding to physical sizes of about 1.6 and 2.3 fm, respectively. Up to the distances probed no sign of string breaking is detectable. We also compute the recently proposed scale $r_0$ defined by $r_0^2 F(r_0) = 1.65$.