---
title: Lagrangian Descriptors in Dynamical Systems
url: https://www.emergentmind.com/topics/lagrangian-descriptors-lds
type: topic
---

# Lagrangian Descriptors in Dynamical Systems

Searching arXiv for recent and foundational papers on Lagrangian descriptors to ground the encyclopedia entry.
arXiv Search Query: "Lagrangian descriptors review invariant manifolds Hamiltonian"
Lagrangian descriptors (LDs) are finite-time, trajectory-based scalar diagnostics designed to expose the geometric skeleton of dynamical systems by accumulating a positive observable along orbits. Across the literature, LDs are used to reveal stable and unstable manifolds, separatrices, transport barriers, invariant tori, slow manifolds, chaotic layers, and bounded-motion regions by identifying sharp, non-smooth, or rapidly varying features in scalar fields over initial conditions [1705.10213]. Their scope now spans autonomous and nonautonomous flows, discrete maps, stochastic differential equations, dissipative systems, chemical reaction dynamics, fluid transport, celestial mechanics, and recent semiclassical and quantum generalizations [1106.1306].

## 1. Definitions and principal formulations

The original continuous-time formulation treats an LD as a finite-time accumulation of trajectory arc length. For a system
\[
\dot{\mathbf{x}}=\mathbf{v}(\mathbf{x},t),
\]
the basic symmetric descriptor is
\[
M(\mathbf{x}_0,t_0,\tau)=\int_{t_0-\tau}^{t_0+\tau}\|\mathbf{v}(\mathbf{x}(t;\mathbf{x}_0),t)\|\,dt,
\]
with forward and backward components \(M^+\) and \(M^-\) obtained by integrating on \([t_0,t_0+\tau]\) and \([t_0-\tau,t_0]\), respectively [1106.1306]. A closely related componentwise \(p\)-norm family, emphasized in the theoretical literature, is
\[
M_p(\mathbf{x}_0,t_0,\tau)=\int_{t_0-\tau}^{t_0+\tau}\sum_{i=1}^{n}\big|\dot{x}_i(t;\mathbf{x}_0)\big|^p\,dt,\qquad 0<p\le 1,
\]
or equivalently, in autonomous notation, by integrating \(\sum_i |f_i(\mathbf{x}(t),t)|^p\) along the orbit [1705.10213]. The choice \(0<p<1\) is repeatedly used because it sharpens singular features.

Discrete-time analogues replace the time integral by orbitwise sums of step increments. For maps, one standard definition is
\[
MD_p(x_0,N)=\sum_{i=-N}^{N-1}\|x_{i+1}-x_i\|_p,
\]
while componentwise forms
\[
\sum_{j=-T}^{T-1}\sum_{i=1}^{N}|z^{(i)}_{j+1}-z^{(i)}_j|^p
\]
are used in higher-dimensional symplectic maps [1705.11057]. This discrete viewpoint is central in the analysis of area-preserving maps and in later LD-based chaos diagnostics [2305.00978].

The observable accumulated by an LD is not unique. Beyond arc length, the literature includes acceleration-based, curvature-based, and higher time-derivative variants, as well as formulations that integrate phase-space coordinates and momenta directly rather than velocities. In the LiCN\(\leftrightarrow\)LiNC application, for example, the descriptor is
\[
M(\mathbf{z}_0,\tau)=\int_{-\tau}^{\tau}\sum_{i=1}^{4}|z_i(t)|^p\,dt,
\]
with \(p=0.4\), and \(\mathbf{z}\) denotes the phase-space coordinates themselves, not the velocity field [2110.09989]. In stochastic settings, the corresponding stochastic LD is defined along sample paths by summing path increments,
\[
MS_p(x_0,t_0,\tau,\omega)=\sum_{i=-N}^{N-1}\big\|X_{t_{i+1}}(\omega)-X_{t_i}(\omega)\big\|_p,
\]
so that randomness enters through the path rather than through an explicit diffusion term in the integrand [1705.11074].

## 2. Geometric mechanism and theoretical basis

