CDCC 계산을 위한 고속 축소 기저 에뮬레이터 개발

초록

본 연구는 연속-이산화 결합채널(CDCC) 계산에 적용할 수 있는 축소 기저(Reduced Basis) 에뮬레이터를 제안한다. Proper Orthogonal Decomposition(POD)으로 얻은 스냅샷 기반 저차원 기저에 갈레르킨 투영을 적용해 18개의 광학 퍼텐셜 파라미터를 동시에 에뮬레이션한다. deuteron‑58Ni 21.6 MeV 산란을 테스트케이스로 사용했으며, 전산 시간은 100배 이상 단축되면서 탄성 산란 단면의 상대오차는 0.1 % 이하로 유지된다. 이 방법은 베이지안 추정과 불확실성 정량화에 필요한 대규모 샘플링을 실현 가능하게 만든다.

상세 분석

본 논문은 CDCC 계산이 갖는 두 가지 근본적인 병목을 정확히 짚어낸다. 첫째, 다중 부분파(Jmax≈100‑200)와 수십수백 개의 채널(Nch≈30‑50)을 포함하는 대규모 선형 시스템(Ntot≈5 000‑10 000)을 매 파동수마다 풀어야 한다는 점이다. 둘째, 광학 퍼텐셜 파라미터가 18차원에 달하는 고차원 파라미터 공간을 베이지안 샘플링(10⁴‑10⁶ 회)으로 탐색하려면 개별 CDCD 계산이 수분수시간 걸리는 현실적 한계가 존재한다. 저자는 이러한 문제를 “해결책은 파라미터가 변함에 따라 파동함수가 부드럽게 변한다는 물리적 직관”에 기반한다는 점을 강조한다. 즉, 몇 개의 대표적인 파라미터 집합(스냅샷)만 완전 계산으로 얻어두면, 그 주변의 해는 선형 결합으로 충분히 근사될 수 있다. 이를 구현하기 위해 Proper Orthogonal Decomposition, 즉 스냅샷 행렬에 대한 특이값 분해(SVD)를 수행한다. 특이값의 급격한 감소는 해 공간이 저차원 초평면에 거의 포함된다는 것을 의미하며, 에너지 기준(Eₙb>1‑εtol, εtol≈10⁻⁶)으로 5‑50개의 기저 벡터만 선택하면 전체 변동의 99.9999 % 이상을 포착한다.

선택된 기저 Xᵣ는 각 부분파 J마다 별도로 구축되며, 이는 J마다 채널 수와 매시 포인트가 다르기 때문이다. 온라인 단계에서는 새로운 파라미터 θ에 대해 전체 시스템 행렬 M(θ)를 그대로 구성하되, 이를 기저에 투사해 차원 n_b×n_b의 축소 시스템 Mʳα= bʳ을 푼다. 여기서 Galerkin 투사는 잔차를 기저와 직교하도록 강제함으로써, 원래 방정식의 물리적 구조를 보존하면서도 연산량을 크게 감소시킨다. 결과적으로, 원래 5 000‑10 000 차원의 선형 방정식을 20‑30 차원으로 축소해 수밀리초 수준의 평가 시간을 달성한다.

정밀도 검증에서는 deuteron‑58Ni 21.6 MeV 산란을 대상으로 18개의 광학 퍼텐셜 파라미터(깊이, 반경, 확산 등)를 동시에 변동시켰다. 에뮬레이터가 예측한 탄성 단면은 전 범위에서 상대오차가 0.1 % 이하이며, 특히 전형적인 실험 오차(≈1 %)보다 훨씬 작은 수준이다. 이는 POD 기반 기저가 파라미터 변동에 대한 비선형 효과까지 충분히 포착함을 의미한다. 또한, 속도 향상은 약 10² 배이며, 이는 베이지안 추정에 필요한 수십만 샘플을 실시간으로 생성할 수 있는 수준이다.

한계점으로는(1) 기저 구축 단계가 여전히 고비용이며, 파라미터 공간이 크게 확장될 경우 스냅샷 수가 증가할 수 있다. (2) 현재 구현은 선형 시스템에만 적용 가능하므로, 비선형 반응 메커니즘이나 다체 상호작용을 포함한 확장에는 추가 연구가 필요하다. 그럼에도 불구하고, CDCC와 같은 복잡한 결합채널 문제에 POD‑Galerkin 접근을 성공적으로 적용한 점은 향후 핵반응 이론 전반에 걸친 불확실성 정량화와 데이터 기반 파라미터 추정에 혁신적인 도구가 될 것이다.