배치상호작용

양자화학 계산화학 마지막 수정: 2026-07-19 04:38:27

1. 개요[편집]

배치상호작용 Configuration Interaction
약칭	CI
앙자츠	슬레이터 행렬식의 선형 결합
성질	변분적(상계 보장)
극한	Full CI = 주어진 기저에서의 엄밀해
치명적 약점	절단 시 크기 확장성 상실

원리는 세상에서 제일 단순하다. 가능한 전자 배치를 전부 늘어놓고 섞으면 된다. 문제는 “전부”가 우주의 원자 수보다 빨리 늘어난다는 것.

배치상호작용(configuration interaction, CI)은 다전자 파동함수를 여러 개의 슬레이터 행렬식(전자 배치, configuration)의 선형 결합으로 쓰고, 그 계수를 변분 원리에 따라 최적화해 전자 상관 에너지를 얻는 post-하트리-폭 방법이다. 보통 HF 기준 행렬식 $|\Phi_0\rangle$ 에서 전자를 점유 궤도에서 비점유 궤도로 들뜨게 해 배치를 생성한다.

|\Psi_{\text{CI}}\rangle = c_0|\Phi_0\rangle + \sum_{ia} c_i^a |\Phi_i^a\rangle + \sum_{i<j,\,a<b} c_{ij}^{ab}|\Phi_{ij}^{ab}\rangle + \cdots

하트리-폭 방법은 전자가 서로의 평균장만 느낀다고 가정한다. 실제 전자는 서로를 순간순간 피해 다니는데, 행렬식 하나로는 이 회피 운동을 담을 수 없다. 여러 배치를 섞으면 “이 전자가 여기 있을 때 저 전자는 저기 있을 확률이 높다”는 상관을 표현할 수 있게 된다. CI는 그 상관을 담는 가장 직관적이고 가장 오래된 방식이며, 모든 후속 방법론의 언어를 제공한 조상 격 이론이다.¹

2. 변분성과 Full CI[편집]

CI 방정식은 배치 공간에서의 해밀토니안 행렬 고유값 문제로 귀결된다.

\mathbf{H}\mathbf{c} = E\,\mathbf{c}, \qquad H_{IJ} = \langle \Phi_I | \hat{H} | \Phi_J \rangle

계수를 변분적으로 최적화하니 CI 에너지는 항상 참값의 상계(upper bound)다. 계산 결과가 실제보다 낮게 나올 일이 없다는 이 성질은 벤치마크 관점에서 대단히 편안하다. 결합 클러스터 방법이 변분적이지 않아 이론적으로는 참값 아래로 내려갈 수 있는 것과 대비된다.

가능한 모든 들뜸을 다 포함하면 Full CI(FCI)가 되고, 이는 주어진 기저함수 안에서 비상대론적 전자 슈뢰딩거 방정식의 엄밀해다. 이론의 최종 답이 정의되어 있다는 점에서 CI는 특별하다. 문제는 규모다. $2K$ 개의 스핀궤도에 $N$ 개 전자를 넣는 경우의 수는 이항계수로 폭증한다.

N_{\text{det}} = \binom{2K}{N}

원자가 전자 몇십 개, 궤도 몇십 개 수준만 되어도 행렬식 수가 $10^{10}$ 을 가볍게 넘긴다. FCI가 실제로 계산 가능한 범위는 전자 20개 안팎, 궤도 20개 안팎의 소규모 계에 한정된다. 그래서 FCI의 실질적 역할은 다른 방법을 채점하는 기준점이다. 새 근사법이 나오면 작은 분자에서 FCI와 비교해 성적표를 받는다.

3. 절단과 브릴루앙 정리[편집]

현실에서는 들뜸 차수를 자른다. 단일(S)까지 자르면 CIS, 이중(D)까지면 CISD, 삼중까지면 CISDT 하는 식이다.

여기서 브릴루앙 정리(Brillouin’s theorem)가 중요한 사실을 알려준다. HF로 최적화된 궤도에서는 기준 행렬식과 단일 들뜸 사이의 해밀토니안 행렬 원소가 0이다.

\langle \Phi_0 | \hat{H} | \Phi_i^a \rangle = 0

즉 단일 들뜸만 넣으면 바닥상태 에너지가 전혀 안 내려간다. CIS의 바닥상태 에너지는 정확히 HF와 같다. 그래서 CIS는 바닥상태 상관법이 아니라 들뜬상태 계산법으로 쓰인다. 값싼 대신 정확도는 대체로 실망스럽지만, 큰 분자의 UV-Vis 스펙트럼을 개괄적으로 훑을 때는 여전히 등장한다.²

바닥상태 상관의 주력은 CISD다. 이중 들뜸이 상관 에너지의 대부분을 담당하고, 단일 들뜸은 궤도 완화(orbital relaxation)를 통해 간접적으로 기여한다. 비용은 대략 $N^6$ 로 CCSD와 같은 급이다.

4. 크기 확장성이라는 사형선고[편집]

그런데 CISD는 오늘날 거의 쓰이지 않는다. 같은 $N^6$ 비용의 CCSD가 모든 면에서 낫기 때문인데, 그 이유가 크기 확장성(size-extensivity)이다.

무한히 멀리 떨어진 동일한 분자 $n$ 개를 생각하자. 물리적으로 총 에너지는 정확히 단분자 에너지의 $n$ 배여야 한다. 그런데 절단된 CI는 이 조건을 어긴다. 이유는 선형 앙자츠에 있다. 분자 A에서 이중 들뜸, 분자 B에서 이중 들뜸이 동시에 일어난 상태는 전체 계에서 보면 사중 들뜸이다. CISD는 사중 들뜸을 포함하지 않으므로 이 상태를 통째로 빠뜨린다. 분자 수가 늘수록 빠뜨리는 비중이 커지고, 상관 에너지는 분자당 기준으로 점점 작아진다.

