Chain statistics as a path integral¶
Consider a single polymer chain, defined by a space-curve representation. The probability distribution (or Green function) G(b|a) of the chain extending from a to b is
where the state a denotes location xa at arclength sa and state b denotes location xb at arclength sb. The summation over emph{all paths} indicates a sum of all possible chain conformations (or paths) that begin at state a and end at state b. This concept of a path integral is extensively developed by Feynman and Hibbs~cite{feynman} to describe the statistical treatment of quantum mechanical particles. For polymer statistics, the weighting of a path (or chain conformation) is governed by statistical mechanics. The statistical weight ϕ[x(s)] of a path x(s) is given by the integration of the Boltzmann weighting of the path. Define the Action S for the particular path to be
Thus, the statistical weight of a path is given by
where A is a normalization constant to be discussed later.
The Hamiltonian \mathcal{H} includes differential elastic deformation energy e_{elas}(x',x'',...;s) and a potential energy per unit length V(x,s). The Green function G(b|a) is written in %the form path-integral form
Currently the derivation of the path model is in one dimension; however, the ideas are easily extended to higher dimensions.
Schematic of path summation in the space-arclength plane¶
Gaussian Chain model Green function—explicit path integration¶
Consider a polymer chain defined by the Gaussian Chain model in 1-D with no external field (V=0). The action is given by
where n=s/b is a dimensionless arclength parameter. This gives a statistical weight
For this treatment, convert the path integral into a discrete-chain representation. This procedure provides a practical method to perform the summation over all possible chain conformations. Break the chain into M segments of length \delta = (n_{b}-n_{a})/M. The position of the m, and the endpoints are x_{0}=x_{a} and x_{M}=x_{b}. The limit M \rightarrow \infty recovers our continuous chain representation
Discrete-chain representation converts the chain into M segments (in this case M=6) of length \delta=(n_{b}-n_{a})/M.¶
The statistical weight is now approximated as
and the Green function is approximated as
The x_{1} integral is given by
where x_{1}^{\star}=\frac{x_{0}+x_{2}}{2}.
The subsequent x_{2} integral is given by
where x_{2}^{\star}=\frac{x_{0}+2x_{3}}{3}. The x_{m} integral after performing integrals x_{1}, x_{2},ldots, x_{m-1} is
This gives the Green function
where A is set to ensure \int dx_{b} G(b|a)=1
Gaussian Chain model Green function—Schrödinger equation¶
Generally, it is more convenient to develop the path integral into a diffusion equation (or Schr"{o}dinger equation) to solve for the chain statistics. Consider the Gaussian Chain model in 1-D with an external potential, giving an action
where we set n_{a}=0 and n_{b}=N (without loss of generality), and assume V=V(x) only. The chain obeys Markovian statistics; thus, the Green function for a chain that starts at A (emph{i.e.} position x_{A} and arclength n_{A}) to state C (emph{i.e.} x_{C}, n_{C}) through intermediate state B (emph{i.e.} x_{B}, n_{B}) is
Schematic of a chain that obeys Markovian statistics¶
Take the state C to differ an infinitesimal length from state B by setting n_{C}=N+\delta and n_{B}=N (n_{A}=0). For simplicity, set x_{A}=x_{0}, x_{B}=y, x_{C}=x, giving
In the limit \delta \rightarrow 0, the Green function G(x, N+\delta | y,N) will be dominated by a single path whose action is
Assume the contribution to the integral is mainly near y = x + \eta where \eta \ll 1. The Green function is thus given by
The integrand decays to zero for \eta > \sqrt{2 b^{2} \delta } thus the integrand must be expanded to linear order in \delta and quadratic order in \eta in order to capture the lowest order behavior of the Green function. The expansion of the Green function yields the solvable integral
The solution for the normalization constant A is determined by solution of the leading term on the right hand side:
thus A = b \sqrt{ 2 \pi \delta}. The two additional integrals are given by
The Green function is given by
Thus, the Green function for the 1-D Gaussian Chain obeys
with initial condition G(x,0|x_{0},0)=\delta(x-x_{0}). For V(x)=0, the solution for the Green function is
which agrees with our previous derivation (Sec.~ref{chap:gaussian-explicit}).
The governing differential equation for the Gaussian Chain in 3D is written as an extension to our 1D equation as
with initial condition G(\vec{R},0|\vec{R}_{0},0)=\delta(\vec{R}-\vec{R}_{0}). The solution of this differential equation subject to a given potential V gives the equilibrium probability distribution for the Gaussian Chain model. For V(\vec{R})=0, the solution for the Green function for a Gaussian Chain in 3-D is