Hard-Disk Equation of State: First-Order Liquid-Hexatic Transition in Two Dimensions with Three Simulation Methods

ORAL

Abstract

We report large-scale computer simulations of the hard-disk system at high densities in the region of the melting transition~[1]. Our simulations reproduce the equation of state, previously obtained using the event-chain Monte Carlo algorithm, with a massively parallel implementation of the local Monte Carlo method~[2] and with event-driven molecular dynamics. We analyze the relative performance of these simulation methods to sample configuration space and approach equilibrium. Phase coexistence is visualized for individual configurations via the local orientations, and positional correlation functions are computed. Our results confirm the first-order nature of the liquid-hexatic phase transition in hard disks. \\[4pt] [1] J.A. Anderson, M. Engel, S.C. Glotzer, M. Isobe, E.P. Bernard, W. Krauth, arXiv:1211.1645. \newline [2] J.A. Anderson, E. Jankowski, T.L. Grubb, M. Engel, S.C. Glotzer, arXiv:1211.1646.

Authors

  • Michael Engel

    Department of Chemical Engineering, University of Michigan, University of Michigan, University of Michigan, Dept. of Chem. Eng., Chemical Engineering Department

  • Joshua Anderson

    University of Michigan, Department of Chemical Engineering, University of Michigan, Ann Arbor

  • Masaharu Isobe

    Nagoya Institute of Technology

  • Etienne P. Bernard

    Massachusetts Institute of Technology

  • Werner Krauth

    \'{E}cole Normale Sup\'{e}rieure

  • Sharon Glotzer

    Department of Chemical Engineering, University of Michigan, University of Michigan, Department of Chemical Engineering, University of Michigan, Ann Arbor, Materials Science and Engineering Department, University of Michigan, Ann Arbor MI