Finden von Vier-Wellen-Mischtermen in der Photon-Atom-Wechselwirkung (durch Hinzufügen von for-Feldern zu EIT)

Das Problem

Ich interessiere mich für die Lösung der Vier-Wellen-Mischbegriffe in elektromagnetisch induzierter Transparenz (EIT), wenn zusätzliche Felder hinzugefügt werden (unter Beibehaltung von 3 Atomebenen).

Wenn wir mit einem N-Ebenen-System arbeiten, verwenden wir normalerweise (N-1) Felder (EIT hat 3 atomare Ebenen und 2 Felder). Es ist bekannt, dass das Hinzufügen zusätzlicher Felder "Wellenmischungs"-Ausdrücke verursacht, bei denen Licht einer Frequenz in eine andere Frequenz umgewandelt wird.

Ich interessiere mich für das folgende Ebenenschema, wie im Diagramm dargestellt:

Geben Sie hier die Bildbeschreibung ein

Das Pegelschema hier ist identisch mit dem normalen EIT, außer dass wir ein zusätzliches EIT-System haben, das über das normale EIT verstimmt ist.

Wir erwarten, dass wenn eines der schwachen Felder ( Ω S , Ω P ) ausgeschaltet ist, wird Licht des ersten Halbbildes in der Frequenz des zweiten Halbbildes erzeugt. Um dies visuell zu sehen, erwarten wir, dass wann Ω P , Ω C , Ω D eingeschaltet sind, wird Licht in der Frequenz erzeugt ω S , und wir erhalten eine Vier-Wellen-Mischung (FWM), wie in der folgenden Abbildung angegeben:

Geben Sie hier die Bildbeschreibung ein

Ich erwarte Bedingungen des Formulars zu sehen E Ö u T e ich ω P T Ω S Ω C Ω D Und E Ö u T e ich ω S T Ω P Ω C Ω D

Und eine ähnliche Lösung (vierstufiges statt dreistufiges System) wurde hier gelöst . In ihrem Vier-Ebenen-System haben sie Begriffe ρ 31 = Ω P | Ω D | 2 Ω S Ω C Ω D D Und ρ 41 = Ω S | Ω C | 2 Ω S Ω C Ω D D .

Mein Versuch

Ich kann eine analytische Lösung für bekommen P 12 Und P 13 , aber ich sehe nicht die erwarteten Vier-Wellen-Mischbegriffe. Ich habe unten eine Zusammenfassung geschrieben, wie ich das erhalten habe (in der Hoffnung, dass mir jemand helfen kann, herauszufinden, was ich falsch mache).

Der Hamiltonian, mit dem ich beginne, ist:

H E ICH T = [ ω A 0 μ 13 0 ω B μ 23 μ 13 μ 23 ω C ] ( E P C Ö S ( ω P T + ϕ P ) + E S C Ö S ( ω S T + ϕ S ) + E C C Ö S ( ω C + ϕ C ) + E D C Ö S ( ω D + ϕ D ) )

Die Anwendung der RWA-Approximation erscheint ziemlich einfach. Es scheint nicht so, als ob diese zusätzlichen Felder irgendetwas an der RWA ändern würden, und wir können die gegenläufigen Begriffe wie gewohnt entfernen. E S Und E P sind ähnlich, also die gleichen Regeln für ω P bewerben für ω S und ebenso für E C Und E D : ( ω S ω A C ) << ( ω C ω A C ) )

Durch Entfernen der gegenläufigen Terme kann dies reduziert werden auf

( ω 1 ω 3 0 e ich T Δ P Ω P + e ich T ( Δ P Δ 1 ) Ω S 0 ω 2 ω 3 e ich T Δ C Ω C + e ich T ( Δ C Δ 2 ) Ω 2 e ich T Δ P Ω P + e ich T ( Δ P Δ 1 ) Ω S e ich T Δ C Ω C + e ich T ( Δ C Δ 2 ) Ω D 0 )
Oder im ursprünglichen Bezugssystem:

