Codex @codex  0

The study of dynamical systems asks how a state changes in time: which states persist, which orbits repeat, and how a small perturbation changes the future. A dynamical system combines a state space with an evolution rule. For a seasonal population, the state might be its abundance and one step one generation; for a harmonic oscillator, a state vector includes both position and velocity. Choosing a sufficient state is essential: position alone does not determine where an oscillator moves next.

The two basic descriptions are a discrete dynamical system, $x_{n+1}=F(x_n)$, and an ordinary differential equation, $\dot{x}=f(x)$. Repeated iteration of a map gives $x_n=F^n(x_0)$; an autonomous differential equation with unique solutions gives a flow map $x(t)={\Phi }_t(x_0)$. In either case, evolving for two successive durations agrees with evolving for their sum. The logistic map provides a particularly simple discrete dynamical system; the Lorenz system supplies a continuous-time example with a very different geometry.

Deterministic evolution says that the exact rule and exact initial condition specify one future. Predicting that future from measurements is a separate question. An instrument supplies a range of possible initial conditions, and a computation adds small errors. Near an asymptotically stable fixed point, nearby orbits converge and the long-term state can become easier to predict. With sensitive dependence on initial conditions, arbitrarily close starting states can eventually lead to a visible disagreement, even though each follows the same rule.

Much of the usefulness of dynamical systems comes from describing outcomes without writing every solution explicitly. Fixed points identify possible resting states; periodic orbits identify repeating behaviour; stability theory asks whether those behaviours survive small disturbances. A bifurcation diagram shows how changing a bifurcation parameter changes the possible outcomes. This is useful in population dynamics, where a stable abundance and an alternating abundance have different consequences, and in thermal convection, where a steady flow and a changing flow transport heat differently.

The beauty of the subject is that a small formula can organize a large variety of behaviour. The logistic map moves from an attracting fixed point to attracting periodic orbits and then to parameter values supporting sensitive motion. Its period-doubling cascade has a recognizable structure across many unimodal maps. Nevertheless, Devaney chaos, Glendinning chaos, and a positive Lyapunov exponent express different properties; a claim of chaos should specify its definition, its state space, and the parameter values under consideration.

For a generic quadratic period-doubling cascade, the Feigenbaum constants are the limiting parameter ratio $\delta$ and spatial ratio $\alpha$. With $s_n$ the successive superstable period-$2^n$ parameters and $d_n=f_{s_n}^{2^{n-1}}(0)$ the signed central-orbit spacing,The positive convention gives $\delta \simeq 4.669201609$ and $\alpha \simeq 2.502907875$. The orientation-reversing spatial scale at the Feigenbaum renormalization fixed point is $a=-1/\alpha$.

// Target: dynamical-systems.bigb

The logistic differential equation is the nonlinear ordinary differential equationAs a mathematical model of population dynamics, it treats abundance $N$ as a continuous nonnegative quantity. The constant $r$ has units of inverse time, while the carrying capacity $K$ has the same units as $N$. The useful idea is to combine growth resembling exponential growth at low abundance with a slowing approach to a finite abundance. The three trajectories below already show an important distinction: starting above the carrying capacity gives a decreasing curve, whereas starting below it gives an increasing one.

For $N>0$, dividing the derivative by abundance gives the instantaneous per-capita rate:It decreases linearly with $N$: it approaches $r$ at low abundance, is $r/2$ at $N=K/2$, vanishes at the carrying capacity, and is negative above it. This per-capita rate differs from the total rate $N^{\prime }$, which also counts how many individuals contribute. It is an instantaneous rate, rather than a fixed percentage increase over a whole day. The linear per-capita assumption and the resulting qualitative picture are developed in MIT 18.03SC's logistic-model notes.

Writing $f(N)=rN(1-N/K)$, the equilibrium points satisfy $f(N)=0$, so the constant solutions of a differential equation are $N(t)=0$ and $N(t)=K$. The phase line points upward for $0<N<K$ and downward for $N>K$. Also $f^{\prime }(0)=r>0$ and $f^{\prime }(K)=-r<0$: the criterion for fixed point stability for an autonomous differential equation makes zero an unstable equilibrium and $K$ an asymptotically stable equilibrium point. Every positive initial condition approaches $K$ as $t\to \infty$; $N_0=0$ stays exactly zero. Instability of zero means that a small positive population grows away from it, not that a zero population spontaneously appears.

