Period-doubling renormalization operator ($\mathcal{T}(f)(x)=f(f(ax))/a$)

For a normalized even unimodal interval map with $f(0)=1$ and $a=f(1)\in (-1,0)$, defineThis returns the second iterate near the critical point to the original spatial and height normalization. Then $\mathcal{T}(f)(0)=f(1)/a=1$. A restrictive-interval condition is also needed for the composition to remain unimodal and map $[-1,1]$ into itself; normalization alone is not sufficient.

// Target: dynamical-systems.bigb

Fix $R>1$ and write $h(z)={\sum }_{n\ge 0}{h}_n((z-1)/R{)}^n$ with absolutely summable real coefficients. The coefficient space is a Banach space under $\Vert h\Vert ={\sum }_n|h_n|$. The normalized even maps $p(x)=1-x^{2}h(x^2)$ form an affine copy of that space and are analytic where $|x^2-1|<R$. On any smaller coefficient disk, the tail beyond degree $N$ is bounded by $q^{N+1}{\sum }_{n>N}|h_n|$ for a radius ratio $q<1$. Analytic composition and Cauchy estimates on smaller domains permit certified derivative and truncation bounds. A small ball about an admissible map must additionally preserve the quadratic critical maximum and restrictive-interval conditions; the coefficient norm alone does not impose them.

// Target: banach-space.bigb

A Feigenbaum fixed point is an even analytic unimodal interval map $g$ with $g(0)=1$, a critical maximum of even order $d$, $a=g(1)\in (-1,0)$ andIn the quadratic class $g^{\prime \prime }(0)<0$. Its spatial scaling uses the Feigenbaum constants convention $\alpha =-1/a$; the values depend on the critical order. The central restrictive intervals are successively reduced by $|a|$ under renormalization. This fixed point describes the limiting shape of repeatedly renormalized maps at a period-doubling accumulation.

// Target: dynamical-systems.bigb

An approximate analytic fixed point of the period-doubling renormalization operator can be certified by a frozen Newton correction for a fixed-point equation. A coefficient Banach space, rigorous residual and derivative bounds, and an a posteriori contraction ball give existence and local uniqueness in the full function space. Interval arithmetic controls finite computations; analytic estimates control the infinite coefficient tail.

// Target: dynamical-systems.bigb

Suppose $g(x)=1+c{x}^d+O(x^{d+2})$, with even $d$ and $c<0$. Expanding the period-doubling renormalization operator exactly at zero gives $a^{d-1}{g}^{\prime }(1)=1$. Approximating globally by $1+c{x}^d$ gives $c\approx a-1$ and $g^{\prime }(1)\approx dc$, hence $d(a-1)a^{d-1}\approx 1$. For $d=2$, $a\approx (1-\sqrt{3})/2$, $g(x)\approx 1-(1+\sqrt{3})x^2/2$ and $\alpha \approx 1+\sqrt{3}$. The coefficient relation is a polynomial-truncation approximation, not an exact identity for the full analytic fixed point.

// Target: numerical-analysis.bigb

At the normalized quadratic fixed point $g$, one expanding eigenvalue $\delta >1$ of $D\mathcal{T}(g)$ and a contracting complement give a stable manifold of finite codimension one. A generic one-parameter family crosses that manifold transversely at its cascade accumulation. Its unstable coordinate is initially proportional to parameter distance and multiplies by $\delta$ per renormalization, yielding $s_{\infty }-s_n\sim C{\delta }^{-n}$. Stable components disappear, while the fixed point's orientation-reversing scale $a$ yields signed spatial ratios $a^{-1}=-\alpha$.

// Target: dynamical-systems.bigb