Topics:

More recent papers are here.

Lattice gauge theory 

    Composite Higgs:   (see more recent work)

T. DeGrand, Y. Liu, E. T. Neil, Y. Shamir, and B. SvetitskySpectroscopy of SU(4) gauge theory with two flavors of sextet fermions, Phys. Rev. D 91 (2015) 114502 (arXiv:1501.05665 [hep-lat], January 2015).
We present a first look at the spectroscopy of SU(4) gauge theory coupled to two flavors of Dirac fermions in the two-index antisymmetric representation, which is a real representation.  We compute meson and diquark masses, the pseudoscalar and vector meson decay constants, and the masses of six-quark baryons. We make comparisons with large-Nc expectations.
TeX source and figures (zip file), PDF
T. DeGrand, Y. Liu, E. T. Neil, Y. Shamir, and B. SvetitskySpectroscopy of SU(4) lattice gauge theory with fermions in the two index anti-symmetric representation, talk given by Y. Liu at Lattice 2014, the 32nd International Symposium on Lattice Field Theory, Columbia University, New York, NY, June 2014, published in proceedings, Proc. Sci. (LATTICE2014) 275 (arXiv:1412.4851 [hep-lat], December 2014).
We present a study of spectroscopy of SU(4) lattice gauge theory coupled to two flavors of Dirac fermions in the anti-symmetric two index representation. The fermion representation is real, and the pattern of chiral symmetry breaking is SU(2Nf ) → SO(2Nf ) with N flavors of Dirac fermions. It is an interesting generalization of QCD, for several reasons: it allows direct exploration of an alternate large Nc expansion, it can be simulated at non-zero chemical potential with no sign problem, and several UV completions of composite Higgs systems are built on it. We present preliminary results on the baryon and meson spectra of the theory and compare them with SU(3) results and with expectations for large-Nc scaling.
TeX source and figures (zip file), PDF; slides of the talk
  Lattice technicolor:

T. DeGrand, Y. Shamir, and B. Svetitsky, Gauge theories with fermions in two-index representations, talk given by Y. Shamir at Lattice 2013, the 31st International Symposium on Lattice Field Theory, Mainz, Germany, July 2013, published in proceedings, Proc. Sci. (LATTICE 2013) 064 (arXiv:1310.2128 [hep-lat], October 2013).
After some introductory comments on the peculiar features of slowly running theories, I will report results obtained using the Schrödinger functional technique for two gauge theories that are believed to lie near the bottom of the conformal window: the SU(3) theory with two adjoint Dirac fermions, and the SU(4) theory with six Dirac fermions in the two-index antisymmetric representation. In both cases we find a small beta function in strong coupling, but we cannot confirm or rule out an infrared fixed point. In both theories the mass anomalous dimension levels off, staying well below 0.5, much like the theories with fermions in the two-index symmetric representation investigated earlier.
TeX source and figures (zip file), PDF; slides of the talk
T. DeGrand, Y. Shamir, and B. Svetitsky, Near the sill of the conformal window: gauge theories with fermions in two-index representations, Phys. Rev. D 88 (2013) 054505 (arXiv:1307.2425 [hep-lat], July 2013).
We apply Schrödinger functional methods to two gauge theories with fermions in two-index representations: the SU(3) theory with Nf = 2 adjoint fermions, and the SU(4) theory with Nf = 6 fermions in the two-index antisymmetric representation. Each theory is believed to lie near the bottom of the conformal window for its respective representation. In the SU(3) theory we find a small beta function in strong coupling but we cannot confirm or rule out an infrared fixed point. In the SU(4) theory we find a hint of walking—a beta function that approaches the axis and then turns away from it. In both theories the mass anomalous dimension remains small even at the strongest couplings, much like the theories with fermions in the two-index symmetric representation investigated earlier.
TeX source and figures (zip file), PDF
B. Svetitsky, Conformal or confining — results from lattice gauge theory for higher-representation gauge theories, talk given at Quark Confinement and the Hadron Spectrum X, Munich, Germany, October 2012, published in proceedings, Proc. Sci. ConfinementX (2013) 271 (arXiv:1301.1877 [hep-lat], January 2013).
We have calculated the running coupling in SU(2), SU(3), and SU(4) gauge theories to see whether they have infrared fixed points. An infrared fixed point means no confinement: It means that the long-distance physics is conformal, without a mass scale and indeed without a particle spectrum. All these theories run slowly in perturbation theory, so a nonperturbative determination is interesting.
TeX source and figures (zip file), PDF; slides of the talk

T. DeGrand, Y. Shamir, and B. Svetitsky, Mass anomalous dimension in sextet QCD, Phys. Rev. D 87 (2013) 074507 (arXiv:1201.0935 [hep-lat], Jan 2012, revised Dec 2012).
We extend our previous study of the SU(3) gauge theory with Nf = 2 flavors of fermions in the sextet representation of color. Our tool is the Schrödinger functional method. By changing the lattice action, we push the bulk transition of the lattice theory to stronger couplings and thus reveal the beta function and the mass anomalous dimension γm over a wider range of coupling, out to g2 ~ 11. Our results are consistent with an infrared fixed point, but walking is not ruled out. Our main result is that γm never exceeds 0.45, making the model unsuitable for walking technicolor.  We use a novel method of extrapolation to the large-volume/continuum limit, tailored to near-conformal theories.
TeX source and figures (zip file), PDF
T. DeGrand, Y. Shamir, and B. Svetitsky, SU(4) lattice gauge theory with decuplet fermions: Schrödinger functional analysisPhys. Rev. D 85 (2012) 074506 (arXiv:1202.2675 [hep-lat], February 2012).
We complete a program of study of SU(N) gauge theories coupled to two flavors of fermions in the two-index symmetric representation by performing numerical simulations in SU(4). The beta function, defined and calculated via the Schrödinger functional, runs more slowly than the two-loop perturbative result. The mass anomalous dimension levels off in strong coupling at a value of about 0.45, rendering this theory unsuitable for walking technicolor. A large-N comparison of these data with results from SU(2) and SU(3) reveals striking regularities.
TeX source and figures (zip file), PDF, Mathematica notebook (PDF)

