Sparse identification of quasipotentials via a combined data-driven method

TL;DR

This paper introduces a novel data-driven method combining neural networks and symbolic regression to discover quasipotentials in nonlinear dynamical systems, aiding in understanding metastability and energy landscapes from data.

Contribution

It presents a new approach for directly deriving parsimonious quasipotential equations from data, integrating neural networks with symbolic regression for multistable systems.

Findings

Successfully applied to systems with known quasipotentials

Effective on noisy data and higher-dimensional systems

Provides analytical forms of quasipotentials for energy landscape analysis

Abstract

The quasipotential function allows for comprehension and prediction of the escape mechanisms from metastable states in nonlinear dynamical systems. This function acts as a natural extension of the potential function for non-gradient systems and it unveils important properties such as the maximum likelihood transition paths, transition rates and expected exit times of the system. Here, we demonstrate how to discover parsimonious equations for the quasipotential directly from data. Leveraging machine learning, we combine two existing data-driven techniques, namely a neural network and a sparse regression algorithm, specifically designed to symbolically describe multistable energy landscapes. First, we employ a vanilla neural network enhanced with a renormalization and rescaling procedure to achieve an orthogonal decomposition of the vector field. Next, we apply symbolic regression to extract the downhill and circulatory components of the decomposition, ensuring consistency with the underlying dynamics. This symbolic reconstruction involves a simultaneous regression that imposes constraints on both the orthogonality condition and the vector field. We implement and benchmark our approach using an archetypal model with a known exact quasipotential, as well as a nanomechanical resonator system. We further demonstrate its applicability to noisy data and to a four-dimensional system. Our model-unbiased analytical forms of the quasipotential is of interest to a wide range of applications aimed at assessing metastability and energy landscapes, serving to parametrically capture the distinctive fingerprint of the fluctuating dynamics.