분자 스캐터링 폼팩터 고속 계산을 위한 GPU 가속 프레임워크 SCarFFF

초록

본 논문은 TD‑DFT와 전이 밀도 행렬을 이용해 분자 수준의 다크 물질‑전자 산란 폼팩터를 3차원 전방향으로 계산하는 GPU 가속 코드 SCarFFF를 소개한다. B3YLP 교환‑상관 함수와 이중‑ζ 가우시안 기저함수를 사용해 10개의 무거운 원자를 가진 분자에 대해 12개의 전이 폼팩터를 약 5초 내에 얻을 수 있음을 시연한다. 반정밀 FFT 방식, 반정밀 반분석식 방식, 고정밀 전수치식 방식 등 세 가지 계산 경로를 제공해 정확도와 속도 사이의 트레이드오프를 자유롭게 선택한다.

상세 분석

이 연구는 다크 물질 탐지에 유망한 방향성 스칸틸링 분자결정의 전자 전이 폼팩터를 효율적으로 계산하는 방법론을 제시한다. 핵심은 전이 밀도 행렬(TDM)을 이용해 전자 전이의 공간적 구조를 직접적으로 얻고, 이를 푸리에 변환하여 동적 구조인자 S(q,ω)와 연결하는 것이다. 저자들은 TD‑DFT의 Casida 방정식을 풀어 X와 Y 행렬을 구하고, 이를 하나의 Tpq 행렬로 결합한다. Tpq는 원자 궤도 기반의 가우시안 함수에 직접 작용하므로, 별도의 파동함수 재구성 없이도 ⟨φ_a|e^{i q·r}|φ_i⟩ 형태의 전이 적분을 빠르게 계산할 수 있다.

코드 구현은 Julia 언어와 CUDA 기반 GPU 가속을 활용한다. 세 가지 계산 모드가 제공되는데, (1) 전통적인 3‑D FFT를 이용해 전이 전하밀도 Φ_{g→s}(r)를 격자화하고 푸리에 변환하는 방법은 메모리와 격자 크기에 민감하지만, 저정밀 요구 상황에서 가장 빠르게 결과를 산출한다. (2) 반정밀 반분석식 방법은 가우시안형 궤도에 대한 알려진 푸리에 변환식을 이용해 정확도를 유지하면서 연산량을 크게 줄인다. (3) 전수치식 방법은 격자와 FFT 없이 직접 적분을 수행해 가장 높은 정밀도를 제공한다.

성능 테스트에서는 10개의 무거운 원자를 포함한 유기 분자(예: p‑x릴렌)에 대해 12개의 전이 폼팩터를 5 초 이내에 계산했으며, 이는 기존 양자화학 패키지 대비 2~3 order of magnitude 빠른 속도다. 또한, 격자 크기와 박스 길이(L)의 선택이 폼팩터의 저주파 정확도에 미치는 영향을 상세히 분석해, 실용적인 ‘scale = 2’ 설정이 대부분의 다크 물질 모델(특히 무거운 매개자 모델)에서 충분히 정확함을 입증한다.

코드의 확장성도 강조한다. 현재는 스핀‑비의존 폼팩터만 구현했지만, 향후 스핀‑의존 상호작용을 위한 스핀 응답 함수도 추가될 예정이다. 또한, 다양한 기저함수(6‑31G*, cc‑pVDZ, cc‑pVTZ 등)와 교환‑상관 함수(B3YLP)를 지원해 사용자가 요구하는 정밀도와 계산 비용 사이를 자유롭게 조정할 수 있다. 이러한 특성은 수천 개 이상의 후보 물질을 고속으로 스크리닝하여 최적의 방향성 다크 물질 탐지 물질을 찾는 데 큰 도움이 될 것으로 기대된다.