B스플라인으로 페시보드 분자와 할로 상태의 고정밀 스펙트럼 계산

초록

**
본 논문은 B‑스플라인 기반의 유한 기저법을 이용해 이원자 분자의 로-진동 고유스펙트럼, 특히 해리 한계 근처의 얇은(halo) 및 페시보드(Feshbach) 분자 상태를 정확하고 효율적으로 계산한다. 기존의 이산 변수 표현(DVR) 방식과 비교했을 때, B‑스플라인은 ‘고스트’ 상태 발생을 억제하고 에너지 정확도를 10⁻⁷ cm⁻¹ 수준까지 향상시킨다. 또한 1/R⁶ 장거리 퍼텐셜에 대해 마지막 결합 상태의 에너지와 산란 길이 사이의 양자 결함 이론(QDT) 관계를 수치적으로 검증한다.

**

상세 분석

**
본 연구는 두 가지 핵심 기술적 진보를 제시한다. 첫째, B‑스플라인을 Galerkin 방식으로 적용해 라디얼 슈뢰딩거 방정식의 일반화 고유값 문제를 구성하고, 경계조건을 만족하도록 첫·마지막 스플라인을 제외한 n‑2개의 기저함수를 사용한다. 스플라인 차수 k=15와 적절히 설계된 비균등 격자(매핑 격자)를 결합함으로써, 짧은 거리의 급격한 퍼텐셜 변화와 장거리 1/R³·1/R⁶ 꼬리 모두를 고해상도로 포착한다. 매핑 격자는 로컬 드브로이 길이에 따라 격자 간격을 자동 조정하므로, 해리 한계 근처에서 파동함수의 지수적 감쇠와 다수의 노드가 존재하는 영역을 효율적으로 샘플링한다. β 파라미터를 0.4~0.7 사이로 조정해 격자 밀도를 최적화하면서도 ‘고스트’ 상태를 최소화한다는 점이 특히 눈에 띈다.

둘째, B‑스플라인 결과를 Morse 퍼텐셜(분석 해석이 가능한 기준)과 실제 Cs₂, Li₂ 등 1/R³·1/R⁶ 장거리 퍼텐셜에 적용해 비교 검증하였다. Morse 경우, B‑스플라인은 상대 오차 10⁻¹¹ cm⁻¹ 수준을 달성했으며, DVR은 10⁻⁷ cm⁻¹ 수준에 머물렀다. 특히 해리 한계 근처에서 DVR은 ‘고스트’ 고유값(v=168 등)을 생성해 결합 상태 수를 잘못 예측했지만, B‑스플라인은 격자 크기와 무관하게 이러한 인공 상태가 전혀 나타나지 않았다. 실제 1/R⁶ 퍼텐셜에 대해서는 마지막 결합 상태 에너지 E₀와 산란 길이 aₛ 사이의 관계를 QDT식
(E₀ = -\frac{\hbar^2}{2\mu (aₛ - \bar a)^2})
(여기서 (\bar a)는 평균 길이)와 비교했을 때, 수치 결과가 QDT 예측을 넓은 파라미터 구간에서 정확히 재현함을 확인했다. 이는 기존 유효 범위 이론보다 훨씬 넓은 적용 범위를 제공한다는 의미다.

또한, B‑스플라인은 전역적인 비선형 최소화가 필요 없는 Galerkin 접근법을 사용함으로써, 대규모 행렬 대각화만으로 고유값을 얻을 수 있다. 메모리 요구량은 기저 함수 수에 비례하지만, 스플라인이 로컬 지원을 갖기 때문에 행렬이 희소(sparse)하고 효율적인 Lanczos·Arnoldi 알고리즘 적용이 가능하다. 결과적으로, 수천 개의 결합 상태를 포함하는 복잡한 퍼텐셜도 수십 분 내에 정확히 계산할 수 있다.

이러한 기술적 장점은 초저온 원자·분자 물리학, 특히 Feshbach 공명에 의한 분자 결합, 광학 라인스펙트럼, 그리고 양자 결함 이론 검증 등에 직접적인 활용 가능성을 열어준다.

**