Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Partial-core returns, complete era maps and successor curves

Definition
BianchiIX_EraReturnData

by trppppp · Nov 6, 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.

Define the core radius, angle, and clock factor by

r(m,v)=(m−1)2+v2,θ(m,v)=arsinh⁡1−mv,d(r,u,w)=1−(u−w)2r2.r(m,v)=\sqrt{(m-1)^2+v^2},\qquad \theta(m,v)=\operatorname{arsinh}\frac{1-m}{v},\qquad d(r,u,w)=1-\frac{(u-w)^2}{r^2}.r(m,v)=(m−1)2+v2​,θ(m,v)=arsinhv1−m​,d(r,u,w)=1−r2(u−w)2​.

For y=(m,u,v,w)y=(m,u,v,w)y=(m,u,v,w), write ry=r(m,v)r_y=r(m,v)ry​=r(m,v), θy=θ(m,v)\theta_y=\theta(m,v)θy​=θ(m,v), and dy=d(ry,u,w)d_y=d(r_y,u,w)dy​=d(ry​,u,w). All quotients use the total convention x/0=0x/0=0x/0=0.

The core state space is R3\mathbb R^3R3 with coordinates Y=(r,u,p)Y=(r,u,p)Y=(r,u,p). Define

η(m,u,v,w)=(r(m,v),u,vw),ιz(r,u,p)=(1−rtanh⁡z,u,rcosh⁡z,pr/cosh⁡z).\eta(m,u,v,w)=(r(m,v),u,vw),\qquad \iota_z(r,u,p)=\left(1-r\tanh z,u,\frac r{\cosh z}, \frac p{r/\cosh z}\right).η(m,u,v,w)=(r(m,v),u,vw),ιz​(r,u,p)=(1−rtanhz,u,coshzr​,r/coshzp​).

For fixed (z,Y)(z,Y)(z,Y) put v=r/cosh⁡zv=r/\cosh zv=r/coshz, w=p/vw=p/vw=p/v, d=d(r,u,w)d=d(r,u,w)d=d(r,u,w), and define

V(z,Y)=(−κ(u,v,w)−tanh⁡z (u−w)2/rd,u[−r/2+tanh⁡z+w/cosh⁡z+(u2−w2)/(2r)]d,−2p(1+κ(u,v,w))rd).\mathcal V(z,Y)=\left( \frac{-\kappa(u,v,w)-\tanh z\,(u-w)^2/r}{d},\quad \frac{u[-r/2+\tanh z+w/\cosh z+(u^2-w^2)/(2r)]}{d},\quad \frac{-2p(1+\kappa(u,v,w))}{rd}\right).V(z,Y)=(d−κ(u,v,w)−tanhz(u−w)2/r​,du[−r/2+tanhz+w/coshz+(u2−w2)/(2r)]​,rd−2p(1+κ(u,v,w))​).

For real kkk, put εk=1/(10000k)\varepsilon_k=1/(10000k)εk​=1/(10000k) and define

Kk={(m,u,v,w):∣E(m,v)−k2∣≤1/8,v≥1,0≤u≤εk,0≤vw≤εk}.\mathcal K_k=\{(m,u,v,w): |E(m,v)-k^2|\leq1/8,\quad v\geq1, \quad0\leq u\leq\varepsilon_k,\quad0\leq vw\leq\varepsilon_k\}.Kk​={(m,u,v,w):∣E(m,v)−k2∣≤1/8,v≥1,0≤u≤εk​,0≤vw≤εk​}.

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​.

Write G(y)=6F(y)/D(y)G(y)=6F(y)/D(y)G(y)=6F(y)/D(y), λ(y)=6rydy/D(y)\lambda(y)=6r_yd_y/D(y)λ(y)=6ry​dy​/D(y), and Sk(R,U,P)=(R,kU,P)S_k(R,U,P)=(R,kU,P)Sk​(R,U,P)=(R,kU,P). Define

Pyξ=dηyξ−dθyξ V(θy,η(y)),Wy,t0,w(t)=PΨt0(y,t)(Jy,t0(t)w).\mathcal P_y\xi=d\eta_y\xi-d\theta_y\xi\,\mathcal V(\theta_y,\eta(y)), \qquad \mathcal W_{y,t_0,w}(t)= \mathcal P_{\Psi_{t_0}(y,t)}(\mathcal J_{y,t_0}(t)w).Py​ξ=dηy​ξ−dθy​ξV(θy​,η(y)),Wy,t0​,w​(t)=PΨt0​​(y,t)​(Jy,t0​​(t)w).

The entrance parametrization and radial exit factor are

e(k,u,w)=(1+k2−1,u,1,w),e(k,u,w)=(1+\sqrt{k^2-1},u,1,w),e(k,u,w)=(1+k2−1​,u,1,w), fy=V(θy,η(y))0,cy=−tanh⁡θy,hy=cy+fy/ry,ρy=cy/hy.f_y=\mathcal V(\theta_y,\eta(y))_0,\qquad c_y=-\tanh\theta_y, \qquad h_y=c_y+f_y/r_y,\qquad \rho_y=c_y/h_y.fy​=V(θy​,η(y))0​,cy​=−tanhθy​,hy​=cy​+fy​/ry​,ρy​=cy​/hy​.

For projective initial data yyy, put

Ay,t0(t)=1+(αΞ(y),t0(t))0.A_{y,t_0}(t)=1+(\alpha_{\Xi(y),t_0}(t))_0.Ay,t0​​(t)=1+(αΞ(y),t0​​(t))0​.

Coordinates of three-vectors are indexed by 0,1,20,1,20,1,2; projective coordinates are indexed by m,u,v,wm,u,v,wm,u,v,w. 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.

Define the core radius, angle, and clock factor by

r(m,v)=(m−1)2+v2,θ(m,v)=arsinh⁡1−mv,d(r,u,w)=1−(u−w)2r2.r(m,v)=\sqrt{(m-1)^2+v^2},\qquad \theta(m,v)=\operatorname{arsinh}\frac{1-m}{v},\qquad d(r,u,w)=1-\frac{(u-w)^2}{r^2}.r(m,v)=(m−1)2+v2​,θ(m,v)=arsinhv1−m​,d(r,u,w)=1−r2(u−w)2​.

For y=(m,u,v,w)y=(m,u,v,w)y=(m,u,v,w), write ry=r(m,v)r_y=r(m,v)ry​=r(m,v), θy=θ(m,v)\theta_y=\theta(m,v)θy​=θ(m,v), and dy=d(ry,u,w)d_y=d(r_y,u,w)dy​=d(ry​,u,w). All quotients use the total convention x/0=0x/0=0x/0=0.

The core state space is R3\mathbb R^3R3 with coordinates Y=(r,u,p)Y=(r,u,p)Y=(r,u,p). Define

η(m,u,v,w)=(r(m,v),u,vw),ιz(r,u,p)=(1−rtanh⁡z,u,rcosh⁡z,pr/cosh⁡z).\eta(m,u,v,w)=(r(m,v),u,vw),\qquad \iota_z(r,u,p)=\left(1-r\tanh z,u,\frac r{\cosh z}, \frac p{r/\cosh z}\right).η(m,u,v,w)=(r(m,v),u,vw),ιz​(r,u,p)=(1−rtanhz,u,coshzr​,r/coshzp​).

