완전기저 극한

양자화학 수치해석 마지막 수정: 2026-07-25 04:45:33

1. 개요[편집]

완전기저 극한 Complete Basis Set limit (CBS)
분야	양자화학 × 수치해석
수렴 파라미터	기저 카디널 수 X (cc-pVXZ)
HF 수렴	지수적
상관에너지 수렴	X^-3
우회로	F12 / 명시적 상관법

기저를 무한대로 키운 답을, 기저를 무한대로 키우지 않고 알아내는 기술.

완전기저 극한(complete basis set limit, CBS limit)은 기저함수 집합의 크기를 무한대로 키웠을 때 도달하는, 기저에 더는 의존하지 않는 이론 고유의 값이다. 실제 계산은 언제나 유한 기저에서 돌아가므로 CBS 값은 직접 계산되지 않고, 크기가 다른 몇 개의 유한 기저 결과를 알려진 수렴 법칙에 맞춰 외삽해서 추정한다.

발상 자체는 CFD·구조해석 쪽의 메시 독립성 연구나 리처드슨 외삽법과 정확히 같다. 격자를 반씩 줄이며 $h \to 0$ 답을 추정하듯, 기저를 D→T→Q로 키우며 $X \to \infty$ 답을 추정한다. 다른 점이 하나 있다면, 양자화학의 수렴이 훨씬 느리고 훨씬 비싸다는 것. 격자 절반은 비용 8배지만, cc-pVXZ에서 X를 하나 올리면 결합 클러스터 방법 비용이 10배 넘게 뛴다.¹

2. 변분 원리가 보장하는 것[편집]

기저를 키운다는 것은 시행 파동함수가 사는 부분공간을 넓힌다는 뜻이다. 변분법의 원리에 따라 하트리-폭 방법이나 배치상호작용처럼 변분적인 방법의 총에너지는 항상 참값보다 위에 있고, 부분공간이 포함관계로 커지면 위에서부터 단조적으로 내려간다.

E_{N_1} \ge E_{N_2} \ge \cdots \ge E_{\text{CBS}} \ge E_{\text{exact}}

이 단조성이 CBS 추정의 뼈대다. 수렴이 한 방향이라는 걸 알면 외삽이 성립하고, 곡선이 위아래로 튀면 기저 선택이나 SCF 수렴부터 의심해야 한다는 진단 기준도 생긴다.²

상자 안 1D 조화진동자를 사인 함수 N개짜리 기저로 실제 변분 계산한 결과. N을 늘리면 바닥에너지 추정값이 위에서부터 단조 감소하며 정확값 0.5로 접근하고, 오차를 로그 축에 그리면 수렴 속도가 직선으로 드러난다. 어디까지나 1D 장난감 모형이며 실제 분자 기저집합(cc-pVXZ)이 아니다 — 특히 여기엔 전자-전자 첨점이 없어서, 실제 상관에너지의 X^-3 수렴보다 훨씬 빠르게 수렴한다.

3. 상관일치 기저의 계통적 수렴[편집]

아무 기저나 키운다고 외삽이 되는 건 아니다. 6-31G에서 6-311G로 가는 식의 포플(Pople) 계열은 크기가 커지긴 해도 어떤 법칙으로 수렴하는지가 정의되어 있지 않다. 던닝의 상관일치 기저(correlation-consistent, cc-pVXZ, X = D·T·Q·5·6)는 애초에 이 문제를 겨냥해 설계됐다. 상관에너지에 비슷한 크기로 기여하는 함수들을 한 껍질(shell)로 묶어 한꺼번에 추가하기 때문에, X를 올릴 때 에너지가 정해진 꼴로 줄어든다. 음이온·약한 상호작용·들뜬 상태를 다룰 때는 확산 함수를 덧댄 aug-cc-pVXZ를 쓴다.

핵심은 HF 부분과 상관 부분의 수렴 속도가 질적으로 다르다는 사실이다.

성분	수렴 거동	실무적 의미
HF(또는 SCF) 에너지	지수적	X = Q~5면 사실상 수렴
전자 상관 에너지	X^-3 멱함수	X = 6에서도 안 끝남

