Generic and Isometric Embeddings in Reservoir Computers
Allen G Hart

TL;DR
This paper proves that reservoir computers can embed input system attractors in a way that preserves their structure, with explicit constructions for linear systems and conditions for high-dimensional reservoirs.
Contribution
It establishes conditions for generalized and isometric embeddings in reservoir systems, linking topological and metric preservation to reservoir dimension and linearity.
Findings
Reservoir systems can admit topological embeddings of input attractors.
High-dimensional reservoirs enable isometric embeddings.
Explicit isometric embeddings are constructed for linear systems.
Abstract
We prove that a generic reservoir system admits a generalized synchronization that is a topological embedding of the input system's attractor. We also prove that for sufficiently high reservoir dimension (given by Nash's embedding theorem) there exists an isometric embedding generalized synchronization. The isometric embedding can be constructed explicitly when the reservoir system and source dynamics are linear.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsNeural Networks and Reservoir Computing · Ferroelectric and Negative Capacitance Devices · Advanced Memory and Neural Computing
Generic and Isometric Embeddings in Reservoir Computers
Allen G Hart
(October 20, 2025)
Abstract
We prove that a generic reservoir system admits a generalized synchronization that is a topological embedding of the input system’s attractor. We also prove that for sufficiently high reservoir dimension (given by Nash’s embedding theorem) there exists an isometric embedding generalized synchronization. The isometric embedding can be constructed explicitly when the reservoir system and source dynamics are linear.
Lead Paragraph
Reservoir computing is a machine learning paradigm that uses a high-dimensional dynamical system to model and predict the behavior of an underlying time-dependent source system that is only partially observed. A key feature of this framework is the generalized synchronization, where the reservoir learns to mirror the evolution of the source system. In the special case that the generalized synchronization is an embedding, the synchronization can not only predict the future observations of the source system but also faithfully reconstructs the hidden topological structure of the underlying dynamics. In this work, we show that such embeddings occur generically in reservoir systems, providing rigorous justification for a phenomenon long observed in practice. Moreover, we establish conditions under which reservoirs can achieve isometric embeddings, preserving not just the topology but also the geometric properties such as lengths and angles. These results demonstrate that reservoirs have the capacity to represent dynamical systems without distortion, opening the way to more reliable learning, analysis, and discovery of invariants in nonlinear time series.
1 Introduction
Reservoir computing is a machine learning paradigm designed to efficiently process temporal data using a high-dimensional dynamical system called the reservoir, which transforms input signals into rich internal states. When a reservoir is driven by a sequence of observations from a dynamical system, the combined system can exhibit generalized synchronization (GS) Hart2021 ; Hart2020 ; HartThesis2022 ; Hart2024 ; Grigoryeva2021 ; Grigoryeva2023 ; HaraKokubu2022 ; Hu2022 ; Kobayashi2021 ; Lu2017 ; Lu2018 ; Lymburn2019Reservoir ; Nazerian2023 ; Wu2024 ; Verzelli2020 ; Chen2022 ; Ohkubo2024 ; Parlitz2024 ; Platt2021 ; Platt2022 ; Ruffini2018 where the reservoir state becomes functionally dependent on the input system’s state via a map , i.e., .
In the special case where the GS map is an embedding, the topological properties of the underlying dynamical system are preserved and faithfully replicated in the reservoir Hart2020 ; Verzelli2020 ; Pathak2017a . This enables not only the prediction of future observations, but also the recovery of geometric information about the original system such as the eigenvalues of fixed points, Lyapunov exponents, and homology groups Hart2020 ; Pathak2017a . Empirically, when generalized synchronization occurs, it often looks like an embedding. So in this paper we will prove that, under suitable conditions (inspired by Whitney Whitney1936 ) a generic reservoir map possesses a generalized synchronization which is an embedding. This result is closely related to the celebrated Takens embedding theorem Takens1981Detecting ; Grigoryeva2023 . We also prove that for sufficiently high reservoir dimension (given by the bounds of Nash’s embedding theorem) there exists an isometric embedding generalized synchronization. These embeddings additionally preserve the angles and lengths of the source dynamics.
2 Generic Embeddings
First, we will define a reservoir system.
Definition 2.1**.**
(Reservoir System) Let be a smooth dimensional manifold and a diffeomorphism on . Let be a scalar observation function. Let denote the reservoir map. Then the reservoir system is defined
[TABLE]
Under certain conditions, the reservoir system has a (not necessarily unique Grigoryeva2021 ,ceni2020echo ) associated generalised synchronisation.
Definition 2.2**.**
(Generalized Synchronization) Let be an open set and a reservoir system. The pair is said to have the generalized synchronization on the -invariant compact set if the graph
[TABLE]
is invariant under the evolution operator i.e. and and the graph is an attractor in the sense that there exists an such that for any
[TABLE]
Definition 2.3**.**
For a given open set and reservoir system let be the interior of the set of all reservoir maps such that has a GS .
Next we introduce Whitney’s 2nd embedding Theorem - a central result stating that if the reservoir dimension is greater than twice the manifold dimension then generic maps are embedding. This lower bound on the dimension will be used to prove the major result of this section; Theorem 2.5.
Theorem 2.4** (Whitney Whitney1936 ; Lee2013 ).**
Let be a smooth manifold of dimension . If , then a generic map in is an embedding. Equivalently, the subset of embeddings is open and dense.
Now we are ready to state and prove the major theorem - that for generic reservoir maps that have a GS, the associated GS is an embedding. The idea behind the proof is to define a continuous open mapping from the reservoir map to the GS and use the genericity of embeddings guaranteed by Whitney’s theorem to obtain the genericity of the maps . We will assume that the source system has the conditions necessary for Takens’ Theorem to hold.
Theorem 2.5**.**
Suppose that has finitely many periodic orbits in . Suppose that for each periodic orbit of with period , the derivative has distinct eigenvalues. Suppose that . Then for a generic observation function and generic reservoir map , the generalized synchronization is an embedding of .
Proof.
Let
[TABLE]
The set of generalised synchronisations is open and nonempty by Lemma 6.1. Furthermore, the set of embeddings is an open and dense subset by Whitney’s embedding Theorem Whitney1936 ; Lee2013 . Now it follows (by Lemma 6.2) that is a dense open subset of . Now define the map as the mapping of the reservoir map to its associated GS . Now using that is continuous and open (Lemma 6.3), and Lemma 6.4, it follows that is a generic subset of , completing the proof. ∎
Remark 2.6**.**
The periodic-orbit non-degeneracy assumption (distinct eigenvalues for each ) mirrors the standard hypotheses in Takens-type embedding theorems, ensuring transversality of delayed coordinates (see Huke’s formulation of Takens’ Theorem). The assumption is not very restrictive - and is in fact satisfied for generic evolution operators . This was established independently by Kupka Kupka1963 Kupka1963 and Smale1963 Smale. In fact the original Takens embedding theorem Takens1981Detecting did not spell out the periodic-orbit non-degeneracy assumptions in detail - and was shown by Takens to hold for a generic pair .
Remark 2.7**.**
Theorem 2.5 is closely related to the result in Grigoryeva2023 , which shows almost all linear reservoir systems reservoir produce a generalized synchronization that is generically an embedding of the attractor into the reservoir state space Grigoryeva2023 . We can view this result as a sort of special case of Theorem 2.5 which holds for generally nonlinear systems.
Remark 2.8**.**
The original Takens theorem Takens1981Detecting is set on a compact manifold , but the proof goes through using the same arguments when we are considering the restriction of the delay map to a compact attractor . Moreover, the dimension of the reservoir space can be bounded below by some integer lower than if the attractor has lower dimension than the manifold dimension. This sharper bound is outside the scope of this paper but is analysed in sauer1991embedology .
3 Discussion
Theorem 2.5 states that if the reservoir dimension is larger than twice the manifold dimension (following Whitney), then for a generic reservoir map , the associated GS is an embedding of the attractor . Hence, if we take a random sample from a non-singular distribution over (or an open subset of ), then we will almost surely obtain a reservoir map whose GS is an embedding.
In the reservoir computing paradigm, we typically generate a reservoir map by randomly drawing weights from a parameterized class. If this class has the universal approximation property—as is the case for Echo State Networks (ESNs)—then it has been observed empirically that, under suitable hyperparameter choices, almost all realizations of the random weights yield an embedding GS Haluszczynski2019 ; Gauthier2021 ; Wikner2021 . This empirical fact is explained precisely by Theorem 2.5, together with the universal approximation property of ESNs GrigoryevaOrtega2018 .
The results in this section are set in discrete time, so it is natural to wonder if they hold in continuous time as well. The arguments in this paper directly use Takens’ embedding theorem, which does not have a direct continuous-time analogue, so the problem is more complicated than simply repeating the same reasoning for vector fields. There is some analysis of embedding GS in continuous time in Hart2024 , wong2024contraction , and it seems plausible that similar techniques can be used to prove the result.
4 Isometric Embeddings
Having established that, generically, an embedding GS exists, we now move on to consider under what conditions an isometric embedding GS exists. An isometric embedding preserves both topological properties and geometrical properties, including lengths, angles, and curvatures.
Definition 4.1**.**
(Isometric embedding) Let be a Riemannian metric on the manifold . A map is an isometric embedding of if it is an embedding of and if
[TABLE]
The celebrated Nash-embedding Theorem states that for large enough any manifold can be isometrically embedded into .
Theorem 4.2**.**
(Nash nash1954 delellis2016masterpieces ) Let be a Riemannian manifold of dimension . Then for there exists an isometric embedding .
We will use Theorem 4.2 to prove that for sufficiently large reservoir dimension there exists a reservoir map that has an isometric embedding GS.
Theorem 4.3**.**
Suppose that has finitely many periodic orbits in . Suppose that for each periodic orbit with period , the derivative has distinct eigenvalues. Let and suppose satisfies the Nash bound of Theorem 4.2; . Then for generic observation function there exists a reservoir map whose generalized synchronization is an isometric embedding of .
Proof.
By Theorem 4.2, there exists an isometric embedding . By Theorem 2.5, there exists a reservoir map with embedding GS . Let be a diffeomorphism defined on a neighborhood of such that on . Define the conjugated reservoir map
[TABLE]
Then the reservoir system
[TABLE]
has the generalized synchronization which is an isometric embedding. ∎
Remark 4.4**.**
In the proof we constructed a reservoir map by conjugating with the (local) diffeomorphism . This is reminiscent of system isomorphism Tsuda1992 GrigoryevaOrtega2021 Grigoryeva2023 .
Definition 4.5** (System isomorphism Tsuda1992 GrigoryevaOrtega2021 Grigoryeva2023 ).**
Two reservoir maps are said to be isomorphic if there exists a diffeomorphism such that
[TABLE]
5 Linear Reservoir Systems
In the special case of the linear reservoir system where
[TABLE]
for and , with linear, and , we can construct reservoir systems with isometric embedding generalized synchronizations explicitly. This simultaneous linearization of both the evolution operator and reservoir system arises in the vicinity of hyperbolic fixed points Hartman1960 Grobman1959 . Furthermore, we will show that any linear reservoir system that has an embedding GS is isomorphic to one that is isometrically embedding.
Definition 5.1** (Linear system isomorphism).**
Let and be linear reservoir maps with and . We say that and are linearly isomorphic if there exists an invertible matrix such that
[TABLE]
Equivalently, and are isomorphic in the sense of Definition 4.5 with .
Theorem 5.2**.**
Let be a linear map, and represent the eigenvalue-eigenvector pairs of . Suppose the eigenvalues of are distinct and that . Suppose that and satisfy
[TABLE]
Furthermore, suppose that the vectors
[TABLE]
are linearly independent. Let be an orthonormal basis with respect to the Riemannian metric . We can express these vectors with respect to the eigenbasis as follows
[TABLE]
for a matrix . Let be a matrix with th column . Then for any rotation matrix the reservoir map
[TABLE]
where
[TABLE]
is linearly isomorphic to and has generalized synchronization
[TABLE]
which is an isometric embedding.
Proof.
To establish isometry we need to show that
[TABLE]
So we start with
[TABLE]
then observe that
[TABLE]
for canonical unit vector . Now introduce an arbitrary rotation matrix , which represents the arbitrary choice of orthonormal basis with respect to . Then
[TABLE]
Thus
[TABLE]
Now let
[TABLE]
for scalars and and observe that
[TABLE]
∎
Remark 5.3**.**
In the case , the choice of completing the columns of is not unique. One may choose any such that
[TABLE]
and then for any block rotation of the form described in the proof,
[TABLE]
Remark 5.4**.**
The full-rank condition on the matrix (ensuring ) is not restrictive in practice. Indeed, Proposition 4.4 in Grigoryeva2023 shows that when and are drawn randomly from any non-degenerate continuous distribution, the columns (for distinct ) are almost surely linearly independent. Hence, the condition we impose holds with probability 1 in typical random reservoir initializations.
Remark 5.5**.**
It is notable that the isometric form constructed above does not depend on the specific values of the eigenvalues beyond the necessary convergence conditions. In particular, isometry is achievable for any admissible set of eigenvalues, as long as the rank condition on is met. This highlights that the spectral data of do not obstruct the construction of an isometric embedding.
6 Conclusion
Generalized synchronization (GS) is the fundamental mechanism that allows reservoir systems to learn and represent dynamical systems. When the GS is an embedding, the reservoir state space provides a faithful reconstruction of the geometry of the underlying attractor, making possible the recovery of dynamical invariants and topological features. Our first main result shows that such embedding GS maps occur generically, formalizing what has been widely observed empirically in reservoir computing.
Going further, we established that isometric embeddings can also arise in reservoir systems. Isometric GS maps not only preserve topology, but also the metric structure of the source dynamics, including angles and lengths. In the linear case, we showed that every reservoir system with embedding GS is isomorphic to a system with isometric GS. The isometry preserves the eigenvalues of the reservoir matrix, proving that many different spectra are consistent with an isometric GS. This reveals a stronger structural property of reservoirs: they can, in principle, represent the geometry of dynamical systems without distortion. Future work could investigate continuous-time analogues of these results and explore how isometric GS maps might improve learning, stability, and invariant discovery in practical reservoir computing applications.
Appendix
Lemma 6.1**.**
Suppose that has finitely many periodic orbits in . Suppose that for each periodic orbit of with period , the derivative has distinct eigenvalues. Then for generic observation functions the set of generalised synchronisations is open and non-empty
Proof.
We let
[TABLE]
where
[TABLE]
so that
[TABLE]
This is the Takens delay map which is an embedding GS under the specified conditions on and for generic . This establishes non-emptiness. Now for any small perturbation of denoted define a reservoir map that satisfies
[TABLE]
and
[TABLE]
for in the normal space. Then we can smoothly continue on the remaining . Then . This establishes openness. ∎
Lemma 6.2**.**
If is open and nonempty and is dense and open (all in a topological space ), then is a dense open subset of (with the subspace topology).
Proof.
Openness. Since and are open in , their intersection is open in , hence open in with the subspace topology.
Density. Using the subspace–closure identity
[TABLE]
which is Proposition 16.4 in Munkres2000Topology , we have
[TABLE]
Because is dense in , , so . Hence is dense in . ∎
Lemma 6.3**.**
* is open and continuous.*
Proof.
First we define an equivalence relation on as if have generalised synchronisations that are equal on . Now the quotient map is open, continuous (in the quotient topology) and surjective. Now we define the map as mapping the equivalence class to the common GS . Now we observe that . Now using that is a homeomorphism by construction, is follows the composition is open and continuous which completes the proof. ∎
Lemma 6.4**.**
Let be topological spaces and let be continuous and open. If is comeager (generic), then is comeager (generic) in . In particular, if is a Baire space, then is dense in .
Proof.
Write with each open and dense in . By continuity, is open in for each . To see density, let be nonempty and open. Since is open, is open in , hence (as is dense). Pick ; then there exists with , so . Thus every nonempty open meets , i.e. is dense in . Therefore is comeager in . If is Baire, every comeager subset is dense, proving the last claim. ∎
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Allen G. Hart, James L. Hook, and Jonathan H. P. Dawes. Echo state networks trained by tikhonov least squares are l 2 ( μ ) l^{2}(\mu) approximators of ergodic dynamical systems. Physica D: Nonlinear Phenomena , 421:132882, 2021.
- 2[2] Allen G. Hart, James L. Hook, and Jonathan H. P. Dawes. Embedding and approximation theorems for echo state networks. Neural Networks , 128:234–247, 2020.
- 3[3] Allen G. Hart. Reservoir Computing with Dynamical Systems . Ph D thesis, University of Bath, 2022.
- 4[4] Allen G. Hart. Generalised synchronisations, embeddings, and approximations for continuous time reservoir computers. Physica D: Nonlinear Phenomena , 458:133956, 2024.
- 5[5] Lyudmila Grigoryeva, Allen G. Hart, and Juan-Pablo Ortega. Chaos on compact manifolds: Differentiable synchronizations beyond the Takens theorem. Physical Review E , 103(6):062204, 2021.
- 6[6] Lyudmila Grigoryeva, Allen G. Hart, and Juan-Pablo Ortega. Learning strange attractors with reservoir systems. Nonlinearity , 36(9):4674–4708, 2023.
- 7[7] Masashi Hara and Hiroshi Kokubu. Learning dynamics by reservoir computing. Journal of Dynamics and Differential Equations , 36(2):515–540, 2022.
- 8[8] Wancheng Hu, Yibin Zhang, Rencai Ma, and Qionglin Dai. Synchronization between two linearly coupled reservoir computers. Chaos, Solitons & Fractals , 157:111882, 2022.
