Excited state properties of formamide in water solution: An ab initio study

We present ab initio quantum calculation of the optical properties of formamide in vapor phase and in water solution. We employ time dependent density functional theory for the isolated molecule and many-body perturbation theory methods for the system in solution. An average over several molecular dynamics snapshots is performed to take into account the disorder of the liquid. We find that the excited stateproperties of the gas-phase formamide are strongly modified by the presence of the water solvent: the geometry of the molecule is distorted and the electronic and optical properties are severely modified. The important interaction among the formamide and the water molecules forces us to use fully quantum methods for the calculation of the excited stateproperties of this system. The excitonic wave function is localized both on the solute and on part of the solvent.

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

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

Towards Accurate and Efficient Description of Excited States

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

Ground and Electronic Excited States from Pairing Matrix Fluctuation and Particle-Particle Random Phase Approximation

The accurate description of ground and electronic excited states is an important and challenging topic in quantum chemistry. The pairing matrix fluctuation, as a counterpart of the density fluctuation, is applied to this topic. From the pairing matrix fluctuation, the exact electron correlation energy as well as two electron addition/removal energies can be extracted. Therefore, both ground state and excited states energies can be obtained and they are in principle exact with a complete knowledge of the pairing matrix fluctuation. In practice, considering the exact pairing matrix fluctuation is unknown, we adopt its simple approximation --- the particle-particle random phase approximation (pp-RPA) --- for ground and excited states calculations. The algorithms for accelerating the pp-RPA calculation, including spin separation, spin adaptation, as well as an iterative Davidson method, are developed. For ground states correlation descriptions, the results obtained from pp-RPA are usually comparable to and can be more accurate than those from traditional particle-hole random phase approximation (ph-RPA). For excited states, the pp-RPA is able to describe double, Rydberg, and charge transfer excitations, which are challenging for conventional time-dependent density functional theory (TDDFT). Although the pp-RPA intrinsically cannot describe those excitations excited from the orbitals below the highest occupied molecular orbital (HOMO), its performances on those single excitations that can be captured are comparable to TDDFT. The pp-RPA for excitation calculation is further applied to challenging diradical problems and is used to unveil the nature of the ground and electronic excited states of higher acenes. The pp-RPA and the corresponding Tamm-Dancoff approximation (pp-TDA) are also applied to conical intersections, an important concept in nonadiabatic dynamics. Their good description of the double-cone feature of conical intersections is in sharp contrast to the failure of TDDFT. All in all, the pairing matrix fluctuation opens up new channel of thinking for quantum chemistry, and the pp-RPA is a promising method in describing ground and electronic excited states.

Ultrafast charge-transfer in organic photovoltaic interfaces: geometrical and functionalization effects

Understanding the microscopic mechanisms of electronic excitation in organic photovoltaic cells is a challenging problem in the design of efficient devices capable of performing sunlight harvesting. Here we develop and apply an ab initio approach based on time-dependent density functional theory and Ehrenfest dynamics to investigate photoinduced charge transfer in small organic molecules. Our calculations include mixed quantum–classical dynamics with ions moving classically and electrons quantum mechanically, where no experimental external parameter other than the material geometry is required. We show that the behavior of photocarriers in zinc phthalocyanine (ZnPc) and C60 systems, an effective prototype system for organic solar cells, is sensitive to the atomic orientation of the donor and the acceptor units as well as the functionalization of covalent molecules at the interface. In particular, configurations with the ZnPc molecules facing on C60 facilitate charge transfer between substrate and molecules that occurs within 200 fs. In contrast, configurations where ZnPc is tilted above C60 present extremely low carrier injection efficiency even at longer times as an effect of the larger interfacial potential level offset and higher energetic barrier between the donor and acceptor molecules. An enhancement of charge injection into C60 at shorter times is observed as binding groups connect ZnPc and C60 in a dyad system. Our results demonstrate a promising way of designing and controlling photoinduced charge transfer on the atomic level in organic devices that would lead to efficient carrier separation and maximize device performance.

Thiolated Au18cluster: preferred Ag sites for doping, structures, and optical and chiroptical properties

