Quantifying Rare Events in Stochastic Reaction-Diffusion Dynamics Using Tensor Networks

Abstract

The interplay between stochastic chemical reactions and diffusion can generate rich spatiotemporal patterns. While the timescale for individual reaction or diffusion events may be very fast, the timescales for organization can be much longer. That separation of timescales makes it particularly challenging to anticipate how the rapid microscopic dynamics gives rise to macroscopic rates in the non-equilibrium dynamics of many reacting and diffusing chemical species. Within the regime of stochastic fluctuations, the standard approach is to employ Monte Carlo sampling to simulate realizations of random trajectories. Here, we present an alternative numerically tractable approach to extract macroscopic rates from the full ensemble evolution of many-body reaction diffusion problems. The approach leverages the Doi-Peliti second-quantized representation of reaction-diffusion master equations along with compression and evolution algorithms from tensor networks. By focusing on a Schlögl model with one-dimensional diffusion between $L$ otherwise well-mixed sites, we illustrate the potential of the tensor network approach to compute rates from many-body systems, here with approximately $3\times 10^{15}$ microstates. Specifically, we compute the rate for switching between metastable macrostates, with the expense for computing those rates growing subexponentially in $L$. Because we directly work with ensemble evolutions, we crucially bypass many of the difficulties encountered by rare event sampling techniques$\unicode{x2013}$detailed balance and reaction coordinates are not needed.

References in corpus (13)

• The density-matrix renormalization group in the age of matrix product states
• Dynamic first-order phase transition in kinetically constrained models of glasses
• First-order dynamical phase transition in models of glasses: an approach based on ensembles of histories
• A Strictly Single-Site DMRG Algorithm with Subspace Expansion
• Studying Rare Events using Forward-Flux Sampling: Recent Breakthroughs and Future Outlook
• Engineering reaction-diffusion networks with properties of neural tissue
• A variational approach to the stochastic aspects of cellular signal transduction
• Finite time large deviations via matrix product states
• A probabilistic framework for particle-based reaction-diffusion dynamics using classical Fock space representations
• Preserving Correlations Between Trajectories for Efficient Path Sampling
• Dynamical mean-field theory: from ecosystems to reaction networks
• Using tensor network states for multi-particle Brownian ratchets
• Computing time-periodic steady-state currents via the time evolution of tensor network states

Cited by in corpus (3)

• Chemical master equation parameter exploration using DMRG
• Redundant basis interpretation of Doi-Peliti method and an application
• From high-dimensional committors to reactive insights