+ All Categories
Home > Documents > Klasick á-kvantová korespondence ve fázovém prostoru

Klasick á-kvantová korespondence ve fázovém prostoru

Date post: 17-Jan-2016
Category:
Upload: tyrell
View: 25 times
Download: 3 times
Share this document with a friend
Description:
Klasick á-kvantová korespondence ve fázovém prostoru. lekce (X). Obsah:. vlnová funkce z Wignerovy distribuce WKB teorie a Wignerova reprezentace autokorelační funkce z časově závislé Wignerovy distribuce - PowerPoint PPT Presentation
45
str. 1 TMF045 letní semestr 2006 TMF045 letní semestr 2006 X X Klasická-kvantová korespondence ve fázovém prostoru lekce (X)
Transcript
Page 1: Klasick á-kvantová korespondence ve fázovém prostoru

str.

1

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Klasická-kvantová korespondence

ve fázovém prostoru

lekce (X)

Page 2: Klasick á-kvantová korespondence ve fázovém prostoru

str.

2

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Obsah:

vlnová funkce z Wignerovy distribuce WKB teorie a Wignerova reprezentace autokorelační funkce z časově závislé

Wignerovy distribuce rozdíl mezi kvantovou a klasicky

propagovanou distribucí – „nebezpečné“ krostermy, Morseho index

Wignerova metoda škálování s velikostí systému zavedení Morseho indexu pro výpočet acf

teplotní závislost WD kvantové korekce difúzní Monte Carlo ve fázovém prostoru

stacionární Schrödingerova rovnice ve fázovém prostoru

Page 3: Klasick á-kvantová korespondence ve fázovém prostoru

str.

3

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Vlnová funkce z Wignerovy distribuce

• motivace: známe (např. klasickou aprox. časového vývoje) Wignerovu distribuci, z níž chceme získat vlnovou funkci

• řešení: využijeme definice středních hodnot pomocí W.d.

• výhoda využití středních hodnot pro klasické aproximace:– při aproximativní W.d., kdy ji aproximujeme

klasickou distribucí, je chyba nejmenší u momentů nízkého řádu. Např. norma, stř. hodnota hybnosti nebo pozice považujeme za momenty nízkého řádu. Naopak stř. hodnota p10 by patřila mezi momenty poměrně vyššího řádu.

– ZDŮVODNĚNÍ: momenty vyššího řádu jsou daleko více ovlivněny přítomností rychlých oscilačních členů na W.d., které nejsou reprodukovány klasickou mechanikou.

Page 4: Klasick á-kvantová korespondence ve fázovém prostoru

str.

4

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Vlnová funkce z Wignerovy distribuce

• využití středních hodnot (opak. minulé látky a nové označení)– Wignerova distribuce funkce f

– Wignerova reprezentace operátoru A

– kvantová stopa – pro B def. jako matice hustoty stavu f z ní vyplývá vztah pro střední hodnoty

Page 5: Klasick á-kvantová korespondence ve fázovém prostoru

str.

5

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Vlnová funkce z Wignerovy distribuce

• definice amplitudy a fáze vlnové funkce pomocí středních hodnot veličin

• amplituda

– korespondující Wignerovy reprezentace obou operátorů hustoty:

– výsledná definice amplitudy vlnové funkce pomocí W.d.:

Page 6: Klasick á-kvantová korespondence ve fázovém prostoru

str.

6

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Vlnová funkce z Wignerovy distribuce

• derivace fáze:– definice fáze vlnové funkce:

– výpočet derivace fáze podle x a t (času) umožňuje využít jednoduché stř. hodnoty:

– dosadíme za derivace vlnové funkce pomocí známých definic:

Page 7: Klasick á-kvantová korespondence ve fázovém prostoru

str.

7

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Vlnová funkce z Wignerovy distribuce

– a jsou dány jako stopy operátorů:

– evaluace stopy DVOU operátorů pomocí W.d. je známa. Nyní je třeba součin tří operátorů výše rozdělit, takže vznikne delta funkce v x a součinový operátor B:

– Wignerův obraz delta funkce v x je delta funkce ve fázovém prostoru (viz výše), zatímco Wignerův obraz operátorů B je dán takto:

ˆˆ

ˆ ˆ

ˆ ˆ

ˆ ˆ ˆ ˆ ˆ

p

p

H

H V T

x x p f f x x B

B f f p

x x H f f x x B

B B B f f V T

