Computational Photochemistry

Size: px
Start display at page:

Download "Computational Photochemistry"

Transcription

1 Computational Photochemistry Filipp Furche Department of Chemistry February 09, 2015 Turbomole Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

2 Acknowledgments Jefferson Bates (Temple) Henk Eshuis (Montclair) CHE Filipp Furche (UC Irvine) DE-SC Computational Photochemistry Enrico Tapavicza (CSULB) Turbomole 02/09/ / 27

3 Introduction Photocatalytic Water Splitting Enrico Berardo, Martijn Zwijnenburg (UCL) Can photochemical water splitting efficiently be catalyzed by TiO nanoclusters? (TiO 2 ) 4 +4 H 2 O model Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

4 Introduction Wish List for Computational Photochemistry 1 Enable calculations on systems with 100s of atoms, high density of excited states, multiple intersections, and near-degeneracies 2 Automatically explore relevant parts of configuration space 3 Account for nuclear quantum effects, such as non-radiative decay 4 Describe bond breaking 5 Require minimal user input ( black box ) Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

5 Introduction Outline 1 Introduction 2 Non-Adiabatic Dynamics 3 Time-Dependent Density Functional Theory 4 Implementation 5 Applications 6 Conclusions Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

6 Non-Adiabatic Dynamics Born-Oppenheimer Approximation Central assumption: Due to mass difference, nuclear kinetic energy is small compared to electronic energy differences. 1 Solve electronic Schrödinger equation for fixed nuclear positions R (adiabatic separation): Ĥ el (x R)Ψ n (x R) = E n (R)Ψ n (x R) 2 Solve nuclear Schrödinger equation using Born-Oppenheimer (BO) potential energy surface: ( ˆT nuc + E n (R))Φ kn (R) = E kn Φ kn (R) Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

7 Non-Adiabatic Dynamics Breakdown of Born-Oppenheimer BO breaks down if electronic excitation energies are comparable to nuclear kinetic energy Close to avoided crossings and conical intersections BO cannot describe non-radiative transitions and excited state decay C. Xie, J. Ma, X. Zhu, D. H. Zhang, D. R. Yarkony, D. Xie, H. Guo, J. Phys. Chem. Lett. 5 (2014), 1055 Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

8 Non-Adiabatic Dynamics Tully Surface Hopping MD Born-Oppenheimer expansion of time-dependent wavefunction: Ψ(t R) = n C n (t) Ψ n (R) Time-dependent Schrödinger equation: iċ = HC iσ(t)c Coupling matrix σ mn = Ṙ τ mn, non-adiabatic coupling τ mn = Ψ m ξ Ψ n Nuclei propagated classically on one surface using Verlet method Surface hop occurs if g mn = 2Re J. C. Tully, J. Chem. Phys. 93 (1990), t+ t t dt C nσ nm C m C n C m > rnd(0, 1) Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

9 Non-Adiabatic Dynamics Building Blocks for Non-Adiabatic Dynamics Excitation energies F. Furche, D. Rappoport in Computational Photochemistry, Ed. M. Olivucci, Elsevier, Amsterdam, 2005, p. 93 Analytical energy gradients F. Furche, R. Ahlrichs, J. Chem. Phys. 117 (2002), 7433 Analytical ground-to-excited state couplings R. Send, F. Furche, J. Chem. Phys. 132 (2010), Analytical excited-to-excited state couplings Q. Ou, G. D. Bellchambers, F. Furche, J. E. Subotnik, J. Chem. Phys., in press (2015) Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

10 Time-Dependent Density Functional Theory Time-Dependent Density Functional Theory Map interacting problem onto non-interacting time-dependent Kohn-Sham (TDKS) system: i t φ j(t) = ĤKS [ρ](t) φ j (t), ρ(t, x) = j φ j (t, x) 2 Time-dependent density uniquely determines exchange correlation (XC) potential: Ĥ KS [ρ](t) = ˆT + ˆV H [ρ](t) + ˆV XC [ρ](t) + ˆV en + V nn + ˆV ext (t) Adiabatic approximation: V XC [ρ](t, x) = δe XC [n] δn(x) n(x)=ρ(t,x) E. Runge, E. K. U. Gross, Phys. Rev. Lett. 52 (1984), 997 Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

