Untersuchung geometrischer Phaseneffekte in Kohlenstoffnanoröhren-Nanodraht-Resonatoren

Zusammenfassung

Die geometrische Phase ist eine zusätzliche Phasenentwicklung in der Wellenfunktion von Schwingungen, die potenziell in einem breiten Spektrum von Wissenschaft und Technologie anwendbar ist. Die Eigenschaften der geometrischen Phase im gequetschten Zustand für einen auf Kohlenstoffnanoröhren basierenden Nanodraht-Resonator wurden mittels der Invarianten-Operatormethode untersucht. Die Einführung eines linearen Invariantenoperators, der nützlich für die Behandlung komplizierter zeitabhängiger Hamilton-Systeme ist, ermöglichte die Herleitung der analytischen Formel der geometrischen Phase. Anhand dieser wurde das Zeitverhalten der geometrischen Phase anhand relevanter Illustrationen analysiert. Der Einfluss der Quetschparameter auf die Entwicklung der geometrischen Phase wurde untersucht. Die geometrische Phase oszilliert im Großen und Ganzen, und die Einhüllende dieser Oszillation nimmt mit der Zeit zu. Die Zunahme-Rate der geometrischen Phase ist groß, wenn Parameter wie die klassische Amplitude der Oszillation, der Dämpfungsfaktor und die Amplitude der Treibkraft groß sind. Wir bestätigten einen sehr starken Anstieg der geometrischen Phase über die Zeit, wenn die Winkel Frequenz des Systems sich der Resonanzwinkelfrequenz nähert. Unsere Entwicklung bezüglich der Eigenschaften der geometrischen Phase ist entscheidend für das Verständnis topologischer Merkmale in Nanodrahtschwingungen.

Einleitung

Mechanische Schwingungen kleinster Resonatoren, wie kohlenstoffnanoröhrenbasierter (CNT-basierter) Nanodrähte, halbleitender Nanodrähte, Graphene und levitierter Partikel, sind seit über einem Jahrzehnt ein Hauptforschungsthema in der Nanowissenschaft. Aktive Forschungen zu elektromechanischen Schwingungen von Nanodraht-Resonatoren, angetrieben durch eine externe periodische Kraft, wurden sowohl theoretisch als auch experimentell durchgeführt. Insbesondere CNT-basierte Nanodraht-Resonatoren haben aufgrund ihrer außergewöhnlichen Empfindlichkeit mit hohen Gütefaktoren gegenüber kleinen Störungen aus der Umgebung erhebliches Interesse als nanoskalige mechanische Geräte geweckt. Freitragende CNT-basierte Nanodraht-Resonatoren sind vielversprechende Kandidaten für Apparate zur Messung vielfältiger physikalischer Größen wie EM-Wellen, kleine Kräfte, Massen, Temperaturen und Rauschen.

Analysen der quantalen Phasenentwicklung in Nanodrahtschwingungen sind erforderlich, um zugrundeliegende Systemmerkmale theoretisch zu erhellen. Bezüglich quantaler Schwingungszustände CNT-basierter Nanodraht-Resonatoren tritt die geometrische Phase sowie die übliche dynamische Phase als ergänzende Phasenentwicklung auf. Die geometrische Phase ist ein anholonomes Quantenzustandsmerkmal, das in diversen Physikbereichen anwendbar ist. Analysen der geometrischen Phase können potenziell zur Charakterisierung von Nanoeigenschaften von Nanodrähten herangezogen werden, wie Resonanzprofile, starke Quantenschwingungen, Spannungsrelaxationsmechanismen, das Auftreten von Dirac-Magnetoplasmonen und die Topologie von Aharonov-Bohm-Schwingungen.