Page 8: Klasick á-kvantová korespondence ve fázovém prostoru

str.

8

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Vlnová funkce z Wignerovy distribuce

– Bp

– BH

– důkaz viz níže

– dosazením dostaneme výrazy pro derivaci fáze vlnové funkce, které lze přímo využít pro zpětnou transformaci z Wignerovy reprezentace do vlnové funkce:

Page 9: Klasick á-kvantová korespondence ve fázovém prostoru

str.

9

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Vlnová funkce z Wignerovy distribuce

– semiklasická aproximace: nemá smysl započítávat korekce energie vyššího než druhého řádu.

– I. řád je klasická energie– II. řád je harmonická korekce na nulové kmity– klasická Liouvillova rovnice je dána kvantovým

rozvojem do druhého řádu, vyšší korekce neobsahuje, nemá tedy cenu je zahrnovat pro zpětnou transformaci.

• využití zpětné transformace:

– odvození WKB teorie z reprezentace kv. mech. ve fázovém prostoru (vysvětlení kvantování energie u vázaných stavů z dynamiky klasické distribuce) – viz níže

– jiné modelové semiklasické aplikace s využitím pro interpretaci kvantové mechaniky (viz seminář o nehermitovské mechanice)

v semiklas. aproximaci

Page 10: Klasick á-kvantová korespondence ve fázovém prostoru

str.

10

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova reprezentace součinových operátorů

• definice součinového operátoru:

• operátory, které jsou funkcí x (např. potenciál)

– definice Wignerovy reprezentace BV

– úprava kvantové amplitudy

– dosazení:

ˆˆAB A f f

ˆ ˆA V x

Page 11: Klasick á-kvantová korespondence ve fázovém prostoru

str.

11

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova reprezentace součinových operátorů

– rozvoj potenciálu do Taylorovy řady:

– dosazení:

• operátory, které jsou funkcí p (např. kinetický operátor)

– odvození je analogické. Použijeme definici W.d. pomocí hybnostní reprezentace, kde vystupuje opačná Fourierova transformace matice hustoty.

ˆ ˆA T p

0

1,

! 2

n n n

T n nn

i f d TB x p

n x dp

Page 12: Klasick á-kvantová korespondence ve fázovém prostoru

str.

12

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova reprezentace součinových operátorů

• Hermitovsky sdružený součinový operátor je zobrazen jako komplexně sdružená funkce ve fázovém prostoru, pokud A=A+

– důkaz:

– v aplikacích se často setkáváme se součty a rozdíly B a B+, což vede na použití pouze Re či Im složky B

ˆˆAB f f A

*

*

ˆ,2 2

.

ˆ2 2

ˆ2 2

,

ip xA

ip x

ip x

A

x xB x p d x e x f f A x

subs x x

x xd x e x f f A x

x xd x e f x x A f

B x p

Page 13: Klasick á-kvantová korespondence ve fázovém prostoru

str.

13

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova reprezentace součinových operátorů

• často používáme pouze reálné příp. pouze imaginární části Wig. reprezentace součinových operátorů:

2

2,4,

2

2,4,

2 2

2

1 2

1,3,

1 2

1Re ,

! 2

1Re ,

! 2

8

1Im ,

! 2

2

1Im ,

! 2

n n n n

V n nn

n n n n

T n nn

n n n n

T n nn

n

V

f d VB x p V f

n p dx

f d TB x p T f

n x dp

fT f

x

f d TB x p

n x dp

f p

x

B x pn

1,3,

1 2

3,5,

1

2 ! 2

n n n

n nn

n n n n

n nn

f d V

p dx

f dV f d V

x dx n p dx

Page 14: Klasick á-kvantová korespondence ve fázovém prostoru

str.

14

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Základní rovnice kvantové mechaniky formulované ve fázovém prostoru

• kvantová Liouvillova rovnice:

– závislost matice hustoty na čase se zobrazí jako závislost W.d. na čase:

– dosadíme za součinový operátor:

ˆ ˆ

ˆ ˆ

i H i Ht t

i it t t

H H

,, , 2 Im ,H H H

x pi B x p B x p i B x p

t

1 2 1

3,5,

,

1

! 2

n n n n

n nn

x p p dV

t x x dx

d V

n p dx

Page 15: Klasick á-kvantová korespondence ve fázovém prostoru

str.

15

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Základní rovnice kvantové mechaniky formulované ve fázovém prostoru