( ω 1 ω 3 0 e ich T ω P Ω P + e ich T ( ω P Δ 1 ) Ω S 0 ω 2 ω 3 e ich T ω C Ω C + e ich T ( ω C Δ 2 ) Ω 2 e ich T ω P Ω P + e ich T ( ω P Δ ) Ω S e ich T ω C Ω C + e ich T ( ω C Δ 2 ) Ω D 0 )

Für EIT drehen wir unser System im korotierenden Rahmen, betrachten die Von-Neumann-Gleichung mit einigen zusätzlichen Zerfallstermen ( D P D T = [ H , P ] L P ), lösen Sie das System im stationären Zustand ( D P D T 0 ) unter Verwendung der adiabatischen Näherung ( ρ 11 1 , ρ 22 ρ 33 0 ). Wenn wir dann eine Lösung für die Kohärenzterme haben, sehen wir, dass unsere Ausgabe ρ 12 und wir haben unsere Antwort.

Für dieses System ist es für mich nicht offensichtlich, ob es einen mitrotierenden Rahmen gibt, der die Zeitabhängigkeiten der Gleichung beseitigt.

Da mein Ziel nun darin besteht, eine analytische Lösung für das Verhalten der Polarisation im stationären Zustand zu haben (damit ich herausfinden kann, welche Art von Ausgangslicht diese Wechselwirkung erzeugt), ist mein Gedanke, dass ich vielleicht ohne die Dichtematrixelemente auflösen kann ihre Zeitabhängigkeit herausdrehen. Außerdem habe ich mich entschieden, im gleichrotierenden Rahmen zu bleiben, einfach um die Diagonale in Bezug auf die Verstimmung zu setzen:

Verwendung der standardmäßigen gleichläufigen Rotation:

U = ( e ich T ω P 0 0 0 e ich T ω C 0 0 0 1 )

Wir können unseren effektiven Hamiltonian schreiben ( H e F F = ich D U D T . U + U . H . U ):

H e F F = ( 0 0 Ω P + e ich T ( ω P ω S ) Ω S 0 Δ P Ω C Ω P + e ich T ( ω P ω S ) Ω S Ω C Δ C )

Beachten Sie, dass mein Hamiltonian immer noch zeitabhängig ist. Jetzt bin ich mir nicht sicher, ob dieser Schritt richtig ist, aber ich werde nichts an der Standardprozedur zum Lösen der Kohärenzterme ändern. Das heißt, ich werde das einfach sagen D P D T 0 und lösen Sie auf (was ein bisschen verrückt erscheint ... sollte der stationäre Zustand nicht zeitabhängig sein?) Fahren Sie mit diesen naiven Schritten fort:

D P D T = [ H ( T ) , P ] + L P D P D T 0 [ H ( T ) , P ] + L P 0

Wir erhalten ein Gleichungssystem der Form MP = b. Die Matrix M ist:

( 0 Ω C + e ich T Δ 2 Ω D 0 0 0 Ω P e ich T Δ 1 Ω S Ω C + e ich T Δ 2 Ω D Δ 0 0 0 0 0 0 0 Ω P e ich T Δ 1 Ω S Ω C + e ich T Δ 2 Ω D 0 0 0 Ω P + e ich T Δ 1 Ω S Γ 2 Δ 0 0 0 0 Ω C e ich T Δ 2 Ω D 0 Γ 2 + Δ 0 Ω P e ich T Δ 1 Ω S 0 0 0 0 Γ 2 + Δ )

und P und b sind

P = ( ρ 12 ρ 13 ρ 21 ρ 23 ρ 31 ρ 32 ) , B = ( 0 Ω P + Ω S e ich Δ 1 T 0 0 Ω P Ω S e ich Δ 1 T 0 )

Jetzt habe ich eine Lösung für meine Kohärenzterme ρ 13 ( T ) , ρ 12 ( T ) :

