Solid-state optical absorption from optimally tuned time-dependent range-separated hybrid density functional theory

We present a framework for obtaining reliable solid-state charge and optical excitations and spectra from optimally tuned range-separated hybrid density functional theory. The approach, which is fully couched within the formal framework of generalized Kohn-Sham theory, allows for the accurate prediction of exciton binding energies. We demonstrate our approach through first principles calculations of one- and two-particle excitations in pentacene, a molecular semiconducting crystal, where our work is in excellent agreement with experiments and prior computations. We further show that with one adjustable parameter, set to produce the known band gap, this method accurately predicts band structures and optical spectra of silicon and lithium fluoride, prototypical covalent and ionic solids. Our findings indicate that for a broad range of extended bulk systems, this method may provide a computationally inexpensive alternative to many-body perturbation theory, opening the door to studies of materials of increasing size and complexity.

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.

Carbon monoxide adsorption on molybdenum phosphides: Fourier transform infrared spectroscopic and density functional theory studies

Molybdenum phosphide (MoP) and supported molybdenum phosphide (MoP/gamma-Al2O3) have been prepared by the temperature-programmed reduction method. The surface sites of the MoP/gamma-Al2O3 catalyst were characterized by carbon monoxide (CO) adsorption with in situ Fourier transform infrared (FT-IR) spectroscopy. A characteristic IR band at 2037 cm(-1) was observed on the MoP/gamma-Al2O3 that was reduced at 973 K. This band is attributed to linearly adsorbed CO on Mo atoms of the MoP surface and is similar to IR bands at 2040-2060 cm(-1), which correspond to CO that has been adsorbed on some noble metals, such as platinum, palladium, and rhodium. Density functional calculations of the structure of molybdenum phosphides, as well as CO chemisorption on the MoP(001) surface, have also been studied on periodic surface models, using the generalized gradient approximation (GGA) for the exchange-correlation functional. The results show that the chemisorption of CO on MoP occurred mainly on top of molybdenum, because the bonding of CO requires a localized mininum potential energy. The adsorption energy obtained is DeltaH(ads) approximate to -2.18 eV, and the vibrational frequency of CO is 2047 cm-1, which is in good agreement with the IR result of CO chernisorption on MoP/gamma-Al2O3.

Effect of secondary ligands' size on energy transfer and electroluminescent efficiencies for a series of europium(III) complexes, a density functional theory study

In this paper, a quantum chemistry method was used to investigate the effect of different sizes of substituted phenanthrolines on absorption, energy transfer, and the electroluminescent performance of a series of Eu(TTA)(3)L (L = [1,10] phenanthroline (Phen), Pyrazino[2,3-f][1,10]phenanthroline (PyPhen), 2-methylprrazino[2,3-f][1,10] phenanthroline(MPP), dipyrido[3,2-a:2',3'-c]phenazine(DPPz), 11-methyldipyrido[3,2-a:2',3'c]phenazine(MDPz), 11.12-dimethyldipyrido[3,2-a:2',3'-c]phenazine(DDPz), and benzo[i]dipyrido[3,2-a:2',3'-c]phenazine (BDPz)) complexes. Absorption spectra calculations show that different sizes of secondary ligands have different effects on transition characters, intensities, and absorption peak positions.

Reformulating time-dependent density functional theory with non-orthogonal localized molecular orbitals.

Time-dependent density functional theory (TDDFT) has broad application in the study of electronic response, excitation and transport. To extend such application to large and complex systems, we develop a reformulation of TDDFT equations in terms of non-orthogonal localized molecular orbitals (NOLMOs). NOLMO is the most localized representation of electronic degrees of freedom and has been used in ground state calculations. In atomic orbital (AO) representation, the sparsity of NOLMO is transferred to the coefficient matrix of molecular orbitals (MOs). Its novel use in TDDFT here leads to a very simple form of time propagation equations which can be solved with linear-scaling effort. We have tested the method for several long-chain saturated and conjugated molecular systems within the self-consistent charge density-functional tight-binding method (SCC-DFTB) and demonstrated its accuracy. This opens up pathways for TDDFT applications to large bio- and nano-systems.

Experimental and density functional theory study of the vibrational properties of 2-mercaptobenzimidazole in interaction with gold