– srovnání s klasickou Liouvillovou rovnicí:– odvozena na základě zachování objemu klasické

distribuce ve fázovém prostoru:

– pro harmonický potenciál je stejná kvantová a klasická pohybová rovnice

– vyšší členy kvantového rozvoje jsou zodpovědné za:

– vznik oscilačních členů při propagaci Wignerovy distribuce

– tzv. hluboké kvantové jevy: tunelování a odraz nad ostrou hranou

0d dx dp

dt t x dt p dt

p dV

t x p dx

Page 16: Klasick á-kvantová korespondence ve fázovém prostoru

str.

16

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Základní rovnice kvantové mechaniky formulované ve fázovém prostoru

• kvantová difúzní rovnice:– závislost matice hustoty na reciproké teplotě:

– formulace ve fázovém prostoru– nyní matice hustoty komutuje s Hamiltoniánem, proto

můžeme psát Blochovu rovnici takto:

– převod do fázového prostoru

ˆ

ˆ

ˆ

ˆ

ˆ ˆ ˆ ˆ

n

n

E Hn n

n

n n n

E Hn n n

n

e e

kde H E

E e He H

ˆ ˆ ˆˆ ˆ2 H H

,2 , , 2Re ,H H H

x pB x p B x p B x p

Page 17: Klasick á-kvantová korespondence ve fázovém prostoru

str.

17

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Základní rovnice kvantové mechaniky formulované ve fázovém prostoru

– dosazení za součinové operátory

– srovnání s klasickou hustotou ve fázovém prostoru:

– klasická matice hustoty pro nulovou teplotu je dána delta funkcí v okolí minimální geometrie

– kvantové členy II. řádu jsou zodpovědné za nulové kmity (jsou to difúzní členy, které lze realizovat podobně jako u metody difúzní Monte Carlo)

– kvantové členy vyšších řádů jsou zodpovědné za oscilační členy při velmi nízkých teplotách (přítomné např. u dvojitého minima)

– pro vysokou teplotu přechází Wignerova distribuce na klasickou

2 2 2

2 2

2

4,6,

1

8

1

! 2

n n n n

n nn

H V xx p

d V

n p dx

He

H

Page 18: Klasick á-kvantová korespondence ve fázovém prostoru

str.

18

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Základní rovnice kvantové mechaniky formulované ve fázovém prostoru

– srovnání s metodou difúzní Monte Carlo– účelem DMC je získat základní stav u větších

systémů propagací vlnové funkce v imaginárním čase, kde vlnová funkce je reprezentována množinou náhodných bodů

– náhodné body se propagují pomocí náhodné procházky:

2 2

2

ˆ

2

i Hti

tt

Vx

i

v každém krokuse zmenšuje potenciál

v každém krokuse zvětší neurčitostdíky difúznímu členu

Page 19: Klasick á-kvantová korespondence ve fázovém prostoru

str.

19

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Základní rovnice kvantové mechaniky formulované ve fázovém prostoru

– realizace propagace v imaginárním čase pro Wignerovy distribuce (tj. pro matici hustoty)

– propagace v imaginárním čase odpovídá vývoji matice hustoty v reciproké teplotě:

– metodu DMC lze realizovat ve fázovém prostoru, kde nedostáváme jen základní stav, ale také teplotní závislost

ˆ ˆ

ˆ ˆ

H H

H H

2 2 2

2 2

1

8H V x

x p

v každém kroku sesnižuje klas. energie

v každém kroku se zvyšujeneurčitost v x a p díkydifúzním členům

Page 20: Klasick á-kvantová korespondence ve fázovém prostoru

str.

20

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Základní rovnice kvantové mechaniky formulované ve fázovém prostoru

• Stacionární Schrödingerova rovnice:

– pro jakékoli X musí platit:

ˆn n nH E

ˆ

ˆtr tr

n n n n

n n n n n

H E

H E

ˆnHB součinový operátor

, ,

, ,

, ,

n

n

H

n

H n

L dxdp x p B x p

P E dxdp x p x p

B x p E x p

Page 21: Klasick á-kvantová korespondence ve fázovém prostoru

str.

21

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Základní rovnice kvantové mechaniky formulované ve fázovém prostoru

– dosadíme

– napíšeme v operátorovém tvaru:

– využití: v principu lze počítat W.d. vlastních stavů Hamiltoniánu přímo v fázovém prostoru

