Source author record

Andrea Pelissetto

Andrea Pelissetto 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

73works
13topics
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

73 published item(s)

preprint2025arXiv

Out-of-equilibrium spinodal-like scaling behaviors across the magnetic first-order transitions of 2D and 3D Ising systems

We study the out-of-equilibrium scaling behavior of two-dimensional and three-dimensional Ising systems, when they are slowly driven across their {\em magnetic} first-order transitions at low temperature $T<T_c$, where $T_c$ is the temperature of their continuous transition. We consider Kibble-Zurek (KZ) protocols in which a spatially homogenous magnetic field $h$ varies as $h(t)=t/t_s$ with a time scale $t_s$. The KZ dynamics starts from negatively-magnetized configurations equilibrated at $h_i<0$ and stops at a positive value of $h$ where the configurations acquire a positive average magnetization. We consider the Metropolis and the heat-bath dynamics, which are two specific examples of a purely relaxational dynamics. We focus on two different dynamic regimes. We consider the out-equilibrium finite-size scaling (OFSS) limit in which the system size $L$ and the time scale $t_s$ diverge simultaneously, keeping an appropriate combination fixed. Then, we analyze the KZ dynamics in the thermodynamic limit (TL), obtained by taking first the $L\to\infty$ limit at fixed $t$ and $t_s$, and then considering the scaling behavior in the large-$t_s$ limit. Our numerical results provide evidence of OFSS, as predicted by general scaling arguments. The results in the TL show the emergence of spinodal-like behaviors: The passage from the negatively-magnetized phase to the positively-magnetized one occurs at positive values $h_*>0$ of the magnetic field, which decrease as $h_* \sim 1/(\ln t_s)^κ$, with $κ= 2$ and $κ=1$ in two and three dimensions, respectively, for $t_s\to\infty$. We identify $σ\equiv t (\ln t)^κ/t_s$ as the relevant scaling variable associated with the KZ dynamics in the TL.

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

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

Scaling properties of the dynamics at first-order quantum transitions when boundary conditions favor one of the two phases

We address the out-of-equilibrium dynamics of a many-body system when one of its Hamiltonian parameters is driven across a first-order quantum transition (FOQT). In particular, we consider systems subject to fixed boundary conditions, favoring one of the two phases separated by the FOQT: more precisely, boundary conditions that favor the same magnetized phase (EFBC) or opposite phases (OFBC) at the two ends of the chain. These issues are investigated within the paradigmatic one-dimensional quantum Ising model, in which FOQTs are driven by the longitudinal magnetic field h. We study the dynamic behavior for an instantaneous quench and for a protocol in which h is slowly varied across the FOQT. We develop a dynamic finite-size scaling theory for both EFBC and OFBC, which displays some remarkable differences with respect to the case of neutral boundary conditions. The corresponding relevant time scale shows a qualitative different size dependence in the two cases: it increases exponentially with the size in the case of EFBC, and as a power of the size in the case of OFBC.

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 monopole-free CP(N-1) models

We investigate the phase diagram, and the nature of the phase transitions, of three-dimensional monopole-free CP(N-1) models, characterized by a global U(N) symmetry and a U(1) gauge symmetry, and the absence of monopoles. We present numerical analyses based on Monte Carlo simulations for N=2,4,10,15, and 25. We observe a finite-temperature transition in all cases, related to the condensation of a local gauge-invariant order parameter. For N=2 we are unable to draw any definite conclusion on the nature of the transition. The results may be interpreted by either a very weak first-order transition or a continuous transition with anomalously large scaling corrections. However, the results allow us to exclude that the system develops the critical behavior of the O(3) vector universality class, as it occurs in the standard three-dimensional CP(1) model without monopole suppression. For N=4,10,15, the transition is of first order, and significantly weaker than that observed in the presence of monopoles. For N=25 the results are consistent with a conventional continuous transition. We compare our results with the existing literature and with the predictions of different field-theory approaches. Our results are consistent with the scenario in which the model undergoes continuous transitions for large values of N, in agreement with analytic large-N calculations.

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

Large-N behavior of three-dimensional lattice CP(N-1) models

We investigate the phase diagram and critical behavior of a three-dimensional lattice CP(N-1) model in the large-N limit. Numerical evidence of first-order transitions is always observed for sufficiently large values of N, i.e. N>2 up to N=100. The transition becomes stronger---both the latent heat and the surface tension increase---as N increases. Moreover, on the high-temperature side, gauge fields decorrelate on distances of the order of one lattice spacing for all values of N considered. Our results are consistent with a simple scenario, in which the transition is of first order for any N, including N=\infty. We critically discuss the analytic large-N calculations that predicted a large-N continuous transition, showing that one crucial assumption made in these computations fails for the model we consider.

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

Off-equilibrium scaling behaviors driven by time-dependent external fields in three-dimensional O(N) vector models

We consider the dynamical off-equilibrium behavior of the three-dimensional O$(N)$ vector model in the presence of a slowly-varying time-dependent spatially-uniform magnetic field ${\bm H}(t) = h(t)\,{\bm e}$, where ${\bm e}$ is a $N$-dimensional constant unit vector, $h(t)=t/t_s$, and $t_s$ is a time scale, at fixed temperature $T\le T_c$, where $T_c$ corresponds to the continuous order-disorder transition. The dynamic evolutions start from equilibrium configurations at $h_i < 0$, correspondingly $t_i < 0$, and end at time $t_f > 0$ with $h(t_f) > 0$, or vice versa. We show that the magnetization displays an off-equilibrium scaling behavior close to the transition line ${\bm H}(t)=0$. It arises from the interplay among the time $t$, the time scale $t_s$, and the finite size $L$. The scaling behavior can be parametrized in terms of the scaling variables $t_s^κ/L$ and $t/t_s^{κ_t}$, where $κ>0$ and $κ_t > 0$ are appropriate universal exponents, which differ at the critical point and for $T < T_c$. In the latter case, $κ$ and $κ_t$ also depend on the shape of the lattice and on the boundary conditions. We present numerical results for the Heisenberg ($N=3$) model under a purely relaxational dynamics. They confirm the predicted off-equilibrium scaling behaviors at and below $T_c$. We also discuss hysteresis phenomena in round-trip protocols for the time dependence of the external field. We define a scaling function for the hysteresis loop area of the magnetization that can be used to quantify how far the system is from equilibrium.

preprint2016arXiv

Phase diagram and multicritical behaviors of mixtures of 3D bosonic gases

