[论文解读] Melting of sodium clusters

Juan A. Reyes Nava, Ignacio L. Garzón+2 more|ArXiv.org|Jul 17, 2002

nanoparticles nucleation surface interactionsEarth and Planetary Sciences参考文献 2被引用 16

一句话总结

本研究利用基于多体Gupta势的微正则系综分子动力学模拟，研究了钠团簇（Naₙ，N = 13–147）的熔化转变过程，该势函数基于紧束缚哈密顿量的二阶矩近似，用于描述金属键合。研究揭示了尺寸依赖的熔化机制、不规则的熔点趋势，以及Na₅₅的计算熔点比实验值低约40%，突显了经验势在定量捕捉团簇熔化行为方面的局限性，尽管其在定性上与实验趋势相符。

ABSTRACT

Thermal stability properties and the melting-like transition of Na_n, n=13-147, clusters are studied through microcanonical molecular dynamics simulations. The metallic bonding in the sodium clusters is mimicked by a many-body Gupta potential based on the second moment approximation of a tight-binding Hamiltonian. The characteristics of the solid-to-liquid transition in the sodium clusters are analyzed by calculating physical quantities like caloric curves, heat capacities, and root-mean-square bond length fluctuations using simulation times of several nanoseconds. Distinct melting mechanisms are obtained for the sodium clusters in the size range investigated. The calculated melting temperatures show an irregular variation with the cluster size, in qualitative agreement with recent experimental results. However, the calculated melting point for the Na_55 cluster is about 40 % lower than the experimental value.

研究动机与目标

• 通过长时间尺度的微正则系综分子动力学模拟，研究钠团簇的固-液转变过程。
• 检验多体Gupta势在长时间模拟（最长达50 ns）下描述熔化转变的性能。
• 将MD模拟获得的熔点与熔化机制，与蒙特卡罗模拟及第一性原理方法的结果进行比较。
• 分析表面熔化与熔化过程中结构演变在Naₙ团簇中的作用。
• 评估长时间模拟下物理可观测量（如热容曲线、比热容和均方根键长涨落）的收敛性。

提出的方法

• 对Naₙ团簇（N = 13, 20, 55, 135, 142, 147）进行了微正则系综分子动力学模拟，采用基于紧束缚哈密顿量二阶矩近似的多体Gupta势。
• 势能函数定义为 Vi = A∑j≠i e^(-p(r_ij/r₀ - 1)) - ξ(∑j≠i e^(-2q(r_ij/r₀ - 1)))^(1/2)，通过有效双体与三体相互作用描述金属键合。
• 在约50纳秒的模拟时间内，计算了热容曲线、比热容和均方根（rms）键长涨落等物理可观测量。
• 采用两种标准确定熔点：比热容的峰值与Lindemann准则（δ ≈ 0.15）。
• 通过追踪原子轨迹分析加热过程中的结构演变，识别出表面熔化与完全扩散运动等阶段。
• 模拟初始结构采用能量最低的异构体，包括二十面体结构及其不完整变体（如由Na₁₄₇衍生的Na₁₃₅、Na₁₄₂）。

实验结果

研究问题

• RQ1 钠团簇的熔点如何随尺寸变化？该变化是否与实验观测一致？
• RQ2 Naₙ团簇中的主导熔化机制是什么？其从表面熔化到体相熔化的演化过程如何？
• RQ3 长时间模拟（最长50 ns）如何影响熔化指标（如比热容与Lindemann参数）的收敛性？
• RQ4 为何Na₅₅团簇的计算熔点显著低于实验值？这对Gupta势的准确性有何启示？
• RQ5 微正则MD模拟的结果与蒙特卡罗模拟及第一性原理方法在定性和定量上如何比较？

主要发现

• Na₅₅团簇的计算熔点比实验值低约40%，表明存在显著的定量偏差。
• 熔点随团簇尺寸表现出不规则依赖关系，其中Na₁₄₂的熔点最高（通过比热容峰值测得为190 K），与实验趋势一致。
• Lindemann准则给出的熔点低于比热容峰值，表明结构软化先于完全扩散运动的出现。
• 熔化始于团簇表面，原子运动性由表面向内传播，证实表面熔化是体相熔化的前驱过程。
• 无论起始于完整还是不完整的二十面体结构，Na₁₃₅、Na₁₄₂与Na₁₄₇的热容曲线、比热容与rms键长涨落几乎完全一致，表明其结构具有强鲁棒性。
• 更长的模拟时间（最长50 ns）导致更低的熔点，表明时间平均可观测量的收敛性得到改善。