Hamilton-optikk

Hamilton-optikk er en matematisk metode for beskrivelse av lysstråler i geometrisk optikk. Den er basert på Fermats prinsipp og ble utviklet i løpet av noen få år rundt 1830 av William Rowan Hamilton. Han etablerte forskjellige, karakteristiske funksjoner som oppfyller ulike differensialligninger avhengig av det medium lyset går gjennom. Ved å angi utgangspunkt og retning av en lysstråle kan man herav beregne hvordan strålen videre vil bevege seg.

Metoden fikk stor oppmerksomhet i 1832 da Hamilton kunne forutsi et spesielt, optisk fenomen i dobbeltbrytende krystaller og som kort tid etterpå ble eksperimentelt verifisert. På slutten av 1800-tallet ble hans beskrivelse av optiske fenomen vist å være i overensstemmelse med Maxwell-teorien for elektromagnetiske bølger når disse forenkles i eikonalapproksimasjonen.

Lignende betraktninger førte til at Hamilton i de følgende årene kunne beskrive bevegelse av materielle partikler på en mer systematisk måte som i dag kalles for Hamilton-mekanikk. Denne fikk større betydning enn hans optiske arbeid og ville senere danne grunnlaget for moderne kvantemekanikk.

Karakteristisk funksjon

Ved bruk av Fermats prinsipp kan man beregne hvordan en lysstråle beveger seg fra et punkt r₁ = Mal:Nowrap gjennom et medium med brytningsindeks n(r) til et annet punkt Mal:Nowrap ved å løse differensialligningen som følger fra dette variasjonsprinsippet. En vilkårlig lysstråle kan beskrives som en kurve r = r(s) når den parametrisere ved buelengden s. Matematisk kan dette prinsippet da formuleres som at den optiske lengden gitt ved integralet

V (𝐫_{1}, 𝐫_{2}) = \int_{𝐫_{1}}^{𝐫_{2}} d s n (𝐫)

langs strålen, skal ha en ekstremalverdi som vanligvis er et minimum. Dette uttrykket er Hamiltons første, karakteristiske funksjon.^[1]

En generell variasjon r → r + δr av lysets gang gjennom mediet resulterer i en tilsvarende variasjon i den optiske veilengden,

δ V = \int_{𝐫_{1}}^{𝐫_{2}} (δ d s) n (𝐫) + \int_{𝐫_{1}}^{𝐫_{2}} d s δ n (𝐫)

hvor δn = ∇n⋅δr. Ved å benytte at Mal:Nowrap slik at Mal:Nowrap, følger at

δ d s = \frac{d 𝐫}{d s} \cdot δ d 𝐫 = 𝐬 \cdot d δ 𝐫

der s = d r/ds er tangentvektoren til lysstrålen med Mal:Nowrap. Vektoren Mal:Nowrap oppfyller derfor Mal:Nowrap og blir ofte kalt for den optiske impulsen til strålen.^[2]

Variasjonen δV kan forenkles ved en partiell integrasjon. Når man lar også endepunktene r₁ og r₂ være frie, finner man dermed

δ V = 𝐩_{2} \cdot δ 𝐫_{2} - 𝐩_{1} \cdot δ 𝐫_{2} + \int_{𝐫_{1}}^{𝐫_{2}} d s (\nabla n - \frac{d 𝐩}{d s}) \cdot δ 𝐫

For en lysstråle mellom to faste punkt som betyr at δr₁ = δr₂ = 0, vil nå Fermats prinsipp Mal:Nowrap bety at integralet må være lik med null. For at dette kravet skal være oppfylt for alle variasjoner δr, må

\frac{d 𝐩}{d s} = \nabla n

være oppfylt. Dette er Euler-Lagrange-ligningen for kurven som bestemmer banen til den virkelige lysstrålen. Den viser hvordan retningen til lysstrålen er gitt ved forandringer i brytningsindeksen.

Eikonalligningen

Det generelle uttrykket for variasjonen av den karakteristiske funksjonen inneholder mer informasjon. Hvis man betrakter virkelige lysstråler mellom to variable endepunkt r₁ og r₂ slik at integralet ikke bidrar til δV, vil derfor ha den differensielle sammenhengen

