Schmidtova hlava s vyznačenou vlnou Poznámky V[O]AF — elektromagnetické vlnění

Známe z elmy: E(r,t),B(r,t),ε=εrε0,μ=μrμ0

Uvažujme prostředí bez volných nábojů a proudů. Známe Maxwellovy rovnice:

E=0(Gaussu˚v zaˊkon) B=0(neexistence magnetickeˊho monopoˊlu) ×E=Bt(Faradayu˚v zaˊkon) ×B=εμEt(Ampeˊr-Maxwellu˚v zaˊkon)

Z toho dostaneme jakési homogenní řešení. Pokud přidáme zpátky náboje a proudy, stačí přičíst partikulární řešení, ale to ať si najde někdo jiný.

Zapůsobme rotací na Faradaye: Rotující Faraday
×(×E)=t(×B)=εμ2Et2

Zároveň ×(×E)=alias baˊc mıˊnus caˊbrotace rotace jegrad div mıˊnus laplace(E)2E=Gauss2E, z toho plyne

2Et2=1εμ2E

Analogicky lze odvodit

2Bt2=1εμ2B

To jsou vlastně vlnové rovnice. Označme si rychlost:

v1εμ c1ε0μ0(ve vakuu) ncv=εμε0μ0(index lomu)

Pokud nechceme být úplně obecní, můžeme hledat řešení rovnic v d'Alembertovském tvaru:

E=E0F(krvt) B=B0F(krvt)

To není úplně obecné, protože předpokládáme stejnou funkci a stejný směr šíření.

To jsou jistě řešení vlnových rovnic, ale ještě jsme nedokázali, že jsou to i řešení původních Maxwellových rovnic. Dosaďme tedy (s hnusným fyzikálním zápisem):

E=iEi=i(E0iF(njxjvt))=E0iniF(nrvt)=maˊ byˊt0 E0n=0E0r

Analogicky

B0n=0B0r

Tedy E,B jsou příčné vlny. (Za tohle hrozí okamžité vyhození u zkoušky, pokud to nevíte!)

Dosaďme dále:

×E=Bt εijkjEk=tBi εijkj(E0kF(nlxlvt))=t(B0iF(nlxlvt)) εijkE0knjF(nrvt)=B0i(v)F(nrvt) n×E0=vB0

Z toho plyne:

Příklad (harmonická elektromagnetická vlna). F(x)exp(ikx) E=E0exp(i(ωtkr)) B=B0exp(i(ωtkr)) ω=v|k|=v|kr|

Pokud je vlna jen jedna, potom bez újmy na obecnosti můžeme volit r=rz, potom vlna vypadá jako stará dobrá vlna na struně:

E=E0exp(i(ωtkz)) B=B0exp(i(ωtkz))

Vyzařování elektromagnetických vln

Mějme v prostoru jeden náboj q v počátku.

E se pohybuje společně s nábojem, tedy v daném čase musíme počítat s posunutým počátkem:

E(r,t)=qr4πε0r3,rrvt

Uvažujme nyní, že v čase t0 je náboj v klidu, v čase 0tτ rovnoměrně zrychluje se zrychlením a a v čase tτ se pohybuje rovnoměrně rychlostí v=aτ.

Jsme ve vakuu, takže elektromagnetické vzruchy se šíří rychlostí c. Tudíž pokud mě zajímá, jaké je u mě elektrické pole, musím vědět, kde byl zdroj v retardovaném čase tr=t|r|c.

Tím vzniknou jakési dvě slupky oddělující tři prostory, ta větší kolem počátku a ta menší kolem bodu, kde náboj přestal zrychlovat. Uvnitř menší slupky a vně větší slupky jsou radiální čáry od jejich středu, mezi nimi propojení konců čar, čímž vzniká jakási divná pavučina.

Předpokládejme τt, tedy „akcelerační slupka“ je daleko od počátku. Vezměme něco jako válec (vlastně spíš komolý kužel), který má podstavy kolmé na radiální směr a jeho pláštěm prochází vnější slupka, a aplikujme na něj Gaussův zákon:

SEdS=Quvnitrˇε0=naˊboj nenıˊve vaˊlci0

Plášť do integrálu nepřispívá, protože normála na něj je kolmá na pole, takže

0=SEdS=podstavyEdS

Pole v bližší podstavě je rovnoběžná složka toho v akcelerační slupce, pole ve vzdálenější slupce je to statické. Tedy

0=EstatSpESp E=Estat

Pomocí nějaké pofidérní geometrie dostaneme i kolmou složku:

cτvt=podobnost EE E=vtcτE=arc2E E=qa4πε0c2r

Z toho plyne:

Eradiacˇnıˊ(r,t)=qa(tr)4πε0c2r,tr=trc

Teď spadne z nebe: BE,B=Ec,Bn,E,B,n pravotocˇiveˊ (prostě Létající špagetové monstrum určilo, že to funguje stejně jako předtím a nepochybuje se o tom).

Příklad (oblíbená harmonická vlna). a(t)a0cos(ωt)z a=asin(θ) a=aa=a(ar0)r0,r0r|r| a(t)=a0cos(ωtkr)z,ω=ck Eradiacˇnıˊ(r,t)=q4πε0c201rsin(θ)cos(ωtkr)

Energetické veličiny v elektromagnetickém poli

Definice. Hustota energie wdEdV
Věta (viz ELMA). w=12(εE2+B2μ) (platí pro každé pole, nejen pro nějaké, jak bylo odvozeno na elmě)

Tok energie je

wt=12t(εEjEj+BjBjμ)=εEEt+BμBt=MaxwellovyrovniceE(×B)B(×E)μ=E×BμPoyntingu˚vvektor

Poslední úprava je v podstatě jiná varianta „bác mínus cáb“.

Tedy energie z daného místa „vytéká“ ve směru, do kterého ukazuje Poyntingův vektor S.

Nyní aplikujeme Gaussovu větu:

