The Notebook Grid: Exact Addresses, Three Reflections, Pell Branches, and Rational Carriers

Author: Jeffery Huckstead · Cerebral Graphix (cerebralgraphix.com) · ORCID 0009-0007-0234-2177 Version 4.1 · Reader revised 2 November 2026 · DOI 10.5281/zenodo.23839329 Companion to Notebook_Grid_v4_1.pdf (54 Pages · 15 Figures · DOI 10.5281/zenodo.23839329) Offline Self-Contained HTML Companion
Quick Section Navigation & Text Filter (Notebook_Grid_v4_1.pdf Companion View):

The Notebook Grid: Exact Addresses, Three Reflections, Pell Branches, and Rational Carriers

with Odd-Square Factorizations, Stack Shadows, and Divisor-Shell Arithmetic

Jeffery Huckstead · Cerebral Graphix (cerebralgraphix.com) · ORCID 0009-0007-0234-2177

Version 4.1 · November 2, 2026

Abstract

Writing the positive integers in successive triangular chambers gives an exact address N(r,c) = Tr+c-2 + r for every cell of the positive integer quadrant. We derive its two inverse formulas, characterize the admissible chamber-offset wedge (m,V), and prove the quadratic address families traced by affine cell paths. Holding a chamber m = r+c-1 fixed, three natural phases—the native row fraction κrow = r/m, the cell-midpoint fraction κmid = (r-1/2)/m, and the nonlinear square-root phase κroot = √(2N)-m+1/2—induce three distinct reflection maps under κ ↦ 1-κ. Only the midpoint reflection Rmid(N) = m2+1-N permutes a whole chamber; the native reflection Rrow(N) = m2-N preserves the proper subset 1 ≤ r ≤ m-1; and the square-root reflection Rroot(N) = (2m-√(2N))2/2 has an integer image if and only if the cell lies on the superdiagonal N = 2a2 (m=2a), which it fixes. For scaled-square paths Nk = A(Qk)2, a floor correction ℓk partitions the path into exact polynomial branches, separates the first branch departure from sparse power samples, and identifies norm-one triangular seams, algebraic returns at k=2P, and complete negative-Pell diagonal squares. Applying the midpoint fold 𝔣(κmid) = r̃c̃/(4m2) defines an upper-semicircle carrier map 𝒞(r,c) whose rational points are parameterized uniquely by odd squares (r̃,c̃) = λ(U2,W2) with primitive first chamber m0 = (U2+W2)/2 ≡ 1 (mod 4). Its least nontrivial chamber is m0 = 5, and its 45∘ Pell record sequence m0,j ∈ {1,5,29,169,985,…} maps bijectively to diagonal chambers m = Gj ∈ {1,7,41,239,1393,…} carrying square addresses N = m0,j2. On the address sheet pN = (-r/(m+1), 1/(m+1)), the identity (2n+1)2 = 8Tn + 1 = 4n(n+1) + 1 unites the eight-sector odd-square partition, the balanced split {n, n+1} and its finite reciprocal product, Euler's cosine generating function, a fixed-gap Gamma ratio, and depth-weighted isometric stack projection, while the same addresses carry divisor-band progressions, k-dimensional Jordan-totient shells, the factorial-residual dichotomy Q(N) = N/gcd(N,m!) ∈ {1}∪ℙ of asymptotic density log 2, and an entire-function realization of the triangular Taylor coefficients. Two self-contained appendices establish the complete rational-carrier chamber-budget theorems and the all-direction union-shadow theory of the three-dimensional odd-square stack Ωn.

Route Through the Paper. The exposition follows a single chain of maps, domains, and declared measures:

  1. Triangular Address & Chamber Coordinates (§1–§2): Positive quadrant cells (r,c) map bijectively to addresses N(r,c)=Tm-1+r in chamber m=r+c-1 with offset V=r-c and odd legs (r̃,c̃)=(2r-1,2c-1), turning affine lattice rays into degree-2 polynomial address sequences with period-4 offset parity.
  2. Chamber Phases, Three Reflections & Scaled-Square Pell Branches (§3–§4): Fixing m distinguishes κrow, κmid, and κroot (with the anti-diagonal pair correspondence (k, p-k)), whose reflections and floor-branch drift ℓk = ⌊1/2+kε⌋ govern scaled-square Pell branches, norm-one triangular seams, and diagonal square returns.
  3. Rational Carriers, Address Sheet & Odd-Square Factorization (§5–§6, Apps. A–B): Folding κmid lifts odd-square cells onto rational points of the upper semicircle 𝒞(r,c), while row normalization embeds the Grid into Q=[-1,0]×[0,1], separating primitive ratio from scale height, active-span rescaling (1/m) from staircase fixed points (1/k), flat screen sampling (2ω) from depth-weighted stack sampling (6ω(1-ω)), and uniting five routes on (2n+1)2 = 8Tn + 1.
  4. Open Questions, Divisor Shells, Jordan Totients & Factorial Residuals (§7–§8): Separating finite and open prime-pair and elliptic-curve questions (§7) from proved arithmetic on the addresses (§8) yields divisor-shell progressions, three separate Jordan-totient zeta limits, the complete factorial-residual classification Q(N)∈{1}∪ℙ with Dickman density log 2, and an entire-function differential lift of Tn/n!.

Scope, Proofs, and Verification. Every address, reflection, floor-branch, rational-carrier, address-sheet, odd-square factorization, Brouncker continuant, Euler-transform, divisor-shell, Jordan-totient, factorial-residual, entire-function coefficient-lift, chamber-budget, and three-dimensional stack-shadow theorem in §§1–6, §8, and Appendices A–B is proved algebraically, analytically, or by positive-integer descent in this paper. Classical completeness classifications on the rails (§4.6) and smooth-number asymptotics (§8.4) cite Nagell [1], Leu and Li [2], Cohn [3], and de Bruijn [13] with explicit hypotheses (with a self-contained 5-adic valuation proof for base 5 on the diagonal). Finite-range numerical and statistical statements carry explicit bounds: the shell-numerator diagnostic table (§8.2, Remark 8.1a) is a deterministic finite calculation at M ∈ {250, 500, 1000, 2000}; the odd-chamber mirror-prime verification (§7.1) covers the finite interval 1403 ≤ m ≤ 20001; and all-index mirror-prime existence (m ≥ 1403), all-index non-intersection of the a=2 and 45∘ record-chamber families beyond m0 ∈ {1,5,29} (§7.2), two-point prime occupancy (§8.5), and the analytic questions of §8.6 remain open. Throughout, the finite odd-side (j=2n+1) atomic screen-maximum laws Prflat(M=i)=(2i-1)/j2 and Prstack(M=i)=(2i-1)[1+⌊(j-i)/2⌋]/C(j+2, 3) are kept distinct from the j→∞ weak-limit densities 2ω and 6ω(1-ω) of M/j; free congruence of the eight half-open max-norm sectors Σ0,…,Σ7 in the lattice-preserving Euclidean group E(2,ℤ) = ℤ2 ⋊ O(2,ℤ) ⊂ E(2) = ℝ2 ⋊ O(2) via F(a,b)=(b+1,a) is kept distinct from origin-fixing O(2) non-equivalence; and the orthochronous boost conjugacy CD LD CD-1 = ΛβD ∈ SO+(1,1) is restricted to the positive Lorentz sheet P1 > 0 (the Cayley images μ(1/√D) are fundamental integer Pell units exactly for D ∈ {2,3}). An optional standard-library Python script (check_notebook_grid_v4_1.py) accompanies the supplementary archive for finite-range computational corroboration; no proof step in the paper depends on executing code.

Summary of Main Results

SectionCore Object & Local CoordinatesPrincipal Theorem / IdentityProof Method & Scope
§1. Addresses & ChambersCell (r,c) ∈ ℤ≥ 12, chamber m=r+c-1, offset V=r-c, odd legs r̃=2r-1, c̃=2c-1N(r,c)=Tm-1+r, 2N=m2+V+1; admissible iff m ≥ 1, | V| ≤ m-1, m+V ≡ 1 (mod 2); m = ⌈(√(8N+1)-1)/2⌉ = ⌊√(2N)+1/2⌋Props. 1.1–1.2, Fig. 1 (Algebraic proof)
§2. Affine Paths & RailsBaseline (r0,c0), increments (Δr,Δc), L=r0+c0-2, D=Δr+ΔcQuadratic address polynomial of degree ≤ 2; five rails; fixed-V rail parity N(r,r-V) ≡ TV+1 (mod 2) with period 4 (1,1,0,0)Props. 2.1–2.2, Fig. 2 (Algebraic proof)
§3. Three ReflectionsPhases κrow=r/m, κmid=(r-1/2)/m, κroot=√(2N)-m+1/2 (ε=κroot-1/2); pair (k,p-k)Rrow(N)=m2-N, Rmid(N)=m2+1-N, Rroot(N)=(2m-√(2N))2/2; Rroot(N)∈ℤ ⇔ N=2a2; 3/(8m) < κroot-κmid ≤ 1/(2m); (k,p-k)↔ RmidThm. 3.1, Props. 3.2–3.3, Rem. 3.4, Fig. 3 (Algebraic proof)
§4. Floor Branches & PellNk=A(Qk)2, P=round(Q√(2A)), ε=Q√(2A)-P, δ=P2-2AQ2, ℓk=⌊1/2+kε⌋mk=Pk+ℓk, Vk=-δ k2-2Pkℓk-ℓk2-1; ℓk=0 ⇔ -P ≤ δ k < P; seams k=P, returns k=2P; diagonal descent m2-2Q2=-1; rail powersProps. 4.1–4.8, Fig. 3b (Algebraic & descent proofs; [1–3] + 5-adic proof in §4.6)
§5. Rational CarriersMidpoint circle 𝒞(r,c)=(V/m,(√(r̃c̃))/m)=(𝒥(κmid),2√(𝔣(κmid)))Rational iff (r̃,c̃)=λ(U2,W2) (U,W,λ odd, gcd(U,W)=1), m0=(U2+W2)/2≡ 1 (mod 4); least nontrivial m0=5; bridge Ndiag(Gj)=m0,j2Thms. 5.1–5.2, Fig. 4 (Algebraic & descent proofs)
§6. Address Sheet & Odd-Square FactorizationSite pN=(-r/(m+1),1/(m+1)), odd square (2n+1)2=8Tn+1, balanced split {n,n+1}, stack SnClosure, hull & 1/√(k+1) lift; 8-sector E(2,ℤ) split; finite 𝒮LΠ△,L(-1/8)=1; PGF G(τ); Gamma PL(D); atomic PMFs → 2ω,6𝔣(ω); Brouncker & S+Props. 6.1–6.5, Lem. 6.5a, Thm. 6.5, Figs. 5–7, 9, 12 (Algebraic & analytic proofs)
§7. Open QuestionsOdd chamber m=2M+1, center Nc=(m2+1)/2; curve E: C2=(A-2)(A2-200)Mirror-prime pairs (Nc-k,Nc+k) verified on 1403 ≤ m ≤ 20001; shared a=2 & 45∘ records {1,5,29} pulled back to E0,E1,E2 ∈ E(ℤ)§7.1–§7.3, Fig. 8 (Finite verification; all-index questions open)
§8. Shells, Jordan & Q(N)Shells Rd(k,h)=d2 k-d+h, Jordan Jk(d), residual Q(N)=N/gcd(N,m!), entire EEi(z)Within-shell uniqueness & φ(d2)=dφ(d); three separate Jordan zeta limits; Q(N)∈{1}∪ℙ, exceptional prime powers {p2}∪{27}, density log 2; entire lift G±(z)Thms. 8.1–8.3, Prop. 8.4, Fig. 10 (Local proofs + [13]; §8.2 table finite; 2-pt occupancy open)
App. A. Carrier BudgetsInteger budget M ≥ 1, excluded target p*=((A2-C2)/D0,(2AC)/D0), leg chart κleg=r/(m+1)All-budget 45∘ Pell optimum; Bézout tail M ≥ Mcut with M E*(M)2 → 2/D0 & (3,4) budget-17 exception (3,5); two-chart m0,j-1 advanceThms. A.1–A.3, Fig. A.1 (Algebraic & descent proofs)
App. B. Stack ΩnΩn=⋃j=0n [-j-1/2,j+1/2]2×[j,j+1], viewing vector v=(A,C,D)Exact union area law; hiding cone max(| A|,| C|)≤| D|; extrema [(n+1)2, (2n+1)√((2n+1)2+2)]; mean (1+√2)n2+(4+χ)n+C*+εn; sharp τ=1/2Thms. B.1–B.3, Fig. B.1 (Geometric & analytic proofs)

1. The Address, Its Chambers, and Phase Centering Conventions

1.1 Enumerating the Positive Quadrant

Let the triangular polynomial be defined for every complex argument z ∈ ℂ, and write its integer values with subscripts:

T(z) = (z(z+1))/2, Tn = T(n) (n ∈ ℤ), N(r,c) = Tr+c-2 + r (r,c ∈ ℤ≥ 1). (1.1)

For each fixed anti-diagonal sum r+c, the term r counts the sites along that anti-diagonal. Define the chamber index m, the row offset V, and the centered odd legs (r̃, c̃) by

m = r + c - 1, V = r - c, r̃ = 2r - 1 = m + V, c̃ = 2c - 1 = m - V. (1.2)

For each fixed chamber m ≥ 1, the rows r = 1, …, m (with c = m + 1 - r) carry the consecutive integer addresses Tm-1 + 1, …, Tm. The first chamber (m=1) contains 1; the second (m=2) contains 2 and 3; the third (m=3) contains 4, 5, and 6. These disjoint consecutive intervals partition ℤ≥ 1 without gaps. We distinguish the integer address N written in a cell, its row and column (r,c), and its chamber m; none is a substitute for the others.

Proposition 1.1 (Exact Chamber Wedge and Parity Lattice). A positive integer cell (r,c) ∈ ℤ≥ 12 of address N = N(r,c) satisfies

2N = m2 + V + 1, r = (m + V + 1)/2, c = (m - V + 1)/2. (1.3)

Conversely, an integer pair (m,V) ∈ ℤ2 represents a positive integer cell (r,c) ∈ ℤ≥ 12 if and only if

m ≥ 1, | V | ≤ m - 1, m + V ≡ 1 (mod 2). (1.4)

Proof. In chamber m, N = Tm-1 + r = (m(m-1))/2 + r. Multiplying by 2 and substituting 2r = m + V + 1 yields 2N = m2 - m + (m + V + 1) = m2 + V + 1, and solving r+c=m+1, r-c=V gives the coordinate recovery formulas in (1.3). Conversely, given (m,V) ∈ ℤ2, the parity condition m+V ≡ 1 (mod 2) makes r and c integers, and the inequalities m ≥ 1, -m+1 ≤ V ≤ m-1 make r ≥ 1 and c ≥ 1. Substituting (r,c) into (1.1) recovers N. Thus (1.4) is both necessary and sufficient. In particular, any algebraic curve in the (m,V) plane represents actual positive grid cells only where it meets both the parity lattice and the chamber wedge (1.4). □

1.2 Two Exact Inverse Formulas and the Origin Convention

Proposition 1.2 (Inverse Address Recovery). For each positive integer address N ∈ ℤ≥ 1, its unique chamber m and cell coordinates (r,c) are given by

Tm-1 < N ≤ Tm, m = ⌈ (√(8N+1) - 1)/2 ⌉, r = N - Tm-1, c = m + 1 - r. (1.5)

Equivalently, the chamber m is recovered by nearest-integer rounding of √(2N):

m - 1/2 < √(2N) < m + 1/2, m = ⌊ √(2N) + 1/2 ⌋. (1.6)

Proof. Because the chamber intervals (Tm-1, Tm] ∩ ℤ partition ℤ≥ 1, the unique chamber m satisfies ((m-1)m)/2 < N ≤ (m(m+1))/2. Completing the square 8Tm + 1 = (2m+1)2 gives 2m - 1 < √(8N+1) ≤ 2m + 1, which is equivalent to (1.5). Within chamber m, integrality of N strengthens the lower inequality N ≥ Tm-1 + 1 to

2N ≥ m2 - m + 2 = (m - 1/2)2 + 7/4 > (m - 1/2)2,

while the upper inequality N ≤ Tm gives

2N ≤ m2 + m = (m + 1/2)2 - 1/4 < (m + 1/2)2.

Taking positive square roots yields the strict inequalities in (1.6) and the floor formula m = ⌊ √(2N) + 1/2 ⌋. A nearest-integer tie √(2N) = m ± 1/2 can never occur for an integer N, because 2N is an integer whereas (m ± 1/2)2 = m2 ± m + 1/4 ∉ ℤ. □

As a worked address, N = 25 lies between T6 = 21 and T7 = 28. By (1.5)–(1.6), it has chamber m = 7, cell (r,c) = (4,4), and offset V = 0. The identity 2 · 25 = 72 + 1 describes its chamber position (m=7, V=0), whereas 25 = 52 describes its square address (Q=5). Keeping the chamber index m and the address root Q separate is essential in §§4–5.

Figure 1: Positive Quadrant Enumeration, Admissible Chamber Wedge, and Inverse Bounds
Figure 1 (Positive quadrant address, chamber wedge, and inverse recovery): A. Positive quadrant cells (r,c)∈ℤ≥ 12 colored by triangular chamber bands m=r+c-1≤ 7 with cumulative addresses N(r,c)=Tm-1+r and worked diagonal cell N=25 at (r,c)=(4,4), m=7, V=0. B. Admissible chamber-offset lattice (m,V) in Proposition 1.1 (2N = m^2+V+1, |V|≤ m-1, m+V≡ 1 (mod 2)), bounded by the first-column edge V=m-1 (c=1, yellow) and the first-row edge V=1-m (r=1, orange), showing the odd-m diagonal center V=0 (purple) and the even-m superdiagonal center V=-1 (teal). C. Nearest-integer root inverse m=⌊√(2N)+1/2⌋ in Proposition 1.2, showing √(2N)∈(m-1/2,m+1/2) strictly inside each chamber band with no half-integer ties.

Origin Convention. The positive grid has its corner cell at (r,c) = (1,1), holding N(1,1) = 1. Translating the drawing's corner to a Cartesian origin (0,0) replaces (r,c) by (r',c') = (r-1, c-1) ∈ ℤ≥ 02:

r' = r - 1, c' = c - 1, N(r'+1, c'+1) = Tr'+c' + r' + 1 (r',c' ≥ 0). (1.7)

This coordinate translation changes where the axes sit in a drawing without creating a cell of address 0 in the positive Grid. Extending T(z) to z=0 or negative arguments (such as the vertex z = -1/2 where T(-1/2) = -1/8) is a separately declared polynomial operation.

1.3 Relation to Rational Triangular Coordinates and Centering Conventions

The cumulative address (1.1) shares its indexing with the rational triangular-fractional coordinate studied in [4, 11], where chamber m ≥ 1 carries the m entries

a(m,r) = m - 1 + r/m (1 ≤ r ≤ m), (1.8)

and the rank of the entry a(m,r) in this enumeration is exactly the address N(r,c) = Tm-1 + r with c = m + 1 - r [7]. Two structural distinctions are maintained throughout this paper:

  1. Native Row Centering vs. Cell-Midpoint Centering: Using the affine centering map 𝒥(w) = 2w - 1 and quadratic fold 𝔣(w) = w(1-w) = (1 - 𝒥(w)2)/4 [5], centering a chamber by its native row phase κrow = r/m yields the integer displacement
jrow = 2r - m = V + 1, 𝒥(κrow) = (V+1)/m. (1.9)

By contrast, cell-transpose symmetry (r,c) ↔ (c,r) centers a chamber by its cell-midpoint phase κmid = (r - 1/2)/m = r̃/(2m), for which

𝒥(κmid) = (r̃ - c̃)/(2m) = V/m, 4𝔣(κmid) = (r̃c̃)/m2 = 1 - (V/m)2 = 1 - 4(κmid - 1/2)2. (1.10)

Consequently, the native half-row center 𝒥(κrow) = 0 is the superdiagonal cell V = -1 (occurring only in even chambers m = 2a, at r=a, c=a+1, N=2a2), whereas the midpoint fixed center 𝒥(κmid) = 0 is the diagonal cell V = 0 (occurring only in odd chambers m = 2k-1, at r=c=k, N=k2+(k-1)2). This unit shift jrow = V + 1 explains why the two centering conventions select opposite chamber parities.

  1. Signed Binomial Arrays vs. Positive Quadrant Addresses: Companion studies of alternating binomial coefficient arrays (-1)jC(k, j) and their row-reversal involutions [4] act on finite polynomial coefficient tables rather than on the positive integer cell addresses N(r,c) and fixed-chamber spatial reflections proved in §3 below.

2. Straight Cell Paths, Five Rails, and Parity

2.1 Every Affine Cell Path Has a Quadratic Address

Proposition 2.1 (Affine Path Polynomial). Fix a baseline cell (r0, c0) ∈ ℤ≥ 12, integer increments (Δr, Δc) ∈ ℤ2, and an integer parameter k ∈ ℤ, retaining precisely the indices k for which r(k) = r0 + Δr k ≥ 1 and c(k) = c0 + Δc k ≥ 1. Set L = r0 + c0 - 2 and D = Δr + Δc. Then

N(r0 + Δr k, c0 + Δc k) = D2/2 k2 + ((D(2L+1))/2 + Δr)k + TL + r0. (2.1)

Proof. Along the path, the triangular index is r(k) + c(k) - 2 = L + Dk. Expanding T(L + Dk) = TL + (D(2L+1))/2 k + D2/2 k2 and adding r(k) = r0 + Δr k yields (2.1). □

Every affine cell path therefore traces an address polynomial in k of degree at most 2. When D = Δr + Δc = 0, the path remains inside a single chamber m = L + 1 and its quadratic term vanishes, leaving a linear function of k.

2.2 The Five Distinguished Rails

Five affine rays carry most of the structure below. We name them by position; the colors are the keys used in the figures (and in the original notebook). Each is indexed by an integer parameter k ≥ 1:

Rail (figure color)Cell (r,c)Chamber m and Offset VExact Address N(r,c)Line in (m,r̃)
First-row edge (orange)(1, k)m = k, V = 1 - kTk-1 + 1r̃ = 1
First-column edge (yellow)(k, 1)m = k, V = k - 1Tkr̃ = 2m - 1
Diagonal (purple)(k, k)m = 2k - 1, V = 0k2 + (k-1)2 = 4Tk-1 + 1r̃ = m
Superdiagonal (teal)(k, k+1)m = 2k, V = -12k2 = m2/2r̃ = m - 1
Subdiagonal (pink)(k+1, k)m = 2k, V = +12k2 + 1 = (m2+2)/2r̃ = m + 1

In chamber and odd-row coordinates (m, r̃), all five families are straight lines:

r̃ = 1, r̃ = 2m - 1, r̃ = m, r̃ = m - 1, r̃ = m + 1. (2.2)

As functions of the chamber m, however, the three interior rails carry quadratic addresses:

N(k,k) = (m2+1)/2 = 4Tk-1 + 1, N(k,k+1) = m2/2, N(k+1,k) = (m2+2)/2. (2.3)

Thus a statement that the diagonal is "linear" requires naming the horizontal axis: in (m, r̃) it is the line r̃ = m; against the triangular input Tk-1 its address is 4Tk-1+1; against the chamber m its address is the quadratic (m2+1)/2.

Figure 2: The Five Rails on the Notebook Grid, Straight Rays in (m, r-tilde), and Fixed-V Parity 4-Cycle
Figure 2 (The five rails, straight odd-row rays vs. quadratic chamber addresses, and fixed-V parity 4-cycle): A. The 8× 8 positive quadrant (r,c) showing the five rails on the Grid: the first-row edge (1,k) (orange, N=Tk-1+1), the first-column edge (k,1) (yellow, N=Tk), the diagonal (k,k) (purple, V=0, N=4Tk-1+1), the superdiagonal (k,k+1) (teal, V=-1, N=2k2), and the subdiagonal (k+1,k) (pink, V=+1, N=2k2+1). B. In chamber and odd-row coordinates (m,r̃), all five rails are straight lines (2.2), whereas against the chamber m the interior rails carry the quadratic addresses (2.3). C. Constant address parity N(r,r-V)≡ TV+1 (mod 2) along each fixed-V rail (Proposition 2.2), repeating with period 4 in the pattern (1,1,0,0) from V=0,1,2,3 and inducing the one-dimensional four-state Gray cycle of Remark 2.3.

2.3 Fixed-Offset Rails, Parity Period Four, and Odd Powers of Two

Proposition 2.2 (Constant Parity of Fixed-V Rails). Fix an integer offset V ∈ ℤ, and let r ∈ ℤ vary over r ≥ max(1, V+1) with c = r - V ≥ 1. Then

N(r, r - V) = 2r2 - 2(V+1)r + TV+1, N(r, r - V) ≡ TV+1 (mod 2). (2.4)

Consequently, every fixed-V rail has constant address parity throughout the positive quadrant, and as V increases by 1 that rail parity repeats with period 4 in the pattern (1, 1, 0, 0) starting at V = 0, 1, 2, 3.

Proof. In (1.1), r + c - 2 = 2r - V - 2, so

N(r, r - V) = ((2r - V - 2)(2r - V - 1))/2 + r = 2r2 - 2(V+1)r + ((V+1)(V+2))/2 = 2r2 - 2(V+1)r + TV+1.

The polynomial identity T(-V-2) = T(V+1) holds for all V ∈ ℤ. Since 2r2 - 2(V+1)r is even for every integer r, N(r, r-V) ≡ TV+1 (mod 2). Moreover, TV+5 - TV+1 = 4V + 14 ≡ 0 (mod 2), giving period 4 with values T1=1, T2=3≡ 1, T3=6≡ 0, T4=10≡ 0 (mod 2) at V = 0, 1, 2, 3. □

Remark 2.3 (Connection to the One-Dimensional Gray Cycle and the Modular Coefficient Rows). In one dimension, consecutive triangular parities form the adjacent pair (Tn-1 mod 2, Tn mod 2), which traverses the four-state Gray cycle 00 → 01 → 11 → 10 → 00 by the quarter-turn F(a,b) = (b, 1-a), while index reflection n ↦ -1-n acts by E(a,b) = (1-a, b), together generating D4. That one-dimensional index-cycle reflection on modular rising-factorial rows Πn(y) = ∏j=1n (1+jy) mod (B, yK+1) [7] acts across consecutive triangular indices n and is distinct from the fixed-chamber spatial reflections in §3 below.

Finally, the odd powers of two lie exactly on the even-parity superdiagonal V = -1:

22j+1 = 2(2j)2 ⇒ (r, c) = (2j, 2j + 1), m = 2j+1, V = -1 (j ∈ ℤ≥ 0). (2.5)

This sparse exponential subsequence lies inside the quadratic rail N = 2k2 (k = 2j) and must not be confused with the piecewise parabolic branches traced by even powers of two (4j = (2j)2) in §4.

3. Three Reflections and Their Surviving Integer Domains

3.1 Three Ways to Center and Reflect a Chamber

Hold an integer chamber m ≥ 1 fixed. Within chamber m, the native row phase κrow, the cell-midpoint phase κmid, and the square-root phase κroot (with local root displacement ε) are defined by

κrow = r/m, κmid = (r - 1/2)/m = r̃/(2m), ε = √(2N) - m, κroot = ε + 1/2 = √(2N) - m + 1/2. (3.1)

On a positive cell (r,c) in chamber m, κrow = r/m ∈ (0, 1] is a right-endpoint step phase (attaining 1 on the first-column edge r=m, c=1, whereas the address-sheet coordinate -p1 = κleg = r/(m+1) = r/(r+c) ∈ [1/(m+1), m/(m+1)] ⊂ (0,1) of §6.1 divides by m+1 = r+c so both coordinates remain strictly positive inside the open unit interval), while both κmid = (r-1/2)/m and κroot = √(2N)-m+1/2 are interior cell phases lying strictly in (0, 1) by Proposition 1.2 (-1/2 < ε < 1/2). Consequently, at the first-column endpoint r=m, c=1 (N=Tm), the involution κrow ↦ 1 - κrow sends κrow=1 to 0, mapping Rrow(Tm) = m2 - Tm = Tm-1 out of chamber m into the preceding chamber m-1 (or to 0 at m=1). Although all three coordinates represent positions across chamber m, applying the scalar involution κ ↦ 1 - κ while holding m fixed induces three different address maps.

Theorem 3.1 (Three Reflection Maps and Their Domains). Fix an integer chamber m ≥ 1. Reflecting κrow, κmid, and κroot by κ ↦ 1 - κ induces the following address transformations on chamber m:

Selected PhaseInduced Address MapSame-Chamber Positive Cell Domain & ActionFixed Positive Cells in Chamber m
κrow ↦ 1 - κrowRrow(N) = m2 - N(r,c) ↦ (c-1, r+1); preserves chamber m iff c > 1 (1 ≤ r ≤ m-1, empty for m=1)Even m=2a: superdiagonal cell (a, a+1), V=-1, N = m2/2 = 2a2
κmid ↦ 1 - κmidRmid(N) = m2 + 1 - N(r,c) ↦ (c, r); permutes every positive cell of chamber mOdd m=2k-1: diagonal cell (k,k), V=0, N = (m2+1)/2
κroot ↦ 1 - κrootRroot(N) = ((2m - √(2N))2)/2Continuous involution on √(2N) ∈ (m-1/2, m+1/2); integer image iff fixed superdiagonal cellEven m=2a: superdiagonal cell (a, a+1), V=-1, N = m2/2 = 2a2

Moreover, in terms of the offset V = 2mε + ε2 - 1, the transformed offsets and the exact nonlinear difference between Rroot and Rrow are

V'row = -V - 2, V'mid = -V, V'root = -V - 2 + 2ε2, Rroot(N) - Rrow(N) = ε2. (3.2)

Proof. Under κrow ↦ 1 - κrow, the row r = mκrow maps to r' = m - r = c - 1, giving address Tm-1 + (m - r) = m2 - N and column c' = m + 1 - r' = r + 1. This image lies in chamber m with positive coordinates iff 1 ≤ r ≤ m-1 (c > 1). At the first-column endpoint r = m, c = 1 (N = Tm), Rrow(Tm) = m2 - Tm = Tm-1, which is a positive address in the preceding chamber m-1 when m ≥ 2 and equals 0 (outside the positive grid) when m = 1.

Under κmid ↦ 1 - κmid, r - 1/2 maps to m - (r - 1/2), so r' = m + 1 - r = c, c' = r, and Rmid(N) = Tm-1 + (m + 1 - r) = m2 + 1 - N, which is a bijection on {1, …, m}.

Under κroot ↦ 1 - κroot, ε = κroot - 1/2 maps to -ε, so √(2N) = m + ε maps to m - ε = 2m - √(2N) ∈ (m - 1/2, m + 1/2), giving Rroot(N) = (2m - √(2N))2/2. Expanding 2N = (m+ε)2 = m2 + 2mε + ε2 and comparing with 2N = m2 + V + 1 gives V + 1 = 2mε + ε2, so

Rroot(N) = ((m - ε)2)/2 = (m2 - 2mε + ε2)/2 = (m2 - (V + 1 - ε2) + ε2)/2 = m2 - N + ε2 = Rrow(N) + ε2,