d V (𝐫_{1}, 𝐫_{2}) = 𝐩_{2} \cdot d 𝐫_{2} - 𝐩_{1} \cdot d 𝐫_{1}

Retningene p₁ og p₂ til lysstrålene i endepunktene er derfor gitt ved gradienten til den karakteristiske funksjonen,

𝐩_{1} = - \frac{\partial}{\partial 𝐫_{1}} V (𝐫_{1}, 𝐫_{2}), 𝐩_{2} = \frac{\partial}{\partial 𝐫_{2}} V (𝐫_{1}, 𝐫_{2})

De to lysretningene står derfor vinkelrett på flatene definert ved V(r₁,r₂) = konstant. De opptrer derfor her på tilsvarende måte som ekvipotensialflater benyttes for å beskrive vektorfelt. Betydningen av slik fiktive flater i optikken ble først innsett av Malus på begynnelsen av 1800-tallet.^[3]

Da man generelt har at den optiske impulsvektoren oppfyller Mal:Nowrap, vil den karakteriskte funksjonen til Hamilton være en løsning av differensialligningen

(\nabla V)^{2} = (\frac{\partial V}{\partial x})^{2} + (\frac{\partial V}{\partial y})^{2} + (\frac{\partial V}{\partial z})^{2} = n^{2}

Dette er akkurat samme ligning som opptrer for eikonalfunksjonen S(r) som opptrer når man beskriver lysstrålene i bølgeteori. Den har et differensial Mal:Nowrap slik at man formelt kan skrive

V (𝐫_{1}, 𝐫_{2}) = S (𝐫_{2}) - S (𝐫_{1})

Dermed er disse to beskrivelsene knyttet sammen.^[4]

Eksempel

Den karakteristiske funksjonen kan formelt finnes ved å løse eikonalligningen. For de fleste situasjoner er det vanskelig eller umulig å finne eksakte løsninger. Ett av de enkleste tilfeller hvor dette er mulig, er refleksjon i et plant speil. For to punkt Mal:Nowrap og Mal:Nowrap som ligger i xy-planet, er den optiske veilengden for en lysstråle som går langs en rett linje fra det ene punktet til det andre gitt som

V_{0} (𝐫_{1}, 𝐫_{2}) = \sqrt{(x_{1} - x_{2})^{2} + (y_{1} - y_{2})^{2}}

når brytningsindeksen er konstant n = 1. Dette er også den karakteristiske funksjonen i dette trivielle tilfellet.

Hvis nå et speil ligger i planet y = 0, vil lysstrålen deles i to linjer ved at den reflekteres i et punkt (x, 0). Den fulle veilengden blir dermed

V (𝐫_{1}, 𝐫_{2}) = \sqrt{(x - x_{1})^{2} + y_{1}^{2}} + \sqrt{(x_{2} - x)^{2} + y_{2}^{2}}

Fermats prinsipp krever nå at denne er et minimum, det vil si at posisjonen x er bestemt ved at den deriverte Mal:Nowrap. Det gir at

x = \frac{x_{1} y_{2} + x_{2} y_{1}}{y_{1} + y_{2}}

slik at x - x₁ = k y₁ og x₂ - x = k y₂ hvor k = (x₂ - x₁)/(y₁ + y₂). Innsatt gir dette for funksjonen

\begin{matrix} V (𝐫_{1}, 𝐫_{2}) & = y_{1} \sqrt{1 + k^{2}} + y_{2} \sqrt{1 + k^{2}} \\ = \sqrt{(x_{1} - x_{2})^{2} + (y_{1} + y_{2})^{2}} \end{matrix}

Gradientene ∇₁ og ∇₂ vil nå gi de trigonometriske funksjonene til vinkelene som lysstrålen danner med aksene i de to endepunktene. Bortsett fra et fortegn, er disse vinklene like store som de skal være ifølge refleksjonsloven.

Andre karakteristiske funksjoner

Da retningen p₁ til lyset sendt ut fra r₁ er gitt ved gradienten av den karakteristiske funksjonen i dette punktet, kan man ved en Legendre-transformasjon definere en ny, karakteristisk funksjon

