Source author record

Claudio Bonati

Claudio Bonati 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

55works
14topics
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

55 published item(s)

preprint2022arXiv

Color-flavor reflection in the continuum limit of two-dimensional lattice gauge theories with scalar fields

We address the interplay between local and global symmetries in determining the continuum limit of two-dimensional lattice scalar theories characterized by $SO(N_c)$ gauge symmetry and non-Abelian $O(N_f)$ global invariance. We argue that, when a quartic interaction is present, the continuum limit of these model corresponds in some cases to the gauged non-linear $σ$ model field theory associated with the real Grassmannian manifold $SO(N_f)/(SO(N_c)\times SO(N_f-N_c)$), which is characterized by the invariance under the color-flavor reflection $N_c\leftrightarrow N_f-N_c$. Monte Carlo simulations and Finite-Size Scaling analyses, performed for $N_f=7$ and several values of $N_c$, confirm the emergence of the color-flavor reflection symmetry in the scaling limit, and support the identification of the continuum limit.

preprint2022arXiv

Phase diagram and Higgs phases of 3D lattice SU(Nc) gauge theories with multiparameter scalar potentials

We consider three-dimensional lattice SU(Nc) gauge theories with degenerate multicomponent (Nf>1) complex scalar fields that transform under the fundamental representation of the gauge SU(Nc) group and of the global U(Nf) invariance group, interacting with the most general quartic potential compatible with the global (flavor) and gauge (color) symmetries. We investigate the phase diagrams, identifying the low-temperature Higgs phases and their global and gauge symmetries, and the critical behaviors along the different transition lines. In particular, we address the role of the quartic scalar potential, which determines the Higgs phases and the corresponding symmetry-breaking patterns. Our study is based on the analysis of the minimum-energy configurations and on numerical Monte Carlo simulations. Moreover, we investigate whether some of the transitions observed in the lattice model can be related to the behavior of the renormalization-group flow of the continuum field theory with the same symmetries and field content around its stable charged fixed points. For Nc=2, numerical results are consistent with the existence of charged critical behaviors for Nf > Nf*, with 20 < Nf* < 40.

preprint2022arXiv

Scalar gauge-Higgs models with discrete Abelian symmetry groups

We investigate the phase diagram and the nature of the phase transitions of three-dimensional lattice gauge-Higgs models obtained by gauging the Z_N subgroup of the global Z_q invariance group of the Z_q clock model (N is a submultiple of q). The phase diagram is generally characterized by the presence of three different phases, separated by three distinct transition lines. We investigate the critical behavior along the two transition lines characterized by the ordering of the scalar field. Along the transition line separating the disordered-confined phase from the ordered-deconfined phase, standard arguments within the Landau-Ginzburg-Wilson framework predict that the behavior is the same as in a generic ferromagnetic model with Z_p global symmetry, p being the ratio q/N. Thus, continuous transitions belong to the Ising and to the O(2) universality class for p=2 and p>3, respectively, while for p=3 only first-order transitions are possible. The results of Monte Carlo simulations confirm these predictions. There is also a second transition line, which separates two phases in which gauge fields are essentially ordered. Along this line we observe the same critical behavior as in the Z_q clock model, as it occurs in the absence of gauge fields.

preprint2021arXiv

Berezinskii-Kosterlitz-Thouless transitions in two-dimensional lattice SO($N_c$) gauge theories with two scalar flavors

We study the phase diagram and critical behavior of a two-dimensional lattice SO($N_c$) gauge theory ($N_c \ge 3$) with two scalar flavors, obtained by partially gauging a maximally O($2N_c$) symmetric scalar model. The model is invariant under local SO($N_c$) and global O(2) transformations. We show that, for any $N_c \ge 3$, it undergoes finite-temperature Berezinskii-Kosterlitz-Thouless (BKT) transitions, associated with the global Abelian O(2) symmetry. The transition separates a high-temperature disordered phase from a low-temperature spin-wave phase where correlations decay algebraically (quasi-long range order). The critical properties at the finite-temperature BKT transition and in the low-temperature spin-wave phase are determined by means of a finite-size scaling analysis of Monte Carlo data.

preprint2021arXiv

Multicritical point of the three-dimensional Z_2 gauge Higgs model

