Collective Variables for Conformational Polymorphism in Molecular Crystals

Abstract

Controlling polymorphism in molecular crystals is crucial in the pharmaceutical, dye, and pesticide industries. However, its theoretical description is extremely challenging, due to the associated long timescales ($>1\mu s$). We present an efficient procedure for identifying collective variables that promote transitions between conformational polymorphs in molecular dynamics simulations. It involves applying a simple dimensionality reduction algorithm to data from short ($\sim ps$) simulations of the isolated conformers that correspond to each polymorph. We demonstrate the utility of our method in the challenging case of the important energetic material, CL-20, which has three anhydrous conformational polymorphs at ambient pressure. Using these collective variables in Metadynamics simulations, we observe transitions between all solid polymorphs in the biased trajectories. We reconstruct the free energy surface and identify previously unknown defect and intermediate forms in the transition from one known polymorph to another. Our method provides insights into complex conformational polymorphic transitions of flexible molecular crystals.

References in corpus (12)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• PLUMED 2: New feathers for an old bird
• Combining Machine Learning and Computational Chemistry for Predictive Insights Into Chemical Systems
• Caliber based spectral gap optimization of order parameters (SGOOP) for sampling complex molecular systems
• Data-Driven Collective Variables for Enhanced Sampling
• Neural networks-based variationally enhanced sampling
• Predicting polymorphism in molecular crystals using orientational entropy
• Navigating at Will on the Water Phase Diagram
• Finding Efficient Collective Variables: The Case of Crystallization
• Making the best of a bad situation: a multiscale approach to free energy calculation
• Blind search for complex chemical pathways using Harmonic Linear Discriminant Analysis
• CO2 packing polymorphism under pressure: mechanism and thermodynamics of the I-III polymorphic transition

Cited by in corpus (1)

• A unified framework for machine learning collective variables for enhanced sampling simulations: $\text{mlcolvar}$