Determination of nuclear quadrupole moments for $^{25}$Mg, $^{87}$Sr, and $^{135,137}$Ba via configuration-interaction combined with a coupled-cluster approach

Questo articolo impiega un approccio di interazione di configurazione più cluster accoppiato per calcolare i gradienti di campo elettrico e le costanti della struttura iperfine del momento di dipolo magnetico per gli stati a bassa energia di Mg, Sr e Ba, consentendo la determinazione accurata dei momenti quadrupolari nucleari per $^{25}$Mg, $^{87}$Sr e $^{135,137}$Ba che rivelano discrepanze significative con i valori precedentemente adottati per stronzio e bario.

L'essenziale

Immaginate il nucleo di un atomo non come una sfera perfetta, ma come un palloncino leggermente schiacciato o allungato. Questa forma non è casuale; è una "firma" specifica di come i protoni e i neutroni sono disposti all'interno. Gli scienziati chiamano questa forma momento quadrupolare nucleare. Conoscere l'esatta forma di questo "palloncino" è fondamentale per comprendere le regole fondamentali della fisica, dal modo in cui gli atomi si legano tra loro a come si comportano i materiali.

Tuttavia, misurare direttamente questa deformazione è incredibilmente difficile. È come cercare di indovinare la forma esatta di un palloncino dentro una scatola sigillata e buia, solo ascoltando il suono che emette quando lo si picchietta.

L'esperimento: Ascoltare l'"umore" atomico

In questo articolo, l'autore, Yong-Bo Tang, agisce come un maestro ingegnere del suono. Si concentra su tre specifiche "famiglie" di atomi: Magnesio (Mg), Stronzio (Sr) e Bario (Ba).

Quando questi atomi vengono eccitati (come una corda di chitarra pizzicata), emettono un "umore" o vibrazione molto specifica chiamata struttura iperfine. Questo umore è causato dall'interazione tra gli elettroni che ruotano all'esterno e il nucleo schiacciato all'interno.

• La parte misurata: Gli scienziati hanno già misurato l'altezza di questo umore in laboratorio con estrema precisione.
• L'anello mancante: Per capire la forma del nucleo (il momento quadrupolare) da quell'altezza, è necessario sapere esattamente come sono disposti gli elettroni attorno al nucleo. Questa disposizione crea un "gradiente del campo elettrico" (pensa alla pendenza di una collina lungo la quale gli elettroni rotolano verso il basso).

Il problema: La collina è troppo ripida per essere calcolata

Calcolare la forma di quella "collina" (il gradiente del campo elettrico) è un incubo per i computer. Gli elettroni non stanno solo fermi; danzano l'uno intorno all'altro, spingendosi e attirandosi in modi complessi chiamati correlazione elettronica.

• Se ignorate queste danze, il vostro calcolo della collina sarà errato.
• Se cercate di calcolare ogni singola danza, il vostro computer si blocca.

I tentativi precedenti di calcolare questo fenomeno erano come cercare di mappare una catena montuosa usando una foto satellitare sfocata. I risultati erano incoerenti. Per lo Stronzio e il Bario, studi diversi davano risposte diverse, con alcuni risultati che differivano fino al 10%.

La soluzione: Un approccio "Coltellino Svizzero" ibrido

Per risolvere questo problema, Tang ha sviluppato un nuovo metodo computazionale che combina due tecniche potenti:

1. Interazione di Configurazione (CI): Questa è come guardare ogni possibile modo in cui gli elettroni potrebbero disporsi, uno alla volta. È esaustiva ma lenta.
2. Cluster Accoppiato (CC): Questa è come usare una scorciatoia sofisticata per prevedere come gli elettroni si influenzano a vicenda in gruppi. È veloce ma a volte perde i dettagli minimi.

Il metodo di Tang, CI+CC, è il meglio dei due mondi. Utilizza la "scorciatoia" per gestire le grandi e pesanti interazioni tra gli elettroni del nucleo (core), e poi utilizza il metodo "esaustivo" per perfezionare i dettagli degli elettroni esterni. È come usare un drone per mappare la forma generale di una foresta, per poi inviare una squadra di escursionisti per misurare l'altezza esatta di ogni singolo albero.

I risultati: Chiarire la confusione

Utilizzando questo "coltellino svizzero" ad alta precisione, Tang ha calcolato i gradienti del campo elettrico per diversi stati a bassa energia di Mg, Sr e Ba. Ha poi combinato i suoi calcoli con gli "umori" sperimentali noti per determinare le forme nucleari.

Ecco cosa ha scoperto:

• Magnesio (25Mg): Il risultato è un match perfetto con gli esperimenti precedenti. È come sintonizzare una radio e trovare la stazione con un segnale cristallino. La forma calcolata concorda con quanto trovato utilizzando esperimenti di "raggi X muonici" (un altro modo, ad alta tecnologia, per misurare).
• Strontio (87Sr): Qui la trama si complica. Il risultato di Tang suggerisce che il nucleo sia circa il 10% più schiacciato rispetto al valore attualmente accettato nei libri di testo. Il vecchio valore derivava dall'osservazione di uno ione dello Stronzio (un atomo che ha perso un elettrone), mentre Tang ha osservato l'atomo neutro. La differenza suggerisce che il vecchio calcolo potrebbe aver trascurato alcune sottili danze elettroniche.
• Bario (135,137Ba): Similmente allo Strontio, i risultati di Tang per il Bario sono circa il 4% differenti dai valori attualmente accettati derivati dagli ioni del Bario.

La conclusione

L'articolo conclude che, sebbene il metodo funzioni magnificamente per il Magnesio, esiste una discrepanza significativa per lo Strontio e il Bario rispetto ai valori considerati "gold standard" attualmente utilizzati dagli scienziati.

Tang suggerisce che la differenza potrebbe derivare dal fatto che gli attuali calcoli "gold standard" hanno trascurato un tipo specifico di interazione elettronica chiamata eccitazione tripla (dove tre elettroni interagiscono simultaneamente). Proprio come un coro suona diversamente se tre cantanti armonizzano in un modo non previsto, queste interazioni triple potrebbero spostare l' "altezza" della forma dell'atomo.

In sintesi: l'autore ha costruito un modello computazionale migliore per misurare la forma dei nuclei atomici. Per il Magnesio, il modello ha confermato ciò che già sapevamo. Per lo Strontio e il Bario, il modello suggerisce che le misurazioni "ufficiali" attuali potrebbero essere leggermente errate, lasciando intendere che dobbiamo guardare più da vicino a come tre elettroni interagiscono per ottenere la vera forma di questi nuclei atomici.

Sommario tecnico

Sintesi Tecnica: Determinazione dei momenti quadrupoli nucleari per $^{25}\text{Mg}$, $^{87}\text{Sr}$ e $^{135,137}\text{Ba}$

Problematica
Il momento quadrupolo elettrico nucleare ($Q$) è un parametro fondamentale che descrive la simmetria della struttura nucleare, con significative implicazioni per la fisica atomica, nucleare e della materia condensata. Sebbene le costanti della struttura iperfine sperimentali ($B$) possano essere misurate con alta precisione, l'estrazione di valori accurati di $Q$ richiede calcoli teorici altrettanto precisi dei gradienti di campo elettrico ($q$). La letteratura esistente rivela discrepanze notevoli nei valori raccomandati per i momenti quadrupoli nucleari di $^{87}\text{Sr}$ e $^{137}\text{Ba}$, in particolare quando si confrontano i valori derivati da atomi neutri rispetto a quelli derivati da ioni. Ad esempio, i valori derivati dall'atomo neutro di Sr sono circa l'8% più grandi di quelli raccomandati dai dati degli ioni, e simili incongruenze esistono per il Ba. Inoltre, gli approcci teorici standard spesso faticano a tenere conto dei complessi effetti di correlazione elettronica necessari per calcoli del gradiente ad alta precisione in sistemi alcalino-terrosi pesanti.

Metodologia
Gli autori impiegano un approccio ibrido relativistico che combina il metodo del Configuration Interaction (CI) con il metodo del Coupled-Cluster (CC), denominato metodo RCI+all-order. Questo framework è progettato per tenere conto in modo esaustivo delle correlazioni elettroniche core-core, core-valence e valence-valence, che sono decisive per calcoli accurati del gradiente di campo elettrico.

Componenti metodologiche chiave includono:

• Hamiltoniana e Interazione: L'interazione iperfine è modellata utilizzando operatori tensoriali sferici per le interazioni di dipolo magnetico ($k=1$) e quadrupolo elettrico ($k=2$). Il nucleo è trattato come una sfera uniformemente magnetizzata e carica.
• Schema di calcolo ibrido: Il metodo utilizza un insieme finito di funzioni di tipo Gaussiano a temperatura uniforme per espandere le funzioni d'onda di Dirac. I potenziali di correlazione sono derivati tramite cinque distinte strategie computazionali per valutare l'incertezza:
  1. CI + Teoria delle perturbazioni molti corpi del secondo ordine (MBPT(2)).
  2. CI + Linearized Coupled-Cluster Singles and Doubles (LCCSD).
  3. CI + Full Coupled-Cluster Singles and Doubles (CCSD).
  4. CI + LCCSD con parametri di riscalamento (CI+LCCSDs).
  5. CI + CCSD con parametri di riscalamento (CI+CCSDs).
