Sensitivity of solid phase stability to the interparticle potential range: studies of a new Lennard-Jones like model

Olivia Sandra Moro,Vincent Ballenegger,Tom Underwood,Nigel B. Wilding
DOI: https://doi.org/10.1039/d3cp05474h
IF: 3.3
2024-02-07
Physical Chemistry Chemical Physics
Abstract:In a recent article, Wang {\em et al} (Phys. Chem. Chem. Phys., 2020, {\bf 22}, 10624) introduced a new class of interparticle potential for molecular simulations. The potential is defined by a single range parameter, eliminating the need to decide how to truncate truly long-range interactions like the Lennard-Jones (LJ) potential. The authors explored the phase diagram for a particular value of the range parameter for which their potential is similar in shape to the LJ 12-6 potential. We have reevaluated the solid phase behaviour of this model using both Lattice Switch Monte Carlo and thermodynamic integration. In addition to finding that the boundary between hexagonal close packed (hcp) and face centred cubic (fcc) phases presented by Wang {\em et al} was calculated incorrectly, we show that owing to its finite range, the new potential exhibits several `artifact' reentrant transitions between hcp and fcc phases. The artifact phases, which do not occur in the full (untruncated) LJ system, are also found for typically adopted forms of the truncated and shifted LJ potential. However, whilst in the latter case one can systematically investigate and correct for the effects of the finite range on the calculated phase behaviour (a correction beyond the standard long-range mean field tail correction being required to eliminate the artifact phases), this is not possible for the new potential because the choice of range parameter affects the entire potential shape. Our results highlight that potentials with finite range may fail to represent the crystalline phase behavior of systems with long-range dispersion interactions, even qualitatively.
chemistry, physical,physics, atomic, molecular & chemical
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to re - evaluate the behavior of a new finite - range interaction potential (LJL potential) proposed by Wang et al. in the solid - state phase diagram. Specifically, the authors of the paper used the lattice - switching Monte Carlo (LSMC) and thermodynamic integration (TI) methods to study the behavior of this new potential function in the solid - state region, especially the transition boundary between the hexagonal close - packed (hcp) and face - centered cubic (fcc) phases. Through these methods, the authors found significant differences from the results previously reported by Wang et al., especially on the hcp - fcc phase - transition boundary. In addition, the paper also explored that when the particle number density is increased, the LJL potential exhibits multiple re - entering hcp - fcc transitions, which are absent in the complete Lennard - Jones (LJ) potential. The key findings of the paper include: 1. **Re - evaluation of the hcp - fcc phase - transition boundary**: The authors found that there were errors in the hcp - fcc phase - transition boundary calculated by Wang et al., and proposed a new phase - transition boundary. 2. **Re - entering hcp - fcc transitions**: As the particle number density increases, the LJL potential exhibits multiple re - entering hcp - fcc transitions, which are absent in the complete LJ potential. 3. **Temperature dependence**: These re - entering hcp - fcc transitions exist in a wide temperature range, even extending to T = 0, which is significantly different from the phase diagram of the complete LJ potential. 4. **Error in free - energy calculation**: The authors believe that there were errors in the free - energy calculation by Wang et al., which led to the phase diagrams they reported being inconsistent with the actual situation. Through these findings, the paper emphasizes the possible limitations of finite - range potential functions when simulating long - range dispersive interaction systems, especially in describing crystal - phase behaviors. These results are of great significance for understanding the applicability and limitations of finite - range potential functions in molecular simulations.