進行表面波下でのグラフェン上粒子の速度ロック現象

概要

熱揺らぎと振動による高速拡散がナノスケールで検出されています。本論文では、進行表面波を持つグラフェン層上の粒子運動を、分子動力学（MD）シミュレーションと理論モデルで研究しました。粒子が特定の前提条件下で波の速度で移動し続ける「速度ロック効果」が証明されました。粒子と波打つ表面間のファンデルワールス（vdW）ポテンシャルを曲率の関数として表現し、相対波フレーム座標系でのポテンシャルの「水たまり」に基づいてメカニズムを解明しました。2つの前提条件を提案します：粒子の初期位置がポテンシャルの水たまり内にあること、初期運動エネルギーが粒子を水たまりから飛び出させないこと。パラメータ解析により、速度ロック領域が波長、振幅、粒子-波間ポアポテンシャルに影響されることを示しました。波長が小さく、振幅が大きく、vdWポテンシャルが強いほど速度ロック領域は拡大します。本研究は、ポテンシャル水たまり理論に基づく粒子の新たな協調運動を明らかにし、ナノスケールでの高速拡散現象の説明となり得ます。

はじめに

最近、マイクロ/ナノスケールで一連の表面波/フォノン誘起高速輸送と拡散現象が検出されています。まず、カーボンナノチューブやグラフェンリボンに沿った熱泳現象が広範に調査されてきました。熱揺らぎがカーボンナノチューブ（CNT）表面に軸方向熱勾配を課すことで連続的な水流を可能にすることが確認されています。非平衡MDシミュレーションで、大きなグラフェン基板上の熱勾配を用いて小さなグラフェンナノフレークの運動制御の可能性が探られました。さらに、グラフェンと六方晶窒化ホウ素（h-BN）表面での熱駆動水滴輸送がMDシミュレーションで研究されています。これらの現象は特定のフォノンモードと相関すると示唆されています。例えば、SchoenらはCNT内の熱泳運動をチューブの呼吸モードに帰属させました。Panizonらはグラフェン上の屈曲進行波/フォノンが吸着物に運動量を渡し輸送を引き起こすことを指摘しました。熱泳現象と同様に、Angelosらはグラフェン上の温度誘起伝播リップルが、液体水中の水分子自己拡散より2-3桁速い水ナノ滴の高速拡散を導くことを示しました。

熱揺らぎに加え、振動もCNT内外での粒子と液滴輸送を可能にすることが確認されています。例えば、線偏光横音響波が線形運動量をナノ滴に渡すと、ナノ滴が約30 nm/nsの速度でナノチューブに沿って輸送されます。Guoらは、振動カンチレバー内の水分子が遠心力で駆動され、CNTの固定端から自由端へ連続流動することをMDシミュレーションで実証しました。単壁カーボンナノチューブ（SWCNT）を用いた新規ナノスケール単方向水分子輸送が、振動電荷と非対称表面エネルギー複合SWCNTで設計されました。Zhouらは、機械的振動で駆動されるSWCNTベースのナノサイズ水ポンプでの電流反転を調査し、水電流が機械的振動周波数に敏感に依存することを確認しました。ChangとGuoはCNTでのドミノ波を発見し、内部分子を最大1 km/sの高速で射出できることを示しました。単壁カーボンナノチューブでの可逆ドミノプロセスも証明されています。

熱揺らぎと振動による高速拡散と輸送現象がナノスケールで多数検出される中、表面の上下運動が拡散と輸送を促進することは確認されています。しかし、波と粒子運動の関係は不明確で統一されていません。主な説明は表面の運動量が粒子や表面外の液滴に輸送されるというものです。しかし、振幅、周波数、粒子-表面相互作用の関係はこの説明から導き出せません。さらに、Angelosらはグラフェン表面での吸着物高速拡散にはグラフェン曲率の一方の符号への明確な選好が必要であることを指摘し、波打つ表面形態による相互作用ポテンシャルが高速拡散と密接に関連することを示唆しました。したがって、波打つ表面と外部粒子の相互作用を探ることは、ナノスケールでの高速輸送と拡散メカニズム理解に不可欠です。

本論文では、レナード・ジョーンズ（L-J）対ポテンシャルで描写されるvdW相互作用に基づき、波打つグラフェン表面外部の粒子を研究することで、波動と粒子速度の協調関係をMDシミュレーションで実証しました。粒子が波打つ表面に落下した際の全体速度が、特定の前提条件下で進行波と同じ速度に保たれる「速度ロック効果」が確認されました。次に、粒子と波表面の相互作用ポテンシャルを曲率の関数として表現し、ポテンシャル水たまり理論を構築しました。この理論に基づき、速度ロック効果の2つの前提条件を提案し、理論で予測される軌道と速度がMDシミュレーション結果とよく一致することを示しました。また、波長と振幅、vdW相互作用パラメータの効果を分析し、グラフェン表面での液滴サーフィン現象で検出された制御則と良く一致することを示しました。波駆動速度ロック効果のメカニズムは、粒子速度と波打つ表面の新たな協調関係を明らかにします。