11 Time-Dependent Density Functional Theory The TDKS Eigenvalue Problem Interacting excitation energies and transition densities are accessible from linear response of TDKS system Sympletic RPA eigenvalue problem/ Casida s Equations: [( ) ( )] ( ) A B 1 0 Xn Ω B A n = RPA pseudo-norm: X T n X n Y T n Y n = 1 Eigenvalues Ω n are excitation energies, eigenvectors (X n, Y n ) T are one-particle transition density matrices Dimension is N ph = N occ N virt N 2 Y n A. D. McLachlan, M. A. Ball, Rev. Mod. Phys. 36 (1964), 884 M. E. Casida, in Recent advances in density functional methods, Vol. 1, Ed. D. P. Chong, World Scientific, Singapore, 1995, 155 Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

12 Time-Dependent Density Functional Theory Orbital Rotation Hessians Hybrid TDDFT Orbital Rotation Hessians: (A + B) iajb = (ɛ a ɛ i )δ ij δ ab + 2(ia jb) + 2f XC++ iajb c x [(ib ja) + (ij ab)] (A B) iajb = (ɛ a ɛ i )δ ij δ ab + 2f XC iajb + c x [(ib ja) (ij ab)] Generalized XC kernels: Density matrix derivatives f XC±± µνκλ = 2 E XC D ± µν D ± κλ Special cases: TDHF (E XC = 0, c x = 1), TDA/CIS (B = 0), non-hybrid TDDFT (c x = 0) J. E. Bates, F. Furche, J. Chem. Phys. 137 (2012), Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

13 Time-Dependent Density Functional Theory Locally Range-Separated Hybrids Locally range-separated PBE: E XC lrsh-pbe = E X sr-pbe + E X lr + E C PBE Long-range exchange: E X lr = σ Brandon Krull d 3 r 1 d 3 r 2 γ σ (r 1, r 2 ) 2 1 e ασ(r) r 1 r 2 2, r 1 r 2 where α σ (R) = c ρ σ (R) /ρ σ (R) Efficient semi-numerical integration A. V. Krukau, G. E. Scuseria, J. P. Perdew, A. Savin, J. Chem. Phys. 129 (2008), Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

14 Time-Dependent Density Functional Theory Locally Range-Separated Hybrids No tuning necessary Correct high density scaling possible Minimal self-interaction error Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

15 Implementation Generalized Symplectic Davidson Method An alternative which has not been tried is to allow B to be non-orthogonal and to solve the generalized eigenvalue problem in step 1. [E. R. Davidson, Comput. Phys. Commun. 53 (1989), 49] Minimize excitation energy on an iteratively expanded Krylov subspace using raw preconditioned residual vectors as basis No orthogonalization necessary Since residual norm decreases, (much) fewer integrals needed for later iterations Makes linear scaling exchange methods competitive Favorably combined with RI-J approximation Numerically stable Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

16 Implementation Performance System Sym. N at Method Property Speedup Rel. Error 1 C B3LYP/SV(P) α iso % 2 T 282 PBE0/SV(P) E(2 1 A) % 3 C 1 67 PBE0/SV(P) E(2 1 A) % Relative speedups of new algorithm using link+ri on 1 Intel X5560 core (2.80GHz) Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

17 Applications SH-TDDFT Simulations on TiO 2 -Water Clusters Enrico Berardo, Martijn Zwijnenburg (UCL) Simulations show heterolytic water splitting in S 1 state B3LYP/SV(P), 278 K Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

18 Applications Vitamin D Photochemistry Provitamin D (Pro) Lumisterol (Lumi) R:=C 9 H 17 (D 2 ) R:=C 8 H 17 (D 3 ) R:=CH 3 (this study) Previtamin D (Pre) Vitamin D (Vita) Toxisterols (Toxi) Tachysterol (Tachy) E. Tapavicza, A. M. Meyer, F. Furche, Phys. Chem. Chem. Phys. 13 (2011), Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

