# Quantum Techniques for Reaction Networks.

1. IntroductionA "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:

[H.sub.2]O [??] [H.sup.+] + O[H.sup.-] [H.sub.2]O + [H.sup.+] [??] [H.sub.3][O.sup.+] (1)

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 [H.sub.2]O, the proton [H.sup.+], the hydroxyl ion O[H.sup.-], and the hydronium ion [H.sub.3][O.sup.+].

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 the rate 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 the master 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 identical classical 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 its 1th energy level, its state is a vector that we can formally write as [z.sup.l] for some formal variable z, or better, [z.sup.l] / [square root of (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

[PSI] = [[infinity].summation over (l=0)] [[psi].sub.l] [z.sup.l]/[square root of (l!)] (2)

where the coefficient [[psi].sub.l] [member of] 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.sup.[dagger]], which multiplies any power series by z:

[a.sup.[dagger]][PSI] = z[PSI]. (3)

We can decrease the energy level using the "lowering operator" a, which differentiates any power series:

a[PSI] = d/dz [PSI]. (4)

If we choose an inner product for which the vectors [z.sup.l] / [square root of (l!)] are orthonormal, these operators are indeed adjoint to each other:

<a[PSI], [PHI]> = <[PSI], [a.sup.[dagger]][PHI]>. (5)

In checking this one sees why the normalization factor is required.

We can think of [z.sup.l] 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 [z.sub.1] and [z.sub.2]:

[mathematical expression not reproducible] (6)

where [mathematical expression not reproducible] is the amplitude for having [l.sub.1] particles of the first kind and [l.sub.2] particles of the second kind. We have operators that create and annihilate particles of each kind:

[a.sup.[dagger].sub.i][PSI] = [z.sub.i][PSI],

[a.sub.i][PSI] = [partial derivative]/[partial derivative][z.sub.i] [PSI]. (7)

Each creation operator [a.sup.[dagger].sub.i] is adjoint to the corresponding annihilation operator [a.sub.i], and these operators obey the "canonical commutation relations":

[[a.sub.i], [a.sub.j]] = 0

[[a.sub.i.sup.[dagger]], [[a.sub.t.sup.[dagger]]] = 0

[[a.sub.i], [a.sup.[dagger].sub.j]] = [[delta].sub.ij]. (8)

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

[mathematical expression not reproducible] (9)

Here, [mathematical expression not reproducible] is the probability of having [l.sub.1] particles of the first kind and [l.sub.2] 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 [PSI] to describe a probability distribution, we require that all the coefficients be nonnegative and sum to 1. We call [PSI] with these properties a "mixed state". A special case is a monomial [mathematical expression not reproducible]. This is called a "pure state": it describes a completely definite situation where there are [l.sub.1] particles of the first kind and [l.sub.2] 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 Theorem 8 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 equation exactly 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]

[mathematical expression not reproducible] (10)

Here 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) [alpha]: the birth of one healthy cell, which has no input and one H as output

(ii) [beta]: the death of a healthy cell, which has one H as input and no output

(iii) [gamma]: the infection of a healthy cell, which has one H and one V as input, and one I as output

(iv) [delta]: "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) [epsilon]: the death of an infected cell, which has one I as input and no output

(vi) [zeta]: the death of a virion, which has one V as input and no output

So, we have a finite set of species:

S = {H, I, V} (11)

but reactions go between "complexes", which are finite linear combinations of species with natural number coefficients, such as H + V or 2 H + 3V or 0. (In the literature on reaction networks the complex 0, corresponding to nothing at all, is often denoted 0.) We can think of complexes as elements of [N.sup.s], 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. A graph consists of a set E of edges, a set V of vertices, and source and target maps s, t : E [right arrow] V saying where each edge starts and ends.

Then we may say the following.

Definition 2. A reaction network (S, K, T, s, t) consists of

(i) a finite set S of species

(ii) a finite set of complexes K [subset or equal to] [N.sup.S]

(iii) a graph s, t : T [right arrow] K with complexes as vertices and some finite set of reactions T as edges