W (𝐩_{1}, 𝐫_{2}) = V (𝐫_{1}, 𝐫_{2}) - 𝐩_{2} \cdot 𝐫_{2}

Ved å benytte differensialet dV(r₁,r₂) = p₂⋅d r₂ - p₁⋅d r₁, blir nå

d W (𝐫_{1}, 𝐩_{2}) = - 𝐩_{1} \cdot d 𝐫_{1} - 𝐫_{2} \cdot d 𝐩_{2}

Dermed er

𝐩_{1} = - \frac{\partial}{\partial 𝐫_{1}} W (𝐫_{1}, 𝐩_{2}), 𝐫_{2} = - \frac{\partial}{\partial 𝐩_{2}} W (𝐫_{1}, 𝐩_{2})

slik at den muliggjør å bestemme retningen til strålen i begynnelsespunktet og dens endepunkt da den selv er en funksjon av begynnelsespunktet og retningen i sluttpunktet. Mens den første funksjonen V(r₁,r₂) kalles for en punkt-til-punkt-funksjon, er denne andre, karakteristiske funksjonen er punkt-til-vinkel-funksjon.^[4]

Den tredje og siste, karakteristiske funksjonen er vinkel-til-vinkel-funksjonen gitt ved enda en Legendre-transformasjon,

T (𝐩_{1}, 𝐩_{2}) = W (𝐩_{1}, 𝐫_{2}) + 𝐩_{1} \cdot 𝐫_{1}

og dermed

d T (𝐩_{1}, 𝐩_{2}) = 𝐫_{1} \cdot d 𝐩_{1} - 𝐫_{2} \cdot d 𝐩_{2}

Nå er derfor

𝐫_{1} = \frac{\partial}{\partial 𝐩_{1}} T (𝐩_{1}, 𝐩_{2}), 𝐫_{2} = - \frac{\partial}{\partial 𝐩_{2}} T (𝐩_{1}, 𝐩_{2})

slik at denne funksjonen gjør det mulig å bestemme begynnelsespunktet og sluttpunktet ut fra retningene til strålen i disse to punktene.

For refleksjon i et sfærisk speil kan denne funksjonen beregnes. Angir man retningen til strålen med de tre koordinatene p = (α,β,γ) med brytningsindeks n = 1, er

T (𝐩_{1}, 𝐩_{2}) = R \sqrt{(α_{1} - α_{2})^{2} + (β_{1} - β_{2})^{2} + (γ_{1} - γ_{2})^{2}}

der R er redius til speilet. Uttrykket kan lett generaliseres til også å gjelde for brytning i en sfærisk grenseflate.^[4]

Lagrange-formalisme

Den optiske impulsen p = n d r/ds er proporsjonal med cosinus til de tre vinklene som angir retningen til lysstrålen. Disse er ikke uavhengige av hverandre, noe som medfører at denne impulsen har en lengde som er fiksert av brytningsindeksen n. Som ved benyttelsen av Fermats prinsipp kan man unngå dette problemet ved å benytte en annen parametrisering av den kurven som lysstrålen beskriver.

Av stor praktisk og teoretisk betydning er å velge en av de kartesiske koordinatene i posisjonsvektoren r = (x,y,z). Hvis lysets gang er hovedsakelig langs z-aksen, er det da naturlig å velge z som uavhengig parameter. Veien til lysstrålen er da gitt ved funksjonene Mal:Nowrap og Mal:Nowrap. Den optiske veilengden mellom to punkt eller eikonalfunksjonen er da gitt ved integralet

S (𝐫) = \int_{𝐫_{0}}^{𝐫} d s n (𝐫) = \int_{𝐫_{0}}^{𝐫} d z n (x, y, z) \sqrt{1 + {\dot{x}}^{2} + {\dot{y}}^{2}}

hvor $\dot{x} = d x / d z$ og $\dot{y} = d y / d z$ . Matematisk er nå dette problemet det samme som Hamiltons virkningsprinsipp for en partikkel som beveger seg i et todimensjonalt rom med koordinater q = (x,y) hvor z spiller rollen som en tidskoordinat.

Historisk var det de optiske arbeidene som førte Hamilton til en lignende beskrivelse av bevegelsen til materielle partikler. Hans Hamilton-mekanikk viste seg på sikt å få større betydning enn hans matematiske formulering av optiske fenomen.^[1]

