pith. sign in

arxiv: 2606.06728 · v2 · pith:236KYGU2new · submitted 2026-06-04 · 🧮 math.DS

Data-driven methods for computation of optimal linear response in high-dimensional dynamical systems

Pith reviewed 2026-06-27 23:01 UTC · model grok-4.3

classification 🧮 math.DS
keywords optimal linear responsetransfer operatorsKoopman operatorsdynamical systemsdata-driven methodsspectrum manipulationalmost-invariant setshigh-dimensional systems
0
0 comments X

The pith

Kernel-smoothed approximations of transfer operators enable data-driven optimization of infinitesimal perturbations to manipulate system spectra.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a method to identify the smallest change to a nonlinear dynamical system that produces a chosen effect on its long-term statistics, such as slowing the loss of memory in almost-cycles or shifting their oscillation rates. It builds kernel-smoothed estimates of the transfer and Koopman operators directly from trajectory data, even when observations are high-dimensional, then uses these estimates inside an optimization problem to locate the desired perturbation. The same construction yields vector fields that show how the perturbation alters the observed behavior. Applications to periodic, chaotic, and climate models produce perturbations that align with the stated goals.

Core claim

Kernel-smoothed approximations of the transfer and Koopman operators, built from trajectory data, are inserted into a spectral optimization problem whose solution is the optimal infinitesimal perturbation that realizes any prescribed manipulation of the spectrum, including increases in frequency or reductions in decay rate for eigenvalues tied to almost-cycles or almost-invariant sets.

What carries the argument

Kernel-smoothed approximations of the transfer and Koopman operators that turn spectral manipulation into a tractable optimization problem solved from data.

If this is right

  • The optimization can be solved to increase frequency or suppress decay of correlations for almost-cycles identified by the approximated operator.
  • Optimal-response vector fields constructed from the same data visualize the physical effect of the perturbation under any chosen observations.
  • The procedure scales to high-dimensional trajectory data, as shown by its application to an Earth system model for El Nino Southern Oscillation.
  • The resulting perturbations are nontrivial yet consistent with the target dynamical objectives in the tested periodic, chaotic, and climate examples.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the transfer from approximation to true system holds, the same workflow could locate targeted adjustments that lengthen or shorten specific climate oscillations.
  • The vector-field construction might be reused to interpret how a chosen perturbation affects statistics under observations different from those used to build the operator.
  • Repeating the procedure on systems whose exact linear response is known analytically would directly test how much approximation error is tolerable before the computed perturbation loses its optimality.

Load-bearing premise

The kernel-smoothed operator approximations remain accurate enough that perturbations optimal for the approximation remain optimal for the true nonlinear system.

What would settle it

Take the perturbation found from the approximated operators, apply it to the original dynamical system, and measure whether the relevant eigenvalues of the true transfer operator move exactly as the optimization predicted.

Figures

Figures reproduced from arXiv: 2606.06728 by Dimitrios Giannakis, Gary Froyland, Nicholas Peters.

