论文部分内容阅读
本文对液态金属快速凝固技术和理论的发展与现状、纳米团簇的结构与能量以及分子动力学（MD）模拟技术的发展历程等方面进行了简要回顾。对MD模拟方法从原理到实现的各个环节、从基本的NVE系综到能够与实验结果进行直接对比分析的NPT系综进行了分析与整理,修正了有关文献中NPH和NPT系综计算方程中的错误,并对随机数、初始位形、元胞链表等MD算法进行了改进。发展完善了以HA-Pair为基础的基本团簇分析法（Basic Cluster Analysis,BCA）:将CTIM（Cluster Type Index Method）扩展为CTIM-2,以表征更多种类的基本团簇;提出用基本团簇类型来标识其中心原子的分类方法（Basic Cluster Label-ing,BCL）,BCL结合三维图形技术可简单清晰地描述相结构和纳米团簇的对称特征;并开发了配套的三维可视化图形分析软件。采用quantum Sutton-Chen（QSC）多体势对液态金属银在不同冷速下的凝固过程进行了分子动力学模拟。结果表明,冷速对凝固的最终结构具有关键性作用,且形成非晶的临界冷速约为1013K/s。高于此冷速凝固,可得到以二十面体为主要特征的金属玻璃。低于临界冷速,能形成hcp和fcc基本团簇按不同比例共存的混合结构甚至hcp、fcc晶体。且混合共存体中fcc和hcp中心原子有随机分布和规则分层两种形式。模拟得到伊辛模型为↑↑↓↓和↓↓↓↑（hcp和fcc中心原子单层和双层交替）的4H和12R分层晶体。随着冷速下降,随机分布先于规则分层形成;结晶温度Tc、fcc结构在凝固体系中的比例、体系结构的对称度都不断升高;但hcp结构和1422键型的比例却是先增后减,增减趋势发生转换的冷速约为5.0×1012K/s。细致模拟和深入分析表明:快冷结晶过程经历了两次相变三种状态:第一次液-固相变,亚稳bcc结构而非稳定的fcc结构首先从液态中形成;而且初始结晶温度Tic随冷速下降而升高。第二次固-固相变,亚稳bcc结构转变成不同的更加稳定的结构。由此可得快凝过程遵循“step rule”规则,结合形成自由能垒和快凝的特点合理解释了“随冷速降低,Cu和Ag结晶体系中hcp结构的比例先增后减”的现象。另外,同一冷速下多次模拟得到的固态结构并不相同,以及不同冷速下得到的凝固结构也并不总是遵循“冷速越慢体系能量最低”的关系,因此运动学因素对凝固过程及最终结构有重要影响。对2.5～10nm的非幻数纳米液滴,在冷速103和102K/ns下的快凝模拟表明:绝大多数凝固过程包含一次一级相变和一次连续相变;而且随着尺寸增大,初始凝固温度Tis平均值升高,凝固得到的纳米颗粒中,原子平均能量降低且趋于块体系统的平均能量。Tis的平均值与尺寸的关系,多次凝固中Tis所表现的混沌特征,以及晶体、二十面体、十面体等类型的纳米颗粒的比例等模拟结果都与实验结果符合得很好。除了规则晶体、十面体和二十面体以外,借助BCL和三维可视化方法,还发现了一类结构新颖的纳米颗粒—表面异构体（surface-isomers）。从几何结构来看,它们是在规则的晶体、十面体或二十面体外面再增加几个特殊排列方式（...ABCB或...ABCBA等）的密排原子层构成。虽然外层结构的改变在表面引入了一些非密排的（1 0 0）面,但外形更接近球体,整体能量仍可降低,特别是那些具有良好的全局（三重或五重）对称性的颗粒,应该足够稳定以致能在实验室观察到。实际上具有Umb-Dh结构的AuPd和AuCu合金颗粒已经在实验室观察到;我们模拟得到的银Fractal-Ih颗粒与Hendy对铅的模拟计算结果也完全吻合。