Design of shallow acceptors in ZnO: First-principles band-structure calculations

p-type doping is a great challenge for the full utilization of ZnO as short-wavelength optoelectronic material. Due to a large electronegative characteristic of oxygen, the ionization energy of acceptors in ZnO is usually too high. By analyzing the defect wave-function character, we propose several approaches to lower the acceptor ionization energy by codoping acceptors with donor or isovalent atoms. Using the first-principles band-structure method, we show that the acceptor transition energies of V-Zn-O-O can be reduced by introducing F-O next to V-Zn to reduce electronic potential, whereas the acceptor transition energy of N-O-nZn(Zn) (n=1-4) can be reduced if we replace Zn by isovalent Mg or Be to reduce the anion and cation kinetic p-d repulsion, as well as the electronic potential.

Metamorphosis of self-organized quantum dots into quantum wires

Spontaneous formation of InAs quantum wires in InAlAs on InP(001) via sequential chain-like coalescence of quantum dots along [1 (1) over bar 0] is realized. Theoretical calculations based on the energetics of interacting steps provide a qualitative explanation for the experimental results. Sequential coalescence of initially isolated dots reduces the total free energy strikingly. Thus the wire-like structure is energetically favorable. (C) 1998 Elsevier Science B.V.

INTERVALLEY-GAMMA-CHI DEFORMATION POTENTIALS IN III-V ZINCBLENDE AND SI SEMICONDUCTORS BY AB-INITIO PSEUDOPOTENTIAL CALCULATIONS

