전자 상관

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

1. 개요[편집]

전자 상관 Electron Correlation
분야	양자화학 × 계산화학
정의	정확한 비상대론적 에너지 − HF 극한 에너지
부호	항상 음수
크기	총에너지의 약 1%, 화학의 거의 전부
두 얼굴	동적 상관 / 정적(비동적) 상관
기저 수렴	X^-3 (느림)

총에너지의 1%. 그런데 화학은 전부 그 1% 안에서 일어난다.

전자 상관(electron correlation)은 전자들이 서로의 평균장이 아니라 순간의 위치에 반응해 서로를 피해 다니는 효과 전체를 가리키며, 계산화학에서는 그 효과가 에너지에 기여하는 몫인 상관에너지(correlation energy)로 정량화된다. 뢰브딘(Löwdin)의 표준 정의는 이렇다.

E_{\text{corr}} = E_{\text{exact}}^{\text{nonrel}} - E_{\text{HF}}^{\text{limit}}

두 항 모두 완전기저 극한에서, 같은 핵 배치, 같은 비상대론적 해밀토니안 기준이라는 단서가 붙는다. 하트리-폭 방법은 변분적이므로 $E_{\text{HF}} \ge E_{\text{exact}}$ 이고, 따라서 $E_{\text{corr}}$ 은 정의상 항상 음수다. “상관을 많이 회수했다”는 말은 곧 “더 음수로 내려갔다”는 뜻이다.¹

크기는 초라하다. 물 분자 하나의 총 전자에너지가 −76 하트리 남짓인데 상관에너지는 −0.3 하트리 정도, 약 0.4%다. 그런데 화학 결합 하나가 대략 0.1~0.2 하트리이고, 우리가 다투는 반응 장벽은 1 kcal/mol = 0.0016 하트리 스케일이다. 총에너지에서 상관은 잔돈이지만, 에너지 차이에서 상관은 원금이다.

2. 정의를 둘러싼 잔소리[편집]

정의가 깔끔해 보여도 실전에서는 함정이 세 개 있다.

어떤 HF인가. 개껍질(open-shell)에서 RHF/ROHF/UHF 중 무엇을 기준으로 잡느냐에 따라 $E_{\text{corr}}$ 의 값이 달라진다. UHF는 스핀 오염을 대가로 상관 일부를 이미 먹고 들어가므로 남은 상관에너지가 작게 나온다. 숫자를 비교하려면 기준부터 맞춰야 한다.
어떤 기저인가. 유한 기저에서 계산한 $E_{\text{FCI}} - E_{\text{HF}}$ 는 엄밀히는 기저 안에서의 상관에너지다. 진짜 $E_{\text{corr}}$ 은 완전기저 극한에서만 정의된다.
상대론은 별도 항목. 무거운 원소에서 스칼라 상대론 효과는 상관보다 클 수도 있는데, 정의상 상관에너지에 포함되지 않는다. 두 항을 따로 계산해 더하는 것이 컴포지트 방법의 기본 문법이다.

3. 동적 상관과 정적 상관[편집]

상관은 물리적으로 성격이 전혀 다른 두 갈래로 나뉜다. 이 구분이 어떤 방법을 쓸지를 결정하므로, 계산화학 실무의 첫 갈림길이 여기다.

동적 상관(dynamic correlation)은 전자 두 개가 가까워질 때 서로 밀쳐내며 순간적으로 회피하는 단거리 효과다. 기여하는 배치 하나하나는 작지만 개수가 어마어마해서, 무수히 많은 작은 계수의 들뜸이 합쳐져 만들어진다. 개별 계수가 작으므로 섭동론이나 결합 클러스터 방법처럼 단일 기준에서 출발하는 방법이 잘 듣는다.

정적 상관(static correlation, 비동적 상관 nondynamic correlation)은 반대다. 애초에 비슷한 에너지의 전자배치가 여러 개 대등하게 섞여 있어서 하나의 슬레이터 행렬식으로는 출발점조차 못 만드는 상황을 말한다. 소수의 배치가 각각 큰 계수를 갖는다.

교과서 예제는 H $_2$ 해리다. 결합이 늘어나면 $\sigma_g^2$ 와 $\sigma_u^2$ 배치가 에너지적으로 붙어버려서, RHF는 원자 두 개가 아니라 이온쌍 성분을 절반 섞은 엉뚱한 극한으로 발산한다. 이 지점에서 CCSD(T)의 (T) 보정이 폭주하며 에너지 곡선을 아래로 찢어놓는 광경은 양자화학 수업의 단골 장면이다.²

