str. 1 TMF045 letní semestr 2006 TMF045 letní semestr 2006 V V Časová propagace vlnové funkce na mřížce IV. - metoda (t,t’) pro časově závislý Hamiltonián - odchozí okrajová podmínka: metoda komplexního škálování a neskalární absorpční potenciál (Lekce V)
Časová propagace vlnové funkce na mřížce IV. - metoda (t,t ’) pro časově závislý Hamiltonián - odchozí okrajová podmínka: metoda komplexního škálování a neskalární absorpční potenciál. (Lekce V ). Metoda (t,t ’). řešení časově závislých Hamiltoniánů: - PowerPoint PPT Presentation
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
str.
1
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Časová propagace vlnové funkce na mřížce IV.
- metoda (t,t’)pro časově závislý Hamiltonián
- odchozí okrajová podmínka: metoda komplexního škálování a neskalární absorpční potenciál
(Lekce V)
str.
2
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Metoda (t,t’)
• řešení časově závislých Hamiltoniánů:
– propagace po velmi krátkých časových krocích, na nichž lze Hamiltonián považovat přibližně za konstantní např. rozděleným propagátorem
– Adams-Moulton prediktor-korektor,…– metoda (t,t’) – velmi efektivní
• princip (t,t’) :– zavedeme novou funkci s dvěma
časovými souřadnicemi. Podmínka: nová funkce se na diagonále rovná propagované funkci.
ˆ
ˆ 0
ˆ ?
ti H t t
tt u t
u t
, , ,t t t
x t t x t
str.
3
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Metoda (t,t’)
– časová závislost na diagonále musí splňovat časovou Schrödingerovu rovnici:
– derivujeme podle času
– derivace diagonály – pomocí substituce
– dosazení ze Schrödingerovy rovnice:
, , ,t t t
x t t x tt t
0
,2
, , ,2 2
,t t t
t tt t t
x t t x tt t
x tt t t
ˆ , ,t t t
i i H x t x tt t
str.
4
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Metoda (t,t’)
– různé zápisy pravé strany (která musí platit pro diagonálu) pomocí dvou rozměrů:
– omezíme se na to, že Hamiltonián závisí na čase jen jako skalární operátor, takže lze napsat také např.
– rovnice pro dva rozměry, kde diagonála splňuje časově závislou Schr. r.
0
ˆ ˆ, , , , ,2 2
H x t x t H x t x t t
0
ˆ , ,
ˆ , , ,2 2 2
ˆ , , ,t t t
H x t x t
H x t x t t
H x t x t t
ˆ , , ,i i H x t x t tt t
str.
5
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Metoda (t,t’)
– řešení rovnice pro dvě časové souřadnice:
– proměnná t’ slouží jako nová souřadnice
– Hamiltonián je zobecněný operátor Floquetova typu
– propagujeme v t“, Floquet Hamiltonián nezávisí na t“, čili můžeme evoluční operátor napsat v obvyklém tvaru
– počáteční stav může být závislý na t’ libovolně
ˆ , , ,i H x t i x t tt t
ˆ ˆ ,FH H x t it
t t
str.
6
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Metoda (t,t’)
– shrnutí rovnic (t,t’) :
ˆ ,
, ;0 ;0
ˆ , , ;
ˆ, ; , ; , ;0
ˆ
; , ;
F
F
iH x t t
t t
x t x
i H x t x t tt
x t t U x t t x t
U t e
x t x t t
str.
7
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Elektromagnetické pole
• náš případ bude Hamiltonián popisující účinek elektromagnetického pole na molekulu
• definice Hamiltoniánu a „gauge invariance“
– Lorentzova síla na náboj q dána el. polem a magn. indukcí
– zavedení elektromagnetického vektorového a skalárního potenciálu:
– klasický Hamiltonián
1 dq
c dt
rF E B
,1, ,
, ,
tt t
c t
t t
A rE r r
B r A r
2
1, , , ,
2
qH t t q t
c
r p p A r r
str.
8
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Elektromagnetické pole
– volnost v definici A a phi – můžeme měnit definici, aniž by se změnila Lorentzova síla takto:
– Hamiltonián není invariantní vůči těmto transformacím (gauges)
– po kvantizaci získáme různé formy Hamiltoniánu pro různé „gauges“
– jejich řešení (vlnové funkce) se navzájem liší ve fázi
– měřitelné veličiny se neliší
, , ,
1, , ,
t t t
t t tc t
A r A r r
r r r
,
, ,i q
tct e t
r
r r
str.
9
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Elektromagnetické pole
• dipólová aproximace:– zanedbáme prostorovou závislost pole– podmínka: vlnová délka mnohem delší než
rozměry molekuly
• Hamiltonián interakce molekuly s elektromagnetickým polem v dipólové aproximaci některých známých „gauges“
– „length gauge“
– „reduced momentum gauge“
0
ˆ ˆ ˆd te
H Hc dt
A
r
0H
, t tA r A
… Hamiltonián molekuly bez pole
0ˆ ˆ ˆ
eH H t
c p A
str.
10
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• po excitaci adiabatickým pulsem (tj. intensita narůstá pomalu ve srovnání s kmitočtem sin(ωt) ) se molekula postupně dostane do kvazistacionárního stavu. Toto platí i pro velmi silná pole indukující generaci vyšších harm. frekvencí řádu TW/cm2.
• Definice el.mag. pole
• Hamiltonián:– length gauge
– red. mom. gauge
• obě varianty reprez. periodickým Hamiltoniánem
1 10 0sin sinA t c t t A t c t
0 0ˆ ˆ ˆ cosH H xe t
00
ˆ ˆ ˆ sine
H H p t
str.
11
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• kvazistacionární stavy periodického Ham.
• kde je periodické řešení Floquetova operátoru HF (Floq. stav)
; ,iEt
E Ex t e x t
,E x t
ˆ ˆ ,
ˆ , ,
F
F E E
H H x t it
H x t E x t
Příklad:
Ukažte, že Floq. stavy pro případ, že nemáme žádné pole, jsou dány vlastními stavy molekuly.
(návod: separace proměnných x a t)
str.
12
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• Řešení Floq. op. ve Fourierově bázi:
• Floq. matice
,1
, in tk k n
k n
x t x e c
n=-1
n=0
n=1
n’ -1 0 1
n -1 0 1
2T
ˆin t in tFe H e ψ ψ
vlastní stavymolekuly
str.
13
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• využijeme, že HF je součet H0, dipólu a d/dt
• část H0
• část dipólu– length gauge
0ˆ ˆ cosFH H x x t i
t
0 , 0 , ,ˆ ˆin t in t
n n n n k k ke H e H E ψ ψ ψ ψ
0
0
0, 1 1,
ˆ cos
ˆ ˆ2
ˆ2
in t in t
in t i t in t in t i t in t
n n n n
e x e t e
ex e e e x e e e
ex
ψ ψ
ψ ψ ψ ψ
ψ ψ
1cos
2i t i tt e e
dipólová matice
str.
14
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
– reduced momentum gauge
• část d/dt
• pro M funkcí psi(x) a N frekvenčních kanálů získáme MxN kvazi-energií, ale jen M je netriviálních, ostatní jsou posunuté o (n hbar ω) oproti centrální Brillouinově zóně
0
0
01, , 1
ˆ sin
ˆ ˆ2
ˆ2
in t in t
in t i t in t in t i t in t
n n n n
ep e t e
iep e e e p e e e
iep
ψ ψ
ψ ψ ψ ψ
ψ ψ
sin2
i t i tit e e
, ,in t in t
n n k ke i e nt
ψ ψ
str.
15
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• řešení: – přímá diagonalizace Floq. matice je nevýhodná
(matice je příliš velká, přičemž se ve výsledku celá informace zbytečně opakuje N-krát )
– Floquetovy stavy pro jeden fixní čas jsou také řešením evolučního operátoru pro jeden optický cyklus: je možné toho využít tak, že se konstruuje příslušný evoluční operátor v bázi stavů psi(x) a diagonalizací se získá M Floq. stavů v daném fixním čase. Jejich další propagací se dopočítá i časová závislost.
– důkaz :
ˆ ; ;
ˆ , ,
ˆ , , . . .
E E
i iEt E t T
E E
iET
E E
u T x t x t T
u T e x t e x t
u T x t e x t c b d
str.
16
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• evoluční operátor v bázi psi(x)
• je potřeba propagovat funkce psi(x) s časově závislým Hamiltoniánem – velmi výhodně lze použít metodu (t,t’), viz. str. 6
– 1. počáteční funkce je konstantní v t´:
– 2. propagujeme v t od 0 do T:
– 3. vyjádříme PSI v bázi pro x a t´:
, ˆk k k ku u T
, ;0 ;0k kx t x
ˆ, ; , , , ;0k kx t T U x t T x t
,1
, ; in tk k n k k
k n
x t T e x c
str.
17
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
– 4. kde koeficienty jsou dány takto (ortogonalita báze)
– 5. dosadíme za ket- (viz výše)
– 6. koeficienty můžeme zapsat jako nultý sloupec velké matice (srovnej Floquetova matice)
– 7. z rovnice 3 získáme matici evolučního operátoru pro jednu periodu pomocí koeficientů c:
, , ;in tn k k k kc e x x t T
,
ˆ , , , ;0
ˆ , , ;0
in tn k k k k
in tk k
c e x U x t T x t
e x U x t T x
, ,0
,ˆ , ,
n k k n
in t in tn k n k k k
k kc c kde
c e x U x t T e x
str.
18
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
– 8. výpočet koeficientů c pro dosazení do rov. 7, viz. definice rov. 6. Definice zobezněného evolučního operátoru viz str. 6. Nejčastěji se používá rozvoj zobecněného evolučního operátoru do Taylorovy řady, možno použít jakékoli propagační metody pro časově nezávislé Hamiltoniány.
,
,
,
,
01
0
,01
ˆ ;
ˆ
,k k t T
in Tk n
k k
in Tnk
k k
k kk n
in Tk k n k k
k n
k
k
k
n
u u T
u T x t T
e x c
eu
e c
c
str.
19
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
– 9. z toho vyplývá pro koeficienty:
– 10. výpočet elementů <Hm>
– aproximace konečné báze – níže uvedené platí přesně jen pro nekonečnou bázi
ˆ
0
1ˆ ˆ, ,!
F
miTH m
Fm
iTU x t T e H
m
,0
1
!
ˆ
m
n k n km
in t m in tk F k
iTc
m
e x H e x
1
ˆ ?
1
ˆ ˆ
...
ˆ
ˆ
mF
mF
m m
FmF
F
H
H
H
H H
F
nk nk
nk nk
nknk nk
nk nk
k k nk
H
n n
str.
20
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
– 11. dosazení do rov. 9
– pomocí malých matic
– 12. jak vypočítat 0-tý sloupec m-té mocniny Floquetovy matice – využití symetrie této matice (viz str.12-14), takže není třeba ukládat celou Floq. matici.
,0,0
0
1
!
mm
nk knk k
m
iTc
m
FH
1
,,0 ,0
,, ,0
1
,0 ,0
1
,0 ,0
m m
n nn nn
n nn n n n
m m
n n
m
n n nn
n
n
F F F
F F
F F
F F
H H H
H H 1
H H
H H
,0 ,00
1
!
mm
n nm
iT
m
Fc H
str.
21
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• otázka počtu frekvenčních kanálů:– Kdybychom diagonalizovali F. matici,
odpovídá počet použitých Fourierových funkcí počtu zahrnutých frekvenčních kanálů. U metody (tt´) se Fourierova báze netýká časové proměnné, ale pomocné proměnné t´. Numerická zkušenost je, že počet zahrnutých frekvenčních kanálů ve Floquetově stavu v metodě (tt´) je několikrát vyšší než počet bází pro t´.
• numerické řešení rov. 7– rov. 7
, ,0in T
k k nk kn
u e c
str.
22
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• postupným sčítáním matic cn,0
• rychlejší způsob představuje využití maticového násobení ve spojení s použitím knihoven BLAS (LAPACK)– předefinujeme matice takto:
– rovnici 7 nahradíme takto:
– což lze napsat jako maticové násobení:
,0in T
nn
e
u c
c-1,0
c0,0
c1,0
,1
,0, 1
1,
k kk M k
nk kn k M k
in Tn
u u
c c
f e
1 , 1in T
k M k n k M kn
u e c
u f c
str.
23
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
– pozn. s využitím knihoven BLAS vede druhý způsob k násobnému zrychlení výpočtu ve srovnání s optimalizovaným C. Podle platformy jde o urychlení 4x (HP workstation) až 8x (PC).
Příklad:
Navrhněte, jak realizovat druhý způsob výpočtu výše pomocí Matlabového příkazu reshape.
str.
24
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• Výpočet spektra vyšších harmonických frekvencí:– genrovaný dipól obsahyje zřetelné
složky vyšších harmonických frekvencí– důvod: elektrony jsou urychlovány
polem o základní frekvenci, ale přitom naráží na kulombický potenciál jader, který je urychluje a zpomaluje oproti „tahu pole“
jádro
kulombický potenciál+ vnější proměnné el.pole pro elektron
e
str.
25
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• vztah dipólu a vyzařované intenzity
• úprava metodou integrace per partes – vztah vhodný pro výpočty v momentové
gauge
– vztah vhodný pro výpočty v length gauge
22
20
1 Tin td t
I n e dtT dt
2
20 0
1 T Tin t in td t d tin
e dt e dtT dt T dt
ˆ
d te p t
dt
22
20 0
1 T Tin t in td t n
e dt t e dtT dt T
ˆt e x t
str.
26
TMF0
45
le
tní s
emes
tr 2
006
V
TMF0
45
le
tní s
emes
tr 2
006
V
Periodický Hamiltonián
• časově závislá střední hodnota dipólu –– opět využijeme metodu (tt´)– konstruujeme propagátory pro menší
kroky (ne jeden pro celou periodu) v bázi vlastních stavů H0 (bez elmag pole) pro
– počáteční funkce je psik propagovaná do v (m-1)tau