컴포지트 방법

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

1. 개요[편집]

컴포지트 방법 Composite Methods
분야	양자화학 × 계산열화학
발상	비싼 계산 하나 대신 싼 계산 여럿을 가법 합성
계열	Gaussian-n, CBS-QB3, ccCA, Weizmann-n, HEAT
목표	화학적 정확도 1 kcal/mol → 서브 kJ/mol
아킬레스건	가법성 가정, 경험적 HLC

CCSD(T)/CBS를 통째로 돌릴 돈이 없다면, 그 값을 여러 조각의 덧셈으로 추정하면 된다. 문제는 세상이 그렇게 덧셈으로 되어 있느냐는 것.

컴포지트 방법(composite method, 복합 방법)은 하나의 고수준·대기저 계산을 직접 수행하는 대신, 서로 다른 이론 수준과 기저에서 얻은 여러 개의 값싼 계산을 정해진 공식으로 합성해 고정밀 에너지를 추정하는 프로토콜의 총칭이다. 목표는 명확하다. 생성열·원자화 에너지·이온화 에너지·반응 장벽을 화학적 정확도(chemical accuracy, 1 kcal/mol ≈ 4.184 kJ/mol) 안에서 예측하는 것.¹

핵심 관찰은 이렇다. 결합 클러스터 방법의 비용은 $N^7$ 이고 기저는 $X^{-3}$ 로 기어가므로 CCSD(T)/CBS는 원자 열 개만 넘어도 손이 안 간다. 그런데 큰 기저 효과는 값싼 MP2에서 대부분 드러나고, 고차 상관 효과는 작은 기저에서도 대체로 옳게 나온다. 그렇다면 각각을 자기가 잘하는 조건에서 따로 계산해 더하면 된다는 것.

E \approx E[\text{저수준}/\text{대기저}] + \Delta E[\text{고수준} - \text{저수준}]_{\text{소기저}}

이 한 줄이 컴포지트 방법 전체의 문법이다. 나머지는 어떤 조각을 몇 개 쌓느냐의 차이일 뿐이다.²

2. 보정 항목 사전[편집]

고급 프로토콜의 최종 에너지는 대체로 다음 항들의 합이다. 각 항이 무엇을 메우는지를 알면 논문의 표를 읽을 수 있다.

항목	무엇을 메우나	대략적 크기(원자화 에너지 기준)
HF/CBS	참조 행렬식의 기저 불완전성	수 kcal/mol, 지수 외삽으로 쉽게 수렴
CCSD(T)/CBS 원자가 상관	전자 상관 본체	수십 kcal/mol, 최대 항목
고차 들뜸 (T→CCSDT, (Q), CCSDTQ)	(T) 근사의 잔여 오차	0.1~수 kcal/mol, 다중참조성 클수록 커짐
코어-밸런스 상관	frozen-core 근사가 버린 내각 상관	1분자당 0.5~수 kcal/mol, 2주기 원소에서 무시 못 함
스칼라 상대론 (DKH, MVD2)	내각 전자의 상대론적 수축	1주기·2주기는 작지만 3주기부터 급증
스핀-궤도 결합	원자·라디칼의 다중항 갈라짐	원자화 에너지에서 특히 중요(예: F, Cl 원자)
DBOC(대각 보른-오펜하이머 보정)	본-오펜하이머 근사 자체의 오차	대개 0.1 kcal/mol 이하, 수소 많으면 무시 못 함
ZPE·비조화 보정	영점 진동 에너지, 조화 근사 오차	수 kcal/mol, 정확도를 좌우하는 숨은 복병

여기서 자주 저평가되는 것이 마지막 두 줄이다. 전자 구조를 아무리 정밀하게 뽑아도 영점 에너지를 조화 진동수에 경험 계수만 곱해 때우면 그 오차가 그대로 최종값에 실린다. 서브 kJ급을 노리는 프로토콜이 진동 계산에 전자 구조만큼의 공을 들이는 이유다.

3. Gaussian-n 계열 — 대중화의 주역[편집]

1989년 포플(Pople) 그룹의 G1을 시작으로 G2(1991), G3(1998), G4(2007)로 이어진 계보다. 골격은 언제나 동일하다.

값싼 수준에서 구조 최적화와 진동수 계산(G4는 B3LYP/6-31G(2df,p)).
중간 수준의 기준 에너지(G3·G4는 MP4 또는 CCSD(T)/작은 기저).
확산 함수·고차 분극 함수·큰 기저의 효과를 각각 따로 계산한 증분으로 더함.
경험적 고준위 보정(higher-level correction, HLC)을 얹음.