19 Applications Provitamin D Ring-Opening Typical trajectory: Dark (S 2 ) state is not involved Consistent with recent experimental detection of fluorescence K.-C. Tang, A. Rury, M. B. Orozco, J. Egendorf, K. G. Spears, R. J. Sension, J. Chem. Phys. 134 (2011), Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

20 Applications Excited State Decay Rates Comparison of cyclohexadiene (black), α-terpinene (orange), provitamin D (red) Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

21 Applications Origin of Biexponential Decay Comparison of lifetimes (fs) to experiment PBE0 Exp. Cyclohexadiene τ τ α-terpinene τ τ Provitamin D τ τ Biexponential decay caused by productive/non-productive pathway At 300 K, τ(pro) 5τ(CHD) K. Kosma, S. A. Trushin, W. Fuss, W. E. Schmid, Phys. Chem. Chem. Phys. 11 (2009), 172 K.-C. Tang, A. Rury, M. B. Orozco, J. Egendorf, K. G. Spears, R. J. Sension, J. Chem. Phys. 134 (2011), Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

22 Applications Quantum Yields Quantum yields for ring-opening (%): PBE0 Exp. Cyclohexadiene (Cyclohexane) α-terpinene 54 Provitamin D (EtOH) Lumisterol (EtOH) K. Kosma, S. A. Trushin, W. Fuss, W. E. Schmid, Phys. Chem. Chem. Phys. 11 (2009), 172 K.-C. Tang, A. Rury, M. B. Orozco, J. Egendorf, K. G. Spears, R. J. Sension, J. Chem. Phys. 134 (2011), Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

23 Applications Previtamin D Photochemistry Tachy formation: Z/E isomerization by Hula-twist mechanism Quantum yields (in %): Product Computed Exp. 303 nm 254 nm Pre 69 Lumi Pro Tachy Toxi Vita Total simulation time 0.2 ns(!), few weeks on 10 cluster nodes Embarrassingly parallel H. J. C. Jacobs, J. W. J. Gielen, E. Havinga, Tetrahedron Lett. 22 (1981), 4013 Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

24 Applications Acetaldehyde Photochemistry Dissociation into open shell fragments is common Conventional wisdom: Multi-reference methods needed K. H. J. Giesbertz, E. J. Baerends, Chem. Phys. Lett. 461 (2008), 338 Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

25 Applications Spin-Unrestricted Non-Adiabatic TDDFT Simulations Typical bond-breaking trajectory: PBE0/SV(P), 10 a.u. timestep, 157 nm excitation Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

26 Conclusions Conclusions and Outlook TDDFT-SH enables non-adiabatic MD simulations of molecules with up to 100 atoms Mechanisms, quantum yields, product distributions are obtained with qualitative accuracy Hybrid functionals are necessary Re-formulation of Davidson method in conjunction with RI-J yields speedup of up to 5 for hybrid functionals Non-adiabatic MD simulations using single-reference TDDFT can describe homolytic bond breaking Next steps: Double excitations, intersystem crossing, massively parallel implementation Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

27 Conclusions Further Reading E. Tapavicza, G. Bellchambers, J. C. Vincent, F. Furche, Phys. Chem. Chem. Phys. 15 (2013), F. Furche, R. Ahlrichs, C. Hättig, M. Klopper, F. Weigend M. Sierka, WIREs Comput. Mol. Sci. 4 (2014), 91 Furche group publication list: Filipp Furche (UC Irvine) Computational Photochemistry 02/09/ / 27

TDDFT as a tool in biophysics

TDDFT 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 information

TDDFT in Chemistry and Biochemistry III

TDDFT 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 information

Exchange 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 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 information

Conical Intersections. Spiridoula Matsika

