Quantum Chemistry

2016
Vlček, V. ; Eisenberg, H. R. ; Steinle-Neumann, G. ; Rabani, E. ; Neuhauser, D. ; Baer, R. Spontaneous charge carrier localization in extended one-dimensional systems. Phys. Rev. Lett. 2016, 116, 186401.Abstract

Charge carrier localization in extended atomic systems has been described previously as being driven by disorder, point defects, or distortions of the ionic lattice. Here we show for the first time by means of first-principles computations that charge carriers can spontaneously localize due to a purely electronic effect in otherwise perfectly ordered structures. Optimally tuned range-separated density functional theory and many-body perturbation calculations within the GW approximation reveal that in trans-polyacetylene and polythiophene the hole density localizes on a length scale of several nanometers. This is due to exchange-induced translational symmetry breaking of the charge density. Ionization potentials, optical absorption peaks, excitonic binding energies, and the optimally tuned range parameter itself all become independent of polymer length as it exceeds the critical localization length. Moreover, we find that lattice disorder and the formation of a polaron result from the charge localization in contrast to the traditional view that lattice distortions precede charge localization. Our results can explain experimental findings that polarons in conjugated polymers form instantaneously after exposure to ultrafast light pulses.

vlcek2016.pdf
2015
Eisenberg, H. R. ; Baer, R. Exothermic Mechanism for the Abstraction of Hydrogen from Methane on Li-Doped MgO. 2015, 119, 196–215. Publisher's VersionAbstract

Despite its importance, it is still not known by which mechanism Li-doped MgO catalyzes the oxidative coupling of methane to ethane. Nevertheless, it is commonly assumed that the mechanism goes through catalytic H abstraction from methane via a Li+O– surface defect. In this paper we use first-principles density functional theory calculations to show that the reaction is significantly more exothermic when the Li+O– defect is situated on a step edge instead of on the flat surface. We find that the reaction on the step is exothermic by 0.25 eV, whereas it is endothermic by 0.3 eV on the flat surface. The presence of the Li dopant in the step edge is crucial for the exothermicity of the reaction. These findings suggest that surface steps which include lithium defects could be responsible for the catalytic behavior of Li/MgO. Following the binding of hydrogen to the Li+O– defect on the step edge the methyl radical can either depart to the gas phase or bind to an adjacent step-edge oxygen atom, increasing the exothermicity of the overall process to 0.8 eV. Activation energies of 0.2 eV for the first pathway and 0.5–0.8 eV for the second were calculated.

Neuhauser, D. ; Rabani, E. ; Cytter, Y. ; Baer, R. Stochastic Optimally Tuned Range-Separated Hybrid Density Functional Theory. J. Phys. Chem. A 2015, 120, 3071–3078.Abstract

We develop a stochastic formulation of the optimally tuned range-separated hybrid density functional theory that enables significant reduction of the computational effort and scaling of the nonlocal exchange operator at the price of introducing a controllable statistical error. Our method is based on stochastic representations of the Coulomb convolution integral and of the generalized Kohn–Sham density matrix. The computational cost of the approach is similar to that of usual Kohn–Sham density functional theory, yet it provides a much more accurate description of the quasiparticle energies for the frontier orbitals. This is illustrated for a series of silicon nanocrystals up to sizes exceeding 3000 electrons. Comparison with the stochastic GW many-body perturbation technique indicates excellent agreement for the fundamental band gap energies, good agreement for the band edge quasiparticle excitations, and very low statistical errors in the total energy for large systems. The present approach has a major advantage over one-shot GW by providing a self-consistent Hamiltonian that is central for additional postprocessing, for example, in the stochastic Bethe–Salpeter approach.

neuhauser2015.pdf
Rabani, E. ; Baer, R. ; Neuhauser, D. Time-dependent stochastic Bethe-Salpeter approach. Phys. Rev. B 2015, 91, 235302.Abstract

