Continuum Core Model: Deriving $Rₚ$

Continuum Core Model: Deriving Rₚ: Agreement with data is a result here, never a reason. The one identification the continuum chain needs, declared and gated. LOCK: fix the identification rule L_q=λ_C and η=L_q/λ_C=1 in canon_lock (SSOT and No-Tuning). Grade [F] forced.

Agreement with data is a result here, never a reason. The one identification the continuum chain needs, declared and gated. LOCK: fix the identification rule L_q=λ_C and η=L_q/λ_C=1 in canon_lock (SSOT and No-Tuning).

Purpose (locking the outputs)

The purpose of this chapter is to construct a continuum (continuous-approximation) core model that determines the core radius Rₚ using internal definitions only, and to derive and lock the following outputs.

  1. The defining expression of the core radius, Rₚ:=R[core state] (including definition and regime).
  2. A formal statement of the core stability condition (locked as one of: rest condition / equilibrium condition / switch condition).
  3. The dimensionless conclusion for the length-selection ratio:
    \begin{equation} \frac{R_p}{L_q}=\alpha=\frac{2}{\pi}, \end{equation}

    where L_q is the core selection length (an internal length scale) and α is the rectification constant whose unique source is fixed in Chapter 5.

  4. The insertion points for downstream derivations (event rates, mass/force scales) and an explicit statement of forbidden retroactive justification.

This chapter does not introduce Rₚ by appealing to external texts. The only admissible inputs are internal definitions and internal conventions locked in canon_lock/analysis_lock. If a comparison to external numerics is required, that comparison is treated only as a validation (Gate) item; it is never used as a premise of the derivation.

Inputs (LOCK) and regime (scope)

The derivation of Rₚ is carried out only when the following input items are locked.

Substrate (Stone) regime and VP axioms

This chapter assumes the VP axiom set locked in §3.1 (infinite rigidity / impenetrability / plenum / local rule / adjacency). In particular, we fix the following two conditions as regime conditions.

  1. Stone regime: the volume invariance and impenetrability of VP are maintained as global constraints.
  2. Plenum regime: no “empty space” is introduced as an independent degree of freedom; “deficits/throats/gaps” at the core boundary appear only as derived structural quantities of placement and adjacency.

If these regime conditions are not locked, the continuum approximation in this chapter loses its scope, and its conclusions are judged INCONCLUSIVE.

Canonical cell (Anchor Cell) and length scales

The canonical cell is locked to CELL-CUBE, and the representative cell length D_anch is locked as the edge length (edge). The half-length scale must be fixed by the following derived definition.

\begin{equation} r_0:=\frac{D_{\mathrm{anch}}}{2}. \end{equation}

The continuum coordinate R in this chapter is used as the radial distance inside the canonical cell; the unit and coordinate convention of R are locked by the coordinate convention in analysis_lock.

Rectification constant α as an input

This chapter uses, as an input, the rectification constant whose unique source is fixed in §5.1:

\begin{equation} \alpha=\frac{2}{\pi}. \end{equation}

α is not re-derived in this chapter; it is inserted only into the dimensionless ratio conclusion (S06_goal_ratio) for core length selection.

Status of the core selection length L_q

L_q is an internal length scale defined as the “selection length” in the core model. The status of L_q must be locked as one of the following.

  1. CANON-INPUT: lock L_q as a canonical input in canon_lock (including meaning/unit/scope).
  2. CANON-DERIVED: define L_q as a canonical derived quantity and lock its derivation rule (from which structural quantity it is derived and how) in analysis_lock.

If the status of L_q is not locked, the dimensionless conclusion Rₚ/L_q is itself undefined, and the derivation in this chapter does not stand.

Internal definition of the continuum core model (a continuous approximation of geo-structure)

In this chapter, “continuum” does not mean importing an external field theory. It means an internal procedure that aggregates discrete structural quantities of VP placement into functions of a radial variable R. The continuous approximation is fixed by the following rules.

Definition of the core region and the core boundary

The core is a region defined with respect to the domain center x_c. The core boundary is defined by the following switch condition.

  1. For radius R, define an indicator χ_core(R)∈0,1 that judges whether the discrete structural index inside radius R satisfies a “rigidity/transfer” condition.
  2. Define the core radius Rₚ as a transition point of χ_core(R). The discrete/continuous decision rule for the transition point (first transition / threshold transition / stable transition) must be locked in analysis_lock.

Therefore Rₚ is not an “arbitrary geometric length” but an operationally defined length determined by a transition of a discrete structural indicator.

Standard form of radial aggregates

A radial aggregate S(R) at radius R is defined only in the following standard form.

\begin{equation} S(R) := \mathcal{A}\Bigl(\{\Omega_i\}_{i\in\mathcal{V}(R)},\ \mathcal{G}_c(R),\ \mathcal{P}(R)\Bigr), \end{equation}

where

All continuum expressions in this chapter must be reducible to the form (S06_radial_aggregate); if the reduction rule is not locked, the continuum expressions are undefined.

Internal definition of the stability condition (selection radius as an intersection of two scaling laws)

This chapter adopts the internal principle that the core radius Rₚ is selected as the intersection of two radial scaling laws. This is not an external justification; it is implemented by defining the scaling laws of the “aggregate cost” at the core boundary under the VP Stone and plenum constraints, and then locking their intersection as the selection rule.

Where the two aggregate cost functions are defined

For radius R, define the following two costs (or pressure-like scalar aggregates):

