Elektrodynamika 1 Elektrické a magnetické veličiny, jednotky SI Elektrický proud I je v systému SI základní veličina, jednotka je 1 Ampere (1 A). Definice: Stejné proudy ve 2 rovnoběžných drátech ve vzdálenosti 1 m mají velikost 1 A, když vzájemná přitahovací síla na 1 m drátu je 2 • 10~7N. Náboj q: Zdroj elektromagnetických sil, pohybuje se ve vodiči, když teče elektrický proud. Jednotka: 1 Coulomb (1 C). Definice: Při elektrickém proudu o 1 A proteče průřezem vodiče 1 C za sekundu. 1 C = 1 As až 6 • 1018 elementárních nábojů. Coulombův zákon Přitažlivá nebo odpudivá síla dvou nábojů q± a q2 v místech x\ a x2 F = k (x[ - x2), F = kqf, n2 n2 (1) kde ri2 = \x± — ô^l a k je konstanta. Dimenze [k] této konstanty vyplývá z (1). AT r As-As Nm2 N = M —-T- ^ [k] nr A2s2' Definice Amperu byla tak zvolena, aby k = I0~7c2 N Ä21 kde c je rychlosti světla. V SI se píše z historických důvodů 47re0' s tzv. „dielektrickou konstantou vakua" e0. Poznámka: Mohli bychom předpokládat k = 1. Tím bychom určovali dimenzi náboje, vyjádřenou mechanickými veličinami. V systému cgs dostáváme z Coulombova zákona gem \qf i s _1 dyn = —— = —-, =3- [q\ = g2 ema s . s2 cm2 To je jednotka náboje v Gaufiově systému. Elektrické pole: Síla souboru nábojů qi,..., q^ v místech další náboj q v bodě x je popsaná působením elektrického pole v bodě x na náboj. Intenzita eletrického pole E(x) je lokální vlastnost prostoru, v němž se nacházejí elektické náboje. i N ~" — ~" Intenzita = síla na jednotkový náboj. Dimenze: [E] ~ č - ^č- Napětí U. Práce = změna potenciální energie při pohybu náboje v elektrickém poli je daná integrálem W = í Fds. V případě pohybu o vzdálenost / podél homogenního elektrického pole (např. v deskovém kondensátoru) platí jednoduše W = qEl =: qU. Elektrické napětí je rozdíl potenciální energie jednotkového náboje na dvou bodech. Dimenze: [U] = ^ Jednotka: 1 Volt = IV = l£. Volt se používá pro běžné označení jednotky elektrického pole, [E] = 1^. Magnetické pole. (Stacionární) elektrické proudy vyvolávají síly na pohybující se náboje (Lorentzova síla). Analogicky k zavedení intenzity elektrického pole uvažujeme sílu souboru eletrických proudů v daném uspořádání vodičů na úsek dl jednoho dalšího drátu s konstantním proudem I. Síla je úměrná I a dl, dF(x) = Idíx B(x). (3) Magnetická indukce B(x) zahrnuje působení všech elektických proudů v bodě x. Dimenze: Z (3) vyplývá * ™ ^ N J Js Vs N = Am B , =>• B = — = — - = — - = -. Am Anr Cm^ nr Jednotka: 1 Tesla = 1T = 1% Příspěvek elektrického proudu v úseku drátu dl na místě X2 k magnetickému poli v bodě x\ (analogicky Coulombovu zákonu) je dB(xx) = kmIdíx , f1 ~ Z2|q. (4) Fl — x2\ V SI se píše km = fiQ = je tzv. „magnetická permeabilita vakua". Síla mezi dvěma infinitesimálními úseky vodiče, umístěnými v bodech s elek- trickými proudy I\ a I2 (analogon Coulombova zákona): F{xx) = ^hl2 dh x (df2 x f 1 " *2 ^ 47T \ \Xi — X2\ 2 Hustota náboje: p(x) lim av^o AV kde q je naboj v objemu AV, který obsahuje bod x. Hustota proudu: kde A A je element průřezu vodiče a je jednotkový vektor ve směru elektrického proudu J. 2 Maxwellovy rovnice, statický případ J. C. Maxwell našel v 19. století dvě vektorové a dvě skalární parciální diferenciální rovnice 1. řadu, které spojují elektrické a magnetické pole navzájem a s elektrickým nábojem. Z těch rovnic lze odvodit celá elektrodynamika. divE P_ 1 a 9E _ — rotB = e0 — + j, jjLq Ol TOtE = div B = 0. dB ~dt'' Vektorové operátory div a rot lze vyjádřit pomoci operátoru Nabla, V - (— — —\ [dx' <9í/' d z I (5) (6) div v = V • v, rot v = V x v. (7) V případě statických polí, když časové derivace jsou nulové, odpojují se elektrické a magnetické pole. Z Maxwellových rovnic plyne /em rovnice kontinuity t. j. zachování náboje. % + diví=0, Elektrostatika Základní rovnice divE P_ eo' rotE = 0. 3 Elektrostatické pole je bezvírové pole se zřídlem ^. Pro takové pole existuje skalární potenciál ^lř <27> Greenova funkce není jednoznačná. G(x,x') plus libovolné řešení homogénne diferenciální rovnice A0 = 0 je také Greenovou funkcí stejného operátoru. Greenova věta Odečtením dvou identit, V(wVw) = u Av + Vu V v a V (v Vit) = v Au + Vv Vit, dostaneme V(it Vv — vVu) = u Av — v Au. Integrovaním a aplikací Gaufiovy věty odvodíme Greenovu větu: / [u Av — v Au) dV = / (u\7v — v\7u) řídS, Jv Jav (28) 6 kde V je oblast a n je normálový vektor na okraj dV. Podle okrajových podmínek má Greenova věta dvě standardní aplikace. Dirichletův problém: Známe p(x) uvnitř nějakého objemu V a potenciál (x')ndS'. (29) Jv e0 Jav Neumannův problém: Známe gradient potenciálu na dV. Vybereme Neumannovu Greenovu funkci G^, jejíž gradient se rovná nule na okraji. (j)(x) = - í GN(x,x')^^-d3x'- í GN(x,x')f{x). Energie dvou nábojů = energie náboje q± v bodě x\ v potenciálu vyvolaném nábojem q2 v bodě x2 plus výraz s přehozenými náboji lomeno dvěma, aby se energie nepočítala dvakrát. 1 U = ~(qi M%i) + 92 01 (x2)), kde (f)i i^X^j =--~-. 47re0 \x — x*i \ Energie spojitého rozložení náboje: U = ^jp(f>dV. (31) Pro rozložení bodových nábojů je potenciál 004) = T— Y,—, rab=\xa-xb\ (32) 47re0 ^ rab a energie U=—Ti^. (33) 8^0 t?a rab 7 5 Multipólový rozklad pole V následujícím budeme hledat rozvoj elektrického potenciálu ve velké vzdálenosti od zdroje. 5.1 Laplaceova rovnice ve sférických souřadnicích Laplaceova rovnice A0 = 0 (34) platí všude, kde se nenachází náboj. Ve sférických souřadnicích má Laplacián následující tvar: 1 d ( 2ď0\ , 1 d ( , qď0\ , 1 ď20 A" ~~7 T7~ I ' TT- ) • ~7> 77 777 777 ' •>:•>, „ (35) Separace proměnných: Pokoušíme se přeměnit parciální diferenciální rovnici v tří obyčejné diferenciální rovnice pro funkce jednotlivých proměnných, R(r), ©(■$) a K tomu předpokládáme řešení ve tvaru součinu 0(r, ů, ip) = R{r) ■ &(ů) ■ $(y>) (36) a dosadíme do Laplaceovy rovnice 1 d ( 2 di? ^ \ 1 <9 / n „ d0 \ 1 „ ^ d2$ rz or \ dr ) rA sin ů oů \ dů J rz sin » d(^ 2.2..... (37) Násobíme r ^q^ů a píšeme část, závisející na íp, na pravou stranu sin2?? d / 9dR\ sinů d / d0\ 1 d2$ --r — H---sin^- =--- R dr\ dr) O dů\ dů) $ dp2' Levá strana teď závisí na r a i), pravá strana na p. Z toho vyplývá, že se obě strany musí rovnat konstantě, kterou nazveme m2. Z pravé strany dostaneme obyčejnou diferenciální rovnici, d2$ „ -— + m2$ = 0. 38 dpz Rovnici vyplývající z levé strany můžeme upravit podobným způsobem, 1 d / , dR\ m2 lid/, d6\ ~ --r — = —2^--777^ sin$ — = A , Rdr \ dr) sin ů sin wtídw \ d» / kde se opět obě strany musí rovnat konstantě, označené A2. Z toho dostaneme dvě další obyčejné diferenciální rovnice, ±(rf)-#R = 0 (39) dr V dr / d / . ndQ\ m 2 si^7777 M ^" 7777*77 © = 0, (40) sin ů dů y dů I y sin ů tak že místo parciální diferenciální rovnice máme rovnice (38), (39) a (40). Řešení: Rovnice (38) má řešení &m(p) = Cm cos mip + Sm sin mip (41) příslušné k parametru m, který musí být celočíselný, aby $ bylo periodické ve ip. Radiální rovnice (39) má řešení Ri{r) = Airl + Bi_ r-l+l : (42) kde A2 = /(/ +1). V rovnici (40) píšeme cosů = x. Řešení, které obsahuje obě separační konstanty mal, označíme Přm. Takovou úpravou dostaneme Legendreovu rovnici (1-x2) d2Přm(:r) dPřm(x dx2 2x- dx 1(1 + 1) m 1-x2 PT(x) = o. (43) 5.2 Legendreovy polynomy Ortogonální bází řešení pro m = 0 jsou Legendreovy polynomy Pi(x), vyhovující jednodušší rovnici d dx 9,dPi(x) + l(l + l)Pl(x) = 0 (44) s nezáporným celočíselným parametrem /. Legendreovy polynomy jsou ortogonální v intervalu (—1,1) J Pk(x)Pi(x)dx = 0 Wk^l. (45) Legendreovy polynomy se objevují jako koeficienty v rozvoji tzv. vytvářející funkce (l_2xí + í2)V2 ]^Pl{x)tl- Použitím Leibnizova pravidla dm[f(x)g(x)} dxm E k=0 m\ dm f(x) d g (x) k\ (m — k)\ dxm~k dxk dostaneme m-násobným derivováním rovnice (44) (1 - x2)f"(x) - 2x(m + l)f'(x) + (/ - m){l + m+ l)f(x) = 0, (46) (47) (48) kde f{x) = dmPi(x)/dxm. Substituce f{x) = (1 — x2) ml2g{x) vede k tomu, že funkce g(x) musí splňovat rovnici (43), je tedy konečně prw = (i--2r/2^1- Legendreovy polynomy lze vyjádřit pomocí Rodriguesova vzorce 1 dř Pi(x) /! 2ř dx1 9 (x2 i r. (49) (50) Využitím tothoto vztahu můžeme rozšířit (49) na oblast záporních m, tedy P™(x) = L—;L (1 - x2)m/2---(l-x2)1, -Km (ů,° = S7T / / 1 1-(3 cos2??- 1) Y£= -\ — sintfcostfeT 16n V 87r Velmi důležitým speciálním případem rozkladu (54) je vztah pro Legendreův polynom obecného úhlu 7 mezi dvěma vektory n = (sinůcos99, sin-$sin99, cosů) a n' = (sin a cos sin a sin cos a), tedy cos 7 = n ■ n' = cos ů cos a + sin $ sin a cos(ip — /?), 47T Pz(cos7) 21 - E ^m*(«,/3)^m(^^)- (56) 10 5.4 Multipólový rozklad rozložení náboje Uvažujeme rozložení náboje uvnitř koule o poloměru R, ^ = {o r>R ^ Potenciál mimo koule je dán vzorcem (25) Ane0J \x-x\ Ane0 J Jx2 + x2 - 2\x\ ■ \x'\cos-f 1 1 /"dV , P(£/) (58) ^0lXlJ yi-2fcos7+(f)2 7 je úhel mezi i a f'. V integrálu se objeví vytvářející funkce Legendreových polynomů, tak lir 00 /V\Z =---/ dVp(:r')£Pz(cos7) - (59) 47re0 r J ^ V r J (psali jsme \x\ = r a \x'\ = r'). Použitím (56) dostaneme rozvoj 1 r 4-ir r'1 m = ^ J d3x'p(x') J2 ^Y^ŽTT Yr^', V'WW, f)- (60) Pomocí multipólových momentů qlm:= í dVp(£')/r(^>') (61) můžeme konečně psat potenciál jako superpozici kulových funkcí 1 Ä ^ lq,. W = Ě0 [=0m=-I V případě bodového náboje víme, že pole je dáno Coulombovým potenciálem. Je-li náboj q umístěn mimo počátek souřadnic, např. na ose z (v bodě z = R), je potenciál dán vztahem q ~ /r (63) r 47re0r ^ 'V ; V r 7 ' Pro r y> R převažuje rotačně souměrná (vzhledem k počátku souřadnic, nikoli poloze náboje) složka / = 0. Umístíme-li však na ose z ještě náboj opačné velikosti do z = —R, vyruší se identické příspěvky členů s / = 0 a pro r >>> R převažuje pak dipólová složka (1 = 1) 2qRP1(cosů) D cosi? 0dip = "j--ž- = ~A--2~' (64) 47re0 rA 47re0 rÁ 11 kde D = 2qR označuje dipólový moment. Podobně, umístíme-li na ose z v z = ±R náboje g a v počátku náboj — 2g, vyruší se identické příspěvky členů s / = 0 a / = 1 a pro r y> R převažuje pak kvadrupólová složka (/ = 2) 2qR2P2(cos$) Q l-3cos2tf . . 0quad =----3- = "--3-, (65) Anen rá Anen rá kde Q = qR2 je kvadrupólový moment. Obecně jsou multipólové momenty závislé na umístění v souřadném systému, s výjimkou nej nižšího nenulového momentu. 6 Magnetostatika 6.1 Analogie mezi elektrostatikou a magnetostatikou Integrální tvar infinitesimální rovnice (4), vyjádřený pomocí hustoty proudu, je f ,3 x-x _ p0 f 3 , , ^ 1 B(x) = — / dVj(i') x ——^t- = -— / dVj(i') x V An J \x — x \6 Au J \x — x | An -V x /d3x' Z{Xih =: rot A{x). (66) J \x — x I Zavedli jsme vektorový potenciál l(f) = ^/,dVT|2lT. (67) v ' An J \x-x\ v ' Vektorový potenciál není jednoznačný, protože můžeme přičíst gradient libovolné funkce, jehož rotace je identicky nulová. Taková transformace, A(x) —> A(x) + grad k(x) se nazývá kalibrační transformace. (Stejně je skalární potenciál jednoznačný jenom až na konstantu.) Coulombova kalibrace div A = 0 vede k značnému zjednodušení: Dosazením do Maxwellovy rovnice dostaneme rotB(x) = iotľotA(x) = grad div A (x) — AA(x) = — AA(x) ^° ľ j3 / r?(->'\ a 1 dV/(f) A—— = p0m, (68) An J \x — x \ tedy vektorovou Poissonovu rovnici pro A. Následující tabulka ukazuje analogie mezi elektrostatikou a magnetostatikou. 12 Elektrostatika Veličina, vztah Magnetostatika dF = dg E definice pole dF = J dfx B dg hustota zřídel I dl = jd3x F = qE síla na náboj F = qv x B div E = p/e0 rot Ě = 0 rovnice pole rot B = pQj div B = 0 = —grad 0 potenciály 5 = rot A A0 = -^ rovnice potenciálu AA = -hoj (f)(x) = —— /d3x' 47T60 J potenciál zřídel Analogické magnetické pole k elektrickému Coulombovu poli je pole lineárního vodiče. Uvažujme nekonečně dlouhý, rovný drát ve směru osy z s elektrickým proudem I j(x)=IÔ(x)ô(y)e3. (69) Podle (66) je magnetické pole - Hol [°° x- (0,0, z') B[x) =- / dz e-i x —---—. 70 Zavedeme polární souřadnice p = \/x2 + y2, ip = arctan^, z a uvedomíme se, že pole musí být konstantní ve směru z, tak že stačí počítat B v rovině (x, y). Výpočet vede k Biotovu-Savartovu zákonu B{p) = f^. (71) Z7T p 6.2 Magnetické pole kruhové smyčky. Do vztahu pro vektorový potenciál (67) dosadíme hustotu proudu j(x') dV = 15{p' - a) ô(z') ey p' dp' dz' dep', kde ey = — sin(p' — dt m- (103) Analogicky dostaneme v statickém případě E = 0 B ■ tds = dS rot B ■ n dS = /iq / f ■ n dS. s Js (104) Integrál na pravé straně představuje celkový elektrický proud J, z toho vyplývá Ampěreův zákon _ as B ■ tds = (105) 8.2 Aplikace — vzájemná indukce a vlastní indukce Uvažujeme dvě geometricky pevné smyčky s proměnným proudem ve smyčce 2. Indukované napětí ve smyčce 1 vyvolané změnou pole buzeného smyčkou 2 je d B2 ■ ňi dSi dt J(i) Po dosazení dostáváme ř/i = M- (i) B2 ■ ňi dSi Ä2.dilt M = ^í ^. (106) (i) 4?r .1(2) r12 dh 12 Mi Ho 12 di2 ■ d£i dí 47T J(l) J(2) \Xi — x2i Pokud by tekl proměnný proud smyčkou 1, bylo by indukované napětí ve smyčce 2 dh (107) U2 = M21-^ dt M12 = M, 21 M. (108) Ale také změna magnetického toku smyčkou 1 vytvoří indukované napětí v této smyčce, stejně platí pro smyčku 2. Obecně tedy můžeme psát r dh d/2 , Ah T dl2 ^ = Mďľ-Í2ďľ Časová změna energie magnetického pole je rovna záporně vzaté práci dW dt -U1I1-U2I2 = L1L (109) (110) takže pro energii magnetického pole je W= ^Lilf + ^Ltlš-Mhh, LľL2>M2. (111) Energii magnetického pole máme ovšem také vyjádřenou jako W = - í B-HdV = - í j-ÄdV. (112) 2 Jv 2 Jv Při odvození obou výrazů v této rovnici je postupně využito vztahů B = rot Ä, H ■ rot Ä - Ä ■ rot H = div (Ä x H) , rot H = j. (113) Vztahu pro energii využijeme pro výpočet vlastní indukčnosti L = —p: B2dV. (114) Polz Jv Uvažujme dvě solenoidální cívky, každou o N závitech, těsně na sobě. Průřez cívek je S* a jejich délka l. Pole první a druhé cívky jsou tedy přibližně (Ampěreův zákon) Bi ~ ——, B2 « —— (115) a pro indukčnosti máme Lx^L^-M^^p^. (116) Pro energii magnetického pole pak W = ^(I1 + I2f. (117) 9 Časově proměnná elektromagnetická pole 9.1 Dynamické potenciály, kalibrace Předpokládáme dynamické potenciály o A0+-divA = dt e0 AA - ĚoVo-QJj- ~ grad \dwA + eo/^o^ I = -fM)J- (119) 19 S využitím kalibrační transformace Ä^Ä+gradý, -►-^ (12°) můžeme mít Lorenzovu kalibraci (Ludwig Valentin Lorenz ^ Hendrik Antoon Lorentz) divl+eo/xo^ = 0 (121) a dostáváme tak pro potenciály nehomogenní rovnice (122) Označili jsme rychlost světla ve vakuu Rovnice (122) jsou Maxwellovy rovnice pro potenciály, spolu s kalibrací (121) jsou ekvivalentní (5). 9.2 Rovinná a kulová vlna V případě volného elektromagnetického pole popisují homogenní rovnice odpovídající (122) šíření vln. Vlnová rovnice v jednorozměrném případě popisuje rovinnou vlnu (ve směru x) d2^(x,t) ld2^(x,t) dx2 c2 dt2 "U- [lZd) Obecné řešení je #M) = /(t-^)+$(t-^). (124) Vezmeme jako príklad Gaufiovu funkci / = exp[— (t — ^)2]. Maximum se nachází při t — | = 0, pohybuje se tedy rychlostí c ve směru rostoucího x. Dalším jednorozměrném příkladem je sféricky symetrická vlnová rovnice, ld2(rijj(r,t)) 1 d2ý(r,t) r dr2 c2 dt2 s obecným řešením 0 (125) ^(r,í) = Í/(í-£)+%(í + £). (126) Na toto řešení se můžeme divat jako na rozbíhavou nebo sbíhavou kulovou vlnu. Tvar řešení také ukazuje, že rychlost šíření je c. V (lineárním) materiálovém prostředí se nahrazuje e0 —► ere0 a /x0 —► prpo- Pak dostaneme z rovnice (122) rychlost šíření c/n, když zavedeme index lomu TI — -\J ďp jJjf. 20 9.3 Obecné řešení nehomogenní rovnice pro potenciály. Pro obecné řešení rovnice (122) ještě chybí partikulární řešení nehomogenní rovnice. Zavedeme Fourierovu transformaci časové závislosti potenciálu a hustoty náboje 1 r°° &t) = 7T du4>(x,u)e-^, (127) Z7T J-oo 1 r°° p(£,t) = — áujp(x,uj)e-lu]t (128) a dosadíme do rovnice (122), áu A + — c/)(x, u) e-wt =--/ duj p(x, uj) e-tujt. (129) d \ cz J e0 J-oo Exponenciální funkce s různými co jsou nezávislé, proto platí (A + £) 4.(2,,) = -«|^>. (130) Hledáme Greenovu funkci diferenciálního operátoru na levé straně, definovanou vztahem (A + k2) G(x, x", k) = ô3{x - x"). (131) Řešení této rovnice závisí jen na absolutní hodnotu r = \x — x'\: péiikr G(x,x",k) = G(r,k) ==-. (132) Aur Důkaz: 1) (A + fc2)— =[~r + k2)^- = 0 (133) j-dr2 platí pro r / 0. 2) Násobíme levou stranu testovací funkcí f{x) a integrujeme přes celý prostor. Díky (133) stačí integrál přes infinitesimální kouli o poloměru e kolem počátku, d3x f(x) (A + A-2)-= / d3x f{x) (A + k2)-. (134) r Jr<€ r Rozvijeme exponenciální funkci a využijeme A 1/r = —47r<53(:? — x') d°xf{x) (A + k") [ 1 " 1'° ■ ^ ľ r.,^í,l . ,,1 /r = -l* ieyr (140) An\x-x'\ An\x-x'\ Potenciál časově závislého rozložení náboje vypapdá formálně jako statický potenciál, s rozdílem, že časový argument závisí na vzdálenost \x — x'\, t. j. potenciál a elektrické pole se šíří konečnou rychlostí. Takto získaný potenciál se nazývá retardovaný potenciál, analogicky to platí pro vektorový potenciál. 0ret(f,í) =-- / d3x'ľK ' _„c \ (141) 47re0 J \x — x Aretfrt) = d3x' nS':!~lZfA\ (142) 47T J \X — X \ 9.4 Pole časově proměnného dipólu Uvažujme bodové náboje, soustředěné všechny kolem počátku souřadnic. Pak můžeme pro vektorový potenciál psát [J=t,. r\Azm, _ r„ „-? (4- r A&t] " S I 'V-' " č J dV = ^eava[t--) (143) neboli pomocí dipólového momentu Ä(x,t) « £L f>(t - r-j , p{t) = J2eaSa{t). (144) Skalární potenciál spočteme integrací vztahu (Lorenzova kalibrace) -c2divA (145) 90 2 dt 22 Jednoduchými úpravami dostaneme Po div A --rot Ä-- Skalární potenciál je tedy 4>(x, t Pro intenzity dostaneme E(x, t) 4irr3 Po 4ti- r 3 x x x 1 x 47rcn r3 H<--) + -p(í-- f\*--) + -p U--c Aneor3 [r2 p[t--) -x x-p [t- - + — c/ ť p [t--l XX x : B(x,t) p [t--x x, p [t-- 4nr3 Dostatečně daleko od dipólu máme É(x,t) ---- D [x,t--] x n, Ane0cz r V c, B(x,t) = —-D (x,t--4ncr V c kde jsme označili D [ x, t — - j = p (t — - ) x n, x Pro hustotu energie a Poyntingúv vektor platí 1 1 jjlq J 16n2c4e0r2 w = \ Le- + ^B> ■D' S Ex B 1 Po lQn2c3e0 r2 Je přirozeně S_ W cn. Příklad: Vezměme rozložení proudu ve tvaru / 71Z \ j(x,t) = I6(x)6(y) sin í — J cos(^í)ě; Podle (143) a (144) spočteme snadno 0 < z < L. p (t) =- sm(cijŕ) ez a podle (150) D [ x, t - - tylo 2LIuj n ( r --sinů sin oj [t-- 7t V c 23 Příklad: Rotující náboj v rovině (x,y) v dipólové aproximaci: P = Qro (cos sm 0), D qr-quj {z sin ujt,—z cos ut, y cos ut — x sin ut). (156) (157) Časově středovaný Poyntingův vektor popisuje výkon záření do elementu prostorového úhlu, q2rQĹj4 1 + cos2,$ (x, y, z) ^ Z 167r2c3e0r2 2 Celková vyzařovaná energie za sekundu je pak p= (S)-ňds q2rQĹú4 q2a2 (158) (159) kde a označuje zrychlení r0uj2 na kruhové dráze. Táto ztráta energie by vedla k rychlému kolapsu atomů v klasické elektrodynamice. 9.5 Liénardův - Wiechertův potenciál Ať se nabitá částice pohybuje po zadané trajektorii x = x*o(t). Hustota náboje je pak p(x,t) = e5{3\x-xQ(t)). (160) Vzorec pro skalární potenciál přepíšeme jako ťU ť-t+V-* ')dťdV 47T60 J \X — x'\ 47T60 J R(ť) kde jsme označili R{ť) = x — x$(ť), R(t') = \R(t')\. S pomocí vztahu (161) ô ť -t R(ť)\ 5{t'-tr) R(tr)-v(tr) R(tr (162) cR(tr) napíšeme výraz pro skalární potenciál jako (163) Výraz pro vektorový potenciál je pak obdobně (164) 24 Vezměme teď jednoduchý případ pohybu s konstantní rychlostí podél minku pro nalezení časového zpoždění přepíšeme na c2(t - trf = {x- vtrf + y2 + z2, odkud r-\ r.r 1 1 - - )tr = t - - c" C c (x - vty + i - - Ui/ + Jmenovatel výrazů (163) a (164) pro potenciály můžeme psát jako V(X - Vtr) c(t — tr)--= C Po malé úpravě pak dostáváme 0(:?,í) vx 1 - -7 Ur 47T6, 0 Wl - v- r 2 pro skalární potenciál a l(í,r) = i. Ur'./mi. (n. pro vektorový potenciál, kde jsme označili Ax(%, t) — — Vektor intenzity elektrického pole je E{x,t) 1 1 47re0 j\_vt r 2 ,y * ^ 3 2 (x - vt,y,z) a vektor indukce magnetického pole je Pro vektor hustoty impulsu pole G = enE x b dostáváme e2/in 1 v G(x,t) 16tt2 1 - ^ r a pro hustotu energie ty = [e^e2 + b2/jjq) ji výraz 1 (x - t;í)2 + (l + S) (ž/2 + ty(f,r) 327T2e0 1 , * 6 25 10 Základy teorie relativity 10.1 Principy Princip relativity: všechny přírodní zákony jsou stejné ve všech inerciálních souřadných soustavách. Inerciální soustavy jsou takové, kde se pohyb dějě s konstantní rychlostí. Interakce částic se v nerelativistcké mechanice popisuje pomocí interakční potenciální energie, která je funkcí polohy interagujících částic. Tento způsob popisu v sobě obsahuje předpoklad o okamžitém působení. Princip konečné rychlosti šíření signálu: Rychlost šíření interakce je konečná. Z principu relativity je tato rychlost ve všech inerciálních soustavách stejná. Z Maxwell-ových rovnic je vidět, že jde o rychlost světla ve vakuu c = 299 792 458 ms"1. (175) Toto je exaktní hodnota, určující tak délkovou jednotku jednotkou času. Sjednocení principu relativity s principem konečné rychlosti šíření signálu je nazýváno Einsteinovým principem relativity. 10.2 Interval, vlastní čas. Uvažujme dvě události: emisi a absorpci fotonu. V soustavě K je (x2 - Xlf + (y2 - yif + (z2 - Zlf - c2(t2 - h)2 = 0, (176) v soustavě K' pak {x>2 - x',)2 + (y>2 - y[)2 + (4 - z[)2 - c2{ť2 - ťx)2 = 0. (177) Zavedeme obecně kvadrát intervalu mezi dvěma událostmi jako s\2 = c2(t2 - h)2 - (x2 - Xl)2 - (y2 - Vl)2 - (z2 - Zl)2, (178) popřípadě pro infinitesimálně blízké události ds2 = c2dr2 — dx2 — dy2 — dz2. (179) Je-li interval roven nule v nějaké inerciální souřadné soustavě K, je roven nule i v libovolné jiné soustavě K'. Potom tedy musí být ds2 = k{v) ds'2. (180) Vzhledem k homogenitě prostoru a času nemůže faktor úměrnosti záviset na souřadnicích, vzhledem k isotropii prostoru může pak tento faktor záviset pouze na velikosti relativní rychlosti uvažovaných inerciálních soustav. Uvažujeme-li tři soustavy K, K\ a K2, dostáváme ds2 = k(vi) ds2, ds2 = k(v2) ds2,, ds2 = k(ví2) ds2, =3- = ^^12^' (181) 26 a protože levá strana poslední rovnice nezávisí na úhlu mezi vektory rychlostí v\ a v2, zatímco pravá strana může, musí být k{v) = 1. (182) Kvadrát intervalu mezi dvěma událostmi (178) nebo mezi dvěma infinitesimálně blízkými událostmi (179) je stejný ve všech inerciálních soustavách. V předešlých úvahách se připojuje čas přirozeným způsobem k prostoru, proto je výhodně definovat čtyřrozměrný prostoročas či Minkowskiho prostor, v němž se počítá kvadrát prostorového intervalu záporně a kvadrát časového intervalu kladně (nebo opačně). Oznámení událost má význam bodu v čtyřrozměrném prostoručase. Označme si v soustavě K tu = Í2-Í1, 4 = (x2-x1)2+(y2-yi)2 + (z2-z1)2 s\2 = c2t\2-í\2. (183) Zkoumejme, existuje-li taková soustava K', kde by se obě události odehrály v jednom bodě prostoru, tedy že platí £'12 = 0. Máme tak podmínku s\2 = c2t\2 — í\2 = c2ť22 > 0,; takový interval se nazývá časupodobný. Naopak požadavek na to, aby existovala soustava, ve které obě události nastanou současně (t'l2 = 0), vede k podmínce s\2 = c2t\2 — í\2 = —£'12 < 0,; interval se pak nazývá prostorupodobný. V soustavě, která se pohybuje s daným hmotným bodem (£'12 = 0), můžeme tedy definovat vlastní čas jako 1 1 rS2 ŕ2 ( v2\2 ť2-t[ = - ds= 1 - — dt. (184) c Jsi Jíi > n2 V případě konstantní rychlosti v dostaneme jednoduchý vztah mezi parametrem s trajektorie tělesa a časovou souřadnicí t, v2 s2-s1 = c^l--(t2-t1). (185) 10.3 Lorentzova transformace Soustava K se pohybuje vůči inerciální soustavě K' rychlostí v podél osy x. Z elementárních úvah je zřejmé, že čtverec intervalu s2 = c2t2 — x2 se nezmění při transformaci ct = x' sinh-í/> + cť cosh-í/>, x = x' cosh-í/> + cť sinh-í/>, (186) podobně jako se nezmění čtverec vzdálenosti l2 = x2 + y2 při transformaci x = x' cos ip + y' sin 99, y = — x'simp + y'cos ip. (187) Pro počátek soustavy K' (bod x' = 0) máme v soustavě K z definice x/t = v, jednak z (186) x/t = ctanh-í/>, máme tedy tanh-í/> = v/c a vztah (186) můžeme zapsat jako Lorentzovu transformaci (188) 27 Vždy jsou uváděny dva klasické příklady na použití vztahu (188). (a) V soustavě K je podél osy x v klidu měřítko, jehož dvě rysky mají v této soustavě souřadnice xl7 x2. Vzdálenost (klidová) rysek je tedy Ax0 = x2 — x\. Vzdálenost v soustavě K' (souřadnice jsou určovány ve stejném čase t[ = ť2) je Ax = x'2 — x[ = Axoyl — ^2. Mluvíme o kontrakci délky. (b) V soustavě K' se v časech t[ a ť2 odehrají dvě události v jediném místě x[ = x2, y[ = y'2, z[ = z'2 (interval mezi událostmi je tedy Aí0 = ť2 —t[). V soustavě K je interval mezi těmito událostmi At = t2 —1\ = At0 j\J\ — ^ ■ Mluvíme pak o dilataci času. Vztah (188) IIlŮZGIIlG ZcipSclt i v diferenciálním tvaru cdť + ^áx' cdt =-c-- dx' + v dť . . dx = — , dy = dy , dz = dz . (189) Pro transformaci složek vektoru rychlosti (w = dx/dt, w' = dx'/dť) dostaneme z (189) vztah w'x + v w'x v ' w'Jl - K ^2 w. w'Jl - K r2 Sledujeme-li šíření světelného paprsku v rovině (wx = ccosů, wy = c sinů, wz w'x = ccosů', w' = c sinů', w'z = 0), dostaneme vztah (aberace světla) sin-ŕ? cos ů' sin ů'. (190) 0 resp. (191) Pro «/cC 1 položíme ů = ů' — Aů a porovnáním nejnižšího členu Taylorova rozvoje dostaneme obvykle uváděný vztah v Aů sinů'. (192) 10.4 Ctyřvektory, čtyřtenzory Nejprve definujeme podstatné tenzory. Metrický tenzor v Minkowskiho prostoru a jednotkový tenzor jsou 9ik = 9 ik 0 o f 1 o 0 -1 0 0-1 V o o o 0 \ o o -1J ^ 1 0 0 0 ^ 0 10 0 0 0 10 V o o o i y (193) Qik se nazývá kovariantní metrika, inversní metrika gtk se nazývá kontravariantní. Dále definujeme kontravariantní a kovariantní úplně antisymetrický tenzor 4. řádu pomocí vztahů Jklm / 0123 ^iklm (e0123 -1). (194) 28 Čtyřvektor souřadnic události (kontravariantní a kovariantní) zapisujeme jako xL = (x°, x1, x2, x3) = (ct,x), Xi = (x0,x1,x2,x3) = (ct,—x). (195) Metrika nám udává invariantní prostoročasovou „délku" vektoru x% s2 = gik x1 xk = gik Xi xk = xl Xi = c2t2 - (x2 + y2 + z2). (196) Přitom platí Xi=gikxk, x^g^Xk (197) (zvednutí a spuštění indexů = transformace mezi kontravariantními a kovariantními indexy). V čtyřrozměrném zápisu můžeme Lorentzovu transformaci (ve směru x) (188) psát ve tvaru xi = Aikx'h (198) s Lorentzovou maticí A* ( 1 ll 0 0 \ -cl 7 0 0 0 0 10 V o o o i ) (199) kde 7 je běžné zkrácení 7 := ,-r- (200) 10.5 Ctyřrychlost a čtyřzrychlení Definujeme čtyřvektor rychlosti přirozeným způsobem jako derivaci čtyřvektoru událostí, ze kterých se skládá světočára (čtyřrozměrná trajektorie) tělesa, podle parametru s (= c krát vlastní čas) «* = ^, U* = I ~T=^, r—, 1 , lťut = l. (201) ds \ A /1_y* fl u% je tedy tečným vektorem světočáry. Obdobně čtyřvektor zrychlení dúl d2x% - d, w "'"• = a (202) Podívejme se na relativistický popis pohybu s konstantním zrychlením. V souřadné soustavě, kde rychlost částice je momentálně nulová (v = 0), máme < = (1,0,0,0), ďK=(0,-,0,0), (203) kde a je obyčejné zrychlení. V souřadné soustavě pohybující se rychlostí v ve směru x je rychlost a zrychlení l v \ ( v_ Av_ dv ul = I —;^=, . :,0,0 , é =\—c3dt 9,-&-o,0,0|. (204) 1 T^?'cx/T^?' j i-a2Wi--^2'' 1 29 Po malé úpravě (z rovnosti ) dostáváme d / v dt \ . Ix -2 = a. (205) S počátečními podmínkami vq = 0, xq = 0 dostáváme řešení pro konstantní zrychlení a v= / ^ = -h/1+(-) "H- (206) ' i /at aŕ cz / /_ fáty \ s2' a \ V " ' W ' 10.6 Relativistický impuls Jak v klasické mechanice, tak existuje i ve speciální teorii relativity princip nejmenšího účinku. Jako invariantní a jednoduchý účinek bodové částice se nabízí integrál délky podél světočáry. Z důvodu dimenze a abychom v nerelativistické limitě dostali pro účinek známý nerelativistický výraz, musíme konstantu úměrnost zvolit rovnu —mc, tedy S = -mc jb ds = -mc2 jí*" Jl - ^ dr. (207) Lagrangeova funkce a impuls jsou v2 ^ OL m v S volbou faktoru —mc dostaneme v přiblížení « tedy nerelativistickou Lagrangovu funkci volné částice minus konstantu mc2, která je bez významu pro pohybové rovnice klasické mechaniky. Hamiltonova funkce je pak 2 _ Tfl C I- H = p ■ v — L = . ■. = Jp2c2 + m2c4; (210) v případě «30 c 1íp'02 + v_ _p/12> 7 )p'12 + l F'02^ 0 p/32 7 íj?'03 + £ ^13 ^ \ 0 } Převedeno do vektorů intenzity a indukce platí Ey=1(E'y + vB'z), Ez = 1{E'z-vB'y) Fx = E'x, Bx = B'x, B,n ^\B'y--2E'zr 32 B? 1\B'Z + -E'y (230) (231) V nerelativistickém přiblížení (v/c —► 0) přechází (231) na Ě = Ě' - v x B' B = B'. (232) Invarianty pole můžeme zkonstruovat z tenzoru pole. Poněvadž je antisymetrický, zúžení nedává nic a máme až kvadratické výrazy im kn rp p p pi. y y 1 tk 1 mn ± %k ± inv, -e 2 ikmn F'ik F'mn Fjh -kF1, mv. (233) Duální tenzor vyjádřený pomocí intenzity elektrického pole a indukce magnetického pole má tvar " BX By ik ( o Bx By y B z 0 c Invarianty mají pak vyjádření FtkFtk Bz \ E, E, c 0 c E,,, E. c Ex c 0 c 4- E-B (234) (235) 12 Synchrotronové záření 12.1 Liénardův-Wiechertův potenciál Počítáme potenciál pole, vytvářeného jedním nábojem, který se pohybuje po trajektorii x = Xo(t), v čase t v bodě P(x, y, z). Potenciál je dán stavem pohybu částice v čase ť, pro který platí (doba potřebná pro šíření světelného signálu) c(t - ť) = R(ť) = \x-x0(ť)\. (236) V souřadné soustavě, ve které je částice v čase ť v klidu, máme právě Coulombův zákon cp(x,t) Ä(x, t) = 0. (237) 47T60 R(ť) ' Podmínku (236) zapíšeme ve čtyřrozměrném (kovariantním) tvaru jako podmínku toho, že interval mezi událostmi „emise fotonu" (cť,Xo(t')) a „absorpce fotonu" (cb,x) leží na světelném kuželu, tedy pro rozdíl čtyřvektorů událostí Rk = (c(r — ť),x — xo(t')) platí Rk Rk = 0. (238) 33 Pomocí tohoto nulového čtyřvektoru a jednotkového čtyřvektoru rychlosti částice u< = I -r^>-r=f) > v = v(ť) = ^1, ^ = 1 (239) se pokusíme zapsat čtyřvektor potenciálu pole tak, aby pro v = 0 (tj. pro čtyřvektor uL = (1,0)) přešel do tvaru (237). Z možných kombinací snadno nalezeme výsledek Aí=(*,ä)= 6 u[k. (240) Pokud nevypisujeme explicitně argumenty, musíme mít na pamětí, že levé strany vztahů jsou uvažovány v čase t, pravé strany v čase ť. V trojrozměrném značení pak má (240) tvar 0 = ---A-^š, A=6-^ * (241) 47T60 R (l _ J|)' 47TÄ(l_JgN Výsledek (241) je přirozeně stejný jako (163) a (164). Při výpočtu polí 3Ä -> E = -V0 - —, B = V x A (242) budeme potřebovat následující triky pro výpočet parciálních derivací: Derivováním vztahu (236) podle t dostáváme dR _dRdť _ Ř-vdť _ ( dť\ dť _ 1 lh~dťljt~ R~~dt ~C[1~~dt)^~dt~i-^É' ^2A3' v 7 1 cR Obdobně derivováním vztahu (236) podle x dostáváme f) R R R -cVť = VR{ť) = ^ Vť + £ Vť =--. . (244) Výraz pro potenciály ve (242) pak budeme chápat jako funkce f(x,ť), a budeme počítat parciální derivace podle x při konstantním ť a podle ť při konstantním x. Porovnáváním diferenciálů df = Vf-dx + ^dt = Vf-dx+^ dť =hf + % Vť) -dät+ZĹ^ dt (245) přepíšeme (242) jako - - ,, _ „ d(/)(x, ť) -> . dÄ(x, ť) dť E = -V(f)(x, ť) - Y\ ' ; Vť v ; - - -r/- ,x ^ , cM(í,ť) 5 = V x A(f, ť) + Vť x —i^-^ 34 Pro intenzitu elektrického pole dostáváme pak E 4?re0 l-v4 71 n x I I n — -c i x -u; i?2 1 c2i? 1 (247) zatímco pro indukci magnetického pole pole B -n x E = - C 47T (■u x rt) R2 1 3--h n x n x I I n — -c ] x -u; ci? 1 (24í Označili jsme jednotkový vektor n = R/R a zrychlení w = dv/dť. Limitní případy pro v/c —► 0 jsou en 47reni?2 5 efj,0(v x n) 4nR2 (249) 12.2 Intenzita záření Poyntingův vektor (energie, procházející jednotkovou plochou za jednotku času, dimenze [Jm_2s-1]) je Š = —Ě x B = e0cE2n (250) Po a intenzitu záření (tj. energie, vyzařovanou za sekundu do elementu prostorového úhlu, [Watt]) spočteme tedy jako dl = lim Š-nR2díl. (251) R^oo Po dosazení z (250) a (247) 2(rt • w)(v ■ w) dl 167r2eoC3 w (n • w)2 cl díl. (252) Pro v/c —► 0 dostáváme s označením n-w = w cos£ pro celkovou vyzařovanou intenzitu e2w2 2tt e2 -u;2 ifi 2 s/ dr] d£ sin£(l - cos2£) V klidové soustavě částice je tedy (s označením I = dE/dt) (253) dE e2w2 Ô7T 6qC3 dŕ, -u* drr* ds (1,0), w* dul ds 0, w (254) Relativisticky invariantní výraz (tj. diferenciál čtyřvektoru impulsu) vytvořený z čtyř-vektorů rychlosti a zrychlení, který v klidové soustavě přejde na výrazy ze vztahu (254), je pak ■E J, P" = [ —,P dp1 e2 duk duk 6ireoC ds ds 35 dx1 e2 duk duk 6ireoC ds ds u% ds. (255) V laboratorní soustavě tedy máme pro celkovou vyzařovanou intenzitu výraz p2 w2 — (w x 6ne0c3 (256) Zde jsme potřebovali vyjádření čtyřvektoru rychlosti i zrychlení v laboratorní soustavě. Abychom nemuseli při výpočtu čtyřvektoru zrychlení užit obecné Lorentzovy transfor- 2 —»// / -i 7^2 mace, vypočteme wl derivováním známého tvaru uL = yí/y 1 — fr,v/ \ c\Jl potom V homogenním magnetickém poli se nabitá částice pohybuje rychlostí v po kružnici poloměru R, její zrychlení w = v2/R je kolmé k rychlosti. Dosazením do vztahu (256) e2 v4 e2c ( p \4 e2c / T . , ( — 1 «-í-1 (258) 67re0c3 (i _ v^_y 6ne0R2 \mcj 6ne0R2 \mc2, V posledním výrazu ve (258) jsme použili aproximace vysokých energií, kde pro kinetickou energii platí T = \/p2c2 + m2c4 — nic2 « pc. Z tohoto výrazu je také zřejmé, že synchrotronové záření je omezujícím faktorem při urychlování lehkých částic (elektronů a positronů). Pro normovači hodnotu Rq 0, 5 km můžeme psát I « (i?0/Ä)2(T/mc2)4eVs-1. Jsou-li rychlost a zrychlení v určitém okamžiku rovnoběžné, dostáváme (n • v = vcosů, rychlost podél osy z) pro úhlové rozložení záření výraz e2w2 sin2-# ^ 16?r eoc3 fl - ^COS^ c Pro hodnoty v/c —► 1 má úhlové rozložení velmi úzké, ale „dvouhrbé" maximum kolem ů = 0. Jsou-li rychlost a zrychlení v určitém okamžiku navzájem kolmé, dostáváme (n ■ v = v cos ů, n • w = w cos p sin ů, rychlost podél osy z a zrychlení podél osy x) pro úhlové rozložení dl e2w2 16n2e0c3 1 - \) sin2 ů cos2 p (l-fcosí?)4 (l-fcosí?)6 díl (260) 13 Maxwellovy rovnice v čtyřrozměrné formulaci 13.1 Čtyřrozměrný vektor proudu, rovnice kontinuity Definujeme čtyřvektor proudu (pro částici: náboj krát čtyřrychlost) dxl f = p— = (cp,pv) = (cp,j). (261) 36 Náboj, který ubude v nějakém objemu, můžeme zapsat dvojím způsobem -ít.ípdV = fJ'-ridS- (262) S pomocí Gaufiovy věty pak z (262) plyne (263) tedy (objem je libovolný) rovnice kontinuity _± _ do df V • j+ — = — = 0. dt dxl (264) 13.2 Homogenní Maxwellovy rovnice Z vyjádření tensoru elektromagnetického pole pomocí potenciálu snadno odvodíme platnost vztahu dFik dFH dFu dxl dx'L dxk 0. (265) Na levé straně je úplně antisymetrický tensor třetího řádu, představuje pouze čtyři různé rovnice. Zřetelněji je to vidět, užijeme-li zápis duálního (pseudo)vektoru 2 dxk dxk (266) Nultá komponenta dává tvrzení o nezřídlovém charakteru magnetického pole, další tři komponenty Faradayův indukční zákon V -B = 0, (267) 13.3 Nehomogenní Maxwellovy rovnice Čtyřrozměrný zápis nehomogenních Maxwellových rovnic, obsahujících hustotu náboje a proudu, je q pík dxk -Poj ■ (265 Nultá komponenta je rovnice pro divergenci intenzity elektrického pole (Gaufiova věta elektrostatiky), zbývající tři pro rotaci magnetického pole (Ampěreův zákon) ^ ň 1 9E V x B = -^^- + PoJ-cz ot (269) 37 13.4 Tensor energie-impulsu Z hustoty energie W=\ (e0Ě2 + —B2] , (270) 1 \ ßo ) z Poyntingova vektoru Š=—ĚxB (271) ßo a z Maxwellova tensoru napětí aaP = e0EaEp + —BaBp - Wôap (272) elektromagnetického pole (ve vakuu) můžeme sestavit čtyřrozměrný tensor energie-impulsu, Tk=[ c . (273) \ -Oa. —CTa/3 ) Pomocí tensoru elektromagnetického pole dostáváme jednoduchý výraz Tík = — (-glmFaFkm + -aíkFlmFlm ßo V 4 (274) Tensor energie-impulsu soustavy částic zapíšeme pomocí analogie s tensorem energie-impulsu elektromagnetického pole. Hustotu impulsu soustavy částic napíšeme jako na = pcua, fi = J2ma^3)(x-xa). (275) a Hustota impulsu je u elektromagnetického pole rovna hustotě toku energie dělené c2. Výraz (275) bude tedy analogicky roven T0a/c. Veličina pc je nultou komponentou (stejně jako hustota náboje u čtyřvektoru proudu) čtyřvektoru toku hmoty pdxl/dt. Tensor energie-impulsu tak můžeme konečně psát jako ^ = = V-T?- (276) Pro tensor energie-impulsu elektromagnetického pole dostaneme s využitím Maxwell-ových rovnic ^r = -^ŕ, eiUm^ = 0 (277) ôx ôx výraz ^-kTFk = -FtkJk. (278) Pro tensor energie-impulsu soustavy částic dostaneme s využitím pohybových rovnic pc ^ = p Fk uk & pc ^ = Fk Jk (279) a rovnice zachování hmotnosti (rovnice kontinuity pro čtyřvektor toku hmotnosti, podobně jako pro čtyřvektor proudu) 0 ( dxk\ p — )= 0 (280) dxk \ dt výraz dxk Spojením (278) a (281) dostáváme zákon zachování d 9 Tk = rkJk^ (2glj n , ,T'f + TPk = 0. (282) oxk v ' Pro tensor energie-impulsu platí (rovnost právě pro elektromagnetické pole) T\ > 0. (283) 14 Elektromagnetické vlny 14.1 Vlnová rovnice Vezmeme nehomogenní Maxwellovy rovnic ve vakuu (p = 0, j = 0) a dosadíme vyjádření pole pomocí potenciálů dxk ' y dxi dxl) ' a 92Ak kl d2A* 9 -:--9 -= 0. dxi dxk dxk dxl Lorenzova kalibrační podmínka (121) nabývá formu čtyřdivergence dAk (284) dxk a zjednoduší (284) na vlnovu rovnici d2 A{ 0 (285) Pomocí d'Alembertova operátoru f w = °- <286) □-A-i^ (287) máme pak ve třírozměrném zápisu 1 dcf) , V-A = 0, □ = 0, LU = 0. (288) c2 ot 39 Vlnové rovnice, spolu s Lorenzovou kalibrační podmínkou, jsou ekvivalentní Maxwell-ovým rovnicím pro volné elektromagnetické pole. Konsistence kalibračních podmínek s rovnicemi pole se dokáže takto: 1) Zvolíme Lorenzovu kalibraci dkAk = 0 na počáteční nadploše t = 0. 2) Řešíme rovnici D Ak = 0. 3) Protože Ak je řešení vlnové rovnice, platí jt dkAk = Ä0 + VÄ= AA0 + VÄ=V (VA0 + AJ = -VÉ = 0. dkAk = 0 a í dkAk = 0 dt pro t = 0. 4) Když Ak je řešení vlnové rovnice, pak je dkAk také řešení: UdkAk = dkUAk = 0. 5) Z toho vyplývá, že dkAk = 0 všude. 14.2 Rovinná monochromatická vlna Řešení hledáme ve tvaru rovinné vlny, tedy konstantní čtyřvektor násobený komplexním fázovým faktorem A* = Re{ať exp(ikjXj)}, h k' = 0, h a' = 0. (289) Poslední vztah ve (289) je dán Lorenzovou kalibrační podmínkou. Vlnový čtyřvektor zapisujeme jako k%=^-,kj, k = ~ň' n=l. (290) Velmi jednoduše popíšeme pomocí charakteristik rovinné monochromatické vlny Dop-plerův jev. Mějme zdroj světla, který je v klidu v soustavě Kq. Soustava Kq se pohybuje vzhledem k laboratorní soustavě K rychlostí v. Ať je úhel mezi směrem pohybu zdroje a směrem šíření světla a. Potom platí 7^0 _ V 1.1 L.0 _ ^ C ^ jLO _ 7,0 _ ^ ^(0) - r, u) "(o) - r,—v£' (o) „2 c ,i k1 -^k° uj(0) u kfo = —i 2 , kfo = —^ cosqí(o), k = — cos a. (291) a odtud r .o " = "(o) -iv (292) 1 — - cos a c Pro rychlosti malé ve srovnání s rychlostí světla máme v lv 2 uú cj(o) 1 H— cos a H---- cos 2a . (293) c ' 2 c2 40 Tensor energie-impulsu je i* + Re {crdi exp (2ikj xj^ j . (294) r2 1 Tk = — W k'kk, W =- uj2 2/xq a a- Ve střední hodnotě podle času je druhý člen ve výrazu pro hustotu energie roven nule. Oba invarianty (235) jsou rovny nule. Se speciální volbou kalibrace (spojené ovšem s jednou určitou inerciální souřadnou soustavou) máme pro rovinnou vlnu ve směru x A1 = (0, A), A = ay cos(ut — kx + ot)ey + az sin(cIqD- (313) k 15 Rozptyl záření volnými náboji 15.1 Thomsonův vzorec Zavedeme pojem účinného průřezu. Ať dl značí intenzitu záření, tj. střední hodnotu energie vyzařované soustavou za jednotku času do elementu prostorového úhlu díl a S je střední hodnota Poyntingova vektoru (střední hodnota toku energie) dopadajícího záření. Potom je diferenciální účinný průřez (účinný průřez rozptylu do elementu prostorového úhlu díl) veličina rozměru elementu plochy dl , . da = (314) Uvažujme teď rozptyl elektromagnetické vlny jedním jednotkovým volným nábojem. Budeme předpokládat, že rychlost získaná nábojem bude malá a že vlnová délka je mnohem větší než amplituda vyvolaných kmitů náboje okolo původní polohy (kam umístíme počátek souřadnic), tedy můžeme psát m——r = eĚo cos (k ■ x — ut + a) « cĚq cosíut — a). (315) dtz v ' Pro intenzitu dipólového záření kmitajícího náboje máme podle (151) ve směru ň 4 4 dl = _ 26 2 AĚQ x n|2cos2(^t - a) díl = qo 6 E20 sm2ůdQ (316) 167Tze0mzcó 62'kzeQmzci a pro střední hodnotu Poyntingova vektoru dopadající vlny S = ce0 E0 cos2(cjr — a) = — ce0 E0, (317) takže pro diferenciální účinný průřez je -,2 2 do- = I---t J sin2?? díl (31* \47reo mcz 43 Celkový účinný průřez je pak dán Thomsonovým vzorcem 2 \ 2 3 \4ne0mc2j 3 e' (7 = ^f^_7| =-nr2e, (319) kde re je klasický poloměr elektronu. 15.2 Modifikace Thomsonova vzorce Uvažujme nyní nikoliv volný náboj, ale tlumený oscilátor, tedy d2x dx dť2 1 dt u m Pro dipólový moment p = ex odsud dostáváme + 7——h uj0 x = — Eq COS Uút. (320) e2 (ĺJq — co2) cos uot + 7cj sin ut m (cJq — u2)2 + 72cj2 P=—-T~2-, -.9.. .9.- ^0- (321) Celkový účinný průřez je v tomto případě Š7T 0 cj4 o = —K —2- . (322) 3 (cjg - cj2)2 + 72w2 16 Index lomu Definujeme polarizovatelnost ol(uj) jako konstantu úměrnosti ve vztahu mezi (lokálním) elektrickým polem E\oc a dipólovým momentem p. Vyjdeme z komplexního zápisu (320) d. d.iř/ o_. g * c * , , , , 377 + T-TT + x = — ^loc = — E0 exp(-iut). (323) dŕ^ dŕ m m Potom e2 1 p = e0a(uj)Eioc, a(uj) =--5-:--. (324) e0m ĺJq — vyui — uz Polarizace je pak P = N p. Musíme ovšem uvážit, jaké pole působí na náboj. Připomeňme z elektrostatiky, že je-li v dielektriku s homogenním polem dutina, je lokálni pole rovno Éloc = É, Éloc = É + -P, Éloc = É + ^-P, (325) eo 3e0 podle toho, jde-li o štěrbinu podél nebo napříč pole nebo o kulovou dutinu. (V případě štěrbiny napříč pole máme E\oc = j^D.) Pro úplnost poznamenejme, že pro magnetické pole máme v podobné situaci Bloc = B-M, Bloc = B, Bloc = B-^M. (326) 44 Pro dielektrika uvažujeme o vázaných nábojích uvnitř kulové dutiny, můžeme tedy psát P = e0 NaĚloc = eoNa(Ě+^-py (327) tak že Na Index lomu dostaneme z relace e0E e0E (za velmi častého předpokladu h(uj) = po): Na •> D P / \ n2 = er = — = ! + — (329) "2 = 1 + T^ps- (330) Obvyklá forma tohoto vztahu je (Clausius - Mosotti) n2 - 1 3^— = Na. (331) nz + 2 Ve vodiči uvažujeme o téměř volných elektronech (nevázaných k atomu, tedy = 0) a dále máme pro konstantu 7 (ze dvou různých vyjádření proudu a zápisu změny impulsu za dobu mezi srážkami) Ne2 j = aE, j = NeVd, mvďy = eE =>• 7 =-. (332) m a (vd je zprůměrovaná rychlost elektronů - drift.) Také lokální pole je rovno vnějšímu, opět díky neustálému pohybu téměř volných elektronů. Odtud máme pro index lomu úl „Ne2 = 1 " -^-^^ < = —■ (333) ujp je tzv. plasmová frekvence. uj2 + íuj uj2 a ' p m en p O" J V plasmatu je 7 zanedbatelné, t. zn. e0/a —► 0, a kvadrát indexu lomu je n2 = 1 - % (334) Když > ujp, je n reálné a plasma propustí e.m. vlny, když uj < ujp, je n imaginárni, což znamená, že plasma odrazuje vlny. Krátké a dlouhé radiové vlny se odrazují od ionosféry a vrátí se k zemi, ultrakrátké propagují do prostoru. Kovy jsou průhledné pro ultrafialové světlo. Při vstupu vesmírních lodí do atmosféry se zahřívá vyduch, tím se zvyšuje počet iontů, t. j. plasmová frekvence v okolí, a komunikace je přechodně přerušena. 45