Recently, the X-ray determined structure of the thiolated Au18 cluster has been reported. In this communication, we addressed a study of structures and chiroptical properties of thiolated Au18 cluster doped with up to ten Ag atoms, which have been calculated by Time Dependent Density Functional Theory (TD-DFT). The number of Ag atoms was steadily varied and more stable isomers showed optical and Circular Dichroism (CD) spectra distinct from that found for the parent Au18 cluster. Doping with more than four Ag atoms results in enhancement of the oscillator strength of the HOMO–LUMO peak and it is expected that this feature can be exploited for photoluminescence applications.

Theoretical Rationalization of the Singlet–Triplet Gap in OLEDs Materials: Impact of Charge-Transfer Character

New materials for OLED applications with low singlet–triplet energy splitting have been recently synthesized in order to allow for the conversion of triplet into singlet excitons (emitting light) via a Thermally Activated Delayed Fluorescence (TADF) process, which involves excited-states with a non-negligible amount of Charge-Transfer (CT). The accurate modeling of these states with Time-Dependent Density Functional Theory (TD-DFT), the most used method so far because of the favorable trade-off between accuracy and computational cost, is however particularly challenging. We carefully address this issue here by considering materials with small (high) singlet–triplet gap acting as emitter (host) in OLEDs and by comparing the accuracy of TD-DFT and the corresponding Tamm-Dancoff Approximation (TDA), which is found to greatly reduce error bars with respect to experiments thanks to better estimates for the lowest singlet–triplet transition. Finally, we quantitatively correlate the singlet–triplet splitting values with the extent of CT, using for it a simple metric extracted from calculations with double-hybrid functionals, that might be applied in further molecular engineering studies.

New indolium fluorofores for two-photon absortion (2PA) bioimaging applications

The development of organic materials with 2PA has attracted intensive attention in the past two decades [1]. In two-photon bio-imaging applications the design of the chromophore requires to have a good cross-section (σ2PA) and good biological compatibility which depends on the molecular volume and polarity [2]. In this work, we present the design, synthesis and characterization of new indolium derivatives. These compounds are easy to achieve with good yields and good photophysical properties. In addition, time-dependent density functional theory (TDDFT) has been carried out to investigate the energy level of the ground and excited state. Their spectral properties and assays performed on cultured cells, demonstrate the potential of these compounds as fluorescent probes with application in two-photon bio-imaging.

Selective Sonogashira Couplings in 1,2,4,5-tetrazine derivatives

1,2,4,5-Tetrazines are six-membered heterocyclic compounds in which the four nitrogen atoms are displayed in a symmetric fashion. Their reactivity is quite different from other heterocyclic aromatic systems due to its unique electron-withdrawing character, comparable to tetra-nitrobenzene. 1 In particular, 1,2,4,5- tetrazines are known to take part in [4+2] inverse-Diels–Alder cycloaddition processes which efficiently lead to the construction of substituted pyridazine systems that are important in drug development and biomarker applications. 2 However, the electronic character of 1,2,4,5-tetrazines hampered the development of 3- ethynyl- and 3,6-diethynyl-1,2,4,5-tetrazine derivatives for molecular electronic applications, proved by the scarcity of examples found in the literature. 3 Herein, we describe the synthesis and characterization of two novel ethynyl-based 1,2,4,5-tetrazine derivatives. Synthesis of 3,6-(4-bromophenyl)-1,2,4,5-tetrazine precursor (1) was achieved in good yield by Pinner’s method, starting from 4-bromobenzonitrile. Despite its low solubility in common organic solvents, this precursor was found to react smoothly under typical Sonogashira coupling conditions to selectively afford the 3-ethynyl (2) and 3,6-diethynyl (3) protected derivatives (Figure 1). Reaction conditions were evaluated in order to provide the best yields and to promote selectivity of the mono- or disubstituted ethynyl derivatives. Finally, deprotection was achieved affording, in the case of compound 3, an unprecedented 3,6- diethynyl-1,2,4,5-tetrazine compound. Time-Dependent Density Functional Theory (TDDFT) calculations for both deprotected ethynyl derivatives were used to simulate electronic spectra. A deep knowledge of the relevant electronic transitions involved and quantitatively satisfactory results of the calculated electronic excitations in comparison with experimental data were obtained.

Spin and density longitudinal response of quantum dots in the time-dependent local-spin-density approximation

