The Spacetime Resonance Theorem and the Origin of Dark Mass
Abstract
This work applies a well-known idea from acoustics and optics—multiple scattering in a medium—to gravity itself. In linearised general relativity the gravitational field is determined by a retarded Green function, but in most applications its “causal tail’’ (the part inside the light cone) is discarded. When this tail is kept and combined with the linear damped-oscillator susceptibility that arises from the kinematic density–sloshing of baryonic matter (mathematically identical to the response used in acoustical and optical scattering theory), spacetime behaves as an effective transport medium for weak curvature waves. From this simple starting point we derive a single theorem—the Spacetime Resonance Theorem—showing how coherent, time-varying baryonic motion generates a small resonant curvature field. Multiple gravitational scattering causes part of this oscillatory field to accumulate as a slow, cycle-averaged component. This stored curvature energy behaves exactly like an effective dark-mass density. The resulting framework naturally reproduces the observed dark-mass profiles of spiral galaxies, ellipticals, clusters, and filaments; predicts “local silence’’ and “homogeneous silence’’ (matching the Bullet Cluster and the CMB); and remains fully consistent with standard GR. No new particles or modifications of gravity are introduced: dark-mass behaviour emerges from the resonant, multiple-scattering response of baryonic structure within ordinary weak-field GR.
Full text
The Spacetime Resonance Theorem and the Origin of Dark Mass Bryan Hennelly ∗1,2,3 1Department of Electronic Engineering, Maynooth University, Maynooth, Co. Kildare, Ireland 2Department of Computer Science, Maynooth University, Maynooth, Co. Kildare, Ireland 3The Hamilton Institute, Maynooth University, Maynooth, Co. Kildare, Ireland 19 November 2025 Abstract In linearised general relativity the gravitational field is determined by a retarded Green function, but in most applications its “causal tail” is discarded. When this tail is kept and combined with the linear dampedoscillator susceptibility that arises from the kinematic density–sloshing of baryonic matter (mathematically identical to the response used in acoustical and optical scattering theory), this yields an effective–medium description of the weak–field gravitational response. Coherent, time–varying kinematic sloshing drives a collective curvature field whose delayed rescattering produces a small but persistent cycle–averaged component. This stored curvature energy appears as an effective dark–mass density and reproduces the observed dark–mass phenomenology. No new fields, particles, or modifications of GR are introduced: all dynamics arise from standard linearised gravity once the usual quasistatic truncation of the retarded kernel is avoided. The resulting causal–resonant framework reproduces the effective mass distributions of spiral galaxies, ellipticals, clusters, and filaments; yields natural “local” and “homogeneous” silence; and remains fully consistent with Solar–System tests and with the CMB. In this picture, dark–mass behaviour emerges as curvature energy stored by the collective resonant response of baryonic structure. ∗Email: [email protected] 1
Contents 1 Introduction 5 2 Spacetime Resonance Theorem 6 2.1 The Fundamental Oscillator–Field Principle . . . . . . . . . . . . 6 2.2 The Spacetime Resonance Theorem . . . . . . . . . . . . . . . . . 9 2.3 Corollary: Density–Controlled Restoring Frequency . . . . . . . . 12 2.3.1 Spectral origin of the density scaling . . . . . . . . . . . . 12 2.3.2 Agreement with known weak-field gravitational oscillators 13 2.4 Physical Interpretation and Preview . . . . . . . . . . . . . . . . . 13 3 Dark Mass 15 4 Resonant Field Dynamics 15 4.1 ChannelSeparation.......................... 15 4.2 Spacetime as an Effective Reversible Medium . . . . . . . . . . . 17 4.3 Collective Mode Dynamics and the Masson . . . . . . . . . . . . . 18 4.4 Linear Spatial Perturbations and the Propagation Law . . . . . . 19 5 Gravitational Scattering Theory 21 5.1 Single–Scatterer Response . . . . . . . . . . . . . . . . . . . . . . 22 5.2 Multiple Scattering and the Foldy–Lax Hierarchy . . . . . . . . . 23 5.3 Effective–Medium and Ensemble Averaging . . . . . . . . . . . . . 24 5.4 Transport and Diffusion Limits . . . . . . . . . . . . . . . . . . . 25 5.5 Unified Interpretation . . . . . . . . . . . . . . . . . . . . . . . . . 26 6 Spatiotemporal Feedback and Resonant Coherence 27 6.1 Constituent Entities & States of the Feedback System . . . . . . . 27 6.2 Feedback Channels and Interactions . . . . . . . . . . . . . . . . . 29 6.3 Temporal Coherence and Mutual Susceptibility . . . . . . . . . . 32 6.3.1 Phase Locking and Mutual Susceptibility . . . . . . . . . . 32 6.3.2 Global Frequency Selection in Inhomogeneous Baryonic Systems ............................ 34 6.3.3 Polarisation of Temporal Coherence . . . . . . . . . . . . . 34 7 Spatial Coherence of the Resonant Channel 36 7.1 Kirchhoff–Sommerfeld Representation . . . . . . . . . . . . . . . . 37 7.2 Effective Potential and Energy Storage in the Resonant Medium . 39 7.3 Single Point Source in a Resonant Channel–1 Background . . . . 40 7.4 Interference of Two Point Sources . . . . . . . . . . . . . . . . . . 42 7.5 Temporal Memory and the Steady–State Reduction . . . . . . . . 46 2
8 Steady–State (Attractor) Solutions 47 8.1 Spiral Galaxies as Coherent m= 1 (Laguerre–Gaussian) Modes . . . . . . . . . . . . . . . . . . . . . 49 8.1.1 Observing the Steady–State Solution . . . . . . . . . . . . 49 8.1.2 Energy Source: Epicyclic Sloshing with Coherent Filtering 51 8.1.3 The First Moment . . . . . . . . . . . . . . . . . . . . . . 52 8.1.4 The Thin–Disk Model: Solving the 2D Helmholtz Equation RevealstheHalo ....................... 55 8.1.5 Extending the Thin–Disk to Three Dimensions . . . . . . . 57 8.1.6 Gain from Incoherent Foldy–Lax Scattering . . . . . . . . 58 8.1.7 Stored Energy, Dark–Mass Equivalence, and EnvironmentalGain ............................ 59 8.1.8 Numerical Results: 3D Thin–Disk Model . . . . . . . . . . 60 8.1.9 Physical Interpretation . . . . . . . . . . . . . . . . . . . . 64 8.2 Ellipticals as Random Resonant Media . . . . . . . . . . . . . . . 68 8.2.1 Observing the Statistical Steady State . . . . . . . . . . . 68 8.2.2 Complex Baryonic Forcing with Memory . . . . . . . . . . 70 8.2.3 Four–dimensional speckle geometry of the gravitational field 72 8.2.4 Ensemble reduction: masson–speckle summation . . . . . . 74 8.2.5 Transport Mean Free Path and Collective Damping . . . . 76 8.2.6 Core Regularization in the Non–Transport Regime . . . . 78 8.2.7 Gradient–Energy Boost and Final Transport–Regularized HaloLaw ........................... 79 8.2.8 Numerical Results . . . . . . . . . . . . . . . . . . . . . . 80 8.2.9 Halo Formation from Coherent Energy Flux . . . . . . . . 81 8.2.10 Physical Interpretation . . . . . . . . . . . . . . . . . . . . 84 8.3 Clusters as Radiative Transport Extensions of Ellipticals . . . . . 88 8.3.1 Observing the Statistical Transport Equilibrium . . . . . . 88 8.3.2 Mathematical Formulation: Overlapping Elliptical TransportEnvelopes ........................ 91 8.3.3 Numerical Results . . . . . . . . . . . . . . . . . . . . . . 95 8.3.4 Physical Interpretation . . . . . . . . . . . . . . . . . . . . 98 9 Local Silence and the Preservation of Classical Gravity 101 10 Prediction: Absence of Gain in Homogeneous Systems 102 10.1 Three Distinct Regimes for Homogeneous Clouds . . . . . . . . . 103 10.2 Theorem: Homogeneous Media Produce No Resonant Gain . . . . 103 10.3 Interpretation and Connection to Perturbation Theory . . . . . . 104 11 Emergence of Multiple Scattering from a Weakly Inhomogeneous Continuous Medium 105 11.1 A Brief Review of Garvitational Multiple Scatter Theory . . . . . 105 11.2 The Resonant Frequency in a Continuous Medium . . . . . . . . . 107 11.3 Weak Random Potential and the Ballistic Regime . . . . . . . . . 107 3
11.4 Nonlinear Growth Toward the Transport Threshold . . . . . . . . 108 11.5 Emergent Resonant Patches and Effective Discreteness . . . . . . 108 11.6 Linear Limit and Connection to the CMB Power Spectrum . . . . 109 12 Conclusion 110 13 Dedication and Acknowledgement 111 A Appendix: Bessel and Hankel Functions Used in the Main Text 113 A.1 Cylindrical Bessel and Hankel Functions . . . . . . . . . . . . . . 113 A.2 Spherical Bessel Functions . . . . . . . . . . . . . . . . . . . . . . 114 A.3 Physical Interpretation . . . . . . . . . . . . . . . . . . . . . . . . 114 B Appendix: Code for Simulating the 3D Thin–Disk Resonant Halo 115 C Appendix: Code for Simulating Elliptical Galaxies - Transport Theory 124 D Appendix: Code for Simulating Cluster Galaxies 134 4
1 Introduction The empirical success of Newtonian and Einsteinian gravity [ 1 , 2 ] on laboratory, Solar–System, and stellar scales is beyond dispute. Yet across galaxies, clusters, and cosmological environments, the observed gravitational field exceeds the response predicted from baryons alone. Rather than postulating a new particulate component, we revisit a well–known but rarely exploited aspect of general relativity: the full retarded Green–function structure of the linearised field. In the weak–field limit, curvature responds through the causal kernel Gret , whose support extends throughout the past light cone. These interior contributions, the “tails”, are genuine features of GR but are normally discarded when the field is approximated by its instantaneous Poisson limit. Restoring the full retarded response exposes a small but physically meaningful effect: curvature reacts to time–varying matter through a delayed, frequency–dependent, finite-memory response that is already present in the linearised theory. A central point—formalised later by the Spacetime Resonance Theorem (Section 2)—is that we do not introduce an oscillator model for the mass distribution. The oscillator structure emerges automatically from two generic ingredients already present in weak–field gravity: (i) baryonic elements undergo small, kinematic density–sloshing driven by their orbital motion, epicyclic drift, and shear, so that each element carries a weak but persistent time–dependent component of the density; and (ii) these elements communicate through the causal retarded gravitational kernel, which propagates curvature perturbations between different locations with finite delay. Whenever spatially distributed masses exhibit such kinematic oscillations and interact through a retarded field—whether acoustic, electromagnetic, elastic, plasma, or gravitational—linearity and finite propagation speed inevitably generate phase lag, delayed mutual forcing, and therefore acollective damped–oscillator susceptibility. Weak–field GR inherits this same universal structure once the Poisson (truly instantaneous) limit is lifted. Baryonic elements therefore behave as gravitational scatterers: each responds to incident curvature with a linear susceptibility χ⋆ ( ω ) and re–emits a delayed perturbation. When many such elements participate, their retarded interactions produce a gravitational analogue of multiple scattering. The resulting ensemble renormalises the propagator in the standard Dyson sense, yielding an effective medium with a collective restoring frequency Ω, damping rate Γ, and propagation speed ceff . The coherent curvature response, Φ res , is purely dynamical. Because the retarded kernel carries long memory, its oscillatory part does not cancel cycle–to–cycle. A slowly accumulated, reversible curvature energy builds up, defining an effective mass density ρeff ∝ ⟨| Φ res|2⟩ that sources the ordinary Poisson equation. Thus all local tests of gravity remain unchanged, while extended systems (spirals, ellipticals, and clusters) develop substantial effective dark–mass envelopes from stored curvature energy. In this multiple–scattering framework, coherent baryonic motion in rotating disks, orbiting cores, or clustered environments drives constructive rescattering 5
of curvature. Depending on spatial contrast and coherence length, this produces two distinct observational regimes: large–scale coherent modes in spirals, and diffusive effective–medium halos in ellipticals and clusters. Homogeneous or weakly perturbed environments, lacking spatial contrast, exhibit no rescattering and remain resonantly silent, a fact consistent with the early Universe and the CMB. The present work develops this causal–resonant picture from first principles. Beginning with the retarded solution of linearised GR and the linear susceptibility of baryonic matter, we derive a self–consistent effective medium that accounts for spiral galaxies, ellipticals, clusters, and cosmological structure without introducing new fields or non–baryonic particles. The resulting Spacetime Resonance Theorem provides the microscopic foundation for the collective resonant mode and the stored curvature energy that behaves observationally as dark mass. 2 Spacetime Resonance Theorem The central result of this work rests upon a structural feature shared by all linear oscillator ensembles embedded in any causal medium described by a retarded Green function. Such behaviour is familiar across physics: acoustic scatterers in fluids, electromagnetic dipoles in dielectrics, plasma oscillators, elastic and seismological modes, collective excitations in condensed matter, and optical systems such as mode–locked lasers, where many cavity modes synchronise through retarded electromagnetic interactions. To make this universality explicit, we begin by formulating a general principle, independent of gravitational assumptions, that codifies the emergence of collective frequencies, phase coherence, renormalised propagation, and the universal damped–oscillator susceptibility produced by causal retardation and spatially distributed sources. Weak–field general relativity will then appear as a direct and natural specialisation of this structure. We emphasise that this Principle is not a theorem in the mathematical sense; it is a structural summary of universal behaviour in linear oscillator–field systems, intended as a conceptual framework for the rigorous results that follow. 2.1 The Fundamental Oscillator–Field Principle Principle (Fundamental Oscillator–Field Principle). Let {qi ( t ) } be an ensemble of linear oscillators (mechanical, electromagnetic, elastic, gravitational, or otherwise) that interact through a mediating field whose propagation is described by a causal retarded Green function Gret ( x, x′ ). Assume only linearity, causality, and that the oscillatory elements are spatially distributed, so that the field generated by one element is rescattered at distinct locations. Under these assumptions the coupled oscillator–field system exhibits the following universal structural features of linear response theory. 1. Collective behaviour and spectral emergence. 6
Retarded field-mediated interactions couple the phases and amplitudes of the oscillators. In the frequency domain the interaction is encoded in the coupling operator Kij(ω) = gi(ω)Gret(ω;xi,xj)gj(ω),(1) where gi ( ω ) denotes the linear coupling strength of oscillator i to the field. When the effective coupling bandwidth is sufficiently large compared to the natural detuning spread, the ensemble can exhibit collective long-time behaviour governed by the spectral structure of K(ω). 2. Damped–oscillator susceptibility. Each isolated oscillator responds linearly to a harmonic component of the mediating field. The susceptibility of an oscillator with natural frequency ω0and local damping coefficient γis the Lorentzian χ(ω) = 1 ω2 0−ω2−iγ ω,(2) the universal analytic structure enforced by linearity, causality, and the second–order local dynamics of the oscillator. This form appears identically across mechanical, electromagnetic, plasma, elastic, and weak-field gravitational systems. 3. Multiple scattering and self–energy renormalisation. Because the oscillators are spatially distributed, the field they source is rescattered from multiple locations. The net effect modifies the propagator through the Dyson relation G−1 eff (ω, k) = G−1 0(ω, k)−Σ(ω, k),(3) where G0 denotes the bare propagator of the mediating field and Σ( ω, k) is the self–energy encoding all multiple-scattering and feedback processes. In the long-wavelength, low-frequency regime, analyticity of the retarded response ensures an expansion Σ(ω, k)≃a0+a1(iω)+a2ω2+b2k2,(4) where a0 , a1 , a2 , and b2 are real coefficients determined by the statistics, spatial distribution, and coupling strengths of the oscillators. Substituting (4) into (3) yields a homogenised second–order effective evolution equation with collective parameters (Ω,Γ, ceff),(5) in which Ω is the emergent collective frequency, Γ the effective damping rate, and ceff an emergent propagation speed. 7
The above features represent a broad and robust structural pattern of linear response theory. They arise whenever spatially distributed oscillatory elements interact through a causal retarded field, independent of the specific physical medium. Spectral viewpoint. Because the field-mediated coupling (1) acts linearly on the vector of oscillator amplitudes, it may be regarded as an integral operator (or as a finite-dimensional matrix in discretised form). For an amplitude vector q= (q1, . . . , qN), the action of the operator is (K(ω)q)i= N X j=1 Kij(ω)qj,(6) or, in the continuous limit, ( Kq )( x ) = RK ( ω ; x, x′ ) q ( x′ ) d3x′ . The long-time dynamics of the coupled system are governed by the spectral properties of this operator. In particular, collective behaviour corresponds to the emergence of a normal mode associated with the dominant eigenvalue of K ( ω ), so that questions of uniqueness, stability, and coherence of the collective mode reduce to the spectral structure of the operator (6). Conditions for a Unique Dominant Collective Mode. For each fixed frequency ω , the coupling operator K ( ω ) acts linearly on the amplitude vector q and therefore possesses a (generally complex) spectrum defined by the eigenvalue equation K(ω)ψ=λ ψ, (7) with λ∈C and nonzero eigenfunction (or eigenvector) ψ . Because K ( ω ) is typically non-Hermitian, its eigenvalues need not be real. In linear response theory the long-time behaviour of the system is controlled by the eigenvalue whose real part ℜ ( λ ) is largest, since ℜ ( λ ) determines the temporal growth or decay rate of the associated collective mode. Thus ordering the spectrum by real part is the natural criterion for dynamical dominance and motivates the conditions below. To identify a unique collective normal mode and relate its frequency directly to the dominant spectral component of K ( ω ), one may further impose three conditions C1–C3: (C1) Irreducibility: no nontrivial partition of { 1 , . . . , N} blocks all cross-terms of Kij(ω). (C2) Dominant eigenvalue: K ( ω ) possesses a single eigenvalue whose real part strictly exceeds that of all others. (C3) Spectral gap: the dominant eigenvalue is separated from the remainder of the spectrum by a positive gap. When these conditions hold, the long-time behaviour is governed by a single collective mode whose emergent frequency is determined directly from this 8
dominant eigenvalue. In the gravitational application that follows, these spectral conditions are satisfied automatically by the positivity and compactness of the long-wavelength weak-field kernel. Principle 2.1 summarises a structural pattern so widespread that it is seldom stated explicitly: whenever spatially distributed oscillatory elements communicate through a causal retarded field, the system generically develops collective modes, a renormalised (effective) propagator, a well-defined self–energy encoding multiple scattering, and the universal Lorentzian susceptibility characteristic of linear second–order response. No reference to mechanical, electromagnetic, elastic, plasma, or gravitational specifics is required for the emergence of this general structure. 2.2 The Spacetime Resonance Theorem Before applying the Principle to gravity, we note that the weak-field Einstein equations satisfy its assumptions: the linearised curvature dynamics are governed by a causal retarded Green function, the equations are strictly linear in the perturbation, and baryonic matter provides spatially distributed oscillatory sources. Thus weak-field GR lies fully within the scope of Principle 2.1. When one does not impose the Newtonian instantaneous limit and instead retains the full retarded gravitational Green function, spacetime curvature plays the role of the causal medium; baryonic mass elements constitute the spatially distributed oscillatory sources; and gravitational rescattering generates the corresponding self–energy. The gravitational case therefore inherits the entire universal structure: collective behaviour, Lorentzian susceptibility, multiple scattering, Dyson renormalisation, and effective medium propagation. Before stating the theorem, we establish that the long-wavelength weak-field gravitational interaction automatically satisfies the spectral conditions (C1)–(C3) defined after Principle 2.1. Statement (Spectral Properties of the Weak–Field GR Kernel). In the long-wavelength, low-frequency regime, the linearised retarded gravitational kernel may be written schematically as the positive integral operator K(x, x′) = 4πG G0(x, x′)ρb(x′),(8) where G0 denotes the ω→ 0 (static, long-wavelength) limit of the retarded Green function of the linearised Einstein operator. In this limit the kernel reduces to the strictly positive Newtonian potential 1 / (4 π|x−x′| ), and it is in this static regime that positivity is invoked. We treat K as an integral operator acting on continuous functions over any bounded spatial domain. The restriction to a bounded spatial domain ensures that the Newtonian kernel 1 / (4 π|x−x′| ) defines a compact integral operator on this function space, as is standard in potential theory. 9
into its static and oscillatory frequency sectors, as established in the Spacetime Resonance Theorem. No approximation is involved: the split follows directly from the linearity of the retarded weak–field equations. The slowly varying component of the source generates the static Newtonian potential, while the phase–coherent oscillatory component excites the resonant curvature field, and the total solution is their exact superposition. Accordingly we decompose both the baryonic density and the curvature potential, ρb=ρb,0+ρb,1,Φ = ΦN+ Φres,(33) where the decomposition is functional rather than physical: ρb,1 denotes the narrow–band coherent part of the source and ρb,0 the slowly varying remainder. This exact frequency–sector split underlies the definitions of Channel 1 (the resonant, dynamic response) and Channel 2 (the quasi–static Newtonian response) used throughout the remainder of this section. Channel 1: The Resonant Channel. A useful intuition comes from a simple mechanical analogue: a slightly deformable water balloon attached to a rotating string. As it moves in a circle, the balloon does not maintain a rigid shape; variations in angular acceleration cause its interior to slosh, producing a coherent, phase–lagged deformation. In a rotating stellar disc, stars behave in an analogous way: because their orbital speed and centrifugal balance vary with radius, each stellar element undergoes small but coherent compressions and rarefactions along its orbit. This orbital sloshing generates the oscillatory density component ρb,1 that drives Channel 1. In pressure–supported systems (ellipticals and spheroids), individual stellar motions are random, but neighbouring stars still share locally correlated accelerations as they respond to the same smooth background potential. These correlated displacements produce a weak but coherent fluctuating density field: local patches “breathe” in and out together, providing the analogue of ρb,1 in systems without ordered rotation. Thus both discs and ellipticals produce the same type of phase–coherent forcing, differing only in the geometry of the susceptibility and coherence domains. Whenever such coherent sloshing is present, the frequency–dependent part of the retarded gravitational propagator becomes active. In the effective–medium limit it drives a collective, phase–coherent curvature response governed by ∂2 tΦres + 2Γ ∂tΦres + Ω2Φres −c2 eff∇2Φres = 4πG ρb,1.(34) Channel 1 therefore describes the resonant gravitational response of the stellar ensemble to these kinematic density fluctuations. It is active only when ρb,1 = 0 and vanishes identically in the static limit. Channel 2: The Newtonian Channel. If the source contains no phase– coherent time dependence, the frequency–dependent terms in Eq. (34) play no role. With ∂t Φ = 0 the retarded solution reduces to the usual instantaneous constraint ∇2ΦN= 4πG ρb,0,(35) 16
the standard Newtonian limit arising from the ω→ 0 behaviour of the retarded propagator. Thus Channel 2 is simply the static limit of the same underlying response. The two channels always coexist in general systems. Channel 1 carries the dynamic, phase–coherent component of the curvature field, while Channel 2 sets the quasi–static geodesic background on which that field propagates. Their separation in Eqs. (33) – (35) is therefore conceptual rather than ontological: it cleanly distinguishes the frequency–dependent part of the retarded response from the instantaneous part. Their slow mutual interaction—dynamic energy accumulation modifying the quasi–static background—will be analysed later in the feedback architecture of Section 6. The remainder of this book is devoted entirely to Channel 1. Section 4.2 reformulates its dynamics in terms of a standard effective–medium wave equation; Section 5 develops the corresponding single– and multiple–scattering theory. Channel 2 requires no further development until the feedback structure is assembled in Section 6. 4.2 Spacetime as an Effective Reversible Medium Channel 1 requires a clear physical interpretation. The aim of this section is not to attribute material properties to spacetime, but to clarify that once the full retarded Green function of linearised GR is retained, the effective response generated by coherent baryonic forcing and multiple gravitational scattering behaves mathematically like a reversible, compressible medium with a single longitudinal mode. The “medium” in this sense is not spacetime itself but the collective gravitational response of the baryonic ensemble. This viewpoint mirrors the standard effective–medium interpretation used in multiple–scattering acoustics (e.g. the experiments of Derode, Tourin & Fink [ 4 , 5 ]), where microscopic scattering produces a macroscopic, reversible wave equation with well–defined resonance, linewidth, and effective propagation speed. Whenever the baryonic source contains a coherent oscillatory component, the frequency–dependent part of the effective propagator (§2) excites a collective mode governed by ∂2 tΦres + 2Γ ∂tΦres + Ω2Φres −c2 eff∇2Φres = 4πG ρb,1,(36) the fundamental evolution law of Channel 1. To expose its structure, compare Eq. (36) with the standard damped acoustic wave equation in a compressible, inviscid medium, ∂2 tu+ 2Γ ∂tu+ Ω2u−c2 s∇2u=fext ρ0 ,(37) where u is the velocity potential. This parallel is purely mathematical: both systems are governed by the same second–order operator with a restoring frequency Ω, a phase–lag coefficient Γ, and a wave speed ( ceff or cs ). The structural 17
correspondence is therefore Φres ↔u, ceff ↔cs,4πG ρb,1↔fext/ρ0.(38) No physical identification is implied: the effective curvature field is not a material compression, but its governing equation shares the same mathematical form as a reversible acoustic medium, exactly as in the effective–medium formulation of multiple–scattering acoustics [4, 5]. This equivalence extends to the instantaneous quadratic energy densities: Eu=1 2ρ0h(∂tu)2+c2 s(∇u)2+Ω2u2i,EΦ=1 8πGh(∂tΦres)2+c2 eff(∇Φres)2+Ω2Φ2 resi, (39) each conserved in the limit Γ → 0. Under harmonic driving at frequency ω , the resonant field attains the steady–state amplitude |Φres|RMS =4πG |ρb,1| p(Ω2−ω2)2+ (2Γω)2, Q =Ω 2Γ,(40) precisely analogous to a driven acoustic resonator. High– Q environments (Γ ≪ Ω) support long–lived, nearly lossless oscillations of the gravitational potential— the collective curvature mode of the baryonic ensemble, directly mirroring the coherent modes observed in multiply scattering acoustic media [4, 5]. Fundamental Equivalence (mathematical). Channel 1—the resonant gravitational response—is governed by the same linear operator that appears in effective reversible acoustic media. The equivalence is structural rather than material: pressure stress in a fluid ↔ curvature stress in the effective gravitational medium. This structural correspondence provides a robust dynamical interpretation of Channel 1 and motivates the use of standard wave–medium tools (normal modes, phase lag, quality factors, effective sound speeds). These concepts form the foundation for the scattering and coherence theory developed in the next section. 4.3 Collective Mode Dynamics and the Masson Channel 1 admits a natural variational formulation, reflecting its character as a genuine wave sector of the effective gravitational response of the baryonic ensemble. In the absence of damping (Γ=0), the resonant field Φ res is governed by the conservative Lagrangian density Lres =1 8πGh(∂tΦres)2−c2 eff(∇Φres)2−Ω2Φ2 resi,(41) whose Euler–Lagrange equation is the undamped form of Eq. (36) . The three terms correspond respectively to inertial response, spatial stiffness, and the collective restoring force determined by the baryonic self–energy. 18
The homogeneous solutions of Eq. (41) form a discrete set of coherent oscillatory wave packets—the normal modes of the effective gravitational medium. In condensed–matter language, a normal mode of a reversible, multiple–scattering medium is called a phonon. By strict structural analogy (established in Section 4.2), the oscillatory curvature mode of Channel 1 is the gravitational counterpart of such a phonon. For clarity, and to distinguish it from the familiar metric gravitational wave of vacuum GR, we refer to this collective excitation as a masson. The term is descriptive rather than ontological: a masson is not a new particle species, but the coherent curvature mode of the effective gravitational propagator generated by the baryonic ensemble itself. Each masson carries phase and reversible energy, propagates with effective speed ceff , and couples to baryonic structure through the coherent susceptibility χ⋆(Ω) introduced in the scattering theory. The reversible energy associated with these modes is the quadratic curvature energy density already defined in Eq. (30) . This energy behaves gravitationally as an effective mass density and later feeds Channel 2 through the dark–mass mechanism developed in Section 6. In the driven system, the coefficient Γ controls linewidth and temporal phase coherence, with the corresponding quality factor Q=Ω 2Γ.(42) High– Q media (Γ ≪ Ω) support long–lived, spatially coherent curvature oscillations; macroscopic signatures of masson dynamics sustained by multiple gravitational scattering. The masson therefore provides the microscopic language for Channel 1: a collective curvature excitation of the effective gravitational medium. Its propagation, linewidth, and energy transport form the basis of the ensemble scattering and coherence theory developed in Section 5. Before constructing that scattering theory, we require the corresponding spatial propagation law for these modes. This follows from linear perturbation of Eq. (36) , expanding the resonant field about a quasi–static background and resolving it into spatial Fourier components. The resulting dispersion relation, group velocity, damping rate, and energy balance are summarised in the next subsection; these constitute the linear inputs for the multiple–scattering framework that follows. 4.4 Linear Spatial Perturbations and the Propagation Law In this section we summarise the spatial perturbation theory of the causal–resonant field. Nothing new is introduced here: this is the standard linearisation used in wave mechanics to extract dispersion, damping, and energy transport. Its purpose is simply to determine how the already–identified collective mode of Channel 1 19
propagates spatially as a weak gravitational wave within the effective medium established by multiple scattering. For completeness, the few algebraic steps that lead to the relations summarised below follow directly from inserting a small perturbation δ Φ into the effective-medium PDE of Channel 1 and performing the standard Fourier decomposition. Linearisation about a quasi–static background yields the scalar wave–type equation ∇2δΦ + 1 c2 eff −∂2 t−2Γ ∂t+ Ω2δΦ = 4πG c2 eff δρb,(43) where δ Φ is the weak spatial perturbation of the resonant field and δρb its small baryonic source. This perturbative split is purely mathematical and should not be confused with the physical temporal split ρb = ρb,0 + ρb,1 used to define the two gravitational channels; here we examine only the spatial wave content of the established collective mode. A Fourier mode δΦ∝ei(k·x−ωt)satisfies the dispersion relation ω2+ 2iΓω−Ω2=−c2 effk2,(44) from which ω2 r= Ω2+c2 effk2,|δΦ|∝e−Γt,(45) and the phase/group velocities follow as vph =ωr k, vg=c2 effk ωr .(46) Thus Channel 1 supports weakly damped spatial waves with natural frequency Ω and linewidth Γ, exactly as in reversible acoustic or optical media—consistent with the effective–medium analogy established in [4, 5]. The corresponding energy balance is ∂tEΦ+∇·SΦ=1 c2 eff δρb∂tδΦ−2Γ(∂tδΦ)2,(47) where the Γ–term represents the reversible phase–lag attenuation associated with the frequency–dependent self–energy of the effective propagator. Cycle–averaging yields the usual quality factor, Q=Ω 2Γ,(48) which measures temporal phase coherence. Under harmonic driving at frequency ωthe RMS response is |δΦ|RMS =4πG |δρb| p(Ω2−ω2)2+ (2Γω)2,(49) peaking strongly at ω≃Ω with bandwidth 2Γ. 20
For high driving frequencies ω≫Ω, |δΦ|RMS ∼4πG |δρb| ω2→0,(ω≫Ω),(50) demonstrating that the collective mode acts as a strict gravitational low–pass filter. Only frequencies near Ω can excite Channel 1; higher–frequency baryonic fluctuations are strongly suppressed. This mirrors the behaviour of ideal compressible media, where frequency response is controlled by the effective bulk modulus: high–frequency density perturbations cannot drive the medium. Together, Eqs. (43) – (49) show that the resonant field behaves as a weakly damped, energy–storing continuum supporting dispersive waves with finite bandwidth ∆ω= Ω/Q. The associated coherence time, τcoh ∼Q Ω,(51) links the microscopic phase–lag parameter Γ to the macroscopic persistence of curvature oscillations. As demonstrated in Section 6, the effective memory of the system extends to τmem ∼ Genvτcoh , providing the basis for sustained coherent amplification in Channel 1. These linear propagation and damping laws supply the foundation for the single– and multiple–scattering theory developed next in Section 5. 5 Gravitational Scattering Theory Having established that the dynamical (Channel 1) sector of gravity is governed by the effective resonant wave equation derived in the previous section, we now construct the corresponding gravitational scattering theory. Throughout this section we work exclusively with the oscillatory baryonic component ρb,1 , since only temporally varying sources can excite and sustain the collective resonant field Φ res . The static component ρb,0 generates the background Newtonian field but does not contribute to scattering. The goal in what follows is to develop the gravitational analogues of the Lippmann– Schwinger, Foldy–Lax, and Dyson equations using the propagator Gret associated with the effective-medium PDE. This provides the microscopic foundation for the collective behaviour analysed later—including coherent spiral modes and diffusive elliptical/cluster regimes (Section 8). In this framework each baryonic mass acts as a gravitational scatterer characterised by a linear susceptibility χ⋆ (Ω), which encodes its phase-lagged response to an oscillatory potential at frequency Ω. More generally, the gravitational response of a localized mass is a polarizability tensor χij ⋆ (Ω), whose symmetry is set by the local spacetime environment. In isotropic settings the tensor reduces to the scalar χ⋆ (Ω) used throughout this section, whereas in anisotropic geometries the effective medium develops preferred directions of response. This tensorial 21
structure will become essential when comparing planar spiral systems (with lateral polarization of coherent gain) to the isotropic coherence cells of ellipticals, but for the present microscopic theory the scalar form suffices. An incident perturbation δ Φ inc produces a delayed re-emission described by the retarded Green function Gret , and the resulting scattered field propagates to all other masses. Repeated re-scattering among many inclusions generates a hierarchy of coupled responses, precisely analogous to classical multiple-scattering theory in acoustics and electromagnetism [6, 7, 8]. Thus baryonic matter plays a dual role: its oscillatory component ρb,1 acts as the continuous driving source of the resonant field, while the same masses also act as discrete gravitational scatterers that re–emit curvature perturbations through their susceptibility χ⋆ (Ω). The iterative network of scattering and re-scattering links the microscopic delayed response of individual baryons to the emergent, large-scale propagation properties of the collective resonant field. It therefore supplies the necessary bridge between the effective wave operator introduced earlier and the macroscopic gravitational behaviour developed in the remainder of the book. 5.1 Single–Scatterer Response A discrete baryonic mass m⋆ with linear gravitational susceptibility χ⋆ (Ω) responds to an incident oscillatory potential Ainc through the standard Lippmann–Schwinger relation Asc(x)=t⋆(Ω) GΩ(x−x⋆)Ainc(x⋆), t⋆(Ω) = 4πG m⋆χ⋆(Ω),(52) where GΩ is the frequency–domain form of the retarded Green function associated with the effective–medium propagation equation. The quantity t⋆ (Ω) defines the gravitational scattering strength of the inclusion, capturing its phase–lagged re–emission of the incident curvature perturbation. In general, the susceptibility of a localized gravitating mass is a polarizability tensor χij ⋆ (Ω). Local geometric anisotropies (for example, planar kinematics or coherent streaming) can preferentially amplify particular tensor components, thereby polarizing the local gravitational response. In the isotropic approximation used here the tensor reduces to its trace, χ⋆ (Ω), but the full tensor structure will reappear in Section 8.1.9 (spiral galaxies) and Section 8.2.10 (elliptical galaxies). Because the resonant sector is weakly damped (Γ ≪ Ω) and the underlying dynamics are those of linearised GR, the interaction is reversible: no curvature energy is irreversibly absorbed by the scatterer. Instead, the inclusion stores curvature perturbations only transiently, re–emitting them with a phase delay determined by χ⋆(Ω) and the causal structure of GΩ. Equation (52) therefore provides the microscopic building block of the gravitational multiple–scattering theory developed below, linking the individual baryonic response to the collective, ensemble–level behaviour of the resonant field. 22
5.2 Multiple Scattering and the Foldy–Lax Hierarchy In classical wave physics, the interaction of many scatterers with a coherent incident field is described by the Foldy–Lax hierarchy [ 6 , 7 , 8 ], which provides the standard self–consistent formulation of multiple scattering among discrete inclusions. Originally developed for acoustics and later generalised to electromagnetism, elasticity, and photonics, the same framework applies directly to the present gravitational problem once the effective propagator GΩ has been specified by the linearised dynamics of Section 4.4. Throughout this subsection the field amplitude A denotes the oscillatory (resonant) curvature perturbation Φres driven by the baryonic component ρb,1. Each baryonic inclusion re–emits the incident curvature fluctuation with a phase delay determined by its susceptibility χ⋆ (Ω). These delayed re–emissions propagate to all other inclusions via the Green function GΩ , and the resulting network of mutual couplings constitutes a gravitational analogue of standard multiple– scattering theory. When the spacing between inclusions is much smaller than the resonant wavelength ( kΩd≪ 1), their phases remain locked to the incident field and the scattered amplitudes add coherently. In this fully coherent Foldy–Lax limit, Acoh ∝N|χ⋆(Ω)|Ainc,(53) so that the ensemble behaves as a single collective radiator with intensity Icoh ∝N2|χ⋆(Ω)|2|Ainc|2,(54) It is convenient to quantify this as an environmental gain Genv ≡|Atot|2 |Ainc|2,Genv =N2|χ⋆(Ω)|2(coherent limit).(55) For an ensemble of N scatterers, the total field satisfies the discrete Lippmann– Schwinger system A(x)=Ainc(x)+t(Ω) X n GΩ(x−xn)A(xn),(56) with t(Ω) = 4πGm⋆χ⋆(Ω). In operator notation this becomes A= (1 −tGΩ)−1Ainc,(57) whose Neumann expansion, A=Ainc +tGΩAinc + (tGΩ)2Ainc +··· ,(58) sums all orders of mutual gravitational re–scattering. In the limit Γ → 0 (purely real GΩ ), the response is strictly coherent and Eq. (57) reduces to the linear scaling of Eq. (53). 23
If instead the relative phases of different inclusions become effectively random—because of geometric delays, finite memory (Γ > 0), or dynamical phase drift—then scattered amplitudes no longer combine coherently. They add in intensity only: Itot =N|χ⋆(Ω)|2|Ainc|2,(59) yielding the standard incoherent scaling Genv =N|χ⋆(Ω)|2,pGenv =√N|χ⋆(Ω)|.(60) The spiral–disk regime discussed in Section 8.1 lies in the long–wavelength, phase– ordered limit ( kΩd≪ 1), where both spatial phasing and delayed temporal response organise coherently into a global m = 1 pattern. As coherence decreases—due to finite Γ, increasing disorder, or large separations—the system transitions smoothly toward the statistical regime treated in the next subsection, where ensemble averaging of Eq. (56) produces an effective–medium description appropriate to elliptical and cluster environments. This subsection has treated multiple scattering for a fixed arrangement of scatterers, with coherence determined purely by the phase relations among their delayed re–emissions. The following subsection extends the framework to statistically disordered or dynamically evolving ensembles, where averaging over many such configurations yields the emergent effective parameters of the resonant gravitational field. 5.3 Effective–Medium and Ensemble Averaging The Foldy–Lax hierarchy describes the fully coherent interaction of sub-wavelength baryonic scatterers with an incident oscillatory gravitational field. In realistic configurations—particularly in irregular, thick, or dynamically evolving systems—phase coherence among scatterers is degraded as separations approach the resonant wavelength or as differential motions introduce random phase delays. Under such conditions the local fields A (x n ) at different inclusions can no longer be treated as perfectly correlated, and the collective behaviour must be described statistically through an ensemble average over possible configurations. Averaging the discrete Lippmann–Schwinger equation (56) over many realisations of scatterer positions yields the standard effective–medium (Dyson) equation for the coherent mean field: (∇2+k2 eff)⟨A⟩= 0, k2 eff =k2+ Σ(Ω),(61) where Σ(Ω) is the ensemble self–energy. The self–energy encapsulates the average phase delay and finite dwell time associated with multiple gravitational rescatterings among baryonic inclusions. Its real part modifies the dispersion relation, while its imaginary part determines the ensemble attenuation rate: Γeff =ceff Im keff,(62) 24
which is reduced relative to the microscopic damping Γ whenever reversible rescattering retains phase coherence over extended domains. The degree of coherent amplification produced by the environment can be summarised using the environmental gain factor Genv , defined earlier in Eq. (60) as the total intensity gain. In the effective–medium limit this quantity can be conveniently written in terms of the coherence volume Vcell of the emergent Bloch-like mode: Genv ∝1 Vcell =1 ℓ2 ⊥ℓ∥ ,(63) where ( ℓ⊥, ℓ∥ ) denote the transverse and longitudinal coherence lengths of the averaged field. From a microscopic perspective, the local intensity gain in a coherent patch scales as the square of its average coupling efficiency divided by its coherence volume: Genv ∝|X|2 Vcell =|X|2 ℓ2 ⊥ℓ∥ ,(64) where X is the effective (ensemble–averaged) coupling amplitude. Large coherence domains (large ℓ⊥ and ℓ∥ ) thus imply strong environmental gain, while shrinking coherence lengths reduce Genv by increasing the number of independent scattering cells. Systems with strong spatial phase order (e.g. thin, rotating disks) possess large coherence volumes and hence large environmental gain. These satisfy Γeff ≃Γ,(65) recovering the coherent Foldy–Lax limit. Increasing disorder or dynamical randomness decreases the coherence lengths, lowers Genv , and drives the system toward regimes where ensemble damping is dominated by phase randomisation rather than the microscopic Γ. Equation (61) therefore provides the statistical extension of the coherent multiple–scattering framework. Ordered systems correspond to real, phase–aligned self–energies with large environmental gain, whereas disordered ensembles are described by complex, stochastic self–energies whose imaginary parts encode diffusive loss of coherence. This effective–medium formulation supplies the macroscopic link between the microscopic scattering interactions and the global intensity structure analysed later in Section 8.2. 5.4 Transport and Diffusion Limits When the relative phases of individual baryonic scatterers become random and time–dependent, the ensemble field correlation function obeys the standard radiative–transfer equation familiar from acoustic, electromagnetic, and elastic multiple–scattering theory [8]: ˆ s·∇I(x,ˆ s)=−I(x,ˆ s) ℓs +1 ℓsZP(ˆ s,ˆ s′)I(x,ˆ s′)dΩ′+S(x,ˆ s),(66) 25
coexistence of fast and slow regimes ensures that resonant dynamics, long–term memory, and the large–scale gravitational structure remain mutually consistent. The ensemble gain Genv clarifies that Φ res represents the coherent superposition of many baryonic emitters. Since | Φ res| ∝ √Genv , the stored component scales as ρeff ∝ Genv . This amplification is purely statistical—reflecting coherent phase alignment across many baryons—and does not violate causality; it follows from standard multiple–scattering theory applied to curvature perturbations within the weak–field limit. The causal relations illustrated schematically in Fig. 1 are expanded in the table below. For compactness, the figure merges the resonant field Φ res and the formation of the stored effective mass ρeff into a single step. In the table this process is separated into two stages: (i) a fast feed, in which the oscillatory effective curvature field deposits energy into the spacetime response encoded by resonant propagation, and (ii) a slow accumulation, in which this energy is integrated over the memory interval τmem = Γ −1 to form the time–averaged effective mass density ρeff . This makes explicit that propagation and rescattering operate on the wave timescale, while the long–term stored component arises only after temporal coarse–graining over the resonant memory kernel. 6.3 Temporal Coherence and Mutual Susceptibility The behaviour described in this subsection is the concrete gravitational realisation of the collective phase-locking and emergent frequency predicted by the Fundamental Oscillator–Medium Theorem (Section 2). 6.3.1 Phase Locking and Mutual Susceptibility The resonant channel linking baryonic oscillations ρb,1 to the effective curvature field Φ res behaves as a coupled oscillator system mediated entirely by causal propagation of curvature perturbations. Baryons drive the resonant field through the Green response GΩ of the effective resonant operator (Section 4), and the resulting Φ res feeds back on the baryons through the same phase–lagged gravitational interaction. No inverse operator acts directly on the mass distribution: feedback proceeds solely through delayed curvature forces generated by Φ res and their influence on baryonic motion. This wave–mediated, causal exchange constitutes Channel 1, the resonant feedback loop of the system. In an ensemble of baryonic oscillators (i.e. baryonic elements, each acting as a weak kinematic oscillator through local density–sloshing) the coupling is collective. 2 Individual elements contribute coherently to the resonant curvature 2 The onset of coherence does not require a finely tuned coincidence of phases. As in all weakly damped oscillator ensembles, each baryonic element carries a small stochastic kinematic oscillatory component ( ρb,1 = 0 at arbitrarily low amplitude), and these fluctuations are continually exchanged through the causal Green response GΩ . Any neighbouring pair whose natural frequencies differ by less than the curvature–mediated locking bandwidth ( |ωi−ωj|≲ Γ) 32
field, so that the resulting amplitude scales as Φres ∝pGenv GΩρb,1,(72) where Genv is the ensemble gain that measures statistical phase alignment across the environment. This factor modifies only the strength of collective reinforcement: it counts how many baryonic emitters share a common phase and does not alter causality or the underlying Green function. At the microscopic level this collective response may be described by an effective susceptibility χ (Ω) of the baryonic ensemble, so that the channel behaves as a weakly damped, polarizable medium whose overall gain is controlled by Genv|χ(Ω)|2. In the high– Q regime ( Q = Ω / 2Γ ≫ 1), the memory interval Γ −1 spans many oscillation periods, so baryons and curvature waves remain nearly phase–locked: their relative phase drifts only slowly compared with a single cycle. Repeated exchange between ρb,1 and Φ res then acts as a conservative mode–locking mechanism. Small phase offsets are continuously corrected by the delayed curvature forces, and the ensemble self–organises into a coherent global resonance at frequency Ω, statistically reinforced by Genv . This self–locking is directly analogous to mode locking in optical and acoustic resonators: reversible, phase–lagged exchange sustains coherence, while ensemble gain sets the effective field strength. The reversible exchange between ρb,1 and Φ res gradually loads energy into the resonant channel, which accumulates over long times as the effective stored component ρeff . Because ρeff evolves on the memory interval τmem ∼ Γ −1 , it acts as a slowly varying background contribution that modulates the effective resonant parameters B = ( ceff, Ω , Γ) entering the Channel 1 dynamics. Thus the coupled curvature–matter system continually re–tunes its own resonant environment as energy circulates among ρb,1↔Φres →ρeff,(73) with ρeff feeding back into the quasi–static Channel 2 potential Φ N via Poisson’s equation. Temporal coherence therefore emerges from the slow, collective evolution of mutual susceptibility within this coupled system. As long as the coherence time τcoh =Q Ω,(74) exceeds the characteristic timescale of baryonic perturbations, the resonant channel remains phase–locked and globally coherent. Under these conditions the effective curvature field behaves as a self–sustaining resonator: its phase is stabilised collectively via √Genv and the effective susceptibility, while its spatial envelope evolves on the slower timescales associated with ρeff and the Newtonian background. undergoes slow phase pulling and converges toward a common phase. Once a small coherence patch forms, its enhanced coherent emission strengthens the local forcing and recruits adjacent elements, leading to progressive domain coalescence. This nucleation-and-merging process is directly analogous to the formation of coherent modes in optical and acoustic resonators and requires no special initial alignment of phases. 33
6.3.2 Global Frequency Selection in Inhomogeneous Baryonic Systems The Spacetime Resonance Theorem establishes that any finite coherence region behaves as a weakly damped, self–gravitating resonant enclosure with a collective natural frequency Ω. The Corollary in Section 2 motivates (by analogy with all known weak–field gravitational oscillators) that Ω 2 should scale with an effective coarse–grained baryonic density inside the coherence volume. This density–controlled scaling is not derived from first principles here; rather, it is a physically informed supposition guided by the universal behaviour of familiar gravitational modes (stellar pulsations, Jeans modes, disc breathing modes). We now examine how such a density–dependent frequency behaves when the baryonic medium is spatially inhomogeneous. Consider an extended region in which the coarse–grained baryonic density varies smoothly with position. Each local coherence patch then possesses a slightly different natural frequency, inherited from the local mean density. Once curvature– mediated exchange is active, these patches no longer oscillate independently: the retarded curvature field sourced by each patch acts back on its neighbours with a finite phase lag, producing mutual frequency pulling. Two adjacent patches whose natural frequencies lie within one another’s locking range will synchronise to a common intermediate value, thereby forming a larger coherent enclosure. If the density variation is mild across the system, this merger process repeats iteratively. Curvature exchange continually transports phase–coherent perturbations between neighbouring regions, increasing the coherence volume and pulling the collective frequency toward the density–weighted mean of the expanding domain. Provided that the curvature–mediated coupling bandwidth ( ∼ Γ) exceeds the small natural detunings between adjacent patches, this inductive coalescence proceeds over the long memory time ( Q≫ 1), ultimately producing a single phase–locked resonant enclosure that spans the full inhomogeneous region and oscillates at one emergent global frequency Ω. This behaviour is the precise gravitational analogue of frequency pulling and domain coalescence in inhomogeneously broadened optical resonators (Siegman [ 10 ]). In such systems, local oscillators with slightly different natural frequencies do not compete; delayed coupling through the common field drives them toward a shared frequency and enforces a common phase. In the present gravitational setting, curvature exchange plays the same role: it continually redistributes phase–coherent perturbations among neighbouring baryonic regions and selects a unique, globally coherent oscillation frequency for the entire inhomogeneous structure. 6.3.3 Polarisation of Temporal Coherence The temporal coherence of the resonant channel is not only a matter of frequency locking; it also possesses a polarisation structure determined by the gravitational susceptibility of the baryonic ensemble. In the general case the linear response 34
of a mass element to an oscillatory curvature perturbation is described by a susceptibility tensor χij(Ω),(75) which maps an incident curvature perturbation Aj into a delayed re–emitted field proportional to χijAj . Here χij (Ω) is the linear response of the local baryonic density field to an oscillatory curvature perturbation, not an intrinsic property of individual stars. The effective gravitational polarisation states are given by the eigenvectors of χij(Ω), obtained by diagonalising χijvj (α)=χ(α)vi (α),(76) so that each eigenmode vi (α) behaves as an independent, polarised oscillator with scalar susceptibility χ(α). In a statistically isotropic medium (e.g. pressure–supported ellipticals) the tensor reduces to a multiple of the identity, χij(Ω) = χ(Ω) δij,(77) so all polarisations respond equally; no preferred temporal coherence direction exists. This isotropy leads naturally to three–dimensional coherence domains and the spherical transport envelopes characteristic of ellipticals. In rotating discs, by contrast, the baryonic mass distribution and its dynamical response are strongly anisotropic. Vertical excursions are weak, thin, and largely incoherent, whereas in–plane motions involve substantial coherent mass with long memory. Accordingly the susceptibility tensor decomposes into two unequal sectors, χij = diagχ∥, χ∥, χ⊥,|χ∥|≫|χ⊥|,(78) so the curvature–matter feedback channel amplifies perturbations in the galactic plane far more efficiently than those perpendicular to it. The eigenmodes of Eq. (76) therefore separate into: •adominant in–plane polarisation with large eigenvalue χ∥, and •aweak vertical polarisation with small eigenvalue χ⊥. The resonant feedback loop selects whichever eigenmode yields the greatest coherent reinforcement. Because the environmental gain satisfies Genv,(α)∝ |χ(α)|2, the in–plane mode overwhelmingly dominates in thin discs, forming the planar ( xy ) breathing pattern that drives global m = 1 spiral modes (Section 8.1). Feedback then locks both the baryonic forcing and the effective curvature field into the same polarisation eigenvector: the oscillatory mass distribution drives curvature primarily in the xy plane, and the effective curvature field in turn acts back on the baryons in that same direction, closing the self–consistent loop. 35
In isotropic systems (ellipticals, clusters) no such preferred eigen–direction exists. The susceptibility tensor has three equal eigenvalues, so the resonant mode is unpolarised in the mean: each polarisation contributes equally to the temporal coherence, producing the spatially isotropic speckle cells and large–scale transport envelopes derived in Section 8.2. In mildly anisotropic systems (triaxial ellipticals, S0 galaxies), the eigenvalues of χij split mildly, and the dominant temporal coherence aligns along whichever principal axis maximises |χ(α)|2, giving rise to systematically flattened but non–rotational transport halos. Thus the polarisation of temporal coherence is not an independent ingredient but a direct consequence of the susceptibility tensor of the ensemble. Geometry and kinematics determine χij (Ω); diagonalisation selects the polarisation eigenmodes; and the feedback loop amplifies whichever mode has the largest eigenvalue. Spiral discs, lenticulars, ellipticals, and triaxial systems thereby emerge as distinct polarisation phases of the same curvature–matter resonator. 7 Spatial Coherence of the Resonant Channel The preceding sections established that the resonant sector (Channel 1) acts as a causal, self–locking oscillator whose feedback dynamics drive the system toward a globally synchronized state. Once temporal coherence has been achieved, all baryonic emitters contributing to the oscillatory component ρb,1 acquire a common phase and oscillate at the same resonant frequency Ω = ωres . The purpose of this section is to describe the corresponding spatial coherence of the Channel 1 effective curvature field. Because Channel 1 is defined entirely by the oscillatory source ρb,1 , all coherent structures discussed here belong strictly to the resonant sector. The spatially coherent field Φ res is not directly observable. Its macroscopic imprint enters only through the slowly accumulated, time–averaged stored component ρeff ∝DGΩ∗ρb,1 2Eτmem ,(79) where the average is taken over the long memory time τmem = Γ −1 . Equation (264) expresses the fact that only the stored curvature energy feeds Channel 2; the coherent spatial patterns within Channel 1 remain internally dynamical and observationally silent. The spatial coherence kernel GΩ (x , x ′ ) implicitly defines a spatial polarizability tensor for the resonant channel: it selects which directions of the baryonic fluctuations couple most efficiently into the coherent field. Disk–like coherence domains amplify lateral components of the response, whereas spherical or ellipsoidal coherence cells weight the components almost isotropically. This geometric dependence sets the relative spatial response weights that later enter the susceptibility tensor χij (Ω); its explicit consequences are developed in the spiral and elliptical sections. 36
As proposed in Section 6.3, curvature-mediated frequency pulling drives neighbouring regions toward a common temporal frequency, so a global carrier Ω is already selected. Spatial coherence therefore concerns the geometry of the resonant envelope at this locked frequency, not the mechanism of temporal locking. With temporal coherence established, the resonant field admits the decomposition Φres(x, t) = ReA(x)e−iΩt,(80) where the complex envelope A (x) satisfies a stationary Helmholtz–type equation. Constructive interference among the phase–locked modes produces long–lived, slowly evolving resonant envelopes, the gravitational analogue of standing–wave patterns in optical or acoustic cavities. These envelopes evolve only through slow baryonic advection and the long curvature–memory time of the resonant sector. The next subsection derives the corresponding Helmholtz envelope equation and its exact Kirchhoff–Sommerfeld representation, identifying the spatial coherence scales that govern the internal structure of the resonant field. 7.1 Kirchhoff–Sommerfeld Representation Assuming global synchronization at frequency Ω, all curvature modes carry the common temporal factor e−iΩt . Substituting the ansatz (80) into the linearized field equation (43) and neglecting slow temporal modulations gives the stationary Helmholtz envelope equation ∇2A(x) + Ω2+i2ΓΩ c2 eff A(x) = 4πG c2 eff δρb(x),(81) which determines the spatial profile of the phase–locked resonant field. The exact solution of Eq. (81) is the standard Kirchhoff–Sommerfeld representation: A(x) = G c2 eff Zd3x′eikΩ|x−x′| |x−x′|δρb(x′),(82) with complex wavenumber kΩ=Ω+iΓ ceff , λres =2πceff Ω.(83) The imaginary part of kΩ produces a slow spatial attenuation over the damping length Ldamp =ceff Γ=ceffτmem,(84) which measures the spatial extent of curvature memory. For weak damping (Γ≪Ω) and for r≪Ldamp, one may take kΩ≃2π/λres as real. On galactic and cluster scales the coherence length of the resonant curvature field is of order the resonant wavelength, ℓcoh ∼λres . In spiral galaxies the baryonic disc occupies only a fraction of this coherent envelope, typically adisk ∼ 0 . 1– 0 . 2 λres , but the entire region still belongs to a single phase–locked resonant 37
domain. Within such a domain the curvature phase varies only slowly, and the envelope A (x) behaves as a quasi–static, phase–coherent standing pattern maintained by curvature memory and delayed rescattering, rather than by geometric boundaries as in optical cavities. This definition of the coherence domain is fully consistent with the polarisation–based coherence cells introduced earlier: its size is determined by the curvature wavelength λres , not by the geometric extent of the baryonic material contained within it. In this representation the source term δρb (x) in Eq. (81) should be understood as the coherent Channel 1 driver: the oscillatory component of the baryonic density ρb,1expressed in the global resonant phase frame. Formally, one may write δρb(x) = δρb(x)eiφd(x),(85) where φd (x) encodes the local phase delay of the baryonic motion relative to the curvature carrier e−iΩt . The complex amplitude in Eq. (81) therefore represents the phase–weighted mass density that actually couples to the resonant mode, not a static geometric overdensity. The structure above is directly familiar from scalar Fourier optics 3 (see Goodman [ 11 ]). In a coherent wave system the physical source is represented by acomplex distribution, and the field generated at any point is obtained by a phase–weighted superposition of all emitter contributions. In this language the oscillatory baryonic source δρb (x) = |δρb|eiφd(x) plays the role of a complex aperture distribution, while the curvature kernel GΩ supplies the propagating phase front. The field envelope A (x) is therefore obtained by the standard coherent integral defined in Eq. (82) , which is exactly the Kirchhoff–Sommerfeld representation of the Helmholtz solution. The complex phase of δρb ensures that each baryonic emitter is automatically “rephased” by the propagation kernel before contributing to the field; the integral is simply the coherent superposition of these phase–weighted contributions. This form makes clear that the effective drive of the resonant channel is not determined by the geometric asymmetry of the mass distribution but by its coherent projection onto the curvature phase front. In general, the Green function GΩ(x,x′) = eikΩ|x−x′| |x−x′|(86) acts as the continuous phase reference against which the oscillatory mass distribution is compared. The integral A(x) = G c2 eff Zd3x′GΩ(x,x′)δρb(x′) (87) 3 Unlike standard free–space optical propagation—which lacks an internal distributed memory—this curvature–mediated regime behaves as a quasi–static, phase–coherent resonator whose geometry and long memory determine the envelope structure. Section 7.2 shows how these coherent envelopes generate stored curvature energy in the Poisson channel. Closer physical analogues exist in acoustical media with delayed multiple scattering; see Section 4.2. 38
can thus be interpreted as the phase–matched projection of the complex source δρb onto the resonant effective curvature field. In regions where the baryonic phases φd (x ′ ) are aligned with the kernel, contributions add constructively; where they decorrelate, they cancel statistically. A nearly symmetric configuration may therefore drive the resonant field strongly if its internal phase pattern is coherent, whereas a highly asymmetric distribution may drive it weakly if its phases are random. This viewpoint applies equally to ordered and disordered systems. In a spiral disc, the phase pattern φd (x) may align with a low–order azimuthal structure (treated later via explicit Fourier moments), while in an elliptical galaxy the same formalism describes a random ensemble of local “speckle” phasors whose net drive is governed by their phase–weighted average under Eq. (87) . In all cases it is the coherently projected mass, rather than the raw geometric density, that controls the strength of the Channel 1 field and hence the rate at which stored curvature energy ρeff builds up in Section 7.2. 7.2 Effective Potential and Energy Storage in the Resonant Medium Within the two–channel framework developed above, the physically observable gravitational field is the Channel 2 potential Φ N . This potential is sourced by the time–averaged baryonic density ρb,0 together with the long–term stored curvature energy ρeff, ∇2ΦN= 4πGρb,0+ρeff,(88) and represents the quasi–static geodesic structure of the system. Crucially, Φ N does not contain the instantaneous oscillatory contribution from the resonant channel. The oscillatory effective curvature field Φ res ( t ), driven solely by the coherent baryonic component ρb,1 , averages to zero over a period and is not directly observable. Only its stored curvature energy, accumulated over the memory time τmem = Γ−1, feeds into Channel 2. Because the effective resonant medium possesses finite causal memory, its response to the oscillatory field is phase–lagged. This delay prevents perfect cancellation of curvature stresses over each cycle: part of the resonant curvature energy injected by ρb,1 remains stored in the medium rather than being returned immediately. The result is a slowly evolving curvature bias—the quantity ρeff —that persists on timescales long compared to 2 π/ Ω and sources the observable Channel 2 potential. The local energy density of the Channel 1 resonant field, expressed through the complex spatial envelope A(x) introduced in Eq. (80), is EΦ(x) = 1 8πGhc2 eff|∇A|2+ Ω2|A|2i,(89) representing the reversible partition between curvature–strain gradients and temporal restoring energy. By the equivalence of mass and energy, this stored 39
curvature energy behaves gravitationally as an effective mass density, ρeff(x) = EΦ(x) c2=1 8πGc2hc2 eff|∇A|2+ Ω2|A|2i.(90) This identification ensures that all curvature energy retained in the effective resonant medium enters Channel 2 in precisely the same manner as conventional matter. For statistically incoherent systems—such as elliptical galaxies or cluster environments composed of many locally uncorrelated oscillators—the relevant observable is the ensemble average of Eq. (90), ρeff(x)=1 8πGc2hc2 eff|∇A|2+ Ω2|A|2i,(91) since both intensity and gradient fluctuations contribute to the macroscopic curvature bias. In either the coherent or statistical limit, the stored curvature energy produces a gravitational field identical in form to that attributed in the ΛCDM framework to a dark–matter component, ∇2ΦN= 4πGρb,0+ρeff⇐⇒ ∇2ΦΛCDM = 4πGρb,0+ρDM.(92) Thus ρeff acts observationally as the “missing mass” required for flat rotation curves and weak–lensing signatures, even though its origin lies entirely in the resonant, energy–storing dynamics of Channel 1 rather than in particulate dark matter. The next section examines how this effective resonant medium responds to a single baryonic source, whose localized oscillation constitutes the fundamental building block of all collective resonant structures. 7.3 Single Point Source in a Resonant Channel–1 Background As in the preceding subsections, all results here apply exclusively to Channel 1. The fields derived are purely resonant, oscillatory curvature responses driven by the coherent baryonic component ρb,1 and are therefore not directly observable. Their only physical imprint arises through the stored curvature energy they deposit into the effective resonant medium, which increments the long–term effective density ρeff that later sources Channel 2. For a resonant point mass mb the Channel–1 driver is represented as a Dirac–delta oscillatory source [11]: δρb,1(x)=mbδ(3) x−xp.(93) Inserting this source into the Kirchhoff–Sommerfeld representation (Eq. (82) ) yields the corresponding Channel–1 resonant field δΦp(r): δΦp(x) = Gmb c2 eff eikΩr r, r =|x−xp|.(94) 40
Equation (94) is simply the Green–function evaluation of the Channel–1 Helmholtz operator with a delta–function source. The real, physical Channel–1 oscillation corresponds to the real part of Eq. (94) , giving δΦp(r) = Gmb c2 eff e−r/Ldamp j02πr λres ,(95) where j0 describes spherical spreading and the exponential factor represents attenuation over the curvature–memory length Ldamp . This oscillation remains confined to Channel 1 and has no direct dynamical effect on matter. The quantity of physical interest is the stored curvature energy associated with this oscillation, obtained from the RMS magnitude |δΦp(r)|21/2=Gmb c2 eff e−r/Ldamp r√2j0(2πr/λres).(96) This RMS field is what enters the resonant energy density and ultimately modifies the effective mass density ρeff . In the general case the effective density contains both gradient and curvature contributions, as in Eq. (91) . For the single, spherically symmetric point–source mode considered here we adopt the long–wavelength, slowly varying–envelope approximation ( |∇A|≪kΩ|A| ), so that the gradient term contributes only a small correction to the total stored energy. Retaining only the dominant curvature term ∝Ω2⟨|A|2⟩yields the illustrative expression below. The contribution of this single oscillatory source to the stored curvature energy is therefore ρeff(r) = Ω2 4πGc2|δΦp(r)|2 =G m2 bΩ2 8π c2c4 eff e−2r/Ldamp r2j2 02πr λres . (97) This stored energy is the only mechanism by which a point source modifies observable dynamics: its contribution adds to the total ρeff appearing in the Channel–2 Poisson equation. The instantaneous oscillatory field never contributes directly to ΦN, but the associated ρeff does. Thus the spherical Channel–1 response to a point mass provides the elementary “building block” from which collective resonant fields and the resulting effective density ρeff of extended systems are constructed. In realistic galactic settings the true ρeff arises from the interference of many such point–source contributions to form the global amplitude envelope A (x). The next subsection therefore considers the interference of two oscillatory point sources as the simplest nontrivial example. 41
introduced in Eq. (80) . In this limit, the governing field equation reduces to the Helmholtz envelope equation (81), whose spatial operator Lres−space[A]≡ ∇2A+k2 ΩA, k2 Ω=Ω2+i2ΓΩ c2 eff ,(115) is the Channel–1 resonant operator. In coherent systems (e.g. spiral discs) the wavenumber kΩ is set by the intrinsic resonant scale of the medium, while in disordered systems (e.g. ellipticals) it may be replaced by an effective value keff obtained through the scattering and transport formalisms of Sections 5.3 and 5.4. In all cases, Lres−space governs the spatial structure of the steady resonant envelope that sources the Channel–2 Newtonian potential through the stored energy density ρeff. A stationary attractor is obtained when the spatial curvature pattern A (x) becomes a fixed point of the feedback loop. Its governing relation may be written in the compact form Lres−space[A] + pGenv A= 4πG ρb,(116) where √Genv represents the scalar spatial environmental gain arising from Foldy– Lax–type multiple scattering, speckle–like phonon transport, and reverberation within the resonant medium. Temporal memory—encoded by the quality factor Q = Ω / 2Γ and discussed in Section 7.5—acts as a separate multiplicative amplitude gain. In practice we often work with the temporally renormalized envelope A→Q A , so that Eq. (116) contains only the spatial amplification √Genv . The normalized eigenfunction A (x) then describes the stable Channel–1 curvature mode sustained by a particular feedback geometry, environmental scattering statistics, and damping state. Physically, these steady–state attractors are self–maintaining resonant states. Curvature energy dissipated by damping is continually replenished through coherent or statistical re–excitation by the baryonic ensemble. In the coherent limit—most clearly realised in spiral disks—the system forms a globally ordered, m = 1 standing–wave–like mode analogous to the Foldy–Lax multiple–scattering regime (Section 5.2), where phase–locked rescattering produces spatially coherent gain. For elliptical galaxies, the kinematic disorder increases, and the global phase order is lost and replaced by a statistical balance of intensity, described by the effective–medium and transport equations of Sections 5.3 and 5.4. In such regimes the effective curvature field no longer forms a standing wave but a dynamically sustained intensity distribution that is steady in its ensemble average. The continually fluctuating speckle field redistributes curvature energy stochastically, maintaining the velocity–dispersion pattern that sustains the mode. Despite these differences, all steady configurations satisfy the same causal energy balance: (i) curvature energy injected through organized baryonic motion must equal that removed by damping, and (ii) the local resonant state B = ( ceff, Ω , Γ) 48
relaxes to a stationary form. Each attractor is therefore simultaneously an eigenmode of Eq. (116) and a fixed point of the causal feedback cycle. Once established, it perpetually re–excites itself through the spacetime memory channel described in Section 7.5, while remaining bounded by its own damping and phase statistics. The sections that follow identify three characteristic attractor classes emerging from this unified feedback framework: (i) the spiral–disc configuration, a globally coherent m = 1 resonant mode sustained by ordered rotational streaming; (ii) the elliptical configuration, a statistical or “random–resonant” state governed by diffusive curvature transport and 4D speckle statistics; and (iii) the cluster–core configuration, a nonlinear extension of the elliptical regime in which overlapping diffusive tails undergo collective curvature amplification, yielding a concentrated yet statistically stable resonant core. 8.1 Spiral Galaxies as Coherent m= 1 (Laguerre–Gaussian) Modes 8.1.1 Observing the Steady–State Solution The flat rotation curves of spiral galaxies provide a natural setting for a long–lived, single–frequency resonance of the Channel 1 effective curvature field. Empirically the rotational profile satisfies v(r)≃v0, ω(r) = v0 r,(117) so that all radii share nearly the same orbital speed while retaining their usual differential angular motion. In the present framework this kinematic structure does not impose any global azimuthal synchrony on the baryons; instead, it provides an ordered, predictable background flow in which the effective curvature field can select and sustain its own coherent spatial mode. The only synchrony required for the resonance is temporal. Because v0≪c , the relativistic factor γ(r) = 1 p1−v2 0/c2,(118) is essentially constant across the disk, giving a uniform proper–time relation dτ =dt γ0 .(119) This near–uniformity equalises the oscillation rate of all emitters in proper time, allowing every baryonic element to drive the medium at a common temporal frequency Ω. It synchronises the drive, not the azimuthal motion. The physical origin of this temporal driving is the familiar epicyclic sloshing of stars in a shearing disk. Each star oscillates radially at the local epicyclic frequency κ ( r ), but only the slow, shear–filtered component of this sloshing remains phase coherent 49
across the flat–velocity annulus. It is this coherent, low–drift component—not the full epicyclic motion itself—that selects the effective driving frequency Ω ∼ω ( r0 ). A detailed discussion of this filtering mechanism is given in the next subsection. Once all emitters oscillate at a common temporal frequency Ω, the curvature medium can support a rotating eigenpattern. The large–scale spiral inherits an m=1 azimuthal phase and rotates at the same frequency: Ωpattern = Ω,(m= 1).(120) This is a property of the curvature mode, not of the baryonic disk: the pattern rotates at Ω because the medium responds at this frequency, while the baryons continue to orbit with their natural rate ω(r) = v0/r. For a flat rotation curve, and for any chosen pattern frequency Ω pattern lying between the inner and outer angular velocities of the disk, there exists a corotation radius r0for which ω(r0) = Ωpattern.(121) The shear in the neighbourhood of r0 is relatively weak on dynamical timescales, so baryons there drift only slowly through the pattern and maintain a strong instantaneous projection onto the m = 1 Fourier harmonic. A broad annulus around r0 therefore contributes coherently to the global mode, while more distant radii contribute with smoothly decreasing weight. Thus a flat rotation curve enhances and stabilises the m = 1 curvature eigenpattern: the ordered, low–shear environment sustains a large region with high m= 1 coherence. Spatial organisation of the mode is maintained by rescattering of the Channel 1 field within the baryonic ensemble (Section 5). Each oscillating mass element emits curvature waves that propagate, scatter, and partially reinforce one another, producing a field whose dominant structure is the single–frequency, m = 1 curvature eigenpattern. The detailed gain mechanism (Foldy–Lax multiple scattering) is discussed later; it suffices here that this rescattering network supplies the spatial reinforcement needed for a persistent, phase–coherent mode. Thus the global spiral morphology is the observable imprint of a spacetime– synchronous eigenstate: a coherent m = 1 mode of the Channel 1 operator Lres−space , temporally aligned by the single proper–time frequency Ω and spatially shaped by rescattering in the baryonic disk. The following theorem establishes the existence and stability of this coherent m= 1 spacetime mode. Theorem (Existence and Stability of a Coherent m = 1 Spacetime Mode). Consider a differentially rotating, finite–thickness disk with an approximately flat rotation curve v ( r ) ≃v0 . Let the curvature response of the surrounding medium be described by the driven, damped Helmholtz envelope equation ∇2+k2 ΩA=4πG c2 eff ρb,(122) 50
where A is the complex Channel–1 curvature–potential amplitude, kΩ = (Ω + iΓ)/ceff with Γ>0, and ceff is the curvature–wave speed. Suppose the baryonic emitters drive the medium at a single temporal frequency Ω. Then Eq. (122) admits nontrivial steady–state solutions whose azimuthal dependence consists of a single Fourier harmonic, A(r, θ, z, t)=ℜne−iΩteiθ ˜ A(r, z)o,(123) for some bounded spatial profile ˜ A ( r, z ). These solutions represent coherent, rotating m= 1 curvature patterns. A self–consistent mode arises when the driving frequency matches the orbital frequency at some radius r0, Ω = ω(r0) = v0 r0 ,(124) so that baryons near r0 experience minimal phase drift relative to the pattern. For flat rotation curves a unique such radius exists for any choice of Ω ∈ (0 , ωmax ), and the weak shear around r0 ensures that a broad annulus contributes coherently to the m= 1 harmonic. Because Γ > 0, the homogeneous part of Eq. (122) is exponentially damped. Thus the driven configuration (123) is a linearly attracting steady state: perturbations decay and the system converges toward the coherent m = 1 pattern. The resulting configuration defines a self–consistent steady state in which the effective dark–mass distribution and the baryonic rotation field remain mutually compatible. Corollary (Resonant Coherence Condition). The corotation relation (124) should be interpreted as a steady–state attractor: for the frequency Ω satisfying Ω = ω ( r0 ), the effective curvature field at each point remains phase coherent with its past, and the driven solution is stable under damping. Consequently the m= 1 Fourier sector dominates the long–time response of the disk. Remark on Vertical structure. For numerical work it is convenient to approximate the vertical dependence of ˜ A ( r, z ) by the lowest vertical eigenfunction Z0 ( z ) of the finite–thickness disk and to assume approximate separability, ˜ A ( r, z ) ≈M1 ( r ) Z0 ( z ). This yields an effective 2D in–plane problem for the radial envelope M1 ( r ); the approximation is useful but not required for the existence of the global, phase–locked m= 1 solution. 8.1.2 Energy Source: Epicyclic Sloshing with Coherent Filtering The temporal driving of Channel 1 ultimately originates in the familiar epicyclic motion of stars in a differentially rotating disk. A star on a nearly circular orbit executes small radial excursions, r(t)=rg+Acosκt +ϕ,(125) 51
where rg is the guiding centre, A the oscillation amplitude, ϕ a phase, and κ ( r ) the epicyclic frequency. For a flat rotation curve the standard kinematic identity gives κ(r) = √2ω(r), ω(r) = v0 r,(126) so the sloshing of the radial coordinate—and therefore of the local surface density—occurs at a frequency higher than the orbital rate by a factor √2. Taken over the full disk, this epicyclic motion produces a broad continuum of sloshing frequencies κ ( r ). However, not all of this power is available to the coherent m = 1 mode. The curvature field is highly selective: only those sloshing components whose phases remain slowly drifting relative to the global pattern can contribute constructively. Over the flat–velocity annulus the orbital frequency ω ( r ) varies only weakly with radius, and shear therefore acts as a sharp temporal filter. The only stars whose radial sloshing remains phase coherent over many dynamical times are those whose epicyclic frequency lies within the linewidth of the resonant response, |κ(r)−Ω|≲Γ,(127) where Γ is the damping rate appearing in the Helmholtz operator. Over the coherence time τcoh ∼ (2Γ) −1 the relative phase drift |κ− Ω |τcoh then remains of order unity or smaller, allowing these stars to contribute to a long–lived coherent drive. This condition selects a narrow sloshing band centred on the pattern frequency Ω ≃ω ( r0 ), where r0 is the corotation radius defined by Eq. (121) . Accordingly, the effective driving frequency is not the full kinematic epicycle κ = √2ω , but the shear–filtered component of the sloshing that survives temporal dephasing. Only this small, shear–protected subset of the stellar population contributes constructively to the coherent first azimuthal moment developed in Section 8.1.3. Stars whose sloshing phases drift rapidly relative to the pattern do not add to the coherent drive, even though they continue to oscillate energetically. Crucially, however, all stars participate in rescattering of curvature waves. Every mass element contributes to the Foldy–Lax environmental amplification described in Section 8.1.6, so the gain scales with the full number of scatterers N⋆ , not merely with the coherent set. In this sense coherence determines the effective source, while the entire baryonic population determines the gain. Epicyclic sloshing therefore supplies the raw oscillatory energy; shear filtering isolates a single temporally coherent frequency Ω; and multiple scattering provides the large–scale spatial reinforcement. Together these mechanisms produce the temporally locked, spatially coherent m = 1 resonant mode that underlies the observed spiral halo. 8.1.3 The First Moment The existence theorem establishes that the global effective curvature field settles into a single, temporally locked m = 1 mode. To determine the strength and radial 52
structure of this mode, we now identify the precise baryonic driver that enters the thin–disk Helmholtz equation of Section 8.1.4. This requires extracting, from the full time–dependent mass distribution, the unique component that survives temporal averaging and couples coherently to the m= 1 eigenpattern. We begin by modelling the baryonic emitters as point masses on circular orbits, each oscillating at the common proper–time frequency Ω, ρb(r, θ, t) = N X n=1 m⋆δ(r−rn)δθ−θn(t)e−i[Ωt+∆n(t)],(128) where the orbital phases evolve according to θn(t) = θn0+ω(rn)t, (129) and differential rotation contributes the drift ∆n(t) = −ω(rn)t. (130) All emitters thus share the same temporal frequency while acquiring radius–dependent phases; perfect temporal locking occurs only at radii satisfying ω(r0) = Ω.(131) To extract the coherent azimuthal structure, expand ρbin a Fourier series, ρb(r, θ, t) = ∞ X m=−∞ ρb,m(r, t)eimθ,(132) with ρb,m(r, t) = 1 2πZ2π 0 ρb(r, θ, t)e−imθ dθ. (133) Using (128) gives ρb,m(r, t) = N X n=1 m⋆δ(r−rn)e−i[Ωt+∆n(t)] e−imθn(t),(134) with total phase Ψm,n(t) = −Ωt−∆n(t)−mθn(t).(135) Substituting the orbital evolution, Ψm,n(t) = −mθn0+t[ω(rn)(1 −m)−Ω] .(136) A mode remains phase–coherent only if ω(rn)(1 −m)−Ω = 0,(137) 53
i.e. ω(rn) = Ω 1−m.(138) Only m= 1 allows a physical solution, ω(r0) = Ω,(139) while all m= 1 dephase under differential rotation. It is convenient to write the source in a form that makes the temporal frequencies explicit, ρb(r, θ, t) = ∞ X m=−∞ ˆρm(r) exp−i[(Ω −ω(r))t−(1 −m)θ],(140) with mode frequency ∆Ωm(r)=Ω−(1 −m)ω(r).(141) Only m = 1 satisfies ∆Ω 1 ( r ) = 0 in a flat–rotation disk, and therefore becomes strictly time–independent, ρb,1(r, θ, t) = ˆρ1(r)eiθ.(142) All other harmonics acquire oscillatory factors e−i∆Ωm(r)t and dephase. In a high– Q medium the effective curvature field responds only to the long–time average, ρb(r, θ, t)T−−−→ T→∞ ˆρ1(r)eiθ.(143) Thus the coherent driver of the effective curvature field is precisely the m = 1 component. On scales large compared with the interstellar spacing, the discrete source is replaced by a smooth surface density Σb(r, θ) with vertical profile Z0(z), ρb(r, θ, z)≃Σb(r, θ)Z0(z).(144) The azimuthal average and first moment are Σ0(r) = 1 2πZ2π 0 Σb(r, θ)dθ, M1(r) =Z2π 0 Σb(r, θ)e−iθ dθ. (145) We parametrize the coherent m= 1 component as Σb,1(r, θ) = ϵ(r) Σ0(r)eiθ,(146) so that M1(r)=2π ϵ(r) Σ0(r), ϵ(r) = M1(r) 2πΣ0(r).(147) This coherent moment M1 ( r ) is the unique long–time driver of the global m = 1 spiral resonance. 54
In Section 8.1.4, the stationary envelope equation is reduced to a 2D Helmholtz problem for the midplane effective curvature field. The source term appearing there is precisely ρb,1 ( r, θ ), or equivalently the radial moment M1 ( r ). The remainder of the analysis therefore proceeds by inserting the coherent first moment into the Helmholtz equation to obtain the radial envelope, the outgoing Green function, and ultimately the full spiral morphology. 8.1.4 The Thin–Disk Model: Solving the 2D Helmholtz Equation Reveals the Halo The coherent first moment M1 ( r ) obtained in the previous subsection is the unique long–time driver of the global m = 1 curvature mode. We now determine the spatial structure of the corresponding field by solving the midplane Helmholtz problem. This constitutes the projection of the full three–dimensional eigenmode established in the Theorem in Section 8.1.1 onto the thin disk. 5 The thin–disk treatment represents the limiting case ℓz/Ldamp → 0, providing an accurate description of the galactic midplane and capturing the dominant lateral structure of the coherent m =1 spacetime mode. The 3D model will be built from this in the section that follows. In this limit ( h≪r ), the stationary envelope equation (81) reduces in the midplane to a two–dimensional Helmholtz form, ∂2 ∂r2+1 r ∂ ∂r +1 r2 ∂2 ∂θ2+k2 ΩA(r, θ) = 4πG c2 eff ρb(r, θ),(148) where ρb(r, θ) denotes the projected baryonic density of the disk. Expanding both the curvature envelope and the baryonic source in azimuthal Fourier modes, A(r, θ) = ∞ X m=−∞ Am(r)eimθ, ρb(r, θ) = ∞ X m=−∞ ρm(r)eimθ,(149) is justified because any square–integrable function on the circle ( θ∈ [0 , 2 π )) admits a complete Fourier representation at each radius r . This decomposition follows from the intrinsic 2π–periodicity in θ. The first azimuthal coefficient, ρ1 ( r ), corresponds to the same phase–weighted driver M1 ( r ) introduced in Eq. (147) . Substituting the expansions (149) into Eq. (148) yields, for each mode m, d2Am dr2+1 r dAm dr +k2 Ω−m2 r2Am=4πG c2 eff ρm(r).(150) The homogeneous solutions of Eq. (150) are the cylindrical Bessel functions Jm ( kΩr ) and Ym ( kΩr ), or equivalently the outgoing Hankel function H(1) m = 5 In this projection the envelope A ( r, θ ) represents the lateral component of the full mode A(r, θ, z)=A(r, θ)Z0(z), where Z0(z) is the lowest vertical eigenfunction of scale height ℓz. 55
Jm + iYm , which represents a radiating, causally damped wave (see Appendix A for definitions and asymptotic forms). The coherent spiral attractor corresponds to the dominant m =1 azimuthal mode. Identifying the stationary envelope as A ( r, θ ) = eiθM1 ( r ) / (Γ + i Ω) in Eq. (148) yields the radial Helmholtz equation for M1(r): d2M1 dr2+1 r dM1 dr +k2 Ω−1 r2M1=4πG(Γ+iΩ) c2 eff ρ1(r).(151) The unique outgoing–wave Green function for the 2D Helmholtz operator, consistent with the 3D normalization (∇2+k2 Ω)G=−4πδ3(r), is G1(r, r′) = iπ J1(kΩr<)H(1) 1(kΩr>). This normalization ensures that the far–field asymptotic behaviour A ( r ) ∝ ( GM1/c2 eff ) eikΩr/r matches the three–dimensional outgoing solution of Eq. (82) . The corresponding Green–function solution of Eq. (151) is A(r, θ) = eiθ 4π2iG c2 eff "J1(kΩr)Z∞ r r′M1(r′)H(1) 1(kΩr′)dr′ +H(1) 1(kΩr)Zr 0 r′M1(r′)J1(kΩr′)dr′#, (152) The first integral represents the interior (source–dominated) region that remains regular at r = 0, while the second describes the exterior response that decays outward and satisfies the outgoing (Sommerfeld) boundary condition. It is convenient to define the complex radial amplitude Ar(r) = 4π2iG c2 eff "J1(kΩr)Z∞ r r′M1(r′)H(1) 1(kΩr′)dr′ +H(1) 1(kΩr)Zr 0 r′M1(r′)J1(kΩr′)dr′#, (153) so that the two–dimensional field can be written compactly as A(r, θ) = eiθ Ar(r).(154) This separation makes explicit that all geometric and resonant information resides in the radial envelope Ar ( r ), while the azimuthal dependence eiθ simply encodes the m =1 rotation. The same definition of Ar ( r ) will be used in Section 8.1.5 to construct the three–dimensional field A ( r, θ, z ) and in Section 8.1.7 to evaluate the stored energy and effective dark–mass distribution. 56
8.1.5 Extending the Thin–Disk to Three Dimensions The two–dimensional Helmholtz formulation derived above captures the dominant radial and azimuthal structure of the resonant halo, but a complete description of the coherent m = 1 mode requires inclusion of the finite vertical extent of the disk. In cylindrical coordinates ( r, θ, z ) the stationary curvature envelope satisfies the full three–dimensional Helmholtz equation, ∂2 ∂r2+1 r ∂ ∂r +1 r2 ∂2 ∂θ2+∂2 ∂z2+k2 ΩA(r, θ, z) = 4πG c2 eff ρb(r, θ, z),(155) where ρb ( r, θ, z ) is the baryonic source confined to the thin midplane and its immediate surroundings. The effective curvature field extends beyond this layer in z , but its vertical profile is limited not by the temporal damping length Ldamp = ceff/ Γ, but by a geometric scale height Lz that reflects the finite thickness of the resonant layer itself. Assuming a separable form A ( r, θ, z ) = Ar ( r, θ ) Z ( z ) leads to the vertical eigenvalue equation d2Z dz2+k2 zZ= 0,(156) in which the vertical wavenumber kz is determined from the local dispersion relation k2 Ω=k2 r+m2 r2+k2 z,so that k2 z=k2 Ω−k2 r−m2 r2.(157) Here kr ( r ) denotes the local radial wavenumber of the midplane solution, which may be estimated in WKB form as k2 r≃ |∂rln Ar|2 . Both kr and kz vary only slowly with r , so the product ansatz A ( r, θ, z ) ≃Ar ( r, θ ) Z0 ( z ; r ) remains accurate in the adiabatic (thin–disk) limit. The lowest–order, even vertical solution consistent with midplane symmetry and exponential decay at large |z|is Z0(z;r) = e−|z|/Lz(r) p2Lz(r),(158) where the vertical scale height follows from the evanescent branch of Eq. (157) as L−1 z(r) = pℜ[−k2 z(r)].(159) The normalization in Eq. (158) ensures R∞ −∞|Z0|2dz = 1, so that R|A|2dV is identical to the midplane two–dimensional integral of |Ar|2 . Typically Lz≫ℓz but Lz≪Ldamp , so the effective curvature field remains vertically stratified rather than filling the full damping length. The same Green–function normalization used in the thin–disk solution - see Eq. (152) - is preserved here, ensuring that the three–dimensional extension retains the correct far–field amplitude and outgoing–wave behavior. 57
morphology is a vertically stratified, wave–supported halo consistent with the expected resonant geometry. Figure 6: Three–dimensional effective density distribution. Normalized slices of ρeff for the (a) horizontal (XY) midplane and (b) vertical (XZ) cross–section, each across a 200 × 200 kpc domain. The vertical structure extends to |z| ∼ Lz≃57.65 kpc, consistent with the resonant wavelength scale λres. Although the separable thin–disk ansatz captures the essential resonant behaviour and global halo morphology, it does not constitute a full numerical solution of the three–dimensional Helmholtz equation. More accurate results could be obtained by solving the complete 3D problem without enforcing separability and using the exact baryonic density distribution as a source term. Such an approach would refine the vertical structure, mode coupling, and amplitude normalization within the stellar disk. 8.1.9 Physical Interpretation The global m = 1 spiral mode arises as the dominant solution of the spatial Helmholtz equation driven by a baryonic disk whose differential rotation imparts a slowly varying azimuthal phase. When the baryonic density is written as a complex, phase–weighted source, its first azimuthal Fourier component ρ1 ( r ) isolates the shear–filtered, temporally coherent part of the stellar sloshing—the only portion of the broad epicyclic spectrum that remains phase–aligned at the global carrier frequency Ω 6 . Because this m =1 projection captures the largest shear–protected fraction of the disk’s organised in–plane motion, it drives a large–scale curvature wave with the same azimuthal structure. The resulting field 6 “Coherent” here refers to temporal coherence at the frequency Ω; the stellar orbits themselves are not spatially phase–locked. 64
is therefore not an imposed perturbation but the natural, symmetry–selected mode excited by a rotating, phase–biased baryonic source. A crucial physical ingredient underlying the emergence of this global pattern is the planar polarisation of the resonant channel. As established in Section 6.3.3, a thin, rotating disk possesses a strongly anisotropic gravitational susceptibility, χzz ≪χxx ≃χyy,(174) because vertical excursions involve little coherent mass whereas in–plane motions involve the full surface density of the disk. The coherence (or “speckle”) domain of a spiral galaxy is therefore intrinsically oblate: it is extended laterally but strongly confined in z , and its spatial polarizability enhances curvature–matter feedback almost exclusively within the galactic plane. Consequently the resonant channel selects the in–plane polarisation eigenvector of χij (Ω), and both the baryonic drivers and the effective curvature field are locked into this ( x, y ) polarisation from the outset. This has two immediate consequences. First, the baryonic forcing that excites the resonant mode is necessarily planar: only the in–plane components of stellar motion couple efficiently into Channel 1, and these components contain the dominant organised kinetic energy of the disk. Second, the curvature response inherits the same polarisation. The global resonant solution of the Helmholtz equation is therefore a planar m = 1 mode supported by the large eigenvalue χ∥ , while the vertical eigenmode associated with χ⊥ remains weak and cannot enter the coherent feedback loop. In this way the geometry of the disk and the anisotropy of its susceptibility jointly enforce the flatness of spiral structure: the resonant channel simply does not amplify out–of–plane perturbations. Once this planar polarisation is established, the first azimuthal Fourier component ρ1 ( r ) of the phase-weighted baryonic density provides the dominant coherent driver. Because ρ1 ( r ) captures the largest organised part of the stellar motion, it naturally excites the m =1 curvature pattern with the same polarisation and the same azimuthal phase. This mode is not chosen arbitrarily: it is the symmetryselected eigen-solution of the anisotropically polarised resonant channel when driven by a rotating, phase-biased, planar mass distribution. One might ask why the steady–state attractor does not drive the stellar disk toward a perfectly synchronous m = 1 azimuthal motion, as would occur in a rigidly rotating system. The reason is dynamical: a rising rotation curve would eliminate the shear needed to propagate non–axisymmetric structure, while a falling (Keplerian) curve introduces excessive shear that destroys any pattern on a timescale much shorter than the oscillation period. Only a flat rotation curve—the characteristic kinematic feature of spiral galaxies—creates an extended near–corotation domain: a radial band over which the orbital frequency varies slowly enough that the epicyclic sloshing of stars remains close to the global driving frequency Ω. Within this band the sloshing phases drift slowly relative to the pattern, surviving the shear filtering described in Section 8.1.2. This 65
broad near–corotation region is the essential dynamical ingredient that allows a low–order azimuthal mode, particularly m =1, to persist over tens of kiloparsecs without being sheared apart. The modest magnitude of the azimuthal projection factor, ϵ≃ 0 . 1–0 . 3, is not an adjustable parameter but a direct consequence of this shear–filtered coherence. Only stars whose epicyclic sloshing frequencies lie within the resonant linewidth— i.e. those satisfying |κ ( r ) − Ω |≲ Γ—maintain temporal phase coherence over an oscillation cycle and therefore contribute to the coherent first moment. Material significantly inside or outside the near–corotation annulus experiences too rapid a phase drift (from differential rotation and the κ/ω gradient) to contribute coherently. As a result, only a fraction of the disk’s organised orbital energy projects into the m = 1 eigenmode, while the remainder contributes only incoherently. Fractions of order 10–30% are typical for global modes in shearing disks, and because the resonant curvature channel has a large quality factor ( Q≫ 1), even this partial coherent projection provides ample power to maintain the global m=1 mode and its associated halo amplitude. Although the baryonic disk selects the m =1 mode coherently through its first azimuthal moment, the subsequent build–up of the curvature field is governed by incoherent Foldy–Lax multiple scattering (see Section 5.2). The crucial point is that the three–dimensional Helmholtz Green function GΩ (x , x ′ ) depends only on the difference θ−θ′ when expressed in cylindrical coordinates. This azimuthal invariance ensures that the scattering operator commutes with ∂/∂θ: Seimθf(r, z)=eimθ S[f(r, z)],(175) for any integer m . Thus the m quantum number is preserved by all scattering events. The global m = 1 pattern is therefore robust: it is selected coherently at the source level (temporally at Ω) but its amplitude is amplified statistically through random rescattering without ever mixing into other azimuthal modes. The geometry is coherent; the gain is incoherent; the polarisation is protected by symmetry. The curvature envelope A ( r, θ, t ) carries a well-defined azimuthal phase eiθ and therefore possesses an associated orbital angular momentum (OAM), in complete analogy with optical vortex modes and OAM-carrying fields used in quantum communications and holographic optical tweezers [ 12 ]. Because the underlying Channel–1 field arises from a Lagrangian density, one expects a corresponding dynamical law to be obtainable by varying the action with respect to the baryonic degrees of freedom. Such a law would describe the torque exerted by the m = 1 curvature pattern on the stellar disk and quantify the transfer of angular momentum between the field and the orbiting baryons, thereby reducing azimuthal phase drift near corotation. Although we do not derive this torque equation here, its existence is physically well motivated: a rotating, phase-coherent field generically entrains material motion toward its own pattern speed. In the present context this implies that the m = 1 curvature mode can dynamically stabilise itself by gradually synchronising the baryons (in the temporal, frequency–locking sense), 66
thereby increasing the effective coherence factor ε ( r ) over secular timescales and promoting long-term persistence of the global spiral resonance. Only four sets of physical parameters enter the global mode: (i) the resonant frequency Ω, (ii) the effective propagation speed ceff , (iii) the single–scatterer coupling parameters ( X, ϵ ), and (iv) the quality factor Q . There are no further normalisation freedoms. The reference radius r0must lie somewhere within the observed flat portion of the rotation curve, but its precise location is not uniquely fixed by theory. In practice we selected a value near the inner edge of the flat profile, r0≃ 8 . 5 kpc , which yields v0≃ 230 km s−1 and therefore fixes the global frequency through Ω = v0/r0 . Although this choice is not a priori unique, it produces a broad near–corotation region and leads to excellent agreement with the observed halo structure. The effective propagation speed ceff is not introduced as a fit parameter but follows directly from the requirement that the resonant wavelength λres = 2 πceff/ Ω match the observed radial extent of spiral–galaxy halos. Given the extremely small global frequency Ω = 8 . 77 × 10 −16 s−1 fixed by the flat rotation curve, a natural wavelength of order λres ∼ 40–60 kpc can be obtained only if ceff ∼ 10 5m s−1 . Larger values would shift the resonant scale to unphysical, megaparsec lengths, while smaller values would confine it to only a few kiloparsecs. Thus the adopted value ceff = 2 × 10 5m s−1 is fixed by the requirement that the global curvature mode possess a wavelength comparable to the observed halo span of Milky– Way–like spirals, rather than by empirical fine-tuning. The geometric projection factor ϵ quantifies the m =1 component of the baryonic disk and is bounded by |ϵ| ≤ 1; a conservative value ϵ≃ 0 . 3 is adopted. Similarly, X≃ 1 corresponds to near-maximal single–scatterer coupling to the global mode. With these quantities fixed, the only remaining adjustable parameter is the quality factor Q = Ω / (2Γ), which controls the amplitude of the resonant mode and the radial extent of the halo through the damping length Ldamp = ceff/ Γ. A value Q≃ 12 . 6 simultaneously reproduces: (i) the observed halo extent ( Ldamp ≃ 192 kpc ), and (ii) the overall resonant amplitude required to obtain the empirical dark–to–baryonic mass ratio of ∼ 20. By contrast, the resonant wavelength ( λres ≃ 46 kpc ) and the vertical confinement scale ( Lz≃ 57 . 6 kpc ) are fixed by the chosen propagation speed ceff and temporal frequency Ω, and therefore do not depend on Q . Given that the correct amplitude and spatial extent emerge from a single tunable parameter, this appears to provide strong internal consistency for the resonant interpretation. It is also useful to emphasise the extreme slowness of the resonant frequency itself. For the adopted parameters Ω = 8 . 77 × 10 −16 s−1 , the corresponding oscillation period TΩ≃2.3×108yr is a substantial fraction of a gigayear. Over the ∼ 1 . 4 × 10 10 yr age of the Universe the global mode executes only O (60) such oscillations. In this regime the quality factor Q≃ 12 . 6 has a transparent interpretation: it measures how many of these ultra–slow cycles contribute coherently to the stored curvature energy, corresponding to a memory interval of roughly one fifth of the cosmic age. Even if the entire stellar disk were to 67
disappear instantaneously, the stored effective density ρeff would decay only at the intrinsic rate 2Γ, with a half-life τhalf ≃ ( ln 2) / (2Γ) ≈ 10 9yr . Thus the dark–mass analogue would remain substantially present long after the removal of its baryonic drivers. The strong in–plane polarisation of the longitudinal m = 1 disk mode implies that its baryonic forcing must be similarly polarised. A natural structural analogue is found in the family of low–order self–gravitating oscillations reviewed in Section 2.3.2 and by Cox [ 3 ]. These modes—including the classical radial and low-degree ( ℓ = 1–2) “breathing” oscillations of extended gravitating systems—possess geometries tightly constrained by the underlying mass distribution (disk-like, spheroidal, or ellipsoidal). Although such stellar and cluster modes are not globally phase coherent, their existence shows that largescale, geometry-aligned gravitational oscillations are dynamically admissible in weak–field gravity. The Spacetime Resonance Theorem provides the rigorous basis for such geometry-aligned collective modes in extended baryonic systems, while numerical realisations of the theorem demonstrate that the resulting resonant curvature field reproduces the observed dark–mass profiles of spiral galaxies. Thus the observed family of geometrically constrained low–order modes provides a physically natural analogue for the baryonic sourcing of the coherent m = 1 curvature resonance in spiral disks. Finally, the resonant curvature mode behaves as a damped oscillator with natural frequency Ω and damping rate Γ, so that microscopic dissipation removes effective gravitating density at the rate Ploss = 2Γ ρeff . In the steady state this loss is exactly replenished by coherent orbital driving from the baryons, and the resulting mode amplitude is fixed by the combined action of four physical factors: the quality factor Q = Ω / (2Γ), the single–scatterer coupling strength |X| , the geometric m =1 projection factor ϵ , and the statistical Foldy–Lax gain √N⋆ from the N⋆ baryonic emitters within a coherence volume. Together these give the scaling Aeff ∝Q|X|pN⋆ϵ, (176) expressing that the halo is a self–maintaining resonant structure whose long-term amplitude is determined by the balance between orbital energy injection and dissipation at rate Γ. 8.2 Ellipticals as Random Resonant Media 8.2.1 Observing the Statistical Steady State Unlike spirals, where shear filtering isolates a narrow coherent driving frequency from the sloshing spectrum, ellipticals possess no such frequency selection. Their epicyclic and orbital motions remain fully phase-mixed, so the curvature field is driven stochastically rather than coherently. Spirals maintain a coherent m =1 curvature mode through organized rotation and phase–locked rescattering; ellipticals, by contrast, lack any global azimuthal order. The random orbital 68
motion of stars continually scrambles relative phases, so that coherence survives only locally while the ensemble as a whole behaves as a diffusive resonant medium. Curvature energy is still stored and re–emitted through the same causal mechanism, but now in a statistically fluctuating rather than phase–locked form. This transition from coherent to diffusive resonance is the same one described in radiative–transport theory for acoustic and elastic wave fields [ 13 , 14 , 9 , 4 , 5 ], where random rescattering replaces global phase locking with ensemble–averaged diffusion. In this limit, the coherent multiple–scattering picture described by the Foldy Lax hierarchy (Section 5.2) must be replaced by its statistical extension—the effective–medium and transport framework developed in Sections 5.3 and 5.4. Each star acts simultaneously as source and scatterer of curvature waves, but its motion introduces stochastic phase drift on the scale of the resonant period. The resulting field is a dense superposition of delayed, phase–drifting phasors whose interference forms a curvature speckle field: a dynamically evolving, granular intensity pattern produced by the mutual scattering of moving masses. This pattern is statistically analogous to optical speckle [ 15 ]. In spiral galaxies the response is dominated by a single, globally coherent m =1 mode with a well–defined complex amplitude, whereas in the elliptical regime no single standing–wave mode dominates: the dynamically relevant quantity is the ensemble–averaged intensity built from the incoherent superposition of many phasors. The continual overlap and randomization of these phasors drive a diffusive redistribution of curvature energy throughout the ensemble, as described by the radiative–transfer and diffusion relations of Section 5.4. The corresponding ensemble parameters (Γ eff, Ldamp,eff, Qeff ) govern the macroscopic propagation and memory of the effective curvature field. Because curvature energy is repeatedly rescattered and reused, the effective damping rate Γ eff becomes much smaller than the microscopic Γ, and both the memory time and damping length ( Qeff, Ldamp,eff ) are strongly enhanced. These extended coherence measures explain the smooth, long–tailed 1/r2intensity profiles characteristic of elliptical and cluster halos. Elliptical galaxies therefore act as statistical curvature resonators: locally coherent, globally diffusive systems in which curvature energy is perpetually recycled through the random rescattering of the very masses that sustain it. Their apparent “dark” halos arise from the ensemble–averaged intensity of this self–generated curvature speckle field—a dynamically sustained, statistical extension of the same causal resonance that governs spirals. Theorem (Statistical Steady–State Theorem for Ellipticals). Consider a pressure–supported stellar system lacking global rotational order, in which the resonant effective curvature field is described by a complex envelope A (x , t )and an associated effective gravitating density ρeff (x , t ). Let ukin (x , t )denote the stellar kinetic–energy density. Then the curvature–feedback loop reaches a statistical steady state when the ensemble–averaged effective density 69
obeys d dt ρeff=G−L= 0,(177) where G and L denote the local generation and loss rates of effective gravitating density. These satisfy G ∝ ukin,L= 2Γeff ρeff ,(178) with Γeff the transport (multiple–scattering) damping rate. In the source–free exterior region, where the effective curvature field becomes statistically isotropic and diffusive, the steady–state solution implies an ensemble– averaged density profile ρeff(r)∝exp[−2r/Ldamp,eff ] r2,(179) with transport damping length Ldamp,eff =ceff Γeff ,(180) where ceff is the effective curvature–wave speed. Interpretation. Equation (177) expresses the condition for statistical stationarity: the effective gravitating density generated by random stellar motion balances the density removed by transport damping. Equations (179) – (180) then give the resulting spatial form—a transport–regulated curvature halo with an r−2 intensity tail and an exponential cutoff set by Ldamp,eff. 8.2.2 Complex Baryonic Forcing with Memory A key result of this book is that the curvature–envelope field A (x , t ) is the phononic excitation of the resonant effective curvature field of spacetime, sourced and scattered by baryons. 7 Throughout what follows we therefore interpret A as a continuous field of curvature phasors whose interference carries stress and energy through the system. In regions where these phasors reinforce coherently, the field forms localized packets of compression and rarefaction—phonon-like curvature wavepackets which we refer to as massons (Section 4.3). These are not material quanta but localized solutions of the same linear, causal PDE derived earlier. The governing dynamics are elastic, obeying the same forced–damped Helmholtz equation and transport identities as in acoustic multiple scattering. One may picture a vast ensemble of randomly distributed resonant scatterers, each acting 7 Recall that A (x , t ) is the slowly varying complex amplitude (envelope) of the resonant effective curvature field, Φ res = A e−iΩt + c.c. . It represents a continuous field of curvature phasors arising from the same forced–damped oscillator equation that governs single-mass susceptibility χ(Ω). 70
as a kinematic source while simultaneously obeying a common forced–damped oscillator response that defines χ (Ω), and therefore both emitting and resonantly scattering curvature phonons. This is directly analogous to the resonant–rod experiments of Derode et al. [ 4 , 5 ], in which multiple scattering of diffuse ultrasound produces coherent amplification and long–range interference within a random elastic network [16, 17, 18]. 8 In such a scattering environment, the finite transport mean free path reduces the spatial coherence of the resonant effective curvature field while simultaneously extending its effective damping length. These two effects (shorter coherence, longer diffusive decay) govern how curvature energy is redistributed through the medium and ultimately define the smooth, anisotropic envelope of the effective density field. As shown below, this transition from coherent to diffusive propagation is precisely what gives rise to the extended, elliptical morphology of the dark–mass halo. In spiral galaxies, the low–entropy m =1 orbital order enforced a nearly uniform geometric phase. The complex source term mneiϕn was therefore a coherent phasor sum: the driving was essentially instantaneous and spatially phase–locked (see Section 8.1). By contrast, elliptical galaxies host a phase–mixed stellar population. The geometric phase ϕn evolves stochastically due to random stellar velocities, destroying long–range coherence. Coherence survives only locally in space (transverse scale ℓ⊥) and in time (memory scale τcoh). In an elliptical galaxy the baryonic forcing term retains the full, phase–mixed frequency content of the underlying orbital and sloshing motions, since no shear filtering exists to isolate a narrow coherent band. The appropriate source of curvature phasors is therefore the stochastic phasor representation of the baryonic density, ρb(x, t) = N⋆ X n=1 m⋆eiϕn(t)δ(3) x−xn(t),(181) which drives the same forced–damped Helmholtz operator, ∇2+k2 ΩA(x, t) = 4πG c2 eff ρb(x, t),(182) where m⋆ is the stellar mass, x n ( t ) its instantaneous position, and ϕn ( t ) the geometric phase accumulated along its random trajectory. The underlying PDE is identical to that in the coherent (spiral) case; only the phase statistics differ: in spirals the phases cluster near a common driving frequency Ω, whereas in ellipticals they drift stochastically and provide a broadband, incoherent forcing spectrum. 8 In the acoustic case the rods are driven by an external source; in the gravitational system each mass element is a self–driven kinematic emitter, while the same causal forced–damped PDE determines how its emission propagates, scatters, and damps. 71
The field amplitude A (x , t ) may be written using the causal (single–scatterer) Green function: A(x, t) = 4πG c2 eff N⋆⋆ X n=1 m⋆eiϕn(t)GΩ x−xn(t),(183) The causal Green function GΩ (r) describes the response to a single resonantly driven stellar source. Its amplitude decays in the far field as |GΩ(r)|2=e−2|r|/Ldamp (4π)2|r|2, Ldamp =ceff Γ,(184) representing the causal damping of an individual curvature phasor. Ensemble averaging over many such phasors will later produce the transport–renormalized kernel. In Eq. (183), N⋆⋆ denotes the full population of past emitters within the causal memory window t−τmem < t′≤t, where τmem = Γ −1 is the memory time of the underlying forced–damped oscillator response and τcoh = (2Γ) −1 its shorter phase–decorrelation time. Equation (183) therefore represents a four–dimensional ensemble of phasors, each carrying its causal phase drift across the memory window. Not all of the N⋆⋆ phasors are statistically independent: correlations persist over finite regions of space and time, defining a four–dimensional speckle scale determined by the coherence of the same underlying PDE. The following subsection therefore defines the characteristic size of a four–dimensional speckle cell, within which phasors remain phase–correlated. Determining this coherence volume and timescale allows the number of statistically independent phasors within N⋆⋆ —and the corresponding baryonic mass— to be quantified. These quantities are then used in Section 8.2.4 to rewrite Eq. (183) in terms of the independent phasor ensemble and to compute ensemble averages such as ⟨|A|2⟩ and ⟨ρeff⟩. 8.2.3 Four–dimensional speckle geometry of the gravitational field To extract statistical predictions from Eq. (183) , we must determine over what regions of space and time the gravitational field A (x , t ) retains phase coherence. Because elliptical galaxies lack global orbital order, coherence is neither disk–wide nor steady in time: it is localized and transient. The field therefore forms a granular, four–dimensional speckle pattern—the gravitational analogue of diffuse ultrasound speckle. Transverse coherence. Random stellar motion produces geometric phase drift. When the phase difference between two emitters reaches one radian, their contributions become statistically uncorrelated. This defines the transverse coherence length ℓ⊥=λres 2π,(185) 72
so phasors launched from within a transverse patch of radius ℓ⊥ interfere coherently. Longitudinal coherence. Along a star’s trajectory, orbital motion also scrambles phase in time. During the coherence interval τcoh the typical orbital excursion is ℓ∥=v τcoh =v 2Γ,(186) so emissions separated by more than ℓ∥along an orbit are decorrelated. Space–time coherence cell. The smallest region over which phasors retain a common phase is a transverse disk of radius ℓ⊥ extended over a longitudinal distance ℓ∥. Its spatial volume is Vcell =α ℓ2 ⊥ℓ∥=αλres 2π2v 2Γ,(187) where α is a geometric factor of order unity. The corresponding temporal extent is ∆tcell =τcoh =1 2Γ.(188) This pair ( Vcell, ∆ tcell ) is the gravitational analogue of the isoplanatic patch and coherence time of diffuse acoustics. Population statistics. Let Vgal denote the active baryonic volume. The number of disjoint coherence cells at any moment is N3D =Vgal Vcell ,(189) and each cell contains on average Ncell =n⋆Vcell =N⋆Vcell Vgal ,(190) stars, where n⋆is the mean stellar number density. Effective phasor mass. Within a single cell all stars share, to leading order, a common field phase. Their masses therefore add coherently as one composite phasor of effective mass Mcell =Ncell m⋆.(191) The original stellar sum in Eq. (183) may thus be re-expressed as a sum over statistically independent four–dimensional coherence cells, each weighted by its complex mass Mcell and the appropriate Green kernel. In summary, coherence in ellipticals is local in space and finite in time. The gravitational field is driven not by a global mode but by a four–dimensional mosaic of statistically independent phasor ensembles. Rewriting Eq. (183) in this statistical basis enables the evaluation of ensemble quantities such as ⟨|A|2⟩ and the mean effective density ⟨ρeff ⟩, as developed in the following subsection. 73
Defining a correlation length ℓcby ℓ−2 c:= |∇A|2 |A|2,(217) Averaging Eq. (216) and using the definition ℓ−2 c=⟨|∇A|2⟩/⟨|A|2⟩gives ρeff(r)=Ω2 8πGc2D|A(r)|2E1 + c2 eff Ω2ℓ2 c≡Ξ∇ Ω2 8πGc2D|A(r)|2E,(218) where the dimensionless gradient–energy boost is Ξ∇= 1 + 2π Ldamp,eff λres 2 ,(219) obtained by setting ℓc≃Ldamp,eff and Ω = 2 πceff/λres . This gradient–energy boost enhances the total curvature energy stored in the coherent mode without changing the spatial shape of the halo. In the ballistic limit, Ξ ∇∼O (1); in the transport limit, it can become very large when Ldamp,eff ≫λres. Multiplying Eq. (210) by the gradient–energy boost Ξ ∇ and applying the small– r regularization from Eq. (215) yields the final, fully normalized transport halo law: ρeff(r)=4πG c2 eff 2Qeff 2πΞ∇Genv,eff r, r < Rmatch, e−2r/Ldamp,eff (4π)2r2, r ≥Rmatch. (220) For convenience, we collect the transport–renormalized parameters: Genv,eff =|X|2N2 ⋆ Vgal αλres 2π2v 2Γeff , Qeff =Ω 2Γeff ,Ξ∇= 1 +2πLdamp,eff λres 2 . (221) Together, these factors determine the amplitude of the regularized curvature–phonon halo: Genv,eff encodes the coherent mass content of the stellar ensemble, Qeff sets its temporal memory, and Ξ ∇ provides the gradient–energy boost associated with transport–extended correlations. The resulting law defines a smooth, finite–core, exponentially damped envelope that links microscopic scattering dynamics to the macroscopic curvature structure of elliptical halos. 8.2.8 Numerical Results The boosted transport–halo law (Eq. 220) predicts a self-consistent, two–scale halo structure for elliptical galaxies: a small-radius linear core (Section 8.2.6) and a 80
large-radius transport tail ∝exp ( − 2 r/Ldamp,eff ) /r2 , with the amplitude controlled solely by the resonant wavelength λres . The numerical results presented here are obtained using the reference model defined in Appendix C, with Ldamp,eff = 636 kpc tuned to match the observed halo extent of massive ellipticals, and λres = 27 pc chosen to reproduce a realistic dark–to–baryon mass ratio. Figure 7 shows radial diagnostics for the resulting ensemble–averaged halo. Panel (a) reveals that the effective dark density becomes comparable to the stellar density at ∼ 5 Re ( ∼ 25 kpc), rises above it, and gradually dominates throughout the transport tail. Panel (b) quantifies this with the radial density ratio ρeff/ρb , which reaches a maximum near r∼ 700 kpc. Panel (c) shows that the dark–induced rotation curve rises rapidly from the core, develops a modest ∼ 200 km,s −1 shoulder extending over several tens of kiloparsecs, and then declines slowly. This behavior is consistent with that observed in luminous ellipticals [ 20 ]. Panel (d) displays the enclosed mass ratio Meff/Mb : unity near 27 kpc, ∼ 3 by 100 kpc, ∼ 7 by 300 kpc, and asymptoting to ∼ 11 at ∼ 1 . 3 Mpc. These values fall squarely within the observed dark-to-luminous mass fractions for giant ellipticals in deep X-ray and weak-lensing surveys. Figure 8 visualizes the spatial structure of the halo in the central ± 50 kpc: panel (a) overlays the stellar and effective components using a log-scaled composite (red: baryons; green: curvature–transport dark mass), highlighting the strong baryon dominance near the center and the smoother, more extended spatial footprint of the resonant component. Panel (b) shows their combined log-normalized density distribution. Finally, Fig. 9 shows XY and XZ slices extending to ± 1 Mpc, demonstrating the near-spherical envelope of the coherent transport tail. The horizontal and vertical cuts are visually indistinguishable, confirming the expected isotropy of the far-field ensemble solution. The two free quantities in the model have clear and distinct roles: Ldamp,eff ⇒halo radial extent, λres ⇒total dark mass fraction. Tuning λres and Ldamp,eff therefore provides a transparent calibration of ellipticals along their observed halo size–mass sequence without any hidden or fine-tuned parameters. The results above constitute a first demonstration that resonant curvature–phonon transport alone can reproduce the canonical dark-mass phenomenology of pressure-supported galaxies. 8.2.9 Halo Formation from Coherent Energy Flux The boosted transport halo of Eq. (220) can equivalently be derived from energy–flux conservation in the source–free exterior region. In this picture, the resonant halo represents a stationary outflow of curvature–phonon energy, whose spatial profile follows directly from flux continuity and transport damping. 81
Figure 7: Radial halo diagnostics for the boosted transport–law model. (a) Effective vs baryonic mass densities. (b) Density ratio. (c) Circular velocities showing a modest ∼ 200 km s −1 shoulder. (d) Enclosed mass ratio rising to ∼ 11 at ∼1.3 Mpc. 82
Figure 8: Midplane maps ( ± 50 kpc). (a) Composite stellar (red) and transport–dark (green) mass distributions. (b) Combined log-scaled density. The effective component spreads well beyond the stellar extent even in the central region. Figure 9: Large-scale XY and XZ slices ( ± 1 Mpc) of the resonant dark component, log-scaled and normalized. The near-spherical symmetry of the transport halo is evident. 83
Outside the baryonic domain ( r > Rmatch ), the ensemble-averaged field behaves as an outgoing, spherically symmetric transport mode F(r)=ceff u(r), ceff du dr +2u r=−2 Γeff u(r),(222) where u ( r ) is the curvature–phonon energy density. The unique outgoing solution of Eq. (222) is u(r)∝exp[−2r/Ldamp,eff ] r2, Ldamp,eff =ceff Γeff .(223) The corresponding effective gravitating density thus takes the form ρeff(r)=Kflux exp[−2r/Ldamp,eff ] r2,(r > Rmatch) (224) with a proportionality constant determined by flux matching to the interior ensemble: Kflux = Ξ∇Genv,eff Qeff ,(225) where Ξ ∇ is the gradient–energy boost, Qeff = Ω / (2Γ eff ) the transport quality factor, and Genv,eff the coherent–mass coupling defined in Eq. (208). This construction shows that the halo profile arises directly from energy–flux conservation: the r−2 decline is a geometric consequence of spherical outflow, the exponential attenuation encodes transport damping over Ldamp,eff , and the overall amplitude reflects the coherent energy stored and recycled within the stellar ensemble through Genv,eff,Qeff , and Ξ∇. 8.2.10 Physical Interpretation Elliptical galaxies represent the self–sustaining limit of the transport–resonant curvature–phonon framework developed above. In these pressure–supported systems the random three–dimensional stellar motions both generate and scatter the effective curvature field. Each star acts as a moving phase perturbation of the local spacetime curvature, and the ensemble of such perturbations forms a dense, multiply–scattering phononic medium governed by the Dyson self–energy Σ(Ω). Unlike spirals, the stellar motions do not select any preferred driving frequency: the sloshing is broadband and incoherent, and the single oscillation frequency that appears in the curvature field is imposed by the resonant Helmholtz operator itself, not by stellar phase synchronisation. Kinetic energy injected by the baryons is continually converted into curvature phonons, which scatter back upon the baryons and help maintain the velocity dispersion that excites them. The galaxy as a whole thus behaves as a gigantic, statistically stationary, self–excited transport resonator in which microscopic generation and ensemble damping remain in dynamic equilibrium. In the diffusive limit, the relevant resonant scale is determined by the density and geometry of the coherent domains—or “4D speckle cells”—that compose the 84
stellar ensemble. Each cell of characteristic size λres encloses a coherent stellar mass Mcell within a coarse–graining kernel of volume Vcell ∼λ3 res . Within such a domain, phase correlations cause curvature amplitudes to add linearly while ensemble–averaged intensities add quadratically, yielding ρeff ∝M2 cell Vcell c2,(226) where ρeff is the local ensemble–averaged effective gravitating density. This local coherence is enforced by the long-wavelength Green function and finite memory time of the resonant PDE—not by any dynamical phase locking among the stars themselves. The resulting profile therefore depends jointly on the number of coherent sources within a speckle cell and the compactness of the cell itself. A central prediction of the Spacetime Resonance Theorem (Section 2) is that the natural restoring frequency of the resonant curvature mode within any coherence domain scales with the square root of the coarse–grained baryonic density, Ω ∝√¯ρb , where ¯ρb denotes the mass density averaged over the local coherence kernel. This mirrors the well–established ω∼√Gρ scaling found in small–amplitude oscillations of self–gravitating systems. Therefore, the density of a finite region determines the resonant frequency of the curvature mode, and with ceff fixed, the resonant length scale follows as λres = ceff/ Ω ∝¯ρ−1/2 b . In ellipticals the coherence kernel is effectively isotropic, set by the short mean free path of curvature phonons in a dense, three– dimensional stellar field, yielding small λres . In spirals, coherence extends anisotropically through the disk, giving larger coherence volumes, lower ¯ρb, and much longer wavelengths. With ceff ≈ 2 × 10 5m s−1 fixed, the observed contrast between systems follows directly from geometry and density. Typical coarse–grained densities of ¯ρb∼ 10 −21 –10 −22 kg m −3 in ellipticals yield λres ∼ 20–30 pc. Disk–averaged densities around 10 −24 kg m −3 in spirals give λres ∼ 30–50 kpc. Cluster–scale cavities with ¯ρb≲ 10 −26 kg m −3 support modes approaching megaparsec scales. Thus the short wavelengths of ellipticals reflect their compact, three–dimensional coherence and high baryonic density, whereas spirals and clusters express progressively larger, lower–density resonant domains. The issue of sustainment concerns not the static ensemble profile but the persistence of the microscopic fluctuations that continually regenerate the effective curvature field. In the diffusive–transport regime each scattering event introduces a small delay and phase randomization, preventing the instantaneous field from ever relaxing into a perfectly smooth configuration. The central regions therefore exhibit a shifting speckle pattern of bright and dim curvature zones—localized fluctuations evolving on the coherence timescale τcoh,eff ∼ 1 / (2Γ eff ) (Section 8.2.5). Although minuscule compared with the stellar kinetic energy, these fluctuations form a microscopic phase–randomizing bath that continuously refreshes the stochastic component of the halo. On macroscopic scales this agitation establishes a steady exchange between baryonic kinetic motion and the effective gravitating density carried by curvature 85
modes. In ensemble average the balance takes the form d⟨ρeff⟩ dt =⟨G⟩−⟨L⟩ ≃ 0,(227) where the generation rate is proportional to the local kinetic–energy density, G ∝ ukin , and the loss rate is L = 2Γ effρeff . Stationarity therefore requires that random stellar motion replenish the effective density at precisely the rate at which ensemble damping removes it. This is the same self–consistent energy closure that fixes the transport damping rate Γ eff in the Dyson effective–medium theory. This balance is directly analogous to the classical Jeans equilibrium. In the Jeans picture random stellar velocities generate a pressure gradient balancing self–gravity. Here, the same random motions generate curvature–phonon energy ( G ∝ ukin ), while the resulting effective density feeds back through the transport damping term L = 2Γ effρeff . The steady state is thus the resonant–transport counterpart of Jeans support: kinetic agitation sustains the effective curvature field, and the effective curvature field confines the very motions that regenerate it. Because the speckle bath is continually refreshed, the ensemble–averaged effective density behaves as a causal moving average of the instantaneous field: ⟨ρeff⟩t+∆t= (1 −β)⟨ρeff ⟩t+β ρinst(t), β = 1 −e−∆t/τcoh,eff ,(228) with τcoh,eff determining the memory time of the transport process. In stationary regions the average converges to a constant value; in evolving environments it governs how the ensemble field tracks transient curvature fluctuations. Through this continual renewal on the coherence timescale, the system maintains a persistent background of stochastic curvature stresses that both mirror and sustain the random stellar motions themselves—thus closing the feedback loop between baryonic agitation and curvature diffusion. Although the mean effective curvature field is smooth and energetically subdominant, its finite variance implies a nonzero probability of localized curvature intensification. If such peaks persist for several coherence times they may weakly bias the nearby stellar distribution, producing transient clustering or faint substructure. These features arise naturally from the feedback loop itself: the very fluctuations that preserve equilibrium occasionally imprint brief, localized curvature enhancements within an otherwise relaxed halo. The polarisation properties of the curvature–matter resonance in ellipticals follow from the susceptibility tensor χij (Ω) introduced in Section 6.3.3. Although individual four–dimensional speckle cells are locally anisotropic—each having transverse and longitudinal scales ( ℓ⊥, ℓ∥ ) set by the stellar velocity distribution—the ensemble of cells is randomly oriented and statistically well mixed. Because the stellar velocity ellipsoid in pressure–supported ellipticals is nearly 86
isotropic, this random orientation washes out the anisotropy of individual cells in the ensemble average. The linear response therefore reduces to the isotropic form χij(Ω) = χ(Ω) δij,(229) so all tensor eigenmodes share the same susceptibility, χ(x)=χ(y)=χ(z). No preferred polarisation axis exists in this limit. Three–dimensional stellar agitation excites all eigenvectors vi (α) with comparable strength, and rapid phase mixing removes any transient anisotropy. Ellipticals are therefore statistically unpolarised: their curvature–phonon field populates all tensor eigenmodes equally, consistent with the nearly spherical transport envelopes of pressure–supported systems. This framework naturally extends to lenticular, triaxial, and irregular systems. Whenever the stellar velocity distribution departs from isotropy, the ensemble– averaged susceptibility acquires a small but finite eigenvalue splitting, χxx =χyy =χzz,(230) with the principal axes aligned with those of the stellar velocity ellipsoid. The coherence kernel remains three–dimensional but becomes mildly oblate or triaxial, imparting a correspondingly weak polarisation bias to the global effective curvature field. Such systems lie between the isotropic (giant elliptical) and strongly planar (spiral) limits. In S0 galaxies, the residual disk component enhances the in–plane susceptibilities, χxx ≃χyy > χzz,(231) producing an oblate coherence tensor and a flattened transport envelope, though still without the strong phase locking or large eigenvalue contrast characteristic of spiral galaxies. Triaxial ellipticals exhibit three distinct susceptibilities, each aligned with a principal inertial axis, naturally yielding the observed E3–E7 sequence as a continuous family of anisotropic transport–coherence tensors. Thus the projected halo ellipticity traces the eigenvalue hierarchy of χij(Ω). Irregular galaxies represent the partially coherent limit. Low baryonic densities and turbulent, stochastic environments prevent long–lived phase locking, yet short–range curvature correlations arise within transient stellar associations. Each system is therefore a fluctuating ensemble of small coherence patches whose cross–correlations are weak but nonzero, 0<⟨AiA∗ j⟩ ⟨|A|2⟩<1,(232) continually forming, dissolving, and merging without establishing a global mode. This picture agrees with their patchy morphology, irregular kinematics, and diffuse halos: irregulars are turbulence–driven curvature ensembles that occupy the transition between the coherent and fully random extremes of the resonant– speckle hierarchy. 87
8.3 Clusters as Radiative Transport Extensions of Ellipticals 8.3.1 Observing the Statistical Transport Equilibrium In this section we extend the radiative transport formalism developed for isolated ellipticals and spirals to the cluster scale. Whereas the elliptical case required a transport treatment to capture long-lived diffusive tails, and the spiral case was fully described by coherent Foldy–Lax scattering, the cluster environment involves the superposition of many independent elliptical transport envelopes. These envelopes overlap and partially interfere, producing a collective intensity field that must be modelled statistically. The primary goal of this section is to demonstrate that the resulting effective transport density reproduces, in both form and amplitude, the observed darkmatter distribution of clusters as described by the Navarro–Frenk–White (NFW) profile [21, 22]. The NFW model, derived empirically from cosmological N -body simulations, describes the dark-matter halo as ρNFW(r) = ρs x(1+x)2, x =r rs ,(233) with rs and ρs as characteristic scale radius and density, respectively. Recent extensions include the introduction of a “splashback” radius Rsp , beyond which ρ(r) rapidly steepens [23]. The theoretical basis for the cluster-scale treatment presented in this section, lies in the long-standing radiative transport framework developed for diffuse acoustic and optical propagation. Weaver [ 13 , 14 ] established that acoustic energy in heterogeneous solids obeys a radiative transfer equation for the intensity field, with distinct ballistic, coherent, and diffusive regimes. These works provided the first quantitative formulation of acoustic diffusion and localization, defining the statistical kernel that our own model generalises to astrophysical scales. The subsequent experimental validation was achieved by Derode, Tourin, and Fink [ 9 , 4 , 5 ], who investigated ultrasonic propagation through random ensembles of resonant scatterers. Their experiments confirmed the radiative transport regime predicted by Weaver, demonstrating the separation of coherent, ballistic, and diffuse components, as well as the emergence of reverberation gain in multiply scattered fields. These resonant-rod systems represent direct analogues of our elliptical galaxies: individual mesoscopic resonators embedded in a diffusive host that collectively generate extended energy-density tails. Later developments by Weaver and Lobkis [ 16 , 17 ] established that diffuse fields from independent random ensembles can exhibit weak cross-coherence, recoverable through cross-correlation measurements. This phenomenon—now widely used in seismology and acoustics—formalised the concept of overlapping 88
diffuse transport fields, providing a rigorous physical basis for the small nonlinear coupling parameter ηused in our cluster model. Finally, Larose et al. [ 18 ] unified these concepts across acoustics, optics, and geophysics, showing that the superposition of diffuse fields from independent sources can reproduce effective Green’s functions through reverberant correlations. This provides the closest experimental and theoretical precedent for our treatment of overlapping elliptical transport envelopes within clusters: random, extended resonators (galaxies) embedded in a diffusive medium (intracluster baryons) that generate partially coherent, overlapping transport tails and a measurable ensemble reverberation gain. The ensemble behaviour demonstrated in the acoustic and optical transport experiments of Weaver and Derode provides the direct methodological foundation for the present cluster model. Here we apply the same radiative–transport formalism—previously shown to describe overlapping diffuse fields in random media—to the collective curvature–transport regime of spacetime itself. Each elliptical galaxy sustains its own diffusive envelope, generated by the random stellar motions that define its internal transport equilibrium, and these envelopes coexist within a shared effective curvature field (the common transport continuum of spacetime) that supports their mutual overlap. The resulting ensemble of overlapping transport intensities forms a statistically smooth superposition whose weak cross–coherence produces the modest reverberation gain Grev =1+Ngalη, (234) which, when applied at the cluster scale, reproduces the NFW–like curvature of the effective density profile without invoking any additional physical assumptions. In the sections that follow, we adopt identical methodologies to those used in the acoustic and optical transport studies, applying them directly to the astrophysical transport regime to construct and normalize the cluster-scale effective density field. Each elliptical or lenticular acts as amesoscale driver of the effective curvature field, ρeff,i(r)∝exp(−r/Lmean),(235) where Lmean =ceffτtr (236) is the mean transport distance determined by the effective transport time τtr . Within a cluster, these envelopes overlap and interfere weakly. Their mutual correlation defines a small but finite reverberation parameter η, ⟨AiA∗ j⟩=ηq⟨|Ai|2⟩⟨|Aj|2⟩,(i=j),(237) yielding the statistical amplification Grev =1+Ngalη, (238) 89
Grev = 1 + Ngalη to incorporate partial coherence between overlapping tails. The complete algorithm and parameter definitions are given in Appendix D, together with the Python code used to generate the figures. For a fiducial 10 15 M⊙ cluster, we adopt Ngal = 300 ellipticals, Lmean = 636 kpc , σln L = 2 . 5, and η = 8 . 3 × 10 −3 , corresponding to Grev ≃ 3 . 5. The resulting one-dimensional transport profiles are presented in Figure 10, which compares the model’s effective density ρeff ( r ) and enclosed mass Meff ( r ) with a standard NFW benchmark. All quantities were computed on a logarithmic grid spanning 0 . 2 kpc ≤r≤ 3 Mpc , using the 2.5D convolution for ρshape ( r ) and normalized by Λ geom as defined in Eq. (253) to preserve total enclosed mass. The reference NFW halo corresponds to M200 = 1015 M⊙and c200 = 4. Figure 10: Cluster density and enclosed mass. (a,b) Effective density in linear and log–log scales. (c,d) Corresponding enclosed-mass profiles. The transport–overlap model (solid) reproduces the NFW benchmark (dashed) over three decades in radius, with minor deviations at small r attributable to the missing central BCG component. Panels (a) and (b) of Figure 10 show the effective density in linear and logarithmic scales, respectively. The transport–overlap profile reproduces the NFW slope 96
Figure 11: Circular velocity profiles. The transport–overlap prediction (solid) agrees closely with the NFW velocity curve (dashed) across the cluster, confirming that the cumulative mesoscopic transport reproduces the macroscopic gravitational potential. Figure 12: Two-dimensional distributions of effective density and velocity. (a) Normalized log10 ρeff intensity map. (b) Linear-scale circular-velocity field vc . Both panels exhibit smooth radial symmetry and NFW-like scaling. 97
and amplitude across three orders of magnitude in radius, maintaining close agreement from ∼ 30 kpc to 3 Mpc . The slightly increased concavity of ρeff ( r ) at small radii ( r≲ 50 kpc ) possibly arises from the omission of the central brightest cluster galaxy (BCG), whose extended stellar envelope would naturally flatten the inner slope. The shallow dip near r→ 0 originates from the discrete angular integration and is purely numerical. Panels (c) and (d) display the corresponding enclosed-mass profiles. The total mass derived from the transport model follows the NFW prediction almost exactly over the full domain, confirming that the normalization and reverberation gain preserve global mass balance. Within the range 1–3 Mpc , the overlap of exponential tails produces an effective power-law transition identical to that of the NFW halo, with the chosen σln L = 2 . 5 being essential to the match. Smaller values yield profiles that are too steep, whereas larger ones over-broaden the outer envelope. Figure 11 shows the corresponding circular-velocity curve, which matches the NFW form both in amplitude and in the location of its characteristic peak near r≃rs . The gradual decline beyond the virial boundary reflects the splashback taper Wsp ( r ), which enforces a finite physical extent without artificial cutoffs. Figure 12 presents the projected two-dimensional maps of normalized logarithmic density and linear velocity. The maps exhibit excellent radial symmetry and smooth gradients, confirming that the 2.5D convolution produces a numerically stable and physically consistent transport field even in projection. Across all figures, the agreement between the transport–overlap model and the NFW benchmark is excellent for r≳ 2 × 10 −2R200 . The model reproduces the NFW slope and amplitude throughout the 1–3 Mpc range, demonstrating that the combined exponential tails of individual galaxies collectively mimic the canonical dark-matter profile. Minor residuals at small radii are expected given the omission of the BCG and the discretization of the angular integration. Overall, these results show that the observed dark-matter behaviour of galaxy clusters can be reproduced by the coherent superposition of mesoscopic transport envelopes. Each galaxy contributes an exponentially decaying intensity field whose cumulative overlap, after geometric normalization and modest reverberation gain, yields a cluster-scale density profile that is nearly indistinguishable from the NFW form. The model thus provides a self-consistent bridge between between the mesoscopic transport physics of individual elliptical galaxies and the macroscopic gravitational phenomenology observed in clusters. 8.3.4 Physical Interpretation Galaxy clusters represent the statistical–transport limit of the curvature–phonon hierarchy. Spiral galaxies support coherent resonances, and ellipticals sustain self–contained diffusive envelopes. Clusters extend this hierarchy by combining many such envelopes within a shared curvature–transport continuum. The cluster therefore behaves as a radiative–transport ensemble: a largeN superposition 98
of weakly coherent resonators whose intensities add statistically to produce a smooth macroscopic curvature field. Phase information is effectively erased, but the residual cross–coherence encoded in the reverberation parameter η generates a finite amplification of the ensemble mean. This amplified statistical field reproduces the canonical NFW profile without invoking any additional mass component. Each elliptical galaxy emits a curvature–phonon field whose intensity decays according to the single–galaxy transport kernel Keff ( r ). Because the transport length Lmean far exceeds the internal coherence scale of the galaxy, its diffusive envelope cannot readjust instantaneously as the galaxy orbits within the cluster. This produces a natural curvature wake: a trailing region of stored effective density that lingers for times of order τtr and extends over hundreds of kiloparsecs. The curvature wake is therefore the single–galaxy manifestation of the more global reverberant field sustained by the entire cluster. At the ensemble level, the wakes and transport envelopes of many ellipticals overlap and merge. The resulting statistical superposition forms a smooth field characterised by a small but finite mutual coherence, ⟨AiA∗ j⟩=ηq⟨|Ai|2⟩⟨|Aj|2⟩,(i=j),(256) which produces the reverberation gain Grev =1+Ngal η. (257) This gain is the cluster-scale analogue of the reinforcement observed in acoustic and optical multiple scattering: weak cross–coherence among diffuse fields yields a statistically enhanced intensity. Physically, the cluster curvature field isareverberant transport medium, maintained by continual redistribution and reprocessing of curvature energy among its member galaxies. Energetically, the cluster reaches a statistical transport equilibrium in which curvature energy circulates diffusively through the ensemble. Emission from one galaxy is scattered and re-emitted by its neighbours with randomized phase, establishing a globally stationary mean intensity: d ρeff dt ≃G−L≃0.(258) This steady-state condition mirrors the equilibrium reached in optically thick radiative media: microscopic fluctuations decorrelate rapidly, while their ensemble average remains constant. The geometric scale of the cluster ensemble is not set by any modification of the intrinsic resonant frequency Ω but by the transport length of the overlapping diffusive envelopes. Each elliptical hosts its own internal resonance governed by λ(ell) res ; the cluster-scale behaviour instead arises from the propagation and 99
attenuation of these internal modes through the intergalactic curvature continuum. The relevant macroscopic length scale is therefore Lmean =ceff τtr,(259) set by the effective transport time τtr that includes diffusive damping and reverberant return. Typically hundreds of kiloparsecs, Lmean determines how far curvature energy can propagate before merging statistically with the envelopes of neighbouring galaxies. Because many individual transport tails overlap within a few Lmean, the cluster behaves as a continuous curvature field whose radial symmetry emerges naturally from the convolution of exponentially decaying envelopes. The NFW-like slope of the effective density profile reflects the statistical superposition of these diffusive fields, while its amplitude is set by the reverberation gain Grev. Lenticular (S0) galaxies and transitional systems occupy an intermediate position in the transport hierarchy. Their flattened stellar distributions generate anisotropic transport envelopes, though without the global phase locking of spiral disks. Within clusters, S0s are more centrally concentrated and frequently exhibit weak radial or cluster–axis alignments, reflecting the morphology–density relation and the influence of the cluster potential. These factors naturally bias their transport tails toward the cluster interior, enhancing reverberation in the core without requiring any specific disk–orientation mechanism. Spiral galaxies, by contrast, contribute negligibly to the cluster–scale field since their coherent m= 1 modes do not sustain extended transport tails. The same radiative–transport principles extend to cosmic filaments, which may be regarded as cylindrical continuations of the cluster ensemble. Here, resonators (galaxies, groups) lie along an extended axis, and their diffusive envelopes overlap primarily in one dimension. The reduced symmetry decreases the effective gain but preserves the transport mechanism, producing elongated curvature envelopes aligned with the filament spine. In summary, a galaxy cluster behaves as a reverberant curvature–transport system: a self-averaging ensemble of overlapping diffusive envelopes whose collective behaviour reproduces the NFW form of the apparent dark-matter halo. Curvature energy is sustained statistically rather than coherently, maintained through continual diffusive exchange among the constituent galaxies. Spirals correspond to coherent resonant emitters, ellipticals to self-sustaining diffusive cavities, lenticulars to anisotropic transitional resonators, clusters to radiative–transport continua, and filaments to linear transport ensembles. Together these systems form a continuous hierarchy of curvature–matter coupling, governed throughout by the same causal transport framework. 100
9 Local Silence and the Preservation of Classical Gravity The two–channel framework developed in Section 4.1 and Section 6.2 immediately explains why all laboratory, Solar–System, and stellar tests of gravity recover the Newton–Einstein law with complete fidelity. Channel 2—the Newtonian Channel—represents the quasi–static stress–response loop and is identical to classical gravity. It determines the slowly varying gravitational potential ∇2ΦN= 4πG (ρb,0+ρeff),(260) where ρb,0 is the coarse–grained baryonic density and ρeff is the cycle–averaged curvature energy stored over cosmological timescales by Channel 1. Since Channel 2 governs all observable accelerations, it reproduces Newtonian and relativistic gravity in every regime accessible to experiment. Channel 1, the Resonant Channel, is entirely different in character. It is driven solely by the coherent component ρb,1 of the baryonic motion, 9 producing an oscillatory curvature perturbation Φ res . The response is given by the causal susceptibility of the effective field, Φres =GΩ∗ρb,1,(261) but this oscillatory field does not enter the directly observable gravitational potential. Instead, its cycle–averaged curvature energy is transferred extremely slowly into the stored component ρeff over the memory time τmem = Γ−1,(262) which exceeds stellar and galactic dynamical times by many orders of magnitude. For all practical and observational purposes, ρeff is therefore constant. The decomposition ρb=ρb,0+ρb,1,Φ = ΦN+ Φres,(263) makes clear that only Φ N governs geodesics. Channel 1 produces rapidly oscillatory curvature, but its only macroscopic imprint is the ultra–slow accumulation of stored curvature energy, ρeff ∝DGΩ∗ρb,1 2Eτmem ,(264) 9 Here “coherent” refers to the portion of ρb,1 that contributes constructively to the resonant curvature channel. In spiral galaxies this corresponds to a genuinely phase–aligned global m = 1 mode. In ellipticals, the baryonic motion is not globally phase coherent; rather, the long-wavelength retarded Green function induces short-range field coherence within each 4D speckle cell. Thus ρb,1 is globally random but locally coherent through the field response, which is the relevant sense for Channel 1. 101
which behaves observationally as a time–independent dark–mass density. This separation of timescales defines the regime of local silence. No laboratory, Solar–System, or stellar experiment can detect Φ res directly, because any change in ρeff unfolds only over the cosmological memory time (262) . Even a dramatic event—such as instantaneously removing every elliptical galaxy from a cluster— would not erase the cluster’s dark envelope: the stored curvature energy would decay only over many billions of years. To any local observer, the gravitational field would appear unchanged, continuing to follow the Newton–Einstein potential ΦN. Hence Channel 1 is not merely dynamically quiet: it is observationally silent. Its influence enters measurable dynamics only through the quasi–static stored component ρeff that feeds Channel 2. This ensures full consistency with all precision tests of Newtonian and relativistic gravity, while permitting resonant curvature to generate effective dark–mass halos only in systems where coherent baryonic motion persists over cosmic times. A useful analogy is supplied by acoustics. The acoustic field oscillates rapidly relative to the slow buildup of time–averaged pressure; what one measures is not the sound field itself but the quasi–static pressure it leaves behind. Similarly, the resonant curvature field Φ res is not observable in isolation. Only its stored cycle– averaged energy density ρeff , acting gravitationally through Channel 2, enters measurable dynamics. This captures precisely why the resonant channel remains locally silent even while shaping the large–scale gravitational environment. 10 Prediction: Absence of Gain in Homogeneous Systems This section establishes a central and falsifiable prediction of the two–channel causal–resonant framework: a perfectly homogeneous baryonic medium cannot generate curvature gain in Channel 1 under any circumstances. Inhomogeneity is required to form the source–rescatter loops, Foldy–Lax coupling, and reverberation processes that define the effective dark density ρeff . A spatially smooth system excites only Channel 2 (the Newtonian Channel) and is therefore resonantly silent. To avoid confusion, we first clarify how the linear perturbation theory of Section 4.4 relates to homogeneous media. The perturbative equation, repeated here for reference, ∇2δΦ + 1 c2 eff −∂2 t−2Γ ∂t+ Ω2δΦ = 4πG c2 eff δρb,(265) presupposes that multiple scattering has already produced an effective medium with renormalised parameters Ω, Γ, and ceff . A homogeneous configuration cannot generate such a medium and therefore cannot satisfy Eq. (265) in isolation. 102
10.1 Three Distinct Regimes for Homogeneous Clouds A homogeneous cloud may arise in three physically distinct settings: (1) Isolated homogeneous cloud. The cloud possesses a natural frequency ωnat ∝√ρ0 , but generates no scattering, no self–energy, and therefore no renormalised parameters (Ω,Γ, ceff). Thus Channel 1 does not form: Φ(isolated) res = 0.(266) (2) Cloud embedded within an externally generated effective medium but not phase–locked. The surrounding medium possesses (Ω , Γ , ceff ), but the cloud oscillates at its own ωnat and cannot couple coherently through the causal Green function. Hence no Channel 1 field is generated: Φ(embedded,free) res = 0.(267) (3) Cloud embedded within an externally generated effective medium and forced at the global resonance Ω.In this artificial situation the homogeneous cloud may oscillate in temporal phase with the background because the driving field imposes the phase externally. The cloud then satisfies the driven Helmholtz relation ∇2+k2 ΩδΦ = 4πG c2 eff δρb, k2 Ω=Ω2 c2 eff ,(268) but the scattering strength remains σ(x)≡0,(269) because temporal coherence does not create spatial contrast. Thus no Foldy– Lax loops form, no self-energy is generated, and the response remains a strictly single–pass, non-amplifying Helmholtz field. The theorem below applies to all three regimes. 10.2 Theorem: Homogeneous Media Produce No Resonant Gain Theorem (Absence of Resonant Gain in Homogeneous Systems). Let ρb (x) = ρ0 be a perfectly homogeneous baryonic distribution with no spatial contrast. Then, in each of the regimes in Section 10.1, the Channel 1 curvature response exhibits no resonant amplification: 1. Isolated case: No effective medium forms, so Channel 1 does not exist: Φres = 0.(270) 103
2. Embedded but free-running: Without phase–locking, the cloud does not couple to the resonant field: Φres = 0.(271) 3. Embedded and phase–locked: Even when driven at the global resonance Ω, the response is single–pass and non–amplifying because σ (x) ≡ 0: A=GΩ∗(σA) =⇒A= 0.(272) The reverberation gain vanishes, Grev =1+Nη −→ 1, η = 0,(273) and the stored curvature energy remains at the single–pass level: ρ(hom) eff ∼1 Nρ(inhom) eff ,(274) for systems with the same total mass, where N counts the discrete substructures available in the inhomogeneous case. Therefore a perfectly homogeneous baryonic medium cannot generate measurable curvature amplification in Channel 1 under any physical configuration. 10.3 Interpretation and Connection to Perturbation Theory The theorem shows that spatial contrast is essential for Channel 1. Without it, no Foldy–Lax loops form, no self–energy accumulates, and the collective resonant mode does not arise. This explains why the perturbative framework of Section 4.4 cannot be applied to a homogeneous system in isolation: the parameters (Ω , Γ , ceff ) exist only after rescattering has generated an effective medium. In the idealised circumstance of regime (3), where a homogeneous cloud is externally forced to oscillate at the global resonance within an already-existing effective medium, the response reduces to the driven Helmholtz convolution δΦ(x) = ZGΩ(|x−x′|)δρb(x′)d3x′,(275) but its amplitude remains faint, being limited to the single-pass level because no spatial contrast exists to sustain coherent feedback. The causal–resonant model therefore predicts that apparent dark structure correlates with baryonic contrast and temporal coherence, not with total baryonic mass: ρeff ∝(baryonic contrast) ×(coherence time).(276) 104
In cluster collisions such as the Bullet Cluster (1E 0657 − 56), where dense, collisionless galaxies retain their spatial contrast while the diffuse plasma is stripped away, the curvature resonance must remain anchored to the galaxy component rather than to the displaced gas. This behaviour precisely matches observations without invoking non–baryonic matter, and yields a decisive test: if lensing peaks consistently coincide with regions of greatest baryonic contrast, the causal– resonant interpretation gains strong empirical support. 11 Emergence of Multiple Scattering from a Weakly Inhomogeneous Continuous Medium 11.1 A Brief Review of Garvitational Multiple Scatter Theory Before attempting to analyse the early Universe, it is useful to summarise the assumptions that underlie the effective–medium framework developed thus far. Throughout previous sections we have considered a baryonic mass distribution ρb(x) = N X j=1 mjδ(x−xj),(277) consisting of spatially separated stellar or galactic constituents. Within a coherent domain these constituents oscillate at a common curvature frequency Ω2∝G¯ρb,(278) where ¯ρbis the root–mean baryonic density over the coherent volume. Each baryonic element responds to the resonant curvature field as a damped oscillator, ¨ Xj+ 2Γ ˙ Xj+ Ω2Xj=αΦres(xj),(279) which reinjects curvature waves back into the medium. Substituting the induced oscillators into the Green–function representation yields the self–consistent Foldy– Lax equation, Φres(x)=Φinc(x) + ZGΩ(|x−x′|)σ(x′) Φres(x′)d3x′,(280) where the contrast–defined scattering strength is σ(x)∝δρb(x).(281) Applying the Foldy–Lax operator to both sides of Eq. (280) produces the renormalised effective–medium wave equation, ∇2Φres +1 c2 eff −∂2 t−2Γ ∂t+ Ω2Φres =4πG c2 eff ρb,1,(282) 105
familiar with when I began this project, yet essential for deriving the behaviours of elliptical galaxies and clusters. I sincerely hope I have not misinterpreted their theories. I gratefully acknowledge the support of Research Ireland / Science Foundation Ireland and the European Union, whose funding over the past twenty years enabled my research in optics and allowed me to deepen my understanding of physics. Last and most importantly, I thank my wife Karen and my daughters for their unwavering support as I wrestled with this theory over these past several months. 112
A Appendix: Bessel and Hankel Functions Used in the Main Text This appendix is provided purely for convenience and reference. It contains no new results or derivations but simply summarizes the definitions and basic properties of the Bessel and Hankel functions employed throughout Section 7.4–Section 8.1.4. Two distinct families appear in the resonant–curvature formalism: 1. the cylindrical Bessel functions Jm ( z ) and Ym ( z ), which solve the two–dimensional Helmholtz equation relevant to thin–disk geometries, and 2. the spherical (or “normal”) Bessel functions jℓ ( z ), which arise in three–dimensional interference or isotropic–propagation problems. Both sets are standard special functions with well–known analytic forms. A.1 Cylindrical Bessel and Hankel Functions The cylindrical Bessel functions of the first and second kind, Jm ( z ) and Ym ( z ), are the two linearly independent solutions of Bessel’s differential equation, d2 dz2+1 z d dz +1−m2 z2um(z)=0.(A.1) Jm ( z ) is finite at z = 0 and represents the regular solution, while Ym ( z ) diverges logarithmically and represents the singular one. For small arguments (z≪1), Jm(z)≃1 m!z 2m, Ym(z)≃ −(m−1)! π2 zm ,(A.2) and for large arguments (z≫1), Jm(z)≃r2 πz cosz−mπ 2−π 4, Ym(z)≃r2 πz sinz−mπ 2−π 4. (A.3) The Hankel functions of the first and second kind, H(1) m(z) = Jm(z) + iYm(z), H(2) m(z)=Jm(z)−iYm(z),(A.4) represent outgoing and incoming cylindrical waves, respectively. They satisfy the same differential equation as Eq. (A.1) and are particularly convenient for enforcing causal (Sommerfeld) boundary conditions. Their large–argument forms, H(1) m(z)≃r2 πz ei(z−mπ 2−π 4), H(2) m(z)≃r2 πz e−i(z−mπ 2−π 4),(A.5) show explicitly that H(1) m corresponds to an outward–propagating (radiative) wave and H(2) m to an inward–propagating one. In the main text, H(1) 1 ( kΩr ) represents the outward, causally damped response of the resonant curvature field in Eq. (152). 113
A.2 Spherical Bessel Functions The spherical Bessel functions jℓ ( z ) arise as the radial eigenfunctions of the three–dimensional Helmholtz equation in spherical coordinates, d2 dz2+2 z d dz +1−ℓ(ℓ+ 1) z2uℓ(z) = 0.(A.6) They are related to the cylindrical functions by jℓ(z) = rπ 2zJℓ+1/2(z),(A.7) and describe standing or radiating waves in isotropic three–dimensional media. The lowest orders are j0(z) = sin z z, j1(z) = sin z z2−cos z z,(A.8) which govern the monopole and dipole components, respectively. In Section 7.4, j0 defines the symmetric interference term between two coherent curvature sources, while j1 captures the first azimuthal modulation responsible for the observed spiral–mode pattern. A.3 Physical Interpretation Bessel and Hankel functions provide the analytic basis for describing wave propagation in both planar and spherical geometries: •Jm ( kΩr ) and H(1) m ( kΩr ) describe oscillatory and radiative behaviour within the thin, nearly–planar baryonic disk. •jℓ ( kΩr ) describes isotropic or quasi–spherical propagation relevant to interference, scattering, and large–scale resonant coupling. Together these functions form the mathematical bridge between the two–dimensional spiral disk and the three–dimensional resonant halo. They are summarized here solely for convenience; the appendix introduces no new analytical content beyond the main text. 114
B Appendix: Code for Simulating the 3D Thin–Disk Resonant Halo This appendix provides the Python implementation used to generate the three dimensional Helmholtz resonant halo simulations presented in the main text (Section 8.1.8; Figs. 4–6). The code evaluates the full integral solution of the driven three dimensional Helmholtz equation for the coherent m =1 curvature mode of a rotating exponential disk, including vertical stratification, damping, and environmental amplification. The stationary curvature envelope A(r, θ, z) satisfies the Helmholtz equation (∇2+k2 Ω)A=4πG c2 eff ρb(r, θ, z), where ρb is the baryonic source density and kΩ = (Ω + i Γ) /ceff is the complex resonant wavenumber. The coherent azimuthal dependence is fixed by the factor eiθ selecting the global m =1 mode. Vertical structure is imposed through an exponential confinement exp ( −|z|/Lz ), reflecting the evanescent vertical branch of the dispersion relation. The resonant wavelength and damping length, expressed in terms of the effective propagation speed ceff and damping rate Γ, are λres =2πceff Ω, Ldamp =ceff Γ, and are implemented directly as physical scales in the simulation. The local resonant frequency Ω = v0/r0 is fixed entirely by the observed rotation curve of the galaxy and is therefore not a free parameter. Likewise, the propagation speed ceff is chosen so that the resulting resonant wavelength lies in the observed 40–60 kpc range; for Ω ≃ 8 . 77 × 10 −16 s−1 this requires ceff ∼ 10 5m s−1 , and we adopt ceff = 2 × 10 5m s−1 . With these fixed, only four physical inputs govern the global mode: 1. Resonant frequency Ω=v0/r0(fixed by the rotation curve). 2. Quality factor Q = Ω / (2Γ) (sets the damping rate and halo extent Ldamp ). 3. Coupling efficiency X with |X| ≤ 1 (sets the mean single–scatterer response). 4. Geometric projection factor ϵ (weights the m = 1 azimuthal component of the disk). Once these parameters are specified, the integral solution determines the complete three–dimensional structure of the resonant halo without further normalisation. The midplane envelope Ar ( r ) is computed directly from the integral kernel of the 115
Helmholtz operator; the vertical profile multiplies this by exp ( −|z|/Lz ) without altering its radial amplitude. The corresponding effective energy density scales as ρeff(r, θ, z)∝Q2|X|2N⋆|A(r, θ, z)|2e−2|z|/Lz, consistent with the amplitude–level Foldy–Lax gain √N⋆ and the Q –dependent buildup of stored curvature energy. Integrated quantities—the effective column density, enclosed effective mass Meff ( r ), and total circular velocity vtot ( r )—are calculated directly from this field. This implementation reproduces the radial profiles, midplane maps, and orthogonal slices shown in Figs. 4–6 and provides a complete numerical realisation of the 3D resonant halo model discussed in Section 8.1.8. 1# ================================================================ 2# 3D Helmholtz Resonant Halo Model m = 1 (LaguerreGaussian Mode) 3# --------------------------------------------------------------- 4# Author: Bryan Hennelly 5# Date: 25 October 2025 6# --------------------------------------------------------------- 7# Implements the full threedimensional causalresonant halo model 8# described in the main text (Sections full-3D-with-gain, 9# stored-energy-and-gain). All symbols, definitions, and physical 10 # meanings correspond directly to those in the book. 11 # 12 # The model solves the BesselHankel form of the 2D Helmholtz equation 13 # for the coherent m = 1 mode of a thin rotating galactic disk, then 14 # extends it vertically as an exponential layer of scale L_z. 15 # 16 # --- PARAMETER TUNING SUMMARY ---------------------------------- 17 # Omega : Angular frequency v0 / r0 fixed by rotation curve scale. 18 # Q : Resonance quality factor adjusts the outer halo tail. 19 # X : Mean coupling efficiency (0 X 1) sets effective gain. 20 # eps : Geometric m = 1 projection factor tunes overall amplitude. 21 # --------------------------------------------------------------- 22 # These four govern amplitude and spatial morphology. 23 # --------------------------------------------------------------- 24 # 25 # --- NORMALIZATION NOTE ---------------------------------------- 26 # The field A_r(r) is *already depth-normalized*: it represents the 27 # midplane envelope of the full 3D mode. No additional normalization 28 # with L_z is required or applied. The vertical extent of the stored 29 # curvature energy is represented by multiplying the midplane energy 30 # density by the effective column depth L_damp (optionally replaced 31 # by L_z in tests). Normalizing with 1/sqrt(L_z) would distort the 32 # A_r amplitude and is intentionally avoided here. 33 # --------------------------------------------------------------- 34 # 35 # --- L_z INTERPRETATION ---------------------------------------- 36 # The vertical scale height L_z is a geometric confinement parameter 37 # (typically L_z ~ 0.3 * L_damp). It sets the half-thickness of the 38 # exponential envelope exp(-|z|/L_z) used for visualizations. This is 39 # *not* a normalization constant but a physical e-folding scale. 40 # --------------------------------------------------------------- 116
41 # ================================================================ 42 43 import numpy as np 44 import mpmath as mp 45 import matplotlib.pyplot as plt 46 47 # ---------- Global plot style ---------- 48 plt.rcParams.update({ 49 ’font.size’: 16, 50 ’axes.titlesize’: 18, 51 ’axes.labelsize’: 16, 52 ’xtick.labelsize’: 14, 53 ’ytick.labelsize’: 14, 54 ’legend.fontsize’: 14 55 }) 56 57 # ---------- Physical constants ---------- 58 G = 6.67430e-11 59 Msun = 1.98847e30 60 pc = 3.085677581491367e16 61 kpc = 1.0e3 * pc 62 c_light = 2.99792458e8 63 64 # ================================================================ 65 # USER-TUNABLE MODEL PARAMETERS 66 # ------------------------------------------------ 67 # These are the only physically meaningful dials: 68 # Omega (fixed by v0, r0), Q, X, and eps. 69 # ================================================================ 70 v0 = 230e3 # flat rotation speed [m/s] 71 r0 = 8.5 * kpc # reference radius [m] 72 Omega = v0 / r0 # angular frequency [s^-1] 73 N_star = 6.0e10 # number of baryonic emitters (stars) 74 75 c_eff = 2.0e5 # effective wave speed [m/s] 76 Q = 12.6 # quality factor (temporal coherence) 77 X = 1 # mean coupling efficiency (|X| 1) 78 G_env_root = np.sqrt(X*X*N_star) # amplitude-level environmental gain 79 Gamma = Omega / (2 * Q) # damping rate 80 eps = 0.3 # geometric m=1 projection factor 81 # ================================================================ 82 83 # ---------- Derived scales ---------- 84 k_complex = (Omega + 1j * Gamma) / c_eff 85 lambda_res = 2 * np.pi * c_eff / Omega 86 L_damp = c_eff / Gamma 87 L_z = 0.3 * L_damp # vertical confinement scale << L_damp 88 89 print(f"Omega= {Omega}, L_damp = {L_damp/kpc:.2f} kpc, _res = ,→{lambda_res/kpc:.2f} kpc, L_z = {L_z/kpc:.2f} kpc") 90 91 # ---------- Exponential baryonic disk ---------- 92 M_disk_total = N_star * Msun 93 R_d = 3.0 * kpc 117
94 h_z = 0.30 * kpc 95 Sigma0 = M_disk_total / (2.0 * np.pi * R_d**2) 96 97 def Sigma_r(r): return Sigma0 * np.exp(-r / R_d) 98 def rho_b_midplane(r): return Sigma_r(r) / (2.0 * h_z) 99 100 # ---------- Radial grid ---------- 101 r_max = 1000.0 * kpc 102 N = 1800 103 r = np.linspace(1.0e-3 * kpc, r_max, N) 104 105 # ---------- Special functions ---------- 106 def J1(z): return mp.besselj(1, z) 107 def H1(z): return mp.hankel1(1, z) 108 109 J1_vals = np.array([complex(J1(k_complex * ri)) for ri in r]) 110 H1_vals = np.array([complex(H1(k_complex * ri)) for ri in r]) 111 112 # ---------- Source moment ---------- 113 M1_vol = 2.0 * np.pi * eps * rho_b_midplane(r) 114 115 # ---------- Complex trapezoidal integration ---------- 116 def cumtrapz_complex(y, x): 117 out = np.zeros_like(y, dtype=complex) 118 for iin range(1, len(x)): 119 out[i] = out[i-1] + 0.5 * (y[i] + y[i-1]) * (x[i] - x[i-1]) 120 return out 121 122 integrand_J = r * M1_vol * J1_vals 123 integrand_H = r * M1_vol * H1_vals 124 125 I_J_cum = cumtrapz_complex(integrand_J, r) 126 I_H_cum = cumtrapz_complex(integrand_H, r) 127 I_H_total = I_H_cum[-1] 128 129 # ---------- Envelope field A(r) ---------- 130 prefac = 4.0 * np.pi**2 * 1j * G / (c_eff**2) 131 A_r = prefac * (J1_vals * (I_H_total - I_H_cum) + H1_vals * I_J_cum) 132 133 # ---------- Apply coherence gain ---------- 134 # Environmental amplification (Q (X N_star)) 135 gain_field = Q * G_env_root 136 A_r *= gain_field 137 138 # ---------- Effective and baryonic densities ---------- 139 # Field energy density |A|, already depth-normalized. 140 rho_eff_mid = (2*Omega**2 / (8.0 * np.pi * G * c_light**2)) * np.abs(A_r)**2 141 rho_b = rho_b_midplane(r) 142 143 # ---------- Vertically averaged surface densities ---------- 144 # The vertical column depth is represented by L_damp (no extra L_z ,→normalization) 145 Sigma_eff = L_damp * rho_eff_mid 146 118
147 # ---------- Enclosed masses ---------- 148 M_b, M_eff = np.zeros_like(r), np.zeros_like(r) 149 for iin range(1, len(r)): 150 dM_b = 2 * np.pi * 0.5 * (r[i]*rho_b[i] + r[i-1]*rho_b[i-1]) * ,→(r[i]-r[i-1]) * (2*h_z) 151 dM_eff = 2 * np.pi * 0.5 * (r[i]*Sigma_eff[i] + r[i-1]*Sigma_eff[i-1]) * ,→(r[i]-r[i-1]) 152 M_b[i] = M_b[i-1] + dM_b 153 M_eff[i] = M_eff[i-1] + dM_eff 154 155 print(f"Check: integrated baryonic mass = {M_b[-1]/Msun:.2e} Msun (target ,→{M_disk_total/Msun:.2e})") 156 157 158 # ================================================================ 159 # FIGURE 1 Radial Profiles and Rotation Curves 160 # ================================================================ 161 162 r_kpc = r / kpc 163 fig, axs = plt.subplots(2, 2, figsize=(11, 9)) 164 plt.subplots_adjust(wspace=0.35, hspace=0.35) 165 166 plt.rcParams.update({ 167 ’font.size’: 21, 168 ’axes.titlesize’: 23, 169 ’axes.labelsize’: 21, 170 ’xtick.labelsize’: 19, 171 ’ytick.labelsize’: 19, 172 ’legend.fontsize’: 19 173 }) 174 175 # (a) Amplitude |A| 176 axs[0,0].plot(r_kpc, np.abs(A_r), color=’black’) 177 axs[0,0].set_title("(a) Amplitude $|A|$") 178 axs[0,0].set_xlabel("r [kpc]") 179 axs[0,0].set_ylabel(r"$|A(r)|$") 180 181 # (b) Mass density 182 Msun_pc3 = (pc**3) / Msun 183 axs[0,1].plot(r_kpc, rho_eff_mid*Msun_pc3, label=r"$\rho_{\rm eff}$", ,→color=’blue’) 184 axs[0,1].plot(r_kpc, rho_b*Msun_pc3, label=r"$\rho_{\rm b}$", ,→color=’orange’, linestyle=’--’) 185 axs[0,1].set_yscale(’log’) 186 axs[0,1].set_ylim(1e-8, 1e2) 187 axs[0,1].set_title("(b) Mass density") 188 axs[0,1].set_xlabel("r [kpc]") 189 axs[0,1].set_ylabel(r"Density [$M_\odot\,{\rm pc}^{-3}$]") 190 axs[0,1].legend() 191 192 # (c) Rotation curves 193 limit = r_kpc <= 50 194 axs[1,0].plot(r_kpc[limit], np.sqrt(G*M_b[limit]/r[limit])/1e3, ’--’, ,→label=r"$v_{\rm b}$", color=’orange’) 119
195 axs[1,0].plot(r_kpc[limit], np.sqrt(G*M_eff[limit]/r[limit])/1e3, ,→label=r"$v_{\rm eff}$", color=’blue’) 196 axs[1,0].plot(r_kpc[limit], ,→np.sqrt(G*(M_b[limit]+M_eff[limit])/r[limit])/1e3, label=r"$v_{\rm ,→tot}$", color=’black’) 197 axs[1,0].set_title("(c) Rotation curves (50 kpc)") 198 axs[1,0].set_xlabel("r [kpc]") 199 axs[1,0].set_ylabel("Velocity [km s$^{-1}$]") 200 axs[1,0].legend(loc=’center right’, frameon=True, facecolor=’white’, ,→framealpha=0.3, edgecolor=’none’, fontsize=16) 201 202 # (d) Mass ratio 203 axs[1,1].plot(r_kpc, M_eff / np.maximum(M_b, 1e-30), color=’purple’) 204 axs[1,1].set_title("(d) Enclosed mass ratio") 205 axs[1,1].set_xlabel("r [kpc]") 206 axs[1,1].set_ylabel(r"$M_{\rm eff}/M_{\rm b}$") 207 208 fig.suptitle("3D Helmholtz Halo Model (m = 1 Mode)", 209 fontsize=26, fontweight=’bold’, y=0.98) 210 plt.tight_layout(rect=[0, 0, 1, 0.96]) 211 plt.savefig("C:/Users/BryanH/Downloads/fig1_spiral_profiles.png", dpi=300, ,→facecolor=’white’) 212 plt.show() 213 214 # ================================================================ 215 # FIGURE 2 2D Midplane Maps (Amplitude, Phase, and Densities) 216 # ================================================================ 217 218 Nmap = 400 219 x = np.linspace(-50, 50, Nmap) * kpc 220 y = np.linspace(-50, 50, Nmap) * kpc 221 X, Y = np.meshgrid(x, y) 222 R = np.sqrt(X**2 + Y**2) 223 TH = np.arctan2(Y, X) 224 225 A_radial = np.interp(R, r, A_r) 226 rho_b_2D = np.interp(R, r, rho_b) 227 # --- corrected factor of 2 for consistency with rho_eff_mid --- 228 rho_eff_2D = (2 * Omega**2 / (8.0 * np.pi * G * c_light**2)) * ,→np.abs(A_radial)**2 229 A_2D = A_radial * np.exp(1j * TH) 230 231 amp_map = np.abs(A_2D) 232 phase_map = (np.angle(A_2D) + np.pi) / (2.0 * np.pi) 233 rho_tot_2D = rho_b_2D + rho_eff_2D 234 235 amp_norm = amp_map / np.max(amp_map) 236 phase_norm = phase_map 237 rho_b_norm = rho_b_2D / np.max(rho_b_2D) 238 rho_eff_norm = rho_eff_2D / np.max(rho_eff_2D) 239 rho_tot_norm = rho_tot_2D / np.max(rho_tot_2D) 240 241 fig, axs = plt.subplots(2, 2, figsize=(10, 10)) 242 120
243 axs[0,0].imshow(amp_norm, extent=[-50, 50, -50, 50], origin=’lower’, ,→cmap=’viridis’) 244 axs[0,0].set_title("(a) |A| amplitude\n(normalized)", fontsize=20) 245 axs[0,0].axis(’off’) 246 247 axs[0,1].imshow(phase_norm, extent=[-50, 50, -50, 50], origin=’lower’, ,→cmap=’gray’) 248 axs[0,1].set_title("(b) Phase arg(A)\n(normalized)", fontsize=20) 249 axs[0,1].axis(’off’) 250 251 rgb = np.zeros((*rho_eff_norm.shape, 3)) 252 rgb[..., 0] = rho_b_norm 253 rgb[..., 1] = rho_eff_norm 254 axs[1,0].imshow(rgb, extent=[-50, 50, -50, 50], origin=’lower’, vmin=0, ,→vmax=1) 255 axs[1,0].set_title("(c) Baryonic (red)\n and resonant (green)", ,→fontsize=20) 256 axs[1,0].axis(’off’) 257 258 axs[1,1].imshow(rho_tot_norm, extent=[-50, 50, -50, 50], origin=’lower’, ,→cmap=’gray’) 259 axs[1,1].set_title("(d) Combined density\n (normalized)", fontsize=20) 260 axs[1,1].axis(’off’) 261 262 fig.suptitle("Spiral Helmholtz Solutions\nNormalized Results (100 kpc 100 ,→kpc)", 263 fontsize=22, fontweight=’bold’, y=0.93) 264 plt.tight_layout(rect=[0, 0, 1, 0.95]) 265 plt.subplots_adjust(wspace=-0.3, hspace=0.25) 266 plt.savefig("C:/Users/BryanH/Downloads/fig2_spiral_maps.png", dpi=300, ,→facecolor=’white’) 267 plt.show() 268 269 # ================================================================ 270 # FIGURE 3 Dark Mass Density Slices (XY and XZ) 271 # ================================================================ 272 # The vertical XZ slice shows exponential decay exp(-2|z|/L_z), 273 # where L_z represents the geometric confinement scale. The vertical 274 # contrast therefore illustrates the halos stratification typically 275 # a few 10 kpc thick and is not rescaled by any normalization factor. 276 # ================================================================ 277 278 Nxy = 600 279 WIDTH = 200 280 extent_xy = [-WIDTH/2, WIDTH/2, -WIDTH/2, WIDTH/2] 281 extent_xz = [-WIDTH/2, WIDTH/2, -WIDTH/2, WIDTH/2] 282 283 x = np.linspace(-WIDTH/2, WIDTH/2, Nxy) * kpc 284 y = np.linspace(-WIDTH/2, WIDTH/2, Nxy) * kpc 285 z = np.linspace(-WIDTH/2, WIDTH/2, Nxy) * kpc 286 287 # --- Horizontal (XY) slice at z = 0 --- 288 X, Y = np.meshgrid(x, y) 289 R_xy = np.sqrt(X**2 + Y**2) 121
127 # STELLAR BARYONIC INPUT 128 # ============================================================ 129 130 M_tot = 1.0e11 * Msun 131 m_star = 1.0 * Msun 132 N_star = M_tot/m_star 133 Re = 5*kpc 134 a_h = Re/1.8153 135 rho0 = M_tot/(2*np.pi*a_h**3) 136 V_gal = 2*np.pi*a_h**3 137 138 v_rms = 200e3 139 alpha = 1.0 140 X = 1.0 # |X| = 1 for now 141 142 print("\n--- Galaxy (baryons) ---") 143 print(f"M_tot = {M_tot/Msun:.2e} Msun") 144 print(f"N_star = {N_star:.2e}") 145 print(f"V_gal = {V_gal/(kpc**3):.2e} kpc^3") 146 print(f"v_rms = {v_rms/1e3:.1f} km/s") 147 print(f"alpha = {alpha:.2f}") 148 149 # ============================================================ 150 # TRANSPORTENHANCED HALO DENSITY (grouped analytic form) 151 # ============================================================ 152 153 # Environmental gain (for diagnostic only; not changing normalization) 154 G_env_eff = (X**2 155 * (N_star**2 / V_gal) 156 * alpha 157 * (lambda_res / (2*pi))**2 158 * (v_rms / (2*Gamma_eff))) 159 160 # Working prefactor (kept identical to original numerics) 161 K_base = (G * alpha * v_rms) / (2.0 * c**2 * c_eff**3) \ 162 * L_damp_eff * Q_eff \ 163 * (N_star**2 * m_star**2 / V_gal) * X**2 164 165 def rho_eff_far(r): 166 return Xi_nabla * K_base * np.exp(-2.0*r/L_damp_eff) / ((4.0*pi)**2 * ,→r**2) 167 168 # ============================================================ 169 # INNER-CORE REGULARIZATION ( r for r < R_match) 170 # ============================================================ 171 172 R_match = 5.0 * kpc 173 rho_match = rho_eff_far(R_match) 174 slope = rho_match / R_match 175 176 def rho_eff_piecewise(r): 177 r = np.asarray(r, dtype=float) 178 return np.where(r < R_match, slope*r, rho_eff_far(r)) 179 128
180 # ============================================================ 181 # OPTIONAL: Hernquist baryons (for comparison) 182 # ============================================================ 183 184 def rho_b_hernquist(r): 185 r = np.asarray(r, dtype=float) 186 return (M_tot / (2.0 * np.pi)) * (a_h / (r * (r + a_h)**3)) 187 188 def M_b_enclosed(r): 189 r = np.asarray(r, dtype=float) 190 return M_tot * (r**2) / (r + a_h)**2 191 192 # ============================================================ 193 # RADIAL GRID + ENSEMBLE MASS INTEGRATION 194 # ============================================================ 195 196 r = np.geomspace(0.05 * kpc, 2.0e3 * kpc, 2000) 197 rho_eff_arr = rho_eff_piecewise(r) 198 rho_b_arr = rho_b_hernquist(r) 199 200 M_eff = np.zeros_like(r) 201 for iin range(1, len(r)): 202 r0, r1 = r[i-1], r[i] 203 y0 = 4.0 * np.pi * r0**2 * rho_eff_arr[i-1] 204 y1 = 4.0 * np.pi * r1**2 * rho_eff_arr[i] 205 M_eff[i] = M_eff[i-1] + 0.5 * (y0 + y1) * (r1 - r0) 206 207 M_b = M_b_enclosed(r) 208 209 ratio_rho = rho_eff_arr / np.maximum(rho_b_arr, 1e-99) 210 ratio_M = M_eff / np.maximum(M_b, 1e-99) 211 212 # ============================================================ 213 # DIAGNOSTICS 214 # ============================================================ 215 216 print("\n--- Diagnostics ---") 217 print(f"G_env_eff = {G_env_eff:.3e}") 218 print(f"Q_eff = {Q_eff:.3e}") 219 print(f"Xi_nabla = {Xi_nabla:.3e}") 220 print(f"K_base = {K_base:.3e} [SI, identical to original]") 221 print(f"rho_eff(R_match)= {rho_eff_piecewise(R_match):.3e} kg m^-3") 222 idx_100kpc = np.searchsorted(r, 100*kpc) 223 print(f"M_eff(100 kpc) = {float(M_eff[idx_100kpc])/Msun:.3e} Msun") 224 225 # ============================================================ 226 # Plots 227 # ============================================================ 228 229 plt.rcParams.update({ 230 ’font.size’: 21, 231 ’axes.titlesize’: 23, 232 ’axes.labelsize’: 21, 233 ’xtick.labelsize’: 19, 129
234 ’ytick.labelsize’: 19, 235 ’legend.fontsize’: 19 236 }) 237 238 r_kpc = r/kpc 239 240 # ================================================================ 241 # FIGURE 1 Radial Profiles and Rotation Curves (Matched Formatting) 242 # ================================================================ 243 244 fig, axs = plt.subplots(2, 2, figsize=(11, 9)) 245 plt.subplots_adjust(wspace=0.35, hspace=0.35) 246 247 # (a) Effective and baryon densities 248 Msun_pc3 = Msun/(pc**3) 249 axs[0,0].loglog(r_kpc, rho_b_arr/Msun_pc3, ’--’, 250 label=r’$\rho_{\rm b}$’, color=’orange’) 251 axs[0,0].loglog(r_kpc, rho_eff_arr/Msun_pc3, 252 label=r’$\langle\rho_{\rm eff}\rangle$’, color=’blue’) 253 axs[0,0].set_title("(a) Mass density") 254 axs[0,0].set_xlabel("r [kpc]") 255 axs[0,0].set_ylabel(r"Density [$M_\odot\,{\rm pc}^{-3}$]") 256 axs[0,0].legend() 257 258 # (b) Density ratio 259 axs[0,1].semilogx(r_kpc, ratio_rho, color=’purple’) 260 axs[0,1].set_title("(b) Density ratio") 261 axs[0,1].set_xlabel("r [kpc]") 262 axs[0,1].set_ylabel(r"$\rho_{\rm eff}/\rho_{\rm b}$") 263 264 # (c) Rotation curves 265 limit = r_kpc <= 100.0 266 axs[1,0].plot(r_kpc[limit], np.sqrt(G*M_b[limit]/r[limit])/1e3, 267 ’--’, label=r’$v_{\rm b}$’, color=’orange’) 268 axs[1,0].plot(r_kpc[limit], np.sqrt(G*M_eff[limit]/r[limit])/1e3, 269 label=r’$v_{\rm eff}$’, color=’blue’) 270 axs[1,0].plot(r_kpc[limit], 271 np.sqrt(G*(M_b[limit]+M_eff[limit])/r[limit])/1e3, 272 label=r’$v_{\rm tot}$’, color=’black’) 273 axs[1,0].set_title("(c) Rotation curves (100 kpc)") 274 axs[1,0].set_xlabel("r [kpc]") 275 axs[1,0].set_ylabel("Velocity [km s$^{-1}$]") 276 axs[1,0].legend() 277 278 # (d) Enclosed mass ratio 279 axs[1,1].semilogx(r_kpc, ratio_M, color=’green’) 280 axs[1,1].set_title("(d) Enclosed mass ratio") 281 axs[1,1].set_xlabel("r [kpc]") 282 axs[1,1].set_ylabel(r"$M_{\rm eff}/M_{\rm b}$") 283 284 fig.suptitle("Elliptical Transport Halo: Radial Dynamics", 285 fontsize=26, fontweight=’bold’, y=0.98) 286 plt.tight_layout(rect=[0, 0, 1, 0.95]) 287 plt.savefig("C:/Users/BryanH/Downloads/elliptical_fig1_radial.png", 130
288 dpi=300, facecolor=’white’) 289 plt.show() 290 291 # ================================================================ 292 # FIGURE 2 2D Midplane: Composite & Total (LOG SCALE, 50 kpc) 293 # ================================================================ 294 295 Nmap = 600 296 extent = 50.0 # kpc 297 x = np.linspace(-extent, extent, Nmap) * kpc 298 y = np.linspace(-extent, extent, Nmap) * kpc 299 X, Y = np.meshgrid(x, y) 300 R = np.sqrt(X**2 + Y**2) 301 302 # Interpolate baryonic and effective densities on the grid 303 rho_b_2D = np.interp(R, r, rho_b_arr) 304 rho_eff_2D = np.interp(R, r, rho_eff_arr) 305 rho_tot_2D = rho_b_2D + rho_eff_2D 306 307 # Avoid zeros for log scaling 308 eps = 1e-40 309 rho_b_2D = np.maximum(rho_b_2D, eps) 310 rho_eff_2D = np.maximum(rho_eff_2D, eps) 311 rho_tot_2D = np.maximum(rho_tot_2D, eps) 312 313 # Shared physical log domain for baryonic and transport-dark maps 314 rho_min = min(rho_b_2D.min(), rho_eff_2D.min()) 315 rho_max = max(rho_b_2D.max(), rho_eff_2D.max()) 316 317 def log_norm(arr): 318 return np.clip( 319 np.log10(arr / rho_min) / np.log10(rho_max / rho_min), 320 0.0, 1.0 321 ) 322 323 Rb = log_norm(rho_b_2D) # Red channel 324 Gd = log_norm(rho_eff_2D) # Green channel 325 326 # Composite RGB 327 RGB = np.dstack((Rb, Gd, np.zeros_like(Rb))) 328 329 # Total density: its own grayscale log stretch 330 log_tot = log_norm(rho_tot_2D) 331 log_tot = log_tot / np.max(log_tot) 332 333 fig, axs = plt.subplots(1, 2, figsize=(11, 6)) 334 335 # (a) Composite map 336 axs[0].imshow(RGB, 337 extent=[-extent, extent, -extent, extent], 338 origin=’lower’, vmin=0, vmax=1) 339 axs[0].set_title("(a) Baryonic Mass (red)\n and Dark Mass (green)", 340 fontsize=20) 341 axs[0].axis(’off’) 131
342 axs[0].set_aspect(’equal’) 343 344 # (b) Total mass density 345 axs[1].imshow(log_tot, 346 extent=[-extent, extent, -extent, extent], 347 origin=’lower’, cmap=’gray’, vmin=0, vmax=1) 348 axs[1].set_title("(b) Combined density\n (normalized)", fontsize=20) 349 axs[1].axis(’off’) 350 axs[1].set_aspect(’equal’) 351 352 fig.suptitle("Elliptical Midplane Density (100 kpc 100 kpc) Log Scale", 353 fontsize=22, fontweight=’bold’, y=0.999) 354 355 fig.subplots_adjust(left=0.02, right=0.98, 356 bottom=0.03, top=0.82, wspace=0.10) 357 358 plt.savefig("C:/Users/BryanH/Downloads/elliptical_fig2_midplane_50kpc.png", 359 dpi=300, facecolor=’white’) 360 plt.show() 361 362 # ================================================================ 363 # FIGURE 3 Resonant Dark Mass Density Slices (LOG SCALE, 1 Mpc) 364 # ================================================================ 365 366 Nxy = 700 367 extent = 1000.0 # kpc (1 Mpc) 368 369 # Grids 370 x = np.linspace(-extent, extent, Nxy) * kpc 371 y = np.linspace(-extent, extent, Nxy) * kpc 372 z = np.linspace(-extent, extent, Nxy) * kpc 373 374 # ---------- XY slice ---------- 375 X, Y = np.meshgrid(x, y) 376 R_xy = np.sqrt(X**2 + Y**2) 377 rho_xy = np.interp(R_xy, r, rho_eff_arr) 378 379 # ---------- XZ slice ---------- 380 Xz, Z = np.meshgrid(x, z) 381 R_xz = np.sqrt(Xz**2 + Z**2) # spherical symmetry for elliptical case 382 rho_xz = np.interp(R_xz, r, rho_eff_arr) 383 384 # ---------- Shared log normalization ---------- 385 eps = 1e-99 386 rho_xy = np.maximum(rho_xy, eps) 387 rho_xz = np.maximum(rho_xz, eps) 388 389 rho_min = min(rho_xy.min(), rho_xz.min()) 390 rho_max = max(rho_xy.max(), rho_xz.max()) 391 392 def log_norm2(arr): 393 return np.clip(np.log10(arr / rho_min) / np.log10(rho_max / rho_min), 394 0.0, 1.0) 395 132
396 img_xy = log_norm2(rho_xy) 397 img_xz = log_norm2(rho_xz) 398 399 fig, axs = plt.subplots(1, 2, figsize=(12, 6)) 400 401 # (a) Horizontal slice XY 402 axs[0].imshow(img_xy, extent=[-extent, extent, -extent, extent], 403 origin=’lower’, cmap=’gray’, vmin=0, vmax=1) 404 axs[0].set_title("(a) Dark mass density\n Horizontal slice XY", 405 fontsize=20) 406 axs[0].axis(’off’) 407 axs[0].set_aspect(’equal’, adjustable=’box’) 408 409 # (b) Vertical slice XZ 410 axs[1].imshow(img_xz, extent=[-extent, extent, -extent, extent], 411 origin=’lower’, cmap=’gray’, vmin=0, vmax=1) 412 axs[1].set_title("(b) Dark mass density\n Vertical slice XZ", 413 fontsize=20) 414 axs[1].axis(’off’) 415 axs[1].set_aspect(’equal’, adjustable=’box’) 416 417 fig.suptitle("Elliptical Transport Halo Log-scaled Density Slices (1 Mpc)", 418 fontsize=22, fontweight=’bold’, y=0.999) 419 420 fig.subplots_adjust(left=0.02, right=0.98, 421 bottom=0.03, top=0.82, wspace=0.10) 422 423 plt.savefig("C:/Users/BryanH/Downloads/elliptical_fig3_slices_1Mpc.png", 424 dpi=300, facecolor=’white’) 425 plt.show() Listing 2: Python code for the ensemble-averaged resonant halo of an elliptical galaxy. — 133
D Appendix: Code for Simulating Cluster Galaxies This appendix provides the complete Python implementation used to generate Figures 10–12. The code follows directly from the radiative–transport formalism developed in Section 8.3.1, implementing the overlapping elliptical–halo convolution in reduced “2.5D” form for numerical stability. Each galaxy is represented by an exponential transport kernel K(s;L) = e−2s/L (4π)2s2,(D.1) ensemble–averaged over a lognormal distribution of damping lengths L with logarithmic dispersion σln L . The angularly averaged 2.5D convolution is evaluated as ρshape(r)=2πZ∞ 0 ngal(r′)r′⟨K(s)⟩µdr′, s2=r2+r′2−2rr′µ, (D.2) using Gauss–Legendre quadrature for the µ–integration. The galaxy number–density profile adopts an NFW-like form with a smooth splashback taper, ngal(r)=n0 Wsp(r) x(1+x)2, Wsp(r) = 1 1 + exp(r−Rsp)/∆sp,(D.3) ensuring a finite outer extent without discontinuities. The reduced convolution result ρshape(r) is normalized by the geometric factor Λgeom =Ngal M(elliptical) DM Z∞ 0 4πr2ρshape(r)dr ,(D.4) and scaled by the nonlinear reverberation gain Grev =1+Ngalη, (D.5) to yield the final effective density ρeff(r)=Grev Λgeom Sbase ρshape(r).(D.6) For the fiducial 1015 M⊙cluster, the adopted parameters are Lmean = 636 kpc, σln L= 2.5, Ngal = 300, η = 8.3×10−3, corresponding to Grev ≃ 3 . 5. The code produces three diagnostic figures: (i) cluster density and enclosed mass, (ii) circular velocity profile, and (iii) 2D distributions of ρeff and vc, reproducing the results presented in Figures 10–12. 134
1 2# ===================================================================== 3# Galaxy Cluster Static Transport-Overlap vs. NFW Benchmark 4# 5# Purpose: 6# Compute the cluster-scale effective transport density _eff(r) 7# arising from the collective overlap of mesoscopic transport halos 8# surrounding individual galaxies, and compare its emergent profile 9# with the canonical NFW dark-matter distribution. 10 # 11 # Model Overview: 12 # 1. K_tail(s) defines the dimensionless single-galaxy transport kernel, 13 # representing the exponential intensity tail of an elliptical halo. 14 # 15 # 2. The elliptical-halo microphysics determine the absolute scaling 16 # through Sbar_base, which sets the normalization of _eff(r) 17 # once the effective damping length L_damp_eff is chosen. 18 # Here, we adopt the same L_damp_eff as used in the single-elliptical 19 # simulation for internal consistency. 20 # 21 # 3. The spatial convolution is performed in a reduced (~2.5D) form 22 # for numerical stability. One geometric length factor is omitted 23 # and later restored by _geom, derived via global mass conservation. 24 # 25 # 4. A nonlinear reverberation gain factor, 26 # G_rev = 1 + N_gal * , 27 # accounts for weak cross-coherence between overlapping transport 28 # fields, where is the fractional coupling coefficient. 29 # 30 # Key Control Parameters: 31 # _lnL Logarithmic dispersion of transport lengths (heterogeneity): 32 # broadens the effective overlap and governs the curvature 33 # of _eff(r). Crucial for matching the NFW slope. 34 # 35 # Reverberation coupling coefficient: 36 # controls the degree of collective amplification among 37 # overlapping galaxy halos. Values in the range 1010 38 # yield physically plausible coherence gains. 39 # 40 # Once the galaxy microphysics and damping scale are fixed, 41 # _lnL and become the two primary tuning parameters determining 42 # the cluster-scale density and velocity profiles. 43 # 44 # Final Expression: 45 # _eff(r) = G_rev * _geom * Sbar_base * _shape(r) 46 # ===================================================================== 47 48 import numpy as np 49 import matplotlib.pyplot as plt 50 51 # --------------------------------------------------------------------- 52 # Fundamental constants and units 53 # --------------------------------------------------------------------- 135
54 G = 6.67430e-11 # m^3 kg^-1 s^-2 55 c = 2.99792458e8 # m/s 56 Msun= 1.98847e30 # kg 57 pc = 3.085677581491367e16 # m 58 kpc = 1.0e3 * pc 59 Mpc = 1.0e6 * pc 60 pi = np.pi 61 62 # --------------------------------------------------------------------- 63 # Cluster-scale parameters 64 # --------------------------------------------------------------------- 65 L_bar = 636.0 * kpc # main transport length [m] 66 sigma_lnL = 2.5 # lognormal scatter in L 67 R_cl = 2.0 * Mpc # cluster radius [m] 68 N_gal = 300 # number of elliptical galaxies 69 70 # Structural (cluster-scale) parameters 71 R_sp = 2.5 * Mpc # splashback radius [m] 72 Delta_sp = 0.5 * Mpc # splashback transition width [m] 73 r_s = 1.0 * Mpc # satellite NFW scale radius [m] 74 n0_sat = 2.0e-63 # satellite number-density scale [m^-3] 75 76 # --------------------------------------------------------------------- 77 # Elliptical-halo microphysics Sbar_base 78 # --------------------------------------------------------------------- 79 L_damp_eff = L_bar 80 c_eff_E = 2.0e5 81 Gamma_eff = c_eff_E / L_damp_eff 82 lambda_res = 27.0 * pc 83 Omega_E = 2.0 * pi * c_eff_E / lambda_res 84 Q_eff = Omega_E / (2.0 * Gamma_eff) 85 Xi_nabla = 1.0 + (2.0 * pi * L_damp_eff / lambda_res)**2 86 87 M_tot_E = 2.0e11 * Msun # total mass of one elliptical 88 m_star_E = Msun 89 N_star_E = M_tot_E / m_star_E 90 Re_E = 5.0 * kpc 91 a_h_E = Re_E / 1.8153 92 V_gal_E = 2.0 * pi * a_h_E**3 93 v_E = 200.0e3 94 alpha_E = 1.0 95 96 K_base = (G * alpha_E * v_E) / (2.0 * c**2 * c_eff_E**3) \ 97 * L_damp_eff * Q_eff * (N_star_E**2 * m_star_E**2 / V_gal_E) 98 Sbar_base = Xi_nabla * K_base 99 100 print("=== Single Elliptical Kernel Parameters ===") 101 print(f"L_bar = {L_bar/kpc:.1f} kpc, Q_eff = {Q_eff:.3e}, Xi_nabla = ,→{Xi_nabla:.3e}") 102 print(f"Sbar_base = {Sbar_base:.3e} [kg/m^3]\n") 103 104 # --------------------------------------------------------------------- 105 # Dimensionless transport kernel (shape only) 106 # --------------------------------------------------------------------- 136
107 def K_tail(s, L): 108 """Dimensionless transport kernel.""" 109 s_eps = 1.0 * pc 110 s_clp = np.maximum(s, s_eps) 111 return np.exp(-2.0 * s_clp / L) / (((4.0 * pi)**2) * s_clp**2) 112 113 def K_tail_ensemble(s, L_mean, sigma_ln): 114 """Lognormal-averaged kernel ensemble.""" 115 if sigma_ln <= 0.0: 116 return K_tail(s, L_mean) 117 xi = np.array([-2, -1, 0, 1, 2]) * sigma_ln 118 w = np.array([1, 4, 10, 4, 1], dtype=float) 119 w /= w.sum() 120 out = np.zeros_like(s) 121 for wi, Li in zip(w, L_mean * np.exp(xi)): 122 out += wi * K_tail(s, Li) 123 return out 124 125 Kfun = lambda s: K_tail_ensemble(s, L_bar, sigma_lnL) 126 127 # --------------------------------------------------------------------- 128 # Galaxy number-density profile (satellite population) 129 # --------------------------------------------------------------------- 130 def W_splash(r): 131 """Smooth splashback taper.""" 132 return 1.0 / (1.0 + np.exp((r - R_sp) / Delta_sp)) 133 134 def n_sat_NFW(r): 135 """Satellite population ~ 1/[x(1+x)^2].""" 136 x = np.maximum(r / r_s, 1e-9) 137 return W_splash(r) * n0_sat / (x * (1.0 + x)**2) 138 139 # --------------------------------------------------------------------- 140 # Reduced spherical convolution (~2.5D form) 141 # --------------------------------------------------------------------- 142 def angular_average_K(r, rp, Kfun): 143 """Angularly averaged kernel integral K(s) d.""" 144 mu, w = np.polynomial.legendre.leggauss(48) 145 r, rp = np.atleast_1d(r).reshape(-1,1), np.atleast_1d(rp).reshape(1,-1) 146 mu3, w3 = mu[:,None,None], w[:,None,None] 147 s = np.sqrt(r**2 + rp**2 - 2.0 * r * rp * mu3) 148 return np.sum(w3 * Kfun(s), axis=0) 149 150 def rho_from_population(n_pop): 151 """Compute reduced convolution _shape(r).""" 152 I = angular_average_K(r, r, Kfun) 153 integrand = n_pop * r 154 return 2.0 * pi * (I @ (integrand * dr)) 155 156 # --------------------------------------------------------------------- 157 # Radial grid and convolution 158 # --------------------------------------------------------------------- 159 r = np.geomspace(0.2 * kpc, 3.0 * Mpc, 480) 160 dr = np.gradient(r) 137
[20] Ortwin Gerhard, Andi Kronawitter, RP Saglia, and Ralf Bender. Dynamical family properties and dark halo scaling relations of giant elliptical galaxies. The Astronomical Journal, 121(4):1936, 2001. [21] Julio F Navarro. The structure of cold dark matter halos. In Symposiuminternational astronomical union, volume 171, pages 255–258. Cambridge University Press, 1996. [22] Julio F Navarro, Carlos S Frenk, and Simon DM White. A universal density profile from hierarchical clustering. The Astrophysical Journal, 490(2):493, 1997. [23] Surhud More, Benedikt Diemer, and Andrey V Kravtsov. The splashback radius as a physical halo boundary and the growth of halo mass. The Astrophysical Journal, 810(1):36, 2015. [24] Sarah M Hansen, Timothy A McKay, Risa H Wechsler, James Annis, Erin Scott Sheldon, and Amy Kimball. Measurement of galaxy cluster sizes, radial profiles, and luminosity functions from sdss photometric data. The Astrophysical Journal, 633(1):122, 2005. 144