Calibrated simulations for dynamic focusing of ultrasound through the temporal window

Technische Samenvatting: Gekalibreerde Simulaties voor Dynamische Focussering van Ultrasone Golven door het Temporale Venster

Probleemstelling
Gefocuste ultrasone neuromodulatie (FUS) komt op als een instrument voor therapeutische klinische toepassingen en causale hersenonderzoeken. Een kritieke uitdaging bij het inzetten van multi-element transducer arrays voor dynamisch sturen is de discrepantie tussen akoestische simulaties en gemeten focale druk. Terwijl single-focus transducers vaak vertrouwen op een enkele kalibratiewaarde, introduceren multi-element arrays die over 3D-ruimte kunnen sturen complexe fouten door mechanische of elektrische cross-talk en onnauwkeurige aannames met betrekking tot schedeltransmissie. Bestaande simulatiemodellen slagen er vaak niet in om deze factoren te accounten, wat leidt tot significante mismatches (groter dan 20%) tussen de voorspelde en werkelijke drukoutputs. Bovendien is de microstructuur van de heterogeniteit van de schedel, met name de verstrooiing door trabeculair bot, moeilijk nauwkeurig te simuleren met de huidige beeldresoluties, wat leidt tot een onderschatting van de attenuatieverliezen. Het artikel adresseert de behoefte aan efficiënte, nauwkeurige kalibratiemethoden die de kloof tussen simulatie en realiteit overbruggen om de veiligheid en precisie van gepersonaliseerde neuromodulatiebehandelingen te waarborgen.

Methodologie
De auteurs ontwikkelden een kalibratiekader dat combineert met ijle experimentele bemonstering (sparse sampling) en geoptimaliseerde akoestische simulaties over een 3D-stuurbare subspace.

• Definitie van de Target Ruimte: Relatieve doelcoördinaten voor vier bilaterale diepe hersenstructuren (nucleus caudatus, amygdala, nucleus centromedianus van de thalamus en de nucleus subthalamica) werden geëxtraheerd uit 157 proefpersonen. Deze coördinaten definieerden een afgekort kegelvormig kalibratieruimte die 95% van de potentiële targets vastlegt.
• Simulatie-Metingskalibratie:
  • Akoestische Holografie: Initiële pogingen om holografie te gebruiken voor brondrukmapping toonden significante mismatches (>20%) door niet-gemodelleerde element cross-talk, wat de noodzaak voor directe meting over de stuurruimte vereiste.
  • IJle Bemonstering (Sparse Sampling): De auteurs stelden vast dat een 2D-kalibratievlak accuraat gekarakteriseerd kon worden door 50 punten te bemonsteren met een even afstand van 6 mm, waarbij de discrepantiegegevens werden gefit met een anisotrope Gaussische functie ($R^2 = 0,83$). Voor de 3D-ruimte gaven Bayesian model selectie en elbow detectie aan dat 20 punten voldoende waren om de discrepantie te modelleren, een bevinding die generaliseerde over vijf transducers met identieke configuraties.
  • Asymmetrie: Hoewel simulaties axiale symmetrische velden suggereerden, toonden experimentele metingen van 8 vlakken in stappen van 45° bescheiden verschillen door fabricagevariabiliteit, wat suggereert dat volledige 3D-karakterisering noodzakelijk is voor niet-axiale arrays.
• Optimalisatie van Schedelattenuatie:
  • Ex Vivo Validatie: Fragmenten van het temporale venster werden geëxtraheerd uit 5 cadaver schedels (waarbij 8 bruikbare fragmenten werden verkregen) en gescand via Zero Echo Time (ZTE) MRI. Deze fragmenten werden gevalideerd tegen virtuele extracties uit de MRI's van 157 levende proefpersonen, waarbij werd bevestigd dat de dikte en invalshoek vergelijkbaar waren.
  • Afleiding van de Attenuatiecoëfficiënt: Door middel van time-reversal simulaties gericht op 40 locaties in de kalibratieruimte voor elk fragment, pasten de auteurs iteratief de homogene botattenuatiecoëfficiënt ($\alpha_0$) aan. De waarde die de gemiddelde kwadratische fout tussen gesimuleerde en gemeten piekdrukken minimaliseerde, werd bepaald op 15,3 dB/cm/MHz.
  • Foutgrenzen: De studie berekende de 5e en 95e percentiel drukfouten, waarbij werd vastgesteld dat de werkelijke drukken 44% lager of 33% hoger konden zijn dan de gesimuleerde waarden bij gebruik van de geoptimaliseerde coëfficiënt.

Belangrijkste Resultaten

• Efficiëntie van Bemonstering: De studie stelde vast dat hoog-resolutie veldbemonstering niet vereist is voor kalibratie; 50 punten voor 2D-vlakken en 20 punten voor 3D-volumes bieden voldoende nauwkeurigheid om de simulatie-realiteit discrepantie te modelleren met anisotrope Gaussische functies.
• Geoptimaliseerde Attenuatie: Een geoptimaliseerde botattenuatiecoëfficiënt van 15,3 dB/cm/MHz werd afgeleid voor het temporale venster, wat overeenkomt met eerdere experimentele schattingen (13,3–16,6 dB/cm/MHz bij 1 MHz).
• Karakterisering van Discrepantie: De mismatch tussen simulatie en meting bleek een geleidelijk, laagfrequent component te bevatten, waarbij de grootste discrepanties optraden bij nabijveld laterale targetpunten.
• Generaliseerbaarheid: De kalibratiebenadering en de geoptimaliseerde attenuatiecoëfficiënt bleken toepasbaar op verschillende array-ontwerpen (concentrische cirkel, Fermat spiraal, rooster) en meerdere transducer-units.

Betekenis en Claims
Het artikel claimt dat dit kalibratiekader en de afgeleide geoptimaliseerde tempore venster coëfficiënten een praktische oplossing bieden om de nauwkeurigheid van rapportage en veiligheidsbeoordelingen voor gepersonaliseerde neuromodulatie te verbeteren. Door af te stappen van enkelvoudige, conservatieve schattingen of ongekalibreerde simulaties, stellen de auteurs een datagestuurde benadering voor die:

1. Veiligheid Verhoogt: Het vastleggen van de werkelijke piekdruk in plaats van de beoogde focale positie zorgt voor conservatieve veiligheidskalibraties voor parameters zoals de Mechanical Index (MI) en de spatial peak pulse average intensity ($I_{SPPA}$).
2. Nauwkeurigheid Verbetert: De geoptimaliseerde attenuatiecoëfficiënt houdt rekening met de verstrooiingsverliezen inherent aan de botmicrostructuur die in simulaties doorgaans worden onderschat.
3. Schaalbaarheid Faciliteert: De ijle bemonsteringsmethode (20–50 punten) maakt kalibratie haalbaar voor individuele apparaten die in het veld worden ingezet, waarmee de logistieke uitdagingen van het karakteriseren van elk punt in een stuurbereik worden aangepakt.

De auteurs benadrukken dat hoewel thermische derivaten ook cruciaal zijn voor de veiligheid, de gekalibreerde drukvelden die hier worden afgeleid, direct kunnen worden geïntegreerd in thermische simulaties (bijv. de Pennes' bioheat equation) door de absorptiecoëfficiënt aan te passen om het niet-verstrooiende deel van de attenuatie te representeren. De studie concludeert dat het duidelijk specificeren van deze attenuatiemethoden essentieel is voor het begeleiden van replicatiepogingen en het waarborgen van de veiligheid van stuurbare ultrasone arrays in menselijke toepassingen.