Comparison of reactivity on step and terrace sites of Pd (3 3 2) surface for the dissociative adsorption of hydrogen: A quantum chemical molecular dynamics study

The notion of “active sites” is fundamental to heterogeneous catalysis. However, the exact nature of the active sites, and hence the mechanism by which they act, are still largely a matter of speculation. In this study, we have presented a systematic quantum chemical molecular dynamics (QCMD) calculations for the interaction of hydrogen on different step and terrace sites of the Pd (3 3 2) surface. Finally the dissociative adsorption of hydrogen on step and terrace as well as the influence of surface hydrogen vacancy for the dissociative adsorption of hydrogen has been investigated through QCMD. This is a state-of-the-art method for calculating the interaction of atoms and molecules with metal surfaces. It is found that fully hydrogen covered (saturated) step sites can dissociate hydrogen moderately and that a monovacancy surface is suitable for significant dissociative adsorption of hydrogen. However in terrace site of the surface we have found that dissociation of hydrogen takes place only on Pd sites where the metal atom is not bound to any pre-adsorbed hydrogen atoms. Furthermore, from the molecular dynamics and electronic structure calculations, we identify a number of consequences for the interpretation and modeling of diffusion experiments demonstrating the coverage and directional dependence of atomic hydrogen diffusion on stepped palladium surface.     

Investigation of the dissociative adsorption for cyclopropane on the copper surface by density functional theory and quantum chemical molecular dynamics method

The dissociative adsorption of cyclopropane on the copper surface was studied using quantum chemical molecular dynamics method with “Colors-Excite” code and density functional theory by Amsterdam Density Functional program (ADF2000). The excited state of cyclopropane was used as adsorbate to simulate the dissociated adsorption under an irradiation energy of ca. 10 eV. One of the C-C bonds in cyclopropane was broken and the two new bonds between cyclopropane and copper surface were formed. The electrons transferred from the copper atoms to cyclopropane with a value of about 0.2e. The shorter distances between the carbons and surface copper atoms showed the existence of strong interaction. Consistently, the results indicated metallacyclopentane was the most possible intermediate species in dissociative adsorption by ADF2000 and “Colors-Excite” method.     

Surface reconstruction on stishovite SiO2,HfO2 and rutile TiO2 （001）

This paper systematically investigates the surface reconstruction processes and patterns on stishovite SiO 2,HfO 2 and rutile TiO 2 (001) by using classical molecular dynamics.It is found that these three surfaces relax instead of reconstruction at 0 K,and have little possibility to reconstruct below 40 K.Above 40 K,surface reconstructions take place as collective atomic motion which can be speeded by higher temperature or compressed strain.Several reconstruction patterns with approximate surface energies are found,and electrostatic potentials on them are also provided in comparison with possible microscopic results.     

立方铂纳米粒子的形状变化与熔化特性的分子动力学研究

采用分子动力学方法结合嵌入原子多体势,对立方铂纳米粒子的热稳定性进行了模拟研究.计算结果表明,立方纳米粒子在升温过程中首先转变为由{111}和{100}面所构成的十四面体,然后再转变为球形,最后熔化为液态.通过计算立方铂钠米粒子的统计半径,发现形状转变温度在1250 K左右.尽管形状不同,立方纳米粒子和球形纳米粒子的熔点是相同的. 关键词： 纳米粒子 热稳定性 分子动力学     

A molecular dynamics study of the structural change differences between Au225 and Au369 clusters on MgO surfaces at low temperature

The differences in structural change between Au₂₂₅ and Au₃₆₉ clusters with their (111) facets supported on MgO(100) surfaces at 5 K are studied by using molecular-dynamics simulations with the atomic interchange potentials of the Au/MgO interface. The parameters are obtained from the ab initio energies using the Chen-Möbius inversion method. Analyses of the pair distribution functions show that the two Au clusters use different deformation processes to adjust the distances between the interface atoms, owing to the misfit between the atom distances among the clusters and the substrates. The local structural changes are identified by atomic density profiles.     

