Simulation of melting of two dimensional Lennard-Jones solids

Abstract

We study the nature of melting of a two dimensional (2D) Lennard-Jones solid using large scale Monte Carlo simulation. We use systems of up to 102,400 particles to capture the decay of the correlation functions associated with translational order (TO) as well as the bond-orientational (BO) order. We study the role of dislocations and disclinations and their distribution functions. We computed the temperature dependence of the second moment of the two order parameters, i.e., of the TO order parameter as well as of the order parameter associated with BO order and using finite-size scaling we determined the corresponding two anomalous dimension critical exponents eta. We also computed the temperature dependent distribution of these two order parameters on the complex plane which support a two stage melting with a hexatic phase as an intermediate phase. The analysis of our results leads to a consistent picture strongly supporting a two stage melting scenario as predicted by the Kosterlitz, Thouless, Halperin, Nelson, and Young (KTHNY) theory. We find that the two correlation lengths associated with translational and bond-orientational order have a distinctly different temperature dependence each diverging at different temperature and that their finite size scaling properties are consistent with the KTHNY theory. We also used the temperature dependence of etas and their theoretical bounds to provide estimates for the two critical temperatures, which are found to be in agreement with those determined from the divergent correlation lengths and from the Binder ratio.

14 two-column revtex pages and 20 figures

Cited by in corpus (23)

• Identification of structure in condensed matter with the topological cluster classification
• Stone-Wales Defects Preserve Hyperuniformity in Amorphous Two-Dimensional Materials
• Specific heat in two-dimensional melting
• Effect of a potential softness on the solid-liquid transition in a two-dimensional core-softened potential system
• Attraction tames two-dimensional melting: from continuous to discontinuous transitions
• Can random pinning change the melting scenario of two-dimensional core-softened potential system?
• How Close to Two Dimensions Does a Lennard-Jones System Need to Be to Produce a Hexatic Phase?
• Multiple Phase Transitions in Extended Hard Core Lattice Gas Models in Two Dimensions
• Possible quantum liquid crystal phases of helium monolayers
• Phase behaviour of Lennard-Jones particles in two dimensions
• Random pinning elucidates the nature of melting transition in two-dimensional core-softened potential system
• Nearly hyperuniform, nonhyperuniform, and antihyperuniform density fluctuations in two-dimensional transition metal dichalcogenides with defects
• Roles of energy dissipation in a liquid-solid transition of out-of-equilibrium systems
• Phase diagram of the system with the repulsive shoulder potential in two dimensions: density functional approach
• Accurate determination of the translational correlation function of two-dimensional solids
• The High Density Phase of the $k$-NN Hard Core Lattice Gas Model
• Liquid-Hextic-Solid Phase Transition of a Hard-Core Lattice Gas with Third Neighbor Exclusion
• On the connection between sliding friction and phonon lifetimes: Thermostat induced thermolubricity effects in molecular dynamics simulations
• The melting of the classical two dimensional Wigner crystal
• Thermal melting of a vortex lattice in a quasi two-dimensional Bose gas
• The Liquid--Hexatic Transition for Soft Disks
• Collective dynamics and shattering of disturbed two-dimensional Lennard-Jones crystals
• Topological Defects, Inherent Structures, and Hyperuniformity