典型固液界面热力学与动力学性质的分子动力学研究.docVIP

典型固液界面热力学与动力学性质的分子动力学研究 【摘要】：研究固液界面的性质对认识液体的润湿、晶体的生长与形核等都具有重要的意义。目前,采用计算模拟的方法对固液界面的性质进行研究已成为可能,而且大量有关固液界面性质的研究成果都来自于计算模拟。本文采用分子动力学方法,对典型体系固液界面的热力学与动力学性质进行分析与研究,主要研究内容可以概括为以下四个部分：一、对固体表面的纳米液滴进行分子动力学模拟,主要研究相互作用的细节对固液界面结构与润湿性质的影响。我们发现当相互作用较强,固液原子尺寸失配不是很大的情况下,固液界面附近的液体原子排列出现了长程有序结构。多体相互作用、固液相互作用的强度与距离,都会影响液滴在固体表面的润湿性质,并且随着固液相互作用强度的增大,液体原子与固体原子尺寸的失配对接触角的影响效果增强。二、对铁液滴在石墨烯表面和单壁碳纳米管中的润湿性质进行了模拟。模拟结果表明,对放在石墨烯表面的纳米铁液滴,接触角随液滴尺寸的增加而减小；但对包裹在碳管中的纳米铁液滴,接触角随管径的增加而减小；但接触角的大小与包裹液滴的长度无关。我们还根据修正后的杨氏方程估算了纳米铁液滴在石墨表面三相接触线的线张力的大小,得到了与理论研究相同的结果。通过对碳纳米管中的液体密度的分析,观察到液体的密度在垂直于管壁和平行于管壁的两个方向上都存在周期的波动,而且在碳管的轴向,密度波动的周期与碳管的周期一致。三、研究降温过程中,碳纳米管包裹的金属液滴的固化过程与固化结构。发现不同体系的降温过程中,能量随温度的变化规律不同,液体固化后的结构也不同。在碳纳米管管壁对内部液体影响较弱的体系中,随温度的降低,能量的变化规律与自由液体相似,且液体固化后形成晶体结构；反之,如果管壁对内部液体影响较强,能量随温度的变化关系复杂,降温过程中甚至观察到亚稳结构的出现,液体最终的固化结构通常为多层结构。四、采用自由固化方法,系统的研究了体系尺寸对固化过程中界面处温度梯度的影响。研究表明,体系在生长方向的长度对界面处温度梯度的影响很强。考虑到界面处温度梯度的存在,采用界面处的真实温度计算了金属铁与金属镁的动力学系数,并研究了铁和镁动力学系数的各向异性性质。采用界面处真实温度计算的结果约为采用热浴温度进行计算的两倍,但无论采用什么样的温度计算动力学系数,对相同的体系而言,其各向异性性质没有发生改变。【关键词】：分子动力学接触角碳纳米管线张力动力学系数 【学位授予单位】：华东师范大学 【学位级别】：博士 【学位授予年份】：2011 【分类号】：TB383.1 【目录】：摘要6-8ABSTRACT8-13第一章绪论13-251.1固液界面的结构13-141.2固液界面的热力学性质14-171.2.1固液界面自由能15-161.2.2三相平衡和表面润湿16-171.3固液界面的动力学性质17-181.4研究现状18-221.5本论文的主要研究内容22-25第二章分子动力学方法25-452.1分子动力学基本概念25-282.1.1物理基础25-262.1.2分子动力学模拟的步骤26-282.2原子间相互作用势28-332.2.1原子间相互作用28-302.2.2原子作用势举例30-332.2.2.1Lennard-Jones势30-312.2.2.2EAM势312.2.2.3LJ-EAM势31-322.2.2.4键序势32-332.3模拟技巧33-372.3.1时间积分算法33-362.3.1.1Verlet算法33-352.3.1.2预估校正法35-362.3.2周期性边界条件36-372.4系综的模拟与宏观物理量的控制37-452.4.1系综的模拟37-392.4.2常见控温方法39-412.4.2.1速度标度方法392.4.2.2Berendsen恒温槽方法39-402.4.2.3Andersen恒温方法402.4.2.4Nose-Hoover热浴40-412.4.3常见控压方法41-452.4.3.1Berdendsen方法41-422.4.3.2Andersen方法422.4.3.3Parrinell-Rahman方法42-432.4.3.4有限体系加压方法43-45第三章纳米液滴接触角的模拟45-613.1引言45-483.2模拟计算的方法48-523.2.1势函数的选取48-503.2.2分子模拟的方法与细节50-523.3计算结果与讨论52-593.3.1对润湿性质的影响52-583.3.2界面附近液体的结构58-593.4小结59-61第四章铁液滴对石墨表面和碳管的润湿61-794.1引言61-644.2计算细节64-674.2.1势函数的选取64-654.2.2模拟体系的选取65-664.2.3分子模拟的实现66-674.3计算结果及讨论67-774.3.1液滴