Conical Intersections. Spiridoula Matsika Conical Intersections Spiridoula Matsika The Born-Oppenheimer approximation Energy TS Nuclear coordinate R ν The study of chemical systems is based on the separation of nuclear and electronic motion The

More information

Key concepts in Density Functional Theory (I) Silvana Botti

Key concepts in Density Functional Theory (I) Silvana Botti From the many body problem to the Kohn-Sham scheme European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre

More information

Linear response time-dependent density functional theory

Linear 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 information

Density Functional Theory for Electrons in Materials

Density 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 information

Time-Dependent Density-Functional Theory

Time-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 information

VALENCE Hilary Term 2018

VALENCE 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 information

Short Course on Density Functional Theory and Applications VIII. Time-dependent DFT

Short 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 information

Adiabatic connection for near degenerate excited states

Adiabatic 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 information

Convergence behavior of RPA renormalized many-body perturbation theory

Convergence behavior of RPA renormalized many-body perturbation theory Convergence behavior of RPA renormalized many-body perturbation theory Understanding why low-order, non-perturbative expansions work Jefferson E. Bates, Jonathon Sensenig, Niladri Sengupta, & Adrienn Ruzsinszky

More information

Introduction to Computational Chemistry

Introduction to Computational Chemistry Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry room B430, Chemicum 4th floor vesa.hanninen@helsinki.fi September 3, 2013 Introduction and theoretical backround September

More information

Théorie de la fonctionnnelle de la densité avec séparation de portée pour les forces de van der Waals

Théorie de la fonctionnnelle de la densité avec séparation de portée pour les forces de van der Waals Théorie de la fonctionnnelle de la densité avec séparation de portée pour les forces de van der Waals Julien Toulouse 1 Iann Gerber 2, Georg Jansen 3, Andreas Savin 1, János Ángyán 4 1 Laboratoire de Chimie

More information

A 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 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 information

Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems

Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems Mitglied der Helmholtz-Gemeinschaft Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems Birkbeck University, London, June the 29th 2012 Edoardo Di Napoli Motivation and Goals

More information

Introduction to Density Functional Theory

Introduction to Density Functional Theory Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary

More information

CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122

CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122 CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 191 THANKS TO MANY COLLABORATORS, INCLUDING SY VOSKO DAVID LANGRETH ALEX

More information

Ab initio Molecular Dynamics Born Oppenheimer and beyond

Ab initio Molecular Dynamics Born Oppenheimer and beyond Ab initio Molecular Dynamics Born Oppenheimer and beyond Reminder, reliability of MD MD trajectories are chaotic (exponential divergence with respect to initial conditions), BUT... With a good integrator

More information

First-principles Molecular Dynamics Simulations

First-principles Molecular Dynamics Simulations First-principles Molecular Dynamics Simulations François Gygi University of California, Davis fgygi@ucdavis.edu http://eslab.ucdavis.edu http://www.quantum-simulation.org MICCoM Computational School, Jul

More information

Density Functional Theory. Martin Lüders Daresbury Laboratory

Density Functional Theory. Martin Lüders Daresbury Laboratory Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei

More information

Nonadiabatic Molecular Dynamics with Kohn-Sham DFT: Modeling Nanoscale Materials. Oleg Prezhdo

Nonadiabatic Molecular Dynamics with Kohn-Sham DFT: Modeling Nanoscale Materials. Oleg Prezhdo Nonadiabatic Molecular Dynamics with Kohn-Sham DFT: Modeling Nanoscale Materials Oleg Prezhdo CSCAMM Mar 11, 2010 Outline Nonadiabatic MD with Kohn-Sham DFT Advantages & Validity Quantum Backreaction &

More information

Rabi oscillations within TDDFT: the example of the 2 site Hubbard model

Rabi oscillations within TDDFT: the example of the 2 site Hubbard model Rabi oscillations within TDDFT: the example of the 2 site Hubbard model Johanna I. Fuks, H. Appel, Mehdi,I.V. Tokatly, A. Rubio Donostia- San Sebastian University of Basque Country (UPV) Outline Rabi oscillations

