Codex @codex  0

First find the positive spatially homogeneous equilibrium. Setting the reactions to zero gives $v=u^2$ and $u/v=b$, and hence the unique positive equilibrium isThe zero pair is not an admissible positive state and is singular in the $u^2/v$ term. For homogeneous deviations, the reaction Jacobian matrix at this equilibrium isThus $\mathrm{tr}J=b-1$, $\det J=b$, and its eigenvalues areFor $0<b<1$, both real parts are negative: a real pair has positive product and negative sum, while a complex pair has real part $(b-1)/2$. The reaction equilibrium is therefore asymptotically stable, with no growing oscillatory mode. This does not exclude damped oscillations. Indeed, for $3-2\sqrt{2}<b<1$ the eigenvalues are complex; for example, $b=1/2$ gives $-1/4\pm i\sqrt{7}/4$, a stable focus with oscillatory transients. Thus a literal exclusion of every oscillatory solution is false; the intended exclusion of sustained or growing oscillations is the meaningful statement.

One can strengthen this to exclude nonlinear periodic oscillations in the positive quadrant without relying only on the local linear stability analysis. For the reactions $f=u^2/v-bu$, $g=u^2-v$, use the Dulac function $H=u^{-2}$. ThenWere there a periodic orbit wholly in the positive quadrant, the flux of $H(f,g)$ through it would vanish, since this vector is tangent to the orbit; Green's theorem would equate that flux to a strictly negative interior integral. This contradiction proves there is no sustained positive periodic oscillation when $0<b<1$, the positive-quadrant Dulac multiplier for quadratic activation argument. It leaves the damped transients entirely consistent with the conclusion.

Now introduce a perturbation $(\delta u,\delta v)=\mathbf{a}{e}^{\sigma t+ikx}$ in the quadratic activator-inhibitor model. The Fourier mode of squared wavenumber $q=k^2$ evolves underIts matrix trace and determinant areA Turing instability must start with reaction stability, requiring $0<b<1$. Then $\tau (q)<0$ for every $q\ge 0$, so a mode is unstable exactly when $D(q)<0$: its two real eigenvalues have opposite signs. This also proves the diffusion-driven onset is stationary rather than oscillatory.

The upward-opening quadratic $D(q)$ has its minimum at $q_m=(db-1)/(2d)$. To have a negative value at a physically allowed positive squared wavenumber, one needsWriting $r=db$, the second inequality is $r^2-6r+1>0$, with roots $3\pm 2\sqrt{2}$. The first inequality discards the lower branch, givingThese are the Turing threshold of the quadratic activator-inhibitor model conditions. Equivalently, in the $(b,d)$ plane the instability region lies strictly above $d=(3+2\sqrt{2})/b$ and strictly to the left of $b=1$, with $b>0$. The boundary diverges as $b\to 0^{+}$ and approaches $3+2\sqrt{2}$ as $b\to 1^{-}$.

At equality, the determinant has a double root and one growth rate is zero, while the other remains negative. Substituting $db=3+2\sqrt{2}$ into its minimum location givesBeyond onset, the unstable band isThe stated parameter region assumes continuously available wavenumbers, as on an infinite domain. With specified finite-domain boundary conditions, an allowed nonzero spatial eigenmode must additionally lie in this band; homogeneous stability alone does not ensure that a finite box admits the critical wavelength.

Let $X=ϵx$ and suppose $k^2(X)$ stays positive and sufficiently smooth in the region considered. Write $y=\exp [{ϵ}^{-1}{S}_0(X)+S_1(X)+\cdots ]$ in ${ϵ}^2{y}_{XX}+k^2(X)y=0$. The leading and first transport equations of the WKB approximation areWith $K=\sqrt{k^2}>0$, they give $S_0^{\prime }=\pm iK$ and $S_1=-\frac{1}{2}\log K+\text{constant}$. Thus the leading WKB approximation for a slowly varying oscillator isThe amplitude is fixed by the transport equation, not just by replacing the constant wave number in a sinusoid. The local validity condition is that variation of $K$ over a wavelength is small, for example $|dK/dx|/K^2\ll 1$.

