# Quantum Techniques for Reaction Networks

Authors: John C. Baez
Journal: Advances in Mathematical Physics
 (2018)

---

## Abstract

Reaction networks are a general formalism for describing collections of classical entities interacting in a random way. While reaction networks are mainly studied by chemists, they are equivalent to Petri nets, which are used for similar purposes in computer science and biology. As noted by Doi and others, techniques from quantum physics, such as second quantization, can be adapted to apply to such systems. Here we use these techniques to study how the “master equation” describing stochastic time evolution for a reaction network is related to the “rate equation” describing the deterministic evolution of the expected number of particles of each species in the large-number limit. We show that the relation is especially strong when a solution of master equation is a “coherent state”, meaning that the numbers of entities of each kind are described by independent Poisson distributions. Remarkably, in this case the rate equation and master equation give the exact same formula for the time derivative of the expected number of particles of each species.


---

## Body

## 1. Introduction

A “reaction network” describes how various kinds of classical particles can interact and turn into other kinds. They are commonly used in chemistry. For example, this reaction network gives a very simplified picture of what is happening in a glass of water:(1)H2O⇄H++OH-H2O+H+⇄H3O+The reactions here are all reversible, but this is not required by the reaction network formalism. The relevant particles here are not atoms, but rather “species” including the water molecule H2O, the proton H+, the hydroxyl ion OH-, and the hydronium ion H3O+.While the dynamics of the reactions are fundamentally quantum-mechanical, reaction networks give a simplified description of what is going on. First, given a “rate constant” for each reaction, one can write down a differential equation called therate equation, which describes the time evolution of the concentration of each species in a deterministic way. This is especially useful in the limit where there are many particles of each species. At a more fine-grained level, one can describe reactions stochastically, using themaster equation. In the rate equation the concentration of each species is treated as a continuous variable, but, in the master equation, the number of particles of each species is an integer.In this paper, we explain how techniques from quantum physics and in particular second quantization can be used to study reaction networks. This work is part of a broader program emphasizing the links between quantum mechanics and a subject we call “stochastic mechanics”, where probabilities replace amplitudes [1–3].Using second quantization methods to model classical stochastic systems goes back at least to 1976, with the work of Doi [4]. It has been further developed by Grassberger, Scheunert [5], and many others. The big surprise one meets at the very beginning of this subject is that the canonical commutation relations between creation and annihilation operators, long viewed as a hallmark of quantum theory, are also perfectly suited to systems of identicalclassical particles interacting in a probabilistic way.In quantum theory, the simplest example of a Fock space is used to describe the states of a quantum harmonic oscillator. When such an oscillator is in itslth energy level, its state is a vector that we can formally write as zl for some formal variable z, or better, zl/l!, where the normalization factor comes in handy later. A general quantum state is a linear combination of such vectors, so we can write it as a power series(2)Ψ=∑l=0∞ψlzll!where the coefficient ψl∈C is the amplitude for the oscillator to be in its lth energy level. We can increase the energy level using the “raising operator” a†, which multiplies any power series by z:(3)a†Ψ=zΨ.We can decrease the energy level using the “lowering operator” a, which differentiates any power series:(4)aΨ=ddzΨ.If we choose an inner product for which the vectors zl/l! are orthonormal, these operators are indeed adjoint to each other:(5)aΨ,Φ=Ψ,a†Φ.In checking this one sees why the normalization factor is required.We can think ofzl as a state containing l quanta of energy. In applications to quantum physics, we think of these quanta as actual particles, for example, photons. If we have k different kinds of particles or a single kind of particle with k different states, we describe the quantum state of a collection of particles using a power series in k variables. For example, if k=2 there is a Fock space consisting of power series in two variables z1 and z2:(6)Ψ=∑l1=0∞∑l2=0∞ψl1,l2z1l1l1!z2l2l2!where ψl1,l2∈C is the amplitude for having l1 particles of the first kind and l2 particles of the second kind. We have operators that create and annihilate particles of each kind:(7)ai†Ψ=ziΨ,aiΨ=∂∂ziΨ.Each creation operator ai† is adjoint to the corresponding annihilation operator ai, and these operators obey the “canonical commutation relations”: (8)ai,aj=0ai†,aj†=0ai,aj†=δij.These rules are the basis of many practical calculations in particle physics, where products of creation and annihilations operators are used to model processes in which particles interact and turn into other particles.Doi noticed that all this has a stochastic analogue where we work with probabilities rather than amplitudes. For example, there is a stochastic Fock space whose elements are power series(9)Ψ=∑l1=0∞∑l2=0∞ψl1,l2zl1zl2.Here, ψl1,l2∈R is theprobability of having l1 particles of the first kind and l2 particles of the second kind. The formalism is somewhat simpler than in the quantum case; for example, the normalization factors are no longer required. For Ψ to describe a probability distribution, we require that all the coefficients be nonnegative and sum to 1. We call Ψ with these properties a “mixed state”. A special case is a monomial z1l1z2l2. This is called a “pure state”: it describes a completely definite situation where there are l1 particles of the first kind and l2 of the second kind.This method of describing a probability distribution using a power series is far from new. It was introduced by DeMoivre [6] and perfected by Laplace [7] in a paper published in 1782. He dubbed these power series “generating functions”, and they have been used in probability theory ever since. Doi’s new realization was that spaces of generating functions are analogous to the Fock spaces in quantum theory. In particular, products of creation and annihilation operators can be used to describe stochastic processes in which objects interact and turn into other objects.A natural example is chemistry, where reactions are often described stochastically using the master equation or deterministically using the rate equation. There is a long line of work on both of these equations, initiated by Horn and Jackson [8] along with Feinberg [9] in the late 1970s and continuing to this day. This line of work does not use quantum techniques, but it might profit from them and also contribute to their development. To illustrate how, here we use these techniques to study how the master equation reduces to the rate equation in the large-number limit. This is formally similar to the “classical limit” in which a quantum field theory reduces to a classical field theory.In Theorem8 we start with the master equation and derive an equation for the rate of change of the expected number of particles of each species. This equation is similar to the rate equation and reduces to it in a certain approximation. However, in Theorem 9 we show that this equation gives the rate equationexactly when the numbers of particles of each species are described by independent Poisson distributions. This relies on a nice fact: the generating function of a product of independent Poisson distributions is an eigenvector of all the annihilation operators.In quantum mechanics, eigenvectors of the annihilation operators are called “coherent states”. Many facts about classical mechanics have simple generalizations to quantum mechanics if we restrict attention to coherent states [10]. Here we are seeing a similar phenomenon in stochastic mechanics: the deterministic dynamics given by the rate equation match the stochastic dynamics given by the master equation in the special case of coherent states. We explore this further in a companion paper [3].The application of quantum techniques to reaction networks is just one of many interesting connections between quantum physics and network theory. For example, the evolution of networks displays phase transitions analogous to Bose–Einstein condensation [11, 12]. For more, see the review by Biamonte, Faccin, and De Domenico [13].