ddtVwdV=AVSdA dEVdt=ΦS
Definice. Hustota hybnosti gdpdV
Věta. Ve vakuu platí: g=Sc2=ε0E×B
Příklad. Vezměme rovinnou elektromagnetickou vlnu a strčme jí do cesty nabitou částici q. Bez újmy na obecnosti uvažujme E=Ex,B=By. ma=αvtrˇenıˊ+FL Nechť je náboj silně tlumený, tedy maαv. Potom vFLα=qα(E+v×B) Elektrická část Lorentzovy síly ho potáhne střídavě nahoru a dolů, takže vlastně není moc zajímavá. Podívejme se na magnetickou část. Když o tom budeme uvažovat s použitím pravé ruky, zjistíme, že vždy bude tlačit náboj ve směru vlny. (Ještě tím, jak ho tlačí, mu dá složku rychlosti ve směru vlny, která trochu změní magnetickou část síly, ale to zanedbáme.) dpdt=FL=q(Exx+vxByzvzByx)=q(0+vxByz+0) dEdt=FLv=q(Ev+v(v×B))=q(Ev+B(v×v))=qExvx (pro maximalizaci zmatení je E tentokrát energie, nikoliv velikost E) Z toho plyne: dpdt=1cdEdtz Matematiko, odpusť mi mé hříchy… dpdV=1cdEdVz g=wcz To sice nevypadá jako ten obecný vztah, ale později se ukáže, že je to vlastně to samé.

Vztahy v rovinné elektromagnetické vlně

E pro změnu značí velikost E… aspoň doufám.

E=vB=Bεμ B=r×Ev w=12(εE2+B2μ)=εE2 S=E×Bμ=1μvE×(r×E)=εμE2r g=cwrc2=wrc g=Sc2(???)
Definice. Intenzita I|S|
Příklad. rz EE0xcos(ωtkz) I=vw=vεE2=εμE02x2cos(ωtkz)2=12εμE02

Tlak záření

Už máme hybnost g=Sc2. Taky máme tok energie S, kde dE=Sn. Kdo ví, jestli je to to samé S… Co bychom ještě chtěli? Co třeba tlak?

Představme si, že máme destičku, na kterou dopadá záření. Již víme, že má nějakou hybnost, takže bude na destičku působit silou a rozpočtením na plochu dostaneme tlak.

p=dFdS=dpdSdt

Pozor: p je tlak, ale dp je (diferenciálně malá) hybnost!

Dopadající hybnost je dp=gdV (z definice hustoty hybnosti). dV je objem záření, které dopadne za čas dt, tedy dV = cdtdS. Z toho máme

p=gcdtdSdSdt=gc=Sc=w

Pozor: S je jiné S než to předchozí S, které je v dS.

Příklad. Co kdybychom chtěli vědět, jak moc sluneční záření tlačí na Zemi? Na internetu asi někde najdeme, že S=1361Wm2, tedy p=4.5×106Pa. (Všimněme si, že jednotka S není m2, protože to není to S, které značí povrch. Chápeme?)

Energie vyzářené vlny pohybujícího se náboje

Již víme:

Erad=qa(tr)4πε0c2r tr=trc a=a(ar0)r0 B=Ec S=E×Bμ0=EBrμ0=E2rμ0c=ε0μ0E2r a=asin(θ)

Už se vzdáme myšlenky, že S je zároveň plocha a Poyntingův vektor, a budeme značit plochu A.

Správného fyzika teď bude zajímat celkový výkon protékající sférou o poloměru r:

dP=SrdA=SdA=ε0μ0E2 P=sfeˊradP=sfeˊraSdA=02π0πε0μ0E2r2sin(θ)dθdφ=ε0μ0r2r2(qa4πε0c2)20πsin(θ)3dθ=43ε0μ0(qa4πε0c2)2

„Přelitím“ konstant dostaneme Larmorův vzorec

P(t,r)=μ0q26πca(tr)2
Příklad (harmonický). xAcos(ωt)z a=Aω2cos(ωt)z P=μ0q26πcA2ω4cos(ωtr)2=μ0q26πcA2ω4cos(ωtkr)2 P=μ0q212πcA2ω4

Index lomu n

Uvažujme elektromagnetickou vlnu šířící se látkou. V té jsou nějaké atomy skládající se z těžkého jádra a lehkých elektronů navázaných na jádro. Vazba je inverzně kvadratická, takže si můžeme přédstavit, že elektron je k jádru přichycený pružinkou. Jelikož jádra jsou mnohem těžší než elektrony, jejich pohyb zanedbáme a budeme se zabývat pohybem elektronů. Ten je popsán pohybovou rovnicí, kde zároveň působí Lorentzova síla od dopadající vlny:

mx¨+mω02x=FL=e(E+v×B)

Obecný postup výpočtu: modelxpPεεrn

Budeme uvažovat nerelativistické rychlosti, tedy můžeme zanedbat magnetickou složku:

FBevB=evEc=βFEFE

Vezměme si konkrétně harmonickou vlnu:

E(t)=E0cos(ωt)

Vznikne buzený harmonický oscilátor. Obecný předpis pro něj je

x¨+2δx˙+ω02x=B^(t) B^=Bexp(iΩt) x^=Aexp(iΩt) A=Bω02Ω2+2iδΩ

V našem případě δ=0,Ω=ω, tedy

A=Bω02ω2

(ovšem B je budící síla, nikoliv magnetické pole!)

x(t)=emE0ω02ω2cos(iωt)=emE(t)ω02ω2

Dipólový moment jednoho atomu je p(t)=ex(t). Zajímá-li nás celý objem, chceme vektor polarizace (viz ELMA), což je objemová hustota dipólového momentu: PdpdV. Pokud zavedeme početní hustotu elektronů N, máme P=Np. Z elmy máme vzorec (vlastně spíš definici):

ε0E+P=εE ε0(1+Ne2mε0(ω02ω2))E=εE ε=(1+Ne2mε0(ω02ω2))ε0 εr=1+Ne2mε0(ω02ω2) n=εrμrεr=1+Ne2mε0(ω02ω2)

Jelikož různé elektrony můžou mít různé ωk,Nk, máme obecně

n1+kNke2mε0(ωk2ω2)

(Aby bylo jasno: ωk je úhlová rychlost pružinky, na které poskakuje atom, ω je úhlová rychlost vlny.)

Vyšlo nám tedy, že n závisí na ω (což například umožňuje rozklad barevného spektra na hranolu).

