Advanced density functional theory methods for materials science

In this work we chiefly deal with two broad classes of problems in computational materials science, determining the doping mechanism in a semiconductor and developing an extreme condition equation of state. While solving certain aspects of these questions is well-trodden ground, both require extending the reach of existing methods to fully answer them. Here we choose to build upon the framework of density functional theory (DFT) which provides an efficient means to investigate a system from a quantum mechanics description.

Zinc Phosphide (Zn₃P₂) could be the basis for cheap and highly efficient solar cells. Its use in this regard is limited by the difficulty in n-type doping the material. In an effort to understand the mechanism behind this, the energetics and electronic structure of intrinsic point defects in zinc phosphide are studied using generalized Kohn-Sham theory and utilizing the Heyd, Scuseria, and Ernzerhof (HSE) hybrid functional for exchange and correlation. Novel 'perturbation extrapolation' is utilized to extend the use of the computationally expensive HSE functional to this large-scale defect system. According to calculations, the formation energy of charged phosphorus interstitial defects are very low in n-type Zn₃P₂ and act as 'electron sinks', nullifying the desired doping and lowering the fermi-level back towards the p-type regime. Going forward, this insight provides clues to fabricating useful zinc phosphide based devices. In addition, the methodology developed for this work can be applied to further doping studies in other systems.

Accurate determination of high pressure and temperature equations of state is fundamental in a variety of fields. However, it is often very difficult to cover a wide range of temperatures and pressures in an laboratory setting. Here we develop methods to determine a multi-phase equation of state for Ta through computation. The typical means of investigating thermodynamic properties is via ’classical’ molecular dynamics where the atomic motion is calculated from Newtonian mechanics with the electronic effects abstracted away into an interatomic potential function. For our purposes, a ’first principles’ approach such as DFT is useful as a classical potential is typically valid for only a portion of the phase diagram (i.e. whatever part it has been fit to). Furthermore, for extremes of temperature and pressure quantum effects become critical to accurately capture an equation of state and are very hard to capture in even complex model potentials. This requires extending the inherently zero temperature DFT to predict the finite temperature response of the system. Statistical modelling and thermodynamic integration is used to extend our results over all phases, as well as phase-coexistence regions which are at the limits of typical DFT validity. We deliver the most comprehensive and accurate equation of state that has been done for Ta. This work also lends insights that can be applied to further equation of state work in many other materials.

Fundamental Studies of Carbon Oxygen Bond Activation in Nickel Diphosphine Ether Complexes. And, Metallomacrocycles as Ligands: Synthesis and Characterization of Aluminum-Bridged Bisglyoximato Complexes of Iron and Cobalt

In order to develop better catalysts for the cleavage of aryl-X bonds fundamental studies of the mechanism and individual steps of the mechanism have been investigated in detail. As the described studies are difficult at best in catalytic systems, model systems are frequently used. To study aryl-oxygen bond activation, a terphenyl diphosphine scaffold containing an ether moiety in the central arene was designed. The first three chapters of this dissertation focus on the studies of the nickel complexes supported by this diphosphine backbone and the research efforts in regards to aryl-oxygen bond activation.

Chapter 2 outlines the synthesis of a variety of diphosphine terphenyl ether ligand scaffolds. The metallation of these scaffolds with nickel is described. The reactivity of these nickel(0) systems is also outlined. The systems were found to typically undergo a reductive cleavage of the aryl oxygen bond. The mechanism was found to be a subsequent oxidative addition, β-H elimination, reductive elimination and (or) decarbonylation.

Chapter 3 presents kinetic studies of the aryl oxygen bond in the systems outlined in Chapter 2. Using a series of nickel(0) diphosphine terphenyl ether complexes the kinetics of aryl oxygen bond activation was studied. The activation parameters of oxidative addition for the model systems were determined. Little variation was observed in the rate and activation parameters of oxidative addition with varying electronics in the model system. The cause of the lack of variation is due to the ground state and oxidative addition transition state being affected similarly. Attempts were made to extend this study to catalytic systems.

Chapter 4 investigates aryl oxygen bond activation in the presence of additives. It was found that the addition of certain metal alkyls to the nickel(0) model system lead to an increase in the rate of aryl oxygen bond activation. The addition of excess Grignard reagent led to an order of magnitude increase in the rate of aryl oxygen bond activation. Similarly the addition of AlMe3 led to a three order of magnitude rate increase. Addition of AlMe₃ at -80 °C led to the formation of an intermediate which was identified by NOESY correlations as a system in which the AlMe3 is coordinated to the ether moiety of the backbone. The rates and activation parameters of aryl oxygen bond activation in the presence of AlMe3 were investigated.

