[论文解读] Statistical Mechanics of Semiflexible Chains: A Meanfield Variational Approach

D. Thirumalai, Bae‐Yeun Ha|ArXiv.org|May 20, 1997

Protein Structure and DynamicsBiochemistry, Genetics and Molecular Biology参考文献 7被引用 21

一句话总结

本文提出一种基于平均场变分法的半柔性聚合物链模型，通过用全局平均约束替代局部切向矢量约束，实现可处理的统计力学计算。该方法能准确预测链端距离分布及拉伸力下的弹性响应，与DNA的模拟和实验结果定量吻合，但在强向列场中低估了持久长度，原因在于未能捕捉对称性破缺效应。

ABSTRACT

We describe a simple meanfield variational approach to study a number of properties of intrinsically stiff chains which are appropriate models for a large class of biopolymers. We present the calculation of the distribution of end-to-end distance and the elastic response of stiff chains under tension using this approach. In the former example we find that the simple expression almost quantitatively fits the results of computer simulation. For the case of the stiff chain under tension we recover analytically all the known limits. We obtain quantitative agreement with recent experiments on the stretching of DNA. The limitations of our approach are also discussed.

研究动机与目标

• 开发适用于DNA和生物大分子等本征刚性聚合物的可处理统计力学框架。
• 解决在蠕虫状链模型中强制执行局部切向矢量约束的不可行性问题。
• 计算构象性质，如链端距离分布和张力下的弹性响应。
• 通过与DNA拉伸的模拟和实验数据对比，检验该方法的准确性。
• 识别在对称性破缺情形（如强向列场）下的局限性。

提出的方法

• 将局部约束 ||u(s)|| = 1 替换为对切向矢量幅值的全局平均约束。
• 使用变分自由能泛函来近似系统的实际统计力学行为。
• 应用驻定相近似推导有效持久长度和关联函数。
• 通过作用量中的二次相互作用项引入外场（如向列场或张力场）。
• 在各种极限（弱场与强场）下解析求解所得方程，以提取物理可观测量。
• 通过与已知解析极限及端到端距离和弹性性能的模拟数据对比，验证结果。

实验结果

研究问题

• RQ1 如何在保留持久长度和链端距离分布等关键物理性质的前提下，简化半柔性链的统计力学模型？
• RQ2 在张力作用下，刚性链的弹性响应如何？该平均场方法在多大程度上能再现已知的解析极限？
• RQ3 向列场如何影响半柔性链的构象性质，特别是持久长度？
• RQ4 为何平均场变分方法无法捕捉强向列场下持久长度的指数增长？
• RQ5 该方法能否仅通过轮廓长度和持久长度作为参数，定量描述实验测得的DNA拉伸数据？

主要发现

• 平均场变分方法导出的端到端距离分布表达式极为简洁，且与模拟结果几乎定量一致。
• 该方法在解析上恢复了张力作用下刚性链弹性响应的所有已知极限。
• 在DNA拉伸实验中，该模型仅使用轮廓长度和持久长度作为可调参数，即可与实验数据达到定量一致。
• 在弱向列场极限下，方法预测持久长度沿场方向线性增加，因子为 (1 + 1/3 gl₀)，与先前工作一致。
• 在强向列场极限下，该方法错误预测持久长度仅有限增加（因子为2），而精确结果表明其随 exp(√(glₚ)) 指数增长，原因在于未能捕捉瞬子效应。
• 该方法在对称性破缺区域失效，因其无法描述有效势能中简并极小值之间的隧穿行为，凸显了其在这些情形下的局限性。