Source author record

Massimo D'Elia

Massimo D'Elia appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

57works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

57 published item(s)

preprint2022arXiv

Localisation of Dirac eigenmodes and confinement in gauge theories: the Roberge-Weiss transition

Ample numerical evidence from lattice calculations shows a strong connection between the confining properties of gauge theories at finite temperature and the localisation properties of the low-lying Dirac eigenmodes. In this contribution we discuss recent progress on this topic, focussing on results for QCD at imaginary chemical potential $μ_I/T = π$ at temperatures above the Roberge-Weiss transition temperature. These confirm the general picture of low modes turning from delocalised to localised at the deconfinement transition, in a previously unexplored setup with a genuine, physical transition in the presence of dynamical fermions. This further supports the use of Dirac eigenmodes as a tool to investigate the mechanisms behind confinement and the deconfinement transition.

preprint2022arXiv

Towards glueball masses of large-$N$ $\mathrm{SU}(N)$ pure-gauge theories without topological freezing

In commonly used Monte Carlo algorithms for lattice gauge theories the integrated autocorrelation time of the topological charge is known to be exponentially-growing as the continuum limit is approached. This $\mathit{topological}\,\,\textit{freezing}$, whose severity increases with the size of the gauge group, can result in potentially large systematics. To provide a direct quantification of the latter, we focus on $\mathrm{SU}(6)$ Yang--Mills theory at a lattice spacing for which conventional methods associated to the decorrelation of the topological charge have an unbearable computational cost. We adopt the recently proposed $\mathit{parallel}\,\,\mathit{tempering}\,\,\mathit{on}\,\,\mathit{boundary}\,\,\mathit{conditions}$ algorithm, which has been shown to remove systematic effects related to topological freezing, and compute glueball masses with a typical accuracy of $2-5\%$. We observe no sizeable systematic effect in the mass of the first lowest-lying glueball states, with respect to calculations performed at nearly-frozen topological sector.

preprint2021arXiv

Confining and chiral properties of QCD in extremely strong magnetic fields

We investigate, by numerical lattice simulations, the static quark-antiquark potential, the flux tube properties and the chiral condensate for $N_f = 2+1$ QCD with physical quark masses in the presence of strong magnetic fields, going up to $eB = 9$ GeV$^2$, with continuum extrapolated results. The string tension for quark-antiquark separations longitudinal to the magnetic field is suppressed by one order of magnitude at the largest explored magnetic field with respect to its value at zero magnetic background, but is still non-vanishing; in the transverse direction, instead, the string tension is enhanced but seems to reach a saturation at around 50 % of its value at $B = 0$. The flux tube shows a consistent suppression/enhancement of the overall amplitude, with mild modifications of its profile. Finally, we observe magnetic catalysis in the whole range of explored fields with a behavior compatible with a lowest Landau level approximation, in particular with a linear dependence of the chiral condensate on $B$ which is in agreement, within errors, with that already observed for $eB \sim 1$ GeV$^2$.

preprint2021arXiv

Localization properties of Dirac modes at the Roberge-Weiss phase transition

We study the localization properties of the low-lying Dirac eigenmodes in QCD at imaginary chemical potential $\hatμ_I=π$ at temperatures above the Roberge-Weiss transition temperature $T_{\rm RW}$. We find that modes are localized up to a temperature-dependent "mobility edge" and delocalized above it, and that the mobility edge extrapolates to zero at a temperature compatible with $T_{\rm RW}$. This supports the existence of a strong connection between localization of the low Dirac modes and deconfinement, studied here for the first time in a model with a genuine deconfinement transition in the continuum limit in the presence of dynamical fermions.

preprint2021arXiv

Phase diagram of QCD in a magnetic background

We provide numerical evidence that the thermal QCD crossover turns into a first order transition in the presence of large enough magnetic background fields. The critical endpoint is found to be located between $eB = 4$ GeV$^2$ (where the pseudocritical temperature is $T_c = (98 \pm 3)$ MeV) and $eB = 9$ GeV$^2$ (where the critical temperature is $T_c = (63 \pm 5)$ MeV). Results are based on the analysis of quark condensates and number susceptibilities, determined by lattice simulations of $N_f = 2+1$ QCD at the physical point, discretized with three different lattice spacings, $a = 0.114, 0.086$ and $0.057$ fm, via rooted stout staggered fermions and a Symanzik tree level improved pure gauge action. We also present preliminary results regarding the confining properties of the thermal theory, suggesting that they could change drastically going across the phase transition

