How to make the Born-Oppenheimer approximation exact: A fresh look at potential energy surfaces and Berry phases

Transcription

1 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)

2 Process of vision Light-induced isomeriztion

3 "Triad molecule": Candidate for photovoltaic applications C.A. ozzi et al, Nature Communications 4, 1602 (2013) S.M. Falke et al, Science 344, 1001 (2014) TDDFT propagation with clamped nuclei

4 "Triad molecule": Candidate for photovoltaic applications C.A. ozzi et al, Nature Communications 4, 1602 (2013) S.M. Falke et al, Science 344, 1001 (2014) Moving nuclei

5 Hamiltonian for the complete system of N e electrons with coordinates r r and N n nuclei with coordinates r 1 N e 1 N n Ĥ Tˆ n () Ŵ nn () Tˆ e (r) Ŵ ee (r) Vˆ en (, r) with Tˆ n N n 1 2 2M Tˆ e N e i1 2 i 2m Ŵ nn 1 2 N n, Z Z Ŵ ee 1 2 N e j j, k k r j 1 r k Vˆ en N e n j1 N r 1 j Z Stationary Schrödinger equation Ĥ r, Er,

6 Hamiltonian for the complete system of N e electrons with coordinates r r and N n nuclei with coordinates r 1 N e 1 N n Ĥ Tˆ n () Ŵ nn () Tˆ e (r) Ŵ ee (r) Vˆ en (, r) with Tˆ n N n 1 2 2M Tˆ e N e i1 2 i 2m Ŵ nn 1 2 N n, Z Z Ŵ ee 1 2 N e j j, k k r j 1 r k Vˆ en N e n j1 N r 1 j Z Time-dependent Schrödinger equation i r,, t Hr, V r,, t laser r,, t t V laser N e N n r,, t rj Z E f tcos t j 1 1

7 Born-Oppenheimer approximation solve ext BO BO Tˆ (r) Ŵ (r) Vˆ (r) Vˆ (r, ) Φ r Φ r BO e ee e en for each fixed nuclear configuration. Make adiabatic ansatz for the complete molecular wave function: Ψ BO BO BO r, r χ and find best χ BO by minimizing <Ψ BO H Ψ BO > w.r.t. χ BO :

8 Born-Oppenheimer approximation solve ext BO BO Tˆ (r) Ŵ (r) Vˆ (r) Vˆ (r, ) Φ r Φ r BO e ee e en for each fixed nuclear configuration. Make adiabatic ansatz for the complete molecular wave function: Ψ BO BO BO r, r χ and find best χ BO by minimizing <Ψ BO H Ψ BO > w.r.t. χ BO :

9 Nuclear equation ext 1 BO BO Tˆ n() Ŵnn () Vˆ n () A (-i ) υ M Φ BO * BO BO BO r Φ r dr χ Eχ Tˆ n Berry connection A BO υ γ BO BO * BO Φ r(-i υ )Φ r BO C A C d dr is a geometric phase In this context, potential energy surfaces and the vector potential follow from an APPOXIMATION (the BO approximation). BO BO A

10 Nuclear equation ext 1 BO BO Tˆ n() Ŵnn () Vˆ n () A (-i ) υ M Φ BO * BO BO BO r Φ r dr χ Eχ Tˆ n Berry connection A BO υ γ BO BO * BO Φ r(-i υ )Φ r BO C A C d dr is a geometric phase In this context, potential energy surfaces and the vector potential follow from an APPOXIMATION (the BO approximation). BO BO A

11 Geometric Phases Concept of geometric phase: Discovered by S. Pancharatnam (1956) Proc. Indian Acad. Sci. A 44: In the context of quantum mechanics: Michael V. Berry (1984) Proc. oyal Society 392 (1802),

12 Whenever the Hamiltonian of a quantum system depends on a vector of parameters,, the Berry phase is defined as: where the line integral is along a closed loop, C, in parameter space. A non-vanishing value of only appears when C encircles some non-analyticity.