Figure 1
Figure 1. Figure 1: Left: Heat map of a 50 × 50 diagonal block of the 1000 × 1000 transition matrix P0. Right: Spectrum of P0, visualized in terms of the magnitudes, |λ (k) |, and frequencies ν (±k) = arg λ (±k)/2π∆t, associated with the eigenvalues λ (±k) . ν (±k) = arg λ (±k)/2π ∆t, are displayed in [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Left: Heat map of a 50×50 diagonal block of the 1000×1000 perturbation matrix P˙ . Right: The optimal-response vector field (red arrows), for the frequency￾increasing perturbation, originating on the corresponding observational data points. The circle rotation direction is anticlockwise. occurs on the geometric circle in R 2 . We emphasise that there is no a priori con￾straint on the optimal-response vecto… view at source ↗
Figure 3
Figure 3. Figure 3: Left: The variable rotation angle α(θ) (black line) plotted as a function of the angle θ. Also plotted are the magnitudes of the orthogonal projections of the optimal-response vector fields onto the tangent spaces of the instantaneous velocities in R 2 , for the frequency-increasing perturbation (red) and magnitude-increasing per￾turbation (green). Centre: Magnitudes and frequencies of the eigenvalues of P… view at source ↗
Figure 4
Figure 4. Figure 4: Left: The observational data is coloured according to the forward differ￾ence α(θ). The arrows represent the optimal-response vector field of the frequency￾increasing perturbation at a selection of points on the circle. Right: The same infor￾mation for the magnitude-increasing perturbation. The scaling factor and number of arrows were chosen independently to aid in visualisation, and are not the same betwe… view at source ↗
Figure 5
Figure 5. Figure 5: Magnitudes and frequencies of the leading eigenvalues of [PITH_FULL_IMAGE:figures/full_fig_p018_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Left: The observational data is coloured according to the argument of v0. The arrows represent the optimal-response vector field of the frequency perturba￾tion. For visibility, arrow lengths have been uniformly scaled and only a random sample of arrows are drawn, sampled with weighting inversely proportional to the density of points. Right: Entries of v0 are plotted in the complex plane, coloured by their … view at source ↗
Figure 7
Figure 7. Figure 7: Left: The observational data is coloured according to the magnitude of v0. The arrows represent the optimal-response vector field of the magnitude perturbation. Arrow selection and length scaling are as in [PITH_FULL_IMAGE:figures/full_fig_p020_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The magnitudes and frequencies of the leading eigenvalues of [PITH_FULL_IMAGE:figures/full_fig_p022_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The observational data is coloured according to the eigenvector [PITH_FULL_IMAGE:figures/full_fig_p023_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Left: The magnitudes and frequencies of the largest 30 eigenvalues of P0, approximating those of the transfer operator for the SST data. Right: The real and imaginary parts of the ENSO eigenvector. The magnitude and angle of each point can be interpreted respectively as the ENSO strength and phase during each month of the data. The cycle is broken up into 8 phases (polar sectors) each spanning π/4 radians… view at source ↗
Figure 11
Figure 11. Figure 11: Left: Phase composites of SST anomalies (colours) and surface wind anomalies (arrows) based on the ENSO eigenvector v0. Each arrow represents the mean wind anomaly in the surrounding 7.5 ◦×7.5 ◦ square. Center: Phase composites of the optimal-response vector fields for the SST field (colours) and surface wind field (arrows), for the ENSO frequency-increasing perturbation. Right: The same information for t… view at source ↗
Figure 12
Figure 12. Figure 12: SST anomalies (black line) and optimal-response SST fields for the [PITH_FULL_IMAGE:figures/full_fig_p029_12.png] view at source ↗
read the original abstract

We develop a data-driven framework for estimating optimal linear response of nonlinear dynamical systems. Our approach is based on kernel-smoothed approximations of the transfer/Koopman operators of the system, built from possibly high-dimensional observations along trajectories. Combining these operator approximations with the theory developed in [Antown et al. (2018), J. Stat. Phys., 170(6), 1051-1087], we formulate a computationally tractable optimization problem for the optimal infinitesimal perturbation realising any desired manipulation of the spectrum. We also introduce a notion of optimal-response vector fields for visualising, and physically interpreting, the system's response to the optimal perturbation under arbitrary observations. Our focus is on finding perturbations that optimally increase the frequency or optimally suppress the decay of correlations of almost-cycles or almost-invariant sets associated with the eigenvalues of the kernel-smoothed transfer operator. We illustrate our approach with applications to low-dimensional periodic and chaotic systems, as well as a high-dimensional example involving the El Nino Southern Oscillation in a comprehensive Earth system model. In these examples our approach discovers nontrivial optimal perturbations of the system, which are post hoc natural and consistent with the desired dynamical objectives.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper develops a data-driven framework for estimating optimal linear response in nonlinear dynamical systems. It uses kernel-smoothed approximations of transfer/Koopman operators constructed from trajectory data (possibly high-dimensional), combines these with the optimization theory of Antown et al. (2018) to pose a tractable problem for infinitesimal perturbations that achieve desired spectral manipulations (e.g., increasing frequency or suppressing correlation decay of almost-cycles), and introduces optimal-response vector fields for interpretation. The method is illustrated on low-dimensional periodic/chaotic systems and a high-dimensional ENSO example in an Earth system model, where the discovered perturbations are reported as nontrivial yet physically natural.

Significance. If the kernel approximations are shown to yield perturbations whose spectral effects transfer to the true nonlinear operator, the framework would offer a practical, observation-based route to optimal control of spectral properties in high-dimensional systems, with direct relevance to climate dynamics and chaos control. The combination of data-driven operator learning with an external optimization theory is a natural extension, and the introduction of optimal-response vector fields provides a useful interpretive tool.

major comments (3)
  1. [Formulation of the optimization problem] Formulation of the optimization problem (section referenced in the abstract): the kernel-smoothed operator is used to define the optimization problem, but no a-priori error bound or Lipschitz-type control is given on how the approximation error affects the non-convex optimizer; consequently it is not shown that a perturbation optimal for the smoothed operator remains near-optimal (or even produces the desired spectral shift) when the true transfer operator is substituted.
  2. [Numerical examples] Numerical examples section: the ENSO and low-dimensional examples report that the discovered perturbations are 'post hoc natural,' yet no quantitative metric (e.g., table of eigenvalue displacements or correlation-decay rates) compares the effect of the optimal vector field under the kernel-smoothed operator versus the true underlying dynamics or a finer-resolution reference operator.
  3. [Optimal-response vector fields] Definition of optimal-response vector fields: the construction appears to rely on the same kernel-smoothed operator used for the optimization; without an accompanying error analysis or sensitivity test, it is unclear whether these vector fields correctly represent the response of the original nonlinear system under arbitrary observations.
minor comments (2)
  1. Notation for the kernel bandwidth/smoothing parameter is introduced without an explicit symbol in the abstract and early sections; a consistent symbol and discussion of its selection would improve readability.
  2. The manuscript cites Antown et al. (2018) for the core optimization theory; a brief self-contained recap of the relevant theorem (or at least the precise statement used) would help readers who have not consulted the reference.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thorough and constructive report. We address each major comment below and outline the revisions we will make.

read point-by-point responses
  1. Referee: [Formulation of the optimization problem] the kernel-smoothed operator is used to define the optimization problem, but no a-priori error bound or Lipschitz-type control is given on how the approximation error affects the non-convex optimizer; consequently it is not shown that a perturbation optimal for the smoothed operator remains near-optimal (or even produces the desired spectral shift) when the true transfer operator is substituted.

    Authors: We agree that rigorous a-priori bounds on the propagation of kernel approximation error through the non-convex spectral optimization would be desirable. Deriving such bounds appears technically difficult because the problem is non-convex and the underlying operators are infinite-dimensional. In the revised manuscript we will add a dedicated discussion of this limitation together with numerical sensitivity experiments that vary kernel bandwidth, sample size, and regularization; these tests will quantify how much the computed optimal perturbations and resulting eigenvalue shifts change under controlled perturbations of the approximated operator. revision: partial

  2. Referee: [Numerical examples] the ENSO and low-dimensional examples report that the discovered perturbations are 'post hoc natural,' yet no quantitative metric (e.g., table of eigenvalue displacements or correlation-decay rates) compares the effect of the optimal vector field under the kernel-smoothed operator versus the true underlying dynamics or a finer-resolution reference operator.

    Authors: The observation is correct. For the low-dimensional periodic and chaotic examples we will insert tables that directly compare the eigenvalue displacements and correlation-decay rates obtained from the kernel-smoothed operator against the exact transfer operator (or a high-resolution reference). For the high-dimensional ENSO example, an exact operator is unavailable; we will instead report results obtained from several independent long trajectories and from cross-validation across different kernel parameters to demonstrate consistency of the discovered perturbations. revision: yes

  3. Referee: [Optimal-response vector fields] the construction appears to rely on the same kernel-smoothed operator used for the optimization; without an accompanying error analysis or sensitivity test, it is unclear whether these vector fields correctly represent the response of the original nonlinear system under arbitrary observations.

    Authors: We will revise the section on optimal-response vector fields to include explicit sensitivity tests with respect to kernel bandwidth and data subsampling. We will also add a clarifying paragraph stating that the vector fields are constructed from the data-driven approximation and therefore inherit its limitations; the tests will illustrate the degree of stability of the visualized fields under these variations. revision: yes

Circularity Check

0 steps flagged

No circularity; derivation combines external theory with data-driven approximations

full rationale

The paper constructs kernel-smoothed transfer/Koopman operator approximations from trajectory data and invokes the optimization framework of the external Antown et al. (2018) reference to set up a tractable problem for optimal perturbations. No step reduces a claimed result to a fitted quantity or self-defined input by construction, no load-bearing self-citation chain exists, and the derivation remains self-contained against the stated external benchmark and data-driven inputs. The transfer of optimality from the smoothed operator to the true system is an assumption separate from any definitional circularity.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The approach rests on the applicability of the 2018 theory to kernel approximations and on the choice of kernel and smoothing parameters, which are not specified as fixed or derived.

free parameters (1)
  • kernel bandwidth / smoothing parameter
    Kernel smoothing requires a bandwidth choice that affects the operator approximation; not stated as derived from first principles.
axioms (1)
  • domain assumption The theory developed in Antown et al. (2018) applies directly to the kernel-smoothed operator approximations.
    The optimization problem is formulated by combining the approximations with that external theory.
invented entities (1)
  • optimal-response vector fields no independent evidence
    purpose: Visualising and physically interpreting the system's response to the optimal perturbation under arbitrary observations.
    New interpretive object introduced in the paper; no independent evidence supplied.

pith-pipeline@v0.9.1-grok · 5738 in / 1298 out tokens · 23799 ms · 2026-06-27T23:01:11.235142+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

94 extracted references · 65 canonical work pages

  1. [1]

    Ghil and V

    M. Ghil and V. Lucarini. The physics of climate variability and climate change. Rev. Mod. Phys., 92:035002, 2020. doi:10.1103/RevModPhys.92.035002

  2. [2]

    Lucarini, F

    V. Lucarini, F. Ragone, and F. Lunkeit. Predicting climate change using re- sponse theory: Global averages and spatial patterns.J. Stat. Phys., 166:1036– 1064, 2017. doi:10.1007/s10955-016-1506-z

  3. [3]

    D. Ruelle. Differentiation of SRB states.Commun. Math. Phys., 187:227–241,

  4. [4]

    doi:10.1007/s002200050134

  5. [6]

    V. Baladi. Linear response, or else, 2014. URLhttps://arxiv.org/abs/1408

  6. [7]

    H. B. Callen and T. A. Welton. Irreversibility and generalized noise.Phys. Rev., 83(1):34–40, 1951. doi:10.1103/PhysRev.83.34

  7. [8]

    R. Kubo. The fluctuation–dissipation theorem.Rep. Prog. Phys., 29:255–284,

  8. [9]

    doi:10.1088/0034-4885/29/1/306

  9. [10]

    R. H. Kraichnan. Classical fluctuation–relaxation theorem.Phys. Rev., 113(5),

  10. [11]

    doi:10.1103/PhysRev.113.1181. 30

  11. [13]

    Gallavotti and E

    G. Gallavotti and E. G. D. Cohen. Dynamic ensembles in stationary states. Phys. Rev. Lett., 74(14):2694–2697, 1995. doi:10.1103/PhysRevLett.74.2694

  12. [14]

    Gallavotti

    G. Gallavotti. Onsager hypothesis: Onsager reciprocity and fluctuation–dissipation theorem.J. Stat. Phys., 84(5/6):899–925, 1996. doi:10.1007/bf02174123

  13. [15]

    L.-S. Young. What are SRB measures, and which dynamical systems have them?J. Stat. Phys., 108:733–754, 2002. doi:10.1023/a:1019762724717

  14. [16]

    Bowen and D

    R. Bowen and D. Ruelle. The ergodic theory of Axiom A flows.Inventiones Mathematicae, 29:181–202, 1975. doi:10.1007/bf01389848

  15. [17]

    Gou¨ ezel and C

    S. Gou¨ ezel and C. Liverani. Banach spaces adapted to Anosov systems.Ergod. Theory Dyn. Syst., 26:189–217, 2006. doi:10.1017/s0143385705000374

  16. [18]

    R. V. Abramov and A. Majda. New approximations and tests of linear fluctuation-response for chaotic nonlinear forced-dissipative dynamical systems. J. Nonlinear Sci., 18:303–341, 2008. doi:10.1007/s00332-007-9011-9

  17. [19]

    Baladi and D

    V. Baladi and D. Smania. Smooth deformations of piecewise expanding uni- modal maps.Discrete and Continuous Dynamical Systems, 23(3):685–703,

  18. [20]

    doi:10.3934/dcds.2009.23.685

  19. [21]

    Hairer and A

    M. Hairer and A. J. Majda. A simple framework to justify linear response theory.Nonlinearity, 23(4):909–922, 2010. doi:10.1088/0951-7715/23/4/008

  20. [22]

    Baladi, T

    V. Baladi, T. Kuna, and V. Lucarini. Linear and fractional response for the srb measure of smooth hyperbolic attractors and discontinuous observables. Nonlinearity, 30(3), 2017. doi:10.1088/1361-6544/aa5b13

  21. [23]

    Baladi and M

    V. Baladi and M. Todd. Linear response for intermittent maps.Communica- tions in Mathematical Physics, 347, 2016. doi:10.1007/s00220-016-2577-z

  22. [24]

    Wormell and G.A

    C.L. Wormell and G.A. Gottwald. On the validity of linear response theory in high-dimensional deterministic dynamical systems.Journal of Statistical Physics, 172, 2018. doi:10.1007/s10955-018-2106-x

  23. [25]

    G. A. Gottwald, C. Wormell, and J Wouters. On spurious detection of linear response and misuse of the fluctuation-dissipation theorem in fi- nite time series.Physica D Nonlinear Phenomena, 331:89–101, 2016. doi:10.1016/j.physd.2016.05.010

  24. [26]

    A rigor- ous computational approach to linear response.Nonlinearity, 31(3):1073–1109, 2018

    Wael Bahsoun, Stefano Galatolo, Isaia Nisoli, and Xiaolong Niu. A rigor- ous computational approach to linear response.Nonlinearity, 31(3):1073–1109, 2018. 31

  25. [27]

    Chandramoorthy and Q

    N. Chandramoorthy and Q. Wang. Efficient computation of linear response of chaotic attractors with one-dimensional unstable manifolds.SIAM J. Appl. Dyn. Syst., 21(2):735–781, 2022. doi:10.1137/21m1405599

  26. [28]

    A. Ni. Fast differentiation of hyperbolic chaos.Archive for Rational Mechanics and Analysis, 250(1), 2026

  27. [29]

    Froyland and M

    G. Froyland and M. Phalempin. Optimal linear response for Anosov diffeomor- phisms, 2025. URLhttps://arxiv.org/abs/2504.16532

  28. [30]

    T. Bell. Climate sensitivity from fluctuation dissipation: Some simple model tests.J. Atmos. Sci, 37(8):1700–1708, 1980. doi:10.1175/1520- 0469(1980)037<1700:csffds>2.0.co;2

  29. [31]

    A. S. Gritsoun. Fluctuation-dissipation theorem on attractors of atmo- spheric models.Russ. J. Numer. Anal. Math. Modelling, 16(2):115–133, 2001. doi:10.1515/rnam-2001-0203

  30. [32]

    A. S. Gritsoun, G. Branstator, and V. P. Dymnikov. Construction of the lin- ear response operator of an atmospheric general circulation model to small external forcing.Russ. J. Numer. Anal. Math. Modelling, 17(5):399–416, 2002. doi:10.1515/rnam-2002-0503

  31. [33]

    A. J. Majda, R. V. Abramov, and M. J. Grote.Information Theory and Stochas- tics for Multiscale Nonlinear Systems, volume 25 ofCRM Monograph Series. Americal Mathematical Society, Providence, 2005

  32. [34]

    A. J. Majda, R. Abramov, and B. Gershgorin. High skill in low-frequency cli- mate response through fluctuation dissipation theorems despite structural insta- bility.Proc. Natl. Acad. Sci, 107(2):581–586, 2010. doi:10.1073/pnas.091299710

  33. [35]

    Dellnitz and O

    M. Dellnitz and O. Junge. On the approximation of complicated dynamical be- havior.SIAM J. Numer. Anal., 36:491, 1999. doi:10.1137/s0036142996313002

  34. [36]

    Sch¨ utte, W

    Ch. Sch¨ utte, W. Huisinga, and P. Deuflhard. Transfer operator approach to conformational dynamics in biomolecular systems. In B. Fiedler, editor,Ergodic Theory, Analysis, and Efficient Simulation of Dynamical Systems, pages 191–

  35. [37]

    doi:10.1007/978-3-642-56589-2 9

    Springer-Verlag, Berlin, 2001. doi:10.1007/978-3-642-56589-2 9

  36. [38]

    I. Mezi´ c. Spectral properties of dynamical systems, model reduction and decom- positions.Nonlinear Dyn., 41:309–325, 2005. doi:10.1007/s11071-005-2824-x

  37. [39]

    Mezi´ c and A

    I. Mezi´ c and A. Banaszuk. Comparison of systems with complex behavior. Phys. D., 197:101–133, 2004. doi:10.1016/j.physd.2004.06.015

  38. [40]

    S. E. Otto and C. W. Rowley. Koopman operators for estimation and control of dynamical systems.Annu. Rev. Control Robot. Auton. Syst., 4:59–87, 2021. doi:10.1146/annurev-control-071020-010108. 32

  39. [41]

    Colbrook

    M. Colbrook. The multiverse of dynamic mode decomposition algorithms. InHandbook of Numerical Analysis, pages 127–230. Amsterdam, 2024. doi:10.1016/bs.hna.2024.05.004

  40. [42]

    S. L. Brunton, M. Budisi´ c, E. Kaiser, and J. N. Kutz. Modern Koop- man theory for dynamical systems.SIAM Rev., 64(2):229–340, 2022. doi:10.1137/21m1401243

  41. [43]

    Baladi.Dynamical Zeta Functions and Dynamical Determi- nants for Hyperbolic Maps: A Functional Approach

    V. Baladi.Dynamical Zeta Functions and Dynamical Determi- nants for Hyperbolic Maps: A Functional Approach. Springer, 2018. doi:https://doi.org/10.1007/978-3-319-77661-3

  42. [44]

    Santos Guti´ errez and V

    M. Santos Guti´ errez and V. Lucarini. On some aspects of the response to stochastic and deterministic forcings.J. Phys. A: Math. Theor., art. 425002,

  43. [45]

    doi:10.1088/1751-8121/ac90fd

  44. [46]

    Zagli, M

    N. Zagli, M. J. Colbrook, V. Lucarini, I. Mezi’c, and J. Moroney. Bridging the gap between Koopmanism and response theory: Using natural variability to predict forced response.SIAM J. Appl. Dyn. Syst., 25(1):196–229, 2026. doi:10.1137/24m1699206

  45. [47]

    Lucarini, M

    V. Lucarini, M. Santos Guti´ errez, J. Moroney, and N. Zagli. A general frame- work for linking free and forced fluctuations via Koopmanism.Chaos Solitons Fractals, 202(Part 1):117540, 2026. doi:10.1016/j.chaos.2025.117540

  46. [48]

    M. O. Williams, I. G. Kevrekidis, and C. W. Rowley. A data-driven approxi- mation of the Koopman operator: Extending dynamic mode decomposition.J. Nonlinear Sci., 25(6):1307–1346, 2015. doi:10.1007/s00332-015-9258-5

  47. [50]

    Froyland, C

    G. Froyland, C. Gonz´ alez-Tokman, and A. Quas. Detecting isolated spectrum of transfer and Koopman operators with Fourier analytic tools.J. Comput. Dyn., 1(2):249–278, 2014. doi:10.3934/jcd.2014.1.249

  48. [51]

    M. Porte. Linear response for Dirac observables of Anosov diffeomor- phisms.Discrete and Continuous Dynamical Systems, 39(4):1799–1819, 2019. doi:10.3934/dcds.2019078

  49. [52]

    Froyland

    G. Froyland. Statistically optimal almost-invariant sets.Phys. D., 200:205–219,

  50. [53]

    doi:10.1016/j.physd.2004.11.008

  51. [54]

    Froyland, D

    G. Froyland, D. Giannakis, B. Lintner, M. Pike, and J. Slawinska. Spectral analysis of climate dynamics with operator-theoretic approaches.Nat. Com- mun., 12:6570, 2021. doi:10.1038/s41467-021-26357-x

  52. [55]

    M. M. Castro and G. Froyland. On the structure of complex spectra and eigenfunctions of transfer and Koopman operators, 2025. URLhttps://arxiv. org/abs/2505.05770. 33

  53. [56]

    Froyland and N

    G. Froyland and N. Santitissadeekorn. Optimal mixing enhancement.SIAM Journal on Applied Mathematics, pages 1444–1470, 2017

  54. [57]

    Froyland, P

    G. Froyland, P. Koltai, and M. Stahn. Computation and optimal perturbation of finite-time coherent sets for aperiodic flows without trajectory integration. SIAM Journal on Applied Dynamical Systems, 19(3):1659–1700, 2020

  55. [58]

    Antown, D

    F. Antown, D. Dragiˇ cevi´ c, and G. Froyland. Optimal linear responses for Markov chains and stochastically perturbed dynamical systems.J. Stat. Phys., 170(6):1051–1087, 2018. doi:10.1007/s10955-018-1985-1

  56. [59]

    Antown, G

    F. Antown, G. Froyland, and S. Galatolo. Optimal linear response for Markov Hilbert–Schmidt integral operators and stochastic dynamical systems.J. Non- linear Sci., 32:79, 2022. doi:10.1007/s00332-022-09839-0

  57. [60]

    Santos Guti´ erez, N

    M. Santos Guti´ erez, N. Zagli, and G. Carigi. Markov matrix perturbations to optimize dynamical and entropy functionals, 2025. URLhttps://arxiv.org/ abs/2507.14040

  58. [61]

    Das and D

    S. Das and D. Giannakis. Delay-coordinate maps and the spectra of Koopman operators.J. Stat. Phys., 175(6):1107–1145, 2019. doi:10.1007/s10955-019- 02272-w

  59. [62]

    Giannakis

    D. Giannakis. Delay-coordinate maps, coherence, and approximate spectra of evolution operators.Res. Math. Sci., 8:8, 2021. doi:10.1007/s40687-020-00239- y

  60. [63]

    F. Takens. Detecting strange attractors in turbulence. InDynamical Systems and Turbulence, volume 898 ofLecture Notes in Mathematics, pages 366–381. Springer, Berlin, 1981. doi:10.1007/bfb0091924

  61. [65]

    P. J. Webster, A. Moore, J. Loschnigg, and M. Leban. Coupled ocean dy- namics in the Indian Ocean during 1997–98.Nature, 401:356–360, 1999. doi:10.1038/43848

  62. [66]

    Sauer, J

    T. Sauer, J. A. Yorke, and M. Casdagli. Embedology.J. Stat. Phys., 65(3–4): 579–616, 1991. doi:10.1007/bf01053745

  63. [67]

    J. C. Robinson. A topological delay embedding theorem for infinite-dimensional dynamical systems.Nonlinearity, 18(5):2135–2143, 2005. doi:10.1088/0951- 7715/18/5/013

  64. [68]

    J. C. Robinson. A topological time-delay embedding theorem for infinite- dimensional cocycle dynamical systems.Discrete Cont. Dyn. Syst. Ser. B, 9 (3&4):731–741, 2008. doi:10.3934/dcdsb.2008.9.731. 34

  65. [69]

    E. R. Deyle and G. Sugihara. Generalized theorems for nonlinear state space re- construction.PLoS ONE, 6(3):e18295, 2011. doi:10.1371/journal.pone.0018295

  66. [70]

    E. A. Nadaraya. On estimating regression.Theory Probab. Appl., 9(1):141–142,

  67. [71]

    G. S. Watson. Smooth regression analysis.Sankhya Ser. A, 26(4):359–372, 1964

  68. [72]

    E. M. Stein.Harmonic Analysis: Real-Variable Methods, Orthogonality, and Oscillatory Integrals, volume 43 ofPrinceton Mathematical Series. Princeton University Press, Princeton, 1993

  69. [73]

    Haj lasz, P

    P. Haj lasz, P. Koskela, and H. Tuominen. Sobolev embeddings, extensions and measure density condition.Journal of Functional Analysis, 254(5):1217–1234,

  70. [74]

    doi:10.1016/j.jfa.2007.11.020

    ISSN 0022-1236. doi:10.1016/j.jfa.2007.11.020

  71. [75]

    Giannakis and M

    D. Giannakis and M. J. Latifi-Jebelli. Kernel smoothing operators on thick open domains.Anal. Math., 2026. In press

  72. [76]

    Bj¨ orn, P

    J. Bj¨ orn, P. MacManus, and N. Shanmugalingam. Fat sets and pointwise bound- ary estimates forp-harmonic functions in metric spaces.J. Anal. Math., 85: 339–369, 2001. doi:10.1007/bf02788087

  73. [77]

    Canto, L

    J. Canto, L. Ihnatsyeva, J. Lehrb¨ ack, and A. V. V¨ ah¨ akangas. Capaci- ties and density conditions in metric spaces.Potential Anal., 62(2), 2025. doi:10.1007/s11118-024-10137-5

  74. [78]

    R. Z. Khas’minskii. The behavior of a self-oscillating system acted upon by slight noise.Journal of Applied Mathematics and Mechanics, 27(4):1035–1044,

  75. [79]

    doi:https://doi.org/10.1016/0021-8928(63)90184-9

    ISSN 0021-8928. doi:https://doi.org/10.1016/0021-8928(63)90184-9

  76. [80]

    Y. I. Kifer. On small random perturbations of some smooth dynamical systems.Math. USSR-Izv., 8:1083–1107, 1974. doi:10.1070/IM1974v008n05ABEH002139

  77. [81]

    S. M. Ulam.Problems in Modern Mathematics. Dover Publications, Mineola, 1964

  78. [82]

    J. G. Kemeny and J. L. Snell.Finite Markov Chains. Springer New York, 1983

  79. [83]

    Froyland, C

    G. Froyland, C. P. Rock, and K. Sakellariou. Sparse eigenbasis approximation: Multiple feature extraction across spatiotemporal scales with application to co- herent set identification.Communications in Nonlinear Science and Numerical Simulation, 77:81–107, 2019. doi:https://doi.org/10.1016/j.cnsns.2019.04.012

  80. [84]

    Froyland, D

    G. Froyland, D. Giannakis, E. Luna, and J. Slawinska. Revealing trends and persistent cycles of non-autonomous systems with operator-theoretic tech- niques: Applications to past and present climate dynamics.Nat. Commun., 15: 4268, 2024. doi:10.1038/s41467-024-48033-6. 35

Showing first 80 references.