대칭적응 섭동론

양자화학 계산화학 마지막 수정: 2026-08-18 04:49:03

1. 개요[편집]

대칭적응 섭동론 Symmetry-Adapted Perturbation Theory
약칭	SAPT
분야	양자화학 × 분자간 상호작용
구하는 것	상호작용 에너지를 직접 (차이가 아니라)
성분	정전기 · 교환 · 유도 · 분산 (+ 각각의 교환 짝)
BSSE	구조적으로 발생하지 않음
계층	SAPT0 → SAPT2+ → SAPT2+3, 그리고 SAPT(DFT)

“이 리간드가 왜 더 세게 붙나요?” — 총 에너지 하나만 들고는 답할 수 없는 질문이다. SAPT는 그 숫자를 물리 성분으로 쪼개서 돌려준다.

대칭적응 섭동론(Symmetry-Adapted Perturbation Theory, SAPT)은 두 분자 A와 B의 상호작용 에너지를 각각의 단량체 파동함수에서 출발해 섭동 전개로 직접 계산하는 방법이며, 그 과정에서 상호작용을 정전기·교환반발·유도·분산이라는 물리적으로 의미 있는 성분으로 자동 분해한다는 점이 결정적 특징이다. 예지오르스키(Jeziorski), 셰프치크(Szalewicz), 모시니스키(Moszynski) 등 폴란드 학파가 1970년대부터 다듬어온 체계다.

이름의 “대칭 적응”은 단량체 사이의 전자 교환에 대한 반대칭성을 섭동 전개 각 차수에서 강제로 회복시킨다는 뜻이다. 이 한 단어가 이 이론의 존재 이유 전부이므로 아래에서 따로 다룬다.

2. 초분자 방식의 문제[편집]

상호작용 에너지를 구하는 표준적인 방법은 뺄셈이다.

E_{\text{int}} = E_{AB} - E_A - E_B

이 초분자(supermolecular) 방식은 개념이 단순하고 어떤 전자구조 방법에도 얹을 수 있다는 장점이 있지만, 실무에서 세 가지 대가를 치른다.

기저중첩오차(BSSE). 이량체 계산에서는 A가 B의 기저함수를 빌려 쓸 수 있어 $E_{AB}$ 가 부당하게 낮아진다. 카운터포이즈 보정(보이스-베르나르디)으로 단량체도 이량체 기저에서 다시 계산해 상쇄해야 하고, 계산이 세 번으로 늘어난다.
큰 수의 작은 차이. 물 이량체를 예로 들면 $E_{AB}$ 는 −152 하트리 스케일인데 우리가 원하는 답은 −0.008 하트리다. 두 큰 수의 다섯 자리 아래 차이를 신뢰하려면 두 계산이 정확히 같은 수준으로 수렴해야 한다.
해석 불가. 답은 숫자 하나뿐이다. 정전기 때문인지 분산 때문인지, 교환반발이 얼마나 상쇄했는지는 알 길이 없다. 사후에 에너지 분해 분석(EDA)을 붙이는 방법들이 있지만 그건 이미 얻은 답을 쪼개는 것이고, 스킴마다 답이 달라진다.

SAPT는 뺄셈을 하지 않는다. 상호작용 항 $V$ (A의 핵·전자와 B의 핵·전자 사이 쿨롱 상호작용 전부)를 섭동으로 두고 상호작용 에너지 자체를 급수로 전개한다. 그러니 큰 수의 차이 문제도, 초분자적 의미의 BSSE도 원리적으로 생기지 않고, 성분 분해는 계산의 부산물이 아니라 정의 그 자체다.

3. 왜 반대칭화가 필요한가[편집]

가장 소박한 발상은 이렇다. 해밀토니안을 $H = H_A + H_B + V$ 로 쓰고, 무섭동 파동함수를 단순곱 $\Phi_A\Phi_B$ 로 놓고, 통상적인 레일리-슈뢰딩거 섭동론을 $V$ 에 대해 돌린다. 이것을 분극 전개(polarization expansion)라 한다. 1차에서 정전기 에너지가, 2차에서 유도와 분산이 나온다. 반데르발스 힘의 $-C_6/r^6$ 이 정확히 이 2차 분산 항이다.

