【要約】Psi4の活用事例（分子間相互作用） SAPT計算の基礎と実践 その１ [Zenn_Python] | Summary by TechDistill
> Source: Zenn_Python
Execute Primary Source
// Problem
分子設計を行う研究者は、分子間の結合エネルギーの物理的な内訳が不明であるという課題に直面している。従来の超分子法では、以下の問題が生じる。
- ・結合エネルギーの総量しか算出できない。
- ・「なぜ結合するのか」という物理的な起源が特定できない。
- ・静電的な引力と分散力のどちらが支配的か判断できない。
- ・見かけ上の引力の裏にある、強い反発成分を見落とすリスクがある。
// Approach
研究者は、相互作用エネルギーを直接摂動理論に基づいて計算するSAPTを採用することで、この問題を解決する。Psi4を用いて、以下のステップで解析を行う。
- ・相互作用を4つの物理的成分に分解する。
1.静電 (Elst): 電荷分布によるクーロン相互作用。
2.交換反発 (Exch): パウリの排他律による反発。
3.誘導 (Ind): 分極による相互作用。
4.分散 (Disp): 電子の揺らぎによる分散力。
- ・Pythonスクリプトにより、座標定義から計算、単位変換までを自動化する。
- ・jun-cc-pvdzなどの適切な基底関数系を選択し、計算精度を確保する。
// Result
ジメトキシベンゼンとメタンの計算事例により、相互作用の「質」を定量的に明らかにした。具体的には以下の成果を得ている。
- ・合計エネルギー -0.679 kcal/mol に対し、分散力が主成分であることを特定した。
- ・静電引力と、無視できない大きさの交換反発（+0.211 kcal/mol）の存在を可視化した。
- ・超分子法では見落としがちな、引力と反発のせめぎ合いを詳細に記述できた。
Senior Engineer Insight
> 計算化学のワークフローにおいて、Pythonによる自動化は極めて重要である。Psi4はAPIが整備されており、計算パイプラインへの組み込みが容易だ。ただし、SAPT計算は超分子法に比べ計算コストが高くなる。実務では、精度（SAPT2+等）と計算リソースのトレードオフを考慮した、適切な計算レベルの選定が鍵となる。