Comparative studies of cation doping of ZnO with Mn, Fe, and Co.

MSINDO calculations were performed to elucidate the effect of doping and defects on the electronic properties of zinc oxide. The cyclic cluster Zn(48)O(48) served as a model for the bulk. Band gaps and stabilities of spin states were determined for doping with Mn, Fe, and Co. It was found that a substantial lowering of the band gap occurs for Mn and Co doping. In contrast, the Fe doping has a rather insignificant effect on the band gap. Additional defects such as oxygen vacancy, zinc interstitial, or zinc vacancy were also studied. In the case of substitution of zinc by two transition metal atoms, various spin states arise which can be classified as antiferromagnetic, ferrimagnetic, or ferromagnetic. We find that the spin state with the lowest multiplicity, antiferromagentic or ferrimagnetic, is more stable than the high-multiplicity ferromagnetic state in most considered cases, but additional zinc vacancies or mixed doping may reverse this trend.

Unusual optical properties of Mn-doped ZnO: the search for a new red pigment-a combined experimental and theoretical study.

A pigment of your imagination: A range of polycrystalline solid solutions of a zinc-rich Zn(x-1)Mn(x)O system (see figure) have been prepared and studied in terms of their colour, diffuse reflectance spectra, Mn valence state and electronic structure. The intense optical absorption arises from Mn(2+) doping and is thought to be due to forbidden or partially forbidden transitions between the valence and the conduction band.We report an investigation of zinc-rich polycrystalline solid solutions of the Zn(1-x)Mn(x)O system concerning the colour, the diffuse reflectance spectra, the valence state of manganese and the electronic structure. Samples were prepared by a chemical-vapour-transport-assisted route and optimized with respect to colour strength. In agreement with previous experimental results, EPR studies showed that manganese is in the divalent charge state. The nature of the very intense optical absorption, which is caused by Mn(2+) doping and determines the colour of the material, is discussed. It is argued that the Mn(2+)-induced optical absorption is due to forbidden or partially forbidden transitions between the valence and the conduction band that involve Mn admixed states. This assignment is also confirmed by quantum chemical calculations using the semiempirical molecular orbital method MSINDO.

Influence of intrinsic defects on the properties of zinc oxide.

The effects of oxygen vacancies and zinc interstitials on the structure and energy of zinc oxide were studied with the semiempirical MO method MSINDO. Cyclic clusters were chosen as model systems. Single and multiple removal of oxygen atoms and zinc interstitials in zinc oxide served to determine the defect formation energy and the band gap. The interaction between two and three oxygen vacancies was investigated. The vacancies cause a decrease of the band gap, which originates from an occupied defect level. This is also found for zinc interstitials under zinc rich conditions. The defect formation energy of such zinc interstitials is found to be lower than that of oxygen vacancies at 0 K but decreases for oxygen vacancies and increases for zinc interstitials with increasing temperature.

Anion substitution in zinc chalcogenides.

Anion substitution effects on the structure and energy of zinc chalcogenides were studied with the semiempirical molecular orbital method MSINDO. Cyclic clusters of different sizes were chosen as model systems. The convergence of the bulk properties of the perfect clusters with increasing cluster size was tested. Single and multiple substitution of oxygen atoms in zinc oxide by sulfur and of sulfur atoms in zinc sulfide by oxygen served to determine the energetics of substitution for these two cases. It was found that the substitution of oxygen by sulfur in ZnO is easier than the substitution of sulfur by oxygen in ZnS in agreement with experimental results. The interaction between two oxygen atoms vs. two selenium atoms in zinc sulfide was investigated. Oscillations of the cluster energy in dependence of the distance between the two doping atoms were observed. These are explained by the relative sites of the doping atoms in the crystal lattice. The magnitude of the oscillations is smaller in ZnS:Se than in ZnS:O, because the difference between the anion radii of S2- and Se2- is smaller than between S2- and O2-. This is also reflected in the band gap.

Molecular dynamics study of the thermal entropy in mixed zinc chalcogenides.

Molecular dynamics (MD) calculations were performed to determine the vibrational contribution to the entropy of mixing and its importance for the mixing of ZnO/ZnS and ZnS/Zn3P2. These systems were modeled by cyclic clusters Zn48O48, Zn32S32, and Zn48P32. The mixed cyclic clusters considered were Zn48O47S, Zn32S31O, Zn33S30P2, and Zn47S2P30. For each of the clusters, the entropy was calculated in the range of the experimental temperature of the mixing process. The convergence of the entropy with respect to the number of MD steps was studied. Finally, the thermal part of the entropy of mixing was determined, and its dependence on the number of MD steps was investigated. It was found that the thermal entropy is important for the Gibbs free energy of mixing near the miscibility gaps.

Molecular dynamics investigation of oxygen vacancy diffusion in rutile.

Oxygen vacancy diffusion in rutile was studied by Born-Oppenheimer molecular dynamics techniques in the framework of the semiempirical molecular orbital method MSINDO. Migration of an oxygen vacancy from the rutile (110) surface towards the bulk was simulated. The metadynamics technique was employed to accelerate the diffusion processes. In this way, transition state structures and activation energies for the diffusion processes were obtained. Rate constants and the time scale of diffusion processes were estimated for different temperatures using the calculated activation energy. It was found that the vacancies in the bulk are less stable than on the surface. The feasibility of oxygen vacancy diffusion under experimental conditions is discussed.

Toward an understanding of the formation of vanadia-titania catalysts.

