MFGSB (ver. 1.0): Computer code for self-consistent mean-field calculations of atomic nuclei using Gaussian expansion method

Il codice informatico MFGSB (versione 1.0), disponibile nel repository dell'Università di Chiba, esegue calcoli di campo medio autoconsistenti per i nuclei atomici utilizzando il metodo di espansione gaussiana.

L'essenziale

Immagina di voler capire come è fatto un LEGO gigante che rappresenta il cuore di un atomo (il nucleo atomico). Questo "cuore" è fatto di migliaia di piccoli pezzi (protoni e neutroni) che si muovono, si attraggono e si respingono in un balletto complesso.

Il documento che hai condiviso presenta un nuovo manuale di istruzioni (un codice informatico chiamato MFGSB) creato da un ricercatore dell'Università di Chiba in Giappone. Questo programma è come un super-architetto virtuale capace di ricostruire e prevedere come si comportano questi nuclei atomici.

Ecco una spiegazione semplice, con qualche metafora per rendere il tutto più chiaro:

1. Cos'è questo "MFGSB"?

Pensa al nucleo atomico come a una folla di persone in una stanza buia che si tengono per mano. Per capire come si muovono, devi calcolare le forze che agiscono su ognuno.
Il codice MFGSB è un programma scritto in un linguaggio informatico classico (FORTRAN77) che fa questi calcoli in modo auto-consistente.

• Cosa significa "auto-consistente"? Immagina di disegnare una mappa del traffico. Se muovi un'auto, il traffico cambia, quindi devi ridisegnare la mappa. Se ridisegni la mappa, le auto si spostano di nuovo. Il programma fa questo ciclo infinite volte finché la mappa e le auto non si "accordano" perfettamente. Non usa approssimazioni facili, ma calcola tutto esattamente.

2. Il trucco magico: Il "Metodo dell'Espansione Gaussiana"

Il programma usa un metodo speciale chiamato GEM (Gaussian Expansion Method).

• L'analogia: Immagina di dover disegnare una montagna complessa. Potresti usare tanti piccoli mattoncini quadrati (come i pixel), ma il risultato sarebbe sgranato. Oppure potresti usare pennellate morbide e curve (come l'acquerello).
• Il metodo GEM usa queste "pennellate curve" (funzioni gaussiane) per descrivere la forma dei nuclei.
• Il vantaggio: È come se avessi un set di pennelli magici che funzionano bene per qualsiasi tipo di montagna, sia che sia una collina bassa o un picco altissimo. Non devi cambiare i pennelli se vuoi disegnare un nucleo diverso; lo stesso set di strumenti serve per quasi tutti i nuclei dell'universo. Questo fa risparmiare tempo e rende il programma molto versatile.

3. Cosa può fare questo programma?

Il codice è una vera e propria "fabbrica di simulazioni". Può fare diversi tipi di calcoli:

• Simulazioni statiche: Guarda il nucleo fermo (come una foto).
• Simulazioni dinamiche: Guarda come le particelle si legano tra loro (come un video).
• Simmetrie: Può guardare il nucleo da diverse angolazioni. Può vederlo come una sfera perfetta, come un pallone da rugby allungato (asse simmetrico) o in forme ancora più strane.

4. Cosa serve per farlo funzionare?

Per far girare questo "super-architetto", il tuo computer ha bisogno di:

• Molta memoria: Per il calcolo di base servono 8 GB di RAM (come se avessi un tavolo di lavoro enorme). Se vuoi fare calcoli più complessi, il tavolo deve raddoppiare di dimensioni ogni volta.
• Librerie speciali: Il programma ha bisogno di "utensili" già pronti (chiamati BLAS e LAPACK) per fare i calcoli matematici velocemente.

5. Come si ottiene?

Il codice è gratuito e disponibile online.

• Il pacchetto completo: È un file zippato di 16,6 GB. È enorme! È come scaricare l'intero archivio di un museo.
• L'opzione leggera: Se non hai spazio, puoi scaricare solo il "motore" del programma e scaricare i dati specifici (le "ricette" delle forze tra le particelle) separatamente, come se scaricassi solo gli ingredienti che ti servono per quella ricetta specifica.

In sintesi

Il documento annuncia la disponibilità di un potente strumento scientifico che permette ai fisici di "vedere" dentro i nuclei atomici con una precisione incredibile, usando un metodo matematico intelligente che non richiede di riscrivere le regole ogni volta che studiano un nuovo elemento. È come avere una macchina del tempo e uno scanner 3D per il mondo subatomico, tutto accessibile gratuitamente da Chiba.

Se qualcuno vuole usare i risultati di questo programma per una ricerca, deve solo citare i quattro articoli scientifici originali dove sono stati sviluppati i metodi matematici alla base.

Sommario tecnico

Riassunto Tecnico: MFGSB (ver. 1.0)

Titolo: Codice informatico per calcoli di campo medio autoconsistente di nuclei atomici utilizzando il metodo di espansione gaussiana.
Autore: H. Nakada (Università di Chiba, Giappone).
Data: 5 febbraio 2026 (preprint arXiv).

1. Il Problema

Il calcolo della struttura dei nuclei atomici richiede la risoluzione di equazioni di campo medio autoconsistente (SCMF - Self-Consistent Mean-Field). Le sfide principali in questo campo includono:

• La necessità di gestire interazioni nucleari efficaci complesse (inclusi termini tensoriali e di scambio di Coulomb) senza approssimazioni eccessive.
• La difficoltà nel descrivere efficientemente l'asintotica dipendente dall'energia delle funzioni d'onda delle particelle singole o dei quasiparticelle.
• La necessità di adattare i parametri delle funzioni di base a diversi nuclidi e deformazioni nucleari, un processo spesso laborioso e soggetto a errori.
• La gestione computazionalmente costosa del moto del centro di massa (c.m.) e dei termini di interazione a due corpi.

2. Metodologia

Il documento presenta MFGSB, un codice scritto in FORTRAN77 progettato per eseguire calcoli SCMF basati sul Metodo di Espansione Gaussiana (GEM - Gaussian Expansion Method).

• Approccio Matematico: Il codice utilizza una serie di funzioni di base gaussiane per espandere le funzioni d'onda. Questo approccio permette di descrivere con precisione moderata ma efficiente l'asintotica delle funzioni d'onda, che è cruciale per i nuclei debolmente legati.
• Interazioni: Il codice è compatibile con varie funzioni per le interazioni efficaci a due nucleoni. Una caratteristica distintiva è la capacità di gestire l'interazione di Yukawa fino al canale tensoriale, nonché il termine di scambio di Coulomb e il termine a due corpi del moto del centro di massa senza approssimazioni aggiuntive.
• Simmetrie: I calcoli possono essere eseguiti sotto diverse ipotesi di simmetria applicate ai campi a un corpo:
  • Simmetria sferica (J) con conservazione della parità (P) e inversione temporale (T).
  • Simmetria assiale (rispetto all'asse z, $J_z$) con conservazione di P, T e simmetria di riflessione rispetto all'asse y (R).
  • Simmetria $J_z$ con conservazione di P e simmetria RT.
• Metodi Numerici: Per raggiungere l'autoconsistenza, sono disponibili due metodi:
  • Metodo Iterativo: L'Hamiltoniana HF o HFB viene diagonalizzata ad ogni iterazione.
  • Metodo del Gradiente Coniugato: La diagonalizzazione avviene solo alla fine del processo.
  • Per i calcoli HF+BCS, è possibile scegliere il metodo per la parte HF, mentre la parte BCS utilizza sempre il metodo iterativo.
• Tipologie di Calcolo: Il codice supporta calcoli non autoconsistenti (potenziale di Woods-Saxon o armonico), Hartree-Fock (HF), Hartree-Fock + BCS, e Hartree-Fock-Bogolyubov (HFB).

3. Contributi Chiave

I principali contributi tecnici del codice MFGSB sono:

1. Parametri di Base Universali: A differenza di altri metodi che richiedono la sintonizzazione dei parametri delle funzioni di base per ogni singolo nuclide o deformazione, i parametri del GEM sono insensibili a questi fattori. Un'unica serie di parametri è applicabile a un'ampia porzione della carta dei nuclidi, semplificando notevolmente l'uso del codice.
2. Gestione Completa delle Interazioni: È l'unico codice SCMF disponibile che gestisce l'interazione di Yukawa fino al canale tensoriale e tratta esattamente i termini di scambio di Coulomb e del moto del centro di massa.
3. Efficienza Computazionale: L'uso del GEM permette di coprire quasi tutti i nuclei con un singolo file di elementi di matrice a due corpi, riducendo la necessità di generare nuovi set di dati per ogni sistema.
4. Flessibilità: Supporto per vincoli sul momento di quadrupolo e scelta tra diversi metodi di convergenza.

4. Risultati e Specifiche Tecniche

• Disponibilità: Il codice (versione 1.0) è stato reso disponibile nel repository dell'Università di Chiba.
• Dimensioni e Struttura: Il pacchetto completo (mfgsb.zip) è di circa 16,6 GB, principalmente a causa dei dati di interazione. È disponibile una versione "leggera" (mfgsb light.zip) che separa i file di interazione per facilitare il download.
• Requisiti di Sistema:
  • La memoria richiesta dipende dal valore di taglio angolare ($\ell_{cut}$). Per il valore predefinito $\ell_{cut} = 7$, il codice richiede 8 GB di RAM. Per ogni incremento unitario di $\ell_{cut}$, la memoria richiesta raddoppia approssimativamente.
  • È necessario l'installazione delle librerie BLAS e LAPACK.
• Applicabilità: Il codice è stato validato per calcoli su un vasto spettro di nuclei, dimostrando la robustezza dei parametri di base universali.

5. Significato

MFGSB rappresenta un avanzamento significativo nella fisica nucleare teorica computazionale. La sua capacità di utilizzare un'unica serie di parametri per una vasta gamma di nuclei elimina l'incertezza e il tempo di calcolo associati alla sintonizzazione dei parametri di base. Inoltre, la trattazione rigorosa delle interazioni tensoriali e del moto del centro di massa senza approssimazioni aggiuntive lo rende uno strumento ideale per studi di precisione sulla struttura nucleare, specialmente per nuclei esotici o in regioni dove gli effetti di correlazione e tensoriali sono dominanti. La disponibilità pubblica del codice favorisce la riproducibilità e l'ulteriore sviluppo nella comunità scientifica.