Adamátor Zápisky Hlášky

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=−∂B→∂t(Faradayu˚v zaˊkon) ∇→×B=εμ∂E→∂t(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→)=−εμ∂2E→∂t2

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

∂2E→∂t2=1εμ∇2E→

Analogicky lze odvodit

∂2B→∂t2=1εμ∇2B→

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

v≔1εμ c≔1ε0μ0(ve vakuu) n≔cv=εμε0μ0(index lomu)

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

E→=E→0⋅F(k→⋅r→−vt) B→=B→0⋅F(k→⋅r→−vt)

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(njxj−vt))=E0iniF′(n→⋅r→−vt)=maˊ byˊt0 E→0⋅n→=0∴E→0⟂r→

Analogicky

B→0⋅n→=0∴B→0⟂r→

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→=−∂B→∂t εijk∂jEk=−∂tBi εijk∂j(E0kF(nlxl−vt))=−∂t(B0iF(nlxl−vt)) εijkE0knjF(n→⋅r→−vt)=−B0i(−v)F(n→r→−vt) n→×E→0=vB→0

Z toho plyne:

Příklad (harmonická elektromagnetická vlna). F(x)≔exp⁡(ikx) E→=E→0exp⁡(−i(ωt−k→⋅r→)) B→=B→0exp⁡(−i(ωt−k→⋅r→)) ω=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→=E→0exp⁡(−i(ωt−kz)) B→=B→0exp⁡(−i(ωt−kz))

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)=qr→′4πε0r′3,r→′≔r→−v→t

Uvažujme nyní, že v čase t≤0 je náboj v klidu, v čase 0≤t≤τ 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:

∮SE→⋅dS→=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=∮SE→⋅dS→=∮podstavyE→⋅dS→

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=EstatSp−E∥Sp E∥=Estat

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

cτv⟂t=podobnost △E∥E⟂ E→⟂=v⟂tcτE∥=a⟂rc2E∥ E⟂=qa⟂4πε0c2r

Z toho plyne:

E→radiacˇnıˊ(r→,t)=−qa→⟂(tr)4πε0c2r,tr=t−rc

Teď spadne z nebe: B→⟂E→,B=Ec,B→⟂n→,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→⟂=a→−a→∥=a→−(a→⋅r→0)r→0,r→0≔r→|r→| a→(t)=a0cos⁡(ωt−kr)z→,ω=ck Eradiacˇnıˊ(r→,t)=q4πε0c2⏟≕ℰ01rsin⁡(θ)cos⁡(ωt−kr)

Energetické veličiny v elektromagnetickém poli

Definice. Hustota energie w≔dEdV
Věta (viz ELMA). w=12(εE→2+B→2μ) (platí pro každé pole, nejen pro nějaké, jak bylo odvozeno na elmě)

Tok energie je

∂w∂t=12∂∂t(εEjEj+BjBjμ)=εE→⋅∂E→∂t+B→μ⋅∂B∂t=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:

ddt∫VwdV=−∮A≔∂VS→⋅dA→ dEVdt=−ΦS→
Definice. Hustota hybnosti g→≔dp→dV
Věta. Ve vakuu platí: g→=S→c2=ε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→=−αv→⏟trˇenıˊ+F→L Nechť je náboj silně tlumený, tedy ma→≪αv→. Potom v→≈F→Lα=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.) ⟨dp→dt⟩=⟨F→L⟩=⟨q(Exx→+vxByz→−vzByx→)⟩=q(0→+⟨vxBy⟩z→+0→) ⟨dEdt⟩=⟨F→L⋅v→⟩=⟨q(E→⋅v→+v→⋅(v→×B→))⟩=⟨q(E→⋅v→+B→⋅(v→×v→))⟩=⟨qExvx⟩ (pro maximalizaci zmatení je E tentokrát energie, nikoliv velikost E→) Z toho plyne: ⟨dp→dt⟩=1c⟨dEdt⟩z→ Matematiko, odpusť mi mé hříchy… ⟨dp→dV⟩=1c⟨dEdV⟩z→ ⟨g→⟩=⟨w⟩cz→ 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=v⋅B=Bεμ B→=r→×E→v w=12(εE→2+B→2μ)=εE2 S→=E→×B→μ=1μvE→×(r→×E→)=εμE2r→ g→=cwr→c2=wr→c g→=S→c2(???)
Definice. Intenzita I≔⟨|S→|⟩
Příklad. r→≔z→ E→≔E0x→cos⁡(ωt−kz) I=⟨vw⟩=⟨vεE2⟩=εμ⟨E02x→2cos⁡(ωt−kz)2⟩=12εμE02

Tlak záření

