Learning Relativistic Geodesics and Chaotic Dynamics via Stabilized Lagrangian Neural Networks
Abdullah Umut Hamzaogullari, Arkadas Ozakin
TL;DR
This work addresses the instability of training Lagrangian Neural Networks (LNNs) by introducing Hessian-based regularization, gradient clipping, physics-aware scaling, and careful activation/initialization choices to enable learning of complex dynamics. The authors develop two regularization strategies—an eigenvalue penalty and a scale-invariant Sylvester-criterion approach—to enforce physically meaningful mass-matrix signatures for both classical and relativistic systems, including AdS$_4$ geodesics. They demonstrate substantial gains in stability and accuracy, achieving, for example, ~96% reductions in validation loss for double pendulum and enabling the learning of geodesic Lagrangians in non-relativistic and relativistic spacetimes, with potential to extract metric components from learned mass matrices. While limitations remain (notably the need for invertible Hessians and handling of coordinate singularities), the approach broadens LNN applicability to automated discovery of geometric structures and spacetime metrics from trajectory data, with broad implications for physics-informed machine learning and scientific discovery.
Abstract
Lagrangian Neural Networks (LNNs) can learn arbitrary Lagrangians from trajectory data, but their unusual optimization objective leads to significant training instabilities that limit their application to complex systems. We propose several improvements that address these fundamental challenges, namely, a Hessian regularization scheme that penalizes unphysical signatures in the Lagrangian's second derivatives with respect to velocities, preventing the network from learning unstable dynamics, activation functions that are better suited to the problem of learning Lagrangians, and a physics-aware coordinate scaling that improves stability. We systematically evaluate these techniques alongside previously proposed methods for improving stability. Our improved architecture successfully trains on systems of unprecedented complexity, including triple pendulums, and achieved 96.6\% lower validation loss value and 90.68\% better stability than baseline LNNs in double pendulum systems. With the improved framework, we show that our LNNs can learn Lagrangians representing geodesic motion in both non-relativistic and general relativistic settings. To deal with the relativistic setting, we extended our regularization to penalize violations of Lorentzian signatures, which allowed us to predict a geodesic Lagrangian under AdS\textsubscript{4} spacetime metric directly from trajectory data, which to our knowledge has not been done in the literature before. This opens new possibilities for automated discovery of geometric structures in physics, including extraction of spacetime metric tensor components from geodesic trajectories. While our approach inherits some limitations of the original LNN framework, particularly the requirement for invertible Hessians, it significantly expands the practical applicability of LNNs for scientific discovery tasks.