T. DeGrand, Y. Shamir, and B. Svetitsky, Gauge theories with fermions in the two-index symmetric representation, talk given by T. DeGrand at Lattice 2011, the 29th International Symposium on Lattice Field Theory, Squaw Valley, Lake Tahoe, California, July 2011, published in proceedings, Proc. Sci. LATTICE2011 (2011) 060 (arXiv:1110.6845 [hep-lat], October 2011).
We summarize our recent work on gauge theories with two flavors of fermions in the two-index symmetric representation: SU(2) gauge theory with adjoint fermions, SU(3) with sextets, and SU(4) with ten-dimensional-representation fermions. All three systems have beta functions smaller than their perturbative value, approaching a fixed point near the expected two-loop zero. In all cases the mass anomalous dimension is small, under 0.5.
TeX source and figures (zip file), PDF; slides of the talk
T. DeGrand, Y. Shamir, and B. Svetitsky, Infrared fixed point in SU(2) gauge theory with adjoint fermions, Phys. Rev. D 83 (2011) 074507 (arXiv:1102.2843 [hep-lat], February 2011).
We apply Schrödinger-functional techniques to the SU(2) lattice gauge theory with Nf = 2 flavors of fermions in the adjoint representation. Our use of hypercubic smearing enables us to work at stronger couplings than did previous studies, before encountering a critical point and a bulk phase boundary.  Measurement of the running coupling constant gives evidence of an infrared fixed point g* where 1/g*2 = 0.20(4)(3). At the fixed point, we find a mass anomalous dimension γm(g*) = 0.31(6).
TeX source and figures (zip file), PDF
B. Svetitsky, Y. Shamir, and T. DeGrand Sextet QCD: slow running and the mass anomalous dimension, talk given at Lattice 2010, the 28th International Symposium on Lattice Field Theory, Villasimius, Italy, June 2010, published in proceedings, Proc. Sci. LATTICE2010 (2010) 072 (arXiv:1010.3396 [hep-lat], October 2010).
I report the results of Schrödinger functional calculations in the SU(3) gauge theory with two flavors of color sextet fermions, defined with the Wilson–clover action using nHYP fat links. While we cannot confirm the infrared fixed point seen with thin links, we find very slow evolution of the coupling constant, so slow that extraction of the mass anomalous dimension is straightforward.
TeX source and figures (zip file),  PDF; slides of the talk

T. DeGrand, Y. Shamir, and B. Svetitsky, Running coupling and mass anomalous dimension of SU(3) gauge theory with two flavors of symmetric-representation fermions, Phys. Rev. D 82 (2010) 054503 (arXiv:1006.0707 [hep-lat], June 2010).
We have measured the running coupling constant of SU(3) gauge theory coupled to Nf = 2 flavors of symmetric representation fermions, using the Schrödinger functional scheme. Our lattice action is defined with hypercubic smeared links which, along with the larger lattice sizes, bring us closer to the continuum limit than in our previous study. We observe that the coupling runs more slowly than predicted by asymptotic freedom, but we are unable to observe fixed point behavior before encountering a first order transition to a strong coupling phase. This indicates that the infrared fixed point found with the thin-link action is a lattice artifact. The slow running of the gauge coupling permits an accurate determination of the mass anomalous dimension for this theory, which we observe to be small, γm<  0.6, over the range of couplings we can reach.  We also study the bulk and finite-temperature phase transitions in the strong coupling region.
TeX source and figures (zip file), PDF

O. Machtey and B. SvetitskyMetastable nonconfining states in SU(3) lattice gauge theory with sextet fermionsPhys. Rev. D 81 (2010) 014501 (arXiv:0911.0886 [hep-lat], November 2009).
We study the SU(3) lattice gauge theory, with two flavors of sextet Wilson–clover fermions, near its finite-temperature phase transition. We find metastable states that have Wilson line expectation values whose complex phases are near 2π/3 or π. The true equilibrium phase at these couplings and temperatures has its Wilson line oriented only towards the positive real axis, in agreement with perturbation theory.
TeX source and figures (zip file), PDF

B. SvetitskyLattice gauge theory in technicolor, talk given at PANIC 2008, the International Conference on Particles and Nuclei, Eilat, Israel, November 2008, published in proceedings, Nucl. Phys. A 827 (2009) 547 (arXiv:0901.2103 [hep-lat], January 2009).
The methods of lattice gauge theory may be applied to gauge theories besides QCD, in fact to any gauge group and any representation of matter fields (as long as the coupling is not chiral).  Such theories are useful for model building beyond the Standard Model, for instance in technicolor models.  We have carried out Monte Carlo simulations of the SU(3) gauge theory with color sextet fermions.  Our result for its discrete beta function indicate an infrared fixed point that makes the theory conformal rather than confining.  The lattice theory's phase diagram shows no separation of chiral and confinement scales, measured when the quark mass is nonzero.
TeX source and figures (zip file), PDF; slides of the talk

T. DeGrand, Y. Shamir, and  B. Svetitsky, Phase structure of SU(3) gauge theory with two flavors of symmetric-representation fermions, Phys. Rev. D 79 (2009) 034501 (arXiv:0812.1427 [hep-lat], December 2008).
We have performed numerical simulations of SU(3) gauge theory coupled to Nf = 2 flavors of symmetric representation fermions. The fermions are discretized with the tadpole-improved clover action. Our simulations are  done on lattices of length L = 6, 8, and 12. In all simulation volumes we observe a crossover from  a strongly coupled confined phase to a weak coupling deconfined phase.  Degeneracies in screening masses, plus the behavior of the pseudoscalar decay constant, indicate that the deconfined phase is also a phase in which chiral symmetry is restored.  The movement of the confinement transition as the volume is changed is consistent with avoidance of the basin of attraction of an infrared fixed point of the massless theory.
TeX source and figures (zip file), PDF

