Reading guide · Proof index

A smooth quadratic reduction with parameters

This companion proves a signed quadratic normal form with smooth parameters, its critical-point Jacobian and the finite-chart consequence for compact parameter sets. Two exact families illustrate curved coordinates, zero coordinate pivots and the limitations of determinant bounds.

Mathematical sources: Jiří Lebl, Basic Analysis II, theorem numbering from version 6.3 (15 May 2026), gives the inverse theorem and proof and the implicit theorem and proof, Theorems 8.5.1 and 8.5.6. Guillemin–Sternberg’s 13 January 2010 author draft, §13.14.3, pp. 455–457, presents the classical Morse lemma and its smooth parameter variant. The proof here uses scalar elimination, the smooth inverse upgrade and the integral Taylor formula.

The programme proof index locates the local arguments; the human-proof references identify external statements and proofs, which retain their own terms. The finite-dimensional companion gives P1–P5; the compact-integral companion gives mixed derivatives and parameter integration. Gaussian and quadratic arguments are in the analytic module.

Original programme exposition: CC0 to the extent rights exist. Cited human works are external mathematical sources and retain their own terms.

M1. Parameter Morse lemma and its Jacobian

Let U⊂RnU\subset\mathbb R^n and S⊂RpS\subset\mathbb R^p be open, and let ϕ∈C∞(U×S;R)\phi\in C^\infty(U\times S;\mathbb R). Suppose

Dxϕ(x0,s0)=0,H0=Dx2ϕ(x0,s0) is invertible. D_x\phi(x_0,s_0)=0,\qquad H_0=D_x^2\phi(x_0,s_0) \text{ is invertible}.

There are neighbourhoods of s0s_0 and of z=0z=0, a smooth critical branch xc(s)x_c(s), and a smooth map F(z,s)F(z,s) such that for each nearby ss, z↦F(z,s)z\mapsto F(z,s) is a local diffeomorphism and

F(0,s)=xc(s),ϕ(F(z,s),s)=c(s)+12∑j=1nσjzj2,c(s)=ϕ(xc(s),s),σj∈{1,−1}.(M1.1) F(0,s)=x_c(s),\qquad \phi(F(z,s),s)=c(s)+\tfrac12\sum_{j=1}^n\sigma_j z_j^2, \qquad c(s)=\phi(x_c(s),s),\quad \sigma_j\in\{1,-1\}. \tag{M1.1}

The signs are independent of ss on this neighbourhood. If H(s)=Dx2ϕ(xc(s),s)H(s)=D_x^2\phi(x_c(s),s), then

∑jσj=sgn⁡H(s),∣det⁡DzF(0,s)∣=∣det⁡H(s)∣−1/2.(M1.2) \sum_j\sigma_j=\operatorname{sgn}H(s),\qquad |\det D_zF(0,s)|=|\det H(s)|^{-1/2}. \tag{M1.2}

Proof inputs. The freely accessible C1C^1 inverse and implicit proofs linked above supply the base theorems. Their smooth parameter upgrade P3, matrix inverse P2, Schur calculation P4 and congruence and determinant argument P5 are written locally. The mixed-partial identity P11.2, Taylor formula P12.6 and joint parameter integration P12.7 give the differential and integral steps. The positive square root P8 supplies existence and smoothness; finite compactness supplies the compact-parameter step.

For n=0n=0, the coordinate space is a singleton. Take xc(s)=x0x_c(s)=x_0, F(0,s)=x0F(0,s)=x_0, and c(s)=ϕ(x0,s)c(s)=\phi(x_0,s). The quadratic sum and signature are 0, while the empty determinant and its Jacobian factor are 1. This proves all assertions in that case. In the proof below, n≥1n\geq1.

