Abstract:Stiff systems of ordinary differential equations (ODEs) arise in a wide range of scientific and engineering disciplines and are traditionally solved using implicit integration methods due to their stability and efficiency. However, these methods are computationally expensive, particularly for applications requiring repeated integration, such as parameter estimation, Bayesian inference, neural ODEs, physics-informed neural networks, and MeshGraphNets. Explicit exponential integration methods have been proposed as a potential alternative, leveraging the matrix exponential to address stiffness without requiring nonlinear solvers. This study evaluates several state-of-the-art explicit single-step exponential schemes against classical implicit methods on benchmark stiff ODE problems, analyzing their accuracy, stability, and scalability with step size. Despite their initial appeal, our results reveal that explicit exponential methods significantly lag behind implicit schemes in accuracy and scalability for stiff ODEs. The backward Euler method consistently outperformed higher-order exponential methods in accuracy at small step sizes, with none surpassing the accuracy of the first-order integrating factor Euler method. Exponential methods fail to improve upon first-order accuracy, revealing the integrating factor Euler method as the only reliable choice for repeated, inexpensive integration in applications such as neural ODEs and parameter estimation. This study exposes the limitations of explicit exponential methods and calls for the development of improved algorithms.

Dual stepsize explicit numerical integration method and applications

A Transient Stability Numerical Integration Algorithm for Variable Step Sizes Based on Virtual Input

Rethinking to Finite Difference Time-Step Integrations

Numerical Derivative-based Flexible Integration Algorithm for Power Electronic Systems Simulation Considering Nonlinear Components

A precise Runge–Kutta integration and its application for solving nonlinear dynamical systems

An Efficient Nondissipative Higher-Order Single-Step Integration Method for Long-Term Dynamics Simulation

On Enhanced Second-Order Explicit Integration Methods with Controllable Algorithmic Dissipation and Adjustable Sub-Step Size for Hyperbolic Problems

A Two-Step Time Integration Method with Desirable Stability for Nonlinear Structural Dynamics

Explicit Multi-Time Step Computing Method in the Form of Non-Integer Step Ratio

New Explicit Integration Algorithms with Controllable Numerical Dissipation for Structural Dynamics

POWER SYSTEM FULL DYNAMIC SIMULATION PARTI:NUMERICAL METHOD

PRECISE TIME-INTEGRATION WITH[1mm］ DIMENSION EXPANDING METHOD

Variable Step Euler Method for Real-Time Simulation

A Three-Stage Explicit Time Integration Method with Controllable Numerical Dissipation

Euler Methods with Different Step-Sizes for Different Equations for Stochastic Differential Equations

An Improved Sub-Step Composite Time Integration Formulation with Enhanced Performance on Linear and Nonlinear Dynamics

Highly Accurate and Efficient Time Integration Methods with Unconditional Stability and Flexible Numerical Dissipation

Highly Precise Time Integration Method For Linear Structural Dynamic Analysis

Numerical Solutions of Energy-Conserving Time Integration Methods

Performance Evaluation of Single-step Explicit Exponential Integration Methods on Stiff Ordinary Differential Equations

A Robust Time-Integration Algorithm for Solving Nonlinear Dynamic Problems with Large Rotations and Displacements