2018/11/30 by Natalie P. Schieber, Michael R. Shirts · 7 citations
Materials Science · Physics and Astronomy · #Block Copolymer Self-Assembly #Cutoff #Dispersion (optics) #Interpolation (computer graphics) #Line (geometry) #Material Dynamics and Properties #Molecular dynamics #Phase (matter) #Phase diagram #Phase space #Solidification and crystal growth phenomena #cond-mat.stat-mech #physics.comp-ph #van der Waals force
paper · pdf · doi:10.1063/1.5080431
published in The Journal of Chemical Physics 150(16), 164112 (American Institute of Physics) · main paper: 11 pages with 7 figures; supplementary material: 3 pages wirth 5 figures
openalex created_date 2018/11/16 · openalex publication_date 2019/04/24 · arxiv created 2019/05/03 · arxiv updated 2019/05/06 · openalex updated_date 2026/08/05
In this study, we incorporate configuration mapping between simulation ensembles into the successive interpolation of multistate reweighting (SIMR) method in order to increase phase space overlap between neighboring simulation ensembles. This significantly increases computational efficiency over the original SIMR method in many situations. We use this approach to determine the coexistence curve of face-centered cubic-hexagonal close-packed Lennard-Jones spheres using direct molecular dynamics and SIMR. As previously noted, the coexistence curve is highly sensitive to the treatment of the van der Waals cutoff. Using a cutoff treatment, the chemical potential difference between phases is moderate and SIMR quickly finds the phase equilibrium lines with good statistical uncertainty. Using a smoothed cutoff results in nonphysical errors in the phase diagram, while the use of particle mesh Ewald for the dispersion term results in a phase equilibrium curve that is comparable with previous results. The drastically closer free energy surfaces for this case test the limits of this configuration mapping approach to phase diagram prediction.