T. DeGrandY. Shamir, and B. Svetitsky,   Exploring the phase diagram of sextet QCD, talk given at Lattice 2008, the 26th International Symposium on Lattice Field Theory, Williamsburg, Virginia, July 2008, published in proceedings, Proc. Sci. LATTICE2008 (2008) 063 (arXiv:0809.2953 [hep-lat], September 2008).
As a follow up to the previous talk about the beta function of SU(3) gauge theory with Nf=2 symmetric representation (clover) fermions [arXiv:0809.2885 [hep-lat]], we describe our explorations of the β–κ plane, away from the massless limit. Our simulations are mostly done on lattices of length L=8 and 12. We observe a phase transition from a strong coupling confined phase to a deconfined, chirally restored phase. The line of transitions avoids (so far) the location of the infrared fixed point discussed in the last talk.
TeX source and figures (gzipped tar file), PDF; slides of the talk

B. Svetitsky, Y. Shamir, and T. DeGrand Nonperturbative infrared fixed point in sextet QCD, talk given at Lattice 2008, the 26th International Symposium on Lattice Field Theory, Williamsburg, Virginia, July 2008, published in proceedings, Proc. Sci. LATTICE2008 (2008) 062 (arXiv:0809.2885 [hep-lat], September 2008).
The SU(3) gauge theory with fermions in the sextet representation is one of several theories of interest for technicolor models. We have carried out a Schrodinger functional (SF) calculation for the lattice theory with two flavors of Wilson fermions. We find that the discrete beta function changes sign when the SF renormalized coupling is in the neighborhood of g2=2.0, showing a breakdown of the perturbative picture even though the coupling is weak. The most straightforward interpretation is an infrared-stable fixed point.
TeX source and figures (zip file),  PDF; slides of the talk

Y. Shamir, B. Svetitsky, and T. DeGrand, Zero of the discrete beta function in SU(3) lattice gauge theory with color sextet fermions, Phys. Rev. D 78 (2008) 031502(R) (arXiv:0803.1707 [hep-th], March 2008).
We have carried out a Schrodinger functional (SF) calculation for the SU(3) lattice gauge theory with two flavors of Wilson fermions in the sextet representation of the gauge group. We find that the discrete beta function, which governs the change in the running coupling under a discrete change of spatial scale, changes sign when the SF renormalized coupling is in the neighborhood of g2=2.0. The simplest explanation is that the theory has an infrared-attractive fixed point, but more complicated possibilities are allowed by the data. While we compare rescalings by factors of 2 and 4/3, we work at a single lattice spacing.
TeX source and figures (zip file), PDF

    Dense baryonic matter:
B. Svetitsky, Y. Shamir, and M. Golterman, Why (staggered fermions)1/4 fail at finite density, talk given at Lattice 2006, the 24th International Symposium on Lattice Field Theory, Tucson, Arizona, July 2006, published in proceedings, Proc. Sci. LAT2006 (2006) 148 (hep-lat/0609051, September 2006).
Because the staggered fermion determinant is complex at nonzero μ, taking its fourth root leads to phase ambiguities. These unphysical effects cause the measure to become discontinuous; the problem becomes acute when Re μ ≳ mπ/2 (when T > 0 this rough bound probably moves towards larger μ). We show how to overcome the problem, but only very close to the continuum limit. This regime may be beyond reach with current resources.
TeX source and figures (zip file), PDF

M. Golterman, Y. Shamir, and B. Svetitsky, Breakdown of staggered fermions at nonzero chemical potential, Phys. Rev. D 74 (2006) 071501(R) (hep-lat/0602026, February 2006 [revised August 2006]).
The staggered fermion determinant is complex when the quark chemical potential μ is nonzero. Its fourth root, used in simulations with dynamical fermions, will have phase ambiguities that become acute when Re μ is sufficiently large. We show how to resolve these ambiguities, but our prescription only works very close to the continuum limit. We argue that this regime is very far from current capabilities. Other procedures require being even closer to the continuum limit, or fail altogether, because of unphysical discontinuities in the measure. At zero temperature the breakdown is expected when Re μ ≳ mπ/2. Estimates of the location of the breakdown at nonzero temperature are less certain.
TeX source and figures (zip file), PDF

B. Bringoltz, Dense baryonic matter in strong coupling lattice gauge theory (Ph.D. thesis), hep-lat/0407018, July 2004.
B. Bringoltz and B. Svetitsky, Anisotropic Goldstone bosons of strong-coupling lattice QCD at high density, Phys. Rev. D 70 (2004) 074512 (hep-lat/0405013, May 2004).
We calculate the spectrum of excitations in strong-coupling lattice QCD in a background of fixed baryon density, at a substantial fraction of the saturation density. We employ a next-nearest-neighbor fermion formulation that possesses the SU(Nf)xSU(Nf) chiral symmetry of the continuum theory. We find two types of massless excitations: type I Goldstone bosons with linear dispersion relations and type II Goldstone bosons with quadratic dispersion relations.  Some of the type I bosons originate as type II bosons of the nearest-neighbor theory. Bosons of either type can develop anisotropic dispersion relations, depending on the value of Nf and the baryon density.
TeX sourcePDF

B. Bringoltz and B. Svetitsky, Spontaneous symmetry breaking in strong-coupling lattice QCD at high density , Phys. Rev. D 69 (2004) 014502 (hep-lat/0310032, October 2003).
We determine the patterns of spontaneous symmetry breaking in strong-coupling lattice QCD in a fixed background baryon density. We employ a next-nearest-neighbor fermion formulation that possesses the SU(Nf)xSU(Nf) chiral symmetry of the continuum theory. We find that the global symmetry of the ground state varies with Nf and with the background baryon density. In all cases the condensate breaks the discrete rotational symmetry of the lattice as well as part of the chiral symmetry group.
TeX sourcePDF

B. Bringoltz, Order from disorder in lattice QCD at high density, Phys. Rev. D 69 (2004) 014508 (hep-lat/0308018, August 2003).
B. Bringoltz and B. Svetitsky, Lattice gauge theory with baryons at strong coupling, Phys. Rev. D 68 (2003) 034501 (hep-lat/0211018, November 2002).
We study the effective Hamiltonian for strong-coupling lattice QCD in the case of non-zero baryon density. In leading order the effective Hamiltonian is a generalized antiferromagnet. For naive fermions, the symmetry is U(4Nf) and the spins belong to a representation that depends on the local baryon number. Next-nearest neighbor (nnn) terms in the Hamiltonian break the symmetry to U(Nf)xU(Nf). We transform the quantum problem to a Euclidean sigma model which we analyze in a 1/Nc expansion. In the vacuum sector we recover spontaneous breaking of chiral symmetry for the nearest-neighbor and nnn theories. For non-zero baryon density we study the nearest-neighbor theory only, and show that the pattern of spontaneous symmetry breaking depends on the baryon density.
TeX source and figures (packed with uufiles), PDF

