time-dependent electron localization function co-workers: tobias burnus miguel marques alberto...
TRANSCRIPT
![Page 1: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/1.jpg)
Time-Dependent Electron Localization FunctionTime-Dependent Electron Localization Function
Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)
![Page 2: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/2.jpg)
Electron dynamics happens on the femto-second time scale.To probe it we need atto-second pulses.
Questions:
• How much time does it take to break a bond in a laser field?
• How long takes an electronic transition from one state to another?
• In a molecular junction, how much time does it take until a steady-state current is reached (after switching on a bias)? Is it reached at all?
Those are questions outside the realm of linear-responsetheory. To study them we have to propagate in time the TDSE or -for larger systems- the TDKS equations.
![Page 3: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/3.jpg)
How can one give a mathematical meaning to intuitive chemical concepts such as
• Single, double, triple bonds• Lone pairs
Note: • Density (r) is not useful!• Orbitals are ambiguous (w.r.t. unitary
transformations)
Electron Localization FunctionElectron Localization Function
![Page 4: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/4.jpg)
D r, r'P r, r' :
r
= conditional probability of finding an electron with spin at if we know with certainty that there is an electron with the same spin at .
= diagonal of two-body density matrix
3 4 N
23 33 N 3 3 N N
...
D r, r' d r ... d r r , r' , r ..., r
= probability of finding an electron with spin at and another electron with the same spin at .
r
r'
r'
r
![Page 5: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/5.jpg)
Coordinate transformation
r'
r
sIf we know there is an electron with spin at , then
is the (conditional) probability of finding another electron at , where is measured from the reference point .
2
s 0
dp r,s 1p r,s p r,0 s C r s
ds 3
0 0
Expand in a Taylor series:
The first two terms vanish.
Spherical average 2
0 0
1p r, s sin d d P r, s , ,
4
r
P ,r r
s
r
s
s
If we know there is an electron with spin at , then is the conditional probability of finding another electron at the distance s from .
p ,r s
r
r
![Page 6: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/6.jpg)
is a measure of electron localization.
small means strong localization at
σC r
Why? , being the s2-coefficient, gives the probability of
finding a second like-spin electron very near the reference
electron. If this probability very near the reference electron is
low then this reference electron must be very localized.
σC r
σC r
r
![Page 7: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/7.jpg)
C is always ≥ 0 (because p is a probability) and is not bounded from above.
Define as a useful visualization of localization(A.D. Becke, K.E. Edgecombe, JCP 92, 5397 (1990))
1 ELF 0 Advantage: ELF is dimensionless and
where
is the kinetic energy density of the uniform gas.
1
1ELF 2
rC
rCuni
rr6
5
3rC uni35322uni
σC r
![Page 8: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/8.jpg)
ELF
A. Savin, R. Nesper, S. Wengert, and T. F. Fässler, Angew. Chem. Int. Ed. 36, 1808 (1997)
![Page 9: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/9.jpg)
12-electron 2D quantum dot with four minima
Density ELF
E. Räsänen, A. Castro and E.K.U. Gross, Phys. Rev. B 77, 115108 (2008).
![Page 10: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/10.jpg)
For a determinantal wave function one obtains
2N
2deti
i 1
r1C r r
4 r
in the static case (i.e. for real-valued orbitals):
(A.D. Becke, K.E. Edgecombe, JCP 92, 5397 (1990))
in the time-dependent case:
T. Burnus, M. Marques, E.K.U.G., PRA (Rapid Comm) 71, 010501 (2005)
2N
2deti
i 1
r, t1C r, t r, t
4 r, t
2
j r, t r, t
![Page 11: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/11.jpg)
Acetylene in a strong laser field (ħω = 17.15 eV, I = 1.21014 W/cm2) [Snapshots of TDELF]
![Page 12: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/12.jpg)
Scattering of a high-energy proton from ethylene (Ekin(proton) = 2 keV) [Snapshots of TDELF]
![Page 13: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/13.jpg)
INFORMATION ACCESSIBLE THROUGH TDELF
How long does it take to break a bond ina laser field?
Which bond breaks first, which second, etc, in a collision process?
Are there intermediary (short-lived) bonds formed during a collision, which are not present any more in the collision products ?
![Page 14: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/14.jpg)
TDELF movies produced from TD Kohn-Sham equations
'rr
t'r'rdrtvrt't'rv 3
KSvxc[ (r’t’)](r t)
rt rtvm2
rtt
i jKS
22
j
propagated numerically on real-space grid using octopus code
octopus: a tool for the application of time-dependent density functional theory, A. Castro, M.A.L. Marques, H. Appel, M. Oliveira, C.A. Rozzi, X. Andrade, F. Lorenzen, E.K.U.G., A. Rubio, Physica Status Solidi 243, 2465 (2006).
![Page 15: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/15.jpg)
+ +
-5 Å +5 Å
0 Å
xy
R
+
– –(1) (2)
MODEL
Nuclei (1) and (2) are heavy: Their positions are fixed
S. Shin, H. Metiu, JCP 102, 9285 (1995), JPC 100, 7867 (1996)
![Page 16: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/16.jpg)
Anti-parallel spins Parallel spins
M. Erdmann, E.K.U.G., V. Engel, JCP 121, 9666 (2004)
![Page 17: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/17.jpg)
Parallel spins
M. Erdmann, E.K.U.G., V. Engel, JCP 121, 9666 (2004)
![Page 18: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/18.jpg)
Anti-parallel spins
TD-ELF is a measureof non-adiabaticity
![Page 19: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/19.jpg)
Most commonly used approximation for trρvxc
Adiabatic Approximation
adiab approxxc xc,stat
n ( r t )v r t : v n
homstat,xcv = xc potential of static homogeneous e-gas
How restrictive is the adiabatic approximation, i.e. the neglect of memory in the functional vxc[ρ(r’,t’)](r,t) ?Can we assess the quality of the exact adiabatic approximation?
e.g. ALDA homxc xc,statv r t : v r t
![Page 20: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/20.jpg)
1D MODEL
1D MODEL
Restrict motion of electrons and nuclei to 1D (along polarization axis of laser)
Replace in Hamiltonian all 3D Coulomb interactions by soft 1D interactions (Eberly et al)
22222
11
zzyx = constant
![Page 21: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/21.jpg)
Two goals of 1D calculations
1. Qualitative understanding of physical processes, such as double ionization
of He
2. Exact reference to test approximate xc functionals of time-dependent density functional theory
![Page 22: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/22.jpg)
How can we assess the quality of the adiabatic approximation?
Solve 1D model for He atom in strong laser fields (numerically) exactly. This yields exact TD density ρ(r,t).
Inversion of one-particle TDSE yields exact TDKS potential
Inversion of one-particle ground-state SE yields exact static KS potential that gives (for each separate t) ρ(r,t) as a ground-state density. This is the exact adiabatic approximation of the TDKS potential.
![Page 23: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/23.jpg)
Solid line: exact
Dashed line: exact adiabatic
E(t) ramped over 27 a.u. (0.65 fs) to the value E=0.14 a.u. and then kept constant
t = 0 t = 21.5 a.u. t = 43 a.u.
M. Thiele, E.K.U.G., S. Kuemmel, Phys. Rev. Lett. 100, 153004 (2008)
![Page 24: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/24.jpg)
M. Thiele, E.K.U.G., S. Kuemmel, Phys. Rev. Lett. 100, 153004 (2008)
4-cycle pulse with λ = 780 nm, I1= 4x1014W/cm2, I2=7x1014W/cm2
Solid line: exact
Dashed line: exact adiabatic
![Page 25: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/25.jpg)
PRIZE QUESTION No 3
For which kind of processes would you expect that the (exact) adiabatic approximation does not work?
![Page 26: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/26.jpg)
By virtue of time-dependent 1-1 correspondence, ALL observables are functionals of the TD density
some observables are easily expressed in terms of the density (no approximations involved)
e.g. TD dipole moment
HHG spectrum obtained from
Other observables are more difficult to express in terms of the density (involving further approximation)
e.g. ionization yields
2
ωd
rzdr, tρd(t) 3
![Page 27: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/27.jpg)
Calculation of ionization yields (for He)
divide |R3 in: a large “analyzing volume” A (where (r t) is actually calculated
AB
and its complement B = |R3 \ A
normalization of many-body wave function
B B
2
B A
2
A
2
2123
A
13 2t r r rdrd1
p(0)(t) p(+1)(t) p(+2)(t)
pair correlation function trtr
trr2:trrg
21
2
2121
M. Petersilka and E.K.U. Gross, Laser Physics 9, 105 (1999).
![Page 28: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/28.jpg)
A A
212123
13
A
32
A A
212123
13
A
31
A A
212123
130
trrρgtrρtrρrdrd2
1+trrρd1tp
trrρgtrρtrρrdrd-trrρdtp
trrρgtrρtrρrdrd2
1tp
![Page 29: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/29.jpg)
x-only limit for g[](r1,r2,t);
2
1t,r,rg 21onlyx
resulting ionization probabilities (mean-field expressions:
P0(t) = N1s(t)2
P+1(t) = 2N1s(t) (1- N1s(t))
P+2(t) = (1- N1s(t))2
where:
N1s(t) := d3r (r, t) = d3r| 1s (r, t)|2
A A
12
two-electron-system:
![Page 30: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/30.jpg)
Correlation Contributions
t,r,rgt,r,rg 21onlyx21 + gc[](r1,r2,t)
exactifies the mean-field expressions:
P0(t) = N1s(t)2 + K(t)
P+1(t) = 2N1s(t) (1- N1s(t)) - 2K(t)
P+2(t) = (1- N1s(t))2 + K(t)
correlation correction:
K(t) := d3r1 d3r2 (r1, t) (r2, t) gc [] (r1, r2, t)
A A
12
![Page 31: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/31.jpg)
The calculation involves two approximate functionals:
1. The xc potential vxc[](r t)
2. The pair correlation function g[](r1r2 t)
Which approximation is more critical?
![Page 32: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/32.jpg)
1D Helium atom (with soft Coulomb interaction) (Lappas, van Leeuwen, J. Phys. B 31, L249 (1998)
P(He+) exact
P(He++) exact
P(He+) with exact density and g=1/2
P(He++) with exact density and g=1/2
![Page 33: Time-Dependent Electron Localization Function Co-workers: Tobias Burnus Miguel Marques Alberto Castro Esa Räsänen Volker Engel (Würzburg)](https://reader035.vdocument.in/reader035/viewer/2022062720/56649f115503460f94c23e98/html5/thumbnails/33.jpg)