Non-adiabatic dynamics with external laser fields
|
|
- Ralph Mosley
- 5 years ago
- Views:
Transcription
1 January 18, 21
2 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 equilateral H 3 Protonated formaldimine, CH 2NH Appl. H + 2 H Appl. LiF LiF 7 Appl. Oxirane Oxirane
3 Mixed Quantum-Classical Dynamics Definitions {r, m} fast (light) degrees of freedom {R,M} slow (heavy) degrees of freedom Mean-field (Ehrenfest-type) dynamics R t t dt E el (t ) Ψ(r, R,t) = Φ(r,t)Ω(R,t)e i Z Z E el (t) = drdr Φ (r,t)ω (R,t)H el Φ(r,t)Ω(R,t) Born-Oppenheimer (BO) and non-adiabatic BO Ψ(r,R, t) = X I Φ I (r,r)ω I (R,t) Φ I (r,r) Solution of the time-independent Schrödinger equation (Only parametrically dependent on R) Other possible Ansätze
4 Quantum Dynamics with Trajectories Molecular Hamiltonian H tot = X γ z } { 2 X 2 2 γ 2 i + V en (r,r) 2M γ 2m i i H Inserting Ψ(r, R,t) = P I Φ I(r,R)Ω I (R,t) in the TDSE Nuclei Ω I (R, t) = A I (R,t)e i S I (R,t) in the limit ( ) ( A J = F(A t J, A J,S J, S J ) + P I J A J S J = P 1 t γ 2M γ ( γs J ) 2 + H II i hṙ d JI + i H IJ e i R t t dτ (H II H JJ ) Motion evolves on each PES with transport of amplitude between PESs governed by the off-diagonal terms d JI and H JI. Electrons are NOT dynamical. Define the PES computed on-the-fly
5 Quantum Dynamics with Trajectories Independent Trajectory Approximation In the ITA each classical trajectory representing the nuclear wavepacket evolve independently from all the other (missing coherence/decoherence effects) and the amplitudes only change through the coupling terms between PESs A J t = X I J A J»Ṙ d JI + i H IJ e i R t t dτ (H II H JJ ) J. Tully 1 has introduced a stochastic algorithm that (in the limit ) decide about the surface that drags the "classical" nuclear dynamics. 1 J.C. Tully, R.K. Preston, JCP, 55, 562 (1971)
6 Quantum Dynamics with Trajectories Tully s Surface Hopping In Tully s dynamics, the classical trajectories evolve adiabatically according to Born-Oppenheimer dynamics until a hop between two potential energy surfaces (H II and H JJ ) occurs according to a Monte Carlo-type procedure. In the "fewest switches" algorithm, the transition probability from state I to state J in the time interval [t, t + dt] is Z t+dt g IJ (t,t + dt) 2 dτ Im[A J(τ)A I H JI (τ)] Re[A J (τ)a I (τ)ξ JI (τ)] t A I (τ)a I (τ), (1) where Ξ JI (τ) = P γ Ṙγ d γ JI (τ) = P γ and only if X R γ K i 1g JK < ζ < X K I where ζ is generated randomly in the interval [, 1] Φ t I R γ Φ J, and a hop occurs if g JK, (2)
7 Quantum Dynamics with Trajectories Tully s Surface Hopping In the ITA a single trajectory follows adiabatically a PES (I) until a transition to a different PES (J) is accepted with a given probability g(a IJ ). The amplidudes A I evolve according to i Ȧ J = X I A I (H JI i Ṙ d JI ) A swarm of independent trajectories will then reproduce the approximate dynamics of the nuclear wavefunction amplitude. Tully s claim The algorithm "guarantee" that at any instant of time for the ENSEMBLE of trajectories the fraction of trajectories assigned to any PES is approximately equal to the relative population of the state A I (t) 2.
8 Quantum Dynamics with Trajectories Protonated formaldimine is a model compound for the study of isomerization in rhodopsin chromophore retinal. Photo-excitation promotes the system mainly into S 2 and therefore the relaxation occurring in the excited states involves at least 3 states: S (GS), S 1 and S 2.
9 Interaction Hamiltonian Interaction Hamiltonian In the Born-Oppenheimer approximation for the separation of electronic and nuclear degrees of freedom (which is assumed in Tully s dynamics), the total (non relativistic) Hamiltonian is given by H tot = H mol + H rad + H int where the interaction Hamiltonian (with no spin-magnetic field contributions) is obtained from the standard prescription p p ea/c, H int = X i» e 2mc (p i A(r i,t) + A(r i,t) p i ) + e2 2mc A(r i,t)a(r 2 i,t). The vector potential is of the form A = A ǫ λ e ik r+iωt
10 Interaction Hamiltonian In a typical molecular transition in the optical region, the wavelength of the interacting photon is much greater than the linear dimension of the molecule, λ = 2ω/k >> r molecule and therefore e ik r = 1 ik r + O(k r) 2 can be replaced by the leading term, 1. The radiation field coupling matrix element becomes (using A = A ǫ λ ) H int JI = iω JI A c µ JIe iωt where µ JI = e Φ J X α r α Φ I is the the transition dipole vector, and ω JI = (E J E I )/.
11 Interaction Hamiltonian Tully s differential equations for the state amplitudes i ȦJ = X I A I (H JI i Ṙ d JI) become i ȦJ = X I A i (H JI i Ṙ d JI + iω JI A c µ JI eiωt ) In addition, a classical electrostatic interaction term of the form E nucl (R α ) = X γ Z γr γ E(t) with E(t) = 1 A(t) is added to the potential energy of the nuclei, which, c t according to Tully s scheme, follow a "classical" trajectory on a single PES with possible surface hops according to the hopping probability g JI.
12 NAC in DFT NAC evaluation For any one-body operator, Ô, a mapping between MBPT and TDDFT quantities gives (we only consider transitions from the ground state Φ ) O S 1/2 e I = ω 1/2 I Φ Ô Φ I where the operator Ô = P iaσ o iaσâ aσâ iσ has components o iaσ = φ iσ Ô ψ aσ 2 with ω I = E I E and S ijσ,klτ = δ ik δ jl δ στ (f kσ f lσ )(ǫ lσ ǫ kσ ), e I are the TDDFT eigenvectors of the pseudoeigenvalue equation, Ωe I = ω 2 Ie I, q Ω ijσ,klτ = δ στ δ ik δ jl (ǫ lτ ǫ kσ ) q(f iσ f jσ )(ǫ jσ ǫ iσ )K ijσ,klτ (f kτ f lτ )(ǫ lτ ǫ kτ ) K is P the matrix form of the TDDFT kernel 2 iaσ stands for P N P σ {α,β}. P i=1 a=1
13 NAC in DFT The operator Ô can be the position operator, ˆr or the gradient R Ĥ ee, or any other one-body operator of interest, and may depend parametrically on the atomic positions, R. 3 For practical purposes we introduce the auxiliary linear-response many-electron wavefunctions 4 as a linear combination of singly excited Slater determinants Φ I [{φ }] = X iaσ c I iaσ â aσâ iσ Φ[{φ }], with c I iaσ s S 1 iaσ ω I e I iaσ where Φ [{φ }] is the Slater determinant of all occupied KS orbitals {φ iσ } N i=1, which, at a turn, are promoted into a virtual (unoccupied) orbitals, ψ aσ. 3 E. Tapavicza, IT, U. Rothlisberger, PRL, 98, 231 (27); C.P. Hu, O. Sugino, JCP, 127, 6413 (27) 4 IT, E. Tapavicza, U. Rothlisberger, JCP, 13, (27); JCP, 131, (29).
14 NAC in DFT The non-adiabatic coupling elements at the mid step t + δt/2 of a TDDFT MD dynamics can therefore be calculated as E Ṙ d I t+δt/2 = D Φ(r;R(t)) R ΦI (r; R(t)) D Φ(r; Ṙ = R(t)) t Φ E I (r;r(t)) 1 h i Φ (r;r(t)) Φ I (r; R(t + δt)) Φ (r; R(t + δt)) Φ I (r;r(t)) 2δt The non-adiabatic coupling vectors Ω I = Φ I (R) R Φ (R) = Φ I (R) R H Φ (R) E I (R) E (R) In second quantization, ˆr = P ia r iaâ aσâ iσ, and r is the vector with elements r ia = φ a ˆr ψ a in the basis of the KS orbitals. The transition dipole matrix elements µ I = eω 1/2 I r S 1/2 e I
15 Sternheimer Implementation Linear response perturbative TDDFT equations (Sternheimer 1951) NX j=1 (H σ δ ij ǫ ij ) φ I,jσ + Q σ δv σscf I φ iσ = ω I φ I,iσ where Q σ = 1 P N i=1 φ iσ φ iσ, H σ is the Kohn-Sham Hamiltonian, ǫ ij is a Lagrangian multiplier ensuring the orthogonality of the ground state orbitals, Z δv σscf I (r, ω) = d 3 r 1 r r + δ 2 E xc δρ σ(r) δρ σ(r )!ρ ρ=ρ Iσ (r, ω) (using ALDA) with ρ Iσ (r, t) = P i φ iσ(r)φ I,iσ (r, t) + φ I,iσ (r, t)φ iσ (r). In this formulation: Φ I [{φ }, {φ I, }] = X iσ ˆr I,iσâiσ Φ [{φ }] where ˆr I,iσ is the creator operator for the linear response orbital φ I,iσ.
16 Implementation The gradient of the Hamiltonian operator in Ω IJ = Φ I (R) R Φ J (R) = Φ I (R) R Ĥ Φ J (R) E I (R) E J (R) is split into the explicit and the implicit terms, R Ĥ[ρ](R) = Ĥ[ρ](R) Z ρ + R dr Ĥ[ρ](R) ρ(r) ρ(r) R, The first one can be computed analytically and describes the contribution from the electrostatic interaction between nuclei and electrons, and from the local and nonlocal parts of the pseudopotentials. The second term instead adds the contribution from the Hartree and exchange-correlation potentials.
17 Implementation There are two possible solutions of the implicit term: Using finite differences. Costly because there are 3 M nuclear displacements (M is the number of atoms). Using Density Functional Perturbation Theory (DFPT) with perturbing functional E pert γ [ρ] = Z drρ(r) V ions,el(r) R γ where the perturbation is described by an infinitesimal nuclear displacement, λ e Aµ, along the cartesian axis centered on each atom. The index A labels the ions and µ the cartesian coordinates along which the displacements occur (µ = 1, 2,3). We use the collective index γ for both indices A and µ. δρ γ(r) = X iσ φ iσ (r)δφ γ,iσ(r) + δφ γ,iσ (r)φ iσ(r)
18 Implementation The linear response-density is then used to compute the matrix elements of the implicit term of R Ĥ[ρ](R) according to Φ I X σ Z dr Ĥ[ρ](R) ρ σ(r) ρ σ(r) R Φ J = Φ I X γ X Z e γ σ dr Ĥ[ρ](R) ρ σ(r) where Φ J is either the ground state wavefunction (J = ) or 5 Φ J = NX X ˆr J,iσâJ,iσ Φ i=1 σ ρ γσ(r) Φ J 5 IT, E. Tapavicza, U. Rothlisberger, JCP, 13, 1241 (29)
19 Nearly equilateral H 3 Energy [ev] Ω 1 component [Bohr -1 ] R [Å] d12 x H 2 H 1 d 3 Upper panel: Adiabatic potential energy surfaces of the ground state and first singlet excited state of H 3 calculated along the path depicted here above. Solid lines: TDDFT/ALDA; Dashed lines: reference. lower panel: Corresponding z-component of the NACVs on the approaching H atom. The ALDA/TDDFT results (circles) are shown together with some reference data. Un-rescaled TDDFT/ALDA z components of the NACVs are also shown (crosses). a a IT, E. Tapavicza, U. Rothlisberger, JCP, 13, 1241 (29); THEOCHEM, 914, 22 (29). χ H 3 z
20 Protonated formaldimine, CH 2 NH Energy [ev] 1 5 Ω Ω 1 [Bohr [Bohr -1 ] ] Time evolution of the potential energy surfaces (upper panel) S (solid line), S 1 (dashed line), S 2 (dash dotted line), and S 3 (dotted line). The state that drives the dynamics is shown in red. The main structural change in S 2 during the first 1 fs is a CN bond elongation. The middle and lower panel show the corresponding time series of the absolute value of the NACVs of the carbon (black) and nitrogen (red) atoms computed for the pair of states S /S 1 and S 1/S 2, respectively fs
21 H + 2 Photoinduced excitation of H A.J. Cohen, P. Mori-Sanchez, W. Yang, Science, 321, 792 (28)
22 H H 2 + Rabi oscillations.8.6 C i fs
23 H + 2 H 2 + DFT BLYP Energy [a.u.] c 2 c r [Å] fs
24 H + 2 H 2 + DFT BLYP Energy [a.u.] c 2 c fs fs
25 H + 2 H 2 + HF Energy [a.u.] C 2 C r [Å] fs
26 H + 2 H 2 + HF Energy [a.u.] C 2 C fs fs
27 LiF Photoinduced excitation of LiF
28 LiF 1 LiF (1,,) - Without hops.8 C(t) GS S1 S2 S Time [fs]
29 LiF 1 LiF Pulse (1,,) - High intensity C(t)^ Li-F [Angstrom] Time [fs]
30 LiF C(t) 2 Li-F [Å] S1 S2 S3 GS LiF Pulse (1,1,1) A(t) Time [fs]
31 LiF LiF Pulse (1,1,1) Energy [hartree] Time [fs]
32 LiF 1 LiF Pusle (1,1,1) - High intensity - Out of phase (omega =.5 a.u.).8 C(t) Li-F [Å] Time [fs]
33 Oxirane Photoinduced excitation of Oxirane
34 Oxirane
35 Oxirane xxx
36 Oxirane Energy [hartree] C(t) S1 S2 S3 S4 GS Oxirane Succession of short pulses A(t) Time [fs]
37 Oxirane Acknowledgments Dr. Enrico Tapavicza (University of California Irvine) Basile Curchod (EPFL) Ursula Roethlisberger (EPFL) Also thanks to N. Maitra for introducing me to the photo-dissociation H + 2 system.
Time-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 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 informationElectronic Structure Methods for molecular excited states
Electronic Structure Methods for molecular excited states Summer School on Time dependent Density-Functional Theory: Prospects and Applications Benasque, Spain September 2008 Ivano Tavernelli EPFL CH-1015
More informationReminder from last lecture: potential energy surfaces
Outline 1 Why Quantum Dynamics? A TDDFT approach Ivano Tavernelli IBM -Research Zurich Summer School of the Max-Planck-EPFL Center for Molecular Nanoscience & Technology 2 Nonadiabatic Bohmian dynamics
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 informationMixed quantum-classical dynamics. Maurizio Persico. Università di Pisa Dipartimento di Chimica e Chimica Industriale
Mixed quantum-classical dynamics. Maurizio Persico Università di Pisa Dipartimento di Chimica e Chimica Industriale Outline of this talk. The nuclear coordinates as parameters in the time-dependent Schroedinger
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 informationModels for Time-Dependent Phenomena
Models for Time-Dependent Phenomena I. Phenomena in laser-matter interaction: atoms II. Phenomena in laser-matter interaction: molecules III. Model systems and TDDFT Manfred Lein p.1 Outline Phenomena
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationOVERVIEW 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 informationAb 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 informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation
More informationAb 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 informationNonadiabatic 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 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 informationMolecular Dynamics. Park City June 2005 Tully
Molecular Dynamics John Lance Natasa Vinod Xiaosong Dufie Priya Sharani Hongzhi Group: August, 2004 Prelude: Classical Mechanics Newton s equations: F = ma = mq = p Force is the gradient of the potential:
More informationMD 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 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 informationWhat 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 informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » To calculate all the properties of a molecule or crystalline system knowing its atomic information: Atomic species Their coordinates The Symmetry
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 LECTURE II Review of ground-state DFT LECTURE
More informationTheoretical 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 informationKey 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 informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
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) Hamiltonian for
More informationAN ACCELERATED SURFACE-HOPPING METHOD FOR COMPUTATIONAL SEMICLASSICAL MOLECULAR DYNAMICS. Laren K. Mortensen
AN ACCELERATED SURFACE-HOPPING METHOD FOR COMPUTATIONAL SEMICLASSICAL MOLECULAR DYNAMICS by Laren K. Mortensen A senior thesis submitted to the faculty of Brigham Young University in partial fulfillment
More informationYingwei 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 information1 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 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 informationRabi 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 informationTDDFT 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 informationElectron 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 informationLecture 10. Born-Oppenheimer approximation LCAO-MO application to H + The potential energy surface MOs for diatomic molecules. NC State University
Chemistry 431 Lecture 10 Diatomic molecules Born-Oppenheimer approximation LCAO-MO application to H + 2 The potential energy surface MOs for diatomic molecules NC State University Born-Oppenheimer approximation
More informationBorn-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 informationSecond quantization. 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 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 informationAndré 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 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 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 informationElectronic Structure of Clusters and. Nanocrystals
Electronic Structure of Clusters and Nanocrystals James R. Chelikowsky and Yousef Saad Department of Chemical Engineering and Materials Science Department of Computer Science Institute for the Theory of
More informationMolecular Simulation I
Molecular Simulation I Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationMixing Quantum and Classical Mechanics: A Partially Miscible Solution
Mixing Quantum and Classical Mechanics: A Partially Miscible Solution R. Kapral S. Nielsen A. Sergi D. Mac Kernan G. Ciccotti quantum dynamics in a classical condensed phase environment how to simulate
More informationAn Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method
Fakultät für Mathematik und Naturwissenschaften - Lehrstuhl für Physikalische Chemie I / Theoretische Chemie An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method Jan-Ole Joswig
More informationTheoretical 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 informationConical 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 information2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
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 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 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 informationTime-Dependent Electron Localization Function! (TD-ELF)! GOAL! Time-resolved visualization of the breaking and formation of chemical bonds.!
Time-Dependent Electron Localization Function! (TD-ELF)! GOAL! Time-resolved visualization of the breaking and formation of chemical bonds.! How can one give a rigorous mathematical meaning to chemical
More informationThe interaction of light and matter
Outline The interaction of light and matter Denise Krol (Atom Optics) Photon physics 014 Lecture February 14, 014 1 / 3 Elementary processes Elementary processes 1 Elementary processes Einstein relations
More informationDensity 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 informationThe frequency-dependent Sternheimer equation in TDDFT
The frequency-dependent Sternheimer equation in TDDFT A new look into an old equation Miguel A. L. Marques 1 Centre for Computational Physics, University of Coimbra, Portugal 2 LPMCN, Université Claude
More informationReactivity 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 informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part I Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationIntroduction 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 informationComputational 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 informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationJack Simons Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
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 informationSolid 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 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 informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 10.1038/NPHYS2210 Femtosecond torsional relaxation Theoretical methodology: J. Clark, S. Tretiak, T. Nelson, G. Cirmi & G. Lanzani To model non-adiabatic excited state dynamics
More informationAN 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 informationExchange-Correlation Functionals
Exchange-Correlation Functionals E.K.U. Gross Max-Planck Institute of Microstructure Physics Halle (Saale) The Hamiltonian of Condensed Matter (atoms, molecules, clusters, solids) Hamiltonian for the complete
More informationChem 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 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 informationTime-dependent density functional theory
Annual Review of Physical Chemistry 2004 1056-8700/97/0610-00 Time-dependent density functional theory M. A. L. Marques and E.K. U. Gross Institut für Theoretische Physik, Freie Universität Berlin, Arnimallee
More informationDynamics Simulation of Excited State Intramolecular Proton Transfer. Felix Plasser Institut für Theoretische Chemie Universität Wien
Dynamics Simulation of Excited State Intramolecular Proton Transfer Felix Plasser Institut für Theoretische Chemie Universität Wien Wien, am 26. 3. 2009 Contents 1 Introduction 1 2 Theory 5 2.1 The Schrödinger
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 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 informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Raman Centre for Atomic, Molecular and Optical
More informationI. CSFs Are Used to Express the Full N-Electron Wavefunction
Chapter 11 One Must be Able to Evaluate the Matrix Elements Among Properly Symmetry Adapted N- Electron Configuration Functions for Any Operator, the Electronic Hamiltonian in Particular. The Slater-Condon
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
More informationAdvanced TDDFT I: Memory and Initial-State Dependence
Advanced TDDFT I: Memory and Initial-State Dependence when the adiabatic approximation commits a crime Where were you at the time the photon was annihilated? I..um I just can t remember! V xc( r,t ) n(r,t)
More informationEnergy and Forces in DFT
Energy and Forces in DFT Total Energy as a function of nuclear positions {R} E tot ({R}) = E DF T ({R}) + E II ({R}) (1) where E DF T ({R}) = DFT energy calculated for the ground-state density charge-density
More information2m 2 Ze2. , where δ. ) 2 l,n is the quantum defect (of order one but larger
PHYS 402, Atomic and Molecular Physics Spring 2017, final exam, solutions 1. Hydrogenic atom energies: Consider a hydrogenic atom or ion with nuclear charge Z and the usual quantum states φ nlm. (a) (2
More informationModeling Ultrafast Deactivation in Oligothiophenes via Nonadiabatic Dynamics
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2015 Supplementary Data for Modeling Ultrafast Deactivation in Oligothiophenes via Nonadiabatic
More informationIntroduction 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 informationAdvanced Electronic Structure Theory Density functional theory. Dr Fred Manby
Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)
More informationAnswers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014
Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde
More informationMagnetism. Eric Bousquet. University of Liège. Abinit School, Lyon, France 16/05/2014
Magnetism University of Liège eric.bousquet@ulg.ac.be Abinit School, Lyon, France 16/05/2014 Outline Origin of magnetism: Inside an atom Between 2 atoms Interaction with ligands Interaction through ligands
More informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More informationTime-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 informationCHEM-UA 127: Advanced General Chemistry I
1 CHEM-UA 127: Advanced General Chemistry I I. OVERVIEW OF MOLECULAR QUANTUM MECHANICS Using quantum mechanics to predict the chemical bonding patterns, optimal geometries, and physical and chemical properties
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 informationLecture 10. Transition probabilities and photoelectric cross sections
Lecture 10 Transition probabilities and photoelectric cross sections TRANSITION PROBABILITIES AND PHOTOELECTRIC CROSS SECTIONS Cross section = = Transition probability per unit time of exciting a single
More informationNMR and IR spectra & vibrational analysis
Lab 5: NMR and IR spectra & vibrational analysis A brief theoretical background 1 Some of the available chemical quantum methods for calculating NMR chemical shifts are based on the Hartree-Fock self-consistent
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 informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationNoncollinear spins in QMC: spiral Spin Density Waves in the HEG
Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they
More informationPseudopotentials 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 informationMolecular energy levels
Molecular energy levels Hierarchy of motions and energies in molecules The different types of motion in a molecule (electronic, vibrational, rotational,: : :) take place on different time scales and are
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 informationChapter 2 Quantum chemistry using auxiliary field Monte Carlo
Chapter 2 Quantum chemistry using auxiliary field Monte Carlo 1. The Hubbard-Stratonovich Transformation 2. Neuhauser s shifted contour 3. Calculation of forces and PESs 4. Multireference AFMC 5. Examples
More informationNonadiabatic Dynamics: Mean-Field and Surface Hopping
John von Neumann Institute for Computing Nonadiabatic Dynamics: Mean-Field and Surface Hopping Nikos L. Doltsinis published in Quantum Simulations of Complex Many-Body Systems: From Theory to Algorithms,
More informationInstitut Néel Institut Laue Langevin. Introduction to electronic structure calculations
Institut Néel Institut Laue Langevin Introduction to electronic structure calculations 1 Institut Néel - 25 rue des Martyrs - Grenoble - France 2 Institut Laue Langevin - 71 avenue des Martyrs - Grenoble
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 informationQuantum Electrodynamical TDDFT: From basic theorems to approximate functionals
Quantum Electrodynamical TDDFT: From basic theorems to approximate functionals Ilya Tokatly NanoBio Spectroscopy Group - UPV/EHU San Sebastiàn - Spain IKERBASQUE, Basque Foundation for Science - Bilbao
More information