// =============================================================================
// Multi-Center Dynamics of Adjacent Causal Riesz Sectors
// =============================================================================
//
// All algebraic identities machine-checked with Z3 via the Kleis language.
// Complete source: https://kleis.io
//
// =============================================================================

import "stdlib/prelude.kleis"
import "stdlib/templates/arxiv_paper.kleis"

// =============================================================================
// Paper Metadata
// =============================================================================

define paper_title = "Multi-Center Dynamics of Adjacent Causal Riesz Sectors"

define paper_authors = [
    Author("Engin Atik", "1")
]

define paper_affiliations = [
    Affiliation(1, "Kleis Research", "https://kleis.io")
]

define paper_abstract = "In Projected Ontology Theory (POT), the gravitational kernel $K(r) = -A\/r + B ln(r\/r_0)$ is a linear combination of adjacent causal Riesz sectors $R_2^+$ (Newtonian) and $R_3^+$ (galactic). Previous work derived flat rotation curves and the baryonic Tully-Fisher relation (BTFR) for individual galaxies. This paper develops the multi-center dynamics.

The naive coupling $U = M_1 M_2 K(r)$ is inconsistent with galaxy-dependent $lambda_i$: reciprocity in the $R_3^+$ sector requires $lambda\/M = \"const\"$, but the BTFR gives $lambda prop sqrt(M)$. We identify the required structure: a symmetric bilinear $R_3^+$ interaction charge $q_i$, with pair potential
$
U_(i j)(r) = -G M_i M_j \/ r + gamma q_i q_j ln(r\/r_0), quad C_(i j) = gamma q_i q_j = C_(j i).
$
Action--reaction symmetry is then structural. The single-source field becomes a derived quantity, and the relationship between the interaction charge $q_i$ and observables ($lambda_i$, $M_i$) becomes a result rather than an assumption. The BTFR constrains the collective galactic charge $Q prop sqrt(M)$, implying a mass-dependent $R_3^+$ response for coherent extended sources, while microscopic test-body universality is preserved.