preprint2020arXiv

$θ$-dependence in the small-$N$ limit of $2d$ $CP^{N-1}$ models

We present a systematic numerical study of $θ$-dependence around $θ=0$ in the small-$N$ limit of $2d$ $CP^{N-1}$ models, aimed at clarifying the possible presence of a divergent topological susceptibility in the continuum limit. We follow a twofold strategy, based on one side on direct simulations for $N = 2$ and $N = 3$ on lattices with correlation lengths up to $O(10^2)$, and on the other side on the small-$N$ extrapolation of results obtained for $N$ up to $9$. Based on that, we provide conclusive evidence for a finite topological susceptibility at $N = 3$, with a continuum estimate $ξ^2 χ= 0.110(5)$. On the other hand, results obtained for $N = 2$ are still inconclusive: they are consistent with a logarithmically divergent continuum extrapolation, but do not yet exclude a finite continuum value, $ξ^2 χ\sim 0.4$, with the divergence taking place for $N$ slightly below 2 in this case. Finally, results obtained for the non-quadratic part of $θ$-dependence, in particular for the so-called $b_2$ coefficient, are consistent with a $θ$-dependence matching that of the Dilute Instanton Gas Approximation at the point where $ξ^2 χ$ diverges.

preprint2020arXiv

Quantum computation of thermal averages in the presence of a sign problem

We illustrate the application of Quantum Computing techniques to the investigation of the thermodynamical properties of a simple system, made up of three quantum spins with frustrated pair interactions and affected by a hard sign problem when treated within classical computational schemes. We show how quantum algorithms completely solve the problem, and discuss how this can apply to more complex systems of physical interest, with emphasis on the possible systematics and on their control.

preprint2020arXiv

Reconfinement, localization and thermal monopoles in $SU(3)$ trace-deformed Yang-Mills theory

We study, by means of numerical lattice simulations, the properties of the reconfinement phase transition taking place in trace deformed $SU(3)$ Yang-Mills theory defined on $\mathbb{R}^3\times S^1$, in which center symmetry is recovered even for small compactification radii. We show, by means of a finite size scaling analysis, that the reconfinement phase transition is first-order, like the usual $SU(3)$ thermal phase transition. We then investigate two different physical phenomena, which are known to characterize the standard confinement/deconfinement phase transition, namely the condensation of thermal magnetic monopoles and the change in the localization properties of the eigenmodes of the Dirac operator. Regarding the latter, we show that the mobility edge signalling the Anderson-like transition in the Dirac spectrum vanishes as one enters the reconfined phase, as it happens in the standard confined phase. Thermal monopoles, instead, show a peculiar behavior: their density decreases going through reconfinement, at odds with the standard thermal theory; nonetheless, they condense at reconfinement, like at the usual confinement transition. The coincidence of monopole condensation and Dirac mode delocalization, even in a framework different from that of the standard confinement transition, suggests the existence of a strict link between them.

preprint2019arXiv

$θ$-dependence and center symmetry in Yang-Mills theories

We investigate the relation between the realization of center symmetry and the dependence on the topological parameter $θ$ in $SU(N)$ Yang-Mills theories, exploiting trace deformations as a tool to regulate center symmetry breaking in a theory with a small compactified direction. We consider, in particular, $SU(4)$ gauge theory, which admits two possible independent deformations, and study, as a first step, its phase diagram in the deformation plane for two values of the inverse compactified radius going up to $L^{-1} \sim 500$ MeV, comparing the predictions of the effective 1-loop potential of the Polyakov loop with lattice results. The $θ$-dependence of the various phases is then addressed, up to the fourth order in $θ$, by numerical simulations: results are found to coincide, within statistical errors, with those of the standard confined phase iff center symmetry is completely restored and independently of the particular way this happens, i.e. either by local suppression of the Polyakov loop traces or by long range disorder.

preprint2019arXiv

Topology via Spectral Projectors with Staggered Fermions

The spectral projectors method is a way to obtain a theoretically well posed definition of the topological susceptibility on the lattice. Up to now this method has been defined and applied only to Wilson fermions. The goal of this work is to extend the method to staggered fermions, giving a definition for the staggered topological susceptibility and testing it in the pure $SU(3)$ gauge theory. Besides, we also generalize the method to higher-order cumulants of the topological charge distribution.

preprint2016arXiv