手法

MDシミュレーションはLAMMPSソフトウェアパッケージで実装しました。波打つ表面を原子数密度ρ = 3.85 × 10¹⁹ m⁻²のグラフェン層と仮定しました。グラフェンシートは初期状態でz = 0 Åの平坦状態で、x方向に6344 Åの長さを持ち、単位格子サイズは6000原子です。y軸方向に周期長12.2 Åで周期境界条件を適用しました。ここでは、モデルを簡略化し波打つ表面の幾何学的効果に焦点を当てるため、質量m = 0.83 × 10⁻²⁵ kgの球形粒子を考慮しました。初期状態で、粒子をz = 7 Å、x = 200 Åに配置しました。z方向に-50 m/s、x方向に約2000 m/sの初期速度を与えました。z方向初期速度の開始時間を設定することで、粒子が波打つ表面に着地する初期位置を制御しました。

グラフェン原子間相互作用には反応性経験的結合次数（REBO）ポテンシャルを採用しました。一方、粒子Pとグラフェン各炭素原子の相互作用にはレナード・ジョーンズポテンシャルを選択しました：

u(R) = ε(σ/R)¹² - ε(σ/R)⁶ (1)

ここでε = 5.92 × 10⁻²¹ J、σ = 4 × 10⁻¹⁰ mです。粒子Pと曲面間の平衡高さはh = 4.2 × 10⁻¹⁰ mとし、垂直抗力がゼロになる条件とシミュレーション結果から決定しました（詳細は補足1: 1参照）。

進行波関数は正弦波形を取ります：

y = A sin(2πx/λ - ωt) (2)

振幅A = 1 × 10⁻⁹ m、波長λ = 21.75 × 10⁻⁹ m（別記除く）。角周波数ω = 2π/10⁻¹²（周期10 psに対応）、波速v_wave = 2175 m/s = λω/2πです。進行波を励起するため、グラフェンの左端10 Å（y ∈ [-10, 0] Å）を上記振幅と周波数でz方向に揺らしました。また、x > 6010 Åの炭素原子をクランプしてグラフェンシートを安定化しました。特に、平坦グラフェンシートをシミュレートする場合、非クランプグラフェン原子もz軸に沿って初期位置に弱いばね定数0.0938 eV/Ųで拘束します（A = 0に設定）。

非固定炭素原子に初期温度5 Kを割り当てました。この温度は、グラフェンの屈曲モードと伸張モードの調和結合による熱活性化リップルを排除し、機械的励起による進行波の効果に焦点を当てるために設定されました。構造はマイクロカノニカルアンサンブル（NVE）でタイムステップ1 fsで時間発展させました。この発展を監視し、シミュレーション全体で温度がほぼ変化しないことを確認しました。

結果と考察

波打つグラフェン表面上と平坦グラフェン表面上での粒子軌道を図1に示します。時間間隔は波打つ表面の周期としています。粒子の相対位置が波の山や谷に対して変化しないことがわかり、粒子が波打つ表面にロックされ、その速度が波速と等しいことを意味します。比較として、同一初期位置での平坦表面上の粒子全体運動は波打つ表面上より明らかに遅く、平坦表面では摩擦により粒子速度が急速に減衰しますが、波打つグラフェン表面では摩擦が働いていないように見えます。異なるシミュレーション温度と波動関数パラメータでのMDシミュレーションケースを補足1: 1に示します。原子Xeと分子C₆₀の波打つ平坦表面上での運動シミュレーションも、現象の一般性を確認するために実施し、補足1: 2に示します。

ナノスケールでの速度ロック効果メカニズムを理解するため、波打つ表面Sと外部粒子Pの相互作用を考慮したモデルを構築しました（図2a, b）。波打つ表面の波長と振幅をλとA、PとSの最近接高さをh、Sの数密度をρ_sとします。MDシミュレーションでは、PとSの相互作用をvdW相互作用とし、L-Jポテンシャルで描写します：

U_L-J = ε[(σ/r)¹² - (σ/r)⁶]

すると、PとSの相互作用はL-J対ポテンシャルに基づき、平均曲率とガウス曲率の関数として記述できることが証明されています：

U_6-12 = (4πρ_sεσ¹²/5h¹⁰)[1 - hH + h²H² + (9h²/16)(H² - K)] - (2πρ_sεσ⁶/h⁴)[1 - hH + h²H² + (3h²/4)(H² - K)] (3)

