Determinazione strutturale ai raggi X: dal cristallo alla struttura

Un cristallo racconta la propria struttura attraverso la luce che diffrange. La cristallografia ai raggi X a cristallo singolo è ancora oggi il metodo più potente per determinare con precisione atomica dove si trova ogni atomo in un reticolo: dalla posizione dei picchi di diffrazione si ricavano le distanze interplanari, dai loro moduli il fattore di struttura, e da quest’ultimo — superando il famoso problema della fase — la mappa della densità elettronica. Questo articolo percorre l’intera catena, dal cristallo montato sul goniometro al file CIF depositato in banca dati.

Dal cristallo alla raccolta dati: la geometria della diffrazione

Il processo comincia con un cristallo di qualità sufficiente (tipicamente lato ≥ 50 μm) montato su un nylon loop o un capillare e irradiato con raggi X monocromatici. La condizione affinché un insieme di piani reticolari (hkl) diffranga è la legge di Bragg:

nλ = 2 dhkl sinθ

dove λ è la lunghezza d’onda, d_hkl la distanza interplanare, θ il mezzo-angolo di diffrazione. In un diffrattometro moderno (a quattro cerchi o con detector 2D), il cristallo viene ruotato mentre il rivelatore raccoglie migliaia di riflessi. L’insieme completo di riflessioni — con posizioni (hkl) e intensità misurate I(hkl) — è il set di dati grezzo.

Un dettaglio raramente sottolineato nei testi di base: la posizione angolare di un riflesso identifica univocamente la famiglia (hkl) e quindi la distanza d_hkl, ma l’intensità non è affatto banale: dipende da quanti elettroni stanno sui piani e come sono disposti, dall’eventuale assorbimento del cristallo, da fattori geometrici (Lorentz-polarizzazione) e da fenomeni di estinzione. Prima ancora di passare alla struttura, i dati grezzi vanno corretti per questi effetti.

Il fattore di struttura F(hkl)

Dall’intensità misurata si ricava il modulo del fattore di struttura |F(hkl)| (dopo le correzioni): esso è la trasformata di Fourier della densità elettronica sulla cella elementare, e dipende dal tipo e dalla posizione di ogni atomo nella cella:

|F(hkl)| = Σj fj exp[2πi(hxj + kyj + lzj)]

Nella formula, f_j è il fattore di scattering atomico del j-esimo atomo (proporzionale al numero atomico Z per piccoli angoli), e x_j, y_j, z_j sono le sue coordinate frazionali. L’esperimento misura |F(hkl)|² ∝ I(hkl), ma non la fase φ(hkl). Ecco il problema della fase: senza la fase, non è possibile eseguire la trasformata di Fourier inversa per ricostruire la densità elettronica ρ(x,y,z) e quindi la struttura.

Metodi per risolvere il problema della fase

Tre approcci dominano:

• Metodi diretti: statistiche probabilistiche sulle fasi (relazioni di tripli, formula di Sayre) permettono di stimare le fasi di un sottoinsieme di riflessi forti e da qui di propagare la soluzione. Efficaci per strutture organiche medio-piccole (fino a ~1000 atomi in cella).
• Metodo di Patterson: la mappa di Patterson è la trasformata di Fourier di |F|²; contiene picchi in corrispondenza di tutti i vettori inter-atomici. Particolarmente utile quando la struttura contiene un atomo pesante (Z elevato) che domina lo scattering.
• Rimpiazzo isomorfo molecolare (MIR) / diffrazione anomala (SAD/MAD): usati in cristallografia proteica, dove i metodi diretti non convergono per strutture con decine di migliaia di atomi.

Un aspetto sottile che non tutti i manuali evidenziano chiaramente: anche dopo aver ottenuto un modello iniziale ragionevole, la struttura è ancora «approssimata». Occorre un ciclo iterativo di raffinamento per far convergere il modello sui dati sperimentali.

Il raffinamento e il fattore R

