[論文レビュー] Bayesian Estimation of the Hydroxyl Radical Diffusion Coefficient at Low Temperature and High Pressure from Atomistic Molecular Dynamics

Carter T. Butts, Rachel W. Martin|arXiv (Cornell University)|Aug 25, 2021

Spectroscopy and Quantum Chemical Studies参考文献 60被引用数 7

ひとこと要約

本研究では、TIP4P/2005水モデルを用いた原子論的分子動力学（MD）シミュレーションを用いて、低温・高圧条件下における水溶液中でのヒドロキシルラジカル（·OH）の拡散係数を推定する階層ベイズ推論フレームワークを提示する。著者らは、298 Kおよび1 atmにおける実験的拡散係数0.23 Ų/psを再現するようにLennard-Jonesパラメータを最適化することで、CHARMM互換の·OH力場を開発し、有限サイズ効果をYehとHummerの解析的補正により補正したlog-log多項式モデルを構築した。

ABSTRACT

The hydroxyl radical is the primary reactive oxygen species produced by the radiolysis of water, and is a significant source of radiation damage to living organisms. Mobility of the hydroxyl radical at low temperatures and/or high pressures is hence a potentially important factor in determining the challenges facing psychrophilic and/or barophilic organisms in high-radiation environments (e.g., ice-interface or undersea environments in which radiative heating is a potential heat and energy source). Here, we estimate the diffusion coefficient for the hydroxyl radical in aqueous solution, using a hierarchical Bayesian model based on atomistic molecular dynamics trajectories in TIP4P/2005 water over a range of temperatures and pressures.

研究の動機と目的

• 低温・高圧条件下における水溶液中でのヒドロキシルラジカル（·OH）の拡散係数を推定すること。この条件下では実験的測定が困難である。
• TIP4P/2005水中で298 Kおよび1 atmにおける既知の拡散係数0.23 Ų/psを正確に再現できるCHARMM互換の·OH力場を構築すること。
• 分子動力学シミュレーションにおける有限サイズ効果をYehとHummerの解析的補正を用いて補正すること。
• ベイズ後確率推論に基づき、温度および圧力の関数としてD·OHを予測するlog-log多項式モデルを構築すること。
• 極限環境に生息する好冷性・好圧性生物における放射線誘発損傷のメカニズムを解明するため、極限条件下での·OH移動度の役割を、その拡散行動の定量的評価を通じて評価すること。

提案手法

• 複数のシミュレーション条件における不確実性を統合するために、階層ベイズモデルを用いて原子論的MD軌道から·OH拡散係数を推定する。
• OおよびH原子のLennard-Jones井戸の深さ（ε）および半径（rmin/2）パラメータを最適化することで、CHARMMベースの·OH力場をパラメータ化し、298 Kおよび1 atmにおける実験的D·OH = 0.23 Ų/psを再現する。
• 4次元パラメータ空間（εO、εH、rminO/2、rminH/2）を探索するため、250セットの準ランダムHalton系列サンプリングを用いる。
• 各シミュレーションでは、1 nsの生産段階をNpT系で実行し、温度制御にはランジュバンダイナミクス、圧力制御にはランジュバン-ノーゼ-フーバーピストンを用いる。
• ·OH軌道からDMDを計算するためにBullerjahnらの共分散推定法を適用し、その後YehとHummerの式を用いて無限大サイズへの補正D∞ = DMD + (ζkBT)/(2πηL) を行う。
• パラメータ選定は水素結合のエネルギー的性質に基づき、O–H供与体とO–O受容体の幾何学的配置が互いに適切に不安定化される構成を優先し、物理的現実性を確保する。

実験結果

リサーチクエスチョン

• RQ1 実験的データが欠落する低温・高圧条件下における水中·OHの拡散係数は何か？
• RQ2 環境に適した物理的正確性を持つCHARMM互換の·OH力場を、常温常圧下での既知の拡散係数を再現できるようにどのようにパラメータ化できるか？
• RQ3 MDシミュレーションにおける有限サイズ効果が推定された拡散係数に与えるバイアスはどの程度で、どのように補正できるか？
• RQ4 ·OH拡散係数は温度および圧力にどのように依存するか？また、シミュレーションデータから予測モデルを導出できるか？
• RQ5 ·OHの移動度は、氷上または深海環境に生息する極限生物における放射線誘発損傷にどのような役割を果たすか？

主な発見

• 最適な·OH力場パラメータは、εO = -0.30 kcal/mol、rminO/2 = 1.5 Å、εH = -0.15 kcal/mol、rminH/2 = 1.0 Åであり、298 Kおよび1 atmで実験的D·OH = 0.23 Ų/psを再現する。
• YehとHummerの補正により、推定された拡散係数は平均で約10%低減され、特に小さなシミュレーションボックスで最大の補正が行われる。
• 最終的なD·OHの後確率分布は、log10(D·OH) = a + b·log10(T) + c·log10(P) + d·log10(T)·log10(P) のlog-log多項式モデルで最もよく記述され、係数はベイズ推論から導出された。
• モデルは、温度低下および圧力上昇に伴いD·OHが顕著に減少することを予測しており、これは溶媒の粘度増加および分子移動度低下と整合的である。
• D·OH推定における絶対相対誤差（ARE）を最小化するパラメータセットは、1.3–1.8 ÅのO–O間隔における反発的相互作用を強く支持しており、空間的排除効果が正確な拡散を実現するために不可欠であることを示している。
• 水素結合解析により、選択されたパラメータセットがO–O対向配置をO–H供与体配置に対して適切に不安定化していることが確認され、現実的な溶媒化挙動が確保されている。