Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
  • Free to Read

Hard-disk equation of state: First-order liquid-hexatic transition in two dimensions with three simulation methods

Michael Engel, Joshua A. Anderson, Sharon C. Glotzer, Masaharu Isobe, Etienne P. Bernard, and Werner Krauth
Phys. Rev. E 87, 042134 – Published 30 April 2013
An article within the collection: Physical Review E 25th Anniversary Milestones
PDFHTMLExport Citation
Abstract
Collection
Authors
Article Text
  • INTRODUCTION
  • SIMULATION METHODS
  • PERFORMANCE COMPARISON
  • RESULTS AND DISCUSSION
  • CONCLUSION
  • ACKNOWLEDGEMENTS
  • APPENDICES
    • Supplemental Material
    References

    Abstract

    We report large-scale computer simulations of the hard-disk system at high densities in the region of the melting transition. 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 and with event-driven molecular dynamics. We analyze the relative performance of these simulation methods to sample configuration space and approach equilibrium. Our results confirm the first-order nature of the melting phase transition in hard disks. Phase coexistence is visualized for individual configurations via the orientational order parameter field. The analysis of positional order confirms the existence of the hexatic phase.

    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    1 More
    • Received 9 November 2012

    DOI:https://doi.org/10.1103/PhysRevE.87.042134

    ©2013 American Physical Society

    Collections

    This article appears in the following collection:

    Physical Review E 25th Anniversary Milestones

    The year 2018 marks the 25th anniversary of Physical Review E. To celebrate the journal’s rich legacy, during the upcoming year we highlight a series of papers that made important contributions to their field. These milestone articles were nominated by members of the Editorial Board of Physical Review E, in collaboration with the journal’s editors. The 25 milestone articles, including an article for each calendar year from 1993 through 2017 and spanning all major subject areas of the journal, will be unveiled in chronological order and will be featured on the journal website.

    Authors & Affiliations

    Michael Engel1, Joshua A. Anderson1, Sharon C. Glotzer1,*, Masaharu Isobe2,3, Etienne P. Bernard4, and Werner Krauth5,†

    • 1Department of Chemical Engineering, University of Michigan, Ann Arbor, Michigan 48109, USA
    • 2Graduate School of Engineering, Nagoya Institute of Technology, Nagoya, 466-8555, Japan
    • 3Department of Chemistry, University of California, Berkeley, California 94720, USA
    • 4Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
    • 5Laboratoire de Physique Statistique, École Normale Supérieure, UPMC, CNRS, 24 Rue Lhomond, 75231 Paris Cedex 05, France

    • *sglotzer@umich.edu
    • †werner.krauth@ens.fr

    Article Text

    Click to Expand

    Supplemental Material

    Click to Expand

    References

    Click to Expand
    Issue

    Vol. 87, Iss. 4 — April 2013

    Reuse & Permissions
    Access Options
    CHORUS

    Article part of CHORUS

    Author publication services for translation and copyediting assistance advertisement

    Authorization Required

    Log In
    Other Options
    ×

    Images

    • Figure 1
      Figure 1
      (a) Pair-correlation function g(r) close to contact for N=5122 at density η=0.698 using LMC. Error bars are computed through 64 independent simulations. g(r) is fitted with a fourth-order polynomial. Difference between g(r) and the polynomial fit with the single-precision data (b) and double-precision data (c) for the histogram.Reuse & Permissions
    • Figure 2
      Figure 2
      Autocorrelation function of the global orientation order parameter Ψ6(t) for N=5122, η=0.698 obtained with LMC, EDMD, ECMC, and MPMC. (a) Time is measured in number of attempted displacements (or collisions) per disk. (b) Time is measured in CPU or GPU hours.Reuse & Permissions
    • Figure 3
      Figure 3
      Equation of state from ECMC, EDMD, and MPMC for N=2562 and N=5122. Error bars are mostly smaller than the symbols. Results agree within one standard deviation. The inset shows the relative pressure difference ΔP/P of EDMD and MPMC with respect to ECMC for N=2562.Reuse & Permissions
    • Figure 4
      Figure 4
      Equation of state from ECMC and MPMC for N=10242. Error bars are mostly smaller than the symbols. Results agree within one standard deviation. The inset shows the relative pressure difference ΔP/P of MPMC with respect to ECMC.Reuse & Permissions
    • Figure 5
      Figure 5
      Orientational order parameter field ψ(x) of configurations obtained with the MPMC algorithm for system size N=10242. With increasing density, (a) pure liquid (η=0.700), (b) a bubble of hexatic phase (η=0.704), and (c) a stripe regime of hexatic phase (η=0.708) are visible. The interface between the liquid and the hexatic phase is extremely rough. (d) A scale bar illustrates the size of the fluctuations. The phase of ψ is represented via the color wheel.Reuse & Permissions
    • Figure 6
      Figure 6
      Positional order parameter field χ(x) of configurations obtained with (a), (b) the MPMC and (c), (d) ECMC algorithms for system size N=10242. (a), (c) In the hexatic phase (η=0.718), positional order is short-range. (b), (d) Toward higher density (η=0.720), fluctuations are much weaker and bounded, as expected for a continuous transition to a solid phase. The scale bar and the color code for the phase of χ are identical to the scale bar and the color code in Fig. 5.Reuse & Permissions
    • Figure 7
      Figure 7
      The positional correlation function Cq0(r) shows exponential decay at density η=0.718 and approaches a power law ∝r−1/3 at η=0.720. (a) System size N=5122: Excellent agreement between ECMC and MPMC. (b) System size N=10242: Excellent agreement in the hexatic phase (η=0.718) and fair agreement at the approach of the solid phase (η=0.720). Our algorithms fall out of strict equilibrium in the solid and long-scale correlations become sensitive to the boundary conditions.Reuse & Permissions
    • Figure 8
      Figure 8
      Influence of the configuration averaging on the positional configuration functions in the hexatic phase. For single configurations (n=1, data incoherently averaged), positional correlations cannot decay below a level given by the square root of the ratio of the correlated domain size to the system size. Coherent averaging over n=2, 4, 8, 16, and 32 configurations reduces the noise level and the correlations. The data corresponds to a single long MPMC run with N=10242 particles at density η=0.718.Reuse & Permissions
    ×

    Sign up to receive regular email alerts from Physical Review E

    Log In

    Cancel
    ×

    Search

    Article Lookup

    Paste a citation or DOI
    Enter a citation
    ×