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 (33)

• The density-matrix renormalization group in the age of matrix product states
• Efficient classical simulation of slightly entangled quantum computations
• A Practical Introduction to Tensor Networks: Matrix Product States and Projected Entangled Pair States
• The ITensor Software Library for Tensor Network Calculations
• Time-dependent variational principle for quantum lattices
• Unifying time evolution and optimization with matrix product states
• Time-evolution methods for matrix-product states
• Mechanisms of noise-resistance in genetic oscillators
• Sampling rare switching events in biochemical networks
• Dynamic first-order phase transition in kinetically constrained models of glasses
• Forward Flux Sampling for rare event simulations
• Simulating rare events in equilibrium or non-equilibrium stochastic systems
• First-order dynamical phase transition in models of glasses: an approach based on ensembles of histories
• Post-Matrix Product State Methods: To tangent space and beyond
• Tangent-space methods for uniform matrix product states
• Time integration of tensor trains
• A Strictly Single-Site DMRG Algorithm with Subspace Expansion
• Studying Rare Events using Forward-Flux Sampling: Recent Breakthroughs and Future Outlook
• Enhancement of the stability of genetic switches by overlapping upstream regulatory domains
• Stochastic Simulation of Reaction-Diffusion Systems: A Fluctuating-Hydrodynamics Approach
• Dynamical phase behavior of the single- and multi-lane asymmetric simple exclusion process via matrix product states
• Engineering reaction-diffusion networks with properties of neural tissue
• A variational approach to the stochastic aspects of cellular signal transduction
• Asymptotic form of the approach to equilibrium in reversible recombination reactions
• Non-meanfield deterministic limits in chemical reaction kinetics far from equilibrium
• 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
• Field Theory of free Active Ornstein-Uhlenbeck Particles
• Dynamical mean-field theory: from ecosystems to reaction networks
• Using tensor network states for multi-particle Brownian ratchets
• Algebraic probability, classical stochastic processes, and counting statistics
• 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