Quantum Monte Carlo study of the CO interaction with a dimer model surface for Cr(110)

The chemisorption of CO on a Cr( 110) surface is investigated using the quantum Monte Carlo method in the diffusion Monte Carlo (DMC) variant and a model Cr2CO cluster. The present results are consistent with the earlier ab initio HF study with this model that showed the tilted/ near-parallel orientation as energetically favoured over the perpendicular arrangement. The DMC energy difference between the two orientations is larger (1.9 eV) than that computed in the previous study. The distribution and reorganization of electrons during CO adsorption on the model surface are analysed using the topological electron localization function method that yields electron populations, charge transfer and clear insight on the chemical bonding that occurs with CO adsorption and dissociation on the model surface.

有分子内态激发和化学反应的流动的数值模拟

本文旨在研究双原子分子化学反应的微观模型、表面复合反应机制和分子 相继碰撞中的能量传递过程等问题。这些问题的本质是如何考虑分子和原子之 间的相互作用，传统的研究途径是从能量观点出发的。本文尝试将力学的观点 和方法，与分子模型和统计平均手段结合起来，在计算机的帮助下，分析和研 究这些问题。 本文所建立的依赖空间取向化学反应模型认为：化学键断裂的充要条件是 作用于化学键的拉伸力超过它的抗拉极限。由此导出的双原子分子的反应微观 判据及反应速率常数Arrhenius形式的表达式，反映了碰撞空间取向以及反应势 能面的影响。对这个模型的检验和应用表明：它能够解释了置换反应速率常数 实验数据拟合公式k_f=aTbexp(-Ts/T)中Ts<ab initio)方法。它的突出特点是能够模拟分子之间的相继碰撞，直接给出了分子转动能或振动能随碰撞次数的变化过程，这对于了解和认识热力学非平衡现象特别有帮助口对同核双原子分子和惰性气体分子碰撞 过程的模拟结果，与实验或理论结果符合，说明单个靶分子方法的基本想法和 实施步骤，以及我们对惰性气体分子和同核双原子分子及其原子相互作用的考 虑，是现实可行的。 此外，本文附录还简要介绍了作者在读博士期间完成的其它工作，包括一 种自然的宏观流动量分解方法及封闭的平均流动方程、过渡领域三维复杂外形 绕流场的直接统计模拟以及过渡领域高超声速圆柱绕流的直接统计模拟等。

Mechanical and chemical bonding properties of ground state BeH2

The crystal structure, mechanical properties and electronic structure of ground state BeH2 are calculated employing the first-principles methods based on the density functional theory. Our calculated structural parameters at equilibrium volume are well consistent with experimental results. Elastic constants, which well obey the mechanical stability criteria, are firstly theoretically acquired. The bulk modulus B, Shear modulus G, Young's modulus E, and Poisson's ratio upsilon are deduced from the elastic constants. The bonding nature in BeH2 is fully interpreted by combining characteristics in band structure, density of states, and charge distribution. The ionicity in the Be-H bond is mainly featured by charge transfer from Be 2s to H 1s atomic orbitals while its covalency is dominated by the hybridization of H 1s and Be 2p states. The Bader analysis of BeH2 and MgH2 are performed to describe the ionic/covalent character quantitatively and we find that about 1.61 (1.6) electrons transfer from each Be (Mg) atom to H atoms.

Effects of residual C and O impurities on photoluminescence in undoped GaN epilayers

Photoluminescence (PL) was investigated in undoped GaN from 4.8 K to room temperature. The 4.8 K spectra exhibited recombinations of free exciton, donor-acceptor pair (DAP), blue and yellow bands (Ybs). The blue band (BB) was also identified to be a DAP recombination. The YB was assigned to a recombination from deep levels. The energy-dispersive X-ray spectroscopy show that C and O are the main residual impurities in undoped GaN and that C concentration is lower in the epilayers with the stronger BB. The electronic structures of native defects, C and O impurities, and their complexes were calculated using ab initio local-density-functional (LDF) methods with linear muffin-tin-orbital and 72-atomic supercell. The theoretical analyses suggest that the electron transitions from O-N states to C-N and to V-Ga states are responsible for DAP and the BB, respectively, and the electron transitions between the inner levels of the C-N-O-N complex may be responsible for the YB in our samples. (C) 2002 Elsevier Science B.V. All rights reserved.

Energy bands and acceptor binding energies of GaN

The energy bands of zinc-blende and wurtzite GaN are calculated with the empirical pseudopotential method, and the pseudopotential parameters for Ga and N atoms are-given. The calculated energy bands are in agreement with those obtained by the ab initio method. The effective-mass theory for the semiconductors of wurtzite structure is established, and the effective-mass parameters of GaN for both structures are given The binding energies of acceptor states are calculated by solving strictly the effective-mass equations. The binding energies of donor and acceptor are 24 and 142 meV for the zinc-blende structure, 20 and 131, and 97 meV for the wurtzite structure, respectively, which are consistent with recent experimental results. It is proposed that there are two kinds of acceptor in wurtzite GaN. One kind is the general acceptor such as C, which substitutes N, which satisfies the effective-mass theory. The other kind of acceptor includes Mg, Zn, Cd, etc., the binding energy of these accepters is deviated from that given by the effective mass theory. In this report, wurtzite GaN is grown by the molecular-beam epitaxy method, and the photoluminescence spectra were measured. Three main peaks are assigned to the donor-acceptor transitions from two kinds of accepters. Some of the transitions were identified as coming from the cubic phase of GaN, which appears randomly within the predominantly hexagonal material. [S0163-1829(99)15915-0].

