Ground and Excited State Intramolecular Proton Transfer in Salicylic Acid: an Ab Initio Electronic Structure Investigation

Energetics of the ground and excited state intramolecular proton transfer in salicylic acid have been studied by ab initio molecular orbital calculations using the 6-31G** basis set at the restricted Hartree-Fock (RHF) and configuration interaction-single excitation (CIS) levels and also using the semiempirical method AM1 at the RHF level as well as with single and pair doubles excitation configuration interaction spanning eight frontier orbitals (PECI = 8). The ab initio potential energy profile for intramolecular proton transfer in the ground state reveals a single minimum corresponding to the primary form, in the first excited singlet state, however, there are two minima corresponding to the primary and tautomeric forms, separated by a barrier of similar to 6 kcal/mol, thus accounting for dual emission in salicylic acid. Electron density changes with electronic excitation and tautomerism indicate no zwitterion formation. Changes in spectral characteristics with change in pH, due to protonation and deprotonation of salicylic acid, are also accounted for, qualitatively. Although the AM1 calculations suggest a substantial barrier for proton transfer in the ground as well as the first excited state of SA, it predicts the transition wavelength in near quantitative accord with the experimental results for salicylic acid and its protonated and deprotonated forms.

Prediction of the formation and stabilities of energetic salts and ionic liquids based on ab initio electronic structure calculations

A computational approach to predict the thermodynamics for forming a variety of imidazolium-based salts and ionic liquids from typical starting materials is described. The gas-phase proton and methyl cation acidities of several protonating and methylating agents, as well as the proton and methyl cation affinities of many important methyl-, nitro-, and cyano- substituted imidazoles, have been calculated reliably by using the computationally feasible DFT (B3LYP) and MP2 (extrapolated to the complete basis set limit) methods. These accurately calculated proton and methyl cation affinities of neutrals and anions are used in conjunction with an empirical approach based on molecular volumes to estimate the lattice enthalpies and entropies of ionic liquids, organic solids, and organic liquids. These quantities were used to construct a thermodynamic cycle for salt formation to reliably predict the ability to synthesize a variety of salts including ones with potentially high energetic densities. An adjustment of the gas phase thermodynamic cycle to account for solid- and liquid-phase chemistries provides the best overall assessment of salt formation and stability. This has been applied to imidazoles (the cation to be formed) with alkyl, nitro, and cyano substituents. The proton and methyl cation donors studied were as follows: HCl, HBr, HI, (HO)(2)SO2, HSO3CF3 (TfOH), and HSO3(C6H4)CH3 (TsOH); CH3Cl, CH3Br, CH3I, (CH3O)(2)SO2, CH3SO3CF3 (TfOCH3) and CH3SO3(C6H4)CH3 (TsOCH3). As substitution of the cation with electron-withdrawing groups increases, the triflate reagents appear to be the best overall choice as protonating and methylating agents. Even stronger alkylating agents should be considered to enhance the chances of synthetic success. When using the enthalpies of reaction for the gas-phase reactants (eq 6) to form a salt, a cutoff value of - 13 kcal mol(-1) or lower (more negative) should be used as the minimum value for predicting whether a salt can be synthesized.

Effective t-J model Hamiltonian parameters of monolayered cuprate superconductors from ab initio electronic structure calculations

The magnetic coupling constant of selected cuprate superconductor parent compounds has been determined by means of embedded cluster model and periodic calculations carried out at the same level of theory. The agreement between both approaches validates the cluster model. This model is subsequently employed in state-of-the-art configuration interaction calculations aimed to obtain accurate values of the magnetic coupling constant and hopping integral for a series of superconducting cuprates. Likewise, a systematic study of the performance of different ab initio explicitly correlated wave function methods and of several density functional approaches is presented. The accurate determination of the parameters of the t-J Hamiltonian has several consequences. First, it suggests that the appearance of high-Tc superconductivity in existing monolayered cuprates occurs with J/t in the 0.20¿0.35 regime. Second, J/t=0.20 is predicted to be the threshold for the existence of superconductivity and, third, a simple and accurate relationship between the critical temperatures at optimum doping and these parameters is found. However, this quantitative electronic structure versus Tc relationship is only found when both J and t are obtained at the most accurate level of theory.

Emergence of half-metallicity in suspended NiO chains: Ab initio electronic structure and quantum transport calculations

