A General Path-Integral Monte Carlo Method for Exact Simulations of Chemical Reaction Networks.

Chemical Reaction Networks (CRNs) are stochastic many-body systems used to model real-world chemical systems through a differential Master Equation (ME); analytical solutions to these equations are only known for the simplest systems. In this paper, we construct a path-integral inspirited framework for studying CRNs. Under this scheme, the time-evolution of a reaction network can be encoded in a Hamiltonian-like operator.