难熔金属钒熔化行为的局域原子结构模拟与分析

蒋元祺

doi:10.7498/aps.69.20200185

摘要

Abstract

作者及机构信息

Authors and contacts

文章全文

1. 引　言
2. 计算方法与参数设置
3. 结果分析与讨论
3.1 体系势能、体积与熔化速率的关系
3.2 径向分布函数与微结构的演化
3.3 不同熔化速率下微结构的分布与演化
4. 结　论

参考文献

施引文献

摘要

采用经典分子动力学(molecular dynamics, MD)方法, 模拟了16000个钒原子在5种不同熔化速率(γ₁ = 1 × 10¹¹ K/s, γ₂ = 1 × 10¹² K/s, γ₃ = 1 × 10¹³ K/s, γ₄ = 1 × 10¹⁴ K/s 与γ₅ = 1 × 10¹⁵ K/s)下原子结构的熔化行为. 结果表明: 不同熔化速率对难熔金属钒的熔点影响明显, 不过随着温度升高, 体系特征原子结构诸如体心立方(BCC)、六角密堆(HCP)、面心立方(FCC)、简单立方(SC)以及二十面体(ICO)的相对分布次序并不随熔化速率的改变而改变, 温度仍然是影响原子结构分布的主要因素. 通过从头算分子动力学(ab initio MD)与热力学分析发现, ICO能够在液态金属区域稳定存在, 一方面是因为其孤立团簇的相对稳定性和团簇寿命要优于晶体型原子团簇, 另一方面是因为其拥有相对较高的团簇熵与相对较低的自由能.

关键词:

分子动力学 /
钒 /
熔化 /
团簇

Abstract

By using large-scale atomic/molecular massively parallel simulator (LAMMPS) code, a molecular dynamics simulation is performed in the NPT ensemble at zero pressure to investigate the influence of melting rates γ on the evolutional characteristics of vanadium atomic structure such as body-centered cubic (BCC), hexagonal close-packed structure (HCP), face centered cubic (FCC), simple cubic (SC) and icosahedra (ICO) during the rapid melting of solid vanadium crystal at five different melting rates (γ₁ = 1 × 10¹¹ K/s, γ₂ = 1 × 10¹² K/s, γ₃ = 1 × 10¹³ K/s, γ₄ = 1 × 10¹⁴ K/s , γ₅ = 1 × 10¹⁵ K/s), in which 16000 atoms in a cubic box under the periodic boundary condition are considered, and their motion equations are solved by Verlet’s algorithm in the velocity form in time steps of 1 fs. Constant pressure P and temperature T are imposed by a modified Nose-Hoover method for both P and T variables, and an embedded-atom model (EAM) potential is utilized. For identifying the local atomic structures of liquid and solid vanadium at different temperatures, a polyhedral template matching method (PTMM) is used by measuring the root-mean square deviation (RMSD), in which clusters are classified as the topology of the local atomic environment without any ambiguity in the classification. Subsequently, the variation of the potential energy, entropy and Gibbs free energy of FCC, HCP, BCC and ICO vanadium clusters are calculated through ab initio MD simulation in the canonical ensemble (NVT) at selected temperatures, and the lowest-energy dynamic structure and its corresponding static heating structure are also shown in this paper. Based on the above calculated results, it is found that the melting point of refractory metal vanadium increases obviously with the increase of heating rate, but the heating rate only presents a limited effect on the population of atomic structure for each of BCC, HCP, FCC, SC and ICO. Namely, the temperature still plays a dominant role in the rapid melting process of V rather than heating rate. Moreover, the ab initio MD simulation and thermodynamics analysis further reveal that lots of ICO clusters of vanadium can exist stably in the liquid region rather than in solid crystal, which is not only due to its higher stability and longer lifetime than those of crystalline atomic clusters, but also because ICO possesses higher entropy and lower Gibbs free energy in high temperature liquid region.

Keywords:

PACS:

36.40.-c	(Atomic and molecular clusters)
36.40.Ei	(Phase transitions in clusters)
61.46.Bc	(Structure of clusters (e.g., metcars; not fragments of crystals; free or loosely aggregated or loosely attached to a substrate))
64.70.dj	(Melting of specific substances)

作者及机构信息

蒋元祺^1,2

1.
南昌师范学院物理系, 南昌　330032

2.
吉林大学物理学院, 长春　130012

通信作者: 蒋元祺, yuanqi325@163.com

基金项目: 江西省自然科学基金(批准号: 20202BAB204004, 20171BAB216001)、江西省教育厅科学技术研究项目(批准号: GJJ191114, GJJ161242)、南昌师范学院首批“青蓝学者”人才计划和国家自然科学基金(批准号: 51871096, 11664028)资助的课题

Authors and contacts

Jiang Yuan-Qi^1,2

1.
Department of Physics, Nanchang Normal University, Nanchang 330032, China

2.
College of Physics, Jilin University, Changchun 130012, China

Corresponding author: Jiang Yuan-Qi, yuanqi325@163.com

Funds: Project supported by the Jiangxi Provincial Natural Science Foundation, China (Grant Nos. 20202BAB204004, 20171BAB216001), the Scientific Research Project of Jiangxi Provincial Education Department, China (Grant Nos. GJJ191114, GJJ161242), the Qinglan Scholars Program of Nanchang Normal University, China, and the National Natural Science Foundation of China (Grant Nos. 51871096, 11664028)

文章全文

1. 引　言

凝固与熔化是自然界中广泛存在的两类自然现象^[1,2], 在工业生产和日常生活中发挥着重要作用, 是研究物质相变和结构演化的有效途径. 通常认为结晶应该只发生在晶化温度附近, 熔化应该只发生在熔点温度处, 但现实中存在“过冷”^[1]与“预熔”^[2]现象, 且在凝固的形核机制与晶体的熔化机制方面仍存在诸多未解之谜^[1,2]. 近年针对液态过渡金属(transition metal, TM)的快凝过程, 发现其中存在着大量的五次对称性原子团与二十面体团簇^[1], 它们不但可以抑制形核^[1-4], 影响金属玻璃的形成能力^[5,6], 而且还与部分畸变的晶体型原子团簇之间存在竞争与演变现象^[7,8]. 2014年Zhong等^[9]利用10¹⁴ K/s的超高冷却速率, 首次在实验上将具有体心立方(BCC)结构的液态难熔单质金属Ta, V, Mo以及W制备成了金属玻璃. 这极大地鼓舞了从事计算模拟的研究人员^[10,11]. 他们以此为基础, 先后从特征原子团簇的竞争与演变^[7]、化学硬度的高低^[8]、近黄金比序^[11]的新型原子结构以及化学序参数的变化^[12]等角度, 相继揭示了液态Cu-Zr合金与单质金属Ag, Ta等在快凝过程中其原子与电子层面的物理机制. 此外还有大量工作深入分析了团簇的遗传^[13,14]、自组织^[15]及其枝状连接^[16], 阐述了团簇结构与玻璃形成能力之间的关系^[17,18], 促进了学界对凝固机理的理解, 进一步深化了对非晶形成的认识^[19,20].