\begin{equation} \Pi_{4}(R) := \frac{C_{4}}{R^{4}}, \qquad \Pi_{5}(R) := \frac{C_{5}}{R^{5}}, \end{equation}

where C₄,C₅ are constants derived from internal aggregation (or locked coefficients), and the exponents R⁻⁴ and R⁻⁵ are derived and locked internally from domain/boundary/counting conventions. This overview fixes only the definition slot and role of (S06_two_scalings); the explicit derivation of C₄,C₅ and the internal derivation of the exponents are completed in later sections of this chapter.

Stability condition for the selection radius (core radius)

The selection (stability) condition for the core radius Rₚ must be locked as one of the following (exclusive choice; forbidden to swap depending on results).

  1. Intersection condition: Π₄(Rₚ)=Π₅(Rₚ).
  2. Rest condition: for some locked energy-like aggregate U(R), dU/dR|_(R=Rₚ)=0, and U is defined so that this condition is equivalent to the intersection in (S06_two_scalings).
  3. Switch condition: define a switch observable so that the transition point of a structural indicator χ_core(R) coincides with the intersection point of (S06_two_scalings).

The canonical choice of this chapter is the intersection condition. Under the intersection condition,

\begin{equation} \Pi_{4}(R_p)=\Pi_{5}(R_p) \ \Longrightarrow\ R_p=\frac{C_{5}}{C_{4}}. \end{equation}

Therefore the dimensionless ratio conclusion (S06_goal_ratio) is completed once the internal derivation of C₅/C₄ is completed.

6.5 Allowed conclusion forms in this chapter (implementation of the ban on external justification)

The conclusions of this chapter are permitted only in the following forms.

  1. Derivations of dimensionless ratios such as Rₚ/L_q (including insertion of the rectification constant α).
  2. Internal length-selection formulas of the form Rₚ=α L_q (only when the status of L_q is locked).
  3. Numerical values may be presented only when REALIZATION or validation (Gate) items are locked; numerical agreement must not be used to justify a definition/axiom/rectification convention.

Therefore this chapter does not use external doctrines (equations of other theories, definitions of other constants, or external justifications) as grounds for the derivation of Rₚ. Correspondence with external texts is handled only in a separate “correspondence table / non-use scope declaration” section and does not affect the conclusion status of this chapter.

LOCK/Gate connections for this section (none if empty)

6.1 Linking L_q and λ_C (L_q=λ_C)

6.1.1 Purpose

This section fixes the two length symbols L_q and λ_C appearing in the continuum core model as internal definitions, and then confirms the identification rule

\begin{equation} L_q=\lambda_C \end{equation}

as a canonical lock (LOCK). The outputs of this section are (i) the definition of L_q, (ii) the definition of λ_C, (iii) the internal basis for the identification (axiom-based), and (iv) the substitution rules after the identification.

6.1.2 Definition: L_q (core selection length)

[D-6.1-1] Dimensionless radial coordinate

Lock the core center x_c inside the canonical cell and define the radial distance as

\begin{equation} R := \|\mathbf{x}-\mathbf{x}_c\|. \end{equation}

The coordinate system of R and the selection convention of x_c are locked in analysis_lock.

The radial description of the core model proceeds in a dimensionless coordinate ξ; define its normalization length as L_q:

\begin{equation} \xi := \frac{R}{L_q}. \end{equation}

[D-6.1-2] Meaning of L_q (selection length)

Define L_q as the unique length scale in the core model that simultaneously satisfies the following two conditions.

  1. (Normalization) It is the normalization reference length such that every radial aggregate of the core model can be expressed not as a function of R but only as a function of ξ=R/L_q.
  2. (Transition) It is the selection length such that the boundary transition of the core (a transition of the core indicator or a transition of a cost function) occurs at a single value ξ=ξₚ.

Therefore L_q is not “an arbitrary length,” but is defined as “the normalization length that fixes the core transition at a single dimensionless value.” This is an internal definition needed for the continuum core model to stand and does not use external-text justification.

6.1.3 Definition: λ_C (core phase-completion length)

[D-6.1-3] Core phase variable

Define a phase variable ψ(R) describing a radial unfolding of an internal core state as

\begin{equation} \psi:\ [0,\infty)\to\mathbb{R}, \qquad \psi(0)=0, \end{equation}

where ψ is an internal variable used to represent a “phase cycle” inside the core. The protocol for generating ψ (which aggregate is used as a phase) is locked in analysis_lock. This section uses only (i) the existence of ψ as an internal variable of the continuum core model and (ii) the fact that “cycle completion” is locked at 2π.

[D-6.1-4] Meaning of λ_C (first cycle-completion radius)

Define λ_C as the first positive length that satisfies

\begin{equation} \lambda_C := \inf\{\, R>0\ |\ \psi(R)=2\pi\,\}. \end{equation}

In (S06_lambdaC_def), 2π is the full-cycle canonical constant locked in Chapter 5 (the reference used by the rectification convention), and λ_C is the length scale at which that full cycle first completes in the radial unfolding of the core. Therefore λ_C is not an “external constant” but an “internal length defined by core phase completion.”

6.1.4 Basis for the identification (axiom-based)

The identification L_q=λ_C is not an extra assumption but a canonical selection (lock) needed to satisfy the global rules of this document (No-Tuning, SSOT, single-transition definition). The basis consists of the following four items.

(G1) Single-transition principle: the core boundary must be unique

