명시적 상관법

양자화학 계산화학 마지막 수정: 2026-08-16 04:23:07

1. 개요[편집]

명시적 상관법 Explicitly Correlated Methods (R12 / F12)
분야	양자화학 × 계산화학
핵심	파동함수에 r₁₂를 직접 넣는다
상관인자	슬레이터형 exp(−γ r₁₂)
기저 수렴	X^-3 → X^-7
대가	보조기저(CABS) + RI 근사
대표	MP2-F12, CCSD(T)-F12a/b

기저를 키워서 뾰족함을 흉내 내지 말고, 뾰족함을 손으로 그려 넣어라.

명시적 상관법(explicitly correlated method)은 다전자 파동함수에 전자 간 거리 $r_{12}$ 를 명시적인 변수로 포함하는 항을 집어넣어, 궤도 곱으로는 표현하기 어려운 전자-전자 첨점 근처의 구조를 직접 기술하는 방법론이다. 통칭 R12(선형 $r_{12}$ 를 쓰던 초기형)와 F12(일반 상관인자 $f(r_{12})$ 를 쓰는 현대형)로 불린다.

동기는 단순하다. 전자 상관 에너지는 기저 크기에 대해 $X^{-3}$ 이라는 굼뜬 속도로 수렴해서, 결합 클러스터 방법을 cc-pV5Z로 돌려야 겨우 소수점을 붙잡는 처지가 된다. F12를 붙이면 VDZ~~VTZ급 기저에서 QZ~~5Z급 값이 나온다. 비용이 $N^7$ 인 세계에서 기저를 두 단계 건너뛴다는 것은 계산 시간을 수백 분의 일로 줄인다는 뜻이라, 오늘날 고정밀 열화학은 F12를 기본값으로 깐다.¹

2. 카토 첨점 조건[편집]

출발점은 카토(Kato, 1957)가 증명한 첨점 조건(cusp condition)이다. 쿨롱 퍼텐셜은 두 입자가 만나는 지점에서 발산하는데, 슈뢰딩거 방정식의 해가 유한한 에너지를 가지려면 그 발산을 파동함수의 운동에너지 항이 정확히 상쇄해야 한다. 그러려면 파동함수가 합류점에서 뾰족한 꺾임을 가져야 한다.

핵-전자 합류점에서는 이렇게 쓴다( $Z$ 는 핵전하).

\left.\frac{\partial \bar{\Psi}}{\partial r_{iA}}\right|_{r_{iA}=0} = -Z\,\bar{\Psi}(r_{iA}=0)

전자-전자 합류점에서는 부호가 반대이고, 결정적으로 스핀에 따라 계수가 다르다.

\left.\frac{\partial \bar{\Psi}}{\partial r_{12}}\right|_{r_{12}=0} = \tfrac{1}{2}\,\bar{\Psi}(0) \;\;(\text{반평행 스핀}), \qquad \tfrac{1}{4}\,\bar{\Psi}(0) \;\;(\text{평행 스핀})

윗줄은 스핀이 다른 전자쌍(싱글렛 쌍)의 조건이다. 아랫줄이 스핀이 같은 쌍인데, 계수가 절반인 이유는 페르미 구멍 때문에 애초에 파동함수가 $r_{12}$ 에 비례해 0으로 사라지기 때문이다. 즉 평행 스핀 쌍은 첨점이 훨씬 약하다. 이 1/2과 1/4이라는 두 숫자는 뒤에 나올 고정 진폭 앙자츠에서 그대로 재활용된다.²

3. 왜 행렬식 전개로는 안 되는가[편집]

문제는 우리가 쓰는 도구다. 슬레이터 행렬식은 한 전자 함수의 곱이고, 곱함수는 어디서나 매끄럽다. 매끄러운 함수를 아무리 유한 개 더해도 $|r_{12}|$ 같은 꺾임이 생기지 않는다. 흉내는 낼 수 있지만 각운동량이 높은 함수를 끝없이 쌓아야 한다.

정량적으로는 부분파 전개로 따진다. 각운동량 $L$ 까지 포함했을 때 상관에너지의 증분이 $(L+\tfrac{1}{2})^{-4}$ 로 줄고, 꼬리를 다 합치면 잔여 오차가 $X^{-3}$ 이 된다. 이것이 완전기저 극한 외삽 공식의 뿌리이자, 동시에 외삽으로는 원인을 없앨 수 없다는 선언이기도 하다. 우리는 잘못된 함수 형태를 쓰고 있고, 외삽은 그 대가를 사후에 추정할 뿐이다.

