Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Terminal variations, reciprocal phases and natural cone maps

Definition
BianchiIX_TerminalConeData

by trppppp · Nov 5, 2026 · Mathlib 0df444a (Lean v4.33.1)

bianchi-ix-era-statisticsdynamical-systems

For real variables p,s,a,b,cp,s,a,b,cp,s,a,b,c, set

K=a2+b2+c2−2(ab+bc+ca),S=p2+s2,C=S+K.K=a^2+b^2+c^2-2(ab+bc+ca),\qquad S=p^2+s^2,\qquad C=S+K.K=a2+b2+c2−2(ab+bc+ca),S=p2+s2,C=S+K.

The shear components and wall rates are

P=Kp−2a2+b2+c2+ab−2bc+ac,P=Kp-2a^2+b^2+c^2+ab-2bc+ac,P=Kp−2a2+b2+c2+ab−2bc+ac, Q=Ks+3(−b2+c2+ab−ac),Q=Ks+\sqrt{3}(-b^2+c^2+ab-ac),Q=Ks+3​(−b2+c2+ab−ac), r1=−S+2p,r2=−S−p+3s,r3=−S−p−3s.r_1=-S+2p,\qquad r_2=-S-p+\sqrt{3}s,\qquad r_3=-S-p-\sqrt{3}s.r1​=−S+2p,r2​=−S−p+3​s,r3​=−S−p−3​s.

These functions define the normalized vacuum Bianchi IX vector field (P,Q,r1a,r2b,r3c)(P,Q,r_1a,r_2b,r_3c)(P,Q,r1​a,r2​b,r3​c). For y=(m,u,v,w)∈R4y=(m,u,v,w)\in\mathbb R^4y=(m,u,v,w)∈R4, define

K(u,v,w)=u2+v2+w2−2(uv+vw+wu),κ=u(v+w−u),K(u,v,w)=u^2+v^2+w^2-2(uv+vw+wu),\qquad \kappa=u(v+w-u),K(u,v,w)=u2+v2+w2−2(uv+vw+wu),κ=u(v+w−u), Fm=−v2+w2+u(v−w)−mκ,Fu=12(1−m2−(v−w)2+u2)u,F_m=-v^2+w^2+u(v-w)-m\kappa,\qquad F_u=\tfrac12(1-m^2-(v-w)^2+u^2)u,Fm​=−v2+w2+u(v−w)−mκ,Fu​=21​(1−m2−(v−w)2+u2)u, Fv=(m−1−κ)v,Fw=(−m−1−κ)w,F=(Fm,Fu,Fv,Fw),F_v=(m-1-\kappa)v,\qquad F_w=(-m-1-\kappa)w,\qquad F=(F_m,F_u,F_v,F_w),Fv​=(m−1−κ)v,Fw​=(−m−1−κ)w,F=(Fm​,Fu​,Fv​,Fw​),

and E(m,v)=(m−1)2+v2E(m,v)=(m-1)^2+v^2E(m,v)=(m−1)2+v2.

For x=(p,s,a,b,c)x=(p,s,a,b,c)x=(p,s,a,b,c), put

ψ(p,z)=3 z1+p,Π(x)=(ψ(p,s),ψ(p,a),ψ(p,b),ψ(p,c)).\psi(p,z)=\frac{\sqrt3\,z}{1+p},\qquad \Pi(x)=(\psi(p,s),\psi(p,a),\psi(p,b),\psi(p,c)).ψ(p,z)=1+p3​z​,Π(x)=(ψ(p,s),ψ(p,a),ψ(p,b),ψ(p,c)).

For y=(m,u,v,w)y=(m,u,v,w)y=(m,u,v,w), set D(y)=3+m2+K(u,v,w)D(y)=3+m^2+K(u,v,w)D(y)=3+m2+K(u,v,w) and

ΞA(y)=(A−1,Am3,Au3,Av3,Aw3),Ξ(y)=Ξ6/D(y)(y).\Xi_A(y)=\left(A-1,\frac{Am}{\sqrt3},\frac{Au}{\sqrt3}, \frac{Av}{\sqrt3},\frac{Aw}{\sqrt3}\right),\qquad \Xi(y)=\Xi_{6/D(y)}(y).ΞA​(y)=(A−1,3​Am​,3​Au​,3​Av​,3​Aw​),Ξ(y)=Ξ6/D(y)​(y).

Write N(x)=(P,Q,r1a,r2b,r3c)N(x)=(P,Q,r_1a,r_2b,r_3c)N(x)=(P,Q,r1​a,r2​b,r3​c) for the normalized vector field and C(x)=C(p,s,a,b,c)\mathcal C(x)=C(p,s,a,b,c)C(x)=C(p,s,a,b,c) for its constraint.

The normalized and projective state spaces are R5\mathbb R^5R5 and R4\mathbb R^4R4 with their finite-product norms. The maps Π,Ξ,D,F\Pi,\Xi,D,FΠ,Ξ,D,F act on these state spaces.

The projective vector field in normalized time is

G(y)=6D(y)F(y).G(y)=\frac6{D(y)}F(y).G(y)=D(y)6​F(y).

Write X=R5X=\mathbb R^5X=R5, Y=R4Y=\mathbb R^4Y=R4, and let αx,t0\alpha_{x,t_0}αx,t0​​ be the maximal solution of α′=N(α)\alpha'=N(\alpha)α′=N(α) with initial value xxx at t0t_0t0​ and domain Imax⁡(x,t0)I_{\max}(x,t_0)Imax​(x,t0​). Put

Jx,t0(t)=dxαx,t0(t),Ψt0(y,t)=Π(αΞ(y),t0(t)),J_{x,t_0}(t)=d_x\alpha_{x,t_0}(t),\qquad \Psi_{t_0}(y,t)=\Pi(\alpha_{\Xi(y),t_0}(t)),Jx,t0​​(t)=dx​αx,t0​​(t),Ψt0​​(y,t)=Π(αΞ(y),t0​​(t)), Jy,t0(t)=dΠαΞ(y),t0(t)∘JΞ(y),t0(t)∘dΞy.\mathcal J_{y,t_0}(t)=d\Pi_{\alpha_{\Xi(y),t_0}(t)} \circ J_{\Xi(y),t_0}(t)\circ d\Xi_y.Jy,t0​​(t)=dΠαΞ(y),t0​​(t)​∘JΞ(y),t0​​(t)∘dΞy​.

For Y=(k,u,w)∈R3Y=(k,u,w)\in\mathbb R^3Y=(k,u,w)∈R3, define

e−(Y)=(1−k2−1,u,1,w),xY=Ξ(e−(Y)),βY(t)=Ψ0(e−(Y),t),g(x)=(Π(x))v.e_-(Y)=(1-\sqrt{k^2-1},u,1,w),\quad x_Y=\Xi(e_-(Y)),\quad \beta_Y(t)=\Psi_0(e_-(Y),t),\quad g(x)=(\Pi(x))_v.e−​(Y)=(1−k2−1​,u,1,w),xY​=Ξ(e−​(Y)),βY​(t)=Ψ0​(e−​(Y),t),g(x)=(Π(x))v​.

For any real height η\etaη, put

Hη,Y={t>0:t∈Imax⁡(xY,0), g(αxY,0(t))=η},H_{\eta,Y}=\{t>0:t\in I_{\max}(x_Y,0),\ g(\alpha_{x_Y,0}(t))=\eta\},Hη,Y​={t>0:t∈Imax​(xY​,0), g(αxY​,0​(t))=η}, Tη(Y)=inf⁡Hη,Y,Pη(Y)=βY(Tη(Y)).T_\eta(Y)=\inf H_{\eta,Y},\qquad P_\eta(Y)=\beta_Y(T_\eta(Y)).Tη​(Y)=infHη,Y​,Pη​(Y)=βY​(Tη​(Y)).

The infimum uses the total real-order convention. Let TailValid⁡η(Y)\operatorname{TailValid}_\eta(Y)TailValidη​(Y) mean that Tη(Y)T_\eta(Y)Tη​(Y) is the least element of Hη,YH_{\eta,Y}Hη,Y​, and, writing T=Tη(Y)T=T_\eta(Y)T=Tη​(Y), that

t∈Imax⁡(xY,0),1+(αxY,0(t))0>0(0≤t≤T),t\in I_{\max}(x_Y,0),\quad 1+(\alpha_{x_Y,0}(t))_0>0 \quad(0\leq t\leq T),t∈Imax​(xY​,0),1+(αxY​,0​(t))0​>0(0≤t≤T), g(αxY,0(t))>η(0≤t<T),g(\alpha_{x_Y,0}(t))>\eta\quad(0\leq t<T),g(αxY​,0​(t))>η(0≤t<T), dgαxY,0(t)N(αxY,0(t))<0,∣(Π(αxY,0(t)))w∣<1/2(0≤t≤T).dg_{\alpha_{x_Y,0}(t)}N(\alpha_{x_Y,0}(t))<0,\quad |(\Pi(\alpha_{x_Y,0}(t)))_w|<1/2\quad(0\leq t\leq T).dgαxY​,0​(t)​N(αxY​,0​(t))<0,∣(Π(αxY​,0​(t)))w​∣<1/2(0≤t≤T).

All vector and operator norms use the finite-product supremum norms.

For y∈R4y\in\mathbb R^4y∈R4, put xy=Ξ(y)x_y=\Xi(y)xy​=Ξ(y),

T+(y)=inf⁡{t>0:t∈Imax⁡(xy,0), (Ψ0(y,t))v=1},P+(y)=Ψ0(y,T+(y)).T_+(y)=\inf\{t>0:t\in I_{\max}(x_y,0),\ (\Psi_0(y,t))_v=1\}, \qquad P_+(y)=\Psi_0(y,T_+(y)).T+​(y)=inf{t>0:t∈Imax​(xy​,0), (Ψ0​(y,t))v​=1},P+​(y)=Ψ0​(y,T+​(y)).