13 Standard representation of the full TD wave function Expand full molecular wave function in complete set of BO states: BO,J J Ψ r,,t Φ r χ,t J and insert expansion in the full Schrödinger equation standard BO non-adiabatic coupling terms from T n acting on Φ r.,j

14 Plug Born-Huang expansion in full TDSE: i, t T, t, t t k n k k k 2 BO BO,k i, j i j, t j M NAC-1 2 j 2M BO 2 BO,k, j j, t NAC-2 The dynamics is "non-adiabatic" when the NAC terms cannot be neglected

15 BO Φ1, r E BO 1 BO E Φ r 0 BO 0, Ψ BO, t,t, t BO r, χ Φ r χ Φ r , 01 1, When only few BO-PES are important, the BO expansion gives a perfectly clear picture of the dynamics

16 Example: NaI femtochemistry Na + + I - Na + I

17 Example: NaI femtochemistry Na + + I - Na + I emitted neutral Na atoms

18 Effect of tuning pump wavelength (exciting to different points on excited surface) λ pump /nm Different periods indicative of anharmonic potential T.S. ose, M.J. osker, A. Zewail, JCP 91, 7415 (1989)

19 For larger systems one would like to (one has to) treat the nuclei classically.

20 Trajectory-based quantum dynamics

21 For larger systems one would like to (one has to) treat the nuclei classically. But what s the classical force when the nuclear wave packet splits??

22 For larger systems one would like to (one has to) treat the nuclei classically. But what s the classical force when the nuclear wave packet splits??

23 For larger systems one would like to (one has to) treat the nuclei classically. But what s the classical force when the nuclear wave packet splits?? There is only one correct answer!

24 Outline Show that the factorisation Ψ r, r χ can be made exact Concept of exact PES and exact Berry phase Concept of exact and unique time-dependent PES Mixed quantum-classical treatment

25 THANKS Axel Schild Ali Abedi Federica Agostini Yasumitsu Suzuki Seung Kyu Min Neepa Maitra (Hunter College, CUNY) yan equist Nikitas Gidopoulo (Durham University, UK)

26 Theorem I The exact solutions of Ĥ r, Er, can be written in the form Ψ r, r χ 2 where dr Φ r 1 for each fixed. N.I. Gidopoulos, E.K.U. Gross, Phil. Trans.. Soc. 372, (2014)

27 Proof of Theorem I: Given the exact electron-nuclear wavefuncion Ψ r, Choose: 2 is : e dr r, with some real-valued funcion S Φ r : Ψ r, / χ Then, by construction, dr Φ 2 r 1

28 Proof of Theorem I: Given the exact electron-nuclear wavefuncion Ψ r, Choose: 2 is : e dr r, with some real-valued funcion S Φ r : Ψ r, / χ Then, by construction, dr Φ 2 r 1 Note: If we want () to be smooth, S() may be discontinuous

29 Immediate consequences of Theorem I: 1. The diagonal of the nuclear N n -body density matrix is identical with χ 2 in this sense, χ 2 2 proof: Γ dr Ψr, dr Φ 2 r χ χ 1 2 can be interpreted as a proper nuclear wavefunction.

30 Theorem II: r and satisfy the following equations: ext Eq. Tˆ Ŵ Vˆ Vˆ i A e ee e Ĥ BO en N n 1 2M N n 1 i χ A M χ N 1 ext nn n 2M 2 i A Φ r Φ r n 2 Eq. i A Ŵ Vˆ χ Eχ * where A i r r dr N.I. Gidopoulos, E.K.U. Gross, Phil. Trans.. Soc. 372, (2014)

31 Theorem II: r and satisfy the following equations: ext Eq. Tˆ Ŵ Vˆ Vˆ i A e ee e Ĥ BO en N n 1 2M N n 1 i χ A M χ N 1 ext nn n 2M 2 i A Φ r Φ r n 2 Eq. i A Ŵ Vˆ χ Eχ * where A i r r dr N.I. Gidopoulos, E.K.U. Gross, Phil. Trans.. Soc. 372, (2014) Exact PES Exact Berry potential

