Au纳米薄膜的熔化机理及其结构演变的分子动力学模拟

利用分子动力学模拟详细研究了不同厚度的Au纳米薄膜的熔化机理和结构演变. 模拟结果表明所有Au纳米薄膜的熔化行为分为两个阶段,即表面预熔和均相熔化. 只有最外层原子出现了预熔化行为, 其他内层原子在均相熔化之前始终保持稳定的固态,这与零维的Au纳米团簇和一维的Au纳米线的预熔化行为是不同的. 同时Au纳米薄膜的熔化温度随着薄膜厚度的增加而升高. 在预熔化过程中,在原子水平上发现了所有的Au纳米薄膜的f100g晶面向f111g晶面转变的表面重建过程. 对于最薄的L2纳米薄膜,当温度低于500 K 时表面应力不能诱导这样的表面重建. 然而一维的Au纳米线在更低温度下就能够观察到了由表面应力诱导的表面重建过程. 这主要是因为Au纳米线具有更高的比表面积所导致的. 另外研究结果还表明当模拟温度达到某一特定值时,由双原子层组成的Au纳米薄膜能够分裂成一维的纳米线.     

铂纳米晶在升温过程中结构演化与熔化特征的原子级模拟研究

采用分子动力学方法结合量子修正Sutton-Chen型多体力场,对由{100}面和{111}面构成的十四面体Pt纳米晶在升温过程中的热稳定性和熔化机制进行了计算机模拟研究,并引入统计半径和Lindemann指数来分析它的结构和形状演化过程. 结果表明：该纳米晶在1500 K时形状开始发生变化,并在1700 K时转变为球形. 铂纳米晶粒在1500 K时开始出现表面预熔,在1650 K时表面完全熔化并开始向内部传播,最终在1730 K时整体熔化为液态粒子. 表面预熔的出现对形状转变的发生是有利的. 关键词： 纳米晶 结构 熔化 分子动力学     

铂纳米粒子热稳定性的原子级模拟研究

本文采用分子动力学结合嵌入原子多体势,模拟了铂纳米粒子在升温过程中的热稳定性和熔化机制,并利用共近邻分析方法分析了它的微结构演化过程。模拟的结果表明：铂纳米粒子的熔点明显低于体材料的熔点;由于表面层原子的结合力较弱,在升温过程中表面会首先出现预熔;纳米粒子的熔化是从表面层开始的,并随着温度的升高,熔化的表面层会逐渐向内部扩展,最终导致纳米粒子整体转变为液态结构;当温度低于表面预熔温度时,纳米粒子保持良好的晶态结构。     

纳米团簇熔化过程的分子动力学模拟

本文采用分子动力学结合嵌入原子多体势,模拟了不同半径的Ni纳米团簇的升温熔化过程,研究团簇尺寸对熔点和表面能的影响.模拟结果表明:团簇的熔点显著低于体材料的熔点.团簇熔化的过程首先是在团簇的表面出现预熔,然后向团簇内部扩展,直到整个团簇完全熔为液态.在模拟的纳米尺度范围内,团簇的熔点与团簇尺寸基本成线性关系.团簇的表面能随着团簇尺寸的增大而减小,而且表面能均高于体材料的表面能.     

分子动力学模拟Cu(010)基体对负载Co-Cu双金属团簇熔化过程的影响

纳米团簇负载到基体上的结构演化和热稳定性是其走向技术应用的关键. 本文用分子动力学结合嵌入原子方法模拟了具有二十面体初始结构的Co₂₈₁Cu₂₈₀ 混合双金属团簇在Cu(010)基体上的熔化过程, 考察了基体的Cu原子可以自由移动(自由基体)和固定(固定基体)两种条件对负载团簇熔化的影响. 发现基体条件对团簇的熔化有明显的影响. 在自由基体上团簇原子的温度-能量曲线存在明显的团簇熔化时的能量突变点, 熔点为1320 K, 低于固定基体上团簇的熔点1630 K. 在升温过程中团簇的二十面体结构会在基体表面发生外延生长. 外延团簇随着温度增加发生表面预熔, 预熔原子会逐渐向基体表面扩散形成薄层, 直至完全熔化. 自由基体上团簇原子的嵌入行为会使原子的分布状态产生不同于固定基体上的演变.     

超细Pt纳米线结构和熔化行为的分子动力学模拟研究

基于EAM原子嵌入势, 对临界尺寸下的自由Pt纳米线的奇异结构和熔化行为进行分子动力学模拟. 模拟结果显示, 超细Pt纳米线的熔点随径向尺寸和结构的不同而发生明显改变; 引入林德曼因子, 令其临界值为0.03, 以此得到对应熔点值大小与通过势能-温度变化曲线找出的一致, 又比较了纳米线各层粒子平均林德曼指数的大小, 对各层纳米结构的热稳定性进行定量标度; 综合分析发现螺旋结构纳米线的熔化从内核开始, 而多边形结构的纳米线的熔化从外壳层开始.     

碳纳米管内金纳米线的结构与热稳定性