Structures of hydrated vanadia species on the TiO2-anatase surfaces were investigated using the semiempirical molecular orbital method MSINDO. The (101), (001), and (100) surfaces of anatase were considered. They were modeled by appropriate two-dimensional cyclic clusters of TiO2. Monomeric and dimeric hydrated vanadia species on the anatase surfaces were simulated by adsorbing VO4H3 and V2O7H4 molecules, respectively. Different adsorption structures were considered, and their stabilities at 300 and 600 K were tested by constant-temperature Born-Oppenheimer molecular dynamics simulations in the framework of MSINDO. Structural features of the vanadia-titania catalysts found in extended X-ray absorption fine structure, secondary ion mass spectrometry, IR, Raman, and NMR spectroscopy and conductivity experiments can be explained by the present calculations.

Miscibility of zinc sulfide and zinc phosphide.

Solid solutions in the system zinc sulfide/zinc phosphide (Zn(2+)(x)S(2-2xP(2x)) were investigated using the cyclic cluster model within the semiempirical MSINDO method. Results of cyclic cluster calculations for binding energies of the perfect ZnS and Zn(3)P(2) are presented and compared with the experimental data. The miscibility of ZnS and Zn(3)P(2) over the whole composition range of 0 < x < 1 was investigated by calculating the Gibbs free energy of mixing Delta(M)G for different values of x. A miscibility gap was found at both ends of the composition range and compared with experimental data.

Models for the treatment of crystalline solids and surfaces.

Crystalline solids and surfaces have become a subject of growing interest. The difficulty of a comprehensive description of a variety of phenomena by a single method has led to the development of many models. These models can be classified as nonperiodic and periodic models. The former include free clusters, saturated clusters, and embedded clusters. The latter two models serve to remove the boundary effects of the free clusters. No perfect avoidance of such effects can be achieved in this way. The cyclic cluster model overcomes this difficulty in a natural way. It is periodic with a finite periodicity. An embedding can take into account a long-range effect in ionic crystals. Previous periodic approaches relied on the large unit cell model, which is related to the supercell approach. For perfect crystals the conventional unit cell approach is a well-known standard. However, its disadvantage is the unphysical periodicity of defects, which is avoided in the cyclic cluster model. The present article presents a description of these models together with selective applications to solid-state systems and surfaces.

Molecular dynamics implementation in MSINDO: study of silicon clusters.

Born-Oppenheimer molecular dynamics is implemented in the semiempirical self-consistent field molecular orbital method MSINDO. The method is employed for the investigation of the structure and dynamics of silicon clusters of various sizes. The reliability of the present parameterization for silicon compounds is demonstrated by a comparison of the results of simulated annealing and of density functional calculations of Si(n) clusters (n = 5-7). The melting behavior of the Si(7) cluster is investigated and the MSINDO results are compared to previous high-level calculations. The efficiency of the present approach for the treatment of large systems is demonstrated by an extensive simulated annealing study of the Si(45) and Si(60) clusters. New Si(45) and Si(60) structures are found and evaluated. The relative stability of various energy minimum structures is compared with density functional calculations and available literature data.

Bond energies for molecules, clusters, and deposit systems.

A new approach is suggested to the assignment of bond energies in molecules and clusters. It uses a shareholder principle for the redistribution of the shifts in atomic energies, which arise in a molecule, on the bonds. The scheme is directly suitable for semiempirical methods, where only one- and two-center terms occur. MSINDO calculations are performed to demonstrate the suitability of the approach for molecules and clusters. As an application the bonding in a deposit system is analyzed for the case of copper on magnesium oxide. It is found that copper atoms do not only bind to the preferred oxygen sites but also substantially to the magnesium sites. The copper-copper bonds are the strongest and will determine the structure of copper clusters on magnesium oxide surfaces.

Mechanism for large first hyperpolarizabilities of phosphonic acid stilbene derivatives.

This paper presents calculations of dipole moments (mu), static polarizabilities (alpha), and first hyperpolarizabilities (beta) of phosphonic acid stilbene derivatives calculated in the framework of density functional theory. These calculations were performed using a finite field approach implemented in the density functional program ALLCHEM and were of an all-electron type using local exchange-correlation functional and specially designed basis sets. The molecular structures have been fully optimized using the semiempirical program MSINDO. Some of the investigated stilbenes have been synthesized very recently while others are described for the first time. Donor and acceptor groups of these analogues have been modified and the influence of these changes on the first hyperpolarizabilities has been investigated. This work demonstrates that the nonlinear optical response beta of these compounds increases dramatically when the acceptor moiety is displaced by analogues containing alkali metal groups. A general mechanism for the design of novel nonlinear optical materials with large first hyperpolarizabilities is described.

On the Alternation Effect in Substituted Indolizines and Their Aza-analogs.

The SINDO1 method is used to calculate charge distribution and total energies in the series of isomeric nitro- and methoxy-substituted indolizines, azaindolizines, and N-methylazaindolizinium cations. The influence of donor and acceptor substituents on the total energy and charges may be treated according to the simple alternation rule. In the series of isomers (substituted or azaindolizines) the lower energies are observed for those molecules, where acceptor substituents are arranged consonant to the polarity of the chain N-C(5)-C(6)-C(7)-C(8)-C(9)-C(1)-C(2)-C(3) of the indolizine. The results are used to treat known peculiarities of ring opening/transformation of indolizine and its aza-derivatives. It is proved that the energies of C(3) and C(5) adducts of isomeric azaindolizines with hydroxyl anion follow the same alternation pattern.