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

Philippe De Forcrand

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

ResearcherUnclaimed source record
hep-lathep-phcond-mat.stat-mechhep-thnucl-thcond-mat.str-elcond-matcond-mat.otherDistributed, Parallel, and Cluster Computinggr-qcphysics.comp-phq-bio

Catalog footprint

What is connected

59works
12topics
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

59 published item(s)

preprint2022arXiv

Diquark properties from full QCD lattice simulations

We study diquarks on the lattice in the background of a static quark, in a gauge-invariant formalism with quark masses down to almost physical $m_π$. We determine mass differences between diquark channels as well as diquark-quark mass differences. The lightest and next-to-lightest diquarks have "good" scalar, $\bar{3}_F$, $\bar{3}_c$, $J^P=0^+$, and "bad" axial vector, $6_F$, $\bar{3}_c$, $J^P=1^+$, quantum numbers, and a bad-good mass difference for $ud$ flavors, $198(4)~\rm{MeV}$, in excellent agreement with phenomenological determinations. Quark-quark attraction is found only in the "good" diquark channel. We extract a corresponding diquark size of $\sim 0.6~\rm{fm}$ and perform a first exploration of the "good" diquark shape, which is shown to be spherical. Our results provide quantitative support for modeling the low-lying baryon spectrum using good light diquark effective degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Good and bad diquark properties and spatial correlations in lattice QCD

We study good, bad and not-even-bad diquarks on the lattice in a gauge-invariant formalism in full QCD. We establish their spectral masses with short extrapolations to the physical point, observing agreement with phenomenological expectations. We find that only the good diquark has attractive quark-quark spatial correlations, with spherical shape and size $\sim0.6$ fm. Our results provide quantitative support for modelling the low-lying baryon spectrum using good light diquark effective degrees of freedom.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Aspects of topological actions on the lattice

We consider a lattice action which forbids large fields, and which remains invariant under smooth deformations of the field. Such a "topological" action depends on one parameter, the field cutoff, but does not have a classical continuum limit as this cutoff approaches zero. We study the properties of such an action in 4d compact U(1) lattice gauge theory, and compare them with those of the Wilson action. In both cases, we find a weakly first-order transition separating a confining phase where monopoles condense, and a Coulomb phase where monopoles are exponentially suppressed. We also find a different, critical value of the field cutoff where monopoles completely disappear. Finally, we show that a topological action simplifies the measurement of the free energy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Higgs-Yukawa model with higher dimension operators via extended mean field theory

Using extended mean field theory (EMFT) on the lattice, we study properties of the Higgs-Yukawa model as an approximation of the standard model Higgs sector, and the effect of higher dimension operators. We note that the discussion of vacuum stability is completely modified in the presence of a $ϕ^6$ term, and that the Higgs mass no longer appears fine tuned. We also study the finite temperature transition. Without higher dimension operators the transition is found to be second order (crossover with gauge fields) for the experimental value of the Higgs mass $M_h=125$ GeV. By taking a $ϕ^6$ interaction in the Higgs potential as a proxy for a UV completion of the standard model, the transition becomes stronger and turns first order if the scale of new physics, i.e. the mass of the lightest mediator particle, is around $1.5$ TeV. This implies that electroweak baryogenesis may be viable in models which introduce new particles around that scale.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Oscillating propagators in heavy-dense QCD

Using Monte Carlo simulations and extended mean field theory calculations we show that the $3$-dimensional $Z_3$ spin model with complex external fields has non-monotonic spatial correlators in some regions of its parameter space. This model serves as a proxy for heavy-dense QCD in $(3+1)$ dimensions. Non-monotonic spatial correlators are intrinsically related to a complex mass spectrum and a liquid-like (or crystalline) behavior. A liquid phase could have implications for heavy-ion experiments, where it could leave detectable signals in the spatial correlations of baryons.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sampling of General Correlators in Worm Algorithm-based Simulations

Using the complex $ϕ^4$-model as a prototype for a system which is simulated by a worm algorithm, we show that not only the charged correlator $<ϕ^{*}(x)ϕ(y)>$, but also more general correlators such as $<|ϕ(x)||ϕ(y)|>$ or $<\text{arg}(ϕ(x))\text{arg}(ϕ(y))>$, as well as condensates like $<|ϕ|>$, can be measured at every step of the Monte Carlo evolution of the worm instead of on closed-worm configurations only. The method generalizes straightforwardly to other systems simulated by worms, such as spin or sigma models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Slab Method to Measure the Topological Susceptibility