A quantum chemical study of ZrO2 atomic layer deposition growth reactions on the SiO2 surface

Zirconium oxide (ZrO₂) is one of the leading candidates to replace silicon oxide (SiO₂) as the gate dielectric for future generation metal-oxide-semiconductor (MOS) based nanoelectronic devices. Experimental studies have shown that a 1–3 monolayer SiO₂ film between the high permittivity metal oxide and the substrate silicon is needed to minimize electrical degradation. This study uses density functional theory (DFT) to investigate the initial growth reactions of ZrO₂ on hydroxylated SiO₂ by atomic layer deposition (ALD). The reactants investigated in this study are zirconium tetrachloride (ZrCl₄) and water (H₂O). Exchange reaction mechanisms for the two reaction half-cycles were investigated. For the first half-reaction, reaction of gaseous ZrCl₄ with the hydroxylated SiO₂ surface was studied. Upon adsorption, ZrCl₄ forms a stable intermediate complex with the surface SiO₂–OH^* site, followed by formation of SiO₂–O–Zr–Cl^* surface sites and HCl. For the second half-reaction, reaction of H₂O on SiO₂–O–Zr–Cl^* surface sites was investigated. The reaction pathway is analogous to that of the first half-reaction; water first forms a stable intermediate complex followed by evolution of HCl through combination of a Cl atom from the surface site and an H atom from H₂O. The results reveal that the stable intermediate complexes formed in both half-reactions can lead to a slow film growth rate unless process parameters are adjusted to lower the stability of the complex. The energetics of the two half-reactions are similar to those of ZrO₂ ALD on ZrO₂ and as well as the energetics of ZrO₂ ALD on hydroxylated silicon. The energetics of the growth reactions with two surface hydroxyl sites are also described.     

熔融Cu55团簇在Cu(010)表面上凝固过程的分子动力学模拟

采用基于嵌入原子方法的分子动力学,模拟了熔融Cu₅₅团簇在Cu衬底（010）表面上以两个不同降温速率降温过程中结构的变化.模拟结果表明,降温速率对团簇结构的变化有很大影响.较快的降温速率使得降温过程中团簇原子具有较低的能量;较慢的降温速率有助于高温时位于衬底内的原子向衬底表面扩散,排列形成面心立方结构. 关键词： 团簇 凝固 分子动力学 表面     

与HO2和HO3相关的复杂反应的量子动力学

Complex-forming reactions widely exist in gas-phase chemical reactions.Various complexforming bimolecular reactions have been investigated and interesting phenomena have been discovered.The complex-forming reactions usually have small or no barrier in the entrance channel,which leads to obvious differences in kinetic and dynamic characteristics compared with direct reactions.Theoretically,quantum state-resolved reaction dynamics can provide the most detailed microscopic dynamic mechanisms and is now feasible for a direct reaction with only one potential barrier.However,it is of great challenge to construct accurate potential energy surfaces and perform accurate quantum dynamics calculations for a complex polyatomic reaction involving deep potential wells and multi-channels.This paper reviews the most recent progress in two prototypical oxyhydrogen complex-forming reaction systems,HO₂ and HO₃,which are significant in combustion,atmospheric,and interstellar chemistry.We will present a brief survey of both computational and experimental work and emphasize on some unsolved problems existing in these systems.     

Adsorption of NO on unreduced and reduced CeO2 surfaces: A plane-wave DFT study

The adsorption of NO on the (1 1 1) and (1 1 0) surfaces of ceria (CeO₂) was studied using projector-augmented wave (PAW) method based density-functional theory within the generalized gradient approximation (GGA). Several adsorption sites for NO on the stoichiometric surfaces are found, all with weak molecule-surface interaction. The adsorption on the reduced surfaces is much stronger. The O-ends of the adsorbed NO molecules fill the oxygen vacancies and the N-O bonds are elongated. If two such adsorbed NO molecules, residing at neighbouring sites, meet, their N-ends will form a strong N-N bond with little or no barrier. This is an intermediate step towards dissociation of free N₂ which is calculated to be strongly thermodynamically driven.     