We investigate the Bose-Einstein condensation patterns, the critical and multicritical behaviors of three-dimensional mixtures of bosonic gases with short-range density-density interactions. These systems have a global U(1)+U(1) symmetry, as the system Hamiltonian is invariant under independent U(1) transformations acting on each species. In particular, we consider the three-dimensional Bose-Hubbard model for two lattice bosonic gases coupled by an on-site inter-species density-density interaction. We study the phase diagram and the critical behaviors along the transition lines characterized by the Bose-Einstein condensation of one or both species. We present mean-field calculations and numerical finite-size scaling analyses of quantum Monte Carlo data. We also consider multicritical points, close to which it is possible to observe the condensation of both gas components. We determine the possible multicritical behaviors by using field-theoretical perturbative methods. We consider the U(1)+U(1)-symmetric Landau-Ginzburg-Wilson Phi4 theory and determine the corresponding stable fixed points of the renormalization-group flow. The analysis predicts that, in all cases, the multicritical behavior is analogous to the one that would be observed in systems of two identical gases, with an additional Z_2 exchange symmetry.

preprint2015arXiv

Bose-Einstein condensation and critical behavior of two-component bosonic gases

We study Bose-Einstein condensation (BEC) in three-dimensional two-component bosonic gases, characterizing the universal behaviors of the critical modes arising at the BEC transitions. For this purpose, we use field-theoretical (FT) renormalization-group (RG) methods and perform mean-field and numerical calculations. The FT RG analysis is based on the Landau-Ginzburg-Wilson Phi4 theory with two complex scalar fields which has the same symmetry as the bosonic system. In particular, for identical bosons with exchange Z_2,e symmetry, coupled by effective density-density interactions, the global symmetry is Z_2e X U(1) X U(1). At the BEC transition it may break into Z_2,e X Z_2 X Z_2 when both components condense simultaneously, or to U(1) X Z_2 when only one component condenses. This implies different universality classes for the corresponding critical behaviors. Numerical simulations of the two-component Bose-Hubbard model in the hard-core limit support the RG prediction: when both components condense simultaneously, the critical behavior is controlled by a decoupled XY fixed point, with unusual slowly-decaying scaling corrections arising from the on-site inter-species interaction.

preprint2015arXiv

Coarse-graining polymer solutions: a critical appraisal of single- and multi-site models

We critically discuss and review the general ideas behind single- and multi-site coarse-grained (CG) models as applied to macromolecular solutions in the dilute and semi-dilute regime. We first consider single-site models with zero-density and density-dependent pair potentials. We highlight advantages and limitations of each option in reproducing the thermodynamic behavior and the large-scale structure of the underlying reference model. As a case study we consider solutions of linear homopolymers in a solvent of variable quality. Secondly, we extend the discussion to multi-component systems presenting, as a test case, results for mixtures of colloids and polymers. Specifically, we found the CG model with zero-density potentials to be unable to predict fluid-fluid demixing in a reasonable range of densities for mixtures of colloids and polymers of equal size. For larger colloids, the polymer volume fractions at which phase separation occurs are largely overestimated. CG models with density-dependent potentials are somewhat less accurate than models with zero-density potentials in reproducing the thermodynamics of the system and, although they presents a phase separation, they significantly underestimate the polymer volume fractions along the binodal. Finally, we discuss a general multi-site strategy, which is thermodynamically consistent and fully transferable with the number of sites, and that allows us to overcome most of the limitations discussed for single-site models.

preprint2015arXiv

Critical mass renormalization in renormalized phi4 theories in two and three dimensions

We consider the O(N)-symmetric phi4 theory in two and three dimensions and determine the nonperturbative mass renormalization needed to obtain the phi4 continuum theory. The required nonperturbative information is obtained by resumming high-order perturbative series in the massive renormalization scheme, taking into account their Borel summability and the known large-order behavior of the coefficients. The results are in good agreement with those obtained in lattice calculations.

preprint2015arXiv

Integral-equation analysis of single-site coarse-grained models for polymer-colloid mixtures

We discuss the reliability of integral-equation methods based on several commonly used closure relations in determining the phase diagram of coarse-grained models of soft-matter systems characterized by mutually interacting soft and hard-core particles. Specifically, we consider a set of potentials appropriate to describe a system of hard-sphere colloids and linear homopolymers in good solvent, and investigate the behavior when the soft particles are smaller than the colloids, which is the regime of validity of the coarse-grained models. Using computer-simulation results as a benchmark, we find that the hypernetted-chain approximation provides accurate estimates of thermodynamics and structure in the colloid-gas phase in which the density of colloids is small. On the other hand, all closures considered appear to be unable to describe the behavior of the mixture in the colloid-liquid phase, as they cease to converge at polymer densities significantly smaller than those at the binodal. As a consequence, integral equations appear to be unable to predict a quantitatively correct phase diagram.

preprint2015arXiv

Off-equilibrium finite-size method for critical behavior analyses

We present a new dynamic off-equilibrium method for the study of continuous transitions, which represents a dynamic generalization of the usual equilibrium cumulant method. Its main advantage is that critical parameters are derived from numerical data obtained much before equilibrium has been attained. Therefore, the method is particularly useful for systems with long equilibration times, like spin glasses. We apply it to the three-dimensional Ising spin-glass model, obtaining accurate estimates of the critical exponents and of the critical temperature with a limited computational effort.

preprint2015arXiv

Phase diagram of mixtures of colloids and polymers in the thermal crossover from good to $θ$ solvent

We determine the phase diagram of mixtures of spherical colloids and neutral nonadsorbing polymers in the thermal crossover region between the $θ$ point and the good-solvent regime. We use the generalized free-volume theory (GFVT), which turns out to be quite accurate as long as $q = R_g/R_c\lesssim 1$ ($R_g$ is the radius of gyration of the polymer and $R_c$ is the colloid radius). Close to the $θ$ point the phase diagram is not very sensitive to solvent quality, while, close to the good-solvent region, changes of the solvent quality modify significantly the position of the critical point and of the binodals. We also analyze the phase behavior of aqueous solutions of charged colloids and polymers, using the extension of GFVT proposed by Fortini et al., J. Chem. Phys. 128, 024904 (2008).

preprint2015arXiv

Quantum Ising chains with boundary fields

We present a detailed study of the finite one-dimensional quantum Ising chain in a transverse field in the presence of boundary magnetic fields coupled with the order-parameter spin operator. We consider two magnetic fields located at the boundaries of the chain that have the same strength and that are aligned in the same or in the opposite direction. We derive analytic expressions for the gap in all phases for large values of the chain length L, as a function of the boundary field strength. We also investigate the behavior of the chain in the quantum ferromagnetic phase for oppositely aligned fields, focusing on the magnet-to-kink transition that occurs at a finite value of the magnetic field strength. At this transition we compute analytically the finite-size crossover functions for the gap, the magnetization profile, the two-point correlation function, and the density of fermionic modes. As the magnet-to-kink transition is equivalent to the wetting transition in two-dimensional classical Ising models, our results provide new analytic predictions for the finite-size behavior of Ising systems in a strip geometry at this transition.

preprint2015arXiv

Quantum transitions driven by one-bond defects in quantum Ising rings