不过遗憾的是, 近二十年来对同样是以温度为主导的自然现象—金属熔化的研究力度远不及对凝固的研究, 绝大部分研究集中在晶体结构的识别^[21]、晶粒生长^[22]、纳米颗粒的表面预熔^[23]、晶格振动^[24]及之后液滴的形成^[25]等, 对其中团簇的研究也主要侧重于结构稳定性^[26,27]、化学序^[28,29]、热力学^[30,31]与动力学行为^[32]、团簇结构竞争^[33]与团簇相转变^[33,34]、不同初始构型^[35]及不同几何尺寸^[36,37]对团簇熔点的影响^[38,39]、金属小团簇熔点与其块体材料熔点的异同^[40,41](包括团簇原子数太少而难以熔化的机理解释^[42])、核-壳^[43]二元反对称Mackay二十面体与扶手型结构之间的转变^[44,45], 以及近两年新兴起的用机器学习^[46]来研究和识别纳米颗粒熔化过程中的原子分布与竞争关系等方面. 而对熔化过程中特征原子团簇的结构演变、分布以及不同熔化速率对微观原子结构的影响涉及较少, 一些具体科学问题仍不清楚^[2,47], 且针对原子团簇的相关学术报道呈逐年下降趋势. 当然这与当前的研究热点有关, 但对熔化现象, 尤其是对其中典型原子团簇的深入研究会更好地促进对一级相变过程的理解. 本文以含有16000个V原子的难熔金属为研究对象, 以300 K为熔化的起始温度, 模拟了其在5种速率下的熔化过程, 初步分析了BCC、面心立方(FCC)、六角密堆(HCP)、二十面体(ICO)以及简单立方(SC)原子团随温度及熔化速率的演变关系, 基本澄清了隐藏在其背后的物理机制.

2. 计算方法与参数设置

将含有16000个钒原子的BCC晶体结构置于一个立方盒子中, 在常温常压(normal pressure and temperature, NPT)系综下, 按周期性边界条件, 采用分子动力学(molecular dynamics, MD)程序LAMMPS^[48]来模拟其熔化过程. 原子间相互作用势采用嵌入原子势(embedded-atom model, EAM)^[49]. 起始熔化温度为300 K, 熔化速率分别为γ₁ = 1 × 10¹¹ K/s, γ₂ = 1 × 10¹² K/s, γ₃ = 1 × 10¹³ K/s, γ₄ = 1 × 10¹⁴ K/s 与γ₅ = 1 × 10¹⁵ K/s, 时间步长为1 fs, 实时输出温度T、势能P、体积V以及原子的位置坐标信息. 在模拟完成后, 利用Ovito程序^[50]与多面体模板匹配法(polyhedral template matching method, PTMM)^[51]进一步可视化分析熔化过程中体系特征原子团簇的构型. 文中涉及标准BCC, FCC, HCP以及ICO的从头算分子动力学(ab initio MD)模拟时, 采用的是正则系综(NVT), 模拟的时间步长为5 fs, 模拟步数为500步, 总模拟时间为2.5 ps.

3. 结果分析与讨论

3.1 体系势能、体积与熔化速率的关系

针对含有16000个钒原子的BCC结构, 图1中首先给出了5种不同熔化速率下势能曲线随温度的演变过程. 从图1(a)与(b)中可以看出, 4种熔化速率(γ₁ = 1 × 10¹¹ K/s, γ₂ = 1 × 10¹² K/s, γ₃ = 1 × 10¹³ K/s, γ₄ = 1 × 10¹⁴ K/s)的能量曲线存在明显的突变, 其熔化转变温度分别为T_m1 = 3246 K, T_m2 = 3309 K, T_m3 = 3384 K以及T_m4 = 3553 K. 而随着熔化速率进一步升高至 $1\!\times \!10^{15}$ K/s, 体系势能曲线的突变消失, 在T = 3250 K附近呈现渐变现象, 直至4236 K与另外4种熔化曲线重合. 即熔化速率越高, T_m也越高. 导致这一现象产生的原因是过高的熔化速率使得原子体系在极短的时间内来不及达到稳定的动力学与热力学平衡. 与此同时, 图2中也同时给出了5种不同熔化速率下体积随温度的演变过程. 从图2知, 在熔化速率γ ≤ γ₄ (1 × 10¹⁴ K/s)时, 不同熔化速率下体系的体积存在明显的突变, 说明有一级相变产生, 其突变所对应的温度点与图1中势能曲线的转变点呈现完美的一一对应关系, 反映出固态钒的体积明显小于液态钒的体积. 当考虑不同熔化速率下体系的压强随温度的波动变化时, 图3中明确显示熔化速率最低时(γ₁ = 1 × 10¹¹ K/s), 压强波动最大, 而当熔化速率越来越高时, 压强的波动幅度则逐渐减小, 当速率升高至1 × 10¹⁵ K/s时, 压强的变化幅度在图3中是最小的. 该现象说明包含16000个钒原子的周期性模拟体系在熔化速率较低时, 研究对象有相对充分的时间来达到相对稳定的亚稳或者稳定平衡状态, 而速率较高时, 体系则经历剧烈的热力学过程.

$图 1 5种不同熔化速率下金属钒体系的势能变化与温度之间的演变关系, (b)是(a)的局部放大图\r\nFig. 1. The potential energy of the vanadium metal as a function of temperature at various melting rates γ, (b) is a partial enlarged view of (a).$

图 1 5种不同熔化速率下金属钒体系的势能变化与温度之间的演变关系, (b)是(a)的局部放大图

Fig. 1. The potential energy of the vanadium metal as a function of temperature at various melting rates γ, (b) is a partial enlarged view of (a).

下载: 全尺寸图片幻灯片

$图 2 5种不同熔化速率下金属钒体系的体积变化与温度之间的演变关系, (b)是(a)的局部放大图\r\nFig. 2. The volume of the vanadium metal as a function of temperature at various melting rates γ, (b) is a partial enlarged view of (a).$

图 2 5种不同熔化速率下金属钒体系的体积变化与温度之间的演变关系, (b)是(a)的局部放大图

Fig. 2. The volume of the vanadium metal as a function of temperature at various melting rates γ, (b) is a partial enlarged view of (a).

下载: 全尺寸图片幻灯片

$图 3 不同熔化速率下体系压强随温度的演变关系\r\nFig. 3. The pressure of the vanadium metal as a function of temperature at various melting rates γ.$

图 3 不同熔化速率下体系压强随温度的演变关系

Fig. 3. The pressure of the vanadium metal as a function of temperature at various melting rates γ.

下载: 全尺寸图片幻灯片

3.2 径向分布函数与微结构的演化

