Predicting molecular crystal growth: Automated surface terminations, interchangeable interaction energy models and finite particles in the CrystalClear protocol
- 1School of Molecular and Life Sciences, Curtin University, PO Box U1987, Perth, Western Australia, Australia, 6845
- 2Department of Chemistry, School of Natural Sciences, University of Manchester, Oxford Road, Manchester, UK, M13 9PL
- 3Crystogen, 46 Grafton Street, Manchester, United Kingdom, M13 9NT
Abstract
The incorporation of crystal growth processes extends the utility of predictive techniques beyond simply crystal geometry and a thermodynamic ranking with the ultimate goal being actionable information to target specific crystallisation. Pairwise model interaction energies with implicit solvation have been previously demonstrated in the CrystalClear protocol to provide a computationally efficient route to incorporate such processes, though some limitations restricted its use. We present several advancements to overcome these limitations, incorporating a native implementation of the semi-empirical GFN2-xTB model with SMD solvation, and validating both GFN2-xTB and the CE-1p model against reference lattice energies and solubilities. Further, automated enumeration of symmetry-unique faces and construction of solvated surface energies as part of the same protocol facilitates direct prediction of solvent-dependent equilibrium crystal shapes while eliminating the tedium and errors frequently associated with manual construction of surface geometries required to evaluate surface free energies. To demonstrate the potential of these advancements we present two case studies, oxindole and ritonavir, evaluating the accuracy of the crystal growth protocol in reproducing experimental results for particle size effects, polymorph stability and more. All methods presented here are freely available and implemented in open source software, ensuring accessibility, extendability and transferability to the broader community.
Related
- Predicting Adsorbate-Solvent Interaction Energies in Zeolite Pores Using Convolutional Neural Networks with Attention Mechanisms.
- Benchmarking a fast energy model for predicting the stability of solvated surface sites for organic molecular crystals: The paracetamol–water interface
- Large integration time-step in molecular dynamics simulation artificially enhances the strength of hydrophobic interaction
- The Molecular Bases of the Interaction between a Saponin from the Roots of Gypsophila paniculata L. and Model Lipid Membranes
- Molecular Docking and Molecular Dynamics Studies: Evaluation of the Potential Interaction between Ciprofloxacin-DAB-PAMAM Dendrimer-Based Prodrug Prototypes and Human Carboxylesterase 1
BibTeX
@misc{brookes2026predicting,
title = {{Predicting molecular crystal growth: Automated surface terminations, interchangeable interaction energy models and finite particles in the CrystalClear protocol}},
author = {James R. Brookes and Peter R. Spackman and Julian D. Gale and Michael W. Anderson and Alvin J. Walisinghe},
howpublished = {ChemRxiv},
year = {2026},
month = oct,
doi = {10.26434/chemrxiv.15011118/v1},
url = {https://chemrxiv.org/doi/full/10.26434/chemrxiv.15011118/v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Predicting molecular crystal growth: Automated surface terminations, interchangeable interaction energy models and finite particles in the CrystalClear protocol
AU - James R. Brookes
AU - Peter R. Spackman
AU - Julian D. Gale
AU - Michael W. Anderson
AU - Alvin J. Walisinghe
PY - 2026
DA - 2026/10/27
PB - ChemRxiv
DO - 10.26434/chemrxiv.15011118/v1
UR - https://chemrxiv.org/doi/full/10.26434/chemrxiv.15011118/v1
AB - The incorporation of crystal growth processes extends the utility of predictive techniques beyond simply crystal geometry and a thermodynamic ranking with the ultimate goal being actionable information to target specific crystallisation. Pairwise model interaction energies with implicit solvation have been previously demonstrated in the CrystalClear protocol to provide a computationally efficient route to incorporate such processes, though some limitations restricted its use. We present several advancements to overcome these limitations, incorporating a native implementation of the semi-empirical GFN2-xTB model with SMD solvation, and validating both GFN2-xTB and the CE-1p model against reference lattice energies and solubilities. Further, automated enumeration of symmetry-unique faces and construction of solvated surface energies as part of the same protocol facilitates direct prediction of solvent-dependent equilibrium crystal shapes while eliminating the tedium and errors frequently associated with manual construction of surface geometries required to evaluate surface free energies. To demonstrate the potential of these advancements we present two case studies, oxindole and ritonavir, evaluating the accuracy of the crystal growth protocol in reproducing experimental results for particle size effects, polymorph stability and more. All methods presented here are freely available and implemented in open source software, ensuring accessibility, extendability and transferability to the broader community.
ER -