The last two chapters involve the study of metalla-macrocycles as ligands. Chapter 5 details the synthesis of a variety of glyoxime backbones and diphenol precursors and their metallation with aluminum. The coordination chemistry of iron on the aluminum scaffolds was investigated. Varying the electronics of the aluminum macrocycle was found to affect the observed electrochemistry of the iron center.

Chapter 6 extends the studies of chapter 5 to cobalt complexes. The synthesis of cobalt dialuminum glyoxime metal complexes is described. The electrochemistry of the cobalt complexes was investigated. The electrochemistry was compared to the observed electrochemistry of a zinc analog to identify the redox activity of the ligand. In the presence of acid the cobalt complexes were found to electrochemically reduce protons to dihydrogen. The electronics of the ancillary aluminum ligands were found to affect the potential of proton reduction in the cobalt complexes. These potentials were compared to other diglyoximate complexes.

Electrocatalytic activity of transition metal substituted heteropolytungstates

The electrochemical and electrocatalytic behavior of a series of heteropolytungstate anions in which a tungsten atom in the well known Keggin structure has been replaced by an iron atom is described. All of the iron substituted ions exhibit a one electron reversible couple associated with the Fe³⁺ center and a pair of two electron waves attributed to electron addition and removal from the tungsten oxo framework. The pH and ionic strength effects upon the various electrochemical processes are discussed and interpreted in terms of a competition between protonation and ion pairing of the highly negatively charged ions.

The anions are efficient catalysts for the electroreduction of hydrogen peroxide. A catalytic mechanism involving a formally Fe(IV) intermediate is proposed. Pulse radiolysis experiments were employed to detect the intermediate and evaluate the rate constants for the reactions in which it is formed and decomposed. A chain mechanism for the catalytic decomposition of hydrogen peroxide in which the Fe center shuttles between the +2, +3, and +4 oxidation states is proposed to explain the non-integral stoichiometry observed for the iron substituted polytungstate catalyzed electroreduction of hydrogen peroxide.

The anions are also efficient electrocatalyst for the electrochemical conversion of nitric oxide to ammonia. The catalyzed reduction does not produce hydroxylamine as an intermediate and appears to depend upon the ability of the multiply reduced heteropolytungstates to deliver several electrons to the bound NO group in a concerted step. A valuable feature of the heteropolytungstates is the ease at which the formal potentials of the several redox couples they exhibit may be shifted by changing the identity of the central heteroatom. Exploitation of this feature provided diagnostic information that was decisive in establishing the mechanism of electrocatalytic reduction.

The iron substituted heteropolytungstates are not degraded by repeated cycling between their oxidized and reduced states. They also show superior activity compared to their unsubstituted analogues, indicating that the Fe center acts as a binding site that facilitates inner-sphere electron transfer processes. The basic electrochemistry of several other transition metal substituted Keggin ions is also described.

A pessoa estomizada e o processo de inclusão no trabalho: contribuição para a enfermagem

Estudo cujo objeto tratou da inclusão do cliente estomizado no mundo do trabalho. Os objetivos foram: identificar as dificuldades e facilidades dos clientes estomizados para inclusão no mundo do trabalho; analisar as possibilidades de inclusão no mundo do trabalho; discutir, a partir do ponto de vista do cliente estomizado, as orientações fornecidas pelos enfermeiros com vistas à inclusão no mundo do trabalho. O referencial teórico baseou-se no campo da Saúde do Trabalhador, enfocando a reabilitação profissional e o capítulo de bases conceituais abordou o conhecimento da estomaterapia, do mundo do trabalho, da deficiência física e aspectos legais que envolvem a reabilitação do estomizado no mundo laboral. O desenho metodológico foi de uma pesquisa descritiva, exploratória, de natureza qualitativa, realizada com 20 clientes estomizados definitivos, aos quais se aplicou uma entrevista semiestruturada. O método de análise dos dados foi a Análise Temática de Conteúdo, a qual fez emergir quatro categorias: a) Sentidos do Trabalho para o Ser Estomizado; b) O Estomizado e Sua Problemática Biopsicossocial; c) Contexto Social e Aspectos Legais Envolvendo a Inclusão do Estomizado no Mundo Laboral; d) O Enfermeiro e Sua Participação na Reabilitação do Cliente Estomizado. Os resultados revelaram que a maioria dos sujeitos trabalhava informalmente e recebia ao mesmo tempo algum auxílio governamental. Ressalta-se que eles reconheciam a ilegalidade desta situação, porém, julgavam-na necessária devido aos baixos valores dos benefícios, enfatizando-se a sensação de utilidade causada pelo fato de trabalharem. Referiram que o retorno ao trabalho era prejudicado devido a empecilhos encontrados nas dimensões psíquica, física e social, as quais estavam articuladas intimamente. Enfatizaram grande dificuldade em encontrarem empregos adequados às suas especificidades, e que não prejudicassem sua condição de saúde, pois há necessidade de banheiros adaptados, de não exposição ao calor na região do estoma e nem a esforços físicos severos, sem contar com a necessidade de um emprego que lhes permita flexibilidade para irem às consultas da equipe multiprofissional. Os maiores empecilhos sociais referiram-se ao desconhecimento e descaso social e governamental a respeito do que é ser estomizado, pois esta problemática não é divulgada, e nem conhecida pela maioria da população. Em relação aos enfermeiros, os sujeitos foram quase unânimes em referirem falta de orientação por parte desses profissionais, acerca de esclarecerem sobre sua inclusão no mundo do trabalho. Este fato caracterizou-se como preocupante, pois os enfermeiros são educadores por excelência e a orientação está intimamente ligada ao processo de reabilitação. Concluiu-se que o retorno ao trabalho foi considerado essencial, mas existem inúmeras dificuldades para que este retorno e manutenção no universo laboral. Estes empecilhos os levam a adquirirem aposentadorias precoces ou auxílios-doença. Há de se rever o processo de reabilitação da pessoa com estoma, especialmente no que se refere a sua inclusão no mundo do trabalho, no sentido de melhor prepará-la para suas potencialidades e limitações, destacando-se que ela não é incapaz e que existem atividades formais em que elas podem ser produtivas e felizes.