문제는 이 전개가 파울리 원리를 위반한다는 것이다. 단순곱 $\Phi_A\Phi_B$ 는 A의 전자와 B의 전자를 맞바꾸는 데 대해 반대칭이 아니다. 그 결과 두 가지가 벌어진다.

교환반발이 통째로 없다. 분극 전개는 인력 성분만 내놓으므로, 두 분자가 가까워지면 에너지가 끝없이 내려간다. 실제 퍼텐셜 곡선의 벽이 생기지 않는다.
전개가 옳은 상태로 수렴하지 않는다. 분극 급수는 진짜 반대칭 바닥상태가 아니라 파울리 원리를 무시한 다른 상태로 수렴하며, 그 차이(교환 결손)는 급수 차수를 올려도 사라지지 않는다.¹

해법이 반대칭화 연산자 $\mathcal{A}$ 를 각 차수에 끼워 넣는 것이다. 즉 섭동 함수를 만들 때마다 A-B 사이 전자 치환을 포함해 반대칭화한다. 그 결과 각 분극 항에 짝이 되는 교환 항이 생긴다.

E_{\text{int}} = \underbrace{E_{\text{elst}} + E_{\text{ind}} + E_{\text{disp}}}_{\text{분극 성분}} + \underbrace{E_{\text{exch}} + E_{\text{exch-ind}} + E_{\text{exch-disp}}}_{\text{교환 짝}} + \delta_{\text{HF}}

교환 항들은 전자밀도가 겹치는 곳에서만 유의미하고 거리에 대해 지수적으로 감쇠한다. 그래서 멀리서는 SAPT가 고전적 다중극 전개와 일치하고, 가까이서는 벽이 제대로 선다.

4. 성분 읽는 법[편집]

성분	부호	물리	거리 의존
정전기 $E_{\text{elst}}$	±	얼어붙은 두 전자밀도 사이의 쿨롱	다중극 전개 꼴, 느림
교환반발 $E_{\text{exch}}$	+	파울리 원리에 의한 밀도 겹침 벌금	지수 감쇠
유도 $E_{\text{ind}}$	−	상대 전기장에 의한 분극 + 전하이동	다중극 꼴
교환유도	+	유도를 되돌리는 겹침 보정	지수 감쇠
분산 $E_{\text{disp}}$	−	요동 상관, 런던 힘	$-C_6/r^6$ 이하
교환분산	+	분산의 겹침 보정	지수 감쇠
$\delta_{\text{HF}}$	−	2차를 넘는 고차 유도(HF 수준)	—

$\delta_{\text{HF}}$ 항은 순수 섭동론적이 아니라 초분자 HF 상호작용 에너지에서 2차까지의 SAPT 유도·교환유도 항을 뺀 나머지로 정의된다. 유도는 극성 계에서 급수 수렴이 느리기 때문에 그 꼬리를 통째로 빌려오는 실용적 타협이며, 이온성이나 강한 수소결합 계에서는 반드시 켜야 한다.

해석할 때 두 가지를 조심해야 한다.

각 성분은 크고 총합은 작다. 수소결합 문서에 있는 물 이량체 분해가 전형적이다 — 정전기 −8, 교환 +7, 유도 −2, 분산 −2 남짓이 모여 총 −5 kcal/mol이 된다. 성분 하나를 10% 틀리면 총합은 20% 틀린다.
전하이동은 독립 성분이 아니다. SAPT에서 전하이동은 유도 안에 들어가 있고, 이를 분리하려는 시도는 기저 의존적이다. 기저를 키우면 “전하이동”이 줄고 “분극”이 커진다. “이 결합의 X%가 전하이동”이라는 문장을 볼 때는 항상 스킴과 기저를 확인해야 한다.

5. 계층과 비용[편집]

SAPT의 두 번째 섭동 축은 **단량체 내부의 전자 상관**이다. 상호작용 $V$ 에 대한 차수와 단량체 상관 $W$ 에 대한 차수를 각각 어디까지 가져가느냐로 계층이 정해진다.