Pro fázovou rychlost máme v=cn. Pokud ω<ω0, potom n>1, tedy světlo má v této látce menší rychlost než ve vakuu. Je-li ω mezi ω0 a ω1ω02+Ne2mε0, vychází index lomu záporný, takže v takové látce se světlo vůbec šířit nemůže. A nakonec pokud ω>ω0, bude n<1, takže fázová rychlost světla bude větší než ve vakuu — ale to nevadí, protože fázovou rychlostí se nemůže šířit informace. (Analogie: Pokud si představíme břeh, na který téměř kolmo dopadá vlna, zdánlivá rychlost vlny na břehu bude vb=vcos(α), což taky může teoreticky přesáhnout rychlost světla.)

Disperzní vztah v plazmatu

Představíme si plazma jako volné elektrony (na kladné ionty se vykašleme, protože jsou skoro nehybné), tedy ω0=0. Dosazením do vzorečku dostaneme

n=1Ne2mε0ω2

Z disperzního vztahu ω=cnk vymlátíme n2=c2k2ω2, tudíž

1Nc2mε0ω2=c2k2ω2 ω2=Nc2mε0+c2k2ωp2+c2k2

Pokud ω>ωp, najdeme přípustné k, tedy prostředí je transparentní. Pokud však ω<ωp, prostředí je reaktivní.

Vlnovod

Chceme nechat šířit elektromagnetické pole bez úbytku amplitudy. K tomu použijeme vlnovod — dutou kovovou trubku.

Pro šíření elektrického pole známe vlnovou rovnici a divergenční rovnici:

2Et2=c22EE=0

Jak z toho uděláme něco úplně jiného? Přece tak, že přidáme okrajové podmínky!

Představme si, že jsou stěny dokonale vodivé, tedy elektrony se v nich vždy postaví tak, aby vyrušily elektrické pole. To můžeme odvodit z diferenciálního Ohmova zákona:

j=σE E=jσσ0

Představme si vlnovod jako dlouhou obdélníkovou trubici ve směru z. Nechť je pro dolní stěnu x=0, horní stěnu x=a, levou stěnu y=0, pravou stěnu y=b. Takže máme okrajové podmínky

E(x{0,a},y,z,t)=E(x,y{0,b},z,t)=0

Nezajímá nás obecné řešení, ale jenom pro vlnu šířící se ve směru z.

Mohli bychom zkusit zanedbat x a y, ale z toho by nevyšlo nic zajímavého, takže se vykašleme jenom na x. To znamená, že Ey(x,y,z,t)=Ey(0,y,z,t)=0, podobně pro z, takže jediná nenulová složka E je x. Dosadíme tedy do rovnic a dostaneme:

E=Exx=0=maˊ byˊt0 2Ext2=maˊ byˊtc2(2Exy2+2Exz2)

Předpokládejme harmoničnost ve směru z: Ex(y,z,t)=Y(y)exp(i(ωtkz))

ω2Ex=c2(Y(y)exp(i(ωtkz))kz2Ex) ω2Y=c2(Yk2Y) Y+(ω2c2kz2)KY=0

Pokud K0, dokážeme to snadno vyřešit a zjistíme, že funguje jenom triviální řeséní. Nechť tedy K>0,kyK. Máme z toho rovnici harmonického oscilátoru, tedy

Y(y)=αcos(kyy)+βsin(kyy)

Pokud se ještě zamyslíme nad okrajovými podmínkami, zjistíme, že Y(0)=Y(b)=0.

A teď už to bude úplně stejný jako pro tu strunu! —Josef Schmidt, 2022
ky(m)mπb Y(m)(y)=sin(ky(m)y) Ex(m)(y)=E0sin(ky(m)y)exp(i(ωtkzz))

Jelikož má platit ky=ω2c2k2 a zároveň ky=mπb, dostáváme disperzní vztah

ω2=(mcπb)2+(ckz)2ωmin(m)2+c2k2

(kde kkz, protože vlna je koneckonců ve směru z)

Pokud si představíme řez vlnovodem, tak v prvním módu vznikne jakási dvourozměrná postupná vlna skládající se z kopců, ve druhém módu budou dva kopce vedle sebe, a tak dále.

Zjistili jsme, že pro každý mód je nějaká minimální frekvence ωmin(m), která se může transparentně šířit. Pokud bude ω<ωmin, můžeme označit kiϰ, tedy ϰωmin2ω2c, a budeme mít tlumenou vlnu (tedy reaktivní mód):

Ex=E0sin(mπyb)exp(κz)exp(iωt)

Pokud si spočítáme fázovou a grupovou rychlost, vyjde nám vφ=ωkz>c,vg=dωdkz<c.

Podmínky napojení elektromagnetického pole na rozhraní dvou látkových prostředí

Mějme rozhraní dvou prostředí. Vlevo máme ε1,μ1,E1(r,t),B1(r,t). Vpravo máme ε2,μ2,E2(r,t),B2(r,t).

Samozřejmě použijeme Maxwellovy rovnice:

(εE)=0 B=0 ×E=Bt ×Bμ=(εE)t

Problém je, že na rozhraní to úplně nefunguje, protože tam funkce nejsou diferencovatelné, takže tam musíme použít integrální tvary:

SεEdS=0 SBdS=0

Dejme kolem počátku válec mířící ve směru, ve kterém se mění prostředí.

S=Sp1+Sp2+Spl

Dělat to pořádně by bylo na dlouho, takže to pojďme dělat nepořádně a rovnou posílat h0,Sp0.

Φp1=ε1E1rSp Φp2=ε2E2rSp Φpl=α2πrh

Z toho nějak

ε1E1=ε2E2 B1=B2

(tím se myslí složky kolmé na rozhraní)

Teď by nás ještě zajímaly rovnoběžné složky. No když ty kolmé jsme odvodili přes divergenční rovnice, tak to takhle tentokrát použít rotační rovnice (opět v integrálním tvaru)?

lBμdl=ddtSεEdS=ddtSεE×ndS=strˇednıˊhodnotaαlh

Jako křivku zvolme obdélník, jehož dvě strany procházejí kolmo rozhraním.

