DOCSLIB.ORG

Critical Behavior of the Extended Hubbard Model with Bond Dimerization

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
Critical Behavior of the Extended Hubbard Model with Bond Dimerization Critical behavior of the extended Hubbard model with bond dimerization Satoshi Ejimaa, Florian Langea,b, Fabian H. L. Esslerc, Holger Fehskea a Institut f¨urPhysik, Ernst-Moritz-Arndt-Universit¨atGreifswald, 17487 Greifswald, Germany b Computational Condensed Matter Physics Laboratory, RIKEN, Wako, Saitama 351-0198, Japan c The Rudolf Peierls Centre for Theoretical Physics, Oxford University, Oxford OX1 3NP, United Kingdom Abstract Exploiting the matrix-product-state based density-matrix renormalization group (DMRG) technique we study the one-dimensional extended (U-V) Hubbard model with explicit bond dimerization in the half-filled band sector. In particular we investigate the nature of the quantum phase transition, taking place with growing ratio V=U between the symmetry-protected-topological and charge-density-wave insulating states. The (weak-coupling) critical line of continuous Ising transitions with central charge c = 1=2 terminates at a tricritical point belonging to the universality class of the dilute Ising model with c = 7=10. We demonstrate that our DMRG data perfectly match with (tricritical) Ising exponents, e.g., for the order parameter β = 1=8 (1/24) and correlation length ν = 1 (5/9). Beyond the tricritical Ising point, in the strong-coupling regime, the quantum phase transition becomes first order. Keywords: Extended Hubbard model, (tricritical) Ising universality class 1. Introduction methods. Half a century has passed since it was proposed, yet At present, quantum phase transitions between topo- the Hubbard model [1] is still a key Hamiltonian for logically trivial and nontrivial states arouse great in- the investigation of strongly correlated electron sys- terest [10, 11, 12]. In this context, extensions of tems. Originally designed to describe the ferromag- the half-filled EHM also attracted attention, mainly netism of transition metals, in successive studies the with regard to the formation of symmetry-protected- Hubbard model has also been used for heavy fermions topological (SPT) states [11]. Including an alternating and high-temperature superconductors. The physics of ferromagnetic spin interaction [13] or an explicit dimer- the model is governed by the competition between the ization [14] in the EHM, the SDW and BOW phases itinerancy of the charge carriers and their local Coulomb are completely replaced by an SPT insulator, whereby a interaction. In one dimension (1D), seen from a theoret- quantum phase transition occurs between the SPT and ical point of view, the Hubbard model is a good starting the CDW, the area of which shrinks. Most interestingly, point to explore, for example, Tomonaga-Luttinger liq- the SPT-CDW continuous Ising transition with central uid behavior (including spin-charge separation). charge c = 1=2 ends at a tricritical point, belonging to While the 1D Hubbard model is exactly solvable by the universality class of the tricritical Ising model, a sec- Bethe Ansatz [2], most of its extensions are no longer ond minimal model with c = 7=10 [15, 16]. Above this integrable. This is even true if only the Coulomb in- point, the quantum phase transition becomes first order. teraction between electrons on nearest-neighbor lattice In Ref. [14] it has been demonstrated that the transi- sites is added. The ground-state phase diagram of this tion region of the EHM with bond dimerization can be so-called extended Hubbard model (EHM) is still a described by the triple sine-Gordon model by extend- arXiv:1706.07486v1 [cond-mat.str-el] 22 Jun 2017 hotly debated issue. At half filling, this relates in partic- ing the former bosonization analysis [17]. The predic- ular to the recently discovered bond-order-wave (BOW) tions of field theory regarding power-law (exponential) state located in between spin-density-wave (SDW) and decay of the density-density (spin-spin) and bond-order charge-density-wave (CDW) phases [3, 4]. To charac- correlation functions are shown to be in excellent ac- terize the BOW state and determine its phase boundaries cordance with the numerical data obtained by a matrix- considerable efforts were undertaken in the last few product-states (MPS) based density-matrix renormal- years, using both analytical [5, 6] and numerical [7, 8, 9] ization group (DMRG) technique [18, 19]. Preprint submitted to Physica B June 26, 2017 6 first order χ = 400 5 100 CDW (∆c > 0, ∆s > 0) 5 U=t = 4 χ ξ 4 tricritical point critical end point 50 4 (a) V=t 3 χ = 200 0 2 continuous (c = 1) BOW (∆c > 0, ∆s > 0) 1:9 2:0 2:1 2:2 V=t SDW (∆c > 0, ∆s = 0) 1 2:1 (b) V = 2:160t (≈ Vc) 0 0 2 4 6 8 10 12 U=t χ 2 ∗ S c ' 0:996 Figure 1: DMRG ground-state phase diagram of the 1D EHM (1) at half filling [9]. The red dotted line gives the continuous SDW-BOW y = a · x + b transition. The bold (thin) blue dashed line marks the continuous 1:9 (first-order) BOW-CDW transition and the green dashed-dotted line 5:6 5:8 6 6:2 6:4 denotes the first-order SDW-CDW transition. ln ξχ The Ising criticality of the EHM with explicit dimer- Figure 2: (a): Correlation length ξχ of the EHM as a function of V=t ization was established in early work [17] that also spec- for U=t = 4 obtained from iDMRG. The dashed line indicates the ifies the critical exponents. The critical exponents at BOW-CDW transition point. (b): von Neumann entropy S χ as a func- tion of logarithm of ξχ at V ≈ Vc for U=t = 4. The iDMRG data the tricritical point should differ from those at the ordi- for ln ξχ > 6 (χ ≥ 1800) provide us the numerically obtained central nary Ising transition because the tricritical Ising quan- charge c∗ ' 0:996 by fitting to Eq. (2). tum phase transition belongs to a different universality class. Simulating the neutral gap and the CDW order pa- system shows fluctuating SDW order. The spin (charge) rameter by DMRG, in this paper we will determine the excitations are gapless (gapped) 8U > 0 [2]. At fi- critical exponents at both Ising and tricritical Ising tran- nite V, for V=U . 1=2, the ground state is still a sitions. The paper is structured as follows. Section 2 SDW. When V=U becomes larger than 1/2 a 2kF-CDW introduces the model Hamiltonians under consideration is formed. As pointed out first by Nakamura [3, 4] and discusses their ground-state properties. The critical and confirmed later by various analytical and numeri- exponents will be derived in Sect. 3. Section 4 summa- cal studies [8, 9, 20, 21], the SDW and CDW phases are rizes our main results. separated by a narrow BOW phase below the critical EHM EHM end point, (Uce ,Vce ) ≈ (9:25t; 4:76t). In the BOW phase translational symmetry is spontaneously broken, 2. Model which implies that the spin gap opens passing the SDW- EHM 2.1. Extended Hubbard model BOW phase boundary at fixed U < Uce . Increasing V further, the system enters the CDW phase with finite The Hamiltonian of the EHM is defined as spin and charge gaps. The BOW-CDW Gaussian tran- X ˆ y sition line with central charge c = 1 terminates at the HEHM = −t (ˆc jσcˆ + H:c:) j+1σ EHM EHM jσ tricritical point, (Utr ; Vtr ) ≈ (5:89t; 3:10t) [9]. For ! ! EHM EHM X 1 1 Utr < U < Uce , the BOW-CDW transition becomes +U nˆ j" − nˆ j# − first order, characterized by a jump in the spin gap (see, 2 2 j Fig. 3 in Ref. [9]). Figure 1 summarizes the rich physics X +V (ˆn j − 1)(ˆn j+1 − 1) ; (1) of the half-filled EHM. j The criticality at the continuous BOW-CDW transi- tion line can be verified numerically by extracting, e.g., wherec ˆy (ˆc ) creates (annihilates) an electron with jσ jσ the central charge from the the correlation length (ξχ) y spin projection σ ="; # at Wannier site j,n ˆ jσ = cˆ jσcˆ jσ, and von Neumann entropy (S χ), where ξχ can be ob- andn ˆ j = nˆ j" + nˆ j#. In the Hubbard model limit (V = 0), tained from the second largest eigenvalue of the trans- at half-filling, no long-range order exists. Instead the fer matrix for some bond dimension χ used in a infinite 2 8 Previous studies of this model have shown that the low δ=t = 0:2 first order lying excitations in the large-U limit are chargeless spin-triplet and spin-singlet excitations [24, 25, 26, 27, 6 28, 29, 30], whereby the dynamics is described by an ef- tricritical Ising fective spin-Peierls Hamiltonian. Moreover, at finite U, point (c = 7=10) the Tomonaga-Luttinger parameters have been explored 4 CDW V=t at and near commensurate fillings by DMRG [31]. Par- continuous Ising ticularly for half filling, it has been proven by pertur- (c = 1=2) PI 2 bative [32, 33] and renormalization group [6, 34, 35] approaches that the system realizes Peierls insulator +δ −δ (PI) and CDW phases in the weak-coupling regime. 0 According to weak-coupling renormalization-group re- 0 2 4 6 8 10 12 U=t sults [6], any finite bond dimerization δ will change the universality class of the continuous BOW-CDW transi- Figure 3: Ground-state phase diagram of the 1D EHM with bond tion (realized in the pure EHM) from Gaussian to Ising dimerization in the half-filled band sector [14]. The red solid line type. Thereby the PI-CDW transition in the weak-to- marks the PI-CDW phase boundaries for δ/t = 0:2.
Load more