UO2 晶体中低密勒指数晶面表面能的分子动力学模拟

采用基于刚性离子势的分子动力学模拟方法初步计算了UO₂晶体中(100), (110)和(111) 3种低密勒指数晶面在300–1500 K范围内的表面能大小. 结果表明, 3种晶面的表面能大小随温度的升高而降低, 与实验结果趋势一致; 原子排列最紧密的(111)晶面具有最低的表面能, 3种晶面的表面能大小从高到低依次为(100), (110)和(111)晶面; 达到平衡状态下的表面层原子相对于体内原子层在表面的法线方向上发生了明显的压缩并且表面层原子的对称性也降低了, 表面原子的弛豫效应一直影响到了第5层. 计算研究结果将有助于深入认识UO₂燃料中裂变气体气泡的聚集长大以及燃料的辐照肿胀开裂行为. 关键词： 分子动力学 2')" href="#">UO₂ 低密勒指数晶面 表面能     

Ab initio path-integral molecular dynamics and the quantum nature of hydrogen bonds

The hydrogen bond(HB) is an important type of intermolecular interaction, which is generally weak, ubiquitous,and essential to life on earth. The small mass of hydrogen means that many properties of HBs are quantum mechanical in nature. In recent years, because of the development of computer simulation methods and computational power, the influence of nuclear quantum effects(NQEs) on the structural and energetic properties of some hydrogen bonded systems has been intensively studied. Here, we present a review of these studies by focussing on the explanation of the principles underlying the simulation methods, i.e., the ab initio path-integral molecular dynamics. Its extension in combination with the thermodynamic integration method for the calculation of free energies will also be introduced. We use two examples to show how this influence of NQEs in realistic systems is simulated in practice.     

Classical and quantum molecular dynamics of cation in (CH3NH3)3 Sb2Br9 polycrystal as studied by 1H NMR.

¹H spin-lattice relaxation times and second moments were determined for polycrystalline (CH₃NH₃)₃Sb₂Br₉ sample in a wide range of temperature (5–200 K) at 24.6 and 55.2 MHz. ²H NMR spectra of (CD₃NH₃)₃Sb₂Br₉ were recorded between 5 K and room temperature. The relaxation time is interpreted as a result of motion of two different non-equivalent types of monomethylammonium cations occurring at the 2:1 proportion in a unit cell. Below 30 K, the relaxation processes via tunneling are suggested to dominate. Above 30 K, only classical behaviour of methylammonium cations is detected. Two monomethylammonium cations relax with the classical correlated C₃ reorientation and the rotational tunnelling mechanism, while the third cation exhibits only the classical correlated reorientation. The dynamic parameters of these motions have been determined.     

Coalescence between Cu57 and Cu58 clusters at a room temperature: molecular dynamics simulations

Three coalescence processes of Cu57-Cu57,Cu57-Cu58,and Cu58-Cu58 clusters at 300 K are investigated by employing molecular dynamics simulations.According to the evolutions of mean square displacement and local atom packing,the coalescence process can be separated into three stages including an approaching stage,a coalescing stage,and a coalesced stage.The simulations show that the coalescence processes and the formed products are sensitive to the respective initial structures of,and the relative configuration between,the two coalescing icosahedron-based clusters.     

Study of reduction processes over cerium oxide surfaces with atomic hydrogen using ultra accelerated quantum chemical molecular dynamics

Ceria plays an important role in catalysis, due to its ability to store and release oxygen depending on the condition present in the catalyst environment. To analyze the role of ceria in catalytic reactions, it is necessary to know the details of the interaction of ceria surface with environmentally sensitive molecules. This study was conducted using ultra accelerated quantum chemical molecular dynamics. Its purpose was to investigate the reduction process of the (1 1 1) and (1 1 0) surfaces of ceria with atomic hydrogen as well as water desorption mechanisms from the surfaces. This simulation demonstrated that when a high-energy colliding hydrogen atoms are adsorbed on the ceria, it pulls up an O atom from the ceria surfaces and results in the formation of a H₂O molecule. This is the first dynamics simulation related to such reduction processes based on quantum chemistry.     