Contrary to the antiferromagnetic and insulating character of bulk NiO, one-dimensional chains of this material can become half metallic due to the lower coordination of their atoms. Here we present ab initio electronic structure and quantum transport calculations of ideal infinitely long NiO chains and of more realistic short ones suspended between Ni electrodes. While infinite chains are insulating, short suspended chains are half-metallic minority-spin conductors that displays very large magnetoresistance and a spin-valve behavior controlled by a single atom.

Magnetic structure of Li2CuO2: From ab initio calculations to macroscopic simulations

The magnetic structure of the edge-sharing cuprate compound Li2CuO2 has been investigated with highly correlated ab initio electronic structure calculations. The first- and second-neighbor in-chain magnetic interactions are calculated to be 142 and -22 K, respectively. The ratio between the two parameters is smaller than suggested previously in the literature. The interchain interactions are antiferromagnetic in nature and of the order of a few K only. Monte Carlo simulations using the ab initio parameters to define the spin model Hamiltonian result in a Nel temperature in good agreement with experiment. Spin population analysis situates the magnetic moment on the copper and oxygen ions between the completely localized picture derived from experiment and the more delocalized picture based on local-density calculations.

Ab initio calculation of electronic coupling in the photosynthetic reaction center

We have carried out an ab initio electronic structure calculations of electron transfer couplings between chromophores in the bacterial photosynthetic reaction center. The couplings agree remarkably well with parameters obtained from recent quantum dynamical modeling of experimental data assuming an explicit intermediate mechanism. We also have computed couplings on the M-side of the reaction center and have found that the interaction of the primary donor to the M-side intermediate bacteriochlorophyll is quite small because of destructive interference of the two localized coupling matrix elements. This may explain the slow rate of electron transfer down the M-side of the reaction center.

Ab initio electronic structures of rhombohedral and cubic HgXO3 (X = Ti, Pb)

First-principles calculations were performed for orthorhombic HgO, rhombohedral and cubic phases of HgTiO3 (HTO) and HgPbO3 (HPO). The calculations show that in the rhombohedral phase HTO is a direct gap insulator with a gap of ~1.6 eV. The rhombohedral phase of HPO, on the other hand, shows a weak metallic character. The results provide an explanation for the electrical properties of these compounds. The cubic phases of HTO and HPO are invariably metallic in nature, thereby suggesting that for HTO the rhombohedral–cubic transition must also be accompanied by a change in the electrical state. Examination of the electronic density of states of these systems revealed no significant on-site mixing of Hg 5d and Hg 6s states in any of these materials.

LDA or GGA? A combined experimental inelastic neutron scattering and ab initio lattice dynamics study of alkali metal hydrides

In a previous work, we carried out inelastic neutron scattering (INS) spectroscopy experiments and preliminary first principles calculations on alkali metal hydrides. The complete series of alkali metal hydrides, LiH, NaH, KH, RbH and CsH was measured in the high-resolution TOSCA INS spectrometer at ISIS. Here, we present the results of ab initio electronic structure calculations of the properties of the alkali metal hydrides using both the local density approximation (LDA) and the generalized gradient approximation (GGA), using the Perdew–Burke–Ernzerhof (PBE) parameterization. Properties calculated were lattice parameters, bulk moduli, dielectric constants, effective charges, electronic densities and inelastic neutron scattering (INS) spectra. We took advantage of the currently available computer power to use full lattice dynamics theory to calculate thermodynamic properties for these materials. For the alkali metal hydrides (LiH, NaH, KH, RbH and CsH) using lattice dynamics, we found that the INS spectra calculated using LDA agreed better with the experimental data than the spectra calculated using GGA. Both zero-point effects and thermal contributions to free energies had an important effect on INS and several thermodynamic properties.

Facilitating Realistic and Efficient Applications of Ab Initio Electronic Dynamics for Polyatomic Systems

Thesis (Ph.D.)--University of Washington, 2016-06

Structure and properties of quaternary and tetragonal Heusler compounds for spintronics and spin transfer torque applications

