Roi Baer

Reliable Prediction of Charge Transfer Excitations in Molecular Complexes Using Time-Dependent Density Functional Theory

We show how charge transfer excitations at molecular complexes can be calculated quantitatively using time-dependent density functional theory. Predictive power is obtained from range-separated hybrid functionals using nonempirical tuning of the range-splitting parameter. Excellent performance of this approach is obtained for a series of complexes composed of various aromatic donors and the tetracyanoethylene acceptor, paving the way to systematic nonempirical quantitative studies of charge-transfer excitations in real systems.

Tuned Range-Separated Hybrids in Density Functional Theory

We review density functional theory (DFT) within the Kohn-Sham (KS) and the generalized KS (GKS) frameworks from a theoretical perspective for both time-independent and time-dependent problems. We focus on the use of range-separated hybrids within a GKS approach as a practical remedy for dealing with the deleterious long-range self-repulsion plaguing many approximate implementations of DFT. This technique enables DFT to be widely relevant in new realms such as charge transfer, radical cation dimers, and Rydberg excitations. Emphasis is put on a new concept of system-specific range-parameter tuning, which introduces predictive power in applications considered until recently too difficult for DFT.

Excitation Gaps of Finite-Sized Systems from Optimally Tuned Range-Separated Hybrid Functionals

Excitation gaps are of considerable significance in electronic structure theory. Two different gaps are of particular interest. The fundamental gap is defined by charged excitations, as the difference between the first ionization potential and the first electron affinity. The optical gap is defined by a neutral excitation, as the difference between the energies of the lowest dipole-allowed excited state and the ground state. Within many-body perturbation theory, the fundamental gap is the difference between the corresponding lowest quasi-hole and quasi-electron excitation energies, and the optical gap is addressed by including the interaction between a quasi-electron and a quasi-hole. A long-standing challenge has been the attainment of a similar description within density functional theory (DFT), with much debate on whether this is an achievable goal even in principle. Recently, we have constructed and applied a new approach to this problem. Anchored in the rigorous theoretical framework of the generalized Kohn-Sham equation, our method is based on a range-split hybrid functional that uses exact long-range exchange. Its main novel feature is that the range-splitting parameter is not a universal constant but rather is determined from first principles, per system, based on satisfaction of the ionization potential theorem. For finite-sized objects, this DFT approach mimics successfully, to the best of our knowledge for the first time, the quasi-particle picture of many-body theory. Specifically, it allows for the extraction of both the fundamental and the optical gap from one underlying functional, based on the HOMO-LUMO gap of a ground-state DFT calculation and the lowest excitation energy of a linear-response time-dependent DFT calculation, respectively. In particular, it produces the correct optical gap for the difficult case of charge-transfer and charge-transfer-like scenarios, where conventional functionals are known to fail. In this perspective, we overview the formal and practical challenges associated with gap calculations, explain our new approach and how it overcomes previous difficulties, and survey its application to a variety of systems.

Density Functional Theory with Correct Long-Range Asymptotic Behavior

We derive an exact representation of the exchange-correlation energy within density functional theory (DFT) which spawns a class of approximations leading to correct long-range asymptotic behavior. Using a simple approximation, we develop an electronic structure theory that combines a new local correlation energy (based on Monte Carlo calculations applied to the homogeneous electron gas) and a combination of local and explicit long-ranged exchange. The theory is applied to several first-row atoms and diatomic molecules where encouraging results are obtained: good description of the chemical bond at the same time allowing for bound anions, reasonably accurate affinity energies, and correct polarizability of an elongated hydrogen chain. Further stringent tests of DFT are passed, concerning ionization potential and charge distribution under large bias.

PREDICTION OF CHARGE-TRANSFER EXCITATIONS IN COUMARIN-BASED DYES USING A RANGE-SEPARATED FUNCTIONAL TUNED FROM FIRST PRINCIPLES

We study the description of charge-transfer excitations in a series of coumarin-based donor-bridge-acceptor dyes. We show that excellent predictive power for the excitation energies and oscillator strengths in these systems is obtained by using a range-separated hybrid functional within the generalized Kohn-Sham approach to time-dependent density functional theory. Key to this success is a step for tuning the range separation parameter from first principles. We explore different methods for this tuning step, which are variants of a recently suggested approach for charge-transfer excitations [T. Stein et al., J. Am. Chem. Soc. 131, 2818 (2009)]. We assess the quality of prediction by comparing to excitation energies previously published for the same systems using the approximate coupled-cluster singles and doubles (CC2) method.

