고정밀 분자 상대론 계산으로 구한 87Sr 핵 전기 사중극자 모멘트

초록

본 연구는 SrO와 SrS의 실험값인 핵 사중극자 결합 상수(NQCC)를 이용하고, X2C‑mmf 기반의 고차 전자 상관 방법(CCSD(T), CCSD‑T, CCSD˜T)으로 계산한 전기장 구배(EFG)를 결합하여 87Sr의 핵 전기 사중극자 모멘트(Q)를 새롭게 추정한다. 최종값 Q = 0.33666 ± 0.00258 b 로, 기존 권고값보다 약 10 % 크게 제시되며 최근 원자 계산 결과와도 일치한다.

상세 분석

이 논문은 핵 전기 사중극자 모멘트(Q)를 분자 수준에서 결정하는 두 단계 접근법을 체계적으로 구현한다. 첫 번째 단계는 SrO와 SrS의 핵 사중극자 결합 상수(NQCC)를 고정밀 마이크로파·라만 분광법 등 실험적 방법으로 확보한 뒤, 두 번째 단계에서 전자 구조 계산을 통해 해당 핵 주변의 전기장 구배(EFG)를 정확히 예측한다. EFG 계산에 있어 저자들은 정확한 두-성분 평균장(X2C‑mmf) 해밀토니안을 채택하고, 스칼라 상대론 효과와 스핀‑오빗 상호작용, 그리고 Gaunt 항까지 포함한 전자 상호작용을 전부 고려하였다.

전자 상관 효과는 핵심적인 오류 원인으로 인식되어, CCSD(T), CCSD‑T, 그리고 CCSD˜T(완전 반복 CCSDT에서 얻은 T 보정) 세 가지 고차 상관 방법을 모두 적용하였다. 특히 CCSD˜T는 작은 기저와 제한된 활성 공간에서 완전 반복 CCSDT 계산을 수행한 뒤, 그 차이를 큰 기저와 넓은 활성 공간에서의 CCSD 결과에 보정으로 더하는 복합 스킴이다. 이 과정에서 저자들은 새로운 병렬 텐서 라이브러리인 Cyclops Tensor Framework(CTF)를 TAPP 인터페이스를 통해 DIRAC 코드에 통합함으로써, 수천 개의 텐서 축소 연산을 수백 노드에 걸쳐 효율적으로 수행할 수 있었다.

기저 집합은 dyall.ae4z를 기본으로 하되, Sr 원자에 두 개의 d‑tight 함수(지수 1.45205990×10⁴, 6.52893664×10³)를 추가한 even‑tempered quadruple‑ζ 수준을 사용하였다. 활성 공간은 가상 궤도 에너지 컷오프 50 Eh로 제한했으며, 이는 EFG 수렴성을 실험적으로 검증한 결과이다. 또한, VIBCAL 유틸리티를 이용해 진동 보정도 포함함으로써 전자적 EFG와 실제 실험값 사이의 차이를 최소화하였다.

계산된 EFG와 실험 NQCC를 Eq.(1) NQCC = 234.9647 × Q × q에 대입하면, 두 분자 각각에 대해 Q 값을 얻을 수 있다. 두 값을 가중 평균하면 최종 Q = 0.33666 b이며, 통계적 불확도와 시스템적 오차를 합산한 총 불확도는 ±0.00258 b이다. 이 결과는 기존 권고값(≈0.305 b)보다 약 10 % 높으며, Lu et al. (2019)의 원자 계산값(0.328 b)과도 차이가 있지만, 최근 Tang et al. (2025)의 하이브리드 CI+CC 결과(0.336 b)와는 거의 일치한다.

결론적으로, 고차 전자 상관, 정확한 상대론 효과, 그리고 대규모 병렬 텐서 연산을 결합한 분자 수준의 EFG 계산이 핵 전기 사중극자 모멘트의 정확한 추정에 매우 유효함을 입증한다. 또한, 기존에 활용되지 않았던 SrO·SrS의 NQCC 데이터를 성공적으로 활용함으로써, 분자 기반 방법이 원자 기반 방법을 보완하거나 대체할 수 있는 가능성을 제시한다.