Adsorption mechanisms of translationally-energetic O2 and N2: direct dissociation versus direct molecular chemisorption

A direct dissociation mechanism has been traditionally assigned to molecular beam data that exhibit an increase in the initial adsorption probability with increasing kinetic energy. Yet, recent experiments of nitrogen and oxygen adsorption provide support for an alternative high kinetic energy pathway in which incident energy assists in surmounting barriers to molecular chemisorption on a surface as the first step to dissociation. Moreover, systems for which the experimental evidence supports such a mechanism also demonstrate that molecularly chemisorbed intermediates can be spectroscopically observed at low temperatures and coverages from exposure to a gas in thermal equilibrium at room temperature. Likewise, such observations have not been measured for systems which are consistent with direct dissociation. A consideration of this trend regarding the existence of molecularly chemisorbed states and the implications for the dominant, dissociative chemisorption pathway at high kinetic energy is presented for a number of gas surface systems.     

O-vacancy and surface on CeO2: A first-principles study

Atomic and electronic structures of CeO₂ (1 1 1), (1 1 0) and (1 0 0) surfaces are investigated using the first-principles density functional theory taking into account the on-site Coulomb interaction. Both the stoichiometric and O-deficient surfaces are examined in order to clarify the overall features. The CeO₂ (1 1 1) is found to be the most stable surface, followed by the (1 1 0) and (1 0 0) surfaces, consistent with experimental observations. Three surfaces exhibit different features of relaxation. Large relaxations are found at the (1 1 0) and (1 0 0) surfaces, while very small changes are observed at the (1 1 1) surface. It is found that the O-vacancy occurs more readily at the (1 1 0) surface as compared with the (1 1 1) surface. Furthermore, the formation energies of the O-vacancy in the surfaces are lower than that in the bulk. The energetically favorable O-vacancy locates in the second O-atomic layer for the (1 1 1) while at the surface layer for the (1 1 0). The excess electrons left with the removal of the O atom are distributed in the first two layers with certain (a considerable) fraction filling the Ce-4f states.     

Kinetic modelling of the H2---O2 reaction on Pd and of its influence on the hydrogen response of a hydrogen sensitive Pd metal-oxide-semiconductor device

We have established a model for the water forming reaction on Pd in the temperature range 350–475 K. Importantly, the model takes into account the possibility that hydrogen may absorb and adsorb at interface sites on supported Pd catalysts. It is shown that already at modest conditions interface adsorption may significantly affect reaction rates. The model may also be used to quantify the response of a hydrogen sensitive Pd-MOS device during hydrogen sensing in oxygen. In the case of Pd supported on SiO₂, the concentration of interface sites is so low that interface hydrogen adsorption will have only a minor influence on a catalytic reaction. The fact that a Pd-MOS device may be used as a very sensitive hydrogen detector at atmospheric oxygen conditions, despite a steric oxygen blocking of hydrogen dissociation sites, is predicted by the model.     

低温下Au959团簇负载于MgO(100)表面后结构变化的分子动力学研究

采用基于Chen-Mbius反演方法,从金属/金属氧化物界面第一性原理计算的粘结能结果中推导出的Au/MgO原子间相互作用势的正则系综（NVT）分子动力学,模拟了在10 K条件下,Au₉₅₉团簇负载于MgO（100）表面后团簇结构的变化.根据原子对分析技术和对分布函数的分析表明,由于团簇界面处原子间距与载体原子间距相匹配,置于载体上的Au团簇经过一个变形过程后,较其孤立自由表面时的团簇体积变大. 关键词： 团簇 分子动力学 计算机模拟 表面