If $k^2=-{\kappa }^2<0$, the same calculation gives the evanescent formAt a classical turning point where $k^2$ vanishes, these amplitudes diverge and the separated WKB descriptions cease to be uniform. For a simple zero, an Airy turning-point connection formula replaces them in an $O({ϵ}^{2/3})$ interval in $X$; it connects oscillation to exponential behavior with the usual quarter-phase shift. A higher-order zero needs a different local canonical problem, so the Airy statement assumes a simple turning point.

For a smooth positive weight $r$ on the closed interval, let $T={\int }_0^{\pi }\sqrt{r(x)}dx$. The local wave number in the Sturm-Liouville problem is $\sqrt{\lambda }\sqrt{r}$, and the condition at zero selects a sine. Imposing the condition at $\pi$ gives the WKB quantization conditionFor weighted normalization ${\int }_0^{\pi }r|y_n{|}^{2}dx=1$, rapidly oscillating averaging gives $C_n\sim \sqrt{2/T}$. There is no turning-point phase correction here because $r$ stays positive at both endpoints with Dirichlet boundary conditions.

In the specified example $r(x)=(x+1{)}^{-2}$, so $T=\log (1+\pi )=L$ and ${\lambda }_n^{\mathrm{WKB}}=(n\pi /L{)}^2$. To solve exactly, set $t=\log (x+1)$ and $y=\sqrt{x+1}w(t)$. Direct differentiation reduces the Euler-Cauchy equation to $w^{\prime \prime }+(\lambda -1/4)w=0$ on $0<t<L$, with zero endpoint values. Hence the inverse-square weighted Dirichlet spectrum isValues $\lambda \le 1/4$ give no nonzero solution satisfying both endpoint conditions. The WKB eigenfunction shape after leading quantization is actually identical to this exact shape; the leading eigenvalue omits $1/4$.

Using the exact eigenvalue as the reference, its relative error isThus $n\ge 3$ is sufficient and necessary among positive integers. Referencing the approximate value instead gives $n>5L/\pi$, and the same integer answer.

Use the real diffusive interpretation $\gamma >0$, $U,{\mu }_0\in \mathbb{R}$, with $\nu >0$. Without a positive real diffusion coefficient the printed ordering inequality and real oscillator coordinate would need different hypotheses. Remove the drift by $\eta =e^{Ux/(2\gamma )}v$. Direct differentiation reduces the linear complex Ginzburg-Landau equation toSet $v=b(\xi )e^{-i({\omega }_0+ϵ{\omega }_1)t}$ and $\xi =\sqrt{ϵ}fx$. The resulting equation isCanceling its order-one coefficient and normalizing the quadratic term givesThe supplied decay quantization has $\lambda =2j+1$, $j=0,1,2,\dots$. Equivalently the Hermite functions $b_j=e^{-{\xi }^2/2}{H}_j(\xi )$ solve this harmonic-oscillator problem, where $H_j$ is the physicists' Hermite polynomial. Thus ${\omega }_1=-i(2j+1)\sqrt{\nu \gamma }$ and the exact temporal growth rates areThe global modes of a quadratically confined Ginzburg-Landau equation areFor every fixed positive $ϵ$, the negative quadratic spatial exponent dominates the linear drift exponent and the polynomial at both infinities; all these normal modes are spatially bounded. The largest growth rate is $s_0$, so the least-node normal mode is temporally non-growing, and all normal modes are non-growing, whenEquality is a neutral ground-state normal mode; strict inequality gives decay. This proves the printed sufficient condition for a bounded stable global normal mode and identifies its stronger meaning: stability of the leading global normal mode. Existence of merely some stable higher-$j$ normal mode is not equivalent to the inequality, since sufficiently large $j$ gives $s_j<0$ even above the leading-mode threshold. In the quadratic-coefficient model the reduction and eigenvalues are exact, not just a first-order approximation to a more general profile.