A time-dependent formulation for electron-hole excitations in extended finite systems, based on the Bethe-Salpeter equation (BSE), is developed using a stochastic wave function approach. The time-dependent formulation builds on the connection between time-dependent Hartree-Fock (TDHF) theory and the configuration-interaction with single substitution (CIS) method. This results in a time-dependent Schrödinger-like equation for the quasiparticle orbital dynamics based on an effective Hamiltonian containing direct Hartree and screened exchange terms, where screening is described within the random-phase approximation (RPA). To solve for the optical-absorption spectrum, we develop a stochastic formulation in which the quasiparticle orbitals are replaced by stochastic orbitals to evaluate the direct and exchange terms in the Hamiltonian as well as the RPA screening. This leads to an overall quadratic scaling, a significant improvement over the equivalent symplectic eigenvalue representation of the BSE. Application of the time-dependent stochastic BSE (TDsBSE) approach to silicon and CdSe nanocrystals up to size of  3000 electrons is presented and discussed.

rabani2015.pdf
Gao, Y. ; Neuhauser, D. ; Baer, R. ; Rabani, E. Sublinear scaling for time-dependent stochastic density functional theory. J. Chem. Phys. 2015, 142, 034106.Abstract

A stochastic approach to time-dependent density functional theory is developed for computing the absorption cross section and the random phase approximation (RPA) correlation energy. The core idea of the approach involves time-propagation of a small set of stochastic orbitals which are first projected on the occupied space and then propagated in time according to the time-dependent Kohn-Sham equations. The evolving electron density is exactly represented when the number of random orbitals is infinite, but even a small number ( 16) of such orbitals is enough to obtain meaningful results for absorption spectrum and the RPA correlation energy per electron. We implement the approach for silicon nanocrystals using real-space grids and find that the overall scaling of the algorithm is sublinear with computational time and memory.

gao2015.pdf
2014
Neuhauser, D. ; Gao, Y. ; Arntsen, C. ; Karshenas, C. ; Rabani, E. ; Baer, R. Breaking the Theoretical Scaling Limit for Predicting Quasiparticle Energies: The Stochastic GW Approach. Phys. Rev. Lett. 2014, 113, 076402.Abstract

We develop a formalism to calculate the quasiparticle energy within the GW many-body perturbation correction to the density functional theory. The occupied and virtual orbitals of the Kohn-Sham Hamiltonian are replaced by stochastic orbitals used to evaluate the Green function G, the polarization potential W, and, thereby, the GW self-energy. The stochastic GW (sGW) formalism relies on novel theoretical concepts such as stochastic time-dependent Hartree propagation, stochastic matrix compression, and spatial or temporal stochastic decoupling techniques. Beyond the theoretical interest, the formalism enables linear scaling GW calculations breaking the theoretical scaling limit for GW as well as circumventing the need for energy cutoff approximations. We illustrate the method for silicon nanocrystals of varying sizes with Ne > 3000 electrons.

neuhauser2014.pdf
Baratz, A. ; White, A. J. ; Galperin, M. ; Baer, R. Effects of Electromagnetic Coupling on Conductance Switching of a Gated Tunnel Junction. The Journal of Physical Chemistry Letters 2014, 5 3545–3550.Abstract

Using a combination of density functional theory and quantum master equations approach, we study the effect of electromagnetic (EM) coupling on the nonequilibrium steady-state behavior of a recently introduced gated molecular junction. This junction was demonstrated in a previous publication to exhibit sharp current switching near a certain critical DC field Ez*, which induces intramolecular charge transfer, and here, we analyze the steady-state population and current when an AC EM field (EMF) is present. The AC EMF at frequency $ømega_0$ produces pronounced population and current features at gate fields Ez = Ez* ± $\hbar ømega_0/ez$ (where $e_z$ is the dipole of the charge-transfer state) and thus allows additional sharp switching capability at lower gate fields. We found that even when EMF is absent, the EM coupling itself changes the overall steady-state population and current distributions because it allows for relaxation via spontaneous emission

baratz2014.pdf
Egger, D. A. ; Weissman, S. ; Refaely-Abramson, S. ; Sharifzadeh, S. ; Dauth, M. ; Baer, R. ; Kümmel, S. ; Neaton, J. B. ; Zojer, E. ; Kronik, L. Outer-valence Electron Spectra of Prototypical Aromatic Heterocycles from an Optimally Tuned Range-Separated Hybrid Functional. J. Chem. Theory Comput. 2014, 10, 1934–1952.Abstract