– klasická aproximace vlastních stavů je dána klasickými orbity:

0

11

! 2

n n n n nn

nn n n nn

i d V d TE

n p dx x dp

0

ˆ

1ˆ 1! 2

n n n

n n n n nn

n n n nn

E

i d V d T

n dx p dp x

H

H =

ˆ ,

, n n n

H x p

H x p E

H

Page 22: Klasick á-kvantová korespondence ve fázovém prostoru

str.

22

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Základní rovnice kvantové mechaniky formulované ve fázovém prostoru

– korespondence s WKB teorií:

– pokud vezmeme v úvahu vždy jen izolované větve klasických orbitů (např. u harmonického oscilátoru pohyb jen jedním směrem), a ty převedeme do vlnové funkce (viz níže) pak získáme WKB aproximaci:

– amplituda ve fázovém prostoru je dána rychlostí ve fázovém prostoru c:

x

p

2 2 21, c x p

c

Page 23: Klasick á-kvantová korespondence ve fázovém prostoru

str.

23

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Základní rovnice kvantové mechaniky formulované ve fázovém prostoru

– amplituda vlnové funkce je dána projekcí hustoty ve fázovém prostoru do osy x:

– derivace fáze podle x je dána průměrnou hybností

– závěr: WKB teorii lze získat na základě analýzy klasické aproximace stacionárních stavů ve (Wignerově) fázovém prostoru

– při převodu klasické W.d.do vl.fce je třeba ji rozdělit do větví oddělených body obratu

0

1exp

x

WKBx

ix p x dx

p x

2

2

1cos

1 1cos

WKB

WKB

c

x cx p

Page 24: Klasick á-kvantová korespondence ve fázovém prostoru

str.

24

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Autokorelační funkce z časově závislé Wignerovy distribuce

• motivace: – předpokládaný systém:

více dimenzionální systém, např. 10-50 st. volnosti

nelineární chováním (anharmonický, případně s chaotickým chováním)

vlnová funkce obsazuje značnou část prostoru a není ji možno numericky vyjádřit pomocí vhodné báze ani na mřížce

– předpokládáme, že časově závislou Wignerovu distribuci LZE aproximovat klasicky vyvíjenou distribucí (níže se věnujeme problematice této aproximace. Předpokládáme využitelnost reprezentace Monte Carlo body ve fázovém prostoru.)

– chceme se vyhnout transformaci klasické distribuce do vlnové funkce (nevešla by se do počítače) výpočet autokorelační funkce je žádoucí provést PŘÍMO ze znalosti mnohadimenzionální distribuce.

– obecněji odvodíme výpočet pro kros-korelační funkci, později ukážeme využití kros-korelačních funkcí pro klasickou aproximaci

Page 25: Klasick á-kvantová korespondence ve fázovém prostoru

str.

25

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Autokorelační funkce z časově závislé Wignerovy distribuce

• definice kros-korelační funkce:

• autokorelační funkce je speciální případ:

• pro výpočet C z W.d. rozdělíme C na amplitudu a fázi, podobně jako u převodního vztahu pro vlnovou funkci. Tím získáme vhodnou podobu vztahů pro klasickou aproximaci.

ˆ

,i

HtC t t e

, 0 0

0

C t t

kde

2

*

, , exp ,

, tr

,2exp ,

,

iC t C t t

C t t t

C tit

C t

Page 26: Klasick á-kvantová korespondence ve fázovém prostoru

str.

26

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Autokorelační funkce z časově závislé Wignerovy distribuce

• výpočet amplitudy jako překryvu dvou Wignerových distribucí:

– dobrý souhlas mezi kvantovou acf a kvaziklasickou aproximací známý od 70-80 let, problém představoval výpočet fáze acf: „Well into the non-linear regime the classical and quantum traces agree. If only the phases can be put back!“ E. Heller 1975

2, tr

, ; , ;

C t t t

d d t

x p x p x p

Úvaha: V literatuře lze u mnoha autorů do nedávné doby pozorovat mylný předpoklad, že klasická distribuce ve fázovém prostoru neobsahuje informaci o kvantové fázi. Bylo by však možné získat propagaci amplitudy aniž by informace o fázi byla v klasické distribuci nějak přítomná?

Page 27: Klasick á-kvantová korespondence ve fázovém prostoru

str.

27

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Autokorelační funkce z časově závislé Wignerovy distribuce

• výpočet kvantové fáze kros-korelační funkce– využijeme takové definice kvantové fáze, která

