Movatterモバイル変換


[0]ホーム

URL:


Jump to content
WikipediaThe Free Encyclopedia
Search

Ideal chain

From Wikipedia, the free encyclopedia
Molecular model for describing polymers
icon
This articleneeds additional citations forverification. Please helpimprove this article byadding citations to reliable sources. Unsourced material may be challenged and removed.
Find sources: "Ideal chain" – news ·newspapers ·books ·scholar ·JSTOR
(April 2025) (Learn how and when to remove this message)

Anideal chain (orfreely-jointed chain) is the simplestmodel inpolymer chemistry to describepolymers, such asnucleic acids andproteins. It assumes that themonomers in a polymer are located at the steps of a hypotheticalrandom walker that does not remember its previous steps. By neglecting interactions amongmonomers, this model assumes that two (or more) monomers can occupy the same location. Although it is simple, its generality gives insight about thephysics of polymers.

In this model, monomers are rigid rods of a fixed lengthl, and their orientation is completely independent of the orientations and positions of neighbouring monomers. In some cases, the monomer has a physical interpretation, such as anamino acid in apolypeptide. 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 theKuhn length. For example,chromatin is modeled as a polymer in which each monomer is a segment approximately14–46 kbp in length.[1]

Model

[edit]

Nmers form the polymer, whose total unfolded length is:L=Nl,{\displaystyle L=N\,l,} whereN 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 atthermodynamic equilibrium, all of its shape configurations are equally likely to occur as the polymer fluctuates in time, according to theMaxwell–Boltzmann distribution.

