下载:
-
金属纳米颗粒的性质不同于同种材料的大块金属玻璃, 如磁性、光学、电子、催化性能和吸收性能等方面[1-8], 因而在生物工程、电子器件、化学化工等领域具有很好的应用前景[5,9-12]. 研究表明, 金属纳米颗粒的性质与其颗粒的形状、尺寸、结构和组成有关[13-15]. 由于实验很难捕捉到冷却过程中任意时刻非晶态体系的结构信息, 因此鲜有快速凝固与非晶纳米颗粒微观结构的相关性的实验研究报道. 分子动力学模拟虽是一个有效的研究工具, 但由于缺乏有效的表征无序体系微观结构的方法, 对纳米液滴的结构研究难以有效开展.
材料的微观结构决定其宏观特性, 但由于非晶没有长程序, 晶体学的方法难以有效描述纳米非晶的微观结构, 最近发展的最大标准团簇分析法能不依赖任何人为设置的参数唯一客观地确定以每个原子为中心的微观结构, 为有效描述非晶材料的微观结构奠定了基础. 纳米液滴凝固是获得纳米颗粒的主要方法, 由于纳米液滴的比表面积随液滴尺寸的增大而减小, 因此纳米液滴的尺寸必然影响纳米颗粒的结构.
近年来, 作为玻璃形成能力较好的二元体系, 非晶CuZr合金纳米颗粒的力学性质分子动力学模拟得到了主要关注[16-20], 但关于CuZr纳米液滴快速凝固过程的微观结构演变的研究未见报道. 本研究选取元素比例相同但尺寸不同的CuZr纳米液滴, 采用分子动力学方法在固定的冷却速率下模拟纳米液滴的凝固过程, 分析凝固过程中基本物理量(如势能和空间分布特征)和微观结构的演化特征, 结果发现纳米液滴的凝固过程包含液-液和液-固两个连续的相变过程, 基于拓扑密堆(topologically close-packed, TCP)局域结构的数量和结构参数随温度的变化很好地体现了凝固过程的宏观特征.
-
本工作利用Lammps开源软件, 基于2009年Mendelev等[21]提出的嵌入原子模型EAM势函数, 对Cu64Zr36纳米液滴的快速凝固过程进行了分子动力学模拟. 总共模拟了6个尺寸的纳米颗粒, 分别包含1000, 2000, 4000, 5000, 7500和10000个原子数. 首先, 将Cu原子和Zr原子按64∶36的比例随机放置在一个立方体的盒子中, 模拟计算程序根据纳米颗粒包含的原子数目自动计算盒子的大小, 原子越多盒子越大. 因此, 包含不同原子数的纳米液滴在不同的模拟盒子中模拟. 为了避免在计算过程中丢失原子, 采用周期性边界条件, 并且在盒子内增加30 Å的真空层以保证纳米颗粒的有限边界特征. 在NVT系综下在2200 K (高于熔点)等温弛豫1 ns, 让系统达到平衡态, 然后以1 K/ps的冷却速率将系统快速凝固至300 K, 每间隔1 K保存一次数据, 以便进行分析.
周期性边界条件的引入可以将立方盒子无限复制填充满整个空间, 当一个原子移出盒子边界时, 其周期性映像从相对的面进入原盒子. 纳米颗粒表面原子能量高, 在实验和模拟过程中均有可能从表面脱离成为孤立原子. 对于模拟来说, 这种情况将破坏系综(原子数目不变)的基本特征. 给纳米颗粒加上周期性边界条件后, 脱离的原子可能重新与颗粒主体结合, 否则这些脱离主体的孤立原子将逐渐远离纳米颗粒主体, 不仅破坏系综特征, 而且严重影响模拟计算的效率. 在周期性边界条件下, 在颗粒表面增加30 Å的真空层, 这个距离超过了原子相互作用势的作用范围, 因此可以保证纳米颗粒的表面原子不会与临近周期性映像得到的原子发生作用, 确保其有限边界特性. 实验制备纳米颗粒的冷速可达到100 K/ps[22], 一般而言, 对于本文所研究尺寸范围内的纳米液滴, 1 K/ps冷速可以达到热力学平衡, 虽然凝固过程会稍稍偏离平衡[23], 但这种偏离对非晶化过程的结构特征没有明显影响.
-
描述非晶合金原子结构的常用方法有Voronoi分析[24]、HA-Pair分析[25]和CTIM 团簇分析[26]等, 但是这些分析方法都需要预先设定一定的参数. 例如 Voronoi 团簇分析方法必须预先确定将多大面积的多边形近似忽略, HA-Pair和CTIM分析方法需要预先设定截断距离Rc. 为了解决结构分析结果对截断距离Rc的依赖问题, 这里使用最大标准团簇分析(largest standard cluster analysis, LaSCA)法[27-30], 以便更有效地对原子结构进行分析. LaSCA法无需预先设置Rc, 每个原子的Rc可以由文献[27]中提出的拓扑准则来确定, 可以客观有效地描述无序体系的局部原子结构.
由LaSCA方法确定的结构分为3个层次, 分别为最大标准团簇 (largest standard cluster, LaSC)、中心近邻子团簇(center neighbor sub-cluster, CNS)和共同近邻(common-near-neighbor, CNN). 一个LaSC包含中心原子及其近邻; 中心原子和任意近邻原子都能组成一个根原子对(root pair or reference pair, RP), CNS由RP及其共有近邻CNN组成. CNS的拓扑结构可以用CNS-index(Sijk)来表示, 其中i表示CNN原子总数目; j表示CNN之间的成键总数目; k为CNN原子之间的最长链的长度. 因此任意一个LaSC都可以用一组n/Sijk来描述, n是LaSC中一种Sijk的数目, 前缀S可以省略. 例如图1(a) 的LaSC可表示为[12/555], 其中包含12个S555; 图1(b) 表示的CNS为S555, 包含5个CNN(图1(c)). 体心立方(body-centered cubic, BCC)晶体的结构单元则用[6/444, 8/666] 来表征; 面心立方(face-centered cubic, FCC)和六角密排(hexagonal close-packed, HCP)晶体的结构单元分别为[12/421]和[6/421, 6/422], 具体如图1(d)—(f)所示.
图 1 LaSC的基本概念及典型结构图示 (a)二十面体及其(b)构成单元 S555; (c) S555的共有近邻之间的连接关系; (d) BCC、(e) FCC和(f) HCP晶体的基本结构单元
Figure 1. Basic concept and typical structure of LaSC: (a) Icosahedron and (b) a constituent unit S555; (c) onnection between CNNs of S555; basic structural units of three typical crystals of (d) BCC, (e) FCC, and (f) HCP.
根据欧拉定理, 只含有 S444, S555和S666的LaSC满足特定的关系: 2 × n4 + n5 = 12, 其中n4和n5分别是S444和S555的数目, 这里S666的数目n6没有拓扑限制. 若中心原子的近邻总数(配位数CN = n4 + n5 + n6)和n4已知, 则此类LaSC的拓扑结构可唯一确定, 并用这两个整数来表征. 由于n4的变化范围小(从0到6), 可用几个符号分别标记. 当n4 = 0时, 用字母Z表示; 当n4为1—6时, 依次用A—F表示. 第二个整数为配位数CN, 这样每一个LaSC便可写成“L + CN ”的形式来表征. 如二十面体LaSC [12/555]记为Z12, 而BCC晶体的结构单元[6/444, 8/666]表征为F14. 至少含有一个S555的此类结构显然不能构成简单晶体, 且只含有四面体间隙, 具有典型的TCP特征, 我们称之为TCP LaSC. 显然, BCC晶体的结构单元F14不是TCP LaSC. LaSCA可以唯一确定以每一个原子为中心的LaSC, 因此体系内的所有原子可以根据LaSC的种类进行分类. 如Z12 LaSC的中心原子为Z12原子, TCP LaSC的中心原子统称为TCP原子, 以此类推.
-
平均原子能量(atomic potential energy, APE)随温度T变化的E-T曲线可以初步反映合金在快速凝固过程中基本性质的变化. 通常, 随着温度的降低, E-T曲线的斜率发生变化表明发生了相变. 对E-T曲线的高温段和低温端进行线性拟合, 若两条拟合线在E-T曲线的同一侧相交, 即曲线的斜率单调连续变化, 是连续相变的基本特征, 则凝固一般为非晶化, 否则可能为结晶(经历一级相变). 图2为所有Cu64Zr36纳米液滴在快速凝固过程中APE随温度的变化曲线. 可以看到, E-T曲线夹在两条拟合线之间, 说明纳米液滴的凝固可能存在复杂相变或者是结晶的过程.
图 2 不同尺寸的Cu64Zr36纳米液滴快速凝固过程中平均原子能量随温度的变化
Figure 2. Evolution of average atomic potential energy of per atom with temperature during rapid solidification of Cu64Zr36 nanodroplets of different sizes.
作为统计体系结构特征的一个重要参数, 双体分布函数(pair distribution function, PDF)可反映系统原子分布的总体特征. 图3(a)是300 K时不同尺寸的Cu64Zr36纳米颗粒的PDF曲线, 可以看出PDF曲线第一峰尖锐狭窄且峰值最高, 说明体系内短程有序度较高; 曲线第二峰宽而平缓, 有轻微的分裂; 而从第三峰开始, PDF曲线不再出现有明显起伏的峰, 说明体系内的原子排列呈长程无序, 表明纳米颗粒具有非晶态的结构特征. 图3(b)为300 K时不同尺寸的Cu64Zr36纳米颗粒原子排列的三维可视化图, 可以看到所有原子分布无序且均匀, 没有晶体的长程有序周期性排列, 所呈现的分布情况具有非晶态固体的原子排列特征, 进一步表明了纳米液滴的凝固过程为玻璃转变, 最终形成了非晶态的纳米颗粒. 既然最终得到的是非晶纳米颗粒, 则E-T曲线体现的类晶化过程特征需要深入研究.
图 3 300 K时Cu64Zr36纳米颗粒的(a) PDF曲线和(b)原子排列的三维可视化图
Figure 3. (a) Pair distribution functions g(r) curves and (b) three-dimensional visualization of atomic arrangement of Cu64Zr36 nanoparticles at 300 K.
从每条E-T曲线的线性拟合线可以看出, 在高温段和低温段平均原子能量以不同的斜率线性减少, 如图2(a)—(f)所示, 黑色虚线为高温段的线性拟合线, 偏离线性变化的临界温度(Ts)意味着某种相变的开始; 随着尺寸的增加, Ts分别达到1055, 1124, 1162, 1145, 1150, 1168 K. 绿色虚线为低温段的线性拟合线, 偏离线性变化的临界温度(Tg)是最后一次相变的结束; 随着尺寸的增加, Tg分别达到644, 665, 694, 677, 670, 710 K. 从高温段的Ts到低温段的Tg凝固区间, E-T曲线的斜率不是简单的线性变化, 需要深入研究.
-
为了明确Ts到Tg凝固区间的斜率变化, 先对E-T曲线进行多次平滑处理, 再用多项式 “y = ax3 + bx2 + cx + d ”对E-T曲线进行拟合, 最后对其求一阶微分和二阶微分. 如图4(a)所示, 一阶微分表明了E-T曲线斜率变化的情况, 每条一阶微分曲线都存在极大值点, E-T曲线的斜率先增加到最大值然后减小. 图4(b)是E-T曲线多项式拟合后的二阶微分, 二阶微分为零, 即E-T曲线的拐点, 为斜率变化最快的点. 随着纳米液滴尺寸的增大, 这个临界温度分别是868, 939, 960, 939, 942, 960 K. 因此, 在Tg < T < Ts区间存在一个斜率变化最快的临界温度, 这个临界温度的存在表明纳米液滴在快速凝固过程中可能存在复杂的相变过程.
图 4 E-T曲线的(a)斜率和(b)二阶导数的变化
Figure 4. (a) Slope and (b) the second derivative of the E-T curves.
至此可以肯定, T > Ts时体系是纳米液滴, 而T < Tg时为纳米固体, 在液态和固态之间可能存在液-液相变和液-固相变. 为了确定相变的类型, 需要研究熔体的基本结构性质. LaSCA不仅不需要设定决定局域结构的截断距离, 而且可以输出体系内所有原子的平均截断半径
$ \bar R_{\rm c} $ [27], 该参数是短程序的最基本特征—短程序特征长度. 在不发生相变的条件下, 短程序的特征长度$ \bar R_{\rm c} $ 与温度呈简单的线性关系, 且液体和固体的变化率存在明显区别, 是区分液体和固体的重要结构参数.图5是各个尺寸的纳米颗粒的
$ \bar R_{\rm c} $ 随凝固温度的变化过程, 可以看到每条$ \bar R_{\rm c} $ -T曲线上不仅有一个斜率转变点, 而且这个转折点的温度恰好与能量曲线的驻点温度一致, 分别是868, 939, 960, 936, 942, 960 K, 将这个临界温度记为Tls. 因此可以确定, T > Tls系统为液态, 在温度区间Ts > T > Tls发生了液-液相变, 而Tls > T > Tg则是液-固相变. 即在[Ts, Tls ]阶段为液-液相变区间, 在[Tls, Tg]阶段为液-固相变区间.
图 5 纳米液滴凝固过程的局域结构平均截断半径
$ \bar R_{\rm c} $ 随温度的变化 (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000Figure 5. Evolution of the
$ \bar R_{\rm c} $ with temperature during the solidification of nanodroplets: (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000. -
研究表明, TCP 结构不仅能成功描述结构复杂的sigma相[31]和A15相[32]的结晶度, 还能体现金属玻璃的基本结构特征[33], 能完美解释PDF曲线第一峰和第二峰分裂出来的次峰的结构起源[34], 固态非晶体系内的TCP含量与其非晶形成能力正相关[33], 且与Tg随压强的变化高度一致[35]. 纳米液滴凝固过程中TCP原子的百分含量(这是个中程序参数, 因为只涉及最多到第二近邻的局域结构特征)随温度的变化如图6所示. 在液-液相变开始之前 (T > Ts), TCP原子的数量呈指数增长; 在液-液相变区间[Ts, Tls], TCP原子数量的增长速率变快; 在液-固相变区间[Tls, Tg], TCP原子数量持续快速增长; 在液-固相变结束之后(T < Tg), TCP原子数量线性增长. 但从图6可以看出, 除了Tg以外, 其他两个临界温度在TCP原子的数量随温度的变化过程中没有明显的趋势变化.
图 6 纳米液滴凝固过程中TCP原子的百分含量随温度的变化 (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000
Figure 6. Evolution of the percentage of TCP atoms with temperature during the solidification of nano-droplets: (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000.
考虑到Tls是连续相变的特征温度, 用基于基本结构(TCP原子)的高阶参数来深入研究结构演化过程的基本特征. 考虑由TCP原子构成的中长程序结构, 这里把互为近邻的TCP原子通过TCP键(S444, S555或者S666)相互连接构成的结构称为TCP团簇, 一个TCP团簇所包含的原子总数称为团簇的尺寸或大小. 接下来分别研究TCP和Z12结构以下两个参数的变化: 1)团簇的数量, 表示为NC; 2)最大TCP团簇尺寸, 用Smax表示. 这两个参数均为长程序参数, 因为至少涉及两个TCP原子的链接(某个原子的第三近邻), 直至整个体系的尺寸.
图7为所有纳米液滴凝固过程中TCP团簇的数量NC随温度的变化. 可以看到, 所有纳米液滴凝固过程中, 随着温度下降, NC -T为一条先增后减的单峰曲线, 3个特征温度(特别是Ts和Tls)均显而易见. 在T > Ts时, NC线性增加; 在Ts > T > Tls时, 虽然TCP原子的数目仍然在快速增加(图6), 但是NC基本不变; 在Tls > T > Tg时, NC快速减少; 在T < Tg时, NC缓慢减少.
图 7 纳米液滴凝固过程中TCP团簇的数量NC随温度的变化 (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000
Figure 7. Evolution of the number of TCP clusters NC with temperature during the solidification of nano-droplets: (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000.
图8是不同尺寸的纳米液滴凝固过程中最大TCP团簇的尺寸Smax随温度的变化. 可以看到, 在T > Ts时, Smax很小且没有明显的变化. 在这个阶段, 虽然TCP原子和TCP团簇的数量都在增加, 但是团簇与团簇之间比较分散和独立, 因此TCP团簇的尺寸都还比较小. 在液-液相变区间[Ts, Tls], Smax开始有所增加, 但是增加的速率比较缓慢. 在这个阶段TCP团簇的数量不变(见图7), 而TCP原子却不断增长, 增长的TCP原子比较平均地分配给各个小团簇, 因此Smax增加缓慢. 在液-固相变区间[Tls, Tg], Smax增加的速率变快, 此时TCP团簇的数目快速减少(见图7), TCP原子仍在增加, TCP团簇快速长大或相互合并使得Smax快速增加; 在T < Tg以后, TCP原子的增速降低, TCP团簇的减少速率变得缓慢, 于是Smax增加的速率也逐渐变得缓慢.
图 8 纳米液滴在凝固过程中最大TCP团簇的尺寸Smax随温度的变化 (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000
Figure 8. Evolution of the size of the maximum TCP cluster (Smax) with temperature during the solidification of nano-droplets: (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000.
以上结果表明, 纳米液滴快速凝固过程中, TCP结构的演化与3个临界温度密切相关. 图9 进一步显示3个临界温度与纳米颗粒内TCP原子的数目随尺寸的变化趋势基本一致: 都随着尺寸的增加呈非线性非单调的变化, 并且都在N = 4000时取得极大值后有所减少, 然后增加. TCP结构与3个临界温度(Ts, Tls和Tg)的强相关性以及固态纳米颗粒内的TCP原子数目随尺寸的变化趋势的一致性, 表明TCP结构体现了纳米液滴和纳米颗粒的本质结构特征.
图 9 临界温度和TCP原子数目与尺寸的相关性 (a)液-液相变的起始温度Ts; (b)液-固相变的起始温度Tls; (c)玻璃转变温度Tg; (d) 300 K时纳米颗粒内TCP原子的百分含量
Figure 9. Correlation of critical temperature and the percentage of TCP atoms with the size of nano-droplets: (a) Initial temperature of the liquid-liquid transformation (Ts); (b) initial temperature of the liquid-solid transformation (Tls); (c) glass transition temperature (Tg); (d) percentage of TCP atoms in nanoparticles at 300 K.
一直以来, 具有五重称性的二十面体通常是研究过冷液体和非晶体结构的重点关注对象[36-40]. 图10(a)—(c)为所有纳米液滴在快速凝固过程中基于Z12的结构参数随温度的变化. 可以看到, 随着温度的降低, Z12原子的数目、Z12团簇的数量 (NC)和最大Z12团簇的尺寸Smax均呈非线性增加, 但这3个参数随温度的变化只能较好体现Tg这一个临界温度, 其他两个临界温度附近没有明显的趋势变化. 图10(d)为300 K是纳米颗粒内Z12原子的百分含量随纳米颗粒尺寸的变化, 与图9(a)—(c)对比可知, Z12原子数量随尺寸的变化趋势与3个临界温度随尺寸的变化趋势不一致. 因此, TCP结构比二十面体能更有效地描述Cu64Zr36纳米液滴的结构特征和凝固过程发生的关键结构变化.
图 10 纳米液滴在凝固过程中基于Z12的结构参数随温度的演化 (a) Z12原子数量; (b) Z12团簇数量(NC); (c) 最大Z12团簇的尺寸Smax; (d) 300 K时纳米颗粒内Z12 原子的百分含量随颗粒尺寸的变化
Figure 10. Evolution of Z12-based structure parameters with temperature during solidification of nano-droplets: (a) Number of Z12 atoms; (b) number of Z12 clusters (NC); (c) size of the largest Z12 cluster (Smax); (d) evolution of the percentage of Z12 atoms with the size of nanoparticles at 300 K.
图6显示, 在液-液相变期间(Tls < T < Ts), TCP原子的数量增加. 而图7表明团簇数量基本不变, 则TCP团簇的平均尺寸增加. 但从图8可知, 最大TCP团簇的尺寸在此温度区间没有明显变化, 所以此阶段的TCP团簇尺寸差别不大, 且缓慢长大. 这种情况类似于结晶过程的形核阶段: 晶核数目快速增加, 但均小于临界尺寸, 平均尺寸缓慢增长. 对于在此阶段的非晶体系, 独立的TCP原子聚集形成大量较小的TCP团簇. 液-固相变区间(Tg < T < Tls), TCP原子继续快速增加(图6), TCP团簇的数目快速减少(图7), TCP团簇的尺寸必然快速增大, 最大TCP团簇的尺寸Smax指数增长, 且纳米液滴的尺寸越大, 增长速率越快(见图8). 此阶段类似于结晶过程的晶核长大阶段, 液-固相变是TCP团簇的长大阶段. 当系统温度小于Tg时, TCP原子继续缓慢增加, TCP团簇的数量缓慢减少, 最大TCP团簇的尺寸缓慢增大, 是TCP团簇的粗化阶段. 而温度高于Ts时, TCP原子数量缓慢增长, 但基本相互独立, 可称为TCP团簇的坯胎阶段.
因此, 纳米液滴凝固得到非晶纳米颗粒的过程经历了高温熔体、液-液相变、液-固相变和非晶固体这4个阶段, 对应着TCP团簇的坯胎、聚集、长大和粗化这4个结构演化阶段. 但从二十面体的角度来看, 这4个阶段并不明显. 因此, 以TCP结构为研究对象, 可以更为全面地解读金属液滴凝固过程的结构演化本质, 推动凝固理论的发展.
综上所述, TCP结构是CuZr合金熔体的内禀结构, 随着温度降低逐渐增加, 并经过液-液相变和液-固相变由液体凝固成非晶态固体. 作为液态的内禀结构, 固态体系中的TCP团簇与液态体系中的同种结构相比, 只是空间结构上更加规则, 这种结构变化能通过局域位置的微调来完成, 因此其能量差异较小, 也没有明显的结构转化能垒, 因此凝固过程中宏观上的能量体积均连续缓慢变化, 体现为连续相变而非具有突变特征的一级相变. 另一方面, 如果高温熔体内晶体结构的微观结构单元也随着温度降低而(快速)增加, 则与TCP结构存在强烈的竞争, 而且晶态微观结构能量更低, 一旦结晶过程开始, 将迅速完成, 形成一级相变. 如果要得到非晶固体, 则必须利用极高的冷却速率抑制结晶将内禀的TCP结构保留至固态. 相反, 如果高温熔体内, 晶态微观结构很少, 则容易形成非晶固体. 因此可以通过调节合金熔体的成分增强其TCP结构的形成能力来增强非晶形成能力.
-
运用分子动力学方法, 对不同尺寸的Cu64Zr36纳米液滴的快速凝固过程进行模拟及分析. 结果表明, 所有纳米液滴在快速凝固过程中均发生玻璃转变, 最后形成非晶态的纳米颗粒. 通过对纳米液滴的凝固过程的结构演化分析得出以下结论:
1) LaSCA输出的短程序特征长度
$ \bar R_{\rm c} $ 随温度的变化特征能有效区分凝固过程金属体系的物态变化.2) 纳米液滴的凝固连续经历了液-液相变和液-固相变两个关键阶段, 液-固相变在Tg结束.
3) TCP结构的演化过程能充分体现纳米液滴的两次相变的基本特征, 以二十二面体为中心的结构演化不能完整体现这些特点.
4) 从TCP结构的角度, 纳米液滴的凝固过程可以分为坯胎、聚集、长大和粗化这4个阶段.
相对于以前备受关注的二十面体, 基于TCP结构的参数更能充分地描述熔体结构演化的细节, 为理解纳米液滴凝固过程的微观结构演变提供了一个新的思路和方法, 对非晶态凝固理论的发展具有重要意义.
-
由于冷却技术和合金非晶形成能力的限制, 实验室难以得到大块非晶, 而纳米液滴的快速冷却要相对容易, 因此纳米液滴的模拟研究更容易得到实验的验证. 本文运用分子动力学方法, 模拟不同尺寸的Cu64Zr36纳米液滴在1.0 × 1012 K/s冷却速率下的凝固过程, 并采用平均原子能量、双体分布函数、三维可视化和最大标准团簇分析等方法分析其微观结构的演化. 对能量曲线和微观结构短程序特征长度的统计分析表明, 所有纳米液滴的凝固过程都经历了液-液相变和液-固相变, 最后形成了非晶态纳米颗粒. 拓扑密堆(topologically close-packed, TCP)结构的演化过程能充分体现纳米液滴两次相变的基本特征, 但二十面体不能. 从TCP团簇的角度, 纳米液滴的整个凝固过程可以分为坯胎、聚集、长大和粗化4个阶段. TCP结构能体现出非晶纳米液滴和颗粒的基本结构特征, 对于完善凝固理论具有重要意义.It is difficult to obtain bulk amorphous alloys experimentally due to the limitation of cooling technology and the ability to form amorphous alloy. However, the rapid cooling of nano-droplets is relatively easy, so the simulation research of nano-droplets is easier to verify experimentally. In this work, the molecular dynamics simulation for the rapid cooling of Cu64Zr36 nano-droplets of different sizes is conducted at a cooling rate of 1.0 × 1012 K/s, and the evolution of microstructure is analyzed in terms of the average potential energy, the pair distribution function, the three-dimensional visualization, and the largest standard cluster analysis. The analysis of the energy curves and the characteristic length for short-range-ordered microstructure show that the solidification process for all nano-droplets undergoes liquid-liquid transition and liquid-solid transition, and finally forms amorphous nanoparticles. Comparing with the icosahedron, the evolution of the topologically close-packed (TCP) structures can reflect the basic characteristics of phase transitions effectively. Based on the evolution of TCP clusters, the entire solidification process of nano-droplets can be divided into four stages: embryo, aggregation, growth and coarsening. The TCP structure embodies the basic structural characteristics of amorphous nano-droplets and particles, which is of great significance in perfecting the solidification theory.
-
Keywords:
- CuZr alloy /
- topologically close-packed structure /
- LaSC /
- molecular dynamics simulation /
- solidification process
[1] Duan S B, Wang R M 2013 Prog. Nat. Sci. 23 113
Google Scholar
[2] Taylor M G, Austin N, Gounaris C E, Mpourmpakis G 2015 ACS Catal. 5 6296
Google Scholar
[3] Yan Y, Warren S C, Fuller P, Grzybowski B A 2016 Nat. Nanotechnol. 11 603
Google Scholar
[4] Mpourmpakis G, Andriotis A N, Vlachos D G 2010 Nano Lett. 10 1041
Google Scholar
[5] Yan Y C, Du J S S, Gilroy K D, Yang D R, Xia Y N, Zhang H 2017 Adv. Mater. 29 1605997
Google Scholar
[6] Cuenya B R, Behafarid F 2015 Surf. Sci. Rep. 70 135
Google Scholar
[7] Wessels J M, Nothofer H G, Ford W E, Wrochem F V, Scholz F, Vossmeyer T, Schroedter A, Weller H, Yasuda A 2004 J. Am. Chem. Soc. 126 3349
Google Scholar
[8] Kelly K L, Coronado E, Zhao L L, Schatz G C 2003 J. Phys. Chem. B 107 668
Google Scholar
[9] De M, Ghosh P S, Rotello V M 2008 Adv. Mater. 20 4225
Google Scholar
[10] Kim M, Lee C, Ko S M, Nam J M 2019 J. Solid State Chem. 270 295
Google Scholar
[11] Ferrando R, Jellinek J, Johnston R L 2008 Chem. Rev. 108 845
Google Scholar
[12] Talapin D V, Lee J S, Kovalenko M V, Shevchenko E V 2010 Chem. Rev. 110 389
Google Scholar
[13] Bratlie K M, Lee H, Komvopoulos K, Yang P and Somorjai G A 2007 Nano Lett. 7 3097
Google Scholar
[14] Cuenya B R 2010 Thin Solid Films 518 3127
Google Scholar
[15] Johnston R L 1998 Philos. Trans. R. Soc. London, Ser. A 356 211
Google Scholar
[16] Adibi S, Branicio P S, Ballarini R 2016 RSC Adv. 6 13548
Google Scholar
[17] Yuan S Y, Branicio P S 2020 Int. J. Plast. 134 102845
Google Scholar
[18] Zheng K, Branicio P S 2020 Phys. Rev. Mater. 4 076001
Google Scholar
[19] Zhang M, Li Q M, Zhang J C, Zheng G P, Wang X Y 2019 J. Alloys Compd. 801 318
Google Scholar
[20] Nandam S H, Ivanisenko Y, Schwaiger R, Śniadecki Z, Mu X, Wang D, Chellali R, Boll T, Kilmametov A, Bergfeldt T, Gleiter H, Hahn H 2017 Acta Mater. 136 181
Google Scholar
[21] Mendelev M I, Kramer M J, Ott R T, Sordelet D J, Yagodin D, Popel P 2009 Philos. Mag. 89 967
Google Scholar
[22] Zhong L, Wang J W, Sheng H W, Zhang Z, Mao S X 2014 Nature 512 177
Google Scholar
[23] Nelli D, Ferrando R 2019 Nanoscale 11 13040
Google Scholar
[24] Mauro N A, Wessels V, Bendert J C, Klein S, Gangopadhyay A K, Kramer M J, Hao S G, Rustan G E, Kreyssig A, Goldman A I, Kelton K F 2011 Phys. Rev. B 83 184109
Google Scholar
[25] Honeycutt J D, Andersen H C 1987 J. Phys. Chem. 91 4950
Google Scholar
[26] Liu R S, Li J Y, Dong K J, Zheng C X, Liu H R 2002 Mater. Sci. Eng. B 94 141
Google Scholar
[27] Tian Z A, Liu R S, Dong K J, Yu A B 2011 EPL 96 36001
Google Scholar
[28] Tian Z A, Dong K J, Yu A B 2015 Ann. Phys. 354 499
Google Scholar
[29] Tian Z A, Dong K J, Yu A B 2014 Phys. Rev. E 89 032202
Google Scholar
[30] Tian Z A, Dong K J, Yu A B 2013 AIP Conf. Proc. 1542 373
Google Scholar
[31] Wu Z Z, Mo Y F, Lang L, Yu A B, Xie Q, Liu R S, Tian Z A 2018 Phys. Chem. Chem. Phys. 20 28088
Google Scholar
[32] Mo Y F, Tian Z A, Lin L, Riu R S, Zhou L L, Hou Z Y, Peng P, Zhang T Y 2019 J. Non-Cryst. Solids 513 111
Google Scholar
[33] Lin L, Deng H Q, Tian Z A, Gao F, Hu W Y, Wen D D, Mo Y F 2019 J. Alloys Compd. 775 1184
Google Scholar
[34] 栗晶晶, 田泽安 2020 低温物理学报 42 81
Google Scholar
Li J J, Tian Z A 2020 Low. Temp. Phys. Lett. 42 81
Google Scholar
[35] Zhou L L, Mo Y F, Tian Z A, Li F Z, Xie X L, Liu R S 2021 J. Mater. Sci. 56 4220
Google Scholar
[36] Sheng H W, Luo W K, Alamgir F M, Bai J M, Ma E 2006 Nature 439 419
Google Scholar
[37] Inoue A, Kimura H 2001 J. Light Metals 1 31
Google Scholar
[38] Jiang H, Wei X, Lu W, Liang D D, Wen Z, Wang Z, Xiang H, Shen J 2019 J. Non-Cryst. Solids 521 119531
Google Scholar
[39] Luo W K, Sheng H W, Alamgir F M, Bai J M, He J H, Ma E 2004 Phys. Rev. Lett. 92 145502
Google Scholar
[40] Lee M, Kim H K, Lee J C 2010 Met. Mater. Int. 16 877
Google Scholar
-
图 1 LaSC的基本概念及典型结构图示 (a)二十面体及其(b)构成单元 S555; (c) S555的共有近邻之间的连接关系; (d) BCC、(e) FCC和(f) HCP晶体的基本结构单元
Fig. 1. Basic concept and typical structure of LaSC: (a) Icosahedron and (b) a constituent unit S555; (c) onnection between CNNs of S555; basic structural units of three typical crystals of (d) BCC, (e) FCC, and (f) HCP.
图 2 不同尺寸的Cu64Zr36纳米液滴快速凝固过程中平均原子能量随温度的变化
Fig. 2. Evolution of average atomic potential energy of per atom with temperature during rapid solidification of Cu64Zr36 nanodroplets of different sizes.
图 3 300 K时Cu64Zr36纳米颗粒的(a) PDF曲线和(b)原子排列的三维可视化图
Fig. 3. (a) Pair distribution functions g(r) curves and (b) three-dimensional visualization of atomic arrangement of Cu64Zr36 nanoparticles at 300 K.
图 4 E-T曲线的(a)斜率和(b)二阶导数的变化
Fig. 4. (a) Slope and (b) the second derivative of the E-T curves.
图 5 纳米液滴凝固过程的局域结构平均截断半径
$ \bar R_{\rm c} $ 随温度的变化 (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000Fig. 5. Evolution of the
$ \bar R_{\rm c} $ with temperature during the solidification of nanodroplets: (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000.
图 6 纳米液滴凝固过程中TCP原子的百分含量随温度的变化 (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000
Fig. 6. Evolution of the percentage of TCP atoms with temperature during the solidification of nano-droplets: (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000.
图 7 纳米液滴凝固过程中TCP团簇的数量NC随温度的变化 (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000
Fig. 7. Evolution of the number of TCP clusters NC with temperature during the solidification of nano-droplets: (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000.
图 8 纳米液滴在凝固过程中最大TCP团簇的尺寸Smax随温度的变化 (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000
Fig. 8. Evolution of the size of the maximum TCP cluster (Smax) with temperature during the solidification of nano-droplets: (a) N = 1000; (b) N = 2000; (c) N = 4000; (d) N = 5000; (e) N = 7500; (f) N = 10000.
图 9 临界温度和TCP原子数目与尺寸的相关性 (a)液-液相变的起始温度Ts; (b)液-固相变的起始温度Tls; (c)玻璃转变温度Tg; (d) 300 K时纳米颗粒内TCP原子的百分含量
Fig. 9. Correlation of critical temperature and the percentage of TCP atoms with the size of nano-droplets: (a) Initial temperature of the liquid-liquid transformation (Ts); (b) initial temperature of the liquid-solid transformation (Tls); (c) glass transition temperature (Tg); (d) percentage of TCP atoms in nanoparticles at 300 K.
图 10 纳米液滴在凝固过程中基于Z12的结构参数随温度的演化 (a) Z12原子数量; (b) Z12团簇数量(NC); (c) 最大Z12团簇的尺寸Smax; (d) 300 K时纳米颗粒内Z12 原子的百分含量随颗粒尺寸的变化
Fig. 10. Evolution of Z12-based structure parameters with temperature during solidification of nano-droplets: (a) Number of Z12 atoms; (b) number of Z12 clusters (NC); (c) size of the largest Z12 cluster (Smax); (d) evolution of the percentage of Z12 atoms with the size of nanoparticles at 300 K.
-
[1] Duan S B, Wang R M 2013 Prog. Nat. Sci. 23 113
Google Scholar
[2] Taylor M G, Austin N, Gounaris C E, Mpourmpakis G 2015 ACS Catal. 5 6296
Google Scholar
[3] Yan Y, Warren S C, Fuller P, Grzybowski B A 2016 Nat. Nanotechnol. 11 603
Google Scholar
[4] Mpourmpakis G, Andriotis A N, Vlachos D G 2010 Nano Lett. 10 1041
Google Scholar
[5] Yan Y C, Du J S S, Gilroy K D, Yang D R, Xia Y N, Zhang H 2017 Adv. Mater. 29 1605997
Google Scholar
[6] Cuenya B R, Behafarid F 2015 Surf. Sci. Rep. 70 135
Google Scholar
[7] Wessels J M, Nothofer H G, Ford W E, Wrochem F V, Scholz F, Vossmeyer T, Schroedter A, Weller H, Yasuda A 2004 J. Am. Chem. Soc. 126 3349
Google Scholar
[8] Kelly K L, Coronado E, Zhao L L, Schatz G C 2003 J. Phys. Chem. B 107 668
Google Scholar
[9] De M, Ghosh P S, Rotello V M 2008 Adv. Mater. 20 4225
Google Scholar
[10] Kim M, Lee C, Ko S M, Nam J M 2019 J. Solid State Chem. 270 295
Google Scholar
[11] Ferrando R, Jellinek J, Johnston R L 2008 Chem. Rev. 108 845
Google Scholar
[12] Talapin D V, Lee J S, Kovalenko M V, Shevchenko E V 2010 Chem. Rev. 110 389
Google Scholar
[13] Bratlie K M, Lee H, Komvopoulos K, Yang P and Somorjai G A 2007 Nano Lett. 7 3097
Google Scholar
[14] Cuenya B R 2010 Thin Solid Films 518 3127
Google Scholar
[15] Johnston R L 1998 Philos. Trans. R. Soc. London, Ser. A 356 211
Google Scholar
[16] Adibi S, Branicio P S, Ballarini R 2016 RSC Adv. 6 13548
Google Scholar
[17] Yuan S Y, Branicio P S 2020 Int. J. Plast. 134 102845
Google Scholar
[18] Zheng K, Branicio P S 2020 Phys. Rev. Mater. 4 076001
Google Scholar
[19] Zhang M, Li Q M, Zhang J C, Zheng G P, Wang X Y 2019 J. Alloys Compd. 801 318
Google Scholar
[20] Nandam S H, Ivanisenko Y, Schwaiger R, Śniadecki Z, Mu X, Wang D, Chellali R, Boll T, Kilmametov A, Bergfeldt T, Gleiter H, Hahn H 2017 Acta Mater. 136 181
Google Scholar
[21] Mendelev M I, Kramer M J, Ott R T, Sordelet D J, Yagodin D, Popel P 2009 Philos. Mag. 89 967
Google Scholar
[22] Zhong L, Wang J W, Sheng H W, Zhang Z, Mao S X 2014 Nature 512 177
Google Scholar
[23] Nelli D, Ferrando R 2019 Nanoscale 11 13040
Google Scholar
[24] Mauro N A, Wessels V, Bendert J C, Klein S, Gangopadhyay A K, Kramer M J, Hao S G, Rustan G E, Kreyssig A, Goldman A I, Kelton K F 2011 Phys. Rev. B 83 184109
Google Scholar
[25] Honeycutt J D, Andersen H C 1987 J. Phys. Chem. 91 4950
Google Scholar
[26] Liu R S, Li J Y, Dong K J, Zheng C X, Liu H R 2002 Mater. Sci. Eng. B 94 141
Google Scholar
[27] Tian Z A, Liu R S, Dong K J, Yu A B 2011 EPL 96 36001
Google Scholar
[28] Tian Z A, Dong K J, Yu A B 2015 Ann. Phys. 354 499
Google Scholar
[29] Tian Z A, Dong K J, Yu A B 2014 Phys. Rev. E 89 032202
Google Scholar
[30] Tian Z A, Dong K J, Yu A B 2013 AIP Conf. Proc. 1542 373
Google Scholar
[31] Wu Z Z, Mo Y F, Lang L, Yu A B, Xie Q, Liu R S, Tian Z A 2018 Phys. Chem. Chem. Phys. 20 28088
Google Scholar
[32] Mo Y F, Tian Z A, Lin L, Riu R S, Zhou L L, Hou Z Y, Peng P, Zhang T Y 2019 J. Non-Cryst. Solids 513 111
Google Scholar
[33] Lin L, Deng H Q, Tian Z A, Gao F, Hu W Y, Wen D D, Mo Y F 2019 J. Alloys Compd. 775 1184
Google Scholar
[34] 栗晶晶, 田泽安 2020 低温物理学报 42 81
Google Scholar
Li J J, Tian Z A 2020 Low. Temp. Phys. Lett. 42 81
Google Scholar
[35] Zhou L L, Mo Y F, Tian Z A, Li F Z, Xie X L, Liu R S 2021 J. Mater. Sci. 56 4220
Google Scholar
[36] Sheng H W, Luo W K, Alamgir F M, Bai J M, Ma E 2006 Nature 439 419
Google Scholar
[37] Inoue A, Kimura H 2001 J. Light Metals 1 31
Google Scholar
[38] Jiang H, Wei X, Lu W, Liang D D, Wen Z, Wang Z, Xiang H, Shen J 2019 J. Non-Cryst. Solids 521 119531
Google Scholar
[39] Luo W K, Sheng H W, Alamgir F M, Bai J M, He J H, Ma E 2004 Phys. Rev. Lett. 92 145502
Google Scholar
[40] Lee M, Kim H K, Lee J C 2010 Met. Mater. Int. 16 877
Google Scholar
计量
- 文章访问数: 1714
- PDF下载量: 46
- 被引次数: 0
