冷速对液态金属Ga凝固过程中微观结构演变影响的模拟研究

采用分子动力学方法对液态金属Ga凝固过程中不同冷却速率对微观结构演变的影响进行了模拟跟踪研究. 运用HA键型指数法和原子成团类型指数法(CTIM)分析了金属原子Ga的成键类型和形成的基本原子团结构. 结果发现，冷却速率对凝固过程中的微观结构起着非常重要的作用. 在以1.0×10¹⁴，1.0×10¹³，1.0×10¹²K/s的速率冷却时，系统形成以与1311，1301键型相关的菱面体结构为主，夹杂着立方体、六角密集等其他团簇结构所构成的非晶态结构；在以1.0×10¹¹K/s的速率冷却时，系统明显发生结晶，其结晶转变温度T_c约为198K，且冷却速率越慢，结晶转变温度T_c越高，形成以与1421键型相关的斜方晶体(经可视化分析确认)为主要构型的晶态结构. 这将为研究液态金属的结晶转变过程提供一种新方法. 关键词： 液态金属Ga 凝固过程 微结构转变 分子动力学模拟     

冷速对液态金属Na凝固过程中微观结构影响的模拟研究

采用分子动力学方法对液态金属Na在四种不同冷速下的快速凝固过程进行了模拟跟踪研究.采用双体分布函数g(r)曲线、Honeycutt-Andersen键型指数法和原子团类型指数法对凝固过程中微观结构的变化进行了分析.结果表明：冷却速率对微结构的转变有决定性影响，当冷速为1.0×10¹⁴和1.0×10¹³K/s时，系统形成以1551和1541键型或以缺陷多面体基本原子团(13 1 10 2)和二十面体基本原子团(12 0 12 0)为主体的非晶态结构；当冷速为1.0×10¹²和1.0×10¹¹K/s时，系统则形成以1441和1661键型或以体心立方基本原子团(14 6 0 8)为主体的晶态结构.同时发现：不同冷速对液态金属Na在液态和过冷态时微观结构的影响甚小；但不同冷速对其固态(非晶态利晶态)时的微观结构有显著的影响，且要在液-固转变点(分别在玻璃转变温度T_g和晶化起始温度T_c)附近或以后才能充分展现出来.根据这一特点，有可能建立另一种确定液态金属T_g和T_c的新方法.原子团类型指数法比键型指数法更有利于研究液态、非晶态等无序体系和一些晶化体系的具体结构特征. 关键词： 液态金属Na 凝固过程 分子动力学模拟 原子团类型指数法     

冷速对液态金属Pb凝固过程中微观团簇结构演变影响的模拟研究

采用分子动力学方法对六种不同冷却速率对液态金属Pb凝固过程中微观团簇结构演变的影响进行了模拟跟踪研究.采用双体分布函数、Honeycutt-Andersen键型指数法、原子团类型指数法(CTIM-2)、平均配位数等方法对凝固过程中微观团簇结构的演变进行了分析.结果表明：系统存在一个形成非晶态或晶态结构的临界冷速(介于5×10¹²与1×10¹²之间),大于这个临界冷速时,系统将形成以1551,1541和1431键型为主的非晶态结构;当小于这个临界冷速时,系统将先形成以1441和1661键型或以bcc基本原子团(14 6 0 8 0 0)为主的晶态结构,并稳定存在一段时间,然后又迅速转变为以1421和1422键型为主或以fcc基本原子团(12 0 0 0 12 0)和hcp基本原子团(12 0 0 0 6 6)以一定比例并存的部分晶态结构;同时发现,冷速对系统中fcc结构和hcp结构的相对比例有明显的影响,冷速越低,fcc结构所占的比例越多,越倾向于形成完美的fcc晶态结构. 关键词： 液态金属铅 凝固过程 团簇结构演变 分子动力学模拟     

不同冷速对金属铜Cu凝固过程中微观团簇结构演变影响的模拟研究

