Using first principle electronic structure calculations within density functional theory and the supercell model, we have investigated the nature and formation energies of defect states associated with Ga and Te vacancies and Ge and Sn substitutional impurities in GaTe. We have also calculated the band structure of pure GaTe for comparison with systems with defects and also to find out the importance of spin–orbit interaction (SOI) on its band structure. We find that the top valence band at the Γ-point shifts up in energy by 0.1 eV due to the mixing of Te p x –p y and p z bands, this splitting being considerably smaller than in atoms where it is 0.8 eV. From an analysis of charge densities and band structures associated with the defect states, we find that most of them are strongly localized and lie deep in the band gap region. The calculated binding energy of the deep defect state and the ε(−1 / −2) transition level associated with the Ga vacancy appears to be in good agreement with experiment. Formation energy calculations suggest that V Ga is the preferred intrinsic defect in GaTe.

VL - 21 ER - TY - JOUR T1 - Thermochemical Kinetics for Multi-reference Systems: Addition Reactions of Ozone JF - Journal of Physical Chemistry A Y1 - 2009 A1 - Piotr Piecuch A1 - Zhao, Y. A1 - Truhlar, D.G. AB -The 1,3-dipolar cycloadditions of ozone to ethyne and ethene provide extreme examples of multireference singlet-state chemistry, and they are examined here to test the applicability of several approaches to thermochemical kinetics of systems with large static correlation. Four different multireference diagnostics are applied to measure the multireference characters of the reactants, products, and transition states; all diagnostics indicate significant multireference character in the reactant portion of the potential energy surfaces. We make a more complete estimation of the effect of quadruple excitations than was previously available, and we use this with CCSDT/CBS estimation of Wheeler et al. (Wheeler, S. E.; Ess, D. H.; Houk, K. N. J. Phys. Chem. A 2008, 112, 1798.) to make new best estimates of the van der Waals association energy, the barrier height, and the reaction energy to form the cycloadduct for both reactions. Comparing with these best estimates, we present comprehensive mean unsigned errors for a variety of coupled cluster, multilevel, and density functional methods. Several computational aspects of multireference reactions are considered: (i) the applicability of multilevel theory, (ii) the convergence of coupled cluster theory for reaction barrier heights, (iii) the applicability of completely renormalized coupled cluster methods to multireference systems, (iv) the treatment by density functional theory, (v) the multireference perturbation theory for multireference reactions, and (vi) the relative accuracy of scaling-type multilevel methods as compared with additive ones. It is found that scaling-type multilevel methods do not perform better than the additive-type multilevel methods. Among the 48 tested density functionals, only M05 reproduces the best estimates within their uncertainty. Multireference perturbation theory based on the complete-active-space reference wave functions constructed using a small number of reaction-specific active orbitals gives accurate forward barrier heights; however, it significantly underestimates reaction energies.

VL - 113 IS - 19 ER - TY - JOUR T1 - Three Modes of Metal-Enriched Star Formation in the Early Universe JF - \apj Y1 - 2009 A1 - Smith, B. D. A1 - Turk, M. J. A1 - Sigurdsson, S. A1 - O'Shea, B. W. A1 - Norman, M. L. KW - early universe KW - stars: formation VL - 691 ER - TY - JOUR T1 - Turbulence-Interface Interactions in a Two-Fluid Homogeneous Flow JF - Physics of Fluids Y1 - 2009 A1 - Li, Z. A1 - Jaberi, F.A. AB -The two-way interactions between the turbulent velocity field and the interface in an incompressible two-fluid homogeneous turbulent flow are studied with a recently developed Lagrangian–Eulerian interfacial particle level-set method. The numerical results confirm that the rate of change of the interface area is directly related to the work done by the surface tension force. While the surface tension damps the surrounding turbulence in the “interface stretching period” to oppose the increase in interface area, it is shown to actually increase the turbulent kinetic energy when the interface experiences compression. Additionally, the surface tension force is found to generate strong vortical motions close to the interface through the baroclinic torque effects. There is also an increase in strain rate and the viscous dissipation rate of turbulent kinetic energy in the interface region. The effect of interface on the surrounding turbulence appears primarily in the direction perpendicular to the interface. Analysis of the vorticity and kinetic energy equations indicates that the turbulence-interface interactions are strongly dependent on the fluids’ density ratio and the Weber number.