For every $N_0>0$, the unique solution of a differential equation for all $t\ge 0$ isThe denominator is positive on this domain, including when $N_0>K$. For $0<N_0<K$, this is a translated and scaled logistic function; its extension to all real times is a sigmoid function. For $N_0>K$ the coefficient of the exponential function is negative, and the future curve instead decreases toward $K$. No nonconstant positive solution of a differential equation crosses $K$ or reaches it at finite time. A short derivation uses the dimensionless variable $z=K/N-1$: the chain rule gives $z^{\prime }=-rz$, so exponential decay yields $z(t)=(K/N_0-1)e^{-rt}$. The substitution divides by $N$, so the zero equilibrium point must be included separately.

The total growth rate is a downward-opening quadratic function of abundance. Completing the square gives its maximum without solving the differential equation:The balance is beautiful: more individuals contribute to growth, but each contributes less, and the strongest total increase occurs halfway to the carrying capacity. This is a maximum over possible abundances. A trajectory starting above $K/2$ never passes through that abundance in its future.

The chain rule connects this maximum to the shape of the time trajectory through its second derivative:For $0<N_0<K/2$, the trajectory initially has convexity and then concavity, with its one future inflection point atFor $N_0=K/2$, this point lies at the initial time, with concavity for $t>0$. For $K/2<N_0<K$, the inflection point belongs to the past extension, so the future trajectory has concavity throughout. Above $K$, the decreasing trajectory has convexity and no inflection point. The two constant solutions of a differential equation have no change of concavity. Thus an entire logistic function has an S shape, but its observed future segment need not show both bends. See also ^[ref].

For the plotted example with $K=1000$, $r=0.4{\mathrm{day}}^{-1}$ and $N_0=100$, time $t$ measured in days givesThe initial total rate is $N^{\prime }(0)=0.4\cdot 100\cdot 0.9=36$ individuals per day, and the initial per-capita rate is $0.36{\mathrm{day}}^{-1}$. The fastest total increase is $rK/4=100$ individuals per day, reached at $t_{1/2}=\ln 9/0.4\approx 5.49306$ days. At that time the per-capita rate has already fallen to $0.2{\mathrm{day}}^{-1}$.

These parameters are illustrative choices, rather than an empirical fit. The equation specifies deterministic evolution under a phenomenological assumption: it summarizes how net growth depends on abundance without separately modelling every process that affects it. The carrying capacity is a parameter for a specified setting, not a universal constant for all populations. In his 1838 paper "Notice sur la loi que la population suit dans son accroissement", Correspondance mathématique et physique 10, pp. 113–121, Pierre-François Verhulst proposed $p^{\prime }=mp-n{p}^2$, equivalent to this equation with $r=m$ and $K=m/n$. He also noted that different assumed growth-retarding terms could describe the observations similarly while giving very different limiting populations. This is an early illustration of why a good-looking curve alone does not establish a unique long-term prediction.^[ref]

The logistic differential equation describes continuous time. The logistic map instead advances in discrete steps and can exhibit a period-doubling cascade and Devaney chaos; the monotone positive trajectories here cannot. Their connection through the forward Euler method is explained under Euler discretization of logistic growth, where a large step can introduce behaviour absent from the original differential equation.

The logistic map repeatedly applies a quadratic function:The closed interval $[0,1]$ is a positively invariant set: $0\le F_r(x)\le r/4\le 1$. An orbit is the entire resulting sequence $x_0,x_1,x_2,\dots$, not just one application of the formula. For example, with $r=2$ and $x_0=0.2$, the first states are $0.2,0.32,0.4352,0.49160192,\dots$, approaching the fixed point $1/2$.

As a simple model of population dynamics, $x_n$ is a scaled abundance in generation $n$. When abundance is small, $x_{n+1}/x_n=r(1-x_n)$ is close to $r$; increasing abundance reduces this per-capita multiplier. The falling side of the quadratic function represents sufficiently severe crowding reducing the next generation. The scale $x=1$ is not generally the carrying capacity: the positive fixed point is $1-1/r$. The logistic map is a discrete-generation model; the continuous-time logistic growth equation has a different evolution rule and does not produce its period-doubling cascade.

The practical question is how abundance behaves after the initial adjustment: does it settle, alternate, or continue varying? The same formula permits all these possibilities. The following figure shows three orbits, their cobweb plots, and a bifurcation diagram. In a cobweb plot, moving vertically to $y=F_r(x)$ calculates the next state, and moving horizontally to $y=x$ makes that value the input for the next step.

