BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Owe Philipsen

Owe Philipsen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
hep-lathep-phnucl-thhep-thDistributed, Parallel, and Cluster Computing

Catalog footprint

What is connected

66works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

66 published item(s)

preprint2022arXiv

Deconfinement critical point of heavy quark effective lattice theories

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 $N_f=1,2$ effective theories up to $\mathcal{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 $\mathcal{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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Lattice gauge theory computation of the static force

We explore a novel approach to compute the force between a static quark and a static antiquark with lattice gauge theory directly. The approach is based on expectation values of Wilson loops or Polyakov loops with chromoelectric field insertions. We discuss theoretical and technical aspects in detail, in particular, how to compensate large discretization errors with a multiplicative renormalization factor and the evaluation using a multilevel algorithm. We also compare numerical results for the static force to corresponding results obtained in the traditional way, i.e., by computing first the static potential and then taking the derivative.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Non-perturbative insights into the spectral properties of QCD at finite temperature

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 ill-posed 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 \, \text{MeV}$, and outline why this data supports the conclusion that there exists a distinct pion state above the chiral pseudo-critical temperature $T_{\!\text{pc}}$.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Deconfinement critical point of lattice QCD with $N_{\rm f}=2$ Wilson fermions

The ${\rm SU}(3)$ pure gauge theory exhibits a first-order thermal deconfinement transition due to spontaneous breaking of its global $Z_3$ 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 $N_f=2$ degenerate standard Wilson fermions on $N_τ\in\{6,8,10\}$ lattices, corresponding to lattice spacings $a=0.12\, {\rm fm}$, $a=0.09\, {\rm fm}$, $a=0.07\, {\rm fm}$, respectively. Significant cut-off 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 $m_π^c\approx 4 {\rm GeV}$ with a remaining systematic error of $\sim 20\%$. Our results allow to assess the accuracy of the LO and NLO 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Finite size and cut-off effects on the Roberge-Weiss transition in $N_\text{f}=2$ QCD with Staggered fermions

In the absence of a genuine solution to the sign problem, lattice studies at imaginary quark chemical potential are an important tool to constrain the QCD phase diagram. We calculate the values of the tricritical quark masses in the Roberge-Weiss plane, $μ=\imathπT/3$, which separate mass regions with chiral and deconfinement phase transitions from the intermediate region, for QCD with $N_\text{f}=2$ unimproved staggered quarks on $N_τ=6$ lattices. A quantitative measure for the quality of finite size scaling plots is developed, which significantly reduces the subjective judgement required for fitting. We observe that larger aspect ratios are necessary to unambiguously determine the order of the transition than at $μ=0$. Comparing with previous results from $N_τ=4$ we find a $\sim50$% reduction in the light tricritical pion mass. The heavy tricritical pion mass stays roughly the same, but is too heavy to be resolved on $N_τ=6$ lattices and thus equally afflicted with cut-off effects. Further comparison with other discretizations suggests that current cut-off effects on the light critical masses are likely to be larger than $\sim100$%, implying a drastic shrinking of the chiral first-order region to possibly zero.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

QCD in the heavy dense regime for general $N_c$: On the existence of quarkyonic matter

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, $N_c$. 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 $N_c$ is made large. Moreover, in the baryon condensed regime, we find the pressure to scale as $p\sim N_c$ through three orders in the hopping expansion. Such a phase differs from a hadron gas with $p\sim N_c^0$, or a quark gluon plasma, $p\sim N_c^2$, 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 $N_c$ limit, which is not yet possible in practice. However, in the controlled range of lattice spacings and $N_c$-values, our results are stable when the limits are approached in this order. We discuss possible implications for physical QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A real-time lattice simulation of the thermalization of a gluon plasma: first results

To achieve an understanding of the thermalization of a quark-gluon plasma, starting from QCD without using model assumptions, is a formidable task. We study the early stage dynamics of a relativistic heavy ion collision in the framework of real time simulations of classical Yang-Mills theory in a static box with the color glass condensate as initial condition. Our study generalizes a previous one by Fukushima and Gelis from SU(2) to the realistic case of SU(3). We calculate the chromo-electric and chromo-magnetic energy densities as well as the ratio of longitudinal and transverse pressure as a function of time as probes for thermalization. Our preliminary results on coarse lattices show the occurrence of Weibel instabilities prior to thermalization.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Computing the static potential using non-string-like trial states

We present a method for computing the static quark-antiquark potential, which is not based on Wilson loops, but where the trial states are formed by eigenvector components of the covariant lattice Laplace operator. We have tested this method in SU(2) Yang-Mills theory and have obtained results with statistical errors of similar magnitude compared to a standard Wilson loop computation. The runtime of the method is, however, significantly smaller, when computing the static potential not only for on-axis, but also for many off-axis quark-antiquark separations, i.e. when a fine spatial resolution is required.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Equation of state for cold and dense heavy QCD

A previously derived three-dimensional effective lattice theory describing the thermodynamics of QCD with heavy quarks in the cold and dense region is extended through order $\sim u^5κ^8$ in the combined character and hopping expansion of the original four-dimensional Wilson action. The systematics of the effective theory is investigated to determine its range of validity in parameter space. We demonstrate the severe cut-off effects due to lattice saturation, which afflict any lattice results at finite baryon density independent of the sign problem or the quality of effective theories, and which have to be removed by continuum extrapolation. We then show how the effective theory can be solved analytically by means of a linked cluster expansion, which is completely unaffected by the sign problem, in quantitative agreement with numerical simulations. As an application, we compute the cold nuclear equation of state of heavy QCD. Our continuum extrapolated result is consistent with a polytropic equation of state for non-relativistic fermions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Roberge-Weiss transition in $N_\text{f}=2$ QCD with Wilson fermions and $N_τ=6$

QCD with imaginary chemical potential is free of the sign problem and exhibits a rich phase structure constraining the phase diagram at real chemical potential. We simulate the critical endpoint of the Roberge-Weiss (RW) transition at imaginary chemical potential for $N_\text{f}=2$ QCD on $N_τ=6$ lattices with standard Wilson fermions. As found on coarser lattices, the RW endpoint is a triple point connecting the deconfinement/chiral transitions in the heavy/light quark mass regions and changes to a second-order endpoint for intermediate masses. These regimes are separated by two tricritical values of the quark mass, which we determine by extracting the critical exponent $ν$ from a systematic finite size scaling analysis of the Binder cumulant of the imaginary part of the Polyakov loop. We are able to explain a previously observed finite size effect afflicting the scaling of the Binder cumulant in the regime of three-phase coexistence. Compared to $N_τ=4$ lattices, the tricritical masses are shifted towards smaller values. Exploratory results on $N_τ=8$ as well as comparison with staggered simulations suggest that significantly finer lattices are needed before a continuum extrapolation becomes feasible.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Roberge-Weiss transition in $N_f=2$ QCD with staggered fermions and $N_τ=6$

The QCD phase diagram at imaginary chemical potential exhibits a rich structure and studying it can constrain the phase diagram at real values of the chemical potential. Moreover, at imaginary chemical potential standard numerical techniques based on importance sampling can be applied, since no sign problem is present. In the last decade, a first understanding of the QCD phase diagram at purely imaginary chemical potential has been developed, but most of it is so far based on investigations on coarse lattices ($N_τ=4$, $a=0.3\:$fm). Considering the $N_f=2$ case, at the Roberge-Weiss critical value of the imaginary chemical potential, the chiral/deconfinement transition is first order for light/heavy quark masses and second order for intermediate values of the mass: there are then two tricritical masses, whose position strongly depends on the lattice spacing and on the discretization. On $N_τ=4$, we have the chiral $m_π^{\text{tric.}}=400\:$MeV with unimproved staggered fermions and $m_π^{\text{tric.}}\gtrsim900\:$MeV with unimproved pure Wilson fermions. Employing finite size scaling we investigate the change of this tricritical point between $N_τ=4$ and $N_τ=6$ as well as between Wilson and staggered discretizations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The $N_f=2$ QCD chiral phase transition with Wilson fermions at zero and imaginary chemical potential

The order of the thermal phase transition in the chiral limit of Quantum Chromodynamics (QCD) with two dynamical flavors of quarks is a long-standing issue and still not known in the continuum limit. Whether the transition is first or second order has important implications for the QCD phase diagram and the existence of a critical endpoint at finite densities. We follow a recently proposed approach to explicitly determine the region of first order chiral transitions at imaginary chemical potential, where it is large enough to be simulated, and extrapolate it to zero chemical potential with known critical exponents. Using unimproved Wilson fermions on coarse $N_t=4$ lattices, the first order region turns out to be so large that no extrapolation is necessary. The critical pion mass $m_π^c\approx 560$ MeV is by nearly a factor 10 larger than the corresponding one using staggered fermions. Our results are in line with investigations of three-flavour QCD using improved Wilson fermions and indicate that the systematic error on the two-flavour chiral transition is still of order 100\%.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The QCD deconfinement critical point for $N_τ=8$ with $N_f=2$ flavours of unimproved Wilson fermions

QCD at zero baryon density in the limit of infinite quark mass undergoes a first order deconfinement phase transition at a critical temperature $T_c$ corresponding to the breaking of the global centre symmetry. In the presence of dynamical quarks this symmetry is explicitly broken. Lowering the quark mass the first order phase transition weakens and terminates in a second order Z(2) point. Beyond this line confined and deconfined regions are analytically connected by a crossover transition. As the continuum limit is approached (i.e. the lattice spacing is decreased) the region of first order transitions expands towards lower masses. We study the deconfinement critical point with standard Wilson fermions and $N_f=2$ flavours. To this end we simulate several kappa values on $N_τ=8$ and various aspect ratios in order to extrapolate to the thermodynamic limit, applying finite size scaling. We estimate if and when a continuuum extrapolation is possible.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Heavy dense QCD from a 3d effective lattice theory

The cold and dense regime of the QCD phase diagram is to this day inaccessible to first principle lattice calculations owing to the sign problem. Here we present progress of an ongoing effort to probe this particularly difficult regime utilising a dimensionally reduced effective lattice theory with a significantly reduced sign problem. The effective theory is derived by combined character and hopping expansion and is valid for heavy quarks near the continuum. We show an extension of the effective theory to order $u^5κ^8$ in the cold regime. A linked cluster expansion is applied to the effective theory resulting in a consistent mechanism for handling the effective theory fully analytically. The new results are consistent with the ones from simulations confirming the viability of analytic methods. Finally we resum the analytical result which doubles the convergence region of the expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical corrections to the strong coupling effective Polyakov-line action for finite T Yang-Mills theory

We consider a three-dimensional effective theory of Polyakov lines derived previously from lattice Yang-Mills theory and QCD by means of a resummed strong coupling expansion. The effective theory is useful for investigations of the phase structure, with a sign problem mild enough to allow simulations also at finite density. In this work we present a numerical method to determine improved values for the effective couplings directly from correlators of the 4d Yang-Mills theory. For values of the gauge coupling up to the vicinity of the phase transition, the dominant short range effective coupling are well described by their corresponding strong coupling series. We provide numerical results also for the longer range interactions, Polyakov lines in higher representations as well as four-point interactions, and discuss the growing significance of non-local contributions as the lattice gets finer. Within this approach the critical Yang-Mills coupling $β_c$ is reproduced to better than one percent from a one-coupling effective theory on $N_τ=4$ lattices while up to five couplings are needed on $N_τ=8$ for the same accuracy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

QCD phase diagram from the lattice at strong coupling

The phase diagram of lattice QCD in the strong coupling limit can be measured in the full $μ$-$T$ plane, also in the chiral limit. In particular, the phase diagram in the chiral limit features a tricritical point at some $(μ_c,T_c)$. This point may be related to the critical end point expected in the QCD phase diagram. We discuss the gauge corrections to the phase diagram at strong coupling and compare our findings with various possible scenarios in continuum QCD. We comment on the possibility that the tricritical point at strong coupling is connected to the tricritical point in the continuum, massless QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The $N_f= 2$ chiral phase transition from imaginary chemical potential with Wilson Fermions

The order of the thermal transition in the chiral limit of QCD with two dynamical flavours of quarks is a long-standing issue. Still, it is not definitely known whether the transition is of first or second order in the continuum limit. Which of the two scenarios is realized has important implications for the QCD phase diagram and the existence of a critical endpoint at finite densities. Settling this issue by simulating at successively decreased pion mass was not conclusive yet. Recently, an alternative approach was proposed, extrapolating the first order phase transition found at imaginary chemical potential to zero chemical potential with known exponents, which are induced by the Roberge-Weiss symmetry. For staggered fermions on $N_t=4$ lattices, this results in a first order transition in the chiral limit. Here we report of $N_t=4$ simulations with Wilson fermions, where the first order region is found to be large.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The nature of the Roberge-Weiss Transition in $N_f=2$ QCD with Wilson Fermions on $N_τ=6$ lattices

The finite temperature chiral and deconfinement phase transitions at zero density for light and heavy quarks, respectively, have analytic continuations to imaginary chemical potential. At some critical imaginary chemical potential, they meet the Roberge-Weiss transition between adjacent $Z(3)$ sectors. For light and heavy quarks, where the chiral and deconfinement transitions are first order, the transition lines meet in a triple point. For intermediate masses chiral or deconfinement transitions are crossover and the Roberge-Weiss transition ends in a second order point. At the boundary between these regimes the junction is a tricritical point, as shown in studies with $N_f=2,3$ flavors of staggered and Wilson quarks on $N_τ=4$ lattices. Employing finite size scaling we investigate the nature of this point as a function of quark mass for $N_f=2$ flavors of Wilson fermions with a temporal lattice extent of $N_τ=6$. In particular we are interested in the change of the location of tricritical points compared to our ealier study on $N_τ=4$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Heavy dense QCD and nuclear matter from an effective lattice theory

A three-dimensional effective lattice theory of Polyakov loops is derived from QCD by expansions in the fundamental character of the gauge action, u, and the hopping parameter, κ, whose action is correct to κ^n u^m with n+m=4. At finite baryon density, the effective theory has a sign problem which meets all criteria to be simulated by complex Langevin as well as by Monte Carlo on small volumes. The theory is valid for the thermodynamics of heavy quarks, where its predictions agree with simulations of full QCD at zero and imaginary chemical potential. In its region of convergence, it is moreover amenable to perturbative calculations in the small effective couplings. In this work we study the challenging cold and dense regime. We find unambiguous evidence for the nuclear liquid gas transition once the baryon chemical potential approaches the baryon mass, and calculate the nuclear equation of state. In particular, we find a negative binding energy per nucleon causing the condensation, whose absolute value decreases exponentially as mesons get heavier. For decreasing meson mass, we observe a first order liquid gas transition with an endpoint at some finite temperature, as well as gap between the onset of isospin and baryon condensation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The chiral phase transition in two-flavor QCD from imaginary chemical potential

We investigate the order of the finite temperature chiral symmetry restoration transition for QCD with two massless fermions, by using a novel method, based on simulating imaginary values of the quark chemical potential $μ=iμ_i,μ_i\in\mathbb{R}$. Our method exploits the fact that, for low enough quark mass $m$ and large enough chemical potential $μ_i$, the chiral transition is decidedly first order, then turning into crossover at a critical mass $m_c(μ)$. It is thus possible to determine the critical line in the $m - μ^2$ plane, which can be safely extrapolated to the chiral limit by taking advantage of the known tricritical indices governing its shape. We test this method with standard staggered fermions and the result of our simulations is that $m_c(μ=0)$ is positive, so that the phase transition at zero density is definitely first order in the chiral limit, on our coarse $N_t=4$ lattices with $a\simeq 0.3\,\mathrm{fm}$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The lattice QCD phase diagram in and away from the strong coupling limit

We study lattice QCD with four flavors of staggered quarks. In the limit of infinite gauge coupling, "dual" variables can be introduced, which render the finite-density sign problem mild and allow a full determination of the $μ-T$ phase diagram by Monte Carlo simulations, also in the chiral limit. However, the continuum limit coincides with the weak coupling limit. We propose a strong-coupling expansion approach towards the continuum limit. We show first results, including the phase diagram and its chiral critical point, from this expansion truncated to next-to-leading order.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The nature of the Roberge-Weiss transition in $N_f=2$ QCD with Wilson fermions

At imaginary values of the quark chemical potential $μ$, Quantum Chromodynamics shows an interesting phase structure due to an exact center, or Roberge-Weiss (RW), symmetry. This can be used to constrain QCD at real $μ$, where the sign problem prevents Monte Carlo simulations of the lattice theory. In previous studies of this region with staggered fermions it was found that the RW endpoint, where the center transition changes from first-order to a crossover, depends non-trivially on the quark mass: for high and low masses, it is a triple point connecting to the deconfinement and chiral transitions, respectively, changing to a second-order endpoint for intermediate mass values. These parameter regions are separated by tricritical points. Here we present a confirmation of these findings using Wilson fermions on $N_τ=4$ lattices. In addition, our results provide a successful quantitative check for a heavy quark effective lattice theory at finite density.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Phase Diagram of Strong Coupling QCD including Gauge Corrections

The strong coupling limit of lattice QCD with staggered fermions has been studied for decades, both via Monte Carlo and via mean field theory. In this model, the finite density sign problem can be made mild and the full phase diagram can be obtained, even in the chiral limit. It is however desirable to understand the effect of a finite lattice gauge coupling $β$ on the phase diagram in the $μ-T$ plane in order to understand how it evolves into the phase diagram of continuum QCD. Here we discuss how to construct a partition function for non-zero lattice coupling, exact to $\mathcal{O}(β)$, and present corresponding Monte Carlo results, in particular for corrections to the chiral susceptibility and to the phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The second order hydrodynamic transport coefficient $κ$ for the gluon plasma from the lattice

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 $120^3\times 6,8$ lattices with $a<0.1$ fm. In the temperature range $2T_\mathrm{c}-10T_\mathrm{c}$ we find $κ=0.36(15)T^2$. 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effective lattice Polyakov loop theory for finite temperature Yang-Mills

Effective Polyakov loop theories are a useful tool for an investigation of pure Yang-Mills theory and full QCD. A systematic derivation of the effective action can be done in a spatial strong coupling expansion. Quite accurate predictions for the deconfinement phase transition of Yang-Mills theory have been obtained in this approach. Besides the critical couplings, further observables can be measured in the effective theory. These provide additional tests for the reliability of the strong coupling approach and the truncation of the effective action. In this contribution we will present recent results for the free energy of the static quark-antiquark pair and the equation of state.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effective lattice Polyakov loop theory vs. full SU(3) Yang-Mills at finite temperature

A three-dimensional effective theory of Polyakov loops has recently been derived from Wilson's Yang-Mills lattice action by means of a strong coupling expansion. It is valid in the confined phase up to the deconfinement phase transition, for which it predicts the correct order and gives quantitative estimates for the critical coupling. In this work we study its predictive power for further observables like correlation functions and the equation of state. We find that the effective theory correctly reproduces qualitative features and symmetries of the full theory as the continuum is approached. Regarding quantitative predictions, we identify two classes of observables by numerical comparison as well as analytic calculations: correlation functions and their associated mass scales cannot be described accurately from a truncated effective theory, due to its inherently non-local nature involving long-range couplings. On the other hand, phase transitions and bulk thermodynamic quantities are accurately reproduced by the leading local part of the effective theory. In particular, the effective theory description is numerically superior when computing the equation of state at low temperatures or the properties of the phase transition.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the definition and interpretation of a static quark anti-quark potential in the colour-adjoint channel

We study possibilities to define a static quark anti-quark pair in a colour-adjoint orientation based on Wilson loops with generator insertions, using both lattice QCD and leading order perturbation theory in various gauges. Non-perturbatively, the only way to obtain non-zero results while maintaining positivity of the Hamiltonian is by some form of temporal gauge. In this case the correlator is equivalent to a gauge invariant correlation function of a static quark anti-quark pair and a static adjoint quark, the resulting three-point potential is attractive. Saturating open colour indices with colour magnetic fields instead also leads to a gauge invariant correlator. However this object is found to couple to the singlet sector only. There appears to be no lattice observable that reproduces the repulsive adjoint potential predicted by perturbation theory in Lorenz or Coulomb gauges.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Onset Transition to Cold Nuclear Matter from Lattice QCD with Heavy Quarks

Lattice QCD at finite density suffers from a severe sign problem, which has so far prohibited simulations of the cold and dense regime. Here we study the onset of nuclear matter employing a three-dimensional effective theory derived by combined strong coupling and hopping expansions, which is valid for heavy but dynamical quarks and has a mild sign problem only. Its numerical evaluations agree between a standard Metropolis and complex Langevin algorithm, where the latter is free of the sign problem. Our continuum extrapolated data clearly show a first order phase transition building up at $μ_B \approx m_B$ as the temperature approaches zero. An excellent description of the data is achieved by an analytic solution in the strong coupling limit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Onset Transition to Cold Nuclear Matter from Lattice QCD with Heavy Quarks to $κ^4$

We present results of our ongoing studies of an effective three-dimensional theory of thermal lattice QCD with heavy Wilson quarks. This is done by combined strong coupling and hopping parameter expansions. The full quark determinant of four dimensional lattice QCD is expanded in orders of the hopping parameter $κ$, the dimensional reduction is achieved by integrating over the spatial links. We present the calculation of the effective theory through order $κ^nu^m$ with $n+m=4$. This theory is then used to simulate heavy quarks near the cold and dense limit. For nonzero chemical potential the theory suffers from a sign problem, wich is avoided by employing stochastical quantisation. Continuum extrapolated results for the onset of nuclear matter are shown and the region of convergence of the effective theory is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

QCD thermodynamics with O(a) improved Wilson fermions at Nf=2

We present an update of our study of the phase diagram of two-flavour QCD at zero baryon density with dynamical $O(a)$ improved Wilson quarks. All simulations are done on lattices with a temporal extent of $N_t=16$ and spatial extent $L=32,48$ and 64, ensuring that discretisation effects are small and finite size effects can be controlled. In the approach to the chiral limit we currently have three scans with pion masses between 540 and 200 MeV. In this proceedings article the focus is on the new scan at $m_π=200$ MeV and the measurement of screening masses. We also present first results concerning a test of scaling in the approach to the chiral limit and the chiral extrapolation of the difference of screening masses in scalar and pseudoscalar channels, which provides a measure for the strength of the anomalous breaking of the $U_A(1)$ symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The chiral phase transition for two-flavour QCD at imaginary and zero chemical potential

The chiral symmetry of QCD with two massless quark flavours gets restored in a non-analytic chiral phase transition at finite temperature and zero density. Whether this is a first-order or a second-order transition has not yet been determined unambiguously, due to the difficulties of simulating light quarks. We investigate the nature of the chiral transition as a function of quark mass and imaginary chemical potential, using staggered fermions on N_t=4 lattices. At sufficiently large imaginary chemical potential, a clear signal for a first-order transition is obtained for small masses, which weakens with decreasing imaginary chemical potential. The second-order critical line m_c(mu_i), which marks the boundary between first-order and crossover behaviour, extrapolates to a finite m_c(mu_i=0) with known critical exponents. This implies a definitely first-order transition in the chiral limit on relatively coarse, N_t=4 lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The QCD equation of state from the lattice

The equation of state of QCD at finite temperatures and baryon densities has a wide range of applications in many fields of modern particle and nuclear physics. It is the main ingredient to describe the dynamics of experimental heavy ion collisions, the expansion of the early universe in the standard model era and the interior of compact stars. On most scales of interest, QCD is strongly coupled and not amenable to perturbative investigations. Over the last decade, first principles calculations using lattice QCD have reached maturity, in the sense that for particular discretisation schemes simulations at the physical point have become possible, finite temperature results near the continuum limit are available and systematic errors begin to be controlled. This review summarises the current theoretical and numerical state of the art based on staggered and Wilson fermions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraints on the two-flavor QCD phase diagram from imaginary chemical potential

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 mu=i pi/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 two-flavour chiral limit.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Lattice QCD based on OpenCL

We present an OpenCL-based Lattice QCD application using a heatbath algorithm for the pure gauge case and Wilson fermions in the twisted mass formulation. The implementation is platform independent and can be used on AMD or NVIDIA GPUs, as well as on classical CPUs. On the AMD Radeon HD 5870 our double precision dslash implementation performs at 60 GFLOPS over a wide range of lattice sizes. The hybrid Monte-Carlo presented reaches a speedup of four over the reference code running on a server CPU.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Phase transitions in heavy-quark QCD from an effective theory

With combined hopping parameter and strong coupling expansions, we calculate a dimensionally reduced Polyakov-loop effective theory valid for heavy quarks at nonzero temperature and arbitrary chemical potential. We numerically compute the critical endpoint of the deconfinement transition as a function of quark masses and number of flavours. We also investigate the applicability of the model to the low-T and high density region, specifically in terms of baryon condensation phenomena.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pseudo-Critical Temperature and Thermal Equation of State from $N_f=2$ Twisted Mass Lattice QCD

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 $T_c$ for three values of the pion mass in the range of $m_\mathrm{PS} \simeq 300$ and $500 \mathrm{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 $N_τ=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 $m_\mathrm{PS} \simeq 400$ and $700 \mathrm{MeV}$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

QCD thermodynamics with two flavours of Wilson fermions on large lattices

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 N_t=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 κ, 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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Centre symmetric 3d effective actions for thermal SU(N) Yang-Mills from strong coupling series

We derive three-dimensional, Z(N)-symmetric effective actions in terms of Polyakov loops by means of strong coupling expansions, starting from thermal SU(N) Yang-Mills theory in four dimensions on the lattice. An earlier action in the literature, corresponding to the (spatial) strong coupling limit, is thus extended by several higher orders, as well as by additional interaction terms. We provide analytic mappings between the couplings of the effective theory and the parameters $N_τ,β$ of the original thermal lattice theory, which can be systematically improved. We then investigate the deconfinement transition for the cases SU(2) and SU(3) by means of Monte Carlo simulations of the effective theory. Our effective models correctly reproduce second order 3d Ising and first order phase transitions, respectively. Furthermore, we calculate the critical couplings $β_c(N_τ)$ and find agreement with results from simulations of the 4d theory at the few percent level for $N_τ=4-16$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Status of the QCD phase diagram from lattice calculations

The present knowledge of the QCD phase diagram based on simulations of lattice QCD is summarised. The main questions are whether there is a critical point in the QCD phase diagram and whether it is related to a chiral phase transition. It is shown that QCD at imaginary chemical potentials has a rich phase structure, which can be determined in a controlled way without sign problem and which severely constrains the phase structure at real chemical potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Strong coupling effective theory with heavy fermions

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, mu/T, T) space over the whole range of chemical potentials from zero up to infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Strong-coupling effective action(s) for SU(3) Yang-Mills

We apply strong-coupling expansion techniques to finite-temperature lattice pure gauge theory, obtaining dimensionally reduced $Z_N$-symmetric effective theories. The analytic mappings between the effective couplings and the original one, viz. $β$, allow to estimate the transition point $β_c$ of the 4D theory for a large range of the imaginary-time extent $N_τ$ of the lattice. We study the models for SU(3) via Monte Carlo simulation, finding satisfactory agreement with the critical point of the original theories especially at low $N_τ$. We have fixed an error in the group measure used in arXiv:1010.0951 and provide here the correct numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The QCD deconfinement transition for heavy quarks and all baryon chemical potentials

Using combined strong coupling and hopping parameter expansions, we derive an effective three-dimensional theory from thermal lattice QCD with heavy Wilson quarks. The theory depends on traced Polyakov loops only and correctly reflects the centre symmetry of the pure gauge sector as well as its breaking by finite mass quarks. It is valid up to certain orders in the lattice gauge coupling and hopping parameter, which can be systematically improved. To its current order it is controlled for lattices up to N_τ\sim 6 at finite temperature. For nonzero quark chemical potentials, the effective theory has a fermionic sign problem which is mild enough to carry out simulations up to large chemical potentials. Moreover, by going to a flux representation of the partition function, the sign problem can be solved. As an application, we determine the deconfinement transition and its critical end point as a function of quark mass and all chemical potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Towards corrections to the strong coupling limit of staggered lattice QCD

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(β)$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraining the QCD phase diagram by tricritical lines at imaginary chemical potential

We present unambiguous evidence from lattice simulations of QCD with three degenerate quark species for two tricritical points in the (T,m) phase diagram at fixed imaginary μ/T=iπ/3 mod 2π/3, one in the light and one in the heavy mass regime. These represent the boundaries of the chiral and deconfinement critical lines continued to imaginary chemical potential, respectively. It is demonstrated that the shape of the deconfinement critical line for real chemical potentials is dictated by tricritical scaling and implies the weakening of the deconfinement transition with real chemical potential. The generalization to non-degenerate and light quark masses is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraints for the QCD phase diagram from imaginary chemical potential

We present unambiguous evidence from lattice simulations of N_f=3 QCD for two tricritical points in the (T,m) phase diagram at fixed imaginary μ/T=iπ/3 mod. 2π/3, one in the light and one in the heavy quark regime. Together with similar results in the literature for N_f=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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cutoff effects of Wilson fermions on the QCD equation of state to O(g^2)

We compute the O(g^2) contribution to the thermodynamic pressure for Wilson fermions in the standard, the twisted mass, and clover improved formulation in lattice perturbation theory, including finite mass effects. We compare the continuum approaches of these discretizations for the massive ideal and interacting gas. In all cases, for N_t \geq 8 cutoff effects of Wilson type fermions are comparable to those of staggered fermions, but asymptotic scaling requires $N_τ>10$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Effective Polyakov-loop theory for pure Yang-Mills from strong coupling expansion

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 $N_τ$ at once.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The deconfinement transition of finite density QCD with heavy quarks from strong coupling series

Starting from Wilson's action, we calculate strong coupling series for the Polyakov loop susceptibility in lattice gauge theories for various small N_τin the thermodynamic limit. Analysing the series with Padé approximants, we estimate critical couplings and exponents for the deconfinement phase transition. For SU(2) pure gauge theory our results agree with those from Monte-Carlo simulations within errors, which for the coarser N_τ=1,2 lattices are at the percent level. For QCD we include dynamical fermions via a hopping parameter expansion. On a N_τ=1 lattice with N_f=1,2,3, we locate the second order critical point where the deconfinement transition turns into a crossover. We furthermore determine the behaviour of the critical parameters with finite chemical potential and find the first order region to shrink with growing μ. Our series moreover correctly reflects the known Z(N) transition at imaginary chemical potential.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The pressure of strong coupling lattice QCD with heavy quarks, the hadron resonance gas model and the large N limit

In this paper we calculate the pressure of pure lattice Yang-Mills theories and lattice QCD with heavy quarks by means of strong coupling expansions. Dynamical fermions are introduced with a hopping parameter expansion, which also allows for the incorporation of finite quark chemical potential. We show that in leading orders the results are in full agreement with expectations from the hadron resonance gas model, thus validating it with a first principles calculation. For pure Yang-Mills theories we obtain the corresponding ideal glueball gas, in QCD with heavy quarks our result equals that of an ideal gas of mesons and baryons. Another finding is that the Yang-Mills pressure in the large N limit is of order $\sim N^0$ to the calculated orders, when the inverse 't Hooft coupling is used as expansion parameter. This property is expected in the confined phase, where our calculations take place.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Towards the N_f=2 deconfinement transition temperature with O(a) improved Wilson fermions

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 $\Ord(a)$ improved Wilson fermions. We start with $N_{t}=12$ and 16 lattices and pion masses of 600 to 450 MeV, aiming at chiral and continuum limits with light quarks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Towards the N_f=2 deconfinement transition temperature with O(a) improved Wilson fermions: An update

We give an update on our current project to determine the transition temperature and the order of the deconfinement transition in the chiral limit of two flavour QCD. We use nonperturbatively O(a) improved Wilson fermions of the Sheikholeslami-Wohlert type, employing the efficient deflation accelerated DDHMC algorithm. We start at lattices with N_t>=12 and pion masses below 600 MeV, aiming at chiral and continuum limits with light quarks.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Heavy Quark Thermalization in Classical Lattice Gauge Theory: Lessons for Strongly-Coupled QCD

Thermalization of a heavy quark near rest is controlled by the correlator of two electric fields along a temporal Wilson line. We address this correlator within real-time, classical lattice Yang-Mills theory, and elaborate on the analogies that exist with the dynamics of hot QCD. In the weak-coupling limit, it can be shown analytically that the dynamics on the two sides are closely related to each other. For intermediate couplings, we carry out non-perturbative simulations within the classical theory, showing that the leading term in the weak-coupling expansion significantly underestimates the heavy quark thermalization rate. Our analytic and numerical results also yield a general understanding concerning the overall shape of the spectral function corresponding to the electric field correlator, which may be helpful in subsequent efforts to reconstruct it from Euclidean lattice Monte Carlo simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Lattice QCD Thermodynamics on the Grid

We describe how we have used simultaneously ${\cal O}(10^3)$ nodes of the EGEE Grid, accumulating ca. 300 CPU-years in 2-3 months, to determine an important property of Quantum Chromodynamics. We explain how Grid resources were exploited efficiently and with ease, using user-level overlay based on Ganga and DIANE tools above standard Grid software stack. Application-specific scheduling and resource selection based on simple but powerful heuristics allowed to improve efficiency of the processing to obtain desired scientific results by a specified deadline. This is also a demonstration of combined use of supercomputers, to calculate the initial state of the QCD system, and Grids, to perform the subsequent massively distributed simulations. The QCD simulation was performed on a $16^3\times 4$ lattice. Keeping the strange quark mass at its physical value, we reduced the masses of the up and down quarks until, under an increase of temperature, the system underwent a second-order phase transition to a quark-gluon plasma. Then we measured the response of this system to an increase in the quark density. We find that the transition is smoothened rather than sharpened. If confirmed on a finer lattice, this finding makes it unlikely for ongoing experimental searches to find a QCD critical point at small chemical potential.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Screened perturbation theory for 3d Yang-Mills theory and the magnetic modes of hot QCD

Perturbation theory for non-abelian gauge theories at finite temperature is plagued by infrared divergences which are caused by magnetic soft modes ~g^2T, 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 ~g^6 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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Towards a determination of the chiral critical surface of QCD

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 {m_{u,d}, m_s,μ} parameter space. The potential critical endpoint of the QCD (T,μ)-phase diagram is widely expected to be part of this surface. Since for μ=0 with physical quark masses QCD is known to exhibit an analytic crossover, this expectation requires the region of chiral transitions to expand with μfor a chiral critical endpoint to exist. Instead, on coarse N_t=4 lattices, we find the area of chiral transitions to shrink with μ, which excludes a chiral critical point for QCD at moderate chemical potentials μ_B < 500 MeV. First results on finer N_t=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 N_f=2 and N_f=2+1 theories and its consequences for physical QCD.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Strong coupling expansion for finite temperature Yang-Mills theory in the confined phase

We perform euclidean strong coupling expansions for Yang Mills theory on the lattice at finite temperature. After setting up the formalism for general SU(N), we compute the first few terms of the series for the free energy density and the lowest screening mass in the case of SU(2). To next-to-leading order the free energy series agrees with that of an ideal gas of glueballs. This demonstrates that in the confined phase the quasi-particles indeed correspond to the T=0 hadron excitations, as commonly assumed in hadron resonance gas models. Our result also fixes the lower integration constant for Monte Carlo calculations of the thermodynamic pressure via the integral method. In accord with Monte Carlo results, we find screening masses to be nearly temperature independent in the confined phase. This and the exponential smallness of the pressure can be understood as genuine strong coupling effects. Finally, we analyse Pade approximants to estimate the critical couplings of the phase transition, which for our short series are only ~25% accurate. However, up to these couplings the equation of state agrees quantitatively with numerical results on N_t=1-4 lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The curvature of the critical surface (m_ud,m_s)^{crit}(mu): a progress report

At zero chemical potential mu, the order of the temperature-driven quark-hadron transition depends on the quark masses m_{u,d} and m_s. Along a critical line bounding the region of first-order chiral transitions in the (m_{u,d},m_s) plane, this transition is second order. When the chemical potential is turned on, this critical line spans a surface, whose curvature at mu=0 can be determined without any sign or overlap problem. Our past measurements on N_t=4 lattices suggest that the region of quark masses for which the transition is first order shrinks when mu is turned on, which makes a QCD chiral critical point at small mu/T unlikely. We present results from two complementary methods, which can be combined to yield information on higher-order terms. It turns out that the O(mu^4) term reinforces the effect of the leading O(mu^2) term, and there is strong evidence that the O(mu^6) and O(mu^8) terms do as well. We also report on simulations underway, where the strange quark is given its physical mass, and where the lattice spacing is reduced.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The chiral critical line of N_f=2+1 QCD at zero and non-zero baryon density

We present numerical results for the location of the chiral critical line at finite temperature and zero and non-zero baryon density for QCD with N_f=2+1 flavours of staggered fermions on lattices with temporal extent N_t=4. For degenerate quark masses, we compare our results obtained with the exact RHMC algorithm with earlier, inexact R-algorithm results and find a reduction of 25% in the critical quark mass, for which the first order phase transition changes to a smooth crossover. Extending our analysis to non-degenerate quark masses, we map out the chiral critical line up to the neighbourhood of the physical point, which we confirm to be in the crossover region. Our data are consistent with a tricritical point at a strange quark mass of ~500 MeV. Finally, we investigate the shift of the critical line with finite baryon density, by simulating with an imaginary chemical potential for which there is no sign problem. We observe this shift to be very small or, conversely, the critical endpoint μ^c(m_{u,d},m_s) to be extremely quark mass sensitive. Moreover, the sign of this shift is opposite to standard expectations. If confirmed on a finer lattice, it implies the absence of a critical endpoint or phase transition for chemical potentials μ_B < 500 MeV. We thus argue that finer lattices are required to settle even the qualitative features of the QCD phase diagram.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Cooling, Physical Scales and Topology

We develop a cooling method controlled by a physical cooling radius that defines a scale below which fluctuations are smoothed out while leaving physics unchanged at all larger scales. We apply this method to study topological properties of lattice gauge theories, in particular the behaviour of instantons, dislocations and instanton--anti-instanton pairs. Monte Carlo results for the SU(2) topology are presented. We find that the method provides a means to prevent instanton--anti-instanton annihilation under cooling. While the instanton sizes are largely independent from the smoothing scale, the density and pair separations are determined by the particular choice made for this quantity. We discuss the questions this raises for the "physicality" of these concepts.

0 citations0 reviewsNo code linkedOpen access