针对熔化过程中系统微结构的演变, 径向分布函数(radial distribution function, RDF)是一种简单有效的方法, 与实验上经常使用的结构因子互为傅里叶变换, 是表征结构的一个重要参数, 数学定义如下^[1]:

${g_{{\rm{\alpha \beta }}}}(r) = \dfrac{N}{{4{\rm{\pi }}\rho {N_{\rm{\alpha }}}{N_{\rm{\beta }}}}}\sum\limits_{i = 1}^{{N_{\rm{a}}}} {\sum\limits_{j = 1}^{{N_{\rm{\beta }}}} {\delta \left( {r - \left| {{{{r}}_{ij}}} \right|} \right)} },$

(1)

其 ${g_{{\rm{\alpha \beta }}}}(r)$ 表示以所有 ${\rm{\alpha }}$ 粒子为中心, 在距中心原子r处的球壳内的单位体积中找到 ${\rm{\beta }}$ 粒子的平均几率. N是粒子总数, ${N_{\rm{\alpha }}}$ 和 ${N_{\rm{\beta }}}$ 分别是 ${\rm{\alpha }}$ 类和 ${\rm{\beta }}$ 类原子的数目, $\left| {{{{r}}_{ij}}} \right|$ 是第i个原子( ${\rm{\alpha }}$ 类)和第j ( ${\rm{\beta }}$ 类)个原子间的距离, $\rho$ 是原子数密度. 对于不区分原子种类的单原子系统, 其表达式简化如下^[1]:

$g(r) = \dfrac{1}{{4{\rm{\pi }}\rho N}}\sum\limits_{i = 1}^N {\sum\limits_{j = 1,j \ne i}^N {\delta (r - \left| {{{{r}}_{ij}}} \right|)} } .$

(2)

结合上述RDF的定义以及文献[1]中的示意图, 图4中给出了金属钒以1 × 10¹² K/s的速率熔化时的RDF曲线图. 可知, RDF曲线中峰和谷的变化过程可以分为如下四个阶段: 首先在低温300 K时, 明锐分立的衍射峰反映出完美的晶体结构特征; 在随后的400—3200 K (温度间隔ΔT = 100 K)的温度区间内, 随着温度的逐渐升高, 分立尖锐的峰逐渐变得模糊, 尤其是在r = 5 Å处的两个较窄的、分立的峰逐渐合并为一个较宽的峰, 即300 K时的第二峰逐渐减小并融入第三峰; 当温度继续升高, 在3200 K ≤ T ≤ 3320 K的温度区间内, 除第一峰外, RDF曲线其余峰变得更加平缓, 但仍未完全消失, 呈现出典型的固态向液态熔化演变的过渡特征; 当温度继续升高至3320 K以上, 在r ≥ 5 Å的区域内RDF系列峰逐渐消失, 呈现液态金属的典型特征. RDF曲线的峰与谷的演变规律与势能曲线与温度之间的演变规律具有很好的对应关系.

$图 4 熔化速率为γ2 = 1 × 1012 K/s时RDF在不同温度区间的演变关系\r\nFig. 4. The RDF for vanadium metal at different temperature obtained by rapid heating rate γ2 = 1 × 1012 K/s.$

图 4 熔化速率为γ₂ = 1 × 10¹² K/s时RDF在不同温度区间的演变关系

Fig. 4. The RDF for vanadium metal at different temperature obtained by rapid heating rate γ₂ = 1 × 10¹² K/s.

下载: 全尺寸图片幻灯片

在此基础上, 图5更加直观地给出了钒晶体熔化过程中不同类型原子团簇(FCC, HCP, BCC, ICO以及SC)的演变示意图. 可知300 K时, 体系的确呈现完美的晶体结构, 原子排列整齐, 然而当温度继续升高至3300 K时, 原子已经明显偏离了其平衡位置, BCC显著减少, 之后当温度每升高10 K其结构变化都是异常剧烈的. 当T升高至4500 K时, 体系的原子排布已经是完全“无序”的排列, 直观地印证了图4中RDF反映的体系结构特征. 与此同时, 也要承认以速率1 × 10¹² K/s进行的熔化所得的熔点温度与实验熔化温度有一定差距, 这主要是由于熔化速率过高所致.

$图 5 晶体钒熔化过程不同类型团簇的演变示意图(γ2 = 1 × 1012 K/s)\r\nFig. 5. Different type atomic schematic diagraph in vanadium system at heating rate γ2 = 1 × 1012 K/s.$

图 5 晶体钒熔化过程不同类型团簇的演变示意图(γ₂ = 1 × 10¹² K/s)

Fig. 5. Different type atomic schematic diagraph in vanadium system at heating rate γ₂ = 1 × 10¹² K/s.

下载: 全尺寸图片幻灯片

3.3 不同熔化速率下微结构的分布与演化

在完成5种不同熔化速率的MD模拟之后, 要进一步获悉熔化过程中晶体型原子团簇之间的演变和彼此之间的相互转化, 以及要发生这种转化所需的基本条件等, 首先就需要对其中承载基本信息的原子团簇进行识别. 本文将运用多面体模板匹配法(polyhedral template matching method, PTMM)^[51]来完成这一任务. 简单来说, PTMM方法是基于均方根偏差(root-mean-square deviation, RMSD)的数学方法来判断两点之间结构的相似性. 对于给定的两个原子A与B, 其RMSD可定义为

${\rm{RMSD}}({{A}},{{B}}) = \sqrt {\frac{1}{N}\sum\limits_{i = 1}^N {\left\| {{{{A}}_i} - {{{B}}_j}} \right\|_{}^2} },$

(3)

而对于两个以上的多原子体系的叠加问题, 主要是要找到一个旋转矢量A与一个缩放矢量B, 使得RMSD最小. 这种情况下RMSD的定义式转变如下^[51]:

$\begin{split} &{\rm{RMSD}}({{A}},{{B}}) = \\ & \mathop {{\rm{min}}}\limits_{s {{Q}}}\sqrt {\frac{1}{N}\sum\limits_{i = 1}^N {{{\Big\| {s\left[ {{{{a}}_i} - {{A}}} \right] - {{\big( {{{Q}}{{[{{{b}}_i}-{{B}}]}^{\rm{T}}}} \big)}^{\rm{T}}}} \Big\|}^2}} },\end{split}$

(4)

其中N是原子数, s是A的最佳缩放因子, Q是右手正交矩阵, ${{A}} = \dfrac{1}{N}\displaystyle\sum\limits_{i = 1}^N {{{ {{a}}}_i}}$ 与 ${{B}} = \dfrac{1}{N}\displaystyle\sum\limits_{i = 1}^N {{{ {{b}}}_i}}$ 分别是A与B原子的重心, ${ {{a}}_i}$ 与 ${ {{b}}_i}$ 分别代表A与B原子的任意变化矢量. 通常可以利用正交矩阵与缩放因子的共同调节使得RMSD可以保持一个测定的不变量, 从而获得变换后各点距原点的距离^[51]:

