Introduction to linear-response, and time-dependent density-functional theory. Heiko Appel
|
|
- Dylan Nelson
- 5 years ago
- Views:
Transcription
1 Introduction to linear-response, and time-dependent density-functional theory Heiko Appel Fritz-Haber-Institut der Max-Planck-Gesellschaft, Berlin Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
2 Motivation Where is electron dynamics important? Electron-hole pair creation and exciton propagation in solar cells Photosynthesis and energy transfer in light-harvesting antenna complexes Quantum computing (e.g. electronic transitions in ultracold atoms) Molecular electronics, quantum transport... Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
3 Motivation - Decoherence in Quantum Mechanics Today s Google frontpage: 126. birthday of Erwin Schrödinger Schrödinger cat state Ψ = 1 2 (a(t) ψ A + d(t) ψ D ) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
4 Outline Linear Response in DFT Response functions Casida equation Sternheimer equation Real-space representation and real-time propagation Real-space representation for wavefunctions and Hamiltonians Time-propagation schemes Optimal control of electronic motion Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
5 Time-dependent density-functional theory One-to-one correspondence of time-dependent densities and potentials v(r, t) 1 1 ρ(r, t) For fixed inital states, the time-dependent density determines uniquely the time-dependent external potential and hence all physical observables. E. Runge, and E.K.U. Gross, Phys. Rev. Lett. 52, 997 (1984). Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
6 Time-dependent density-functional theory One-to-one correspondence of time-dependent densities and potentials v(r, t) 1 1 ρ(r, t) For fixed inital states, the time-dependent density determines uniquely the time-dependent external potential and hence all physical observables. E. Runge, and E.K.U. Gross, Phys. Rev. Lett. 52, 997 (1984). Time-dependent Kohn-Sham system The time-dependent density of an interacting many-electron system can be calculated as density N ρ(r, t) = ϕ j(r, r) 2 of an auxiliary non-interacting Kohn-Sham system ) i tϕ j(r, t) = ( vs[ρ](r, t) ϕ j(r, t) 2m with a local multiplicative potential v S[ρ(r, t )](r, t) = v ext(r, t) + j=1 ρ(r, t) r r d3 r + v xc[ρ(r, t )](r, t) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
7 Linear Response Theory Hamiltonian Ĥ(t) = Ĥ0 + Θ(t t0)v1(r, t) Initial condition: for times t < t 0 the system is in the ground-state of the unperturbed Hamiltonian Ĥ0 with potential v0 and density ρ0(r) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
8 Linear Response Theory Hamiltonian Ĥ(t) = Ĥ0 + Θ(t t0)v1(r, t) Initial condition: for times t < t 0 the system is in the ground-state of the unperturbed Hamiltonian Ĥ0 with potential v0 and density ρ0(r) For times t > t 0, switch on perturbation v 1(r, t): leads to time-dependent density ρ(r, t) = ρ 0(r) + δρ(rt) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
9 Linear Response Theory Hamiltonian Ĥ(t) = Ĥ0 + Θ(t t0)v1(r, t) Initial condition: for times t < t 0 the system is in the ground-state of the unperturbed Hamiltonian Ĥ0 with potential v0 and density ρ0(r) For times t > t 0, switch on perturbation v 1(r, t): leads to time-dependent density ρ(r, t) = ρ 0(r) + δρ(rt) Functional Taylor expansion of ρ[v](r, t) around v 0: ρ[v](r, t) = ρ[v 0 + v 1](r, t) = ρ[v 0](r, t) δρ[v](rt) + v 1(r t )d 3 r dt δv(r t ) v0 δ 2 ρ[v](rt) + v 1(r t )v 1(r t )d 3 r dt d 3 r dt δv(r t )δv(r t ) v Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
10 Computing Linear Response Different ways to compute first order response in DFT Response functions, Casida equation (frequency-dependent) perturbation theory, Sternheimer equation real-time propagation with weak external perturbation Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
11 Response functions Functional Taylor expansion of ρ[v](r, t) around external potential v 0: δρ[v](rt) ρ[v 0 + v 1](r, t) = ρ[v 0](r) + v 1(r t )d 3 r dt +... δv(r t ) v0 Density-density response function of interacting system χ(rt, r t ) := δρ[v](rt) δv(r t ) v0 Θ(t t ) 0 [ˆρ(r, t) H, ˆρ(r, t ) H] 0 Response of non-interacting Kohn-Sham system: δρ[vs](rt) ρ[v S,0 + v S,1](r, t) = ρ[v S,0](r) + v S(r t )d 3 r dt +... δv S(r t ) v0 Density-density response function of time-dependent Kohn-Sham system χ S(rt, r t ) := δρs[vs](rt) δv S(r t ) vs,0 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
12 Derivation of response equation Definition of time-dependent xc potential v xc(rt) = v KS(rt) v ext(rt) v H(rt) Take functional derivative δv xc(rt) δρ(r t ) = δvks(rt) δρ(r t ) δvext(rt) δρ(r t ) δ(t t ) r r f xc(rt, r t ) := χ 1 S (rt, r t ) χ 1 (rt, r t ) W c(rt, r t ) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
13 Derivation of response equation Definition of time-dependent xc potential v xc(rt) = v KS(rt) v ext(rt) v H(rt) Take functional derivative δv xc(rt) δρ(r t ) = δvks(rt) δρ(r t ) δvext(rt) δρ(r t ) δ(t t ) r r f xc(rt, r t ) := χ 1 S (rt, r t ) χ 1 (rt, r t ) W c(rt, r t ) Act with reponse functions from left and right χ S W c + f xc = χ 1 S χ 1 χ χ S(W c + f xc)χ = χ χ S Dyson-type equation for response functions χ = χ S + χ S(W c + f xc)χ Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
14 First order density response Exact density response to first order In integral notation ρ 1(rt) = + ρ 1 = χv 1 = χ Sv 1 + χ S(W c + f xc)ρ 1 [ d 3 r dt χ S(rt, r t ) v 1(r t ) ] d 3 r dt (W c(r t, r t ) + f xc(r t, r t ))ρ 1(r t ) For practical application: iterative solution with approximate kernel f xc f xc(r t, r t ) = δvxc[ρ](r t ) δρ(r t ) ρ0 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
15 Lehmann representation of linear response function Exact many-body eigenstates Ĥ(t = t 0) m = E m m Lehmann representation of linear density-density response function: χ(r, t;, r, t ) = Θ(t t ) 0 [ˆρ(r, t) H, ˆρ(r, t ) H] 0 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
16 Lehmann representation of linear response function Exact many-body eigenstates Ĥ(t = t 0) m = E m m Lehmann representation of linear density-density response function: χ(r, t;, r, t ) = Θ(t t ) 0 [ˆρ(r, t) H, ˆρ(r, t ) H] 0 Neutral excitation energies are poles of the linear response function! ( 0 ˆρ(r)H m m ˆρ(r χ(r, r ) H 0 ; ω) = lim 0 ˆρ(r ) ) H m m ˆρ(r) H 0 η 0 + ω (E m E m 0) + iη ω + (E m E 0) + iη Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
17 Lehmann representation of linear response function Exact many-body eigenstates Ĥ(t = t 0) m = E m m Lehmann representation of linear density-density response function: χ(r, t;, r, t ) = Θ(t t ) 0 [ˆρ(r, t) H, ˆρ(r, t ) H] 0 Neutral excitation energies are poles of the linear response function! ( 0 ˆρ(r)H m m ˆρ(r χ(r, r ) H 0 ; ω) = lim 0 ˆρ(r ) ) H m m ˆρ(r) H 0 η 0 + ω (E m E m 0) + iη ω + (E m E 0) + iη Exact linear density response to perturbation v 1(ω) ρ 1(ω) = ˆχ(ω)v 1(ω) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
18 Lehmann representation of linear response function Exact many-body eigenstates Ĥ(t = t 0) m = E m m Lehmann representation of linear density-density response function: χ(r, t;, r, t ) = Θ(t t ) 0 [ˆρ(r, t) H, ˆρ(r, t ) H] 0 Neutral excitation energies are poles of the linear response function! ( 0 ˆρ(r)H m m ˆρ(r χ(r, r ) H 0 ; ω) = lim 0 ˆρ(r ) ) H m m ˆρ(r) H 0 η 0 + ω (E m E m 0) + iη ω + (E m E 0) + iη Exact linear density response to perturbation v 1(ω) ρ 1(ω) = ˆχ(ω)v 1(ω) Relation to two-body Green s function i 2 G (2) (r, t;, r, t, r, t;, r, t ) = χ(r, t;, r, t ) + ρ(r)ρ(r ) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
19 Lehmann representation of linear response function Exact many-body eigenstates Ĥ(t = t 0) m = E m m Lehmann representation of linear density-density response function: χ(r, t;, r, t ) = Θ(t t ) 0 [ˆρ(r, t) H, ˆρ(r, t ) H] 0 Neutral excitation energies are poles of the linear response function! ( 0 ˆρ(r)H m m ˆρ(r χ(r, r ) H 0 ; ω) = lim 0 ˆρ(r ) ) H m m ˆρ(r) H 0 η 0 + ω (E m E m 0) + iη ω + (E m E 0) + iη Exact linear density response to perturbation v 1(ω) ρ 1(ω) = ˆχ(ω)v 1(ω) Relation to two-body Green s function i 2 G (2) (r, t;, r, t, r, t;, r, t ) = χ(r, t;, r, t ) + ρ(r)ρ(r ) Current-current response function: Π α,β (r, t;, r, t ) = Θ(t t ) 0 [ĵ α(r, t) H, ĵ β (r, t ) H] 0 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
20 Excitation energies Dyson-type equation for response functions in frequency space [ˆ1 ˆχ S(ω)(Ŵc + ˆf xc(ω))]ρ 1(ω) = χ Sv 1(ω) ρ 1(ω) has poles for exact excitation energies Ω j ρ 1(ω) for ω Ω j On the other hand, rhs χ Sv 1(ω) stays finite for ω Ω j hence the eigenvalues of the integral operator [ˆ1 ˆχ S(ω)(Ŵc + ˆf xc(ω))]ξ(ω) = λ(ω)ξ(ω) vanish, λ(ω) 0 for ω Ω j. Determines rigorously the exact excitation energies [ˆ1 ˆχ S(Ω j)(ŵc + ˆf xc(ω j))]ξ(ω j) = 0 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
21 Casida equation (Non-linear) eigenvalue equation for excitation energies ΩF j = ω 2 j F j with Ω iaσ,jbτ = δ σ,τ δ i,jδ a,b (ɛ a ɛ i) (ɛ a ɛ i)k iaσ,jbτ (ɛb ɛ j) and K iaσ,jbτ (ω) = d 3 r [ d 3 r 1 ] φ iσ(r)φ jσ(r) r r + fxc(r, r, ω) φ kτ (r)φ lτ (r) Eigenvalues ω j are exact vertical excitation energies Eigenvectors can be used to compute oscillator strength Drawback: need occupied and unoccupied orbitals Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
22 Adiabatic approximation Adiabatic approximation: evaluate static Kohn-Sham potential at time-dependent density v adiab xc [ρ](rt) := v static DFT xc [ρ(t)](rt) Example: adiabatic LDA vxc ALDA [ρ](rt) := vxc LDA (ρ(t)) = αρ(r, t) 1/ Exchange-correlation kernel f ALDA xc (rt, r t ) = δvalda xc [ρ](rt) δρ(r t ) = δ(t t )δ(r r ) valda xc ρ(r) = δ(t t )δ(r r ) 2 e hom xc n 2 ρ0 (r) ρ0 (r) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
23 Failures of the adiabatic approximation in linear response H 2 dissociation is incorrect E( 1 Σ + u ) E( 1 Σ + g ) R 0 (in ALDA) Gritsenko, van Gisbergen, Grling, Baerends, JCP 113, 8478 (2000). sometimes problematic close to conical intersections response of long chains strongly overestimated Champagne et al., JCP 109, (1998) and 110, (1999). in periodic solids f xc(q, ω, ρ) = c(ρ), whereas for insulators, q 0 1/q 2 divergent f exact xc charge transfer excitations not properly described Dreuw et al., JCP 119, 2943 (2003). Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
24 Relation of TDHF eigenvalue problem to RPA, CIS and drccd energies RPA equation ( A B B A ) ( X Y ) = ( X Y ) ω RPA correlation energy Ec RPA = 1 T r(ω A) 2 CIS correlation energy from Tamm-Dancoff approximation TDA: B = 0 Ec CIS = 1 T r( ω A) 2 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
25 Relation of TDHF eigenvalue problem to RPA, CIS and drccd energies RPA equation ( A B B A ) ( X Y ) = ( X Y ) ω RPA correlation energy Ec RPA = 1 T r(ω A) 2 CIS correlation energy from Tamm-Dancoff approximation TDA: B = 0 Ec CIS = 1 T r( ω A) 2 Keeping only particle-hole ring contractions, yiels matrix Ricatti equation for CCD cluster amplitudes B + AT + TA + TBT = 0, t ab ij = T ia,jb Correlation energy in direct ring Coupled Cluster Doubles (drccd) E drccd c = 1 2 T r(bt) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
26 Relation of TDHF eigenvalue problem to RPA, CIS and drccd energies RPA equation ( A B B A ) ( X Y ) = ( X Y ) ω Multiplication of RPA equation with X 1 from right yields A + BT = XωX 1, where T := YX 1 Taking trace yields correlation energies 2E drccd c = T r(bt) = T r(xωx 1 A) = T r(ω A) = 2E RPA c Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
27 Relation of TDHF eigenvalue problem to RPA, CIS and drccd energies RPA equation ( A B B A ) ( X Y ) = ( X Y ) ω Multiplication of RPA equation with X 1 from right yields A + BT = XωX 1, where T := YX 1 Taking trace yields correlation energies 2E drccd c = T r(bt) = T r(xωx 1 A) = T r(ω A) = 2E RPA c Multiplication of RPA equation with (T, 1) from left and X 1 from right ( ) ( ) ( ) A B 1 1 (T, 1) B A YX 1 = (T, 1) YX 1 XωX 1 Expanding yields drccd Ricatti equation B + AT + TA + TBT = 0 T := YX 1 satisfies drccd amplitude equation G. Scuseria, et. al. J. Chem. Phys. 129, (2008). Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
28 Computing Linear Response Different ways to compute first order response in DFT Response functions, Casida equation (frequency-dependent) perturbation theory, Sternheimer equation real-time propagation with weak external perturbation Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
29 Sternheimer equation Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
30 Sternheimer equation Perturbed Hamiltonian and states (zero frequency) (Ĥ0 + λh1 +...)(ψ0 + λψ1 +...) = (E0 + λe1 +...)(ψ0 + λψ1 +...) Expand and keep terms to first order in λ Ĥ 0ψ 0 + λh 1ψ 0 + λh 0ψ 1 = E 0ψ 0 + λe 0ψ 1 + λe 1ψ 0 + O(λ 2 ) Use Ĥ0ψ0 = E0ψ0 (Ĥ0 E0)ψ1 = (Ĥ1 E1)ψ0, Sternheimer equation Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
31 Sternheimer equation in TDDFT (Weak) monochromatic perturbation v 1(r, t) = λr i cos(ωt) Expand time-dependent Kohn-Sham wavefunctions in powers of λ ψ m(r, t) = exp( i(ɛ (0) m + λɛ (1) m )t) { ψ m (0) (r) + 1 } 2 λ[exp(iωt)ψ(1) m (r, ω) + exp( iωt)ψ m (1) (r, ω)] Insert in time-dependent Kohn-Sham equation and keep terms up to first order in λ Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
32 Sternheimer equation in DFT Frequency-dependent response (self-consistent solution!) [Ĥ(0) ɛ j ± ω + iη] ψ (1) (r, ±ω) = Ĥ(1) (±ω)ψ (0) (r), with first-oder frequency-dependent perturbation ρ1(r, Ĥ (1) ω) (ω) = v(r) + r r d3 r + f xc(r, r, ω)ρ 1(r, ω)d 3 r and first-order density response occ. { } ρ 1(r, ±ω) = [ψ (0) (r)] ψ (1) (r, ω) + [ψ (1) (r, ω)] ψ (0) (r) m Main advantages Only occupied states need to be considered Scales as N 2, where N is the number of atoms (Non-)Linear system of equations. Can be solved with standard solvers Disadvantage Converges slowly close to a resonance Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
33 Different types of perturbations The response equations can be used for different types of perturbations Electric perturbations v(r) = r i Response contains information about polarizabilities, absorption, fluoresence, etc. Magnetic perturbations v(r) = L i Response contains e.g. NMR signals, etc. Atomic displacements v(r) = v(r) R i Response contains e.g. phonons, etc. Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
34 Outline Linear Response in DFT Response functions Casida equation Sternheimer equation Real-space representation and real-time propagation Real-space representation for wavefunctions and Hamiltonians Time-propagation schemes Optimal control of electronic motion Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
35 Real-space grids Simulation volumes: sphere, cylinder, parallelepiped Minimal mesh: spheres around atoms, filled with uniform mesh of grid points Typically zero boundary condition, absorbing boundary, optical potential Finite-difference representation ( stencils ) for the Laplacian/kinetic energy Pseudopotentials Domain-parallelization Domain A Domain B Ghost Points of Domain B Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
36 Real-space grids Example: five-point finite difference Laplacian in 2D 1 2 ψ 2m x 1 1 [ 2 2m h ψ 2m y 1 1 [ 2 2m h 2 Stencil notation for kinetic energy 1 1 2m h 2 Leads to sparse matrices ] ψ(i 1, j) + 2ψ(i, j) ψ(i + 1, j) ] ψ(i, j 1) + 2ψ(i, j) ψ(i, j + 1) ψ(i, j) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
37 Real-space grids Size of Hamiltonian matrix can easily reach 10 7 x10 7 Basic operation Ĥψ sparse matrix vector operations Sparse solvers Conjugate gradients Krylov subspace/lanczos methods Davidson or Jacobi-Davidson algorithm Multigrid methods Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
38 Real-time evolution for the time-dependent Kohn-Sham system Time-dependent Kohn-Sham equations i tϕ j(r, t) = ( 2 2 2m v S[ρ(r, t )](r, t) = v(r, t) + ρ(r, t) = N ϕ j(r, r) 2 j=1 ) + vs[ρ](r, t) ϕ j(r, t) ρ(r, t) r r d3 r + v xc[ρ(r, t )](r, t) Initial value problem ϕ j(r, t) = ϕ (0) j (r) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
39 Real-time evolution for the time-dependent Kohn-Sham system Time-dependent Kohn-Sham equations i tϕ j(r, t) = ( 2 2 2m v S[ρ(r, t )](r, t) = v(r, t) + ρ(r, t) = N ϕ j(r, r) 2 j=1 ) + vs[ρ](r, t) ϕ j(r, t) ρ(r, t) r r d3 r + v xc[ρ(r, t )](r, t) Initial value problem ϕ j(r, t) = ϕ (0) j (r) Time-evolution operator Û(t, t0) ϕ j(r, t) = Û(t, t0)ϕj(r, t0) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
40 Properties of Û(t, t 0) Û(t, t 0) is a non-linear operator The propagator is unitary Û = Û 1 In the absence of magnetic fields the propagator is time-reversal symmetric Û 1 (t, t 0) = Û(t0, t) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
41 Properties of Û(t, t 0) Û(t, t 0) is a non-linear operator The propagator is unitary Û = Û 1 In the absence of magnetic fields the propagator is time-reversal symmetric Û 1 (t, t 0) = Û(t0, t) Equation of motion for the propagator i t Û(t, t 0) = Ĥ(t)Û(t, t0), Û(t 0, t 0 ) = ˆ1 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
42 Properties of Û(t, t 0) Û(t, t 0) is a non-linear operator The propagator is unitary Û = Û 1 In the absence of magnetic fields the propagator is time-reversal symmetric Û 1 (t, t 0) = Û(t0, t) Equation of motion for the propagator i t Û(t, t 0) = Ĥ(t)Û(t, t0), Û(t 0, t 0 ) = ˆ1 Representation in integral form Û(t, t 0) = ˆ1 i t t 0 dτĥ(τ)û(τ, t0) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
43 Properties of Û(t, t 0) Û(t, t 0) is a non-linear operator The propagator is unitary Û = Û 1 In the absence of magnetic fields the propagator is time-reversal symmetric Û 1 (t, t 0) = Û(t0, t) Equation of motion for the propagator i t Û(t, t 0) = Ĥ(t)Û(t, t0), Û(t 0, t 0 ) = ˆ1 Representation in integral form Û(t, t 0) = ˆ1 i t t 0 dτĥ(τ)û(τ, t0) Iterated solution of integral equation - time-ordered exponential ( i) n t t t Û(t, t 0) = dt 1 dt 2... dt n ˆT [ Ĥ(t 1)Ĥ(t2)... Ĥ(tn)] n! t 0 t 0 t 0 n=0 = ˆT exp( i t t 0 dτĥ(τ)) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
44 Real-time evolution - Short-time propagation Group property of exact propagator Û(t 1, t 2) = Û(t1, t3)û(t3, t2) Split propagation step in small short-time propagation intervals Û(t, t 0) = N 1 j=1 Û(t j, t j + t j) Why is this a good idea? If we want to resolve frequencies up to ω max, the time-step should be no larger than 1/ω max The time-dependence of the Hamiltonian is small over a short-time interval The norm of the time-ordered exponential is proportional to t. Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
45 Real-time evolution - Magnus expansion Time-ordered evolution operator Û(t, t 0) = ( i) n n=0 n! = ˆT exp( i t dt 1 t 0 t t t 0 dτĥ(τ)) t 0 dt 2... t dt n ˆT [ Ĥ(t 1)Ĥ(t2)... Ĥ(tn)] t 0 Magnus expansion Û(t + t, t) = exp (ˆΩ1 + ˆΩ 2 + ˆΩ ) 3 + Magnus operators ˆΩ 1 = i ˆΩ 2 = t+ t t t+ t τ1. t t Ĥ(τ)dτ [Ĥ(τ1), Ĥ(τ2)]dτ2dτ1 Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
46 Real-time evolution - Magnus expansion Second-order Magnus propagator - Exponential midpoint rule ) Û (2) (t + t, t) = exp (ˆΩ1 + O( t 3 ) ˆΩ 1 = iĥ(t + t/2) + O( t3 ). Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
47 Real-time evolution - Magnus expansion Second-order Magnus propagator - Exponential midpoint rule ) Û (2) (t + t, t) = exp (ˆΩ1 + O( t 3 ) ˆΩ 1 = iĥ(t + t/2) + O( t3 ). Fourth-order Magnus propagator Û (4) (t + t, t) = exp ˆΩ 1 = (ˆΩ1 + Ω 2 ) + O( t 5 ) i(ĥ(τ1) + Ĥ(τ2)) t 2 + O( t5 ). 3 t 2 ˆΩ 2 = i[ĥ(τ1), Ĥ(τ2)] + O( t 5 ). 12 τ 1,2 = t + ( ± 6 ) t Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
48 Real-time evolution - Crank-Nicholson/Cayley propagator Padé approximation of exponential, e.g. lowest order (Crank-Nicholson) 1 iĥ t/2 exp( iĥ t) 1 + iĥ t/2 Need only action of operator on a state vector Ψ(t + t) = 1 iĥ t/2 Ψ(t) 1 + iĥ t/2 (Non-)Linear system of equations at each time-step (1 + iĥ t/2) Ψ(t + t) = (1 iĥ t/2) Ψ(t) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
49 Real-time evolution - Operator splitting methods Typically, the Hamiltonian has the form Ĥ = ˆT + ˆV ˆT is diagonal in momentum space, ˆV in position space Baker-Campbell-Hausdorff relation eâe ˆB = exp(â + ˆB [Â, ˆB] +...) Split-Operator exp( i t( ˆT + ˆV )) exp( i t ˆT /2) exp( i t ˆV ) exp( i t ˆT /2) Use FFT to switch between momentum space and real-space. Higher-order splittings possible, but require more FFTs Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
50 Real-time evolution - Enforced time reversal symmetry Enforced time-reversal symmetry exp(+i t t Ĥ(t + t)) Ψ(t + t) = exp( i 2 2 Ĥ(t)) Ψ(t) Propagator with time-reversal symmetry Û ETRS (t + t, t) = exp( i t t Ĥ(t + t)) exp( i 2 2 Ĥ(t)) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
51 Real-time evolution - Matrix exponential Û CN (t + t, t) = 1 iĥ t/2 1 + ( i tĥ(t iĥ t/2 ) Û EM (t + t, t) = exp + t/2) Û SO (t + t, t) = exp( i t ˆT /2)exp( i t ˆV )exp( i t ˆT /2) Û ETRS (t + t, t) = exp( i t t Ĥ(t + t))exp( i 2 2 Ĥ(t))... Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
52 Real-time evolution - Matrix exponential Û CN (t + t, t) = 1 iĥ t/2 1 + ( i tĥ(t iĥ t/2 ) Û EM (t + t, t) = exp + t/2) Û SO (t + t, t) = exp( i t ˆT /2)exp( i t ˆV )exp( i t ˆT /2) Û ETRS (t + t, t) = exp( i t t Ĥ(t + t))exp( i 2 2 Ĥ(t))... Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
53 Real-time evolution - Matrix exponential C. Moler and C. Van Loan, Nineteen Dubious Ways to Compute the Exponential of A Matrix, SIAM Review 20, 801 (1978) C. Moler and C. Van Loan, Nineteen Dubious Ways to Compute the Exponential of A Matrix, Twenty-Five Years Later, SIAM Review 45, 3 (2003) Task: Compute exponential of operator/matrix Taylor series Chebyshev polynomials Padé approximations Scaling and squaring Ordinary differential equation methods Matrix decomposition methods Splitting methods Task: Compute eâv for given v Taylor series Chebyshev rational approximation Lanczos-Krylov subspace projection Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
54 Real-time evolution - Movie time Proton scattering of fast proton with ethene Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
55 Octopus code Octopus: real-space, real-time TDDFT code, available under GPL Page (Parsec: real-space, real-time code using similar concepts) libxc: Exchange-Correlation library, available under LGPL (used by many codes: Abinit, APE, AtomPAW, Atomistix ToolKit, BigDFT, DP, ERKALE, GPAW, Elk, exciting, octopus, Yambo) Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
56 Optimal control theory Control of ring current in a quantum ring external potential density profile Optimal Control of Quantum Rings by Terahertz Laser Pulses, E. Räsänen, et. al, Phys. Rev. Lett. 98, (2007). Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
57 Optimal control theory Goal: find optimal laser pulse ɛ(t) that drives the system to a desired state Φ f maximize overlap functional J 1[Ψ] = Ψ(T ) Φ f 2. constrain laser intensity T J 2[ɛ] = α 0 ɛ 2 (t) dt. 0 Lagrange multiplier density to ensure evolution with TDSE J 3[Ψ, χ, ɛ] = 2 Im T 0 (i ) χ(t) t Ĥ(t) Ψ(t) dt, Find maximum of J 1[Ψ] + J 2[ɛ] + J 3[Ψ, χ, ɛ] Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
58 Optimal control theory First variation of the functional δj = δ ΨJ + δ χj + δ ɛj = 0 Control equations δ ΨJ = 0 : δ χj = 0 : ) (i t Ĥ(t) χ(t) = 0, χ(t ) = Φ f Φ f Ψ(T ) ) (i t Ĥ(t) Ψ(t) = 0, Ψ(0) = Φ i, δ ɛj = 0 : α 0 ɛ(t) = Im χ(t) ˆµ Ψ(t). Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
59 Optimal control theory Optimal laser pulse and level population Optimal Control of Quantum Rings by Terahertz Laser Pulses, E. Räsänen, et. al, Phys. Rev. Lett. 98, (2007). Heiko Appel (Fritz-Haber-Institut der MPG) Intro to real-space, linear-response, and TD methods August 12, / 44
Propagators for TDDFT
Propagators for TDDFT A. Castro 1 M. A. L. Marques 2 A. Rubio 3 1 Institut für Theoretische Physik, Fachbereich Physik der Freie Universität Berlin 2 Departamento de Física, Universidade de Coimbra, Portugal.
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 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 informationLinear-response excitations. Silvana Botti
from finite to extended systems 1 LSI, CNRS-CEA-École Polytechnique, Palaiseau, France 2 LPMCN, CNRS-Université Lyon 1, France 3 European Theoretical Spectroscopy Facility September 3, 2008 Benasque, TDDFT
More informationTime-dependent density functional theory
Time-dependent density functional theory E.K.U. Gross Max-Planck Institute of Microstructure Physics Halle (Saale) OUTLINE Phenomena to be described by TDDFT Basic theorems of TDDFT Approximate xc functionals:
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 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 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 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 European Theoretical Spectroscopy
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 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 informationTime-dependent density functional theory
Time-dependent density functional theory E.K.U. Gross Max-Planck Institute of Microstructure Physics Halle (Saale) OUTLINE Phenomena to be described by TDDFT: Matter in weak and strong laser fields Basic
More informationNumerical Simulation of Spin Dynamics
Numerical Simulation of Spin Dynamics Marie Kubinova MATH 789R: Advanced Numerical Linear Algebra Methods with Applications November 18, 2014 Introduction Discretization in time Computing the subpropagators
More information1 Time-Dependent Two-State Systems: Rabi Oscillations
Advanced kinetics Solution 7 April, 16 1 Time-Dependent Two-State Systems: Rabi Oscillations a In order to show how Ĥintt affects a bound state system in first-order time-dependent perturbation theory
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 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 informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationExact exchange (spin current) density-functional methods for magnetic properties, spin-orbit effects, and excited states
Exact exchange (spin current) density-functional methods for magnetic properties, spin-orbit effects, and excited states DFT with orbital-dependent functionals A. Görling Universität Erlangen-Nürnberg
More informationLinear response theory and TDDFT
Linear response theory and TDDFT Claudio Attaccalite http://abineel.grenoble.cnrs.fr/ CECAM Yambo School 2013 (Lausanne) Motivations: +- hν Absorption Spectroscopy Many Body Effects!!! Motivations(II):Absorption
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 informationTime Dependent Density-Functional Theory - Linear Response
University of Tennessee, Knoxville Trace: Tennessee Research and Creative Exchange Masters Theses Graduate School 12-2013 Time Dependent Density-Functional Theory - Linear Response Bryan Edman Sundahl
More informationAdaptively Compressed Polarizability Operator
Adaptively Compressed Polarizability Operator For Accelerating Large Scale ab initio Phonon Calculations Ze Xu 1 Lin Lin 2 Lexing Ying 3 1,2 UC Berkeley 3 ICME, Stanford University BASCD, December 3rd
More information2.5 Time dependent density functional theory
.5 Time dependent density functional theory The main theorems of Kohn-Sham DFT state that: 1. Every observable is a functional of the density (Hohenger-Kohn theorem).. The density can be obtained within
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 informationTime-dependent Density Functional Theory
Time-dependent Density Functional Theory Miguel A. L. Marques and E. K. U. Gross 1 Introduction Time-dependent density-functional theory (TDDFT) extends the basic ideas of ground-state density-functional
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 informationTime Evolving Block Decimation Algorithm
Time Evolving Block Decimation Algorithm Application to bosons on a lattice Jakub Zakrzewski Marian Smoluchowski Institute of Physics and Mark Kac Complex Systems Research Center, Jagiellonian University,
More informationCorrelated Light-Matter Interactions in Cavity Quantum Electrodynamics
Correlated in Cavity Quantum Electrodynamics [1] Fritz-Haber-Institut der Max-Planck-Gesellschaft, Berlin [2] NanoBio Spectroscopy group and ETSF, Universidad del País Vasco, San Sebastián, Spain DPG Frühjahrstagung
More informationTime-Dependent Density Functional Perturbation Theory New algorithms with applications to molecular spectra
Scuola Internazionale Superiore di Studi Avanzati - ma per seguir virtute e conoscenza - Time-Dependent Density Functional Perturbation Theory New algorithms with applications to molecular spectra Thesis
More informationKey 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 information2.1 Green Functions in Quantum Mechanics
Chapter 2 Green Functions and Observables 2.1 Green Functions in Quantum Mechanics We will be interested in studying the properties of the ground state of a quantum mechanical many particle system. We
More informationThe GW Approximation. Manish Jain 1. July 8, Department of Physics Indian Institute of Science Bangalore 1/36
1/36 The GW Approximation Manish Jain 1 Department of Physics Indian Institute of Science Bangalore July 8, 2014 Ground-state properties 2/36 Properties that are intrinsic to a system with all its electrons
More informationNeutral Electronic Excitations:
Neutral Electronic Excitations: a Many-body approach to the optical absorption spectra Claudio Attaccalite http://abineel.grenoble.cnrs.fr/ Second Les Houches school in computational physics: ab-initio
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 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 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 : 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 informationFaddeev Random Phase Approximation (FRPA) Application to Molecules
Faddeev Random Phase Approximation (FRPA) Application to Molecules Matthias Degroote Center for Molecular Modeling (CMM) Ghent University INT 2011 Spring Program Fermions from Cold Atoms to Neutron Stars:
More informationGW Many-Body Theory for Electronic Structure. Rex Godby
GW Many-Body Theory for Electronic Structure Rex Godby Outline Lecture 1 (Monday) Introduction to MBPT The GW approximation (non-sc and SC) Implementation of GW Spectral properties Lecture 2 (Tuesday)
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 informationExcited states and tp propagator for fermions
Excited states and tp propagator for fermions So far sp propagator gave access to E 0 Ground-state energy and all expectation values of 1-body operators E +1 Energies in +1 relative to ground state n E0
More informationNANOSCALE SCIENCE & TECHNOLOGY
. NANOSCALE SCIENCE & TECHNOLOGY V Two-Level Quantum Systems (Qubits) Lecture notes 5 5. Qubit description Quantum bit (qubit) is an elementary unit of a quantum computer. Similar to classical computers,
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 informationAll electron optimized effective potential method for solids
All electron optimized effective potential method for solids Institut für Theoretische Physik Freie Universität Berlin, Germany and Fritz Haber Institute of the Max Planck Society, Berlin, Germany. 22
More informationScattering Theory. In quantum mechanics the basic observable is the probability
Scattering Theory In quantum mechanics the basic observable is the probability P = ψ + t ψ t 2, for a transition from and initial state, ψ t, to a final state, ψ + t. Since time evolution is unitary this
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 informationOptimal Controlled Phasegates for Trapped Neutral Atoms at the Quantum Speed Limit
with Ultracold Trapped Atoms at the Quantum Speed Limit Michael Goerz May 31, 2011 with Ultracold Trapped Atoms Prologue: with Ultracold Trapped Atoms Classical Computing: 4-Bit Full Adder Inside the CPU:
More informationNon-linear optics, k p perturbation theory, and the Sternheimer equation
Non-linear optics, k p perturbation theory, and the Sternheimer equation David A. Strubbe Department of Materials Science and Engineering Massachusetts Institute of Technology, Cambridge, MA Formerly Department
More informationLecture 1: The Equilibrium Green Function Method
Lecture 1: The Equilibrium Green Function Method Mark Jarrell April 27, 2011 Contents 1 Why Green functions? 2 2 Different types of Green functions 4 2.1 Retarded, advanced, time ordered and Matsubara
More informationQuantum Physics III (8.06) Spring 2016 Assignment 3
Quantum Physics III (8.6) Spring 6 Assignment 3 Readings Griffiths Chapter 9 on time-dependent perturbation theory Shankar Chapter 8 Cohen-Tannoudji, Chapter XIII. Problem Set 3. Semi-classical approximation
More informationTime Dependent Density Functional Theory. Francesco Sottile
An Introduction Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France European Theoretical Spectroscopy Facility (ETSF) Belfast, 29 Jun 2007 Outline 1 Introduction: why TD-DFT? 2 (Just)
More informationElectronic Structure: Density Functional Theory
Electronic Structure: Density Functional Theory S. Kurth, M. A. L. Marques, and E. K. U. Gross Institut für Theoretische Physik, Freie Universität Berlin, Arnimallee 14, D-14195 Berlin, Germany (Dated:
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 informationIntroduction 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 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 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 informationThe 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 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 information1 Mathematical preliminaries
1 Mathematical preliminaries The mathematical language of quantum mechanics is that of vector spaces and linear algebra. In this preliminary section, we will collect the various definitions and mathematical
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
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 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 informationPseudopotential 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 informationElectron Correlation - Methods beyond Hartree-Fock
Electron Correlation - Methods beyond Hartree-Fock how to approach chemical accuracy Alexander A. Auer Max-Planck-Institute for Chemical Energy Conversion, Mülheim September 4, 2014 MMER Summerschool 2014
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 informationBlock 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 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 informationPhys460.nb Back to our example. on the same quantum state. i.e., if we have initial condition (5.241) ψ(t = 0) = χ n (t = 0)
Phys46.nb 89 on the same quantum state. i.e., if we have initial condition ψ(t ) χ n (t ) (5.41) then at later time ψ(t) e i ϕ(t) χ n (t) (5.4) This phase ϕ contains two parts ϕ(t) - E n(t) t + ϕ B (t)
More informationMany-Body Perturbation Theory. Lucia Reining, Fabien Bruneval
, Fabien Bruneval Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France European Theoretical Spectroscopy Facility (ETSF) Belfast, 27.6.2007 Outline 1 Reminder 2 Perturbation Theory
More informationKey concepts in Density Functional Theory (II) Silvana Botti
Kohn-Sham scheme, band structure and optical spectra European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address:
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 informationEfficient Modeling Techniques for Time- Dependent Quantum System with Applications to Carbon Nanotubes
University of Massachusetts Amherst ScholarWorks@UMass Amherst Masters Theses 1911 - February 214 21 Efficient Modeling Techniques for Time- Dependent Quantum System with Applications to Carbon Nanotubes
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 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 informationA Brief Introduction to Linear Response Theory with Examples in Electromagnetic Response
A Brief Introduction to Linear Response Theory with Examples in Electromagnetic Response Robert Van Wesep May 3, 2013 In order to gain information about any physical system, it is necessary to probe the
More informationNumerical Integration of the Wavefunction in FEM
Numerical Integration of the Wavefunction in FEM Ian Korey Eaves ike26@drexel.edu December 14, 2013 Abstract Although numerous techniques for calculating stationary states of the schrödinger equation are
More informationIntroduction to Path Integral Monte Carlo. Part I.
Introduction to Path Integral Monte Carlo. Part I. Alexey Filinov, Jens Böning, Michael Bonitz Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany
More informationComputational Photochemistry
Computational Photochemistry Filipp Furche Department of Chemistry February 09, 2015 Turbomole Filipp Furche (UC Irvine) Computational Photochemistry 02/09/2015 1 / 27 Acknowledgments Jefferson Bates (Temple)
More informationTime Dependent Density Functional Theory. Francesco Sottile
Applications, limitations and... new frontiers Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France European Theoretical Spectroscopy Facility (ETSF) Vienna, 19 January 2007 /55 Outline
More informationApproximation Methods in QM
Chapter 3 Approximation Methods in QM Contents 3.1 Time independent PT (nondegenerate)............... 5 3. Degenerate perturbation theory (PT)................. 59 3.3 Time dependent PT and Fermi s golden
More informationOptimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics
Optimized Effective Potential method for non-collinear Spin-DFT: view to spin-dynamics Sangeeta Sharma 1,2, J. K. Dewhurst 3, C. Ambrosch-Draxl 4, S. Pittalis 2, S. Kurth 2, N. Helbig 2, S. Shallcross
More informationUltrafast laser-driven spin dynamics in antiferromagnetic and ferromagnetic solids: Ab-initio studies with real-time TDDFT
Ultrafast laser-driven spin dynamics in antiferromagnetic and ferromagnetic solids: Ab-initio studies with real-time TDDFT E.K.U. Gross Max-Planck Institute of Microstructure Physics Halle (Saale) First
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 informationTDDFT for excitation energies
TDDFT for ultrafast electronic dynamics Ivano Tavernelli and Basile Curchod Laboratory of Computational Chemistry and Biochemistry, EPFL, Lausanne TDDFT SCHOOL BENASQUE 2014 Outline 1 Overview 2 TDDFT
More informationQuantum Light-Matter Interactions
Quantum Light-Matter Interactions QIC 895: Theory of Quantum Optics David Layden June 8, 2015 Outline Background Review Jaynes-Cummings Model Vacuum Rabi Oscillations, Collapse & Revival Spontaneous Emission
More informationCP2K: the gaussian plane wave (GPW) method
CP2K: the gaussian plane wave (GPW) method Basis sets and Kohn-Sham energy calculation R. Vuilleumier Département de chimie Ecole normale supérieure Paris Tutorial CPMD-CP2K CPMD and CP2K CPMD CP2K http://www.cpmd.org
More informationThe Overhauser Instability
The Overhauser Instability Zoltán Radnai and Richard Needs TCM Group ESDG Talk 14th February 2007 Typeset by FoilTEX Introduction Hartree-Fock theory and Homogeneous Electron Gas Noncollinear spins and
More informationElectronic 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 informationTime-dependent density functional perturbation theory
Time-dependent density functional perturbation theory Iurii Timrov and Tommaso Gorni SISSA Scuola Internazionale Superiore di Studi Avanzati, Trieste, Italy Advanced Quantum ESPRESSO developers' meeting:
More informationQuantum Mechanics Solutions. λ i λ j v j v j v i v i.
Quantum Mechanics Solutions 1. (a) If H has an orthonormal basis consisting of the eigenvectors { v i } of A with eigenvalues λ i C, then A can be written in terms of its spectral decomposition as A =
More informationRandom-phase approximation correlation methods for molecules and solids
Random-phase approximation correlation methods for molecules and solids Andreas Goerling, Andreas Hesselmann To cite this version: Andreas Goerling, Andreas Hesselmann. Random-phase approximation correlation
More informationTime-dependent linear response theory: exact and approximate formulations. Emmanuel Fromager
June 215 ISTPC 215 summer school, Aussois, France Page 1 Time-dependent linear response theory: exact and approximate formulations Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie
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 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 informationErrata 1. p. 5 The third line from the end should read one of the four rows... not one of the three rows.
Errata 1 Front inside cover: e 2 /4πɛ 0 should be 1.44 10 7 ev-cm. h/e 2 should be 25800 Ω. p. 5 The third line from the end should read one of the four rows... not one of the three rows. p. 8 The eigenstate
More informationMP463 QUANTUM MECHANICS
MP463 QUANTUM MECHANICS Introduction Quantum theory of angular momentum Quantum theory of a particle in a central potential - Hydrogen atom - Three-dimensional isotropic harmonic oscillator (a model of
More informationPhys 622 Problems Chapter 5
1 Phys 622 Problems Chapter 5 Problem 1 The correct basis set of perturbation theory Consider the relativistic correction to the electron-nucleus interaction H LS = α L S, also known as the spin-orbit
More informationExcitation 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 information2 The Density Operator
In this chapter we introduce the density operator, which provides an alternative way to describe the state of a quantum mechanical system. So far we have only dealt with situations where the state of a
More informationThe quantum state as a vector
The quantum state as a vector February 6, 27 Wave mechanics In our review of the development of wave mechanics, we have established several basic properties of the quantum description of nature:. A particle
More information