Let us callR{\displaystyle {\vec {R}}} the total end to end vector of an ideal chain andr1,,rN{\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 mersN is large, so that thecentral 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:R=i=1Nri=0 {\displaystyle \left\langle {\vec {R}}\right\rangle =\sum _{i=1}^{N}\left\langle {\vec {r}}_{i}\right\rangle ={\vec {0}}~}

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

Sincer1,,rN{\displaystyle {\vec {r}}_{1},\ldots ,{\vec {r}}_{N}} areindependent, it follows from thecentral limit theorem thatR{\displaystyle {\vec {R}}} is distributed according to anormal distribution (or gaussian distribution): precisely, in 3D,Rx,Ry,{\displaystyle R_{x},R_{y},} andRz{\displaystyle R_{z}} are distributed according to anormal distribution ofmean 0 and ofvariance:σ2=Rx2Rx2=Rx20{\displaystyle \sigma ^{2}=\langle R_{x}^{2}\rangle -\langle R_{x}\rangle ^{2}=\langle R_{x}^{2}\rangle -0}Rx2=Ry2=Rz2=Nl23{\displaystyle \left\langle R_{x}^{2}\right\rangle =\left\langle R_{y}^{2}\right\rangle =\left\langle R_{z}^{2}\right\rangle =N\,{\frac {l^{2}}{3}}}

So thatR2=Nl2=Ll {\displaystyle \langle {R^{2}}\rangle =N\,l^{2}=L\,l~}. The end to end vector of the chain is distributed according to the followingprobability density function:P(R)=(32πNl2)3/2e3|R|22Nl2{\displaystyle P({\vec {R}})=\left({\frac {3}{2\pi Nl^{2}}}\right)^{3/2}e^{-{\frac {3\left|{\vec {R}}\right|^{2}}{2Nl^{2}}}}}

The average end-to-end distance of the polymer is:R2=Nl=Ll {\displaystyle {\sqrt {\left\langle {R^{2}}\right\rangle }}={\sqrt {N}}\,l={\sqrt {L\,l}}~}

A quantity frequently used in polymer physics is theradius of gyration:RG=Nl6 {\displaystyle \langle {\mathit {R}}_{G}\rangle ={\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 polymerNl{\displaystyle N\,l} at thethermodynamic limit. This result is a general property of statistical systems.

Mathematical remark: the rigorous demonstration of the expression of the density of probabilityP(R){\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 thatRx{\displaystyle R_{x}},Ry{\displaystyle R_{y}} andRz{\displaystyle R_{z}} are distributed according to a centerednormal distribution of varianceNl2/3{\displaystyle N\,l^{2}/3}. Then, the expression given above forP(R){\displaystyle P({\vec {R}})} is not the only one that is compatible with such distribution forRx{\displaystyle R_{x}},Ry{\displaystyle R_{y}} andRz{\displaystyle R_{z}}. However, since the components of the vectorsr1,,rN{\displaystyle {\vec {r}}_{1},\ldots ,{\vec {r}}_{N}} areuncorrelated for the random walk we are considering, it follows thatRx{\displaystyle R_{x}},Ry{\displaystyle R_{y}} andRz{\displaystyle R_{z}} are alsouncorrelated. This additional condition can only be fulfilled ifR{\displaystyle {\vec {R}}} is distributed according toP(R){\displaystyle P({\vec {R}})}. Alternatively, this result can also be demonstrated by applying a multidimensional generalization of thecentral limit theorem, or throughsymmetry arguments.

Generality of the model

[edit]

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 neglectssteric effects). Since theN mers are of a rigid, fixed length, the model also does not consider bond stretching, though it can be extended to do so.[2]

Other fluctuating polymer models that consider no interaction between monomers and no excluded volume, like theworm-like chain model, are all asymptotically convergent toward this model at thethermodynamic limit. For purpose of this analogy aKuhn 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

[edit]

If the two free ends of an ideal chain are pulled apart by some sort of device, then the device experiences a force exerted by the polymer. As the ideal chain is stretched, its energy remains constant, and its time-average, orinternal energy, also remains constant, which means that this force necessarily stems from a purelyentropic effect.

Thisentropic force is very similar to the pressure experienced by the walls of a box containing anideal gas. Theinternal energy of anideal gas depends only on its temperature, and not on the volume of its containing box, so it is not anenergy effect that tends to increase the volume of the box like gaspressure does. This implies that thepressure of an ideal gas has a purelyentropic origin.

What is the microscopic origin of such anentropic force or pressure? The most general answer is that the effect of thermal fluctuations tends to bring athermodynamic 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 maximumentropy.

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 statesri{\displaystyle {\vec {r}}_{i}} of each individual monomer (withi varying from1 toN). 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 itsentropy 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, themean of this force will be derived. The generality of the expression obtained at thethermodynamic limit will then be discussed.

Ideal chain under length constraint

[edit]

The case of an ideal chain whose two ends are attached to fixed points will be considered in this sub-section. The vectorR{\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Ω(R){\displaystyle \Omega ({\vec {R}})} (micro-states are defined in the introduction to this section). Since the ideal chain'senergy is constant, each of these micro-states is equally likely to occur. Theentropy associated to a macro-state is thus equal to:S(R)=kBlog(Ω(R)),{\displaystyle S({\vec {R}})=k_{\text{B}}\log(\Omega ({\vec {R}})),}wherekB{\displaystyle k_{\text{B}}} is theBoltzmann constant

The above expression gives the absolute (quantum)entropy of the system. A precise determination ofΩ(R){\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 densityP(R){\displaystyle P({\vec {R}})} associated with theend-to-end vector of theunconstrained ideal chain, above. Since all micro-states of the ideal chain are equally likely to occur,P(R){\displaystyle P({\vec {R}})} is proportional toΩ(R){\displaystyle \Omega ({\vec {R}})}. This leads to the following expression for the classical (relative)entropy of the ideal chain:S(R)=kBlog(P(R))+Cst,{\displaystyle S({\vec {R}})=k_{\text{B}}\log(P({\vec {R}}))+C_{st},}whereCst{\displaystyle C_{st}} is a fixed constant. Let us callF{\displaystyle {\vec {F}}} the force exerted by the chain on the point to which its end is attached. From the above expression of theentropy, 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 vectorR{\displaystyle {\vec {R}}}. If the operator changesR{\displaystyle {\vec {R}}} by a tiny amountdR{\displaystyle d{\vec {R}}}, then the variation ofinternal energy of the chain is zero, since theenergy of the chain is constant. This condition can be written as:0=dU=δW+δQ {\displaystyle 0=dU=\delta W+\delta Q~}

δW{\displaystyle \delta W} is defined as the elementary amount ofmechanical work transferred by the operator to the ideal chain, andδQ{\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-stateR{\displaystyle {\vec {R}}} to macro-stateR+dR{\displaystyle {\vec {R}}+d{\vec {R}}}, the system passes through a series ofthermodynamic equilibrium macro-states. This has two consequences:

  1. first, the amount ofheat received by the system during the transformation can be tied to the variation of itsentropy:δQ=TdS,{\displaystyle \delta Q=T\,dS,} whereT is the temperature of the chain.
  2. second, in order for the transformation to remain infinitely slow, themean force exerted by the operator on the end points of the chain must balance themean force exerted by the chain on its end points. Callingfop{\displaystyle {\vec {f}}_{\text{op}}} the force exerted by the operator andf{\displaystyle {\vec {f}}} the force exerted by the chain, we have:δW=fopdR=fdR{\displaystyle \delta W=\langle {\vec {f}}_{\text{op}}\rangle \cdot d{\vec {R}}=-\langle {\vec {f}}\rangle \cdot d{\vec {R}}}

We are thus led to:f=TdSdR=kBTP(R)dP(R)dR{\displaystyle \langle {\vec {f}}\rangle =T{\frac {dS}{d{\vec {R}}}}={\frac {k_{\text{B}}T}{P({\vec {R}})}}{\frac {dP({\vec {R}})}{d{\vec {R}}}}}f=kBT3RNl2{\displaystyle \langle {\vec {f}}\rangle =-k_{\text{B}}T{\frac {3{\vec {R}}}{Nl^{2}}}}

The above equation is theequation of state of the ideal chain. Since the expression depends on thecentral limit theorem, it is only exact in the limit of polymers containing a large number of monomers (that is, thethermodynamic limit). It is also only valid for small end-to-end distances, relative to the overall polymercontour length, where the behavior is like a hookean spring. Behavior over larger force ranges can be modeled using acanonical ensemble treatment identical to magnetization of paramagnetic spins. For the arbitrary forces the extension-force dependence will be given byLangevin functionL{\displaystyle {\mathcal {L}}}:RNl=coth(flkBT)kBTfl=L(flkBT),{\displaystyle {\frac {R}{Nl}}=\coth \left({\frac {fl}{k_{\text{B}}T}}\right)-{\frac {k_{\text{B}}T}{fl}}={\mathcal {L}}\left({\frac {fl}{k_{\text{B}}T}}\right),}where the extension isR=|R|{\displaystyle R=|{\vec {R}}|}.

For the arbitrary extensions the force-extension dependence can be approximated by:[3]flkBT=L1(RNl)3RNl+15(RNl)2sin(7R2Nl)+(RNl)31RNl,{\displaystyle {\frac {fl}{k_{\text{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}}}},}whereL1{\displaystyle {\mathcal {L}}^{-1}} is the inverseLangevin function,N is the number of bonds[4] in the molecule (therefore if the molecule hasN bonds it hasN + 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.[5]

Ideal polymer exchanging length with a reservoir

[edit]

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 forcefop{\displaystyle {\vec {f}}_{\text{op}}} on the ideal chain. In this case the two ends of the polymer fluctuate around amean positionR{\displaystyle \langle {\vec {R}}\rangle }. The ideal chain reacts with a constant opposite forcef=fop{\displaystyle {\vec {f}}=-{\vec {f}}_{\text{op}}}.

For an ideal chain exchanging length with a reservoir, a macro-state of the system is characterized by the vectorf{\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 theStatistical mechanics article about this).[6] 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 thethermodynamic limit. Theequation of state of the ideal chain remains the same, except thatR{\displaystyle {\vec {R}}} is now subject to fluctuations:f=kBT3RNl2 .{\displaystyle {\vec {f}}=-k_{\text{B}}T{\frac {3\langle {\vec {R}}\rangle }{Nl^{2}}}~.}

Ideal chain under a constant force constraint – calculation

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

Consider a freely jointed chain ofN{\displaystyle N} bonds of lengthl{\displaystyle l} subject to a constant elongational forcef{\displaystyle f} applied to its ends along thez{\displaystyle z} axis and an environment temperatureT{\displaystyle T}. An example could be a chain with two opposite charges+q{\displaystyle +q} andq{\displaystyle -q} at its ends in a constantelectric fieldE{\displaystyle {\vec {E}}} applied along thez{\displaystyle z} axis as sketched in the figure on the right. If the directCoulomb interaction between the charges is ignored, then there is a constant forcef{\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.U=qER=fR=fRz{\displaystyle U=-q{\vec {E}}\cdot {\vec {R}}=-{\vec {f}}\cdot {\vec {R}}=-fR_{z}}

Thus, different chain conformation have different statisticalBoltzmann factorsexp(U/kBT){\displaystyle \exp(-U/k_{\text{B}}T)}.[4]

Thepartition function is:Z=statesexp(U/kBT)=statesexp(fRzkBT){\displaystyle Z=\sum _{\text{states}}\exp(-U/k_{\text{B}}T)=\sum _{\text{states}}\exp \left({\frac {fR_{z}}{k_{\text{B}}T}}\right)}

Everymonomer connection in the chain is characterized by a vectorri{\displaystyle {\vec {r}}_{i}} of lengthl{\displaystyle l} and anglesθi,φi{\displaystyle \theta _{i},\varphi _{i}} in thespherical coordinate system. The end-to-end vector can be represented as:Rz=i=1Nlcosθi{\textstyle R_{z}=\sum _{i=1}^{N}l\cos \theta _{i}}. Therefore:

Z=exp(flkBTi=1Ncosθi)i=1Nsinθidθidφi=[0π2π sinθi exp(flkBTcosθi)dθi]N=[2πfl/(kBT)(exp(flkBT)exp(flkBT))]N=[4πsinh(fl/(kBT))fl/(kBT)]N{\displaystyle {\begin{aligned}Z&=\int \exp \left({\frac {fl}{k_{\text{B}}T}}\sum _{i=1}^{N}\cos \theta _{i}\right)\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 \left({\frac {fl}{k_{\text{B}}T}}\cos \theta _{i}\right)\,d\theta _{i}\right]^{N}\\&=\left[{\frac {2\pi }{fl/(k_{\text{B}}T)}}\left(\exp \left({fl \over k_{\text{B}}T}\right)-\exp \left(-{\frac {fl}{k_{\text{B}}T}}\right)\right)\right]^{N}\\&=\left[{4\pi \sinh(fl/(k_{\text{B}}T)) \over fl/(k_{\text{B}}T)}\right]^{N}\end{aligned}}}

TheGibbs free energy G can be directly calculated from the partition function:G(T,f,N)=kBTlnZ(T,f,N)=NkBT[ln(4πsinh(flkBT))ln(flkBT)]{\displaystyle G(T,f,N)=-k_{\text{B}}T\,\ln Z(T,f,N)=-Nk_{\text{B}}T\left[\ln \left(4\pi \sinh \left({\frac {fl}{k_{\text{B}}T}}\right)\right)-\ln \left({\frac {fl}{k_{\text{B}}T}}\right)\right]}

The Gibbs free energy is used here because the ensemble of chains corresponds to constant temperatureT{\displaystyle T} and constant forcef{\displaystyle f} (analogous to theisothermal–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:R=Gf=Nl[coth(flkBT)kBTfl]{\displaystyle \langle R\rangle =-{\frac {\partial G}{\partial f}}=Nl\left[\coth \left({\frac {fl}{k_{\text{B}}T}}\right)-{\frac {k_{\text{B}}T}{fl}}\right]}

This expression is theLangevin functionL{\displaystyle {\mathcal {L}}}, also mentioned in previous paragraphs:

The average distanceR{\displaystyle \langle {\vec {R}}\rangle } of the chain as a function ofα{\displaystyle \alpha }.

L(α)=coth(α)1α{\displaystyle {\mathcal {L}}(\alpha )=\coth(\alpha )-{1 \over \alpha }}whereα=flkBT{\displaystyle \alpha ={\frac {fl}{k_{\text{B}}T}}}.

For small relative elongations (RRmax=lN{\displaystyle \langle R\rangle \ll R_{\text{max}}=lN}) the dependence is approximately linear,L(α)α3 for α1{\displaystyle {\mathcal {L}}(\alpha )\cong {\frac {\alpha }{3}}\qquad {\text{ for }}\alpha \ll 1}and followsHooke's law as shown in previous paragraphs:f=kBT3RNl2{\displaystyle {\vec {f}}=k_{\text{B}}T{\frac {3\langle {\vec {R}}\rangle }{Nl^{2}}}}

See also

[edit]

References

[edit]
  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. ^Buche, M.R.; Silberstein, M.N.; Grutzik, S.J. (2022). "Freely jointed chains with extensible links".Phys. Rev. E.106 024502.arXiv:2203.05421.doi:10.1103/PhysRevE.106.024502.
  3. ^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.
  4. ^abPolymer PhysicsISBN 019852059-X, 76, Rubinstein
  5. ^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.
  6. ^Buche, M.R.; Silberstein, M.N. (2020). "Statistical mechanical constitutive theory of polymer networks: The inextricable links between distribution, behavior, and ensemble".Phys. Rev. E.102 012501.arXiv:2004.07874.doi:10.1103/PhysRevE.102.012501.
Retrieved from "https://en.wikipedia.org/w/index.php?title=Ideal_chain&oldid=1322916005"
Categories:
Hidden categories:

[8]ページ先頭

©2009-2025 Movatter.jp