Calibrated simulations for dynamic focusing of ultrasound through the temporal window

기술 요약: 시간적 창(Temporal Window)을 통한 초음파의 동적 집속을 위한 보정된 시뮬레이션

문제 정의
집속 초음파(FUS) 신경 조절은 치료적 임상 적용 및 인과적 뇌 연구를 위한 도구로 부상하고 있다. 동적 조향(dynamic steering)을 위해 다중 요소 트랜스듀서 어레이를 배치하는 데 있어 핵심적인 과제는, 음향 시뮬레이션과 측정된 초점 압력 사이의 불일치이다. 단일 초점 트랜스듀서는 종종 단일 보정값에 의존하지만, 3D 공간 전체를 조향할 수 있는 다중 요소 어레이는 기계적 또는 전기적 크로스토크(cross-talk) 및 두개골 투과에 대한 부정확한 가정으로 인해 복잡한 오류를 유발한다. 기존의 시뮬레이션 모델은 이러한 요인들을 반영하지 못해 예측된 압력 출력과 실제 출력 사이에 상당한 차이(20% 초과)를 보이는 경우가 많다. 또한, 두개골의 미세구조적 불균질성, 특히 골수(trabecular bone) 산란은 현재의 영상 해상도로 정확하게 시뮬레이션하기 어려우며, 이는 감쇠 손실의 과소평가로 이어진다. 본 논문은 개인 맞춤형 신경 조절 치료의 안전성과 정밀도를 보장하기 위해 시뮬레이션과 현실 사이의 간극을 메울 수 있는 효율적이고 정확한 보정 방법의 필요성을 다룬다.

방법론
저자들은 3D 조향 가능한 하위 공간에서 희소 실험 샘플링(sparse experimental sampling)과 최적화된 음향 시뮬레이션을 결합한 보정 프레임워크를 개발하였다.

• 대상 공간 정의: 157명의 피험자로부터 추출한 네 가지 양측 심부 뇌 구조(미상핵, 편도체, 시상 중심핵, 시상하핵)의 상대적 목표 좌표를 사용하였다. 이 좌표들은 잠재적 타겟의 95%를 포착하는 절단 원뿔형(truncated cone) 보정 공간을 정의하였다.
• 시뮬레이션-측정 보정:
  • 음향 홀로그래피: 소스 압력 매핑을 위해 홀로그래피를 사용하려는 초기 시도는 모델링되지 않은 요소 간 크로스토크로 인해 상당한 불일치(>20%)를 보였으며, 이에 따라 조향 공간 전체에 걸친 직접 측정이 필요하게 되었다.
  • 희소 샘플링: 저자들은 50개의 지점을 6mm 간격으로 등간격 샘ша플링함으로써 2D 보정 평면을 정확하게 특성화할 수 있음을 확인하였고, 이 불일치 데이터를 이방성 가우시안 함수(anisotropic Gaussian function, $R^2 = 0.83$)로 피팅하였다. 3D 공간의 경우, 베이지안 모델 선택(Bayesian model selection) 및 엘보우 탐지(elbow detection)를 통해 20개의 지점만으로도 불일치를 모델링하기에 충분하다는 것을 확인하였으며, 이 결과는 동일한 구성을 가진 5개의 트랜스듀서에 걸쳐 일반화되었다.
  • 축대칭성: 시뮬레이션은 축대칭적 장(field)을 시사했으나, 45° 단계로 8개 평면을 측정한 실험 결과 제조 변동성으로 인한 미세한 차이가 발견되었으며, 이는 비축대칭 어레이의 경우 완전한 3D 특성 파악이 필요함을 시사한다.
