Running head: CATALYSIS IN DIELS–ALDER REACTION

Computational Investigation of Catalytic Effects in the Diels–Alder Reaction Using ORCA and Avogadro

LJH

Hansung Science High School, Department of Chemistry

본 연구는 한성과학고등학교 화학 자율연구의 일환으로 수행되었으며, 외부 지도교수의 직접적인 지도 없이 진행되었다. Correspondence to bigrirason@gmail.com.

Chemistry Route

Room: dft · Type: paper · Status: active
Parent: Chemistry MOC · Study track: Study MOC
Why here: 계산화학 방법, 전이상태 탐색, 활성화 자유에너지 비교가 핵심이라 DFT room에 둔다. 유기 반응 해석은 related_rooms: [organic]으로 연결한다.
Next action: 계산 수준과 reference state를 통일해 ΔG‡를 재계산하고, 각 TS를 IRC로 검증한다.

Abstract

본 연구는 cyclopentadiene과 acrolein 사이의 Diels–Alder 반응에 대해 여러 촉매 조건의 전이상태와 겉보기 활성화 자유에너지(ΔG‡)를 ORCA와 Avogadro를 이용해 비교하였다. 무촉매, 단순화한 iminium 활성화, AlCl₃ Lewis acid, η² Cu⁺–diene 및 η⁴ Fe(CO)₃–diene 조건에서 OptTS와 진동수 계산을 수행했다. 최종 출력 기준 겉보기 ΔG‡는 무촉매 31.87, iminium 6.04, AlCl₃ 9.15, Cu⁺–diene 21.26, Fe(CO)₃–diene 61.54 kcal/mol이었다. diene 결합형 계에서는 η² 조건의 장벽이 가장 낮고 η⁴ 조건의 장벽이 가장 높아, hapticity 증가에 따른 단조적 활성 증가 가설과 다른 비단조 경향을 보였다. 다만 금속 종류·전하·리간드와 reference state가 동시에 달라 이 결과를 hapticity만의 인과 효과나 촉매 간 절대적 우열로 단정할 수 없다. 현재 결과는 dienophile의 frontier orbital 안정화, 금속–diene 결합 양식 및 접근 기하를 함께 검토해야 함을 시사한다.

Keywords: Diels–Alder, ORCA, 전이상태 이론, Lewis acid 촉매, iminium 유기촉매, hapticity, HOMO–LUMO 상호작용

---

Introduction

Diels–Alder 반응은 conjugated diene과 dienophile 사이의 [4+2] cycloaddition으로, frontier molecular orbital(FMO) 이론, endo/exo 선택성, 그리고 thermodynamic·kinetic product의 경쟁 등 유기화학의 핵심 개념들이 집약된 모델 반응이다. 특히 MacMillan @MacMillan2000 의 first highly enantioselective organocatalytic Diels–Alder 연구 이후, iminium 활성화 기반의 유기촉매와 Lewis acid 촉매가 dienophile의 LUMO를 안정화시켜 반응 속도를 향상시킨다는 사실이 광범위하게 입증되어 왔다.

그러나 대부분의 교과서적 논의는 dienophile 활성화에 집중되어 있으며, 반대로 diene에 금속이 결합한 경우에는 diene의 π 전자 분포와 접근 기하가 어떻게 바뀌는지를 함께 고려해야 한다. Hapticity(ηⁿ)는 금속에 연속적으로 배위하는 리간드 원자의 수를 뜻하며, 금속–diene 착물에서는 η²와 η⁴ 같은 결합 양식이 가능하다 @hapticity_review . 본 연구는 이 결합 양식의 차이가 계산된 전이상태와 활성화 자유에너지에 어떤 경향을 만드는지 탐색했다.

본 연구는 cyclopentadiene과 acrolein의 표준 Diels–Alder 반응을 기준으로, 다음의 다섯 조건을 ORCA를 이용해 계산하고 비교한다: (1) 무촉매, (2) iminium 유기촉매, (3) AlCl₃ Lewis acid (acrolein 활성화), (4) Fe(CO)₃–cyclopentadiene 착물 (η⁴ diene 활성화), (5) Cu⁺–cyclopentadiene 착물 (η² diene 활성화).

연구 질문

1. RQ1: 온도 변화에 따른 열보정은 Diels–Alder 반응의 ΔG‡에 어떤 경향을 만드는가?
2. RQ2: acrolein 대신 cyclopentadiene 측에 Lewis acid (또는 금속 착물)를 결합시키는 전략이 dienophile 활성화 전략에 필적하는 촉매 효과를 낼 수 있는가?

가설

• H1: Diels–Alder 반응은 두 분자가 하나의 고리형 전이상태를 형성하므로 병진·회전 자유도 감소로 인해 엔트로피적으로 불리하며, 온도가 높아질수록 ΔG‡의 열보정 항이 증가할 것이다.
• H2: Orbital overlap의 재배열과 reaction plane 고정 효과를 고려할 때, 금속-diene 착물의 hapticity가 클수록 ΔG‡가 더 크게 감소할 것이다.

