Wie kann ich die Achsen der Polarisationsellipse aus dem Jones-Vektor des Lichts erhalten?

Ich analysiere den Polarisationszustand einer monochromatischen, kohärenten Lichtquelle, für die ich den Jones-Vektor der Polarisation kenne.

E = ( E. x E. y ) = ( | E. x | e i φ x | E. y | e i φ y ) , E. = ( E. x E. y ) = ( | E. x | e ich φ x | E. y | e ich φ y ) ,
und ich möchte es in Bezug auf eine Haupt- und eine Nebenachse der Elliptizität erweitern, dh in der Form
E = e i φ ( A u ^ + i B v ^ ) = e i φ ( A ( cos ( θ ) Sünde ( θ ) ) + i B ( - sin ( θ ) cos ( θ ) ) ) , E. = e ich φ ( EIN u ^ + ich B. v ^ ) = e ich φ ( EIN ( cos ( θ ) Sünde ( θ ) ) + ich B. ( - - Sünde ( θ ) cos ( θ ) ) ) ,
oder wie grafisch wie folgt dargestellt:

Bildquelle

Wikipedia bietet ein mehrstufiges Verfahren, das die Stokes-Parameter durchläuft, aber ich denke, es gibt sicherlich einen saubereren und direkteren Weg, dies zu erreichen EIN EIN , B. B. , u ^ u ^ , v ^ v ^ , θ θ und die Komponenten A u ^ EIN u ^ und B v ^ B. v ^ , von E. x E. x und E. y E. y und es ist nicht besonders offensichtlich aus den Suchergebnissen, die ich finden kann. Was ist der sauberste Weg, dies zu tun?


Um es klar auszudrücken: Was meiner Meinung nach in den vorhandenen Ressourcen fehlt und was die Frage direkt verlangt, ist eine explizite Reihe von möglichst einfachen Verbindungen für die genannten Parameter (alle EIN EIN , B. B. , u ^ u ^ , v ^ v ^ , θ θ und die Komponenten A u ^ EIN u ^ und B v ^ B. v ^ ) in Bezug auf die kartesischen Komponenten E. x E. x und E. y E. y . Schemata, die einfach an andere komplexe Manipulationen gesendet werden, sind bereits bei Wikipedia erhältlich und entsprechen nicht den Anforderungen der Frage.

Um ganz klar zu sein, dieses Diagramm, das Sie aus Wikipedia streichen, scheint so etwas wie die Punkte zu sein
( p ( ψ ) , q ( ψ ) ) = ( | E. x | cos ( ψ + ϕ x ) ,   | E. y | cos ( ψ + ϕ y ) . ( p ( ψ ) , q ( ψ ) ) = ( | E. x | cos ( ψ + ϕ x ) , | E. y | cos ( ψ + ϕ y ) .
Was den unternehmungslustigen Schüler davon abhält, beispielsweise zu minimieren / zu maximieren p 2 + q 2 p 2 + q 2 etwas zu finden wie (rektale Extraktionsmethode)
bräunen ( 2 ψ c ) = - ( | E. x | 2 Sünde ( 2 ϕ x ) + | E. y | 2 Sünde ( 2 ϕ y ) ) / ( | E. x | 2 cos ( 2 ϕ x ) + | E. y | 2 cos ( 2 ϕ y ) ) bräunen ( 2 ψ c ) = - - ( | E. x | 2 Sünde ( 2 ϕ x ) + | E. y | 2 Sünde ( 2 ϕ y ) ) /. ( | E. x | 2 cos ( 2 ϕ x ) + | E. y | 2 cos ( 2 ϕ y ) )
und dann geben Sie eine Reihe von Antworten in Bezug auf ψ c ψ c ? Nur dass es nicht elegant genug ist, um diesen unordentlichen Arkustangens in der Mitte zu haben?
Um das Feld im Koordinatensystem der Ellipse auszudrücken, müssen Sie die Eigenachsen dieser Ellipse aus dem Jones-Polarisationsvektor ermitteln. Der einzige Weg, dies zu tun, den ich kenne, ist die Diagonalisierung von 2 × 2 2 × 2 Matrix. Die dabei zu diagonalisierende Matrix ist die Dichtematrix des Lichts. Die Zerlegung dieser Matrix in die Pauli-Matrizen gibt Ihnen die Stokes-Parameter. Das ist also nicht die Frage nach der Physik, sondern die Frage, wie man die Berechnung vereinfacht, nicht wahr?
Vielleicht könnten Sie das sechste Kapitel [in diesem Buch] [1] lesen. [1]: optics.byu.edu/BYUOpticsBook_2015.pdf

Antworten (3)

Lassen Sie mich ein zweites Mal versuchen. Ich verwende https://math.stackexchange.com/questions/1204131/converting-a-rotated-ellipse-in-parametric-form-to-cartesian-form als Ressource.

Vorherige Manipulationen

Das physikalische elektrische Feld ist

E. p h y s = R e [ E e i ω t ] = Re [ ( | E. x | e i φ x | E. y | e i φ y ) e i ω t ] = ( | E. x | cos ( ω t + φ x ) | E. y | cos ( ω t + φ y ) ) = ( | E. x | [ cos ( ω t ) cos ( φ x ) - Sünde ( ω t ) sin ( φ x ) ] | E. y | [ cos ( ω t ) cos ( φ y ) - Sünde ( ω t ) sin ( φ y ) ] ) E. p h y s = R. e [ E. e ich ω t ]] = R. e [ ( | E. x | e ich φ x | E. y | e ich φ y ) e ich ω t ]] = ( | E. x | cos ( ω t + φ x ) | E. y | cos ( ω t + φ y ) ) = ( | E. x | [ cos ( ω t ) cos ( φ x ) - - Sünde ( ω t ) Sünde ( φ x ) ]] | E. y | [ cos ( ω t ) cos ( φ y ) - - Sünde ( ω t ) Sünde ( φ y ) ]] )

