On the relation between orbital-localization and self-interaction errors in the density functional theory treatment of organic semiconductors

Description of intermolecular charge transfer with subsystem density-functional theory

Efficient quantum-chemical methods that are able to describe intermolecular charge transfer are crucial for modeling organic semiconductors. However, the correct description of intermolecular charge transfer with density-functional theory (DFT) is hampered by the fractional charge error of approximate exchange-correlation (xc) functionals. Here, we investigate the charge transfer induced by an external electric field in a tetrathiafulvalene-tetracyanoquinodimethane (TTF-TCNQ) complex as a test case. For this seemingly simple model system, a supermolecular DFT treatment fails with most conventional xc functionals. Here, we present an extension of subsystem DFT to subsystems with a fractional number of electrons. We show that within such a framework, it becomes possible to overcome the fractional charge error by enforcing the correct dependence of each subsystem's total energy on the subsystem's fractional charge. Such a subsystem DFT approach allows for a correct description of the intermolecular charge transfer in the TTF-TCNQ model complex. The approach presented here can be generalized to larger molecular aggregates and will thus allow for modeling organic semiconductor materials accurately and efficiently.

Electronic structure of BN-aromatics: Choice of reliable computational tools

The importance of having reliable calculation tools to interpret and predict the electronic properties of BN-aromatics is directly linked to the growing interest for these very promising new systems in the field of materials science, biomedical research, or energy sustainability. Ionization energy (IE) is one of the most important parameters to approach the electronic structure of molecules. It can be theoretically estimated, but in order to evaluate their persistence and propose the most reliable tools for the evaluation of different electronic properties of existent or only imagined BN-containing compounds, we took as reference experimental values of ionization energies provided by ultra-violet photoelectron spectroscopy (UV-PES) in gas phase-the only technique giving access to the energy levels of filled molecular orbitals. Thus, a set of 21 aromatic molecules containing B-N bonds and B-N-B patterns has been merged for a comparison between experimental IEs obtained by UV-PES and various theoretical approaches for their estimation. Time-Dependent Density Functional Theory (TD-DFT) methods using B3LYP and long-range corrected CAM-B3LYP functionals are used, combined with the ΔSCF approach, and compared with electron propagator theory such as outer valence Green's function (OVGF, P3) and symmetry adapted cluster-configuration interaction ab initio methods. Direct Kohn-Sham estimation and "corrected" Kohn-Sham estimation are also given. The deviation between experimental and theoretical values is computed for each molecule, and a statistical study is performed over the average and the root mean square for the whole set and sub-sets of molecules. It is shown that (i) ΔSCF+TDDFT(CAM-B3LYP), OVGF, and P3 are the most efficient way for a good agreement with UV-PES values, (ii) a CAM-B3LYP range-separated hybrid functional is significantly better than B3LYP for the purpose, especially for extended conjugated systems, and (iii) the "corrected" Kohn-Sham result is a fast and simple way to predict IEs.

Accurate description of the electronic structure of organic semiconductors by GW methods

Electronic properties associated with charged excitations, such as the ionization potential (IP), the electron affinity (EA), and the energy level alignment at interfaces, are critical parameters for the performance of organic electronic devices. To computationally design organic semiconductors and functional interfaces with tailored properties for target applications it is necessary to accurately predict these properties from first principles. Many-body perturbation theory is often used for this purpose within the GW approximation, where G is the one particle Green's function and W is the dynamically screened Coulomb interaction. Here, the formalism of GW methods at different levels of self-consistency is briefly introduced and some recent applications to organic semiconductors and interfaces are reviewed.

Investigating the electronic structure of a supported metal nanoparticle: Pd in SiCN

A supporting matrix of SiCN does not significantly change the electronic properties of catalytically active Pd nanoparticles.

CuPc/Au(1 1 0): Determination of the azimuthal alignment by a combination of angle-resolved photoemission and density functional theory

Here we report on a combined experimental and theoretical study on the structural and electronic properties of a monolayer of Copper-Phthalocyanine (CuPc) on the Au(1 1 0) surface. Low-energy electron diffraction reveals a commensurate overlayer unit cell containing one adsorbate species. The azimuthal alignment of the CuPc molecule is revealed by comparing experimental constant binding energy (kxky )-maps using angle-resolved photoelectron spectroscopy with theoretical momentum maps of the free molecule's highest occupied molecular orbital (HOMO). This structural information is confirmed by total energy calculations within the framework of van-der-Waals corrected density functional theory. The electronic structure is further analyzed by computing the molecule-projected density of states, using both a semi-local and a hybrid exchange-correlation functional. In agreement with experiment, the HOMO is located about 1.2 eV below the Fermi-level, while there is no significant charge transfer into the molecule and the CuPc LUMO remains unoccupied on the Au(1 1 0) surface.

Size Effects in the Interface Level Alignment of Dye-Sensitized TiO2 Clusters

