DLPNO-CCSD(T) · ORCA 사용자 가이드

이 페이지에서 다룹니다

왜 DLPNO-CCSD(T)인가
기본 사용법
정밀도 단계 — Loose / Normal / Tight
완전 기저집합 외삽 (CBS)
F12 — 기저집합 한계에 빠르게 접근
열린 껍질 — UHF 레퍼런스
권장 워크플로우
한 걸음 더 — 다중참조 (CASSCF · NEVPT2)

왜 DLPNO-CCSD(T)인가

CCSD(T)는 단일참조 양자화학의 “황금 표준”이다. 웬만한 닫힌 껍질 화학종에서 상대 에너지를 ±1 kcal/mol(이른바 화학 정확도) 안으로 재현한다. 다만 표준 CCSD(T)는 비용이 시스템 크기에 대해 O(N⁷)로 늘어서, 원자 20개를 넘어가면 비현실적으로 비싸진다.

DLPNO(Domain-based Local Pair Natural Orbital) 근사가 이 비용을 확 줄여 준다. 국소 궤도와 도메인 기반 페어 자연 궤도를 써서 점근적으로 선형 확장(linear scaling)에 가까워지면서도, 기준 CCSD(T) 에너지의 약 99.9%를 재현한다. 매뉴얼은 DLPNO-CCSD(T)를 두고 “DLPNO-MP2가 되면 DLPNO-CCSD(T)도 된다”고 적어 뒀다. 실제로 원자 수백 개짜리 분자에도 쓸 수 있다.

기본 사용법

키워드 라인에 DLPNO-CCSD(T)와 보조 기저집합을 같이 적으면 된다.

# 닫힌 껍질 단일점 — 일반적인 정밀 에너지
! DLPNO-CCSD(T) cc-pVTZ cc-pVTZ/C TightSCF

* xyzfile 0 1 mol.xyz

보조 기저집합(/C 접미사)는 RI-MP2와 같은 형식이다. DLPNO-CCSD(T)는 RI 적분이 필수라 보조 기저집합을 반드시 같이 지정해야 한다.

주 기저집합	권장 보조 기저집합
`cc-pVDZ`	`cc-pVDZ/C`
`cc-pVTZ`	`cc-pVTZ/C`
`cc-pVQZ`	`cc-pVQZ/C`
`def2-SVP`	`def2-SVP/C`
`def2-TZVP`	`def2-TZVP/C`
`def2-TZVPP`	`def2-TZVPP/C`

정밀도 단계 — Loose / Normal / Tight

DLPNO의 컷오프 임계값은 세 단계로 묶여 있다.

키워드	전형적 비용	상대 정확도	권장 사용처
`LoosePNO`	0.5×	약 95~98%	대규모 스크리닝, 빠른 탐색.
`NormalPNO`	1× (기본)	약 99.5%	일반 연구. 기본 추천.
`TightPNO`	약 2~3×	약 99.9%	벤치마크, 발표 직전의 정밀 에너지.

# 벤치마크 수준 정확도
! DLPNO-CCSD(T) TightPNO cc-pVTZ cc-pVTZ/C TightSCF

상대 에너지 비교에서는 동일한 임계값을

반응 에너지나 결합 해리 에너지를 비교할 때는 모든 종에 같은 PNO 임계값을 써야 오차가 상쇄된다. 한 종은 NormalPNO로, 다른 종은 TightPNO로 계산하면 그 차이만 0.5 kcal/mol쯤 어긋날 수 있다.

완전 기저집합 외삽 (CBS)

CCSD(T)는 완전 기저집합 한계(CBS, Complete Basis Set limit)로 외삽할 때 가장 정확하다. ORCA의 Extrapolate 키워드가 이걸 자동으로 해 준다.

# cc-pVTZ와 cc-pVQZ로 CCSD(T) 단일점을 자동 외삽
! DLPNO-CCSD(T) Extrapolate(3/4,cc) AutoAux TightSCF

