[論文レビュー] To Pair or Not to Pair? Machine-Learned Explicitly-Correlated Electronic Structure for NaCl in Water

Niamh O’Neill, Benjamin X. Shi+3 more|arXiv (Cornell University)|Nov 2, 2023

Spectroscopy and Quantum Chemical Studies被引用数 2

ひとこと要約

本研究では、NaClを水に溶解した系をモデル化するために、RPAおよびMP2法を用いた機械学習波動関数ポテンシャルを開発した。これにより、溶媒化したイオン間の平均力のポテンシャル（PMF）を正確に捉えることができた。RPAおよびMP2ともに、イオン対状態と溶媒共有状態のエネルギーがほぼ簡退しており（0.2 kcal/mol以内）、長年の電解質シミュレーションにおける矛盾を解消した。

ABSTRACT

The extent of ion pairing in solution is an important phenomenon to rationalise transport and thermodynamic properties of electrolytes. A fundamental measure of this pairing is the potential of mean force (PMF) between solvated ions. The relative stabilities of the paired and solvent shared states in the PMF and the barrier between them are highly sensitive to the underlying potential energy surface. However direct application of accurate electronic structure methods is challenging, since long simulations are required. We develop wavefunction based machine learning potentials with the Random Phase Approximation (RPA) and second order Moller-Plesset (MP2) perturbation theory for the prototypical system of Na and Cl ions in water. We show both methods in agreement, predicting the paired and solvent shared states to have similar energies (within 0.2 kcal/mol). We also provide the same benchmarks for different DFT functionals as well as insight into the PMF based on simple analyses of the interactions in the system.

研究の動機と目的

• 計算モデルにおける水酸化ナトリウム塩（NaCl）のイオン対形成挙動についての合意が得られていない問題に対処すること。
• 電子相関効果を正確に捉えるために、イオン-溶媒相互作用に不可欠な波動関数に基づく機械学習ポテンシャルを構築すること。
• 高水準な電子構造計算手法を用いて、水中NaClの平均力のポテンシャル（PMF）をベンチマークすること。
• RPA、MP2およびさまざまなDFT関数を用いた、溶媒化構造およびイオン対形成のエネルギー的性質の予測性能を比較すること。
• イオン-水およびイオン-イオン相互作用エネルギーの分析を通じて、イオン対形成のエネルギーに信頼性のある基準を提供すること。

提案手法

• NaClを水に溶解した系に対して、RPAおよびMP2の電子構造計算に基づく機械学習ポテンシャルのトレーニング。
• 周期的境界条件を用いた95個の水分子を含むボックスモデルを用いて、溶媒化したNa⁺およびCl⁻イオンをシミュレート。
• イオン間距離座標に沿ったアンブルマッピングを用いて、平均力のポテンシャル（PMF）を計算。
• 真空中におけるイオン-イオン結合曲線を、基準となる電子構造計算手法を用いて計算し、固有の相互作用強度を評価。
• エネルギー分解を用いてイオン-水相互作用エネルギーを定量：EIW = Esys − EW − EG。
• 溶媒化構造を解明するために、水-水、イオン-水、および水-水素間の径方向分布関数（RDF）を分析。

実験結果

リサーチクエスチョン

• RQ1 RPAおよびMP2法は、水中NaClにおけるペア状態と溶媒共有状態の相対的安定性を予測する上で、どのように比較されるか？
• RQ2 高水準な波動関数法とDFT関数を用いて計算された場合の、Na⁺とCl⁻間のPMFに定量的な差異はどの程度か？
• RQ3 revPBE-D3、r2SCAN、optB88-vdWなどの異なるDFT関数は、RPAおよびMP2と比較して、イオン-水およびイオン-イオン相互作用をどの程度正確に記述できるか？
• RQ4 核量子効果は、NaCl/水系におけるPMFおよび溶媒化構造にどの程度の影響を及ぼすか？
• RQ5 イオン-水およびイオン-イオン相互作用エネルギーの分解から、イオン対形成の本質に関するどのような知見が得られるか？

主な発見

• RPAおよびMP2ポテンシャルは、イオン対状態（CIP）と溶媒共有状態（SSIP）のエネルギーをほぼ簡退しており、差はわずか0.2 kcal/molにとどまる。
• CIPの井戸の深さは、RPAで27.9 kcal/mol、MP2で27.8 kcal/molであり、浅いが安定したイオン対形成井戸であることが示された。
• Na⁺のイオン-水相互作用エネルギーは、RPAで約-173.8 kcal/mol、MP2で-172.3 kcal/molであり、Cl⁻の相互作用はわずかに安定化的で約+4.8 kcal/molであった。
• delocalization誤差のため、revPBE-D3などのDFT関数はNa⁺およびCl⁻を水に過剰に結合させ、Cl⁻-水相互作用エネルギーが+12.68 kcal/molに達する。
• RPAおよびMP2モデルにおける水の自己拡散係数は、実験値とよく一致しており、古典的力場やDeepMDポテンシャルを上回る性能を示した。
• RDF解析により、RPAおよびMP2は現実的な溶媒殻を再現しており、実験的期待値と整合する明確な第一殻の協調数が得られた。