SAPT0 — 단량체를 하트리-폭 수준으로 놓는다. 즉 $W$ 차수 0. 지배 비용은 분산 항이며 대략 $O(N^5)$ . 밀도적합(density fitting, DF)을 붙인 DF-SAPT0가 사실상의 실무 기본값이고, 100원자급 계도 감당한다. 정확도는 오차 상쇄에 상당히 의존하는데, 특히 jun-cc-pVDZ 같은 특정 기저와 짝지었을 때 성능이 좋다는 경험칙이 널리 쓰인다. 우연한 궁합이라는 점을 알고 써야 한다.²
SAPT2, SAPT2+, SAPT2+(3), SAPT2+3 — 단량체 상관을 차례로 더 넣는다. 비용이 $O(N^6)$ , $O(N^7)$ 으로 뛰며, 정확도는 CCSD(T)/CBS 참조값에 체계적으로 접근한다. 괄호 안 숫자는 $V$ 에 대한 3차 항을 포함하는지를 뜻한다. 여기에 결합 클러스터 이중여기(CCD)로 분산을 다듬은 변종들이 붙는다.
SAPT(DFT) / DFT-SAPT — 단량체를 밀도범함수이론으로 기술하고, 분산은 결합 콘-샴(coupled-KS) 방식으로 얻은 진동수 의존 밀도 감수율에서 계산한다. 단량체 상관을 DFT가 싸게 담아주므로 SAPT2+ 급 정확도를 훨씬 낮은 비용에 얻는 경우가 많다. 다만 콘-샴 궤도의 점근 거동이 나쁘면 유도·분산이 어긋나므로, 점근 보정된 범함수를 쓰고 이온화 퍼텐셜을 맞춰주는 작업이 필수다. 이 준비를 생략하면 조용히 틀린 답이 나온다.

즉 분산 보정이 DFT에 분산을 집어넣는 문제라면, SAPT(DFT)는 DFT를 단량체 기술에만 쓰고 분산은 처음부터 따로 정확히 계산하는 반대 방향의 전략이다.

6. 확장[편집]

밀도적합·자연궤도 — 4중심 적분을 3중심으로 근사해 SAPT0를 실용 영역으로 끌어내렸다. 오늘날 SAPT를 돌린다는 것은 대개 DF 구현을 돌린다는 뜻이다.
F-SAPT(functional-group SAPT) — 성분 에너지를 작용기 단위로 다시 분할한다. “이 결합 에너지의 −1.2 kcal/mol은 이 페닐기와 저 카르보닐 사이의 분산에서 왔다” 같은 문장을 만들어준다. 신약 프로젝트에서 구조-활성 관계를 물리적으로 설명할 때 강력하다.³
I-SAPT(intramolecular SAPT) — A와 B가 같은 분자의 두 조각이고 그 사이를 잇는 링커 C가 있을 때, 링커의 영향을 분리해 분자 내 상호작용을 SAPT로 다룬다. 분자 내 수소결합의 세기를 논할 때 쓰인다.
A-SAPT — 원자 단위 분할. 지도(map)처럼 시각화하기 좋다.
다체 SAPT 확장으로 삼체 항을 다루는 갈래도 있고, 물 모형의 다체 퍼텐셜이나 분극성 힘장 개발에서 참조값 생산에 쓰인다.

7. 무엇에 쓰나[편집]

힘장 파라미터화가 가장 직접적인 응용이다. 고전 힘장의 함수 형태 — 점전하 정전기, 레너드-존스 퍼텐셜의 반발과 $-r^{-6}$ , 유도 쌍극자 — 는 SAPT 성분과 일대일로 대응한다. 그러니 총 상호작용 에너지 하나에 전체를 맞추는 대신, 성분별로 따로 맞출 수 있다. 정전기 항은 SAPT 정전기 성분에, 반발 항은 교환 성분에, 분산 계수는 분산 성분에. 이렇게 하면 파라미터가 물리적 의미를 유지하고, 전이 가능성(transferability)이 좋아진다. AMOEBA류 다중극·분극 힘장과 최근의 물리 기반 힘장들이 이 경로를 따랐다.

