🤖 AI Summary
This study addresses the challenge of predicting irreversible phase transitions in real-world systems by proposing an equation-free dynamical systems analysis framework based on Local Neural Operators (LNOs). Methodologically, it introduces the first integration of spatiotemporal LNOs into Krylov subspace iterative methods, enabling multiscale fixed-point computation, stability assessment, and bifurcation analysis without requiring explicit governing differential equations. Unlike conventional surrogate models, the framework supports rigorous system-level analysis and directly identifies critical phenomena. Evaluated on three canonical nonlinear PDE benchmarks—Allen–Cahn, Liouville–Bratu–Gelfand, and FitzHugh–Nagumo—the approach successfully detects multiple pitchfork, Hopf, and saddle-node bifurcations. It significantly accelerates computer-aided modeling and analysis of spatiotemporal dynamics while preserving mathematical fidelity and interpretability.
📝 Abstract
Neural Operators (NOs) provide a powerful framework for computations involving physical laws that can be modelled by (integro-) partial differential equations (PDEs), directly learning maps between infinite-dimensional function spaces that bypass both the explicit equation identification and their subsequent numerical solving. Still, NOs have so far primarily been employed to explore the dynamical behavior as surrogates of brute-force temporal simulations/predictions. Their potential for systematic rigorous numerical system-level tasks, such as fixed-point, stability, and bifurcation analysis - crucial for predicting irreversible transitions in real-world phenomena - remains largely unexplored. Toward this aim, inspired by the Equation-Free multiscale framework, we propose and implement a framework that integrates (local) NOs with advanced iterative numerical methods in the Krylov subspace, so as to perform efficient system-level stability and bifurcation analysis of large-scale dynamical systems. Beyond fixed point, stability, and bifurcation analysis enabled by local in time NOs, we also demonstrate the usefulness of local in space as well as in space-time ("patch") NOs in accelerating the computer-aided analysis of spatiotemporal dynamics. We illustrate our framework via three nonlinear PDE benchmarks: the 1D Allen-Cahn equation, which undergoes multiple concatenated pitchfork bifurcations; the Liouville-Bratu-Gelfand PDE, which features a saddle-node tipping point; and the FitzHugh-Nagumo (FHN) model, consisting of two coupled PDEs that exhibit both Hopf and saddle-node bifurcations.