l=levaˊ+hornıˊ+pravaˊ+dolnıˊh0levaˊ+pravaˊ=levaˊB1μ1tdlpravaˊB2μ2tdl=maˊ byˊt0 B1μ1=B2μ2 E1=E2

Tím jsme dostali informaci jen o jednom rovnoběžném směru, tedy musíme obdélník ještě pootočit a udělat to znova!

Polarizace

Elektromagnetické vlnění je kolmé na směr šíření, tedy máme k dispozici dva nezávislé směry, kam může směřovat. Nechť směr šíření je z, potom báze elektromagnetického vlnění je (x,y). Předpis pro obecnou harmonickou elektrickou vlnu je tedy

E=Ex0xexp(i(ωtkz+φ1))+Ey0yexp(i(ωtkz+φ2))
Věta. Fázový rozdíl nezávisí na z,t.

Tedy fáze se mění hodně rychle, ale amplitudy a fázový rozdíl zůstávají pořád stejné.

Použijme takový trik. Posuňme si obě fáze o stejnou hodnotu, konkrétně φ1, tedy místo φ1,φ2 budeme mít 0,φ2φ1, což jsou konstanty.

Zvolme nyní fixní z a zahrňme podobným způsobem do fáze kz. Tím dostaneme

E(t)=Ex0xexp(i(ωt+φ1))+Ey0yexp(i(ωt+φ2))

Tohle je parametrická rovnice elipsy! Tedy vektor v každém kolmém řezu vlny opisuje elipsu. Říkáme, že vlna je elipticky polarizovaná.

Jaká je intenzita?

I=S=ε0μ0E2=ε0μ0Ex02x2cos(ωt+φ1)2+Ey02y2cos(ωt+φ1)2+2Ex0Ey0(xy)cos(ωt+φ1)cos(ωt+φ2)=ε0μ0Ex02cos(ωt+φ1)2+Ey02cos(ωt+φ1)2=ε0μ0Ex02+Ey022

Můžeme rozlišit dva speciální případy:

Pravo- a levotočivost samozřejmě funguje i u elipsy. Jenom v případě lineární polarizace se nerozlišuje.

Lineární polarizátor

Lineární polarizátor mění světlo na lineárně polarizované. Má nějakou osu propustnosti n a propouští pouze složku vlny rovnoběžnou s tímto směrem. Ein=E+E Eout=E=(Einn)n

Vsuvka: Jonesův kalkulus

E(t)=Ex0xexp(i(ωt+φ1))+Ey0yexp(i(ωt+φ2))=(Ex0xexp(iφ1)+Ey0yexp(iφ2))E^exp(iωt) E^(Ex0exp(iφ1)Ey0exp(iφ2))
Jonesův vektor pro lineární polarizátor bude E^out=𝐏nE^in (kde 𝐏 je projektor). Eout=(Einxnx+Einxny)(nxny)=(Einxnx2+EinynxnyEinxnxny+Einyny2)=(nx2nxnynxnyny2)𝐏Ein

Malusův zákon

Mějme dva lineární polarizátory za sebou, chceme zjistit intenzitu výstupu z druhého vůči vstupu do něj (jímž je vlna lineárně polarizovaná od prvního polarizátoru).

Nechť n1x,n2n(cos(θ)sin(θ)).

U druhého polarizátoru máme:

Ein=E0xcos(ωt) Iin=Ein2=E022 Eout=(Einn)n=E0cos(θ)ncos(ωt) Iout=Eout2=E022cos(θ)2=Iincos(θ)2

Vlnová destička

Fázový rozdíl δφ se normálně nemění, ale vlnová destička ho dokáže změnit o nějaké Δφ.

Nechť d je tloušťka destičky. Kdyby byla z normálního materiálu, bylo by:

Eout=E0xxcos(ωt+φ1kd)+E0yycos(ωt+φ2kd)

Tedy by fázový rozdíl neměnila. Ale my ji můžeme vyrobit z dvojlomného krystalu, což je materiál, který má v rovnoběžném směru jiný index lomu než v kolmém směru. Potom:

Eout=E0xxcos(ωt+φ1k1d)+E0yycos(ωt+φ2k2d) δφ=(ωt+φ1k1d)+(ωt+φ2k2d)=φ1φ2+(k1k2)dδφ+Δφ Δφ=(k2k1)d=ωc(n1n2)d=2πλ0(n2n1)d

Pokud Δφ=π2, jde o čtvrtvlnovou destičku. Pokud Δφ=π, jde o půlvlnovou destičku.

Pro zjednodušení můžeme o oměma fázím přičíst k2d, čímž dostaneme

Eout=E0xxcos(ωt+φ1(k2k1)d)+E0yycos(ωt+φ2)

Z toho vzniká takové pravidlo:

Máš-li vlnovou destičku s Δφ>0, přičti Δφ ve složce odpovídající menšímu indexu lomu.

Nyní aplikujme Jonesův kalkulus:

E^out=𝐃ΔφE^in E^in=E^n1+E^n2=𝐈E^in=(𝐏n1+𝐏n2)E^in=𝐏n1E^in+𝐏n2E^in E^out=(exp(iΔφ)𝐏n1+𝐏n2)E^in 𝐃=exp(iΔφ)𝐏n1+𝐏n2

Optická aktivita

Tahle podsekce vůbec nedává smysl a ani není ve skriptech, takže ji možná radši zanedbejte.

Vlnová destička přidává fázi do x-ové složky. Co kdyby nějaký optický prvek přidával fázi do obou složek?

Ein=Exxcos(ωt)=E02(xcos(ωt)+ysin(ωt))+E02(xcos(ωt)ysin(ωt)) Eout=Exxcos(ωt+Δφ)=E02(xcos(ωt+Δφ)+ysin(ωt+Δφ))+E02(xcos(ωt+Δφ)ysin(ωt+Δφ))==(E0cos(Δφ2)x+E0sin(Δφ2)y)cos(ωt+Δφ2)=E0ncos(ωt+Δφ2)

Intuice za tímhle je taková, že když se sečtou dva vektory točící se v opačném směru, vznikne tím vektor mířící na konstantní stranu a měnící sinusově délku, přičemž směr, kam míří, je polovina fázového rozdílu.

Nepolarizované světlo

Zatím máme pouze popis polarizovaného světla.

