Explicit treatment of hydrogen bonds in the universal force field: Validation and application for metal-organic frameworks, hydrates, and host-guest complexes

Abstract

A straightforward means to include explicit hydrogen bonds within the Universal Force Field is presented. Instead of treating hydrogen bonds as non-bonded interaction subjected to electrostatic and Lennard-Jones potentials, we introduce an explicit bond with negligible bond order, thus maintaining the structural integrity of the H-bonded complexes and avoiding the necessity to assign arbitrary charges to the system. The explicit hydrogen bond changes the coordination number of the acceptor site and the approach is thus most suitable for systems with under-coordinated atoms, such as many metal-organic frameworks, however, it also shows excellent performance for other systems involving a hydrogen-bonded framework. In particular, it is an excellent means for creating starting structures for molecular dynamics and for investigations employing more sophisticated methods.The approach is validated for the hydrogen bonded complexes in the S22 dataset and then employed for a set of metal-organic frameworks from the Computation-Ready Experimental (CoRE) database and several hydrogen bonded crystals including water ice and clathrates. We show that direct inclusion of hydrogen bonds reduces the maximum error in predicted cell parameters from 66% to only 14% and the mean unsigned error is similarly reduced from 14% to only 4%. We posit that with the inclusion of hydrogen bonding, the solvent-mediated breathing of frameworks such as MIL-53 is now accessible to rapid UFF calculations, which will further the aim of rapid computational scanning of metal-organic frameworks while providing better starting points for electronic structure calculations.

Proximity Effect in Crystalline Framework Materials: Stacking‐Induced Functionality in MOFs and COFs

“…In this work, the universal force field (UFF) [121] extended for MOFs [122,123] and hydrogen bonds were employed. [124] There are various other forcefield alternatives available for MOFs, which typically perform superior for a particular MOF subclass, but do not provide the almost complete range of elements that are covered by extended UFF. Some [125] can, however, produce superior structures that are comparable to those of quantum chemical approaches.…”

A Computational Protocol Combining DFT and Cheminformatics for Prediction of pH-Dependent Redox Potentials

“…This is due to a different treatment of electrostatic interactions (the main component of hydrogen bonding) in the two force fields. 27,28 The relationship between intramolecular hydrogen bonding and stability has been reported before. 10,29 Calculation of the redox potential at pH 0.…”

Proximity Effect in Crystalline Framework Materials: Stacking-Induced Functionality in MOFs and COFs

“…In this work we employed the Universal Force Field (UFF) [123] extended for MOFs [124,125] and hydrogen bonds. [126] There are various other forcefield alternatives available for MOFs, which typically perform superior for a particular MOF subclass, but do not provide the almost complete range of elements that are covered by extended UFF. Some [127] can, however, produce superior structures that are comparable to those of quantum chemical approaches.…”