[論文レビュー] A simple function for calculating the interaction between a molecule and a graphene sheet

Xiongce Zhao|arXiv (Cornell University)|Sep 19, 2012

Quantum Electrodynamics and Casimir Effect参考文献 2被引用数 4

ひとこと要約

本稿では、Lennard-Jones、誘導、双極子-四極子、四極子-四極子力のすべてを、垂直方向の分離距離のみに依存する解析的に導出された統合ポテンシャル関数のセットを提示する。この方法は、原子的総和と比較してエネルギー最小値の誤差が2%未満であり、計算コストを無視できるほど低く抑えつつ、吸着等温線および密度プロファイルを高い忠実度で再現する。

ABSTRACT

We present a novel potential model for calculating the interaction between a molecule and a single graphene sheet. The dispersion/repulsion, induction, dipole-quadrupole, quadrupole-quadrupole interactions between a fluid molecule and a graphene sheet are described by integrated functions that are only dependent on the separation between the molecule and the graphene along its planar normal. The derived potential functions are in excellent agreement with the computationally demanding atom-explicit summation method. Typical errors of the integrated potential are less than 2% in the energy minimum compared with the exact atom-explicit summation. To examine the practical effectiveness of the newly developed functions, Monte Carlo simulations were performed to model the adsorption of two representative gases in graphene sheets using both the integrated and atom-explicit potentials. The integrated potential results in same adsorption isotherms and density profiles for the adsorbed phase while it only requires negligible computing time compared with that using the atom-explicit method. The newly developed potential functions provide a simple and accurate approach to calculating the physical interaction between molecules and graphene sheets.

研究の動機と目的

• 分子-グラフェン相互作用を計算的に効率的かつ正確に計算するための手法を開発すること。
• グラフェン-流体系に対するシンプルで効果的なポテンシャルモデルの欠如に応えること。
• 高価な原子的総和を置き換えることで、グラフェンにおける吸着の大型分子シミュレーションを可能にすること。
• 代表的な気体（メタンおよび水）を用いて、新しいポテンシャルを原子的シミュレーションと比較して検証すること。
• 統合ポテンシャルが層状転移や密度プロファイルといった重要な物理的性質を正確に保持していることを示すこと。

提案手法

• 無限に広がるグラフェンシート全体にわたる積分を用いて、分散力、反発力、誘導力、双極子-四極子力、四極子-四極子力の解析的表現を導出する。
• グラフェン格子の対称性と周期性を活用して、全ポテンシャルエネルギーを垂直方向の分離距離zのみの関数として表現する。
• Lennard-Jonesおよび多極子に基づく相互作用項を、炭素原子の2次元平面にわたって統合する。
• さまざまなz値において、直接的な原子的総和との比較を通じて、統合ポテンシャルのエネルギー値を検証する。
• 統合ポテンシャルと原子的ポテンシャルの両方を用いてモンテカルロシミュレーションを実施し、吸着等温線および密度プロファイルを比較する。
• 次元なし化学ポテンシャルμ*を用いて、吸着相における相転移および層状構造の挙動を分析する。

実験結果

リサーチクエスチョン

• RQ1 原子的総和と比較して、解析的に統合されたポテンシャルが分子とグラフェンシート間の相互作用エネルギーをどれほど正確に再現できるか？
• RQ2 分子シミュレーションにおいて、統合ポテンシャルの計算コストは原子的総和法と比べてどの程度低いのか？
• RQ3 メタンおよび水の単層グラフェン内における吸着等温線および密度プロファイルについて、統合ポテンシャルと原子的ポテンシャルがどれほど同じ結果を再現するか？
• RQ4 統合ポテンシャルは、第一層および第二層の開始といった吸着相における層状転移を正しく捉えているか？
• RQ5 双極子-四極子および四極子-四極子項を含む多極子相互作用をモデル化する際、統合ポテンシャルはどの程度頑健か？

主な発見

• 統合ポテンシャルは、全テスト対象の相互作用タイプにおいて、原子的総和法と比較してエネルギー最小値の誤差が2%未満である。
• 統合ポテンシャルを用いたモンテカルロシミュレーションは、原子的ポテンシャルを用いた場合とほぼ区別できない吸着等温線および密度プロファイルを生成した。
• メタンの層状転移はμ* ≈ -8.6（第一層）およびμ* ≈ -7.3（第二層）で発生し、両ポテンシャルともに同一の転移位置を示した。
• 水では、μ* ≈ -33.4で連続的な濡れ転移が発生し、密度プロファイルはz = ±4.4, ±10.6 Åおよびz = ±1.5, ±13.5 Åに一致する第一層および第二層の位置を示した。
• 統合ポテンシャルを用いたシミュレーションの計算時間は、原子的総和法に比べて1%未満であった。
• この手法は、吸着相の構造的詳細、特に分子層に対応する密度プロファイルの鋭いピークを正しく捉えており、物理的正確性を裏付けた。