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.

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.

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.

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.

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.

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

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

The catalytic role of an isolated-Ti atom in the hydrogenation of Ti-doped Al(001) surface : An ab initio density functional theory calculation

Recent work [S. Chaudhuri, J.T. Muckerman, J. Phys. Chem. B 109 (2005) 6952] reported that two Ti-substituted atoms on an Al(0 0 1) surface can form a catalytically active site for the dissociation of H2, but the diffusion barrier of atomic H away from Ti site is as high as 1.57 eV. By using ab initio density functional calculations, we found that two hydrogen molecules can dissociate on isolated-Ti atom doped Al(0 0 1) surface with small activation barriers (0.21 and 0.235 eV for first and second H2, respectively). Additionally, the diffusion barrier of atomic H away from Ti site is also moderate (0.47 eV). These results contribute further towards understanding the improved kinetics observed in recycling of hydrogen with Ti-doped NaAlH4.

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.

CO oxidation on Pt(111):An ab initio density functional theory study

CO oxidation on Pt(111) is studied with ab initio density functional theory. The low energy pathway and transition state for the reaction are identified. The key event is the breaking of an O-metal bond prior to the formation of a chemisorbed CO2 molecule. The pathway can be rationalized in terms of competition of the O and C atoms for bonding with the underlying surface, and the predominant energetic barrier is the strength of the O-metal bond.

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.

The Role of Anharmonicity in Hydrogen-Bonded Systems: The Case of Water Clusters

The nature of vibrational anharmonicity has been examined for the case of small water clusters using second-order vibrational perturbation theory (VPT2) applied on second-order Møller–Plesset perturbation theory (MP2) potential energy surfaces. Using a training set of 16 water clusters (H2O)n=2–6,8,9 with a total of 723 vibrational modes, we determined scaling factors that map the harmonic frequencies onto anharmonic ones. The intermolecular modes were found to be substantially more anharmonic than intramolecular bending and stretching modes. Due to the varying levels of anharmonicity of the intermolecular and intramolecular modes, different frequency scaling factors for each region were necessary to achieve the highest accuracy. Furthermore, new scaling factors for zero-point vibrational energies (ZPVE) and vibrational corrections to the enthalpy (ΔHvib) and the entropy (Svib) have been determined. All the scaling factors reported in this study are different from previous works in that they are intended for hydrogen-bonded systems, while others were built using experimental frequencies of covalently bonded systems. An application of our scaling factors to the vibrational frequencies of water dimer and thermodynamic functions of 11 larger water clusters highlights the importance of anharmonic effects in hydrogen-bonded systems.

Benchmark Structures and Binding Energies of Small Water Clusters with Anharmonicity Corrections

For (H2O)n where n = 1–10, we used a scheme combining molecular dynamics sampling with high level ab initio calculations to locate the global and many low lying local minima for each cluster. For each isomer, we extrapolated the RI-MP2 energies to their complete basis set limit, included a CCSD(T) correction using a smaller basis set and added finite temperature corrections within the rigid-rotor-harmonic-oscillator (RRHO) model using scaled and unscaled harmonic vibrational frequencies. The vibrational scaling factors were determined specifically for water clusters by comparing harmonic frequencies with VPT2 fundamental frequencies. We find the CCSD(T) correction to the RI-MP2 binding energy to be small (<1%) but still important in determining accurate conformational energies. Anharmonic corrections are found to be non-negligble; they do not alter the energetic ordering of isomers, but they do lower the free energies of formation of the water clusters by as much as 4 kcal/mol at 298.15 K.

Accurate Predictions of Water Cluster Formation, (H2O)n=2-10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.