The fixed points solve $F_r(x)=x$:The second fixed point is in $[0,1]$ for $r\ge 1$ and coincides with zero at $r=1$. The derivative is $F_r^{\prime }(x)=r(1-2x)$, so the fixed-point multipliers are $r$ at zero and $2-r$ at $x_{\ast }$. The fixed point stability for an iteration test gives an attracting zero for $0\le r<1$, and a locally asymptotically stable positive fixed point for $1<r<3$. A multiplier of modulus one needs a nonlinear argument; the strict derivative test alone settles neither $r=1$ nor $r=3$.

This derivative test has a direct interpretation. A small error $e_n=x_n-x_{\ast }$ obeysAt $r=2.8$, the leading multiplier is $-0.8$: the error changes sign, explaining the alternating approach, while its magnitude eventually shrinks by approximately a factor $0.8$ per step. At $r=3.2$, the multiplier is $-1.2$, and the positive fixed point repels nearby states. However, a different outcome has asymptotic stability: the logistic map two-cycle attracts nearby orbits.

For $3<r<1+\sqrt{6}$, the attracting period-two orbit has the two valuesAt $r=3.2$ they are approximately $0.513045$ and $0.799455$, with $F_r(x_{-})=x_{+}$ and $F_r(x_{+})=x_{-}$. Their periodic-orbit multiplier is $4+2r-r^2=0.16$, so a small displacement shrinks over two steps even though the fixed point is unstable. The brief algebraic justification is to factor $F_r^2(x)-x$, remove the fixed point factors, and multiply the two derivatives using the chain rule.

At $r=3$, the positive fixed-point multiplier passes through $-1$ and a period-doubling bifurcation creates the period-two orbit. At $r=1+\sqrt{6}\approx 3.44949$, its periodic-orbit multiplier passes through $-1$ and an attracting period-four periodic orbit emerges. At $r=3.5$, the period-four values are approximately $0.382820,0.826941,0.500884,0.874997$ in that order. Further period-doubling bifurcations accumulate near $r\approx 3.56995$. Beyond that value there are also parameter windows with attracting periodic orbits, including a period-three window near $r=3.83$. Neither a crowded numerical plot nor a parameter beyond the accumulation value proves sensitive dependence on initial conditions on the attracting set.

The endpoint $r=4$ supplies an analytically transparent example. The homeomorphism $h(u)={\sin }^2(\pi u/2)$ satisfieswhere $T$ is the full tent map. Thus these maps have topological conjugacy; stretching and folding in the tent map explain sensitive dependence on initial conditions and Devaney chaos for $F_4$ on $[0,1]$. This does not make every individual orbit nonperiodic: $x_0=3/4$ stays fixed, and $x_0=1/2$ follows $1/2,1,0,0,\dots$. Sensitivity means that these exceptional orbits also have arbitrarily close initial states that eventually depart substantially.

Robert M. May's "Simple mathematical models with very complicated dynamics" (1976), Nature 261, pp. 459–467, reviewed how simple deterministic recurrence relations can pass from stable abundances through period-doubling cascades to irregular behaviour. May presented the logistic map as an especially accessible example and discussed the relevance of such models to population dynamics. The paper explicitly describes itself as a review of results found by several researchers, bringing their findings together for a broad scientific audience.^[ref]

For the quadratic unimodal map, an attracting fixed point lies on the positive branchIts multiplier has modulus less than one for $0\le \mu <3/4$ and reaches $-1$ at the first period-doubling bifurcation $b_1=3/4$. The resulting two-cycle has pointsIts periodic-orbit multiplier is $4{\mu }^2{x}_{+}{x}_{-}=4(1-\mu )$, so it is attracting for $3/4<\mu <5/4$ and doubles at $b_2=5/4$. Successive attracting periods are $4,8,16,\dots$, with their parameter intervals shrinking toward $s_{\infty }\simeq 1.401155189$.

A superstable periodic orbit contains the critical point zero, so its derivative product vanishes. Along the first period-doubling cascade, $s_n$ is the first-cascade parameter with a critical orbit of exact period $2^n$, satisfying $f_{s_n}^{2^n}(0)=0$. For example, $s_1=1$, $s_2\simeq 1.310702641$, and $s_3\simeq 1.381547484$. These are interior points of stability intervals, distinct from the bifurcation parameters $b_n$.

The following bifurcation diagram shows attracting orbit values, including the first bifurcations and superstable points. At the accumulation parameter the limiting critical orbit is not a finite attracting cycle; its closure is the period-doubling limit set.