$θ$ dependence in $SU(3)$ Yang-Mills theory from analytic continuation

We investigate the topological properties of the $SU(3)$ pure gauge theory by performing numerical simulations at imaginary values of the $θ$ parameter. By monitoring the dependence of various cumulants of the topological charge distribution on the imaginary part of $θ$ and exploiting analytic continuation, we determine the free energy density up to the sixth order order in $θ$, $f(θ,T) = f(0,T) + {1\over 2} χ(T) θ^2 (1 + b_2(T) θ^2 + b_4(T) θ^4 + O(θ^6))$. That permits us to achieve determinations with improved accuracy, in particular for the higher order terms, with control over the continuum and the infinite volume extrapolations. We obtain $b_2=-0.0216(15)$ and $|b_4|\lesssim 4\times 10^{-4}$.

preprint2016arXiv

$θ$ dependence of 4D $SU(N)$ gauge theories in the large-$N$ limit

We study the large-$N$ scaling behavior of the $θ$ dependence of the ground-state energy density $E(θ)$ of four-dimensional (4D) $SU(N)$ gauge theories and two-dimensional (2D) $CP^{N-1}$ models, where $θ$ is the parameter associated with the Lagrangian topological term. We consider its $θ$ expansion around $θ=0$, $E(θ)-E(0) = {1\over 2}χ\,θ^2 ( 1 + b_2 θ^2 + b_4θ^4 +\cdots)$ where $χ$ is the topological susceptibility and $b_{2n}$ are dimensionless coefficients. We focus on the first few coefficients $b_{2n}$, which parametrize the deviation from a simple Gaussian distribution of the topological charge at $θ=0$. We present a numerical analysis of Monte Carlo simulations of 4D $SU(N)$ lattice gauge theories for $N=3,\,4,\,6$ in the presence of an imaginary $θ$ term. The results provide a robust evidence of the large-$N$ behavior predicted by standard large-$N$ scaling arguments, i.e. $b_{2n}= O(N^{-2n})$. In particular, we obtain $b_2=\bar{b}_2/N^2 + O(1/N^4)$ with $\bar{b}_2=-0.23(3)$. We also show that the large-$N$ scaling scenario applies to 2D $CP^{N-1}$ models as well, by an analytic computation of the leading large-$N$ dependence.

preprint2016arXiv

Axion phenomenology and $θ$-dependence from $N_f = 2+1$ lattice QCD

We investigate the topological properties of $N_f = 2+1$ QCD with physical quark masses, both at zero and finite temperature. We adopt stout improved staggered fermions and explore a range of lattice spacings $a \sim 0.05 - 0.12$ fm. At zero temperature we estimate both finite size and finite cut-off effects, comparing our continuum extrapolated results for the topological susceptibility $χ$ with predictions from chiral perturbation theory. At finite temperature, we explore a region going from $T_c$ up to around $4\, T_c$, where we provide continuum extrapolated results for the topological susceptibility and for the fourth moment of the topological charge distribution. While the latter converges to the dilute instanton gas prediction the former differs strongly both in the size and in the temperature dependence. This results in a shift of the axion dark matter window of almost one order of magnitude with respect to the instanton computation.

preprint2016arXiv

Field-strength correlators for QCD in a magnetic background

We present the results of an exploratory study (by means of Monte Carlo simulations on the lattice) of the properties of the gauge-invariant two-point correlation functions of the gauge-field strengths for $N_f=2$ QCD at zero temperature and in the presence of a magnetic background field: the analysis provides evidence for the emergence of anisotropies in the nonperturbative part of the correlators and for an increase of the gluon condensate as a function of the external magnetic field.

preprint2016arXiv

Gauge-invariant field-strength correlators for QCD in a magnetic background

We consider the properties of the gauge-invariant two-point correlation functions of the gauge-field strengths for QCD in the presence of a magnetic background field. We discuss the general structure of the correlators in this case and provide the results of an exploratory lattice study for $N_f = 2$ QCD discretized with unimproved staggered fermions. Our analysis provides evidence for the emergence of anisotropies in the non-perturbative part of the correlators and for an increase of the gluon condensate as a function of the external magnetic field.

preprint2016arXiv

Magnetic field effects on the static quark potential at zero and finite temperature