Máme dvě časové škály:

Měření polarizace

E=Ex0xcos(ωt+φ1)+Ey0ycos(ωt+φ2)

Měřicí přístroj má nějakou rozlišovací dobu trozT

I=E2troz

Definujme další čtyři intenzity:

IxEx2=Ex02cos(ωt+φ1)2=Ex022 IyEy2=Ey02cos(ωt+φ1)2=Ey022 IxyExEy=Ex0Ey0cos(ωt+φ1)cos(ωt+φ2)==Ex0Ey0cos(φ1φ2)2 IxyEx(ωtπ2)Ey=Ex0Ey0cos(ωt+φ1)cos(ωt+φ2)==Ex0Ey0sin(φ1φ2)2

Tedy měřením Ix dokážeme zjistit Ex, analogicky Ey, a měřením Exy,Exy zjistíme δφ. Ale jak je změříme?

U Ix,Iy je to snadné. Pokud dáme světlu do cesty polarizátor pod úhlem π4, bude výsledná intenzita Ix+Iy2+Ixy, z čehož snadno dopočteme Ixy. Pokud ještě před něj přidáme čtvrtvlnovou destičku, dostaneme podobně Ixy.

A co nepolarizované světlo? Teď máme už tři časové škály: periodu T1014s, koherenční dobu tkoh108s a rozlišovací dobu troz. Pokud troztkoh, jde o rychlý přístroj a můžeme prostě sledovat změny polarizace. Horší případ je troztkoh — „pomalý přístroj“? Jelikož Ttkoh, budou již střední hodnoty závislé na čase:

Ix(t)=Ex2T=Ex0(t)22 Iy(t)=Ey2T=Ey0(t)22 Ixy(t)=T=Ex0(t)Ey0(t)cos(δφ(t))2 Ixy(t)=T=Ex0(t)Ey0(t)sin(δφ(t))2

Ovšem intenzity, které naměříme, jsou Ixprˇıˊstroj=Ix(t)troz,

Obecně se Ex0,Ey0,φ1,φ2 mění tak, že rovnoměrně pokrývají polarizační stavy, ale celková intenzita je konstantní. Ze symetrie plyne, že Ix=Iy=I2. Také ze symetrie máme cos(δφ)=sin(δφ)=0, z čehož takovým fyzikálním podvodem můžeme vyvodit, že Ixy=Ixy=0.

Zřejmě pro nepolarizované světlo P=0. Pro polarizované světlo:

|P|=(IxIy)2+2Ixy2+2Ixy2(Ix+Iy)2=Ex042Ex02Ey02+Ey04+4Ex02Ey02cos(δφ)2+4Ex02Ey02sin(δφ)2Ex04+2Ex02Ey02+Ey04=1

Pro částečně polarizované světlo vyjde něco mezi. Tedy Stokesovy parametry tvoří jakousi kouli s nepolarizovaným světlem ve středu a polarizovaným světlem na okraji.

Z hlediska měření jsou tedy všechna nepolarizovaná světla ekvivalentní.

Nepolarizované světlo nemusí nutně být úplně náhodné. Pokud bychom například pravidelně střídali lineární polarizaci vodorovně a svisle, taky by všechno vycházelo, tedy Ix=Iy=12I0,Ixy=Ixy=0. Nulovost Ixy vyplývá z toho, že v každém okamžiku je Ex0=0Ey0=0. Nulovost Ixy vyplývá z toho, že δφ=0, protože jde o lineární polarizaci. To ale neznamená, že nepolarizované světlo můžeme vytvořit superpozicí polarizovaných světel, musíme je skutečně střídat!

Odraz elektromagnetické vlny na rovinném rozhraní

Mějme rovinu ve směru xy. Na rovinu dopadá Ed, odrazí se Er, projde Et. Uvažujme harmonickou vlnu E=Ed0exp(i(ωtkdr)),kd=kdnd,ω=cn1kd.

Předpokládejme, že Er=exp(i(ωrtkrr)),kr=krnr,ωr=cn1kr a Et=exp(i(ωttktt)),kt=ktnt,ωt=cn2kt.

Natočíme osy tak, aby kd=(0,kdy,kdz). Tím zajistíme, že Ed nezávisí na x, tedy vlna má translační symetrii podél osy x. Pokud nedokážeme poznat složku x dopadající vlny, nedokážeme ji poznat ani u odražené a prošlé vlny, tedy i ty nebudou záviset na x: kr=(0,kry,krz),kt=(0,kty,ktz). Na x tedy vůbec nezáleží — jde o rovinný problém.

Obrázek situace
kd=(0,kdsin(ϑd),kdcos(ϑd)) kr=(0,krsin(ϑr),krcos(ϑr)) kt=(0,ktsin(ϑt),ktcos(ϑt))

Z Faradayova zákona vyplývá první podmínka napojení:

E1(z=0)=E2(z=0) Ed(z=0)+Er(z=0)=Et(z=0) Ed0exp(i(ωtkdsin(ϑd)y))+Er0exp(i(ωtkrsin(ϑr)y))=Et0exp(i(ωtktsin(ϑt)y))

Exponenciály se musí rovnat ve všech bodech. Tedy kdyby byly lineárně nezávislé, jediné řešení by byly nulové amplitudy, takže by tam žádná vlna nebyla, což nechceme! Aby byly lineárně závislé, musí se rovnat, tedy ω=ωr=ωt,kdsin(ϑd)=krsin(ϑr)=ktsin(ϑt).

Jelikož ωd=ωr, z disperzního vztahu plyne kd=kr, tedy ϑd=ϑr. To je zákon odrazu.

Poku naopak použijeme ωr=ωt, z disperzního vztahu máme krn1=ktn2, tedy n1sin(ϑd)=n2sin(ϑr). To je Snellův zákon lomu.

Aby vůbec existoval vhodný úhel průchodu, musí být úhel dopadu menší než kritický úhel ϑcarcsin(n2n1). Pokud bude větší, náš předpoklad již nefunguje! Mohli bychom předpokládat, že se potom všechno odrazí, ale to je tak trochu lež.

Nechť ϑd>ϑc, tedy jde o totální odraz. To znamená, že druhé prostředí se chová jako reaktivní.