The theoretical rationale for LDs is that invariant manifolds separate trajectories with qualitatively different forward or backward histories. When a positive quantity is accumulated along those trajectories, the resulting scalar field changes abruptly across the separating set. In the rigorous framework developed for \(M_p\), a “singular feature” is defined by nonexistence of a normal derivative of the LD across a surface, and such singular features are proven to coincide with stable and unstable manifolds in several two-dimensional hyperbolic settings [1705.10213].

For the linear autonomous saddle
\[
\dot{x}=\lambda x,\qquad \dot{y}=-\lambda y,\qquad \lambda>0,
\]
the descriptor takes the explicit form
\[
M_p((x_0,y_0),0,\tau)=2\,\frac{\lambda^{p-1}}{p}\sinh(\lambda p\tau)\,\big(|x_0|^p+|y_0|^p\big).
\]
For \(0<p<1\), derivatives such as \(\partial M_p/\partial x_0\propto \operatorname{sign}(x_0)|x_0|^{p-1}\) become unbounded at \(x_0=0\), and likewise at \(y_0=0\). The stable and unstable manifolds are therefore singular sets of \(M_p\) for any finite \(\tau\) [1705.10213]. Analogous results are established for nonlinear autonomous saddles via Moser normal form, and for linear and nonlinear nonautonomous hyperbolic systems via exact calculations or conjugacy arguments [1705.10213].

Forward and backward integrations isolate different invariant objects. In the standard hyperbolic picture, forward LDs emphasize stable manifolds and backward LDs emphasize unstable manifolds, while symmetric windows reveal both simultaneously [1705.10213]. In time-dependent chemical dynamics this same principle is used operationally: minima of \(M^+\) identify \(W^s\), minima of \(M^-\) identify \(W^u\), and their intersection defines the anchor for a non-recrossing dividing surface [1705.00248].

A complementary 1-DoF viewpoint removes time altogether. For planar Hamiltonian systems with one degree of freedom, the geometric LD is defined as the Euclidean arc length of an energy level curve,
\[
\mathrm{LD}_{\mathrm{geo}}(E)\equiv \ell(E)=\int_{\Gamma(E)} ds.
\]
In the pendulum, Duffing, and fish-tail separatrix problems, \(\ell(E)\) attains a cusp-like maximum at the separatrix energy, and its derivative diverges as
\[
\left|\frac{d\ell}{dE}\right|=\mathcal{O}\big(|E-E_{\mathrm{sep}}|^{-1/2}\big),
\]
which isolates the geometric mechanism behind LD singularities without reference to a finite time window [2112.05857].

## 3. Variants across deterministic, stochastic, and dissipative dynamics

The map-based extension preserves the same structural idea. For two-dimensional area-preserving autonomous and nonautonomous maps, explicit formulae show that discrete LDs with \(p<1\) are non-differentiable across stable and unstable manifolds of hyperbolic saddles. In linear autonomous maps of the form
\[
x_{n+1}=\lambda x_n,\qquad y_{n+1}=\lambda^{-1}y_n,\qquad \lambda>1,
\]
the resulting discrete LD is proportional to \(|x_0|^p+|y_0|^p\), so directional derivatives become unbounded along the invariant axes. Rotated saddles and nonlinear maps exhibit the same phenomenon, with singular directions converging to the true invariant manifolds as the iteration horizon increases [1705.11057].

For stochastic differential equations, LDs are interpreted within the random dynamical systems framework. The stochastic extension
\[
dX_t=b(X_t,t)\,dt+\sigma(X_t,t)\,dW_t
\]
uses sample-path increment sums \(MS_p\), optionally ensemble averaged over Wiener realizations. In the noisy saddle, stochastic Duffing oscillator, and stochastic double gyre, the resulting fields reveal random fixed points, random stable and unstable manifolds, and transport barriers. The mechanism remains the same: for \(0<p\le 1\), non-differentiability appears when the dominant forward or backward increment terms vanish on the relevant random manifolds [1705.11074].

