Density functional theory embedding for correlated wavefunctions

Methods that exploit the intrinsic locality of molecular interactions show significant promise in making tractable the electronic structure calculation of large-scale systems. In particular, embedded density functional theory (e-DFT) offers a formally exact approach to electronic structure calculations in which the interactions between subsystems are evaluated in terms of their electronic density. In the following dissertation, methodological advances of embedded density functional theory are described, numerically tested, and applied to real chemical systems.

First, we describe an e-DFT protocol in which the non-additive kinetic energy component of the embedding potential is treated exactly. Then, we present a general implementation of the exact calculation of the non-additive kinetic potential (NAKP) and apply it to molecular systems. We demonstrate that the implementation using the exact NAKP is in excellent agreement with reference Kohn-Sham calculations, whereas the approximate functionals lead to qualitative failures in the calculated energies and equilibrium structures.

Next, we introduce density-embedding techniques to enable the accurate and stable calculation of correlated wavefunction (CW) in complex environments. Embedding potentials calculated using e-DFT introduce the effect of the environment on a subsystem for CW calculations (WFT-in-DFT). We demonstrate that WFT-in-DFT calculations are in good agreement with CW calculations performed on the full complex.

We significantly improve the numerics of the algorithm by enforcing orthogonality between subsystems by introduction of a projection operator. Utilizing the projection-based embedding scheme, we rigorously analyze the sources of error in quantum embedding calculations in which an active subsystem is treated using CWs, and the remainder using density functional theory. We show that the embedding potential felt by the electrons in the active subsystem makes only a small contribution to the error of the method, whereas the error in the nonadditive exchange-correlation energy dominates. We develop an algorithm which corrects this term and demonstrate the accuracy of this corrected embedding scheme.

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.

A variational framework for spectral discretization of the density matrix in Kohn Sham density functional theory

Kohn-Sham density functional theory (KSDFT) is currently the main work-horse of quantum mechanical calculations in physics, chemistry, and materials science. From a mechanical engineering perspective, we are interested in studying the role of defects in the mechanical properties in materials. In real materials, defects are typically found at very small concentrations e.g., vacancies occur at parts per million, dislocation density in metals ranges from $10^{10} m^{-2}$ to $10^{15} m^{-2}$, and grain sizes vary from nanometers to micrometers in polycrystalline materials, etc. In order to model materials at realistic defect concentrations using DFT, we would need to work with system sizes beyond millions of atoms. Due to the cubic-scaling computational cost with respect to the number of atoms in conventional DFT implementations, such system sizes are unreachable. Since the early 1990s, there has been a huge interest in developing DFT implementations that have linear-scaling computational cost. A promising approach to achieving linear-scaling cost is to approximate the density matrix in KSDFT. The focus of this thesis is to provide a firm mathematical framework to study the convergence of these approximations. We reformulate the Kohn-Sham density functional theory as a nested variational problem in the density matrix, the electrostatic potential, and a field dual to the electron density. The corresponding functional is linear in the density matrix and thus amenable to spectral representation. Based on this reformulation, we introduce a new approximation scheme, called spectral binning, which does not require smoothing of the occupancy function and thus applies at arbitrarily low temperatures. We proof convergence of the approximate solutions with respect to spectral binning and with respect to an additional spatial discretization of the domain. For a standard one-dimensional benchmark problem, we present numerical experiments for which spectral binning exhibits excellent convergence characteristics and outperforms other linear-scaling methods.

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.

Structure, chemistry, and energetics of organic and inorganic adsorbates on Ga-rich GaAs and GaP(001) surfaces

The work described in this dissertation includes fundamental investigations into three surface processes, namely inorganic film growth, water-induced oxidation, and organic functionalization/passivation, on the GaP and GaAs(001) surfaces. The techniques used to carry out this work include scanning tunneling microscopy (STM), X-ray photoelectron spectroscopy (XPS), and density functional theory (DFT) calculations. Atomic structure, electronic structure, reaction mechanisms, and energetics related to these surface processes are discussed at atomic or molecular levels.

First, we investigate epitaxial Zn₃P₂ films grown on the Ga-rich GaAs(001)(6×6) surface. The film growth mechanism, electronic properties, and atomic structure of the Zn₃P₂/GaAs(001) system are discussed based on experimental and theoretical observations. We discover that a P-rich amorphous layer covers the crystalline Zn3P2 film during and after growth. We also propose more accurate picture of the GaP interfacial layer between Zn₃P₂ and GaAs, based on the atomic structure, chemical bonding, band diagram, and P-replacement energetics, than was previously anticipated.

