对《关于〈一维圆环上双原子链的马德隆常数〉的解析解》的讨论

本文绘制了一维等距离子晶体的累加法、埃夫琴法和平均值法的各级马德隆常数的曲线.本文推导了二聚化直线离子晶体马德隆常数的求和式,比较了两种解析式的优劣.本文详细地推导了圆环上等距和二聚化离子晶体的马德隆常数的解析式,使读者更易于理解.本文修正了《关于〈一维圆环上双原子链的马德隆常数〉的解析解》的两个笔误,还建议修正二聚化马德隆常数的定义式,使其物理意义更明显.     

熔体初始温度对液态金属Ni凝固过程中微观结构演变影响的模拟研究

采用量子 Sutton-Chen多体势, 对熔体初始温度热历史条件对液态金属Ni快速凝固过程中微观结构演变的影响进行了分子动力学模拟研究. 采用双体分布函数g(r)曲线、键型指数法、原子团类型指数法和三维可视化等分析方法对凝固过程中微观结构的演变进行了分析. 结果表明: 熔体初始温度对凝固微结构有显著影响, 但在液态和过冷态时的影响并不明显, 只有在结晶转变温度T_c附近才开始充分显现出来. 体系在1×10¹² K/s的冷速下, 最终均形成以1421和1422键型或面心立方(12 0 0 0 12 0)与六角密集(12 0 0 0 6 6) 基本原子团为主的晶态结构. 末态时, 不同初始温度体系中的主要键型和团簇的数目有很大的变化范围, 且与熔体初始温度的高低呈非线性变化关系. 然而, 体系能量随初始温度呈线性变化关系, 初始温度越高, 末态能量越低, 其晶化程度越高. 通过三维可视化分析进一步发现, 在初始温度较高的体系中, 同类团簇结构的原子出现明显的分层聚集现象, 随着初始温度的下降, 这种分层现象将被弥散开去. 可视化分析将更有助于对凝固过程中微观结构演变进行更为深入的研究. 关键词： 液态金属Ni 熔体初始温度 微观结构 分子动力学模拟     

非晶态金属Al中微团簇多面体结构及其演变规律的模拟研究

对50 000个金属Al原子从液态急冷形成非晶态结构的过程进行了分子动力学模拟研究.采用原子键型和原子团簇类型指数法,发现在其微观结构的转变过程中,与不断增加的1551键型密切相关的二十面体及其缺陷结构的形成起了非常突出的作用.通过各种键型原子团簇的计算,获得团簇多面体结构的分布数据.在可视化的显示下,得到一幅十分清晰的关于团簇多面体结构的特征及其演变的图景.建立了一种简单而有效的团簇多面体结构研究方法.同时也对模拟计算结果的微观机理给出了相应的讨论,这对于深入理解非晶态结构的形成机制及其微观过程,将有重     

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

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

二元液态金属Cu-Ni和Ag-Cu凝固过程的分子动力学模拟研究

用Quantum Sutton-Chen多体势对Ag6Cu4和CuNi液态金属凝固过程进行了分子动力学模拟研究.在冷却速率2×1012到2×1014K/s范围内,CuNi总是形成fcc晶体结构,而Ag6Cu4总是形成非晶态结构.考虑到CuNi及AgCu中原子半径之比分别为1.025和1.13,那么模拟结果证实了原子的尺寸差别是非晶态合金形成的一个主要影响因素.此外采用键对及原子多面体类型指数法对凝固过程中微观结构组态变化的分析,不但能说明二十面体结构在非晶态合金形成和稳定性中所起的关键作用,又有助于对液态金属的凝固过程、非晶态结构特征的深入理解.     

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.     

冷速对液态合金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晶态结构．     

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

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

Polymorph selection of magnesium under different pressures: A simulation study

Molecular dynamics simulations were used to investigate the influence of pressure on the structural properties and dynamics of magnesium(Mg)during rapid solidification.The dynamics analysis revealed that,with an increase in pressure,the dynamics of Mg melt slowed down sharply and the dynamical heterogeneities increased,leading to a denser structure.Atom-level structural analysis using the cluster-type index method suggested that the predominant structure transformed from hexagonal closed-packed to face-centered cubic with increasing pressure from 0 GPa to 5 GPa,and then transformed to the A15 complex crystal structure as the pressure increased above 10 GPa.In addition,the nature of polymorph selection was investigated by analyzing the phonon dispersion of Mg under different pressures.These findings provide a novel insight into polymorphic transitions of Mg under pressure and guide the selection of Mg polymorphs for practical applications.     

冷速对液态金属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之间的确存在一个形成非晶态结构的临界冷速.