The relative abundances of Sn, Te, Xe, Ba, and Ce in the cosmic radiation

Elements with even atomic number (Z) in the interval 50 ≤ Z ≤ 58 have been resolved in the cosmic radiation using the Heavy Nuclei Experiment on the HEAO-3 satellite. Their relative abundances have been compared with the results expected from pure r-process material, pure s-process material, and solar system material, both with and without a modification due to possible first ionization potential effects. Such effects may be the result of the preferential acceleration, and hence enhancement in the cosmic rays, or those elements having low first ionization potentials. We find that our measurements are inconsistent with pure r-process material at the greater than 98% confidence level whether or not the first ionization potential adjustments are made.

In addition, we have compared our results with mixtures having varying ratios of pure r-process material to pure s-process material. We find that, if no first ionization potential effects are included,

(r/s)CRS/(r/s)SS = 0.20^+0.18_-0.14

where CRS refers to the cosmic ray source and SS refers to the solar system, consistent with having an almost pure s-process source. If the first ionization potential adjustments are applied

(r/s)CRS/(r/s)SS = 1.5^+1.1_-0.7

consistent with a solar system mixture.

Methods for melting temperature calculation

Melting temperature calculation has important applications in the theoretical study of phase diagrams and computational materials screenings. In this thesis, we present two new methods, i.e., the improved Widom's particle insertion method and the small-cell coexistence method, which we developed in order to capture melting temperatures both accurately and quickly.

We propose a scheme that drastically improves the efficiency of Widom's particle insertion method by efficiently sampling cavities while calculating the integrals providing the chemical potentials of a physical system. This idea enables us to calculate chemical potentials of liquids directly from first-principles without the help of any reference system, which is necessary in the commonly used thermodynamic integration method. As an example, we apply our scheme, combined with the density functional formalism, to the calculation of the chemical potential of liquid copper. The calculated chemical potential is further used to locate the melting temperature. The calculated results closely agree with experiments.

We propose the small-cell coexistence method based on the statistical analysis of small-size coexistence MD simulations. It eliminates the risk of a metastable superheated solid in the fast-heating method, while also significantly reducing the computer cost relative to the traditional large-scale coexistence method. Using empirical potentials, we validate the method and systematically study the finite-size effect on the calculated melting points. The method converges to the exact result in the limit of a large system size. An accuracy within 100 K in melting temperature is usually achieved when the simulation contains more than 100 atoms. DFT examples of Tantalum, high-pressure Sodium, and ionic material NaCl are shown to demonstrate the accuracy and flexibility of the method in its practical applications. The method serves as a promising approach for large-scale automated material screening in which the melting temperature is a design criterion.

We present in detail two examples of refractory materials. First, we demonstrate how key material properties that provide guidance in the design of refractory materials can be accurately determined via ab initio thermodynamic calculations in conjunction with experimental techniques based on synchrotron X-ray diffraction and thermal analysis under laser-heated aerodynamic levitation. The properties considered include melting point, heat of fusion, heat capacity, thermal expansion coefficients, thermal stability, and sublattice disordering, as illustrated in a motivating example of lanthanum zirconate (La₂Zr₂O₇). The close agreement with experiment in the known but structurally complex compound La₂Zr₂O₇ provides good indication that the computation methods described can be used within a computational screening framework to identify novel refractory materials. Second, we report an extensive investigation into the melting temperatures of the Hf-C and Hf-Ta-C systems using ab initio calculations. With melting points above 4000 K, hafnium carbide (HfC) and tantalum carbide (TaC) are among the most refractory binary compounds known to date. Their mixture, with a general formula Ta_xHf_1-xC_y, is known to have a melting point of 4215 K at the composition Ta₄HfC₅, which has long been considered as the highest melting temperature for any solid. Very few measurements of melting point in tantalum and hafnium carbides have been documented, because of the obvious experimental difficulties at extreme temperatures. The investigation lets us identify three major chemical factors that contribute to the high melting temperatures. Based on these three factors, we propose and explore a new class of materials, which, according to our ab initio calculations, may possess even higher melting temperatures than Ta-Hf-C. This example also demonstrates the feasibility of materials screening and discovery via ab initio calculations for the optimization of "higher-level" properties whose determination requires extensive sampling of atomic configuration space.

