[논문 리뷰] The effect of non-ionizing excitations on the diffusion of ion species and inter-track correlations in FLASH ultra-high dose rate radiotherapy

Ramin Abolfath, Alex Baikalov외 3명|arXiv (Cornell University)|2022. 01. 19.

Radiation Therapy and Dosimetry참고 문헌 26인용 수 26

한 줄 요약

이 논문은 플래시 초고선량률(UHDR) 방사선 치료에서 비이온화 상태의 진동(Non-ionizing excitations, NIS)이 양자화된 철학 포논을 통해 냉각된 이온 종류(IS)에 열 에너지를 전달하여 펌토세컨드 수준의 외향적 IS 폭발을 일으키는 미세구조적 메커니즘을 제안한다. 이는 확산 계수의 급격한 증가와 트랙 간 상관관계를 유도하며, 현재의 몬테카를로 코드가 트랙 간 반응을 과소평가하고 관측된 플래시 효과를 예측하지 못하는 이유를 설명한다.

ABSTRACT

We present a microscopic mechanism that accounts for the outward burst of "cold" ion species (IS) in a high-energy particle track due to coupling with ''hot" non-ion species (NIS). IS refers to radiolysis products of ionized molecules, whereas NIS refers to non-ionized excitations of molecules in a medium. The interaction is mediated by a quantized field of acoustic phonons, a channel that allows conversion of thermal energy of NIS to kinetic energy of IS, a flow of heat from the outer to the inner core of the track structure. We demonstrate the coexistence of "hot" NIS with "cold" IS in the radiation track structures right after their generation. NIS, concentrated within nano-scales volumes wrapping around IS, are the main source of intensive heat-waves and the outward burst of IS due to femto-second time scale IS-NIS coupling. By comparing the transport of IS coupled to NIS with identical configurations of non-interacting IS in thermal equilibrium at room temperature, we demonstrate that the energy gain of IS due to the surrounding hot nanoscopic volumes of NIS significantly increases their effective diffusion constants. The much higher diffusion constants predicted in the present model suggest higher inter-track chemical reaction rates at FLASH-UHDR, as well as lower intra-track reaction rates.

연구 동기 및 목표

• 플래시 UHDR 방사선 치료에서 증가된 트랙 간 상관관계의 물리적 기원을 설명하는 것.
• 현재 몬테카를로 코드가 관측된 생물학적 보존 효과를 예측하지 못하는 이유를 해결하는 것.
• 비이온화 상태의 진동(NIS)이 이온 종류(IS)의 확산 및 트랙 구조 역학에 미치는 영향을 정량화하는 것.
• NIS에서 IS로의 비가역적 열 흐름을 이끄는, 기본 원리에 기반한 메커니즘을 제공하는 것.
• 기존 몬테카를로 시뮬레이션 프레임워크에 IS-NIS 결합을 통합하여 정확한 플래시 모델링을 가능하게 하는 것을 제안하는 것.

제안 방법

• 아토세컨드 스케일에서 단계별 몬테카를로 시뮬레이션을 수행하여 IS 및 NIS의 공간 좌표와 에너지 침착량을 측정한다.
• 몬테카를로 트랙 시뮬레이션을 통해 온도 프로파일을 계산하고, 루더퍼드 산란 공식을 통해 분석적으로 검증한다.
• 온도 프로파일을 경계 조건으로 사용하여 다중 스케일 분자 동역학(MD) 시뮬레이션을 수행하여 급격한 트랙 팽창을 모델링한다.
• 기본 원리에 기반한 MD에서 유도된 확률적 군집화 랑주뱅 방정식을 통해 IS의 운동 방정식을 유도하며, 이는 포논을 통한 NIS에서 IS로의 비가역적 에너지 흐름을 기술한다.
• MD 시뮬레이션에서 고전적 레너드존 포텐셜을 사용하여 결과를 검증한다.
• 비평형 에너지 전달의 영향을 분리하기 위해 NIS 결합 유무에 따라 IS의 확산을 비교한다.

실험 결과

연구 질문

• RQ1 비이온화 상태의 진동(NIS)은 초고선량률에서 방사선 트랙 내 이온 종류(IS)의 확산에 어떻게 영향을 미치는가?
• RQ2 수소 방사선 분해에서 뜨거운 NIS에서 냉각된 IS로의 에너지 전달에 철학 포논이 수행하는 역할은 무엇인가?
• RQ3 왜 현재의 몬테카를로 코드는 플래시-UHDR 실험에서 트랙 간 상관관계를 예측하지 못하는가?
• RQ4 NIS 유도 에너지 전달은 실온에서의 확산과 비교해 IS의 효과적 확산 계수에 어떤 영향을 미치는가?
• RQ5 IS-NIS 결합은 플래시-UHDR에서 트랙 간 화학 반응 속도를 어느 정도 향상시키는가?

주요 결과

• 방사선 트랙 형성 직후에 뜨거운 NIS와 냉각된 IS의 공존이 확인되며, NIS는 IS를 감싸는 나노스케일 부피에 집중되어 있다.
• 비평형 NIS와 IS 간의 결합으로 인해 IS의 폭발적 운동이 펌토세컨드 수준에서 발생하며, 이는 실온에서의 값보다 훨씬 높은 효과적 확산 계수를 유도한다.
• 10 Gy 선량 및 40 Gy/s 선량률 조건에서, 모델은 현재의 몬테카를로 코드보다 더 큰 트랙 팽창과 강한 트랙 간 상관관계를 예측하며, 최근의 활성산소종 집합 및 비활성산소종(NROS) 형성 가설과 일치한다.
• NIS에서 IS로의 에너지 전달은 내부 열 흐름으로 작용하며, 포논을 통한 결합을 통해 IS의 비가역적 운동 에너지 증가를 이끈다.
• 이전 몬테카를로 시뮬레이션에서 트랙 간 겹침을 과소평가한 이유는 고정된 실온 확산 계수 때문임을 모델이 설명한다.
• 저자들은 모든 주요 몬테카를로 코드(Geant4, TOPAS-nBio, PARTRAC 등)에 IS-NIS 결합 및 내부 트랙 온도 프로파일을 통합하여 플래시-UHDR에서의 반응-확산을 정확히 모델링할 것을 권고한다.