which establishes all relations in (3.2). □

Single-Cell Distinguishing Witness. For address N = 5 in chamber m = 3 (cell (2,2), V = 0), the three reflections give

Rrow(5) = 4, Rmid(5) = 5, Rroot(5) = 23 - 6√10 ≈ 4.026334, (3.3)

with 4 < 23 - 6√10 < 5. This single cell rules out any identification of the three reflections.

3.2 Integer-Image Rigidity and Root–Midpoint Discrepancy

Proposition 3.2 (Integer-Image Rigidity of Root Reflection). Among positive integer cells, Rroot(N) is an integer if and only if the cell lies on the superdiagonal:

m = 2a, N = 2a2, (r,c) = (a, a+1), ε = 0 (a ∈ ℤ≥ 1), (3.4)

and every such cell is fixed by Rroot.

Proof. Expanding Rroot(N) = 2m2 + N - 2m√(2N) shows that if Rroot(N) ∈ ℤ, then √(2N) ∈ ℚ. Because 2N ∈ ℤ, a rational square root of 2N must be an integer. By Proposition 1.2, m - 1/2 < √(2N) < m + 1/2, so the only integer in this open interval is m. Hence √(2N) = m, so 2N = m2 and ε = 0. Integrality of N = m2/2 forces m = 2a even, giving N = 2a2 and (r,c) = (a, a+1). Conversely, ε = 0 gives Rroot(2a2) = 2a2. □

Proposition 3.2 proves without any unbounded search that Rroot never exchanges distinct integer Pell families: off the fixed superdiagonal, Rroot(N) is always irrational. Even the affine map Rrow need not preserve squares (in chamber m=7, Rrow(25) = 49 - 25 = 24).

Proposition 3.3 (Root–Midpoint Phase Discrepancy Bounds). On every positive integer cell (r,c) in chamber m ≥ 1,

m(κroot - κmid) = (1 - ε2)/2, 3/(8m) < κroot - κmid ≤ 1/(2m). (3.5)

Equality at the upper bound 1/(2m) holds if and only if the cell is a superdiagonal cell (3.4), and the lower bound 3/(8m) is strict on every positive cell.

Proof. Since V = 2mε + ε2 - 1, we have

κmid = 1/2 + V/(2m) = 1/2 + ε + (ε2 - 1)/(2m) = κroot - (1 - ε2)/(2m),

which gives m(κroot - κmid) = (1 - ε2)/2. Proposition 1.2 gives 0 ≤ ε2 < 1/4, which yields 3/(8m) < κroot - κmid ≤ 1/(2m), with upper equality iff ε = 0. □

Remark 3.4 (Anti-Diagonal Pair Correspondence and Even-Chamber Midpoint Reflection). For any integer p ≥ 2 and step 1 ≤ k ≤ p - 1 (no primality condition on p is required), the anti-diagonal pair (r, c) = (k, p - k) studied in the companion sum–difference geometry [15] lies on the Notebook Grid at

(r, c) = (k, p - k), m = p - 1, V = 2k - p, N(k, p - k) = Tp-2 + k. (3.6)

Exchanging the pair (k, p-k) ↦ (p-k, k) sends N = Tp-2 + k to N' = Tp-2 + (p - k). Since 2Tp-2 + p = (p-2)(p-1) + p = (p-1)2 + 1 = m2 + 1, we have N' = m2 + 1 - N = Rmid(N): coordinate exchange on (k, p-k) is identically the midpoint reflection Rmid on chamber m = p - 1. When p is odd, the chamber m = p - 1 is even and every offset V = 2k - p is odd, so the diagonal center V = 0 is absent from chamber m and the central adjacent pair ((p-1)/2, (p+1)/2), ((p+1)/2, (p-1)/2) occupies the superdiagonal and subdiagonal offsets V = -1 and V = +1. Equation (3.6) records the discrete pair correspondence; the companion polar scalar of [15], the exchange-commuting 2× 2 Gram family Ba = I + aE (characterized in §6.6, Remark 6.6b(1)), and Gaussian-integer lattice squares live in distinct ambient spaces and are not identified as a single point set, nor is this even chamber m = p - 1 (for odd p) identified with the odd-chamber (m = 2M + 1) midpoint prime-pair setting of §7.1.

Figure 3: Three Chamber Reflections, Single-Cell N=5 Witness, and Phase Discrepancy Band
Figure 3 (Three chamber reflections, distinguishing witness N=5, and root–midpoint phase discrepancy): A. In chamber m=3 (addresses {4,5,6}), the diagonal cell N=5 at (2,2) separates all three reflections (3.3): Rrow(5)=4 < Rroot(5)=23-6√10≈ 4.026334 < Rmid(5)=5, while at the first-column endpoint N=6 (c=1) the native map Rrow(6)=3=T2 drops into chamber m=2. B. Fixed-cell parity split (Theorem 3.1) and root-reflection integrality rigidity (Proposition 3.2): Rmid fixes the diagonal V=0 (purple) in odd chambers m=2k-1, whereas Rrow fixes the superdiagonal V=-1, N=2a2 (teal) in even chambers m=2a, which is also the unique positive locus where Rroot(N)∈ℤ. C. Scaled phase discrepancy m(κroot-κmid)=(1-ε2)/2∈(3/8,1/2] across ε=√(2N)-m∈(-1/2,1/2) (Proposition 3.3).

4. Scaled-Square Branches, Pell Seams, and Classical Powers

4.1 Exact Floor-Selected Polynomial Branches

Proposition 4.1 (Floor-Branch Dictionary). Fix positive integers A, Q ∈ ℤ≥ 1, and let P ∈ ℤ≥ 1 be the unique nearest integer to Q√(2A). For each integer k ≥ 1, define

Nk = A(Qk)2, ε = Q√(2A) - P, δ = P2 - 2AQ2, ℓk = ⌊ 1/2 + kε ⌋. (4.1)

Then the chamber mk, offset Vk, and cell coordinates (rk, ck) of Nk are given exactly by

mk = Pk + ℓk, Vk = -δ k2 - 2Pkℓk - ℓk2 - 1, (rk, ck) = ((mk + Vk + 1)/2, (mk - Vk + 1)/2). (4.2)

Proof. By Proposition 1.2, mk = ⌊ √(2Nk) + 1/2 ⌋ = ⌊ k(P + ε) + 1/2 ⌋ = Pk + ℓk because Pk ∈ ℤ. Substituting mk = Pk + ℓk and 2AQ2 = P2 - δ into Vk = 2Nk - mk2 - 1 = 2AQ2 k2 - (Pk + ℓk)2 - 1 yields (4.2). □

On any interval of integer indices k where the integer floor ℓk = L remains constant, mk is linear in k and Vk is quadratic in k. That polynomial formula applies only on the floor-selected domain L - 1/2 ≤ kε < L + 1/2.

Proposition 4.2 (Asymmetric First-Branch Criterion). For the scaled-square path (4.1), a positive integer index k ≥ 1 lies on the initial branch ℓk = 0 if and only if

ℓk = 0 ⇔ -P ≤ δ k < P. (4.3)

Proof. The condition ℓk = 0, i.e. mk = Pk, holds iff TPk-1 < A(Qk)2 ≤ TPk. Multiplying by 2 gives (Pk)2 - Pk < (P2 - δ)k2 ≤ (Pk)2 + Pk. Subtracting P2 k2 and dividing by k > 0 yields -P < -δ k ≤ P, or equivalently -P ≤ δ k < P. □

The strict upper bound δ k < P (for δ > 0) versus the closed lower bound | δ | k ≤ P (for δ < 0) comes directly from the half-open chamber interval Tm-1 < N ≤ Tm. On the first branch ℓk = 0, the offset simplifies to Vk + 1 = -δ k2.

4.2 The Zero-Norm Domain

Proposition 4.3 (Zero-Norm Characterization). For positive integers A, Q ∈ ℤ≥ 1 and P = round(Q√(2A)),

δ = 0 ⇔ A = 2b2 and P = 2bQ for some b ∈ ℤ≥ 1, Nk = 2(bQk)2. (4.4)

In that case ℓk = 0, mk = 2bQk, and Vk = -1 for all k ≥ 1 (the straight superdiagonal).

Proof. If δ = P2 - 2AQ2 = 0, then (P/Q)2 = 2A ∈ ℤ, so P/Q is an integer whose square is even, hence P/Q = 2b for some b ∈ ℤ≥ 1. Thus A = 2b2 and P = 2bQ. The converse is immediate. □

In particular, (A,P,Q) = (8, 4, 1) has δ = 0. Stating that δ = 0 occurs "only for A = 2" requires restricting A to be squarefree.

4.3 The Worked 181/128 Square Branch, First Departure, and Phase-Linearity Shift

Take A = 1, P = 181, Q = 128, for which δ = 1812 - 2 · 1282 = 32761 - 32768 = -7. By Proposition 4.2, the first branch ℓk = 0 holds iff 7k ≤ 181, so its exact integer domain is 1 ≤ k ≤ 25:

Nk = (128k)2, mk = 181k, Vk = 7k2 - 1 (1 ≤ k ≤ 25). (4.5)

The powers of four 47, 48, 49, 410, 411 appear along this branch at k = 1, 2, 4, 8, 16 at cells (94, 88), (195, 168), (418, 307), (948, 501), and (2344, 553). The first integer departure of the full branch Nk = (128k)2 occurs at k = 26, where ℓ26 = 1, m26 = 4707, and (r26, c26) = (13, 4695). By contrast, the next sampled power of four (412) occurs later at k = 32, where ℓ32 = 1, m32 = 5793, V32 = -4418, and (r32, c32) = (688, 5106). Distinguishing the full quadratic sequence from its sparse exponential subsequence prevents confusing k=26 with k=32.

Exact Phase-Linearity Shift on the First Branch. On 1 ≤ k ≤ 25, substituting mk = 181k and Vk = 7k2 - 1 into 2rk = mk + Vk + 1 and r̃k = mk + Vk yields

2κrow(k) = (2rk)/mk = 1 + 7/1812 mk, 2κmid(k) = r̃k/mk = 1 + 7/1812 mk - 1/mk (1 ≤ k ≤ 25). (4.6)

Thus the native row normalization 2κrow is strictly linear in the chamber mk across the entire first branch, whereas the centered odd-row normalization 2κmid = r̃k/mk carries the exact reciprocal offset -1/mk.

Remark 4.3a (Cross-Reference: The (7, 180, 181) Reciprocal Bridge and the Ramanujan–Nagell Solutions). The arithmetic parameters of the 181/128 floor branch in (4.5)–(4.6) align with two structural bridges developed in [4, 16]:

  1. Grid 181/128 Branch Parameters (δ = -7, m1 = 181 = 180 + 1): At the base step k=1, the Ramanujan–Nagell relation 1812 - 2· 1282 = -7 (1812 + 7 = 215) places N1 = 1282 = 16384 = 47 in chamber m1 = 181 = 180 + 1 immediately above T180 = (180· 181)/2 = 16290, at (r1, c1) = (94, 88) and offset V1 = 6. In (4.6), the phase-linearity slope is 7/181 per step k (7/1812 per chamber unit mk), governed by |δ| = 7, while 2κmid(k) - 2κrow(k) = -1/mk = -1/(181k). Because P = 181 = 180 + 1 is prime, all φ(181) = 180 proper interior row phases κrow = r/181 (1 ≤ r ≤ 180) in chamber m1 = 181 are reduced fractions, and at n = 180 the rising-factorial interior gcd of Theorem 8.3(5) evaluates to g180 = 181.
  2. Reciprocal, Catalan, and Degree-to-Radian π-Phase Bridge at n = 7 [4]: At the third Mersenne index n = 7 = 23 - 1, the triangular number is T7 = 28 (2T7 = 56 = 7· 8). Because n+1 = 8 is composite, T7/7! = 1/180 with reduced denominator b7 = (7!)/T7 = 5040/28 = 180 = (3!)2 C3 (the j=3 odd-slice Catalan denominator in T2j+1/(2j+1)! = 1/[(j!)2 Cj]), giving
1/7 →÷ 4 1/T7 = 1/28 →÷ 180 1/(7!) = 1/5040, (1/T7)∘ = 1/7·1/4·π/180 rad = π/(180 T7) rad = π/(7!) rad = π/5040 rad,

while across dyadic halvings i ≥ 0, 1/(7· 2i) = 100-i(Ai + 1/7) with Ai = (50i - 1)/7 ∈ ℤ (50 ≡ 1 (mod 7)).

  1. The Value u = 1/8 and the Rational π Bracket [4]: In the genus-zero triangular cosine product Π△(u) = ∏n≥ 1(1 - u/Tn) = -(cos(π/2 √(1+8u)))/(2π u) of [4] (Remark 6.6a), the value u = 1/8 gives the first-factor ratio u/(T1 - u) = 1/(8T1 - 1) = 1/7 and rational bounds 1/8 < π - 3 < 1/7 ⇔ 25/8 < π < 22/7, with 11/4 - (7π)/8 = 7/8 (22/7 - π). Scope: On the 181/128 floor branch, (7, 180, 181) arises from δ = 1812 - 2· 1282 = -7 and prime chamber m1 = 181, whereas at n = 7 it arises from 7!/T7 = 180 and degree-to-radian conversion π/180; spatial chamber phase κrow = r/m and angular radian conversion are kept distinct.
  2. The Complete Ramanujan–Nagell Solution Set [1, 16]: Writing the odd chamber index as m1 = 2n + 1 so m12 = 8Tn + 1, the binary-square condition m12 + 7 = 8(Tn + 1) = 2μ is Nagell's Mersenne-triangular equation M2 + 7 = 2μ ⇔ Tn = 2μ-3 - 1 [1], whose complete nonnegative integer solution set is
(n, M=2n+1, Tn, μ) ∈ {(0, 1, 0, 3), (1, 3, 1, 4), (2, 5, 3, 5), (5, 11, 15, 7), (90, 181, 4095, 15)}.
  • The largest solution n = 90 (2n = 180, M = 181, T90 = 4095 = 212-1, μ = 15): Here 1812 + 7 = 32768 = 215, so on 1 ≤ k ≤ 25, 2N(rk,ck) = mk2 + Vk + 1 = (1812 + 7)k2 = 215k2 ⇔ N(rk,ck) = (128k)2, and 1 + 7/1812 = 215/1812 = ((128√2)/181)2.
  • The solution n = 5 (M = 11, T5 = 15 = 24-1, μ = 7): The solution M = 11 (112 + 7 = 27) yields the penultimate V=6 binary square N = 64 = 43 in chamber m = 11 (Eq. 4.13), the unique nontrivial iterated triangular-factorial identity TT5 = T15 = 120 = 5! (since TTn = n! fails for n ∈ {1,2,3,4} and for n ≥ 6 Stirling's bound gives n! > 1/8 n4 ≥ TTn for all n ≥ 7 with T21=231 ≠ 720 = 6!), and the largest pure power of 3 on the subdiagonal at k = 11, where 2(11)2 + 1 = 243 = 35 in chamber m = 22 (Eq. 4.14).
  • The reciprocal-π product: In Π△(-1/8) = ∏n≥ 1((2n+1)2)/(8Tn) = 4/π (§6.6, Theorem 6.5), the four positive solutions n ∈ {1, 2, 5, 90} (M ∈ {3, 5, 11, 181}) are the unique stages where numerator (2n+1)2 = 2μ - 7 and denominator 8Tn = 2μ - 8 are simultaneously shifted powers of two.

4.4 Norm-One Seams, k = 2P Returns, and the Square-Triangular Converse

Proposition 4.4 (Norm-One Launch and Triangular Seams). Assume δ = P2 - 2AQ2 ∈ {-1, +1}. Then the per-step root drift is

ε = -δ/(P + Q√(2A)), 1/(| ε |) = P + Q√(2A). (4.7)

At k = 1, the path launches at m1 = P, V1 = -δ - 1. At k = P, the path lands on an exact triangular seam:

NP = A(PQ)2 = TP2 if δ = -1 (in chamber mP = P2, ℓP = 0), ; TP2 - 1 if δ = +1 (in chamber mP = P2 - 1, ℓP = -1). (4.8)

Proof. Multiplying (P - Q√(2A))(P + Q√(2A)) = δ gives (4.7). At k = P, 2NP = P2(2AQ2) = P2(P2 - δ), so NP = (P2(P2+1))/2 = TP2 when δ = -1 and NP = ((P2-1)P2)/2 = TP2-1 when δ = +1. By Proposition 4.2, k=P satisfies -P ≤ δ P < P for δ = -1 (ℓP = 0) and violates it by one chamber for δ = +1 (ℓP = -1). □

Proposition 4.5 (Exact Algebraic Return at k = 2P). Under δ = P2 - 2AQ2 = ± 1, set M = P2 + 2AQ2. At k = 2P,

m2P = M = P2 + 2AQ2, V2P = -2, (r2P, c2P) = ((M-1)/2, (M+3)/2), (4.9)

corresponding to the quadratic-unit squaring identity (P + Q√(2A))2 = M + 2PQ√(2A).

Proof. Expanding M2 - 2A(2PQ)2 = (P2 - 2AQ2)2 = δ2 = 1 gives 2N2P = 2A(2PQ)2 = M2 - 1. Since M = 2P2 - δ ≥ 3 is odd, (M - 1/2)2 < M2 - 1 < (M + 1/2)2, so Proposition 1.2 gives m2P = M and Proposition 1.1 gives V2P = 2N2P - M2 - 1 = -2. □

Proposition 4.6 (Complete Square-Triangular Converse). Every positive square-triangular number Tn = L2 (n, L ∈ ℤ≥ 1) arises as an A = 1 norm-one seam (4.8): there exist unique positive integers P, Q ∈ ℤ≥ 1 with P2 - 2Q2 = δ ∈ {-1, +1} such that L = PQ and n = P2 (δ = -1) or n = P2 - 1 (δ = +1).

Proof. In n(n+1) = 2L2, gcd(n, n+1) = 1, so the odd member of {n, n+1} is a square P2 and the even member is twice a square 2Q2. Then L = PQ and P2 - 2Q2 = ∓ 1. Since | Q√2 - P | = 1/(P + Q√2) < 1/2, P = round(Q√2), placing (P,Q) in the A=1 branch family (4.8). □

4.5 Square Addresses on the Diagonal and the Subdiagonal

Proposition 4.7 (Diagonal-Square Pell Equivalence). For any square address N = Q2 (Q ∈ ℤ≥ 1) in chamber m = ⌊ Q√2 + 1/2 ⌋, its offset is V = 2Q2 - m2 - 1. Consequently:

  1. N = Q2 lies on the diagonal (V = 0) if and only if
m2 - 2Q2 = -1, m = 2k - 1, Q2 = k2 + (k-1)2 (k = r = c). (4.10)
  1. N = Q2 lies on the subdiagonal (V = +1) if and only if
m2 - 2Q2 = -2, m = 2k, Q2 = 2k2 + 1 (r = k+1, c = k), (4.11)

whose first positive instances are (m, Q) = (4, 3), (24, 17), (140, 99) with square addresses N = 9, 289, 9801.

Proposition 4.8 (Complete Negative-Pell Descent for Diagonal Squares). Every positive integer solution of m2 - 2Q2 = -1 is generated from the base solution (m0, Q0) = (1, 1) by iterating

(m, Q) ↦ (3m + 4Q, 2m + 3Q). (4.12)

Proof. The map (4.12) preserves m2 - 2Q2 = -1 and positivity. Conversely, let (m,Q) be a positive solution with Q > 1. Testing Q ∈ {2,3,4} gives 2Q2 - 1 ∈ {7, 17, 31}, none of which is a square, so Q ≥ 5. Hence 4/3 < m/Q = √(2 - 1/Q2) < 3/2. Applying the inverse integral matrix gives

(m', Q') = (3m - 4Q, 3Q - 2m),

which satisfies m' > 0, Q' > 0, (m')2 - 2(Q')2 = -1, and Q' < Q (since m > Q ⇒ 3Q - 2m < Q). Iterating this strict positive-integer descent must terminate at Q = 1, m = 1. Reversing the steps yields (4.12). □

The first five diagonal-square solutions are (m, Q) = (1,1), (7,5), (41,29), (239,169), (1393,985), with diagonal row indices k = (m+1)/2 ∈ {1, 4, 21, 120, 697} and square addresses N = Q2 ∈ {1, 25, 841, 28561, 970225}.

Figure 3b: Scaled-Square 181/128 Floor Branch, Phase-Linearity Shift, and Pell Seams / Descent
Figure 3b (Scaled-square 181/128 floor branch, phase-linearity shift, and Pell seams/descent): A. Piecewise sample across 1≤ k≤ 34 along Nk=(128k)2 (A=1,P=181,Q=128,δ=-7): the first polynomial branch ℓk=0 (mk=181k, Vk=7k2-1, Eq. (4.6)) holds strictly on 1≤ k≤ 25 (carrying 47,…,411 at k=1,2,4,8,16), departs first at k=26 (ℓ26=1, m26=4707, (13,4695)), and reaches the next power sample 412 at k=32 (m32=5793, V32=-4418). B. On the first branch 1≤ k≤ 25, native row phase 2κrow(k)=1+7/1812 mk is strictly linear in mk, whereas 2κmid(k)=2κrow(k)-1/mk carries the reciprocal offset -1/mk (4.6). C. Logarithmic chamber growth of diagonal squares (m2-2Q2=-1, Proposition 4.8), subdiagonal squares (m2-2Q2=-2), and norm-one first-column-edge triangular seams Tn=(PQ)2 at k=P (Propositions 4.4–4.6).

4.6 Classical Classifications of Powers on the Rails

Several visually prominent power intersections on the rails reduce to classical exponential Diophantine equations:

  • Powers of 2 on the First-Row Edge and Offset V = 6: On the first-row edge (1, c), N = 2j (j ≥ 0) is equivalent to Tc-1 + 1 = 2j, or
(2c - 1)2 + 7 = 2j+3. (4.13)

By Nagell's theorem [1], the positive odd solutions 2c-1 are 1, 3, 5, 11, 181, giving j ∈ {0, 1, 2, 4, 12} and addresses N ∈ {1, 2, 4, 16, 4096}. Similarly, since odd powers of 2 sit on V = -1 (2.5), an even power 4k sits at offset V = 6 iff P2 + 7 = 22k+1 with P = round(2k√2); Nagell's list intersected with the nearest-integer bound leaves P ∈ {11, 181}, so only 64 (m=11) and 16384 (m=181) are powers of 2 with V = 6.

  • Powers of 3 on the Subdiagonal: On the subdiagonal (k+1, k), N = 3j (k, j ≥ 1) is
2k2 + 1 = 3j. (4.14)

By Leu and Li [2], the only positive solutions are (k, j) ∈ {(1,1), (2,2), (11,5)}, giving addresses N ∈ {3, 9, 243}.

  • Powers of 5 on the Diagonal: On the diagonal (k, k), N = 5j (k, j ≥ 1) is
(2k - 1)2 + 1 = 2 · 5j. (4.15)

Cohn's perfect-Pell-power lemma [3] excludes j > 2 for base 5, leaving (k, j) ∈ {(2,1), (4,2)}, diagonal chambers m = 2k - 1 ∈ {3, 7}, and diagonal addresses N = N(k,k) = 5j ∈ {5, 25} (where (2k-1)2 + 1 = 2N ∈ {10, 50}; Cohn's exceptional higher-power solution for an arbitrary base gives (2392+1)/2 = 28561 = 134 at k=120). Moreover, for the fixed base 5 in (4.15), no appeal to Ljunggren's k=4 theorem is required, using two short self-contained prerequisites: (i) in ℤ[√10], any integral unit ω = A + C√10 > 1 of norm A2 - 10C2 = ± 1 has positive integers A, C ≥ 1, where C = 1 gives A = 3 (ε10 = 3 + √10 of norm -1) and C ≥ 2 gives A + C√10 ≥ 1 + 2√10 > ε10, so choosing a ≥ 1 with ε10a ≤ ω < ε10a+1 places ωε10-a ∈ [1, ε10) and forces ω = ε10a (with a odd whenever Norm(ω) = -1); and (ii) for the Pell recurrence matrix G = ([ 2 1 ; 1 0 ]), G3 = ([ 12 5 ; 5 2 ]) ≡ 2I (mod 5), so P3a+i ≡ 2a Pi (mod 5) (i ∈ {0,1,2}) with P0 = 0, P1 = 1, P2 = 2, whence 5 | Pn ⇔ 3 | n. Consequently:

  • When j = 2b + 1 is odd, (4.15) is (2k-1)2 - 10(5b)2 = -1, so Prerequisite (i) gives 5b = C2a+1 = ∑i=0a C(2a+1, 2i+1) 32(a-i) 10i in (3+√10)2a+1 (a ≥ 0). Using C(2a+1, 2i+1) = (2a+1)/(2i+1) C(2a, 2i) and v5(2i+1) ≤ i - 1 < i for all i ≥ 1, the i-th summand has valuation ≥ v5(2a+1) + i - v5(2i+1) ≥ v5(2a+1) + 1 > v5(2a+1) (a sharp one-unit gap, attained at a=1, i=1 where (3+√10)3 = 117 + 37√10 has 37 = 27 + 10 with v5(10)=1). Thus the i=0 term (2a+1)32a strictly dominates the 5-adic valuation, forcing v5(C2a+1) = v5(2a+1) and hence 5b | (2a+1). Since C2a+1 ≥ (2a+1)32a > 2a+1 for all a ≥ 1, this holds only at a=0, b=0, j=1 (k=2, chamber m=3, diagonal address N=5).
  • When j = 2b is even (b ≥ 1), (4.15) is (2k-1)2 - 2(5b)2 = -1, so Proposition 4.8 and Prerequisite (ii) give (2k-1) + 5b√2 = (1+√2)2n+1 with 5b = P2n+1 (P0=0, P1=1, Pi+2=2Pi+1+Pi) and 3 | (2n+1). Writing 2n+1 = 3(2a+1) (a ≥ 0) and expanding (1+√2)3(2a+1) = (7+5√2)2a+1 gives P3(2a+1) = 5∑i=0a C(2a+1, 2i+1) 49a-i 50i. Because 50i = 2i· 52i and v5(2i+1) ≤ i - 1 implies 2i - v5(2i+1) ≥ i + 1 ≥ 2 for all i ≥ 1, each i ≥ 1 summand has valuation ≥ 1 + v5(2a+1) + 2i - v5(2i+1) ≥ 1 + v5(2a+1) + 2, so v5(P3(2a+1)) = 1 + v5(2a+1), forcing 5b-1 | (2a+1). Since P3(2a+1)/5 ≥ (2a+1)49a > 2a+1 for all a ≥ 1, this holds only at a=0, b=1, j=2 (k=4, chamber m=7, diagonal address N=25).

Sources. The completeness statements (4.13) and (4.14) are taken from [1, 2]; (4.15) for base 5 follows both from [3] and from the self-contained 5-adic argument above. In (4.14), factoring (1 + k√(-2))(1 - k√(-2)) = (1 + √(-2))j(1 - √(-2))j in the Euclidean domain ℤ[√(-2)] reduces 2k2 + 1 = 3j to Re((1+√(-2))j) = ± 1 (j ∈ {1, 2, 5}).

5. Rational Carriers, the Route to 5, and the Diagonal Bridge

5.1 The Midpoint Fold Circle Map and Unique Rational Carrier Parameterization

In [5, §2.1], the analytic fold operator 𝔣(w) = w(1-w) = (1 - 𝒥(w)2)/4 is applied to the complex argument s of ξ(s) to produce q = 𝔣(s). On the discrete address grid, applying the same operators (𝒥, 𝔣) to the cell-midpoint phase κmid = r̃/(2m) defines the centered carrier map on the open upper unit semicircle:

𝒞(r,c) = (𝒥(κmid), 2√(𝔣(κmid))) = ((r̃ - c̃)/(r̃ + c̃), (2√(r̃c̃))/(r̃ + c̃)) = (V/m, (√(r̃c̃))/m). (5.1)

Because r, c ≥ 1, both odd legs r̃ = 2r-1 and c̃ = 2c-1 are positive, so the second coordinate of 𝒞(r,c) is strictly positive and the boundary endpoints (± 1, 0) are excluded. This is a specified algebraic map on positive grid cells, distinct from polar-shell embeddings Z = √d e2π i h/d and from κroot.

Theorem 5.1 (Unique Rational Carrier Classification). For a positive integer cell (r,c) ∈ ℤ≥ 12, the circle point 𝒞(r,c) in (5.1) belongs to ℚ2 if and only if there exist unique positive odd integers U, W, λ ∈ ℤ≥ 1 with gcd(U, W) = 1 such that

(r̃, c̃) = λ(U2, W2). (5.2)

Equivalently, the cell (r,c), its chamber m, its primitive first chamber m0, its carrier point p(U,W) = 𝒞(r,c), and its midpoint phase κmid are given by

(r, c) = ((λ U2 + 1)/2, (λ W2 + 1)/2), m = λ m0, m0 = (U2 + W2)/2, (5.3)
p(U,W) = ((U2 - W2)/(2m0), (UW)/m0), κmid = U2/(U2 + W2). (5.4)

Here λ = gcd(r̃, c̃) is the odd ray-repeat scale, the reduced denominator of κmid is 2m0, and the primitive first chamber satisfies

m0 ≡ 1 (mod 4), ((U2 - W2)/2, UW, m0) is a primitive Pythagorean triple (5.5)

(with signed first leg, reducing to the degenerate triple (0, 1, 1) at (U,W) = (1,1)).

Proof. Since V/m ∈ ℚ for every cell, 𝒞(r,c) ∈ ℚ2 iff √(r̃c̃) ∈ ℚ, which holds iff the odd integer r̃c̃ is a perfect square. Let λ = gcd(r̃, c̃), which is a positive odd integer. Then r̃/λ and c̃/λ are coprime positive odd integers whose product (r̃c̃)/λ2 is a rational square, hence an integer square. By coprimality, r̃/λ = U2 and c̃/λ = W2 for unique positive odd coprime integers U, W. Conversely, (5.2) makes √(r̃c̃) = λ UW ∈ ℤ. For a fixed oriented point p = (p1, p2), p2/(1 - p1) = U/W determines the coprime positive pair (U,W) uniquely, while λ = m/m0 runs over all positive odd integers. Since any odd square is ≡ 1 (mod 8), U2 + W2 ≡ 2 (mod 8), so m0 = (U2+W2)/2 ≡ 1 (mod 4). □

5.2 First Chambers 5, 13, 17, the Route to 5, and Multi-Ray Chambers

The only primitive diagonal carrier has U = W = 1, giving p(1,1) = (0,1) at first chamber m0 = 1, cell (1,1). For any nontrivial primitive carrier (U ≠ W), the smallest unordered pair of distinct positive odd coprime integers is {3, 1}. Therefore, the least nontrivial rational carrier chamber in the centered chart is m0 = (32 + 12)/2 = 5, realized at cell (5,1) with p(3,1) = (4/5, 3/5) and its transpose (1,5) with p(1,3) = (-4/5, 3/5).