## 2. Reaction Networks

A reaction network, first defined by Aris [14], consists of a set of different kinds of particles—each kind being called a “species”—together with a set of processes, called “reactions”, that turn collections of particles of various species into collections of particles of other species. As the name suggests, it is convenient to draw a reaction network as a graph.While reaction networks first arose in chemistry, their uses are not limited to this context. Here is an example that arose in work on HIV: the human immunodeficiency virus is [15](10)0⇄βαHH+V→γI→δI+VI→ϵ0V→ζ0Here we have three species: (i)
H: healthy white blood cells(ii)
I: infected white blood cells(iii)
V: virions (that is, individual virus particles)We also have six reactions:(i)
α: the birth of one healthy cell, which has no input and one H as output(ii)
β: the death of a healthy cell, which has one H as input and no output(iii)
γ: the infection of a healthy cell, which has one H and one V as input, and one I as output(iv)
δ: “production”, or the reproduction of the virus in an infected cell, which has one I as input and one I and one V as output(v)
ϵ: the death of an infected cell, which has one I as input and no output(vi)
ζ: the death of a virion, which has one V as input and no outputSo, we have a finite set of species:(11)S=H,I,Vbut reactions go between “complexes”, which are finite linear combinations of species with natural number coefficients, such as H+V or 2H+3V or 0. (In the literature on reaction networks the complex 0, corresponding to nothing at all, is often denoted ∅.) We can think of complexes as elements of NS, the set of functions from S to the natural numbers. The reaction network is a graph with certain complexes as vertices and reactions as edges.We can formalize this as follows. There are various kinds of graphs; the kind we want is sometimes called “directed multigraphs” or “quivers”, but, to keep our terminology simple, we make the following definition.Definition 1.
Agraph consists of a set E ofedges, a set V ofvertices, andsource andtarget maps s,t:E→V saying where each edge starts and ends.Then we may say the following.Definition 2.
Areaction network  (S,K,T,s,t) consists of (i)
a finite setS ofspecies(ii)
a finite set ofcomplexes  K⊆NS(iii)
a graphs,t:T→K with complexes as vertices and some finite set ofreactions  T as edgesIn the chemistry literature this sort of thing is called a “chemical reaction network”, but we want to emphasize that they are useful more generally, so we drop the adjective. Petri nets are an alternative formalism that is entirely equivalent to reaction networks, often used in subjects other than chemistry [16–18]. In the Petri net literature the species are called “places” and the reactions are called “transitions”. For a more complete translation guide between reaction networks and Petri nets, see [3].For convenience we shall often writek=|S| for the number of species present in a reaction network and identify the set S with the set {1,…,k}. This lets us write any complex as a k-tuple of natural numbers. In particular, we write the source and target of any reaction τ as(12)sτ=s1τ,…,skτ∈Nk,tτ=t1τ,…,tkτ∈Nk.


