We introduce Hebbian-Oscillatory Co-Learning (HOC-L), a unified two-timescale dynamical framework for joint structural plasticity and phase synchronization in bio-inspired sparse neural architectures. HOC-L couples two recent frameworks: the hyperbolic sparse geometry of Resonant Sparse Geometry Networks (RSGN), which employs Poincaré ball embeddings with Hebbian-driven dynamic sparsity, and the oscillator-based attention of Selective Synchronization Attention (SSA), which replaces dot-product attention with Kuramoto-type phase-locking dynamics. The key mechanism is synchronization-gated plasticity: the macroscopic order parameter r(t) of the oscillator ensemble gates Hebbian structural updates, so that connectivity consolidation occurs only when sufficient phase coherence signals a meaningful computational pattern. We prove convergence of the joint system to a stable equilibrium via a composite Lyapunov function and derive explicit timescale separation bounds. The resulting architecture achieves O(n⋅k) complexity with k≪n, preserving the sparsity of both parent frameworks. Numerical simulations confirm the theoretical predictions, demonstrating emergent cluster-aligned connectivity and monotonic Lyapunov decrease.
Figures & tables
Figure 1: Overall HOC-L architecture. The left stream (blue) implements the RSGN [ 18 ] structural pathway: inputs are embedded in the Poincaré ball, sparse hyperbolic neighborhoods are constructed, and Hebbian plasticity modifies structural connectivity. The right stream (orange) implements the SSA [ 21 ] oscillatory pathway: tokens receive oscillator encodings, phase-locking attention computes synchronization-based weights, and the macroscopic order parameter r(t) is extracted. The central coupling bridge (green) implements synchronization-gated plasticity: the order parameter gates Hebbian updates, and structural changes modulate future synchronization.
Figure 2: Simulation results: Two-timescale dynamics in HOC-L. A Kuramoto oscillator system ( N=50 , K=2.0 ) is simulated jointly with synchronization-gated Hebbian weight updates for 1000 time steps ( dt=0.05 ), using the smooth sigmoid gate G(r)=σ(β(r−rc)) with β=20 . The fast curve (red) shows the running average of SSA phase update magnitudes ∣dθi/dt∣ , exhibiting rapid oscillations as oscillators seek synchronization. The slow curve (blue) shows the running average of Hebbian weight update magnitudes ∣ΔWij∣ , which are continuously modulated by G(r) : negligible when r≪rc and substantial once r exceeds rc . The vertical dashed line (green) marks the moment when G(r) crosses 0.5 (i.e., r crosses rc=0.5 ). The orange curve (right axis) traces the macroscopic order parameter r(t) . The timescale separation τfast≪τslow (Section III ), where τfast is the characteristic timescale of phase updates and τslow is that of Hebbian structural updates, is clearly visible in the distinct frequency content of the two signals.
Figure 3: Simulation results: Phase synchronization to structural plasticity coupling mechanism in a Kuramoto system with N=8 oscillators ( K=3.0 , 2000 time steps). (A) Final oscillator phases on the unit circle: orange markers denote phase-locked (synchronized) oscillators that have formed a coherent cluster, while grey markers denote desynchronized oscillators with distinct intrinsic frequencies. (B) Order parameter r(t) trajectory over time, with the critical threshold rc=0.5 shown as a dashed line; the shaded region indicates intervals where the plasticity gate is open ( r>rc ). (C) Emergent weight matrix Wij after simulation: strong positive weights (warm colors) develop between synchronized oscillators, while connections involving desynchronized oscillators remain weak or near zero, demonstrating that synchronization-gated Hebbian plasticity produces sparse, cluster-aligned connectivity.
Figure 4: Simulation results: Convergence basin of the HOC-L Lyapunov function. The 3D surface shows the composite Lyapunov function V(W,θ)=−2NK∑i,jcos(θi−θj)+2λ∥W∥F2 evaluated on a grid over the projected weight norm ( ∥W∥F -axis) and mean phase deviation ( ϕˉ -axis). Colored trajectories trace the joint dynamics of 8 random initializations through the (∥W∥F,ϕˉ,V) space, each computed from actual numerical integration of the coupled Kuramoto–Hebbian ODE system ( N=20 oscillators, 600 steps per trajectory). All trajectories converge toward the basin minimum (red star) where oscillators are synchronized ( ϕˉ≈0 ) and weights are regularized, confirming the monotonic decrease of V along system trajectories as guaranteed by Theorem 3 . Contour lines at the base reveal the basin geometry projected onto the (∥W∥F,ϕˉ) plane.