We investigate the multicritical behavior of the three-dimensional Z_2 gauge Higgs model, at the multicritical point (MCP) of its phase diagram, where one first-order transition line and two continuous Ising-like transition lines meet. The duality properties of the model determine some features of the multicritical behavior at the MCP located along the self-dual line. Moreover, we argue that the system develops a multicritical XY behavior at the MCP, which is controlled by the stable XY fixed point of the three-dimensional multicritical Landau-Ginzburg-Wilson field theory with two competing scalar fields associated with the continuous Z_2 transition lines meeting at the MCP. This implies an effective enlargement of the symmetry of the multicritical modes at the MCP, to the continuous group O(2). We also provide some numerical results to support the multicritical XY scenario.

preprint2020arXiv

Asymptotic low-temperature behavior of two-dimensional RP$^{N-1}$ models

We investigate the low-temperature behavior of two-dimensional (2D) RP$^{N-1}$ models, characterized by a global O($N$) symmetry and a local ${\mathbb Z}_2$ symmetry. For $N=3$ we perform large-scale simulations of four different 2D lattice models: two standard lattice models and two different constrained models. We also consider a constrained mixed O(3)-RP$^2$ model for values of the parameters such that vector correlations are always disordered. We find that all these models show the same finite-size scaling (FSS) behavior, and therefore belong to the same universality class. However, these FSS curves differ from those computed in the 2D O(3) $σ$ model, suggesting the existence of a distinct 2D RP$^2$ universality class. We also performed simulations for $N=4$, and the corresponding FSS results also support the existence of an RP$^3$ universality class, different from the O(4) one.

preprint2020arXiv

Asymptotic low-temperature critical behavior of two-dimensional multiflavor lattice SO(Nc) gauge theories

We address the interplay between global and local gauge nonabelian symmetries in lattice gauge theories with multicomponent scalar fields. We consider two-dimensional lattice scalar nonabelian gauge theories with a local SO(Nc) (Nc >= 3) and a global O(Nf) invariance, obtained by partially gauging a maximally O(Nf x Nc)-symmetric multicomponent scalar model. Correspondingly, the scalar fields belong to the coset S(Nf Nc-1)/SO(Nc), where S(N) is the N-dimensional sphere. In agreement with the Mermin-Wagner theorem, these lattice SO(Nc) gauge models with Nf >= 3 do not have finite-temperature transitions related to the breaking of the global nonabelian O(Nf) symmetry. However, in the zero-temperature limit they show a critical behavior characterized by a correlation length that increases exponentially with the inverse temperature, similarly to nonlinear O(N) sigma models. Their universal features are investigated by numerical finite-size scaling methods. The results show that the asymptotic low-temperature behavior belongs to the universality class of the two-dimensional RP(Nf-1) model.

preprint2020arXiv

Flux tubes and string breaking in three dimensional SU(2) Yang-Mills theory

We consider the three dimensional SU(2) Yang-Mills theory with adjoint static color sources, studying by lattice simulations how the shape of the flux tube changes when increasing the distance between them. The disappearance of the flux tube at string breaking is quite abrupt, but precursors of this phenomenon are present already when the separation between the sources is smaller than its critical value, a fact that influences also some details of the static potential.

preprint2020arXiv

Higher-charge three-dimensional compact lattice Abelian-Higgs models

We consider three-dimensional higher-charge multicomponent lattice Abelian-Higgs (AH) models, in which a compact U(1) gauge field is coupled to an N-component complex scalar field with integer charge q, so that they have local U(1) and global SU(N) symmetries. We discuss the dependence of the phase diagram, and the nature of the phase transitions, on the charge q of the scalar field and the number N>1 of components. We argue that the phase diagram of higher-charge models presents three different phases, related to the condensation of gauge-invariant bilinear scalar fields breaking the global SU(N) symmetry, and to the confinement/deconfinement of external charge-one particles. The transition lines separating the different phases show different features, which also depend on the number N of components. Therefore, the phase diagram of higher-charge models substantially differs from that of unit-charge models, which undergo only transitions driven by the breaking of the global SU(N) symmetry, while the gauge correlations do not play any relevant role. We support the conjectured scenario with numerical results, based on finite-size scaling analyses of Monte Carlo simuations for doubly-charged unit-length scalar fields with small and large number of components, i.e. N=2 and N=25.

preprint2020arXiv

Lattice Abelian-Higgs model with noncompact gauge fields