The structures and transformations on Si(113) surface

The surface structures of the Si(113)-(1 X 1), Si(113)-(3 X 1) and Si(113)-(3 X 2) have been studied theoretically by means of an ab initio quantum chemical CNDO method. We address not only the importance of the surface energy but also the energy minimization and the barrier height in the different structural conversion. We found that (1) the relaxed Si(113)-(1 X 1) structure. (2) the Si(113)-(3 X 1) close to the Si(113) Ranke (3 X 1)-2 model; (3) the atomic positions of Si(113)-(3 X 2) corrugated arrangement. (C) 1997 Elsevier Science B.V.

Effects of residual C and O impurities on photoluminescence in undoped GaN epilayers

Photoluminescence (PL) was investigated in undoped GaN from 4.8 K to room temperature. The 4.8 K spectra exhibited recombinations of free exciton, donor-acceptor pair (DAP), blue and yellow bands (Ybs). The blue band (BB) was also identified to be a DAP recombination. The YB was assigned to a recombination from deep levels. The energy-dispersive X-ray spectroscopy show that C and O are the main residual impurities in undoped GaN and that C concentration is lower in the epilayers with the stronger BB. The electronic structures of native defects, C and O impurities, and their complexes were calculated using ab initio local-density-functional (LDF) methods with linear muffin-tin-orbital and 72-atomic supercell. The theoretical analyses suggest that the electron transitions from O-N states to C-N and to V-Ga states are responsible for DAP and the BB, respectively, and the electron transitions between the inner levels of the C-N-O-N complex may be responsible for the YB in our samples. (C) 2002 Elsevier Science B.V. All rights reserved.

Theoretical study on analytical potential function and spectroscopic parameters for CaF molecule

The equilibrium properties and potential energy curves of the ground electronic state of CaF have been calculated using the Brueckner Doubles calculation with a triples contribution added [BD(T)] and the gradient-corrected density functional theory with three-parameter exact exchange mixing (B3LY-P) method, with 6-311 + G*,6-311 + G(2df,2pd) and 6-311 + G(3df,3pd) basis sets. All the computational PECs are fitted to analytical potential energy functions using Murrell-Sorbie, Huxley and Tang-Toennies potentials. Based on this, the spectroscopic parameters are calculated, and then compared with some other theoretical and experimental data. (C) 2004 Elsevier B.V. All rights reserved.

Golden-rule treatment on the ClO/ClO+ electron-transfer system

The structures, properties and electron transfer reactivity of the ClO/ClO+ coupling system are studied in this paper at ab initio (HF and MP2) levels and the density functional theory (DFT: B3LYP, B3P86, B3PW91) levels employing 6311 + G(3df) basis set and on the basis of the golden-rule of the time-dependent perturbation theory. Investigations indicate that the results got from the B3LYP method employing 6-311 + G(3df) basis set is in excellent agreement with the experiment. The activation energies, the stabilization energies and the electronic coupling matrix elements have also been calculated by using the B3LYP/6-311 + G(3df) method, and then the electron transfer rates are determined at this level. The electronic coupling matrix element of EC.6 is very small, only 0.03 kcal/mol, while that of EC.7 is the biggest, being 12.41 kcal/mol, the corresponding electron transfer rate is also the fastest among these seven encounter complexes. The averaged electron transfer rate is about 1.672 X 10(11) M-1 s(-1). It is indicated that the structures optimized by B3LYP method are more reliable than the results got from the other four methods. It also testified that the electronic coupling matrix element is the vital factor that significantly affects the electron transfer rate. (C) 2003 Elsevier B.V. All rights reserved.

Theoretical study on the ClO/ClO- system electron-transfer reactivity by the golden-rule

The structures, properties and electron transfer reactivity of the ClO/ClO- coupling system are studied in this paper at ab initio (UHF and UMP2) levels and the Density Functional Theory (DFT: UB3LYP, UB3P86, UB3PW91) levels employing 6311 + G(3df) basis set and on the basis of the Golden-rule of the time-dependent perturbation theory. Investigations indicate that the results obtained using the UB3LYP method employing 6-311 + G(3df) basis set is in excellent agreement with the experiment. For this coupling system, six stable coupling modes have been found which correspond to six different encounter complexes and denote six different electron transfer mechanism: four O-O directly linked structures (one collinear: D-h, one anti-parallel: C-s, two twist: C-2) and two Cl-O linked structures (cis- and anti- C-s structures). The activation energies, the stabilization energies and the electronic coupling matrix elements have also been calculated for the electron transfer reactions via these six different mechanism at the UB3LYP/6-311 + G(3df) level, and then the electron transfer rates are determined at the same level. The most favorable coupling mode to the electron transfer is the anti-parallel mechanism. The averaged electron transfer rate is about 5.58 X 10(11) M-1 s(-1). It is also implied that the B3LYP method can give more reasonable results for the electron transfer reactivity of this system. (C) 2003 Elsevier B.V. All rights reserved.

The properties and possible transformation path for C12B24N24

The structure and frequencies of C12B24N24 have been calculated by means of an ab initio method. By comparing the average bond energies with C-60, the calculated results predict that the cage C12B24N24 is a stable molecule. The calculated results indicate that the cage molecule C12B24N24 has a relative large HOMO-LUMO energy gap and a low rigidity The structures and stability of six possible isomers of C2B4N4 are used to suggest a possible transformation path from the pentagon CB2N2 to the C12B24N24 materials. (C) 2001 John Wiley & Sons, Inc.