The core radius Rₚ is defined as a single boundary indicating the transition “core interior→core exterior” (a transition of the core indicator or a transition of a cost function). If L_q and λ_C exist as independent lengths, the length that normalizes the core transition (L_q) and the length that defines phase completion (λ_C) can point to different transitions. Then the definition of the core boundary becomes duplicated, Rₚ cannot be fixed as a single length, and the core boundary becomes an ambiguous object depending on procedures. Therefore, to preserve the definition that the core boundary is unique, the two lengths must be identified so that they point to the same boundary.

(G2) SSOT principle: keep only one length for the same meaning

L_q and λ_C are defined by (S06_xi_def) and (S06_lambdaC_def), respectively, but the physical meaning they indicate must converge to the same slot (the core selection length). In the core model, “radial scale selection” and “phase-cycle completion” both indicate the point where internal organization ends and a boundary forms. If two lengths that indicate the same slot are kept simultaneously, the same meaning is split across multiple symbols, violating SSOT. To satisfy SSOT, the two symbols must be unified into one canonical entry; the unification rule is (S06_Lq_eq_lC).

(G3) No-Tuning principle: do not leave L_q/λ_C as a tunable degree of freedom

If L_q and λ_C are treated as independent inputs, an additional dimensionless degree of freedom appears:

\begin{equation} \eta := \frac{L_q}{\lambda_C}. \end{equation}

This ratio η can directly enter the core-radius ratios Rₚ/L_q, Rₚ/λ_C and subsequent derivations (event rate/mass/force), opening a path to adjust η to match outcomes. However, this document forbids adjusting degrees of freedom after seeing results (No-Tuning). To close the core selection without introducing a free knob, η must be fixed by a lock. This section declares

\begin{equation} \eta := 1, \end{equation}

which is equivalent to (S06_Lq_eq_lC).

(G4) Compatibility with the unique source of rectification constants: do not introduce a new universal constant beyond π

In Chapter 5, the rectification constants α=2/π and δ=1/π² are locked with a unique source. If one separates L_q and λ_C, one needs an additional constant or an additional convention to determine η in (S06_eta_def). This document does not introduce a new universal constant beyond the rectification-constant system. Therefore the only canonical choice that removes redundancy in the core selection length is η=1, i.e., L_q=λ_C.

6.1.5 Substitution rules after the identification

After the identification (S06_Lq_eq_lC) is locked, fix the following reuse rules globally.

  1. L_q and λ_C are two notations for the same entry, and canon_lock records them as a single entry (single value / single meaning under one lock_id).
  2. In formula manipulations, any occurrence of L_q may be replaced by λ_C, and vice versa. The replacement is only a notational substitution following a defined identification; it is not a new derivation or a new premise.
  3. Within the same version, it is forbidden to treat L_q and λ_C as different values or to split them into different object assignments/geometric meanings.
  4. To release the identification or to change to a different ratio (η≠ 1), one must version-bump canon_lock and perform full re-derivation / full re-validation.

LOCK/Gate connections for this section (none if empty)

6.2 Balance of 1/R⁴ vs 1/R⁵ → Rₚ/L_q=2/π

Concept links: the ratio equals the rectification α=2/π of §5.1; the same balance reappears in the gravity cap (§17.4).

6.2.1 Radial aggregation geometry and standard symbols

Lock the core center as x_c and define the radial distance by

\begin{equation} R:=\|\mathbf{x}-\mathbf{x}_c\| \end{equation}

R is an aggregation coordinate inside the canonical cell (CELL-CUBE); the cell geometry is not replaced by a sphere. Define the radial aggregation surface (a level set) at radius R by

\begin{equation} \mathbb{S}_R := \{\mathbf{x}\ |\ \|\mathbf{x}-\mathbf{x}_c\|=R\}, \qquad A(R):=\mathrm{Area}(\mathbb{S}_R)=4\pi R^2. \end{equation}

Here A(R) is the area of the radial aggregation surface; it is a definition of the aggregation surface and is independent of the canonical cell definition (cube).

The core selection length L_q is the length scale locked in §6.1. In this section, we use L_q in two roles.

  1. Unit length for radial aggregation: the reference length for making R dimensionless.
  2. Unit patch length on the radial aggregation surface: the minimal linear scale used to partition the aggregation surface.

Fix the unit patch linear size as L_q and define the unit patch area as

\begin{equation} A_0 := L_q^2. \end{equation}

A₀ is the minimal aggregation unit on the radial aggregation surface and must not be changed after seeing results.

6.2.2 Insertion point of the direction rectification factor α

When a directional component is aggregated into a scalar in radial aggregation, we use the direction rectification factor α as an input from the unique source in §5.1:

\begin{equation} \alpha=\left\langle |\cos\theta| \right\rangle=\frac{2}{\pi}. \end{equation}

Here θ is an angular variable aggregated with the uniform measure over one full cycle [0,2π); the averaging convention is locked in §5.1. In this section, α is used only as the rectification factor meaning “the effective contribution of the radial component survives only by a factor α on average”; α itself is not re-derived here.

6.2.3 Definition and expansion of the 1/R⁴ collapse term Π₄(R)

This section defines the “collapse term” (a pressure-like constraint indicator) Π₄(R) at the core boundary by two rounds of geometric dilution and an inverse correction for rectification loss.

6.2.3.1 Geometric dilution 1: area dilution (number of patches)

Partition the aggregation surface S_R at radius R by the unit patch area A₀. Define the number of patches as

\begin{equation} N_A(R) := \frac{A(R)}{A_0}=\frac{4\pi R^2}{L_q^2}. \end{equation}