$\begin{split} &{\rm{RMSD}}({{A}},{{B}}) = \\ &\mathop {{\rm{min}}}\limits_{s{{Q}}} \sqrt {\dfrac{1}{N}\sum\limits_{i = 1}^N {{{\Big\| {s\left[ {{{{a}}_i} - {{A}}} \right] - \dfrac{1}{{l({{B}})}}{{\big( {{{Q}}{{[{b_i} - {{B}}]}^{\rm{T}}}} \big)}^{\rm{T}}}} \Big\|}^2}} } , \end{split}$

(5)

其中 $l( {{B}}) = \dfrac{1}{N}\displaystyle\sum\nolimits_{i = 1}^N {\left\| {{b_i} - {{B}}} \right\|}$ . 之后确定一个中心原子和其近邻原子之间的结构, 就只需先给出标准模板团簇构型(如图5中的ICO, BCC, HCP, FCC及SC多面体), 然后通过RMSD来判断该标准团簇模板与目标团簇构型之间的结构相似性. RMSD数值越小, 条件越发苛刻, 目标结构与模板的相似度也就越高, 所表征的团簇结构也就越接近于模板结构, 该方法称之为PTMM, 目前已被广泛应用于对非晶与玻璃态结构的表征, 接下来将用其分析不同类型团簇的原子结构.

众所周知, 金属钒在常温下是BCC结构, 图5中已经初步给出了FCC, HCP, BCC, ICO以及SC团簇原子随温度的变化示意图, 在此基础上, 图6给出了5种不同熔化速率(γ₁ = 1 × 10¹¹ K/s, γ₂ = 1 × 10¹² K/s, γ₃ = 1 × 10¹³ K/s, γ₄ = 1 × 10¹⁴ K/s以及γ₅ = 1 × 10¹⁵ K/s)下FCC, HCP, BCC, ICO以及SC团簇原子数随温度变化的具体情况. 从图6总的来看, 在低温区间内, 不论是高熔化速率还是低熔化速率, 当温度由室温逐渐升高时, BCC团簇在起始阶段呈现缓慢减少的趋势, 但当温度逐渐逼近T_m时, BCC原子急剧减少, HCP则迅速增加, 之后ICO, FCC以及SC等原子团簇也在T = T_m这一瞬间陡然增加. 而一旦当温度T > T_m时, HCP与ICO则由最大值开始急剧下降, 同时FCC也开始减少. 与此同时虽然SC的数量在高温阶段较少, 但其数目并未随温度的增加而发生明显变化. 简言之, 图6中隐藏着如下的物理规律: 即在BCC急剧消失的时候, 在熔化温度附近有部分原子团簇快速转化为HCP, ICO, FCC以及SC团簇, 这一熔化规律并不随熔化速率的升高而发生改变. 之后为了进一步了解晶体钒在熔化过程中, 不同熔化速率对同一类型团簇原子数目的影响, 图7分别给出了FCC, HCP, BCC, ICO, SC以及其他类型团簇的原子分布随温度的演变关系. 从图7可见, 不论是FCC, HCP, ICO还是SC类型的原子团簇, 其数目都在熔化温度T_m附近明显增加, 与此同时伴随着BCC的急剧减少(图7(c)). 对同一团簇而言, 考虑熔化速率效应, 发现在速率γ ≤ 1 × 10¹² K/s时, 速率的改变对团簇数目的影响并不显著. 当γ大于1 × 10¹² K/s时, 熔化速率越高, BCC类型团簇减少得越慢, 而FCC, HCP, ICO以及SC则增长得越慢. 之后随着温度逐渐升高, FCC, HCP以及ICO的数目在达到顶峰后才逐渐减少.

$图 6 利用多面体模板匹配法分析得到的5种不同熔化速率下FCC, HCP, BCC, ICO以及SC类型团簇原子随温度的演变关系\r\nFig. 6. The population of the FCC, HCP, BCC, ICO and SC types atoms in V system as a function of temperature obtained from the polyhedral template matching at various rates, respectively.$

图 6 利用多面体模板匹配法分析得到的5种不同熔化速率下FCC, HCP, BCC, ICO以及SC类型团簇原子随温度的演变关系

Fig. 6. The population of the FCC, HCP, BCC, ICO and SC types atoms in V system as a function of temperature obtained from the polyhedral template matching at various rates, respectively.

下载: 全尺寸图片幻灯片

$图 7 利用多面体模板匹配法分析得到的5种不同熔化速率下各种类型团簇原子分布随温度的演变关系　(a) FCC; (b) HCP; (c) BCC; (d) ICO; (e) SC; (f)其他类型团簇\r\nFig. 7. The fraction of the various types atoms in V system as a function of temperature obtained from the polyhedral template matching at five various rates: (a) FCC; (b) HCP; (c) BCC; (d) ICO; (e) SC; (f) other types atoms counts.$

图 7 利用多面体模板匹配法分析得到的5种不同熔化速率下各种类型团簇原子分布随温度的演变关系　(a) FCC; (b) HCP; (c) BCC; (d) ICO; (e) SC; (f)其他类型团簇

Fig. 7. The fraction of the various types atoms in V system as a function of temperature obtained from the polyhedral template matching at five various rates: (a) FCC; (b) HCP; (c) BCC; (d) ICO; (e) SC; (f) other types atoms counts.

下载: 全尺寸图片幻灯片

针对图6与图7中 ICO原子团簇在高温阶段持续存在的物理现象, 图8进一步利用ab initio MD模拟了BCC(15个原子)、FCC(13个原子)、HCP(13个原子)以及ICO(13个原子)分别在300, 500 K和1100 K时的结构稳定性与团簇寿命. 结果显示在300 K和500 K时, FCC与HCP均在模拟结束后转变为了ICO, BCC也在模拟结束后转变为了含有ICO团簇核心的几何结构, 而ICO则几乎保持其初始结构而未发生变化. 当温度升高至1100 K时, BCC演变为一个含有畸变ICO的结构, 而FCC与HCP则经过一系列不同的亚稳结构最终演变为非规则构型. 唯独ICO仍然保持其基本构型, 且未发生明显改变, 这说明在高温阶段ICO的稳定性的确高于晶体型原子团簇, 这也为ICO在温度高于熔化温度后的大量出现提供了团簇层面的解释.

$图 8 FCC, HCP, BCC以及ICO 团簇平均每原子势能随模拟步长的演变趋势　(a) 300 K; (b) 500 K; (c) 1100 K\r\nFig. 8. The variation of the potential energy of per atom of FCC, HCP, BCC and ICO cluster as a function of time step, respectively: (a) 300 K; (b) 500 K; (c) 1100 K$

图 8 FCC, HCP, BCC以及ICO 团簇平均每原子势能随模拟步长的演变趋势　(a) 300 K; (b) 500 K; (c) 1100 K