More information

Density-Fitting Approximations to the Electron Repulsion Integrals

Density-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 information

XYZ of ground-state DFT

XYZ 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 information

Pseudopotentials for hybrid density functionals and SCAN

Pseudopotentials for hybrid density functionals and SCAN Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania Why do we need pseudopotentials?

More information

Theoretical Photochemistry WiSe 2016/17

Theoretical Photochemistry WiSe 2016/17 Theoretical Photochemistry WiSe 2016/17 Lecture 8 Irene Burghardt burghardt@chemie.uni-frankfurt.de) http://www.theochem.uni-frankfurt.de/teaching/ Theoretical Photochemistry 1 Topics 1. Photophysical

More information

(1/2) M α 2 α, ˆTe = i. 1 r i r j, ˆV NN = α>β

(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 information

Time-dependent density functional theory (TDDFT)

Time-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 information

Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory

Electronic 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 information

Post Hartree-Fock: MP2 and RPA in CP2K

Post Hartree-Fock: MP2 and RPA in CP2K Post Hartree-Fock: MP2 and RPA in CP2K A tutorial Jan Wilhelm jan.wilhelm@chem.uzh.ch 4 September 2015 Further reading MP2 and RPA by Mauro Del Ben, Jürg Hutter, Joost VandeVondele Del Ben, M; Hutter,

More information

CHEM6085: Density Functional Theory Lecture 10

CHEM6085: Density Functional Theory Lecture 10 CHEM6085: Density Functional Theory Lecture 10 1) Spin-polarised calculations 2) Geometry optimisation C.-K. Skylaris 1 Unpaired electrons So far we have developed Kohn-Sham DFT for the case of paired

More information

Key concepts in Density Functional Theory

Key concepts in Density Functional Theory From the many body problem to the Kohn-Sham scheme ILM (LPMCN) CNRS, Université Lyon 1 - France European Theoretical Spectroscopy Facility (ETSF) December 12, 2012 Lyon Outline 1 The many-body problem

More information

Intermediate 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 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 information

Electron States of Diatomic Molecules

Electron 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 information

Module 6 1. Density functional theory

Module 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 information

Two-Component Ehrenfest Dynamics. Li Research Group. University of Washington, Seattle

Two-Component Ehrenfest Dynamics. Li Research Group. University of Washington, Seattle Two-Component Ehrenfest Dynamics Li Research Group University of Washington, Seattle Outline Quantum Electronic Dynamics Background Decay of Coherent-Exciton Time-Dependent Two-Component HF/DFT Theory

More information

Assessment of range-separated time-dependent density-functional theory for calculating C 6 dispersion coefficients

Assessment 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 information

Time-dependent density functional theory (TDDFT)

Time-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 information

Dispersion energies from the random-phase approximation with range separation

Dispersion energies from the random-phase approximation with range separation 1/34 Dispersion energies from the random-phase approximation with range separation Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Web page: www.lct.jussieu.fr/pagesperso/toulouse/

More information

Non-adiabatic dynamics with external laser fields

Non-adiabatic dynamics with external laser fields January 18, 21 1 MQCD Mixed Quantum-Classical Dynamics Quantum Dynamics with Trajectories 2 Radiation Field Interaction Hamiltonian 3 NAC NAC in DFT 4 NAC computation Sternheimer Implementation Nearly

More information

What s the correct classical force on the nuclei

What s the correct classical force on the nuclei What s the correct classical force on the nuclei E.K.U. Gross Max-Planck Institute of Microstructure Physics Halle (Saale) What s the correct classical force on the nuclei or: How to make the Born-Oppenheimer

More information

How 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 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 information

Fractional-charge and Fractional-spin errors in Range-Separated Density-Functional Theory with Random Phase Approximation.

Fractional-charge and Fractional-spin errors in Range-Separated Density-Functional Theory with Random Phase Approximation. Fractional-charge and Fractional-spin errors in Range-Separated Density-Functional Theory with Random Phase Approximation. Bastien Mussard, Julien Toulouse Laboratoire de Chimie Théorique Institut du Calcul