Therefore the area fraction of a single patch on the full aggregation surface (the area-dilution factor) is

\begin{equation} f_A(R) := \frac{1}{N_A(R)}=\frac{A_0}{A(R)}=\frac{L_q^2}{4\pi R^2}. \end{equation}

We fix f_A(R) as a purely geometric factor representing “the average share received by a unit patch” on the aggregation surface.

6.2.3.2 Geometric dilution 2: angular dilution (directional window to hit a fixed patch)

Define the degree of directional alignment toward a fixed patch (linear size L_q) as the “directional window” to hit that patch. At radius R, viewing a patch of linear size L_q on the unit sphere, define a representative polar-angle width by

\begin{equation} \Delta\vartheta(R) := \frac{L_q}{R} \end{equation}

and define the corresponding representative solid angle by

\begin{equation} \Delta\Omega(R) := \bigl(\Delta\vartheta(R)\bigr)^2 = \left(\frac{L_q}{R}\right)^2. \end{equation}

Normalize the full directional space by the total solid angle 4π. The angular dilution factor is then

\begin{equation} f_\Omega(R) := \frac{\Delta\Omega(R)}{4\pi} = \frac{1}{4\pi}\left(\frac{L_q}{R}\right)^2. \end{equation}

We fix f_Ω(R) as the geometric factor representing “the directional window required by a unit patch.”

6.2.3.3 Rectification-loss correction (inverse factor 1/α)

Even when the direction is aligned, the radial component is a signed projection so cancellations occur on average. Since the effective contribution of the radial component survives only by α ((S06_02_alpha)), an inverse correction factor 1/α is required to secure the same “effective contribution.” This correction is not an ad hoc fit but a direct consequence of the rectification convention (§5.1).

6.2.3.4 Definition of the collapse term Π₄(R)

Define the collapse term (pressure-like constraint indicator) Π₄(R) at the core boundary as

\begin{equation} \Pi_{4}(R) := \Pi_\star\, \frac{1}{\alpha}\, f_A(R)\, f_\Omega(R), \end{equation}

where Π_* is the “unit constraint strength” locked inside the regime (dimensionless or internal unit) and is a reference value independent of R. Π_* is a common factor that cancels in this section; its magnitude does not affect the ratio conclusion.

Substitute (S06_02_fA) and (S06_02_fOmega) into (S06_02_Pi4_def) to fully expand:

\begin{align} \Pi_{4}(R) &= \Pi_\star\, \frac{1}{\alpha}\, \left(\frac{L_q^2}{4\pi R^2}\right) \left(\frac{1}{4\pi}\left(\frac{L_q}{R}\right)^2\right) \notag\\ &= \Pi_\star\, \frac{1}{\alpha}\, \left(\frac{L_q^2}{4\pi R^2}\right) \left(\frac{L_q^2}{4\pi R^2}\right) \notag\\ &= \Pi_\star\, \frac{1}{\alpha}\, \frac{L_q^4}{(4\pi)^2 R^4} \notag\\ &= \Pi_\star\, \frac{1}{\alpha}\, \frac{L_q^4}{16\pi^2}\, \frac{1}{R^4}. \end{align}

Therefore the collapse term is fixed to the 1/R⁴ scaling, and its coefficient includes the inverse factor 1/α.

6.2.4 Definition and expansion of the 1/R⁵ rigidity term Π₅(R)

This section defines the “rigidity term” (a pressure-like constraint indicator) Π₅(R) in the substrate (Stone) regime by adding a radial-chain lock factor to the same two rounds of geometric dilution.

6.2.4.1 Radial-chain lock factor η_R(R)

Define the number of unit lengths L_q stacked along the radial direction up to radius R (the number of radial layers) by

\begin{equation} N_R(R) := \frac{R}{L_q}. \end{equation}

Under the local rule (§3.1), radial transfer is represented as a chain of unit layers, and the average contribution of the chain is fixed by a rectification factor inversely proportional to the layer count. Define this as the “radial-chain lock factor”

\begin{equation} \eta_R(R) := \frac{1}{N_R(R)}=\frac{L_q}{R}. \end{equation}

η_R(R) is a geometric factor determined by the definition (S06_02_NR) and must not be adjusted after seeing results.

6.2.4.2 Definition of the rigidity term Π₅(R)

Define the rigidity term Π₅(R) as

\begin{equation} \Pi_{5}(R) := \Pi_\star\, f_A(R)\, f_\Omega(R)\,\eta_R(R). \end{equation}

Unlike the collapse term, the rigidity term does not include 1/α. The attenuation of the rigidity term is determined not by “rectification loss of a directional projection” but only by “radial-chain locking.”

Substitute (S06_02_fA), (S06_02_fOmega), and (S06_02_etaR) into (S06_02_Pi5_def) to fully expand:

\begin{align} \Pi_{5}(R) &= \Pi_\star\, \left(\frac{L_q^2}{4\pi R^2}\right) \left(\frac{1}{4\pi}\left(\frac{L_q}{R}\right)^2\right) \left(\frac{L_q}{R}\right) \notag\\ &= \Pi_\star\, \left(\frac{L_q^2}{4\pi R^2}\right) \left(\frac{L_q^2}{4\pi R^2}\right) \left(\frac{L_q}{R}\right) \notag\\ &= \Pi_\star\, \frac{L_q^5}{(4\pi)^2 R^5} \notag\\ &= \Pi_\star\, \frac{L_q^5}{16\pi^2}\, \frac{1}{R^5}. \end{align}