Density functional theory with optimally tuned range-separated hybrid (OT-RSH) functionals has been recently suggested [Refaely-Abramson et al. Phys. Rev. Lett. 2012, 109, 226405] as a nonempirical approach to predict the outer-valence electronic structure of molecules with the same accuracy as many-body perturbation theory. Here, we provide a quantitative evaluation of the OT-RSH approach by examining its performance in predicting the outer-valence electron spectra of several prototypical gas-phase molecules, from aromatic rings (benzene, pyridine, and pyrimidine) to more complex organic systems (terpyrimidinethiol and copper phthalocyanine). For a range up to several electronvolts away from the frontier orbital energies, we find that the outer-valence electronic structure obtained from the OT-RSH method agrees very well (typically within  0.1–0.2 eV) with both experimental photoemission and theoretical many-body perturbation theory data in the GW approximation. In particular, we find that with new strategies for an optimal choice of the short-range fraction of Fock exchange, the OT-RSH approach offers a balanced description of localized and delocalized states. We discuss in detail the sole exception found—a high-symmetry orbital, particular to small aromatic rings, which is relatively deep inside the valence state manifold. Overall, the OT-RSH method is an accurate DFT-based method for outer-valence electronic structure prediction for such systems and is of essentially the same level of accuracy as contemporary GW approaches, at a reduced computational cost.

egger2014.pdf
2013
Zohar, G. ; Baer, R. ; Rabani, E. Multiexciton generation in IV–VI nanocrystals: the role of carrier effective mass, band mixing, and phonon emission. The journal of physical chemistry letters 2013, 4 317–322.Abstract

We study the role of the effective mass, band mixing, and phonon emission on multiexciton generation in IV–VI nanocrystals. A four-band k · p effective mass model, which allows for an independent variation of these parameters, is adopted to describe the electronic structure of the nanocrystals. Multiexciton generation efficiencies are calculated using a Green’s function formalism, providing results that are numerically similar to impact excitation. We find that multiexciton generation efficiencies are maximized when the effective mass of the electron and hole are small and similar. Contact with recent experimental results for multiexciton generation in PbS and PbSe is made.

zohar2013multiexciton.pdf
Neuhauser, D. ; Rabani, E. ; Baer, R. Expeditious Stochastic Approach for MP2 Energies in Large Electronic Systems. J. Chem. Theory Comput. 2013, 9 24–27.Abstract

A fast stochastic method for calculating the second order Møller-Plesset (MP2) correction to the correlation energy of large systems of electrons is presented. The approach is based on reducing the exact summation over occupied and unoccupied states to a time-dependent trace formula amenable to stochastic sampling. We demonstrate the abilities of the method to treat systems with thousands of electrons using hydrogen passivated silicon spherical nanocrystals represented on a real space grid, much beyond the capabilities of present day MP2 implementations.

neuhauser2013.pdf
Neuhauser, D. ; Rabani, E. ; Baer, R. Expeditious Stochastic Calculation of Random-Phase Approximation Energies for Thousands of Electrons in Three Dimensions. J. Phys. Chem. Lett. 2013, 4 1172–1176.Abstract

A fast method is developed for calculating the random phase approximation (RPA) correlation energy for density functional theory. The correlation energy is given by a trace over a projected RPA response matrix, and the trace is taken by a stochastic approach using random perturbation vectors. For a fixed statistical error in the total energy per electron, the method scales, at most, quadratically with the system size; however, in practice, due to self-averaging, it requires less statistical sampling as the system grows, and the performance is close to linear scaling. We demonstrate the method by calculating the RPA correlation energy for cadmium selenide and silicon nanocrystals with over 1500 electrons. We find that the RPA correlation energies per electron are largely independent of the nanocrystal size. In addition, we show that a correlated sampling technique enables calculation of the energy difference between two slightly distorted configurations with scaling and a statistical error similar to that of the total energy per electron.

neuhauser2013a.pdf
Ge, Q. ; Gao, Y. ; Baer, R. ; Rabani, E. ; Neuhauser, D. A Guided Stochastic Energy-Domain Formulation of the Second Order Møller–Plesset Perturbation Theory. J. Phys. Chem. Lett. 2013, 5 185–189.Abstract