Fig. 8. The variation of the potential energy of per atom of FCC, HCP, BCC and ICO cluster as a function of time step, respectively: (a) 300 K; (b) 500 K; (c) 1100 K

下载: 全尺寸图片幻灯片

此外, 考虑到熔化过程中的MD模拟^[35]与实验^[37,40]均是在有限温度下进行的, 所以本文通过计算团簇的原子振动频率^[52]后, 进一步从熵与吉布斯自由能的角度来研究各团簇的热力学稳定性^[52], 期望以此找到具有5次对称性而缺乏平移对称性的ICO能够在高温区稳定存在的原始因素. 因为熵与吉布斯自由能均是描述热力学系统的重要态函数, 其大小均可清晰反映系统所处状态的稳定情况, 指明热力学过程进行的方向, 可以为研究对象提供定量表述, 即熵越高, 吉布斯自由能越低, 孤立体系越稳定. 基于此, 图9中分别给出了FCC, HCP, ICO及BCC(为了与其他三种标准团簇具有可比性，此处是仅含13原子的标准BCC碎片)原子团簇的熵与吉布斯自由能随温度的演变关系. 可知, 虽然各团簇的熵与自由能随温度的演变趋势恰好相反, 但其所指示的物理意义^[52]却完全一致, 即与FCC, HCP, BCC相比, ICO的熵最高, 吉布斯自由能最低, 而且随着温度逐渐升高, 该现象愈加显著. 这充分表明高温液态金属区域ICO的热力学稳定性远远优于各晶体型原子团, 为液态金属中ICO原子团的大量存在提供了热力学层面的理论解释.

$图 9 FCC, HCP, BCC及ICO 原子团簇的熵与吉布斯自由能随温度的演变关系　(a) 熵; (b) 吉布斯自由能\r\nFig. 9. The variation of the entropy and Gibbs free energy of FCC, HCP, BCC and ICO cluster as a function of temperature, respectively: (a) Entropy; (b) Gibbs free energy.$

图 9 FCC, HCP, BCC及ICO 原子团簇的熵与吉布斯自由能随温度的演变关系　(a) 熵; (b) 吉布斯自由能

Fig. 9. The variation of the entropy and Gibbs free energy of FCC, HCP, BCC and ICO cluster as a function of temperature, respectively: (a) Entropy; (b) Gibbs free energy.

下载: 全尺寸图片幻灯片

4. 结　论

本文采用MD方法, 通过LAMMPS模拟了5种不同熔化速率(γ₁ = 1 × 10¹¹ K/s, γ₂ = 1 × 10¹² K/s, γ₃ = 1 × 10¹³ K/s, γ₄ = 1 × 10¹⁴ K/s与γ₅ = 1 × 10¹⁵ K/s)下16000个钒原子在同一初始几何结构下局域特征原子结构的熔化行为. 研究表明: 在不同速率的熔化过程中, 钒的熔点随着速率的升高而依次升高, 熔化速率对不同类型原子团簇的绝对数目影响显著, 但对不同类型原子团簇数目分布的相对次序无显著影响. 针对特征原子团簇进行的ab initio MD模拟和热力学分析表明, ICO原子团簇在高温区之所以能够持续存在, 不仅是因为ICO原子团簇的稳定性和团簇寿命明显优于FCC, HCP以及BCC, 而且是因为ICO原子团簇在高温区具有相对较高的团簇熵与相对较低的吉布斯自由能.

参考文献