Let V+(y)V_+(y)V+​(y) mean that this infimum is a least element, the full interval [0,T+(y)][0,T_+(y)][0,T+​(y)] lies in the maximal domain with 1+(αxy,0(t))0>01+(\alpha_{x_y,0}(t))_0>01+(αxy​,0​(t))0​>0, the active wall is below one before the endpoint, and its normalized-time derivative is strictly positive throughout the interval. The infimum has its total real-order meaning.

For y=(m,u,v,w)y=(m,u,v,w)y=(m,u,v,w) define the cyclic chart by

H(y)=(m+1)2+K(u,v,w),A(y)=(1−4(m+1)/H,4v/H,4w/H,4u/H).H(y)=(m+1)^2+K(u,v,w),\qquad \mathcal A(y)=(1-4(m+1)/H,4v/H,4w/H,4u/H).H(y)=(m+1)2+K(u,v,w),A(y)=(1−4(m+1)/H,4v/H,4w/H,4u/H).

Its normalized symmetry is

A(p,s,a,b,c)=((−p+3s)/2,(−3p−s)/2,b,c,a).A(p,s,a,b,c)=((-p+\sqrt3s)/2,(-\sqrt3p-s)/2,b,c,a).A(p,s,a,b,c)=((−p+3​s)/2,(−3​p−s)/2,b,c,a).

Put zh(Y)=A(Ph(Y))z_h(Y)=\mathcal A(P_h(Y))zh​(Y)=A(Ph​(Y)) and

T(h,Y)=Th(Y)+T+(zh(Y)),P(h,Y)=P+(zh(Y)),xYc=AΞ(e−(Y)).\mathcal T(h,Y)=T_h(Y)+T_+(z_h(Y)),\qquad \mathcal P(h,Y)=P_+(z_h(Y)),\qquad x_Y^c=A\Xi(e_-(Y)).T(h,Y)=Th​(Y)+T+​(zh​(Y)),P(h,Y)=P+​(zh​(Y)),xYc​=AΞ(e−​(Y)).

Write r(m,v)=(m−1)2+v2r(m,v)=\sqrt{(m-1)^2+v^2}r(m,v)=(m−1)2+v2​ and g∗(x)=(Π(x))vg_*(x)=(\Pi(x))_vg∗​(x)=(Π(x))v​.

On the invariant plane with zero active wall, use y=(m,u,w)∈R3y=(m,u,w)\in\mathbb R^3y=(m,u,w)∈R3, the inclusion ι(m,u,w)=(m,u,0,w)\iota(m,u,w)=(m,u,0,w)ι(m,u,w)=(m,u,0,w) and the restriction πs(m,u,v,w)=(m,u,w)\pi_s(m,u,v,w)=(m,u,w)πs​(m,u,v,w)=(m,u,w). Define

f(y)=(mu2−(m+1)uw+w2,1−m2−w2+u22u,−(m+1+uw−u2)w),f(y)=\left(mu^2-(m+1)uw+w^2, \frac{1-m^2-w^2+u^2}{2}u,-(m+1+uw-u^2)w\right),f(y)=(mu2−(m+1)uw+w2,21−m2−w2+u2​u,−(m+1+uw−u2)w), Δ(y)=3+m2+(u−w)2,fo(y)=6f(y)/Δ(y).\Delta(y)=3+m^2+(u-w)^2,\qquad f_o(y)=6f(y)/\Delta(y).Δ(y)=3+m2+(u−w)2,fo​(y)=6f(y)/Δ(y).

Write J=dfJ=dfJ=df, H=d2fH=d^2fH=d2f, Jo=dfoJ_o=df_oJo​=dfo​, Ho=d2foH_o=d^2f_oHo​=d2fo​ for their explicit Jacobian and Hessian expressions. The denominator derivatives are dΔ=(2m,2(u−w),2(w−u))d\Delta=(2m,2(u-w),2(w-u))dΔ=(2m,2(u−w),2(w−u)) and d2Δ(v,z)=2v0z0+2(v1−v2)(z1−z2)d^2\Delta(v,z)=2v_0z_0+2(v_1-v_2)(z_1-z_2)d2Δ(v,z)=2v0​z0​+2(v1​−v2​)(z1​−z2​). The exact original-time boundary flow and its projected variation are

Φ(y,t)=πsΨ0(ιy,t),L(y,t)=πsJιy,0(t)ι.\Phi(y,t)=\pi_s\Psi_0(\iota y,t),\qquad L(y,t)=\pi_s\mathcal J_{\iota y,0}(t)\iota.Φ(y,t)=πs​Ψ0​(ιy,t),L(y,t)=πs​Jιy,0​(t)ι.

For a fixed scale KKK, set ΦK(y,s)=Φ(y,Ks)\Phi_K(y,s)=\Phi(y,Ks)ΦK​(y,s)=Φ(y,Ks), LK(y,s)=L(y,Ks)L_K(y,s)=L(y,Ks)LK​(y,s)=L(y,Ks) and QK(y,s)=dyLK(y,s)Q_K(y,s)=d_yL_K(y,s)QK​(y,s)=dy​LK​(y,s). For directions v,zv,zv,z, define

SK(y,s;v,z)=K[Ho(ΦK)(LKv,LKz)+Jo(ΦK)QK(v,z)].S_K(y,s;v,z)=K\left[H_o(\Phi_K)(L_Kv,L_Kz)+J_o(\Phi_K)Q_K(v,z)\right].SK​(y,s;v,z)=K[Ho​(ΦK​)(LK​v,LK​z)+Jo​(ΦK​)QK​(v,z)].

All vector norms are sup norms, and operator norms are induced by them. Write b(y)=∣u∣+∣w∣b(y)=|u|+|w|b(y)=∣u∣+∣w∣, E(y)=e512b(y)E(y)=e^{512b(y)}E(y)=e512b(y), vs=(v1,v2)v_s=(v_1,v_2)vs​=(v1​,v2​), and

H(y;v,z)=b(y)∥v∥∥z∥+∥vs∥∥z∥+∥v∥∥zs∥,\mathcal H(y;v,z)=b(y)\|v\|\|z\|+\|v_s\|\|z\|+\|v\|\|z_s\|,H(y;v,z)=b(y)∥v∥∥z∥+∥vs​∥∥z∥+∥v∥∥zs​∥, B2(y)=192b(y)(40000E(y)3)+10000E(y)2.B_2(y)=192b(y)(40000E(y)^3)+10000E(y)^2.B2​(y)=192b(y)(40000E(y)3)+10000E(y)2.

The closed boundary region is Sc={(m,u,0,w):∣m−c∣≤1/10, ∣u∣,∣w∣≤1/100}\mathcal S_c=\{(m,u,0,w):|m-c|\leq1/10,\ |u|,|w|\leq1/100\}Sc​={(m,u,0,w):∣m−c∣≤1/10, ∣u∣,∣w∣≤1/100}. The boundary projection P(m,u,w)P(m,u,w)P(m,u,w) is the limit of Φ(m,u,w;t)0\Phi(m,u,w;t)_0Φ(m,u,w;t)0​ as t→∞t\to\inftyt→∞, using the total limit convention of the formal definitions.

Set p(y)=P(y0,y1,y3)p(y)=P(y_0,y_1,y_3)p(y)=P(y0​,y1​,y3​), ρ(y)=(p(y)−1)2+y22\rho(y)=\sqrt{(p(y)-1)^2+y_2^2}ρ(y)=(p(y)−1)2+y22​​, and χ(y)=(ρ(y),y1,y2,y3)\chi(y)=(\rho(y),y_1,y_2,y_3)χ(y)=(ρ(y),y1​,y2​,y3​). An admissible inverse ψ\psiψ at yyy is C2C^2C2 at χ(y)\chi(y)χ(y), maps χ(y)\chi(y)χ(y) to yyy, and is a two-sided local inverse of χ\chiχ; its derivative is the inverse of the triangular derivative of χ\chiχ. In the following notation put yh=A∘Phy_h=\mathcal A\circ P_hyh​=A∘Ph​ and zh=χ∘yhz_h=\chi\circ y_hzh​=χ∘yh​. Define

Fh(Y)=r(P(h,Y)0,P(h,Y)2),Dh=Fh−(zh)0,Qh=2/Fh−2/(zh)0.F_h(Y)=r(\mathcal P(h,Y)_0,\mathcal P(h,Y)_2),\quad D_h=F_h-(z_h)_0,\quad Q_h=2/F_h-2/(z_h)_0.Fh​(Y)=r(P(h,Y)0​,P(h,Y)2​),Dh​=Fh​−(zh​)0​,Qh​=2/Fh​−2/(zh​)0​.

For I(m,u)=(1+m2+u2)/mI(m,u)=(1+m^2+u^2)/mI(m,u)=(1+m2+u2)/m let

λ+(I)=(I+I2−4)/2,s(m,u)=2/(λ+(I(m,u))−1).\lambda_+(I)=(I+\sqrt{I^2-4})/2,\qquad s(m,u)=2/(\lambda_+(I(m,u))-1).λ+​(I)=(I+I2−4​)/2,s(m,u)=2/(λ+​(I(m,u))−1).

