不同热浴下Co团簇熔化行为的分子动力学模拟

运用分子动力学方法,采用Berendsen热浴和Nose-Hoover热浴分别研究了Co_n(n=13,55,147)团簇的熔化特性,模型采用Gupta相互作用势.模拟结果表明:两种热浴对钴团簇熔点及预熔化区间给出了基本一致的描述.所研究团簇体系在给定温度下长时间内各Co团簇中单个原子的速率(速度)分布与麦克斯韦速率(速度)分布曲线符合很好.     

大尺寸贵金属纳米团簇热力学性质的分子动力学模拟研究

采用分子动力学方法和原子嵌入势模拟了大尺寸金(n=1136～1556)、银(n=1088～1724)、铜(n=1000～1600)、铂(n=1004～1800)原子纳米团簇的熔化过程,得出了相应纳米团簇的势能随温度的变化曲线以及热容量随温度的变化曲线,研究了各种原子纳米团簇熔点与其团簇尺寸的关系.模拟结果表明团簇的熔点随团簇尺寸增大而升高,并逐渐向大块晶体靠拢.所有纳米团簇在熔化过程中在熔点附近都出现负热容现象,通过对团簇熔化前后结构的比较,分析了导致这种现象的原因.     

仲氢及正氘团簇的量子定域和超流动性研究

本文采用路径积分Monte Carlo (PIMC)方法研究仲氢团簇和正氘团簇的结构和超流体性质.比较仲氢团簇(分子数N≤40)在T=0.4 K,T=1.6 K的临近能量差和超流比例(superfluid fractions),发现仲氢团簇在T=0.4 K时出现"量子熔化"(quantum melting);研究了仲氢团簇的量子定域(quantum localization)和 "量子熔化"的关系,对T=0.4 K时仲氢及正氘团簇的量子定域与超流动性 关键词： 团簇 量子定域 量子熔化 超流比例     

Aln(n=13–32)团簇熔化行为的分子动力学模拟研究

基于半经验的Gupta原子间多体相互作用势, 采用分子动力学方法并结合模拟退火及淬火技术, 系统研究了小尺寸铝团簇Al_n (n=13–32)的熔化行为. 模拟结果表明: 除个别尺寸(Al₁₃ 和Al₁₉) 外, 团簇熔化过程比热曲线普遍呈现杂乱无规(无明显单峰)现象, 这与实验观测小Al团簇比热普遍无规的结果完全一致. 通过分析不同温度点上团簇淬火构型的势能分布图给出了小Al团簇比热呈现无规或有规现象的成因. 对于比热无规团簇, 可以利用原子等价指数判断给出团簇熔点, 所得团簇熔点随团簇尺寸的变化趋势与实验观测结果完全一致. 关键词： Gupta势 团簇 分子动力学 熔化     

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

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

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

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

支撑Co团簇熔化行为的分子动力学模拟

运用分子动力学结合退火方法,模拟研究了所含原子数介于400～20000之间的自由及支撑Co团簇的熔化行为.模拟中团簇Co原子间相互作用采用半经验的Gupta多体相互作用势.讨论了不同衬底势(Lennard-Jones势与Morse势)对支撑团簇熔化行为的影响.结果表明:对于给定原子数目的支撑Co团簇,在衬底势强度相同时,两种衬底势下对支撑团簇熔点及预熔化区间给出了一致的描述;随衬底势强度增加,支撑团簇熔点升高,且两种衬底势所描述的支撑Co团簇也类同于自由Co团簇都表现出较好的熔点与尺寸依赖的线性关系.     

铜、银和铂原子纳米团簇负热容现象的分子动力学模拟研究

本文采用微正则分子动力学方法模拟研究了铂、铜和银原子纳米团族从固态到液态的熔化过程,得到热容量随温度变化关系,结果表明这三种金属纳米团簇在熔化过程中均出现了负热容现象,并通过对团簇热能随温度的变化关系以及团簇原子数径向分布的分析,探讨了产生负热容现象的微观机制.     

ConCu(38-n)(n=0～38)团簇基态结构和熔化行为

基于半经验的Gupta多体势采用遗传算法和分子动力学方法并结合模淬火技术,系统研究了ConCu(38-n)(n=0～38)团簇的基态结构与熔化行为．结果表明：除Co7Cu31团簇基态为类Ih结构外,其它ConCu(38-n)混合团簇的基态结构均是在单质Cu38（Co38）的Oh基态结构基础上的畸变;对于混合团簇的基态结构,随Co原子增加均表现出Co原子先占中心后占表面的特征;通过分析基态团簇二阶差分能和混合能表明Co7Cu31和Co14Cu24为幻数结构团簇;ConCu(38-n)混合团簇在熔化过程中均表现出热容曲线无明显宽峰的异常熔化行为;通过对团簇原子等价指数和淬火结构势能分布图对团簇异常熔化行为进行了分析表征,指出在团簇熔化过程中两种动力学稳定结构（类Oh结构与类Ih结构）之间的相互竞争对团簇熔化行为的重要影响．     