Dissipative systems require additional care because forward and backward dynamics are strongly asymmetric and backward integrations may blow up in finite time. In that setting LDs have been used to detect hyperbolic saddles with unequal expansion and contraction rates, limit cycles, slow manifolds, strange attractors, and transition ellipsoids in damped Hamiltonian systems. The reported practice includes asymmetric forward and backward times to balance stable and unstable signatures, and variable-time stopping when trajectories escape a bounded computational region [2103.14719]. In the Hopf normal form, LDs distinguish the case \(\beta<0\), where the origin is a stable focus, from \(\beta>0\), where a stable limit cycle of radius \(r^*=\sqrt{\beta/\sigma}\) appears as a ring in the LD field [2103.14719]. In slow–fast settings such as the van der Pol oscillator and the rotating-hoop model, LD ridges align with slow invariant manifolds [2103.14719].

## 4. Computation, parameter choice, and numerical practice

The numerical workflow is usually minimal: choose a grid of initial conditions, integrate trajectories over a finite horizon, accumulate the chosen observable, and visualize the resulting scalar field. In continuous-time problems, the LD can be accumulated by augmenting the ODE with a scalar equation
\[
\dot{M}=g(\mathbf{x}(t),t),\qquad M(t_0)=0,
\]
and integrating the extended system [1106.1306]. In discrete-time problems, the procedure is an orbit sum of step lengths or increment norms [1705.11057].

Parameter selection is problem dependent, but several quantitative prescriptions recur. For hyperbolic structures, subunit exponents are favored because they sharpen LD singularities. Reported choices include \(p=0.4\) in LiCN isomerization, \(p=1/2\) in the double pendulum and in several dissipative examples, and \(p\in(0.1,0.5)\) in stochastic computations [2110.09989]. Integration time is crucial. In the LiCN study, the practical criterion
\[
\tau \gtrsim C\,|\lambda_u|^{-1},\qquad C=\mathcal{O}(5\text{–}10),
\]
is used to resolve manifolds associated with a periodic orbit whose unstable Floquet exponent is \(\lambda_u\). At \(E=4000\ \mathrm{cm}^{-1}\), the reported values are \(|\lambda_u|^{-1}\approx 1163\) a.u. for the transition-state periodic orbit and \(|\lambda_u|^{-1}\approx 2170\) a.u. for the unstable saddle-node branch; production maps use \(\tau=2\times 10^4\) a.u. [2110.09989].

Nonautonomous and stochastic applications impose additional numerical requirements. The stochastic LD paper uses Euler–Maruyama integration, two-sided Brownian motion when backward integration is needed, and ensemble averaging when single-realization structures are noisy [1705.11074]. In the double pendulum, forward LDs are computed with adaptive Dormand–Prince \(8(9)\), tolerance \(10^{-8}\), and \(\tau=700\) in dimensionless time units [2403.07000]. Recent work on discrete maps replaces finite-difference estimates of LD second derivatives by differential algebra, obtaining machine-precision values of
\[
\|\Delta \mathrm{LD}(x_0;N)\|=\sum_{i=1}^{d}\left|\frac{\partial^2}{\partial x_{0,i}^2}\mathrm{LD}(x_0;N)\right|,
\]
thereby improving threshold robustness in thin resonant webs and depleted meshes [2410.19435].

## 5. Applications in fluids, chemistry, astrodynamics, and transition-state dynamics

In geophysical fluid dynamics, LDs were introduced as global descriptors for aperiodically time-dependent flows and applied to altimetric ocean data. In the Kuroshio region, the symmetric arc-length LD identified distinguished trajectories, hyperbolic manifolds, and non-hyperbolic eddy cores on the same scalar map. For \(\tau=15\) and \(30\) days, singular features aligned with independently computed stable and unstable manifolds, while the loss of smoothness inside an eddy core between \(\tau=30\) and \(\tau=72\) days was interpreted as a finite confinement-time indicator, with \(2\tau\approx 60\) days giving an upper bound on core trapping [1003.0368].

