Elektrodynamika 1 Elektrické a magnetické veličiny, jednotky SI Elektricky proud I je v systému SI základní veličina, jednotka je 1 Ampere (1 A). Definice: Stejne proudy ve 2 rovnobežných drátech ve vzdálenosti 1 m mají velikost 1 A, když vzájemna pritahovací síla na 1 m dratu 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: Pri elektrickem proudu o 1 A proteče prurežem vodiče 1 C ža sekundu. 1 C = 1 As ~ 6 • 1018 elementárních naboju. CoulombUv zákon Pritažliva nebo odpudivá síla dvou nýboju q1 a q2 v místech x1 a X2 _ F = k^Y2 (x\ — x2), F = k^Y2, r12 r12 kde r12 = \X1 — X2\ a k je konstanta. Dimenze [k] teto konstanty výplyva ž (1). AT rilAs • As Nm2 N=[k]-— [k] = TY2 ■ m2 A2s2 Použitím rychlosti svetla c platí pro experimentalne urcenou konstantu k = 10-7c2— ■ A2 V SI se píse ž historickych duvodu k =: 1 4neo' kde e0 je tzv. „dielektricka konstanta vakua". (1) Požnamka: Mohli bychom predpoklídat k = 1. Tím bychom urcovali dimenži naboje, vyjídrenou mechanickymi velicinami. V systemu cgs dostavíme ž Coulombova žaíkona , gcm [q]2 i, _1 dyn =—— = —2, [q] = g2 cm2 s ■ s2 cm2 To je jednotka naboje v Gaufiove systemu. Elektrické pole: Síla souboru naboju q1}. ■ ■ ,qN v místech X1, ■ ■ ■, x n na dalsí naboj q v bode X je popsaní pusobením elektrickeho pole v bode X na naboj. Intenzita eletrického pole E(X) je lokílní vlastnost prostoru, v nemž se nachažejí elektickíe naíboje. F (X) = j—y. q q* t^ts =: q E (X)■ (2) 4neo i=1 \x — Xi\s 1 Intenzita = síla na jednotkový náboj. Dimenze: [E] — N — J c = mc ■ Napětí U. Prýce — zmena potenciální energie při pohybu nýboje v elektrickém poli je daná integrýlem W — j F ds. V prípade pohybu o vzdalenost l podel homogenního elektrickeho pole (např. v deskovem kondensatoru) platí jednoduše W — qEl —: qU. Električke nýpetí je rozdíl potencialní energie jednotkoveho naboje na dvou bodech. Dimenze: [U] — ^ Jednotka: 1 Volt — 1V — 1 ±. Volt se pouzívý pro bezne oznýcení jednotky elektrickeho pole, [E] — 1 m. Magnětickě pole. (Stacionírní) elektricke proudy vyvolavají síly na pohybující se naboje (Lorentzova síla). Analogicky k zavedení intenzity elektrickeho pole uvazujeme sílu souboru eletřických proudu v danem usporadíní vodicu na usek dl jednoho dalsího drítu s konstantním proudem I. Síla je umerní I a dl, dF (x) — I dl x B (x). (3) Magnetická indukce B (x) zahrnuje pusobení vsech elektickych proudu v bode x. Dimenze: Z (3) vyplyva N J Js Vs N — Am[B], [B] — — — — — — — —. Am Am2 Cm2 m2 Jednotka: 1 Tesla — 1T — 1 . Prispevek elektrickeho proudu v íseku drítu dl na míste x2 k magnetickemu poli v bode xi (analogicky Coulombovu zakonu) je dB (xi) — km I df x .f"1 ~ x' . (4) \x1 — x2\3 V SI se píse km — ^, fi0 — -0^? je tzv. „magneticka permeabilita vakua". Síla mezi dvema infinitesimalními íseky vodice, umísteními v bodech s elek- trickími proudy I1 a I2 (analogon Coulombova zakona): 4n \ \xi — x2\3 J 2 Hustota náboje: p(x) = lim ' v / AT/" .f Q AV^0 A V kde q je naboj v objemu AV, který obsahuje bod x. Hustota proudu: V (x) = Ä XÄ dl' kde AA je element prurezu vodiče a ^ je jednotkový vektor ve smeru elektrickeho proudu I. 2 Maxwellovy rovnice, statický prípad J. C. Maxwell nasel v 19. století dve vektorove a dve skalýrní parciální diferenciální rovnice 1. radu, které spojují električke a magneticke pole navzajem a s elektrickým nabojem. Z tech rovnic lze odvodit cela elektrodynamika. (5) Vektorove operítorý div a rot lze výjadrit pomoci operatoru Nabla, V ío_ d 1\ \dx' ôy ôzj div v = V • v, rot v = V x v. (6) (7) V prípade statickích polí, kdýž casove derivace jsou nulove, odpojují se elektricke a magnetickíe pole. Elektrostatika Zakladní rovnice divV = P, rotE = 0. (8) eo Elektrostaticke pole je bezvírove pole se zrídlem p. Pro takove pole existuje skalarní potenciíl 0 tak ze E = -grad 0 = -V(. (9) Dosazením do Maxwellový rovnice dostaneme —div grad 0 = —V20 = P. 3 Definujeme Laplaceův operátor A- V2 = — — — (10) ' dx2 dy2 d z2' nabýva rovnice pro potenciál následující tvar, P I (11) A0 = — P en Tato rovnice se nazáva Poissonova rovnice. Oproti Maxwellove rovnici divE = —p ma výhodu, ze jejím rešením je jedna skalární funkce 0(x) místo vektoroveho pole E (x). Magnetostatika rotB = tm?, divB = 0. (12) Magnetostaticke pole je bezzrídlove, silocarý jsou uzavrene; existuje vektorový potenciál A, tak Ze B = rot?, (13) protože div rot A = 0. Dosazená do Maxwellový rovnice vede k rotB = rot rotA = grad divA — AA = ttnj, nebo AA — V (VA) = —tm?. (14) Potencialý nejsou jednoznacne urcený svámi rovnicemi (11) a (14), 0 je jednoznacne až na libovolnou konstantu, A? az na gradient libovolne funkce, ponevadz rot grad/ = 0. Potenciálý lze kalibrovat, jsou kalibracná, nikoliv fýzikalná veliciný. 3 Greenova funkce Poissonovy rovnice, ô -distribuce Poissonova rovnice (11) je eliptickaá lineáarná parciáalná diferenciáalná rovnice druháeho radu. Pro takove rovnice existuje rozsahlá teorie, mý budeme ale zkonstruovat resem na záakláad e znáamýách fýzikáalnáích skute cnostáí. Elektricke pole bodoveho naboje v pocatku je E = 7^4., (15) 4nen \x\- K ' p ráíslu snýá potenciáal je 0 =T^-- (16) 4nen r 4 Pro hustotu bodového náboje platí p{x) = 0 y x = 0 a J d3xp(x) = q. Přestože taková funkce v klasickém smyslu neexistuje, existuje matematický objekt, který popisuje idealižaci bodoveho naboje a pomoci toho můžeme skladat električke pole a potencial libovolneho rozložení naboje ž výražu (15) a (16). Hustota bodoveho nýboje je žkonstruovana pomocí „Diracovy 5-funkce", který lže definovat jako limita funkci' s integrýlem od minus do plus nekonecna rovným jedne, ktere jsou „vľce a vľce soustredene na jeden bod", napr. 1 1 _*2 5(x) := limge(c) = —^= lim-e e . (17) Takove limity majý vlastnost, že /OO ŕ'CO f (x) 5(x)dx = hm/ f (x) gĚ (x) dx = f (0), (18) -OO e >0 J — OO nebo obecneji, /OO f (x') 5 (x - x')dx' = f (x) (19) O pro libovolnou spojitou funkci f. Integral s 5-funkcľ tedy vybŕra hodnotu funkce v urcitem bode. „Výber urcite funkcný hodnoty" je exaktný výžnam 5 v rýmci distribucý t.j. jako žo-bražený vhodneho prostoru funkcý jako nekonecne diferencovatelne funkce s omeženmi nosicem, do realných nebo komplexných ďseL V jažyce distribucý se pýse ekvivalent rovnice (19) 5x[f] = f (x). (20) Distribuce 5x priražuje funkci f funkcní hodnotu f (x); pro její definici stací, že f je spojitýa funkce. Hustota bodoveho naboje v pocýtku souradne soustavy se vyjýdruje trojrožmernou 5-funkcý p(x) = q53(x) = q5(x) 5(y) 5(z). (21) Dosadýme potenciýl a hustotu nýboje v bode x' do Poissonovy rovnice, A0(x) = -S- A—^— = -q 53(x - x'), (22) 47re0 \x — x '\ e0 odvodýme ž toho A = 53(x - x'). (23) 4n \x - x'\ G(x,x '):=—-^— (24) se nažýva Greenova funkce Laplaceova operatoru (=potenciýl jednoho bodoveho naboje.) 5 Poněvadž Poissonova rovnice je lineární, můžeme zkonstruovat potenciál rozložení náboje p(x) jako superpozici potenciálů bodových nábojů, tedy ((x) = -L /d3x'p(x')—^, (25) poprípáde ((x) = — íd3x' G(x, x') P^. (26) J eo To je resení Poissonovy rovnice pro rozložení náboje p(x). Důkáz: —Ax0(x) = /d3x AxG(x, x') ^ = /d3x 53(x — x') ^ = p(X). Index x ždůrážnuje působení Lápláciánu ná necárkováne sourádnice. Greenová fůnkce není jednožnácná. G(x,x') plůs libovolne resení homogenne diferenciální rovnice A( = 0 je táke Greenovoů fůnkcí stejneho operátorů. Greenova věta, Greenův vzorec Odectením dvoů identit, V (uVv) = u Av + Vu V v á V(vVu) = v Au + Vv Vu, dostáneme V (u Vv — vVu) = u Av — v Au. Integrováním á áplikácí Gáůfiovy vety odvodíme Greenovů vetů: / (u Av — v Au) dV = (u V v — v Vu)n dS, Jy Jav (27) kde V je objem á n je normálni vektor ná okráj dV. Podle okrájovích podmínek má Greenová vetá dve stándárdní áplikáce. Dirichletův problem: Zníme p(x) ůvnitr nejákeho objemů V á potenciíl ( ná okráji dV. Dosádíme u = ( á v = G do Greenová vžorce. / [((x') Ax>G(x, x') — G(x, x') Ax>((x')] d3x' = Jv = I \(f>(x') V x' G(x, x') — G(x, x') Vx> ((x')] n dS'. av Podle predpokládů je první clen ná práve stráne žnámí, drůhy, kterí obsáhůje grádi-ent potenciálů, nikoliv. Ták vyůžíme možnost volby Greenovy fůnkce á žkonstrůůjeme 6 takovou, která se rovná nule na okraji dV. Greenově funkci s tou vlastností říkáme Dirichletovu Greenovu funkci, GD. Potenciál ve V je pak dán vzorcem 0(f) = - f G d (x, x') dV + / VGD (x, x') <(x') n dS'. (28) JV t0 JdV Neumannův problém: Známe gradient potenciálu na dV. Vybereme Neumannovu Greenovu funkci GN, jejíž gradient se rovná nule na okraji. <(x) = -f GN (x,x') p(x^ d3x' -f GN (x,x') V
(x). Energie dvou naboju = energie naboje qi v bode x1 v potencialu vyvolanem nabojem q2 v bode x2 plus váraž s prehoženámi naboji lomeno dvema, aby se energie nepočítala dvakrat. 1 U = 2(q1 <2(x1) + q2 01(x2)), kde 1 4>i(x) q 4nt0 \x — x, Energie spojiteho rozložení náboje: 1 U =2| p< dV. (30) Pro rozložení bodovách naboju je potencial 1 v- qb 9((xa) = ^ a energie 0(x«) = ~.- V]-, rab = \Xa — £b\ (31) 4vreo rab U = — T, —. (32) 5 MultipOlový rozklad pole V nasledujícím budeme hledat rožvoj elektrickeho potenciálu ve velke vždalenosti od ždroje. 7 5.1 Laplaceova rovnice ve sférických souřadnicích Ve sférických souřadnicích mí Laplacián následující tvar: aó = - d (r2 dó) + 1 d (sin ůdó) + 1 cľ± (33) r2 dr \ dr J r2 sin ů dů \ dů J r2 sin2 ůdp2 Separace proménných: Pokousíme se přeménit parciální diferencialní rovnici do tří obýcejnách diferencialních rovnic pro funkce jednotlivách promenných, R(r), 0(9) a $(p). K tomu predpokladame resení ve tvaru soucinu Ó(r, ů, p) = R(r) • O(ů) • $(p) (34) a dosadíme do Laplaceový rovnice 1 d ( 2 dR ^ \ 1 ô ( n „de \ 1 „ ^ d2$ AÓ = r2 — • O • $+-—— — sinů • R • — • $ + 2 . O — = 0. r2 y dr y r2 sin ů dů \ dů J r2 sin2 ů dp2 2.2..... (35) Nísobíme rRQ1^ 9 a píseme cast, zívisející na p, na pravou stranu sin2ů_d (2 dRA sinů á_ (. ? d0) = 1 d2$ Leví strana ted' zavisí na r a ů, prava strana na p. Z toho výplíva, Ze se obe straný musí rovnat konstante, kterou nazveme m2. Z prave straný dostaneme obýcejnou difer-enciíalni rovnici, d2 $ — + m2$ = 0. (36) dp2 Rovnici výplývající z leve straný muZeme upravit podobním zpusobem, IA (r2 dR A = Jll___LI A (sin ů d0 A = A2, R dr V dr y sin2 ů sin ů O dů y dů y ' kde se opet obe straný musí rovnat konstante, oznacene A2. Z toho dostaneme dve dalsí obýcejne diferencialní rovnice, dr(r2 f) - a2r=o (37) " sinů dů (sin ů S) + (A2 - sŘ) e = 0' (38) tak ze místo parciílní diferencialní rovnice mame rovnice (36), (37) a (38). Resení: Rovnice (36) ma resení $m(p) = Cm cos mp + Sm sin mp (39) 8 príslůsne k párámetrů m, ktery můsí byt celocíselní, áby $ bylo periodicke ve p. Rádiílní rovnice (37) má resení Ri(r) = Ai r1 + ^, (40) kde A2 = l(l + 1). V rovnici (38) píseme cosů = x. Resení, ktere obsáhůje obe integrácní konstánty m á l, ožnáďme Pm. Tákovoů ůprávoů dostáneme Legendreovů rovnici (i — x2)d!pjmrx) — 2x ^ + ['(' + i) — i^] Pľ'(x) = 0. (41) 5.2 Legendreovy polynomy Ortogonální báží resení pro m = 0 jsoů Legendreovy polynomy Pi(x), vyhovůjící jednodůssí rovnici dd- [(1 — x2)dPj^]+ l(l + 1) Pi(x) = 0 (42) s nežáporním celocíselnym párámetrem l. Legendreovy polynomy jsoů ortogonální v interválů (—1,1) J Pk(x) Pi(x)dx = 0 V k = l. (43) Legendreovy polynomy se objevůjí jáko koeficienty v rožvoji tžv. vytvárející fůnkce (1 — 2xt + t2)./2 = g P(x) *■ (44) Poůžitím Leibnižová právidlá dm[/(x) g(x)] = m m! dm-k / (x) dkg(x) (45) dxm k=0 k!(m — k)! dxm-k dxk ( ) dostáneme m-níásobníym derivováíníím rovnice (42) (1 — x2)f" (x) — 2x(m + 1)/'(x) + (l — m)(l + m + 1)/(x) = 0, (46) kde /(x) = dmPl(x)/dxm. Sůbstitůce /(x) = (1 — x2)-m/2g(x) vede k tomů, že fůnkce g(x) můsí splnovát rovnici (41), je tedy konecne Piľ(x) = (1 — x2)m/2 dmPlM. (47) Legendreovy polynomy lže vyjádrít pomocí Rodrigůesová vžorce Pi (x) = ^dx(x2—1)1■ (48) Vyůžitím tothoto vžtáhů můžeme rožsírit (47) ná oblást žáporních m, tedy Pim(x) = (1 — x2)m/2^(1 — x2)i, —l < m < l. (49) Polynomy Pm se nážívájí pridrůžene Legendreovy polynomy. Námi definováne Pm(x) nebo Pi(x) nejsoů ná interválů (-1,1) normováne ná jednicků. Ostátne růžne drobne i vetsí odchylky v definicích speciálních fůnkcí jsoů díky historickemů vívoji bohůžel žcelá bežne. 9 5.3 Kulove funkce Pomocí přidružených Legendreových polynomů definujeme úplný ortonormální soubor kulových funkcí (t.j. káždou funkci úhlových promenných ve sferických sourádnicích můžeme nápsát pomocí (nekonecne) rády techto funkcí) 77X^7 P™(cos ů) exp(rn^). (50) Plátí tedy Jo dV Jo dů sinůY^^ů,
R (55)
10
(57)
Potencial mimo koule je dan vžorcem (25)
*(x) = ~A- / d3X' i,. =-- / d3X' . K ' -
AneoJ \x — x \ 4n^oj 2 + x2 — 2 \X\ • \x' \ cos 7
= /d3x' ; P("') ■ (56)
4n£0 \X 1 ^ — 2 X cos 7 + (f)2
7 je uhel meži X a XV integrílu se objeví vytvarející funkce Legendreovích polynomu, tak
11 r °° rr'\1
*(X) = — -\ d3x' p(X') £ Pi (cos 7) 7J
(psali jsme \X\ = r a \X'\ = r'). Použítím (54) dostaneme rožvoj
1 r 4n r'1
*(X) = —j d3x'p(XOEgT+TrjV1 Y^(V',v')Yr(V,v)■ (58)
Pomocí multipílovích momentu
qm := / d3x' p(X') r'1 Y™*(ůV) (59) mužeme konecne psat potencial jako superpožici kulovích funkcí
*(x) = E 2ľ-T^ Yr(^. (60)
V prípade bodoveho naboje víme, že pole je dano Coulombovím potencialem. Je-li naboj q umísten mimo pocatek souradnic, napr. na ose z (v bode z = R), je potencial dan vžtahem
*=4^žfl(c°s'')(R)'' r sR
'=0 1 (61)
q 0 R
Pro r R prevažuje rotacne soumerna (vžhledem k pocatku souradnic, nikoli polože naboje) složka l = 0. Umístíme-li vsak na ose z jeste naboj opacne velikosti do z = — R, vyrusí se identicke príspevky clenu s l = 0 a pro r >> R prevažuje pak dipóloví složka
(l=1)
2qRP1(cos ů) D cos ů *dip = -.--2- = -.--ž~ , (62)
kde D = 2qR ožnacuje dipílovy moment. Podobne, umístíme-li na ose z v z = ±R naboje q a v pocatku naboj — 2q, vyrusí se identicke príspevky clenu s l = 0 a l = 1 a pro r R prevažuje pak kvadrupolova složka (l = 2)
. _ 2qR2 P2(cosů) Q 1 — 3 cos2 ů
*quad = — --3- = "j--3-, (63)
kde Q = qR2 je kvadrupolovy moment. Obecne jsou multipolove momenty žívisle na umístení v souradnem systemu, s vyjimkou nejnížsího nenuloveho momentu.
11
6 Magnetostatika
6.1 Analogie mezi elektrostatikou a magnetostatikou
Integrální tvar infinitesimální rovnice (4), vyjádřený pomocí hustoty proudu, je B (x) = ^ í dV f(x') x ,! ~ t' 3 = - ^ í d3x' f(x') x V—
= ^V x í d3x' f(xl, =: rot f(x). (64) 4n J \x - x \
Zavedli jsme vektorový potenciál
A(x) = ^ I dV-P^. (65)
Vektorový potenciál není jednoznaCný, protože mUzeme priCíst gradient libovolne funkce, jehoz rotace je identicky nulová. Taková transformace,
A(x) —> A(x) + grad A(x)
se nazívá kalibracní transformace. (Stejne je skalární potenciál jednoznacní jenom až na konstantu.) Volba A = 0 (Coulombova kalibrace) vede k vlastnosti div A = 0, protoze V\x - x= -V '\x - xa z toho dostaneme integrací per partes V 'f(x'), coz se rovná nule podle staticke rovnice kontinuity. Dosazením do Maxwellovy rovnice dostaneme tak
rot B (x) = rotrotA (x) = graddivA (x) - AA (x) = -AA (x)
= -£° / d3x' j(x') A—^ = m (x), (66)
tedy vektorovou Poissonovu rovnici pro A.
Následující tabulka ukazuje analogie mezi elektrostatikou a magnetostatikou.
Elektrostatika Velicina, vztah Magnetostatika
dF = dq E definice pole dF = I df x B
dq hustota zrídel I dl = f d3x
F = qE síla na níaboj F = qv x B
div E = p/t0 rot E = 0 rovnice pole rot B = p,0f div B = 0
E = -grad 0 potenciíaly Bf = rot Af
= - p to rovnice potenciaílu AA = -^f
0(x) = ~ě^l d3x' p(x') potenciál zrídel Af(x) = |°/ d3x' f (x')
12
Analogicke magneticke pole k elektrickemu Coulombovu poli je pole linearního vodice. Uvažujme nekonecne dlouhý, rovní drat ve smeru osý z s elektrickím proudem
I
j (X) = Iô(x) ó(y) e3. (67)
Podle (64) je magneticke pole
4n J-oo \x — (0, 0, z')\3
Zavedeme polírní souradnice p = \Jx2 + y2, p = arctan ^, z a uvedomíme se, ze pole musíí býít konstantníí ve sm eru z, tak ze sta cíí po cíítat BV v rovin e (x, y). Víýpo cet vede k Biotovu-Savartovu zakonu
B (P) = ^. (69)
6.2 Magnetické pole kruhové smyCky.
Do vztahu pro vektoroví potencial (65) dosadíme hustotu proudu J(x') d3x' = 15(p' — a) ô(z') ev> p' dp' dz' dtp', kde ěy = — sin(p' — p) ep + cos(p' — p)ev, a dostaneme
.4 (p,z) = ^ (p'z)ěV' Av(p,z) = — —-2-2—ň-^. (70)
2n jo (a2 + p2 + z2 — 2ap cos p) 2
Integral lze výjadrit uplními eliptickými integralý E a K,
A,(p, z) = Q* [ (l — K (k) — E (k)] , (71)
Aip(p,z) =
kde
k2 = ( +7 + 2 , K (k) = í * d4 2 , E (k) = / V1 — k2 sin2 e de.
(a+p)2 +z2 j° ^1 — k2 sin2 e Jo
(72)
Pri výpoctu indukce potrebujeme identitý
dE (k) = E (k) — K (k) dK (k) = E (k) K (k)
dk = k , dk = k(1 — k2) k~. ( )
Potom mame pro slozký indukce (B^ = 0)
BP(p,z) = —^ = ^ , z —K(k) + a2 + p22 + z22 E(k)l , (74)
^ ' ôz 2vr (a + p)2 + z2 L (a — p)2 + z2 ^J' 1 '
Bz (p, z) =--tt^- = t;--/ = K (k) +-- 2 E(k) . (75)
p ôp 2n ^(a + p)2 + z2 L (a — p)2 + z2
13
7 Maxwellovy rovnice v materiálovém prostředí
7.1 Polarizace, magnetizace
Když se materiálové prostředí nachází pod vlivem vnějšího elektromagnetického pole, např. mezi deskami kondensátoru nebo v magnetickem poli cívky, rozlišujeme vnejší náboje a proudy, pext a pxt, a naboje p a proudy jindukovane vnejsím polem v materiálu. Naboje a proudy uvnitr atomu, molekul, nebo elementarních buňek krystalu, ktere vyvolavají mnohem vetsí pole než jsou makroskopicka, mužeme zanedbat, když stňredujeme pňres prostorovíe oblasti podstatnňe vňetňsí neňz vzdalenost elementíarních jednotek materialu. Stredovane Maxwellovy rovnice jsou
divE = {P) + "ext rotE = - d4-
C\ (76)
1 d E
—rot B = eg — + {j) + Jext divB = 0.
p,0 dt
E a B jsou stredovaní pole, {) oznacuje stredovaní zdroju.
Celkový níboj vazaní na prostredí, ktere je plne uzavreno uvnitr oblasti V, je roven nule
0 = {p) = -divp, (77)
/ {p) dV
pricemz P = 0 vne materialu. Potom je totiz
/ {p) dV = - í divP dV = í P • n dS = 0, Jy Jy Js
kde S = dV. Uvazujme dipolovy moment
J x{p) dV = - J x divP dV = - J x {n • P) dS + J (P -V)x dV = J P dV. (78)
Uvaňzujme nyníí uzavňrenou plochu S uvnitňr materiíalu. Celkovyí proud touto plochou vazaní na prostredí je dín celkovou hodnotou casove zmeny prumetu vektoru polarizace P
f {])■ n ds = -f d^ dV = / d divP dV = f dPJ • n dS
J s J v dt J v dt J s dt
dPP
= =rot M + —, (79)
pricemz M = 0 vne materialu.
/rot M • n dS = 0 protoze dS = ®. s
V statickem prípade je dP/dt = 0. Uvazujme magneticky moment
1 í x x {j) dV = 2 ľ x x rot M dV = (80)
2 J V 2 J V
2 ľ x x (n x M) dS - 1 í (M x V) x x dV = í M dV.
2 J S 2 J V J V
14
Definice vektoru polarižace P a magnetižace M pomoc! momentu je duležita pro jed-nožnacnost, jinak by vyhovovaly take P + rot/ a M + grad f. Jako vedlejsý výsledek dostaneme že spojený rovnic (77) a (79) diferenciýalný (= mikroskopickou) rovnici kontinuity
^ + div(J) = 0. (81)
7.2 Makroskopické Maxwellovy rovnice
Zavedeme vektory indukce elektrickeho pole a intenžity magnetickeho pole jako
D = t0Ě + Ě, H =— B - M. (82)
Použitým polarižace (77) a magnetižace (79) dostaneme makroskopicke Maxwellovy rovnice,
, á á dB
divD = pext, rotĚ = - — ,
dt (83)
á dD _ á rot H = — + /ext, div B = 0.
dt
Jak Maxwellovy rovnice ve vseobecne forme, tyto rovnice jsou konsistentný s rovnicľ kontinuity
^ + div áext = 0. (84)
V homogenným isotropným linearnmi prostredý bež disperse jsou vžtahy meži D a Ě a meži H a B obžvlast jednoduche, pole jsou ýmerna,
D = tr e0Ě, H =^— B, (85)
kde tr je dielektrický konstanta a konstanta fir je magnetický permeabilita materiýlu. („Lineýrm prostredý" žnamený, že tyto veliciny jsou konstantný.)
7.3 Energie a impuls elektromagnetického pole
Mejme spojite rožložení naboje p. e ožnacuje energii nýboje obsaženeho v objemu A V, Ae žmena energie pri pohybu v elektromagnetickem poli. p a / ožnacují ted' externí veliciny.
- - 1 Ae
Ae = F ■ Ax, F = pĚ AV + j x B AV = — — = j ■ Ě. (86)
AV At
S využitým vžtahu
ĚĚ ■ (ý x H) - H -ý x Ě) = V ■ [H x Ě) (87) odvodýme ž Maxwellovyých rovnic vyýraž
H ■ ddB+ě ■ddD=-j-ě - V ■ (ě x h) . (88)
dt dt
15
Na pravé straně vystupuje vykonaná práce a divergence toku, výraz na levé straně
můžeme tedy interpretovat jako Časovou zmenu hustoty energie. Po zavedení veliCin hustoty energie W a Poyntingova vektoru S
W = 1 (Ě - D + B - H), S = E x H (89)
muzeme (88) psat jako
díw dV + / j-E dV + / S - n dS = 0. (90)
dt j v j v jav K '
S má tedy význam hustoty toku energie. Obdobnou uvahu mužeme provest pro impuls. Pri prechodu ke spojitemu rozložení naboje je
Ap = FAt, F = pE AV + j x B AV = = pE + j x B. (91)
Z Maxwellových rovnic odvodíme výraz Bfí f)f)
d x—+-Qt~ x B = E (ý - ď) - b x (ý x Hr) +H (ý - b) - d x (ý x E)-jx b - pE.
\ t . o (92)
Poslední dva cleny na prave strane popisují Lorentzovu sílu, muzeme tedy casovou
derivaci na leve strane interpretovat jako casovou zmenu hustoty impulsu
G = D x B = trnreo^o E x H = ^ S, (93)
c2
pokud první c ty ri c leny na prave strane lze psát jako divergence toku impulsu. Zavedli jsme
c2 :=-, n2 := er fir. (94)
c je rychlost sírení elektromagnetických vln ve vakuu, tedy rychlost svetla, jak budeme brzo videt, a n je index lomu, jehoz vyznam bude take z rejmý v souvislosti s s írení elektromagnetickych vln. Po uprave, kdy p redpokladýme, z e permitivita a permeabilita nezavisí na prostorovych souradnicích, muzeme psat
'e(v - d)- d x (vx e)] .=j:—d3 - 2 e - ď
j—1 3
3 _0 „ 1 )
3 ) 1
(95)
'H(ý - B) - B x (V x H)}i = ^ — (Ht Bj - 2 bij H - É)
j—1 j
a zakon zachovaní ma tvar
d! Gi dV + f \pEt + (j x É)] dV W ]T Ttj nj dS = 0. (96) )t v v i av j—1
16
Definovali jsme Maxwellův tensor napětí Tý jako
Ta = -(Ei D, + Hi B,) + - ótJ(E • D + H • B).
2
(97)
Takto definovaný Maxwellův tensor určuje tok impulsu z uvažovaného objemu. Jeho stopa je rovna hustote energie
w -E Ta = 0.
(98)
i=1
8 Integrální formy Maxwellových rovnic, indukce
8.1 Integrovane rovnice
Uvažujeme prostorovou oblast V s okrajem dV a pločhu S s okrajem dS. n označuje normainí jednotkový vektor na pločhu dV, popr. S, t tečný vektor na dS. Integru-
jeme rovnici diví? elektrostatiký.
pres objem V, dostaneme podle Gaufiový vetý Gaufiuv zakon
Jav
E-n dS
Q
(99)
Když integrujeme rotaci elektrického pole přes plochu S, můžeme uplatňovat Stokes-ovu vetu (t je teCny vektor okraje dS):
é E] •t ds = rot E] •n dS = Jas Js
d
_ B ■ n dS. ôtJs
(100)
Integra^ B-n dS definuje magnetičký tok $TO pločhou S. Odvodili jsme tedý Faradaýuv indukční zakon: Električke nýpetí v uzavrene dratene smýčče se rovna minus časove derivači magnetičkeho toku,
- - d E] •t ds = -— $m.
as dt
(101)
Analogičký dostaneme v statičkem prípade E = 0
é B - t ds = rot B - n dS = [io j - n dS. Jas Js Js
(102)
Integral na prave strane predstavuje čelkový električky proud I, z toho výplýva Ampereuv zýkon
as
as
IB •t ds = hq I.
(103)
17
8.2 Aplikace — vzajemna indukce a vlastní indukce
Uvažujeme dve geometricky pevne smycky s promenným proudem ve smycce 2. In-dukovane napetý ve smycce 1 vyvolane žmenou pole buženeho smyckou 2 je
U1 = - f B2 ■ ň1 dS1, / B2 ■ ň1 dS1 =1 Ä2 ■ dľi, Ä2 = —. (104)
at J (i) J (i) J(i) 4n J (2) r12
Po dosazení dostáváme
Ui = Mi2 ^, Mi2 = - f i i . (105)
dt 4n i(i^(2) \xi — x2\
Pokud by tekl proměnná proud smyčkou 1, bylo by indukované nápetí ve smyčce 2
d J
U2 = M2idjji, Mi2 = M2i = M. (106) dt
Ale take zmena magnetickeho toku smyčkou 1 vytvorí indukovane nápetí v teto smyčce, stejne platí pro smyčku 2. Obecne tedy muzeme psat
Ui =—Li dt+M § U=M §—Lf §• (107)
Casova zmena energie magnetičkeho pole je rovna zaporne vzate príči
f = 'i — U2J2 = L ^ + — M(Jif + J<£) . (108)
takze pro energii magnetičkeho pole je
W =2 Li I2 + 1 L2 I| — M li I2, Li L2 > M2. (109) Energii magnetičkeho pole mame ovsem take vyjídrenou jako
W = U B ■H dV 4 U-ÄdV. (110)
Pri odvození obou vírazu v teto rovniči je postupne využito vztahu
B = rotÄ, H ■ rot Ä — Ä ■ rot H = div [Ä x Hřj , rot H = Ä- (111)
Vztahu pro energii vyuzijeme pro vípočet vlastní indukčnosti
1
L
B2 dV. (112)
Jy
fi012 J V
Uvažujme dve solenoidalní cívky, každou o N žavitech, tesne na sobe. Prurež cívek je S a jejich delka í. Pole prvný a druhe cŕvky jsou tedy približne (Ampereuv žakon)
Bi « tí^EL, B2 « ^ (113)
a pro indukcnosti mame
Li « L2 ~ M « 0(°N£ S. (114) Pro energii magnetickýeho pole pak
W = ^ (Ii + h?. (115)
18
9 Casove promenna elektromagneticka pole 9.1 Dynamicke potencialy, kalibrace
Predpokladame dýnamicke potenciálý $(x, t) a A?(x,t) a B = rotA. Pak výplává z
Maxwellový rovnice rotE = —B, ze i casova derivace vektoroveho potencialu prispíva k elektrickáemu poli
Dosazenáí do dal sáích Maxwellovýách rovnic vede k
d A
B = rotA, E = -gradó - —■ (116)
Aó + — div A = - — dt en
AA - t0ß0dp - grad I divA + e0ß0) = -ßoj-
(117)
S využitím kalibrační transformace
Ä — Ä +gradt, 0 — 0 - (US)
můžeme mít Lorenžovu kalibraci (Ludwig Valentin Lorenz = Hendrik Antoon Lorentž)
divÄ + t0ßo 90 = 0 (119) dt
a dostáváme tak pro potencialy nehomogenní rovnice
d2 Ä I
(120)
Oznacili jsme rýchlost svetla ve vakuu c
= 1
Rovnice (120) jsou Maxwellový rovnice pro potencialý, spolu s kalibraci' (119) jsou ekvivalentná (5).
9.2 Rovinna a kulova vlna
V prápade volneho elektromagnetickeho pole popisujá homogenná rovnice odpovídajíCí (120) sŕrem vln. Vlnova rovnice v jednorozmernem prápade popisuje rovinnou vlnu (ve smeru x)
d 2^(x,t) 1 d2 #r, t) =0 (121)
dx2 c2 dt2 Obecné rešení je
#M) = f(t - X) + g(t + (122)
19
Vezmeme jako príklad Gaufiovu funkci f = exp[-(t - x)2]. Maximum se nachazí pri t - x = 0, pohybuje se tedy rychlostí c ve smeru rostoucího x.
Dalsím jednorozmernem príkladem je sfericky symetrickí vlnoví rovnice v troj-rozmňerníem prípadňe,
1 d 2(r//(r,t)) 1 d2^(r,t) =0
r dr2 c2 dt2
s obecnym resením
/(r,t) = lf(t - ri) +1 g(t + ri). (124)
Na toto resení se muzeme divat jako na rozbíhavou nebo sbíhavou kulovou vlnu. Tvar ňreňsení takíe ukazuje, ňze rychlost ňsíňrení je c.
V (lineírním) materiílovem prostredí se nahrazuje e0 — ere0 a fi0 — (irfi0. Pak dostaneme z rovnice (120) rychlost sírení c/n, kdyz zavedeme index lomu
9.3 Obecne řešení nehomogenní rovnice pro potencialy.
Pro obecníe ňreňseníí rovnice (120) jeňstňe chybíí partikulíarníí ňreňseníí nehomogenníe rovnice. Zavedeme Fourierovu transformaci casove zavislosti potencialu a hustoty níboje
1 r°°
Vť =--7-^
(241)
Vyraz pro potenciály ve (239) pak budeme chípat jako funkce f (x,t'), a budeme pocítat parciální derivace podle x pri konstantním t' a podle t' pri konstantním x. Porovnáváním diferenciálu
Of dt
Of Oť
df = Vf • dx + Oť dt = Vf • dx + OL dt = V + ^ Vť) • dx + Oť O- dt (242)
prepíšeme (239) jako
Oj_Oť_
E = - V jj(x,ť)
O(j)(x, ť) f , OA (x, ť) Oť
t
Oť Ot
B = V x A(x,ť) + Vť x OA(fx,t')
Pro intenzitu elektrickíeho pole dostaívíame pak
Ef
4neo
1 _ v_
1_c_
i)
R2 1
c
+
n x
n
"t
x w
c2 R (1 - v-f
zatímco pro indukci magnetickeho pole pole
B = - n x E =-
c 4n
1 _ v_
1 2
c2
(v x n)
R2 1
1
n x
n x ((n -cR Í1 -
x W
3
(243)
(244)
(245)
Oznacili jsme jednotkoví vektor n = R/R a zrychlení w = dv/dt'. Limitní prípady pro v/c — 0 jsou
Ef
efn
4vreo R2
Bf
e^0(v x n) 4nR2
(246)
12.2 Intenzita zaření
Poyntinguv vektor (energie, procházející jednotkovou plochou za jednotku casu, dimenze [Jm-2s-1]) je
f = — E x B = tocE2n (247) Vo
1
1
e
c
c
34
a intenžitu žáŤem (tj. energie, vyžarovanou ža sekundu do elementu prostoroveho uhlu, [Watt]) spocteme tedy jako
dl = lim Š • n R2 díl (248)
R—>oo
Po dosaženíí ž (247) a (244)
2(ř? • w)(v^w) w2 (1 — (n • w)2
e2
dl
16n2e0c3
C2
dí (249)
Pro v/c — 0 dostavíme s ožnacením n • w = w cos £ pro celkovou vyžarovanou intenžitu
e2w2 /■ 2n ľ n e2 w2 I =T6Tw-J dW d£ sin £(1 — cos2 £) = 6-- ■ (250)
V klidove soustave castice je tedy (s ožnacením I = dE/dt)
, ^ e2w2 , * dx* , ->. * du* í w\ ,
dE =6^60^ dť d'P = 0, u* = — = (1,0), w* = — = [0,^j ■ (251)
Relativisticky invariantní vyraž (tj. diferencial ctyrvektoru impulžu) vytvoreny ž ctyr-vektoru rychlosti a žrychlení, kterí v klidove soustave prejde na víražy že vžtahu (251), je pak (m = 1)
* E * e2 duk duk * e2 duk duk *
p* = -,p), dp* = —--— —k dx* = —-— —k u* ds. (252)
V c J 6ne0c ds ds 6ne0c ds ds
V laboratorníí soustav e tedy míame pro celkovou vyža rovanou intenžitu vyíraž
e2 w2 — (w x f) 2
6ne0c3 ^1 — v^A
(253)
Zde jsme pot rebovali vyjíad reníí cty rvektoru rychlosti i žrychleníí v laboratorníí soustav e. Abychom nemuseli pri vípoctu ctyrvektoru žrychlení užit obecne Lorentžovy transformace, vypocteme w* derivovíním žnameho tvaru u* = 1 — v^,v/ {c\J 1 — ^
potom (
(
w*
v ^ w__w + v(v ^ w) i (254)
cS i1 — S) c211 — ^) c4 í1 — $
\C ^ — J) ° V " ^ ^ — C2
V homogenním magnetickem poli se nabita castice pohybuje rychlostí v po kružnici polomeru R, její žrychlení w = v2/R je kolme k rychlosti. Dosažením do vžtahu (253)
I = v4 = ( JL Y~ (JLf 4 (255)
6ne0c3 R2 fi — v£j2 6ne0R2 \mcj ~ 6ne0R2 Vmc2/ ' V ;
V posledním víražu ve (255) jsme použili aproximace vysokích energií, kde pro kinetickou energii platí T = \Jp2c2 + m2c4 — mc2 ~ pc. Z tohoto vyražu je take žrejme, že
I
35
synchrotronové záření je omezujícím faktorem při urychlování lehkých částic (elektronů a positronů). Pro normovací hodnotu R0 ~ 0, 5 km můZeme psát I « (Ro/R)2(T/mc2)4 eVs-1.
Jsou-li rychlost a zrychlení v urcitem okamZiku rovnobeZne, dostáváme (n -v = v cos ů, rychlost podel osy z) pro uhlove rozloZená zárem vyraz
dI
e2w2 16n2e0c3
sin2 ů
(l - v cos ů) 6
díl
(256)
Pro hodnoty v/c — 1 ma uhlove rozlozená velmi uzke, ale „dvouhrbe" maximum kolem ů = 0. Jsou-li rychlost a zrychlená v urcitem okamziku navzajem kolme, dostaváme (n - v = v cos ů, n - w = w cos ip sin ů, rychlost podel osy z a zrychlená podel osy x) pro áhlove rozlozená
dI
e2w2 16n2eoC3
1
1 - v cos ů
c
(1 - g) sin2 ů cos2 p - f cos ů)6
(257)
13 Maxwellovy rovnice v čtyřrozměrné formulaci
13.1 Čtyrrozměrný vektor proudu, rovnice kontinuity
Definujeme ctyrvektor proudu (pro castici: naboj krát ctyrrychlost)
j1 = P
(cp,pv) = (cp,J).
Naboj, ktery ubude v nejakem objemu, muzeme zapsat dvojím zpusobem
_3_ - dt
J p dV = j> j - n dS.
S pomocí Gaufiovy vety pak z (259) plyne
dV = 0,
tedy (objem je libovolnyá) rovnice kontinuity
(258)
(259)
(260)
^ _ dp dj
V - j + — = — = 0.
dt 3x%
(261)
13.2 Homogenní Maxwellovy rovnice
Z vyjádření tensoru elektromagnetického pole pomocí potenciálu snadno odvodíme platnost vztahu
dFtk + dFja + dFu dxl dx% dxk
0.
(262)
36
Na levé straně je úplně antisymetrický tensor třetího řádu, představuje pouze čtyři různé rovnice. Zřetelneji je to videt, uZijeme-li zápis duálního (pseudo)vektoru
1 eiklm dFlm = d*F ik = o
2 dxk dxk
(263)
Nultá komponenta davá tvrzení o nezrídlověm charakteru magnetickěho pole, dalsí tri komponenty Faradayuv indukcná zakon
V • B = 0, V x Ě
dB "dt'
(264)
13.3 Nehomogenní Maxwellovy rovnice
Ctyřrozmerný zápis nehomogenních Maxwellových rovnic, obsahujících hustotu náboje a proudu, je
dF
ik
dxk
-ßoj i.
(265)
Nultá komponenta je rovnice pro divergenci intenzity elektrickeho pole (Gaufiova věta elektrostatiky), zbávající tri pro rotaci magnetickěho pole (Ampěreuv zákon)
V ^Ě = í V x B = +«■
(266)
13.4 Tensor energie-impulsu
Z hustoty energie
W = \ (t0Ě2 + — £2)
2
z Poyntingova vektoru
l
a z Maxwellova tensoru nápětí
S = — Ě x B ßo
l
&aß = Í0ĚaĚa +--BaBß - W&aß
ßo
(267)
(268)
(269)
elektromagnetickeho pole (ve vakuu) muzeme sestavit ctyrrozměrná tensor energie-impulsu,
iS " ) . (270)
Tik =
Pomocí tensoru elektromagnetického pole dostáváme jednoduchý výraz
T
ik
l
ßo
(
-gimFilFkm +
4 gik FimF lm^j
(27l)
37
Tensor energie-impulsu soustavy částic zapíšeme pomocí analogie s tensorem energie-impulsu elektromagnetického pole. Hustotu impulsu soustavy částic napíseme jako
na = /icua, / = J2 maS{3) (x — xa). (272)
a
Hustota impulsu je u elektromagnetickeho pole rovna hustote toku energie delene c2. Vyraz (272) bude tedy analogicky roven T0a/c. Velicina /ic je nultou komponentou (stejne jako hustota naboje u ctyrvektoru proudu) ctyrvektoru toku hmoty i dxl/dt. Tensor energie-impulsu tak muzeme konecne psít jako
Tik = /c ^= icuiuk dS, Tik = Tki. (273) ds dt dt
Pro tensor energie-impulsu elektromagnetickeho pole dostaneme s vyuzitím Maxwell-ovíych rovnic
= — lo J \ eiklm ^ = 0 (274) d
víyraz
qI^f = —Fkjk. (275) Pro tensor energie-impulsu soustavy cíastic dostaneme s vyu zitíím pohybovyích rovnic
ic ddUÍ = p Fik uk <* ic dutÍ = Fik jk (276)
ds dt
a rovnice zachovíaníí hmotnosti (rovnice kontinuity pro cty rvektor toku hmotnosti, podobn e jako pro cty rvektor proudu)
dxk \ dt
víyraz
d(l dí) = 0 (277)
d Tk = Fik jk. (278)
Spojeníím (275) a (278) dostíavaíme zíakon zachovaíníí
^ (t* + tpk) = 0. (279)
Pro tensor energie-impulsu platí (rovnost prave pro elektromagneticke pole)
ti > 0. (280)
38
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ů
dF'k -o, r» = ^>'aAl aAj
axk ' v axj axl)'
%1 a2 Ak ki a2 a
9 —:--g -= 0.
Lorenzova kalibrační podmínka (119) nabývá formu čtyřdivergence
aAk
(281)
axk
a zjednoduší (281) na vlnovu rovnici
0 (282)
ki a 2 Ai
9 axk axl ( )
Pomocí d'Alembertova operátoru
máme pak ve trírozmernem zapisu
+ V • AA =0, U