小Cu团簇熔化过程中的预熔化行为

用分子动力学方法研究了不同大小的Cu原子团簇(Cu13,Cu38,Cu55)的熔化行为.发现对于小的团簇(Cu13),在接近熔化温度时,其内部存在着一种特殊的预熔化行为:内部原子的相对扭转,这种扭转使得Cu13团簇在保持基态结构不变的同时,其内部的原子又具有一定的流动性,对于较大的团簇(Cu38,Cu55)则不存在这种现象.     

Dynamical consequences of a constraint on the Langevin thermostat in molecular cluster simulation

We investigate some unusual behaviour observed while performing molecular dynamics simulations of small molecular clusters using a constrained Langevin thermostat. Atoms appear to be thermalised to different temperatures that depend on their mass and on the total number of particles in the system. The deviation from the zeroth law of thermodynamics can be considerable for small systems of heavy and light particles. We trace this behaviour to the absence of thermal noise acting on the centre of mass of the system. This is demonstrated by solving the stochastic dynamics for the constrained thermostat and comparing the results with simulation data. By removing the constraint, the Langevin thermostat may be restored to its intended behaviour. We also investigate a Langevin thermostat constrained to have zero total force acting on its centre of mass, and find similar deficiencies.     

Structural properties of complex (dusty) plasma upon crystallization and melting

A change in the local order of a bounded complex (dusty) plasma in the process of its crystallization and melting has been examined by molecular dynamics simulations. The dynamics of microparticles is considered in the framework of a Langevin thermostat, the pair interaction between charged particles is described by a screened Coulomb potential (Yukawa potential) with the hard wall potential as a confinement. It has been shown that the beginning of the crystallization of such a system is accompanied by the formation of clusters with the hexagonal close packed (hcp) structure; a noticeable number of these clusters are then transformed to the face centered cubic (fcc) phase. A plasma crystal formed after crystallization consists of the metastable hcp phase, fcc clusters, and a small number of clusters with a body centered cubic (bcc) crystal lattice. Beginning with a certain threshold value of the thermostat temperature, the number of fcc/bcc clusters decreases sharply with increasing temperature, which is an important signature of the beginning of the melting of the plasma crystal.     

Size dependence of multipolar plasmon resonance frequencies and damping rates in simple metal spherical nanoparticles

We have simulated the non-equilibrium dynamics of methanol adsorbed in FAU zeolite driven by external microwave (MW) radiation. We have modelled steady states produced by augmenting this MW-driven system with a thermostat that acts as a balancing heat sink. We have compared results from an implicit thermostat (Andersen velocity replacement) and an explicit thermostat (helium atoms subjected to Andersen velocity replacement). We find very good agreement between the implicit and explicit thermostats for energy distributions and diffusion coefficients produced under MW-heated steady-state conditions. This augurs well for the continued use of implicit thermostats, which are computationally more efficient.     

Impact of thermostat on interfacial thermal conductance prediction from non-equilibrium molecular dynamics simulations

The knowledge of interfacial thermal conductance (ITC) is key to understand thermal transport in nanostructures. The non-equilibrium molecular dynamics (NEMD) simulation is a useful tool to calculate the ITC. In this study, we investigate the impact of thermostat on the prediction of the ITC. The Langevin thermostat is found to result in larger ITC than the Nose-Hoover thermostat. In addition, the results from NEMD simulations with the Nose-Hoover thermostat exhibit strong size effect of thermal reservoirs. Detailed spectral heat flux decomposition and modal temperature calculation reveal that the acoustic phonons in hot and cold thermal reservoirs are of smaller temperature difference than optical phonons when using the Nose-Hoover thermostat, while phonons in the Langevin thermostat are of identical temperatures. Such a non-equilibrium state of phonons in the case of the Nose-Hoover thermostat reduces the heat flux of low-to-middle-frequency phonons. We also discuss how enlarging the reservoirs or adding an epitaxial rough wall to the reservoirs affects the predicted ITC, and find that these attempts could help to thermalize the phonons, but still underestimate the heat flux from low-frequency phonons.     

Numerical distortion and effects of thermostat in molecular dynamics simulations of single-walled carbon nanotubes