Už máme hybnost g→=S→c2. Taky máme tok energie S→, kde dE=S→⋅n→. 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=1361W⋅m−2, tedy p=4.5×10−6Pa. (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:

E→rad=−qa→⟂(tr)4πε0c2r tr=t−rc a→⟂=a→−(a→⋅r→0)r→0 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=S→⋅r→dA=SdA=ε0μ0E2 P=∫sfeˊradP=∫sfeˊraSdA=∫02π∫0πε0μ0E2r2sin⁡(θ)dθdφ=ε0μ0r2r2(qa4πε0c2)2∫0π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ý). x≔Acos⁡(ωt)z→ a→=−Aω2cos⁡(ωt)z→ P=μ0q26πcA2ω4cos⁡(ωtr)2=μ0q26πcA2ω4cos⁡(ωt−kr)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→=F→L=−e(E→+v→×B→)

Obecný postup výpočtu: model→x→→p→→P→→ε→εr→n

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

FB≤evB=evEc=βFE≪FE

Vezměme si konkrétně harmonickou vlnu:

E→(t)=E→0cos⁡(ω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)=−emE→0ω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: P→≔dp→dV. 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ě

n≈1+∑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=1−Ne2mε0ω2

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

1−Nc2mε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:

∂2E→∂t2=c2∇2E→∇→⋅E→=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→=∂Ex∂x=0=maˊ byˊt0 ∂2Ex∂t2=maˊ byˊtc2(∂2Ex∂y2+∂2Ex∂z2)

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

−ω2Ex=c2(Y′′(y)exp⁡(i(ωt−kz))−kz2Ex) −ω2Y=c2(Y′′−k2Y) Y′′+(ω2c2−kz2)⏟≕KY=0

Pokud K≤0, dokážeme to snadno vyřešit a zjistíme, že funguje jenom triviální řeséní. Nechť tedy K>0,ky≔K. 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(ωt−kzz))

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

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

