[論文レビュー] Multiscale modeling of a rectifying bipolar nanopore: explicit-water versus implicit-water simulations

Zoltán Ható, Mónika Valiskó+3 more|Repository of the Academy's Library (Library of the Hungarian Academy of Sciences)|Nov 30, 2017

Nanopore and Nanochannel Transport Studies参考文献 116被引用数 13

ひとこと要約

本研究では、整流性バイポーラナノポアを対象として、明示的水分子を用いた分子動力学（MD）シミュレーションと、陽性水を用いたNernst-Planck/局所平衡モンテカルロ（NP+LEMC）シミュレーションを比較した。分子スケールの詳細に顕著な差異があるにもかかわらず、両モデルは定量的に類似した電流-電圧応答および軸方向イオン濃度プロファイルを生成し、分子的リアリズムを犠牲にしても、低次元モデルが本質的なデバイススケールの物理学を捉えられることを示している。

ABSTRACT

In a multiscale modeling approach, we present computer simulation results for a rectifying bipolar nanopore on two modeling levels. In an all-atom model, we use explicit water to simulate ion transport directly with the molecular dynamics technique. In a reduced model, we use implicit water and apply the Local Equilibrium Monte Carlo method together with the Nernst-Planck transport equation. This hybrid method makes the fast calculation of ion transport possible at the price of lost details. We show that the implicit-water model is an appropriate representation of the explicit-water model when we look at the system at the device (i.e., input vs.\ output) level. The two models produce qualitatively similar behavior of the electrical current for different voltages and model parameters. Looking at details of concentration and potential profiles, we find profound differences between the two models. These differences, however, do not influence the basic behavior of the model as a device because they do not influence the $z$-dependence of the concentration profiles which are the main determinants of current. These results then address an old paradox: how do reduced models, whose assumptions should break down in a nanoscale device, predict experimental data? Our simulations show that reduced models can still capture the overall device physics correctly, even though they get some important aspects of the molecular-scale physics quite wrong; reduced models work because they include the physics that is necessary from the point of view of device function. Therefore, reduced models can suffice for general device understanding and device design, but more detailed models might be needed for molecular level understanding.

研究の動機と目的

• 整流性バイポーラナノポアを通過するイオン輸送をシミュレートするにあたり、陽性水を用いた粗粒化モデルと、全原子的明示的水分子を用いた分子動力学（MD）モデルの比較を目的とする。
• 簡略化モデルが分子的仮定に欠陥があるにもかかわらず、実験的電流-電圧特性を再現できる理由という、長年のパラドックスを解明することを目的とする。
• 分子スケールの物理的不正確さがあるにもかかわらず、デバイススケールの予測に有効な条件を同定することを目的とする。
• ナノスケールイオン輸送デバイスのマルチスケールモデリングにおいて、計算効率と分子的詳細さのトレードオフを評価することを目的とする。

提案手法

• GROMACSを用いた全原子的分子動力学（MD）シミュレーション。明示的水分子を用い、CHARMM27力場と周期的境界条件を採用した。
• 陽性水を用いた低次元モデルと、Nernst-Planck（NP）輸送方程式を局所平衡モンテカルロ（LEMC）と結合させ、イオンフラックスおよび電気化学ポテンシャルを計算した。
• NP+LEMCモデルでは、イオンを誘電率ε = 78.5の連続体中での荷電硬球として扱い、拡散係数はバルク値から導出した。
• 電流はMDではイオン数のカウント、NP+LEMCではフラックスの統合により計算し、両モデル間で軸方向濃度および電位プロファイルを比較した。
• 境界条件の違い：MDは周期的境界条件と外部電場を用いるが、NP+LEMCは円筒型コンパートメントのバルク濃度および電位を固定した。
• 電位プロファイルは2通りの方法で検証した：テストチャージサンプリングとポアソン方程式の解法。NP+LEMCでは良好な一致が得られ、MDでも妥当な一致が確認された。

実験結果

リサーチクエスチョン

• RQ1 整流性バイポーラナノポアの電流-電圧応答を予測するにあたり、明示的水分子を用いたMDと陽性水を用いたNP+LEMCシミュレーションは、どのように比較されるか？
• RQ2 2つのモデル間のイオンおよび水の構造に関する分子スケールの差異が、デバイススケールの電気的特性にどの程度影響を及えるか？
• RQ3 点イオンや陽性溶媒といった簡略化仮定を用いる低次元モデルが、既知の物理的近似があるにもかかわらず、実験的整流挙動を再現できるのはなぜか？
• RQ4 明示的水モデルに比べ、陽性水モデルにおいてイオン濃度および電位プロファイルのどの側面が保持され、どの側面が歪められるか？
• RQ5 分子現象の詳細を捉えられないにもかかわらず、陽性水モデルはデバイス設計に信頼できるか？

主な発見

• 2つのモデルは、さまざまな電圧、表面電荷、イオン濃度条件下で、定性的に同一の電流-電圧応答を示した。
• 明示的水と陽性水モデルの間で、軸方向（z依存）濃度プロファイルが顕著に類似しており、局所的構造や溶媒和の差にもかかわらず、一致が見られた。
• 陽性水モデルは、電流の主要因たるz依存イオン濃度プロファイルを正確に捉えているため、その予測的成功が説明できる。
• 特に狭窄部領域において、局所的電位およびイオン分布プロファイルに顕著な差が認められ、低次元モデルが分子スケールの物理を正しく再現していないことが示された。
• NP+LEMCモデルの計算コストは著しく低く（数時間）、MD（数日～数週間）に比べてパラメータースキャンに効率的である。
• NP+LEMCモデルでは、2通りの電位計算手法（テストチャージとポアソン方程式の解法）の結果が良好に一致し、数値的安定性が裏付けられた。