Method

Materials (Computational Setup)

모든 계산은 다음 도구로 수행되었다:

• 분자 구조 생성·시각화: Avogadro
• 양자화학 계산: ORCA 5.0.2
• 사용 장비: 삼성 갤럭시북 4 (14인치 노트북)

구조 최적화와 진동수 계산의 기본 수준은 B3LYP/def2-SVP였다. 촉매 반응 복합체와 금속 착물 계산 다수에는 D3BJ dispersion correction, def2/J 및 RIJCOSX가 사용되었다. 반면 acrolein·cyclopentadiene 단독 구조와 무촉매 TS에는 D3BJ가 포함되지 않았다. 또한 입력 파일에서 연속체 용매 모형은 확인되지 않았으므로 결과는 기체상 계산으로 해석했다. 이 불일치는 현재 결과를 탐색적·정성적 비교로 제한하는 요인이다.

Procedure

1. 기준 반응계 선정

계산 자원 한계를 고려하여 비교적 작은 분자로 구성된 Diels–Alder 반응계를 선정하였다. Diene으로 cyclopentadiene, dienophile로 acrolein을 사용하였다.

2. 반응물 단독 최적화 및 기본 계산

Avogadro에서 각 반응물의 초기 구조를 제작한 뒤 ORCA로 구조 최적화 및 진동수 계산을 수행하여 Gibbs 자유에너지, 엔탈피, 엔트로피 보정값, HOMO/LUMO 에너지를 얻었다. 활성화 자유에너지의 기본 정의는 다음과 같다:

$$\Delta G^{‡}=G_{\text{TS}}-(G_{\text{cyclopentadiene}}+G_{\text{acrolein}})$$

3. 무촉매 전이상태 탐색 (OptTS)

다음 기준으로 전이상태를 검증하였다:

1. Imaginary frequency가 정확히 1개일 것.
2. 해당 imaginary mode가 두 C–C 결합 형성 방향과 일치할 것.
3. 분자 회전·병진·메틸 회전 등은 거짓 TS로 분류.
4. 형성 중인 두 C–C 결합 거리가 합리적 범위에 있을 것.

4. Iminium 촉매 모델

MacMillan의 유기촉매 전략에서 핵심이 되는 iminium 활성화 원리를 따라 acrolein-derived iminium 모델을 구성하였다. 이 모델은 실제 키랄 MacMillan 촉매, 반대이온 및 용매를 포함한 실험계를 그대로 재현한 것이 아니라, dienophile의 전자받개 오비탈이 안정화되는 원리를 조사하기 위한 단순화 모델이다. 활성화 자유에너지는 reactant complex를 기준으로 계산하였다:

$$\Delta G_{\text{iminium}}^{‡}=G_{\text{TS,iminium}}-G_{\text{reactant complex}}$$

5. Lewis Acid 촉매 모델

Carbonyl group이 hard Lewis base임을 고려하여 hard Lewis acid 위주로 후보를 선정하였다.

촉매 후보	선정 이유	예상 효과
BF₃	대표적 Lewis acid 기준 모델	acrolein LUMO 안정화
BH₃	비교군 (약한 Lewis acid)	모델 비교
AlCl₃	강한 Lewis acid	가장 큰 LUMO 안정화 예상

실제 계산은 AlCl₃–acrolein 착물 중심으로 수행하였다:

$$\Delta G_{{\text{AlCl}}_3}^{‡}=G_{{\text{TS,AlCl}}_3}-G_{{\text{AlCl}}_3\text{-acrolein complex}}$$

6. Diene-결합형 금속 촉매 모델

Diene에 결합한 금속이 반응 평면 고정과 orbital overlap 재배열을 통해 촉매 효과를 낼 수 있다는 가설 하에 다음 후보를 설정하였다.

금속 촉매 후보	예상 hapticity	선정 이유
Cu⁺	η²	diene 한쪽 π결합과 상호작용
Ag⁺	η²	Cu⁺와 유사하나 더 soft한 금속 중심
Fe(CO)₃	η⁴	cyclopentadiene 전체 π계와 결합

Fe(CO)₃–cyclopentadiene 착물에서 OptTS 계산 시 초기에는 imaginary frequency가 매우 작고 형성 중 C–C 결합이 3 Å 이상으로 나타나 진정한 TS로 보기 어려웠다. Acrolein의 접근 방향과 형성될 두 C–C 결합 거리의 비대칭성을 조정하여 더 큰 imaginary frequency를 가지는 TS 후보를 확보하였다.

Cu⁺ 착물의 경우 생성물에 지나치게 유사한 구조가 형성되거나 수렴 실패가 발생하여, 형성될 C–C 결합 거리를 2.6–2.7 Å 수준으로 늘린 constrained 구조를 먼저 만든 뒤 OptTS를 수행하였다.