Hamilton-funksjon

Fra virkningsintegralet kan man nå identifisere den optiske Lagrange-funksjonen med

L (q, \dot{q}, z) = n (q, z) \sqrt{1 + {\dot{q}}^{2}}

Den kanoniske impulsen

p = \frac{\partial L}{\partial \dot{q}} = \frac{n \dot{q}}{\sqrt{1 + {\dot{q}}^{2}}}

har to komponenter (p_x, p_y) som angir retningen til lysstrålen i forhold til z-aksen. Herfra finner man at $\dot{q} = n / (1 + {\dot{q}}^{2})^{1 / 2}$ slik at den tilsvarende Hamilton-funksjonen følger fra en Legendre-transformasjonen som

H = \dot{q} p - L = - \sqrt{n^{2} - p^{2}} = - p_{z}

Lysets gang beskrives nå av de to Hamilton-ligningene som begge er av første grad,

\begin{matrix} \dot{q} & = \frac{\partial H}{\partial p} = \frac{p}{\sqrt{n^{2} - p^{2}}} \\ \dot{p} & = - \frac{\partial H}{\partial q} = \frac{n}{\sqrt{n^{2} - p^{2}}} \frac{\partial n}{\partial q} \end{matrix}

De er ekvivalent med Euler-Lagrange-ligningen som følger direkte fra Fermats prinsipp og er av andre grad.^[2]

Paraksial approksimasjon

Hvis lysstrålen beveger seg hele veien tett langs z-aksen, kan bevegelsesligningene forenkles. Da er de optiske impulsene eller retningsvinklene små slik at Mal:Nowrap slik at Hamilton-funksjon tilnærmet blir

H = - n + \frac{p^{2}}{2 n}

som for en ikke-relativistisk partikkel. Det betyr at med samme nøyaktighet er $\dot{q} = p / n$ slik at den andre Hamilton-ligningen git

\frac{d}{d z} (n \dot{q}) = \frac{\partial n}{\partial q}

Da koordinaten q = (x,y), er dette ekvivalent med de samme paraksiale ligninger som også kan finnes direkte fra Fermats prinsipp i denne approksimasjonen.

Referanser

↑ ^1,0 ^1,1 M. Born and E. Wolf, Principles of Optics, Pergamon Press, Oxford (1980). ISBN 0-08-026482-4.
↑ ^2,0 ^2,1 H.A. Buchdahl, An Introduction to Hamiltonian Optics, Dover Publications, New York (1970). ISBN 0-486-67597-1.
↑ O. Darrigol, A History of Optics, Oxford University Press, Oxford (2012). ISBN 978-0-19-876695-7.
↑ ^4,0 ^4,1 ^4,2 J.L. Synge, Geometrical Optics: An Introduction to Hamilton's Method, Cambridge University Press (1962). ISBN 978-0-521-06590-0.

[BW-1] 1,0 ^1,1 M. Born and E. Wolf, Principles of Optics, Pergamon Press, Oxford (1980). ISBN 0-08-026482-4.

[Buchdahl-2] 2,0 ^2,1 H.A. Buchdahl, An Introduction to Hamiltonian Optics, Dover Publications, New York (1970). ISBN 0-486-67597-1.

[Darrigol-3] O. Darrigol, A History of Optics, Oxford University Press, Oxford (2012). ISBN 978-0-19-876695-7.

[Synge-4] 4,0 ^4,1 ^4,2 J.L. Synge, Geometrical Optics: An Introduction to Hamilton's Method, Cambridge University Press (1962). ISBN 978-0-521-06590-0.

[1]

[2]

[3]

[4]

Hamilton-optikk

Innhold

Karakteristisk funksjon

Eikonalligningen

Eksempel

Andre karakteristiske funksjoner

Lagrange-formalisme

Hamilton-funksjon

Paraksial approksimasjon

Referanser

Navigasjonsmeny

Hamilton-optikk

Karakteristisk funksjon

Eikonalligningen

Eksempel

Andre karakteristiske funksjoner

Lagrange-formalisme

Hamilton-funksjon

Paraksial approksimasjon

Referanser

Navigasjonsmeny

Søk