32 How do the exact PES look like?

33 MODEL S. Shin, H. Metiu, JCP 102, 9285 (1995), JPC 100, 7867 (1996) (1) (2) Å -5 Å +5 Å x Nuclei (1) and (2) are heavy: Their positions are fixed

34

35

36 Exact Berry connection * A dr r i r Insert: r r, / i : e * 2 A Im dr r, r, / A J / 2 with the exact nuclear current density J

37 Another way of reading this equation: 2 J {A } Conclusion: The nuclear Schrödinger equation Nn 1 ext nn n 2M 2 i A Ŵ Vˆ χ Eχ yields both the exact nuclear N-body density and the exact nucler N-body current density A. Abedi, N.T. Maitra, E.K.U. Gross, JCP 137, 22A530 (2012)

38 Question: Can the true vector potential be gauged away, i.e. is the true Berry phase zero?

39 Question: Can the true vector potential be gauged away, i.e. is the true Berry phase zero? Look at Shin-Metiu model in 2D: (1) (2) + + +

40 BO-PES of 2D Shin-Metiu model

41 BO-PES of 2D Shin-Metiu model conical intersection with Berry phase

42 Non-vanishing Berry phase results from a non-analyticity in the electronic wave function BO r as function of. Such non-analyticity is found in BO approximation.

43 Non-vanishing Berry phase results from a non-analyticity in the electronic wave function BO r as function of. Such non-analyticity is found in BO approximation. Does the exact electronic wave function show such non-analyticity as well (in 2D Shin-Metiu model)? Look at D rφ rdr as function of nuclear mass M. S.K. Min, A. Abedi, K.S. Kim, E.K.U. Gross, PL 113, (2014)

44 D() M =

45 D() M =

46 Question: Can one prove in general that the exact molecular Berry phase vanishes?

47 Question: Can one prove in general that the exact molecular Berry phase vanishes? Answer: No! There are cases where a nontrivial Berry phase appears in the exact treatment.. equist, F. Tandetzky, EKU Gross, Phys. ev. A 93, (2016).

48 Time-dependent case

49 Theorem T-I The exact solution of t i r,, t H r,, t r,, t can be written in the form r,, t r, t, t where 2 dr r, t 1, t for any fixed. A. Abedi, N.T. Maitra, E.K.U.G., PL 105, (2010) JCP 137, 22A530 (2012)

50 Eq. r, t and, t Theorem T-II satisfy the following equations Tˆ e N n Eq. Ŵ 1 M ee Vˆ ext e N n r, t Vˆ r, i A, t t Ĥ BO i χ, t χ, t en A 1 2M, t i A, tφ r i Φ r, t 2 t N n 1 2M 2 ext i A, t Ŵ Vˆ, t, tχ, t i χ, t nn n t A. Abedi, N.T. Maitra, E.K.U.G., PL 105, (2010) JCP 137, 22A530 (2012)

51 Eq. r, t and, t Theorem T-II satisfy the following equations Tˆ e N n Ŵ 1 M ee Vˆ ext e N n r, t Vˆ r, i A, t t Ĥ BO i χ, t χ, t en A 1 2M, t i A, tφ r i Φ r, t 2 t Eq. Exact Berry potential Exact TDPES N n 1 2M 2 ext i A, t Ŵ Vˆ, t, tχ, t i χ, t nn n t A. Abedi, N.T. Maitra, E.K.U.G., PL 105, (2010) JCP 137, 22A530 (2012)

52 How does the exact time-dependent PES look like?

53 Example: Nuclear wave packet going through an avoided crossing (Zewail experiment) A. Abedi, F. Agostini, Y. Suzuki, E.K.U.Gross, PL 110, (2013) F. Agostini, A. Abedi, Y. Suzuki, E.K.U. Gross, Mol. Phys. 111, 3625 (2013)

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75 New MD scheme: Perform classical limit of the nuclear equation, but retain the quantum treatment of the electronic degrees of freedom. A. Abedi, F. Agostini, E.K.U.Gross, EPL 106, (2014)