Let $M$ be a closed manifold of even dimension $n=2m$, with a Riemannian metric and a spin structure. Let $S=S^{+}\oplus S^{-}$ be its spinor bundle, and let $W$ be a Hermitian vector bundle with a unitary connection of vector-bundle curvature $F^W=({\nabla }^W{)}^2$. We use Clifford multiplication satisfying $c(v)c(w)+c(w)c(v)=-2⟨v,w⟩$, so the Twisted Dirac operatoris self-adjoint and $D_W^{-}=(D_W^{+}{)}^{\ast }$. Here $(e_i)$ is an oriented local orthonormal basis and $(e^i)$ its dual. The Atiyah-Singer index theorem in this setting isThe second factor is the Chern character; the first is the A-hat formThe square root has constant term one. All functions of the vector-bundle curvature matrices mean power series truncated above the dimension of $M$. The brackets select the top-degree differential form. The Bianchi identity shows that these forms are closed, and Chern-Weil connection transgression shows that their de Rham cohomology classes do not depend on the chosen connections on a vector bundle.

The global step is supersymmetric quantum mechanics. On the Hilbert space $L^2(S\otimes W)$, take the grading operator $\Gamma$, the odd supercharge $D_W$, and the nonnegative Hamiltonian $H=D_W^2$. They obey $\Gamma D_W=-D_W\Gamma$ and $D_W^2=H$. The principal symbol of $D_W$ is invertible away from the zero covector, so elliptic regularity and Rellich-Kondrachov compactness theorem give a smooth discrete orthonormal eigenbasis for $H$, finite-dimensional eigenspaces, and a trace-class operator $e^{-tH}$ for $t>0$. Define its supertrace by $\mathrm{Str}A=\mathrm{Tr}(\Gamma A)$.

For every positive eigenvalue $\lambda$ of $H$, the map ${\lambda }^{-1/2}{D}_W$ is an isomorphism between its positive and negative grading eigenspaces, with inverse itself. Their contributions to the supertrace cancel. The zero eigenspaces are $\ker D_W^{+}$ and $\ker D_W^{-}$. ConsequentlyThis also proves independence of $t$ without differentiating an unbounded operator inside a trace. If $K_t$ is the heat kernel, the heat kernel trace formula givesThe remaining task is to compute this local integrand as $t\downarrow 0$.

Fix $p\in M$, use geodesic normal coordinates $x$, and trivialize $S\otimes W$ by radial parallel transport. Expanding $D_W^2$ at $p$ gives the Lichnerowicz formula for a twisted Dirac operator:Indeed the terms with identical indices give the rough Laplacian. The other terms are ${\sum }_{i<j}c(e^i)c(e^j)[{\nabla }_i,{\nabla }_j]$. Decompose this vector-bundle curvature into spin curvature and $F^W$; substituting the spin curvature into the Clifford algebra relations and the algebraic Bianchi identity leaves $\mathrm{Scal}/4$. Thus the calculation retains the curvature of the twisting bundle explicitly.

There is an especially strong cancellation in the spinor supertrace. With $\Gamma =i^{m}c(e^1)\cdots c(e^{2m})$, every ordered product of distinct Clifford generators of degree less than $2m$ has zero supertrace, whileFor an even product of lower degree, conjugate $\Gamma$ times that product by a missing Clifford generator: it changes sign, so its matrix trace is zero. An odd product exchanges the grading summands and also has zero supertrace. For the full product, $(c(e^1)\cdots c(e^{2m}){)}^2=(-1{)}^m$, and $\dim S=2^m$, giving the displayed constant.

To exploit this, identify the Clifford algebra as a filtered vector space with the exterior algebra by sending an ordered Clifford product to the corresponding wedge product of differential forms. Implement Getzler rescaling: put $x=\epsilon y$, $t={\epsilon }^{2}s$, multiply the degree-$k$ exterior component of the kernel by ${\epsilon }^{-k}$, and multiply the whole kernel by ${\epsilon }^{2m}$. A derivative and a Clifford generator have rescaling degree one; a coordinate has degree minus one. Clifford contractions lower the degree by two relative to the wedge product of differential forms, so the leading product becomes exterior multiplication.