In simulations of a model with topological sectors, algorithms which proceed in small update steps tend to get stuck in one sector, especially on fine lattices. This distorts the numerical results; in particular it is not straightforward to measure the topological susceptibility chi_t. We test a method to measure chi_t even if configurations from only one sector are available. It is based on the topological charges in sub-volumes, which we denote as "slab". This enables the evaluation of chi_t, as we demonstrate with numerical results for non-linear sigma-models and for 2-flavour QCD. In the latter case, the gradient flow is applied for the smoothing of the gauge configurations, and the slab method results for chi_t are stable over a broad range of flow times.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Two-Flavor Lattice QCD with a Finite Density of Heavy Quarks: Heavy-Dense Limit and "Particle-Hole" Symmetry

We investigate the properties of the half-filling point in lattice QCD (LQCD), in particular the disappearance of the sign problem and the emergence of an apparent particle-hole symmetry, and try to understand where these properties come from by studying the heavy-dense fermion determinant and the corresponding strong-coupling partition function (which can be integrated analytically). We then add in a first step an effective Polyakov loop gauge action (which reproduces the leading terms in the character expansion of the Wilson gauge action) to the heavy-dense partition function and try to analyze how some of the properties of the half-filling point change when leaving the strong coupling limit. In a second step, we take also the leading nearest-neighbor fermion hopping terms into account (including gauge interactions in the fundamental representation) and mention how the method could be improved further to incorporate the full set of nearest-neighbor fermion hoppings. Using our mean-field method, we also obtain an approximate ($μ$,T) phase diagram for heavy-dense LQCD at finite inverse gauge coupling $β$. Finally, we propose a simple criterion to identify the chemical potential beyond which lattice artifacts become dominant.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Condensation phenomena in two-flavor scalar QED at finite chemical potential

We study condensation in two-flavored, scalar QED with non-degenerate masses at finite chemical potential. The conventional formulation of the theory has a sign problem at finite density which can be solved using an exact reformulation of the theory in terms of dual variables. We perform a Monte Carlo simulation in the dual representation and observe a condensation at a critical chemical potential $μ_c$. After determining the low-energy spectrum of the theory we try to establish a connection between $μ_c$ and the mass of the lightest excitation of the system, which are naively expected to be equal. It turns out, however, that the relation of the critical chemical potential to the mass spectrum in this case is non-trivial: Taking into account the form of the condensate and making some simplifying assumptions we suggest an adequate explanation which is supported by numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effects of higher dimension operators on the Standard Model Higgs sector

We study the effect of higher dimension operators on the electroweak finite temperature phase transition in two sectors of the Standard Model. Firstly, the Higgs-Yukawa sector, consisting of the Higgs doublet and the massive Standard Model fermions, is studied with an approximate method, Extended Mean Field Theory. Secondly, the gauge-Higgs sector, consisting of the Higgs doublet and the gauge fields of the weak interaction, is studied using Monte Carlo simulations. In both cases we find that a cutoff scale of around 1.5 TeV is needed to make the electroweak phase transition first order at the experimental value of the Higgs boson mass, which is a requirement for making electroweak baryogenesis viable.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gauge-invariant signatures of spontaneous gauge symmetry breaking by the Hosotani mechanism

The Hosotani mechanism claims to achieve gauge-symmetry breaking, for instance $SU(3) \to SU(2)\times U(1)$. To verify this claim, we propose to monitor the stability of a topological defect stable under a gauge subgroup but not under the whole gauge group, like a $U(1)$ flux state or monopole in the case above. We use gauge invariant operators to probe the presence of the topological defect to avoid any ambiguity introduced by gauge fixing. Our method also applies to an ordinary gauge-Higgs system.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Integrating out lattice gauge fields

The sign problem is a major obstacle to our understanding of the phase diagram of QCD at finite baryon density. Several numerical methods have been proposed to tackle this problem, but a full solution to the sign problem is still elusive. Motivated by this problem and by recent advances in diagrammatic Monte Carlo methods, we find a new exact representation of the partition function of pure lattice gauge theory that contains no link variables. This approach can be easily extended to include staggered fermions, and results in a diagrammatic representation of fermionic states as arrangements of monomers, dimers, and fermionic loops saturating the spacetime lattice. Our representations are exact for any value of the lattice coupling, and extend previous representations that are only valid in the strong coupling limit and at $O(β)$. As a concrete example, we construct a monomer-dimer-loop representation of compact lattice QED.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lattice simulation of the SU(2) chiral model at zero and non-zero pion density