Proof: locate and translate the critical point. Apply the implicit theorem, upgraded by P3, to the smooth map Dxϕ(x,s)D_x\phi(x,s). Its derivative in xx at (x0,s0)(x_0,s_0) is H0H_0, so there is a unique smooth branch xc(s)x_c(s) in a sufficiently small fixed product neighbourhood. Replace xx by xc(s)+ux_c(s)+u and subtract c(s)c(s). The resulting smooth function f(u,s)f(u,s) has f(0,s)=0f(0,s)=0, Duf(0,s)=0D_uf(0,s)=0 and an invertible Hessian for every nearby ss. The last assertion follows by continuity of the determinant and its nonzero value at s0s_0.

Choose a nonzero scalar pivot. A nonzero symmetric matrix has a vector with nonzero quadratic value. To verify this without assuming a normal form, if a diagonal entry is nonzero, use its coordinate vector. If every diagonal entry is zero, some off-diagonal entry hijh_{ij} is nonzero, and (ei+ej)TH(ei+ej)=2hij≠0(e_i+e_j)^TH(e_i+e_j)=2h_{ij}\ne0. Make a fixed invertible linear change of variables taking that vector to the first coordinate direction. Such a change can be constructed by adjoining all coordinate vectors except one in which the chosen vector has a nonzero component; expansion of the resulting matrix in that row gives a nonzero determinant. This change is fixed at s0s_0, so it is smooth in the parameters. After shrinking the neighbourhood, ∂12f(0,s)\partial_1^2 f(0,s) stays nonzero with one fixed sign.

Eliminate one variable. Write u=(u1,u′)u=(u_1,u'). P3 applied to ∂1f\partial_1 f gives a smooth function u1=q(u′,s)u_1=q(u',s) for which ∂1f(q(u′,s),u′,s)=0\partial_1 f(q(u',s),u',s)=0. Since u=0u=0 is critical and the solution is locally unique, q(0,s)=0q(0,s)=0. Put g(u′,s)=f(q(u′,s),u′,s)g(u',s)=f(q(u',s),u',s). With v=u1−q(u′,s)v=u_1-q(u',s), the integral Taylor formula proved in P12.6 gives