We investigate the static $Q\bar{Q}$ potential at zero and finite temperature in the presence of a constant and uniform external magnetic field $\vec{B}$, for several values of the lattice spacing and for different orientations with respect to $\vec{B}$. As a byproduct, we provide continuum limit extrapolated results for the string tension, the Coulomb coupling and the Sommer parameter at $T = 0$ and $B = 0$. We confirm the presence in the continuum of a $B$-induced anisotropy, regarding essentially the string tension, for which it is of the order of 15\% at $|e| B \sim 1~{\rm GeV}^2$ and would suggest, if extrapolated to larger fields, a vanishing string tension along the magnetic field for $|e| B \gtrsim 4$ GeV$^2$. The angular dependence for $|e| B \lesssim 1$ GeV$^2$ can be nicely parametrized by the first allowed term in an angular Fourier expansion, corresponding to a quadrupole deformation. Finally, for $T \neq 0$, the main effect of the magnetic field is a general suppression of the string tension, leading to a precocious loss of the confining properties: this happens even before the appearance of inverse magnetic catalysis in the chiral condensate, supporting the idea that the influence of the magnetic field on the confining properties is the leading effect originating the decrease of $T_c$ as a function of $B$.

preprint2016arXiv

Roberge-Weiss endpoint at the physical point of $N_f = 2+1$ QCD

We study the phase diagram of $N_f = 2+1$ QCD in the $T - μ_B$ plane and investigate the critical point corresponding to the onset of the Roberge-Weiss transition, which is found for imaginary values of $μ_B$. We make use of stout improved staggered fermions and of the tree level Symanzik gauge action, and explore four different sets of lattice spacings, corresponding to $N_t = 4,6,8,10$, and different spatial sizes, in order to assess the universality class of the critical point. The continuum extrapolated value of the endpoint temperature is found to be $T_{\rm RW} = 208(5)$ MeV, i.e. $T_{\rm RW}/T_c \sim 1.34(7)$, where $T_c$ is the chiral pseudocritical temperature at zero chemical potential, while our finite size scaling analysis, performed on $N_t = 4$ and $N_t = 6$ lattices, provides evidence for a critical point in the $3d$ Ising universality class.

preprint2015arXiv

Curvature of the chiral pseudo-critical line in QCD

We determine the curvature of the pseudo-critical line of strong interactions by means of numerical simulations at imaginary chemical potentials. We consider $N_f=2+1$ stout improved staggered fermions with physical quark masses and the tree level Symanzik gauge action, and explore two different sets of lattice spacings, corresponding to temporal extensions $N_t = 6$ and $N_t = 8$. Both the renormalized chiral condensate and the renormalized chiral susceptibility are used to locate the transition. The determinations obtained from the two quantities are in good agreement, a preliminary continuum extrapolation yields $κ= 0.013(2)(1)$. We also investigate the impact of a non-zero strange quark chemical potential and compare our results to previous determinations in the literature, discussing the possible sources of systematic errors affecting the various procedures.

preprint2015arXiv

Curvature of the chiral pseudo-critical line in QCD: continuum extrapolated results

We determine the curvature of the pseudo-critical line of strong interactions by means of numerical simulations at imaginary chemical potentials. We consider $N_f=2+1$ stout improved staggered fermions with physical quark masses and the tree level Symanzik gauge action, and explore four different sets of lattice spacings, corresponding to $N_t = 6,8,10,12$, in order to extrapolate results to the continuum limit. Our final estimate is $κ= 0.0135(20)$.

preprint2015arXiv

Heavy quarkonia in strong magnetic fields

We investigate the influence of a homogeneus and constant strong external magnetic field on the heavy-meson spectrum. Quarkonium states $c\bar{c}$ and $b\bar{b}$ are described within a non-relativistic framework and by means of a suitable potential model based on the Cornell parametrization. In particular, in this work we propose a model which takes into account the possible anisotropies emerging at the level of the static quark-antiquark potential, as observed in recent lattice studies. The investigation is perfomed both with and without taking into account the anisotropy of the static potential, in order to better clarify its effects

preprint2014arXiv

Anisotropy of the quark-antiquark potential in a magnetic field

We investigate the static $\overline{Q}Q$-potential for $N_f = 2+1$ QCD at the physical point in the presence of a constant and uniform external magnetic field. The potential is found to be anisotropic and steeper in the directions transverse to the magnetic field than in the longitudinal one. In particular, when compared to the standard case with zero background field, the string tension increases (decreases) in the transverse (longitudinal) direction, while the absolute value of the Coulomb coupling and the Sommer parameter show an opposite behavior.

