Neural Network Matrix Product Operator: A Multi-Dimensionally Integrable Machine Learning Potential
Kentaro Hino, Yuki Kurashige
TL;DR
The paper presents NN-MPO, a neural-network-based PES model that leverages matrix product operator and tensor-train representations to enable accurate high-dimensional energy surfaces with a small training set. By mapping inputs to a latent space with a trainable Coordinator and expressing the weight tensor in TT form, NN-MPO achieves spectroscopic accuracy ($3.03\ \mathrm{cm}^{-1}$ MAE on a 6D H$_2$CO PES trained with 625 points) while permitting efficient evaluation of multi-dimensional integrals essential for quantum dynamics. The approach demonstrates compatibility with DMRG-based vibrational calculations, achieving $0.325\ \mathrm{cm}^{-1}$ MAE against exact MPO results for up to twenty excited states, underscoring its potential for quantum simulations beyond classical ML PES methods. The work also discusses practical considerations for scalability, symmetry, and coordinate choices, outlining pathways to extend NN-MPO to larger systems through localized coordinates and improved regularization. Overall, NN-MPO bridges traditional SOP/HDMR tensor representations and modern neural networks, offering a scalable, accurate framework for quantum molecular simulations and related high-dimensional integral tasks.
Abstract
A neural network-based machine learning potential energy surface (PES) expressed in a matrix product operator (NN-MPO) is proposed. The MPO form enables efficient evaluation of high-dimensional integrals that arise in solving the time-dependent and time-independent Schrödinger equation and effectively overcomes the so-called curse of dimensionality. This starkly contrasts with other neural network-based machine learning PES methods, such as multi-layer perceptrons (MLPs), where evaluating high-dimensional integrals is not straightforward due to the fully connected topology in their backbone architecture. Nevertheless, the NN-MPO retains the high representational capacity of neural networks. NN-MPO can achieve spectroscopic accuracy with a test mean absolute error (MAE) of 3.03 cm$^{-1}$ for a fully coupled six-dimensional ab initio PES, using only 625 training points distributed across a 0 to 17,000 cm$^{-1}$ energy range. Our Python implementation is available at https://github.com/KenHino/Pompon.