ω2=c2n22|k1|2=c2n22(kty2+ktz2)=c2n22(kdy2+ktz2)=c2n22(kd2sin(ϑd)2+ktz2)=c2n22(ω2c2n12sin(ϑ1)2+ktz2) ϰ2ktz2=ω2c2(n22n12sin(ϑd)2) Et=Et0exp(i(ωtktr))=Et0exp(i(ωtktyy))exp(ϰz)

Vznikne vlna, která se šíří ve směru y a exponenciálně ubývá ve směru z.

Tak tedy víme, kterými směrem vlna půjde, ale jaká bude intenzita?

Vlna dopadá po nějaké rovině dopadu kolmé na rozhraní. Rozdělme její polarizaci na složku v rovině dopadu () a složku kolmou na rovinu dopadu ().

Postup je podobný, jako když jsme něco podobného dělali pro strunu, akorát složitější, protože struna je lineární a elektromagnetická vlna je rovinná.

Nějaké technické výpočty, které si nemáme zapisovat, protože je Schmidt stejně opisuje ze skript.

V nějakou chvíli se ϑd,ϑt přejmenovaly na ϑ1,ϑ2.

Vyjde nějaký přesný výsledek, který se nijak nejmenuje, ten se dá zaokrouhlit pomocí aproximace μ1,2μ0 na tzv. Fresnelovy vzorce:

Rn2cos(ϑ1)n1cos(ϑ2)n2cos(ϑ1)+n1cos(ϑ2)Rn1cos(ϑ1)n2cos(ϑ2)n1cos(ϑ1)+n2cos(ϑ2)

(kde R je definováno podobně jako u struny)

Když použijeme vztah n1n2=sin(ϑ2)sin(ϑ1), dostaneme zjednodušené Fresnelovy vzorce:

Rtan(ϑ2ϑ1)tan(ϑ2+ϑ1)Rsin(ϑ2ϑ1)sin(ϑ2+ϑ1)

Při odrazu se světlo částečně polarizuje, protože jedna polarizace se odráží víc. Při tzv. Brewsterově úhlu se polarizuje úplně. Jak ho určíme? Chceme, aby bylo R=0. Toho můžeme u obou zlomků docílit tím, že čitatel bude 0 nebo jmenovatel bude . Nastavit čitatele na 0 nám nepomůže, protože potom bude ϑ1=ϑ2 a po dosazení do původního výrazu zjistíme, že vlastně nevyjde nula. Sínus nemůže jít do nekonečna, takže jediná naděje je tan(ϑ2ϑ1)=0ϑ2ϑ1=±π2. Tedy odražený a prošlý paprsek mají svírat pravý úhel a ϑB=arctan(n2n1).

Brewsterův úhel se dá využít k polarizaci odrazem nebo k tzv. Brewsterovu oknu (pokud pošleme správně polarizované světlo pod správným úhlem, nic se neodrazí, tedy materiál se bude chovat jako dokonale průhledný).

Zkouška

V pátek 16. 12. bude možná předtermín

Druhý zápočtový test se píše ve středu 14. 12. místo přednášky

V zápočtovém testu příklady na 15, 10, 10, 5 bodů, nutno získat 20

Na limity klasické fyziky už se nezkouší

Nejdůležitější téma je polarizace

Interference

Maxwellovy rovnice jsou lineární, tedy splňují princip superpozice. Ovšem v intenzitě už se vyskytují druhé mocniny, takže lineární není.

I=E2=(E1+E2)2=E12+E22+2E1E2interferencˇnıˊcˇlen

(Teď už ve vzorci pro intenzitu kašleme na ε0μ0, protože je k ničemu.)

Rozdělíme paprsek rozdělovačem paprsků a získáme interferometr, který normálně vypadá takto (Z = zdroj, B = rozdělovač, M = zrcadlo, D = stínítko):

D B Z M M

Pro zjednodušení můžeme uvažovat tzv. pedagogický interferometr, který je jednodušší, ale nedá se reálně sestrojit, protože nejsme schopni zařídit, aby vysílané paprsky byly přesně stejné:

D B Z2 M Z1

Budeme považovat elektrické pole za skalární, protože vektorovost teď není důležitá. Pole budeme uvažovat se stejnou frekvencí, ale jinými amplitudami. Elektrické pole ze zdroje 1, 2 je EZ1,2(t)E1,2cos(ωt+φ1,2).

To, co z jednotlivých zdrojů dopadne, je to samé, ale v retardovaném čase, tedy:

ED1,2(t)EZ1,2(tr)=E1,2cos(ωtkl1,2+φ1,2)

Co vlastně chceme zjistit, je dopadající intenzita:

ID=(E1D+E2D)2=E12cos()2+E22cos()2+2E1E2cos(ωtkl1+φ1)cos(ωtkl2+φ2)=E122+E222+E1E2(cos(2ωtk(l1+l2)+φ1+φ2)+cos(k(l2l1)+φ1φ2))=I1+I2+2I1I2cos(k(l2l1)+φ1φ2)I1+I2+Ii

Fáze interference je tedy posunuta o Δφ=k(l2l1)=2πclλ. Všechno kromě λ známe, takže tímto způsobem můžeme měřit vlnovou délku.

Jelikož cos je na intervalu [1,1], máme:

Id[I1+I22I1I2,I1+I2+2I1I2]=[(I1I2)2,(I1+I2)2]

Pokud by vysílané intenzity byly stejné, tedy I1=I2I0, potom ID[0,4I0], konkrétně I0=2I0(1+cos(Δφ)).

Tohle je všechno celkem jednoduché, ale vlastně to není úplně pravda. Reálně nejsou φ1,2 konstanty.

Uvažujme podobně jako dříve, že fáze φ1,2(t) jsou konstantní po nějakou koherenční dobu a potom se náhodně změní. Jak to ovlivní náš výpočet? Zaprvé musíme do výrazu pro ED přidat závislost na retardovaném čase:

ED1,2(t)=E1,2cos(ωtkl1,2+φ1,2(tr1,2))

Zadruhé už při středování nemůžeme uvažovat periodičnost. To se dá vyřešit tím, že budeme středovat jen přes jednu periodu, protože ta je mnohem kratší než koherenční doba.