Recommended publications
  • 978-1-62808-450-4 Ch1.Pdf
    In: Recent Advances in Magnetism Research ISBN: 978-1-62808-450-4 Editor: Keith Pace © 2013 Nova Science Publishers, Inc. No part of this digital document may be reproduced, stored in a retrieval system or transmitted commercially in any form or by any means. The publisher has taken reasonable care in the preparation of this digital document, but makes no expressed or implied warranty of any kind and assumes no responsibility for any errors or omissions. No liability is assumed for incidental or consequential damages in connection with or arising out of information contained herein. This digital document is sold with the clear understanding that the publisher is not engaged in rendering legal, medical or any other professional services. Chapter 1 TOWARDS A FIELD THEORY OF MAGNETISM Ulrich Köbler FZ-Jülich, Institut für Festkörperforschung, Jülich, Germany ABSTRACT Experimental tests of conventional spin wave theory reveal two serious shortcomings of this typical local or atomistic theory: first, spin wave theory is unable to explain the universal temperature dependence of the thermodynamic observables and, second, spin wave theory does not distinguish between the dynamics of magnets with integer and half- integer spin quantum number. As experiments show, universality is not restricted to the critical dynamics but holds for all temperatures in the long range ordered state including the low temperature regime where spin wave theory is widely believed to give an adequate description. However, a rigorous and convincing test that spin wave theory gives correct description of the thermodynamics of ordered magnets has never been delivered. Quite generally, observation of universality unequivocally calls for a continuum or field theory instead of a local theory.
    [Show full text]
  • 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]
  • Arxiv:1012.0653V3
    Quantum phase transitions in transverse field spin models: From Statistical Physics to Quantum Information Amit Dutta∗ Department of Physics, Indian Institute of Technology, Kanpur 208 016, India Uma Divakaran† Theoretische Physik, Universit¨at des Saarlandes, 66041 Saarbr¨ucken, Germany Diptiman Sen‡ Centre for High Energy Physics, Indian Institute of Science, Bangalore 560 012, India Bikas K. Chakrabarti§ Theoretical Condensed Matter Physics Division and Centre for Applied Mathematics and Computational Science, Saha Institute of Nuclear Physics, 1/AF Bidhan Nagar, Kolkata 700 064, India Thomas F. Rosenbaum¶ Department of Physics and the James Franck Institute, University of Chicago, Chicago, IL 60637 Gabriel Aeppli∗∗ London Centre for Nanotechnology and Department of Physics and Astronomy, UCL, London, WC1E 6BT, United Kingdom We review quantum phase transitions of spin systems in transverse magnetic fields taking the examples of some paradigmatic models, namely, the spin-1/2 Ising and XY models in a transverse field in one and higher spatial dimensions. Beginning with a brief overview of quantum phase transitions, we introduce the model Hamiltonians and discuss the equivalence between the quantum phase transition in such a model and the finite temperature phase transition in a higher dimensional classical model. We then provide exact solutions in one spatial dimension connecting them to conformal field theoretical studies when possible. We also discuss Kitaev models, and some other exactly solvable spin systems in this context. Studies of quantum phase transitions in the presence of quenched randomness and with frustrating interactions are presented in details. We discuss novel phenomena like Griffiths-McCoy singularities associated with quantum phase transitions of low- dimensional transverse Ising models with random interactions and transverse fields.
    [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]
  • Bootstrapping Hypercubic and Hypertetrahedral Theories in Three Dimensions
    CERN-TH-2018-012 Bootstrapping hypercubic and hypertetrahedral theories in three dimensions Andreas Stergiou Theoretical Physics Department, CERN, Geneva, Switzerland There are three generalizations of the Platonic solids that exist in all dimensions, namely the hypertetrahedron, the hypercube, and the hyperoctahedron, with the latter two being dual. Conformal field theories with the associated symmetry groups as global symmetries can be argued to exist in d = 3 spacetime dimensions if the " = 4 − d expansion is valid when " ! 1. In this paper hypercubic and hypertetrahedral theories are studied with the non-perturbative numerical conformal bootstrap. In the N = 3 cubic case it is found that a bound with a kink is saturated by a solution with properties that cannot be reconciled with the " expansion of the cubic theory. Possible implications for cubic magnets and structural phase transitions are discussed. For the hypertetrahedral theory evidence is found that the non-conformal window that is seen with the " expansion exists in d = 3 as well, and a rough estimate of its extent is given. arXiv:1801.07127v4 [hep-th] 28 Apr 2018 January 2018 1. Introduction The numerical conformal bootstrap [1] in three dimensions has produced impressive results, especially concerning critical theories in universality classes that also contain scalar conformal field theories (CFTs). For the Ising universality class the bootstrap is in fact the state-of-the-art method for the determination of critical exponents [2]. For theories with continuous global symmetries there has also been significant progress, with important developments for the Heisenberg universality class [3{5]. In this work we study critical theories with discrete global symmetries in three dimensions, focusing on the hypercubic and hypertetrahedral symmetry groups.
    [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]