Stochastic Evolution of Galactic Star Formation with Halo Coupling, AGN Quenching and Hopf Bifurcation Dynamics

Abstract

We present a computational framework for galactic evolution based on a coupled stochastic nonlinear oscillator, implemented with the \textbf{Stochastic Hopf Engine}. Gas density ($G$) and star formation rate ($S$) co-evolve through a supercritical Hopf bifurcation, capturing the transition from quiescent stability to merger-driven starbursts. Scatter in dark matter halo properties, modeled as multiplicative noise via the \textbf{Euler--Maruyama method}, broadens the bifurcation into a regime where noise-induced bursts occur below the deterministic threshold. Simulations reveal a periodic signature, the \textbf{Galactic Heartbeat}, emerging as a deterministic limit cycle validated by the \textbf{data3} resonance peak in the star-formation spectrum. A radial reduction yields an effective \textbf{Fokker--Planck equation} for burst amplitude; its stationary solution matches numerical PDFs, providing statistical closure. Including differential shear $Ω(r)$ and spatially varying bifurcation fields reproduces spiral morphologies and AGN-driven quenching. Driving the growth parameter sub-critical ($r_{agn} < 0$) yields ``Red and Dead'' cores via attractor collapse. Dark matter halo scatter suppresses mean star formation while enhancing intermittency, offering a minimal yet interpretable framework linking local feedback and global potentials to macroscopic galactic evolution.

Figures (17)

• Figure 1: Noise-induced starburst dynamics driven by dark matter halo scatter. Left: Phase-space trajectory showing a noise-broadened Hopf structure. Middle: Time series of the starburst amplitude revealing intermittent bursts. Right: Stationary amplitude distribution compared with the analytical Fokker--Planck prediction, demonstrating excellent agreement.
• Figure 2: Noise-induced bifurcation diagram of starburst amplitude as a function of halo scatter $\sigma_k$. Increasing scatter leads to a continuous broadening of amplitude distributions, indicating a stochastic transition from quiescent to burst-dominated star formation.
• Figure 3: Noise-induced starburst cycles driven by dark matter halo scatter.Top row: (a) Stochastic bifurcation diagram showing the radial amplitude $A=\sqrt{x^2+y^2}$ as a function of halo scatter $\sigma_k$. The deterministic Hopf point at $r=\bar{k}_{\rm DM}$ broadens into a finite-width bifurcation band, indicating noise-induced transitions. (b) Mean starburst amplitude $\langle A\rangle$ versus $\sigma_k$, demonstrating suppression of coherent oscillations with increasing halo variability. (c) Variance of the amplitude, $\mathrm{Var}(A)$, revealing enhanced burstiness and intermittency at large $\sigma_k$. Bottom row: Phase-space trajectories $(x,y)$ for increasing feedback strength $r$, illustrating quiescent regimes ($r\le\bar{k}_{\rm DM}$) and noise-sustained starburst cycles for $r>\bar{k}_{\rm DM}$. These results validate the stochastic Hopf framework and demonstrate that halo mass--concentration scatter alone can induce bursty star formation even in otherwise marginally stable systems.
• Figure 4: Multi-panel diagnostic of galactic evolution. Top: SFR history with marked merger window. Middle: Phase space transition from stable focus to limit cycle. Bottom: Spectral analysis identifying the characteristic pulse (data3) at 0.157 Hz.
• Figure 5: Diagnostic of Galactic Death. Top: SFR light curve showing the collapse at $t=300$. Middle: Phase space showing the limit cycle collapsing into the "Red and Dead" sink. Bottom: Spectral analysis confirming the disappearance of the galactic heartbeat.