固液接触現象の分子動力学

固液接触现象的分子动力学

• 批准号: 01J05754
• 负责人: 木村 達人
• 金额: $ 1.28万
• 依托单位: The University of Tokyo
• 依托单位国家: 日本
• 项目类别: Grant-in-Aid for JSPS Fellows
• 财政年份: 2001
• 资助国家: 日本
• 起止时间: 2001 至 2002
• 项目状态: 已结题
• 来源: https://kaken.nii.ac.jp/grant/KAKENHI-PROJECT-01J05754/
• 关键词: 分子動力学法; 水; カーボンナノチューブ; 相変化

分子動力学法を用いて単層カーボンナノチューブ内の水から氷への相変化シミュレーションを行った.中央に炭素原子3120個からなる(10,10)のSWNT(直径1.4nm,長さ19.5nm)を配置し,その内部に192個の水分子を配置した計算系を用いた.水分子には単純でありながらも,表面張力などをよく再現できるSPC/Eモデルを適用し,炭素問共有結合ポテンシャルとしてはBrennerポテンシャルを用いた.水-炭素間のポテンシャルに関しては,van der Waals力を表す酸素原子-炭素原子間のLennard-Jonesポテンシャルと,水分子の電荷と炭素原子の四重極との間のクーロンポテンシャルの和で表現した.計算開始から最初の50psは系全体に対して速度スケーリング法による温度制御を行い,300Kに系の温度を保った.その後,系を冷却していくことにより水の相変化シミュレーションを行ったが,SWNT内の水を直接温度制御することは避け,外側のSWNTの温度制御のみによって熱流速一定条件で冷却を行った.本研究で用いた系は小さく凝固熱自体も小さいため,通常の熱流速一定条件の相変化で見られる温度一定領域は観測されなかった.冷却速度を変えた2通りの計算を行い,酸素原子間距離分布の時間変化から結晶化の判定をしたところ,220K,200Kという異なる相変化温度が観測された.220Kでの相変化では八角形のチューブ状の氷結晶が得られたが,200Kでの相変化では中心に一列水分子が並ぶことによって八角形構造の一部に欠陥が生じている結晶が生成された.

Molecular dynamics method を with い て 単 layer カ ー ボ ン ナ ノ チ ュ ー ブ の water in か ら 氷 へ の phase - change シ ミ ュ レ ー シ ョ ン を line っ た. Central に carbon atoms 3120 か ら な る の (10, 10) SWNT (19.5 nm) 1.4 nm in diameter, long さ を configuration し そ の internal に を configuration 192 の water molecules し た department を calculated using い た. Water molecules に は 単 pure で あ り な が ら も, surface tension な ど を よ く reappearance で き る SPC/E モ デ ル を applicable し, carbon asked a total of combining ポ テ ン シ ャ ル と し て は Brenner ポ テ ン シ ャ ル を with い た. Between water and carbon の ポ テ ン シ ャ ル に masato し て は, van der Waals force between を table す acid element - carbon atoms の Lennard - Jones ポ テ ン シ ャ ル と, water molecules の charge と carbon atoms の quad extremely と の between の ク ー ロ ン ポ テ ン シ ャ ル の and で し た. Computing started か ら 50 ps は の is initially all に し seaborne て speed ス ケ ー リ ン グ method に よ る temperature suppression を い, 300 k に の temperature を bartender っ た. そ の after department を cooling し て い く こ と に よ り water の phase - change シ ミ ュ レ ー シ ョ ン を line っ た が, SWNT を の water in the direct temperature suppression す る こ と は け avoidance, lateral の SWNT <s:1> temperature control みによって heat flow rate certain conditions で cooling を rows った. This research で い た department は small さ く freezing autologous も small さ い た め, usually の - の phase change heat flow condition of certain で see ら れ る temperature field は 観 measuring さ れ な か っ た. Cooling speed を - え た 2 tong り の computing を い, acid element interatomic distance distribution の time variations change か ら crystallization の determine を し た と こ ろ, 220 k, 200 k と い う different な - る phase transition temperature が 観 measuring さ れ た. 220 k で の phase - change で は octagonal の チ ュ ー ブ shape の 氷 crystallization が must ら れ た が, 200 k で の phase - change で に a は center が column water molecules and ぶ こ と に よ っ て octagonal structure の a に owe 陥 が raw じ て い る crystallization が generated さ れ た.

项目成果

S.Maruyama, T.Kimura, M.-C.Lu: "Molecular Scale Aspects of Liquid Contact on a Solid Surface"Thermal Science and Engineering. vol.10, no.6. 23-29 (2002)

S.Maruyama、T.Kimura、M.-C.Lu：“固体表面上液体接触的分子尺度方面”热科学与工程。