We develop an alternative formulation in the energy-domain to calculate the second order Møller–Plesset (MP2) perturbation energies. The approach is based on repeatedly choosing four random energies using a nonseparable guiding function, filtering four random orbitals at these energies, and averaging the resulting Coulomb matrix elements to obtain a statistical estimate of the MP2 correlation energy. In contrast to our time-domain formulation, the present approach is useful for both quantum chemistry and real-space/plane wave basis sets. The scaling of the MP2 calculation is roughly linear with system size, providing a useful tool to study dispersion energies in large systems. This is demonstrated on a structure of 64 fullerenes within the SZ basis as well as on silicon nanocrystals using real-space grids.

ge2013.pdf
2012
Baer, R. ; Rabani, E. Expeditious stochastic calculation of multiexciton generation rates in semiconductor nanocrystals. Nano Lett. 2012, 12, 2123–2128.Abstract

A stochastic method is developed to calculate the multiexciton generation (MEG) rates in semiconductor nanocrystals (NCs). The numerical effort scales near-linearly with system size allowing the study of MEG rates up to diameters and exciton energies previously unattainable using atomistic calculations. Illustrations are given for CdSe NCs of sizes and energies relevant to current experimental setups, where direct methods require treatment of over 1011 states. The approach is not limited to the study of MEG and can be applied to calculate other correlated electronic processes.

baer2012expeditious.pdf
2011
Karolewski, A. ; Stein, T. ; Baer, R. ; Kümmel, S. Communication: Tailoring the optical gap in light-harvesting molecules. J. Chem. Phys. 2011, 134, 151101.Abstract

Systematically varying the optical gap that is associated with charge-transfer excitations is an important step in the design of light-harvesting molecules. So far the guidance that time-dependent density functional theory could give in this process was limited by the traditional functionals' inability to describe charge-transfer excitations. We show that a nonempirical range-separated hybrid approach allows to reliably predict charge-transfer excitations for molecules of practically relevant complexity. Calculated absorption energies agree with measured ones. We predict from theory that by varying the number of thiophenes in donor-acceptor-donor molecules, the energy of the lowest optical absorption can be tuned to the lower end of the visible spectrum. Saturation sets in at about five thiophene rings. (C) 2011 American Institute of Physics. [doi:10.1063/1.3581788]

karolewski2011.pdf
2010
Stein, T. ; Eisenberg, H. R. ; Kronik, L. ; Baer, R. Fundamental gaps in finite systems from eigenvalues of a generalized Kohn-Sham method. Phys. Rev. Lett. 2010, 105, 266802.Abstract

We present a broadly applicable, physically motivated, first-principles approach to determining the fundamental gap of finite systems from single-electron orbital energies. The approach is based on using a range-separated hybrid functional within the generalized Kohn-Sham approach to density functional theory. Its key element is the choice of a range-separation parameter such that Koopmans’ theorem for both neutral and anion is obeyed as closely as possible. We demonstrate the validity, accuracy, and advantages of this approach on first, second and third row atoms, the oligoacene family of molecules, and a set of hydrogen-passivated silicon nanocrystals. This extends the quantitative usage of density functional theory to an area long believed to be outside its reach.

stein2010.pdf
Rabani, E. ; Baer, R. Theory of multiexciton generation in semiconductor nanocrystals. Chem. Phys. Lett. 2010, 496, 227–235.Abstract

We develop a generalized framework based on a Green’s function formalism to calculate the efficiency of multiexciton generation in nanocrystal quantum dots. The direct/indirect absorption and coherent/incoherent impact ionization mechanisms, often used to describe multiexciton generation in nanocrystals, are reviewed and rederived from the unified theory as certain approximations. In addition, two new limits are described systematically – the weak Coulomb coupling limit and the semi-wide band limit. We show that the description of multiexciton generation in nanocrystals can be described as incoherent process and we discuss the scaling of multiexciton generation with respect to the photon energy and nanocrystal size. Illustrations are given for three prototype systems: CdSe, InAs and silicon quantum dots.

rabani2010.pdf
Baer, R. ; Rabani, E. Can Impact Excitation Explain Efficient Carrier Multiplication in Carbon Nanotube Photodiodes?. Nano Lett. 2010, 10, 3277–3282.Abstract