preprint2014arXiv

Comparison of the gradient flow with cooling in $SU(3)$ pure gauge theory

The gradient (Wilson) flow has been introduced recently in order to provide a solid theoretical framework for the smoothing of ultraviolet noise in lattice gauge configurations. It is interesting to ask how it compares with other, more heuristic and numerically cheaper smoothing techniques, such as standard cooling. In this study we perform such a comparison, focusing on observables related to topology. We show that, already for moderately small lattice spacings, standard cooling and the gradient flow lead to equivalent results, both for average quantities and configuration by configuration.

preprint2014arXiv

Magnetic susceptibility and equation of state of N_f = 2+1 QCD with physical quark masses

We determine the free energy of strongly interacting matter as a function of an applied constant and uniform magnetic field. We consider N_f = 2+1 QCD with physical quark masses, discretized on a lattice by stout improved staggered fermions and a tree level improved Symanzik pure gauge action, and explore three different lattice spacings. For magnetic fields of the order of those produced in non-central heavy ion collisions (eB ~ 0.1 GeV^2) strongly interacting matter behaves like a medium with a linear response, and is paramagnetic both above and below the deconfinement transition, with a susceptibility which steeply rises in the deconfined phase. We compute the equation of state, showing that the relative increase in the pressure due to the magnetic field gets larger around the transition, and of the order of 10 % for eB ~ 0.1 GeV^2.

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}$.

preprint2014arXiv

Universal scaling effects of a temperature gradient at first-order transitions

We study the effects of smooth inhomogeneities at first-order transitions. We show that a temperature gradient at a thermally-driven first-order transition gives rise to nontrivial universal scaling behaviors with respect to the length scale of the variation of the local temperature T(x). We propose a scaling ansatz to describe the crossover region at the surface where T(x)=Tc, where the typical discontinuities of a first-order transition are smoothed out. The predictions of this scaling theory are checked, and get strongly supported, by numerical results for the 2D Potts models, for a sufficiently large number of q-states to have first-order transitions. Comparing with analogous results at the 2D Ising transition, we note that the scaling behaviors induced by a smooth inhomogeneity appear quite similar in first-order and continuous transitions.

preprint2013arXiv

Change of theta dependence in 4D SU(N) gauge theories across the deconfinement transition

We investigate the dependence of four-dimensional SU(N) gauge theories on the topological theta term at finite temperature and, in particular, across the deconfinement transition. For this purpose, we exploit the lattice formulation of the theory and present numerical results for the expansion of the free energy up to O(theta^6), for N=3 and N=6. Our numerical analysis shows that the theta dependence of 4D SU(N) gauge theory experiences a drastic change across the deconfinement transition: the low-temperature phase is characterized by a large-N scaling with theta/N as relevant variable, while in the high-temperature phase the scaling variable is just theta and the free energy is essentially determined by the instanton-gas approximation. The crossover between the two different behaviours gets sharper with increasing N, suggesting that the instanton-gas regime sets in just above Tc at large N.

preprint2013arXiv

Phase diagram of the 4D U(1) model at finite temperature

We explore the phase diagram of the 4D compact U(1) gauge theory at finite temperature as a function of the gauge coupling and of the compactified Euclidean time dimension L_t. We show that the strong-to-weak coupling transition, which is first order at T=0 (L_t=\infty), becomes second order for high temperatures, i.e. for small values of L_t, with a tricritical temporal size \bar{L_t} located between 5 and 6. The critical behavior around the tricritical point explains and reconciles previous contradictory evidences found in the literature.

preprint2013arXiv

Phase diagram of Yang-Mills theories in the presence of a theta term

We study the phase diagram of non-Abelian pure gauge theories in the presence of a topological theta term. The dependence of the deconfinement temperature on theta is determined on the lattice both by analytic continuation and by reweighting, obtaining consistent results. The general structure of the diagram is discussed on the basis of large-N considerations and of the possible analogies and dualities existing with the phase diagram of QCD in presence of an imaginary baryon chemical potential.

preprint2013arXiv

Susceptibility of the QCD vacuum to CP-odd electromagnetic background fields

We investigate two flavor QCD in presence of CP-odd electromagnetic background fields and determine, by means of lattice QCD simulations, the induced effective theta term to the first order in the scalar product of E and B. We employ a rooted staggered discretization and study lattice spacings down to 0.1 fm and Goldstone pion masses around 480 MeV. In order to deal with a positive measure, we consider purely imaginary electric fields and real magnetic fields, then exploiting analytic continuation. Our results are relevant to a description of the effective pseudoscalar QED-QCD interactions.

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.