[1]	Cheng Y Q, Ma E 2011 Prog. Mater. Sci. 56 379 Google Scholar
[2]	Baletto F 2005 Rev. Mod. Phys. 77 371 Google Scholar
[3]	Turci F, Tarjus G, Royall C P 2017 Phys. Rev. Lett. 118 215501 Google Scholar
[4]	Shen Y T, Kim T H, Gangopadhyay A K, Kelton K F 2009 Phys. Rev. Lett. 102 057801 Google Scholar
[5]	Sprakel J, Zaccone A, Spaepen F, Schall P, Weitz D 2017 Phys. Rev. Lett. 118 088003 Google Scholar
[6]	Şopu D, Stukowski A, Stoica M, Scudino S 2017 Phys. Rev. Lett. 119 195503 Google Scholar
[7]	Jiang Y Q, Wen D D, Peng P 2017 J. Mole. Liquids 230 271 Google Scholar
[8]	Jiang Y Q, Wen D D, Peng P, Han S C, Hou Z Y 2015 Comput. Mater. Sci. 99 156 Google Scholar
[9]	Zhong L, Wang J, Sheng H, Zhang Z, Mao S X 2014 Nature 512 177 Google Scholar
[10]	蒋元祺, 彭平 2018 物理学报 67 132101 Google Scholar Jiang Y Q, Peng P 2018 Acta Phys. Sin. 67 132101 Google Scholar
[11]	Jiang Y Q, Peng P 2020 Chin. Phys. B. 29 046105 Google Scholar
[12]	Jiang Y Q, Wen D D, He W X, Peng P 2018 Mol. Simulat. 44 1183 Google Scholar
[13]	武振伟, 李茂枝, 徐莉梅, 汪卫华 2017 物理学报 66 176405 Google Scholar Wu Z W, Li M Z, Xu L M, Wang W H 2017 Acta Phys. Sin. 66 176405 Google Scholar
[14]	邓永和, 文大东, 彭超, 韦彦丁, 赵瑞, 彭平 2016 物理学报 65 066401 Google Scholar Deng Y H, Wen D D, Peng C, Wei Y D, Zhao R, Peng P 2016 Acta Phys. Sin. 65 066401 Google Scholar
[15]	孙保安, 王利峰, 邵建华 2017 物理学报 66 178103 Google Scholar Sun B A, Wang L F, Shao J H 2017 Acta Phys. Sin. 66 178103 Google Scholar
[16]	Yang M H, Cai B, Sun Y, Zhang F, Wang Y F, Wang C Z, Ho K M 2019 Phys. Rev. Mater. 3 125602 Google Scholar
[17]	孙奕韬, 王超, 吕玉苗, 胡远超, 罗鹏, 刘明, 咸海杰, 赵德乾, 丁大伟, 孙保安, 潘明祥, 闻平, 白海洋, 柳延辉, 汪卫华 2018 物理学报 67 126101 Google Scholar Sun Y T, Wang C, Lü Y M, Hu Y C, Luo P, Liu M, Xian H J, Zhao D Q, Ding D W, Sun B A, Pan M X, Wen P, Bai H Y, Liu Y H, Wang W H 2018 Acta Phys. Sin. 67 126101 Google Scholar
[18]	Li M X, Zhao S F, Lu Z, Hirata A, Wen P, Bai H Y, Chen M W, Schroers J, Liu Y H, Wang W H 2019 Nature 569 99 Google Scholar
[19]	Sheng H W, Luo W K, Alamgir F M, Bai J M, Ma E 2006 Nature 439 419 Google Scholar
[20]	Hirata A, Kang L J, Fujita T, Klumov B, Matsue K, Kotani M, Yavari A R, Chen M W 2013 Science 341 376 Google Scholar
[21]	Wang H Q, Li H F 2015 RSC Adv. 5 94685 Google Scholar
[22]	Sosso C C, Chen J, Stephen J, Fitzner M, Pedevilla P, Zen A, Michaelides A 2016 Chem. Rev. 116 7078 Google Scholar
[23]	Taherkhani F, Akbarzadeh H, Rezania H 2014 J. Alloys. Compd. 617 746 Google Scholar
[24]	Chen L, Wang Q, Xiong L 2017 J. Nanopart. Res. 19 300 Google Scholar
[25]	Zhang Z, Hu W, Xiao S 2006 Phys. Rev. B 73 125443 Google Scholar
[26]	Nunez S, Lopez J M, Aguado A 2012 Nanoscale 4 6481 Google Scholar
[27]	Steenbergen K G, Gaston N 2013 Phys. Chem. Chem. Phys. 15 15325 Google Scholar
[28]	Taran S, Garip A K, Arslan H 2020 J. Clus. Sci. 10876 Google Scholar
[29]	Taran S, Garip A K, Arslan H 2020 Chin. Phys. B 29 077801 Google Scholar
[30]	Elatresh S F, Bonev S A, Gregoryanz E, Ashcroft N W 2016 Phys. Rev. B 94 104107 Google Scholar
[31]	Zhao C Y, Tao Y, Yu Y S 2020 Int. J. Mass. Tran. 150 119382 Google Scholar
[32]	Feng D, Feng Y, Yuan S, Zhang X, Wang G 2017 Appl. Therm. Eng. 111 1457 Google Scholar
[33]	Fang X W, Wang C Z, Yao Y X, Ding Z J, Ho K M 2011 Phys. Rev. B 83 224203 Google Scholar
[34]	Aguado A, Jarrold M F 2011 Annu. Rev. Phys. Chem. 62 151 Google Scholar
[35]	Lekka C E, Papaconstantopoulos D A 2007 Surf. Sci. 601 3937 Google Scholar
[36]	Pyfer K L, Kafader J O, Yalamanchali A, Jarrold M F 2014 J. Phys. Chem. A 118 4900 Google Scholar
[37]	Haberland H, Hippler T, Donges J, Kostko O, Schmidt M, Issendorff B 2005 Phys. Rev. Lett. 94 035701 Google Scholar
[38]	Zhang L, Sun H 2010 Phys. Status Solidi 207 1178 Google Scholar
[39]	Jena P, Sun Q 2018 Chem. Rev. 118 5755 Google Scholar
[40]	Chacko S, Joshi K, Kanhere D G 2004 Phys. Rev. Lett. 92 135506 Google Scholar
[41]	Rapacioli M, Tarrat N, Spiegelman F 2018 J. Phys. Chem. A 122 4092 Google Scholar
[42]	Breaux G.A, Benirschke R C, Sugai T, Kinnear B S, Jarrold M F 2003 Phys. Rev. Lett. 91 215508 Google Scholar
[43]	Nelli D, Ferrando R 2019 Nanoscale 11 13040 Google Scholar
[44]	Settem M, Kanjarla A K 2020 Sci. Rep. 10 3296 Google Scholar
[45]	Wu J, Qi L, You H, Gross A, Li J, Yang H 2012 J. Am. Chem. Soc. 134 11880 Google Scholar
[46]	Zhang J, Chen J, Hu P, Wang H 2020 Chin. Chem. Lett. 31 890 Google Scholar
[47]	Mao H K, Chen X J, Ding Y, Li B, Wang B 2018 Rev. Mod. Phys. 90 015007 Google Scholar
[48]	Plimpton S 1995 J. Comput. Phys. 117 1 Google Scholar
[49]	Cai J, Ye Y 1996 Phys. Rev. B 54 8398 Google Scholar
[50]	Stukowski A 2010 Model. Simul. Mater. Sc. 18 015012 Google Scholar
[51]	Larsen P M, Schmidt S, Schiøtz J 2016 Model. Simul. Mater. Sc. 24 055007 Google Scholar
[52]	蒋元祺 2015 博士学位论文(长沙: 湖南大学) Jiang Y Q 2015 Ph. D. Dissertation (Changsha: Hunan University) (in Chinese)

施引文献

期刊类型引用(1)

王海鹏，王庆. 表征液态金属结构双体分布函数的讲授方法. 物理与工程. 2022(01): 155-161 .

百度学术

其他类型引用(0)

图 1 5种不同熔化速率下金属钒体系的势能变化与温度之间的演变关系, (b)是(a)的局部放大图

Fig. 1. The potential energy of the vanadium metal as a function of temperature at various melting rates γ, (b) is a partial enlarged view of (a).

下载: 全尺寸图片幻灯片

图 2 5种不同熔化速率下金属钒体系的体积变化与温度之间的演变关系, (b)是(a)的局部放大图

Fig. 2. The volume of the vanadium metal as a function of temperature at various melting rates γ, (b) is a partial enlarged view of (a).

下载: 全尺寸图片幻灯片

图 3 不同熔化速率下体系压强随温度的演变关系

Fig. 3. The pressure of the vanadium metal as a function of temperature at various melting rates γ.

下载: 全尺寸图片幻灯片

图 4 熔化速率为γ₂ = 1 × 10¹² K/s时RDF在不同温度区间的演变关系

Fig. 4. The RDF for vanadium metal at different temperature obtained by rapid heating rate γ₂ = 1 × 10¹² K/s.

下载: 全尺寸图片幻灯片

图 5 晶体钒熔化过程不同类型团簇的演变示意图(γ₂ = 1 × 10¹² K/s)

Fig. 5. Different type atomic schematic diagraph in vanadium system at heating rate γ₂ = 1 × 10¹² K/s.

下载: 全尺寸图片幻灯片

图 6 利用多面体模板匹配法分析得到的5种不同熔化速率下FCC, HCP, BCC, ICO以及SC类型团簇原子随温度的演变关系

Fig. 6. The population of the FCC, HCP, BCC, ICO and SC types atoms in V system as a function of temperature obtained from the polyhedral template matching at various rates, respectively.

下载: 全尺寸图片幻灯片

图 7 利用多面体模板匹配法分析得到的5种不同熔化速率下各种类型团簇原子分布随温度的演变关系　(a) FCC; (b) HCP; (c) BCC; (d) ICO; (e) SC; (f)其他类型团簇

Fig. 7. The fraction of the various types atoms in V system as a function of temperature obtained from the polyhedral template matching at five various rates: (a) FCC; (b) HCP; (c) BCC; (d) ICO; (e) SC; (f) other types atoms counts.

下载: 全尺寸图片幻灯片

图 8 FCC, HCP, BCC以及ICO 团簇平均每原子势能随模拟步长的演变趋势　(a) 300 K; (b) 500 K; (c) 1100 K

