Ideal chain

An ideal chain (or freely-jointed chain) is the simplest model to describe polymers, such as nucleic acids and proteins. It only assumes a polymer as a random walk and neglects any kind of interactions among monomers. Although it is simple, its generality gives insight about the physics of polymers.

In this model, monomers are rigid rods of a fixed length l, and their orientation is completely independent of the orientations and positions of neighbouring monomers, to the extent that two monomers can co-exist at the same place. In some cases, the monomer has a physical interpretation, such as an amino acid in a polypeptide. In other cases, a monomer is simply a segment of the polymer that can be modeled as behaving as a discrete, freely jointed unit. If so, l is the Kuhn length. For example, chromatin is modeled as a polymer in which each monomer is a segment approximately 14-46 kbp in length.[1]

The model

N mers form the polymer, whose total unfolded length is:

${\displaystyle L=N\,l}$, where N is the number of mers.

In this very simple approach where no interactions between mers are considered, the energy of the polymer is taken to be independent of its shape, which means that at thermodynamic equilibrium, all of its shape configurations are equally likely to occur as the polymer fluctuates in time, according to the Maxwell–Boltzmann distribution.

Let us call ${\displaystyle {\vec {R}}}$ the total end to end vector of an ideal chain and ${\displaystyle {\vec {r}}_{1},\ldots ,{\vec {r}}_{N}}$ the vectors corresponding to individual mers. Those random vectors have components in the three directions of space. Most of the expressions given in this article assume that the number of mers N is large, so that the central limit theorem applies. The figure below shows a sketch of a (short) ideal chain.

The two ends of the chain are not coincident, but they fluctuate around each other, so that of course:

${\displaystyle \langle {\vec {R}}\rangle =\Sigma _{i=1}^{N}\langle {\vec {r}}_{i}\rangle ={\vec {0}}~}$

Throughout the article the ${\displaystyle \langle \rangle }$ brackets will be used to denote the mean (of values taken over time) of a random variable or a random vector, as above.

Since ${\displaystyle {\vec {r}}_{1},\ldots ,{\vec {r}}_{N}}$ are independent, it follows from the Central limit theorem that ${\displaystyle {\vec {R}}}$ is distributed according to a normal distribution (or gaussian distribution): precisely, in 3D, ${\displaystyle R_{x},R_{y},}$ and ${\displaystyle R_{z}}$ are distributed according to a normal distribution of mean 0 and of variance:

${\displaystyle \sigma ^{2}=\langle R_{x}^{2}\rangle -\langle R_{x}\rangle ^{2}=\langle R_{x}^{2}\rangle -0}$

${\displaystyle \langle R_{x}^{2}\rangle =\langle R_{y}^{2}\rangle =\langle R_{z}^{2}\rangle =N\,{\frac {l^{2}}{3}}}$

So that ${\displaystyle \langle {\vec {R^{2}}}\rangle =N\,l^{2}=L\,l~}$. The end to end vector of the chain is distributed according to the following probability density function:

${\displaystyle P({\vec {R}})=\left({\frac {3}{2\pi Nl^{2}}}\right)^{3/2}e^{-{\frac {3{\vec {R}}^{2}}{2Nl^{2}}}}}$

The average end-to-end distance of the polymer is:

${\displaystyle {\sqrt {\langle {\vec {R^{2}}}\rangle }}={\sqrt {N}}\,l={\sqrt {L\,l}}~}$

A quantity frequently used in polymer physics is the radius of gyration:

${\displaystyle {\mathit {R}}_{G}={\frac {{\sqrt {N}}\,l}{{\sqrt {6}}\ }}}$

It is worth noting that the above average end-to-end distance, which in the case of this simple model is also the typical amplitude of the system's fluctuations, becomes negligible compared to the total unfolded length of the polymer ${\displaystyle N\,l}$ at the thermodynamic limit. This result is a general property of statistical systems.

