Search papers, labs, and topics across Lattice.
This paper introduces a computational framework combining correlated wavefunction theory (cWFT) and the many-body expansion to predict the lattice energies of complex molecular crystals with sub-chemical accuracy. The framework achieves accuracy within experimental uncertainties for 23 organic molecules and 13 ice polymorphs at a cost comparable to hybrid DFT, significantly improving efficiency compared to previous cWFT methods. The authors demonstrate the framework's applicability to large, drug-like molecules like axitinib and ROY, and use the generated reference data to parametrize a highly accurate DFT functional, both of which are made openly available.
Achieve sub-chemical accuracy in predicting the lattice energies of complex molecular crystals, including drug-like molecules, at DFT-level computational cost, unlocking high-throughput polymorph screening.
Molecules can form myriad crystalline polymorphs, each with distinct properties affecting their performance across diverse applications, from pharmaceuticals to functional materials and more. Predicting the thermodynamically most stable polymorph from first principles remains a formidable challenge. It requires methods that scale to large, technologically-relevant molecules while achieving very high accuracy (below 1 kJ/mol) on relative lattice energies. Such accuracy, often termed sub-chemical accuracy, is generally beyond the reach of the workhorse density functional theory (DFT). In this work, we introduce a framework, combining advances in correlated wavefunction theory (cWFT) and the many-body expansion, to deliver accurate, cost-effective predictions of complex molecular crystals. For 23 organic molecules and 13 ice polymorphs, we predict crystal lattice energies to within experimental uncertainties at costs comparable to hybrid DFT, while being several orders of magnitude more efficient than previous cWFT approaches. We extend this approach to a set of large, drug-like molecules including axitinib and ROY, previously inaccessible to cWFT and where DFT is insufficient, achieving sub-chemical accuracy on the relative energies between challenging polymorphs. With the reference data generated throughout this work, we have been able to further parametrize a DFT functional with unprecedented accuracy aligning with our predictions. This cWFT framework as well as DFT functional are made openly available, providing new ranking tools to facilitate efficient high-throughput screening of molecular crystal polymorphs.
