I'm trying to understand how the Doi-Peliti (DP) action is constructed, and specifically how they compute expectation values. To this end, I've been using the book by Taüber as a reference (Critical Dynamics: A Field Theory Approach to Equilibrium and Non-Equilibrium Scaling Behaviour).
The point I seem to be missing is during the discretized → path integral procedure of the expectation value of an observable. For a single chemical species with lattice occupation numbers {ni} this is defined as ⟨A(t)⟩=∑{ni}A({ni})P({ni},t)
where
P({ni},t) denotes the probability of observing a configuration
{ni} and follows a master-type equation. In the DP formalism, the chemical reactant is assigned a site-specific bossonic ladder algebra
ai,
a†i and one finds that the expectation value above may be expressed in this language by
⟨A(t)⟩=⟨P|A({a†iai})|Φ(t)⟩
Here the projection operator
⟨P|=⟨0|∏ieai, and the state vector
|Φ(t)⟩=∑{ni}P({ni},t)|{ni}⟩ satisfies the imaginary time Schrödinger equation
∂t|Φ(t)⟩=−H({a†i},{ai})|Φ(0)⟩
(
H({a†i},{ai}) is meant to indicate that
H is normal-ordered).
By shifting the operator ∏ieai in the above expression for ⟨A(t)⟩ over to the right, one obtains ⟨A(t)⟩=⟨0|˜A({a†i→1},{ai})e−H({a†i→1+a†i},{ai})t|˜Φ(0)⟩
in which
˜A({1},{ai}) is obtained from
A by normal ordering and replacing
ai by
1 (e.g.
a†iaia†jaj→aiδij+aiaj), and
|˜Φ(0)⟩=∏ieai|Φ(0)⟩
Here comes the part I seem to fail to understand. If we denote by U(t2,t1)=e−H({a†i→1+a†i},{ai})(t2−t1)
then clearly
U(t2,t1)=U(t2,t′)U(t′,t1). We may thus split the time-evolution operator
U in
⟨A(t)⟩ up into many pieces and insert the completeness relation
1=∫∏idϕ∗idϕi2πie∑iϕ∗iϕi|ϕ⟩⟨ϕ|
(
i in the denominator is the imaginary unit and
|ϕ⟩ is a coherent state) inbetween each time-step to obtain
⟨A(t)⟩=∫(∏i,kdϕ∗i(tk)dϕi(tk)2πi)⟨0|˜A({1},{ai})|ϕ(tf)⟩(∏j⟨ϕ(tj)|U(tj,tj−1)|ϕ(tj−1)⟩)×⟨ϕ(t0)|˜Φ(0)⟩
The matrix elements
⟨ϕ(tj)|U(tj,tj−1)|ϕ(tj−1)⟩
are easily calculated. However, to me it seems that
˜A({1},{ai})|ϕ(tf)⟩=˜A({1},{ϕi(tf)})|ϕ(tf)⟩
since
ai|ϕ⟩=ϕi|ϕ⟩. In particular, I don't find it obvious how the above tends to the path integral
⟨A(t)⟩=∫∏iD[ϕ∗i,ϕi]˜A({1},{ϕi(t)})e−A[ϕ∗i,ϕi]
(for some action
A I leave unspecified), as it seems as though the observable
˜A should only be evaluated at the final point
ϕ(tf).
Sorry for the very long message. Any help would be greatly appreciated!
This post imported from StackExchange Physics at 2017-08-11 12:47 (UTC), posted by SE-user john