## 3. The Rate Equation

The amount of each species is represented by a “state” of the chemical reaction network. There are various options here. Apure state is a vector x∈Nk with ith entry xi specifying the number of instances of the ith species. This generalizes in two ways. First, we define aclassical state to be any vector x∈[0,∞)k of nonnegative real numbers and think of such a state as specifying the expected number or concentration of the species present. Second, we define amixed state to be a probability distribution over the pure states. This assigns a probability ψl to each l∈Nk. The rate equation describes the time evolution of a classical state, while the master equation describes the time evolution of a mixed state.To write down either of these equations, we must first specify the system’s dynamics. In this paper we consider systems that evolve by the law of mass action [14]. Very roughly, this law states that the rate at which reactions occur is proportional to the product of the concentrations of all its input species. We call the constant of proportionality for each reaction itsrate constant.Definition 3.
Astochastic reaction network  (S,K,T,s,t,r) consists of a reaction network together with a map r:T→(0,∞) assigning arate constant to each reaction.This concept is equivalent to another concept in the literature, namely, that of a stochastic Petri net [3].Given a stochastic reaction network, the rate equation says that for each reactionτ the time derivative of a classical state x is the product of the following:(i)
the vectort(τ)-s(τ) whose ith component is the change in the number of the ith species due to the reaction τ(ii)
the concentration of each input speciesi of τ raised to the power given by the number of times it appears as an input, namely, si(τ)(iii)
the rate constantr(τ) of τDefinition 4.
Therate equation for a stochastic reaction network (S,K,T,s,t,r) is(13)ddtxt=∑τ∈Trτtτ-sτxtsτwhere x:R→[0,∞)k and we have used multi-index notation to define(14)xsτ=x1s1τ⋯xkskτ.For the HIV reaction network in (10), if we also use the Greek letter names of the reactions as names for their rate constants, we get this rate equation:(15)dHdt=α-βH-γHVdIdt=γHV-ϵIdVdt=-γHV+δI-ζV.This example involves only complexes where the coefficient of each species is at most 1. It is worthwhile comparing an example with larger coefficients, for example, the dissociation of a diatomic gas into atoms and the recombination of these atoms into a molecule:(16)Cl2⇄βα2Cl.This reaction network gives the following rate equation:(17)dCl2dt=-αCl2+βCl2dCldt=2αCl2-2βCl2.The recombination of two atoms of chlorine into a diatomic molecule occurs at a rate proportional to the square of the concentration of atomic chlorine, so both terms arising from reaction β are proportional to Cl2. The reaction α produces two atoms of Cl while the reaction β uses up two atoms of Cl, so both terms in the second equation have a coefficient of 2. All this follows from (13).


## 4. The Master Equation