TWO-DIMENSIONAL ARRAY OF ION TRAPS ON A PLANAR CHIP

We investigate a planar ion chip design with a two-dimensional array of linear ion traps for the scalable quantum information processor. The segmented electrodes reside in a single plane on a substrate and a grounded metal plate, a combination of appropriate rf and DC potentials are applied to them for stable ion confinement, and the trap axes are located above the surface at a distance controlled by the electrodes' lateral extent and the substrate's height as discussed. The potential distributions are calculated using static electric field qualitatively. This architecture is conceptually simple and many current microfabrication techniques are feasible for the basic structure. It may provide a promising route for scalable quantum computers.

High pressure states in condensed matter: I. High pressure behavior of the iron-sulphur system with applications to the earth's core. II. Empirical equation of state for organic compounds at high pressures

Part I:

The earth's core is generally accepted to be composed primarily of iron, with an admixture of other elements. Because the outer core is observed not to transmit shear waves at seismic frequencies, it is known to be liquid or primarily liquid. A new equation of state is presented for liquid iron, in the form of parameters for the 4th order Birch-Murnaghan and Mie-Grüneisen equations of state. The parameters were constrained by a set of values for numerous properties compiled from the literature. A detailed theoretical model is used to constrain the P-T behavior of the heat capacity, based on recent advances in the understanding of the interatomic potentials for transition metals. At the reference pressure of 10⁵ Pa and temperature of 1811 K (the normal melting point of Fe), the parameters are: ρ = 7037 kg/m³, K_S0 = 110 GPa, K_S' = 4.53, K_S" = -.0337 GPa-1, and γ = 2.8, with γ α ρ^-1.17. Comparison of the properties predicted by this model with the earth model PREM indicates that the outer core is 8 to 10 % less dense than pure liquid Fe at the same conditions. The inner core is also found to be 3 to 5% less dense than pure liquid Fe, supporting the idea of a partially molten inner core. The density deficit of the outer core implies that the elements dissolved in the liquid Fe are predominantly of lower atomic weight than Fe. Of the candidate light elements favored by researchers, only sulfur readily dissolves into Fe at low pressure, which means that this element was almost certainly concentrated in the core at early times. New melting data are presented for FeS and FeS₂ which indicate that the FeS₂ is the S-hearing liquidus solid phase at inner core pressures. Consideration of the requirement that the inner core boundary be observable by seismological means and the freezing behavior of solutions leads to the possibility that the outer core may contain a significant fraction of solid material. It is found that convection in the outer core is not hindered if the solid particles are entrained in the fluid flow. This model for a core of Fe and S admits temperatures in the range 3450K to 4200K at the top of the core. An all liquid Fe-S outer core would require a temperature of about 4900 K at the top of the core.

Part II.

The abundance of uses for organic compounds in the modern world results in many applications in which these materials are subjected to high pressures. This leads to the desire to be able to describe the behavior of these materials under such conditions. Unfortunately, the number of compounds is much greater than the number of experimental data available for many of the important properties. In the past, one approach that has worked well is the calculation of appropriate properties by summing the contributions from the organic functional groups making up molecules of the compounds in question. A new set of group contributions for the molar volume, volume thermal expansivity, heat capacity, and the Rao function is presented for functional groups containing C, H, and O. This set is, in most cases, limited in application to low molecular liquids. A new technique for the calculation of the pressure derivative of the bulk modulus is also presented. Comparison with data indicates that the presented technique works very well for most low molecular hydrocarbon liquids and somewhat less well for oxygen-bearing compounds. A similar comparison of previous results for polymers indicates that the existing tabulations of group contributions for this class of materials is in need of revision. There is also evidence that the Rao function contributions for polymers and low molecular compounds are somewhat different.

Electron transfer in chemically and genetically modified myoglobins