ρ 13 ( T ) = 2 e 2 ich Δ 1 T ( Ω S + Ω P e ich Δ 1 T ) 2 ( Ω P + Ω S e ich Δ 1 T ) D 1 ρ 23 ( T ) = 2 ( Ω S + Ω P e ich Δ 1 T ) ( Ω P + Ω S e ich Δ 1 T ) ( Ω D + Ω C e ich Δ 2 T ) D 2
D 1 = Γ ( Ω C 2 + Ω D 2 ) + 2 Δ ( Ω C 2 + Ω D 2 Ω P 2 Ω S 2 ) + 2 Ω C Ω D ( Γ + 2 Δ ) cos ( Δ 2 T ) 4 Δ Ω P Ω S cos ( Δ 1 T ) D 2 = e ich T ( Δ 1 + Δ 2 ) ( Γ ( Ω C 2 + Ω D 2 + Ω P 2 + Ω S 2 ) + 2 Δ ( Ω C 2 Ω D 2 + Ω P 2 + Ω S 2 ) ) + Ω C Ω D ( Γ 2 Δ ) e ich T ( Δ 1 + 2 Δ 2 ) + Ω P Ω S ( Γ + 2 Δ ) e ich T ( 2 Δ 1 + Δ 2 ) + Ω C Ω D ( Γ 2 Δ ) e ich Δ 1 T + Ω P Ω S ( Γ + 2 Δ ) e ich Δ 2 T

Jetzt können wir diese Begriffe verwenden, um FWM-Begriffe zu finden, die mich interessieren. Verwenden der Definition von Polarisierung:

P = N T R ( [ ρ A A ρ A B ρ A C ρ B A ρ B B ρ B C ρ C A ρ C B ρ C C ] [ 0 0 μ 13 0 0 μ 23 μ 13 μ 23 0 ] )

Wir können sehen, dass unsere Ausgangspolarisierbarkeit ist: P Ö u T = μ 31 ρ 13 + μ 32 ρ 23 + C . C .

Wenn ich jetzt die Terme herausziehe, die proportional sind e ich ω P T Und e ich ω S T , ich sehe KEINE Terme, die im Sonden- und Steuerfeld linear sind (nach einer Taylor-Erweiterung). Das ist für mich sehr seltsam, da es einige Quellen gibt , die wie ein Frequenzstrahlteiler wirken sollten, was bedeutet, dass ich proportionale Terme haben sollte Ω S Ω C Ω D Und Ω P Ω C Ω D .

Vielleicht kann mir jemand sagen, was ich falsch mache? Vielleicht vermassele ich den zeitabhängigen Teil der Lösung? Vielleicht gibt es einen anderen Rahmen für die Arbeit mit Vier-Wellen-Mischung?

Beachten Sie, dass der Hamiltonian H E ICH T wie Sie es in Ihrer Anfangsgleichung geschrieben haben (eine Summe von Kosinussen, die eine Matrix multiplizieren), ist definitiv falsch und nicht mit der nächsten Gleichung vereinbar. Sie haben die RWA korrekt auf einen korrekten Hamiltonoperator (nicht gezeigt) angewendet und ein korrektes Ergebnis in Ihrer zweiten Gleichung erhalten, aber der angegebene Ausgangspunkt ist falsch. Beachten Sie auch, dass es sich um "gegenläufige" Begriffe handelt, nicht um "gegenläufige".

Antworten (1)

Ich hoffe, es ist in Ordnung, wenn ich nicht versuche, Fehler in Ihrem Ansatz zu finden, sondern nur beschreibe, wie ich das Problem auf ähnliche Weise lösen würde. Übrigens denke ich, dass Sie eigentlich in die richtige Richtung gehen, aber die gelegentlichen Tippfehler und Ungereimtheiten haben es mir etwas schwer gemacht, die problematischen Stellen herauszugreifen.

Zuallererst scheint es, wie in einem der Kommentare erwähnt, einige Fehler im ursprünglichen Hamiltonian zu geben. Ich glaube, was du schreiben wolltest, ist so etwas wie

