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$ の鎖を長さ $a = L/N$ の $N$ 個のセグメントに離散化し、配向角 $\theta_i$ を導入することでポリマー問題を再定式化した。曲げエネルギーは、隣接セグメント間の角度差（$\Delta \theta_i$）に対して二次的に定義される。

手法上の核心的な革新は、対角分布関数（$C_{i,i+1}^{rs}$）と単一サイト角密度（$\rho_i^r$）を結びつける厳密な閉関係式の導出である。これは、関連するポリマー配置間のエネルギー差の体系的な分析を通じて達成される。4 つの補助分布関数を定義し、特定の曲げ項を分離する固有のハミルトニアンを構築することで、著者らは分配関数の比とエネルギー差を関連付ける恒等式を導出した。これにより、因数分解を仮定することなく、$C_{i,i+1}^{rs}$ を $\rho_i^r$ および補助項に結びつける厳密な代数関係（式 18）が得られる。

この閉関係式を用いて、著者らは平均場近似を回避し、エントロピー汎関数を厳密に積分した。曲げエネルギー、厳密なエントロピー項、および外部ポテンシャルからなる全自由エネルギー汎関数 $\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 塩基対）の全原子 MD シミュレーションに対してベンチマークされた。平均曲げ角に関する DWLC の予測は、MD 結果と「完璧な」一致（0.1 度以内）を示した。さらに、モデルは WLC 座標における曲げ角分布の線形化（$\ln P$ 対 $1-\cos \theta$）を決定係数 $R^2 = 0.999$ で再現し、一方、代替のサブエラスティックチェーン（SEC）モデルの適合には失敗した。
• 相挙動: この枠組みは濃縮ポリマー溶液に適用され、持久長およびアスペクト比に依存する臨界濃度を含む、等方性 - ネマティック相転移の予測に成功した。

意義
本論文は、ポリマー統計力学における長年の課題を解決する、離散ワームライクチェーンモデルに対する初の厳密解析解を提供すると主張している。計算の扱いやすさを維持しつつ歴史的な近似を除去することで、このフレームワークは原子論的シミュレーションと連続体力学理論の間の統合された架け橋を確立する。著者らは、このアプローチが生物学的および材料的な文脈における半柔軟性ポリマーの定量的研究のための多用途なツールを提供し、それ以外では計算的に実行不可能な高スループットなパラメータ探索および逆推論問題を可能にすると論じている。この研究は、二次曲げエネルギーを持つ離散定式化が、非常に短い長さスケール（DNA の単一ヘリックスターン）においても連続体ポリマー物理学を正確に捉えることを実証している。