The master equation describes the time evolution of mixed states. For the rest of this paper we fix a stochastic reaction network(S,K,T,s,t,r). Recall that, for each l∈Nk, a “mixed state” gives the probability ψl that, for each i, exactly li∈N instances of the ith species are present. Themaster equation says that(18)ddtψl′=∑l∈NkHl′lψlfor some matrix of numbers H determined by the stochastic reaction network. The matrix element Hl′l is the probability per time for the pure state l to evolve to the pure state l′. This probability is a sum over reactions:(19)Hl′l=∑τ∈THτl′l.To describe the matrixH(τ), we need “falling powers”. For any natural numbers n and p, we define the pthfalling power of n by (20)np_=nn-1⋯n-p+1.This is the number of ways of choosing an ordered p-tuple of distinct elements from a set with n elements. Note that np_=0 if p>n. More generally, for any multi-indices l and m, we define(21)lm_=l1m1_⋯lkmk_.This is the number of ways to choose, for each species i=1,…,k, an ordered list of mi distinct things from a collection of li things of that species. Thus, we expect a factor of this sort to appear in the master equation.Using this notation, we define the matrixH(τ) as follows:(22)Hτl′l=rτlsτ_δl′,l+tτ-sτ-δl′l.Here δ is the Kronecker delta, which equals 1 if its subscripts are equal and 0 otherwise. The first term describes the rate for the complex l to become the complex l′ via the reaction τ: it equals the rate constant for this reaction multiplied by the number of ways this reaction can occur. The second term describes the rate at which l goes away due to this reaction. Thus, the size of the second term is equal to that of the first, but its sign is opposite.Putting together (19) and (22), we obtain this formula for the Hamiltonian:(23)Hl′l=∑τ∈Trτlsτ_δl′,l+tτ-sτ-δl′l.


## 5. The Stochastic Fock Space