Set Sh=s((yh)0,(yh)1)S_h=s((y_h)_0,(y_h)_1)Sh​=s((yh​)0​,(yh​)1​), Qhs=2/(zh)0−ShQ_h^s=2/(z_h)_0-S_hQhs​=2/(zh​)0​−Sh​, Eh=2/Fh−ShE_h=2/F_h-S_hEh​=2/Fh​−Sh​, and Bh=2/FhB_h=2/F_hBh​=2/Fh​. The positive wall multipliers m0,m1m_0,m_1m0​,m1​ are defined by integrating the two multiplicative wall rates along the exact exit flow, so that (Ph)1=um0(P_h)_1=u m_0(Ph​)1​=um0​ and (Ph)3=wm1(P_h)_3=w m_1(Ph​)3​=wm1​. Write Bih=4mi/H(Ph)B_i^h=4m_i/H(P_h)Bih​=4mi​/H(Ph​) for i=0,1i=0,1i=0,1. For Y=(2+2α,u,w)Y=(2+2\alpha,u,w)Y=(2+2α,u,w) let h=ηαJh=\eta\alpha^Jh=ηαJ. All derivatives of maps with subscript hhh hold this height fixed. A cone direction is d=(1,d1,d2)d=(1,d_1,d_2)d=(1,d1​,d2​) with ∣d1∣≤u|d_1|\leq\sqrt u∣d1​∣≤u​, ∣d2∣≤w|d_2|\leq\sqrt w∣d2​∣≤w​. For y=(m,u,v,w)∈R4y=(m,u,v,w)\in\mathbb R^4y=(m,u,v,w)∈R4, define

K(u,v,w)=u2+v2+w2−2(uv+vw+wu),κ=u(v+w−u),K(u,v,w)=u^2+v^2+w^2-2(uv+vw+wu),\qquad \kappa=u(v+w-u),K(u,v,w)=u2+v2+w2−2(uv+vw+wu),κ=u(v+w−u), Fm=−v2+w2+u(v−w)−mκ,Fu=12(1−m2−(v−w)2+u2)u,F_m=-v^2+w^2+u(v-w)-m\kappa,\qquad F_u=\tfrac12(1-m^2-(v-w)^2+u^2)u,Fm​=−v2+w2+u(v−w)−mκ,Fu​=21​(1−m2−(v−w)2+u2)u, Fv=(m−1−κ)v,Fw=(−m−1−κ)w,F=(Fm,Fu,Fv,Fw),F_v=(m-1-\kappa)v,\qquad F_w=(-m-1-\kappa)w,\qquad F=(F_m,F_u,F_v,F_w),Fv​=(m−1−κ)v,Fw​=(−m−1−κ)w,F=(Fm​,Fu​,Fv​,Fw​),

and E(m,v)=(m−1)2+v2E(m,v)=(m-1)^2+v^2E(m,v)=(m−1)2+v2.

Write r(m,v)=(m−1)2+v2r(m,v)=\sqrt{(m-1)^2+v^2}r(m,v)=(m−1)2+v2​ and define

g(R,v)=R2−v2,d(R,u,v,w)=g(R,v)+u(v+w−u),g(R,v)=\sqrt{R^2-v^2},\qquad d(R,u,v,w)=g(R,v)+u(v+w-u),g(R,v)=R2−v2​,d(R,u,v,w)=g(R,v)+u(v+w−u), fR=−u(v+w−u)R−g(R,v)(u−w)2/Rd,f_R=\frac{-u(v+w-u)R-g(R,v)(u-w)^2/R}{d},fR​=d−u(v+w−u)R−g(R,v)(u−w)2/R​, as=2g(R,v)−R2+2vw−w2+u22d,aw=g(R,v)−2−u(v+w−u)d.a_s=\frac{2g(R,v)-R^2+2vw-w^2+u^2}{2d},\qquad a_w=\frac{g(R,v)-2-u(v+w-u)}d.as​=2d2g(R,v)−R2+2vw−w2+u2​,aw​=dg(R,v)−2−u(v+w−u)​.

The exit clock rate is ω(y)=1−m+u(v+w−u)\omega(y)=1-m+u(v+w-u)ω(y)=1−m+u(v+w−u), and its field is Flog⁡(y)=F(y)/ω(y)F_{\log}(y)=F(y)/\omega(y)Flog​(y)=F(y)/ω(y). Division and square roots have their total real-valued Lean meanings.

For Y=(R,u,w)Y=(R,u,w)Y=(R,u,w) define Ev(Y)=(fR,asu,aww)E_v(Y)=(f_R,a_su,a_ww)Ev​(Y)=(fR​,as​u,aw​w), evaluated at (R,u,v,w)(R,u,v,w)(R,u,v,w). The coordinate maps are Q(m,u,v,w)=(r(m,v),u,w)Q(m,u,v,w)=(r(m,v),u,w)Q(m,u,v,w)=(r(m,v),u,w) and Qv−1(R,u,w)=(1−g(R,v),u,v,w)Q_v^{-1}(R,u,w)=(1-g(R,v),u,v,w)Qv−1​(R,u,w)=(1−g(R,v),u,v,w). For a curve Y(s)Y(s)Y(s) set

Mi(Y,t)=exp⁡∫0tai(Y0(s),Y1(s),e−s,Y2(s)) ds(i=s,w).M_i(Y,t)=\exp\int_0^t a_i(Y_0(s),Y_1(s),e^{-s},Y_2(s))\,ds \quad(i=s,w).Mi​(Y,t)=exp∫0t​ai​(Y0​(s),Y1​(s),e−s,Y2​(s))ds(i=s,w).

Let Rα,b,A,D\mathcal R_{\alpha,b,A,D}Rα,b,A,D​ consist of the states satisfying

7/4≤R≤17/4,∣u∣,∣w∣≤1/100,∣u∣+∣w∣≤Ab,∣R−(2+2α)∣≤Db.7/4\leq R\leq17/4,\quad |u|,|w|\leq1/100,\quad |u|+|w|\leq Ab,\quad |R-(2+2\alpha)|\leq Db.7/4≤R≤17/4,∣u∣,∣w∣≤1/100,∣u∣+∣w∣≤Ab,∣R−(2+2α)∣≤Db.

Put

B(α,b,A,D,T)=[10(α+Db)+100(Ab)2]T+100Ab+5,B(\alpha,b,A,D,T)=[10(\alpha+Db)+100(Ab)^2]T+100Ab+5,B(α,b,A,D,T)=[10(α+Db)+100(Ab)2]T+100Ab+5, L(η,α,J)=log⁡(1/η)+Jlog⁡(1/α).L(\eta,\alpha,J)=\log(1/\eta)+J\log(1/\alpha).L(η,α,J)=log(1/η)+Jlog(1/α).

For real variables p,s,a,b,cp,s,a,b,cp,s,a,b,c, set

K=a2+b2+c2−2(ab+bc+ca),S=p2+s2,C=S+K.K=a^2+b^2+c^2-2(ab+bc+ca),\qquad S=p^2+s^2,\qquad C=S+K.K=a2+b2+c2−2(ab+bc+ca),S=p2+s2,C=S+K.

The shear components and wall rates are

P=Kp−2a2+b2+c2+ab−2bc+ac,P=Kp-2a^2+b^2+c^2+ab-2bc+ac,P=Kp−2a2+b2+c2+ab−2bc+ac, Q=Ks+3(−b2+c2+ab−ac),Q=Ks+\sqrt{3}(-b^2+c^2+ab-ac),Q=Ks+3​(−b2+c2+ab−ac), r1=−S+2p,r2=−S−p+3s,r3=−S−p−3s.r_1=-S+2p,\qquad r_2=-S-p+\sqrt{3}s,\qquad r_3=-S-p-\sqrt{3}s.r1​=−S+2p,r2​=−S−p+3​s,r3​=−S−p−3​s.

These functions define the normalized vacuum Bianchi IX vector field (P,Q,r1a,r2b,r3c)(P,Q,r_1a,r_2b,r_3c)(P,Q,r1​a,r2​b,r3​c). For y=(m,u,v,w)∈R4y=(m,u,v,w)\in\mathbb R^4y=(m,u,v,w)∈R4, define

K(u,v,w)=u2+v2+w2−2(uv+vw+wu),κ=u(v+w−u),K(u,v,w)=u^2+v^2+w^2-2(uv+vw+wu),\qquad \kappa=u(v+w-u),K(u,v,w)=u2+v2+w2−2(uv+vw+wu),κ=u(v+w−u), Fm=−v2+w2+u(v−w)−mκ,Fu=12(1−m2−(v−w)2+u2)u,F_m=-v^2+w^2+u(v-w)-m\kappa,\qquad F_u=\tfrac12(1-m^2-(v-w)^2+u^2)u,Fm​=−v2+w2+u(v−w)−mκ,Fu​=21​(1−m2−(v−w)2+u2)u, Fv=(m−1−κ)v,Fw=(−m−1−κ)w,F=(Fm,Fu,Fv,Fw),F_v=(m-1-\kappa)v,\qquad F_w=(-m-1-\kappa)w,\qquad F=(F_m,F_u,F_v,F_w),Fv​=(m−1−κ)v,Fw​=(−m−1−κ)w,F=(Fm​,Fu​,Fv​,Fw​),

and E(m,v)=(m−1)2+v2E(m,v)=(m-1)^2+v^2E(m,v)=(m−1)2+v2.

For x=(p,s,a,b,c)x=(p,s,a,b,c)x=(p,s,a,b,c), put

ψ(p,z)=3 z1+p,Π(x)=(ψ(p,s),ψ(p,a),ψ(p,b),ψ(p,c)).\psi(p,z)=\frac{\sqrt3\,z}{1+p},\qquad \Pi(x)=(\psi(p,s),\psi(p,a),\psi(p,b),\psi(p,c)).ψ(p,z)=1+p3​z​,Π(x)=(ψ(p,s),ψ(p,a),ψ(p,b),ψ(p,c)).

For y=(m,u,v,w)y=(m,u,v,w)y=(m,u,v,w), set D(y)=3+m2+K(u,v,w)D(y)=3+m^2+K(u,v,w)D(y)=3+m2+K(u,v,w) and