We investigate quantum scaling phenomena driven by lower-dimensional defects in quantum Ising-like models. We consider quantum Ising rings in the presence of a bond defect. In the ordered phase, the system undergoes a quantum transition driven by the bond defect between a magnet phase, in which the gap decreases exponentially with increasing size, and a kink phase, in which the gap decreases instead with a power of the size. Close to the transition, the system shows a universal scaling behavior, which we characterize by computing, either analytically or numerically, scaling functions for the gap, the susceptibility, and the two-point correlation function. We discuss the implications of these results for the nonequilibrium dynamics in the presence of a slowly-varying parallel magnetic field h, when going across the first-order quantum transition at h=0.

preprint2015arXiv

Relevance of the axial anomaly at the finite-temperature chiral transition in QCD

We investigate the nature of the finite-temperature chiral transition in QCD with two light flavors, in the case of an effective suppression of the the U(1)_A symmetry breaking induced by the axial anomaly, which implies the symmetry breaking U(2)_L X U(2)_R -> U(2)_V, instead of SU(2)_L X SU(2)_R -> SU(2)_V. For this purpose, we perform a high-order field-theoretical perturbative study of the renormalization-group (RG) flow of the corresponding three-dimensional multiparameter Landau-Ginzburg-Wilson Phi4 theory with the same symmetry-breaking pattern. We confirm the existence of a stable fixed point (FP), and determine its attraction domain in the space of the bare quartic parameters. Therefore, the chiral QCD transition might be continuous also if the U(1)_A symmetry is effectively restored at Tc. However, the corresponding universality class differs from the O(4) vector universality class which would describe a continuous transition in the presence of a substantial U(1)_A symmetry breaking at Tc. We estimate the critical exponents of the U(2)_L X U(2)_R -> U(2)_V universality class by computing and analyzing their high-order perturbative expansions. These results are important to discriminate among the different scenarios for the scaling behavior of QCD with two light flavors close to the chiral transition.

preprint2015arXiv

Three-dimensional antiferromagnetic CP(N-1) models

We investigate the critical behavior of three-dimensional antiferromagnetic CP(N-1) [ACP(N-1)] models in cubic lattices, which are characterized by a global U(N) symmetry and a local U(1) gauge symmetry. Assuming that critical fluctuations are associated with a staggered gauge-invariant (hermitian traceless matrix) order parameter, we determine the corresponding Landau-Ginzburg-Wilson (LGW) model. For N=3 this mapping allows us to conclude that the three-component ACP(2) model undergoes a continuous transition that belongs to the O(8) vector universality class, with an effective enlargement of the symmetry at the critical point. This prediction is confirmed by a detailed numerical comparison of finite-size data for the ACP(2) and the O(8) vector models. We also present a renormalization-group (RG) analysis of the LGW theories for N>3. We compute perturbative series in two different renormalization schemes and analyze the corresponding RG flow. We do not find stable fixed points that can be associated with continuous transitions.

preprint2014arXiv

Accurate coarse-grained models for mixtures of colloids and linear polymers under good-solvent conditions

A coarse-graining strategy, previously developed for polymer solutions, is extended here to mixtures of linear polymers and hard-sphere colloids. In this approach groups of monomers are mapped onto a single pseudoatom (a blob) and the effective blob-blob interactions are obtained by requiring the model to reproduce some large-scale structural properties in the zero-density limit. We show that an accurate parametrization of the polymer-colloid interactions is obtained by simply introducing pair potentials between blobs and colloids. For the coarse-grained model in which polymers are modelled as four-blob chains (tetramers), the pair potentials are determined by means of the iterative Boltzmann inversion scheme, taking full-monomer pair correlation functions at zero-density as targets. For a larger number $n$ of blobs, pair potentials are determined by using a simple transferability assumption based on the polymer self-similarity. We validate the model by comparing its predictions with full-monomer results for the interfacial properties of polymer solutions in the presence of a single colloid and for thermodynamic and structural properties in the homogeneous phase at finite polymer and colloid density. The tetramer model is quite accurate for $q\lesssim 1$ ($q=\hat{R}_g/R_c$, where $\hat{R}_g$ is the zero-density polymer radius of gyration and $R_c$ is the colloid radius) and reasonably good also for $q=2$. For $q=2$ an accurate coarse-grained description is obtained by using the $n=10$ blob model. We also compare our results with those obtained by using single-blob models with state-dependent potentials.

preprint2014arXiv

Depletion effects in colloid-polymer solutions

The surface tension, the adsorption, and the depletion thickness of polymers close to a single nonadsorbing colloidal sphere are computed by means of Monte Carlo simulations. We consider polymers under good-solvent conditions and in the thermal crossover region between good-solvent and $θ$ behavior. In the dilute regime we consider a wide range of values of $q$, from $q = 0$ (planar surface) up to $q\approx 30$-50, while in the semidilute regime, for $ρ_p/ρ_p^*\le 4$ ($ρ_p$ is the polymer concentration and $ρ_p^*$ is its value at overlap), we only consider $q = 0,0.5,1$ and 2. The results are compared with the available theoretical predictions, verifying the existing scaling arguments. Field-theoretical results, both in the dilute and in the semidilute regime, are in good agreement with the numerical estimates for polymers under good-solvent conditions.

preprint2014arXiv

Finite-size scaling at first-order quantum transitions

We study finite-size effects at first-order quantum transitions (FOQTs). We show that the low-energy properties show a finite-size scaling (FSS) behavior, the relevant scaling variable being the ratio of the energy associated with the perturbation driving the transition and the finite-size energy gap at the FOQT point. The size dependence of the scaling variable is therefore essentially determined by the size dependence of the gap at the transition, which in turn depends on the boundary conditions. Our results have broad validity and, in particular, apply to any FOQT characterized by the degeneracy and crossing of the two lowest-energy states in the infinite-volume limit. In this case, a phenomenological two-level theory provides exact expressions for the scaling functions. Numerical results for the quantum Ising chain in transverse and parallel magnetic fields support the FSS ansatzes.

preprint2014arXiv

Finite-size scaling at quantum transitions

We develop the finite-size scaling (FSS) theory at quantum transitions, considering generic boundary conditions, such as open and periodic boundary conditions, and also the corrections to the leading FSS behaviors. Using renormalization-group (RG) theory, we generalize Wegner's scaling Ansatz to the quantum case, classifying the different sources of scaling corrections. We identify nonanalytic corrections due to irrelevant (bulk and boundary) RG perturbations and analytic contributions due to regular backgrounds and analytic expansions of the nonlinear scaling fields. To check the general predictions, we consider the quantum XY chain in a transverse field. For this model exact or numerically accurate results can be obtained by exploiting its fermionic quadratic representation. We study the FSS of several observables, such as the free energy, the energy differences between low-energy levels, correlation functions of the order parameter, etc., confirming the general predictions in all cases. Moreover, we consider bipartite entanglement entropies, which are characterized by the presence of additional scaling corrections, as predicted by conformal field theory.

preprint2014arXiv

Finite-size scaling at the first-order quantum transitions of quantum Potts chains