More information

Ab initio molecular dynamics : BO

Ab initio molecular dynamics : BO School on First Principles Simulations, JNCASR, 2010 Ab initio molecular dynamics : BO Vardha Srinivasan IISER Bhopal Why ab initio MD Free of parametrization and only fundamental constants required. Bond

More information

Yingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2

Yingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2 Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer

More information

Electron Correlation

Electron Correlation Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation

More information

DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY

DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES

More information

Pseudopotential generation and test by the ld1.x atomic code: an introduction

Pseudopotential generation and test by the ld1.x atomic code: an introduction and test by the ld1.x atomic code: an introduction SISSA and DEMOCRITOS Trieste (Italy) Outline 1 2 3 Spherical symmetry - I The Kohn and Sham (KS) equation is (in atomic units): [ 1 ] 2 2 + V ext (r)

More information

Chem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.

Chem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components. Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative

More information

Theoretical Photochemistry SoSe 2014

Theoretical Photochemistry SoSe 2014 Theoretical Photochemistry SoSe 2014 Lecture 9 Irene Burghardt (burghardt@chemie.uni-frankfurt.de) http://www.theochem.uni-frankfurt.de/teaching/ Theoretical Photochemistry 1 Topics 1. Photophysical Processes

More information

Fast Eigenvalue Solutions

Fast Eigenvalue Solutions Fast Eigenvalue Solutions Techniques! Steepest Descent/Conjugate Gradient! Davidson/Lanczos! Carr-Parrinello PDF Files will be available Where HF/DFT calculations spend time Guess ρ Form H Diagonalize

More information

Range-separated density-functional theory with long-range random phase approximation

Range-separated density-functional theory with long-range random phase approximation Range-separated density-functional theory with long-range random phase approximation Julien Toulouse 1 Wuming Zhu 1, Andreas Savin 1, János Ángyán2 1 Laboratoire de Chimie Théorique, UPMC Univ Paris 6

More information

1 Density functional theory (DFT)

1 Density functional theory (DFT) 1 Density functional theory (DFT) 1.1 Introduction Density functional theory is an alternative to ab initio methods for solving the nonrelativistic, time-independent Schrödinger equation H Φ = E Φ. The

More information

CHEM6085: Density Functional Theory

CHEM6085: Density Functional Theory Lecture 5 CHEM6085: Density Functional Theory Orbital-free (or pure ) DFT C.-K. Skylaris 1 Consists of three terms The electronic Hamiltonian operator Electronic kinetic energy operator Electron-Electron

More information

Advanced TDDFT II. Double Excitations in TDDFT

Advanced 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 information

Density Functional Theory

Density Functional Theory Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction

More information

Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory

Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National

More information

3/23/2010 More basics of DFT Kieron Burke and friends UC Irvine Physics and Chemistry References for ground-state DFT ABC of DFT, by KB and Rudy Magyar, http://dft.uci.edu A Primer in Density Functional

More information

André Schleife Department of Materials Science and Engineering

André Schleife Department of Materials Science and Engineering André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic

More information

AN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University

AN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,

More information

Computational Methods. Chem 561

Computational Methods. Chem 561 Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational

More information

Born-Oppenheimer Approximation

Born-Oppenheimer Approximation Born-Oppenheimer Approximation Adiabatic Assumption: Nuclei move so much more slowly than electron that the electrons that the electrons are assumed to be obtained if the nuclear kinetic energy is ignored,

More information

Quantum Molecular Dynamics Basics

Quantum Molecular Dynamics Basics Quantum Molecular Dynamics Basics Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological

More information

Introduction to DFTB. Marcus Elstner. July 28, 2006

Introduction to DFTB. Marcus Elstner. July 28, 2006 Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD

More information

OVERVIEW OF QUANTUM CHEMISTRY METHODS