vede na využití středních hodnot – takže místo fáze počítáme její časovou derivaci

– časová derivace acf

*

*

*

*

2

*

2

2, log , log ,

, ,2,

, ,

, , . .

,

Im , ,,

,

it C t C t

C t C tit

C t C t

C t C t c c

C t

C t C tt

C t

ˆ ˆˆ,

ˆ

i iHt Hti

C t e Het t

iH t

Page 28: Klasick á-kvantová korespondence ve fázovém prostoru

str.

28

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Autokorelační funkce z časově závislé Wignerovy distribuce

– dosazení do rovnice pro fázi kros-korelační funkce podle času výše:

– jmenovatel je amplituda kros-korelační funkce, kterou získáme jako překryv Wignerových distribucí – viz výše

– čitatel je dán kv. stopou tří operátorů, kterou je třeba vyjádřit jako stopu DVOU operátorů, abychom ji mohli vyjádřit pomocí Wignerovy reprezentace:

*

*

ˆ, ,

ˆtr

1 ˆIm , , Re tr

iC t C t t H t

iH t t

C t C t H t t

2

ˆRe tr,

,

H t tt

C t

Page 29: Klasick á-kvantová korespondence ve fázovém prostoru

str.

29

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Autokorelační funkce z časově závislé Wignerovy distribuce

– dosadíme za Re B:

ˆRe tr

ˆRe tr

Re , ; , ;

H

H

H t t

B t t

d d B t

x p x p x p

2 2 2

2 2

2

4,6,

1Re

8

1

! 2

H

n n n n

n nn

B H V xx p

d V

n p dx

2 2 2

2 2

ˆRe tr

, ; , ;

8

H t t

d d H t

d d V

x p x p x p

1x p x

x p

Page 30: Klasick á-kvantová korespondence ve fázovém prostoru

str.

30

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Autokorelační funkce z časově závislé Wignerovy distribuce

– v prvním přiblížení je derivace fáze kros-korelační funkce podle času dána průměrnou energií překryvu mezi distribucemi

– kvantové členy II. řádu jsou zodpovědné za posuv spekrálních čar díky nulovým kmitům (Fourierova transformace acf dává spektrum vlastních energií)

Page 31: Klasick á-kvantová korespondence ve fázovém prostoru

str.

31

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

• počáteční vlnovou funkci (často Gaussián v mnoha dimenzích) převedeme na Wignerovu distribuci a tu pak reprezentujeme jako sadu Monte Carlo bodů s jednotkovou váhou.

• počítáme časově závislou klasickou distribuci s počáteční podmínkou, která je dána počáteční Wignerovou distribucí. To uděláme tak, že propagujeme klasické trajektorie, jejichž počáteční podmínky byly určeny výše. Tím získáme nové body v čase t, které reprezentují požadovanou distribuci jako Monte Carlo body s jednotkovou váhou.

, ,k kx x p x p

, , ;k k clx t p t x p t

Page 32: Klasick á-kvantová korespondence ve fázovém prostoru

str.

32

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

• výpočet amplitudy acf pomocí MC reprezentace:

• výpočet fáze pomocí MC reprezentace:

– u klasické Wignerovy metody nemá cenu přidávat vyšší korekce než II. řádu, poněvadž již nejsou obsaženy v samotné klasické propagaci

2

1

, , ; , ;

, ;N

k kk

C t d d t

x p x p x p

x p

1

1

2 2

2 2 21 1

1

, , ;,

, ;

, ; , ;

8 , ;

N

k k k kk

N

k kk

N N

k k k kk k

N

k kk

Ht

V

x p x p

x p

1x p x x p

x p

x p

Page 33: Klasick á-kvantová korespondence ve fázovém prostoru

str.

33

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

• problém klasické propagace:

– počáteční stav: počáteční stav by měl být reprezentován pozitívně definitní Wignerovou distribucí ze dvou důvodů:

– a. kros-termy lze velmi obtížně vzorkovat malým množstvím MC bodů, protože tvoří jemnou oscilační strukturu

– b. časová závislost kros-termů nelze (ani přibližně) aproximovat klasickou Liouvillovou rovnicí: hbar je obsaženo v jejich definici ve jmenovateli vyšších derivací, takže kvantová Liouvillova rovnice ve fázovém prostoru diverguje s hbar!

– příklad (Heller):

kvantově klasick

y

