A fully decoupled and structure-preserving relaxation Crank--Nicolson finite element method for Gross--Pitaevskii--Poisson model

Probleemstelling

Het artikel richt zich op de ontwikkeling van efficiënte numerieke methoden voor het oplossen van het gekoppelde Gross-Pitaevskii-Poisson (GPP) systeem. Dit model beschrijft de dynamica van ultrakoude plasma's en omvat twee interagerende kwantumgolffuncties ($\psi_+$ en $\psi_-$, respectievelijk voor cationen en anionen) die gekoppeld zijn via een zelfconsistent Poisson-potentiaal ($\phi$).

De belangrijkste uitdagingen bij het numeriek oplossen van dit systeem zijn:

1. Niet-lineariteit: Het systeem bevat complexe niet-lineaire interactietermen (zelfinteractie en kruisinteractie) en niet-lineaire ladingsdichtheidstermen in de Poisson-vergelijking.
2. Koppeling: De golffuncties en het potentiaal veld zijn sterk aan elkaar gekoppeld, wat traditioneel leidt tot gekoppelde, niet-lineaire algebraïsche systemen die moeilijk en duur op te lossen zijn.
3. Behoudswetten: Het fysieke model behoudt specifieke invarianten, namelijk de massa van elke component en de totale energie. Een goede numerieke methode moet deze behoudswetten ook op het discrete niveau garanderen om fysisch correcte resultaten te verkrijgen op lange tijdschalen.

Bestaande methoden gebruiken vaak iteratieve procedures voor de niet-lineariteiten of splitsingsmethoden die meerdere lineaire systemen per tijdstap vereisen, wat de rekentijd aanzienlijk verhoogt.

Methodologie

De auteurs stellen een volledig ontkoppelde, structuurbehoudende relaxatie Crank-Nicolson eindige-elementenmethode (FEM) voor. De kern van de aanpak bestaat uit de volgende stappen:

1. Herschrijving met hulpparameters:
   Er worden twee nieuwe hulpparameters geïntroduceerd: $Z_+ = |\psi_+|^2$ en $Z_- = |\psi_-|^2$. Hierdoor worden de vier niet-lineaire interactietermen en de twee niet-lineaire ladingsdichtheidstermen herschreven in een equivalente vorm. Dit transformeert het oorspronkelijke niet-lineaire probleem in een systeem dat lineair is in de onbekenden op het nieuwe tijdstip.

2. Relaxatie Crank-Nicolson Schema:
   Voor de tijdsdiscretisatie wordt een relaxatie Crank-Nicolson schema toegepast. In plaats van de niet-lineariteiten op het huidige tijdstip direct te evalueren, worden ze benaderd door een lineaire combinatie van waarden op het huidige en vorige tijdstip (via de $Z$-variabelen). Dit maakt het mogelijk om het systeem lineair te houden zonder iteraties.

3. Volledige Ontkoppeling:
   Het resulterende numerieke schema is volledig ontkoppeld. Op elke tijdstap moeten slechts twee onafhankelijke lineaire systemen worden opgelost voor de golffuncties $\psi_+$ en $\psi_-$, en één Poisson-vergelijking voor het potentiaal $\phi$. Er is geen noodzaak om een groot, gekoppeld niet-lineair systeem op te lossen of iteratieve solvers te gebruiken.

4. Ruimtelijke Discretisatie:
   De ruimtelijke discretisatie gebeurt met de eindige-elementenmethode (FEM) met polynomen van graad $k$.

Belangrijkste Bijdragen

De paper introduceert een aantal significante innovaties:

• Efficiëntie door ontkoppeling: In tegenstelling tot bestaande splitsingsmethoden (zoals time-splitting sine spectral methods) die tot zeven lineaire systemen per tijdstap vereisen, vereist deze methode slechts twee lineaire systemen voor de golffuncties en één Poisson-vergelijking. Dit leidt tot een aanzienlijke reductie in rekentijd.
• Structuurbehoud: De methode behoudt strikt de discrete massa van elke component en een gemodificeerde discrete energie op het volledig discrete niveau. Dit is cruciaal voor het langdurig stabiel houden van de simulatie van fysieke systemen.
• Rigoureuze Foutanalyse: De auteurs bewijzen de welgesteldheid (well-posedness) en uniciteit van het schema. Ze leiden optimale foutschattingen af, bewijzend dat de methode tweede-orde nauwkeurigheid in de tijd en optimale $(k+1)$-de orde convergentie in de ruimte (in de $L^2$-norm) bereikt.
• Uniciteit en Lineariteit: Het bewijs toont aan dat het lineaire systeem op elke stap een unieke oplossing heeft, mits de initiële massa's compatibel zijn.

Resultaten

De theoretische resultaten worden gevalideerd door uitgebreide numerieke experimenten:

1. Convergentie: Numerieke tests in zowel 1D als 2D tonen aan dat de methoden de voorspelde convergentieorden halen: tweede orde in de tijd en $(k+1)$-de orde in de ruimte voor $P_k$ elementen.
2. Behoudswetten: Grafieken tonen aan dat de veranderingen in totale massa en energie op het niveau van machine-precisie blijven (rond-off fouten), wat bevestigt dat de behoudswetten perfect worden behouden.
3. Dynamisch Gedrag: De methode slaagt erin complexe dynamische fenomenen nauwkeurig weer te geven, waaronder:
   • Densiteitsgolven in 1D.
   • Continue golven in 2D.
   • Instabiliteitsverschijnselen (zoals modulatie-instabiliteit) die optreden bij specifieke parameterinstellingen, wat de robuustheid van de methode onderstreept.

Significantie

Deze studie is van groot belang voor de numerieke analyse van niet-lineaire kwantumsystemen met langeafstandsinteracties.

• Computational Efficiency: Door de noodzaak van iteraties en gekoppelde systemen te elimineren, maakt de methode het mogelijk om grotere en langdurigere simulaties van ultrakoude plasma's en multi-component Bose-Einstein condensaten (BECs) uit te voeren.
• Fysische Betrouwbaarheid: Het garanderen van behoudswetten op het discrete niveau voorkomt artefacten die vaak voorkomen bij standaard methoden, waardoor de resultaten fysisch consistenter zijn.
• Algemene Toepasbaarheid: Hoewel de methode specifiek is ontwikkeld voor het GPP-systeem, kan de relaxatie-techniek met hulpparameters ook worden toegepast op andere gekoppelde niet-lineaire Schrödinger-systemen, zoals het gekoppelde GPE-systeem zonder Poisson-koppeling.

Kortom, de auteurs bieden een krachtig, snel en fysiek correct alternatief voor bestaande numerieke methoden voor een complexe klasse van partiële differentiaalvergelijkingen.