Lanczos Method for QRPA Strength Functions in Atomic Nuclei

Technische Samenvatting: Symmetrische Lanczos-methode voor QRPA-sterktefuncties

Probleemstelling
De berekening van ladingsveranderende nucleaire responsfuncties is essentieel voor het beschrijven van $\beta$-verval, dubbel-$\beta$-verval en collectieve excitaties in atoomkernen. Voor open-schil kernen met een gemiddelde massa en zware kernen, in het bijzonder vervormde systemen, biedt de Quasiparticle Random-Phase Approximation (QRPA) een noodzakelijk microscopisch kader. De dimensie van de QRPA-matrix groeit echter snel met de grootte van de quasiparticle-ruimte, waardoor expliciete constructie en volledige eigendecompositie onbetaalbaar duur worden.

Hoewel de Finite Amplitude Method (FAM) een matrix-vrije alternatief biedt door responsen te berekenen via geïnduceerde velden, lossen conventionele FAM-gebaseerde benaderingen het lineaire-respons-probleem frequentie-voor-frequentie op. Dit vereist het oplossen van een nieuw lineair systeem voor elke complexe frequentie op de gekozen contour. Bijgevolg wordt het reconstrueren van een sterktefunctie over een breed energie-interval computationeel kostbaar. Het artikel adresseert de behoefte aan een methode die spectrale informatie globaal over een uitgebreid frequentiebereik kan vastleggen zonder de kosten van herhaalde frequentie-specifieke oplossingen.

Methodologie
De auteurs ontwikkelen een symmetrisch Lanczos-raamwerk om charge-changing QRPA-sterktefuncties te benaderen. De methodologie verloopt via de volgende stappen:

1. Spectrale Herformulering: Vertrekkend vanuit de FAM lineaire-responsvergelijkingen, leiden de auteurs equivalente spectrale representaties af. In het geval dat QRPA-matrices $A$ en $B$ reëel zijn, formuleren zij een gereduceerd eigenwaardeprobleem met betrekking tot de matrixproducten $MK$ en $KM$, waarbij $M = A + B$ en $K = A - B$.
2. Symmetrische Lanczos-projectie: In tegenstelling tot de standaard QRPA-matrix die niet-Hermitisch is in de Euclidische inwendige product, is de gereduceerde operator $MK$ zelf-geadjungeerd met betrekking tot de $K$-inwendige product ($\langle x, y \rangle_K = x^T K y$). Deze eigenschap maakt de toepassing van een symmetrisch Lanczos-proces mogelijk.
3. Single-Run Benadering: Het Lanczos-algoritme genereert een tridiagonale matrix $T_m$ uit een enkele Krylov-subruimte run. De eigenwaarden van $T_m$ benaderen de gekwadrateerde excitatie-energieën ($\Omega_i^2$), en de bijbehorende Ritz-vectoren benaderen de QRPA-transitieamplituden.
4. Reconstructie van de Sterktefunctie: Met behulp van de Ritz-paren reconstrueert de methode de met Lorentziaanse breedte gesmeerde (Lorentzian-smeared) sterktefunctie over een breed energie-interval. Dit staat in contrast met frequentie-voor-frequentie solvers (FAM of GMRES) die afzonderlijk moeten worden uitgevoerd voor elk frequentiepunt.
5. Initialisatie: De Lanczos-recursie wordt geïnitialiseerd met een vector afgeleid van de externe veldcombinatie $f_+ = F^{20} + F^{02}$, wat ervoor zorgt dat de projectie de symmetrie-sectoren van het axiaal vervormde systeem respecteert.

Belangrijkste Bijdragen

• Afleiding van het Gereduceerde Eigenwaardeprobleem: Het artikel stelt de equivalentie vast tussen het volledige QRPA lineaire-respons-probleem en een gereduceerd eigenwaardeprobleem met betrekking tot $MK$ en $KM$, wat het gebruik van symmetrische Lanczos-projectie vergemakkelijkt.
• Matrix-vrije Implementatie: De methode vermijdt de assemblage van de volledige QRPA-matrix. Matrix-vector producten (matvecs) worden berekend met behulp van de geïnduceerde veldformalisme, vergelijkbaar met FAM, maar toegepast op de gereduceerde operatorstructuur.
• Efficiëntie in Brede Spectrale Scans: De benadering toont aan dat een enkele Krylov-run de dominante spectrale informatie kan vastleggen die nodig is voor een volledig sterkteprofiel, wat een significant alternatief biedt voor puntgewijze frequentieberekeningen.

Numerieke Resultaten
De methode werd getest op twee kernen: de medium-zware $^{112}\text{Sn}$ en de zware zeldzame aarde $^{150}\text{Nd}$, met gebruik van de Skyrme energie-densiteitsfunctionaal SkM*.

• Referentievalidatie: De auteurs hebben eerst de Generalized Minimal Residual (GMRES) methode als een efficiënte frequentie-voor-frequentie referentie vastgesteld. GMRES reproduceerde de geconvergeerde iteratieve FAM (IFAM) sterkteprofielen, maar vereiste minder iteraties (bijv. 2347 vs. 3781 voor $^{112}\text{Sn}$).
• Lanczos Nauwkeurigheid: Vergelijking van de Lanczos-benadering met de GMRES-referentie:
  • Voor $^{112}\text{Sn}$ verbeterde het verhogen van de Krylov-dimensie van $m=50$ naar $m=100$ de overeenstemming met de referentie aanzienlijk, waarbij zowel brede distributies als lokale piekstructuren werden gevangen.
  • Voor $^{150}\text{Nd}$ verminderde het verhogen van $m$ van 100 naar 200 eveneens de discrepanties.
  • Kullback-Leibler (KL) divergentie-metrieken bevestigden dat het verhogen van de Lanczos-dimensie de fout op distributieniveau substantieel verminderde (bijv. de KL-divergentie voor $^{112}\text{Sn}$ daalde van $2.89 \times 10^{-2}$ bij $m=50$ naar $2.34 \times 10^{-3}$ bij $m=100$).
• Computationele Kosten: De Lanczos-methode bleek efficiënter dan frequentie-voor-frequentie solvers voor brede spectrale scans. Een 100-staps Lanczos-run voor $^{112}\text{Sn}$ vereiste ongeveer 600 matvec-equivalente operaties, vergeleken met 2347 voor de GMRES-referentie. Vergelijkbare besparingen werden waargenomen voor $^{150}\text{Nd}$.

Betekenis en Claims
Het artikel claimt dat de symmetrische Lanczos-projectie een efficiënte en nauwkeurige benadering biedt voor QRPA-sterktefunctieberekeningen wanneer spectrale informatie vereist is over een uitgebreid frequentiebereik. De methode behoudt de matrix-vrije efficiëntievoordelen van FAM-achtige berekeningen, terwijl het profiteert van de bredere spectrale informatie die beschikbaar is via Lanczos-projectie.

De auteurs positioneren dit werk als een praktisch alternatief voor zware en vervormde kernen waar expliciete QRPA-matrices te groot zijn om te bouwen en herhaalde frequentie-voor-frequentie oplossingen computationeel onhaalbaar worden. Het artikel claimt niet de frequentie-specifieke solvers te vervangen in gevallen waar slechts enkele discrete frequenties nodig zijn, maar richt zich op het optimaliseren van de reconstructie van volledige respons-profielen. Toekomstig werk suggereert de uitbreiding van de methode naar aanvullende externe velden en grotere systematische berekeningen.