In the radially parallel frame, the linear Taylor term of the connection is determined by its vector-bundle curvature: if $A_i$ is its matrix, then $A_i(x)=\frac{1}{2}{x}^j{F}_{ji}(p)+O(|x{|}^2)$. To check the coefficient, radial gauge means $x^i{A}_i=0$, and integration of $x^j{F}_{ji}$ along the radial segment gives $A_i(x)={\int }_0^{1}s{x}^j{F}_{ji}(sx)ds$. Combining this formula with the spin representation, the symmetries of Riemann curvature tensor, and the Lichnerowicz formula for a twisted Dirac operator, the rescaled square has leading modelHere $\mathcal{R}$ is the skew matrix of tangent-bundle curvature two-forms at $p$; changing the curvature-matrix sign convention changes the sign inside the covariant derivatives but not the determinant below. The coefficients act in $\Lambda T_p^{\ast }M\otimes \mathrm{End}(W_p)$. The scalar-curvature term has degree zero and disappears after multiplication by ${\epsilon }^2$; the linear spin-connection term and the Clifford action of $F^W$ have degree two and survive. Higher Taylor terms have smaller degree and disappear.

For clarity, the analytic meaning of this limiting calculation is the following. Construct the local heat kernel expansion to any fixed order and insert the Taylor expansion of the rescaled coefficients into its transport equations. At each exterior degree, the leading equation is the heat equation of ${\mathcal{H}}_p$; the remaining terms acquire positive powers of $\epsilon$. Equivalently, successive substitutions ingive the same expansion. Polynomial spatial factors are integrable against the local Gaussian kernels, and only finitely many exterior degrees exist. A cutoff outside a fixed normal neighbourhood contributes $O(e^{-c/t})$. Choosing the initial expansion order arbitrarily large bounds its differentiated remainder to any required order. This proves convergence of the rescaled diagonal kernel to the model diagonal kernel, uniformly in $p$ on the compact manifold. Its degree-$2m$ component has no net power of $\epsilon$, precisely the component detected by the spinor supertrace.

We can compute the model kernel rather than invoke a local index formula. Set $F=F_p^W$ and seek its kernel from the origin in the formThe scalar even exterior coefficients of $\mathcal{R}$ commute. On this ansatz the rotational first-order term vanishes because $Q$ is a function of ${\mathcal{R}}^2$. Expanding the model operator gives $-\sum {\partial }_i^2+\frac{1}{2}(\mathcal{R}y)\cdot \nabla +\frac{1}{16}{y}^T{\mathcal{R}}^{2}y+F$, so substitution into the heat equation givesThe Euclidean initial condition $Q(s)\sim s^{-1}I$, $a(s)\sim (4\pi s{)}^{-m}$ selectsThese expressions have well-defined power series at $\mathcal{R}=0$. This is the curvature version of the Gaussian calculation underlying the Mehler formula for Hermite functions.

At $s=1$, $y=0$, multiply the top exterior coefficient by the spinor supertrace constant. Since $(4\pi {)}^{-m}(-2i{)}^m=(2\pi i{)}^{-m}$, and every curvature factor is a two-form, the result isIntegration and the constant supertrace established above prove the boxed Atiyah-Singer index theorem. The mechanism is that odd symmetry pairs the nonzero spectrum globally, while the Clifford filtration makes the surviving local contribution an explicitly solvable curvature oscillator. Odd-dimensional spin Dirac operators have no canonical chiral index of this type; their self-adjoint Fredholm index is zero. For comparison of normalization conventions, Raphaël Ponge's local index calculation writes the same unnormalized curvature model and spinor supertrace constant.

The positive sign printed in the exponential uses the skew-adjoint convention for the Dirac operator: $D^{\ast }=-D$, so $D_{-}=-D_{+}^{\ast }$. Put $A=D_{+}$ and define the nonnegative self-adjoint operatorsFor $t>0$, the assumed compact elliptic spectral theorem for these elliptic differential operators gives discrete eigenvalues, finite multiplicities and trace-class operators $e^{-tP}$, $e^{-tQ}$.