B. Bringoltz and B. Svetitsky, Towards a strong-coupling theory of QCD at finite density, talk and poster presented at Lattice 2002, the 20th International Symposium on Lattice Field Theory, Cambridge, Massachusetts, June 2002, published in Nucl. Phys. B (Proc. Suppl.) 119 (2003) 565 (hep-lat/0209005, September 2002).
We apply strong-coupling perturbation theory to the QCD lattice Hamiltonian. We begin with naive, nearest-neighbor fermions and subsequently break the doubling symmetry with next-nearest-neighbor terms. The effective Hamiltonian is that of an antiferromagnet with an added kinetic term for baryonic "impurities," reminiscent of the t–J model of high-Tc superconductivity. As a first step, we fix the locations of the baryons and make them static. Following analyses of the t–J model, we apply large-N methods to obtain a phase diagram in the (Nc,Nf) plane at zero temperature and baryon density. Next we study a simplified U(3) toy model, in which we add baryons to the vacuum. We use a coherent state formalism to write a path integral which we analyze with mean field theory, obtaining a phase diagram in the (nB,T) plane.
TeX source and figures (packed with uufiles), PDF

    Microscopic spectra:
K. Splittorff and B. Svetitsky, The sign problem via imaginary chemical potential, Phys. Rev. D 75 (2007) 114504 (hep-lat/0703004, March 2007).
We calculate an analogue of the average phase factor of the staggered fermion determinant at imaginary chemical potential. Our results from the lattice agree well with the analytical predictions in the microscopic regime for both quenched and phase-quenched QCD. We demonstrate that the average phase factor in the microscopic domain is dominated by the lowest-lying Dirac eigenvalues.
TeX source and figures (zip file), PDF

P. H. Damgaard, U. M. Heller, K. SplittorffB. Svetitsky, and D. Toublan, Microscopic eigenvalue correlations in QCD with imaginary isospin chemical potential, Phys. Rev. D 73 (2006) 105016 (hep-th/0604054, April 2006).
We consider the chiral limit of QCD subjected to an imaginary isospin chemical potential. In the ε-regime of the theory we can perform precise analytical calculations based on the zero-momentum Goldstone modes in the low-energy effective theory. We present results for the spectral correlation functions of the associated Dirac operators.
TeX source and figures (zip file), PDF

P. H. Damgaard, U. M. Heller, K. SplittorffB. Svetitsky, and D. Toublan, Extracting Fπ from small lattices: unquenched results, Phys. Rev. D 73 (2006) 074023 (hep-lat/0602030, February 2006).
We calculate the response of the microscopic Dirac spectrum to an imaginary isospin chemical potential for QCD with two dynamical flavors in the chiral limit. This extends our previous calculation from the quenched to the unquenched theory. The resulting spectral correlation function in the ε-regime provides here, too, a new and efficient way to measure Fπ on the lattice. We test the method in a hybrid Monte Carlo simulation of the theory with two staggered quarks.
TeX source and figures (zip file), PDF

P. H. Damgaard, U. M. Heller, K. Splittorff, and B. Svetitsky, A new method for determining Fπ on the lattice, Phys. Rev. D 72 (2005) 091501(R) (hep-lat/0508029, August 2005)
We derive the two-point spectral correlation function of the Dirac operator with a specific external source in the ε-regime of QCD. This correlation function has a unique and strong dependence on Fπ, and thus provides an novel way to extract Fπ from lattice simulations. We test the method in a quenched lattice simulation with staggered fermions.
TeX source and figures (packed with uufiles), PDF

P. H. Damgaard, U. M. Heller, R. Narayanan, and B. Svetitsky, Divergent chiral condensate in the quenched Schwinger model, Phys. Rev. D 71 (2005) 114503 (hep-lat/0504012, April 2005).
We calculate numerically the eigenvalue distribution of the overlap Dirac operator in the quenched Schwinger model on a lattice. The distribution does not fit any of the three universality classes of spontaneous chiral symmetry breaking, and its strong volume dependence indicates that the chiral condensate in the quenched theory is an ill-defined and divergent quantity. When we reweight configurations with the Dirac determinant to study the theory with Nf = 1, we obtain a distribution of eigenvalues that is well-behaved and consistent with the theory of explicit symmetry breaking due to the anomaly.
TeX source and figures (packed with uufiles), PDF

P. H. Damgaard, U. M. Heller, R. Niclasen, and B. Svetitsky, Patterns of spontaneous chiral symmetry breaking in vectorlike gauge theories, Nucl. Phys. B 633 (2002) 97–113 (hep-lat/0110028, October 2001).
It has been conjectured that spontaneous chiral symmetry breaking in strongly coupled vectorlike gauge theories falls into only three different depending on the gauge group and the representations carried by the fermions. We test this proposal by studying SU(2), SU(3) and SU(4) lattice gauge theories with staggered fermions in different irreducible representations. Staggered fermions away from the continuum limit should, for all complex representations, still belong to the continuum class of spontaneous symmetry breaking. But for all real and pseudo-real representations we show that staggered fermions should belong to incorrect symmetry breaking classes away from the continuum, thus generalizing previous results. As an unambiguous signal for whether chiral symmetry breaks, and which breaking pattern it follows, we look at the smallest Dirac eigenvalue distributions. We find that the patterns of symmetry breaking are precisely those conjectured.
TeX source and figures (packed with uufiles), PDF

    Lattice fermions and chiral symmetry:
B. Svetitsky, Localization and lattice fermions, talk given at the Workshop on Computational Hadron Physics, Nicosia, Cyprus, September 14–17, 2005, published in Nucl. Phys. B (Proc. Suppl.) 153 (2006) 314 (hep-lat/0511003, November 2005)
I review how the phenomenology of localization applies to fermions in lattice gauge theory, present measurements of the localization length and other quantities, and discuss the consequences for things like the overlap kernel.
TeX source and figures (zip file), PDF