The efficiency of dye-sensitized solar cells (DSCs) depends critically on the electronic structure of the interfaces in the active region. We employ recently developed dispersion-inclusive density functional theory (DFT) and GW methods to study the electronic structure of TiO2 clusters sensitized with catechol molecules. We show that the energy level alignment at the dye-TiO2 interface is the result of an intricate interplay of quantum size effects and dynamic screening effects and that it may be manipulated by nanostructuring and functionalizing the TiO2. We demonstrate that the energy difference between the catechol LUMO and the TiO2 LUMO, which is associated with the injection loss in DSCs, may be reduced significantly by reducing the dimensions of nanostructured TiO2 and by functionalizing the TiO2 with wide-gap moieties, which contribute additional screening but do not interact strongly with the frontier orbitals of the TiO2 and the dye. Precise control of the electronic structure may be achieved via "interface engineering" in functional nanostructures.

Efficient Parameter Estimation of Generalizable Coarse-Grained Protein Force Fields Using Contrastive Divergence: A Maximum Likelihood Approach

Maximum Likelihood (ML) optimization schemes are widely used for parameter inference. They maximize the likelihood of some experimentally observed data, with respect to the model parameters iteratively, following the gradient of the logarithm of the likelihood. Here, we employ a ML inference scheme to infer a generalizable, physics-based coarse-grained protein model (which includes Go̅-like biasing terms to stabilize secondary structure elements in room-temperature simulations), using native conformations of a training set of proteins as the observed data. Contrastive divergence, a novel statistical machine learning technique, is used to efficiently approximate the direction of the gradient ascent, which enables the use of a large training set of proteins. Unlike previous work, the generalizability of the protein model allows the folding of peptides and a protein (protein G) which are not part of the training set. We compare the same force field with different van der Waals (vdW) potential forms: a hard cutoff model, and a Lennard-Jones (LJ) potential with vdW parameters inferred or adopted from the CHARMM or AMBER force fields. Simulations of peptides and protein G show that the LJ model with inferred parameters outperforms the hard cutoff potential, which is consistent with previous observations. Simulations using the LJ potential with inferred vdW parameters also outperforms the protein models with adopted vdW parameter values, demonstrating that model parameters generally cannot be used with force fields with different energy functions. The software is available at https://sites.google.com/site/crankite/.

Negative differential resistance in a hybrid silicon-molecular system: resonance between the intrinsic surface-states and the molecular orbital

It has been a long-term desire to fabricate hybrid silicon-molecular devices by taking advantages of organic molecules and the existing silicon-based technology. However, one of the challenging tasks is to design applicable functions on the basis of the intrinsic properties of the molecules, as well as the silicon substrates. Here we demonstrate a silicon-molecular system that produces negative differential resistance (NDR) by making use of the well-defined intrinsic surface-states of the Si (111)-√3 × √3-Ag (R3-Ag/Si) surface and the molecular orbital of cobalt(II)-phthalocyanine (CoPc) molecules. From our experimental results obtained using scanning tunneling microscopy/spectroscopy, we find that NDR robustly appears at the Co(2+) ion centers of the CoPc molecules, independent of the adsorption configuration of the CoPc molecules and irrespective of doping type and doping concentration of the silicon substrates. Joint with first principle calculations, we conclude that NDR is originated from the resonance between the intrinsic surface-state band S(1) of the R3-Ag/Si surface and the localized unoccupied Co(2+)d(z(2)) orbital of the adsorbed CoPc molecules. We expect that such a mechanism can be generally used in other silicon-molecular systems.

Collectively induced quantum-confined Stark effect in monolayers of molecules consisting of polar repeating units

The electronic structure of terpyrimidinethiols is investigated by means of density-functional theory calculations for isolated molecules and monolayers. In the transition from molecule to self-assembled monolayer (SAM), we observe that the band gap is substantially reduced, frontier states increasingly localize on opposite sides of the SAM, and this polarization in several instances is in the direction opposite to the polarization of the overall charge density. This behavior can be analyzed by analogy to inorganic semiconductor quantum-wells, which, as the SAMs studied here, can be regarded as semiperiodic systems. There, similar observations are made under the influence of a, typically external, electric field and are known as the quantum-confined Stark effect. Without any external perturbation, in oligopyrimidine SAMs one encounters an energy gradient that is generated by the dipole moments of the pyrimidine repeat units. It is particularly strong, reaching values of about 1.6 eV/nm, which corresponds to a substantial electric field of 1.6 × 10(7) V/cm. Close-lying σ- and π-states turn out to be a particular complication for a reliable description of the present systems, as their order is influenced not only by the docking groups and bonding to the metal, but also by the chosen computational approach. In the latter context we demonstrate that deliberately picking a hybrid functional allows avoiding pitfalls due to the infamous self-interaction error. Our results show that when aiming to build a monolayer with a specific electronic structure one can not only resort to the traditional technique of modifying the molecular structure of the constituents, but also try to exploit collective electronic effects.