We propose a flux representation based lattice formulation of the partition function corresponding to the SU(2) principal chiral Lagrangian, including a chemical potential and scalar/pseudo-scalar source terms. Lattice simulations are then used to obtain non-perturbative properties of the theory, in particular its mass spectrum at zero and non-zero pion density. We also sketch a method to efficiently measure general one- and two-point functions during the worm updates.

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

Topological Susceptibility from Slabs

In quantum field theories with topological sectors, a non-perturbative quantity of interest is the topological susceptibility chi_t. In principle it seems straightforward to measure chi_t by means of Monte Carlo simulations. However, for local update algorithms and fine lattice spacings, this tends to be difficult, since the Monte Carlo history rarely changes the topological sector. Here we test a method to measure chi_t even if data from only one sector are available. It is based on the topological charges in sub-volumes, which we denote as slabs. Assuming a Gaussian distribution of these charges, this method enables the evaluation of chi_t, as we demonstrate with numerical results for non-linear sigma-models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

U(1) lattice gauge theory with a topological action

We investigate the phase diagram of the compact $U(1)$ lattice gauge theory in four dimensions using a non-standard action which is invariant under continuous deformations of the plaquette angles. Just as for the Wilson action, we find a weakly first order transition, separating a confining phase where magnetic monopoles condense, and a Coulomb phase where monopoles are dilute. We also find a third phase where monopoles are completely absent. The topological action offers an algorithmic advantage for the computation of the free energy.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Deformations of infrared-conformal theories in two dimensions

We study two exactly solvable two-dimensional conformal models, the critical Ising model and the Sommerfield model, on the lattice. We show that finite-size effects are important and depend on the aspect ratio of the lattice. In particular, we demonstrate how to obtain the correct massless behavior from an infinite tower of finite-size-induced masses and show that it is necessary to first take the cylindrical geometry limit in order to get correct results. In the Sommerfield model we also introduce a mass deformation to measure the mass anomalous dimension, $γ_m$. We find that the explicit scale breaking of the lattice setup induces corrections which must be taken into account in order to reproduce $γ_m$ at the infrared fixed point. These results can be used to improve the methodology in the search for the conformal window in QCD-like theories with many flavors.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Extended Mean Field study of complex $ϕ^4$-theory at finite density and temperature

We review the Extended Mean Field Theory (EMFT) approximation and apply it to complex, scalar $ϕ^4$-theory on the lattice. We study the critical properties of the Bose condensation driven by a nonzero chemical potential $μ$ at both zero and nonzero temperature and determine the $(T,μ)$ phase diagram. The results are in very good agreement with recent Monte Carlo data for all parameter values considered. EMFT can be formulated directly in the thermodynamic limit which allows us to study lattice spacings for which Monte Carlo studies are not feasible with present techniques. We find that the EMFT approximation accurately reproduces many known phenomena of the exact solution, like the "Silver Blaze" behavior at zero temperature and dimensional reduction at finite temperature.

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 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
preprint2013arXiv

Complex $φ^4$ Theory at Finite Temperature and Density via Extended Mean Field Theory

We review the Extended Mean Field (EMFT) approximation and apply it to complex, scalar $φ^4$ theory on the lattice. We determine the ($T,μ$) phase diagram and study the critical properties of the Bose condensation transition, at zero and finite temperature. In both cases we obtain results which agree very well with recent Monte Carlo studies. Within our approximation we can reach the thermodynamic limit and can thus obtain results at lattice spacings unreachable by present Monte Carlo simulations. We find that our approximation reproduces accurately many phenomena of the model, like the "Silver Blaze" behavior at zero temperature and dimensional reduction at finite temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Conformality in many-flavour lattice QCD at strong coupling