$$\Delta G_{{\text{Fe(CO)}}_3}^{‡}=G_{{\text{TS,Fe(CO)}}_3}-(G_{{\text{Fe(CO)}}_3\text{-Cp complex}}+G_{\text{acrolein}})$$ $$\Delta G_{{\text{Cu}}^{+}}^{‡}=G_{{\text{TS,Cu}}^{+}}-(G_{{\text{Cu}}^{+}\text{-Cp complex}}+G_{\text{acrolein}})$$

Data Analysis

모든 ΔG‡ 값은 Hartree(Eh) 단위에서 계산한 뒤 1 Eh = 627.509 kcal/mol을 이용하여 환산하였다. 298.15 K의 값은 각 최종 ORCA 출력의 Gibbs 자유에너지로 계산했다. 온도 의존성은 273.15–1073.15 K, 1 atm 조건의 thermal correction 파일을 이용해 별도로 검토했다.

변인 통제 주의

본 연구에서 무촉매 조건은 분리 반응물 합 기준, iminium·Lewis acid 조건은 reactant complex 기준으로 계산되었다. Cu⁺·Fe(CO)₃ 조건은 metal–diene complex와 분리된 acrolein의 합을 기준으로 했다. 따라서 절대값 비교에는 reference state mismatch가 존재하며, 후속 연구에서는 모든 조건을 동일한 분리 반응물 기준으로 재계산할 필요가 있다.

온도 스윕 CSV의 범위

thermo_sweep_complex_to_TS.csv와 thermo_sweep_acrolein_only.csv의 Eel_* 열은 모두 0이다. 현재 CSV의 dG_activation은 전체 활성화 자유에너지가 아니라 열보정 차이이므로, 전자에너지 차이를 더하지 않은 값을 완전한 ΔG‡로 해석해서는 안 된다.

Results

4.1 반응물 및 촉매 Complex의 계산 결과

Table 1. 주요 반응물 및 complex의 ORCA 계산 결과

구조	Gibbs 자유에너지 / Eh	HOMO / Eh	LUMO / Eh	HOMO–LUMO gap / eV	비고
Acrolein	−191.62627218	−0.259295	−0.069046	5.177	dienophile
Cyclopentadiene	−193.76577644	−0.217834	−0.018402	5.427	diene
Iminium (acrolein-derived)	−250.57600236	−0.457898	−0.264469	5.263	단순화한 활성화 모델
Iminium–Cp reactant complex	−444.37849800	−0.360523	−0.231852	3.501	반응 복합체
AlCl₃–acrolein–Cp reactant complex	−2008.02977059	−0.253753	−0.120536	3.625	Lewis acid 반응 복합체
Fe(CO)₃–Cp complex	−1796.91550073	−0.217460	−0.028285	5.148	η⁴ diene complex
Cu⁺–Cp complex	−1833.77229864	−0.420072	−0.272939	4.004	η² diene complex

주. Frontier orbital 값은 각 최종 ORCA 출력의 마지막 ORBITAL ENERGIES 블록에서 읽었다. 특히 이전 표의 Cu⁺–Cp HOMO/LUMO 값은 Fe(CO)₃–Cp 값이 중복 입력되어 있었으므로 수정했다.

Iminium 및 AlCl₃ 반응 복합체의 전체계 gap은 각각 3.501, 3.625 eV로 Fe(CO)₃–Cp 복합체의 5.148 eV보다 작다. 이는 dienophile을 직접 활성화한 두 조건의 낮은 내부 장벽과 정성적으로 일치한다. 그러나 전체계 HOMO와 LUMO가 실제 diene π와 dienophile π*에 얼마나 국소화되어 있는지는 이 수치만으로 알 수 없다. 따라서 gap은 보조 지표로 사용하고, 오비탈 등위면과 결합 분석을 함께 확인해야 한다.

4.2 전이상태 및 활성화 자유에너지 비교

Table 2. 촉매 조건에 따른 전이상태 Gibbs 자유에너지와 활성화 자유에너지

조건	G_TS / Eh	기준 반응물 G / Eh	ΔG‡ / Eh	ΔG‡ / kcal·mol⁻¹	해석
무촉매 DA	−385.34126249	−385.39204862	0.05078613	31.87	기준 반응
Iminium 촉매	−444.36886745	−444.37849800	0.00963055	6.04	가장 큰 촉매 효과
AlCl₃–acrolein	−2008.01519377	−2008.02977059	0.01457682	9.15	LUMO 활성화
Fe(CO)₃–diene	−1988.44369654	−1988.54177291	0.09807637	61.54	오히려 불리
Cu⁺–diene	−2025.36469866	−2025.39857082	0.03387216	21.26	무촉매보다 유리, 그러나 약함

