液态金属银凝固过程及微观结构演变特征的模拟研究

论文摘要

本文对液态金属快速凝固技术和理论的发展与现状、纳米团簇的结构与能量以及分子动力学（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对铅的模拟计算结果也完全吻合。

论文目录

• 摘要
• Abstract
• 第1章 绪论
• 1.1 引言
• 1.2 凝固理论的形成与发展
• 1.3 快速凝固基本原理
• 1.3.1 经典形核理论
• 1.3.2 单相合金与共晶合金的快速凝固
• 1.3.3 亚稳相与非晶态的形成
• 1.4 非平衡溶质分配理论
• 1.4.1 溶质“捕获”
• 1.4.2 非平衡溶质分配基本原理
• 1.5 快速凝固的组织特征
• 1.6 纳米团簇结构与能量
• 1.6.1 纳米团簇的形成机制
• 1.6.2 纳米团簇的能量特征
• 1.6.3 纳米团簇的结构特征
• 1.7 快速凝固的计算机模拟与发展
• 1.8 本研究工作的目的及主要内容
• 第2章 分子动力学模拟方法及改进
• 2.1 引言
• 2.2 基本原理
• 2.3 边界条件
• 2.4 运动积分
• 2.5 积分步长
• 2.6 初始条件
• 2.7 计算流程
• 2.8 不同系综的分子动力学实现
• 2.8.1 NVE系综
• 2.8.2 NVT系综
• 2.8.3 NPH系综
• 2.8.4 NPT系综
• 2.9 算法改进
• 2.9.1 公式修正
• 2.9.2 高质量随机数
• 2.9.3 均匀分布的初始位形
• 2.9.4 最佳元胞分割
• 2.10 本章小结
• 第3章 微观结构分析方法
• 3.1 引言
• 3.2 双体分布函数（PDF）
• 3.3 键型指数法（HA-Pair）
• 3.4 基本团簇分析法（BCA）
• 3.5 角分布函数（ADF）
• 3.6 三维可视化分析
• 3.7 分析方法综述
• 3.8 本研究使用的势函数和模拟方法
• 3.9 本章小结
• 第4章 凝固结构与冷却速率的关系
• 4.1 引言
• 4.2 模拟条件和方法
• 4.3 结果与讨论
• 4.3.1 双体分布函数
• 4.3.2 平均能量
• 4.3.3 键型分析
• 4.3.4 基本团簇分析
• 4.3.5 角分布函数分析
• 4.3.6 可视化分析
• 4.3.7 讨论
• 4.4 本章小结
• 第5章 亚稳BCC相在快速凝固结晶过程中的形成与转变
• 5.1 引言
• 5.2 模拟条件与方法
• 5.3 结果与讨论
• 5.3.1 两次相变与三种状态
• 5.3.2 过渡相的结构特征
• 5.3.3 过渡bcc相的形成与持续
• 5.3.4 过度bcc相的亚稳特性及其解体
• 5.3.5 最终结构与冷却速率的关系
• 5.4 讨论
• 5.5 本章小结
• 第6章 自由纳米液滴的凝固结构
• 6.1 引言
• 6.2 模拟和分析方法
• 6.3 结果与讨论
• 6.3.1 共同特征
• 6.3.2 常规结构
• 6.3.3 新型结构
• 6.3.4 讨论
• 6.4 本章小结
• 结论
• 参考文献
• 致谢
• 附录