Therefore the rigidity term is fixed to the 1/R⁵ scaling. It shares the same geometric coefficient (4π)⁻² with the collapse term and has one additional order of attenuation by the factor L_q/R.

6.2.5 Balance condition and the complete derivation of Rₚ/L_q=2/π

Define the core radius Rₚ as the transition point where the “collapse term” and the “rigidity term” balance. Lock the balance condition as

\begin{equation} \Pi_{4}(R_p)=\Pi_{5}(R_p). \end{equation}

Substitute (S06_02_Pi4_final) and (S06_02_Pi5_final) into (S06_02_balance_condition):

\begin{align} \Pi_\star\, \frac{1}{\alpha}\, \frac{L_q^4}{16\pi^2}\, \frac{1}{R_p^4} &= \Pi_\star\, \frac{L_q^5}{16\pi^2}\, \frac{1}{R_p^5}. \end{align}

Cancel the common factors Π_* and 16π² on both sides:

\begin{equation} \frac{1}{\alpha}\,L_q^4\,\frac{1}{R_p^4} = L_q^5\,\frac{1}{R_p^5}. \end{equation}

Divide both sides by L_q⁴/Rₚ⁴:

\begin{equation} \frac{1}{\alpha} = \frac{L_q}{R_p}. \end{equation}

Therefore the dimensionless ratio of the core radius is

\begin{equation} \frac{R_p}{L_q}=\alpha. \end{equation}

Since α is fixed in §5.1 as α=2/π,

\begin{equation} \frac{R_p}{L_q}=\frac{2}{\pi}. \end{equation}

Moreover, since L_q=λ_C is locked in §6.1, the following substitution is allowed within the same locked version:

\begin{equation} R_p=\frac{2}{\pi}\,L_q=\frac{2}{\pi}\,\lambda_C. \end{equation}

Equations (S06_02_Rp_over_Lq_final) and (S06_02_Rp_lambdaC) are fixed as the conclusions of this section.

6.2.6 Validity conditions (regime prerequisites) and handling violations

For the conclusion (S06_02_Rp_over_Lq_final) to hold, the following items must be locked.

  1. The aggregation convention that adopts L_q as the unit patch length must be locked (A₀=L_q²).
  2. The direction rectification factor α must be locked by the canonical convention in §5.1 (α=2/π).
  3. The radial-chain lock factor η_R(L_q/R) must be locked by the regime convention.

If any of the above items is not locked, or if they are mixed within the same output, the definitions (S06_02_Pi4_def)(S06_02_Pi5_def) collapse and the corresponding conclusion loses conclusion status.

LOCK/Gate connections for this section (none if empty)

6.3 Numeric substitution for Rₚ and a summary of invariants

6.3.1 Inputs (LOCK) and reference equations (starting point)

In this section, lock the following inputs and reference equations.

  1. Rectification constant:
    \begin{equation} \alpha=\frac{2}{\pi}. \end{equation}
  2. Length identification:
    \begin{equation} L_q=\lambda_C. \end{equation}
  3. Ratio conclusion for the core radius (derived in §6.2):
    \begin{equation} \frac{R_p}{L_q}=\alpha. \end{equation}

The three equations above are not re-derived here. This section combines them to (i) compute a numerical value of Rₚ, (ii) define and compute the geometric cross section σ_geom, and (iii) summarize the invariant 4/π.

6.3.2 Full expansion of the Rₚ formula

From (S06_03_Rp_over_Lq),

\begin{equation} \frac{R_p}{L_q}=\alpha \quad\Longrightarrow\quad R_p=\alpha\,L_q. \end{equation}

Substitute (S06_03_Lq_eq_lC) into (S06_03_Rp_alpha_Lq):

\begin{equation} R_p=\alpha\,\lambda_C. \end{equation}

Substitute (S06_03_alpha) into (S06_03_Rp_alpha_lC):

\begin{equation} R_p=\frac{2}{\pi}\,\lambda_C. \end{equation}

Therefore, the numerical substitution for Rₚ in this section is fixed as the procedure “evaluate (S06_03_Rp_final_form) using the locked value of λ_C.”

6.3.3 Substituting the locked value of λ_C → Rₚ=0.8412fm

λ_C is the “core phase-completion length” defined in §6.1. In this section, assume that the following value is locked by canon_lock:

\begin{equation} \lambda_C = 1.3213538700998668\ \mathrm{fm}. \end{equation}

Substitute (S06_03_lC_value_lock) into (S06_03_Rp_final_form):

\begin{align} R_p &=\frac{2}{\pi}\,\lambda_C =\frac{2}{\pi}\times 1.3213538700998668\ \mathrm{fm} \notag\\ &= 0.8412\ \mathrm{fm}. \end{align}

With the unit-conversion convention 1fm=10⁻¹⁵m,

\begin{equation} R_p = 0.8412\times 10^{-15}\ \mathrm{m} = 8.412\times 10^{-16}\ \mathrm{m}. \end{equation}

6.3.4 Back-calculation (consistency): reconstructing λ_C from Rₚ=0.8412fm

Equation (S06_03_Rp_final_form) is equivalent to

\begin{equation} \lambda_C=\frac{\pi}{2}\,R_p. \end{equation}

Substituting (S06_03_Rp_value_fm) into (S06_03_lC_from_Rp) gives