For fixed (z,Y)(z,Y)(z,Y) put v=r/cosh⁡zv=r/\cosh zv=r/coshz, w=p/vw=p/vw=p/v, d=d(r,u,w)d=d(r,u,w)d=d(r,u,w), and define

V(z,Y)=(−κ(u,v,w)−tanh⁡z (u−w)2/rd,u[−r/2+tanh⁡z+w/cosh⁡z+(u2−w2)/(2r)]d,−2p(1+κ(u,v,w))rd).\mathcal V(z,Y)=\left( \frac{-\kappa(u,v,w)-\tanh z\,(u-w)^2/r}{d},\quad \frac{u[-r/2+\tanh z+w/\cosh z+(u^2-w^2)/(2r)]}{d},\quad \frac{-2p(1+\kappa(u,v,w))}{rd}\right).V(z,Y)=(d−κ(u,v,w)−tanhz(u−w)2/r​,du[−r/2+tanhz+w/coshz+(u2−w2)/(2r)]​,rd−2p(1+κ(u,v,w))​).

At (z,Y)(z,Y)(z,Y) use the abbreviations

s=1/cosh⁡z,c=−tanh⁡z,δ=u−w,s=1/\cosh z,\quad c=-\tanh z,\quad \delta=u-w,s=1/coshz,c=−tanhz,δ=u−w, C=u(v+w−u),a=−r/2+tanh⁡z+w/cosh⁡z+(u2−w2)/(2r),g=−C+cδ2/r.C=u(v+w-u),\quad a=-r/2+\tanh z+w/\cosh z+(u^2-w^2)/(2r), \quad g=-C+c\delta^2/r.C=u(v+w−u),a=−r/2+tanhz+w/coshz+(u2−w2)/(2r),g=−C+cδ2/r.

Let Q(z,Y)=(C,d,a,g)\mathcal Q(z,Y)=(C,d,a,g)Q(z,Y)=(C,d,a,g). Define the 4×34\times34×3 matrix H(z,Y)H(z,Y)H(z,Y) by its rows, in the order (C,d,a,g)(C,d,a,g)(C,d,a,g) and columns (r,u,p)(r,u,p)(r,u,p):

HC=(u(s−w/r), v−2u+w, u/v),H_C=(u(s-w/r),\ v-2u+w,\ u/v),HC​=(u(s−w/r), v−2u+w, u/v), Hd=(2δ(u−2w)/r3, −2δ/r2, 2δ/(vr2)),H_d=(2\delta(u-2w)/r^3,\ -2\delta/r^2,\ 2\delta/(vr^2)),Hd​=(2δ(u−2w)/r3, −2δ/r2, 2δ/(vr2)), Ha=(−1/2−sw/r+(3w2−u2)/(2r2), u/r, (1−w/v)/r),H_a=(-1/2-sw/r+(3w^2-u^2)/(2r^2),\ u/r,\ (1-w/v)/r),Ha​=(−1/2−sw/r+(3w2−u2)/(2r2), u/r, (1−w/v)/r), Hg=(−u(s−w/r)+cδ(3w−u)/r2, −(v−2u+w)+2cδ/r, −u/v−2cδ/(rv)).H_g=(-u(s-w/r)+c\delta(3w-u)/r^2,\ -(v-2u+w)+2c\delta/r,\ -u/v-2c\delta/(rv)).Hg​=(−u(s−w/r)+cδ(3w−u)/r2, −(v−2u+w)+2cδ/r, −u/v−2cδ/(rv)).

Define the two rates and the 2×32\times32×3 matrix JJJ by

L0(z,Y)=a/d,L1(z,Y)=−2(1+C)/(rd),L_0(z,Y)=a/d,\qquad L_1(z,Y)=-2(1+C)/(rd),L0​(z,Y)=a/d,L1​(z,Y)=−2(1+C)/(rd), J0j=Haj/d−aHdj/d2,J_{0j}=H_{aj}/d-aH_{dj}/d^2,J0j​=Haj​/d−aHdj​/d2, J1j=δj02(1+C)r2d−2HCjrd+2(1+C)Hdjrd2,j=0,1,2.J_{1j}=\delta_{j0}\frac{2(1+C)}{r^2d} -\frac{2H_{Cj}}{rd}+\frac{2(1+C)H_{dj}}{rd^2},\qquad j=0,1,2.J1j​=δj0​r2d2(1+C)​−rd2HCj​​+rd22(1+C)Hdj​​,j=0,1,2.

Here δj0\delta_{j0}δj0​ is the Kronecker symbol. The total formulas are defined also when a wall amplitude vanishes, and use x/0=0x/0=0x/0=0.

Define the rates in normalized time and the positive flow factors by

R(y)=λ(y)L(θy,η(y)),Mi(y,t0,a,b)=exp⁡(∫abRi(Ψt0(y,t)) dt),i=0,1.R(y)=\lambda(y)L(\theta_y,\eta(y)),\qquad M_i(y,t_0,a,b)=\exp\left(\int_a^b R_i(\Psi_{t_0}(y,t))\,dt\right), \quad i=0,1.R(y)=λ(y)L(θy​,η(y)),Mi​(y,t0​,a,b)=exp(∫ab​Ri​(Ψt0​​(y,t))dt),i=0,1.

The interval integral is oriented and uses its total Bochner-integral definition for every choice of endpoints and initial data.

For real k,zk,zk,z, define

Γk(z)=2arctan⁡(sinh⁡z)+2410000kz,C∗=exp⁡(2π+4810000).\Gamma_k(z)=2\arctan(\sinh z)+\frac{24}{10000k}z, \qquad C_* = \exp\left(2\pi+\frac{48}{10000}\right).Γk​(z)=2arctan(sinhz)+10000k24​z,C∗​=exp(2π+1000048​).

Division at zero has its total real-field meaning.

Write r(y)=r(y0,y2)\mathfrak r(y)=r(y_0,y_2)r(y)=r(y0​,y2​). For a time function τ:R4→R\tau:\mathbb R^4\to\mathbb Rτ:R4→R put miτ(y)=Mi(y,0,0,τ(y))m_i^\tau(y)=M_i(y,0,0,\tau(y))miτ​(y)=Mi​(y,0,0,τ(y)). Let τ↓(y)\tau_\downarrow(y)τ↓​(y) be the infimum of positive times in the normalized maximal domain at which (Ψ0(y,t))v=1(\Psi_0(y,t))_v=1(Ψ0​(y,t))v​=1 with negative event rate, and set F↓(y)=η(Ψ0(y,τ↓(y)))F_\downarrow(y)=\eta(\Psi_0(y,\tau_\downarrow(y)))F↓​(y)=η(Ψ0​(y,τ↓​(y))). The incoming condition means yv>1y_v>1yv​>1, or yv=1y_v=1yv​=1 and ym≥1y_m\geq1ym​≥1. Define

Qk(u,p;v,dz)=∣v0∣+k∣v1∣+∣v2∣+(k2u+p/k)∣dz∣.Q_k(u,p;v,dz)=|v_0|+k|v_1|+|v_2|+(k^2u+p/k)|dz|.Qk​(u,p;v,dz)=∣v0​∣+k∣v1​∣+∣v2​∣+(k2u+p/k)∣dz∣.

All vector norms are supremum norms and operator norms are induced by them. 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.

Define the core radius, angle, and clock factor by