采用分子动力学模拟方法, 研究了填充在(8,8)单壁碳纳米管内的Au纳米线的结构和热稳定性. 研究表明, 经高温退火至室温, Au在碳纳米管内能生成多样而稳定的结构上明显区别于自由状态Au纳米线的壳层螺旋结构Au纳米线, 其螺旋结构会随着温度的变化而转变. 束缚在碳纳米管内的壳层螺旋结构Au纳米线有非常好的热稳定性, 稳定温度高于块体Au晶体的熔化温度. 关键词： 纳米线 碳纳米管 热稳定性 分子动力学模拟     

温度对金属纳米线势能分布的影响

采用三维分子动力学模拟方法,以面心立方金属银为研究对象,基于Finnis-Sinclair型嵌入原子法(EAM)多体势,模拟研究了纳米线势能分布特征在常温下及其在不同温度直到熔化过程中的变化,给出了常温及不同温度银纳米线势能分布比例和势能分布函数.结果表明:常温下,纳米线高势能原子比例随纳米线横截面尺寸的减小而增大,势能分布函数曲线各峰位几乎与纳米线横截面尺寸无关;纳米线熔化前的势能分布函数曲线具有多个波峰,随着温度增加,峰数减少且峰位右移;熔化后,多峰特征消失,只有一个宽化的峰.     

量子局域效应和应力对GaSb纳米线电子结构影响的第一性原理研究

采用基于密度泛函理论的第一性原理计算方法, 研究了不同晶体结构和尺寸的GaSb纳米线能带结构特性和载流子的有效质量, 以及单轴应力对GaSb纳米线能带结构的调控. 研究结果表明: 闪锌矿结构[111]方向和纤锌矿结构[0001]方向的小尺寸GaSb纳米线均出现间接带隙的能带结构, 并可通过单轴应力来实现纳米线能带结构由间接带隙到直接带隙的转变, 其中, 闪锌矿结构[111]方向GaSb纳米线仅在受到单轴拉伸应力时才发生能带由间接带隙到直接带隙的转变, 而纤锌矿结构[0001]方向GaSb纳米线无论受单轴拉伸还是压缩应力的作用均可实现能带由间接带隙到直接带隙的转变; [111]和[0001]方向GaSb纳米线的带隙和载流子有效质量与纳米线直径呈非线性关系, 并随纳米线直径的减小而增大; 同一方向和尺寸的GaSb纳米线, 其空穴有效质量要小于电子有效质量, 这表明小尺寸GaSb纳米线有利于空穴载流子输运.     

Pt-Ag合金纳米线热力学性质的分子动力学模拟研究

基于EAM原子嵌入势,采用分子动力学方法,对临界尺寸下的Pt0.95Ag0.05合金纳米线多边形结构的熔化行为进行了计算模拟.结果表明:径向尺寸对Pt0.95Ag0.05合金纳米线的熔点影响较为显著,而长度对其影响较小;引入林德曼因子得到的熔点和用势能-温度变化曲线找到的熔点基本一致;合金纳米线的染色原子由外向内运动;综合分析发现Pt0.95Ag0.05合金纳米线以先外后内的模式进行熔化.     

硅团簇熔化行为的紧束缚分子动力学研究

利用紧束缚分子动力学方法研究了硅团簇Si_n(n=5—10)的熔化行为.给出了团簇 熔化潜热 和熔点随团簇尺寸的变化关系，表明团簇熔化潜热和熔点强烈依赖于团簇的原子数.计算结 果表明硅团簇熔化机理与金属团簇熔化有很大不同，金属小团簇的熔化是一个从低温类固态 向高温类固态转变的过程，在转变温区，类固态和类液态处于动力学共存，而硅团簇在转变 温区则是处于一种中间态，这种中间态既不是类固态又不是类液态.比较了用不同计算方法 和定义方法所得硅团簇熔点. 关键词： 紧束缚 硅团簇 熔化潜热     

Thermal instability of decahedral structures in platinum nanoparticles

We conduct molecular dynamics simulations of 887 and 1389-atom decahedral platinum nanoparticles using an embedded atom potential. By constructing microcanonical caloric curves, we identify structural transitions from decahedral to fcc in the particles prior to melting. The transitions take place during phase coexistence and appear to occur via melting of the decahedral structure and subsequent recrystallisation into the fcc structure.     

Zn-Al二元合金熔化过程的低频内耗研究

利用强迫振动扭摆方法对ZnAl二元合金熔化过程的低频内耗进行了研究.结果表明,ZnAl二元合金熔化过程的内耗峰与其固态相变内耗峰的特征有较大差异.结合该合金熔化过程的微观结构变化,初步分析了内耗峰的形成机理. 关键词： Zn-Al合金 低频内耗 熔化     

Quantum mechanical aspect of first order phase transition of crystals

