Collective Variables for Conformational Polymorphism in Molecular Crystals

Oren Elishav,Roy Podgaetsky,Olga Meikler,Barak Hirshberg

DOI: https://doi.org/10.1021/acs.jpclett.2c03491

2022-11-15

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.

Chemical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the challenge of controlling conformational polymorphism in molecular crystals. Specifically, the authors focus on the polymorphism of molecular crystals, which is very important in the pharmaceutical, dye, and pesticide industries. The polymorphism of these materials has a significant impact on their chemical and physical properties. Therefore, being able to effectively control and understand the transformation mechanisms between these polymorphs is crucial for practical applications.

However, it is extremely difficult to theoretically describe these polymorphic phenomena, mainly because these processes usually occur on a very long time scale (> 1 microsecond), which is beyond the capabilities of standard simulation methods. To solve this problem, the authors propose an efficient method to identify collective variables (CVs) that can promote the transformation between conformational polymorphs. This method determines these CVs by applying a simple dimensionality reduction algorithm to the short - time (~ picosecond) simulation data of the isolated conformations corresponding to each polymorph.

The authors pay special attention to CL - 20 (Hexanitrohexaazaisowurtzitane), which is an important energetic material and has three anhydrous conformational polymorphs at normal pressure. Using these CVs for metadynamics simulations, the authors observed the transformation between all solid - state polymorphs and reconstructed the free - energy surface, and discovered previously unknown defects and intermediate forms during the transformation from one known polymorph to another. This study not only provides new insights into the complex conformational polymorph transformation but also provides a theoretical basis for designing experimental methods to isolate the desired polymorphs.