r(m,v)=(m−1)2+v2,θ(m,v)=arsinh⁡1−mv,d(r,u,w)=1−(u−w)2r2.r(m,v)=\sqrt{(m-1)^2+v^2},\qquad \theta(m,v)=\operatorname{arsinh}\frac{1-m}{v},\qquad d(r,u,w)=1-\frac{(u-w)^2}{r^2}.r(m,v)=(m−1)2+v2​,θ(m,v)=arsinhv1−m​,d(r,u,w)=1−r2(u−w)2​.

For y=(m,u,v,w)y=(m,u,v,w)y=(m,u,v,w), write ry=r(m,v)r_y=r(m,v)ry​=r(m,v), θy=θ(m,v)\theta_y=\theta(m,v)θy​=θ(m,v), and dy=d(ry,u,w)d_y=d(r_y,u,w)dy​=d(ry​,u,w). All quotients use the total convention x/0=0x/0=0x/0=0.

The core state space is R3\mathbb R^3R3 with coordinates Y=(r,u,p)Y=(r,u,p)Y=(r,u,p). Define

η(m,u,v,w)=(r(m,v),u,vw),ιz(r,u,p)=(1−rtanh⁡z,u,rcosh⁡z,pr/cosh⁡z).\eta(m,u,v,w)=(r(m,v),u,vw),\qquad \iota_z(r,u,p)=\left(1-r\tanh z,u,\frac r{\cosh z}, \frac p{r/\cosh z}\right).η(m,u,v,w)=(r(m,v),u,vw),ιz​(r,u,p)=(1−rtanhz,u,coshzr​,r/coshzp​).

For fixed (z,Y)(z,Y)(z,Y) put v=r/cosh⁡zv=r/\cosh zv=r/coshz, w=p/vw=p/vw=p/v, d=d(r,u,w)d=d(r,u,w)d=d(r,u,w), and define

V(z,Y)=(−κ(u,v,w)−tanh⁡z (u−w)2/rd,u[−r/2+tanh⁡z+w/cosh⁡z+(u2−w2)/(2r)]d,−2p(1+κ(u,v,w))rd).\mathcal V(z,Y)=\left( \frac{-\kappa(u,v,w)-\tanh z\,(u-w)^2/r}{d},\quad \frac{u[-r/2+\tanh z+w/\cosh z+(u^2-w^2)/(2r)]}{d},\quad \frac{-2p(1+\kappa(u,v,w))}{rd}\right).V(z,Y)=(d−κ(u,v,w)−tanhz(u−w)2/r​,du[−r/2+tanhz+w/coshz+(u2−w2)/(2r)]​,rd−2p(1+κ(u,v,w))​).

For real kkk, put εk=1/(10000k)\varepsilon_k=1/(10000k)εk​=1/(10000k) and define

Kk={(m,u,v,w):∣E(m,v)−k2∣≤1/8,v≥1,0≤u≤εk,0≤vw≤εk}.\mathcal K_k=\{(m,u,v,w): |E(m,v)-k^2|\leq1/8,\quad v\geq1, \quad0\leq u\leq\varepsilon_k,\quad0\leq vw\leq\varepsilon_k\}.Kk​={(m,u,v,w):∣E(m,v)−k2∣≤1/8,v≥1,0≤u≤εk​,0≤vw≤εk​}.

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​.

Write G(y)=6F(y)/D(y)G(y)=6F(y)/D(y)G(y)=6F(y)/D(y), λ(y)=6rydy/D(y)\lambda(y)=6r_yd_y/D(y)λ(y)=6ry​dy​/D(y), and Sk(R,U,P)=(R,kU,P)S_k(R,U,P)=(R,kU,P)Sk​(R,U,P)=(R,kU,P). Define

Pyξ=dηyξ−dθyξ V(θy,η(y)),Wy,t0,w(t)=PΨt0(y,t)(Jy,t0(t)w).\mathcal P_y\xi=d\eta_y\xi-d\theta_y\xi\,\mathcal V(\theta_y,\eta(y)), \qquad \mathcal W_{y,t_0,w}(t)= \mathcal P_{\Psi_{t_0}(y,t)}(\mathcal J_{y,t_0}(t)w).Py​ξ=dηy​ξ−dθy​ξV(θy​,η(y)),Wy,t0​,w​(t)=PΨt0​​(y,t)​(Jy,t0​​(t)w).

The entrance parametrization and radial exit factor are

e(k,u,w)=(1+k2−1,u,1,w),e(k,u,w)=(1+\sqrt{k^2-1},u,1,w),e(k,u,w)=(1+k2−1​,u,1,w), fy=V(θy,η(y))0,cy=−tanh⁡θy,hy=cy+fy/ry,ρy=cy/hy.f_y=\mathcal V(\theta_y,\eta(y))_0,\qquad c_y=-\tanh\theta_y, \qquad h_y=c_y+f_y/r_y,\qquad \rho_y=c_y/h_y.fy​=V(θy​,η(y))0​,cy​=−tanhθy​,hy​=cy​+fy​/ry​,ρy​=cy​/hy​.

For projective initial data yyy, put

Ay,t0(t)=1+(αΞ(y),t0(t))0.A_{y,t_0}(t)=1+(\alpha_{\Xi(y),t_0}(t))_0.Ay,t0​​(t)=1+(αΞ(y),t0​​(t))0​.

Coordinates of three-vectors are indexed by 0,1,20,1,20,1,2; projective coordinates are indexed by m,u,v,wm,u,v,wm,u,v,w.

For Y∈R3Y\in\mathbb R^3Y∈R3, put xY=Ξ(e(Y))x_Y=\Xi(e(Y))xY​=Ξ(e(Y)), βY(t)=Ψ0(e(Y),t)\beta_Y(t)=\Psi_0(e(Y),t)βY​(t)=Ψ0​(e(Y),t), and

HY={t>0:t∈Imax⁡(xY,0), (βY(t))v=1},T(Y)=inf⁡HY,F(Y)=η(βY(T(Y))).H_Y=\{t>0:t\in I_{\max}(x_Y,0),\ (\beta_Y(t))_v=1\}, \qquad T(Y)=\inf H_Y,\qquad F(Y)=\eta(\beta_Y(T(Y))).HY​={t>0:t∈Imax​(xY​,0), (βY​(t))v​=1},T(Y)=infHY​,F(Y)=η(βY​(T(Y))).

The infimum uses the total real-order convention. Write Valid⁡(Y)\operatorname{Valid}(Y)Valid(Y) when T(Y)T(Y)T(Y) is the least element of HYH_YHY​, every t∈[0,T(Y)]t\in[0,T(Y)]t∈[0,T(Y)] lies in Imax⁡(xY,0)I_{\max}(x_Y,0)Imax​(xY​,0) with Ae(Y),0(t)>0A_{e(Y),0}(t)>0Ae(Y),0​(t)>0, the active wall is strictly above one for 0<t<T(Y)0<t<T(Y)0<t<T(Y), and

dgαxY,0(T(Y))N(αxY,0(T(Y)))<0,g(x)=(Π(x))v.dg_{\alpha_{x_Y,0}(T(Y))}N(\alpha_{x_Y,0}(T(Y)))<0, \qquad g(x)=(\Pi(x))_v.dgαxY​,0​(T(Y))​N(αxY​,0​(T(Y)))<0,g(x)=(Π(x))v​.

All vector and operator norms below use the finite-product supremum norm. 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​.

