High-dimensional neural network potentials for solvation: The case of protonated water clusters in helium

Abstract

The design of accurate helium-solute interaction potentials for the simulation of chemically complex molecules solvated in superfluid helium has long been a cumbersome task due to the rather weak but strongly anisotropic nature of the interactions. We show that this challenge can be met by using a combination of an effective pair potential for the He-He interactions and a flexible high-dimensional neural network potential (NNP) for describing the complex interaction between helium and the solute in a pairwise additive manner. This approach yields an excellent agreement with a mean absolute deviation as small as 0.04 kJ/mol for the interaction energy between helium and both, hydronium and Zundel cations compared to CCSD(T) reference calculations with an energetically converged basis set. The construction and improvement of the potential can be performed in a highly automated way, which opens the door for applications to a variety of reactive molecules to study the effect of solvation on the solute as well as the solute-induced structuring of the solvent. Furthermore, we show that this NNP approach yields very convincing agreement with the CCSD(T) reference for properties like many-body spatial and radial distribution functions. This holds for the microsolvation of the protonated water monomer and dimer by a few helium atoms up to their solvation in bulk helium as obtained from path integral simulations at about 1 K.

References in corpus (4)

• Separable Dual Space Gaussian Pseudo-potentials
• How van der Waals interactions determine the unique properties of water
• Efficient first-principles calculation of the quantum kinetic energy and momentum distribution of nuclei
• High Order Path Integrals Made Easy

Cited by in corpus (15)

• Towards Exact Molecular Dynamics Simulations with Machine-Learned Force Fields
• Committee neural network potentials control generalization errors and enable active learning
• Automated Fitting of Neural Network Potentials at Coupled Cluster Accuracy: Protonated Water Clusters as Testing Ground
• PiNN: A Python Library for Building Atomic Neural Networks of Molecules and Materials
• Path Integral Molecular Dynamics for Fermions: Alleviating the Sign Problem with the Bogoliubov Inequality
• Converged Colored Noise Path Integral Molecular Dynamics Study of the Zundel Cation down to Ultra-low Temperatures at Coupled Cluster Accuracy
• Reduced Rovibrational Coupling Cartesian Dynamics for Semiclassical Calculations: Application to the Spectrum of the Zundel Cation
• Perspective: Non-conventional Force Fields for Applications in Spectroscopy and Chemical Reaction Dynamics
• Nuclear Quantum Effects in Liquid Water Are Negligible for Structure but Significant for Dynamics
• Random sampling versus active learning algorithms for machine learning potentials of quantum liquid water
• Neural Network Interaction Potentials for para-Hydrogen with Flexible Molecules
• Machine-learning interatomic potentials achieving CCSD(T) accuracy for systems with extended covalent networks and van der Waals interactions
• When Theory Meets Experiment: What Does it Take to Accurately Predict ${}^1$H NMR Dipolar Relaxation Rates in Neat Liquid Water from Theory?
• The Homunculus Brain and Categorical Logic
• Machine-Learned Potentials for Solvation Modeling