OVERVIEW OF QUANTUM CHEMISTRY METHODS OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density

More information

Basics of DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA

Basics of DFT. Kieron Burke and Lucas Wagner. Departments of Physics and of Chemistry, University of California, Irvine, CA 92697, USA Basics of 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) Basics of DFT Lausanne12 1

More information

Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC

Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Structure of Cement Phases from ab initio Modeling Crystalline C-S-HC Sergey V. Churakov sergey.churakov@psi.ch Paul Scherrer Institute Switzerland Cement Phase Composition C-S-H H Solid Solution Model

More information

The potential of Potential Functional Theory

The potential of Potential Functional Theory The potential of Potential Functional Theory IPAM DFT School 2016 Attila Cangi August, 22 2016 Max Planck Institute of Microstructure Physics, Germany P. Elliott, AC, S. Pittalis, E.K.U. Gross, K. Burke,

More information

Ab initio molecular dynamics: Basic Theory and Advanced Met

Ab initio molecular dynamics: Basic Theory and Advanced Met Ab initio molecular dynamics: Basic Theory and Advanced Methods Uni Mainz October 30, 2016 Bio-inspired catalyst for hydrogen production Ab-initio MD simulations are used to learn how the active site

More information

Large Scale Electronic Structure Calculations

Large Scale Electronic Structure Calculations Large Scale Electronic Structure Calculations Jürg Hutter University of Zurich 8. September, 2008 / Speedup08 CP2K Program System GNU General Public License Community Developers Platform on "Berlios" (cp2k.berlios.de)

More information

ABC of ground-state DFT

ABC of ground-state DFT ABC 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) ABC of ground-state

More information

Density matrix functional theory vis-á-vis density functional theory

Density 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 information

Introduction to density-functional theory. Emmanuel Fromager

Introduction to density-functional theory. 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. Institut

More information

DFT calculations of NMR indirect spin spin coupling constants

DFT calculations of NMR indirect spin spin coupling constants DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling

More information

Quantum Mechanical Simulations

Quantum 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 information

Quantum Physics III (8.06) Spring 2008 Final Exam Solutions

Quantum Physics III (8.06) Spring 2008 Final Exam Solutions Quantum Physics III (8.6) Spring 8 Final Exam Solutions May 19, 8 1. Short answer questions (35 points) (a) ( points) α 4 mc (b) ( points) µ B B, where µ B = e m (c) (3 points) In the variational ansatz,

More information

Supporting Information. I. A refined two-state diabatic potential matrix

Supporting Information. I. A refined two-state diabatic potential matrix Signatures of a Conical Intersection in Adiabatic Dissociation on the Ground Electronic State Changjian Xie, Christopher L. Malbon, # David R. Yarkony, #,* Daiqian Xie,,%,* and Hua Guo,* Department of

More information

Pseudo-Hermitian eigenvalue equations in linear-response electronic-structure theory

Pseudo-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 information

Walter 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 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 information

Electrochemistry project, Chemistry Department, November Ab-initio Molecular Dynamics Simulation

Electrochemistry project, Chemistry Department, November Ab-initio Molecular Dynamics Simulation Electrochemistry project, Chemistry Department, November 2006 Ab-initio Molecular Dynamics Simulation Outline Introduction Ab-initio concepts Total energy concepts Adsorption energy calculation Project

More information

Excited State Processes