\begin{align} \lambda_C &=\frac{\pi}{2}\times 0.8412\ \mathrm{fm} \notag\\ &= 1.3213538700998668\ \mathrm{fm}, \end{align}

which matches (S06_03_lC_value_lock). Therefore the locks Rₚ/L_q=α and L_q=λ_C are numerically consistent within the same version (this check is used only as a Gate input and not as a justification).

6.3.5 Definition and numerics of the geometric cross section σ_geom

[D-6.3-1] Definition of geometric cross section

Once the core radius Rₚ is defined, define the geometric cross section (geometric area) as

\begin{equation} \sigma_{\mathrm{geom}} := \pi R_p^2. \end{equation}

Definition (S06_03_sigma_def) is purely geometric and must not be reinterpreted as a different notion (effective cross section, etc.). If an effective cross section is needed, it must be introduced as a separate symbol with a separate definition.

Numerical substitution

Substitute (S06_03_Rp_value_fm) into (S06_03_sigma_def):

\begin{align} \sigma_{\mathrm{geom}} &=\pi\,(0.8412\ \mathrm{fm})^2 \notag\\ &=\pi\times 0.70761744\ \mathrm{fm}^2 \notag\\ &=2.223045751056016\ \mathrm{fm}^2. \end{align}

With the unit-conversion convention (1fm)²=10⁻³⁰m²,

\begin{equation} \sigma_{\mathrm{geom}} = 2.223045751056016\times 10^{-30}\ \mathrm{m}^2. \end{equation}

6.3.6 Summary of the 4/π invariant (definition–expansion–conclusion)

In this section, an “invariant” means a dimensionless or normalized combination that remains valid inside the regime regardless of unit choices as long as the lock on the ratio Rₚ/L_q is maintained.

6.3.6.1 Invariant I: σ_geom/L_q²=4/π

From the conclusion (S06_03_Rp_over_Lq) of §6.2,

\begin{equation} R_p=\alpha L_q \end{equation}

and since α=2/π,

\begin{equation} R_p=\frac{2}{\pi}L_q. \end{equation}

Substitute (S06_03_Rp_2overpi_Lq) into (S06_03_sigma_def):

\begin{align} \sigma_{\mathrm{geom}} &=\pi R_p^2 =\pi\left(\frac{2}{\pi}L_q\right)^2 \notag\\ &=\pi\left(\frac{4}{\pi^2}\right)L_q^2 \notag\\ &=\frac{4}{\pi}\,L_q^2. \end{align}

Therefore the following invariant holds.

\begin{equation} \boxed{ \frac{\sigma_{\mathrm{geom}}}{L_q^2} = \frac{4}{\pi} } \qquad (\text{Invariant I}). \end{equation}

Applying the lock L_q=λ_C from §6.1 gives the equivalent form

\begin{equation} \boxed{ \frac{\sigma_{\mathrm{geom}}}{\lambda_C^2} = \frac{4}{\pi} } \qquad (\text{Invariant I, equivalent form}). \end{equation}

6.3.6.2 Numerical check (Invariant I)

From (S06_03_lC_value_lock), since L_q=λ_C, we have

\begin{equation} L_q^2 = (1.3213538700998668\ \mathrm{fm})^2 = 1.7459760500278958\ \mathrm{fm}^2. \end{equation}

Using (S06_03_sigma_value_fm2) and (S06_03_Lq2_value),

\begin{align} \frac{\sigma_{\mathrm{geom}}}{L_q^2} &= \frac{2.223045751056016}{1.7459760500278958} \notag\\ &= 1.2732395447351628, \end{align}

and since

\begin{equation} \frac{4}{\pi}=1.2732395447351628, \end{equation}

the numerical ratio agrees with (S06_03_invariant_4overpi). This check is a validation computation; it does not use external texts as a basis for the invariant.

6.3.7 Where the invariant is used (locking insertion points)

Invariant I can be used as a geometric normalization in derivations of the following form, and its use must have its insertion point locked.

  1. In every expression where σ_geom appears, when normalizing σ_geom by L_q², use (S06_03_invariant_4overpi) to fix the constant as 4/π.
  2. For expressions written in terms of λ_C, use the equivalent form (S06_03_invariant_4overpi_lC).

This use is permitted only as a geometric consequence of the already-locked ratio Rₚ/L_q=α and α=2/π, not as an “ad hoc constant correction.”

LOCK/Gate connections for this section (none if empty)

6.4 Continuum→Discrete (82+7) stability conditions

6.4.1 Purpose

This section declares a list of minimal stability conditions required when the core radius Rₚ and related invariants derived in Chapter 6 (e.g., Rₚ/L_q=2/π, σ_geom/L_q²=4/π) descend to the discrete structure “core 82 + shell 7” in Chapter 8. This section does not derive the detailed coordinates or coupling conventions of the discrete structure. It fixes only the necessary conditions imposed by the continuum results on the discrete structure, in the form of (i) condition identifier, (ii) condition content, and (iii) handling on violation.

6.4.2 Premises (locked continuum-side results)

The continuum-side locked results referenced in this section are as follows.

  1. Core radius ratio:
    \begin{equation} \frac{R_p}{L_q}=\frac{2}{\pi}. \end{equation}
  2. Length identification:
    \begin{equation} L_q=\lambda_C. \end{equation}
  3. Geometric cross-section invariant:
    \begin{equation} \frac{\sigma_{\mathrm{geom}}}{L_q^2}=\frac{4}{\pi}, \qquad \sigma_{\mathrm{geom}}:=\pi R_p^2. \end{equation}

