전이 상태 · NEB

전이 상태 탐색의 세 가지 길

ORCA에서 전이 상태(TS, saddle point)를 찾는 길은 크게 셋이다.

NEB-TS: 처음 시도할 때

NEB-TS는 Nudged Elastic Band와 OptTS를 합친 방법이다. 반응물과 생성물 구조만 주면 알아서 둘을 잇는 경로를 만들고, 그 위에서 가장 높은 점을 골라 정밀하게 TS로 수렴시킨다.

아래는 아세트산 분자 내부 양성자 이동의 NEB-TS 입력 예다(매뉴얼 6.3.16절 예제).

# reactant.inp — 반응물 구조에서 시작
! XTB NEB-TS

%neb
   neb_end_xyzfile "product.xyz"    # 생성물 좌표 파일
end

* xyz 0 1
  C    0.416   0.039  -0.014
  C    0.042   0.012   1.440
  O    1.524   0.177  -0.454
  O   -0.654  -0.128  -0.804
  H   -0.391  -0.126  -1.737
  H   -0.913   0.507   1.585
  H   -0.058  -1.026   1.751
  H    0.820   0.485   2.030
*

product.xyz에는 같은 분자에서 “양성자가 다른 산소로 옮겨간” 구조의 좌표가 들어 있어야 한다. 원자 순서가 반응물과 똑같아야 한다는 점만 주의하면 된다.

NEB-TS의 결과는 이렇게 요약된다.

---------------------------------------------------------------
                         PATH SUMMARY
---------------------------------------------------------------
Image Dist.(Ang.)    E(Eh)        dE(kcal/mol)     max(|Fp|)
  0     0.000     -14.45993           0.00          0.00011
  1     0.426     -14.44891           6.91          0.00092
  2     0.652     -14.42864          19.63          0.00084
  3     0.805     -14.41132          30.50          0.00075
  4     0.932     -14.40562          34.08          0.00057   <= CI (전이 상태)
  5     1.044     -14.41047          31.03          0.00057
  ...
  9     1.869     -14.45988           0.03          0.00013

“CI”(Climbing Image)로 표시된 이미지가 수렴된 전이 상태이고, 그 옆 dE 값(kcal/mol)이 활성화 에너지의 1차 추정치다. NEB-TS는 이 점을 OptTS로 한 번 더 정밀 수렴시키니, 최종 결과는 출력 파일 끝의 FINAL SINGLE POINT ENERGY를 보면 된다.

NEB의 여러 변형

NEB 변형은 비용과 강건성이 조금씩 다르다.

OptTS (초기 추정이 좋을 때)

이미 TS에 가까운 초기 구조와 정밀한 헤시안이 있으면 OptTS가 가장 효율적이다. 단, 헤시안은 반드시 같은 메소드의 진동수 계산으로 먼저 만들어 둬야 한다.

# 1단계: 초기 추정 구조에서 헤시안 계산
! B3LYP D4 def2-SVP NumFreq
* xyzfile 0 1 ts_guess.xyz

$new_job

# 2단계: OptTS
! B3LYP D4 def2-SVP OptTS Freq
%geom
   InHess     Read
   InHessName "job1.hess"
   Calc_Hess  true          # 매 5 step마다 헤시안 재계산
   Recalc_Hess 5
end
* xyzfile 0 1

ScanTS: 스캔과 TS 검색 자동화

반응 좌표가 결합 길이나 각도 하나로 분명하면 ScanTS가 직관적이다. ORCA가 표면을 스캔하다가 최대값을 지나면, 그 점을 초기 추정 삼아 OptTS를 시작한다.

! B3LYP D4 def2-SVP ScanTS

%geom Scan
   B 2 5 = 2.5, 1.0, 16      # 2-5 결합을 2.5 → 1.0 Å로 스캔
end
end

* xyzfile 0 1 reactant.xyz

스캔이 최대를 지나는 순간 자동으로 OptTS로 넘어간다. 스캔 중간에 헤시안을 한 번 계산해 두면 정확도가 올라간다. 전체 스캔이 끝날 때까지 기다리게 하고 싶으면 fullScan true를 더한다.

IRC로 전이 상태 검증

TS를 얻었으면 그게 정말 의도한 반응을 잇는 TS인지 확인해야 한다. IRC(Intrinsic Reaction Coordinate)는 TS에서 양쪽으로 가장 가파른 내리막을 따라 내려가 반응물과 생성물에 닿는 경로를 추적한다.

! B3LYP D4 def2-SVP IRC

%irc
   MaxIter   100
   PrintLevel 1
   Direction  both    # forward / backward / both
   InitHess   read
   Hess_Filename "ts.hess"
end

* xyzfile 0 1 ts.xyz

IRC는 같은 입력 파일에서 !OptTS, !ScanTS, !NEB-TS, !AnFreq, !NumFreq와 같이 쓸 수 있어, “TS 검색 → 진동수 확인 → IRC”를 한 번에 끝낼 수도 있다.

# NEB-TS → AnFreq → IRC를 한 번에
! B3LYP D4 def2-SVP NEB-TS AnFreq IRC

%neb
   neb_end_xyzfile "product.xyz"
end

* xyzfile 0 1 reactant.xyz

권장 워크플로우

처음 보는 반응의 TS를 찾을 때는 보통 이 순서로 간다.

1. 반응물·생성물 구조 최적화: 가벼운 방법(예: r2SCAN-3c 또는 XTB)으로.
2. XTB-NEB-TS: 대략적인 TS를 빠르게 잡는다. 보통 몇 분이면 끝난다.
3. DFT OptTS: 위에서 얻은 TS 구조를 초기값으로 정밀 수렴. 헤시안은 함께 계산.
4. AnFreq: 허수 진동수가 정확히 하나인지 확인.
5. IRC: 양방향 IRC로 반응물·생성물에 도달하는지 검증.
6. 단일점 정밀화: 필요하면 DLPNO-CCSD(T)/cc-pVTZ 등으로 최종 에너지를 다듬는다.

IRC를 빼먹으면 “전이 상태처럼 보이지만 사실은 다른 반응을 잇는” 가짜 TS를 그대로 발표하게 될 수 있다.