Primitive Pair (U,W)First Chamber m0 = (U2+W2)/2First Cell (r,c)Carrier Point p(U,W)Midpoint Phase κmidRole in Classification
(1, 1)1(1, 1)(0, 1)1/2Unique diagonal carrier (45∘ record j=0)
(3, 1)5(5, 1)(4/5, 3/5)9/10Least nontrivial chamber (45∘ record j=1)
(5, 1)13(13, 1)(12/13, 5/13)25/26Second nontrivial chamber (not a 45∘ record)
(5, 3)17(13, 5)(8/17, 15/17)25/34Third nontrivial chamber; its transpose (3,5) wins for (A,C)=(3,4) at M=17
(7, 1)25(25, 1)(24/25, 7/25)49/50First composite primitive carrier chamber
(7, 3)29(25, 5)(20/29, 21/29)49/58Second nontrivial 45∘ Pell record (j=2)

Three scope clarifications govern this classification:

  1. Prime vs. Composite First Chambers: By (5.5) and Fermat's two-square theorem (p = A2 + C2 ⇒ U = A+C, W = | A-C |), an odd prime p occurs as a primitive first chamber m0 if and only if p ≡ 1 (mod 4). However, m0 is not a prime generator: (U,W) = (7,1) gives the composite first chamber m0 = 25.
  2. Non-Injectivity of m0 Across Distinct Rays: A primitive first chamber m0 need not determine a unique carrier ray, even up to transpose. At m0 = 65 = 5 · 13, the two coprime odd pairs (U,W) = (11,3) and (9,7) both satisfy (112+32)/2 = (92+72)/2 = 65, yielding first cells (61,5) and (41,25) and distinct non-mirror points (56/65, 33/65) and (16/65, 63/65).
  3. Exclusion of Opposite-Parity Targets from 𝒞(r,c): In the continuous extension κmid = a2/(a2+1), setting a = 2 gives (3/5, 4/5). Because rational centered cells require a = U/W with U, W both odd, (3/5, 4/5) is excluded from the centered rational carrier chart (whereas (4/5, 3/5) is admitted at m0 = 5 via a = 3/1). Appendix A gives the complete budget approximation theory for excluded targets and shows how adjoining the leg-sum phase κleg = r/(m+1) realizes (3/5, 4/5) at cell (4,1), chamber 4.

5.3 The Bidirectional Carrier-to-Diagonal Square Bridge

The 45∘ Pell carrier record sequence begins at (U0, W0) = (1,1) and advances by

(Uj+1, Wj+1) = (2Uj + Wj, Uj) (j ≥ 0). (5.6)

Appendix A (Theorem A.1) proves that for every chamber budget M ≥ 1, the last available ray in (5.6) with m0,j ≤ M uniquely minimizes Euclidean chord distance on 𝒞(r,c) to the 45∘ point p* = (1/√2, 1/√2). Here we prove that this carrier sequence maps bijectively to the diagonal-square Pell sequence of Proposition 4.8.

Theorem 5.2 (Bidirectional Bridge Between Carrier First Chambers and Diagonal Squares). For each primitive carrier pair (U,W), define its triangle legs and their sum by

m0 = (U2 + W2)/2, ℓ1 = (U2 - W2)/2, ℓ2 = UW, G = ℓ1 + ℓ2, δleg = ℓ1 - ℓ2. (5.7)

Then G2 + δleg2 = 2m02. In particular, | δleg | = 1 holds if and only if (U,W) = (Uj, Wj) (or its transpose for j=0) in the Pell carrier recurrence (5.6), and in that case Gj2 + 1 = 2m0,j2. The grid diagonal cell in diagonal chamber mdiag = Gj at row and column rdiag = cdiag = (Gj + 1)/2 then carries the exact square address

Ndiag(Gj) = (Gj2 + 1)/2 = m0,j2. (5.8)

Conversely, every positive integer solution (m, Q) of the diagonal-square equation m2 - 2Q2 = -1 in (4.10) determines a unique primitive carrier pair (U,W) in (5.6) with first carrier chamber m0 = Q and leg sum G = m.

Proof. Direct expansion gives

G2 + δleg2 = (ℓ1 + ℓ2)2 + (ℓ1 - ℓ2)2 = 2(ℓ12 + ℓ22) = 2(((U2 - W2)2)/4 + U2 W2) = ((U2 + W2)2)/2 = 2m02.

Along (5.6), δleg, j+1 = ((2Uj+Wj)2 - Uj2)/2 - (2Uj+Wj)Uj = -((Uj2 - Wj2)/2 - Uj Wj) = -δleg, j, so δleg, j = (-1)j+1 and Gj2 + 1 = 2m0,j2. By Proposition 4.7, the diagonal cell in chamber mdiag = Gj has address Ndiag = (Gj2 + 1)/2 = m0,j2.

Conversely, let (m, Q) be any positive integer solution of m2 + 1 = 2Q2. Set G = m and m0 = Q (with (U,W)=(1,1) for G=m0=1). For G > 1, choose δleg ∈ {-1, +1} so that ℓ1 = (G + δleg)/2 is even and ℓ2 = (G - δleg)/2 is odd. Then ℓ12 + ℓ22 = (G2 + 1)/2 = m02 and gcd(ℓ1, ℓ2) = 1. Consequently m0 + ℓ1 and m0 - ℓ1 are coprime positive odd integers whose product is m02 - ℓ12 = ℓ22, forcing m0 + ℓ1 = U2 and m0 - ℓ1 = W2 for unique coprime positive odd integers U > W with | ℓ1 - ℓ2 | = 1. By the descent in Theorem A.1, (U,W) = (Uj, Wj). □

Record Index jPrimitive Carrier Pair (Uj,Wj)Carrier First Chamber m0,jFirst Carrier Cell (r,c)Right-Triangle Legs (ℓ1,j,ℓ2,j)Grid Diagonal Chamber mdiag=GjGrid Diagonal Cell (r,r)Diagonal Square Address Ndiag=m0,j2
0(1, 1)1(1, 1)(0, 1) (degenerate)1(1, 1)1 = 12
1(3, 1)5(5, 1)(4, 3)7(4, 4)25 = 52
2(7, 3)29(25, 5)(20, 21)41(21, 21)841 = 292
3(17, 7)169(145, 25)(120, 119)239(120, 120)28,561 = 1692
4(41, 17)985(841, 145)(696, 697)1393(697, 697)970,225 = 9852

The first nondegenerate bridge is the 3–4–5 triangle (j=1): its hypotenuse is the carrier first chamber m0,1 = 5 (at carrier cell (5,1)), its leg sum is the grid diagonal chamber mdiag = 4 + 3 = 7 (at diagonal cell (4,4)), and the diagonal square address is Ndiag = 52 = 25. Naming these three distinct roles resolves the index ambiguity between the sequence 1, 5, 29, 169, 985, … and the diagonal chamber sequence 1, 7, 41, 239, 1393, ….

Figure 4: Bidirectional Carrier-to-Diagonal Square Bridge
Figure 4 (Rational carriers and the carrier-to-diagonal bridge): The centered carrier map 𝒞(r,c) of (5.1) places the odd-square cells (r̃,c̃) = λ(U2,W2) on rational points of the upper unit semicircle. The 45∘ Pell records (Uj,Wj) = (1,1),(3,1),(7,3),(17,7),… have first chambers m0,j = 1, 5, 29, 169,…, and Theorem 5.2 sends each to the diagonal chamber Gj = 1, 7, 41, 239,… whose diagonal cell carries the square address N = m0,j2.

6. A Bounded View of the Address Set

6.1 The Bounded Address Sheet, Closure, and Convex Hull

Dividing each row index by its chamber scale maps the infinite quadrant into a bounded planar region while preserving address order. Each address N, in chamber m with m + 1 = r + c ≥ 2, maps to a planar site pN in the second quadrant [7]:

m + 1 = ⌊ √(2N - 3/4) + 3/2 ⌋ ≥ 2, r = N - Tm-1, pN = (p1, p2) = (-r/(m+1), 1/(m+1)), (6.1)

where (p1, p2) ∈ ℝ2 are the coordinates of the address sheet. For each chamber m ≥ 1, the address N satisfies Tm-1 < N ≤ Tm and 1 ≤ r ≤ m. At height p2 = 1/(m+1) ∈ (0, 1/2], the row's extreme sites are (-p2, p2) (at r=1) and (p2 - 1, p2) (at r=m).

Proposition 6.1 (Discrete Sites, Topological Closure, and Isosceles Right Hull). Every site pN in (6.1) satisfies 0 < p2 ≤ 1/2 and p2 - 1 ≤ p1 ≤ -p2. The topological closure of the infinite site set {pN : N ≥ 1} in ℝ2 is precisely

{pN : N ≥ 1} = {pN : N ≥ 1} ∪ ([-1, 0] × {0}), (6.2)

so the sites are not dense in the interior of the triangle. The closed convex hull of the site set is the isosceles right triangle conv(A, B, C) with vertices

A = (-1, 0), B = (0, 0), C = (-1/2, 1/2), (6.3)

whose right angle is at the apex C = p1 (AB2 = 1, AC2 = BC2 = 1/2, (A-C)·(B-C) = 0).

Proof. Any sequence of sites converging to a point with limiting height p2 > 0 has bounded m, hence lies in a finite subset of {pN} and is eventually constant. For p2 = 0 and any a ∈ [0, 1], choosing rm ∈ {1, …, m} with rm/(m+1) → a as m → ∞ produces the limiting base point (-a, 0). Since every pN lies in conv(A,B,C) and A, B, C belong to the closure (6.2), the closed convex hull of the infinite site set is conv(A,B,C). (A finite prefix m ≤ mmax has a smaller hull; at mmax=1 it is the singleton {C}.) □

Example (Primitive columns of the address sheet). In the finite crop 1≤ r,c≤ L, write (r,c)=k(a,b), where k=gcd(r,c) and gcd(a,b)=1. Equations (1.1) and (6.1) give

N=Tk(a+b)-2+ka, p=(-a/(a+b),1/(k(a+b))).

The primitive ordered pair fixes the vertical column; the common scale fixes the height. For L=6, the cells (2,1),(4,2),(6,3) have addresses 3,14,34 and heights 1/3,1/6,1/9, all at p1=-2/3. Keeping the full point retains the pair. Discarding height identifies its common multiples. The conserved fraction is the leg-sum phase κleg = r/(r+c), distinct from the row and midpoint phases below and from the odd-leg gcd in Theorem 5.1.

Figure 5: Primitive Columns on the Address Sheet and Exact Cell Recovery
Figure 5 (Source map & coordinate recovery): On the finite crop 1≤ r,c≤ 6, writing (r,c)=k(a,b) with k=gcd(r,c) and gcd(a,b)=1 places all common multiples of a primitive cell (a,b) on the vertical column p1 = -a/(a+b) at heights p2 = 1/[k(a+b)]. For (a,b)=(2,1), the cells (2,1),(4,2),(6,3) carrying addresses N=3,14,34 lie on p1=-2/3 at heights 1/3,1/6,1/9; the open base marker (-2/3,0) is a limit point in the closure (6.2), not a discrete Grid cell.

Remark 6.1b (Euler Totient φ(m+1) Chamber Bridge, Denominator-20 Partition, and Square Coprime Count). Let Euler's totient φ(M) = #{1 ≤ j ≤ M : gcd(j,M) = 1} = M∏p | M(1 - 1/p) count the positive integers up to M that are coprime to M (φ(1)=1; φ counts coprime numerators, composite ones included, not primes). In chamber m, every cell satisfies r+c = m+1 with 1 ≤ r ≤ m and c = m+1-r, so gcd(r,c) = gcd(r, m+1). Because the omitted endpoint r = m+1 (c=0) is not coprime to m+1 for m ≥ 1, the number of primitive cells (r,c) at height 1/(m+1) of the address sheet (6.1) is exactly

#{(r,c) ∈ ℤ≥ 12 : r+c = m+1, gcd(r,c) = 1} = φ(m+1), ∑d | (m+1), d > 1 φ(d) = m (m ≥ 1).

For m=8 (r+c=9), the φ(9) = 6 primitive cells are (1,8), (2,7), (4,5), (5,4), (7,2), (8,1) (including the composite coordinates 4 and 8 because gcd(4,9)=gcd(8,9)=1). For m=19 (r+c=20), reducing the twenty fractions 1/20, …, 20/20 partitions them by reduced denominator d | 20: d=20 retains the φ(20) = 20(1-1/2)(1-1/5) = 8 primitive cells (1,19), (3,17), (7,13), (9,11), (11,9), (13,7), (17,3), (19,1) at height p2 = 1/20 (paired by r ↔ 20-r so complementary native abscissae sum to -1), while d=10,5,4,2 contribute φ(10)+φ(5)+φ(4)+φ(2) = 4+4+2+1 = 11 scaled cells lifted from earlier primitive columns by k = 20/d, and d=1 (20/20, i.e., (20,0)) lies outside the positive quadrant, giving 8 + 11 = 19 = m. Keep the predicates distinct: primitive cells in chamber m use the leg-sum denominator m+1 and count φ(m+1), whereas asking which native row fractions κrow = r/m are already in lowest terms requires gcd(r,m)=1 and counts φ(m) (and the odd-square carrier condition in Theorem 5.1 imposes the further restriction (r̃,c̃)=λ(U2,W2)). Separately, in the L × L square 1 ≤ a,b ≤ L, the lower triangle 1 ≤ b ≤ a ≤ L contains ∑k=1L φ(k) coprime pairs and shares only (1,1) with the upper triangle, so under independent uniform sampling on {1,…,L}2,

C2(L) = #{1 ≤ a,b ≤ L : gcd(a,b)=1} = 2∑k=1L φ(k) - 1, Pr(gcd(a,b)=1) = (C2(L))/L2 →L→∞ 6/π2 = 1/(ζ(2)),

by standard Möbius inversion ∑1≤ a,b≤ L, gcd(a,b)=1 1 = ∑i=1L μ(i)⌊ L/i⌋2 = 6/π2 L2 + O(L log L) [17].

Remark 6.1a (Unit-Square Frame, Half-Diagonal Length 1/√2, Row vs. Midpoint Half-Phase Pencils, and Reflected Extension C-). The address sheet (6.1) sits inside the unit square Q = [-1, 0] × [0, 1] with corners O = B = (0,0), A = (-1,0), Bul = (-1,1), and D = (0,1). Its first positive site (m = r = c = N = 1) is the unit-square center C+ = C = (-1/2, 1/2). Holding r = 1 fixes the square diagonal p1 = -p2 (OBul, of length √2), which the native hull condition p1 ≥ p2 - 1 restricts to 0 ≤ p2 ≤ 1/2; hence the native fixed-row segment OC+ has exact Euclidean length √2/2 = 1/√2 (which also equals the orthogonal projection A·e = D·e = 1/√2 of the unit edges onto the unit diagonal e = (-1,1)/√2, while Area(conv(A,O,C+)) = 1/4 = 1/2 Area(conv(A,O,D))). On this same chart, three distinct lines meet the square's geometry differently:

  1. Fixed-row ray (r = 1): p1 = -p2 through O = (0,0) and C+ = (-1/2, 1/2) (approaching O as m → ∞).
  2. Row half-phase (κrow = r/m = 1/2): p1 = (p2 - 1)/2 through the upper apex D = (0,1) and the native site P = (-1/3, 1/3) ((m,r,c) = (2,1,2), N = 2), crossing the fixed-row ray at P ≠ C+.
  3. Midpoint half-phase (κmid = (r - 1/2)/m = 1/2): the vertical line p1 = -1/2, which contains every positive diagonal site r = c (m = 2r - 1, p2r-1,r = (-1/2, 1/(2r))) including C+ = (-1/2, 1/2).

The row and midpoint half-phase lines share the boundary foot F = (-1/2, 0) (where their phase difference p2/[2(1-p2)] vanishes as p2 → 0+), while the vertical line p1 = -1/2 continues downward through F to the horizontal-axis reflection C- = (-1/2, -1/2)—an explicit reflected extension below the positive sheet with | OC+ | = | OC- | = 1/√2 and complex diagonal coincidence iω+ = ω+ = ω- for ω± = (-1 ± i)/2 (whereas the row-half-phase line p1 = (p2 - 1)/2 passes through (-3/4, -1/2) ≠ C-). This unit-square half-diagonal scale 1/√2 is distinct from the 1/√5 inner-square side arising in five-part equal-area dissections of a unit square [17].

Hull side lengths, angle bisector and the Pell target. In the unit-base normalization of (6.3), the hull 𝒦=conv(A,B,C) has side lengths 1,1/√2,1/√2. Therefore

L△=1+√2, Area(𝒦)=1/4, Rin=(2Area(𝒦))/L△ =1/(2L△)=(√2-1)/2.

The base-angle bisector gives tan(π/8)=2Rin=L△-1. Write the centered rational carrier in (5.4) as p(U,W)=(C1,C2)=(cos ψ,sin ψ), reserving (p1,p2) for the address sheet here. Then

C2/(1-C1)=U/W=cot(ψ/2).

Thus the 45∘ carrier target has U/W=cot(π/8)=L△. Under the existing Pell step (U,W)↦(2U+W,U), the signed error E=U-L△ W obeys E'=-E/L△. This connects the hull's side-length constant to the carrier's half-angle parameter and its Pell approximation. The native first site (-1/2,1/2) corresponds to cell (1,1), whose carrier (5.1) is the apex (0,1); the 45∘ point is the separate approximation target.

6.2 The Sloping Sides and Apex Circle vs. Circumcircle

Proposition 6.2 (Triangular Polynomial Difference and Sum on the Address Sheet). In planar coordinates (p1, p2) ∈ ℝ2, the difference and sum of T(p1) and T(p2 - 1) satisfy

2(T(p1) - T(p2 - 1)) = (p1 + p2)(p1 - p2 + 1) = (E12 - E22)/4, (6.4)
T(p1) + T(p2 - 1) = 0 ⇔ (p1 + 1/2)2 + (p2 - 1/2)2 = 1/2 ⇔ E12 + E22 = 2, (6.5)

where E1 = 2p1 + 1 and E2 = 2p2 - 1 center the apex C = (-1/2, 1/2) at the origin.

The locus T(p1) - T(p2 - 1) = 0 factors into the two sloping sides p1 = -p2 (BC) and p1 = p2 - 1 (AC). The circle T(p1) + T(p2 - 1) = 0 has center at the apex C = (-1/2, 1/2) and radius 1/√2; it passes through A=(-1,0) and B=(0,0), and its closed disk contains conv(A,B,C). It is not the triangle's circumcircle: the actual circumcircle passing through A, B, and C has center (-1/2, 0) and radius 1/2, with equation (p1 + 1/2)2 + p22 = 1/4.

6.3 The Scalene Affine Variant

Replacing the numerator r in (6.1) by r - (m - 1) = N - ((m+2)(m-1))/2 defines the variant sites

p̃N = ((r - m + 1)/(m + 1), 1/(m + 1)) (2 - m ≤ r - m + 1 ≤ 1). (6.6)

Proposition 6.3 (Affine Involution and Scalene Hull). The two bounded address sheets are related site-by-site by the area-preserving affine involution

(p1, p2) ↦ (-p1 + 2p2 - 1, p2) (det = -1). (6.7)

The closed convex hull of {p̃N : N ≥ 1} is the scalene triangle with vertices (-1, 0), (0, 0), (1/2, 1/2), bounding edges p2 = 0, p1 = p2, and p1 = 3p2 - 1, and squared side lengths 1, 1/2, 5/2.

6.4 Equilateral Lift and Regular k-Simplex Orthogonal Projection

The isosceles right triangle conv(A,B,C) in (6.3) is the orthogonal projection onto p3 = 0 of the spatial equilateral triangle of side 1 with vertices  = (-1, 0, 0), B̂ = (0, 0, 0), and Ĉ = (-1/2, 1/2, 1/(√2)), tilted from the horizontal plane by ψ2 = arccos(1/√3) ≈ 54.7356∘.

Proposition 6.4 (Regular k-Simplex Orthogonal Projection Factor). In ℝk+1 (k ≥ 1), let Δk = conv(e0, e1, …, ek) be the standard regular k-simplex of edge length √2, and let projΔk = conv(0, e1, …, ek) be its orthogonal projection onto the coordinate hyperplane orthogonal to e0. Then the edge vectors ei - e0 (1 ≤ i ≤ k) have Gram matrix G = Ik + 11T with det G = k + 1, and the intrinsic k-volumes satisfy

Volk(Δk) = (√(k+1))/(k!), Volk(projΔk) = 1/(k!), (Volk(projΔk))/(Volk(Δk)) = cos ψk = 1/(√(k+1)). (6.8)

Distinction Between Two Meanings of "Simplex Shadow". Proposition 6.4 projects a k-dimensional simplex orthogonally onto a k-dimensional hyperplane and tracks its intrinsic k-volume ratio 1/√(k+1) (with angles ψk = 45∘, 54.7356∘, 60∘, 63.4349∘ for k = 1, 2, 3, 4). By contrast, in Reading the Fold v1.3 (§3.2), a "simplex shadow" is the pushforward of uniform probability measure on an n-simplex to a one-dimensional line, producing the normalized Curry–Schoenberg B-spline density Mn on [-n, 0] [6]. For n = 2, the graph of M2 has vertices (-2,0), (-1,1), (0,0) and area 1; scaling both coordinates by 1/2 produces the triangle conv(A,C,B) of area 1/4, which is a geometric similarity of plane figures rather than a probability-preserving rescaling.

Remark 6.5 (Discrete Triangular Chamber Law on a Finite L × L Grid Crop and Weak Convergence to M2). While the infinite positive quadrant ℤ≥ 12 has chamber sizes #{(r,c) ∈ ℤ≥ 12 : r+c-1 = m} = m that grow linearly without bound, restricting the grid to a finite origin-anchored square crop

𝒬L = {1, 2, …, L}2 ⊂ ℤ≥ 12 (L ∈ ℤ≥ 1)

equipped with uniform discrete probability weight 1/L2 on each of its L2 cells produces an exact discrete triangular distribution for the chamber m = r + c - 1 ∈ {1, 2, …, 2L - 1} of a uniformly chosen cell (r,c):

