Before recombination, photons and baryons formed a tightly coupled plasma. Sound waves propagated through this plasma, driven by the competition between radiation pressure and gravity. At recombination, the photons decoupled and these sound waves were frozen — preserved as the acoustic peak pattern we observe today in the CMB angular power spectrum, and as the preferred clustering scale of approximately 490 million light years known as the Baryon Acoustic Oscillation.
Both signatures are treated as fossil evidence of a one-time event. The standard model asserts they could only have been produced this way.
The BFUT challenge is precise and limited. This paper does not deny that the observations are real. It denies that they uniquely require a singular Big Bang origin. The mechanism producing acoustic oscillations — Thomson scattering coupling of photons and baryons against diffusion-like damping — operates continuously wherever ionised plasma exists. It is not unique to the pre-recombination era. And the BAO scale of approximately 490 million light years is not a raw observation: it is a number derived by assuming the very expansion history it is cited to support.
The Big Flare-Up Theory proposes that ongoing astrophysical processes — expanding HII shells around star-forming regions, AGN-driven cavities in cluster plasmas, supernova superbubbles, ripple-like pressure disturbances in intracluster media — continuously inject shell-like and ripple-like acoustic power into the baryon-photon system. This ongoing injection is balanced by scale-dependent damping. The result is a steady-state power spectrum at every scale.
The master equation governing this balance is:
dP(k,t)/dt = I(k,t) − D(k,t) × P(k,t)
In statistical steady state this gives P(k) = I(k) / D(k), where I(k) is the ongoing injection power at wavenumber k from shell-like source populations, and D(k) = D₀ + D₂k² is the effective damping law with the k² term physically motivated by the Silk-like diffusion coefficient D_diff = c/(3 n_e σ_T).
Shell geometry alone generates oscillatory Fourier structure: a shell of radius R produces a power contribution proportional to [sin(kR)/(kR)]². This is not a property of a Big Bang. It is a property of any shell-like structure in three-dimensional space.
To test whether this mechanism is computationally viable — not merely a theoretical analogy — a numerical simulation programme was conducted across six progressively refined versions. The simulation implements the master equation through discrete time-stepping toward steady state. Shell source populations are modelled as Gaussian-weighted distributions of shell radii using compensated finite-thickness sinc-squared kernels. The real-space correlation function is obtained by numerical Fourier transform. The angular power spectrum is obtained by a Limber-like projection.
The question asked of each simulation version was simple: does the mechanism produce recognisable CMB-like acoustic peaks and a BAO-like real-space correlation bump? And does the result improve systematically as the physics becomes more realistic?
| Version | Peaks | Peak Positions (ℓ) | BAO-like bump | Status |
|---|---|---|---|---|
| v2 | 1 | ~210 | ~50 Mpc | First directional signal |
| v3-lite | 2 | ~210, ~387 | ~80.6 Mpc | First encouraging result |
| v4-fast | 2 | 222, 415 | ~100 Mpc | Clear convergence |
| v5-fast | 3 | 222, 414, 602 | ~120 Mpc | ✓ First major breakthrough |
| v6-fast | 3 | 218, 411, 595 | ~135 Mpc | ✓ Strong proof-of-principle |
The monotonic improvement is the scientific result. Every refinement of the physics — larger source scale, increased compensation offset, broader projection depth, two-band source mixture — moved the simulation closer to the observed values. This is the signature of a mechanism that is physically real rather than a curve-fitting exercise. A model that improves consistently under more realistic physics is making a genuine scientific claim.
| Parameter | Value | Physical meaning |
|---|---|---|
| Source band 1 | R* = 300 Mpc, σ = 16 Mpc, weight 0.70 | Dominant shell scale and spread |
| Source band 2 | R* = 420 Mpc, σ = 30 Mpc, weight 0.30 | Secondary shell population |
| Compensation offset ΔR | 42 Mpc | Finite shell thickness correction |
| Compensation strength β | 0.98 | Inner shell subtraction weight |
| Shell thickness σ_shell | 10 Mpc | Physical shell wall width |
| Damping D₂ | 1.5 | Silk-type diffusion coefficient |
| Damping D₄ | 0 | Higher-order suppression (unused) |
| Projection centre χ* | 16,500 Mpc | Effective comoving depth |
| Projection width σ_χ | 2,600 Mpc | Depth window width |