For the reduced two-body problem with effective potential $U_(\"eff\")(r) = -a\/r + C ln(r\/r_0) + L^2\/(2 mu r^2)$, the circular-orbit condition is quadratic ($C r^2 + a r - L^2\/mu = 0$), yielding exactly one positive circular radius. Every such orbit is necessarily stable: $U_(\"eff\")''(r_(\"circ\")) = a\/r^3 + 2C\/r^2 > 0$. The strongest result: for $C = gamma q_i q_j > 0$ at arbitrarily large separation, $U_(\"eff\")(r) arrow +infinity$ as $r arrow infinity$, so no finite-energy trajectory can escape to infinity ($L >= 0$, covering radial and nonradial motion).

This identifies a large-distance domain-of-validity constraint: if asymptotically separating states occur, the isolated attractive $C ln(r\/r_0)$ interaction cannot extend unchanged to arbitrarily large separation. Possible completions include finite $R_3^+$ coherence, a sector transition, or environmental modification at large distances --- consistent with the observed convergence steepening beyond the virial radius.

The Bullet Cluster (1E 0657-56) is treated as an application: the GR-inferred convergence slopes ($kappa prop 1\/R$) bracket the observed values ($n = 0.9$--$1.2$), the mass-gas offset follows kinematically, and all descriptions use one parameter per component ($B = sigma^2$). The algebraic identities underlying the construction were machine-checked in Kleis with Z3."

define paper_keywords = "Riesz distributions, multi-center dynamics, interaction charge, galaxy clusters, Bullet Cluster, sector transition, baryonic Tully-Fisher relation, Projected Ontology Theory"


// =============================================================================
// Section 1: Introduction
// =============================================================================

define sec_intro = ArxivSection("Introduction",
"In Paper 1 (Atik 2026a), the gravitational kernel
$
K(r) = -A / r + B ln(r\/r_0)
$
was derived within Projected Ontology Theory (POT) as a linear combination of adjacent causal Riesz sectors $R_2^+$ (Newtonian) and $R_3^+$ (galactic). Flat rotation curves ($v^2 = lambda$) and the baryonic Tully-Fisher relation ($M_(\"bar\") = kappa v^4$) were established as Z3-verified theorems. Paper 1 identified multi-galaxy systems as an open test arena.

The present paper closes that gap by developing the *multi-center dynamics* of the $R_3^+$ sector. The central structural fact is *linearity*: admissible POT kernels are linear operators on flows (Atik 2026b, `kernel_lin_add`). When two galaxies interact, their potentials add. But two questions arise immediately:

*First, the coupling question.* If each galaxy generates a field $Phi_j = -mu_j\/r + lambda_j ln(r\/r_0)$, and the force on center $i$ is $bold(F)_(i j) = -M_i nabla Phi_j$, then the $R_3^+$ force components are $F_(i j)^((\"log\")) = M_i lambda_j \/ r$ and $F_(j i)^((\"log\")) = M_j lambda_i \/ r$ (forces depend only on $d\/d r [ln(r\/r_0)] = 1\/r$). These are equal and opposite only if $lambda_i \/ M_i = lambda_j \/ M_j$. But the BTFR gives $lambda prop sqrt(M)$, so $lambda \/ M prop 1\/sqrt(M)$, which varies. The naive pairwise potential $U = M_1 M_2 K(r)$ is *inconsistent* with observed galaxy-dependent $lambda_i$.

*Second, the escape question.* For an attractive logarithmic potential $U(r) tilde C ln(r\/r_0)$ with $C > 0$, $U(r) arrow +infinity$ as $r arrow infinity$. The usual Newtonian bound/unbound classification does not apply. The escape energy is infinite.

Reciprocity and the BTFR together impose requirements on any viable $R_3^+$ interaction law. The minimal factorized bilinear realization is a symmetric $R_3^+$ interaction charge $q_i$ with coupling $C_(i j) = gamma q_i q_j$, making action--reaction symmetry structural rather than patched. The second issue, rather than being a nuisance, provides an independent theoretical constraint: the logarithmic sector must have a finite domain of validity.

The paper proceeds as follows. Section 2 states assumptions. Sections 3--4 develop multi-center superposition and the phantom density profile. Section 5 introduces the $R_3^+$ interaction charge and derives the symmetric pair force law. Section 6 derives conservation laws from interaction symmetry. Section 7 develops the two-body dynamics, proving existence and stability of circular orbits and the impossibility of asymptotic escape. Section 8 derives the finite-domain constraint. Section 9 discusses phantom density and lensing convergence as GR-observer inferences. Section 10 applies the theory to the Bullet Cluster. Sections 11--12 compare with alternatives and discuss limitations.

All algebraic identities underlying the construction were machine-checked in Kleis with Z3. A complete verification table appears in the Appendix.")


// =============================================================================
// Section 2: Assumptions
// =============================================================================

define sec_assumptions = ArxivSection("Assumptions and Applicability",
"Every encoding decision and physical idealization is stated explicitly. If a conclusion depends on an assumption, it is listed here.

=== A1. Spherical symmetry of the Green kernel

The kernel $K(r) = -A\/r + B ln(r\/r_0)$ is the spherically symmetric ($ell = 0$ multipole) Green's function. Extended, aspherical sources produce aspherical fields by convolution, but the kernel itself has no preferred direction. The superposition theorems depend *only* on linearity, not on spherical symmetry of the source.

=== A2. Galaxy clusters vs. individual galaxies

The $R_3^+$ sector produces radial acceleration $a = -B\/r$ regardless of system type. For clusters, $B = sigma^2$ rather than $v_c^2$. The resulting phantom density $rho prop 1\/r^2$ is the isothermal sphere --- derived, not assumed.

=== A3. No morphological restriction

The kernel does not depend on galaxy morphology. The Bullet Cluster comprises 291 spectroscopic members (Cha et~al.~2025), predominantly elliptical galaxies.

=== A4. Bullet Cluster geometry

System 1E 0657-56 at $z = 0.296$, observed $tilde 10 degree$ from the plane of the sky. JWST + DECam mass ratio $q = 10.14^(+3.22)_(-2.47)$ (Cho et~al.~2025). Halo relative velocity $tilde 2600$ km/s. Time since core passage $tilde 100$--$120$ Myr.

=== A5. Static kernel approximation

The superposition uses the $omega = 0$ restriction of the retarded Riesz distributions, valid because the light-crossing time ($tilde 3$ Myr) is much shorter than the dynamical timescale ($tilde 100$ Myr).

=== A6. Three-body structure

JWST reveals the main cluster is bimodal (two halos at ratio $tilde 1.3:1$, separated by 170 kpc). POT handles $N$ centers by linear superposition.

=== A7. No kernel self-interaction

Admissible kernels are linear operators (`kernel_lin_add`). Experimental bounds on dark matter self-interaction $sigma\/m$ are not applicable.

=== A8. Gas dynamics are external

Gas dynamics (ram pressure, shocks, stripping) are external to the kernel algebra. POT uses the same gas physics as $Lambda$CDM and MOND.")


// =============================================================================
// Section 3: Multi-Center Field Superposition
// =============================================================================

define sec_superposition = ArxivSection("Multi-Center Field Superposition",
"=== Two-center potential (GI-1)

By kernel linearity:
$
Phi_(\"total\")(bold(x)) = [-mu_A / r_A + lambda_A ln(r_A \/ r_0)] + [-mu_B / r_B + lambda_B ln(r_B \/ r_0)]
$
where $r_i = |bold(x) - bold(x)_i|$ and $r_0$ is an arbitrary reference scale (cancels from forces). When centers coincide ($r_A = r_B = r$), $lambda_A ln(r\/r_0) + lambda_B ln(r\/r_0) = (lambda_A + lambda_B) ln(r\/r_0)$ (GI-2). The _instantaneous_ test-particle circular velocity is $v_(\"coinc\")^2 = lambda_A + lambda_B$ (GI-3).

This is a coincident-source superposition, not necessarily the equilibrium source strength of a recohered merged object. If the coherence functional $Q = cal(C)[rho, \"structure\"]$ is nonlinear, the relaxed post-merger source may have
$
lambda_(\"eq\") = sqrt(lambda_A^2 + lambda_B^2) quad (\"BTFR equilibrium with additive mass\"),
$
which for identical progenitors gives $lambda_(\"eq\") = sqrt(2) dot lambda < 2 lambda = lambda_(\"coinc\")$.

The relaxation sequence is therefore:
$
lambda_(\"coinc\") = lambda_A + lambda_B quad arrow.r quad lambda_(\"eq\") = sqrt(lambda_A^2 + lambda_B^2).
$
The overshoot $eta = lambda_(\"coinc\") \/ lambda_(\"eq\")$ decays as the merged system recoheres through $cal(C)$.

=== Unequal-mass generalization (GI-8)

For logarithmic strength ratio $q = lambda_A \/ lambda_B$ (not mass ratio):
- Instantaneous velocity amplification: $v_(\"coinc\")^2 \/ v_B^2 = q + 1$.
- BTFR overshoot factor at coincidence: $f = (q+1)^2 \/ (q^2+1)$, bounded in $(1, 2]$.
- For identical mergers ($q = 1$): $f = 2$, a $100%$ excess in $v^4$.

Note: $q$ is the *logarithmic strength* ratio, not the total mass ratio. Under BTFR, $M prop lambda^2$, so a mass ratio of 10 corresponds to $q = sqrt(10) approx 3.16$. For clusters, whether the galaxy BTFR applies directly is an open question (Section 12).

#rect(stroke: 0.5pt, inset: 8pt, width: 100%)[*Three dynamical regimes.* (1) The *approach trajectory* of two coherent centers is governed by the pair interaction $U_(i j)$ (Section 5), involving receiver-dependent coupling $C_(i j)$. (2) At *coincidence*, the instantaneous superposed field gives $v_(\"coinc\")^2 = lambda_A + lambda_B$ and the overshoot factor $f$ (GI-2, GI-3, GI-8). (3) After *recoherence*, the relaxed object's field is determined by $cal(C)[rho_(\"merged\")]$ and must satisfy the BTFR with the total mass, yielding $v_(\"eq\")^2 = sqrt(lambda_A^2 + lambda_B^2)$. The overshoot $f$ could in principle become observable as a transient elevated rotation velocity in recently merged galaxies.]")


// =============================================================================
// Section 4: Phantom Density
// =============================================================================

define sec_phantom = ArxivSection("Phantom Density Profile",
"The $R_3^+$ acceleration $a_(\"log\") = -B\/r$ is interpreted by a Newtonian observer as a phantom mass:
$
M_(\"phantom\")(r) = B r \/ G, quad rho_(\"phantom\")(r) = B / (4 pi G r^2).
$
This is the isothermal sphere --- universal $1\/r^2$, one parameter ($B = sigma^2$). It replaces the NFW profile's two fitted parameters ($M_(200)$, concentration $c$). Z3-verified properties:
- $M_(\"phantom\")(2r) = 2 M_(\"phantom\")(r)$ (linear growth, GI-9b).
- $rho(2r) = rho(r)\/4$ ($1\/r^2$ law, GI-9c).
- $4 pi r^2 rho = B\/G$ (constant mass per shell, GI-9d).")


// =============================================================================
// Section 5: The R₃⁺ Interaction Charge
// =============================================================================

define sec_charge = ArxivSection("The $R_3^+$ Interaction Charge",
"=== The coupling problem

If $Phi_j = -mu_j\/r + lambda_j ln(r\/r_0)$ and $bold(F)_(i j) = -M_i nabla Phi_j$, then
$
F_(i j)^((\"log\")) = M_i lambda_j / r, quad F_(j i)^((\"log\")) = M_j lambda_i / r.
$
Newton's third law requires $M_i lambda_j = M_j lambda_i$, i.e., $lambda\/M = \"const\"$. But the BTFR gives $lambda prop sqrt(M)$, so $lambda\/M prop 1\/sqrt(M)$ --- not constant. The naive $U = M_1 M_2 K(r)$ is inconsistent.

=== The required structure: symmetric bilinear coupling

Reciprocity rules out directly identifying the phenomenological single-source coefficient $lambda_i$ with an ordinary mass-proportional pair coupling. A symmetric pair interaction therefore requires an additional coupling structure. The minimal *factorized bilinear* realization of the required symmetric coupling is a sector charge $q_i$ with $C_(i j) = gamma q_i q_j$. Other symmetric functions $C_(i j) = F(Q_i, Q_j)$ with $F(x,y) = F(y,x)$ would also satisfy reciprocity; we study the bilinear case as the simplest. We distinguish three objects:
- $M_i$: Newtonian/gravitational charge ($R_2^+$ coupling).
- $q_i$: $R_3^+$ interaction charge (new).
- $lambda_i$: observed single-source logarithmic response (phenomenological).

The pairwise interaction potential is
$
U_(i j)(r) = -G M_i M_j / r + gamma q_i q_j ln(r \/ r_0),
$
where $C_(i j) = gamma q_i q_j = C_(j i)$ is structurally symmetric and $r_0$ is an arbitrary reference scale. The reference scale adds a constant $C_(i j) ln r_0$ to the potential; it cancels from all forces ($d\/d r ln(r\/r_0) = 1\/r$), from the circular-orbit equation, and from the no-escape theorem (which depends only on $C ln r arrow +infinity$). Action--reaction symmetry is automatic (GI-12a, GI-13a). The force decomposes as a Newtonian $1\/r^2$ term plus a logarithmic $1\/r$ term (GI-12c). The pair-dependent crossover radius is
$
r_(c, i j) = G M_i M_j / (gamma q_i q_j)
$
(GI-12d, GI-12f).

=== Universality structure

At the microscopic level, universal free fall in $R_3^+$ requires $q\/M = \"const\"$ for test bodies. But galactic source strengths show $lambda prop sqrt(M)$ (BTFR), which constrains coherent extended sources to acquire a collective effective charge $Q_(\"galaxy\") = cal(C)[rho, \"structure\", \"coherence\"]$ distinct from the elementary probe coupling:
$
q_(\"micro\") / M = \"const\" quad (\"universality\"), quad Q_(\"galaxy\") prop sqrt(M) quad (\"BTFR constraint\").
$

The consequence is that coherent galaxies do *not* experience universal free fall in $R_3^+$:
$
a_i^((3)) = gamma Q_i Q_j / (M_i r) prop sqrt(M_j) / (sqrt(M_i) dot r),
$
which depends on $M_i$. Microscopic equivalence and collective coherent response are *different levels of description*. This is testable physics: the $R_3^+$ sector distinguishes test particles from extended coherent sources.

=== Relation to the source field (bridge)

The single-source field $Phi_j = -mu_j\/r + lambda_j ln(r\/r_0)$ (Sections 3--4) and the pair interaction $U_(i j) = -G M_i M_j\/r + gamma q_i q_j ln(r\/r_0)$ (this section) are different effective descriptions:

- *Test-particle field.* A microscopic test body ($q_i = alpha M_i$, mass $M_i arrow 0$) in the field of source $j$ experiences acceleration $a_i = gamma alpha q_j \/ r$. Identifying $lambda_j = gamma alpha q_j$ recovers the source-only field $Phi_j$. This is the regime of Sections 3--4.

- *Coherent-center interaction.* Two galaxies with effective charges $Q_i$, $Q_j$ interact through $C_(i j) = gamma Q_i Q_j$. The acceleration of galaxy $i$ is $a_i = gamma Q_i Q_j \/ (M_i r)$, which depends on the receiver. This is the regime of Sections 5--7.

The source field $lambda_j$ is therefore a derived quantity: it is the acceleration per unit $q\/M$ induced by source $j$, valid for probes with universal $q\/M = alpha$. When both objects are coherent extended sources, the pair coupling $C_(i j)$ replaces the single-source parametrization.

=== Differential acceleration of coherent sources

The non-universal free fall prediction has a definite mass scaling. If the BTFR constraint gives $Q_i = k sqrt(M_i)$ for coherent extended sources, then the $R_3^+$ acceleration of source $i$ in the field of source $j$ is
$
a_(i arrow.l j)^((3)) = (gamma Q_i Q_j) / (M_i r) = (gamma k^2 sqrt(M_j)) / (sqrt(M_i) dot r).
$
For two coherent bodies $A$, $B$ at the same distance $r$ from a common source $j$:
$
|a_A^((3))| / |a_B^((3))| = sqrt(M_B) / sqrt(M_A).
$
The heavier coherent body accelerates *less* in the $R_3^+$ field. By contrast, a microscopic test body ($q = alpha M$) experiences acceleration $lambda_j \/ r$ independent of its mass.

This is not merely qualitative. The mass scaling $a prop 1\/sqrt(M)$ for coherent receivers is a definite functional form, testable in principle via satellite galaxy dynamics or cluster infall: a dwarf galaxy ($M tilde 10^9 M_sun$) falling toward a giant ($M tilde 10^(12) M_sun$) would experience $R_3^+$ acceleration $tilde 30$ times larger than the giant experiences from the dwarf, whereas the test-particle (Newtonian equivalence) prediction gives equal $R_3^+$ acceleration per unit mass. The absolute normalization of the effect remains tied to the unresolved coherence functional $cal(C)$, but the *scaling* is determined.")


// =============================================================================
// Section 6: Conservation Laws
// =============================================================================

define sec_conservation = ArxivSection("Conservation Laws from Interaction Symmetry",
"Because $C_(i j) = gamma q_i q_j = C_(j i)$, action--reaction symmetry follows directly from the pair potential:
$
bold(F)_(i j) = -bold(F)_(j i) quad (\"GI-13a\").
$
Thus the analogue of Newton's third law is a derived property of the interaction, not an independent postulate. Since $U_(i j)$ depends only on $|bold(x)_i - bold(x)_j|$, the force between any pair is central (along the separation vector). The standard conservation chain follows:
- *Momentum*: $d bold(P)_(\"total\")\/d t = 0$ from translational invariance of $U(|bold(x)_i - bold(x)_j|)$. Internal forces cancel pairwise by the derived action--reaction symmetry (GI-13a).
- *Angular momentum*: $d bold(L)_(\"total\")\/d t = 0$ from rotational invariance. Since the force is central ($bold(r)_(i j) times bold(F)_(i j) = 0$), no torque is generated.
- *Center of mass*: $bold(v)_(\"COM\") = \"const\"$ (GI-13b). For unequal masses, the lighter object accelerates more (GI-13c), but the center of mass remains unaccelerated.

For $N$ centers, there are $N(N-1)\/2$ distinct pairs (GI-13d), each with a symmetric interaction. These conservation laws follow from the symmetry and invariances of the pair potential; they are not imported as Newtonian axioms.")


// =============================================================================
// Section 7: Two-Body Dynamics
// =============================================================================

define sec_twobody = ArxivSection("Two-Body Dynamics with Logarithmic Potential",
"=== Effective potential

The reduced two-body problem uses $mu_(\"red\") = M_A M_B \/ (M_A + M_B)$, relative coordinate $r = |bold(x)_A - bold(x)_B|$, and effective potential
$
U_(\"eff\")(r) = -a / r + C ln r + L^2 / (2 mu r^2), quad a = G M_A M_B, quad C = gamma q_A q_B > 0.
$

=== Unique positive circular radius (GI-14b)

Setting $d U_(\"eff\") \/ d r = 0$:
$
a / r^2 + C / r = L^2 / (mu r^3).
$
Both the Newtonian and logarithmic forces are attractive (same side), balanced by the centrifugal term. Multiplying by $r^3$:
$
C r^2 + a r - L^2 / mu = 0.
$
This is *quadratic*, not cubic. The discriminant $Delta = a^2 + 4 C L^2 \/ mu > 0$ always (GI-14b2), giving exactly one positive root:
$
r_(\"circ\") = (-a + sqrt(a^2 + 4 C L^2 \/ mu)) / (2 C).
$

=== Necessary stability (GI-14c)

The second derivative at the circular orbit is
$
U_(\"eff\")''(r) = -2a / r^3 - C / r^2 + 3 L^2 / (mu r^4).
$
Using $L^2\/mu = a r + C r^2$ from the orbit equation:
$
U_(\"eff\")''(r_(\"circ\")) = a / r_(\"circ\")^3 + 2 C / r_(\"circ\")^2 > 0.
$
Since $a > 0$, $C > 0$, and $r > 0$, *every circular orbit of the $R_2^+ + R_3^+$ pair potential is stable*. This is a clean theorem, not a case analysis.

=== No finite-energy asymptotic escape (GI-14d)

The machine-checked statement is:
$
forall E in RR, exists R > 0, forall r > R: quad U_(\"eff\")(r) > E.
$
The proof decomposes into four verified lemmas: (L1) for $r > a$, $-a\/r > -1$; (L2) the centrifugal barrier $L^2\/(2 mu r^2) >= 0$ (valid for all $L$ including $L = 0$, covering radial motion); (L3) for any target, $exists R$ with $C ln R > E + 1$; (L4) $ln$ is strictly increasing on $(0, infinity)$, axiomatized as a property of the real logarithm. The assembly is then machine-verified *relative to this axiom*. The assembly combines all three terms of $U_(\"eff\")$ and quantifies over *all* $r > R$, not merely one witness:
$
U_(\"eff\")(r) = underbrace(-a\/r, > -1) + underbrace(C ln r, > E + 1) + underbrace(L^2\/(2 mu r^2), >= 0) > E.
$
Therefore *finite-energy asymptotic escape $r arrow infinity$ is impossible*.

For $L = 0$, the Newtonian $-a\/r$ term drives $U arrow -infinity$ as $r arrow 0$, so radial collapse remains possible in the point-source idealization. The theorem concerns escape, not capture.

=== Thermalization and merger (GI-14e)

A two-point-mass Hamiltonian cannot irreversibly dissipate relative orbital energy (Liouville theorem). But a galaxy has $tilde 10^(11)$ internal degrees of freedom. Even in a fully Hamiltonian system:
$
E_(\"relative\") arrow E_(\"internal\") + E_(\"tidal\") + E_(\"ejecta\").
$
Individual stars are not the two reduced coordinates --- they are third bodies in a time-dependent potential. Some gain energy and are ejected; others become more tightly bound. Stars carry away positive energy and angular momentum while the remaining cores bind further. No fundamental dissipation need be inserted.

If the $R_3^+$ field extends globally, even ejected stars face the same infinite escape-energy barrier: tidal ejection is local (to large but finite $r$), not truly asymptotic, unless $R_3^+$ has finite domain.

=== Relation to violent relaxation

The redistribution in GI-14e has a well-established collisionless analogue in galactic dynamics. In *violent relaxation* (Lynden-Bell 1967), a rapidly varying collective gravitational potential changes individual stellar energies even though the full many-body dynamics remains Hamiltonian. Subsequent phase mixing can drive the system toward a quasi-stationary state without requiring microscopic dissipation. Post-merger simulations measure a rapid ($tilde 1$ Gyr) phase in which the energy distributions of the two progenitors converge by $tilde 80%$, followed by a slower ($tilde 8$ Gyr) evolution toward a fully mixed state (Young et~al.~2018). This provides a possible dynamical setting for the post-merger recoherence introduced in Section 3: the transition from the coincident-source strength $lambda_(\"coinc\")$ to the BTFR-constrained equilibrium strength $lambda_(\"eq\")$ may occur during collective phase-space relaxation. Importantly, changing the gravitational law changes the relaxation and merging behavior: MOND $N$-body simulations show significantly longer violent relaxation and merging timescales than Newtonian gravity with dark matter (Nipoti et~al.~2007, 2008). We do not derive this identification here. In particular, the functional $cal(C)[rho, \"structure\", \"coherence\"]$ and its relaxation timescale remain unknown. Determining whether collisionless evolution under the $R_2^+ + R_3^+$ interaction dynamically produces the required recoherence is an $N$-body problem.

*Epistemological note:* GI-14a through GI-14d (including the three lemmas) are verified algebraic identities checked by Z3. GI-14e is a physical consequence of Hamiltonian mechanics stated analytically; the Liouville theorem is not encoded in Kleis.")


// =============================================================================
// Section 8: The Finite-Domain Constraint
// =============================================================================

define sec_transition = ArxivSection("The Finite-Domain Constraint on $R_3^+$",
"=== The mathematical theorem (GI-15a)

For $C = gamma q_i q_j > 0$ at arbitrarily large separation, no finite-energy asymptotic two-body scattering state exists. This is established by GI-14d.

=== The physics inference (GI-15b)

We do not literally observe an astronomical pair at $t arrow infinity$, $r arrow infinity$. We observe systems whose kinematics, environments, and reconstructed trajectories are conventionally interpreted as flybys, escaping systems, or hyperbolic encounters.

The logical structure is:
$
C = gamma q_i q_j > 0 \"for arbitrarily large\" r arrow.r.double \"no finite-energy scattering at infinity\" quad [\"theorem\"].
$
If astrophysical systems genuinely possess asymptotically separating states, as standard dynamical reconstructions indicate, then the isolated $C ln(r\/r_0)$ interaction cannot extend unchanged to arbitrarily large separation. The two-body dynamics therefore identifies a previously undetermined domain-of-validity condition, requiring one of:
- Finite $R_3^+$ coherence length.
- Transition to a less confining asymptotic regime --- for example, reversion to $R_2^+$, deactivation of the $R_3^+$ contribution, or participation of higher sectors only with signs\/cancellations that avoid stronger confinement.
- Environmental\/cosmological modification of the isolated pair approximation.
- Some equivalent breakdown of the global logarithmic approximation.

In particular, a naive transition to an attractive higher sector such as $R_4^+ tilde r$ or $R_5^+ tilde r^2$ would not resolve the escape constraint, since those asymptotics are more confining than $ln r$.

=== Scale universality (GI-15c)

The finite-domain issue is not peculiar to galaxy-center scattering. It arises whenever *any* object --- star, galaxy, cluster --- attempts to escape a globally logarithmic well. Even stellar tidal ejection during mergers is only local, not truly asymptotic, unless $R_3^+$ terminates. This makes the finite-domain constraint a universal consequence of the theory, not a special feature of the two-galaxy problem.

=== Supporting evidence (GI-15d)

The observed convergence profile steepens beyond the virial region ($kappa$ departs from $1\/R$) rather than remaining $1\/R$ indefinitely (Umetsu et~al.~2011). This independently supports the finite-domain interpretation.

=== Summary

#quote(block: true)[The adjacent $R_3^+$ sector admits a conservative multi-center dynamics, but consistency with scattering forces breakdown of the globally isolated $C ln r$ approximation.]

The multi-center dynamics thus determines a constraint on the previously unspecified large-distance completion of the single-center $R_3^+$ theory.")


// =============================================================================
// Section 9: Lensing as GR-Observer Inference
// =============================================================================

define sec_lensing = ArxivSection("Phantom Density, GR-Inferred Convergence, and Mass-Gas Offset",
"=== Phantom density and GR-inferred convergence (GI-10)

The $R_3^+$ acceleration produces phantom mass $M_(\"phantom\")(r) = B r\/G$ (Section 4). The quantity $rho_(\"phantom\") prop 1\/r^2$ is the mass distribution a *Newtonian or GR dynamical interpretation* assigns to the POT acceleration field --- it is not a geometric lensing prediction. POT's actual photon deflection depends on the metric representative, which remains an open problem.

The projected surface density is
$
Sigma(R) = B / (4 G R), quad kappa(R) = B / (4 G Sigma_(\"cr\") R).
$
This is what a GR observer *infers* from the effective mass distribution. The $kappa prop 1\/R$ scaling is therefore an interesting correspondence with the GR-inferred effective density, not evidence for POT lensing per se. Whether the profile shape survives at the conformal level requires a separate argument that has not been completed.

Properties of the GR-inferred convergence (verified algebraic identities):
- Two-center convergence superposes linearly (GI-10c).
- Convergence ratio at equal distance = $B_A\/B_B$ (GI-10d).

=== Mass-gas offset (GI-11)

The kernel center tracks the collisionless component via momentum conservation. Gas is retarded by hydrodynamic friction. After gas halts:
$
Delta x = v_0 (t - tau)
$
where $tau$ is the gas friction timescale. This is a *kinematic consequence conditional on the gas-retardation timescale* $tau$ and is not a distinctive POT prediction --- any theory with collisionless gravitational centers produces the same kinematics. The timescale $tau approx 40$ Myr is an external input from gas dynamics, not derived from POT. For the Bullet Cluster: $Delta x approx 186$ kpc, consistent with the observed $approx 150$ kpc subcluster displacement.

=== The source-location problem

A linear kernel acts on its baryonic source distribution. The Bullet Cluster ICM gas constitutes $10$--$16%$ of the baryonic mass within 250 kpc apertures. If $R_3^+$ responds to ordinary mass density, the displaced gas should shift the effective gravitational center toward the gas, not the galaxies.

This is precisely where the coherence charge $Q = cal(C)[rho, \"structure\", \"coherence\"]$ becomes structurally relevant. If the $R_3^+$ source is the *coherent* charge rather than the mass density, then a gravitationally coherent stellar system (galaxy, BCG) carries a different effective $R_3^+$ charge density from shocked, thermalized ICM plasma. The collisionless component dominates the $R_3^+$ field not because it has more mass, but because it maintains the phase-space coherence that sources $Q$.

Rather than a single prediction, this gives the coherence functional $cal(C)$ a third independent boundary condition. The three constraints are:
$
cases(
  Q_(\"gal\") prop sqrt(M) & \"(BTFR)\",
  q_(\"micro\") \/ M = \"const\" & \"(microscopic universality)\",
  Q_(\"shocked gas\") < Q_(\"coherent stellar\") & \"(Bullet morphology)\".
)
$
A future construction of $cal(C)$ that cannot simultaneously satisfy all three kills this version of the theory. The constraints are therefore falsifiable mathematical development, not an unconstrained rescue parameter. Until $cal(C)$ is derived, the Bullet Cluster offset demonstrates consistency with the coherence interpretation and sharpens the requirements on $cal(C)$, but does not independently confirm it.")


// =============================================================================
// Section 10: Application — The Bullet Cluster
// =============================================================================

define sec_bullet = ArxivSection("Application: The Bullet Cluster",
"The following comparisons are *qualitative consistency checks*, not precision fits. Observational uncertainties (especially in convergence slope) are substantial, and the comparisons depend on the GR-inferred convergence framework discussed in Section 9.

=== GR-inferred convergence profile slope

#figure(
  table(columns: 3, stroke: 0.5pt, inset: 6pt,
    table.header([*Component*], [*Observed $n$*], [*GR-inferred POT $n$*]),
    [Main cluster], [$1.2$], [$1$],
    [Subcluster], [$0.9$], [$1$],
    [Mean], [$1.05$], [$1$],
  ),
  caption: [GR-inferred convergence profile slope $kappa prop R^(-n)$. The $R_3^+$ phantom density yields $n = 1$. Model-derived slopes from lensing reconstructions bracket this value.],
  kind: table,
)

=== BTFR overshoot (coincident-source transient)

The coincident-source overshoot factor $f = (q+1)^2\/(q^2+1)$ applies at the instant of center alignment for any pair with $q = lambda_A\/lambda_B$. For identical galaxies ($q = 1$), $f = 2$; for $q = 3$, $f = 8\/5$. This is a *transient* state: the recohered merged object relaxes toward $lambda_(\"eq\") = sqrt(lambda_A^2 + lambda_B^2)$ to restore BTFR consistency, with the overshoot decaying on the recoherence timescale of $cal(C)$.

Whether this formula applies to the Bullet Cluster *as a cluster merger* requires establishing a cluster analogue of the galaxy BTFR. The observed 10:1 mass ratio is $M_A\/M_B$, not $lambda_A\/lambda_B$; under BTFR, these differ by a square root. We therefore present the overshoot as a general algebraic prediction for galaxy mergers rather than a specific Bullet Cluster number.

=== Mass-gas offset

Free-streaming upper bound: $v_0 dot t approx 293$ kpc. With $tau approx 40$ Myr: $Delta x approx 186$ kpc. Observed: $approx 150$ kpc subcluster displacement (Cha et~al.~2025). As noted in Section 9, the offset is a generic kinematic consequence of any collisionless gravitational center; the POT-specific content is that no particulate dark component is needed.

=== GR-inferred convergence ratio

Observed (Clowe et~al.~2006): $overline(kappa)_(\"main\") = 0.36$, $overline(kappa)_(\"sub\") = 0.20$, ratio $= 1.80$. GR-inferred POT single-halo value: $49\/36 approx 1.36$. The discrepancy is qualitatively explained by the main cluster's bimodal structure (two BCGs contribute within the 100 kpc aperture).")


// =============================================================================
// Section 11: Comparison
// =============================================================================

define sec_comparison = ArxivSection("Comparison with Other Approaches",
"#figure(
  table(columns: 4, stroke: 0.5pt, inset: 6pt,
    table.header([], [$Lambda$CDM], [MOND], [POT]),
    [Method], [Simulate + fit], [Nonlinear PDE], [Derive algebraic identities],
    [Single-source params], [2 ($M_(200)$, $c$)], [1--3 ($mu$, $a_0$, $m_nu$)], [1 ($B = sigma^2$)],
    [Profile shape], [NFW (fitted)], [Depends on $mu$], [$1\/r^2$ (derived)],
    [Superposition], [Linear], [Fails], [Linear],
    [$N$-body force], [Simulated], [Expensive PDE], [Closed-form pairwise],
    [Mass-gas offset], [From simulation], [Not predicted], [Kinematic (cond. on $tau$)],
    [Strong lensing], [Full prediction], [Full prediction], [Open],
    [Escape/scattering], [Standard], [Standard], [Finite-domain constraint],
    [External field effect], [Ordinary tidal], [Nonlinear EFE], [Ordinary tidal],
    [Formal consistency], [Numerical simulation], [Numerical PDE], [SMT machine-checked],
  ),
  caption: [Structural comparison of approaches.],
  kind: table,
)

=== External field effect

In MOND, the internal dynamics of a subsystem can be modified by the magnitude of the external field it is embedded in --- the *external field effect* (EFE) --- because MOND is a nonlinear theory. POT's kernel is linear. The total potential is $Phi_(\"total\") = Phi_(\"internal\") + Phi_(\"external\")$, and both contribute acceleration: $bold(a) = -nabla Phi_(\"internal\") - nabla Phi_(\"external\")$. A spatially uniform external acceleration disappears from relative internal motion in a freely falling frame, but an external gradient (tidal field) does not:
$
Delta a_(\"ext\") tilde (nabla nabla Phi_(\"ext\")) dot Delta bold(x).
$
The distinction is that POT has *no nonlinear modification* of the internal field law caused by the magnitude of an external field. Ordinary superposed external and tidal fields remain. This is in principle testable: POT predicts ordinary tidal effects on subclusters, while MOND predicts the nonlinear EFE.")


// =============================================================================
// Section 12: Limitations
// =============================================================================

define sec_limitations = ArxivSection("Limitations and Open Problems",
"=== What POT can derive

Profile shape ($1\/r^2$), GR-inferred convergence scaling ($kappa prop 1\/R$ within the $R_3^+$ domain), mass peak locations, mass-gas offset (conditional on $tau$), multi-center superposition, coincident-source BTFR overshoot (transient), pair-dependent crossover, unique stable circular orbits, and the finite-domain constraint.

=== What POT cannot yet predict

*Strong lensing and actual photon deflection.* Requires the metric representative (open). The GR-inferred convergence is a derived correspondence, not a lensing prediction.

*Convergence at large $R$.* The GR-inferred $1\/R$ profile holds within the $R_3^+$ domain. The finite-domain constraint (GI-15) *predicts* that $kappa$ must steepen at some scale, consistent with observations. The transition scale is not yet derived.

*The coherence functional $cal(C)$.* Reciprocity and the BTFR together constrain $Q prop sqrt(M)$ for coherent extended sources. The detailed functional form of $cal(C)[rho, \"structure\"]$ is not determined --- this is a key open problem for future work. The bilinear sector charge $C_(i j) = gamma q_i q_j$ is the minimal factorized bilinear realization of the required symmetric coupling; other symmetric constructions $C_(i j) = F(Q_i, Q_j)$ are not excluded.

*Cluster BTFR analogue.* The coincident-source overshoot formula $f = (q+1)^2\/(q^2+1)$ is derived for *instantaneous* superposed fields with additive logarithmic strengths. Whether this applies to cluster-scale mergers (where $B = sigma^2$ replaces $v_c^2$) requires establishing a cluster analogue of the galaxy BTFR. Until then, the overshoot is a general algebraic prediction for galaxy mergers, not a specific Bullet Cluster number.

*Post-merger recoherence.* Kernel linearity guarantees $Phi[rho_A + rho_B] = Phi[rho_A] + Phi[rho_B]$, so the coincident-source field is exact. But if the coherence functional $cal(C)$ is nonlinear, $cal(C)[rho_A + rho_B] eq.not cal(C)[rho_A] + cal(C)[rho_B]$, and the relaxed merged source may have $lambda_(\"eq\") eq.not lambda_A + lambda_B$. The BTFR-consistent equilibrium $lambda_(\"eq\") = sqrt(lambda_A^2 + lambda_B^2)$ is a constraint, not a derivation, until $cal(C)$ is constructed. Collisionless violent relaxation (Section 7) provides a candidate physical process for this recoherence, but existing merger simulations use either Newtonian gravity with dark matter or MOND --- not the $R_2^+ + R_3^+$ pair potential. An $N$-body calculation with $U = -a\/r + C ln(r\/r_0)$ is therefore required to determine whether the remnant approaches the BTFR-constrained state and on what timescale.

*Non-universal free fall.* The scaling $a prop 1\/sqrt(M)$ for coherent sources (Section 5.0.3) is a testable consequence of $Q prop sqrt(M)$, but quantitative predictions require knowing the absolute normalization of $gamma k^2$, which depends on the coherence functional $cal(C)$.

*Scale dependence and sector selection.* The finite-domain result raises a broader scale question. The $R_2^+$ contribution reproduces the familiar inverse-square regime where it dominates, while the adjacent $R_3^+$ contribution is relevant to the galactic regime considered here. GI-14d shows that the latter cannot simply be extrapolated unchanged to arbitrarily large separation if asymptotic scattering states are admitted. This motivates, but does not determine, the question of whether the physical sector composition changes again on cluster or cosmological scales. No particular higher-sector assignment is assumed here. The formal hierarchy $R_3^+ tilde ln r$, $R_4^+ tilde r$, $R_5^+ tilde r^2, ...$ does not by itself answer this question, since attractive participation of the higher sectors would increase rather than remove the asymptotic confinement.

=== Verification hierarchy

The paper's results have distinct epistemological status:
- *GI-1 through GI-14d* (including lemmas): verified algebraic identities and dynamical consequences, machine-checked with Z3. Lemma L4 (monotonicity of $ln$) is axiomatized as a property of the real logarithm; the assembly is verified relative to this axiom. The GI-13 identities verify algebraic preconditions for conservation (force symmetry, COM identity); the dynamical statements $dot(P) = 0$ and $dot(L) = 0$ follow analytically but are not separately encoded. The no-escape theorem (GI-14d) covers both radial ($L = 0$) and nonradial ($L > 0$) motion.
- *GI-14e*: physical consequence of Hamiltonian mechanics (analytical, not Z3). The Liouville theorem is not encoded in Kleis.
- *GI-15a*: mathematical asymptotic result (= GI-14d, Z3-verified).
- *GI-15b--d*: physics inference and observational evidence (not Z3).")


// =============================================================================
// Section 13: Conclusion
// =============================================================================

define sec_conclusion = ArxivSection("Conclusion",
"The POT gravitational kernel $K(r) = -A\/r + B ln(r\/r_0)$ admits a well-defined multi-center mechanics once the $R_3^+$ interaction charge is properly introduced:

#enum(
  [The naive coupling $U = M_1 M_2 K(r)$ is inconsistent with the BTFR. Reciprocity and the BTFR together require an additional coupling structure; a symmetric bilinear $R_3^+$ charge $C_(i j) = gamma q_i q_j$ is the minimal factorized bilinear realization and makes action--reaction symmetry structural.],
  [Conservation of momentum and angular momentum follow from the symmetry of the pair interaction.],
  [The circular-orbit condition is quadratic, yielding exactly one positive circular radius. Every such orbit is necessarily stable ($U_(\"eff\")'' = a\/r^3 + 2C\/r^2 > 0$).],
  [No finite-energy trajectory can escape to infinity under a globally valid $C ln r$ interaction. Observed scattering systems therefore constrain the $R_3^+$ domain.],
  [Coincident-source superposition is an exact instantaneous identity; the post-merger equilibrium source strength depends on the (not yet derived) coherence functional $cal(C)$. The BTFR overshoot would therefore constitute a potentially observable transient if the required recoherence dynamics occurs.],
  [The Bullet Cluster is an application, not the conceptual center. POT describes it with one parameter per component, deriving the isothermal sphere rather than assuming it.]
)

The paper's strongest result is not that POT accommodates the Bullet Cluster, but that *multi-center dynamics feeds a nontrivial constraint back into the foundational theory*: the isolated logarithmic $R_3^+$ interaction cannot extend unchanged to arbitrarily large separation. The interaction theory thus determines a constraint on the previously unspecified large-distance completion of the single-center theory.")


// =============================================================================
// Acknowledgments and References
// =============================================================================

define ack = ArxivAcknowledgments("The formal verification infrastructure was built using Kleis (https://kleis.io) with Z3 as the backend SMT solver.")

define ref_paper1 = ArxivReference("atik2026a", "Atik, E. (2026). Adjacent Causal Riesz Sectors and Galactic Gravity in Projected Ontology Theory. Preprint, Kleis Research.")
define ref_paper0 = ArxivReference("atik2026frc", "Atik, E. (2026). Flat Galactic Rotation Curves as a Theorem of Projected Ontology. Preprint, Kleis Research.")
define ref_paper2 = ArxivReference("atik2026b", "Atik, E. (2026). Emergent Conformal Geometry from Projection Kernels. Kleis Research.")
define ref_paper3 = ArxivReference("atik2026c", "Atik, E. (2026). Projective Metric Reconstruction from Coherent Kernel Families. Kleis Research.")
define ref_clowe2006 = ArxivReference("clowe2006", "Clowe, D., et al. (2006). A Direct Empirical Proof of the Existence of Dark Matter. ApJ, 648, L109.")
define ref_bradac2006 = ArxivReference("bradac2006", "Bradac, M., et al. (2006). Strong and Weak Lensing United. III. ApJ, 652, 937.")
define ref_cha2025 = ArxivReference("cha2025", "Cha, S., et al. (2025). A High-Caliber View of the Bullet Cluster through JWST. ApJ, 987, L15.")
define ref_cho2025 = ArxivReference("cho2025", "Cho, B. Y., et al. (2025). Joint JWST-DECam Lensing Reveals That the Bullet Cluster Is a Minor Merger. arXiv:2512.03150.")
define ref_springel2007 = ArxivReference("springel2007", "Springel, V., & Farrar, G. R. (2007). The speed of the 'bullet' in 1E0657-56. MNRAS, 380, 911.")
define ref_lage2014 = ArxivReference("lage2014", "Lage, C., & Farrar, G. (2014). Constrained Simulation of the Bullet Cluster. ApJ, 787, 144.")
define ref_randall2008 = ArxivReference("randall2008", "Randall, S. W., et al. (2008). Constraints on Self-Interaction of Dark Matter. ApJ, 679, 1173.")
define ref_angus2007 = ArxivReference("angus2007", "Angus, G. W., et al. (2007). On the Proof of Dark Matter, the Law of Gravity, and the Mass of Neutrinos. ApJ, 654, L13.")
define ref_zhang2026 = ArxivReference("zhang2026", "Zhang, X., et al. (2026). MOND with the IGIMF: Revisiting the Bullet Cluster with JWST. Phys. Rev. D.")
define ref_mcgaugh2000 = ArxivReference("mcgaugh2000", "McGaugh, S. S., et al. (2000). The baryonic Tully-Fisher relation. ApJL, 533, L99.")
define ref_kneib2011 = ArxivReference("kneib2011", "Kneib, J.-P., & Natarajan, P. (2011). Cluster lenses. A&A Rev., 19, 47.")
define ref_umetsu2011 = ArxivReference("umetsu2011", "Umetsu, K., et al. (2011). Cluster Mass Profiles from a Bayesian Analysis of Weak-Lensing Distortion and Magnification. ApJ, 729, 127.")
define ref_lyndbell = ArxivReference("lyndenbell1967", "Lynden-Bell, D. (1967). Statistical mechanics of violent relaxation in stellar systems. MNRAS, 136, 101.")
define ref_young2018 = ArxivReference("young2018", "Young, L. M., Williams, L. L. R., & Hjorth, J. (2018). Dynamics of merging: post-merger mixing and relaxation of an Illustris galaxy. JCAP, 2018(02), 033.")
define ref_nipoti2007 = ArxivReference("nipoti2007", "Nipoti, C., Londrillo, P., & Ciotti, L. (2007). Galaxy merging in MOND. ApJ, 660, 256.")
define ref_moura = ArxivReference("moura2008", "de Moura, L., & Bjorner, N. (2008). Z3: An efficient SMT solver. TACAS 2008.")


// =============================================================================
// Appendix
// =============================================================================

define sec_appendix = ArxivSection("Appendix: Verified Result Summary",
"#figure(
  table(columns: 3, stroke: 0.5pt, inset: 6pt,
    table.header([*Theory file*], [*Verified identities*], [*Primary content*]),
    [pot_galactic_interaction], [62], [Superposition, charge, conservation, two-body dynamics],
    [pot_galactic_interaction_numerical], [20], [Bullet Cluster dimensionless comparisons],
  ),
  caption: [Machine-checked algebraic identities (Z3 backend).],
  kind: table,
)

=== Verification hierarchy

#figure(
  table(columns: 3, stroke: 0.5pt, inset: 6pt,
    table.header([*Result*], [*Status*], [*Content*]),
    [GI-1 to GI-11], [Verified identity], [Superposition (coincident-source), BTFR overshoot (transient), phantom density, GR-inferred convergence, offset],
    [GI-12a to GI-12f], [Verified identity], [Symmetric pair interaction, crossover, force law],
    [GI-13a to GI-13d], [Verified identity], [Force symmetry, COM, acceleration inequality, pair count (algebraic preconditions for conservation; dynamical $dot(P) = 0$, $dot(L) = 0$ follow analytically)],
    [GI-14a to GI-14c], [Verified identity], [Reduced mass, unique positive circular radius (product of roots $< 0$), necessary stability],
    [GI-14d (L1--L4 + assembly)], [Verified (axiom: ln monotone)], [$forall E exists R forall r > R: U_(\"eff\")(r) > E$ ($L >= 0$)],
    [GI-14e], [Analytical], [Thermalization requires many-body (Liouville)],
    [GI-15a (= GI-14d)], [Verified identity], [No scattering at infinity for global $C > 0$],
    [GI-15b--d], [Physics/observational], [Finite-domain inference and evidence],
  ),
  caption: [Epistemological status of each result.],
  kind: table,
)

=== Theoretical spine

$
\"Riesz linearity\" arrow \"symmetric sector charges\" (M_i, q_i) arrow \"central pair force\" arrow \"conservation laws\"
$
$
arrow \"unique stable circular orbit\" arrow \"no asymptotic escape for global\" R_3^+ arrow \"finite-domain constraint.\"
$")


// =============================================================================
// Assemble and Compile
// =============================================================================

define all_elements = [
    sec_intro,
    sec_assumptions,
    sec_superposition,
    sec_phantom,
    sec_charge,
    sec_conservation,
    sec_twobody,
    sec_transition,
    sec_lensing,
    sec_bullet,
    sec_comparison,
    sec_limitations,
    sec_conclusion,
    ack,
    sec_appendix,
    ref_paper1, ref_paper0, ref_paper2, ref_paper3,
    ref_clowe2006, ref_bradac2006, ref_cha2025,
    ref_cho2025, ref_springel2007, ref_lage2014, ref_randall2008,
    ref_angus2007, ref_zhang2026, ref_mcgaugh2000, ref_kneib2011,
    ref_umetsu2011, ref_lyndbell, ref_young2018, ref_nipoti2007, ref_moura
]

define galactic_interaction_paper = arxiv_paper(
    paper_title,
    paper_authors,
    paper_affiliations,
    paper_abstract,
    paper_keywords,
    all_elements
)

// =============================================================================
// Paper Compilation
// =============================================================================

example "compile" {
    let typst_output = compile_arxiv_paper(galactic_interaction_paper) in
    out(typst_raw(typst_output))
}

example "validate" {
    assert(valid_arxiv_paper(galactic_interaction_paper) = true)
    out("Multi-center dynamics paper is valid!")
}