It is widely believed that chiral symmetry is spontaneously broken at zero temperature in the strong coupling limit of staggered fermions, for any number of colors and flavors. Using Monte Carlo simulations, we show that this conventional wisdom, based on a mean-field analysis, is wrong. For sufficiently many fundamental flavors, chiral symmetry is restored via a bulk, first-order transition. This chirally symmetric phase appears to be analytically connected with the expected conformal window of many-flavor continuum QCD. We perform simulations in the chirally symmetric phase at zero quark mass for various system sizes L, and measure the torelon mass, the Dirac spectrum and the hadron spectrum. All masses go to zero with 1/L. L is hence the only infrared length scale. Thus, the strong-coupling chirally restored phase appears as a convenient laboratory to study IR-conformality. Finally, we present a conjecture for the phase diagram of lattice QCD as a function of the bare coupling and the number of quark flavors.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical Mean Field Approximation Applied to Quantum Field Theory

We apply the Dynamical Mean Field (DMFT) approximation to the real, scalar phi^4 quantum field theory. By comparing to lattice Monte Carlo calculations, perturbation theory and standard mean field theory, we test the quality of the approximation in two, three, four and five dimensions. The quantities considered in these tests are the critical coupling for the transition to the ordered phase and the associated critical exponents nu and beta. We also map out the phase diagram in four dimensions. In two and three dimensions, DMFT incorrectly predicts a first order phase transition for all bare quartic couplings, which is problematic, because the second order nature of the phase transition of lattice phi^4-theory is crucial for taking the continuum limit. Nevertheless, by extrapolating the behaviour away from the phase transition, one can obtain critical couplings and critical exponents. They differ from those of mean field theory and are much closer to the correct values. In four dimensions the transition is second order for small quartic couplings and turns weakly first order as the coupling increases beyond a tricritical value. In dimensions five and higher, DMFT gives qualitatively correct results, predicts reasonable values for the critical exponents and considerably more accurate critical couplings than standard mean field theory. The approximation works best for small values of the quartic coupling. We investigate the change from first to second order transition in the local limit of DMFT which is computationally much less intensive. We also discuss technical issues related to the convergence of the non-linear self-consistency equation solver and the solution of the effective single-site model using Fourier-space Monte Carlo updates in the presence of a phi^4-interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Euclidean Dynamical Triangulation revisited: is the phase transition really first order?

The transition between the two phases of 4D Euclidean Dynamical Triangulation [1] was long believed to be of second order until in 1996 first order behavior was found for sufficiently large systems [3,4]. However, one may wonder if this finding was affected by the numerical methods used: to control volume fluctuations, in both studies [3,4] an artificial harmonic potential was added to the action; in [4] measurements were taken after a fixed number of accepted instead of attempted moves which introduces an additional error. Finally the simulations suffer from strong critical slowing down which may have been underestimated. In the present work, we address the above weaknesses: we allow the volume to fluctuate freely within a fixed interval; we take measurements after a fixed number of attempted moves; and we overcome critical slowing down by using an optimized parallel tempering algorithm [6]. With these improved methods, on systems of size up to 64k 4-simplices, we confirm that the phase transition is first order.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Scale hierarchy in high-temperature QCD

Because of asymptotic freedom, QCD becomes weakly interacting at high temperature: this is the reason for the transition to a deconfined phase in Yang-Mills theory at temperature $T_c$. At high temperature $T \gg T_c$, the smallness of the running coupling $g$ induces a hierachy betwen the "hard", "soft" and "ultrasoft" energy scales $T$, $g T$ and $g^2 T$. This hierarchy allows for a very successful effective treatment where the "hard" and the "soft" modes are successively integrated out. However, it is not clear how high a temperature is necessary to achieve such a scale hierarchy. By numerical simulations, we show that the required temperatures are extremely high. Thus, the quantitative success of the effective theory down to temperatures of a few $T_c$ appears surprising a posteriori.

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
preprint2012arXiv

A surprise with many-flavor staggered fermions in the strong coupling limit

It is widely believed that chiral symmetry is spontaneously broken at zero temperature in the strong coupling limit of staggered fermions, for any number of colors and flavors. Using Monte Carlo simulations, we show that this conventional wisdom, based on a mean-field analysis, is wrong. For sufficiently many fundamental flavors, chiral symmetry is restored via a bulk, first-order transition. This chirally symmetric phase appears to be analytically connected with the expected conformal window of many-flavor continuum QCD. We perform simulations in the chirally symmetric phase at zero quark mass for various system sizes L, and measure the torelon mass and the Dirac spectrum. We find that all observables scale with L, which is hence the only infrared length scale. Thus, the strong-coupling chirally restored phase appears as a convenient laboratory to study IR-conformality. Finally, we present a conjecture for the phase diagram of lattice QCD as a function of the bare coupling and the number of quark flavors.

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