Die Untersuchung der geometrischen Phase im Zusammenhang mit nicht-adiabatischer Dynamik kann Einblicke für nanomechanische Systeme liefern, die für die Weiterentwicklung präziser Simulationsverfahren notwendig sind. Die Präparation, Manipulation und Detektion von Quantenzuständen sind wichtige Faktoren in Quantentechnologien. Ziel der vorliegenden Forschung ist es, Licht auf das Zeitverhalten der geometrischen Phase zu werfen, das in Quantenzuständen von Nanodrahtschwingungen auftritt. Um den Mechanismus CNT-basierter Nanodrahtschwingungen zu verstehen, werden wir die Zeitentwicklung der geometrischen Phase im gequetschten Zustand untersuchen, welcher ein klassisch-ähnlicher Quantenzustand wie der kohärente Zustand ist. Der Vorteil des gequetschten Zustands besteht darin, dass die Unschärfe eines Quadratur in diesem Zustand auf Kosten der Erhöhung der Unschärfe des anderen Quadratur erheblich reduziert werden kann, während eine solche Unschärfemodulation im kohärenten Zustand unmöglich ist. Insbesondere werden wir die Auswirkungen der Resonanz auf die geometrische Phase analysieren. Da die Resonanzenergie signifikant von der Energie des nicht-resonanten Zustands abweicht, ist das topologische Verhalten der Wellenfunktion nicht-trivial und kann erheblich von normalen Situationen abweichen. Der Einfluss der Änderung physikalischer Parameter und der Quetschparameter auf die Entwicklung der geometrischen Phase wird ebenfalls rigoros analysiert. Geometrische Phasen sind in dynamischen Systemen allgegenwärtig und können auf verschiedene moderne Technologien angewendet werden, wie Quantenberechnung, Intensitätsinterferometrien, photonisches Multitasking, Quanten-Sensing-Protokolle und Wellenstabilitätsmessungen.

Der Hamiltonian des Systems beinhaltet Zeitfunktionen, die mit der Dämpfung des Systems und der externen Treibkraft assoziiert sind. Daher handelt es sich um eine Art zeitabhängiger Hamilton-Systeme (TDHS), deren quantenmechanische Probleme bis kürzlich umfangreich untersucht wurden. Die Zeitfunktion im Hamiltonian eines TDHS kann in den meisten Fällen nicht von der Funktion der kanonischen Variablen getrennt werden, was die konventionelle Variablentrennungsmethode zur Lösung der Schrödinger-Gleichung unanwendbar macht. Eine alternative leistungsstarke Methode, die zur Überwindung dieser Schwierigkeit entwickelt wurde, ist die Invarianten-Operatormethode, eingeführt von Lewis und Riesenfeld. Diese Methode ist ein sehr nützliches mathematisches Werkzeug, wenn wir Quantenslösungen eines TDHS ableiten. Viele quantenmechanische Probleme, die durch TDHS beschrieben werden, werden basierend auf dieser Methode untersucht. Beispiele sind chaotische Teilchenstreuung, Lichtausbreitung in zeitabhängigen Medien, Kontrolle getriebener gefangener Elektronen und Nichtklassikalität von Quanten-Nanoelektronikkreisen. Es gibt eine Vielzahl weiterer Methoden für quantenmechanische Behandlungen von TDHS, die die unitäre Transformationsmethode, die Lie-algebraische Methode und die Hamilton-Schätzmethode umfassen.

Da das System ein TDHS ist, verwenden wir die Invarianten-Operatormethode, um Quantenslösungen des Systems zu erhalten. Ein linearer Invarianter Operator, der in Terms des Annihilationsoperators dargestellt wird, wird eingeführt. Während Annihilations- und Erzeugungsoperatoren aufgrund der Zeitabhängigkeit des Systems zeitabhängig dargestellt werden, können sowohl kohärente als auch gequetschte Zustände mit diesen Leiteroperatoren erhalten werden. Die geometrische Phase des Systems wird analytisch durch Nutzung der Wellenfunktion im gequetschten Zustand evaluiert. Die Zeitentwicklung der geometrischen Phase wird detailliert auf Basis ihrer Illustrationen mit diversen Parameterwahlen analysiert.

Methoden

Um die geometrische Phase zu untersuchen, müssen wir zunächst die klassische Bewegungsgleichung der Nanodrahtspitze aufstellen. Da die geometrische Phase in der quantenmechanischen Wellenentwicklung eines TDHS erscheint, ist es notwendig, Wellenfunktionen in einem spezifischen Quantenzustand abzuleiten, den wir handhaben. Wir betrachten den gequetschten Zustand wie im einführenden Teil erwähnt. Wellenfunktionen in diversen Quantenzuständen eines TDHS, einschließlich des gequetschten Zustands, können aus der Invarianten-Operatormethode gewonnen werden.

Die Bewegungsgleichung für die zeitabhängige Amplitude x für einen Biegemodus eines freitragenden Kohlenstoffnanoröhrchens mit effektiver Masse m ist gegeben durch:

ẍ + (ω₀/Q + ηx²)ẋ + (ω₀² + βx²)x = f_d cos(ωt)