I(t)(E1D2+E2D2)T=I1+I2+2I1I2cos(k(l2l1)+φ1(tr1)φ2(tr2)) ID=I(t)troz
Definice. Pole je prostorově koherentní v bodech A,B, pokud ze znalosti pole v A lze předpovědět pole v B ve stejném čase.
Definice. Pole je časově koherentní v časech t1,2, pokud ze znalosti pole v t1 lze předpovědět pole v t2 ve stejném bodě.

Uvažujme několik podpřípadů:

Reálně máme jen pomalé přístroje, z čehož plyne poučka: „Nekoherentní vlny spolu neinterferují“. (Tím se myslí, že po vystředování v přístroji je interferenční člen nulový.) Proto normálně nevidíme všude interferenci.

Pokud vlny nedopadnou přesně ve stejném směru, ale trochu našikmo, vzniknou tím známé interferenční proužky. (Prý se to počítalo na cvičení.)

Difrakce

Již víme, jak spolu interferují dvě vlny, ale co kdyby jich bylo nekonečně mnoho?

Mějme zdroj světla. Umístěme světlu do cesty překážku ve tvaru roviny s malým otvorem. Následovat bude další rovina, tentokrát bez otvoru – stínítko. Zajímá nás rozložení intenzity na stínítku.

Z

Provedeme sérii jednoduchých úvah.

Babinetův princip

Uvažujme neprůhlednou přepážku. Na ni bude dopadat pole Edop a ona bude „indukovat“ pole Eind. Z definice neprůhlednosti musí být Edop+Eind=0.

Z

Přepážku si rozdělíme na části A,B. Zřejmě Eind=EindA+EindB.

Z

Nyní dejme pryč B. Za přepážkou již bude nějaké pole EA. Jelikož předtím tam nebylo nic, musí být EA=EindB. Tedy to, co projde, je opak toho, co by se předtím indukovalo na B.

Tato úvaha není úplně korektní, protože ve skutečnosti atomy obou částí jsou nějak propojené, takže odstraněním B může dojít ke změně EindA. Ale pokud je otvor dostatečně velký (Dλ), můžeme tyto efekty zanedbat.

Dalšimi hrátkami můžeme dostat vzorce:

EA+EB=(EindA+Eindb)(Babinetu˚v princip) EA=EdopEB

Difrakční integrál

Předpokládejme:

Každý malý úsek přepážky vyzařuje vlnu E0rexp(i(ωtkr)). Celkově tedy máme

EA=EindB=BE0rexp(i(ωtkr))dS

Tím bychom mohli říct, že je úloha vyřešena, ale výraz uvnitř integrálu je vlastně dost hnusný, takže ještě provedeme nějaké aproximace.

Zaveďme v rovině přépážky souřadnice X,Y a v rovině stínítka souřadnice x,y. Nechť L je vzdálenost mezi přepážkou a stínítkem. Pokud si vybereme bod na přepážce a bod na stínítku, vzdálenost mezi nimi bude

r=L2+(Xx)2+(Yy)2

Integrál tedy bude vypadat takto:

EA(x,y)=SE0r(x,y,X,Y)exp(i(ωtkr(x,y,X,Y)))dXdY

Navíc jsme zanedbali, že vyzařování je anizotropní, což můžeme vyřešit tím, že uvažujeme jen malou část stínítka. Také z našich úvah nelze určit E0, tedy jak moc vlna rozkmitá atomy; na to by byl potřeba složitější model. To nám ale nevadí, protože chceme stanovit pouze relativní intenzity.

Historicky to vychází s Huygensova-Fresnelova principu, kdy ještě neměli model atomů a jen na základě experimentů postulovali, jak to funguje. Nejdříve Huygens zjistil, že z otvoru jdou sférické vlnoplochy, poté Fresnel navrhl, že u nich funguje princip superpozice.

Zaveďme vzdálenost RL2+x2+y2, která bude ze středu díry (pro jednoduchost v počátku) do bodu na stínítku. Na rozdíl od r nezávisí na integračních proměnných, takže bychom ji chtěli nějak nacpat do integrálu. Zjistěme, jaký je mezi nimi rozdíl:

r2=R2+(Xx)2+(Yy)2x2y2=R2+X2+Y22xX2yY

Použijme obvyklý fyzikální postup – vezmeme nejhrubší aproximaci, při které jev úplně nezmizí. Vybereme tyto aproximace:

r=R1+X2+Y22xX2yYR2R(1xX+yYR2)

Ještě hrubší aproximace je rR nebo dokonce rL. Ovšem tyto dvě nemůžeme dosadit do integrálu na obou místech, protože potom by to byl integrál z konstanty a nevyšlo by nic zajímavého. Rozhodneme se tak, že do jmenovatele za r dosadíme R a do kosínu lineární vztah.

EA(x,y)SE0Rexp(i(ωtk(RxX+yYR)))dXdY=E0Rexp(i(ωtkR))Sexp(ikR(xX+yY))dXdY

Pokud definujeme f(X,Y)[(X,y)B], můžeme integrovat přes celé 2, akorát s tím, že to ještě vynásobíme f(X,Y). To nás povede na dvourozměrnou Fourierovu transformaci.

Kdybychom vrátili zpět ten kvadratický člen, který jsme zanedbali, přibyl by do fáze a způsobil by fázový rozdíl. Aby tato aproximace byla dobrá, musí tento rozdíl být malý, čímž vzniká podmínka Fraunhoferovy difrakce:

ΔφkX2+Y22R1

Díra v přepážce je podmnožina koule o průměru D, takže platí X2+Y2D24. Také z definice máme RL. Z toho plyne

Δφ2πλD2412L=πD24λLD2λLmaˊ byˊt1

Aproximaci tedy můžeme provést, pokud platí zjednodušená podmínka:

D2λL

Jaký je geometrický význam kxX+yYR? Z každého bodu díry vychází vlna a dopadá do bodu na stínítku. Zajímá nás, jakou dráhu při tom urazí. Zaveďme úhel θ na spojnici mezi středem díry a bodem na stínítku, tedy sin(θ)=xR. Pokud si nakreslíme obrázek s paprskem od jiného bodu díry, zjistíme, že zároveň sin(θ)=lX, kde l je kolmý průmět bodu díry na spojnici mezi středem díry a bodem stínítka. Z toho plyne l=xXR, tedy našli jsme geometrickou podstatu výrazu. Máme lr0rX, kde r0,rX jsou délky spojnic.