갈림길 신호	해석	처방
$T_1$ 진단값이 0.02 초과(닫힌껍질 기준)	단일 기준 가정이 흔들림	다중기준 검토
자연 궤도 점유수가 0이나 2에서 크게 벗어남	정적 상관 존재	CASSCF 활성 공간 설정
(T) 보정이 CCSD 상관의 10%를 넘음	섭동적 삼중항 붕괴 조짐	CCSDT(Q)나 다중기준 방법
결합 해리·전이금속·라디칼·들뜬상태	구조적으로 다중참조	CASSCF+NEVPT2, MRCI, DMRG

물론 이 구분은 연속적이고 정의가 애매하다. “정적 상관 30%“처럼 딱 잘라 말할 방법이 없고, 그래서 이 바닥에서 가장 자주 벌어지는 리뷰어-저자 싸움의 소재이기도 하다.

4. 쿨롱 구멍과 페르미 구멍[편집]

상관을 에너지가 아니라 공간 분포로 보면 그림이 훨씬 선명해진다. 전자 하나를 어딘가에 고정했을 때 다른 전자를 발견할 확률이 주변에서 얼마나 파이는지를 보면 된다.

페르미 구멍(Fermi hole) — 반대칭성에서 온다. 같은 스핀 전자는 파울리 배타 원리 때문에 같은 자리에 있을 수 없으므로 확률이 정확히 0으로 파인다. HF 행렬식은 이 구멍을 이미 정확히 갖고 있다. 교환에너지가 곧 이 구멍의 값이다.
쿨롱 구멍(Coulomb hole) — 쿨롱 반발에서 온다. 스핀에 상관없이 두 전자는 서로를 피하므로 $r_{12} \to 0$ 근처의 확률이 더 파인다. HF 행렬식에는 이 구멍이 없다. 다른 스핀 전자 사이에는 회피 장치가 아예 없기 때문이다.

즉 상관에너지의 정체는 쿨롱 구멍을 파느라 얻는 에너지 이득이다. 이 구멍을 실제 파동함수에서 뽑아내 그려 보인 것이 쿨슨-닐슨(Coulson–Neilson, 1961)의 헬륨 계산으로, 정확한 파동함수와 HF 파동함수의 전자쌍 분포 차이를 $r_{12}$ 의 함수로 그리자 원점 근처가 움푹 파이고 바깥이 부풀어 오른 모양이 나왔다. “상관 = 구멍”이라는 직관이 여기서 굳어졌다.

형식적으로는 밀도행렬 언어가 가장 깔끔하다. 전자 해밀토니안은 한 몸 항과 두 몸 항뿐이므로 총에너지는 2차 축소 밀도행렬(2-RDM) $\Gamma_2$ 만 알면 정확히 결정된다.

E = \operatorname{Tr}(\hat{h}\,\gamma_1) + \tfrac{1}{2}\operatorname{Tr}(\hat{v}_{12}\,\Gamma_2)

HF는 $\Gamma_2$ 를 1차 밀도행렬 $\gamma_1$ 의 반대칭 곱으로 인수분해해 버린다(교환항까지만). 상관이란 정확히 그 인수분해에서 벗어난 나머지다. 그렇다면 $\Gamma_2$ 만 변분적으로 최적화하면 되지 않나 싶은데, 실제로 그렇게 하면 **어떤 파동함수로도 만들어질 수 없는 $\Gamma_2$ **가 튀어나온다. 이것이 악명 높은 $N$ -표현가능성(N-representability) 문제이고, 2-RDM 직접 최적화가 아직 주류가 못 된 이유다.³

5. 누가 얼마나 긁어오는가[편집]

주어진 기저 안에서 Full CI가 100%라고 두고, 얌전한 닫힌껍질 분자에서의 감각값은 대략 이렇다. 계에 따라 크게 흔들리니 숫자를 외우지는 말 것.

방법	비용	상관 회수(감각값)	성격
MP2	$N^5$	80~90%	싸고 빠름, 분산은 과대평가
CCSD	$N^6$	95% 내외	단일·이중 완전 포함
CCSD(T)	$N^7$	99% 내외	금본위, 화학적 정확도
CISD	$N^6$	CCSD보다 낮음	크기 확장성 위반
Full CI	지수적	100%(정의)	채점 기준