wobei ω₀ die Resonanzwinkelfrequenz, Q der Gütefaktor, f_d die elektrostatische Treibkraft geteilt durch m, η der nichtlineare Dämpfungskoeffizient und β der Duffing-Parameter ist. Nehmen wir der Einfachheit halber an, dass die Verschiebung der Spitze hinreichend klein im Vergleich zur CNT-Drahtlänge ist. Dann können wir die nichtlinearen Terme in Gl. (1) vernachlässigen, was zu folgendem führt:

ẍ + (ω₀/Q)ẋ + ω₀²x = f_d cos(ωt)

Der Hamiltonian des Systems, der Gl. (2) liefert, ist gegeben durch:

Ĥ = e^(-γt) p̂²/(2m) + ½ m e^(γt) [ω₀² x̂² - 2f_d cos(ωt) x̂]

wobei γ = ω₀/Q. Die klassische Lösung von Gl. (2) setzt sich aus einer Komplementärfunktion X_c(t) und einer Partikulärlösung X_p(t) zusammen, gegeben durch:

X_c(t) = X_c,0 e^(-γt/2) cos(Ωt + φ)

X_p(t) = X_p,0 cos(ωt - δ)

wobei X_c,0 eine Konstante ist, Ω = √(ω₀² - γ²/4), φ eine beliebige Phase, und

X_p,0 = f_d / √((ω₀² - ω²)² + γ²ω²)

δ = tan⁻¹(γω / (ω₀² - ω²))

Die klassische Lösung im Impulsraum ist analog gegeben, wobei die Komplementärfunktion P_c(t) = m e^(γt) Ẋ_c(t) und die Partikulärlösung P_p(t) = m e^(γt) Ẋ_p(t) ist. Um die geometrische Phase des Systems zu untersuchen, müssen wir zunächst Quantenslösungen ableiten. Beachten Sie, dass der Hamiltonian des Systems in Gl. (3) explizit zeitabhängig ist. Um Quantenslösungen des Systems abzuleiten, verwenden wir die Invarianten-Operatormethode, welche ein nützliches Werkzeug bei der Behandlung solcher zeitvariabler Systeme ist. Ein Invarianter Operator Î des Systems kann aus der Liouville-von-Neumann-Gleichung abgeleitet werden, gegeben durch dÎ/dt = ∂Î/∂t + [Î, Ĥ]/(iℏ) = 0. Daher erhalten wir nach strenger Auswertung durch Einsetzen von Gl. (3) in diese Gleichung einen linearen Invariantenoperator der Form:

Î = Â e^(iΩt)

wobei  der Annihilationsoperator ist, gegeben durch:

 = (2ℏmΩ)^(-1/2) [ m(Ω + iγ/2) e^(γt/2) [x̂ - X_p(t)] + i e^(-γt/2) [p̂ - P_p(t)] ]

Die hermitesche Adjungierte von Gl. (9), †, ist der Erzeugungsoperator.

Wir können die Eigenwertgleichung von  als  |A⟩ = A |A⟩ ausdrücken. Durch Auswertung erhalten wir den Eigenwertausdruck:

A(t) = A(0) e^(-iΩt)

wobei A(0) = A₀ e^(-iφ) mit A₀ = [mΩ/(2ℏ)]^(1/2) X_c,0.

Während der kohärente Zustand |A⟩ der Eigenzustand von  ist, ist der gequetschte Zustand der Eigenzustand eines Operators B̂, gegeben durch:

B̂ = μ + ν†

wobei μ und ν komplexe Variablen sind, die die Gleichung |μ|² - |ν|² = 1 erfüllen. Wenn wir die Eigenwertgleichung von B̂ in der Form B̂ |B⟩ = B |B⟩ schreiben, ist |B⟩ der gequetschte Zustand. Durch Lösen dieser Gleichung im Konfigurationsraum erhalten wir die Wellenfunktion im gequetschten Zustand (Gleichung 16). Für μ=1 und ν=0 reduziert sich Gl. (16) auf die Wellenfunktion im kohärenten Zustand. Diese Wellenfunktion wird im nächsten Abschnitt verwendet, um die geometrische Phase im gequetschten Zustand abzuleiten.

Ergebnisse und Diskussion