ΞA(y)=(A−1,Am3,Au3,Av3,Aw3),Ξ(y)=Ξ6/D(y)(y).\Xi_A(y)=\left(A-1,\frac{Am}{\sqrt3},\frac{Au}{\sqrt3}, \frac{Av}{\sqrt3},\frac{Aw}{\sqrt3}\right),\qquad \Xi(y)=\Xi_{6/D(y)}(y).ΞA​(y)=(A−1,3​Am​,3​Au​,3​Av​,3​Aw​),Ξ(y)=Ξ6/D(y)​(y).

Write N(x)=(P,Q,r1a,r2b,r3c)N(x)=(P,Q,r_1a,r_2b,r_3c)N(x)=(P,Q,r1​a,r2​b,r3​c) for the normalized vector field and C(x)=C(p,s,a,b,c)\mathcal C(x)=C(p,s,a,b,c)C(x)=C(p,s,a,b,c) for its constraint.

The normalized and projective state spaces are R5\mathbb R^5R5 and R4\mathbb R^4R4 with their finite-product norms. The maps Π,Ξ,D,F\Pi,\Xi,D,FΠ,Ξ,D,F act on these state spaces.

The projective vector field in normalized time is

G(y)=6D(y)F(y).G(y)=\frac6{D(y)}F(y).G(y)=D(y)6​F(y).

Write X=R5X=\mathbb R^5X=R5, Y=R4Y=\mathbb R^4Y=R4, and let αx,t0\alpha_{x,t_0}αx,t0​​ be the maximal solution of α′=N(α)\alpha'=N(\alpha)α′=N(α) with initial value xxx at t0t_0t0​ and domain Imax⁡(x,t0)I_{\max}(x,t_0)Imax​(x,t0​). Put

Jx,t0(t)=dxαx,t0(t),Ψt0(y,t)=Π(αΞ(y),t0(t)),J_{x,t_0}(t)=d_x\alpha_{x,t_0}(t),\qquad \Psi_{t_0}(y,t)=\Pi(\alpha_{\Xi(y),t_0}(t)),Jx,t0​​(t)=dx​αx,t0​​(t),Ψt0​​(y,t)=Π(αΞ(y),t0​​(t)), Jy,t0(t)=dΠαΞ(y),t0(t)∘JΞ(y),t0(t)∘dΞy.\mathcal J_{y,t_0}(t)=d\Pi_{\alpha_{\Xi(y),t_0}(t)} \circ J_{\Xi(y),t_0}(t)\circ d\Xi_y.Jy,t0​​(t)=dΠαΞ(y),t0​​(t)​∘JΞ(y),t0​​(t)∘dΞy​.

For Y=(k,u,w)∈R3Y=(k,u,w)\in\mathbb R^3Y=(k,u,w)∈R3, define

e−(Y)=(1−k2−1,u,1,w),xY=Ξ(e−(Y)),βY(t)=Ψ0(e−(Y),t),g(x)=(Π(x))v.e_-(Y)=(1-\sqrt{k^2-1},u,1,w),\quad x_Y=\Xi(e_-(Y)),\quad \beta_Y(t)=\Psi_0(e_-(Y),t),\quad g(x)=(\Pi(x))_v.e−​(Y)=(1−k2−1​,u,1,w),xY​=Ξ(e−​(Y)),βY​(t)=Ψ0​(e−​(Y),t),g(x)=(Π(x))v​.

For any real height η\etaη, put

Hη,Y={t>0:t∈Imax⁡(xY,0), g(αxY,0(t))=η},H_{\eta,Y}=\{t>0:t\in I_{\max}(x_Y,0),\ g(\alpha_{x_Y,0}(t))=\eta\},Hη,Y​={t>0:t∈Imax​(xY​,0), g(αxY​,0​(t))=η}, Tη(Y)=inf⁡Hη,Y,Pη(Y)=βY(Tη(Y)).T_\eta(Y)=\inf H_{\eta,Y},\qquad P_\eta(Y)=\beta_Y(T_\eta(Y)).Tη​(Y)=infHη,Y​,Pη​(Y)=βY​(Tη​(Y)).

The infimum uses the total real-order convention. Let TailValid⁡η(Y)\operatorname{TailValid}_\eta(Y)TailValidη​(Y) mean that Tη(Y)T_\eta(Y)Tη​(Y) is the least element of Hη,YH_{\eta,Y}Hη,Y​, and, writing T=Tη(Y)T=T_\eta(Y)T=Tη​(Y), that

t∈Imax⁡(xY,0),1+(αxY,0(t))0>0(0≤t≤T),t\in I_{\max}(x_Y,0),\quad 1+(\alpha_{x_Y,0}(t))_0>0 \quad(0\leq t\leq T),t∈Imax​(xY​,0),1+(αxY​,0​(t))0​>0(0≤t≤T), g(αxY,0(t))>η(0≤t<T),g(\alpha_{x_Y,0}(t))>\eta\quad(0\leq t<T),g(αxY​,0​(t))>η(0≤t<T), dgαxY,0(t)N(αxY,0(t))<0,∣(Π(αxY,0(t)))w∣<1/2(0≤t≤T).dg_{\alpha_{x_Y,0}(t)}N(\alpha_{x_Y,0}(t))<0,\quad |(\Pi(\alpha_{x_Y,0}(t)))_w|<1/2\quad(0\leq t\leq T).dgαxY​,0​(t)​N(αxY​,0​(t))<0,∣(Π(αxY​,0​(t)))w​∣<1/2(0≤t≤T).

All vector and operator norms use the finite-product supremum norms.

For y∈R4y\in\mathbb R^4y∈R4, put xy=Ξ(y)x_y=\Xi(y)xy​=Ξ(y),

T+(y)=inf⁡{t>0:t∈Imax⁡(xy,0), (Ψ0(y,t))v=1},P+(y)=Ψ0(y,T+(y)).T_+(y)=\inf\{t>0:t\in I_{\max}(x_y,0),\ (\Psi_0(y,t))_v=1\}, \qquad P_+(y)=\Psi_0(y,T_+(y)).T+​(y)=inf{t>0:t∈Imax​(xy​,0), (Ψ0​(y,t))v​=1},P+​(y)=Ψ0​(y,T+​(y)).

Let V+(y)V_+(y)V+​(y) mean that this infimum is a least element, the full interval [0,T+(y)][0,T_+(y)][0,T+​(y)] lies in the maximal domain with 1+(αxy,0(t))0>01+(\alpha_{x_y,0}(t))_0>01+(αxy​,0​(t))0​>0, the active wall is below one before the endpoint, and its normalized-time derivative is strictly positive throughout the interval. The infimum has its total real-order meaning.

For y=(m,u,v,w)y=(m,u,v,w)y=(m,u,v,w) define the cyclic chart by

H(y)=(m+1)2+K(u,v,w),A(y)=(1−4(m+1)/H,4v/H,4w/H,4u/H).H(y)=(m+1)^2+K(u,v,w),\qquad \mathcal A(y)=(1-4(m+1)/H,4v/H,4w/H,4u/H).H(y)=(m+1)2+K(u,v,w),A(y)=(1−4(m+1)/H,4v/H,4w/H,4u/H).

Its normalized symmetry is

A(p,s,a,b,c)=((−p+3s)/2,(−3p−s)/2,b,c,a).A(p,s,a,b,c)=((-p+\sqrt3s)/2,(-\sqrt3p-s)/2,b,c,a).A(p,s,a,b,c)=((−p+3​s)/2,(−3​p−s)/2,b,c,a).

Put zh(Y)=A(Ph(Y))z_h(Y)=\mathcal A(P_h(Y))zh​(Y)=A(Ph​(Y)) and

T(h,Y)=Th(Y)+T+(zh(Y)),P(h,Y)=P+(zh(Y)),xYc=AΞ(e−(Y)).\mathcal T(h,Y)=T_h(Y)+T_+(z_h(Y)),\qquad \mathcal P(h,Y)=P_+(z_h(Y)),\qquad x_Y^c=A\Xi(e_-(Y)).T(h,Y)=Th​(Y)+T+​(zh​(Y)),P(h,Y)=P+​(zh​(Y)),xYc​=AΞ(e−​(Y)).

Write r(m,v)=(m−1)2+v2r(m,v)=\sqrt{(m-1)^2+v^2}r(m,v)=(m−1)2+v2​ and g∗(x)=(Π(x))vg_*(x)=(\Pi(x))_vg∗​(x)=(Π(x))v​.

Let β(Y,t)=Ψ0(e−(Y),t)\beta(Y,t)=\Psi_0(e_-(Y),t)β(Y,t)=Ψ0​(e−​(Y),t) and c(Y,t)=−log⁡β(Y,t)2c(Y,t)=-\log\beta(Y,t)_2c(Y,t)=−logβ(Y,t)2​. Call (Y,t)(Y,t)(Y,t) regular when Y0>1Y_0>1Y0​>1, D(e−(Y))>0D(e_-(Y))>0D(e−​(Y))>0, the normalized flow exists at ttt with positive chart clock, and β2>0\beta_2>0β2​>0, β0<1\beta_0<1β0​<1, ω(β)>0\omega(\beta)>0ω(β)>0. For an interval III, let τI(Y,s)\tau_I(Y,s)τI​(Y,s) be a solution in III of c(Y,t)=sc(Y,t)=sc(Y,t)=s, selected by choice, or zero if there is none. Define

XI(Y,s)=Q(β(Y,τI(Y,s))),ℓi(Y,s)=∫0sai(XI(Y,r)0,XI(Y,r)1,e−r,XI(Y,r)2) dr.X_I(Y,s)=Q(\beta(Y,\tau_I(Y,s))),\qquad \ell_i(Y,s)=\int_0^s a_i(X_I(Y,r)_0,X_I(Y,r)_1,e^{-r},X_I(Y,r)_2)\,dr.XI​(Y,s)=Q(β(Y,τI​(Y,s))),ℓi​(Y,s)=∫0s​ai​(XI​(Y,r)0​,XI​(Y,r)1​,e−r,XI​(Y,r)2​)dr.

