Morphogen-Controlled Topology-Preserving Sparse Growth
Implementation & benchmark of arXiv:2607.24619 — Control of morphology and topology in a lattice model of branching morphogenesis
Source paper: Control of morphology and topology in a lattice model of branching morphogenesis arXiv:2607.24619 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Replace fixed sparse masks with a stochastic birth-death process for neural connections or spatial units. A diffusing morphogen-like utility field controls where connections are added or removed, while a local simple-point test rejects removals or additions that would disconnect a layer or alter a prescribed computational topology. This creates an adaptive sparse architecture with a tunable compact-to-branched transition rather than unconstrained magnitude pruning.
Formulas
Mathematical statement
Let $n_r\in\{0,1\}$ denote occupancy of lattice site or candidate connection $r$, let $c_r\geq 0$ be its morphogen or utility concentration, and let $e$ range over local neighbors. The paper uses local growth and degrowth rates $\nu_r^+$ and $\nu_r^-$: $\nu_r^+=\bar\nu^+\sum_e f_r(1-n_r)n_{r+e}\bar\partial_e^+c_r$ and $\nu_r^-=\bar\nu^-\sum_e f_r n_r(1-n_{r+e})\bar\partial_e^-c_r$. Here $f_r$ is a local control or topology factor, $\bar\nu^+$ and $\bar\nu^-$ are base rates, and $\bar\partial_e^\pm c_r$ are nonnegative directional morphogen-response operators. In the neural adaptation, $r$ indexes candidate edges, $n_r$ is the binary mask, $c_r$ is a smoothed estimate of edge utility, and $f_r=0$ whenever a proposed local update violates connectivity or a prescribed topology constraint. A continuous-time approximation is $\dot p_r=(1-p_r)\nu_r^+-p_r\nu_r^-$ for edge occupancy probabilities $p_r\in[0,1]$. The key mechanism is multiplicative gating by both local neighborhood occupancy and the morphogen gradient: growth occurs near active structure and high utility, while pruning is suppressed when it would break topology.
Implementation notes
Integrate the module into a spatial CNN, graph neural network, or Transformer linear layer represented by a sparse bipartite graph. Maintain a binary mask $n_r$ for candidate edges and a scalar field $c_r$ on edges or nodes. Every $K$ optimizer steps, compute edge utility $s_r$ from the exponential moving average of $|w_r\,\partial L/\partial w_r|$, diffuse it with $c\leftarrow c+(\Delta t/\tau_c)(s-c+D_c\Delta_{graph}c)$, and estimate nonnegative directional responses using $\bar\partial_e^+c_r=\max(c_{r+e}-c_r,0)$ and $\bar\partial_e^-c_r=\max(c_r-c_{r+e},0)$. For each minibatch of proposals, sample additions with probability $1-\exp(-\Delta t\nu_r^+)$ and removals with probability $1-\exp(-\Delta t\nu_r^-)$. Before accepting a removal, temporarily delete the edge and run a cheap local connectivity or simple-point test; set $f_r=0$ if an output unit loses all inputs, an input loses all outputs, or the graph's monitored connected-component count changes. Keep the parameter budget fixed by pairing accepted additions and removals. The paper supplies the birth-death construction and topology-preserving local gating; the utility source, diffusion discretization, and graph topology test are engineering adaptations. First run a CIFAR-10 experiment with a 90%-sparse 3-layer CNN, comparing magnitude pruning, RigL-style rewiring, and this method at equal FLOPs. Sweep $\rho=\bar\nu^+/\bar\nu^-$. Measure graph radius of gyration versus active-edge count, validation accuracy, and connected components. The predicted signature is a compact-to-branched transition as $\rho D_c$ increases, while the topology-preserving regime keeps component count constant and yields $R_g\propto N^\alpha$ with $\alpha$ approaching approximately $0.67$ for sufficiently large clusters. A failure of this exponent, absence of a reproducible transition, or a large accuracy loss at fixed FLOPs falsifies the transfer.
Verification
This idea has not been verified yet.
Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.
Artifacts
Artifacts unavailable.