Charge transfer, double and bond-breaking excitations with time-dependent density matrix functional theory
|
|
- Laura Gertrude Tucker
- 6 years ago
- Views:
Transcription
1 Version Charge transfer, double and bond-breaking excitations with time-dependent density matrix functional theory K.J.H. Giesbertz, E. J. Baerends, and... Afdeling Theoretische Chemie, Vrije Universiteit, De Boelelaan 183, 181 HV Amsterdam, The Netherlands (Dated: April 3, 28) Time-dependent density functional theory in its current adiabatic implementations exhibits three remarkable failures: a) totally wrong behavior of the excited state surface along a bond-breaking coordinate; b) lack of doubly excited configurations, affecting again excited state surfaces; c) much too low charge transfer excitation energies. These TDDFT failure cases are all strikingly exhibited by prototype two-electron systems such as dissociating H 2 (zero σ g σ u excitation energy at R, missing lowest excited state (a double excitation) at R > 2.5 bohr) and HeH + (charge transfer excitations). The exact density matrix functional being known in two-electron systems, the performance of time-dependent density matrix theory can be tested in those systems. We find for TDDMFT that: a) The first order (linear response) formulation does not prohibit the correct representation of doubly excited states; b) Within previously formulated simple adiabatic approximations the bonding-to-antibonding excited state surface as well as charge transfer excitations are described without problems, but not the double excitations; c) An adiabatic approximation is formulated in which also the double excitations are fully accounted for. I. INTRODUCTION Dissociation of molecular systems poses a serious challenge within the DFT framework. This is already true for the ground state energy curve, cf. e.g. Refs [1, 2]. For excited states the problems are even worse, there are various kinds of excitations which are not correctly represented in TDDFT calculations. Fig. I displays the total failure of TDDFT (in its adiabatic GGA (BP86) variant) for the potential energy surface (PES) of the first excited state of H 2, 1 1 Σ + u. As highlighted in Ref. [3], this PES goes to zero instead of going through a minimum and asymptotically to ca. 1 ev. B3LYP doesn t improve the situation the least. This state corresponds to the simple σ g σ u orbital excitation in H 2. The prototypical photochemical bond breaking event resulting from bonding to antibonding orbital excitation takes place on the triplet surface (1 3 Σ + u ) corresponding to this singlet excited state, hence the great importance of this type of excitation. In Fig. I we show the TDDFT Σ + g excitation energies of dissociating H 2. It is evident that there is an exact state (the third at R e = 1.4 Bohr) which is totally missing in the TDDFT calculations. In the accurate calculations it becomes lower with increasing R and has avoided crossings with the second and first state, becoming the lowest excited 1 Σ + g state (2 1 Σ + g ) from ca. 2.5 Bohr onwards. This is a doubly excited state, (σ g ) 2 (σ u ) 2. It is totally missing in the TDDFT calculations. After some initial optimism that such doubly excited states states would be accurately calculated in TDDFT [4 6], it has become clear this is not the case [7 1]. We note in passing that also in the higher part of the spectrum - above.6 a.u. - TDDFT performs very poorly, there is hardly any agreement with the exact spectrum. In Fig. I the excitation energies of HeH + are shown along the dissociation coordinate for the lowest three 1 Σ excited states. The TDDFT excitation energies exhibit the well-known severe underestimation at long distance, where these excitations have strong charge transfer character, from He(1s) to H(1s) in 2Σ and from He(1s) to H(2s, 2p z ) in 3Σ and 4Σ. The hybrid functional B3LYP improves only little on the pure GGA BP86. We note that the excited state energy does not have a 1/R asymptotic behavior, the fragments being He + and neutral H. A remedy for these deficiencies of time-dependent density functional theory might be sought in density matrix functional theory (DMFT) [11 18]. For the ground state total energy curve along the dissociation coordinate, successes have already been reported [16] with approximate DM functionals that represent improvements on the forms formulated long ago by Müller [11] and Buijse and Baerends [12, 13]. It is the purpose of the present paper to address the abovementioned problem cases for TDDFT excitation energy calculation with DMFT, using a further development of linear response based timedependent density matrix functional theory [19, 2]. We will restrict ourselves to the two-electron systems H 2 and HeH +, which exhibit all the problems, but which have the advantage that we do not have to introduce an approximate functional since the exact one is known for two-electron systems. The spatial part of a singlet twoelectron wave function is symmetric at all times t, so the matrix C(t) of expansion coeeficients in a two-electron product basis χ p (r)χ q (r ) is symmetric and can be diagonalised at each t, as originally noted by Löwdin and Shull [21] Ψ(r, r ; t) = pq = µ C pq (t)χ p (r)χ q (r ) c µ (t)φ µ (rt)φ µ (r t). (1) The orbitals {φ µ } are time-dependent, since the wave function will have a different spectral representation
2 Σ + u 3Σ + u.5 3Σ + 4Σ + 1Σ + u 2Σ + R He-H.8 4Σ + g FIG. 2: Σ + TDDFT excitation energies for HeH +. Solid lines: exact in cc-pvtz basis; dashed: TDDFT-BP86; dotted: TDDFT-B3LYP Σ + g 3Σ + g FIG. 1: Σ + u and Σ + g TDDFT excitation energies for H 2. Solid lines: exact in cc-pvtz basis; dashed: TDDFT-BP86; dotted: TDDFT-B3LYP. matrix reduces to the particularly simple form Γ(r 1, r 2, r 2, r 1; t) = pqrs Γ pqrs (t)χ s(r 1)χ r(r 2)χ q (r 2 )χ p (r 1 ) = µνρσ[ 2cµ (t)c σ(t)δ µν δ στ ] φ σ(r 1t)φ ρ(r 2t)φ ν (r 2 t)φ µ (r 1 t). (3) This shows that in a singlet two-electron system the twomatrix can be reconstructed easily from the one-matrix. It is this property which gives two-electron systems their special place in DMFT. The electron-electron interaction energy W (t) = Ψ(t) Ŵ Ψ(t) can be written in terms of the matrix W W µν (t) ρστ Γ µρστ (t)w τσρν (t) at different times. To distinguish this special timedependent basis {φ µ (t)} from the stationary one, {χ p }, we use Greek letters as indices for the time-dependent basis and Latin for the time-independent functions. The one-particle density matrix can be written = 2 c σ(t)c µ (t)w σσµν (t) (4) σ w pqrs = dr 1 dr 2 χ p(r 1 )χ q(r 2 )w(r 1, r 2 ) χ r (r 2 )χ s (r 1 ) (5) γ(r, r ; t) = pq = µ γ pq (t)χ p (r)χ q(r ) n µ (t)φ µ (rt)φ µ(r t) (2) as W (t) = 1 2Tr{W }. Note that the two-electron integrals are related to the ones in conventional notation as w pqrs = pq sr. The equation of motion of the onematrix can now be written in the form that was derived in [19] and used in [2] which shows the {φ µ (t)} are the time-dependent NOs and the wavefunction coefficients c µ (t) of Eq. 1 are related to the NO occupation numbers, n µ (t) = 2 c µ (t) 2. The two- i γ kl (t) = ( hks (t)γ sl (t) γ ks (t)h sl (t) ) + s ( W kl (t) W kl(t) ), (6)
3 3 where W kl (t) is defined as the contraction of the twomatrix in the stationary basis with the two-electron integrals in the same basis, W kl (t) = rst Γ krst(t)w tsrl (t). The simplifications afforded by the simple representation of the one-matrix and two-matrix in terms of the diagonal wavefunction coefficients {c µ (t)}, can now be exploited to derive (with considerable algebra, see Ref. ) exact equations for the excitation energies. We first define the matrices A kl,sr (δ kr δ ls + δ ks δ lr )E (w klsr + w klrs ) (δ kr h ls + δ ks h lr + δ lr h ks + δ ls h kr ) (7a) C kl,r A kl,rr (7b) E k,r A kk,rr, (7c) where E is the energy of the ground state. Note that A is hermitian, and C and E are just submatrices of A. In order to obtain a symmetric matrix problem, we define quantities related to the real and imaginary parts of the lower triangle density matrix elements δγ kl, l < k and the real diagonal elements δγ kk = δn κ : δ γ R/I kl () δγr/i kl () 2(c l ± c k ) (8) We then obtain the equations CC T A CE A C δ γr () δ γ I () = δ γr () δ γ I () (9) EC T C T E2 δñ() δñ(). The values for which this system gives non-trivial solutions represent excitation energies. The equations have to be solved iteratively, e.g. by diagonalizing the matrix for trial values of, and search for those values of where an eigenvalue ɛ i () becomes equal to the input. Beyond the linear response approach, no approximations have yet been introduced, and indeed the exact excitation energies are found, equal to those form a full CI calculation in the same basis set. This demonstrates that the linear response formulation does not preclude obtaining double excitations. Het berekenen van excitatie enrgieën met lineaire response is geen benadering, dus ik zou gewoon zeggen dat er geen benaderingen zijn gemaakt. For practical applications (systems with more than two electrons), we need adiabatic approximations, which yield an -independent matrix so that simple diagonalisation will yield the excitation energies. The limit of the exact TDDMFT matrix of Eq. 9 is relevant for the static limit of frequency dependent polarizability calculations. It is, however, immediately clear that the 1/ terms in the exact Eq. 9 will diverge for small. Not surprisingly, the so-called static approximation (SA) introduced in Ref. [19] yields a matrix in which those terms have been eliminated: A A C δ γr () δ γ I () = δ γr () δ γ I () (1) C T δñ() δñ() Σ + u Σ + g FIG. 3: Σ + u and Σ + g excitations for H 2. Solid lines: exact and AA2. Dotted: SA. Dashed: AA1. Unfortunately, the SA approximation turns out to be deficient in the sense that it yields m excitation energies equal zero. As we will see, it is precisely the doubly excited states that turn out to be totally wrong (there are m diagonal excitations (ψ 1 ) 2 (ψ i ) 2 ); ψ 1 is the occupied 1σ g or 1σ orbital in our systems). In an attempt to improve upon the SA, the so-called AA has been proposed [2], in which the third set of equations of the SA has been replaced by the static perturbation equation for the energy under a (static) potential perturbation, C T δ γ R () + Eδñ() δṽ D () = c δe(), (11) Again, this approximation, which we will call AA1, does not solve the problem with the double excitations, which are actually totally absent in this approximation. We propose an adiabatic approximation in which the static perturbation equation 11 is incorporated in the third set of the equations for the (δ γ R, δ γ I, δñ) vectors, in such
4 Σ + R He-H FIG. 4: Σ + excitations for HeH +. Solid lines: exact and AA2. Dotted: SA. Dashed: AA1. a way that in the limit these equations reduce to the static perturbation equation. In that way the full set of equations becomes consistent, because imposing the equation 11 is actually a requirement for the elimination of the 1/ terms in the first two sets of equations. We obtain the AA2 set of equations A A C δ γr () δ γ I () = δ γr () δ γ I () (12) C T C T E δñ() δñ(). We compare in Fig. I the adiabatic TDDMFT excitation energies with the exact ones for the Σ + u excited states. Interestingly, in this symmetry the SA, AA1 and AA2 all coincide with the exact excitation energies along the whole binding distance. This is related to the symmetry: we find that for different symmetries than the ground state, all adiabatic approximations yield perfect excitation energies. This also holds for the Π u and Π g excited states. The significant observation, however, is that the adiabatic TDDMFT, in contrast to adiabatic TDDFT, has no problem with the σ g σ u excitation along the whole bond-breaking coordinate. Results for the Σ + g excitations are shown in Fig. I. Now the adiabatic approximations are not identical to the exact ones. However, the performance of the adiabatic approximations SA and AA1 is rather good, except for the fact that they are completely missing the state with predominantly doubly excited character, that is high-lying at short distance and then on its way down crosses the other states (see the crosses indicating the doubly excited electronic character on the different curves). This deficiency of the SA and AA1 adiabatic approximations in TDDMFT is actually very similar to the adiabatic TDDFT deficiency discussed in the Introduction. On the other hand, the AA2 approximation coincides with the exact results. It brings in the coupling to the diagonal doubles very effectively, in fact perfectly. Finally, the HeH + excitations in the Σ + symmetry are shown in Fig I. The low lying Σ + excitations are, similar to the Σ + g excitations in H 2, not perfect in the SA and AA1 approximations, but they are clearly rather accurate. In this case the AA1 performs better than the SA, especially for the lowest Σ + excitation and also for the other two. The AA2 is exact for all the Σ (and Π) excited states. The important point is that at long distance the lowest excited state, 2Σ + at ca..4 a.u., is a pure CT state, He(1s) H(1s). On the other hand, 3Σ + at ca..74 a.u. is at long distance a mixed local (He(1s) He(2s)) and charge transfer (He(1s) H(2s, 2p)) excitation. The 4Σ + at ca..82 a.u. is mostly CT, He(1s) H(2s, 2p). The different characters of these states do not seem to influence the accuracy of the adiabatic TD- DMFT approximations. We conclude that the well-known failure cases of TDDFT can in principle be solved without difficulty in adiabatic TDDMFT. Whether adiabatic TDDMFT will prove to be widely applicable, also in N-electron systems with N > 2, will depend on the quality of the approximate density matrix functionals that are needed in those cases. Good progress is being made with such functionals for the ground state, both at equilibrium geometry [22] and along the dissociation coordinate,[16] which offers good prospects for the response calculations. [1] E. J. Baerends, Phys. Rev. Lett. 87, 1334 (21). [2] M. Grüning, O. V. Gritsenko and E.J. Baerends, J. Chem. Phys. 118, 7183 (23). [3] O. Gritsenko, S. J. A. van Gisbergen, A. Görling and E.J. Baerends, J. Chem. Phys. 113, 8478 (2). [4] S. Hirata and M. Head-Gordon, Chem. Phys. Lett. 314, 291 (1999). [5] C-P Hsu, S. Hirata and M. Head-Gordon, J. Phys. Chem. A 15, 451 (21). [6] D.J. Tozer, R.D. Amos, N.C. Handy, B.J. Roos and L. Serrano-Andres, Mol. Phys. 97, 859 (1999). [7] N.T. Maitra, F. Zhang, R.J. Cave and K. Burke, J. Chem. Phys. 12, 5932 (24). [8] R.J. Cave, F. Zhang, N.T. Maitra and K. Burke, Chem. Phys. Lett. 389, 39 (24). [9] J. Neugebauer and E.J. Baerends, J. Chem. Phys. 121, 6155 (24). [1] I.A. Mikhailov, Phys. Rev. A 77, 1251 (28). [11] A.M.K. Müller, Phys. Lett. A 19, 446 (1984). [12] M. Buijse, Ph.D. thesis, Vrije Universiteit, Amsterdam (1991). [13] M. Buijse and E.J. Baerends, Mol. Phys. 1, 41 (22).
5 5 [14] C. Kollmar and B.A. Heß, J. Chem. Phys. 12, 3158 (24). [15] C. Kollmar, J. Chem. Phys. 121, (24). [16] O.V. Gritsenko, K. Pernal and E.J. Baerends, J. Chem. Phys. 122, 2412 (25). [17] M. Piris, Int. J. Quant. Chem. 16, 193 (26). [18] M. Piris, X. Lopez and J.M. Ugalde, J. Chem. Phys. 126, (27). [19] K. Pernal, O. Gritsenko and E.J. Baerends, Phys. Rev. A 75, 1256 (27). [2] K. Pernal, K. Giesbertz, O. Gritsenko and E.J. Baerends, J. Chem. Phys. 127, (27). [21] P.-O. Löwdin and H. Shull, Phys. Rev. 11, 173 (1956). [22] P. Leiva and M. Piris, J. Chem. Phys. 123, (25).
Excitation energies of dissociating H 2 : A problematic case for the adiabatic approximation of time-dependent density functional theory
JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 19 15 NOVEMBER 2000 Eitation energies of dissociating H 2 : A problematic case for the adiabatic approximation of time-dependent density functional theory
More informationTDDFT Explorations in Phase-Space. Neepa T. Maitra Hunter College and the Graduate Center of the City University of New York
TDDFT Explorations in Phase-Space Neepa T. Maitra Hunter College and the Graduate Center of the City University of New York TDDFT in Phase-Space Outline Why phase-space DFT? Some motivating examples and
More informationIntroduction and Overview of the Reduced Density Matrix Functional Theory
Introduction and Overview of the Reduced Density Matrix Functional Theory N. N. Lathiotakis Theoretical and Physical Chemistry Institute, National Hellenic Research Foundation, Athens April 13, 2016 Outline
More informationExtending Kohn-Sham density-functional theory
Extending Kohn-Sham density-functional theory Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Email: julien.toulouse@upmc.fr Web page: www.lct.jussieu.fr/pagesperso/toulouse/ January
More informationCoupled-Cluster Perturbative Triples for Bond Breaking
Coupled-Cluster Perturbative Triples for Bond Breaking Andrew G. Taube and Rodney J. Bartlett Quantum Theory Project University of Florida INT CC Meeting Seattle July 8, 2008 Why does chemistry need triples?
More informationIntroduction to multiconfigurational quantum chemistry. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Notations
More informationLinear response time-dependent density functional theory
Linear response time-dependent density functional theory Emmanuel Fromager Laboratoire de Chimie Quantique, Université de Strasbourg, France fromagere@unistra.fr Emmanuel Fromager (UdS) Cours RFCT, Strasbourg,
More informationDensity-functional-theory response-property calculations with accurate exchange-correlation potentials
PHYSICAL REVIEW A VOLUME 57, NUMBER 4 APRIL 1998 Density-functional-theory response-property calculations with accurate exchange-correlation potentials S. J. A. van Gisbergen, 1 F. Kootstra, 2 P. R. T.
More informationAdiabatic connection for near degenerate excited states
PHYSICAL REVIEW A 69, 052510 (2004) Adiabatic connection for near degenerate excited states Fan Zhang Department of Physics and Astronomy, Rutgers University, 136 Frelinghuysen Road, Piscataway, New Jersey
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationTDDFT in Chemistry and Biochemistry III
TDDFT in Chemistry and Biochemistry III Dmitrij Rappoport Department of Chemistry and Chemical Biology Harvard University TDDFT Winter School Benasque, January 2010 Dmitrij Rappoport (Harvard U.) TDDFT
More informationTime-Dependent Density-Functional Theory
Summer School on First Principles Calculations for Condensed Matter and Nanoscience August 21 September 3, 2005 Santa Barbara, California Time-Dependent Density-Functional Theory X. Gonze, Université Catholique
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationSupporting Information
Supporting Information The Journal of Physical Chemistry A Determination of Binding Strengths of a Host-Guest Complex Using Resonance Raman Scattering Edward H. Witlicki, Stinne W. Hansen, Martin Christensen,
More informationDensity matrix functional theory vis-á-vis density functional theory
Density matrix functional theory vis-á-vis density functional theory 16.4.007 Ryan Requist Oleg Pankratov 1 Introduction Recently, there has been renewed interest in density matrix functional theory (DMFT)
More informationAccurate atomic quantum defects from particle-particle random phase approximation
Accurate atomic quantum defects from particle-particle random phase approximation Yang Yang Department of Chemistry, Duke University, Durham, NC 27708, U.S. Kieron Burke Department of Chemistry, University
More informationOrbital dependent correlation potentials in ab initio density functional theory
Orbital dependent correlation potentials in ab initio density functional theory noniterative - one step - calculations Ireneusz Grabowski Institute of Physics Nicolaus Copernicus University Toruń, Poland
More informationDiscontinuity of the chemical potential in reduced-density-matrix-functional theory for open-shell systems
PHYSICAL REVIEW A 79, 022504 2009 Discontinuity of the chemical potential in reduced-density-matrix-functional theory for open-shell systems N. Helbig, 1,2,3 N. N. Lathiotakis, 1,3,4 and E. K. U. Gross
More informationf xc Advanced TDDFT PDF Created with deskpdf PDF Writer - Trial :: II. Frequency-Dependent Kernels: Double Excitations
Advanced TDDFT II. Frequency-Dependent Kernels: Double Excitations f xc Neepa T. Maitra Hunter College and the Graduate Center of the City University of New York First, quick recall of how we get excitations
More informationAdvanced TDDFT II. Double Excitations in TDDFT
Advanced TDDFT II. Double Excitations in TDDFT f xc Neepa T. Maitra Hunter College and the Graduate Center of the City University of New York First, quick recall of how we get excitations in TDDFT: Linear
More informationQuantum Mechanics Solutions
Quantum Mechanics Solutions (a (i f A and B are Hermitian, since (AB = B A = BA, operator AB is Hermitian if and only if A and B commute So, we know that [A,B] = 0, which means that the Hilbert space H
More informationDensity functional calculation of nuclear magnetic resonance chemical shifts
Density functional calculation of nuclear magnetic resonance chemical shifts Christoph van Wüllen Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany Received 22 August
More informationAsymptotic nodal planes in the electron density and the potential in the effective equation for the square root of the density
Eur. Phys. J. B 018) 91: 160 https://doi.org/10.1140/epjb/e018-905-3 Regular Article THE EUROPEAN PHYSICAL JOURNAL B Asymptotic nodal planes in the electron density and the potential in the effective equation
More informationDensity matrix functional theory. Rob Klooster
Density matrix functional theory Rob Klooster March 26, 2007 Summary In this report, we will give an overview of Density Matrix Functional Theory (DMFT). In the first part we will discuss the extended
More informationQUANTUM CHEMISTRY FOR TRANSITION METALS
QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline
More informationAssessment of range-separated time-dependent density-functional theory for calculating C 6 dispersion coefficients
1/10 Assessment of range-separated time-dependent density-functional theory for calculating C 6 dispersion coefficients Julien Toulouse 1,2, Elisa Rebolini 1, Tim Gould 3, John F. Dobson 3, Prasenjit Seal
More informationJULIEN TOULOUSE, PAOLA GORI-GIORGI, ANDREAS SAVIN
Scaling Relations, Virial Theorem, and Energy Densities for Long-Range and Short-Range Density Functionals JULIEN TOULOUSE, PAOLA GORI-GIORGI, ANDREAS SAVIN Laboratoire de Chimie Théorique, CNRS et Université
More informationElectric properties of molecules
Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the
More informationand states of CH 2 Non-adiabatic effects in the à 2 B 2 through coupling vibrational modes Journal of Physics: Conference Series Related content
Journal of Physics: Conference Series PAPER OPEN ACCESS Non-adiabatic effects in the à B and states of CH F + through coupling vibrational modes To cite this article: Rudraditya Sarkar 016 J. Phys.: Conf.
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry Chemicum 4th floor vesa.hanninen@helsinki.fi September 10, 2013 Lecture 3. Electron correlation methods September
More informationElectronic Excitation Energies of Molecular Systems from the Bethe Salpeter Equation
8 Electronic Excitation Energies of Molecular Systems from the Bethe Salpeter Equation Example of the H Molecule Elisa Rebolini, Julien Toulouse, and Andreas Savin CONTENTS 8. Introduction... 67 8. Review
More informationarxiv: v2 [physics.chem-ph] 14 Jan 2013
Multi-configuration time-dependent density-functional theory based on range separation Emmanuel Fromager a, Stefan Knecht b and Hans Jørgen Aa. Jensen b a Laboratoire de Chimie Quantique, arxiv:1211.4829v2
More informationMulticonfigurational Quantum Chemistry. Björn O. Roos as told by RL Department of Theoretical Chemistry Chemical Center Lund University Sweden
Multiconfigurational Quantum Chemistry Björn O. Roos as told by RL Department of Theoretical Chemistry Chemical Center Lund University Sweden April 20, 2009 1 The Slater determinant Using the spin-orbitals,
More informationStudy of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal
Study of Ozone in Tribhuvan University, Kathmandu, Nepal Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal 1 Country of the Mt Everest 2 View of the Mt Everest 3 4 5
More informationOn the Uniqueness of Molecular Orbitals and limitations of the MO-model.
On the Uniqueness of Molecular Orbitals and limitations of the MO-model. The purpose of these notes is to make clear that molecular orbitals are a particular way to represent many-electron wave functions.
More informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
More informationHow to make the Born-Oppenheimer approximation exact: A fresh look at potential energy surfaces and Berry phases
How to make the Born-Oppenheimer approximation exact: A fresh look at potential energy surfaces and Berry phases E.K.U. Gross Max-Planck Institute of Microstructure Physics Halle (Saale) Process of vision
More informationCorrelated two-electron quantum dynamics in intense laser fields
Correlated two-electron quantum dynamics in intense laser fields Dieter Bauer University of Rostock, Germany 7th Workshop on Time-Dependent Density Functional Theory Benasque, Spain, September 20 23, 2016
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationTime-dependent density functional theory
Time-dependent density functional theory E.K.U. Gross Max-Planck Institute for Microstructure Physics OUTLINE LECTURE I Phenomena to be described by TDDFT Some generalities on functional theories LECTURE
More informationDensity Functional Theory (DFT) Time Dependent Density Functional Theory The working equations Performance Properties.
DFT and TDDFT Contents 1 Density Functional Theory (DFT) 2 3 4 4 Basic literature: Jensen: Introduction to Computational Chemistry, Wiley and Sons. Koch, Holthausen: A Chemist s Guide to Density Functional
More informationSpeed of light c = m/s. x n e a x d x = 1. 2 n+1 a n π a. He Li Ne Na Ar K Ni 58.
Physical Chemistry II Test Name: KEY CHEM 464 Spring 18 Chapters 7-11 Average = 1. / 16 6 questions worth a total of 16 points Planck's constant h = 6.63 1-34 J s Speed of light c = 3. 1 8 m/s ħ = h π
More informationMulti-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
More informationElectron States of Diatomic Molecules
IISER Pune March 2018 Hamiltonian for a Diatomic Molecule The hamiltonian for a diatomic molecule can be considered to be made up of three terms Ĥ = ˆT N + ˆT el + ˆV where ˆT N is the kinetic energy operator
More informationarxiv: v2 [cond-mat.other] 27 Jun 2008
Self-interaction correction in a simple model arxiv:0804.0684v [cond-mat.other] 7 Jun 008 P. M. Dinh a, J. Messud a, P.-G. Reinhard b, and E. Suraud a Abstract a Laboratoire de Physique Théorique, Université
More informationTDDFT as a tool in biophysics
TDDFT as a tool in biophysics The primary event in vision Robert Send Universität Karlsruhe 09.09.08 Robert Send TDDFT as a tool in biophysics 09.09.08 1 / 28 Outline 1 Human vision 2 The methods 3 The
More informationLecture 3: Quantum Satis*
Lecture 3: Quantum Satis* Last remarks about many-electron quantum mechanics. Everything re-quantized! * As much as needed, enough. Electron correlation Pauli principle Fermi correlation Correlation energy
More informationTheoretical Photochemistry WiSe 2017/18
Theoretical Photochemistry WiSe 2017/18 Lecture 7 Irene Burghardt (burghardt@chemie.uni-frankfurt.de) http://www.theochem.uni-frankfurt.de/teaching/ Theoretical Photochemistry 1 Topics 1. Photophysical
More informationSelf-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT
Self-Consistent Implementation of Self-Interaction Corrected DFT and of the Exact Exchange Functionals in Plane-Wave DFT Kiril Tsemekhman (a), Eric Bylaska (b), Hannes Jonsson (a,c) (a) Department of Chemistry,
More informationBuilding a wavefunction within the Complete-Active. Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy
Building a wavefunction within the Complete-Active Active-Space Coupled-Cluster Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy Dmitry I. Lyakh (Karazin Kharkiv
More informationIntermediate DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA
Intermediate DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA October 10-19th, 2012 Kieron (UC Irvine) Intermediate DFT Lausanne12
More informationXYZ of ground-state DFT
XYZ of ground-state DFT Kieron Burke and Lucas Wagner Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA January 5-9th, 2014 Kieron (UC Irvine) XYZ of ground-state
More informationFrom the Bethe-Salpeter equation to non-relativistic approaches with effective two-body interactions
INT, Seattle - 15.11.2004 p. 1/28 From the Bethe-Salpeter equation to non-relativistic approaches with effective two-body interactions Lukas THEUSSL theussl@triumf.ca Triumf, 4004 Wesbrook Mall, Vancouver,
More informationAdvanced TDDFT II. II. Frequency-Dependent Kernels: Double Excitations. f xc
Advanced TDDFT II II. Frequency-Dependent Kernels: Double Excitations f xc Neepa T. Maitra Hunter College and the Graduate Center of the City University of New York Plan -- Double-Excitations in TDDFT
More informationWalter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.
Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm
More informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for
More informationShort Course on Density Functional Theory and Applications VIII. Time-dependent DFT
Short Course on Density Functional Theory and Applications VIII. Time-dependent DFT Samuel B. Trickey Sept. 2008 Quantum Theory Project Dept. of Physics and Dept. of Chemistry trickey@qtp.ufl.edu Time-dependent
More informationarxiv: v1 [cond-mat.str-el] 10 Aug 2010
ROHF Theory Made Simple Takashi Tsuchimochi and Gustavo E. Scuseria Department of Chemistry, Rice University, Houston, TX 775-1892 Department of Physics and Astronomy, Rice University, Houston, TX 775-1892
More informationVALENCE Hilary Term 2018
VALENCE Hilary Term 2018 8 Lectures Prof M. Brouard Valence is the theory of the chemical bond Outline plan 1. The Born-Oppenheimer approximation 2. Bonding in H + 2 the LCAO approximation 3. Many electron
More informationCombining density-functional theory and many-body methods
Combining density-functional theory and many-body methods Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Vrije Universiteit Amsterdam, Netherlands November 2017 Outline 2/23 1
More informationPopulation Analysis. Mulliken Population Analysis APPENDIX S
APPENDIX S Population Analysis On p. 665, electronic density ρ is defined. If the wave function is a Slater determinant p. 397) and assuming the double occupancy of orbitals ϕ i, we have see 11.7) ρ r
More informationSCF calculation on HeH +
SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the
More informationOptimization of quantum Monte Carlo wave functions by energy minimization
Optimization of quantum Monte Carlo wave functions by energy minimization Julien Toulouse, Roland Assaraf, Cyrus J. Umrigar Laboratoire de Chimie Théorique, Université Pierre et Marie Curie and CNRS, Paris,
More informationPractical Issues on the Use of the CASPT2/CASSCF Method in Modeling Photochemistry: the Selection and Protection of an Active Space
Practical Issues on the Use of the CASPT2/CASSCF Method in Modeling Photochemistry: the Selection and Protection of an Active Space Roland Lindh Dept. of Chemistry Ångström The Theoretical Chemistry Programme
More information3: Many electrons. Orbital symmetries. l =2 1. m l
3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience
More informationModule 6 1. Density functional theory
Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft
More informationNWChem: Coupled Cluster Method (Tensor Contraction Engine)
NWChem: Coupled Cluster Method (ensor Contraction Engine) What we want to solve H Ψ = E Ψ Many Particle Systems Molecular/Atomic Physics, Quantum Chemistry (electronic Schrödinger equations) Solid State
More informationBeyond the Hartree-Fock Approximation: Configuration Interaction
Beyond the Hartree-Fock Approximation: Configuration Interaction The Hartree-Fock (HF) method uses a single determinant (single electronic configuration) description of the electronic wavefunction. For
More informationMulti-configuration time-dependent density-functional theory based on range separation
Multi-configuration time-dependent density-functional theory based on range separation Emmanuel Fromager, Stefan Knecht, and Hans Jørgen Aa. Jensen Citation: J. Chem. Phys. 138, 8411 213); doi: 1.163/1.4792199
More informationA new generation of density-functional methods based on the adiabatic-connection fluctuation-dissipation theorem
A new generation of density-functional methods based on the adiabatic-connection fluctuation-dissipation theorem Andreas Görling, Patrick Bleiziffer, Daniel Schmidtel, and Andreas Heßelmann Lehrstuhl für
More informationDivergence in Møller Plesset theory: A simple explanation based on a two-state model
JOURNAL OF CHEMICAL PHYSICS VOLUME 112, NUMBER 22 8 JUNE 2000 Divergence in Møller Plesset theory: A simple explanation based on a two-state model Jeppe Olsen and Poul Jørgensen a) Department of Chemistry,
More informationNo. 2 lectronic state and potential energy function for UH where ρ = r r e, r being the interatomic distance and r e its equilibrium value. How
Vol 12 No 2, February 2003 cfl 2003 Chin. Phys. Soc. 1009-1963/2003/12(02)/0154-05 Chinese Physics and IOP Publishing Ltd lectronic state and potential energy function for UH 2+* Wang Hong-Yan( Ψ) a)y,
More informationBridging Scales Through Wavefunction Analysis
Bridging Scales Through Wavefunction Analysis Felix Plasser Institute for Theoretical Chemistry, University of Vienna Excited States Bridging Scales Marseille, November 7 10, 2016 F. Plasser Wavefunction
More informationTime-dependent density functional theory (TDDFT)
Advanced Workshop on High-Performance & High-Throughput Materials Simulations using Quantum ESPRESSO ICTP, Trieste, Italy, January 16 to 27, 2017 Time-dependent density functional theory (TDDFT) Ralph
More informationExact factorization of the electron-nuclear wave function and the concept of exact forces in MD
Exact factorization of the electron-nuclear wave function and the concept of exact forces in MD E.K.U. Gross Max-Planck Institute for Microstructure Physics Halle (Saale) OUTLINE Thanks Exact factorisation
More informationDensity-Fitting Approximations to the Electron Repulsion Integrals
Density-Fitting Approximations to the Electron Repulsion Integrals C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Created on 19 July 2010 Here we will follow some
More informationLecture 9 Electronic Spectroscopy
Lecture 9 Electronic Spectroscopy Molecular Orbital Theory: A Review - LCAO approximaton & AO overlap - Variation Principle & Secular Determinant - Homonuclear Diatomic MOs - Energy Levels, Bond Order
More informationE x E x HF, while E c of DFT is approximated with the difference E c
Exchange and correlation energy in density functional theory: Comparison of accurate density functional theory quantities with traditional Hartree Fock based ones and generalized gradient approximations
More informationTheoretical study of the low-lying excited singlet states of furan
JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 2 8 JULY 2003 Theoretical study of the low-lying excited singlet states of furan E. V. Gromov, A. B. Trofimov, and N. M. Vitkovskaya Laboratory of Quantum
More information4 Post-Hartree Fock Methods: MPn and Configuration Interaction
4 Post-Hartree Fock Methods: MPn and Configuration Interaction In the limit of a complete basis, the Hartree-Fock (HF) energy in the complete basis set limit (ECBS HF ) yields an upper boundary to the
More informationMolecular Physics. Attraction between the ions causes the chemical bond.
Molecular Physics A molecule is a stable configuration of electron(s) and more than one nucleus. Two types of bonds: covalent and ionic (two extremes of same process) Covalent Bond Electron is in a molecular
More informationSize-extensive wave functions for QMC A linear-scaling GVB approach
Size-extensive wave functions for QMC A linear-scaling GVB approach Claudia Filippi, University of Twente, The Netherlands Francesco Fracchia, University of Pisa, Italy Claudio Amovilli, University of
More informationSpecial 5: Wavefunction Analysis and Visualization
Special 5: Wavefunction Analysis and Visualization Felix Plasser Institute for Theoretical Chemistry, University of Vienna COLUMBUS in China Tianjin, October 10 14, 2016 F. Plasser Wavefunction Analysis
More informationStructure of diatomic molecules
Structure of diatomic molecules January 8, 00 1 Nature of molecules; energies of molecular motions Molecules are of course atoms that are held together by shared valence electrons. That is, most of each
More informationQ-Chem Workshop Examples Part 2
Marek Freindorf Q-Chem Workshop Examples Part 2 Louisville, KY March 2010 CO, MP2 Calculations Open Avogadro using the "Build" button of QUI Page 2 CO, MP2 Calculations Using Avogadro to create coordinates
More informationCHAPTER 11 MOLECULAR ORBITAL THEORY
CHAPTER 11 MOLECULAR ORBITAL THEORY Molecular orbital theory is a conceptual extension of the orbital model, which was so successfully applied to atomic structure. As was once playfuly remarked, a molecue
More informationExcitation energies from time-dependent density-functional theory beyond the adiabatic approximation
JOURNAL OF CHEMICAL PHYSICS VOLUME 11, NUMBER 1 1 JULY 004 Excitation energies from time-dependent density-functional theory beyond the adiabatic approximation C. A. Ullrich a) Department of Physics, University
More informationVariational Methods for Electronic Structure
Variational Methods for Electronic Structure The hydrogen atom is a two-body system consisting of a proton and an electron. If spin and relativistic effects are ignored, then the Schrödinger equation for
More informationOrbital entanglement in bond-formation processes
Orbital entanglement in bond-formation processes Katharina Boguslawski a, Pawel Tecmer a, Gergely Barcza b, Örs Legezab,1, and Markus Reiher a,2 arxiv:10.20v2 [physics.chem-ph] 11 May 201 a ETH Zurich,
More informationExchange-Correlation Functional
Exchange-Correlation Functional Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
More information(1/2) M α 2 α, ˆTe = i. 1 r i r j, ˆV NN = α>β
Chemistry 26 Spectroscopy Week # The Born-Oppenheimer Approximation, H + 2. Born-Oppenheimer approximation As for atoms, all information about a molecule is contained in the wave function Ψ, which is the
More informationTime-dependent density functional theory (TDDFT)
04/05/16 Hands-on workshop and Humboldt-Kolleg: Density-Functional Theory and Beyond - Basic Principles and Modern Insights Isfahan University of Technology, Isfahan, Iran, May 2 to 13, 2016 Time-dependent
More informationElectronic structure of correlated electron systems. Lecture 2
Electronic structure of correlated electron systems Lecture 2 Band Structure approach vs atomic Band structure Delocalized Bloch states Fill up states with electrons starting from the lowest energy No
More informationwould represent a 1s orbital centered on the H atom and φ 2px )+ 1 r 2 sinθ
Physical Chemistry for Engineers CHEM 4521 Homework: Molecular Structure (1) Consider the cation, HeH +. (a) Write the Hamiltonian for this system (there should be 10 terms). Indicate the physical meaning
More informationTime-dependent linear-response variational Monte Carlo.
Time-dependent linear-response variational Monte Carlo. Bastien Mussard bastien.mussard@colorado.edu https://mussard.github.io/ Julien Toulouse julien.toulouse@upmc.fr Sorbonne University, Paris (web)
More informationSupplemental Material: Experimental and Theoretical Investigations of the Electronic Band Structure of Metal-Organic Framework of HKUST-1 Type
Supplemental Material: Experimental and Theoretical Investigations of the Electronic Band Structure of Metal-Organic Framework of HKUST-1 Type Zhigang Gu, a Lars Heinke, a,* Christof Wöll a, Tobias Neumann,
More informationPseudo-Hermitian eigenvalue equations in linear-response electronic-structure theory
1/11 Pseudo-Hermitian eigenvalue equations in linear-response electronic-structure theory Julien Toulouse Université Pierre & Marie Curie and CNRS, 4 place Jussieu, Paris, France Web page: www.lct.jussieu.fr/pagesperso/toulouse/
More informationAdiabatic-connection fluctuation-dissipation density-functional theory based on range separation
Adiabatic-connection fluctuation-dissipation density-functional theory based on range separation Julien Toulouse 1 I. Gerber 2, G. Jansen 3, A. Savin 1, W. Zhu 1, J. Ángyán 4 1 Laboratoire de Chimie Théorique,
More information