We address recent experiments (Science 2009, 325, 1367) reporting on highly efficient multiplication of electron?hole pairs in carbon nanotube photodiodes at photon energies near the carrier multiplication threshold (twice the quasi-particle band gap). This result is surprising in light of recent experimental and theoretical work on multiexciton generation in other confined materials, such as semiconducting nanocrystals. We propose a detailed mechanism based on carrier dynamics and impact excitation resulting in highly efficient multiplication of electron?hole pairs. We discuss the important time and energy scales of the problem and provide analysis of the role of temperature and the length of the diode.

baer2010b.pdf
2009
Livshits, E. ; Baer, R. ; Kosloff, R. Deleterious Effects of Long-Range Self-Repulsion on the Density Functional Description of O-2 Sticking on Aluminum. J. Phys. Chem. A 2009, 113, 7521–7527. Publisher's VersionAbstract

Density functional theory (DFT) with semilocal functionals such as the local-density and generalized gradients approximations predicts that the dissociative adsorption of oxygen on Al (111) goes through without a barrier in stark contradiction to experimental findings. This problem motivated our study of the reaction of oxygen colliding with a small aluminum cluster Al-5. We found semilocal functionals predict a minute barrier to sticking, associated with smeared long-range charge transfer from the metal to the oxygen. Hybrid B3LYP predicts a larger barrier while the range-separated the Baer-Neuhauser-Livshits (BNL, Phys. Chem. Chem. Phys. 2007, 9, 2932.) functional finds a more prominent barrier. BNL predicts short-ranged and more abrupt charge transfer from the surface to the oxygen. We conclude that spurious self-repulsion inherent in semilocal functionals causes early electron-transfer, long-range attraction toward the surface and low reaction barriers for these systems. The results indicate that the missing DFT barrier for O-2 sticking on Al (111) may be due to Spurious self-repulsion.

livshits2009.pdf
Eisenberg, H. R. ; Baer, R. A new generalized Kohn-Sham method for fundamental band-gaps in solids. Phys. Chem. Chem. Phys. 2009, 11, 4674–4680. Publisher's VersionAbstract

We developed a method for calculating the ground-state properties and fundamental band-gaps of solids, using a generalized Kohn-Sham approach combining a local density approximation (LDA) functional with a long-range explicit exchange orbital functional. We found that when the range parameter is selected according to the formula gamma = A/(epsilon(infinity) (epsilon) over tilde) where epsilon(infinity) is the optical dielectric constant of the solid and (epsilon) over tilde = 0.84 and A = 0.216 a(0)(-1), predictions of the fundamental band-gap close to the experimental values are obtained for a variety of solids of different types. For most solids the range parameter g is small (i.e. explicit exchange is needed only at long distances) so the predicted values for lattice constants and bulk moduli are similar to those based on conventional LDA calculations. Preliminary calculations on silicon give a general band structure in good agreement with experiment.

eisenberg2009.pdf
2008
Livshits, E. ; Baer, R. A Density Functional Theory for Symmetric Radical Cations from Bonding to Dissociation. J. Phys. Chem. A 2008, 112, 12789–12791. Publisher's VersionAbstract

It has been known for quite some time that approximate density functional (ADF) theories fail disastrously when describing the dissociative symmetric radical cations R2+. By considering this dissociation limit, previous work has shown that Hartree-Fock (HF) theory favors the R+1-R-0 charge distribution, whereas DF approximations favor the R+(0.5)-R+0.5. Yet, general quantum mechanical principles indicate that both these (as well as all intermediate) average charge distributions are asymptotically energy degenerate. Thus, HF and ADF theories mistakenly break the symmetry but in a contradicting way. In this letter, we show how to construct system-dependent long-range corrected (LC) density functionals that can successfully treat this class of molecules, avoiding the spurious symmetry breaking. Examples and comparisons to experimental data is given for R = H, He, and Ne, and it is shown that the new LC theory improves considerably the theoretical description of the R-2(+) bond properties, the long-range form of the asymptotic potential curve, and the atomic polarizability. The broader impact of this finding is discussed as well, and it is argued that the widespread semiempirical approach which advocates treating the LC parameter as a system-independent parameter is in fact inappropriate under general circumstances.

livshits2008.pdf

Pages