preprint2013arXiv

The magnetic susceptibility in QCD

Recently much work has been devoted to the study of QCD coupled to a background magnetic field. Strongly interacting matter acts as a magnetic medium and it is natural to study the properties of this medium, in particular to understand if it behaves like a diamagnetic or a paramagnetic material. A serious difficulty in studying these properties by means of LQCD simulations is the quantization of the magnetic field in a toroidal geometry. We will expose a method to overcome this difficulty and we will present data obtained for the $N_f=2$ theory that show that the QCD medium is paramagnetic.

preprint2013arXiv

The Maximal Abelian Gauge in SU(N) Gauge Theories and Thermal Monopoles for N = 3

We discuss and propose a proper extension of the Abelian projection based on the Maximal Abelian Gauge to SU(N) gauge theories. Based, on that, we investigate the properties of thermal Abelian monopoles in the deconfined phase of the SU(3) pure gauge theory. Such properties are very similar to those already found for SU(2), confirming the relevance of the magnetic component close to Tc and the possible condensation of thermal monopoles as the deconfinement temperature is crossed from above. Moreover, we study the correlation functions among monopoles related to different U(1) subgroups, which show interesting features and reveal the presence of non-trivial interactions.

preprint2013arXiv

Theta dependence of 4D SU(N) gauge theories at finite temperature

We report a study of the dependence of 4D SU(N) gauge theories on the topological theta term at finite temperature, and in particular in the large-N limit. We show that the theta dependence drastically changes across the deconfinement transition. The low-temperature phase is characterized by a large-N scaling with theta/N as relevant variable, while in the high-temperature phase the free energy is essentially determined by the dilute instanton-gas approximation, with a simple theta dependence of the free-energy density proportional to 1-Cos(theta).

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.

preprint2012arXiv

Lattice QCD Simulations in External Background Fields

We discuss recent results and future prospects regarding the investigation, by lattice simulations, of the non-perturbative properties of QCD and of its phase diagram in presence of magnetic or chromomagnetic background fields. After a brief introduction to the formulation of lattice QCD in presence of external fields, we focus on studies regarding the effects of external fields on chiral symmetry breaking, on its restoration at finite temperature and on deconfinement. We conclude with a few comments regarding the effects of electromagnetic background fields on gluodynamics.

preprint2012arXiv

The disorder parameter of dual superconductivity in QCD revisited

We discover the origin of the pathologies of the disorder parameter used in previous papers to detect dual superconductivity of QCD vacuum, and we remove them by defining an improved disorder parameter. A check of the approach is made by numerical simulations of SU(2) gauge theory, which demonstrate that the approach is consistent and with it that deconfinement is a transition from dual superconductor to normal.

preprint2012arXiv

Theta dependence of the deconfinement temperature in Yang-Mills theories

We determine the theta dependence of the deconfinement temperature of SU(3) pure gauge theory, finding that it decreases in presence of a topological theta term. We do that by performing lattice simulations at imaginary theta, then exploiting analytic continuation. We also give an estimate of such dependence in the limit of a large number of colors N, and compare it with our numerical results.

preprint2011arXiv

Chiral Properties of Strong Interactions in a Magnetic Background

We investigate the chiral properties of QCD in presence of a magnetic background field and in the low temperature regime, by lattice numerical simulations of N_f = 2 QCD. We adopt a standard staggered discretization, with a pion mass around 200 MeV, and explore a range of magnetic fields (180 MeV)^2 \leq|e|B \leq (700 MeV)^2, in which we study magnetic catalysis, i.e. the increase of chiral symmetry breaking induced by the background field. We determine the dependence of the chiral condensate on the external field, compare our results with existing model predictions and show that a substantial contribution to magnetic catalysis comes from the modified distribution of non-Abelian gauge fields, induced by the magnetic field via dynamical quark loop effects.

preprint2011arXiv

QCD Phase Transition in a Strong Magnetic Background

We investigate the properties of the deconfining/chiral restoring transition for two flavor QCD in presence of a uniform background magnetic field. We adopt a standard staggered discretization of the fermion action and a lattice spacing of the order of 0.3 fm. We explore different values of the bare quark mass, corresponding to pion masses in the range 200 - 480 MeV, and magnetic fields up to |e|B ~ 0.75 GeV^2. The deconfinement and chiral symmetry restoration temperatures remain compatible with each other and rise very slightly (< 2 % for our largest magnetic field) as a function of the magnetic field. On the other hand, the transition seems to become sharper as the magnetic field increases.