Put c=u(v+w−u)c=u(v+w-u)c=u(v+w−u), s(m,u,v,w)=m−1−cs(m,u,v,w)=m-1-cs(m,u,v,w)=m−1−c, b(m,u,v,w)=(1−m2−(v−w)2+u2)/2b(m,u,v,w)=(1-m^2-(v-w)^2+u^2)/2b(m,u,v,w)=(1−m2−(v−w)2+u2)/2, D(m,u,v,w)=3+m2+K(u,v,w)D(m,u,v,w)=3+m^2+K(u,v,w)D(m,u,v,w)=3+m2+K(u,v,w), and r(m,v)=(m−1)2+v2r(m,v)=\sqrt{(m-1)^2+v^2}r(m,v)=(m−1)2+v2​.

Define

Cμ={(m,u,v,w):∣−m−μ∣≤1/2, 0≤u≤1/10000, 0≤v,w≤1},\mathcal C_\mu=\{(m,u,v,w):|-m-\mu|\leq1/2, \ 0\leq u\leq1/10000,\ 0\leq v,w\leq1\},Cμ​={(m,u,v,w):∣−m−μ∣≤1/2, 0≤u≤1/10000, 0≤v,w≤1}, ω−(y)=s(−ym,yu,yv,yw),H(y)=F(y)/ω−(y).\omega_-(y)=s(-y_m,y_u,y_v,y_w),\qquad H(y)=F(y)/\omega_-(y).ω−​(y)=s(−ym​,yu​,yv​,yw​),H(y)=F(y)/ω−​(y).

For a normalized trajectory write

β(t)=Π(αx,t0(t)),A(t)=1+(αx,t0(t))0,I=Imax⁡(x,t0).\beta(t)=\Pi(\alpha_{x,t_0}(t)),\quad A(t)=1+(\alpha_{x,t_0}(t))_0, \quad I=I_{\max}(x,t_0).β(t)=Π(αx,t0​​(t)),A(t)=1+(αx,t0​​(t))0​,I=Imax​(x,t0​).

For p=(μ,u,w)∈R3p=(\mu,u,w)\in\mathbb R^3p=(μ,u,w)∈R3, set e(p)=(−μ,u,1,w)e(p)=(-\mu,u,1,w)e(p)=(−μ,u,1,w), xp=Ξ(e(p))x_p=\Xi(e(p))xp​=Ξ(e(p)), and

T−(p)=inf⁡{t>0:t∈Imax⁡(xp,0), (Ψ0(e(p),t))w=1},P−(p)=Ψ0(e(p),T−(p)).T_-(p)=\inf\{t>0:t\in I_{\max}(x_p,0),\ (\Psi_0(e(p),t))_w=1\}, \qquad P_-(p)=\Psi_0(e(p),T_-(p)).T−​(p)=inf{t>0:t∈Imax​(xp​,0), (Ψ0​(e(p),t))w​=1},P−​(p)=Ψ0​(e(p),T−​(p)).

Let V−(p)V_-(p)V−​(p) mean that this infimum is a least element, the full interval [0,T−(p)][0,T_-(p)][0,T−​(p)] lies in the maximal domain with positive AAA, the growing wall is below one before the endpoint, and its normalized-time derivative is strictly positive at the endpoint. The infimum has its total real-order meaning.

Write e+(R,u,w)=(1+R2−1,u,1,w)e_+(R,u,w)=(1+\sqrt{R^2-1},u,1,w)e+​(R,u,w)=(1+R2−1​,u,1,w), e−(R,u,w)=(1−R2−1,u,1,w)e_-(R,u,w)=(1-\sqrt{R^2-1},u,1,w)e−​(R,u,w)=(1−R2−1​,u,1,w), p−(R,u,w)=(R2−1−1,u,w)p_-(R,u,w)=(\sqrt{R^2-1}-1,u,w)p−​(R,u,w)=(R2−1​−1,u,w), and R(m,u,v,w)=(−m,u,w,v)\mathcal R(m,u,v,w)=(-m,u,w,v)R(m,u,v,w)=(−m,u,w,v). The physical connection coordinates are

C(p)=((−P−(p)0−1)2+1,P−(p)1,P−(p)2).\mathscr C(p)=\left(\sqrt{(-P_-(p)_0-1)^2+1},P_-(p)_1,P_-(p)_2\right).C(p)=((−P−​(p)0​−1)2+1​,P−​(p)1​,P−​(p)2​).

Below TTT denotes the canonical core return time from e+e_+e+​ and F(Y)=η(Ψ0(e+(Y),T(Y)))F(Y)=\eta(\Psi_0(e_+(Y),T(Y)))F(Y)=η(Ψ0​(e+​(Y),T(Y))) its radial and wall coordinates. Let R(Y)=C(p−(Y))\mathscr R(Y)=\mathscr C(p_-(Y))R(Y)=C(p−​(Y)), tr(Y)=T−(p−(Y))t_r(Y)=T_-(p_-(Y))tr​(Y)=T−​(p−​(Y)), G=R∘FG=\mathscr R\circ FG=R∘F, and tG(Y)=T(Y)+tr(F(Y))t_G(Y)=T(Y)+t_r(F(Y))tG​(Y)=T(Y)+tr​(F(Y)). A cone direction at YYY has d0=1d_0=1d0​=1, ∣d1∣≤Y1|d_1|\leq\sqrt{Y_1}∣d1​∣≤Y1​​, ∣d2∣≤Y2|d_2|\leq\sqrt{Y_2}∣d2​∣≤Y2​​. For a map HHH with positive output walls, its positive cone property means dH(Y)ddH(Y)ddH(Y)d has positive radial component and ∣(dH(Y)d)i∣≤H(Y)i(dH(Y)d)0|(dH(Y)d)_i|\leq\sqrt{H(Y)_i}(dH(Y)d)_0∣(dH(Y)d)i​∣≤H(Y)i​​(dH(Y)d)0​ for i=1,2i=1,2i=1,2. 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​.

In the era formulas below, FFF denotes the core return map and G=R∘FG=\mathscr R\circ FG=R∘F its composition with reflection. For Z=(R,u,w)Z=(R,u,w)Z=(R,u,w) use the natural height h∗(Z)=η((R−2)/2)Jh_*(Z)=\eta((R-2)/2)^Jh∗​(Z)=η((R−2)/2)J and define S∗(Z)=P(h∗(Z),Z)S_*(Z)=\mathcal P(h_*(Z),Z)S∗​(Z)=P(h∗​(Z),Z), R∗(Z)=r(S∗(Z)0,S∗(Z)2)R_*(Z)=r(S_*(Z)_0,S_*(Z)_2)R∗​(Z)=r(S∗​(Z)0​,S∗​(Z)2​), Q∗(Z)=2/R∗(Z)Q_*(Z)=2/R_*(Z)Q∗​(Z)=2/R∗​(Z) and t∗(Z)=T(h∗(Z),Z)t_*(Z)=\mathcal T(h_*(Z),Z)t∗​(Z)=T(h∗​(Z),Z). The complete era maps are

Zm=F∘Gm,Sm=S∗∘Zm,Rm=R∗∘Zm,Pm=2/Rm.Z_m=F\circ G^m,\qquad S_m=S_*\circ Z_m,\qquad R_m=R_*\circ Z_m,\qquad P_m=2/R_m.Zm​=F∘Gm,Sm​=S∗​∘Zm​,Rm​=R∗​∘Zm​,Pm​=2/Rm​.

