BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Simon Hands appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We present the most recent results from the FASTSUM collaboration for hadron properties at high temperature. This includes the temperature dependence of the light and charmed meson and baryon spectrum, as well as properties of heavy quarkonia. The results are obtained using anisotropic lattices with a fixed scale approach. We also present the status of our next generation gauge ensembles.
We employ the domain wall fermion (DWF) formulation of the Thirring model on a lattice in 2+1+1 dimensions and perform $N=1$ flavor Monte Carlo simulations. At a critical interaction strength the model features a spontaneous $\mathrm{U}(2)\rightarrow\mathrm{U}(1)\otimes \mathrm{U}(1)$ symmetry breaking; we analyse the induced spin-0 mesons, both Goldstone and non-Goldstone, as well as the correlator of the fermion quasiparticles, in both resulting phases. Crucially, we determine the anomalous dimension $η_ψ\approx 3$ at the critical point, in stark contrast with the Gross-Neveu model in 3$d$ and with results obtained with staggered fermions. Our numerical simulations are complemented by an analytical treatment of the free fermion correlator, which exhibits large early-time artifacts due to branch cuts in the propagator stemming from unbound interactions of the fermion with its heavy doublers. These artifacts are generalisable beyond the Thirring model, being an intrinsic property of DWF, or more generally Ginsparg-Wilson fermions.
The infamous sign problem makes it impossible to probe dense (baryon density $μ_B>0$) QCD at temperatures near or below the deconfinement threshold. As a workaround, one can explore QCD-like theories such as two-colour QCD (QC2D) which don't suffer from this sign problem but are qualitively similar to real QCD. Previous studies on smaller lattice volumes have investigated deconfinement and colour superfluid to normal matter transitions. In this study we look at a larger lattice volume $N_s=24$ in an attempt to disentangle finite volume and finite temperature effects. We also fit to a larger number of diquark sources to better allow for extrapolation to zero diquark source.
The central idea of this review is to consider quantum field theory models relevant for particle physics and replace the fermionic matter in these models by a bosonic one. This is mostly motivated by the fact that bosons are more ``accessible'' and easier to manipulate for experimentalists, but this ``substitution'' also leads to new physics and novel phenomena. It allows us to gain new information about among other things confinement and the dynamics of the deconfinement transition. We will thus consider bosons in dynamical lattices corresponding to the bosonic Schwinger or Z$_2$ Bose-Hubbard models. Another central idea of this review concerns atomic simulators of paradigmatic models of particle physics theory such as the Creutz-Hubbard ladder, or Gross-Neveu-Wilson and Wilson-Hubbard models. Finally, we will briefly describe our efforts to design experimentally friendly simulators of these and other models relevant for particle physics.
We study two-color QCD with two flavors of Wilson fermion as a function of quark chemical potential mu and temperature T, for two different lattice spacings and two different quark masses. We find that the quarkyonic region, where the behaviour of the quark number density and the diquark condensate are described by a Fermi sphere of almost free quarks distorted by a BCS gap, extends to larger chemical potentials with decreasing lattice spacing or quark mass. In both cases, the quark number density also approaches its non-interacting value. The pressure at low temperature is found to approach the Stefan-Boltzmann limit from below.
Baryonic correlation functions provide an ideal tool to study parity doubling and chiral symmetry using lattice simulations. We present a study using $2+1$ flavors of anisotropic Wilson clover fermions on the FASTSUM ensembles and find clear evidence that parity doubling emerges in the quark-gluon plasma. This result is confirmed on the level of spectral functions, which are obtained using a MEM reconstruction. We further highlight the importance of Gaussian smearing in this study.
At zero temperature the negative-parity ground states of the nucleon and delta baryons are non-degenerate with the positive-parity partners due to spontaneous breaking of chiral symmetry. However, chiral symmetry is expected to be restored at sufficiently high temperature, in particular when going from the hadronic to the quark-gluon plasma (QGP) phase. This would imply that channels with opposite parity become degenerate. We study the nucleon (spin $1/2$) and $Δ$ (spin $3/2$) baryons in both parity sectors using lattice QCD. The range of temperatures spans both the hadronic and QGP phases. Using the FASTSUM anisotropic $N_f = 2 + 1$ ensembles, we analyze the correlation functions and the spectral functions using respectively exponential fits and the Maximum Entropy Method. We find clear evidence of in-medium effects in the hadronic phase, especially for the negative-parity ground state, and of parity doubling in the QGP phase.
Interacting theories of N relativistic fermion flavors in reducible spinor representations in 2+1 spacetime dimensions are formulated on a lattice using domain wall fermions (DWF), for which a U(2N) global symmetry is recovered in the limit that the wall separation $L_s$ is made large. The Gross-Neveu (GN) model is studied in the large-N limit and an exponential acceleration of convergence to the large-$L_s$ limit is demonstrated if the usual parity-invariant mass $m\barψψ$ is replaced by the U(2N)-equivalent $im_3\barψγ_3ψ$. The GN model and two lattice variants of the Thirring model are simulated for N = 2 using a hybrid Monte Carlo algorithm, and studies made of the symmetry-breaking bilinear condensate and its associated susceptibility, the axial Ward identity, and the mass spectrum of both fermion and meson excitations. Comparisons are made with existing results obtained using staggered fermions. For the GN model a symmetry-breaking phase transition is observed, the Ward identity is recovered, and the spectrum found to be consistent with large-N expectations. There appears to be no obstruction to the study of critical UV fixed-point physics using DWF. For the Thirring model the Ward identity is not recovered, the spectroscopy measurements are inconclusive, and no symmetry breaking is observed all the way up to the effective strong coupling limit. This is consistent with a critical Thirring flavor number $N_c<2$, contradicting earlier staggered fermion results.
In 2+1 dimensions, Dirac fermions in reducible, i.e. four-component representations of the spinor algebra form the basis of many interesting model field theories and effective descriptions of condensed matter phenomena. This paper explores lattice formulations which preserve the global U(2N ) symmetry present in the massless limit, and its breakdown to U(N)xU(N) implemented by three independent and parity-invariant fermion mass terms. I set out generalisations of the Ginsparg-Wilson relation, leading to a formulation of an overlap operator, and explore the remnants of the global symmetries which depart from the continuum form by terms of order of the lattice spacing. I also define a domain wall formulation in 2+1+1d, and present numerical evidence, in the form of bilinear condensate and meson correlator calculations in quenched non-compact QED using reformulations of all three mass terms, to show that U(2N) symmetry is recovered in the limit that the domain-wall separation tends to infinity. The possibility that overlap and domain wall formulations of reducible fermions may coincide only in the continuum limit is discussed.
We present a lattice QCD calculation of the charge diffusion coefficient, the electrical conductivity and various susceptibilities of conserved charges, for a range of temperatures below and above the deconfinement crossover. The calculations include the contributions from up, down and strange quarks. We find that the diffusion coefficient is of the order of 1/(2πT) and has a dip around the crossover temperature. Our results are obtained with lattice simulations containing 2+1 dynamical flavours on anisotropic lattices. The Maximum Entropy Method is used to construct spectral functions from correlators of the conserved vector current.
The equivalence of domain wall and overlap fermion formulations is demonstrated for lattice gauge theories in 2+1 spacetime dimensions with parity-invariant mass terms. Even though the domain wall approach distinguishes propagation along a third direction with projectors ${1\over2}(1\pmγ_3)$, the truncated overlap operator obtained for finite wall separation $L_s$ is invariant under interchange of $γ_3$ and $γ_5$. In the limit $L_s\to\infty$ the resulting Ginsparg-Wilson relations recover the expected U($2N_f$) global symmetry up to O($a$) corrections. Finally it is shown that finite-$L_s$ corrections to bilinear condensates associated with dynamical mass generation are characterised by whether even powers of the symmetry-breaking mass are present; such terms are absent for antihermitian bilinears such as $i\barψγ_3ψ$, markedly improving the approach to the large-$L_s$ limit.
We introduce effective field theories for the electronic properties of graphene in terms of relativistic fermions propagating in 2+1 dimensions, and outline how strong inter-electron interactions may be modelled by numerical simulation of a lattice field theory. For strong enough coupling an insulating state can form via condensation of particle-hole pairs, and it is demonstrated that this is a theoretical possibility for monolayer graphene. For bilayer graphene the effect of an interlayer bias voltage can be modelled by the introduction of a chemical potential (akin to isopsin chemical potential in QCD) with no accompanying sign problem; simulations reveal the presence of strong interactions among the residual degrees of freedom at the resulting Fermi surface, which is disrupted by an excitonic condensate. We also present preliminary results for the quasiparticle dispersion, which permit direct estimates of both the Fermi momentum and the induced gap.
The properties of the ground state of two-color QCD at non-zero baryon chemical potential $μ$ present an interesting problem in strongly-interacting gauge theory; in particular the nature of the physically-relevant degrees of freedom in the superfluid phase in the post-onset regime $μ>m_π/2$ still needs clarification. In this study we present evidence for in-medium effects at high $μ$ by studying the wavefunctions of mesonic and diquark states using orthodox lattice simulation techniques, made possible by the absence of a Sign Problem for the model with $N_f=2$. Our results show that beyond onset the spatial extent of hadrons decreases as $μ$ grows, and that the wavefunction profiles are consistent with the existence of a dynamically-gapped Fermi surface in this regime.
It is expected that nucleons and their parity partners become degenerate when chiral symmetry is restored. We investigate this question in the context of the thermal transition from the hadronic phase to the quark-gluon plasma, using lattice QCD simulations with N_f=2+1 flavours. We observe a clear sign of parity doubling in the quark-gluon plasma. Besides, we find that the nucleon ground state is, within the uncertainty, largely independent of the temperature, whereas temperature effects are substantial in the negative-parity (N^*) channel, already in the confined phase.
The spectrum of nucleons and their parity partners is studied as a function of temperature spanning the deconfinement transition. We analyse our results using the correlation functions directly, exponential fits in the hadronic phase, and the Maximum Entropy Method. These techniques all indicate that there is degeneracy in the parity partners' channels in the deconfined phase. This is in accordance with the expectation that there is parity doubling and chiral symmetry in the deconfined phase. In the hadronic phase, we also find that the nucleon ground state is largely independent of temperature, whereas there are substantial temperature effects in the negative parity channel. All results are obtained using our FASTSUM 2+1 flavour ensembles.
Monte Carlo simulation of a 2+1 dimensional model of voltage-biased bilayer graphene, consisting of relativistic fermions with chemical potential mu coupled to charged excitations with opposite sign on each layer, has exposed non-canonical scaling of bulk observables near a quantum critical point found at strong coupling. We present a calculation of the quasiparticle dispersion relation E(k) as a function of exciton source j in the same system, employing partially twisted boundary conditions to boost the number of available momentum modes. The Fermi momentum k_F and superfluid gap Delta are extracted in the limit j tends to zero for three different values of mu, and support a strongly interacting scenario at the Fermi surface with Delta of order O(mu). We propose an explanation for the observation mu < k_F in terms of a dynamical critical exponent z < 1.
QCD at high chemical potential has interesting properties such as deconfinement of quarks. Two-color QCD, which enables numerical simulations on the lattice, constitutes a laboratory to study QCD at high chemical potential. Among the interesting properties of two-color QCD at high density is the diquark condensation, for which we present recent results obtained on a finer lattice compared to previous studies. The quark propagator in two-color QCD at non-zero chemical potential is referred to as the Gor'kov propagator. We express the Gor'kov propagator in terms of form factors and present recent lattice simulation results.
The FASTSUM collaboration has initiated a detailed study of thermal QCD using 2+1 flavours of improved Wilson quarks on anisotropic lattices. Spatial volumes of (3fm)^3 and (4fm)^3 are used at fixed cut-off with temperatures ranging from 40 to 350 MeV (corresponding to temporal lattice extents of 128 to 16 lattice units). Results presented here include the deconfinement temperature and a study of the restoration of chiral symmetry, together with a brief summary of our collaboration's other results from these ensembles.
The FASTSUM collaboration has been carrying out simulations of N_f=2+1 QCD at nonzero temperature in the fixed-scale approach using anisotropic lattices. Here we present the status of these studies, including recent results for electrical conductivity and charge diffusion, and heavy quarkonium (charm and beauty) physics.
We give an update of results from two-colour, two-flavour QCD. Using a Wilson fermion action we calculate thermodynamic quantities as a function of chemical potential μ. Calculating the Karsch Coefficients non-perturbatively gives us access to the derivative method. Compared to our previously published results, we have improved our analysis leading to revised and more accurate estimates for the renormalised energy density, pressure and the trace anomaly.
The FASTSUM Collaboration presents its first results for the electric charge susceptibility in QCD using 2+1 dynamical flavours of Wilson quark on anisotropic lattices. Spatial volumes of (2.94 fm)^3 are used at fixed cut-off with temperatures ranging from below T_c to ~2 T_c.
A lattice calculation is presented for the electrical conductivity of the QCD plasma with 2+1 dynamical flavours at nonzero temperature. We employ the conserved lattice current on anisotropic lattices using a tadpole-improved clover action and study the behaviour of the conductivity over a wide range of temperatures, both below and above the deconfining transition. The conductivity is extracted from a spectral-function analysis using the Maximal Entropy Method and a discussion of its systematics is provided.
We formulate a model of N_f=4 flavors of relativistic fermion in 2+1d in the presence of a chemical potential mu coupled to two flavor doublets with opposite sign, akin to isopsin chemical potential in QCD. This is argued to be an effective theory for low energy electronic excitations in bilayer graphene, in which an applied voltage between the layers ensures equal populations of particles on one layer and holes on the other. The model is then reformulated on a spacetime lattice using staggered fermions, and in the absence of a sign problem, simulated using an orthodox hybrid Monte Carlo algorithm. With the coupling strength chosen to be close to a quantum critical point believed to exist for N_f<N_fc\approx4.8, it is found that there is a region below saturation where both the carrier density and a particle-hole "excitonic" condensate scale anomalously with increasing mu, much more rapidly that the corresponding quantities in free field theory, while the conventional chiral condensate is strongly suppressed. The corresponding ground state is speculated to be a strongly-correlated degenerate fermion system, with a remnant Fermi surface distorted by a superfluid excitonic condensate. The model thus shows qualitatively different behaviour to any model with mu=/=0 previously studied by lattice simulation.
The FASTSUM Collaboration has calculated several quantities relevant for QCD studies at non-zero temperature using the lattice technique. We report here our results for the (i) interquark potential in charmonium; (ii) bottomonium spectral functions; and (iii) electrical conductivity. All results were obtained with 2+1 flavours of dynamical fermions on an anisotropic lattice which allows greater resolution in the temporal direction.
We study two-color QCD with two flavors of Wilson fermion as a function of quark chemical potential mu and temperature T. We find evidence of a superfluid phase at intermediate mu and low T where the quark number density and diquark condensate are both very well described by a Fermi sphere of nearly-free quarks disrupted by a BCS condensate. Our results suggest that the quark contribution to the energy density is negative (and balanced by a positive gauge contribution), although this result is highly sensitive to details of the energy renormalisation. We also find evidence that the chiral condensate in this region vanishes in the massless limit. This region gives way to a region of deconfined quark matter at higher T and mu, with the deconfinement temperature, determined from the renormalised Polyakov loop, decreasing only very slowly with increasing chemical potential. The quark number susceptibility chi_q does not exhibit any qualitative change at the deconfinement transition. We argue that this is because chi_q is not an appropriate measure of deconfinement for 2-color QCD at high density.
The FASTSUM collaboration presents a study on the temperature dependence of the electrical conductivity $σ$ in the quark-gluon plasma, using the methods of lattice QCD. Correlators of the exactly conserved vector current are measured at different temperatures across the deconfinement transition, using ensembles of $2+1$ flavours of dynamical fermions on anisotropic lattices. We then employ bayesian methods (MEM) to extract the relevant spectral functions, which are found to be consistent with $σ/T$ rising as a function of $T$. The robustness of the results is verified by a detailed analysis of the systematics involved in the bayesian reconstruction of the spectral functions.
The heavy quarkonium spectrum of Two Color QCD (QC$_2$D) at non-zero quark chemical potential $μ$ and temperature $T$ with $μ/T\gg1$ has been calculated in both $S$- and $P$-wave channels using a lattice non-relativistic formulation of QC$_2$D. As $μ$ is varied, the quarkonium spectra reveal three separate regions, corroborating previous findings that there are three distinct physical regimes of QC$_2$D at low temperature and high baryon density: hadronic matter, quark/quarkyonic matter, and deconfined matter. The results are interpreted in terms of the formation of heavy-light $Qq$ states in the two-color baryonic medium.
We study two-color QCD with two flavors of Wilson fermion as a function of quark chemical potential mu and temperature T. We find evidence of a superfluid phase at intermediate mu and low T where the quark number density and diquark condensate are both very well described by a Fermi sphere of nearly-free quarks disrupted by a BCS condensate. This gives way to a region of deconfined quark matter at higher T and mu, with the deconfinement temperature decreasing only very slowly with increasing chemical potential. We find that heavy quarkonium bound states persist in the S-wave channels at all T and mu, with an energy reflecting the phase structure. P-wave states appear not to survive in the quarkyonic region.
I present recent results from lattice simulations of SU(2) gauge theory with Nf=2 Wilson quark flavors, at non-zero quark chemical potential. The thermodynamic equation of state is discussed along with the nature of the high density matter which forms. It is conjectured that deconfinement may mean different things for bulk and Fermi surface phenomena.
We present results from a simulation of SU(2) lattice gauge theory with N_f=4 flavors of Wilson fermion and non-zero quark chemical potential mu, using the same 12^3x24 lattice, bare gauge coupling, and pion mass in cut-off units as a previous study with N_f=2. The string tension for N_f=4 is found to be considerably smaller implying smoother gauge field configurations. Thermodynamic observables and order parameters for superfluidity and color deconfinement are studied, and comparisons drawn between the two theories. Results for quark density and pressure as functions of mu are qualitatively similar for N_f=2 and N_f=4; in both cases there is evidence for a phase in which baryonic matter is simultaneously degenerate and confined. Results for the stress-energy tensor, however, suggest that while N_f=2 has a regime where dilute matter is non-relativistic and weakly-interacting, N_f=4 matter is relativistic and strongly-interacting for all values of mu above onset.
We present the first evidence from lattice simulations that the magnetic monopoles in three dimensional compact quantum electrodynamics (cQED3) with N_f=2 and N_f= 4 four-component fermion flavors are in a plasma phase. The evidence is based mainly on the divergence of the monopole susceptibility (polarizability) with the lattice size at weak gauge couplings. A weak four-Fermi term added to the cQED3 action enabled simulations with massless fermions. The exact chiral symmetry of the interaction terms forbids symmetry breaking lattice discretization counterterms to appear in the theory's effective action. It is also shown that the scenario of a monopole plasma does not depend on the strength of the four-Fermi coupling. Other observables such as the densities of "isolated" dipoles and monopoles and the so-called specific heat show that a crossover from a dense monopole plasma to a dilute monopole gas occurs at strong couplings. The implications of our results on the stability of U(1) spin liquids in two spatial dimensions are also discussed.
We present results from lattice simulations of a monolayer graphene model at non-zero temperature. At low temperatures for sufficiently strong coupling the model develops an excitonic condensate of particle-hole pairs corresponding to an insulating phase. The Berezinskii-Kosterlitz-Thouless phase transition temperature is associated with the value of the coupling where the critical exponent delta governing the response of the order parameter at criticality to an external source has a value close to 15. The critical coupling on a lattice with temporal extent N_t=32 (T=1/(N_t a_t) where a_t is the temporal lattice spacing) and spatial extent N_s=64 is very close to infinite coupling. The value of the transition temperature normalized with the zero temperature fermion mass gap Delta_0 is given by T_BKT/Delta_0=0.055(2). This value provides an upper bound on the transition temperature, because simulations closer to the continuum limit where the full U(4) symmetry is restored may result in an even lower value. In addition, we measured the helicity modulus Upsilon and the fermion thermal mass Delta_T(T), the later providing evidence for a pseudogap phase with Delta_T>0 extending to arbitrarily high T.
We study the quark number susceptibility in SU(2) lattice gauge theory with two Wilson quark flavours at non-zero chemical potential and low temperature. We present some technical aspects of the issue and numerical results obtained at different lattices and different parameters. We discuss what this observable can teach us about the phase diagram of the model and in particular about the relation between susceptibility and Polyakov loop
We study the topological charge fluctuations of an SU(2) lattice gauge theory containing both N_f=2 and 4 flavors of Wilson fermion, at low temperature with non-zero chemical potential $μ$. The topological susceptibility, chi_T, is used to characterize differing physical regimes as mu is varied between the onset of matter at mu_o and and color deconfinement at mu_d. Suppression of instantons by matter via Debye screening is also investigated, revealing effects not captured by perturbative predictions. In particular, the breaking of scale invariance leads to the mean instanton size rho becoming mu-dependent in the regime between onset and deconfinement, with a scaling rho~1/mu^2 over the range mu_o<mu<mu_d, resulting in an enhancement of chi_T immediately above onset.
We present results from simulations of Two Color QCD with two Wilson quark flavors in the presence of a quark chemical potential mu at two different lattice spacings. The equation of state, conformal anomaly, superfluid order parameter and Polyakov line are all discussed. Our results suggest that the transition from hadronic to quark matter, and that from confined to deconfined matter occur at distinct values of mu, consistent with the existence of a quarkyonic phase in this model.
It has long been known that QCD undergoes a deconfining phase transition at high temperature. One of the consequent features of this new, quark-gluon phase is that hadrons become unbounded. In this talk meson correlation functions at non-zero momentum are studied in the deconfined phase using the Maximum Entropy Method.
We consider QCD at very low temperatures and non-zero quark chemical potential from lattice Monte Carlo simulations of the two-color theory in a very small spatial volume (the attoscale). In this regime the quark number rises in discrete levels in qualitative agreement with what is found analytically at one loop on S3xS1 with radius R_S3 << 1/Λ_QCD. The detailed level degeneracy, however, cannot be accounted for using weak coupling arguments. At each rise in the quark number there is a corresponding spike in the Polyakov line, also in agreement with the perturbative results. In addition the quark number susceptibility shows a similar behaviour to the Polyakov line and appears to be a good indicator of a confinement-deconfinement type of transition.
To leading order in perturbation theory, we solve QCD, defined on a small three sphere in the large N and Nf limit, at finite chemical potential and map out the phase diagram in the (mu,T) plane. The action of QCD is complex in the presence of a non-zero quark chemical potential which results in the sign problem for lattice simulations. In the large N theory, which at low temperatures becomes a conventional unitary matrix model with a complex action, we find that the dominant contribution to the functional integral comes from complexified gauge field configurations. For this reason the eigenvalues of the Polyakov line lie off the unit circle on a contour in the complex plane. We find at low temperatures that as mu passes one of the quark energy levels there is a third-order Gross-Witten transition from a confined to a deconfined phase and back again giving rise to a rich phase structure. We compare a range of physical observables in the large N theory to those calculated numerically in the theory with N=3. In the latter case there are no genuine phase transitions in a finite volume but nevertheless the observables are remarkably similar to the large N theory.
In this proceedings we summarize our calculation of the phase diagram of QCD at non-zero temperature and chemical potential on S^1 x S^3 from one-loop perturbation theory [1], which is valid in the limit R << 1/Lambda, where R is the radius of S^3. We calculate several observables including the Polyakov lines and the quark number, for large number of colors N and large number of quark flavors Nf, on S^1 x S^3, and compare with results for the same system with N = 3, and with results for N=2 lattice QCD. For N > 2 the action is complex and the dominant contributions to the path integral occur in the space of complexified gauge field configurations. This results in the expectation values of the eigenvalues of the Polyakov line lying off the unit circle and out in the complex plane. This is an important issue for the lattice, and also for the calculation on S^1 x S^3 in the large N limit where we obtain analytical results using the saddle point approximation. It is thus necessary to adapt available techniques to locate the stationary solutions in the complexified gauge field configuration space.
A 2+1 dimensional fermion field theory is proposed as a model for the low-energy electronic excitations in monolayer graphene. The model consists of N=2 four-component Dirac fermions moving in the plane and interacting via a contact interaction between charge densities. For strong couplings there is a continuous transition to a Mott insulting phase. We present results of an extensive numerical study of the model's critical region, including the order parameter, its associated susceptibility, and for the first time the quasiparticle propagator. The data enables an extraction of the critical exponents at the transition, including the dynamical critical exponent, which are hypothesised to be universal features of a quantum critical point. The relation of our model with others in the literature is discussed, along with the implications for physical graphene following from our value of the critical coupling.
We formulate a model of relativistic fermions moving in two Euclidean dimensions based on a tight-binding model of graphene. The eigenvalue spectrum of the resulting Dirac operator is solved numerically in smooth U(1) gauge field backgrounds carrying an integer-valued topological charge Q, and it is demonstrated that the resulting number of zero-eigenvalue modes is in accord with the Atiyah-Singer index theorem applied to two continuum flavors. A bilinear but gauge non-invariant chirality operator appropriate for distinguishing the topological zero modes is identified. When this operator is used to calculate Q, it is found that the maximum topological charge capable of being measured in this fashion scales with the perimeter of the lattice. Some concluding remarks compare these results to what is known for staggered lattice fermions.
QED3 with fermi and gap anisotropies is considered to be a candidate effective field theory of high temperature superconductors. Simulations of a variant of the theory have demonstrated that there is evidence consistent with a phase transition to a chirally restored phase as the velocity anisotropy increases, and that the correlation lengths of pions in spatial directions behave in a manner consistent with the anisotropy being a relevant parameter of the theory. We present the first measurements of the fermion propagator in Landau gauge for this theory and suggest that the structure of this theory is more complex than previously thought.
In the high-temperature phase of QCD meson spectral functions at nonzero momentum are expected to have a nontrivial and interesting structure. In order to provide a reference point for lattice studies employing e.g. the Maximal Entropy Method, we discuss several characteristics of meson spectral functions in the infinite-temperature limit. We report on ongoing work in quenched QCD with staggered fermions.
We report on a lattice study of two-color QCD with adjoint staggered fermions at high density. We find that the model has no early onset and we report on results for diquark condensation, from simulations with and without a diquark source term.
I review recent developments in the study of strongly interacting field theories with non-zero chemical potential mu. In particular I focus on (a) the determination of the QCD critical endpoint in the (mu,T) plane; (b) superfluid condensates in Two Color QCD; and (c) Fermi surface effects in the NJL model. Some remarks are made concerning the relation of superconductivity with the sign problem.
Four-fermi models in dimensionality $2<d<4$ exhibit an ultra-violet stable renormalization group fixed point at a strong value of the coupling constant where chiral symmetry is spontaneously broken. The resulting field theory describes relativistic fermions interacting non-trivially via exchange of scalar bound states. We calculate the $O(1/N_f)$ corrections to this picture, where $N_f$ is the number of fermion species, for a variety of models and confirm their renormalizability to this order. A connection between renormalizability and the hyperscaling relations between the theory's critical exponents is made explicit. We present results of extensive numerical simulations of the simplest model for $d=3$, performed using the hybrid Monte Carlo algorithm on lattice sizes ranging from $8^3$ to $24^3$. For $N_f=12$ species of massless fermions we confirm the existence of a second order phase transition where chiral symmetry is spontaneously broken. Using both direct measurement and finite size scaling arguments we estimate the critical exponents $β$, $γ$, $ν$ and $δ$. We also investigate symmetry restoration at non-zero temperature, and the scalar two-point correlation function in the vicinity of the bulk transition. All our results are in excellent agreement with analytic predictions, and support the contention that the $1/N_f$ expansion is accurate for this class of models.