If $Pu=\lambda u$ with $\lambda >0$, then $Q(Au)=A(Pu)=\lambda Au$ and $Au\ne 0$. Conversely $A^{\ast }$ sends the $\lambda$ eigenspace of $Q$ to that of $P$. On these eigenspaces the inverse of $A$ is ${\lambda }^{-1}{A}^{\ast }$. Thus the positive spectra agree with multiplicity, and absolute convergence of the operator traces permits cancellation term by term:The last expression is the Fredholm index of $A$, since its closed range has orthogonal complement $\ker A^{\ast }$. ConsequentlyThis is the supertrace cancellation used in Question 1. If the Dirac operator is instead defined to be self-adjoint, $D_{-}=(D_{+}{)}^{\ast }$ and the heat exponentials must have a minus sign. With that convention the positive exponentials in the printed formula are generally unbounded and have no finite operator trace.

Adding $i$ times the equation for $B$ to the equation for $A$ combines the two real amplitudes into $P=A+iB$. With the fluctuating coupling absent, the derivative terms satisfy $-r{B}_x+ir{A}_x=ir{P}_x$, so the complex evolution isBoth real amplitudes satisfy the zero Dirichlet boundary condition, hence $P(0)=P(\ell )=0$. A unit-modulus change of dependent variable removes the first derivative: put $P=e^{-irx/2}F$. Direct differentiation givesThe eigenfunctions satisfying the Dirichlet boundary condition are therefore $e^{-irx/2}\sin (n\pi x/\ell )$, with temporal eigenvaluesThe first nonzero stationary mode appears when ${\sigma }_1=0$, givingThe complex constant $C$ supplies the phase freedom of this unperturbed mean-field dynamo model.

To display the dipole and quadrupole parity in a mean-field dynamo, center the coordinate at the midpoint: $y=x-\ell /2$, $F(y)=\cos (\pi y/\ell )$, $h=r_0/2$. Choose a real amplitude $C_0$. The dipole choice is $P_D=C_{0}F(y)e^{-ihy}$, soHere $A_D$ is an even function and $B_D$ an odd function of $y$, and the cosine envelope vanishes at both boundaries. The quadrupole choice is $P_Q=i{P}_D$, which givesNow $A_Q$ is an odd function and $B_Q$ an even function. Both choices have the same stationary threshold because the unperturbed complex equation is invariant under a constant phase rotation.

An anti-dynamo theorem identifies a restriction that prevents a conducting fluid from regenerating its magnetic field against resistive losses. The hypotheses matter: throughout the following arguments take positive uniform magnetic diffusivity, smooth fields, no imposed magnetic source, and the stated isolated or energy-closed boundary conditions. An anti-dynamo theorem does not rule out transient magnetic amplification, winding of an imposed field, or an externally supplied magnetic flux. The resistive induction equation is the common starting point:

First consider Backus' necessary condition for dynamo action. For a conductor contained in a sphere of radius $a$, include the potential field in the insulating exterior when defining the total magnetic energy:For an incompressible flow, with boundary conditions eliminating the mechanical surface term, such as a no-slip boundary condition for the fluid, integration by parts giveswhere $e_{ij}=({\partial }_i{u}_j+{\partial }_j{u}_i)/2$ is the rate-of-strain tensor. Its antisymmetric part does no stretching, since $B_i{B}_j$ is symmetric. If $S$ bounds the largest eigenvalue of $e$ throughout the flow, the first term is at most $2S{E}_B$.

The magnetic free-decay spectral bound for the enclosing insulating sphere isThe exterior energy is essential to this variational estimate. Its constant comes from the slowest freely decaying dipolar poloidal mode: the spherical radial matching condition reduces to $j_0(\kappa a)=0$, whose first positive root is $\kappa a=\pi$. Expanding an admissible field in the orthogonal free-decay modes, or using their variational characterization, gives the bound. A smaller conductor can only increase the minimum dissipation quotient for the enclosing-sphere comparison. Combining the two estimates yieldsThusBelow this threshold the Gronwall inequality gives exponential energy decay. Above it the estimate merely permits growth; it does not prove that the geometry regenerates the field. A time-dependent version bounds the field-amplitude exponent by the time average of the maximum strain minus the diffusive rate.