Here i=0,1i=0,1i=0,1 denotes as,awa_s,a_was​,aw​ respectively. Write

V(α,b,A,D,T)=2B(α,b,A,D,T)+2048(1+AbT).V(\alpha,b,A,D,T)=2B(\alpha,b,A,D,T)+2048(1+AbT).V(α,b,A,D,T)=2B(α,b,A,D,T)+2048(1+AbT).

The tail box is 7/4≤R≤17/47/4\leq R\leq17/47/4≤R≤17/4, ∣u∣,∣w∣≤1/100|u|,|w|\leq1/100∣u∣,∣w∣≤1/100, 0≤v≤10\leq v\leq10≤v≤1. Write GR=∇R,u,wfRG_R=\nabla_{R,u,w}f_RGR​=∇R,u,w​fR​, Gs=∇R,u,wasG_s=\nabla_{R,u,w}a_sGs​=∇R,u,w​as​, Gw=∇R,u,wawG_w=\nabla_{R,u,w}a_wGw​=∇R,u,w​aw​ for their explicit algebraic gradients there. For real variables p,s,a,b,cp,s,a,b,cp,s,a,b,c, set

K=a2+b2+c2−2(ab+bc+ca),S=p2+s2,C=S+K.K=a^2+b^2+c^2-2(ab+bc+ca),\qquad S=p^2+s^2,\qquad C=S+K.K=a2+b2+c2−2(ab+bc+ca),S=p2+s2,C=S+K.

The shear components and wall rates are

P=Kp−2a2+b2+c2+ab−2bc+ac,P=Kp-2a^2+b^2+c^2+ab-2bc+ac,P=Kp−2a2+b2+c2+ab−2bc+ac, Q=Ks+3(−b2+c2+ab−ac),Q=Ks+\sqrt{3}(-b^2+c^2+ab-ac),Q=Ks+3​(−b2+c2+ab−ac), r1=−S+2p,r2=−S−p+3s,r3=−S−p−3s.r_1=-S+2p,\qquad r_2=-S-p+\sqrt{3}s,\qquad r_3=-S-p-\sqrt{3}s.r1​=−S+2p,r2​=−S−p+3​s,r3​=−S−p−3​s.

These functions define the normalized vacuum Bianchi IX vector field (P,Q,r1a,r2b,r3c)(P,Q,r_1a,r_2b,r_3c)(P,Q,r1​a,r2​b,r3​c). For y=(m,u,v,w)∈R4y=(m,u,v,w)\in\mathbb R^4y=(m,u,v,w)∈R4, define

K(u,v,w)=u2+v2+w2−2(uv+vw+wu),κ=u(v+w−u),K(u,v,w)=u^2+v^2+w^2-2(uv+vw+wu),\qquad \kappa=u(v+w-u),K(u,v,w)=u2+v2+w2−2(uv+vw+wu),κ=u(v+w−u), Fm=−v2+w2+u(v−w)−mκ,Fu=12(1−m2−(v−w)2+u2)u,F_m=-v^2+w^2+u(v-w)-m\kappa,\qquad F_u=\tfrac12(1-m^2-(v-w)^2+u^2)u,Fm​=−v2+w2+u(v−w)−mκ,Fu​=21​(1−m2−(v−w)2+u2)u, Fv=(m−1−κ)v,Fw=(−m−1−κ)w,F=(Fm,Fu,Fv,Fw),F_v=(m-1-\kappa)v,\qquad F_w=(-m-1-\kappa)w,\qquad F=(F_m,F_u,F_v,F_w),Fv​=(m−1−κ)v,Fw​=(−m−1−κ)w,F=(Fm​,Fu​,Fv​,Fw​),

and E(m,v)=(m−1)2+v2E(m,v)=(m-1)^2+v^2E(m,v)=(m−1)2+v2.

For x=(p,s,a,b,c)x=(p,s,a,b,c)x=(p,s,a,b,c), put

ψ(p,z)=3 z1+p,Π(x)=(ψ(p,s),ψ(p,a),ψ(p,b),ψ(p,c)).\psi(p,z)=\frac{\sqrt3\,z}{1+p},\qquad \Pi(x)=(\psi(p,s),\psi(p,a),\psi(p,b),\psi(p,c)).ψ(p,z)=1+p3​z​,Π(x)=(ψ(p,s),ψ(p,a),ψ(p,b),ψ(p,c)).

For y=(m,u,v,w)y=(m,u,v,w)y=(m,u,v,w), set D(y)=3+m2+K(u,v,w)D(y)=3+m^2+K(u,v,w)D(y)=3+m2+K(u,v,w) and

ΞA(y)=(A−1,Am3,Au3,Av3,Aw3),Ξ(y)=Ξ6/D(y)(y).\Xi_A(y)=\left(A-1,\frac{Am}{\sqrt3},\frac{Au}{\sqrt3}, \frac{Av}{\sqrt3},\frac{Aw}{\sqrt3}\right),\qquad \Xi(y)=\Xi_{6/D(y)}(y).ΞA​(y)=(A−1,3​Am​,3​Au​,3​Av​,3​Aw​),Ξ(y)=Ξ6/D(y)​(y).

Write N(x)=(P,Q,r1a,r2b,r3c)N(x)=(P,Q,r_1a,r_2b,r_3c)N(x)=(P,Q,r1​a,r2​b,r3​c) for the normalized vector field and C(x)=C(p,s,a,b,c)\mathcal C(x)=C(p,s,a,b,c)C(x)=C(p,s,a,b,c) for its constraint.

The normalized and projective state spaces are R5\mathbb R^5R5 and R4\mathbb R^4R4 with their finite-product norms. The maps Π,Ξ,D,F\Pi,\Xi,D,FΠ,Ξ,D,F act on these state spaces.

The projective vector field in normalized time is

G(y)=6D(y)F(y).G(y)=\frac6{D(y)}F(y).G(y)=D(y)6​F(y).

Write X=R5X=\mathbb R^5X=R5, Y=R4Y=\mathbb R^4Y=R4, and let αx,t0\alpha_{x,t_0}αx,t0​​ be the maximal solution of α′=N(α)\alpha'=N(\alpha)α′=N(α) with initial value xxx at t0t_0t0​ and domain Imax⁡(x,t0)I_{\max}(x,t_0)Imax​(x,t0​). Put

Jx,t0(t)=dxαx,t0(t),Ψt0(y,t)=Π(αΞ(y),t0(t)),J_{x,t_0}(t)=d_x\alpha_{x,t_0}(t),\qquad \Psi_{t_0}(y,t)=\Pi(\alpha_{\Xi(y),t_0}(t)),Jx,t0​​(t)=dx​αx,t0​​(t),Ψt0​​(y,t)=Π(αΞ(y),t0​​(t)), Jy,t0(t)=dΠαΞ(y),t0(t)∘JΞ(y),t0(t)∘dΞy.\mathcal J_{y,t_0}(t)=d\Pi_{\alpha_{\Xi(y),t_0}(t)} \circ J_{\Xi(y),t_0}(t)\circ d\Xi_y.Jy,t0​​(t)=dΠαΞ(y),t0​​(t)​∘JΞ(y),t0​​(t)∘dΞy​.

On the invariant plane with zero active wall, use y=(m,u,w)∈R3y=(m,u,w)\in\mathbb R^3y=(m,u,w)∈R3, the inclusion ι(m,u,w)=(m,u,0,w)\iota(m,u,w)=(m,u,0,w)ι(m,u,w)=(m,u,0,w) and the restriction πs(m,u,v,w)=(m,u,w)\pi_s(m,u,v,w)=(m,u,w)πs​(m,u,v,w)=(m,u,w). Define

f(y)=(mu2−(m+1)uw+w2,1−m2−w2+u22u,−(m+1+uw−u2)w),f(y)=\left(mu^2-(m+1)uw+w^2, \frac{1-m^2-w^2+u^2}{2}u,-(m+1+uw-u^2)w\right),f(y)=(mu2−(m+1)uw+w2,21−m2−w2+u2​u,−(m+1+uw−u2)w), Δ(y)=3+m2+(u−w)2,fo(y)=6f(y)/Δ(y).\Delta(y)=3+m^2+(u-w)^2,\qquad f_o(y)=6f(y)/\Delta(y).Δ(y)=3+m2+(u−w)2,fo​(y)=6f(y)/Δ(y).

Write J=dfJ=dfJ=df, H=d2fH=d^2fH=d2f, Jo=dfoJ_o=df_oJo​=dfo​, Ho=d2foH_o=d^2f_oHo​=d2fo​ for their explicit Jacobian and Hessian expressions. The denominator derivatives are dΔ=(2m,2(u−w),2(w−u))d\Delta=(2m,2(u-w),2(w-u))dΔ=(2m,2(u−w),2(w−u)) and d2Δ(v,z)=2v0z0+2(v1−v2)(z1−z2)d^2\Delta(v,z)=2v_0z_0+2(v_1-v_2)(z_1-z_2)d2Δ(v,z)=2v0​z0​+2(v1​−v2​)(z1​−z2​). The exact original-time boundary flow and its projected variation are

Φ(y,t)=πsΨ0(ιy,t),L(y,t)=πsJιy,0(t)ι.\Phi(y,t)=\pi_s\Psi_0(\iota y,t),\qquad L(y,t)=\pi_s\mathcal J_{\iota y,0}(t)\iota.Φ(y,t)=πs​Ψ0​(ιy,t),L(y,t)=πs​Jιy,0​(t)ι.

For a fixed scale KKK, set ΦK(y,s)=Φ(y,Ks)\Phi_K(y,s)=\Phi(y,Ks)ΦK​(y,s)=Φ(y,Ks), LK(y,s)=L(y,Ks)L_K(y,s)=L(y,Ks)LK​(y,s)=L(y,Ks) and QK(y,s)=dyLK(y,s)Q_K(y,s)=d_yL_K(y,s)QK​(y,s)=dy​LK​(y,s). For directions v,zv,zv,z, define