ここで、点P₁は表面S上で粒子Pに最も近い点で、HとKは点P₁での平均曲率とガウス曲率です（図2a参照）。この曲率ベースポテンシャル[式(3)]はマイクロ/ナノスケールの多くの異常現象説明に用いられてきましたが、本ケースでの信頼性は、上記パラメータでMDシミュレーションの表面ポテンシャルと比較検証し、図2cに示しました。

波打つ表面が粒子Pに与える影響を分析する前に、摩擦を調査し考慮すべきです。粒子と波表面の摩擦はナノスケールで非常に複雑です。原始的な摩擦推定として、MDで平坦グラフェン層上の粒子運動をシミュレートし（補足1: 3詳細）、便宜上波打つ表面の代わりに平坦表面を用いました。この近似は補足1: 3でポテンシャル水たまりメカニズムと合わせて推定されています。上記パラメータで摩擦はf = -5.2 × 10⁻¹³ Nと推定されました。

次に、摩擦を考慮して表面Sと粒子Pの相対ポテンシャルを調べます。まず、図2bの赤色で示す相対波フレーム座標PXYを構築します。これは波速で移動し、進行波に対して静止しています。したがって、進行波はPXYで「凍結」されます。粒子がグラフェンに対して右方向に移動し続けるため、それに作用する摩擦は表面に沿って常に左方向です。その結果、相対波フレームポテンシャルは曲率ベースポテンシャルから摩擦による仕事を引いたものになります：

P = U_n + f·x (4)

式(4)に曲率ベースポテンシャルU_nと摩擦を代入し、相対波フレームポテンシャルを評価・プロットしました（図2d）。

波フレーム座標PXYが進行波と共に移動するため、ポテンシャル内の粒子Pの初期位置が粒子軌道を決定します。粒子Pの初期速度をv₀、波速をv_waveとすると、図2dに基づき2つの前提条件を提案できます：粒子Pの初期位置が赤色ゾーンλ₁のポテンシャル水たまり内にあること；初期波フレーム運動エネルギー½m(v₀ - v_wave)² ≤ ΔUを満たすこと。すると、粒子は水たまりから飛び出せず、代わりに水たまり内で捕捉され揺動します。絶対座標の観点から、粒子Pはポテンシャル水たまり内で振動しながら進行波と共に伝播し、全体速度が波速付近にロックされる、すなわち速度ロック効果が生じます。一方、粒子Pの初期位置が青色ゾーンλ₂内にあるか、相対初期運動エネルギー½m(v₀ - v_wave)² > ΔUの場合、粒子Pは単一水たまりに留まらず、波フレームポテンシャル面に沿って左方向のより低い水たまりへホップします。絶対座標では、粒子は進行波に遅れ、別の力の平衡が満たされるまで後退します。そのような平衡の一つは、粒子がグラフェン上で停止し摩擦が消失することです。興味深いことに、文献[21]でPanizonらは速度差があると進行波が粒子によって散乱され推進力を与えることを明らかにしており、粒子の最終速度がゼロより大きくなる可能性を示唆しています。

理論を定式化しより良く説明するため、ニュートンの運動法則により粒子Pの運動方程式をさらに確立します。粒子Pに作用する駆動力は、垂直抗力F_nと接線力F_tの2部分を含みます（図2b）：

F_n = ∂U_6-12/∂h; F_t = ∂U_6-12/∂H ∇H + ∂U_6-12/∂K ∇K (5)

L-Jポテンシャルでは原子間に引力と斥力の両方が存在し、外部粒子Pは垂直抗力F_nがゼロになる高さhに留まります。高さhの決定計算は補足1: 2に記載。粒子Pのx方向運動方程式は：

mẍ = F_x - f (6)

ここで、F_xは接線力F_tのx方向成分です（図2b）。式(6)を計算すると粒子軌道が得られます。正弦波表面ではガウス曲率はゼロで平均曲率はOzx面での曲線曲率に等しく、K = 0、H = κ [52]です。式(5)を(6)に代入すると、粒子Pの運動軌道を数値的に解けます。

ロックとアンロックの例を図3に示します。初期位置（図3a）が図2dのロック領域λ₁に対応する場合、理論とMDシミュレーション結果の軌道を図3bで比較しました。平坦グラフェン表面では粒子が摩擦により非常に短時間で停止する一方、波打つ表面では粒子が右方向に移動し続けます。理論軌道はMDシミュレーション結果に近似します。この傾向は図3cの粒子速度（シミュレーション時間10倍）でさらに確認されます。粒子がロックゾーンに落ち初期速度が波速と等しいため、ポテンシャル水たまり内で振動し全体速度が波速と等しくなり、推測と一致します。初期位置（図3d）が図2dのアンロック領域λ₂に落ちる粒子の場合、波打つ表面上の粒子軌道は定数に収束します（図3e）かつ速度分布でさらに確認されます。興味深いことに、速度アンロック領域に落ちても、進行波は平坦グラフェン表面と比較して粒子運動を促進します。図3fは、シミュレーション時間より長い時間で速度がゼロに減少することを示します。更なる例を補足1: 3に示します。

