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‡ 재계산.

Abstract

본 연구는 cyclopentadiene과 acrolein 사이의 Diels–Alder 반응에 대해 다양한 촉매 조건에서의 활성화 자유에너지(ΔG‡)를 ORCA와 Avogadro를 이용한 양자화학 계산으로 비교하였다. 무촉매 반응, iminium 유기촉매, AlCl₃ Lewis acid 촉매, 그리고 diene에 결합한 두 종류의 금속 착물(Fe(CO)₃, Cu⁺) 조건을 설정하고 각 조건에서 전이상태를 OptTS 계산으로 탐색하였다. 무촉매 반응은 31.86 kcal/mol의 ΔG‡를 나타냈으며, iminium 촉매(6.04 kcal/mol)와 AlCl₃–acrolein 착물(9.15 kcal/mol)이 가장 큰 활성화 에너지 감소 효과를 보였다. 반면 Fe(CO)₃–diene 착물은 오히려 ΔG‡를 증가시켰고(61.54 kcal/mol), Cu⁺–diene 착물은 무촉매보다 유리했으나(21.26 kcal/mol) 유기촉매·Lewis acid 방식보다 약했다. 결과는 dienophile의 LUMO를 직접 안정화하는 촉매 전략이 diene–metal 착물 형성을 통한 접근 방식보다 효과적임을 시사한다.

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 착물(예: cyclopentadiene-Fe(CO)₃)은 cyclopentadiene의 π 전자계가 금속과 hapticity(η²–η⁴)로 결합하여 반응 평면이 고정되고, orbital overlap 분포가 변화한다 @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‡가 낮게 유지되지만, 고온에서는 촉매 기능이 약화되어 무촉매 반응의 ΔG‡로 점근할 것이다.
• H2: Orbital overlap의 재배열과 reaction plane 고정 효과를 고려할 때, 금속-diene 착물의 hapticity가 클수록 ΔG‡가 더 크게 감소할 것이다.

Method

Materials (Computational Setup)

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

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

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의 유기촉매 전략을 따라 acrolein이 iminium ion 형태로 활성화되는 모델을 구성하였다. 섭동 이론(perturbation theory)에 의하면 iminium 형성은 dienophile의 LUMO를 낮추어 diene HOMO–dienophile LUMO 간격을 감소시킨다. 활성화 자유에너지는 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을 이용하여 환산하였다. 온도 의존성은 270–1070 K 범위에서 표준 압력(1 atm) 조건으로 ORCA의 thermal correction을 통해 산출하였다.

변인 통제 주의

본 연구에서 무촉매 조건은 분리 반응물 합 기준, iminium·Lewis acid 조건은 reactant complex 기준으로 계산되었다. 따라서 절대값 비교에는 reference state mismatch 문제가 존재하며, 후속 연구에서는 모든 조건을 동일한 reference state(예: 모든 분리 반응물 합)로 재계산할 필요가 있다.

Results

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

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

구조	Gibbs 자유에너지 / Eh	HOMO / Eh	LUMO / Eh	비고
Acrolein	−191.62627218	−0.25848	−0.06911	dienophile
Cyclopentadiene	−193.76577644	−0.21783	−0.01650	diene
Iminium (acrolein-derived)	−250.57600236	−0.45790	−0.25873	acrolein 활성화 모델
AlCl₃–acrolein complex	−2008.02977059	−0.25913	−0.11486	Lewis acid–dienophile
Fe(CO)₃–Cp complex	−1796.91550073	−0.21746	−0.02829	η⁴ diene complex
Cu⁺–Cp complex	−1833.77229864	−0.21746	−0.02829	η² diene complex

Iminium의 LUMO는 −0.25873 Eh로 acrolein 단독(−0.06911 Eh) 대비 약 0.19 Eh 낮아졌으며, AlCl₃ 착물 또한 acrolein의 LUMO를 −0.11486 Eh로 안정화시켰다. 반면 Fe(CO)₃·Cu⁺와 cyclopentadiene의 착물에서 LUMO 변화는 미미하였다.

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

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

조건	G_TS / Eh	기준 반응물 G / Eh	ΔG‡ / Eh	ΔG‡ / kcal·mol⁻¹	해석
무촉매 DA	−385.34127499	−385.39204862	0.05077363	31.86	기준 반응
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로 환산하였다.

4.3 온도 의존성 (RQ1)

무촉매 반응

