Molecular dynamics simulation of hydrogen enhancing dislocation emission

The interactive pair potential between Al and H is obtained based on the ab initio calculation and the Chen-Mobius 3D lattice inversion formula. By utilizing the pair potentials calculated, the effects of hydrogen on the dislocation emission from crack tip have been studied. The simulated result shows that hydrogen can reduce the cohesive strength for Al single crystal, and then the critical stress intensity factor for partial dislocation emission decreases from 0.11 MPa root m (C-H = 0) to 0.075 MPa root m (C-H=0.72%) and 0.06 MPa root m (C-H = 1.44%). This indicates thar hydrogen can enhance the dislocation emission. The simulation also shows that atoms of hydrogen can gather and turn into small bubbles, resulting in enhancement of the equilibrium vacancy concentration.

Theoretical investigation of epitaxial deformation and the hcp-bcc transition of alkali metals

The deformation of alkali metals K, Rb, and Cs under epitaxial deformation is studied via the ab initio pseudopotential plane wave method using the local-density approximation. Under loading from the stable fee phase, metastable stares along directions [001], [111], and [201] are identified. One metastable state, presented at direction [201], has a very low symmetry in contrast to the planes [001] and [201]. Our results show that the softening direction and sequences of growth is significantly affected by the existence of the metastable states and magnitude of the energy barrier. The resulting softening sequences from soft to hard are [201], [110], [001], and [111] under biaxial compression and [001], [111], [201], and [110] under biaxial tension. An orthorhombic deformation path is used to investigate the fact, that the structure of the alkali films K and Cs evolve from the quasihexagonal structure into the (110)-oriented bcc structure, observed by experiments.

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.