Quantifying surface free energy of molecular crystal β-HMX using non-equilibrium molecular dynamics methods

Abstract

Granular molecular crystals show a large variation in the orientation, length, shape, and cohesive interactions of interfaces when embedded in a polymer matrix. But the determination of the associated excess of free energy is not straightforward, especially due to the entropic contribution which is not negligible for molecular compounds. The surface free energy (SFE) is also important to predict crystal shapes and growth or to prove crucial insights into the mesoscale interfaces in the granular composites often with interfacial defects, voids, and pores. In this paper, we use Molecular Dynamics (MD) and assess and adapt two non-equilibrium methods, namely, non-equilibrium thermodynamic integration (NETI) and modified steered MD (SMD), to determine the SFE of various facets of a high-energy molecular crystal of interest: β-HMX. Starting with defect-free surfaces, both methods agree well with experiment concerning the most stable facets, whose energy is further lowered by a large entropic contribution. For some facets, surface creation (debonding process) and annihilation (bonding process) show fundamentally different paths, henceforth irreversible. Then, we discuss some sources of discrepancy between theoretical and experimental SFE, considering non-reversible dissipative events, surface roughness, and differences between various experimental techniques. These non-equilibrium methods can be applied to a wide range of molecular crystals and to study interface stability in polycrystals or with binding polymers, including the effects of various defects.