In this paper, single-walled carbon nanotubes （SWCNTs） are studied through molecular dynamics （MD） simulation. The simulations are performed at temperatures of 1 and 300K separately, with atomic interactions characterized by the second Reactive Empirical Bond Order （REBO） potential, and temperature controlled by a certain thermostat, i.e. by separately using the velocity scaling, the Berendsen scheme, the Nose-Hoover scheme, and the generalized Langevin scheme. Results for a （5,5） SWCNT with a length of 24.5 nm show apparent distortions in nanotube configuration, which can further enter into periodic vibrations, except in simulations using the generalized Langevin thermostat, which is ascribed to periodic boundary conditions used in simulation. The periodic boundary conditions may implicitly be applied in the form of an inconsistent constraint along the axis of the nanotube. The combination of the inconsistent constraint with the cumulative errors in calculation causes the distortions of nanotubes. When the generalized Langevin thermostat is applied, inconsistently distributed errors are dispersed by the random forces, and so the distortions and vibrations disappear. This speculation is confirmed by simulation in the case without periodic boundary conditions, where no apparent distortion and vibration occur. It is also revealed that numerically induced distortions and vibrations occur only in simulation of nanotubes with a small diameter and a large length-to-diameter ratio. When MD simulation is applied to a system with a particular geometry, attention should be paid to avoiding the numerical distortion and the result infidelity.     

Monte Carlo simulations of copper clustering in Fe–Cu alloys under irradiation

We present the computational approach for studying the microstructures of Cu clusters in Fe–Cu alloys by combining the molecular dynamics (MD) simulation and Monte Carlo methods. The MD simulation is used to characterize the primary damage resulting from the displacement cascade in Fe. Then, using the Metropolis Monte Carlo methods, the microstructure of the Cu clusters is predicted under the assumption that the system will evolve towards the equilibrium state. The formation of the Cu clusters is apparent for Fe–Cu alloys of a higher Cu content (1.0 w/o), whereas the degree of Cu clustering is not significant for the lower Cu content (0.1 w/o) alloys. The atomic configuration of the Cu–vacancy complex under irradiation, produced by this simulation, is in a fair agreement with the experiments. The simulation is expected to provide important information on the Cu-cluster morphology, which is useful for experimental data analysis.     

Dopant dependent stability of Co n TM+ (TM = Ti,V, Cr,and Mn) clusters

In this paper, we present a photofragmentation study of mass-selected transition metal-doped cobalt cluster cations Co _n TM⁺  (n = 8–18, TM = Ti, V, Cr, and Mn). Time-of-flight spectra recorded after laser excitation of mass-selected clusters in the gas phase show that the evaporation of a cobalt atom is the most facile dissociation channel for clusters with TM = Ti and V, suggesting an enhanced stability of the doped clusters compared to the bare ones. In contrast, for Co _n TM⁺ with TM = Cr and Mn, the loss of the dopant atom is found to be the preferred dissociation channel. Co₁₃Cr⁺ is a notable exception and favors dissociation by loss of a neutral Co atom. It is implied that substituting Mn and Cr generally destabilizes the cobalt clusters with the exception of Co₁₂Cr⁺, which is relatively more stable than Co ₁₃ ⁺ . Additional measurements of V _n Co⁺ (n = 9–16) show that the loss of a Co atom is still the most facile dissociation channel, which is in agreement with the predicted stronger V?V bond compared to the V?Co one.     

Kinetic Monte Carlo simulation of the growth of metal clusters on regular array of defects on insulator

A Kinetic Monte Carlo simulation of the nucleation and growth of Pd clusters on a nanostructured alumina substrate is presented. The new Monte Carlo simulation program allows to derive the 3D shape of the growing clusters without performing a full all atoms simulation. The simulation shows, like in previous pure 2D simulations, that clusters nucleate exclusively on the defects of the nanostructure in a limited range of substrate temperature. Around 300 K, the clusters have a compact faceted shape and they grow, at not too large coverage, layer by layer. These results are in agreement with previous studies of the nucleation and growth of Pd clusters on an ultrathin alumina film on Ni₃Al (1 1 1).     

X-Ray Photoelectron and Raman Spectroscopies,Comparative Studies of Supported Molybdenum Cobalt and Cobalt Molybdenum Catalysts

Laser Raman Spectroscopy is used to study catalysts : Co, Mo or Co-Mo supported on alumina. Well defined salts spectra are compared to those of these solids, which permits to establish superficial structures of cobalt, molybdenum or both.

For molybdenum catalysts, we found two structures depending upon the concentration

isolated MoO^2- ₄ molecular ions at low concentration < 4%

polymolybdate phase at higher concentration

For cobalt catalysts two structures are also observed: - tetrahedrally surrounded co ²+ ions at low concentration Co³O⁴ clusters at higher concentrations

For cobalt-molybdenum (8%) catalysts three cases occur no effect is observable at very low concentration of cobalt (0.5%)

a well defined compound appears for Co/Mo atomic ratio 0.25

when the cobalt concentration increases a third structure is observed but not identified

All these results are compared to X Ray photoelectron spectroscopic studies and are in good agreement with them.