Fattore di struttura e intensità di diffrazione

Sapere dove appaiono i riflessi di diffrazione dipende dalla geometria del reticolo. Sapere quanto sono intensi dipende dagli atomi che stanno dentro la cella unitaria. Il legame tra i due è il fattore di struttura F(hkl): una somma di onde complesse, una per ogni atomo nella cella, che codifica sia le posizioni atomiche sia la capacità di ogni atomo di diffrangere i raggi X.

Comprendere il fattore di struttura è essenziale per capire le estinzioni dovute alla simmetria, il problema della fase e i metodi per risolvere strutture cristalline.

Il fattore di scattering atomico f_j

Ogni atomo diffrange i raggi X in proporzione alla sua nuvola elettronica: un atomo con più elettroni diffrange di più. Il fattore di scattering atomico f_j (o fattore di forma) è una funzione di sinθ/λ: vale il numero atomico Z per θ = 0 (diffusione coerente da tutti gli elettroni in fase) e cala rapidamente all’aumentare di θ perché le parti diverse della nuvola elettronica non sono più in fase tra loro. Questo calo è più rapido per gli atomi leggeri (piccola nuvola) che per quelli pesanti. Di conseguenza, gli atomi pesanti (Fe, Pb, U) dominano le intensità dei riflessi ad alto angolo, rendendo difficile localizzare gli atomi leggeri (H, Li) in loro presenza.

La formula del fattore di struttura

Per il riflesso con indici hkl, il fattore di struttura è la somma sugli n atomi della cella:

F(hkl) = Σ_j f_j exp[2πi(hx_j + ky_j + lz_j)]

dove f_j è il fattore di scattering atomico dell’atomo j e (x_j, y_j, z_j) sono le sue coordinate frazionali nella cella (compresi tra 0 e 1). L’esponenziale complesso rappresenta il contributo di fase dell’atomo j: la sua posizione determina di quanto la sua onda è sfasata rispetto all’origine. La somma di tutte queste onde sfasate dà l’ampiezza complessa totale F(hkl).

L’intensità misurata nel diffrattogramma è proporzionale al quadrato del modulo:

I(hkl) ∝ |F(hkl)|² = [Σ_j f_j cos(2π(hx_j+ky_j+lz_j))]² + [Σ_j f_j sin(2π(hx_j+ky_j+lz_j))]²

Estinzioni dovute al reticolo centrato

L’esempio più chiaro del potere del fattore di struttura è la spiegazione delle estinzioni sistematiche per reticoli centrati. Per un reticolo I (centrato al corpo), ogni atomo all’origine è accompagnato da un atomo identico a (1/2, 1/2, 1/2). Il contributo del secondo atomo ha fase 2π(h/2 + k/2 + l/2) = π(h+k+l). Quando h+k+l è dispari, questa fase vale π (rad) e le due onde si cancellano: F(hkl) = 0. Il riflesso è sistematicamente assente. Per reticolo F (fcc), gli atomi equivalenti per centratura aggiungono fasi che si cancellano a meno che h, k, l siano tutti pari o tutti dispari (stessa parità). Queste regole si derivano in modo rigoroso inserendo le posizioni degli atomi equivalenti nella formula del fattore di struttura.

Il problema della fase

Il diffrattografo misura l’intensità I(hkl) ∝ |F(hkl)|²: perde l’informazione sulla fase di F(hkl), cioè l’argomento del numero complesso. Conoscendo solo i moduli |F|, non è possibile calcolare direttamente la trasformata di Fourier inversa per ottenere la densità elettronica. Questo è il problema della fase (o phase problem) della cristallografia.

Tre famiglie di metodi lo affrontano: i metodi di Patterson (mappa della convoluzione delle densità, utile per strutture con pochi atomi pesanti), i metodi diretti (sfruttano la distribuzione statistica delle fasi per stimarle, efficaci per strutture con atomi simili come i composti organici) e il metodo dell’atomo pesante con raffinamento di Fourier. Il raffinamento Rietveld per polveri parte da un modello strutturale iniziale e minimizza iterativamente la differenza tra profilo calcolato e osservato, aggiustando sia le posizioni atomiche sia i fattori di temperatura.

Esempi di calcolo: Cu (fcc) e NaCl

Per il rame (fcc, atomi identici), F(hkl) = 4f_Cu se h,k,l stessa parità, altrimenti 0. Per NaCl (fcc con Na a (0,0,0) e Cl a (1/2,0,0) e siti equivalenti), F(hkl) = 4(f_Na ± f_Cl): il segno dipende da h+k+l. Se h+k+l pari, F = 4(f_Na + f_Cl); se dispari, F = 4(f_Na − f_Cl). Per KCl, f_K+ ≈ f_Cl− (stessa configurazione elettronica), quindi i riflessi con h+k+l dispari sono quasi nulli: il diffrattogramma di KCl sembra quello di un reticolo semplice (sc), non fcc — una delle curiosità storiche più citate in cristallografia.

Questo esempio storico (KCl che «sembra» cubico primitivo) mostrò quanto il contenuto atomico della cella — non solo la sua geometria — determini il diffrattogramma. È anche un monito: un’analisi superficiale degli indici di KCl potrebbe portare a concludere erroneamente che il reticolo è cubico semplice.

Domande frequenti

Cosa misura esattamente il diffrattometro?

Misura l’intensità dell’onda diffusa in ciascuna direzione: I(hkl) ∝ |F(hkl)|². Il quadrato del modulo del fattore di struttura contiene le informazioni sulle posizioni e sulla natura degli atomi, ma non la fase di F. È questa perdita della fase che rende il problema della determinazione strutturale non banale.

Perché gli atomi leggeri sono difficili da localizzare con i raggi X?

Perché il fattore di scattering f_j è proporzionale al numero di elettroni. L’idrogeno (1 elettrone) contribuisce meno dell’1% all’intensità totale in presenza di atomi pesanti. La diffrazione di neutroni risolve il problema perché i neutroni diffrangono sui nuclei (non sugli elettroni) e l’idrogeno ha lunghezza di scattering comparabile ai metalli pesanti (anche se con segno negativo, il che causa altre complicazioni).

Cos’è il raffinamento Rietveld?

È un metodo per determinare i parametri strutturali (coordinate atomiche, parametri di cella, fattori di temperatura) da un diffrattogramma su polveri, minimizzando la differenza tra il profilo di diffrazione calcolato e quello osservato. Non richiede monocristalli e è diventato il metodo standard per la caratterizzazione strutturale di materiali policristallini in campo accademico e industriale.

Perché le estinzioni sistematiche rivelano il gruppo spaziale?

Perché ogni elemento di simmetria con traslazione parziale (asse elicoidale, piano di scorrimento, reticolo centrato) produce un gruppo di riflessi con F = 0 esatto, per cancellazione delle fasi. Identificare quali classi di riflessi sono sistematicamente assenti permette di determinare quali elementi di simmetria sono presenti e quindi di identificare il gruppo spaziale.

Cosa sono i metodi diretti in cristallografia?

Sono algoritmi statistici per stimare le fasi di F(hkl) sfruttando il fatto che la densità elettronica deve essere positiva e localizzata sugli atomi. Le relazioni tra le fasi di riflessi diversi (relazioni di tripla, relazione di Sayre) permettono di determinare le fasi iterativamente a partire dai riflessi più intensi. I metodi diretti funzionano meglio quando gli atomi sono tutti simili (strutture organiche con soli C, N, O, H) e sono oggi implementati in tutti i software di cristallografia automatica.