Calibrated simulations for dynamic focusing of ultrasound through the temporal window

Sintesi Tecnica: Simulazioni Calibrate per la Focalizzazione Dinamica degli Ultrasuoni attraverso la Finestra Temporale

Definizione del Problema
La neuromodulazione mediante ultrasuoni focalizzati (FUS) sta emergendo come uno strumento per applicazioni cliniche terapeutiche e per l'indagine causale del cervello. Una sfida critica nell'impiego di array di trasduttori multi-elemento per la sterzata dinamica è la discrepanza tra le simulazioni acustiche e le pressioni focali misurate. Spesso, i trasduttori a fuoco singolo si affidano a un singolo valore di calibrazione, mentre gli array multi-elemento, capaci di sterzare nello spazio 3D, introducono errori complessi dovuti al cross-talk meccanico o elettrico e ad assunzioni imprecise sulla trasmissione cranica. I modelli di simulazione esistenti spesso non riescono a tenere conto di questi fattori, portando a discrepanze significative (superiori al 20%) tra la pressione di uscita prevista e quella reale. Inoltre, l'eterogeneità microstrutturale del cranio, in particolare lo scattering dell'osso trabecolare, è difficile da simulare accuratamente con le attuali risoluzioni di imaging, portando a una sottostima delle perdite di attenuazione. Il documento affronta la necessità di metodi di calibrazione efficienti e accurati che colmino il divario tra simulazione e realtà per garantire la sicurezza e la precisione dei trattamenti di neuromodulazione personalizzati.

Metodologia
Gli autori hanno sviluppato un framework di calibrazione che combina il campionamento sperimentale sparso con simulazioni acustiche ottimizzate attraverso un sottospazio sterzabile 3D.

• Definizione dello Spazio Target: Le coordinate relative dei target per quattro strutture cerebrali profonde bilaterali (nucleo caudato, amigdala, nucleo centromediano del talamo e nucleo subtalamico) sono state estratte da 157 soggetti. Queste coordinate hanno definito uno spazio di calibrazione a cono tronco che cattura il 95% dei potenziali target.
• Calibrazione Simulazione-Misurazione:
  • Olografia Acustica: I tentativi iniziali di utilizzare l'olografia per la mappatura della pressione della sorgente hanno rivelato discrepanze significative (>20%) dovute al cross-talk non modellato degli elementi, rendendo necessaria la misurazione diretta attraverso lo spazio di sterzata.
  • Campionamento Sparso: Gli autori hanno determinato che un piano di calibrazione 2D può essere caratterizzato accuratamente campionando 50 punti con una spaziatura equidistante di 6 mm, adattando i dati di discrepanza con una funzione gaussiana anisotropa ($R^2 = 0,83$). Per lo spazio 3D, la selezione del modello bayesiano e il rilevamento del "gomito" (elbow detection) hanno indicato che 20 punti sono sufficienti per modellare la discrepanza, un risultato che si è generalizzato su cinque trasduttori con configurazioni identiche.
  • Assiometria: Sebbene le simulazioni suggerissero campi assimetrici, le misurazioni sperimentali di 8 piani a passi di 45° hanno rivelato modeste differenze dovute alla variabilità di produzione, suggerendo che la caratterizzazione 3D completa è necessaria per array non assimetrici.
• Ottimizzazione dell'Attenuazione Cranica:
  • Validazione Ex Vivo: Frammenti della finestra temporale sono stati estratti da 5 crani cadaverici (ottenendo 8 frammenti utilizzabili) e scansionati tramite MRI Zero Echo Time (ZTE). Questi frammenti sono stati validati rispetto a estrazioni virtuali dai MRI di 157 soggetti vivi, confermando la somiglianza in termini di spessore e angolo di incidenza.
  • Derivazione del Coefficiente di Attenuazione: Utilizzando simulazioni di time-reversal per mirare a 40 posizioni nello spazio di calibrazione per ogni frammento, gli autori hanno modificato iterativamente il coefficiente di attenuazione ossea omogeneo ($\alpha_0$). Il valore che minimizzava l'errore quadratico medio tra le pressioni di picco simulate e misurate è stato determinato in 15,3 dB/cm/MHz.
  • Confini di Errore: Lo studio ha calcolato gli errori di pressione al 5° e 95° percentile, riscontrando che le pressioni reali potevano essere inferiori del 44% o superiori del 33% rispetto ai valori simulati utilizzando il coefficiente ottimizzato.

Risultati Chiave

• Efficienza di Campionamento: Lo studio ha stabilito che non è richiesto un campionamento del campo ad alta risoluzione per la calibrazione; 50 punti per i piani 2D e 20 punti per i volumi 3D sono sufficienti per modellare la discrepanza tra simulazione e realtà utilizzando funzioni gaussiane anisotrope.
• Attenuazione Ottimizzata: È stato derivato un coefficiente di attenuazione ossea ottimizzato di 15,3 dB/cm/MHz per la finestra temporale, in linea con le stime sperimentali precedenti (13,3–16,6 dB/cm/MHz a 1 MHz).
• Caratterizzazione della Discrepanza: La discrepanza tra simulazione e misurazione è risultata contenere una componente graduale a bassa frequenza, con le maggiori discrepanze che si verificano nei punti target laterali in campo vicino.
• Generalizzabilità: L'approccio di calibrazione e il coefficiente di attenuazione ottimizzato si sono dimostrati applicabili a diversi design di array (cerchio concentrico, spirale di Fermat, griglia) e a molteplici unità di trasduttori.

Significato e Rivendicazioni
Il documento sostiene che questo framework di calibrazione e i coefficienti ottimizzati della finestra temporale derivati offrono una soluzione pratica per migliorare l'accuratezza della reportistica e delle valutazioni di sicurezza per la neuromodulazione personalizzata. Allontanandosi da stime singole e conservative o da simulazioni non calibrate, gli autori propongono un approccio basato sui dati che:

1. Migliora la Sicurezza: Catturare la pressione di picco effettiva anziché la posizione focale prevista garantisce calibrazioni di sicurezza conservative per parametri come l'Indice Meccanico (MI) e l'intensità di picco spaziale media pulsata ($I_{SPPA}$).
2. Migliora l'Accuratezza: Il coefficiente di attenuazione ottimizzato tiene conto delle perdite di scattering inerenti alla microstruttura ossea che vengono tipicamente sottostimate nelle simulazioni.
3. Facilita la Scalabilità: Il metodo di campionamento sparso (20–50 punti) rende la calibrazione fattibile per i singoli dispositivi distribuiti sul campo, affrontando le sfide logistiche della caratterizzazione di ogni punto in un intervallo di sterzata.

Gli autori sottolineano che, sebbene anche le derivate termiche siano critiche per la sicurezza, i campi di pressione calibrati derivati qui possono essere integrati direttamente nelle simulazioni termiche (ad esempio, l'equazione del calore biologico di Pennes) regolando il coefficiente di assorbimento per rappresentare la frazione di attenuazione non diffusa. Lo studio conclude che specificare chiaramente questi metodi di attenuazione è essenziale per guidare gli sforzi di replicazione e garantire la sicurezza degli array di ultrasuoni sterzabili nelle applicazioni umane.