PRAHA 2001/2010/2012 STUDIJNÍ TEXT PRO ZVÍDAVÉ TTTFFF222::: KKKVVVAAANNNTTTOOOVVVÁÁÁ TTTEEEOOORRRIIIEEE PETR KULHÁNEK FEL ČVUT PŘEDMLUVA Teoretická mechanika vychází ze zobecněných zkušeností člověka, z toho, jak vnímáme svět kolem sebe v našich měřítkách – v tzv. makrosvětě. Snažíme-li se zákony teoretické mechaniky aplikovat na tělesa malých rozměrů (atomy, částice) – tzv. mikrosvět, nebudou již předpovědi ve shodě s experimentem. V mikrosvětě platí jiné zákony. Například samotný akt měření může ovlivnit objekty mikrosvěta. Chceme-li určit polohu fotbalového míče, zachytíme okem fotony odražené od míče a informaci zpracujeme. Chceme-li určit polohu elektronu, odražený foton, z kterého na polohu usuzujeme, udělí elektronu nezanedbatelný impuls a změní jeho stav. Asi největší rozdíl mezi jevy v makrosvětě a mikrosvětě souvisí s komutativností. V makrosvětě jsme si zvykli na to, že jevy, které pozorujeme, jsou komutativní – nezáleží na pořadí. Je jedno, zda nejprve provedeme měření A a poté měření B nebo naopak. Zkrátka AB = BA. V mikrosvětě tomu tak ale není. Akt měření ovlivňuje stav objektů a záleží na tom, které měření provedeme jako první. To je také hlavním důvodem selhání teoretické mechaniky při popisu mikrosvěta. Teoretická mechanika je založena na komutujících matematických objektech. Jedinou nekomutující strukturou jsou Poissonovy závorky, a to navíc ještě pomocnou. První jevy v mikrosvětě, které byly v příkrém rozporu s teoretickou mechanikou, byly objeveny na počátku 20. století. Jejich analýza vedla ke zrodu kvantové teorie – jedné ze dvou nejúspěšnějších teorií v dějinách lidstva (kvantová teorie, obecná teorie relativity). Základní rovnice a vztahy zůstávají shodné s teoretickou mechanikou, platí však pro zcela jiné objekty. Například Lieova algebra Poissonových závorek je aplikována na jisté operátory představující dynamické proměnné. Předpovědi dnešní kvantové teorie se shodují s experimentem na mnoho platných cifer. Uveďme nyní základní rozdíly světa malých rozměrů – mikrosvěta – oproti situacím, na které jsme zvyklí z našeho okolí – makrosvěta: 1) diskrétní hladiny některých dynamických proměnných (například energie, moment hybnosti ...) – v dané situaci můžeme naměřit jen určité hodnoty u sledované veličiny a žádné jiné. V makrosvětě jsou měřené hodnoty spojité. 2) dualismus vln a částic – objekty mikrosvěta se mohou chovat jako vlny i jako částice. 3) nekomutativnost aktu měření – při měření hodnot dvou dynamických proměnných (například polohy a rychlosti) může výsledek záležet na pořadí provedení měření. Akt měření totiž ovlivňuje stav systému, po měření se systém obecně nachází v jiném stavu než před měřením. 4) relace neurčitosti – zvýšení přesnosti měření jedné dynamické proměnné v některých případech sníží přesnost měření jiné dynamické proměnné. Tato měření se navzájem ovlivňují a jsou nekomutativní. 5) nedeterminismus kvantové teorie – dva experimenty připravené za stejných podmínek mohou dopadnout různě. Při provedení mnoha pokusů zjistíme, že výsledky mají pravděpodobnostní charakter. Jsme tedy schopni předpovědět jen to, s jakou pravděpodobností naměříme ten či onen možný jev, nikoli který jev konkrétně nastane. Fyzika se tak dostala před úlohu vytvořit takovou teorii, která by souhlasila s experimenty v mikrosvětě a v makrosvětě přecházela v klasickou teoretickou mechaniku. Konstrukcí kvantové teorie se budeme zabývat v této části sylabu. Aktuální verzi sylabu naleznete na serveru www.aldebaran.cz v sekci Studium. Nahrávky přednášek z roku 2015 naleznete na adrese https://www.youtube.com/playlist?list=PLBjIROWPicvWvL1yhucO4cQaL0jApTxir. Petr Kulhánek OBSAH 2. KVANTOVÁ TEORIE 5 2.1 VZNIK A VÝVOJ KVANTOVÉ TEORIE 5 2.2 (M) OPERÁTORY V KVANTOVÉ TEORII 9 2.2.1 UNITÁRNÍ PROSTORY (PROSTORY SE SKALÁRNÍM SOUČINEM) 9 2.2.2 OPERÁTORY 12 2.2.3 PROJEKČNÍ OPERÁTORY 17 2.2.4 ROZVOJ PRVKU DO BÁZE 19 2.2.5 SPEKTRÁLNÍ TEORIE 21 2.3 ZÁKLADNÍ PRINCIPY KVANTOVÉ TEORIE 26 2.3.1 ZÁKLADNÍ AXIOMY A DEFINICE 26 2.3.2 KOMPATIBILITA MĚŘENÍ A HEISENBERGOVY RELACE 30 2.3.3 VLASTNÍ STAVY ENERGIE, SCHRÖDINGEROVA ROVNICE 33 2.4 HARMONICKÝ OSCILÁTOR 36 2.4.1 ŘEŠENÍ POMOCÍ VLNOVÉ MECHANIKY (SCHRÖDINGER) 36 2.4.2 ŘEŠENÍ BEZ VOLBY REPREZENTACE (DIRAC) 40 2.4.3 ŘEŠENÍ POMOCÍ MATICOVÉ MECHANIKY (HEISENBERG) 43 2.5 SFÉRICKY SYMETRICKÝ POTENCIÁL 45 2.5.1 MOMENT HYBNOSTI 47 2.5.2 ŘEŠENÍ V X REPREZENTACI, KULOVÉ FUNKCE 50 2.5.3 JEDNODUCHÉ SYSTÉMY: OSCILÁTOR, VODÍK, JÁMA 52 2.6 ČASOVÝ VÝVOJ 54 2.6.1 EVOLUČNÍ OPERÁTOR 54 2.6.2 ČASOVÁ SCHRÖDINGEROVA ROVNICE 56 2.6.3 OSCILACE NEUTRIN 57 2.6.4 DVOUŠTĚRBINOVÝ EXPERIMENT 58 2.6.5 EHRENFESTOVY TEORÉMY, VIRIÁLOVÝ TEORÉM 59 2.7 RELATIVISTICKÁ KVANTOVÁ TEORIE, SPIN 62 2.7.1 PROSTOROVÁ ROTACE A LORENTZOVA TRANSFORMACE 62 2.7.2 SPIN 63 2.7.3 KLEINOVA-GORDONOVA ROVNICE 65 2.7.4 DIRACOVA ROVNICE 69 2.7.5 POZITRON, C SYMETRIE 79 2.7.6 ELEKTRON A JEHO POLE, U(1) SYMETRIE 81 2.8 SOUSTAVA STEJNÝCH ČÁSTIC 85 2.8.1 OPERÁTOR VÝMĚNY DVOU ČÁSTIC 85 2.8.2 BOSONY A FERMIONY, PAULIHO PRINCIP 86 2.8.3 DRUHÉ KVANTOVÁNÍ 87 2.8.4 UKÁZKA DRUHÉHO KVANTOVÁNÍ PRO KLEINOVO-GORDONOVO POLE 89 PŘÍLOHA – ZOBECNĚNÉ FUNKCE 92 P1. DIRACOVA DISTRIBUCE 92 P2. KONVOLUCE 94 P3. GREENŮV OPERÁTOR A GREENOVA FUNKCE 95 Kvantová teorie Vznik a vývoj 5 2. KVANTOVÁ TEORIE 2.1 VZNIK A VÝVOJ KVANTOVÉ TEORIE Shrňme nyní základní experimentální fakta, která vedla ke zrodu kvantové teorie: Záření absolutně černého tělesa: V absolutně černém tělese (lze za ně považovat například každou hvězdu) je v rovnováze látka a záření při nějaké konkrétní teplotě T. Sledujeme-li vyzařování absolutně černého tělesa, zjistíme, že na různých frekvencích vyzařuje s různou intenzitou. Experimentálně pozorovaný průběh energie vyzářené na jednotkovou frekvenci je na obrázku. Teoretické výpočty křivky záření absolutně černého tělesa, které prováděli Rayleigh, Jeans a Wien, vedly k odlišným závislostem. Buď divergovaly v infračervené (IR) nebo v ultrafialové (UV) oblasti spektra. Správnou formuli uhodl až Max Planck v srpnu 1900 tím, že zkoušel porovnávat různé funkce s naměřenými údaji. Jeho výsledek zněl: dI/d ~ 3 exp[– const /T ]. Za další dva měsíce odvodil Planck tuto závislost i teoreticky za předpokladu, že energie světla o určité frekvenci  se nemění spojitě, ale je celistvým násobkem základního energetického kvanta .1005.1; 34 Js  E (2.1) Veličina  se nazývá redukovaná Planckova konstanta. Planck původně použil předpoklad o kvantování energie pro zjednodušení matematických výpočtů. Později se ukázalo, že energie elektromagnetického záření určité frekvence je skutečně kvantována, tj. její pozorované hodnoty nejsou spojité, ale mění se skokem o základní energetické kvantum  . Fotoelektrický jev (fotoefekt): Při dopadu světla (elektromagnetického záření) na povrch kovu může být z kovu vytržen elektron, který opustí povrch kovu. K uvolňování elektronů z kovu dochází při frekvencích světla vyšších než prahová frekvence 0, která je pro daný kov charakteristická. Máme-li k dispozici světlo s frekvencí nižší než prahovou, emise elektronů nenastane, byť bychom použili světlo se sebevětší intenzitou. Tento experiment je v rozporu s představou o světle jako elektromagnetickém vlnění. K fotoefektu by mělo docházet při každé frekvenci a dostatečnou energii k emisi by mělo jít získat zvýšením intenzity dopadajícího světla. Řešení podal A. Einstein v roce 1905. Elektromagnetické vlnění se chová při fotoefektu jako částice. Tyto částice nazval fotony. Energie jednoho fotonu záření o frekvenci  je právě energie jednoho energetického kvanta (2.1). Vysvětlení fotoelektrického jevu je nyní velice jednoduché. Na povrchu kovu dochází ke srážce fotonu s elektronem. Aby foton vyrazil elektron, musí mít vyšší energii než je vazbová energie elektronu v kovu: iE  . Prahová frekvence zřejmě je 0 iE   . Celková energetická bilance 2 i e 1 2 E m   v se nazývá Einsteinova rovnice pro fotoefekt. Energie dopadlého fotonu se spotřebuje na vytržení elektronu z kovu a na kinetickou energii vylétávajícího elektronu. IR UV dI d/   záření elektron kov Kvantová teorie Vznik a vývoj 6 Elektromagnetické vlnění tedy můžeme považovat za soubor fotonů. Proto i při záření absolutně černého tělesa se mění energie záření o dané frekvenci skokem – tento skok představuje přírůstek nebo úbytek jednoho fotonu. Comptonův jev A. H. Compton v roce 1923 zjistil, že rentgenové paprsky odražené od povrchu grafitu mění svoji vlnovou délku. Podle klasických představ by vlny měly rozkmitat povrchové elektrony a ty generovat vlnu se stejnou frekvencí. Vysvětlení: Fotony se opět chovají jako částice, srážejí se s elektrony a při srážce ztrácí část energie, a proto mění svou vlnovou délku. Ohyb elektronů: Fotoelektrický jev ukázal, že vlnění se může chovat v určitých situacích jako částice. Naopak, někdy se částice chovají jako vlny. Například svazek elektronů procházející štěrbinou nebo dvouštěrbinou po dopadu na stínítko vytvoří typický ohybový obrazec. Nemůžeme předem říci, kam který elektron dopadne, ale při velkém množství elektronů můžeme určit pravděpodobnosti dopadu do konkrétního místa na stínítku. Vzniklý ohybový obrazec je tedy typickým statistickým jevem. Dnes jsou vlnové vlastnosti elektronů využívány například v elektronových mikroskopech. Elektrony mají výrazně kratší vlnovou délku než viditelné světlo a proto je rozlišovací schopnost elektronového mikroskopu podstatně vyšší než optického. Poprvé byly vlnové vlastnosti elektronu pozorovány C. J. Davissonem a L. H. Germerem v roce 1927. Zkoumali odraz elektronů od povrchu niklu. Po vyžíhání niklu došlo k rekrystalizaci a odražené elektrony začaly vykazovat na přesných velkých krystalech ohybový obrazec. Poznámka: Částice popisujeme čtveřicí veličin (E, p). Definice energie E a hybnosti p souvisí se symetriemi při posunutí v čase a v prostoru (teorém Noetherové). Vlnění popisujeme čtveřicí veličin (, k). Úhlová frekvence  je definována jako změna fáze vlnění s časem  = /t a vlnový vektor k je změna fáze vlnění s prostorovými souřadnicemi k = /x. Při periodickém ději s konstantní periodou T v čase a  v prostoru (vlnová délka) lze psát  = 2/T, k = 2/Louis de Broglie vyslovil hypotézu, že objekty mikrosvěta se chovají jako vlny i jako částice (dualismus vln a částic). Převodní vztah má tvar: ., kp   E (2.2) Často nás zajímá vlnová délka vlnění odpovídajícího konkrétní částici, například elektronu v elektronovém mikroskopu. Ze vztahu (2.2) máme  2vm a tedy . 2 vm    (2.3) Existence atomu: Podle klasického planetárního modelu atomu obíhají záporně nabité elektrony kolem kladně nabitého jádra tak, jako ve Sluneční soustavě obíhají planety kolem Slunce. Odstředivá síla je vyrovnána přitažlivou Coulombovou silou. Mezi gravitačními a elektromagnetickými jevy je ale podstatný rozdíl. Z Maxwellovy teorie elektromagnetického pole plyne, že každá nabitá částice, která se pohybuje se zrychlením, vyzařuje elektromagnetické vlnění a ztrácí tak elektrony štěrbina počet elektronů stínítko Kvantová teorie Vznik a vývoj 7 energii. Při kruhovém pohybu elektronu kolem jádra se mění směr rychlosti, zrychlení dv/dt je nenulové (míří do centra atomu, jde o dostředivé zrychlení) a elektron ztrácí energii zářením. Pohybuje se po spirále, až dopadne na jádro atomu. Tento proces trvá například pro vodík 10–11 s. Podle klasické teorie by tedy za velice krátkou dobu neměly žádné atomy existovat!! Na tento paradox upozornil poprvé dánský fyzik Niels Bohr. Niels Bohr vytvořil tzv. Bohrův model atomu na základě tří umělých postulátů, které přidal ke klasické teorii: 1) elektrony se pohybují jen po tzv. stacionárních drahách – tj. po takových drahách, ve kterých je odpovídající de Broglieho vlnová délka ze vztahu (2.3) „namotána“ na oběžnou dráhu tj. obvod dráhy je n-násobkem vlnové délky. ! Tato dráha není možná Tato dráha je možná . 2 ;2 n n m nr v    Index n čísluje možné stavy elektronu v atomu (rn možný poloměr dráhy, vn rychlost na n-té dráze, En odpovídající energie) podle počtu vlnových délek elektronu na jeho oběžné dráze. 2) na stacionární dráze elektron nezáří. 3) při přeskoku elektronu mezi dvěma stacionárními hladinami dojde k vyzáření fotonu o energii odpovídající rozdílu energií těchto hladin. Tento jednoduchý Bohrův model atomu není řešením výše uvedeného paradoxu, jde spíše o postulování nebo konstatování experimentálně známých skutečností. Navíc je tento model aplikovatelný jen na nejjednodušší atomy s jediným elektronem v obalu (H, He+). Tento model ale poprvé správně určil hladiny energie elektronu v atomu vodíku a vysvětlil spektrum atomu vodíku. Heisenbergovy relace neurčitosti: Při měření polohy a hybnosti objektu mikrosvěta budou nepřesnosti měření x, p splňovat relaci (přes k se nesčítá) .3,2,1; 2  kpx kk  (2.4) Čím přesněji určíme polohu objektu, tím méně přesně určíme jeho hybnost a naopak. Samotný akt měření ovlivňuje náš objekt, ale relace (2.4) je splněna i tehdy, neprovedeme-li měření vůbec. Jde o principiální hranici danou přírodou, za kterou nelze nahlédnout. Například obyčejný ohyb světla na štěrbině lze chápat jako důsledek relací neurčitosti pro fotony. Průchod fotonů štěrbinou není nic jiného než pokus o určení jejich polohy y s přesností y (velikost štěrbiny). Fotony, které prošly štěrbinou, určitě měly v okamžiku průchodu souřadnici y rovnou souřadnici y štěrbiny. Zmenšíme-li šířku štěrbiny y, zvýšíme přesnost měření y; podle relací (2.4) se ale zvýší nepřesnost py určení odpovídající komponenty hybnosti. Výsledkem je známý ohybový jev – fotony za štěrbinou vyletují s danou pravděpodobností do různých směrů se střední kvadratickou fluktuací hybnosti py danou Heisenbergovými relacemi neurčitosti. Kvantová teorie Vznik a vývoj 8 Výčet experimentálních faktů, které jsme uvedli výše, není zdaleka úplný. Všechny ale přispěly ke zrodu kvantové teorie, popisující pro nás nezvyklý svět atomů a elementárních částic. Podejme nyní stručný přehled jejího vývoje. V roce 1925 formuloval Werner Heisenberg ve svých 25 letech maticovou mechaniku – každé dynamické proměnné přiřadil čtvercovou matici (zpravidla nekonečnou), jejíž vlastní čísla byly měřitelné hodnoty příslušné veličiny. Šlo o teorii pramenící z vynikající intuice, na základě které bylo možné určit například energetická spektra různých atomů (nejen vodíku). V roce 1926 Erwin Schrödinger formuloval vlnovou kvantovou mechaniku. Řešením rovnice  EV m            2 2  (2.5) pro vlnovou funkci  bylo opět možné určit hodnoty energie E pro objekt v potenciálním poli V(x, y, z). Obě konstrukce – Heisenbergova i Schrödingerova – poskytovaly shodné výsledky. Spíše než o ucelenou teorii šlo v té době o návod, jak určit energetické spektrum. Obecnou konstrukci kvantové teorie na Hilbertových prostorech provedl P. A. M. Dirac. Ukázalo se, že Heisenbergova a Schrödingerova mechanika se liší jen jinou volbou příslušného Hilbertova prostoru. Až doposud byla budována nerelativistická kvantová teorie. Zobecnění na relativistický případ provedli Klein a Gordon pro spin částice s = 0 a Dirac pro spin částice s = 1/2 (Kleinova-Gordonova rovnice, Diracova rovnice). S relativistickou kvantovou teorií byla objasněna podstata spinu, Dirac předpověděl existenci pozitronu, ale především byl postaven základ pro vybudování kvantové elektrodynamiky (Dirac – 1949). Odsud byl již jen krůček ke vzniku kvantové teorie elektromagnetického pole (Dirac, Feynman), ve které dochází i ke kvantování samotného elektromagnetického pole (tzv. druhé kvantování). Výsledky kvantové teorie pole lze přehledně zapisovat pomocí tzv. Feynmanových diagramů. Na základě různých symetrií v přírodě se od 60. let bouřlivě vyvíjí kalibrační teorie, například Weinbergova-Salamova teorie elektroslabé interakce, která sjednocuje teoretický pohled na interakci elektromagnetickou a slabou, rozvíjí se kvantová chromodynamika – teorie silné interakce, teorie GUT sjednocující elektroslabou a silnou interakci a probíhají intenzívní pokusy o formulaci Einsteinových-Diracových rovnic supersymetrických teorií SUSY pokoušejících se o jednotný popis všech čtyř interakcí. Lidstvo stále více poznává svět elementárních částic a jeho zákonitosti. Následující kapitolu budeme věnovat matematice, kterou je třeba znát pro pochopení kvantové teorie. Vlastní stavbou kvantové teorie se budeme zabývat až v kapitole 2.3 a následujících. Kvantová teorie Operátory 9 2.2 (M) OPERÁTORY V KVANTOVÉ TEORII V této kapitole se budeme zabývat nejdůležitější matematikou potřebnou v kvantové teorii. Veškeré úvahy jsou z důvodu jednoduchosti provedeny pro případ, kdy vlastní čísla operátorů jsou navzájem různá a tvoří spočetnou množinu. Obecnější případy vícenásobných vlastních čísel a spojitého spektra jsou krátce diskutovány v závěru kapitoly. 2.2.1 Unitární prostory (prostory se skalárním součinem) V kapitole 1.4.1 jsme rozšířili pojem vektoru na obecnější objekty než jsou uspořádané trojice a zavedli lineární vektorový prostor. Pozorně si znovu tuto pasáž přečtěte! Nyní analogicky rozšíříme pojem skalárního součinu pro různé lineární vektorové prostory. Budeme důsledně používat Diracovu symboliku, ve které jsou prvky lineárních vektorových prostorů značeny symboly  axf |,|,| a skalární součiny  bayxgf |,|,| , atd. R 3 prostor reálných trojic .| ,),,(|,),,(| 332211 321321 součinskalárníkk gfgfgfgf gggfff   gf gf Norma vektoru (velikost) se definuje vztahem .||||| 2 3 2 2 2 1 pro 3 fff  R fff (2.6) Pro reálné trojice znázorněné jako úsečky opatřené šipkami je norma vektoru rovna délce úsečky a platí cos|||||||||  gfgf , kde  je úhel sevřený oběma vektory. Z tohoto vztahu plyne okamžitě Schwartzovo lemma: | || || || || .   f g f g (2.7) Výsledkem operace skalárního součinu je číslo, v případě lineárního vektorového prostoru R3 reálné číslo, v obecném případě bude výhodné uvažovat i o čísle komplexním. Norma vektoru (velikost) musí ale vždy být nezáporné reálné číslo. R N prostor reálných N-tic .| ,;),,(|,),,(| 1 11 11 kk N k kkNN llNN gfgfgfgf gfggff       gf gf ,R V platnosti zůstávají definice normy i Schwartzovo lemma. C N prostor komplexních N-tic .| ,;),,(|,),,(| * 1 ** 1 * 1 11 kk N k kkNN llNN gfgfgfgf gfggff       gf gf ,C Skalární součin definujeme v jednom z argumentů komplexně sdružený (dohodou v levém). Pro komplexní číslo z = a + i b je velikost (norma) čísla dána vztahem * z z z Kvantová teorie Operátory 10 Právě proto, aby pro komplexní čísla zůstalo v platnosti, že norma vektoru je odmocnina skalárního součinu vektoru se sebou samým, je v definici skalárního součinu komplexní sdružení v jednom z argumentů. Při výše uvedené definici skalárního součinu bude výsledkem sice komplexní číslo, ale norma vektoru zůstane reálná nezáporná: * * 2 2 1 1 1 0N N Nf f f f f f f f f         . Opět platí Schwartzovo lemma. l 2 prostor komplexních posloupností (N-tice s N ) 1 1 1 * * * * 1 1 1 | { , , , } { } , | { } ; , | . n l l l l l l n n k k k k k f f f g f g f g f g f g f g                     f g f g     C , Takto definovaný skalární součin má smysl jen pro konvergentní posloupnosti. Do prostoru l 2 můžeme zahrnout jen takové prvky f| , pro které je |||| f , tj. požadujeme .|pro 2 1 * lff k kk    f Potom je ,|,|pro||||||||||||| 2 1 * lgf k kk     gfgfgf neboť Schwartzovo lemma platí i v případě nekonečných posloupností. L2 () prostor komplexních funkcí reálné proměnné Při dalším zobecnění prostoru l 2 si můžeme index k představit spojitý. Místo k budeme psát x : fx. Výraz fx není ale nic jiného než komplexní funkce reálné proměnné (spojitého indexu), kterou je zvykem zapisovat ve tvaru f (x) , tj. .)()(| ,,;),(|),(| *      dxxgxf gfxxggxff xx gf gf ,CR Analogicky jako v l 2 je třeba do prostoru L2 zahrnout jen prvky s |||| f , tj. .)(pro)()( 2* L   xfdxxfxf Potom je 2* |,|pro|||||||||)()(|||| L    gfgfgf dxxgxf a skalární součin má smysl. Schwartzovo lemma platí i pro integrály. L2 se někdy nazývá prostor funkcí integrovatelných s kvadrátem. Lze ho definovat i pro jiný definiční obor než (   ), potom píšeme L2 (M), kde M je definiční obor funkcí f (x)  L2 (M). Nyní můžeme přistoupit k obecné definici prostorů se skalárním součinem. Kvantová teorie Operátory 11 UNITÁRNÍ PROSTOR (prostor se skalárním součinem) – unitárním prostorem nazveme lineární vektorový prostor V (s operací + : V  V V a operací  : C   V  V ), na kterém je definována další operace C VV:| (tzv. skalární součin) s vlastnostmi 0|0;0)4 ,)||(|)3 ,||)2 ,|||)1 **     fffff gfgfgffg gfgf hfgfhgf   Poznámky: 1) Přidáním operace [ , ] z lineárního vektorového prostoru získáme Lieovu algebru, přidáním operace < | > získáme unitární prostor. 2) První dvě operace v definici znamenají linearitu v pravém argumentu. Z třetí operace plyne antilinearita v levém argumentu (aditivnost + vytknutí komplexně sdružené konstanty). 3) Symbolika zápisu pochází od P. A. M. Diraca. Nazývá se také braketová symbolika nebo brakety (z anglického bracket = závorka). < | > „bracket“ < | „bra“ (lze matematicky definovat, duál, naznačená operace skalárního součinu) | > „ket“ (vektor z V ) 4) Pro komplexní N-tice lze interpretovat f| jako sloupcovou matici, |f jako transponovanou komplexně sdruženou matici:   .; ** 1 1 N N ff f f              ff Potom je skalární součin   kk N N gf g g ff * 1 ** 1               gf definován za pomoci maticového násobení. Pro jiné prostory než n-tice není pro naše účely třeba jednotlivé části skalárního součinu  gf nějak interpretovat. 5) Pro L2 lze chápat     dxxf )(| * f jako naznačenou operaci skalárního součinu. Je jen třeba doplnit patřičnou funkci, na kterou operace působí. Podobná situace je u derivování, napíšeme-li jen d /dx. Paprsek – nechť V je unitární prostor, f| jeho nenulový prvek. Paprskem nataženým na f| nazveme množinu prvků }|,|},0{\;||;|{ VC  gffgg  . | f > Hilbertův prostor – úplný unitární prostor (hranice prostoru je jeho součástí). Separabilní Hilbertův prostor – Hilbertův prostor se spočetnou bází. Kvantová teorie Operátory 12 2.2.2 Operátory Operátorem rozumíme zobrazení ,:ˆ VV A které prvku f| prostoru V přiřazuje prvek g| tohoto prostoru: .||ˆ  gfA V platnosti zůstává běžné názvosloví používané pro zobrazení (vzor, obraz, definiční obor, obor hodnot,…). Příklad 1: R3 Operátorem na R3 může být libovolná matice 33, například . 110 100 001 |ˆ obecně,| 3 1 1 1 2 1 110 100 001 |ˆ 1 2 1 |; 110 100 001 ˆ 32 3 1 3 2 1                                                                                            ff f f f f f f gf f A A A Příklad 2: L 2 ()   .|e)1(e|ˆ e|;ˆ     gf f xx x xx dx d x dx d D D Jednotkový operátor:  ff ||ˆ1 . Pro n-tice je jednotkovým operátorem diagonální matice s jednotkami na diagonále (jednotková matice) – ověřte! Kvadrát operátoru: Druhou mocninu operátoru můžeme definovat, je-li obor funkčních hodnot operátoru podmnožinou jeho definičního oboru, potom   .|ˆˆ|ˆ 2  ff AAA Příklad 3: Operátor derivace .e)42(e2e|ˆ e|;ˆ 222 2 22 xxx x xx dx d dx d dx d dx d                f f D D Mocnina operátoru: Analogicky definujeme indukcí obecnou mocninu operátoru   .|ˆˆ|ˆ 1  ff -nn AAA Funkce operátoru: Nechť f (x) je analytická funkce s Taylorovým rozvojem .)( 0     k k k xcxf Kvantová teorie Operátory 13 Potom můžeme definovat funkci operátoru .ˆ)ˆ( 0     k k kcf AA (2.8) Připomeňme si zde rozvoje některých důležitých funkcí: . !8!6!4!2 1ch , !9!7!5!3 sh , !8!6!4!2 1cos , !9!7!5!3 sin , !4!3!2 1e ,1 1 1 8642 9753 8642 9753 432 5432              xxxx x xxxx xx xxxx x xxxx xx xxx x xxxxx x x Příklad 4: Na prostoru C 2 je zadán maticový operátor 0 iˆ i 0        A . Určete exp( ˆA ). 2 3 2 4 3 2 2 1 2 0 i 0 i 1 0ˆ ˆ , i 0 i 0 0 1 ˆ ˆ ˆ ˆ ˆ ˆ , ˆ ˆ ˆ ˆ ˆ ˆ ˆ , ˆ ˆ ˆ ˆ, , 1,2, .n n n                                       A 1 A A A A 1 A A A A A A A 1 A A A 1 Nyní již snadno nalezneme hledanou funkci matice: . 1ch1shi 1shi1chˆ)1sh(ˆ)1ch( ˆ !7 1 !5 1 !3 1 1ˆ !6 1 !4 1 !2 1 1 !4 ˆ !3 ˆ !2 ˆ ˆ1)ˆ(exp 432                       A1 A1 AAA AA   Takto získají smysl například i výrazy typu sin(d/dx ) a podobně. Později se naučíme funkci operátoru nalézt pomocí spektrálního rozvoje operátoru. Jde o efektivnější způsob než je Taylorův rozvoj. Inverzní operátor: Inversním operátorem k Aˆ nazveme takový operátor 1ˆ  A , že .ˆˆˆˆˆ 11 1AAAA   Kvantová teorie Operátory 14 K danému operátoru Aˆ může být nalezení inversního operátoru značně obtížné, někdy inversní operátor neexistuje vůbec. Sdružený operátor: Sdruženým operátorem k Aˆ nazveme takový operátor †ˆA , že †ˆ ˆ| | .    f g f gA A Působení operátoru Aˆ v pravé straně skalárního součinu dopadne stejně jako působení k němu sdruženého operátoru v levé části skalárního součinu. Sdružený operátor k Aˆ nemusí vždy existovat. Příklad 5: 1 †1 2i 1 2 1 0ˆ ˆ ˆ; ; . 0 i 0 i 2i i                        A A A Skutečně 1 1 1 1 2 2 2 1 1† 2 1 2 1 1 2i 1 2 1 0ˆ ˆ ˆ , 0 i 0 i 0 1 1 2 1 2i 1 0ˆ ˆ ˆ ; 0 i 0 i 0 1 2i1 2iˆ | , i0 i 1 0ˆ | 2i i2i i ˆ| g g g g g f f f f f f                                                                                           g f f g A A 1 A A 1 A A A     1 2* * * * * 2 1 1 1 2 2 2 2 1† * * * * * * 1 1 2 1 1 1 2 2 2 2 2i 2i i , i ˆ | 2i i 2i i . g g f f g f g f g g g f f f f g f g f g g                       f gA Poznámka: Nalézt sdružený operátor pro matice je velmi snadné, původní matici stačí komplexně sdružit a transponovat (překlopit kolem diagonály), tj. † *ˆ ˆ( )T A A . Uveďme nyní velmi užitečné vztahy pro výpočet inversního a sdruženého operátoru pro součin dvou operátorů: ,ˆˆ)ˆˆ( 111   ABBA (2.9) † † †ˆ ˆ ˆ ˆ( ) .AB B A (2.10) Jejich důkaz je triviální přímo z definice inversního a sdruženého operátoru: ,ˆ/ˆˆ)ˆˆ( ,ˆ/ˆˆˆ)ˆˆ( 111 11 zprava zprava     ABABA B1BABA .ˆˆ)ˆˆ( 111   ABBA Analogicky postupujeme i pro sdružený operátor: † † † †ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ( ) | | | | .          f g f g f g f gAB AB A B B A Kvantová teorie Operátory 15 Komutativita operátorů: Pro operátory je obecně ABBA ˆˆˆˆ  . Říkáme, že operátory nekomutují. Míru nekomutativnosti můžeme posoudit za pomoci tzv. komutátoru ˆ ˆ ˆ ˆ ˆ ˆ[ ] .,  A B AB BA (2.11) Je-li komutátor operátorů Aˆ , Bˆ nulový, operátory komutují, je-li různý od nuly nekomutují. Výsledkem komutátoru je opět operátor. Uveďme nejdůležitější vlastnosti komutátorů (dokažte!) .0ˆˆ[ˆ[ˆˆ[ˆ[ˆˆ[ˆ[)4 ,ˆˆ[ˆˆ[)3 ,ˆˆ[ˆˆ[ˆˆˆ[2) ,ˆ,ˆ[ˆˆ[)1     ]]]]]] ]] ]]] ]] BACACBCBA BABA CBCACBA ABBA ,,,,,, ,, ,,, ,  To jsou ale právě definiční vlastnosti Lieovy algebry (1.22) až (1.25). Komutátory tvoří Lieovu algebru na prostoru operátorů. Příklad 6: Mějme na L2 dva operátory: xdxd  XD ˆˆ a , například na prvek 5 |x působí takto .|ˆ,5|ˆ 655455 xxxxxx dx d x  XD Určeme jejich komutátor  ˆ ˆ ˆ ˆ ˆ ˆ[ , ] | ( ) | ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) | d d d d x x f x x f x x f x dx dx dx dx f x x f x x f x f x                  f f f D X DX XD .1XD XD ˆ]ˆ,ˆ[ |pro||]ˆ,ˆ[ 2   Lfff Podobně můžeme určovat i další komutační relace. V celém tomto textu se budeme zabývat lineárními operátory, tj. operátory s lineární odezvou:  gfgf |ˆ|ˆ)||(ˆ AAA  . Všechny dosud uvedené operátory byly lineární. V kvantové teorii se setkáme především se dvěma druhy lineárních operátorů - operátory unitárními a operátory Hermitovými. Uveďme nyní definice těchto operátorů. UNITÁRNÍ OPERÁTORY Definice: unitární operátor zachovává skalární součin, tj. .ˆ|ˆ||ˆ|,|ˆ|  gfgfggff UUUU (2.12) Skalární součin se před a po působení unitárního operátoru nezmění. Věta: Pro unitární operátory je sdružený a inversní operátor totožný, tj. † 1ˆ ˆ . U U (2.13) Důkaz: Z definice sdruženého operátoru víme, že †ˆ ˆ ˆ ˆ| | .    f g f gU U U U Kvantová teorie Operátory 16 Pro zachování skalárního součinu (definice unitárního operátoru) je nutné, aby †ˆ ˆ ˆ ,U U 1 to ale podle definice inversního operátoru právě znamená, že † 1ˆ ˆ  U U . Příklad 7: operátor xi eˆ U na prostoru L 2 je unitární   .|)()( )(e)(e)(e)(eˆ|ˆ ,)(e|ˆ,)(| ,)(e|ˆ,)(| * i*ii*i i i        gf gf gg ff dxxgxf dxxgxfdxxgxf xgxg xfxf xxxx x x UU U U HERMITOVY OPERÁTORY Definice: Hermitův operátor působí v obou částech skalárního součinu stejně, tj. .ˆ||ˆ  gfgf AA (2.14) Věta: Pro Hermitův operátor je sdružený operátor roven původnímu, je samosdružený: †ˆ ˆ .A A (2.15) Důkaz: Plyne okamžitě z definice sdruženého operátoru. Poznámka: V přesné matematice se definice samosdruženého a Hermitova operátoru nepatrně liší požadavky na definiční obor, pro naše účely nebudeme samosdružené a Hermitovy operátory rozlišovat. Vzhledem k tomu, že Hermitův operátor působí v obou částech skalárního součinu stejně, často se píše .ˆ|ˆ||ˆ  gfgfgf |AAA Centrální pozice Aˆ naznačuje, že podle vlastního uvážení můžeme operátorem zapůsobit vlevo či vpravo. Tato struktura se někdy nazývá sendvič. Příklad 8: operátor xd d iˆ B na prostoru L 2 (  ) je hermitovský   .ˆ| )( i)( )( )(i)()(i )( )( i)()(i|ˆ *** .**                           gf gf B B dx xd xgd xfdx xd xgd xfxgxf dxxg xd xfd dxxgxf xd d partesp Výraz v hranaté závorce je nulový, neboť funkce z L 2 (  ) jsou integrovatelné s kvadrátem na (  ) a proto musí být .,pro0)(lim)(lim 2 L  gfxgxf xx Samotný operátor derivace dxdDˆ hermitovský není, při provedení integrace per partes by se zaměnilo znaménko a  gfgf DD ˆ||ˆ . Kvantová teorie Operátory 17 2.2.3 Projekční operátory Cílem této kapitoly bude naučit se nacházet projekce vektorů do předem zadaného paprsku. Jde o úlohu, která má prvořadý význam nejen pro kvantovou teorii. Rozvoje do různých typů řad (například Fourierova řada) nejsou nic jiného než hledání projekcí zadané funkce do vektorů nějaké báze, které reprezentují paprsek v prostoru. Z celého paprsku stačí vzít jediný vektor, který paprsek zcela popíše. Vybereme tento „reprezentativní“ vektor jednotkový, tj. tak, aby 1tj  ||||.,1| aaa . Snadno představitelná je situace v R2. Na obrázku vidíme jednotkový vektor a| , který representuje paprsek a vektory  hgf |,|,| , které do tohoto paprsku budeme promítat. | >g P | >f | >f | >h P | >h P| >g | >a paprsek ^ ^ ^  Projekce libovolného vektoru f| má velikost cos|||| f a směr ||||| aa  . Znaménko funkce cos ve velikosti projekce určuje, zda promítaný vektor míří ve směru a| nebo ve směru opačném. Projekci lze napsat jako součin její velikosti a směru: .| | | | |||||||| | |||| | |||||||| | |||| |||| | cos|||||ˆ               a aa fa a aa fa a a fa fa f a a ff P Povšimněte si, že při úpravách výrazů v Diracově symbolice můžeme libovolně stěhovat čísla (normy vektorů a skalární součiny). Výraz pro projekci se skládá z koeficientu, který určuje velikost projekce a vektoru a| . Zapíšeme-li formálně koeficient až za vektor, získáme operátorový tvar: . | || | | ||ˆ       aa faa aa fa afP Označme . | ||ˆ    aa aa P (2.16) Tento výraz se nazývá projekčním operátorem. Sám o sobě význam nemá, jde o naznačenou operaci skalárního součinu, která musí být vykonána. Teprve zapůsobením Pˆ na nějaký vektor f| dostaneme smysluplný výraz – projekci vektoru danou koeficientem  aafa || a směrem vektoru a| . Situace je podobná operátoru d/dx, také jde jen o naznačenou derivaci, která musí být vykonána na konkrétní funkci. Bude-li vektor a| , do kterého promítáme jednotkový, výrazy se ještě zjednoduší: .|||ˆ ,||ˆ   faaf aa P P (2.17) Koeficient projekce je  fa | a směr je a| . Kvantová teorie Operátory 18 Příklad 9: Nalezněte projekci vektoru f| do vektorů a| a b| . Vektory jsou dány takto: . 1 1 |;|;|                      baf 1 1 3 1 Nejprve nalezneme projekční operátory:         . 2121 2121 2 11 11 1 1 11 11 1 1 | ||ˆ , 2121 2121 2 11 11 1 1 11 11 1 1 | ||ˆ                                                                       bb bb aa aa b a P P Nyní snadno nalezneme hledané projekce: . 1 2121 2121 |ˆ , 2121 2121 |ˆ                                           13 1 2 2 3 1 f f b a P P Využijme nyní tohoto příkladu a vyzkoušejte si, že platí jednoduché a užitečné relace. Výpočty jsou natolik snadné, že zde uvedu jen výsledky: 0 1 1 2 3 2 | >a| >b | >f 1. † †ˆ ˆ ˆ ˆ; .a a bb P P P P 2. .ˆˆ;ˆˆ 22 bbaa PPPP  3. .ˆˆˆ 1PP  ba První relace znamená, že projekční operátory jsou hermitovské. Pro matice je význam jednoduchý: Matice se po překlopení kolem hlavní diagonály a následném komplexním sdružení nezmění. Druhá relace má také velmi jednoduchý význam: Projekce dvakrát provedená po sobě (kvadrát operátoru) je shodná s projekcí provedenou jednou. Obě vlastnosti jsou pro projekční operátory charakteristické a většinou se považují za definici projekčního operátoru: Definice: Projekčním operátorem nazveme takový lineární operátor, který je hermitovský a platí PP ˆˆ 2  . Poznámka: Snadno lze ukázat, že obě vlastnosti jsou splněny pro definici (2.16), například pro druhou vlastnost: .ˆ | || || ||| | || | ||ˆ 2 PP              aa aa aaaa aaaa aa aa aa aa V prostředním výrazu jsme zkrátili skalární součin v čitateli (uprostřed) s jedním ze skalárních součinů ve jmenovateli. Jde o prostá komplexní čísla, která lze vytknout před výrazy a krátit. Kvantová teorie Operátory 19 Význam třetí relace je také snadno pochopitelný: Vektory a| a b| jsou navzájem kolmé a v rovině tvoří ortogonální bázi (bázi složenou z kolmých vektorů). Projekce do těchto vektorů nejsou ničím jiným než rozkladem původního vektoru do této báze. Zkuste si obě projekce sečíst. Dostanete původní vektor. Právě matematickým vyjádřením faktu, že součet všech projekcí dá původní vektor je třetí relace: .||ˆ|ˆˆˆˆ  fff baba PP1PP Někdy se této relaci říká relace úplnosti. Pokud je daná báze úplná (nechybí v ní žádný vektor), potom je součet všech projekčních operátorů roven jednotkovému operátoru. To znamená, že součet všech projekcí libovolného vektoru dá původní vektor. 2.2.4 Rozvoj prvku do báze Mějme v unitárním prostoru spočetnou bázi (maximální množinu lineárně nezávislých vektorů) }{| ke . Vhodnou lineární kombinací jednotlivých prvků můžeme vždy zajistit, aby prvky báze byly navzájem kolmé a měly jednotkovou velikost. Takové prvky budeme označovat jen pořadovým číslem: k| . Báze složená z vektorů k| má dvě základní vlastnosti: ,| kllk  (2.18) .ˆ|| 1 k kk (2.19) Vlastnost (2.18) je vyjádřením ortonormality. Skalární součin dvou různých vektorů je nulový (jsou navzájem kolmé) a dvou stejných je roven jedné (vektory báze mají jednotkovou velikost). Vlastnost (2.19) je relací úplnosti báze. Součet všech projekčních operátorů dá jednotkový operátor, tj. součet všech projekcí libovolného vektoru dá původní vektor. V relaci úplnosti je možné vynechat sumaci a využít Einsteinovu sumační konvenci. V případě neseparabilních prostorů s nespočetnou bází mají obě relace tvar: | ( ) ,x y y x    (2.20) .ˆ|| 1 dxxx (2.21) Na pravé straně relace ortonormality je místo Kroneckerova symbolu Diracova  distribuce a v relaci úplnosti je místo sumace projekčních operátorů integrace. Rozvoj prvku do báze je v Diracově symbolice mimořádně jednoduchý. Stačí před prvek vsunout relaci úplnosti v podobě jednotkového operátoru: .|,||||ˆ|   ffff kckckk k k k k 1 (2.22) Koeficienty rozvoje jsou dány skalárním součinem rozvíjeného vektoru s prvkem báze. Příklad 10: Fourierova řada Uvažujme prostor L2 (0, 2) periodických funkcí integrovatelných s kvadrátem. Díky požadavku f(0) = f(2) tvoří v tomto prostoru úplnou bázi soustava funkcí (viz elementární učebnice z matematiky nebo látka z prvních semestrů): .,2,1,0,e| i  kf kx k Tyto funkce jsou sice navzájem kolmé, ale nejsou jednotkové: Kvantová teorie Operátory 20 .2 0 eee| 2 0 )(i 2 0 ii lk xklxlxk lk kl kl dxdxff             pro2 pro Nulovost skalárního součinu pro l  k plyne z periodičnosti trigonometrických funkcí na intervalu <0, 2>. Vydělíme-li prvky báze jejich velikostí, získáme ortonormální bázi ,,2,1,0,e 2 1 | i  kk kx  pro kterou platí relace ortonormality (2.18) a relace úplnosti (2.19). Rozvoj libovolné funkce z našeho prostoru je potom ,|,||   ff kckc k k k neboli ,)(e 2 1 ,e 2 1 )( 2 0 ii      dxxfccxf kx k k kx k což jsou známé vztahy pro Fourierovu řadu. Nezapomeňte, že skalární součin je v levém argumentu komplexně sdružený, proto je v koeficientu ck minus. Reprezentace V dané bázi můžeme snadno přepsat operátorovou rovnici  gf ||ˆA . Vložíme před vektor jednotkový operátor.    zleva|/|ˆ ||ˆˆ k|ll| l gf gf A 1A .|ˆ  gf k|ll||k l A (2.23) Získaný výraz není ale nic jiného než maticové násobení 11 1 1 1 1 ˆ ˆ , , , , n k l n n nn n n n A A f g A k| |l | f n| A A f g | g n|                                                 f f g g          A A neboli .k l lkl gfA  (2.24) Jestliže operátoru přiřadíme čtvercovou matici  l||kA lk AA ˆˆ a vektoru sloupcovou matici  ff |lf| l , můžeme s operátorovými výrazy zacházet jako s obyčejným násobením matic. Hovoříme o tom, že jsme zvolili reprezentaci daného prostoru. Ve skutečnosti nejde o nic jiného než o volbu konkrétní báze. Má-li báze nekonečný spočetný počet členů, budou vektorům odpovídat posloupnosti a operátorům nekonečné matice. Vidíme, že v libovolném separabilním Hilbertově prostoru existuje jednoznačné zobrazení prvků na prostor posloupností l 2 (izomorfismus). V případě neseparabilních prostorů s nespočetnou bází získáme obdobnou relaci: Kvantová teorie Operátory 21 ,|ˆ  gf xdy|yy||x A (2.25) což není nic jiného než integrální transformace.   gf x|xgy|yfy|x|yxA xgdyyfyxA )(,)(,ˆ),( ;)()(),( A (2.26) s jádrem A(x, y). Veličiny x, y hrají úlohu spojitého indexu. Přechod mezi bázemi Máme-li k dispozici dvě sady bázových vektorů }{| k a }{| k , bude mezi koeficienty rozvojů platit jednoduchý vztah, který opět odvodíme jen vložením relace úplnosti: .,|     |kkSfSk|kk|kf kk k kkk k k ff Matice S se nazývá matice přechodu. 2.2.5 Spektrální teorie V teorii operátorů patří k základním úlohám nalézt směry, ve kterých se působení daného operátoru projevuje jako komplexní natahování: .;||ˆ C  ffA (2.27) Vektor f| se nazývá vlastním vektorem (charakteristickým vektorem) operátoru Aˆ a koeficient natahování  vlastním číslem (charakteristickým číslem). Například u tenzoru setrvačnosti leží vlastní vektory ve směru os, ve kterých těleso při rotaci nehází. U lineárních operátorů je i každý násobek vlastního vektoru vlastním vektorem se stejným vlastním číslem. Jde tedy o celý paprsek v Hilbertově prostoru, neboli vlastní směr. Takových vlastních směrů a čísel může existovat u lineárních operátorů celá řada, jejich maximální počet je roven dimenzi prostoru (počtu prvků báze). U separabilních prostorů můžeme tedy úlohu nalezení vlastních čísel a vektorů formulovat rovnicemi: .;,2,1;||ˆ C kkkk k  ffA (2.28) Množina všech vlastních čísel },,,{ 1  n se nazývá spektrum operátoru Aˆ . Nalezneme-li spektrum operátoru a jeho vlastní vektory, můžeme relativně snadno řešit rovnice obsahující tento operátor. Pomocí vlastních čísel a vektorů lze například řešit soustavy obyčejných lineárních diferenciálních rovnic (viz kapitola 1.5 první části sylabu). Věta: Hermitovský operátor má reálná vlastní čísla a vlastní vektory dvou různých vlastních čísel jsou navzájem kolmé. Důkaz: Při výpočtu skalárního součinu využijeme hermitovosti a ponecháme operátor nejprve působit v pravé části skalárního součinu a poté v levé. Výsledek musí být stejný: . ||ˆ || ˆ| * * R        kkk kkkkk kkkkkk kk    ffff ffff ff A A Vlastní čísla jsou tedy reálná. V druhé části budeme postupovat obdobně. Pro vlastní číslo z levé části skalárního součinu využijeme první části důkazu: Kvantová teorie Operátory 22 * | | ˆ| ˆ | | | k l l l k l k l k l k l k k lk                       f f f f f f f f f f f f A A | | ( ) | 0 | 0 pro . l k l k k l l k k l k l k l                    f f f f f f f f Věta: Vlastní čísla unitárního operátoru leží na komplexní jednotkové kružnici a vlastní vektory dvou různých vlastních čísel jsou navzájem kolmé. Důkaz: Při důkazu vyjdeme s definice unitárního operátoru: .|ˆ|ˆ| *  kkkkkkkk ffffff UU Porovnáním prvního a posledního výrazu je zřejmé, že .1||1*  kkk  Zbývá dokázat kolmost vektorů. K tomu budeme potřebovat pomocnou větu (lemma), která říká, že inverzní unitární operátor má vlastní čísla 1/λk. 1 * 1ˆ ˆ ˆ ˆ k k k k k k k k k k k k k k         f f f f f f f f f f f fU U U U . Z porovnání prvního a posledního výrazu plyne 1 1ˆ k k k  f fU . Nyní již budeme postupovat analogicky jako v případě hermitovského operátoru. V levé části skalárního součinu využijeme opět první části důkazu: 1 * | | ˆ 1| ˆ | | | | | ( ) | 0 | 0 pro . k l l l k l k l k l k l k k l k l k l k k l l k k l k l k l                                             f f f f f f f f f f f f f f f f f f f f U U Poznámky: 1. Reálné vlastní hodnoty hermitovských operátorů budou v kvantové teorii využity jako možné výsledky měření dynamické proměnné, které odpovídá operátor Aˆ . Ani kolmost vlastních vektorů není bez užitku. Vhodný hermitovský operátor nám může v podobě svých vlastních vektorů „porodit“ výhodnou ortonormální bázi v Hilbertově prostoru. 2. Unitární operátory se v kvantové teorii využívají k popisu časového vývoje stavu objektu. Kvantová teorie Operátory 23 Příklad 11: Určeme vlastní čísla a vektory matice Aˆ z příkladu 4: , 0i i0ˆ        A . 0 0 i i 0i i0 ||ˆ 2 1 2 1 2 1                                         f f f f f f    ffA Rovnice má netriviální řešení jen pokud je determinant roven nule. Z této podmínky plynou dvě možné hodnoty vlastního čísla . Pro každou z nich potom již snadno učíme příslušný vlastní vektor. Pozor! Podmínka pro determinant činí rovnice pro složky vlastního vektoru závislé. To je ale v pořádku, řešení rovnic musí mít jeden volný parametr, tak aby představovalo celý paprsek v prostoru. K vlastním vektorům můžeme najít normované vlastní vektory a příslušné projekční operátory. Výsledek je: . 1i i1 2 1 |22|ˆ, 1 i 2 1 2|, 1 i |,1 , 1i i1 2 1 |11|ˆ, 1 i 2 1 1|, 1 i |,1 222 111                                                   P P c c f f   Povšimněte si, že matice A byla hermitovská (matice překlopená kolem diagonály a komplexně sdružená je shodná s původní maticí). Proto má reálná vlastní čísla a vlastní vektory tvoří ortonormální soustavu: .01|2,12|2  2|11|1 Tato soustava vektorů je úplná, tvoří bázi v prostoru komplexních dvojic: .1PP ˆˆˆ|22||11| 21  Věta o spektrálním rozvoji: Nechť Aˆ je lineární operátor s množinou vlastních vektorů, která tvoří úplnou ortonormální bázi v Hilbertově prostoru. Potom můžeme pro analytickou funkci operátoru definovanou Taylorovým rozvojem (2.8) psát .ˆ)(||)()ˆ(   k kk k k fkkff PA  (2.29) Důkaz: Nejprve si povšimněme působení mocnin operátoru na vlastní vektory: .|)(||ˆ|)ˆ( ,|ˆ ,ˆ|ˆ ,|ˆ 22       kkk n k n nk n n nk k n kk n kkkkk kkk fccf ffff ff fff ff     AA A AA A | || |  V další fázi důkazu budeme funkcí operátoru působit již na obecný vektor. Provedeme ale jeho rozklad do báze složené z vlastních vektorů, kde působení známe z poslední rovnosti: .||)(||)ˆ(|)ˆ(    fff k k k kkfkkff AA To je ale přesně rovnost, kterou jsme chtěli dokázat. Vynecháme-li ve výrazech libovolný vektor f| , na který operátory působí, dostáváme větu o spektrálním rozvoji. Kvantová teorie Operátory 24 Poznámky: 1) Má-li operátor vícenásobné vlastní číslo stupně N, není to na závadu. Vlastní vektory k vícenásobnému vlastnímu číslu tvoří celý podprostor P dimenze N a lze zvolit N nezávislých kolmých vlastních vektorů odpovídajících tomuto vícenásobnému číslu. 2) Je-li prostor neseparabilní, bude za některých dalších předpokladů možné větu o spektrálním rozvoji modifikovat do tvaru: .||)()ˆ(   dxxxff xA 3) Pro vyjádření funkce operátoru je často mnohem jednodušší použít místo Taylorova rozvoje větu o spektrálním rozvoji. Příslušná řada probíhá jen přes vlastní čísla operátoru. 4) Známe-li spektrum operátoru a jeho vlastní vektory, můžeme snadno napsat jeho libovolnou funkci a řešit tak rovnice, ve kterých se tato funkce operátoru vyskytuje. Speciálně inverzní operátor je dán formulí 1 1 1ˆ ˆ| | .k k kk k k k       A P Vidíme, že pro jeho existenci je kromě předpokladů věty nutná nenulovost všech vlastních čísel. Příklad 12: Na prostoru C 2 je zadán maticový operátor 0 iˆ i 0        A . Určete exp( Aˆ ). Tento příklad jsme již řešili pomocí Taylorova rozvoje jako příklad 4. Z příkladu 11 známe vlastní čísla, vektory i projekční operátory tohoto maticového operátoru. Z věty o spektrálním rozvoji můžeme proto napsat: . 1ch1shi 1shi1ch 1i i1 2 1 e 1i i1 2 1 eˆeˆee 11 21 ˆ 21                          PPA λλ Na rozdíl od Taylorova rozvoje má řada nyní jen dva členy. Příklad 13: Nalezněte časový vývoj průběhu teploty tyče délky L, jejíž oba konce jsou udržovány na nulové teplotě. Počáteční teplota tyče je dána funkcí )(0 xT . Úkolem je najít řešení rovnice teplotní difúze ),(, 2 2 xtTT x T t T        s počátečními a okrajovými podmínkami .0),()0,(,)(),( 00  LtTtTxTxtT Přeformulujme nyní úlohu do Diracovy symboliky. Zaveďme nejprve Hilbertův prostor  0)()0(),0(:)( 2  LffLLfxfH . Jde o prostor funkcí definovaných na intervalu <0, L>, periodických a integrovatelných s kvadrátem. Okrajová podmínka původní rovnice je přesunuta do definice prostoru. Kdyby byly okrajové hodnoty na obou stranách různé, mohli bychom teplotní funkci zrcadlově rozšířit (tím bychom získali stejnou hodnotu teploty na obou koncích tyče). Jestliže by hodnota teploty tyče na obou koncích byla nenulová, můžeme posunout počátek teplotní stupnice. Vzhledem k derivacím v rovnici teplotní difúze to na tvar rovnice nemá vliv. Požadavek nulové teploty na obou koncích tyče tedy není na újmu obecnosti řešení. Úloha má nyní tvar: .ˆ,|,|ˆ| 2 2 dx d TTκ td Td   AA H (2.30) Kvantová teorie Operátory 25 Z příkladu 8 víme, že operátor xd d iˆ B je hermitovský (má reálná vlastní čísla a kolmou soustavu vlastních vektorů). Proto i kvadrát tohoto operátoru 2ˆˆ BA  je hermitovský. Znaménko minus zde není podstatné, jen zajistí nezáporná vlastní čísla operátoru Aˆ . Řešení úlohy můžeme formálně napsat okamžitě (!!) : .)(|e)(| 0 )(ˆ 0   tTtT ttA (2.31) Skutečně: dosadíme-li počáteční čas, dá řešení počáteční podmínku. Derivujeme-li řešení podle času, zjistíme, že řešení (2.31) splňuje výchozí rovnici (2.30): .)(|ˆ)(|eˆ)(|e)(| 0 )(ˆ 0 )(ˆ 00   tTκtTκtT dt d tT dt d ttκttκ AA AA Jediný problém je, že v nalezeném řešení (2.31) vystupuje funkce operátoru Aˆ . Abychom ji mohli určit, musíme znát spektrum a vlastní vektory operátoru Aˆ na prostoru H. Řešme tedy nejprve úlohu .0)()0(,0 |,||ˆ   Lffff fff   HA Řešení této obyčejné lineární diferenciální rovnice s okrajovými podmínkami je: .,2,1,sin 2 |,sin|, 2 22              kx L k L kx L k cf L k kk   Vlastní čísla jsou reálná ( Aˆ je hermitovský operátor), vlastní vektory jsou navzájem kolmé a tvoří přirozenou bázi v H. Napsat řešení (2.31) je nyní již jen jednoduchou aplikací věty o spektrálním rozvoji: ;|e)(||e)(|e)(| 1 )( 1 0 )( 0 )(ˆ 02 22 00        k tt L κk k k ttκλttκ kctTkktTtT k  A .)(|,sin 2 | 0        tTkcx L k L k k  Řešení můžeme samozřejmě zapsat i standardně, bez použití Diracovy symboliky: .)(sin 2 ;sine 2 ),( 0 0 1 )( 02 22                  L k k tt L κk k dxxTx L k L cx L k c L xtT   Pro t = t0 máme Fourierovu řadu pro počáteční podmínku. Pro t  t0 nejde o nic jiného než o rozvoj do jednotlivých Fourierových modů. Každá Fourierova komponenta ubývá exponenciálně s časem. Vidíme, že věta o spektrálním rozvoji nám může být užitečná i při řešení parciálních diferenciálních rovnic. Kvantová teorie Základní principy 26 2.3 ZÁKLADNÍ PRINCIPY KVANTOVÉ TEORIE Jak už víme, klasická mechanika selhala při popisu dějů mikrosvěta zejména proto, že je postavena na komutujících objektech. V mikrosvětě děje ale nekomutují. Základním cílem bude tedy místo dynamických proměnných používat nekomutující objekty (operátory) a nalézt vztah mez operátory a reálnými měřitelnými veličinami. 2.3.1 Základní axiomy a definice I. Redefinice stavu V klasické mechanice je stav částice určen polohou a hybností. Vzhledem k tomu, že v mikrosvěte nelze současně tyto veličiny měřit a měření jedné ovlivní měření druhé, je nutné pojem stavu definovat novým způsobem. Fázové trajektorie již nelze v mikrosvětě popsat křivkami. Vnímáme je s přesností danou relacemi neurčitosti 2 px . Můžeme si představit, že fázovou trajektorii vidíme jako rozmazanou čáru s rozlišením daným čtverečkem o ploše 2 (sledujeme-li jednu souřadnici a jí odpovídající zobecněnou hybnost). Kompatibilita: Řekneme, že dvě dynamické proměnné jsou kompatibilní, jestliže měření jedné veličiny neovlivní měření veličiny druhé. Příkladem kompatibilních proměnných jsou (x, y), příkladem nekompatibilních proměnných jsou (x, px). Kompatibilita je symetrická vlastnost, ale není transitivní. (A komp B)  (B komp A) , (A komp B)  (B komp C) / (A komp C) . Příklad 14: (x komp y)  (y komp px) / (x komp px). Úplná množina pozorovatelných (ÚMP): Maximální nezávislá množina vzájemně kompatibilních dynamických proměnných. Jakákoli další dynamická proměnná už s nimi není kompatibilní. Například v nerelativistické teorii jsou nejznámější ÚMP (x, y, z), (px, py, pz) a pro centrální pole ještě energie, kvadrát momentu hybnosti a jeho jedna komponenta (E, L2 , L3). Současně lze tedy změřit všechny tři souřadnice nebo všechny tři složky hybností. Nelze již ale současně změřit všechny tři složky momentu hybnosti. Stav systému: Řekneme, že známe stav systému, známe-li výsledek měření některé ÚMP. Stavem tedy nazveme jen to, co lze ve skutečnosti současně změřit. Základním rysem nové teorie musí být nekomutující objekty – operátory. Místo dynamických proměnných z klasické mechaniky (souřadnice, hybnost, energie, moment hybnosti,…) budeme používat operátory (operátor souřadnice, hybnosti, ...). Nekomutativnost těchto operátorů bude vyjadřovat nekomutativnost aktu měření dynamických proměnných v mikrosvětě. Veličiny naměřené přístrojem v mikrosvětě jsou reálná čísla, někdy spojitá (poloha částice), někdy diskrétní (například jednotlivé hladiny energie elektronu vázaného v atomu). Jak získat z operátoru dynamické proměnné sadu reálných čísel spojitého nebo diskrétního charakteru? Takovou sadou je právě spektrum hermitovských operátorů. Dynamickým proměnným budeme tedy přiřazovat hermitovské operátory. Každý operátor působí na prvky nějakého Hilbertova prostoru H. Musíme se tedy ptát, jaký význam bude v naší teorii mít sám Hilbertův prostor a také prvky, na které operátory působí. Později uvidíme, že příliš nezáleží na volbě Hilbertova prostoru. Podstatné jsou spíše vztahy mezi dynamickými proměnnými, nyní operátory. Kvantová mechanika založená na prostoru L 2 funkcí integrovatelných s kvadrátem je známá Schrödingerova vlnová mechanika x p x p Kvantová teorie Základní principy 27 vedoucí na Schrödingerovu rovnici a vlnové funkce. Kvantová teorie založená na prostoru l 2 posloupností sčitatelných s kvadrátem je Heisenbergova maticová mechanika. Obě teorie se na první pohled zdají naprosto odlišné. Přesto vlastní čísla operátorů v obou teoriích jsou stejná a obě teorie tak dávají stejné předpovědi. Hilbertův prostor se všemi svými prvky a s operátory, které na prvky působí, koresponduje s vlastnostmi celého systému z klasické mechaniky. Místo systému budeme proto v kvantové teorii hovořit o Hilbertově prostoru daného systému (například Hilbertův prostor elektronu). Zbývá rozluštit poslední hádanku – k čemu jsou prvky Hilbertova prostoru? Již v úvodu jsme si řekli, že v mikrosvětě sám akt měření ovlivní stav systému. Před měřením je systém v jiném stavu než po měření. Akt měření dynamické proměnné zastupuje v kvantové teorii hermitovský operátor této proměnné. Působením tohoto operátoru na prvek prostoru dostáváme jiný prvek tohoto prostoru. A to je přesně to, co hledáme. Prvky (vektory) prostoru tedy představují stav systému. Akt měření je působení příslušného operátoru na stav (prvek prostoru) a nový stav je prvek, který vznikl působením operátoru. Vlastní číslo operátoru prezentuje naměřenou hodnotu a vypovídá tak o stavu systému. Víme už, že násobky každého vlastního vektoru jsou opět vlastním vektorem. V H tedy existuje k danému vlastnímu číslu celý vlastní směr (paprsek). Stavu systému proto musí odpovídat celý paprsek v H, nikoli jen jeden jediný vektor. Přicházíme tak ke třem základním axiomům kvantové teorie:  systému přiřadíme Hilbertův prostor H systém  H  stavu systému přiřadíme paprsek |  > v H stav  |  dynamické proměnné přiřadíme hermitovský operátor A na H A  Aˆ S linearitou budované teorie se pojí velmi důležitý princip:  Princip superpozice: Nechť H| a H| reprezentují dva různé stavy systému. Potom je každý vektor  |ψα|α 21  je také fyzikálně realizovatelný stav. Bez tohoto požadavku by nebylo možné vybudovat lineární teorii. II. Měření v kvantové teorii Akt měření dynamické proměnné A v nějakém stavu znamená aplikaci operátoru Aˆ této dynamické proměnné na daný stav | . Operátorem Aˆ a stavem | musí tedy být zcela jednoznačně dáno, co je a co není možné na systému naměřit. Odpověď na tuto otázku poskytují tzv. interpretační postuláty:  Postulát A. Měříme-li dynamickou proměnnou A, mohu na systému naměřit jen některou z vlastních hodnot {aj} operátoru Aˆ této dynamické proměnné:  j|aj| jAˆ .  Postulát B. Pozorování dynamické proměnné A na systému, který byl připraven ve vlastním stavu j| operátoru Aˆ , vede zcela jistě k naměření vlastní hodnoty aj.  Postulát C. Je-li systém připraven v obecném stavu H| , vede opakované měření veličiny A k různým výsledkům aj. Střední hodnota těchto opakovaných měření bude rovna   |ˆ| A . Kvantová teorie Základní principy 28 Poznámky: 1) Opakovaná měření si nemusíme představovat tak, že bychom na stejném systému opakovali neustále tatáž měření. V praxi by to nebylo proveditelné. Těžko můžeme na jednom jediném elektronu zopakovat nějaké měření. Musíme mít připraveno velké množství systémů ve stejném stavu (například svazek elektronů) a opakovat měření na různých elektronech (systémech). 2) Výraz pro střední hodnotu je nejjednodušším možným výrazem složeným z operátoru Aˆ a stavu | , který dá jako výsledek reálné číslo. Střední hodnotu bývá zvykem označovat nebo A . 3) Automaticky předpokládáme, že stavové vektory jsou normovány k jedné. Není-li stavový vektor normován, musíme výraz pro střední hodnotu vydělit ještě kvadrátem normy stavového vektoru: . | |ˆ|      A A 4) Výraz pro střední hodnotu rozepsaný v prostoru L2 (R 3 ) dá: . )()( )(ˆ)( 3* 3*   xxx xxx d d A   A 5) Jsou-li všechny systémy připraveny ve vlastním stavu j| operátoru Aˆ , dá střední hodnota daná postulátem C samozřejmě příslušnou vlastní hodnotu podle postulátu B a všechna měření dají v tomto výjimečném případě stejný výsledek (přes j se nesčítá): . | | | |ˆ| j j a jj jja jj jj A        A III. Statistická interpretace stavového vektoru Rozvineme-li stavový vektor do ortonormální spočetné báze | n  nebo nespočetné báze | x  | | | |nn n n        , resp. | | | ( )|x x dx x x dx        , nazýváme koeficienty rozvoje n , respektive ( )x , amplitudou pravděpodobnosti, že systém nalezneme ve stavu | n , respektive | x . K tomu nás opravňuje fakt, že jde o projekce stavového vektoru do patřičného prvku báze. Z ortonormality báze a normovanosti stavu k jedné okamžitě plyne * 1n n   , resp. * ( ) ( ) 1x x dx   a výrazy * n n nw   , resp. * ( ) ( ) ( )w x x x  proto chápeme jako pravděpodobnost realizace stavu | n  resp. hustotu pravděpodobnosti nalezení systému ve stavu | x  . Pravděpodobnosti jsou automaticky normalizovány k jedné. IV. Princip korespondence Posledním ze základních principů kvantové teorie je princip korespondence. Vymezuje, které části z teoretické mechaniky je možné převzít v kvantové teorii.  Princip korespondence pro základní relace. Základní relace mezi dynamickými proměnnými v teoretické mechanice a příslušnými operátory v kvantové mechanice se mohou lišit jen pořadím operátorů.  Princip korespondence pro algebru Poissonových závorek. Struktura Poissonových závorek v teoretické mechanice je identická se strukturou komutátorů v kvantové teorii: Kvantová teorie Základní principy 29 .] CBA C B A ˆˆ,ˆ[},{ ˆ ˆ ˆ kCBA C B A     První část principu korespondence platí pro jednoduché relace mezi dynamickými proměnnými, které neobsahují derivace. Například definice Hamiltonovy funkce v potenciálním poli V ),,( 2 222 zyxV m ppp H zyx    přejde v definici Hamiltonova operátoru .)ˆ,ˆ,ˆ( 2 ˆˆˆ ˆ 222 ZYX PPP H V m zyx    (2.32) Výraz pro potenciální energii je typickou funkcí operátoru, kterou jsme probírali v kapitole o operátorech. Pro výrazy typu A = xp nelze kvantový analog jednoznačně určit. Může jím být buď PXA ˆˆˆ  nebo XPA ˆˆˆ  . Operátory nekomutují, a proto záleží na jejich pořadí. Správná varianta z obou možných musí být vybrána na základě experimentu. Stejně tak můžeme z různých Lagrangeových funkcí téhož systému obdržet různé kvantové teorie a správnou variantu je třeba opět vybrat na základě experimentu. Druhá část principu korespondence se týká Poissonových závorek – výrazů, které v klasické mechanice obsahují derivace dynamických proměnných. Poissonovým závorkám v kvantové teorii odpovídají komutátory dynamických proměnných. Nelze však položit rovnost mezi komutační relací a Poissonovou závorkou. Důvody jsou hned dva: 1) rozměrový: Poissonova závorka obsahuje derivace, které do výrazů vnášejí fyzikální rozměr, komutátory nikoli  je třeba použít rozměrový převodní koeficient k. 2) principiální: Dynamickým proměnným v kvantové teorii můžeme přiřazovat jen hermitovské operátory (mají reálná vlastní čísla, která interpretujeme jako naměřitelné hodnoty). Jsou-li operátory odpovídající A a B hermitovské, musí být operátor odpovídající C také hermitovský. To lze opět zajistit pomocí konstanty k. Určeme nyní podmínku na konstantu k, která plyne z požadavku hermitovosti operátorů: † † † † † * † † * * ˆˆ ˆ[ , ] ˆˆ ˆ ˆ ˆ / ˆ ˆ ˆˆ ˆ ˆ ˆ / ˆˆ ˆ ˆ ˆ ˆˆ ˆ[ , ] k k k k k .           A B C AB BA C B A A B C O O BA AB C A B C Porovnáme-li počátečný a koncový výraz, musí platit k* = – k. To ale splňují jen ryze imaginární čísla. Převodní konstanta tedy musí mít tvar: i .k   (2.33) Konstanta  je nějaké reálné číslo a je jedinou fundamentální konstantou kvantové teorie. Tato konstanta se bude vyskytovat ve všech předpovědích kvantové teorie (například v energetickém spektru elektronu vázaného v atomech, ve vztazích pro záření absolutně černého tělesa, v Heisenbergových relacích neurčitosti, atd.). Její hodnotu je možné změřit experimentálně na základě těchto předpovědí a je rovna: .100545887.1 34 sJ  (2.34) Kvantová teorie Základní principy 30 Princip korespondence pro Poissonovy závorky můžeme stručně zapsat jako .]BA ˆ,ˆ[ i 1 },{  BA (2.35) Tím jsme zakončili přehled základních principů kvantové teorie. Jelikož jde o základní neodvoditelné principy, na kterých teorii stavíme, bylo by možné jen stroze vypsat axiomy, postuláty a principy označené v této kapitole čtverečkem. Doplňující texty se snaží jen poukázat na to, že právě tato volba základních axiomů je přirozená a povede k cíli. O správnosti základních principů však mohou rozhodnout jedině experimenty ověřující výpovědi z těchto principů plynoucí. 2.3.2 Kompatibilita měření a Heisenbergovy relace Rozhodnout o tom, zda se měření dvou dynamických proměnných ovlivňují či nikoli, je jednoduché. Stačí znát komutátor operátorů těchto proměnných. Je-li tento komutátor nulový, je ABBA ˆˆˆˆ  a měření se neovlivňují. Základní komutátory pro souřadnice a hybnosti můžeme odvodit z principu korespondence, ostatní už pak z vlastností komutátorů. Pro Poissonovy závorky mezi souřadnicemi a hybnostmi platí (1.55): .},{,0},{},{ kllklklk pxppxx  Tomu odpovídají podle principu korespondence komutační relace: ˆ ˆ ˆ ˆ ˆ ˆ ˆ[ , ] [ , ] 0 , [ , ] i .k l k l k l kl   X X P P X P 1 (2.36) Z nich je zřejmé, že současně lze u objektu změřit všechny tři souřadnice nebo hybnosti. Také se vzájemně neovlivní měření například souřadnice x a hybnosti py. Jediná měření, která se vzájemně ovlivňují a u kterých záleží na pořadí měření (nenulový komutátor) je měření zobecněné souřadnice a jí příslušné zobecněné hybnosti. Rovnice (2.36) jsou základními komutačními relacemi v kvantové teorii. Bylo by samozřejmě možné hledat ostatní složitější komutační relace také z Poissonových závorek. Výhodnější je odvozovat je ze základních relací (2.36) a vlastností Lieovy algebry komutátorů. Tím se oprostíme od klasické mechaniky a nemusíme se k ní při každé komutační relaci vracet. Kvantová mechanika začíná „žít vlastním životem“. To, co přebrala z klasické mechaniky prostřednictvím principu korespondence, jsou jen relace (2.36). Odvoďme nyní komutační relaci mezi první a druhou komponentou momentu hybnosti:     2313312321331323 3312313231321332 3123132331321332 3113233221 ˆ]ˆˆ,ˆ[]ˆˆ,ˆ[ˆˆ]ˆˆ,ˆ[]ˆˆ,ˆ[ˆ ˆ]ˆˆ,ˆ[]ˆˆ,ˆ[ˆˆ]ˆˆ,ˆ[]ˆˆ,ˆ[ˆ ]ˆˆ,ˆˆ[]ˆˆ,ˆˆ[]ˆˆ,ˆˆ[]ˆˆ,ˆˆ[ ]ˆˆˆˆ,ˆˆˆˆ[]ˆ,ˆ[ PPXXPXPXPPXXPXPX PPXXPXPXPPXXPXPX PXPXPXPXPXPXPXPX PXPXPXPXLL 2 3 3 1 2 3 3 1 3 2 1 3 2 3 1 3 2 1 3 3 2 3 1 3 1 2 3 3 2 1 3 3 3 3 2 1 3 2 3 1 3 3 1 2 3 3 1 2 3 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ[ , ] [ , ] [ , ] [ , ] ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ[ , ] [ , ] [ , ] [ , ] ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ[ , ] [ , ] [ , ] [ , ] ˆ ˆ                 X X P P X P X P X X P P X X PP X X P P X P X P X X P P X X P P X X P P X P X P X X P P X X PP X 1 2 3 3 2 1 3 1 3 3 2 3 1 3 2 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ[ , ] [ , ] [ , ] [ , ]   X P P X P X P X X P P X X P P .312212112 23311332 ˆi)ˆˆˆˆ(iˆˆˆiˆˆˆi 0ˆ]ˆ,ˆ[ˆ000000000000ˆ]ˆ,ˆ[ˆ0 LPXPXP1XP1X PPXXPPXX    Kvantová teorie Základní principy 31 Je jasné, že postup je velmi zdlouhavý, ale přímočarý. Hledanou komutační relaci postupně „rozmělňujeme“ podle pravidel Lieovy algebry až na elementární relace mezi souřadnicemi a hybnostmi. Prakticky všechny symbolicky orientované programy či programovací jazyky bez problémů tuto úlohu řeší za nás a obsahují balíky pro výpočet komutačních relací. Analogickým postupem můžeme nalézt komutační relace pro ostatní komponenty momentu hybnosti. Není to ale nutné, stačí je získat cyklickou záměnou souřadnicových os. Kompletní komutační relace pro moment hybnosti potom jsou: .ˆi]ˆ,ˆ[ ,ˆi]ˆ,ˆ[ ,ˆi]ˆ,ˆ[ 213 132 321 LLL LLL LLL       (2.37) Výsledkem je, že současně není možné změřit žádné dvě komponenty momentu hybnosti. Měření každé komponenty ovlivní měření kterékoli jiné komponenty. Zaveďme operátor kvadrátu velikosti momentu hybnosti .ˆˆˆˆ 2 3 2 2 2 1 2 LLLL  (2.38) Stejným postupem jako dříve dopočteme komutační relace kvadrátu momentu s jednotlivými komponentami. Tentokráte při „rozmělňování“ komutační relace postačí dostat se jen k relacím (2.37) pro moment hybnosti. Jejich výsledek už známe. Po výpočtu dostaneme: .3,2,1,0]ˆ,ˆ[ 2  kkLL (2.39) Není tedy možné současně změřit dvě různé komponenty momentu hybnosti. Vždy je ale možné změřit kvadrát velikosti momentu hybnosti a jednu z jeho libovolných komponent, zpravidla se používá třetí komponenta. Ze zatím provedených úvah je zřejmé, že současně můžeme měřit dynamické proměnné {x, y, z} nebo {px, py, pz} nebo {L2 , L3}. V kapitole 2.5 uvidíme, že v případě sféricky symetrického potenciálu je s poslední množinou kompatibilní ještě energie. Jde o základní tři úplné množiny pozorovatelných (ÚMP) pro nerelativistickou částici. Nalezli jsme tedy jednoduchý postup, pomocí kterého zjistíme, které veličiny lze současně měřit a které ne. Postačí nalézt komutátor odpovídajících operátorů. Tento postup nám ale umožní odpověď typu ano/ne. V případě, že dynamické proměnné spolu současně měřit nelze, se musíme ptát, jak moc naruší akt měření jedné proměnné akt měření proměnné druhé. Na tuto otázku odpovídají Heisenbergovy relace neurčitosti, které si nyní odvodíme. Předtím si uveďme přehled základních statistických pojmů a jejich operátorových analogií v kvantové teorii: statistika kvantová teorie a střední hodnota   |ˆ| Aa střední hodnota aaa  odchylka 1AA ˆˆˆ a operátor odchylky 0a vztah pro a 0|ˆ|   A vztah pro a Kvantová teorie Základní principy 32   222 aaa  variance 2 2 2 |ˆ| |ˆ| |ˆ|       A A A vztah pro  2 a (variance)   222 aaaakv  střední kvadr. odchylka   |ˆ| 2 Akva střední kvadr. odchylka Zkuste si dokázat oba dva statistické vztahy v klasické statistice i v kvantové teorii. V obou případech stačí jen dosadit z příslušných definic. Nyní již můžeme přistoupit k odvození relací neurčitosti. Předpokládejme, že máme dvě nekompatibilní proměnné: .ˆ]ˆ,ˆ[;ˆ,ˆ, CBABA BA Nalezněme střední kvadratické chyby měření:     (*1) (*2) (*3) 2 2 2 2 2 2 2 2 ˆ ˆ ˆ ˆ ˆ ˆ| ( ) | | ( ) | | | ˆ ˆ ˆ ˆ ˆ ˆ|| || || || | | | | | | kv kva b |                                               A B A A B B A B A B A B (*4) (*5) 2 2 2 1 1ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ| | ( ) ( ) | 2 2 1 1ˆ ˆ ˆ ˆ ˆ ˆ| | ( ) | |[ ] | 2 2 | | , |                                  A B B A A B B A A B B A A B 2 21 1 ˆˆ ˆ| |[ ] | | | | . 2 2 , | |        A B C Po odmocnění dostáváme konečný tvar Heisenbergových relací: .|ˆ|| 2 1   |ba kvkv C (2.40) * Poznámky k odvození: (1) využití hermitovosti operátorů; (2) Schwartzovo lemma (2.7); (3) rozdělení na symetrickou (S) a antisymetrickou (D) část; (4) symetrická část S je reálná (je součtem dvou navzájem komplexně sdružených čísel), antisymetrická část D je naopak ryze imaginární (je rozdílem dvou navzájem komplexně sdružených čísel) a tvoří dohromady komplexní číslo, pro které platí |S+D| = |x+iy| = (x2 +y2 )1/2 ≥ |y| = |D|. (5) jednotkový operátor v definici Aˆ komutuje s čímkoli. Známe-li výsledek komutační relace operátorů příslušících dvěma dynamickým proměnným, můžeme z Heisenbergových relací určit míru ovlivnění jednoho měření druhým. Toto vzájemné ovlivnění výsledků měření závisí na stavu, ve kterém je systém připraven. Jen jsouli obě dynamické proměnné ve vztahu zobecněná souřadnice – zobecněná hybnost, nezávisí vzájemné ovlivnění na stavu systému: . 2 ||| 2 |ˆ|| 2 ˆi]ˆ,ˆ[     |px 11PX To je nejznámější tvar relací neurčitosti . 2   px (2.41) Kvantová teorie Základní principy 33 Jde o první konkrétní měřitelný výsledek z námi budované teorie, který obsahuje jedinou konstantu teorie – Planckovu konstantu. Poznámka: Operátory kinetické a potenciální energie zpravidla vzájemně nekomutují. To má za následek, že není současně přesně zjistitelná kinetická i potenciální energie a částice se může (na rozdíl od klasické fyziky) „přehoupnout“ přes potenciálovou bariéru (tzv. tunelový jev). 2.3.3 Vlastní stavy energie, Schrödingerova rovnice V minulé kapitole jsme se naučili rozhodnout, které dynamické proměnné lze společně měřit a které ne. Pomocí Heisenbergových relací neurčitosti jsme schopni i kvalitativně postihnout míru narušení jednoho měření měřením druhým. Nyní se budeme věnovat druhé základní úloze kvantové mechaniky: Nalézt spektrum operátoru energie – hodnoty energie, které je možné na systému naměřit. Úlohu můžeme zformulovat například takto: ,||ˆ  nEn nH (2.42) ,)ˆ,ˆ,ˆ( 2 ˆˆˆ ) ˆ ( 2 ˆ ˆ 2222 ZYX PPP X P H V m V m zyx      (2.43) ˆ ˆ ˆ ˆ ˆ ˆ ˆ, , 0 , , i .k l k l k l kl              X X P P X P 1 (2.44) Budeme hledat vlastní hodnoty operátoru energie (Hamiltonova operátoru) ze vztahu (2.42). Tato rovnice pro vlastní hodnoty Hamiltonova operátoru se nazývá Schrödingerova rovnice. Operátor energie (Hamiltonův operátor) je dán vztahem (2.43). Základní operátory, ze kterých je složen Hamiltonův operátor, podléhají komutačním relacím (2.36), resp. (2.44). Nezáleží příliš na tom, jaký Hilbertův prostor zvolíme. V příští kapitole uvidíme řešení harmonického oscilátoru v různých prostorech H, vždy dostaneme stejné spektrum Hamiltonova operátoru. Na daném prostoru je nejpodstatnější zvolit Hermitovy operátory zobecněných souřadnic a hybností tak, aby splňovaly komutační relace(2.44). Ukažme si nyní přepis Schrödingerovy rovnice v prostoru L2 (R3 ) funkcí integrovatelných s kvadrátem na celém prostoru R 3 . Nejjednodušším operátorem na tomto prostoru je operátor násobení souřadnicí. Tento operátor ztotožníme s operátorem souřadnice: .ˆ;ˆ;ˆ zyx  ZYX (2.45) Nyní zbývá nalézt hermitovské operátory pro hybnost tak, aby splňovaly komutační relace (2.44). V příkladu 6 jsme si ukázali, že operátor derivace a operátor násobení souřadnicí splňují v jedné dimenzi komutační relaci 1XD ˆ]ˆ,ˆ[  , neboli 1DX ˆ]ˆ,ˆ[  . Je zřejmé, že operátor dxdiˆ P splňuje relaci 1PX ˆi]ˆ,ˆ[  a v jedné dimenzi plní úlohu operátoru hybnosti. Sám operátor derivace není hermitovský, ale operátor derivace vynásobený ryze imaginární konstantou již hermitovský je (viz příklad 8). Aby ve třech dimenzích operátor hybnosti splňoval komutační relace (2.44) a platila volba (2.45) pro operátor souřadnice, musí operátor hybnosti mít tvar ,iˆ;iˆ;iˆ zyx zyx           PPP neboli (2.46) .i ˆ     P (2.47) Schrödingerova rovnice (2.42) s operátorem energie (2.43), volbou prostoru H = L2 (R3 ) a operátory (2.45) a (2.46) vede potom na slavnou Schrödingerovu rovnici v x reprezentaci: Kvantová teorie Základní principy 34 .)()()( 2 2 xxx nnn EV m          (2.48) Řešení této rovnice pro konkrétní potenciál V poskytne spektrum operátoru energie {En} jakožto množinu možných měřitelných hodnot energie. Poznámka 1: Řešení rovnice (2.48) lze nalézt pro každou hodnotu energie. Ne vždy je však toto řešení z prostoru L2 (R3 ). Je proto vždy třeba vybrat z možných řešení jen ta, která jsou integrovatelná s kvadrátem, tj. donekonečna ubývají dostatečně rychle, aby zajistila integrovatelnost. Poznámka 2: Existuje jednoduchý způsob, jak odhadnout typ spektra pro daný potenciál. Může-li se v klasické mechanice částice vzdálit do nekonečna, je spektrum operátoru energie spojité. Nemůže-li se ani v jednom směru vzdálit do nekonečna je spektrum operátoru energie diskrétní. V V0 V1 x Připomeňme, že v klasické mechanice se částice může pohybovat tam, kde je její celková energie větší než potenciální. To plyne ze vztahu )(2/2 xVmE  v . Jde vlastně o podmínku nezápornosti kinetické energie. V nakreslené situaci je pro E < V0 spektrum energie diskrétní, pro E > V0 spojité. Příklad 15: x r r r r x x x V V0 V0 V0 V0 V1 V0 V V V V V V V1 2 3 4 5 6 7 8 1D potenciály; ),( x 1. Symetrická pravoúhlá jáma. Diskrétní spektrum pro E < V0, spojité pro E > V0. Pro nekonečnou jámu ( 0V ) je spektrum jen diskrétní. 2. Bariéra. Spojité spektrum. 3. Nesymetrická pravoúhlá jáma. Diskrétní spektrum pro E < V0, spojité pro E > V0. 4. Harmonický oscilátor. Diskrétní spektrum. Kvantová teorie Základní principy 35 3D potenciály; r > 0 5. Sférická jáma. Diskrétní spektrum pro E < V0, spojité pro E > V0. V podobném potenciálu se pohybuje například neutron v atomovém jádře. 6. Coulombova bariéra (sférická jáma kombinovaná s Coulombovým odpuzováním). Diskrétní spektrum pro E < V0, spojité pro E > V0. Nenulová pravděpodobnost průniku bariérou‚ tunelový jev (operátor kinetické a potenciální energie nekomutuje). V podobném potenciálu se pohybuje proton nebo  částice v atomovém jádře. 7. Coulombův přitažlivý potenciál. Diskrétní spektrum pro E < 0, spojité pro E > 0. V podobném potenciálu se pohybuje elektron v atomárním obalu. Stavy se zápornou energií jsou vázané stavy, stavy s kladnou energií jsou volné (elektron není vázán k atomovému jádru). 8. Sférický harmonický oscilátor. Diskrétní stavy energie. Systém je při vychýlení do kteréhokoli směru vracen do počátku podle předpisu V(r) = 1/2 kr2 . Celá kapitola 2.4 bude věnována řešení Schrödingerovy rovnice pro vlastní hodnoty operátoru energie na příkladu harmonického oscilátoru. Harmonický oscilátor je velice důležitým systémem. Jedná se o parabolické přiblížení k jakékoli potenciální energii s minimem a mnoho systémů lze považovat za harmonický oscilátor alespoň v prvním přiblížení. Nakonec i samo elektromagnetické pole je soustavou fotonů – elementárních harmonických oscilátorů. Mohli jsme jistě zvolit příklad jednodušší – pravoúhlou jámu či bariéru. V těchto příkladech je potenciální energie po částech konstantní a řešení Schrödingerovy rovnice je víceméně triviální. Jen je třeba jednotlivé části řešení v místech změny předpisu potenciální energie správně navázat. Tyto elementární úlohy jsou řešeny v každé úvodní učebnici základního kursu fyziky a zvídavý čtenář si je tam jistě najde. Kvantová teorie Harmonický oscilátor 36 2.4 HARMONICKÝ OSCILÁTOR Na příkladu harmonického oscilátoru si ukážeme typické řešení úlohy pro vlastní hodnoty operátoru energie. Naše úloha je   .ˆiˆ,ˆ ,ˆ 2 1 2 ˆ ˆ ,||ˆ 22 2 1PX X P H H    m m nEn n (2.49) Jde o problém vlastních hodnot Hamiltonova operátoru s konkrétním průběhem potenciální energie a zadanými základními komutačními relacemi mezi operátorem polohy a hybnosti. V kapitole 2.4.1 úlohu vyřešíme v rámci klasické Schrödingerovy vlnové mechaniky. Za Hilbertův prostor zvolíme prostor L2 (R ), volba operátorů (2.45) a (2.46) povede na diferenciální rovnici (2.48) v jedné dimenzi. Řešení této rovnice se provádí rozvojem do nekonečných řad, které je třeba „oříznout“ tak, aby řešení bylo z prostoru L2 (R ), tj. integrovatelné s kvadrátem. Odsud získáme spektrum operátoru energie. V kapitole 2.4.2 si ukážeme řešení úlohy (2.49) bez volby reprezentace. Nebudeme vůbec volit konkrétní podobu Hilbertova prostoru. Řešení nalezneme jen z formulace úlohy (2.49). Uvidíme tak, že konkrétní volba Hilbertova prostoru není podstatná. Při tomto přístupu si zavedeme kreační a anihilační operátory, které svým působením posouvají energetické hladiny o jednu výše či níže. Tyto operátory jsou v kvantové teorii velmi užitečné, a proto se s nimi seznámíme již nyní u jednoduchého příkladu harmonických oscilací. V kapitole 2.4.3 si ukážeme řešení úlohy (2.49) na prostoru l 2 nekonečných posloupností sčitatelných s kvadrátem (tzv. Heisenbergova maticová mechanika). Operátory zde budou nekonečné matice. Možná se vám zdá obtížná úloha najít vlastní čísla nekonečných matic. Problém ale není tak složitý. Jestliže za vektory báze zvolíme vlastní vektory příslušného operátoru, bude matice odpovídající tomuto operátoru diagonální. Vlastní čísla diagonálních matic se hledají snadno  . Jsou to právě prvky na diagonále. Třemi různými způsoby tak uvidíte řešení jednoho a téhož problému. V kvantové teorii jde totiž o vnitřní strukturu teorie, nikoli o konkrétní reprezentaci, ve které výpočet provádíme. 2.4.1 Řešení pomocí vlnové mechaniky (Schrödinger) Hamiltonova funkce jednodimenzionálního harmonického oscilátoru je dána součtem kinetické a potenciální energie (1.44) . 2 1 2 ),( 22 2 xm m p pxH  (2.50) Hamiltonův operátor je v prostoru L2 (–, +) potom dán jednoduchou relací . 2 1 2 ˆ 22 2 22 xm dx d m   H (2.51) Odpovídající Schrödingerova rovnice pro vlastní funkci )(x z prostoru L2 (–, +) má tvar Kvantová teorie Harmonický oscilátor 37 )()( 2 1 2 22 2 22 xExxm dx d m             (2.52) Jde o obyčejnou lineární diferenciální rovnici druhého řádu s nelineárním koeficientem u nulté derivace. Standardní tvar této rovnice (s jednotkovým koeficientem u nejvyšší derivace) je: 2 2 2 2 2 2 2 2 0 . d mE m x dx              (2.53) Rovnici budeme řešit ve čtyřech krocích: 1. substituce ve vnitřní (nezávislé) proměnné V nezávislé proměnné budeme volit takovou substituci, která „zbezrozměrní“ rovnici. Přesuňme koeficienty tak, aby byly symetrické u proměnné x 0 22 2 2          E x m dx m d (2.54) a proveďme substituci ,x m     (2.55) po které Schrödingerova rovnice získá bezrozměrný tvar 2 2 2 2 0 , ( ) . d E E d             (2.56) 2. substituce ve vnější (závislé) proměnné V závislé proměnné budeme volit takovou substituci, které zohlední chování vlnové funkce pro   . Pro velká  můžeme zanedbat poslední člen v rovnici (2.56) a přibližně platí .e0 22 2 2 2        d d (řešení stačí dosadit do původní rovnice a zanedbat členy s nižšími mocninami ). Kladné z nalezených řešení evidentně není z prostoru L2 , integrál z kvadrátu přes celý prostor by byl nekonečný. Vlnová funkce se tedy pro velká  musí chovat jako ]2exp[ 2  . To nás přivádí k substituci pro závislou proměnnou 2 /2 ( ) e ( ) ,u     (2.57) po jejímž provedení dostaneme rovnici .0)1(2  uuu  (2.58) Derivace se automaticky rozumí podle nové proměnné . V principu by z matematického hlediska bylo v pořádku říci „v rovnici (2.53) provedeme substituce (2.55) a (2.57), výsledná rovnice je (2.58)“. V bodech 1 a 2 jsme si jen ukázali, jaké pohnutky nás k těmto substitucím vedou, protože postup je obdobný i u jiných průběhů potenciálu. 3. rozvoj řešení do mocninné řady Řešení rovnice (2.58) budeme hledat ve tvaru mocninné řady Kvantová teorie Harmonický oscilátor 38 .)( 0     k k kcu  Snadno nalezneme první a druhou derivaci .)1()(;)( 0 2 0 1         k k k k k k ckkucku  Výrazy pro u a její derivace dosadíme do rovnice (2.58): .0)1(2)1( 000 2          k k k k k k k k k cckckk  Jednotlivé členy upravíme tak, aby mocniny  byly stejné (v prvním členu položíme k–2 = l): .0)1(2)2)(1( 002 2          l l l l l l l l l cclcll  První dva členy prvního součtu jsou nulové, a proto můžeme spodní hranici posunout na l = 0:   .0)12()2)(1( 0 2     l l ll clcll  Má-li být polynomiální výraz identicky nulový pro každou hodnotu , musí být nulové všechny koeficienty, tj. výrazy v hranaté závorce. Získáváme tak rekurentní relaci pro koeficienty cl naší řady: . )2)(1( )12( 2 ll c ll l c     (2.59) Budeme-li znát koeficienty c0 a c1, budeme znát celé řešení, protože z rekurentní relace můžeme spočítat .,,, ,,, 7531 6420   cccc cccc   Koeficienty c0 a c1 tak hrají roli dvou integračních konstant řešení diferenciální rovnice (2.58) druhého řádu. Sudá část řady se počítá z c0 a lichá z c1. 4. oříznutí řady Nalezené řešení je ve tvaru nekonečné mocninné řady. Řeší sice původní rovnici, ale není z prostoru L2 . Aby bylo řešení z L2 (integrovatelné s kvadrátem), musí být řada konečná, tedy polynomiální. Prakticky to znamená, že koeficienty řady musí být od určitého l = n nulové. V rekurentní relaci (2.59) bude čitatel pro toto l = n nulový a veškeré odvozené koeficienty cl s l  n nulové. Vidíme, že nebude možné takto „oříznout“ současně sudé i liché členy řady. Proto jsou možná jen sudá (c0  0, c1 = 0) nebo jen lichá řešení (c0 = 0, c1  0) představující sudý nebo lichý polynom stupně n. Podmínka oříznutí (nulovost čitatele) v (2.59) je 2n +1−  = 0 a plyne z ní po vyjádření spektrum energie harmonického oscilátoru: .)21(  nEn (2.60) Poznámky: 1) Nezapomínejte, že energie E (vlastní hodnota operátoru Hˆ ) je po celou dobu výpočtu schována v bezrozměrné konstantě (vlastním číslu) . 2) Sama Schrödingerova rovnice má řešení pro každou hodnotu energie. Tato řešení ale nejsou integrovatelná s kvadrátem, až výběr integrovatelných funkcí (oříznutí řady) vede k diskrétnímu spektru operátoru energie (jen pro Kvantová teorie Harmonický oscilátor 39 některé vybrané hodnoty energie ubývá řešení v  dostatečně rychle, aby bylo integrovatelné s kvadrátem). Tato situace je typická pro spojité průběhy potenciální energie s minimem. 3) Základní hladina energie 20 E je nenulová! Ani při nulové absolutní teplotě není harmonický oscilátor v klidu a vykonává tzv. nulové kmity (například oscilace krystalové mříže). Při absolutní nule se hmota nachází ve stavu s nejnižší možnou energií, nikoli však v klidu. 4) Spektrum operátoru energie je ekvidistantní, rozdíl dvou libovolných sousedních energetických hladin je   nn EEE 1 : To je právě známý Planckův vztah z počátku století. Energie jakýchkoli kmitů se nemůže měnit spojitě, ale po skocích (energetických kvantech) .E (2.61) 5) Zde se také nachází jedna z prvních možností experimentálního určení Planckovy konstanty měřením energetických kvant (například při fotoelektrickém jevu: vyrážení elektronů z povrchu kovu kvanty energie elektromagnetického záření – fotony). Zatím byla Planckova konstanta jediným neurčeným parametrem základních postulátů kvantové teorie. 6) Polynomiální řešení pro funkci u se nazývají Hermitovy polynomy a označujeme je Hn(). Pro dané n nejprve určíme bezrozměrné vlastní číslo n 2 2( 1 2) 2 1n n E n n             a z rekurentní formule (2.59) určíme pomocí c0 nebo c1 (podle toho zda jde o sudý či lichý polynom) ostatní koeficienty rozvoje. Pro c0 ≠ 0, c1 = 0 nebo c0 = 0, c1 ≠ 0 se nalezené polynomy nazývají Hermitovy. Prvních několik Hermitových polynomů vychází: )15434()(,)21()( ,)3441()(,)( ,)32()(,)( 53 15 2 02 42 0411 3 1300       cHcH cHcH cHcH Koeficienty c0 a c1 jsme volili rovny jedné. Stupeň polynomu n udává současně počet nulových bodů polynomu (počet průsečíků s osou ). 7) Hermitovy polynomy se prakticky snadno počítají nenormované z rekurentní formule )(2)(2)( 11    nnn HnHH . Pro první polynomy vychází:    12016032)(,24)( ,124816)(,2)( ,128)(,1)( 35 5 2 2 24 41 3 30    HH HH HH Normovací koeficienty (i s exponenciální vahou exp[−2 ]) jsou dány vztahem n n n 2! 1 2/1    . 8) Celkové řešení spektrálního problému je 2 /21 , | ( ) ( )e ; 0,1, 2, 2 n n n nE n n H n                  (2.62) Vlastní funkce )(n tvoří přirozený úplný ortonormální systém na prostoru L2 (–, +), který pro    „dosti rychle“ ubývá k nule. 9) Hustota pravděpodobnosti, že částice kmitající s energií En (oscilátor ve stavu |n>) se nachází v poloze x (resp. bezrozměrné poloze ), je dána výrazem nnnw  * . Pro několik prvních stavů je vykreslena na obrázku. Pravděpodobnost má oscilující charakter a existuje malá nenulová pravděpodobnost výskytu oscilátoru i za klasickými body obratu. Tento obraz nastává pro systémy s nízkou teplotou a je zcela odlišný od klasického řešení. Pro velké energie (vysoká n) by se měla křivka blížit klasické pravděpodobnosti výskytu oscilátoru (1.49). Vidíme však, že Kvantová teorie Harmonický oscilátor 40 oscilace jsou sice velmi husté, ale existuje značné množství bodů, ve kterých je kvantová pravděpodobnost nulová. Nic takového však u makroskopických systémů neměříme. Proč? To je dáno rozlišovací schopností makroskopických přístrojů. Žádný přístroj nebude měřit polohu s takovou přesností, aby registroval jednotlivá minima pravděpodobnosti u vysokých energetických stavů. Přístroj ve skutečnosti určuje polohu s konečnou přesností, do které se vejde řada minim a registruje jen střední hodnotu hustoty pravděpodobnosti. A tou je právě klasická křivka. – A – A – AA A A w w wwklas wklas wklas w0 w1 w11 2.4.2 Řešení bez volby reprezentace (Dirac) Úlohu (2.49) budeme nyní řešit obecně. Hamiltonův operátor nejprve přepíšeme do bezrozměrného tvaru:   .ˆ 2 1ˆ 2 ˆ 2 ˆ ˆ 2 1ˆˆˆ 22 2 22 PX HP XPXH      m m m m , (2.63) Převedení do bezrozměrného tvaru naprosto není nutné, veškeré další úvahy by bylo možné provádět i s rozměrovým hamiltoniánem a všechny následující vztahy by se lišily o konstantu  , kterou jsme hamiltonián vydělili. Důvodem je to, že vztahy získané z bezrozměrného hamiltoniánu jsou poněkud názornější. Pro komutující čísla je možné součet kvadrátů „odmocnit“ pomocí vztahu )i()i(22 bababa  . U nekomutujících objektů není situace tak jednoduchá. Zaveďme operátory: .; PXaPXa ˆ 2 1 iˆ 2 ˆˆ 2 1 iˆ 2 ˆ      m m m m   (2.64) Oba tyto operátory jsou pro kvantovou teorii velmi důležité. Nazývají se anihilační a kreační operátory (smysl tohoto názvu uvidíme za chvíli). Kreační a anihilační operátory, jako jedny z mála v kvantové teorii, nejsou hermitovské a nepůsobí tedy v obou částech skalárního součinu stejně. Platí pro ně některé důležité relace, například:     ˆ 1ˆ ˆ(1) , 2 ˆ 1ˆ ˆ(2) , 2 ˆ ˆ ˆ(3) , 2 ˆ ˆ ˆ(4) i , 2 ˆ ˆ ˆ(5) , , ˆ ˆ ˆ(6) , , ˆˆ ˆ(7) , . m m                                             H a a H a a X a a P a a H a a H a a a a 1 (2.65) Kvantová teorie Harmonický oscilátor 41 Důkaz všech relací je triviální. Stačí jen dosadit z definice kreačních a anihilačních operátorů  aˆ , aˆ (2.64) a využít základní komutační relace 1PX ˆiˆ,ˆ ][ . Relace (1) a (2) jsou zobecněním relace )i()i(22 bababa  pro nekomutující objekty a představují formální odmocnění hamiltoniánu. Kreační a anihilační operátory jsou lineární kombinací operátoru souřadnice a operátoru hybnosti. Proto je možné naopak operátor souřadnice a hybnosti vyjádřit jako lineární kombinaci kreačních a anihilačních operátorů – viz relace (3) a (4). Známe-li kreační a anihilační operátor, můžeme z relací (1) až (4) zpětně zrekonstruovat hamiltonián, operátor souřadnice a operátor hybnosti. Komutační relace (5) až (7) vyjadřují základní vlastnosti kreačních a anihilačních operátorů: Uvidíme, že relace (5) znamená, že anihilační operátor posouvá stavy systému o energetickou hladinu  dolů a relace (6) znamená, že kreační operátor posouvá stav o energetickou hladinu  vzhůru. Relace (7) je potom vzájemnou relací mezi anihilačním a kreačním operátorem. V následující větě dokážeme, že operátor  aˆ je kreačním operátorem, tj. posouvá energetické stavy o jednotku vzhůru (kreuje, vytváří energetické kvantum). Věta: Nechť  nEn n ||ˆH . Potom  1|~|ˆ nna . Důkaz: .1|~|ˆ |ˆ)(|ˆˆ |ˆ)(|)ˆˆ(|)ˆˆˆ(|ˆˆ )6(         nn nEn nEnEnn n nn a aaH aaaaHaaH     Zcela analogicky můžeme z relace (5) ukázat, že pro anihilační operátor platí  1|~|ˆ nna . Zavedeme-li normovací konstanty (požadujeme, aby všechny stavy byly normovány k jedné, tj. ortonormální bázi z vlastních vektorů operátoru energie), můžeme posouvání v energetickém spektru prováděné kreačním a anihilačním operátorem jednoduše zapsat jako .1||ˆ ,1||ˆ     nn nn n n   a a (2.66) Normovací konstanty  určíme později. Nyní naše úsilí zaměříme na nalezení spektra Hamiltonova operátoru pro harmonický oscilátor. Hamiltonův operátor je součtem kvadrátů dvou Hermitových operátorů a je proto pozitivně definitní, tj. jeho vlastní čísla jsou nezáporná. Kreační a anihilační operátory posouvají ve spektru energie o konstantní hodnotu (energetické kvantum). Musí tedy existovat stav s nejnižší možnou energií, která je nezáporná. Tento stav nazýváme základní stav a označujeme ho |ZS>. Zapůsobíme-li na základní stav anihilačním operátorem, musíme dostat nulový vektor |0>, jehož velikost je nulová a netvoří paprsek v Hilbertově prostoru, tedy nejde o reálný fyzikální stav (není již co anihilovat, jsme v základním stavu s nejnižší možnou energií). Pro základní stav tedy platí: 0 ˆ ˆ| ZS | ZS ; | ZS | .E     0H a Nalezněme kvadrát velikosti poslední relace (skalární součin prvku se sebou samým): Kvantová teorie Harmonický oscilátor 42 (2.65.1) 0 0 0 ˆ ˆ ˆ ˆZS| | ZS 0 ZS| | ZS 0 2 1 1 1ˆ ˆZS| | ZS ZS| | ZS 0 ZS| ZS 0 2 2 1 0 . 2 2 E E E                                          H 1 a a H 1 Známe-li hodnotu základního energetického stavu, můžeme další hodnoty energií získat působením kreačního operátoru, ten posouvá v energii o konstantu  , je tedy jasné, že    , 2 5 2 , 2 3 02 01     EE EE 0 1 ; 0,1,2, . 2 nE E n n n               Spektrum harmonického oscilátoru jsme získali jen z vlastností Hamiltonova operátoru, resp. jen z formulace úlohy (2.49). Nikde jsme nevolili konkrétní reprezentaci, konkrétní Hilbertův prostor. Kreační a anihilační operátory, se kterými jsme se zde poprvé setkali, mají značný význam v kvantové teorii pole, kde pomocí podobných operátorů kreujeme a anihilujeme jednotlivé částice přítomné v systému. Zde u harmonického oscilátoru jen kreujeme či anihilujeme energetické kvantum a tím se dostáváme o jednu hladinu výše nebo níže. Aby naše odvození bylo úplné, určíme na závěr normovací konstanty ve výrazu (2.66). Vytvořme nejprve kvadrát normy obou relací: 2 (2.65.1,2) 2 ˆ ˆ| | 1| 1ˆ | | 1 ˆ | | 1 ˆ ˆ| | 1| 1 nn n n n n n nn n n n n n n n                                 aaa a a a 22 22 ˆ 11 | 1| 1| | 1| 1 22 . ˆ 11 | 1| 1| | 1| 1 22 n nn n nn E n n n nn n n n E n n n nn n n n                                                    H H Z požadavku normovanosti vlastních vektorů operátoru energie k jedné máme: 2 2 1 ( 1 2) 1 1 2 2 1 ( 1 2) 1 . 2 2 n n n n E n n E n n                                              Fázový faktor při odmocňování komplexního čísla není podstatný (jednotkovou velikost vektoru neovlivní). Výsledné působení kreačního a anihilačního operátoru (2.66) i s normovací konstantou tedy je: .1||ˆ ,1|1|ˆ   nnn nnn a a (2.67) Kvantová teorie Harmonický oscilátor 43 Tento výsledek si můžete snadno zapamatovat: Pod odmocninou je vždy pořadové číslo vyššího energetického stavu z obou stran rovnice. Zajímavé vlastnosti má ještě jeden operátor: .ˆˆˆ aaN   (2.68) Zapůsobme s tímto operátorem na n-tý energetický stav: (2.67) ˆ ˆ ˆ ˆ| | | 1 | | .n n n n n n n n n           N a a a Vlastním číslem tohoto operátoru je počet kvant přítomných v daném energetickém stavu. V kvantové teorii pole má tento operátor význam operátoru počtu částic. 2.4.3 Řešení pomocí maticové mechaniky (Heisenberg) Řešme nyní ještě jednou úlohu (2.49) o harmonickém oscilátoru na prostoru nekonečných posloupností l 2 sčitatelných s kvadrátem. Na prostoru n-tic jsou operátory matice nn. Na prostoru nekonečných posloupností (n  ) budou operátory nekonečně rozměrné matice. Úkol tedy je: najít nekonečně rozměrné matice X, P, H, které vyhovují úloze (2.49). Tyto matice nemusíme hledat „na zelené louce“. S tím co víme o kreačních a anihilačních operátorech je snadno zkonstruujeme. Tyto matice najdeme v energetické reprezentaci, to znamená, že určíme maticové elementy operátorů HPX ˆ,ˆ,ˆ v bázi vytvořené z vlastních vektorů Hamiltonova operátoru. Všechny tři operátory umíme zkonstruovat pomocí kreačních a anihilačních operátorů podle relace (2.65). A působení kreačních a anihilačních operátorů na zvolenou bázi také známe – viz relace (2.67). Konstrukce elementů příslušných matic je tedy víceméně triviální záležitostí:       (2.65) (2.67) , 1 , 1 ˆ ˆ ˆ| | | | 2 1 | 1 | 1 2 1 , , 0,1,2. 2 k l k l k l X k l k l m l k l l k l m l l k l m                                   X a a       (2.65) (2.67) , 1 , 1 ˆ ˆ ˆ| | i | | 2 i 1 | 1 | 1 2 i 1 , , 0,1,2. 2 k l k l k l m ω P k l k l m ω l k l l k l m ω l l k l                               P a a (2.65) (2.67)1ˆ ˆ ˆ| | | | 2 1 , , 0,1,2. 2 k l kl H k l ω k l l ω k l                           H a a Kvantová teorie Harmonický oscilátor 44 Napišme si tyto matice: , 0300 3020 0201 0010 2                       m X , 0300 3020 0201 0010 2 i                       m P . 2 5 00 0 2 3 0 00 2                              H Ověřte si, že skutečně  , i X P 1 a H = P2 /2m + m2 X2 /2 podle požadavků úlohy (2.49). Ze znalosti matic X a P jsme již mohli Hamiltonovu matici určit přímo z této relace. Poslední co zbývá je nalézt vlastní čísla matice H. Tato úloha je mimořádně jednoduchá. U diagonální matice jsou vlastní čísla právě prvky na diagonále. Pokud tento fakt nevíte, výpočet je jednoduchý: | | ( ) | 0 det ( ) 0 3 5 0 2 2 2 1 , 0, 1, 2, . 2 n Ε E E Ε Ε Ε E n n                                                         H H 1 H 1 Opět tedy máme vztah (2.60) pro spektrum harmonického oscilátoru. Kvantová teorie Sférický potenciál 45 2.5 SFÉRICKY SYMETRICKÝ POTENCIÁL Sféricky symetrickým (centrálním) nazýváme potenciál, který závisí jen na vzdálenosti od určitého centrálního bodu. Pro popis pohybu těles v sféricky symetrickém potenciálu je velmi výhodný sférický souřadnicový systém. Mezi nejznámější sférické potenciály patří sférický harmonický oscilátor, sférická jáma a Coulombův potenciál. Sférický oscilátor si můžete představit jako tělísko v počátku souřadnic, od kterého vedou pružiny na všechny strany. Kdykoli ho vychýlíme, bude působit vratná síla směrem do centra. Průběh potenciální energie je kvadratický. Sférická jáma přibližně odpovídá potenciálu, který pociťuje neutron zachycený v atomovém jádře. Jaderné síly na hranici jámy (r = a) jsou značné (v idealizaci (2.69) dokonce nekonečné) a v jiných oblastech velmi slabé. Coulombův potenciál se uplatní například ve vodíkovém atomu, kdy osamocený elektron podléhá působení jediného protonu v atomovém jádře. Nezapomínejte, že r  (0, ). Průběhy těchto známých potenciálů jsou: . 1 4 )()3( , 0 )()2( , 2 1 )()1( 0 0 2 rr qQ rV arV ar rV krrV            (2.69) V klasické mechanice bude popsán systém Lagrangeovou funkcí, zobecněnými hybnostmi a energií a Hamiltonovou funkcí ve sférickém souřadnicovém systému takto: 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 222 2 2 2 2 2 2 1 1 1 sin ( ) , 2 2 2 , sin , , 1 1 1 sin ( ) , 2 2 2 ( ) ( ) . 2 22 2 sin 2 r r r L mr mr mr V r p mr p mr p mr E mr mr mr V r ppp p H V r V r m mmr mr mr                                L      (2.70) Již v klasické mechanice jsme si ukázali, že zobecněné hybnosti odpovídající úhlovým proměnným jsou komponenty momentu hybnosti. Druhá část Hamiltonovy funkce odpovídá rotačním stupňům volnosti a lze ji zapsat pomocí vektoru momentu hybnosti L vzhledem k ose z, od které je odvozen sférický souřadnicový systém. Z předchozího již víme, že jednotlivé komponenty momentu hybnosti nejsou současně měřitelné a nekomutují spolu (2.37). Současně ale můžeme měřit kvadrát momentu hybnosti (2.38) a libovolnou z jeho komponent (2.39). U sféricky symetrického problému budeme preferovat třetí osu a třetí komponentu. Osa z má preferované postavení při budování sférického souřadnicového systému, ve skutečnosti je však lhostejné, kterou z komponent V r (1) (2) (3) Kvantová teorie Sférický potenciál 46 momentu hybnosti zvolíme do úplné množiny pozorovatelných. Je-li v systému přítomno vnější magnetické pole, volíme zpravidla souřadnicový systém tak, aby třetí osa mířila ve směru tohoto pole, osa z je potom současně směrem vnějšího magnetického pole. Je-li systém sféricky symetrický, potom s operátory 2ˆL a 3 ˆL ještě komutuje Hamiltonův operátor Hˆ . To je vidět již z klasického rozpisu (2.70). Víme totiž, že zobecněné souřadnice nekomutují jedině se svými zobecněnými hybnostmi. V komutátoru ],[ HL ˆˆ2 mohou tedy být jediné nenulové členy s úhlovou částí hamiltoniánu, tou je ale právě násobek 2ˆL . Operátor sám se sebou komutuje, takže výsledek může být jedině nulový. Podobně komutátor ],[ HL ˆˆ 3 může mít jediné nenulové části s úhlovou částí hamiltoniánu, tj. ~ ],[ 2 3 ˆˆ LL . Tento komutátor je ale opět nulový podle(2.39). Nalezli jsme tak trojici nezávislých komutujících operátorů, která tvoří úplnou množinu pozorovatelných u nerelativistického sféricky symetrického problému (v relativistické úloze k těmto proměnným ještě přibude spin): .0ˆˆˆˆˆˆ 3 2 3 2  ],[],[],[ HLHLLL (2.71) U soustavy nezávislých vzájemně komutujících operátorů je možné hledat společné vlastní vektory ke všem operátorům. U sféricky symetrického problému budeme tedy řešit soustavu tří rovnic pro vlastní vektory .,,|,,|ˆ ,,,|,,|ˆ ,,,|,,|ˆ 3 2    mlνmlν mlνmlν mlνEmlν m l ν   L L H (2.72) Index  čísluje energetické stavy, index l stavy kvadrátu momentu hybnosti a index m stavy projekce momentu hybnosti do libovolné osy (zvolili jsme třetí). Vlastní čísla jsme označili E, , . Tuto soustavu je třeba řešit současně. Co by se stalo, kdybychom například řešili jen rovnici pro energii? Nalezená vlastní čísla E by samozřejmě byla v pořádku, ale ke každému vlastnímu číslu (každé hodnotě energie) by existovalo více nezávislých vlastních vektorů (ve skutečnosti se od sebe liší čísly l a m, to ale nevíme, protože řešíme jen první rovnici). Tomuto typu spektra říkáme degenerované spektrum. Znamená to jen to, že k danému vlastnímu číslu existuje více vlastních vektorů. Odlišili bychom je od sebe až pomocí dalších operátorů, které komutují s operátorem, jehož spektrum právě hledáme. V následujících dvou kapitolách se budeme zabývat momentem hybnosti, tedy druhou a třetí rovnicí v (2.72). Řešení pro moment hybnosti je stejné pro všechny průběhy potenciální energie. V kapitole 2.5.1 nalezneme řešení bez použití konkrétní reprezentace a v kapitole 2.5.2 naznačíme, jak by se při řešení postupovalo v x reprezentaci. První rovnicí v (2.72) se budeme zabývat v kapitole 2.5.3. Řešení pro energii (energetické spektrum) již samozřejmě závisí na průběhu potenciální energie a je jiné například pro vodík a jiné pro sférický oscilátor. Navíc řešení pro energii závisí na číslech l a m. To je logické: moment hybnosti souvisí s rotačními stavy systému a ty k energii přispívají. Vidíme to konec konců i v hamiltoniánu (2.70), kde je právě rotační část energie vyjádřena přes kvadrát momentu hybnosti. Kvantová teorie Sférický potenciál 47 2.5.1 Moment hybnosti Základními komutačními relacemi pro moment hybnosti jsou vztahy (2.37) a (2.39): .0]ˆ,ˆ[ ˆi]ˆ,ˆ[ 3 2 321   LL LLL záměny,cyklické Zaveďme nyní tzv. posuvné operátory .ˆiˆˆ 21 LLL  (2.73) Tyto operátory budou mít podobný význam jako kreační a anihilační operátory u energie harmonického oscilátoru. Budou nás totiž posouvat ve spektru momentu hybnosti. Napišme přehledně jejich důležité vlastnosti (všechny lze snadno odvodit z definice posuvných operátorů a z komutačních relací momentu hybnosti):     1 2 2 2 3 3 2 2 3 3 3 3 2 1ˆ ˆ ˆ(1) , 2 1ˆ ˆ ˆ(2) , 2i ˆ ˆ(3) , ˆ ˆ ˆ ˆ ˆ(4) , ˆ ˆ ˆ ˆ ˆ(5) , ˆ ˆ ˆ(6) , 2 , ˆ ˆ ˆ(7) , , ˆ ˆ(8) , 0 .                                                L L L L L L L L L L L L L L L L L L L L L L L L L L (2.74) Známe-li posuvné operátory, můžeme z relací (1), (2) a (6) zrekonstruovat celý moment hybnosti. Úlohu, kterou budeme nyní řešit, lze zformulovat takto: .,|,|ˆ ,,|,|ˆ 3 2     L L Dokažme nejprve, že posuvné operátory posouvají vlastní vektory ve třetí komponentě momentu hybnosti o Planckovu konstantu: Lemma 1: .,|~,|ˆ  L Důkaz: Označme    ,|ˆ| L . Aplikujme operátory 3 ˆL a 2ˆL na tento vektor: (2.74.7) 3 3 3 ˆ ˆ ˆ ˆ ˆ ˆ ˆ| | , ( ) | , ( ) | , ( )|                         L L L L L L L , (2.74.8) 2 2 2ˆ ˆ ˆ ˆ ˆ ˆ| | , | , | , |                   L L L L L L . Vidíme, že posuvné operátory Lˆ posouvají ve spektru operátoru 3 ˆL o konstantu  . Ve spektru operátoru 2ˆL nedělají posuvné operátory nic. Lˆ tedy mění jen hodnotu projekce momentu hybnosti do zvolené osy. Kvantová teorie Sférický potenciál 48 Lemma 2: Při daném λ je spektrum operátoru 3 ˆL omezené, tj. existuje μmin a μmax . Důkaz: V relacích     2 2 2 2 2 2 1 2 3 ˆ | , | , , ˆ ˆ ˆ ˆ| , | , ( ) | ,                        L L L L L jsou operátory na levých stranách pozitivně definitní. Proto musí platit 0,0 2   . Zřejmě tedy musí být 02   , a proto   , a existuje min a max. Nyní již spektrum momentu hybnosti odvodíme standardním způsobem. Podobně jako u harmonického oscilátoru zapůsobíme posuvným operátorem na první (resp. poslední stav). Výsledek působení musí být nulový, protože další stav již neexistuje: .0,|ˆ0,|ˆ minmax    LL vytvořme kvadrát normy těchto vektorů: .0,|ˆˆ|,0,|ˆˆ|, minminmaxmax    LLLL Součiny operátorů vyjádříme z (2.64.4) a (2.64.5). .0,|ˆˆˆ|,0,|ˆˆˆ|, min3 2 3 2 minmax3 2 3 2 max   LLLLLL  Po zapůsobení operátorů máme: .0||,||)(0||,||)( 2 minmin 2 min 2 maxmax 2 max    vynulováním koeficientů u obou relací dostáváme: ,min 2 minmax 2 max    neboli .)()( minminmaxmax    (*) Posuvné operátory posouvají ve spektru třetí komponenty momentu hybnosti o Planckovu konstantu, proto musí také současně platit: maxmaxminminminmin ,,,3,2,,    . zaveďme bezrozměrné číslo /m . Potom maxmaxminminminmin ,1,,3,2,1, mmmmmmm   . Lemma 3: Označíme-li mmax  l, potom je mmin = – l. Důkaz: Z relací (*) snadno zjistíme, že 2 min min min min min min ( 1) ( 1) ( 1) 0 11 1 4 ( 1) 1 (1 2 ) . 2 2 l l m m m m l l ll l l m m l                     První řešení je ve sporu s předpokladem, druhé dokazuje uvedené tvrzení. Číslo m tedy může nabývat celkem 2l+1 různých hodnot z množiny },1,,2,1,{ lllllm   . Počet hodnot 2l+1 musí být nezáporné celé číslo a proto samo číslo l může nabývat jen poločíselných hodnot Kvantová teorie Sférický potenciál 49 }, 2 5 ,2, 2 3 ,1, 2 1 ,0{ l . Vlastní číslo )1()()( 2 maxmax  llll  . Závěr: Výsledky celého odvození můžeme zformulovat takto: .1,|~,|ˆ ;},1,,2,1,{,,|,|ˆ ;}, 2 5 ,2 2 3 ,1, 2 1 ,0{,,|)1(,|ˆ 3 22     mlml lllllmmlmml lmlllml L L L   (2.75) Poznámky k řešení: (velmi důležité, čtěte pozorněji než samo řešení!!!) 1) Číslo l čísluje velikost momentu hybnosti a nazývá se vedlejší kvantové číslo (hlavní kvantové číslo čísluje energii). Číslo m čísluje projekci momentu hybnosti do libovolné osy. Vzhledem k tomu, že nabitá rotující částice má nenulový magnetický moment, a toto číslo bylo poprvé zavedeno pro elektron v atomárním obalu vodíku, nazývá se magnetické kvantové číslo. 2) Možné hodnoty velikosti momentu hybnosti a jeho projekce do třetí osy jsou: .,,1,, ;...,3,2,1,0,)1(|| 3 lllmmL lllL     (2.76) 3) Poločíselné hodnoty, které jsme odvodili pro číslo l jsou skutečně také možné. Realizují se u spinu, jehož operátor má stejnou komutační strukturu jako moment hybnosti. V Schrödingerovské X reprezentaci (následující kapitola) tyto hodnoty nedostaneme. Volba reprezentace zde znamená ztrátu části řešení. To, že poločíselné hodnoty l jsou již součástí komutačních relací (2.37) bylo objeveno až relativně pozdě (v roce 1968 Kaufmannem) postupem podobným našemu odvození. 4) Z výsledku (2.76), respektive (2.75) plyne skutečný význam Planckovy konstanty. Jedná se o elementární kvantum momentu hybnosti. Při měření momentu hybnosti budeme vždy měřit projekci momentu do určité osy, dané měřícím zařízením. Tato projekce je vždy násobkem Planckovy konstanty. 5) Vidíme, že stavy s konkrétním vedlejším kvantovým číslem l jsou degenerovány – existuje více vlastních vektorů | l, m>, které přísluší stejnému kvantovému číslu l. Tyto vektory se od sebe liší kvantovým číslem m a jejich počet je 2l+1 (tzv. stupeň degenerace, který označujeme #). 6) Historicky byly označovány kvantové stavy velikosti momentu hybnosti elektronu v obalu atomu vodíku písmeny s, p, d, f podle následující tabulky: 7#stav 5#stav 3#stav 1#stav     3,2,1,0,1,2,33 2,1,0,1,22 1,0,11 00 mfl mdl mpl msl 7) Vztah pro velikost kvadrátu momentu hybnosti lze dostat také jako aritmetický průměr všech možných hodnot. Například pro l = 2 jsou možné hodnoty projekcí Lx, Ly nebo Lz rovny –2, –, 0, , 2. Průměrná hodnota kvadrátu je proto dána vztahem .65)404(33 2222222222   zzyx LLLLL Velikost 6|| L přesně podle vztahu (2.76). 8) Není obtížné napočítat maticové elementy  mlml k ,|ˆ|, L operátoru momentu hybnosti ve vlastní reprezentaci pomocí posuvných operátorů podobně jako u harmonického oscilátoru v kapitole (2.4.3). Pro l = 0 může být m i m jen 0 a proto jde o jediný prvek. Tato matice působí 2 -2 0 Lz l = 2   - Kvantová teorie Sférický potenciál 50 na skalární veličiny, hovoříme o skalární reprezentaci. Pro l = 1/2 může nabývat m i m hodnot –1/2 a +1/2. Jde o matice 2×2 působící na uspořádané dvojice, které nazýváme spinory. Jedná se o tzv. spinorovou reprezentaci. Pro l = 1 může nabývat m i m hodnot –1, 0 a +1. Jde o matice 3×3 působící na uspořádané trojice, které nazýváme vektory. Jedná se o tzv. vektorovou reprezentaci. Všimněte si, že matice L3 jsou diagonální s vlastními čísly na diagonále. Spinorová reprezentace ( l = 1/2) 1 2 3 0 1 0 i 1 0 ; ; . 1 0 i 0 0 12 2 2                            L L L (2.77) Vektorová reprezentace ( l = 1) 1 2 3 0 1 0 0 i 0 1 0 0 1 0 1 ; i 0 i ; 0 0 0 . 0 1 0 0 i 0 0 0 1                                    L L L (2.78) Matice pro l = 1/2 se nazývají Pauliho matice (bez násobících koeficientů). 9) Známé tvrzení Bohrova modelu, že na obvod dráhy elektronu v atomárním obalu musí připadnout celistvý násobek vlnových délek, je možné s pomocí vztahu (2.3) přepsat takto: 2 2 2n n n n n n r n r m r n m          v v a nejde tedy o nic jiného než o kvantování projekce momentu hybnosti. 2.5.2 Řešení v x reprezentaci, kulové funkce V x reprezentaci budeme problém sférického potenciálu řešit ve sférických souřadnicích (jsou nejbližší symetrii potenciální energie). Je třeba řešit soustavu rovnic (2.72), která bude mít nyní tvar: 2 3 ˆ ( , , ) ( , , ) , ˆ ( , , ) ( , , ) , ˆ ( , , ) ( , , ) . ν l m r E r r r r r                        H L L (2.79) Operátory zapsané ve sférických souřadnicích mají tvar: 2 2 2 2 2 2 2 3 2 2 2 2 2 2 2 2 2 ˆˆˆ ( ) ( ) , 2 2 2 2 ˆ , ˆ i ; ˆ1 , 1 , 1 1 sin . sin sin r r r r r V r V r m I m mr r r r r rr                                                               p L H L L L (2.80) Kinetická energie v Hamiltonově operátoru vede v Schrödingerově rovnici na člen Kvantová teorie Sférický potenciál 51   . 1 222 2 222         rmmm rzyx  V kartézských souřadnicích se Laplaceův operátor štěpí na součet druhých derivací podle jednotlivých os, tomu odpovídá rozklad kinetické energie na složky Tx, Ty a Tz. Ve sférických souřadnicích se Laplaceův operátor dělí na radiální a úhlovou část, tomu odpovídá rozklad kinetické energie na radiální a úhlovou část. Právě úhlová část kinetické energie je rotační energie spojená s momentem hybnosti a proto kvadrátu momentu hybnosti odpovídá úhlová část Laplaceova operátoru. Hledané řešení ),,(  r samozřejmě závisí na kvantových číslech , l, m. Řešení budeme hledat v separovaném tvaru )()()(),,(  hgrfr  . (2.81) Nejdříve řešme poslední z rovnic (2.79):    i ( ) ( ) ( ) ( ) i ( ) exp i . m m m f r g h f r g h dg g g c d                              Nalezené řešení musí být periodické v úhlu :  ,2,1,0;)2()0(  mmgg m V x reprezentaci jsme opět odvodili kvantování projekce momentu hybnosti. Projekce momentu hybnosti může nabývat jen celistvých násobků Planckovy konstanty. Poločíselná řešení nelze v x reprezentaci nalézt. Přechodem ke konkrétní reprezentaci přicházíme o část řešení. Hledané řešení má nyní tvar: i1 ( , , ) ( ) e ( ) ; 0, 1, 2, 2 m r f r h m          (2.82) Konstantu c jsme zvolili tak, aby nalezené řešení bylo normováno k jedné. Jako další krok dosadíme toto řešení do druhé rovnice (2.80) a budeme ji řešit 2 2 i i 2 2 2 2 2 1 1 sin e ( ) e ( ) sin sin 1 sin 0 . sin sin m m l l h h d dh m h d d                                                  Výsledkem je obyčejná diferenciální rovnice pro funkci h(), která se řeší standardními matematickými postupy přesahujícími rámec tohoto sylabu. Výsledkem jsou polynomiální funkce v cos  a sin , které se nazývají přidružené Legendreovy polynomy Plm(cos ) a jsou definované vztahem  ,1,0;||;,2,1,0;)1( !2 )1( )( 2 22      mlmlx dx d l x xP l ml ml l m lm (2.83) Pro m = 0 se tyto polynomy nazývají Legendreovy polynomy. Příslušné vlastní číslo je 2 )1(  lll (2.84) Celá úhlová část řešení se nazývá kulová funkce a označuje se Kvantová teorie Sférický potenciál 52 i1 ( , ) e (cos ) 2 .m lm lmY P      (2.85) Celkové řešení druhých dvou rovnic soustavy (2.79) tedy je lmmm lll PrfYrfr m l lm m lm    ||;,1,0;; ,2,1,0;)1( ;)(cose)( 2 1 ),()(),,( 2 i         (2.86). Odvozené kvantování momentu hybnosti je až na absenci poločíselných hodnot shodné se vztahy odvozenými jinou cestou v předchozí kapitole. Pro radiální funkci f(r) lze řešení získat z první rovnice (2.79). Toto řešení závisí na tvaru potenciální energie. Pro některé základní tvary potenciální energie bude řešení diskutováno v příští kapitole. Na závěr uveďme příklady některých kulových funkcí: i i 00 10 11 1, 1 2 i 20 21 1 3 3 3 ; cos ; e sin ; e sin ; 4 8 82 5 15 (1 3cos ); e cos sin ; 16 8 Y Y Y Y Y Y                            2.5.3 Jednoduché systémy: oscilátor, vodík, jáma Nyní zbývá řešit první z rovnic (2.79) – rovnici pro energii. Tato rovnice nám poskytne energetické spektrum a radiální část celého řešení (r, , ). Jak energetické spektrum, tak radiální část mohou záviset na kvantových číslech l a m z předchozího řešení a budou závislé na konkrétním tvaru potenciální energie V(r). V poslední rovnici (2.79) známe působení rotační části kinetické energie Hamiltonova operátoru na celkovou vlnovou funkci. To je dáno působením kvadrátu momentu hybnosti podle druhé z rovnic (2.79). Známe již i vlastní číslo l podle vztahu (2.84). Po zapůsobení rotační části zkrátíme úhlové části g() a h() na obou stranách rovnice a získáme rovnici pro radiální část řešení: 2 2 2 2 2 1 ( 1) ( ) ( ) ( ) 2 2 l l d d l l r V r f r E f r m dr drr mr               . (2.87) Povšimněte si, že v rovnici vystupuje vedlejší kvantové číslo l a energetické spektrum proto nezávisí jen na radiálním čísle , které čísluje energii, ale i na vedlejším kvantovém čísle l. Řešení rovnice (2.87) se provádí standardními metodami (rozvoj do řady, hledání asymptotického chování, oříznutí nekonečné řady). Uvedeme výsledky výpočtů pro potenciální energii sférického harmonického oscilátoru, prostorové jámy a Coulombův potenciál (2.69). Harmonický oscilátor Pro potenciální energii harmonického oscilátoru vychází energetické spektrum    )23()232( 2 1 )( 22  nlErmrV l . (2.88) Nejmenší možná hodnota energie (nulové kmity) je 23 . Radiální kvantové číslo  čísluje pořadí radiálních stavů a zpravidla také počet průsečíků radiálního řešení s osou x. Většinou se zavádí tzv. hlavní kvantové číslo n, které skutečně čísluje stavy energie: Kvantová teorie Sférický potenciál 53 nlnln  1,0,,2,1,0;2   . (2.89) Spektrum oscilátoru je degenerované (ke každé hodnotě energie přísluší více stavů, každé n lze složit z více kombinací  a l). Snadno určíme stupeň degenerace, uvědomíme-li si, že ke každému vedlejšímu kvantovému číslu existuje 2l + 1 hodnot magnetických čísel m: . 2 )2)(1( 1421)2(212# 2/ 0 2/ 0     nn nnl nn l n   (2.90) Řadu (2.90) jsme sečetli jako aritmetickou řadu. Každá energetická slupka n obsahuje (n + 1)(n + 2)/2 stavů. Coulombický potenciál Pro Coulombickou potenciální energii vychází energetické spektrum 2 2 e e 2 2 2 2 0 1 ( ) 4 2 ( 1) 2 l m mqQ V r E r r l n                 . (2.91) Hlavní kvantové číslo n číslující stavy energie jsme zavedli vztahem 1 ; 1, 2, , 0,1 1 .n l n l n       (2.92) Stupeň degenerace bude .1221)1(212# 2 1 0 1 0 nnnl nn l n         (2.93) Jde-li o atom vodíku, může mít každý elektron ještě dva spinové stupně volnosti ms = 1/2 a celkový počet stavů v jedné energetické slupce je proto 2n2 . Tyto stavy se liší hodnotou kvantových čísel l, m, ms. Kvantová jáma Sférická konečná kvantová jáma s potenciálem        arV ar rV 0 0 )( (2.94) nemá analytické řešení. Problém lze řešit jen numericky nebo graficky. Kvantová teorie Časový vývoj 54 2.6 ČASOVÝ VÝVOJ Prozatím jsme se v kvantové teorii zabývali stacionárními stavy, tj. stavy systému, které se v čase nevyvíjí. Skutečné kvantové stavy jsou lineárními kombinacemi stacionárních stavů (prvků báze) a koeficienty těchto kombinací se mění s časem. Přechod stavu z jednoho času do času pozdějšího provádí tzv. evoluční operátor (operátor časového vývoje). 2.6.1 Evoluční operátor Evoluční operátor převádí známý stav čase t0 na stav, do kterého se vyvine v čase t:  )(|),(ˆ)(| 00 tttt  U . (2.95) Evoluční operátor musí splňovat některé podmínky a požadavky: 1) počáteční podmínka: vývoj z počátečního času do počátečního času nemění stav 1U ˆ),(ˆ 00 tt . 2) semigrupová podmínka: vývoj ze stavu t1 do t2 dopadne stejně, je-li proveden naráz nebo přes mezičas t:   )(|),(ˆ),(ˆ)(|)(|),(ˆ)(| 11221122 2121 tttttttttt ttttt  UUU Porovnáním obou postupů získáme semigrupovou podmínku ),(ˆ),(ˆ),(ˆ 1212 tttttt UUU  . 3) unitarita: časový vývoj nemění normování stavu: 0 0 † 0 0 0 0 ( ) ( ) ( ) ( ) , ˆ ˆ , t t t t          U U †ˆ ˆ ˆU U 1. 4) inverze: inverzní evoluční operátor má obrácené pořadí argumentů. Odvodíme ze semigrupové podmínky:  1UUU ˆ),(ˆ),(ˆ),(ˆ ! 00 tttttt ),(ˆ),(ˆ 00 1 tttt UU  . 5) spojitost: samovolný vývoj stavu (bez aktu měření), který popisuje evoluční operátor, musí být spojitý: 0 0 0 ˆ ( , ) ( ) je spojité pro at t t t    U H . Nyní odvodíme základní rovnici pro evoluční operátor. Vyjdeme z definice střední hodnoty dynamické proměnné (viz tabulka v kapitole 2.3.2) a tuto střední hodnotu budeme derivovat podle času: † † † 0 0 0 0 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ . da d d d d dt dt dt dt dt          U U A U AU AU U A Jinou možností je přímo zavést operátor časové derivace dynamické proměnné vztahem † 0 0 ˆ ˆ ˆ ˆ . da dt        A U AU Porovnáním obou postupů získáme rovnici Kvantová teorie Časový vývoj 55 † † † ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆd d dt dt   U U AU U A U AU, (*) ve které za časovou derivaci operátoru dynamické proměnné dosadíme časový vývoj dynamické proměnné zapsaný v Poissonových závorkách (1.53) převedený do kvantové podoby pomocí principu korespondence (2.35):   1ˆ ˆ ˆ, , i A A H         A A H . (2.96) Získáme tak rovnici, ze které se budeme snažit získat rovnici pro evoluční operátor: † † † ˆ ˆ 1ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, i d d dt dt       U U AU U A U A H U   † † † † ˆ ˆ 1ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ i d d dt dt     U U AU U A U AHU U HAU (**) Ve všech následujících úpravách využíváme unitaritu UU † = U † U = 1. Z rovnice (**) je třeba vyloučit operátor U † a jeho derivaci podle času, kterou získáme derivováním definice unitarity podle času a násobením výsledku operátorem U † zprava: † † † † † † ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ0 0 d d d d dt dt dt dt         U U U U U U 1 U U U U † † † ˆ ˆ ˆ ˆ . d d dt dt   U U U U Výsledek dosadíme do rovnice (**) a vynásobíme ji operátorem U zleva a U † zprava:  † † † † † ˆ ˆ 1ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ i d d dt dt       U U U U AU U A U AHU U HAU  † † ˆ ˆ 1ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ i d d dt dt       U U U A A U AH HA † ˆ 1ˆ ˆ ˆ ˆ, , i d dt            U A U A H † ˆ 1ˆ ˆ i d dt    U U H UH U ˆˆ ˆ i  dt d  . (2.97) Právě odvozená rovnice se nazývá rovnice časového vývoje. Zapůsobíme-li touto operátorovou rovnicí na počáteční stav |0>, provede evoluční operátor vývoj stavu do času t a získaná rovnice pro |(t)> se nazývá časová Schrödingerova rovnice: )(ˆ )( i t dt td   H . (2.98) Kvantová teorie Časový vývoj 56 2.6.2 Časová Schrödingerova rovnice Řešení časového vývoje lze najít relativně snadno, není-li Hamiltonův operátor explicitní funkcí času, tj. závisí jen na operátorech zobecněných souřadnic a hybností. V takovém případě je výhodné volit v Hilbertově prostoru popisovaného systému bázi generovanou vlastními vektory Hamiltonova operátoru (2.42): 1H ˆ;;||ˆ   n nmn nnnmnEn  . Do těchto vektorů rozvineme hledaný stav, koeficienty rozvoje budou funkcemi času:   n n ntat |)()(| . Řešení v tomto tvaru dosadíme do časové Schrödingerovy rovnice a získáme lineární rovnici pro koeficienty an(t). .e)( ;)(i zleva|/;|)(|i ;|)(ˆ|i )( i 0ttE mm mm m n nn n n n n n n m cta Eta dt ad mnEtan dt ad ntan dt ad            H Řešení časového vývoje tedy je:    n ttE n nct n |e)(| )( i 0  . (2.99) Poněkud elegantnější řešení je nalézt přímo evoluční operátor jako superpozici projektorů generovaných Hamiltonovým operátorem pomocí věty o spektrálním rozvoji. Řešení rovnice pro evoluční operátor lze formálně zapsat jako 0 0 1 1ˆ ( ) ( ) i i 0 0 ˆ ˆ ˆˆ ˆi ( , ) e ( , ) e | | nt t E t t n d t t t t n n dt           HU U UHU . Nyní zapůsobíme nalezeným evolučním operátorem na počáteční stav | 0 >: 0 1 ( ) i 0| ( ) e | | nE t t n t n n        a získáme tak okamžitě řešení časové Schrödingerovy rovnice: 0 1 ( ) i 0 0 | ( ) e | ; | | ; | . nE t t n n n n n t c n c n c n                 (2.100) Kvantová teorie Časový vývoj 57 Příklad 16: Nalezněte vývoj pravděpodobnosti systému, jehož počáteční stav je zadán jako reálná lineární kombinace dvou reálných vlastních funkcí Hamiltonova operátoru. Například může jít o dva stavy harmonického oscilátoru nebo kvantové jámy či o dvoustavový systém. Požadavek reálnosti vlastních funkcí a koeficientů je jen z důvodu jednoduchosti výpočtu. Řešení: Počáteční stav je kombinací dvou vlastních stavů 1 a 2 Hamiltonova operátoru )()()( 22110 xcxcx   , časový vývoj je )(e)(e),( 2 i 21 i 1 21 xcxcxt tEtE     a výsledná hustota pravděpodobnosti pro reálné vlastní funkce vychází ][ )( i )( i 2121 2 22 2 11 2 1212 ee)()(),( tEEtEE ccccxtw    . Celková pravděpodobnost je součtem pravděpodobnosti, že se systém nachází ve stavu 1, ve stavu 2 a interferenčního členu, který je pro kvantové procesy typický. Výsledek lze jednoduše zapsat takto: .;)cos()()()(),( 21 Etxfxwxwxtw   Frekvence časových oscilací pravděpodobnosti odpovídá Planckovu kvantování  E . ■ 2.6.3 Oscilace neutrin Neutrina (elektronové, mionové a tauonové) jsou ve skutečnosti lineární kombinací vlastních stavů hmoty k kV   , kde , ,e   a 1,2,3k  . Index α popisuje generace neutrin a index k vlastní hmotnostní stavy. Transformační matice je unitární a poprvé ji zavedli Ziro Maki, Masami Nakagawa a Shoichi Sakata v roce 1962, aby vysvětlili oscilace neutrin předpovězené Brunem Pontecorvem. Pro pochopení principu oscilací předpokládejme jen existenci dvou generací neutrin a mixáž ve tvaru e 1 2 μ 1 2 cos sin , sin cos .                 Unitární matici jsme zapsali jako běžnou rotační matici za pomoci úhlu θ. Za letu neutrin se budou hmotnostní stavy vyvíjet a mixážní poměry měnit. Spíše než časový vývoj nás ale bude zajímat vývoj stavu podél letící částice. Vzhledem k tomu, že pro rovinnou vlnoplochu platí i ( ) i( ) e e Et t      p x k x  , budeme moci vývoj hmotnostních stavů podél letu neutrina zapsat takto: i ( ) e (0) kp x k kx   . Například stav elektronového neutrina se za letu bude měnit podle formule 1 2 i i e 1 2( ) e cos (0) e sin (0) . p x p x x        Amplituda pravděpodobnosti, že se elektronové neutrino bude za letu jevit pozorovateli jako čistě mionové neutrino (dané svou počáteční kombinací), bude Kvantová teorie Časový vývoj 58 e μ μ e(0) ( )x    A . Po provedení skalárního součinu (hmotnostní stavy tvoří ortonormální bázi) máme e μ 2 1 i i cos sin exp expp x p x                     A . Neutrina mají velmi malou hmotnost a relativistické energie, a proto lze využít rozvoj   2 322 2 2 2 ( / ) 1 / 2 k k k k m cE E p E c m c m c E c c E       Amplitudu pravděpodobnosti nyní snadno upravíme e μ 2 3 1 2 3 2 3 1 2 3 2 3i 2 1 i 2 3 2 i 2 3 2 4 cos sin e exp i exp i 2 2 cos sin e exp i 1 2 cos sin e exp i 4 Ex c E m c x c E E m c m c x c E E m c m c x x E E m c x E m c x E                                                                              A 2 3 2 3 1 2 3 i 2 3 2 4 exp i 4 2icos sin e sin 4 E m c m c x c E E m c x E m c x E                                          Pokud jsou vlastní hmotnosti různé (stačí jedna nenulová), dojde k oscilacím neutrin (poprvé byla pozorována v roce 1998 na detektoru Super-Kamiokande). Pravděpodobnost, že původní elektronové neutrino nalezneme jako mionové je periodickou funkcí vzdálenosti od zdroje e μ 2 3 2 * 2 2 2 2 2 2 1sin 2 sin ; 4 m c x m m m E                P A AA . Z různých experimentů je možné určit úhel mixáže a střední vzdálenost přeměny neutrin 2 3 4 E L m c     , ze které plyne pouze rozdíl kvadrátů hmotností neutrin. Skutečná neutrina mají tři generace a transformační matice je 3×3 a obsahuje tři úhly. Princip se ale nemění. Z měření plyne, že platí Δm12 2 = (7,59 ± 0,21)×10−5 eV2 (KamLAND, 2005) a Δm23 2 = (2,43 ± 0,13)×10−3 eV2 (MINOS, 2006). Pro mixážní úhly máme přibližně θ12 ~ 33°, θ23 ~ 45°, θ13 ≤ 9°. Mixážní matice se tak rozhodně nepodobá diagonální matici, jako je tomu v případě mixáží kvarků. 2.6.4 Dvouštěrbinový experiment Představme si, že na dvě štěrbiny dopadá proud částic. Po průchodu štěrbinami se na stínítku zaznamenává, kam která dopadla. Výsledkem je klasický interferenční vlnový obrazec s maximem dopadů paradoxně mezi oběma štěrbinami. Podobně jako v předchozí kapitole se sčítají amplitudy pravděpodobností obou možností, nikoli samotné pravděpodobnosti. Na výsledku nic nezmění ani počet přítomných částic: bude-li tok zleva velmi slabý a v průměru se bude vyskytovat v oblasti experimentu jediná částice, nikdy nezjistíme, kterým otvorem prošla. Po dosti dlouhé době získáme statistický obraz dopadu částic na stínítko podle obrázku. Můžeme si třeba myslet, že část částice prošla jedním otvorem a část druhým, nebo že interferovala sama se sebou. Takové úvahy nemají reálný smysl. Pro Kvantová teorie Časový vývoj 59 posouzení statistického výsledku mnoha opakovaných dopadů je důležitý jen souhlas experimentálního výsledku s předpovědí danou teorií. počet částic stínítko ? Jiný obraz se nám naskytne, pokusíme-li se zjistit, kudy částice prolétla. Zakryjeme-li jeden z otvorů, bude maximum dopadajících částic proti otevřenému otvoru. Můžeme vymyslet rafinovanější postup. Budeme sledovat například pomocí částic světla – fotonů, kudy částice prolétla. Bude-li foton málo energetický, bude mít příliš dlouhou vlnovou délku na to, aby určil, kudy částice prolétla. Bude-li ale foton mít pro detekci dosti krátkou vlnovou délku, můžeme skutečně rozhodnout, kudy prolétla částice. Ale něco za něco: krátkovlnný foton má značnou energii a silně ovlivní stav prolétající částice. Dokonce natolik, že interferenční obrazec zcela vymizí. Obecně platí: nepokusíme-li se o detekci, sčítají se amplitudy pravděpodobnosti a statistika dopadů má charakter interferenčního jevu. Pokusíme-li se o detekci, interference zaniká a sčítají se klasicky samotné pravděpodobnosti. Těžko se nám tento fakt přijímá. Je to vlastnost mikrosvěta, která se nám zdá velmi podivná. Naše zkušenosti z makrosvěta jsou založeny na komutujících objektech. Právě nekomutativnost jevů v mikrosvětě vede ke skládání amplitud pravděpodobností možností, které jsou k dispozici, a k interferenčnímu jevu. počet částic stínítko ? detekce 2.6.5 Ehrenfestovy teorémy, viriálový teorém V této kapitole si probereme tři základní teorémy týkající se časového vývoje. První Ehrenfestův teorém První teorém se týká časového vývoje operátoru souřadnice. Pro jednoduchost ho odvodíme v jednorozměrném případě, vyjdeme z principu korespondence a časového vývoje (1.53):              . ˆ ˆiˆi i2 1ˆˆ,ˆˆ,ˆˆ i2 1 )ˆ(,ˆ1ˆ,ˆ i2 1 )ˆ( 2 ˆ ,ˆ i 1ˆ,ˆ i 1ˆ 2 2 mmm V im V mdt d P PPPPXPXP XXPXX P XHX X               První Ehrenfestův teorém je tak analogií definice hybnosti z klasické mechaniky: Kvantová teorie Časový vývoj 60 mdt d PX ˆˆ  . (2.101) Druhý Ehrenfestův teorém Druhý Ehrenfestův teorém se týká časového vývoje operátoru hybnosti. Budeme postupovat podobně jako v předchozím případě:         .)ˆ(,ˆ1 )ˆ(,ˆ1ˆ,ˆ i2 1 )ˆ( 2 ˆ ,ˆ i 1ˆ,ˆ i 1ˆ 2 2 XPXPPPX P PHP P V i V im V mdt d            Hodnotu posledního komutátoru určíme takto: Nejprve nalezneme komutátor operátoru hybnosti s libovolnou mocninou operátoru souřadnice (indukcí) a výsledek budeme člen po členu aplikovat na operátor potenciálu rozvinutý do mocninného Taylorova rozvoje: 2 1 1 ˆ ˆ ˆ ˆ ˆ, , i , ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, , , i 2 ˆ ˆ ˆ, i , ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ, , , i ( 1) , ˆ ˆ, ( ) i . ˆ n n n n n n n n V V                                                              P X X P 1 P X X P X P X X X P X X P X X P X P X X X P X X Základním předpokladem těchto úvah je samozřejmě rozvinutelnost potenciální energie do Taylorovy řady. Po dosazení za vypočtený komutátor druhý Ehrenfestův teorém vychází: X P ˆ ˆ    V dt d , (2.102) což je vlastně kvantovou analogií Newtonových pohybových rovnic (záporně vzatý gradient potenciální energie je působící sílou). Viriálový teorém Viriálový teorém je velmi užitečný nejen v kvantové teorii, ale i ve statistické fyzice. Určuje střední hodnotu kinetické energie obsažené v systému z tvaru energie potenciální. Určeme nejprve maticové elementy komutátoru dynamické proměnné A s Hamiltonovým operátorem v energetické reprezentaci: .)(|ˆ|)(|ˆˆˆˆ||]ˆ,ˆ[| nmnmnm AEEmnEEmnmn  AAHHAHA Pro n = m máme .0|]ˆ,ˆ[|  nn HA Za operátor dynamické proměnné A budeme nyní volit součin souřadnice a hybnosti: .0| ˆ ˆ||ˆ ˆ | ,0|]ˆ,ˆ[ˆ||ˆ]ˆ,ˆ[| ,0|]ˆ,ˆˆ[|    n dt d nn dt d n nnnn nn P XP X HPXPHX HPX Za časový vývoj souřadnice a hybnosti dosadíme z Ehrenfestových teorémů: Kvantová teorie Časový vývoj 61 .| ˆ ˆ 2 1 || 2 ˆ | 2     n V nn m n X X P Ve třech dimenzích je výsledek součtem příspěvků v jednotlivých osách. Na levé straně stojí střední hodnota kinetické energie systému, napravo tzv. operátor viriálu:     n V nnn k k | ˆ ˆ 2 1 ||ˆ| X XT . (2.103) Pro jednorozměrný harmonický oscilátor je operátor viriálu přímo roven potenciální energii: .ˆ 2 1 ˆ ˆ 2 1ˆ 2 1 )ˆ( 22 X X XXX k V kV     Střední hodnoty kinetické a potenciální energie jsou si proto v každém stavu rovny. Poznámka: Již v roce 1933 upozornil F. Zwicky, že v kupě galaxií ve Vlasech Bereniky je pohyb galaxií větší, než by odpovídalo viriálovému teorému pro gravitační potenciální energii. Řešením je existence další neviditelné (temné) hmoty v této kupě. Později byl podobný problém zjištěn Verou Rubinovou i pro oběžné rychlosti hvězd v periferních oblastech samotných galaxií. Řešením je opět existence haló z temné hmoty v okolí galaxie. Viriálový teorém může být proto velmi užitečný i pro makroskopické nekvantové systémy. Svítící (registrované) hmoty v galaxiích je jen asi 1 %. V roce 2000 se pomocí HST ukázalo, že až 50 procent hmoty Galaxie může být soustředěno ve velmi starých a málo svítících bílých trpaslících, které doposud nebyly pozorovatelné. Patřily pravděpodobně k první generaci hvězd před cca 12 miliardami let a vyplňují celé haló Galaxie. Obdobně tomu bude asi i u ostatních galaxií. K řešení problému temné hmoty ale bílí trpaslíci zdaleka nestačí. S největší pravděpodobností jde o neznámou formu hmoty nebaryonové povahy. Kvantová teorie Relativita a spin 62 2.7 RELATIVISTICKÁ KVANTOVÁ TEORIE, SPIN 2.7.1 Prostorová rotace a Lorentzova transformace Prostorová rotace Pootočíme-li souřadnicovým systémem kolem osy z o úhel , lze transformaci zapsat jako . ,cossin ,sincos , zz yxy yxx tt       Časovou souřadnici budeme dávat na nultou pozici, při prostorové rotaci se čas nemění. Celou transformaci popíšeme pomocí rotační matice Rz, Podobně můžeme popsat rotace kolem ostatních souřadnicových os (stačí cyklicky zaměnit x  y  z  x): . 1000 0cossin0 0sincos0 0001 , cos0sin0 0100 sin0cos0 0001 , cossin00 sincos00 0010 0001                                                                         zyx RRR Rotace patří mezi unitární transformace. Připomeňme si, že unitární operátory zachovávají skalární součin, proto platí    * 2 det det 1 det det 1 |det | 1.        U U 1 U U U U U Pro reálné matice může být determinant všech unitárních transformací roven buď +1 (rotace) nebo –1 (zrcadlení). Snadno se přesvědčíme, že determinant všech tří rotačních matic je roven jedné. S rotační symetrií se pojí zachování veličiny, kterou nazýváme moment hybnosti. Tato veličina je danou symetrií definována (viz teorém Noetherové, kap. 1.3.1). Lorentzova transformace Velmi příbuznou transformací k rotacím je Lorentzova transformace popisující přechod mezi dvěma vzájemně se rovnoměrně pohybujícími inerciálními souřadnicovými systémy, předpokládejme, že v ose x: . , , /1 , /1 / 22 22 2 zz yy c x x c cxt t         v vt v v Tuto známou transformaci lze zapsat podstatně elegantněji v maticové podobě. Zavedeme-li relativistické proměnné zxyxxxctx  3210 ;;, a relativistické koeficienty 2 1/1;/   cv , Kvantová teorie Relativita a spin 63 budou matice Lorentzovy transformace (v ostatních osách matice získáme cyklickou záměnou) mít tvar . 00 0100 0010 00 , 1000 00 0010 00 , 1000 0100 00 00                                                          zyx  Determinant transformačních matic je roven 1)1(det 22222   a jde tedy opět o rotace, tentokrát v rovině dané časovou a jednou prostorovou osou. Charakter rotací lépe vynikne, zapíšeme-li Lorentzovy matice pomocí tzv. rapidity (je definována vztahem )/(arcth cu v ): ch sh 0 0 ch 0 sh 0 ch 0 0 sh sh ch 0 0 0 1 0 0 0 1 0 0 , , . 0 0 1 0 sh 0 ch 0 0 0 1 0 0 0 0 1 0 0 0 1 sh 0 0 ch x y z u u u u u u u u u u u u                                            S Lorentzovou symetrií (experiment dopadne stejně ve dvou inerciálních soustavách, které se navzájem pohybují rovnoměrně přímočaře) se pojí existence nové zachovávající se veličiny, která se nazývá spin. 2.7.2 Spin V minulé kapitole jsme viděli, že podobnou úlohu, jakou má prostorová rotace má i Lorentzova transformace. Jde také o rotaci, ale v rovině dané časovou a jednou prostorovou souřadnicí o imaginární úhel nazývaný rapidita. Rotační symetrie odpovídá symetrii systému vzhledem k pootočení, Lorentzova symetrie odpovídá stejnému chování systému v různých, navzájem se rovnoměrně pohybujících, inerciálních souřadnicových systémech. S oběma symetriemi se pojí odpovídající zákony zachování: rotační symetrie  moment hybnosti L Lorentzova symetrie  spin S Spin má velmi podobné vlastnosti jako moment hybnosti, lze si ho však jen velmi těžko představit. Značně nepřesné, ale přesto ilustrativní, je představit si částici obíhající kolem centra a současně rotující kolem vlastní osy. V této klasické představě odpovídá momentu hybnosti orbitální rotace a spinu vlastní rotace. Skutečné částice ani neobíhají kolem centra, ani nerotují kolem vlastní osy. Jejich celkový rotační stav je dán dvěma veličinami – momentem hybnosti (orbitálním momentem) a spinem (vnitřním momentem). Obě veličiny se mohou skládat, potom hovoříme o spinorbitální interakci, neboli LS interakci či LS vazbě. Operátor spinu má stejné komutační relace jako moment hybnosti (2.37), (2.39) .0]ˆ,ˆ[ ˆi]ˆ,ˆ[ 3 2 321   SS SSS ,záměnycyklické (2.104) Stejně tak jako u momentu hybnosti zavádíme dvě kvantová čísla popisující spin: spinové číslo neboli spin s určující velikost a magnetické spinové číslo ms určující projekci spinu do L S Kvantová teorie Relativita a spin 64 třetí osy. Pro spin lze pomocí posuvných operátorů odvodit stejně jako pro moment hybnosti vztah (2.76) .,,1,, ;...,23,1,21,0,)1(|| 3 sssmmS sssS ss     (2.105) Tentokrát se ale realizují i poločíselné hodnoty, které jsme pro komutační strukturu (2.104) respektive (2.37) odvodili dříve. Hodnota spinu s je pro elementární částice neměnnou charakteristikou, stejně tak jako hodnota elektrického náboje Q nebo klidové hmotnosti m0. Spin některých částic leptony (elektron, tauon, mion, neutrina) 1/2 kvarky (d,u,s,c,b,t) 1/2 skalární mezony ( kaony) 0 vektorové mezony (,kaony) 1 hadrony (neutron, proton,  hyperon) 1/2 hadrony (, ) 3/2 intermediální bosony (, W , Z0 , gluony) 1 gravitony 2 Přítomnost spinu zvyšuje stupeň degenerace energetických stavů. Například elektron v atomárním obalu, který má energetický stav určený hlavním kvantovým číslem, již nemá stupeň degenerace n2 , ale 2n2 . Elektron má totiž spin 1/2 a jeho stavy jsou určeny čtveřicí čísel n, l, m, ms. Projekce spinu ms může nabývat dvou hodnot 1/2 a počet stavů se zdvojnásobuje. Částice s nenulovým spinem vykazují magnetický moment, aniž by měly orbitální moment hybnosti. Magnetické vlastnosti částic proto nemusí souviset jen se skutečným rotačním pohybem částic, ale i s „vlastním momentem“ – spinem. V přítomnosti nehomogenního magnetického pole reagují částice na toto pole. Stavy, které původně odpovídaly jediné energii, se štěpí na multiplety blízkých energetických podhladin. Stupeň degenerace se snižuje, stavy s různým m a ms mají různou energii. Hovoříme o tzv. sejmutí degenerace v přítomnosti magnetického pole. Pec Magnet StínítkoKolimátor Spin byl poprvé pozorován ve Sternově-Gerlachově experimentu (1925). Atomy stříbra odpařující se z pícky byly kolimovány do svazku procházejícího nehomogenním magnetickým polem. Na tyto elementární magnetické momenty v nehomogenním poli působí síla (1.123) F = B. Magnetický moment jednotlivých stavů je různý a proto je různá i výsledná působící síla a energie daného stavu. Kdyby neexistoval spin, nebude se stav l = 0 štěpit vůbec (m = 0), stav l = 1 se bude štěpit na tři různé podstavy (m = 0, 1) a na stínítku se vytvoří jedna nebo tři stříbrné skvrny (i ve vyšších stavech l půjde vždy o lichý počet skvrn). Na stínítku však byly pozorovány dvě stříbrné skvrny, což svědčí o elektronu s orbitálním stavem l = 0 a spinovým stavem s = 1/2 (magnetické vlastnosti jsou určeny dvěma projekcemi Kvantová teorie Relativita a spin 65 ms = 1/2). Sudý počet projekcí znamená poločíselné řešení komutačních relací (2.104) respektive (2.37). Hypotézu o existenci vlastního momentu elektronu, který má podobné vlastnosti jako orbitální moment, podali ještě před teoretickým objasněním spinu Uhlenbeck a Goudsmit v roce 1925. Na následujícím obrázku je numerická simulace (Yamanashi University) Stern Gerlachova experimentu. Stavy s projekcí ms = +1/2 jsou označeny modře, stavy ms = 1/2 červeně. V malé vzdálenosti se na stínítku objeví dvě výrazné stříbrné skvrny, ve větší vzdálenosti nejsou pravděpodobnosti dopadu atomů v jednotlivých stavech výrazně prostorově oddělené. 2.7.3 Kleinova-Gordonova rovnice Schrödingerova rovnice není relativistická, a proto nemůže správně popsat spin. Při jejím odvození jsme používali nerelativistický tvar Hamiltonovy funkce. Výsledkem byla Schrödingerova časová rovnice (2.98), která má v x reprezentaci tvar 0 2 i )( 2      V mt   . V rovnici se nachází první časová derivace a druhé prostorové derivace, čas a prostor není rovnoprávný, rovnice zjevně není relativistická. Relativistickou konstrukci lze vytvořit jak ve druhých (Kleinova-Gordonova rovnice), tak v prvních (Diracova rovnice) derivacích. V této kapitole se budeme zabývat konstrukcí správné rovnice ve druhých derivacích. ■ Kleinova-Gordonova rovnice Rovnici poprvé odvodili Oskar Klein a Walter Gordon. Předpokládejme, že hledáme lineární rovnici, která limitně při malých rychlostech přejde v Schrödingerovu rovnici. U lineárních rovnic platí princip superpozice a obecné řešení lze vždy složit z rovinných vln i[ ]i[ ] i[ ] ( ) ( )e ( )e ( , )ek xk x t k x a k a k a           k x k . (2.106) Třírozměrné vektory jsou označeny tučně. Složky vlnového vektoru kα musí být nutně závislé, neboť i parciální vlny (2.106) musí splňovat hledanou rovnici. Taková závislost se nazývá disperzní relace a můžeme ji zapsat v implicitním tvaru ( , ) 0  k . (2.107) V některých případech je možné nalézt explicitní závislost ω = ω(k). Obecná vlnová funkce bude superpozicí i[ ] 4 i[ ] 3 ( ) ( )e ( ) ( , )ek x t x a k d k a d          k x k k . (2.108) Diracova distribuce zajišťuje automatické splnění disperzní relace (2.107). Parciální (rovinné) vlny lze snadno derivovat: ( ) i ( )k kx k x     (2.109) a parciálním derivacím odpovídají algebraické výrazy Kvantová teorie Relativita a spin 66 ik    (2.110) S využitím duality (2.2) máme i p    . (2.111) Nejpřirozenějším přechodem od komutujícího k nekomutujícímu popisu je tedy zavedení operátorů na L2 předpisem ˆ i ; ˆ . p x x          . (2.112) Poznámka 1: Snadno dopočteme, že takto zavedené operátory splňují komutační relace, které jsou ve shodě s principem korespondence mezi Poissonovými závorkami a komutátory ˆ ˆ ˆ ˆ, , 0 , ˆ ˆ, i . p p x x x p g                       (2.113) Poznámka 2: V (3+1) D formalizmu lze první ze vztahů (2.112) zapsat jako ˆ i , ˆ i . E t       p x   (2.114) Odlišné znaménko u časové proměnné souvisí s relativistickými transformačními vlastnostmi čtyřvektorů. Druhou relaci jsme již používali v x reprezentaci operátoru hybnosti, viz (2.47). Najděme nyní velikost čtyřhybnosti 2 2 2 / . E c E p p p c              p p (2.115) tato hodnota musí být ve všech souřadnicových soustavách stejná a můžeme ji určit v klidové soustavě částice, kde je E = m0c2 , p = 0: 2 2 0p p m c    . (2.116) V (3+1) formalizmu jde o známou Pythagorovu větu pro energii 2 2 2 2 4 .E c m c p Tento vztah je správným relativistickým vztahem pro energii volné částice, a proto se o něho musí opírat odvození relativistické varianty Schrödingerovy rovnice. Přepišme proto (2.116) do operátorové podoby:  2 2 0 0ˆ ˆ ˆ0 ; ip p m c p         . (2.117) Rovnice (2.117) je Kleinova-Gordonova rovnice pro volnou částici. Po dosazení za operátory získáme jiný často používaný tvar Kleinovy-Gordonovy rovnice  2 0 0 ; m c        . (2.118) Kleinova-Gordonova rovnice je relativistickou analogií Schrödingerovy rovnice pro volnou částici. Při malých rychlostech limitně přechází v nerelativistickou Schrödingerovu rovnici. Jde o lineární rovnici a každé její „rozumné“ řešení je možné zapsat pomocí Fourierovy transformace jako superpozici rovinných vln. Konstanta κ je v normální soustavě jednotek (c = 1, ħ = 1) rovna klidové hmotnosti částice. ■ Nerelativistická limita Kleinovu Gordonovu rovnici (2.117) můžeme v operátorovém tvaru zapsat jako Kvantová teorie Relativita a spin 67 2 2 2 2 4 0 ˆˆ ˆE c m c p 1. (2.119) Obě strany formálně odmocníme. Odmocninu chápeme jako funkci operátoru ve smyslu (2.8) nebo (2.29): 2 2 2 2 2 4 2 2 0 0 02 2 2 2 0 0 2 2 0 0 ˆ ˆˆ ˆ ˆˆ ˆ 2 ˆˆˆ 2 E c m c m c m c m c m c E m c m                  p p p 1 1 1 p 1  . Záporné znaménko před odmocninou jsme zatím vynechali jako nefyzikální a budeme se jím zabývat až v kapitole věnované Diracově rovnici. První člen můžeme chápat jako konstantní/nulovou potenciální energii (posunutím o konstantu se potenciální energie nezmění) a druhý je běžná kinetická energie částice. Po dosazení za operátory z (2.114) získáme časovou Schrödingerovu rovnici (2.98) s nulovou, resp. konstantní potenciální energií. Pro malé rychlosti (hybnosti) Kleinova-Gordonova rovnice přechází ve Schrödingerovu rovnici. ■ Pravděpodobnostní interpretace Hustota ρ a tok pravděpodobnosti j výskytu částice by měly splňovat rovnici kontinuity (zákon zachování pravděpodobnosti výskytu částice) ve tvaru / 0 ; c j j           j . (2.120) Ukažme, že takový zákon zachování je v Kleinově-Gordonově rovnici obsažen. Nalezněme kombinaci ψ* (2.118)–ψ(2.118)* :    * 2 2 * * * * * 0 0 0 .                                    Nyní v obou výrazech využijeme identitu f ∂αg = ∂α(f g)– (∂αf )g:          * * * * 0                          . Pokud v posledním výrazu zvýším první index a snížíme druhý, vyruší se s druhým výrazem:    * * 0             * * 0 ;j j              . (2.121) Čtyřvektor jα reprezentuje nenormovanou pravděpodobnost výskytu částice. Hustota pravděpodobnosti j0 (v SI j0 /c) není bohužel pozitivně definitní a Kleinova Gordonova rovnice připouští i záporné hustoty pravděpodobnosti. Řešením tohoto problému (vyústí v existenci antičástic) se budeme zabývat v kapitole věnované Diracově rovnici. ■ Disperzní relace Po dosazení rovinné vlny (2.106) do Kleinovy-Gordonovy rovnice získáme disperzní relaci 2 2 2 2 2 2 2 2 2 .c k c c k c         (2.122) Záporné řešení (odpovídá záporné energii ħω) budeme opět považovat za nefyzikální. Standardním postupem určíme fázovou a grupovou rychlost: Kvantová teorie Relativita a spin 68 2 2 2 f 2 2 g 2 2 2 2 2 1 1 , 4 . 1 1 4 c c k k c c k k                       v v Na první pohled je zřejmé, že grupová rychlost je vždy podsvětelná. Z Hamiltonových rovnic mechaniky g H x k k p               v plyne, že grupová rychlost vlnového balíku je analogem mechanické rychlosti pohybující se částice. Oproti tomu fázová rychlost je vždy nadsvětelná a nemá význam přenosu informace. Mezi oběma rychlostmi je jednoduchý vztah vfvg = c2 . Obě rychlosti závisí na vlnové délce parciální vlny, tj. dochází k disperzi. ■ Kleinova-Gordonova rovnice pro nabitou částici v elektromagnetickém poli V přítomnosti elektromagnetického pole se v Hamiltonově funkci (1.90) vyskytovala kanonická (zobecněná) hybnost v kombinaci p – QA. Obdobně tomu musí být i v Kleinově Gordonově rovnici (2.117), která má pro nabitou částici v elektromagnetickém poli tvar    2 2 0 0ˆ ˆ ˆ0 ; ip QA p QA m c p                 . (2.123) Po dosazení za operátor hybnosti a roznásobení všech členů máme 2 2 2 2 0 i 2i 0Q A A m c Q A QA                   . (2.124) Využijeme-li kalibrační podmínku (1.188), tj. položíme-li ∂αAα = 0, získáme výslednou rovnici 2 2 0 2 2i 0 ; m cQ Q A A A                    (2.125) pro popis nabité částice v elektromagnetickém poli. ■ Vodíkový atom Naznačme nyní, jak by se postupovalo při hledání spektra vodíkového atomu z KleinovyGordonovy rovnice. Elektron s nábojem Q = –e je v poli jádra, které lze vyjádřit vztahy 0 0 , 4 0 . Ze A c rc      A (2.126) Kleinova-Gordonova rovnice získá tvar 22 2 2 2 1 2i ( , , , ) 0 . e e t r c c c tc t                         Laplaceův operátor rozložíme na radiální a úhlovou část stejně jako v nerelativistickém případě (2.80). Budeme hledat stacionární řešení, tj. časovou část vlnové funkce budeme předpokládat ve tvaru exp(–iωt) = exp(–iEt/ħ), prostorovou část zapíšeme jako součin radiální a úhlové části, kterou již známe z dřívějška, viz (2.85): Kvantová teorie Relativita a spin 69 2 2 2 2 2 2 2 2 2 2 ˆ 1 2i e ( ) ( , ) 0 , i Et r lm Z Z R r Y r c tr c t r                      L kde jsme označili 04 e c     (2.127) tzv. konstantu jemné struktury. Po provedení časových derivací získáme rovnici 22 2 42 0 2 2 2 2 ˆ 2 ( ) ( , ) 0 .r lm E m c Z Z E R r Y r r cr c                      L Nyní zapůsobíme operátorem 2ˆL na úhlovou část vlnové funkce podle vztahu (2.75) a vyjádříme radiální část Laplaceova operátoru 2 2 42 2 2 0 2 2 2 2 1 ( 1) 2 ( ) 0 . E m cd d l l Z Z E r R r dr dr r cr r c              Jde o obyčejnou diferenciální rovnici, která se řeší standardními postupy (asymptotické chování, rozvoj do řady, oříznutí). Výsledkem jsou tzv. Laguerrovy polynomy a energetické spektrum   2 2 4 4 2 2 6 0 0 2 4 3 (2 1) ( ) 82 (2 1) nl Z Z E m c m c n l Z n l n                    O . (2.128) Hlavní kvantové číslo je definováno stejně jako v nerelativistickém případě, druhý člen v hranaté závorce reprezentuje první relativistickou korekci a současně sejmutí degenerace spektrálních čar. ■ Problémy Kleinova-Gordonova rovnice má tři základní problémy: 1. Druhé časové derivace znamenají zadání počáteční podmínky nejen na vlnovou funkci (reprezentuje stav systému), ale i na první časovou derivaci vlnové funkce, což je fyzikálně jen obtížně interpretovatelné. 2. Hustota pravděpodobnosti není pozitivně definitní. 3. Kleinova-Gordonova rovnice poskytuje i záporné energetické stavy. 2.7.4 Diracova rovnice Správnou relativistickou kvantovou rovnici pro nabitou částici, ve které jsou obsaženy jen první derivace, odvodil Paul Adrien Maurice Dirac (1902–1984) v roce 1928. Ukázalo se, že jde o mnohem vhodnější rovnici pro elektron, než je Kleinova-Gordonova rovnice. Tím, že rovnice je jen v prvních derivacích, postačí zadat počáteční hodnotu vlnové funkce a automaticky odpadá nutnost zadávat první derivaci vlnové funkce. U Diracovy rovnice je hustota pravděpodobnosti pozitivně definitní a tak odpadá i druhý základní problém Kleinovy-Gordonovy rovnice. Problém záporných energetických stavů nicméně přetrvává a Dirac tyto stavy interpretoval jako stavy příslušející antičástici k elektronu – pozitronu. Ten byl objeven až v roce 1932 Carlem Andersonem. ■ Diracova rovnice Hledejme rovnici, která má stejný tvar jako Schrödingerova rovnice, ale Hamiltonův operátor je lineární funkcí prostorových derivací: Kvantová teorie Relativita a spin 70 1 2 3 1 2 3 ˆ ˆ; .i H H a a a b t              (2.129) Z rozměrových důvodů budeme namísto koeficientů ak a b hledat koeficienty αk a β, které jsou bezrozměrné:  1 2 3 2 1 2 3 0 ˆ iH c m c           . (2.130) Na koeficienty máme dvě základní podmínky 1. Kvadrát Hamiltonovy funkce musí dát pravou stranu (2.119), tj. 2 2 2 2 4 0 ˆˆ ˆ ;H c m c p 1 (2.131) tím bude každé řešení Diracovy rovnice řešením Kleinovy-Gordonovy rovnice (nikoli naopak, druhé derivace některá řešení přidají). 2. Nová rovnice musí být relativisticky kovariantní (tj. její tvar se nesmí změnit po provedení Lorentzovy transformace souřadnic a polí). Za chvíli uvidíme, že tyto podmínky nesplňují žádné číselné koeficienty a hledaná čísla αk a β musí být matice. Vyjděme z podmínky (2.131), do které dosadíme hamiltonián (2.130) a za operátor hybnosti z (2.114):    2 2 2 2 4 0 2 2 2 2 2 4 0 0 0 ˆˆ ˆ ; ˆˆi i ;k l k l H c m c c m c c m c c m c              p 1 p 1   2 2 2 2 2 4 2 2 2 4 0 0 0 ˆik l k k k l kc cm c m c c m c                1   Porovnáním členů na levé a pravé straně máme na koeficienty : 2 , 0 , ˆ . k l k l k k              1 (2.132) První relaci upravíme snadno na tvar   0 pro ,1 2 pro .2 k l l k k l l k k l k l k l                    (2.133) Požadavky (2.132), resp. (2.133) nesplňují žádná reálná ani komplexní čísla. Budeme proto hledat soustavu čtyř matic, jejichž zajímavé vlastnosti nejprve přehledně sepíšeme a vzápětí dokážeme 1. Matice αk a β antikomutují (každá s každou):    , , 0 ;k l k k l      . (2.134) 2. Kvadráty matic αk a β dají jednotkovou matici:         2 2 2 21 2 3 ˆ .       1 (2.135) 3. Matice αk a β jsou hermitovské:   † † , .k k      (2.136) 4. Vlastní čísla matic αk a β mohou nabývat jen hodnot +1 a –1. 5. Stopa matic αk a β je nulová. 6. Matice αk a β jsou nezávislé. Kvantová teorie Relativita a spin 71 Dokažme nyní jednotlivá tvrzení Ad 1) Antikomutační relace matic αk a β plynou okamžitě z relací (2.132) a (2.133). Poznamenejme, že antikomutátor dvou objektů je definován jako {A,B}≡AB+BA. Ad 2) Tvrzení opět plyne okamžitě z relací (2.132) a (2.133). Ad 3) Hermitovost matic αk a β plyne z požadavku na hermitovost operátoru energie (2.130). Ad 4) Z podmínky (2.135) plyne, že vlastní čísla matic αk a β leží na jednotkové kružnici v komplexní rovině, tj. | λ | = 1. Hermitovské matice ale mají reálná vlastní čísla, tedy připadají v úvahu pouze hodnoty λ = ±1. Ad 5) Stopou matice nazýváme součet diagonálních členů Tr( ) k kAA . (2.137) Tr je zkratkou z anglického trace. Stopa matice se nezmění při cyklické záměně matic: 1 2 2 1Tr( ) Tr( )N NA A A A A A  , (2.138) tj. první matici můžeme odstěhovat na poslední místo v součinu (nebo poslední na první). Nyní již snadno dokážeme, že stopa hledaných matic je nulová:            2 Tr Tr Tr Tr Tr Trk k k k k k               . Nejprve jsme přidali β2 , což je ale jednotková matice. Poté jsme jednu matici β odstěhovali na konec za pomoci cyklické záměny a vrátili ji zpět na původní pozici s využitím antikomutativnosti matic α a β. Přečteme-li si začátek a konec, máme        Tr Tr 2Tr 0 Tr 0k k k k          . Obdobně můžeme postupovat u matice β:            2 Tr Tr ( ) Tr Tr Tr Tr 0k k k k k                  . Ad 6) Předpokládejme závislost matic, tj. například matici β bude možné vyjádřit jako lineární kombinaci ostatních: k kc   . Vynásobme relaci zleva maticí β: 1 1 Tr( ) Tr( ) Tr( ) Tr(0) 0 . 2 2 k k k k k k k k c c c c                 1 1 Jde o spor, neboť stopa jednotkové matice nalevo je nenulová. Matice tedy musí být nezávislé. Tím jsou všechna tvrzení (1 až 6) dokázána. ♦ Stopa matic je invariantem, tj. ve všech bázích/souřadnicových soustavách je stejná. Pokud u hermitovské matice za bázi zvolíme její vektory, bude matice diagonální a na diagonále budou její vlastní čísla. Stopa matice je proto součtem vlastních čísle matice. V našem případě jsou vlastní čísla +1 nebo –1, stopa matice je nulová, a proto musí mít hledané matice sudou dimenzi (aby součet čísel +1 a –1 mohl dát nulu). Kvantová teorie Relativita a spin 72 N = 2 V kapitole věnované momentu hybnosti jsme odvodili tzv. spinorovou reprezentaci momentu hybnosti (2.77). Matice spinu bez příslušných koeficientů se nazývají Pauliho matice: 1 2 30 1 0 i 1 0 ; ; . 1 0 i 0 0 1                            (2.139) Pauliho matice mají námi hledané vlastnosti. Jsou hermitovské, antikomutují mezi sebou, jejich kvadráty jsou jednotkové matice, vlastní čísla jsou +1 a –1, součet členů na diagonále je nulový. Jejich jedinou nevýhodou je, že jsou jen tři. My hledáme soustavu čtyř nezávislých antikomutujících matic. Ve dvou dimenzích taková soustava ale neexistuje. Další nezávislou maticí k Pauliho maticím je jednotková matice, ale ta s nimi komutuje, nikoli antikomutuje. Navíc u ní není součet diagonálních členů nulový. N = 4 Ve čtyřech dimenzích existuje celkem 16 nezávislých matic a skutečně z nich lze vybrat 4 antikomutující matice požadovaných vlastností. Jde o nejmenší počet dimenzí, ve kterých lze vyřešit Diracovu úlohu. Existuje více způsobů, jak vybrat hledanou soustavu antikomutujících matic. Dirac je blokově skládal z Pauliho matic a nalezl řešení 3 1 ; k k k k                         1 0 0 1 0 1 0 . (2.140) Každý prvek matice znamená blok 2×2. Výsledné Diracovy matice tedy jsou: 1 2 3 1 0 0 0 0 0 0 1 0 1 0 0 0 0 1 0 ; ; 0 0 –1 0 0 1 0 0 0 0 0 –1 1 0 0 0 0 0 0 i 0 0 1 0 0 0 i 0 0 0 0 –1 ; . 0 i 0 0 1 0 0 0 i 0 0 0 0 –1 0 0                                                           (2.141) Ověřte si, že všechny matice jsou hermitovské, mají vlastní čísla +1 a –1, součet prvků na diagonále je 0, v kvadrátu dají jednotkovou matici a každá matice antikomutuje s každou. Diracova rovnice pro volnou částici má nyní jednoduchý tvar:   1 22 0 3 4 ( , ) ( , ) i i ; ( , ) ( , ) k k t t c m c tt t                          x x x x   . (2.142) Koeficienty rovnice jsou matice 4×4, vlnovou funkci proto tvoří čtveřice funkcí (nejde o čtyřvektor!). Jiná volba čtveřice Diracových matic by vedla na tatáž fyzikální řešení. ■ Operátor rychlosti, záporné energie Určeme operátor rychlosti částice jako operátor časového vývoje polohy podle principu korespondence (2.96): 2 0 1 1 1ˆˆ , , i , i i i i , i , i . i i i k k k k l k l l l l k l k l k k l l l dx x H x c m c x c dt c c c x x p c                                                     v Kvantová teorie Relativita a spin 73 Matice αk tak mají (až na konstantu c) význam operátoru rychlosti: ˆk k cv . (2.143) formálně lze zapsat všechny tři relace naráz ˆˆ cv   . (2.144) Za pomoci operátorů rychlosti a hybnosti získá Diracova rovnice (2.142) jednoduchý tvar:  2 0 ˆˆi ; ˆ ˆˆ , i . m c t c             v p v p      (2.145) Řešme nyní Diracovu rovnici pro částici v klidu, tj. s nulovým operátorem rychlosti 1 1 2 22 0 3 3 4 4 i .m c t                                   Řešení je: 2 0 1 1 2 0 2 2 2 0 3 3 2 0 4 4 ( , ) ( )exp i , ( , ) ( )exp i , ( , ) ( )exp i , ( , ) ( )exp i . m c t A t m c t A t m c t A t m c t A t                                     x x x x x x x x     Porovnáme-li řešení s časovou částí rovinné vlny exp[–iωt] = exp[–i(E/ħ)t], je zřejmé,že první dvě řešení odpovídají kladné energii E = m0c2 a druhá dvě záporné energii E = – m0c2 . Problém záporných energetických stavů tak Diracova rovnice nevyřešila. Ukázalo se, že Diracova rovnice popisuje chování částic se spinem 1/2 (například elektron). Čtveřice  se nazývá bispinor. Má speciální transformační vlastnosti. Horní dvě komponenty bispinoru popisují stavy částice s projekcí spinu +1/2 a –1/2 a mají kladnou energii. Dolní dvě komponenty mají zápornou energii a Dirac je interpretoval jako stavy antičástice s projekcí spinu +1/2 a –1/2. Diracova rovnice je v jistém smyslu „odmocněním“ KleinovyGordonovy rovnice postavené na vztahu E2 = p2 c2 + m0 2 c4 . Proto stavy se zápornou energií nejsou překvapením. Elegantní však bylo Diracovo vysvětlení: Všechny záporné stavy jsou zaplněny (Diracovo moře elektronů se zápornou energií). Nezaplněný stav se chová jako „díra“, kterou Dirac interpretoval jako antičástici s kladnou energií. Rozborem tvaru rovnice a jejích řešení teoreticky Dirac v roce 1928 předpověděl existenci pozitronu ještě před jeho experimentálním objevem v roce 1932 (Carl Anderson). ■ Pravděpodobnostní interpretace Při odvození rovnice kontinuity pro pravděpodobnost budeme postupovat stejně jako u Kleinovy-Gordonovy rovnice, jen namísto komplexního sdružení budeme využívat hermitovské sdružení jednotlivých matic i základního bispinoru, který tvoří vlnovou funkci. Hermitovsky sdružený bispinor má tvar Kvantová teorie Relativita a spin 74  † * * * * 1 2 3 4     . (2.146) Nalezněme nyní kombinaci ψ† (2.142) – ψ (2.142)† :       † † †† † 2 2 0 0i ( i ) i ik k k kc m c c m c t t                                   ,   † † † 2 † † 2 † 0 0i i i ik k k kc m c c m c t t                               ,   † † † †k k k kc c t t                     ,    † † k k c t          ,    † † 0 .k kc c ct           Získali jsme tak rovnici kontinuity ve tvaru 0 † † 0 ; ˆˆ ˆ, , . j j c c           j v v     (2.147) Hustota pravděpodobnosti je dána vztahem 0 † * * * * P 1 1 2 2 3 3 4 4/ 0j c                 (2.148) a je tedy pozitivně definitní. Tok pravděpodobnosti je zobecněním vztahu pro klasický tok (hustota × rychlost), rychlost nahrazuje operátor rychlosti. Výsledný tok pravděpodobnosti je ale obyčejným vektorem, neboť každá z jeho komponent je součinem řádkové, čtvercové a sloupcové matice, tj. dá obyčejné číslo. ■ Diracova rovnice pro nabitou částici v elektromagnetickém poli Zobecnění z volné částice na částici v poli provedeme stejně jako u Kleinovy-Gordonovy rovnice, tj. nahradíme ˆ ˆp p QA     . (2.149) V (3+1) symbolice máme ; . i i Q ct ct c i i Q               A x x     (2.150) Diracova rovnice (2.145) získá nyní tvar   2 0 ˆˆi Q Q m c t              v p A   neboli 2 0 ˆˆ ˆi m c Q Q t                  v p A v    . (2.151) Oproti volné částici přibyl napravo interakční hamiltonián podle vztahu 0 I 2 0 0 I ˆ ˆi ; ˆˆ ˆˆ ˆ, . H H t H m c H Q Q               v p A v    (2.152) Kvantová teorie Relativita a spin 75 I v tomto případě jde jen o přímé zobecnění interakčního členu známého z Lagrangeovy funkce v klasické mechanice. ■ Kovariantní tvar Diracovy rovnice Přenásobme Diracovu rovnici (2.142) zleva maticí β a poté převeďme všechny členy doleva:  2 0 0 i i i i 0 k k k k c m c t m c ct                             . Získali jsme tak nejznámější tvar Diracovy rovnice 0 0 i 0 ; , ( ) k k m c             , (2.153) ve které jsou koeficienty derivací tzv. gama matice 0 1 2 3 1 0 0 0 0 0 0 1 0 1 0 0 0 0 1 0 ; ; 0 0 1 0 0 1 0 0 0 0 0 1 1 0 0 0 0 0 0 i 0 0 1 0 0 0 i 0 0 0 0 1 ; ; 0 i 0 0 1 0 0 0 i 0 0 0 0 1 0 0                                                         (2.154) Mezi původními maticemi a maticemi gama existují jednoduché transformační vztahy: 0 , , . k k k k          (2.155) Důkaz třetího vztahu je triviální, druhý vztah stačí násobit maticí β zleva. Pro prostorovou část obou sad matic tak platí jednoduché pravidlo: násobením maticí β zleva dostaneme odpovídající matici z druhé sady. Matice gama opět antikomutují, nejsou již ale hermitovské a kvadráty prostorových matic nedají jednotkovou matici, ale minus jednotkovou matici:   21 1 1 1 1 1 1 1 1                1. Obdobně bychom postupovali u ostatních matic, platí tedy   2 ; 1,2,3k k   1 . (2.156) Zaveďme nyní dvě užitečné a často používané operace. První z nich je Diracovo sdružení: † 0 A A . (2.157) Jde o hermitovské sdružení doplněné násobením maticí γ0 zprava. Druhou operací je Feynmanovo zúžení (Feynman slash): K K   . (2.158) Za pomoci těchto operací lze elegantně zapsat složky čtyřtoku pravděpodobnosti (2.147) Kvantová teorie Relativita a spin 76 0 † † 0 0 0 † † 0 , .k k k k j c c c j c c c                        Jednotně tedy můžeme psát 0 ;j j c        . (2.159) Diracovu rovnice (2.153) lze přepsat do „úsporného“ tvaru  0ˆ 0p m c    , neboli (2.160)  0ˆ 0p m c   . (2.161) Pro nabitou částici v elektromagnetickém poli bude mít Diracova rovnice nyní velmi jednoduchý tvar  0ˆ A 0p Q m c    (2.162) Spin se stal automatickou součástí relativistických rovnic kvantové teorie. Rovnice KleinovaGordonova se nakonec ukázala vhodnou rovnicí pro skalární částice (se spinem nula), rovnice Diracova pro částice se spinem ½ (elektrony, neutrina, kvarky). Právě na ní je postavena dnešní kvantová elektrodynamika. Příklad 17:     0 ; , 2; 0 , 2 2; 1,2,3 g                           (2.163) Kvadrát matic α a β je roven jednotkové matici. U matic γ tomu tak není, (γ0 )2 =1, ale (γk )2 = −1 pro k = 1,2,3. Je to přirozené, v Minkowského metrice se prostorová část chová vždy jinak než časová část. Ve vztahu (2.163) je výsledek na pravé straně vždy násoben jednotkovou maticí, tu ale nebývá zvykem psát. Příklad 18:   2 2 0 1 , 2 KK K K K K K K g K K K                            K . Příklad 19:  1 , 2 g                                   Příklad 20: 0 0 ( / ) ( / ) kk kk k k kk k E c pE p p p p p c p E c                               11 0 0 0 1 10 . Příklad 21: Dokažme, že platí relace (užitečná při výpočtu sdružených matic γμ† ) † 0 0      . (2.164) Relaci (2.164) zleva a zprava vynásobíme γ0 a dokážeme platnost vztahu γ0 γμ† γ0 = γμ : Sdruženou matici γμ† rozepíšeme z (2.155) za pomocí sady hermitovských matic α a β:   †0 0† 0 0 0 0 0 0 0 0 0 † 0 0 0 0 0 0 0 0 ; .k k k k k                               Kvantová teorie Relativita a spin 77 ■ Matice C Existuje řada dalších zajímavých matic, které lze odvodit ze základní sady matic γμ . Zaveďme nejprve C matici, která bude užitečná při nábojovém sdružení (přechodu od částic k antičásticím a také při výpočtu transponovaných matic γμT : 2 0 0 0 0 1 0 0 1 0 i 0 1 0 0 1 0 0 0                C . (2.165) Matice má prvky jen na vedlejší diagonále, a to střídající se hodnoty +1 a –1. Stopa matice je nulová. Na první pohled je zřejmé, že pro C platí zajímavá vlastnost: T 1     C C C C. (2.166) Chceme-li tedy najít inverzní matici, transponovanou či hermitovsky sdruženou, stačí jen změnit znaménko matice (vyměnit pořadí +1 a –1 na vedlejší diagonále). Pokud potřebujeme nalézt transponovanou matici γμT , lze k tomu využít matice C: T T 1 ; resp.          C C C C (2.167) Tvrzení se dokáže pouhým převedením matic γμ na sadu matic α a β, které se po transpozici nezmění. ■ Matice γ5 Další důležitou maticí, která má využití při popisu levopravé symetrie, je matice 5 0 1 2 3 0 0 1 0 0 0 0 1 i . 1 0 0 0 0 1 0 0                        0 1 1 0 (2.168) Tato matice je hermitovská, její kvadrát je roven jednotkové matici, je lineárně nezávislá na ostatních γ maticích a antikomutuje s nimi:     5† 5 2 , , , 0 ; 0,1,2,3.             1 (2.169) ■ Matice Σ a báze Γk Z matic γμ můžeme zkonstruovat ještě matice Σ, které jsou užitečné při definici báze na prostoru matic a při hledání transformačních vlastností antisymetrických tenzorů. Matice definujeme jako komutátory , 2 i         . (2.170) Zřejmě platí 0 ; , ; .              (2.171) Je zjevné, že existuje celkem šest nezávislých matic Σαβ , například Σ01 , Σ02 , Σ03 , Σ12 , Σ13 , Σ23 . Ostatní prvky jsou buď nulové nebo je lze dopočíst z antisymetrie. Zajímavou bází na prostoru matic 4×4 je následujících 16 matic: Kvantová teorie Relativita a spin 78  5 5 , , , ,k          1 . (2.172) Celkem snadno lze ukázat, že matice Γk jsou nezávislé. Protože jich je 16, tvoří bázi na prostoru matic 4×4. Jejich kvadrát je +1 nebo –1, vhodným vynásobením ±i by bylo možné docílit, aby kvadrát byl vždy roven jednotkové matici, je to však zbytečné. Stopa všech, s výjimkou jednotkové matice Γ1 , je nulová. Násobek libovolných dvou různých matic Γk je až na znaménko roven některé další matici Γ:   2 1) jsou lineárně nezávislé, 2) 1, 4 1; 3) Tr( ) 0 1, 4) , ; , : . k k k k l m k k k l k l m k l                   (2.173) Z rovnice kontinuity (2.159) je zřejmé, že veličina    je čtyřvektor. Obdobně lze pomocí členů báze Γk zkonstruovat i další veličiny, které se transformují charakteristickým způsobem: 5 5 skalár pseudoskalár vektor pseudovektor antisymetrický tenzor              (2.174) Kvantová teorie Relativita a spin 79 2.7.5 Pozitron, C symetrie Nejprve popíšeme úvahy, které vedly Paula Diraca k předpovědi existence pozitronu a poté matematickou transformaci (nábojové sdružení), která převede Diracovu rovnici pro elektron na rovnici pro pozitron. Dirac předpověděl existenci pozitronu v roce 1928, objeven byl Carlem Andersonem v kosmickém záření v roce 1932. ■ Diracovo moře Jak jsme viděli, z Diracovy rovnice vycházejí záporné energetické stavy. Záporné energie se ale v přírodě nevyskytují a tak Diraca napadlo, že tyto stavy jsou všechny zaplněny elektrony a žádný z nich není volný, proto je nepozorujeme. Vakuum je podle této představy tvořeno mořem elektronů v záporných energetických stavech, tzv. Diracovým mořem. 0 +m c0 2 +m c0 2 +m c0 2 –m c0 2 –m c0 2–m c0 2 0 E E EE Diracovo moře záporných stavů přilétá foton s dostateč ou energiín foton vyráží elektron a zůstává díra díra se vůči okolí chová jako kladně nabitá částice s kladnou energií 0 Představme si, že do tohoto moře vletí foton s energií větší než je dvojnásobek klidové energie elektronu. Potom může z Diracova moře vytrhnout elektron a má-li dostatečnou energii, převést ho do některého energetického stavu s kladnou energií. V záporném Diracově moři zůstane díra – prázdný energetický stav, který se vůči okolí jeví jako kladně nabitá oblast s kladnou energií (hmotností). Dirac tuto díru v roce 1928 interpretoval jako kladně nabitou částici, která má jinak shodné vlastnosti s elektronem, a nazval ji pozitron. V roce 1929 tento koncept rozšířil na všechny částice a zavedl pojem antičástice – objektu, který má opačné hodnoty všech kvantových nábojů oproti původní částici. Navenek se tedy zdá, jakoby se původní foton rozpadl na elektron-pozitronový pár. Pohyb volného elektronu by měl být Diracovým mořem ovlivněn. Elektron interaguje s blízkými elektrony v záporných stavech, odtlačuje je od sebe a z dálky vypadá, jako by měl menší náboj, než skutečně má. Z větší vzdálenosti proto nevidíme skutečný náboj elektronu, ale náboj odstíněný Diracovým mořem. Čím více se přibližujeme k letícímu elektronu, tím více vnímáme jeho skutečný, holý náboj. ■ Nábojové sdružení Diracova rovnice pro elektron ve vnějším poli má tvar (2.162)       0 0 0 ˆ A 0 0 0 p Q m c i Q A m c i Q A m c                                  0 0i e A m c         (2.175) Rovnice pro pozitron by měla mít tvar Kvantová teorie Relativita a spin 80  0 C 0i e A m c         , (2.176) kde ψC je vlnová funkce pozitronového řešení. Diracovu rovnici nejprve hermitovsky sdružíme a poté ji transponujeme. Po těchto dvou operacích přejde rovnice pro elektron v rovnici pro pozitron. Vzpomeňte si, že operace hermitovského sdružení a transpozice splňují vlastnost (2.10):    † T† † T T ; . AB B A AB B A (2.177) Proveďme tedy Hermitovo sdružení Diracovy rovnice pro elektron (2.175): † † 0 0i e A m c                  . Derivace nyní působí samozřejmě vlevo, tj. na vlnovou funkci. Hermitovo sdružení nyní aplikujme na jednotlivé členy v závorce † † † 0i 0e A m c                   a sdružené matice gama vyjádříme ze vztahu (2.164) † 0 0 0 0 0i 0e A m c                       . Rovnici vynásobíme zprava maticí γ0 : † 0 0 0 i 0 i 0 e A m c e A m c                                           . Po operaci hermitovského sdružení se změnilo znaménko prvého členu, vlnová funkce je nalevo a má tvar diracovsky sdruženého bispinoru. Nyní vrátíme vlnovou funkci doprava za pomoci operace transpozice     T T 0i 0e A m c          . Provedeme transpozice všech členů v závorce   TT T 0i 0e A m c          a transponované matice gama vyjádříme ze vztahu (2.167)   T 0i 0e A m c         C C C C 1 . Celou rovnici vynásobíme maticí C−1 zleva   T1 0i 0e A m c          C C C a využijeme relace (2.166) pro inverzní matici C−1 = −C:    T 0i 0e A m c         C Transpozice tedy změnila znaménko posledního členu rovnice a výsledek je Kvantová teorie Relativita a spin 81     0 T C i 0 ;C . e A m c             C  (2.178) Získali jsme hledanou rovnici pro pozitron, který má opačný náboj a nezměněnou hmotnost. Pokud známe v dané situaci řešení ψ pro elektron, bude ve stejné situaci odpovídajícím řešením pro pozitron vlnová funkce ψC = C( )T . Řešení pro pozitron tedy není novým řešením, je obsaženo v řešení pro elektron. Diracova rovnice nezmění svůj tvar (tzv. kovariance), provedeme-li nábojové sdružení neboli C transformaci  T ; . A A       C (2.179) V původní Diracově rovnici odpovídají dvě řešení s kladnou energií elektronu s projekcí spinu +1/2 a −1/2 (proto je řešení zdvojené) a řešení se zápornou energií pozitronu s projekcí spinu +1/2 a −1/2. Dvě dvojice spojené do jedné čtveřice vlnových funkcí se nazývají bispinor. Po provedení transformace (2.179) mají pozitronová řešení naopak kladnou energii a elektronová zápornou, takže se na situaci můžeme dívat obráceně a elektron interpretovat jako díru v Diracově moři pozitronů obsazujících záporné energetické stavy. 2.7.6 Elektron a jeho pole, U(1) symetrie Nabité částice v přírodě generují elektromagnetická pole popsaná Maxwellovými rovnicemi (resp. kvantovou teorií elektromagnetického pole) a samy se v těchto polích pohybují ve shodě s Lorentzovou pohybovou rovnicí (resp. Diracovou rovnicí). V této kapitole se nejprve zaměříme na kompletní Lagrangeův popis soustavy pole + elektron a poté se budeme věnovat U(1) symetrii, ze které přímo plyne nutnost existence pole v okolí elektronu. ■ Lagrangeův popis Hustota Lagrangeovy funkce pro interakci nabité částice a elektromagnetického pole má tvar field int part  L L L L (2.180) Polní část známe z teoretické mechaniky, viz vztah (1.194) field 0 1 ; 4 F F F A A             L . (2.181) Také interakční část jsme poznali v teoretické mechanice. Pokud jsou částice popsány čtyřtokem jμ a pole čtyřpotenciálem Aμ je nejjednodušším skalárem kombinace jμ Aμ: int j A L . (2.182) V klasické fyzice je tok náboje částic dán výrazem jμ = (ρQc, j), v kvantové teorii musí tok částic sledovat pravděpodobnost jejich výskytu a tak musí být úměrný toku pravděpodobnosti (2.159). V případě náboje bude koeficientem úměrnosti samotný náboj: j Q    . (2.183) Interakční člen tak získá tvar int j A Q A     L . (2.184) Zbývá nám tedy najít hustotu Lagrangeovy funkce pro samotné částice (elektrony), ze které plyne Diracova rovnice. Velmi jednoduchým skalárem sestaveným přímo za pomoci Diracovy rovnice je výraz Kvantová teorie Relativita a spin 82  part 0i m c     L . (2.185) Pokud budeme interpretovat pole  a ψ jako nezávislá, je hustota Lagrangeovy funkce part part ( , , )   L L a příslušné Lagrangeovy rovnice dají part part 0 part part 0 0 0 ; ( ) 0 0 . ( ) i m c i m c                                                     L L L L První rovnice je Diracova rovnice pro nabitou částici, druhá rovnice je v tuto chvíli jen pomocnou rovnicí pro diracovsky sdružené pole. Povšimněte si, že hmotový člen změnil znaménko. U Feynmanových diagramů to bude odpovídat přítoku hmoty (či odtoku hmoty) do (z) daného vrcholu. U derivací je v obou případech naznačen směr jejich působení. Nyní již můžeme zapsat kompletní hustotu Lagrangeovy funkce pro částici a pole:  0 0 1 ; 4 ; . F F j A i m c F A A j Q                               L (2.186) Poznámky: 1) První člen odpovídá volnému poli, druhý interakci pole s částicí a třetí volné částici. 2) Lagrangeova funkce je funkcí polí , ,A   a jejich derivací. Lagrangeovy rovnice pro pole A dají Maxwellovy rovnice, Lagrangeovy rovnice pro pole  dají Diracovu rovnici. 3) Pokud ponecháme jen první člen, získáme z hustoty Lagrangeovy funkce Maxwellovy rovnice ve vakuu. Pokud ponecháme jen poslední člen, dostaneme Diracovu rovnici volné částice. První a druhý člen dají Maxwellovy rovnice se zdrojovými členy (pole interaguje s částicemi), druhý a třetí člen dají Diracovu rovnici pro částici v přítomnosti elektromagnetického pole (částice interaguje s polem). 4) Interakční člen (druhý) spolu s částicovým členem (třetím) lze sloučit do podoby, která vede na Diracovu rovnici s elektromagnetickým polem:  Dir 0i Q A m c         L (2.187) ■ U(1) symetrie Hustota Lagrangeovy funkce Diracovy částice a čtyřtok reprezentující tok náboje  Dir 0 ;i Q A m c j Q                  L se nezmění při transformaci i i e , e , .A A A                   (2.188) Touto transformací není změněn ani tenzor pole Fμν , a tím ani Lagrangeova funkce Lfield elektromagnetického pole. Celá teorie reprezentovaná lagranžiánem (2.186) je kovariantní vzhledem k transformaci (2.188). Transformace představuje otočení vlnové funkce v každém Kvantová teorie Relativita a spin 83 bodě časoprostoru o stejný úhel α. Jde o unitární operaci (nemění skalární součin) s jedním parametrem (úhlem α), proto se tato transformace označuje U(1). Představuje tzv. vnitřní symetrii teoretického popisu interakce pole-částice. Důsledkem této symetrie je existence elektrického náboje, který se zachovává. U(1) symetrie v jiných teoriích (lagranžiánech) vede na existenci obdobných kvantových nábojů, jako je elektrický, které se v daných procesech zachovávají. ■ U(1)loc symetrie Prozkoumejme nyní, jak by se změnil lagranžián volné částice, pokud bychom připustili, aby úhel α potočení vlnové funkce byl v každém bodě časoprostoru jiný. Představte si nekonečnou prostorovou mříž, v jejímž každém bodě je malý míček. U(1) symetrie v našem modelu odpovídá tomu, že otočíme všechny míčky naráz o stejný úhel. Po této transformaci bude mříž vypadat stejně jako před ní. Nyní si představme, že budeme otáčet v různých časech různé míčky o různé úhly. Výsledek? V libovolném čase naše prostorová mříž vypadá stejně jako před začátkem otáčení. A právě takovou symetrii nazýváme U(1)loc symetrií: i ( ) i ( ) ( ) ( )e , ( ) ( )e . x x x x x x              Proměnná x symbolizuje celou událost x = (t,x).Tato transformace opět nezmění čtyřtok (2.183). Jak se ale změní Lagrangeova funkce volné částice? Dosaďme čárkované veličiny do Lagrangeovy funkce (2.185) a proveďme derivace vlnové funkce a exponenciely:       i ( ) i ( ) part 0 0 , 0 part e e . x x i m c i m c i m c                                         L L Hustota Lagrangeovy funkce po transformaci změnila svůj tvar. Přibyl člen ħγμ α,μ s derivacemi úhlu pootočení, který kopíruje polní člen v hustotě Lagrangeovy funkce nabité částice v přítomnosti pole (2.187). Výsledek je velmi zajímavý. Pokud bychom trvali na tom, aby hustota Lagrangeovy funkce pro volnou částici splňovala symetrii U(1)loc, musíme do teorie přidat elektromagnetické pole. Požadavek, aby byla Diracova rovnice kovariantní vzhledem k U(1)loc symetrii, vede na požadavek existence elektromagnetického pole v okolí částice! Samo elektromagnetické pole se při transformaci změní tak, aby kompenzovalo nově vzniklý člen ħγμ α,μ. Uvažujme tedy transformaci i ( ) i ( ) ( )e , ( )e , ( ) ( ) x x x x A A A x A x                      a opět proveďme výpočet čárkovaného lagranžiánu, tentokrát s elektromagnetickým polem:       part 0 i ( ) i ( ) 0 , 0 e ( ) e ( ) . x x i Q A m c i Q A A m c i Q A Q A m c                                                            L hustota Lagrangeovy funkce se nezmění, pokud bude vnitřní kulatá závorka nulová, tj. ,A Q      . Tím jsme získali návod pro správnou transformaci elektromagnetického pole. Celková transformace U(1)loc bude tedy mít tvar: Kvantová teorie Relativita a spin 84 i ( ) i ( ) ( )e , ( )e , ( ) . x x x x A A A x Q                        (2.189) Jak jsme ukázali, nezmění se při této transformaci součet Lint + LDir. Snadno ověříme, že transformace U(1)loc nemá vliv ani na tenzor elektromagnetického pole a tím na polní část lagranžiánu Lfield. Tok pravděpodobnosti (2.159) transformace také neovlivní. Celá teorie je tak kovariantní vzhledem k U(1)loc transformaci. O kvantové teorii elektromagnetického pole se proto většinou hovoří jako o U(1)loc teorii. Symetrie U(1)loc zajišťuje provázanost nabité částice (elektronu) a elektromagnetického pole. Pokud budeme aplikovat Kleinovu-Gordonovu, resp. Diracovu rovnici na soustavu částic, zjistíme, že statistické chování více částic je u každé z rovnic odlišné. Kleinova-Gordonova rovnice je vhodnou rovnicí pro částice se spinem 0, které nesplňují Pauliho vylučovací princip. Diracova rovnice je naopak vhodná pro částice se spinem ½, které Pauliho vylučovací princip splňují. Chováním soustavy stejných částic v kvantové teorii se budeme zabývat v následující kapitole. Kvantová teorie Stejné částice 85 2.8 SOUSTAVA STEJNÝCH ČÁSTIC Stejnými částicemi nazýváme dvě částice se shodnými parametry (hmotou, nábojem, spinem,…). Z hlediska teoretické mechaniky je trajektorie těchto částic dána Hamiltonovými rovnicemi a známe-li počáteční polohy a rychlosti částic, lze přesně predikovat budoucí polohy částic a v každém okamžiku říci která je která. V kvantové teorii můžeme předpovědět jen pravděpodobnost výskytu částice v nějakém místě a čase. Tato pravděpodobnost má maximum v místě klasické trajektorie a se vzdáleností od ní zpravidla exponenciálně ubývá a dosti daleko od klasické trajektorie je sice velmi malá, nikoli však nulová. Máme-li dvě stejné částice, nikdy si nemůžeme být jisti, která částice je která. Pravděpodobnost výskytu jedné částice v místě druhé je nenulová. Hovoříme o tom, že stejné částice jsou v kvantové teorii nerozlišitelné. Hamiltonův operátor se při záměně dvou stejných částic nezmění: 2112 ˆˆ HH  . (2.190) 2.8.1 Operátor výměny dvou částic Pro jednoduchost budeme uvažovat jen dvě částice, u kterých sledujeme dynamickou proměnnou A (nejlépe celou úplnou množinu pozorovatelných). Stav, ve kterém má první částice hodnotu a1 a druhá částice hodnotu a2 označíme  21,|| aa . Opačnou situaci, kdy první částice má hodnotu a2 a druhá a1, označíme  12,|| aa . Díky nerozlišitelnosti identických částic v kvantové mechanice musí být oba stavy závislé (vyjadřují ve skutečnosti jeden a tentýž kvantový stav), proto  2112 ,|,| aaaa  . (2.191) Zaveďme nyní operátor vzájemné výměny částic vztahem  211212 ,|,|ˆ aaaaP (2.192) a prozkoumejme jeho vlastnosti: .0ˆˆ)3( 1)2( ˆˆ)1( ][ 2,1 2    H,P , ,1P  (2.193) Důkaz (1): Dvojnásobná záměna částic vede na původní konfiguraci. Důkaz (2): Vlastními vektory jsou vektory | a | definované výše:  21122112 ,|,|,|ˆ aaaaaa P . (2.194) Číslo  je vlastním číslem operátoru výměny. Proveďme nyní dvojnásobnou výměnu jednak pomocí prvního vztahu (2.193) a jednak podle (2.194): .11 ,| ,| ,| 2 21 2 21 21 2           aa aa aaP Kvantová teorie Stejné částice 86 Hodnota vlastních čísel operátoru výměny je zřejmá již z prvního vztahu (2.193). Jde o unitární a hermitovský operátor. Vlastní čísla musí ležet na jednotkové kružnici v komplexní rovině a současně být reálná. Jediné takové hodnoty jsou  1. Důkaz (3): V důkazu využijeme časovou Schrödingerovu rovnici (2.98): 2 1 1 2 12 12 1 2 12 2 1 21 2 1 12 1 2 12 12 12 1 2 | , | ,ˆ ˆ ˆ ˆ ˆ| , | , | , i i | ,ˆ ˆ ˆi | , . d a a d a a a a a a a a dt dt d a a a a dt                  H P H H P P P H 2.8.2 Bosony a fermiony, Pauliho princip Z předchozího rozboru je zřejmé, že .,|,|ˆ,| 212112  aaaaaa P (2.195) Vlnová funkce dvou částic může být jen symetrická nebo antisymetrická. Neexistuje nic mezitím. Částice mohou být jen dvojího druhu: se symetrickými vlnovými funkcemi (bosony) nebo s antisymetrickými vlnovými funkcemi (fermiony) vzhledem ke vzájemné záměně. Tuto vlastnost nelze změnit ani časovým vývojem, protože operátor výměny částic podle třetího vztahu (2.193) komutuje s Hamiltonovým operátorem a jeho časový vývoj je proto nulový. Vznikne-li částice jako fermion či boson, zůstává takovou až do svého zániku. ■ Bosony Bosony mají symetrickou vlnovou funkci  2112 ,|,| aaaa . (2.196) Budou li oba stavy stejné, tj. a1 = a2 = a, získáme relaci  aaaa ,|,| , která je vždy splněna a proto může existovat více bosonů ve stejném kvantovém stavu. Při nízkých teplotách mají bosony dokonce snahu kumulovat se v nejnižším možném energetickém stavu a vytvářet tzv. bosonový kondenzát. Ten je známý zejména v supratekutosti a supravodivosti. Statistika, které podléhá soustava bosonů, se nazývá Boseho-Einsteinova statistika a zabýváme se jí v části TF3 (Statistická fyzika). Z dalšího vývoje kvantové mechaniky se ukázalo, že bosony jsou vždy částice s celočíselným spinem (0, 1, 2,…) a pro tyto částice lze zavést kreační operátory splňující jednoduché komutační relace (viz následující kapitola). Nejtypičtějšími představiteli této rodiny jsou skalární (s = 0) a vektorové (s = 1) mezony, dále všechny intermediální částice (foton, W , Z0 a gluony se spinem 0 a graviton se spinem 2). ■ Fermiony Fermiony mají antisymetrickou vlnovou funkci  2112 ,|,| aaaa . (2.197) Budou li oba stavy stejné, tj. a1 = a2 = a, získáme relaci  aaaa ,|,| , která není nikdy splněna, a proto nemůže existovat více fermionů ve stejném kvantovém stavu. Tomuto faktu se říká Pauliho vylučovací princip. Při nízkých teplotách obsazují fermiony postupně jednotlivé energetické hladiny, například v atomárním obalu může být na každé hladině jen tolik elektronů, kolik kvantových stavů tato hladina představuje (to je dáno stupněm degenerace). V atomárním obalu tedy nemohou existovat dva elektrony se stejnými kvantovými čísly n, l, m, ms.. Statistika, které podléhá soustava fermionů, se nazývá FermihoDiracova statistika a budeme se jí zabývat v dalším díle (TF3) tohoto sylabu. Fermiony jsou vždy částice s poločíselným spinem (1/2, 3/2,…) a pro tyto částice lze zavést kreační operátory splňující jednoduché antikomutační relace (viz následující kapitola). Nejtypičtějšími představiteli této rodiny částic jsou leptony (elektron, mion, tauon a neutrina Kvantová teorie Stejné částice 87 se spinem 1/2), kvarky (d, u, s, c, b, t se spinem 1/2), částice složené ze tří kvarků, neboli baryony (neutron, proton,  hyperon se spinem 1/2 a například  baryony se spinem 3/2). BOSONY FERMIONY spin celočíselný poločíselný vlnová funkce symetrická antisymetrická statistika Boseho-Einsteinova Fermiho-Diracova Pauliho princip nesplňují splňují kreační operátory splňují komutační relace splňují antikomutační relace 2.8.3 Druhé kvantování Představme si, že máme N stejných částic, které obsazují stavy nějaké dynamické proměnné. N1 částic je v prvním stavu (hodnota a1), N2 částic je ve druhém stavu (hodnota a2), atd. Čísla Nk nazýváme obsazovací čísla stavu k. Součet všech obsazovacích čísel je roven počtu částic: NN k k  . (2.198) Pro bosony je ,3,2,1,0kN Pro fermiony je situace jednodušší. V daném stavu může být nejvýše jeden fermion, tj. 1,0kN . Příslušný stav soustavy N stejných částic s danými obsazovacími čísly označíme   ,,,,|| 21 kNNN . (2.199) Tomuto zápisu říkáme reprezentace obsazovacích čísel a příslušné stavy nazýváme Fockovy stavy. Dále se situace bude lišit pro bosony a pro fermiony. ■ Bosony Zaveďme podobně jako u harmonického oscilátoru kreační a anihilační operátory do stavu k definičními vztahy (normovací konstanty ponecháme stejné jako u harmonického oscilátoru): .,1,,,|,,,,|ˆ ,,1,,,|1,,,,|ˆ 2121 2121     kkkk kkkk NNNNNNNa NNNNNNNa (2.200) Přímo z těchto definičních relací (pouhým zapůsobením na stavový vektor (2.199) snadno spočteme komutační relace kreačních a anihilačních operátorů: .]ˆ,ˆ[ ,0]ˆ,ˆ[ ,0]ˆ,ˆ[ lklk lk lk aa aa aa      (2.201) Zaveďme další operátor kkk aaN ˆˆˆ   . (2.202) Tento operátor se nazývá (analogicky jako u harmonického oscilátoru) operátor počtu částic ve stavu k, protože zapůsobením na stavový vektor získáme počet částic ve stavu k: Kvantová teorie Stejné částice 88 .,,,,|,,,,| ,1,,,|ˆ,,,,|ˆˆ 2121 2121      kkkkk kkkkkk NNNNNNNNN NNNaNNNNaa Operátor celkového počtu částic potom je    k kk aaN ˆˆˆ . (2.203) Je-li úplná množina pozorovatelných spojitá můžeme celý postup zopakovat pro spojité proměnné. Například v x reprezentaci lze zavést )(ˆ x  kreační operátor do polohy x, )(ˆ x anihilační operátor z polohy x. Komutační relace budou obdobné, jen místo Kroneckerova symbolu vystupuje na pravé straně Diracova  distribuce (viz příloha): .)()](ˆ,)(ˆ[ ,0)](ˆ,)(ˆ[ ,0)](ˆ,)(ˆ[ yxyx yx yx         (2.204) Operátor hustoty počtu částic se zavádí vztahem )(ˆ)(ˆˆ xxk   N , (2.205) operátor počtu částic vyskytujících se v intervalu je    b a dxxxbaN )(ˆ)(ˆ),(ˆ  (2.206) a operátor celkového počtu částic je ˆ ˆ ˆ( ) ( )N x x dx       . (2.207) Obdobně by se postupovalo ve třech dimenzích. Celý přechod od fyziky jedné částice k fyzice mnoha stejných částic lze formálně provézt nahrazením vlnové funkce kreačními a anihilačními operátory a nahrazením hustoty pravděpodobnosti operátorem hustoty počtu částic: )(ˆ)(ˆ)(ˆ)()()( )(ˆ)( * xxxxxxw xx      N . (2.208) Tomuto postupu se říká druhé kvantování, vlnové funkce popisující systém se stávají operátory a kvantová teorie přechází v kvantovou teorii pole, ve které jsou právě veličiny popisující klasická spojitá pole nahrazovány operátory. Druhý řádek přiřazení (2.208) má ještě jeden důležitý význam: U soustavy stejných částic vyjadřujeme pravděpodobnost děje operátorem hustoty počtu částic, tak jak to bývá u skutečných systémů (například svazku stejných částic v experimentu). U jedné částice můžeme hovořit o hustotě pravděpodobnosti jejího výskytu )()(* xx  . Celková pravděpodobnost je rovna jedné, tak, jak to odpovídá normování stavového vektoru. Kvantová teorie Stejné částice 89 ■ Fermiony U fermionů probíhá druhé kvantování obdobně. Opět zavádíme kreační a anihilační operátory lk bb ˆ,ˆ do stavů k a l. Vzhledem k antisymetrii vlnových funkcí musí tyto operátory splňovat antikomutační relace: .0ˆˆˆˆ0,|,|,|,|   kllk bbbbkllkkllk Antikomutátory značíme složenými závorkami a relace (2.201) platná pro bosony, získá pro fermiony tvar: .}ˆ,ˆ{ ,0}ˆ,ˆ{ ,0}ˆ,ˆ{ lklk lk lk bb bb bb      (2.209) Definice spojitých operátorů i operátoru hustoty počtu částic zůstávají shodné. U fermionů jsou všude nahrazeny relace komutační relacemi antikomutačními. V mnoha situacích se chování fermionů a bosonů liší pouze znaménkem (symetrie vlnové funkce; komutační a antikomutační relace; BE a FD statistika). 2.8.4 Ukázka druhého kvantování pro Kleinovo-Gordonovo pole Uvažujme nejjednodušší variantu reálného Kleinova-Gordonova pole pro volnou částici s Lagrangeovou funkcí    2 21 1 2 2        L (2.210) a polní rovnicí  2 0 .   (2.211) Přejdeme-li k soustavě identických částic, změní se pole na operátor ˆ  (2.212) s vlastnostmi † † † ˆ ˆ ˆ ˆ( ), ( ) ( ), ( ) 0 ; ˆ ˆ( ), ( ) ( ) . x y x y x y x y                       (2.213) Veličiny x a y reprezentují celou událost (čas a prostor). Rozviňme nyní polní operátor do rovinných vln (zvlášť označíme kladně frekvenční a zvlášť záporně frekvenční část): † i( ) i( ) 3 ˆ ˆ ˆ( ) ( ) ( )e ( )et t x C a a d            k x k x k k k k . (2.214) V uvedeném vztahu je časová část čtyřvektoru kμ provázána s prostorovou částí přes disperzní relaci ω = ω(k), takže integrace ve skutečnosti probíhá přes všechny čtyři složky. Konstanta C(k) je normovací konstanta, která zajišťuje, aby koeficienty rozvoje (operátory ˆa , ˆa † ) splňovali relace kreačních a anihilačních operátorů. Dosadíme-li rozvoj polního operátoru (2.214) do komutačních relací (2.213), získáme ihned (správná volba C zajistí koeficient 1 u delta funkce v druhé relaci) † † † ˆ ˆ ˆ ˆ( ), ( ) ( ), ( ) 0 ; ˆ ˆ( ), ( ) ( ) . a a a a a a                  k k k k k k k k (2.215) Kvantová teorie Stejné částice 90 Dosadíme-li rozvoj polního operátoru (2.214) do definice Hamiltonovy funkce (1.165) a hybnosti (1.170) máme po elementárních úpravách     † † 3 † † 3 1ˆ ˆ ˆ ˆ ˆ( ) ; 2 1ˆ ˆ ˆ ˆ ˆ( ) , 2 H E aa a a d aa a a d       k k P p k k (2.216) kde jsme označili 2 2 2 2 2 2 2 4 0 ( ) ; ( ) ( ) .E k k c k c p c m c        p k k   (2.217) Vztah ω(k) je dán disperzní relací (2.122). S využitím komutačních relací (2.215) získáme přímo dosazením relace † † † † ˆ ˆ ˆ, ( ) ( ) ( ) ; ˆ ˆ ˆ, ( ) ( ) ( ) ; ˆ ˆ ˆ, ( ) ( ) ( ) ; ˆ ˆ ˆ, ( ) ( ) ( ) . H a E k a H a E k a a a a a                       k k k k P k p k k P k p k k (2.218) Z těchto relací je zjevný význam operátorů ˆa , ˆa + : Pole φ je kvantováno a operátor ˆa † kreuje kvantum pole s energií E(k) a hybností p(k), operátor ˆa stejné kvantum anihiluje. Toto kvantum můžeme interpretovat jako boson s nulovým spinem (pole má jedinou vlnovou funkci odpovídající jediné projekci spinu). Kleinova-Gordonova rovnice získává po druhém kvantování názornou interpretaci. Jde o pole, které můžeme chápat jako soustavu excitací – bosonů s nulovým spinem. ■ Komplexní pole Pokud by Kleinovo–Gordonovo pole bylo komplexní, tj. † 1 2 1 2ˆ ˆ ˆ ˆ ˆ ˆi ; i ,         lze ukázat, že excitace takového pole odpovídají dvěma druhům skalárních bosonů (se spinem 0), které jsou sobě navzájem antičásticemi. ■ Normální uspořádání operátorů Ve vztazích pro energii a hybnost (2.216) se skrývá jeden problém. Pokud bychom hledali střední hodnotu energie a hybnosti ve vakuovém stavu, dostaneme nekonečné hodnoty. Za to může první člen, ve kterém je kreační operátor napravo a na vakuový stav dá nenulovou hodnotu:  † † 3 † 3 3 3 1ˆ ˆ ˆ ˆ ˆ0 0 ( ) 0 0 0 0 2 1 ˆ ˆ( ) 0 0 ( ) 0 0 ( ) . 2 H E aa a a d E aa d E d E d           k k k k k k k k  Tento problém je důsledkem principu korespondence, který neřeší správné pořadí operátorů, které jsou v součinu. Máme-li dvě dynamické proměnné A a B, můžeme součinu AB v kvantové teorii přiřadit dvě možná pořadí operátorů: ˆ ˆ , ˆˆ . AB AB BA     (2.219) Kvantová teorie Stejné částice 91 Jen jedno pořadí bude ovšem odpovídat dějům v přírodě. Toto pořadí není dáno principem korespondence, ale musíme ho vybrat tak, aby výsledky byly v souladu s pozorováním. Správné pořadí operátorů se nazývá normální uspořádání a označujeme ho dvojtečkou, tedy ˆ ˆ: :AB Tento zápis znamená, že pořadí operátorů mezi dvojtečkami není určeno jednoznačně a musíme ho volit ve shodě s experimentem. Použít lze například následující postup: 1) operátory vyjádříme za pomoci příslušných kreačních a anihilačních operátorů (bosonové splňují komutační relace a fermionové antikomutační relace), 2) v součinech budeme anihilační operátory přesouvat doprava podle následujících pravidel: a. pokud vedle sebe jsou dva bosonové operátory, přesuneme anihilační doprava, b. pokud vedle sebe je jeden bosonový a jeden fermionový operátor, přesuneme anihilační doprava, c. pokud vedle sebe jsou dva fermionové operátory, přesuneme anihilační doprava a zaměníme znaménko daného členu. Příklad 22: Nalezněte správné pořadí operátorů v hamiltoniánu Kleinova-Gordonova pole:       † † 3 † † 3 † † 3 † 3 3 1 1ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ: ( ) : ( ) : : 2 2 1 ˆˆ ˆ ˆ ˆ ˆ ˆ( ) ( ) ( ) ( ) 2 H E aa a a d E aa a a d E a a a a d E a a d E d               k k k k k k k k k k kN Střední hodnota hamiltoniánu ve vakuovém stavu již nediverguje, navíc je struktura Hamiltonova operátoru zcela zřejmá, ˆ( )kN je operátor hustoty počtu částic s vlnovým vektorem k. Obdobně musíme upravit i vztah (2.216) pro hybnost. Správné relace tedy jsou: † 3 † 3 ˆ ˆ ˆ( ) ; ˆ ˆ ˆ( ) . H E a a d a a d     k k P p k k (2.220) Příklad 23: Určete správné pořadí operátorů ve výrazu (a jsou bosonové operátory, d fermionové)  † † † † † †1 ˆ ˆ ˆ ˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ: : 2 A aa a a dd d d ad d a      Aplikací výše uvedených pravidel získáme snadno výsledek  † † † † † † † †1 ˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆˆ ˆ ˆˆ ˆ ˆ ˆ ˆ ˆ ˆ ˆ 2 a dA a a a a d d d d d a d a a a d d N N          . Kvantová teorie Zobecněné funkce 92 PŘÍLOHA – ZOBECNĚNÉ FUNKCE P1. Diracova distribuce Ve fyzice se velmi často setkáváme s nutností popsat bodový náboj nebo hmotný bod. Náboj či hmotnost částice si představujeme lokalizované v jediném místě, což s sebou nese problém nekonečné hustoty náboje či hmoty v tomto místě. Řešením je zavedení tzv. zobecněných funkcí, zejména Diracovy distribuce. Ukažme si problém na lineární hustotě náboje lokalizovaného v místě x = 0: 0 ; 0 ( ) 0 0. x x x       (2.221) Integrál z hustoty ale musí dát celkový náboj Q: ( ) .x dx Q    (2.222) Je jasné, že hustota náboje není „normální“ funkcí. Má nenulovou hodnotu v jediném bodě a integrál z ní by přesto měl dát konečné číslo. Takové funkce ale neexistují, můžeme je zavádět jako limitu posloupností funkcí a jejich význam je jen ve skalárním součinu s jinou, tzv. testovací funkcí. f ( )x x  1/ f ( )x x f ( )x x Posloupnost obdélníků Zaveďme si obdélníkové funkce 1/ , /2, /2 ; ( ) 0, /2, /2 . x f x x               (2.223) Všechny obdélníky mají stejnou plochu rovnou jedné a funkce mají zajímavé vlastnosti: 0 pro 01 ( ) 1; (0) ; lim ( ) . 0 pro 0 x f x dx f f x x                (2.224) Diracovu distribuci můžeme formálně zavést jako limitu těchto obdélníkových funkcí 0 ( ) lim ( ) .x f x     (2.225) Posloupnost kopečků Obdélníky z předchozí ukázky nejsou hladké funkce. To ale není nepřekonatelný problém, místo obdélníků můžeme použít funkce spojité se všemi svými derivacemi podle vztahu 2 2 1 ( ) .f x x       (2.226) Plocha pod těmito funkcemi je rovna jedné pro každé , protože Kvantová teorie Zobecněné funkce 93 2 2 1 1 1 ( ) arctan 1. 2 2 x f x x                               (2.227) Pro malá  se „kopce“ zužují a přitom se zvětšuje jejich výška: 0 pro 01 ( ) 1; (0) ; lim ( ) . 0 pro 0 x f x dx f f x x                (2.228) Opět můžeme zavést Diracovu distribuci jako limitu těchto spojitých funkcí: 0 ( ) lim ( ) .x f x     (2.229) Posloupnost Dirichletových jader Diracovu Distribuci můžeme zavést také pomocí jednoduché funkce sin ( ) ; (0) 1; ( ) . x f x f f x dx x       Zaveďme posloupnost sin ( ) ,k k kx f x kx  (2.230) která má jednoduché vlastnosti pro 0 ( ) 1; (0) ; lim ( ) . 0 pro 0 k k k k xk f x dx f f x x            Diracovu distribuci lze zavést jako limitu funkcí ( ) lim ( ) .k k x f x   Poznamenejme, že funkce ( )kf x jsou známé z důkazu věty o Fourierově rozvoji do řady a nazývají se Dirichletovo jádro. Fourierův obraz jednotkové funkce Spočtěme nejprve následující integrál: 1 e e sin e e 2 . kk ikx ikx ikx ikx kk kx dk k ix ix kx              Integrál až na koeficient /2 dává Dirichletovo jádro. Diracovu distribuci lze proto napsat jako 1 1 ( ) lim e e 2 2 k ikx ikx k k x dk dk           . (2.231) Integrál v nevlastních mezích chápeme právě ve smyslu uvedené limity. Diracova distribuce je tak Fourierovým obrazem jednotkové funkce. Diracova distribuce nemá vlastnosti běžných funkcí. Přestože je její hodnota nenulová v jediném bodě, dá integrál z ní nenulovou hodnotu. To plyne z limitního charakteru zavedení této distribuce. K jejím základním vlastnostem patří: ( ) ( ) ( ) (0) (0) ( ) (0).x f x dx x f dx f x dx f              (2.232) Kvantová teorie Zobecněné funkce 94 Důvod je snad zřejmý. Distribuce  je všude nulová kromě jediného bodu x = 0. Proto výsledek integrálu může ovlivnit jedině hodnota funkce f v počátku. Tu však můžeme vytknout před integrál a dostaneme jako výsledek hodnotu funkce v počátku. Poznámka 1: Distribuci lze také chápat jako velmi jednoduché zobrazení, které přiřadí funkci její hodnotu v počátku (zobrazení, které přiřadí funkci číslo se nazývá funkcionál). ˆ ( ) (0) ; resp. ( ) (0) T T f x f f x f     . Poznámka 2: Distribuci lze chápat jako funkcionál daný skalárním součinem ˆ ( ) ;gT f x g f Skalární součin působí na libovolnou funkci f z tzv. prostoru testovacích funkcí. Funkce g je pevně daná, definuje toto zobrazení a nazývá se temperovaná distribuce. Čím hezčí vlastnosti budou mít funkce z testovacího prostoru (například budou dostatečně rychle konvergovat k nule na hranicích oblasti), tím horší vlastnosti může mít funkce g definující zobrazení. Za prostor testovacích funkcí může posloužit například Schwarřtzův (Sobolevův) prostor. Poznámka 3: Často se hledají řešení celých rovnic „ve smyslu skalárního součinu“. Například místo rovnice f  řešíme rovnici 0f    , kde  je hledané řešení a  je libovolná funkce z prostoru testovacích funkcí. Tato řešení se nazývají slabá řešení. Jejich třída je mnohem bohatší než byla třída řešení původní rovnice. Nacházená řešení mohou mít „divočejší“ charakter a jsou bližší fyzikální realitě. Jejich hledáním se zabývala vynikající matematička Ladyženská. P2. Konvoluce Na separabilních prostorech (se spočetnou bází) můžeme zobrazení ˆA f g psát v konkrétní reprezentaci v maticovém tvaru .kl l k l A f g (2.233) Jednotkové zobrazení ˆ f f1 je dáno jednotkovou maticí, jejíž prvky tvoří Kroneckerův symbol: .kl l k l f f  (2.234) V případě neseparabilních prostorů je zobrazení dáno funkcí dvou proměnných ( , ) ( ) ( )A x y f y dy g x   . (2.235) Integrál (2.235) se nazývá konvoluce a označuje se * ( , ) ( )A f A x y f y dy    . (2.236) Konvoluce je analogií maticového násobení na neseparabilních prostorech. Roli indexů přebírají spojité proměnné x a y. Roli matice přebírá tzv. jádro konvoluce ( , )A x y . Speciálním případem konvolucí jsou různé integrální transformace (Laplaceova, Fourierova, Abelova, atd.). Jádrem jednotkového operátoru je Diracova distribuce (je nenulová jen pro x y ): ( ) ( ) ( )x y f y dy f x   Diracova distribuce tak na neseparabilních prostorech přebírá úlohu Kroneckerova symbolu. Kvantová teorie Zobecněné funkce 95 P3. Greenův operátor a Greenova funkce Napišme maticové elementy jednotkového operátoru v x reprezentaci (maticové elementy jednotkového operátoru jsou právě Diracovou distribucí): *ˆ( ) 1 ( ) ( ) .n n n n x y y x y n n x f y f x      Ve spojitých prostorech * ( ) ( ) ( )k k k x y f y f x dk    (2.237) Distribuci lze tak napsat pomocí libovolných bázových funkcí, například pomocí báze 1 e 2 ikx k   , dostaneme ( )1 1 1 ( ) e e e ( ) e 2 2 2 iky ikx ik x y ikx k k k x y dk dk x dk              , (2.238) což je výše odvozený vztah (2.231). Greenův operátor Hledejme řešení lineární operátorové rovnice s pravou stranou ˆL f  . (2.239) Z věty o spektrálním rozvoji víme, že řešení je možné zapsat pomocí vlastních vektorů (tvoříli ortonormální bázi) a vlastních ve tvaru 1 ll l l f    . Přepišme řešení takto 1ˆ ˆ; ll G f G l l     . (2.240) Operátor ˆG se nazývá Greenův operátor a je inverzním operátorem k operátoru ˆL . V případě operátoru se spojitým spektrem přejde sumace v integraci. Greenova funkce Zabývejme se nyní speciálním případem - rovnicí s lineárním operátorem a nenulovou pravou stranou na prostoru L2 ˆL f  . (2.241) Hledejme nejprve řešení pro jednotkový impuls na pravé straně (bude reprezentovaný Diracovou distribucí): ˆ ( ) ( )LG x x Toto řešení se nazývá Greenova funkce. Obecné řešení rovnice (2.241) je konvolucí Greenovy funkce a pravé strany rovnice ( ) * ( ) ( )x G f G x y f y dy    . Důkaz je velmi jednoduchý. Ukážeme, že působením operátoru ˆL na nalezené řešení dostaneme pravou stranu původní rovnice: ˆ ˆ( ) ( ) ( ) ( ) ( ) ( ) .L x LG x y f y dy x y f y dy f x      