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:
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 ( ) ausgeschaltet ist, wird Licht des ersten Halbbildes in der Frequenz des zweiten Halbbildes erzeugt. Um dies visuell zu sehen, erwarten wir, dass wann eingeschaltet sind, wird Licht in der Frequenz erzeugt , und wir erhalten eine Vier-Wellen-Mischung (FWM), wie in der folgenden Abbildung angegeben:
Ich erwarte Bedingungen des Formulars zu sehen Und
Und eine ähnliche Lösung (vierstufiges statt dreistufiges System) wurde hier gelöst . In ihrem Vier-Ebenen-System haben sie Begriffe Und .
Mein Versuch
Ich kann eine analytische Lösung für bekommen Und , 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:
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. Und sind ähnlich, also die gleichen Regeln für bewerben für und ebenso für Und : )
Durch Entfernen der gegenläufigen Terme kann dies reduziert werden auf
Für EIT drehen wir unser System im korotierenden Rahmen, betrachten die Von-Neumann-Gleichung mit einigen zusätzlichen Zerfallstermen ( ), lösen Sie das System im stationären Zustand ( ) unter Verwendung der adiabatischen Näherung ( ). Wenn wir dann eine Lösung für die Kohärenzterme haben, sehen wir, dass unsere Ausgabe 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:
Wir können unseren effektiven Hamiltonian schreiben ):
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 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:
Wir erhalten ein Gleichungssystem der Form MP = b. Die Matrix M ist:
und P und b sind
Jetzt habe ich eine Lösung für meine Kohärenzterme :
Jetzt können wir diese Begriffe verwenden, um FWM-Begriffe zu finden, die mich interessieren. Verwenden der Definition von Polarisierung:
Wir können sehen, dass unsere Ausgangspolarisierbarkeit ist:
Wenn ich jetzt die Terme herausziehe, die proportional sind Und , 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 Und .
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?
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
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 Und in beide Richtungen passieren kann, werde ich den einfacheren Fall betrachten, wo wir nur den anwenden pumpen und nach Ausgang bei Frequenz suchen (dh 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 , können wir dies fast als das Standard-EIT-Problem bei Pumpen betrachten Und außer jetzt gibt es ein zusätzliches Signal mit Amplitude die leicht verstimmt ist (um ) aus . Wenn diese Verstimmung klein ist, können wir uns die Kombination von anschauen Und fährt als etwas mit Frequenz aber mit einer langsam variierenden Amplitude und Phase. Die langsame Variation moduliert den Transmissionskoeffizienten des EIT-Systems für den schwachen Sondenton . Im Ausgabefeld sollten Sie also grundsätzlich mit einem Begriff rechnen , Wo ist periodisch mit der Frequenz . Durch Multiplizieren von Termen, die um oszillieren Und , erhalten Sie natürlich eine Komponente bei . Und voila, das ist Ihr Frequenzmischprozess.
Nun, um dies tatsächlich abzuleiten. Zuerst führen wir die Frame-Transformation durch
Auf diese Weise ist der effektive Hamiltonoperator
Wir wissen das Und . Die einzigen Terme, die langsam veränderliche Komponenten aufweisen, die die Annäherung an rotierende Wellen überleben, sind diejenigen, die proportional zu sind , Und (wir nehmen die Verstimmung zwischen Und ist nicht groß genug, um das fallen zu lassen Begriff). Nach der Annäherung erhalten wir
Der Einfachheit halber werden wir die Kopplungsterme mit bezeichnen Und , was uns gibt
Jetzt glaube ich, dass es ein wichtiger Schritt ist, das anzunehmen ist viel langsamer als die Zeitskala, auf der sich das System dem Gleichgewicht nähert. Auf diese Weise können Sie die Dichtematrix annehmen des Systems zu jeder Zeit ist ungefähr der stationäre Zustand, der dem Hamilton-Operator entspricht . Genau wie Sie gehe ich davon aus , und schreiben Sie das System stationärer Gleichungen für , Und (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):
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 ( ), 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 und lass . Das kannst du zeigen hat eine endliche Grenze, die davon abhängt . Dies bedeutet, dass die Dichtematrixelemente als Funktionen der beiden Variablen Und zeigen ein ziemlich fieses Verhalten . Am wichtigsten ist, dass sie nicht zu einer Taylor-Reihe erweitert werden können Und . 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
Jetzt erinnern wir uns daran, dass wir rückgängig machen müssen, um uns wieder in den Laborrahmen zu verwandeln . Das bedeutet multiplizieren von Und von . Wir müssen uns auch ausdrücken Und bezüglich , Und . Das sagt uns dann eine einfache Rechnung besteht aus Termen, die mit Frequenzen oszillieren Und . Keines davon entspricht dem gewünschten . Der Term erster Ordnung im Matrixelement schwingt bei und ist nicht interessant. Der zweite Term hat andererseits eine Komponente bei der Frequenz . Dieser Begriff steht für den Mischprozess, den Sie anstreben. Explizit können wir es schreiben als
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.
Emilio Pisanty