[论文解读] Curvature Distribution of Worm-like Chains in Two and Three Dimensions

Shay M. Rappaport, Shlomi Medalion+1 more|ArXiv.org|Jan 21, 2008

Dynamics and Control of Mechanical Systems参考文献 1被引用 5

一句话总结

该论文表明，在三维空间中，蠕虫状链（WLCs）的曲率分布峰值出现在一个有限且非零的值，这与直观预期（零曲率为最可能）相矛盾。通过几何论证和蒙特卡洛模拟，论文证明了由于统计系综中曲率相关的测度，三维空间中恰好观测到零曲率的概率为零，而这一现象在二维空间中并不成立，此时零曲率仍是概率最高的状态。

ABSTRACT

Bending of worm-like polymers carries an energy penalty which results in the appearance of a persistence length l such that the polymer is straight on length scales smaller than l and bends only on length scales larger than this length. Intuitively, this leads us to expect that the most probable value of the local curvature of a worm-like polymer undergoing thermal fluctuations in a solvent, is zero. We use simple geometric arguments and Monte Carlo simulations to show that while this expectation is indeed true for polymers on surfaces (in two dimensions), in three dimensions the probability of observing zero curvature anywhere along the worm-like chain, vanishes.

研究动机与目标

• 研究作为蠕虫状链建模的半柔性聚合物的局部曲率统计分布。
• 解决直观预期（零曲率最可能）与三维空间中实际统计行为之间的明显矛盾。
• 阐明维度在塑造热涨落聚合物曲率分布中的作用。
• 为解释DNA及其他生物聚合物的实验曲率测量提供理论和模拟基础。
• 将模型扩展以包含本征曲率，从而实现对双链DNA序列依赖性效应的分析。

提出的方法

• 研究采用长度为Δs的离散蠕虫状链模型，通过角度增量Δθ和Δφ更新切向量。
• 利用微分几何推导曲率构型的统计测度，表明在三维空间中测度按κ dκ缩放，而在二维空间中按dκ缩放。
• 曲率通过相邻切向量之间的夹角定义：在小角度极限下，Δsκ ≈ √(2(1−cos(Δθ)))。
• 执行蒙特卡洛模拟，以WLC哈密顿量采样聚合物构象，能量项通过E = ½b∫κ² ds对曲率进行惩罚。
• 从模拟数据中计算曲率分布，并与基于福克-普朗克方程和路径积分形式的解析预测进行比较。
• 通过将弯曲能修改为∝(κ−κ₀)²，将模型推广以包含本征曲率，解析推导出相应的曲率分布，并通过模拟验证。

实验结果

研究问题

• RQ1 为何蠕虫状链在二维空间中曲率最可能值为零，而在三维空间中却不是？
• RQ2 空间维度如何影响热涨落聚合物曲率构型的统计测度？
• RQ3 三维蠕虫状链的曲率分布函数形式为何？为何其峰值出现在非零处？
• RQ4 在表面进行的DNA曲率实验测量如何与聚合物真实的三维统计行为相关联？
• RQ5 能否利用推广的WLC模型，从曲率分布测量中提取双链DNA中的本征曲率？

主要发现

• 在三维空间中，由于统计系综中的κ dκ测度，沿蠕虫状链恰好观测到零曲率的概率恒为零。
• 三维空间中的曲率分布为非高斯分布，且在有限非零曲率值处出现峰值，峰值位置按⟨κ⟩ = 1/√(lₚΔs)缩放。
• 在二维空间中，曲率分布对称且在零曲率处达到峰值，与直观预期一致。
• 蒙特卡洛模拟与基于路径积分形式推导出的解析曲率分布完全一致。
• 引入本征曲率后，分布变为p(κ|κ₀) ∝ κ exp[−½l(κ−κ₀)²]，归一化因子依赖于κ₀。
• 对具有分布本征曲率的随机共聚物的模拟结果，证实了观测曲率分布的预测卷积形式。