SK(y,s;v,z)=K[Ho(ΦK)(LKv,LKz)+Jo(ΦK)QK(v,z)].S_K(y,s;v,z)=K\left[H_o(\Phi_K)(L_Kv,L_Kz)+J_o(\Phi_K)Q_K(v,z)\right].SK​(y,s;v,z)=K[Ho​(ΦK​)(LK​v,LK​z)+Jo​(ΦK​)QK​(v,z)].

All vector norms are sup norms, and operator norms are induced by them. Write b(y)=∣u∣+∣w∣b(y)=|u|+|w|b(y)=∣u∣+∣w∣, E(y)=e512b(y)E(y)=e^{512b(y)}E(y)=e512b(y), vs=(v1,v2)v_s=(v_1,v_2)vs​=(v1​,v2​), and

H(y;v,z)=b(y)∥v∥∥z∥+∥vs∥∥z∥+∥v∥∥zs∥,\mathcal H(y;v,z)=b(y)\|v\|\|z\|+\|v_s\|\|z\|+\|v\|\|z_s\|,H(y;v,z)=b(y)∥v∥∥z∥+∥vs​∥∥z∥+∥v∥∥zs​∥, B2(y)=192b(y)(40000E(y)3)+10000E(y)2.B_2(y)=192b(y)(40000E(y)^3)+10000E(y)^2.B2​(y)=192b(y)(40000E(y)3)+10000E(y)2.

The closed boundary region is Sc={(m,u,0,w):∣m−c∣≤1/10, ∣u∣,∣w∣≤1/100}\mathcal S_c=\{(m,u,0,w):|m-c|\leq1/10,\ |u|,|w|\leq1/100\}Sc​={(m,u,0,w):∣m−c∣≤1/10, ∣u∣,∣w∣≤1/100}. The boundary projection P(m,u,w)P(m,u,w)P(m,u,w) is the limit of Φ(m,u,w;t)0\Phi(m,u,w;t)_0Φ(m,u,w;t)0​ as t→∞t\to\inftyt→∞, using the total limit convention of the formal definitions.

Write πs(m,u,v,w)=(m,u,w)\pi_s(m,u,v,w)=(m,u,w)πs​(m,u,v,w)=(m,u,w), Rs(m,u,w)=P(m,u,w)−P(m,u,0)R_s(m,u,w)=P(m,u,w)-P(m,u,0)Rs​(m,u,w)=P(m,u,w)−P(m,u,0), p(y)=P(πsy)p(y)=P(\pi_s y)p(y)=P(πs​y), ρs(y)=P(y0,y1,0)−1\rho_s(y)=P(y_0,y_1,0)-1ρs​(y)=P(y0​,y1​,0)−1, ρ(y)=(p(y)−1)2+y22\rho(y)=\sqrt{(p(y)-1)^2+y_2^2}ρ(y)=(p(y)−1)2+y22​​ and χ(y)=(ρ(y),y1,y2,y3)\chi(y)=(\rho(y),y_1,y_2,y_3)χ(y)=(ρ(y),y1​,y2​,y3​). Put Es=ρ−ρsE_s=\rho-\rho_sEs​=ρ−ρs​ and Qs=2/ρ−2/ρsQ_s=2/\rho-2/\rho_sQs​=2/ρ−2/ρs​. For y=(m,u,v,w)∈R4y=(m,u,v,w)\in\mathbb R^4y=(m,u,v,w)∈R4, define

K(u,v,w)=u2+v2+w2−2(uv+vw+wu),κ=u(v+w−u),K(u,v,w)=u^2+v^2+w^2-2(uv+vw+wu),\qquad \kappa=u(v+w-u),K(u,v,w)=u2+v2+w2−2(uv+vw+wu),κ=u(v+w−u), Fm=−v2+w2+u(v−w)−mκ,Fu=12(1−m2−(v−w)2+u2)u,F_m=-v^2+w^2+u(v-w)-m\kappa,\qquad F_u=\tfrac12(1-m^2-(v-w)^2+u^2)u,Fm​=−v2+w2+u(v−w)−mκ,Fu​=21​(1−m2−(v−w)2+u2)u, Fv=(m−1−κ)v,Fw=(−m−1−κ)w,F=(Fm,Fu,Fv,Fw),F_v=(m-1-\kappa)v,\qquad F_w=(-m-1-\kappa)w,\qquad F=(F_m,F_u,F_v,F_w),Fv​=(m−1−κ)v,Fw​=(−m−1−κ)w,F=(Fm​,Fu​,Fv​,Fw​),

and E(m,v)=(m−1)2+v2E(m,v)=(m-1)^2+v^2E(m,v)=(m−1)2+v2.

Write r(m,v)=(m−1)2+v2r(m,v)=\sqrt{(m-1)^2+v^2}r(m,v)=(m−1)2+v2​ and define

g(R,v)=R2−v2,d(R,u,v,w)=g(R,v)+u(v+w−u),g(R,v)=\sqrt{R^2-v^2},\qquad d(R,u,v,w)=g(R,v)+u(v+w-u),g(R,v)=R2−v2​,d(R,u,v,w)=g(R,v)+u(v+w−u), fR=−u(v+w−u)R−g(R,v)(u−w)2/Rd,f_R=\frac{-u(v+w-u)R-g(R,v)(u-w)^2/R}{d},fR​=d−u(v+w−u)R−g(R,v)(u−w)2/R​, as=2g(R,v)−R2+2vw−w2+u22d,aw=g(R,v)−2−u(v+w−u)d.a_s=\frac{2g(R,v)-R^2+2vw-w^2+u^2}{2d},\qquad a_w=\frac{g(R,v)-2-u(v+w-u)}d.as​=2d2g(R,v)−R2+2vw−w2+u2​,aw​=dg(R,v)−2−u(v+w−u)​.

The exit clock rate is ω(y)=1−m+u(v+w−u)\omega(y)=1-m+u(v+w-u)ω(y)=1−m+u(v+w−u), and its field is Flog⁡(y)=F(y)/ω(y)F_{\log}(y)=F(y)/\omega(y)Flog​(y)=F(y)/ω(y). Division and square roots have their total real-valued Lean meanings.

For Y=(R,u,w)Y=(R,u,w)Y=(R,u,w) define Ev(Y)=(fR,asu,aww)E_v(Y)=(f_R,a_su,a_ww)Ev​(Y)=(fR​,as​u,aw​w), evaluated at (R,u,v,w)(R,u,v,w)(R,u,v,w). The coordinate maps are Q(m,u,v,w)=(r(m,v),u,w)Q(m,u,v,w)=(r(m,v),u,w)Q(m,u,v,w)=(r(m,v),u,w) and Qv−1(R,u,w)=(1−g(R,v),u,v,w)Q_v^{-1}(R,u,w)=(1-g(R,v),u,v,w)Qv−1​(R,u,w)=(1−g(R,v),u,v,w). For a curve Y(s)Y(s)Y(s) set

Mi(Y,t)=exp⁡∫0tai(Y0(s),Y1(s),e−s,Y2(s)) ds(i=s,w).M_i(Y,t)=\exp\int_0^t a_i(Y_0(s),Y_1(s),e^{-s},Y_2(s))\,ds \quad(i=s,w).Mi​(Y,t)=exp∫0t​ai​(Y0​(s),Y1​(s),e−s,Y2​(s))ds(i=s,w).

Let Rα,b,A,D\mathcal R_{\alpha,b,A,D}Rα,b,A,D​ consist of the states satisfying

7/4≤R≤17/4,∣u∣,∣w∣≤1/100,∣u∣+∣w∣≤Ab,∣R−(2+2α)∣≤Db.7/4\leq R\leq17/4,\quad |u|,|w|\leq1/100,\quad |u|+|w|\leq Ab,\quad |R-(2+2\alpha)|\leq Db.7/4≤R≤17/4,∣u∣,∣w∣≤1/100,∣u∣+∣w∣≤Ab,∣R−(2+2α)∣≤Db.

Put

B(α,b,A,D,T)=[10(α+Db)+100(Ab)2]T+100Ab+5,B(\alpha,b,A,D,T)=[10(\alpha+Db)+100(Ab)^2]T+100Ab+5,B(α,b,A,D,T)=[10(α+Db)+100(Ab)2]T+100Ab+5, L(η,α,J)=log⁡(1/η)+Jlog⁡(1/α).L(\eta,\alpha,J)=\log(1/\eta)+J\log(1/\alpha).L(η,α,J)=log(1/η)+Jlog(1/α).

Define

H0=(2−g(R,v))2+v2,d0=(u−w)2−2v(u+w),H_0=(2-g(R,v))^2+v^2,\quad d_0=(u-w)^2-2v(u+w),H0​=(2−g(R,v))2+v2,d0​=(u−w)2−2v(u+w), e=d0[1+(R2−4+d0)/H0],σv(R,u,w)=(R2+e−2)/2.e=d_0[1+(R^2-4+d_0)/H_0],\qquad \sigma_v(R,u,w)=(\sqrt{R^2+e}-2)/2.e=d0​[1+(R2−4+d0​)/H0​],σv​(R,u,w)=(R2+e​−2)/2.

Put

H0,R=2(g−2)R/g,eR=d0[2R/H0−(R2−4+d0)H0,R/H02],H_{0,R}=2(g-2)R/g,\quad e_R=d_0[2R/H_0-(R^2-4+d_0)H_{0,R}/H_0^2],H0,R​=2(g−2)R/g,eR​=d0​[2R/H0​−(R2−4+d0​)H0,R​/H02​], ewall=1+(R2−4+2d0)/H0.e_{\rm wall}=1+(R^2-4+2d_0)/H_0.ewall​=1+(R2−4+2d0​)/H0​.

