搜索

文章查询

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

非晶Ag晶化过程中不同类型晶核结构的识别与跟踪

李媛 彭平

非晶Ag晶化过程中不同类型晶核结构的识别与跟踪

李媛, 彭平
PDF
HTML
导出引用
导出核心图
  • 采用分子动力学模拟研究了非晶Ag的等温晶化过程, 通过原子轨迹逆向追踪法分析了不同类型晶体团簇的结构遗传与组态演化. 在团簇类型指数法的基础上, 根据基本团簇种类与联结方式不同, 提出了一种可区分fcc单晶、多晶与混晶团簇的分析方法. 在非晶Ag等温晶化过程中, 基于团簇结构的连续遗传性特征, 发展了一种可区分fcc单晶、多晶与混晶晶胚与晶核的结构分析技术. 结果发现: 不论临界尺寸还是几何构型, 不同类型的晶核结构都存在差异, 其中fcc单晶临界尺寸最小, 多晶次之, 混晶最大; fcc单晶与多晶壳层原子中有少量hcp和bcc原子, 而混晶壳层则全部为非晶类原子, 并且fcc单晶、多晶与混晶的临界晶核都不是球型结构.
      通信作者: 彭平, ppeng@hnu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 51871096)资助的课题.
    [1]

    Jiang Y H, Zou J T, Xiao P, Liang S H, Liu F, Gu B 2018 J. Non-Cryst. Solids 481 608

    [2]

    CelikA F, Kazanc S 2013 Physica B 63

    [3]

    Hoang V V, Long N T, Son D N 2014 Comp. Mater. Sci. 95 491

    [4]

    Zhou K, Wang H P, Chang J, Wei B 2015 Chem. Phys. Lett. 639 105

    [5]

    Gasser U, Weeks E R, Schofield A, Pusey P N, Weitz D A 2001 Science 292 258

    [6]

    Lee B S, Burr G W, Shelby R M, Raoux S, Rettner C T, Bogle S N, Darmawikarta K, Bishop S G, Abelson J R 2009 Science 326 980

    [7]

    Lechner W, Dellago C 2008 J. Chem. Phys. 129 114707

    [8]

    Brostow W, Chybicki M, Laskowski R, Rybicki J 1998 Phys. Rev. B 57 13448

    [9]

    文大东, 彭平, 蒋元祺, 田泽安, 刘让苏 2013 物理学报 62 196101

    Wen D D, Peng P, Jiang Y Q, Tian Z A, Liu R S 2013 Acta Phys. Sin. 62 196101

    [10]

    Kelton K F 1991 Solid. State. Phys. 45 75

    [11]

    Liu J, Zhao J Z, Hu Z Q 2007 Mat. Sci. Eng. A 452 103

    [12]

    Honeycutt J D, Andersen H C 1987 J. Phys. Chem. 91 4950

    [13]

    Wu Z W, Li M Z, Wang W H, Song W J, Liu K X 2013 J. Chem. Phys. 138 620

    [14]

    Li R, Wu Y, Xiao J 2014 J. Chem. Phys. 140 046102

    [15]

    Zhou L L, Liu R S, Tian Z A, Liu H R, Hou Z Y, Peng P 2016 Sci. Rep. 6 31653

    [16]

    Yu D Q, Chen M, Yang H, Lü Y J 2009 Phil. Mag. Lett. 89 44

    [17]

    Zhang S L, Wang L M, Zhang X Y, Qi L, Zhang S H, Ma M Z, Liu R P 2015 Sci. Rep. 5 8590

    [18]

    Bai Y W, Bian X F, Qin J Y, Hu L N, Yang J F, Zhang K, Zhao X L, Zhang S, Huang Y Y, Yang C C 2014 J. Appl. Phys. 115 063106

    [19]

    Kalikmanov V I 2013 Nucleation Theory (Vol. 860)(Dordrecht: Springe) pp17-41

    [20]

    Plimpton S 1995 Fast Parallel Algorithms for Short-range Molecular Dynamics (Academic Press Professional, Inc.) pp1−19

    [21]

    Yang Z Y, Lu Z X, Zhao Y P 2009 Comp. Mater. Sci. 46 142

    [22]

    Lei Y W, Sun X R, Zhou R L, Bo Z 2018 Comp. Mater. Sci. 150 1

    [23]

    Foiles S M, Baskes M I, Daw M S 1986 Phys. Rev. B 33 7983

    [24]

    Martyna G J, Tobias D J, Klein M L 1994 J. Chem. Phys. 101 4177

    [25]

    Verlet L 1967 Health Phys. 22 79

    [26]

    邓永和, 文大东, 彭超, 韦彦丁, 赵瑞, 彭平 2016 物理学报 65 066401

    Deng Y H, Wen D D, Peng C, Wei Y D, Zhao R, Peng P 2016 Acta Phys. Sin. 65 066401

    [27]

    Zhou L L, Yang R Y, Tian Z A, Mo Y F, Liu R S 2017 J. Alloy. Compd. 690 633

    [28]

    Sun B, Sun Z W, Ouyang W Z, Xu S H 2014 J.Chem. Phys. 140 134904

    [29]

    Zhang H T, Mo Y F, Liu R S, Tian Z, Liu H R, Hou Z Y, Zhou L L, Liang Y C, Peng P 2018 Mater. Res. Express 5 036507

    [30]

    Hou Z Y, Liu L X, Liu R S, Tian Z A, Wang J G 2010 J. Appl. Phys. 107 629

    [31]

    Lad K N, Jakse N, Pasturel A 2012 J. Chem. Phys. 136 3313

    [32]

    Wen D D, Peng P, Jiang Y Q, Liu R S 2013 J. Non-Crys. Solids 378 61

    [33]

    Hou Z Y, Liu L X, Liu R S, Tian Z A, Wang J G 2010 Chem. Phys. Lett. 491 172

    [34]

    Wen D D, Peng P, Jiang Y Q, Tian Z A, Liu R S, Dong K J 2014 J. Non-Cryst. Solids 388 75

    [35]

    Streitz F H, Glosli J N, Patel M V 2006 Phys. Rev. Lett. 96 225701

  • 图 1  (a)体系平均每个原子能量E与(b)双体分布函数g(r)随弛豫时间t的变化

    Fig. 1.  (a) Systemic total energy per atom E and (b) pair distribution function g(r) as a function of relaxation time t.

    图 2  (a)典型晶体原子数目Ntotal随弛豫时间t的变化; (b)特征时间节点三类晶体原子的空间分布(红色、绿色和蓝色球分别表示fcc, hcp和bcc原子)

    Fig. 2.  (a) Variation of number Ntotal of various crystalline atoms with relaxation time t; (b) their spatial distributions at several special times, where red, green and blue balls represent fcc, hcp and bcc atoms, respectively.

    图 3  fcc基本原子团的四种连接方式示意图 (红色、灰色和橙色球分别表示中心原子、表面原子和共享原子)

    Fig. 3.  Schematic diagram of VS-, ES-, FS- and IS-linkages between two basic fcc clusters, where red, gray and orange balls represent core atoms, shell atoms and shared atoms, respectively.

    图 4  三种fcc扩展晶体团簇及其芯部原子示意图(红色、绿色和灰色球分别表示fcc原子、hcp原子和其他表面原子, 芯原子间实线表示IS联结, 虚线表示FS联结) (a)单晶团簇; (b)多晶团簇; (c)混晶团簇

    Fig. 4.  Schematic diagram of three kinds of fcc extended clusters and their core atoms. Red, green and gray balls represent fcc, hcp, and other shell atoms, respectively. The solid line and dot line between core atoms represent IS-linkages and FS-linkages between fcc basic clusters, respectively. (a) Single-crystal cluster; (b) poly-crystal cluster; (c) hydrid-crystal cluster.

    图 5  被追踪fcc团簇芯原子数目n(t)及其可遗传分数f(t)随弛豫时间t的变化

    Fig. 5.  Number of core atoms n(t) and heritable fraction f(t) of traced fcc extended clusters as a function of relaxation time t.

    图 6  fcc单晶团簇芯原子的遗传与演化示意图. 左边为瞬态遗传, 右边为连续遗传, 深红色与浅红色球分别表示在两种遗传模式中可遗传和不具有遗传性的原子

    Fig. 6.  Schematic diagram of heredity and evolution of core atoms in the traced single-crystal cluster. Dark and light red balls represent heritable atoms and other atoms without heredity in transient and continuous heredities, respectively.

    图 7  三种类型fcc临界晶核示意图 (a)单晶核; (b)多晶核; (c)混晶核(红色、绿色、蓝色和灰色球分别表示fcc, hcp, bcc和其他表面原子)

    Fig. 7.  Schematic diagram of critical crystalline nuclei and their core atoms. Red, green, blue and gray balls denote fcc, hcp, bcc and other shell atoms, respectively. (a) Single-crystal; (b) poly-crystal; (c) hydrid-crystal.

    表 1  被追踪团簇在不同弛豫时间t的芯原子总数n(t)、可遗传至下一个时刻的原子数目ns(t)、可连续遗传的原子数目nc(t)及其中的IS和FS联结数lISc(t)和lFSc(t)

    Table 1.  The number of core atoms n(t), the number ns(t) and nc(t) of heritable core atoms in transient and continuous heredities at relaxation time t as well as the number lISc(t) and lFSc(t) of IS- and FS-linkages between heritable core atoms in the continuous heredity mode.

    Single-crystalPoly-crystalHydrid-crystal
    t/psnnsnclIScnnsnclISclFScnnsnclISc
    2005355134
    190383131824131318201175656116
    1803825223839252238066423352
    1703619132036191320078271721
    1603624111539241115069351315
    1502720882820881522698
    1403021773021771703487
    1303018653018651653465
    12014931169311371131
    11014710177100291110
    10022910249100371310
    902113102313100522010
    801812102012100351610
    702711103113100411810
    602271024710032910
    5018910189100201010
    401370013700014700
    3021007200011200
    下载: 导出CSV
  • [1]

    Jiang Y H, Zou J T, Xiao P, Liang S H, Liu F, Gu B 2018 J. Non-Cryst. Solids 481 608

    [2]

    CelikA F, Kazanc S 2013 Physica B 63

    [3]

    Hoang V V, Long N T, Son D N 2014 Comp. Mater. Sci. 95 491

    [4]

    Zhou K, Wang H P, Chang J, Wei B 2015 Chem. Phys. Lett. 639 105

    [5]

    Gasser U, Weeks E R, Schofield A, Pusey P N, Weitz D A 2001 Science 292 258

    [6]

    Lee B S, Burr G W, Shelby R M, Raoux S, Rettner C T, Bogle S N, Darmawikarta K, Bishop S G, Abelson J R 2009 Science 326 980

    [7]

    Lechner W, Dellago C 2008 J. Chem. Phys. 129 114707

    [8]

    Brostow W, Chybicki M, Laskowski R, Rybicki J 1998 Phys. Rev. B 57 13448

    [9]

    文大东, 彭平, 蒋元祺, 田泽安, 刘让苏 2013 物理学报 62 196101

    Wen D D, Peng P, Jiang Y Q, Tian Z A, Liu R S 2013 Acta Phys. Sin. 62 196101

    [10]

    Kelton K F 1991 Solid. State. Phys. 45 75

    [11]

    Liu J, Zhao J Z, Hu Z Q 2007 Mat. Sci. Eng. A 452 103

    [12]

    Honeycutt J D, Andersen H C 1987 J. Phys. Chem. 91 4950

    [13]

    Wu Z W, Li M Z, Wang W H, Song W J, Liu K X 2013 J. Chem. Phys. 138 620

    [14]

    Li R, Wu Y, Xiao J 2014 J. Chem. Phys. 140 046102

    [15]

    Zhou L L, Liu R S, Tian Z A, Liu H R, Hou Z Y, Peng P 2016 Sci. Rep. 6 31653

    [16]

    Yu D Q, Chen M, Yang H, Lü Y J 2009 Phil. Mag. Lett. 89 44

    [17]

    Zhang S L, Wang L M, Zhang X Y, Qi L, Zhang S H, Ma M Z, Liu R P 2015 Sci. Rep. 5 8590

    [18]

    Bai Y W, Bian X F, Qin J Y, Hu L N, Yang J F, Zhang K, Zhao X L, Zhang S, Huang Y Y, Yang C C 2014 J. Appl. Phys. 115 063106

    [19]

    Kalikmanov V I 2013 Nucleation Theory (Vol. 860)(Dordrecht: Springe) pp17-41

    [20]

    Plimpton S 1995 Fast Parallel Algorithms for Short-range Molecular Dynamics (Academic Press Professional, Inc.) pp1−19

    [21]

    Yang Z Y, Lu Z X, Zhao Y P 2009 Comp. Mater. Sci. 46 142

    [22]

    Lei Y W, Sun X R, Zhou R L, Bo Z 2018 Comp. Mater. Sci. 150 1

    [23]

    Foiles S M, Baskes M I, Daw M S 1986 Phys. Rev. B 33 7983

    [24]

    Martyna G J, Tobias D J, Klein M L 1994 J. Chem. Phys. 101 4177

    [25]

    Verlet L 1967 Health Phys. 22 79

    [26]

    邓永和, 文大东, 彭超, 韦彦丁, 赵瑞, 彭平 2016 物理学报 65 066401

    Deng Y H, Wen D D, Peng C, Wei Y D, Zhao R, Peng P 2016 Acta Phys. Sin. 65 066401

    [27]

    Zhou L L, Yang R Y, Tian Z A, Mo Y F, Liu R S 2017 J. Alloy. Compd. 690 633

    [28]

    Sun B, Sun Z W, Ouyang W Z, Xu S H 2014 J.Chem. Phys. 140 134904

    [29]

    Zhang H T, Mo Y F, Liu R S, Tian Z, Liu H R, Hou Z Y, Zhou L L, Liang Y C, Peng P 2018 Mater. Res. Express 5 036507

    [30]

    Hou Z Y, Liu L X, Liu R S, Tian Z A, Wang J G 2010 J. Appl. Phys. 107 629

    [31]

    Lad K N, Jakse N, Pasturel A 2012 J. Chem. Phys. 136 3313

    [32]

    Wen D D, Peng P, Jiang Y Q, Liu R S 2013 J. Non-Crys. Solids 378 61

    [33]

    Hou Z Y, Liu L X, Liu R S, Tian Z A, Wang J G 2010 Chem. Phys. Lett. 491 172

    [34]

    Wen D D, Peng P, Jiang Y Q, Tian Z A, Liu R S, Dong K J 2014 J. Non-Cryst. Solids 388 75

    [35]

    Streitz F H, Glosli J N, Patel M V 2006 Phys. Rev. Lett. 96 225701

  • [1] 李春丽, 段海明, 买力坦, 开来木. Aln(n=13–32)团簇熔化行为的分子动力学模拟研究. 物理学报, 2013, 62(19): 193104. doi: 10.7498/aps.62.193104
    [2] 刘海涛, 王晓路, 袁 喆, 何春龙, 李家明. 团簇的第一原理分子动力学计算研究:价键优选法. 物理学报, 2005, 54(2): 628-635. doi: 10.7498/aps.54.628
    [3] 张宗宁, 刘美林, 李蔚, 耿长建, 张林, 赵骞. 熔融Cu55团簇在Cu(010)表面上凝固过程的分子动力学模拟. 物理学报, 2009, 58(13): 67-S71. doi: 10.7498/aps.58.67
    [4] 张林, 张彩碚, 祁阳, 徐送宁. 熔融Cu55团簇在铜块体中凝固过程的分子动力学模拟. 物理学报, 2009, 58(13): 40-S46. doi: 10.7498/aps.58.40
    [5] 坚增运, 高阿红, 常芳娥, 唐博博, 张龙, 李娜. Ni熔体凝固过程中临界晶核和亚临界晶核的分子动力学模拟. 物理学报, 2013, 62(5): 056102. doi: 10.7498/aps.62.056102
    [6] 周诗韵, 王 音, 宁西京. 一种寻找团簇异构体的准动力学方法. 物理学报, 2008, 57(1): 387-391. doi: 10.7498/aps.57.387
    [7] 李儒新, 朱频频, 徐至展, 刘晶儒, 刘建胜. 大尺寸团簇在超短超强激光场中的动力学行为. 物理学报, 2001, 50(6): 1121-1127. doi: 10.7498/aps.50.1121
    [8] 张林, 张彩碚, 祁阳. 低温下Au959团簇负载于MgO(100)表面后结构变化的分子动力学研究. 物理学报, 2009, 58(13): 53-S57. doi: 10.7498/aps.58.53
    [9] 张林, 祁阳, 张宗宁, 赵骞. 低温下Cu13团簇负载于Cu(001)表面上结构变化的分子动力学研究. 物理学报, 2009, 58(13): 47-S52. doi: 10.7498/aps.58.47
    [10] 樊沁娜, 李蔚, 张林. 熔融Cu57团簇在急冷过程中弛豫和局域结构转变的分子动力学研究. 物理学报, 2010, 59(4): 2428-2433. doi: 10.7498/aps.59.2428
    [11] 张英杰, 肖绪洋, 李永强, 颜云辉. 分子动力学模拟Cu(010)基体对负载Co-Cu双金属团簇熔化过程的影响. 物理学报, 2012, 61(9): 093602. doi: 10.7498/aps.61.093602
    [12] 司丽娜, 郭丹, 雒建斌. 氧化硅团簇切削单晶硅粗糙峰的分子动力学模拟研究. 物理学报, 2012, 61(16): 168103. doi: 10.7498/aps.61.168103
    [13] 芦鹏飞, 李勇, 杨向东, 刘锦超. 从常温常压到超临界乙醇的分子动力学模拟. 物理学报, 2010, 59(7): 4880-4887. doi: 10.7498/aps.59.4880
    [14] 谢 朝, 侯 氢, 汪 俊, 孙铁英, 龙兴贵, 罗顺忠. 金属钛中氦团簇融合的分子动力学模拟. 物理学报, 2008, 57(8): 5159-5164. doi: 10.7498/aps.57.5159
    [15] 叶子燕, 张庆瑜. 低能Pt原子团簇沉积过程的分子动力学模拟. 物理学报, 2002, 51(12): 2798-2803. doi: 10.7498/aps.51.2798
    [16] 李 鹏, 宁西京, 王 音. C36团簇自组装的分子动力学研究. 物理学报, 2005, 54(6): 2847-2852. doi: 10.7498/aps.54.2847
    [17] 刘丽霞, 侯兆阳, 刘让苏. 过冷液体钾形核初期微观动力学机理的模拟研究. 物理学报, 2012, 61(5): 056101. doi: 10.7498/aps.61.056101
    [18] 张宪刚, 宗亚平, 吴艳. 相场再结晶储能释放模型与显微组织演变的模拟研究. 物理学报, 2012, 61(8): 088104. doi: 10.7498/aps.61.088104
    [19] 张春艳, 刘显明. 氢团簇在飞秒强激光场中的动力学行为. 物理学报, 2015, 64(16): 163601. doi: 10.7498/aps.64.163601
    [20] 张华伟, 李言祥. 金属熔体中气泡形核的理论分析. 物理学报, 2007, 56(8): 4864-4871. doi: 10.7498/aps.56.4864
  • 引用本文:
    Citation:
计量
  • 文章访问数:  328
  • PDF下载量:  6
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-12-13
  • 修回日期:  2019-02-14
  • 上网日期:  2019-03-23
  • 刊出日期:  2019-04-01

非晶Ag晶化过程中不同类型晶核结构的识别与跟踪

  • 湖南大学材料科学与工程学院, 长沙 410082
  • 通信作者: 彭平, ppeng@hnu.edu.cn
    基金项目: 国家自然科学基金(批准号: 51871096)资助的课题.

摘要: 采用分子动力学模拟研究了非晶Ag的等温晶化过程, 通过原子轨迹逆向追踪法分析了不同类型晶体团簇的结构遗传与组态演化. 在团簇类型指数法的基础上, 根据基本团簇种类与联结方式不同, 提出了一种可区分fcc单晶、多晶与混晶团簇的分析方法. 在非晶Ag等温晶化过程中, 基于团簇结构的连续遗传性特征, 发展了一种可区分fcc单晶、多晶与混晶晶胚与晶核的结构分析技术. 结果发现: 不论临界尺寸还是几何构型, 不同类型的晶核结构都存在差异, 其中fcc单晶临界尺寸最小, 多晶次之, 混晶最大; fcc单晶与多晶壳层原子中有少量hcp和bcc原子, 而混晶壳层则全部为非晶类原子, 并且fcc单晶、多晶与混晶的临界晶核都不是球型结构.