• Correzioni: Gli elementi di matrice di transizione includono l'Approssimazione del Campo Medio di Random Phase (RPA), effetti di Brueckner del core, radiazione strutturale e correzioni di normalizzazione (collettivamente denominate correzioni di ordine superiore o HO), nonché interazioni a due particelle (TP) fino al secondo ordine.
• Estrazione di Q: Il momento quadrupolo nucleare viene estratto utilizzando la relazione $Q = B / (234.9648867 \cdot q)$, dove $B$ è la costante della struttura iperfine misurata sperimentalmente e $q$ è il gradiente di campo elettrico calcolato teoricamente.

Contributi Chiave e Risultati
Lo studio calcola sistematicamente i gradienti di campo elettrico e le costanti della struttura iperfine di dipolo magnetico ($A$) per gli stati a bassa energia di $^{25}\text{Mg}$, $^{87}\text{Sr}$ e $^{135,137}\text{Ba}$ neutri. Gli stati specifici analizzati sono:

• Mg: $3s3p \ ^3P_{1,2}$
• Sr: $5s4d \ ^1D_2$ e $5s5p \ ^3P_{1,2}$
• Ba: $6s5d \ ^3D_{1,2,3}$, $6s5d \ ^1D_2$, e $6s6p \ ^1P_1$

Energia e Costanti Iperfine:
Le energie calcolate e le costanti di dipolo magnetico ($A$) mostrano un eccellente accordo con i dati sperimentali (NIST) e altri risultati teorici. Il metodo CI+LCCSD è identificato come il più affidabile, con i valori finali raccomandati per $A$ che corrispondono alle misurazioni sperimentali entro le loro incertezze (tipicamente entro il 2% per la maggior parte degli stati). Ciò valida l'approccio ibrido CI+CC per i sistemi alcalino-terrosi.

Momenti Quadrupoli Nucleari ($Q$):
Combinando i gradienti calcolati con i valori sperimentali di $B$, gli autori determinano i seguenti momenti quadrupoli nucleari:

• $^{25}\text{Mg}$: $0.203(2)$ b. Questo risultato è in perfetto accordo con l'esperimento di raggi X mesonici ($0.201(3)$ b) e con le precedenti determinazioni basate sull'iperfine.
• $^{87}\text{Sr}$: $0.336(4)$ b. Questo valore è circa il 10% più grande del valore attualmente adottato di $0.305(2)$ b (derivato dallo stato $4d_{5/2}$ di Sr$^+$), ma si allinea bene con i valori derivati dallo stato $5p_{3/2}$ di Sr$^+$ e dagli stati neutri di Sr ($5s5p \ ^3P_{1,2}$).
• $^{137}\text{Ba}$: $0.246(4)$ b. È circa il 4% più grande del valore raccomandato di $0.236(3)$ b (derivato dagli stati $5d_{3/2,5/2}$ di Ba$^+$). Tuttavia, concorda eccellentemente con i valori derivati dallo stato $6p_{3/2}$ di Ba$^+$.
• $^{135}\text{Ba}$: $0.161(2)$ b.

Significato e Rivendicazioni
L'articolo sostiene che l'approccio di interazione di configurazione più coupled-cluster fornisce un framework robusto e accurato per la determinazione dei momenti quadrupoli nucleari, particolarmente per nuclei instabili o pesanti dove la misurazione diretta è difficile.

Gli autori evidenziano una discrepanza significativa tra i valori derivati dagli atomi neutri e quelli dagli ioni per Sr e Ba. Nello specifico, i risultati per $^{87}\text{Sr}$ e $^{137}\text{Ba}$ derivati dagli atomi neutri in questo lavoro differiscono di circa il 10% e il 4%, rispettivamente, dai valori attualmente adottati in letteratura (che si basano su dati ionici). Gli autori suggeriscono che queste discrepanze possano derivare dall'omissione dei contributi delle eccitazioni triple nei calcoli standard di coupled-cluster, notando che lavori recenti sugli ioni indicano che le eccitazioni triple contribuiscono per l'1–2% alle costanti iperfini.

Lo studio conclude che, sebbene i risultati per $^{25}\text{Mg}$ siano coerenti con i dati stabiliti, le scoperte per $^{87}\text{Sr}$ e $^{137}\text{Ba}$ mettono in discussione i valori attualmente adottati derivati dagli ioni. Gli autori affermano che le loro estrazioni basate su atomi neutri sono affidabili, supportate dall'alta precisione dei calcoli delle costanti di dipolo magnetico, ma riconoscono che è necessario ulteriore lavoro teorico che includa le eccitazioni triple e le interazioni a tre corpi per risolvere pienamente le differenze tra le determinazioni basate su atomi neutri e ioni.