명시적 상관법의 발상은 정공법이다. 참조 파동함수에 다음과 같은 게미널 항을 더한다.

|\Psi\rangle = |\Psi_{\text{conv}}\rangle + \sum_{i<j}\sum_{kl} t_{ij}^{kl}\,\hat{Q}_{12}\, f(r_{12})\,|kl\rangle

여기서 $\hat{Q}_{12}$ 는 새 항이 이미 궤도 기저로 표현되는 공간과 겹치지 않도록 잘라내는 강직교 사영연산자다. 이 항 하나 덕분에 잔여 오차가 $X^{-3}$ 에서 $X^{-7}$ 급으로 급감한다.

4. R12에서 F12로[편집]

1985, 쿠첼니크(Kutzelnigg) — $f(r_{12}) = r_{12}$ , 즉 선형항을 쓴 MP2-R12를 제안한다. 첨점 조건을 정확히 만족시키는 가장 단순한 선택이지만, 선형 함수는 $r_{12}$ 가 커져도 무한정 자라서 중·장거리 영역을 망친다. 그래서 아주 큰 기저에서만 잘 듣는다는 한계가 있었다.
1991, 쿠첼니크·클로퍼(Klopper) — 3전자·4전자 적분을 항등연산자 분해(RI)로 우회하는 기법을 확립한다. R12가 “이론적으로 아름답지만 못 돌리는” 신세에서 벗어난 결정적 계기.
2004, 텐노(Ten-no) — $f(r_{12})$ 를 슬레이터형 게미널(Slater-type geminal, STG)로 바꾼다.

f(r_{12}) = -\frac{1}{\gamma}\,e^{-\gamma r_{12}}

이 함수는 $r_{12} \to 0$ 에서 기울기가 정확히 1이라 첨점 조건을 맞추면서, $r_{12}$ 가 커지면 상수로 포화해 장거리를 오염시키지 않는다. 쿨롱 구멍의 실제 모양과도 훨씬 닮았다. $\gamma$ 는 보통 $1.0\ a_0^{-1}$ 안팎으로 잡되 기저에 맞춰 조정하며(작은 기저일수록 작게), 결과가 $\gamma$ 에 그리 민감하지 않다는 점이 실무적으로 큰 미덕이다. 실제 구현에서는 STG를 가우시안 6~10개의 합으로 적합해 기존 적분 코드를 재활용한다.

5. RI와 CABS — 공짜가 아닌 부분[편집]

게미널 항을 넣으면 해밀토니안 행렬 원소에 3전자·4전자 적분이 등장한다. 이걸 그대로 계산하면 F12는 태어나자마자 죽는다. 해법이 항등연산자 분해(resolution of the identity, RI)다.

\hat{1} \approx \sum_{p}^{\text{RI 기저}} |p\rangle\langle p|

이 근사를 끼워 넣으면 다전자 적분이 2전자 적분의 곱으로 쪼개진다. 대신 RI가 성립하려면 삽입되는 기저가 충분히 완전해야 한다. 그래서 궤도 기저와 별개로 큰 보조기저가 필요해지는데, 여기서 발레예프(Valeev, 2004)의 CABS(complementary auxiliary basis set)가 등장한다. 보조기저 전체를 쓰는 대신 궤도 기저가 이미 덮고 있는 공간을 직교로 빼낸 나머지만 쓰는 방식으로, 수치 안정성과 정확도를 동시에 잡았다. 오늘날 F12 계산의 입력에는 궤도 기저(cc-pVXZ-F12)·MP2 피팅 기저·CABS·밀도 피팅 기저가 줄줄이 딸려 오며, 이 짝이 안 맞으면 조용히 틀린 답이 나온다. F12의 가장 흔한 사고가 여기서 난다.³

6. 앙자츠와 근사 계층[편집]

논문 제목의 MP2-F12/3*C(FIX) 같은 암호는 세 개의 독립적인 선택이 붙은 것이다.

선택지	뜻	실무 기본값
앙자츠 1 / 2	사영연산자 $\hat{Q}_{12}$ 를 어떻게 정의할지	2 (점유·가상 공간 모두 배제)
근사 A / A’ / B / C	해밀토니안 교환자 항을 어디까지 유지할지	B 또는 C
최적화 / 고정(FIX, SP)	게미널 진폭 $t_{ij}^{kl}$ 을 변분 최적화할지 첨점 조건으로 못 박을지	고정