preprint2011arXiv

QCD simulations with staggered fermions on GPUs

We report on our implementation of the RHMC algorithm for the simulation of lattice QCD with two staggered flavors on Graphics Processing Units, using the NVIDIA CUDA programming language. The main feature of our code is that the GPU is not used just as an accelerator, but instead the whole Molecular Dynamics trajectory is performed on it. After pointing out the main bottlenecks and how to circumvent them, we discuss the obtained performances. We present some preliminary results regarding OpenCL and multiGPU extensions of our code and discuss future perspectives.

preprint2011arXiv

The Roberge-Weiss endpoint in N_f = 2 QCD

We present the results of extensive simulations regarding the critical behavior at the endpoint of the Roberge-Weiss transition for N_f = 2 QCD. We confirm early evidence, presented in arXiv:0909.0254, according to which the Roberge-Weiss endpoint is first order in the limit of large or small quark masses, and second order for intermediate masses. A systematic study of the transition strength as a function of the quark mass in the first order regions, permits us to estimate the tricritical values of the quark mass separating the second order region from the first order ones.

preprint2011arXiv

Toward a well defined monopole creation operator

The lattice implementation of monopole creation operator proves to have problems related to bulk transitions that can possibly affect the interpretation of its mean value as an order parameter for monopole condensation. Preliminary evidence is presented that these unexpected behaviours are in fact only due to lattice artefacts and do not spoil the physical interpretation of the monopole operator.

preprint2010arXiv

Detecting monopoles on the lattice

We address the issue why the number and the location of magnetic monopoles detected on lattice configurations are gauge dependent, in contrast with the physical expectation that monopoles have a gauge invariant status. By use of the Non-Abelian Bianchi Identities we show that monopoles are gauge invariant, but the efficiency of the technique usually adopted to detect them depends on the choice of the gauge in a well understood way. In particular we have studied a class of gauges which interpolates between the Maximal Abelian gauge, where all monopoles are observed, and the Landau gauge, where all monopoles escape detection.

preprint2010arXiv

Non abelian Bianchi identities, monopoles and gauge invariance

A direct connection is proved between the Non-Abelian Bianchi Identities and the Abelian Bianchi identities for the 't Hooft tensor in a generic gauge; the existence of a magnetic current is related to the violation of NABI's. Using this relation it is shown that not all gauges are equivalent to detect monopoles on the lattice, that e.g. the Maximal Abelian Gauge is a legitimate choice while the Landau gauge is not. Nevertheless monopole condensation is found to be a gauge invariant property.

preprint2010arXiv

On the analytic continuation of the critical line

We perform a numerical study of the systematic effects involved in the determination of the critical line at real baryon chemical potential by analytic continuation from results obtained at imaginary chemical potentials. We present results obtained in a theory free of the sign problem, three-color QCD with finite isospin chemical potential, and comment on general features which could be relevant also to the continuation of the critical line in real QCD at finite baryon density.

preprint2010arXiv

The order of the Roberge-Weiss endpoint (finite size transition) in QCD

We consider the endpoint of the Roberge-Weiss (RW) first order transition line present for imaginary baryon chemical potentials. We remark that it coincides with the finite size transition relevant in the context of large $N_c$ QCD and study its order in the theory with two degenerate flavors. The RW endpoint is first order in the limit of large and small quark masses, while it weakens for intermediate masses where it is likely in the Ising 3d universality class. Phenomenological implications and further speculations about the QCD phase diagram are discussed.

preprint2010arXiv

The phase diagram of QCD with four degenerate quarks

We revisit the determination of the pseudo-critical line of QCD with four degenerate quarks at non-zero temperature and baryon density by the method of analytic continuation. We determine the pseudo-critical couplings at imaginary chemical potentials by high-statistics Monte Carlo simulations and reveal deviations from the simple quadratic dependence on the chemical potential visible in earlier works on the same subject. Finally, we discuss the implications of our findings for the shape of the pseudo-critical line at real chemical potential, comparing different possible extrapolations.

preprint2010arXiv

The three-dimensional, three state Potts model in a negative external field