주. kcal/mol 값은 1 Eh = 627.509 kcal/mol로 환산하였다. Iminium과 AlCl₃ 값은 reactant complex 기준의 내부 장벽이며, 다른 행과 reference state가 같지 않다.

4.3 전이상태 검증

모든 최종 TS 출력은 정상 종료·수렴했으며, 최종 진동수 계산에서 imaginary frequency가 정확히 하나였다.

조건	최종 imaginary frequency / cm⁻¹	형성 중인 두 C···C 거리 / Å(근삿값)
무촉매	−455.45	2.07, 2.45
Iminium	−378.15	2.00, 약 2.85
AlCl₃	−400.73	2.00, 약 2.74
Fe(CO)₃–diene	−378.58	1.69, 2.53
Cu⁺–diene	−386.31	2.15, 2.37

Avogadro에서 저장된 normal-mode 궤적을 확인한 결과, 허수 진동 모드는 Diels–Alder 결합 형성 좌표와 관련된 움직임을 보였다. 다만 imaginary frequency 하나와 시각적 mode 확인만으로 양쪽 끝점이 완전히 검증되는 것은 아니므로, 최종 판정에는 양방향 IRC가 필요하다.

4.4 온도 의존성 (RQ1)

273.15–1073.15 K 범위에서 CSV에 기록된 TS–reference 열보정 차이는 온도에 따라 증가했다. 그러나 두 온도 스윕 CSV의 전자에너지 열이 모두 0이므로, 현재 기록된 298.15 K 값인 14.16 kcal/mol(무촉매)과 1.59 kcal/mol(iminium)은 완전한 ΔG‡가 아니다. 전자에너지를 포함한 최종 출력 기반 298.15 K 장벽은 각각 31.87 및 6.04 kcal/mol이다. 따라서 H1의 열보정 경향은 잠정적으로 지지되지만, 완전한 G(T)를 재구성하기 전에는 촉매 조건별 온도 의존성이나 고온에서의 촉매 기능 약화를 결론내릴 수 없다.

4.5 Hapticity 효과 (RQ2)

최적화 구조에서 Fe–C(Cp) 거리는 2.671, 2.126, 2.066, 2.067, 2.130 Å였다. 연속한 네 탄소가 약 2.07–2.13 Å로 가깝고 나머지 한 탄소가 2.67 Å로 멀어 η⁴ 배정과 잘 맞는다. Cu–C(Cp) 거리는 3.353, 2.997, 2.185, 2.069, 2.955 Å로, 연속한 두 탄소만 2.07–2.19 Å로 가까워 η² 배정과 잘 맞는다. 즉, hapticity를 단순 명칭이 아니라 연속한 배위 탄소 수와 금속–탄소 거리로 조작적으로 구분했다.

본 연구의 MO 기반 hapticity 정의

별도 연구 노트 hapticity.pdf에서 제안한 정의에 따라, complex $F=[L-M]$를 금속 fragment $L$과 ligand fragment $M$으로 분해한다. Complex의 $i$번째 MO는 fragment MO의 선형결합으로 쓴다.

$${\Psi }_i=\sum _{a\in N_L}{c}_{ai}{\psi }_a^L+\sum _{b\in N_M}{c}_{bi}{\psi }_b^M$$

각 MO에서 두 fragment의 weight는 다음과 같다.

$$w_i^L=\sum _a|c_{ai}{|}^2,w_i^M=\sum _b|c_{bi}{|}^2,w_i^L+w_i^M=1$$

$w_i^L{w}_i^M\ge \kappa$를 만족해 두 fragment가 충분히 섞이고, fragment orbital의 가중 에너지에 비해 complex MO가 안정화되어

$$\Delta E_i={\epsilon }_i-\sum _a|c_{ai}{|}^2{\epsilon }_a^L-\sum _b|c_{bi}{|}^2{\epsilon }_b^M<0$$

인 MO의 집합을 bonding states $\mathcal{B}$로 선택한다. 금속 원자 $A$의 AO $\mu$와 ligand 원자 $v$의 AO $\nu$ 사이에서 $i$번째 MO가 만드는 overlap contribution을

$$q_{iAv}=\sum _{\mu \in A{O}_A}\sum _{\nu \in A{O}_v}{c}_{\mu i}^{\ast }{c}_{\nu i}{S}_{\mu \nu }$$

로 두면, 본 연구에서 정의한 연속형 hapticity $n_A$는 점유수 $n_i$로 가중한 금속–ligand overlap population이다.

$$n_A=\sum _{i\in \mathcal{B}}{n}_i\sum _{v\in M}{q}_{iAv}$$

이 정의에서 표준 η²·η⁴는 연속적으로 결합한 원자의 위상적 개수이고, $n_A$는 그 결합에 참여하는 MO의 혼합·안정화·점유수·AO overlap을 반영하는 연속적 전자 hapticity다. 따라서 같은 η⁴ 구조라도 네 M–C 접촉의 세기가 약하면 $n_A$는 작아질 수 있고, 반응 도중 ring slippage가 일어나면 $n_A$가 연속적으로 변할 수 있다.

