Spitze eines sich ausbreitenden Wellenpakets: Asymptotik jenseits aller Ordnungen einer Sattelpunktentwicklung

Dies ist eine technische Frage, die sich aus der Abbildung eines nicht verwandten Problems auf die Dynamik eines nicht-relativistischen massiven Teilchens in 1 + 1-Dimensionen ergibt. Dieses Problem wird bei der Asymptotik von einem Term jenseits aller Ordnungen einer Sattelpunktentwicklung (singuläre Terme einer asymptotischen Reihe) dominiert, wie beim Problem der Lebensdauer eines gebundenen Zustands in einer negativen 1+0-Kopplung ϕ 4 Spielzeugmodell.

Betrachten Sie ein Teilchen mit einer anfänglichen (normalisierten) Wellenfunktion

ψ 0 ( X ) = e ( X + e X ) / 2
Diese spezifische Form definiert eine natürliche Einheit für X , beachte die doppelexponentielle Asymptotik von ψ 0 ( X ) als X .

Zeitentwicklung unter dem Hamiltonian H = 1 2 X 2 transformiert die Wellenfunktion zu (mit dem Lehrbuchpropagator )

ψ ( X , T ) = ( 2 π ich T ) 1 / 2 e ich ( X X ' ) 2 / ( 2 T ) ψ 0 ( X ' ) D X '

Meine Frage betrifft die Asymptotik dieses Integrals , insbesondere die führende Front, die sich nach links ausbreitet. Hier bin ich an die Wand gefahren:

Die Sattelpunkterweiterung in T 1 gibt

T | ψ ( X , T ) | 2 e e X [ 1 + ( e 3 X 2 e 2 X ) T 1 / 8 + Ö ( T 2 ) ]
was gut konvergiert (numerisch überprüft) für X 1 , kann aber die Bestellbedingungen nicht erfassen e X / T die bei Negtavie über die doppelte Exponentialfunktion dominieren X .

Für T + die Lösung wird symmetrisch,

| ψ ( X , T ) | 2 = 1 T cosch ( π X / T )

Irgendwelche Ideen/Hinweise werden geschätzt.

EDIT: Das Ergebnis wurde in einer Veröffentlichung verwendet: Phys. Rev. Lett. 109, 216801 (2012) ; Open-Access- Version , siehe Gl. (9).

Kleiner Tippfehler in der Frage (v3): Die Quadratwurzel ( 2 π ich T ) 1 / 2 in der zweiten Formel sollte im Nenner stehen.
@Qmechanic Danke, das wurde behoben! Ich neige zu Tippfehlern :(

Antworten (1)

Nach einigem Ringen und einem nützlichen Hinweis eines Kollegen löste sich das Problem schließlich:

  1. Gehen zu Schwung Raum gibt

    ψ 0 k = ψ 0 ( X ) e ich k X D X = 2 ich k + 1 / 2 Γ ( ich k + 1 / 2 )

  2. Anwendung der Zeitentwicklung mit H = k 2 / 2 gibt

    ψ ( X , T ) = 1 2 π e ich k 2 T / 2 ich k X ψ 0 k D k

  3. Mit Groß- z asymptotische Entwicklung für Γ ( z + 1 ) und Identifizieren eines stationären Phasenpunkts in der Nähe k X / T Geben Sie die führende Ordnung an, die mit der Asymptotik von zusammenfällt T 1 Lösung:

    ψ ( X , T ) ( 2 / T ) e π X / T

  4. Schließlich wird der Vorfaktor wiederhergestellt, indem alle führenden Protokolle in der Erweiterung von aufbewahrt werden Γ ( z + 1 ) und das Auflösen nach dem stationären Punkt unter Verwendung der Lambert-W-Funktion (für die Mathematica die Asymptotik gut handhabt ) gibt die endgültige Antwort

    ψ ( X , T ) = ( 2 / T ) e π X / T ( 2 X / T ) π / T + Ö ( T 1 + ϵ )
    mit ϵ > 0 , Gültig für X T für beide T Groß und klein.

Danke an alle, die aufmerksam waren und mit Ratschlägen geholfen haben.