New algorithms and new results for strong coupling LQCD

We present and compare new types of algorithms for lattice QCD with staggered fermions in the limit of infinite gauge coupling. These algorithms are formulated on a discrete spatial lattice but with continuous Euclidean time. They make use of the exact Hamiltonian, with the inverse temperature beta as the only input parameter. This formulation turns out to be analogous to that of a quantum spin system. The sign problem is completely absent, at zero and non-zero baryon density. We compare the performance of a continuous-time worm algorithm and of a Stochastic Series Expansion algorithm (SSE), which operates on equivalence classes of time-ordered interactions. Finally, we apply the SSE algorithm to a first exploratory study of two-flavor strong coupling lattice QCD, which is manageable in the Hamiltonian formulation because the sign problem can be controlled.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Numerical properties of staggered quarks with a taste-dependent mass term

The numerical properties of staggered Dirac operators with a taste-dependent mass term proposed by Adams [1,2] and by Hoelbling [3] are compared with those of ordinary staggered and Wilson Dirac operators. In the free limit and on (quenched) interacting configurations, we consider their topological properties, their spectrum, and the resulting pion mass. Although we also consider the spectral structure, topological properties, locality, and computational cost of an overlap operator with a staggered kernel, we call attention to the possibility of using the Adams and Hoelbling operators without the overlap construction. In particular, the Hoelbling operator could be used to simulate two degenerate flavors without additive mass renormalization, and thus without fine-tuning in the chiral limit.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Walking near a Conformal Fixed Point: the 2-d O(3) Model at theta near pi as a Test Case

Slowly walking technicolor models provide a mechanism for electroweak symmetry breaking whose nonperturbative lattice investigation is rather challenging. Here we demonstrate walking near a conformal fixed point considering the 2-d lattice O(3) model at vacuum angle $θ\approx π$. The essential features of walking technicolor models are shared by this toy model and can be accurately investigated by numerical simulations. We show results for the running coupling and the beta-function and we perform a finite size scaling analysis of the massgap close to the conformal point.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Continuous Time Monte Carlo for Lattice QCD in the Strong Coupling Limit

We present results for lattice QCD in the limit of infinite gauge coupling, obtained from a worm-type Monte Carlo algorithm on a discrete spatial lattice but with continuous Euclidean time. This is obtained by sending both the anisotropy parameter gamma^2 \sim a/a_t and the number of time-slices N_τto infinity, keeping the ratio γ^2/N_τ\sim aT fixed. The obvious gain is that no continuum extrapolation N_τ-> \infty has to be carried out. Moreover, the algorithm is faster and the sign problem disappears. We compare our computations with those on discrete lattices. We determine the phase diagram as a function of temperature and baryon chemical potential.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Continuous Time Monte Carlo for Lattice QCD in the Strong Coupling Limit

We present results for lattice QCD with staggered fermions in the limit of infinite gauge coupling, obtained from a worm-type Monte Carlo algorithm on a discrete spatial lattice but with continuous Euclidean time. This is achieved by sending both the anisotropy parameter $γ^2\simeq a/\at$ and the number of time-slices $N_τ$ to infinity, keeping the ratio $γ^2/N_τ\simeq aT$ fixed. In this limit, ambiguities arising from the anisotropy parameter $γ$ are eliminated and discretization errors usually introduced by a finite temporal lattice extent $\Nt$ are absent. The obvious gain is that no continuum extrapolation $N_τ\rightarrow \infty$ has to be carried out. Moreover, the algorithm is faster and the sign problem disappears completely. As a first application, we determine the phase diagram as a function of temperature and real and imaginary baryon chemical potential. We compare our computations with those on lattices with discrete Euclidean time. Discretization errors due to finite $\Nt$ in previous studies turn out to be large at low temperatures.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical properties of staggered overlap fermions

We report the results of a numerical study of staggered overlap fermions, following the construction of Adams which reduces the number of tastes from 4 to 2 without fine-tuning. We study the sensitivity of the operator to the topology of the gauge field, its locality and its robustness to fluctuations of the gauge field. We make a first estimate of the computing cost of a quark propagator calculation, and compare with Neuberger's overlap.

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