• 두개골 감쇠 최적화:
  • Ex Vivo 검증: 5개의 사체 두개골에서 추출된 시간적 창 파편(총 8개의 사용 가능한 파편)을 Zero Echo Time (ZTE) MRI로 스캔하였다. 이 파편들은 157명의 생존 피험자 MRI로부터 가상 추출된 데이터와 두께 및 입사각 측면에서 유사함을 확인하였다.
  • 감쇠 계수 도출: 각 파편에 대해 보정 공간 내 40개 위치를 대상으로 시뮬레이션 역전(time-reversal) 시뮬레이션을 수행하여, 균질한 뼈 감쇠 계수($\alpha_0$)를 반복적으로 수정하였다. 시뮬레이션된 피크 압력과 측정된 피크 압력 사이의 평균 제곱 오차를 최소화하는 값은 15.3 dB/cm/MHz로 결정되었다.
  • 오차 경계: 본 연구는 5번째 및 95번째 백분위수 압력 오차를 계산하였으며, 최적화된 계수를 사용할 때 실제 압력이 시뮬레이션 값보다 44% 낮거나 33% 높을 수 있음을 발견하였다.

주요 결과

• 샘플링 효율성: 본 연구는 보정을 위해 고해상도 장(field) 샘플링이 필수적인 것이 아니라, 2D 평면의 경우 50개 지점, 3D 볼륨의 경우 20개의 지점만으로도 이방성 가우시안 함수를 사용하여 시뮬레이션-현실 불일치를 모델링하기에 충분하다는 것을 입증하였다.
• 최적화된 감쇠: 시간적 창에 대해 최적화된 뼈 감쇠 계수인 15.3 dB/cm/MHz가 도출되었으며, 이는 이전의 실험적 추정치(1 MHz에서 13.3–16.6 dB/cm/MHz)와 일치한다.
• 불일치 특성: 시뮬레이션과 측정 사이의 불일치는 점진적이고 저주파적인 성분을 포함하고 있으며, 근접장(near-field) 측면 타겟 지점에서 가장 큰 불일치가 발생하는 것으로 나타났다.
• 일반화 가능성: 이 보정 접근 방식과 최적화된 감쇠 계수는 다양한 어레이 설계(동심원, 페르마 나선, 격자) 및 다수의 트랜스듀서 유닛에 적용 가능함을 보여주었다.

의의 및 주장
본 논문은 이 보정 프레임워크와 도출된 최적화된 시간적 창 계수가 개인 맞춤형 신경 조절을 위한 보고 및 안전성 평가의 정확도를 높이는 실질적인 해결책을 제공한다고 주장한다. 단일하고 보수적인 추정치나 보정되지 않은 시뮬레이션에서 벗어나, 저자들은 다음과 같은 데이터 기반 접근 방식을 제안한다:

1. 안전성 강화: 의도된 초점 위치뿐만 아니라 실제 피크 압력을 포착함으로써 기계적 지수(MI) 및 공간 피크 펄스 평균 강도($I_{SPPA}$)와 같은 파라미터에 대해 보수적인 안전 보정을 보장한다.
2. 정확도 향상: 최적화된 뼈 감쇠 계수는 시뮬레이션에서 흔히 과소평가되는 뼈의 미세구조에 의한 산란 손실을 반영한다.
3. 확장성 용이: 희소 샘플링 방법(20~50개 지점)은 현장에 배치된 개별 장치의 보정을 가능하게 하여, 모든 조향 범위의 지점을 특성화해야 하는 물류적 어려움을 해결한다.

저자들은 열적 유도(thermal derivatives) 또한 안전을 위해 중요하지만, 여기서 도출된 보정된 압력장은 흡수 계수를 비산란(non-scattering) 감쇠 분율을 나타내도록 조정함으로써 열 시뮬레이션(예: Pennes의 생체열 방정식)에 직접 통합될 수 있다고 강조한다. 본 연구는 결론적으로, 재현 노력을 안내하고 인체 적용을 위한 조향 가능한 초음파 어레이의 안전성을 보장하기 위해 이러한 감쇠 방법을 명확히 규정하는 것이 필수적이라고 밝히고 있다.