추가 ORCA fragment-MO 계산과 overlap hapticity

기존 최적화 구조를 고정하고 같은 수준(B3LYP-D3BJ/def2-SVP, def2/J, RIJCOSX)에서 Fe 및 Cu의 reactant complex와 TS에 대해 fragment-resolved single-point 계산을 추가했다. ORCA의 LOEWDIN FRAGMENT POPULATIONS PER MO에서 각 MO의 fragment weight를, MULLIKEN FRAGMENT-FRAGMENT OVERLAP PER MO에서 금속–diene AO overlap population을 추출했다. ORCA가 인쇄하는 fragment-pair 값 $O_i^{A-M}$은 두 AO 순서를 합한 cross term이므로, 위에서 정의한 한쪽 방향의 $q_{iAv}$와는 다음 관계를 가진다.

$$q_i^{A-M}=\sum _{v\in M}{q}_{iAv},O_i^{A-M}=2q_i^{A-M}$$

에 직접 대응한다. Cu 계에서는 $L={\mathrm{Cu}}^{+}$, $M=\mathrm{Cp}$로 나누었다. Fe 계의 MO 혼합도는 보조 리간드의 전자 효과를 포함하도록 $L=\mathrm{Fe}(\mathrm{CO}{)}_3$, $M=\mathrm{Cp}$로 계산하되, $q_i^{A-M}$는 별도의 원자 fragment 출력에서 Fe 원자–Cp 사이의 overlap만 취해 CO–Cp 중첩이 hapticity에 섞이지 않게 했다. TS에서는 acrolein을 세 번째 spectator fragment로 분리했다.

hapticity.pdf에는 $\kappa$의 수치가 지정되어 있지 않으므로 결과를 본 뒤 조정하지 않고 $\kappa =0.05$로 고정했다. 또한 ORCA의 population 출력은 complex MO와 isolated-fragment MO 사이의 일대일 energy mapping을 제공하지 않으므로, 원래의 $\Delta E_i<0$ 조건을 임의의 fragment energy로 대체하지 않았다. 대신 혼합 조건을 통과한 occupied MO 중 $O_i^{A-M}>0$인 궤도를 결합성-overlap 집합 ${\mathcal{B}}_{+}(\kappa )$로 정의하고, 다음의 overlap hapticity를 계산했다.

$$n_A^{(+)}(\kappa )=\sum _{i\in {\mathcal{B}}_{+}(\kappa )}{n}_i{q}_i^{A-M}=\sum _{i\in {\mathcal{B}}_{+}(\kappa )}\frac{n_i}{2}{O}_i^{A-M},{\mathcal{B}}_{+}(\kappa )=\{i\mid n_i>0,w_i^L{w}_i^M\ge \kappa ,O_i^{A-M}>0\}$$

모든 분석 대상 MO는 closed-shell occupied orbital로 $n_i=2$이므로, 각 MO의 $n_i{q}_i^{A-M}$는 ORCA가 인쇄한 $O_i^{A-M}$와 수치적으로 같다. 이 값은 금속–diene AO overlap을 MO별로 분해하여 다시 합한 수치이므로 결합차수의 대용값이 아니다. 다만 $\Delta E_i<0$를 isolated-fragment projection으로 판정한 원래 $n_A$와는 구분하여 $n_A^{(+)}$로 표기한다.

계	표준 hapticity	$n_A^{(+)}$, reactant	$n_A^{(+)}$, TS	TS retention	겉보기 ΔG‡ / kcal·mol⁻¹
Cu⁺–diene	η²	0.1024	0.0723	70.6%	21.26
Fe(CO)₃–diene	η⁴	0.5106	0.2112	41.4%	61.54

추가 계산 파일

입력과 출력은 C:\DFT_work\orcaout\sp-bestreaction\hapticity_mo_analysis에 저장했다. 분석에 사용한 출력은 cu_complex_frag.out, cu_ts_frag.out, fe_complex_frag.out, fe_ts_frag.out, fe_complex_ligandweights.out, fe_ts_ligandweights.out이며 모두 ORCA TERMINATED NORMALLY를 확인했다. Fragment label만 바꾼 Fe 계산 사이의 최종 에너지 차이는 $3\times 10^{-8}$ Eh 이하로, 동일한 수렴 해를 재현했다.

Reactant complex에서 Fe의 $n_A^{(+)}$는 Cu의 약 4.99배로, η⁴ Fe 계가 η² Cu 계보다 훨씬 넓고 강한 금속–diene orbital overlap을 가진다는 구조적 배정과 일치한다. 그러나 Fe 계의 장벽은 오히려 40.28 kcal/mol 더 높다. 또한 TS로 이동할 때 overlap hapticity는 Fe에서 0.2994(58.6%) 감소하지만 Cu에서는 0.0301(29.4%)만 감소한다. 따라서 Fe 계의 높은 장벽은 큰 초기 hapticity 자체가 촉매 활성을 높이기보다, 새 C–C 결합을 만들기 위해 강한 다점 overlap을 해제하고 재조직하는 비용을 증가시킨 결과로 해석할 수 있다.