76 Eq. Theorem T-II Tˆ e N n Eq. Ŵ 1 M ee Vˆ ext e N n r, t Vˆ r, i A, t t Ĥ BO i χ, t χ, t en A 1 2M, t i A, tφ r i Φ r, t 2 t N n 1 2M 2 ext i A, t Ŵ Vˆ, t, tχ, t i χ, t nn n t

77 Eq. Theorem T-II Tˆ e N n Eq. Ŵ 1 M ee Vˆ ext e N n r, t Vˆ r, i A, t t Ĥ BO i χ, t χ, t en A 1 2M, t i A, tφ r i Φ r, t 2 t N n 1 2M 2 ext i A, t Ŵ Vˆ, t, tχ, t i χ, t nn n t

78 Shin-Metiu model

79 Propagation of classical nuclei on exact TDPES

80 Measure of decoherence: Quantum: d C 1, t 2 C2, t 2 χ, t 2 Trajectories 1 N traj I C 1 I t 2 C 2 I (t) 2

81 Algorithm implemented in:

82 The "right" electron-phonon interaction

83 electron-phonon interaction?? Hˆ M k, q c? c b b eph kq k q q k,q, k-q q + k-q -q k k

84 electron-phonon interaction?? Hˆ M k, q c? c b b eph kq k q q k,q, k-q q k M k, q + k-q -q k M k, q In a genuine ab-initio description, what is the exact coupling M k, q?

85 LITEATUE on M k, q : What everybody uses: n p V n ' p' q KS q p p g, p' pq 12 M n,n ' ' 2MN c q obert van Leeuwen, PB 69, (2004): 12 Z 1 i 0, M, 2MN c / r q r q d r1 e r, r1 ; q, e Many textbooks neglect ϵ -1 r1 0, completely (no screening). Higher-order terms (Marini, Ponce, Gonze, PB 91, (2015) (using DFPT): (2)?? eph q, q' ' k k-q- q' -q q - q' ' q' ' k, q, q' '?? H ˆ (k)c? c b b b b

86 Theorem II: r and satisfy the following equations: ext Eq. Tˆ Ŵ Vˆ Vˆ i A e ee e Ĥ BO en N n 1 2M N n 1 i χ A M χ N 1 ext nn n 2M 2 i A Φ r Φ r n 2 Eq. i A Ŵ Vˆ χ Eχ Exact phonons Expand ϵ() around equilibrium positions (to second order): Ĥph? q k bqbq q Eq. 1 2

87 Theorem II: r and satisfy the following equations: ext Eq. Tˆ Ŵ Vˆ Vˆ i A e ee e Ĥ BO en N n 1 2M N n 1 i χ A M χ Exact el-ph interaction 2 i A Φ r Φ r N 1 ext nn n 2M n 2 Eq. i A Ŵ Vˆ χ Eχ Exact phonons Expand ϵ() around equilibrium positions (to second order): Ĥph? q k bqbq q Eq. 1 2

88 ˆ q n p u VKS n ' p' Mq n p,n ' p ' g, p' pq n p,n ' p' fhxc n p,n ' p' 2MN n q Traditional term

89 Summary on exact factorisation t t t Ψ r,, r, χ, is exact A. Abedi, N.T. Maitra, E.K.U. Gross, PL 105, (2010) Exact Berry phase vanishes S.K. Min, A. Abedi, K.S. Kim, E.K.U. Gross, PL 113, (2014) TD-PES shows jumps resembling surface hopping A. Abedi, F. Agostini, Y. Suzuki, E.K.U.Gross, PL 110, (2013) mixed quantum classical algorithms S.K. Min, F Agostini, E.K.U. Gross, PL 115, (2015) correct electron-phonon interaction shows new terms (in addition to standard expression)

90 SFB 450 SFB 685 SFB 762 SPP 1145