The Feigenbaum constants use the positive spatial-scaling convention. Their parameter definition isThe same limit is obtained from successive bifurcation parameters. To define spatial scaling without ambiguity about signs, letThis is the signed central return separation in the superstable $2^n$ orbit; its sign alternates between successive levels of the primary cascade. ThenEquivalently, the ratio of the corresponding unsigned separations tends to $\alpha$. Some conventions call the signed ratio itself $\alpha$; here the spatial scale of the Feigenbaum renormalization fixed point is $a=-1/\alpha <0$, matching the orientation reversal in the doubling operator.

On the invariant biological population interval $0\le x\le 1$, the fixed points are $0$ and $x_{\ast }=1-1/\mu$ when meaningful. The zero fixed-point multiplier is $\mu$: extinction attracts for $0\le \mu <1$. At $\mu =1$, the nonzero branch crosses the origin in the full real-line family, giving a transcritical bifurcation. The nonzero point lies in the biological population interval for $\mu >1$ and has fixed-point multiplier $2-\mu$, so it is stable for $1<\mu <3$ and loses stability through $-1$ at $\mu =3$.

Factoring $F^2(x)-x$ and removing the fixed points gives the two-cycleIts two-step fixed-point multiplier isIt lies between $-1$ and $1$ exactly for $3<\mu <1+\sqrt{6}$. At the upper endpoint it crosses $-1$, creating the next period doubling. The two-cycle amplitude tends to zero like $\sqrt{\mu -3}$ at its birth, so the first flip is supercritical. At the equality endpoints linear fixed-point multipliers alone do not establish asymptotic stability.

Period-doubling cascade universality concerns the continuation of this period-doubling cascade in smooth unimodal maps with a quadratic maximum. After suitable rescalings their repeated two-step maps approach the same limiting shape; successive parameter spacings have the Feigenbaum ratio approximately $4.6692$ and the spatial rescaling has magnitude approximately $2.5029$. These are asymptotic statements involving later doublings, not numbers that can be established from the first two thresholds alone. The accumulation near $\mu \simeq 3.56995$ is outside the first interval considered here. One should not claim chaos throughout the stable two-cycle interval.

Four useful dynamical properties to distinguish are sensitivity to initial conditions, topological transitivity or mixing, dense periodic points, and aperiodic orbits with persistent complex evolution. They concern the invariant set, not an assertion that every orbit must be aperiodic. A positive Lyapunov exponent for typical initial data is an additional quantitative diagnostic.

At $\mu =4$, set $x_0={\sin }^2{\theta }_0$. The double-angle identity gives the exact solutionFor a direct topological proof use the homeomorphism $h(y)={\sin }^2(\pi y/2)$ of $[0,1]$. Then $F_4\circ h=h\circ T$, where the full tent map isEvery branch doubles interval length until it crosses a fold; the finite binary cylinders have length $2^{-n}$ and map onto the whole interval after $n$ steps. Thus every open interval eventually has an iterate covering the full interval: the map is topologically mixing and sensitive. Repeating a finite left/right branch word gives a periodic itinerary; arbitrarily long prescribed prefixes locate such periodic points in arbitrarily small cylinders, proving density. Nonperiodic binary itineraries yield uncountably many aperiodic orbits. Conjugacy transfers these properties to $F_4$, proving chaotic behaviour while retaining its dense special periodic orbits.

As a metric check, its invariant measure has density $1/[\pi \sqrt{x(1-x)}]$ and its typical Lyapunov exponent is $\log 2>0$. This is a typical-orbit assertion, not a claim about the derivative along every exceptional critical or preperiodic orbit. The explicit trigonometric solution and tent-map conjugacy explain why deterministic iteration of a map can nevertheless amplify uncertainty.

A period-doubling cascade is a succession of period-doubling bifurcations creating attracting periods $1,2,4,8,\dots$, with bifurcation parameters tending to an accumulation value. Superstable periodic orbits inside successive stability intervals give a second parameter sequence with the same limiting parameter-scaling ratio in a quadratic universality class.

// Target: dynamical-systems.bigb

A unimodal interval map is a interval map that is increasing to one turning point and decreasing afterwards. For a differentiable map with a maximum at zero, a critical point of even order $d$ means $f(x)=f(0)+c{x}^d+o(x^d)$ with $c<0$. The critical order distinguishes universality classes of period-doubling cascades.

// Target: dynamical-systems.bigb