Katedra teoretické fyziky Přírodovědecká fakulta Univerzita Palackého VYBRANÉ PARTIE SOUČASNÉ FYZIKY Tomáš Opatrný a Lukáš Richterek Olomouc 2005 Abstrakt Tento text si klade za cíl seznámit čtenáře s vybranými partiemi současné fyziky. První část pojednává o základních představách moderní kosmologie, o nejdůležitěšjích vlastnostech Vesmíru a o možných scénářích jeho dalšího vývoje. Seznámení se základnímí principů by mělo čtenáře připravit a motivovat ke studiu další odborné literatury. Druhá oblast se týká kvantového zpracování informace. Dozvíme se tam o tom, proč nás velmi zajímá možnost sestrojení kvantového počítače, jak by takový počítač pracoval a k čemu by mohl být užitečný, co to je kvantová kryptografie a teleportace kvantových stavů. Náš text v žádném případě nepokrývá zvolená témata vyčerpávájícím způsobem, měl by spíše sloužit k základní orientaci a motivovat k hlubšímu studiu problémů, jimž se ostatní fyzikální kurzy na PřF UP věnují jen okrajově nebo vůbec. Od čtenáře – studenta předpokládáme základní znalost diferenciálního i integrálního počtu a fyzikálních pojmů na úrovni základního kurzu fyziky na PřF UP a též základní znalosti z kvantové teorie. Na několika místech se odvoláváme teorii relativity a na statistickou fyziku; znalost těchto disciplin popř. jejich matematického aparátu je pro čtenáře bezesporu výhodou. Cílová skupina Text je primárně určen jako doplněk k volitelnému předmětu „Vybrané kapitoly z fyziky pro 5. ročník učitelství všeobecně vzdělávacích předmětů kombinací, pro něž je fyzika jedním z aprobačních předmětů. Doufáme, že by mohl oslovit i ostatní zájemce o fyziku ať už z řad studentů UP či veřejnosti. Obsah 1 Friedmannovy kosmologické modely . . . . . . . . . . . . . . . . 5 1.1 Základní východiska a principy relativistické kosmologie . . . . . . . . . . . . . . 6 1.2 Friedmannova rovnice . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12 1.3 Robertsonova-Walkerova metrika . . . . . . . . . . . . . . . . . . . . . . . . . . . 17 1.4 Hustota a tlak . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 1.5 Hubbleův zákon . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22 1.6 Šíření světla ve Friedmannových modelech . . . . . . . . . . . . . . . . . . . . . . 25 1.7 Observační parametry vesmíru . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29 1.8 Budoucnost vesmíru . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32 1.9 Stáří pozorovaných objektů . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35 1.10 Fotometrická vzdálenost a Hubbleovy diagramy . . . . . . . . . . . . . . . . . . . 38 1.11 Závěr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44 Literatura ke kapitole 1 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45 2 Základy kvantové informace . . . . . . . . . . . . . . . . . . . . . 47 2.1 Úvod a něco z historie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47 2.1.1 Moorův zákon . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47 2.1.2 Limity výpočetních možností . . . . . . . . . . . . . . . . . . . . . . . . . 48 2.1.3 Ať počítá kvantový systém! . . . . . . . . . . . . . . . . . . . . . . . . . . 49 2.1.4 Problémy konstrukce kvantového počítače . . . . . . . . . . . . . . . . . . 52 2.1.5 Jiné aplikace kvantové informatiky . . . . . . . . . . . . . . . . . . . . . . 53 2.2 Kvantové bity neboli qubity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56 2.2.1 Co je to qubit . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56 2.2.2 Měření kvantových bitů . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58 2.3 Kvantová hradla . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60 3 2.3.1 Reverzibilita kvantového počítání . . . . . . . . . . . . . . . . . . . . . . . 60 2.3.2 Jednobitová hradla . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62 2.3.3 Dvoubitová hradla . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64 2.3.4 Kombinace kvantových hradel . . . . . . . . . . . . . . . . . . . . . . . . . 66 2.3.5 Kvantová Fourierova transformace . . . . . . . . . . . . . . . . . . . . . . . 68 2.4 Shorův algoritmus: nalezení periody diskrétní funkce . . . . . . . . . . . . . . . . 70 2.5 Faktorizace čísel . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73 2.6 Groverův algoritmus: vyhledávání v neuspořádaných seznamech . . . . . . . . . . 75 2.6.1 Popis algoritmu . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75 2.6.2 Překlápění kolem střední hodnoty amplitud . . . . . . . . . . . . . . . . . 76 2.6.3 Překlopení znaménka amplitudy . . . . . . . . . . . . . . . . . . . . . . . . 77 2.7 Kvantová kryptografie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 78 2.7.1 Distribuce klíče . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 79 2.7.2 Odhalení narušitele . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80 2.7.3 Rozvoj kvantové kryptografie . . . . . . . . . . . . . . . . . . . . . . . . . 81 2.8 Kvantová provázanost . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83 2.9 Kvantová teleportace . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86 2.10 Závěr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88 Literatura ke kapitole 2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88 Přílohy A Fyzikální konstanty a jednotky použité v textu . . . . . . . . . 90 B RSA kryptografický kód . . . . . . . . . . . . . . . . . . . . . . . 91 Kapitola 1 Friedmannovy kosmologické modely Studijní cíle: V této kapitole odvodíme základní rovnice popisující rozpínání (i případné smršťování) vesmíru. Ukážeme si vliv základních kosmologických parametrů na možný budoucí osud vesmíru, budeme diskutovat jejich měření i nejpravděpodobnější hodnoty. Klíčová slova: Friedmannova rovnice, Hubbleův a decelerační parametr, hustotní parametr, kosmologická konstanta Λ, Hubbleovy diagramy, inflace. Potřebný čas: 300 minut. Otázky spojené s existencí, vlastnostmi našeho vesmíru i s postavením člověka v něm se vždy promítaly do základů lidské kultury, náboženských a filozofických směrů. Od dob Galileových a Newtonových se rozvíjí kosmologické teorie vycházející z fundamentálních fyzikálních zákonů, jež kladou důraz jak na přesné výpočty a předpovědi, tak především na konfrontaci se stále se zpřesňujícími pozorováními. V tomto smyslu je fyzikální kosmologie vědeckou, tj. vyvratitelnou teorií, v níž jsou naše představy a modely neustále korigovány, aby bylo dosaženo stále větší a lepší shody s tím, co v našem vesmíru pozorujeme. Průvodce studiem Naše chápaní vesmíru a jeho vývoje s neustále vyvíjí, jak po stránce teoretické, tak získáváním a zpracováním stále přesnějších experimentálních dat. V moderní kosmologii se doslova a do písmene uplatňuje celá fyzika – od fyziky mikrosvěta a kvantové teorie pole až po teorii relativity. Zde se budeme zabývat pouze základy relativistické kosmologie, základními pozorovanými vlastnostmi vesmíru a z nich vyplývajícími nejjednoduššími modely. Otázky týkající se raného a velmi raného vesmíru v prvních minutách po velkém třesku, inflační teorie a problematiky formování struktur najde čtenář v podrobnějších monografiích, např. [1.17, 1.23, 1.27, 1.30, 1.32]. 5 Obr. 1.1: Galaxie v souhvězdí Andromedy známá pod označením z katalogu francouzského astronoma Charlese Messiera (1730-1817) jako M31, v NGC katalogu J.L.E. Dreyera (1852–1926) je uvedena pod číslem NGC 224 (oba katalogy jsou dnes k dispozici na internetových adresách http://www.seds.org/messier/ resp. http://www.seds.org/ ~spider/ngc/ngc.html). Jedná se o superobří spirální galaxii typu Sb ve vzdálenosti více jak 2·106 ly, spolu s naší Galaxií je největší v naší místní skupině. Na obrázku vidíme také dvě satelitní trpasličí eliptické galaxie NGC 205 a NGC 221 (M32). Galaxie M31 představuje zvětšenou obdobu naší vlastní Galaxie s hmotností okolo 3·109 M . Do místní skupiny Galaxií patří celkem asi 30 galaxií, z nichž nejznámější jsou satelity naší Mléčné dráhy – nepravidelné galaxie Velké a Malé Magellanovo mračno (LMC neboli Large Magellanic Cloud a SMC neboli Small Magellanic Cloud). Galaxie v místní skupině tvoří gravitačně vázaný systém a díky tomu se od sebe nevzdalují v důsledku rozpínání vesmíru (viz část 1.5), naopak lze předpokládat, že ve velmi vzdálené budoucnosti může dojít ke srážkám např. naší Galaxie a M31. 1.1 Základní východiska a principy relativistické kosmologie V tomto textu se budeme zabývat společnými základy moderních kosmologických teorií, které můžeme shrnout označením standardní kosmologický model. Vychází ze dvou základních principů, zobecněného Koperníkova principu a principu uniformity. Jádro prvního tvoří tvrzení, že se nenalézáme ve středu vesmíru ani v jeho význačMikoláš Koperník (1473–1543) ném bodě. V historickém kontextu jde o velmi důležitý myšlenkový zlom, neboť ve starověku i středověku byla právě Země považována za střed všehomíra, okolo něhož obíhala všechna ostatní tělesa. Průlom představoval heliocentrický model sluneční soustavy Mikoláše Koperníka z r. 1543, jenž zbavil Zemi „výsadního postavení. 6 Dnes víme, že z kosmologického hlediska není význačná ani naše Galaxie nebo kupa galaxií, dokonce i hmota, z níž jsou složena naše těla představuje pouze malou část v porovnání s převažujícím typem hmoty ve vesmíru. Homogenita a izotropie na velkých vzdálenostech představují klíčové vlastnosti vesmíru! „Princip uniformity nebo také „princip homogenity a izotropie tvrdí, že vesmír je ve velkých měřítkách homogenní a izotropní, tj. stejný ve všech místech a ve všech směrech. Je zřejmé, že např. v měřítkách naší sluneční soustavy, popř. naší Galaxie uvedený princip neplatí; ve Slunci popř. v jádru Galaxie je soustředěno mnohem více hmoty než na okraji. Homogenita a izotropie se projevuje až na velkých škálách o rozměrech řádově 103 Mpc, jak o tom svědčí výsledky řady pozorování (viz obr. 1.2– 1.4) O vlastnostech vesmíru vypovídají i zdánlivě triviální otázky; příkladem může být tzv. Olbersův paradox. Otázka „Jak to, že je noční obloha temná? byla poprvé zformulována Keplerem v roce 1610, diskutována Halleym a Cheseauxem v 18. století a zpopularizována německým lékařem a astronomem Heinrichem Wilhelmem Matthäusem Olbersem (objevitelem planetek Pallas a Vesta i komety 13P/Olbers) v r. 1826. Význam otázky vynikne k kontextu dobových filozofických představ, užitečné je i srovnání s odpovědí, které dává současná kosmologie. Pokud by vesmír byl nekonečný a obsahoval nekonečné množství rovnoměrně rozmístěných hvězd, pak by noční obloha měla být stejně jasná jako povrch Slunce. Zjednodušeně řečeno, v každém místě oblohy bychom nalezli svítící hvězdu.1 H. W. M. Olbers (1758–1840)Možná řešení paradoxu jsou následující: Olbersův paradox: i jednoduchá otázka vede k hlubšímu zamyšlení! 1. mezihvězdný prach nám brání vidět vzdálené hvězdy; 2. vesmír obsahuje pouze konečný počet hvězd; 3. hvězdy nejsou rozmístěny rovnoměrně 4. vesmír se rozpíná, světlo od nejvzdálenějších hvězd má takový rudý posuv, že je mimo oblast viditelného světla; 5. vesmír má konečné stáří, světlo od nejvzdálenějších objektů k nám ještě nedora- zilo. Z hlediska moderní kosmologie se nejvýznamněji se uplatňuje právě poslední možnost – vesmír existuje pouze konečnou dobu, podle nejnovějších odhadů asi (13,7 ± 0,2) miliard světelných let [1.5]. Shrnutí Vesmír je ve velkých měřítkách homogenní a izotropní; svědčí o tom rozbory kosmického mikrovlnného záření a statistické sledování rozmístění galaxií v něm. Nejen z Olbersova paradoxu vyplývá, že vesmír netrvá „odjakživa , ale pouze konečnou dobu. 1 Přesnější formulace vychází s poznatku, že světelný tok od vzdálených zdrojů klesá se čtvercem jejich vzdálenosti od nás r, tj. s 1/r2 , ale objem kulové vrstvy, v níž se hvězdy nacházejí je roven 4pr2 dr, výsledný světelný tok z takové vrstvy tak nezávisí na její vzdálenosti; podrobnější diskusi lze najít např. v [1.40]. 7 T = 2, 728 K ∆T = 3, 353 mK ∆T = 18 µK Obr. 1.2: Teplota kosmického mikrovlnného záření podle výsledků družice COBE, konkrétně její součásti DMR (Differential Microwave Radiometer). Tři mapy teploty záření pro různé rozsahy jsou vyneseny v galaktických souřadnicích, rovina naší Galaxie se nachází uprostřed. Horní mapa dokládá, že záření je vysoce homogenní s teplotou T ≈ 2,73 K. Prostřední část ukazuje rozdíly v teplotě řádu mK. Zřetelná dipólová anizotropie odpovídá pohybu sluneční soustavy vůči kosmickému mikrovlnnému záření (jeho klidové soustavě). Pokud tuto anizotropii odečteme, získáme zbývající anizotropii řádu µK na dolní mapě, přičemž výrazný červený pruh uprostřed odpovídá záření z naší Galaxie; více obrázků a informací lze nalézt na http://lambda.gsfc.nasa.gov/product/cobe/ a v [1.4]. Studium anizotropií mikrovlnného záření, jež jsou podle současných představ „otiskem fluktuací hustoty vesmíru v době vzniku mikrovlnného záření (tj. asi 300- 400 000 let po velkém třesku), je považováno za klíč k pochopení vzniku pozorovaných struktur (galaxií). Poskytuje dosud nejpřesnější určení stáří vesmíru i dalších kosmologických parametrů – viz stránky projektu WMAP (Wilkinson Microwave Anisotropy Probe) http://map.gsfc.nasa.gov/ nebo např. [1.37] 8 Obr. 1.3: Úhlové rozdělení galaxií podle projektu APM (Automated Plate Measuring, domovská stránka http://www-astro.physics.ox.ac.uk/~wjs/apm_survey.html) [1.24]. Zahrnuta je asi 1/4 oblohy rozdělená do menších čtverečku, jejichž barva odráží počet galaxií v každém z nich (čím světlejší, tím více). Vidíme, že na velkých škálách se počet galaxií v jednotlivých směrech neliší. Pojmy k zapamatování • Koperníkův princip • princip homogenity • princip uniformity • Olbersův paradox Kontrolní otázky 1. Co rozumíme pod Koperníkovým principem? 2. Jaké nejstarší objekty v současné době ve vesmíru pozorujeme Hubbleovým te- leskopem? 3. Jaké je přibližně stáří vesmíru? Cvičení 1. Naše Galaxie s hmotností 2,5·1011 M a galaxie M 31 v souhvězdí Andromedy o hmotnosti 3,6·1011 M jsou dvě největší galaxie v tzv. místní soustavě galaxií. Předpokládejme, že tvoří dvojnou soustavu a obíhají kolem společného hmotného středu po kruhových drahách. Určete velikost oběžné doby mezi nimi, jestliže vzdálenost mezi nimi je asi 700 kpc [1.38]. 9 rud´y posuv miliardy svˇeteln´ych let Obr. 1.4: 2dF Galaxy redshift survey ukazuje polohu 62 559 z celkového počtu 220 929 galaxií, u nichž byly určeny úhlové souřadnice a rudý posuv [1.9]. V souladu s Hubbleovým zákonem (1.5.1) je radiální vzdálenost galaxií určena rudým posuvem z. Zdánlivě galaxií v prostoru s rostoucí vzdáleností od nás ubývá, ale jde jen o výběrový efekt – ve velkých vzdálenostech můžeme pozorovat pouze nejzářivější objekty. Pro malé rudé posuvy, kde je počet galaxií zmapován důkladněji, vidíme zřetelně vláknovitou strukturu s uzly a mezerami. Homogenita vesmíru předpokládá, že rozmístění galaxií ve vzdálenějších oblastech musí být stejné. Více obrázků a informací lze nalézt na stránkách projektu http://www.mso.anu.edu.au/2dFGRS/. Řešení 1. Podle 3. Keplerova zákona platí a3 T2 = G 4p2 (M1 + M2) , odkud T = 4p2a3 G (M1 + M2) . Po dosazení v jednotkách SI G = 6,67·10−11 m3 ·kg·−1 s−2 , a = 700 kpc ≈ 2,16·1022 m, M1 ≈ 4,975·1041 kg a M1 ≈ 7,161·1041 kg vychází T ≈ 2,22·1018 s ≈ 7,03·1010 let. Galaxie oběhnou okolo společného hmotného středu asi za 70 miliard let. Průvodce studiem K plnému porozumění moderní kosmologii a korektnímu odvození rovnic, popisujících expanzi vesmíru je nezbytná znalost obecné teorie relativity. Avšak 10 Obr. 1.5: Pohled na velmi vzdálené galaxie Hubbleovým teleskopem označovaný jako „the Hubble deep field . Snímek pokrývá malou část oblohy asi 30× menší než Měsíc, tak malou, že se v ní nachází pouze několik jasných hvězd z naší Galaxie. Naopak vidíme řadu galaxií v různých stadiích vývoje – od spirálních a eliptických, až po slabě zářící nepravidelné objekty, jejichž světlo bylo vysláno méně jak miliardu let po velkém třesku, tj. zhruba před 12–13 miliardami let, a které patří k nejstarším dosud pozorovaným objektům ve vesmíru. Vzhledem k homogenitě a izotropii lze i tento malý vzorek oblohy považovat za reprezentativní, typické rozmístění galaxií. i čtenář, jenž se studiem Einsteinovy teorie dosud nezabýval, může řadu jevů odvodit a popsat na základě elementárnějších úvah vyžadujících znalosti na úrovni základního vysokoškolského kurzu fyziky. Potřebuje k tomu především 11 Friedmannovu rovnici (1.2.3) a rovnici pro práci vykonanou tlakem při rozpínání vesmíru. Zájemce o nástin plně relativistického odvození rovnic odkazujeme na následující podkapitolu 1.3. 1.2 Friedmannova rovnice V této části vyjdeme z modelu Newtonovského vesmíru, v němž uvažujeme pouze nerelativistické pohyby a gravitace je chápána jako síla působící mezi hmotnými částicemi [1.19, 1.23]. Friedmannova rovnice (1.2.3) v tomto případě popisuje zachování celkové mechanické energie částice (galaxie) v průběhu vesmírné expanze. V relativistickém odvození namísto klasické hustoty hmotnosti musíme započítat hustotu energie v souladu se známou rovnicí ekvivalence hmotnosti a energie E = mc2 a namísto gravitace jako síly uvažovat zakřivení prostoročasu, nicméně výsledek je formálně tentýž, bez jakýchkoli korekčních členů či přiblížení. Do hustoty energie však musíme započítat všechny její formy. K pochopení a správné interpretaci šíření světla v zakřiveném vesmíru (Hubbleovy diagramy, fotometrická vzdálenost apod. – viz část 1.10) se pak bez teorie relativity neobejdeme! Rozpínání vesmíru pozorujeme jako vzájemné vzdalování galaxií objevené Edwinem Hubblem (viz část 1.5). Náš newtonovský model si lze představit jako rozpína- Východiskem našich úvah je zákon zachování energie. jící se plyn, jehož „molekuly jsou celé galaxie (ty samy se nerozpínají). Gravitace naopak galaxie přitahuje k sobě a rozpínání zpomaluje. Zvolíme-li počátek souřadnic v „místě některé galaxie (např. naší), potom na galaxii ve vzdálenosti R působí pouze hmota obsažená v kouli o poloměru R (působení vnějších vrstev se vyruší, to je důležitá vlastnost gravitační i coulombovské síly v klasické fyzice), jejíž hmotnost M = 4pR3 /3, kde je hustota vesmírné „tekutiny . Mechanická energie uvažované galaxie bude potom E = 1 2 m dR dt 2 − GMm R = 1 2 m dR dt 2 − 4p 3 m R2 . Předpokládáme-li, že na galaxii působí ještě izotropní síla F = Λmc2 R/3, jíž odpovídá potenciální energie2 Podle nejnovějších pozorování musíme počítat s kosmologickou konstantou. U = − 1 6 Λmc2 R2 , vychází celková energie dané galaxie E = 1 2 m dR dt 2 − 1 6 Λmc2 R2 − 4p 3 m R2 . (1.2.1) Parametr Λ o rozměru m−2 nazýváme kosmologickou konstantou a standardně bývá interpretována jako vliv energie samotného vakua. 2 V daném kontextu jde o ad hoc předpoklad; zavedení kosmologické konstanty přirozeně plyne z Einsteinových rovnic obecné teorie relativity (1.3.2) a v důsledku novějších pozorování (viz např. [1.5, 1.20]) o nutnosti zahrnout i tento člen dnes málokdo pochybuje. 12 ˇcas Obr. 1.6: K zavedení „comoving souřadnic: souřadnicový systém je unášen expanzí, takže uvažovaný objekt zůstává stále ve zvoleném bodě souřadnicové sítě, ostatní body se od něj v důsledku rozpínání vzdalují. Vzhledem k homogenitě a izotropii rozpínání je výhodné zavést tzv. „comoving souřadnice, které jsou „unášeny rozpínáním a pro zvolenou galaxii zůstávají po celou dobu konstantní. Zvětšování vzájemných vzdáleností je pak popsáno expanzním faktorem a = a(t) závisejícím na čase (viz obr. 1.6). Pro skutečnou vzdálenost proto platí Pro jednodušší popis zavádíme „comoving souřadnice a expanzní faktor. R = a(t)x (1.2.2) a po dosazení do (1.2.1) získáváme 2E mx2 = da dt 2 − 1 3 Λc2 a2 − 8p 3 a2 . Jak energie E, hmotnost m i „comoving souřadnice x jsou v rámci našeho newtonovského modelu pro uvažovanou galaxii konstantami, použijeme-li ekvivalenci hmotnosti a energie E = mc2 , bude 2E/(mx2 ) = 2c2 /x2 . Výraz 2/x2 odráží geometrii vesmíru v „comoving souřadnicích a obvykle klademe 2/x2 = −k. Po dosazení do předcházející rovnice dospějeme k základní, tzv. Friedmannově rovnici popisující rozpínání vesmíru pomocí expanzního faktoru Friedmannova rovnice 1 a2 da dt 2 = 8pG 3 + 1 3 Λc2 − kc2 a2 . (1.2.3) Nazvána je po ruském matematikovi A. A. Friedmannovi3 , který v roce 1922 odvodil a předpověděl rozpínání vesmíru z Einsteinových rovnic obecné teorie relativity. Výše uvedené odvození dokazuje platnost rovnice (1.2.3) pro newtonovský model vesmíru, v rámci obecné teorie relativity lze dokázat (viz část 1.3) její platnost pro všechny homogenní izotropní kosmologické modely popsané metrickým tenzorem (1.3.1). 3 V české literatuře se setkáme se dvěma různými přepisy jeho jména: s fonetickým ruským transkriptem „Fridman a nebo s podobou obvyklou v anglosaské literatuře, kterou používáme v textu. Friedmannova maminka, profesionální pianistka, se za svobodna jmenovala Ludmila Vojáčková [1.13]. 13 k > 0 k < 0 k = 0 Obr. 1.7: Dvourozměrné analogie různých geometrií vesmíru v závislosti na hodnotách parametru k (opraveno podle http://map.gsfc.nasa.gov/m_uni/uni_101shape.html). Průvodce studiem Friedmannova rovnice určuje změnu expanzního faktoru s časem, tj. „rychlost rozpínání. Každodenní zkušenost nám napovídá, že k popisu běžných pohybů potřebujeme znát i zrychlení, tj. druhou derivaci expanzního faktoru podle času – pouze s Friedmannovou rovnicí zcela určitě nevystačíme. I když zavedení konstant k a Lambda ve Friedmannově rovnici můžeme chápat formálně, mají svůj fyzikální význam. Těmto otázkám se budeme věnovat ve zbytku podkapi- toly. A. A. Friedmann (1888–1925)Parametr k s rozměrem m−2 charakterizuje konstantní křivost vesmíru a po vhodném přeškálování souřadnic lze uvažovat pouze tři hodnoty k = 0,±1; potom už jde o parametr bezrozměrný a s takovým budeme nadále pracovat [1.17]. Hodnota k = 0 odpovídá nekonečnému plochému vesmíru, v němž platí zákony eukleidovské geometrie (součet úhlu v trojúhelníku je roven 180◦ , obvod kruhu je roven 2pr). Hodnota k = 1 odpovídá uzavřenému konečnému vesmíru, jehož geometrie je podobná geometrii kulové plochy (lze jej v principu „obejít dokola , součet úhlů v trojúhelníku je větší než 180◦ , obvod kruhu menší než 2pr). Konečně hodnota k = −1 odpovídá Parametr k určuje geometrii vesmíru: buď je otevřený, plochý nebo uzavřený. nekonečnému otevřenému vesmíru, jehož geometrie je podobná geometrii sedlové plochy (součet úhlů v trojúhelníku je menší než 180◦ , obvod kruhu větší než 2pr); dvourozměrné analogie zmíněných geometrií vesmíru jsou znázorněny na obr. 1.7. V současné době nelze s jistotou říci, jakou geometrii náš vesmír má, zda je přesně plochý, jak jej vnímáme z běžné každodenní zkušenosti (kdy běžně používáme uvedené vzorce eukleidovské geometrie), nebo je na velkých škálách zakřiven. Víme, že 14 se jeho geometrie od Euklidovy liší jen velmi málo, což teoreticky vysvětlujeme pomocí mechanismu inflace (viz část 1.7). Podrobnější rozbor topologie našeho vesmíru v závislosti na parametru k lze nalézt např. v [1.15, 1.17, 1.27, 1.32]. Druhou rovnici popisující expanzi vesmíru získáme z práce, kterou přitom hmota – „plyn galaxií – vykoná. V části 1.4 ukážeme, že různé formy energie, jež mohou být ve vesmíru obsaženy, přispívají k tlaku „plynu rozdílně a vedou k odlišné závislosti hustoty energie na čase. Uvažujme kouli o objemu V , jejíž objem se změní o dV . Energie obsažená v tomto objemu bude E = 4p 3 c2 R3 . Považujeme-li rozpínání za adiabatické (nedodáváme žádnou vnější energii), musí podle 1. věty termodynamické platit dE + p dV = 0 neboli dE dt + p dV dt = 0. Dosadíme-li postupně dE dt = 4p c2 R2 dR dt + 4p 3 c2 R3 d dt , dV dt = 4pR2 dR dt a přejdeme ke „comoving souřadnicím (1.2.2), dospějeme k další důležité rovnici pro změnu hustoty energie s časem Další důležitá rovnice určuje změnu hustoty vesmíru s časem. d dt + 3 a da dt + p c2 = 0. (1.2.4) Ke změně hustoty energie přispívá jednak zvětšení objemu (člen da/dt, jednak práce vykonaná při expanzi tlakem p, jež se v našem newtonovském modelu mění na gravitační potenciální energii. Derivujeme-li Friedmannovu rovnici (1.2.3) podle času, získáme 2 1 a da dt 1 a2  a d2 a dt2 − da dt 2   = 8pG 3 d dt + 2kc2 a3 da dt , po dosazení za d /dt z (1.2.4) dále 1 a d2 a dt2 − 1 a da dt 2 = − 4pG + p c2 + kc2 a2 a odečtením Friedmanovy rovnice (1.2.3) druhou základní rovnici popisující „zrychlení rozpínání vesmíru Třetí a poslední ze stěžejních rovnic určuje „zrychlení kosmologické expanze. 1 a d2 a dt2 = − 4pG 3 + 3p c2 + 1 3 Λc2 . (1.2.5) 15 Vidíme, že v rozporu s intuitivním předpokladem, že tlak „plynu galaxií se podílí na rozpínání vesmíru, přispívá naopak ke zpomalení expanze podobně jako Λ < 0. Naopak Λ > 0 přispívá ke zrychlení expanze. Tlak p skutečně nemůže vyvolávat rozpínání vesmíru, neboť v homogenním vesmíru nemohou existovat tlakové síly (vyžadovaly by gradient tlaku. Při konkrétních výpočtech se neobejdeme bez stavové rovnice. K řešení získaných rovnic je ještě zapotřebí znát závislost tlaku na hustotě p = = p( ). Závisí na typu materiálu (energie) a nazýváme ji stavovou rovnicí. S nejdůležitějšími případy (z kosmologického hlediska) stavových rovnic se setkáme v části 1.4. Shrnutí Friedmannova rovnice (1.2.3), rovnice pro časovou změnu hustoty energie (1.2.4) a z nich vyplývající rovnice (1.2.5) popisují rozpínání vesmíru prostřednictvím změny expanzního faktoru a(t) s časem. K jejich řešení potřebujeme znát ještě stavovou rovnici, tj, závislost tlaku na hustotě charakterizující typ materiálu (energie) vyplňující vesmír. Geometrie vesmíru se liší v závislosti na parametru k a mohou nastat 3 případy: plochý, otevřený a uzavřený vesmír. Pojmy k zapamatování • comoving souřadnice • kosmologická konstanta • Friedmannova rovnice • stavová rovnice • plochý, uzavřený a otevřený vesmír Kontrolní otázky 1. Jak závisí geometrie vesmíru na parametru k? 2. Uveďte příklady vzorců eukleidovské geometrie, které v zakřiveném prostoru ne- platí. 3. Jakou geometrii má náš vesmír? 4. Jakou závislost udává stavová rovnice? 5. Čím se vyznačují „comoving souřadnice? 6. Jaká vliv má kosmologická konstanta Λ na rozpínání vesmíru? Průvodce studiem Pro čtenáře obeznámeného se základy obecné teorie relativity v následující kapitole naznačíme, jak lze rovnice (1.2.3) a (1.2.5) odvodit z Einsteinových rovnic aplikovaných na homogenní izotropní vesmír. Ostatní čtenáři mohou tuto část při prvním čtení přeskočit, neboť její prostudování není nezbytně nutné k pochopení dalšího výkladu. Dále se budeme odvolávat pouze na Robertsonovu-Walkerovu metriku (1.3.1). 16 1.3 Robertsonova-Walkerova metrika V obecné teorii relativity je základní veličinou popisující geometrii prostoročasu metrický tenzor gµν. Lze ukázat (viz např. [1.15, 1.27]), že předpoklad homogenity a izotropie splňuje Robertsonova-Walkerova metrika Robertsonova- Walkerova metrika popisuje prostoročas s konstantní křivostí. ds2 = −c2 dt2 + a(t)2 dr2 1 − kr2 + r2 dϑ2 + sin2 ϑ dϕ2 , (1.3.1) pro kterou v „comoving souřadnicích x0 = ct, x1 = r, x2 = ϑ a x3 = ϕ má metrický tenzor tvar gµν =      −1 0 0 0 0 a(t)2 /(1 − kr2 ) 0 0 0 0 a(t)2 r2 0 0 0 0 a(t)2 r2 sin2 ϑ      . Parametr k = ±1,0 zavedený v předcházející podkapitole 1.2 charakterizuje křivost prostoročasu, jež v důsledku homogenity a izotropie musí být konstantní. Obsah vesmíru modelujeme ideální tekutinou s tenzorem energie hybnosti (viz např. [1.16, 1.17, 1.26, 1.30]) Tµν = + p c2 uµuν + gµνp. V klidové soustavě tekutiny pro složky čtyřrychlosti platí u0 = −c, ur = uϑ = uϕ = = 0, takže Tµν =      c2 0 0 0 0 pa(t)2 /(1 − kr2 ) 0 0 0 0 pa2 r2 0 0 0 0 pa2 r2 sin2 ϕ      . Po dosazení do Einsteinových rovnic obecné teorie relativity (viz např. [1.6, 1.15, 1.25, 1.36]) Albert Einstein (1879–1955)Gµν + Λgµν = 8pG c4 Tµν (1.3.2) získáme dvě nezávislé rovnice Řešení Einsteinových rovnic vypočtené pro Robertsonovu- Walkerovu metriku. 3 a2   da dt 2 + kc2   − Λc2 = 8pG , − 1 a2  kc2 + 2a d2 a dt2 + da dt 2   + Λc2 = 8pG p c2 . První z nich dává přímo Friedmanovu rovnici (1.2.3), z druhé odečtením (1.2.3) pak získáme (1.2.5). Konečně z relativistického zákona lokálního zachování energie- hybnosti4 ν νTµν = 0, 4 Zákon musíme skutečně chápat v lokálním smyslu, protože pokud se zakřivení prostoročasu mění, energie nezůstává konstantní, ale odpovídajícím způsobem se mění; příkladem může být energie mikrovlnného záření, která v důsledku rozpínání vesmíru klesá (viz vztah (1.4.4)). V plochém prostoročase speciální teorie relativity ovšem zákon zachování platí globálně [1.15]. 17 kde ν značí kovariantní derivaci, obdržíme přímo (1.2.4). Zdůrazněme ještě jednou, že obecná teorie relativity dává hlubší a konzistentnější odvození rovnic, které jsme v části 1.2 získali pomocí newtonovského kosmologického modelu. Moderní relativistická kosmologie představuje jednu z nejúspěšnějších a nejzajímavějších aplikací Einsteinovy teorie relativity. Shrnutí Homogenní izotropní vesmír s konstantní křivostí je v „comoving sférických souřadnicích popsán Robertsonovou-Walkerovou metrikou (1.3.1). Zakřivení prostoru ovlivňuje šíření světelných signálů. Pojmy k zapamatování • Robertsonova-Walkerova metrika Kontrolní otázky 1. Jak lze jednoduše zdůvodnit, že křivost vesmíru musí být konstantní? Průvodce studiem Na konci části 1.2 jsme se zmínili o významu stavové rovnice: charakterizuje typ materiálu, kterým je v našem modelu vesmír vyplněn a který určuje jeho rozpínání. V následující podkapitole probereme nejdůležitější příklady stavových rovnic a jim odpovídajících druhů energie. 1.4 Hustota a tlak Tlak tekutiny v kosmologickém modelu závisí na typu materiálu, jímž je vesmír zaplněn. V zásadě může jít o běžnou hmotu, záření, vakuum, popř. o jejich kombinaci (což je samozřejmě nejrealističtější model); s každým z těchto typů energie je spojen jiný tlak. Vztah p = p( ) mezi tlakem a hustotou popisuje stavová rovnice. Dosazením za p do (1.2.4) pak pro každý typ materiálu získáme předpis, jak se jeho hustota mění s časem. Stavová rovnice má pro uvažované typy energie tvar Stavová rovnice v kosmologii není vůbec složitá! p = w c2 . (1.4.1) Dosazením do rovnice (1.2.4) a integrací získáváme závislost ∝ 1 a3+3w . (1.4.2) Hmotu, tj. materiál z něhož se skládají i naše těla, popisujeme v kosmologickém měřítku jako tzv. prach. Srážky galaxií jsou odlišné od srážek molekul v plynu a 18 nevytvářejí žádný tlak, proto wm = 0.5 Pro závislost hustoty energie hmoty m na čase pak podle (1.4.2) vychází Závislost hustoty hmoty na čase odpovídá změně objemu. m ∝ 1 a3 ; (1.4.3) hustota energie v kouli o poloměru a(t) zůstává konstantní, neboť d d ma3 = 0. Vidíme, že hustota energie hmoty – prachu klesá pouze v důsledku rozpínání vesmíru. Stavovou rovnici pro záření (ideální „plyn tvořený z fotonů) lze odvodit z následující úvahy: uvažujme fotony s hybností p, jež se odrážejí od rovinné stěny s plochou A postavenou kolmo na osu x. Při odrazu od stěny se změní složka hybnosti fotonu px na −px, ostatní složky zůstanou stejné. Průměrná síla, kterou působí fotony na stěnu, je dána změnou jejich hybnosti v uvažovaném časovém intervalu, tj. 2px násobek počtu fotonů, které v tomto intervalu do stěny narazí. Uvažujeme-li fotony s různými hybnostmi p a uvědomíme-li si, že v průměru se pouze polovina pohybuje v kladném směru osy x (polovina v opačném) a tlak je určen podílem síly a plochy, na kterou síla působí, pro tlak na stěnu A vychází pγ = 1 2 2 pxvx nγ, kde nγ je počet fotonů v jednotkovém objemu (koncentrace) a čárkou „nad značíme střední časovou hodnotu příslušného výrazu. Fotony se pohybují rychlostí c a pohyby ve směru os x, y i z jsou stejně pravděpodobné, takže v průměru se ve směru osy x pohybuje pouze 1/3 z jejich celkového počtu, takže pγ = 1 3 c |p| nγ. Připomeňme, že energie Eγ fotonu o vlnové délce λ je rovna Eγ = h c λ = c |p| , kde h je Planckova konstanta. Pro záření tak dostáváme stavovou rovnici6 S zářením je spojen nenulový tlak!pγ = 1 3 γc2 tj. wγ = 1/3. Dosazením do (1.4.2) obdržíme Hustota záření s časem klesá rychleji než hustota prachu. γ ∝ 1 a4 . (1.4.4) 5 Srážky galaxií znamenají v podstatě jejich prostoupení, nikoli odraz pružných „kuliček ; výsledkem jsou – v závislosti na hmotnosti interagujících galaxií - většinou různé nepravidelné galaxie nebo galaktický kanibalismus – pohlcení jedné galaxie druhou. 6 Přesné odvození stavové rovnice ideálního fotonového plynu lze nalézt ve většině učebnic statistické fyziky, např. [1.10, 1.18, 1.21, 1.22] 19 Hustota energie záření tak klesá rychleji než hustota hmoty podle rovnice (1.4.3), protože se nejenom mění počet fotonů v jednotkovém objemu v důsledku rozpínání vesmíru, ale navíc se mění i jejich energie v důsledku rudého posuvu popsaného v části 1.6. Ze statistického rozboru záření černého tělesa popsaného Planckovým zákonem a Wienovým posunovacím zákonem (viz např. [1.10, 1.18, 1.21, 1.22]) plyne, že nejvíce energie připadá na fotony s energií E ≈ 2,8 kBT a střední energie fotonů záření při teplotě T je E ≈ 2,7 kBT. Při rozpínání vesmíru pak vlnová délka λ, energie fotonu E a rovnovážná teplota záření T závisejí na expanzním faktoru a(t) λ ∝ a(t), E ∝ 1 λ ∝ 1 a(t) , T ∝ E ∝ 1 a . Počet fotonů se přitom zachovává, takže jejich počet nγ v objemové jednotce musí při rozpínání klesat nγ ∝ 1 a(t)3 , zatímco energie fotonů v této jednotce se mění podle vztahu (1.4.4) γ = nγE ∝ T4 v souladu se Stefanovým-Boltzmannovým zákonem. Konečně hustota energie vakua se s časem nemění, tj. Hustoty vakua na čase nezávisí: vakuum je prostě vakuum! d v dt = 0, (1.4.5) tudíž d dt va3 = v da3 dt . Z rovnice (1.2.4) pak vychází Stavová rovnice vakua nám ukazuje, že energie vakua je spojena s tlakem. pv = − vc2 neboli wv = −1. Výsledek je v souladu s požadavkem, že hustota energie vakua musí být stejná pro všechny pozorovatele, tj. invariantní pro vzhledem k Lorentzově transformaci speciální teorie relativity. V homogenním vesmíru musí být konstantní v prostoru, a kombinace časové i prostorové derivace v Lorentzových transformacích vyžaduje, aby nezávisela ani na čase (důkaz invariantnosti a odpovídající tvar tenzoru energie-hybnosti vakua lze nalézt např. v [1.19]). V poslední době byla předložena řada hypotéz, předpokládajících, že vesmír obsahuje speciální látku zvanou „quintessence , pro niž w < −1/3, popř. závisí na Zda existuje „quintessence nebo jde o pouhou spekulaci nelze dnes ještě rozhodnout. čase. Důvodem je snaha vysvětlit zrychlující se rozpínání vesmíru, jež by existence takové – dosud experimentálně neprokázané – formy hmoty vysvětlovala. Ze závislostí (1.4.3), (1.4.4) a (1.4.5) vyplývá, že vzájemný poměr jednotlivých druhů energie se s časem mění. Hustota energie záření klesá rychleji, než hustota energie hmoty; v raných stádiích vývoje vesmíru proto bylo dominantní záření. Je-li vakuum nositelem energie, bude v někdy v budoucnu dominantní právě vakuum. V současné době jsme v etapě rozpínání vesmíru, kde je dominantní nebo alespoň velmi významnou složkou „prach (chladná hmota). 20 Jiná dominantní energie má za následek jinou závislost expanzního faktoru na a(t) na čase. Kvalitativní rozdíl můžeme demonstrovat na nejjednodušším modelu s k = 0 i Λ = 0. Friedmannova rovnice (1.2.3) se potom zjednoduší na da dt 2 ∝ a2 . Pro prach po dosazení z (1.4.3) s počáteční podmínkou a(t = 0 = 0 vychází da dt ∝ 1 a1/2 =⇒ a ∝ t2/3 ; (1.4.6) podobně pro záření podle (1.4.4) da dt ∝ 1 a =⇒ a ∝ t1/2 (1.4.7) a pro vakuum v souladu s (1.4.5) da dt ∝ a =⇒ a ∝ exp (Ht) , (1.4.8) kde H je Hubbleův parametr definovaný vztahem (1.5.2), v tomto případě je kon- Inflací rozumíme rychlé a obrovské rozpínání vesmíru způsobené energií vakua. stantní. Exponenciální rozpínání způsobené energií vakua nazýváme inflací.. Ve skutečném vesmíru se uplatňují všechny zmíněné složky, ale v různých fázích vývoje různou měrou. Na počátku vývoje prošel fází dominantního záření, nyní dominuje prach a v daleké budoucnosti bude pravděpodobně rozhodující energie vakua. Přechod mezi érami dominantního záření a dominantního prachu pro model s k = 0 a Λ = 0 je znázorněn na obr. 1.8. Shrnutí ln t ln γ m γ > m m > γ Obr. 1.8: Éry dominantního záření a prachu Vesmír je naplněn různými typy energie, z nichž nejvýznamnější vliv na jeho expanzi mají chladná hmota (prach), záření a energie vakua. Jejich vzájemné zastoupení se s časem mění; v dávné minulosti prošel vesmír érou dominantního záření, nyní se nachází ve stavu dominantního prachu a ve vzdálené budoucnosti pravděpodobně bude dominovat energie vakua. Předpokládá se, že energie vakua dominovala i v krátkém časovém intervalu (řádově t = 10−34 s po velkém třesku). Tuto fázi označujeme jako inflace a vesmír se během ní exponenciálně zvětšil. Pojmy k zapamatování • quintessence • inflace 21 Kontrolní otázky 1. Co se kromě rozpínání vesmíru podílí na změně energie záření? 2. Jaké typy energie dominovaly v různých fázích vývoje vesmíru? 3. K čemu dochází ve fázi vývoje zvané inflace? Úkoly k textu 1. Dokažte, že pokud ve vesmíru dominuje energie popsaná obecnou stavovou rovnicí (1.4.1), bude při k = 0 a Λ = 0 expanzní faktor záviset na čase podle vztahu a ∝ t2/(3+3w) . Ověřte, že rovnice (1.4.6), (1.4.7) i (1.4.8) jsou speciálními případy získané závislosti. 2. Ukažte, že pokud k = 0 a rozpínání vesmíru u ovlivňuje pouze prach s hustotou energie m > 0 a vakuum s hustotou energie v > 0, potom (viz [1.31]) a3/2 ∝ sinh 6pG vt . Průvodce studiem Po odvození základních rovnic můžeme nyní přistoupit k podrobnější diskusi vlastností našich matematických modelů a jejich konfrontaci s výsledky pozorování reálného vesmíru. V následující podkapitole se zaměříme na ústřední pozorovanou skutečnost, která ve své době překvapila i samotného Einsteina: náš vesmír se rozpíná. . . 1.5 Hubbleův zákon Označíme-li v rychlost, jíž se od nás vzdaluje galaxie ve vzdálenosti d = a(t)r, kde r reprezentuje „comoving souřadnice, lze psát v = dd dt = da dt r = 1 a da dt d neboli Hubbleův zákon: jednoduchý vztah mezi rychlostí a vzdáleností. v = Hd, (1.5.1) kde jsme zavedli časově proměnný parametr H = H(t) = 1 a da dt . (1.5.2) 22 Vztah (1.5.1) nazýváme Hubbleovým zákonem a parametr H(t) definovaný rovnicí (1.5.2) Hubbleovým parametrem. Současná hodnotu parametru H0 = H(t0) představuje jeden z důležitých observačních údajů – Hubbleovu konstantu7 . Jednoduchý vztah (1.5.1) vyjadřuje skutečnost, že čím je od nás objekt dále, tím se od nás více vzdaluje a tím je světlo, které od objektu přichází, posunuto k delším vlnovým délkám; vzhledem k proměnnosti Hubbleova parametru H(t) s časem však nejde o přímou úměru. Hubbleův zákon představuje přímý experimentální důkaz o rozpíE. P. Hubble (1889 - 1953) nání našeho vesmíru a poprvé byl ověřen Edwinem Powellem Hubblem v roce 1929 pomocí 100-palcového dalekohledu observatoře na Mt. Wilsonu v USA. Většinou počítáme s hodnotou Hubbleovy konstanty H0 = 100 h km·s−1 ·Mpc,−1 kde h ∈ 0,55 − 0,75. (1.5.3) Bezrozměrný parametr h, odrážející neurčitost v současném měření Hubbleovy Hubbleova konstanta je současnou hodnotou Hubbleova parametru. konstanty, slouží jako „měřítko mapy pro kosmologické vzdálenosti, proto bývá výhodné pracovat s poměrnými veličinami, v nichž se vykrátí (Hubbleův parametr (1.5.2), decelerační parametr (1.7.7)). Ze studia anizotropií kosmického mikrovlnného záření vychází hodnota h = 0,71 ± 0,04 [1.5], k podobným výsledkům konvergují i další nezávislá pozorování [1.12]. I když se Hubbleova konstanta tradičně udává v km·s−1 ·M−1 pc, v soustavě SI má rozměr s−1 . Její převrácenou hodnotu nazýváme Hubbleovým časem tH . V miliardách let (Gyr) vychází tH = 1 H0 = 9,78 h−1 Gyr. (1.5.4) Hubbleův čas představuje hrubý odhad stáří vesmíru – tak dlouho by trval, pokud Hubbleův čas představuje hrubý odhad stáří vesmíru. by se po celou dobu rozpínal rovnoměrně (srovnejte se vztahem (1.7.9)). Shrnutí Rychlost, s jakou se od nás vzdalují galaxie, se zvětšuje s jejich vzdáleností od nás podle Hubbleova zákona (1.5.1). Hubbleův parametr vystupující v této rovnici jako koeficient úměrnosti závisí obecně na čase. Jeho současnou hodnotu nazýváme Hubbleovou konstantou a patří k nejdůležitějším experimentálním údajům o našem vesmíru. Její převrácená hodnota se nazývá Hubbleovým časem a slouží jako hrubý odhad stáří vesmíru. Pojmy k zapamatování • Hubbleův zákon • Hubbleův parametr • Hubbleova konstanta • Hubbleův čas 7 V kosmologii obvykle indexem „0 značíme současné hodnoty veličin. 23 Kontrolní otázky 1. Jaká je podle současných pozorování hodnota Hubbleovy konstanty? 2. Za jakých podmínek by stáří vesmíru bylo rovno Hubbleovu času? Cvičení 1. Najděte závislost H = H(t) Hubbleova parametru na čase pro éru dominantního záření a dominantního prachu s využitím rovnic (1.4.7) a (1.4.6). Uvažujte nejjednodušší případ plochého vesmíru k = 0 bez energie vakua Λ = 0. Řešení 1. Z rovnic (1.2.3) a (1.4.3) pro éru dominantního prachu plyne H2 = 1 a2 da dt 2 ∝ 1 a3 neboli √ a da ∝ dt. Integrací s počáteční podmínkou a(t = 0) = 0 pak vychází a3/2 ∝ t =⇒ a ∝ t2/3 . Označíme-li konstantu úměrnosti K, bude a = Kt2/3 , a0 = Kt2/3 , =⇒ a = a0 t t0 2/3 , kde indexem 0 značíme současné hodnoty. Dosazením posledního výrazu do (1.5.2) vychází H = 2 3t , konkrétně H0 = 2/(3t0). Obráceně pak pro stáří takového vesmíru vychází t0 = 2 3H0 = 2 3 tH, kde tH = 1/H0 je Hubbleův čas. Pro model s dominantním zářením se výchozí rovnice změní na H2 = 1 a2 da dt 2 ∝ 1 a4 resp. a da ∝ dt. Analogickým způsobem pak odvodíme a = a0 t t0 1/2 , H = 1 2t . 24 Průvodce studiem V části 1.4 jsme viděli, že hustota energie záření při expanzi vesmíru klesá rychleji než hustota prachu a tudíž i rychleji, než by odpovídalo pouhému zvětšení objemu. Proč tomu tak je? V následující podkapitole odhalíme příčinu: při rozpínání vesmíru se mění vlnová délka (a tím samozřejmě i frekvence) světelných signálů. V kombinaci s výše popsaným Hubbleovým zákonem to znamená, že čím je pozorovaná galaxie od nás dále, tím je změna vlnové délky větší. Připomeňme, že pozorování galaxií je zároveň výletem do jejich minulosti, neboť na překonání vzdálenosti od nich k nám potřebuje světlo jistý čas (viz část 1.9). 1.6 Šíření světla ve Friedmannových modelech Při popisu šíření světelných signálů musíme vzít v úvahu geometrii vesmíru charakterizovanou metrickým tenzorem. Studujme pro jednoduchost radiální šíření světelných signálů v Robertsonově-Walkerově prostoročase (1.3.1). Světlo se podle obecné teorie relativity šíří po nulových geodetických čarách, pro něž platí ds2 = 0, pro čistě radiální pohyb ze zdroje o souřadnici r = re k pozorovateli o souřadnici r = r0 jsou navíc ϑ i ϕ konstantní, tj. dϑ = dϕ = 0. Dosazením do (1.3.1) získáme Pro šíření světelných signálů položíme v (1.3.1) ds = 0. cdt a(t) = ± dr √ 1 − kr2 , kde kladné resp záporné znaménko odpovídá fotonům šířícím se ve směru rostoucího resp. klesajícího r. Po integraci dostáváme t0 te cdt a(t) = ± r0 re dr √ 1 − kr2 . (1.6.1) Označíme-li periody v místě vyslání a zachycení signálu po řadě Te a T0 a jim odpovídající vlnové délky λe = cTe a λ0 = cT0, bude pro signál emitovaný v čase te + Te analogicky platit t0+T0 te+Te cdt a(t) = ± r0 re dr √ 1 − kr2 . (1.6.2) Vzhledem k tomu, že r je „comoving souřadnicí, která se nemění s časem, budou na čase nezávislé pravé strany rovnic (1.6.1) a (1.6.2), musí být proto rovny i jejich levé strany t0+T0 te+Te cdt a(t) = t0 te cdt a(t) . S využitím vlastností mezí určitých integrálů po úpravě dostáváme te te+Te cdt a(t) + t0 te cdt a(t) + t0+T0 t0 cdt a(t) = t0 te cdt a(t) , 25 t0+T0 t0 cdt a(t) = te+Te te cdt a(t) . V malých časových intervalech odpovídajících periodě světelných vln Te a T0 lze pokládat expanzní faktor a(t) za konstantní, tudíž Vlnová délka λ je úměrná expanznímu faktoru a. cT0 a(t0) = cTe a(te) , neboli λ0 a(t0) = λe a(te) . Vlnová délka světla se mění úměrně s expanzním faktorem. Protože a(t) s časem roste, zvětšuje se i vlnová délka směrem k červenému konci spektra. V praxi charakterizujeme změnu vlnové délky nejčastěji tzv. rudým posuvem z definovaným vztahem Změnu vlnové délky vyjadřuje rudý posuv. z = λ0 − λe λe (1.6.3) neboli 1 + z = λ0 λe = a(t0) a(te) . (1.6.4) Příklady spekter galaxií s různými rudými posuvy jsou na obr. 1.9. Naměření rudého posuvu vzdálených galaxií je přímým důkazem rozpínání našeho vesmíru. Např. signál, pro který naměříme z = 1, byl vyslán v čase te, kdy a(te) = a(t0)/2, tj. vesmír měl poloviční velikost než dnes. Nejjednodušší vysvětlení rudého posuvu předpokládá, že se vlnová délka světla zvětšuje spolu s vesmírem (obr. 1.10). Nabízí se i vysvětlení, že jde vlastně o dopplerovský posun způsobený vzdalováním zdrojů světla v důsledku rozpínání vesmíru. Odvozený vzorec pro rudý posuv z (1.6.4) je v souladu se vztahem pro Dopplerův posuv ve speciální teorii relativity pouze v prázdném prostoru, prakticky pro dostatečně blízké zdroje (z = = v/c Ć 1). V zakřiveném prostoročase je korektní zavedení relativní rychlosti zdroje a pozorovatele komplikovanější. Obecně relativistický výpočet v [1.28] dokazuje, že kosmologický rudý posuv lze skutečně interpretovat jako dopplerovský, i když v literatuře není na tuto otázku jednotný názor. Shrnutí Šíření světla odráží geometrii prostoročasu. V rozpínajícím se homogenním a izotropním vesmíru roste vlnová délka světla úměrně velikosti vesmíru. Mírou změny vlnové délky je rudý posuv spektrálních čar. Pojmy k zapamatování • rudý posuv Kontrolní otázky 1. Dokážete na základě definice rudého posuvu vysvětlit modrý posuv? 2. Pro které objekty ve vesmíru můžeme naměřit rudý a modrý posuv? 3. Jaké je vysvětlení kosmologického rudého posuvu? 26 λ z = 0,067 λ z = 0,264 Obr. 1.9: Příklady spekter dvou galaxií s různými rudými posuvy. Vlnová délka je udávána v angströmech, 1 ˚A = 0,1 nm. Vyznačeny jsou i typické spektrální čáry; všimněme si různé polohy těchto čar jež odráží různý rudý posuv. Svislá osa znázorňuje relativní světelný tok připadající na jednotlivé části spektra. Upraveno podle [1.8]. ˇcas Obr. 1.10: K vysvětlení rudého posuvu: vlnová délka fotonu se „rozpíná s vesmírem jako u fotonu v uzavřené krabici Cvičení 1. Jaký rudý posuv naměříme u galaxií patřících do kupy galaxií v souhvězdí Panny, vzdalují-li se od nás rychlostí asi 1 000 km·s−1 ? Rychlost světla je přibližně 300 000 km·s−1 [1.39]. 27 2. Které z emisních čar v následující tabulce můžeme z povrchu Země pozorovat v optickém oboru spektra u kvasaru s následujícím rudým posuvem: a) z = 0,1, b) z = 1, c) z = 4. Tabulka hlavních emisních čar u aktivních galaxií a kvasarů [1.38]: Lα 121,6 nm N V 124,0 nm C IV 154,9 nm C III 190,9 nm Mg II 279,8 nm O II 372,7 nm Ne III 386,8 nm Hδ 410,2 nm Hγ 434,1 nm Řešení 1. Pro dané hodnoty v = 1 000 km·s−1 , c = 300 000 km·s−1 , tedy v/c = 1/300 ≈ 0,0033 Ć 1. Rychlost je dostatečně malá, takže můžeme použít přibližný vztah z ≈ v c ≈ 0,0033. 2. Označme za = 0,1, zb = 1 a zc = 4 rudé posuvy. Viditelné světlo má vlnové délky v intervalu I0 = 400 nm,780 nm . Podle (1.6.3) platí z = λ0 λe − 1, kde λ0 je vlnová délka pozorovaná na zemi a λe vlnová délka světla vyzářeného kvazarem. Zřejmě λe = λ0 1 + z . Přepočítejme podle posledního vztahu intervaly vlnových délek, v nichž může být při daných rudých posuvech vyzářeno světlo, aby jeho vlnová délka λ0 na Zemi zapadala do intervalu viditelného světla: a) Iza = I0/(1 + za) = 363,6 nm,709,1 nm , b) Izb = I0/(1 + zb) = 200,0 nm,390,0 nm , c) Izc = I0/(1 + zc) = 80,0 nm,156,0 nm . Porovnáním získaných intervalů se zadanou tabulkou zjistíme, že viditelné budou spektrální čáry a) O II, Ne III, Hδ a Hγ, b) Mg II, O II a Ne III, c) Lα, N V a C IV. 28 Průvodce studiem Známe již klíčové rovnice (1.2.3)–(1.2.5) a přirozeně se ptáme: jaká byla minulost a jaká může být daleká budoucnost právě našeho jedinečného vesmíru? Zdaleka ne všechny parametry vystupující ve zmíněných rovnicích umíme snadno určit – vždyť už v části 1.5 jsme viděli, že Hubbleova konstanta není známa s velkou přesností. Pro řadu úloh je výhodnější pracovat s bezrozměrnými, tzv. observačními parametry, jejichž hodnoty lze určit přímo z pozorování. Zavedení těchto veličin je věnována následující podkapitola. 1.7 Observační parametry vesmíru Ze všech Friedmannových modelů je nejjednodušší případ s k = 0 a Λ = 0. Odpovídá mu nekonečný vesmír s eukleidovskou geometrií, který se bude neustále rozpínat. Z Friedmannovy rovnice (1.2.3) pro hustotu vesmíru v tomto modelu dostáváme Kritická hustota určuje hustotu plochého vesmíru obsahujícího pouze prach. c = 3H2 8pG (1.7.1) označovanou jako kritická hustota. Vzhledem k tomu, že Hubbleův parametr H se mění s časem, není ani hodnota kritické hustoty neměnná. Pro současnou hodnotu (1.5.3) dostáváme c0 = 3H2 0 8pG = 1,88 h2 ·10−26 kg·m−3 = 2,78 h−1 ·1011 M⊗/ h−1 Mpc . (1.7.2) Vidíme, že ačkoli v běžně užívaných jednotkách jde o hustotu velmi malou, je průměrná hustota našeho vesmíru této hodnotě velmi blízká, i když lokálně (např. v našem těle, na Zemi, ve Sluneční soustavě) může být mnohem větší. Galaxie totiž obsahují miliardy hvězd a jejich hmotnost se pohybuje okolo 1011 –1012 hmotností Slunce M a 1 Mpc je zároveň typickou mezigalaktickou vzdáleností. Jsou i další argumenty, např. inflační modely, podle nichž se hustota našeho vesmíru od c0 liší jen velmi málo. Hustotní parametry udávají poměr hustot vůči c. Pro jednotlivé druhy energie (prach, záření i vakuum) můžeme zavést hustotní parametr prachu Ωm (t) = m c , (1.7.3) záření Ωγ (t) = γ c (1.7.4) nebo vakua Ωv (t) = v c = 8pG v 3H2 = Λc2 3H2 . (1.7.5) Protože prach i záření hrají v dynamice kosmologických modelů obdobnou úlohu a jednu z těchto složek považujeme většinou za dominantní (po většinu doby expanze 29 je to právě prach), budeme jejich součet nadále označovat jako Ωm. Dosazením do Friedmannovy rovnice (1.2.3) vychází kc2 H2a2 = Ωm + Ωv − 1. (1.7.6) Vidíme, že pokud má mít vesmír plochou eukleidovskou geometrii k = 0, musí být splněno Ωm + Ωv = 1. Pokud se rozpínání vesmíru zpomaluje, potom se součin H2 a2 = (da/dt)2 zmen- Inflace je jednou z možností jak vysvětlit, proč je náš vesmír (téměř) plochý. šuje a zlomek na levé straně rovnice (1.7.6) zvětšuje. Je-li odchylka Ωm + Ωv od 1 nyní malá, musela být v minulosti ještě mnohem, mnohem menší, hovoříme o tzv. problému plochosti vesmíru: jako to, že je geometrie vesmíru vyladěna tak blízko k eukleidovské? Tato skutečnost je jednou z motivací k předpokladu, že vesmír ve velmi raném stadiu prošel fází tzv. inflace, prudkého rozpínání s H ≈ konst., kdy se hodnota zlomku zmenšila úměrně 1/a2 a vesmír se stal (téměř) plochým. Potvrzuje to i rozbor anizotropií mikrovlnného záření s výsledkem Ωm +Ωv = 1,02±0,02 [1.5]. Decelerační parametr charakterizuje rozpínání vesmíru na větších vzdálenostech. Důležitou charakteristikou expanze vesmíru je kromě rychlosti rozpínání, jejíž mírou je Hubbleova konstanta, také „zrychlení . Až do 90-tých let minulého století se předpokládalo, že expanze se bude vždy zpomalovat, namísto akcelerace bylo proto přirozenější zavést veličinu zvanou decelerační parametr q definovaný vztahem q = − d2 a dt2 a da dt 2 . (1.7.7) Dosazením do (1.2.5) s využitím (1.7.3)–(1.7.5) pak vychází q = Ωm 2 − Ωv. (1.7.8) Zřejmě pokud Ωm/2 − Ωv = 0, rozpínání probíhá s konstantní rychlostí, pro q > 0 se zpomaluje a pro q < 0 se zrychluje. Rozvineme-li vztah pro expanzní faktor a(t) v Taylorovu řadu okolo současného stáří vesmíru t = t0, dostáváme a(t) = a(t0) + da dt t0 (t − t0) + d2 a dt2 t0 (t − t0)2 + . . . a pro rudý posuv objektu, od něhož bylo světlo vysláno z čase t podle (1.6.4) vychází 1 1 + z = a(t) a(t0) = 1 + H0 (t − t0) − q0 2 H2 0 (t − t0)2 + . . . , kde H0 = H(t0), q0 = q(t0) představují současné hodnoty Hubbleova a deceleračního parametru. 30 Z Friedmannovy rovnice (1.7.6) pak lze vyjádřit současnou hodnotu expanzního faktoru a0 = c H0 k Ωm0 + Ωv0 − 1 . V případě plochého vesmíru k = 0 vychází neurčitý výraz 0/0, pro který z limitního chování obdržíme tzv. Hubbleovu délku Hubbleova délka představuje hrubý odhad velikosti pozorovaného vesmíru. dH = ctH = c H0 ≈ 2 998 h−1 Mpc, (1.7.9) která odpovídá dráze, kterou urazí světlo za Hubbleův čas (1.5.4) v plochém vesmíru. Přibližně tak charakterizuje velikost pozorovaného vesmíru, tj. oblasti, z níž můžeme zachytit světelné (i jakékoli jiné) signály. Shrnutí Kritickou hustotou nazýváme hustotu plochého homogenního vesmíru s nulovou energií vakua. Hustotu všech typů energie pak s výhodou vyjadřujeme poměrem k této kritické hustotě pomocí hustotních parametrů. Spolu s deceleračním parametrem, popisujícím zrychlování resp. zpomalování expanze, představují další důležité měřitelné charakteristiky vesmíru. Pojmy k zapamatování • hustotní parametr • decelerační parametr • Hubbleova délka Kontrolní otázky 1. Jaký vztah musí splňovat hustotní parametry pro plochý vesmír? 2. Jak moc se liší celková hustota vesmíru od kritické a jakým procesem to umí moderní kosmologie vysvětlit? 3. Jaký je význam deceleračního parametru? 4. Jaký má význam Hubbleova délka? Cvičení 1. Kolika nukleonům v m3 odpovídá kritická hustota vesmíru pro hodnotu Hubbleovy konstanty H0 = 71 km·s−1 ·Mpc [1.39]? Úkoly k textu 1. Ukažte, že je-li v nějakém čase Ω = Ωm + Ωv = 1, potom se Ω s časem nemění. 2. Jaká podmínka musí být splněna, aby hustotní parametr vakua Ωv byl konstantní a nezávisel na čase? 3. Předpokládejme, že v určitém čase je v rovnici (1.7.8) q = 0. Co musí být splněno, aby decelerační parametr zůstal nulový? 31 Řešení 1. Kritickou hustotu určíme ze vztahu (1.7.1) dosazením za příslušné konstanty H0 = 71 km·s−1 ·Mpc = 2,30·10−18 s−1 a G = 6,67·10−11 m3 ·kg·−1 s2−1 . Je-li hmotnost jednoho nukleonu (protonu nebo neutronu) mu ≈ 1,66·10−27 kg (tzv. hmotnostní atomová jednotka), bude odpovídající počet nukleonů v m3−1 n = 3H2 0 8pGmu ≈ 5,7. Kritická hustota odpovídá necelým šesti nukleonům v 1 m3 . Průvodce studiem Jedna z ústředních otázek kosmologie zní: bude se náš vesmír stále rozpínat nebo se někdy v budoucnu smrští zpět do stavu podobného jako na počátku? Odpověď dávají hodnoty hustotních parametrů, jejichž přesné určení představuje výzvu mnoha experimentálním projektům (Boomerang, 2dF Galaxy Redshift Survey, Sloan Digital Sky Survey, Supernova Cosmology Project, Wilkinson Anisotropy Probe). Konktrétní možnosti vývoje vesmíru v závislosti na obsahu množství hmoty a energie vakua diskutujeme v následující části. 1.8 Budoucnost vesmíru Moderní kosmologie předpokládá, že vesmír se rozpíná z velmi horkého, superhustého stavu, označovaného jako velký třesk nebo anglicky „Big Bang Jednou ze zajímavých otázek, které můžeme zkoumat je, zda rozpínání vesmíru bude trvat věčně nebo zda se zastaví, vesmír dosáhne maximálního rozměru a poté se začne znovu smršťovat a skončí „velkým křachem (anglicky „Big Crunch ). Různé případy, jež mohou nastat, přehledně znázorňujeme v rovině Ωm a Ωv (záření, které dominovalo krátkou dobu na počátku vesmíru neuvažujeme). Rozpínání se zřejmě zastaví, pokud d2 a/dt2 < 0 (neboli q > 0) a tato podmínka zůstává splněna, dokud da/dt = 0 (resp. H = 0). Podle (1.7.8) tato situace vždy nastane pro Ωv < 0, neboť potom Je-li Ωv < 0, expanze se jednou změní ve smršťování. q = Ωm 2 − Ωv > 0 a ve Friedmannově rovnici (1.7.6) lze nalézt a, pro něž H = 0. Jestliže Ωm 5 1, potom se rozpínání zastaví pro Ωv < 0, ale nikoli pro Ωv = 0. Protože kladná hodnota Ωv zvětšuje jak d2 a/dt2 , tak da/dt, vidíme, že v tomto případě se rozpínání ve smršťování také nezmění. Konečně při Ωm = 1 také existuje maximální hodnota Ωv, při které se rozpínání ještě zastaví a přejde v kontrakci (viz obr. 1.11). V tomto limitním případě dospěje vesmír do stavu, kdy q = 0 i H = 0, který odpovídá tzv. Einsteinovu statickému 32 vesmíru8 (v němž by tyto podmínky platily po celou dobu). Úpravou Friedmanovy Při Ωv > 0 je teoretických možností více, dokonce i nerozpínající se Einsteinův statický vesmír. rovnice (1.7.6) se započtením závislostí hustoty energie prachu na expanzním faktoru (1.4.3)–(1.4.5) lze psát H2 = H2 0 Ωm0 a3 0 a3 + H2 0 Ωv0 − kc2 a2 . Z této rovnice lze pro t = t0 vyjádřit kc2 = H2 0 a2 0 (Ωm0 + Ωv0 − 1) a dosadit zpět, čímž dostaneme H2 H2 0 = Ωm0 a3 0 a3 − a2 0 a2 + Ωv0 1 − a2 0 a2 + a2 0 a2 = 0. (1.8.1) Podobně podmínka q = 0 podle (1.7.8) dává Ωv0 = 4Ωm0 a0 2a 3 . Dosazením do (1.8.1) a zavedením substituce a = a0/(2x) dospějeme ke kubické rovnici 4x3 − 3x + 1 − 1 Ωm0 = 0. (1.8.2) Podle definice x < 1, rovnici proto lze řešit pomocí goniometrické substituce x = = cos β. S využitím vzorce (viz např. [1.35]) cos 3β = 4 cos3 β − 3 cos β obdržíme z (1.8.2) postupně cos 3β = − 1 − 1 Ωm0 , cos (3β − p) = 1 − 1 Ωm0 , x = cos 1 3 arccos 1 − 1 Ωm0 + p 3 , Ωv0 = 4Ωm0x3 = 4Ωm0 cos3 1 3 arccos 1 − 1 Ωm0 + p 3 . (1.8.3) Poslední závislost (ve zvětšeném měřítku) je také vykreslena na 1.11. Podrobnější rozbor této problematiky lze nalézt např. v [1.11]. 8 Kosmologická konstanta Λ byla Einsteinem zavedena právě proto, aby získal statický, nerozpínající se model vesmíru. Po Hubbleových pozorováních, která rozpínání vesmíru prokázala, označil Einstein zavedení Λ jako jeden ze svých největších omylů, dnes se však zdá být omylem velmi inspirujícím. . . 33 1 20 1 2 3 -1 0 1 2 3 2 3 CMB Clusters Knop et al. (2003) Spergel et al. (2003) Allen et al. (2002) Supernova Cosmology Project Ω Ωv m bez velk´eho tˇresku Pozorov´an´ı supernov st´ale se rozp´ın´a ploch´y uzavˇren´y otevˇren´y historie Vesm ´ıru kr´atk´a pˇr´ıliˇs m˚uˇze se smrˇst’ovat rozp´ın´an´ı se zrychluje rozp´ın´an´ı se zpomaluje Obr. 1.11: Znázornění výsledků tří různých nezávislých měření v rovině parametrů Ωm0 a Ωv0: pozorování vzdálených supernov [1.20], kosmického mikrovlnného záření (CMB) [1.37] a studia rozložení hmoty ve skupinách galaxií [1.1]. Vidíme, že pozorování konvergují k hodnotám Ωm0 = 0,3 a Ωv0 = = 0,7. Upraveno podle [1.33]. Shrnutí Budoucnost vesmíru – věčné rozpínání versus budoucí smršťování, zrychlování versus zpomalování expanze – lze jednoznačně předpovědět v závislosti na hodnotách hustotních parametrů hmoty a vakua. Pro jejich současné nejpravděpodobnější hodnoty se vesmír bude stále rychleji rozpínat. Pojmy k zapamatování • velký třesk • velký křach • Einsteinův statický vesmír Kontrolní otázky 1. Jaké jsou nejpravděpodobnější hodnoty hustotních parametrů Ωm0 a Ωv0? 2. Co z těchto hodnot plyne pro budoucnost Vesmíru? Úkoly k textu 1. Ověřte, že pokud Ωm0 = 54/28, potom z rovnice (1.8.2) vychází Ωv0/Ωm0 = = 1/54. 34 2. Ukažte, že pokud Ωm0 − 1 Ć 1, potom aproximací rovnice (1.8.3) získáme Ωv0 = 4Ωm0 27 1 − 1 Ωm0 3 , tj. stejný výsledek, jako kdybychom v rovnici (1.8.2) zanedbali kubický člen. Průvodce studiem Již v souvislosti s rudým posuvem jsme připomínali, že pozorování vesmíru je cestou do minulosti. Zejména u velmi vzdálených objektů je důležité určit dobu, kterou jejich světlo potřebovalo k uražení vzdálenosti k nám. Pozorujemeli dnes objekty, jejichž světlo bylo vyzářeno před 10 a více miliardami let (viz obr. 1.5), získáváme konkrétnější představu o formování galaxií, formování prvních hvězd apod. Díky rozpínání vesmíru a možné křivosti jeho geometrie však určování časů a vzdáleností není jednoduché. Nejpraktičtější je najít jejich vztah k experimentálně dobře definovaným veličinám: rudému posuvu z a hustotním parametrům. 1.9 Stáří pozorovaných objektů V této části odvodíme vztah mezi rudým posuvem z světla vyslaného vzdálenou galaxií a časem t, který světelný signál potřeboval k překonání vzdálenosti k nám. Určování vzdálenosti pozorovaných objektů se pak budeme věnovat v následující podkapitole 1.10. Obě veličiny musí odrážet geometrii vesmíru a nutně proto závisejí na hustotě hmoty Ωm (prachu popř. záření) i hustotě vakuové energie Ωv. Vyjdeme ze zavedení Hubbleova parametru (1.5.2), jež přepíšeme do tvaru Vyjádříme změnu času v závislosti na expanzním faktoru... dt = 1 H da a . (1.9.1) Ze vztahu (1.8.1) dále plyne H = H0 Ωm0 a3 0 a3 − a2 0 a2 + Ωv0 1 − a2 0 a2 + a2 0 a2 a z definice rudého posuvu (1.6.4) také ...a expanzní faktor v závislosti na rudém posuvu. a = a0 1 + z , da = − a0 (1 + z)2 dz, kde a0 značí současnou hodnotu expanzního faktoru. Označíme-li rudý posuv pozorované galaxie ze a expanzní faktor v čase vyslání světla ae, po dosazení do (1.9.1) 35 obdržíme [1.7] ∆t = t0 − te = 1 H0 a0 ae da a Ωm0 a3 0 a3 − a2 0 a2 + Ωv0 1 − a2 0 a2 + a2 0 a2 = (1.9.2) = 1 H0 ze 0 dz (1 + z) Ωm0 (1 + z)3 + Ωv0 + (1 − Ωm0 − Ωv0) (1 + z)2 . Pokud bychom započítali i hustotu záření, která klesá s a podle vztahu (1.4.4) a je popsána hustotním faktorem Ωγ, analogickými úvahami odvodíme ∆t = 1 H0 ze 0 dz (1 + z) Ωm0 (1 + z)3 + Ωγ0 (1 + z)4 + Ωv0 + (1 − Ω0) (1 + z)2 , (1.9.3) kde jsme označili Ω0 = Ωm0 + Ωγ0 + Ωv0 současnou hodnotu celkového hustotního faktoru zahrnujícího všechny formy energie [1.30]. Protože 1 + z > 0, bude Ωm0 (1 + z)3 > Ωm0 (1 + z)2 a Ωv0 < Ωv0 (1 + z)2 , větší Ωm0 vede ke zkrácení ∆t, naopak větší kladné Ωv0 má za následek větší ∆t. Kvalitativně tak můžeme posoudit vliv jednotlivých forem energie na stáří samotného vesmíru, jež je extrapolací uvedeného vztahu. Numerický model vývoje expanzního faktoru v závislosti na čase pro plochý vesmír (tj. Ω0 = 1) je znázorněn na obr. 1.12. Shrnutí Čas, za který k nám doputovalo světlo ze vzdálených objektů, stejně jako stáří samotného vesmíru závisejí na zastoupení jednotlivých forem energie. Kontrolní otázky 1. Jaký má vliv přítomnost hmoty na stáří vesmíru? 2. Jak může stáří vesmíru ovlivnit energie kosmologická konstanta? Cvičení 1. Ukažte, že pokud by vesmír obsahoval pouze energii vakua 0 < Ωv0 < 1, potom z rovnice (1.9.2) s počáteční podmínkou a(t = 0) = 0 získáme rovnici a = a0 1 − Ωv0 Ωv0 sinh Ωv0H0t . Vyšetřete limitní případy Ωv0 → 0 a Ωv0 → 1. 36 0,0 0,5 1,0 1,5 –20 –10 0 10 1 0,1 0,01 0,001 0,0001 0 0,5 1 1.5 2 3 budoucnostdnesminulost Miliardy let vzhledem k souˇcasnosti VelikostVesm´ıruvpomˇerukdneˇsn´ı relativn´ı jasnost rud´yposuv nejprve zpomaluje, pak zrychluje ...nebost´alezpomaluje Po inflaˇcn´ı f´azi Vesm´ır... zkolabuje st´ale se rozp´ın´a Obr. 1.12: Rekonstrukce a možné předpovědi rozpínání (popř. budoucího smršťování) vesmíru z výsledků měření vzdálených supernov (data znázorněná černými body) pro plochý vesmír s eukleidovskou geometrií. Současná hodnota expanzního faktoru je přeškálována na hodnotu a0 = 1, takže a = 1/(1 + z). Křivky v oblasti s modrým pozadím představují modely, u nichž převáží vliv kosmologické konstanty a rozpínání se bude zrychlovat. Křivky v oblasti se žlutým pozadím odpovídají modelům, u nichž se rozpínání bude zpomalovat, u posledních dvou modelů nakonec dojde ke zpětnému smršťování. Upraveno podle [1.33]. Řešení 1. Z rovnice (1.9.2) při Ωm0 = 0 plyne H0 dt = da a Ωv0 1 − a2 0 a2 + a2 0 a2 nebo – zavedeme-li proměnnou u = a/a0 H0 dt = du √ 1 − Ωv0 + Ωv0u2 . Vzhledem k počátečním podmínkám v zadání úlohy můžeme psát t 0 H0 dt = u 0 du √ 1 − Ωv0 + Ωv0u2 , 37 což vede k výsledku H0 t = 1 √ Ωv0 argsinh √ Ωv0 u √ 1 − Ωv0 = 1 √ Ωv0 argsinh √ Ωv0 √ 1 − Ωv0 a a0 . Obrácená závislost potom má tvar ze zadání úlohy a = a0 1 − Ωv0 Ωv0 sinh Ωv0H0t . Zřejmě pokud Ωv0 → 0, sinh Ωv0H0t ≈ Ωv0H0t a a → a0H0t, naopak pokud Ωv0 → 1, potom t ∝ argsinh √ Ωv0 √ 1 − Ωv0 a a0 → ∞. Takový model vesmíru by existoval nekonečně dlouho a nezačínal by velkým třeskem (viz levá horní část obr. 1.11). Průvodce studiem Pracujeme-li v astronomii s jednotkou světelný rok (1 ly), chápeme ji jako vzdálenost uraženou světlem ve vakuu za jeden rok. Na první pohled by mělo stačit vynásobit pravou stranu rovnice (1.9.3) rychlostí světla ve vakuu c a zjistili bychom vzdálenost objektu s příslušným rudým posuvem. Jenže přesně vzato tak jednoduché to není! Protože náš vesmír nemá geometrii příliš odlišnou od Euklidovy, pro bližší galaxie a přibližně to sice platí a někdy je dokonce takový odhad i postačující. Chceme-li však studovat Hubbleovy diagramy a zkoumat, jakým hustotním parametrům odpovídají, musíme se tímto problémem zabývat detailněji. 1.10 Fotometrická vzdálenost a Hubbleovy diagramy Hubbleovy diagramy znázorňují závislost vzdálenosti pozorovaného objektu, resp. jeho hvězdné velikosti, na rudém posuvu z a poskytují detailnější obraz vlastností vesmíru na velkých vzdálenostech. Zpočátku, ve 30-ých letech 20. století přesně odpovídaly Hubbleovu zákonu (1.5.1) s konstantním Hubbleovým parametrem H. Novější data získaná např. v rámci „Supernova cosmology project [1.20, 1.33] zahrnují pozorování mnohem vzdálenějších objektů a ukazují, že rozpínání vesmíru se skutečně s časem mění. 38 Rychlost, z jakou se od nás vzdalují daleké galaxie určujeme z rudého posuvu spektrálních čar (viz část 1.6). Zatímco rudý posuv je možné změřit s poměrně vy- Standardní svíčky jsou „kalibrované zdroje, u nichž známe (nebo předpokládáme) jejich skutečný zářivý výkon. sokou přesností, mezigalaktické vzdálenosti musíme určovat nepřímo. Nejčastěji se využívá fotometrických zákonů: skutečnosti, že intenzita osvětlení (obecně ozáření) nějaké plochy klesá se čtvercem její vzdálenosti od zdroje. Pokud bychom znali svítivost vzdáleného zdroje, potom bychom změřením světelného výkonu dopadajícího do našich měřicích přístrojů mohli určit jeho vzdálenost. K takovým měřením se používají tzv. standardní svíčky („standard candles ), tj. objekty, jejichž absolutní svítivost určujeme z jiných, nezávislých fyzikálních vlastností (viz obr. 1.13). Vlivem rozpínání vesmíru a možnému zakřivení jeho geometrie se takto určená, tzv. fotometrická vzdálenost, liší od pouhého součinu c∆t, kde ∆t je dáno vztahem (1.9.3) V praxi nejčastěji měříme fotometrickou vzdálenost. nalezeným v předcházející podkapitole 1.9. Pro šíření signálu z (1.6.1) získáváme [1.15, 1.30] t0 te cdt a(t) = re 0 dr √ 1 − kr2 = χ. (1.10.1) Analogickými úvahami jako při odvození (1.9.3) vypočítáme integrál na levé straně χ = c a0H0 ze 0 dz Ωm0 (1 + z)3 + Ωγ0 (1 + z)4 + Ωv0 + (1 − Ω0) (1 + z)2 (1.10.2) Integrál na pravé straně (1.10.1) dává různé výsledky pro různé hodnoty k, což odráží vliv geometrie vesmíru na šíření světelných signálů [1.15, 1.19] re = sin χ pro k = 1, re = χ pro k = 0, re = sinh χ pro k = −1. (1.10.3) Označíme-li L absolutní svítivost zdroje (celkový vyzářený výkon) a F jeho zdánlivou svítivost (zářivý tok, tj. výkon záření dopadající v místě pozorovatele na jednotkovou plochou kolmou na směr šíření záření), potom v plochém prostoru pro fotometrickou vzdálenost dL platí [1.15, 1.17] Zavedení fotometrické vzdálenosti odráží pokles zářivého toku se čtvercem vzdálenosti od zdroje. F = L 4pd2 L (1.10.4) Nachází-li se zdroj v místě s „comoving souřadnicí re a vyšle signál v čase te, dopadne světlo v současnosti (tj. v čase t0) na plochu S = 4pa2 0r2 e ; re je přitom obecně dáno vztahy (1.10.3). V důsledku kosmologického rudého posuvu (1.6.4) však bude energie každého fotonu namísto hc/λe rovna hc/λ0 = = hc/[λe (1 + z)]. Díky poklesu frekvence záření fotony vyzářené v intervalu ∆te dopadnou do místa pozorování v intervalu ∆t0 = ∆te (1 + z) 39 Obr. 1.13: Supernova 1994Dna okraji galaxie NGC 4526 je příkladem supernovy typu Ia, jejichž kulminační jasnost převýší svítivost celého jádra galaxie; proto lze supernovy pozorovat i ve velkých vzdálenostech. V tomto případě se z kosmologického hlediska díváme „za humna přibližně 20 Mpc. Supernovy typu Ia jsou příkladem standardních svíček, u nichž byla empiricky zjištěna závislost mezi kulminační svítivostí a dobou poklesu jasnosti [1.33]. Pro zářivý tok v místě pozorování pak vychází F = L 4pa2 0r2 e (1 + z)2 a srovnáním s (1.10.4) pro fotometrickou vzdálenost odvodíme Fotometrická vzdálenost závisí na současné velikosti vesmíru a rudém posuvu pozorované galaxie. dL = rea0 (1 + z) , (1.10.5) kde za re opět musíme dosadit z (1.10.3) a (1.10.2). Pro malé rudé posuvy přibližně do hodnoty z ≈ 0,4 vystačíme s užitečnou aproximací. Rozvineme-li expanzní faktor a(t) v Taylorovu řadu v okolí současné hodnoty 40 a(t0) a(t) ≈ a(t0) + da dt 0 (t − t0) + d2 a dt2 0 (t − t0)2 + . . . a dosadíme podle (1.5.2) a (1.7.7) a(t) ≈ a(t0) 1 + H0 (t − t0) − 1 2 q0H2 0 (t − t0)2 + . . . S použitím (1.10.5) pak dospějeme ke vztahům [1.17] dL = c H0 z + 1 2 (1 − q0) z2 + . . . respektive z = H0 c dL + 1 2 (1 + q0) H0 c d2 L + . . . . U členů vyšších řádů se již projevuje zakřivení geometrie a situace je mnohem komplikovanější [1.17]. Pro malé z ≈ v/c pak v 1. přiblížení dostáváme Hubbleův zákon (1.5.1). Dodejme, že Edwin Hubble při svých měřeních z roku 1929 pozoroval galaxie se z ≈ 0,003. V praktické astronomii se častěji užívají historicky zavedené magnitudy. Při samotných pozorováních pracujeme obvykle namísto fotometrické vzdálenosti s historicky zavedenými hvězdnými velikostmi neboli magnitudami, i když v případě kosmologie nejde o pozorování hvězd, nýbrž galaxií. Magnitudy byly zavedeny již ve starověku Ptolemaiem a Hipparchem a odrážejí zkušenost, že naše fyziologické vjemy jsou úměrné logaritmu intenzity podnětu (dnes se však už pouze na lidské smysly zdaleka nespoléháme). Absolutní hvězdná velikost (magnituda) M je definována vztahem M = −2,5 log10 L L + 4,74, (1.10.6) kde L = 3,85·1026 W je zářivý výkon Slunce, relativní hvězdná velikost (magnituda) pak m = −2,5 log10 F F v 10 pc + 4,74, kde F v 10 pc = 3,21·10−10 W·m−2 je zářivý výkon Slunce, který by dopadal na jednotku plochy, pokud by se nacházelo ve vzdálenosti 10 pc. Obecně je proto zdánlivá magnituda objektu rovna jeho absolutní magnitudě, kterou by měl ve vzdálenosti 10 pc od nás. Rozdíl Fotometrickou vzdálenost pak nahrazujeme modulem vzdálenosti. m − M = 5 log10 dL 10 pc nazývaný modulem vzdálenosti je astronomickým vyjádřením vztahu (1.10.4). Příklad závislosti magnitudy na rudém posuvu pro vzdálené supernovy je na obr. 1.14. 41 20 21 22 23 24 25 0,01 0,02 0,04 0,1 14 16 18 20 22 24 26 0,40,2 0,6 1,0 0,40,2 0,6 1,0 magnituda rudý posuv Type Ia Supernovae Calan/Tololo Supernova Survey High-Z Supernova Search Supernova Cosmology Project slabźí Zpomalující se rozpínání Zrychlující se rozpínání prázdný hustota hmoty 0 1 0,1 1 0,01 0,001 0,0001 Relativníjasnost 0,70,8 0,6 0,5 Velikost Vesmíru [vzhledem k dnešní] Obr. 1.14: Příklad Hubbleova diagramu, tj. závislosti pozorované magnitudy na rudém posuvu z, odpovídající datům získaným pozorováním supernov typu Ia. Pro z = 0,1 (odpovídá vzdálenosti asi 109 ly se začínají křivky lišit v závislosti na předpokládaných hodnotách hustotních parametrů Ωm0 a Ωv0. Červené křivky odpovídají modelům bez kosmologické konstanty Ωv0 = 0, a hodnotám Ωm0 = 1 až Ωm0 = 0 (prázdný vesmír). Naměřeným datům nejlépe odpovídá modrá křivka s hodnotami Ωm0 ≈ 0,3 a Ωv0 ≈ 2Ωm0, což podle rovnice (1.7.8) znamená, že se rozpínání vesmíru v současné době zrychluje. Upraveno podle [1.33]. Shrnutí Vlastnosti vesmíru na velkých vzdálenostech dobře charakterizují Hubbleovy diagramy, tj. závislosti vzdálenosti galaxií na rudém posuvu. Měření vzdáleností představuje v kosmologii velký problém, většinou pracujeme s tzv. fotometrickou vzdáleností, při samotných pozorováních i s hvězdnou velikostí neboli magnitudou. Také z Hubbleových diagramů pro supernovy ve vzdálených galaxiích vyplývá, že rozpínání vesmíru se v současné době zrychluje. Pojmy k zapamatování • Hubbleův diagram • standardní svíčky • fotometrická vzdálenost • hvězdná velikost (magnituda) • modul vzdálenosti 42 Kontrolní otázky 1. Jakou závislost představují Hubbleovy diagramy? 2. Jak zavádíme absolutní a relativní magnitudu? Cvičení 1. Ve spektru kvasaru 3C 273 je emisní čára vodíku Hβ o laboratorní vlnové délce 486,1 nm posunuta o 77,8 nm směrem k dlouhovlnnému konci spektra. Určete a) vzdálenost kvasaru, b) lineární rozměru kvasaru, jestliže jeho úhlový průměr činí 0,24 , c) velikost výtrysku z kvasaru, jehož úhlová velikost je 19,5 , d) zářivý výkon kvasaru, jestliže absolutní bolometrická hvězdná velikost je −25 mag [1.38]. 2. Zářivý výkon kvasarů dosahuje 1040 W. Fyzikální podstata procesů umožnujících uvolňování tak obrovské energie není dosud definitivně objasněna. Vypočtěte, jaké množství hmoty v jednotkách M za rok by se muselo přeměnit, aby pokrylo tento výkon a) při termonukleární reakci s účinností η = 0,01, b) při akreci na černou díru s účinností η = 0,2 [1.38]. 3. Šířka čáry Hβ ve spektru jádra seyfertovské galaxie je zhruba 3 nm. jaké jsou charakteristické rychlosti pohybu mračen plynu v jádře takové galaxie [1.38]? 4. Rádiový zdroj v jádře aktivní galaxie má úhlovou velikost 0,001 , kosmologický rudý posuv je z = 0,5. Odhadněte lineární rozměry zdroje v pc [1.38]. Uvažujte Hubbleovu konstantu H0 = 71 km·s−1 ·M−1 pc. Řešení 1. a) Ze zadání λ = 486,1 nm, ∆λ = 77,8 nm, z = ∆λ/λ ≈ 0,16. Můžeme proto použít přibližný vztah z ≈ v/c, neboli v = cz. Z Hubbleova zákona (1.5.1) v ≈ H0d a konečně d = v H0 = cz H0 = ∆λc λH0 , což pro H0 ≈ 71 km·s−1 ·Mpc −1 dává d ≈ 674 Mpc. b) Pro úhel α = 0,24 ≈ 1,16·10−6 rad vychází lineární rozměr lα = dα ≈ 782 pc. c) Analogicky pro úhel β = 19,5 ≈ 9,45·10−5 rad vyjde lβ = dβ ≈ 63,7 kpc. d) Z rovnice (1.10.6), kde položíme M = −25 odvodíme L = 10(4,74−M)/2,5 L = 101,9−0,4M L po dosazení zářivého výkonu Slunce L = 3,85·1026 W vychází zářivý výkon kvasaru L = 3,03·1038 W. 43 2. Vyjdeme z Einsteinova vztahu mezi hmotností a energií E = Lt = η∆mc2 , odkud dostáváme a) ∆m1 = Lt ηc2 = 3,53·1032 kg ≈ 177 M , b) ∆m2 = Lt ηc2 = 1,76·1031 kg ≈ 9 M . 3. Ze zadání plyne 2∆λ = 3 nm, vlnová délka čáry Hβ = 486,1 nm. Koeficient 2 u ∆λ vyjadřuje skutečnost, že plyn se pohybuje jak směrem k nám, tak od nás, pohyb v jednom směru tak statisticky zodpovídá za polovinu rozmazání spektrální čáry. Podle vztahu (1.6.3) dále platí z = ∆λ λ ≈ 0,0062 Ć 1. Můžeme proto opět použít přibližný vztah z ≈ v/c a pro c ≈ 3·108 m·s−1 vychází v = c z = c ∆λ λ ≈ 9,23·105 m·s−1 . 4. Vzhledem k velké hodnotě rudého posuvu musíme použít relativistický vztah 1 + z = 1 + v c 1 − v c . V kombinaci s Hubbleovým zákonem (1.5.1) v ≈ H0r a vlastnosti úhlů v radiánech d ≈ αr (přesná rovnost by platila v plochém prostoru) vychází d ≈ α v H0 = αc H0 (1 + z)2 − 1 (1 + z)2 + 1 . Pro zadané hodnoty α = 0,001 = 4,85·10−9 rad, c = 2,99·105 km·s−1 , z = 0,5 a H0 = 71 km·s−1 ·M−1 pc pod dosazení obdržíme d ≈ 7,86·10−6 Mpc = 7,86 pc. 1.11 Závěr V této kapitole jsme nastínili základní představy moderní relativistické kosmologie, základní scénáře možného vývoje vesmíru i výsledky nejnovějších měření jeho základních observačních parametrů. V následujících letech můžeme očekávat další, přesnější a důkladnější pozorování, která bezpochyby naše porozumění vesmíru prohloubí a přinesou i nejedno překvapení. Čtenáře s hlubším zájmem o tuto problematiku můžeme odkázat na mnoho zajímavých monografií: od populárních [1.2, 1.3, 1.13, 1.14, 1.34, 1.39] až po odborné [1.17, 1.27, 1.29, 1.30, 1.32, 1.38]. Řečeno slovy amerického nositele Nobelovy ceny S. Weinberga: „úsilí pochopit vesmír je jednou z velmi mála věcí, která zvedá lidský život trochu nad úroveň frašky a dává mu něco z krásy tragédie [1.39]. 44 Literatura ke kapitole 1 [1.1] Allen S.W., Schmidt R.W., Fabian A.C.: „Cosmological constraints from the Xray gas mass fraction in relaxed lensing clusters observed with Chandra . Mon. Not. Roy. Astron. Soc. 334 (2002), L11. ArXiv: http://xxx.lanl.gov/abs/ astro-ph/0205007. [1.2] Barrow J.D.: Teorie všeho. Mladá fronta, Praha 1999. [1.3] Barrow J.D.: Teorie ničeho. Mladá fronta, Praha 2005. [1.4] Bennett C.L., Banday A.J., Gorski K.M. a kol.: „Four-Year COBE DMR Cosmic Microwave Background Observations: Maps and Basic Results . Astrophys. J. 464 (1996), L1–L4. URL: http://adsabs.harvard.edu/cgi-bin/nph-bib_ query?bibcode=1996ApJ...464L...1B&db_key=AST. [1.5] Bennett C.L. a kol.: „First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Preliminary Maps and Basic Results . Astrophys. J. Suppl. 148 (2003), 1–27. ArXiv: http://xxx.lanl.gov/abs/astro-ph/ 0302207. [1.6] Carroll S.M.: Spacetime and Geometry: An Introduction to General Relativity. Addison-Wesley 2003. ArXiv: http://xxx.lanl.gov/abs/gr-qc/9712019. [1.7] Carroll S.M., Press W.H., Turner E.L.: „The cosmological constant . Annual Rev. Astron. Astrophys. 30 (1992), 499–542. URL: http://adsabs.harvard.edu/ cgi-bin/nph-bib_query?bibcode=1992ARA%26A..30..499C&db_key=AST. [1.8] Colless M., Boyle B.: „Redshift Surveys with 2dF . (1997). ArXiv: http://xxx. lanl.gov/abs/astro-ph/9710268. [1.9] Colless M. a kol. (The 2DFGRS): „The 2dF Galaxy Redshift Survey: Spectra and redshifts . Mon. Not. Roy. Astron. Soc. 328 (2001), 1039. ArXiv: http:// xxx.lanl.gov/abs/astro-ph/0106498 URL: http://www.mso.anu.edu.au/ 2dFGRS/. [1.10] Čulík F., Noga M.: Úvod do štatistickej fyziky a termodynamiky. Alfa, Bratislava 1982. [1.11] Felten J.E., Isaacman R.: „Scale factors R(t) and critical values of the cosmological constant Λ in Friedman universes . Rev. Mod. Phys. 58(3) (1986), 689–698. [1.12] Freedman W.L. a kol.: „Final Results from the Hubble Space Telescope Key Project to Measure the Hubble Constant . Astrophys. J. 553 (2001), 47–72. ArXiv: http://xxx.lanl.gov/abs/astro-ph/0012376. [1.13] Gygar J.: Vesmírná zastavení. Pyramida, Praha 1990. [1.14] Gygar J.: Vesmír jaký je. Mladá fronta, Praha 1997. [1.15] Hartle J.B.: Gravity: An Introduction to Einstein’s Relativity. Addison Wesley, San Francisco 2003. [1.16] Horský J., Bartoň S.: Relativistický vesmír. Ando Publishing, Brno 1997. [1.17] Horský J., Novotný J., Štefaník M.: Úvod do fyzikální kosmologie. Academia, Praha 2004. [1.18] Janků V.: Základy statistické fyziky. UP, Olomouc 1983. [1.19] Jordan T.F.: „Cosmology calculation almost without general realtivity . Am. J. Phys. 73(7) (2005), 653–662. ArXiv: http://xxx.lanl.gov/abs/astro-ph/ 0309756. 45 [1.20] Knop R.A. a kol. (The Supernova Cosmology Project): „New Constraints on ΩM , ΩΛ, and w from an Independent Set of Eleven High-Redshift Supernovae Observed with HST . Astrophys. J. 598 (2003), 102–137. ArXiv: http://xxx. lanl.gov/abs/astro-ph/0309368. [1.21] Kvasnica J.: Statistická fyzika. Academia, Praha 1983 a 1998. [1.22] Landau, L.D., Lifxic, E.M.: Statistiqeska fizika Qast 1. Nauka, Moskva 1976. [1.23] Liddle A.R.: An Introduction to modern Cosmology. John Wiley & Sons Ltd, Chichester 1999. [1.24] Maddox S.J., Efstathiou G., Sutherland W.J.: „The APM Galaxy Survey - Part Two - Photometric Corrections . Mon. Not. Roy. Astron. Soc. 246 (1990), 433– 457. URL: http://adsabs.harvard.edu/cgi-bin/nph-bib_query?bibcode= 1990MNRAS.246..433M&db_key=AST. [1.25] Misner C.W., Thorne K.S., Wheeler J.A.: Gravitation. W. H. Freeman &, San Francisco 1973. [1.26] Møller C.: The Theory of Relativity. Clarendon Press, Oxford 1962. [1.27] Narlikar J.V.: Introduction to Cosmology. Cambridge University Press, Cambridge 1993. [1.28] Narlikar J.V.: „Spectral shifts in general relativity . Am. J. Phys. 62(10) (1994), 903–907. [1.29] Padmanabhan T.: Cosmology and Astrophysics Through Problems. Cambridge University Press, Cambridge 1996. [1.30] Peacock J.A.: Cosmological Physics. Cambridge University Press, Cambridge 1999. [1.31] Peebles P.J.E.: „Tests of cosmological models constrained by inflation . Astrophys. J. 284 (1984), 439–444. URL: http://adsabs.harvard.edu/cgi-bin/ nph-bib_query?bibcode=1984ApJ...284..439P&db_key=AST. [1.32] Peebles P.J.E.: Principles of physical cosmology. Princeton University Press, New Jersey 1993. [1.33] Perlmutter S.: „Supernovae, Dark Energy, and the Accelerating Universe . Physics Today 56(4) (2003), 53–60. URL: http://panisse.lbl.gov/. [1.34] Rees M.: Náš neobyčejný vesmír. Dokořán, Praha 2002. [1.35] Rektorys K. a kol.: Přehled užité matematiky I, II. SNTL, Praha 1988. [1.36] Schutz B.: Gravity from the ground up. Cambridge University Press, Cambridge 2003. URL: http://www.gravityfromthegroundup.org/. [1.37] Spergel D.N. a kol. (WMAP): „First Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Determination of Cosmological Parameters . Astrophys. J. Suppl. 148 (2003), 175. ArXiv: http://xxx.lanl.gov/abs/ astro-ph/0302209. [1.38] Štefl V., Krtička J.: Úlohy z astrofyziky. PřF MU, Brno 2000. [1.39] Weinberg S.: První tři minuty. Mladá fronta, Praha 1998. [1.40] Wesson P.S.: „Olbers’s paradox and the spectral intensity of the extragalactic background light . Astrophys. J. 367 (1991), 399–406. URL: http://adsabs.harvard.edu/cgi-bin/nph-bib_query?bibcode=1991ApJ. ..367..399W&db_key=AST. 46 Kapitola 2 Základy kvantové informace Studijní cíle: Dozvíme se proč se lidé začali zabývat otázkami kvantové informace a komunikace, co to jsou qubity a kvantová hradla, jaký je princip kvantového počítače a k čemu by mohl být dobrý, jak funguje kvantová kryptografie či kvantová teleportace. Klíčová slova: Kvantový bit, qubit, kvantové hradlo, kvantový počítač, kvantově provázané stavy. Potřebný čas: 450 minut. 2.1 Úvod a něco z historie 2.1.1 Moorův zákon Naše výpočetní možnosti neustále rostou. Kvantitativně to vystihuje tzv. Moorův zákon: množství tranzistorů které dokážeme umístit na jeden čip se zdvojnásobí zhruba každý rok a půl. Toto pravidlo v ještě optimističtější podobě - zdvojnásobení každý rok - formuloval Gordon Moore, spoluzakladatel firmy Intel roku 1965. Napsal tehdy, že tento trend bylo možné pozorovat několik dosavadních roků a předpověděl, že ještě pár let by to tak mělo jít dál. Za deset let bychom se prý měli dostat z tehdejších šedesáti tranzistorů na neuvěřitelných šedesát tisíc. Moore svůj odhad později zmírnil na zdvojnásobení každé dva roky, ale sám neočekával platnost vyřčeného pravidla po delší dobu než nějakých deset či patnáct let. Uplynulo čtyřicet let a jak si můžeme ověřit, exponenciální nárůst počtu tranzistorů se nezastavil. Například roku 1989 obsahoval mikroprocesor Intel486 1,2 milionů tranzistorů, roku 2004 udává Intel pro svůj Itanium mikroprocesor 590 milionů tranzistorů. To je téměř devět zdvojnásobení během patnácti let, čili zdvojnásobení každý rok a osm měsíců. Jednoho dne se ale rostoucí exponenciála musí zastavit. Velikost atomu je jistě principiální překážka - tranzistor může být sotva menší než jeden atom. Nejspíše ale ještě dříve narazíme na jiné těžkosti - snad to bude naše schopnost odvádět odpadní teplo, možná to budou ekonomické bariéry. Kdy se dnešní exploze zastaví? Bude trvat ještě deset nebo patnáct let? Ať už ji zastaví cokoliv, ještě několit roků bude vývoj v různých oblastech fyziky pomáhat v udržení platnosti Moorova zákona. Ale na věky to zřejmě nepůjde. 47 2.1.2 Limity výpočetních možností Kvantově mechanické modely I když je růst možností výpočetní techniky impozantní, existují problémy, se kterými si naše počítače sotva kdy poradí. Jedním z nich je numerické modelování větších kvantově mechanických systémů. Problémem je již pouhé zadání stavu do paměti počítače. Zatímco v klasické fyzice je stav popsán zadáním polohy a hybnosti (tedy dvou vektorů) pro každou částici, v kvantové mechanice se stav popisuje pomocí vlnové funkce. Tato funkce má tolik argumentů, kolik stupnů volnosti má celý systém. Představme si jednu částici, pohybující se v jednom rozměru. Klasická fyzika potřebuje na popis jejího stavu dvě reálná čísla (polohu a hybnost), kvantová fyzika musí zadat komplexní číslo pro každý bod. Rekněme, že pro numerickou simulaci si vystačíme s rozdělením úsečky, podél které se částice pohybuje, na 100 dílků. Ve sto bodech musíme zadat hodnotu vlnové funkce - tedy její reálnou a imaginární část. Předpokládejme, že pro každou z nich nám stačí přesnost daná pamětí jednoho bytu - čili 8 bitů, tedy přesnost omezená na zhruba tři platné číslice. Pro jednu dimenzi tedy potřebujeme zhruba 200 bytů. Buďme velkorysí a řekněme, že se spokojíme jen se 100 byty. Pokud uvažujeme dvourozměrný pohyb částice, musíme zadat její vlnovou funkci na síti 100×100 bodů. Při požadované přesnosti to znamená 10 tisíc bytů, což je stále hravě zvládnutelné. Pro třírozměrný pohyb je to 1003 = 106 , tedy milion bytů. Stále je to řešitelné i na běžném notebooku. Uvažujme ale nyní dvě částice v třírozměrném prostoru. Konfigurační prostor tohoto systému je šestirozměrný a na jeho popis už potřebujeme 1006 = 1012 bytů, tedy zhruba terabyte. Naše počítače už dostávají zabrat. Zajímá nás kvantověmechanický problém tří interagujících částic ve třech rozměrech (třeba elektrony v atomu lithia)? Pro numerickou simulaci potřebujeme 1009 = 1018 bytů a to je již nad naše možnosti. Podobně se do problémů dostaneme i při studování mnohem jednodušších systémů - například spinů elektronů. Spin elektronu (jeho vnitřní moment hybnosti) může nabývat jen dvou projekcí do každého směru. Rekněme, že měříme průmět spinu do osy z: můžeme dostat hodnotu buď +¯h nebo −¯h. Každá z těchto dvou hodnot odpovídá možnému stavu spinu. Tyto stavy si můžeme označit jako | ↑ a | ↓ . Kromě toho se však podle zákonů kvantové mechaniky každý spin může nacházet v libovolné superpozici těchto dvou stavů, tedy ve stavu |ψ = α| ↑ + β| ↓ , kde α a β jsou dvě komplexní čísla splnující |α|2 + |β2 | = 1 (o významu těchto čísel a jejich vlivu na výsledky měření se dozvíme více v následujících kapitolách). Zadání těchto dvou čísel definuje stav jednoho spinu. Pokud budeme mít za úkol popsat stav dvou spinů, musíme zadat čtyři koeficienty v superpozici |ψ = α1| ↑↑ + α2| ↑↓ + α3| ↓↑ + α4| ↓↓ . (2.1.1) Při celkovém počtu N spinů to představuje 2N koeficientů. Možnost numerického modelování interagujících spinů (což je důležité například pro studium magnetických vlastností látek) tak velice rychle klesá s počtem uvažovaných spinů. 48 Průvodce studiem Trochu to připomíná pověst o vynálezci šachu. Když se ho král zeptal, jakou odměnu si žádá za autorství tak úžasné hry, vynálezce si prý řekl o pár zrnek obilí. Na první políčko šachovnice nechť král položí jedno zrnko, na druhé dvě, na třetí čtyři a na každé další dvojnásobný počet toho, co bylo na předchozím políčku. Na celou šachovnici by to představovalo 264 − 1 zrnek. Ačkoliv se na první pohled mohlo zdát, že se jedná o směšně malou odměnu, mohl ji král stěží splnit. Výsledné množství řádově tisíckrát přesahuje současnou roční světovou produkci pšenice. K podobné beznaději přijdeme při snahách modelovat kvantově mechanický vývoj systémů složených z více částic. Sice můžeme dostat řadu užitečných výsledků použitím různých aproximací a zjednodušujících předpokladů, na sledování libovolného stavu ale naše možnosti nestačí. „Složité problémy Kromě modelování kvantově mechanických systémů zná matematika řadu úloh, které se s rostoucím vstupem stávají stále hůře řešitelné tak, že velmi rychle přesáhnou možnosti i těch nejlepších počítačů. Typickým příkladem je problém obchodního cestujícícho: máme zadáno N měst, která se mají navštívit. V jakém pořadí městy projedeme, aby byla uražená cesta nejkratší? Další příklad by se dal popsat jako úkol ubytovat N studentů ve čtyřlůžkových pokojích na kolejích s tím, že každý student uvede několik požadavků s kým nesnese, aby byl na pokoji. Důležitou úlohou, jejíž složitosti vděčíme za bezpečný přenos dat, je faktorizace velkých čísel. Pokud budete mít za úkol rozložit na prvočinitele číslo 119, budete za chvilku hotovi, ale u čísla 16637 už to půjde déle. S rostoucím počtem číslic se tento problém rychle stane neřešitelným i pro naše počítače. Kódování pomocí velkých čísel, kde klíč představuje znalost jejich prvočinitelů, je základem některých dnešních kryptografických systémů (o principu tzv. RSA kódu se můžeme dočíst v příloze B). 2.1.3 Ať počítá kvantový systém! Kvantové superpozice Skepse nad nemožností využít naše počítače k simulaci kvantových systémů se ale dá obrátit. V osmdesátých letech s tímto nápadem přišel Richard Feynman (nositel Nobelovy ceny za fyziku z roku 1965 za kvantovou elektrodynamiku a autor vynikajících „Feynmanových přednášek z fyziky [2.1]). Každý počítač je konec konců fyzikální systém a počítání je fyzikální proces. Všechny dosavadní počítače však k reprezentaci číslic a jejich zpracování využívají pouze zákony klasické fyziky (to, že děje uvnitř tranzistoru atd. jsou ve své podstatě kvantové, tu nehraje roli). Pokud 49 bychom ale využili všech možností, které nám dává kvantová mechanika, vypadalo by počítání podstatně jinak. V „klasickém počítači (tedy v tom, který pracuje na základě klasické fyziky) může být každý bit ve stavu 0 nebo 1. Soubor N bitů odpovídá 2N možnostem, v jakých stavech tento soubor může být. Pokud bychom však využili kvantové mechaniky, víme, že podle principu superpozice se každý systém může nacházet v libovolné superpozici jakýchkoliv jiných přípustných stavů. S N kvantovými bity můžeme tedy pracovat s mnohem více než s 2N stavy—máme k dispozici i jakoukoliv jejich super- pozici. Průvodce studiem Princip superpozice je jedním ze základních principů kvantové fyziky: pokud může být systém v jednom stavu nebo ve druhém, může být i v jakékoliv superpozici těchto stavů. Erwin Schrödinger, jeden ze zakladatelů kvantové fyziky dospěl k názoru, že z tohoto principu vyplývají důsledky, které se příčí zdavému rozumu—a už málem chtěl kvantovou fyziku zavrhnout. Například radioaktivní jádro může dojít po určité době do superpozice stavů „nerozpadlé jádro a „nové jádro a alfa částice . Jenže alfa částice může spustit řetězec reakcí v Geiger-Müllerově počítači, proudový impuls pak může spustit pekelný stroj, který rozbije ampuli s jedem a jed usmrtí kočku. Pokud může být atom v superpozici stavů (tj. zároveň v jednom i druhém stavu), může být v takové superpozici i sama kočka—živá a mrtvá zároveň—a to je podezřelé! Uvažujme například pouhé dva bity. Jejich hodnoty si můžeme označit jako (0,0), (0,1), (1,0) a (1,1), případně v desítkové soustavě jako 0, 1, 2 a 3. V kvantové fyzice by těmto hodnotám odpovídaly například orientace dvojic spinů, které odpovídají stavům označeným jako |0 , |1 , |2 a |3 . Kromě těchto stavů však kvantově mechanický systém může být v mnoha jiných stavech, například |ψ1 = 1 √ 2 |0 − 1 √ 2 |2 , |ψ2 = 1 2 |1 + i √ 3 2 |2 , |ψ3 = 1 2 (0 + |1 + i|2 − |3 ) , |ψ4 = 1 2 (0 − |1 + |2 − |3 ) , . . . Snadno si představíme, že možností, do jakého stavu náš systém umístit, s rostoucím počtem bitů rychle přibývá. Můžeme toho ale nějak využít pro naše počítání? Deutschova úloha Nebylo příliš jednoduché najít konkrétní příklad, kdy by vývoj kvantového systému mohl představovat výpočetní algoritmus, který by byl efektivnější, než stan- 50 dardní algoritmy. Jedním z nich je tzv. Deutschova úloha. Předpokládejme, že kolega má v počítači zadanou funkci proměnné x ∈ {0,1}, která může nabývat hodnot 0 nebo 1. Celkem existují čtyři takové funkce: f1(x) = 0, f2(x) = 1, f3(x) = x a f4(x) = 1−x. My zatím nevíme, o kterou funkci se jedná, můžeme ale do ní dosazovat a ptát se na výsledek. Řekněme, že máme odpovědět ano či ne na otázku, zda je zadaná funkce konstantní. Pokud bychom měli takovouto funkci naprogramovanou v běžném počítači, museli bychom do ní dosazovat celkem dvakrát. Když však budeme kódovat informaci pomocí kvantových stavů a využijeme principu superpozice, stačí pouze jediné dosazení. Shorova faktorizace Ačkoliv po Deutschově úloze přišlo ještě několik dalších jednoduchých schémat, kde byl kvantový přístup efektivnější než s klasickými algoritmy, nějakou dobu to byly spíše jen teoretické hříčky, u nichž nelze očekávat nějakou praktickou aplikaci. Přelom nastal roku 1994, kdy Peter Shor z telekomunikačních laboratoří AT&T publikoval algoritmus, který by na kvantovém počítači dokázal efektivně faktorizovat velká čísla (internetově přístupná verze viz [2.4]). Realizace takového algoritmu by měla zásadní dopad na bezpečnost přenosu kódovaných informací, protože běžně užívané kódy jsou založeny na složitosti faktorizace velkých čísel. Jak takový algoritmus funguje se dozvíme v kapitolách 2.4 a 2.5. Groverovo vyhledávání Jiným příkladem efektivity kvantového zpracování informace je Groverův vyhledávací algoritmus. Ve svém článku z roku 1997 nazvaném „Kvantová mechanika pomáhá hledat jehlu v kupce sena [2.6] přišel L. K. Grover z Bellových laboratoří s efektivním postupem vyhledávání v rozsáhlých databázích. Představme si, že neznámá kráska nám zanechala své telefonní číslo a my si nemůžeme vzpomenout na její jméno. Můžeme vzít telefonní seznam a probírat jedno číslo po druhém. Kolik záznamů musíme přečíst, než padneme na ten správný? Sice můžeme mít štěstí a být úspěšní hned napoprvé, můžeme však být smolaři a z N záznamů projít všechny. V průměru však budeme muset číst N/2-krát. Groverův kvantový algoritmus umožňuje, aby se počat čtení omezil na řádově √ N. Tedy pokud má seznam milión položek, nemusíme jich přečíst 500 000, ale stačí jen zhruba tisíc. O principu Groverova algoritmu se dozvíme v kapitole 2.6. Modelování kvantových systémů Zajímavým výsledkem je i nalezený způsob, jak kvantovým počítačem modelovat evoluci kvantových systémů. Přesněji řečeno, jakým způsobem z elementárních kvantových „procesorů sestavit předem zadaný evoluční operátor popisující zkoumaný kvantový systém. 51 2.1.4 Problémy konstrukce kvantového počítače Křehkost kvantových bitů Jakkoliv nalezené kvantové algoritmy vypadají slibně, zásadní problém spočívá v tom, jak takový počítač fyzicky sestrojit. Kvantové bity (zvané qubity) jsou buď velmi citlivé na rušivé vlivy okolí, anebo je značně obtížné je přinutit, aby spolu navzájem „mluvily . První problém lze ilustrovat na tomto příkladu: bit může být realizován třeba jako jeden ze dvou spinových stavů nějakého atomu drženého v pasti (může jít například o iont udržovaný na svém místě vhodným elektrickým polem, nebo neutrální atom zachycený v tzv. optické mřížce, stojaté vlně tvořené interferujícími laserovými svazky). Řekněme, že pokud spin valenčního elektronu míří dolů, znamená to logickou nulu, zatímco spin mířící vzhůru znamená logickou jedničku. Potřebujeme připravit a udržet qubit ve stavu 1√ 2 (|0 +|1 ). Elektron však mění svůj spin pod vlivem vnějšího magnetického pole—a to může být způseobeno třeba okolními atomy. Bez naší kontroly tak může dojít třeba k tomu, že se nula překlopí na jedničku a naopak. Pokud takovémuto překlápění nějak zabráníme, stále ještě snadno může docházet k fázovým posuvům: vstupní superpozice se změní na 1√ 2 (|0 + eiϕ |1 ), kde ϕ je nějaká neznámá fáze, daná konkrétní interakcí atomu s okolím. Výsledný stav qubitu je tedy odlišný od vstupního a na výstup z počítače se pak v žádném případě nemůžeme spolehnout. Musíme tedy hledat způsob jak systém qubitů co nejlépe izolovat od okolí, případně jak korigovat drobné chyby, které se v průběhu času vyskytly. Oba přístupy jsou dnes v centru intenzivního vědeckého výzkumu. Průvodce studiem V křehkosti kvantových superpozic spočívá jedna z odpovědí, proč nepozorujeme „Schrödingerovy kočky , tedy superpozice různých stavů (živý a zároveň mrtvý) u makroskopických objektů. Vlivem i nepatrné interakce s okolím (například místnost s kočkou opustí jediný foton) se systém dostává do jednoho z mnoha „normálních stavů. Uchování neporušených superpozic mnohaqubitových stavů (mnoho qubitů je skoro makroskopický objekt) je pak jednou z hlavních výzev pro realizaci kvantového počítače. Nesnadná interakce mezi qubity Některé systémy kvantových bitů jsou naopak poměrně robustní. Ač se to nezdá, patří k nim třeba světlo. Qubit můžeme realizovat pomocí polarizace jediného fotonu. Pokud bude elektrické pole horizontálně letícího fotonu kmitat vodorovně (↔), bude to znamenat logickou nulu, svislé kmitání ( ) bude znamenat jedničku. V průhledném prostředí se takový foton může šířit na veliké vzdálenosti aniž by se jeho polarizační stav nějak narušil. Problém však je přinutit fotony interagovat: jak dosáhnout toho, aby polarizační stav jednoho fotonu ovlivnil polarizaci druhého fotonu? Něco takového je ale nezbytně nutné, pokud bychom chtěli postavit „fotonový 52 kvantový počítač s polarizačními qubity. I zde je spousta zajímavých otázek, které čekají na své vyřešení. 2.1.5 Jiné aplikace kvantové informatiky Kvantová kryptografie Ačkoliv by kvantové počítače mohly představovat určitou hrozbu pro naše soukromí (rychlá faktorizace velkých čísel by znamenala snadné rozluštění kódovaných zpráv), přináší jiné odvětví kvantové informatiky naopak možnost jak dosáhnout bezpečného přenosu utajovaných dat. Využívá se tu principiální nemožnosti spolehlivě rozlišovat mezi neortogonálními kvantovými stavy a toho, že měření kvantový stav obecně narušuje. Základem kvantové kryptografie je distribuce tajného klíče mezi dvěma vzdálenými partnery tak, aby nikdo jiný klíč nemohl znát. Jakmile dvě strany klíč znají (je to v podstatě náhodná posloupnost nul a jedniček), stačí jej přičíst k binárně kódované zprávě. Takto zašifrovaný vzkaz se pak vnějšímu pozorovateli jeví jako zcela chaotická změť nul a jedniček bez jakýchkoliv vnitřních korelací, které by se daly využít k získání nějaké užitečné informace. Tajný klíč se dá předat například posláním řady qubitů. Jak uvidíme později, je přitom zapotřebí náhodně měnit bázi, v jaké se nuly a jedničky kódují do kvantových stavů. Pokud by nějaký špión chtěl předávaný klíč odposlechnout, musel by přitom do systému vnést takový šum, že by byl snadno odhalitelný. Kvantová kryptografie je z oblastí kvantové informatiky pravděpodobně nejblíže komerčnímu využití. Jako ukázku jejích možností provedl vědecký tým prof. Antona Zeilingera z vídeňské univerzity první skutečnou bankovní transakci na jaře roku 2004: slabounkým proudem fotonů proběhla platba mezi vídeňskou radnicí a pobočkou místní banky. O kvantové kryptografii bude pojednávat kapitola 2.7. Kvantová teleportace Ve sci-fi příbězích funguje teleportace tak, že cestující je ve startovní stanici rozebrán na své nejmenší součástky - řekněme atomy, informace o těchto atomech je poslána světelným či jiným elektromagnetickým zářením do cílové stanice, kde je pak cestující z odpovídajících součástek opět poskládán dohromady. Kvantová teleportace znamená podobnou přepravu kvantových stavů. Systém ve startovní stanici se nachází v nějakém neznámém kvantovém stavu. Pokud bychom tento stav změřili a zjištěnou informaci poslali do cílové stanice, mohli bychom tam připravit jiný systém ve stejném stavu. Problém je ovšem v tom, že nelze určit stav kvantového systému měřením na jediném exempláři. Pokud bychom totiž změřili jednu fyzikální veličinu, řekněme polohu nějaké částice, zničíme tím informaci o konjugované veličině, v našem případě o její hybnosti. Mohli bychom se snažit operovat v rámci Heisenbergových relací neurčitosti a měřit obě veličiny současně za cenu snížení přesnosti. Představme si například harmonický oscilátor o (úhlové) frekvenci ω připravený v nějakém (pro nás neznámém) koherentním stavu. Pokud bychom se 53 pokusili změřit jeho polohu a hybnost a výsledné hodnoty použít k přípravě stavu jiného harmonického oscilátoru, byl by takto vytvořený stav podobný originálu, ale obsahoval by přidaný šum. Množství šumu by odpovídalo fluktuacím termálního stavu o střední energii jednoho energetického kvanta ¯hω, kde ¯h = 1,055×10−34 Js je Planckova konstanta h dělená 2p. Čistší výsledek nám kvantová mechanika s tímto protokolem nedovolí. Tento problém se však dá překonat, pokud by výchozí a cílová stanice sdílely předem připravený kvantově provázaný stav. To je zvláštní stav systému, který se skládá ze dvou (či více) odlehlých podsystémů. Můžeme pak měřit určité kombinované veličiny systému, jehož stav se má teleportovat a provázaného podsystému ve výchozí stanici. Mohou to být například součet poloh a rozdíl hybností - operátory těchto veličin navzájem komutují a lze je tedy měřit současně s libovolnou přesností. Výsledné hodnoty pošleme do cílové stanice, kde se na jejich základě působí na tamní podsystém - operátor jej přesune a „šťouchne do něj tak, aby se přesným způsobem změnila jeho poloha a hybnost. Na konci celé této operace bude stav vstupního systému zničen, ale v cílové stanici bude připraven stav totožný s tím, jaký byl původně ve výchozí stanici. Takováto „kvantová teleportace byla nedávno experimentálně realizována na optických pulsech, přičemž teleportované veličiny byly jednak proměnné odpovídající intenzitě elektrického pole světelného svazku (analogické poloze a hybnosti mechanického systému), ale i polarizační proměnné analogické spinu. V tomto druhém případě dává měření kombinovaného systému ve výchozí stanici dva bity informace (jeden ze čtyř stavů systému složeného ze dvou spinů 1/2), které se pošlou do cílové stanice tak, aby tam mohl být zrekonstruován teleportovaný qubit. O teleportaci kvantových stavů qubitu se dozvíme více v kapitole 2.9 Husté kódování Teleportací „naruby by se mohlo nazývat husté kódování, kdy dvě komunikující strany (řekněme Alice a Bob) na počátku sdílejí kvantově provázané dvojice systémů. Pokud chce Alice poslat nějakou zprávu Bobovi, bude mu předávat qubity například fotony s vhodnou polarizací. Při vhodné manipulaci pak s každým qubitem předá dva bity informace, tedy dvakrát tolik, než co by umožňovalo posílání qubitů bez sdílených provázaných stavů. Klonování Operací, kterou nám kvantová mechanika nedovolí provádět dokonale, je přesné kopírování kvantových stavů, někdy nazývané klonování. To, co je v počítačích naprosto běžnou a nezbytnou součástí mnoha algoritmů, tedy zkopírování dat z jednoho registru do druhého, s qubity učinit nelze. Plyne to z principu superpozice: představme si zařízení, které do prázdného registru x kopíruje informaci o hodnotách 0, 1 z registru y, tedy |0 x|0 y → |0 x|0 y a |0 x|1 y → |1 x|1 y . Pokud by však qubit v registru y byl nějakou superpozicí stavů |0 a |1 , třeba a|0 y + b|1 y, byl by pak výsledný stav obou registrů a|0 x|0 y + b|1 x|1 y, což je něco jiného, než dvě kopie původního stavu, tedy (a|0 x + b|1 x)(a|0 y + b|1 y). Nemožnost přesně 54 kopírovat kvantovou informaci nám dovoluje spolehnout se na kvantovou kryptografii. Nicméně můžeme kvantové stavy klonovat alespoň přibližně s tím, že kopírování vnáší do soustavy jistý šum. Jak zvolit optimální klonovací strategii? Jaké ohrožení soukromí by pak takovéto klonování představovalo? Jak skombinovat teleportaci s klonováním pro „kvantové faxování či „teleklonování ? To jsou některé z mnoha otázek, které dnes kvantová informatika řeší a které jsou podstatné pro případnou kvantovou komunikaci. Motivací je umožnit výměnu dat mezi kvantovými počítači i mezi vzdálenými stranami, které chtějí sdílet tajný kryptografický klíč a potřebují zamezit vlivu ztrát, které jsou podstatné při přenosu na větší vzdálenosti. Shrnutí Kvantová informatika využívá kódování informace pomocí kvantových stavů různých fyzikálních systémů. Důležitým aspektem je zde hlavně princip superpozice stavů. Kvantový počítač by měl být schopen efektivně řešit některé problémy, které by na klasických počítačích trvaly příliš dlouho. Jedná se například o faktorizaci velkých čísel, vyhledávání v rozsáhlých databázích nebo modelování časového vývoje kvantových systémů. Kvantová kryptografie slouží k bezpečnému předávání tajného klíče. K dalším aplikacím kvantové informatiky patří husté kódování či teleportace kvantových stavů. Pojmy k zapamatování • Moorův zákon, • superpozice kvantových stavů, • qubit, • kvantový počítač, • kvantová kryptografie, • kvantová teleportace. Kontrolní otázky 1. Jakým způsobem rostou nároky na paměť počítače v závislosti na velikosti systému při modelování kvantových stavů? 2. Proč je zajímavá otázka jak efektivně faktorizovat velká čísla? 3. Jaké jsou hlavní překážky pro postavení kvantového počítače? 4. Proč lidé považují kvantovou kryptografii za bezpečnou? Úkoly k textu 1. Předpokládejme, že Moorův zákon bude platit ještě neomezenou dobu a že například každý rok a půl se zdvojnásobí velikost paměti, kterou mají naše počítače k dispozici. Řekněme, že jsme nyní schopni uložit 10 GB v paměti počítače pro popis fyzikálního stavu. Za jak dlouho bychom pak byli schopni uložit informaci o klasickém mikrostavu jednoho molu plynu, jestliže pro popis stavu jedné molekuly potřebujeme 5 bytů informace? 55 2. Řekněme, že na zadání jedné amplitudy pravděpodobnosti potřebujeme 5 bytů informace, a že chceme popsat kvantový stav soustavy spinů, z nichž každý má dvě možné projekce. Kolikaspinový systém můžeme popsat za dnešních podmínek (jako v předchozí úloze) a kolikaspinový systém bude možné popsat poté, co půjde do paměti našeho počítače uložit informaci o klasickém mikrostavu jednoho molu částic? Řešení 1. Jeden mol je ∼ 6 × 1023 částic, pro popis stavu potřebujeme ∼ 3 × 1024 bytů informace, tedy asi 3×1014 krát víc, než máme k dispozici. Protože 3×1014 ≈ 248 , musíme celkem 48krát zdvojnásobit naši paměťovou kapacitu, na což bychom při stálé platnosti Moorova zákona potřebovali asi 72 let. 2. Na kvantový popis N-spinového systému musíme zadat 2N amplitud a tedy potřebujeme 5 × 2N bytů paměti. K dispozici máme zatím 1010 bytů, což je zhruba 5×231 , takže můžeme popsat stav systému s asi 31 spiny. Z předchozí úlohy víme, že pro popis klasického mikrostavu jednoho molu potřebujeme 48krát zdvojnásobit naši paměťovou kapacitu. Každé zdvojnásobení paměťové kapacity nám umožní zvýšit počet popisovaných spinů o jeden, takže nakonec bychom mohli pracovat s asi 79 spiny. 2.2 Kvantové bity neboli qubity 2.2.1 Co je to qubit Za kvantový bit můžeme považovat libovolný fyzikální systém, u nějž jsou ve hře dva vzájemně ortogonální kvantové stavy. Z matematického hlediska tyto stavy odpovídají dvěma ortogonálním vektorům v Hilbertově prostoru. Fyzikálně to jsou například dvě ortogonální polarizace fotonu, dva vibrační stavy harmonického oscilátoru (třeba iontu v iontové pasti), dva spinové stavy atomu nebo iontu, případně stav fotonu, který se může nacházet ve dvou cestách nějakého interferometru. Mohou to být i dvě různé hodnoty proudu v supravodivém interferenčním zařízení SQUID (Superconducting QUantum Interference Device), případně dvě hodnoty náboje na supravodivém ostrůvku odděleném od jiných supravodičů tenkou bariérou z nesupravodivého materiálu—tzv. Josephsonovým spojem. Při volbě dvou bázových ortogonálních stavů dvoustavového systému máme velkou volnost: pokud jde o polarizace fotonu, může to být vertikální a horizontální polarizace, nebo polarizace 45◦ šikmo vzhůru a šikmo dolů, nebo levotočivá a pravotočivá kruhová polarizace, případně jakákoliv jiná dvojice lineárně či elipticky polarizovaných stavů. Je dobré si připomenout, že pokud si zvolíme nějakou bázovou dvojici stavů, jakýkoliv jiný stav je možné vyjádřit jako jejich superpozici. Napříplad 56 pravotočivě polarizovaný foton ve stavu |R a levotočivě polarizovaný foton ve stavu |L lze vyjádřit pomocí vertikální polarizace | a horizontální polarizace | ↔ jako |R = 1 √ 2 (| + i| ↔ ) , (2.2.1) |L = 1 √ 2 (i| + | ↔ ) . (2.2.2) Pro bázové stavy přitom platí podmínka ortonormality, tedy skalární součin vektoru se sebou samým je jednička a skalární součin s druhým (ortogonálním) vektorem je nula, | = 1, (2.2.3) | ↔ = 0, (2.2.4) ↔ | = 0, (2.2.5) ↔ | ↔ = 1. (2.2.6) Průvodce studiem Volba vhodné báze je důležitým krokem: může nám velice usnadnit některé výpočty, nebo může být vhodná pro fyzikální implementaci. Pokud pracujeme například se spiny atomů, které se nacházejí v magnetickém poli, je vhodné volit za bázové stavy ty, které mají definovanou projekci spinu do směru magnetického pole: bázové stavy se pak „nehýbou . Pokud bychom zvolili jiný směr než magnetické pole, systém by v daném bázovém stavu nezůstal a osciloval by mezi různými bázovými stavy. Jakmile máme zvolenou bázi, můžeme si její dva vektory označit jako |0 a |1 a každý stav systému vyjádřit jako jejich superpozici |ψ = α|0 + β|1 , kde α a β jsou komplexní čísla splňující |α|2 + |β|2 = 1. Těmto komplexním číslů říkáme často amplitudy. Pokud bychom chtěli vyjádřit stav |ψ pomocí jiné báze, dané vektory |u a |v , můžeme psát |ψ = ˜α|u + ˜β|v . Vztah mezi dvojicí amplitud ˜α, ˜β a α, β je dán skalárními součiny mezi vektory báze |0 , |1 a vektory báze |u , |v . Vynásobíme-li vektor |ψ zleva bázovými vektory |u a |v , dostaneme u|ψ = ˜α u|u + ˜β u|v = ˜α = α u|0 + β u|1 , (2.2.7) v|ψ = ˜α v|u + ˜β v|v = ˜β = α v|0 + β v|1 , (2.2.8) což lze psát v maticovém tvaru ˜α ˜β = U11 U12 U21 U22 α β , (2.2.9) kde matice U daná vztahy U = U11 U12 U21 U22 = u|0 u|1 v|0 v|1 , (2.2.10) se nazývá matice přechodu od báze |0 , |1 k bázi |u , |v . 57 2.2.2 Měření kvantových bitů Měření je v kvantové fyzice svou povahou nedeterministické, to znamená, že často není principiálně možné ze známého stavu předpovědět výsledek měření, ale jen pravděpodobnosti, s jakými ten který výsledek může nastat. Tím, jakou veličinu budeme měřit (např. průmět spinu elektronu do osy x, nebo kruhovou polarizaci fotonu) určujeme bázi, do které se bude stavový vektor promítat. Druhé mocniny absolutních hodnot amplitud u jednotlivých bázových stavů jsou rovny pravděpodobnostem, s jakými zjistíme příslušný výsledek. Průvodce studiem Za tuto tzv. „statistickou interpretaci vlnové funkce dostal Max Born v roce 1954 Nobelovu cenu. Poprvé se objevila v jeho článku z roku 1926, kde se zabývá myšlenkou, že nově vzniklá kvantová mechanika je vhodná i pro popis nestacionárních dějů, například srážek částic. V článku je věta říkající, že nalezená vlnová funkce Φ (zkoumal tam srážku elektronu s atomem) může mít pouze jednu interpretaci, a to pravděpodobnost, že se odražený elektron šíří daným směrem. Pod čarou se pak objevuje poznámka doplněná při korektuře: „Přesnější úvaha ukazuje, že pravděpodobnost je úměrná kvadrátu Φ. Dva řádečky a Nobelova cena z toho—tomu říkám efektivita práce! Pokud uvažujeme například foton ve stavu |ψ = 1 2 |R − i √ 3 2 |L a měříme kruhovou polarizaci fotonu, zjistíme s pravděpodobností |1 2 |2 = 1/4 pravotočivou a s pravděpodobností |i √ 3 2 |2 3/4 levotočivou polarizaci. Pokud bychom měli stav fotonu zadán jako lineárně polarizovaný ve svislém směru | , převedením do báze kruhových polarizací zjistíme, že | = 1√ 2 |R − i√ 2 |L a tím pádem každou z kruhových polarizací bychom naměřili s pravděpodobností 1/2. Je třeba zdůraznit, že ačkoliv qubit může být v nekonečně mnoha různých stavech, jeho měřením nikdy nedostaneme víc než jeden bit informace: jednu ze dvou možných odpovědí. Pokud uvažujeme systém o více qubitech, měření na jednom qubitu může ovlivnit stav zbývajícího systému. Představme si například dva qubity připravené ve stavu |ψ = 1 2 |00 + 1 2 |01 + 1 2 |10 − 1 2 |11 (2.2.11) a měříme stav prvního qubitu v bázi |0 , |1 . Stav daný rovnicí (2.2.11) si můžeme přepsat do tvaru |ψ = 1 √ 2 |0 ⊗ 1 √ 2 (|0 + |1 ) + 1 √ 2 |1 ⊗ 1 √ 2 (|0 − |1 ). (2.2.12) Měření na prvním qubitu nám dá výsledek |0 nebo |1 , každý s pravděpodobností 1/2. Pokud bude výsledek 0, bude druhý qubit ve stavu |+ ≡ 1√ 2 (|0 + |1 ), pokud bude výsledek 1, změní se stav druhého qubitu na |− ≡ 1√ 2 (|0 − |1 ). Pokud 58 bychom tedy nyní měřili stav druhého qubitu v bázi dané stavy |+ a |− , dostali bychom jednoznačný výsledek určený tím, co bylo naměřeno na prvním qubitu. Na druhé straně, pokud bychom druhý qubit měřili v bázi |0 a |1 , byl by výsledek zcela náhodný: s pravděpodobností 1/2 bychom dostali 0 a s pravděpodobností 1/2 bychom dostali 1. Stav z rovnice (2.2.11) však lze přepsat i ve tvaru |ψ = 1 √ 2 1 √ 2 (|0 + |1 ) ⊗ |0 + 1 √ 2 1 √ 2 (|0 − |1 ) ⊗ |1 = 1 √ 2 |+ ⊗ |0 + 1 √ 2 |− ⊗ |1 . (2.2.13) To znamená, že pokud budeme měřit první qubit v bázi |+ a |− , dostane se podle výsledku druhý qubit buď do stavu |0 nebo |1 . Pokud pak měříme stav druhého qubitu v bázi |0 , |1 , bude výsledek jednoznačně určený z měření na prvním qubitu. Pokud bychom však měřili druhý qubit v bázi |+ a |− , bude výsledek zcela náhodný. Stav z rovnice(2.2.11) má zajímavé korelační vlastnostmi mezi dvěma qubity; takovému stavu se říká kvantově provázaný, neboli entanglovaný. Takovéto stavy mají velký význam pro přenos kvantové informace, více se o nich dozvíme v kapitole 2.8. Shrnutí Qubit čili kvantový bit je libovolný systém s dvěma ortogonálními stavy, například spin elektronu či polarizace fotonu. Qubit může nést jeden bit informace, může však být připraven v libovolné superpozici bázových stavů. Měření na qubitu nese všechny rysy kvantového měření: výsledek může být buď předpověditelný nebo náhodný, podle zvolené báze. Měřením však původní kvantový stav zaniká a systém (pokud není též detekcí zničen jako foton při fotodetekci) je připraven ve zcela novém kvantovém stavu. Pojmy k zapamatování • Qubit, • báze, • kvantové měření. Kontrolní otázky 1. Kolik bitů informace můžete předat kolegovi, jestliže mu dáte jeden qubit? 2. Kolik bázových stavů má systém složený z pěti qubitů? 3. S jakými pravděpodobnostmi bychom dostali jednotlivé výsledky, pokud bychom stav z předchozí úlohy měřili v bázi |+ a |− ? 59 Úkoly k textu 1. Qubit je připraven ve stavu i√ 3 |0 − 2 3 |1 . Provedeme na něm měření v bázi |0 , |1 . S jakými pravděpodobnostmi dostaneme jednotlivé výsledky? 2. S jakými pravděpodobnostmi bychom dostali jednotlivé výsledky, pokud bychom stav z předchozí úlohy měřili v bázi |+ a |− ? Řešení 1. |0 s pravděpodobností 1/3 a |1 s pravděpodobností 2/3. 2. |+ i |− s pravděpodobností 1/2. 2.3 Kvantová hradla Pro funkci kvantového počítače je nezbytné, abychom dokázali s qubity manipulovat. Transformace qubitů probíhá v hradlech—to jsou zařízení, ve kterých dochází k nějaké jasně definované změně kvantového stavu qubitu. 2.3.1 Reverzibilita kvantového počítání Každá transformace qubitů odpovídá evoluci kvantového stavu nějakého fyzikálního systému. Pokud je celý systém izolovaný od svého okolí, je takováto transformace popsána unitárním operátorem, tedy operátorem U, pro který platí UU† = = U† U = I, kde I je jednotkový operátor. Pro maximální využití výhod kvantového počítání je zapotřebí pracovat právě s unitárními transformacemi. Protože každý unitární operátor lze invertovat (inverzním operátorem k U je jeho hermitovsky sdružený operátor U† ), je možné každou operaci v kvantovém počítači obrátit a nechat ji běžet nazpět. To je velký rozdíl oproti práci klasických počítačů, který však s sebou přináší určité komplikace. Průvodce studiem Reverzibilním počítáním se zabývali teoretici i na základě klasické fyziky: Fredkin a Toffoli v osmdesátých letech přišli s balistickým modelem reverzibilního počítače. Ten je tvořen soustavou kulečníkových koulí (jejichž přítomnost či nepřítomnost kóduje binární čísla) a překážek tvořících hradla. Koule se pohybují beze ztrát energie podle Newtonových zákonů a narážejí do překážek a do sebe navzájem. Výsledek „výpočtu pak odečteme z poloh koulí na konci jejich dráhy. Protože srážkové děje v klasické fyzice mohou probíhat oběma směry, můžeme nechat takovýto počítač „běžet zpět , aby z výsledných dat udělal data vstupní. 60 FANAOUT: kopírování Operace FANOUT (tedy „rozfouknutí ) vytváří dvě kopie vstupního bitu. Aby tato operace mohla být unitární, znamená to, že velikost vstupu musí být rovna velikosti výstupu (měřeno dimenzí Hilbertova prostoru vstupního a výstupního stavu). Protože na výstupu jsou dva bity, musí být vstupní bit doplněn dalším pomocným bitem, který je nastaven na nějakou pevnou hodnotu. Příklad realizace FANOUT ukážeme v odstavci 2.3.4. Operaci FANOUT můžeme pro kvantové počítání sestrojit, nicméně je třeba si uvědomit, že nám dovoluje pouze vytvořit z bitu 0 dva nulové bity, případně z bitu 1 dva bity o hodnotě 1. Nemůžeme však „klonovat kvantové stavy, tedy vytvářet kopie libovolných superpozic—viz odstavec 2.1.5. Nulování Operace, která by z libovolné vstupní hodnoty vytvořila nulu, neboli jakýkoliv vstupní stav transformovala na |0 , není unitární. Znamená to tedy, že obecně nemůžeme v kvantovém počítači vynulovat nepotřebné bity a volné registry použít pro ukládání užitečných dat. Pomoci si můžeme pouze tehdy, pokud máme dvě (či více) kopií nějakého bitu. V tomto případě můžeme použít inverzi k operaci FANOUT a ze dvou kopií vytvořit jedinou s tím, že jeden bit se vynuluje. Pro vynulování poslední kopie daného bitu (o neznámé hodnotě) však nemůžeme použít unitární transformaci. Takovéto operaci by odpovídalo „stlačování fázového prostoru nebo snižování dimenze Hilbertova prostoru našeho fyzikálního systému. Z termodynamického hlediska by to odpovídalo poklesu entropie systému. Podle druhého termodynamického zákona však entropie izolovaného systému klesat nemůže, entropie se můžeme zbavit pouze tím, že ji předáme nějakému jinému fyzikálnímu systému. Pokud budeme předávat tuto entropii okolí, které je v termodynamické rovnováze při teplotě T, může se jeho entropie zvýšit o ∆S pouze tehdy, pokud mu dodáme energii (ve formě tepla) T∆S. Při nulování neznámého bitu se entropie daného registru sníží o kB ln 2, kde kB = 1,38×10−23 J/K je Boltzmannova konstanta. Proč? Neznámý bit nabývá hodnoty 0 s pravděpodobností 1/2 a hodnoty 1 také s pravděpodobností 1/2. Entropie takového registru je tedy S1 = −kB n pn ln pn = − −kB[1 2 ln 1 2 + 1 2 ln 1 2 ] = kB ln 2. Bit s hodnotou 0 odpovídá systému s nulovou entropií: p0 = 1 a p1 = 0, tedy S1 = 0. Celkem tedy entropie registru poklesne o ∆S = kB ln 2. Předáme-li takovou entropii do okolí o teplotě T, musíme disipovat teplo kBT ln 2. Rolf Landauer (1927 - 1999) z IBM, největší přínos má ve zkoumání fyzikálních mezí pro naše výpočetní možnosti. Tomuto tvrzení se říká Landauerův princip: pro vynulování neznámého bitu je nutno disipovat alespoň kBT ln 2 tepla. V současných počítačích se disipuje mnohem více tepla (bity v nich nejsou kódovány pomocí dvou kvantových stavů elementárních „dvoustavových systémů, ale pomocí stavů mnohem větších objektů). Landauerův princip udává určitou mez pro minimální disipaci tepla, ke které se můžeme blížit, pokud naše počítače budou v maximální míře pracovat s reverzibilními kroky. 61 Průvodce studiem Landauerův princip se ukazuje také jako konečné rozluštění záhady Maxwellova démona. Maxwell uvažoval inteligentní bytůstku, která by oddělovala v rozdělené nádobě s plynem rychlé molekuly od pomalejších: po čase by tak v jedné části byla vyšší teplota než ve druhé. To by nám pak umožňovalo zapojit mezi obě části tepelný stroj, který by část tepla přeměnit na užitečnou práci. To je ale proti druhému termodynamickému zákonu. Někde asi musí být chyba - a přes stovku let a řadu falešných stop trvalo, než fyzikové zjistili kde: při každém měření rychlosti molekul Maxwellův démon zaplňuje svoji paměť. Až svou paměť zaplní úplně, bude ji muset vynulovat. Ale na to musí podle Landauerova principu disipovat energii (znehodnotit ji na teplo). Znehodnocené energie bude nakonec přinejmenším tolik, kolik jsme získali užitečné práce. Pro kvantové počítání to znamená, že pokud ireverzibilní kroky nepřicházejí v úvahu (například pro ztrátu koherence), musíme mít v zásobě dostatek volných bitů, do kterých se v průběhu výpočtů odkládá „odpad . 2.3.2 Jednobitová hradla Hradla I, X, Y a Z Nejjednoduššímm typem hradel jsou taková, která mění stav jediného qubitu. Ta zcela nejjednodušší, označená jako operátory I, X, Y a Z transformují qubity následujícím způsobem: I : |0 → |0 , |1 → |1 , (2.3.1) X : |0 → |1 , |1 → |0 , (2.3.2) Y : |0 → −|1 , |1 → |0 , (2.3.3) Z : |0 → |0 , |1 → −|1 . (2.3.4) Tyto transformace lze také vyjádřit pomocí matic jako I = 1 0 0 1 , X = 0 1 1 0 , Y = 0 1 −1 0 , Z = 1 0 0 −1 , Transformace I je identita a „hradlo tu má triviální funkci, hradlo X je negace, hradlo Z představuje fázový posuv o p radiánů a hradlo Y = ZX je kombinací 62 posledních dvou. Uvedená hradla se někdy znázorňují pomocí symbolů na obr. 2.1. Na schématu může proměnná x nabývat hodnot 0 a 1, přičemž součet ⊕ se definuje jako 0 ⊕ 0 = 1 ⊕ 1 = 0 a 0 ⊕ 1 = 1 ⊕ 0 = 1. I Y |x |x |x |x |x (−1)x+1 |x⊕1 (−1)x |x |x ⊕ 1 X Z Obr. 2.1: Schéma jednobitových hradel I, X, Y a Z. Hadamardovo hradlo Tuto čtveřici jednobitových transformací je třeba doplnit o hradlo, připravující superpozici bázových stavů podle předpisu H : |0 → 1 √ 2 (|0 + |1 ) |1 → 1 √ 2 (|0 − |1 ) , (2.3.5) což lze vyjádřit pomocí matice H = 1 √ 2 1 1 1 −1 . (2.3.6) Toto je takzvaná Hadamardova transformace. Její schéma, spolu s působením na qubit ve stavu |0 je na obr. 2.2. Algebraicky můžeme Hadamardovu transformaci vyjádřit též jako H|x = 1√ 2 [(−1)x |x + |x ⊕ 1 ). 1√ 2 (|0 + |1 )|0 H Obr. 2.2: Působení Hadamardovy transformace na nulový qubit. Obzvláště důležitou aplikací této transformace je vytváření mnohabitových stavů, pokud operátor H působí na každý z N vynulovaných bitů: (H ⊗ H ⊗ . . . ⊗ H)|00 . . . 0 = 1 √ 2N (|0 + |1 ) ⊗ (|0 + |1 ) ⊗ . . . ⊗ (|0 + |1 ) = 1 √ 2N 2N −1 x=0 |x , (2.3.7) kde |x jsou N-bitové stavy a x probíhá všechny hodnoty čísel od 0 do 2N − 1 v binárním tvaru. Takovéto Hadamardově transformaci aplikované na Nqubitů se říká Walshova, případně Walshova-Hadamardova transformace. Aplikací Walshovy transformace připravujeme vstupní registr v superpozici všech možných čísel od 0 do 2N − 1, přičemž každé číslo tu vystupuje se stejnou vahou. 63 Fázový posuv Zobecněním hradla Z na libovolné úhly je fázový posuv o hodnotu φ, který můžeme vyjádřit transformací Z(φ) : |0 → |0 , |1 → eiφ |1 , (2.3.8) případně pomocí matice jako Z(φ) = 1 0 0 eiφ . Rotace qubitu Zobecněním hradel I a Y je rotace qubitu o úhel θ, daná maticí U(θ) = cos θ 2 sin θ 2 − sin θ 2 cos θ 2 . (2.3.9) Jak si můžeme ověřit, pro θ = 0 dostáváme U(0) = I a pro θ = p dostáváme U(p) = Y . 2.3.3 Dvoubitová hradla Řízená negace CNOT Aby mohl kvantový počítač fungovat, je nutné, aby mohl jeden bit ovlivňovat druhý. Operací, která patří k těm nejdůležitějším, je kontrolovaná (řízená) negace, zvaná CNOT (controlled not): podle toho, v jakém stavu je první qubit, se druhý qubit buď nezmění, nebo se neguje. Formálně ji můžeme popsat jako CNOT : |00 → |00 |01 → |01 |10 → |11 |11 → |10 , (2.3.10) případně pomocí matice CNOT =      1 0 0 0 0 1 0 0 0 0 0 1 0 0 1 0      . (2.3.11) Schéma hradla CNOT je na obrázku 2.3. 64 |x |x |y |x ⊕ y Obr. 2.3: Schéma hradla CNOT. Řízený fázový posuv Jiným příkladem dvoubitového hradla je řízený fázový posuv: fáze stavu se změní pouze tehdy, když oba qubity budou nabývat hodnotu 1: B(φ) : |00 → |00 |01 → |01 |10 → |10 |11 → eiφ |11 , (2.3.12) případně pomocí matice B(φ) =      1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 eiφ      , (2.3.13) což odpovídá schématu na obr. 2.4. |x |y eixyφ |x eixyφ |yφ Obr. 2.4: Schéma řízeného fázového posuvu. Pro kvantové počítání nejdůležitějším případem řízeného fázového posuvu je takový, který odpovídá hodnotě φ = p, tedy matici B(p) =      1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 −1      . (2.3.14) Jak se můžeme přesvědčit, lze z tohoto hradla a pomocí jednobitových Hadamardových hradel sestavit hradlo CNOT, což odpovídá obr. 2.5. 65 H H =π Obr. 2.5: Řízená negace sestavená ze dvou Hadamardových hradel a řízeného fázového posuvu. 2.3.4 Kombinace kvantových hradel Hradlo SWAP Řetězením jednobitových a dvoubitových kvantových hradel můžeme získávat další kvantové obvody—některé velmi jednoduché, jiné komplikovanější. Například ze tří hradel CNOT můžeme sestrojit hradlo SWAP, tedy výměnu bitu mezi dvěma registry, viz obr. 2.6. |x |x|x|y |y|x |y |y = Obr. 2.6: Výměna bitu mezi registry pomocí tří hradel CNOT. Toffoliho hradlo Toffoliho hradlo je tříbitové hradlo a patří k těm nejdůležitějším pro kvantové výpočty. Jedná se o negaci řízenou dvěma vstupními bity: právě tehdy, když oba řídící bity nabývají hodnotu 1, překlopí se hodnota řízeného bitu. Schématicky znázorňujeme Toffoliho hradlo jako na obrázku 2.7. Symbolicky můžeme zapsat funkci |y |x |y |x |z |z ⊕ (xy) Obr. 2.7: Toffoliho hradlo překlápí bit z právě tehdy když oba bity x a y nabývají hodnoty 1. Toffoliho hradla vztahem (x,y,z) → (x,y,z⊗(xy)). Pokud na jeden vstupní bit přivádíme pevně danou hodnotu, chová se Toffoliho hradlo jako dvoubitové hradlo, pokud na dva vstupní bity přivedeme pevné hodnoty, můžeme dostat hradla transformující 66 jediný bit. Přitom o jaké hradlo půjde, můžeme zvolit výběrem vstupních bitů: z ⊗ (xy) =    xy pro z = 0 (AND) x ⊕ z pro y = 1 (CNOT) x pro y = z = 1 (NOT) x pro y = 1, z = 0 (FANOUT) (2.3.15) Kvantové Toffoliho hradlo lze sestrojit pomocí řízených negací a jednobitových hradel, příklad je na obrázku 2.8. |z |y |x |y |z ⊕ (xy) |x U(π/4) U(π/4) U(−π/4) U(−π/4) Obr. 2.8: Toffoliho hradlo sestrojené pomocí tří hradel CNOT a čtyř hradel rotujících qubit o úhel p/4. (Ve skutečnosti přidává tento obvod v některých případech k řízenému bitu fázový posun, ten je však v případě potřeby snadno odstranitelný.) Z Toffoliho hradel a z řízené negace lze sestavit například sčítací obvod, jak ukazuje obr. 2.9. V něm zůstává první bit x nezměněn, druhý bit nese na výstupu |x |y |x |0 |x ⊕ y |xy Obr. 2.9: Kvantový sčítací obvod pro dva qubity sestavený z Toffoliho hradla a řízené negace. hodnotu součtu bitů x a y modulo 2, x ⊕ y, a třetí bit slouží jako přenos: jsou-li oba vstupní bity rovny jedné, je jejich součet (modulo 2) 0 a „jedna jde dál . Řetězením dvoubitových sčítacích hradel pak můžeme sestavit hradlo sčítající libovolně velká vstupní čísla. Podobným způsobem pak můžeme konstruovat řadu dalších matematických operací. Průvodce studiem Poznamenejme, že pro sestavení univerzálního klasického počítače potřebujeme jako stavebí kameny několik typů elementárních hradel, vždy však mezi nimi bude alespoň jedno tříbitové hradlo—např. Toffoliho. U kvantových obvodů však pro univerzální obvod vystačíme s jedno- a dvoubitovými hradly: díky možnosti vytvářet superpozice stavů z nich lze sestrojit Toffoliho hradlo, jak jsme ukázali na obr. 2.8. Něco takového u klasických obvodů není možné. 67 2.3.5 Kvantová Fourierova transformace Mezi ostatními funkcemi vytvořenými z elementárních hradel patří k nejdůležitějším Fourierova transformace, která z N-bitového kvantového stavu |x vytváří superpozici F : |x → 1 2N/2 2N −1 y=0 exp i 2pxy 2N |y . (2.3.16) Všimněme si, že tato transformace přeměňuje každý bázový stav (tedy každý stav odpovídající nějakému binárnímu číslu) v superpozici všech bázových stavů, zastoupených se stejnou vahou. Chová se tedy podobně jako Walshova-Hadamardova transformace, kromě triviálního případu N = 1 však působí odlišně. Konkrétně pro N = 2 dostáváme pro Fourierovu transformaci F : |00 → 1 2 (|00 + |01 + |10 + |11 ) , |01 → 1 2 (|00 + i|01 − |10 − i|11 ) , |10 → 1 2 (|00 − |01 + |10 − |11 ) , |11 → 1 2 (|00 − i|01 − |10 + i|11 ) , (2.3.17) zatímco Walshova-Hadamardova transformace působí jako H : |00 → 1 2 (|00 + |01 + |10 + |11 ) , |01 → 1 2 (|00 − |01 + |10 − |11 ) , |10 → 1 2 (|00 + |01 − |10 − |11 ) , |11 → 1 2 (|00 − |01 − |10 + |11 ) . (2.3.18) Průvodce studiem Ve fyzice hraje Fourierova transformace velmi důležitou roli, zjednodušeně můžeme říci, že nás informuje o tom, jak jsou v dané funkci (nebo signálu) zasoupeny různé frekvence. Rozklad světelného svazku hranolem do spektrálních komponent je svým způsobem fyzikální implementace Fourierovy transformace. V řadě případů je jednodušší pracovat „ve frekvenční oblasti spíše než přímo s časovým průběhem. Kvantová Fourierova transformace je v jistém ohledu analogií takovéhoto spektrálního rozkladu pro rozklad nějakého stavu do různých bází. Vypišme ještě explicitně kvantovou Fourierovu transformaci pro tříbitové stavy: F : 68 |000 → 1 √ 8 (|000 + |001 + |010 + |011 + |100 + |101 + |110 + |111 ) , |001 → 1 √ 8 |000 + eip 4 |001 + i|010 + ei3p 4 |011 − |100 + e−i3p 4 |101 − i|110 + e−ip 4 |111 , |010 → 1 √ 8 (|000 + i|001 − |010 − i|011 + |100 + i|101 − |110 − i|111 ) , |011 → 1 √ 8 |000 + ei3p 4 |001 − i|010 + eip 4 |011 − |100 + e−ip 4 |101 + i|110 + e−i3p 4 |111 , |100 → 1 √ 8 (|000 − |001 + |010 − |011 + |100 − |101 + |110 − |111 ) , |101 → 1 √ 8 |000 + e−i3p 4 |001 + i|010 + eip 4 |011 − |100 + eip 4 |101 − i|110 + ei3p 4 |111 , |110 → 1 √ 8 (|000 − i|001 − |010 + i|011 + |100 − i|101 − |110 + i|111 ) , |111 → 1 √ 8 |000 + ei−p 4 |001 − i|010 + e−i3p 4 |011 − |100 + e−i3p 4 |101 + i|110 + eip 4 |111 . (2.3.19) Hradlo vytvářející takovouto transformaci lze sestrojit z jednobitových Hadamardových hradel a z dvoubitových hradel řízeného fázového posuvu o úhly p/2 a p/4, jak je vidět z obrázku 2.10. Pro sestrojení vícebitové Fourierovy transformace bychom H H H 0 1 2 π/2 π/4 π/2 Obr. 2.10: Kvantový obvod vytvářející Fourierovu transformaci na třech vstupních bitech. Konečný výsledek je třeba ještě bitově invertovat (nebo bity číst v opačném pořadí), což je triviální operace. postupovali induktivně: pro přidávání dalšího bitu bychom napojili příslušný registr hradlem řízeného fázového posuvu na všechny předchozí registry a na závěr připojili hradlo Hadamardovy transformace. Přitom fázový posuv mezi k-tým a j-tým bitem je φ = p/2|k−j| . Shrnutí Při zpracování kvantové informace se musí brát v potaz reverzibilita kvantových transformací. Některé operace, které běžně provádějí naše počítače, jsou ireverzibilní a při kvantovém zpracování informace je třeba je vhodně modifikovat. Kvantová hradla provádějí transformace stavů qubitů. Nejjednodušší hradla transformují jednobitové stavy, u vícebitových hradel hodnoty jednoho či více bitů ovlivňují jiné bity. K nejdůležitějším jednobitovým hradlům patří ty, které realizují fázový posuv, rotace a Hadamardovu transformaci. Nejdůležitější dvoubitová hradla realizují řízenou negaci a řízený fázový posuv. Kombinací jednobitových a dvoubitových hradel 69 lze získat složitější kvantové obvody. Toffoliho hradlo, které může posloužit jako základní kámen pro sčítací obvody, lze sestrojit kombinací hradel řízené negace a jednobitových rotací. Ze sčítacích obvodů pak můžeme, podobně jako v klasické informatice sestavit obvod realizující komplikovanější funkce. Kromě toho patří k nejdůležitějším mnohobitovým obvodům ten, kterým se uskutečňuje kvantová Fourierova transformace. Pojmy k zapamatování • Reverzibilita, operace FANOUT, mazání informace, Landauerův princip, • jednobitová hradla: fázový posuv, rotace, Hadamardova transformace, • Walshova-Hadamardova transformace, • dvoubitová hradla: řízená negace a řízený fázový posuv, • vícebitová hradla: Toffoliho hradlo, sčítací obvod, kvantová Fourierova transfor- mace Kontrolní otázky 1. Jak vypadají matice popisující Hadamardovu transformaci, fázový posuv o p/4 a negaci jednoho qubitu? 2. Jak vypadají matice popisující dvoubitovou Walshovu-Hadamardovu transformaci, řízenou negaci, řízený fázový posuv a Fourierovu transformaci? Úkoly k textu 1. Uvažujte dvoubitové hradlo řízeného fázového posuvu B(p). Na jeho vstup přivedeme stav Ψ− = 1√ 2 (|10 − |01 ). Jak bude vypadat výstupní stav? Jak by se změnil tímto hradlem stav stav Φ+ = 1√ 2 (|00 + |11 )? 2. Uvažujme dvoubitové hradlo, které je sestavené takto: první bit projde Hadamardovou transformací, následuje řízený fázový posuv B(p) a poté druhý bit projde Hadamardovou transformací. Jak se změní vstupní stav Ψ− = = 1√ 2 (|10 − |01 )? Řešení 1. 1√ 2 (|10 − |01 ) → 1√ 2 (|10 − |01 ) (stav Ψ− se nezmění), 1√ 2 (|00 + |11 ) → 1√ 2 (|00 − |11 ) (stav Φ+ se změní na Φ− ). 2. Výsledný stav bude 1√ 2 (|0 − |1 ) ⊗ |1 . 2.4 Shorův algoritmus: nalezení periody diskrétní funkce Peter W. Shor z AT&T laboratoří, objevitel kvantového algoritmu pro faktorizaci čísel. Jsme-li schopni pomocí jednoduchých kvantových hradel realizovat výpočet nějaké funkce (například pomocí sčítacích obvodů a podobně) a pokud máme sestro- 70 jený obvod pro kvantovou Fourierovu transformaci, můžeme tímto zařízením zjistit, zda je zadaná funkce periodická a jakou má periodu. Postup je následující. Nejprve připravíme N vstupních bitů v superpozici všech stavů, |0 → 1 2N/2 2N −1 x=0 |x . (2.4.1) Toho nejsnáze dosáhneme Walshovou-Hadamardovou transformací (viz rovnice (2.3.7)), tedy použitím jednoho Hadamardova hradla pro každý vstupní, původně vynulovaný bit. Poté aplikujeme obvod „počítající funkci f, |x; 0 → |x; f(x) : 1 2N/2 2N −1 x=0 |x; 0 → 1 2N/2 2N −1 x=0 |x; f(x) . (2.4.2) Protože jsme měli na vstupu superpozici všech 2N bázových stavů, je jediným průchodem přes posloupnost hradel vypočítána funkce f pro všech 2N hodnot jejího argumentu. Toto je krok, který nám klasické počítače neumožní. Nyní na vstupní bity x aplikujeme kvantovou Fourierovu transformaci (viz odst. 2.3.5). Stav se tak změní 1 2N/2 2N −1 x=0 |x; f(x) → 1 2N/2 2N −1 x=0 2N −1 y=0 exp i 2pxy 2N |y; f(x) . (2.4.3) Předpokládejme nyní, že funkce f má periodu r, tedy f(x + r) = f(x). Tím pádem můžeme stav |f(x) „vytknout a výraz 2N −1 x=0 exp i2pxy 2N |f(x) psát jako 2N −1 x=0 exp i 2pxy 2N |f(x) = r−1 x=0 |f(x) n exp i 2p(x + nr)y 2N = r−1 x=0 exp i 2pxy 2N |f(x) n exp i 2pnry 2N ,(2.4.4) kde sčítání přes n probíhá od 0 do [(2N − 1 − x)/r] a výraz [z] znamená celou část čísla z. Výsledný stav z (2.4.3) lze tedy zapsat ve tvaru 1 2N/2 2N −1 y=0 |y r−1 x=0 exp i 2pxy 2N |f(x) n exp i 2pnry 2N . (2.4.5) Suma na konci tohoto výrazu nabývá velkých hodnot pouze tehdy, když je příslušná exponenciála přibližně rovna 1 nezávisle na n, tedy pokud y nabývá hodnot blízkých y ≈ 2N r m, (2.4.6) kde m je nějaké celé číslo. Pro jiné hodnoty mají příspěvky exp i2pnry 2N pro různé hodnoty n různé fáze a interferencí se buď úplně nebo částečně vyruší. 71 0 5 10 15 20 25 30 0 0.05 0.1 0.15 0.2 y P(y) r=5 N=5 0 5 10 15 20 25 30 0 0.02 0.04 0.06 0.08 0.1 0.12 y P(y) N=5 r=10 Obr. 2.11: Pravděpodobnosti detekce pětibitového (N = 5) stavu |y v případě, že funkce f(x) má periodu r = 5 a r = 10. Pokud nyní provedeme měření na vstupním registru, výsledkem bude nějaká hodnota y z intervalu mezi 0 a 2N −1. Ovšem hodnoty y, pro které se příspěvky exponenciál v sumě v (2.4.5) vyruší, nastanou s téměř nulovou pravděpodobností. S největší pravděpodobností naměříme některou z hodnot y, pro kterou je splňen vztah (2.4.6). Průběh těchto pravděpodobností je ilustrován na obr. 2.11 a 2.12. Z naměřené hodnoty y a vztahu (2.4.6) pak můžeme tipovat jaká je perioda r; dosazením několika málo argumentů do funkce f si pak můžeme správnost tipu ověřit. 0 50 100 150 200 250 0 0.05 0.1 0.15 0.2 y P(y) r=5 N=8 0 50 100 150 200 250 0 0.02 0.04 0.06 0.08 0.1 y P(y) r=10 N=8 Obr. 2.12: Pravděpodobnosti detekce osmibitového (N = 8) stavu |y v případě, že funkce f(x) má periodu r = 5 a r = 10. Průvodce studiem Je dobré si uvědomit několik skutečností: 1. Kvantové počítání je „pravděpodobnostní : výsledek měření je náhodný, i tak nám však poskytuje cennou informaci. Pokud bychom se snažili zjistit periodu funkce f na klasickém počítači, museli bychom funkční hodnotu spočítat mnohokrát, řádově ∼ 2N -krát. Při kvantovém počítání stačí několik málo dosazení do funkce realizované pomocí kvantových hradel a poté dosazení vytipovaných hodnot do funkce počítané na klasickém počítači. 72 2. Je zajímavé, že poté, co jsme kvantově napočítali hodnoty funkce f, x |x,0 → x |x,f(x) , již registr s funkčními hodnotami |f(x) nepotřebujeme a můžeme jej vymazat. Vše zajímavé se nyní odehrává již jen v „nedotčeném registru argumentů |x ! To je důsledek použití superpozice bázových stavů na vstupu. Shrnutí Pro nalezení periody nějaké diskrétní funkce nejprve připravíme superpozici stavů odpovídajících všem vstupním argumentům, pak tento stav použijeme jako argument pro výpočet funkce f a poté provedeme výpočet kvantové Fourierovy transformace na registru agrumentů. Nakonec změříme hodnotu v registru argumentů (zatímco stav v registru funkčních hodnot můžeme zcela ignorovat) a z ní se pokusíme vytipovat hledanou periodu. Správnost tipu pak ověříme přímým dosazením do funkce. Protože těchto úkonů je mnohem méně, než kolik bychom museli učinit při ∼ 2N -násobném dosazování do funkce, představuje tento algoritmus podstatnou úsporu při hledání periody. Pojmy k zapamatování • registr argumentů, • registr funkčních hodnot, • kvantová Fourierova transformace, • měření stavu 2.5 Faktorizace čísel Předešlé výsledky mají velký význam pro efektivní faktorizaci velkých čísel. To by nám umožnilo efektivně luštit tzv RSA kryptografický kód, (viz příloha B). Předpokládejme, že hodláme faktorizovat číslo N. Bude nám stačit nalézt jediný faktor, což redukuje úlohu na jednodušší případ. Můžeme postupovat následujícím způsobem: 1. Zvolíme číslo z a ověříme u něj, že je nesoudělné s N. To lze velmi snadno tzv. Euklidovým algoritmem. Pokud je číslo z soudělné s N, máme první faktor čísla N—největší společný dělitel N a z. 2. Na kvantovém obvodu sestrojíme funkci f(x) = zx (modN). Tuto funkci si můžeme zapsat jako posloupnost x 0 1 2 . . . r − 1 r r + 1 r + 2 . . . 2r 2r + 1 . . . f(x) 1 z z2 . . . zr−1 zr zr+1 zr+2 . . . z2r z2r+1 . . . f(x) 1 z z2 . . . zr−1 1 z z2 . . . 1 z . . . Druhý řádek je tu psán jako posloupnost mocnin, které je však třeba brát modulo N. Ve třetím řádku jsou již zapsány funkční hodnoty tak, že je zřejmé, že funkce je periodická s periodou r: číslo r je prvním netriviálním exponentem, pro který platí zr (modN) = 1. 73 3. Pomocí kvantového algoritmu (viz předchozí kapitola) nalezneme hodnotu periody r. Pokud je toto číslo sudé, pokračujeme k bodu 4, pokud není, vrátíme se k bodu 1 a volíme novou hodnotu z. (Při každém opakování máme 50% šanci dostat sudý výsledek, takže takovýchto opakování nebude mnoho.) 4. Po nalezení sudého r můžeme vztah zr (modN) = 1 přepsat ve tvaru zr/2 2 − 1 = 0(modN), zr/2 + 1 zr/2 − 1 = 0(modN). (2.5.1) To znamená, že součin na levé straně rovnice (2.5.1) je násobkem čísla N a tudíž buď zr/2 + 1 nebo zr/2 − 1 (případně obě čísla) má společného dělitele s N. 5. Nakonec Euklidovým nalezneme společného dělitele zr/2 + 1 a N, případně společného dělitele zr/2 − 1 a N. Toto číslo tedy řeší náš problém. Průvodce studiem V roce 1977 byla vypsána cena 100 dolarů pro toho, kdo dokáže faktorizovat číslo N = 114 381 625 757 888 867 669 235 779 976 146 612 010 218 296 721 242 362 562 561 842 935 706 935 245 733 897 830 597 123 563 958 705 058 989 075 147 599 290 026 879 543 541. Řešení přišlo až v roce 1994, kdy 600-členný tým, který spojil síly svých výpočetních prostředků ohlásil výsledek: N = pq, kde q = 3 490 529 510 847 650 949 147 849 619 903 898 133 417 764 638 493 387 843 990 820 577 a p = 32 769 132 993 266 709 549 961 988 190 834 461 413 177 642 967 992 942 539 798 288 533. Pro ilustraci si zkusme tímto postupem faktorizovat číslo 15. Jako číslo z si zvolme třeba 7. Jako posloupnost funkčních hodnot f(x) = 7x (mod 15) pak dostáváme x 0 1 2 3 4 5 6 7 8 9 10 . . . f(x) 1 7 4 13 1 7 4 13 1 7 4 . . . Je vidět, že perioda funkce f(x) je r = 4, což je sudé číslo. Můžeme tedy psát 74/2 + 1 74/2 − 1 = 0(mod 15), 50 × 48 = 0(mod 15). (2.5.2) Tedy buď 50 nebo 48 musí mít s číslem 15 společného dělitele. Jak vidíme (nebo ověříme Euklidovým algoritmem), největší společný dělitel čísel 50 a 15 je 5 a tedy 5 je dělitelem čísla 15. Zkusme si ještě faktorizovat číslo 391. Za číslo z si zvolme například 4. Pro posloupnost funkčních hodnot f(x) = 4x (mod 391) je pak již výhodné využít počítače. Dostáváme pak x 0 1 2 3 4 5 6 . . . 42 43 44 45 . . . f(x) 1 4 16 64 256 242 186 . . . 220 98 1 4 . . . 74 Perioda funkce f(x) je tedy r = 44, což je opět sudé číslo. Můžeme tedy psát 422 + 1 422 − 1 = 0(mod 391), 17 592 186 044 417 × 17 592 186 044 415 = 0(mod 391). (2.5.3) Euklidovým algoritmem zjistíme, že čísla 17 592 186 044 417 a 391 jsou nesoudělná, kdežto čísla 17 592 186 044 415 a 391 mají společného dělitele 23. Číslo 23 je tedy jedním z hledaných dělitelů a číslo 391 lze faktorizovat jako 391 = 23 × 17. Shrnutí Efektivní faktorizace velkých čísel je podstatná pro luštění kryptografického RSA kódu. Pro nalezení dělitelů čísla N zvolíme s ním nesoudělné číslo z a nalezneme periodu funkce f(x) = zx (modN) jako nějaké sudé číslo r. Kvantový počítač je v tomto kroku mnohem efektivnější než klasické počítače. Z periody r pak snadno zjistíme alespoň jednoho dělitele čísla N. 2.6 Groverův algoritmus: vyhledávání v neuspořádaných seznamech Lov Grover z Bellových laboratoří, autor kvantového vyhledávacího algoritmu. 2.6.1 Popis algoritmu Groverův algoritmus pomáhá vyhledat v nesetříděném seznamu o délce N číslo x, které splňuje předem zadanou podmínku. Předpokládejme, že N je číslo splňující 2N = N a nechť Up je kvantový obvod, který dává jako výstup odpověď, zda vstup splňuje naši podmínku, Up : |x; 0 → |x,P(x) , (2.6.1) kde P(x) = 1 pokud x splňuje podmínku 0 pokud x nesplňuje podmínku (2.6.2) Registr pro funkční hodnoty tedy může být pouze jednobitový. Kvantové vyhledávání probíhá takto: 1. Vstupní registr připravíme v superpozici všech bázových stavů (například z vynulovaného registru pomocí Walshovy-Hadamardovy transformace ), |0 → 1 2N/2 2N −1 x=0 |x . (2.6.3) 2. Na takto připravený vstupní stav aplikujeme kvantový obvod Up, získáme tak stav 1 2N/2 2N −1 x=0 |x; P(x) . (2.6.4) 75 Znamená to, že pokud je na výstupu registru pro funkční hodnotu jednička, musí být v registru pro argument funkce hodnota, která splňuje naši podmínku, nazvěme ji x0. Pokud bychom tedy provedli měření na registru pro funkční hodnotu a získali hodnotu 1, měření na registru pro x by nám dalo hledané číslo x0. Problém však je, že výsledek 1 na registru pro funkční hodnotu získáme s mizivě malou pravděpodobností ∼ 1/2N . Groverův trik spočívá v transformaci kvantového stavu z rovnice(2.6.4) tak, aby vzrostla velikost amplitudy pravděpodobnosti pro stav |x0; 1 . To se děje následujícím způsobem: 3. Provedeme transformaci stavu, která změní amplitudu aj na −aj u těch bázových stavů |xj , pro které P(xj) = 1. 4. Transformujeme stav tak, že se amplitudy pravděpodobnosti „překlápí kolem střední hodnoty amplitud. Tyto dva kroky jsou znázorněny na obr. 2.13. střední hodnota střední hodnota střední hodnota x x x x x x 0 0 0 Obr. 2.13: Překlápění amplitud v Groverově vyhledávacím algoritmu. 5. Posloupnost operací - aplikace operátoru Up → překlopení amplitudy aj na −aj → překlopení amplitud kolem jejich střední hodnoty opakujeme p 4 2N/2 -krát. 6. V registru pro argumenty funkce provedeme měření a odečteme výsledek. 2.6.2 Překlápění kolem střední hodnoty amplitud Abychom mohli aplikovat Groverův algoritmus, musíme dokázat pomocí základních kvantových hradel sestavit obvody překlápějící amplitudy pravděpodobností. Pro překlápění amplitud kolem jejich střední hodnoty potřebujeme zajistit transfor- maci 2N −1 k=0 ak|xk → 2N −1 k=0 (2A − ak)|xk , (2.6.5) 76 kde A označuje střední hodnotu amplitud ak. Tato operace odpovídá matici 2N ×2N D =      21−N − 1 21−N . . . 21−N 21−N 21−N − 1 . . . 21−N . . . . . . . . . . . . 21−N 21−N . . . 21−N − 1      . (2.6.6) Protože platí, že DD† = I, je matice D unitární a může odpovídat nějakému vývoji kvantových stavů. Naším úkolem je nyní ukázat, že tuto matici lze efektivně implementovat pomocí nepříliš vysokého počtu elementárních hradel (jejich počet by měl být úměrný N a ne dimenzi matice 2N ). Grover ukázal, že matici D lze zapsat ve tvaru D = WRW, kde W je WalshovaHadamardova transformace a R je fázový posuv o p u jediného bitu, R =      1 0 . . . 0 0 −1 . . . 0 . . . . . . . . . . . . 0 0 . . . −1      . (2.6.7) Je užitečné si napsat matici R jako R = R − I, kde R =         2 0 0 . . . 0 0 0 0 . . . 0 0 0 2 . . . 0 . . . . . . . . . . . . . . . 0 0 0 . . . 0         . (2.6.8) Protože WRW je W(R − I)W = WR W − I a jak si můžeme ověřit, platí WR W =      21−N 21−N . . . 21−N 21−N 21−N . . . 21−N . . . . . . . . . . . . 21−N 21−N . . . 21−N      . (2.6.9) Tím pádem je WR W −I = D. Překlápění amplitud kolem střední hodnoty lze tedy dosáhnout aplikací dvou Walshových transformací a fázovým posuvem u jednoho bitu. 2.6.3 Překlopení znaménka amplitudy Abychom dokázali překlopit znaménko u té amplitudy, která přísluší stavu splňujícímu podmínku P(x) = 1, uvažujeme operátor Up působící obecně Up : |x; b → |x; b ⊕ P(x) . Aplikujeme nyní Up na superpozici |ψ = 1 2N/2 2N −1 x=0 αx|x a qubit b ve stavu |b = 1√ 2 (|0 − |1 ). Rozdělme nyní množinu všech x na dvě podmnožiny: X0 = {x|P(x) = 0} a X1 = {x|P(x) = 1} a nechme působit operátor Up na |ψ; b : Up|ψ; b = 1 2(N+1)/2 Up   x∈X0 αx|x,0 + x∈X1 αx|x,0 − x∈X0 αx|x,1 − x∈X1 αx|x,1   77 = 1 2(N+1)/2   x∈X0 αx|x,0 ⊕ 0 + x∈X1 αx|x,0 ⊕ 1 − x∈X0 αx|x,1 ⊕ 0 − x∈X1 αx|x,1 ⊕ 1   = 1 2(N+1)/2   x∈X0 αx|x,0 + x∈X1 αx|x,1 − x∈X0 αx|x,1 − x∈X1 αx|x,0   = 1 2N/2   x∈X0 αx|x − x∈X0 αx|x,   |b . (2.6.10) Vidíme, že stav |b zůstane nezměněn, ale v registru pro argument funkce dojde k překlopení znaménka u těch amplitud, pro které je P(x) = 1. Shrnutí Groverův algoritmus umožňuje efektivně zjistit, který argument x splňuje podmínku P(x) = 1. Pokud tímto argumentem může být libovolné číslo z 2N , museli bychom při klasickém vyhledávání dosazovat do funkčního vztahu průměrně 2N−1 -krát. Groverův algoritmus vyžaduje pouze p 4 2N/2 -násobné dosazení. Tento algoritmus je možné implementovat pomocí kombinací Walshových-Hadamardových transformací, jednobitových fázových posuvů a kvantového obvodu realizujícího testování podmínky P(x). Pojmy k zapamatování • Překlápění amplitud kolem střední hodnoty, • překlápění znaménka amplitudy hledaného stavu, • efektivní vyhledávání. 2.7 Kvantová kryptografie Charles Bennett z IBM, jeden z objevitelů kvantové kryptografie. Ačkoliv by kvantové počítače mohly znamenat určitou hrozbu pro naše soukromí, přináší jiné odvětví kvantové informatiky naopak možnost, jak dosáhnout bezpečného přenosu utajovaných zpráv. Využívá se tu principiální nemožnosti spolehlivě rozlišovat mezi neortogonálními kvantovými stavy a toho, že měření kvantový stav obecně narušuje. Základem kvantové kryptografie je distribuce tajného klíče mezi dvěma vzdálenými partnery tak, aby nikdo jiný klíč nemohl znát. Jakmile dvě strany klíč znají (je to v podstatě náhodná posloupnost nul a jedniček), stačí jej přičíst k binárně kódované zprávě. Takto zašifrovaná zpráva se pak vnějšímu pozorovateli jeví jako zcela náhodná posloupnost nul a jedniček bez jakýchkoliv vnitřních korelací, které by se daly využít k získání nějaké užitečné informace. Průvodce studiem Sdílejí-li dvě strany zcela náhodnou, tajnou posloupnost čísel, může ji vysílající strana přičíst ke své zprávě (zapsané pomocí číslic) a přijímající strana pak 78 od šifrované zprávy odečíst. Tímto způsobem lze předat tajnou zprávu, aniž by z ní někdo nepovolaný mohl zjistit cokoliv rozumného. Problém je, že takovýto klíč můžeme použít pouze jednou: při opakovaném použití lze v předávaných zprávách najít korelace, které by mohly vést k „prosáknutí citlivých informací ven. Možnost předávat tajný klíč jako velmi dlouhou posloupnost náhodných čísel je tedy zcela podstatná. 2.7.1 Distribuce klíče 1 0 ALICE BOB 1 0 1 Obr. 2.14: Schéma kvantové kryptografie, protokol BB84. Jak se takovýto bezpečný klíč distribuuje? Uveďme si pro příklad princip kódu nazvaného BB84 podle jeho autorů Charlese Bennetta a Gilese Brassarda a roku, ve kterém byl publikován. Dvě komunikující strany, zvané tradičně Alice a Bob chtějí sdílet tajný klíč, viz obr. 2.14. Alice posílá Bobovi posloupnost qubitů, které kódují bity ve dvou možných bázích. Pro jednoduchost si představme různé lineární polarizace fotonů: ty mohou být orientovány buď vodorovně (stav | ↔ , což bude označovat hodnotu bitu 0) nebo svisle (stav | , což označuje hodnotu bitu 1). Toto odpovídá „červené bázi „⊕ na obr. 2.14. Alice však využije i kódování v jiné bázi: foton může být lineárně polarizován ve směru +45 stupňů (stav | = 1√ 2 (| + | ↔ ) což bude znamenat bit 0) nebo ve směru −45 stupňů (stav | = 1√ 2 (| − | ↔ ), což znamená hodnotu bitu 1). Toto odpovídá „modré bázi „⊗ na obr. 2.14. Alice náhodně přepíná báze, ve kterých bity kóduje a také hodnoty bitů jsou zcela náhodné. Fotony zaslané Bobovi tedy mohou tvořit posloupnost jako: 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. . . . | | | | ↔ | | | ↔ | | ↔ | | | ↔ . . . 1 1 1 0 0 1 0 0 0 1 1 0 . . . Bob měří polarizace fotonů, které k němu přichází. Protože však neví, v jaké bázi 79 Alice kódovala, náhodně přepíná svůj detektor tak, aby rozlišoval buď mezi horizontální a vertikální polarizací (báze „⊕ ), nebo mezi polarizací ve směru +45 stupňů a −45 stupňů (báze „⊗ ). Připomeňme, že podle zákonů kvantové fyziky je výsledek měření zcela náhodný, pokud Bob měří v jiné bázi, než v jaké Alice kódovala svůj bit. Jestliže tedy Alice poslala foton | a Bob měřil v bázi „⊕ , dostane se stejnou pravděpodobností výsledek | jako | ↔ . Předpokládejme, že posloupnost bází, ve kterých Bob měřil je (ve vztahu k výše uvedeným Aliciným stavům) 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. . . . | | | | ↔ | | | ↔ | | ↔ | | | ↔ . . . 1 1 1 0 0 1 0 0 0 1 1 0 . . . ⊗ ⊗ ⊕ ⊗ ⊕ ⊕ ⊗ ⊕ ⊕ ⊗ ⊗ ⊕ . . . Bobovy naměřené hodnoty pak mohou být například 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. . . . | | | ↔ | | | | | ↔ | ↔ | | | ↔ . . . 1 1 0 1 1 1 1 0 0 0 1 0 . . . Ve druhém, šestém, devátém, jedenáctém a dvanáctém případě měl Bob stejnou bázi jako Alice a musí tedy získat stejnou polarizaci, jakou Alice posílala. V ostatních případech je naproti tomu výsledek Bobova měření zcela náhodný. Bob nyní Alici sdělí, jaké báze pro svá měření používal (ne však výsledky měření) a Alice mu odpoví, ve kterých případech použila stejnou bázi jako Bob. Komunikující strany nyní vědí, kdy mohou očekávat stejné polarizace fotonů, přiřadí jim odpovídající hodnoty bitů a ostatní data zahodí. Pro náš příklad jsou tedy ponechané polarizace ty, které odpovídají fotonům číslo 2, 6, 9, 11 a 12, tedy | , | , | ↔ , | a | ↔ , což představuje hodnoty bitů „1,1,0,1,0 . Tato čísla mohou představovat začátek klíče, sdíleného Alicí a Bobem. Giles Brassard z univerzity v Montrealu, jeden z objevitelů kvantové kryptografie. 2.7.2 Odhalení narušitele Co když se však někdo pokoušel tuto komunikaci odposlechnout a získat nějakou znalost o klíči? Nějaký narušitel (nebo narušitelka - eavesdropper, většinou pojmenovaná Eva) zachytává fotony vyslané Alicí, proměří je a pošle Bobovi. Řekněme, že Eva zná všechny detaily komunikace mezi Alicí a Bobem a ví, že fotony mohou přilétat s některou ze čtyř lineárních polarizací | , | ↔ , | a | . Nemůže však vědět, pro kterou polarizaci se rozhodla Alice v každém jednotlivém případě. Co kdyby se Eva rozhodla detekovat polarizace a náhodně při tom střídat báze? Představme si například Evinu posloupnost bází 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. . . . ⊗ ⊕ ⊗ ⊗ ⊕ ⊕ ⊕ ⊗ ⊗ ⊕ ⊕ ⊕ . . . Její naměřené hodnoty mohou být například 80 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. . . . ⊗ ⊕ ⊗ ⊗ ⊕ ⊕ ⊕ ⊗ ⊗ ⊕ ⊕ ⊕ . . . | | | | | ↔ | | ↔ | | | | | ↔ . . . Ve třetím, šestém, sedmém, osmém, desátém a dvanáctém případě to byla stejná báze jako u Alice a musí tedy získat stejnou polarizaci, ostatní hodnoty jsou náhodné. Fotony s polarizacemi odpovídajícími Eviným datům nyní putují k Bobovi. U šestého a dvanáctého fotonu mají Alice, Eva i Bob stejnou bázi, takže budou mít i stejné naměřené hodnoty polarizace. U druhého, devátého a jedenáctého fotonu jsou báze Alice a Boba stejné, Eva má však báze odlišné. Foton, který přijde Bobovi, bude dávat zcela náhodnou hodnotu polarizace. Ač měli Bob a Alice stejné báze, mohou u daného fotonu pozorovat různé polarizace - a to je podezřelé! Nesoulad mezi vysílanými a naměřenými hodnotami v případech, kdy Alice i Bob měli stejné báze je znakem toho, že narušitel mohl být na drátě. V praxi to znamená, že komunikující strany obětují část dat (řekněme tisíc bitů z celkového počtu deset tisíc) na jejich porovnání. V případě shody mohou Alice a Bob předpokládat, že jejich komunikaci nikdo nenarušil a ty bity, u nichž měli stejnou bázi, použijí jako tajný klíč. To, že by Eva měla štěstí a podařilo se jí odposlechout výměnu dat, aniž by Alice a Bob zjistili nějakou neshodu, je velice malá, v našem případě 2−250 . Pokud se budou některé z kontrolních bitů lišit, znamená to možnost narušení soukromí a Alice a Bob přestanou komunikovat. Poznamenejme, že bezpečnost kvantové kryptografie je mimo jiné založena na nemožnosti klonování kvantových stavů. 2.7.3 Rozvoj kvantové kryptografie Náš příklad ilustroval kryptografický princip na základě polarizovaných fotonů. Existuje řada jiných návrhů, využívajících dalších vlastností světla. Je možné kódovat qubity pomocí interferometrické fáze, pomocí intenzity světelných pulsů, či pomocí prostorových módů, v nichž se pulsy šíří. Výzkum v této oblasti je velice intenzivní a kvantová kryptografie má zatím ze všech odvětví kvantové informatiky nejblíž ke komerčnímu využití. Z řady mezinárodních projektů, které v této oblasti působí, jmenujme například IST QuComm, financovaný z Pátého rámcového programu EU, jehož řešitel prof. Anders Karlsson získal roku 2004 Descartovu cenu (udělovala se v Praze), projekt SECOQC, jehož cílem je vyvinout funkční prototyp komerčně využitelného kryptografického spojení, či projekt COVAQIAL, soustředící se na využití spojitých proměnných pro kvantovou komunikaci. Na posledních dvou jmenovaných projektech se podílí i Univerzita Palackého v Olomouci. K jejím důležitým výsledkům patří i jedno z prvních experimentálně fungujících kryptografických zařízení, které může sloužit i jako systém pro vzájemnou identifikaci komunikujících stran. To je v podstatě kombinace kryptografického systému na distribuci tajného klíče, spojeného s identifikační procedurou pro rozpoznání komunikujících stran [2.7]. Jako qubity tu fungovaly slabé nanosekundové pulsy, cestující optickým vláknem, přičemž informace se nekódovala do jejich polarizace, ale do jejich načasování: Alice svůj puls vláknovým děličem rozdělí na dvě části, z nichž jedna projde delší trasu než druhá, a poté je opět vláknovým děličem spojí (snímek 81 Obr. 2.15: Kvantová kryptografie v Olomouci: Dr. Haderka dolaďuje „Alici , vysílající část kryptografického vláknového interferometru. z laboratoře je na obr. 2.15). Tyto dvě trasy tak fungují jako dvě ramena Machova - Zehnderova interferometru, jejichž vzájemnou fázi můžeme nastavovat. Pokud je puls slabý tak, že obsahuje maximálně jediný foton, bude po průchodu takovýmto interferometrem ve stavu odpovídajícím superpozici dvou možných stavů: „foton vepředu (prošel kratším ramenem) a „foton vzadu (zdržel se v delším rameni). Tyto dvě možnosti spolu s jejich libovolnými superpozicemi představují qubit, s jehož detekcí se musí vypořádat Bob - tedy přijímací stanice. K detekci mu poslouží podobný interferometr, jako má Alice, přičemž nastavení fáze tu hraje roli výběru báze. Obě stanice spojovalo vlákno o celkové délce 0,5 km a zařízení bylo schopné vygenerovat zhruba 6 kilobitů klíče za sekundu. Průvodce studiem V dnešní době již fungují vláknové kryptografy na vzdálenost několika desítek kilometrů se současným rekordem cca 120 km. Důvodem, proč je těžké dostat se dále, jsou absorpční ztráty: narozdíl od běžné optické komunikace je nelze kompenzovat zesilovači, protože ty by do systému vnesly podobný šum, jako případný narušitel. Nicméně i tato vzdálenost může přinést zajímavé aplikace, např. pro komunikaci mezi bankami větší aglomerace. Připomeňme, že první bankovní transakce na bázi kvantové kryptografie proběhla na jaře roku 2004 mezi vídeňskou radnicí a pobočkou místní banky - své schopnosti tu demonstrovala skupina prof. Antona Zeilingera z vídeňské univerzity. Druhou možností je posílat fotony volným prostorem. Za vhodných atmosférických podmínek tak lze komunikovat pomocí fotonů šířících se vzduchem několik desítek kilometrů. Na delší vzdálenosti se však můžeme dostat, pokud foton část své cesty urazí vzduchoprázdným prostorem - pro kryptografickou komunikaci pak bude nutné využít družicového systému. 82 Shrnutí Kvantová kryptografie umožňuje dvěma stranám sdílet tajný klíč. Je založena na výměně kvantových stavů a jejich měření. U protokolu BB84 Alice posílá Bobovi bity kódované pomocí stavů ve dvou různých bázích, Bob na těchto stavech provádí měření. Jak Alice, tak i Bob své báze náhodně mění. Po zaslání všech stavů se strany vzájemně informují o bázích, které použily pro každý stav. Jako součást klíče použijí bity z těch případů, u kterých se jejich báze shodovaly. Pro kontrolu se u části takovýchto bitů informují, zda se skutečně hodnoty bitů shodují. Pokud by se nějaký narušitel pokoušel vysílané stavy zachytit a proměřit, projevilo by se to tím, že některé z kontrolních bitů by se neshodovaly. Pojmy k zapamatování • Protokol BB84, • kryptografický klíč • střídání bází • kvantové měření 2.8 Kvantová provázanost Důležitým aspektem kvantové informatiky je existence tzv. entanglovaných, neboli kvantově provázaných stavů. (Jako kuriozitu je poprvé zmínil Erwin Schrödinger, který je pojmenoval „verschraenkte Zustände , v angličtině se objevují pod názvem „entangled states .) Jedná se o stavy systému, skládajícího se ze dvou či více navzájem odlehlých podsystémů. Vlastnosti těchto podsystémů (představme si například systém dvou fotonů, šířících se opačnými směry) mohou být vzájemně korelované. Pokud je celý systém v kvantově provázaném stavu, budou vzájemně korelované různé páry veličin, a to tak, že korelací je více, než by nám dovolovaly úvahy na základě klasické fyziky (přesněji řečeno, korelace mezi měřeními na dvou kvantových podsystémech mohou narušovat tzv. Bellovy nerovnosti). Pro ilustraci uvažujme například dva fotony A a B v kvantovém stavu |ψ = 1 √ 2 (| A| ↔ B − | ↔ A| B) . (2.8.1) Pokud bychom měřili polarizaci jediného fotonu, dostali bychom zcela náhodnou hodnotu. Polarizace každého fotonu je ovšem vždy přesně opačná, než má jeho partner: pokud u fotonu A naměříme vertikální polarizaci, bude mít foton B s jistotou horizontální polarizaci a naopak. A neplatí to jen pro měření v bázi vertikální a horizontální polarizace. Zkusme pro ilustraci převést stav |ψ do báze kruhových polarizací. Ze vztahů (2.2.1) a (2.2.2) zjistíme, že stavy | ↔ a | mohou být zapsány ve tvaru | = 1 √ 2 (|R − i|L ) , (2.8.2) | ↔ = 1 √ 2 (−i|R + |L ) . (2.8.3) 83 Stav |ψ tedy můžeme psát jako |ψ = 1 2 √ 2 [(|R A − i|L A) (−i|R B + |L B) − (−i|R A + |L A) (|R B − i|L B)] = 1 √ 2 (|R A|L B − |L A|R B) . (2.8.4) Znamená to tedy, že pokud je foton A pravotočivě polarizovaný, bude foton B levotočivý a naopak. Kdykoliv ale zjistíme, že foton A měl horizontální polarizaci, bude foton B vertikálně polarizovaný. Jak je to možné, když přece lineární polarizace fotonu znamená, že se foton bude chovat zcela náhodně při průchodu analyzátorem detekujícím kruhovou polarizaci? Dá se uvažovat, že si každý foton s sebou nese jakousi ukrytou informaci, která rozhoduje o tom, jak na který analyzátor reagovat a je to jen naše neznalost, díky které považujeme chování fotonu za „náhodné ? John S. Bell roku 1964 ukázal, že pokud by takovéto „lokální skryté proměnné existovaly, musely by statistiky měření různých veličin splňovat určité nerovnosti. Kvantová mechanika však předpovídá možnost narušení Bellových nerovností a řada experimentů jí dává zapravdu. Plyne z toho, že existenci lokálních skrytých proměnných musíme nejspíš ze svých úvah o přírodě vyloučit. Stav |ψ , nazývaný singlet, je typickým příkladem „entanglovaného , tedy kvantově provázaného stavu. Je ale nutné si uvědomit, že ne všechny korelace odpovídají kvantové provázanosti. Představme si, že nám nějaký zdroj posílá dvojice fotonů, které budou s 50% pravděpodobností ve stavu | A| ↔ B a s 50% pravděpodobností ve stavu | ↔ A| B. V kvantovém formalismu by to odpovídalo matici hustoty = 1 2 | A | ⊗ | ↔ B ↔ | + 1 2 | ↔ A ↔ | ⊗ | B |. (2.8.5) I v tomto případě jsou lineární polarizace obou fotonů korelované: pokud u fotonu A zjistíme vertikální polarizaci, bude foton B s jistotou polarizován horizontálně a naopak. Narozdíl od předchozího případu se stavem |ψ , u stavu však nebudeme pozorovat žádné korelace mezi kruhovými polarizacemi obou fotonů. Matice hustoty z rovnice (2.8.5) neodpovídá kvantově provázanému stavu. Přesněji se kvantová provázanost definuje takto: jakýkoliv kvantový stav dvou podsystémů A a B lze popsat maticí hustoty AB. Pokud je možné tuto matici hustoty zapsat jako součet AB = k pk (k) A (k) B , (2.8.6) kde 0 < pk 5 1 a (k) A jsou matice hustoty podsystému A a (k) B jsou matice hustoty podsystému B, nazývá se stav AB faktorizovatelný. Pokud neexistuje žádná možnost, jak matici hustoty AB zapsat ve tvaru (2.8.6), je příslušný stav kvantově provázaný. V kvantové informatice se často pracuje s kvantově provázanými stavy dvou qubitů. Je proto velmi užitečné zavést pro dvoubitové systémy bázi, jejíž všechny 84 stavy jsou kvantově provázané. Jedná se o takzvanou Bellovu bázi: |Φ+ = 1 √ 2 (|00 + |11 ) , (2.8.7) |Φ− = 1 √ 2 (|00 − |11 ) , (2.8.8) |Ψ+ = 1 √ 2 (|01 + |10 ) , (2.8.9) |Ψ− = 1 √ 2 (|01 − |10 ) . (2.8.10) Jak se můžeme přesvědčit, singletní stav z našeho příkladu odpovídá Bellovu stavu |Ψ− . Průvodce studiem Příroda nám někdy poskytuje zdroj entanglovaných stavů. Příkladem může být frekvenční konverze fotonů v nelineárních krystalech: z jednoho fotonu o frekvenci ω se mohou stát dva fotony s nižšími frekvencemi ω1 a ω2, pro něž platí ω1 + ω2 = ω. Tyto dva nové fotony jsou v kvantově provázaném stavu. Často ale bývá problém, jak dva původně nezávislé fotony (například přilétající ze dvou různých zdrojů) přivést do vzájemně kvantově provázaného stavu. Korelace u kvantově provázaných stavů svědčí o tom, že příroda se chová podstatně jinak, než by odpovídalo naší běžné intuici, z určitého pohledu to vypadá, jakoby výsledek měření na jednom podsystému okamžitě (nadsvětelnou rychlostí) ovlivnil chování jiných částí systému. Zároveň však toto „nadsvětelné ovlivnění má takový charakter, že ho nelze bezprostředně využít k předání informace: informace sama se může šířit maximálně rychlostí světla ve vakuu. Tyto vlastnosti kvantových systémů daly příčinu k mnoha filosofickým diskusím o tom, co vlastně znamená kvantový stav, zda odpovídá nějaké realitě nebo spíš naší subjektivní znalosti a co vlastně „za tím vším je. Ale v poslední době se spíše lidé snaží zjistit, jak by se tyto podivuhodné vlastnosti přírody daly prakticky využít. Například Artur Ekert roku 1991 zjistil, že provázané stavy by byly vynikajícím prostředkem pro kryptografii (jeho protokol E91 je odlišný od výše uvedeného protokolu BB84). Kromě toho budou kvantově provázané stavy důležité pro chod kvantových počítačů a již dnes jsou základem dalších zajímavých procesů: kvantové teleportace, hustého kódování a dalších. Shrnutí Kvantová provázanost (entanglement) je vlastnost stavů systémů, které jsou složené ze dvou (nebo více) podsystémů a měření na nich vykazují určitý typ korelací. Kvantově provázané stavy jsou takové, které nejsou faktorizovatelné—jejich matice hustoty nelze přepsat do tvaru konvexní sumy součinů matic hustoty jednotlivých podsystémů. Kvantově provázané stavy mohou narušovat tzv. Bellovy nerovnosti, 85 zatímco faktorizovatelné stavy nikoliv. U stavů dvou qubitů se často používají kvantově provázané stavy tvořící tzv. Bellovu bázi. Pojmy k zapamatování • Faktorizovatelný stav, • Bellova báze, • singlet, • kvantová provázanost. Cvičení 1. Dokažte, že Bellovy stavy tvoří ortonormální systém stavů dvou qubitů. 2. Najděte unitární matici, která permutuje Bellovy stavy |Φ+ → |Φ− → |Ψ+ → |Ψ− → |Φ+ 2.9 Kvantová teleportace Jednou z aplikací kvantové provázanosti je kvantová teleportace. Předpokládejme, že Alice má nějaký systém—pro jednoduchost předpokládejme qubit—v neznámém kvantovém stavu. Sama tento stav nepotřebuje, ale ráda by ho neporušený předala Bobovi, který je na vzdáleném místě a nelze mu qubit předat fyzicky. Alice si však s Bobem může telefonovat. Jak tento problém vyřešit? Jako první myšlenka nás může napadnout, že Alice jednoduše změří svůj qubit, výsledek zatelefonuje Bobovi a ten si svůj qubit připraví v tomto stavu. Ovšem něco takového nám kvantová mechanika nedovolí. Alice může měřit svůj qubit jen v jediné bázi a pokud takové měření provede, nutně musí přijít o veškerou informaci, která mohla být zakódovaná v jiné bázi. Předpokládejme například, že Alice má svůj qubit ve stavu 1 2 |0 + i √ 3 2 |1 . Pokud by Alice svůj qubit změřila v bázi |0 , |1 , dostala by buď jedničku nebo (s menší pravděpodobností) nulu, ale nemohla by zjistit žádnou informaci o fázi mezi těmito dvěma stavy. Předpokládejme ale navíc, že Alice s Bobem již sdílejí nějaký kvantově provázaný stav, například singlet |Ψ− . V tomto případě může Alice Bobovi svůj stav teleportovat. Pro popis teleportace nám poslouží obr. 2.16. Pro odlišení si označme Alicin vstupní systém 0 a systémy, které nesou provázaný stav jako 1 (tento systém bude mít Alice) a 2 (bude mít Bob). Předpokládejme, že Alicin qubit je ve stavu α|0 0 + β|1 0. Na počátku je tedy celý tříčásticový systém ve stavu |ψ 012 = (α|0 0 + β|1 0) 1 √ 2 (|0 1|1 2 − |1 1|0 2) . (2.9.1) Teleportace pak probíhá následujícím způsobem: 1. Alice provede měření na kombinovaném systému 0 a 1 v Bellově bázi. Tím získá výsledek, který je buď Φ+ , Φ− , Ψ+ , nebo Ψ− . Její vstupní stav je tímto měřením zlikvidován a navíc nemůže o tomto stavu již nic zjistit ani z výsledku měření: v každém z Bellových stavů je totiž stejně zastoupen jak stav |0 0 tak i stav |1 0. 86 |ψin |ψout |Ψ− Měření v bázi Φ+ , Φ− , Ψ+ , Ψ− Manipulace Alice Bob ENT 0 1 2 Obr. 2.16: Kvantová teleportace jednoho qubitu. Alice má qubit v nějakém neznámém stavu |ψin a chce ho předat Bobovi, aniž by ho mohla fyzicky poslat. Může k tomu využít sdílený kvantově provázaný stav. 2. Alice výsledek měření zatelefonuje Bobovi (v podstatě mu předá jen dva bity informace). 3. Bob podle výsledku měření provede jednu ze čtyř možných transformací na svém systému 2. Poté je systém 2 ve stejném kvantovém stavu, v jakém byl systém 0 na počátku. Proč tomu tak je zjistíme, když budeme působit na stav (2.9.1) projekčními operátory Bellových stavů, což odpovídá měření v Bellově bázi: |Φ+ 01 Φ+ |ψ 012 = |Φ+ 01(α|1 2 − β|0 2), (2.9.2) |Φ− 01 Φ− |ψ 012 = |Φ+ 01(α|1 2 + β|0 2), (2.9.3) |Ψ+ 01 Ψ+ |ψ 012 = |Ψ+ 01(β|1 2 − α|0 2), (2.9.4) |Ψ− 01 Ψ− |ψ 012 = |Ψ− 01(β|1 2 + α|0 2). (2.9.5) Znamená to tedy, že pokud Alice zjistí stav Ψ− , nebude Bob dělat se svým qubitem nic a jeho stav bude roven vstupnímu stavu na Alicině qubitu. Pokud Alice zjistí stav Ψ+ , bude muset Bob fázově posunout amplitudu u stavu |0 2 o p. Pokud Alice zjistí některý ze stavů Φ, bude muset Bob rotovat svůj qubit mezi bázovými stavy |0 2 a |1 2 a případně i provést fázový posuv. Každopádně po některé z těchto jednoduchých operací bude výstupní stav Bobova qubitu roven vstupnímu stavu qubitu u Alice, a to aniž se kdokoliv mohl tuto hodnotu stavu dozvědět. Shrnutí Kvantová teleportace je předávání kvantového stavu z jednoho systému na druhý— vzdálený. Využívá se při ní sdíleného kvantově provázaného stavu a klasické komunikace. Alice provede měření v Bellově bázi na vstupním stavu a její částí sdíleného provázaného stavu. Tím se ovšem vstupní stav zničí. Výsledek měření předá Alice Bobovi a ten na jeho základě provede transformaci na své části sdíleného provázaného stavu. Tím původně zničený vstupní stav nově vznikne na Bobově podsystému. Kontrolní otázky 1. Jaké hradlo by měl použít Bob na svůj qubit, pokud Alice naměřila stav |Ψ+ 01? 87 2.10 Závěr Kvantová informatika si získává místo jak ve fyzikálních, tak v matematických oblastech. Po matematické stránce jsou nejzajímavější možné nové algoritmy a protokoly zpracování kvantové informace. Z fyzikálního hlediska je pak nejdůležitějším problémem jejich konkrétní implementace. Tento text by měl sloužit jako základní orientace pro ty, kdo by se chtěli některým otázkám kvantové informatiky věnovat. Případný zájemce najde další informace buď v původní literatuře (např. [2.4, 2.6]), v přehledových článcích (k vynikajícím patří Braunsteinův [2.5]) nebo v monografiích (např.[2.3, 2.2]). Literatura ke kapitole 2 [2.1] Feynman R.P., Leighton R.B., Sands M.: Feynmanovy přednášky z fyziky. Fragment, Havlíčkův Brod 2002. [2.2] Dušek M.: Koncepční otázky kvantové teorie. Vydavatelství Univerzity Palackého, Olomouc 2002. [2.3] Gruska J.: Quantum computing. McGraw-Hill, London 1999. [2.4] Shor P.W.: Polynomial-time algorithms for prime factorization and discrete logarithms on a quantum computer. http://www.arxiv.org/abs/quant-ph/9508027 [2.5] Braunstein S.L.: Quantum computation. http://www-users.cs.york.ac.uk/ ~schmuel/comp/comp.html [2.6] Grover, L.K.: Quantum mechanics helps in searching for a needle in a haystack. Phys. Rev. Letters 79, 325-328, 1997. http://www.arxiv.org/abs/quant-ph/ 9706033 [2.7] Dušek, M., Haderka, O., Hendrych, M., Myška, R.: Quantum identification system. Physical Review A 60, p. 149, 1999. http://optics.upol.cz/dusek/clanky/ pra60-149.pdf 88 Přílohy 89 Příloha A Fyzikální konstanty a jednotky použité v textu V přehledné tabulce uvádíme fyzikální konstanty a jednotky použité v textu spolu s jejich hodnotami v soustavě SI. Boltzmannova konstanta kB = 1,380 662 ·10−23 J·K−1 hmotnost Slunce M = 1,989 ·1030 kg gravitační konstanta G = 6,672 0 ·10−11 N·m2 ·kg−2 parsek 1 pc = 3,085 977 56·1016 m Planckova konstanta h = 6,626 176·10−34 J·s rychlost světla ve vakuu c = 2,997 924 580·108 m·s−1 světelný rok 1 ly = 9,460 528 3·1015 m 90 Příloha B RSA kryptografický kód RSA kód je nazván podle svých objevitelů, kterými jsou Ronald Rivest, Adi Shamir a Leonard Adleman. Funguje následujícím způsobem. Předpokládejme, že dva účastníci—Alice a Bob spolu chtějí bezpečným způsobem komunikovat. Alice bude přijímat od Boba zprávy, které nikdo kromě Alice nedokáže rozluštit. 1. Alice zvolí dvě velká prvočísla p a q a vynásobí je, takže dostane číslo N = pq. Pro ilustraci tu ale pracujme s malými prvočísly, zvolme si třeba p = 17 a q = 11, takže N = 17 × 11 = 187. 2. Alice dále vybere číslo e, které je nesoudělné jak s (p − 1) tak s (q − 1). V našem příkladu tedy hledáme číslo nesoudělné s 16 a 10. Zvolme e = 7. 3. Alice dále najde číslo d, které splňuje podmínku e × d = 1(mod (p − 1)(q − 1)). (2.0.1) V našem případě tedy hledáme číslo d splňující 7d = 1(mod 16 × 10), tedy 7d = = 1(mod 160). Jak si můžeme ověřit, hledané číslo je d = 23. 4. Alice může komukoliv sdělit čísla N a e, která budou sloužit jako její veřejný klíč, ale čísla p, q a d si ponechá v tajnosti.Čísla N a d Alici slouží jako její soukromý klíč. 5. Pro zakódování zprávy musí Bob nejprve konvertovat posílaný text na číslo— například pomocí ASCII kódu. Zprávě (zatím nezašifrované) tak odpovídá číslo M. Šifrování se děje předpisem C = Me (mod N). (2.0.2) Šifrovanou zprávou je tedy číslo C. Pro ilustraci předpokládejme, že Bob chce Alici poslat písmeno X, které v ASCII kódu odpovídá číslu 88. Tedy M = 88 a C = 887 (mod 187) = 894 432(mod 187) = 11. (2.0.3) 6. Bob pošle šifrovanou zprávu C Alici. 7. Alice může zprávu rozšifrovat pomocí svého soukromého klíče vztahem M = Cd (mod N). (2.0.4) V našem případě to tedy znamená M = 1123 (mod 187) = 88. (2.0.5) 91 Veřejným klíčem může kdokoliv zašifrovat zprávu, kterou chce poslat Alici. Zašifrovanou zprávu ale může rozluštit jen ten, kdo má soukromý klíč. Aby někdo mohl ze znalosti veřejného klíče získat soukromý klíč, musel by dokázat faktorizovat číslo N, tedy najít jeho činitele p a q. To je ovšem s rostoucí velikostí N výpočetně velice náročné, na klasických počítačích doba výpočtu roste exponenciálně s počtem číslic čísla N. Kvantový počítač by však dokázal číslo N faktorizovat v čase, který roste jen polynomicky s počtem číslic N. 92 Rejstřík Albert Einstein 17 algoritmus – Euklidův 73 – Groverův 51, 75 – Shorův 51, 70 báze 56–58 – Bellova 85, 86 Bell John 84 Bennett Charles 78 Big Bang viz velký třesk Big Crunch viz velký křach Brassard Giles 80 comoving souřadnice 13 čas Hubbleův 23 decelerační parametr viz parametr de- celerační délka Hubbleova 31 Einsteinovy rovnice viz rovnice Ein- steinovy Einsteinův statický vesmír viz vesmír Einsteinův statický Ekert Artur 85 entanglement 83 entropie 61 expanzní faktor 13 faktorizace 73 fázový posuv 64, 77 – řízený 65, 69 Feynman Richard 49 fotometrická vzdálenost viz vzdálenost fotometrická Friedmann Alexandr Alexandrovič 14 Friedmannova rovnice viz rovnice Friedmannova Grover Lov 75 hradlo – CNOT 64–66 – dvoubitové 64 – Hadamardovo 63, 65, 69 – jednobitové 62 – kvantové 60 – sčítací 67 – SWAP 66 – Toffoliho 66 Hubble Edwin Powell 23 Hubbleova délka viz délka Hubbleova Hubbleův čas viz čas Hubbleův hustota kritická 29 hustotní parametr viz parametr hus- totní hvězdná velikost – absolutní 41 – relativní 41 inflace 21, 30 klíč 53, 78, 82 klonování 54, 61, 81 konstanta – Hubbleova 23 – kosmologická 12 Koperník Mikoláš 6 kosmologická konstanta viz konstanta kosmologická kritická hustota viz hustota kritická kryptografie – kvantová 53, 78, 85 kvantová provázanost 83 magnituda – absolutní 41 – relativní 41 93 měření – kvantové 58, 76, 78, 80, 86 modul vzdálenosti 41 nulování 61 Olbers Heinrich Wilhelm 7 Olbersův paradox viz paradox Olber- sův paradox Olbersův 7 parametr – decelerační 30 – Hubbleův 23 – hustotní 29 polarizace 56, 58, 79, 83 prach 18 princip – Koperníkův 6 – Landauerův 61 – superpozice 48, 50 problém plochosti vesmíru 30 qubit 52, 56, 58, 86 quintessence 20 reverzibilita 60 rotace qubitu 64 rovnice – Einsteinovy 17 – Friedmannova 12–16 – stavová 16, 18 – – vakua 20 – – záření 19 RSA 73, 91 rudý posuv 26 řízená negace 64 Shor Peter 70 Schrödinger Erwin 83 spin 48 standardní kosmologický model 6 standardní svíčky 39 stavová rovnice viz rovnice stavová superpozice 48, 71, 75 teleportace 53, 86 transformace – CNOT 64 – fázového posuvu 64 – Hadamardova 63, 65 – Hadamardova-Walshova 63 – kvantová Fourierova 68, 71 – rotace qubitu 64 – Walshova-Hadamardova 68, 71, 75, 77 úloha – Deutschova 50 velký křach 32 velký třesk 32 vesmír – Einsteinův statický 33 – otevřený 14 – plochý 14 – uzavřený 14 vzdálenost fotometrická 39 zákon – druhý termodynamický 61 – Hubbleův 23 – Moorův 47 94