We investigate finite-size effects in quantum systems at first-order quantum transitions. For this purpose we consider the one-dimensional q-state Potts models which undergo a first-order quantum transition for any q>4, separating the quantum disordered and ordered phases with a discontinuity in the energy density of the ground state. The low-energy properties around the transition show finite-size scaling, described by general scaling ansatzes with respect to appropriate scaling variables. The size dependence of the scaling variables presents a particular sensitiveness to boundary conditions, which may be considered as a peculiar feature of first-order quantum transitions.

preprint2014arXiv

Scaling phenomena driven by inhomogeneous conditions at first-order quantum transitions

We investigate the effects of smooth inhomogeneities at first-order quantum transitions (FOQT), such as those arising from the presence of a space-dependent external field, which smooths out the typical discontinuities of the low-energy properties. We argue that scaling phenomena develop at the transition region where the external field takes the value corresponding to the FOQT of the homogenous system. We present numerical evidence of such scaling phenomena at the FOQTs of quantum Ising chains, driven by a parallel magnetic field when the system is in the ferromagnetic phase, and at the FOQT of the q-state Potts chain for q>4, driven by an even temperature-like parameter giving rise to a discontinuity of the ground-state energy density.

preprint2013arXiv

Comparing different coarse-grained potentials for star polymers

We compare different coarse-grained models for star polymers. We find that phenomenological models inspired by the Daoud-Cotton model reproduce quite poorly the thermodynamics of these systems, even if the potential is assumed to be density dependent, as done in the analysis of experimental results. We also determine the minumum value fc of the functionality of the star polymer for which a fluid-solid transition occurs. By applying the Hansen-Verlet criterion we find 35 < fc < 40. This result is confirmed by an analysis based on the modified (reference) hypernetted chain method and is qualitatively consistent with previous work.

preprint2013arXiv

Consistent coarse-graining strategy for polymer solutions in the thermal crossover from Good to Theta solvent