These results are chapter outputs of the continuum model and belong to canon_lock. If the discrete structure fails to satisfy them, the continuum→discrete link does not hold.

6.4.3 Minimal stability conditions for the discrete structure (82+7)

The discrete structure consists of the combination “core 82” and “shell 7.” The continuum results require the following conditions as mandatory. Each condition is independent; violating any one of them judges the continuum→discrete link as FAIL.

[C-82/7-01] Core-boundary radius consistency (radius lock condition)

Given a discrete core (82) coordinate set x_i_(i=1)⁸² and a center x_c, the radial aggregate of the key boundary layer (boundary-candidate set B₈₂⊆1,…,82) must be consistent with Rₚ. Declare the consistency in the following form:

\begin{equation} R_{82} := \mathrm{Agg}\Bigl(\{\|\mathbf{x}_i-\mathbf{x}_c\|\}_{i\in\mathcal{B}_{82}}\Bigr) \equiv R_p, \end{equation}

where the aggregation operator Agg (e.g., median, mean, center of a min–max band) must be locked in analysis_lock. If Agg is not locked, R₈₂ is non-unique; the condition is undefined and judged INCONCLUSIVE. The tolerance (allowable error) ε_R must be registered in advance in gate_lock, and

\begin{equation} |R_{82}-R_p|>\varepsilon_{R} \quad\Longrightarrow\quad \texttt{FAIL-CORE82-RADIUS}. \end{equation}

[C-82/7-02] Cross-section invariant consistency (geometric cross-section condition)

In the discrete core (82), it must be possible to define a “discrete cross section” through the boundary-candidate set B₈₂, and its normalization must be consistent with the invariant (S06_04_invariant_sigma). Define the discrete cross section σ₈₂ in the following form:

\begin{equation} \sigma_{82} := \mathrm{ProjArea}\Bigl(\{\mathbf{x}_i\}_{i\in\mathcal{B}_{82}};\ \mathbf{n}_\sigma\Bigr), \end{equation}

where n_σ is the projection axis and is locked in analysis_lock. The definition of the projection-area operator ProjArea (convex hull / lattice-cell count / pixelization, etc.) must also be locked in analysis_lock.

Define the normalized cross-section ratio as

\begin{equation} I_{\sigma,82}:=\frac{\sigma_{82}}{L_q^2}. \end{equation}

The requirement from the continuum invariant is declared as

\begin{equation} \left|I_{\sigma,82}-\frac{4}{\pi}\right|>\varepsilon_{\sigma} \quad\Longrightarrow\quad \texttt{FAIL-CORE82-SIGMA}. \end{equation}

where ε_(σ) is a threshold locked in gate_lock.

[C-82/7-03] Uniqueness of boundary transition (single core-boundary condition)

The continuum model locks the core boundary as a single transition point. The discrete structure must also have a single transition. Define a “core indicator” χ₈₂(r) for radius r and require that there is only one transition point:

\begin{equation} \chi_{82}(r)\in\{0,1\}, \qquad r\mapsto \chi_{82}(r)\ \text{has exactly one }0\to 1\text{ transition}. \end{equation}

The decision rule for the transition (which structural quantity defines χ₈₂, and the thresholds that judge “one transition”) must be locked in analysis_lock. If two or more transitions occur, or if there is no transition, treat it as a failure of the continuum→discrete link:

\begin{equation} \text{transition count}\neq 1 \quad\Longrightarrow\quad \texttt{FAIL-CORE82-TRANSITION}. \end{equation}

[C-82/7-04] Locality of shell (7) attachment (preservation of the local rule)

Shell (7) must form attachment/cancellation/survival structures locally near the core boundary. This is the requirement that the discrete construction is compatible with the local-rule axiom in §3.1. Declare the following condition:

\begin{equation} \forall k\in\{1,\ldots,7\},\ \exists i(k)\in\mathcal{B}_{82}\ \text{s.t.}\ \|\mathbf{s}_k-\mathbf{x}_{i(k)}\|\le \rho_{\mathrm{attach}}. \end{equation}

Here sₖ is a shell vector or shell marker point (locked in Chapter 8), and ρ_attach is the attachment-radius threshold. Lock ρ_attach by normalizing it with the length scale L_q:

\begin{equation} \rho_{\mathrm{attach}} := \eta_{\mathrm{attach}}\,L_q, \qquad \eta_{\mathrm{attach}}>0\ \text{is locked}. \end{equation}

Handle violation as

\begin{equation} \exists k\ \text{s.t.}\ \min_{i\in\mathcal{B}_{82}}\|\mathbf{s}_k-\mathbf{x}_i\|>\rho_{\mathrm{attach}} \quad\Longrightarrow\quad \texttt{FAIL-SHELL7-LOCAL}. \end{equation}

[C-82/7-05] Non-degeneracy of the cancellation–survival convention (existence of a survival vector)

When the continuum results descend to the discrete structure, the cancellation–survival convention of shell (7) must remain non-degenerate. This means that a survival vector V is definable (not non-unique/ambiguous) and must not collapse to the zero vector.

\begin{equation} \mathbf{V}:=\sum_{k=1}^{7}\mathbf{s}_k, \qquad \|\mathbf{V}\|\ge V_{\min}. \end{equation}

V_(min) is a threshold locked in gate_lock in internal units or normalized by L_q. If |V|

\begin{equation} \|\mathbf{V}\|<V_{\min} \quad\Longrightarrow\quad \texttt{FAIL-SHELL7-DEGEN}. \end{equation}