In 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 write k = [absolute value of (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 [tau] as

s ([tau]) = ([s.sub.1] ([tau]), ..., [s.sub.k] ([tau])) [member of] [N.sup.k],

t ([tau]) = ([ta.sub.1] ([tau]), ..., [t.sub.k] ([tau])) [member of] [N.sup.k]. (12)

3. The Rate Equation

The amount of each species is represented by a "state" of the chemical reaction network. There are various options here. A pure state is a vector x [member of] [N.sup.k] with ith entry [x.sub.i] specifying the number of instances of the ith species. This generalizes in two ways. First, we define a classical state to be any vector x [member of] [[0, [infinity]).sup.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 a mixed state to be a probability distribution over the pure states. This assigns a probability [[psi].sub.l] to each l [member of] [N.sup.k]. 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 its rate constant.

Definition 3. A stochastic reaction network (S, K, T, s, t, r) consists of a reaction network together with a map r : T [right arrow] (0, [infinity]) assigning a rate 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 t the time derivative of a classical state x is the product of the following:

(i) the vector t([tau]) - s([tau]) whose ith component is the change in the number of the ith species due to the reaction [tau]

(ii) the concentration of each input species i of [tau] raised to the power given by the number of times it appears as an input, namely, [s.sub.i]([tau])

(iii) the rate constant r([tau]) of [tau]

Definition 4. The rate equation for a stochastic reaction network (S, K, T, s, t, r) is

[mathematical expression not reproducible] (13)

where x : R [right arrow] [[0, [infinity]).sup.k] and we have used multi-index notation to define

[mathematical expression not reproducible] (14)

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:

dH/dt = [alpha] - [beta]H - [gamma]HV

dI/dt = [gamma]HV - [epsilon]I

dV/dt = -[gamma]HV + [delta]I - [zeta]V. (15)

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:

[mathematical expression not reproducible] (16)

This reaction network gives the following rate equation:

d[Cl.sub.2]/dt = -[alpha][Cl.sub.2] + [beta][Cl.sup.2]

dCl/dt = 2[alpha][Cl.sub.2] - 2[beta][Cl.sup.2]. (17)

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 [beta] are proportional to [Cl.sup.2]. The reaction a produces two atoms of Cl while the reaction [beta] 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 [member of] [N.sup.k], a "mixed state" gives the probability [[psi].sub.l] that, for each i, exactly [l.sub.i] [member of] N instances of the ith species are present. The master equation says that

[mathematical expression not reproducible] (18)

for some matrix of numbers H determined by the stochastic reaction network. The matrix element [H.sub.l'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:

[mathematical expression not reproducible] (19)

To describe the matrix H([tau]), we need "falling powers". For any natural numbers n and p, we define the _pth falling power of n by

[n.sup.p]- = n{n - 1) ... (n - p + 1). (20)

This is the number of ways of choosing an ordered p-tuple of distinct elements from a set with 11 elements. Note that [n.sup.p]- = 0 if p > n. More generally, for any multi-indices l and m, we define

[mathematical expression not reproducible] (21)

This is the number of ways to choose, for each species i = l, ..., k, an ordered list of m; distinct things from a collection of [l.sub.i] things of that species. Thus, we expect a factor of this sort to appear in the master equation.

Using this notation, we define the matrix H(t) as follows:

[mathematical expression not reproducible] (22)

Here [delta] 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 [tau]: 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:

[mathematical expression not reproducible] (23)

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 variables [z.sub.1], ..., [z.sub.k], with the coefficient of

[mathematical expression not reproducible] (24)

being the probability [[psi].sub.l]. That is, we write any mixed state as

[mathematical expression not reproducible] (25)

Because a mixed state is a probability distribution, the coefficients [[psi].sub.l] must be nonnegative and sum to 1. Indeed, in this formalism mixed states are precisely the formal power series [PSI] with

[mathematical expression not reproducible] (26)

The simplest mixed states are the monomials [z.sup.l]; these are the pure 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 the stochastic Fock space to be the vector space R [[[z.sub.1], ..., [z.sub.k]]] of all real formal power series in the variables [z.sub.1], ..., [z.sub.k]. We treat the Hamiltonian H for the master equation as the operator on stochastic Fock space with

[mathematical expression not reproducible] (27)

for all

[mathematical expression not reproducible] (28)

where the matrix entries [H.sub.l'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. Let 1 [less than or equal to] i [less than or equal to] k. The creation operator [a.sup.[dagger].sub.i] is given by

[a.sup.[dagger].sub.i][PSI] = [z.sub.i][PSI]. (29)

This takes any pure state to the pure state with one additional instance of the ith species:

[mathematical expression not reproducible] (30)

The corresponding annihilation operator is given by formal differentiation:

[a.sub.i][PSI] = [partial derivative]/[partial derivative][z.sub.i] [PSI]. (31)

This takes any pure state to the pure state with one fewer instance of the ith species, but multiplied by a coefficient [l.sub.i]:

[mathematical expression not reproducible] (32)

This represents the fact that there are [l.sub.i] ways to annihilate one of the instances of the ith species present.

In what follows we use multi-index notation to define, for any n [member of] [N.sup.k],

[mathematical expression not reproducible] (33)

and

[mathematical expression not reproducible] (34)

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

[mathematical expression not reproducible] (35)

Proof. We start with a basic result on annihilation and creation operators.

Lemma 6. For any multi-indices l, m we have

[mathematical expression not reproducible] (36)

Proof. The first equation follows inductively from the definition of the creation operators, which implies

[mathematical expression not reproducible] (37)

The second follows inductively from the definition of annihilation operators, which implies

[mathematical expression not reproducible] (38)

Note that [l.sup.[m.bar]] = 0 if [m.sub.i] > [l.sub.i] for any i, and in this case we also have [a.sup.m][z.sup.l] = 0, so the equation [a.sup.m][z.sup.l] = [l.sup.[m.bar]][z.sup.l-m] holds in this case because both sides vanish.

We now prove the theorem. For any multi-index l, (23) gives

[mathematical expression not reproducible] (39)

Using Lemma 6 we obtain

[mathematical expression not reproducible] (40)

and thus

[mathematical expression not reproducible] (41)

The last step requires a little justification. Technically speaking, the monomials [z.sup.l] are not a basis of the stochastic Fock space R[[[z.sub.1], ..., [z.sub.k]]], because not every formal power series is a finite linear combination of monomials. However, these monomials form a "topological basis" of R[[[z.sub.1], ..., [z.sub.k]]], in the sense that every element of this vector space can be expressed as an convergent 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 [z.sup.l].

Next we investigate what the master equation says about the time evolution of expected values. For this, we start by defining

[mathematical expression not reproducible] (42)

for any formal power series

[mathematical expression not reproducible] (43)

In general the sum defining <[PSI]> may not converge, but it converges and equals 1 when [PSI] is a mixed state. Suppose O is an operator on the space of formal power series in the variables [z.sub.1], ..., [z.sub.k]. We define the expected value of O in the mixed state [PSI] to be <O[PSI]>, 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 <[PSI], O[PSI]>, where the brackets denote an inner product, in stochastic mechanics we use <O[PSI]>. For more details see [1, 2].

We are especially interested in the expected value of the number operators

[N.sub.i] = [a.sup.[dagger].sub.i][a.sub.i]. (44)

Note that

[mathematical expression not reproducible] (45)

Since [[psi].sub.l] is the probability of being in a pure state where there are [l.sub.i] 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

d/dt <[N.sub.i]{PSI](t)> (46)

assuming that [PSI](t) obeys the master equation. For this, we need to define the falling powers of an operator A as follows:

[A.sub.[p.bar]] = A(A - l) ... (A - p + 1) (47)

for any p [member of] N. If m is a multi-index we define

[mathematical expression not reproducible] (48)

This is well-defined because the number operators [N.sub.i] commute.

Lemma 7. For any multi-indices l and m we have

[mathematical expression not reproducible] (49)

Proof. By Lemma 6 and the definition of the number operators we have

[mathematical expression not reproducible] (50)

for all 1 [less than or equal to] i [less than or equal to] 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 write N for the vector, or list, of number operators ([N.sub.1], ..., [N.sub.k]). Then, given a mixed state [PSI], we set

<N[PSI]> = (<[N.sub.1][PSI]>, ..., <[N.sub.k][PSI]>). (51)

If the expected values here are well-defined, <N[PSI]> [member of] [[0, [infinity]i).sup.k] is a classical state. It is a simplified description of the mixed state [PSI], 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 [PSI](t) evolving in time according to the master equation, we have

[mathematical expression not reproducible] (52)

assuming the expected values and their derivatives exist.

Proof. First note using Lemmas 6 and 7 that for any multi-indices l, m, n and any natural number i we have

[mathematical expression not reproducible] (53)

Thus if at any given time

[mathematical expression not reproducible] (54)

and [PSI](t) evolves in time according to the master equation, we have

[mathematical expression not reproducible] (55)

Recall that N is our notation for the vector of number operators ([N.sub.1], ..., [N.sub.k]). Thus, the above equation is equivalent to

[mathematical expression not reproducible] (56)

The above theorem makes it clear that the rate equation is closely related to the master equation, since the former says

[mathematical expression not reproducible] (57)

while the latter implies

[mathematical expression not reproducible] (58)

Suppose we let x(t) be the classical state coming from the mixed state [PSI](t):

x (t) = <N[PSI](t)>. (59)

The rate equation for x(t) would then follow from the master equation for [PSI](t) if we had

[<N[PSI] (t)>.sup.m] = <[N.sup.[m.bar]][PSI](t)> (60)

for every multi-index m. This equation is not true in general. However, it should hold approximately in a suitable limit of large numbers. And surprisingly, it holds exactly for coherent states.

6. Coherent States

For any c [member of] [[0, [infinity]).sup.k], we define the coherent state with expected value c to be the mixed state

[[PSI].sub.c] = [e.sup.cxz]/[e.sup.c] (61)

where c x z is the dot product of the vectors c and z, and we set [mathematical expression not reproducible]. Equivalently,

[mathematical expression not reproducible] (62)

where [c.sup.n] and [z.sup.n] are defined as products in our usual way, and n! = [n.sub.1]! ... [n.sub.k]!. The name "coherent state" comes from quantum mechanics [10 ], where we think of the coherent state [[PSI].sub.c] as the quantum state that best approximates the classical state c. In the state [[PSI].sub.c], the probability of having [n.sub.i] things of the ith species is equal to

[mathematical expression not reproducible] (63)

This is precisely the definition of a Poisson distribution with mean equal to [c.sub.i]. The state [[PSI].sub.c] is thus a product of independent Poisson distributions.

Theorem 9. Given any stochastic reaction network, let [PSI](t) be a mixed state evolving in time according to the master equation. If [PSI](t) is a coherent state when t = [t.sub.0], then

x(t) = <N[PSI](t)> (64)

obeys the rate equation when t = [t.sub.0].

Proof. We prove this using a series of lemmas.

Lemma 10. For any multi-index m and any coherent state [[PSI].sub.c] we have

[a.sup.m][[PSI].sub.c] = [c.sup.m][[PSI].sub.c] (65)

and

[a.sup.[dagger]m][[PSI]sub.c] = [z.sup.m][[PSI].sub.c]. (66)

Proof. The second equation is immediate from the definition, while the first follows from

[mathematical expression not reproducible] (67)

Lemma 11. For any multi-index m we have

[N.sub.[m.bar]] = [a.sup.[dagger]m] [a.sup.m]. (68)

Proof. For any multi-index l, Lemmas 6 and 7 imply

[mathematical expression not reproducible] (69)

Since the states [z.sup.l] form a topological basis of R[[[z.sub.1], ..., [z.sub.k]]] and all the operators in question are continuous, the lemma follows.

Lemma 12. For any multi-index m and any [PSI] [member of] R[[[z.sub.1], ..., [z.sub.k]]] we have

<[a.sup.[dagger]m][PSI]> = <[PSI]>. (70)

Proof. We have

[mathematical expression not reproducible] (71)

Lemma 13. For any multi-index m and any coherent state [[PSI].sub.c] we have

<[N.sub.[m.bar]][[PSI].sub.c]> = [<N[[PSI].sub.c]>.sup.m]. (72)

Proof. By Lemmas 11 and 10 we have

[N.sub.[m.bar]][[PSI].sub.c] = [a.sup.[dagger]m][a.sup.m][[PSI].sub.c] = [c.sup.m][a.sup.[dagger]m][[PSI].sub.c] (73)

so by Lemma 12

<[N.sub.[m.bar]][[PSI].sub.c]> = <[c.sup.m] [a.sup.[dagger]m][[PSI].sub.c]> = [c.sup.m] <[[PSI].sub.c]> = [c.sup.m]. (74)

As a special case we have <[N.sub.i][[PSI].sub.c]> = [c.sub.i], so we also have

[<N[[PSI].sub.c]>.sup.m] = [c.sup.m] (75)

and it follows that <[N.sup.[m.bar]][[PSI].sub.c]> = [<N[PSI]>.sup.m].

To prove Theorem 9, we start by setting

x(t) = <N[PSI](t)> (76)

where [PSI](t) is a solution of the master equation. Theorem 8 implies that the time derivative of x equals

[mathematical expression not reproducible] (77)

If [PSI](t) is a coherent state at time t = [t.sub.0], Lemma 13 implies that

[mathematical expression not reproducible] (78)

so we have

[mathematical expression not reproducible] (79)

which is the rate equation at time [t.sub.0], 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, Theorem 8 implies that, as long as the expected values <[N.sup.[m.bar]][PSI](t)> are close to the powers [<N[PSI](t)>.sup.m], at least when m is the source of some reaction, the expected values x(t) = <N[PSI](t)> will approximately 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 [PSI](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].

https://doi.org/10.1155/2018/7676309

Data Availability

No data were used to support this study.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The author thank Jacob Biamonte, Brendan Fong, and the students at U. C. Riverside who helped come up with a preliminary version of the proof of Theorem 8, notably Daniel Estrada, Reeve Garrett, Michael Knap, Tu Pham, Blake Pollard, and Franciscus Rebro and thank John Rowlands and Blake Stacey for catching typos and other mistakes. This work was done at the Centre of Quantum Technology, in Singapore, and the Mathematics Department of U. C. Riverside.

References

[1] J. C. Baez and J. Biamonte, Quantum Techniques for Stochastic Mechanics, World Scientific, Singapore, 2018.

[2] J. C. Baez and B. Fong, "A Noether theorem for Markov processes," Journal of Mathematical Physics, vol. 54, no. 1, p. 013301, 2013.

[3] J. C. Baez and B. Fong, "Quantum techniques for studying equilibrium in reaction networks," Journal of Complex Networks, vol. 3, no. 1, pp. 22-34, 2015.

[4] M. Doi, "Second quantization representation for classical many-particle system," Journal of Physics A: Mathematical and General, vol. 9, no. 9, pp. 1465-1477, 1976.

[5] P. Grassberger and M. Scheunert, "Fock-space methods for identical classical objects," Fortschritte der Physik, vol. 28, no. 10, pp. 547-578, 1980.

[6] A. DeMoivre, in Miscellanea Analytica de Seriebus et Quadraturis, Book VII, pp. 194-197, Touson & Watts, London, UK, 1730.

[7] P.-S. Laplace, "Memoire sur les suites," Memoires de l'Academie Royale des Sciences de Paris, pp. 207-309, 1782.

[8] F. Horn and R. Jackson, "General mass action kinetics," Archive for Rational Mechanics and Analysis, vol. 47, pp. 81-116, 1972.

[9] M. Feinberg, Lectures On Reaction Networks, 1979, Available at https://crnt.osu.edu/LecturesOnReactionNetworks.

[10] J. R. Klauder and B. Skagerstam, Coherent States, World Scientific, Singapore, 1985.

[11] G. Bianconi and A.-L. Barabasi, "Bose-Einstein condensation in complex networks," Physical Review Letters, vol. 86, no. 24, pp. 5632-5635, 2001.

[12] M. Alberto Javarone and G. Armano, "Quantum-classical transitions in complex networks," Journal of Statistical Mechanics: Theory and Experiment, no. 4, p. P04019, 2013.

[13] J. Biamonte, M. Faccin, and M. De Domenico, Complex networks: from classical to quantum, Available as https://arxiv.org/abs/1702.08459.

[14] R. Aris, "Prolegomena to the rational analysis of systems of chemical reactions," Archive for Rational Mechanics and Analysis, vol. 19, pp. 81-99, 1965.

[15] A. Korobeinikov, "Global properties of basic virus dynamics models," Bulletin of Mathematical Biology, vol. 66, no. 4, pp. 879-883, 2004.

[16] I. Koch, "Petri nets - A mathematical formalism to analyze chemical reaction networks," Molecular Informatics, vol. 29, no. 12, pp. 838-843, 2010.

[17] I. Koch, W. Reisig, and F. Schreiber, Modeling in Systems Biology, Springer, 2011.

[18] D. J. Wilkinson, Stochastic Modelling for Systems Biology, Taylor and Francis, New York, NY, USA, 2006.

[19] D. F. Anderson and T. G. Kurtz, "Continuous time Markov chain models for chemical reaction networks," in Design and Analysis of Biomolecular Circuits: Engineering Approaches to Systems and Synthetic Biology, H. Koeppl et al., Ed., pp. 3-42, Springer, 2011.

[20] D. F. Anderson, G. Craciun, and T. G. Kurtz, "Product-form stationary distributions for deficiency zero chemical reaction networks," Bulletin of Mathematical Biology, vol. 72, no. 8, pp. 1947-1970, 2010.

John C. Baez (iD) (1,2)

(1) Department of Mathematics, University of California, Riverside, CA 92521, USA

(2) Centre for Quantum Technologies, National University of Singapore, Singapore 117543

Correspondence should be addressed to John C. Baez; baez@math.ucr.edu

Received 28 June 2018; Accepted 13 September 2018; Published 1 November 2018

Academic Editor: Giampaolo Cristadoro

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Baez, John C. |

Publication: | Advances in Mathematical Physics |

Article Type: | Report |

Geographic Code: | 1USA |

Date: | Jan 1, 2018 |

Words: | 5993 |

Previous Article: | Decay Rate and Energy Gap for the Singularity Solution of the Inhomogeneous Landau-Lifshitz Equation on [S.sup.2]. |

Next Article: | Stability and Energy-Casimir Mapping for Integrable Deformations of the Kermack-McKendrick System. |

Topics: |