Second, DFT calculations are carried out in order to understand water-induced oxidation mechanisms on the Ga-rich GaP(001)(2×4) surface. Structural and energetic information of every step in the gaseous water-induced GaP oxidation reactions are elucidated at the atomic level in great detail. We explore all reasonable ground states involved in most of the possible adsorption and decomposition pathways. We also investigate structures and energies of the transition states in the first hydrogen dissociation of a water molecule on the (2×4) surface.

Finally, adsorption structures and thermal decomposition reactions of 1-propanethiol on the Ga-rich GaP(001)(2×4) surface are investigated using high resolution STM, XPS, and DFT simulations. We elucidate adsorption locations and their associated atomic structures of a single 1-propanethiol molecule on the (2×4) surface as a function of annealing temperature. DFT calculations are carried out to optimize ground state structures and search transition states. XPS is used to investigate variations of the chemical bonding nature and coverage of the adsorbate species.

Chain-forming zintl antimonides as novel thermoelectric materials

Zintl phases, a subset of intermetallic compounds characterized by covalently-bonded "sub-structures," surrounded by highly electropositive cations, exhibit precisely the characteristics desired for thermoelectric applications. The requirement that Zintl compounds satisfy the valence of anions through the formation of covalent substructures leads to many unique, complex crystal structures. Such complexity often leads to exceptionally low lattice thermal conductivity due to the containment of heat in low velocity optical modes in the phonon dispersion. To date, excellent thermoelectric properties have been demonstrated in several Zintl compounds. However, compared with the large number of known Zintl phases, very few have been investigated as thermoelectric materials.

From this pool of uninvestigated compounds, we selected a class of Zintl antimonides that share a common structural motif: anionic moieties resembling infinite chains of linked MSb₄ tetrahedra, where $M$ is a triel element. The compounds discussed in this thesis (A₅M₂Sb₆ and A₃MSb₃, where A = Ca or Sr and M = Al, Ga and In) crystallize as four distinct, but closely related "chain-forming" structure types. This thesis describes the thermoelectric characterization and optimization of these phases, and explores the influence of their chemistry and structure on the thermal and electronic transport properties. Due to their large unit cells, each compound exhibits exceptionally low lattice thermal conductivity (0.4 - 0.6 W/mK at 1000 K), approaching the predicted glassy minimum at high temperatures. A combination of Density Functional calculations and classical transport models were used to explain the experimentally observed electronic transport properties of each compound. Consistent with the Zintl electron counting formalism, A₅M₂Sb₆ and A₃MSb₃ phases were found to have filled valence bands and exhibit intrinsic electronic properties. Doping with divalent transition metals (Zn²⁺ and Mn²⁺) on the M³⁺ site, or Na¹⁺ on the A³⁺ site allowed for rational control of the carrier concentration and a transition towards degenerate semiconducting behavior. In optimally-doped samples, promising peak zT values between 0.4 and 0.9 were obtained, highlighting the value of continued investigations of complex Zintl phases.

Experimental chemical physics of droplets and clusters

Time-of-flight measurements of energetic He atoms, field ionization of cryogenic liquid helium clusters, and time-of-flight and REMPI spectroscopy of radical salt clusters were investigated experimentally. The excited He atoms were generated in a corona discharge. Two strong neutral peaks were observed, accompanied by a prompt photon peak and a charged peak. All peaks were correlated with the pulsing of the discharge. The neutral hyperthermal and metastable atoms were formed by different mechanisms at different stages of the corona discharge. Positively charged helium droplets were produced by ionization of liquid helium in an electrostatic spraying experiment. The fluid emerging from a thin glass capillary was ionized by a high voltage applied to a needle inside the capillary. Fine droplets (less than 10 µm in diameter) were produced in showers with currents as high as 0.4 µA at 2-4 kV. The high currents resulting from field ionization in helium and the low surface tension of He I, led to charge densities that greatly exceeded the Rayleigh limit, thus resulting in coulombic explosion of the liquid. In contrast, liquid nitrogen formed a well-defined Taylor cone with droplets having diameters comparable to the jet (≈100 µm) at lower currents (10 nA) and higher voltages (8 kV). The metal-halide clusters of calcium and chlorine were generated by laser ablation of calcium metal in a Ar/CCl₄ expansion. A visible spectrum of the Ca₂Cl₃ cluster was observed from 651 to 630 nm by 1 +1' REMPI. The spectra were composed of a strong origin band at 15 350.8 cm^-1 and several weak vibronic bands. Density functional calculations predicted three minimum energy isomers. The spectrum was assigned to the ²B₂ ← X ²A₁ transition of a planar C_2V structure having a ring of two Cl and two Ca atoms and a terminal Cl atom. The ring isomer of Ca₂Cl₃ has the unpaired electron localized on one Ca²⁺ ion to form a Ca⁺ chromophore. A second electronic band of Ca₂Cl₃ was observed at 720 nm. The band is sharply different from the 650 nm band and likely due to a different isomer.

