Resonance-based integrators for stochastic Schrödinger equations. Convergence and long-time error bounds

Technische Samenvatting: Resonantie-gebaseerde Integratoren voor Stochastische Schrödingervergelijkingen

Probleemstelling
Dit werk behandelt de numerieke benadering van de stochastische Schrödingervergelijking (SSE) gedreven door additieve $Q$-Wienerruis. De studie richt zich op twee specifieke gevallen: de lineaire vergelijking met een ruwe potentiaal ($F(\xi) = V(x)\xi$) en de kubisch niet-lineaire vergelijking ($F(\xi) = |\xi|^2\xi$). Een primaire uitdaging in dit domein is de lage temporele regulariteit van oplossingen die wordt veroorzaakt door stochastische convoluties, wat doorgaans de convergentieorde van standaard numerieke schema's beperkt. Bovendien zijn bestaande theoretische analyses voor dergelijke vergelijkingen grotendeels beperkt tot korte tijdschalen. Het artikel heeft tot doel resonantie-gebaseerde integratoren met lage regulariteit te ontwikkelen en te analyseren die hogere convergentieordes bereiken onder verminderde ruimtelijke regulariteitsaannames en, cruciaal, uniforme lange-tijd foutgrenzen vaststellen in regimes met kleine potentialen en ruis.

Methodologie
De voorgestelde methodologie is gebaseerd op Fourier-analyse en de constructie van resonantie-gebaseerde numerieke integratoren. De kernstrategie omvat:

1. Gedraaide Variabelen: Het transformeren van de oplossing $\xi(t)$ naar een gedraaide variabele $v(t) = e^{-it\Delta}\xi(t)$ om de dominante oscillerende termen te isoleren.
2. Exacte Integratie van Resonante Termen: De integratoren worden geconstrueerd door de dominante termen die voortkomen uit de Fourier-expansie van de niet-lineaire of potentiaalinteracties exact te integreren. Dit houdt in dat tijdintegralen van oscillerende functies (bijvoorbeeld $\int_0^\tau e^{is(\ell_1^2 + 2\ell_1\ell_2)} ds$) worden benaderd door de exacte integraal voor resonante gevallen (waar de fase verdwijnt) te behouden en niet-resonante gevallen te benaderen.
3. Regulariteits-compensatie Oscillatie (RCO) Techniek: Een nieuwe uitbreiding van de RCO-techniek, die eerder in deterministische settings werd gebruikt, wordt aangepast aan het stochastische kader. Deze techniek is essentieel om het verlies aan temporele regulariteit (doorgaans $O(\sqrt{\tau})$) veroorzaakt door stochastische convoluties te overwinnen. Door gebruik te maken van de compensatie van oscillaties in niet-resonante termen en gebruik te maken van een afsnijparameter voor hoge frequenties, leiden de auteurs verbeterde foutschattingen af.
4. Afsnijstrategie: Voor het kubisch niet-lineaire geval, waar sterke convergentie wordt gehinderd door een gebrek aan exponentiële integreerbaarheid, passen de auteurs een op afsnijden gebaseerde strategie toe in combinatie met technieken voor lage regulariteit om convergentieresultaten per pad af te leiden.

Belangrijkste Bijdragen en Resultaten

• Lineaire Vergelijking met Ruwe Potentiaal:

  • Convergentie: Een resonante integrator met lage regulariteit (SLR1) wordt geconstrueerd. Deze bereikt sterke convergentie van orde $\gamma \in (0, 1]$ in $L^2(\Omega; H^\sigma)$ voor oplossingen in $H^{\sigma+\gamma} Opmerkelijk is dat voor convergentie van de eerste orde in$H^\sigma$, de methode slechts $H^{\sigma+1}$ regulariteit vereist voor de beginwaarde en potentiaal, wat een verbetering is ten opzichte van de klassieke $H^{\sigma+2}$-vereiste.
  • Lange-tijd Grenzen: In een regime met een potentiaal van $O(\varepsilon^2)$ en ruis van $O(\varepsilon)$, wordt aangetoond dat de exacte oplossing uniforme momentgrenzen heeft tot tijden van $O(\varepsilon^{-2})$. Voor dit regime wordt een niet-resonante integrator (SNRLR1) ontworpen, die uniforme lange-tijd foutgrenzen van grootte $O(\varepsilon^2 \tau)$ in $L^2(\Omega; L^2)$ vertoont, een aanzienlijke verbetering ten opzichte van klassieke methoden.

• Kubisch Niet-lineaire Vergelijking:

  • Convergentie: Zowel resonante (SLR2) als niet-resonante (SNRLR2) integratoren worden geïntroduceerd. Convergentieresultaten per pad worden afgeleid bij lage regulariteit, waarbij orde $\gamma \in (0, 1)$ in $H^\sigma$ wordt bereikt voor oplossingen in $H^{\sigma+\gamma}$.
  • Lange-tijd Grenzen: In het zwak niet-lineaire regime ($O(\varepsilon^2)$ niet-lineariteit, $O(\varepsilon)$ ruis), wordt bewezen dat de niet-resonante integrator (SNRLR2) uniforme lange-tijd foutgrenzen per pad van grootte $O(\varepsilon^2 \tau^\delta)$ voldoet voor elke $\delta < 1$ tot tijden van $O(\varepsilon^{-2})$.

• Theoretische Innovatie:

  • Het artikel presenteert de eerste uitbreiding van de RCO-techniek naar stochastische dispersieve vergelijkingen. Deze uitbreiding stelt de auteurs in staat verbeterde foutschattingen te herstellen die de standaard $\sqrt{\tau}$-beperking opgelegd door stochastische ruis overwinnen, waardoor lange-termijnfouten met een factor $\varepsilon^2$ worden gereduceerd in vergelijking met andere integratoren met lage regulariteit.

Betekenis en Beweringen
De auteurs beweren dat dit werk de eerste vaststelling is van lange-tijd foutgrenzen voor integratoren met lage regulariteit die worden toegepast op stochastische dispersieve vergelijkingen. De betekenis ligt in:

1. Verminderde Regulariteitsvereisten: De methoden vereisen minder ruimtelijke regulariteit ($H^{\sigma+1}$ versus $H^{\sigma+2}$) om convergentie van de eerste orde te bereiken, waardoor ze toepasbaar zijn op ruwere potentialen en beginwaarden.
2. Lange-tijd Stabiliteit: De afleiding van uniforme foutgrenzen over tijdschalen van $O(\varepsilon^{-2})$ adresseert een gat in de literatuur, waarbij eerdere resultaten beperkt waren tot korte tijdschalen.
3. Stochastische RCO-uitbreiding: De succesvolle aanpassing van de RCO-techniek aan het stochastische setting biedt een nieuw analytisch hulpmiddel voor het behandelen van de temporele onregelmatigheid die inherent is aan stochastische convoluties.

Numerieke Verificatie
De theoretische resultaten worden bevestigd door middel van numerieke experimenten. Deze experimenten demonstreren:

• Sterke en per pad convergentiesnelheden van orde één voor het lineaire geval onder $H^1$-regulariteit.
• De begrensdheid van fouten voor de niet-resonante integrator (SNRLR1) over lange tijdsintervallen ($O(\varepsilon^{-2})$), in tegenstelling tot de foutgroei die wordt waargenomen bij de resonante (SLR1) en klassieke exponentiële methoden.
• De $O(\varepsilon^2)$-reductie in foutgrootte voor het niet-resonante schema in het regime met kleine potentiaal/ruis.
• Vergelijkbare convergentie en lange-tijdgedrag voor het kubisch niet-lineaire geval met behulp van de voorgestelde integratoren.