성적표는 준수하다. G3가 G2/97 시험집합에서 평균절대편차 약 1.0 kcal/mol, G4가 더 큰 G3/05 시험집합(454개 에너지)에서 약 0.8 kcal/mol이다. 유기 소분자 열화학이라면 이 정도로 충분한 경우가 많고, 무엇보다 상용 코드에서 키워드 한 줄이라는 접근성이 압도적이었다.

3.1. HLC 논쟁[편집]

문제는 4번이다. HLC는 물리적 항이 아니라 원자가 전자 수만의 함수인 경험 파라미터다. 대략 이런 꼴이다.

E_{\text{HLC}} = -A\,n_\beta - B\,(n_\alpha - n_\beta)

$n_\alpha, n_\beta$ 는 짝수·홀수 스핀 원자가 전자 수이고, $A$ · $B$ 등의 상수는 실험 데이터에 맞춰 피팅한다. 비판은 세 갈래다.

ab initio가 아니다. 실험값에 맞춘 상수가 들어간 순간 “제1원리”라는 간판을 내려야 한다는 지적. Gn 지지자들은 “그럼 교환-상관 범함수에 파라미터 넣는 DFT는 뭐냐”고 반박하고, 논쟁은 30년째 무승부다.
구조에 무감각하다. HLC는 전자 개수만 보므로 이성질체 두 개에 똑같은 값이 들어가 정확히 상쇄된다. 즉 이성질화 에너지에는 아무 도움이 안 된다. 도움이 필요한 곳에서 도움이 안 되는 보정.
전자 짝지음이 바뀌면 어긋난다. 결합이 끊어져 라디칼 두 개가 되는 반응처럼 $n_\alpha - n_\beta$ 가 변하는 과정에서는, 복합체의 HLC와 조각들의 HLC 합이 일치하지 않아 크기 일관성 논란이 생긴다.

이 비판을 정면으로 받아 HLC를 아예 없앤 것이 윌슨 그룹의 ccCA(correlation consistent Composite Approach, 2006)다. 던닝의 상관일치 기저와 정직한 외삽만으로 구성해 경험 파라미터 없이 1 kcal/mol 근처를 찍었고, “가법성만으로도 여기까지 된다”는 것을 보여준 상징적 결과가 됐다.

4. CBS-QB3 — 다른 각도의 절약[편집]

피터슨(Petersson) 그룹의 CBS(complete basis set) 계열, 그중 CBS-QB3(1999)는 절약 지점을 다른 데서 찾았다. B3LYP로 구조·진동수를 잡고, MP2 쌍 에너지의 점근 거동을 이용해 기저 극한을 외삽한 뒤, 그 위에 MP4(SDQ)와 CCSD(T) 증분을 소기저에서 얹는다. 정확도는 G3급(평균절대편차 1 kcal/mol 남짓)이면서 계산이 가볍고, 특히 전이상태와 반응 장벽에서 평판이 좋아 반응속도론·전이상태 이론 쪽에서 오래 애용됐다.

5. Weizmann-n과 HEAT — 정밀도의 끝[편집]

경험 파라미터를 버리고 항을 끝까지 쌓아 올리는 노선이다. 마틴(Martin)·데올리베이라의 W1/W2(1999)와 그 후속 W3/W4, 그리고 스탠턴·게이지 그룹 등이 함께 만든 HEAT(High-accuracy Extrapolated Ab initio Thermochemistry, 2004)가 대표다.

W1 — HF·CCSD·(T) 성분을 각각 다른 기저 계열로 따로 외삽하고 코어-밸런스·상대론·스핀-궤도를 더한다. 원자화 에너지 평균절대편차 0.3~0.5 kcal/mol 수준.
W2 — 더 큰 기저, 더 엄격한 외삽. 경험 파라미터를 사실상 제거.
W4 — CCSDT, CCSDT(Q), 나아가 CCSDTQ5까지 쌓는다. 목표는 0.1 kcal/mol 대, 즉 서브 kJ/mol. 이 수준이 되면 실험 열화학 데이터를 거꾸로 교정하는 기준이 된다.
HEAT — 조성이 비슷하되 DBOC와 비조화 ZPE를 처음부터 정식 항목으로 포함시켰다. 목표 정확도 1 kJ/mol.

