In the theory of discrete dynamical systems, the loss of stability of a fixed point can occur in several ways. While the Fold bifurcation deals with an eigenvalue crossing $+1$ and the Flip bifurcation with an eigenvalue crossing $-1$ (period-doubling), the **Neimark-Sacker bifurcation** describes the scenario where a pair of complex conjugate eigenvalues crosses the unit circle. This is the discrete-time analogue of the [Hopf bifurcation](/pages/1047) in [continuous](/page/Continuity) systems: the unit circle $|z| = 1$ plays the role of the imaginary axis, and the invariant closed curve that emerges plays the role of the [limit](/page/Limit) cycle.

The question the theory answers is: given a family of maps depending on a parameter $\alpha$, with a fixed point whose eigenvalues cross the unit circle as $\alpha$ crosses a critical value, does an invariant closed curve branch off from the fixed point? If so, is the curve attracting (nearby orbits spiral onto it) or repelling (nearby orbits are pushed away)? The answer depends on the interplay between the linear spectral crossing and the leading-order nonlinear terms, distilled into a single computable scalar — the First Lyapunov Coefficient — exactly as in the continuous case.

## Setup

We consider a smooth family of maps on $\mathbb{R}^2$ with a real parameter $\alpha \in \mathbb{R}$:

\begin{align*} x \mapsto f(x, \alpha), \quad x \in \mathbb{R}^2, \quad \alpha \in \mathbb{R}, \end{align*}

where $f \in C^k(U \times V, \mathbb{R}^2)$ with $k \ge 5$, $U \subseteq \mathbb{R}^2$ a neighbourhood of the origin, and $V \subseteq \mathbb{R}$ a parameter interval. We assume that $x = 0$ is a fixed point for all parameter values:

\begin{align*} f(0, \alpha) = 0 \quad \text{for every } \alpha \text{ near } 0. \end{align*}

This is not restrictive: if a family of fixed points $x^*(\alpha)$ depends smoothly on $\alpha$, the coordinate change $y = x - x^*(\alpha)$ reduces to this form. The linearisation at the origin is the $\alpha$-dependent Jacobian matrix $A(\alpha) := D_x f(0, \alpha)$, whose eigenvalues govern the local stability of the fixed point.

When the system is higher-dimensional ($x \in \mathbb{R}^n$ with $n > 2$), the [center manifold theorem](/page/Center%20Manifold%20Theorem) reduces the bifurcation analysis to a two-dimensional center manifold, just as in the continuous case. The essential dynamics are captured by the restriction of the map to this invariant surface.

## Spectral Hypothesis

The bifurcation requires the eigenvalues of $A(\alpha)$ to cross the unit circle through a conjugate pair.

[definition: Neimark-Sacker Spectral Hypothesis] The map $x \mapsto f(x, \alpha)$ satisfies the **Neimark-Sacker spectral hypothesis** at $\alpha = 0$ if the Jacobian $A(0) = D_x f(0, 0)$ has eigenvalues: \begin{align*} \lambda(0), \overline{\lambda(0)} = e^{\pm i\theta_0}, \quad \theta_0 \in (0, \pi), \end{align*} lying on the unit circle with non-zero imaginary part. [/definition]

The condition $\theta_0 \in (0, \pi)$ excludes real eigenvalues: $\theta_0 = 0$ gives a double eigenvalue at $+1$ (Fold), and $\theta_0 = \pi$ gives a double eigenvalue at $-1$ (Flip). The eigenvalues $e^{\pm i\theta_0}$ generate a discrete rotation of $\mathbb{R}^2$ by angle $\theta_0$ per iterate, which is the linear mechanism behind the closed invariant curve.

The key difference from the continuous Hopf setting is that the stability [boundary](/page/Boundary) is the unit circle $|\lambda| = 1$ rather than the imaginary axis $\mathrm{Re}(\lambda) = 0$. A fixed point of a map is stable when all eigenvalues satisfy $|\lambda| < 1$ (inside the unit disc) and unstable when any eigenvalue satisfies $|\lambda| > 1$.

## Complexification and Normal Form

To analyse the nonlinear dynamics near the bifurcation, we pass to complex coordinates on $\mathbb{R}^2$. This requires two steps: first, a real change of basis that puts the linear part in rotation form; second, a complexification that diagonalises the rotation.

### Putting the Linear Part in Rotation Form

The Jacobian $A(0)$ at the bifurcation has complex eigenvalues $\lambda, \bar{\lambda} = e^{\pm i\theta_0}$. A real $2 \times 2$ matrix with complex eigenvalues cannot be diagonalised over $\mathbb{R}$, but it can be brought to **rotation form** — the real [Jordan normal form](/theorems/864) for a conjugate pair.

The procedure uses the eigenvectors. Since $A(0)$ is real but $\lambda$ is complex, the eigenvector $w$ satisfying $A(0)w = \lambda w$ is necessarily complex. Write it as $w = p + iq$ where $p, q \in \mathbb{R}^2$ are the real and imaginary parts. Define the change-of-basis matrix:

\begin{align*} Q = \begin{pmatrix} | & | \\ p & q \\ | & | \end{pmatrix}. \end{align*}