Classical molecular dynamics and Metropolis Monte Carlo simulations were carried out to investigate the thermal stability and melting behaviors of free-standing Pd-Pt bimetallic nanowires (NWs) with pentagonal multi-shell-type (PMS-type) structure in the whole composition range. Equilibrium configurations at 100 K are predicted in the semi-grand canonical ensemble. Pd-Pt PMS-type NWs are stable with a multishell structure of alternating Pd and Pt compositions and Pd segregating systematically to the surface. On thermal heating, an interesting composition-dependent structural transformation from the PMS-type to face-centred-cubic (FCC) by overcoming a high energy barrier is observed for Pd-Pt bimetallic NWs before the melting. Consequently, the system energy is decreased. The FCC structure is found more stable than PMS-type over the whole range of composition. The melting of Pd-Pt bimetallic NWs is also studied. It is found to start at the edges, then propagate over the whole surface, and next to the interior. It occurs in a composition-dependent range of temperature.     

Simulation of the processes of structuring of copper nanoclusters in terms of the tight-binding potential

The processes of melting and crystallization of copper nanoclusters with a radius ranging from 0.69 to 3.05 nm have been investigated using the molecular dynamics simulation. The performed simulation has shown that the melting begins with the surface of the cluster. Another feature of this phase transition is that it occurs in a temperature range where the liquid and solid phases can coexist. However, it is found that, for small copper clusters, the melting and crystallization temperatures coincide with each other. Moreover, it is established that the parent face-centered cubic structure of these small clusters (N < 150 atoms) transforms into a structure with fivefold symmetry even at temperatures of the order of 150–170 K. The behavior of some thermodynamic characteristics of copper nanoclusters is investigated in the vicinity of the solid-liquid phase transition. Analysis of the data obtained has revealed a number of regularities that are in agreement with the results of analytical calculations. In particular, the melting and crystallization temperatures of copper nanoparticles are linear functions of N ^?1/3. However, the melting heat ΔH _m and the melting entropy ΔS _m vary in a more complex manner. It is noted that the formation of a cluster structure depends on the conditions used for cooling from the liquid phase. Slow cooling results predominantly in the formation of a face-centered cubic phase, whereas rapid cooling in the majority of cases leads to the formation of an icosahedral modification. Therefore, the simulation performed has demonstrated the possibility of controlling the formation of a structure of copper nanoclusters during crystallization.     

Melting of lithium hydride under pressure

We have computed the melting line of lithium hydride up to 200 GPa using the two-phase simulation technique coupled with first-principles molecular dynamics. Our predicted melting temperature at high pressures varies slowly with compression, ranging from 2000 to 2450 K at 50-200 GPa pressures. The compressed fluid close to the melting line retains the ionic character of the low pressure molten state, while at higher temperatures dynamical hydrogen clustering processes are observed, which are accompanied by changes in the electronic structure.     

Neural Network Based on Quantum Chemistry for Predicting Melting Point of Organic Compounds

采用反向传播神经网络与定量结构-性质关系中量子化学参数相结合的方法测定有机化合物的熔点．11个反映分子间作用力和分子对称性的描述符作为输入变量,通过分子模型和PM3半经验分子轨道理论计算量子化学参数．用260个化合物训练由MatLab方法建立的神经网络,预测了73个化合物的熔点,并与文献中的实验数据进行比较,结果表明这种人工神经网络与定量结构-性质关系结合的方法可以预测有机化合物的熔点,平均绝对偏差5%．     

Two-dimensional mesoscopic vortex clusters in a superconducting ring: Shell structure and melting

The structure and phase transitions in the mesoscopic system of vortices in a quasi-two-dimensional superconducting ring are investigated. The shell structure of the mesoscopic system of vortices is studied, and its variation with the number of vortices and the parameters of the superconducting ring is analyzed. Two mechanisms of formation of new shells in vortex clusters with an increasing number of vortices in an increasing magnetic field are discovered: the generation of a new shell in a cluster and the splitting of the internal shell into two shells. The melting of vortex clusters and their thermodynamic parameters are analyzed using the Monte Carlo method. It is found that the melting of shell-type clusters occurs in two stages, orientation melting taking place at the lower temperature (during which nearly crystalline adjacent shells start rotating relative to each other) and blurring of the vortex structure occurring at the higher temperature. The shells obtained by splitting upon an increase in the number of vortices do not participate in orientational melting. The two-stage form of melting is associated with the smaller height of potential barriers being surmounted during the rotation of shells relative to one another as compared to the barrier for vortices jumping from one shell to another.     

Characteristic features of the melting of two-dimensional mesoscopic Wigner clusters

Two-dimensional Wigner microclusters in a semiconductor dot are studied. Their melting is investigated in detail and it is shown that, for typical mesoscopic clusters possessing a shell structure, melting occurs in two stages: orientational melting (rotation of the shells relative to one another) and total melting, where the shells start to overlap with one another and exchange particles. An example of a “magic” microstructure which has a triangular structure and melts in a single stage is presented. For this, the temperature dependences of various quantities characterizing cluster structure are investigated. The change in the distribution of cluster configurations over local minima of the potential energy with increasing temperature is investigated. At temperatures below the temperature of total melting, a cluster is always located near the configuration of a global minimum and, at temperatures above the temperature of complete melting, a cluster can be located with finite probability near configurations corresponding to various local minima of the potential energy. Fiz. Tverd. Tela (St. Petersburg) 41, 1499–1504 (August 1999)