A trajectory-driven algorithm for differentiating SRB measures on unstable manifolds

Abstract

SRB measures are limiting stationary distributions describing the statistical behavior of chaotic dynamical systems. Directional derivatives of SRB measure densities conditioned on unstable manifolds are critical in the sensitivity analysis of hyperbolic chaos. These derivatives, known as the SRB density gradients, are by-products of the regularization of Lebesgue integrals appearing in the original linear response expression. In this paper, we propose a novel trajectory-driven algorithm for computing the SRB density gradient defined for systems with high-dimensional unstable manifolds. We apply the concept of measure preservation together with the chain rule on smooth manifolds. Due to the recursive one-step nature of our derivations, the proposed procedure is memory-efficient and can be naturally integrated with existing Monte Carlo schemes widely used in computational chaotic dynamics. We numerically show the exponential convergence of our scheme, analyze the computational cost, and present its use in the context of Monte Carlo integration.

arXiv comment25 pages, 9 figures, submitted to journal, under review

References in corpus (14)

• Dynamical ensembles in stationary states
• A review of linear response theory for general differentiable dynamical systems
• Theory and computation of covariant Lyapunov vectors
• Fast numerical test of hyperbolic chaos
• Relaxation of finite perturbations: Beyond the Fluctuation-Response relation
• Linear response for macroscopic observables in high-dimensional systems
• Linear response, or else
• Non-stationary non-uniform hyperbolicity: SRB measures for dissipative maps
• A computable realization of Ruelle's formula for linear response of statistics in chaotic systems
• Computational assessment of smooth and rough parameter dependence of statistics in chaotic dynamical systems
• Linear response theory for diffeomorphisms with tangencies of stable and unstable manifolds. [A contribution to the Gallavotti-Cohen chaotic hypothesis.]
• Ergodic Sensitivity Analysis of One-Dimensional Chaotic Maps
• An ergodic averaging method to differentiate covariant Lyapunov vectors
• Differentiating densities on smooth manifolds

Cited by in corpus (3)

• On Some Aspects of the Response to Stochastic and Deterministic Forcings
• Heterogeneity of the Attractor of the Lorenz '96 Model: Lyapunov Analysis, Unstable Periodic Orbits, and Shadowing Properties
• Space-split algorithm for sensitivity analysis of discrete chaotic systems with unstable manifolds of arbitrary dimension