B. Svetitsky, Y. Shamir, and M. Golterman, Localization of lattice fermions: lessons for overlap, talk given at Lattice 2005, the 23rd International Symposium on Lattice Field Theory, Dublin, July 2005, published in proceedings, Proc. Sci. LAT2005 (2005) 129 (hep-lat/0508015, August 2005).
Lattice fermions in a fluctuating gauge field can show localization, much like electrons in a disordered potential. We study the spectrum of localized and extended states of supercritical Wilson fermions in gauge ensembles generated with plaquette and improved actions. When the Wilson fermion operator is used to construct the overlap kernel, the mobility edge, that is the boundary between the localized and extended states, determines the range of the kernel.
TeX source and figures (packed with uufiles), PDF

M. Golterman, Y. Shamir, and B. Svetitsky, Localization properties of lattice fermions with plaquette and improved gauge actions, Phys. Rev. D 72 (2005) 034501 (hep-lat/0503037, March 2005).
We determine the location λc of the mobility edge in the spectrum of the hermitian Wilson operator in pure-gauge ensembles with plaquette, Iwasaki, and DBW2 gauge actions. The results allow mapping a portion of the (quenched) Aoki phase diagram. We use Green function techniques to study the localized and extended modes. Where λc>0 we characterize the localized modes in terms of an average support length and an average localization length, the latter determined from the asymptotic decay rate of the mode density. We argue that, since the overlap operator is commonly constructed from the Wilson operator, its range is set by the value of λc−1 for the Wilson operator. It follows from our numerical results that overlap simulations carried out with a cutoff of 1 GeV, even with improved gauge actions, could be afflicted by unphysical degrees of freedom as light as 250 MeV.
TeX source and figures (packed with uufiles), PDF

M. Golterman, Y. Shamir, and B. Svetitsky, Mobility edge in lattice QCD, Phys. Rev. D 71 (2005) 071502(R) (hep-lat/0407021, July 2004).
We determine the location λc of the mobility edge in the spectrum of the hermitian Wilson operator on quenched ensembles. We confirm a theoretical picture of localization proposed for the Aoki phase diagram. When λc>0 we also determine some key properties of the localized eigenmodes with eigenvalues |λ|<λc. Our results lead to simple tests for the validity of simulations with overlap and domain-wall fermions.
TeX source and figure (packed with uufiles), PDF

Federico Berruto, Richard C. Brower, and Benjamin Svetitsky, Effective Lagrangian for strongly coupled domain wall fermions, Phys. Rev. D 64 (2001) 114504 (hep-lat/0105016, May 2001).
We derive the effective Lagrangian for mesons in lattice gauge theory with domain-wall fermions in the strong-coupling and large-Nc limits. We use the formalism of supergroups to deal with the Pauli–Villars fields, needed to regulate the contributions of the heavy fermions. We calculate the spectrum of pseudo-Goldstone bosons and show that domain wall fermions are doubled and massive in this regime. Since we take the extent and lattice spacing of the fifth dimension to infinity and zero respectively, our conclusions apply also to overlap fermions.
TeX source and figures (packed with uufiles), compressed Postscript , PDF
Richard C. Brower and Benjamin Svetitsky, Hamiltonian domain wall fermions at strong coupling, Phys. Rev. D 61 (2000) 114511 (hep-lat/9912019, December 1999). TeX sourcePDF
    Other lattice papers:

B. Svetitsky, O. Raviv, and Y. Shamir, Beta function of three-dimensional QED, talk given at Lattice 2014, the 32nd International Symposium on Lattice Field Theory, Columbia University, New York, NY, June 2014, published in proceedings, Proc. Sci. (LATTICE2014) 051 (arXiv:1410.0118 [hep-lat], October 2014).
We have carried out a Schrödinger-functional calculation for the Abelian gauge theory with Nf = 2 four-component fermions in three dimensions. We find no fixed point in the beta function, meaning that the theory is confining rather than conformal.
TeX source and figures (zip file), PDF; slides of the talk
T. DeGrand, Y. Shamir, and B. Svetitsky, Suppressing dislocations in normalized hypercubic smearing, Phys. Rev. D 90 (2014) 054501 (arXiv:1407.4201 [hep-lat], July 2014).
Normalized hypercubic smearing improves the behavior of dynamical Wilson-clover fermions, but has the unwanted side effect that it can occasionally produce spikes in the fermion force. These spikes originate in the chain rule connecting the derivative with respect to the smeared links to the derivative with respect to the dynamical links, and are associated with the presence of dislocations in the dynamical gauge field. We propose and study an action designed to suppress these dislocations. We present evidence for improved performance of the hybrid Monte Carlo algorithm. A side benefit is improvement in the properties of valence chiral fermions.
TeX source and figures (zip file), PDF
O. Raviv, Y. Shamir, and B. Svetitsky, Non-perturbative beta function in three-dimensional electrodynamics, Phys. Rev. D 90 (2014) 014512 (arXiv:1405.6916 [hep-lat], May 2014).
We apply the Schrödinger functional method to the Abelian gauge theory in three dimensions with Nf = 2 four-component fermions. We find that the calculated beta function does not cross zero in the range of coupling we study. This implies that the theory exhibits confinement and mass generation, rather than a conformal infared regime.
TeX source and figures (zip file), PDF
Y. Shamir, B. Svetitsky, and E. Yurkovsky,  Improvement via hypercubic smearing in triplet and sextet QCD, Phys. Rev. D 83 (2011) 097502 (arXiv:1012.2819 [hep-lat], December 2010).
We study non-perturbative improvement in SU(3) lattice gauge theory coupled to fermions in the fundamental and two-index symmetric representations. Our lattice action is defined with hypercubic smeared links incorporated into the Wilson–clover fermion kernel. Using standard Schrödinger-functional techniques we estimate the clover coefficient cSW and find that discretization errors are much smaller than in thin-link theories.
TeX source and figures (zip file),  PDF

