[논문 리뷰] The effect of ion-induced shock waves on the transport of reacting species around energetic ion tracks

Pablo de Vera, Eugene Surdutovich외 3명|arXiv (Cornell University)|2017. 03. 27.

Ion-surface interactions and analysis참고 문헌 48인용 수 11

한 줄 요약

이 연구는 고에너지 이온 궤적 주변 수소 이온에 의해 유도된 충격파가 수소 과산화물 라디칼과 같은 반응성 종의 반경 방향 이동을 확산보다 80배 이상 향상시키는 것으로 나타내며, 반응성 종의 재결합을 억제하여 브라그 피크 영역에서 화학적 동역학을 크게 변화시킨다. 이는 방사선 치료에서 방사화학적 서명을 통한 간접적 실험적 탐지에 새로운 메커니즘을 제안한다.

ABSTRACT

The passage of energetic ions through tissue initiates a series of physico-chemical events, which lead to biodamage. The study of this scenario using a multiscale approach brought about the theoretical prediction of shock waves initiated by energy deposited within ion tracks. These waves are being explored in this letter in different aspects. The radial dose that sets their initial conditions is calculated using diffusion equations extended to include the effect of energetic $δ$-electrons. The resulting shock waves are simulated by means of reactive classical molecular dynamics. The simulations predict a characteristic distribution of reactive species which may have a significant contribution to biodamage, and also suggests experimental means to detect the shock waves.

연구 동기 및 목표

• 고에너지 이온 궤적 주변에서 이온에 의해 유도된 충격파가 반응성 종의 이동 및 화학적 성질에 미치는 영향을 조사하기 위해.
• 충격파가 브라그 피크 영역에서 OH 라디칼과 H₂O₂의 생존 및 분포에 미치는 영향을 평가하기 위해.
• 충격파가 방사선 손상의 화학적 단계에 미치는 영향, 특히 반응 경로 및 속도 상수의 변화를 평가하기 위해.
• 측정 가능한 반응성 종의 생성률(G-값) 변화를 통해 충격파를 실험적으로 추론할 수 있는 메커니즘을 제안하기 위해.
• 확산 방정정식에서 유도된 반경 방향 선량 프로파일을 반응성 MD 시뮬레이션과 융합하여 현실적인 충격파 동역학을 모델링하기 위해.

제안 방법

• 확장된 확산 방정식을 사용하여 브라그 피크 영역에서 저에너지 δ-전자 기여를 포함한 반경 방향 선량 프로파일을 계산하였다.
• MBN Explorer 소프트웨어를 사용하여 반응성 힘장 모델을 적용한 반응성 고전적 분자 동역학 시뮬레이션을 수행하여 화학 반응 및 충격파 전파를 모델링하였다.
• 이온 궤적에 의한 국소적 에너지 침착 기반으로 충격파를 모의하였으며, 압력파는 약 10⁻¹³ 초의 시간 스케일에서 형성되었다.
• 충격파 유무에 따라 라디칼 이동을 비교 분석하였으며, 시간에 따른 평균 제곱 이동 및 반응 수율을 분석하였다.
• 화학 동역학은 OH 라디칼 및 H₂O₂ 분자의 수를 모니터링하여 추적하였으며, 재결합 반응은 명시적으로 포함되었다.
• 확산 전용 케이스에 아인슈타인 관계를 적용하여 OH 확산 계수(D_OH = 0.166 Ų/ps)를 추출하였으며, 이는 시뮬레이션 설정의 타당성을 검증하였다.

실험 결과

연구 질문

• RQ1 이온에 의해 유도된 충격파는 브라그 피크 영역에서 OH 라디칼의 반경 방향 이동에 어떤 영향을 미치는가?
• RQ2 충격파는 OH 라디칼이 H₂O₂로 재결합하는 것을 어느 정도 억제하는가?
• RQ3 측정 가능한 반응성 종 생성률(G-값)의 변화를 통해 충격파 존재 여부를 실험적으로 추론할 수 있는가?
• RQ4 선형 에너지 전달(LET)에 따라 충격파 강도가 어떻게 변화하는가? 특히 브라그 피크 영역과 피크 외부 영역에서의 차이를 고려하여.
• RQ5 반응성 종의 이동에 있어 충격파에 의한 집단적 유동과 확산 중 어느 것이 더 기여하는가?

주요 결과

• 충격파는 확산 전용 경우에 비해 OH 라디칼의 평균 제곱 이동을 80배 이상 증가시켜 이동 효율성이 크게 향상됨을 시사한다.
• 확산 전용 케이스에서 유도된 OH 확산 계수(D_OH = 0.166 Ų/ps)는 기존의 알려진 값(예: GEANT4-DNA에서 0.28 Ų/ps)과 양호한 일치를 보였다.
• 16 ps 시뮬레이션 후, 충격파 존재 조건에서 생존한 OH 라디칼 수는 확산 전용 조건보다 40% 더 많았다.
• 충격파는 공간적 분리와 극한의 반응 조건으로 인해 H₂O₂ 형성을 억제하여 OH 재결합을 방지한다.
• 충격파 효과는 고에너지 침착이 일어나 강한 압력파가 형성되는 브라그 피크 영역에서 가장 강력하며, 이 영역 외부에서는 크게 약화된다.
• 이러한 결과는 충격파가 방사선 손상의 물리적 단계와 화학적 단계를 연결할 수 있음을 시사하며, 방사선 모델에서 이 두 단계의 시간적 분리 개념을 재고할 필요가 있음을 시사한다.