We can understand the Hamiltonian for the master equation at a deeper level using techniques borrowed from quantum physics. The first step is to introduce a stochastic version of Fock space. To do this, we write any mixed state as a formal power series in some variablesz1,…,zk, with the coefficient of(24)zl=z1l1⋯zklkbeing the probability ψl. That is, we write any mixed state as(25)Ψ=∑l∈Nkψlzl.Because a mixed state is a probability distribution, the coefficients ψl must be nonnegative and sum to 1. Indeed, in this formalismmixed states are precisely the formal power series Ψ with(26)∑l∈Nkψl=1,ψl≥0.The simplest mixed states are the monomials zl; these are thepure states, where there is a definite number of things of each species.In quantum mechanics we use a similar sort of Fock space, but with power series having complex rather than real coefficients. To obtain a Hilbert space, we often restrict attention to power series obeying for which a certain norm is finite. However, power series not obeying this condition are also useful, for example, to ensure that operators of interest are everywhere defined, instead of merely densely defined.Thus, for the purposes of studying the master equation, we define thestochastic Fock space to be the vector space R[[z1,…,zk]] of all real formal power series in the variables z1,…,zk. We treat the Hamiltonian H for the master equation as the operator on stochastic Fock space with(27)HΨ=∑l,l′∈NkHl′lψlzl′for all(28)Ψ=∑l∈Nkψlzl,where the matrix entries Hl′l are given by (23).This allows us to express the Hamiltonian in terms of certain special operators on the stochastic Fock space: the creation and annihilation operators. Our notation here follows that used in quantum physics, where the creation and annihilation operators are adjoints of each other. Let1≤i≤k. Thecreation operator  ai† is given by(29)ai†Ψ=ziΨ.This takes any pure state to the pure state with one additional instance of the ith species:(30)ai†z1l1⋯zili⋯zklk=z1l1⋯zili+1⋯zklk.The correspondingannihilation operator is given by formal differentiation:(31)aiΨ=∂∂ziΨ.This takes any pure state to the pure state with one fewer instance of the ith species, but multiplied by a coefficient li:(32)aiz1l1⋯zili⋯zklk=liz1l1⋯zili-1⋯zklk.This represents the fact that there are li ways to annihilate one of the instances of the ith species present.In what follows we use multi-index notation to define, for anyn∈Nk,(33)an=a1n1⋯aknkand(34)a†n=a1†n1⋯ak†nk.These expressions are well-defined because all the annihilation operators commute with each other and similarly for the creation operators.Theorem 5.
For any stochastic reaction network, the Hamiltonian is given by(35)H=∑τ∈Trτa†tτ-a†sτasτ.Proof. We start with a basic result on annihilation and creation operators.Lemma 6.
For any multi-indicesl,m we have(36)a†mzl=zl+m,amzl=lm_zl-m.Proof.
The first equation follows inductively from the definition of the creation operators, which implies(37)ai†zl=z1l1⋯zili+1⋯zklk.The second follows inductively from the definition of annihilation operators, which implies(38)aizl=liz1l1⋯zili-1⋯zklk.Note that lm_=0 if mi>li for any i, and in this case we also have amzl=0, so the equation amzl=lm_zl-m holds in this case because both sides vanish.We now prove the theorem. For any multi-indexl, (23) gives(39)Hzl=∑l′∈NkHl′lzl′=∑τ∈Trτlsτ_zl+tτ-sτ-zl.Using Lemma 6 we obtain(40)Hzl=∑τ∈Trτa†tτ-a†sτasτzland thus(41)H=∑τ∈Trτa†tτ-a†sτasτ.The last step requires a little justification. Technically speaking, the monomials zl are not a basis of the stochastic Fock space R[[z1,…,zk]], because not every formal power series is afinite linear combination of monomials. However, these monomials form a “topological basis” of R[[z1,…,zk]], in the sense that every element of this vector space can be expressed as anconvergent infinite linear combination of these monomials, with respect to a certain topology. (In this topology, a sequence of formal power series converges if each coefficient converges.) The annihilation and creation operators and indeed all the operators discussed in this paper are continuous in this topology. So, to check equations between these operators, it suffices to check them on the monomials zl.Next we investigate what the master equation says about the time evolution of expected values. For this, we start by defining(42)Ψ=∑l∈Nkψlfor any formal power series(43)Ψ=∑l∈Nkψlzl.In general the sum defining 〈Ψ〉 may not converge, but it converges and equals 1 when Ψ is a mixed state. Suppose O is an operator on the space of formal power series in the variables z1,…,zk. We define theexpected value of O in the mixed state Ψ to be 〈OΨ〉, assuming this converges. This generalizes the usual definition of the expected value of a random variable, in a way that emphasizes the analogy between stochastic mechanics and quantum mechanics. While in quantum mechanics we compute expected values of observables using expressions of the form 〈Ψ,OΨ〉, where the brackets denote an inner product, in stochastic mechanics we use 〈OΨ〉. For more details see [1, 2].We are especially interested in the expected value of thenumber operators(44)Ni=ai†ai.Note that(45)NiΨ=∑l∈Nkliψl.Since ψl is the probability of being in a pure state where there are li things of the ith species, the above formula indeed gives the expected value of the number of things of this species.Our goal is to compute(46)ddtNiΨtassuming that Ψ(t) obeys the master equation. For this, we need to define the falling powers of an operator A as follows:(47)Ap_=AA-1⋯A-p+1for any p∈N. If m is a multi-index we define(48)Nm_=N1m1_⋯Nkmk_.This is well-defined because the number operators Ni commute.Lemma 7.
For any multi-indicesl and m we have(49)Nm_zl=lm_zl.Proof.
By Lemma6 and the definition of the number operators we have(50)Nizl=lizlfor all 1≤i≤k. The lemma follows directly from this.There is a way to extract a classical state from a mixed state, which lets us connect the rate equation to the master equation. To do this, we writeN for the vector, or list, of number operators (N1,…,Nk). Then, given a mixed state Ψ, we set(51)NΨ=N1Ψ,…,NkΨ.If the expected values here are well-defined, 〈NΨ〉∈[0,∞)k is a classical state. It is a simplified description of the mixed state Ψ, which only records the expected number of instances of each species.The next result says how this sort of classical state evolves in time, assuming that the mixed state it comes from obeys the master equation.Theorem 8.
For any stochastic reaction network and any mixed stateΨ(t) evolving in time according to the master equation, we have(52)ddtNΨt=∑τ∈Trτsτ-tτNsτ_Ψtassuming the expected values and their derivatives exist.Proof.
First note using Lemmas6 and 7 that for any multi-indices l,m,n and any natural number i we have(53)Nia†manzl=ln_Nia†mzl-n=ln_Nizl-n+m=ln_li-ni+mi.Thus if at any given time(54)Ψt=∑l∈Nkψltzland Ψ(t) evolves in time according to the master equation, we have(55)ddtNiΨt=NiHΨt=∑τ∈T∑l∈NkrτNia†tτ-a†sτasτψltzl=∑τ∈T∑l∈Nkrτlsτ_siτ-tiτψlt=∑τ∈T∑l∈Nkrτsiτ-tiτNsτ_zlψlt=∑τ∈Trτsiτ-tiτNsτ_Ψt.Recall that N is our notation for the vector of number operators (N1,…,Nk). Thus, the above equation is equivalent to(56)ddtNΨt=∑τ∈Trτsτ-tτNsτ_Ψt.The above theorem makes it clear that the rate equation is closely related to the master equation, since the former says(57)ddtxt=∑τ∈Trτtτ-sτxtsτwhile the latter implies(58)ddtNΨt=∑τ∈Trτsτ-tτNsτ_Ψt.Suppose we let x(t) be the classical state coming from the mixed state Ψ(t):(59)xt=NΨt.The rate equation for x(t) would then follow from the master equation for Ψ(t) if we had(60)NΨtm=Nm_Ψtfor every multi-index m. This equation is not true in general. However, it should holdapproximately in a suitable limit of large numbers. And surprisingly, it holdsexactly for coherent states.