The vibrational properties of the 2-mercaptobenzimidazole (MBI) molecule in interaction with gold were examined by a combined approach of FTIR measurements and density functional theory (DFT). A complete assignment of the 42 normal modes of MBI has been performed on the basis of DFT calculations at the B3PW91 level in complement to the Raman and FTIR spectra. Calculations demonstrated that, on the deprotonated MBI molecule, the negative charge is localized on the sulfur atom, favoring the formation of a gold-sulfur bond upon reaction of MBI with gold. This was confirmed by the very good agreement between the calculated spectrum and the experimental spectra of different gold-MBI compounds, indicating that the vibrational properties of adsorbed MBI are chiefly determined by the coordination through the sulfur atom. © 2006 American Chemical Society.

General rules for predicting where a catalytic reaction should occur on metal surfaces: A density functional theory study of C-H and C-O bond breaking/making on flat, stepped, and kinked metal surfaces

To predict where a catalytic reaction should occur is a fundamental issue scientifically. Technologically, it is also important because it can facilitate the catalyst's design. However, to date, the understanding of this issue is rather limited. In this work, two types of reactions, CH4 CH3 + H and CO C + 0 on two transition metal surfaces, were chosen as model systems aiming to address in general where a catalytic reaction should occur. The dissociations of CH4 - CH3 + H and CO --> C + O and their reverse reactions on flat, stepped, and kinked Rh and Pd surfaces were studied in detail. We find the following: First, for the CH4 Ch(3) + H reaction, the dissociation barrier is reduced by similar to0.3 eV on steps and kinks as compared to that on flat surfaces. On the other hand, there is essentially no difference in barrier for the association reaction of CH3 + H on the flat surfaces and the defects. Second, for the CO C + 0 reaction, the dissociation barrier decreases dramatically (more than 0.8 eV on Rh and Pd) on steps and kinks as compared to that on flat surfaces. In contrast to the CH3 + H reaction, the C + 0 association reaction also preferentially occurs on steps and kinks. We also present a detailed analysis of the reaction barriers in which each barrier is decomposed quantitatively into a local electronic effect and a geometrical effect. Our DFT calculations show that surface defects such as steps and kinks can largely facilitate bond breaking, while whether the surface defects could promote bond formation depends on the individual reaction as well as the particular metal. The physical origin of these trends is identified and discussed. On the basis of our results, we arrive at some simple rules with respect to where a reaction should occur: (i) defects such as steps are always favored for dissociation reactions as compared to flat surfaces; and (ii) the reaction site of the association reactions is largely related to the magnitude of the bonding competition effect, which is determined by the reactant and metal valency. Reactions with high valency reactants are more likely to occur on defects (more structure-sensitive), as compared to reactions with low valency reactants. Moreover, the reactions on late transition metals are more likely to proceed on defects than those on the early transition metals.

A Density Functional Theory Study on the Active Center of Fe-Only Hydrogenase: Characterization and Electronic Structure of the Redox States

We have carried out extensive density functional theory (DFT) calculations for possible redox states of the active center in Fe-only hydrogenases. The active center is modeled by [(H(CH(3))S)(CO)(CN(-))Fe(p)(mu-DTN)(mu-CO)Fe(d)(CO)(CN(-))(L)](z) (z is the net charge in the complex; Fe(p)= the proximal Fe, Fe(d) = the distal Fe, DTN = (-SCH(2)NHCH(2)S-), L is the ligand that bonds with the Fed at the trans position to the bridging CO). Structures of possible redox states are optimized, and CO stretching frequencies are calculated. By a detailed comparison of all the calculated structures and the vibrational frequencies with the available experimental data, we find that (i) the fully oxidized, inactive state is an Fe(II)-Fe(II) state with a hydroxyl (OH(-)) group bonded at the Fe(d), (ii) the oxidized, active state is an Fe(II)-Fe(l) complex which is consistent with the assignment of Cao and Hall (J. Am. Chem. Soc. 2001, 123, 3734), and (iii) the fully reduced state is a mixture with the major component being a protonated Fe(l)-Fe(l) complex and the other component being its self-arranged form, Fe(II)-Fe(II) hydride, Our calculations also show that the exogenous CO can strongly bond with the Fe(II)-Fe(l) species, but cannot bond with the Fe(l)-Fe(l) complex. This result is consistent with experiments that CO tends to inhibit the oxidized, active state, but not the fully reduced state. The electronic structures of all the redox states have been analyzed. It is found that a frontier orbital which is a mixing state between the e(g) of Fe and the 2pi of the bridging CO plays a key role concerning the reactivity of Fe-only hydrogenases: (1) it is unoccupied in the fully oxidized, inactive state, half-occupied in the oxidized, active state, and fully occupied in the fully reduced state; (ii) the e(g)-2pi orbital is a bonding state, and this is the key reason for stability of the low oxidation states, such as Fe(l)-Fe(l) complexes; and (iii) in the e(g)-2pi orbital more charge accumulates between the bridging CO and the Fe(d) than between the bridging CO and the Fe(p), and the occupation increase in this orbital will enhance the bonding between the bridging CO and the Fe(d), leading to the bridging-CO shift toward the Fe(d).