Fig. 8. The variation of the potential energy of per atom of FCC, HCP, BCC and ICO cluster as a function of time step, respectively: (a) 300 K; (b) 500 K; (c) 1100 K

下载: 全尺寸图片幻灯片

图 9 FCC, HCP, BCC及ICO 原子团簇的熵与吉布斯自由能随温度的演变关系　(a) 熵; (b) 吉布斯自由能

Fig. 9. The variation of the entropy and Gibbs free energy of FCC, HCP, BCC and ICO cluster as a function of temperature, respectively: (a) Entropy; (b) Gibbs free energy.

下载: 全尺寸图片幻灯片

[1]	Cheng Y Q, Ma E 2011 Prog. Mater. Sci. 56 379 Google Scholar
[2]	Baletto F 2005 Rev. Mod. Phys. 77 371 Google Scholar
[3]	Turci F, Tarjus G, Royall C P 2017 Phys. Rev. Lett. 118 215501 Google Scholar
[4]	Shen Y T, Kim T H, Gangopadhyay A K, Kelton K F 2009 Phys. Rev. Lett. 102 057801 Google Scholar
[5]	Sprakel J, Zaccone A, Spaepen F, Schall P, Weitz D 2017 Phys. Rev. Lett. 118 088003 Google Scholar
[6]	Şopu D, Stukowski A, Stoica M, Scudino S 2017 Phys. Rev. Lett. 119 195503 Google Scholar
[7]	Jiang Y Q, Wen D D, Peng P 2017 J. Mole. Liquids 230 271 Google Scholar
[8]	Jiang Y Q, Wen D D, Peng P, Han S C, Hou Z Y 2015 Comput. Mater. Sci. 99 156 Google Scholar
[9]	Zhong L, Wang J, Sheng H, Zhang Z, Mao S X 2014 Nature 512 177 Google Scholar
[10]	蒋元祺, 彭平 2018 物理学报 67 132101 Google Scholar Jiang Y Q, Peng P 2018 Acta Phys. Sin. 67 132101 Google Scholar
[11]	Jiang Y Q, Peng P 2020 Chin. Phys. B. 29 046105 Google Scholar
[12]	Jiang Y Q, Wen D D, He W X, Peng P 2018 Mol. Simulat. 44 1183 Google Scholar
[13]	武振伟, 李茂枝, 徐莉梅, 汪卫华 2017 物理学报 66 176405 Google Scholar Wu Z W, Li M Z, Xu L M, Wang W H 2017 Acta Phys. Sin. 66 176405 Google Scholar
[14]	邓永和, 文大东, 彭超, 韦彦丁, 赵瑞, 彭平 2016 物理学报 65 066401 Google Scholar Deng Y H, Wen D D, Peng C, Wei Y D, Zhao R, Peng P 2016 Acta Phys. Sin. 65 066401 Google Scholar
[15]	孙保安, 王利峰, 邵建华 2017 物理学报 66 178103 Google Scholar Sun B A, Wang L F, Shao J H 2017 Acta Phys. Sin. 66 178103 Google Scholar
[16]	Yang M H, Cai B, Sun Y, Zhang F, Wang Y F, Wang C Z, Ho K M 2019 Phys. Rev. Mater. 3 125602 Google Scholar
[17]	孙奕韬, 王超, 吕玉苗, 胡远超, 罗鹏, 刘明, 咸海杰, 赵德乾, 丁大伟, 孙保安, 潘明祥, 闻平, 白海洋, 柳延辉, 汪卫华 2018 物理学报 67 126101 Google Scholar Sun Y T, Wang C, Lü Y M, Hu Y C, Luo P, Liu M, Xian H J, Zhao D Q, Ding D W, Sun B A, Pan M X, Wen P, Bai H Y, Liu Y H, Wang W H 2018 Acta Phys. Sin. 67 126101 Google Scholar
[18]	Li M X, Zhao S F, Lu Z, Hirata A, Wen P, Bai H Y, Chen M W, Schroers J, Liu Y H, Wang W H 2019 Nature 569 99 Google Scholar
[19]	Sheng H W, Luo W K, Alamgir F M, Bai J M, Ma E 2006 Nature 439 419 Google Scholar
[20]	Hirata A, Kang L J, Fujita T, Klumov B, Matsue K, Kotani M, Yavari A R, Chen M W 2013 Science 341 376 Google Scholar
[21]	Wang H Q, Li H F 2015 RSC Adv. 5 94685 Google Scholar
[22]	Sosso C C, Chen J, Stephen J, Fitzner M, Pedevilla P, Zen A, Michaelides A 2016 Chem. Rev. 116 7078 Google Scholar
[23]	Taherkhani F, Akbarzadeh H, Rezania H 2014 J. Alloys. Compd. 617 746 Google Scholar
[24]	Chen L, Wang Q, Xiong L 2017 J. Nanopart. Res. 19 300 Google Scholar
[25]	Zhang Z, Hu W, Xiao S 2006 Phys. Rev. B 73 125443 Google Scholar
[26]	Nunez S, Lopez J M, Aguado A 2012 Nanoscale 4 6481 Google Scholar
[27]	Steenbergen K G, Gaston N 2013 Phys. Chem. Chem. Phys. 15 15325 Google Scholar
[28]	Taran S, Garip A K, Arslan H 2020 J. Clus. Sci. 10876 Google Scholar
[29]	Taran S, Garip A K, Arslan H 2020 Chin. Phys. B 29 077801 Google Scholar
[30]	Elatresh S F, Bonev S A, Gregoryanz E, Ashcroft N W 2016 Phys. Rev. B 94 104107 Google Scholar
[31]	Zhao C Y, Tao Y, Yu Y S 2020 Int. J. Mass. Tran. 150 119382 Google Scholar
[32]	Feng D, Feng Y, Yuan S, Zhang X, Wang G 2017 Appl. Therm. Eng. 111 1457 Google Scholar
[33]	Fang X W, Wang C Z, Yao Y X, Ding Z J, Ho K M 2011 Phys. Rev. B 83 224203 Google Scholar
[34]	Aguado A, Jarrold M F 2011 Annu. Rev. Phys. Chem. 62 151 Google Scholar
[35]	Lekka C E, Papaconstantopoulos D A 2007 Surf. Sci. 601 3937 Google Scholar
[36]	Pyfer K L, Kafader J O, Yalamanchali A, Jarrold M F 2014 J. Phys. Chem. A 118 4900 Google Scholar
[37]	Haberland H, Hippler T, Donges J, Kostko O, Schmidt M, Issendorff B 2005 Phys. Rev. Lett. 94 035701 Google Scholar
[38]	Zhang L, Sun H 2010 Phys. Status Solidi 207 1178 Google Scholar
[39]	Jena P, Sun Q 2018 Chem. Rev. 118 5755 Google Scholar
[40]	Chacko S, Joshi K, Kanhere D G 2004 Phys. Rev. Lett. 92 135506 Google Scholar
[41]	Rapacioli M, Tarrat N, Spiegelman F 2018 J. Phys. Chem. A 122 4092 Google Scholar
[42]	Breaux G.A, Benirschke R C, Sugai T, Kinnear B S, Jarrold M F 2003 Phys. Rev. Lett. 91 215508 Google Scholar
[43]	Nelli D, Ferrando R 2019 Nanoscale 11 13040 Google Scholar
[44]	Settem M, Kanjarla A K 2020 Sci. Rep. 10 3296 Google Scholar
[45]	Wu J, Qi L, You H, Gross A, Li J, Yang H 2012 J. Am. Chem. Soc. 134 11880 Google Scholar
[46]	Zhang J, Chen J, Hu P, Wang H 2020 Chin. Chem. Lett. 31 890 Google Scholar
[47]	Mao H K, Chen X J, Ding Y, Li B, Wang B 2018 Rev. Mod. Phys. 90 015007 Google Scholar
[48]	Plimpton S 1995 J. Comput. Phys. 117 1 Google Scholar
[49]	Cai J, Ye Y 1996 Phys. Rev. B 54 8398 Google Scholar
[50]	Stukowski A 2010 Model. Simul. Mater. Sc. 18 015012 Google Scholar
[51]	Larsen P M, Schmidt S, Schiøtz J 2016 Model. Simul. Mater. Sc. 24 055007 Google Scholar
[52]	蒋元祺 2015 博士学位论文(长沙: 湖南大学) Jiang Y Q 2015 Ph. D. Dissertation (Changsha: Hunan University) (in Chinese)