## 6. Coherent States

For anyc∈[0,∞)k, we define thecoherent state with expected value c to be the mixed state(61)Ψc=ec·zecwhere c·z is the dot product of the vectors c and z, and we set ec=ec1+…+cn. Equivalently,(62)Ψc=1ec∑n∈Nkcnn!zn,where cn and zn are defined as products in our usual way, and n!=n1!⋯nk!. The name “coherent state” comes from quantum mechanics [10], where we think of the coherent state Ψc as the quantum state that best approximates the classical state c. In the state Ψc, the probability of having ni things of the ith species is equal to(63)e-cicinini!.This is precisely the definition of a Poisson distribution with mean equal to ci. The state Ψc is thus a product of independent Poisson distributions.Theorem 9.
Given any stochastic reaction network, letΨ(t) be a mixed state evolving in time according to the master equation. If Ψ(t) is a coherent state when t=t0, then(64)xt=NΨtobeys the rate equation when t=t0.Proof. We prove this using a series of lemmas.Lemma 10.
For any multi-indexm and any coherent state Ψc we have(65)amΨc=cmΨcand(66)a†mΨc=zmΨc.Proof.
The second equation is immediate from the definition, while the first follows from(67)aiΨc=∂∂ziec·zec=ciΨc.Lemma 11.
For any multi-indexm we have(68)Nm_=a†mam.Proof.
For any multi-indexl, Lemmas 6 and 7 imply(69)a†mamzl=lm_a†mzl-m=lm_zl=Nm_zl.Since the states zl form a topological basis of R[[z1,…,zk]] and all the operators in question are continuous, the lemma follows.Lemma 12.
For any multi-indexm and any Ψ∈R[[z1,…,zk]] we have(70)a†mΨ=Ψ.Proof.
We have(71)a†mΨ=∑l∈Nkψlzl+m=∑l∈Nkψl=Ψ.Lemma 13.
For any multi-indexm and any coherent state Ψc we have(72)Nm_Ψc=NΨcm.Proof.
By Lemmas11 and 10 we have(73)Nm_Ψc=a†mamΨc=cma†mΨcso by Lemma 12(74)Nm_Ψc=cma†mΨc=cmΨc=cm.As a special case we have 〈NiΨc〉=ci, so we also have(75)NΨcm=cmand it follows that 〈Nm_Ψc〉=〈NΨ〉m.
⁡
To prove Theorem9, we start by setting(76)xt=NΨtwhere Ψ(t) is a solution of the master equation. Theorem 8 implies that the time derivative of x equals(77)x˙t=∑τ∈Trτsτ-tτNsτ_Ψt.If Ψ(t) is a coherent state at time t=t0, Lemma 13 implies that(78)Nsτ_Ψt0=NΨt0sτ=xt0sτso we have(79)x˙t0=∑τ∈Trτsτ-tτxt0sτwhich is the rate equation at time t0, as desired.In general, the above result applies at only one moment in time. The reason is that if a solution of the master equation is a coherent state at some time, it need not be a coherent state at later (or earlier) times. Still, Theorem8 implies that, as long as the expected values 〈Nm_Ψ(t)〉 areclose to the powers 〈NΨ(t)〉m, at least when m is the source of some reaction, the expected values x(t)=〈NΨ(t)〉 willapproximately obey the rate equation. This could be studied in detail using either the techniques introduced here or those discussed by Anderson and Kurtz [19].Furthermore, in some cases whenΨ(t) is initially a coherent state it continues to be so for all times. In this case, x(t) continues to exactly obey the rate equation. For example, this is true when all the complexes in the reaction network consist of a single species. It is also true for a large class of equilibrium solutions of the master equation, as shown by Anderson, Craciun, and Kurtz [20]. For more on how quantum techniques apply to this situation, see the companion to this paper [3].



---

*Source: 7676309-2018-11-01.xml*