Quasiparticle Spectra from a Nonempirical Optimally Tuned Range-Separated Hybrid Density Functional

We present a method for obtaining outer-valence quasiparticle excitation energies from a density-functional-theory-based calculation, with an accuracy that is comparable to that of many-body perturbation theory within the GW approximation. The approach uses a range-separated hybrid density functional, with an asymptotically exact and short-range fractional Fock exchange. The functional contains two parameters, the range separation and the short-range Fock fraction. Both are determined nonempirically, per system, on the basis of the satisfaction of exact physical constraints for the ionization potential and frontier-orbital many-electron self-interaction, respectively. The accuracy of the method is demonstrated on four important benchmark organic molecules: perylene, pentacene, 3,4,9,10-perylene-tetracarboxylic-dianydride (PTCDA), and 1,4,5,8-naphthalene-tetracarboxylic-dianhydride (NTCDA). We envision that for the outer-valence excitation spectra of finite systems the approach could provide an inexpensive alternative to GW, opening the door to the study of presently out of reach large-scale systems.

A study of conical intersection effects on scattering processes: The validity of adiabatic single‐surface approximations within a quasi‐Jahn–Teller model

Conical intersections between Born–Oppenheimer potential energy surfaces create singularities which are known to have a direct effect on the symmetry of the nuclear wave functions. In this article is presented a quasi‐Jahn–Teller model to study the symmetry effects of these singularities on nonreactive and reactive scattering processes. Applying this model, we were able to determine in what way and to what extent the conical intersection affects the relevant S‐matrix elements. Having the results of this study available, conclusions concerning more realistic systems were derived.

Density matrix description of laser-induced hot electron mediated photodesorption of NO from Pt(111)

Abstract Based on the numerical solution of the Liouville-von Neumann equation for dissipative systems, the photodesorption dynamics of NO/Pt(111) are studied. Dissipative terms are used to describe the quenching of electronically excited states on the metal, electronic dephasing and the indirect (hot-electron mediated) excitation processes in the DIMET and DIET limits. Norm and energy flow, desorption probabilities and density time-of-flight spectra are computed.

Improved Fermi operator expansion methods for fast electronic structure calculations

Linear scaling algorithms based on Fermi operator expansions (FOE) have been considered significantly slower than other alternative approaches in evaluating the density matrix in Kohn–Sham density functional theory, despite their attractive simplicity. In this work, two new improvements to the FOE method are introduced. First, novel fast summation methods are employed to evaluate a matrix polynomial or Chebyshev matrix polynomial with matrix multiplications totalling roughly twice the square root of the degree of the polynomial. Second, six different representations of the Fermi operators are compared to assess the smallest possible degree of polynomial expansion for a given target precision. The optimal choice appears to be the complementary error function. Together, these advances make the FOE method competitive with the best existing alternatives.

Ab initio study of the alternating current impedance of a molecular junction

The small-bias conductance of the C6 molecule, stretched between two metallic leads, is studied using time-dependent density functional theory within the adiabatic local density approximation. The leads are modeled by jellium slabs, the electronic density and the current density are described on a grid, whereas the core electrons and the highly oscillating valence orbitals are approximated using standard norm-conserving pseudopotentials. The jellium leads are supplemented by a complex absorbing potential that serves to absorb charge reaching the edge of the electrodes and hence mimic irreversible flow into the macroscopic metal. The system is rapidly exposed to a ramp potential directed along the C6 axis, which gives rise to the onset of charge and current oscillations. As time progresses, a fast redistribution of the molecular charge is observed, which translates into a direct current response. Accompanying the dc signal, alternating current fluctuations of charge and currents within the molecule and the metallic leads are observed. These form the complex impedance of the molecule and are especially strong at the plasmon frequency of the leads and the lowest excitation peak of C6. We study the molecular conductance in two limits: the strong coupling limit, where the edge atoms of the chain are submerged in the jellium and the weak coupling case, where the carbon atoms and the leads do not overlap spatially.