The exit linearization is ξ↦(GR⋅ξ,uGs⋅ξ+asξ1,wGw⋅ξ+awξ2)\xi\mapsto(G_R\cdot\xi,uG_s\cdot\xi+a_s\xi_1, wG_w\cdot\xi+a_w\xi_2)ξ↦(GR​⋅ξ,uGs​⋅ξ+as​ξ1​,wGw​⋅ξ+aw​ξ2​); its time-dependent second-variation source is obtained by differentiating this expression along the exit clock flow. The physical-time exit clock rate is 6ω/D6\omega/D6ω/D. For the cyclic map write N=(−4(m+1),4v,4w,4u)N=(-4(m+1),4v,4w,4u)N=(−4(m+1),4v,4w,4u), so A=(1,0,0,0)+N/H\mathcal A=(1,0,0,0)+N/HA=(1,0,0,0)+N/H; its derivative and Hessian are the quotient-rule expressions using N′N'N′, dHdHdH, and d2Hd^2Hd2H. The natural height is h∗(Y)=η((Y0−2)/2)Jh_*(Y)=\eta((Y_0-2)/2)^Jh∗​(Y)=η((Y0​−2)/2)J. The natural state and phase are P(h∗(Y),Y)\mathcal P(h_*(Y),Y)P(h∗​(Y),Y) and Fh∗(Y)(Y)F_{h_*(Y)}(Y)Fh∗​(Y)​(Y). These definitions also include λ−(I)=(I−I2−4)/2\lambda_-(I)=(I-\sqrt{I^2-4})/2λ−​(I)=(I−I2−4​)/2.

Definition code
import Mathlib.Analysis.Calculus.ImplicitContDiff
import Mathlib.Tactic.NormNum
import Mathlib.Analysis.Calculus.Deriv.Pow
import Mathlib.Analysis.Calculus.Deriv.Mul
import Mathlib.Analysis.Calculus.Deriv.Add
import Mathlib.Analysis.Real.Sqrt
import Mathlib.Tactic.Linarith
import Mathlib.Tactic.Ring
import Mathlib.Analysis.Calculus.Deriv.Inv
import Mathlib.Tactic.FieldSimp
import Mathlib.Tactic.LinearCombination
import Mathlib.Analysis.SpecialFunctions.Sqrt
import Mathlib.Analysis.SpecialFunctions.Pow.Deriv
import Definitions.Def_BianchiIX_TerminalConnectionData
import Mathlib.Analysis.SpecialFunctions.Arsinh
import Definitions.Def_BianchiIX_Algebra_Projective
import Mathlib.Analysis.Calculus.ContDiff.Operations
import Mathlib.Analysis.Calculus.Deriv.Prod
import Mathlib.Analysis.Calculus.MeanValue
import Mathlib.Data.Fin.VecNotation
import Mathlib.Tactic.FinCases
import Mathlib.Tactic.FunProp
import Mathlib.Analysis.SpecialFunctions.Arcosh
import Mathlib.MeasureTheory.Integral.IntervalIntegral.FundThmCalculus
import Mathlib.Analysis.Calculus.InverseFunctionTheorem.Deriv
import Mathlib.Analysis.Calculus.Deriv.MeanValue
import Mathlib.Analysis.SpecialFunctions.Integrals.Basic
import Mathlib.Analysis.ODE.ExistUnique
import Mathlib.Analysis.Calculus.FDeriv.Partial
import Mathlib.Analysis.SpecialFunctions.Log.Deriv
import Definitions.Def_BianchiIX_Dynamics_ProjectiveState
import Mathlib.Analysis.Normed.Module.FiniteDimension
import Mathlib.LinearAlgebra.Determinant
import Mathlib.Topology.Instances.Matrix
import Mathlib.LinearAlgebra.Matrix.Trace
import Mathlib.Analysis.Calculus.FDeriv.Analytic
import Mathlib.Topology.Algebra.Module.Determinant
import Mathlib.Analysis.Calculus.Deriv.Pi
import Mathlib.Data.Real.Basic
import Mathlib.Analysis.Calculus.ContDiff.Comp
import Mathlib.Analysis.Calculus.Deriv.Comp
import Mathlib.Analysis.ODE.Gronwall
import Mathlib.Data.Matrix.Mul
import Mathlib.Topology.Order.ProjIcc
import Mathlib.LinearAlgebra.Matrix.ToLin
import Mathlib.Topology.Order.Compact
import Mathlib.MeasureTheory.Integral.IntervalIntegral.DistLEIntegral
import Definitions.Def_BianchiIX_Dynamics_MaximalFlow
import Definitions.Def_BianchiIX_Dynamics_ProjectiveFlowVariation
import Definitions.Def_BianchiIX_Dynamics_CoreTail
import Mathlib.Analysis.Calculus.FDeriv.Symmetric
import Mathlib.Tactic.Module
import Mathlib.Topology.Order.IntermediateValue
import Mathlib.Analysis.SpecialFunctions.Log.Basic
import Mathlib.Tactic.Positivity
import Mathlib.Topology.Instances.Real.Lemmas
import Mathlib.Tactic.GCongr
import Mathlib.Analysis.Calculus.Deriv.Slope
import Mathlib.Topology.Order.OrderClosed
import Mathlib.Algebra.BigOperators.Fin
import Mathlib.Analysis.SpecialFunctions.Trigonometric.ArctanDeriv
import Mathlib.Analysis.SpecialFunctions.Trigonometric.DerivHyp
import Mathlib.Topology.MetricSpace.Thickening
import Mathlib.Analysis.Calculus.ParametricIntervalIntegral
import Mathlib.MeasureTheory.Integral.DominatedConvergence
import Definitions.Def_BianchiIX_EntrancePhaseResponseData
import Mathlib.Topology.MetricSpace.Cauchy
import Mathlib.Analysis.Calculus.UniformLimitsDeriv
import Definitions.Def_BianchiIX_Dynamics_CoreTailMultipliers
import Definitions.Def_BianchiIX_Algebra_CoreAngle
import Mathlib.Analysis.Calculus.InverseFunctionTheorem.ContDiff
import Definitions.Def_BianchiIX_CompactPhaseExpansionData
import Mathlib.Analysis.SpecialFunctions.ImproperIntegrals
import Mathlib.Analysis.Calculus.ParametricIntegral
import Mathlib.MeasureTheory.Integral.Bochner.ContinuousLinearMap
import Mathlib.Analysis.Calculus.FDeriv.Mul
import Mathlib.Analysis.SpecialFunctions.ExpDeriv
import Mathlib.MeasureTheory.Integral.IntervalIntegral.IntegrationByParts
import Mathlib.Analysis.SpecialFunctions.Pow.Real
import Mathlib.Analysis.Normed.Operator.Basic
import Mathlib.LinearAlgebra.Matrix.Notation
import Mathlib.Analysis.Calculus.ContDiff.Deriv
import Mathlib.Analysis.SpecialFunctions.Pow.Continuity

section

noncomputable section

namespace BianchiIX

open Set Filter
open scoped Topology

def intervalScalarInverse {E : Type*} (Q : E × ℝ → ℝ) (I : Set ℝ) (p : E × ℝ) : ℝ := by
  classical
  exact if h : ∃ t ∈ I, Q (p.1, t) = p.2 then Classical.choose h else 0









end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def singleStrongInvariant (m u : ℝ) : ℝ := (1 + m ^ 2 + u ^ 2) / m

def singleStrongUpperRoot (I : ℝ) : ℝ := (I + Real.sqrt (I ^ 2 - 4)) / 2

def singleStrongLowerRoot (I : ℝ) : ℝ := (I - Real.sqrt (I ^ 2 - 4)) / 2

def singleStrongPhase (m u : ℝ) : ℝ := 2 / (singleStrongUpperRoot (singleStrongInvariant m u) - 1)





















end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def exitSingleBase (R v : ℝ) : ℝ := (2 - exitTailGap R v) ^ 2 + v ^ 2

def exitSingleDefect (u v w : ℝ) : ℝ := (u - w) ^ 2 - 2 * v * (u + w)

def exitSingleError (R u v w : ℝ) : ℝ :=
  exitSingleDefect u v w * (1 + (R ^ 2 - 4 + exitSingleDefect u v w) / exitSingleBase R v)

def exitSinglePhase (R u v w : ℝ) : ℝ := (Real.sqrt (R ^ 2 + exitSingleError R u v w) - 2) / 2















end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def exitTailRadialGradient (R u v w : ℝ) : Fin 3 → ℝ :=
  let t := exitTailGap R v
  let d := exitTailSpeed R u v w
  let g := exitTailRadial R u v w
  ![(-coupling u v w - (1 / t - t / R ^ 2) * (u - w) ^ 2 - g * R / t) / d,
    (-(v + w - 2 * u) * R - 2 * t * (u - w) / R - g * (v + w - 2 * u)) / d,
    (-u * R + 2 * t * (u - w) / R - g * u) / d]

def exitTailStrongRateGradient (R u v w : ℝ) : Fin 3 → ℝ :=
  let t := exitTailGap R v
  let d := exitTailSpeed R u v w
  let a := exitTailStrongRate R u v w
  ![(R / t - R - a * R / t) / d,
    (u - a * (v + w - 2 * u)) / d, (v - w - a * u) / d]

def exitTailWeakRateGradient (R u v w : ℝ) : Fin 3 → ℝ :=
  let t := exitTailGap R v
  let d := exitTailSpeed R u v w
  let a := exitTailWeakRate R u v w
  ![((1 - a) * R / t) / d,
    (-(1 + a) * (v + w - 2 * u)) / d, (-(1 + a) * u) / d]