In chemical reaction dynamics, LDs have become a diagnostic for invariant manifolds associated with reaction barriers. For the LiCN\(\leftrightarrow\)LiNC isomerization, the rotationless molecule is modeled by a 2-DoF Hamiltonian with two minima at \(\vartheta=0\) and \(\vartheta=\pi\), energies \(2281\ \mathrm{cm}^{-1}\) and \(0\ \mathrm{cm}^{-1}\), and a barrier top at \(E_{\mathrm{SP}}=3455\ \mathrm{cm}^{-1}\). On a Poincaré section built along the minimum-energy path, singular features of the LD locate the invariant manifolds of the transition-state periodic orbit and of the unstable orbit born at the saddle-node bifurcation \(E_{\mathrm{SN}_1}=3440.6\ \mathrm{cm}^{-1}\). The same framework supports a Morse-based equivalent 2-DoF potential and an adiabatic 1-DoF reduction that captures the principal separatrices and the emergence of a secondary dynamical barrier as radial excitation \(n_R\) increases [2110.09989].

In nonautonomous reaction dynamics, LDs are used constructively rather than just diagnostically. For a driven Gaussian barrier coupled to a bath mode, minima of forward and backward LDs on reactive-coordinate sections define \(W^s\) and \(W^u\); their intersection over all bath coordinates forms an anchor surface \(T(t)=W^s(t)\cap W^u(t)\). The corresponding time-dependent dividing surface is empirically non-recrossing. In the reported 2-DoF model, the LD-based dividing surface yields zero recrossings for all 160,000 test trajectories, and for a thermal ensemble at \(k_B T=5.13\) gives a rate \(k=3.03\), with Arrhenius parameters \(k_\infty=3.14\) and \(\Delta E_{\mathrm{eff}}=0.135\) [1705.00248].

Astrodynamics provides a different use case: finite-time boundedness rather than asymptotic invariant sets. In the perturbed planar bi-elliptic restricted four-body problem for the Didymos–Dimorphos system, forward LDs in the state \(x=[x,y,\dot{x},\dot{y}]^\top\) reveal phase-space structures that delimit bounded motion around Dimorphos over a finite horizon. Without solar radiation pressure, regions of bounded motion are visually identifiable; with SRP, the majority of those structures break down, leaving a large region of unstable motion with rare exceptions [2307.14943].

## 6. LD-based chaos diagnostics and comparison with FTLE-, SALI-, and Poincaré-based methods

Although LDs were introduced as geometric visualization tools, a substantial literature now uses them as quantitative chaos indicators. One route computes finite differences of LD values for neighboring initial conditions. In low-dimensional conservative systems, the Difference of Neighboring orbits’ Lagrangian Descriptors and Ratio of Neighboring orbits’ Lagrangian Descriptors,
\[
D_L^n,\qquad R_L^n,
\]
classify regular and chaotic motion with better than \(90\%\) agreement with SALI in both the Hénon–Heiles flow and the standard map, while using only short integration windows and forward LDs [2208.09853]. In the 4D symplectic standard map, the related indices \(D_L^n\), \(R_L^n\), and
\[
S_L^n(z)=\frac{1}{n}\sum_{i=1}^{n}\left|\frac{LD(y_i^+)-2LD(z)+LD(y_i^-)}{(\sigma^{(i)})^2}\right|
\]
achieve \(P_A\gtrsim 90\%\) against SALI for a balanced choice \(T=1000\), \(n=2\), \(\sigma=10^{-3}\), with \(D_L^n\) and \(S_L^n\) generally outperforming \(R_L^n\) [2305.00978].

A simpler recent proposal uses only one neighboring orbit:
\[
\Delta L(x_0;\beta)=|L(x_0)-L(x_0+\beta)|.
\]
The heuristic argument is that \(\Delta L\) grows linearly with the integration horizon in regular regions and essentially exponentially in chaotic ones, reflecting the underlying separation law. In the reported benchmarks, this “difference LD” matches or surpasses more elaborate LD-based diagnostics. For Hénon–Heiles at \(H=1/15\), the F1-score of \(\Delta L\) against SALI is \(1.0000\), while several multi-neighbor LD indicators fail badly; for the standard map at \(K=1.5\), \(\Delta L\) attains an F1-score of \(0.9933\) [2506.16660].