특히 세 번째가 재미있다. 진폭을 자유롭게 최적화하면 작은 기저에서 기저중첩오차와 비슷한 종류의 과결합이 생기고 수치적으로 불안정해진다. 그래서 텐노가 제안한 고정 진폭 앙자츠(SP, fixed-amplitude)는 대각 진폭을 앞서 본 첨점 계수 1/2과 1/4로 그냥 못 박아 버린다. 최적화를 포기했는데 결과가 더 좋고 기저중첩오차에도 강하며 크기 확장성까지 지켜지는, 보기 드문 공짜 점심이라 지금은 사실상 표준이다.

결합 클러스터 쪽에서는 아들러·크니지아·베르너(2007)의 CCSD(T)-F12a / F12b가 주력이다. 완전한 F12 이론을 CC에 이식하면 방정식이 감당 못 하게 불어나므로, 게미널 기여를 CCSD 방정식에 근사적으로 얹는 실용 노선이다. 삼중 들뜸 (T)에는 F12 보정이 직접 들어가지 않아서, 관례적으로 MP2-F12 쪽 기저 개선 비율을 곱해 스케일링한다 — 정직하게 말하면 땜빵이고, 그래서 (T) 부분은 여전히 기저에 민감하다.

7. 실측 감각과 비용[편집]

조합	대략 대응하는 종래 계산	메모
MP2-F12 / VDZ-F12	MP2 / aug-cc-pVQZ~5Z	상관에너지 기준
CCSD(T)-F12b / VTZ-F12	CCSD(T) / V5Z 급	열화학 실무의 국룰
CCSD(T)-F12b / VQZ-F12	사실상 CBS 근처	벤치마크급

비용은 공짜가 아니다. F12 항 때문에 MP2 기준 2~4배, CCSD 기준으로도 상당한 추가 시간이 붙고 메모리도 늘어난다. 그러나 기저를 한 단계 올리는 비용이 $1.6^7 \approx 27$ 배임을 생각하면 압도적으로 남는 장사다. 남은 함정만 짚자면 이렇다.

F12 전용 궤도 기저(cc-pVXZ-F12)는 일반 cc-pVXZ와 다르다. 이름이 비슷하다고 섞어 쓰면 안 된다. F12용 기저에 다시 CBS 외삽을 거는 것도 권장되지 않는다 — 수렴 법칙이 이미 바뀌었기 때문이다.
코어 상관을 다룰 거면 코어까지 상관시키는 기저(cc-pCVXZ-F12)와 frozen-core 해제가 함께 가야 한다.
F12는 동적 상관의 기저 수렴만 고친다. 정적 상관, 즉 다중참조 문제는 조금도 해결해주지 않는다. CCSD(T)-F12가 무너지는 계에서는 CCSD(T)도 똑같이 무너진다.
밀도 피팅과 CABS 기저가 없는 원소(주로 무거운 원소·유효 코어 퍼텐셜 조합)에서는 그냥 못 돌린다. 주기율표 아래쪽으로 갈수록 지원이 얇아진다.

이 모든 조각이 최종적으로 조립되는 곳이 컴포지트 방법이다. 현대의 W $n$ -F12나 HEAT 계열 프로토콜은 참조 상관 에너지를 F12로 뽑고, 그 위에 고차 들뜸·코어-밸런스·상대론 보정을 층층이 얹는다.

8. 관련 문서[편집]

재미있는 역사적 사실 하나. 1929년 힐레라스(Hylleraas)가 헬륨 바닥상태를 계산할 때 이미 $r_{12}$ 를 파동함수에 직접 넣어 당시로선 경이적인 정확도를 냈다. 즉 명시적 상관은 궤도 전개보다 먼저 나온 아이디어였는데, 다전자 적분이 감당이 안 돼 반세기 동안 창고에 처박혀 있었다. RI가 창고 문을 열어준 셈. ↩
여기서 “평행 스핀은 1/2, 반평행은 1/4”로 뒤집어 외우는 사고가 흔하다. 물리로 되짚으면 헷갈릴 일이 없다 — 페르미 구멍이 이미 파여 있는 평행 스핀 쪽이 추가로 팔 일이 적으므로 계수가 작다. ↩
증상이 고약하다. 잘못된 보조기저를 써도 계산은 멀쩡히 끝나고 에너지도 그럴듯하게 나온다. 그저 몇 mHartree 틀릴 뿐인데, 하필 그 mHartree를 보겠다고 F12를 켠 것이다. 출력 파일에서 CABS 싱글 보정값이 비정상적으로 크면 일단 기저 조합부터 의심하자. ↩