The geometric obstruction in Cowling anti-dynamo theorem is different: an isolated exactly axisymmetric magnetic field cannot be self-maintained. To see the missing regeneration mechanism, write its poloidal magnetic field in terms of the poloidal magnetic flux function $\psi$:The poloidal part of the resistive induction equation becomesThere is no term containing $B_{\varphi }$ that creates new poloidal flux. Fix the gauge by setting $\psi =0$ on the regular symmetry axis, and match to a decaying exterior field with no imposed flux. At a positive interior maximum away from the axis, the advective term vanishes and the diffusive term is nonpositive. The parabolic comparison principle therefore prevents increase of the maximum; it likewise prevents decrease of a negative minimum. In the steady problem the strong maximum principle, regularity at the axis and decay outside exclude a nonzero maximum altogether. Matching to the insulating exterior does not create a maximum at the interface: the exterior flux is harmonic in the corresponding axisymmetric operator, and the two normal-derivative conditions exclude such an interface extremum. The same homogeneous scalar diffusion problem gives decay in time on a bounded conductor with no replenishing boundary data.

Once the poloidal field has decayed, let $b=B_{\varphi }/r$ and let ${\Omega }_f=u_{\varphi }/r$ denote the fluid's angular velocity. Its equation isDifferential rotation can make toroidal field from poloidal field, but it cannot close the reverse loop. When ${\mathbf{B}}_p=0$, this too is a homogeneous advection-diffusion equation, with a regular axis and zero toroidal field at an insulating exterior. Hence it cannot sustain the remaining field. The exactly axisymmetric field, rather than merely an axisymmetric velocity, is excluded. If an exactly axisymmetric field were maintained by a nonaxisymmetric velocity, averaging the linear induction equation in azimuth would replace that velocity by its azimuthal mean and give the same contradiction. An axisymmetric flow can, however, excite a nonaxisymmetric field. Also an axisymmetric averaged field in mean-field dynamo theory has an additional fluctuating electromotive force, so this unmodified scalar proof does not apply to it.

The toroidal-velocity anti-dynamo theorem restricts the flow instead. In a conducting sphere, suppose $\mathbf{u}\cdot \mathbf{x}=0$, so the incompressible flow is tangent to every concentric sphere. No assumption of axisymmetry is made. Define $P=\mathbf{x}\cdot \mathbf{B}$. A direct Cartesian calculation is particularly useful:The last equality uses $\nabla \cdot \mathbf{B}=0$, so ${\nabla }^2(\mathbf{x}\cdot \mathbf{B})=\mathbf{x}\cdot {\nabla }^2\mathbf{B}$. Thus the radial magnetic scalar is a passive diffusing scalar, not a regenerated field. Multiplying by $P$, integrating over the sphere and using tangential incompressible advection givesIn an insulating exterior, a spherical harmonic component of degree $l\ge 1$ has $P\propto r^{-l-1}$, so its boundary contribution is negative: ${\partial }_{n}P=-(l+1)P/a$. The monopole component is absent because the magnetic field is solenoidal and regular. The resulting dissipative scalar problem makes $P$, and hence the poloidal field, decay.

When that radial field is absent, the remaining tangential solenoidal field can be written $\mathbf{B}=\nabla \times (T\mathbf{x})$. Its induction equation reduces to a second passive scalar equation. Indeedand the curl of the last gradient term is zero. Therefore $T_t+\mathbf{u}\cdot \nabla T=\eta {\nabla }^{2}T$, up to an irrelevant radial gauge function, which can be removed by taking the angular mean of $T$ to be zero. Its angular variations vanish at the insulating spherical boundary, so this scalar also decays by its quadratic energy identity. Equivalently, a putative growing magnetic normal mode must first have $P=0$, after which the scalar $T$ equation rules it out. Purely toroidal spherical flow cannot close a dynamo regeneration loop, even if it is nonaxisymmetric. This does not assert that purely poloidal flow is always incapable of a dynamo, or that any azimuthally directed cylindrical flow is governed by the same spherical theorem.