금속 없음(형식적 η⁰), Cu⁺, Fe(CO)₃ 순서로 reactant의 $n_A^{(+)}$는 0, 0.1024, 0.5106이지만 겉보기 ΔG‡는 31.87, 21.26, 61.54 kcal/mol이다. 따라서 H2(hapticity ↑ → ΔG‡ ↓)는 지지되지 않았고, 현재 세 조건에서는 중간 orbital overlap에서 최소 장벽이 나타나는 비단조 경향이 확인되었다. 다만 Cu⁺와 Fe(CO)₃는 금속 종류, 전하, 보조 리간드 및 전자 수가 동시에 다르므로 이를 hapticity만의 인과 효과로 단정할 수 없다.

4.6 전자상태와 리간드 결합 분석

최적화가 끝난 구조의 전하·orbital energy를 다시 확인하고, 추가 single-point 출력의 fragment population과 MO별 fragment overlap을 결합해 분석했다. 원자 전하와 d population은 basis set에 의존하므로 산화수를 직접 뜻하지 않으며, 같은 계산 수준에서 reactant complex와 TS 사이의 상대적 변화만 해석했다.

4.6.1 금속 중심의 전자상태

계	ORCA 입력 전하	다중도	형식적 전자상태	금속 Mulliken charge: complex → TS	금속 d population: complex → TS
Fe(CO)₃–diene	0	1	Fe(0), d⁸, 18-electron closed shell	−0.794 → −0.418	6.676 → 6.594
Cu⁺–diene	+1	1	Cu(I), d¹⁰ closed shell	+0.574 → +0.565	9.828 → 9.845
AlCl₃–acrolein	0	1	Al(III) Lewis acid, closed shell	+0.726 → +0.738	해당 없음
Iminium	+1	1	closed-shell iminium cation	금속 없음	해당 없음

Fe(CO)₃–diene은 중성 singlet이며, 중성 diene을 4-electron donor, 세 CO를 각각 2-electron donor로 세는 형식적 18-electron Fe(0) 착물이다. ORCA의 Fe 원자 Mulliken charge가 음수인 것은 형식적 산화수가 음수라는 뜻이 아니라, Fe–CO 및 Fe–diene 결합의 강한 공유성과 ligand-to-metal donation/metal-to-ligand back-donation이 Mulliken 분할에 반영된 결과다. 반면 Cu 계는 +1 singlet이고 Cu d population이 9.83–9.84로 거의 채워져 있어 형식적 Cu(I) d¹⁰ 배정과 잘 맞는다.

4.6.2 MO별 orbital mixing과 금속–diene overlap

각 MO에서 ${\Omega }_i=w_i^L{w}_i^M$를 금속 조각–diene 조각의 mixing index로 두었다. ${\Omega }_i$가 0이면 한 fragment에 국소화된 궤도이고, 두 fragment의 weight가 같을수록 이론적 최댓값 0.25에 가까워진다. 아래 표는 $\kappa =0.05$를 통과한 결합성-overlap MO 중 $n_i{q}_i^{A-M}$ 기여가 큰 궤도들이다. 모든 궤도는 closed-shell occupied MO이므로 $n_i=2$이다.

구조	MO	${\epsilon }_i$ / Eh	${\Omega }_i$	ORCA $O_i^{A-M}$	$q_i^{A-M}=O_i/2$	$n_i{q}_i^{A-M}$
Cu reactant	25	−0.5194	0.1895	0.0608	0.0304	0.0608
Cu reactant	26	−0.4950	0.1261	0.0367	0.0184	0.0367
Cu TS	38	−0.5105	0.0595	0.0353	0.0177	0.0353
Cu TS	39	−0.4671	0.1430	0.0204	0.0102	0.0204
Cu TS	45	−0.4112	0.1148	0.0140	0.0070	0.0140
Fe reactant	23	−0.8812	0.0711	0.1122	0.0561	0.1122
Fe reactant	47	−0.2850	0.2426	0.0995	0.0498	0.0995
Fe reactant	51	−0.2175	0.2406	0.0882	0.0441	0.0882
Fe reactant	49	−0.2562	0.1463	0.0572	0.0286	0.0572
Fe TS	60	−0.3165	0.2025	0.0611	0.0306	0.0611
Fe TS	62	−0.2800	0.1738	0.0610	0.0305	0.0610
Fe TS	64	−0.2438	0.1368	0.0301	0.0151	0.0301