The temperature dependences of the reduction potentials (E^o') of wildtype human myoglobin (Mb) and three site-directed mutants have been measured by using thin-layer spectroelectrochemistry. Residue Val68, which is in van der Waals contact with the heme in Mb, has been replaced by Glu, Asp, and Asn. At pH 7.0, reduction of the heme iron (III) in the former two proteins is accompanied by uptake of a proton by the protein. The changes in E^o', and the standard entropy (ΔS^o') and enthalpy (ΔH^o') of reduction in the mutant proteins were determined relative to values for wild-type; the change in E^o' at 25°C was about -200 millivolts for the Glu and Asp mutants, and about -80 millivolts for the Asn mutant. Reduction of Fe(III) to Fe(II) in the Glu and Asp mutants is accompanied by uptake of a proton. These studies demonstrate that Mb can tolerate substitution of a buried hydrophobic group by potentially charged and polar residues, and that such amino acid replacements can lead to substantial changes in the redox thermodynamics of the protein.

Through analysis of the temperature dependence and shapes of NMR dispersion signals, it is determined that a water molecule is bound to the sixth coordination site of the ferric heme in the Val68Asp and in the Val68Asn recombinant proteins while the carboxyl group of the sidechain of Glu68 occupies this position in Val68Glu. The relative rhombic distortions in the ESR spectra of these mutant proteins combined with H₂¹⁷O and spin interconversion experiments performed on them confirm the conclusions of the NMRD study.

The rates of intramolecular electron transfer (ET) of (NH₃)₅Ru-His48 (Val68Asp, His81GIn, Cys110AIa)Mb and (NH₃)₅Ru-His48 (Val68GIu,His81GIn,Cys110Ala)Mb were measured to be .85(3)s^-1 and .30(2)s^-1, respectively. This data supports the hypothesis that entropy of 111 reduction and reorganization energy of ET are inversely related. The rates of forward and reverse ET for (NH₃)₅ Ru-His48 (Val68GIu, His81 GIn, Cys110AIa)ZnMb -7.2(5)•10⁴s^-1and 1.4(2)•10⁵s^-1, respectively- demonstrate that the placement of a highly polar residue nearby does not significantly change the reorganization energy of the photoactive Zn porphyrin.

The distal histidine imidazoles of (NH₃)₄isnRu-His48 SWMb and (NH₃)₅Ru-His48 SWMb were cyanated with BrCN. The intramolecular ET rates of these BrCN-modified Mb derivatives are 5.5(6)s^-1 and 3.2(5)s^-1, respectively. These respective rates are 20 and 10 times faster than those of their noncyanated counterparts after the differences in ET rate from driving force are scaled according to the Marcus equation. This increase in ET rate of the cyanated Mb derivatives is attributed to lower reorganization energy since the cyanated Mb heme is pentacoordinate in both oxidation states; whereas, the native Mb heme loses a water molecule upon reduction so that it changes from six to five coordinate. The reorganization energy from Fe-OH₂ dissociation is estimated to be .2eV. This conclusion is used to reconcile data from previous experiments in our lab. ET in photoactive porphyrin-substituted myoglobins proceed faster than predicted by Marcus Theory when it is assumed that the only difference in ET parameters between photoactive porphyrins and native heme systems is driving force. However, the data can be consistently fit to Marcus Theory if one corrects for the smaller reorganization in the photoactive porphyrin systems since they do not undergo a coordination change upon ET.

Finally, the intramolecular ET rate of (NH₃)₄isnRu-His48 SWMb was measured to be 3.0(4)s^-1. This rate is within experimental error of that for (NH₃)₄pyrRu-His48 SWMb even though the former has 80mV more driving force. One likely possibility for this observation is that the tetraamminepyridineruthenium group undergoes less reorganization upon ET than the tetraammineisonicotinamideruthenium group. Moreover, analysis of the (NH₃)₄isnRu-His48 SWMb experimental system gives a likely explanation of why ET was not observed previously in (NH₃)₄isnRu-Cytochrome C.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

Development and Application of Embedding Methods for the Simulation of Large Chemical Systems

The high computational cost of correlated wavefunction theory (WFT) calculations has motivated the development of numerous methods to partition the description of large chemical systems into smaller subsystem calculations. For example, WFT-in-DFT embedding methods facilitate the partitioning of a system into two subsystems: a subsystem A that is treated using an accurate WFT method, and a subsystem B that is treated using a more efficient Kohn-Sham density functional theory (KS-DFT) method. Representation of the interactions between subsystems is non-trivial, and often requires the use of approximate kinetic energy functionals or computationally challenging optimized effective potential calculations; however, it has recently been shown that these challenges can be eliminated through the use of a projection operator. This dissertation describes the development and application of embedding methods that enable accurate and efficient calculation of the properties of large chemical systems.

Chapter 1 introduces a method for efficiently performing projection-based WFT-in-DFT embedding calculations on large systems. This is accomplished by using a truncated basis set representation of the subsystem A wavefunction. We show that naive truncation of the basis set associated with subsystem A can lead to large numerical artifacts, and present an approach for systematically controlling these artifacts.

Chapter 2 describes the application of the projection-based embedding method to investigate the oxidative stability of lithium-ion batteries. We study the oxidation potentials of mixtures of ethylene carbonate (EC) and dimethyl carbonate (DMC) by using the projection-based embedding method to calculate the vertical ionization energy (IE) of individual molecules at the CCSD(T) level of theory, while explicitly accounting for the solvent using DFT. Interestingly, we reveal that large contributions to the solvation properties of DMC originate from quadrupolar interactions, resulting in a much larger solvent reorganization energy than that predicted using simple dielectric continuum models. Demonstration that the solvation properties of EC and DMC are governed by fundamentally different intermolecular interactions provides insight into key aspects of lithium-ion batteries, with relevance to electrolyte decomposition processes, solid-electrolyte interphase formation, and the local solvation environment of lithium cations.

Synthesis, Oxidation and Photophysics of Perfluoroborated Tetrakis(pyrophosphito)diplatinate (II) and Density Functional Theory (DFT) Study of Electrochemical CO2 Reduction by Mn Catalysts

In the first part of this thesis (Chapters I and II), the synthesis, characterization, reactivity and photophysics of per(difluoroborated) tetrakis(pyrophosphito)diplatinate(II) (Pt(POPBF₂)) are discussed. Pt(POP-BF₂) was obtained by reaction of [Pt₂(POP)₄]^4- with neat boron trifluoride diethyl etherate (BF₃·Et₂O). While Pt(POP-BF₂) and [Pt₂(POP)₄]^4- have similar structures and absorption spectra, they differ in significant ways. Firstly, as discussed in Chapter I, the former is less susceptible to oxidation, as evidenced by the reversibility of its oxidation by I2. Secondly, while the first excited triplet states (T₁) of both Pt(POP-BF₂) and [Pt₂(POP)₄]^4- exhibit long lifetimes (ca. 0.01 ms at room temperature) and substantial zero-field splitting (40 cm^-1), Pt(POP-BF₂) also has a remarkably long-lived (1.6 ns at room temperature) singlet excited state (S₁), indicating slow intersystem crossing (ISC). Fluorescence lifetime and quantum yield (QY) of Pt(POP-BF₂) were measured over a range of temperatures, providing insight into the slow ISC process. The remarkable spectroscopic and photophysical properties of Pt(POP-BF₂), both in solution and as a microcrystalline powder, form the theme of Chapter II.

In the second part of the thesis (Chapters III and IV), the electrochemical reduction of CO₂ to CO by [(L)Mn(CO)₃]^- catalysts is investigated using density functional theory (DFT). As discussed in Chapter III, the turnover frequency (TOF)-limiting step is the dehydroxylation of [(bpy)Mn(CO)₃(CO₂H)]^0/- (bpy = bipyridine) by trifluoroethanol (TFEH) to form [(bpy)Mn(CO)₄]^+/0. Because the dehydroxylation of [(bpy)Mn(CO)₃(CO₂H)]^- is faster, maximum TOF (TOF_max) is achieved at potentials sufficient to completely reduce [(bpy)Mn(CO)₃(CO₂H)]⁰ to [(bpy)Mn(CO)₃(CO₂H)]^-. Substitution of bipyridine with bipyrimidine reduces the overpotential needed, but at the expense of TOF_max. In Chapter IV, the decoration of the bipyrimidine ligand with a pendant alcohol is discussed as a strategy to increase CO₂ reduction activity. Our calculations predict that the pendant alcohol acts in concert with an external TFEH molecule, the latter acidifying the former, resulting in a ~ 80,000-fold improvement in the rate of TOF-limiting dehydroxylation of [(L)Mn(CO)₃(CO₂H)]^-.

An interesting strategy for the co-upgrading of light olefins and alkanes into heavier alkanes is the subject of Appendix B. The proposed scheme involves dimerization of the light olefin, operating in tandem with transfer hydrogenation between the olefin dimer and the light alkane. The work presented therein involved a Ta olefin dimerization catalyst and a silica-supported Ir transfer hydrogenation catalyst. Olefin dimer was formed under reaction conditions; however, this did not undergo transfer hydrogenation with the light alkane. A significant challenge is that the Ta catalyst selectively produces highly branched dimers, which are unable to undergo transfer hydrogenation.

Topological strings, double affine Hecke algebras, and exceptional knot homology

In this thesis, we consider two main subjects: refined, composite invariants and exceptional knot homologies of torus knots. The main technical tools are double affine Hecke algebras ("DAHA") and various insights from topological string theory.

In particular, we define and study the composite DAHA-superpolynomials of torus knots, which depend on pairs of Young diagrams and generalize the composite HOMFLY-PT polynomials from the full HOMFLY-PT skein of the annulus. We also describe a rich structure of differentials that act on homological knot invariants for exceptional groups. These follow from the physics of BPS states and the adjacencies/spectra of singularities associated with Landau-Ginzburg potentials. At the end, we construct two DAHA-hyperpolynomials which are closely related to the Deligne-Gross exceptional series of root systems.

In addition to these main themes, we also provide new results connecting DAHA-Jones polynomials to quantum torus knot invariants for Cartan types A and D, as well as the first appearance of quantum E6 knot invariants in the literature.

Investigations of DNA-mediated protein oxidation

DNA charge transport (CT) involves the efficient transfer of electrons or electron holes through the DNA π-stack over long molecular distances of at least 100 base-pairs. Despite this shallow distance dependence, DNA CT is sensitive to mismatches or lesions that disrupt π-stacking and is critically dependent on proper electronic coupling of the donor and acceptor moieties into the base stack. Favorable DNA CT is very rapid, occurring on the picosecond timescale. Because of this speed, electron holes equilibrate along the DNA π-stack, forming a characteristic pattern of DNA damage at low oxidation potential guanine multiplets. Furthermore, DNA CT may be used in a biological context. DNA processing enzymes with 4Fe4S clusters can perform DNA-mediated electron transfer (ET) self-exchange reactions with other 4Fe4S cluster proteins, even if the proteins are quite dissimilar, as long as the DNA-bound [4Fe4S]^3+/2+ redox potentials are conserved. This mechanism would allow low copy number DNA repair proteins to find their lesions efficiently within the cell. DNA CT may also be used biologically for the long-range, selective activation of redox-active transcription factors. Within this work, we pursue other proteins that may utilize DNA CT within the cell and further elucidate aspects of the DNA-mediated ET self-exchange reaction of 4Fe4S cluster proteins.

Dps proteins, bacterial mini-ferritins that protect DNA from oxidative stress, are implicated in the survival and virulence of pathogenic bacteria. One aspect of their protection involves ferroxidase activity, whereby ferrous iron is bound and oxidized selectively by hydrogen peroxide, thereby preventing formation of damaging hydroxyl radicals via Fenton chemistry. Understanding the specific mechanism by which Dps proteins protect the bacterial genome could inform the development of new antibiotics. We investigate whether DNA-binding E. coli Dps can utilize DNA CT to protect the genome from a distance. An intercalating ruthenium photooxidant was employed to generate oxidative DNA damage via the flash-quench technique, which localizes to a low potential guanine triplet. We find that Dps loaded with ferrous iron, in contrast to Apo-Dps and ferric iron-loaded Dps which lack available reducing equivalents, significantly attenuates the yield of oxidative DNA damage at the guanine triplet. These data demonstrate that ferrous iron-loaded Dps is selectively oxidized to fill guanine radical holes, thereby restoring the integrity of the DNA. Luminescence studies indicate no direct interaction between the ruthenium photooxidant and Dps, supporting the DNA-mediated oxidation of ferrous iron-loaded Dps. Thus DNA CT may be a mechanism by which Dps efficiently protects the genome of pathogenic bacteria from a distance.

Further work focused on spectroscopic characterization of the DNA-mediated oxidation of ferrous iron-loaded Dps. X-band EPR was used to monitor the oxidation of DNA-bound Dps after DNA photooxidation via the flash-quench technique. Upon irradiation with poly(dGdC)₂, a signal arises with g = 4.3, consistent with the formation of mononuclear high-spin Fe(III) sites of low symmetry, the expected oxidation product of Dps with one iron bound at each ferroxidase site. When poly(dGdC)₂ is substituted with poly(dAdT)₂, the yield of Dps oxidation is decreased significantly, indicating that guanine radicals facilitate Dps oxidation. The more favorable oxidation of Dps by guanine radicals supports the feasibility of a long-distance protection mechanism via DNA CT where Dps is oxidized to fill guanine radical holes in the bacterial genome produced by reactive oxygen species.

We have also explored possible electron transfer intermediates in the DNA-mediated oxidation of ferrous iron-loaded Dps. Dps proteins contain a conserved tryptophan residue in close proximity to the ferroxidase site (W52 in E. coli Dps). In comparison to WT Dps, in EPR studies of the oxidation of ferrous iron-loaded Dps following DNA photooxidation, W52Y and W52A mutants were deficient in forming the characteristic EPR signal at g = 4.3, with a larger deficiency for W52A compared to W52Y. In addition to EPR, we also probed the role of W52 Dps in cells using a hydrogen peroxide survival assay. Bacteria containing W52Y Dps survived the hydrogen peroxide challenge more similarly to those containing WT Dps, whereas cells with W52A Dps died off as quickly as cells without Dps. Overall, these results suggest the possibility of W52 as a CT hopping intermediate.

DNA-modified electrodes have become an essential tool for the study of the redox chemistry of DNA processing enzymes with 4Fe4S clusters. In many cases, it is necessary to investigate different complex samples and substrates in parallel in order to elucidate this chemistry. Therefore, we optimized and characterized a multiplexed electrochemical platform with the 4Fe4S cluster base excision repair glycosylase Endonuclease III (EndoIII). Closely packed DNA films, where the protein has limited surface accessibility, produce EndoIII electrochemical signals sensitive to an intervening mismatch, indicating a DNA-mediated process. Multiplexed analysis allowed more robust characterization of the CT-deficient Y82A EndoIII mutant, as well as comparison of a new family of mutations altering the electrostatics surrounding the 4Fe4S cluster in an effort to shift the reduction potential of the cluster. While little change in the DNA-bound midpoint potential was found for this family of mutants, likely indicating the dominant effect of DNA-binding on establishing the protein redox potential, significant variations in the efficiency of DNA-mediated electron transfer were apparent. On the basis of the stability of these proteins, examined by circular dichroism, we proposed that the electron transfer pathway in EndoIII can be perturbed not only by the removal of aromatic residues but also through changes in solvation near the cluster.

While the 4Fe4S cluster of EndoIII is relatively insensitive to oxidation and reduction in solution, we have found that upon DNA binding, the reduction potential of the [4Fe4S]^3+/2+ couple shifts negatively by approximately 200 mV, bringing this couple into a physiologically relevant range. Demonstrated using electrochemistry experiments in the presence and absence of DNA, these studies do not provide direct molecular evidence for the species being observed. Sulfur K-edge X-ray absorbance spectroscopy (XAS) can be used to probe directly the covalency of iron-sulfur clusters, which is correlated to their reduction potential. We have shown that the Fe-S covalency of the 4Fe4S cluster of EndoIII increases upon DNA binding, stabilizing the oxidized [4Fe4S]³⁺ cluster, consistent with a negative shift in reduction potential. The 7% increase in Fe-S covalency corresponds to an approximately 150 mV shift, remarkably similar to DNA electrochemistry results. Therefore we have obtained direct molecular evidence for the shift in 4Fe4S reduction potential of EndoIII upon DNA binding, supporting the feasibility of our model whereby these proteins can utilize DNA CT to cooperate in order to efficiently find DNA lesions inside cells.

In conclusion, in this work we have explored the biological applications of DNA CT. We discovered that the DNA-binding bacterial ferritin Dps can protect the bacterial genome from a distance via DNA CT, perhaps contributing to pathogen survival and virulence. Furthermore, we optimized a multiplexed electrochemical platform for the study of the redox chemistry of DNA-bound 4Fe4S cluster proteins. Finally, we have used sulfur K-edge XAS to obtain direct molecular evidence for the negative shift in 4Fe4S cluster reduction potential of EndoIII upon DNA binding. These studies contribute to the understanding of DNA-mediated protein oxidation within cells.

Electrocatalytic reduction of dioxygen at chemically modified electrodes

The kinetics of the reduction of O₂ by Ru(NH₃)₆⁺² as catalyzed by cobalt(II) tetrakis(4-N-methylpyridyl)porphyrin are described both in homogeneous solution and when the reactants are confined to Nafion coatings on graphite electrodes. The catalytic mechanism is determined and the factors that can control the total reduction currents at Nafion-coated electrodes are specified. A kinetic zone diagram for analyzing the behavior of catalyst-mediator-substrate systems at polymer coated electrodes is presented and utilized in identifying the current-limiting processes. Good agreement is demonstrated between calculated and measured reduction currents at rotating disk electrodes. The experimental conditions that will yield the optimum performance of coated electrodes are discussed, and a relationship is derived for the optimal coating thickness.

The relation between the reduction potentials of adsorbed and unadsorbed cobalt(III) tetrakis(4-N-methylpyridyl)porphyrin and those where it catalyzes the electroreduction of dioxygen is described. There is an unusually large change in the formal potential of the Co(III) couple upon the adsorption of the porphyrin on the graphite electrode surface. The mechanism in which the (inevitably) adsorbed porphyrin catalyzes the reduction of O₂ is in accord with a general mechanistic scheme proposed for most monomeric cobalt porphyrins.

Four new dimeric metalloporphyrins (prepared in the laboratory of Professor C. K. Chang) have the two porphyrin rings linked by an anthracene bridge attached to meso positions. The electrocatalytic behavior of the diporphyrins towards the reduction of O₂ at graphite electrodes has been examined for the following combination of metal centers: Co-Cu, Co-Fe, Fe-Fe, Fe-H₂. The Co-Cu diporphyrin catalyzes the reduction of O₂ to H₂O₂ but no further. The other three catalysts all exhibit mixed reduction pathways leading to both H₂O₂ and H₂O. However, the pathways that lead to H₂O do not involve H₂O₂ as an intermediate. A possible mechanistic scheme is offered to account for the observed behavior.