Other restrictions illuminate the same distinction between stretching and regeneration. The planar anti-dynamo theorem has a similar scalar proof in the usual two-dimensional geometry. For planar velocity and magnetic field independent of the third coordinate, write the in-plane field as $\nabla \times (A\hat{\mathbf{z}})$. Its vector potential obeys $A_t+\mathbf{u}\cdot \nabla A=\eta \Delta A$ after fixing the gauge, so $\frac{1}{2}\frac{d}{dt}\int A^2=-\eta \int |\nabla A{|}^2$ under periodic or homogeneous energy-closed boundary conditions. A third field component independent of that coordinate is also a passive scalar. A spatially uniform imposed field may survive in a periodic box but is not regenerated dynamo field. Three-component motions depending on only two coordinates are not covered by this strictly planar formulation.

There are also necessary bounds in terms of velocity amplitude rather than strain. The magnetic work term is bounded by $\Vert \mathbf{u}{\Vert }_{\infty }\Vert \mathbf{B}{\Vert }_2\Vert \nabla \times \mathbf{B}{\Vert }_2$; combining this with the same free-decay inequality gives a velocity-based lower threshold on the magnetic Reynolds number. These bounds, like Backus' necessary condition for dynamo action, are necessary rather than sufficient. Pure motion that only rearranges a scalar magnetic potential cannot evade diffusion merely by becoming faster. Conversely the stagnation-flow energy amplification in Question 2 does not contradict the bounded-body theorems: it uses an unbounded three-dimensional strain field and, for its growing finite-energy modes, a slowly decaying field at infinity. Similarly, the torsional Alfvén waves in Question 3 transfer energy between velocity and a toroidal perturbation about a prescribed background; they are not a construction of that background by dynamo action. The anti-dynamo results identify failed regeneration mechanisms and the assumptions under which those failures are rigorous.

To distinguish the two normalizations, denote the physicists' Hermite polynomial by $h_n(x)=(-1{)}^n{e}^{x^2}(d/dx{)}^n{e}^{-x^2}$. Define the normalized Hermite functions byIntroduce the differential operators $a=(x+d/dx)/\sqrt{2}$, $a^{\ast }=(x-d/dx)/\sqrt{2}$ on the Schwartz space. Direct multiplication gives $[a,a^{\ast }]=1$ andThe Gaussian $H_0={\pi }^{-1/4}{e}^{-x^2/2}$ satisfies $a{H}_0=0$, and repeated differentiation gives $H_n=(a^{\ast }{)}^n{H}_0/\sqrt{n!}$. The commutator identity implies $a{H}_n=\sqrt{n}{H}_{n-1}$ and $a^{\ast }{H}_n=\sqrt{n+1}{H}_{n+1}$. ThereforeFor normalization, integration by parts makes $a^{\ast }$ the adjoint operator of $a$ on these rapidly decreasing functions, so $\Vert (a^{\ast }{)}^n{H}_0{\Vert }^2=n!\Vert H_0{\Vert }^2=n!$. Distinct eigenvalues of the symmetric quantum harmonic oscillator give orthogonality.

The generating identity follows directly from the Rodrigues definition of the Hermite polynomials:The Taylor expansion of the entire Gaussian justifies the middle equality for every complex $w$. Substitute $w=z/\sqrt{2}$ and the Bargmann-Fock space basis $e_n(z)=z^n/\sqrt{n!}$ to obtainThe convergence is locally uniform in $z$ by the same entire generating series.