VL - 21 IS - 9 ER - TY - JOUR T1 - Theoretical studies of defect states in GaTe JF - Journal of Physics: Condensed Matter Y1 - 2008 A1 - Rak, Zs A1 - S. D. Mahanti A1 - Mandal, Krishna C A1 - Fernelius, N.C. VL - 21 IS - 1 ER - TY - JOUR T1 - Transport Model Simulations of Projectile Fragmentation Reactions at 140 {MeV/nucleon} JF - Physical Review C Y1 - 2008 A1 - M. Mocko A1 - M. B. Tsang A1 - D. Lacroix A1 - A. Ono A1 - P. Danielewicz A1 - W. G. Lynch A1 - R. J Charity KW - Nuclear Experiment KW - Nuclear Theory AB -The collisions in four different reaction systems using {\$ˆ{40,48}\$Ca} and {\$ˆ{58,64}\$Ni} isotope beams and a Be target have been simulated using the Heavy Ion Phase Space Exploration and the Antisymmetrized Molecular Dynamics models. The present study mainly focuses on the model predictions for the excitation energies of the hot fragments and the cross sections of the final fragments produced in these reactions. The effects of various factors influencing the final fragment cross sections, such as the choice of the statistical decay code and its parameters have been explored. The predicted fragment cross sections are compared to the projectile fragmentation cross sections measured with the A1900 mass separator. At {\$E/A=140\$} {MeV,} reaction dynamics can significantly modify the detection efficiencies for the fragments and make them different from the efficiencies applied to the measured data reported in the previous work. The effects of efficiency corrections on the validation of event generator codes are discussed in the context of the two models.

VL - 78 IS - 2 N1 - {Phys.Rev.C78:024612,2008} ER - TY - JOUR T1 - Theoretical Characterization of End-on and Side-on Peroxide Coordination in Ligated Cu2O2 Models JF - The Journal of Physical Chemistry A Y1 - 2006 A1 - C.J. Cramer A1 - A. Kinal A1 - M. Wloch A1 - Piotr Piecuch A1 - L. Gagliardi AB -The relative energetics of μ-η1:η1 (trans end-on) and μ-η2:η2 (side-on) peroxo isomers of Cu2O2 fragments supported by 0, 2, 4, and 6 ammonia ligands have been computed with various density functional, coupled-cluster, and multiconfigurational protocols. There is substantial disagreement between the different levels for most cases, although completely renormalized coupled-cluster methods appear to offer the most reliable predictions. The significant biradical character of the end-on peroxo isomer proves problematic for the density functionals, while the demands on active space size and the need to account for interactions between different states in second-order perturbation theory prove challenging for the multireference treatments. In the latter case, it proved impossible to achieve any convincing convergence.

VL - 110 IS - 40 ER - TY - JOUR T1 - Theoretical Models on the Cu2O2 Torture Track. Mechanistic Implications for Oxytyrosinase and Small-molecule Analogs JF - The Journal of Physical Chemistry A Y1 - 2006 A1 - C.J. Cramer A1 - M. Wloch A1 - Piotr Piecuch A1 - Cristina Puzzarini A1 - L. Gagliardi AB -Accurately describing the relative energetics of alternative bis(mu-oxo) and mu-eta2:eta2 peroxo isomers of Cu2O2 cores supported by 0, 2, 4, and 6 ammonia ligands is remarkably challenging for a wide variety of theoretical models, primarily owing to the difficulty of maintaining a balanced description of rapidly changing dynamical and nondynamical electron correlation effects and a varying degree of biradical character along the isomerization coordinate. The completely renormalized coupled-cluster level of theory including triple excitations and extremely efficient pure density functional levels of theory quantitatively agree with one another and also agree qualitatively with experimental results for Cu2O2 cores supported by analogous but larger ligands. Standard coupled-cluster methods, such as CCSD(T), are in most cases considerably less accurate and exhibit poor convergence in predicted relative energies. Hybrid density functionals significantly underestimate the stability of the bis(mu-oxo) form, with the magnitude of the error being directly proportional to the percentage Hartree-Fock exchange in the functional. Single-root CASPT2 multireference second-order perturbation theory, by contrast, significantly overestimates the stability of bis(mu-oxo) isomers. Implications of these results for modeling the mechanism of C-H bond activation by supported Cu2O2 cores, like that found in the active site of oxytyrosinase, are discussed.

VL - 110 ER - TY - JOUR T1 - Two New Classes of Non-Iterative Coupled-Cluster Methods Derived from the Method of Moments of Coupled-Cluster Equations JF - Molecular Physics Y1 - 2006 A1 - Piotr Piecuch A1 - M. Wloch A1 - M.D. Logriguito A1 - J. R Gour VL - 104 IS - 13 & 14 ER -