최적 방사선 치료 계획을 위한 최소 엔트로피 모델

초록

시간 의존 볼츠만 전송 방정식을 기반으로 방사선 입자의 조직 내 상호작용을 모델링하고, 두 가지 목표(처방 선량 추적 및 선형‑이차 세포 사멸 모델)를 최소화하는 최적 제어 문제를 설정한다. 최소 엔트로피(M₁) 근사로 상태·대입 방정식을 해결하고, 최적성 조건을 도출한 뒤 2‑차원 수치 실험을 통해 방법의 효율성을 확인한다.

상세 분석

본 논문은 방사선 치료 계획을 수학적 최적 제어 문제로 정형화하고, 기존의 Fermi‑Eyges 기반 선량 계산이 갖는 비균질 조직에 대한 한계를 극복하기 위해 시간‑의존 볼츠만 전송 방정식을 채택하였다. 입자 밀도 ψ(t,x,Ω)는 위치 x와 방향 Ω에 대한 함수이며, 흡수·산란 단면 σₜ, σₛ와 Henyey‑Greenstein 산란 커널을 이용해 조직의 이방성 및 비균질성을 정확히 기술한다. 제어 변수는 공간·시간·방향에 의존하는 외부·내부 방사원 q(t,x,Ω)이며, 물리적 제한으로 q≥0 및 각 위치에서의 총 방사량이 상한 U(x) 이하가 되도록 제약한다.

두 가지 목적함수를 고려한다. 첫 번째는 처방 선량 D*(x)와 실제 총 선량 Dψ(x)=∫₀ᵀ∫_{S²}ψ dΩdt 사이의 L² 차이를 최소화하는 quadratic tracking functional이며, 가중치 c₁(x)·c₂(x)로 종양·위험 조직·정상 조직을 구분한다. 두 번째는 선형‑이차(LQ) 모델 S=exp(−αD−βD²)를 이용해 종양 세포의 생존률을 최소화하고, 동시에 정상 조직 손상을 억제하는 비선형 목적함수이다. LQ 파라미터 α,β는 조직별로 다르게 설정될 수 있어 실제 임상 데이터와 연계가 가능하다.

볼츠만 방정식은 고차원(시간·3D 위치·2D 방향) PDE로 직접 해석이 어려우므로, 최소 엔트로피(M₁) 근사를 도입한다. M₁ 폐쇄는 0차·1차 모멘트(에너지 밀도와 플럭스)만을 보존하면서도 비음향성, 비균질성, 비정상적인 방사선 흐름을 물리적으로 타당하게 재현한다. 이 근사는 연산량을 크게 감소시키면서도 양의성, 에너지 보존, 최대 전파 속도 제한 등 중요한 물리적 제약을 유지한다.

최적성 조건은 라그랑주 승수법과 연산자 이론을 이용해 도출된다. 상태 방정식 −Aψ+Σψ=Kψ+q와 그에 대한 대입 방정식 Aλ+Σλ=Kλ+∂J/∂ψ를 결합한 최적성 시스템을 얻으며, 제어는 투영 연산 proj_{