Filtern
Erscheinungsjahr
Dokumenttyp
- Konferenzveröffentlichung (35) (entfernen)
Sprache
- Englisch (35)
Volltext vorhanden
- ja (35)
Gehört zur Bibliographie
- nein (35) (entfernen)
Schlagworte
Institut
- Physik (35) (entfernen)
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).
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.
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.
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 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).
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
LatticeQCD using OpenCL
(2011)
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.
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.
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.
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.
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 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.
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.
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.
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.
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.