Contrasting conformational behaviors of molecules XXXI and XXXII in the seventh blind test of crystal structure prediction

Accurate modeling of conformational energies is key to the crystal structure prediction of conformational polymorphs. Focusing on molecules XXXI and XXXII from the seventh blind test of crystal structure prediction, this study employs various electronic structure methods up to the level of domain-local pair natural orbital coupled cluster singles and doubles with perturbative triples [DLPNO-CCSD(T1)] to benchmark the conformational energies and to assess their impact on the crystal energy landscapes. Molecule XXXI proves to be a relatively straightforward case, with the conformational energies from generalized gradient approximation (GGA) functional B86bPBE-XDM changing only modestly when using more advanced density functionals such as PBE0-D4, ωB97M-V, and revDSD-PBEP86-D4, dispersion-corrected second-order Møller-Plesset perturbation theory (SCS-MP2D), or DLPNO-CCSD(T1). In contrast, the conformational energies of molecule XXXII prove difficult to determine reliably, and variations in the computed conformational energies appreciably impact the crystal energy landscape. Even high-level methods such as revDSD-PBEP86-D4 and SCS-MP2D exhibit significant disagreements with the DLPNO-CCSD(T1) benchmarks for molecule XXXII, highlighting the difficulty of predicting conformational energies for complex, drug-like molecules. The best-converged predicted crystal energy landscape obtained here for molecule XXXII disagrees significantly with what has been inferred about the solid-form landscape experimentally. The identified limitations of the calculations are probably insufficient to account for the discrepancies between theory and experiment on molecule XXXII, and further investigation of the experimental solid-form landscape would be valuable. Finally, assessment of several semi-empirical methods finds r²SCAN-3c to be the most promising, with conformational energy accuracy intermediate between the GGA and hybrid functionals and a low computational cost.