In mathematics, and in particular the study of dynamical systems, the idea of stable and unstable sets or stable and unstable manifolds give a formal mathematical definition to the general notions embodied in the idea of an attractor or repellor. In the case of hyperbolic dynamics, the corresponding notion is that of the hyperbolic set.
The gravitational tidal forces acting on the rings of Saturn provide an easy-to-visualize physical example. The tidal forces flatten the ring into the equatorial plane, even as they stretch it out in the radial direction. Imagining the rings to be sand or gravel particles ("dust") in orbit around Saturn, the tidal forces are such that any perturbations that push particles above or below the equatorial plane results in that particle feeling a restoring force, pushing it back into the plane. Particles effectively oscillate in a harmonic well, damped by collisions. The stable direction is perpendicular to the ring. The unstable direction is along any radius, where forces stretch and pull particles apart. Two particles that start very near each other in phase space will experience radial forces causing them to diverge, radially. These forces have a positive Lyapunov exponent; the trajectories lie on a hyperbolic manifold, and the movement of particles is essentially chaotic, wandering through the rings. The center manifold is tangential to the rings, with particles experiencing neither compression nor stretching. This allows second-order gravitational forces to dominate, and so particles can be entrained by moons or moonlets in the rings, phase locking to them. The gravitational forces of the moons effectively provide a regularly repeating small kick, each time around the orbit, akin to a kicked rotor, such as found in a phase-locked loop.
The discrete-time motion of particles in the ring can be approximated by the Poincaré map. The map effectively provides the transfer matrix of the system. The eigenvector associated with the largest eigenvalue of the matrix is the Frobenius–Perron eigenvector, which is also the invariant measure, i.e the actual density of the particles in the ring. All other eigenvectors of the transfer matrix have smaller eigenvalues, and correspond to decaying modes.
The following provides a definition for the case of a system that is either an iterated function or has discrete-time dynamics. Similar notions apply for systems whose time evolution is given by a flow.
Let
X
f\colonX\toX
p
f
p
Ws(f,p)=\{q\inX:fn(q)\topasn\toinfty\}
p
Wu(f,p)=\{q\inX:f-n(q)\topasn\toinfty\}.
Here,
f-1
f
f\circf-1=f-1\circf=idX
idX
X
If
p
k
fk
p
Ws(f,p)=Ws(fk,p)
Wu(f,p)=Wu(fk,p).
U
p
p
s | |
W | |
loc |
(f,p,U)=\{q\inU:fn(q)\inUforeachn\geq0\}
u | |
W | |
loc |
(f,p,U)=
s | |
W | |
loc |
(f-1,p,U).
If
X
Ws(f,p)=\{q\inX:d(fn(q),fn(p))\to0forn\toinfty\}
Wu(f,p)=Ws(f-1,p),
where
d
X
p
Suppose now that
X
f
l{C}k
k\geq1
p
U
p
l{C}k
p
Es
Eu
Df(p)
f
l{C}k
Diffk(X)
l{C}k
X
l{C}k
If
X
f