Equation of State at Finite Density from Imaginary Chemical Potential

We perform two flavor QCD simulations with an imaginary chemical potential and measure derivatives of the pressure up to 4th order as a function of the imaginary chemical potential and the temperature $T \in [0.83 T_c, 2 T_c]$. For temperatures $T \geq T_c$, these derivatives are fitted by a Taylor series in $μ/T$ about $μ=0$. A fit limited to 4th order describes the data poorly at all temperatures, showing that we are sensitive to 6th order contributions. Similarly, a 6th order fit fails for temperatures $T_c \leq T \leq 1.05 T_c$, showing the need for 8th order terms. Thus, our method may offer a computational advantage over the direct measurement of Taylor coefficients at $μ=0$. At temperatures $T \leq T_c$, we fit our data with a hadron resonance gas ansatz. The fit starts to fail at $T \gtrsim 0.95 T_c$. Using our fits, we also reconstruct the equation of state as a function of real quark and isospin chemical potentials.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nuclear Physics from lattice QCD at strong coupling

We study numerically the strong coupling limit of lattice QCD with one flavor of massless staggered quarks. We determine the complete phase diagram as a function of temperature and chemical potential, including a tricritical point. We clarify the nature of the low temperature dense phase, which is strongly bound nuclear matter. This strong binding is explained by the nuclear potential, which we measure. Finally, we determine, from this first-principle limiting case of QCD, the masses of atomic nuclei up to A=12 "carbon".

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Simulating QCD at finite density

In this review, I recall the nature and the inevitability of the "sign problem" which plagues attempts to simulate lattice QCD at finite baryon density. I present the main approaches used to circumvent the sign problem at small chemical potential. I sketch how one can predict analytically the severity of the sign problem, as well as the numerically accessible range of baryon densities. I review progress towards the determination of the pseudo-critical temperature T_c(mu), and towards the identification of a possible QCD critical point. Some promising advances with non-standard approaches are reviewed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The phase diagram of Yang-Mills theory with a compact extra dimension

We present a non-perturbative study of the phase diagram of SU(2) Yang-Mills theory in a five-dimensional spacetime with a compact extra dimension. The non-renormalizable theory is regularized on an anisotropic lattice and investigated through numerical simulations in a regime characterized by a hierarchy between the scale of low-energy physics, the inverse compactification radius, and the cutoff scale. We map out the structure of the phase diagram and the pattern of lines corresponding to fixed values of the ratio between the mass of the fifth component of the gauge field and the non-perturbative mass gap of the four-dimensional modes. We discuss different limits of the model, and comment on the implications of our findings.

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

Nuclear physics from strong coupling QCD

The strong coupling limit (beta_gauge = 0) of QCD offers a number of remarkable research possibilities, of course at the price of large lattice artifacts. Here, we determine the complete phase diagram as a function of temperature T and baryon chemical potential mu_B, for one flavor of staggered fermions in the chiral limit, with emphasis on the determination of a tricritical point and on the T ~ 0 transition to nuclear matter. The latter is known to happen for mu_B substantially below the baryon mass, indicating strong nuclear interactions in QCD at infinite gauge coupling. This leads us to studying the properties of nuclear matter from first principles. We determine the nucleon-nucleon potential in the strong coupling limit, as well as masses m_A of nuclei as a function of their atomic number A. Finally, we clarify the origin of nuclear interactions at strong coupling, which turns out to be a steric effect.

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
preprint2006arXiv

Testing and tuning symplectic integrators for Hybrid Monte Carlo algorithm in lattice QCD

We examine a new 2nd order integrator recently found by Omelyan et al. The integration error of the new integrator measured in the root mean square of the energy difference, $\braΔH^2\ket^{1/2}$, is about 10 times smaller than that of the standard 2nd order leapfrog (2LF) integrator. As a result, the step size of the new integrator can be made about three times larger. Taking into account a factor 2 increase in cost, the new integrator is about 50% more efficient than the 2LF integrator. Integrating over positions first, then momenta, is slightly more advantageous than the reverse. Further parameter tuning is possible. We find that the optimal parameter for the new integrator is slightly different from the value obtained by Omelyan et al., and depends on the simulation parameters. This integrator could also be advantageous for the Trotter-Suzuki decomposition in Quantum Monte Carlo.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Contribution of disconnected diagrams to the hyperfine splitting of charmonium

