Addressing Open Problem O4 — the Whitney A₁ fold radius for galactic attractors
The dm³ contact ODE on the manifold $(\mathbb{R}^3,\, \alpha = dz - r^2\,d\theta)$ is (Grossi 2026, Vol IV §2):
with coupling parameter $\varepsilon = 2$ for the electromagnetic (macro) realization and attractor $\Gamma = \{r = 1,\, \dot\theta = 1,\, \dot z = 1\}$. The inner basin boundary is:
The value $r^*(2, 0) = 0.77594058$ is certified by certify_rstar.py (DOP853, rtol=$10^{-12}$, bisection tolerance $10^{-7}$). No closed-form expression for $r^*(\varepsilon, z_0)$ is known for any $\varepsilon \neq 0$ [Conjecture 2.2 of rstar_compute.py].
Open Problem O4 (WP57, §7) asks: what is the analogous fold radius at galactic scale ($g^{96}$), and can it be derived from the contact geometry?
We compute $r^*(\varepsilon, z_0 = 0)$ for $\varepsilon \in \{0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 4.0\}$ by bisection on the dm³ ODE (DOP853, rtol=$10^{-13}$, atol=$10^{-15}$, $T=1000$, bisection tolerance $5 \times 10^{-10}$).
| ε | r*(ε) | Note |
|---|---|---|
| 0.5 | 0.3779748572 | — |
| 1.0 | 0.5722354585 | near 1/√3 = 0.5774 (Δ = 5.2×10⁻³, not equal) |
| 1.5 | 0.6931303713 | near ln 2 = 0.6931472 (Δ = 1.7×10⁻⁵ — tantalizing, not proved equal) |
| 2.0 | 0.7759405754 | electromagnetic realization — certified [VERIFIED] |
| 2.5 | 0.8353301711 | — |
| 3.0 | 0.8789150842 | — |
| 4.0 | 0.9352791770 | — |
$r^*(\varepsilon)$ is monotone increasing: larger coupling pushes the basin boundary outward. The function is smooth and bounded in $(0,1)$. The limiting behaviors:
Fig. 1 — r*(ε): basin boundary as a function of coupling strength ε = ν/κ · bisection on dm³ ODE (DOP853) · green: MW at ε ≈ 1.61 · blue: 2:1 resonance at ε = 2
At $g^{96}$, the G-chain $G = U \circ F \circ K \circ C$ acts on a single galactic disc modeled by $(\mathbb{R}^3, \alpha = dz - r^2\,d\theta)$ at the disc scale (Grossi 2026, chGravity-scales.html). The coordinate identifications at galactic scale:
| Variable | Galactic meaning | Normalized to |
|---|---|---|
| r | Cylindrical radius from galactic centre | $R_{\rm disc}$ — disc scale radius; attractor at $r=1$ |
| θ | Orbital phase (azimuth) | $2\pi$ — one full orbit |
| z | Normalized gravitational potential along the orbit | Contact condition $\dot z = r^2$ sets the scale |
Under these identifications, the inner basin $\{r < r^*\}$ is the bulge — stars that orbit within the fold and are drawn toward the galactic center — and the outer basin $\{r > r^*\}$ is the disc, attracted toward the attractor at $r = 1$ (the exponential disc scale). This matches the geometric description in chGravity-scales.html: "the bulge is the inner basin $(r < r^*)$, the disc is the outer basin $(r > r^*)$."
The coupling parameter $\varepsilon$ in the galactic ODE controls the strength of the $z$-coupling: how strongly vertical position affects radial dynamics. In a galactic disc, vertical and radial motions are coupled through the potential. The natural dimensionless ratio is:
Using the galactic identification, we compute $r^*(\varepsilon_{\rm gal})$ for the Milky Way. Frequencies are evaluated at multiple radii to show the radial dependence:
| R (kpc) | κ (km/s/kpc) | ν (km/s/kpc) | ε = ν/κ | r*(ε) | R*_phys (kpc) |
|---|---|---|---|---|---|
| 2.0 | 155.6 | 177 | 1.138 | 0.610 | 2.14 |
| 3.5 (R_disc) | 88.9 | 143 | 1.607 | 0.713 | 2.50 |
| 5.0 | 62.2 | 115 | 1.853 | 0.755 | 2.64 |
| 8.2 (R_☉) | 37.9 | 73 | 1.924 | 0.765 | 2.68 |
| — | any κ | 2κ (exact) | 2.000 | 0.77594058 | 0.77594058 × R_disc |
κ computed as √2 × v_c/R for flat rotation curve v_c = 220 km/s. ν scaled from solar-circle value ν(R_☉ = 8.2 kpc) = 73 km/s/kpc (Bienaymé et al. 2006) using ν ∝ exp[(R_☉ − R)/(2h_R)] with h_R = 3.5 kpc. R*_phys = r*(ε) × R_disc = r*(ε) × 3.5 kpc.
Observationally, the MW bulge stellar population ends at $\sim 2$–$3\,\rm kpc$ (Portail et al. 2015; Bovy et al. 2019). The predicted $R^* \approx 2.5\,\rm kpc$ sits within this range.
The exact value $r^*_{\rm gal} = r^*_{\rm macro} = 0.77594058$ — the universality conjectured in chGravity-scales.html — holds when $\varepsilon_{\rm gal} = 2$ exactly. Under the proposed identification, this requires:
The condition $\nu = 2\kappa$ is the 2:1 vertical resonance: stars complete two vertical oscillations for every one epicyclic (radial) oscillation. This resonance is well-studied in galactic dynamics. It is the mechanism by which bar-driven orbital heating creates boxy/peanut-shaped bulges — stars captured in the 2:1 resonance family form the X-shaped structure observed edge-on (Combes & Sanders 1981; Bureau & Athanassoula 2005; Quillen 2002). The Milky Way's boxy/peanut bulge extends to $\sim 2.5$–$3\,\rm kpc$, and the 2:1 resonance radius has been located at $R_{2:1} \approx 2$–$3\,\rm kpc$ for the MW bar (Portail et al. 2017).
For a flat rotation curve, $\kappa = \sqrt{2}\,\Omega$ exactly. The 2:1 resonance condition $\nu = 2\kappa = 2\sqrt{2}\,\Omega$ then gives:
This is a constraint on the midplane density $\rho_{\rm mid}$ — satisfied at a specific radius $R_{2:1}$ in each galaxy. For the G-chain to give $r^*_{\rm gal} = 0.77594058$ universally, the disc scale radius $R_{\rm disc}$ must coincide with $R_{2:1}$. Whether this is generic, or specific to $\Lambda$CDM-typical disc galaxies, is an open question.
The analysis above is a derivation attempt, not a theorem. The logical structure is:
Step (3) is now proved (Lemma 9.1) — this eliminates the entire class of contact-geometric derivations as the proof route. The remaining open steps are (4) and (5): showing the galactic ODE has the dm³ form and identifying ε_gal from the orbital mechanics. These are claims about galactic physics, not contact geometry.
The CONJECTURE in chGravity-scales.html concerns the merger of two galactic $g^{96}$ attractors (Milkomeda), not the single-galaxy fold. The merger fold is a distinct object. The relevant boundary is the Hill sphere — the tidal radius of each galaxy in the field of the other — or equivalently the L₁ Lagrange point of the two-body gravitational system.
For the MW–Andromeda system ($M_{\rm MW} = 1.5\times10^{12}\,M_\odot$, $M_{\rm And} = 1.0\times10^{12}\,M_\odot$, separation $d_0 = 770\,\rm kpc$):
In normalized merger coordinates ($r = d/d_0$), the merger fold is at $r_{\rm merger}^* \approx 0.490$ — not $0.77594058$. The single-galaxy fold and the merger fold are distinct: the former is the bulge–disc boundary within each galaxy ($r$ in units of $R_{\rm disc}$); the latter is the tidal stripping boundary between the two galaxies ($r$ in units of $d_0$). The conjecture in chGravity-scales.html identifies both with the "same geometric object" magnified by scale — they are Whitney $A_1$ singularities, but in different coordinate systems and with different numerical values of $r^*/r_{\rm ref}$.
rstar_compute.py: no closed-form expression for $r^*(\varepsilon, z_0)$ is known for any $\varepsilon \neq 0$. The dm³ ODE is a non-autonomous system with no obvious integrating factor. Possible routes: (a) Lyapunov function bounding the basin; (b) matched asymptotic expansion in $\varepsilon$; (c) Poincaré section analysis reducing to a 1D map whose fixed point is $r^*$. A closed-form $r^*(\varepsilon)$ would complete O.4.2 and give a fully explicit expression for the galactic fold radius.Path 1 sought $r^*_{\rm gal}$ by deriving $\varepsilon$ from the contact structure of the galactic potential. §8 showed all three candidate routes (O.4.2a–c) fail or stall. The failure has a structural cause.
What Lemma 9.1 rules out. Any argument of the form "derive $r^*_{\rm gal}$ from the contact structure of $\Phi_{\rm eff}$" — including all three routes in O.4.2a–c — cannot close the proof, because $r^*$ does not live in the contact structure. Changing $\varepsilon$ moves $r^*$; replacing $\alpha$ entirely does not. This is why path 1 stalled: it was solving the wrong equation.
Path 2 — the correct task. Identify the ODE governing stellar orbit amplitudes at galactic scale in the normalized coordinates $(r = R/R_{\rm disc},\; z_c)$. Match its coefficients to the dm³ form. Read off $\varepsilon_{\rm gal}$. The contact form plays no role in this step. The three structural elements of the dm³ ODE and their galactic meaning:
| Term | Role in dm³ | Galactic interpretation | Status |
|---|---|---|---|
| r(1 − r²) | Conservative restoring force toward $r = 1$; linearizes to $-2(r-1)$ at the attractor | Centrifugal + gravitational balance near the circular orbit; linearizes to $-\kappa^2(r-1)$ at $R_{\rm disc}$. The nonlinear form differs ($-\kappa^2(r-1) \neq r(1-r^2)$ globally) — matching requires the full effective potential. | [MODEL] |
| e^{−z_c} | Exponential envelope — coupling decays as orbit winds up ($z_c$ grows monotonically on $\Gamma$) | Disc density or relaxation envelope: both decay exponentially in normalized orbital time $z_c = \int_0^t r^2\,d\theta$. Either interpretation gives $e^{-z_c}$. | [MODEL] |
| ε | Coupling amplitude — the only parameter that determines $r^*$ | The dimensionless ratio of vertical to radial oscillation strength at $R_{\rm disc}$. The candidate $\varepsilon_{\rm gal} = \nu/\kappa$ is the unique such ratio intrinsic to the disc potential — but confirming it requires deriving the ODE coefficient explicitly (O.4.2, restated below). | [OPEN] |
Elements (I) and (II) establish the existence of a fold. Element (III) alone determines its location. The contact form contributes nothing to any of the three.
The following derives the identification $\varepsilon_{\rm gal} = \nu/\kappa$ from galactic orbital mechanics. The contact form is not used — consistent with Lemma 9.1. All approximations are labelled.
Use azimuth $\theta = \phi$ as the independent variable and define $z_c = \Omega_c t$ (the guiding-centre clock, $\Omega_c = v_c/R_{\rm disc}$). Angular momentum conservation gives $\dot\phi = L/R^2 = \Omega_c/r^2$, so:
The contact form $\alpha = dz_c - r^2\,d\theta = 0$ is Kepler's second law written as a 1-form. It is a conclusion of the dynamics, not an input.
For a flat rotation curve ($v_c = \text{const}$, $\partial\Phi/\partial R = v_c^2/R$), the conservative radial force per unit azimuth, in normalized coordinates $r = R/R_{\rm disc}$, is exactly:
This is not a normal-form ansatz. It is the flat-rotation-curve effective force with the $r^4$ Jacobian of the $t\to\theta$ change of variable. The linearization at $r = 1$ gives $\partial_r[r(1-r^2)]|_1 = -2$; keep this number. Converting to first order requires:
So $e^{-z} = J_z/J_z(0)$: the fraction of initial vertical action still in the vertical mode. $z = 0$ is a full reservoir; $z \to \infty$ is vertically cold.
At the resonance $\ell\nu = m\kappa$ (resonant angle $\psi = \ell\theta_z - m\theta_r$), averaging the perturbation Hamiltonian over all non-resonant angles leaves only $H_{\rm res} = A(J_r, J_z)\cos\psi$. Hamilton's equations give:
One unit of vertical action destroyed creates $\nu/\kappa$ units of radial action. Therefore $\varepsilon \equiv \nu/\kappa$. The coupling enters as $\varepsilon(r-1)e^{-z}$ in $\dot r$ (vanishes on the circular orbit) and $-\varepsilon(r-1)^2 e^{-z}$ in $\dot z$ (back-reaction).
Linearizing the assembled ODE about $r = 1$ gives:
The circular orbit is marginally stable when $\varepsilon = 2$. This $2$ came from the radial force term in Step 11.2 ($\partial_r[r(1-r^2)]|_1 = -2$), which came from the flat-rotation-curve force $L^2/R^3 - v_c^2/R$. It has no relation to resonance ratios. Yet the 2:1 vertical resonance also gives $\nu/\kappa = 2$. Two completely independent inputs both land on 2 — this is the main internal check on $\varepsilon = \nu/\kappa$. Any other candidate ($\nu^2/\kappa^2$, $\kappa/\nu$, $\nu/(2\kappa)$, …) would place the marginal-stability threshold at a different resonance ratio.
Setting both $\dot r = 0$ and $\dot z = 0$ simultaneously, $\varepsilon$ cancels identically and the fixed-point $r$-coordinates satisfy:
These are heptagonal constants, independent of $\varepsilon$. But $r^* = 0.77594058$ is not a fixed point. It is the $z = 0$ crossing of the stable manifold of the saddle at $(r_1, z_{\rm saddle})$, where
The saddle's $r$-coordinate is $\varepsilon$-independent; its $z$-coordinate grows as $\ln\varepsilon$. As $\varepsilon$ increases, the saddle rises in $z$, the stable manifold tilts, and its $z = 0$ crossing $r^*(\varepsilon)$ increases — exactly matching the bisection table (§2). At $\varepsilon = 2$: $z_{\rm saddle} \approx 1.135$, and the manifold descends from $(0.445, 1.135)$ to cross $z = 0$ at $r^* \approx 0.776$.
| Approximation | Error order | Affects ε? |
|---|---|---|
| Flat rotation curve (input assumption) | — | sets κ and ν |
| ż = r² (angular momentum conservation) | exact | no |
| [A1] Epicyclic averaging ⟨r²⟩ → r² | O(a²) | no |
| [A2] Overdamped limit Γ ≫ Ω_c | O(Ω_c/Γ) | no |
| [A5] Reservoir closure J_z ∝ e^{−z}, single rate Γ | definition + one-process | no |
| [A6] Linear-response contact exchange | O((r−1)², z) | yes — weakest link |
| Single isolated resonance ℓν = mκ | resonance overlap → chaos | yes if violated |
Fig. 2 — Interactive phase portrait of the dm³ ODE · drag ε slider to move r* · Blue: disc basin (converges to attractor r = 1) · Red: escape basin (r → 0) · Yellow dashes: ṙ = 0 nullcline · Orange: saddle at r₁ = 2cos(3π/7) whose stable manifold crosses z = 0 at r* · Red dot: r* basin boundary
| Claim | Status |
|---|---|
| r* is a dynamical invariant of the ODE, independent of the contact form (Lemma 9.1) | [VERIFIED — α absent from ODE by inspection] |
| r*_gal = r*(ε_gal) given the galactic ODE has dm³ form (Theorem G2) | [MODEL | VERIFIED for ε=2] |
| dm³ ODE with ε=2 has r*(ε=2, z₀=0) = 0.77594058 | [VERIFIED — certify_rstar.py] |
| r*(ε) monotone in ε; table above to 10 s.f. | [VERIFIED — numerical, DOP853] |
| ε_gal = ν/κ at R_disc | [MODEL, DERIVED — §11: resonant action transfer J̇_r/J̇_z = −ν/κ + marginal-stability check; conditional on [A2], [A6]] |
| r* is z=0 crossing of stable manifold of saddle at r₁ = 2cos(3π/7) (Theorem G3) | [VERIFIED — phase-portrait analysis; r₁ from r³−r²−2r+1=0] |
| r*_MW ≈ 0.713, R*_phys ≈ 2.5 kpc | [MODEL] |
| Universality r*_gal = 0.77594058 ⟺ ν = 2κ at R_disc | [MODEL — follows from identification] |
| 2:1 resonance ⟺ boxy/peanut bulge boundary | [MODEL + lit support (Combes & Sanders 1981)] |
| r*(3/2) = ln 2 exactly | [OPEN — O.4.1] |
| Galactic ODE derivation from Φ_eff | [OPEN — O.4.2, key step] |
| Observational test R* ≈ 2.5 kpc | [OPEN — O.4.3] |
| Closed form r*(ε) for general ε | [OPEN — O.4.4] |
📄 Nota sobre o artigo PT/BR —
O PDF abaixo (Geometria de Contato e a Dobra Galáctica, norma ABNT NBR 6022)
é um artigo de síntese independente redigido para periódicos brasileiros.
Não é um capítulo deste livro nem um duplicado deste artigo de trabalho:
apresenta os resultados centrais de WP58 em português e propõe a identificação
de r* com a cáustica Whitney A1 de lentes gravitacionais —
conteúdo que, no corpo da série, está desenvolvido em dois lugares distintos:
·
ch-einstein.html · Einstein e a Forma de Contato
— § "A Dobra Whitney e as Lentes" (a identificação geométrica, uma página)
·
WP59 · Lentes de Matéria Escura
— tratamento completo: Teoremas T1–T3, previsões F1–F5, comparação com CDM/SIDM/FDM,
excesso Natarajan et al. (2025)
Open Problem O.4.4 of this paper lists “Lyapunov function bounding the basin” as one of three untaken routes. WP70 promotes it to a stated open problem (L1–L3) and marks every use of r* as a basin boundary elsewhere in the series as unlicensed until it is answered.
WP70 also adopts this paper’s lemma that r* is independent of the contact form, and draws the consequence: the missing Lyapunov function is a question in planar ODE theory, not contact geometry. And it generalises O.4.1 — since r*(ε) is monotone onto (0,1), it passes arbitrarily close to any target, so the ln 2 near-coincidence at ε = 1.5 is not evidence unless ε is fixed independently.