A. Kovner and B. Svetitsky, Interaction potential in compact three-dimensional QED with mixed action, Phys. Rev D 60 (1999) 105032 (hep-lat/9811015, November 1998).
We use a variational wave function to calculate the energy of the interaction between external charges in the compact Abelian gauge theory in 2+1 dimensions with mixed action. Our variational wave functions preserve the compact gauge invariance of the theory both in the vacuum and in the charged sectors. We find that a good estimate of the interaction energy is obtained only when we allow more variational parameters in the charged sector than in the vacuum sector. These extra parameters are the profile of an induced electric field. We find that the theory has a two-phase structure: When the charge-2 coupling is large and negative there is no mass gap in the theory and no confinement, while otherwise a mass gap is generated dynamically and the theory confines charges. The pure Wilson theory is in the confining phase.
TeX sourcePDF
B. Svetitsky, Effective spin models for the confinement phase transition, lecture given at Lattice '98, the 16th International Symposium on Lattice Field Theory, Boulder, Colorado, 13–18 July 1998, published in Nucl. Phys. B (Proc. Suppl) 73 (1999) 432 (hep-lat/9809008, September 1998).
Spatial correlations—bubbles, domain walls, etc.—can best be studied by concentrating on the degrees of freedom most relevant to the problem. For the finite temperature confinement transition, I integrate out all gauge degrees of freedom, leaving only spins—Ising or Potts—related to the Wilson line. I present problems that arise in the course of this transformation and some results for the effective spin action.
TeX sourcePDF
B. Svetitsky, The Ising model and bubbles in the quark–gluon plasma, lecture given at Strong and Electroweak Matter, Eger, Hungary, 21–25 May 1997, published in proceedings, ed. F. Csikor and Z. Fodor (World Scientific, Singapore, 1998), pp. 342–6 (hep-lat/9707031, July 1997).
I review evidence for the stability of bubbles in the quark–gluon plasma near the confinement phase transition. In analogy with the much-studied oil–water emulsions, this raises the possibility that there are many phases between the pure plasma and the pure hadron gas, characterized by spontaneous inhomogeneity and modulation. In studying emulsions, statistical physicists have reproduced many of their phases with microscopic models based on Ising-like theories with competing interactions. Hence we seek an effective Ising Hamiltonian for the SU(3) gauge theory near its transition.
TeX source and figures (packed with uufiles),  PDF
B. Svetitsky and N. Weiss, Ising description of the transition region in SU(3) gauge theory at finite temperature, Phys. Rev. D 56 (1997) 5395 (hep-lat/9705007).
We attempt the numerical construction of an effective action in three dimensions for Ising spins which represent the Wilson lines in the four-dimensional SU(3) gauge theory at finite temperature. For each configuration of the gauge theory, each spin is determined by averaging the Wilson lines over a small neighborhood and then projecting the average to ±1 according to whether the neighborhood is ordered or disordered. The effective Ising action, determined via the lattice Schwinger–Dyson equations, contains even (two-spin) and odd (one- and three-spin) terms with short range. We find that the truncation to Ising degrees of freedom produces an effective action which is discontinuous across the gauge theory's phase transition. This discontinuity may disappear if the effective action is made more elaborate.
TeX source and figures (packed with uufiles), PDF
G. G. Batrouni and B. Svetitsky, Accelerated dynamics in simulations of first order phase transitions, Phys. Rev. B 36 (1987) 5647.
We show that the method of Fourier acceleration can be used to speed up tunnelling between degenerate and near-degenerate minima of the effective action in Langevin simulations of statistical systems. Acceleration factors of two orders of magnitude are attained. Such improvement is important for the study of weak first-order phase transitions, for simulations of spin glasses, and for the method of simulated annealing in optimization problems.
TeX source (packed with uufiles), compressed Postscript — figures not available 
B. Svetitsky, Status of lattice gauge theory, Invited talk at Quark Matter '86, the Fifth International Conference on Ultra-Relativistic Nucleus–Nucleus Collisions, Asilomar, CA, April 1986. Published in proceedings, Nucl. Phys. A 461 (1987) 71c.
I review recent results in lattice gauge theory which pertain to the physics of high temperature QCD and quark matter. Topics covered include: the transition temperature in the pure gauge theory, scaling and finite temperature behavior in QCD with quarks, problems with non-zero chemical potential, the Debye length in the high temperature plasma, and quantities relevant to thermodynamics and hydrodynamics in relativistic heavy ion collisions.
TeX source (gzipped), compressed Postscript — figures not available 

Lattice Skyrme model


Alec J. Schramm and Benjamin Svetitsky, Exact topological density in the lattice Skyrme model, lecture given at Strong and Electroweak Matter (SEWM 2000), Marseille, France, 14–17 June 2000, published in proceedings, ed. C. P. Korthals-Altes (World Scientific, Singapore, 2001), pp. 359–63 (hep-lat/0009014, September 2000). TeX source and figures (packed with uufiles), PDF
Alec J. Schramm and Benjamin Svetitsky, Topology and metastability in the lattice Skyrme model, Phys. Rev. D 62 (2000) 114020 (hep-lat/0008003, August 2000). TeX source and figures (packed with uufiles), PDF

Hot nuclear matter


G. Kälbermann, J. M. Eisenberg, and B. Svetitsky, Hot nuclear matter with dilatons, Nucl. Phys. A600 (1996) 436 (nucl-th/9510011).
We study hot nuclear matter in a model based on nucleon interactions deriving from the exchange of scalar and vector mesons. The main new feature of our work is the treatment of the scale breaking of quantum chromodynamics through the introduction of a dilaton field. Although the dilaton effects are quite small quantitatively, they affect the high-temperature phase transition appreciably. We find that inclusion of the dilaton leads to a metastable high-density state at zero pressure, similar to that found by Glendenning who considered instead the admixture of higher baryon resonances.
TeX source and figures (packed with uufiles), PDF

Particle pair production in strong electric fields; the flux tube model