* xyzfile 0 1 mol.xyz

Extrapolate(X/Y,basis) 형식이다(콤마 뒤에 공백을 넣으면 파싱 오류로 죽는다). X, Y는 기저집합의 카디널 수(2 = DZ, 3 = TZ, 4 = QZ, 5 = 5Z)이고, basis는 cc(Dunning), def2, aug-cc 같은 기저집합 패밀리다. DLPNO와 함께 쓸 때는 AutoAux(보조 기저집합 자동 생성)를 반드시 붙여야 한다.

입력	실제 수행 작업
`Extrapolate(2/3,cc)`	cc-pVDZ와 cc-pVTZ의 두 계산을 자동 외삽.
`Extrapolate(3/4,cc)`	cc-pVTZ와 cc-pVQZ로 외삽. 가장 일반적.
`Extrapolate(2/3,def2)`	def2-SVP와 def2-TZVPP로 외삽.
`Extrapolate(3,cc)`	cc-pVDZ, cc-pVTZ, cc-pVQZ의 세 계산 + 두 쌍의 외삽.

EP2 · EP3 — 더 큰 기저집합을 위한 효율적 외삽

큰 기저집합(예: cc-pV5Z)에서 CCSD(T)를 직접 돌리기는 비싸도, MP2 수준에서는 되는 경우가 많다. EP2/EP3 외삽은 “큰 기저집합의 MP2 + 작은 기저집합의 CCSD(T) − MP2 보정”을 묶어 효율을 끌어올린다.

# EP2: 작은 기저집합에서 CC, 큰 기저집합에서 MP2 외삽
! DLPNO-CCSD(T) ExtrapolateEP2(2/3,cc) AutoAux TightSCF

F12 — 기저집합 한계에 빠르게 접근

F12 명시적 상관 방법은 더 작은 기저집합으로도 CBS 한계 근처 정확도를 내게 해 준다. DLPNO-CCSD(T)-F12는 이 둘을 합친 것이다.

# DLPNO-CCSD(T)-F12 with approach D — F12로 큰 기저집합 효과 대체
! DLPNO-CCSD(T)-F12D cc-pVDZ-F12 cc-pVDZ-F12-CABS cc-pVDZ-F12/C TightSCF

F12에는 전용 기저집합과 CABS(Complementary Auxiliary Basis Set)가 필요하다. 기저집합 이름 끝의 -F12, -F12-CABS를 빠뜨리지 않도록 한다. F12를 쓰면 cc-pVDZ-F12 한 번 계산이 cc-pVQZ나 cc-pV5Z CBS 외삽과 비슷한 정확도를 낸다.

열린 껍질 — UHF 레퍼런스

라디칼이나 다중항 시스템에서는 UHF 레퍼런스를 쓰는 “UHF-DLPNO-CCSD(T)”를 쓴다.

! UHF DLPNO-CCSD(T) cc-pVTZ cc-pVTZ/C TightSCF

* xyz 0 2          # 도플렛 라디칼
  O   0.0   0.0   0.000
  H   0.0   0.0   0.969
*

스핀 오염 확인 필수

UHF 레퍼런스를 쓸 때는 SCF 결과의 ⟨S²⟩ 값을 꼭 확인한다. 도플렛(S = 1/2)의 이상적인 ⟨S²⟩는 0.75인데, 스핀 오염이 심하면 1.0을 넘기기도 한다. 그러면 뒤따르는 CCSD(T) 결과도 믿기 어렵다. 오염이 심하면 ROHF 레퍼런스를 시도하거나, 다중참조 방법(CASSCF/NEVPT2)을 고려해야 한다.

권장 워크플로우

DLPNO-CCSD(T)는 보통 단일점 에너지로만 쓴다. 그래디언트도 구현돼 있지만 비싸서, 구조 최적화는 가벼운 방법(DFT)으로 하고 그 위에 정밀 에너지를 얹는 “composite” 워크플로우가 표준이다.