Dies ist eine parametrische Gleichung für eine Ellipse, die vom elektrischen Feld verfolgt wird.

Hauptachsenwinkel

Lassen

a = | E. x | cos ( φ x ) ein = | E. x | cos ( φ x )

b = | E. x | Sünde ( φ x ) b = | E. x | Sünde ( φ x )

c = | E. y | cos ( φ y ) c = | E. y | cos ( φ y )

d = | E. y | Sünde ( φ y ) d = | E. y | Sünde ( φ y )

Im Vergleich mit der verknüpften math.SE-Frage gibt die akzeptierte Antwort an, dass die Haupt- und Nebenachse auf der Ellipse (die auf den Ursprung zentriert ist) zeigen

ω t = 1 2 Arctan 2 ( a b + c d ) ( a 2 + c 2 ) - ( b 2 + d 2 ) + k π 2 ( 0 k 3 )   . ω t = 1 2 Arctan 2 ( ein b + c d ) ( ein 2 + c 2 ) - - ( b 2 + d 2 ) + k π 2 ( 0 k 3 ) .

Vom OP geforderte Erweiterung

In der Tat ist diese Größe der Winkel θ θ in der Erweiterung in der OP-Frage erscheinen. Allerdings abhängig von welchem ​​Wert für k k wird in Mai gewählt π / 2 - θ π /. 2 - - θ Eine Unterscheidung der Fälle in Abhängigkeit vom Sektor des Arktans ist erforderlich.

Dies ergibt natürlich auch u ^ u ^ und v ^ v ^ Die Erweiterung kann nun leicht erhalten werden, indem der Jones-Vektor auf diese Basis projiziert wird.

Die Formeln können mitgehen, stellen jedoch eine enge Formlösung für das Problem dar, bis hin zur Fallunterscheidung der Auswahl k k . Ich sehe nicht ein, wie eine einfachere Lösung existieren kann, da die Formel für den angegebenen Winkel nicht algebraisch reduzierbar erscheint.

Entschuldigung, dass ich das Kopfgeld nicht vergeben habe, aber das war ziemlich kurz von dem, wonach ich gesucht habe. Ich habe jetzt die sauberen Ausdrücke gefunden, nach denen ich gesucht habe, und sie als Antwort hinzugefügt.
@EmilioPisanty keine Sorge, ich verstehe jetzt was du meinst (siehe Kommentar unter deiner Antwort).

Der sauberste Weg, dies zu tun, wird von Michael Berry in der Zeitung angeboten

