Calibrated simulations for dynamic focusing of ultrasound through the temporal window

Resumen Técnico: Simulaciones Calibradas para el Enfoque Dinámico de Ultrasonido a Través de la Ventana Temporal

Planteamiento del Problema
La neuromodulación mediante ultrasonido focalizado (FUS) está emergiendo como una herramienta para aplicaciones clínicas terapéuticas e investigación causal del cerebro. Un desafío crítico al desplegar arreglos de transductores de múltiples elementos para el direccionamiento dinámico es la discrepancia entre las simulaciones acústicas y las presiones focales medidas. Mientras que los transductores de foco único a menudo dependen de un único valor de calibración, los arreglos de múltiples elementos capaces de direccionarse en el espacio 3D introducen errores complejos debido al cross-talk mecánico o eléctrico y a suposiciones inexactas respecto a la transmisión del cráneo. Los modelos de simulación existentes a menudo no logran contabilizar estos factores, lo que conduce a desajustes signific efectivos (superiores al 20%) entre la salida de presión predicha y la real. Además, la heterogeneidad microestructural del cráneo, particularmente la dispersión del hueso trabecular, es difícil de simular con precisión utilizando las resoluciones de imagen actuales, lo que conduce a una subestimación de las pérdidas por atenuación. El artículo aborda la necesidad de métodos de calibración eficientes y precisos que cierren la brecha entre la simulación y la realidad para garantizar la seguridad y la precisión de los tratamientos de neuromodulación personalizados.

Metodología
Los autores desarrollaron un marco de calibración que combina el muestreo experimental disperso con simulaciones acústicas optimizadas a través de un subespacio direccionable en 3D.

• Definición del Espacio Objetivo: Las coordenadas relativas de los objetivos para cuatro estructuras cerebrales profundas bilaterales (núcleo caudado, amígdala, núcleo centromediano del tálamo y núcleo subtalámico) fueron extraídas de 157 sujetos. Estas coordenadas definieron un espacio de calibración de cono truncado que captura el 95% de los objetivos potenciales.
• Calibración Simulación-Medición:
  • Holografía Acústica: Los intentos iniciales de utilizar la holografía para el mapeo de la presión de la fuente revelaron desajustes significativos (>20%) debido al cross-talk de los elementos no modelado, lo que hizo necesario el la medición directa a través del espacio de direccionamiento.
  • Muestreo Disperso: Los autores determinaron que un plano de calibración 2D podía caracterizarse con precisión mediante el muestreo de 50 puntos con un espaciamiento equidistante de 6 mm, ajustando los datos de discrepancia con una función gaussiana anisotrópica ($R^2 = 0.83$). Para el espacio 3D, la selección de modelos bayesianos y la detección del codo indicaron que 20 puntos eran suficientes para modelar la discrepancia, un hallazgo que se generalizó a cinco transductores con configuraciones idénticas.
  • Acimetría: Aunque las simulaciones sugerían campos axisimétricos, las mediciones experimentales de 8 planos a pasos de 45° revelaron diferencias modestas debido a la variabilidad de fabricación, lo que sugiere que es necesaria una caracterización 3D completa para arreglos no aximétricos.
• Optimización de la Atenuación del Cráneo:
  • Validación Ex Vivo: Se extrajeron fragmentos de la ventana temporal de 5 cráneos de cadáveres (obteniendo 8 fragmentos utilizables) y se escanearon mediante resonancia magnética de Tiempo de Eco Cero (ZTE). Estos fragmentos fueron validados contra extracciones virtuales de los MRIs de 157 sujetos vivos, confirmando la similitud en el espesor y el ángulo de incidencia.
  • Derivación del Coeficiente de Atenuación: Utilizando simulaciones de reversión de tiempo para apuntar a 40 ubicaciones dentro del espacio de calibración para cada fragmento, los autores modificaron iterativamente el coeficiente de atenuación ósea homogéneo ($\alpha_0$). El valor que minimizó el error cuadrático medio entre las presiones pico simuladas y las medidas fue de 15.3 dB/cm/MHz.
  • Límites de Error: El estudio calculó los percentiles 5 y 95 de los errores de presión, encontrando que las presiones reales podían ser un 44% menores o un 33% mayores que los valores simulados utilizando el coeficiente optimizado.

Resultos Clave

• Eficiencia de Muestreo: El estudio estableció que no se requiere un muestreo de campo de alta resolución para la calibración; 50 puntos para planos 2D y 20 puntos para volúmenes 3D son suficientes para modelar la discrepancia entre simulación y realidad utilizando funciones gaussianas anisotrópicas.
• Atenuación Optimizada: Se derivó un coeficiente de atenuación ósea optimizado de 15.3 dB/cm/MHz, alineándose con estimaciones experimentales previas (13.3–16.6 dB/cm/MHz a 1 MHz).
• Caracterización de la Discrepancia: Se encontró que el desajuste entre la simulación y la medición contiene un componente gradual de baja frecuencia, con las mayores discrepancias ocurriendo en los puntos laterales de la región cercana (near-field).
• Generalizabilidad: El enfoque de calibración y el coeficiente de atenuación optimizado mostraron ser aplicables a diferentes diseños de arreglos (círculo concéntrico, espiral de Fermat, rejilla) y a múltiples unidades de transductores.

Significancia y Reivindicaciones
El artículo afirma que este marco de calibración y los coeficientes optimizados de la ventana temporal ofrecen una solución práctica para mejorar la precisión de los informes y las evaluaciones de seguridad para la neuromodulación personalizada. Al alejarse de las estimaciones únicas y conservadoras o de las simulaciones no calibradas, los autores proponen un enfoque basado en datos que:

1. Mejora la Seguridad: Capturar la presión pico real en lugar de la posición focal prevista asegura calibraciones de seguridad conservadoras para parámetros como el Índice Mecánico (MI) y la intensidad de pico espacial promedio de pulso ($I_{SPPA}$).
2. Mejora la Precisión: El coeficiente de atenuación optimizado tiene en cuenta las pérdidas por dispersión inherentes a la microestructura ósea que suelen subestimarse en las simulaciones.
3. Facilita la Escalabilidad: El método de muestreo disperso (20–50 puntos) hace que la calibración sea factible para dispositivos individuales desplegados en el campo, abordando los desafíos logísticos de caracterizar cada punto en un rango de direccionamiento.

Los autores enfatizan que, si bien los derivados térmicos también son críticos para la seguridad, los campos de presión calibrados derivados aquí pueden integrarse directamente en simulaciones térmicas (por ejemplo, la ecuación de bioheat de Pennes) ajustando el coeficiente de absorción para representar la fracción de atenuación no dispersiva. El estudio concluye que especificar claramente estos métodos de atenuación es esencial para guiar los esfuerzos de replicación y garantizar la seguridad de los arreglos de ultrasonido direccionables en aplicaciones humanas.