Validation of an Ab Initio-Informed Electronic Stopping Model for Large-Scale Atomistic Simulations

Quando particelle ad alta energia colpiscono un materiale solido, non rimbalzano semplicemente come palle da biliardo; esse squarciano il reticolo atomico, scagliando via gli atomi e innescando una caotica cascata di collisioni. Questo processo, noto come danno da radiazione, è il nemico silenzioso dei materiali nei reattori nucleari e nelle astronavi, dove l'intensa radiazione può gradualmente degradare l'integrità strutturale. Per prevedere quanto durerà un materiale, gli scienziati devono comprendere esattamente come queste particelle veloci perdano la loro energia mentre si fanno strada attraverso il solido. Per decenni, il modo standard per modellare questo fenomeno è stato quello di trattare gli atomi come particelle classiche che si muovono sotto l'influenza dei propri vicini, ignorando gli elettroni che orbitano attorno ad essi. Questa semplificazione funziona bene per molte cose, ma fallisce quando una particella si muove abbastanza velocemente da strappare energia a quegli elettroni, un fenomeno chiamato arresto elettronico. Senza un modo preciso per calcolare questa perdita di energia, le previsioni su quanto profondamente una particella penetrerà o quanto danno causerà rimangono inaffidabili, basandosi spesso su ipotesi arbitrarie piuttosto che sulla fisica fondamentale.

Un team di ricercatori dell'Università di Aalto, in Finlandia, ha compiuto ora un passo significativo verso la risoluzione di questo vuoto, testando un nuovo modello più sofisticato rispetto a dati reali. Si sono concentrati su un approccio specifico chiamato modello unificato a due temperature, che tenta di descrivere come l'energia fluisca tra gli atomi in movimento e la nuvola elettronica senza utilizzare numeri aggiustabili o inventati. Invece di indovinare quanta frizione avverta uno ione, questo modello deriva il proprio comportamento direttamente dalle leggi fondamentali della meccanica quantistica, in particolare osservando come cambia la densità elettronica al passaggio dello ione. I ricercatori volevano sapere se questo modello teoricamente perfetto potesse effettivamente prevedere ciò che accade in un esperimento reale, o se sarebbe crollato di fronte alla realtà disordinata di un materiale fisico. Per scoprirlo, hanno simulato un fascio di ioni di silicio, ciascuno dotato di 100 kiloelettronvolt di energia, che attraversa una sottile fetta di cristallo di silicio spessa solo 53 nanometri. Hanno poi confrontato i risultati generati dal computer con misurazioni effettuate in laboratorio, dove gli ioni passavano attraverso un foglio identico e venivano catturati da un rilevatore posto a 29 centimetri di distanza.

L'esperimento è stato progettato per essere incredibilmente sensibile ai dettagli della simulazione. Ruotando il cristallo di silicio, i ricercatori potevano costringere gli ioni a percorrere diverse traiettorie: alcuni scivolavano fluidamente lungo i canali aperti tra i file di atomi, mentre altri si scontravano casualmente con il reticolo atomico. Questa variazione è cruciale perché la quantità di energia persa da uno ione dipende fortemente dalla sua traiettoria; un percorso fluido attraverso i canali del cristallo causa una minore perdita di energia rispetto a un percorso caotico attraverso gli atomi. I modelli precedenti spesso trattavano tutti i percorsi allo stesso modo o utilizzavano regole semplici che cambiavano bruscamente in base alla velocità, mancando questi sottili effetti direzionali. Il nuovo modello, tuttavia, tiene conto del fatto che la nuvola elettronica risponde in modo diverso a seconda di dove si trova lo ione e di quanto velocemente si muove rispetto ai suoi vicini. Nella simulazione, i ricercatori hanno anche dovuto tenere conto di un sottile strato di ossido che si forma naturalmente sulla superficie del foglio di silicio. Hanno scoperto che l'inclusione di uno strato di ossido di 22 angstrom su entrambi i lati del foglio era essenziale per riprodurre i pattern di scattering osservati nell'esperimento reale, poiché questo strato sposta gentilmente gli ioni, disperdendoli quanto basta per corrispondere ai dati osservati.

Quando il team ha confrontato i risultati delle simulazioni con le misurazioni sperimentali, l'accordo è stato sorprendente. Il nuovo modello ha previsto con successo non solo la quantità media di energia persa dagli ioni, ma anche la specifica forma della distribuzione dell'energia. Per gli ioni che viaggiano dritti lungo i canali del cristallo, il modello ha previsto un intervallo ristretto di perdita di energia, mentre gli ioni che seguivano percorsi casuali mostravano una diffusione molto più ampia. Cosa più importante, la simulazione ha catturato una distinta struttura a due picchi nei dati per certi angoli, una caratteristica che deriva dalla miscela di ioni che rimanevano su percorsi fluidi e di quelli che venivano deviati. I vecchi modelli più semplici non riuscivano a riprodurre questa complessità, prevedendo spesso perdite di energia che erano o troppo alte o troppo basse, o mancando completamente della dispersione dei dati. I ricercatori hanno anche verificato le prestazioni del modello per ioni che si muovevano in direzioni casuali, dove la struttura cristallina non offre canali fluidi. Anche in questo caso, il modello ha retto, corrispondendo ai dati sperimentali senza la necessità di essere modificato o ricalibrato per l'angolo specifico del fascio.

Lo studio conferma che il modello unificato a due temperature è uno strumento affidabile per prevedere come le radiazioni danneggino i materiali. A differenza dei metodi precedenti che facevano affidamento su soglie arbitrarie o aggiustamenti empirici, questo approccio utilizza una descrizione dell'attrito elettronico che è radicata in calcoli basati sui primi principi, il che significa che è derivata dalle proprietà fondamentali del materiale stesso. I ricercatori hanno dimostrato che questo metodo può tracciare accuratamente la dissipazione di energia degli ioni che si muovono ad alte velocità attraverso un cristallo, catturando le sottili differenze tra percorsi fluidi e caotici. Questo livello di accuratezza è vitale per progettare materiali in grado di resistere agli ambienti ostili dello spazio e della tecnologia nucleare. Dimostrando che il modello funziona nel regime di arresto elettronico — la fase specifica in cui gli ioni perdono energia verso gli elettroni — il team ha fornito una solida base per future simulazioni su larga scala. Queste simulazioni possono ora essere utilizzate per prevedere il danno da radiazione con una fiducia precedentemente irraggiungibile, offrendo una via più chiara per comprendere come i materiali evolvono sotto il bombardamento incessante delle radiazioni.