テーマ概要
潤滑油は、自動車・建機・家電など身の回りの様々な製品に用いられており、機械をスムーズに動かし、故障を防ぐ上で欠かせない存在です。
潤滑油には、性能向上のため添加剤が配合されています。中でも、ZnDTP（ジアルキルジチオリン酸亜鉛）は有名な添加剤です。古くからその作用が研究され、化学反応により複雑な物質変化を伴うことがわかっています。詳細なメカニズム解明のため、分子シミュレーションも適用されてきました。しかし、ZnDTPが複数種の元素で構成されることから、高精度な力場の作成自体が困難でした。
今回、化学反応を観察可能かつ各種元素の組合せに対応するMatlantisを用いて、ZnDTPの化学反応ダイナミクスを計算しました。
計算モデルと計算方法
Matlantis上でZnDTPの分子構造を作成し、機械のしゅう動面を模した①α-Fe2O3(0001)表面に配置しました。ZnDTPを構成する4つのアルキル基Rは、プロピル基またはブチル基としました。また、しゅう動面の自然酸化膜が剥落し、新生面が露出した状況を模擬して、②α-Fe(110)表面上にZnDTPが配置されたモデルも作成しました。これら2つのモデルに対して、Matlantisによる分子動力学（MD）シミュレーションを行いました。
計算結果と展望
MDシミュレーションの結果、①α-Fe2O3(0001)表面のZnDTPは吸着するのみで、分解反応は確認されませんでした。一方、②α-Fe(110)表面では、時間の経過に伴ってDTPやZnが解離するなどの分解反応、Pの配位状態変化などが観察されました。電荷の変化を解析した結果、①ではZnDTPと表面との間でほとんど電荷移動しないのに対し、②では顕著な電荷移動が生じました。このように表面性状がZnDTPの電荷状態に影響を及ぼすこと、それによって反応性に相違が生じることを明らかにしました。また、①の温度を673 Kとしたところ、2つのDTPに分解する反応が生じ、実験による想定[1][2]と一致する結果を得ました。
前述の通り、本計算で扱ったモデルは、ZnDTPを構成する亜鉛・硫黄・リン・酸素・炭素・水素と表面の鉄を加えると、7種類の元素を含みます。例えば従来の古典MDによる検討では、複数元素の組合せを網羅した力場の作成が必要であり、シミュレーションを始めるまでに多大な時間を費やしていました。一方、Matlantisは元素の組合せを気にする必要がなく、ユーザーによる力場作成が不要であり、「気軽に」複雑系の化学反応シミュレーションを始めることができました。
このような強みを活かし、例えばモリブデンを含んだ添加剤とZnDTPとの相乗作用や、最適な反応条件を明らかにするなど、これまで手が出しにくかった「かゆいところに手が届く」シミュレーションが期待されます。