Comparison with FTLEs is one of the oldest themes in the LD literature. In the periodically forced double gyre, forward and backward FTLE ridges approximate stable and unstable manifolds, but the paper reports that FTLE ridge locations can shift with \(\tau\), depend on neighbor selection, and cannot be computed at boundaries, whereas LD singular contours preserve manifold location as \(\tau\) increases and reveal KAM islands as smooth low-\(M\) basins [2103.15557]. In the ocean-flow study, LDs recover both hyperbolic and non-hyperbolic regions in a single field, while FTLEs require separate forward and backward computations and may exhibit “ghost” structures in transient settings [1003.0368]. In deterministic-chaos applications, the second-derivative LD index \(DLD_L\) is presented as a non-variational alternative to FLI, MEGNO, and SALI, with about \(90\%\) agreement against SALI on 4D coupled standard maps [2309.10643]. The differential-algebra refinement of \(\|\Delta \mathrm{LD}\|\) goes further: finite-difference implementations can misclassify up to \(20\%\) of initial conditions in thin resonant webs, whereas differential algebra yields clean bimodal distributions and sensitivity comparable to MEGNO, FLI, and SALI across regimes [2410.19435].

## 7. Conceptual issues, controversies, and emerging directions

A recurrent conceptual issue is objectivity. The theoretical framework explicitly argues that LDs are not objective under time-dependent Euclidean frame changes and that this is not necessarily a defect: in nonautonomous dynamics, different frames can carry genuinely different phase portraits, so a strictly objective scalar would be unable to distinguish them [1705.10213]. This stands in contrast to diagnostics designed to be frame invariant. A related point is that some applied papers still describe LDs as heuristic and note that their precise mathematical connection to invariant manifolds remains debated outside the classes where rigorous results are available [2307.14943]. The present literature therefore supports two statements simultaneously: rigorous manifold-detection theorems exist for several hyperbolic settings, but general interpretive claims should still be tied to the dynamical context and to validation against independent structure calculations where possible [1705.10213].

Recent work extends LDs beyond classical deterministic mechanics. In the path-integral formulation, a quantum LD is defined as the fluctuation average of the classical LD around extremal trajectories. For the 1-DoF Hamiltonian saddle
\[
H(q,p)=\frac{\lambda}{2}(p^2-q^2),
\]
stable and unstable manifolds \(p=\pm q\) broaden into finite-width structures under quantum fluctuations. With a mode cutoff \(N\), time window \(T\), and Planck constant restored, the reported width estimate is
\[
\sigma_{\mathrm{rms}}=\sqrt{\frac{N\hbar}{4T\lambda}},
\]
and Monte Carlo sampling agrees with this prediction within \(1\%\) across the explored range [2604.04128]. A related Bohmian construction moves from a single wavefunction to a two-parameter preparation space of Gaussian wavepackets labeled by \((q_0,p_0)\). For the inverted harmonic oscillator, the wavepacket-center dynamics are exactly classical, and the preparation-space stability matrix yields an \(\mathcal{O}(e^{\omega T})\) bound on LD sensitivity, suggesting a geometric counterpart to the exponential growth familiar from semiclassical OTOC analyses [2603.20803].

These developments suggest that LDs have evolved from a finite-time flow-visualization tool into a broader geometric framework for transport, instability, and phase-space organization. The common thread remains unchanged: a scalar field built from trajectories can encode qualitative dynamical transitions as singularity patterns, whether the setting is a Kuroshio eddy, a chemical transition state, a stochastic transport barrier, a resonance web in a symplectic map, or a broadened quantum separatrix [1106.1306].

Source: https://www.emergentmind.com/topics/lagrangian-descriptors-lds