방사선 전달 물리 불확실성 베이지안 추정

본 논문은 관성 구속 핵융합(I​CF) 실험에서 사용되는 방사선‑수소역학 시뮬레이션의 미세 물리 모델, 특히 방사선 전달에 대한 정확성을 정량적으로 평가하기 위해 베이지안 추론 방법을 제시한다. 저자들은 실험 데이터와 기존 이론(사전) 사이의 불일치를 해결하려면 수많은 물리 모델 파라미터와 실험적 잡음 파라미터를 동시에 고려해야 한다는 점을 강조한다. 이를 위해 ‘관심 파라미터’와 ‘잡음 파라미터’를 구분하고, 잡음 파라미터의 선형 응답을 매트릭스 A로 표현해 전체 확률 모델을 구축한다. 이 모델은 수정 χ² 함수(식 1)로 정리되며, 첫 번째 항은 전통적인 데이터‑모델 차이, 두 번째 항은 잡음 파라미터가 초래하는 정보 손실, 세 번째 항은 정규화, 네 번째 항은 사전 분포를 반영한다. 이렇게 함으로써 사전 지식과 실험 데이터를 동시에 고려한 최대 사후 확률(MAP) 추정이 가능해진다. 실험 대상은 NIF의 단일 샷(N110625)으로, 게르마늄 도핑 플라스틱 캡슐을 사용해 방사선 사진법으로 implosion 속도, 남은 어블레이터 질량, 그리고 연료 쉘의 ρR을 시간별로 측정하였다. 저자들은 다섯 개의 수정 인자(드라이브 강도, 드라이브 타이밍, Au M‑쉘 강도, C K‑엣지 흡수, Ge L‑엣지 흡수)를 정의하고, 각각에 대해 사전 평균값(대부분 1, 타이밍은 0 ns)과 사전 오차를 설정하였다. 잡음 파라미터는 캡슐의 물리적 치수·밀도·조성 29개를 1 % 오차로 가정하였다. 최적화는 유전 알고리즘을 사용해 수행되었다. 25세대, 각 세대 92개의 개체를 통해 약 2 300번의 HYDRA 1D 시뮬레이션을 실행했으며, 이는 200 CPU시간 이하로 완료되었다. 두 가지 경우(사전 포함 vs 사전 제외)를 비교함으로써 실험 데이터가 사전 지식에 비해 어느 정도의 정보를 제공하는지 정량화하였다. 사전을 포함한 경우 최적화된 수정 인자는 드라이브 강도 0.57±0.90, 타이밍 –0.45 ns±0.01 ns, Au M‑쉘 1.84±0.97, C K‑엣지 0.92±1.0, Ge L‑쉘 1.15±1.16으로 나타났다. 사전을 제외하면 모든 인자가 1에 가까운 값으로 수렴해, 실험 데이터만으로는 충분한 제약을 제공하지 못함을 보여준다. 결과 해석에 따르면, X‑ray 구동 강도와 Ge L‑쉘 흡수에 대한 기존 모델에 오차가 존재함을 확인했으며, 이는 사전 오차 범위 내에서 조정 가능하다. 반면, 방사선 전달 자체(예: 흡수 계수)의 수정만으로는 1D HYDRA 시뮬레이션과 실험 사이의 큰 차이를 해소할 수 없으며, 다차원 불안정성, 플라즈마 흐름, 혹은 기타 미세 물리 효과가 추가로 고려되어야 함을 시사한다. 또한, 베이지안 프레임워크와 GA 기반 최적화가 수십 개의 파라미터를 동시에 다루면서도 사전 지식과 실험 데이터를 일관되게 통합할 수 있는 강력한 도구임을 입증한다. 이 방법은 향후 ICF 실험 설계와 해석, 그리고 고에너지 밀도 물리 전반에 걸친 모델 검증에 활용될 수 있다.