Exact Solution of the Discrete Wormlike Chain Model

技术摘要：离散蠕虫状链模型的精确解

问题陈述
半柔性聚合物，如双链 DNA、F-肌动蛋白和微管，是生物功能和材料设计的基础。虽然连续蠕虫状链（WLC）模型成功描述了它们的力学性质，但它依赖于缺乏简单闭式解的路径积分，因此需要诸如微扰理论或计算昂贵的数值模拟等近似方法。相反，离散蠕虫状链（DWLC）模型将聚合物离散化为具有角度变量的刚性段，将问题简化为离散概率分布。然而，获取 DWLC 的精确解一直难以实现。以往的方法通常依赖于平均场近似（假设成对相关性可分解为 $C_{i,i+1}^{rs} \approx \rho_i^r \rho_{i+1}^s$）或微扰展开，这会引入系统误差或将有效性限制在特定区域。这项工作解决了在任意外力下推导 DWLC 模型无近似精确解的挑战。

方法论
作者通过将总长度为 $L$ 的链离散化为 $N$ 个长度为 $a = L/N$ 的段，并引入取向角 $\theta_i$，重新表述了聚合物问题。弯曲能定义为相邻段之间角度差（$\Delta \theta_i$）的二次函数。

核心方法论创新在于推导了一个精确闭合关系，将成对角度分布函数（$C_{i,i+1}^{rs}$）与单点角度密度（$\rho_i^r$）联系起来。这是通过对相关聚合物构型之间的能量差异进行系统分析实现的。通过定义四个辅助分布函数并构建隔离特定弯曲项的不同哈密顿量，作者推导出了一个恒等式，将配分函数之比与能量差异联系起来。这产生了一个精确的代数关系（公式 18），将 $C_{i,i+1}^{rs}$ 与 $\rho_i^r$ 及辅助项联系起来，而无需假设因子分解。

利用该闭合关系，作者精确积分了熵泛函，避免了平均场近似。构建了包含弯曲能、精确熵项和外势的总自由能泛函 $\Omega[\rho, C]$。热力学平衡通过最小化该自由能确定，从而得到一组关于 $\rho_i^r$ 和 $C_{i,i+1}^{rs}$ 的耦合自洽方程。这些方程使用定点迭代算法数值求解，作者报告该算法收敛迅速（通常在 1–10 次迭代内）。

主要贡献

1. 精确闭合关系：推导了成对相关性与单点密度之间的严格关系，消除了对平均场因子分解假设的需求。
2. 精确熵泛函：利用精确闭合关系积分熵泛函，提供了角度相关性的完整统计力学描述。
3. 计算框架：构建了数值可解的耦合自洽方程，提供了一种计算高效的替代分子动力学（MD）模拟的方案，同时在离散模型内保持精确性。
4. 统一框架：一个理论结构，自然地将从单链力学扩展到多链系统及相变（例如各向同性 - 向列相转变）。

结果与验证
该理论框架针对已知物理极限及现有的理论/计算结果进行了验证：

• 极限情况：该模型在持久长度趋于零（$l_p \to 0$）的极限下正确恢复了随机线团行为（按 $\sqrt{N}$ 标度），并在持久长度趋于无穷大（$l_p \to \infty$）时恢复了刚性棒极限（$\langle R \rangle/Na \to 1$）。
• 与连续理论的比较：DWLC 模型生成的力 - 延伸曲线与 Marantan & Mahadevan (2018) 提供的连续 WLC 模型的解析近似表现出极好的一致性（当 $f l_p > 0.1$ 时误差 < 4%）。
• 与分子动力学的比较：该模型针对 Mazur (2007) 对短 DNA 片段（4 到 24 个碱基对）的全原子分子动力学模拟进行了基准测试。DWLC 对平均弯曲角的预测与分子动力学结果“完美”吻合（误差在 0.1°以内）。此外，该模型在 WLC 坐标（$\ln P$ 对 $1-\cos \theta$）中重现了弯曲角分布的线性化，决定系数 $R^2 = 0.999$，而未能拟合替代的次弹性链（SEC）模型。
• 相行为：该框架应用于浓聚合物溶液，成功预测了各向同性 - 向列相转变，包括临界浓度对持久长度和长径比的依赖性。

意义
该论文声称提供了离散蠕虫状链模型的第一个精确解析解，解决了聚合物统计力学中长期存在的问题。通过去除历史近似同时保持计算可行性，该框架在原子模拟和连续弹性理论之间建立了一座统一的桥梁。作者认为，这种方法为生物和材料背景下半柔性聚合物的定量研究提供了一种多功能工具，使得高通量参数探索和逆推断问题成为可能，而这些任务在以往因计算量过大而难以实现。该工作表明，具有二次弯曲能的离散形式即使在极短的长度尺度（如 DNA 的单螺旋圈）上也能准确捕捉连续聚合物物理。