그래서 총에너지를 통째로 외삽하는 것은 두 개의 다른 함수를 하나의 식에 욱여넣는 짓이다. 실무의 국룰은 분리 외삽이다.

4. 두 점·세 점 외삽 공식[편집]

상관에너지는 헬가커(Helgaker) 두 점 공식이 표준이다. $E_X^{\text{corr}} = E_\infty^{\text{corr}} + A X^{-3}$ 를 연속한 두 카디널 수에 대해 풀면 $A$ 가 소거된다.

E_\infty^{\text{corr}} = \frac{X^3 E_X^{\text{corr}} - Y^3 E_Y^{\text{corr}}}{X^3 - Y^3}, \qquad Y = X-1

HF는 지수 꼴이므로 별도로 처리한다. 카튼-마틴(Karton–Martin) 두 점 공식이 흔히 쓰인다.

E_X^{\text{HF}} = E_\infty^{\text{HF}} + A\,(X+1)\,e^{-9\sqrt{X}}

세 점을 쓸 수 있으면 $E_X = E_\infty + A e^{-bX}$ 를 세 값으로 비선형 적합하는 펠러(Feller) 방식도 쓴다. 다만 점이 늘어난다고 항상 안전해지지는 않는다. DZ(X=2)는 수렴 점근 영역에 들어와 있지 않아서 외삽에 넣으면 오히려 결과를 망치는 일이 잦다 — (T,Q) 조합이 사실상 표준이고, DZ는 넣지 않는다는 것이 축적된 경험칙이다.³

무엇을 외삽할지도 중요하다. 원자화 에너지처럼 절대 총에너지 차이가 그대로 드러나는 양은 CBS가 필수지만, 반응 엔탈피·장벽 높이·상호작용 에너지처럼 양쪽에서 비슷한 오차가 상쇄되는 상대량은 훨씬 작은 기저에서도 쓸 만한 값이 나온다. 반대로 평형 구조나 진동 주파수는 대개 TZ급에서 이미 안정되므로, 구조·주파수는 싼 수준에서 잡고 그 위에서 단일점 에너지만 CBS로 올리는 계층적 프로토콜이 표준이 됐다. 계산 조건을 CCSD(T)/CBS//B3LYP/def2-TZVP처럼 이중 슬래시로 적는 표기가 바로 이 분업을 박제한 것이다.

5. 왜 이렇게 느린가 — 전자 첨점과 F12[편집]

$X^{-3}$ 이라는 굼뜬 수렴에는 명확한 물리적 이유가 있다. 참 파동함수는 두 전자가 겹치는 지점 $r_{12} \to 0$ 에서 뾰족한 꺾임, 즉 전자-전자 첨점(cusp)을 갖는다. 카토(Kato)의 첨점 조건은 이렇게 쓴다.

\left.\frac{\partial \bar{\Psi}}{\partial r_{12}}\right|_{r_{12}=0} = \frac{1}{2}\,\bar{\Psi}(r_{12}=0)

문제는 우리가 쓰는 기저가 전부 한 전자 함수의 곱이라는 점이다. 매끄러운 곱함수들로 $|r_{12}|$ 같은 꺾임을 흉내 내려면 각운동량이 높은 함수를 끝없이 넣어야 한다. 부분파 전개에서 각운동량 $L$ 까지 넣었을 때 남는 증분이 $(L+\tfrac{1}{2})^{-4}$ 로 줄고, 이를 꼬리까지 합하면 정확히 $X^{-3}$ 이 나온다. 즉 이 느림은 기저 설계의 게으름이 아니라 곱함수 전개라는 형식 자체의 한계다.

우회로가 F12 / 명시적 상관법(명시적 상관법)이다. 파동함수에 $r_{12}$ 를 명시적으로 담은 항(현대 F12는 슬레이터형 게미널 $-\gamma^{-1}e^{-\gamma r_{12}}$ )을 직접 집어넣어 첨점을 손으로 그려 넣는다. 그러면 남은 부분파 증분이 $(L+\tfrac{1}{2})^{-7}$ 수준으로 급감하고, 실무적으로 CCSD(T)-F12/VTZ-F12가 통상 CCSD(T)/V5Z급 값을 준다. 기저 두 단계를 공짜로 건너뛰는 셈이라, 오늘날 고정밀 열화학은 F12를 기본값으로 깐다.

