Calibrated simulations for dynamic focusing of ultrasound through the temporal window

Resumo Técnico: Simulações Calibradas para o Foco Dinâmico de Ultrassom através da Janela Temporal

Definição do Problema
A neuromodulação por ultrassom focalizado (FUS) está emergindo como uma ferramenta para aplicações clínicas terapêuticas e investigação causal do cérebro. Um desafio crítico na implementação de matrizes de transdutores de múltiplos elementos para direcionamento dinâmico é a discrepância entre as simulações acústicas e as pressões focais medidas. Enquanto transdutores de foco único frequentemente dependem de um único valor de calibração, matrizes de múltiplos elementos capazes de direcionamento em 3D introduzem erros complexos devido ao cross-talk mecânico ou elétrico e suposições imprecisas sobre a transmissão do crânio. Os modelos de simulação existentes frequentemente falham em considerar esses fatores, levando a incompatibilidades significativas (superiores a 20%) entre os resultados de pressão previstos e reais. Além disso, a heterogeneidade microestrutural do crânio, particularmente a dispersão do osso trabecular, é difícil de simular com precisão usando as resoluções de imagem atuais, levando à subestimação das perdas por atenuação. O artigo aborda a necessidade de métodos de calibração eficientes e precisos que aproximem a lacuna entre a simulação e a realidade para garantir a segurança e a precisão de tratamentos de neuromodulação personalizados.

Metodologia
Os autores desenvolveram uma estrutura de calibração combinando amostragem experimental esparsa com simulações acústicas otimizadas através de um subespaço direcionável em 3D.

• Definição do Espaço Alvo: Coordenadas relativas de alvos para quatro estruturas cerebrais profundas bilaterais (caudado, amígdala, núcleo centromediando do tálamo e núcleo subtalâmico) foram extraídas de 157 sujeitos. Essas coordenadas definiram um espaço de calibração de cone truncado capturando 95% dos potenciais alvos.
• Calibração Simulação-Medição:
  • Holografia Acústica: Tentativas iniciais de usar holografia para mapeamento de pressão de fonte revelaram discrepâncias significativas (>20%) devido ao cross-talk de elementos não modelado, necessitando de medição direta através do espaço de direcionamento.
  • Amostragem Esparsa: Os autores determinaram que um plano de calibração 2D pode ser caracterizado com precisão amostrando 50 pontos com um espaçamento equidistante de 6 mm, ajustando os dados de discrepância com uma função gaussiana anisotrópica ($R^2 = 0,83$). Para o espaço 3D, a seleção de modelo bayesiano e a detecção de cotovelo (elbow detection) indicaram que 20 pontos eram suficientes para modelar a discrepância, um achado que se generalizou para cinco transdutores com configurações idênticas.
  • Axisimetria: Embora as simulações tenham sugerido campos axisimétricos, as medições experimentais de 8 planos em passos de 45° revelaram diferenças modestas devido à variabilidade de fabricação, sugerindo que a caracterização 3D completa é necessária para matrizes não axisimétricas.
• Otimização da Atenuação do Crânio:
  • Validação Ex Vivo: Fragmentos da janela temporal foram extraídos de 5 crânios de cadáveres (resultando em 8 fragmentos utilizáveis) e escaneados via Ressonância Magnética de Tempo de Eco Zero (ZTE). Esses fragmentos foram validados contra extrações virtuais de MRIs de 157 sujeitos vivos, confirmando a semelhança em espessura e ângulo de incidência.
  • Derivação do Coeficiente de Atenuação: Usando simulações de reversão de tempo para atingir 40 localizações no espaço de calibração para cada fragmento, os autores modificaram iterativamente o coeficiente de atenuação óssea homogêneo ($\alpha_0$). O valor que minimizou o erro quadrático médio entre as pressões de pico simuladas e medidas foi determinado como 15,3 dB/cm/MHz.
  • Limites de Erro: O estudo calculou os erros de pressão do 5º e 95º percentil, encontrando que as pressões reais poderiam ser 44% menores ou 33% maiores que os valores simulados usando o coeficiente otimizado.

Resultos Principais

• Eficiência de Amostragem: O estudo estabeleceu que não é necessária uma amostragem de campo de alta resolução para calibração; 50 pontos para planos 2D e 20 pontos para volumes 3D fornecem precisão suficiente para modelar a discrepância entre simulação e realidade usando funções gaussianas anisotrópicas.
• Atenuação Otimizada: Um coeficiente de atenuação óssea otimizado de 15,3 dB/cm/MHz foi derivado, alinhando-se com estimativas experimentais anteriores (13,3–16,6 dB/cm/MHz a 1 MHz).
• Caracterização da Discrepância: A incompatibilidade entre simulação e medição revelou conter um componente gradual de baixa frequência, com as maiores discrepâncias ocorrendo em pontos alvo laterais de campo próximo.
• Generalizabilidade: A abordagem de calibração e o coeficiente de atenuação otimizado mostraram-se aplicáveis a diferentes designs de matriz (círculo concêntrico, espiral de Fermat, grade) e múltiplas unidades de transdutor.

Significância e Alegações
O artigo afirma que esta estrutura de calibração e os coeficientes otimizados da janela temporal derivados oferecem uma solução prática para melhorar a precisar do relato e das avaliações de segurança para neuromodulação personalizada. Ao se afastar de estimativas únicas e conservadoras ou de simulações não calibradas, os autores propõem uma abordagem baseada em dados que:

1. Aumenta a Segurança: Capturar a pressão de pico real, em vez da posição focal pretendida, garante calibrações de segurança conservadoras para parâmetros como o Índice Mecânico (MI) e a intensidade de pico espacial média de pulso ($I_{SPPA}$).
2. Melhora a Precisão: O coeficiente de atenuação otimizado contabiliza as perdas por dispersão inerentes à microestrutura óssea que são tipicamente subestimadas em simulações.
3. Facilita a Escalabilidade: O método de amostragem esparsa (20–50 pontos) torna a calibração viável para dispositivos individuais implantados em campo, abordando os desafios logísticos de caracterizar cada ponto em um intervalo de direcionamento.

Os autores enfatizam que, embora as derivadas térmicas também sejam críticas para a segurança, os campos de pressão calibrados derivados aqui podem ser diretamente integrados em simulações térmicas (por exemplo, equação de bioaquecimento de Pennes) ao ajustar o coeficiente de absorção para representar a fração de atenuação não dispersiva. O estudo conclui que especificar claramente esses métodos de atenuação é essencial para orientar esforços de replicação e garantir a segurança de matrizes de ultrassom direcionáveis em aplicações humanas.