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.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the explicit treatment of hydrogen bonds in the Universal Force Field (UFF). Specifically, the author proposes a method to handle hydrogen bonds by introducing an explicit bond with an almost negligible bond order, rather than treating them as non - bonding interactions (such as van der Waals forces and Coulomb forces). This method aims to maintain the structural integrity of hydrogen - bonded complexes and avoid the need to assign arbitrary charges to the system. The paper mainly focuses on the application effect of this new method in Metal - Organic Frameworks (MOFs), hydrates and host - guest complexes. ### Main problems and solutions 1. **Problem**: When dealing with hydrogen bonds, traditional UFF and other classical force fields usually consider hydrogen bonds as non - bonding interactions, which requires the definition of atomic charges. However, defining charges is not practical in some cases, especially when dealing with host - guest systems that are completely described by force fields. In addition, non - bonding interactions may lead to large forces during the initial geometric optimization, making it difficult to converge to the expected local minimum. 2. **Solutions**: - **Introducing explicit hydrogen bonds**: The author proposes a new method to handle hydrogen bonds by introducing an explicit bond with a very small bond order (e.g., 0.001). This method increases the coordination number of the receptor site, but does not have a significant impact on fully coordinated atoms. - **Avoiding charge definition**: This method does not require the definition of atomic charges, thus avoiding the problem of unclear charge definition and reducing the computational cost. - **Improving the accuracy of structure prediction**: By explicitly handling hydrogen bonds, the author shows that the maximum error in predicting unit parameters is reduced from 66% to 14%, and the mean absolute error is reduced from 14% to 4%. ### Verification and application - **Verification**: The author uses the hydrogen - bonded complexes in the S22 dataset for verification, and the results show that the explicit hydrogen bond method performs well in predicting hydrogen bond distances. - **Application**: This method has been applied to a series of MOFs, hydrates and host - guest complexes, including the narrow - pore structure of MIL - 53(Al), water ice, methane hydrate and some crystal structures containing hydrogen bonds. The results show that the explicit hydrogen bond method performs excellently in the structure prediction of these systems and significantly improves the prediction accuracy. ### Conclusion The method proposed by the author provides a simple and effective way to handle hydrogen bonds, which not only avoids the problem of charge definition, but also significantly improves the accuracy of structure prediction. This is of great significance for quickly generating high - quality initial structures, especially when performing molecular dynamics simulations and electronic structure calculations.

Explicit treatment of hydrogen bonds in the Universal Force Field: validation and application for Metal-Organic Frameworks, hydrates and host-guest complexes

Accurate and Compatible Force Fields for Molecular Oxygen, Nitrogen, and Hydrogen to Simulate Gases, Electrolytes, and Heterogeneous Interfaces

Explicit implementation of hydrogen transport in metals

Chapter 5 Treatment of Hydrogen Bonds in Protein Simulations

Hydrogen storage in MOF-5: A van der Waals density functional theory study

Density Functionals for Hydrogen Storage: Defining the H2Bind275 Test Set with Ab Initio Benchmarks and Assessment of 55 Functionals

Explicitly correlated intermolecular distances and interaction energies of hydrogen bonded complexes

Automatic learning of hydrogen-bond fixes in an AMBER RNA force field

Role of Kinetic Exchange and Coulomb Interaction in Bonding of Hydrogen Molecular Systems and Excited States

An analytical bond-order potential for the copper–hydrogen binary system

Universal Hydrogen Bond Symmetrisation Dynamics Under Extreme Conditions

First-principles study of solid hydrogen: Comparison among four exchange-correlation functionals

Theoretical and experimental analysis of H2 binding in a prototype metal organic framework material

Hydrogen-bonded frameworks for molecular structure determination

Quantum Monte Carlo Calculations of Dihydrogen Binding Energetics on Ca Cations: an Assessment of Errors in Density Functionals for Weakly Bonded Systems

Improved association in a classical density functional theory for water

Predicting Hydrogen Storage Performances in Porous Aromatic Frameworks Containing Carboxylate Functional Groups with Divalent Metallic Cations

Flexible Force Field Parameterization through Fitting on the Ab Initio-Derived Elastic Tensor

Proposed Hydrogen-Bonding Index of Donor or Acceptor Reflecting Its Intrinsic Contribution to Hydrogen-Bonding Strength

Predicting Hydrogen Adsorption Uptake in Metal-Organic Framework (mof) Containing Unsaturated Coordinate Coppers Using Molecular Simulation

Ab Initio Study of Structural and Mechanical Property of Solid Molecular Hydrogens