采用分子动力学方法和Quantum Sutton-Chen(Q-SC)多体势对由5万个液态金属铜(Cu)原子构成的系统在三个不同冷却速率下的凝固过程中微观团簇结构转变的影响进行了模拟研究.运用双体分布函数、Honeycutt-Andersen(HA)键型指数法、原子团类型指数法(CTIM-2)和可视化分析等方法,对凝固过程中微观团簇结构的演变特性进行了分析研究.结果发现:由非晶体向晶体转变的临界速度约为1.0×1013 K/s,在此冷速下系统形成非晶体和晶体以一定比例并存的混合结构;在冷速为1.0×1014 K/s冷却时系统形成以1551、1541、1532、1431键型为主的非晶体结构,非晶转化温度约为673K;在以4.0×1012 K/s速度冷却时,系统从673K就开始结晶,并形成以1421和1422二种键型或由这二种键型构成的面心立方(FCC)(12 0 0 0 12 0)和六角立方(HCP)(12 0 0 0 6 6)基本原子团为主的晶体结构,尤其是由1421键型构成的面心立方(12 0 0 0 12 0)基本原子团在晶体生长和微观团簇结构形成过程中占主导地位.同时发现,冷速对金属Cu系统中的FCC结构和HCP结构的相对比例有显著的影响,冷速越低,FCC基本原子团以及由其构成的团簇结构越多.     

冷速对液态金属Mg凝固过程中微观结构演变的影响

采用分子动力学方法对不同冷速下液态金属镁(Mg)快速凝固过程中的微观结构演变进行了模拟研究.并采用能量-温度(E-T)曲线、双体分布函数、Honeycutt-Andersen键型指数法、原子团簇类型指数法(CTIM-3)以及三维可视化等方法系统地考察了凝固过程中微观结构演变与相转变过程.结果发现:在以冷速为1×10~(11)K/s的凝固过程中,亚稳态bcc相优先形成,随后大量解体,其变化规律符合Ostwald规则,系统最终形成以hcp结构为主体与fcc结构共存,中间还夹杂部分bcc结构的致密晶体结构.在1×10~(12)K/s冷速下,结晶过程呈现迟缓现象,形成bcc结构的初始温度降低,系统形成以hcp居多、与bcc和fcc三相共存的结构,且因相互竞争、相互制约而导致不易形成粗大的晶粒结构.而在1×10~(13)K/s冷速下,系统则形成以1551,1541,1431键型为主的多种非晶态基本原子团组成的非晶态结构.此外,在冷速1×10~(12)与1×10~(13)K/s之间的确存在一个形成非晶态结构的临界冷速.     

不同冷速对金属铜Cu凝固过程中微观团簇结构演变影响的模拟研究

采用分子动力学方法和Quantum Sutton-Chen（Q-SC）多体势对由5万个液态金属铜（Cu）原子构成的系统在三个不同冷却速率下的凝固过程中微观团簇结构转变的影响进行了模拟研究．运用双体分布函数、Honeycutt-Andersen（HA）键型指数法、原子团类型指数法（CTIM-2）和可视化分析等方法,对凝固过程中微观团簇结构的演变特性进行了分析研究.结果发现：由非晶体向晶体转变的临界速度约为1.0×1013K/s,在此冷速下系统形成非晶体和晶体以一定比例并存的混合结构;在冷速为1.0×1014K/s冷却时系统形成以1551、1541、1532、1431键型为主的非晶体结构,非晶转化温度约为673K;在以4.0×1012K/s速度冷却时,系统从673K就开始结晶,并形成以1421和1422二种键型或由这二种键型构成的面心立方（FCC）（12 0 0 0 12 0）和六角立方（HCP）（12 0 0 0 6 6）基本原子团为主的晶体结构,尤其是由1421键型构成的面心立方（12 0 0 0 12 0）基本原子团在晶体生长和微观团簇结构形成过程中占主导地位. 同时发现,冷速对金属Cu系统中的FCC结构和HCP结构的相对比例有显著的影响,冷速越低,FCC基本原子团以及由其构成的团簇结构越多。     

Ti75Al25合金微观结构演变的分子动力学模拟

本文利用分子动力学模拟方法,研究了液态Ti₇₅Al₂₅合金在不同冷却条件下形成晶体及非晶的过程(Q1:1.0×10¹³K Ks^-1,Q2:1.0×10¹¹Ks^-1).利用平均原子体积、双体分布函数、键角分布函数、键对分析和Voronoi多面体方法研究了微观局域结构随温度的变化关系.研究发现:在Q1冷却过程中,液态Ti₇₅A1₂₅合金在1000 K发生玻璃化转变,形成非晶结构;而在Q2冷却过程中,液态Ti₇₅Al₂₅合金发生结晶,并最终形成hcp晶体结构.     