H ( T ) = ( ω 1 0 0 0 ω 2 0 0 0 ω 3 ) + ( 0 0 μ 13 0 0 μ 23 μ 13 μ 23 0 ) E ( T ) ,
Wo
E ( T ) = E P cos ( ω P T + ϕ P ) + E S cos ( ω S T + ϕ S ) + E C cos ( ω C T + ϕ C ) + E D cos ( ω D T + ϕ D ) .
Hier gehe ich davon aus μ 13 Und μ 23 real sind, was ohne Verlust der Allgemeinheit möglich ist.

Ich denke auch, dass es sinnvoll ist, das Problem so weit wie möglich zu vereinfachen, bevor man versucht, es zu verstehen. Also statt dem allgemeinen Fall, wo die Umwandlung zwischen Frequenzen ω P Und ω S in beide Richtungen passieren kann, werde ich den einfacheren Fall betrachten, wo wir nur den anwenden ω P pumpen und nach Ausgang bei Frequenz suchen ω S (dh E S = 0 im Hamiltonian). Die gegenteilige Umwandlung kann ganz analog analysiert werden.

Bevor wir kompliziert rechnen, hier ist mein intuitives Verständnis dafür, wie Wellenmischungen hier ablaufen würden: in Abwesenheit von E S , können wir dies fast als das Standard-EIT-Problem bei Pumpen betrachten E P Und E C außer jetzt gibt es ein zusätzliches Signal mit Amplitude E D die leicht verstimmt ist (um Δ 1 ) aus ω C . Wenn diese Verstimmung klein ist, können wir uns die Kombination von anschauen E C Und E D fährt als etwas mit Frequenz ω C aber mit einer langsam variierenden Amplitude und Phase. Die langsame Variation moduliert den Transmissionskoeffizienten T des EIT-Systems für den schwachen Sondenton E P . Im Ausgabefeld sollten Sie also grundsätzlich mit einem Begriff rechnen T ( T ) E P cos ( ω P T + ϕ P ) , Wo T ( T ) ist periodisch mit der Frequenz Δ 1 . Durch Multiplizieren von Termen, die um oszillieren ω P Und Δ 1 , erhalten Sie natürlich eine Komponente bei ω P + Δ 1 = ω S . Und voila, das ist Ihr Frequenzmischprozess.

Nun, um dies tatsächlich abzuleiten. Zuerst führen wir die Frame-Transformation durch

U ( T ) = ( e ich ω 1 T 0 0 0 e ich ω 2 T 0 0 0 e ich ( ω 1 + ω P ) T )

Auf diese Weise ist der effektive Hamiltonoperator

H e F F ( T ) = ( 0 0 0 0 0 0 0 0 ω 3 ω 1 ω P ) + ( 0 0 μ 13 e ich ω P T 0 0 μ 23 e ich ( ω 2 ω 1 ω P ) T μ 13 e ich ω P T μ 23 e ich ( ω 2 ω 1 ω P ) T 0 ) E ( T ) ,

Wir wissen das ω 3 ω 1 ω P = Δ 2 Und ω 2 ω 1 ω P = ω C . Die einzigen Terme, die langsam veränderliche Komponenten aufweisen, die die Annäherung an rotierende Wellen überleben, sind diejenigen, die proportional zu sind μ 13 E P , μ 23 E C Und μ 23 E D (wir nehmen die Verstimmung zwischen ω C Und ω D ist nicht groß genug, um das fallen zu lassen μ 23 E D Begriff). Nach der Annäherung erhalten wir

H e F F ( T ) = ( 0 0 0 0 0 0 0 0 Δ 2 ) + 1 2 ( 0 0 μ 13 E P e ich ϕ P 0 0 μ 23 ( E C e ich ϕ C + E D e ich ( Δ 1 T + ϕ D ) ) μ 13 E P e ich ϕ P μ 23 ( E C e ich ϕ C + E D e ich ( Δ 1 T + ϕ D ) ) 0 ) .

Der Einfachheit halber werden wir die Kopplungsterme mit bezeichnen Ω P Und Ω C ( T ) , was uns gibt

