Monte Carlo Investigations of Phase Transitions: Status and Perspectives Kurt Bindera , Erik Luijtena; ∗, Marcus Muller a , Nigel B

Total Page:16

File Type:pdf, Size:1020Kb

Monte Carlo Investigations of Phase Transitions: Status and Perspectives Kurt Bindera , Erik Luijtena; ∗, Marcus Muller a , Nigel B Physica A 281 (2000) 112–128 www.elsevier.com/locate/physa Monte Carlo investigations of phase transitions: status and perspectives Kurt Bindera , Erik Luijtena; ∗, Marcus Muller a , Nigel B. Wildinga; 1, Henk W.J. Bloteb aInstitut fur Physik, Johannes Gutenberg-Universitat, Staudinger Weg 7, D-55099 Mainz, Germany bDepartment of Physics, Delft University of Technology, P.O. Box 5046, 2600 GA Delft, The Netherlands Abstract Using the concept of ÿnite-size scaling, Monte Carlo calculations of various models have become a very useful tool for the study of critical phenomena, with the system linear dimension as a variable. As an example, several recent studies of Ising models are discussed, as well as the extension to models of polymer mixtures and solutions. It is shown that using appropriate cluster algorithms, even the scaling functions describing the crossover from the Ising universality class to the mean-ÿeld behavior with increasing interaction range can be described. Additionally, the issue of ÿnite-size scaling in Ising models above the marginal dimension (d∗ = 4) is discussed. c 2000 Elsevier Science B.V. All rights reserved. Keywords: Critical phenomena; Ising model; Crossover scaling; Polymers; Finite-size scaling 1. Introduction It is a common belief that at the present time, about 30 years after the renormalization- group theory of critical phenomena was invented [1], static critical behavior of sys- tems in thermal equilibrium is rather well understood. In particular, this is expected to be true for the most intensively studied case, the Ising universality class [2,3], to which systems such as uniaxial ferromagnets, binary alloys, simple uids, uid mix- tures, polymer solutions and polymer blends belong [4]. However, in the present work, we shall draw attention to some aspects of critical behavior in Ising-like spin systems which are, even today, still incompletely understood. The ÿrst of these concerns the ∗ Correspondence address: Max-Planck-Institut fur Polymerforschung, Postfach 3148, 55021 Mainz, Germany. E-mail address: [email protected] (E. Luijten) 1 Present address: Department of Physics and Astronomy, The University of Edinburgh, Edinburgh EH9 3JZ, UK. 0378-4371/00/$ - see front matter c 2000 Elsevier Science B.V. All rights reserved. PII: S 0378-4371(00)00025-X K. Binder et al. / Physica A 281 (2000) 112–128 113 problem of crossover between the Ising universality class and mean-ÿeld critical be- havior. This crossover occurs, for instance, when the interaction range (and hence the “Ginzburg number” G entering the Ginzburg criterion [5]) is varied [6–12]. A closely related crossover is found for symmetrical polymer mixtures when the chain length N of the polymers is varied [4,13–23]. A part of this crossover (though typically not the full extent of the crossover scaling function) can be probed experimentally near the critical point of uids and uid binary mixtures [24–27]. While the Ginzburg criteria [5,14–16] provide a qualitative understanding of this crossover, the quantitatively ac- curate theoretical prediction of the crossover scaling function is a challenging problem [28–36], and hence Monte Carlo studies [6–12,17–19] are of great potential beneÿt. In particular, the question as to what extent (if at all) such crossover scaling functions are universal is an intriguing one [10–12,24–26,36]. Another very interesting crossover which can also be studied is that which occurs near the critical point of unmixing for polymer solutions in a bad solvent [13,37–46]. For chain length N →∞ the critical temperature Tc(N) moves towards the -temperature, where a single coil undergoes a transition from a swollen coil to a collapsed globule. This limit corresponds to a tricritical point [13]. Monte Carlo analyses of critical phenomena typically apply ÿnite-size scaling con- cepts [47–52]. However, care is necessary in the proper application of these methods in the mean-ÿeld limit. In fact, the standard formulation of ÿnite-size scaling (“linear dimensions L scale with the correlation length ”) implies that the hyperscaling rela- tion [2,3] between critical exponents should hold [48,53], which is not the case for mean-ÿeld exponents (apart from d = d∗ = 4 dimensions). This problem already arises for Ising models with short-range interactions for d¿d∗ [54–66], and some disagree- ments between Monte Carlo results [54,55,57] and theoretical predictions [56,65] have stimulated a long-standing debate (see Ref. [66] for a detailed review). 2. Mean-ÿeld to Ising crossover We consider the Hamiltonian [6,7] X X X H=kBT = − K(ri − rj)sisj − h0 si ; (1) i j¿i i −d with si =1 and an interaction K(r) ≡ cR for |r|6R and zero elsewhere. The critical behavior of this model on d-dimensional lattices can be studied eciently with a new cluster algorithm adapted for long-range interactions [67]. To analyze the crossover it is instructive to consider the associated Ginzburg–Landau ÿeld theory in continuous space, Z ( Z 1 h c i 1 H()=k T = − dr dr (r)(r0) − v2(r) B d V 2 |r−r0|6R R 2 ) 4 −u0 (r)+h0(r) ; (2) 114 K. Binder et al. / Physica A 281 (2000) 112–128 Fig. 1. Qualitative picture of the renormalization trajectory describing the crossover from the Gaussian ÿxed ∗ ∗ ∗ ∗ point 0 =(0; 0) to the Ising ÿxed point =(r0 ;u ). where (r) is the single-component order-parameter ÿeld, v is a temperature-like parameter and u0 is a constant. After Fourier transformation and suitable rescaling this can be rewritten as (here N is the total number of lattice sites) 1 X h r i H =k T = k2 + 0 B 2 R2 k −k k r u X X X h N + − ; (3) 4R4N k1 k2 k3 −k1−k2−k3 R 2 k=0 k1 k2 k3 where u is related to u0 and h to h0 [7], and r0 in mean-ÿeld theory is the deviation of the temperature from its critical-point value. We are now interested in identifying the crossover scaling variable associated with the crossover from the Gaussian ÿxed point u = 0 and r0 = 0 to the nontrivial Ising ÿxed point (Fig. 1). Because of the trivial character of the Gaussian ÿxed point and the fact the crossover scaling description should hold all the way from the Ising ÿxed 4=(4−d) point to the Gaussian ÿxed point, one can infer the crossover length scale l0 = R exactly! This is done by considering a renormalization by a length scale l, such that the wave number changes from k to k0 = kl, the number of degrees of freedom is 0 −d 0 −1 reduced from N to N = Nl , and k changes into k0 = l k to leave H invariant. From inspection of the terms in the Hamiltonian one can conclude that the singular part of the free energy must satisfy the scaling relation r u h r u h f˜ 0 ; ; = l−df˜ 0 l2; l4−d; l1+d=2 : (4) s R2 R4 R s R2 R4 R ˜ We see that a ÿnite and nonzero value for the second argument of fs is retained exactly when l takes the value of the crossover scale l0. Thus, we conclude that the singular K. Binder et al. / Physica A 281 (2000) 112–128 115 Fig. 2. E ective susceptibility exponent ∗ above T for the three-dimensional Ising model with variable e c interaction range R (numbers in the key) plotted versus t=G, along with three theoretical calculations for this quantity; due to Refs. [34] (BK), [30] (BB), and [28,32] (SF), respectively. From Ref. [10]. part of the free energy scales with R as follows: ˜ −4d=(4−d) ˆ 2d=(4−d) 3d=(4−d) fs = R fs(˜r0R ; u;˜ hR ) ; (5) where a natural choice of coordinates (Fig. 1) is to measurer ˜0 andu ˜ 0 as distances from the Ising ÿxed point, unlike in the original Hamiltonian, where r0 and u0 are distances from the Gaussian ÿxed point. Eq. (5) describes how the temperature distancer ˜0 from criticality and the magnetic ÿeld h scale with the range of interaction R: Note that the crossover exponent is known exactly (unlike other cases of crossover, e.g., between the Ising and Heisenberg universality class in isotropic magnets with varying uniaxial anisotropy [68]). The same result for the crossover exponent follows [6], of course, from simple-minded arguments using the Ginzburg criterion. However, the location of the nontrivial ÿxed point u∗ (Fig. 1), the associated other exponents, and the explicit form of the scaling function ˜ fs cannot be obtained exactly. The calculation of the scaling function for the free-energy density or its derivatives, such as the susceptibility, is a nontrivial task for both renormalization-group and Monte Carlo calculations. This is demonstrated in Fig. 2 where the e ective critical exponent + e of the susceptibility for T¿Tc is plotted versus the thermal crossover scaling −6 variable t=G, with t =(T − Tc)=Tc being the reduced temperature and G = G0R the Ginzburg number in d = 3, for which G0 ≈ 0:277. Note that e ective exponents are deÿned as e ≡−dlnˆ=dln|t|; ˆ ≡ kBTc(R)(@M=@h)T ; (6) 116 K. Binder et al. / Physica A 281 (2000) 112–128 where refers to T ? Tc, respectively, M = hsiT;h, and the range R is deÿned from the second moment of the interaction (z being the e ective coordination number) , X X 2 2 R = |ri − rj| K(ri − rj) K(ri − rj) j6=i j6=i 1 X = |r − r |2 with |r − r |6R : (7) z i j i j m j6=i 2 Here the second equality holds only for a square-well potential and values Rm = 1; 2; 3; 4; 5; 6; 8; 12; 18; 28; 60; 100, and 160 were studied.
Recommended publications
  • Universal Scaling Behavior of Non-Equilibrium Phase Transitions
    Universal scaling behavior of non-equilibrium phase transitions HABILITATIONSSCHRIFT der FakultÄat furÄ Naturwissenschaften der UniversitÄat Duisburg-Essen vorgelegt von Sven LubÄ eck aus Duisburg Duisburg, im Juni 2004 Zusammenfassung Kritische PhÄanomene im Nichtgleichgewicht sind seit Jahrzehnten Gegenstand inten- siver Forschungen. In Analogie zur Gleichgewichtsthermodynamik erlaubt das Konzept der UniversalitÄat, die verschiedenen NichtgleichgewichtsphasenubÄ ergÄange in unterschied- liche Klassen einzuordnen. Alle Systeme einer UniversalitÄatsklasse sind durch die glei- chen kritischen Exponenten gekennzeichnet, und die entsprechenden Skalenfunktionen werden in der NÄahe des kritischen Punktes identisch. WÄahrend aber die Exponenten zwischen verschiedenen UniversalitÄatsklassen sich hÄau¯g nur geringfugigÄ unterscheiden, weisen die Skalenfunktionen signi¯kante Unterschiede auf. Daher ermÄoglichen die uni- versellen Skalenfunktionen einerseits einen emp¯ndlichen und genauen Nachweis der UniversalitÄatsklasse eines Systems, demonstrieren aber andererseits in ubÄ erzeugender- weise die UniversalitÄat selbst. Bedauerlicherweise wird in der Literatur die Betrachtung der universellen Skalenfunktionen gegenubÄ er der Bestimmung der kritischen Exponen- ten hÄau¯g vernachlÄassigt. Im Mittelpunkt dieser Arbeit steht eine bestimmte Klasse von Nichtgleichgewichts- phasenubÄ ergÄangen, die sogenannten absorbierenden PhasenubÄ ergÄange. Absorbierende PhasenubÄ ergÄange beschreiben das kritische Verhalten von physikalischen, chemischen sowie biologischen
    [Show full text]
  • Three-Dimensional Phase Transitions in Multiflavor Lattice Scalar SO (Nc) Gauge Theories
    Three-dimensional phase transitions in multiflavor lattice scalar SO(Nc) gauge theories Claudio Bonati,1 Andrea Pelissetto,2 and Ettore Vicari1 1Dipartimento di Fisica dell’Universit`adi Pisa and INFN Largo Pontecorvo 3, I-56127 Pisa, Italy 2Dipartimento di Fisica dell’Universit`adi Roma Sapienza and INFN Sezione di Roma I, I-00185 Roma, Italy (Dated: January 1, 2021) We investigate the phase diagram and finite-temperature transitions of three-dimensional scalar SO(Nc) gauge theories with Nf ≥ 2 scalar flavors. These models are constructed starting from a maximally O(N)-symmetric multicomponent scalar model (N = NcNf ), whose symmetry is par- tially gauged to obtain an SO(Nc) gauge theory, with O(Nf ) or U(Nf ) global symmetry for Nc ≥ 3 or Nc = 2, respectively. These systems undergo finite-temperature transitions, where the global sym- metry 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 ≥ 3. 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 ≥ 3. We perform numerical simulations for Nc = 3 and Nf = 2, 3. The numerical results are in agreement with the LGW predictions. I. INTRODUCTION The global O(N) symmetry is partially gauged, obtain- ing a nonabelian gauge model, in which the fields belong to the coset SN /SO(N ), where SN = SO(N)/SO(N 1) Global and local gauge symmetries play a crucial role c − in theories describing fundamental interactions [1] and is the N-dimensional sphere.
    [Show full text]
  • Universal Scaling Behavior of Non-Equilibrium Phase Transitions
    Universal scaling behavior of non-equilibrium phase transitions Sven L¨ubeck Theoretische Physik, Univerit¨at Duisburg-Essen, 47048 Duisburg, Germany, [email protected] December 2004 arXiv:cond-mat/0501259v1 [cond-mat.stat-mech] 11 Jan 2005 Summary Non-equilibrium critical phenomena have attracted a lot of research interest in the recent decades. Similar to equilibrium critical phenomena, the concept of universality remains the major tool to order the great variety of non-equilibrium phase transitions systematically. All systems belonging to a given universality class share the same set of critical exponents, and certain scaling functions become identical near the critical point. It is known that the scaling functions vary more widely between different uni- versality classes than the exponents. Thus, universal scaling functions offer a sensitive and accurate test for a system’s universality class. On the other hand, universal scaling functions demonstrate the robustness of a given universality class impressively. Unfor- tunately, most studies focus on the determination of the critical exponents, neglecting the universal scaling functions. In this work a particular class of non-equilibrium critical phenomena is considered, the so-called absorbing phase transitions. Absorbing phase transitions are expected to occur in physical, chemical as well as biological systems, and a detailed introduc- tion is presented. The universal scaling behavior of two different universality classes is analyzed in detail, namely the directed percolation and the Manna universality class. Especially, directed percolation is the most common universality class of absorbing phase transitions. The presented picture gallery of universal scaling functions includes steady state, dynamical as well as finite size scaling functions.
    [Show full text]
  • Sankar Das Sarma 3/11/19 1 Curriculum Vitae
    Sankar Das Sarma 3/11/19 Curriculum Vitae Sankar Das Sarma Richard E. Prange Chair in Physics and Distinguished University Professor Director, Condensed Matter Theory Center Fellow, Joint Quantum Institute University of Maryland Department of Physics College Park, Maryland 20742-4111 Email: [email protected] Web page: www.physics.umd.edu/cmtc Fax: (301) 314-9465 Telephone: (301) 405-6145 Published articles in APS journals I. Physical Review Letters 1. Theory for the Polarizability Function of an Electron Layer in the Presence of Collisional Broadening Effects and its Experimental Implications (S. Das Sarma) Phys. Rev. Lett. 50, 211 (1983). 2. Theory of Two Dimensional Magneto-Polarons (S. Das Sarma), Phys. Rev. Lett. 52, 859 (1984); erratum: Phys. Rev. Lett. 52, 1570 (1984). 3. Proposed Experimental Realization of Anderson Localization in Random and Incommensurate Artificial Structures (S. Das Sarma, A. Kobayashi, and R.E. Prange) Phys. Rev. Lett. 56, 1280 (1986). 4. Frequency-Shifted Polaron Coupling in GaInAs Heterojunctions (S. Das Sarma), Phys. Rev. Lett. 57, 651 (1986). 5. Many-Body Effects in a Non-Equilibrium Electron-Lattice System: Coupling of Quasiparticle Excitations and LO-Phonons (J.K. Jain, R. Jalabert, and S. Das Sarma), Phys. Rev. Lett. 60, 353 (1988). 6. Extended Electronic States in One Dimensional Fibonacci Superlattice (X.C. Xie and S. Das Sarma), Phys. Rev. Lett. 60, 1585 (1988). 1 Sankar Das Sarma 7. Strong-Field Density of States in Weakly Disordered Two Dimensional Electron Systems (S. Das Sarma and X.C. Xie), Phys. Rev. Lett. 61, 738 (1988). 8. Mobility Edge is a Model One Dimensional Potential (S.
    [Show full text]
  • Arxiv:1908.10990V1 [Cond-Mat.Stat-Mech] 29 Aug 2019
    High-precision Monte Carlo study of several models in the three-dimensional U(1) universality class Wanwan Xu,1 Yanan Sun,1 Jian-Ping Lv,1, ∗ and Youjin Deng2, 3 1Anhui Key Laboratory of Optoelectric Materials Science and Technology, Key Laboratory of Functional Molecular Solids, Ministry of Education, Anhui Normal University, Wuhu, Anhui 241000, China 2Hefei National Laboratory for Physical Sciences at Microscale and Department of Modern Physics, University of Science and Technology of China, Hefei, Anhui 230026, China 3CAS Center for Excellence and Synergetic Innovation Center in Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China We present a worm-type Monte Carlo study of several typical models in the three-dimensional (3D) U(1) universality class, which include the classical 3D XY model in the directed flow representation and its Vil- lain version, as well as the 2D quantum Bose-Hubbard (BH) model with unitary filling in the imaginary-time world-line representation. From the topology of the configurations on a torus, we sample the superfluid stiff- ness ρs and the dimensionless wrapping probability R. From the finite-size scaling analyses of ρs and of R, we determine the critical points as Tc(XY) = 2.201 844 1(5) and Tc(Villain) = 0.333 067 04(7) and (t/U)c(BH) = 0.059 729 1(8), where T is the temperature for the classical models, and t and U are respec- tively the hopping and on-site interaction strength for the BH model. The precision of our estimates improves significantly over that of the existing results.
    [Show full text]
  • Crossover Between Mean-Field and Ising Critical Behavior in a Lattice-Gas Reaction-Diffusion Model Da-Jiang Liu the Ames Laboratory
    Physics and Astronomy Publications Physics and Astronomy 2004 Crossover Between Mean-Field and Ising Critical Behavior in a Lattice-Gas Reaction-Diffusion Model Da-Jiang Liu The Ames Laboratory N. Pavlenko Universitat Hannover James W. Evans Iowa State University, [email protected] Follow this and additional works at: http://lib.dr.iastate.edu/physastro_pubs Part of the Chemistry Commons, and the Physics Commons The ompc lete bibliographic information for this item can be found at http://lib.dr.iastate.edu/physastro_pubs/456. For information on how to cite this item, please visit http://lib.dr.iastate.edu/howtocite.html. This Article is brought to you for free and open access by the Physics and Astronomy at Iowa State University Digital Repository. It has been accepted for inclusion in Physics and Astronomy Publications by an authorized administrator of Iowa State University Digital Repository. For more information, please contact [email protected]. Crossover Between Mean-Field and Ising Critical Behavior in a Lattice- Gas Reaction-Diffusion Model Abstract Lattice-gas models for CO oxidation can exhibit a discontinuous nonequilibrium transition between reactive and inactive states, which disappears above a critical CO-desorption rate. Using finite-size-scaling analysis, we demonstrate a crossover from Ising to mean-field behavior at the critical point, with increasing surface mobility of adsorbed CO or with decreasing system size. This behavior is elucidated by analogy with that of equilibrium Ising-type systems with long-range interactions. Keywords critical behavior, lattice-gas reaction-diffusion model, Ising, mean field Disciplines Chemistry | Physics Comments This article is published as Liu, Da-Jiang, N.
    [Show full text]
  • Directed Percolation and Turbulence
    Emergence of collective modes, ecological collapse and directed percolation at the laminar-turbulence transition in pipe flow Hong-Yan Shih, Tsung-Lin Hsieh, Nigel Goldenfeld University of Illinois at Urbana-Champaign Partially supported by NSF-DMR-1044901 H.-Y. Shih, T.-L. Hsieh and N. Goldenfeld, Nature Physics 12, 245 (2016) N. Goldenfeld and H.-Y. Shih, J. Stat. Phys. 167, 575-594 (2017) Deterministic classical mechanics of many particles in a box statistical mechanics Deterministic classical mechanics of infinite number of particles in a box = Navier-Stokes equations for a fluid statistical mechanics Deterministic classical mechanics of infinite number of particles in a box = Navier-Stokes equations for a fluid statistical mechanics Transitional turbulence: puffs • Reynolds’ original pipe turbulence (1883) reports on the transition Univ. of Manchester Univ. of Manchester “Flashes” of turbulence: Precision measurement of turbulent transition Q: will a puff survive to the end of the pipe? Many repetitions survival probability = P(Re, t) Hof et al., PRL 101, 214501 (2008) Pipe flow turbulence Decaying single puff metastable spatiotemporal expanding laminar puffs intermittency slugs Re 1775 2050 2500 푡−푡 − 0 Survival probability 푃 Re, 푡 = 푒 휏(Re) ) Re,t Puff P( lifetime to N-S Avila et al., (2009) Avila et al., Science 333, 192 (2011) Hof et al., PRL 101, 214501 (2008) 6 Pipe flow turbulence Decaying single puff Splitting puffs metastable spatiotemporal expanding laminar puffs intermittency slugs Re 1775 2050 2500 푡−푡 − 0 Splitting
    [Show full text]
  • Block Scaling in the Directed Percolation Universality Class
    OPEN ACCESS Recent citations Block scaling in the directed percolation - 25 Years of Self-organized Criticality: Numerical Detection Methods universality class R. T. James McAteer et al - The Abelian Manna model on various To cite this article: Gunnar Pruessner 2008 New J. Phys. 10 113003 lattices in one and two dimensions Hoai Nguyen Huynh et al View the article online for updates and enhancements. This content was downloaded from IP address 170.106.40.139 on 26/09/2021 at 04:54 New Journal of Physics The open–access journal for physics Block scaling in the directed percolation universality class Gunnar Pruessner1 Mathematics Institute, University of Warwick, Gibbet Hill Road, Coventry CV4 7AL, UK E-mail: [email protected] New Journal of Physics 10 (2008) 113003 (13pp) Received 23 July 2008 Published 7 November 2008 Online at http://www.njp.org/ doi:10.1088/1367-2630/10/11/113003 Abstract. The universal behaviour of the directed percolation universality class is well understood—both the critical scaling and the finite size scaling. This paper focuses on the block (finite size) scaling of the order parameter and its fluctuations, considering (sub-)blocks of linear size l in systems of linear size L. The scaling depends on the choice of the ensemble, as only the conditional ensemble produces the block-scaling behaviour as established in equilibrium critical phenomena. The dependence on the ensemble can be understood by an additional symmetry present in the unconditional ensemble. The unconventional scaling found in the unconditional ensemble is a reminder of the possibility that scaling functions themselves have a power-law dependence on their arguments.
    [Show full text]
  • QCD Critical Point: Fluctuations, Hydrodynamics and Universality Class
    QCD critical point: fluctuations, hydrodynamics and universality class M. Stephanov U. of Illinois at Chicago and RIKEN-BNL QCD critical point: fluctuations, hydrodynamics and universality class – p.1/21 QCD critical point hot QGP T, GeV 3 1 E critical 0.1 point 2 nuclear CFL vacuum matter quark matter 0 1 µB , GeV 1 Lattice at µ =0 → crossover. 2 Sign problem. Models: 1st order transition. 1 + 2 = 3 : critical point E. (As in water at p = 221bar, T = 373◦C – critical opalescence.) Where is point E? Challenge to theory and experiment. QCD critical point: fluctuations, hydrodynamics and universality class – p.2/21 Location of the CP (theory) Source (T, µB ), MeV Comments Label MIT Bag/QGP none only 1st order — Asakawa,Yazaki ’89 (40, 1050) NJL, CASE I NJL/I “ (55, 1440) NJL, CASE II NJL/II Barducci, et al ’89-94 (75, 273)TCP composite operator CO Berges, Rajagopal ’98 (101, 633)TCP instanton NJL NJL/inst Halasz, et al ’98 (120, 700)TCP random matrix RM Scavenius, et al ’01 (93,645) linear σ-model LSM “ (46,996) NJL NJL Fodor, Katz ’01 (160, 725) lattice reweighting I Hatta, Ikeda, ’02 (95, 837) effective potential (CJT) CJT Antoniou, Kapoyannis ’02 (171, 385) hadronic bootstrap HB Ejiri, et al ’03 (?,420) lattice Taylor expansion Fodor, Katz ’04 (162, 360) lattice reweighting II QCD critical point: fluctuations, hydrodynamics and universality class – p.3/21 Location of the CP (theory) 200 T HB lat. Taylor exp. lat. reweighting I 150 lat. reweighting II NJL/inst RM 100 CO LSM CJT NJL/II NJL 50 NJL/I 0 0 200 400 600 800 1000 1200 1400 1600 µB Slope: Allton et al ’02 (Taylor exp.), de Forcrand, Philipsen ’02 (Imµ): 2 T (µ ) µ B ≈ 1 − .00(5 − 8) B T (0) T (0) Allton et al ’03: .013.
    [Show full text]
  • [Cond-Mat.Stat-Mech] 26 Jul 1999
    Ordered phase and scaling in Zn models and the three-state antiferromagnetic Potts model in three dimensions Masaki Oshikawa Department of Physics, Tokyo Institute of Technology, Oh-okayama, Meguro-ku, Tokyo 152-8551, Japan (July 20, 1999) the spontaneously broken Zn symmetry and the high- Based on a Renormalization-Group picture of Zn symmet- temperature disordered phase. The intermediate phase is ric models in three dimensions, we derive a scaling law for the O(2) symmetric and corresponds to the low-temperature Zn order parameter in the ordered phase. An existing Monte phase of the XY model. Carlo calculation on the three-state antiferromagnetic Potts On three dimensional (3D) case, Blankschtein et al.3 model, which has the effective Z6 symmetry, is shown to be in 1984 proposed an RG picture of the Z models, to consistent with the proposed scaling law. It strongly supports 6 the Renormalization-Group picture that there is a single mas- discuss the STI model. They suggested that the transi- sive ordered phase, although an apparently rotationally sym- tion between the ordered and disordered phases belongs metric region in the intermediate temperature was observed to the (3D) XY universality class, and that the ordered numerically. phase reflects the symmetry breaking to Z6 in a large enough system. It means that there is no finite region of 75.10.Hk, 05.50+q rotationally symmetric phase which is similar to the or- dered phase of the XY model. Unfortunately, their paper is apparently not widely known in the related fields. It might be partly because their discussion was very brief I.
    [Show full text]
  • Disorder Driven Itinerant Quantum Criticality of Three Dimensional
    Disorder driven itinerant quantum criticality of three dimensional massless Dirac fermions J. H. Pixley, Pallab Goswami, and S. Das Sarma Condensed Matter Theory Center and Joint Quantum Institute, Department of Physics, University of Maryland, College Park, Maryland 20742- 4111 USA (Dated: August 31, 2018) Progress in the understanding of quantum critical properties of itinerant electrons has been hin- dered by the lack of effective models which are amenable to controlled analytical and numerically exact calculations. Here we establish that the disorder driven semimetal to metal quantum phase transition of three dimensional massless Dirac fermions could serve as a paradigmatic toy model for studying itinerant quantum criticality, which is solved in this work by exact numerical and approx- imate field theoretic calculations. As a result, we establish the robust existence of a non-Gaussian universality class, and also construct the relevant low energy effective field theory that could guide the understanding of quantum critical scaling for many strange metals. Using the kernel polynomial method (KPM), we provide numerical results for the calculated dynamical exponent (z) and corre- lation length exponent (ν) for the disorder-driven semimetal (SM) to diffusive metal (DM) quantum phase transition at the Dirac point for several types of disorder, establishing its universal nature and obtaining the numerical scaling functions in agreement with our field theoretical analysis. PACS numbers: 71.10.Hf,72.80.Ey,73.43.Nq,72.15.Rn I. INTRODUCTION lengths (ξl and ξt) are infinite. Even though the QCP occurs strictly at zero temperature, its effects on physical Phase transitions are ubiquitous in the natural world properties manifest over a broad range of temperatures ranging from the solidification of water to the thermal (T ) and energies (E), such that suppression of magnetism in iron or superconductivity in 1 E, k T ξ− , metals.
    [Show full text]
  • Non-Equilibrium Phase Transitions an Introduction
    Directed Percolation Other Universality Classes Non-equilibrium phase transitions An Introduction Lecture III Haye Hinrichsen University of Würzburg, Germany March 2006 Directed Percolation Other Universality Classes Third Lecture: Outline 1 Directed Percolation Scaling Theory Langevin Equation 2 Other Universality Classes Parity-conserving Class Voter Universality Class DP with long-range interactions Directed Percolation Other Universality Classes Outline 1 Directed Percolation Scaling Theory Langevin Equation 2 Other Universality Classes Parity-conserving Class Voter Universality Class DP with long-range interactions Directed Percolation Other Universality Classes Phenomenological Scaling Theory – Scale Invariance Scaling hypothesis: In the scaling regime the large-scale properties of absorbing phase transitions are invariant under scale transformations (zoom in – zoom out) ∆ → λ ∆ (∆ = p − pc) ρ → λβρ 0 P → λβ P −ν⊥ ξ⊥ → λ ξ⊥ −νk ξk → λ ξk 0 with the four exponents β, β , ν⊥, νk. This is called simple scaling (opposed to multiscaling). Directed Percolation Other Universality Classes Using Scale Invariance If you have an expression... multiply the distance from criticality p − pc by λ. multiply all quantities that measure local activity by λβ 0 multiply all quantities that activate sites locally by λβ multiply all lengths by λ−ν⊥ multiply all times by λ−νk ...then this expression should be asymptotically invariant. Directed Percolation Other Universality Classes Using Scale Invariance 1st Example: Decay of the density at criticality: −α β −νk −α ρ ∼ t ⇒ λ ρ ∼ (λ t) ⇒ α = β/νk The same more formally: ρ(t) = λ−βρ(λ−νk t) Set λ = t1/νk . ⇒ ρ(t) = t−β/ν⊥ ρ(1) Directed Percolation Other Universality Classes Using Scale Invariance 2nd Example: Decay of the density for ∆ 6= 0: ρ(∆, t) = λ−βρ(λ∆, λ−νk t) Set again λ = t1/νk .
    [Show full text]