冷速对液态合金Ca50Zn50快速凝固过程中微观结构演变的影响

采用分子动力学方法对六种不同冷速对原子尺寸相差较大的液态合金Ca50Znso凝固过程中微观结构演变的影响进行了模拟研究,并采用双体分布函数、Honeycutt—Andersen（HA）键型指数法、原子团类型指数法（CTIM一2）、可视化等方法进行了深入分析,结果表明：系统存在一个临界冷速,介于和5×10^11K／s与1×10^11K／s之间,在临界冷速以上（如1×10^11K／s,1×10^11K／s,1×10^11K／s和5×10^11K／s）时,系统形成以1551,1541,1431键型或二十面体基本原子团（12012000）等为主体的非晶态结构;在临界冷速以下时,系统形成以1441和1661键型或bcc基本原子团（1460800）为主体（含有少量的hcp（1200066）和fcc（12000120）基本原子团）的部分晶态结构．在非晶形成的冷速范围内,其总双体分布函数的第一峰明显分裂成与近邻分别为Zn—Zn,Ca—Zn,Ca—Ca相对应的三个次峰;且随着冷速的下降,同类原子近邻的次峰峰值升高、异类原子近邻的次峰峰值下降;Zn原子容易偏聚,随着冷速降低,二十面体的数量增多,非晶态结构也越稳定．在晶态形成的冷速范围内,Zn原子己大量偏聚形成大块bcc晶态结构,Ca原子也部分形成hcp和fcc晶态结构．     

金属铜Cu凝固过程中晶体生长与团簇微观结构演变特性的模拟研究

采用分子动力学方法和Quantum Sutton-Chen（QS-C）多体势,对液态金属铜（Cu）凝固过程中的晶体生长规律及纳米团簇微观结构转变特性进行了模拟跟踪研究．运用Honeycutt-Andersen（HA）键型指数法和新的原子团类型指数法（CTIM-2）分析了金属Cu原子的成键类型和原子团簇结构演变特性.结果发现：当以1.0×1013K/s速率凝固时,系统最终形成晶体和非晶体混合共存结构;在以4.0×1012K/s速度冷却时,系统从673K就开始结晶,并形成以1421和1422二种键型为主的晶体结构;面心立方（FCC）和六角密集（HCP）结构在形成晶体铜时起着非常重要的作用,尤其是由1421键型构成的面心立方（12 0 0 0 12 0）基本原子团在晶体生长和纳米团簇结构形成过程中占主导地位.     

Ni熔体凝固过程中临界晶核和亚临界晶核的分子动力学模拟

本文用分子动力学模拟研究了Ni熔体以不同冷速凝固后微观结构的演变规律, 并通过理论计算确定出了Ni熔体凝固后获得理想非晶的临界条件. 模拟结果发现冷速小于10¹¹ K/s时, Ni 熔体凝固后形成晶态组织; 冷速在10¹¹ K/s到10^14.5 K/s之间时, Ni熔体凝固后形成由晶态结构与非晶态结构组成的混合组织. 冷速小于10¹⁰ K/s, Ni 熔体凝固后形成的晶态组织具有fcc结构; 冷速在10¹⁰ K/s到10^14.5 K/s之间时, Ni熔体凝固后组织中的晶态由fcc和hcp结构层状镶嵌排列构成. 通过分析模拟结果和计算结果, 确定出了Ni熔体凝固后形成理想非晶的临界冷速为10^14.5 K/s. 并发现Ni熔体中临界晶核(冷速等于10^14.5 K/s)和亚临界晶核(冷速大于10^14.5 K/s) 均由fcc和hcp组成层状偏聚结构, 这表明Ni熔体中生长的晶体、临界晶核和晶胚的结构是相同的. 关键词： 分子动力学模拟 晶体团簇 临界冷速 结构     

A simulation study of microstructure evolution during solidification process of liquid metal Ni

