🔬 physics

Assessment of the neutron fluence during operation of CEMHTI Cyclotron using Monte-Carlo codes - Simβ-AD study

Lo studio Simβ-AD ha stabilito una metodologia per la valutazione dell'attivazione dei rifiuti radioattivi presso il ciclotrone del CEMHTI, validando le misurazioni della fluenza neutronica da fogli di attivazione rispetto a quattro codici Monte-Carlo (FLUKA, MCNP6, PHITS e RayXpert®), dimostrando un buon accordo tra i dati sperimentali e le simulazioni utilizzando modelli semplificati.

Autori originali: Jean-Michel HORODYNSKI, Abir HASSANI, Frédéric CHAPELLE, Nicolas ARBOR, Djokhar BETELGUERIEV, Lucia Victoria GARCIA GARCIA, Stéphane HIGUERET, Léo NOWAK, Cédric DOSSAT, Sébastien BOUILLON, David CHAUL

Pubblicato 2026-07-07

📖 4 min di lettura☕ Lettura da pausa caffè

CC BY 4.0

Autori originali: Jean-Michel HORODYNSKI, Abir HASSANI, Frédéric CHAPELLE, Nicolas ARBOR, Djokhar BETELGUERIEV, Lucia Victoria GARCIA GARCIA, Stéphane HIGUERET, Léo NOWAK, Cédric DOSSAT, Sébastien BOUILLON, David CHAULIN, William HATE, Patrice RIFARD, Thierry SAUVAGE

Articolo originale sotto licenza CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ Questa è una spiegazione generata dall'IA dell'articolo qui sotto. Non è stata scritta né approvata dagli autori. Per precisione tecnica, consulta l'articolo originale. Leggi il disclaimer completo

La visione d'insieme: Prevedere il "fantasma" lasciato alle spalle

Immaginate di avere una macchina molto potente (un ciclotrone) che spara particelle minuscole come protoni e deuteroni contro un bersaglio. Quando queste particelle colpiscono, non si fermano semplicemente; creano una pioggia caotica di "fantasmi" invisibili chiamati neutroni. Questi fantasmi rimbalzano nella stanza, colpendo pareti, magneti e attrezzature, rendendoli leggermente radioattivi.

L'obiettivo principale di questo studio era capire: I nostri programmi di simulazione al computer possono prevedere con precisione quanto la stanza diventerà "radioattiva" a causa di questi fantasmi?

I ricercatori volevano costruire una ricetta affidabile (metodologia) per misurare questa radioattività senza dover smontare tutto in seguito. Questo è fondamentale per gestire i rifiuti radioattivi in modo sicuro.

I personaggi della storia

Il Ciclotrone: Un gigantesco acceleratore di particelle presso la struttura CEMHTI in Francia. È come una pista da corsa ad alta velocità per particelle minuscole. Recentemente ha cessato l'attività, rendendolo il luogo perfetto per testare le idee prima che la pista scompaia per sempre.
I "Fantasmi" (Neutroni): Particelle invisibili create quando il fascio colpisce il bersaglio. Sono i principali responsabili che rendono radioattivi altri materiali.
Le "Trappole Appiccicose" (Fogli di attivazione): I ricercatori hanno posizionato piccoli fogli sottili di metallo (Oro, Scandio, Tantalio e Terbio) in diversi punti della stanza. Pensateli come delle trappole adesive per mosche. Quando i fantasmi dei neutroni volano di passaggio, si attaccano al metallo, cambiandolo leggermente. In seguito, i ricercatori hanno misurato quanta "appiccicosità" (radioattività) c'era su ogni foglio.
Le "Sfere di Cristallo" (Codici Monte-Carlo): Questi sono quattro diversi programmi per supercomputer (FLUKA, MCNP6, PHITS e RayXpert). Sono come previsori del tempo. Voi inserite la mappa della stanza e la velocità delle particelle, e loro cercano di prevedere esattamente quanti fantasmi colpiranno le trappole adesive.

L'esperimento: Una verifica della realtà

Il team ha fatto funzionare il ciclotrone con due tipi di fasci (protoni e deuteroni) e ha posizionato le sue "trappole appiccicose" in varie posizioni: vicino all'uscita del fascio, vicino ai magneti e contro le pareti.

Dopo che la macchina è stata spenta, hanno portato i fogli metallici in un laboratorio e hanno misurato esattamente quanto fossero diventati radioattivi. Successivamente, hanno confrontato questi numeri del mondo reale con le previsioni fatte dai quattro programmi informatici.

I risultati: Come se sono comportate le Sfere di Cristallo?

I ricercatori hanno scoperto che i programmi per computer erano in realtà piuttosto bravi a indovinare il risultato, ma non perfetti.