구조 최적화 + 진동수: B3LYP-D4/def2-TZVP 또는 r2SCAN-3c로 수행. 깁스 보정 ΔG_therm 확보.
DLPNO-CCSD(T) 단일점: 위 구조에서 CBS 외삽 또는 F12로 정밀 전자 에너지 E_elec^CC 산출.
용매 보정: 필요 시 (DFT/SMD) − (DFT/gas) 차이로 ΔG_solv 계산.
최종 자유에너지: G_final = E_elec^CC + ΔG_therm^DFT + ΔG_solv.

이 흐름은 입력 파일 세 개로 나눠 돌린다(한 입력의 $new_job으로 묶으면 안 되는 이유는 코드 아래 경고에).

# ── step1.inp ── DFT 최적화 + 진동수 (D4 분산 켬) → step1.xyz 생성
! B3LYP D4 def2-TZVP RIJCOSX def2/J TightOpt Freq
* xyzfile 0 1 mol.xyz

# ── step2.inp ── DLPNO-CCSD(T)/CBS 정밀 단일점 (분산 빼기 — CC가 이미 포함)
! DLPNO-CCSD(T) Extrapolate(3/4,cc) AutoAux TightSCF
* xyzfile 0 1 step1.xyz

# ── step3.inp ── 용매 보정 (SMD), 같은 최적화 구조에서
! B3LYP D4 def2-TZVP RIJCOSX def2/J CPCM(water) TightSCF
%cpcm
   SMD true
   SMDSolvent "water"
end
* xyzfile 0 1 step1.xyz

컴포지트에서 자주 죽는 두 지점 (ORCA 6.1.1 실측)

① 세 단계를 한 입력의 $new_job으로 묶지 말 것. DFT 최적화의 D4 분산 설정이 CC 단계로 넘어가 “dispersion + correlated method” 오류로 죽는다 (CCSD(T)는 분산을 이미 포함해 ORCA가 막는다). 그래서 파일을 나누고 CC 단계엔 D4를 넣지 않는다.
② DLPNO + Extrapolate에는 AutoAux가 필수이고(없으면 “보조 기저집합 필요” 오류), Extrapolate(3/4,cc)는 콤마 뒤 공백이 있으면 파싱 실패한다.

반응 평형이나 자유에너지 장벽을 보고할 때 쓰는 “DFT 최적화 + CCSD(T) 정밀점 + SMD 용매” 조합은 요즘 양자화학에서 가장 표준적인 절차이고, 학술지에서도 무난히 받아들여진다.

DLPNO를 다른 메소드와 결합

DLPNO 근사는 CCSD(T) 말고도 CCSD, MP2, B2PLYP(이중 하이브리드)와도 결합된다. 예를 들어 ! DLPNO-B2PLYP D4 def2-TZVP def2-TZVP/C로 이중 하이브리드 비용을 크게 낮출 수 있어, 중대형 분자의 정밀 단일점에 아주 쓸모 있다.

한 걸음 더 — 다중참조 (CASSCF · NEVPT2)

DLPNO-CCSD(T)는 단일참조 이론의 정점이다. 그런데 결합이 끊어지는 도중의 분자, 전이금속·란타나이드의 비공선 자기 상태, 단일항 디라디칼처럼 슬레이터 행렬식이 둘 이상 비슷한 가중치를 가지는 영역에서는 단일참조 가정 자체가 깨진다. 거기서부터는 활성 공간을 직접 푸는 CASSCF가 출발점이고, 그 위에 NEVPT2로 동적 상관을 얹는다.

활성 공간 잡는 법, state-averaging, NEVPT2, 스핀-궤도 결합·SINGLE_ANISO, 금속·리간드 조각 병합까지는 전용 챕터로 옮겼다 → 13. CASSCF · 다중참조.