We consider a noncompact lattice formulation of the three-dimensional electrodynamics with $N$-component complex scalar fields, i.e., the lattice Abelian-Higgs model with noncompact gauge fields. For any $N\ge 2$, the phase diagram shows three phases differing for the behavior of the scalar-field and gauge-field correlations: the Coulomb phase (short-ranged scalar and long-ranged gauge correlations), the Higgs phase (condensed scalar-field and gapped gauge correlations), and the molecular phase (condensed scalar-field and long-ranged gauge correlations). They are separated by three transition lines meeting at a multicritical point. Their nature depends on the coexisting phases and on the number $N$ of components of the scalar field. In particular, the Coulomb-to-molecular transition line (where gauge correlations are irrelevant) is associated with the Landau-Ginzburg-Wilson $Φ^4$ theory sharing the same SU($N$) global symmetry but without explicit gauge fields. On the other hand, the Coulomb-to-Higgs transition line (where gauge correlations are relevant) turns out to be described by the continuum Abelian-Higgs field theory with explicit gauge fields. Our numerical study is based on finite-size scaling analyses of Monte Carlo simulations with $C^*$ boundary conditions (appropriate for lattice systems with noncompact gauge variables, unlike periodic boundary conditions), for several values of $N$, i.e., $N=2, 4, 10, 15$, and $25$. The numerical results agree with the renormalization-group predictions of the continuum field theories. In particular, the Coulomb-to-Higgs transitions are continuous for $N\gtrsim 10$, in agreement with the predictions of the Abelian-Higgs field theory.

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.

preprint2020arXiv

Three-dimensional lattice multiflavor scalar chromodynamics: interplay between global and gauge symmetries

We study the nature of the finite-temperature transition of the three-dimensional scalar chromodynamics with N_f flavors. These models are constructed by considering maximally O(M)-symmetric multicomponent scalar models, whose symmetry is partially gauged to obtain SU(N_c) gauge theories, with a residual nonabelian global symmetry given by U(N_f) for N_c>2 and Sp(N_f) for N_c=2, so that M = 2 N_c N_f. We find that their finite-temperature transition is continuous for N_f=2 and for all values of Nc we investigated, N_c=2,3,4. Such continuous transitions belong to universality classes related to the global symmetry group of the theory. For N_c=2 it belongs to the SO(5)=Sp(2)/Z_2 universality class, while for N_c>2 it belongs to the SO(3)=SU(2)/Z_2 universality class. For N_f>2, the transition is always of first order. These results match the predictions obtained by using the effective Landau-Ginzburg-Wilson approach in terms of a gauge-invariant order parameter. Our results indicate that the nonabelian gauge degrees of freedom are irrelevant at the transition. These conclusions are supported by an analysis of gauge-field dependent correlation functions, that are always short-ranged, even at the transition.

preprint2020arXiv

Three-dimensional phase transitions in multiflavor lattice scalar SO(Nc) gauge theories

We investigate the phase diagram and finite-temperature transitions of three-dimensional scalar SO(Nc) gauge theories with Nf scalar flavors. These models are constructed starting from a maximally O(N)-symmetric multicomponent scalar model (N = Nc Nf), whose symmetry is partially gauged to obtain an SO(Nc) gauge theory, with O(Nf) or U(Nf) global symmetry for Nc > 2 or Nc = 2, respectively. These systems undergo finite-temperature transitions, where the global symmetry is broken. Their nature is discussed using the Landau-Ginzburg-Wilson (LGW) approach, based on a gauge-invariant order parameter, and the continuum scalar SO(Nc) gauge theory. The LGW approach predicts that the transition is of first order for Nf > 2. For Nf = 2 the transition is predicted to be continuous: it belongs to the O(3) vector universality class for Nc=2 and to the XY universality class for any Nc > 2. We perform numerical simulations for Nc = 3 and Nf = 2,3. The numerical results are in agreement with the LGW predictions.

preprint2020arXiv

Universal low-temperature behavior of two-dimensional lattice scalar chromodynamics

We study the role that global and local nonabelian symmetries play in two-dimensional lattice gauge theories with multicomponent scalar fields. We start from a maximally O($M$)-symmetric multicomponent scalar model, Its symmetry is partially gauged to obtain an SU($N_c$) gauge theory (scalar chromodynamics) with global U$(N_f)$ (for $N_c\ge 3$) or Sp($N_f$) symmetry (for $N_c=2$), where $N_f>1$ is the number of flavors. Correspondingly, the fields belong to the coset $S^M$/SU($N_c$) where $S^M$ is the $M$-dimensional sphere and $M=2 N_f N_c$. In agreement with the Mermin-Wagner theorem, the system is always disordered at finite temperature and a critical behavior only develops in the zero-temperature limit. Its universal features are investigated by numerical finite-size scaling methods. The results show that the asymptotic low-temperature behavior belongs to the universality class of the 2D CP$^{N_f-1}$ field theory for $N_c>2$, and to that of the 2D Sp($N_f$) field theory for $N_c=2$. These universality classes correspond to 2D statistical field theories associated with symmetric spaces that are invariant under Sp($N_f$) transformations for $N_c=2$ and under SU($N_f$) for $N_c > 2$. These symmetry groups are the same invariance groups of scalar chromodynamics, apart from a U(1) flavor symmetry that is present for $N_f \ge N_c > 2$, which does not play any role in determining the asymptotic behavior of the model.

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