Mathematical remark: the rigorous demonstration of the expression of the density of probability ${\displaystyle P({\vec {R}})}$ is not as direct as it appears above: from the application of the usual (1D) central limit theorem one can deduce that ${\displaystyle R_{x}}$, ${\displaystyle R_{y}}$ and ${\displaystyle R_{z}}$ are distributed according to a centered normal distribution of variance ${\displaystyle N\,l^{2}/3}$. Then, the expression given above for ${\displaystyle P({\vec {R}})}$ is not the only one that is compatible with such distribution for ${\displaystyle R_{x}}$, ${\displaystyle R_{y}}$ and ${\displaystyle R_{z}}$. However, since the components of the vectors ${\displaystyle {\vec {r}}_{1},\ldots ,{\vec {r}}_{N}}$ are uncorrelated for the random walk we are considering, it follows that ${\displaystyle R_{x}}$, ${\displaystyle R_{y}}$ and ${\displaystyle R_{z}}$ are also uncorrelated. This additional condition can only be fulfilled if ${\displaystyle {\vec {R}}}$ is distributed according to ${\displaystyle P({\vec {R}})}$. Alternatively, this result can also be demonstrated by applying a multidimensional generalization of the central limit theorem, or through symmetry arguments.

Generality of the model

While the elementary model described above is totally unadapted to the description of real-world polymers at the microscopic scale, it does show some relevance at the macroscopic scale in the case of a polymer in solution whose monomers form an ideal mix with the solvent (in which case, the interactions between monomer and monomer, solvent molecule and solvent molecule, and between monomer and solvent are identical, and the system's energy can be considered constant, validating the hypotheses of the model).

The relevancy of the model is, however, limited, even at the macroscopic scale, by the fact that it does not consider any excluded volume for monomers (or, to speak in chemical terms, that it neglects steric effects).

Other fluctuating polymer models that consider no interaction between monomers and no excluded volume, like the worm-like chain model, are all asymptotically convergent toward this model at the thermodynamic limit. For purpose of this analogy a Kuhn segment is introduced, corresponding to the equivalent monomer length to be considered in the analogous ideal chain. The number of Kuhn segments to be considered in the analogous ideal chain is equal to the total unfolded length of the polymer divided by the length of a Kuhn segment.

Entropic elasticity of an ideal chain

If the two free ends of an ideal chain are attached to some kind of micro-manipulation device, then the device experiences a force exerted by the polymer. The ideal chain's energy is constant, and thus its time-average, the internal energy, is also constant, which means that this force necessarily stems from a purely entropic effect.

This entropic force is very similar to the pressure experienced by the walls of a box containing an ideal gas. The internal energy of an ideal gas depends only on its temperature, and not on the volume of its containing box, so it is not an energy effect that tends to increase the volume of the box like gas pressure does. This implies that the pressure of an ideal gas has a purely entropic origin.

What is the microscopic origin of such an entropic force or pressure? The most general answer is that the effect of thermal fluctuations tends to bring a thermodynamic system toward a macroscopic state that corresponds to a maximum in the number of microscopic states (or micro-states) that are compatible with this macroscopic state. In other words, thermal fluctuations tend to bring a system toward its macroscopic state of maximum entropy.

What does this mean in the case of the ideal chain? First, for our ideal chain, a microscopic state is characterized by the superposition of the states ${\displaystyle {\vec {r}}_{i}}$ of each individual monomer (with i varying from 1 to N). In its solvent, the ideal chain is constantly subject to shocks from moving solvent molecules, and each of these shocks sends the system from its current microscopic state to another, very similar microscopic state. For an ideal polymer, as will be shown below, there are more microscopic states compatible with a short end-to-end distance than there are microscopic states compatible with a large end-to-end distance. Thus, for an ideal chain, maximizing its entropy means reducing the distance between its two free ends. Consequently, a force that tends to collapse the chain is exerted by the ideal chain between its two free ends.

In this section, the mean of this force will be derived. The generality of the expression obtained at the thermodynamic limit will then be discussed.

Ideal chain under length constraint

The case of an ideal chain whose two ends are attached to fixed points will be considered in this sub-section. The vector ${\displaystyle {\vec {R}}}$ joining these two points characterizes the macroscopic state (or macro-state) of the ideal chain. Each macro-state corresponds a certain number of micro-states, that we will call ${\displaystyle \Omega ({\vec {R}})}$ (micro-states are defined in the introduction to this section). Since the ideal chain's energy is constant, each of these micro-states is equally likely to occur. The entropy associated to a macro-state is thus equal to:

${\displaystyle S({\vec {R}})=k_{B}\log(\Omega ({\vec {R}}))}$, where ${\displaystyle k_{B}}$ is Boltzmann's constant

The above expression gives the absolute (quantum) entropy of the system. A precise determination of ${\displaystyle \Omega ({\vec {R}})}$ would require a quantum model for the ideal chain, which is beyond the scope of this article. However, we have already calculated the probability density ${\displaystyle P({\vec {R}})}$ associated with the end-to-end vector of the unconstrained ideal chain, above. Since all micro-states of the ideal chain are equally likely to occur, ${\displaystyle P({\vec {R}})}$ is proportional to ${\displaystyle \Omega ({\vec {R}})}$. This leads to the following expression for the classical (relative) entropy of the ideal chain:

${\displaystyle S({\vec {R}})=k_{B}\log(P({\vec {R}}))+C_{st}}$,

where ${\displaystyle C_{st}}$ is a fixed constant. Let us call ${\displaystyle {\vec {F}}}$ the force exerted by the chain on the point to which its end is attached. From the above expression of the entropy, we can deduce an expression of this force. Suppose that, instead of being fixed, the positions of the two ends of the ideal chain are now controlled by an operator. The operator controls the evolution of the end to end vector ${\displaystyle {\vec {R}}}$. If the operator changes ${\displaystyle {\vec {R}}}$ by a tiny amount ${\displaystyle {\vec {dR}}}$, then the variation of internal energy of the chain is zero, since the energy of the chain is constant. This condition can be written as:

${\displaystyle 0=dU=\delta W+\delta Q~}$

${\displaystyle \delta W}$ is defined as the elementary amount of mechanical work transferred by the operator to the ideal chain, and ${\displaystyle \delta Q}$ is defined as the elementary amount of heat transferred by the solvent to the ideal chain. Now, if we assume that the transformation imposed by the operator on the system is quasistatic (i.e., infinitely slow), then the system's transformation will be time-reversible, and we can assume that during its passage from macro-state ${\displaystyle {\vec {R}}}$ to macro-state ${\displaystyle {\vec {R}}+{\vec {dR}}}$, the system passes through a series of thermodynamic equilibrium macro-states. This has two consequences:

• first, the amount of heat received by the system during the transformation can be tied to the variation of its entropy:

${\displaystyle \delta Q=TdS~}$, where T is the temperature of the chain.

• second, in order for the transformation to remain infinitely slow, the mean force exerted by the operator on the end points of the chain must balance the mean force exerted by the chain on its end points. Calling ${\displaystyle {\vec {f}}_{op}}$ the force exerted by the operator and ${\displaystyle {\vec {f}}}$ the force exerted by the chain, we have:

${\displaystyle \delta W=\langle {\vec {f}}_{op}\rangle \cdot {\vec {dR}}=-\langle {\vec {f}}\rangle \cdot {\vec {dR}}~}$

We are thus led to:

${\displaystyle \langle {\vec {f}}\rangle =T{\frac {dS}{\vec {dR}}}={\frac {k_{B}T}{P({\vec {R}})}}{\frac {dP({\vec {R}})}{\vec {dR}}}~}$

${\displaystyle \langle {\vec {f}}\rangle =-k_{B}T{\frac {3{\vec {R}}}{Nl^{2}}}~}$

The above equation is the equation of state of the ideal chain. Since the expression depends on the central limit theorem, it is only exact in the limit of polymers containing a large number of monomers (that is, the thermodynamic limit). It is also only valid for small end-to-end distances, relative to the overall polymer contour length, where the behavior is like a hookean spring. Behavior over larger force ranges can be modeled using a canonical ensemble treatment identical to magnetization of paramagnetic spins. For the arbitrary forces the extension-force dependence will be given by Langevin function ${\displaystyle {\mathcal {L}}}$:

${\displaystyle {\frac {R}{Nl}}=\coth \left({\frac {fl}{k_{B}T}}\right)-{\frac {k_{B}T}{fl}}={\mathcal {L}}\left({\frac {fl}{k_{B}T}}\right),}$

where the extension is ${\displaystyle R=|{\vec {R}}|}$.

For the arbitrary extensions the force-extension dependence can be approximated by:[2]

${\displaystyle {\frac {fl}{k_{B}T}}={\mathcal {L}}^{-1}\left({\frac {R}{Nl}}\right)\approx 3{\frac {R}{Nl}}+{\frac {1}{5}}\left({\frac {R}{Nl}}\right)^{2}\sin \left({\frac {7R}{2Nl}}\right)+{\frac {\left({\frac {R}{Nl}}\right)^{3}}{1-{\frac {R}{Nl}}}}}$,

where ${\displaystyle {\mathcal {L}}^{-1}}$ is the inverse Langevin function, N is the number of bonds[3] in the molecule (therefore if the molecule has N bonds it has N+1 monomers making up the molecule.).

Finally, the model can be extended to even larger force ranges by inclusion of a stretch modulus along the polymer contour length. That is, by allowing the length of each unit of the chain to respond elastically to the applied force.[4]

Ideal polymer exchanging length with a reservoir

Throughout this sub-section, as in the previous one, the two ends of the polymer are attached to a micro-manipulation device. This time, however, the device does not maintain the two ends of the ideal chain in a fixed position, but rather it maintains a constant pulling force ${\displaystyle {\vec {f}}_{op}}$ on the ideal chain. In this case the two ends of the polymer fluctuate around a mean position ${\displaystyle \langle {\vec {R}}\rangle }$. The ideal chain reacts with a constant opposite force ${\displaystyle {\vec {f}}=-{\vec {f}}_{op}}$

For an ideal chain exchanging length with a reservoir, a macro-state of the system is characterized by the vector ${\displaystyle {\vec {f}}}$.

The change between an ideal chain of fixed length and an ideal chain in contact with a length reservoir is very much akin to the change between the micro-canonical ensemble and the canonical ensemble (see the Statistical mechanics article about this). The change is from a state where a fixed value is imposed on a certain parameter, to a state where the system is left free to exchange this parameter with the outside. The parameter in question is energy for the microcanonical and canonical descriptions, whereas in the case of the ideal chain the parameter is the length of the ideal chain.

As in the micro-canonical and canonical ensembles, the two descriptions of the ideal chain differ only in the way they treat the system's fluctuations. They are thus equivalent at the thermodynamic limit. The equation of state of the ideal chain remains the same, except that ${\displaystyle {\vec {R}}}$ is now subject to fluctuations:

${\displaystyle {\vec {f}}=-k_{B}T{\frac {3\langle {\vec {R}}\rangle }{Nl^{2}}}~}$.

Ideal chain under a constant force constraint - calculation

A diagram of an ideal chain constrained by a constant force.

Consider a freely jointed chain of N bonds of length ${\displaystyle l}$ subject to a constant elongational force f applied to its ends along the z axis and an environment temperature ${\displaystyle T}$. An example could be a chain with two opposite charges +q and -q at its ends in a constant electric field ${\displaystyle {\vec {E}}}$ applied along the ${\displaystyle z}$ axis as sketched in the figure on the right. If the direct Coulomb interaction between the charges is ignored, then there is a constant force ${\displaystyle {\vec {f}}}$ at the two ends.

Different chain conformations are not equally likely, because they correspond to different energy of the chain in the external electric field.

${\displaystyle U=-q{\vec {E}}\cdot {\vec {R}}=-{\vec {f}}\cdot {\vec {R}}=-fR_{z}}$

Thus, different chain conformation have different statistical Boltzmann factors ${\displaystyle exp(-U/k_{B}T)}$.[3]

The partition function is:

${\displaystyle Z=\sum _{states}exp(-U/k_{B}T)=\sum _{states}exp({fR_{z} \over k_{B}T})}$

Every monomer connection in the chain is characterize by a vector ${\displaystyle {\vec {r_{i}}}}$ of length ${\displaystyle l}$ and angles ${\displaystyle \theta _{i},\varphi _{i}}$ in the spherical coordinate system. The end-to-end vector can be represented as: ${\displaystyle R_{z}=\sum _{i=1}^{N}l{\text{ }}\cos \theta _{i}}$. Therefore:

${\displaystyle {\begin{aligned}Z=&\int \exp({fl \over k_{B}T}\sum _{i=1}^{N}\cos \theta _{i})\prod _{i=1}^{N}\sin \theta _{i}d\theta _{i}d\varphi _{i}=\\&\left[\int _{0}^{\pi }2\pi {\text{ }}\sin \theta _{i}{\text{ }}\exp({fl \over k_{B}T}\cos \theta _{i})d\theta _{i}\right]^{N}=\\&\left[{2\pi \over fl/(k_{B}T)}[\exp({fl \over k_{B}T})-\exp(-{fl \over k_{B}T})]\right]^{N}=\\&\left[{4\pi {\text{ }}\sinh(fl/(K_{B}T)) \over fl/(k_{B}T)}\right]^{N}\end{aligned}}}$

The Gibbs free energy G can be directly calculated from the partition function:

${\displaystyle G(T,f,N)=-k_{B}T{\text{ }}ln{\text{ }}Z(T,f,N)=-Nk_{B}T[ln(4\pi {\text{ }}sinh({fl \over k_{B}T}))-ln({fl \over k_{B}T})]}$

The Gibbs free energy is used here because the ensemble of chains corresponds to constant temperature ${\displaystyle T}$ and constant force ${\displaystyle f}$ (analogous to the isothermal-isobaric ensemble, which has constant temperature and pressure).

The average end-to-end distance corresponding to a given force can be obtained as the derivative of the free energy:

${\displaystyle <R>=-{\partial G \over \partial f}=Nl[coth({fl \over k_{B}T})-{1 \over fl/(k_{B}T)}]}$

This expression is the Langevin function ${\displaystyle {\mathcal {L}}}$, also mentioned in previous paragraphs:

The average distance ${\displaystyle \left\langle {\vec {R}}\right\rangle }$ of the chain as a function of ${\displaystyle \alpha }$.

${\displaystyle {\mathcal {L}}(\alpha )=coth(\alpha )-{1 \over \alpha }}$ where, ${\displaystyle \alpha ={fl \over k_{B}T}}$.

For small relative elongations (${\displaystyle \left\langle R\right\rangle \ll R_{max}=lN}$) the dependence is approximately linear,

${\displaystyle {\mathcal {L}}(\alpha )\cong {\alpha \over 3}}$ for ${\displaystyle \alpha \ll 1}$

and follows Hooke's law as shown in previous paragraphs:

${\displaystyle {\vec {f}}=K_{B}T{3\left\langle {\vec {R}}\right\rangle \over Nl^{2}}}$

See also

• Polymer
• Worm-like chain, a more complex polymer model
• Kuhn length
• Coil-globule transition

References

1. Rippe, Karsten (2001). "Making contacts on a nucleic acid polymer". Trends in Biochemical Sciences. 26 (12): 733–740. doi:10.1016/S0968-0004(01)01978-8.
2. Petrosyan, R. (2016). "Improved approximations for some polymer extension models". Rehol Acta. 56: 21–26. arXiv:1606.02519. doi:10.1007/s00397-016-0977-9.
3. Polymer Physics ISBN 019852059-X, 76, Rubinstein
4. Smith, SB; Finzi, L; Bustamante, C (1992). "Direct mechanical measurements of the elasticity of single DNA molecules by using magnetic beads". Science. 258 (5085): 1122–6. Bibcode:1992Sci...258.1122S. doi:10.1126/science.1439819. PMID 1439819.