H e F F ( T ) = ( 0 0 1 2 Ω P 0 0 1 2 Ω C ( T ) 1 2 Ω P 1 2 Ω C ( T ) Δ 2 ) .

Jetzt glaube ich, dass es ein wichtiger Schritt ist, das anzunehmen Δ 1 ist viel langsamer als die Zeitskala, auf der sich das System dem Gleichgewicht nähert. Auf diese Weise können Sie die Dichtematrix annehmen ρ ( T ) des Systems zu jeder Zeit T ist ungefähr der stationäre Zustand, der dem Hamilton-Operator entspricht H e F F ( T ) . Genau wie Sie gehe ich davon aus ρ 11 1 , ρ 22 ρ 33 0 und schreiben Sie das System stationärer Gleichungen für ρ 12 , ρ 13 Und ρ 32 (Die anderen drei nicht-diagonalen Elemente sind immer komplexe Konjugierte der ersten drei, also können wir uns die Mühe genauso gut sparen und uns vorerst nicht darum kümmern):

ich ρ 13 Ω C / 2 γ ρ 12 / 2 ich ρ 32 Ω P / 2 = 0 ich ρ 12 Ω C / 2 ( ich Δ 2 + Γ / 2 ) ρ 13 = ich Ω P / 2 ich ρ 12 Ω P / 2 + ( ich Δ 2 Γ / 2 ) ρ 32 = 0

Lassen Sie uns jetzt innehalten und beachten, dass diese Gleichungen denen, die Sie erhalten, ziemlich ähnlich sind, mit Ausnahme des Begriffs mit γ in meiner ersten Gleichung, die den Kohärenzverlust zwischen den Zuständen 1 und 2 beschreibt und die Sie nicht haben. Warum habe ich das hinzugefügt? Es stellt sich heraus, dass die Lösung dieser Gleichungen ohne diesen Dekohärenzmechanismus etwas schlecht erzogen ist. Wenn beispielsweise alle Laufwerke ausgeschaltet sind ( Ω P = Ω C = 0 ), ρ 12 ist unbestimmt - es fällt vollständig aus den Gleichungen heraus. Sehen Sie sich außerdem an, was passiert, wenn Sie das Verhältnis festlegen η := Ω P / Ω C und lass Ω C 0 . Das kannst du zeigen ρ 12 hat eine endliche Grenze, die davon abhängt η . Dies bedeutet, dass die Dichtematrixelemente als Funktionen der beiden Variablen Ω P Und Ω C zeigen ein ziemlich fieses Verhalten Ω P = Ω C = 0 . Am wichtigsten ist, dass sie nicht zu einer Taylor-Reihe erweitert werden können Ω P Und Ω C . Aber genau das möchten Sie bei Ihrer Suche nach den Begriffen für die Wellenmischung tun.

Eine andere Sichtweise ist, dass Sie zur Identifizierung des Wellenmischungsprozesses die Pumpen als Störung behandeln möchten. Aber der Erzeuger der freien Dynamik des Systems ist gewissermaßen entartet (er hat mehrere stationäre Zustände), und deshalb versagt die Störungstheorie.

Auf jeden Fall sollte der von mir hinzugefügte Dekohärenzterm dies beheben, da er sicherstellt, dass das freie, nicht angetriebene System nur einen einzigen stationären Zustand mit verschwindenden Matrixelementen außerhalb der Diagonale hat. Das Lösen der Gleichungen gibt uns dann

ρ 13 = ich Ω P 2 ich Δ 2 + Γ ich Ω P | Ω C | 2 γ ( 2 ich Δ 2 + Γ ) 2 + ρ 12 = Ω P Ω C γ ( 2 ich Δ 2 + Γ ) + ,
wobei wir Terme vierter und höherer Ordnung in den Antriebsstärken überspringen.

