Refine
Year of publication
Document Type
- Conference Proceeding (35)
- Article (17)
Language
- English (52)
Has Fulltext
- yes (52)
Is part of the Bibliography
- no (52)
Keywords
- Lattice QCD (3)
- Lattice Quantum Field Theory (2)
- Effective Field Theories (1)
- FOS: Physical sciences (1)
- Finite baryon density (1)
- High Energy Physics - Lattice (hep-lat) (1)
- High Energy Physics - Phenomenology (hep-ph) (1)
- High Energy Physics - Theory (hep-th) (1)
- Hybrid Monte Carlo algorithm (1)
- Lattice field theory (1)
Institute
- Physik (52)
- Frankfurt Institute for Advanced Studies (FIAS) (4)
- ELEMENTS (2)
- Informatik (1)
Lattice Yang-Mills theories at finite temperature can be mapped onto effective 3d spin systems, thus facilitating their numerical investigation. Using strong-coupling expansions we derive effective actions for Polyakov loops in the SU(2) and SU(3) cases and investigate the effect of higher order corrections. Once a formulation is obtained which allows for Monte Carlo analysis, the nature of the phase transition in both classes of models is investigated numerically, and the results are then used to predict – with an accuracy within a few percent – the deconfinement point in the original 4d Yang-Mills pure gauge theories, for a series of values of Nt at once.
Perturbation theory for non-abelian gauge theories at finite temperature is plagued by infrared
divergences which are caused by magnetic soft modes ~ g2T, corresponding to gluon fields of
a 3d Yang-Mills theory. While the divergences can be regulated by a dynamically generated
magnetic mass on that scale, the gauge coupling drops out of the effective expansion parameter
requiring summation of all loop orders for the calculation of observables. Some gauge invariant
possibilities to implement such infrared-safe resummations are reviewed. We use a scheme based
on the non-linear sigma model to estimate some of the contributions ~ g6 of the soft magnetic
modes to the QCD pressure through two loops. The NLO contribution amounts to ~ 10% of the
LO, suggestive of a reasonable convergence of the series.
The so-called sign problem of lattice QCD prohibits Monte Carlo simulations at finite baryon
density by means of importance sampling. Over the last few years, methods have been developed
which are able to circumvent this problem as long as the quark chemical potential is m=T <~1.
After a brief review of these methods, their application to a first principles determination of the
QCD phase diagram for small baryon densities is summarised. The location and curvature of the
pseudo-critical line of the quark hardon transition is under control and extrapolations to physical
quark masses and the continuum are feasible in the near future. No definite conclusions can as
yet be drawn regarding the existence of a critical end point, which turns out to be extremely quark
mass and cut-off sensitive. Investigations with different methods on coarse lattices show the lightmass
chiral phase transition to weaken when a chemical potential is switched on. If persisting on
finer lattices, this would imply that there is no chiral critical point or phase transition for physical
QCD. Any critical structure would then be related to physics other than chiral symmetry breaking.
The chiral critical surface is a surface of second order phase transitions bounding the region of
first order chiral phase transitions for small quark masses in the fmu;d;ms;mg parameter space.
The potential critical endpoint of the QCD (T;m)-phase diagram is widely expected to be part of
this surface. Since for m = 0 with physical quark masses QCD is known to exhibit an analytic
crossover, this expectation requires the region of chiral transitions to expand with m for a chiral
critical endpoint to exist. Instead, on coarse Nt = 4 lattices, we find the area of chiral transitions
to shrink with m, which excludes a chiral critical point for QCD at moderate chemical potentials
mB < 500 MeV. First results on finer Nt = 6 lattices indicate a curvature of the critical surface
consistent with zero and unchanged conclusions. We also comment on the interplay of phase
diagrams between the Nf = 2 and Nf = 2+1 theories and its consequences for physical QCD.
We report progress in our exploration of the finite-temperature phase structure of two-flavour lattice
QCD with twisted-mass Wilson fermions and a tree-level Symanzik-improved gauge action
for a temporal lattice size Nt = 8. Extending our investigations to a wider region of parameter
space we gain a global view of the rich phase structure. We identify the finite temperature transition/
crossover for a non-vanishing twisted-mass parameter in the neighbourhood of the zerotemperature
critical line at sufficiently high b . Our findings are consistent with Creutz’s conjecture
of a conical shape of the finite temperature transition surface. Comparing with NLO lattice
cPT we achieve an improved understanding of this shape.
We discuss the use of Wilson fermions with twisted mass for simulations of QCD thermodynamics.
As a prerequisite for a future analysis of the finite-temperature transition making use
of automatic O(a) improvement, we investigate the phase structure in the space spanned by the
hopping parameter k , the coupling b , and the twisted mass parameter m. We present results for
Nf = 2 degenerate quarks on a 163×8 lattice, for which we investigate the possibility of an Aoki
phase existing at strong coupling and vanishing m, as well as of a thermal phase transition at
moderate gauge couplings and non-vanishing m.
Pseudo-Critical Temperature and Thermal Equation of State from Nf = 2 Twisted Mass Lattice QCD
(2012)
We report about the current status of our ongoing study of the chiral limit of two-flavor QCD at finite temperature with twisted mass quarks. We estimate the pseudo-critical temperature Tc for three values of the pion mass in the range of mPS ~ 300 and 500 MeV and discuss different chiral scenarios. Furthermore, we present first preliminary results for the trace anomaly, pressure and energy density. We have studied several discretizations of Euclidean time up to Nt = 12 in order to assess the continuum limit of the trace anomaly. From its interpolation we evaluate the pressure and energy density employing the integral method. Here, we have focussed on two pion masses with mPS ~ 400 and 700 MeV.
We review our knowledge of the phase diagram of QCD as a function of temperature, chemical potential and quark masses. The presence of tricritical lines at imaginary chemical potential m = i p 3 T, with known scaling behaviour in their vicinity, puts constraints on this phase diagram, especially in the case of two light flavors. We show first results in our project to determine the finite-temperature behaviour in the Nf = 2 chiral limit.
We compute hybrid static potentials in SU(3) lattice gauge theory. We present a method to automatically generate a large set of suitable creation operators with defined quantum numbers from elementary building blocks. We show preliminary results for several channels and discuss, which structures of the gluonic flux tube seem to be realized by the ground states in these channels.
The properties of matter at finite baryon densities play an important role for the astrophysics of compact stars as well as for heavy ion collisions or the description of nuclear matter. Because of the sign problem of the quark determinant, lattice QCD cannot be simulated by standard Monte Carlo at finite baryon densities. I review alternative attempts to treat dense QCD with an effective lattice theory derived by analytic strong coupling and hopping expansions, which close to the continuum is valid for heavy quarks only, but shows all qualitative features of nuclear physics emerging from QCD. In particular, the nuclear liquid gas transition and an equation of state for baryons can be calculated directly from QCD. A second effective theory based on strong coupling methods permits studies of the phase diagram in the chiral limit on coarse lattices.
The quark gluon plasma produced in heavy ion collisions behaves like an almost ideal fluid described by viscous hydrodynamics with a number of transport coefficients. The second order coefficient κ is related to a Euclidean correlator of the energy-momentum tensor at vanishing frequency and low momentum. This allows for a lattice determination without maximum entropy methods or modelling, but the required lattice sizes represent a formidable challenge. We calculate κ in leading order lattice perturbation theory and simulations on 1203 × 6, 8 lattices with a < 0.1 fm. In the temperature range 2Tc − 10Tc we find κ = 0.36(15)T2. The error covers both a suitably rescaled AdS/CFT prediction as well as, remarkably, the result of leading order perturbation theory. This suggests that appropriate noise reduction methods on the lattice and NLO perturbative calculations could provide an accurate QCD prediction in the near future.
Lattice simulations employing reweighting and Taylor expansion techniques have predicted a (m;T)-phase diagram according to general expectations, with an analytic quark-hadron crossover at m =0 turning into a first order transition at some critical chemical potential mE. By contrast, recent simulations using imgainary m followed by analytic continuation obtained a critical structure in the fmu;d;ms;T;mg parameter space favouring the absence of a critical point and first order line. I review the evidence for the latter scenario, arguing that the various raw data are not inconsistent with each other. Rather, the discrepancy appears when attempting to extract continuum results from the coarse (Nt =4) lattices simulated so far, and can be explained by cut-off effects. New (as yet unpublished) data are presented, which for Nf = 3 and on Nt = 4 confirm the scenario without a critical point. Moreover, simulations on finer Nt = 6 lattices show that even if there is a critical point, continuum extrapolation moves it to significantly larger values of mE than anticipated on coarse lattices.
Euclidean strong coupling expansion of the partition function is applied to lattice Yang-Mills theory
at finite temperature, i.e. for lattices with a compactified temporal direction. The expansions
have a finite radius of convergence and thus are valid only for b <bc, where bc denotes the nearest
singularity of the free energy on the real axis. The accessible temperature range is thus the
confined regime up to the deconfinement transition. We have calculated the first few orders of
these expansions of the free energy density as well as the screening masses for the gauge groups
SU(2) and SU(3). The resulting free energy series can be summed up and corresponds to a glueball
gas of the lowest mass glueballs up to the calculated order. Our result can be used to fix
the lower integration constant for Monte Carlo calculations of the thermodynamic pressure via
the integral method, and shows from first principles that in the confined phase this constant is
indeed exponentially small. Similarly, our results also explain the weak temperature dependence
of glueball screening masses below Tc, as observed in Monte Carlo simulations. Possibilities and
difficulties in extracting bc from the series are discussed.
We report on the first steps of an ongoing project to add gauge observables and gauge corrections
to the well-studied strong coupling limit of staggered lattice QCD, which has been shown earlier
to be amenable to numerical simulations by the worm algorithm in the chiral limit and at finite
density. Here we show how to evaluate the expectation value of the Polyakov loop in the framework
of the strong coupling limit at finite temperature, allowing to study confinement properties
along with those of chiral symmetry breaking. We find the Polyakov loop to rise smoothly, thus
signalling deconfinement. The non-analytic nature of the chiral phase transition is reflected in the
derivative of the Polyakov loop. We also discuss how to construct an effective theory for non-zero
lattice coupling, which is valid to O(b).
We study a random matrix model for QCD at finite density via complex Langevin dynamics. This model has a phase transition to a phase with nonzero baryon density. We study the convergence of the algorithm as a function of the quark mass and the chemical potential and focus on two main observables: the baryon density and the chiral condensate. For simulations close to the chiral limit, the algorithm has wrong convergence properties when the quark mass is in the spectral domain of the Dirac operator. A possible solution of this problem is discussed.
I review recent developments in determining the QCD phase diagram by means of lattice simulations.
Since the invention of methods to side-step the sign problem a few years ago, a number
of additional variants have been proposed, and progress has been made towards understanding
some of the systematics involved. All available techniques agree on the transition temperature
as a function of density in the regime mq/T <~1. There are by now four calculations with signals
for a critical point, two of them at similar parameter values and with consistent results. However,
it also emerges that the location of the critical point is exceedingly quark mass sensitive. At the
same time sizeable finite volume, cut-off and step size effects have been uncovered, demanding
additional investigations with exact algorithms on larger and finer lattices before quantitative conclusions
can be drawn. Depending on the sign of these corrections, there is ample room for the
eventual phase diagram to look as expected or also quite different, with no critical point at all.
QCD at finite temperature and denisty remains intractable by Monte Carlo simulations for quark
chemical potentials m >∼T. It has been a long standing problem to derive effective theories from
QCD which describe the phase structure of the former with controlled errors. We propose a
solution to this problem by a combination of analytical and numerical methods. Starting from
lattice QCD with in Wilson’s formulation, we derive an effective action in terms of Polyakov
loops by means of combined strong coupling and hopping expansions. The theory correctly
reflects the centre-symmetry in the pure gauge limit and its breaking through quarks. It is valid
for heavy quarks and lattices up to Nt ∼ 6. Its sign problem can be solved and we are able to
calculate the deconfinement transition of QCD with heavy quarks for all chemical potentials.
We extend the recently developed strong coupling, dimensionally reduced Polyakov-loop effective theory from finite-temperature pure Yang-Mills to include heavy fermions and nonzero chemical
potential by means of a hopping parameter expansion. Numerical simulation is employed to investigate the weakening of the deconfinement transition as a function of the quark mass. The
tractability of the sign problem in this model is exploited to locate the critical surface in the (M/T,m/T,T) space over the whole range of chemical potentials from zero up to infinity.
LatticeQCD using OpenCL
(2011)
We report on the status of ongoing investigations aiming at locating the deconfinement critical point with standard Wilson fermions and Nf = 2 flavors towards the continuum limit (standard Columbia plot); locating the tricritical masses at imaginary chemical potential with unimproved staggered fermions at Nf = 2 (extended Columbia plot); identifying the order of the chiral phase transition at μ = 0 for Nf = 2 via extrapolation from non integer Nf (alternative Columbia plot).
The QCD phase diagram as a function of temperature, T, and chemical potential for baryon
number, mB, is still unknown today, due to the sign problem, which prohibits direct Monte Carlo
simulations for non-vanishing baryon density. Investigations in models sharing chiral symmetry
with QCD predict a phase diagram, in which the transition corresponds to a smooth crossover at
zero density, but which is strengthened by chemical potential to turn into a first order transition
beyond some second order critical point. This contribution reviews the lattice evidence in favour
and against the existence of a critical point.
From the colour glass condensate to filamentation: systematics of classical Yang–Mills theory
(2019)
The non-equilibrium early time evolution of an ultra-relativistic heavy ion collision is often described by classical lattice Yang–Mills theory, starting from the colour glass condensate (CGC) effective theory with an anisotropic energy momentum tensor as initial condition. In this work we investigate the systematics associated with such studies and their dependence on various model parameters (IR, UV cutoffs and the amplitude of quantum fluctuations) which are not yet fixed by experiment. We perform calculations for SU() and SU(), both in a static box and in an expanding geometry. Generally, the dependence on model parameters is found to be much larger than that on technical parameters like the number of colours, boundary conditions or the lattice spacing. In a static box, all setups lead to isotropisation through chromo-Weibel instabilities, which is illustrated by the accompanying filamentation of the energy density. However, the associated time scale depends strongly on the model parameters and in all cases is longer than the phenomenologically expected one. In the expanding system, no isotropisation is observed for any parameter choice. We show how investigations at fixed initial energy density can be used to better constrain some of the model parameters.
For a long time, strong coupling expansions have not been applied systematically in lattice QCD thermodynamics, in view of the success of numerical Monte Carlo studies. The persistent sign problem at finite baryo-chemical potential, however, has motivated investigations using these methods, either by themselves or combined with numerical evaluations, as a route to finite density physics. This article reviews the strategies, by which a number of qualitative insights have been attained, notably the emergence of the hadron resonance gas or the identification of the onset transition to baryon matter in specific regions of the QCD parameter space. For the simpler case of Yang–Mills theory, the deconfinement transition can be determined quantitatively even in the scaling region, showing possible prospects for continuum physics.
Attempts to extract the order of the chiral transition of QCD at zero chemical potential, with two dynamical flavors of massless quarks, from simulations with progressively decreasing pion mass, have remained inconclusive because of their increasing numerical cost. In an alternative approach to this problem, we consider the path integral as a function of continuous number Nf of degenerate quarks. If the transition in the chiral limit is first order for Nf≥3, a second-order transition for Nf=2 then requires a tricritical point in between. This, in turn, implies tricritical scaling of the critical boundary line between the first-order and crossover regions as the chiral limit is approached. Noninteger numbers of fermion flavors are easily implemented within the staggered fermion discretization. Exploratory simulations at μ=0 and Nf=2.8, 2.6, 2.4, 2.2, 2.1, on coarse Nτ=4 lattices, indeed show a smooth variation of the critical mass mapping out a critical line in the (m, Nf) plane. For the smallest masses, the line appears consistent with tricritical scaling, allowing for an extrapolation to the chiral limit.
Thermal leptogenesis, in the framework of the standard model with three additional heavy Majorana neutrinos, provides an attractive scenario to explain the observed baryon asymmetry in the universe. It is based on the out-of-equilibrium decay of Majorana neutrinos in a thermal bath of standard model particles, which in a fully quantum field theoretical formalism is obtained by solving Kadanoff-Baym equations. So far, the leading two-loop contributions from leptons and Higgs particles are included, but not yet gauge corrections. These enter at three-loop level but, in certain kinematical regimes, require a resummation to infinite loop order for a result to leading order in the gauge coupling. In this work, we apply such a resummation to the calculation of the lepton number density. The full result for the simplest “vanilla leptogenesis” scenario is by O(1) increased compared to that of quantum Boltzmann equations, and for the first time permits an estimate of all theoretical uncertainties. This step completes the quantum theory of leptogenesis and forms the basis for quantitative evaluations, as well as extensions to other scenarios.
Lattice QCD and functional methods are making significant progress in constraining the QCD phase diagram. As an important milestone, the chiral phase transition with massless u, d-quarks at zero density is now understood to be of second order for all strange quark masses, and a smooth crossover as soon as mu,d, ≠ 0. Together with information on fluctuations and refined reweighted simulations, this bounds a possible critical point to be at µB/T ≲3. On the other hand, an approximately chiral-spin symmetric temperature window has been discovered above the chiral crossover, Tch<T ≳3Tch, with distinct correlator multiplet patterns and a pion spectral function suggesting resonance-like degrees of freedom, which dissolve graduallly with temperature.
The QCD equation of state is not often discussed in cosmology. However, the relic density of
weakly interacting massive particles (WIMPs) depends on the entropy and the expansion rate of
the Universe when they freeze out, at a temperature in the range 400 MeV – 40GeV, where QCD
corrections are still important. We use recent analytic and lattice calculations of the QCD pressure
to produce a new equation of state suitable for use in relic density calculations. As an example,
we show that relic densities calculated by the dark matter package DarkSUSY receive corrections
of several per cent, within the observational accuracy of the Planck CMB mission, due for launch
in 2007.
We present unambiguous evidence from lattice simulations of Nf = 3 QCD for two tricritical points in the (T;m) phase diagram at fixed imaginary m=T = ip=3 mod. 2p=3, one in the light and one in the heavy quark regime. Together with similar results in the literature for Nf = 2 this implies the existence of a chiral and of a deconfinement tricritical line at those values of imaginary chemical potentials. These tricritical lines represent the boundaries of the analytically continued chiral and deconfinement critical surfaces, respectively, which delimit the parameter space with first order phase transitions. It is demonstrated that the shape of the deconfinement critical surface is dictated by tricritical scaling and implies the weakening of the deconfinement transition with real chemical potential. A qualitatively similar effect holds for the chiral critical surface.
For the exploration of the phase diagram of QCD, effective Polyakov loop theories derived from lattice QCD provide a valuable tool in the heavy quark mass regime. In practice, the evaluation of these theories is complicated by the appearance of long-range and multipoint interaction terms. On the other hand, it is well known that for theories with such kind of interactions mean field approximations can be expected to yield reliable results. Here, we apply this framework to the critical endpoint of the deconfinement transition and results are compared to the literature. This treatment can also be used to investigate the phase diagram at non-zero baryon and isospin chemical potential.
The SU(3) spin model with chemical potential corresponds to a simplified version of QCD with static quarks in the strong coupling regime. It has been studied previously as a testing ground for new methods aiming to overcome the sign problem of lattice QCD. In this work we show that the equation of state and the phase structure of the model can be fully determined to reasonable accuracy by a linked cluster expansion. In particular, we compute the free energy to 14-th order in the nearest neighbour coupling. The resulting predictions for the equation of state and the location of the critical end points agree with numerical determinations to O(1%) and O(10%), respectively. While the accuracy for the critical couplings is still limited at the current series depth, the approach is equally applicable at zero and non-zero imaginary or real chemical potential, as well as to effective QCD Hamiltonians obtained by strong coupling and hopping expansions.
Lattice QCD with heavy quarks reduces to a three-dimensional effective theory of Polyakov loops, which is amenable to series expansion methods. We analyse the effective theory in the cold and dense regime for a general number of colours, Nc. In particular, we investigate the transition from a hadron gas to baryon condensation. For any finite lattice spacing, we find the transition to become stronger, i.e. ultimately first-order, as Nc is made large. Moreover, in the baryon condensed regime, we find the pressure to scale as p ∼ Nc through three orders in the hopping expansion. Such a phase differs from a hadron gas with p ∼ N0c, or a quark gluon plasma, p ∼ N2c, and was termed quarkyonic in the literature, since it shows both baryon-like and quark-like aspects. A lattice filling with baryon number shows a rapid and smooth transition from condensing baryons to a crystal of saturated quark matter, due to the Pauli principle, and is consistent with this picture. For continuum physics, the continuum limit needs to be taken before the large Nc limit, which is not yet possible in practice. However, in the controlled range of lattice spacings and Nc-values, our results are stable when the limits are approached in this order. We discuss possible implications for physical QCD.
Lattice constraints on the QCD chiral phase transition at finite temperature and baryon density
(2021)
The thermal restoration of chiral symmetry in QCD is known to proceed by an analytic crossover, which is widely expected to turn into a phase transition with a critical endpoint as the baryon density is increased. In the absence of a genuine solution to the sign problem of lattice QCD, simulations at zero and imaginary baryon chemical potential in a parameter space enlarged by a variable number of quark flavours and quark masses constitute a viable way to constrain the location of a possible non-analytic phase transition and its critical endpoint. In this article I review recent progress towards an understanding of the nature of the transition in the massless limit, and its critical temperature at zero density. Combined with increasingly detailed studies of the physical crossover region, current data bound a possible critical point to μB ≳ 3T.
The nature of the QCD chiral phase transition in the limit of vanishing quark masses has remained elusive for a long time, since it cannot be simulated directly on the lattice and is strongly cutoff-dependent. We report on a comprehensive ongoing study using unimproved staggered fermions with Nf ∈ [2, 8] mass-degenerate flavours on Nτ ∈ {4, 6, 8} lattices, in which we locate the chiral critical surface separating regions with first-order transitions from crossover regions in the bare parameter space of the lattice theory. Employing the fact that it terminates in a tricritical line, this surface can be extrapolated to the chiral limit using tricritical scaling with known exponents. Knowing the order of the transitions in the lattice parameter space, conclusions for approaching the continuum chiral limit in the proper order can be drawn. While a narrow first-order region cannot be ruled out, we find initial evidence consistent with a second-order chiral transition in all massless theories with Nf ≤ 6, and possibly up to the onset of the conformal window at 9 ≲ N∗f ≲ 12. A reanalysis of already published O(a)-improved Nf = 3 Wilson data on Nτ ∈ [4, 12] is also consistent with tricritical scaling, and the associated change from first to second-order on the way to the continuum chiral limit. We discuss a modified Columbia plot and a phase diagram for many-flavour QCD that reflect these possible features.
Recently, an approximate SU(4) chiral spin-flavour symmetry was observed in multiplet patterns of QCD meson correlation functions, in a temperature range above the chiral crossover. This symmetry is larger than the chiral symmetry of massless QCD, and can only arise effectively when colour-electric quark-gluon interactions dynamically dominate the quantum effective action. At temperatures about three times the crossover temperature, these patterns disappear again, indicating the screening of colour-electric interactions, and the expected chiral symmetry is recovered. In this contribution we collect independent evidence for such an intermediate temperature range, based on screening masses and the pion spectral function. Both kinds of observables behave non-perturbatively in this window, with resonance-like peaks for the pion and its first excitation disappearing gradually with temperature. Using symmetry arguments and the known behaviour of screening masses at small densities, we discuss how this chiral spin symmetric band continues into the QCD phase diagram.
In the strong coupling and heavy quark mass regime, lattice QCD dimensionally reduces to effective theories of Polyakov loops depending on the parameters of the original Wilson action β,κ and Nτ. We apply coarse graining techniques to such theories in 1d and 2d, corresponding to lattice QCD at finite temperature and non-zero chemical potential in 1+1d and 2+1d, respectively. In 1d the method is applied to the effective theories up to O(κ4). Using the transfer matrix, the recursion relations are solved analytically. The thermodynamic limit is taken for some observables. Afterwards, continuum extrapolation is performed numerically and results are discussed. In 2d the coarse graining method is applied in the pure gauge and static quark limit. Running couplings are obtained and the fixed points of the transformations are discussed. Finally, the critical coupling of the deconfinement transition is determined in both limits. Agreement to about 12% with Monte Carlo results of 2+1d Yang-Mills theory from the literature is observed.
The order of the chiral phase transition of lattice QCD with unimproved staggered fermions is known to depend on the number of quark flavours, their masses and the lattice spacing. Previous studies in the literature for Nf∈{3,4} show first-order transitions, which weaken with decreasing lattice spacing. Here we investigate what happens when lattices are made coarser to establish contact to the strong coupling region. For Nf∈{4,8} we find a drastic weakening of the transition when going from Nτ=4 to Nτ=2, which is consistent with a second-order chiral transition reported in the literature for Nf=4 in the strong coupling limit. This implies a non-monotonic behaviour of the critical quark or pseudo-scalar meson mass, which separates first-order transitions from crossover behaviour, as a function of lattice spacing.
A lot of effort in lattice simulations over the last years has been devoted to studies of the QCD deconfinement transition. Most state-of-the-art simulations use rooted staggered fermions, while Wilson fermions are affected by large systematic uncertainties, such as coarse lattices or heavy sea quarks. Here we report on an ongoing study of the transition, using two degenerate flavours of nonperturbatively O(a) improved Wilson fermions. We start with Nt = 12 and 16 lattices and pion masses of 600 to 450 MeV, aiming at chiral and continuum limits with light quarks.
We explore the phase diagram of two flavour QCD at vanishing chemical potential using dynamical O(a)-improved Wilson quarks. In the approach to the chiral limit we use lattices with a temporal extent of Nt = 16 and spatial extent L = 32;48 and 64 to enable the extrapolation to the thermodynamic limit with small discretisation effects. In addition to an update on the scans at constant k, reported earlier, we present first results from scans along lines of constant physics at a pion mass of 290 MeV.We probe the transition using the Polyakov loop and the chiral condensate, as well as spectroscopic observables such as screening masses.
Spectral functions encode a wealth of information about the dynamics of any given system, and the determination of their non-perturbative characteristics is a long-standing problem in quantum field theory. Whilst numerical simulations of lattice QCD provide ample data for various Euclidean correlation functions, the inversion required to extract spectral functions is an ill-posed problem. In this work, we pursue previously established constraints imposed by field locality at finite temperature T, namely that spectral functions possess a non-perturbative representation which generalises the well-known Källén-Lehmann spectral form to T > 0. Using this representation, we analyse lattice QCD data of the spatial pseudo-scalar correlator in the temperature range 220–960 MeV, and obtain an analytic expression for the corresponding spectral function, with parameters fixed by the data. From the structure of this spectral function we find evidence for the existence of a distinct pion state above the chiral pseudo-critical temperature Tpc, and contributions from its first excitation, which gradually melt as the temperature increases. As a non-trivial test, we find that the extracted spectral function reproduces the corresponding temporal lattice correlator data for T = 220 MeV.
In quantum field theories at finite temperature spectral functions describe how particle systems behave in the presence of a thermal medium. Although data from lattice simulations can in principle be used to determine spectral function characteristics, existing methods rely on the extraction of these quantities from temporal correlators, which requires one to circumvent an illposed inverse problem. In these proceedings we report on a recent approach that instead utilises the non-perturbative constraints imposed by field locality to extract spectral function information directly from spatial correlators. In particular, we focus on the application of this approach to lattice QCD data of the spatial pseudo-scalar meson correlator in the temperature range 220−960 MeV, and outline why this data supports the conclusion that there exists a distinct pion state above the chiral pseudo-critical temperature Tpc.
The SU(3) pure gauge theory exhibits a first-order thermal deconfinement transition due to spontaneous breaking of its global Z3 center symmetry. When heavy dynamical quarks are added, this symmetry is broken explicitly and the transition weakens with decreasing quark mass until it disappears at a critical point. We compute the critical hopping parameter and the associated pion mass for lattice QCD with Nf=2 degenerate standard Wilson fermions on Nτ∈{6,8,10} lattices, corresponding to lattice spacings a=0.12 fm, a=0.09 fm, a=0.07 fm, respectively. Significant cutoff effects are observed, with the first-order region growing as the lattice gets finer. While current lattices are still too coarse for a continuum extrapolation, we estimate mcπ≈4 GeV with a remaining systematic error of ∼20%. Our results allow us to assess the accuracy of the leading-order and next-to-leading-order hopping expanded fermion determinant used in the literature for various purposes. We also provide a detailed investigation of the statistics required for this type of calculation, which is useful for similar investigations of the chiral transition.
The QCD equation of state at finite baryon density is studied in the framework of a Cluster Expansion Model (CEM), which is based on the fugacity expansion of the net baryon density. The CEM uses the two leading Fourier coefficients, obtained from lattice simulations at imaginary μB, as the only model input and permits a closed analytic form. Excellent description of the available lattice data at both μB = 0 and at imaginary μB is obtained. We also demonstrate how the Fourier coefficients can be reconstructed from baryon number susceptibilities.
Phase transitions in a non-perturbative regime can be studied by ab initio Lattice Field Theory methods. The status and future research directions for LFT investigations of Quantum Chromo-Dynamics under extreme conditions are reviewed, including properties of hadrons and of the hypothesized QCD axion as inferred from QCD topology in different phases. We discuss phase transitions in strong interactions in an extended parameter space, and the possibility of model building for Dark Matter and Electro-Weak Symmetry Breaking. Methodological challenges are addressed as well, including new developments in Artificial Intelligence geared towards the identification of different phases and transitions.
The global center symmetry of quenched QCD at zero baryonic chemical potential is broken spontaneously at a critical temperature Tc leading to a first-order phase transition. Including heavy dynamical quarks breaks the center symmetry explicitly and weakens the first-order phase transition for decreasing quark masses until it turns into a smooth crossover at a Z(2)-critical point. We investigate the Z(2)-critical quark mass value towards the continuum limit for Nf=2 flavors using lattice QCD in the staggered formulation. As part of a continued study, we present results from Monte-Carlo simulations on Nτ=8,10 lattices. Several aspect ratios and quark mass values were simulated in order to obtain the critical mass from a fit of the Polyakov loop to a kurtosis finite size scaling formula. Moreover, the possibility to develop a Ginzburg-Landau effective theory around the Z(2)-critical point is explored.
The thermodynamics of QCD with sufficiently heavy dynamical quarks can be described by a three-dimensional Polyakov loop effective theory, obtained after a truncated character and hopping expansion. We investigate the resulting phase diagram for low temperatures by mean field methods. Taking into account chemical potentials for both baryon number and isospin, we obtain clear signals for a liquid-gas type transition to baryon matter at μI=0 and a Bose-Einstein condensation transition at μB=0, as well as for their connection when both chemical potentials are non-zero.
Effective three-dimensional Polyakov loop theories derived from QCD by strong coupling and hopping expansions are valid for heavy quarks and can also be applied to finite chemical potential μ, due to their considerably milder sign problem. We apply the Monte-Carlo method to the Nf=1,2 effective theories up to O(κ4) in the hopping parameter at μ=0 to determine the critical quark mass, at which the first-order deconfinement phase transition terminates. The critical end point obtained from the effective theory to order O(κ2) agrees well with 4-dimensional QCD simulations with a hopping expanded determinant by the WHOT-QCD collaboration. We also compare with full QCD simulations and thus obtain a measure for the validity of both the strong coupling and the hopping expansion in this regime.
Approaching the continuum limit of the deconfinement critical point for Nf=2 staggered fermions
(2022)
Quenched QCD at zero baryonic chemical potential undergoes a first-order deconfinement phase transition at a critical temperature Tc, which is related to the spontaneous breaking of the global center symmetry. The center symmetry is broken explicitly by including dynamical quarks, which weaken the first-order phase transition for decreasing quark masses. At a certain critical quark mass, which corresponds to the Z(2)-critical point, the first-order phase transition turns into a smooth crossover. We investigate the Z(2)-critical quark mass for Nf=2 staggered fermions on Nτ=8,10 lattices, where larger Nτ correspond to finer lattices. Monte-Carlo simulations are performed for several quark mass values and aspect ratios in order to extrapolate to the thermodynamic limit. We present final results for Nτ=8 and preliminary results for Nτ=10 for the critical mass, which are obtained from fitting to a kurtosis finite size scaling formula of the absolute value of the Polyakov loop.
Quenched QCD at zero baryonic chemical potential undergoes a first-order deconfinement phase transition at a critical temperature Tc, which is related to the spontaneous breaking of the global center symmetry. Including heavy, dynamical quarks breaks the center symmetry explicitly and weakens the first-order phase transition. For decreasing quark masses the first-order phase transition turns into a smooth crossover at a Z2-critical point. The critical quark mass corresponding to this point has been examined with Nf=2 Wilson fermions for several Nτ in a recent study within our group. For comparison, we also locate the critical point with Nf=2 staggered fermions on Nτ=8 lattices. For this purpose we perform Monte Carlo simulations for several quark mass values and various aspect ratios in order to extrapolate to the thermodynamic limit. The critical mass is obtained by fitting to a finite size scaling formula of the kurtosis of the Polyakov loop. Our results indicate large discretization effects, requiring simulations on lattices with Nτ>8.
The so-called Columbia plot summarises the order of the QCD thermal transition as a function of the number of quark flavours and their masses. Recently, it was demonstrated that the first-order chiral transition region, as seen for Nf∈[3,6] on coarse lattices, exhibits tricritical scaling while extrapolating to zero on sufficiently fine lattices. Here we extend these studies to imaginary baryon chemical potential. A similar shrinking of the first-order region is observed with decreasing lattice spacing, which again appears compatible with a tricritical extrapolation to zero.