Zirconocenes as models for homogeneous Ziegler-Natta olefin polymerization catalysts

Using density functional theory, we studied the fundamental steps of olefin polymerization for zwitterionic and cationic Group IV ansa-zirconocenes and a neutral ansa- yttrocene. Complexes [H₂E(C₅H₄)₂ZrMe]ⁿ (n = 0: E = BH₂ (1), BF₂ (2), AlH₂(3); n = +: E = CH₂(4), SiH₂(5)) and H₂Si(C₅H₄)₂YMe were used as computational models. The largest differences among these three classes of compounds were the strength of olefin binding and the stability of the β-agostic alkyl intermediate towards β-hydrogen elimination. We investigated the effect of solvent on the reaction energetics for land 5. We found that in benzene the energetics became very similar except that a higher olefin insertion barrier was calculated for 1. The calculated anion affinity of [CH₃BF₃]^- was weaker towards 1 than 5. The calculated olefin binding depended primarily on the charge of the ansa linker, and the olefin insertion barrier was found to decrease steadily in the following order: [H₂C(C₅H₄)₂ZrMe]⁺ > [F₂B(C₅H₄)₂ZrMe] ≈ [H₂B(C₅H₄)₂ZrMe] > [H₂Si(C₅H₄)₂ZrMe]⁺ > [H₂Al(C₅H₄)₂ZrMe].

We prepared ansa-zirconocene dicarbonyl complexes Me₂ECp₂Zr(CO)₂ (E = Si, C), and t-butyl substituted complexes (t-BuCp)₂Zr(CO)₂, Me₂E(t-BuCp)₂Zr(CO)₂ (E = Si, C), (Me₂Si)₂(t-BuCp)₂Zr(CO)₂ as well as analogous zirconocene complexes. Both the reduction potentials and carbonyl stretching frequencies follow the same order: Me₂SiCp₂ZrCl₂> Me₂CCp₂ZrCl₂> Cp₂ZrCl₂> (Me₂Si)₂Cp₂ZrCl₂. This ordering is a result of both the donating abilities of the cyclopentadienyl substituents and the orientation of the cyclopentadiene rings. Additionally, we prepared a series of analogous cationic zirconocene complexes [LZrOCMe₃][MeB(C₆F₅)₃] (L = CP₂, Me₂SiCp₂, Me₂CCP₂, (Me₂Si)₂Cp₂) and studied the kinetics of anion dissociation. We found that the enthalpy of anion dissociation increased from 10.3 kcal•mol^-1 to 17.6 kcal•mol^-1 as exposure of the zirconium center increased.

We also prepared series of zirconocene complexes bearing 2,2-dimethyl-2-sila-4-pentenyl substituents (and methyl-substituted olefin variants). Methide abstraction with B(C₆F₅) results in reversible coordination of the tethered olefin to the cationic zirconium center. The kinetics of olefin dissociation have been examined using NMR methods, and the effects of ligand variation for unlinked, singly [SiMe₂]-linked and doubly [SiMe₂]-linked bis(cyclopentadienyl) arrangements has been compared (ΔG‡ for olefin dissociation varies from 12.8 to 15.6 kcal•mol^-1). Methide abstraction from 1,2-(SiMe₂)₂(η⁵-C₅H₃)₂Zr(CH₃)-(CH₂CMe₂CH₂CH = CH₂) results in rapid β-allyl elimination with loss of isobutene yielding the allyl cation [{1,2-(SiMe₂)₂(η⁵-C₅H₃)₂Zr(η³-CH₂CH=CH₂)]⁺.

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.

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.