(kde k≔kz, 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 k≕iϰ, 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,E→1(r→,t),B→1(r→,t). Vpravo máme ε2,μ2,E→2(r→,t),B→2(r→,t).

Samozřejmě použijeme Maxwellovy rovnice:

∇→⋅(εE→)=0 ∇→⋅B→=0 ∇→×E→=−∂B→∂t ∇→×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εE→⋅dS→=0 ∮SB→⋅dS→=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 h→0,Sp→0.

Φp1=ε1E→1⋅r→⋅Sp Φp2=ε2E→2⋅−r→⋅Sp Φpl=α⋅2πrh

Z toho nějak

ε1E→1⟂=ε2E→2⟂ B→1⟂=B→2⟂

(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→=ddt∫SεE→⋅dS→=ddt∫SεE→×n→dS=strˇednıˊhodnotaαlh

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

∮l=∫levaˊ+∫hornıˊ+∫pravaˊ+∫dolnıˊ→h→0∫levaˊ+∫pravaˊ=∫levaˊB→1μ1⋅t→dl−∫pravaˊB→2μ2⋅t→dl=maˊ byˊt0 B→1∥μ1=B→2∥μ2 E→1∥=E→2∥

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→=Ex0x→exp⁡(i(ωt−kz+φ1))+Ey0y→exp⁡(i(ωt−kz+φ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)=Ex0x→exp⁡(i(ωt+φ1))+Ey0y→exp⁡(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μ0⟨E2⟩=ε0μ0⟨Ex02x→2cos⁡(ωt+φ1)2+Ey02y→2cos⁡(ωt+φ1)2+2Ex0Ey0(x→⋅y→)cos⁡(ωt+φ1)cos⁡(ωt+φ2)⟩=ε0μ0⟨Ex02cos⁡(ω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. E→in=E→∥+E→⟂ E→out=E→∥=(E→in⋅n→)n→

Vsuvka: Jonesův kalkulus

E→(t)=Ex0x→exp⁡(i(ωt+φ1))+Ey0y→exp⁡(i(ωt+φ2))=(Ex0x→exp⁡(iφ1)+Ey0y→exp⁡(iφ2))⏟≕E→^exp⁡(iωt) E→^≔(Ex0exp⁡(iφ1)Ey0exp⁡(iφ2))
Jonesův vektor pro lineární polarizátor bude E→^out=𝐏n→E→^in (kde 𝐏 je projektor). E→out=(Einxnx+Einxny)(nxny)=(Einxnx2+EinynxnyEinxnxny+Einyny2)=(nx2nxnynxnyny2)⏟𝐏E→in

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ť n→1≔x→,n→2≔n→≔(cos⁡(θ)sin⁡(θ)).

U druhého polarizátoru máme:

E→in=E0x→cos⁡(ωt) Iin=⟨E→in2⟩=E022 E→out=(E→in⋅n→)n→=E0cos⁡(θ)n→cos⁡(ωt) Iout=⟨E→out2⟩=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:

E→out=E0xx→cos⁡(ωt+φ1−kd)+E0yy→cos⁡(ωt+φ2−kd)

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:

E→out=E0xx→cos⁡(ωt+φ1−k1d)+E0yy→cos⁡(ωt+φ2−k2d) δ′φ=(ωt+φ1−k1d)+(ωt+φ2−k2d)=φ1−φ2+(k1−k2)d≕δφ+Δφ Δφ=(k2−k1)d=ωc(n1−n2)d=2πλ0(n2−n1)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

E→out=E0xx→cos⁡(ωt+φ1−(k2−k1)d)+E0yy→cos⁡(ω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→^n→1+E→^n→2=𝐈E→^in=(𝐏n→1+𝐏n→2)E→^in=𝐏n→1E→^in+𝐏n→2E→^in E→^out=(exp⁡(iΔφ)𝐏n→1+𝐏n→2)E→^in 𝐃=exp⁡(iΔφ)𝐏n→1+𝐏n→2

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?

E→in=Exx→cos⁡(ωt)=E02(x→cos⁡(ωt)+y→sin⁡(ωt))+E02(x→cos⁡(ωt)−y→sin⁡(ωt)) E→out=Exx→cos⁡(ωt+Δφ)=E02(x→cos⁡(ωt+Δφ)+y→sin⁡(ωt+Δφ))+E02(x→cos⁡(ωt+Δφ)−y→sin⁡(ωt+Δφ))=⋯=(E0cos⁡(Δφ2)x→+E0sin⁡(Δφ2)y→)cos⁡(ωt+Δφ2)=E0n→cos⁡(ω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→=Ex0x→cos⁡(ωt+φ1)+Ey0y→cos⁡(ωt+φ2)

Měřicí přístroj má nějakou rozlišovací dobu troz⋙T

I=⟨E→2⟩troz

Definujme další čtyři intenzity:

Ix≔⟨Ex2⟩=⟨Ex02cos⁡(ωt+φ1)2⟩=Ex022 Iy≔⟨Ey2⟩=⟨Ey02cos⁡(ωt+φ1)2⟩=Ey022 Ixy≔⟨ExEy⟩=⟨Ex0Ey0cos⁡(ωt+φ1)cos⁡(ωt+φ2)⟩=⋯=Ex0Ey0cos⁡(φ1−φ2)2 Ixy≔⟨Ex(ω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 T≈10−14s, koherenční dobu tkoh≈10−8s a rozlišovací dobu troz. Pokud troz≪tkoh, jde o rychlý přístroj a můžeme prostě sledovat změny polarizace. Horší případ je troz≫tkoh — „pomalý přístroj“? Jelikož T≪tkoh, budou již střední hodnoty závislé na čase:

Ix(t)=⟨Ex2⟩T=Ex0(t)22 Iy(t)=⟨Ey2⟩T=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→|=(Ix−Iy)2+2Ixy2+2Ixy2(Ix+Iy)2=Ex04−2Ex02Ey02+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=0∨Ey0=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á E→d, odrazí se E→r, projde Et. Uvažujme harmonickou vlnu E→=E→d0exp⁡(i(ωt−k→d⋅r→)),k→d=kdn→d,ω=cn1kd.

Předpokládejme, že E→r=exp⁡(i(ωrt−k→r⋅r→)),k→r=krn→r,ωr=cn1kr a E→t=exp⁡(i(ωtt−k→t⋅t→)),k→t=ktn→t,ωt=cn2kt.

Natočíme osy tak, aby k→d=(0,kdy,kdz). Tím zajistíme, že E→d 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: k→r=(0,kry,krz),k→t=(0,kty,ktz). Na x tedy vůbec nezáleží — jde o rovinný problém.

Obrázek situace
k→d=(0,kdsin⁡(ϑd),kdcos⁡(ϑd)) k→r=(0,krsin⁡(ϑr),−krcos⁡(ϑr)) k→t=(0,ktsin⁡(ϑt),ktcos⁡(ϑt))

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

E→1∥(z=0)=E→2∥(z=0) E→d∥(z=0)+E→r∥(z=0)=E→t∥(z=0) E→d0∥exp⁡(i(ωt−kdsin⁡(ϑd)y))+E→r0∥exp⁡(i(ωt−krsin⁡(ϑr)y))=E→t0∥exp⁡(i(ωt−ktsin⁡(ϑ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 ϑc≔arcsin⁡(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|k→1|2=c2n22(kty2+ktz2)=c2n22(kdy2+ktz2)=c2n22(kd2sin⁡(ϑd)2+ktz2)=c2n22(ω2c2n12sin⁡(ϑ1)2+ktz2) −ϰ2≔ktz2=ω2c2(n22−n12sin⁡(ϑd)2) E→t=E→t0exp⁡(i(ωt−k→t⋅r→))=E→t0exp⁡(i(ωt−ktyy))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:

R∥≈n2cos⁡(ϑ1)−n1cos⁡(ϑ2)n2cos⁡(ϑ1)+n1cos⁡(ϑ2)R⟂≈n1cos⁡(ϑ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:

R∥≈tan⁡(ϑ2−ϑ1)tan⁡(ϑ2+ϑ1)R⟂≈sin⁡(ϑ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⟩=⟨(E→1+E→2)2⟩=⟨E12⟩+⟨E22⟩+2⟨E→1⋅E→2⟩⏟interferencˇ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⁡(ωt−kl1,2+φ1,2)

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

ID=⟨(E1D+E2D)2⟩=⟨E12cos⁡(⋯)2+E22cos⁡(⋯)2+2E1E2cos⁡(ωt−kl1+φ1)cos⁡(ωt−kl2+φ2)⟩=E122+E222+⟨E1E2(cos⁡(2ωt−k(l1+l2)+φ1+φ2)+cos⁡(k(l2−l1)+φ1−φ2))⟩=I1+I2+2I1I2cos⁡(k(l2−l1)+φ1−φ2)≕I1+I2+Ii

Fáze interference je tedy posunuta o Δφ=k(l2−l1)=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+I2−2I1I2,I1+I2+2I1I2]=[(I1−I2)2,(I1+I2)2]

Pokud by vysílané intenzity byly stejné, tedy I1=I2≕I0, 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⁡(ωt−kl1,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(l2−l1)+φ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 E→dop a ona bude „indukovat“ pole E→ind. Z definice neprůhlednosti musí být E→dop+E→ind=0.

Z

Přepážku si rozdělíme na části A,B. Zřejmě E→ind=E→indA+E→indB.

Z B  A 

Nyní dejme pryč B. Za přepážkou již bude nějaké pole E→A. Jelikož předtím tam nebylo nic, musí být E→A=−E→indB. 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ě E→indA. Ale pokud je otvor dostatečně velký (D≫λ), můžeme tyto efekty zanedbat.

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

E→A+E→B=−(E→indA+E→indb)(Babinetu˚v princip) E→A=E→dop−E→B

Difrakční integrál

Předpokládejme:

Každý malý úsek přepážky vyzařuje vlnu E→0rexp⁡(i(ωt−kr)). Celkově tedy máme

E→A=−E→indB=∫BE→0rexp⁡(i(ωt−kr))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+(X−x)2+(Y−y)2

Integrál tedy bude vypadat takto:

E→A(x,y)=∬SE→0r(x,y,X,Y)exp⁡(i(ωt−kr(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 R≔L2+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+(X−x)2+(Y−y)2−x2−y2=R2+X2+Y2−2xX−2yY

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

r=R1+X2+Y2−2xX−2yYR2≈R(1−xX+yYR2)

Ještě hrubší aproximace je r≈R nebo dokonce r≈L. 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.

E→A(x,y)≈∬SE→0Rexp⁡(i(ωt−k(R−xX+yYR)))dXdY=E→0Rexp⁡(i(ωt−kR))∬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+Y22R≪1

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

Δφ≤2πλD2412L=πD24λL≤D2λL≪maˊ 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 l≈r0−rX, 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→=E→0Rexp⁡(i(ωt−kR))(exp⁡(ikRxd2)+exp⁡(−ikRxd2))=2E→0Rexp⁡(i(ωt−kR))cos⁡(kdxR)=2E→0Rexp⁡(i(ωt−kR))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˚vvzorecE→0Rexp⁡(i(ωt−kR))∑j=0N−1exp⁡(ikRxjd)=geometrickaˊrˇadaE→0Rexp⁡(i(ωt−kR))exp⁡(ikRxdN)−1exp⁡(ikRxd)−1=vynaˊsobenıˊjednicˇkouE→0Rexp⁡(i(ωt−kR))exp⁡(ikRxdN)−1exp⁡(ikRxd)−1⋅exp⁡(−ikRxdN2)exp⁡(−ikRxdN2)⋅exp⁡(−ikRxd2)exp⁡(−ikRxd2)=sin⁡(x)=exp⁡(ix)−exp⁡(−ix)2E→0Rexp⁡(i(ωt−kR−kRxdN−12))sin⁡(kdNx2R)sin⁡(kdx2R) I=⟨E2⟩=E→022R2sin⁡(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

E→∝absin⁡(kax2R)kax2Rsin⁡(kby2R)kby2R=y=0absin⁡(kasin⁡(θ)2)kasin⁡(θ)2 I∝sin⁡(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⁡(nu−xsin⁡(u))du

Maxima teď budou u0=0,u1≈1.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 tkoh≈10−9s, tedy koherenční délka je lkoh=ctkoh≈30cm. 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.