Calculation of phase diagrams in the multithermal-multibaric ensemble

Abstract

From the Ising model and the Lennard-Jones fluid, to water and the iron-carbon system, phase diagrams are an indispensable tool to understand phase equilibria. In spite of the effort of the simulation community the calculation of a large portion of a phase diagram using computer simulation is still today a significant challenge. Here we propose a method to calculate phase diagrams involving liquid and solid phases by the reversible transformation of the liquid and the solid. To this end we introduce an order parameter that breaks the rotational symmetry and we leverage our recently introduced method to sample the multithermal-multibaric ensemble. In this way in a single molecular dynamics simulation we are able to compute the liquid-solid coexistence line for entire regions of the temperature and pressure phase diagram. We apply our approach to the bcc-liquid phase diagram of sodium and the fcc-bcc-liquid phase diagram of aluminum.

arXiv comment9 pages, 4 figures

References in corpus (12)

• Canonical sampling through velocity-rescaling
• PLUMED 2: New feathers for an old bird
• On representing chemical environments
• Accurate determination of crystal structures based on averaged local bond order parameters
• Comparing molecules and solids across structural and alchemical space
• A Variational Approach to Enhanced Sampling and Free Energy Calculations
• Predicting polymorphism in molecular crystals using orientational entropy
• Enhancing entropy and enthalpy fluctuations to drive crystallization in atomistic simulations
• Solid-liquid interface free energy through metadynamics simulations
• Determining pressure-temperature phase diagrams of materials
• A local order metric for condensed phase environments
• Multithermal-multibaric molecular simulations from a variational principle

Cited by in corpus (17)

• Enhanced sampling methods for molecular dynamics simulations
• A Unified Approach to Enhanced Sampling
• Phase equilibrium of water with hexagonal and cubic ice using the SCAN functional
• Homogeneous ice nucleation in an ab initio machine learning model of water
• Liquid-Liquid Critical Point in Phosphorus
• Driving and characterizing nucleation of urea and glycine polymorphs in water
• Collective Variables for the Study of Crystallization
• Computer Simulation of the Early Stages of Self-Assembly and Thermal Decomposition of ZIF-8
• Phase equilibrium of liquid water and hexagonal ice from enhanced sampling molecular dynamics simulations
• Nucleating a Different Coordination in a Crystal under Pressure: A Study of the B1-B2 Transition in NaCl by Metadynamics
• Deep learning the slow modes for rare events sampling
• Enhancing the formation of ionic defects to study the ice Ih/XI transition with molecular dynamics simulations
• Phase Diagram of ZIF-4 Computed via Well-tempered Metadynamics
• Inverse Design Method with Enhanced Sampling for Complex Open Crystals: Application to Novel Zeolite Self-Assembly in a Coarse-Grained Model
• Melting curves of ice polymorphs in the vicinity of the liquid-liquid critical point
• Determination of melting temperature of hexagonal ice using Lee-Yang phase transition theory
• The Phase Diagram of Carbon Dioxide from Correlation Functions and a Many-body Potential