[1]	邢赫威, 陈占秀, 杨历, 苏瑶, 李源华, 呼和仓. 不凝性气体对纳米通道内水分子流动传热影响的分子动力学模拟. 物理学报, 2024, 73(9): 094701. doi: 10.7498/aps.73.20240192
[2]	文大东, 祁青华, 黄欣欣, 易洲, 邓永和, 田泽安, 彭平. 液态Ta快凝过程中团簇的遗传及其与局域对称性的关联. 物理学报, 2023, 72(24): 246101. doi: 10.7498/aps.72.20231153
[3]	华颖鑫, 陈小辉, 李俊, 郝龙, 孙毅, 王玉峰, 耿华运. 钒的冲击熔化原位X射线衍射测量研究. 物理学报, 2022, 71(7): 076201. doi: 10.7498/aps.71.20212065
[4]	郑治秀, 张林. Fe基体中包含Cu团簇的Fe-Cu二元体系在升温过程中结构变化的原子尺度计算. 物理学报, 2017, 66(8): 086301. doi: 10.7498/aps.66.086301
[5]	汪志刚, 黄娆, 文玉华. Pt-Au核-壳结构纳米粒子热稳定性的分子动力学研究. 物理学报, 2013, 62(12): 126101. doi: 10.7498/aps.62.126101
[6]	周耐根, 胡秋发, 许文祥, 李克, 周浪. 硅熔化特性的分子动力学模拟–-不同势函数的对比研究. 物理学报, 2013, 62(14): 146401. doi: 10.7498/aps.62.146401
[7]	李春丽, 段海明, 买力坦, 开来木. Aln(n=13–32)团簇熔化行为的分子动力学模拟研究. 物理学报, 2013, 62(19): 193104. doi: 10.7498/aps.62.193104
[8]	汪志刚, 黄娆, 文玉华. Au-Pd共晶纳米粒子熔化行为的分子动力学研究. 物理学报, 2012, 61(16): 166102. doi: 10.7498/aps.61.166102
[9]	张英杰, 肖绪洋, 李永强, 颜云辉. 分子动力学模拟Cu(010)基体对负载Co-Cu双金属团簇熔化过程的影响. 物理学报, 2012, 61(9): 093602. doi: 10.7498/aps.61.093602
[10]	韩小静, 王音, 林正喆, 张文献, 庄军, 宁西京. 团簇异构体生长概率的理论预测. 物理学报, 2010, 59(5): 3445-3449. doi: 10.7498/aps.59.3445
[11]	卢敏, 许卫兵, 刘维清, 侯春菊, 刘志勇. 银纳米杆高温熔化断裂弛豫性能的原子级模拟研究. 物理学报, 2010, 59(9): 6377-6383. doi: 10.7498/aps.59.6377
[12]	樊沁娜, 李蔚, 张林. 熔融Cu57团簇在急冷过程中弛豫和局域结构转变的分子动力学研究. 物理学报, 2010, 59(4): 2428-2433. doi: 10.7498/aps.59.2428
[13]	张林, 徐送宁, 李蔚, 孙海霞, 张彩碚. 小尺寸铜团簇冷却与并合过程中结构变化的原子尺度研究. 物理学报, 2009, 58(13): 58-S66. doi: 10.7498/aps.58.58
[14]	文玉华, 孙世刚, 张杨, 朱梓忠. 铂纳米晶在升温过程中结构演化与熔化特征的原子级模拟研究. 物理学报, 2009, 58(4): 2585-2589. doi: 10.7498/aps.58.2585
[15]	赵骞, 张林, 祁阳, 张宗宁. 低温下Cu13团簇负载于Cu(001)表面上结构变化的分子动力学研究. 物理学报, 2009, 58(13): 47-S52. doi: 10.7498/aps.58.47
[16]	张林, 张彩碚, 祁阳. 低温下Au959团簇负载于MgO(100)表面后结构变化的分子动力学研究. 物理学报, 2009, 58(13): 53-S57. doi: 10.7498/aps.58.53
[17]	张宗宁, 刘美林, 李蔚, 耿长建, 赵骞, 张林. 熔融Cu55团簇在Cu(010)表面上凝固过程的分子动力学模拟. 物理学报, 2009, 58(13): 67-S71. doi: 10.7498/aps.58.67
[18]	徐送宁, 张林, 张彩碚, 祁阳. 熔融Cu55团簇在铜块体中凝固过程的分子动力学模拟. 物理学报, 2009, 58(13): 40-S46. doi: 10.7498/aps.58.40
[19]	刘建廷, 段海明. 不同势下铱团簇结构和熔化行为的分子动力学模拟. 物理学报, 2009, 58(7): 4826-4834. doi: 10.7498/aps.58.4826
[20]	张林, 王绍青, 叶恒强. 大角度Cu晶界在升温、急冷条件下晶界结构的分子动力学研究. 物理学报, 2004, 53(8): 2497-2502. doi: 10.7498/aps.53.2497

期刊类型引用(1)
1. 王海鹏，王庆. 表征液态金属结构双体分布函数的讲授方法. 物理与工程. 2022(01): 155-161 . 百度学术
其他类型引用(0)

计量

文章访问数: 8954
PDF下载量: 143
被引次数: 1

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

搜索

留言板

难熔金属钒熔化行为的局域原子结构模拟与分析