Using Implicit Solvent With Neural Network Potentials

by Corin Wagen · Oct 7, 2025

Neural network potentials (NNPs) are typically trained to reproduce the results of gas-phase density-functional-theory calculations. Unfortunately, chemistry is weird in the gas phase and many reactions don't really exist on the gas-phase potential-energy surface, making it tough to use NNPs to study many important chemical processes.

For instance, here's a fairly iconic reaction: the addition of methanethiolate to methyl vinyl ketone. These sorts of thia-Michael additions are fairly common in covalent inhibitors, and we might hope that a simple ground-state two-electron process would be well-described by conventional computational methods. Unfortunately, running a scan for this process in the gas phase shows a confused potential-energy surface without a particularly obvious barrier to thiolate addition.

A thia-Michael reaction run with DFT.

View this calculation on Rowan.

The fundamental problem is not with the specific computational method but with the fact that this modeling is done in the gas phase: the unsolvated anions are so unstable and reactive that this addition becomes almost barrierless. Without the ability to add solvent corrections, even the best NNP will be incapable of usefully modeling this chemistry.

One solution to this problem is to train a separate NNP on implicit-solvent DFT calculations. This has been done before—Roman Zubatyuk and co-workers used transfer learning to build a version of AIMNet that predicted energies with implicit water—but necessitates training a new model for every solvent under study, which isn't ideal.

Instead, we've built machinery here at Rowan to separately compute the energy of solvation using semiempirical quantum chemistry. Using the analytical linearized Poisson–Boltzmann (ALPB) implicit-solvent model in combination with Grimme's GFN2-xTB semiempirical method, we're able to quickly add solvent corrections to calculations from any NNP. While this makes calculations slightly slower (since we have to run xTB calculations in addition to standard NNP calculations), the results are much improved for many reactive systems.

Here's the same thia-Michael addition with ALPB solvation. As the strength of the solvent increases from gas-phase (purple) to hexane (blue) to dichloromethane (green) to water (red), we can see that the barrier becomes more defined and the reaction becomes more clearly downhill, as expected.

A thia-Michael reaction run with different implicit-solvent models.

View these scans on Rowan.

Adding implicit solvent to NNPs can also be useful even in cases where the underlying PES is described qualitatively correctly in the gas phase. The prototypical SN2 reaction between chloride and methyl chloride has a low barrier of around 14 kcal/mol, but addition of an implicit-solvent water model produces a more reasonable barrier of about 23 kcal/mol. (Note that these barriers are estimated from the scan graph, and that actual thermochemical calculations should be performed to get quantitative barrier-height estimates.)

A Sn2 reaction run with different implicit-solvent models.

View these scans on Rowan.

A few caveats are worth noting here. This approach relies on GFN2-xTB/ALPB getting the underlying chemistry approximately correct, and so it might not work in every case. Also, the requisite xTB calculations require matrix diagonalization and thus can become slow for large systems—this isn't as scalable as a pure NNP-based approach. Still, we've found that the ALPB-based approach is quite effective and allows us to model a lot of properties and reactions with much-improved accuracy.

If you want to try combining NNPs with ALPB yourself, it's very simple to do through Rowan. Simply select any NNP from Rowan and then choose ALPB from the solvent menu:

Selecting ALPB implicit water for an NNP in Rowan.

How to add ALPB to an NNP in Rowan.

That's it! We automatically take care of adding implicit-solvent corrections for geometry optimization, single-point energies, and frequency calculations. If this is of interest to your research, you can try this out for free today by making a Rowan account.

Banner background image

What to Read Next

Using Implicit Solvent With Neural Network Potentials

Using Implicit Solvent With Neural Network Potentials

Modeling polar two-electron reactivity accurately with neural network potentials trained on gas-phase DFT.
Oct 7, 2025 · Corin Wagen
Preparing SMILES for Downstream Applications

Preparing SMILES for Downstream Applications

How to quickly use Rowan to predict the correct protomer and tautomer for a given SMILES.
Oct 3, 2025 · Corin Wagen
Better Search and Filtering

Better Search and Filtering

the problem of too many calculations; new ways to search, filter, and sort; how to access these tools; future directions
Sep 30, 2025 · Ari Wagen and Spencer Schneider
Boltz-2 Constraints, Implicit Solvent for NNPs, and More

Boltz-2 Constraints, Implicit Solvent for NNPs, and More

new terms of service; comparing IRCs and conformer searches; contact and pocket constraints for Boltz-2; MOL2 download; implicit-solvent NNPs; draft workflows; optimizing docking efficiency
Sep 22, 2025 · Corin Wagen, Ari Wagen, Jonathon Vandezande, Eli Mann, and Spencer Schneider
Controlling the Speed of Rowan's Docking

Controlling the Speed of Rowan's Docking

Some notes on how docking can be tuned for different applications.
Sep 22, 2025 · Corin Wagen
Studying Scaling in Electron-Affinity Predictions

Studying Scaling in Electron-Affinity Predictions

Testing low-cost computational methods to see if they get the expected scaling effects right.
Sep 10, 2025 · Corin Wagen
Open-Source Projects We Wish Existed

Open-Source Projects We Wish Existed

The lacunæ we've identified in computational chemistry and suggestions for future work.
Sep 9, 2025 · Corin Wagen, Jonathon Vandezande, Ari Wagen, and Eli Mann
How to Make a Great Open-Source Scientific Project

How to Make a Great Open-Source Scientific Project

Guidelines for building great open-source scientific-software projects.
Sep 9, 2025 · Jonathon Vandezande
ML Models for Aqueous Solubility, NNP-Predicted Redox Potentials, and More

ML Models for Aqueous Solubility, NNP-Predicted Redox Potentials, and More

the promise & peril of solubility prediction; our approach and models; pH-dependent solubility; testing NNPs for redox potentials; benchmarking opt. methods + NNPs; an FSM case study; intern farewell
Sep 5, 2025 · Eli Mann, Corin Wagen, and Ari Wagen
Machine-Learning Methods for pH-Dependent Aqueous-Solubility Prediction

Machine-Learning Methods for pH-Dependent Aqueous-Solubility Prediction

Prediction of aqueous solubility for unseen organic molecules remains an outstanding and important challenge in computational drug design.
Sep 5, 2025 · Elias L. Mann, Corin C. Wagen