Diese Dissertation ist in zwei Teile aufgeteilt: Teil 1 befasst sich mit der Vorhersage von Halb-Metallizität in quarternären Heuslerverbindungen und deren Potential für Spintronik-Anwendungen. Teil 2 befasst sich mit den strukturellen Eigenschaften der Mn2-basierenden Heuslerverbindungen und dem Tuning von ihrer magnetischen Eigenschaften bzgl. Koerzitivfeldstärke und Remanenz. Diese Verbindungen sind geeignet für Spin-Transfer Torque-Anwendungen.rnrnIn Teil 1 wurden die folgenden drei Probenserien quarternärer Heuslerverbindungen untersucht: XX´MnGa (X = Cu, Ni und X´ = Fe, Co), CoFeMnZ (Z = Al, Ga, Si, Ge) und Co2−xRhxMnZ (Z = Ga, Sn, Sb). Abgesehen von CuCoMnGa wurden alle diese Verbindungen mittels ab-initio Bandstrukturrechnungen als halbmetallische Ferromagnete prognostiziert. In der XX´MnGa-Verbindungsklasse besitzt NiFeMnGa zwar eine zu niedrige Curie-Temperatur für technologische Anwendungen, jedoch NiCoMnGa mit seiner hohen Spinpolarisation, einem hohen magnetischen Moment und einer hohen Curie-Temperatur stellt ein neues Material für Spintronik-Anwendungen dar. Alle CoFeMnZ-Verbindungen kristallisieren in der kubischen Heuslerstruktur und ihre magnetischen Momente folgen der Slater-Pauling-Regel, was Halbmetalizität und eine hohe Spinpolarisation impliziert. Die ebenfalls hohen Curie-Temperaturen ermöglichen einen Einsatz weit über Raumtemperatur hinaus. In der strukturellen Charakterisierung wurde festgestellt, dass sämtliche Co2−xRhxMnZ abgesehen von CoRhMnSn verschiedene Typen von Unordnung aufweisen; daher war die ermittelte Abweichung von der Slater-Pauling-Regel sowie von der 100%-igen Spinpolarisation dieser Verbindungen zu erwarten. Die Halbmetallizität der geordneten CoRhMnSn-Verbindung sollte nach den durchgeführten magnetischen Messungen vorhanden sein.rnrnIm zweiten Teil wurden Mn3−xCoxGa und Mn2−xRh1+xSn synthetisiert und charakterisiert. Es wurde gezeigt, dass Mn3−xCoxGa im Bereich x = 0.1 − 0.4 in einer tetragonal verzerrten inversen Heuslerstruktur kristallisiert und im Bereich x = 0.6−1 in einer kubisch inversen Heuslerstruktur. Während die tetragonalen Materialien hartmagnetisch sind und Charakeristika aufweisen, die typischerweise für Spin-Transfer Torque-Anwengungen attraktiv sind, repräsentieren die weichmagnetischen kubischen Vertreter die 100% spinpolarisierten Materialien, die der Slater-Pauling-Regel folgen. Mn2RhSn kristallisiert in der inversen tetragonal verzerrten Heuslerstruktur, weist einernhartmagnetische Hystereseschleife auf und folgt nicht der Slater-Pauling-Regel. Bei hohen Rh-Gehalt wird die kubische inverse Heuslerstruktur gebildet. Alle kubischen Proben sind weichmagnetisch und folgen der Slater-Pauling-Regel.

Experimental and theoretical investigation of the structural, chemical, electronic, and high frequency dielectric properties of barium cadmium tantalate-based ceramics

Single-phase Ba(Cd1/3Ta2/3)O-3 powder was produced using conventional solid state reaction methods. Ba(Cd1/3Ta2/3)O-3 ceramics with 2 wt % ZnO as sintering additive sintered at 1550 degreesC exhibited a dielectric constant of similar to32 and loss tangent of 5x10(-5) at 2 GHz. X-ray diffraction and thermogravimetric measurements were used to characterize the structural and thermodynamic properties of the material. Ab initio electronic structure calculations were used to give insight into the unusual properties of Ba(Cd1/3Ta2/3)O-3, as well as a similar and more widely used material Ba(Zn1/3Ta2/3)O-3. While both compounds have a hexagonal Bravais lattice, the P321 space group of Ba(Cd1/3Ta2/3)O-3 is reduced from P (3) under bar m1 of Ba(Zn1/3Ta2/3)O-3 as a result of a distortion of oxygen away from the symmetric position between the Ta and Cd ions. Both of the compounds have a conduction band minimum and valence band maximum composed of mostly weakly itinerant Ta 5d and Zn 3d/Cd 4d levels, respectively. The covalent nature of the directional d-electron bonding in these high-Z oxides plays an important role in producing a more rigid lattice with higher melting points and enhanced phonon energies, and is suggested to play an important role in producing materials with a high dielectric constant and low microwave loss. (C) 2005 American Institute of Physics.

Negative-ion photoelectron spectroscopy, gas-phase acidity, and thermochemistry of the peroxyl radicals CH3OO and CH3CH2OO