The longitudinal dipole response of a quantum dot has been calculated in the far-infrared regime using local-spin-density-functional theory. We have studied the coupling between the collective spin and density modes as a function of the magnetic field. We have found that the spin dipole mode and single-particle excitations have a sizable overlap, and that the magnetoplasmon modes can be excited by the dipole spin operator if the dot is spin polarized. The frequency of the dipole spin edge mode presents an oscillation which is clearly filling factor (v) related. We have found that the spin dipole mode is especially soft for even-n values. Results for selected numbers of electrons and confining potentials are discussed.

Spin and density longitudinal response of quantum dots in the time-dependent local-spin-density approximation

The longitudinal dipole response of a quantum dot has been calculated in the far-infrared regime using local-spin-density-functional theory. We have studied the coupling between the collective spin and density modes as a function of the magnetic field. We have found that the spin dipole mode and single-particle excitations have a sizable overlap, and that the magnetoplasmon modes can be excited by the dipole spin operator if the dot is spin polarized. The frequency of the dipole spin edge mode presents an oscillation which is clearly filling factor (v) related. We have found that the spin dipole mode is especially soft for even-n values. Results for selected numbers of electrons and confining potentials are discussed.

Gradient-dependent density functionals of the Perdew-Burke-Ernzerhof type for atoms, molecules, and solids

One of the standard generalized-gradient approximations (GGAs) in use in modern electronic-structure theory [Perdew-Burke-Ernzerhof (PBE) GGA] and a recently proposed modification designed specifically for solids (PBEsol) are identified as particular members of a family of functionals taking their parameters from different properties of homogeneous or inhomogeneous electron liquids. Three further members of this family are constructed and tested, together with the original PBE and PBEsol, for atoms, molecules, and solids. We find that PBE, in spite of its popularity in solid-state physics and quantum chemistry, is not always the best performing member of the family and that PBEsol, in spite of having been constructed specifically for solids, is not the best for solids. The performance of GGAs for finite systems is found to sensitively depend on the choice of constraints stemming from infinite systems. Guidelines both for users and for developers of density functionals emerge from this work.

Molecular dynamics, density functional, ADMET predictions, virtual screening, and molecular interaction field studies for identification and evaluation of novel potential CDK2 inhibitors in cancer therapy

In this work, we have used molecular dynamics, density functional theory, virtual screening, ADMET predictions, and molecular interaction field studies to design and propose eight novel potential inhibitors of CDK2. The eight molecules proposed showed interesting structural characteristics that are required for inhibiting the CDK2 activity and show potential as drug candidates for the treatment of cancer. The parameters related to the Rule of Five were calculated, and only one of the molecules violated more than one parameter. One of the proposals and one of the drug-like compounds selected by virtual screening indicated to be promising candidates for CDK2-based cancer therapy.

Density matrix functional theory that includes pairing correlations

The extension of density functional theory (DFT) to include pairing correlations without formal violation of the particle-number conservation condition is described. This version of the theory can be considered as a foundation of the application of existing DFT plus pairing approaches to atoms, molecules, ultracooled and magnetically trapped atomic Fermi gases, and atomic nuclei where the number of particles is conserved exactly. The connection with Hartree-Fock-Bogoliubov (HFB) theory is discussed, and the method of quasilocal reduction of the nonlocal theory is also described. This quasilocal reduction allows equations of motion to be obtained which are much simpler for numerical solution than the equations corresponding to the nonlocal case. Our theory is applied to the study of some even Sn isotopes, and the results are compared with those obtained in the standard HFB theory and with the experimental ones.

Density matrix functional theory that includes pairing correlations

The extension of density functional theory (DFT) to include pairing correlations without formal violation of the particle-number conservation condition is described. This version of the theory can be considered as a foundation of the application of existing DFT plus pairing approaches to atoms, molecules, ultracooled and magnetically trapped atomic Fermi gases, and atomic nuclei where the number of particles is conserved exactly. The connection with Hartree-Fock-Bogoliubov (HFB) theory is discussed, and the method of quasilocal reduction of the nonlocal theory is also described. This quasilocal reduction allows equations of motion to be obtained which are much simpler for numerical solution than the equations corresponding to the nonlocal case. Our theory is applied to the study of some even Sn isotopes, and the results are compared with those obtained in the standard HFB theory and with the experimental ones.