For a graph Γ=(U,W)\Gamma=(U,W)Γ=(U,W) write Γ[r]=(r,U(r),W(r))\Gamma[r]=(r,U(r),W(r))Γ[r]=(r,U(r),W(r)), ϕm(r)=Pm(Γ[r])\phi_m(r)=P_m(\Gamma[r])ϕm​(r)=Pm​(Γ[r]), km(r)=Rm(Γ[r])k_m(r)=R_m(\Gamma[r])km​(r)=Rm​(Γ[r]), um(r)=Sm(Γ[r])1u_m(r)=S_m(\Gamma[r])_1um​(r)=Sm​(Γ[r])1​, wm(r)=Sm(Γ[r])3w_m(r)=S_m(\Gamma[r])_3wm​(r)=Sm​(Γ[r])3​. Let Inv⁡a,b(f)\operatorname{Inv}_{a,b}(f)Inva,b​(f) denote Lean's total choice inverse invFunOn⁡f [a,b]\operatorname{invFunOn}f\,[a,b]invFunOnf[a,b]. For a=2(m+1+β)a=2(m+1+\beta)a=2(m+1+β), b=2(m+2−β)b=2(m+2-\beta)b=2(m+2−β), define

gm=Inv⁡a,b(ϕm),hm(y)=2/gm(y).g_m=\operatorname{Inv}_{a,b}(\phi_m),\qquad h_m(y)=2/g_m(y).gm​=Inva,b​(ϕm​),hm​(y)=2/gm​(y).

Write γn(y)=1/(n+y)\gamma_n(y)=1/(n+y)γn​(y)=1/(n+y) for the Gauss inverse branch. The retained phase set is HN=⋃m<N(hm(1−2β),hm(2β))\mathcal H_N=\bigcup_{m<N}(h_m(1-2\beta),h_m(2\beta))HN​=⋃m<N​(hm​(1−2β),hm​(2β)).

A curve Γ=(U,W)\Gamma=(U,W)Γ=(U,W) is admissible at budget BBB if, for every r>2r>2r>2, U(r),W(r)>0U(r),W(r)>0U(r),W(r)>0, U(r)+W(r)≤BU(r)+W(r)\leq BU(r)+W(r)≤B, both functions are C1C^1C1 at rrr, and ∣U′(r)∣≤U(r)|U'(r)|\leq\sqrt{U(r)}∣U′(r)∣≤U(r)​, ∣W′(r)∣≤W(r)|W'(r)|\leq\sqrt{W(r)}∣W′(r)∣≤W(r)​. For an arbitrary graph, define its successor walls by

Um+=Ekm,uma,b; 2/(1−2β), 1/β,Wm+=Ekm,wma,b; 2/(1−2β), 1/β,U_m^+=\mathfrak E_{k_m,u_m}^{a,b;\,2/(1-2\beta),\,1/\beta},\qquad W_m^+=\mathfrak E_{k_m,w_m}^{a,b;\,2/(1-2\beta),\,1/\beta},Um+​=Ekm​,um​a,b;2/(1−2β),1/β​,Wm+​=Ekm​,wm​a,b;2/(1−2β),1/β​,

where E\mathfrak EE is the canonical amplitude extension after composing with the choice inverse of kmk_mkm​, with slope um′/km′u_m'/k_m'um′​/km′​ or wm′/km′w_m'/k_m'wm′​/km′​. The one-step property at (η,J,σ,B,C)(\eta,J,\sigma,B,C)(η,J,σ,B,C) consists of the following. Put β=Bσ\beta=B^\sigmaβ=Bσ, N=⌊1/β⌋NN=\lfloor1/\beta\rfloor_{\mathbb N}N=⌊1/β⌋N​, Kβ=[2β,1−2β]K_\beta=[2\beta,1-2\beta]Kβ​=[2β,1−2β]. One has 0<β<1/1600<\beta<1/1600<β<1/160. For every m<Nm<Nm<N, hmh_mhm​ is strictly decreasing on KβK_\betaKβ​, and for y∈Kβy\in K_\betay∈Kβ​,

