전이 상태 · 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의 여러 변형

매뉴얼에 정리된 변형들은 비용·강건성에서 조금씩 다르다.

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를 그대로 발표하게 될 수 있다.