이 급으로 가면 계산 한 건에 원자 5~~10개가 한계이고, 실행 시간은 수백~~수천 코어시간을 잡아먹는다. 대신 ATcT(Active Thermochemical Tables) 같은 열화학 네트워크의 기준값 공급원 역할을 한다. 실험과 계산이 서로를 채점하는 관계가 여기서 역전되는 셈이다.³

6. 가법성 가정은 어디서 깨지는가[편집]

이 모든 레시피가 서 있는 유일한 가정은 “보정 항들이 서로 독립이다”, 즉

\Delta_{\text{core}}\!\left(\text{큰 기저}\right) \approx \Delta_{\text{core}}\!\left(\text{작은 기저}\right)

같은 식이 성립한다는 것이다. 교차항이 작다는 이 가정은 얌전한 닫힌껍질 분자에서는 잘 맞지만, 다음 경우에 깨진다.

다중참조성이 강한 계. 오존이 전설적인 사례다. 고차 들뜸 보정이 작기는커녕 수 kcal/mol로 커지고, 소기저에서 뽑은 증분이 대기저 값과 전혀 다르다. 컴포지트 방법이 “안 되는 계는 화끈하게 안 되는” 대표 지점.
전이금속·무거운 원소. 코어-밸런스와 스칼라 상대론이 커지고 서로 얽힌다. 애초에 Gn 계열의 파라미터가 주기율표 위쪽 유기 분자로 피팅됐다는 점도 겹친다.
구조가 저수준에서 이미 틀린 경우. 구조 최적화를 B3LYP로 하는 프로토콜에서 B3LYP가 엉뚱한 극소점을 찾으면, 그 위에 아무리 정밀한 단일점 에너지를 얹어도 소용없다. // 오른쪽이 무너지면 왼쪽은 의미가 없다.

그래서 실무 지침은 단순하다. 오차 막대를 보정 항의 크기로 가늠하라. 고차 들뜸 보정이 유난히 크게 나오면 그 계는 가법성이 흔들리고 있다는 신호이고, 그때는 컴포지트 값을 믿는 대신 다중기준 방법으로 갈아타야 한다. 반대로 각 보정이 얌전히 작으면 그 값은 대체로 광고된 정확도를 지킨다.

7. 현대적 변주[편집]

F12화. 참조 상관 에너지를 명시적 상관법으로 뽑는 W1-F12 / W2-F12 계열이 나오면서, 같은 정확도를 훨씬 작은 기저로 얻게 됐다. 컴포지트 방법의 가장 비싼 항이 반값이 된 셈.
초점점 분석(focal-point analysis). 이론 수준과 기저 크기를 2차원 표로 깔고 각 축 방향 증분을 따로 외삽하는 방식. 어느 항이 수렴했고 어느 항이 안 했는지가 표에서 눈으로 보인다는 점에서, 가법성 가정을 가정이 아니라 검증 대상으로 다룬다.
기계학습 대체. 대량의 컴포지트 계산 결과를 학습 데이터로 삼아 회귀 모형으로 생성열을 예측하는 시도가 활발하다. 다만 훈련 분포 밖으로 나가는 순간의 거동은 여전히 검증 및 확인의 숙제로 남아 있다.

8. 관련 문서[편집]

왜 하필 1 kcal/mol인가에 대한 물리적 근거는 사실 없다. 상온 아레니우스 인자로 환산하면 반응 속도상수가 대략 다섯 배쯤 어긋나는 수준인데, 실험 열화학 데이터의 전형적 오차 막대와 얼추 비슷하다는 실용적 이유가 크다. 즉 “이 정도면 실험과 싸울 자격이 있다”는 선이지, 자연이 그어준 선은 아니다. ↩
그래서 컴포지트 방법을 “이론”이 아니라 레시피라고 부르는 사람이 많다. 논문에 실린 것도 실제로 요리법에 가깝다 — 어떤 기저로, 어떤 순서로, 어떤 계수를 곱해 더하라. 재현성 하나는 확실하다는 게 미덕이자, 물리적 통찰이 별로 없다는 게 흠이다. ↩
실제로 W4·HEAT급 계산이 오래된 실험 생성열의 오류를 잡아낸 사례가 여럿 있다. “계산이 실험보다 정확하다”는 문장은 대개 허세지만, 원자 대여섯 개짜리 기체상 소분자의 원자화 에너지라는 아주 좁은 영역에서는 진심으로 참이다. ↩