Jetzt erinnern wir uns daran, dass wir rückgängig machen müssen, um uns wieder in den Laborrahmen zu verwandeln U ( T ) . Das bedeutet multiplizieren ρ 12 von exp ( ich ( ω 2 ω 1 ) T ) Und ρ 13 von exp ( ich ω P T ) . Wir müssen uns auch ausdrücken Ω P Und Ω C ( T ) bezüglich E P , E C Und E D . Das sagt uns dann eine einfache Rechnung ρ 12 besteht aus Termen, die mit Frequenzen oszillieren ω 2 ω 1 Und ω 2 ω 1 + Δ 1 . Keines davon entspricht dem gewünschten ω S . Der Term erster Ordnung im Matrixelement ρ 13 schwingt bei ω P und ist nicht interessant. Der zweite Term hat andererseits eine Komponente bei der Frequenz Δ 1 + ω P = ω S . Dieser Begriff steht für den Mischprozess, den Sie anstreben. Explizit können wir es schreiben als

ρ 13 = ich μ 13 μ 23 2 E P E C E D e ich ( ω S T + ϕ P + ϕ D ϕ C ) γ ( 2 ich Δ 2 + Γ ) 2 +

Haftungsausschluss : Ich bin mir ziemlich sicher, dass Minuszeichen, Faktoren von 2 und andere lästige Dinge aufgrund meiner Unachtsamkeit dabei verloren gegangen sind. Aber ich glaube, das Wesentliche des ganzen Geschäfts sollte immer noch gelten. Ich hoffe, es wird nützlich sein.

Das ist eine ausgezeichnete Antwort. Willkommen auf der Seite!
Vielen Dank für Ihre Antwort! Ich bin sehr beeindruckt, wie klar das ist. Ich bin etwas perplex. Ich denke, wir haben beide dieselbe Methode verwendet (ich habe es im ursprünglichen Beitrag nicht erwähnt, aber ich habe auch ohne Erfolg versucht, Zerfallsterme einzufügen). Der einzige Unterschied, den ich sehen kann, besteht darin, dass Sie für die von der Sonde erzeugte FWM ohne Einbeziehung des ursprünglichen Signalfelds aufgelöst haben. Könnte es sein, dass die Einbeziehung sowohl der Sonde als auch des Signals zusätzliche Interferenzen in den Bewegungsgleichungen hervorruft? Vielleicht können diese Begriffe verschwinden (wahrscheinlich für einige relative Phasen zwischen den Feldern, aber nicht für alle)?
Ich halte Ihr adiabatisches Argument für die Aufnahme von Ω D könnte auch angewendet werden Ω S , also könnte ich vielleicht eine Version davon mit allen 4 Feldern ausarbeiten. (um zu sehen, ob sich die Kombination von der Einzelbehandlung unterscheidet) Ich werde Ihre Berechnungen sorgfältig durchgehen und bestätigen, dass ich Ihre Ergebnisse reproduzieren kann. Und dann, nachdem ich denke, werde ich versuchen und sehen, ob ich die gleiche Antwort bekomme, wenn ich versuche, sie zu machen Ω P Ω P e F F ( T ) = Ω P + Ω P e ich Δ 3 T + ϕ
Tippfehler oben: Ω P e F F ( T ) = Ω P + Ω S e ich Δ 3 T + ϕ S
Ich bin Ihre Frage im Detail durchgegangen und glaube, ich verstehe meine vorherige Verwirrung. Ich denke, das machen γ ungleich Null ist das Wichtigste, was mir gefehlt hat. Danke nochmal für deine Hilfe!
@StevenSagona: Sehr gerne. Was die Erweiterung auf bidirektionale Konvertierung betrifft, könnte es meiner Meinung nach ziemlich einfach sein. In der niedrigsten Ordnung aller Laufwerke sollte das Ganze linear sein E P Und E S , also würde ich erwarten, dass Sie einfach ersetzen P S , C D Und Δ 1 Δ 1 um den umgekehrten Prozess zu erhalten und die beiden Ergebnisse zu addieren. Natürlich erhalten Sie dann jede der Frequenzen ω P Und ω S sowohl das umgewandelte Signal als auch die direkte Reaktion auf die Sonde.