We discuss the contribution of OZI-suppressed diagrams to the hyperfine splitting of charmonium in lattice QCD. We study valence quark mass regions from strange to charm quark masses. No contribution of the disconnected diagram is seen in the vector meson channel. In the pseudo-scalar channel and for valence quark masses around the strange quark, the disconnected contribution induces a considerable increase of the meson mass. This contribution quickly decreases as the quark mass increases. For charmonium the effect is very small although a decrease of the pseudoscalar mass induced by the disconnected contribution cannot be ruled out.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Comparison of SO(3) and SU(2) lattice gauge theory

The Villain form of SO(3) lattice gauge theory is studied and compared to Wilson's SU(2) theory. The topological invariants in SO(3) which correspond to twisted boundary conditions in SU(2) are discussed and lattice observables are introduced for them. An apparent SO(3) phase with negative adjoint Polyakov loop is explained in terms of these observables. The electric twist free energy, an order parameter for the confinement-deconfinement transition, is measured in both theories to calibrate the temperature. The results indicate that lattices with about 700^4 sites or larger will be needed to study the SO(3) confined phase. Alternative actions are discussed and an analytic path connecting SO(3) and SU(2) lattice gauge theory at weak coupling is exhibited. The relevance for confinement of the centre of the gauge group is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2003arXiv

Observing string breaking with Wilson loops

An uncontroversial observation of adjoint string breaking is proposed, while measuring the static potential from Wilson loops only. The overlap of the Wilson loop with the broken-string state is small, but non-vanishing, so that the broken-string groundstate can be seen if the Wilson loop is long enough. We demonstrate this in the context of the (2+1)d SU(2) adjoint static potential, using an improved version of the Luscher-Weisz exponential variance reduction. To complete the picture we perform the more usual multichannel analysis with two basis states, the unbroken-string state and the broken-string state (two so-called gluelumps). As by-products, we obtain the temperature-dependent static potential measured from Polyakov loop correlations, and the fundamental SU(2) static potential with improved accuracy. Comparing the latter with the adjoint potential, we see clear deviations from Casimir scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

't Hooft Loops, Electric Flux Sectors and Confinement in SU(2) Yang-Mills Theory

We use 't Hooft loops of maximal size on finite lattices to calculate the free energy in the sectors of SU(2) Yang-Mills theory with fixed electric flux as a function of temperature and (spatial) volume. Our results provide evidence for the mass gap. The confinement of electric fluxes in the low temperature phase and their condensation in the high temperature phase are demonstrated. In a surprisingly large scaling window around criticality, the transition is quantitatively well described by universal exponents and amplitude ratios relating the properties of the two phases.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Testing Dimensional Reduction in SU(2) Gauge Theory

At high temperature, every $(d+1)$-dimensional theory can be reformulated as an effective theory in $d$ dimensions. We test the numerical accuracy of this Dimensional Reduction for (3+1)-dimensional SU(2) by comparing perturbatively determined effective couplings with lattice results as the temperature is progressively lowered. We observe an increasing disagreement between numerical and perturbative values from $T=4 T_c$ downwards, which may however be due to somewhat different implementations of dimensional reduction in the two cases.

0 citations0 reviewsNo code linkedOpen access
preprint1997arXiv

A simple ansatz to describe thermodynamic quantities of peptides and proteins at low temperatures

We describe a simple ansatz to approximate the low temperature behavior of proteins and peptides by a mean-field-like model which is analytically solvable. For a small peptide some thermodynamic quantities are calculated and compared with numerical results of an all-atoms simulation. Our approach can be used to determine the weights for a multicanonical simulation of the molecule under consideration.

0 citations0 reviewsNo code linkedOpen access
preprint1992arXiv

The fermion sign problem: A new decoupling transformation, and a new simulation algorithm

We discuss the Fermion sign problem and, by examining a very general Hubbard-Stratonovich (HS) transformation, argue that the sign problem cannot be solved with such methods. We propose a different kind of transformation which, while not solving the sign problem, shows more detailed information about the system. With our transformation it is {\it trivial} to tell which auxiliary field configurations give a positive sign and which give a negative sign. We then discuss briefly various properties of this transformation and construct a new algorithm which with one simulation gives results for a whole range of particle densities and Hubbard $U$ values, positive and negative. Our approach is in excellent agreement with exact calculations.

0 citations0 reviewsNo code linkedOpen access