Two-dimensional multicomponent Abelian-Higgs lattice models

We study the two-dimensional lattice multicomponent Abelian-Higgs model, which is a lattice compact U(1) gauge theory coupled with an N-component complex scalar field, characterized by a global SU(N) symmetry. In agreement with the Mermin-Wagner theorem, the model has only a disordered phase at finite temperature and a critical behavior is only observed in the zero-temperature limit. The universal features are investigated by numerical analyses of the finite-size scaling behavior in the zero-temperature limit. The results show that the renormalization-group flow of the 2D lattice N-component Abelian-Higgs model is asymptotically controlled by the infinite gauge-coupling fixed point, associated with the universality class of the 2D CP(N-1) field theory.

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

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

Modeling of strain-induced Pockels effect in Silicon

We propose a theoretical model to describe the strain-induced linear electro-optic (Pockels) effect in centro-symmetric crystals. The general formulation is presented and the specific case of the strained silicon is investigated in detail because of its attractive properties for integrated optics. The outcome of this analysis is a linear relation between the second order susceptibility tensor and the strain gradient tensor, depending generically on fifteen coefficients. The proposed model greatly simplifies the description of the electro-optic effect in strained silicon waveguides, providing a powerful and effective tool for design and optimization of optical devices.

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

preprint2015arXiv

Topology and $θ$ dependence in finite temperature $G_2$ lattice gauge theory

In this work we study the topological properties of the $G_2$ lattice gauge theory by means of Monte Carlo simulations. We focus on the behaviour of topological quantities across the deconfinement transition and investigate observables related to the $θ$ dependence of the free energy. As in $SU(N)$ gauge theories, an abrupt change happens at deconfinement and an instanton gas behaviour rapidly sets in for $T>T_c$.

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

High-performance solution of the transport problem in a graphene armchair structure with a generic potential

We propose an efficient numerical method to study the transport properties of armchair graphene ribbons in the presence of a generic external potential. The method is based on a continuum envelope-function description with physical boundary conditions. The envelope functions are computed in the reciprocal space, and the transmission is then obtained with a recursive scattering matrix approach. This allows a significant reduction of the computational time with respect to finite difference simulations.

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

The Peierls argument for higher dimensional Ising models

The Peierls argument is a mathematically rigorous and intuitive method to show the presence of a non-vanishing spontaneous magnetization in some lattice models. This argument is typically explained for the $D=2$ Ising model in a way which cannot be easily generalized to higher dimension. The aim of this paper is to present an elementary discussion of the Peierls argument for the general $D$-dimensional Ising model.

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

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

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.

preprint2011arXiv

Armchair graphene nanoribbons: PT-symmetry breaking and exceptional points without dissipation

We consider a single layer graphene nanoribbon with armchair edges in a longitudinally constant external potential and point out that its transport properties can be described by means of an effective non-Hermitian Hamiltonian. We show that this system has some features typical of dissipative systems, namely the presence of exceptional points and of PT-symmetry breaking, although it is not dissipative.

preprint2011arXiv

Finite temperature effective string corrections in (3+1)D SU(2) lattice gauge theory

We study the effective string corrections to the inter-quark potential at finite temperature by simulating the SU(2) lattice gauge theory in four dimensions. We provide the first numerical evidence that the logarithmic correction to the potential, which was recently proposed to be a signature of the effective string at finite temperature, is universal also in (3+1)D gauge theory, thus extending previous results limited to the (2+1)D case.

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

Monopoles, abelian projection, and gauge invariance

A direct connection is proved between the Non-Abelian Bianchi Identities(NABI), and the abelian Bianchi identities for the 't Hooft tensor. As a consequence the existence of a non-zero magnetic current is related to the violation of the NABI's and is a gauge-invariant property. The construction allows to show that not all abelian projections can be used to expose monopoles in lattice configurations: each field configuration with non-zero magnetic charge identifies its natural projection, up to gauge transformations which tend to unity at large distances. It is shown that the so-called maximal-abelian gauge is a legitimate choice. It is also proved, starting from the NABI, that monopole condensation is a physical gauge invariant phenomenon, independent of the choice of the abelian projection.

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

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.