Excited State Processes Excited State Processes Photophysics Fluorescence (singlet state emission) Phosphorescence (triplet state emission) Internal conversion (transition to singlet gr. state) Intersystem crossing (transition

More information

Theoretical Photochemistry WiSe 2017/18

Theoretical 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 information

TDDFT II. Alberto Castro

TDDFT II. Alberto Castro Problem set Alberto Castro Institute for Biocomputation and Physics of Complex Systems (BIFI) and Zaragoza Scientific Center for Advanced Modeling (ZCAM), University of Zaragoza, Spain ELK Summit, Lausanne

More information

MD simulation: output

MD simulation: output Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square

More information

What is the course about?

What is the course about? What is the course about? Part 1 Organic Chemistry: organocatalysis 10-10 (1 Å) - 10-9 m Part 2 Structural Biology: peptide conformations Part 3 10-9 - 10-8 m Photochemistry and Photobiology: olefins &

More information

Reactivity and Organocatalysis. (Adalgisa Sinicropi and Massimo Olivucci)

Reactivity and Organocatalysis. (Adalgisa Sinicropi and Massimo Olivucci) Reactivity and Organocatalysis (Adalgisa Sinicropi and Massimo Olivucci) The Aldol Reaction - O R 1 O R 1 + - O O OH * * H R 2 R 1 R 2 The (1957) Zimmerman-Traxler Transition State Counterion (e.g. Li

More information

Time-dependent density functional theory : direct computation of excitation energies

Time-dependent density functional theory : direct computation of excitation energies CECAM Tutorial Electronic excitations and spectroscopies : Theory and Codes Lyon 007 Time-dependent density functional theory : direct computation of excitation energies X. Gonze Université Catholique

More information

Lecture 9: Molecular Orbital theory for hydrogen molecule ion

Lecture 9: Molecular Orbital theory for hydrogen molecule ion Lecture 9: Molecular Orbital theory for hydrogen molecule ion Molecular Orbital Theory for Hydrogen Molecule Ion We have seen that the Schrödinger equation cannot be solved for many electron systems. The

More information

Advanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University

Advanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University Advanced Quantum Chemistry III: Part 3 Haruyuki Nakano Kyushu University 2013 Winter Term 1. Hartree-Fock theory Density Functional Theory 2. Hohenberg-Kohn theorem 3. Kohn-Sham method 4. Exchange-correlation

More information

Introduction to Computational Quantum Chemistry: Theory

Introduction to Computational Quantum Chemistry: Theory Introduction to Computational Quantum Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC 3108 Course Lectures 2007 Introduction Hartree Fock Theory Configuration Interaction Lectures 1 Introduction

More information

Oslo node. Highly accurate calculations benchmarking and extrapolations

Oslo node. Highly accurate calculations benchmarking and extrapolations Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration

More information

Solid State Theory: Band Structure Methods

Solid State Theory: Band Structure Methods Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.

More information

PHOTOCHEMISTRY NOTES - 1 -

PHOTOCHEMISTRY NOTES - 1 - - 1 - PHOTOCHEMISTRY NOTES 1 st Law (Grotthus-Draper Law) Only absorbed radiation produces chemical change. Exception inelastic scattering of X- or γ-rays (electronic Raman effect). 2 nd Law (Star-Einstein

More information

Calculations of electron-molecule scattering cross sections using the Rmatrix method

Calculations of electron-molecule scattering cross sections using the Rmatrix method Calculations of electron-molecule scattering cross sections using the Rmatrix method Jimena Gorfinkiel Department of Physical Sciences The Open University Milton Keynes, UK Jimena.Gorfinkiel@open.ac.uk

More information

Orbital Density Dependent Functionals

Orbital Density Dependent Functionals Orbital Density Dependent Functionals S. Kluepfel1, P. Kluepfel1, Hildur Guðmundsdóttir1 and Hannes Jónsson1,2 1. Univ. of Iceland; 2. Aalto University Outline: Problems with GGA approximation (PBE, RPBE,...)

More information

Adiabatic-connection fluctuation-dissipation density-functional theory based on range separation

Adiabatic-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

The COLUMBUS Project - General Purpose Ab Initio Quantum Chemistry. I. Background and Overview

The COLUMBUS Project - General Purpose Ab Initio Quantum Chemistry. I. Background and Overview The COLUMBUS Project - General Purpose Ab Initio Quantum Chemistry I. Background and Overview Ron Shepard Chemistry Division Argonne National Laboratory CScADS Workshop, Snowbird, Utah, July 23, 2007 Quantum

More information

Extending Kohn-Sham density-functional theory

Extending 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 information