James C. Osborn

Lattice QCD ensembles with four flavors of highly improved staggered quarks

We present results from our simulations of quantum chromodynamics with four flavors of quarks: u, d, s, and c. These simulations are performed with a one-loop Symanzik improved gauge action and the highly improved staggered quark action. We are generating gauge configurations with four values of the lattice spacing ranging from 0.06 to 0.15 fm, and three values of the light quark mass, including the value for which the Goldstone pion mass is equal to the physical pion mass. We discuss simulation algorithms, scale setting, taste symmetry breaking, and the autocorrelations of various quantities. We also present results for the topological susceptibility that demonstrate the improvement of the highly improved staggered quark configurations relative to those generated earlier with the asqtad improved staggered action.

Scaling studies of QCD with the dynamical highly improved staggered quark action

We study the lattice spacing dependence, or scaling, of physical quantities using the highly improved staggered quark (HISQ) action introduced by the HPQCD/UKQCD collaboration, comparing our results to similar simulations with the asqtad fermion action. Results are based on calculations with lattice spacings approximately 0.15, 0.12 and 0.09 fm, using four flavors of dynamical HISQ quarks. The strange and charm quark masses are near their physical values, and the light-quark mass is set to 0.2 times the strange-quark mass. We look at the lattice spacing dependence of hadron masses, pseudoscalar meson decay constants, and the topological susceptibility. In addition to the commonly used determination of the lattice spacing through the static quark potential, we examine a determination proposed by the HPQCD collaboration that uses the decay constant of a fictitious “unmixed ss̄” pseudoscalar meson. We find that the lattice artifacts in the HISQ simulations are much smaller than those in the asqtad simulations at the same lattice spacings and quark masses. PACS numbers: 12.38.Gc,14.20.Dh

Exploring strange nucleon form factors on the lattice

We discuss techniques for evaluating sea quark contributions to hadronic form factors on the lattice and apply these to an exploratory calculation of the strange electromagnetic, axial, and scalar form factors of the nucleon. We employ the Wilson gauge and fermion actions on an anisotropic 24^3 x 64 lattice, probing a range of momentum transfer with Q^2 _0. We discuss the unique systematic uncertainties affecting the latter quantity relative to the continuum, as well as prospects for improving future determinations with Wilson-like fermions.

Adaptive Multigrid Algorithm for the Lattice Wilson-Dirac Operator

We present an adaptive multigrid solver for application to the non-Hermitian Wilson-Dirac system of QCD. The key components leading to the success of our proposed algorithm are the use of an adaptive projection onto coarse grids that preserves the near null space of the system matrix together with a simplified form of the correction based on the so-called γ5-Hermitian symmetry of the Dirac operator. We demonstrate that the algorithm nearly eliminates critical slowing down in the chiral limit and that it has weak dependence on the lattice volume.

Chiral symmetry breaking and the dirac spectrum at nonzero chemical potential

The relation between the spectral density of the QCD Dirac operator at nonzero baryon chemical potential and the chiral condensate is investigated. We use the analytical result for the eigenvalue density in the microscopic regime which shows oscillations with a period that scales as 1/V and an amplitude that diverges exponentially with the volume V = L4. We find that the discontinuity of the chiral condensate is due to the whole oscillating region rather than to an accumulation of eigenvalues at the origin. These results also extend beyond the microscopic regime to chemical potentials mu approximately 1/L.

Unquenched QCD Dirac Operator Spectra at Nonzero Baryon Chemical Potential

Abstract The microscopic spectral density of the QCD Dirac operator at nonzero baryon chemical potential for an arbitrary number of quark flavors was derived recently from a random matrix model with the global symmetries of QCD. In this paper we show that these results and extensions thereof can be obtained from the replica limit of a Toda lattice equation. This naturally leads to a factorized form into bosonic and fermionic QCD-like partition functions. In the microscopic limit these partition functions are given by the static limit of a chiral Lagrangian that follows from the symmetry breaking pattern. In particular, we elucidate the role of the singularity of the bosonic partition function in the orthogonal polynomials approach. A detailed discussion of the spectral density for one and two flavors is given.

Adaptive multigrid algorithm for lattice QCD.

We present a new multigrid solver that is suitable for the Dirac operator in the presence of disordered gauge fields. The key behind the success of the algorithm is an adaptive projection onto the coarse grids that preserves the near null space. The resulting algorithm has weak dependence on the gauge coupling and exhibits very little critical slowing down in the chiral limit. Results are presented for the Wilson-Dirac operator of the 2D U(1) Schwinger model.

Strongly interacting dynamics and the search for new physics at the LHC

We present results for the spectrum of a strongly interacting SU(3) gauge theory with Nf = 8 light fermions in the fundamental representation. Carrying out nonperturbative lattice calculations at the lightest masses and largest volumes considered to date, we confirm the existence of a remarkably light singlet scalar particle. We explore the rich resonance spectrum of the 8-flavor theory in the context of the search for new physics beyond the standard model at the Large Hadron Collider (LHC). Lastly, connecting our results to models of dynamical electroweak symmetry breaking, we estimate the vector resonance mass to be about 2 TeV with a width of roughly 450 GeV, and predict additional resonances with masses below ~3 TeV.

Toward TeV Conformality

We study the chiral properties of an SU(3) gauge theory with N{f} massless Dirac fermions in the fundamental representation when N{f} is increased from 2 to 6. For N{f}=2, our lattice simulations lead to a value of psi psi/F{3}, where F is the Nambu-Goldstone-boson decay constant and psi psi is the chiral condensate, which agrees with the measured QCD value. For N{f}=6, this ratio shows significant enhancement, presaging an even larger enhancement anticipated as N{f} increases further, toward the critical value for transition from confinement to infrared conformality.

Parity Doubling and the S Parameter below the Conformal Window

We describe a lattice simulation of the masses and decay constants of the lowest-lying vector and axial resonances, and the electroweak S parameter, in an SU(3) gauge theory with N(f)=2 and 6 fermions in the fundamental representation. The spectrum becomes more parity doubled and the S parameter per electroweak doublet decreases when N(f) is increased from 2 to 6, motivating study of these trends as N(f) is increased further, toward the critical value for transition from confinement to infrared conformality.