Indexformeln für singuläre Polarisationslinien. MV Berry, J. Opt. A: Reine Appl. Opt. 6 , 675–678 (2004) , Autor eprint .

In Berrys Notation kann das elektrische Feld wie folgt geschrieben werden

E = P + i Q = e i γ ( A + i B ) , E. = P. + ich Q. = e ich γ ( EIN + ich B. ) ,
wo P. P. , Q. Q. , EIN EIN und B. B. sind reelle Vektoren, EIN EIN und B. B. sind jeweils die Haupt- und Nebenachse der Polarisationsellipse, und diese beiden werden bis zu einem Vorzeichen durch definiert A B = 0 EIN B. = 0 und | A | | B | | EIN | | B. | . Mit dieser Notation werden die Polarisationsachsen und die Phase definiert als
γ = 1 2 arg ( E E ) und A + i B = E. E. - - - - - - - - - - - - E. E. - - - - - - - - - - - - E , γ = 1 2 arg ( E. E. ) und EIN + ich B. = E. E. | E. E. | E. ,
oder mit anderen Worten
A = 1 E. E. - - - - - - - - - - - - R e [ E. E. - - - - - - - - - - - - E ] und B = 1 E. E. - - - - - - - - - - - - Ich bin [ E. E. - - - - - - - - - - - - E ] . EIN = 1 | E. E. | R. e [ E. E. E. ]] und B. = 1 | E. E. | ich m [ E. E. E. ]] .
Es gibt eine offensichtliche Mehrdeutigkeit der Zeichen, wenn beide umgedreht werden EIN EIN und B. B. und hinzufügen π π zu γ γ ändert nichts (dh das Drehen der Polarisationsellipse um 180 ° entspricht dem Hinzufügen einer Phase), was sich in den Verzweigungsschnitten sowohl der Argument- als auch der Quadratwurzelfunktion widerspiegelt. Diese greifen natürlich ineinander, solange beide Astschnitte auf demselben Schnitt ausgeführt werden, idealerweise entlang der negativen realen Achse.

Als weiterer kniffliger Punkt sollte man beachten, dass diese Formeln nicht definiert sind, wenn E E = 0 E. E. = 0 , was einer zirkularen Polarisation entspricht; in diesem Fall sowohl die Polarisationsachsen als auch die Phase γ γ an der Hauptachse sind schlecht definiert, so dass dies kein Problem ist.

Als zusätzlichen Bonus gibt dieser Ansatz natürlich auch die Richtung der Normalen zur Ebene der Polarisationsellipse in der Form an

C = 1 2 Ich bin ( E. × E ) = P × Q = A × B , C. = 1 2 ich m ( E. × E. ) = P. × Q. = EIN × B. ,
wo das Kreuzprodukt E. × E. E. × E. ist natürlich imaginär, da sein Konjugat minus sich selbst ist. Dies wird natürlich verschwinden, wenn E. E. und E. E. (oder P. P. und Q. Q. ) sind linear abhängig, was einer linearen Polarisation entspricht; in diesem Fall, B. B. wird verschwinden, weil E. E. - - - - - - - - - - - - E. E. E. E. ist natürlich echt.

Berry Credits

Polarisationssingularitäten in paraxialen Vektorfeldern: Morphologie und Statistik. MR Dennis, Opt. Kommun. 213 , 201–21 (2002) , eprint .

für diese Form und diese Referenz enthält einen umfassenderen Beweis dafür, wie und warum die Zerlegung funktioniert.

Dies ist in der Tat ziemlich einfach, sobald Sie erkennen, dass die Zerlegung als E = e i γ ( A + i B ) E. = e ich γ ( EIN + ich B. ) muss wie oben vorhanden sein, da sich unter diesen Bedingungen das Punktprodukt auf reduziert