𝐿 𝑋  𝑟₀  𝑥  θ 𝑟ₓ 𝑙  

Youngův pokus

Mějme dvě obdélníkové štěrbiny. Zjednodušíme si je na dva bodové otvory, čímž se z integrálu stane jednoduše součet dvou čísel. Z Fraunhoferova integrálu máme

E=E0Rexp(i(ωtkR))(exp(ikRxd2)+exp(ikRxd2))=2E0Rexp(i(ωtkR))cos(kdxR)=2E0Rexp(i(ωtkR))cos(kdsin(θ))

V reálné části se z exp stane cos, což nám umožní spočítat intenzitu:

I=E2=E02R242cos(kdsin(θ)2)2=2E02R2cos(πdsin(θ)λ)2

Intenzita tedy bude mít maxima v bodech, kde πdsin(θ)λ=mπ,m. Pokud bude stínítko daleko od přepážky, budou tím vznikat přibližně rovnoměrně rozmístěné kopce ve vzdálenosti Lλd. Každý kopec má nějaký řád m. Aby bylo |sin(θ)|<1, musí být |m|<dλ, což je tedy maximální řád maxim.

Difrakční mřížka

Jak bychom to mohli zesložitit? Co takhle N bodových štěrbin? Nechť jsou rovnoměrně rozmístěny na úsečce ve vzdálenostech 0,d,2d,,Nd.

E=Fraunhoferu˚vvzorecE0Rexp(i(ωtkR))j=0N1exp(ikRxjd)=geometrickaˊrˇadaE0Rexp(i(ωtkR))exp(ikRxdN)1exp(ikRxd)1=vynaˊsobenıˊjednicˇkouE0Rexp(i(ωtkR))exp(ikRxdN)1exp(ikRxd)1exp(ikRxdN2)exp(ikRxdN2)exp(ikRxd2)exp(ikRxd2)=sin(x)=exp(ix)exp(ix)2E0Rexp(i(ωtkRkRxdN12))sin(kdNx2R)sin(kdx2R) I=E2=E022R2sin(kdNsin(θ)2)2sin(kdsin(θ)2)2

Opět chceme najít maxima. Ta máme v bodech, kde je nulový jmenovatel. Takže zjistíme, že to vychází úplně stejně jako u Youngova pokusu. Co se ovšem změnilo, je šířka maxim. Ta závisí na tom, kde je čitatel nulový, což jsou násobky λNd. Tedy čím více máme štěrbin, tím užší jsou maxima. Pokud definujeme šířku maxim jako vzdálenost mezi nulami čitatele, vyjde 2λNd.

To můžeme použít k vytvoření spektrometru. Ovšem pokud budou dvě vlnové délky příliš blízké, jejich kopce se slijí, takže nepůjdou moc dobře rozlišit. Aby k tomu nedošlo, požadujeme:

|sin(θ1)sin(θ2)|2λNd m|λ1λ2|λ1+λ2N Nλ1+λ2m|λ1λ2|

To nám určuje, kolik děr musíme mít, abychom rozpoznali spektrum.

Obdélníkový otvor

Co zkusit něco složitějšího? Mějme otvor ve tvaru obdélníku. Ve směru x má šířku a, ve směru y výšku b. Na cvičení zjistíme, že vyjde

Eabsin(kax2R)kax2Rsin(kby2R)kby2R=y=0absin(kasin(θ)2)kasin(θ)2 Isin(kasin(θ)2)2(kasin(θ)2)2(sin(α)α)2

Máme nějakou úhlovou šířku maxima Δθ=2λd, ta se nazývá úhlová rozbíhavost svazku.

Kruhový otvor

Pro jednoduchost použijeme polární souřadnice: dXdYρdρdφ. Ze symetrie můžeme očekávat, že indenzita bude záviset pouze na r (což naní to samé r jako předtím, ale vlastně něco jako x, akorát v polárních souřadnicích). Bude z toho extrémně hnusný výpočet, z něhož vyjde tohle:

I(θ)=(J1(kdsin(θ)2)kdsin(θ)2)2

Jn je Besselova funkce, kterou si lze představit jako jakýsi ubývající sínus.

Jn(x)1π0πcos(nuxsin(u))du

Maxima teď budou u0=0,u11.22π.

Skutečný Youngův pokus, tedy dvě obdélníkové štěrbiny

Nechť štěrbiny mají velikost d a vzdálenost mezi nimi je d. Potom budou vznikat jakési dvě difrakce, které se navzájem vynásobí, takže výsledkem bude vlna s vysokou frekvencí zabalená ve vlně s nízkou frekvencí (něco jako u vlnových balíků).

Kritérium rozmazanosti

Proč v běžném životě nepozorujeme difrakci? No protože je to zase interferenční jev, takže tam záleží na koherenci. Koherenční doba je tkoh109s, tedy koherenční délka je lkoh=ctkoh30cm. Interference je viditelná, pokud dráhový rozdíl je menší než koherenční délka, a dráhový rozdíl je dsin(θ). U difrakce máme díry mnohem menší než 30 centimetrů, takže časová koherence nám nevadí.

Horší je ale prostorová koherence. Pokud nám na díru svítí více zdrojů světla, které jsou spolu vzájemně nekoherentní, jejich difrakční obrazce se superponují. Pokud budou všechny na ose, tak nám to moc nevadí, protože obrazce vzniknou stejné. Ale co když budou mimo osu? Tím vznikne jakýsi našikmený obrazec.

Pokud se maximum šikmého obrazce zrovna náhodou trefí do prvního maxima rovného obrazce, pořád se navzájem posílí. Tedy pokud jsou zdroje světla navzájem pod úhlem ϑ, máme kritérium rozmazanosti: ΔϑΔθ. Nebo taky můžou být blízko sebe a trefit se do stejného maxima, to je za podmínky ΔϑΔθ.

Geometrická optika

O té nemá cenu mluvit, protože se to všechno dá odvodit z vlnové optiky.