Chemical master equation parameter exploration using DMRG

Abstract

Well-mixed chemical reaction networks (CRNs) contain many distinct chemical species with copy numbers that fluctuate in correlated ways. While those correlations are typically monitored via Monte Carlo sampling of stochastic trajectories, there is interest in systematically approximating the joint distribution over the exponentially large number of possible microstates using tensor networks or tensor trains. We exploit the tensor network strategy to determine when the steady state of a seven-species gene toggle switch CRN model supports bistability as a function of two decomposition rates, both parameters of the kinetic model. We highlight how the tensor network solution captures the effects of stochastic fluctuations, going beyond mean field and indeed deviating meaningfully from a mean-field analysis. The work furthermore develops and demonstrates several technical advances that will allow steady-states of broad classes of CRNs to be computed in a manner conducive to parameter exploration. We show that the steady-state distributions can be computed via the ordinary density matrix renormalization group (DMRG) algorithm despite having a non-Hermitian rate operator with a small spectral gap, we illustrate how that steady-state distribution can be efficiently projected to an order parameter that identifies bimodality, and we employ excited-state DMRG to calculate a relaxation timescale for the bistability.

References in corpus (56)

• The density-matrix renormalization group in the age of matrix product states
• Entanglement in quantum critical phenomena
• A Practical Introduction to Tensor Networks: Matrix Product States and Projected Entangled Pair States
• Matrix Product Density Operators: Simulation of finite-T and dissipative systems
• The ITensor Software Library for Tensor Network Calculations
• Unifying time evolution and optimization with matrix product states
• Matrix product states represent ground states faithfully
• Tensor networks for complex quantum systems
• Ab Initio Quantum Chemistry using the Density Matrix Renormalization Group
• Sampling rare switching events in biochemical networks
• Forward Flux Sampling for rare event simulations
• Simulating rare events in equilibrium or non-equilibrium stochastic systems
• Entropy, entanglement, and area: analytical results for harmonic lattice systems
• Studying Two Dimensional Systems With the Density Matrix Renormalization Group
• The density matrix renormalization group for ab initio quantum chemistry
• Stochastic gene expression as a many body problem
• A spin-adapted Density Matrix Renormalization Group algorithm for quantum chemistry
• Alternating minimal energy methods for linear systems in higher dimensions. Part I: SPD systems
• Alternating minimal energy methods for linear systems in higher dimensions. Part II: Faster algorithm and application to nonsymmetric systems
• Density matrix renormalization group algorithms with a single center site
• An entanglement-area law for general bosonic harmonic lattice systems
• Variational Matrix Product Operators for the Steady State of Dissipative Quantum Systems
• A Strictly Single-Site DMRG Algorithm with Subspace Expansion
• Efficient Tree Tensor Network States (TTNS) for Quantum Chemistry: Generalizations of the Density Matrix Renormalization Group Algorithm
• Generic Construction of Efficient Matrix Product Operators
• Efficient tensor network simulation of IBM's Eagle kicked Ising experiment
• Wilson renormalization of a reaction-diffusion process
• Skin superfluid, topological Mott insulators, and asymmetric dynamics in interacting non-Hermitian Aubry-Andre-Harper models
• Master equations and the theory of stochastic path integrals
• Markovian Dynamics on Complex Reaction Networks
• Enhancement of the stability of genetic switches by overlapping upstream regulatory domains
• Gauge fixing, canonical forms and optimal truncations in tensor networks with closed loops
• A stochastic spectral analysis of transcriptional regulatory cascades
• Noise-induced multistability in chemical systems: Discrete vs Continuum modeling
• Dynamical phase behavior of the single- and multi-lane asymmetric simple exclusion process via matrix product states
• Classical stochastic dynamics and continuous matrix product states: gauge transformations, conditioned and driven processes, and equivalence of trajectory ensembles
• Conversion of projected entangled pair states into a canonical form
• A variational approach to the stochastic aspects of cellular signal transduction
• Dynamics and large deviation transitions of the XOR-Fredrickson-Andersen kinetically constrained model
• Gauging tensor networks with belief propagation
• Non-Hermitian Fermi-Dirac Distribution in Persistent Current Transport
• Tensor Networks contraction and the Belief Propagation algorithm
• Optimal sampling of dynamical large deviations via matrix product states
• Preserving Correlations Between Trajectories for Efficient Path Sampling
• Optimal sampling of dynamical large deviations in two dimensions via tensor networks
• Tensor-train approximation of the chemical master equation and its application for parameter inference
• Analysis of a density matrix renormalization group approach for transport in open quantum systems
• One-parameter extension of the Doi-Peliti formalism and relation with orthogonal polynomials
• Quantifying Rare Events in Stochastic Reaction-Diffusion Dynamics Using Tensor Networks
• Mean-field vs. stochastic models for transcriptional regulation
• Tensor product approach to modelling epidemics on networks
• Biorthonormal Matrix-Product-State Analysis for Non-Hermitian Transfer-Matrix Renormalization-Group in the Thermodynamic Limit
• Loop Series Expansions for Tensor Networks
• Density-matrix renormalization group algorithm for non-Hermitian systems
• Variational Matrix Product State Approach for Non-Hermitian System Based on a Companion Hermitian Hamiltonian
• From high-dimensional committors to reactive insights

Cited by in corpus (1)

• Self-interacting processes via Doob conditioning