참고로 밀도범함수이론은 이 저주에서 상당히 자유롭다. 콘-샴 행렬식은 첨점을 재현할 의무가 없어서, DFT 총에너지는 대개 TZ~QZ에서 기저 수렴이 사실상 끝난다. “DFT는 TZ면 되는데 CCSD(T)는 왜 QZ에 외삽까지 하냐”는 물음의 답이 여기 있다.

6. 평면파 쪽 대응물과 실무 프로토콜[편집]

평면파 기저에서의 대응물은 컷오프 에너지 수렴 시험이다. 노브가 $E_{\text{cut}}$ 하나뿐이고 역시 변분적·단조적으로 완전기저에 접근하므로, 컷오프를 400·500·600 eV로 올려가며 관심 물리량이 목표 허용오차 안으로 들어오는 지점을 잡는다. 여기서는 외삽 공식을 쓰기보다 그냥 충분히 큰 컷오프를 쓰는 쪽이 관례인데, 계산이 한 전자 이론(DFT)이라 수렴이 빠르고 컷오프를 올리는 비용이 $X$ 를 올리는 비용보다 훨씬 완만하기 때문이다. 대신 유사퍼텐셜 데이터셋이 권장 컷오프를 지정하므로 그 값이 사실상 출발선이 된다.

파동함수 쪽 실무는 컴포지트 방법(컴포지트 방법)으로 굳어졌다. 값싼 계산 여러 개를 조합해 비싼 계산 하나를 흉내 내는 레시피다.

G4 — 여러 이론/기저 조합의 가법적 보정 + 경험적 고준위 보정. G3/05 시험집합에서 평균절대편차 약 0.8 kcal/mol로 화학적 정확도(1 kcal/mol) 언저리.
W1/W2 — 경험 파라미터를 거의 없애고 HF·CCSD·(T) 성분을 각각 다른 기저 계열로 분리 외삽. W1이 0.5 kcal/mol 수준, 상위 W4는 kJ/mol 이하까지 내려간다.
초점점 분석(focal-point analysis) — 이론 수준과 기저 크기를 2차원 표로 깔고, 각 축 방향의 증분을 따로 외삽하는 방식.

공통 철학은 하나다. 기저 수렴과 이론 수렴을 분리해서 각각 외삽한다. 그리고 이 값에는 기저중첩오차가 남아 있을 수 있는데, BSSE는 정의상 CBS 극한에서 0으로 사라지므로 외삽이 잘 되면 그만큼 함께 정리된다 — 자세한 처방은 해당 문서 참고.

7. 관련 문서[편집]

CCSD(T)의 형식 비용은 기저 크기 $N$ 에 대해 $N^7$ 이다. cc-pVTZ→cc-pVQZ에서 함수 개수가 대략 1.5~1.7배 늘어나니, 1.6의 7제곱 ≈ 27배. “Q 한 번 돌려보자”는 말이 왜 사수의 표정을 굳게 만드는지 설명이 된다. ↩
엄밀히 말하면 cc-pVXZ 계열은 서로 완전한 포함관계(nested)가 아니라서 단조성이 수학적으로 보장되지는 않는다. 다만 실전에서 어긋나는 경우는 드물고, 어긋났다면 십중팔구 SCF가 다른 해로 수렴했거나 선형 종속이 발생한 것이다. MP2·CCSD(T)의 상관에너지도 비변분적이라 원리적 보장은 없지만 경험적으로는 얌전히 한 방향으로 수렴한다. ↩
“점을 많이 넣을수록 좋다”는 직관이 외삽에서는 통하지 않는다. 점근 법칙이 성립하지 않는 영역의 데이터는 정보가 아니라 잡음이다. 격자 수렴 연구에서 너무 성긴 격자를 GCI 계산에 넣으면 안 되는 것과 완전히 같은 이야기. ↩