격자 위 원자핵 최초 계산을 위한 NuLattice

초록

NuLattice는 파이온‑없는 유효장 이론의 두‑·세‑입자 접촉 상호작용을 이용해 격자 기반 하트리‑포크, 전자‑상관(FCI), 결합군집(CC) 및 IMSRG 계산을 수행하는 파이썬 패키지이다. 짧은 거리 상호작용의 희소성을 활용해 노트북 수준의 컴퓨팅 자원으로 2H, 3,4He, 8Be, 12C, 16O까지의 바인딩 에너지를 얻는다.

상세 분석

본 논문은 전통적인 조화진동자 기반 핵 구조 계산이 갖는 비효율성을 격자 기반 단일입자 기저로 대체함으로써 해결책을 제시한다. 격자에서는 핵력의 짧은 거리 특성이 그대로 유지되어 두‑·세‑입자 상호작용의 비제로 행렬 원소가 전체 기저 차원 D에 비례하는 O(D) 수준으로 축소된다. 반면 조화진동자 기저에서는 삼입자 상호작용이 D⁵에 달하는 비현실적인 메모리를 요구한다. 저자들은 4 × 4 × 4 격자(L=4)에서 D=4L³=256을 사용해 4He의 해밀토니안 차원을 약 2.5 × 10⁷까지 줄였으며, 실제 메모리 요구량은 20 GB 이하로 노트북에서도 가능하도록 설계했다.

핵 Hamiltonian은 전이 연산자 a†, a 로 구성된 일‑·이‑·삼입자 항을 포함하고, 전자는 라플라시안 근사로 구현된 동역학 항과 접촉형 두‑·세‑입자 상수 v, w 로 정의된다. v와 w는 2H와 4He의 실험 결합 에너지에 맞춰 세 개의 격자 간격(a=2.5, 2.0, 1.7 fm)에서 보정되었다. 이 간단한 LO pion‑less EFT 모델은 SU(4) 대칭을 만족해 스핀·동등성 보존 및 두‑·세‑입자-두‑·세‑입자, 삼‑·세‑입자-삼‑·세‑입자 전이 항이 사라지는 특성을 가진다. 따라서 Hartree‑Fock 계산만으로도 대부분의 바인딩 에너지를 포착하고, CCSD에서 두‑입자‑두‑홀(2p‑2h) 진폭이 작아 높은 정확도를 유지한다.

알고리즘적으로는 모든 비제로 행렬 원소를 파이썬 리스트 형태로 저장하고, scipy.sparse의 CSR 포맷을 이용해 효율적인 행렬‑벡터 연산을 수행한다. Hartree‑Fock 단계에서는 밀도 행렬 ρ와 V, W 리스트를 순회하며 O(N_nonzero) 복잡도로 Fock 행렬을 구성한다. CCSD와 IMSRG는 일반적인 비선형 방정식 풀이(뉴턴‑유사 방법)와 BCH 전개를 사용하지만, 격자 기반 희소성 덕분에 텐서 축소 연산이 크게 간소화된다. 특히, 정상 순서화(two‑body truncation)가 정확히 성립하므로 세‑입자 상호작용을 두‑입자 형태로 완전하게 포함할 수 있다.

결과적으로 NuLattice는 2H, 3He, 4He, 8Be, 12C, 16O에 대해 FCI와 CCSD/IMSRG를 적용해 실험값과 비교했으며, 특히 12C와 16O에서 평균적인 오차는 1‑2 MeV 수준에 머문다. 이는 LO pion‑less EFT의 근본적인 제한(≈30 % 불확실성)보다 작은 차이이며, 격자 크기와 a를 조정하면 수렴성을 확인할 수 있다. 또한, 코드가 BSD‑3 라이선스로 공개돼 교육용 및 연구용으로 자유롭게 활용 가능하며, GPU 가속이나 MPI 병렬화를 추가하면 중간 규모 핵( A≈40)까지 확장 가능함을 시사한다.