1m+2<hm(y)<1m+1,ϕm(2/hm(y))=y,hm is C1 at y,hm′(y)<0,\frac1{m+2}<h_m(y)<\frac1{m+1},\quad \phi_m(2/h_m(y))=y, \quad h_m\text{ is }C^1\text{ at }y,\quad h_m'(y)<0,m+21​<hm​(y)<m+11​,ϕm​(2/hm​(y))=y,hm​ is C1 at y,hm′​(y)<0, ∣hm(y)−γm+1(y)∣≤CB/(m+1)2,∣log⁡∣hm′(y)∣∣γm+1′(y)∣∣≤5CB.|h_m(y)-\gamma_{m+1}(y)|\leq CB/(m+1)^2,\quad \left|\log\frac{|h_m'(y)|}{|\gamma_{m+1}'(y)|}\right|\leq5C\sqrt B.∣hm​(y)−γm+1​(y)∣≤CB/(m+1)2,​log∣γm+1′​(y)∣∣hm′​(y)∣​​≤5CB​.

The successor (Um+,Wm+)(U_m^+,W_m^+)(Um+​,Wm+​) is admissible at budget B4/3B^{4/3}B4/3 and, with r=2/hm(y)r=2/h_m(y)r=2/hm​(y), satisfies km(r)=2/yk_m(r)=2/ykm​(r)=2/y, Um+(2/y)=um(r)U_m^+(2/y)=u_m(r)Um+​(2/y)=um​(r), Wm+(2/y)=wm(r)W_m^+(2/y)=w_m(r)Wm+​(2/y)=wm​(r). The retained branch intervals are pairwise disjoint and λ((0,1)∖HN)≤10Bσ\lambda((0,1)\setminus\mathcal H_N)\leq10B^\sigmaλ((0,1)∖HN​)≤10Bσ. Define

χ(t)=max⁡(1−t,0),P(t)=∫0tχ(u) du,Rs,d,e(x)=s+edsP(ex/s).\chi(t)=\max(1-t,0),\qquad P(t)=\int_0^t\chi(u)\,du,\qquad \mathcal R_{s,d,e}(x)=s+edsP(ex/s).χ(t)=max(1−t,0),P(t)=∫0t​χ(u)du,Rs,d,e​(x)=s+edsP(ex/s).

The integral is oriented, and division uses the total reciprocal 0−1=00^{-1}=00−1=0.

For real functions f,dff,d_ff,df​ and real endpoints l,rl,rl,r, define FFF and VVV by

F(x)={Rf(l),df(l),−1(x−l),x<l,Rf(r),df(r),1(x−r),x≥l, r<x,f(x),otherwise,F(x)=\begin{cases} \mathcal R_{f(l),d_f(l),-1}(x-l),&x<l,\\ \mathcal R_{f(r),d_f(r),1}(x-r),&x\geq l,\ r<x,\\ f(x),&\text{otherwise}, \end{cases}F(x)=⎩⎨⎧​Rf(l),df​(l),−1​(x−l),Rf(r),df​(r),1​(x−r),f(x),​x<l,x≥l, r<x,otherwise,​ V(x)={df(l)χ(−(x−l)/f(l)),x<l,df(r)χ((x−r)/f(r)),x≥l, r<x,df(x),otherwise.V(x)=\begin{cases} d_f(l)\chi(-(x-l)/f(l)),&x<l,\\ d_f(r)\chi((x-r)/f(r)),&x\geq l,\ r<x,\\ d_f(x),&\text{otherwise}. \end{cases}V(x)=⎩⎨⎧​df​(l)χ(−(x−l)/f(l)),df​(r)χ((x−r)/f(r)),df​(x),​x<l,x≥l, r<x,otherwise.​

Write Ff,VfF_f,V_fFf​,Vf​ when comparing several pairs of input functions.

For real functions U,dUU,d_UU,dU​, set

rU(x)=U(x),bU(x)=dU(x)2U(x).r_U(x)=\sqrt{U(x)},\qquad b_U(x)=\frac{d_U(x)}{2\sqrt{U(x)}}.rU​(x)=U(x)​,bU​(x)=2U(x)​dU​(x)​.

Let FU,VUF_U,V_UFU​,VU​ be the root-curve construction applied to (rU,bU)(r_U,b_U)(rU​,bU​), and set

U~(x)=FU(x)2.\widetilde U(x)=F_U(x)^2.U(x)=FU​(x)2.

For another pair W,dWW,d_WW,dW​, use the analogous notation.

Let g=invFunOn⁡k [a,b]g=\operatorname{invFunOn}k\,[a,b]g=invFunOnk[a,b], A(t)=U(g(t))A(t)=U(g(t))A(t)=U(g(t)), and V(t)=U′(g(t))/k′(g(t))V(t)=U'(g(t))/k'(g(t))V(t)=U′(g(t))/k′(g(t)). Denote by Ek,Ua,b;l,r\mathfrak E_{k,U}^{a,b;l,r}Ek,Ua,b;l,r​ the canonical amplitude extension of A,VA,VA,V from [l,r][l,r][l,r], and by Ek,Ua,b\mathfrak E_{k,U}^{a,b}Ek,Ua,b​ the extension from [k(b),k(a)][k(b),k(a)][k(b),k(a)].

The data includes the scalar local-inverse predicate: a C1C^1C1 inverse at f(a)f(a)f(a) taking f(a)f(a)f(a) to aaa, with derivative 1/f′(a)1/f'(a)1/f′(a) and two-sided inverse germs. It also includes the diagonal map Sk(R,U,P)=(R,kU,P)S_k(R,U,P)=(R,kU,P)Sk​(R,U,P)=(R,kU,P), the graph restrictions, the retained phase unions, and the complete one-step property of the era curve.

Definition code
import Mathlib.Analysis.Calculus.InverseFunctionTheorem.ContDiff
import Mathlib.Analysis.Calculus.InverseFunctionTheorem.Deriv
import Mathlib.Tactic.NormNum
import Mathlib.Analysis.SpecialFunctions.Log.Basic
import Mathlib.Tactic.FieldSimp
import Mathlib.Tactic.Linarith
import Mathlib.Tactic.Positivity
import Mathlib.Tactic.Ring
import Mathlib.Analysis.Calculus.Deriv.Prod
import Mathlib.LinearAlgebra.Matrix.Notation
import Mathlib.Tactic.FinCases
import Mathlib.Analysis.Calculus.Deriv.MeanValue
import Mathlib.Analysis.Calculus.ContDiff.Deriv
import Mathlib.Analysis.Calculus.Deriv.Inv
import Mathlib.Analysis.Calculus.Deriv.Add
import Mathlib.Analysis.SpecificLimits.Basic
import Mathlib.Data.List.Induction
import Mathlib.Algebra.BigOperators.Ring.Finset
import Mathlib.MeasureTheory.Measure.Lebesgue.Basic
import Definitions.Def_BianchiIX_Analysis_AmplitudeExtension
import Mathlib.MeasureTheory.Integral.IntervalIntegral.FundThmCalculus
import Mathlib.Analysis.Calculus.Deriv.Pow
import Mathlib.Analysis.Calculus.Deriv.Mul
import Mathlib.Analysis.Real.Sqrt
import Mathlib.Analysis.Calculus.Deriv.Slope
import Mathlib.MeasureTheory.Constructions.BorelSpace.Metrizable
import Mathlib.Analysis.SpecialFunctions.Pow.Continuity
import Mathlib.Analysis.SpecialFunctions.Sqrt
import Definitions.Def_BianchiIX_Dynamics_CoreAngle
import Mathlib.Analysis.Calculus.FDeriv.Symmetric
import Mathlib.Tactic.Module
import Mathlib.Analysis.SpecialFunctions.Arsinh
import Mathlib.Analysis.Calculus.ContDiff.Operations
import Mathlib.Analysis.Calculus.MeanValue
import Mathlib.Data.Fin.VecNotation
import Mathlib.Tactic.FunProp
import Mathlib.Tactic.LinearCombination
import Mathlib.Analysis.SpecialFunctions.Arcosh
import Mathlib.Analysis.SpecialFunctions.Pow.Deriv
import Mathlib.Analysis.SpecialFunctions.Integrals.Basic
import Mathlib.Analysis.ODE.ExistUnique
import Mathlib.Data.Matrix.Mul
import Mathlib.Algebra.BigOperators.Fin
import Mathlib.Analysis.Calculus.FDeriv.Partial
import Mathlib.Analysis.Calculus.ImplicitContDiff
import Mathlib.Analysis.ODE.Gronwall
import Mathlib.Analysis.SpecialFunctions.Trigonometric.ArctanDeriv
import Mathlib.Analysis.SpecialFunctions.Trigonometric.DerivHyp
import Mathlib.Topology.Instances.Real.Lemmas
import Mathlib.MeasureTheory.Integral.IntervalIntegral.DistLEIntegral
import Mathlib.Tactic.GCongr
import Mathlib.Topology.Order.OrderClosed
import Mathlib.Topology.Order.IntermediateValue
import Definitions.Def_BianchiIX_Dynamics_ProjectiveState
import Definitions.Def_BianchiIX_Dynamics_CoreMultipliers
import Mathlib.Analysis.Calculus.ParametricIntervalIntegral
import Mathlib.Topology.Order.Compact
import Mathlib.MeasureTheory.Integral.DominatedConvergence
import Mathlib.Data.Real.Basic
import Definitions.Def_BianchiIX_Dynamics_ProjectiveRegularity
import Definitions.Def_BianchiIX_Dynamics_MaximalFlow
import Definitions.Def_BianchiIX_Dynamics_Events
import Definitions.Def_BianchiIX_Dynamics_ProjectiveFlowVariation
import Mathlib.MeasureTheory.Integral.IntervalIntegral.IntegrationByParts
import Definitions.Def_BianchiIX_Dynamics_ConnectionReturn
import Definitions.Def_BianchiIX_Dynamics_ConnectionPhysicalReturn
import Definitions.Def_BianchiIX_Dynamics_CoreReturnDifferential
import Definitions.Def_BianchiIX_Dynamics_CoreConnectionChart
import Mathlib.Analysis.SpecialFunctions.Log.Deriv
import Mathlib.Analysis.SpecialFunctions.ExpDeriv
import Mathlib.Topology.Order.ProjIcc
import Mathlib.Analysis.Normed.Operator.Basic
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 Definitions.Def_BianchiIX_Dynamics_CompactConnectionReturn
import Definitions.Def_BianchiIX_EntrancePhaseResponseData
import Mathlib.Topology.MetricSpace.Thickening
import Mathlib.Topology.MetricSpace.Cauchy
import Mathlib.Analysis.Calculus.ContDiff.Comp
import Mathlib.Analysis.Calculus.Deriv.Comp
import Mathlib.Analysis.Calculus.UniformLimitsDeriv
import Definitions.Def_BianchiIX_Dynamics_CompactCoreReturn
import Mathlib.Analysis.SpecialFunctions.ImproperIntegrals
import Mathlib.Analysis.Calculus.ParametricIntegral
import Mathlib.MeasureTheory.Integral.Bochner.ContinuousLinearMap
import Mathlib.Analysis.SpecialFunctions.Pow.Real
import Definitions.Def_BianchiIX_TerminalConnectionData
import Definitions.Def_BianchiIX_TerminalConeData
import Mathlib.Analysis.Calculus.FDeriv.Mul
import Mathlib.LinearAlgebra.Matrix.ToLin
import Mathlib.Analysis.Calculus.FDeriv.Comp
import Mathlib.Algebra.Order.Floor.Ring
import Mathlib.Algebra.Order.Floor.Semiring
import Definitions.Def_BianchiIX_Statistics_GaussInverse
import Definitions.Def_BianchiIX_Analysis_AmplitudeSequence
import Mathlib.Topology.Algebra.InfiniteSum.Order
import Mathlib.Data.PNat.Basic
import Mathlib.MeasureTheory.Constructions.BorelSpace.Basic
import Mathlib.MeasureTheory.Function.Floor

section

noncomputable section

namespace BianchiIX

open Filter
open scoped Topology

structure ScalarLocalInverseAt (f : ℝ → ℝ) (a : ℝ) (g : ℝ → ℝ) : Prop where
  apply_base : g (f a) = a
  contDiffAt : ContDiffAt ℝ 1 g (f a)
  hasDerivAt : HasDerivAt g (deriv f a)⁻¹ (f a)
  left_inverse : ∀ᶠ x in 𝓝 a, g (f x) = x
  right_inverse : ∀ᶠ y in 𝓝 (f a), f (g y) = y







end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def graphPhase (F : (Fin 3 → ℝ) → ℝ) (U W : ℝ → ℝ) (r : ℝ) : ℝ := F ![r, U r, W r]





end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set Filter
open scoped Topology

def intervalRadiusInverse (k : ℝ → ℝ) (a b : ℝ) : ℝ → ℝ := Function.invFunOn k (Icc a b)



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set MeasureTheory
open scoped ENNReal

def retainedBranchPhases (N : ℕ) (L R : ℕ → ℝ) : Set ℝ :=
  ⋃ m ∈ Finset.range N, Ioo (R m) (L m)



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set

def radiusAmplitude (k U : ℝ → ℝ) (a b x : ℝ) : ℝ := U (intervalRadiusInverse k a b x)

def radiusAmplitudeSlope (k U : ℝ → ℝ) (a b x : ℝ) : ℝ :=
  deriv U (intervalRadiusInverse k a b x) / deriv k (intervalRadiusInverse k a b x)

def radiusAmplitudeExtension (k U : ℝ → ℝ) (a b : ℝ) : ℝ → ℝ :=
  amplitudeCurveExtension (k b) (k a) (radiusAmplitude k U a b) (radiusAmplitudeSlope k U a b)





end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set

def radiusAmplitudeWindowExtension (k U : ℝ → ℝ) (a b l r : ℝ) : ℝ → ℝ :=
  amplitudeCurveExtension l r (radiusAmplitude k U a b) (radiusAmplitudeSlope k U a b)



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set Filter
open scoped Topology

def scaleCoreVariationLinear (k : ℝ) : CoreAngleState →L[ℝ] CoreAngleState :=
  ContinuousLinearMap.pi ![ContinuousLinearMap.proj 0,
    k • ContinuousLinearMap.proj 1, ContinuousLinearMap.proj 2]



















end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set MeasureTheory

/-- Positive multipliers from an arbitrary projective initial point to its moving exit. -/
def partialCoreReturnMultiplier (τ : ProjectiveState → ℝ) (y : ProjectiveState) (i : Fin 2) : ℝ :=
  coreFlowMultiplier y 0 0 (τ y) i







end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set

def partialCoreVariationBudget (k u p : ℝ) (v : CoreAngleState) (dz : ℝ) : ℝ :=
  |v 0| + k * |v 1| + |v 2| + (k ^ 2 * u + p / k) * |dz|













end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set Filter
open scoped Topology

def partialCoreDownwardReturnTime (y : ProjectiveState) : ℝ :=
  sInf {t : ℝ | 0 < t ∧ t ∈ normalizedMaxDomain (fromProjectiveState y) 0 ∧
    projectiveOriginalFlow y 0 t 2 = 1 ∧
    normalizedEventRate projectiveActiveWall (fromProjectiveState y) 0 t < 0}

def partialCoreDownwardReturnState (y : ProjectiveState) : CoreAngleState :=
  toCoreAngleState (projectiveOriginalFlow y 0 (partialCoreDownwardReturnTime y))









end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set

def coreConnectionParameters (τ : CoreAngleState → ℝ) (Y : CoreAngleState) : ConnectionParameters :=
  coreExitConnectionParameters (coreReturnState τ Y)

def coreReflectionState (τ : CoreAngleState → ℝ) (Y : CoreAngleState) : CoreAngleState :=
  connectionReturnPhysicalState (coreConnectionParameters τ Y)





end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set Filter
open scoped Topology

def reflectionReturnTime (Y : CoreAngleState) : ℝ :=
  connectionReturnTime (coreExitConnectionParameters Y)

def reflectionReturnState (Y : CoreAngleState) : CoreAngleState :=
  connectionReturnPhysicalState (coreExitConnectionParameters Y)

def compactConnectionPhysicalState (η : ℝ) (Y : CoreAngleState) : CoreAngleState :=
  ![compactConnectionPhase η Y, compactConnectionState η Y 1, compactConnectionState η Y 3]











end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def coreReflectionTime (Y : CoreAngleState) : ℝ :=
  reflectionReturnTime (coreReturnState coreRadiusReturnTime Y) + coreRadiusReturnTime Y



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set Filter
open scoped Topology

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

def terminalNaturalBoundedPhase (η : ℝ) (J : ℕ) (Y : ExitTailState) : ℝ :=
  2 / terminalNaturalPhase η J Y





end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set

/-- The last downward core state after `m` nonterminal reflections. -/
def eraCoreExit (m : ℕ) : CoreAngleState → CoreAngleState :=
  coreReturnState coreRadiusReturnTime ∘ (coreReflectionState coreRadiusReturnTime)^[m]

def eraBoundedPhase (η : ℝ) (J m : ℕ) : CoreAngleState → ℝ :=
  terminalNaturalBoundedPhase η J ∘ eraCoreExit m



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

def eraOutputState (η : ℝ) (J m : ℕ) : CoreAngleState → ProjectiveState :=
  terminalNaturalState η J ∘ eraCoreExit m

def eraOutputRadius (η : ℝ) (J m : ℕ) : CoreAngleState → ℝ :=
  terminalNaturalPhase η J ∘ eraCoreExit m



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set

def eraGraphRadius (η : ℝ) (J m : ℕ) (U W : ℝ → ℝ) : ℝ → ℝ :=
  graphPhase (eraOutputRadius η J m) U W

def eraGraphWall (η : ℝ) (J m : ℕ) (U W : ℝ → ℝ) (i : Fin 4) : ℝ → ℝ :=
  graphPhase (fun X => eraOutputState η J m X i) U W



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set

def eraInverseBranch (η : ℝ) (J m : ℕ) (β : ℝ) (U W : ℝ → ℝ) : ℝ → ℝ :=
  fun y => 2 / intervalRadiusInverse (graphPhase (eraBoundedPhase η J m) U W)
    (2 * ((m : ℝ) + 1 + β)) (2 * ((m : ℝ) + 2 - β)) y



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set

def eraAmplitudeExtension (η : ℝ) (J m : ℕ) (β : ℝ) (U W : ℝ → ℝ) (i : Fin 4) : ℝ → ℝ :=
  radiusAmplitudeWindowExtension (eraGraphRadius η J m U W) (eraGraphWall η J m U W i)
    (2 * ((m : ℝ) + 1 + β)) (2 * ((m : ℝ) + 2 - β)) (2 / (1 - 2 * β)) (1 / β)



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set MeasureTheory
open scoped ENNReal

def eraRetainedPhases (η : ℝ) (J N : ℕ) (β : ℝ) (U W : ℝ → ℝ) : Set ℝ :=
  retainedBranchPhases N (fun m => eraInverseBranch η J m β U W (2 * β))
    (fun m => eraInverseBranch η J m β U W (1 - 2 * β))



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set MeasureTheory
open scoped ENNReal

structure EraCurve where
  U : ℝ → ℝ
  W : ℝ → ℝ

namespace EraCurve

def state (Γ : EraCurve) (r : ℝ) : CoreAngleState := ![r, Γ.U r, Γ.W r]

def Admissible (Γ : EraCurve) (B : ℝ) : Prop :=
  ∀ r : ℝ, 2 < r → 0 < Γ.U r ∧ 0 < Γ.W r ∧ Γ.U r + Γ.W r ≤ B ∧
    ContDiffAt ℝ 1 Γ.U r ∧ ContDiffAt ℝ 1 Γ.W r ∧
    |deriv Γ.U r| ≤ Real.sqrt (Γ.U r) ∧ |deriv Γ.W r| ≤ Real.sqrt (Γ.W r)

def successor (Γ : EraCurve) (η : ℝ) (J m : ℕ) (β : ℝ) : EraCurve :=
  ⟨eraAmplitudeExtension η J m β Γ.U Γ.W 1, eraAmplitudeExtension η J m β Γ.U Γ.W 3⟩

/-- Uniform data for one actual era, its selected inverse branches, and its successor curves. -/
structure StepEstimates (Γ : EraCurve) (η : ℝ) (J : ℕ) (σ B C : ℝ) : Prop where
  beta_pos : 0 < B ^ σ
  beta_small : B ^ σ < 1 / 160
  branch : ∀ m : ℕ, m < ⌊1 / B ^ σ⌋₊ →
    let h := eraInverseBranch η J m (B ^ σ) Γ.U Γ.W
    let K := Icc (2 * B ^ σ) (1 - 2 * B ^ σ)
    StrictAntiOn h K ∧ ∀ y ∈ K,
      h y ∈ Ioo (1 / ((m : ℝ) + 2)) (1 / ((m : ℝ) + 1)) ∧
      graphPhase (eraBoundedPhase η J m) Γ.U Γ.W (2 / h y) = y ∧
      ContDiffAt ℝ 1 h y ∧ deriv h y < 0 ∧
      |h y - gaussInv ((m : ℝ) + 1) y| ≤ C * B / ((m : ℝ) + 1) ^ 2 ∧
      |Real.log (|deriv h y| / |deriv (gaussInv ((m : ℝ) + 1)) y|)| ≤ 5 * C * Real.sqrt B
  successor_admissible : ∀ m : ℕ, m < ⌊1 / B ^ σ⌋₊ →
    (Γ.successor η J m (B ^ σ)).Admissible (B ^ (4 / 3 : ℝ))
  successor_match : ∀ m : ℕ, m < ⌊1 / B ^ σ⌋₊ → ∀ y ∈ Icc (2 * B ^ σ) (1 - 2 * B ^ σ),
    let x := 2 / eraInverseBranch η J m (B ^ σ) Γ.U Γ.W y
    eraOutputRadius η J m (Γ.state x) = 2 / y ∧
      (Γ.successor η J m (B ^ σ)).U (2 / y) = eraOutputState η J m (Γ.state x) 1 ∧
      (Γ.successor η J m (B ^ σ)).W (2 / y) = eraOutputState η J m (Γ.state x) 3
  branch_disjoint : (↑(Finset.range ⌊1 / B ^ σ⌋₊) : Set ℕ).PairwiseDisjoint
    (fun m => Ioo (eraInverseBranch η J m (B ^ σ) Γ.U Γ.W (1 - 2 * B ^ σ))
      (eraInverseBranch η J m (B ^ σ) Γ.U Γ.W (2 * B ^ σ)))
  stopped_volume : volume (Ioo (0 : ℝ) 1 \ eraRetainedPhases η J ⌊1 / B ^ σ⌋₊ (B ^ σ) Γ.U Γ.W) ≤
    ENNReal.ofReal (10 * B ^ σ)



end EraCurve



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set

namespace EraCurve

def constant (B : ℝ) : EraCurve := ⟨fun _ => B / 2, fun _ => B / 2⟩



/-- Unretained digits have empty branch domains; a constant curve makes their unused descendants total. -/
def nextCurve (Γ : EraCurve) (η : ℝ) (J : ℕ) (σ B : ℝ) (a : ℕ+) : EraCurve :=
  if a.natPred < ⌊1 / B ^ σ⌋₊ then Γ.successor η J a.natPred (B ^ σ)
  else constant (B ^ (4 / 3 : ℝ))

def family (Γ : EraCurve) (η : ℝ) (J : ℕ) (σ B : ℝ) (w : List ℕ+) : EraCurve :=
  w.reverseRecOn Γ (fun v a γ => γ.nextCurve η J σ (eraAmplitude B v.length) a)









end EraCurve



end BianchiIX

end
end

section

noncomputable section

namespace BianchiIX

open Set

def eraNextState (η : ℝ) (J m : ℕ) (X : CoreAngleState) : CoreAngleState :=
  ![eraOutputRadius η J m X, eraOutputState η J m X 1, eraOutputState η J m X 3]

def eraBudgetReturn (η : ℝ) (J m : ℕ) (p : ℝ × CoreAngleState) : ℝ × CoreAngleState :=
  (p.1 ^ (4 / 3 : ℝ), eraNextState η J m p.2)

def eraInputSelection (σ Bstar : ℝ) (m : ℕ) : Set (ℝ × CoreAngleState) :=
  {p | p.1 ∈ Ioc 0 Bstar ∧ m < ⌊1 / p.1 ^ σ⌋₊ ∧
    0 < p.2 1 ∧ 0 < p.2 2 ∧ p.2 1 + p.2 2 ≤ p.1 ∧
    p.2 0 ∈ Icc (2 * ((m : ℝ) + 1 + p.1 ^ σ)) (2 * ((m : ℝ) + 2 - p.1 ^ σ))}

def eraSelectedInputs (η : ℝ) (J : ℕ) (σ Bstar : ℝ) (m : ℕ) : Set (ℝ × CoreAngleState) :=
  {p | p ∈ eraInputSelection σ Bstar m ∧
    2 / eraOutputRadius η J m p.2 ∈ Ioo (2 * p.1 ^ σ) (1 - 2 * p.1 ^ σ)}











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