English Abstract

    • 非晶晶化初期的形核及其相分离、相选择和形核竞争一直是材料与凝聚态物理研究的热点[1-3]. 长期以来, 人们一直致力于发明一种设备或建立一种方法以检测和表征形核期间晶核的形成与长大[4,5], 但受限于目前衍射、光谱与影像测试技术在高温与原位测试方面的困难, 一直没有取得大的进展[6]. 计算机模拟作为有别于理论研究和实验测试的第三种方法提供了一种可能, 但在如何识别与表征晶胚与晶核的结构特征时却遇到了问题. 目前, 关于无序体系的微结构表征, 虽然已发展了局域序参数(Q4Q6)[7]、Voronoi多面体指数 <m3,m4, m5, m6…>[8]和原子团簇类型指数$Z = \sum {n/\left( {ijkl} \right)} $[9]三种局域原子结构表征方法, 但对于晶胚与晶核这种几个纳米量级的大尺寸晶体团簇[10], 却没有给出明确的表征方法. 先前对晶体团簇的定义主要基于基本原子团间的共享原子数, 如以至少共享3个原子来定义基本晶体团簇的连通性, 赵九洲课题组[11]从Honeycutt-Andersen (H-A)键型指数[12]出发, 标度和跟踪了金属Cu体系中一个由fcc与hcp组成的300个原子左右的等温弛豫晶核; 以至少共享1个邻近原子来区分两个bcc 团簇是否连通, 李茂枝课题组[13]基于Voronoi多面体表征技术[8], 分析了Zr85Cu15合金过冷液体中bcc晶核的形状、尺寸与分布及其演化等. 而对于临界晶核的识别则主要依据模拟体系能量和晶体团簇数量以及被追踪团簇尺寸的突变来推定, 如吴永全课题组[14]基于局域序参数技术, 分析和标度了过冷液体Fe中一个内部结构为fcc/hcp混合体但表面存在许多bcc原子的fcc-Fe临界晶核(N ≈ 327); 刘让苏课题组[15]通过原子团类型指数法(CTIM), 表征了过冷液态Zn等温晶化过程中一个临界晶粒尺寸N约为90—150个原子的fcc/hcp混晶晶核等.

      考虑到先前对形核初期晶胚与晶核的识别主要基于体系热力学参数的变化[11,14-16], 加上在晶核结构表征方面, 没有检测其原子组态结构是松散的扩展团簇还是致密的中程序结构, 即严格区分单晶、多晶, 还是混晶? 为此, 本文拟从相变的一个动力学特征参数(即团簇结构遗传性[9,17,18])出发, 根据经典形核理论[19]关于晶核(可连续长大)与晶胚(有可能中途消失)的概念, 重新定义一种可严格区分晶胚与晶核的团簇分析方法, 并应用于非晶Ag等温晶化时临界晶核的识别与标度.

    • 分子动力学模拟采用LAMMPS (large-scale atomic/molecular massively parallel simulator)程序[20]. 嵌入原子势(embedded atom method, EAM)是金属模拟过程常用的原子间相互作用势[21,22], 本文采用Foiles 等[23]发展的EAM势. 首先使10000个Ag原子随机分布在一个立方盒中, 并给每个原子编号(设定ID). 采用等温等压NPT系综, 恒温恒压(P = 0)调节采用修正的Nose-Hoover法[24], 求解运动方程采用Verlet算法[25], 采用三维周期性边界条件, 时间步长设置为1 fs. 模拟过程为: 首先让体系在1600 K(Ag的熔点约为1235 K)等温弛豫100 ps使之处于平衡态, 接着以1 × 1014 K/s的冷速冷却到300 K以形成非晶, 然后升温至500 K并等温驰豫10000 ps, 期间每隔10 ps收集一次数据, 以记录体系能量和每个原子的空间坐标.

    • 图1(a)示出了500 K等温弛豫过程中体系平均每个原子能量E(t)随时间t的变化. 对应着0 p, 100, 500, 2000和3000 ps, 其双体分布函数g(r)分别如图1(b)所示. 从图1可见, 0 ps时, g(r)曲线的第二峰劈裂, 其峰形呈典型的无定形结构特征, 表明体系属于非晶态; 但在100 ps时, E(t)明显减少, 并且g(r)曲线第一峰与第二峰之间出现了一个微弱次峰, 第二峰的左侧也开始变得尖锐, 说明体系此时已开始晶化, 出现了中程序. 当t ≥ 2000 ps时, E(t)骤减, 并且g(r)曲线逐渐呈现出fcc谱线特征, 表明体系中已大量存在fcc长程有序结构.

      图  1  (a)体系平均每个原子能量E与(b)双体分布函数g(r)随弛豫时间t的变化

      Figure 1.  (a) Systemic total energy per atom E and (b) pair distribution function g(r) as a function of relaxation time t.

    • 模拟体系中各原子的团簇属性及该中心原子与其近邻原子组成的局域原子组态(即短程序)采用基于H-A键型指数[12]的扩展原子团类型指数(Z, n/(ijkl),…)[26]表征. 在这里, $Z = \sum {n/\left( {ijkl} \right)} $为与该中心原子组成基本原子团的近邻原子总数(即配位数), n/(ijkl)表示中心原子与近邻原子所形成的各种H-A键对(ijkl)数目. 如由1个中心原子与Z = 12个近邻原子形成12个1421键对的局域原子组态可表示为(12 12/1421), 此时, 该中心原子的团簇属性即为fcc原子, 其局域原子结构则为fcc基本团簇. 相应地, 模拟体系中的hcp与bcc原子, 可分别表示为(12 6/1421 6/1422)和(14 6/1441 8/1661). 至于非晶二十面体团簇(记作ico), 其CTIM指数则为(12 12/1551).

      图2(a)示出了fcc, hcp, bcc与ico这4种特征原子总数Ntotal随弛豫时间t的变化, 图2(b)则进一步给出了fcc, hcp与bcc晶体原子在t = 0, 30, 100, 500, 2000和3000 ps的空间分布. 从图2(a)可见, 在弛豫的初始阶段10—100 ps区间, 上述各类型团簇的数目变化都不明显; 但从100 ps开始, 则明显可见ico数目减少, 同时fcc, hcp及bcc数目迅速增加, 并且在t ≈ 800 ps时bcc数目达到最大, 而hcp数目则在t ≈ 2000 ps 时迅速减少, 之后趋于稳定. 图2(b)进一步显示: 在t = 0 ps时, 体系中存在少量弥散分布的fcc, hcp和bcc晶体类原子; 从t ≈ 30 ps开始, 晶体类原子逐渐增多, 并且部分聚集在一起; 其间fcc, hcp和bcc相互竞争共存[27,28], 至500 ps时, 依稀可见fcc和hcp的分层排布; 在t = 2000 ps之后, 模拟体系即成为以fcc原子为主、同时存在大量hcp原子的层状结构[29].

      图  2  (a)典型晶体原子数目Ntotal随弛豫时间t的变化; (b)特征时间节点三类晶体原子的空间分布(红色、绿色和蓝色球分别表示fcc, hcp和bcc原子)

      Figure 2.  (a) Variation of number Ntotal of various crystalline atoms with relaxation time t; (b) their spatial distributions at several special times, where red, green and blue balls represent fcc, hcp and bcc atoms, respectively.

    • CTIM只能描述以一个原子为中心的基本团簇(即短程序), 然而图2(b)显示这些基本原子团在系统中并非都是孤立存在的. 部分邻近的基本团簇间会相互联结, 形成尺寸更大的扩展团簇[9,30-32]. 以fcc基本团簇为例, 图3示出了2个基本团簇分别通过顶点共享(VS)、边共享(ES)、面共享(FS)和交叉共享(IS)形成的4个最小扩展团簇. 不难发现: 只有在IS联结模式下, fcc基本原子团的中心原子才能与共享原子形成一个H-A键(即1421键). 此时, 相对于VS, ES和FS联结, IS联结的扩展团簇共享原子最多、原子堆积致密、结构稳定性高[9]. 因此, 对于fcc晶体类团簇而言, 只有以fcc基本键型1421键联结的fcc基本原子团才能被称为fcc中程序. 也就是说, 基本团簇间是否存在H-A键, 将是区分中程序与扩展团簇的唯一判据[9]. 对于如hcp这种双键型构成的基本团簇来说, 其IS联结则存在1421与1422二种键型.

      图  3  fcc基本原子团的四种连接方式示意图 (红色、灰色和橙色球分别表示中心原子、表面原子和共享原子)

      Figure 3.  Schematic diagram of VS-, ES-, FS- and IS-linkages between two basic fcc clusters, where red, gray and orange balls represent core atoms, shell atoms and shared atoms, respectively.

      除了fcc原子外, 模拟体系中还存在hcp与bcc晶态原子, 这些晶态原子的局域原子组态结构间无疑也会存在IS联接, 只是不同基本团簇间的H-A键型不同而已. 这样, 为了区别IS联结时基本原子团种类与其间H-A键型的不同, 可定义: 以IS联接(即存在H-A键对)的同种团簇即为单晶团簇, 而异种原子团的IS联接则称为混晶团簇. 对于以fcc原子为主的晶体类团簇, fcc单晶团簇无疑就是fcc中程序, 而fcc混晶团簇则存在fcc/hcp, fcc/bcc与fcc/hcp/bcc三种类型. 进一步, 考虑到FS也是一种比较致密的联结方式, 不仅常见于非晶合金中[30-32], 而且数量往往比VS和ES高得多[32], 因此, 以FS连接的fcc短程序或中程序本文将定义为fcc多晶团簇. 这样, 针对图2(b)所示各种晶态原子(特别是fcc原子)的分布, 本文提供了一种可区分其集聚状态的结构分析方法.

      作为示例, 本文选取图2(b)t = 30 ps时虚线环绕的fcc原子集聚区, 分别给出了一个fcc单晶和fcc/hcp混晶以及fcc多晶的团簇结构示意图, 参见图4. 在这里, 红色与绿色球分别表示fcc与hcp原子, 通过实线或虚线连接, 构成了团簇的芯部; 灰色球表示非fcc, hcp和bcc等晶态原子, 共同构成了团簇的壳层结构. 其中, 图4(a)表示的是中心原子编号分别为3665和5104的两个fcc基本团簇, 通过IS联结, 形成了一个包含20个原子的fcc单晶团簇; 图4(b)表示的是图4(a)所示fcc单晶团簇与中心原子编号分别为3176和8228的另一个fcc单晶团簇, 通过FS联结(注: 3个共享原子编号分别为4737, 5668和6387), 形成了一个由37个原子组成的fcc多晶团簇; 而图4(c)表示的则是当计及不同类型晶体团簇的IS联结时, 图4(a)所示fcc单晶团簇中编号为5294的hcp原子从fcc壳层原子变成了hcp中心原子, 这样, 以3665和5104两个fcc原子为芯的单晶团簇就变成了一个以两个fcc原子(即3665和5104)和一个hcp原子(即5294)为芯的包含27个原子的fcc/hcp混晶团簇.

      图  4  三种fcc扩展晶体团簇及其芯部原子示意图(红色、绿色和灰色球分别表示fcc原子、hcp原子和其他表面原子, 芯原子间实线表示IS联结, 虚线表示FS联结) (a)单晶团簇; (b)多晶团簇; (c)混晶团簇

      Figure 4.  Schematic diagram of three kinds of fcc extended clusters and their core atoms. Red, green and gray balls represent fcc, hcp, and other shell atoms, respectively. The solid line and dot line between core atoms represent IS-linkages and FS-linkages between fcc basic clusters, respectively. (a) Single-crystal cluster; (b) poly-crystal cluster; (c) hydrid-crystal cluster.

    • 采用原子轨迹逆向追踪法[32,33], 本文进一步分析了体系弛豫过程中fcc单晶、多晶与混晶团簇的遗传与演化. 在这里, 仅仅团簇的芯部(即中心原子组合)被追踪, 并且遗传被定义为相对于前一个状态保持部分中心原子不变的结构演化. 具体即: 在t1时刻选取一个扩展团簇, 通过原子编号ID, 追踪其前一个时刻t2 < t1的结构形态, 定义与t1时刻具有最多相同ID原子的同类型扩展团簇为前驱体, 进而以f(t2)=ns(t2)/n(t1)表征被追踪团簇t2时刻能遗传到t1的遗传分数, 其中, n(t1)为t1时刻团簇的中心原子(即芯原子)总数, ns(t2)为遗传自上一个时刻t2的中心原子数目. 显然, f(t)可表征被追踪团簇在t时刻的瞬态遗传分数. 以t = 3000 ps时体系中最大的fcc单晶、多晶与混晶团簇为例, 图5示出了它们在不同弛豫时间的芯原子数目n(t)和瞬态遗传分数f(t). 由图5可见, 即使在体系能量E(t)很少变化的t < 100 ps区间, 这些团簇也具有结构遗传性, 但f(t)波动大, 说明体系正经历剧烈的原子重排. 在200 ps之后, n(t)和f(t)随弛豫时间t延长而逐渐增大, 并且相对于fcc多晶与混晶团簇, fcc单晶团簇f(t)更高, 同时明显可见, 混晶n(t)最多, 多晶次之, 单晶最少. 表明此阶段fcc单晶、多晶与混晶团簇的演化主要以生长模式为主. 至2000 ps时, 对应着n(t)的迅速升高, f(t)急剧变化, 尤其是fcc单晶团簇. 之后n(t)和f(t)趋于稳定. 表明t = 2000 ps前后, 伴随着大量其他结构属性的原子变成fcc原子, fcc单晶团簇经历了一个从生长到合并再到生长的结构演变过程.

      图  5  被追踪fcc团簇芯原子数目n(t)及其可遗传分数f(t)随弛豫时间t的变化

      Figure 5.  Number of core atoms n(t) and heritable fraction f(t) of traced fcc extended clusters as a function of relaxation time t.

    • 鉴于在t > 200 ps以后, 上述被追踪团簇的演化主要以生长为主, 在瞬态遗传的基础上, 本文进一步考察了它们在200 ps之前的连续遗传性[34]. 在这里, 连续遗传被定义为: 在弛豫历经的各态中一直能检测到某个特征晶型的团簇碎片. 具体即: 从t0 = 200 ps时刻开始, 以$\Delta t$为间隔, 连续逆向追踪t0之前每一个时刻ti = (t0i·$\Delta t$)(其中i为记录的步数)该特征团簇的前驱体, 并以nc(ti)和lc(ti)分别表征tit0每个时刻都存在的具有相同ID的中心原子数和这些具有相同ID中心原子间的IS或FS联结数目(lISclFSc), 直到lc(tn$\Delta t$) = 0. 则在tn时刻, 可连续遗传中心原子仅构成最小fcc单晶、多晶或混晶的团簇碎片. 标记该时刻为连续遗传的起始时刻tn = tonset, 其对应的前驱体即为可连续遗传团簇的初始结构. 按经典形核理论[19], 本文认为可连续遗传的扩展团簇即为晶核, 不具备连续遗传性的扩展团簇则是晶胚, 而与tonset时刻对应的晶核无疑就是临界晶核.

      $\Delta t$ = 10 ps为例, 表1列出了图5所示fcc单晶、多晶与混晶团簇在不同弛豫时间t的芯原子总数n(t)、可遗传至下一个时刻的原子数目ns(t)、可连续遗传的原子数目nc(t)及其中的IS和FS联结数lISc(t)和lFSc(t). 从表1可见, 对于同样的弛豫时间t, 可连续遗传的原子数nc(t)总是小于遗传至下一个时刻的原子数ns(t), 并且当nc(t) ≤ 1时, ns(t)仍然存在一个不小的数值, 表明瞬态遗传与连续遗传反映的是不同中心原子组态结构的演化特征.

      Single-crystalPoly-crystalHydrid-crystal
      t/psnnsnclIScnnsnclISclFScnnsnclISc
      2005355134
      190383131824131318201175656116
      1803825223839252238066423352
      1703619132036191320078271721
      1603624111539241115069351315
      1502720882820881522698
      1403021773021771703487
      1303018653018651653465
      12014931169311371131
      11014710177100291110
      10022910249100371310
      902113102313100522010
      801812102012100351610
      702711103113100411810
      602271024710032910
      5018910189100201010
      401370013700014700
      3021007200011200

      表 1  被追踪团簇在不同弛豫时间t的芯原子总数n(t)、可遗传至下一个时刻的原子数目ns(t)、可连续遗传的原子数目nc(t)及其中的IS和FS联结数lISc(t)和lFSc(t)

      Table 1.  The number of core atoms n(t), the number ns(t) and nc(t) of heritable core atoms in transient and continuous heredities at relaxation time t as well as the number lISc(t) and lFSc(t) of IS- and FS-linkages between heritable core atoms in the continuous heredity mode.

      以fcc单晶团簇为例, 图6进一步示出了它们芯部原子的演化路径与结构形态变化. 从图6可见, 虽然在90—110 ps区间一直存在不少瞬态可遗传原子, 但能连续遗传到200 ps的原子却只有1个孤立fcc原子, 即1785, 此时lISc(t)始终为0, 且团簇时大时小, 说明这个阶段fcc单晶团簇为不稳定晶胚. 在t = 120 ps时, fcc单晶团簇中出现了3个可连续遗传的fcc原子, 即8339, 1785和6873, 并且8339和1785间存在IS联结, 即lISc(120 ps) ≥1, 表明此时可连续遗传的fcc原子成为了一个建造fcc单晶团簇的中程序碎片. 当t ≥ 120 ps时, 这个中程序碎片一直存在, 且团簇尺寸持续增大, 因此, t = 120 ps可认为是连续遗传的起点tonset, 此时的fcc单晶团簇即为临界晶核.

      图  6  fcc单晶团簇芯原子的遗传与演化示意图. 左边为瞬态遗传, 右边为连续遗传, 深红色与浅红色球分别表示在两种遗传模式中可遗传和不具有遗传性的原子

      Figure 6.  Schematic diagram of heredity and evolution of core atoms in the traced single-crystal cluster. Dark and light red balls represent heritable atoms and other atoms without heredity in transient and continuous heredities, respectively.

      为清晰起见, 图7进一步给出了fcc单晶临界晶核的团簇结构与几何构型. 作为比较, 图7同时给出了fcc多晶与混晶临界晶核的芯部与团簇结构形态. 从图7可见, fcc单晶晶核临界最小, 芯部仅有n = 14原子, 加上64个壳层原子, 总共不过N = 78原子. 芯部原子呈近平面分布, 团簇结构呈近半球型, 尺寸约为1.6 nm × 1.0 nm × 1.1 nm, 并且壳层(即表面)还存在少量hcp和bcc原子. 与fcc单晶临界晶核比较, fcc多晶临界晶核稍大(n = 16, N = 95, 尺寸约为1.7 nm × 1.0 nm × 1.6 nm), 团簇结构呈近块体型, 但fcc混晶临界晶核则大得多(n = 37, N = 201, 尺寸约为2.3 nm × 1.9 nm × 2.4 nm), 团簇结构呈山字型, 并且其壳层全部为非晶体类原子. 明显地, fcc单晶、多晶与混晶的形核发生在体系能量突降之后与扩展晶体团簇尺寸陡增之前(参见图1图5), 并且, 不同类型的晶核结构不论临界尺寸还是几何构型都存在差异.

      图  7  三种类型fcc临界晶核示意图 (a)单晶核; (b)多晶核; (c)混晶核(红色、绿色、蓝色和灰色球分别表示fcc, hcp, bcc和其他表面原子)

      Figure 7.  Schematic diagram of critical crystalline nuclei and their core atoms. Red, green, blue and gray balls denote fcc, hcp, bcc and other shell atoms, respectively. (a) Single-crystal; (b) poly-crystal; (c) hydrid-crystal.

      值得指出的是, 上述结果虽然仅从一个团簇逆向追踪导出, 不能认为是非晶Ag等温晶化fcc形核时其临界晶核的真实结果. 但在今后对百万级大体系[35]的模拟研究中, 采用本文提出的团簇表征与跟踪方法, 通过对不同晶化阶段各种特征团簇的逆向追踪, 依据其统计分析结果, 应该可更加准确地给出非晶晶化形核的临界晶粒尺寸和结构形态.

    • 1)根据基本团簇种类与联结方式可将扩展晶体团簇分为三类: 以IS方式联结的fcc团簇为fcc单晶团簇, 以FS和IS方式联结的fcc团簇为fcc多晶团簇, 而以IS方式联结的fcc, hcp和bcc团簇则为混晶团簇.

      2)依据团簇结构遗传性特征可识别和表征fcc单晶、多晶与混晶的晶胚与晶核结构. 只有瞬态遗传性而不具有连续遗传性的为晶胚, 具有连续遗传性的称为晶核, 而与连续遗传起始点对应的扩展晶体团簇则为临界晶核.

      3)非晶Ag等温晶化过程中, fcc单晶、多晶与混晶的形核发生在体系能量突降之后与扩展晶体团簇尺寸陡增之前. 不同类型的晶核结构不论临界尺寸还是几何构型都存在差异, 其中fcc单晶临界尺寸最小, 多晶次之, 混晶最大; fcc单晶与多晶壳层原子中有少量hcp和bcc原子, 而混晶壳层则全部为非晶体类原子, 并且fcc单晶、多晶与混晶的临界晶核都不是球型结构.

参考文献 (35)

目录

    /

    返回文章
    返回