결합 클러스터 방법의 지수 앙자츠 $e^{\hat{T}}$ 는 $\frac{1}{2}\hat{T}_2^2$ 항이 자동으로 “독립적인 두 이중 들뜸”을 만들어내므로 이 문제가 없다. 선형 vs 지수 — 이 한 끗 차이가 CI를 실무 무대에서 밀어냈다.

응급처방으로 나온 것이 데이비드슨 보정(Davidson correction, CISD+Q)이다. 빠뜨린 사중 들뜸의 기여를 근사적으로 되돌려준다.

\Delta E_{Q} = \left(1 - c_0^2\right)\left(E_{\text{CISD}} - E_{\text{HF}}\right)

여기서 $c_0$ 는 기준 행렬식의 계수. 값싸고 꽤 효과적이지만 어디까지나 경험적 땜빵이며, 이 보정을 붙이는 순간 변분성이라는 CI의 유일한 자랑거리도 사라진다.

5. 거대 행렬을 다루는 기술[편집]

CI 행렬은 크기가 $10^9$ 이상으로 가는 대신 극도로 희소하다. 슬레이터-콘돈 규칙에 따라, 두 배치가 세 개 이상의 궤도에서 다르면 행렬 원소가 그냥 0이기 때문이다. 그래서 희소행렬 구조를 전제로 한 특수 기법이 발달했다.

데이비드슨 알고리즘: 최저 몇 개 고유값만 필요하다는 점을 이용한 부분공간 반복법. 크리로프 부분공간법의 사촌이며, 대각 우세한 행렬에 특화된 전처리를 쓴다. 이름 그대로 데이비드슨 보정과 같은 사람이다.
직접 CI(direct CI, Roos 1972): 행렬을 아예 저장하지 않고, 반복법이 필요로 하는 $\mathbf{H}\mathbf{c}$ 곱만 분자적분으로부터 그때그때 계산한다. 메모리 벽을 우회하는 고전적 해법으로, 고유값 문제를 푸는 현대적 관행의 원형이다.

6. MCSCF, MRCI, 그리고 부활[편집]

CI가 완전히 은퇴한 것은 아니다. 오히려 결합 클러스터 방법이 손대지 못하는 영역에서 대체 불가능하다. 결합이 끊어지는 전이상태, 이중결합 비틀림, 전이금속 착물, 들뜬상태처럼 여러 배치가 대등한 비중으로 섞이는 강상관(multireference) 계에서는, 애초에 “좋은 단일 기준 행렬식”이 존재하지 않는다.

여기서 등장하는 것이 MCSCF(multiconfigurational SCF)와 그 대표 격인 CASSCF(complete active space SCF)다. 화학적으로 중요한 궤도 $m$ 개와 전자 $n$ 개를 골라 활성 공간(active space, 흔히 $(n,m)$ 으로 표기)을 정의하고, 그 안에서는 Full CI를 하면서 동시에 궤도까지 최적화한다. 계수와 궤도를 같이 흔들기 때문에 수렴이 까다롭고, 활성 공간 선택은 여전히 사람의 화학적 직관에 의존하는 예술의 영역이다.³

CASSCF는 정적 상관(static correlation)만 담으므로, 동적 상관을 얹으려면 위에 한 층을 더 쌓는다. 섭동론으로 얹으면 CASPT2·NEVPT2, CI로 얹으면 MRCI(multireference CI)다. MRCI는 정확도가 매우 높지만 비용이 살벌하고 크기 확장성 문제도 여전해, 소분자 분광학의 정밀 계산에 주로 쓰인다. 이 계열 전체가 다중기준 방법이라 불린다.

최근에는 CI 자체를 되살리려는 시도도 활발하다.

선택 CI(selected CI): CIPSI, Heat-bath CI 등. 중요한 행렬식만 골라 담고 나머지는 섭동적으로 보정해, 사실상 FCI 품질을 훨씬 싸게 얻는다.
DMRG: 활성 공간을 수십 궤도까지 확장할 수 있게 해, CASSCF의 규모 한계를 크게 밀어냈다.
FCIQMC: 행렬식 공간을 확률적으로 표본하는 양자 몬테카를로 계열. 메모리 한계를 통계로 바꿔치기한다.

결국 CI는 “실무 주력”에서는 물러났지만, 정확한 답의 정의를 제공하고 강상관 영역을 지키는 자리에 여전히 앉아 있다. 계산화학에서 무엇이 옳은지 판정할 때는 결국 이쪽으로 돌아오게 된다.

7. 관련 문서[편집]

“배치상호작용”이라는 번역어가 어색하다는 지적이 많지만, 원어 configuration interaction도 오해를 부르기는 마찬가지다. 배치들이 물리적으로 “상호작용”하는 게 아니라, 해밀토니안 행렬의 비대각 원소를 통해 섞이는 것뿐이다. 용어가 만들어진 1930년대의 어법이 그대로 화석처럼 남았다. ↩
CIS의 정확도는 여기 에너지를 흔히 0.5~1 eV 이상 과대평가하는 수준이라, 정량적 예측용으로 쓰기는 어렵다. 그래도 TD-DFT가 실패하는 전하이동 상태에서는 CIS가 오히려 정성적으로 덜 틀리는 경우가 있어, 아직 완전히 퇴출되지 않았다. ↩
활성 공간을 고르는 일은 논문 심사에서 가장 자주 공격받는 지점이다. “왜 (12,12)를 골랐나”에 대한 정답은 대개 “(14,14)는 안 돌아가서”인데, 그렇게 쓸 수는 없으니 다들 화학적 정당화를 열심히 지어낸다. ↩