f(q(u′,s)+v,u′,s)=g(u′,s)+12v2A(v,u′,s), f(q(u',s)+v,u',s)=g(u',s)+\tfrac12v^2 A(v,u',s),
A(v,u′,s)=2∫01(1−t)∂12f(q(u′,s)+tv,u′,s) dt.(M1.3) A(v,u',s)=2\int_0^1(1-t) \partial_1^2 f(q(u',s)+tv,u',s)\,dt. \tag{M1.3}

Shrink the product neighbourhood so that q(u′,s)q(u',s) and q(u′,s)+vq(u',s)+v lie in a smaller coordinate interval whose closure stays inside the phase domain, with the other coordinates fixed. The whole segment between them then stays in that interval. Specifically apply P12.6 to the function t↦f(q(u′,s)+tv,u′,s)t\mapsto f(q(u',s)+tv,u',s) on [0,1][0,1]. Its first derivative at 0 is zero and its second derivative is v2∂12f(q(u′,s)+tv,u′,s)v^2\partial_1^2f(q(u',s)+tv,u',s). This proves (M1.3) for every sign of vv, including zero. The parameter integration is over a fixed compact interval, so P12.7 proves that AA is jointly smooth. No global improper-integral input from Q1 is needed for this step. At (v,u′)=0(v,u')=0, A(0,0,s)=∂12f(0,s)≠0A(0,0,s)=\partial_1^2 f(0,s)\ne0.

Shrink once more so that AA is nonzero throughout the neighbourhood with a single sign σ1\sigma_1. Define

z1=v∣A(v,u′,s)∣. z_1=v\sqrt{|A(v,u',s)|}.

Its vv-derivative at (v,u′)=(0,0)(v,u')=(0,0) is ∣A(0,0,s)∣>0\sqrt{|A(0,0,s)|}>0. Apply P3 to the map (v,u′,s)↦(z1,u′,s)(v,u',s)\mapsto(z_1,u',s): its derivative is block triangular, with this nonzero entry and identity blocks, hence invertible. We obtain a smooth local inverse, jointly in z1,u′,sz_1,u',s. In these coordinates,

f=g(u′,s)+σ1z12/2.(M1.4) f=g(u',s)+\sigma_1z_1^2/2.\tag{M1.4}

In particular the residual function is independent of z1z_1.

Induct with an invertible residual Hessian. The function gg has g(0,s)=0g(0,s)=0, and its first derivative vanishes at zero: differentiating the composite contributes only first derivatives of ff there. P4 computes its Hessian as the Schur complement of the chosen scalar pivot. That complement is invertible because the original Hessian is invertible. For n=1n=1, (M1.4) already completes the construction. For n>1n>1, apply the same construction to g(u′,s)g(u',s), in n−1n-1 variables and with the same parameters. Finite induction yields a signed quadratic form in all nn coordinates. At each of the finitely many steps, the nonzero pivot and the inverse map persist on an open neighbourhood of the relevant base point. Take preimages of these neighbourhoods under the preceding continuous coordinate maps, so they are neighbourhoods of the same original base point. Their finite intersection is still a neighbourhood; choose smaller product neighbourhoods inside it. Composing the inverse coordinate maps and the initial translation gives the smooth map FF in (M1.1).

Identify its signature and volume factor. Differentiate (M1.1) twice in zz at zero. Terms containing second derivatives of FF are multiplied by Dxϕ(xc(s),s)=0D_x\phi(x_c(s),s)=0, so the chain rule reduces to

[DzF(0,s)]TH(s)DzF(0,s)=diag⁡(σ1,…,σn). [D_zF(0,s)]^T H(s)D_zF(0,s) =\operatorname{diag}(\sigma_1,\ldots,\sigma_n).

P5 gives both the signature identity and the determinant identity in (M1.2). This also proves local constancy of the Hessian signature for this critical family; it has not been assumed as an unproved spectral-continuity assertion. □\square

M2. The precise compact-parameter consequence

Suppose a smooth critical branch is defined on an open neighbourhood of a compact parameter set S0S_0, and its Hessian is invertible on S0S_0. For each point of S0S_0, M1 supplies a local parameter neighbourhood and a signed Morse chart. Before selecting the cover, choose around each parameter point a smaller open parameter neighbourhood with compact closure in its chart domain, and a smaller coordinate ball with compact closure in the coordinate domain. These smaller parameter neighbourhoods still cover S0S_0. P1 and the earlier compactness proofs supply a finite subcover. On each selected product closure, every fixed finite list of derivatives of the coordinate map is bounded by the earlier extreme-value theorem. Its image under (z,s)↦(F(z,s),s)(z,s)\mapsto(F(z,s),s) is compact, by the proved continuous-image theorem, and is contained in the inverse map's open domain. Derivatives of the inverse are bounded on that compact image by the same extreme-value theorem. The largest of finitely many bounds is finite.

Thus later stationary-phase estimates can be proved in finitely many local charts with uniform finite constants. This conclusion does not assert a single global signed frame, nor a uniform coordinate size for a collection of phases whose Hessians approach singularity. It also does not itself construct the partition of unity or identify the clean quotient density; those are additional steps in the full lesson.

M3. A curved family with an exact signed quadratic chart

For s>−1s>-1 put

ϕs(x,y)=12(x−s+y2)2−12(1+s)y2,Fs(z1,z2)=(s+z1−z221+s,z21+s).(M3.1) \phi_s(x,y)=\tfrac12(x-s+y^2)^2-\tfrac12(1+s)y^2, \qquad F_s(z_1,z_2)=\left(s+z_1-\frac{z_2^2}{1+s}, \frac{z_2}{\sqrt{1+s}}\right). \tag{M3.1}

This is an exact example of M1, not a truncation. The inverse coordinates are z1=x−s+y2z_1=x-s+y^2, z2=1+s yz_2=\sqrt{1+s}\,y, so FsF_s is a global smooth diffeomorphism and

ϕs(Fs(z))=(z12−z22)/2,det⁡DFs=(1+s)−1/2.(M3.2) \phi_s(F_s(z))=(z_1^2-z_2^2)/2,\qquad \det DF_s=(1+s)^{-1/2}.\tag{M3.2}

To check the critical point and Hessian, write u=x−s+y2u=x-s+y^2. Then ∂xϕs=u\partial_x\phi_s=u, ∂yϕs=2yu−(1+s)y\partial_y\phi_s=2yu-(1+s)y. Both vanish only at (x,y)=(s,0)(x,y)=(s,0). Differentiating once more gives

D2ϕs=(12y2y2u+4y2−(1+s)),D2ϕs(s,0)=diag⁡(1,−(1+s)). D^2\phi_s= \begin{pmatrix}1&2y\\2y&2u+4y^2-(1+s)\end{pmatrix},\qquad D^2\phi_s(s,0)=\operatorname{diag}(1,-(1+s)).

Its signature is zero and its determinant is −(1+s)-(1+s), which verifies the Jacobian predicted by (M1.2). The direct derivative of FsF_s is upper triangular with diagonal entries 1,(1+s)−1/21,(1+s)^{-1/2}, independently checking that factor. On −1/2≤s≤1/2-1/2\leq s\leq1/2 these maps have bounded derivatives of every fixed order on each fixed compact zz-set. As s↓−1s\downarrow-1, the Hessian loses invertibility and its determinant factor becomes unbounded, illustrating the limitation in M2. Indeed, for every B>0B>0, the inequality 0<1+s<B−20<1+s<B^{-2} gives (1+s)−1/2>B(1+s)^{-1/2}>B, by the positive-root comparison proved in P8.

An exact parameter Morse chart maps a quadratic saddle to curved coordinates

Figure 2. The left panel shows the zz-coordinate rectangle [−1.4,1.4]×[−1.1,1.1][-1.4,1.4]\times[-1.1,1.1]. The other panels show its exact images under FsF_s for s=−1/2s=-1/2 and s=1/2s=1/2. Light curves are images of the same coordinate grid. The three labelled phase levels are mapped by (M3.1), so they have exactly the same phase values by (M3.2). The black point is the critical point, (0,0)(0,0) in the first panel and (s,0)(s,0) in the others. Curves are numerical samples of these exact formulae, not numerical proofs. Reproducible source: figures/draw_parameter_morse.py. The example is constructed here from the parameter normal-form mechanism in M1, with the free human sources identified at the start of this module.

M4. A zero-diagonal pivot and a determinant bound

Exercise. For every real parameter ss, consider

ψs(x,y)=xy+s2y2,Hs=(011s),det⁡Hs=−1.(M4.1) \psi_s(x,y)=xy+\tfrac{s}{2}y^2,\qquad H_s=\begin{pmatrix}0&1\\1&s\end{pmatrix},\qquad \det H_s=-1. \tag{M4.1}

Find its critical point and a signed Morse chart. Explain why the two coordinate directions cannot supply the first pivot when s=0s=0, and why a uniform determinant bound does not give uniform chart derivatives on the whole parameter line. Throughout this exercise, ∥D∥2=sup⁡∣z∣=1∣Dz∣\|D\|_2=\sup_{|z|=1}|Dz| is the Euclidean operator norm.

Solution. The gradient is (y,x+sy)(y,x+sy), so the only critical point is the origin. At s=0s=0 both diagonal Hessian entries vanish, but (e1+e2)TH0(e1+e2)=2(e_1+e_2)^TH_0(e_1+e_2)=2. This is the nonzero direction used in M1. Its quadratic value for general ss is 2+s2+s; that particular pivot is suitable near zero and vanishes at s=−2s=-2, although the full Hessian is invertible there too.

Put u=x+sy/2u=x+sy/2 and v=yv=y. Then ψs=uv\psi_s=uv. Taking u=(z1+z2)/2u=(z_1+z_2)/\sqrt2 and v=(z1−z2)/2v=(z_1-z_2)/\sqrt2 gives the global, jointly smooth chart

Fs(z)=12((1−s/2)z1+(1+s/2)z2, z1−z2),ψs(Fs(z))=12(z12−z22).(M4.2) F_s(z)=\frac1{\sqrt2}\bigl((1-s/2)z_1+(1+s/2)z_2,\ z_1-z_2\bigr), \qquad \psi_s(F_s(z))=\tfrac12(z_1^2-z_2^2). \tag{M4.2}

Its inverse is z1=(x+(s/2+1)y)/2z_1=(x+(s/2+1)y)/\sqrt2, z2=(x+(s/2−1)y)/2z_2=(x+(s/2-1)y)/\sqrt2. Direct expansion gives det⁡DzFs=−1\det D_zF_s=-1, so its absolute Jacobian agrees with M1.2 for every ss. The Hessian has signature zero, since it is congruent to J=diag⁡(1,−1)J=\operatorname{diag}(1,-1).

There cannot be a uniform first-derivative bound for any such choice of Morse charts over all real ss. Indeed, let DD be the derivative at zero of any chart with target JJ. The differentiated normal form gives DTHsD=JD^TH_sD=J. Inverting and multiplying yields

Hs−1=(−s110)=DJDT. H_s^{-1}=\begin{pmatrix}-s&1\\1&0\end{pmatrix}=DJD^T.

For r=DTe1r=D^Te_1, the upper-left entry says −s=rTJr-s=r^TJr, whence ∣s∣≤∣r∣2|s|\leq |r|^2. The row rTr^T is nonzero because DD is invertible. Testing DD on the unit vector r/∣r∣r/|r| gives a first coordinate equal to ∣r∣|r|, and hence

∥D∥2≥∣r∣≥∣s∣.(M4.3) \|D\|_2\geq |r|\geq\sqrt{|s|}. \tag{M4.3}

For every proposed finite bound B>0B>0, choosing ∣s∣>B2|s|>B^2 contradicts it. Thus even the constant value ∣det⁡Hs∣=1|\det H_s|=1 is insufficient on a noncompact parameter set. On each compact parameter interval the explicit chart and its inverse do have bounded derivatives of each fixed order on compact coordinate sets, as M2 states.

For the displayed chart write Ds=DzFs=AsRD_s=D_zF_s=A_sR, where As=(1−s/201)A_s=\left(\begin{smallmatrix}1&-s/2\\0&1\end{smallmatrix}\right) and R=2−1/2(111−1)R=2^{-1/2}\left(\begin{smallmatrix}1&1\\1&-1\end{smallmatrix}\right) is orthogonal. The matrix AsTAsA_s^TA_s has determinant 11 and trace 2+s2/42+s^2/4. Its larger eigenvalue is 1+s2/8+∣s∣s2+16/81+s^2/8+|s|\sqrt{s^2+16}/8, by the quadratic formula. Orthogonal diagonalization Q5 shows that the square of the operator norm is this larger eigenvalue. Taking its positive square root gives

∥Ds∥2=s2+16+∣s∣4.(M4.4) \|D_s\|_2=\frac{\sqrt{s^2+16}+|s|}{4}. \tag{M4.4}

This calculation also fixes the meaning of the plotted norm; it is not an entrywise or Frobenius norm. □\square

A unit circle, its Morse-coordinate image and exact chart-norm bounds

Figure 3. The first two panels show the unit circle and its image under F4F_4, with the coordinate vectors and their images labelled. The last panel plots the exact norm M4.4 of the displayed chart, the lower bound M4.3 valid for every such chart, and the constant absolute Jacobian 11. Its parameter range is [−8,8][-8,8]. Curves are numerical samples of the exact formulae; the proofs are above. Reproducible source: parameter-morse-pivot-norms.py; vector illustration.

Worked exercise and Figure 3: GPT-6.1 Sol (OpenAI). Original programme contribution: CC0 to the extent rights exist. Human sources and plotting software retain their own terms.