Intervalley GAMMA - X deformation potential constants (IVDP's) have been calculated by first principle pseudopotential method for the III-V zincblende semiconductors AlP, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs and InSb. As a prototype crystal we have also carried out calculations on Si. When comparing the calculated IVDP's of LA phonon for GaP, InP and InAs and LO phonon for AlAs, AlSb, GaAs, GaSb and InSb with a previous calculation by EPM in rigid approximation, good agreements are found. However, our ab initio pseudopotential results of LA phonon for AlAs, AlSb, GaAs, GaSb and InSb and LO phonon for GaP, InP and InAs are about one order of magnitude smaller than those obtained by EPM calculations, which indicate that the electron redistributions upon the phonon deformations may be important in affecting GAMMA - X intervalley shatterings for these phonon modes when the anions are being displaced. In our calculations the phonon modes of LA and LO at X point have been evaluated in frozen phonon approximation. We have obtained, at the same time, the LAX and LOX phonon frequencies for these materials from total energy calculations. The calculated phonon frequencies agree very well with experimental values for these semiconductors.

INTERVALLEY GAMMA-CHI DEFORMATION POTENTIALS IN III-V ZINCBLENDE SEMICONDUCTORS BY ABINITIO PSEUDOPOTENTIAL CALCULATIONS

Intervalley GAMMA-X deformation-potential constants (IVDP's) have been calculated by use of a first-principles pseudopotential method for the III-V zinc-blende semiconductors AlP, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs, and InSb. When the calculated IVDP's of LA phonons for GaP, InP, and InAs and of LO phonons for AlAs, AlSb, GaAs, GaSb, and InSb are compared with results of a previous calculation that used the empirical pseudopotential method (EPM) and a rigid-ion approximation, good agreement is found. However, our ab initio pseudopotential results on IVDP's of LA phonons for AlAs, AlSb, GaAs, GaSb, and InSb and of LO phonons for GaP, InP, and InAs are about one order of magnitude smaller than those obtained by use of EPM calculations, indicating that the electron redistribution accompanying crystal-lattice deformation has a significant effect on GAMMA-X intervalley scattering for these phonon modes when the anions are being displaced. In our calculations the LA- and LO-phonon modes at the X point have been evaluated in the frozen-phonon approximation. We have also obtained the LAX- and LOX-phonon frequencies for these materials from total-energy calculations, which agree very well with experimental values for these semiconductors. We have also calculated GAMMA-X hole-phonon scattering matrix elements for the top valence bands in these nine semiconductors, from which the GAMMA-X IVDP's of the top valence bands for the longitudinal phonons and transverse phonons are evaluated, respectively.

CALCULATION OF SPIN-REORIENTATION TEMPERATURE IN ND2FE14B

The spin-reorientation phenomenon in Nd2Fe14B has been investigated using an angular dependent free energy approach. A magnetic Hamiltonian which includes the crystal electric field term and the exchange term has been established using realistic band structure results. The temperature dependence of the molecular field is accounted for by introducing the Brillouin function and the magnetic Hamiltonian is diagonalized within the ground state multiplet of the Nd ion. The eigenstates are then used to form the partition function for the free energy. At each temperature, the direction of the molecular field is obtained by searching for the minimum in the angular parameter space of the free energy. Our calculations show that for Nd2Fe14B, the net magnetic anisotropy direction is canted away from the c axis at a temperature close to the experimentally reported spin-reorientation temperature of 150 K. The temperature dependence of the magnetic structure is found to be very sensitive to the size of the second order crystal field parameter B20.

1ST PRINCIPLES PSEUDOPOTENTIAL CALCULATIONS OF ZONE BOUNDARY PHONON FREQUENCIES OF III-V-ZINCBLENDE SEMICONDUCTORS

Longitudinal zone boundary X phonon frequencies have been calculated by a first principles pseudopotential method for III-V zincblende semiconductors AlP, AlAs, AlSb, GaP, GaAs, GaSb, InP, InAs and InSb. The phonon frequencies have been evaluated from total energy calculations in the frozen phonon approximation. The calculated phonon frequencies agree very well with the experimental values.

Ab initio calculations for the neutral and charged O vacancy in sapphire

The energetics, lattice relaxation, and the defect-induced states of st single O vacancy in alpha-Al2O3 are studied by means of supercell total-energy calculations using a first-principles method based on density-functional theory. The supercell model with 120 atoms in a hexagonal lattice is sufficiently large to give realistic results for an isolated single vacancy (square). Self-consistent calculations are performed for each assumed configuration of lattice relaxation involving the nearest-neighbor Al atoms and the next-nearest-neighbor O atoms of the vacancy site. Total-energy data thus accumulated are used to construct an energy hypersurface. A theoretical zero-temperature vacancy formation energy of 5.83 eV is obtained. Our results show a large relaxation of Al (O) atoms away from the vacancy site by about 16% (8%) of the original Al-square (O-square) distances. The relaxation of the neighboring Al atoms has a much weaker energy dependence than the O atoms. The O vacancy introduces a deep and doubly occupied defect level, or an F center in the gap, and three unoccupied defect levels near the conduction band edge, the positions of the latter are sensitive to the degree of relaxation. The defect state wave functions are found to be not so localized, but extend up to the boundary of the supercell. Defect-induced levels are also found in the valence-band region below the O 2s and the O 2p bands. Also investigated is the case of a singly occupied defect level (an F+ center). This is done by reducing both the total number of electrons in the supercell and the background positive charge by one electron in the self-consistent electronic structure calculations. The optical transitions between the occupied and excited states of the: F and F+ centers are also investigated and found to be anisotropic in agreement with optical data.

THE ISOSPIN DEPENDENCE OF THE NUCLEAR PHASE TRANSITION NEAR THE CRITICAL POINT

The experimental results reveal the isospin dependence of the nuclear phase transition in terms of the Landau Free Energy description of critical phenomena. Near the critical point, different ratios of the neutron to proton concentrations lead to different critical points for the phase transition which is analogous to the phase transitions in He-4-He-3 liquid mixtures. The antisymmetrized molecular dynamics (AMD) and GEMINI models calculations were also performed and the results will be discussed as well.

Computational study on the addition of HCN to methanimine catalyzed by formamidine and formamide

In this paper, we have carried out a theoretical study on the addition of HCN to methanimine with formamidine or formamide using second-order Moller-Plesset perturbation (MP2) method with 6-31 + G(d,p) basis sets. At MP2 level. a high-energy, intermediate has been located for each pathway. The addition of HCN to methanimine with formamidine has the lowest free energy barrier according to the calculations at MP2 level. (C) 2004 Elsevier B.V. All rights reserved.

Crystallization behavior of poly(epsilon-caprolactone) in poly(epsilon-caprolactone) and poly(vinyl methyl ether) mixtures

The morphological development and crystallization behavior of poly(epsilon-caprolactone) (PCL) in miscible mixtures of PCL and poly(vinyl methyl ether) (PVME) were investigated by optical microscopy as a function of the mixture composition and crystallization temperature. The results indicated that the degree of crystallinity of PCL was independent of the mixture composition upon melt crystallization because the glass-transition temperatures of the mixtures were much lower than the crystallization temperature of PCL. The radii of the PCL spherulites increased linearly with time at crystallization temperatures ranging from 42 to 49 degrees C. The isothermal growth rates of PCL spherulites decreased with the amount of the amorphous PVME components in the mixtures. Accounting for the miscibility of PCL/PVME mixtures, the radial growth rates of PCL spherulites were well described by a kinetic equation involving the Flory-Huggins interaction parameter and the free energy for the nuclei formation in such a way that the theoretical calculations were in good agreement with the experimental data. From the analysis of the equilibrium melting point depression, the interaction energy density of the PVME/PCL system was calculated to be -3.95 J/cm(3).

Study on the origin of inverted phase in drying solution-cast block copolymer films

In a previous study, we reported observation of the novel inverted phase (the minority blocks comprising the continuum phase) in kinetically controlled phase separating solution-cast poly(styrene-b-butadiene-b-styrene) (SBS) triblock copolymer films [Zhang et al. Macromolecules 2000, 33, 9561-7]. In this study, we adopt the same approach to investigate the formation of inverted phase in a series of solution-cast poly(styrene-b-butadiene) (SB) asymmetric diblock copolymers having nearly equal polystyrene (PS) weight fraction (about 30 wt %) but different molecular weights. The microstructure of the solution-cast block copolymer films resulting from different solvent evaporation rates, R, was inspected, from which the kinetically frozen-in phase structures at qualitatively different block copolymer concentrations and correspondingly different effective interaction parameter, chieff, can be deduced. Our result shows that there is a threshold molecular weight or range of molecular weight below which the unusual inverted phase is accessible by controlling the solvent evaporation rate. In comparing the present result with that of our previous study on the SBS triblock copolymer, we find that the formation of the inverted phase has little bearing on the chain architecture. We performed numerical calculations for the free energy of block copolymer cylinders and found that the normal phase is always preferred irrespective of the interaction parameter and molecular weight, which suggests the formation of the inverted phase to have a kinetic origin.

Docking and molecular dynamics studies toward the binding of new natural phenolic marine inhibitors and aldose reductase

Phenolic marine natural product is a kind of new potential aldose reductase inhibitors (ARIs). In order to investigate the binding mode and inhibition mechanism, molecular docking and dynamics studies were performed to explore the interactions of six phenolic inhibitors with human aldose reductase (hALR2). Considering physiological environment, all the neutral and other two ionized states of each phenolic inhibitor were adopted in the simulation. The calculations indicate that all the inhibitors are able to form stable hydrogen bonds with the hALR2 active pocket which is mainly constructed by residues TYR48, HIS110 and TRP111, and they impose the inhibition effect by occupying the active space. In all inhibitors, only La and its two ionized derivatives La_ion1 and La_ion2, in which neither of the ortho-hydrogens of 3-hydroxyl is substituted by Br, bind with hALR2 active residues using the terminal 3-hydroxyl. While, all the other inhibitors, at least one of whose ortho-sites of 3- and 6-hydroxyls are substituted by Br substituent which take much electron-withdrawing effect and steric hindrance, bind with hALR2 through the lactone group. This means that the Br substituent can effectively regulate the binding modes of phenolic inhibitors. Although the lactone bound inhibitors have relatively high RMSD values, our dynamics study shows that both binding modes are of high stability. For each inhibitor molecule, the ionization does not change its original binding mode, but it does gradually increase the binding free energy, which reveals that besides hydrogen bonds, the electrostatic effect is also important to the inhibitor–hALR2 interaction.

Theory of the near K-edge structure in electron energy loss spectroscopy

Arguments are given that lead to a formalism for calculating near K-edge structure in electron energy loss spectroscopy (EELS). This is essentially a one electron picture, while many body effects may be introduced at different levels, such as the local density approximation to density functional theory or the GW approximation to the electron self-energy. Calculations are made within the all electron LMTO scheme in crystals with complex atomic and electronic structures, and these are compared with experiment. (c) 2004 Elsevier B.V. All rights reserved.

Reaction Mechanisms of Crotonaldehyde Hydrogenation on Pt(111): Density Functional Theory and Microkinetic Modeling

The microkinetics based on density function theory (DFT) calculations is utilized to investigate the reaction mechanism of crotonaldehyde hydrogenation on Pt(111) in the free energy landscape. The dominant reaction channel of each hydrogenation product is identified. Each of them begins with the first surface hydrogenation of the carbonyl oxygen of crotonaldehyde on the surface. A new mechanism, 1,4-addition mechanism generating enols (butenol), which readily tautomerize to saturated aldehydes (butanal), is identified as a primary mechanism to yield saturated aldehydes instead of the 3,4-addition via direct hydrogenation of the ethylenic bond. The calculation results also show that the full hydrogenation product, butylalcohol, mainly stems from the deep hydrogenation of surface open-shell dihydrogenation intermediates. It is found that the apparent barriers of the dominant pathways to yield three final products are similar on P(111), which makes it difficult to achieve a high selectivity to the desired crotyl alcohol (COL).

Mechanistic Study of 1,3-Butadiene Formation in Acetylene Hydrogenation over the Pd-Based Catalysts Using Density Functional Calculations

Green oil, which leads to the deactivation of the catalysts used for the selective hydrogenation of acetylene, has long been observed but its formation mechanism is not fully understood. In this work, the formation of 1,3-butadiene, known to be the precursor of green oil, on both Pd(111) and Pd(211) surfaces is examined using density functional theory calculations. The pathways containing C-2 + C-2 coupling reactions as well as the corresponding hydrogenation reactions are studied in detail. Three pathways for 1,3-butadiene production, namely coupling plus hydrogenation and further hydrogenation, hydrogenation plus coupling plus hydrogenation, and a two step hydrogenation followed by coupling, are determined. By comparing the effective barriers, we identify the favored pathway on both surfaces. A general understanding toward the deactivation process of the industrial catalysts is also provided. In addition, the effects of the formation of subsurface carbon atoms as well as the Ag alloying on the 1,3-butadiene formation on Pd-based catalysts are also investigated and compared with experimental results.