Methyl, methyl-d(3), and ethyl hydroperoxide anions (CH3OO-, CD3OO-, and CH3CH2OO-) have been prepared by deprotonation of their respective hydroperoxides in a stream of helium buffer, gas. Photodetachment with 364 nm (3.408 eV) radiation was used to measure the adiabatic electron affinities: EA[CH3OO, (X) over tilde (2)A"] = 1.161 +/- 0.005 eV, EA[CD3OO, (X) over tilde (2)A"] = 1.154 +/- 0.004 eV, and EA[CH3CH2OO, (X) over tilde (2)A"] = 1.186 +/- 0.004 eV. The photoelectron spectra yield values for the term energies: DeltaE((X) over tilde 2A"-(A) over tilde 2A')[CH3OO] = 0.914 +/- 0.005 eV, DeltaE((X) over tilde (2)A"-(A) over tilde 2A') [CD3OO] = 0.913 +/- 0.004 eV, and DeltaE((X) over tilde (2)A"-(A) over tilde (2)A')[CH3CH2OO] = 0.938 +/- 0.004 eV. A localized RO-O stretching mode was observed near 1100 cm(-1) for the ground state of all three radicals, and low-frequency R-O-O bending modes are also reported. Proton-transfer kinetics of the hydroperoxides have been measured in a tandem flowing afterglow-selected ion flow tube k(FA-SIFT) to determine the gas-phase acidity of the parent hydroperoxides: Delta (acid)G(298)(CH3OOH) = 367.6 +/- 0.7 kcal mol(-1), Delta (acid)G(298)(CD3OOH) = 367.9 +/- 0.9 kcal mol(-1), and Delta (acid)G(298)(CH3CH2OOH) = 363.9 +/- 2.0 kcal mol(-1). From these acidities we have derived the enthalpies of deprotonation: Delta H-acid(298)(CH3OOH) = 374.6 +/- 1.0 kcal mol(-1), Delta H-acid(298)(CD3OOH) = 374.9 +/- 1.1 kcal mol(-1), and Delta H-acid(298)(CH2CH3OOH) = 371.0 +/- 2.2 kcal mol(-1). Use of the negative-ion acidity/EA cycle provides the ROO-H bond enthalpies: DH298(CH3OO-H) 87.8 +/- 1.0 kcal mol(-1), DH298(CD3OO-H) = 87.9 +/- 1.1 kcal mol(-1), and DH298(CH3CH2OO-H) = 84.8 +/- 2.2 kcal mol(-1). We review the thermochemistry of the peroxyl radicals, CH3OO and CH3CH2OO. Using experimental bond enthalpies, DH298(ROO-H), and CBS/APNO ab initio electronic structure calculations for the energies of the corresponding hydroperoxides, we derive the heats of formation of the peroxyl radicals. The "electron affinity/acidity/CBS" cycle yields Delta H-f(298)[CH3OO] = 4.8 +/- 1.2 kcal mol(-1) and Delta H-f(298)[CH3CH2OO] = -6.8 +/- 2.3 kcal mol(-1).

Theoretical analysis of photoinduced H-atom elimination in thiophenol

The photoinduced hydrogen elimination reaction in thiophenol via the conical intersections of the dissociative (1)pi sigma* excited state with the bound (1)pi pi* excited state and the electronic ground state has been investigated with ab initio electronic-structure calculations and time-dependent quantum wave-packet calculations. A screening of the coupling constants of the symmetry-allowed coupling modes at the (1)pi pi*-(1)pi sigma* and (1)pi sigma*-S-0 conical intersection shows that the SH torsional mode is by far the most important coupling mode at both conical intersections. A model including three intersecting potential-energy surfaces (S-0, (1)pi pi*, (1)pi sigma*) and two nuclear degrees of freedom (SH stretch and SH torsion) has been constructed on the basis of ab initio complete-active-space self-consistent field and multireference second-order perturbation theory calculations. The nonadiabatic quantum wave-packet dynamics initiated by optical excitation of the (1)pi pi* and (1)pi sigma* states has been explored for this three-state two-coordinate model. The photodissociation dynamics is characterized in terms of snapshots of time-dependent wave packets, time-dependent electronic population probabilities, and the branching ratio of the (2)sigma/(2)pi electronic states of the thiophenoxyl radical. The dependence of the timescale of the photodissociation process and the branching ratio on the initial excitation of the SH stretching and SH torsional vibrations has been analyzed. It is shown that the node structure, which is imposed on the nuclear wave packets by the initial vibrational preparation as well as by the transitions through the conical intersections, has a profound effect on the photodissociation dynamics. The effect of additional weak coupling modes of CC twist (nu(16a)) and ring-distortion (nu(16b)) character has been investigated with three-dimensional and four-dimensional time-dependent wave-packet calculations, and has been found to be minor. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4709608]