Il punteggio: Le previsioni del computer erano solitamente comprese tra il 50% e il 150% delle misurazioni reali. Nel mondo della fisica complessa, questo è considerato un "buon accordo". È come un previsore del tempo che dice "pioverà tra 1 e 3 pollici" e in realtà piovono 2 pollici.
Il "glitch" dell'Oro: Uno dei programmi per computer (FLUKA) ha avuto un problema specifico con i fogli d'Oro. Sembrava sbagliare la quantità di radioattività perché non sapeva come gestire uno stato specifico degli atomi d'oro (come confondere un gatto che dorme con uno sveglio). Quando hanno corretto questa logica, le previsioni sono migliorate molto.
Il "problema dell'angolo": Vicino all'uscita del fascio (dove le particelle colpiscono per prime), le previsioni erano un po' più disperse. È come cercare di prevedere come l'acqua schizza quando si getta una pietra in uno stagno; le increspature vicino allo schizzo sono caotiche e difficili da modellare perfettamente. I computer a volte sottostimavano i "fantasmi" a bassa energia in questi angoli stretti.

Perché questo è importante (secondo l'articolo)

Lo studio dimostra che è possibile utilizzare questi codici informatici per stimare quanto rifiuto radioattivo produrrà un impianto a ciclotrone.

Il problema "Beta": Alcuni materiali radioattivi emettono solo un tipo di energia chiamata radiazione "beta", che è difficile da rilevare senza distruggere l'oggetto. L'articolo suggerisce che, se riusciamo a prevedere accuratamente l'attività dei neutroni usando questi codici, possiamo stimare il rifiuto "solo beta" senza dover tagliare i macchinari.
La conclusione: Combinando i modelli informatici con alcune misurazioni del mondo reale (le trappole appiccicose), possiamo creare un metodo affidabile per gestire i rifiuti radioattivi provenienti dagli acceleratori di particelle.

Riassunto in breve

I ricercatori hanno testato quattro diversi programmi per computer per vedere se potevano prevedere quanto una stanza sarebbe diventata radioattiva dopo il funzionamento di un acceleratore di particelle. Hanno usato veri fogli di metallo come "sensori" per catturare le particelle invisibili. I computer hanno fatto un lavoro discreto (entro un fattore di 1,5), dimostrando che possiamo usare queste simulazioni digitali per pianificare la gestione dei rifiuti radioattivi, a patto di conoscere esattamente come si comporta la macchina e correggere alcuni piccoli errori logici nel software.

Sintesi Tecnica: Valutazione della fluenza neutronica durante l'operatività del Ciclotrone CEMHTI utilizzando codici Monte-Carlo

Problematica
La gestione efficace dei rifiuti radioattivi per le strutture ciclotroniche richiede una caratterizzazione radiologica accurata, in particolare per i radionuclidi che emettono solo $\beta$ (ad es. $^{36}$ Cl, $^{49}$ V, $^{55}$ Fe, $^{63}$ Ni), i quali richiedono metodi di misurazione distruttivi. Lo studio Sim $\beta$ -AD è stato lanciato per stabilire una metodologia che combini i codici Monte-Carlo (MC) e rilevatori di neutroni attivi per valutare l'attivazione dei materiali all'interno delle strutture ciclotroniche. La sfida specifica affrontata in questo articolo è il benchmarking dei codici MC rispetto ai dati sperimentali per validare i calcoli della fluenza neutronica e definire il rapporto tra i radionuclidi che emettono solo $\beta$ e quelli che emettono $\gamma$ . Questo lavoro si concentra sulla struttura del Ciclotrone CEMHTI, che ha operato dal 1976 fino alla sua chiusura nel 2023, fornendo un'opportunità unica per validare i modelli rispetto a dati operativi reali prima dello smantellamento dell'impianto.

Metodologia
Lo studio ha impiegato un approccio comparativo tra misurazioni sperimentali e simulazioni numeriche utilizzando quattro distinti codici Monte-Carlo: FLUKA, MCNP6, PHITS e RayXpert.

Modellazione Geometrica e Fisica:
- È stato costruito un modello ad alta fedeltà del bunker CEMHTI, delle aree sperimentali, dei componenti del ciclotrone (magneti, giogo, iniettore) e dei target basandosi su progetti e misurazioni in situ.
- Sono stati modellati due scenari primari del fascio:
  1. Fascio di Protoni: 16 MeV, 25 µA fermi in un blocco di Tantalio, con una perdita secondaria sul deflettore.
  2. Fascio di Deuteroni: 17,5 MeV, 30 µA fermi in un target di Berillio.
- Le impostazioni fisiche specifiche sono state ottimizzate per ogni codice per migliorare la statistica della produzione di neutroni, includendo l'uso di librerie di sezioni d'urto specifiche (JEFF 3.3, TENDL-2019/2021/2024, ENDF/B-VIII.0) e tecniche di riduzione della varianza (ad es. collisioni forzate, campionamento per importanza).
Configurazione Sperimentale:
- Lamine di attivazione passiva composte da Oro (Au), Scandio (Sc), Tantalio (Ta) e Terbio (Tb) sono state posizionate in varie località (estrazione, tubo, giogo, deflettore, parete, porta, magnete) durante due campagne di irraggiamento nel 2021.
- Le lamine sono state irradiate per 5 ore (protoni) e 2 ore (deuteroni).
- Post-irraggiamento, è stata eseguita la $\gamma$ -spettroscopia utilizzando un detector ISOCS con efficienza del 40% per quantificare l'attività dei radionuclidi emettitori di $\gamma$ ( $^{198}$ Au, $^{46}$ Sc, $^{182}$ Ta, $^{160}$ Tb).
Calcoli di Attivazione:
- Le fluenze neutroniche derivate dai codici MC sono state utilizzate come input per calcoli di attivazione offline utilizzando FISPACT-II (con librerie TENDL-2019) e una pre-versione del codice RayActive.
- FLUKA ha inoltre eseguito calcoli di attivazione online utilizzando la sua card RESNUCLEI.

Risultati Chiave
Lo studio ha confrontato le attività misurate con i valori calcolati, ottenendo i seguenti risultati:

Accordo Complessivo: È stato raggiunto un buon accordo tra misurazioni e calcoli in tutti e quattro i codici, con rapporti tra attività misurata e calcolata compresi nell'intervallo 0,5–1,5.
Fascio di Protoni:
- I calcoli online di FLUKA hanno mostrato un buon accordo, fatta eccezione per le lamine d'Oro (Au), dove un rapporto arbitrario 50/50 tra stato fondamentale e isomero nella card RESNUCLEI ha portato a discrepanze. Ciò è stato risolto utilizzando FISPACT-II.
- I codici hanno mostrato maggiori discrepanze per le lamine situate vicino ai punti di perdita del fascio (deflettore, tubo), in particolare per Au, Ta e Tb, che sono sensibili ai neutroni a bassa energia. Ciò suggerisce variazioni nel calcolo della fluenza neutronica a bassa energia a piccoli angoli rispetto al fascio.
- RayXpert e MCNP6 hanno mostrato tendenze simili, mentre PHITS ha seguito la tendenza di FLUKA nella maggior parte dei casi.
Fascio di Deuteroni:
- I calcoli utilizzando FISPACT-II accoppiato con le fluenze di FLUKA, MCNP6 e PHITS hanno mostrato un buon accordo.
- I calcoli online di FLUKA hanno sottostimato significativamente l'attività della lamina d'Au a causa del problema dell'equilibrio dell'isomero.
- Altri codici tendevano a sovrastimare l'attivazione, probabilmente a causa di una sovrastima dei neutroni a bassa energia dalla reazione Be(d,n) e della successiva termalizzazione nelle pareti di cemento.
Performance dei Codici: RayXpert e MCNP6 hanno mostrato maggiori discrepanze rispetto agli altri codici vicino ai punti di perdita del fascio, attribuite alla sottostima della fluenza neutronica a bassa energia a piccoli angoli.

Significatività e Rivendicazioni
Il documento afferma che la stima della fluenza neutronica tramite codici Monte-Carlo, se abbinata a modelli di impianto relativamente accurati, fornisce un metodo affidabile per valutare l'attivazione dei rilevatori passivi in ambienti ciclotronici. Lo studio dimostra che:

Validazione della Metodologia: La combinazione di codici MC e strumenti di attivazione offline (FISPACT-II, RayActive) è una metodologia percorribile per caratterizzare l'attivazione dei rifiuti radioattivi, ottenendo risultati entro un intervallo di performance accettabile (rapporto 0,5–1,5).
Fonti di Incertezza: La principale fonte di incertezza per lo scenario del fascio di protoni è identificata nella mancanza di dati precisi riguardanti l'intensità e la forma della perdita del fascio all'interno dell'acceleratore.
Applicazione Futura: Con il ciclotrone CEMHTI ora spento, la metodologia validata sarà applicata per confrontare le simulazioni numeriche con le misurazioni effettuate su parti smantellate del ciclotrone e sulle pareti di cemento. Ciò faciliterà ulteriori studi sulle differenze specifiche tra i codici MC, le librerie di sezioni d'urto e i comportamenti di termalizzazione neutronica, supportando infine l'obiettivo di Sim $\beta$ -AD di definire i rapporti tra i radionuclidi che emettono solo $\beta$ e quelli che emettono $\gamma$ per la gestione dei rifiuti.

Gli autori concludono che, sebbene siano necessari miglioramenti nella conoscenza dei parametri del fascio, i risultati attuali confermano l'utilità di questi strumenti di simulazione per la caratterizzazione radiologica nelle strutture acceleratorie.