ポテンシャル水たまりメカニズムによれば、粒子の速度ロック効果はポテンシャル波表面によって支配されます。パラメータ効果はポテンシャル水たまり理論に基づいて議論できます。明らかにこれらには波長λ、振幅A、周波数ω、L-Jポテンシャルパラメータが含まれます。以下の解析では、異なるパラメータに関わらず摩擦は同じと仮定します。異なる波長A、振幅λ、L-Jポテンシャルパラメータεでのポテンシャル分布を図4にそれぞれ示します。図4aは、ポテンシャル水たまり深度がλの増加とともに減少し、波長が臨界値を超えると速度ロック領域が消失することを示します。また、低周波数ωは大きなλに関連するため、速度ロック領域はω増加とともに減少します。図4bは、ポテンシャル水たまり深度がA増加とともに増大し、振幅が小さすぎると速度ロック効果が消失することを示します。A/λ比が大きすぎて損傷を防ぐ必要があることに注意です。通常、波や粒子のスケール増加とともにλとAの両方が増加します。スケール効果を研究するため、λ/A比を固定しλまたはAの変化影響を調べました。図4cは、ポテンシャル水たまり深度がλまたはA増加とともに急速に減少することを示します。これは曲率ベース駆動力がスケール増加とともに急速に減少することを示し、大規模波を持つ表面では粒子の速度ロック効果が消失することを意味します。L-Jポテンシャルパラメータεについては、対相互作用ポテンシャルが強いほど速度ロック領域が広く、弱いと速度ロック効果が消失することが確認されました（図4d）。

他の2Dナノ材料では剛性とL-Jポテンシャルパラメータが異なり、異なる周波数と波速をもたらします。パラメータ解析によれば、波打つ表面の適切な波長と振幅を選択することでポテンシャル水たまりが出現します。ポテンシャル水たまりは粒子が波打つ表面と共に移動するための必須条件であるため、この速度ロック効果は短距離相互作用下で多くの2Dナノ材料層でも確立されます。

本論文では1粒子の運動を議論しましたが、依然として熱環境の枠組み内です。ポテンシャル水たまりは、粒子と表面の結合運動の必須条件です。複数粒子の場合、全てがポテンシャル水たまり領域に位置し前提条件を満たせば、それらは捕捉され波打つ表面と共に移動します。パラメータ効果によれば、波長と振幅を調整することで粒子運動を制御可能です。波長が小さく、振幅が大きく、周波数が高い表面波では速度ロック領域が拡大し、波打つ表面での高速拡散も促進されます。パラメトリック解析は、他の多くの文献で検出された高速拡散制御則とも一致します。例えば、Angelosらはグラフェン表面のリップル振幅増加とともに拡散係数が増加することを指摘しました。振幅が十分増加すると速度ロック領域が波長全体を覆い拡散を促進する可能性があり、これは図4bで説明できます。また、谷のポテンシャルが山より常に小さいことを指摘しています（図4）、これは図4に示される山領域の低ポテンシャルに対応します。Caoらは進行表面波存在下でのナノチャネル内流体流れを研究し、速度が振幅と周波数増加とともに増加することを確認しました。これもパラメトリック解析と一致します。

MDシミュレーションは非常に短時間の特性しか反映できませんが、ポテンシャル水たまりメカニズムからこの速度ロック効果の更なる潜在応用を導き出せます。例えば、振幅と周波数を調整することで、ほぼロックまたはアンロック領域を実現し、粒子を移動または停止させることが可能です。表面の垂直方向波動運動が粒子の横方向運動に変換されることは、一種のラチェット運動に類似し、ナノ電気機械システムに応用可能かもしれません。さらに、粒子-表面相互作用が運動に影響するため、異なる対ポテンシャルを持つ粒子では波打つ表面による軌道強化が異なり、相分離につながる可能性があります。

結論

結論として、粒子と進行表面波を持つグラフェン層の独特な関係、すなわち速度ロック現象を実証しました。MDシミュレーションにより、特定条件下で粒子速度が波速付近に保たれることを確認しました。理論モデルを構築しメカニズムを解明しました。ポテンシャル面の水たまりがロック効果を支配します。ロック条件はこのモデルに基づき提案されました：粒子初期位置がポテンシャル水たまり内にあり、初期運動エネルギーが粒子を水たまりから飛び出させないこと。理論予測による粒子軌道はMDシミュレーション結果とよく一致します。波長、振幅、L-Jポテンシャルパラメータの効果を議論しました。本研究は、波打つ表面での高速拡散と輸送理解の新たな視点と、相分離への潜在応用を提供します。