Es ist bekannt, dass die Phase in der quantenmechanischen Wellenentwicklung die geometrische Phase sowie die dynamische Phase umfasst. Die geometrische Phase wurde erstmals 1984 von Berry für ein System entdeckt, das sich zyklisch mit adiabatischer Änderung entwickelt. Gemäß dem adiabatischen Theorem der Quantenmechanik bleibt ein momentaner Eigenzustand eines Quantenzustands in einer zyklischen Entwicklung im Parameterraum später auf demselben Zustand, während eine zusätzliche Akkumulation der Quantenphase auftritt, die Berry-Phase. Eine Verallgemeinerung der Berry-Phase, die nicht-adiabatische, nicht-zyklische und/oder nicht-unitäre Entwicklung des Quantensystems einschließt, ist die geometrische Phase.

Die geometrische Phase im gequetschten Zustand ist gegeben durch:

γ_G(t) = ∫₀ᵗ ⟨B(t')| i ∂/∂t' |B(t')⟩ dt' + γ_G(0)

Die Ableitung der Wellenfunktion nach der Zeit im Konfigurationsraum wird zu einer Ausdrückung mit den Funktionen f₁(t), f₂(t), f₃(t). Nach weiterer Auswertung durch Einsetzen in die obige Gleichung erhalten wir die vollständige geometrische Phase im gequetschten Zustand, gegeben durch Gleichung (28) mit den Gleichungen (23), (24) und (29)–(32).

Die Zeitentwicklung der geometrischen Phase wurde in Abbildungen 1–4 illustriert. Aus Abb. 1 sehen wir, dass die geometrische Phase oszilliert und die Einhüllende dieser Oszillation mit der Zeit zunimmt. Die Zunahme der Einhüllende ist größer, wenn A₀ groß ist. Das Oszillationsmuster wird allmählich unregelmäßiger, je größer die Werte von μ und ν werden. Zudem wird die Amplitude der Oszillation mit der Zeit größer.

Zeitentwicklung der geometrischen Phase für verschiedene Werte von A₀. Die Werte von (μ, ν) in den Grafiken sind (1, 0) für a, (√2, 1) für b und (√3, √2) für c. Verwendete Parameter: m=1, ω₀=1, ω=5, γ=0,35, f_d=1, ℏ=1, φ=0, γ_G(0)=0. Alle Parameter sind dimensionslos. Da A₀ in Abhängigkeit von der klassischen Amplitude X_c,0 der Komplementärfunktion gegeben ist, bestätigen die Grafiken, dass die geometrische Phase groß ist, wenn die Schwingungsamplitude hoch ist. Wir sehen auch, dass die Fluktuation von γ_G(t) groß wird, wenn die Werte μ und ν unter der Bedingung von Gl. (14) zunehmen.

Zeitentwicklung der geometrischen Phase für verschiedene Werte von γ. Der Wert von ω in den Grafiken ist 0,3 für a, 0,99 für b und 5 für c. Die gewählten Quetschparameter sind μ=√2 und ν=1; dies ergibt einen q-gequetschten Zustand zum Anfangszeitpunkt. Weitere Größen: m=1, ω₀=1, A₀=1, f_d=1, ℏ=1, φ=0, γ_G(0)=0. Wir bestätigen, dass die geometrische Phase groß ist, wenn der Dämpfungsfaktor γ in den meisten Fällen groß ist, aber nicht allen. Die Frequenz von Fall b liegt nahe der Resonanzfrequenz, während die von a und c weit davon entfernt sind. Die geometrische Phase für den resonanten Fall (b) nimmt sehr schnell über die Zeit zu.

a–c Diese Grafik entspricht Abb. 2, aber für den Fall, dass die gewählten Quetschparameter μ=√2 und ν=-1 sind, was einen p-gequetschten Zustand zum Anfangszeitpunkt ergibt. Da die Gesamtdarstellungen in diesem Fall nicht stark von den entsprechenden von Abb. 2 abweichen, können wir bestätigen, dass die Entwicklung von γ_G(t) nahezu unabhängig von den Arten der Quetschung ist, solange die Beträge von μ und ν sich nicht ändern.

Zeitentwicklung der geometrischen Phase für verschiedene Werte von f_d. Der Wert von ω in den Grafiken ist 0,3 für a und 5 für b. Verwendet: μ=√2, ν=1, m=1, ω₀=1, γ=0,3, A₀=1, ℏ=1, φ=0, γ_G(0)=0. Mit zunehmender Amplitude (f_d) der Treibkraft wird die geometrische Phase groß.

Die Quetscheffekte im gequetschten Zustand in Abhängigkeit vom Quetschparameter c = μ/ν wurden in Ref. [39] untersucht. Gemäß der Analyse in Ref. [39] entspricht der in Abb. 2 dargestellte gequetschte Zustand (c=√2) dem q-gequetschten Zustand zum Anfangszeitpunkt, während der in Abb. 3 (c=-√2) dem p-gequetschten Zustand entspricht. Durch Vergleich von Abb. 2 und 3 können wir schließen, dass die geometrische Phase im q-gequetschten Zustand nahezu dieselbe ist wie im p-gequetschten Zustand.

Die Effekte von γ auf die Entwicklung der geometrischen Phase lassen sich aus Abb. 2 und 3 bestätigen. Die geometrische Phase nimmt schneller zu, wenn γ groß ist. Durch Vergleich von Abb. 2a und 3a mit Abb. 2c und 3c bestätigen wir, dass die geometrische Phase etwas schneller variiert, wenn ω größer als die Resonanzwinkelfrequenz ist.

Das Zeitverhalten der geometrischen Phase am oder nahe dem Resonanzzustand des Systems könnte von großem Interesse sein. Abb. 2b und 3b zeigen, dass die geometrische Phase sehr schnell zunimmt, wenn ω nahe der Resonanzwinkelfrequenz liegt. Dies bedeutet, dass die Wellenfunktion in dieser Situation signifikant über die Zeit variiert, da die Größe der geometrischen Phase mit der Zeitvariation der Wellenfunktion zusammenhängt. Tatsächlich wird die Amplitude der Drahtschwingung im Resonanzzustand bemerkenswert verstärkt. Übrigens sind die Resonanzwinkelfrequenzen freitragender CNT-basierter Nanodraht-Resonatoren nicht nur hoch, sondern auch mit sehr hohen Gütefaktoren weit abstimmbar. Aus diesem Grund werden die Schwingungsmoden des Systems lange Zeit aufrechterhalten, bis sie vollständig gedämpft sind.

Abb. 4 zeigt, dass die geometrische Phase auch von der Amplitude der Treibkraft f_d beeinflusst wird. Mit zunehmendem f_d ist der Anstieg der geometrischen Phase in der Zeit schnell.

Schlussfolgerung

Wir haben die geometrische Phase im gequetschten Zustand für das System auf Basis der Quantendynamik mit der Schrödinger-Gleichung untersucht. Angesichts der Zeitabhängigkeit des Hamiltonians, der das System beschreibt, wurde die Invarianten-Operatormethode eingeführt, die ein potenzielles Werkzeug zur Ableitung von Quantenslösungen ist, wenn der Hamiltonian in Terms der Zeit beschrieben wird. Mittels dieser Methode wurde die analytische Formel der geometrischen Phase für die CNT-basierte Nanodrahtschwingung erhalten.

Eine detaillierte Analyse der Phaseneffekte, die für ein theoretisches Verständnis der mechanischen Schwingungen notwendig ist, wurde durchgeführt. Unsere Entwicklung der geometrischen Phase ist eine vollständig quantenbasierte mit rigorosen mathematischen Auswertungen. Die geometrische Phase ist empfindlich gegenüber Änderungen mechanischer Parameter und zeigt in einer großen Anzahl von Fällen Oszillationen. Der Einfluss der Quetschparameter auf die Entwicklung der geometrischen Phase wurde ebenfalls analysiert. Wir bestätigten eine starke Zunahme der geometrischen Phasenakkumulation über die Zeit nahe der Resonanzwinkelfrequenz.

Unsere Ergebnisse illustrieren das Zeitverhalten der geometrischen Phase, das in der Schwingung eines CNT-basierten Nanodrahts erscheint. Die Analyse der geometrischen Phase in dieser Arbeit ist wichtig für das Verständnis nicht nur topologischer Merkmale des Systems, sondern auch dynamischer Schwingungen anderer nanodrahtbasierter mechanischer Oszillatoren. Insbesondere haben wir Phaseneigenschaften des Resonanzzustands entwickelt, deren Klärung notwendig ist für die Anwendung des Systems in Quanteninformationstechnologien und anderen quantenbasierten Industrien. Die ähnliche Methode und das Framework, die in dieser Forschung verwendet wurden, können auch auf andere Nanosysteme erweitert werden, wie supraleitende Fabry-Pérot-Resonatoren, Nanocantilever und Qubit-Resonator-Atom-Hybridsysteme.

Abkürzungen

CNT:

Kohlenstoffnanoröhre

EM waves:

Elektromagnetische Wellen

TDHS:

Zeitabhängiges Hamilton-System