E E = e 2 i γ ( A + i B ) ( A + i B ) = e 2 i γ ( A. 2 - B. 2 ) , E. E. = e 2 ich γ ( EIN + ich B. ) ( EIN + ich B. ) = e 2 ich γ ( EIN 2 - - B. 2 ) ,
wo EIN 2 - B. 2 EIN 2 - - B. 2 ist real und positiv, daher ergibt die Argumentation beider Seiten natürlich die Phase als 2 γ = arg ( E E ) . 2 γ = arg ( E. E. ) . In ähnlicher Weise kehrt die Annahme des Moduls dieser Gleichung zurück EIN 2 - B. 2 = | E E | EIN 2 - - B. 2 = | E. E. | , so können wir einfach den Phasenfaktor als erhalten
e 2 i γ = E E. EIN 2 - B. 2 = E E. | E E | ,   damit   e i γ = E E. - - - - - - - - | E E. - - - - - - - - | ,   und deshalb   e - i γ = E. E. - - - - - - - - - - - - | E. E. - - - - - - - - - - - - | ;; e 2 ich γ = E. E. EIN 2 - - B. 2 = E. E. | E. E. | , damit e ich γ = E. E. | E. E. | , und deshalb e - - ich γ = E. E. | E. E. | ;;
die Charakterisierung für A + i B. EIN + ich B. dann folgt aus E = e i γ ( A + i B ) E. = e ich γ ( EIN + ich B. ) durch einfaches Teilen durch e i γ e ich γ .

Das ist fair, ich kann sehen, wie ordentlich das ist. Anfangs schien es mir nur ein geometrisches Problem zu sein, weshalb ich mit einer ziemlich kurzen Antwort antwortete. Ich denke auch, dass Ihre Lösung hier durch den intelligenten Einsatz von Parametrisierung und komplexen Operationen mehr physische Einblicke bietet. Deshalb +1 von mir

Ich denke, Sie können dies mit einer Singularwertzerlegung erreichen. Ich werde mit dem Schreiben beginnen E. E. in der folgenden Form

E = E. x x ^ + E. y y ^ = E. u u ^ + E. v v ^ E. = E. x x ^ + E. y y ^ = E. u u ^ + E. v v ^

wo E. x E. x und E. y E. y sind komplexe Zahlen bekannt. Im Allgemeinen E. x E. y 0 E. x E. y 0 , aber wir wollen E. u E. v = 0 E. u E. v = 0 . Wir können dies in Matrixform schreiben als

( x ^ y ^ ) ( E. x E. y ) = ( u ^ v ^ ) ( E. u E. v ) ( x ^ y ^ ) ( E. x E. y ) = ( u ^ v ^ ) ( E. u E. v )

Wir sollten schreiben können

( x ^ y ^ ) = ( u ^ v ^ ) Λ T. ( x ^ y ^ ) = ( u ^ v ^ ) Λ T.

wo Λ Λ ist eine reale orthogonale Matrix, die eine Drehung von Koordinatensystemen beschreibt. Es ist derzeit unbekannt. Da nehmen wir das an x ^ x ^ und y ^ y ^ sind also orthogonal u ^ u ^ und v ^ v ^ wird auch orthogonal sein. Und wir können auch schreiben

( E. x E. y ) = A ( 1 ich ) ( E. x E. y ) = EIN ( 1 ich )

wo A = ( a c b d ) EIN = ( ein b c d ) ist eine bekannte reale Matrix. Wir können uns zersetzen EIN EIN unter Verwendung einer Singularwertzerlegung

A = U. S. V. T. EIN = U. S. V. T.

wo U. U. , S. S. , und V. V. sind echte Matrizen. U. U. und V. V. sind orthogonale Matrizen und S. = ( A. 0 0 B. ) S. = ( EIN 0 0 B. ) wo EIN EIN und B. B. sind positiv mit A B. EIN B. .

Das alles ergibt

Λ T. U. S. V. T. ( 1 ich ) = ( E. u E. v ) Λ T. U. S. V. T. ( 1 ich ) = ( E. u E. v )

Wenn wir uns entscheiden Λ = U. Λ = U. , was gibt u ^ u ^ und v ^ v ^ , dann reduziert sich dies auf

S. V. T. ( 1 ich ) = ( E. u E. v ) S. V. T. ( 1 ich ) = ( E. u E. v )

Schon seit V. V. ist eine orthogonale Matrix, so ist es auch V. T. V. T. , und S. S. ist diagonal und real, also E. u E. v = 0 E. u E. v = 0 .

Ich habe das in Eile aufgeschrieben, also habe ich vielleicht etwas verpasst oder war ungenau. Ich entschuldige mich, wenn das der Fall ist.