Authors: Patrick Butler, Simon D. Rihm, Sebastian Mosbach, Jethro Akroyd, and Markus Kraft*

• Computational methods benchmarked against 85 experimental MOPs
• Machine learning methods accurately and efficiently optimise MOP geometries
• Optimised MOP structures and properties integrated into The World Avatar
• Demonstrated host-guest screening of optimised MOPs for urea binding

Metal–organic polyhedra (MOPs) have exceptional potential for host–guest chemistry, but their discovery is hindered by the large combinatorial space of their building units. Computational screening offers a powerful workflow for efficiently screening large sets of MOPs; however, reliable results depend on accurate modelling of the geometries and properties of the MOPs. In this paper, we extend our previous workflow, which assembled computation-ready MOPs using purely geometric operations, by incorporating post-assembly computational modelling. We benchmark a range of methods, including machine learning interatomic potentials (MLIPs) and tight-binding DFT, against 85 experimentally resolved MOP structures. The results show that geometry optimisation significantly refines the initial assembled MOP structures in terms of cavity and pore properties when compared against experimental structures, with MLIPs found to achieve good accuracy and excellent convergence rates. We applied the most reliable method to the entire dataset, integrating the resulting data within The World Avatar through the OntoMOPs ontology and enabling natural language querying. Lastly, we demonstrate the utility of this refined dataset by screening for MOPs with the potential to act as hosts for a urea guest molecule.

• This paper draws from preprint 341: Computational Modelling Workflows for Metal-Organic Polyhedra in The World Avatar
• Access the article at the publisher: DOI: 10.1039/D5CP04871K

Associated Themes: