[论文解读] Hartree-Fock Model of a Self-Avoiding Flexible Polymer

C. Bouchiat|arXiv (Cornell University)|Jan 24, 2002

Electrostatics and Colloid Interactions参考文献 1被引用 3

一句话总结

本文提出一种哈特ree-福克（HF）场论方法，用于模拟具有短程排斥单体-单体相互作用的连续自由连接链模型下自避柔性聚合物的力-拉伸行为。关键成果是求解了一个自洽积分方程的解析解，该解能准确预测拉伸曲线，当排斥强度经调节后，与N=100个单体的蒙特卡罗模拟结果具有良好一致性，表明若d ≤ b，则力-拉伸形状对势能范围不敏感。

ABSTRACT

Recent measurements of the force versus extension curves in stretched single stranded DNA, under conditions where the hydrogen bonding between complementary bases is inhibited, provide a new handle for the study of self-avoiding effects in a flexible polymer. We report in this paper upon analytic computations of the force versus extension curves within a continuous version of the freely joining chain model, with monomer-monomer repulsive interactions. The problem is formulated as a Statistical Field Theory model, endowed with a fixed cutoff associated with the curvature of the extension versus force curve, in the free polymer limit. Using a Field Theory version of the Hartree-Fock approximation, the self-avoiding single polymer problem reduces to the solving of a one-dimension self-consistent integral equation. Taking the short range potential limit, we obtain an explicit analytical solution, involving the roots of an algebraic equation. In the low force regime, we find that the slope of the relative extension at zero force increases steadily with the total monomer number N, while it stays constant for a free polymer. We compare, for N=100, our analytic results with those obtained by Monte Carlo simulations of a discrete freely joining chain, in presence of excluded volume effects. Despite a difference between the repulsive potential ranges, an appropriate choice of our single free parameter leads to a good agreement between the two computations.

研究动机与目标

• 开发一种场论框架，用于模拟柔性聚合物中的自避效应，尤其适用于氢键被抑制的单链DNA条件。
• 通过构建具有单体-单体排斥作用的连续自由连接链模型，解决柔性聚合物链中的排斥体积相互作用挑战。
• 应用哈特ree-福克近似，将多体问题简化为可解的一维自洽积分方程。
• 将解析HF结果与蒙特卡罗模拟进行比较，以验证该模型对真实聚合物行为的预测能力。
• 探究力-拉伸曲线形状对排斥势能范围变化的鲁棒性。

提出的方法

• 构建自由连接链模型的连续版本，将聚合物配分函数映射为虚时下的量子力学问题。
• 通过零范围极限引入短程排斥势，使场论表述中出现狄拉克δ函数相互作用。
• 应用哈特ree-福克近似解耦多体相互作用，将问题简化为一维自洽积分方程。
• 在短程势极限下解析求解所得积分方程，得到涉及代数方程根的显式解。
• 利用该解计算作为力函数的相对拉伸度，并将曲线外推至低力区之外，以与模拟结果比较。
• 通过将解析曲线与N=100时的蒙特卡罗模拟数据匹配，校准单个自由参数（耦合常数λ），其中d=b。

实验结果

研究问题

• RQ1 与理想自由连接链相比，自避效应的引入如何改变柔性聚合物的力-拉伸响应？
• RQ2 在d ≤ b的条件下，力-拉伸曲线的形状在多大程度上对单体-单体排斥势能范围d不敏感？
• RQ3 哈特ree-福克近似能否提供对柔性聚合物中排斥体积效应的准确且可解析处理的描述？
• RQ4 解析HF解在多大程度上能准确再现具有排斥体积约束的同一系统中蒙特卡罗模拟的结果？
• RQ5 该模型能否扩展以包含更现实的单体-单体势能形式，如杨-古拉或指数形式，以更好地描述生物聚合物（如ssDNA）？

主要发现

• 在自避情况下，零力时的相对拉伸度随单体总数N的增加而增大，而理想聚合物情况下则保持不变。
• 对于N=100，当λ=0.7时，解析哈特ree-福克解与使用硬核排斥势（d=b）的蒙特卡罗模拟结果表现出良好的定量一致性。
• 当d ≤ b时，力-拉伸曲线形状对势能范围d的变化具有鲁棒性，表明对短程排斥的具体形式不敏感。
• 该模型预测在零力时仍存在非零拉伸，表明系统即使在低力下也保持伸展状态，与实验观测一致。
• 尽管解析解与模拟采用的势能极限不同（零范围 vs. 有限范围），但两者结果高度一致，表明物理行为主要由排斥体积相互作用的强度决定，而非其范围。
• 这些结果鼓励进一步探索哈特ree-福克方法在生物聚合物中的应用，包括引入非平面校正和更现实的相互作用势能。