Pr𝒬L(m = k) = (#{(r, c) ∈ {1,…,L}2 : r + c - 1 = k})/L2 = (L - | k - L |)/L2 (1 ≤ k ≤ 2L - 1). (6.8a)

Indeed, for each k ∈ {1, …, 2L-1}, the admissible rows r ∈ {1, …, L} satisfying 1 ≤ k + 1 - r ≤ L form the interval max(1, k + 1 - L) ≤ r ≤ min(L, k), which contains min(L, k) - max(1, k + 1 - L) + 1 = L - | k - L | cells, and ∑k=12L-1 (L - | k - L |) = L2.

Mapping m to the negative-support convention of [5, §3.2] via the rescaled chamber variable

AL = -(m - 1)/L = -(r - 1)/L - (c - 1)/L ∈ {0, -1/L, …, -2 + 2/L } ⊂ [-2, 0], (6.8b)

the independent discrete uniform coordinates EL = -(r-1)/L and FL = -(c-1)/L converge in distribution as L → ∞ to independent Unif([-1, 0]) variables, and AL = EL + FL converges weakly (and uniformly on the lattice mesh zk = -(k-1)/L with L Pr(AL = zk) = M2(zk) + O(1/L)) to the continuous two-fold Curry–Schoenberg tent density

M2(z) = (M1 * M1)(z) = (1 - | z + 1 |)+ (z ∈ [-2, 0]). (6.8c)

This finite-crop identity requires explicitly declaring the L × L square crop 𝒬L, uniform 1/L2 cell weighting, and the additive chamber variable m = r + c - 1; the finite distribution (6.8a) is discrete, and it provides a concrete n = 2 bridge between the Notebook Grid and [5, §3.2] without implying that empirical Gaussian or moving-average bell plots in other parts of the archive belong to the Mn spline family. Moreover, on the same uniform finite square crop 𝒬L, sorting (r, c) into the pair Y = (min(r,c), max(r,c)) (equivalently (m, | V |)) leaves the chamber m = r + c - 1 and its tent limit unchanged while mapping the L diagonal cells (r = c) to singleton fibers of mass 1/L2 and the L2 - L off-diagonal cells (r ≠ c) to equiprobable two-element fibers of mass 2/L2, forming TL-1 = L(L-1)/2 off-diagonal sorted pairs, so the exact conditional orientation information removed by sorting is H(r, c | Y) = (1 - 1/L) ln 2 (recovered by the single orientation bit I = 1{r > c}). More generally, for an exchangeable pair (A,C) taking values in a finite or countable ordered set, define YAC=(min(A,C),max(A,C)), IAC=1{A>C}, and p==Pr(A=C). Direct conditional entropy on the at-most-two-point fibers gives H(A,C | YAC)=H(IAC | YAC)=(1-p=) ln 2; this remains well-defined even when H(A,C) or H(YAC) is infinite. Moreover, because the sorted pair Y is equivalent to (m, | V |) and the address is N = (m2 + 1 + (2I - 1)| V |)/2, retaining only the chamber m discards | V | as well: with Lm = L - | m - L |, the additional conditional offset entropy lost in passing from Y to m is H(Y | m) = 1/L2 ∑m=12L-1 Lm ln Lm - (1 - 1/L) ln 2, which vanishes at L ∈ {1, 2} and is strictly positive for all L ≥ 3 (at L = 3, H(Y | m) = 1/3 ln 3 - 2/9 ln 2 > 0, witnessed by chamber m = 3 with cells (1,3), (2,2), (3,1) at addresses 4, 5, 6), while H(m) = ln L + 1/2 + o(1) matches the continuous tent differential entropy 1/2 after subtracting the ln L lattice scale.

Remark 6.4a (Tangent-Sphere Cone Ratio, Matched 2/π Inscribed Ratio, and SO+(1,1) Double-Pell Conjugacy [18]). Three three-dimensional geometric identities connect the simplex tilt cos ψ2 = 1/√3 in Proposition 6.4 to the Pell recurrences of §4:

  1. Axial Tangent-Sphere Cone Ratio: In a right circular cone of half-apex angle α ∈ (0,π/2), two spheres of radii R1 > R2 > 0 tangent to the cone wall at apex distances δi = Ri/sin α and externally tangent to each other satisfy δ1 - δ2 = R1 + R2, whence (R1 - R2)/sin α = R1 + R2 and Λcone(α) := R2/R1 = (1 - sin α)/(1 + sin α) = tan2(π/4 - α/2). At α = 45∘ (sin α = 1/√2), Λcone(45∘) = 3 - 2√2 = (√2-1)2, whose reciprocal (1+√2)2 = 3 + 2√2 generates the negative-Pell step (m,Q)↦(3m+4Q, 2m+3Q) of Proposition 4.8; at sin α = 1/√3, Λcone = 2 - √3 = (2+√3)-1.
  2. Cube 1:2:3 Squared-Radius Law and Matched Inscribed Ratio 2/π: A cube [-a,a]3 of half-side a > 0 has insphere, midsphere, and circumsphere radii (Rin, Rmid, Rout) = (a, a√2, a√3) in squared ratio 1:2:3. Under matched inscribed orientation, a square Q2 inscribed in a disk B2(R) (a = R/√2, area 2R2) and a cube Q3 inscribed in a sphere B3(R) (a = R/√3, surface area 6(2a)2 = 8R2) share the exact ratio
(Area(Q2))/(Area(B2(R))) = (2R2)/(π R2) = 2/π = (8R2)/(4π R2) = (Area(∂ Q3))/(Area(∂ B3(R))),

whereas for the dual orientation (an inscribed disk in a square) Area(B2)/Area(Q2) = π/4.

  1. Double-Pell Integrality and Positive-Sheet SO+(1,1) Boost Conjugacy: In ℝk (k > 2), mapping the body-diagonal seeds 1/√k and 1/√(k-1) by the Cayley map μ(β) = (1+β)/(1-β) gives μ(1/√D) = (D+1)/(D-1) + 2/(D-1) √D on P12 - DQ12 = 1. Both coordinates (D+1)/(D-1) and 2/(D-1) are integers iff (D-1) | 2 ⇔ D ∈ {2,3}, singling out k = 3 as the unique dimension k > 2 where both μ(1/√3) = 2+√3 (D=3) and μ(1/√2) = 3+2√2 (D=2) land on fundamental integer Pell units. Normalizing Aβ = ([1 β ; β 1]) (|β| < 1) to Λβ = (1-β2)-1/2Aβ ∈ SO+(1,1) gives Λβ Λγ = Λ(β+γ)/(1+βγ). For any Pell solution (P1,Q1) of P12 - DQ12 = 1 on the positive Lorentz sheet P1 > 0, setting CD = diag(1,√D), LD = ([P1 DQ1 ; Q1 P1]), and βD = (√D Q1)/P1 ∈ (-1,1) gives (1 - βD2)-1/2 = P1 > 0 and the exact orthochronous conjugacy
CD LD CD-1 = [ P1 √D Q1 ; √D Q1 P1 ] = ΛβD ∈ SO+(1,1)

(whereas for P1 ≠ 0 in O(1,1), (1-βD2)-1/2 = | P1| gives sgn(P1)ΛβD, witnessed at (P1,Q1) = (-1,0) where CD LD CD-1 = -I ≠ +I = Λ0). For D ∈ {2,3}, where 1/√D rationalizes to the primitive integer fundamental Pell solution (P1,Q1) = ((D+1)/(D-1),2/(D-1)) ∈ {(3,2), (2,1)}, this gives the exact seed square CD LD CD-1 = Λ1/√D2 = ΛβD. Moreover, for β ≠ 0, AβTAβ = ([1+β2 2β ; 2β 1+β2]) is not a scalar multiple of I, so Aβ cannot lie in the Euclidean similarity group ℝ×·O(2).

Remark 6.4b (Tetrahedral-to-Hexagonal Projection, Adjacent-Pair Conic Basis, D8 Fiber Centralizer, and Interstitial Volume Partition [18]):

  1. Tetrahedral Midpoints, Homothetic Hexagons, and Conic Basis: Partitioning {± 1}3 by the sign of the product of the three coordinates into dual regular tetrahedra Δ±, the edge-midpoint map M({vi,vj}) = 1/2 (vi+vj) sends the 6 edges of Δ+ bijectively onto {± e1,± e2,± e3}. Orthogonal projection P(v) = v - (v· e*)e* along e* = 1/(√3) (1,1,1) sends ±(1,1,1)↦ 0, maps the remaining 6 cube vertices to a regular hexagon of squared radius ‖ P(1,1,-1)‖2 = 3 - 1/3 = 8/3, and maps {± ei} to a homothetic regular hexagon of squared radius ‖ P(e1)‖2 = 1 - 1/3 = 2/3. Writing Σ(z) = T(z) + T(z-1) = z2, δ(z) = T(z) - T(z-1) = z, and β(a,b) = Σ(a+b)-Σ(a)-Σ(b) = 2ab gives a triangular basis (Σ(a), β(a,b), Σ(b), δ(a), δ(b), 1) for all bivariate quadratics, in which {Σ(a)+Σ(b)=1} ∩ {β(a,b)=± b} reduces to a2+b2=1 and b(2a∓ 1)=0, cutting out the six sixth roots of unity (± 1,0) and (±1/2, ±(√3)/2).
  2. Intrinsic D8 Fiber Envelope: On each 4-state fiber Qa = {a, a+6, a+12, a+18} (mod 24) over the 6 projected rays (μ24→z↦ z4μ6), cycling L(i) ≡ i+1 (mod 4) and odd reflection A(i) ≡ b-i (mod 4) (b ∈ {1,3}) satisfy A L A = L-1 and L2 ≡ (·)+2, generating ⟨ A, L⟩ ≅ D8 containing both the tetrahedral stabilizer ⟨ A, L2⟩ (a Klein four-group) and the cube rotation stabilizer C4 = ⟨ L⟩.
  3. Ball–Cube–Ball Volume Tuple: Alternating inward in ℝk between a k-ball Bn(Rn), its inscribed cube Qn (side 2Rn/√k), and the inscribed ball Bn+1(Rn/√k) partitions unit ball volume as 1 = k-k/2 + (2k/(kk/2ωk) - k-k/2) + (1 - 2k/(kk/2ωk)); at k = 3 (ω3 = (4π)/3), the step ratios form ((Vol(Qn))/(Vol(Bn)), (Vol(Bn+1))/(Vol(Qn)), (Vol(Bn+1))/(Vol(Bn))) = (2/(π√3), π/6, 1/(3√3)).

6.5 The Diagonal Staircase, Odd-Square Sectors, and Depth-Weighted Stack Projection

On the address sheet (6.1), where p(r,c) = (-r/(m+1), 1/(m+1)) satisfies p1 < 0 < p2 and the closed convex hull 𝒦 meets the upper-left diagonal boundary at p1 = -p2, the affine line p2 = (k-1)/k + p1 (k ∈ ℤ≥ 1) crosses p1 = -p2 at (p1, p2) = (-(k-1)/(2k), (k-1)/(2k)). Equivalently,

-p1 = p2 = (k-1)/(2k) = Tk-1/k2 = 1/2 - 1/(2k), Tk/k2 + Tk-1/k2 = 1 (k = 3, 4, 5, 6 ⇒ 1/3, 3/8, 2/5, 5/12). (6.9)

(Note: Since a discrete site p(r,c) = (-r/(m+1), 1/(m+1)) on p1 = -p2 requires r = 1 and m + 1 = 2k/(k-1) = 2 + 2/(k-1), for integer k ≥ 2 the crossing is a discrete site if and only if k ∈ {2, 3}, giving cell (1,3), m=3, N=4 at k=2 and cell (1,2), m=2, N=2 at k=3.) For k > 1, setting κstair = Tk/k2 = 1/2 + 1/(2k) gives 𝒥(κstair) = 1/k and 𝔣(κstair) = (Tk Tk-1)/k4 = 1/4 - 1/(4k2), climbing to the fold vertex 1/4 as k → ∞ [5, §2.3].

Square split, 1/k normalization, quadratic fixed-point family, and the far edge of the address sheet. Because Tk + Tk-1 = k2 and Tk - Tk-1 = k (k ∈ ℤ≥ 1), dividing by k2 gives the normalized gap χk+ - χk- = 1/k between χk± = 1/2 ± 1/(2k), with successive differences 1/k - 1/(k+1) = 1/(2Tk). Their common fold value 𝔣(χk±) = 1/4 - 1/(4k2) factors z2 - z + 𝔣(χk±) = (z - χk-)(z - χk+), so χk± are the real fixed points of the quadratic map z ↦ z2 + 𝔣(χk±), whose parameter approaches the parabolic value 1/4 (with multipliers 2χk± = 1 ± 1/k and deficits 1/4 - 𝔣(χk±) = 1/(4k2) summing to ζ(2)/4 = π2/24). On the address sheet (6.1), the far-edge site pm,m = (-1 + 1/(m+1), 1/(m+1)) → (-1, 0) along p1 = p2 - 1, and its raw horizontal step 1/(m+1) rescales by the active m-step span m/(m+1) to 1/m [17].

Figure 6: Diagonal Staircase Crossings, Active Span Rescaling, and Mandelbrot Cusp Fixed Points
Figure 6 (Geometric locus vs. rescaled fixed points): Left: on the address sheet (6.1), the far-edge chamber site pm,m=(-1+1/(m+1),1/(m+1)) approaches (-1,0) along p1=p2-1, while the diagonal staircase crossings (-χk-,χk-) with χk-=(k-1)/(2k) approach the apex (-1/2,1/2) along p1=-p2. Right: active-span normalization of 1/(m+1) by m/(m+1) gives 1/m, whereas square normalization gives χk+-χk-=1/k and the fixed-point pair χk±=1/2± 1/(2k) of z↦ z2+𝔣(χk±) with 𝔣(χk±)=1/4-1/(4k2).

Proposition 6.5 (Eight-Sector Lattice Congruence, Isometric Stack Fiber Multiplicity, and Atomic-to-Weak Screen Laws). Let n ≥ 1 and put j = 2n + 1.

  1. Eight-Sector Decomposition in E(2,ℤ) = ℤ2 ⋊ O(2,ℤ): Centering {1,…,j}2 at (a,b) = (A'-n-1, B'-n-1) ∈ {-n,…,n}2 partitions the j2 = 1 + 8Tn lattice points into the origin (0,0) and eight half-open max-norm sectors Σ0,…,Σ7 obtained from Σ0 = ⋃k=1n {(k,b) : 0 ≤ b < k} and Σ1 = ⋃k=1n {(a,k) : 0 < a ≤ k} under the 4-fold rotation group C4 ⊂ O(2,ℤ) ≅ D4. Each sector has cardinality Tn (indexed locally by N = Tk-1 + r ∈ {1,…,Tn}). In the lattice-preserving Euclidean semidirect product E(2,ℤ) = ℤ2 ⋊ O(2,ℤ) ⊂ E(2) = ℝ2 ⋊ O(2) with group law (b1, M1)(b2, M2) = (b1 + M1 b2, M1 M2) and action (b, M)· v = Mv + b, the lattice reflection-shift F = ((1,0), E) ∈ E(2,ℤ), with E the coordinate exchange, F(a,b) = (b+1, a), maps Σ0 →∼ Σ1 bijectively, so all eight sectors are pairwise congruent as free Euclidean subsets while remaining non-equivalent under origin-fixing maps in O(2) (at n=2, Σ0 has origin squared-norm multiset {1,4,5} whereas Σ1 has {2,5,8}).
  2. Oblique Fiber Multiplicity of the Bottom-Widest Point Stack Sn: Define the discrete point stack Sn = {(p1,p2,ℓ) ∈ ℤ3 : 0 ≤ ℓ ≤ n, | p1|, | p2| ≤ n - ℓ} of total cardinality | Sn| = ∑ℓ=0n (2(n-ℓ)+1)2 = ((n+1)(2n+1)(2n+3))/3 = C(2n+3, 3) = C(j+2, 3) (distinct from the solid cube union Ωn of Appendix B). Along e* = (1,1,1)/√3, orthogonal projection P = I - e*e*T maps (p1,p2,ℓ) to (p1-ℓ)Pe1 + (p2-ℓ)Pe2, indexed injectively by the oblique screen coordinate (a,b) = (p1-ℓ, p2-ℓ) ∈ {-n,…,n}2 with exact fiber multiplicity
wn(a,b) := #{(p1,p2,ℓ) ∈ Sn : (p1-ℓ, p2-ℓ) = (a,b)} = 1 + ⌊ (n - max(a,b))/2 ⌋, ∑a=-nn ∑b=-nn wn(a,b) = C(j+2, 3). (6.10)
  1. Finite Atomic Screen-Maximum PMFs vs. Rescaled Weak-Limit Densities: Set M = max(a+n+1, b+n+1) = max(a,b) + n + 1 ∈ {1,…,j}. At every finite odd side j = 2n+1, M obeys the exact discrete atomic probability mass functions
Prflat(M = i) = (2i - 1)/j2, Prstack(M = i) = ((2i - 1)(1 + ⌊(j - i)/2⌋))/(C(j+2, 3)) (1 ≤ i ≤ j) (6.11)

under uniform 2D screen sampling on {-n,…,n}2 and uniform 3D stack sampling on Sn, respectively. As j → ∞, the rescaled random variable Yj = M/j ∈ (0,1] satisfies the uniform local-limit bounds (ω = i/j)

max1 ≤ i ≤ j | j Prflat(M = i) - 2ω | = 1/j, max1 ≤ i ≤ j | j Prstack(M = i) - 6ω(1 - ω) | ≤ 18/j,

and hence converges weakly on [0,1] to continuous laws with densities pflat(ω) = 2ω and pstack(ω) = 6ω(1-ω) = 6 𝔣(ω).

Proof. (1) For (k,b) ∈ Σ0, 0 ≤ b < k ≤ n implies 1 ≤ b+1 ≤ k ≤ n, so F(k,b) = (b+1, k) ∈ Σ1 with inverse (a,k) ↦ (k, a-1), and F ∈ ℤ2 ⋊ O(2,ℤ) preserves Euclidean distances (‖ F(v)-F(v')‖2 = ‖ E(v-v')‖2 = ‖ v-v'‖2). Arbitrary rotations in O(2) do not preserve ℤ2, so E(2,ℤ) is written ℤ2 ⋊ O(2,ℤ) with normal translation subgroup ℤ2.

(2) Since Pe3 = -Pe1 - Pe2 and {Pe1, Pe2} is a basis of e*⊥, P(p1,p2,ℓ) = (p1-ℓ)Pe1 + (p2-ℓ)Pe2. Given (a,b) ∈ {-n,…,n}2 and an integer ℓ ≥ 0, the lifted point (a+ℓ, b+ℓ, ℓ) belongs to Sn iff ℓ ≤ n, -(n-ℓ) ≤ a+ℓ ≤ n-ℓ, and -(n-ℓ) ≤ b+ℓ ≤ n-ℓ. The lower bounds a ≥ -n and b ≥ -n already hold on {-n,…,n}2, while the upper bounds reduce to 2ℓ ≤ n - a and 2ℓ ≤ n - b, i.e., 0 ≤ ℓ ≤ ⌊(n - max(a,b))/2⌋ ≤ n, which has 1 + ⌊(n - max(a,b))/2⌋ integer solutions ℓ.

(3) Shifting (A', B') = (a+n+1, b+n+1) ∈ {1,…,j}2 gives n - max(a,b) = j - M. Exactly i2 - (i-1)2 = 2i - 1 pairs (A', B') satisfy max(A', B') = i, and wn is constant on {M=i} with value 1 + ⌊(j-i)/2⌋, establishing both PMFs in (6.11). Writing ω = i/j ∈ (0,1], j Prflat(M=i) - 2ω = -1/j. For Prstack, note (j-i+1)/2 ≤ 1 + ⌊(j-i)/2⌋ ≤ (j-i+2)/2 and C(j+2, 3) = (j(j+1)(j+2))/6. Multiplying by 6j/[j(j+1)(j+2)] gives

| j Prstack(M=i) - 6ω(1-ω) | ≤ (6| (2i-1)(1+⌊(j-i)/2⌋) - i(j-i) |)/((j+1)(j+2)) + 6ω(1-ω)(3j+2)/((j+1)(j+2)).

Since 0 ≤ 1 + ⌊(j-i)/2⌋ - (j-i)/2 ≤ 1 and 2i - 1 ≤ 2j, the first numerator has magnitude at most 6[(2j)(1) + j/2] = 15j, while 6ω(1-ω) ≤ 3/2 bounds the second term by (9j+6)/(2(j+1)(j+2)) < 3/j, yielding total error ≤ 15/j + 3/j = 18/j. Riemann-sum convergence over bounded continuous test functions on [0,1] then yields weak convergence of Yj = M/j to 2ω dω and 6ω(1-ω) dω. □

For independent uniform stages (An, Cn) ∈ {1,…,jn}2 (jn = 2n+1, center-hit event En = {An = Cn = n+1} of probability jn-2), writing 𝒮L = ∏n=1L (1 - jn-2) for the probability of no center hit in the first L stages and Π△,L(u) = ∏n=1L (1 - u/Tn) gives at every L ≥ 1 [8]:

𝒮L = ∏n=1L (8Tn)/(8Tn+1) ↓ π/4, Π△,L(-1/8) = ∏n=1L (8Tn+1)/(8Tn) ↑ 4/π, 𝒮L Π△,L(-1/8) = 1, (6.11a)

and since ∑n≥ 1Pr(En) < ∞, the number of center hits I = ∑n≥ 11En is finite almost surely, with entire probability generating function

G(τ) = 𝔼[τI] = ∏n=1∞ (1 - (1-τ)/((2n+1)2)) = π/4 Π△(-τ/8) = (cos(π/2 √(1-τ)))/τ (τ ≠ 0), G(0) = π/4, G(1) = 1. (6.11b)
Figure 7: Discrete Stack Fiber Multiplicity and Finite Atomic PMFs vs. Rescaled Weak-Limit Densities
Figure 7 (Stack projection & declared sampling measures): A. Oblique fiber depth w4(a,b) = 1 + ⌊(4 - max(a,b))/2⌋ across (a,b)∈{-4,…,4}2 summing to | S4| = C(11, 3) = 165. B. Finite odd-side discrete atomic scaled masses j Prflat(M=i) and j Prstack(M=i) (j=9, 25) converging weakly as j→∞ to 2ω and 6ω(1-ω) = 6𝔣(ω).

6.6 Common Triangular / Odd-Square Factorization, Central-Binomial Closed Form, Euler Cosine PGF, and Fixed-Gap Gamma Family

The odd-square address geometry of §§1–5, the diagonal staircase value 𝔣(χ2n+1±) = 1/4 - 1/(4(2n+1)2), and the eight-sector / stack decompositions of Proposition 6.5 all instantiate a single elementary integer identity, (2n+1)2 = 8Tn + 1 = 4n(n+1) + 1, whose finite and asymptotic consequences we prove below [8, 10].

Lemma 6.5a (Rapidity Telescoping, Central-Binomial Product, and Polynomial Rigidity of FG). For each integer n ≥ 1, let (Fn, Gn) = (n, n+1) be the balanced binary split of jn = 2n+1, with share χn = n/(2n+1) = χ2n+1-, contrast δn = 1 - 2χn = 1/(2n+1) ∈ (0,1), odds ratio On = Gn/Fn = (n+1)/n = (1+δn)/(1-δn), and rapidity ϱn := artanh(δn) = 1/2 ln((n+1)/n). Then:

  1. Multiplicative Odds & Additive Rapidity Telescoping: ∏n=1L On = L + 1 and ∑n=1L ϱn = 1/2 ln(L+1).
  2. Variance-Fold Central-Binomial Product: Writing 4𝔣(χn) = 4χn(1-χn) = 1 - δn2 = (4n(n+1))/((2n+1)2) = cosh-2(ϱn), at every finite L ≥ 1,
𝒮L = ∏n=1L 4𝔣(χn) = ∏n=1L cosh-2(ϱn) = (L+1)/((2L+1)2) 16L/(C(2L, L)2). (6.12a)
  1. Polynomial Rigidity of FG: Every symmetric polynomial Ψ(F,G) = Ψ(G,F) over a commutative ring Λ satisfying Ψ(F,0) = Ψ(0,G) = 0 factors uniquely as Ψ(F,G) = FG· Q(F+G, FG) for a polynomial Q over Λ.

Proof. (1) ∏n=1L (n+1)/n = L+1, and artanh(1/(2n+1)) = 1/2 ln((n+1)/n) telescopes to 1/2 ln(L+1). (2) Using ∏n=1L 4n(n+1) = 4L(L+1)(L!)2 and ∏n=1L (2n+1) = ((2L+1)!)/(2L L!) = (2L+1)/2L C(2L, L)L!, squaring the denominator and cancelling (L!)2 gives (4L(L+1))/((2L+1)2 4-LC(2L, L)2) = (L+1)/((2L+1)2) 16L/(C(2L, L)2). (3) Viewing Ψ ∈ (Λ[F])[G], Ψ(F,0) = 0 implies G | Ψ, so Ψ(F,G) = G· Ψ1(F,G); evaluating at F = 0 gives 0 = Ψ(0,G) = G· Ψ1(0,G) in Λ[G], whence Ψ1(0,G) = 0 and F | Ψ1, giving Ψ(F,G) = FG· Q̃(F,G) with Q̃ uniquely determined because FG is a monic regular element in each variable. Symmetry Ψ(F,G) = Ψ(G,F) implies FG(Q̃(F,G)-Q̃(G,F)) = 0, so Q̃ is symmetric and by the fundamental theorem of symmetric polynomials equals a unique Q(F+G, FG). □

Theorem 6.5 (Common Triangular / Odd-Square Factorization and Multi-Route Consequences). For each integer n ≥ 1, let jn = 2n+1 and Tn = n(n+1)/2. All three routes below rest on (2n+1)2 = 8Tn + 1 = 4n(n+1) + 1 together with their stated route hypotheses:

  1. Route 1 — Eight-Sector Partition, Balanced Split, and Finite Reciprocal Identity: The centered square {-n,…,n}2 partitions into (0,0) and eight E(2,ℤ)-congruent half-open max-norm sectors Σ0,…,Σ7 of size Tn (Proposition 6.5(1)). In any discrete binary partition F + G = 2n+1 (F,G ∈ ℤ≥ 0), writing (F,G) = (n-i, n+1+i) (or its swap) with imbalance i ∈ {0,…,n} yields the exact triangular defect sequence FG = (n-i)(n+1+i) = n(n+1) - i(i+1) = 2(Tn - Ti), which is uniquely maximized at i = 0 by the balanced split {n, n+1} (two oriented ordered states (n,n+1), (n+1,n), and C(2n+1, n) labeled cuts of K2n+1, all with FG = 2Tn = ex(2n+1; K3)). At χn = n/(2n+1),
𝔣(χn) = (n(n+1))/((2n+1)2) = (2Tn)/(8Tn+1) = 1/4 - 1/(4(2n+1)2), 1/(4𝔣(χn)) = 1 + 1/(8Tn) = (1 - 1/((2n+1)2))-1,

and at every finite L ≥ 1,

𝒮L = ∏n=1L (8Tn)/(8Tn+1) = (L+1)/((2L+1)2) 16L/(C(2L, L)2) ↓ π/4, Π△,L(-1/8) = ∏n=1L 1/(4𝔣(χn)) = ((2L+1)2)/(L+1) (C(2L, L)2)/16L = 1/𝒮L ↑ 4/π.
  1. Route 2 — Euler Cosine Product and ζ(2) PGF Slope (Under Stage Independence): Under mutually independent Bernoulli center-hit indicators 1En with Pr(En) = (2n+1)-2, I = ∑n=1∞ 1En has entire generating function G(τ) = 𝔼[τI] = ∏n=1∞ (1 - (1-τ)/((2n+1)2)), which by Euler's product cos((π z)/2) = (1-z2)∏n=1∞(1 - z2/((2n+1)2)) at z = √(1-τ) satisfies
G(τ) = (cos(π/2 √(1-τ)))/τ (τ ≠ 0), G(0) = limz→ 1(cos(π z/2))/(1-z2) = π/4 = Pr(I=0), G'(1) = ∑n=1∞ 1/((2n+1)2) = (1 - 2-2)ζ(2) - 1 = π2/8 - 1.
  1. Route 3 — Fixed-Gap Gamma Ratio and Isometric Stack Pushforward: For any real gap parameter D = k - 2 > 0, setting νn(k) = 1/2 n(n+D) (νn(3) = Tn) gives at every finite L ≥ 1
PL(D) := ∏n=1L (1 + D2/(8νn(k))) = ∏n=1L ((n + D/2)2)/(n(n+D)) = (Γ(D+1))/(Γ(1+D/2)2) (Γ(L+1+D/2)2)/(Γ(L+1) Γ(L+1+D)) →L→∞ (Γ(D+1))/(Γ(1+D/2)2),

evaluating to 4/π at D=1, (2(L+1))/(L+2) → 2 at D=2, and 32/(3π) at D=3. Separately, by Proposition 6.5(2)–(3), projecting Sn along e* = (1,1,1)/√3 weights screen coordinates by wn(a,b) = 1 + ⌊(n-max(a,b))/2⌋, carrying uniform 2D flat sampling and uniform 3D stack sampling from their exact finite-j atomic laws (6.11) to the weak-limit densities pflat(ω) = 2ω and pstack(ω) = 6ω(1-ω) = 6 𝔣(ω) on [0,1].

Proof. Route 1 follows from Proposition 6.5(1), Lemma 6.5a(2), and Wallis's asymptotic C(2L, L) ∼ 4L/(√(π L)) (so 𝒮L ∼ L/(4L2)· π L = π/4), which also follows independently from G(0) = π/4 in Route 2. In Route 2, ∑n≥ 1(| 1-τ|)/((2n+1)2) converges locally uniformly on ℂ, so G(τ) is entire; L'Hôpital's rule at z=1 gives limz→ 1(cos(π z/2))/(1-z2) = (-π/2 sin(π/2))/(-2) = π/4, and termwise logarithmic differentiation at τ=1 gives G'(1) = ∑n≥ 1(2n+1)-2 = 3/4 π2/6 - 1 = π2/8 - 1. In Route 3, iterating Γ(z+1) = zΓ(z) gives ∏n=1L (n+a) = (Γ(L+1+a))/(Γ(1+a)) for a ∈ {0, D/2, D}, establishing the exact finite Gamma identity for PL(D), and Stirling's asymptotic (Γ(L+1+D/2)2)/(Γ(L+1)Γ(L+1+D)) = 1 + O(L-1) → 1 yields the limit Γ(D+1)/Γ(1+D/2)2; the stack pushforward is proved in Proposition 6.5(2)–(3). □

Table 6.1 — Hypotheses and Proved Scope of Each Route in Theorem 6.5
Component / RouteRequired Hypotheses & InputsExact Conclusion ProvedExplicit Nonclaim / Scope Boundary
Common Algebraic CoreInteger identity (2n+1)2 = 8Tn + 1 = 4n(n+1)+1.Termwise reciprocity (8Tn)/((2n+1)2) (1+1/(8Tn))=1, finite identity 𝒮LΠ△,L(-1/8)=1, and closed form 𝒮L = (L+1)/((2L+1)2) 16L/(C(2L, L)2).Evaluating the transcendental limit π/4 requires Wallis/Stirling or Euler's cosine product.
Route 1 (Balanced Split & 8-Sector Split)Binary split F+G=2n+1 or centered square {-n,…,n}2.Balanced split {n,n+1} maximizes FG=2Tn, defect sequence 2(Tn-Ti), rigidity of FG, and 8-sector split (2n+1)2 = 1 + 8Tn in ℤ2 ⋊ O(2,ℤ).Sectors Σ0,…,Σ7 are free-congruent in ℤ2 ⋊ O(2,ℤ) ⊂ ℝ2 ⋊ O(2), not O(2)-equivalent about (0,0); unlabeled {n,n+1} ≠ C(2n+1, n) labeled cuts; forces D=1.
Route 2 (Euler Product & PGF G(τ))Euler's cosine product, removable limit at τ=0, and mutual independence of 1En.G(τ) = cos((π/2)√(1-τ))/τ (τ≠ 0), removable G(0)=π/4, and G'(1) = ∑n≥ 1(2n+1)-2 = π2/8-1.Probabilistic interpretation G(τ)=𝔼[τI] requires independent Bernoulli stages.
Route 3a (Fixed-Gap Gamma PL(D))Scalar sequence νn(k) = 1/2 n(n+D) (D=k-2>0) and Gamma recurrence.Exact finite ratio PL(D) and limit Γ(D+1)/Γ(1+D/2)2 for all D > 0.For D>1, matches scalar radial eigenvalues νn(k); does not claim an operator intertwiner or full Sk-1 spectral determinant with multiplicities.
Route 3b (Stack Pushforward)Oblique projection (p1,p2,ℓ)↦(p1-ℓ,p2-ℓ) of Sn, fiber depth wn(a,b), and two source measures at j=2n+1.Finite-j atomic PMFs Prflat, Prstack (6.11); rescaled Yj=M/j converges weakly (bounds 1/j, 18/j) to 2ω and 6ω(1-ω)=6𝔣(ω).Viewing cosine cos ψ2=1/√3 (Prop. 6.4) is an area foreshortening factor; finite-j laws are atomic, and 6𝔣(ω) requires wn(a,b).
Figure 9: Theorem 6.5 Synthesis — Common Core, Sector Congruence, Euler PGF, Gamma Family, and Stack Pushforward
Figure 9 (Theorem 6.5 across four panels): A. Common algebraic core and finite closed form 𝒮L = (L+1)/((2L+1)2) 16L/(C(2L, L)2) ↓ π/4 paired with Π△,L(-1/8) = 1/𝒮L ↑ 4/π. B. Euler cosine product and entire PGF G(τ) = cos((π/2)√(1-τ))/τ (τ≠ 0) with removable value G(0) = π/4 and tangent slope G'(1) = π2/8 - 1 at τ=1. C. Fixed-gap Gamma family PL(D) ↑ Γ(D+1)/Γ(1+D/2)2. D. Finite-j atomic masses vs. weak-limit densities pflat(ω) = 2ω and pstack(ω) = 6ω(1-ω) = 6𝔣(ω).

Remark 6.6a (Triangular Canonical Product Π△(u), Self-Contained Brouncker Continuant Defect Proof, and Completed-Zeta Scope). Because ∑n≥ 1 Tn-1 = 2, Π△(u) = ∏n=1∞ (1 - u/Tn) is entire of genus 0, order 1/2, and type π√2; writing a(a+1) = 2u gives 1 - u/Tn = ((n-a)(n+a+1))/(n(n+1)) and Euler's reflection closed form

Π△(u) = 1/(Γ(1-a) Γ(a+2)) = (sin(π a))/(π a(a+1)), Π△(-1/8) = 1/(Γ(3/2)2) = 4/π (u = -1/8 ⇔ a = -1/2).

For Brouncker's generalized continued fraction B(s) = s + Kn=1∞ ((2n-1)2)/(2s) = 4[(Γ((s+3)/4))/(Γ((s+1)/4))]2 (Re s > 0, where B(1) = 4/π and B(s)B(s+2) = (s+1)2) [16], define the continuant polynomials PLBr(s), QLBr(s) ∈ ℤ[s] on all s ∈ ℂ by the initial conditions

(P-1Br(s), Q-1Br(s)) = (1, 0), (P0Br(s), Q0Br(s)) = (s, 1),

and the three-term recurrence (L ≥ 1)

PLBr(s) = 2s PL-1Br(s) + (2L - 1)2 PL-2Br(s), QLBr(s) = 2s QL-1Br(s) + (2L - 1)2 QL-2Br(s).

For any a, b ≥ -1 and s ∈ ℂ, put ⟨ a, b⟩s := (s+1)2 QaBr(s) QbBr(s+2) - PaBr(s) PbBr(s+2) and write DLBr(s) := ⟨ L, L⟩s, CL-(s) := ⟨ L, L-1⟩s, and CL+(s) := ⟨ L-1, L⟩s (with (-1)!! := 1). At L = 0, D0Br(s) = (s+1)2 - s(s+2) = 1, C0-(s) = -s, and C0+(s) = -(s+2); at L = 1, P1Br(s) = 2s2+1 and Q1Br(s) = 2s give C1-(s) = s - 2, C1+(s) = s + 4, and D1Br(s) = 4s(s+2)(s+1)2 - (2s2+1)(2(s+2)2+1) = -9 = -(3!!)2. By induction on L ≥ 1, expanding the first or second slot via the recurrence yields

CL-(s) = 2s DL-1Br(s) + (2L-1)2 CL-1+(s) = (-1)L-1((2L-1)!!)2[2s - (s + 2L)] = (-1)L ((2L-1)!!)2 (2L - s), ; CL+(s) = 2(s+2) DL-1Br(s) + (2L-1)2 CL-1-(s) = (-1)L-1((2L-1)!!)2[2(s+2) + (2L - 2 - s)] = (-1)L-1((2L-1)!!)2 (s + 2L + 2),

and for L ≥ 2, expanding DLBr(s) = 2s CL+(s) + (2L-1)2[2(s+2)CL-1+(s) + (2L-1)2 DL-2Br(s)] gives

DLBr(s) = (-1)L ((2L-1)!!)2 [-2s(s + 2L + 2) + 2(s+2)(s + 2L) + (2L - 1)2] = (-1)L ((2L-1)!!)2 (2L + 1)2 = (-1)L ((2L+1)!!)2.

Consequently, for every L ≥ 1 and every s ∈ ℂ, DLBr(s) is independent of s and satisfies the exact finite-stage scalar normalization

Π△,L(-1/8) = ∏n=1L ((2n+1)2)/(4n(n+1)) = (((2L+1)!!)2)/(22LL!(L+1)!) = (| DLBr(s)|)/(22LL!(L+1)!) = 1/𝒮L.

(Scope: This identity is an exact finite scalar normalization of ((2N+1)!!)2, not an operator intertwiner between the continuant recurrence and triangular multiplication. Likewise, under u(s) = (s(s-1))/2 = -1/8 + ((s-1/2)2)/2, the branch point u(1/2) = -1/8 and the Riemann–Siegel phase constant -π/8 in θ(t) = t/2 log(t/(2π)) - t/2 - π/8 + 1/(48t) + O(t-3)—giving 1 - 1/8 = 7/8 in N0(t) = 1 + (θ(t))/π + S(t)—share the half-shift s - 1/2 as distinct polynomial and Gamma-phase operations. Writing Ξ△(u) := ξ(s) = 1/2 s(s-1)π-s/2Γ(s/2)ζ(s) in u = s(s-1)/2, negative integers s = -k (k ≥ 1) evaluate at u(-k) = Tk to Ξ△(Tk) = ξ(k+1), which does not uniquely determine Ξ△(u) inside its order-1/2 infinite-type class because Ξ△(u) + b Π△(u) shares all triangular values u = Tk.)*

Remark 6.6b (The Commutant of the Exchange E and the Family Ba = I + aE, Projective Cayley / Möbius Group ℳ ≅ ℝ×, and Self-Contained Euler-Transform Proof for S+ [18, 19]):

  1. The Commutant of the Exchange E and the Family Ba = I + aE: A 2× 2 real matrix L commutes with the coordinate exchange E = ([0 1 ; 1 0]) iff L ∈ 𝒜 = ℝI ⊕ ℝE = ℝP+ ⊕ ℝP- (P± = 1/2 (I ± E)), on which L ↦ L e1 is a linear isomorphism onto ℝ2 intersecting the affine family ([1 a ; a' b]) at the unique matrix Ba = I + aE = (1+a)P+ + (1-a)P- (a'=a, b=1).
  2. Projective Cayley Conjugacies & Completed Möbius Group: On ℝℙ1, the projective Cayley involution C[a:b] = [b-a : b+a] (C(a) = (1-a)/(1+a)) satisfies C(1/a) = -C(a) and C(-a) = 1/C(a), intertwining the Klein-four group {id, -a, 1/a, -1/a} ≅ C2× C2 (with ℝℙ1 fixed sets {0,∞}, {-1,1}, and ∅). Adjoining M∞(a) = 1/a to the open family Mβ(a) = (a-β)/(1-β a) (β ∈ ℝ∖{± 1}, where M2 M-1/2 = M∞) yields the completed abelian group ℳ = {Mβ : β ∈ ℝℙ1∖{-1,1}} ≅ ℝ× via the homomorphism μβ = (1+β)/(1-β) (μ∞ = -1, μβ⋆γ = μβμγ), in which M∞(a) = 1/a is the unique nonidentity involution.
  3. Telescoping Zeta Identity and Exact-Rational Euler Enclosure for S+: Define the positive telescoping series S+ := ∑n=1∞ (ln n)/(n(n+1)) [19]. Using 1/(n(n+1)) = 1/n - 1/(n+1), partial summation on 1 ≤ n ≤ M gives ∑n=1M ln n(1/n - 1/(n+1)) = ∑i=2M (ln i - ln(i-1))/i - (ln M)/(M+1). Letting M → ∞ ((ln M)/(M+1) → 0), expanding -ln(1 - 1/i) = ∑k=1∞ 1/(k ik) (all terms positive for i ≥ 2), and interchanging sums by Tonelli's theorem yields
S+ = ∑i=2∞ (-ln(1 - 1/i))/i = ∑k=1∞ (ζ(k+1) - 1)/k = ∫01 (1 - γ - ψ(2 - z))/z dz ≈ 0.7885305659115090.

For rigorous rational enclosure without Euler–Maclaurin remainders, Euler's binomial transform of the alternating zeta function, (1 - 21-s)ζ(s) = ∑n=0∞ 2-(n+1)Δn fs(1) of the completely monotone sequence fs(j) = j-s (s = k+1 ≥ 2) satisfies 0 < (-1)n Δn fs(1) ≤ 1, so its L-term truncation has positive error bounded strictly by ∑n=L∞ 2-(n+1) = 2-L, while the outer tail beyond K ≥ 1 satisfies 0 < ∑k=K+1∞ (ζ(k+1)-1)/k ≤ ∑k=K+1∞ (k+2)/(k2 2k+1) ≤ (K+3)/((K+1)2 2K+1) (at L=240, K=200, rational enclosure width < 1.57× 10-63; by contrast, S- := ∑k=2∞ (ln k)/(k(k-1)) ≈ 1.2577468869443696 gives S+ + S- ≈ 2.0462774528558786 nats). Moreover, S+ is strictly distinct from three elementary closed-form series whose values hold analytically by telescoping: (i) Stel := ∑n=2∞[ln(1 + 1/(n2-1)) - 1/n2] = 1 + ln 2 - π2/6 ≈ 0.0482131137; (ii) ∑i=1∞2/((4i-1)(4i+1)) = 1 - π/4 ≈ 0.2146018366; and (iii) ∑i=1∞1/((2i-1)(2i)(2i+1)) = ln 2 - 1/2 ≈ 0.1931471806.

Figure 12: Positive Lorentz Sheet, Balanced-Split Defects, Odd-Chamber Midpoint vs. Landau's n^2+1, and the Euler-Transform S_+ Series
Figure 12 (Geometric and analytic scope separations across §§6.4–6.6 & §8.2): A. Positive Lorentz sheet P1 > 0 for CD LD CD-1 = ΛβD ∈ SO+(1,1) vs. the negative sheet (-1,0) ↦ -I ≠ Λ0 = +I, and the D∈{2,3} fundamental seed square. B. Balanced-split defect sequence FG = 2(Tn - Ti) at n=6 (Lemma 6.5a, Theorem 6.5). C. Odd-chamber diagonal midpoint Nc = (m2+1)/2 = 2n2-2n+1 (m=2n-1) vs. the superdiagonal address 2n2 and the standalone polynomial n2+1 (odd primes p ≥ 3). D. Exact-rational Euler-transform series for S+ separated from the elementary closed-form series Stel, 1-π/4, and ln 2 - 1/2.

7. Scope Boundaries and Open Questions

7.1 Mirror-Prime Pairs in Odd Chambers (Finite Verification: 1403 ≤ m ≤ 20001)

For an odd chamber m = 2M + 1, the midpoint reflection Rmid(N) = m2 + 1 - N fixes the diagonal address Nc = (m2 + 1)/2, exchanges the ordered symmetric endpoints Nc - k ↔ Nc + k, and fixes the unordered pair {Nc - k, Nc + k} setwise across the diagonal:

Nc = (m2 + 1)/2, M = (m - 1)/2, (a, b) = (Nc - k, Nc + k) (1 ≤ k ≤ M). (7.1)

Here a + b = m2 + 1 and 2 ≤ b - a = 2k ≤ m - 1. Using distinct unordered prime pairs (a < b, i.e., k ≥ 1), direct sieve computation [9] verifies that every odd chamber in the finite interval 1403 ≤ m ≤ 20001 contains at least one mirror-prime pair (beginning at m = 1403, Nc = 984,205 with kmin = 6, pair {984,199, 984,211}). Proving existence for all odd m ≥ 1403 requires controlling prime pairs in the short moving interval [Nc - M, Nc + M] of length m - 1 ≈ √(2Nc) around Nc, which remains open.

7.2 Specific Record-Set Intersection Problem (Open Beyond {1, 5, 29})

As shown in §5.2, generic distinct primitive carrier rays can share a first chamber (such as (11,3) and (9,7) at m0 = 65). A narrower question concerns the two specific approximation record families in Appendix A: the 45∘ Pell record rays (Uj, Wj) (with | U2 - W2 - 2UW | = 2) and the a = 2 target (3/5, 4/5) record rays (with | U - 2W | = 1). The only primitive rays belonging to both record families are (1,1), (3,1), and (7,3), with first chambers m0 ∈ {1, 5, 29}. Setting τ = 2U + W ≡ 3 or 7 (mod 10) and ε = U - 2W = ± 1 on the a = 2 record chambers mk- = 10k2 + 6k + 1 and mk+ = 10k2 + 14k + 5 gives 10m0 - 1 = τ2, while Pell completeness (G2 - 2m02 = -1 for m0 = P2j+1, where P0=0, P1=1, Pn+2=2Pn+1+Pn and Qn denotes the half-Pell–Lucas sequence Q0=Q1=1, Qn+2=2Qn+1+Qn with backward extension Q-1=-1) reduces the value-set intersection P2j+1 ∈ {mk-, mk+} to the quartic Diophantine equation 50G2 = τ4 + 2τ2 - 49, or equivalently via the companion identity 10P2j+1 - 1 = Qj+22 + Qj-12 - 1 = τ2 (and via the generically degree-two unramified map (τ, G) ↦ (A,C) = (2(τ2+1), 20Gτ) = (20m0, 20Gτ)—which is injective on τ, G > 0—to positive integral points (A,C) ∈ E(ℤ) on the smooth Weierstrass elliptic curve E : C2 = (A - 2)(A2 - 200) of discriminant Δ(E) = 16 · 4 · 200 · 1962 = 491,724,800 = 213 · 52 · 74, where (A,C) ∈ E(ℤ) represents a positive intersection iff A > 0, C > 0, 20 | A, (A-2)/2 = τ2, and (A2-200)/200 = G2 for positive integers τ, G > 0 with τ ≡ 3 or 7 (mod 10), recovering m0 = A/20, ε = -1 if τ ≡ 3 (mod 10) else +1, and U = (2τ+ε)/5, W = (τ-2ε)/5). Because Qj-12 - 1 = 0 iff j ∈ {0, 1, 2}, the three intersections (m0, G, τ) ∈ {(1,1,3), (5,7,7), (29,41,17)} correspond to the particular square root τ = Qj+2 for j = 0, 1, 2 (E0=(20,60), E1=(100,980), E2=(580,13940)). Direct finite enumeration confirms that the two record-chamber sets have no further intersection up to m0 ≤ 107 (and modular quadratic-residue sieving excludes 0 ≤ j ≤ 1000 beyond j ∈ {0,1,2}), while an all-index exclusion proof—either elementary or via a complete integral-point determination on E(ℤ) with pullback filters—remains open.

Figure 8: Odd-Chamber Mirror-Prime Symmetry Across the Diagonal and Record-Set Intersection on E: C^2 = (A-2)(A^2-200)
Figure 8 (Odd-chamber mirror-prime pairs and record-set intersection on E): A. Minimal mirror-prime half-gap kmin(m)≤ M=(m-1)/2 where midpoint reflection Rmid exchanges the prime endpoints Nc-k ↔ Nc+k across Nc=(m2+1)/2 and fixes {Nc-k,Nc+k} setwise in odd chambers m=2M+1 (§7.1; verified across 1403≤ m≤ 20001, while existence for all odd m≥ 1403 remains open). B. Pullback of shared a=2 and 45∘ Pell record chambers m0∈{1,5,29} to integral points E0=(20,60), E1=(100,980), E2=(580,13940) on E: C2=(A-2)(A2-200) (A=20m0, C=20Gτ, §7.2; no further intersections occur in finite checks m0≤ 107 and 0≤ j≤ 1000, while the all-index exclusion question remains open).

7.3 Separating Equations from Exploratory Notebook Sketches

Hand-drawn annotations on early working cards (such as reciprocal ratios 9/56 ↔ 56/9 and Moessner striking diagrams) record heuristic motivation during the initial construction of the Grid [9]. Every theorem in this paper rests solely on the explicitly declared definitions, domains, and proofs given in the text.

8. Divisor Shells, Jordan Totients, and Factorial-Residual Arithmetic on the Grid Addresses

8.1 Relation to the Coordinate Layers and the Shared Identity Criterion

Because the address N(r,c) = Tm-1 + r in chamber m = r+c-1 is exactly the rank of the rational triangular entry a(m,r) = m - 1 + r/m (§1.3), the positive quadrant of §§1–7 serves directly as the coordinate system for three arithmetic structures [11, 12, 14]: (i) divisor-band and shell-numerator progressions (§8.2), (ii) k-dimensional Jordan-totient shells and their three zeta limits (§8.3), and (iii) the triangular factorial-residual calculus and its entire-function Taylor lift (§§8.4–8.4a). Table 8.0 summarizes how the coordinate layers, the odd-square factorization of §6.6, and the arithmetic theorems below fit together logically:

Table 8.0 — Directional Implication Map Across the Coordinate, Factorization, and Arithmetic Layers
Structural LayerExact Mathematical StatementDirectional Logical Relation
Coordinate Layers (§§1–5)Bijection N(r,c) = Tm-1+r, parity wedge 2N = m2+V+1 (r̃+c̃=2m), phase reflections (Rmid, Rrow, Rroot), floor branches ℓk, and odd-square semicircle carriers 𝒞(r,c) with (r̃,c̃)=λ(U2,W2).Supplies the address and phase coordinates on which both the odd-square factorization and the divisor-shell / factorial-residual operators (Rd(k,h), dm, Q(N)) act.
Central Factorization (§6.6, Thm. 6.5)Triangular / Odd-Square Factorization and Reciprocal Identity: (2n+1)2 = 8Tn + 1 = 4n(n+1)+1 ⇔ 𝔣(n/(2n+1)) = (2Tn)/(8Tn+1) ⇔ (1 - jn-2)(1 + 1/(8Tn)) = 1, tying 𝒮L = (L+1)/((2L+1)2) 16L/(C(2L, L)2) ↓ π/4 to Π△,L(-1/8) ↑ 4/π.Unites the 8-sector odd-square partition, the balanced split {n,n+1} (two oriented states (n,n+1),(n+1,n)), the staircase value 𝔣(χ2n+1±), and the Wallis product.
Analytic & Measure Pushforwards (§§6.5–6.6)(a) Entire center-hit PGF G(τ) = (cos((π/2)√(1-τ)))/τ (G(0)=π/4, G'(1)=π2/8-1); (b) Fixed-gap Gamma ratio PL(D) → (Γ(D+1))/(Γ(1+D/2)2); (c) stack fibers wn(a,b)=1+⌊(n-max(a,b))/2⌋ pushing uniform stack measure forward to pstack(ω)=6ω(1-ω)=6𝔣(ω) (vs. flat 2ω).Implied by Theorem 6.5 plus declared route hypotheses: (a) requires stage independence + Euler's cosine product; (b) replaces D=1 by scalar gap D>0; (c) requires the stack fiber multiplicity wn(a,b) along (1,1,1)/√3.
Arithmetic on the Addresses (§§8.2–8.5)(d) Shell numerators Rd(k,h)=d2 k-d+h with the square-gap value Rm(1,m-1)=m2-1 (Thm. 8.1); (e) Jordan conservation ∑d | mJk(d)=mk and three zeta limits (Thm. 8.2); (f) Factorial residual Q(N)=N/gcd(N,m!)∈{1}∪ℙ with density log 2 and entire lift G±(z) (Thm. 8.3, Prop. 8.4).Instantiates N = Tm-1+r across divisor shells and factorials, separating proved one-point / smooth density laws (6/π2, log 2) from open two-point prime occupancy (§7.1, §8.5).

Across all comparisons between geometric and arithmetic constructions, we enforce the following criterion [12]:

Principle 8.1 (Shared Identity Criterion). A repeated numerical constant, exponent, or shape across two constructions is structurally meaningful if and only if both appearances instantiate an explicit common algebraic identity, coordinate transformation, generating function, or measure pushforward. Where no shared identity exists (such as a rational shell defect -1/12 vs. ζ(-1) = -B2/2 = -1/12, a 45∘ cone denominator 84 vs. ζ(-13) = -7/84, or finite parameter counts vs. exponential-sum exponents), the coincidence is classified as a small-denominator collision and kept strictly separate. (All symbols are defined locally and summarized in Appendix C; the notation follows the corpus notation register [7].)

8.2 Divisor-Band Bijection, Shell Numerator Progressions, and Finite Shell Diagnostics

In the rational triangular grid [11], chamber m ≥ 1 consists of the m rational numbers

a(m,r) = m - 1 + r/m = (m(m-1) + r)/m, 1 ≤ r ≤ m,

whose rank is the address N(r,c) = Tm-1 + r and whose fractional part r/m is the native row phase κrow. Grouping entries by their reduced denominator d = m/gcd(m,r) gives the following exact progression structure [11, 14]:

Theorem 8.1 (Divisor-Band Bijection, Shell Numerator Progressions, and the Square-Gap Factorization). For every integer m ≥ 1 and row 1 ≤ r ≤ m, set k = gcd(m,r), d = m/k, and h = r/k so that gcd(h,d) = 1 and 1 ≤ h ≤ d. Then m = dk and r = hk, and the reduced form of a(m,r) has denominator d and numerator

Rd(k,h) = (m(m-1) + r)/(gcd(m,r)) = d2 k - d + h ≡ h - d (mod d2), gcd(h - d, d2) = 1.

Consequently:

  1. Divisor-Band Decomposition & Window Discrepancy: For each divisor d | m, exactly φ(d) entries in chamber m reduce to denominator d, realizing ∑d | mφ(d) = m. For any divisor subset D ⊆ Div(m) and β ∈ [0,1], the coprime window discrepancy Δm,D(β) := ∑d∈ D#{1≤ a≤ dβ : gcd(a,d)=1} - β∑d∈ Dφ(d) satisfies
Δm,D(β) = ∑d∈ D∑i | dμ(i)(⌊(dβ)/i⌋ - (dβ)/i), sup0≤ β≤ 1|Δm,D(β)| ≤ ∑d∈ D 2ω(d), sup0≤ β≤ 1|Δpa,{pa}(β)| = 1 - 1/p.
  1. Strict Within-Shell Uniqueness & Progression Multiplier: For each fixed shell d ≥ 1, (k,h) ↦ Rd(k,h) is injective on ℤ≥ 1 × {1 ≤ h ≤ d : gcd(h,d)=1} (so no integer or prime value ever repeats within a single shell d, although cross-shell values such as R2(2,1) = R3(1,1) = 7 can coincide). As k ≥ 1 varies, shell d traces φ(d) residue classes h - d (mod d2), each coprime to d2, and since d and d2 have the same prime factors, φ(d2) = dφ(d), giving
d2/(φ(d2)) = d2/(d φ(d)) = d/(φ(d)).
  1. Square Gap, the Factorization m2 - 1 = (m-1)(m+1), and the Midpoint vs. n2+1: At r = m - 1 (gcd(m,m-1)=1, d=m, k=1, h=m-1), m2 - a(m,m-1)2 = 2 - 1/m2 = 4(1/4 + 𝔣(χm±)), and Rm(1,m-1) = m2 - 1 = (m-1)(m+1) is composite for every m ≥ 3. On the Notebook Grid, the diagonal and superdiagonal addresses are N(n,n) = 2n2-2n+1 (diagonal, 5 at n=2) and N(n,n+1) = 2n2 (superdiagonal, 8 at n=2), not n2 and n2+1; the Landau polynomial n2+1 (where for an odd prime p there exists n ≥ 1 with p | (n2+1) iff ((-1)/p)=+1 ⇔ p ≡ 1 (mod 4)) connects to the Grid via the odd-chamber midpoint sum N + Rmid(N) = 2Nc = m2+1 (§7.1 and §8.5(3)), not via N(n,n+1).

Proof. Since gcd(h,d) = 1, gcd(Rd(k,h), d) = gcd(d2 k - d + h, d) = gcd(h, d) = 1, so Rd(k,h)/d is reduced and gcd(h-d, d2) = 1. (1) Each r ∈ {1,…,m} writes uniquely as r = (m/d)h with d | m and gcd(h,d)=1 (1≤ h≤ d), giving φ(d) entries per divisor d | m. By 1gcd(a,d)=1 = ∑i | gcd(a,d)μ(i), #{1≤ a≤ dβ : gcd(a,d)=1} = ∑i | dμ(i)⌊ dβ/i⌋; subtracting βφ(d) = ∑i | dμ(i)(dβ/i) gives the formula for Δm,D(β). Since each bracket ⌊·⌋ - (·) lies in (-1, 0] and ∑i | d|μ(i)| = 2ω(d), summing over d ∈ D bounds |Δm,D(β)| ≤ ∑d∈ D 2ω(d); for D = {pa}, Δpa,{pa}(β) = (⌊ pa β⌋ - pa β) - (⌊ pa-1β⌋ - pa-1β), which approaches -(1 - 1/p) as β ↑ 1/pa-1 and has supremum 1 - 1/p. (2) Writing Rd(k,h) - 1 = d2(k-1) + (d(d-1) + h - 1) with 0 ≤ h - 1 ≤ d - 1 places the remainder d(d-1) + h - 1 in [d(d-1), d2 - 1] ⊂ [0, d2 - 1], so Euclidean division by d2 uniquely recovers k - 1 and h - 1; φ(d2) = d2∏p | d2(1-1/p) = d· d∏p | d(1-1/p) = dφ(d). (3) Direct substitution at r=m-1 gives a(m,m-1) = m - 1/m, m2 - (m-1/m)2 = 2 - 1/m2, and Rm(1,m-1) = m2 - 1 = (m-1)(m+1) ≥ 8 composite for m ≥ 3; for an odd prime p, a2 ≡ -1 (mod p) is solvable in ℤ iff (-1)(p-1)/2 = 1 ⇔ p ≡ 1 (mod 4) by Euler's criterion. □

Remark 8.1a (Finite Shell Diagnostics, Explicit Cohort, Like-for-Like CV Weightings, and Two Variance Conventions [14]). To compare finite prime counts across shells 1 ≤ d ≤ M against the Dirichlet–PNT asymptotic heuristic without dividing by log 1 = 0, define the active shell cells and cell weights by

SM*(d) := {(k,h) : 1 ≤ k ≤ ⌊ M/d⌋, 1 ≤ h ≤ d, gcd(h,d)=1, Rd(k,h) ≥ 2}, wd := | SM*(d)|.

Because Rd(k,h) = d2 k - d + h ≥ 2 for every cell except (d,k,h) = (1,1,1) where R1(1,1) = 1, this cutoff excludes only (1,1,1), giving w1 = M - 1 and wd = φ(d)⌊ M/d⌋ for d ≥ 2. Selecting the shell cohort

𝒞M := {1 ≤ d ≤ M : wd ≥ 100}

gives |𝒞M| = 153, 479, 1000, 2000 shells at M = 250, 500, 1000, 2000. Every integer threshold from 97 to 100 gives the same four cohorts; we fix wd ≥ 100. Define the observed active-cell prime density and the cellwise PNT-in-progressions predictor

δM(d) := 1/wd ∑(k,h)∈ SM*(d) 1ℙ(Rd(k,h)), δ̂M(d) := 1/wd ∑(k,h)∈ SM*(d) d/(φ(d) log Rd(k,h)),

and compare the size-blind baseline AM(d) := φ(d)δM(d) against the size-corrected ratio QM(d) := δM(d)/δ̂M(d) on 𝒞M. Here d2/φ(d2) = d/φ(d) is exact, whereas 1/log Rd(k,h) is an asymptotic prime-density heuristic (no uniform theorem for the growing moduli d2 ≤ M2 follows from a finite table). Using average ranks for Spearman ties and population variance, write CVeq(Y) = √(1/(|𝒞M|) ∑d∈𝒞M(Y(d)-Yeq)2)/Yeq for equal-shell CV and CVw(Y) = √(∑d∈𝒞M wd(Y(d)-Yw)2 / ∑d∈𝒞M wd)/Yw for cell-weighted CV. The table reports both like-for-like weightings:

M|𝒞M|Spearman (d, AM)Spearman (d, QM)Equal-Shell CVeq(AM)Equal-Shell CVeq(QM)Cell-Weighted CVw(AM)Cell-Weighted CVw(QM)
250153+0.959+0.0020.5930.1370.6110.138
500479+0.979-0.0220.5460.0950.5840.094
10001000+0.984+0.0400.5390.0790.5790.076
20002000+0.991+0.0240.5410.0590.5830.058

Under either identical weighting scheme (CVeq: 0.539 → 0.079 or CVw: 0.579 → 0.076 at M=1000), dividing by δ̂M(d) reduces shell dispersion substantially. Similarly, at M = 1000, the two binomial-style standard-error conventions

SEdpred := √((δ̂M(d)(1 - δ̂M(d)))/wd), SEdobs := √((δM(d)(1 - δM(d)))/wd)

give 996/1000 shells satisfying |δM(d) - δ̂M(d)| ≤ 2 SEdpred (excluding d ∈ {450, 464, 680, 956}) and 994/1000 shells satisfying |δM(d) - δ̂M(d)| ≤ 2 SEdobs (additionally excluding d ∈ {328, 550}). These bands are descriptive diagnostics on deterministic prime data, not independent Bernoulli trials or confidence-coverage guarantees: higher shells reach larger numerators (Rd(k,h) ∼ d2 k ≤ dM), while AM(d) = φ(d)δM(d) omits the factor d in d/φ(d), explaining why shell stratification acts as a classical PNT-in-progressions coordinate lens rather than a new prime-distribution law.

8.3 The k-Dimensional Jordan Chamber Lift and Three Separate Zeta Limits

Lifting the row phases {r/m : 1≤ r≤ m} to the k-dimensional cube Ck(m) = {(r1/m,…,rk/m) : 1≤ ri≤ m} and assigning each tuple its reduced common denominator dm(r1,…,rk) = m/gcd(m,r1,…,rk) connects Remark 6.1b directly to Jordan's totient Jk(d) = dk∏p | d(1 - p-k) [12]:

Theorem 8.2 (Jordan Shell Conservation, Body-Shell Möbius Formula, and Three Separate Zeta Limits). For every dimension k ≥ 1 and chamber m ≥ 1:

  1. Jordan Shell Conservation, Separated Primitive Counts, and Punctured-Stack Möbius Sieve: For each divisor d | m, #{(r1,…,rk)∈{1,…,m}k : dm(r1,…,rk) = d} = Jk(d), so ∑d | m Jk(d) = mk. More generally, for any subset 𝒲 ⊆ [0,1]k with L𝒲(a) = #{r∈{1,…,a}k : r/a ∈ 𝒲}, A𝒲(d) := #{r∈{1,…,d}k : r/d∈ 𝒲, gcd(d,r1,…,rk)=1} = ∑i | dμ(i)L𝒲(d/i). Row primitivity, common-denominator primitivity, and punctured-stack visibility carry distinct predicates: (a) on the address-sheet row of chamber m, with sites (-r/(m+1), 1/(m+1)) (1 ≤ r ≤ m), the second coordinate forces gcd(m+1, r, 1) = 1 for all m sites, whereas primitive columns require gcd(r, m+1) = 1 (counting φ(m+1)); (b) on the punctured stack Sn* = Sn ∖ {(0,0,0)} of size | Sn| - 1 = C(2n+3, 3) - 1,
ASn*prim := #{(p1,p2,ℓ)∈ Sn* : gcd(| p1|,| p2|,ℓ)=1} = ∑i=1n μ(i)[C(2⌊ n/i⌋+3, 3) - 1] (AS1*prim=9, AS2*prim=25, AS3*prim=65, AS4*prim=121).
  1. Three Separate Divisor-Support Limits: Setting Fk(n) := (Jk(n))/nk = ∑i | n(μ(i))/ik = ∏p | n(1 - p-k) (and primorial Pβ = ∏p≤ β p) separates three distinct asymptotic objects:
(i) Fixedn: Fk(n) = ∑i | n(μ(i))/ik, (ii) Primorial (k≥ 2): limβ→∞ Fk(Pβ) = 1/(ζ(k)), (iii) Modulus-averaged (k≥ 1): limM→∞1/M ∑n≤ M Fk(n) = 1/(ζ(k+1)).
  1. Reconciliation of the Chamber Fraction (φ(n))/(n-1) and the One-Dimensional Density (φ(n))/n: In chamber m ≥ 1, put n = m+1; the m address-sheet sites p1 = -r/n (1 ≤ r ≤ n-1) include φ(n) with gcd(r,n) = 1, so
fnative(m) = (φ(m+1))/m = (φ(n))/(n - 1) = F1(n) + (φ(n))/(n(n-1)), 1/M ∑n=1M F1(n) = 6/π2 + O((log M)/M), 1/(M-1) ∑n=2M (φ(n))/(n-1) = 6/π2 + O((log M)/M).

Proof. (1) Putting a = gcd(m, r1, …, rk) = m/d sets (r1,…,rk) = a(r'1,…,r'k) with r' ∈ {1,…,d}k and gcd(d, r'1, …, r'k) = 1. Since 1gcd(d,r1,…,rk)=1 = ∑i | gcd(d,r1,…,rk)μ(i), writing r = iv (v ∈ {1,…,d/i}k) gives A𝒲(d) = ∑i | dμ(i)L𝒲(d/i), which for 𝒲 = [0,1]k (L𝒲(a) = ak) evaluates to ∑i | dμ(i)(d/i)k = Jk(d). Similarly, (p1,p2,ℓ) ∈ Sn* is divisible by i ≥ 1 iff (p1,p2,ℓ) = i(p'1,p'2,ℓ') with 0 ≤ iℓ' ≤ n ⇔ ℓ' ≤ ⌊ n/i⌋ and | p'1|,| p'2| ≤ ⌊(n-iℓ')/i⌋ = ⌊ n/i⌋ - ℓ', bijecting onto S⌊ n/i⌋* of size | S⌊ n/i⌋| - 1. (2) For k ≥ 2, ∑i=1∞ i-k < ∞, so Fk(Pβ) = ∏p≤ β(1 - p-k) → ∏p(1 - p-k) = 1/ζ(k). For k ≥ 1, 1/M ∑n≤ M Fk(n) = 1/M ∑i≤ M(μ(i))/ik ⌊ M/i⌋ = ∑i≤ M(μ(i))/ik+1 + O(1/M ∑i≤ M i-k) → 1/(ζ(k+1)), with error O((log M)/M) at k = 1 and O(M-1) for k ≥ 2. (3) Since 0 < (φ(n))/(n(n-1)) ≤ 1/(n-1), 1/(M-1) ∑n=2M (φ(n))/(n(n-1)) = O((log M)/M), giving the same Cesàro limit 6/π2. □

8.4 The Triangular Factorial-Residual Calculus on the Grid Addresses

On chamber Im = (Tm-1, Tm], with m(N) = ⌊√(2N)+1/2⌋ the chamber of N, the factorial-residual operator Q(N) := N/(gcd(N, m!)) acts directly on the Notebook Grid address N = N(r,c) = Tm-1 + r [12]:

Theorem 8.3 (Complete Factorial-Residual Classification, the Exceptional Prime Powers {p2}∪{27}, the Successor-Prime Rule, the Rising-Factorial gcd, and the Dickman Density log 2). For every integer N ≥ 1 in chamber m = m(N), let Q(N) = N/gcd(N, m!). Then:

  1. Prime-Residual Dichotomy & Complete Branch Classification: For every N ≥ 1, Q(N) ∈ {1} ∪ ℙ, and any prime p occurs as Q(N) = p at most once in a chamber Im. Moreover, Q(N) > 1 if and only if N belongs to one of three disjoint branches:
  • (i) Prime squares: N = p2 for a prime p (lying in chamber p ≤ m ≤ 2p-1 where vp(m!) = 1, giving Q(p2) = p);
  • (ii) The exceptional cube: N = 27 = 33 in chamber m = 7 ((r,c) = (6,2), v3(7!) = 2, Q(27) = 3; every other prime power pk with k ≥ 3 satisfies Q(pk) = 1);
  • (iii) High-prime multiples: N = ap where p > m is prime and 1 ≤ a < (m+1)/2 (giving Q(ap) = p).
  1. The Successor-Prime Rule on the First-Column Edge: On the first-column edge c = 1 (r = m, N = Tm), Q(Tm) = den(m!/Tm) = am, the reduced numerator of Tm/m! [4], satisfies
am = Q(Tm) = m + 1, if m + 1 is an odd prime, ; 1, otherwise.
  1. Pooled Dickman Density log 2 & Abscissa σc = 2: Let Fres(m) = #{N ∈ Im : Q(N) > 1} and Ares(M) = ∑m=1M Fres(m). For M ≥ 7, the m-smooth count Ssm(M) = #{N ≤ TM : P+(N) ≤ m(N)} and closed count Cres(M) = TM - Ares(M) satisfy Ssm(M) - Cres(M) = π(√TM) + 1, and by de Bruijn's smooth-number theorem [13] (Dickman's function ϱ has ϱ(2) = 1 - log 2),
limM→∞ (Ares(M))/TM = log 2 ≈ 0.693147, ℛ(s) := ∑m=1∞ (Fres(m))/ms ∼ (log 2)/(s - 2) (s ↓ 2).
  1. Half-Index Identity, Reciprocal-Triangular Sum, and the Golden Ratio at Address 1/2: For odd j = 2n+1, 8T(j-1)/2 + 1 = 8Tn + 1 = j2; -1/2 + ∑i=1M 1/Ti = 3/2 - 2/(M+1) → 3/2 = T2/(2!); and the continuous inverse T-1(a) = (√(8a+1)-1)/2 evaluates at a = 1/2 to T-1(1/2) = (√5-1)/2, so T((√5-1)/2) = 1/2.
  2. Successor Split, Rising-Factorial Interior gcd, and Denominator-Jump Trichotomy: For n ≥ 2, write Tn/n! = an/bn in lowest terms, where an ∈ {n+1, 1} (n+1 if n+1 is prime, 1 if n+1 is composite) and the cancelled part dn := (n+1)/an = gcd(n+1, 2(n-1)!) ∈ {1, n+1} satisfy an dn = n+1 and an + dn = n+2. In the rising product Πn(y) = ∏j=1n (1 + jy) = ∑k=0n ek(n) yk (e0(n) = 1, e1(n) = Tn, en(n) = n!), the interior coefficients satisfy gn := gcd(e1(n),…,en-1(n)) = an for all n ≥ 2. For n ≥ 3, the denominator jump Dn := bn / bn-1 = (n-1)dn-1/dn obeys
Dn = n(n-1) = 2Tn-1, n+1 prime, ; (n-1)/(n+1) < 1, n prime, ; (n(n-1))/(n+1) = n - 2 + 2/(n+1) ∉ ℤ, n, n+1 both composite,

so n+1 is prime iff Dn ∈ ℤ, in which case (2Tn+Dn)/2 = n2, (2Tn-Dn)/2 = n, and (2Tn)2 - Dn2 = 4n3. (On the forward ratio bn/bn+1, the single composite-to-prime index n=11 evaluates to 1/(11· 12) = 1/(2T11) = 1/132, which coincides numerically at that single index with |ζ(-9)| = B10/10 = 1/132 under Principle 8.1.)

  1. Primitive Triple Rank Bijection and the Unique Kummer Carry Excess at (j,p)=(5,3): Writing (m,r) = k(d,h) with k = gcd(m,r) and gcd(d,h) = 1, the map (k,d,h) ↦ N(kh, k(d-h)+1) = (kd(kd-1))/2 + kh is a bijection onto ℤ≥ 1, and for a uniformly chosen address N ≤ TM the scale satisfies Pr(k=a) → 6/(π2 a2) as M → ∞. For Λj := ∏i=j2j-1 i = ((2j)!)/(2 j!), δp(j) := vp(Λj) - 2vp(j!) ≤ 0 holds for every prime p ≤ j (j ≥ 2) except at the unique pair (j,p) = (5,3), where δ3(5) = +1 and Λ5/((5!)2) = 21/20.

Proof. (1) Let N ∈ Im = (Tm-1, Tm], so N ≤ (m(m+1))/2 < ((m+1)2)/2, and suppose a prime p divides Q(N), i.e., k := vp(N) > vp(m!).

  • If p > m: Then vp(m!) = 0 and N = ap with 1 ≤ a = N/p ≤ (m(m+1))/(2(m+1)) = m/2 < (m+1)/2 ≤ m (requiring m ≥ 2). Every prime factor p' of a satisfies p' ≤ a < m and vp'(a) ≤ vp'(m!) (since a | m!), while gcd(p, m!) = 1; thus gcd(ap, m!) = a and Q(N) = p (Branch iii). Because | Im| = m < p, at most one multiple of p lies in Im.
  • If p ≤ m: Write N = b pk with p ∤ b; since vp(m!) ≥ 1, we have k ≥ 2. If k = 2, then vp(m!) = 1 ⇒ p ≤ m ≤ 2p - 1, whence N ≤ T2p-1 = p(2p-1) < 2p2, forcing b = 1 and N = p2 with Q(p2) = p (Branch i). If k = 3: for p = 2, N ≥ 8 > T3 = 6 ⇒ m ≥ 4 ⇒ v2(m!) ≥ 3 = k, impossible; for p ≥ 5, m + 1 > √2 p3/2 ≥ √10 p > 3p ⇒ m ≥ 3p ⇒ vp(m!) ≥ ⌊ m/p⌋ ≥ 3, impossible; for p = 3, N = 27b > T6 = 21 ⇒ m ≥ 7, where m ≥ 9 has v3(m!) ≥ 4 > 3 and m ∈ {7,8} has N ≤ T8 = 36 < 54, forcing b = 1, N = 27 ∈ I7 (v3(7!) = 2, Q(27) = 3, Branch ii; note 32 = 9 ∉ I7). If k ≥ 4: for p = 2, N ≥ 16 > T5 ⇒ m ≥ 6 (v2(6!)=4), N ≥ 32 > T7 ⇒ m ≥ 8 (v2(8!)=7≥ 5), N ≥ 64 > T10 ⇒ m ≥ 11 (v2(11!)=8≥ 6), and for k ≥ 7, m + 1 > 2(k+1)/2 > 2k ⇒ v2(m!) ≥ ⌊ m/2⌋ + ⌊ m/4⌋ > k; for p ≥ 3, k ≥ 4, N ≥ 81 > T12 ⇒ m ≥ 13 (v3(13!)=5≥ 4), and in general m + 1 > √2 pk/2 implies ⌊ m/p⌋ ≥ ⌈√2 pk/2 - 1⌉/p - 1 ≥ k whenever (p,k) ≠ (3,4) (checked directly at (3,5): m ≥ 22, v3(22!)=9, and growing exponentially via pk/2-1 ≥ 3k/2-1). Because Branches (i)–(ii) are pure prime powers p2, 33 and Branch (iii) has a | m!, no second prime can divide Q(N), proving Q(N) ∈ {1} ∪ ℙ and within-chamber uniqueness.

(2) Since m!/Tm = (2(m-1)!)/(m+1), Q(Tm) = (m+1)/(gcd(m+1, 2(m-1)!)). At m=1 (m+1=2) and m=3 (m+1=4), m+1 divides 2(m-1)!, so Q(Tm) = 1. If m+1 is an odd prime (m ≥ 2), gcd(m+1, 2(m-1)!) = 1 so Q(Tm) = m+1. If m+1 ≥ 6 is composite, either m+1 = ab with 2 ≤ a < b ≤ (m+1)/2 ≤ m-1 (so ab | (m-1)!) or m+1 = p2 for an odd prime p ≥ 3, where 2p ≤ p2 - 3 < p2 - 1 = m for p ≥ 5 and at p = 3 (m = 8) 3 · 6 | 7!, so p2 | (m-1)! and Q(Tm) = 1.

(3) By (1), Q(N) = 1 iff P+(N) ≤ m(N) and N ∉ {p2 : p∈ℙ} ∪ {27}. Every p2 has m(p2) ≥ p = P+(p2) and 27 has m(27) = 7 > 3, so for M ≥ 7, Ssm(M) - Cres(M) = π(√TM) + 1 = O(M/log M). Because | m(N) - √(2N)| < 1/2, replacing P+(N) ≤ m(N) by P+(N) ≤ √(2TM) alters at most O(M2/log M) = o(TM) integers (since N ≤ ε TM contributes ≤ ε TM and on ε TM < N ≤ TM primes in [√(2ε TM)-1, √(2TM)+1] contribute ≪ TM/log M near the boundary). Applying de Bruijn's theorem Ψ(A, A1/2) = (1 - log 2)A + O(A/log A) [13] at A = TM gives Cres(M)/TM → 1 - log 2, so Ares(M)/TM → log 2. Partial summation on ∑m=1M Fres(m) = ((log 2)/2 + o(1))M2 yields (s-2)ℛ(s) → log 2 as s ↓ 2.

(4) ∑i=1M 1/Ti = ∑i=1M (2/i - 2/(i+1)) = 2 - 2/(M+1), and substituting a = 1/2 in T-1(a) = (√(8a+1)-1)/2 gives (√5-1)/2.

(5) Let gn = gcd(e1(n),…,en-1(n)) for n ≥ 2, and suppose a prime p divides gn. If p ≤ n, then p | n! = en(n) as well, so Πn(y) = 1 + ∑k=1n-1 ek(n)yk + n! yn ≡ 1 (mod p) as polynomials; evaluating at y = -1 gives 0 = Πn(-1) ≡ 1 (mod p) (the factor j = 1 vanishes), a contradiction. Hence no prime p ≤ n divides gn. Since gn | e1(n) = Tn = (n(n+1))/2 and every prime factor of n (and 2, since n ≥ 2) is ≤ n, gn must divide n+1 and be odd. If n+1 is composite, every prime factor of n+1 is ≤ (n+1)/2 < n, forcing gn = 1 = an. If n+1 = p is an odd prime (n = p-1), ∏j=1p-1(1+jy) ≡ 1 - yp-1 (mod p) in 𝔽p[y], so p | ek(p-1) for 1 ≤ k ≤ p-2; since vp(Tp-1) = 1, gp-1 = p = an. Next, Tn/n! = (n+1)/(2(n-1)!) = an/(2(n-1)!/dn) gives bn = (2(n-1)!)/dn, so for n ≥ 3 (where n and n+1 cannot both be prime), Dn = bn/bn-1 = (n-1)dn-1/dn yields the three stated cases.

(6) Bijectivity of (k,d,h) ↦ N(kh, k(d-h)+1) follows from Proposition 1.2 and Theorem 8.1. Since Λj/((j!)2) = (C(2j, j))/(2 j!), Kummer's theorem (vpC(2j, j) = Cp(j), the base-p carry count of j + j) gives δ2(j) = C2(j) - 1 - v2(j!) = 2 wt2(j) - 1 - j ≤ 0 for all j ≥ 2, where wt2(j) counts the ones in the binary expansion of j (δ2(2)=δ2(3)=0, and 2 wt2(j) ≤ 2⌊log2 j⌋ + 2 ≤ j + 1 for j ≥ 4), and δp(j) = Cp(j) - vp(j!) for odd primes p ≤ j. Let L = ⌊logp j⌋ + 1 be the number of base-p digits of j, so Cp(j) ≤ L. If j ≥ p2 (L ≥ 3), Legendre's sum gives vp(j!) ≥ ∑i=1L-1 pL-1-i = (pL-1-1)/(p-1) ≥ (3L-1-1)/2 > L ≥ Cp(j), so δp(j) < 0. If p ≤ j < p2 (L = 2), write j = ap + b (1 ≤ a ≤ p-1, 0 ≤ b ≤ p-1): then vp(j!) = a and Cp(j) ≤ 2, so δp(j) > 0 forces a = 1 and Cp(j) = 2. Adding (1, b)p + (1, b)p produces two base-p carries iff 2b ≥ p and 1 + 1 + 1 ≥ p, forcing p = 3, b = 2, and j = 5 (δ3(5) = 2 - 1 = +1, Λ5/((5!)2) = 21/20). □

8.4a Entire-Function Realization of the Triangular Coefficients

The factorial-weighted triangular sequence Tn/n! of Theorem 8.3(5) admits an exact entire-function realization via the Euler differential operator [15]. Define on all z ∈ ℂ the entire series and differential operators

EEi(z) := ∑n=1∞ zn/(n n!), Θz := zd/(dz), 𝒰 := 1/2 Θz2(Θz + 1), ℒ := 1/2 Θz2(Θz - 1). (8.1)

Here z ∈ ℂ is an analytic variable and n ≥ 1 a Taylor coefficient index. In the notation of [20, §6.2], EEi(z) = -Ein(-z); on any domain in ℂ∖{0} with a consistent branch of the logarithm, EEi(z) = Ei(z) - γ - Log z, and the power series in (8.1) is the branch-independent entire definition used throughout.

Proposition 8.4 (Triangular Coefficient Lift and Centered ODE). For every z ∈ ℂ,

G+(z) := 𝒰 EEi(z) = ∑n=1∞ Tn/(n!) zn = (z(z+2))/2 ez, G-(z) := ℒ EEi(z) = ∑n=1∞ Tn-1/(n!) zn = z2/2 ez. (8.2)

Consequently, G+(z) + G-(z) = z(z+1)ez and G+(z) - G-(z) = z ez, so n![zn](G+ + G-) = n2 and n![zn](G+ - G-) = n for all n ≥ 1. For every integer n ≥ 2, Theorem 8.3(5) gives

num([zn]G+(z)) = n + 1, n + 1 prime, ; 1, n + 1 composite, [zn]G-(z) = 1/(2(n - 2)!). (8.3)

Moreover, the centered entire function Ê(z) := e-z/2EEi(z) satisfies Ê(0) = 0 and the linear differential equation

Ê'(z) + 1/2 Ê(z) = (2 sinh(z/2))/z (z ∈ ℂ), (8.4)

where the right-hand side takes its removable value 1 at z = 0.

Proof. Because EEi(z) has infinite radius of convergence, termwise differentiation holds on ℂ. Since Θz(zn) = n zn, applying 𝒰 and ℒ multiplies zn/(n n!) by 1/2 n2(n± 1), producing (n(n+1))/(2 n!) zn = Tn/(n!) zn and (n(n-1))/(2 n!) zn = Tn-1/(n!) zn. Using ∑n≥ 1n/(n!) zn = z ez and ∑n≥ 1(n(n-1))/(n!) zn = z2 ez gives ∑n≥ 1Tn/(n!) zn = 1/2 (z2 ez + z ez) = (z(z+2))/2 ez and ∑n≥ 1Tn-1/(n!) zn = z2/2 ez, from which G+ ± G- and the n ≥ 2 coefficient reduction (via Theorem 8.3(5), since [zn]G+ = Tn/n!) follow immediately (noting that at n = 1, [z1]G+ = 1 with n+1 = 2, so the reduced-numerator dichotomy requires n ≥ 2). Finally, EEi'(z) = ∑n≥ 1zn-1/(n!) = (ez - 1)/z (equal to 1 at z = 0), so Ê'(z) + 1/2 Ê(z) = e-z/2EEi'(z) = (ez/2 - e-z/2)/z = (2 sinh(z/2))/z. □

Scope. Proposition 8.4 is an exact analytic realization of the triangular Taylor coefficients Tn/n! and the sums n2 and n. Passing from a rational Taylor coefficient [zn]G+ to its lowest-terms numerator num(Tn/n!) is a nonlinear arithmetic operation on the discrete index n ≥ 2; no interchange with infinite zeta-zero explicit-formula sums, prime-distribution theorem, or statement concerning the Riemann hypothesis is implied.

8.5 Residue Sieves, Odd-Prime Wheel Counts, and the Quadratic Characters of the Diagonal and the Subdiagonal

Comparing primitive lattice visibility (§8.3), shell-numerator progressions (§8.2), and symmetric prime-pair configurations (§7.1) isolates the distinction between local congruence obstruction and global prime occupancy [12]:

  1. One-Point Visibility vs. Two-Point Odd-Prime Wheels (3 ≤ p ≤ β): Removing νp forbidden residue classes modulo each prime p leaves local survival factor 1 - νp/p (or 1 - p-k on ℤk). For one-point primitive lattice visibility (k ≥ 2), inclusion-exclusion converges unconditionally to 1/ζ(k) (6/π2 at k=2). For two-point prime patterns—twin pairs {n, n+2} (νp = 2 for odd primes p ≥ 3) and Goldbach midpoint pairs {M-k, M+k} (νℓ = 2 for odd primes ℓ ∤ M, collapsing to νℓ = 1 when ℓ | M)—the Chinese Remainder Theorem gives the exact survivor counts modulo the odd-prime wheel Poddβ = ∏3≤ p≤ β p:
Aβ = ∏3 ≤ p ≤ β (p - 2), Aβ(2M) = (∏3 ≤ ℓ ≤ β, ℓ ∤ M (ℓ - 2))(∏3 ≤ ℓ ≤ β, ℓ | M (ℓ - 1)) = Aβ ∏3 ≤ ℓ ≤ β, ℓ | M (ℓ - 1)/(ℓ - 2),

recovering the exact local factor (ℓ-1)/(ℓ-2) of the Hardy–Littlewood singular series 𝔖(2M) = 2C2∏ℓ | M,ℓ≥ 3(ℓ-1)/(ℓ-2), while global prime occupancy remains obstructed by the classical sieve parity barrier.

  1. Twin-Prime Wheel-Exception Rule: Because every twin prime pair (n, n+2) with n ≥ 3 already has n odd, the parity modulus p = 2 excludes no twin prime pair; (3,5) and (5,7) fail the odd-prime wheel (3 ≤ p ≤ β) only when p ∈ {3, 5} divides one of the primes themselves, leaving n ≡ 5 (mod 6) for all twin pairs beyond (3,5).
  2. The Diagonal Midpoint Sum 2Nc = m2+1 vs. the Subdiagonal 2a2+1: In an odd chamber m, a mirror-prime pair {Nc-k, Nc+k} around Nc = (m2+1)/2 (§7.1) is a Goldbach midpoint representation of 2Nc = m2 + 1 = (Nc-k) + (Nc+k) restricted to 1 ≤ k ≤ (m-1)/2 = O(√Nc). The quadratic residue conditions on the diagonal and the subdiagonal carry distinct Legendre characters:
  • The diagonal midpoint sum 2Nc = m2 + 1 and n2 + 1: Odd prime divisors p | (n2+1) satisfy ((-1)/p) = +1 ⇔ p ≡ 1 (mod 4) (νp = 2 for p ≡ 1 (mod 4), νp = 0 for p ≡ 3 (mod 4)).
  • The subdiagonal Nsub(a) = 2a2 + 1 at (a+1, a) in even chamber m = 2a: Odd prime divisors p | (2a2+1) satisfy (2a)2 ≡ -2 (mod p), so ((-2)/p) = +1 ⇔ p ≡ 1 or 3 (mod 8) (νp = 2 for p ≡ 1,3 (mod 8), νp = 0 for p ≡ 5,7 (mod 8), with a = 1 ↦ Nsub(1) = 3 prime, a = 2 ↦ 9 = 32, a = 3 ↦ 19 prime, and a = 11 ↦ 243 = 35).
Figure 10: Divisor-Shell Numerator Diagnostics, Three Jordan-Totient Zeta Limits, the Factorial-Residual Dickman Density, and Odd-Prime Wheel Midpoint Counts
Figure 10 (Arithmetic structures on the Grid addresses across four panels): A. Active shell cohort wd ≥ 100 (§8.2, Remark 8.1a): dividing by the cellwise PNT-in-progressions predictor δ̂M(d) flattens the size-blind climb AM(d) = φ(d)δM(d) to QM(d) ≈ 1. B. The separate limits of Theorem 8.2: primorial density F2(Pβ) ↓ 1/ζ(2) = 6/π2 and modulus-averaged density 1/M ∑n≤ M F2(n) → 1/ζ(3) ≈ 0.8319. C. Cumulative factorial-residual fraction Ares(M)/TM converging to the Dickman limit 1 - ϱ(2) = log 2 ≈ 0.69315 (Theorem 8.3). D. Odd-prime wheel survivor counts A7(2M) with the midpoint factor (ℓ-1)/(ℓ-2) vs. true midpoint prime pairs D(2M) (§8.5).

8.6 Mathematical Synthesis and Remaining Open Problems

Taken together, §§1–8 establish a clean structural boundary across the Notebook Grid:

  1. Unconditionally Proved Coordinate, Carrier, Measure, and One-Point / Smooth Laws: The quadrant address bijection and its parity wedge (§1), affine quadratic paths and period-4 offset parity (§2), three chamber reflections and the anti-diagonal pair correspondence (§3), floor-selected Pell branches and diagonal square descent (§4), odd-square rational semicircle carriers and chamber budgets (§5, Appendix A), address-sheet geometry, three-dimensional directional stack shadows, and the common (2n+1)2 = 8Tn + 1 balanced-split / PGF / Gamma / stack synthesis (§6, Appendix B), divisor-shell uniqueness and the three Jordan-totient zeta limits (§§8.2–8.3), and the factorial-residual dichotomy Q(N) ∈ {1}∪ℙ with Dickman density log 2 and entire coefficient lift G±(z) (§§8.4–8.4a) all hold unconditionally with self-contained proofs.
  2. Explicitly Open All-Index and Analytic Questions: Five mathematical problems encountered along the boundary of these theorems remain open and are recorded without speculative promotion:
  • (a) All-Index Odd-Chamber Mirror Primes (§7.1): Proving that every odd chamber m ≥ 1403 contains a prime pair {Nc-k, Nc+k} (1 ≤ k ≤ (m-1)/2) around Nc = (m2+1)/2, beyond the verified finite range 1403 ≤ m ≤ 20001.
  • (b) All-Index Record-Chamber Exclusion on E(ℤ) (§7.2): Proving over all indices j ≥ 0 that m0 ∈ {1, 5, 29} (E0, E1, E2) are the only positive intersections between the 45∘ Pell record chambers and the a = 2 excluded-target record chambers on E : C2 = (A - 2)(A2 - 200).
  • (c) Two-Point Prime Occupancy vs. Sieve Parity (§8.5): Passing from exact finite odd-prime wheel survivor counts Aβ, Aβ(2M) and quadratic characters ((-1)/p), ((-2)/p) to infinite prime occupancy along n2+1, 2a2+1, or short chamber intervals.
  • (d) Optimal Multi-Prime Divisor-Window Constant (§8.2): Determining the sharp multi-prime supremum constant refining sup0≤ β≤ 1|Δm,D(β)| ≤ ∑d∈ D 2ω(d) in Theorem 8.1(1).
  • (e) Local and Meromorphic Structure of the Residual Dirichlet Series ℛ(s) (§8.4): Determining whether ℛ(s) = ∑m≥ 1 Fres(m) m-s admits meromorphic continuation left of Re s = 2 beyond its leading asymptotic (log 2)/(s - 2), and characterizing the chamber-by-chamber variance and sign changes of Fres(m) - m log 2.

Appendix A. Complete Chamber-Budget Theory for Rational Carriers

This appendix proves the complete chamber-budget approximation theorems for centered and leg-sum rational carriers on the upper unit semicircle under an integer chamber bound m ≤ M [9].

A.1 The 45∘ Centered Chamber-Budget Theorem and Stepwise Error

Fix an integer chamber budget M ≥ 1, and restrict the admissible set to positive integer cells (r,c) with m = r+c-1 ≤ M whose centered carrier point 𝒞(r,c) in (5.1) is rational. We minimize Euclidean chord distance to the 45∘ target p* = (1/√2, 1/√2). With m0, ℓ1, ℓ2, G = ℓ1 + ℓ2, and δleg = ℓ1 - ℓ2 defined as in (5.7),

G2 + δleg2 = 2m02, ‖ p(U,W) - p* ‖22 = 2 - (√2 G)/m0. (A.1)

Theorem A.1 (All-Budget 45∘ Optimum and Complete Minimizing Cell List). Let (Uj, Wj) be the Pell carrier sequence (5.6) with first chambers m0,j = (Uj2 + Wj2)/2 ∈ {1, 5, 29, 169, 985, …}. For any integer budget M ≥ 1, let j ≥ 0 be the unique index such that m0,j ≤ M < m0,j+1. Then:

  1. The unique minimizing circle point among all rational centered carriers with m ≤ M is p(Uj, Wj).
  2. The complete set of minimizing cells in the grid with m ≤ M consists of its odd-scale repeats:
(r, c) = ((λ Uj2 + 1)/2, (λ Wj2 + 1)/2), λ ∈ ℤ≥ 1 odd, λ m0,j ≤ M, (A.2)

and the number of such minimizing cells is ⌊ (⌊ M/m0,j ⌋ + 1)/2 ⌋.

  1. Let E45(M) = minm ≤ M ‖ p(U,W) - p* ‖2. Along record chambers, m0,j E45(m0,j) → 1/√2, while over all integer budgets M ≥ 1,
liminfM → ∞ M E45(M) = 1/(√2), limsupM → ∞ M E45(M) = (3 + 2√2)/(√2). (A.3)

Proof. Since δleg = (U2 - W2)/2 - UW is odd and G2 + δleg2 = 2m02, any odd prime p | gcd(δleg, m0) would divide G and hence divide ℓ1 = (G+δleg)/2 and ℓ2 = UW, contradicting gcd(U,W)=1. Moreover, for any odd prime p | δleg, coprimality gives p ∤ m0, and G2 + δleg2 = 2m02 reduces to G2 ≡ 2m02 (mod p). Thus 2 must be a quadratic residue modulo p, which excludes p = 3 and p = 5. Consequently every primitive carrier ray has either | δleg | = 1 or | δleg | ≥ 7.

Every | δleg | = 1 pair belongs to (5.6): the base is (1,1), and if U < W then δleg < -UW ≤ -3, while U = W is primitive only at (1,1). For any nonbase solution, U > W and (U - W)2 - 2W2 = ± 2, which forces 2W < U ≤ 3W (for W ≥ 3, U ≤ 2W makes the left side ≤ -W2 < -2, while at W = 1 the unique positive odd nonbase solution is U = 3; conversely U > 3W makes the left side > 2W2 ≥ 2). Hence (U', W') = (W, U - 2W) is positive, odd, and coprime, preserves | δleg | = 1, and strictly decreases U + W until (1,1) is reached. The record chambers satisfy m0,j+1 = 6m0,j - m0,j-1 for j ≥ 1, and the growth bound m0,j+1 < 6m0,j holds for all j ≥ 0 (using m0,1 = 5 < 6m0,0 = 6 at j = 0). Any non-record primitive ray available at budget M < m0,j+1 has | δleg | ≥ 7 and m0 ≤ M < m0,j+1 < 6m0,j, so

(| δleg |)/m0 ≥ 7/m0 > 7/m0,j+1 > 7/(6m0,j) > 1/m0,j = (| δleg, j |)/m0,j.

Since G/m0 = √(2 - (δleg/m0)2) for G > 0 (and G ≤ 0 cannot beat j=0), p(Uj, Wj) strictly beats every non-record ray and every earlier record i < j. Finally, √(2 - Gj/m0,j · √2) ∼ 1/(√2 m0,j) and limj→∞ m0,j+1/m0,j = 3 + 2√2 give (A.3). □

A.2 Excluded Opposite-Parity Rational Targets and the Finite-Prefix Exception

Let A, C ∈ ℤ≥ 1 be coprime integers of opposite parity, and set

D0 = A2 + C2, p* = ((A2 - C2)/D0, (2AC)/D0). (A.4)

Because A, C have opposite parity, p* is never attained on the centered rational chart 𝒞(r,c). For any primitive centered carrier pair (U,W), define the nonzero odd determinant ε = CU - AW. Direct expansion gives the exact chord-distance identity

‖ p(U,W) - p* ‖22 = (2(CU - AW)2)/(D0 m0) = (2ε2)/(D0 m0). (A.5)

Let b+ = (U+, W+) be the unique positive odd pair in 0 < U+ < 2A, 0 < W+ < 2C with CU+ - AW+ = +1, and set b- = (2A - U+, 2C - W+), which satisfies C(2A-U+) - A(2C-W+) = -1. Then all positive odd pairs with | ε | = 1 form the two canonical Bézout progressions b± + 2k(A,C) (k ≥ 0). Define the explicit tail cutoff

Mcut = max{(‖ b+ ‖22)/2, (‖ b- ‖22)/2 } ≤ 2D0 - 2(A + C) + 1. (A.6)

Theorem A.2 (Bézout Tail Theorem and Necessary Finite-Prefix Exception). For every integer budget M ≥ Mcut, the unique minimizing point for p* on the centered chart is the last available | ε | = 1 Bézout ray with m0 ≤ M, and E*(M) = minm ≤ M ‖ p(U,W) - p* ‖2 satisfies

limM → ∞ M E*(M)2 = 2/D0. (A.7)

On the finite initial prefix 1 ≤ M < Mcut, however, non-Bézout rays with | ε | ≥ 3 can win: for (A,C) = (3,4) (D0 = 25), the initial | ε | = 1 rays are b+ = (1,1) (m0 = 1) and b- = (5,7) (m0 = 37 = Mcut), and at budget M = 17 the unique winner is the non-Bézout pair (U,W) = (3,5) (m0 = 17, ε = -3), which achieves ε2/m0 = 9/17 < 1/1 and squared distance 18/425 < 2/25.

Proof. Write p(U,W) = ((U2-W2)/(U2+W2), (2UW)/(U2+W2)) and p* = ((A2-C2)/D0, (2AC)/D0) with D0 = A2+C2. Since both points lie on the unit circle,

‖ p(U,W) - p* ‖22 = 2(1 - p(U,W)· p*) = (2[(A2+C2)(U2+W2) - (A2-C2)(U2-W2) - 4ACUW])/(D0(U2+W2)) ; = (4(CU - AW)2)/(D0(U2+W2)) = (2ε2)/(D0 m0),

where ε = CU - AW and m0 = (U2+W2)/2. Because U, W are both odd and A, C have opposite parity, CU and AW have opposite parity, so ε is always an odd integer (| ε | ∈ {1, 3, 5, …}).

  1. Unique Canonical Odd Bézout Pair Construction: To construct the unique positive odd pair b+ = (U+, W+) in the rectangle 0 < U+ < 2A, 0 < W+ < 2C with CU+ - AW+ = +1, start with any integer Bézout solution (a0, b0) to Ca0 - Ab0 = +1. Since A, C have opposite parity, one coordinate of (a0, b0) is already forced odd by Ca0 - Ab0 = 1; adding (A, C) if necessary makes the other coordinate odd while preserving the determinant. All odd integer solutions then differ by integer multiples of 2(A, C). Reducing the first coordinate modulo 2A by such shifts yields a unique odd representative U+ ∈ (0, 2A). The corresponding second coordinate W+ = (CU+ - 1)/A is an odd integer; since CU+ - 1 ≥ 0, W+ cannot be negative, and because W+ is odd it cannot be zero, while U+ < 2A implies W+ < 2C. This proves existence and uniqueness of b+ = (U+, W+). The complementary pair b- = (2A - U+, 2C - W+) is likewise positive and odd in (0, 2A) × (0, 2C) with determinant C(2A-U+) - A(2C-W+) = -1. All positive odd solutions with | ε | = 1 are therefore b+ + 2k(A, C) and b- + 2k(A, C) (k ≥ 0, since k < 0 makes both coordinates negative), and each is automatically coprime because any common divisor divides | ε | = 1. Since U+ ≤ 2A - 1 and W+ ≤ 2C - 1, both initial chambers satisfy m±, 0 ≤ ((2A-1)2 + (2C-1)2)/2 = 2D0 - 2(A + C) + 1, giving (A.6).
  1. Complementary Projection Interleaving and Uniform Ratio Bound < 9: Set g = (A, C) and let ν = (g· b+)/D0 ∈ (0, 2); the complementary projection of b- is (g· b-)/D0 = 2 - ν. Because | det(g, b±) | = 1, orthogonal decomposition along and perpendicular to g gives
(‖ b+ + 2kg ‖22)/2 = D0/2 ((ν + 2k)2 + D0-2), (‖ b- + 2kg ‖22)/2 = D0/2 ((2 - ν + 2k)2 + D0-2).

The value ν = 1 is impossible: if ν = 1, then b+ - g would be a nonzero integer vector orthogonal to the primitive vector g, hence an integer multiple of (C, -A), whose determinant with g is a multiple of D0 ≥ 5 rather than ± 1. Define β = min(ν, 2 - ν) ∈ (0, 1). Sorting the two complementary progressions by chamber size (without assuming the +1 ray is shorter), their positive projections interleave strictly as

β < 2 - β < 2 + β < 4 - β < 4 + β < ⋯

with corresponding Bézout-ray chambers m̂0 < m̂1 < m̂2 < ⋯, where m̂1 = Mcut = max(m+,0, m-,0) has projection 2 - β > 1. For every index j ≥ 1, the current projection νj ≥ 2 - β > 1 and the next projection satisfies νj+1 < νj + 2, so the successive chamber ratio is strictly bounded by 9:

m̂j+1/m̂j < ((νj + 2)2 + D0-2)/(νj2 + D0-2) < 9 (j ≥ 1),

because 9(ν2 + D0-2) - ((ν+2)2 + D0-2) = 8ν2 - 4ν - 4 + 8D0-2 = 4(ν-1)(2ν+1) + 8D0-2 > 0 for all ν ≥ 1. Moreover, νj+1/νj → 1 implies m̂j+1/m̂j → 1.

  1. All-Budget Tail Optimality, Asymptotic Law, and Finite-Prefix Counterexample: For any integer budget M ≥ Mcut = m̂1, choose the last available Bézout ray j ≥ 1 so that m̂j ≤ M < m̂j+1. Any primitive odd pair outside the two Bézout progressions has odd determinant | ε | ≥ 3; if it is available at budget M, its chamber satisfies m0 ≤ M < m̂j+1, so
ε2/m0 ≥ 9/m0 ≥ 9/M > 9/m̂j+1 > 1/m̂j,

and it loses strictly to the j-th Bézout ray. Earlier Bézout rays (i < j) also lose strictly since m̂i < m̂j and ε2 = 1. Because m̂j+1/m̂j → 1 forces m̂j/M → 1, the minimum squared chord error E*(M)2 = 2/(D0 m̂j) satisfies M E*(M)2 = 2/D0 M/m̂j → 2/D0, proving (A.7). Finally, at (A,C) = (3,4) (D0 = 25), b+ = (1,1) (m̂0 = 1) and b- = (5,7) (m̂1 = Mcut = 37); enumerating the seven primitive ordered odd pairs with m0 ≤ 17—namely (1,1), (1,3), (3,1), (1,5), (5,1), (3,5), (5,3)—shows that at budget M = 17, (U,W) = (3,5) (m0 = 17, ε = -3) achieves ε2/m0 = 9/17 < 1 = 1/m̂0 and squared distance 18/425 < 2/25. □

Remark A.2a (Sharp Dichotomy for the Minimal Tail Onset Mtail(A,C)). Put D0 = A2 + C2, ν = (AU+CW)/D0, ε = CU - AW, and m̂ = m0 = (U2+W2)/2. Orthogonal decomposition D0(U2+W2) = (AU+CW)2 + ε2 gives the exact normalized chamber formula m̂ = (D0 ν2 + ε2/D0)/2, so the two initial | ε | = 1 rays r0 = (U0, W0) (ν = β = p/D0 with integer p = AU0 + CW0 ∈ [A+C, D0-1]) and r1 = (2A-U0, 2C-W0) (ν = 2 - β) have chambers m̂0 = (D0β2 + D0-1)/2 and Mcut = m̂1 = (D0(2-β)2 + D0-1)/2. On [m̂0, Mcut - 1], only r0 is available from | ε | = 1, while the minimum error among available | ε | ≥ 3 rays is non-increasing in M. A primitive positive odd pair (U,W) with | ε | ≥ 3 ties or beats r0 before Mcut iff m̂ = (D0 ν2 + ε2/D0)/2 ≥ ε2 m̂0 = (ε2 D0 β2 + ε2/D0)/2 ⇔ ν2 ≥ ε2β2 ⇔ ν ≥ | ε |β and m̂ < Mcut ⇒ ν < 2 - β. Since ν ≡ ± | ε |β (mod 2) and ν = | ε |β is the non-primitive multiple | ε |r0, any primitive intruder has ν = 2 - | ε |β and (U,W) = w| ε | = (2A - | ε | U0, 2C - | ε | W0), whose positivity and 2 - | ε |β ≥ | ε |β ⇔ β ≤ 1/| ε | imply the same conditions for | ε | = 3 (with 3β ≠ 1 since gcd(A,C)=1 ⇒ 3 ∤ D0). Conversely, when w3 = (2A - 3U0, 2C - 3W0) ∈ ℤ>02 and 3p < D0, w3 is primitive (as gcd(2A-3U0, 2C-3W0) divides both C(2A-3U0)-A(2C-3W0)=± 3 and gcd(2A,2C)=2 if 3 divided both coordinates) and its chamber m̂(w3) = (D0(2-3β)2 + 9/D0)/2 satisfies the strict bounds m̂(w3) - 9m̂0 = 2D0(1-3β) = 2(D0-3p) > 0 and Mcut - m̂(w3) = (4[p(D0-p)-1])/D0 > 0 (since p ≥ A+C ≥ 3 and D0 - p ≥ 1 are positive integers). Thus w3 enters at m̂(w3) ≤ Mcut - 1 and strictly beats r0 across [m̂(w3), Mcut - 1], proving the universal two-value dichotomy: Mtail(A,C) = Mcut if 3U0 < 2A, 3W0 < 2C, 3(AU0+CW0) < D0, and Mtail(A,C) = m̂0 otherwise.

A.3 Adjoining the Positive Leg-Sum Chart κleg and the Explicit Leg-Ratio Involution

For any cell (r,c) ∈ ℤ≥ 12, the positive leg ratio a := r/c ∈ ℚ>0 determines the leg-sum phase through the map

ψ(a) := a/(1 + a) = r/(r + c) = r/(m + 1) = κleg(r,c) ∈ (0,1).

Multiplicative inversion a ↦ 1/a = c/r (exchanging legs (r,c) ↦ (c,r) = Rmid(r,c) in chamber m = r+c-1) conjugates to phase reflection across 1/2:

ψ(1/a) = (1/a)/(1 + 1/a) = 1/(1 + a) = 1 - ψ(a), pleg(c,r) = (-(r-c)/(r+c), (2√(rc))/(r+c)),

which is the algebraic form of the complementary fraction pairing r/n ↔ (n-r)/n = c/n (n = r+c = m+1).

Using the leg-normalized phase κleg = r/(r+c) = r/(m+1) defines the alternative leg-sum carrier map

pleg(r,c) = ((r - c)/(r + c), (2√(rc))/(r + c)). (A.8)

Its rational points correspond uniquely to (r,c) = k(A2, C2) for k ∈ ℤ≥ 1 and coprime A, C ∈ ℤ≥ 1, with first chamber m = A2 + C2 - 1 and repeats at m = k(A2 + C2) - 1. In particular, every opposite-parity target (A.4) is realized exactly on the leg-sum chart at cell (A2, C2) in chamber D0 - 1 (for example, (3/5, 4/5) appears exactly at (4,1) in chamber m = 4).

Theorem A.3 (One-Chamber Advance Under the Combined Two-Chart Budget). When rational points from both 𝒞(r,c) and pleg(r,c) with chamber m ≤ M are admitted to approximate p* = (1/√2, 1/√2), let j ≥ 0 be chosen by m0,j ≤ M + 1 < m0,j+1. For j ≥ 1:

  • At M = m0,j - 1 (e.g., M = 4, 28, 168), the unique minimizing point is the coordinate swap of p(Uj, Wj), realized on the leg-sum chart at (Aj2, Cj2) where Aj = (Uj+Wj)/2 and Cj = (Uj-Wj)/2.
  • For m0,j ≤ M < m0,j+1 - 1, the minimizing points are precisely that leg-sum swapped point and the centered point p(Uj, Wj), which are distinct and equidistant from p*.

Base Budgets and All Cell Realizations. For M ∈ {1, 2, 3} (j = 0), the unique minimizing point is p0 = (0, 1), realized in either chart by the diagonal cells (r, r) with 2r - 1 ≤ M (the swapped boundary point (1, 0) is excluded because both charts require positive input cell parameters r, c ≥ 1, which forbids a zero leg). For j ≥ 1, the centered minimizing cells are the odd-scale repeats in (A.2), the leg-sum swapped point τ(pj) is realized by (r, c) = k(Aj2, Cj2) with k ∈ ℤ≥ 1 and km0,j - 1 ≤ M, and the centered point pj also has leg-sum realizations (r, c) = k(Uj2, Wj2) whenever 2km0,j - 1 ≤ M, giving ⌊(⌊ M/m0,j⌋+1)/2⌋ + ⌊(M+1)/m0,j⌋ + ⌊(M+1)/(2m0,j)⌋ chart-labeled minimizing cells in total.

Proof. First, if any candidate cell (r, c) on either chart has r < c, its carrier point (p1, p2) has p1 < 0 and p2 > 0; transposing the cell (r, c) ↦ (c, r) on the same chart and at the exact same chamber m = r+c-1 produces (-p1, p2), which has strictly larger inner product (-p1 + p2)/√2 > (p1 + p2)/√2 with the target p* = (1/√2, 1/√2). Thus every minimizing cell on either chart satisfies r ≥ c and p1 ≥ 0.

On the positive leg-sum chart (A.8), pleg(r,c) = ((r-c)/(r+c), (2√(rc))/(r+c)) is rational iff rc is a square, i.e., (r, c) = k(A2, C2) for a unique primitive pair of coprime positive integers A, C ≥ 1 and scale k ∈ ℤ≥ 1, with first leg-sum chamber m = A2 + C2 - 1, repeats at k(A2 + C2) - 1, and circle point ((A2-C2)/(A2+C2), (2AC)/(A2+C2)) (where for arbitrary coprime A, C ≥ 1, setting U = A+C and W = | A-C | gives the signed identity A2 - C2 = sgn(A-C) UW). Restricting to the non-negative half-quadrant r ≥ c ⇔ A ≥ C ≥ 1 (p1 ≥ 0), there are two parity cases:

  1. Opposite Parity (A + C odd with A > C ≥ 1): The reduced denominator is m0 = A2 + C2 ≡ 1 (mod 4), the first reduced numerator A2 - C2 > 0 is odd, and the second numerator 2AC > 0 is even. Setting U = A + C and W = A - C > 0 (coprime positive odd integers with (U2+W2)/2 = A2+C2 = m0), we have A2 - C2 = UW and 2AC = (U2 - W2) / 2. Thus this point is τ(p(U,W)) ∈ τ(ℋ) (where τ(p1, p2) = (p2, p1)) and first appears on the leg-sum chart at chamber m = m0 - 1.
  2. Both Odd (A ≥ C ≥ 1 with A ≡ C ≡ 1 (mod 2)): Dividing numerators and denominator by 2, the reduced denominator is m0 = (A2 + C2)/2 ≡ 1 (mod 4), with reduced point (((A2-C2)/2)/m0, (AC)/m0) = p(A, C) already in the centered parity class ℋ. On the centered chart it is already available at chamber m0 (via cell ((A2+1)/2, (C2+1)/2)), whereas its first realization on the leg-sum chart is (r,c) = (A2, C2) at chamber A2 + C2 - 1 = 2m0 - 1 ≥ m0, with equality 2m0 - 1 = m0 if and only if m0 = 1 at the diagonal base (A,C) = (1,1) (where both charts realize p0 = (0,1) at (1,1) in chamber 1) and strict inequality 2m0 - 1 > m0 for all m0 ≥ 5 (for example, (r,c) = (9,1) with (A,C) = (3,1) gives (4/5, 3/5) at leg-sum chamber 9 = 2(5) - 1 > 5, whereas the centered chart already realizes (4/5, 3/5) at (5,1), chamber 5).

By Cases 1 and 2, every admitted point with p1 ≥ 0 at combined budget m ≤ M is either a non-negative centered point p(U,W) ∈ ℋ (U ≥ W) with m0 ≤ M or the coordinate swap τ(p(U,W)) ∈ τ(ℋ) of a positive centered point (U > W) with m0 - 1 ≤ M ⇔ m0 ≤ M + 1. Since τ preserves the linear score p1 + p2, the maximum score p1 + p2 at combined budget M is bounded above by the centered-only optimum score Sj = Gj / m0,j at hypothetical centered budget M + 1, where j ≥ 0 is defined by m0,j ≤ M + 1 < m0,j+1 (Theorem A.1).

For j = 0 (M ∈ {1, 2, 3}), p0 = (0, 1) is admitted from chamber 1 on both charts at (r, r) (2r - 1 ≤ M), while its swap (1, 0) lies outside the positive-parameter domain r, c ≥ 1. For j ≥ 1, the integers Aj = (Uj + Wj)/2 and Cj = (Uj - Wj)/2 are positive, coprime, and of opposite parity with Aj2 + Cj2 = m0,j, so the leg-sum cell (Aj2, Cj2) has chamber m0,j - 1 ≤ M and realizes τ(p(Uj, Wj)), attaining the upper bound Sj. The unswapped point p(Uj, Wj) is admitted iff m0,j ≤ M (entering on the centered chart at m0,j and on the leg-sum chart only at 2m0,j - 1). Because a line p1 + p2 = Sj intersects the unit circle in at most the two points p(Uj, Wj) and τ(p(Uj, Wj))—which are distinct for all j ≥ 1 since Uj > Wj > 0 implies ℓ1,j ≠ ℓ2,j—no other point can tie. □

Figure A.1: All-Budget 45-Degree Pell Staircase, Excluded Target (3,4) Budget-17 Exception, and Two-Chart One-Chamber Advance
Figure A.1 (Rational-carrier chamber budgets and two-chart advance): A. Stepwise scaled chord error M· E45(M) to the 45∘ target p*=(1/√2,1/√2) on the centered chart (Theorem A.1), where finite record values m0,jE45(m0,j) approach liminf M· E45(M) = 1/√2≈ 0.707107 as j→∞ (M=1↦ 0.7654, M=5↦ 0.7089, M=29↦ 0.7072, M=169↦ 0.7071) and pre-record peaks at M=m0,j-1 approach limsup M· E45(M) = (3+2√2)/√2≈ 4.1213. B. Normalized squared error ε2/m0 = D0/2 E*(M)2 for the excluded opposite-parity target (A,C)=(3,4) (D0=25, Theorem A.2 & Remark A.2a): at budget M=17, the non-Bézout ray (U,W)=(3,5) (|ε|=3, ε2/m0=9/17<1) strictly beats b+=(1,1) before the canonical |ε|=1 Bézout tail takes over for all M≥ Mcut=37. C. Adjoining the positive leg-sum chart pleg(r,c) (Theorem A.3) realizes the coordinate swap τ(p(Uj,Wj)) one chamber earlier at M=m0,j-1∈{4,28,168,…}.

Appendix B. The Odd-Square Stack Ωn and Its Directional Shadows

This appendix establishes the all-direction orthographic union-shadow, convex-hull, hiding-cone, and isotropic-mean theorems for the three-dimensional unit-height odd-square stack Ωn, extending the tetrahedral sum ∑j=0n (2j+1)2 = C(2n+3, 3) to continuous solid projections [9].

B.1 Explicit Solid Model and Exact Union-Shadow Area Law

For each integer n ≥ 0, define the unit-height odd-square stack Ωn ⊂ ℝ3 by

Ωn = ⋃j=0n [-j - 1/2, j + 1/2]2 × [j, j+1], (B.1)

and let πv : ℝ3 → v⊥ denote orthogonal projection along a nonzero viewing vector v = (A, C, D) ≠ 0, where D is explicitly declared as the vertical third-axis component. (Orientation, Slab-Center Translation & Affine-Point vs. Linear-Camera Reflection Between Top-Widest Ωn and Bottom-Widest Ωnbot / Sn: In (B.1), slab j ∈ {0,…,n} has side 2j+1 on height [j,j+1], so Ωn is top-widest with widest slab j=n at [n,n+1]. Applying the affine point reflection R3aff(p1,p2,p3) = (p1,p2,n+1-p3) maps Ωn isometrically onto the bottom-widest solid stack Ωnbot = R3aff(Ωn) = ⋃ℓ=0n [-n+ℓ-1/2, n-ℓ+1/2]2 × [ℓ,ℓ+1] (ℓ = n-j). Each bottom slab ℓ ∈ {0,…,n} is a union of unit cubes centered at half-integer heights (p1, p2, ℓ + 1/2) with | p1|,| p2| ≤ n-ℓ; translating those slab centers downward by (0,0,1/2) via (p1,p2,ℓ+1/2) ↦ (p1,p2,ℓ) identifies them bijectively with the discrete integer point stack Sn = {(p1,p2,ℓ)∈ℤ3 : 0≤ℓ≤ n, | p1|,| p2|≤ n-ℓ} of §6.5 Proposition 6.5. On viewing directions v=(A,C,D), the induced linear directional reflection is the homogeneous differential R3lin(v) = dR3aff(v) = (A,C,-D) (distinct from the affine point map R3aff), which conjugates orthographic projection along v on Ωn to projection along R3lin(v)=(A,C,-D) on Ωnbot and preserves 𝒜n(v) and ℋn(v) identically since (B.2)–(B.8) depend only on | A|,| C|,| D|.)* Denote by 𝒜n(v) = Area2(πv(Ωn)) the actual union-shadow area and by ℋn(v) = Area2(πv(convΩn)) the projected convex-hull area.

Theorem B.1 (All-Direction Actual Union-Shadow Area Law). For every integer n ≥ 0 and v = (A, C, D) ≠ 0:

  1. Horizontal Views (D = 0): Because the open slab height bands (j, j+1) are disjoint, their projected strips have disjoint interiors and sum to
𝒜n(A, C, 0) = (n + 1)2 (| A | + | C |)/(√(A2 + C2)). (B.2)
  1. Non-Horizontal Views (D ≠ 0): Set a = min(| A |, | C |)/| D | and b = max(| A |, | C |)/| D |. Then 𝒜n(v) = (| D |)/(‖ v ‖2) 𝒜̃n(a, b), where 𝒜̃n(a, 0) = (2n+1)2 and for b > 0, with ν = b/max(1, a) and (τ)+ = max(0, τ),
𝒜̃n(a, b) = (2n+1)2 + (2n+1)(a + b) + n2[(a - 1)+ + (b - 1)+] + a/(2b) ∑j=0n-1[(ν - 1)+2 - (ν - 2j - 2)+2]. (B.3)

Writing k = min(n, ⌊ ν/2 ⌋), the first term (ν-1)+2 in (B.3) contributes across all n summands whenever ν > 1 while (ν - 2j - 2)+2 contributes only for 0 ≤ j < k, so the corner-notch correction reduces in O(1) arithmetic operations to

a/(2b) [n(ν - 1)+2 - k(ν - k - 1)2 - (k(k2 - 1))/3] = a/(2b) [n(ν - 1)+2 - kν2 + 2ν k(k + 1) - 2/3 k(k + 1)(2k + 1)],

which holds identically across all regimes ν ≤ 1, 1 < ν < 2, and ν ≥ 2.

Negative Control Against Naive Positive-Part Extrapolation. At n = 1 and v = (1/3, 3/2, 1), formula (B.3) gives quotient union area 𝒜̃1(1/3, 3/2) = 541/36 (Euclidean area 541/66), whereas the projected convex hull has quotient area 31/2 (Euclidean area 93/11) and omitting the final corner-notch sum in (B.3) gives the false value 15.

Proof.

  1. Horizontal Views (D = 0): When D = 0, orthogonal projection onto (A, C, 0)⊥ preserves the vertical coordinate p3. For each slab interior j < p3 < j + 1, the horizontal cross-section [-Rj, Rj]2 (Rj = j + 1/2) projects to a segment of length (2j + 1)(| A | + | C |)/√(A2 + C2). Because distinct slabs occupy disjoint open p3-intervals (j, j+1) whose shared horizontal boundaries have planar measure zero, the projected union area is the sum of the individual strip areas:
𝒜n(A, C, 0) = (| A | + | C |)/(√(A2 + C2)) ∑j=0n (2j + 1) = (n + 1)2 (| A | + | C |)/(√(A2 + C2)),

establishing (B.2) directly for all n ≥ 0 without taking a limit as D → 0.

  1. Quotient Projection and All-n Markov Intersection Nesting (D ≠ 0): For D ≠ 0, map each point (p1, p2, p3) ∈ ℝ3 to the horizontal plane p3 = 0 along v via (p1, p2, p3) ↦ (p1 - (A/D)p3, p2 - (C/D)p3). A planar region of quotient area 𝒜̃ in p3 = 0 lifts to an orthogonal projection of Euclidean area (| D |)/(‖ v ‖2) 𝒜̃ in v⊥ because the unit normal to p3 = 0 makes cosine | D | / ‖ v ‖2 with v / ‖ v ‖2. By sign reflection and coordinate swap on the square base [-Rj, Rj]2, we may take sorted nonnegative quotient parameters 0 ≤ a ≤ b with a = min(| A |, | C |)/| D | and b = max(| A |, | C |)/| D |. The j-th slab Ω(j) = [-Rj, Rj]2 × [j, j+1] then has quotient shadow Pj = {(p1 - ap3, p2 - bp3) : (p1, p2, p3) ∈ Ω(j)}, which is the Minkowski sum of [-Rj, Rj]2 and the segment [-(j+1)(a,b), -j(a,b)], described explicitly by
-Rj - a(j+1) ≤ p1 ≤ 1/2 + j(1 - a), -Rj - b(j+1) ≤ p2 ≤ 1/2 + j(1 - b), | bp1 - ap2 | ≤ Rj(a + b).

(For a, b > 0 these inequalities cut out a convex hexagon; if a = 0 < b they reduce to a rectangle, and if a = b = 0 they give the square [-Rj, Rj]2.) For any higher slab k > j, the lower bounds -Rk - a(k+1) < -Rj - a(j+1) and -Rk - b(k+1) < -Rj - b(j+1) strictly relax, and the diagonal strip | bp1 - ap2 | ≤ Rk(a + b) strictly widens (or is identically 0 ≤ 0 when a = b = 0). Thus the only constraints of Pk (k > j) that can cut Pj are the upper coordinate bounds

p1 ≤ p̄1(k) = 1/2 + k(1 - a), p2 ≤ p̄2(k) = 1/2 + k(1 - b).

If a ≤ 1, then p̄1(k) is nondecreasing in k, so p1 ≤ p̄1(k) is already implied on Pj by p1 ≤ p̄1(j); if a > 1, then p̄1(k) strictly decreases with k. The identical dichotomy holds for b and p̄2(k). Consequently, every active constraint from Pk (k > j) tightens monotonically with k, which proves the Markov intersection-nesting property

Pj ∩ Pn ⊆ Pj ∩ Pn-1 ⊆ ⋯ ⊆ Pj ∩ Pj+1.

Therefore, the set of shadow points whose highest covering slab is j simplifies to Δj = Pj ∖ ⋃k > j Pk = Pj ∖ Pj+1 for 0 ≤ j < n, with Δn = Pn. Since Δ0, …, Δn partition ⋃j=0n Pj without any convexity assumption on the union,

𝒜̃n(a, b) = Area2(⋃j=0n Pj) = Area2(Pn) + ∑j=0n-1 Area2(Pj ∖ Pj+1).
  1. Local Difference Area Across All (a, b) Branches and O(1) Sum Reduction: Fix j and translate Pj by (Rj + a(j+1), Rj + b(j+1)) into local nonnegative coordinates (p'1, p'2) ∈ [0, 2j + 1 + a] × [0, 2j + 1 + b]. In (p'1, p'2), Pj is the bounding rectangle [0, 2j + 1 + a] × [0, 2j + 1 + b] minus the two disjoint open corner triangles p'2 > 2j + 1 + (b/a)p'1 and p'1 > 2j + 1 + (a/b)p'2 (for a > 0; each has area ab/2), giving Area2(Pj) = (2j + 1 + a)(2j + 1 + b) - ab = (2j + 1)2 + (2j + 1)(a + b) (and Area2(Pn) = (2n+1)2 + (2n+1)(a+b)). For Pj ∖ Pj+1:
  • If b ≤ 1 (inside the hiding cone), both p̄1(k) and p̄2(k) are nondecreasing, so Pj ⊆ Pj+1 and Area2(Pj ∖ Pj+1) = 0.
  • If a = 0 < b and b > 1 (axis tilt outside the cone), the rectangle Pj = [0, 2j + 1] × [0, 2j + 1 + b] is cut at p'2 ≤ 2j + 2, leaving a protruding strip of width 2j + 1 and height b - 1, so Area2(Pj ∖ Pj+1) = (2j + 1)(b - 1).
  • If 0 < a ≤ b and b > 1, put â = min(1, a) (covering both 0 < a ≤ 1 < b, where â = a, and 1 < a ≤ b, where â = 1). In the translated (p'1, p'2) coordinates, the upper bounds p1 ≤ p̄1(j+1) and p2 ≤ p̄2(j+1) become p'1 ≤ 2j + 1 + â and p'2 ≤ 2j + 2. Thus Pj ∩ Pj+1 is the clipped rectangle [0, 2j + 1 + â] × [0, 2j + 2] minus the portions of the two corner triangles lying inside it:
  • The northwest triangle p'2 > 2j + 1 + (b/a)p'1 meets [0, 2j + 1 + â] × [0, 2j + 2] in a right triangle of legs a/b ≤ â and 1, with area a/(2b).
  • The southeast triangle p'1 > 2j + 1 + (a/b)p'2 has horizontal width (2j + 1 + â) - (2j + 1 + (a/b)p'2) = â - (a/b)p'2 inside p'1 ≤ 2j + 1 + â for 0 ≤ p'2 ≤ min(bâ/a, 2j + 2); integrating ∫0min(bâ/a, 2j+2) (â - (a/b)p'2) dp'2 gives area
(bâ2)/(2a) - a/(2b) ((bâ)/a - 2j - 2)+2.

Subtracting Area2(Pj ∩ Pj+1) = (2j + 1 + â)(2j + 2) - a/(2b) - (bâ2)/(2a) + a/(2b) ((bâ)/a - 2j - 2)+2 from Area2(Pj) = (2j + 1)2 + (2j + 1)(a + b) and using a + b - 1 - â = (a - 1)+ + (b - 1)+ together with the algebraic identity

a/(2b) + (bâ2)/(2a) - â = a/(2b) ((bâ)/a - 1)2 = a/(2b) (ν - 1)+2 (ν = b/(max(1, a)) = (bâ)/a)

gives Area2(Pj ∖ Pj+1) = (2j + 1)[(a - 1)+ + (b - 1)+] + a/(2b) [(ν - 1)+2 - (ν - 2j - 2)+2]. Summing ∑j=0n-1(2j+1) = n2 proves (B.3), and indexing i = j + 1 ∈ {1, …, k} with k = min(n, ⌊ ν/2 ⌋) evaluates ∑i=1k (ν - 2i)2 = kν2 - 2ν k(k+1) + 2/3 k(k+1)(2k+1) in O(1) operations. □

B.2 The Square Hiding Cone, Convex Hull Gap, and Global Extrema

Theorem B.2 (Square Hiding Cone, Hull Area, and Extrema). Let δ+ = (| A | - | D |)+ + (| C | - | D |)+. Then for every v = (A, C, D) ≠ 0:

  1. Projected Convex Hull Area:
ℋn(v) = (| D | (2n+1)2 + (2n+1)(| A | + | C |) + n(n+1)δ+)/(‖ v ‖2). (B.4)
  1. Exact Hiding Cone and Convexity Equivalence (n ≥ 1): The following four conditions are equivalent:
  • All lower slabs hide inside the widest top slab's shadow: πv(Ωn) = πv([-n-1/2, n+1/2]2 × [n, n+1]).
  • The actual shadow πv(Ωn) is a convex polygon.
  • The hull gap vanishes: ℋn(v) - 𝒜n(v) = 0.
  • The viewing direction lies inside the square cone
max(| A |, | C |) ≤ | D |. (B.5)

In polar angle ω from the vertical axis at azimuth ψ, the hiding boundary is ωmax(ψ) = arctan(1/max(|cos ψ|, |sin ψ|)), which equals 45∘ in the coordinate vertical planes and reaches arctan √2 ≈ 54.7356∘ along the isometric cube diagonals (± 1, ± 1, ± 1).

  1. Global Minimum, Isometric View, Global Maximum, and Sharp Hull Gaps:
minv ≠ 0 𝒜n(v) = (n + 1)2 at v = (1,0,0) or (0,1,0), 𝒜n(1,1,1) = ((2n+1)(2n+3))/(√3), (B.6)
maxv ≠ 0 𝒜n(v) = (2n+1)√((2n+1)2 + 2) at (| A| : | C| : | D|) = (1 : 1 : 2n+1), (B.7)
maxv ≠ 0(ℋn(v) - 𝒜n(v)) = √2 n at (1,1,0), maxv ≠ 0(ℋn(v) - 𝒜n(v))/(𝒜n(v)) = n/((n+1)2). (B.8)

Proof.

  1. Three-Dimensional Convex Hull and Facet Projection Formula (B.4): Write 𝒫n = {(p1, p2, p3) : 0 ≤ p3 ≤ n, | p1 |, | p2 | ≤ p3}. Each slab Ω(j) is the Minkowski sum of the base unit cube Ω0 = [-1/2, 1/2]2 × [0, 1] and the horizontal square [-j, j]2 × {j}. Since convexification commutes with Minkowski addition and the convex hull of ⋃j=0n ([-j,j]2 × {j}) is the inverted square pyramid 𝒫n, we have
convΩn = Ω0 + 𝒫n = {(p1, p2, p3) ∈ ℝ3 : 0 ≤ p3 ≤ n + 1, | p1 |, | p2 | ≤ min(p3 + 1/2, n + 1/2)}.

For n > 0, convΩn is a square frustum of height n topped by a square prism of height 1: its boundary consists of a bottom square of area 1 (normal (0,0,-1)), a top square of area (2n+1)2 (normal (0,0,1)), four vertical rectangles of area 2n + 1 (normals (± 1,0,0), (0,± 1,0)), and four slanted trapezoids of area √2 n(n+1) (normals (± 1,0,-1)/√2 and (0,± 1,-1)/√2). For any convex polyhedron, almost every line parallel to v piercing the shadow meets the boundary twice, so ℋn(v) = 1/(2‖ v ‖2) ∑F Area(F) | nF · v |. Summing the ten facets and using 1/2 (| p - | D | | + | p + | D | |) = max(| p |, | D |) = | D | + (| p | - | D |)+ together with (1 + (2n + 1)2)/2 + 2n(n+1) = (2n + 1)2 yields (B.4).

  1. Strict Gap Inequality Outside the Cone and Convexity Equivalence (n ≥ 1): For D ≠ 0, in normalized sorted coordinates 0 ≤ a ≤ b, put D+ = (a - 1)+ + (b - 1)+ and 𝒞n = ∑j=0n-1 Cj(a, b) where Cj(a, b) = a/(2b) [(ν - 1)+2 - (ν - 2j - 2)+2]. Comparing (B.4) and (B.3) gives
ℋn(v) - 𝒜n(v) = (| D |)/(‖ v ‖2) (n D+ - 𝒞n).

Inside the cone (b ≤ 1), D+ = 𝒞n = 0. Outside the cone (b > 1, so D+ > 0), each summand Cj(a,b) satisfies the strict bound 0 ≤ Cj(a,b) < D+/2:

  • If a = 0, then Cj = 0 < (b - 1)/2 = D+/2.
  • If 0 < a ≤ 1 < b, then ν = b and Cj ≤ (a(b-1)2)/(2b) < (b-1)/2 = D+/2 since a(b-1) ≤ b - 1 < b.
  • If 1 < a ≤ b, then ν = b/a and Cj ≤ (a(b/a - 1)2)/(2b) = ((b-a)2)/(2ab) < (a + b - 2)/2 = D+/2 because ab(a + b - 2) - (b - a)2 = a2(b - 1) + b2(a - 1) > 0.

Summing over 0 ≤ j < n (n ≥ 1) proves that for every non-horizontal view outside the cone,

(nδ+)/(2‖ v ‖2) < ℋn(v) - 𝒜n(v) ≤ (nδ+)/(‖ v ‖2),

while for horizontal views (D = 0, n ≥ 1), ℋn(A,C,0) - 𝒜n(A,C,0) = (n(| A | + | C |))/(√(A2+C2)) > 0. Because orthogonal projection commutes with convex hull (πv(convΩn) = convπv(Ωn)), πv(Ωn) is convex iff ℋn(v) - 𝒜n(v) = 0, which for n ≥ 1 holds iff max(| A |, | C |) ≤ | D |, where πv(Ωn) = πv(Ωn,top). (At n = 0, Ω0 is a single cube and every gap vanishes.)

  1. Global Minimum, Global Maximum Across All Four Linear Forms, and Sharp Gaps (B.6)–(B.8): Let M = (n+1)2 = 2n + 1 + n2 and write the naive lower-bound numerator as ℬ = | D | (2n + 1)2 + (2n + 1)(| A | + | C |) + n2δ+. By (B.3) and (B.2), ‖ v ‖2 𝒜n(v) ≥ ℬ in every direction. Since | p | - (| p | - | D |)+ = min(| p |, | D |),
ℬ - M(| A | + | C | + | D |) = (3n2 + 2n)| D | - n2[min(| A |, | D |) + min(| C |, | D |)] ≥ n(n + 2)| D | ≥ 0,

whence 𝒜n(v) ≥ (ℬ)/(‖ v ‖2) ≥ M (| A | + | C | + | D |)/(√(A2 + C2 + D2)) ≥ M = (n+1)2, with equality for n ≥ 1 iff D = 0 and AC = 0 (the horizontal coordinate axes).

For the global maximum, 𝒜n(v) ≤ ℋn(v) everywhere. On the nonnegative octant (A, C, D) ≥ 0 with ‖ v ‖2 = 1, ℋn(v) is the maximum of four linear forms wi · (A, C, D) according to which positive parts in (A-D)+ + (C-D)+ are active:

w0 = (2n + 1, 2n + 1, (2n + 1)2), w1 = (2n + 1 + n(n+1), 2n + 1, (2n + 1)2 - n(n+1)), w1' = (2n + 1, 2n + 1 + n(n+1), (2n + 1)2 - n(n+1)), w2 = (2n + 1 + n(n+1), 2n + 1 + n(n+1), (2n + 1)2 - 2n(n+1)).

Since ‖ w0 ‖22 = (2n + 1)4 + 2(2n + 1)2, direct expansion gives

‖ w1 ‖22 - ‖ w0 ‖22 = ‖ w1' ‖22 - ‖ w0 ‖22 = -2n2(n+1)(3n+1) < 0, ‖ w2 ‖22 - ‖ w0 ‖22 = -2n2(n+1)(5n+1) < 0 (n ≥ 1).

By Cauchy–Schwarz, every linear form is at most ‖ w0 ‖2 = (2n + 1)√((2n + 1)2 + 2), so 𝒜n(v) ≤ ℋn(v) ≤ (2n + 1)√((2n + 1)2+2) over the entire sphere S2, not only inside the cone. Equality holds precisely when (| A | : | C | : | D |) = (2n + 1 : 2n + 1 : (2n + 1)2) = (1 : 1 : 2n+1), which lies inside the hiding cone (a = b = 1/(2n + 1) ≤ 1), where 𝒜n = ℋn = (2n + 1)√((2n + 1)2+2). Finally, ℋn - 𝒜n ≤ (nδ+)/(‖ v ‖2) ≤ (n(| A | + | C |))/(√(A2+C2+D2)) ≤ √2 n (attained at (1,1,0)), and for | A | + | C | > 0, using ‖ v ‖2 𝒜n ≥ ℬ ≥ M(| A | + | C |) and δ+ ≤ | A | + | C | gives

(ℋn(v) - 𝒜n(v))/(𝒜n(v)) ≤ (nδ+)/(M(| A | + | C |)) ≤ n/M = n/((n+1)2),

with equality at every horizontal view (D = 0). □

B.3 Isotropic Mean Shadow, Silver-Ratio Leading Term, and O(n-5) Tail Brackets

Averaging 𝒜n(v) over the unit sphere S2 with respect to the uniform probability measure yields the isotropic mean shadow 𝒜n. Because Ωn is non-convex for n ≥ 1, its total exposed surface area is Area2(∂Ωn) = 2(2n + 1)2 + 4∑j=0n (2j+1) = 12n2 + 16n + 6, and Cauchy's surface average Area2(∂Ωn)/4 = 3n2 + 4n + 3/2 counts ray-piercing multiplicity, exceeding 𝒜n by (2-√2)n2 - χ n + 2/π ∑i=1n I(2i) > 0. By contrast, the convex-hull mean is ℋn = ((2n + 1)2+2(2n + 1))/2 + n(n+1)(√2-1) = (1+√2)n2 + (3+√2)n + 3/2 (using 𝔼| A | = 1/2 and 𝔼 max(| A |, | D |) = 1/(√2) on S2), which shares the quadratic leading term (1+√2)n2 of 𝒜n and differs only at linear order (√2-1-χ)n + 2/π ∑i=1n I(2i).

Theorem B.3 (Exact Isotropic Mean Shadow, Closed Constant C*, and Tail Asymptotics). For every integer n ≥ 0, the isotropic mean actual shadow area of Ωn is

𝒜n = (1 + √2)n2 + (4 + χ)n + 3/2 - 2/π ∑i=1n I(2i) = (1 + √2)n2 + (4 + χ)n + C* + εn, (B.9)

where I(β) = ∫β∞ ((p-β)2)/(p(p2+1)(p2+2)) dp (β ≥ 1), εn = 2/π ∑i > n I(2i), and the exact constants are

χ = (ln 3)/(2π) - 1 + (2√2)/π arctan √2, C* = 4/3 + (ln 2)/(2π) + (ζ(3) + Li3(e-√2 π) - 2Li3(e-π))/π3. (B.10)

To 50 rounded decimal places:

  • 4 + χ = 4.03493671467655295067734553073116140549324173926145,
  • C* = 1.47999653090760042442100878465142982068794763080876,
  • 𝒜1 = 7.93754352744295952118684565150792784372616830715013 (vs. facet average Area2(∂Ω1)/4 = 8.5).

Moreover, for every integer n ≥ 1, writing P3(n) = 1/(24π n) - 1/(48π n2) + 7/(1440π n3) and P4(n) = P3(n) + 1/(320π n4), the tail εn satisfies the global two-sided bounds

P4(n) - 1/(288π n5) < εn < P4(n) - 19/(16128π n5) (n ≥ 1). (B.11)

The shifted sequence (n + τ)εn strictly increases to 1/(24π) for all n ≥ 1 for every real shift τ ≤ 1/2 (with τ = 1/2 the largest shift having this property: 24π(n + 1/2)εn = 1 - 2/(15n2) + O(n-3)). By contrast, the third-order remainder extractor Gn = n3(εn - 1/(24π n) + 1/(48π n2)) is not monotone, satisfying G1 < G2 > G3.

Proof.

  1. Equal Quotient-Chart Integrals and Exact Finite Mean Assembly: On the sorted domain 0 ≤ a ≤ b, define Wβ(a, b) = a/(2b) (ν - β)+2 for β ≥ 1 (with ν = b/max(1,a) and Wβ = 0 at b = 0), so that ∑j=0n-1 Cj(a,b) = n W1(a,b) - ∑i=1n W2i(a,b). Combining the spherical surface Jacobian (1+a2+b2)-3/2 with the orthogonal area factor (1+a2+b2)-1/2 and the 16-fold sign/swap symmetry on S2 gives the averaging operator 4/π ∫0∞ ∫0b (·)/((1+a2+b2)2) da db. Splitting at a = 1 by Tonelli's theorem:
  • Chart 1 (0 ≤ a ≤ 1): Here ν = b ≥ β ≥ 1, leaving the full inner range 0 ≤ a ≤ 1:
∫β∞ ((b - β)2)/(2b) [∫01 a/((1 + b2 + a2)2) da]db = ∫β∞ ((b - β)2)/(4b(b2 + 1)(b2 + 2)) db = 1/4 I(β).
  • Chart 2 (a ≥ 1): Substitute b = aν (ν ≥ β, db = a dν, a/(2b) = 1/(2ν)):
∫β∞ ((ν - β)2)/(2ν) [∫1∞ a/((1 + (1 + ν2)a2)2) da]dν = ∫β∞ ((ν - β)2)/(4ν(ν2 + 1)(ν2 + 2)) dν = 1/4 I(β).

Adding the two equal chart contributions and multiplying by 4/π proves 𝔼S2[(| D |)/(‖ v ‖2) Wβ(a,b)] = 2/π I(β). Partial fractions ((p-β)2)/(p(p2+1)(p2+2)) = β2/(2p) + ((1-β2)p - 2β)/(p2+1) + ((β2/2 - 1)p + 2β)/(p2+2) integrate directly to

I(β) = (β2 - 1)/2 ln(1 + β-2) + (2 - β2)/4 ln(1 + 2β-2) - 2β arctan(β-1) + √2 β arctan((√2)/β),

giving I(1) = (ln 3)/4 - π/2 + √2 arctan √2 and χ = 2/π I(1). Since 𝔼S2[δ+/‖ v ‖2] = 2(𝔼 max(| A |, | D |) - 𝔼| D |) = √2 - 1, subtracting the mean gap 𝔼[ℋn - 𝒜n] = n(√2-1) - n2/π I(1) + 2/π ∑i=1n I(2i) from ℋn = (1+√2)n2 + (3+√2)n + 3/2 establishes (B.9).

  1. Closed Form for C* via the Weighted Log-Sinh Integral: For β ≥ 2, expanding 1/((1+p-2)(1+2p-2)) = ∑k=0∞ (-1)k(2k+1-1)p-2k and integrating termwise via ∫β∞ ((p-β)2)/p2j+3 dp = (2β-2j)/((2j)(2j+1)(2j+2)) gives the absolutely convergent series
I(β) = ∑j=1∞ ((-1)j-1 2(2j - 1)β-2j)/((2j)(2j+1)(2j+2)) (β ≥ 2).

Setting β = 2i and summing over i ≥ 1 gives R∞ = 2/π ∑i=1∞ I(2i) = 4/π ∑j=1∞ ((-1)j-1(2j-1)ζ(2j))/((2j)(2j+1)(2j+2)4j) = 4/π [F(1/√2) - F(1/2)], where F(z) = ∑j=1∞ ((-1)j-1ζ(2j)z2j)/((2j)(2j+1)(2j+2)) = 1/2 ∫01 (1-ψ) ln((sinh(π z ψ))/(π z ψ)) dψ by the Euler product for sinh(π z ψ)/(π z ψ). Writing ln((sinh(kψ))/(kψ)) = kψ - ln(2k) - ln ψ + ln(1 - e-2kψ) and integrating against (1-ψ) yields

L(k) = ∫01 (1 - ψ) ln((sinh(kψ))/(kψ))dψ = k/6 - 1/2 ln(2k) + 3/4 - (ζ(2))/(2k) + (ζ(3) - Li3(e-2k))/(4k2).

Evaluating 2/π [L(π/√2) - L(π/2)] with k1 = π/√2, k2 = π/2, and ζ(2) = π2/6 combines the π-linear and 1/π terms to (k1-k2)/6 - (ζ(2))/2 (k1-1-k2-1) = (π(√2-1))/12 - (π(√2-2))/12 = π/12, which multiplies by 2/π to contribute 1/6 to R∞ (while -1/2 ln(2k1) + 1/2 ln(2k2) = -(ln 2)/4 contributes -(ln 2)/(2π)), yielding R∞ = 1/6 - (ln 2)/(2π) + (2Li3(e-π) - Li3(e-√2 π) - ζ(3))/π3 and C* = 3/2 - R∞ in (B.10).

  1. Signed Power-Tail Bounds for (B.11): Since the alternating terms of I(2i) decrease in magnitude for i ≥ 1, 1/(48i2) - 1/(320i4) < I(2i) < 1/(48i2) - 1/(320i4) + 1/(1536i6), so with the Hurwitz zeta function [20, §25.11] ζ(k, n+1) = ∑i > n i-k,
1/(24π) ζ(2, n+1) - 1/(160π) ζ(4, n+1) < εn < 1/(24π) ζ(2, n+1) - 1/(160π) ζ(4, n+1) + 1/(768π) ζ(6, n+1).

In the Laplace representation ζ(k, n+1) = 1/((k-1)!) ∫0∞ (ψk-1e-nψ)/(eψ - 1) dψ, the partial-fraction series 1/(eψ - 1) = 1/ψ - 1/2 + 2ψ∑j=1∞ 1/(ψ2 + 4π2 j2) has alternating geometric truncations for all ψ > 0:

1/ψ - 1/2 + ψ/12 - ψ3/720 < 1/(eψ - 1) < 1/ψ - 1/2 + ψ/12 - ψ3/720 + ψ5/30240.

Integrating against ψk-1e-nψ/(k-1)! yields 1/n - 1/(2n2) + 1/(6n3) - 1/(30n5) < ζ(2, n+1) < 1/n - 1/(2n2) + 1/(6n3) - 1/(30n5) + 1/(42n7), 1/(3n3) - 1/(2n4) + 1/(3n5) - 1/(6n7) < ζ(4, n+1) < 1/(3n3) - 1/(2n4) + 1/(3n5), and ζ(6, n+1) < ∫n∞ ω-6 dω = 1/(5n5). Substituting these signed brackets gives εn > P4(n) - 1/(288π n5) and εn < P4(n) - 37/(11520π n5) + 41/(20160π n7) ≤ P4(n) - 19/(16128π n5) for all n ≥ 1 (since 37/11520 - 41/20160 = 19/16128), proving (B.11).

  1. Exact Signed-Tail Proof of Shifted Monotonicity: Put Qi = i2 I(2i). Substituting p = βψ in I(β) gives
β2 I(β) = ∫1∞ ((ψ - 1)2)/(ψ(ψ2 + β-2)(ψ2 + 2β-2)) dψ,

whose integrand is positive, strictly increases with β > 0, and is dominated by (ψ-1)2/ψ5 (of integral 1/12). By dominated convergence, Qi = 1/4 (2i)2 I(2i) strictly increases with i ≥ 1 toward limi→∞ Qi = 1/48. For M = n + 1 and any shift τ ≤ 1/2, taking the forward difference gives

π/2 [(n + 1 + τ)εn+1 - (n + τ)εn] = ∑i > M Qi/i2 - ((M + τ - 1)QM)/M2.

Because Qi > QM for every i > M and the trapezoidal comparison on [i-1, i] (i ≥ M+1) gives the strict lower bound 1/M = ∫M∞ ω-2 dω < 1/(2M2) + ∑i > M i-2 ⇒ ∑i > M i-2 > 1/M - 1/(2M2) (whereas the midpoint rule on [i-1/2, i+1/2] gives the upper bound ∑i > M i-2 < 1/(M+1/2)), we obtain

π/2 [(n + 1 + τ)εn+1 - (n + τ)εn] > QM(1/M - 1/(2M2) - (M + τ - 1)/M2) = (QM(1/2 - τ))/M2 ≥ 0 (τ ≤ 1/2).

Since (n+τ)εn → 1/(24π), strict increase implies (n+τ)εn < 1/(24π) for all n ≥ 1 whenever τ ≤ 1/2 (yielding the sharp enclosure 0 < εn < 1/(24π(n+1/2)) at τ = 1/2). Conversely, for any τ > 1/2, the expansion (n+τ)εn = 1/(24π) + (τ - 1/2)/(24π n) + O(n-2) lies strictly above its limit 1/(24π) for all large n, so it cannot be globally increasing. (Note that differentiating the asymptotic truncation (ω+τ)P4(ω) cannot substitute for this exact discrete proof: at τ = 1/2, Φ(ω) = 24π(ω+1/2)P4(ω) = 1 - 2/(15ω2) + 2/(15ω3) + 3/(80ω4) actually has Φ'(1) = -17/60 < 0 and Φ(1) = 83/80 > 757/768 = Φ(2), whereas the exact sequence (n+1/2)εn strictly increases for all n ≥ 1.) Finally, high-precision evaluation of I(2), I(4), I(6) and C* from (B.10) gives G1 ≈ 0.00176526 < G2 ≈ 0.00180513 > G3 ≈ 0.00176694. □

Remark B.3a (Slab-Center Sets, the Constant C*, and Working Precision). In §B.1, let 𝒵ntop = {(p1,p2,n+1/2-ℓ) : 0≤ℓ≤ n, | p1|,| p2|≤ n-ℓ} denote the cube-center set of the top-widest solid stack Ωn. The vertical affine reflection R3aff(p1,p2,p3) = (p1,p2,n+1-p3) maps 𝒵ntop onto the bottom-widest slab-center set 𝒵nbot = {(p1,p2,ℓ+1/2)}, and the translate 𝒵nbot - (0,0,1/2) is exactly the lattice point stack Sn of Proposition 6.5. Throughout Theorem B.3 and Figure B.1(D), the isotropic mean constant C* and the tail sum R∞ = 3/2 - C* are defined by the closed trilogarithm formula (B.10):

C* = 4/3 + (ln 2)/(2π) + (ζ(3) + Li3(e-√2 π) - 2Li3(e-π))/π3 ≈ 1.47999653090760042442100878465142982068794763080876

(the numerical plots in Figure B.1(D) evaluate εn = R∞ - 2/π ∑i=1n I(2i) at 90-digit working precision across 1 ≤ n ≤ 15, while monotonicity for all n when τ ≤ 1/2 rests on the proof of Theorem B.3(4)).

Figure B.1: Solid Odd-Square Stack Quotient Shadow, Square Hiding Cone, and Shifted Tail Monotonicity
Figure B.1 (Three-dimensional odd-square stack Ωn shadows, square hiding cone, and shifted tail monotonicity): A. Quotient shadow P0∪ P1 in the plane p3=0 for the negative-control witness n=1, v=(1/3,3/2,1) (Theorem B.1), showing the top-slab hexagon P1, base protrusion P0∖ P1, and corner-notch gap separating 𝒜̃1(1/3,3/2)=541/36 from the projected convex hull 31/2 and naive extrapolation 15. B. Orthographic shadow area along the diagonal a=b=| A/D|=| C/D| at n=2 (2n+1=5, Theorem B.2), showing ℋ2=𝒜2 on a≤ 1, isometric value 35/√3 at a=1, and global maximum 5√27 at a=1/5. C. Slope-plane map (a,b)=(A/D,C/D) of the square hiding cone max(| a|,| b|)≤ 1 (Eq. B.5). D. Shifted tail product 24π(n+τ)εn from Eq. (B.10), increasing strictly to 1 for τ=0 and at the sharp critical shift τ=1/2, while decreasing from above 1 across 1 ≤ n ≤ 15 at τ=0.65 > 1/2.

Appendix C. Summary of Coordinates, Notation, and Static Geometric Constructions

C.1 Local Symbol and Coordinate Summary

All notation used in §§1–8 and Appendices A–B is defined locally at first appearance and collected below; the letters follow the notation register [7]:

Mathematical ObjectSymbol & Defining FormulaDomain & RoleSection & Equation
Triangular polynomial & integersT(z) = z(z+1)/2, Tn = T(n)z ∈ ℂ, n ∈ ℤ§1.1 (1.1), §6.2 (6.4)–(6.5)
Cell address, chamber & offsetN(r,c) = Tm-1 + r, m = r+c-1, V = r-c(r,c) ∈ ℤ≥ 12; 2N = m2+V+1§1.1 (1.1)–(1.6)
Centered odd legs & anti-diagonal pairr̃ = 2r-1 = m+V, c̃ = 2c-1 = m-V; (k, p-k)r̃+c̃=2m; k↔ p-k induces Rmid§1.1 (1.2), §3.2 (3.6)
Affine path incrementsL = r0+c0-2, D = Δr+ΔcDegree-≤ 2 polynomial along (r0,c0)+k(Δr,Δc)§2.1 (2.1)
Chamber phases & root displacementκrow=r/m, κmid=r̃/(2m), κroot=ε+1/2, κleg=r/(m+1)ε = √(2N) - m ∈ (-1/2, 1/2)§3.1 (3.1), App. A.3 (A.8)
Three chamber reflectionsRrow(N)=m2-N, Rmid(N)=m2+1-N, Rroot(N)=((2m-√(2N))2)/2Induced by κ ↦ 1-κ at fixed m§3.1 Thm. 3.1, (3.2)
Scaled-square floor-branch parametersNk = A(Qk)2, P = round(Q√(2A)), δ = P2-2AQ2, ℓk = ⌊1/2+kε⌋mk = Pk+ℓk; ℓk = 0 ⇔ -P ≤ δ k < P§4.1 (4.1)–(4.3)
Centered rational carrier parameters𝒞(r,c)=(V/m,(√(r̃c̃))/m), (r̃,c̃)=λ(U2,W2), m0=(U2+W2)/2U,W,λ ≥ 1 odd, gcd(U,W)=1, m0 ≡ 1 (mod 4)§5.1 (5.1)–(5.5), App. A
Address sheetpN = (-r/(m+1), 1/(m+1)), p̃N = ((r-m+1)/(m+1), 1/(m+1))Chamber-m row at height p2 = 1/(m+1) in Q = [-1,0]×[0,1]§§6.1–6.3 (6.1)–(6.7)
Discrete stack Sn & solid stack ΩnSn ⊂ ℤ3 (| Sn| = C(2n+3, 3)); Ωn ⊂ ℝ3, viewing vector v=(A,C,D)Fiber depth wn(a,b); orthographic shadow 𝒜n(v)§6.5 (6.10)–(6.11), App. B
Shells, Jordan totient & residualRd(k,h)=d2 k-d+h, Jk(d)=dk∏p | d(1-p-k), Q(N)=N/gcd(N,m!)d=m/gcd(m,r), h=r/gcd(m,r); Q(N)∈{1}∪ℙ§§8.2–8.4 Thms. 8.1–8.3
Rising product & reduced quotientΠn(y)=∏j=1n(1+jy)=∑k ek(n) yk, Tn/n!=an/bn, gn=gcd(e1(n),…,en-1(n))gn=an∈{1,n+1} for n≥ 2§8.4 Thm. 8.3(5)
Entire triangular coefficient liftEEi(z)=∑n≥ 1zn/(n n!), Θz = zd/(dz)G+(z)=(z(z+2))/2 ez, G-(z)=z2/2 ez on ℂ§8.4a (8.1)–(8.4)

C.2 Rail Definitions and Static Geometric Summary

  1. Rails (§2.2): The five reference rays are the first-row edge (1,k) (orange; N = Tk-1+1, r̃=1), the first-column edge (k,1) (yellow; N = Tk, r̃=2m-1), the diagonal (k,k) (purple; V=0, N=(m2+1)/2), the superdiagonal (k,k+1) (teal; V=-1, N=2k2), and the subdiagonal (k+1,k) (pink; V=+1, N=2k2+1). The colors are figure keys only.
  2. Ideal Gapless Solid Ωn (gap = 0) vs. Schematic Display-Gap Render (gap > 0): Theorems B.1–B.3 apply strictly to the gapless union of unit-height slabs Ωn (gap = 0, where each quotient shadow Pj = [-Rj,Rj]2 + [-(j+1)(a,b), -j(a,b)] is a convex hexagon of area (2j+1)2 + (2j+1)(a+b), smaller by ab than its axis-aligned bounding box). Shrinking displayed cells to side (1-gap)×pitch (gap = 0.14 in Figure 11(D)) is solely a visual aid for exposing interior layers and alters the projected silhouette.
Structural DomainGeometric Map & InvariantDistinguishing Witness & Scope
Quadrant, Wedge & Rails (§§1–2)(r,c)↦ N=Tm-1+r ↦ (m,V) straightens quadratic curves N(m) into affine lines in (m,r̃) and preserves period-4 parity TV+1 mod 2.N = 25 ↔ (4,4), m = 7, V = 0; parity cycle (1,1,0,0) on V = 0,1,2,3.
Reflections & Floor Drift (§§3–4)Fixed-m phase flip κ↦ 1-κ separates Rrow, Rmid, Rroot; ℓk = ⌊1/2+kε⌋ selects polynomial branches of Nk = A(Qk)2.N = 5 (m=3): 4 < 23-6√10 < 5; 181/128 jumps at k=26 (m=4707) before 412 (k=32).
Semicircle Budgets (§5 & App. A)Centered 𝒞(r,c)=(+V/m,√(m2-V2)/m) vs. leg-sum pleg(r,c) under budget m ≤ M.Least nontrivial m0=5 at (5,1)↦(4/5,3/5); target (3,4) exception (3,5) at M=17 < Mcut=37.
Address Sheet & Balanced Split (§6)Row map pN=(-r/(m+1),1/(m+1)) into hull 𝒦 (Rin=(√2-1)/2); 8-sector ℤ2⋊O(2,ℤ) congruence and the finite reciprocal identity 𝒮LΠ△,L(-1/8)=1.(2,1)↦ N=3,14,34 on p1=-2/3; m=19↦φ(20)=8; L=1: 𝒮1=8/9, Π△,1=9/8.
3D Point Stack Sn vs. Solid Ωn (§6.5 & App. B)Discrete Sn along (1,1,1)/√3 has fiber depth wn(a,b) and atomic PMFs → 2ω, 6𝔣(ω); solid Ωn (gap = 0) has hiding cone max(| A|,| C|)≤| D|.n=2 stack size | S2|=35; n=1, v=(1/3,3/2,1) gives 𝒜̃1=541/36 < 31/2 and 𝒜1=541/66.
Figure 11: Static Geometric Summary — Straightened Parity Wedge, Address-Sheet Primitive Columns & Hull Incircle, Oblique Fiber Depth, and Solid vs. Display-Gap Domain
Figure 11 (Static geometric summary across four panels): A. Admissible parity wedge (m,V) straightening the five rails into affine rays. B. Address-sheet primitive columns in Q (φ(20)=8 primitive cells in red at m=19) and hull incircle radius Rin(𝒦) = (√2-1)/2. C. Oblique fiber multiplicity w2(a,b) = 1 + ⌊(2-max(a,b))/2⌋ summing to | S2| = 35. D. Comparison between the gapless solid domain (gap = 0) of Appendix B and a schematic display-gap render (gap = 0.14).

References

  1. T. Nagell, The Diophantine equation x2 + 7 = 2n, Arkiv för Matematik 4 (1961), 185–187. https://doi.org/10.1007/BF02592006
  2. M.-G. Leu and G.-W. Li, The Diophantine equation 2x2 + 1 = 3n, Proceedings of the American Mathematical Society 131 (2003), 3643–3645. https://doi.org/10.1090/S0002-9939-03-07212-5
  3. J. H. E. Cohn, Perfect Pell powers, Glasgow Mathematical Journal 38 (1996), 19–20. https://doi.org/10.1017/S0017089500031207
  4. J. Huckstead, Triangular-n Postmaster: From the Triangular Primitive to the Completed Source Frontier, Volume I, Version 6.3, Cerebral Graphix / Zenodo, October 2026. https://doi.org/10.5281/zenodo.23510962
  5. J. Huckstead, Reading the Fold: A Global Newton Series for the Riemann Xi Function, Version 1.3, Cerebral Graphix / Zenodo, October 2026. https://doi.org/10.5281/zenodo.23511096
  6. H. B. Curry and I. J. Schoenberg, On Pólya frequency functions IV: the fundamental spline functions and their limits, Journal d'Analyse Mathématique 17 (1966), 71–107. https://doi.org/10.1007/BF02788653
  7. J. Huckstead, Triangular Index Notation Register and Modular Coefficient-State Notes (Version 1.1) and Cross-Project Symbol Concordance (Version 0.1), Cerebral Graphix working manuscripts, October 2026.
  8. J. Huckstead, Cayley–Casimir Bridge to Exact Semidiscrete Flux Scattering, Version 0.8, Cerebral Graphix / Zenodo, October 2026. https://doi.org/10.5281/zenodo.23307973
  9. J. Huckstead, The Notebook Grid: Exact Addresses, Three Reflections, Pell Branches, and Rational Carriers, Version 1.5 (earlier working-paper deposit), Zenodo, October 2026. https://doi.org/10.5281/zenodo.23459201 (Historical concept record: https://doi.org/10.5281/zenodo.23459200)
  10. J. Huckstead, The Phase-Complement Collar: An Exact Odd-Cell Bridge from Finite-Observer Phase Accounting to the Triangular / Casimir Coordinate (Versions 0.1–0.2), Cerebral Graphix working manuscripts, August 2026.
  11. J. Huckstead, Divisor-Band Bijections and Totient Shells in the Triangular-Fractional Grid, Version 5, Cerebral Graphix / Zenodo, July 2026. https://doi.org/10.5281/zenodo.21440877
  12. J. Huckstead, The Chamber Lift Series: Farey Chambers, Jordan Denominator-Shell Collapse, Zeta-Exponent Three Ledgers, and the Triangular Factorial-Residual Calculus (Versions 4–9), Cerebral Graphix working manuscripts, June–August 2026.
  13. K. Dickman, On the frequency of numbers containing prime factors of a certain relative magnitude, Arkiv för Matematik, Astronomi och Fysik 22A (1930), 1–14; N. G. de Bruijn, On the number of positive integers ≤ x and free of prime factors > y, Indagationes Mathematicae 13 (1951), 50–60.
  14. J. Huckstead, Prime Numerators and Shell Stratification in the Triangular-Fractional Grid (Version 0.1), Cerebral Graphix working manuscript, July 2026.
  15. J. Huckstead, Diagonal-Center Sum–Difference Geometry and Arithmetic (Version 1.1) and Exponential-Integral Triangular Coefficient Investigation (Version 0.1), Cerebral Graphix working manuscripts, 2026.
  16. J. Huckstead, Mersenne-Triangular Stations, Simplex-Factorial Closure, and Triangular Canonical Product / Brouncker Continuant Notes, Cerebral Graphix working manuscripts, August 2026.
  17. J. Huckstead, Bounded Address-Sheet Totient Partitions, Staircase Fixed Points, and Chamber-Offset Entropy Notes, Cerebral Graphix working manuscripts, October 2026.
  18. J. Huckstead, Three-Dimensional Cone/Cube Ratios, Orthochronous Double-Pell Boosts, Tetrahedral-Hexagonal Projections, and Projective Cayley / Möbius Notes, Cerebral Graphix working manuscripts, March–August 2026.
  19. J. Huckstead, Exact-Rational Euler-Transform Enclosure for the Telescoping Logarithmic Series S+, Cerebral Graphix working manuscript, October 2026.
  20. F. W. J. Olver, A. B. Olde Daalhuis, D. W. Lozier, B. I. Schneider, R. F. Boisvert, C. W. Clark, B. R. Miller, B. V. Saunders, H. S. Cohl, and M. A. McClain (eds.), NIST Digital Library of Mathematical Functions, §6.2 (exponential integrals Ei, Ein) and §25.11 (Hurwitz zeta function). https://dlmf.nist.gov/

Acknowledgments and Computational Reproducibility

Author & Collaboration Note. Authored by Jeffery Huckstead (Cerebral Graphix, cerebralgraphix.com, ORCID 0009-0007-0234-2177). AI-assisted mathematical drafting, symbolic reorganization, proof checking, and computational verification tools were used under the author's direction during the preparation of Version 4.1; responsibility for all mathematical definitions, scope boundaries, and conclusions rests with the author.

Optional Computational & Interactive Supplements. Every theorem and proof in this paper is self-contained in print. For readers wishing to replay finite-range algebraic, rational, decimal, and sieve checks, the supplementary archive provides the standalone Python 3 standard-library script check_notebook_grid_v4_1.py alongside an optional offline single-file interactive HTML reader (Notebook_Grid_v4_1_Interactive_Reader.html).