Search

Article

x

留言板

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

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

Phase transitions of palladium under dynamic shock compression

Liu Ze-Tao Chen Bo Ling Wei-Dong Bao Nan-Yun Kang Dong-Dong Dai Jia-Yu

Liu Ze-Tao, Chen Bo, Ling Wei-Dong, Bao Nan-Yun, Kang Dong-Dong, Dai Jia-Yu. Phase transitions of palladium under dynamic shock compression. Acta Phys. Sin., 2022, 71(3): 037102. doi: 10.7498/aps.71.20211511
Citation: Liu Ze-Tao, Chen Bo, Ling Wei-Dong, Bao Nan-Yun, Kang Dong-Dong, Dai Jia-Yu. Phase transitions of palladium under dynamic shock compression. Acta Phys. Sin., 2022, 71(3): 037102. doi: 10.7498/aps.71.20211511

Phase transitions of palladium under dynamic shock compression

Liu Ze-Tao, Chen Bo, Ling Wei-Dong, Bao Nan-Yun, Kang Dong-Dong, Dai Jia-Yu
Article Text (iFLYTEK Translation)
PDF
HTML
Get Citation
  • For palladium (Pd) as a typical high-pressure standard material, studying its structural changes and thermodynamic properties under extreme conditions is widely demanded and challenging. Particularly, the solid-solid phase transition process of Pd under shock loading is understood still scarcely. In this paper, using the classical molecular dynamics simulations with embedded atom method (EAM) based on the interatomic potential, we investigate the phase transition of single crystal Pd from atomic scale under shock loading. A series of structural features is observed in a pressure range of 0–375 GPa, revealing that the structure feature transforms from the initial face-centered cubic (FCC) structure to the stacking faults body-centered cubic (BCC) structure with hexagonal close-packed (HCP) structure, and finally complete melting. Under shock loading of 100 oriented bulk Pd, we find the transformation to BCC structure can take place almost at 70.0 GPa, which is much lower than the previous static calculation result. In addition, we find that the phase transition depends on the direction initially impacting crystal. Under impacting along the 110 direction and the 111 direction, the FCC-BCC phase transition pressures increase to 135.8 GPa and 165.4 GPa, respectively. Also, the introduction of defects will increase the phase transition pressure of FCC-BCC by 20–30 GPa in comparison with perfect crystals, which is verified by the distribution of the potential energy. An interesting phenomenon that FCC-BCC transition pressure of Pd decreases under shock loading is found in this work, which provides a new theoretical insight into the application of high pressure experiments in the future.
      PACS:
      71.15.Pd(Molecular dynamics calculations (Car-Parrinello) and other numerical simulations)
      62.50.Ef(Shock wave effects in solids and liquids)
      61.50.Ks(Crystallographic aspects of phase transformations; pressure effects)
      71.20.Eh(Rare earth metals and alloys)
      Corresponding author: Chen Bo, chenbochain@nudt.edu.cn ; Dai Jia-Yu, jydai@nudt.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11774429, 12104507), the NSAF Joint Fund, China (Grant No. U1830206), the National Key R&D Program of China (Grant No. 2017YFA0403200), and the China Postdoctoral Science Foundation (Grant No. 2019M664024).

    极端条件下物质的结构及其相变动力学是行星物理、冲击波物理、凝聚态物理以及核物理中具有重要研究价值的课题之一. 近年来随着高压实验设备与技术的快速发展, 常压下无法观测的新物理、新现象和新规律在高压研究中逐渐被发现. 快速变化的外部载荷可以产生超高的温度压力环境[1,2], 从而形成与常压下不同的热力学状态[3-5]、独特的物理性能和亚稳态结构[6-10], 并且材料在高温高压下引发相变后的物理、力学性质相比初相有着本质的不同[11], 因此研究极端条件下材料的热力学性质及结构变化动力学过程对于科学研究和工业应用都具有重要的意义[12,13]. 钯是铂族金属中熔点最低, 密度最小的金属元素, 因其卓越的吸氢能力, 钯在电子信息制造业及医疗器械制备等领域[14-16]均有重要应用, 并且成为20世纪80年代末极具关注的冷聚变实验的关键成分[17]. 得益于优异的热力学性能与独特的物理特性, 钯吸引了大量的关注[18-21], 然而目前大多数研究集中于常态下, 对于极端条件下钯的性质则鲜有所见. 常压下钯为典型的面心立方(face-centered cubic, FCC)结构金属, 在高压下仍能维持稳定的FCC结构, 因而在高压实验中常常被用作压力校准器. 然而高温高压下钯的热力学性质以及结构变化过程仍不明晰, 只有少数实验和理论研究在温度-压强同时作用下的相变演化, 这使得研究极端条件下钯的相变过程具有重要意义.

    前人对钯在高温高压条件下的相变研究主要集中在熔化等固-液相变过程. Folies和Adams[22]使用由他们团队[23]早先开发的多体嵌入原子势函数(embedded atom method, EAM)计算金属钯的固相和液相的热力学性质, 通过蒙特卡罗模拟和准简谐近似(quasi harmonic approximation, QHA)方法开展固相和液相的自由能计算, 并利用固相和液相的吉布斯自由能曲线的交汇点来确定钯的熔点. 2011年刘中利等[24]通过基于EAM势的分子动力学模拟计算了钯的熔化曲线, 利用8000个原子构造固-液共存系, 计算得到了钯的熔点, 同时还使用25万个钯原子, 沿着100晶向利用非平衡分子动力学方法(nonequilibrium molecular dynamics, NEMD)模拟冲击波诱导熔化, 相比两相共存系的模拟结果在226 GPa下产生了18.3 %的过热熔化. 2013年Errandonea[25]利用激光加热金刚石顶砧(diamond-anvil cell, DAC)的实验方法首次测量了钯的高压熔化曲线, 与之前刘中利等[24]的理论计算结果有明显的分歧, 实验发现在30 GPa下与理论计算存在约为500 K的差距, 但与1999年Jeong和Chang[26]基于分子动力学的计算结果在20 GPa以上的高压范围内拟合良好. 2015年刘中利等[27]提出一种冲击熔化的模拟方法并证明能够精确有效地确定金属的熔化曲线, 仅用756个钯原子测量了熔化曲线, 且与两相法模拟结果符合, 在理论研究方面显著提升了计算效率, 但仍未解决理论与实验在熔化曲线上的差别. 对于钯的结构性质方面刘中利等[24]利用第一性原理对钯在高压下的相变过程做了详细的分析, 通过比较FCC结构、密排六方结构(hexagonal close packed, HCP) 和体心立方结构(body-centered cubic, BCC)在QHA下的吉布斯自由能, 发现在0—500 GPa的范围内不同温度条件下的FCC和HCP结构的自由能差值曲线没有交点, 得出FCC结构的钯原子在高达500 GPa和5000 K的条件下是结构稳定的. 至今为止, 关于钯在极端条件下的结构相变规律研究仍然匮乏, 并且对于钯相变的理论研究大都停留在准静态压缩. 而物质在冲击压缩这种快速载荷作用下, 其内部动态行为表现与静态[28]或准静态载荷作用下截然不同. 阐明冲击压缩下钯材料结构转变的微观动力学过程并探索其相变的物理机理, 不论是从理论上进一步理解钯的物理性质, 还是将其更好地应用于高温高压环境下, 都具有重要意义.

    实时观测冲击下材料结构的响应极具挑战性, 伴随着计算机水平的快速发展, 至今为止已经发展出多种计算材料学的模拟方法, 其中分子动力学模拟已成为在微纳米尺度上理解材料动态响应的重要研究方法[29-32], 分子动力学能够模拟材料中原子的实时运动轨迹并获得动态演变细节, 在常规实验手段难以达到的时间和空间尺度上对材料的结构特征和相变现象进行模拟研究, 对于理解材料在原子尺度上的冲击特性是极具优势的. 本文利用分子动力学模拟研究了钯在高温高压条件下的热学性质, 用数十万原子通过多尺度冲击压缩技术模拟了不同条件下的FCC理想晶体钯在冲击速度范围为4.2—9.0 km/s下的动态响应结构演化. 研究发现, 当沿着100晶向冲击时, 压强达到35.6 GPa时FCC结构的钯原子含量呈现下降趋势, 同时HCP结构在FCC结构中以层错的形式出现; 当压强达到约为57.5 GPa时BCC结构开始明显增加, HCP结构消失; 当压强约为70.0 GPa时, BCC的含量快速增加, 并在94.9 GPa时达到最大值, 这意味着FCC结构完成了向BCC结构的转变过程. 相较于前人静态压缩模拟的计算结果[23], 本文采用更为贴近实验的动态压缩结果, 揭示了理想晶体钯在1382.3 K和94.9 GPa时完成了FCC-BCC的结构转变现象. 此外, 当沿着110晶向进行冲击压缩时, 在冲击压强约为135.8 GPa时发生了FCC-BCC的相变过程, 而沿着111晶向时对应相变压强为165.4 GPa, 显示出强烈的晶向依赖性. 另外在初始的晶体结构中引入缺陷后, 相较完美晶体的相变压强大约有20—30 GPa的增幅. 本文揭示了单晶钯在动态冲击压缩下FCC-BCC相变规律, 通过模拟原子尺度层面的动力学过程为实验研究提供了理论支撑, 更加深入地理解钯的热力学性质及结构特征, 从而达到拓展应用领域、挖掘其潜在特性的目的.

    冲击波的模拟采用LAMMPS (large-scale atomic/molecular massively parallel simulator)[33]程序包里的多尺度冲击技术(multi-scale shock technique, MSST)[34]模块. 系统中所有原子都遵循修正的拉格朗日关系来更新位置和速度, 以将整个体系约束在冲击雨贡纽条件下. 当受到冲击速度为us的冲击加载时, 物质遵循以下的雨贡纽关系:

    up=us(1ρ0/ρ),
    (1)
    pp0=u2sρ0(1ρ0/ρ),
    (2)
    ee0=p0(1ρ01ρ)+u2s2(1ρ0ρ)2,
    (3)

    这里up, ρ, pe分别是粒子速度、密度、压强和单位质量的能量, 下标为0表示冲击波前的物理量. 将守恒方程(1)—(3)作为额外的约束加入原子的拉格朗日量中, 使整个体系中的原子在达到平衡时处于冲击速度为us的冲击波扫过后的准稳态. 其拉格朗日量可以表示为

    L=T({˙ri})V({ri})+12Q˙υ2+12u2sυ20(υ0υ)2+p0(υ0υ),
    (4)
    Q¨υ=TυVυp0u2sυ20(υ0υ),
    (5)

    这里T, V, Qυ分别是单位质量的动能、势能、类质量参数以及比体积. 通过调整盒子的单轴长度来实现应力的加载, 并在雨贡纽方程的约束下进行演化模拟冲击加载的过程, 直到系统达到平衡(˙υ = 0), 整个体系处于波后的准稳态. 相较于动量镜等基于大体系的非平衡方法, MSST方法通过在运动方程中额外加入雨贡纽方程作为约束, 能将体系约束在冲击状态附近, 因而本工作中所统计的热力学状态都是体系处于冲击波扫过后的准稳态, 并且经过长时间的平衡弛豫过程以达到体系内部的应力与温度均在空间保持一致. 在模拟中通过设置冲击波的速度来体现体系不同的冲击状态, 每一个冲击速度对应于冲击波传播后达到准稳态的压力值和温度值, P-T-V的状态方程变化符合雨贡纽关系. 具体的描述可参考Reed 等[34]的文章.

    本文的MSST计算中, 使用Foiles等[23]开发的多体嵌入原子势对钯原子间的相互作用进行描述. 通过调整模拟盒子单一方向的长度变化来对体系施加单轴方向的应力, 冲击压缩的方向确定为z轴方向, 钯金属体系的初始结构为FCC结构, 计算晶胞尺寸为30 × 30 × 80, 晶格常数为3.8907 Å, 由288000个原子组成, 对3个方向均采用周期性边界条件. 冲击加载前先使用NPT系综在300 K和1 bar (1 bar = 105 Pa)下对整个体系进行50 ps的弛豫. 平衡后对体系加载冲击速度为4.2—9.0 km/s的冲击波, 间隔为0.2 km/s. 设置类质量参数为3600, 人工黏度系数为0.0903, 转换因子为0.01以使得系统更快达到波后的准稳态. 整个MSST模拟持续500 ps, 待体系充分平衡后统计轨迹及热力学数据. 在实现冲击波沿不同晶向进行传递时, 通过调整沿冲击波方向的晶格基矢来确定晶向, 沿不同晶向体系加载冲击之后, 待冲击波扫过之后的体系达到准稳态再获取波后结构信息; 模拟构建含缺陷晶体类型时, 通过随机删除指定占比的原子来设置初始晶格体系的缺陷率来模拟点缺陷的情况, 用来进行本次工作中的对比计算.

    为对冲击过程中的结构进行分析, 采取了两种结构分析方法, 借助开源可视化工具OVITO (open visualization tool)[35]实现. 其一为径向分布函数(radial distribution function, RDF) g(r), RDF定义为以一个粒子为中心, 在半径rr + dr的空间范围内发现另一个粒子的概率, 用以表征结构的无序程度. 其定义为

    g(r)=VN2ijiδ(rrij)
    (6)

    其中V为系统的体积, N为系统中的原子数, riji原子到j原子的位矢. RDF已被用作揭示材料结构相变的有效方法[36-38], 曲线中的峰表示相邻原子相互靠近的概率. RDF曲线中第二峰和第三峰出现融合趋势并形成“短波浪”形状, 该特征的出现可被看作是FCC结构向BCC结构转变的一个判据[39].

    另一种结构分析方法是由Honeycutt和Andersen[40]提出的公共近邻分析(common neighbor analysis, CNA)方法, 该方法已经广泛地运用在固-固及固-液相变等结构分析当中, 并被证实是行之有效的方法[36,41]. 该方法是根据原子周围相邻原子键的拓扑结构进行分类, 具体可以通过区分近邻原子的配位数来分辨结构类型, 例如对于BCC结构其第一近邻配位数为8, 而FCC结构和HCP结构的第一近邻配位数为12. 近邻配位数可以通过截断半径r处的RDF积分得到:

    Nr=4π ρNr0r2g(r)dr
    (7)

    其中ρN为配位原子数密度, 即粒子个数与体系空间体积之比. 对于传统CNA方法, 往往通过选取全局统一的截断半径对所有原子晶体结构进行表征. 但对于本次工作中涉及到的多相系统, 以及考虑到高温或者原子应力分布不均引起的位置波动, 传统CNA方法中全局固定的截断半径难以描述不同环境下的原子结构. 因此使用自适应公共近邻分析(adaptive CNA, a-CNA)的方法[42]对处于不同原子环境中的粒子分别计算截断半径, 该方法是传统CNA的扩展, 为多相系统的结构分析提供更加精准的判断方法. 文献[43, 44]详细介绍了这种方法.

    在分子动力学模拟中, 势函数的准确性对于模拟结果的合理性至关重要[45]. 为了检验所选势函数的准确性和模拟方法的有效性, 首先计算了钯的状态方程(equation of state, EOS). 图1给出了冲击波在沿100晶向传递过程中体积比随压强的变化关系, 同时给出了以往的冲击压缩实验以及理论计算[24,46-49]中钯的EOS关系, 结果表明在0—375 GPa的范围内使用EAM势函数能得到与实验数据一致的EOS预测结果. 对不同大小体系的EOS进行了收敛性测试, 分别选用24000 (图中蓝色方形线条)和288000 (图中红色圆形线条)个原子的体系在相同条件下进行模拟, 两支线条基本符合, 因此证明在EOS上24000个原子体系已达到收敛标准.

    图 1 理论计算得到的雨贡纽P-V曲线与冲击实验数据的比较\r\nFig. 1. Hugoniot P-V curves compared with the shock experimental data and theoretical results.
    图 1  理论计算得到的雨贡纽P-V曲线与冲击实验数据的比较
    Fig. 1.  Hugoniot P-V curves compared with the shock experimental data and theoretical results.

    利用a-CNA方法分析了不同冲击速度加载下体系的结构变化. 不同结构原子的占有率随冲击压力的变化如图2(a)所示, 随着冲击速度的逐渐增大, 压强大于35 GPa时, FCC结构的钯原子在整个体系中所占比例开始下降, HCP结构的原子开始显现, 在图2(b)中观察到HCP结构的钯原子在FCC结构中以层错的形式存在. 当冲击压强增加到57.5 GPa时, BCC结构开始增加, HCP结构随之消失, 同样从图2(b)可以看出, BCC结构也是以层错的形式出现在FCC结构中, 没有出现BCC的整体块状结构. 在70.0—94.9 GPa的压强范围中BCC结构的含量快速增长直到其含量达到最大值, 此时FCC和HCP的原子比例几乎降至0, 意味着FCC结构向BCC结构转换完成. 当压强超过232.0 GPa时, 晶体的结构特征全部消失, 与图3(b)所示的P-T图中对应的间断点(228.3 GPa, 5488.8 K)相对应, 表明此刻出现了熔化, 其微观表现为相对体积和原子平均能量都随着温度升高而线性增大, 达到某一临界温度时体积骤降, 原子平均能量骤升, 之后两者又继续保持线性增大. 该压强间断点与刘中利等[24]100晶向的冲击模拟获得的数据(压强和温度于5855 K和226 GPa时出现熔化)一致. 总之, 本文的结果准确地揭示了沿着100晶向冲击压缩过程中的相变过程及热力学响应. 与小体系静态计算[24]中钯FCC-BCC相变压力大于500 GPa的结果不同, 我们发现沿100晶向加载冲击波后, 在较低的相变压力70.0 GPa时便能出现FCC-BCC的结构转变. 可见在冲击压缩这种载荷快速加载的过程中, 温度与压力等热力学状态的耦合变化, 以及原子剧烈位移的动态过程会降低钯的固-固相变压力.

    图 2 (a)沿$\left\langle {100} \right\rangle $晶向冲击后不同结构特征的钯的原子分数与纵向应力的关系; (b)不同结构特征随纵向应力变化的可视化表示; (c)温度与纵向应力的关系, 其中红色点线图为刘中利等[24]冲击熔化的分子动力学模拟数据\r\nFig. 2. (a) Fraction of atoms with different structural features versus longitudinal stress along the crystallographic direction $\left\langle {100} \right\rangle $; (b) visual representation of the variation of different structural features with longitudinal stress; (c) temperature versus longitudinal stress, the red dotted line plot shows the MD simulation data of Liu et al.[24] shock melting.
    图 2  (a)沿100晶向冲击后不同结构特征的钯的原子分数与纵向应力的关系; (b)不同结构特征随纵向应力变化的可视化表示; (c)温度与纵向应力的关系, 其中红色点线图为刘中利等[24]冲击熔化的分子动力学模拟数据
    Fig. 2.  (a) Fraction of atoms with different structural features versus longitudinal stress along the crystallographic direction 100; (b) visual representation of the variation of different structural features with longitudinal stress; (c) temperature versus longitudinal stress, the red dotted line plot shows the MD simulation data of Liu et al.[24] shock melting.
    图 3 沿不同晶向冲击后的结果 (a) $\left\langle {100} \right\rangle $晶向; (b) $\left\langle {110} \right\rangle $晶向; (c) $\left\langle {100} \right\rangle $晶向. 其中实线表示冲击速度在5.0—7.8 km/s范围内结构的RDF图像, 点划线表示与之对应的配位数曲线, 黑色虚线位置指向初始状态下的第一峰的位置\r\nFig. 3. Results along different crystallographic directions: (a) $\left\langle {100} \right\rangle $; (b) $\left\langle {110} \right\rangle $; (c) $\left\langle {111} \right\rangle $. The solid lines indicate radial distribution functions of deformed micro-structure shocked at a range of shock velocity of 5.0–7.8 km/s. The dotted lines indicate the coordination number curves that correspond to the RDF, and the dashed lines points to the position of the first peak in the initial state.
    图 3  沿不同晶向冲击后的结果 (a) 100晶向; (b) 110晶向; (c) 100晶向. 其中实线表示冲击速度在5.0—7.8 km/s范围内结构的RDF图像, 点划线表示与之对应的配位数曲线, 黑色虚线位置指向初始状态下的第一峰的位置
    Fig. 3.  Results along different crystallographic directions: (a) 100; (b) 110; (c) 111. The solid lines indicate radial distribution functions of deformed micro-structure shocked at a range of shock velocity of 5.0–7.8 km/s. The dotted lines indicate the coordination number curves that correspond to the RDF, and the dashed lines points to the position of the first peak in the initial state.

    为了探究钯在冲击压缩过程中影响其FCC-BCC相变压强的因素, 首先模拟了冲击波沿不同晶向进行传播的过程. 并通过RDF和a-CNA方法对冲击压缩后的结构进行了分析. 图3分别给出了沿100, 110111三个晶向加载冲击的RDF以及配位数曲线. 对于100向, 冲击速度为5.2—7.2 km/s时可以观察到RDF的第二峰与第三峰的融合趋势, 呈现为“短波状”[39]的特征, 在r = 4.25 Å 处(对应于RDF的波谷)附近的配位数呈现出“平台”状特征; 而当相变发生后, 该结构特征消失并且配位数的变化曲线趋于平滑, 表明此径向处的原子配位发生改变, 说明发生了相变, 对应的温度和压强范围在1382—1699 K和95—108 GPa, 这与a-CNA的判断基本一致. 另外, 同样的波状特征在110111起始出现的冲击速度分别为6.6及7.2 km/s, 分别对应冲击压力135.8和182.1 GPa, 意味着相比111晶向, 沿110晶向冲击会在较低压力下发生FCC-BCC的相变, 并且均比100晶向需要更高的冲击压力. 除此之外, 随着冲击速度的增加, 在沿着各个晶向冲击下, 第一峰出现左移(图3中虚线位置), 并有逐渐拓宽的趋势, 峰值高度减小以及整体的展宽, 在低速冲击(5.0 km/s < us < 5.5 km/s)区域, 结合图2(a)中的a-CNA诊断方法及原子结构图发现, 沿100晶向冲击会使体系出现带HCP的层错结构, 沿110晶向会出现FCC-BCC-HCP多相共存的纳米结构, 这些微观结构的出现会导致冲击后RDF峰展宽; 在高速冲击(us > 7.2 km/s)区域, RDF峰的展宽与强度降低是由于熔化或接近熔化导致的非谐效应造成的, 表明随压力升高体系的非谐效应增大, 表现出向液相转变的趋势. 为进一步分析沿不同晶向的冲击结构, 利用a-CNA方法对不同冲击方向的波后体系进行了对比. 图4给出了沿110111晶向在冲击压缩后的结构随冲击压强的变化情况, 将a-CNA采样点中BCC结构增长首次趋于稳定的压力值作为FCC-BCC结构的相变点, 综合图2(a)中沿100晶向冲击的结果, 沿100, 110111这3个晶向加载冲击对应的相变压力依次为70.0, 135.8和165.4 GPa, 这与上述的RDF结果表现基本一致, 相比其他两个晶向, 沿111晶向冲击更不容易发生相的转变.

    图 4 利用a-CNA方法得到的沿不同晶向冲击的不同结构比例随压强的变化 (a) $\left\langle {110} \right\rangle $晶向; (b) $\left\langle {111} \right\rangle $晶向\r\nFig. 4. Fraction of atoms with different structural versus compressing stress along different crystallographic directions obtained by a-CNA: (a) $\left\langle {110} \right\rangle $; (b) $\left\langle {111} \right\rangle $.
    图 4  利用a-CNA方法得到的沿不同晶向冲击的不同结构比例随压强的变化 (a) 110晶向; (b) 111晶向
    Fig. 4.  Fraction of atoms with different structural versus compressing stress along different crystallographic directions obtained by a-CNA: (a) 110; (b) 111.

    另外, 由于在实际材料生长及应用过程中, 很有可能会引入一些缺陷, 因此还探讨了包含不同缺陷率的初始结构对钯在冲击压缩下FCC-BCC相变压力的影响, 结果如图5(a)所示. 可以看出, 缺陷的引入给状态方程带来了影响, 在同等压力下相比完美晶体, 含缺陷的体系会获得更高温度, 在3%缺陷率时达到最大值. 同样, 通过a-CNA跟踪了不同缺陷率体系沿100晶向在冲击后结构特征随着压力的变化情况, 如图5所示. 相比无缺陷晶体(图2(a))可以发现, 随着缺陷的引入相变压力整体呈现增大趋势, 但与缺陷的占比并不遵循正相关关系, 相变压强在3%缺陷率时达到最大值99.1 GPa, 在5%缺陷率时又下降至84.8 GPa. 为进一步分析引入缺陷增大相变压力的原因, 分别计算了完美晶体、3%缺陷率以及5%缺陷率的缺陷晶体在冲击前原始状态及冲击后发生相变时的原子势能, 冲击速度均设置为6 km/s, 图6给出了冲击前后的势能差(ΔEp)在体系中的分布情况, 用以表征体系中FCC转变为BCC结构所需要克服的能垒. 通过图6(a)的颜色分布以及图6(b)的概率分布图可看出, 无缺陷晶体在冲击后的势能相比初始状态是均匀增加的, 引入3%和5%的点缺陷后势能增幅变大且分布不均匀, 因此相比完美晶体, 含有缺陷的体系需要克服更大的能垒用于完成结构相变, 而在缺陷为5%时所需要克服的能量有明显减少, 这与本文a-CNA结果一致, 因此引入点缺陷会增加钯在冲击压缩下FCC-BCC的转变压力.

    图 5 (a)不同缺陷率的晶体结构在冲击压缩后温度与压力的关系; (b), (c), (d)分别表示冲击压缩下含不同缺陷率体系的原子类别占比随冲击压强的变化\r\nFig. 5. (a) Comparison of Hugoniot P-T curve among crystal structures with different porosities defects; (b), (c), (d) comparison of phase transition processes of atoms of system with variation porosities defects under shock compression.
    图 5  (a)不同缺陷率的晶体结构在冲击压缩后温度与压力的关系; (b), (c), (d)分别表示冲击压缩下含不同缺陷率体系的原子类别占比随冲击压强的变化
    Fig. 5.  (a) Comparison of Hugoniot P-T curve among crystal structures with different porosities defects; (b), (c), (d) comparison of phase transition processes of atoms of system with variation porosities defects under shock compression.
    图 6 (a)完美晶体与含缺陷晶体在冲击前后的势能差值分布图; (b)不同体系中势能差值的概率分布曲线\r\nFig. 6. (a) Distribution diagram of potential energy difference between perfect crystal and defective crystal before and after compression; (b) probability distribution curves of potential energy differences in perfect crystal and defective crystal.
    图 6  (a)完美晶体与含缺陷晶体在冲击前后的势能差值分布图; (b)不同体系中势能差值的概率分布曲线
    Fig. 6.  (a) Distribution diagram of potential energy difference between perfect crystal and defective crystal before and after compression; (b) probability distribution curves of potential energy differences in perfect crystal and defective crystal.

    基于分子动力学的方法, 模拟了冲击压缩下钯的结构相变过程, 并探讨了沿不同晶向冲击及不同缺陷率的初始构型对相变过程的影响. 首先得到了与实验一致的雨贡纽状态方程, 并结合RDF图谱的特征峰和配位数曲线以及a-CNA的结构分析手段, 在0—375 GPa的压力区间观察到了一系列结构变化特征, 从初始的FCC结构至出现HCP层错, 到FCC转换成BCC结构, 最后直至熔化. 在沿100晶向的冲击加载下, 发现了低至70.0 GPa的FCC-BCC相变压力, 本文的状态方程与以往的动态实验及NEMD的分子动力学模拟结果拟合良好, 低于之前静态计算中所报道的结果, 由此表明了冲击压缩过程会降低钯的固-固相变压力. 此外还发现, 沿其他晶向加载冲击波会增大FCC-BCC的相变压力, 沿110111方向体系的FCC-BCC相变压力能分别增大至135.8和165.4 GPa. 最后, 在初始结构中引入不同缺陷率的点缺陷, 发现对体系的状态方程产生了影响, 相变压力会明显提升20—30 GPa, 通过比较冲击前后体系的势能差值情况可以看出, 含有缺陷的结构在发生相变的过程中需要克服更大的能垒, 因此点缺陷的引入会增加FCC-BCC的转变压力. 通过此次工作系统地研究了钯在高压动态载荷加载下的热力学性质与相变过程, 发现了冲击压缩加载下降低固-固相变压力的特殊现象, 不仅为高速冲击下的相变过程提供参考方向, 也能为钯的高压实验等极端条件下的应用提供参考和理论依据.

    [1]

    Graziani F, Desjarlais M P, Redmer R, Trickey S B 2014 Frontiers and Challenges in Warm Dense Matter (Vol. 1) (Switzerland: Springer) pp1–3

    [2]

    Moses E I, Boyd R N, Remington B A, Keane C J, Al-Ayat R 2009 Phys. Plasmas 16 041006Google Scholar

    [3]

    Mitchell A C, Nellis W J, Moriarty J A, Heinle R A, Holmes N C, Tipton R E, Repp G W 1991 J. Appl. Phys. 69 2981Google Scholar

    [4]

    Nagao H, Nakamura K J, Kondo K, Ozaki N, Takamatsu K, Ono T, Shiota T, Ichinose D, Tanaka K A, Wakabayashi K 2006 Phys. Plasmas 13 052705Google Scholar

    [5]

    Yokoo M, Kawai N, Nakamura K G, Kondo K I 2008 Appl. Phys. Lett. 92 051901Google Scholar

    [6]

    Nellis W J, Weir S T, Mitchell A C 1996 Science 273 936Google Scholar

    [7]

    Weir S T, Mitchell A C, Nellis W J 1996 Phys. Rev. Lett. 76 1860Google Scholar

    [8]

    Hartley N J, Brown S, Cowan T E, Cunningham E, Döppner T, Falcone R W, Fletcher L B, Frydrych S, Galtier E, Gamboa E J 2019 Sci. Rep. 9 1

    [9]

    Briggs R, Coppari F, Gorman M G, Smith R F, Tracy S J, Coleman A L, Fernandez-Panella A, Millot M, Eggert J H, Fratanduono D E 2019 Phys. Rev. Lett. 123 045701Google Scholar

    [10]

    Sharma S M, Turneaure S J, Winey J M, Li Y, Rigg P, Schuman A, Sinclair N, Toyoda Y, Wang X, Weir N, Zhang J, Gupta Y M 2019 Phys. Rev. Lett. 123 045702Google Scholar

    [11]

    Li Q J, Liu B B 2016 Chin. Phys. B 25 076107Google Scholar

    [12]

    Guo Y L, Wei J H, Liu X, Ke X Z, Jiao Z Y 2021 Chin. Phys. B 30 016101Google Scholar

    [13]

    Lu M Y, Huang Y P, Tian F B, Li D, Duan D F, Zhou Q, Cui T 2020 Chin. Phys. B 29 053104Google Scholar

    [14]

    Hesse R W 2007 Jewelry-Making Through History: an Encyclopedia (Westport: Greenwood Publishing Group) p146

    [15]

    Colon P, Pradelle-Plasse N, Galland J 2003 Dent. Mater. 19 232Google Scholar

    [16]

    Harper C A 1997 Passive Electronic Component Handbook. McGraw-Hill Professional (New York: McGraw-Hill) pp580

    [17]

    Fleischmann M, Pons S, Hawkins M 1989 J. Electroanal. Chem. 261 301Google Scholar

    [18]

    Mutschele T, Kirchheim R 1987 Scr. Mater. 21 1351101

    [19]

    Pundt A, Sachs C, Winter M, Reetz M T, Fritsch D, Kirchheim R 1999 J. Alloys Compd. 293-295 480

    [20]

    Bonivardi A L, Baltanas M A 1994 React. Kinet. Catal. Lett. 52 95

    [21]

    Eastman J A, Thompson L J, Kestel B J 1993 Phys. Rev. 48 84Google Scholar

    [22]

    Foiles S M, Adams J B 1989 Phys. Rev. B 40 9

    [23]

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

    [24]

    Liu Z L, Yang J H, Cai L C, Jing F Q, Alfѐ D 2011 Phys. Rev. B 83 144113Google Scholar

    [25]

    Errandonea D 2013 Phys. Rev. B 87 054108Google Scholar

    [26]

    Jeong J W, Chang K 1999 J. Phys. Condens. Matter 11 3799Google Scholar

    [27]

    Liu Z L, Zhang X L, Cai L C 2015 J. Chem. Phys. 11 637

    [28]

    Zheng Z Y, Zhao J J 2016 Chin. Phys. B 25 076202Google Scholar

    [29]

    孙小伟, 褚衍东, 刘子江, 刘玉孝, 王成伟, 刘维民 2005 物理学报 54 5830Google Scholar

    Sun X W, Zhe Y D, Liu Z J, Liu Y X, Wang C W, Liu W M 2005 Acta Phys. Sin. 54 5830Google Scholar

    [30]

    段芳莉, 王家序, 雒建斌, 温诗铸 2007 物理学报 56 6552Google Scholar

    Duan F L, Wang J X, Luo J B, Weng S Z 2007 Acta Phys. Sin. 56 6552Google Scholar

    [31]

    邵建立, 王裴, 秦承森, 周洪强 2007 物理学报 56 5389Google Scholar

    Shao J L, Wang F, Qing C S, Zhou H Q 2007 Acta Phys. Sin. 56 5389Google Scholar

    [32]

    周化光, 林鑫, 王猛, 黄卫东 2013 物理学报 62 056803Google Scholar

    Zhou H G, Lin X, Wang M, Huang W D 2013 Acta Phys. Sin. 62 056803Google Scholar

    [33]

    Plimpton S J 1995 J. Comput. Phys. 117 1Google Scholar

    [34]

    Reed E J, Fried L E, Joannopoulos J D 2003 Phys. Rev. Lett. 90 235503Google Scholar

    [35]

    Stukowski A 2009 Modell. Simul. Mater. Sci. Eng. 18 015012

    [36]

    韦昭召, 马骁, 柯常波, 张新平 2020 物理学报 69 136102Google Scholar

    Wei Z Z, Ma X, Ke C B, Zhang X P 2020 Acta Phys. Sin. 69 136102Google Scholar

    [37]

    Kien P H 2014 Int. Scholarly Res. Not. 2014 253627

    [38]

    Sun H Y, Kang D D, Hou Y, Dai J Y 2017 Matter Radiat. at Extremes 2 287Google Scholar

    [39]

    Neogi A, Mitra N 2017 Sci. Rep. 7 1Google Scholar

    [40]

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

    [41]

    Schiøtz J, Tolla F D D, Jacobsen K W 1998 Nature 391 561Google Scholar

    [42]

    Stukowski A 2012 Model. Simul. Mater. Sci. Eng. 20 045021Google Scholar

    [43]

    Faken D, Jónsson H 1994 Comput. Mater. Sci. 2 279Google Scholar

    [44]

    Tsuzuki H, Branicio P, Rino J 2007 Comput. Phys. Commun. 177 518Google Scholar

    [45]

    Zeng Q Y, Dai J Y 2020 Sci. China Phys. Mech. Astron. 63 263011Google Scholar

    [46]

    Walsh J M, Rice M H, Mcqueen R G, Yarger F L 1957 Phys. Rev. 108 196Google Scholar

    [47]

    Marsh S P 1979 Los Alamos Shock Hugoniot Data (Berkeley: University of California Press)

    [48]

    Thiel M V 1997 Compendium of Shock Wave Data (Livermore: Lawrence Livermore National Laboratory) Report UCRL 50108

    [49]

    Mcqueen R G, Marsh, S P, Taylor J W, Fritz J N, Carter W J 1970 The Equation of State of Solids from Shock Wave Studies (New York: Academic) pp348–350

  • 图 1  理论计算得到的雨贡纽P-V曲线与冲击实验数据的比较

    Figure 1.  Hugoniot P-V curves compared with the shock experimental data and theoretical results.

    图 2  (a)沿100晶向冲击后不同结构特征的钯的原子分数与纵向应力的关系; (b)不同结构特征随纵向应力变化的可视化表示; (c)温度与纵向应力的关系, 其中红色点线图为刘中利等[24]冲击熔化的分子动力学模拟数据

    Figure 2.  (a) Fraction of atoms with different structural features versus longitudinal stress along the crystallographic direction 100; (b) visual representation of the variation of different structural features with longitudinal stress; (c) temperature versus longitudinal stress, the red dotted line plot shows the MD simulation data of Liu et al.[24] shock melting.

    图 3  沿不同晶向冲击后的结果 (a) 100晶向; (b) 110晶向; (c) 100晶向. 其中实线表示冲击速度在5.0—7.8 km/s范围内结构的RDF图像, 点划线表示与之对应的配位数曲线, 黑色虚线位置指向初始状态下的第一峰的位置

    Figure 3.  Results along different crystallographic directions: (a) 100; (b) 110; (c) 111. The solid lines indicate radial distribution functions of deformed micro-structure shocked at a range of shock velocity of 5.0–7.8 km/s. The dotted lines indicate the coordination number curves that correspond to the RDF, and the dashed lines points to the position of the first peak in the initial state.

    图 4  利用a-CNA方法得到的沿不同晶向冲击的不同结构比例随压强的变化 (a) 110晶向; (b) 111晶向

    Figure 4.  Fraction of atoms with different structural versus compressing stress along different crystallographic directions obtained by a-CNA: (a) 110; (b) 111.

    图 5  (a)不同缺陷率的晶体结构在冲击压缩后温度与压力的关系; (b), (c), (d)分别表示冲击压缩下含不同缺陷率体系的原子类别占比随冲击压强的变化

    Figure 5.  (a) Comparison of Hugoniot P-T curve among crystal structures with different porosities defects; (b), (c), (d) comparison of phase transition processes of atoms of system with variation porosities defects under shock compression.

    图 6  (a)完美晶体与含缺陷晶体在冲击前后的势能差值分布图; (b)不同体系中势能差值的概率分布曲线

    Figure 6.  (a) Distribution diagram of potential energy difference between perfect crystal and defective crystal before and after compression; (b) probability distribution curves of potential energy differences in perfect crystal and defective crystal.

  • [1]

    Graziani F, Desjarlais M P, Redmer R, Trickey S B 2014 Frontiers and Challenges in Warm Dense Matter (Vol. 1) (Switzerland: Springer) pp1–3

    [2]

    Moses E I, Boyd R N, Remington B A, Keane C J, Al-Ayat R 2009 Phys. Plasmas 16 041006Google Scholar

    [3]

    Mitchell A C, Nellis W J, Moriarty J A, Heinle R A, Holmes N C, Tipton R E, Repp G W 1991 J. Appl. Phys. 69 2981Google Scholar

    [4]

    Nagao H, Nakamura K J, Kondo K, Ozaki N, Takamatsu K, Ono T, Shiota T, Ichinose D, Tanaka K A, Wakabayashi K 2006 Phys. Plasmas 13 052705Google Scholar

    [5]

    Yokoo M, Kawai N, Nakamura K G, Kondo K I 2008 Appl. Phys. Lett. 92 051901Google Scholar

    [6]

    Nellis W J, Weir S T, Mitchell A C 1996 Science 273 936Google Scholar

    [7]

    Weir S T, Mitchell A C, Nellis W J 1996 Phys. Rev. Lett. 76 1860Google Scholar

    [8]

    Hartley N J, Brown S, Cowan T E, Cunningham E, Döppner T, Falcone R W, Fletcher L B, Frydrych S, Galtier E, Gamboa E J 2019 Sci. Rep. 9 1

    [9]

    Briggs R, Coppari F, Gorman M G, Smith R F, Tracy S J, Coleman A L, Fernandez-Panella A, Millot M, Eggert J H, Fratanduono D E 2019 Phys. Rev. Lett. 123 045701Google Scholar

    [10]

    Sharma S M, Turneaure S J, Winey J M, Li Y, Rigg P, Schuman A, Sinclair N, Toyoda Y, Wang X, Weir N, Zhang J, Gupta Y M 2019 Phys. Rev. Lett. 123 045702Google Scholar

    [11]

    Li Q J, Liu B B 2016 Chin. Phys. B 25 076107Google Scholar

    [12]

    Guo Y L, Wei J H, Liu X, Ke X Z, Jiao Z Y 2021 Chin. Phys. B 30 016101Google Scholar

    [13]

    Lu M Y, Huang Y P, Tian F B, Li D, Duan D F, Zhou Q, Cui T 2020 Chin. Phys. B 29 053104Google Scholar

    [14]

    Hesse R W 2007 Jewelry-Making Through History: an Encyclopedia (Westport: Greenwood Publishing Group) p146

    [15]

    Colon P, Pradelle-Plasse N, Galland J 2003 Dent. Mater. 19 232Google Scholar

    [16]

    Harper C A 1997 Passive Electronic Component Handbook. McGraw-Hill Professional (New York: McGraw-Hill) pp580

    [17]

    Fleischmann M, Pons S, Hawkins M 1989 J. Electroanal. Chem. 261 301Google Scholar

    [18]

    Mutschele T, Kirchheim R 1987 Scr. Mater. 21 1351101

    [19]

    Pundt A, Sachs C, Winter M, Reetz M T, Fritsch D, Kirchheim R 1999 J. Alloys Compd. 293-295 480

    [20]

    Bonivardi A L, Baltanas M A 1994 React. Kinet. Catal. Lett. 52 95

    [21]

    Eastman J A, Thompson L J, Kestel B J 1993 Phys. Rev. 48 84Google Scholar

    [22]

    Foiles S M, Adams J B 1989 Phys. Rev. B 40 9

    [23]

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

    [24]

    Liu Z L, Yang J H, Cai L C, Jing F Q, Alfѐ D 2011 Phys. Rev. B 83 144113Google Scholar

    [25]

    Errandonea D 2013 Phys. Rev. B 87 054108Google Scholar

    [26]

    Jeong J W, Chang K 1999 J. Phys. Condens. Matter 11 3799Google Scholar

    [27]

    Liu Z L, Zhang X L, Cai L C 2015 J. Chem. Phys. 11 637

    [28]

    Zheng Z Y, Zhao J J 2016 Chin. Phys. B 25 076202Google Scholar

    [29]

    孙小伟, 褚衍东, 刘子江, 刘玉孝, 王成伟, 刘维民 2005 物理学报 54 5830Google Scholar

    Sun X W, Zhe Y D, Liu Z J, Liu Y X, Wang C W, Liu W M 2005 Acta Phys. Sin. 54 5830Google Scholar

    [30]

    段芳莉, 王家序, 雒建斌, 温诗铸 2007 物理学报 56 6552Google Scholar

    Duan F L, Wang J X, Luo J B, Weng S Z 2007 Acta Phys. Sin. 56 6552Google Scholar

    [31]

    邵建立, 王裴, 秦承森, 周洪强 2007 物理学报 56 5389Google Scholar

    Shao J L, Wang F, Qing C S, Zhou H Q 2007 Acta Phys. Sin. 56 5389Google Scholar

    [32]

    周化光, 林鑫, 王猛, 黄卫东 2013 物理学报 62 056803Google Scholar

    Zhou H G, Lin X, Wang M, Huang W D 2013 Acta Phys. Sin. 62 056803Google Scholar

    [33]

    Plimpton S J 1995 J. Comput. Phys. 117 1Google Scholar

    [34]

    Reed E J, Fried L E, Joannopoulos J D 2003 Phys. Rev. Lett. 90 235503Google Scholar

    [35]

    Stukowski A 2009 Modell. Simul. Mater. Sci. Eng. 18 015012

    [36]

    韦昭召, 马骁, 柯常波, 张新平 2020 物理学报 69 136102Google Scholar

    Wei Z Z, Ma X, Ke C B, Zhang X P 2020 Acta Phys. Sin. 69 136102Google Scholar

    [37]

    Kien P H 2014 Int. Scholarly Res. Not. 2014 253627

    [38]

    Sun H Y, Kang D D, Hou Y, Dai J Y 2017 Matter Radiat. at Extremes 2 287Google Scholar

    [39]

    Neogi A, Mitra N 2017 Sci. Rep. 7 1Google Scholar

    [40]

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

    [41]

    Schiøtz J, Tolla F D D, Jacobsen K W 1998 Nature 391 561Google Scholar

    [42]

    Stukowski A 2012 Model. Simul. Mater. Sci. Eng. 20 045021Google Scholar

    [43]

    Faken D, Jónsson H 1994 Comput. Mater. Sci. 2 279Google Scholar

    [44]

    Tsuzuki H, Branicio P, Rino J 2007 Comput. Phys. Commun. 177 518Google Scholar

    [45]

    Zeng Q Y, Dai J Y 2020 Sci. China Phys. Mech. Astron. 63 263011Google Scholar

    [46]

    Walsh J M, Rice M H, Mcqueen R G, Yarger F L 1957 Phys. Rev. 108 196Google Scholar

    [47]

    Marsh S P 1979 Los Alamos Shock Hugoniot Data (Berkeley: University of California Press)

    [48]

    Thiel M V 1997 Compendium of Shock Wave Data (Livermore: Lawrence Livermore National Laboratory) Report UCRL 50108

    [49]

    Mcqueen R G, Marsh, S P, Taylor J W, Fritz J N, Carter W J 1970 The Equation of State of Solids from Shock Wave Studies (New York: Academic) pp348–350

  • [1] WEI Zhaozhao. Molecular dynamics simulation of bending behavior of B2-FeAl alloy nanowires with different crystallographic orientations. Acta Physica Sinica, 2025, 74(3): 036201. doi: 10.7498/aps.74.20241030
    [2] Phase Transitions of Palladium under Dynamic Shock Compression. Acta Physica Sinica, 2021, (): . doi: 10.7498/aps.70.20211511
    [3] Chen Chao, Duan Fang-Li. Effect of functional groups on crumpling behavior and structure of graphene oxide. Acta Physica Sinica, 2020, 69(19): 193102. doi: 10.7498/aps.69.20200651
    [4] Zhou Bian, Yang Liang. Molecular dynamics simulation of effect of cooling rate on the microstructures and deformation behaviors in metallic glasses. Acta Physica Sinica, 2020, 69(11): 116101. doi: 10.7498/aps.69.20191781
    [5] Sun Xiao-Wei, Song Ting, Liu Zi-Jiang, Wan Gui-Xin, Zhang Lei, Chang Wen-Li. Numerical prediction of structural stability and thermodynamic properties for MgF2 with fluorite- type structure under high pressure. Acta Physica Sinica, 2020, 69(15): 156202. doi: 10.7498/aps.69.20200289
    [6] Feng Xiao-Wei, Li Jun-Cheng, Wang Hong-Bo, Chang Jing-Zhen. Mesomechanism of elastic precursor decay in alumina under plate impact loading. Acta Physica Sinica, 2016, 65(16): 166201. doi: 10.7498/aps.65.166201
    [7] Hu Yong-Jin, Wu Yun-Pei, Liu Guo-Ying, Luo Shi-Jun, He Kai-Hua. Structural phase transition, electronic structures and optical properties of ZnTe. Acta Physica Sinica, 2015, 64(22): 227802. doi: 10.7498/aps.64.227802
    [8] Lin Jia-Qi, Li Xiao-Kang, Yang Wen-Long, Sun Hong-Guo, Xie Zhi-Bin, Xiu Han-jiang, Lei Qing-Quan. Molecular dynamics simulation study on the structure and mechanical properties of polyimide/KTa0.5Nb0.5O3 nanoparticle composites. Acta Physica Sinica, 2015, 64(12): 126202. doi: 10.7498/aps.64.126202
    [9] Zheng Xiao-Qing, Yang Yang, Sun De-Yan. Atomistic characterization of a modeled binary ordered alloy solid-liquid interface. Acta Physica Sinica, 2013, 62(1): 017101. doi: 10.7498/aps.62.017101
    [10] Dong Lei, Wang Wei-Guo. Molecular dynamics simulation on the struatural stability of [0 1 1] tilt incoherent 3 gain boundaries in pure copper. Acta Physica Sinica, 2013, 62(15): 156102. doi: 10.7498/aps.62.156102
    [11] Yan Xiao, Xin Zi-Hua, Zhang Jiao-Jiao. Molecular dynamics study on the structure and properties of silicon-graphdiyne. Acta Physica Sinica, 2013, 62(23): 238101. doi: 10.7498/aps.62.238101
    [12] Xu Chun-Long, Hou Zhao-Yang, Liu Rang-Su. Simulation study on thermodynamic, dynamic and structural transition mechanisms during the formation of Ca70Mg30 metallic glass. Acta Physica Sinica, 2012, 61(13): 136401. doi: 10.7498/aps.61.136401
    [13] Xia Dong, Wang Xin-Qiang. Structures and melting behaviors of ultrathin platinum nanowires. Acta Physica Sinica, 2012, 61(13): 130510. doi: 10.7498/aps.61.130510
    [14] Kai Hua, Li Yun-Chao, Guo De-Cheng, Li Shuang, Li Zhi-Jie. Molecular dynamics simulation of the structure characteristic of diamond-like carbon films influence by oblique incidence ion-beam-assisted deposition. Acta Physica Sinica, 2009, 58(7): 4888-4894. doi: 10.7498/aps.58.4888
    [15] Chen Yu-Xiang, Xie Guo-Feng, Ma Ying, Zhou Yi-Chun. Molecular-dynamics simulation of the structure and elastic constants of barium titanium. Acta Physica Sinica, 2009, 58(6): 4085-4089. doi: 10.7498/aps.58.4085
    [16] Cao Li-Xia, Shang Jia-Xiang, Zhang Yue. Molecular dynamics simulation of stress-induced martensitic phase transformation in NiAl. Acta Physica Sinica, 2009, 58(10): 7307-7312. doi: 10.7498/aps.58.7307
    [17] Ma Tian-Bao, Hu Yuan-Zhong, Wang Hui. Molecular dynamics simulation of the growth and structural properties of ultra-thin diamond-like carbon films. Acta Physica Sinica, 2006, 55(6): 2922-2927. doi: 10.7498/aps.55.2922
    [18] Sun Xiao-Wei, Chu Yan-Dong, Liu Zi-Jiang, Liu Yu-Xiao, Wang Cheng-Wei, Liu Wei-Min. Molecular dynamics study on the structural and thermodynamic properties of the zinc-blende phase of GaN at high pressures and high temperatures. Acta Physica Sinica, 2005, 54(12): 5830-5836. doi: 10.7498/aps.54.5830
    [19] Wang Chang-Qing, Jia Yu, Ma Bing-Xian, Wang Song-You, Qin Zhen, Wang Fei, Wu Le-Ke, Li Xin-Jian. Molecular dynamics simulations of various metastable structures on Si(001) at different temperatures. Acta Physica Sinica, 2005, 54(9): 4313-4318. doi: 10.7498/aps.54.4313
    [20] Hu Lin-Hua, Dai Song-Yuan, Wang Kong-Jia. Structural transformation of nanocrystalline titania grown by sol-gel technique and the growth kinetics of crystallites. Acta Physica Sinica, 2003, 52(9): 2135-2139. doi: 10.7498/aps.52.2135
Metrics
  • Abstract views:  5662
  • PDF Downloads:  122
Publishing process
  • Received Date:  16 August 2021
  • Accepted Date:  16 September 2021
  • Available Online:  23 January 2022
  • Published Online:  05 February 2022

/

返回文章
返回