We extend our previously developed coarse-graining strategy for linear polymers with a tunable number n of effective atoms (blobs) per chain [D'Adamo et al., J. Chem. Phys. 137, 4901 (2012)] to polymer systems in thermal crossover between the good-solvent and the Theta regimes. We consider the thermal crossover in the region in which tricritical effects can be neglected, i.e. not too close to the Theta point, for a wide range of chain volume fractions Phi=c/c* (c* is the overlap concentration), up to Phi=30. Scaling crossover functions for global properties of the solution are obtained by Monte-Carlo simulations of the Domb-Joyce model. They provide the input data to develop a minimal coarse-grained model with four blobs per chain. As in the good-solvent case, the coarse-grained model potentials are derived at zero density, thus avoiding the inconsistencies related to the use of state-dependent potentials. We find that the coarse-grained model reproduces the properties of the underlying system up to some reduced density which increases when lowering the temperature towards the Theta state. Close to the lower-temperature crossover boundary, the tetramer model is accurate at least up to Phi<10, while near the good-solvent regime reasonably accurate results are obtained up to Phi<2. The density region in which the coarse-grained model is predictive can be enlarged by developing coarse-grained models with more blobs per chain. We extend the strategy used in the good-solvent case to the crossover regime. This requires a proper treatment of the length rescalings as before, but also a proper temperature redefinition as the number of blobs is increased. The case n=10 is investigated. Comparison with full-monomer results shows that the density region in which accurate predictions can be obtained is significantly wider than that corresponding to the n=4 case.

preprint2013arXiv

Predicting the thermodynamics by using state-dependent interactions

We reconsider the structure-based route to coarse graining in which the coarse-grained model is defined in such a way to reproduce some distributions functions of the original system as accurately as possible. We consider standard expressions for pressure and chemical potential applied to this family of coarse-grained models with density-dependent interactions and show that they only provide approximations to the pressure and chemical potential of the underlying original system. These approximations are then carefully compared in two cases: we consider a generic microscopic system in the low-density regime and polymer solutions under good-solvent conditions. Moreover, we show that the state-dependent potentials depend on the ensemble in which they have been derived. Therefore, care must be used in applying canonical state-dependent potentials to predict phase lines, which is typically performed in other ensembles.

preprint2012arXiv

Coarse-graining strategies in polymer solutions

We review a coarse-graining strategy (multiblob approach) for polymer solutions in which groups of monomers are mapped onto a single atom (a blob) and effective blob-blob interactions are obtained by requiring the coarse-grained model to reproduce some coarse-grained features of the zero-density isolated-chain structure. By tuning the level of coarse graining, i.e. the number of monomers to be mapped onto a single blob, the model should be adequate to explore the semidilute regime above the collapse transition, since in this case the monomer density is very small if chains are long enough. The implementation of these ideas has been previously based on a transferability hypothesis, which was not completely tested against full-monomer results (Pierleoni et al., J. Chem. Phys, 127, 171102 (2007)). We study different models proposed in the past and we compare their predictions to full-monomer results for the chain structure and the thermodynamics in the range of polymer volume fractions Φbetween 0 and 8. We find that the transferability assumption has a limited predictive power if a thermodynamically consistent model is required. We introduce a new tetramer model parametrized in such a way to reproduce not only zero-density intramolecular and intermolecular two-body probabilities, but also some intramolecular three-body and four-body distributions. We find that such a model correctly predicts three-chain effects, the structure and the thermodynamics up to Φ~ 2, a range considerably larger than that obtained with previous simpler models using zero-density potentials. Our results show the correctness of the ideas behind the multiblob approach but also that more work is needed to understand how to develop models with more effective monomers which would allow us to explore the semidilute regime at larger chain volume fractions.

preprint2012arXiv

Colloids and polymers in random colloidal matrices: demixing under good-solvent conditions

We consider a simplified coarse-grained model for colloid-polymer mixtures, in which polymers are represented as monoatomic molecules interacting by means of pair potentials. We use it to study polymer-colloid segregation in the presence of a quenched matrix of colloidal hard spheres. We fix the polymer-to-colloid size ratio to 0.8 and consider matrices such that the fraction f of the volume that is not accessible to the colloids due to the matrix is equal to 40%. As in the Asakura-Oosawa-Vrij (AOV) case, we find that binodal curves in the polymer and colloid volume-fraction plane have a small dependence on disorder. As for the position of the critical point, the behavior is different from that observed in the AOV case: while the critical colloid volume fraction is essentially the same in the bulk and in the presence of the matrix, the polymer volume fraction at criticality increases as f increases. At variance with the AOV case, no capillary colloid condensation or evaporation is generically observed.

preprint2012arXiv

Consistent and transferrable coarse-grained model for semidilute polymer solutions in good solvent

We present a coarse-grained model for linear polymers with a tunable number of effective atoms (blobs) per chain interacting by intra- and inter-molecular potentials obtained at zero density. We show how this model is able to accurately reproduce the universal properties of the underlying solution of athermal linear chains at various levels of coarse-graining and in a range of chain densities which can be widened by increasing the spatial resolution of the multiblob representation, i.e., the number of blobs per chain. The present model is unique in its ability to quantitatively predict thermodynamic and large scale structural properties of polymer solutions deep in the semidilute regime with a very limited computational effort, overcoming most of the problems related to the simulations of semidilute polymer solutions in good solvent conditions.

preprint2012arXiv

Polymers as compressible soft spheres

We consider a coarse-grained model in which polymers under good-solvent conditions are represented by soft spheres whose radii, which should be identified with the polymer radii of gyrations, are allowed to fluctuate. The corresponding pair potential depends on the sphere radii. This model is a single-sphere version of the one proposed in Vettorel et al., Soft Matter 6, 2282 (2010), and it is sufficiently simple to allow us to determine all potentials accurately from full-monomer simulations of two isolated polymers (zero-density potentials). We find that in the dilute regime (which is the expected validity range of single-sphere coarse-grained models based on zero-density potentials) this model correctly reproduces the density dependence of the radius of gyration. However, for the thermodynamics and the intermolecular structure, the model is largely equivalent to the simpler one in which the sphere radii are fixed to the average value of the radius of gyration and radiiindependent potentials are used: for the thermodynamics there is no advantage in considering a fluctuating sphere size.

preprint2012arXiv

Renormalization-group flow and asymptotic behaviors at the Berezinskii-Kosterlitz-Thouless transitions

We investigate the general features of the renormalization-group flow at the Berezinskii-Kosterlitz-Thouless (BKT) transition, providing a thorough quantitative description of the asymptotc critical behavior, including the multiplicative and subleading logarithmic corrections. For this purpose, we consider the RG flow of the sine-Gordon model around the renormalizable point which describes the BKT transition. We reduce the corresponding beta-functions to a universal canonical form, valid to all perturbative orders. Then, we determine the asymptotic solutions of the RG equations in various critical regimes: the infinite-volume critical behavior in the disordered phase, the finite-size scaling limit for homogeneous systems of finite size, and the trap-size scaling limit occurring in 2D bosonic particle systems trapped by an external space-dependent potential.

preprint2011arXiv

Ferromagnetic-glassy transitions in three-dimensional Ising spin glasses

We investigate the ferromagnetic-glassy transitions which separate the low-temperature ferromagnetic and spin-glass phases in the temperature-disorder phase diagram of three-dimensional Ising spin-glass models. For this purpose, we consider the cubic-lattice +-J (Edwards-Anderson) Ising model with bond distribution $P(J) = p δ(J - 1) + (1-p) δ(J + 1)$, and present a numerical Monte Carlo study of the critical behavior along the line that marks the onset of ferromagnetism. The finite-size scaling analysis of the Monte Carlo data shows that the ferromagnetic-glassy transition line is slightly reentrant. As a consequence, for an interval of the disorder parameter p, around p=0.77, the system presents a low-temperature glassy phase, an intermediate ferromagnetic phase, and a high-temperature paramagnetic phase. Along the ferromagnetic-glassy transition line magnetic correlations show a universal critical behavior with critical exponents nu=0.96(2) and eta=-0.39(2). The hyperscaling relation beta/nu = (1 + eta)/2 is satisfied at the transitions, so that beta/nu = 0.305(10). This magnetic critical behavior represents a new universality class for ferromagnetic transitions in Ising-like disordered systems. Overlap correlations are apparently not critical and show a smooth behavior across the transition.

preprint2011arXiv

Finite-size scaling in two-dimensional Ising spin glass models

We study the finite-size behavior of two-dimensional spin-glass models. We consider the +-J model for two different values of the probability of the antiferromagnetic bonds and the model with Gaussian distributed couplings. The analysis of renormalization-group invariant quantities, the overlap susceptibility, and the two-point correlation function confirms that they belong to the same universality class. We analyze in detail the standard finite-size scaling limit in terms of TL^(1/nu) in the +-J model. We find that it holds asymptotically. This result is consistent with the low-temperature crossover scenario in which the crossover temperature, which separates the universal high-temperature region from the discrete low-temperature regime, scales as T_c(L) ~ L^(-theta_S) with theta_S \approx 0.5.

preprint2011arXiv

Fluid-fluid demixing curves for colloid-polymer mixtures in a random colloidal matrix

We study fluid-fluid phase separation in a colloid-polymer mixture adsorbed in a colloidal porous matrix close to the θ-point. For this purpose we consider the Asakura-Oosawa model in the presence of a quenched matrix of colloidal hard spheres. We study the dependence of the demixing curve on the parameters that characterize the quenched matrix, fixing the polymer-to-colloid size ratio to 0.8. We find that, to a large extent, demixing curves depend only on a single parameter f, which represents the volume fraction which is unavailable to the colloids. We perform Monte Carlo simulations for volume fractions f equal to 40% and 70%, finding that the binodal curves in the polymer and colloid packing-fraction plane have a small dependence on disorder. The critical point instead changes significantly: for instance, the colloid packing fraction at criticality increases with increasing f. Finally, we observe for some values of the parameters capillary condensation of the colloids: a bulk colloid-poor phase is in chemical equilibrium with a colloid-rich phase in the matrix.

preprint2010arXiv

Geometrical Properties of Two-Dimensional Interacting Self-Avoiding Walks at the Theta-Point

We perform a Monte Carlo simulation of two-dimensional N-step interacting self-avoiding walks at the theta point, with lengths up to N=3200. We compute the critical exponents, verifying the Coulomb-gas predictions, the theta-point temperature T_theta = 1.4986(11), and several invariant size ratios. Then, we focus on the geometrical features of the walks, computing the instantaneous shape ratios, the average asphericity, and the end-to-end distribution function. For the latter quantity, we verify in detail the theoretical predictions for its small- and large-distance behavior.

preprint2010arXiv

Magnetic and glassy transitions in the square-lattice XY model with random phase shifts

We investigate the magnetic and glassy transitions of the square-lattice XY model in the presence of random phase shifts. We consider two different random-shift distributions: the Gaussian distribution and a slightly different distribution (cosine distribution) which allows the exact determination of the Nishimori line where magnetic and overlap correlation functions are equal. We perform Monte Carlo simulations for several values of the temperature and of the variance of the disorder distribution, in the paramagnetic phase close to the magnetic and glassy transition lines. We find that, along the transition line separating the paramagnetic and the quasi-long-range order phases, magnetic correlation functions show a universal Kosterlitz-Thouless behavior as in the pure XY model, while overlap correlations show a disorder-dependent critical behavior. This behavior is observed up to a multicritical point which, in the cosine model, lies on the Nishimori line. Finally, for large values of the disorder variance, we observe a universal zero-temperature glassy critical transition, which is in the same universality class as that occurring in the gauge-glass model.

preprint2010arXiv

Universality of the glassy transitions in the two-dimensional +- J Ising model

We investigate the zero-temperature glassy transitions in the square-lattice +- J Ising model, with bond distribution $P(J_{xy}) = p δ(J_{xy} - J) + (1-p) δ(J_{xy} + J)$; p=1 and p=1/2 correspond to the pure Ising model and to the Ising spin glass with symmetric bimodal distribution, respectively. We present finite-temperature Monte Carlo simulations at p=4/5, which is close to the low-temperature paramagnetic-ferromagnetic transition line located at p=0.89, and at p=1/2. Their comparison provides a strong evidence that the glassy critical behavior that occurs for $1-p_0<p<p_0$, $p_0=0.897$, is universal, i.e., independent of p. Moreover, we show that glassy and magnetic modes are not coupled at the multicritical zero-temperature point where the paramagnetic-ferromagnetic transition line and the T=0 glassy transition line meet. On the theoretical side we discuss the validity of finite-size scaling in glassy systems with a zero-temperature transition and a discrete Hamiltonian spectrum. Because of a freezing phenomenon which occurs in a finite volume at sufficiently low temperatures, the standard finite-size scaling limit in terms of $TL^{1/ν}$ does not exist: the renormalization-group invariant quantity $ξ/L$ should be used instead as basic variable.

preprint2009arXiv

Quasi-long-range order in the 2D XY model with random phase shifts

We study the square-lattice XY model in the presence of random phase shifts. We consider two different disorder distributions with zero average shift and investigate the low-temperature quasi-long-range order phase which occurs for sufficiently low disorder. By means of Monte Carlo simulations we determine several universal quantities which are then compared with the analytic predictions of the random spin-wave theory. We observe a very good agreement which indicates that the universal long-distance behavior in the whole low-disorder low-temperature phase is fully described by the random spin-wave theory.

preprint2005arXiv

Correction-to-scaling exponents for two-dimensional self-avoiding walks

We study the correction-to-scaling exponents for the two-dimensional self-avoiding walk, using a combination of series-extrapolation and Monte Carlo methods. We enumerate all self-avoiding walks up to 59 steps on the square lattice, and up to 40 steps on the triangular lattice, measuring the mean-square end-to-end distance, the mean-square radius of gyration and the mean-square distance of a monomer from the endpoints. The complete endpoint distribution is also calculated for self-avoiding walks up to 32 steps (square) and up to 22 steps (triangular). We also generate self-avoiding walks on the square lattice by Monte Carlo, using the pivot algorithm, obtaining the mean-square radii to ~0.01% accuracy up to N = 4000. We give compelling evidence that the first non-analytic correction term for two-dimensional self-avoiding walks is Delta_1 = 3/2. We compute several moments of the endpoint distribution function, finding good agreement with the field-theoretic predictions. Finally, we study a particular invariant ratio that can be shown, by conformal-field-theory arguments, to vanish asymptotically, and we find the cancellation of the leading analytic correction.

preprint2005arXiv

Critical behavior of two-dimensional fully frustrated XY systems

We study the phase diagram of the two-dimensional fully frustrated XY model (FFXY) and of two related models, a lattice discretization of the Landau-Ginzburg-Wilson Hamiltonian for the critical modes of the FFXY model, and a coupled Ising-XY model. We present Monte Carlo simulations on square lattices $L\times L$, $L\lesssim 10^3$. We show that the low-temperature phase of these models is controlled by the same line of Gaussian fixed points as in the standard XY model. We find that, if a model undergoes a unique transition by varying temperature, then the transition is of first order. In the opposite case we observe two very close transitions: a transition associated with the spin degrees of freedom and, as temperature increases, a transition where chiral modes become critical. If they are continuous, they belong to the Kosterlitz-Thouless and to the Ising universality class, respectively. Ising and Kosterlitz-Thouless behavior is observed only after a preasymptotic regime, which is universal to some extent. In the chiral case, the approach is nonmonotonic for most observables, and there is a wide region in which finite-size scaling is controlled by an effective exponent $ν_{\rm eff} \approx 0.8$. This explains the result $ν\approx 0.8$ of many previous studies using smaller lattices.

preprint2005arXiv

Transitions and crossover phenomena in fully frustrated XY systems

We study the two-dimensional fully frustrated XY (FFXY) model and two related models, a discretization of the Landau-Ginzburg-Wilson Hamiltonian for the critical modes of the FFXY model and a coupled Ising-XY model, by means of Monte Carlo simulations on square lattices L x L, L=O(10^3). We show that their phase diagram is characterized by two very close chiral and spin transitions, at T_ch > T_sp respectively, of the Ising and Kosterlitz-Thouless type. At T_ch the Ising regime sets in only after a preasymptotic regime, which appears universal to some extent. The approach is nonmonotonic for most observables, with a wide region controlled by an effective exponent nu_eff=0.8.

preprint2004arXiv

Critical behavior of O(2)xO(N) symmetric models

We investigate the controversial issue of the existence of universality classes describing critical phenomena in three-dimensional statistical systems characterized by a matrix order parameter with symmetry O(2)xO(N) and symmetry-breaking pattern O(2)xO(N) -> O(2)xO(N-2). Physical realizations of these systems are, for example, frustrated spin models with noncollinear order. Starting from the field-theoretical Landau-Ginzburg-Wilson Hamiltonian, we consider the massless critical theory and the minimal-subtraction scheme without epsilon expansion. The three-dimensional analysis of the corresponding five-loop expansions shows the existence of a stable fixed point for N=2 and N=3, confirming recent field-theoretical results based on a six-loop expansion in the alternative zero-momentum renormalization scheme defined in the massive disordered phase. In addition, we report numerical Monte Carlo simulations of a class of three-dimensional O(2)xO(2)-symmetric lattice models. The results provide further support to the existence of the O(2)xO(2) universality class predicted by the field-theoretical analyses.

preprint2004arXiv

Crossover behavior in three-dimensional dilute spin systems

We study the crossover behaviors that can be observed in the high-temperature phase of three-dimensional dilute spin systems, using a field-theoretical approach. In particular, for randomly dilute Ising systems we consider the Gaussian-to-random and the pure-Ising-to-random crossover, determining the corresponding crossover functions for the magnetic susceptibility and the correlation length. Moreover, for the physically interesting cases of dilute Ising, XY, and Heisenberg systems, we estimate several universal ratios of scaling-correction amplitudes entering the high-temperature Wegner expansion of the magnetic susceptibility, of the correlation length, and of the zero-momentum quartic couplings.

preprint2004arXiv

Multicritical behavior in frustrated spin systems with noncollinear order

We investigate the phase diagram and, in particular, the nature of the the multicritical point in three-dimensional frustrated $N$-component spin models with noncollinear order in the presence of an external field, for instance easy-axis stacked triangular antiferromagnets in the presence of a magnetic field along the easy axis. For this purpose we study the renormalization-group flow in a Landau-Ginzburg-Wilson ϕ^4 theory with symmetry O(2)x[Z_2 +O(N-1)] that is expected to describe the multicritical behavior. We compute its MS βfunctions to five loops. For N\ge 4, their analysis does not support the hypothesis of an effective enlargement of the symmetry at the multicritical point, from O(2) x [Z_2+O(N-1)] to O(2)xO(N). For the physically interesting case N=3, the analysis does not allow us to exclude the corresponding symmetry enlargement controlled by the O(2)xO(3) fixed point. Moreover, it does not provide evidence for any other stable fixed point. Thus, on the basis of our field-theoretical results, the transition at the multicritical point is expected to be either continuous and controlled by the O(2)xO(3) fixed point or to be of first order.

preprint2004arXiv

Two-Dimensional Heisenberg Model with Nonlinear Interactions: 1/N Corrections

We investigate a two-dimensional classical $-vector model with a generic nearest-neighbor interaction $W(\bsigma_i\cdot \bsigma_j)$ in the large-N limit, focusing on the finite-temperature transition point at which energy-energy correlations become critical. We show that this transition belongs to the Ising universality class. However, the width of the region in which Ising behavior is observed scales as $1/N^{3/2}$ along the magnetic direction and as 1/N in the thermal direction; outside a crossover to mean-field behavior occurs. This explains why only mean-field behavior is observed for $N=\infty$

preprint2003arXiv

Comment on ``Dynamic behavior of anisotropic non-equilibrium driving lattice gases''

In a recent Letter Albano and Saracco study the dynamic critical behavior of some anisotropic driven lattice gases by Monte Carlo (MC) simulations. In this Comment we point out that the Ansätze they use to relate the measured scaling exponents with the critical exponents analytically computed within different field-theoretical approaches do not take properly into account the strongly anisotropic nature of the phase transition, by implicitly assuming $z = z_{\bot} = z_{\parallel}$. As a consequence, at variance with the claims by the authors, their MC data are not conclusive to determine which one of the field theories proposed in the literature correctly describes the universal properties of the phase transition in these lattice gases.

preprint2003arXiv

On the nature of the finite-temperature transition in QCD

We discuss the nature of the finite-temperature transition in QCD with N_f massless flavors. Universality arguments show that a continuous (second-order) transition must be related to a 3-D universality class characterized by a complex N_f X N_f matrix order parameter and by the symmetry-breaking pattern [SU(N_f)_L X SU(N_f)_R]/Z(N_f)_V -> SU(N_f)_V/Z(N_f)_V, or [U(N_f)_L X U(N_f)_R]/U(1)_V -> U(N_f)_V/U(1)_V if the U(1)_A symmetry is effectively restored at T_c. The existence of any of these universality classes requires the presence of a stable fixed point in the corresponding 3-D Phi^4 theory with the expected symmetry-breaking pattern. Otherwise, the transition is of first order. In order to search for stable fixed points in these Phi^4 theories, we exploit a 3-D perturbative approach in which physical quantities are expanded in powers of appropriate renormalized quartic couplings. We compute the corresponding Callan-Symanzik beta-functions to six loops. We also determine the large-order behavior to further constrain the analysis. No stable fixed point is found, except for N_f=2, corresponding to the symmetry-breaking pattern [SU(2)_L X SU(2)_R]/Z(2)_V -> SU(2)_V/Z(2)_V equivalent to O(4) -> O(3). Our results confirm and put on a firmer ground earlier analyses performed close to four dimensions, based on first-order calculations in the framework of the epsilon=4-d expansion. These results indicate that the finite-temperature phase transition in QCD is of first order for N_f>2. A continuous transition is allowed only for N_f=2. But, since the theory with symmetry-breaking pattern [U(2)_L X U(2)_R]/U(1)_V -> U(2)_V/U(1)_V does not have stable fixed points, the transition can be continuous only if the effective breaking of the U(1)_A symmetry is sufficiently large.

preprint2003arXiv

Spin models with random anisotropy and reflection symmetry

We study the critical behavior of a general class of cubic-symmetric spin systems in which disorder preserves the reflection symmetry $s_a\to -s_a$, $s_b\to s_b$ for $b\not= a$. This includes spin models in the presence of random cubic-symmetric anisotropy with probability distribution vanishing outside the lattice axes. Using nonperturbative arguments we show the existence of a stable fixed point corresponding to the random-exchange Ising universality class. The field-theoretical renormalization-group flow is investigated in the framework of a fixed-dimension expansion in powers of appropriate quartic couplings, computing the corresponding $β$-functions to five loops. This analysis shows that the random Ising fixed point is the only stable fixed point that is accessible from the relevant parameter region. Therefore, if the system undergoes a continuous transition, it belongs to the random-exchange Ising universality class. The approach to the asymptotic critical behavior is controlled by scaling corrections with exponent $Δ= - α_r$, where $α_r\simeq -0.05$ is the specific-heat exponent of the random-exchange Ising model.

preprint2003arXiv

The 3-D O(4) universality class and the phase transition in two-flavor QCD

We determine the critical equation of state of the three-dimensional O(4) universality class. We first consider the small-field expansion of the effective potential (Helmholtz free energy). Then, we apply a systematic approximation scheme based on polynomial parametric representations that are valid in the whole critical regime, satisfy the correct analytic properties (Griffiths' analyticity), take into account the Goldstone singularities at the coexistence curve, and match the small-field expansion of the effective potential. From the approximate representations of the equation of state, we obtain estimates of several universal amplitude ratios. The three-dimensional O(4) universality class is expected to describe the finite-temperature chiral transition of quantum chromodynamics with two light flavors. Within this picture, the O(4) critical equation of state relates the reduced temperature, the quark masses, and the condensates around T_c in the limit of vanishing quark masses.

preprint2002arXiv

Two-Dimensional Heisenberg Model with Nonlinear Interactions

We investigate a two-dimensional classical $N$-vector model with a nonlinear interaction $(1 + \bsigma_i\cdot \bsigma_j)^p$ in the large-N limit. As observed for N=3 by Blöte {\em et al.} [Phys. Rev. Lett. {\bf 88}, 047203 (2002)], we find a first-order transition for $p>p_c$ and no finite-temperature phase transitions for $p < p_c$. For $p>p_c$, both phases have short-range order, the correlation length showing a finite discontinuity at the transition. For $p=p_c$, there is a peculiar transition, where the spin-spin correlation length is finite while the energy-energy correlation length diverges.

preprint2001arXiv

Critical structure factors of bilinear fields in O(N)-vector models

We compute the two-point correlation functions of general quadratic operators in the high-temperature phase of the three-dimensional O(N) vector model by using field-theoretical methods. In particular, we study the small- and large-momentum behavior of the corresponding scaling functions, and give general interpolation formulae based on a dispersive approach. Moreover, we determine the crossover exponent $ϕ_T$ associated with the traceless tensorial quadratic field, by computing and analyzing its six-loop perturbative expansion in fixed dimension. We find: $ϕ_T=1.184(12)$, $ϕ_T=1.271(21)$, and $ϕ_T=1.40(4)$ for $N=2,3,5$ respectively.

preprint2000arXiv

Finite-Size Correlation Length and Violations of Finite-Size Scaling

We address the problem of the definition of the finite-volume correlation length. First, we study the large-N limit of the N-vector model, and we show the existence of several constraints on the definition if regularity of the finite-size scaling functions and correct anomalous behaviour above the upper critical dimension are required. Then, we study in detail a model in which the zero mode in prohibited. Such a model is a generalization of the fixed-magnetization Ising model which is equivalent to the lattice gas. Also in this case, we find that the finite-volume correlation length must satisfy appropriate constraints in order to obtain regular finite-size scaling functions, and, above the upper critical dimension, an anomalous scaling behaviour. The large-N results are confirmed by a one-loop calculation in the lattice $ϕ^4$ theory.

preprint2000arXiv

Operator Product Expansion on the Lattice: a Numerical Test in the Two-Dimensional Non-Linear Sigma-Model

We consider the short-distance behaviour of the product of the Noether O(N) currents in the lattice nonlinear sigma-model. We compare the numerical results with the predictions of the operator product expansion, using one-loop perturbative renormalization-group improved Wilson coefficients. We find that, even on quite small lattices (m a \approx 1/6), the perturbative operator product expansion describes that data with an error of 5-10% in a large window 2a \ltapprox x \ltapprox m^{-1}. We present a detailed discussion of the possible systematic errors.

preprint2000arXiv

The critical equation of state of three-dimensional XY systems

We address the problem of determining the critical equation of state of three-dimensional XY systems. For this purpose we first consider the small-field expansion of the effective potential (Helmholtz free energy) in the high-temperature phase. We compute the first few nontrivial zero-momentum n-point renormalized couplings, which parametrize such expansion, by analyzing the high-temperature expansion of an improved lattice Hamiltonian with suppressed leading scaling corrections. These results are then used to construct parametric representations of the critical equation of state which are valid in the whole critical regime, satisfy the correct analytic properties (Griffith's analyticity), and take into account the Goldstone singularities at the coexistence curve. A systematic approximation scheme is introduced, which is limited essentially by the number of known terms in the small-field expansion of the effective potential. From our approximate representations of the equation of state, we derive estimates of universal ratios of amplitudes. For the specific-heat amplitude ratio we obtain A^+/A^-=1.055(3), to be compared with the best experimental estimate A^+/A^-=1.054(1).

preprint1998arXiv

Corrections to Finite-Size Scaling in the Lattice N-Vector Model for Infinite N

We compute the corrections to finite-size scaling for the N-vector model on the square lattice in the large-N limit. We find that corrections behave as log L/L^2. For tree-level improved hamiltonians corrections behave as 1/L^2. In general l-loop improvement is expected to reduce this behaviour to 1/(L^2 \log^l L). We show that the finite-size-scaling and the perturbative limit do not commute in the calculation of the corrections to finite-size scaling. We present also a detailed study of the corrections for the RP^N-model.

preprint1996arXiv

Multi-Grid Monte Carlo via $XY$ Embedding I. General Theory and Two-Dimensional $O(N)$-Symmetric Nonlinear $σ$-Models

We introduce a variant of the multi-grid Monte Carlo (MGMC) method, based on the embedding of an $XY$ model into the target model, and we study its mathematical properties for a variety of nonlinear $σ$-models. We then apply the method to the two-dimensional $O(N)$-symmetric nonlinear $σ$-models (also called $N$-vector models) with $N=3,4,8$ and study its dynamic critical behavior. Using lattices up to $256 \times 256$, we find dynamic critical exponents $z_{int,{\cal M}^2} \approx 0.70 \pm 0.08$, $0.60 \pm 0.07$, $0.52 \pm 0.10$ for $N=3,4,8$, respectively (subjective 68\% confidence intervals). Thus, for these asymptotically free models, critical slowing-down is greatly reduced compared to local algorithms, but not completely eliminated; and the dynamic critical exponent does apparently vary with $N$. We also analyze the static data for $N=8$ using a finite-size-scaling extrapolation method. The correlation length $ξ$ agrees with the four-loop asymptotic-freedom prediction to within $\approx 1\%$ over the interval $12 < ξ< 650$.

preprint1996arXiv

Multi-Grid Monte Carlo via $XY$ Embedding. II. Two-Dimensional $SU(3)$ Principal Chiral Model

We carry out a high-precision simulation of the two-dimensional $SU(3)$ principal chiral model at correlation lengths $ξ$ up to $\sim 4 \times 10^5$, using a multi-grid Monte Carlo (MGMC) algorithm and approximately one year of Cray C-90 CPU time. We extrapolate the finite-volume Monte Carlo data to infinite volume using finite-size-scaling theory, and we discuss carefully the systematic and statistical errors in this extrapolation. We then compare the extrapolated data to the renormalization-group predictions. The deviation from asymptotic scaling, which is $\approx 12%$ at $ξ\sim 25$, decreases to $\approx 2%$ at $ξ\sim 4 \times 10^5$. We also analyze the dynamic critical behavior of the MGMC algorithm using lattices up to $256 \times 256$, finding the dynamic critical exponent $z_{int,{\cal M}^2} \approx 0.45 \pm 0.02$ (subjective 68% confidence interval). Thus, for this asymptotically free model, critical slowing-down is greatly reduced compared to local algorithms, but not completely eliminated.

preprint1995arXiv

Continuum Limits and Exact Finite-Size-Scaling Functions for One-Dimensional $O(N)$-Invariant Spin Models

We solve exactly the general one-dimensional $O(N)$-invariant spin model taking values in the sphere $S^{N-1}$, with nearest-neighbor interactions, in finite volume with periodic boundary conditions, by an expansion in hyperspherical harmonics. The possible continuum limits are discussed for a general one-parameter family of interactions, and an infinite number of universality classes is found. For these classes we compute the finite-size-scaling functions and the leading corrections to finite-size scaling. A special two-parameter family of interactions (which includes the mixed isovector/isotensor model) is also treated, and no additional universality classes appear. In the appendices we give new formulae for the Clebsch-Gordan coefficients and 6--$j$ symbols of the $O(N)$ group, and some new generalizations of the Poisson summation formula; these may be of independent interest.

preprint1993arXiv

A General Limitation on Monte Carlo Algorithms of Metropolis Type

We prove that for any Monte Carlo algorithm of Metropolis type, the autocorrelation time of a suitable ``energy''-like observable is bounded below by a multiple of the corresponding ``specific heat''. This bound does not depend on whether the proposed moves are local or non-local; it depends only on the distance between the desired probability distribution $π$ and the probability distribution $π^{(0)}$ for which the proposal matrix satisfies detailed balance. We show, with several examples, that this result is particularly powerful when applied to non-local algorithms.