Cu reactant의 overlap hapticity 0.1024 중 MO 25와 26이 0.0975를 담당해 전체의 95.2%가 두 개의 지배적 궤도에 집중된다. 반면 Fe reactant에서는 표에 든 네 MO만으로 0.3571이고, 그 밖의 여러 occupied MO에도 양의 Fe–diene overlap이 분산되어 전체가 0.5106이 된다. 이는 Cu의 η² 결합이 소수의 국소적 상호작용에 집중되는 반면, Fe의 η⁴ 결합은 여러 궤도에 걸친 다중심 상호작용이라는 해석과 일치한다.

TS에서 Cu의 지배적 overlap은 유지되지만 전체 합은 70.6%만 남는다. Fe에서는 reactant의 큰 기여를 만들던 MO 23, 47, 51의 패턴이 TS의 MO 60, 62, 64 등으로 재편되고 전체 합은 41.4%만 남는다. 이 변화는 단순한 금속–탄소 거리 변화보다 직접적으로, 반응 좌표를 따라 금속–diene orbital overlap이 크게 재조직됨을 보여 준다.

4.6.3 Frontier orbital과 fragment charge 이동

Reactant complex	HOMO / Eh	LUMO / Eh	HOMO–LUMO gap / eV	겉보기 ΔG‡ / kcal·mol⁻¹
Iminium 계	−0.3605	−0.2319	3.501	6.04
AlCl₃ 계	−0.2538	−0.1205	3.625	9.15
Cu⁺–diene	−0.4201	−0.2729	4.004	21.26
Fe(CO)₃–diene	−0.2175	−0.0283	5.148	61.54

Iminium 및 AlCl₃ reactant complex는 금속–diene 계보다 작은 전역 HOMO–LUMO gap을 보이며 낮은 장벽과 같은 순서를 보였다. 그러나 이 gap은 복합체 전체의 frontier orbital 간격이므로, 실제 Diels–Alder donor–acceptor 상호작용 하나를 직접 측정한 값으로 해석해서는 안 된다. 특히 Cu와 Fe 비교에서는 금속 중심 궤도가 frontier level에 섞이므로, 위의 fragment-resolved overlap 분석을 함께 사용해야 한다.

Iminium 계에서 diene fragment의 Mulliken charge는 +0.186에서 +0.365로 증가하고 iminium fragment는 +0.814에서 +0.635로 감소한다. AlCl₃ 계에서도 diene fragment는 +0.097에서 +0.276으로 증가하고 활성화된 acrolein fragment는 +0.133에서 −0.007로 감소한다. 두 경우 모두 TS로 갈수록 diene에서 dienophile 쪽으로 약 0.18e와 0.14e의 전자밀도가 이동하는 방향이며, 낮아진 dienophile acceptor orbital을 향한 전하 이동이라는 해석과 일치한다.

4.6.4 전자구조 분석의 종합

추가 ORCA 출력은 촉매 효과를 다음 두 항의 경쟁으로 설명한다.

1. Dienophile activation: Iminium과 AlCl₃는 작은 reactant-complex frontier gap과 diene→dienophile 방향의 fragment charge 이동을 보이며, 전이상태의 전하 재분배를 수용한다.
2. Diene decoordination penalty: 금속–diene 계에서는 새 C–C 결합 형성과 동시에 기존 금속–diene orbital overlap을 재편해야 한다. $n_A^{(+)}$ 감소량은 Fe가 0.2994, Cu가 0.0301로 약 9.95배 차이 난다.

따라서 η² Cu 조건의 중간 장벽과 η⁴ Fe 조건의 높은 장벽은 “hapticity가 클수록 반응 평면을 더 잘 고정한다”는 기하학적 이점만으로는 설명되지 않는다. Fe 계는 reactant에서 훨씬 큰 다중심 overlap을 얻지만, TS에 도달하려면 그 overlap의 더 큰 비율을 잃어야 한다. 반응 평면 고정이 지나치게 강해지면 오히려 diene π 전자계의 재구성과 금속에서의 부분적 이탈 비용이 커질 수 있다는 결론이다.

Discussion

핵심 발견

현재 사용한 겉보기·내부 활성화 자유에너지 기준의 순서는 다음과 같다:

$$\text{iminium}>{\text{AlCl}}_3\text{–acrolein}>{\text{Cu}}^{+}\text{–diene}>\text{uncatalyzed}>{\text{Fe(CO)}}_3\text{–diene}$$

이 순서는 dienophile의 frontier orbital을 직접 안정화하는 전략(iminium 또는 Lewis acid)이 현재 모델의 diene 결합형 조건보다 낮은 내부 장벽을 보였음을 뜻한다. 다만 iminium·AlCl₃는 사전 반응 복합체를 기준으로 했으므로, 이를 촉매 간 절대적 성능 순위로 해석해서는 안 된다.

이론적 함의

