Thomas Körzdörfer

Long-range corrected hybrid functionals for π-conjugated systems: Dependence of the range-separation parameter on conjugation length

Long-range corrected (range-separated hybrid) functionals represent a relatively new class of functionals for generalized Kohn-Sham theory that have proven to be very successful, for instance, when it comes to predicting ionization potentials and energy gaps for a wide range of molecules and solids. The results obtained from long-range corrected density functional theory approaches can be improved dramatically, if the range-separation parameter (ω) is optimized for each system separately. In this work, we have optimized ω for a series of π-conjugated molecular systems of increasing length by forcing the resulting functionals to obey the ionization potential-theorem, i.e., that their highest occupied eigenvalue be equal to the ΔSCF ionization potential. The optimized ω values are observed to vary substantially from their default values for the functionals. For highly conjugated chains such as oligoacenes and polyenes, we find that the characteristic length scale of the range-separation, i.e., 1/ω, grows almost linearly with the number of repeat units, for saturated alkane chains, however, 1/ω quickly saturates after 5-6 repeat units. For oligothiophenes, we find that 1/ω grows linearly for the shorter oligomers but then saturates at around 10 repeat units. Our results point to a close relation between the optimal range-separation parameter and the degree of conjugation in the system.

Organic Electronic Materials: Recent Advances in the DFT Description of the Ground and Excited States Using Tuned Range-Separated Hybrid Functionals

CONSPECTUS: Density functional theory (DFT) and its time-dependent extension (TD-DFT) are powerful tools enabling the theoretical prediction of the ground- and excited-state properties of organic electronic materials with reasonable accuracy at affordable computational costs. Due to their excellent accuracy-to-numerical-costs ratio, semilocal and global hybrid functionals such as B3LYP have become the workhorse for geometry optimizations and the prediction of vibrational spectra in modern theoretical organic chemistry. Despite the overwhelming success of these out-of-the-box functionals for such applications, the computational treatment of electronic and structural properties that are of particular interest in organic electronic materials sometimes reveals severe and qualitative failures of such functionals. Important examples include the overestimation of conjugation, torsional barriers, and electronic coupling as well as the underestimation of bond-length alternations or excited-state energies in low-band-gap polymers. In this Account, we highlight how these failures can be traced back to the delocalization error inherent to semilocal and global hybrid functionals, which leads to the spurious delocalization of electron densities and an overestimation of conjugation. The delocalization error for systems and functionals of interest can be quantified by allowing for fractional occupation of the highest occupied molecular orbital. It can be minimized by using long-range corrected hybrid functionals and a nonempirical tuning procedure for the range-separation parameter. We then review the benefits and drawbacks of using tuned long-range corrected hybrid functionals for the description of the ground and excited states of π-conjugated systems. In particular, we show that this approach provides for robust and efficient means of characterizing the electronic couplings in organic mixed-valence systems, for the calculation of accurate torsional barriers at the polymer limit, and for the reliable prediction of the optical absorption spectrum of low-band-gap polymers. We also explain why the use of standard, out-of-the-box range-separation parameters is not recommended for the DFT and/or TD-DFT description of the ground and excited states of extended, pi-conjugated systems. Finally, we highlight a severe drawback of tuned range-separated hybrid functionals by discussing the example of the calculation of bond-length alternation in polyacetylene, which leads us to point out the challenges for future developments in this field.

On the relationship between bond-length alternation and many-electron self-interaction error

Predicting accurate bond-length alternations (BLAs) in long conjugated molecular chains has been a major challenge for electronic-structure theory for many decades. While Hartree-Fock (HF) overestimates BLA significantly, second-order perturbation theory and commonly used density functional theory (DFT) approaches typically underestimate it. Here, we discuss how this failure is related to the many-electron self-interaction error (MSIE), which is inherent to both HF and DFT approaches. We use tuned long-range corrected hybrids to minimize the MSIE for a series of polyenes. The key result is that the minimization of the MSIE alone does not yield accurate BLAs. On the other hand, if the range-separation parameter is tuned to yield accurate BLAs, we obtain a significant MSIE that grows with chain length. Our findings demonstrate that reducing the MSIE is one but not the only important aspect necessary to obtain accurate BLAs from density functional theory.

Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules II: Non-Empirically Tuned Long-Range Corrected Hybrid Functionals