270–1070 K 범위에서 ΔG_TS는 약 0.020 Eh에서 0.079 Eh로 거의 선형 증가하였다 (Fig. 1). 이는 H1의 예측과 일관되며, TS 형성이 엔트로피적으로 불리하다는 사실을 반영한다.

Iminium 촉매 조건

같은 온도 범위에서 ΔG_TS가 0.002 Eh에서 0.011 Eh로 증가하였으며, 무촉매 대비 절대값 자체가 한 자리수 작았다 (Fig. 2). 다만 이 비교는 reactant complex를 기준으로 계산되었으므로 변인 통제가 완전하지 않다.

4.4 Hapticity 효과 (RQ2)

가설 H2(hapticity ↑ → ΔG‡ ↓)는 본 계산에서 기각되었다. η⁴ 결합(Fe(CO)₃)이 오히려 ΔG‡를 61.54 kcal/mol까지 증가시켰고, η² 결합(Cu⁺)은 무촉매보다 약간만 유리한 21.26 kcal/mol를 보였다.

Discussion

핵심 발견

활성화 자유에너지 기준 촉매 효과는 다음 순서로 나타났다:

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

즉, dienophile의 LUMO를 직접 낮추는 전략(iminium 또는 Lewis acid)이 diene 측 활성화보다 우수하다는 결론을 얻었다.

이론적 함의

1. FMO 이론적 해석: dienophile LUMO 안정화는 diene HOMO–dienophile LUMO gap을 직접적으로 줄이므로, perturbation 이론의 분자 적분 ⟨ψ_HOMO|V̂|ψ_LUMO⟩가 가지는 에너지 분모의 영향을 직접 제어한다. 반면 diene 측에 금속이 결합하면 HOMO 분포는 재배열되지만 dienophile LUMO는 그대로이므로 gap 감소 효과가 제한적이다.
2. Hapticity 가설의 반증: η⁴ Fe(CO)₃ 착물에서 ΔG‡가 오히려 증가한 것은, 강한 metal–π 결합이 diene HOMO를 안정화시켜 dienophile LUMO와의 gap을 벌려놓는 부작용을 시사한다. 또한 acrolein의 접근 시 입체 장애가 크게 작용했을 가능성도 있다.
3. 온도 의존성: 두 조건 모두에서 ΔG_TS가 온도에 따라 선형 증가하는 양상은 TS 형성의 엔트로피적 불리함을 확인시켜준다. 다만 본 계산에서는 촉매 분해·기능 상실의 임계 온도 효과는 명시적으로 모델링되지 않았다.

한계

• 변인 통제 문제: 무촉매 조건과 촉매 조건의 reference state가 일치하지 않음 (분리 반응물 합 vs. reactant complex). 본 연구의 비교 결과는 정성적 경향성으로만 해석되어야 한다.
• 계산 수준: 사용한 functional, basis set, dispersion correction의 명시가 부족하여 절대 ΔG‡ 값의 신뢰도에 한계가 있다. 특히 전이금속이 포함된 시스템에서는 multi-reference 효과 검토가 필요하다.
• 단일 conformer 가정: 각 조건에서 가장 안정한 TS 하나만을 비교하였으며, endo/exo 선택성과 conformer ensemble을 충분히 다루지 못했다.
• 하드웨어 한계: 노트북 환경의 계산 자원 제약으로 더 큰 basis set이나 frequency convergence를 엄밀히 추구하지 못했다.

후속 연구

• 모든 조건을 동일한 reference state로 통일하여 ΔG‡ 절대값 비교를 다시 수행한다.
• DFT functional·basis set 의존성 평가 (B3LYP, M06-2X, ωB97X-D 등 교차 검증).
• Soft Lewis acid (예: Cu⁺, Ag⁺)와 hard Lewis acid (AlCl₃)의 dienophile/diene 활성화 효과를 HSAB 관점에서 체계적으로 비교한다.
• Endo/exo 선택성 및 chirality 발생 메커니즘을 iminium 촉매의 conformer 분포로 추적한다.
• TS 검증을 위해 IRC(Intrinsic Reaction Coordinate) 계산을 추가하여 양방향 mode following을 명시적으로 수행한다.

References

APA 7판 형식, 알파벳순. 모든 인용은 본문에서 다뤄져야 함.

• 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

---

상위 노트: Study MOC / Chemistry MOC / Computational Chemistry MOC