A molecular dynamics simulation study has been performed for the microstructure evolution in a liquid metal Ni system during crystallization process at two cooling rates by adopting the embedded atom method (EAM) model potential. The bond-type index method of Honeycutt--Andersen (HA) and a new cluster-type index method (CTIM-2) have been used to detect and analyse the microstructures in this system. It is demonstrated that the cooling rate plays a critical role in the microstructure evolution: below the crystallization temperature $T_{\rm c}$, the effects of cooling rate are very remarkable and can be fully displayed. At different cooling rates of $2.0\times10^{13}$\,K\,$\cdot$\,s$^{-1}$ and $1.0\times10^{12}$\,K\,$\cdot$\,s$^{-1}$, two different kinds of crystal structures are obtained in the system. The first one is the coexistence of the hcp (expressed by (12 0 0 0 6 6) in CTIM-2) and the fcc (12 0 0 0 12 0) basic clusters consisting of 1421 and 1422 bond-types, and the hcp basic cluster becomes the dominant one with decreasing temperature, the second one is mainly the fcc (12 0 0 0 12 0) basic clusters consisting of 1421 bond-type, and their crystallization temperatures $T_{\rm c}$ would be 1073 and 1173\,K, respectively.     

冷速对Cu10Ag90合金凝固结构的影响

通过分子动力学对液态Cu10Ag90合金在四种冷速条件下进行快速凝固模拟。结果显示,1 × 1010和1 × 1011 K/s下系统的平均原子能量分别在750 K和650 K发生突变,冷速越低最终平均原子能量越低;1 × 1012和1 × 1013 K/s的双体分布函数第二峰出现分裂,表明结构处于非晶态。从1 × 1011K/s开始出现尖锐小峰,表明此冷速开始出现晶化现象,1 × 1010 K/s下分裂的峰更加尖锐明显,说明体系形成结晶度较高的晶体结构;1 × 1010 和1 × 1011 K/s下系统凝固后晶体结构含量由高到低分别为fcc, hcp, bcc。冷速越低晶体结构数目越多,系统的有序度更高,结构熵越低。     

Phase and structural transformations in a molecular dynamics model of iron under ultrafast heating and cooling

It is shown that a system of classical particles considered in a molecular dynamics model with Pak-Doyama pairwise interatomic potential adequately describes not only the various structural states of iron (melt, bcc crystal, metal glass) but also the complex self-organization processes occurring in first-and second-order phase transitions (crystallization and vitrification, respectively). When the temperature is varied at a constant rate of 6.6×10¹¹ K/s, crystallization sets in from both the amorphous and the liquid state; at a rate of 1.9×10¹² K/s, crystallization is observed only in the amorphous state; and when heated at a rate of 4.4×10¹² K/s, the model amorphous iron transfers to the liquid state without crystallization. The energy of homogeneous formation of a crystal nucleus in the bulk of the amorphous phase of iron is calculated to be ～0.71 eV under the assumption that there is a spectrum of activation energies.     

液态Mg7Zn3合金快速凝固过程中微观结构演变机理的模拟研究

采用分子动力学方法对液态Mg₇Zn₃合金的快速凝固过程进行了模拟研究,并采用双体分布函数、Honeycutt-Andersen键型指数法、原子团类型指数法(CTIM)以及遗传跟踪等方法对凝固过程中团簇结构的形成、演变特性进行了分析.结果表明:在以冷速为1×10¹² K ·s^-1的凝固过程中, Zn-Zn原子之间的相互成键的概率明显增加,形成以1551,1541,1431键型为主的非晶态结构.其特征键型1551随温度变化 关键词： 液态MgZn合金 快速凝固过程 微观结构演变 分子动力学模拟     

液态金属Ga急冷凝固中微观结构转变的模拟研究

采用分子动力学方法对液态金属Ga的快速凝固过程进行计算机模拟跟踪研究,运用HA键型指数法和原子团类型指数法分析了金属原子Ga的成键类型和形成的基本原子团结构.发现:与通常的液态金属(如Al)相反,随着温度的降低,二十面体及与二十面体相关的1551, 1541键型数目明显下降;与立方体(fcc,bcc),六角密集结构相关的1421,1422键型数目明显增加;而与菱面体相关的1321,1311键型的数目却显著增加,逐渐占据优势.最后形成一种新型的以菱面体结构为主、夹杂着立方体(fcc,bcc)、六角密集(hcp)等团簇结构所组成的非晶态结构.这正是非晶态金属Ga的g(r)曲线分裂的第二峰的顺序为前低后高,而与非晶态金属Al的g(r)曲线(其分裂的第二峰的顺序为前高后低)明显不同的微观结构上的物理根源.