Page 34: Klasick á-kvantová korespondence ve fázovém prostoru

str.

34

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

– pokud je počátečním stavem základní vibrační stav, je podmínka splněna. Pokud je třeba propagovat excitované vibrační stavy, je třeba je vyjádřit jako lineární kombinace Gaussiánů, které mají pozitívně definitní W.d. a autokorelační funkci počítat jako součet.

– rozdělení acf na jednotlivé rekurence:

– kros-termy na kvantové Wignerově distribuci představují informaci o relativním fázovém posuvu mezi „větvemi“ distribuce, které jsou odděleny body obratu. Tato informace chybí v klasicky propagované distribuci. Proto je třeba s jednotlivými větvemi pracovat odděleně, jakoby šlo o nezávislé Wignerovy distribuce.

– rekurence na autokorelační funkci vznikají díky navracejícím se částem klasické distribuce. Po určité době propagace se u neharmonických systémů navracejí trajektorie, které prošly různý počet bodů obratu – jejich příspěvek k acf je potřeba při klasické analýze započítat odděleně.

Page 35: Klasick á-kvantová korespondence ve fázovém prostoru

str.

35

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

• nákres vývoje balíku v neharmonickém potenciálu:

Page 36: Klasick á-kvantová korespondence ve fázovém prostoru

str.

36

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

tyto trajektorie se otočily 2x

tyto trajektorie se otočily 1x

příklad klasické propagace v dvojitém minimu:

Page 37: Klasick á-kvantová korespondence ve fázovém prostoru

str.

37

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

I.Horenko, M.Weiser, J Comp Chem 24 (2003) 1921

Page 38: Klasick á-kvantová korespondence ve fázovém prostoru

str.

38

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

– trajektorie, které se liší počtem bodů obratu se také zpravidla liší energií (perioda oscilace závisí na energii). Proto lze často rozdělit příspěvky k překryvu distribucí při přidání souřadnice energie – viz níže...

• u kvantové propagace naopak vznikají oscilační kros-termy. Pak není třeba oddělovat jenotlivé příspěvky (ani to není možné, pojem klasických trajektorií zde ztrácí smysl) Příklad propagace v Morseho potenciálu:

Page 39: Klasick á-kvantová korespondence ve fázovém prostoru

str.

39

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

M.Gronager, N.E.Henriksen, JCP 102 (1995) 5387.

Page 40: Klasick á-kvantová korespondence ve fázovém prostoru

str.

40

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

• překryv počáteční a konečné W.d. s rozlišenou energií (příklad pro 2D Pullen-Edmondsův potenciál)

současný příspěvektrajektorií, které se otočily5x, 6x a 7x

Page 41: Klasick á-kvantová korespondence ve fázovém prostoru

str.

41

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

• netradiční řešení problému: každou rekurenci je třeba brát zvlášť, spočítat její fázi a amplitudu, a sečíst rekurence jako komplexní funkce:

Page 42: Klasick á-kvantová korespondence ve fázovém prostoru

str.

42

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

– srovnání řešení běžného a netradiční metodou s referenční kvantovou autokorelační funkcí:

• Problém netradičního rozdělení na rekurence: známe pouze amplitudu a derivaci fáze pro každou rekurenci, to znamená, že chybí znalost fázového posuvu mezi rekurencemi

• řešení:

Page 43: Klasick á-kvantová korespondence ve fázovém prostoru

str.

43

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

– najdeme hlavní proudy trajektorií. Nejjednodušeji existuje jediný proud daný střední trajektorií počátečního Gaussiánu. U chaotických systémů však dochází k větvení (bifurkacím). Pak je třeba počítat s více proudy.

– pro každý proud vytvoříme jednu pseudo-trajektorii, která sleduje vývoj proudu ve fázovém prostoru: tato pseudotrajektorie má tu vlastnost, že je periodická (což platí pouze přibližně o proudu skutečných trajektorií)

– podél pseudotrajektorie umístíme hypotetické balíky. Jejich umístění podél pseudo-trajektorie je definováno časem tau, který určuje

Page 44: Klasick á-kvantová korespondence ve fázovém prostoru

str.

44

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“

Page 45: Klasick á-kvantová korespondence ve fázovém prostoru

str.

45

TMF0

45

le

tní s

emes

tr 2

006

X

TMF0

45

le

tní s

emes

tr 2

006

X

Wignerova metoda – „dynamické Monte Carlo“


Recommended