Path integrals for stochastic hybrid reaction-diffusion processes Paul C. Bressloff Department of Mathematics, University of Utah, Salt Lake City, UT, USA We construct a functional path integral for a stochastic hybrid reaction-diffusion (RD) equation, in which the reaction term depends on the discrete state of a randomly switching environment. We proceed by spatially discretizing the RD system and using operator methods and coherent spin states to derive a path integral representation of the lattice model. The path integral specifies the distribution of trajectories in a state-space consisting of the set of local concentrations and the environmental state. Taking the continuum limit yields a path integral for the stochastic hybrid RD equation. We show how in the semi-classical limit, the stochastic dynamics reduces to a system of Langevin equations. The latter consists of an RD equation coupled to a stochastic probability function z that evolved according to an Ornstein-Uhlenbeck process with multiplicative noise. On eliminating the variable z, we obtain an RD equation driven by nonlocal spatiotemporal noise. This is used to determine leading order corrections to the mean-field RD equation. We then apply the theory to the example of diffusion with stochastically gated adsorption under dimerization. Com- parison with the classical problem of diffusion-limited pair-annihilation suggests that the adsorption model exhibits anomalous large-time behavior in one and two spatial dimensions. This motivates an alternative path integral formulation based on an RD master equation that treats diffusion as a set of hopping reactions and incorporates the effects of a randomly switching environment. A small noise expansion now yields a stochastic RD equation with both real and complex noise sources. Finally, diagrammatic and renormalization group methods are used to determine the anomalous large-time behavior of the mean concentration at the upper critical dimension dc = 2. I. INTRODUCTION in the weak noise (adiabatic) limit. In particular, we have shown how to derive a hybrid path integral using two alternative methods: (i) integral representations of the Diffusion processes in randomly switching environ- Dirac delta function [32, 33], which is analogous to the ments are finding increasing applications in cell biology construction of path integrals for stochastic differential and biophysics. Examples include diffusion in domains equations (SDEs) [34{36]; (ii) bra-kets and \quantum- with stochastically-gated boundaries [1{6], diffusion over mechanical" operators [37], similar in spirit to the Doi- fluctuating barriers [7{10], and stochastic gap junctions Peliti formalism for master equations [38{41]. In both [11]. Mathematically speaking, diffusion in a randomly cases the Hamiltonian of the resulting action functional switching environment is an infinite-dimensional version corresponds to the principal eigenvalue of a linear opera- of a so-called stochastic hybrid system. Stochastic hybrid tor, which combines the generator of the discrete Markov systems involve a coupling between a discrete Markov process and the vector fields of the piecewise determin- chain and a continuous stochastic process [12]. If the lat- istic dynamics. This is consistent with more rigorous re- ter evolves deterministically between jumps in the dis- sults obtained using large deviation theory [42{45]. Re- crete state, then the system reduces to a piecewise de- cently, we have derived a different hybrid path integral terministic Markov process (PDMP) [13]. Well known representation that holds outside the weak noise regime examples of finite-dimensional hybrid systems include [46]. Rather than using the \bosonic" annihilation and stochastic gene expression [14{19], voltage fluctuations in creation operators of the Doi-Peliti formalism, we de- neurons [20{27], and motor-driven intracellular transport composed the discrete states in terms of coherent \spin" [28{31]. One method for analyzing the diffusion equation states [47, 48]. Such a decomposition has recently been with switching boundary conditions is to discretize space used to study stochastic gene expression in the presence and construct the Chapman-Kolmogorov (CK) equation of promoter noise and low protein copy numbers [49{51]. for the resulting finite-dimensional stochastic hybrid sys- For example, one can effectively map the stochastic dy- tem [5]. The CK equation can then be used to derive mo- namics of a single genetic switch to a quantum spin-boson ment equations for the stochastic concentration. In the system. continuum limit, this yields a hierarchy of moment equa- arXiv:2103.07759v1 [cond-mat.stat-mech] 13 Mar 2021 tions, with the equation at r-th order taking the form In this paper, we use coherent spin states to construct a of an r-dimensional parabolic partial differential equa- path-integral formulation of a stochastic hybrid reaction- tion (PDE) that couples to lower order moments at the diffusion (RD) equation, under the assumption that the boundaries. Although the diffusing particles are non- switching environment affects the reaction term rather interacting, statistical correlations arise at the popula- than the boundary conditions. A simple example is diffu- tion level due to the fact that they all move in the same sion with stochastically-gated adsorption, where the ad- randomly switching environment. sorption rate depends on the current state of the environ- We have previously developed path integral methods ment [52]. We proceed by first constructing a path in- for studying finite-dimensional stochastic hybrid systems tegral representation of the spatially discretized RD sys- 2 tem and then taking a continuum limit. We show how a Note that U(x; t) is a stochastic concentration field. In semi-classical approximation of the resulting continuum the case of a finite interval, Ω = [0;L], the PDE has to path integral action yields the action for a stochastic RD be supplemented by boundary conditions at x = 0;L. equation with real, spatially colored noise. The latter is In previous work [5], we took the boundary conditions used to determine leading order corrections to the mean- themselves to depend on the environmental state (with field RD equation. We illustrate the theory by consid- Fn = 0). That is, ering the example of diffusion with stochastically gated adsorption. Comparison with the classical problem of bnU(0; t) + cn@xU(0; t) = dn; (2.3a) diffusion-limited pair-annihilation suggests that the ad- 0 0 0 bnU(L; t) + cn@xU(L; t) = dn (2.3b) sorption model exhibits anomalous large-time behavior in one and two spatial dimensions [53]. Therefore, we 0 0 for N(t) = n and constant coefficients bn, bn, cn, cn, derive an alternative path-integral representation based 0 dn, dn. We showed that one way to analyze the effects on an RD master equation that includes the effects of a of a random environment is to spatially discretize the switching environment, and treats diffusion as a set of PDE so that it is converted to a stochastic hybrid ODE hopping reactions that supplement the reaction associ- [5]. Introducing the lattice spacing h and setting uj = ated with absorption. A small noise expansion now yields u(jh), j 2 J = f0; 1;:::; N g with N a = L, leads to the a stochastic RD equation with both real and complex piecewise deterministic ODE noise sources. Finally, we use diagrammatic and renor- malization group (RG) methods in statistical field theory dUi X = ∆n U + F (U ); i 2 J ; (2.4) [54{58] to determine the anomalous large-time behavior dt ij j n i of the mean concentration at the upper critical dimension j2J d = 2. c n where U = (Uj; j 2 J ) and ∆ij is the discrete Laplacian The structure of the paper is as follows. In Sect. II we n introduce the stochastic hybrid RD equation for a sin- for N(t) = n. Away from the boundaries, ∆ij = ∆ij gle chemical species, and derive moment equations for with the concentration. We introduce a two-state model for D stochastically-gated adsorption, and indicate how a non- ∆ij = [δi;j+1 + δi;j−1 − 2δi;j]; 0 < i < N ; (2.5) h2 linear adsorption rate can occur under dimerization, lead- ing to a moment closure problem. In Sect. III, we present n whereas ∆ij is modified at the boundaries i = 0; N in a detailed derivation of the path integral for the spatially order to be consistent with the boundary conditions (see discretized hybrid RD model. We perform the small noise [5] for details). The corresponding probability density expansion of the continuum path integral action in Sect. IV and derive leading order corrections to the determin- ProbfU(t) 2 (u; u + du);N(t) = ng = Pn(u; t)du (2.6) istic mean-field RD equation. In Sect. V we apply the theory to the example of stochastically-gated adsorption evolves according to a differential Chapman-Kolmogorov under dimerization and construct the alternative path in- (CK) equation for the stochastic hybrid system (2.4): tegral representation based on an RD master equation. Finally, the RG analysis is developed in Sect. VI. 20 1 3 @Pn X @ X n = − 4@ ∆ijuj + Fn(ui)A Pn(u; t)5 @t @ui i2J j2J II. STOCHASTIC HYBRID PDE X + QnmPm(u; t): (2.7) Consider a one-dimensional RD process in which the m2I reaction term depends on the current discrete state of It is important to note that the above spatial dis- the environment, N(t). The latter is assumed to evolve cretization is distinct from the scheme used to derive according to a continuous time Markov chain on I ⊂ an RD master equation [39{41]. In the derivation of Z+. Setting P (t) = [N(t) = n], we have the master n P Eq. (2.7) we assumed that h is sufficiently large so that equation the expected number of molecules per lattice site is large enough to write down a local concentration.