여기서 배치상호작용이 왜 밀려났는지가 다시 드러난다. 같은 $N^6$ 을 쓰고도 회수율이 낮고, 계가 커질수록 분자당 회수율이 떨어지는 크기 확장성 위반까지 안고 있다. 반면 밀도범함수이론은 상관을 아예 다른 방식으로 — 파동함수 대신 교환-상관 범함수로 — 근사한다. 그래서 DFT에는 “상관 회수율 몇 %“라는 표현이 성립하지 않는다. 회수하는 게 아니라 모형으로 대체하는 것이기 때문이다.

6. 왜 기저가 그렇게 안 수렴하는가[편집]

상관 계산의 진짜 고통은 방법이 아니라 기저에서 온다. HF 에너지는 기저 크기에 대해 지수적으로 수렴해 cc-pVQZ쯤이면 사실상 끝나는데, 상관에너지는 $X^{-3}$ 이라는 굼뜬 멱함수로 기어간다. cc-pV6Z를 써도 안 끝난다.

이유는 쿨롱 구멍의 모양에 있다. 정확한 파동함수는 두 전자가 겹치는 $r_{12} \to 0$ 에서 뾰족한 꺾임(cusp)을 갖는데, 카토(Kato)의 첨점 조건이 그 기울기를 못 박아 놓았다. 그런데 우리가 쓰는 기저는 전부 한 전자 함수의 곱이라, 매끄러운 곱함수들로 $|r_{12}|$ 꼴의 꺾임을 흉내 내려면 높은 각운동량 함수를 끝없이 쌓아야 한다. 이 느림은 기저 설계의 게으름이 아니라 행렬식 전개라는 형식 자체의 한계다.

대응책은 두 갈래다. 하나는 굼뜬 수렴 법칙을 알고 있으니 외삽하자는 노선(완전기저 극한), 다른 하나는 아예 파동함수에 $r_{12}$ 를 손으로 집어넣어 첨점을 직접 그려주는 노선(명시적 상관법)이다. 오늘날 고정밀 열화학은 둘을 함께 쓰고, 거기에 코어-밸런스·상대론 보정까지 층층이 쌓아 컴포지트 방법이라는 레시피로 굳혔다.

7. 현업 감각[편집]

“상관 없이 결합 에너지를 논하지 마라”는 말은 과장이 아니다. HF는 결합을 체계적으로 과소평가하고, 그 오차가 수십 kcal/mol에 이른다.
반대로 상관을 넣는다고 다 좋아지는 것도 아니다. MP2는 $\pi$ -스택 같은 분산 지배 상호작용을 20~50% 과대평가하는 것으로 악명이 높아, SCS-MP2처럼 스핀 성분별로 계수를 다르게 주는 처방이 따로 발달했다.
국소성이 열쇠다. 동적 상관은 본질적으로 단거리라, 멀리 떨어진 전자쌍의 기여를 잘라내는 국소 상관법(DLPNO-CCSD(T) 등)이 정확도를 거의 잃지 않고 $N^7$ 을 준선형까지 끌어내린다.
결국 실무의 질문은 “상관을 넣을까”가 아니라 **“이 계가 단일참조인가”**다. 여기서 오판하면 아무리 비싼 계산을 돌려도 답이 틀린다.

8. 관련 문서[편집]

그래서 논문에서 상관에너지를 양수로 적고 “회수량”이라 부르는 경우가 있는데, 부호 규약을 명시하지 않으면 표 하나로 독자를 30분간 괴롭힐 수 있다. 정식 정의는 음수다. ↩
RHF의 H $_2$ 해리 곡선이 원자 두 개의 합보다 한참 위로 올라가는 것, 그리고 UHF가 어느 결합 길이(Coulson–Fischer 점)에서 RHF 해에서 갈라져 나오며 그 오류를 스핀 오염과 맞바꾸는 것 — 이 두 그림은 정적 상관을 처음 배울 때 반드시 한 번은 직접 그려봐야 한다. ↩
$N$ -표현가능성 조건을 전부 적어내는 것은 계산 복잡도상 QMA-난해로 알려져 있다. 즉 “2-RDM만 최적화하면 되잖아”라는 아름다운 발상은, 어려움을 없앤 게 아니라 다른 주머니로 옮긴 것이었다. 실무에서는 2-양(positivity) 조건 몇 개만 걸고 반정부호 계획법으로 푸는 근사가 쓰인다. ↩