def exitTailLinearization (v : ℝ) (Y ξ : ExitTailState) : ExitTailState :=
  ![∑ i : Fin 3, exitTailRadialGradient (Y 0) (Y 1) v (Y 2) i * ξ i,
    Y 1 * (∑ i : Fin 3, exitTailStrongRateGradient (Y 0) (Y 1) v (Y 2) i * ξ i) +
      exitTailStrongRate (Y 0) (Y 1) v (Y 2) * ξ 1,
    Y 2 * (∑ i : Fin 3, exitTailWeakRateGradient (Y 0) (Y 1) v (Y 2) i * ξ i) +
      exitTailWeakRate (Y 0) (Y 1) v (Y 2) * ξ 2]





end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def exitLogClock (y : ProjectiveState) : ℝ := -Real.log (y 2)

def exitOriginalClockRate (y : ProjectiveState) : ℝ :=
  (6 / projectiveStateDenom y) * exitClockRate y















end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def exitSingleBaseRadial (R v : ℝ) : ℝ := 2 * (exitTailGap R v - 2) * R / exitTailGap R v

def exitSingleErrorRadial (R u v w : ℝ) : ℝ :=
  exitSingleDefect u v w * (2 * R / exitSingleBase R v -
    (R ^ 2 - 4 + exitSingleDefect u v w) * exitSingleBaseRadial R v / (exitSingleBase R v) ^ 2)

def exitSingleErrorWallFactor (R u v w : ℝ) : ℝ :=
  1 + (R ^ 2 - 4 + 2 * exitSingleDefect u v w) / exitSingleBase R v



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open scoped Topology

def cycleNumerator (i : Fin 4) (y : ProjectiveState) : ℝ :=
  ![-4 * (y 0 + 1), 4 * y 2, 4 * y 3, 4 * y 1] i

def cycleNumeratorSlope (i : Fin 4) (v : ProjectiveState) : ℝ :=
  ![-4 * v 0, 4 * v 2, 4 * v 3, 4 * v 1] i

def cycleDenomDifferential (y v : ProjectiveState) : ℝ :=
  2 * (y 0 + 1) * v 0 + 2 * (y 1 - y 2 - y 3) * v 1 +
    2 * (y 2 - y 1 - y 3) * v 2 + 2 * (y 3 - y 1 - y 2) * v 3

def cycleDenomHessian (v w : ProjectiveState) : ℝ :=
  2 * (v 0 * w 0 + (v 1 - v 2 - v 3) * w 1 +
    (v 2 - v 1 - v 3) * w 2 + (v 3 - v 1 - v 2) * w 3)

def cycleDifferential (i : Fin 4) (y v : ProjectiveState) : ℝ :=
  cycleNumeratorSlope i v / projectiveCycleDenom y -
    cycleNumerator i y * cycleDenomDifferential y v / projectiveCycleDenom y ^ 2

def cycleHessian (i : Fin 4) (y v w : ProjectiveState) : ℝ :=
  -(cycleNumeratorSlope i v * cycleDenomDifferential y w +
    cycleNumeratorSlope i w * cycleDenomDifferential y v +
    cycleNumerator i y * cycleDenomHessian v w) / projectiveCycleDenom y ^ 2 +
    2 * cycleNumerator i y * cycleDenomDifferential y v * cycleDenomDifferential y w /
      projectiveCycleDenom y ^ 3



















end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def exitTailRateGradient (v : ℝ) (Y : ExitTailState) (i : Fin 2) : ExitTailState :=
  if i = 0 then exitTailStrongRateGradient (Y 0) (Y 1) v (Y 2)
  else exitTailWeakRateGradient (Y 0) (Y 1) v (Y 2)









end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set MeasureTheory

def exitTailVariationBudget (α b A D T : ℝ) : ℝ :=
  2 * exitTailRateBudget α b A D T + 2048 * (1 + A * b * T)







end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set Filter
open scoped Topology

def exitEntryFlow (p : ExitTailState × ℝ) : ProjectiveState :=
  projectiveOriginalFlow (coreTailEntry p.1) 0 p.2

def exitEntryClock (p : ExitTailState × ℝ) : ℝ := exitLogClock (exitEntryFlow p)

def exitEntryRegular (Y : ExitTailState) (t : ℝ) : Prop :=
  1 < Y 0 ∧ 0 < projectiveStateDenom (coreTailEntry Y) ∧
    t ∈ normalizedMaxDomain (fromProjectiveState (coreTailEntry Y)) 0 ∧
    0 < 1 + normalizedMaxCurve (fromProjectiveState (coreTailEntry Y)) 0 t 0 ∧
    0 < exitEntryFlow (Y, t) 2 ∧ exitEntryFlow (Y, t) 0 < 1 ∧
    0 < exitClockRate (exitEntryFlow (Y, t))















end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set Filter
open scoped Topology

def exitClockTime (I : Set ℝ) : ExitTailState × ℝ → ℝ := intervalScalarInverse exitEntryClock I

def exitClockFlow (I : Set ℝ) (p : ExitTailState × ℝ) : ExitTailState :=
  toExitTailState (exitEntryFlow (p.1, exitClockTime I p))



















end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set Filter
open scoped Topology

def exitClockVariationRate (I : Set ℝ) (w : ExitTailState) (p : ExitTailState × ℝ) : ExitTailState :=
  fderiv ℝ (exitTailField (Real.exp (-p.2))) (exitClockFlow I p)
    (fderiv ℝ (exitClockFlow I) p (w, 0))







end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set Filter MeasureTheory
open scoped Topology

def exitClockRateIntegrand (I : Set ℝ) (i : Fin 2) (p : ExitTailState × ℝ) : ℝ :=
  exitTailWallRate (Real.exp (-p.2)) (exitClockFlow I p) i

def exitClockLogMultiplier (I : Set ℝ) (i : Fin 2) (Y : ExitTailState) (s : ℝ) : ℝ :=
  ∫ r in (0 : ℝ)..s, exitClockRateIntegrand I i (Y, r)











end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def terminalCyclicInput (h : ℝ) (Y : ExitTailState) : ProjectiveState :=
  projectiveCycle (coreTailState h Y)



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Filter
open scoped Topology

def terminalSinglePhase (h : ℝ) (Y : ExitTailState) : ℝ :=
  singleStrongPhase (terminalCyclicInput h Y 0) (terminalCyclicInput h Y 1)





end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def terminalEntranceInput (h : ℝ) (Y : ExitTailState) : ProjectiveState :=
  entranceRadialCoordinates (terminalCyclicInput h Y)





end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Filter
open scoped Topology

/-- Indices 0 and 1 transport the old U and W walls to the new w and v walls, respectively. -/
def terminalWallFactor (h : ℝ) (Y : ExitTailState) (i : Fin 2) : ℝ :=
  4 * coreTailMultiplier h Y i / projectiveCycleDenom (coreTailState h Y)











end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Filter
open scoped Topology

def terminalConnectionPhase (h : ℝ) (Y : ExitTailState) : ℝ :=
  coreRadius (terminalConnectionState h Y 0) (terminalConnectionState h Y 2)

def terminalConnectionError (h : ℝ) (Y : ExitTailState) : ℝ :=
  terminalConnectionPhase h Y - terminalEntranceInput h Y 0







end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def terminalReciprocalError (h : ℝ) (Y : ExitTailState) : ℝ :=
  2 / terminalConnectionPhase h Y - 2 / terminalEntranceInput h Y 0



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set Filter
open scoped Topology

def stableStrongBase : StableState →L[ℝ] StableState :=
  ContinuousLinearMap.pi ![ContinuousLinearMap.proj 0, ContinuousLinearMap.proj 1, 0]





def stableStrongProjectionError (y : StableState) : ℝ :=
  stableBoundaryProjection (y 0) (y 1) (y 2) - stableBoundaryProjection (y 0) (y 1) 0











end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def entranceSinglePhaseError (y : ProjectiveState) : ℝ :=
  entranceStableRadius y - (stableBoundaryProjection (y 0) (y 1) 0 - 1)











end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Filter
open scoped Topology

def entranceSingleRadius (y : ProjectiveState) : ℝ := stableBoundaryProjection (y 0) (y 1) 0 - 1

def entranceSingleReciprocalError (y : ProjectiveState) : ℝ :=
  2 / entranceStableRadius y - 2 / entranceSingleRadius y







end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Filter
open scoped Topology

def terminalSingleReciprocalError (h : ℝ) (Y : ExitTailState) : ℝ :=
  2 / terminalEntranceInput h Y 0 -
    singleStrongPhase (terminalCyclicInput h Y 0) (terminalCyclicInput h Y 1)





end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Filter Set
open scoped Topology

def terminalConnectionLogWalls (h : ℝ) (Y : ExitTailState) : Fin 2 → ℝ :=
  entranceLogWalls 1 (terminalConnectionState h Y)











end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def terminalBoundedPhaseError (h : ℝ) (Y : ExitTailState) : ℝ :=
  2 / terminalConnectionPhase h Y -
    singleStrongPhase (terminalCyclicInput h Y 0) (terminalCyclicInput h Y 1)







end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def terminalBoundedPhase (h : ℝ) (Y : ExitTailState) : ℝ := 2 / terminalConnectionPhase h Y









end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Filter
open scoped Topology

def terminalNaturalPhase (η : ℝ) (J : ℕ) (Y : ExitTailState) : ℝ :=
  terminalConnectionPhase (η * ((Y 0 - 2) / 2) ^ J) Y



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Filter
open scoped Topology

def terminalNaturalState (η : ℝ) (J : ℕ) (Y : ExitTailState) : ProjectiveState :=
  terminalConnectionState (η * ((Y 0 - 2) / 2) ^ J) Y



end BianchiIX

end
end
Source
https://github.com/l46430640-del/ERA-STATISTICS-FORMAL/tree/d74c241819f896ab8d050b21ace80c8e90eabd9e/BianchiIX

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me