분자 인식의 해석이 두 번째다. 분자 도킹 점수가 왜 그렇게 나왔는지, 리간드에 메틸기 하나를 붙였을 때 이득이 분산에서 왔는지 정전기에서 왔는지 — SAPT는 이 질문에 숫자로 답한다. 특히 F-SAPT를 쓰면 치환기 단위로 기여를 뽑을 수 있어, 의약화학자가 다음 합성을 정할 때 근거가 된다. 결합 자유에너지를 예측하는 도구는 아니지만(엔트로피와 용매화가 빠져 있다), 왜 세냐에 답하는 도구로는 대체재가 마땅치 않다.

벤치마크 참조값 생산도 빼놓을 수 없다. S22·S66 같은 집합은 CCSD(T)/CBS로 총 에너지를 정하지만, 그 집합을 성분별로 분석해 “이 부분집합은 분산 지배, 저 부분집합은 정전기 지배”라고 분류한 것이 SAPT다. 분산 보정 방법이 어떤 계에서 왜 실패하는지 진단할 때 이 분류가 지도 노릇을 한다.

8. 한계[편집]

급수는 점근적이다. SAPT 전개가 수렴한다는 보장은 없고, 실제로 짧은 거리에서는 발산 조짐이 나타난다. 그래서 SAPT는 평형 거리 부근과 그보다 먼 영역에서 신뢰하고, 강하게 겹치는 배치에는 적용하지 않는다.
단량체가 잘 정의되어야 한다. 화학 결합이 형성·절단되거나 전하가 크게 이동하는 계, 강한 정적 상관이 있는 계, 열린 껍질 계에서는 “두 단량체 + 약한 섭동”이라는 전제 자체가 무너진다. 공유결합을 SAPT로 분해하려는 시도는 무의미하다.
기저 문제는 사라지지 않는다. 초분자적 BSSE는 없지만, 각 성분의 기저 수렴은 여전히 느리다. 특히 분산은 완전기저 극한 수렴이 악명 높게 느려서, 확산 함수를 포함한 큰 기저나 명시적 상관법 계열의 도움이 필요하다.
비용. SAPT0조차 분자동역학 궤적에 얹기엔 비싸다. SAPT는 궤적을 돌리는 도구가 아니라 몇 개의 대표 배치를 정밀하게 해부하는 도구다.
성분 이름에 과몰입 금지. 분해는 스킴에 의존하는 개념적 구성물이지 관측량이 아니다. 같은 이량체를 다른 EDA로 쪼개면 다른 숫자가 나오고, 어느 쪽도 “틀린” 것이 아니다. 비교는 항상 같은 스킴 안에서.

9. 관련 문서[편집]

반데르발스 힘 · 분산 보정 · 수소결합 · 물 모형
전자 상관 · 하트리-폭 방법 · 결합 클러스터 방법 · 뫼러-플레세 섭동론
기저중첩오차 · 완전기저 극한 · 기저함수 · 명시적 상관법
밀도범함수이론 · 교환-상관 범함수 · 다중극 전개
힘장 · 레너드-존스 퍼텐셜 · 분자 도킹 · 분자동역학
에너지 분해 분석 · 분극성 힘장 · 결합 자유에너지

교과서적 시험대가 수소 원자 두 개와 헬륨 이량체다. 반대칭화 없이 분극 전개만 돌리면 반발 벽이 아예 서지 않아 두 원자가 가까워질수록 에너지가 계속 내려간다. 물질에 부피를 주는 것이 인력이 아니라 파울리 원리라는 사실을, 급수 하나를 잘못 짜는 것만으로 체험할 수 있다. ↩
SAPT0 + jun-cc-pVDZ 조합은 “오차 상쇄가 잘 되는 지점”으로 널리 쓰이지만, 기저를 키우면 오히려 나빠지는 구간이 있다. 기저를 키웠는데 결과가 나빠지면 보통은 뭔가 잘못된 것이지만, 여기서는 그게 정상이다. 이런 조합을 쓸 때는 “왜 그 기저인가”를 설명할 수 있어야 한다. ↩
F-SAPT 결과 그림은 리간드 위에 작용기별 기여가 색으로 칠해져 나오는데, 발표 자료로 대단히 잘 팔린다. 다만 그 색이 자유에너지가 아니라 고정 구조에서의 전자 상호작용 에너지라는 점은 슬라이드 아래 작은 글씨로만 남는 경우가 많다. ↩