The performance of non-empirically tuned long-range corrected hybrid functionals for the prediction of vertical ionization potentials (IPs) and electron affinities (EAs) is assessed for a set of 24 organic acceptor molecules. Basis set-extrapolated coupled cluster singles, doubles, and perturbative triples [CCSD(T)] calculations serve as a reference for this study. Compared to standard exchange-correlation functionals, tuned long-range corrected hybrid functionals produce highly reliable results for vertical IPs and EAs, yielding mean absolute errors on par with computationally more demanding GW calculations. In particular, it is demonstrated that long-range corrected hybrid functionals serve as ideal starting points for non-self-consistent GW calculations.

A highly K(+)-selective phenylaza-[18]crown-6-lariat-ether-based fluoroionophore and its application in the sensing of K+ ions with an optical sensor film and in cells.

Herein, we report the synthesis of two phenylaza-[18]crown-6 lariat ethers with a coumarin fluorophore (1 and 2) and we reveal that compound 1 is an excellent probe for K(+) ions under simulated physiological conditions. The presence of a 2-methoxyethoxy lariat group at the ortho position of the anilino moiety is crucial to the substantially increased stability of compounds 1 and 2 over their lariat-free phenylaza-[18]crown-6 ether analogues. Probe 1 shows a high K(+)/Na(+) selectivity and a 2.5-fold fluorescence enhancement was observed in the presence of 100 mM K(+) ions. A fluorescent membrane sensor, which was prepared by incorporating probe 1 into a hydrogel, showed a fully reversible response, a response time of 150 s, and a signal change of 7.8% per 1 mM K(+) within the range 1-10 mM K(+). The membrane was easily fabricated (only a single sensing layer on a solid polyester support), yet no leaching was observed. Moreover, compound 1 rapidly permeated into cells, was cytocompatible, and was suitable for the fluorescent imaging of K(+) ions on both the extracellular and intracellular levels.

Electrical response of molecular systems: the power of self-interaction corrected kohn-sham theory.

The accurate prediction of electronic response properties of extended molecular systems has been a challenge for conventional, explicit density functionals. We demonstrate that a self-interaction correction (SIC) implemented rigorously within Kohn-Sham theory via the optimized effective potential (OEP) yields polarizabilities close to the ones from highly accurate wave-function-based calculations and exceeding the quality of exact-change OEP. The orbital structure obtained with the OEP-SIC functional and approximations to it are discussed.

Ionization Energies, Electron Affinities, and Polarization Energies of Organic Molecular Crystals: Quantitative Estimations from a Polarizable Continuum Model (PCM)–Tuned Range-Separated Density Functional Approach

We propose a new methodology for the first-principles description of the electronic properties relevant for charge transport in organic molecular crystals. This methodology, which is based on the combination of a nonempirical, optimally tuned range-separated hybrid functional with the polarizable continuum model, is applied to a series of eight representative molecular semiconductor crystals. We show that it provides ionization energies, electron affinities, and transport gaps in very good agreement with experimental values, as well as with the results of many-body perturbation theory within the GW approximation at a fraction of the computational costs. Hence, this approach represents an easily applicable and computationally efficient tool to estimate the gas-to-crystal phase shifts of the frontier-orbital quasiparticle energies in organic electronic materials.

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 × 107 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.

Accurate description of torsion potentials in conjugated polymers using density functionals with reduced self-interaction error

We investigate the torsion potentials in two prototypical π-conjugated polymers, polyacetylene and polydiacetylene, as a function of chain length using different flavors of density functional theory. Our study provides a quantitative analysis of the delocalization error in standard semilocal and hybrid density functionals and demonstrates how it can influence structural and thermodynamic properties. The delocalization error is quantified by evaluating the many-electron self-interaction error (MESIE) for fractional electron numbers, which allows us to establish a direct connection between the MESIE and the error in the torsion barriers. The use of non-empirically tuned long-range corrected hybrid functionals results in a very significant reduction of the MESIE and leads to an improved description of torsion barrier heights. In addition, we demonstrate how our analysis allows the determination of the effective conjugation length in polyacetylene and polydiacetylene chains.

Benchmark of GW methods for azabenzenes

Many-body perturbation theory in the GW approximation is a useful method for describing electronic properties associated with charged excitations. A hierarchy of GW methods exists, starting from non-self-consistent G0W0, through partial self-consistency in the eigenvalues (ev-scGW) and in the Green function (scGW0), to fully self-consistent GW (scGW). Here, we assess the performance of these methods for benzene, pyridine, and the diazines. The quasiparticle spectra are compared to photoemission spectroscopy (PES) experiments with respect to all measured particle removal energies and the ordering of the frontier orbitals. We find that the accuracy of the calculated spectra does not match the expectations based on their level of self-consistency. In particular, for certain starting points G0W0 and scGW0 provide spectra in better agreement with the PES than scGW.