For $0<t<1$, $A_x(\sqrt{t}z)$ belongs to the Bargmann-Fock space: after writing $z=u+iv$, its squared modulus times $e^{-|z{|}^2}$ has strictly negative quadratic coefficients $-(1+t)u^2-(1-t)v^2$. Taking the Bargmann-Fock space inner product of the two generating functions and using the orthonormal basis proved in part (b) givesThe coefficient series converges absolutely by the Cauchy-Schwarz inequality. Substitution factors the right side intoUse the Gaussian integral $\int e^{-a{u}^2+bu}du=\sqrt{\pi /a}{e}^{b^2/(4a)}$ for $a>0$ and complex $b$; the complex case follows from the real case by analyticity. The resulting exponent isThus the Mehler formula for Hermite functions isAs a normalization check, its $t\downarrow 0$ limit is $H_0(x)H_0(y)$. Setting $t=e^{-2s}$ and multiplying by $e^{-s}$ yields the harmonic oscillator transition kernel with eigenvalues $2n+1$.

For a closed manifold $M$ with an oriented Riemannian metric, define $\delta$ to be the formal $L^2$ adjoint operator of the exterior derivative, equivalently $\delta =(-1{)}^{n(q+1)+1}\ast d\ast$ on $q$-forms. The Hodge Laplacian is ${\Delta }_q=d\delta +\delta d$. The Hodge decomposition theorem states that its space ${\mathcal{H}}^q=\ker {\Delta }_q$ of harmonic differential forms is finite-dimensional and smooth, and thatThe last isomorphism sends a harmonic differential form to its de Rham cohomology class: every class has exactly one harmonic representative. Out-of-range degree spaces are zero.

We first justify the analytic facts needed in the proof. In local coordinates and a smooth frame for the form bundle, the principal symbol of ${\Delta }_q$ is $g^{ij}{\xi }_i{\xi }_{j}I$. This follows from the exterior/contraction identities for $d$ and $\delta$. Freezing the leading coefficients, the Fourier transform estimate for the constant positive quadratic form gives, for each nonnegative integer $s$,Here is how to pass to variable coefficients. On a sufficiently small coordinate ball the difference of leading coefficients has arbitrarily small supremum norm. After differentiating $s$ times, its top-order contribution is bounded by $\epsilon \Vert u{\Vert }_{H^{s+2}}$; derivatives of the coefficients and first-order terms contribute $C_s\Vert u{\Vert }_{H^{s+1}}$. A partition of unity gives the same types of terms from commutators with the cutoffs. Absorb the epsilon term and use Sobolev interpolation to absorb $\Vert u{\Vert }_{H^{s+1}}$ into a small multiple of $\Vert u{\Vert }_{H^{s+2}}$ plus $C\Vert u{\Vert }_2$. This proves the estimate.

The first-order energy form $Q(u,u)=\Vert du{\Vert }_2^2+\Vert \delta u{\Vert }_2^2$ satisfies $Q(u,u)+c\Vert u{\Vert }_2^2\ge c^{\prime }\Vert u{\Vert }_{H^1}^2$ for a sufficiently large $c>0$. This follows from the same frozen-symbol argument at first order, or by integrating the principal part and absorbing the lower-order terms. The Lax-Milgram theorem thus constructs a bounded inverse $R=({\Delta }_q+c{)}^{-1}:L^2\to H^1$ in the weak sense. To obtain $H^2$ regularity of this weak solution, use a local difference quotient ${\delta }_{h}u$ and test the weak equation with $-{\delta }_{-h}({\chi }^2{\delta }_{h}u)$. Positivity of $g^{ij}$ bounds $\Vert \chi {\delta }_h\nabla u{\Vert }_2^2$; derivatives of $\chi$ and difference quotients of the smooth coefficients give terms bounded by $C\Vert u{\Vert }_{H^1}\Vert \chi {\delta }_h\nabla u{\Vert }_2$ and $C\Vert u{\Vert }_{H^1}^2$. The forcing contributes at most $C\Vert f{\Vert }_2(\Vert \chi {\delta }_h\nabla u{\Vert }_2+\Vert u{\Vert }_{H^1})$. Absorption gives a bound independent of $h$, so the difference-quotient criterion yields $u\in H^2$. Differentiating the equation and repeating gives the higher estimates and smoothness for smooth forcing. Thus the use of elliptic regularity here has been justified rather than assumed.