Skip to main content
Article
100%
HD

Imagine teaching a computer to predict the evolution of a cell or the spread of a virus. These systems are "stochastic": their future depends on random events, like rolling a die. Scientists use Gillespie-type algorithms to simulate this, picking the next event based on precise probabilities. But if we want to use these simulations to train an AI with deep learning, a problem arises: the rigid, categorical choice of the event (the "die roll") blocks gradient flow, stopping the network from learning.

A new arXiv study proposes an elegant solution called **Propensity Straight-Through (PST)**. Instead of "softening" randomness with tricks like Gumbel-Softmax (which adds noise and needs tuning), PST differentiates the exact conditional mean of the next step directly. It calculates the exact sensitivity by looking at normalized propensitiesβ€”the reaction probabilitiesβ€”without altering the simulation.

The key lies in the affine state update of continuous-time Markov chains: PST exploits this linearity to get exact per-step gradients, proving that multi-step accumulated error vanishes when downstream dependencies are linear. Results are striking: on benchmarks like genetic oscillators or ion channels, PST matches relaxed methods' accuracy but converges up to 3x faster. The ultimate test? They trained a 200,000+ parameter network to classify MNIST digits using exact "hard" sampling, hitting 98.22% accuracy. PST paves the way for scalable deep learning through exact stochastic trajectoriesβ€”no temperature, no Gumbel samples. A breakthrough for modeling biology, chemistry, and more with AI.

Bonus
100%
HD

The paper addresses the bottleneck of integrating Continuous-Time Markov Chains (CTMCs) with gradient-based optimization. The Gillespie algorithm (SSA) samples discrete events $k$ with probability proportional to propensities $a_k(x)$. This "hard" sampling makes the trajectory non-differentiable w.r.t. parameters $\theta$ governing propensities.

PST solves this by differentiating the exact conditional mean of the state update $\mathbb{E}[\Delta x | x] = \sum_k \nu_k \frac{a_k(x)}{a_0(x)}$, where $\nu_k$ are state-change vectors. Since the update is affine in the state, this mean is a differentiable rational function of propensities. The straight-through gradient uses this exact sensitivity in the backward pass while keeping hard sampling in the forward pass.

Theoretically, the authors derive the per-step discrepancy between the composed one-step sensitivity and the true multi-step sensitivity $\frac{d \mathbb{E}[x_T]}{d\theta}$. They prove this discrepancy vanishes identically if the downstream loss depends affinely on the state, ensuring unbiased gradients in many practical cases (e.g., linear observables).

Empirically, PST matches Gumbel-Softmax Straight-Through (GS-ST) accuracy on reversible dimerization (0.06% error), genetic oscillator (1.7% error), 50-task repressilator suite (0.17% median error), and patch-clamp recordings ($R^2=0.988$). Crucially, PST converges 3.0x faster on the oscillator and 2.1x faster on ion channels due to the absence of temperature and Gumbel relaxation variance. On MNIST, a 203,796-parameter Stochastic Reaction Network (SRN) trained with PST reaches 98.22% accuracy, demonstrating deep-learning scalability. PST thus offers a temperature-free, Gumbel-free path to scalable gradient-based learning through exact stochastic trajectories.

Ask AI

Explore key aspects of this story with 3 quick insights

Source

Authors:
Jose M. G. VilarLeonor Saiz
Date Published:
Original Article arxiv.org
Fact-check verified on 29/10/2026 Primary Study / Original Source
Read our Editorial Guidelines & AI Ethics Policy
Creative Commons License CC BY-NC-SA 4.0

Keep exploring

Discover the full archive of science bites or head back to home

Riconnessione al server... Aggiornamento della sessione in corso
Connessione interrotta Impossibile stabilire il collegamento
Sessione aggiornata Ricarica per visualizzare i contenuti