Geant4 기반 X선 광학 시뮬레이션 확장

초록

본 논문은 Geant4에 X선 반사와 굴절을 위한 물리 모델을 추가하여 일반적인 X선 트레이싱 툴로 활용할 수 있게 한 확장 모듈을 소개한다. 구현된 모듈을 XMM‑Newton의 Wolter‑I 거울 시스템에 적용해 유효 면적을 계산하고, 공개된 교정 데이터와 비교함으로써 정확성을 검증하였다.

상세 분석

Geant4는 입자 물리 시뮬레이션에 최적화된 프레임워크로, 기본 제공되는 전자·광자 상호작용 모델은 주로 고에너지 영역에 초점을 맞추고 있다. 따라서 X선 파장(0.1–10 keV)에서 발생하는 미세한 표면 반사, 거울 코팅의 복합 굴절률, 그리고 작은 입사각에서의 전반사 현상을 정확히 기술하기 위해서는 별도의 물리 프로세스가 필요하다. 저자들은 이러한 요구를 충족시키기 위해 Geant4의 G4VProcess 인터페이스를 상속한 “G4XrayReflectivity” 클래스를 구현하였다. 이 클래스는 실험적으로 측정된 재료의 복소 굴절률 데이터를 입력으로 받아, 입사각에 따른 반사율을 Fresnel 방정식과 표면 거칠기(σ) 보정(Guoy‑phase 포함)으로 계산한다. 또한, 다층 코팅 구조를 지원하도록 스택형 모델을 도입해 각 층마다 다른 n, k 값을 적용할 수 있게 하였다.

거울 표면의 미세 거칠기는 X선 반사 효율에 큰 영향을 미치므로, 저자들은 Gaussian 분포를 가정한 미세 구조 모델을 통해 스펙트럼 전반에 걸친 반사 손실을 정량화하였다. 구현된 프로세스는 Geant4의 표준 트래킹 루프에 자연스럽게 통합되며, 입사광선이 거울 표면에 도달하면 자동으로 반사/투과 판단이 이루어진다. 이와 동시에 입사광선의 파장, 입사각, 그리고 재료 특성을 실시간으로 조회해 물리량을 업데이트한다.

시뮬레이션 검증 단계에서는 XMM‑Newton의 Wolter‑I 설계(두 개의 연속된 포물선·쌍곡선 거울) 모델을 구축하였다. 거울의 기하학적 파라미터와 코팅 재료(주로 Au와 Ni) 정보를 정확히 입력하고, 광원은 무한거리 평면파로 가정해 다양한 에너지(0.5–10 keV)에서 입사시켰다. 결과적으로 유효 면적(Effective Area, EA)은 에너지에 따라 감소하는 경향을 보였으며, 이는 공식 교정 곡선과 3 % 이내의 차이만을 나타냈다. 특히, 저에너지 영역에서 표면 거칠기 보정이 없을 경우 EA가 과대평가되는 현상이 확인되어, 미세 거칠기 모델의 중요성을 입증하였다.

이 확장은 기존 Geant4 기반 방사선 방호·핵물리 시뮬레이션과 통합 가능하다는 장점이 있다. 따라서 복합 시스템(예: X선 검출기와 방사선 차폐가 동시에 존재하는 실험 장치)에서도 일관된 물리 모델링이 가능해진다. 다만, 현재 구현은 정적 거울 표면만을 다루며, 온도·스트레스에 의한 변형이나 동적 변형(진동 등)은 고려되지 않는다. 향후에는 이러한 동적 효과와 다중 반사 경로(예: 다중 거울 시스템)까지 확장할 여지가 있다.