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

Titolo: Un metodo agli elementi finiti di Crank-Nicolson rilassato, completamente disaccoppiato e che preserva la struttura, per il modello Gross-Pitaevskii-Poisson.

1. Il Problema

Il lavoro si concentra sullo sviluppo di metodi numerici efficienti per risolvere il sistema accoppiato Gross-Pitaevskii-Poisson (GPP), utilizzato per modellare la dinamica dei plasmi ultrafreddi. Il sistema è descritto dalle seguenti equazioni in $\mathbb{R}^d \times (0, T]$:

$$\begin{aligned} i\partial_t\psi_{\pm} &= \left[ -\frac{1}{2}\nabla^2 + (g|\psi_{\pm}|^2 + G|\psi_{\mp}|^2) \pm q\phi \right] \psi_{\pm}, \\ \nabla^2\phi &= -4\pi q (|\psi_{+}|^2 - |\psi_{-}|^2). \end{aligned}$$

Dove:

• $\psi_{\pm}$ sono le funzioni d'onda quantistiche (cationi e anioni).
• $\phi$ è il potenziale di Poisson auto-consistente.
• $g$ e $G$ sono i coefficienti di interazione intra-specie e inter-specie.
• $q$ è il parametro di carica.

Sfide Numeriche:
Il sistema presenta diverse difficoltà computazionali:

1. Non linearità multiple: Coesistenza di termini di auto-interazione cubica ($|\psi|^2\psi$), interazione incrociata e termini di densità di carica ($|\psi|^2$) nell'equazione di Poisson.
2. Accoppiamento: La forte interazione tra le equazioni di Schrödinger non lineari e l'equazione di Poisson rende difficile la risoluzione diretta.
3. Preservazione delle proprietà fisiche: È cruciale che lo schema numerico preservi le quantità conservate fisicamente, come la massa di ciascuna componente e l'energia totale (o una sua versione modificata), specialmente su lunghi intervalli temporali.
4. Efficienza: I metodi esistenti (come i metodi di splitting o gli schemi iterativi impliciti) spesso richiedono la risoluzione di sistemi accoppiati complessi o multiple equazioni lineari per passo temporale, risultando costosi.

2. Metodologia

Gli autori propongono un metodo agli elementi finiti (FEM) di Crank-Nicolson rilassato, completamente disaccoppiato e che preserva la struttura.

• Riformulazione con Variabili Ausiliarie:
  Per gestire le non linearità, vengono introdotte due variabili ausiliarie reali: $Z_{\pm} = |\psi_{\pm}|^2$. Questo permette di riscrivere il sistema originale in una forma equivalente dove i termini non lineari sono linearizzati rispetto alle variabili ausiliarie.
  Il sistema riformulato diventa:
  $$\begin{aligned} i\partial_t\psi_{\pm} &= -\frac{1}{2}\nabla^2\psi_{\pm} + gZ_{\pm}\psi_{\pm} + GZ_{\mp}\psi_{\pm} \pm q\phi\psi_{\pm}, \\ \nabla^2\phi &= -4\pi q (Z_{+} - Z_{-}), \\ Z_{\pm} &= |\psi_{\pm}|^2. \end{aligned}$$

• Schema Temporale e Spaziale:

  • Tempo: Viene utilizzato uno schema di Crank-Nicolson con tecnica di rilassamento. Le variabili non lineari ($Z$) e il potenziale ($\phi$) sono valutati in punti medi o aggiornati in modo iterativo ma lineare, evitando la risoluzione di sistemi non lineari.
  • Spazio: Discretizzazione mediante Elementi Finiti ($P_k$) su un dominio limitato con condizioni al contorno periodiche.

• Algoritmo Disaccoppiato:
  A ogni passo temporale, il metodo richiede la risoluzione di:

  1. Due sistemi lineari indipendenti per le funzioni d'onda $\psi_{+}$ e $\psi_{-}$.
  2. Una singola equazione di Poisson per $\phi$.
  3. Due equazioni lineari semplici per aggiornare le variabili ausiliarie $Z_{\pm}$.
     Non sono necessarie iterazioni interne (come nel metodo di Newton) né la risoluzione di sistemi algebrici accoppiati di grandi dimensioni.

3. Contributi Chiave

1. Completamente Disaccoppiato e Lineare: A differenza dei metodi di splitting esistenti (che possono richiedere fino a 7 sistemi lineari per passo) o degli schemi impliciti non lineari, questo metodo risolve solo 3 sistemi lineari indipendenti per passo, garantendo un'efficienza computazionale superiore.
2. Preservazione della Struttura (Structure-Preserving): Lo schema conserva esattamente, a livello discreto:
   • La massa di ciascuna specie ($M(\psi_{\pm})$).
   • Un'energia totale modificata ($E$), che include i termini di interazione e il potenziale elettrostatico.
3. Analisi di Errore Rigorosa: Gli autori dimostrano la ben-postezza e l'unicità della soluzione discreta. Derivano stime di errore ottimali:
   • Ordine 2 in tempo ($O(\tau^2)$).
   • Ordine $(k+1)$ nello spazio ($O(h^{k+1})$) nella norma $L^2$, dove $k$ è il grado dei polinomi degli elementi finiti.
4. Generalizzabilità: Il metodo e l'analisi si estendono naturalmente al sistema GPE accoppiato (senza termine di Poisson, $q=0$).

4. Risultati

• Verifica Teorica: Le stime di errore confermate sperimentalmente mostrano una convergenza ottimale sia in tempo che nello spazio, in accordo con la teoria.
• Conservazione: I test numerici dimostrano che le fluttuazioni di massa ed energia rimangono al livello dell'errore di arrotondamento della macchina (machine precision), confermando la capacità del metodo di preservare le invarianti fisiche su lunghi tempi di simulazione.
• Dinamiche Complesse: Il metodo è stato testato su scenari fisici rilevanti, tra cui:
  • Onde di densità unidimensionali (dimostrando la capacità di catturare instabilità intrinseche del sistema).
  • Onde continue bidimensionali (caso GPE, $q=0$).
  • Onde di densità bidimensionali con accoppiamento Poisson.
    In tutti i casi, il metodo ha mostrato robustezza e accuratezza nel catturare i comportamenti dinamici, inclusi solitoni e pattern di densità.

5. Significato

Questo lavoro rappresenta un avanzamento significativo nella simulazione numerica dei plasmi ultrafreddi e dei condensati di Bose-Einstein multi-componente.

• Efficienza: Riduce drasticamente il costo computazionale rispetto ai metodi di splitting o agli schemi impliciti iterativi, rendendo possibili simulazioni più lunghe e con risoluzioni più fini.
• Affidabilità Fisica: La garanzia della conservazione della massa e dell'energia è fondamentale per studi a lungo termine dove gli errori di dissipazione numerica potrebbero portare a conclusioni fisiche errate.
• Versatilità: La strategia di riformulazione tramite variabili ausiliarie offre un modello promettente per affrontare altri sistemi accoppiati non lineari con termini di Poisson, superando le limitazioni dei metodi tradizionali.

In sintesi, il paper propone un metodo numerico che combina alta efficienza computazionale, accuratezza teorica e fedeltà fisica, risolvendo efficacemente le sfide poste dalla complessa interazione non lineare e dall'accoppiamento di Poisson nel modello GPP.