Insight into the solvent effect: A density functional theory study of cisplatin hydrolysis

The solvent effect on reactions in solutions is crucial for many systems. In this study, the reaction barrier with respect to the number of solvent molecules included in the system is systematically studied using density function theory calculations. Our results show that the barriers rapidly converge with respect to the number of solvent molecules. The solvent effect is investigated by calculating cisplatin hydrolysis in several types of solvents. The results are analyzed and a linear relationship between the reaction barrier and the interaction strength of solvent-reactants is found. Insight into the general solvent effect is obtained. (c) 2006 American Institute of Physics.

A density functional theory study of the a-olefin selectivity in Fischer-Tropsch synthesis

We studied the alpha-olefin selectivity in Fischer-Tropsch (FT) synthesis using density functional theory (131717) calculations. We calculated the relevant elementary steps from C-2 to C-6 species. Our results showed that the barriers of hydrogenation and dehydrogenation reactions were constant with different chain lengths, and the chemisorption energies of alpha-olefins from DFT calculations also were very similar, except for C-2 species. A simple expression of the paraffin/olefin ratio was obtained based on a kinetic model. Combining the expression of the paraffin/olefin ratio and our calculation results, experimental findings are satisfactorily explained. We found that the physical origin of the chain length dependence of paraffin/olefin ratio is the chain length dependence of both the van der Waals interaction between adsorbed alpha-olefins and metal surfaces and the entropy difference between adsorbed and gaseous alpha-olefins, and that the greater chemisorption energy of ethylene is the main reason for the abnormal ethane/ethylene ratio. (c) 2008 Elsevier Inc. All rights reserved.

A density functional theory study of the surface relaxation and reactivity of Cu2O(100)

Density functional theory has been used to investigate the surface relaxation of Cu2O(100) and the adsorption of NO. The calculations indicate the formation of surface copper dimers on relaxation coupled with a large contraction of the spacing between the first and second layers. Local density of states for atoms in the top three layers shows that the third layer copper atoms have the greatest change in bonding character. Adsorption energies have been calculated for the N-down and O-down adsorption of NO on the Cu2O(100) surface. These indicate that N-down adsorption is favoured and that in this case NO-lattice oxygen interactions dominate the adsorbate structure. (C) 2000 Elsevier Science B.V. All rights reserved.

Reaction Mechanisms of Crotonaldehyde Hydrogenation on Pt(111): Density Functional Theory and Microkinetic Modeling

The microkinetics based on density function theory (DFT) calculations is utilized to investigate the reaction mechanism of crotonaldehyde hydrogenation on Pt(111) in the free energy landscape. The dominant reaction channel of each hydrogenation product is identified. Each of them begins with the first surface hydrogenation of the carbonyl oxygen of crotonaldehyde on the surface. A new mechanism, 1,4-addition mechanism generating enols (butenol), which readily tautomerize to saturated aldehydes (butanal), is identified as a primary mechanism to yield saturated aldehydes instead of the 3,4-addition via direct hydrogenation of the ethylenic bond. The calculation results also show that the full hydrogenation product, butylalcohol, mainly stems from the deep hydrogenation of surface open-shell dihydrogenation intermediates. It is found that the apparent barriers of the dominant pathways to yield three final products are similar on P(111), which makes it difficult to achieve a high selectivity to the desired crotyl alcohol (COL).