Benjamin Svetitsky, The flux-tube model of particle creation in nuclear collisions, invited lecture at Nuclear Matter, Hot and Cold: a Symposium in Memory of Judah M. Eisenberg, Tel Aviv, Israel, April 14–16, 1999, published in proceedings, ed. by J. Alster and D. Ashery (Tel Aviv University, Tel Aviv, 2000), pp. 33–40 (hep-ph/9907278, July 1999). TeX source and figures (packed with uufiles), PDF
Melissa A. Lampert and Benjamin Svetitsky, Flux tube dynamics in the dual superconductor, talk given at Quark Matter '99, the 14th International Conference on Ultra-Relativistic Nucleus–Nucleus Collisions, Torino, Italy, May 10–14, 1999, published in Nucl. Phys. A 661 (1999) 617c (hep-ph/9907246, July 1999). TeX source and figures (packed with uufiles), PDF
Melissa A. Lampert and Benjamin Svetitsky, Flux tube dynamics in the dual superconductor, Phys. Rev. D 61 (2000) 034011 (hep-ph/9905455, May 1999). TeX source and figures (gzipped tar file), PDF
F. Cooper, J. M. Eisenberg, Y. Kluger, E. Mottola, and B. Svetitsky, Particle production in the central rapidity region, Phys. Rev. D 48 (1993) 190 (hep-ph/9212206).
We study pair production from a strong electric field in boost-invariant coordinates as a simple model for the central rapidity region of a heavy-ion collision. We derive and solve the renormalized equations for the time evolution of the mean electric field and current of the produced particles, when the field is taken to be a function only of the fluid proper time. We find that a relativistic transport theory with a Schwinger source term modified to take Pauli blocking (or Bose enhancement) into account gives a good description of the numerical solution to the field equations. We also compute the renormalized energy–momentum tensor of the produced particles and compare the effective pressure, energy and entropy density to that expected from hydrodynamic models of energy and momentum flow of the plasma.
Text: TeX source (gzipped), PDF
Figures: uufiles, compressed Postscript 
Y. Kluger, J. M. Eisenberg, and B. Svetitsky, Pair production in a strong electric field: an initial value problem in quantum field theory, Int. J. Mod. Phys. E 2 (1993) 333 (hep-ph/0311293).
We review recent achievements in the solution of the initial-value problem for quantum back-reaction in scalar and spinor QED. The problem is formulated and solved in the semiclassical mean-field approximation for a homogeneous, time-dependent electric field. Our primary motivation in examining back-reaction has to do with applications to theoretical models of production of the quark–gluon plasma, though we here address practicable solutions for back-reaction in general. We review the application of the method of adiabatic regularization to the Klein–Gordon and Dirac fields in order to renormalize the expectation value of the current and derive a finite coupled set of ordinary differential equations for the time evolution of the system. Three time scales are involved in the problem and therefore caution is needed to achieve numerical stability for this system. Several physical features, like plasma oscillations and plateaus in the current, appear in the solution. From the plateau of the electric current one can estimate the number of pairs before the onset of plasma oscillations, while the plasma oscillations themselves yield the number of particles from the plasma frequency.

We compare the field-theory solution to a simple model based on a relativistic Boltzmann–Vlasov equation, with a particle production source term inferred from the Schwinger particle creation rate and a Pauli-blocking (or Bose-enhancement) factor. This model reproduces very well the time behavior of the electric field and the creation rate of charged pairs of the semiclassical calculation. It therefore provides a simple intuitive understanding of the nature of the solution since nearly all the physical features can be expressed in terms of the classical distribution function.

TeX source (packed with uufiles), PDF
Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, and E. Mottola, Fermion pair production in a strong electric field, Phys. Rev. D 45 (1992) 4659.
The initial-value problem for the quantum back-reaction in spinor QED is formulated and solved in the semiclassical mean field approximation, for a homogeneous but time-dependent electric field E(t). We apply the method of adiabatic regularization to the Dirac equation in order to renormalize the expectation value of the current and derive a finite coupled set of ordinary differential equations for the time evolution of the system. We solve this system in (1+1) dimensions numerically and compare the solution to a simple model based on a relativistic Boltzmann–Vlasov equation, with a particle production source term inferred from the Schwinger particle creation rate and a Pauli-blocking factor. This model reproduces very well the time behavior of the electric field and the creation rate of electron–positron pairs of the semiclassical calculation.
TeX source (gzipped), compressed Postscript — figures not available 
J. M. Eisenberg, Y. Kluger, and B. Svetitsky, Pair production in a strong electric field with back-reaction—an interim summary, Acta Phys. Polon. B23 (1992) 577.

Dedicated to Wiesław Czyż on the occasion of his sixty-fifth birthday

We present a summary of the present status of efforts to solve the problem in which pairs are produced in a strong electric field, are accelerated by it, and then react back on it through the counter-field produced by their current. This picture has been used by Bialas and Czyz and others as a model for effects that may possibly arise in the study of the quark–gluon plasma. We here give a didactic review of recent developments in this back-reaction problem. We first present a simple version of the theory of pair tunneling from a fixed electric field, and then sketch how this has been applied to the quark–gluon plasma. Then we turn to a field formulation of the problem for charged bosons, which leads to the need to carry out a renormalization program, outlined again in simple terms. Numerical results for this program are presented for one spatial dimension, the corresponding physical behavior of the system is discussed, and the implications for three spatial dimensions are considered. We exhibit a phenomenological transport equation embodying physics that is essentially identical to that of the field formulation, thus helping to tie the model of Bialas and Czyz for the quark–gluon plasma to a field-theory formulation. Last, we note the status of extensions to (i) the problem with three space dimensions; (ii) the fermion case; (iii) the formulation in terms of boost-invariant variables (as desirable for the quark–gluon plasma); and (iv) transport equations derived in a fundamental and consistent fashion.
TeX source and figures (packed with uufiles)
Text: PDF
Figures: PDF
Y. Kluger, J. M. Eisenberg, B. Svetitsky, F. Cooper, and E. Mottola, Pair production in a strong electric field, Phys. Rev. Lett. 67 (1991) 2427.
We investigate the mechanism of pair creation in scalar QED from spatially homogeneous strong electric fields in 1+1 dimensions. Solution of the semiclassical field equations shows particle creation followed by plasma oscillations. We compare our results with a model based on a relativistic Boltzmann–Vlasov equation with a pair-creation source term related to the Schwinger mechanism. The time evolution of the electric field and the current obtained from the Boltzmann–Vlasov model is surprisingly similar to that found in the semiclassical calculation.
TeX source (gzipped),compressed Postscript — figures not available 

Charm in the quark–gluon plasma


