Hydrogen-Bond-Basicity Prediction

Understanding and controlling the strength of hydrogen bonds is key to modulating a molecule's interactions and physicochemical properties. Rowan's hydrogen-bond-basicity prediction workflow uses neural network potentials, low-cost quantum chemical methods, and molecular electrostatic potential calculations to efficiently and robustly predict the strength of different hydrogen-bond acceptors.

How it works

Rowan's hydrogen-bond-basicity workflow begins with an openconf search on the input molecule. The resulting conformers are optimized and scored with the AIMNet2 neural network potential, then geometrically clustered to remove redundant structures. The lowest-energy retained conformer is used for the electrostatic-potential calculation.

We then locate electrostatic potential minima in the region of each lone pair of the lowest energy conformer, using the low-cost r²SCAN-3c composite density-functional-theory method. pK_BHX values are predicted by linearly scaling V_min, the value of the molecular electrostatic potential at the minimum.

To find values for the scaling constants, we performed Levenberg–Marquardt least-squares fitting to a database of experimentally measured pK_BHX values, predicting the measured per-molecule pK_BHX values by combining the predicted pK_BHX values from each distinct electrostatic potential minimum (up to three per hydrogen-bond acceptor). Here are the scaling constants found through this procedure:

For a more detailed discussion of outliers and potential sources of error, see the preprint.