[C-82/7-06] Regime consistency (Stone/plenum/jamming prerequisites)

The continuum core model assumes the Stone/plenum regime. The discrete core (82+7) must also be defined under the same regime. Declare the following conditions.

  1. The discrete placement does not violate impenetrability.
  2. The discrete placement does not conflict with the plenum convention (do not treat void as an independent object).
  3. The contact graph and backbone decision are consistent with the locked regime coordinate axes (§4.3).

If any of the above is violated, treat it as a regime mismatch:

\begin{equation} \text{regime mismatch} \quad\Longrightarrow\quad \texttt{FAIL-REG-MISMATCH}. \end{equation}

[C-82/7-07] Scale-normalization consistency (single length unit)

All coordinates and length judgments of the discrete structure must be normalized by the same length scale L_q (or the identified λ_C). The following mixes are forbidden within the same output.

  1. Using L_q and λ_C as different values.
  2. Re-introducing Rₚ via a separate definition to evade the R₈₂ matching.
  3. Mixing the cell geometry (cube) and a visualization sphere when computing cross sections or radii.

If a mix is detected, it is an immediate FAIL as a lock_id mix or a meaning conflict:

\begin{equation} \text{normalization mix} \quad\Longrightarrow\quad \texttt{FAIL-LOCK-MIX}\ \text{or}\ \texttt{FAIL-GEO-CONF}. \end{equation}

[C-82/7-08] Existence of tetrahedral locking (a 4-point rigidity certificate)

This condition locks the minimal structural certificate needed to state that “core (82) has full jamming rigidity” (i.e., that one can logically claim the c² stiffness corresponding to Ψ_(rm yield) in the unjamming trigger (unjamming_trigger)).

In 3D, four points are the minimal non-coplanar simplex. If the contact network does not contain this tetrahedral (simplex) backbone, shear rigidity can remain zero or regime-out. Therefore the contact graph of core 82, G_c=(V_c,E_c), must satisfy

\begin{equation} \exists\ (i_1,i_2,i_3,i_4)\subset \mathcal{B}_{82} \ \text{s.t.}\ (i_a,i_b)\in \mathcal{E}_c\ \forall\ a<b. \end{equation}

That is, there must exist a 4-point backbone corresponding to the complete graph K₄ within B₈₂. Denote the corresponding 4-point set by T₄.

To claim “rigidity in full,” this tetrahedron must also be non-degenerate (not coplanar or near-coplanar). In coordinate form, judge this by a threshold on tetrahedral volume (mixed product):

\begin{equation} V_{\mathrm{tet}}(i_1,i_2,i_3,i_4) :=\frac{1}{6}\left|\det\left[\mathbf{x}_{i_2}-\mathbf{x}_{i_1},\ \mathbf{x}_{i_3}-\mathbf{x}_{i_1},\ \mathbf{x}_{i_4}-\mathbf{x}_{i_1}\right]\right| \ge V_{\min}^{\mathrm{tet}}. \end{equation}

Here V_(min)^tet is a threshold locked in gate_lock after normalization by L_q³.

If the condition holds, the fact that core (82) has a jamming backbone containing “tetrahedral locking” is secured, and it becomes the minimal basis to claim the jamming regime (ξ≈ 1) in the dynamic rigidity-mixing narrative (ξ∈[0,1]) of this document. Conversely, if this condition is violated, sentences such as “core rigidity = c²” or “no further contraction” lose conclusion status; only limit claims (CT-LIM) are permitted.

\begin{equation} \neg\exists\ \mathcal{T}_4\ \ \text{or}\ \ V_{\mathrm{tet}}<V_{\min}^{\mathrm{tet}} \quad\Longrightarrow\quad \texttt{FAIL-CORE82-TETRA}. \end{equation}

6.4.4 Handling violations (conclusion status and limit statements)

If any of the conditions [C-82/7-01]~[C-82/7-08] is violated, the continuum→discrete link loses conclusion status. In that case, only “limit statements (CT-LIM)” are permitted, and one must record the FAIL label together with the causal condition ID. Patching violations by “interpretation,” or changing condition definitions (aggregation operator, thresholds, boundary-candidate set, etc.) after seeing results, violates No-Tuning and is forbidden. If changes are necessary, the only allowed path is a version bump followed by full re-validation.

LOCK/Gate connections for this section (none if empty)

  \item LOCK: Fix the continuum outputs ($R_p/L_q=2/\pi$, $L_q=\lambda_C$, $\sigma_{\mathrm{geom}}/L_q^2=4/\pi$) as upstream inputs.
  \item LOCK: Lock the definitions of the discrete stability-condition list [C-82/7-01]\textasciitilde[C-82/7-08] (aggregation operators, projection axis, thresholds, normalization conventions) in \texttt{analysis\_lock}/\texttt{gate\_lock}.
  \item Gate: On any condition violation, immediately assign \texttt{FAIL-CORE82-*}/\texttt{FAIL-SHELL7-*}/\texttt{FAIL-REG-*} labels and revoke conclusion status.
  \item Gate: If condition definitions are not locked (aggregation/threshold/boundary-candidate missing), judge \texttt{INCONCLUSIVE}; if conditions are modified post hoc based on outcomes, judge \texttt{FAIL} in G-NT.
  \item Gate: Mixing different \texttt{lock\_id} combinations or confusing cell geometry is immediate \texttt{FAIL} (G-LOCK/G-SYM).