Theoretical studies of novel organic systems : biradicals, polyradicals, and conducting polymers

Chapter 1

Cyclobutanediyl has been studied in both its singlet and triplet states by ab initio electronic structure theory. The triplet, which is the ground state of the molecule, exists in both C_(2h) and C_(2v) forms, which interconvert via a C_s transition state. For the singlet, only a C_(2h) form is found. It passes, via a C_s transition state, onto the C_(2v) surface on which bicyclobutane is the only minimum. The ring-flipping (inversion) process in bicyclobutane includes the singlet biradical as an intermediate, and involves a novel, nonleast motion pathway. Semiclassical periodic orbit theory indicates that the various minima on both the singlet and triplet surfaces can interconvert via quantum mechanical tunneling.

Chapter 2

The dimethylenepolycyclobutadienes (n) are the non-Kekulé analogues of the classical acenes. Application of a variety of theoretical methods reveals several novel features of such structures. Most interesting is the emergence of a parity rule. When n is even, n is predicted to be a singlet, with n disjoint NBMOs. When n is odd, theory predicts a triplet ground state with (n+1) NBMOs that are not fully disjoint.

Chapter 3

Bi(cyclobutadienyl) (2), the cyclobutadiene analogue of biphenyl, and its homologues tri- (3) and tetra(cyclobutadienyl) (4) have been studied using electronic structure theory. Ab initio calculations on 2 reveal that the central bond is a true double bond, and that the structure is best thought of as two allyl radicals plus an ethylene. The singlet and triplet states are essentially degenerate. Trimer 3 is two allyls plus a dimethylenecyclobutanediyl, while 4 is two coplanar bi(cyclobutadienyl) units connected by a single bond. For both 3 and 4, the quintet, triplet, and singlet states are essentially degenerate, indicating that they are tetraradicals. The infinite polymer, polycyclobutadiene, has been studied by HMO, EHCO, and VEH methods. Several geometries based on the structures of 3 and 4 have been studied, and the band structures are quite intriguing. A novel crossing between the valence and conduction bands produces a small band gap and a high density of states at the Fermi level.

Chapter 4

At the level of Hückel theory, polyfulvene has a HOCO-LUCO degeneracy much like that seen in polyacetylene. Higher levels of theory remove the degeneracy, but the band gap (E_g) is predicted to be significantly smaller than analogous structures such as polythiophene and polypyrrole at the fulvenoid geometry. An alternative geometry, which we have termed quinoid, is also conceivable for polyfulvene, and it is predicted to have a much larger E_g. The effects of benzannelation to produce analogues of polyisothianaphthene have been evaluated. We propose a new model for such structures based on conventional orbital mixing arguments. Several of the proposed structures have quite interesting properties, which suggest that they are excellent candidates for conducting polymers.

Chapter 5

Theoretical studies of polydimethylenecyclobutene and polydiisopropylidene- cyclobutene reveal that, because of steric crowding, they cannot achieve a planar, fully conjugated structure in either their undoped or doped states. Rather, the structure consists of essentially orthogonal hexatriene units. Such a structure is incompatible with conventional conduction mechanisms involving polarons and bipolarons.

First-principles study of transition metal impurities in Si

By using ab initio electronic structure calculations within density functional theory, we study the structural, electronic, and magnetic properties of Si doped with a transition metal impurity. We consider the transition metals of the 3d series V, Cr, Mn, Fe, Co, and Ni. To get insight into the level filling mechanism and the magnetization saturation, we first investigate the transition metal-Si alloys in the zinc-blende structure. Next, we investigate the doping of bulk Si with a transition metal atom, in which it occupies the substitutional site, the interstitial site with tetrahedral symmetry, and the interstitial site with hexagonal symmetry. It is found that all of these transition metal impurities prefer an interstitial position in Si. Furthermore, we show that it is possible to interpret the electronic and magnetic properties by using a simple level filling picture and a comparison is made to Ge doped with the same transition metal atoms. In order to get insight into the effect of a strained environment, we calculate the formation energy as a function of an applied homogeneous pressure and we show that an applied pressure can stabilize the substitutional position of transition metal impurities in Si. Finally, the energies of the ferromagnetic states are compared to those of the antiferromagnetic states. It is shown that the interstitial site of the Mn dopant helps us to stabilize the nearest neighbor substitutional site to realize the ferromagnetic state. For doping of Si with Cr, a ferrimagnetic behavior is predicted.