B. Svetitsky and A. Uziel, Charm—a thermometer of the mixed phase, talk delivered by B.S. at RHIC '97, The Brookhaven Theory Workshop on Relativistic Heavy Ions, July 6–16, 1997 (hep-ph/9709228, September 1997).
A charmed quark experiences drag and diffusion in the quark–gluon plasma, as well as strong interaction with the plasma surface. Our simulations indicate that charmed quarks created in heavy ion collisions will be trapped in the mixed phase and will come to equilibrium in it. Their momentum distribution will thus reflect the temperature at the confinement phase transition.
TeX source and figures (packed with uufiles), PDF

Also given (under the title Charm transverse momentum as a thermometer of the quark–gluon plasma) at the International Europhysics Conference on High Energy Physics, 19–26 August 1997, Jerusalem, Israel, published in proceedings, ed. by D. Lellouch, G. Mikenberg, and E. Rabinovici (Springer-Verlag, Berlin, 1999).

See PDF of this version. 


B. Svetitsky and A. Uziel, Passage of charmed particles through the mixed phase in high-energy heavy-ion collisions, Phys. Rev. D 55 (1997) 2616 (hep-ph/9606284).
We employ a modified cascade hydrodynamics code to simulate the phase transition of an expanding quark–gluon plasma and the passage of a charmed particle through it. When inside the plasma droplets, the charmed quark experiences drag and diffusion forces. When outside the plasma, the quark travels as a D meson and experiences collisions with pions. Additional energy transfer takes place when the quark enters or leaves a droplet. We find that the transverse momentum of D mesons provides a rough thermometer of the phase transition.
TeX source and figures (packed with uufiles),  PDF
D. Levin-Plotnik and B. Svetitsky, Charmonium in a weakly coupled quark–gluon plasma, Phys. Rev. D 52 (1995) 4248 (hep-ph/9503305).
We present a model of charmonium as two heavy quarks propagating classically in a weakly coupled quark–gluon plasma. The quarks interact via a static, color-dependent potential and also suffer collisions with the plasma particles. We calculate the radiation width of the color octet state (for fixed, classical q–q separation) and find that it is long-lived provided a finite gluon mass is used to provide a threshold energy.
TeX source and figures (packed with uufiles), compressed Postscript, PDF (text, figures)

B. Svetitsky, Multiple hadronization and the momentum of charmed paricles: A signal of the quark–gluon plasma, Phys. Lett. 227B (1989) 450.
The transverse momentum of charmed particles should be degraded by their passage through the mixed-phase fluid existing during the course of the quark–hadron phase transition. The mechanism is multiple hadronization, the repeated conversion of quark to meson and back as phase boundaries are crossed. The effect can be reduced by thermalization via ordinary collisions with particles of the plasma or the hadron gas, but passage through a homogeneous dense medium cannot duplicate the momentum degradation. The effect, if observed, thus constitutes a simple signal of the confinement phase transition in ultrarelativistic heavy ion collisions.
TeX source (gzipped), compressed Postscript — figures not available 
B. Svetitsky, Charm—a signal of the confinement phase transition, lecture given at the NATO Advanced Study Institute on the Nuclear Equation of State, Peñiscola, Spain, May 21 – June 3, 1989. Published in proceedings, The Nuclear Equation of State, eds. W. Greiner and H. Stocker (Plenum, New York, 1989), pp. 257–267.
I review some features of proposed signals of the quark–gluon plasma in order to motivate a new one: measuring the transverse momentum of charmed particles. A sharp reduction in the mean transverse momentum will show that the matter produced in a collision has cooled through a confinement phase transition. I argue that other mechanisms will not duplicate the effect.
TeX source (packed with uufiles),compressed Postscript — figures not available 

Perturbative field theory at high temperature


G. Baym, J.-P. Blaizot, and B. Svetitsky, Emergence of new quasiparticles in quantum electrodynamics at finite temperature, Phys. Rev. D 46 (1992) 4043.
We study the spectrum of fermionic excitations in a hot relativistic electron plasma. Numerical analysis of the one-loop electron propagator shows the appearance of new quasiparticle modes, as already studied in the massless fermion limit, as the temperature is raised to exceed the electron mass m. We calculate the relevant spectral densities and show that one mode, whose splitting in energy from the original electron is of the order of m, moves as the temperature is raised towards higher frequency along with the electron while retaining a narrow width. We compare results derived in Feynman and Coulomb gauge. The role of the zero temperature counterterms is discussed.
TeX source and figures(packed with uufiles), compressed Postscript

Bubbles in the quark–gluon plasma


G. Lana and B. Svetitsky, Instability of bubbles in the quark–gluon plasma, Phys. Lett. B285 (1992) 251.
The small surface tension of vacuum bubbles in the quark–gluon plasma, along with a negative curvature tension, implies that the uniform plasma is unstable against the formation of such bubbles. We show that spherical bubbles are in turn unstable against the growth of irregularities.
TeX source (gzipped), PDF
B. Svetitsky, Swiss cheese instability of the quark–gluon plasma, presented at Workshop on QCD at Finite Temperature and Density, Upton, N.Y., Aug 1991. Published in Hot Summer Daze, ed. by A. Gocksch and R. Pisarski (World Scientific, Singapore, 1992) p. 204.
I argue that the quark–gluon plasma is unstable against the formation of bubbles of confining vacuum.
TeX source (gzipped), compressed Postscript — figures not available 
I. Mardor and B. Svetitsky, Bubble free energy at the quark–hadron phase transition, Phys. Rev. D 44 (1991) 878.
We calculate the free energy of finite droplets of quark–gluon plasma, and of finite hadronic bubbles in the bulk plasma, near the confinement phase transition. We sum over free quark and gluon energy levels in the presence of MIT bag boundary conditions. We find that the curvature term in the free energy, proportional to the radius of the droplet/bubble, is far more important than the contribution of the surface tension, proportional to the radius squared. This affects the critical radius for nucleation of plasma droplets in the superheated hadron gas, and seems to lead to instability of the plasma (even when not supercooled) against nucleation of hadron gas bubbles.
TeX source and figures (packed with uufiles), compressed Postscript

Judah


D. S. Koltun and B. Svetitsky, Judah Moshe Eisenberg (obituary), published in Physics Today, October 1998.

HTML, Postscript, PDF