7 1 Skalární vlna a její matematický popis 1.1 Vlna a vlnová rovnice 1.2 Rovinné vlny 1.3 Kulové vlny 1.4 Harmonické vlny 1.5 Komplexní notace harmonických vln 1.6 Intenzita vlnění a výpočet intenzity harmonických vln 1.7 Helmholtzova rovnice 1.8 Poznámka o polárním tvaru řešení Helmholtzovy rovnice 1.9 Paraxiální vlna a paraxiální Helmholtzova rovnice 1.10 Fresnelova aproximace kulové harmonické vlny 1.11 Důležitá věta o řešení Helmholtzovy rovnice Pojmů vlna a vlnění se v různých oblastech fyziky a techniky používá k označení velmi odlišných jevů. Připomeňme jen pojmy jako akustická vlna, povrchová vlna, rázová vlna, elektromagnetická vlna, gravitační vlna. Je proto dosti obtížné najít něco, co je všem těmto jevům společné, a je asi nemožné nabídnout nějakou definici, která by zahrnovala všechny jevy označované jako vlna nebo vlnění. Někdy se setkáváme s tím, že se vlněním označuje proces, při němž se prostorem šíří energie, aniž se současně přenáší hmotnost. Toto vymezení skutečně obstojí pro velkou oblast fyziky. Těžko však obstojí, chceme-li používat pojmů jako de Broglieova vlna, elektronová vlna apod. V elektronové a neutronové difraktografii však tyto pojmy zdomácněly do té míry, že si lze jen obtížně představit, že by mohly být něčím nahrazeny. Učebnice fyziky většinou postupují tak, že uvedou konkrétní příklad nějakého vlnění a ukáží, že vyhovuje jisté rovnici, které se říká vlnová rovnice. Potom, kdykoliv se setkají s procesem, který vyhovuje obdobné rovnici, nazývají tento proces vlněním. Tohoto postupu se přidržíme i my. 1.1 Vlna a vlnová rovnice Vystačíme s nej jednodušším tvarem vlnové rovnice, který charakterizuje vlnění v bodech homogenního, izotropního, nedisperzního a neabsorbujícího prostředí neobsahujících zdroje vlnění. Jak je známo (viz např. [1], str. 301), jde v tomto případě o lineární parciální diferenciální rovnici druhého řádu Funkci \P(r, ť) se říká vlnová funkce a její fyzikální význam závisí na tom, o jaké vlnění jde. U akustického vlnění je to výchylka z rovnovážné polohy, v případě elektromagnetického vlnění je to elektrická intenzita nebo magnetická indukce nebo některá složka těchto vektorů atd. Ve skalární teorii difrakce se pro vlnovou funkci \P často používá názvu rozruch (disturbance). r značí polohový vektor a t čas. Laplaceův operátor V2 je diferenciální operátor (viz např. [1], str. 144), jenž obsahuje nejvýše druhé derivace podle prostorových souřadnic a jeho konkrétní tvar závisí na dimenzi prostoru a na použité soustavě souřadnic. V £3 má v kartézských souřadnicích x\,X2,x$ tvar y2_ #1 ,2) dx\ dx\ dx\ Protože na levé straně vlnové rovnice (1) jsou druhé derivace podle souřadnic a na pravé straně druhá derivace podle času, je zřejmé, že veličina v má rozměr rychlosti a má význam tzv. fázové rychlosti. (V našem případě lze mluvit o konstantě v, neboť prostředí, v němž se vlnění šíří, je homogenní, izotropní a nedisperzní.) Plocha \P(f, t o) = konst., v jejíchž bodech má v určitém okamžiku tQ vlnová funkce touž hodnotu, se nazývá vlnoplocha. (Poznámka k pojmu vlnoplocha: V teorii vlnění jsou velmi významné tzv. harmonické vlny. Počínaje odst. 1.4 se budeme zabývat pouze tímto typem vln. U harmonických vln se však pojmu vlnoplocha používá v jiném smyslu, než jak jsme jej právě zavedli - viz odst. 1.8.) 8 1 SKALÁRNÍ VLNA A JEJÍ MATEMATICKÝ POPIS V závislosti na počátečních a okrajových podmínkách existuje velmi mnoho různých řešení vlnové rovnice. V následujících dvou odstavcích odvodíme dva typy řešení: Tzv. rovinné vlny, když vlnoplochy jsou roviny a kulové vlny, když vlnoplochy jsou kulové plochy. 1.2 Rovinné vlny Uvažujme o vlnění, jehož vlnoplochy jsou roviny kolmé na jednotkový vektor ň (rii, ri2, n^). Vlnoplochy mají tedy rovnici r ■ n = konst. (srov. obr. 1) a vlnová funkce takového vlnění musí mít tvar V{r,t) = V(f ■ ň,t). (1) Takové vlnění má jednorozměrnou povahu. Formálně to nahlédneme, zavedeme-li novou soustavu souřadnic O, yi , y2 , V3 s osou yi ve směru vektoru ň. Rovnice vlnoplochy pak je r-ň = yi (2) a hodnota souřadnice y\ má význam vzdálenosti vlnoplochy od počátku O (viz obr. 1). Derivováním vlnové funkce \P(r- ň,t) dostaneme Obrázek 1: Soustavy souřadnic při odvozování řešení vlnové rovnice ve tvaru rovinné vlny. Podobně Takže dW _ d(f-n) _ dx\ d{f ■ n) dx\ dy\ ' d2V _ d2V 2 dx\ dy\ 1 dx2, dy\ 2' dx\ dy\ 3 „2 T d2W V =- dy\ a vlnová rovnice 1.1(1) přechází v jednorozměrnou vlnovou rovnici (říká se jí též rovnice kmitů struny) 1.3 Kulové vlny 9 dy\ v2 Abychom nalezli její řešení, zavedeme nové proměnné dt2 (3) zi=yi- vt, z2 = yi+ vt. (4) (Nepřehlédněme, že touto substitucí se do řešení vnáší vztah mezi prostorovými souřadnicemi a časem.) Derivováním vyjádříme druhé derivace v rovnici (3) derivacemi podle nových proměnných: dyi d2>P d& dzx d& dz2 d& d& --- H---- =--1--, dzi dyi dz2 dy\ dz\ dz2' d2& d2V dz\ dz\dz2 d2W dž*' ~dt d2V Dosazením do rovnice (3) dostaneme d& dzx d& dz2 dz\ dt dz2 dt 'd2V d2V dz -2 dzi d2V dz\dz2 d z2 d>P\ dz2 d2W dz\dz2 = 0. Obecným řešením této rovnice je součet 9 = #i(zi) + V2{z2) = #1(2/1 - vt) + V2(Vi + vt), (5) (6) kde SP\ a ty2 jsou libovolné funkce, mající ovšem požadované derivace. S pomocí vztahu (2) se vrátíme k původním proměnným r ■ ň, í a dostaneme \P(r ■ ň, t) = \Pi(r ■ n — vt) + ^2(f ■ ň + vt). (7) I když SP\ a ty2 jsou libovolné funkce, můžeme určit směr šíření vln, které tyto funkce představují. Ukážeme, že funkce í\(r • ň — vt) představuje vlnu, jež se šíří rychlostí v ve směru ň a podobně ^2(f ■ ň + vt) vlnu, šířící se rychlostí v ve směru —ň. Zvolme nějaký časový interval Aí. Nahradíme-li ve funkci SP\{f ■ ň — vt) proměnné r, t výrazy r + vňAt, t + Aí, nezmění se hodnota funkce SP\. Jinými slovy, za čas Aí se vlna posune o vňAt, tj. šíří se rychlostí vň. Podobně vlnová funkce ty2(r ■ ň-\-vť) nezmění svou hodnotu, když se proměnné r, t nahradí výrazy r — vňAt, t + Aí. Funkce tedy představuje vlnu, jež se šíří rychlostí —vň. 1.3 Kulové vlny Má-li mít vlnění kulové vlnoplochy se středem v počátku souřadnic, musí mít vlnová funkce tvar \P(f, t) = SP(r,ť), kde b -j^ I Ju 2 I Ju g ■ Laplaceův operátor má v tomto případě (v £3) tvar d2 + 2d_ V2 = (1) (2) Lze se o tom přesvědčit derivováním: 10 1 skalární vlna a její matematicky popis dW ďP dr xx d>P dxi dr dxi r dr ' d2V 1 dW x\ dW x\ d2W dx\ r dr r3 dr + r2 dr2 Podobně pro derivace podle x2 a x3. Po dosazení do 1.1(2) pak použitím (1) dostaneme vlnovou rovnici pro kulové vlny ve tvaru d2W 2dW 1 d2W (3) dr2 r dr v2 dt2 Po vynásobení obou stran této rovnice faktorem r =/= 0 nahlédneme, že ji lze přepsat do tvaru d2{rW) _ 1 d2{rW) dr2 v2 dt2 Dostáváme opět rovnici kmitů struny (srov. 1.2(3)), nyní však pro funkci r!řr(r, t). Takže řešení rovnice (3) má tvar (srov. 1.2(6)) r\P(r, t) = SP1(r - vt) + \P2(r + vt), tj- \P(r,t) = -\P1(r-vt) + -\P2(r + vt), (4) kde W\ a \P2 jsou libovolné funkce mající požadované derivace. Funkce -\P(r — vt) představuje tzv. divergentní kulovou vlnu, tj. vlnu rozbíhající se z počátku rychlostí v. Naproti tomu funkce ^^(r + vt) je konvergentní kulová vlna, tj. vlna sbíhající se rychlostí v do počátku. 1.4 Harmonické vlny Rovinné a kulové vlny ve tvaru 1.2(7) a 1.3(4) jsou velmi obecné. Mohou představovat vlnu tvořenou jen nějakým pulsem, nebo ustálené vlnění nebo něco mezi tím. Nelze tedy s těmito výrazy spojovat nějakou periodicitu, např. vlnovou délku nebo frekvenci. K tomu je třeba vlnění dále specifikovat. Pro teorii vlnění jsou důležité tzv. harmonické vlny. Jsou definovány tím, že vlnová funkce se ve všech bodech r mění s časem t harmonicky, tj. jako funkce coscjí nebo sincjí nebo jejich lineární kombinace. Tím do teorie vlnění vstupují veličiny vztahující se k časové periodičnosti vlnění: Nejmenší perioda T funkcí coscijí, sincijí zřejmě je T = 2tt/lu, takže úhlová frekvence lu a frekvence v jsou určeny vztahy lu=— v=- (1) Důležitost harmonických vln je založena na tom, že řešení vlnové rovnice metodou separace proměnných vede na harmonické vlny a vyjadřuje obecné řešení vlnové rovnice superpozicí harmonických vln (viz odst. 1.6 v dalším textu). Význam představy harmonických vln nesnižuje ani to, že ve skutečnosti neexistují - předpokládají totiž neohraničené trvání a jsou pouze matematickým objektem. Podstatné je, že každou existující vlnu lze vyjádřit jako superpozici harmonických vln (tzv. vlnové klubko), což má nesmírný teoretický i praktický význam. Požadavkem harmoničnosti jsou do značné míry určeny funkce charakterizující harmonické rovinné a harmonické kulové vlny. Konkrétně, rovinnou harmonickou vlnu popisuje podle 1.2(7) funkce >P(r ■ ň =F vt) = a cos r • n u t =f- (2) v níž a a. a jsou konstanty, znaménko mínus platí pro vlnu šířící se ve směru ň a znaménko plus pro vlnu šířící se ve směru —ň. Hodnota vlnové funkce (2) v bodech r a r + An je zřejmě táž, pokud 1.4 Harmonické vlny 11 Veličina A se nazývá vlnová délka. Vlnovou délkou tedy rozumíme nejmenší prostorovou periodu rovinné vlny ve směru šíření vlnění. To je podstatné, neboť to vysvětluje, proč veličina, jež ve vztahu k prostorovým proměnným odpovídá úhlové frekvenci lu, je vektor. Nazývá se vlnový vektor a z výrazu (2) je zřejmé, že je účelné ji definovat výrazem k = —ň. (4) v Dosazením do vztahu (4) za lu/v ze vztahu (3) získáme pro vlnový vektor k výraz k = ^n. (5) A Časová perioda T a vlnová délka A spolu souvisejí vztahem, který dostaneme dosazením do vztahu (3) za lu ze vztahu (1): A = vT. (6) Můžeme tedy s použitím vlnového vektoru k přepsat rovinnou harmonickou vlnu (2) do tvaru W(f, ť) = a cos(ciJÍ =p k ■ r + a). (7) Podobně lze postupovat při odvozování výrazů pro kulové harmonické vlny. Musí mít podle 1.3(4) tvar a r r v W{r,ť) = - cos w(í=f- \+a , (8) v němž a a a jsou konstanty, znaménko mínus platí pro divergentní kulovou harmonickou vlnu a znaménko plus pro konvergentní kulovou harmonickou vlnu. Bez ohledu na faktor 1/r, pouze na základě periodicity funkce kosinus, se zavádí vlnová délka výrazem (3) a platí pro ni vztah (6). Vlnový vektor k se u kulových harmonických vln nezavádí, pouze jeho velikost w 2"7t . . k = - = —. 9 v A (Bylo by ovšem možné i u kulových harmonických vln definovat k = ale nebývá to zvykem a je to zbytečné.) Kulové harmonické vlny tedy charakterizuje výraz \P(r, ť) = - cos(wí =F kr + a). (10) r Konstantě a ve výrazech (2), (7), (8) a (10) se říká amplituda. Je však podstatný rozdíl ve fyzikálním významu této konstanty u rovinných vln a u kulových vln. Představuje-li vlnová funkce \P nějakou fyzikální veličinu (např. u mechanických vln je to výchylka z rovnovážné polohy), představuje amplituda a u rovinných vln (2) resp. (7) touž veličinu a má samozřejmě týž rozměr (neboť funkce kosinus je bezrozměrná). Naproti tomu v případě kulových vln (8) a (10) má amplituda a rozměr oné fyzikální veličiny, o kterou při vlnění jde, násobený délkou (neboť funkce cos(- • -)/r má rozměr délky na mínus prvou). Proto se u kulových vln amplituda a interpretuje jako amplituda harmonických kmitů v jednotkové vzdálenosti od zdroje kulových vln. To je nepraktické, často se na to zapomíná a může to působit komplikace (např. když se opakovaně používá Huygensova-Fresnelova principu při vlnovém popisu optického zobrazení). Chceme-li se tomu vyhnout, musíme charakterizovat kulové harmonické vlny výrazem \P(r, ť) = — cos(wí =F kr + a). (11) kr v němž funkce cos(- • -)/(kr) je bezrozměrná a amplituda a má - stejně jako u rovinné vlny - rozměr fyzikální veličiny, kterou vlnění představuje. Kromě toho má ve výrazu (11) amplituda a význam amplitudy harmonických kmitů ve vzdálenosti od zdroje, pro níž je kr = 1, tj. podle (9) r = X/(2ir), tedy ve vzdálenosti vztažené k samotnému vlnění, a ne k jednotce délky, jež je ovšem zvolena nezávisle na vlnění. 12 1 SKALÁRNÍ VLNA A JEJÍ MATEMATICKÝ POPIS 1.5 Komplexní notace harmonických vln K popisu periodických dějů se v teorii vlnění užívá místo funkcí kosinus resp. sinus komplexní exponenciální funkce. Tato tzv. komplexní notace má výhodu v tom, že některé matematické operace (jmenovitě násobení, derivování a integrování) se provádějí snadněji a zápis výpočtů a odvozování je mnohem přehlednější. Pokud je to možné (omezení jsou naznačena v závěru tohoto odstavce) - nahrazujeme tedy reálné vlnové funkce komplexními, s nimi provádíme potřebné výpočty a teprve nakonec se vracíme opět k reálným funkcím. Komplexní notace se používá zhruba od konce 19. století a během doby vznikla i terminologie, která z této notace vychází. Bylo by sice i dnes možné se obejít bez komplexní notace, bylo by to však těžkopádným anachronismem. Určitou nevýhodou, kterou sebou nesou výhody komplexní notace, je, že komplexní notace výpočty poněkud formalizuje a nenutí člověka ustavičně se zamýšlet nad fyzikálním významem jednotlivých matematických obratů. Proto se v tomto a v následujícím odstavci zaměříme zejména na objasnění a zdůvodnění podmínek, které formalismus komplexní notace umožňují. Základem komplexní notace je Eulerova formule (viz např. [1], str. 15, 169) exp(±iP(f, t) = a exp ±i(wŕ =p k ■ f+ a) . (3) Podobně kulové harmonické vlny 1.4(11) píšeme ve tvaru \P(r, í) = -p- exp [±i(wr =p kr + a)] . (4) kr Při této úmluvě je z Eulerovy formule (1) zřejmé, že ve výrazech (3) a (4) nezáleží na znaménku u imaginární jednotky. Např. o tom, zda kulová vlna je divergentní nebo konvergentní, rozhoduje pouze to, zda znaménka u Lut a kr v argumentu fázoru výrazu (4) jsou opačná nebo stejná, a není důležité, jaké je znaménko u imaginární jednotky. Je tedy vhodné učinit další úmluvu a důsledně používat jedné z obou možností. Rozhodneme se pro znaménko mínus. (Důvod pro tuto volbu vyplyne z textu za vztahem (13) tohoto odstavce.) Rovinné resp. kulové harmonické vlny (3) resp. (4) pak charakterizují výrazy resp. SP{f, t) = a exp i(± k ■ r — u>t — a) = A exp(±i/j • r) exp(—iwí), a A SP(r, ť) = -— exp [i(±kr — Lut — a)] = -— exp(±i/jr) exp(—ÍLut), kr kr (5) (6) kde 1.5 Komplexní notace harmonických vln 13 A = a exp(—ia) (7) je komplexní amplituda obsahující informaci o počáteční fázi a. Podle definice harmonických vln a s použitím komplexní notace musí harmonickou vlnu charakterizovat vlnová funkce ve tvaru součinu &(f, ť) = ij){ř) exp(-iwŕ), (8) v němž ip(r) je funkcí jen prostorových proměnných. Rovinné a kulové harmonické vlny (5) a (6) jsou toho příkladem. Tato faktorizace výrazu pro harmonické vlny do tvaru součinu funkce polohy a funkce času je velkou předností komplexní notace. Vyplývá z ní například, že lineární kombinace harmonických vln \Pj(r,ť) určité frekvence lu je touž lineární kombinací funkcí ipj(f) závislých jen na souřadnicích násobenou fázorem exp(—iujt): A1V1(ŕ,t) + ---+AnVn(ŕ,t) = [A^ir) + ■ ■ ■ + Ani>n(r)} exp(-iwt). (9) Při difrakci vlnění nedochází ke změnám frekvence. Můžeme tedy při výpočtu difrakčních jevů harmonických vln počítat jen s faktorem vlnové funkce, jenž závisí pouze na prostorových proměnných, a považovat za samozřejmé, že časovou závislost charakterizuje fázor exp(—iujt) a vůbec jej nepsat. Podle těchto úmluv charakterizuje tedy výraz tjj(ř) = Aexp(ik ■ ř) (10) rovinnou harmonickou vlnu postupující ve směru vlnového vektoru k, kdežto výraz tp(r) = Aexp(—ik ■ r) (11) představuje rovinnou harmonickou vlnu postupující proti směru vektoru k. Podobně výraz exp^r) kr charakterizuje divergentní kulovou harmonickou vlnu, zatímco výraz exp(-ifcr) ip(r) = A--- (13) kr představuje konvergentní kulovou harmonickou vlnu. Je zřejmé, že uvedená interpretace výrazů (10) až (13) je důsledkem toho, že jsme zvolili záporné znaménko u imaginární jednotky ve fázorech výrazů (3) a (4). Udělali jsme to proto, že při této volbě je kladné znaménko v argumentu rovinné harmonické vlny (10) postupující ve směru k a v argumentu fázoru divergentní kulové harmonické vlny (12). Většina současných autorů volí rovněž záporné znaménko u imaginární jednotky v (3) a (4). Bohužel, ne všichni. Zejména v prvních desetiletích dvacátého století se autoři často rozhodovali pro znaménko plus a i dnes se s touto volbou nezřídka setkáváme (viz např. [2], str. 8, [3], str. 47, [4], str. 51). Představme si nyní, že máme vypočítat nějakou fyzikální veličinu, jež je lineární kombinací (konečného či nekonečného počtu) vlnových funkcí !řj. Uvažme, že lineární kombinace reálných částí funkcí \Pj je reálnou částí lineární kombinace funkcí SPf. ĎjRe (!řj) = Re^^ bj\Pj, b j = reálné koeficienty. (14) i i Použijeme tedy komplexní vyjádření vlnových funkcí \Pj a vytvoříme jejich lineární kombinaci. Ta představuje komplexní vyjádření počítané veličiny. Z ní pak vezmeme jen její reálnou část, neboť jen ta odpovídá fyzikální realitě. Tohoto postupu budeme v dalším mnohokrát používat při matematickém popisu interferencí (součet konečného počtu vlnových funkcí) a difrakčních jevů (součet nekonečně mnoha vlnových funkcí, tj. difrakční integrál). Jinak je tomu, počítáme-li fyzikální veličinu, jež souvisí s jednou nebo několika vlnovými funkcemi jinak než lineárně, např. druhou mocninou nebo součinem. V těchto případech ovšem neplatí, že např. součin reálných částí vlnových funkcí je reálnou částí součinu vlnových funkcí. Proto nelze jednoduše používat komplexní vlnové funkce, a musíme počítat s reálnými vlnovými funkcemi (byť třeba vyjádřenými pomocí komplexních funkcí, jako např. ve výrazech (2)). Výpočty pak bývají mnohem komplikovanější. 14 1 SKALÁRNÍ VLNA A JEJÍ MATEMATICKÝ POPIS Naštěstí jedinou veličinou tohoto typu, s níž se budeme setkávat, je intenzita vlnění. Pojednáme o ní hned v následujícím odstavci a podáme tak příklad, jak se počítají veličiny, jež s vlnovou funkcí souvisejí nelineárně. 1.6 Intenzita vlnění a výpočet intenzity harmonických vln Při pozorování a registraci difrakčních jevů je daleko nejdůležitější veličinou tzv. intenzita světla nebo obecně intenzita vlnění /. Její vyjádření prostřednictvím komplexní vlnové funkce je - aspoň v případě harmonických vln - velmi jednoduché. Intenzita je totiž úměrná čtverci modulu prostorové části ip(f) vlnové funkce. Konstanta této úměrnosti není pro difrakci významná, a proto klademe I{r) = \ip(ř)\2 = i){ř)i)*{ř). (1) Tuto skutečnost nyní zdůvodníme. Především uvedeme, co je to intenzita vlnění. Jak je známo z fyziky (viz např. [5], §34.5, [6], §7.2), je intenzita vlnění skalární veličina představující časovou střední hodnotu hustoty výkonu vlnění. Jinými slovy, je to časová střední hodnota velikosti hustoty toku energie vlnění. U elektromagnetického vlnění je to časová střední hodnota velikosti Poyntingova vektoru. Z praktického hlediska je podstatné, že intenzita vlnění je v naprosto převažujícím počtu případů tou veličinou, kterou detekujeme. Konkrétně, měříme-li výkon vlnění nějakým detektorem s malou vstupní aperturou a je-li plocha detektoru rovinná a kolmá na směr šíření vlnění, je intenzita vlnění v místě detektoru časovou střední hodnotou výkonu dělenou velikostí vstupní plochy. Při pozorování difrakce světla je intenzita světla tou veličinou, která je úměrná jasu matnice. Při fotografování difrakčních jevů souvisí zčernání fotografické emulze s intenzitou světla prostřednictvím tzv. osvitu (expozice), což je součin intenzity světla a expoziční doby. Pro nás je důležité, že u všech druhů vlnění je intenzita vlnění úměrná kvadrátu (reálné) vlnové funkce. U harmonických vln má podle 1.5(8) tento kvadrát tvar jíle^^í)] j2 = |Re[i/i(r)exp(-iwí)]|2 = |Re{[Re[i/i(r)]+ilm[i/i(r)]] (coswí - isinwí)}}2 coswŕ +Im[V>(ŕ)] sinwŕj = = JRe[V>(r)]j cos2wí+|lm[i/i(r)]| sin2 tut + 2Re[i/>(r)]lm[i/>(r)] sincijícos u>t. Vyjádřením trigonometrických funkcí funkcemi dvojnásobného argumentu dostaneme [ne[v(f,t)]y {Re[i/.(r)]}2 + {lm[i/.(f)]} JRe[i/.(f)]}2 - {lm[V>(r)]} cos 2uot + Re[V>(ř)]lm[V>(ř)] sin2wí. (2) Ze vztahu (2) je vidět, že kvadrát (reálné) vlnové funkce - a tedy i hustota výkonu harmonických vln - má složku na čase nezávislou a složku periodicky se měnící s časem. U většiny typů vlnění není možné tuto periodicky se měnící složku sledovat ani okem ani žádnými přístroji. (Frekvence viditelného světla je téměř 1015 Hz, rentgenového záření dokonce 1018 Hz.) Jsme schopni detekovat pouze časovou střední hodnotu hustoty výkonu, tj. veličinu úměrnou výrazu JRefiř^í)] j2 = lim - í JRe[!ř(r,r)] j2dr. Snadno nahlédneme, že [Re[V(f, t)] j' = \ {Re[^(r)] }' + {lm[^(r)] }' Střední hodnota zbývajících dvou sčítanců ve (2) je totiž nulová, neboť (3) (4) 1.7 Helmholtzova rovnice 15 1 [T i. sin2ciJT hm — / coszcijídí = hm - = 0, t 1 /"T • r> i i- 1 —cos2cijt hm — / sinzwŕdŕ = hm--- = 0. t^oo T Jq t^oo 2(jJT Je tedy intenzita harmonického vlnění úměrná čtverci modulu prostorové části ip(f) komplexní vlnové funkce. Konstanta úměrnosti není při popisu difrakčních jevů důležitá, a proto se pod pojmem intenzita I vlnění rozumí pouze čtverec modulu prostorové části vlnové funkce, jak to odpovídá vztahu (1). Možná toto vypuštění konstanty úměrnosti spolu s častým používáním slova intenzita v různých oborech vedlo některé autory učebnic k tomu, aby pro veličinu I určenou vztahem (1) volili jiné názvy. Používá se např. jas (brightness), ozáření (irradiance), osvětlení (illumination) atd. My však budeme veličinu I důsledně nazývat intenzitou. Z výrazů 1.5(10) a 1.5(11) je vidět, že intenzita rovinné harmonické vlny je nezávislá na poloze bodu pozorování: I = AA* = konst. Naproti tomu z 1.5(12) a 1.5(13) vyplývá, že intenzita kulové harmonické vlny klesá se čtvercem vzdálenosti od zdroje resp. od bodu konvergence: I = AA*/{kr)2 . Tak tomu také musí být podle zákona zachování energie. 1.7 Helmholtzova rovnice Budeme nyní hledat řešení vlnové rovnice 1.1(1) metodou separace proměnných. Konkrétně, vyjádříme řešení \P(r,ť) ve tvaru superpozice partikulárních řešení \Pm(r,ť) a o těchto partikulárních řešeních budeme předpokládat, že jsou součinem dvou funkcí, z nichž jedna závisí pouze na prostorových souřadnicích r a druhá pouze na čase t. Symbolicky zapíšeme tuto superpozici součtem ^(r,í)=]>>m(r,í), (1) m kde *m(r,í)=Vm(íOrm(í). (2) Výraz (1) je vskutku jen symbolický zápis superpozice partikulárních řešení \Pm(r,ť). V tuto chvíli nespecifikujeme, zda jde o součet nebo o integrál, ani meze sčítání či integrace. To vše se vyjasní až v průběhu hledání partikulárních řešení ve faktorizovaném tvaru (2). Je hodno pozoru, že požadavek faktorizace partikulárních řešení implikuje rozklad vlnové funkce \P(r,ť) do harmonických vln. Dosadíme-li totiž do vlnové rovnice 1.1(1) partikulární řešení ve tvaru (2), dostaneme V2i/>m(r)Tm(í) = ^ipm(r) d T^!r ■ Podělíme-li tuto rovnici partikulárním řešením (2), dostaneme rovnici V^m(r)_ i %^ 1pm(f) V2 tm(t) jejíž levá strana závisí pouze na souřadnicích, kdežto pravá strana pouze na čase. Protože tato rovnice musí platit pro všechny polohy r a v každém okamžiku t, je zřejmé, že ji lze splnit jen tak, že její pravá i levá strana jsou rovny téže konstantě. Označme tuto konstantu n. Pravá strana rovnice (3) tedy vede na obyčejnou diferenciální rovnici d tm(t) dt2 = V2ktt Át), (4) z níž je zřejmé, že separační konstanta n má fyzikální rozměr délky na méně druhou. Řešeními rovnice (4) jsou funkce expiv^/Jít) a exp(—v^fňt). Protože máme co činit s vlněním v neabsorbujícím prostředí, je smysluplné předpokládat, že tato řešení jsou periodickými funkcemi času. Tomu odpovídá reálná a záporná hodnota separační konstanty k. Položíme tedy k = —k2, tj. yTt = ±ik, kde k je veličina s rozměrem délky na mínus prvou a s nezápornou hodnotou. Každé takové hodnotě k přísluší výše uvedená 16 1 SKALÁRNÍ VLNA A JEJÍ MATEMATICKÝ POPIS dvojice řešení. Je to tedy právě hodnota k, která partikulární řešení (2) specifikuje, a proto ji použijeme místo indexu m. Časově závislý faktor partikulárního řešení (2) má tedy tvar T~k(t) = Ci exp(ivkt) + C*2 exp(-ivkt), (5) kde Ci a Ci jsou konstanty. Je zřejmé, že nezáporná konstanta k souvisí s kruhovou frekvencí lu vztahem vk = uj, tj. k = —. (6) v Také levá strana rovnice (3) je rovna konstantě —k2, čímž dostáváme tzv. Helmholtzovu rovnici nazývanou též stacionární vlnová rovnice: V2^fe(rO + k2Mr) = O- (7) V bodech homogenního, izotropního a neabsorbujícího prostředí neobsahujících zdroje vlnění musí tedy faktor ipk(f) partikulárního řešení fc(r) Helmholtzovy rovnice (7) lze vždy vyjádřit ve tvaru Mrf)=^kr)=^(p)1 (8) kde ip(p) je řešení tzv. normalizované Helmholtzovy rovnice V}iKp)+iKp) = 0. (9) Označíme-li totiž v (9) p = kr, je takže Vp-^(p) + 1>(p) = -^V2^(kf) + ij}{kr) = [V}^{kr) + k?i/>(křj\ = 0. Z toho je vidět, že funkce ipk{f) i ip(kr) jsou řešeními téže Helmholtzovy rovnice (stejné k2), takže platí (8). Při specifikovaných okrajových podmínkách má Helmholtzova rovnice (7) dvě řešení ij}{kr) a ip(—kr). Každé hodnotě k > 0 přísluší tedy čtyři partikulární řešení (2) vlnové rovnice 1.1(1), a sice ip(kr) exp(ikvt), ij}{kr) exp(—ikvt), ip(-k'r) exp(ikvt) a ip(—kr)exp(—ikvť). Řešení (1) vlnové rovnice 1.1(1) má tedy tvar superpozice harmonických vln vyjádřené integrály: \P(r, ŕ) = / tjj(kr) exp(ikvt) dk + / tjj(kr) exp(—ikvt) dk + Jo poo tp(—kr) exp(ikvť) dk + / tp(—kr) exp(—ikvť) dk. (10) o V dalším se budeme vždy zabývat monochromatickými vlnami. Nebudeme tedy potřebovat integrály (10) a pro řešení Helmholtzovy rovnice (7) budeme používat symbolu ip(f) (tj. budeme psát ip(f) místo ip(kř) = ipk(r)). 1.8 Poznámka o polárním tvaru řešení Helmholtzovy rovnice 17 1.8 Poznámka o polárním tvaru řešení Helmholtzovy rovnice V monografiích se často vyskytují formulace vzbuzující dojem, že reálná část výrazu V{r,t) = a(ŕ)expj-i[W - ip(f)] j, (1) v němž a(r) > 0 a (p(ý) značí reálné funkce, vždy charakterizuje harmonickou vlnu (viz např. [6], vztah (23) v odst. 1.3.3, [4], vztah 2.2-1). Funkce ip(r) = a(f) exp[iy(ř)] (2) však nemůže - bez dalších podmínek - být prostorovou částí harmonické vlny, neboť jde o obecnou komplexní funkci reálných proměnných zapsanou v polárním tvaru. Má-li být prostorovou částí harmonické vlny, musí funkce (2) vyhovovat Helmholtzově rovnici. Dosadíme ji tedy do Helmholtzovy rovnice a dostaneme podmínky, které musí splňovat reálné funkce a(r) a ip(r), aby funkce (2) představovala prostorovou část harmonické vlny. Derivováním vypočteme: Vi/>(r) = exp[iy(ŕ)] Va(ŕ) + ŕ0(ŕ) Vy(r), V2i/>(r) = exp [iip(ř)] V2a(r) - tp(r)[Vip(ř)] 2 + iJ2 exp [iip(r)]Va(r)Vip(r) + tP(r) V V(r)j. (3) Dosazením výrazu (3) do Helmholtzovy rovnice dělené funkcí (2), tj. do rovnice V2i/>(r)/i/>(r) + k2 = 0 a oddělením reálné a imaginární části dostaneme hledané podmínky: V2fl(r1-[V^(r)]2 + fc2 = 0, (4) a (ŕ) Vip{r) + Vžip{r) = 0. (5) a(f) Při odvozování základního vztahu geometrické optiky, rovnice o eikonalu (viz např. [4], str. 58), je důležitý vztah vyplývající ze (4) |V^(r)| = Jk2 + ^^. (6) Přesvědčíme se, že vztahy (5) a (6) jsou splněny pro rovinnou a kulovou harmonickou vlnu. V případě rovinné harmonické vlny je (viz 1.5(10) a (11)) a(r) = A = konst. a (p(ý) = ±k ■ r. Je tedy Vo(r) = 0, V2a(r) = 0, Vy(f) = ±k, W2if(r) = 0, takže vztahy (5) a (6) jsou splněny. V případě kulové harmonické vlny je (viz 1.5(12) a (13)) a(f) = A/{kr), f(f) = ±kr. Je tedy Vfl(r) = ^fAr/(kr3), V2a(r) = 0, Wif(r) = ±kr/r, W2if(r) = ±2k/r, takže vztahy (5) a (6) jsou splněny. Rovinné a kulové harmonické vlny bývají uváděny jako příklady tzv. homogenních harmonických vln, u nichž jsou místa konstantní fáze totožná s místy konstantní amplitudy (viz např. [6], §§ 1.3.3, 11.4.2). U rovinné harmonické vlny je fáze konstantní v rovině kolmé na vlnový vektor a amplituda je konstantní v celém prostoru. U kulové harmonické vlny je amplituda i fáze konstantní na každé kulové ploše se středem ve zdroji kulové vlny. Výraz (1) resp. (2) za splnění podmínek (5) a (6) však může představovat obecnější typ vlny, kdy plochy konstantní amplitudy a(f) = konst. nemusí být totožné s plochami konstantní fáze tp(r) = konst. Takové vlny se nazývají nehomogenní harmonické vlny. Ve smyslu definice z odst. 1.1 vlnoplocha v případě nehomogenních harmonických vln neexistuje. V literatuře se však v případě harmonických vln pod pojmem vlnoplocha rozumí plocha konstantní fáze (viz např. [4], str. 51) a v tomto významu budeme nadále pojmu vlnoplocha používat. Příkladem nehomogenních vln jsou tzv. evanescentní vlny. Jsou to např. vlny, které při totálním odrazu „prosakují" do opticky řidšího prostředí. Ve skalární teorii difrakce se s evanescentními vlnami setkáme při rozkladu vlnové funkce do rovinných vln (viz závěr odst. 7.3). O evanescentních vlnách a jejich aplikacích pojednává např. stať [7] a monografie [8]. Zmínka o evanescentních vlnách by mohla vzbudit dojem, že nehomogenní vlny jsou něčím mimořádným a ojedinělým a že většinu vlnových jevů popisuje homogenní vlna. Není tomu tak, spíše opak 18 1 SKALÁRNÍ VLNA A JEJÍ MATEMATICKÝ POPIS je pravdou. Vezměme si například superpozici dvou rovinných harmonických vln o stejné amplitudě a vlnové délce šířících se ve směrech určených jednotkovými vektory ň\ a ň2: ip(f) = a exp(ikňi ■ r) + a exp(ikň2 ■ r). Vyjádřeme tuto komplexní vlnovou funkci tp v polárním tvaru. Vytknutím faktoru aexp (ifc"1"^""2 • r) dostaneme tjj(r) = aexp ik ni + n2 exp ik ni — n2 exp —ik Tli — Tl2 tp(r) = 2a cos í k Tli - Tl2 exp ik ňi + ň2 tj- Z tohoto polárního tvaru vlnové funkce je vidět, že plochy konstantní amplitudy jsou roviny kolmé k vektoru ňi — ň2, zatímco plochy konstantní fáze jsou roviny kolmé k vektoru ňi+ň2. Obě soustavy rovin nejen že nejsou totožné, ale jsou dokonce na sebe kolmé. Tento příklad superpozice dvou rovinných vln naznačuje jednak že interferenční a difrakční jevy představují vždy nehomogenní vlnu, jednak že dělení harmonických vln na homogenní a nehomogenní vlny může být jen „školská moudrost" postrádající praktického významu. 1.9 Paraxiální vlna a paraxiální Helmholtzova rovnice Uvažujme o harmonické vlně, která se šíří v nějakém směru, např. o vlně, která odpovídá světelnému svazku šířícímu se ve směru souřadnicové osy x3. Ve velké vzdálenosti od zdrojů (tj. ve vzdálenosti mnoha vlnových délek) se může taková vlna podobat rovinné vlně. (Výjimkou jsou případy typu průchod vlnění ohniskem.) Proto se často taková vlna charakterizuje výrazem ip(r) = ip0(r) exp(ikx3), (1) v němž „silnou" závislost na x3 vyjadřuje rovinná vlna exp(ikx3) a komplexní funkce ipo(f) tuto rovinnou vlnu pouze moduluje. To znamená, že se předpokládá, že ipo(f) je pomalu se měnící funkcí souřadnice x3. Konkrétně se předpokládá, že v rozmezí vlnové délky je modul změny funkce ipo(f) malý ve srovnání s IV'oC^OI a rovněž modul změny derivace dipo/dx3 je v rozmezí vlnové délky malý ve srovnání s \dipo/dx3\: \ipo(xi,x2,x3 + A) -ipo(xi,x2,x3)\ dip0(ř) dx3 « \M?)\: d2tp0(ř) dip0(ř) dx3 Častěji se tyto podmínky zapisují ve tvaru dip0(ř) dx3 (r) = exp(i/j2ľ3) dx3 (3) kde jsme použili toho, že Vi/>o ■ *3 = dipo/dx3. Dosadíme nyní výraz (3) do Helmholtzovy rovnice: V2i/>(r) + k2tp(ř) = exp(ikx3) V2^0(y) + 2i;T""v/ =0. L vx3 Má-li tedy být funkce (1) řešením Helmholtzovy rovnice, musí komplexní funkce ipo(f) splňovat rovnici dx\ 02 dx2 2ik- dx3 0. (4) Funkce ipo(f) zvolené intuitivně nebo na základě nějaké aproximace však rovnici (4) nesplňují. Přísně vzato paraxiální vlna (1) pak vůbec žádnou vlnou není. Aby bylo možné považovat paraxiální vlnu (1) za smysluplnou aproximaci nějaké vlny, požaduje se většinou, aby funkce ipo(f) vyhovovala tzv. paraxiální Helmholtzově rovnici dx\ dx\ dx3 Vyjasníme nyní, v jakém smyslu je paraxiální Helmholtzova rovnice (5) aproximací rovnice (4). Je zřejmé, že rovnice (5) se získá ze (4) zanedbáním sčítance d2ipo/dx2. Jak se však dá takové zanedbání zdůvodnit? Nic přece není malé ve srovnání s nulou, tedy ani \d2ipo/dx2\. Zanedbání sčítance d2ipo/dx2 zdůvodňujeme takto: Levou stranu rovnice (4) tvoří tři sčítanci d2ip0(r) d2ip0(ř) (0 ^2 + —572—. (6) dx\ d2ip0(r) b2 (ii) -i^, (7) b3 (iii) 2ikd-^l. (8) ox3 Z druhé nerovnosti (2) vyplývá, že modul sčítance (7) je velmi malý ve srovnání s modulem sčítance (8). Má-li však platit rovnice (4), musí být modul sčítance (7) také velmi malý ve srovnání s modulem sčítance (6). (Obě tyto skutečnosti je třeba při intuitivní volbě funkce ipo(f) ověřit.) Součet sčítanců (6) a (8) tvoří levou stranu paraxiální Helmholtzovy rovnice (5). Na její pravé straně by měl být výraz —d2ipo{ř)/dx3, jehož modul je velmi malý ve srovnání s moduly obou sčítanců na levé straně a v tomto smyslu je paraxiální Helmholtzova rovnice (5) smysluplnou aproximací rovnice (4). Má-li tedy výraz (1) představovat paraxiální vlnu, musí funkce ipo(f) splňovat nejen podmínky (2), ale také podmínku dx\ dx\ dx3 a tím - s přihlédnutím k podmínce (2) - také d2Mr) . d2Mr) 0., d^0(f) -— + ^2 ~ ~2lk (9) d2tp0(ř) dx2 d2ip0(f) d2ip0(f) dx\ dx2 (10) V následujícím odstavci tyto skutečnosti ozřejmíme na konkrétním příkladě. 20 1 SKALÁRNÍ VLNA A JEJÍ MATEMATICKÝ POPIS 1.10 Fresnelova aproximace kulové harmonické vlny Uvažujme o kulové vlně vycházející z počátku soustavy souřadnic (obr. 2), tj. o vlně m = e*p(ifcr) (1) Obrázek 2: K odvození Fresnelovy aproximace kulové vlny. Chceme-li vyjádřit kulovou vlnu v bodech P{x\, xi, x3) blízkých ose x3, tj. v bodech, jejichž souřadnice splňují podmínku ^3 e < 1, můžeme délku r aproximovat prvními několika členy mocninného rozvoje 13 1+7-t' ^3 2\2 8x3 (2) 2 8 ) 2x3 Abychom si udělali představu, na kolik členů mocninného rozvoje se můžeme omezit, zvolme hodnoty \Jx\ + x\ = 1 • 10~2 m, 13 = 1 m, k = 1 • 107 m-1, tj. A = 6, 3 • 10~7 m. Ve jmenovateli kulové vlny (1) můžeme s relativní přesností e2/2 položit r « x3. (3) Při zvolených hodnotách je relativní přesnost této aproximace 5 • 10~5. V čitateli kulové vlny je třeba použít přesnější aproximaci. Fázor exp(ikr) je totiž periodickou funkcí, a proto můžeme zanedbat jen ty členy rozvoje (2), jejichž součet A splňuje podmínku kA < 2tt, tj. A < A. Druhý člen rozvoje (2) má při zvolených hodnotách velikost (4) 2x3 = 5•IQ"5 m, a proto jej nelze zanedbat. Třetí člen rozvoje (2) má však hodnotu v 2 {xi + 4y 8x1 = 1, 25 • 10~9 m < 6, 3 • ÍO"7 m = A, takže jej už zanedbat lze. Budeme tedy fázor ve výrazu (1) aproximovat takto: exp(i/jr) « exp ik [ x3 2x3 = exp(ikx3) exp ik 2x~?, (xj + . (5) 1.10 Fresnelova aproximace kulové harmonické vlny 21 Kulovou vlnu (1) pak v bodech blízkých ose x3 můžeme aproximovat výrazem exp(ikr) exp(ikx3) X3 ■ exp (6) Ekvivalentní formy této aproximace použil Fresnel v r. 1818 při interpretaci Fresnelových difrakčních jevů, a proto se jí říká Fresnelova aproximace. V teorii difrakce se jí používá velmi často. Výraz na pravé straně je paraxiální vlnou, v níž funkce ipo(f) modulující rovinnou vlnu exp(ikx3) má tvar ip0(r) = — exp 2^3 2x- (x\ + x\) (7) Přesvědčme se, že tato funkce splňuje požadavky potřebné k tomu, aby pravá strana rovnice (6) mohla být považována za paraxiální vlnu. Za tím účelem počítejme dxi d2ipo(r) dx\ ik—ip0(r), X3 x2 1 = -k2ý0(r)-^ +ikipo(r) — %3 a podobně pro dipo/dx2 a d2ipo/dx2,. Je tedy d2ip0(f) d2ip0(f) dx\ 02 = -k2^ — + 2i^0(r)— = -k2Mr) — x3 í 1 - 2i 2^3 ki^X-^ ~\~ x2^ (8) Dále počítejme derivaci podle x3: dx3 d2ip0(r) dx2 -ikMr^+4 2x2 1 - 2i 2^3 2 | 2 \ ^ X -i | X n ki^x-^ ~\~ x2) 8 1 kx3 x\ + x\ (9) (10) Nerovnosti 1.9(2) přepíšeme do tvaru dx3 k\tp0(ř) «1, č>2V>o(r) dx% dipojr) dx3 < 1 a dosadíme do nich výrazy (9) a (10) a výše uvedené numerické hodnoty. Dostaneme &0o(f) dx 3 k\ip0(r) d2tf>o(ř) dxl dipojr) dx3 2x2 '1+4 k2(x\ + xl)2 = 5 • lO-Vl + 4-10"6 < 1, 64^ 2x2 Vl + 4- 10-6 (H) (12) Vztahy (11) a (12) ukazují, že podmínky 1.9(2) kladené na funkci ipo(r) jsou v případě Fresnelovy aproximace kulové vlny a při zvolených numerických hodnotách splněny. Z výrazů (8) a (9) je zřejmé, že funkce ipo(f) ve tvaru (7) je řešením paraxiální Helmholtzovy rovnice 1.9(5), a tím splňuje podmínky 1.9(9) a 1.9(10). Helmholtzovu rovnici 1.9(4) funkce (7) ovšem nesplňuje. Fresnelova aproximace (6) kulové vlny je tedy dobrou paraxiální aproximací, přísně vzato však nereprezentuje vlnu, neboť pravá strana rovnice (6) nevyhovuje Helmholtzově rovnici. 22 1 skalární vlna a její matematicky popis 1.11 Důležitá věta o řešení Helmholtzovy rovnice Pro teorii difrakce, konkrétně pro formulaci okrajových podmínek, má velký význam tato věta: Buď ip (ŕ) řešením Helmholtzovy rovnice V2i/>(r) + k2iP(r) = 0 (1) v £3. Nechť dále platí, že ij){ř) = 0, Vi/>(r)-n = 0 (2) v bodech libovolně malé, avšak konečné části S" nějaké plochy v £3. Pak je ip(f) = 0 v celém prostoru £3. Důkaz vychází z Greenovy formule (viz Dodatek A) (iŕV2!ŕi-iíiV2ií) dv= -V>iVV>) -ňdS, v J J s (3) Obrázek 3: Objem v v Greenově větě (6). v níž v je konečný objem uzavřený po částech hladkou plochou S, ň je vnější normála k ploše S a ip(f) a ipi(f) jsou libovolné funkce spojité se všemi svými prvními derivacemi uvnitř a na ploše S a se všemi druhými derivacemi uvnitř plochy S. Do rovnice (3) dosadíme za ip řešení Helmholtzovy rovnice (1) a za lpi funkci jež je v £3 řešením Laplaceovy rovnice ýi(ř) = 4 - -> R r V2i/>i(f) = 0. (4) (5) Přitom R značí poloměr kulové plochy k(0,R) specifikované tak (viz obr. 3), že část S" této kulové plochy spolu s částí plochy S', na níž jsou splněny podmínky (2), vymezuje objem V, v němž a na jehož okraji S = S' U S" řešení Helmholtzovy rovnice ip(r) nemění znaménko. Toho lze vždy dosáhnout tím, že zvolíme objem V dostatečně malý. Vzhledem k (1), (4) a (5) tím nabude Greenova formule (3) tvaru 1-1 U(ŕ)dV = R r J j.,s V>(r)V 1 V-0 ídS. (6) Plošný integrál je roven nule na ploše S", na níž jsou splněny podmínky (2). Na části S" kulové plochy HO, R) je REFERENCE 23 r=R r=R dr Vr 1 r r=R r=R n Ř2' (7) Je tedy plošný integrál na pravé straně rovnice (6) roven (l/i?2) JJS„ ip(r) dS". Greenova formule (6) se tím redukuje do tvaru V objemu V je R > r, tj. lR--^*V=Jě tp(f) dS" (8) (9) Předpokládali jsme však, že řešení ip(f) Helmholtzovy rovnice (1) nemění uvnitř a na okraji objemu V znaménko. Rovnici (8) lze tedy splnit jen tak, že ip(f) = 0 uvnitř a na okraji objemu V. Objem V, pro nějž je věta dokázána, můžeme nyní dále zvětšovat tak, že připojíme další objem V' ohraničený plochou ležící uvnitř dříve uvažovaného objemu V a částí nové koule k'(0', R'). Tak posléze dospějeme k tomu, že ip(f) = 0 v celém prostoru E%. Poznámky 1. V důkazu se předpokládá, že řešení ip(f) Helmholtzovy rovnice nemění uvnitř a na okraji objemu V znaménko. To znamená, že se předpokládá, že funkce ip(r) je reálná. Často však pracujeme s řešeními Helmholtzovy rovnice, která jsou komplexními funkcemi reálné proměnné. Věta ovšem platí také v těchto případech, neboť Helmholtzova rovnice je lineární a komplexní řešení lze považovat za lineární kombinaci ip(r) = Reip(r) + ilmip(f) dvou reálných řešení. 2. Analogická věta platí v prostoru E^ libovolné dimenze N > 2. V důkazu se volí řešení Laplaceovy rovnice ve tvaru ipi(r) = (1/RN~2) — (l/rN~2), pro N > 3 a ipi(r) = ln(r/R), pro N = 2. Uvedený důkaz podal pro E2 H. Weber [9] v r. 1869 a reprodukoval F. Pockels v monografii [10] z r. 1891. Od těch dob se v literatuře i přes nespornou důležitost věty důkaz neuvádí. 3. Zajímavým, i když možná evidentně ekvivalentním tvrzením s dokázanou větou, je závěr, že je-li řešení Helmholtzovy rovnice rovno nule uvnitř libovolně malého, avšak konečného objemu V, je identicky rovno nule v celém prostoru. 4. Z linearity Helmholtzovy rovnice vyplývají tyto důsledky dokázané věty: (i) Jsou-li si dvě řešení Helmholtzovy rovnice i jejich derivace ve směru normály rovna v bodech nějaké konečné plošky, jsou si rovna v celém prostoru. (ii) Jsou-li si dvě řešení Helmholtzovy rovnice rovna v bodech nějakého konečného objemu, jsou si rovna v celém prostoru. 5. Tzv. Kirchhoífovy okrajové podmínky (viz odst. 7.4.1) předpokládají, že prochází-li vlnění otvorem v nepropustném stínítku, je v bodech stínítka na odvrácené straně od zdroje vlnová funkce ip = 0 a také derivace ve směru normály Vi/> • ň = 0. V bodech otvoru ve stínítku se naopak předpokládá, že vlnová funkce i derivace ve směru normály je nenulová a táž jako v případě volného šíření, tj. jako za nepřítomnosti jakéhokoli stínítka. Takové okrajové podmínky jsou často výborným fyzikálním modelem, jsou však zřejmě matematicky rozporné. V praxi se jich nicméně používá a budeme jim v dalších kapitolách věnovat mnoho pozornosti. Reference [1] Kvasnica J.: Matematický aparát fyziky. Academia 1989. [2] Cowley J. M.: Diffraction Physics. North-Holland, Amsterdam 1975. [3] Laue M. v.: Materiewellen und ihre Interferenzen. 2. Auflage. Akademische Verlagsgesellschaft, Geest & Portig, Leipzig 1948. 24 REFERENCE [4] Saleh B. E. A., Teich M. C: Základy fotoniky 1. Matfyzpress, Praha 1994. [5] Halliday D., Resnick R., Walker J.: Fyzika. Vysoké učení technické v Brně - Nakladatelství VUTIUM a PROMETHEUS, Praha 2000. [6] Born M., Wolf E.: Principles of Optics. 7th ed. Cambridge University Press 1999. [7] Bryngdahl O.: Evanescent Waves in Optical Imaging. In Progress in Optics XI. (E. Wolf, ed.), North-Holland Publ. Co., Amsterdam 1973, 167-221. [8] Fornel F. de: Evanescent Waves. From Newtonian Optics to Atomic Optics. Springer, Berlin 2001. [9] Weber H.: Ueber die Integration der partiellen Differentialgleichung ||4r + j^ß- + k2u = 0. Mathematische Annalen 1 (1869), 1-36. [10] Pockels F.: Uber die partielle Differentialgleichung Au + k2u = 0 und deren Auftreten in der mathematischen Physik. B. G. Teubner, Leipzig 1891, 212.