Il raffinamento di struttura regola posizioni atomiche, parametri di temperatura e altri parametri di scala per minimizzare la differenza tra i fattori di struttura osservati e quelli calcolati dal modello. La qualità del raffinamento è misurata dal fattore R:

R = Σ||Fobs| − |Fcalc||Σ|Fobs|

Un R del 5–7 % è considerato buono per una struttura organica a piccole molecole; sotto il 5 % è eccellente. Un R eccessivamente basso (es. 1 %) deve far sospettare un errore di modellazione o dati con pseudosimmetria. Affianco a R si usa il R_w (pesato) e il goodness of fit GooF, che idealmente vale 1.

Dalla struttura raffinata al file CIF

Il risultato finale è il CIF (Crystallographic Information File): un formato testo standardizzato dall’IUCr che contiene gruppo spaziale, parametri di cella, coordinate atomiche, fattori di temperatura e numerosi indicatori di qualità. I CIF vengono depositati nelle banche dati: CCDC per organici/metallo-organici, ICSD per inorganici, PDB per proteine. Il CIF è il «documento notarile» della struttura: chiunque nel mondo può scaricarlo, visualizzarlo e riprodurre il calcolo.

Limiti e casi speciali

La tecnica richiede un cristallo singolo di qualità: non è applicabile a polveri, materiali amorfi o campioni in soluzione (per questi esistono altri metodi). Cristalli con simmetria non meriedrica (geminati) complicano il raffinamento. Per strutture con disordine posizionale o molecole dinamicamente disorientate, il modello cristallografico è una media sul tempo e sullo spazio: la struttura «media» può non corrispondere a nessuna configurazione istantanea reale.

Parametri di qualità in cristallografia: riepilogo

Domande frequenti

Che cos’è il fattore di struttura F(hkl)?

È la grandezza complessa che descrive come tutti gli atomi della cella elementare contribuiscono, per un dato riflesso (hkl), all’onda diffratta. Il suo modulo è ricavabile dall’intensità misurata; la fase, invece, non è accessibile direttamente — da qui il problema della fase. La struttura cristallina si ricava solo invertendo F(hkl) nel dominio reale via trasformata di Fourier, ma ciò richiede sia modulo che fase.

Come si risolve il problema della fase?

I tre principali metodi sono: metodi diretti (sfruttano la positività della densità elettronica e relazioni statistiche tra le fasi), metodo di Patterson (mappa dei vettori inter-atomici, utile con atomi pesanti) e tecniche di diffrazione anomala (SAD/MAD), fondamentali in cristallografia macromolecolare. Una volta ottenuto un modello iniziale, il raffinamento iterativo corregge le imprecisioni.

Che cosa misura il fattore R?

Il fattore R (residual index) misura l’accordo tra i fattori di struttura osservati sperimentalmente e quelli calcolati dal modello: R = Σ||F_obs|−|F_calc||/Σ|F_obs|. Un R del 5–7 % è buono per piccole molecole. Valori molto più alti indicano un modello sbagliato o dati di scarsa qualità; valori sospettosamente bassi segnalano possibili errori di modellazione.

Cosa contiene un file CIF?

Il CIF (Crystallographic Information File) è il formato standard per comunicare una struttura cristallina: include il gruppo spaziale, i parametri di cella (a, b, c, α, β, γ), le coordinate frazionali di tutti gli atomi nella unità asimmetrica, i fattori di temperatura e tutti gli indicatori di qualità del raffinamento. Le banche dati internazionali (CCDC, ICSD, PDB) raccolgono centinaia di migliaia di CIF, liberamente scaricabili.

Quando serve la cristallografia ai raggi X in ambito industriale/normativo?

In campo farmaceutico, la forma cristallina (polimorfo) di un principio attivo influenza solubiltà, biodisponibilità e shelf-life: la normativa ICH Q6A richiede la caratterizzazione della forma solida. In chimica industriale e per i dossier REACH, la struttura cristallina è richiesta per identificare univocamente una sostanza, specialmente in caso di isomeri o polimorfismi. I pigmenti inorganici e alcuni additivi alimentari hanno forme cristalline normate.