1. FMO 이론적 해석: Iminium과 AlCl₃ 반응 복합체의 전체계 HOMO–LUMO gap은 각각 3.501, 3.625 eV로 작았고, 낮은 내부 장벽과 정성적으로 일치했다. 그러나 전체계 frontier orbital이 반응에 직접 참여하는 diene π와 dienophile π*인지 확인하지 않았으므로, gap 하나만으로 장벽 순서를 인과적으로 설명할 수 없다.
2. Hapticity 가설과 다른 결과: η⁴ Fe(CO)₃ 착물의 높은 장벽은 강한 metal–π 상호작용에 따른 diene 전자구조 재배열, 반응 평면의 고정, acrolein 접근의 입체·기하학적 제약과 양립한다. η² Cu⁺ 계는 diene을 정렬하면서도 전체 π 계를 η⁴만큼 강하게 묶지 않는 중간 조건일 수 있다. 이 해석은 가능한 메커니즘이며, 결합지수와 오비탈 등위면 분석으로 추가 검증해야 한다.
3. 온도 의존성: 열보정 차이가 온도에 따라 증가한 것은 결합 형성 과정의 엔트로피적 불리함과 양립한다. 그러나 현재 CSV에는 전자에너지가 빠져 있어 완전한 ΔG‡(T) 비교는 아직 이루어지지 않았다.

한계

• Reference state 불일치: 무촉매는 분리 반응물 합, iminium·AlCl₃는 reactant complex, 금속–diene 계는 metal–diene complex와 분리 acrolein의 합을 기준으로 했다. 현재 비교는 정성적 경향으로만 해석해야 한다.
• 계산 수준 불일치: 기본 수준은 B3LYP/def2-SVP지만, D3BJ가 촉매 계 다수에는 포함되고 acrolein·cyclopentadiene 단독 및 무촉매 TS에는 빠져 있다.
• Hapticity 변인의 교란: η²와 η⁴ 비교에서 금속 종류, 전하, 리간드 및 전자 수가 함께 변한다.
• 환경 모형 생략: 계산은 기체상이며, iminium 계의 반대이온과 실제 용매를 포함하지 않았다.
• TS 연결성: 한 개의 imaginary frequency와 mode 시각화는 확인했지만 IRC를 수행하지 않았다.
• 온도 데이터 처리: 온도 스윕 CSV의 Eel이 0이므로 기록된 dG_activation은 전체 장벽이 아니다.
• 단일 conformer 가정: 각 조건에서 가장 안정한 TS 하나만을 비교하였으며, endo/exo 선택성과 conformer ensemble을 충분히 다루지 못했다.
• 하드웨어 한계: 노트북 환경의 계산 자원 제약으로 더 큰 basis set이나 frequency convergence를 엄밀히 추구하지 못했다.

후속 연구

• 모든 구조를 동일한 B3LYP-D3BJ/def2-SVP 수준에서 다시 계산하고, 분리된 catalyst + cyclopentadiene + acrolein 기준으로 전체 ΔG‡를 구한다.
• DFT functional·basis set 의존성을 평가한다(B3LYP, M06-2X, ωB97X-D 등 교차 검증).
• 동일한 금속·리간드·전하 조건에서 η²와 η⁴ 구조를 비교하고, MO별 fragment overlap population 및 반응성 π 오비탈 등위면을 분석한다.
• Soft Lewis acid (예: Cu⁺, Ag⁺)와 hard Lewis acid (AlCl₃)의 dienophile/diene 활성화 효과를 HSAB 관점에서 체계적으로 비교한다.
• Endo/exo 선택성 및 chirality 발생 메커니즘을 iminium 촉매의 conformer 분포로 추적한다.
• TS 검증을 위해 IRC(Intrinsic Reaction Coordinate) 계산을 추가하여 양방향 mode following을 명시적으로 수행한다.
• 온도 스윕 스크립트에서 각 구조의 전자에너지를 복원해 완전한 G(T)를 재구성한다.

References

• Ahrendt, K. A., Borths, C. J., & MacMillan, D. W. C. (2000). New strategies for organic catalysis: The first highly enantioselective organocatalytic Diels–Alder reaction. Journal of the American Chemical Society, 122(17), 4243–4244. https://doi.org/10.1021/ja000092s
• Neese, F. (2022). Software update: The ORCA program system—Version 5.0. WIREs Computational Molecular Science, 12(5), e1606. https://doi.org/10.1002/wcms.1606
• Hanwell, M. D., Curtis, D. E., Lonie, D. C., Vandermeersch, T., Zurek, E., & Hutchison, G. R. (2012). Avogadro: An advanced semantic chemical editor, visualization, and analysis platform. Journal of Cheminformatics, 4, 17. https://doi.org/10.1186/1758-2946-4-17
• Fleming, I. (2010). Molecular orbitals and organic chemical reactions (Reference ed.). Wiley.
• Pearson, R. G. (1963). Hard and soft acids and bases. Journal of the American Chemical Society, 85(22), 3533–3539. https://doi.org/10.1021/ja00905a001

---