We investigate the critical behaviour of the three-dimensional, three state Potts model in presence of a negative external field h, i.e. disfavouring one of the three states. A genuine phase transition is present for all values of |h|, corresponding to the spontaneous breaking of a residual Z_2 symmetry. The transition is first/second order respectively for small/large values of |h|, with a tricritical field h_tric separating the two regimes. We provide, using different and consistent approaches, a precise determination of h_tric; we also compare with previous studies and discuss the relevance of our investigation to analogous studies of the QCD phase diagram in presence of an imaginary chemical potential.

preprint2010arXiv

Thermal Monopole Condensation and Confinement in finite temperature Yang-Mills Theories

We investigate the connection between Color Confinement and thermal Abelian monopoles populating the deconfined phase of SU(2) Yang-Mills theory, by studying how the statistical properties of the monopole ensemble change as the confinement/deconfinement temperature is approached from above. In particular we study the distribution of monopole currents with multiple wrappings in the Euclidean time direction, corresponding to two or more particle permutations, and show that multiple wrappings increase as the deconfinement temperature is approached from above, in a way compatible with a condensation of such objects happening right at the deconfining transition. We also address the question of the thermal monopole mass, showing that different definitions give consistent results only around the transition, where the monopole mass goes down and becomes of the order of the critical temperature itself.

preprint2009arXiv

Thermodynamics of two flavor QCD from imaginary chemical potentials

We study QCD thermodynamics in presence of two independent imaginary chemical potentials coupled to two degenerate flavors of staggered quarks. Analytic continuation is used to determine non-linear susceptibilities, to test the Hadron Resonance Gas (HRG) model below the zero density critical temperature, Tc, and to determine the average phase factor of the fermion determinant. Deviations from HRG predictions, of the order of a few percent, are clearly visible for temperatures T > 0.95 Tc. The determination of non-linear susceptibilities, using different interpolating functions for analytic continuation, gives consistent results and in agreement with Taylor expansion computations, apart from some systematic effects at or right above Tc. Results for the average phase factor are compared with the predictions of Chiral Perturbation Theory; below Tc we are able to distinguish the contribution of different hadron states, which is positive (i.e. tends to mitigate the sign problem) in the case of baryons.

preprint2007arXiv

Analytic continuation from imaginary to real chemical potential in two-color QCD

The method of analytic continuation from imaginary to real chemical potential is one of the most powerful tools to circumvent the sign problem in lattice QCD. Here we test this method in a theory, 2-color QCD, which is free from the sign problem. We find that the method gives reliable results, within appropriate ranges of the chemical potential, and that a considerable improvement can be achieved if suitable functions are used to interpolate data with imaginary chemical potential.

preprint2004arXiv

QCD thermodynamics from an imaginary mu_B: results on the four flavor lattice model

We study four flavor QCD at nonzero temperature and density by analytic continuation from an imaginary chemical potential. The explored region is T = 0.95 T_c < T < 3.5 T_c, and the baryochemical potentials range from 0 to approx. 500 MeV. Observables include the number density, the order parameter for chiral symmetry, and the pressure, which is calculated via an integral method at fixed temperature and quark mass. The simulations are carried out on a 16^3 X 4 lattice, and the mass dependence of the results is estimated by exploiting the Maxwell relations. In the hadronic region we confirm that the results are consistent with a simple resonance hadron gas model, and we estimate the critical density by combining the results for the number density with those for the critical line. In the hot phase, above the endpoint of the Roberge-Weiss transition T_E approx 1.1 T_c the results are consistent with a free lattice model with a fixed effective number of flavor slightly different from four. We confirm that confinement and chiral symmetry are coincident by a further analysis of the critical line, and we discuss the interrelation between thermodynamics and critical behavior. We comment on the strength and weakness of the method, and propose further developments.

preprint2003arXiv

Field theoretical approach to the study of theta dependence in Yang-Mills theories on the lattice

We discuss the extension of the field theoretical approach, already used in the lattice determination of the topological susceptibility, to the computation of further terms in the expansion of the ground state energy $F(θ)$ around $θ= 0$ in SU(N) Yang-Mills theories. In particular we determine the fourth order term in the expansion for SU(3) pure gauge theory and compare our results with previous cooling determinations. In the last part of the paper we make some considerations about the nature of the ultraviolet fluctuations responsible for the renormalization of the lattice topological charge correlation functions; in particular we propose and test an ansatz which leads to improved estimates of the fourth and higher order terms in the expansion of F(θ).