被拉伸链的末端距

原链接不详。原作者：Tim St Pierre

自由连接链被拉伸，末端距与力之间的关系为朗之万函数。

考虑自由连接链，链节长度为$a$，矢量$\vec{a}$表示链节的取向和长度。取向与$x$轴夹角在$\theta$到$\theta+\delta \theta$之间的链节数目正比于球上面积

\begin{equation*} \delta A=2\pi a^2 \sin \theta \delta \theta \end{equation*}

即链节取向与$x$轴夹角为$\theta$的概率为

\begin{equation*} P(\theta)d\theta=C2\pi a^2 \sin \theta d\theta \end{equation*}

其中$C$为归一化常数，由下式定出

\begin{equation*} \int_0^{\theta} P(\theta)d\theta=1 \end{equation*}

于是有

\begin{equation*} P(\theta)=\frac{1}{2} \sin \theta d\theta \end{equation*}

键的取向

如果链被拉伸，不失一般性，假设力沿$x$轴方向。链节不同的取向对应的势能不相等。平均需要$2aF$的势能使沿$x$方向的链节变得沿$x$负方向。

链被拉伸

链节的取向势能为

\begin{equation*} V=-Fa\cos \theta=\vec{F} \cdot \vec{a} \end{equation*}

链节取向与与$x$轴夹角为$\theta$的概率正比于$\exp\left (-\frac{V}{k_BT}\right )$，即

\begin{equation*} P(\theta)d\theta=\left (\frac{1}{2} \sin \theta d\theta\right )\exp\left (-\frac{V}{k_BT}\right ) \end{equation*}

所以$\vec{a}$的$x$分量的平均值为

\begin{equation*} \begin{split} \langle a_x\rangle &=\frac{\int_0^{\pi}a\cos \theta \left (\frac{1}{2} \sin \theta d\theta\right )\exp\left (-\frac{V}{k_BT}\right )}{\int_0^{\pi}\left (\frac{1}{2} \sin \theta d\theta\right )\exp\left (-\frac{V}{k_BT}\right )}\\ &=a\left [\coth\left (\frac{Fa}{k_BT}\right )- \frac{k_BT}{Fa} \right ] \end{split} \end{equation*}

如果链长为$N$，则平均末端距的$x$分量为

\begin{equation*} \langle r_x\rangle=Na\left [\coth\left (\frac{Fa}{k_BT}\right )- \frac{k_BT}{Fa} \right ] \end{equation*}

链在$y$和$z$方向上没有被拉伸，因此

\begin{equation*} \langle r_y\rangle=\langle r_z\rangle=0 \end{equation*}

于是链的平均末端距为

\begin{equation*} \langle r\rangle=Na\left [\coth\left (\frac{Fa}{k_BT}\right )- \frac{k_BT}{Fa} \right ]=Na\mathcal {L}\left(\frac{Fa}{k_BT}\right ) \end{equation*}

其中$\mathcal {L}\left(\frac{Fa}{k_BT}\right )$为朗之万函数

对于双曲余切函数$\coth x=\frac{e^x+e^{-x}}{e^x+e^{-x}}$，展开成级数，为如下形式

\begin{equation*} \coth x=\frac{e^x+e^{-x}}{e^x+e^{-x}}=\frac{1}{x}+\frac{x}{3}+\cdots \end{equation*}

当拉伸力比较小时，$Fa\ll k_BT$，链末端距为

\begin{equation*} \langle r\rangle=\frac{NFa}{3k_BT} \end{equation*}

上式可写为

\begin{equation*} F=\frac{3 k_BT}{Na}\langle r\rangle \end{equation*}

此正是胡克定律。

如果$Fa\gt k_BT$，末端距必须采用朗之万函数形式。

末端距与拉伸力

交换上图坐标系，有下图

力与末端距

下图为实验结果，突变膜蛋白G241C Mutant的力谱（Science 2000, 288, 143）

力谱