A Unified Geometric Model of Ring Strain: From Carbon to Nitrogen and the Role of Lone-Pair Coulomb Repulsion

기술 요약: 고리 변형의 통합 기하학적 모델: 탄소에서 질소까지, 그리고 비공유 전자쌍 쿨롱 반발의 역할

문제 정의
고리 변형 에너지(ring strain energy)는 고리 화합물의 안정성과 반응성을 결정하는 근본적인 요소이다. 유기 화학에서 그 중요성에도 불구하고, 서로 다른 주족 원소 전반에 걸친 변형을 설명할 수 있는 통합된 정량적 이론은 여전히 미해결 상태로 남아 있다. 기존의 접근 방식들은 한계를 지닌다: 분자 역학 힘의 장(force fields)은 양자 역학적 관측값과의 직접적인 연결이 부족하며, 초기 전자 모델(예: "바나나 결합")은 통찰력을 제공했으나 일반화 가능한 예측 기술자가 부족했고, 현대의 머신러닝 모델은 물리적 해석 가능성 없이 "블랙박스"로 기능하는 경우가 많다. 또한, 문헌에는 어떤 2주기 원소가 안정적인 단일 결합 호모원자 고리(homoatomic rings)를 형성할 수 있는지, 그리고 그 이유가 무엇인지에 대한 체계적인 탐구가 부족하다. 본 연구는 안정적인 원소들에 대한 고리 변형을 정량적으로 예측하기 위한 최소한의 물리적으로 해석 가능한 기술자 세트를 찾고, 2주기 호모원자 고리의 안정성 지형을 매핑함으로써 이러한 공백을 메우고자 한다.

방법론
본 연구는 고수준 전자 구조 이론과 고급 통계 회귀 분석을 결합한 엄격한 다단계 계산 워크플로우를 채택한다:

• 안정성 스크리닝: Be, B, C, N, O 고리($n=3–10$)의 초기 구조를 B3LYP-D3(BJ)/def2-SVP 수준에서 최적화하였다. 안정성은 해리 없이 순환적인 국소 최솟값으로 수렴하는 것으로 정의하였다.
• 컨포머 샘플링: 탄소(C3–C14) 및 질소(N3–N14) 고리에 대해 CREST/GFN2-xTB를 사용하여 컨포머 앙상블을 생성하였다. 생성된 컨포머들은 다단계 DFT 퍼널(SVP $\to$ TZVPP $\to$ QZVPP)을 통해 정제되어 앙상블 평균 깁스 자유 에너지를 얻었다.
• 고정밀 에너지 계산: 최종 깁스 자유 에너지는 DLPNO-CCSD(T)/def2-TZVPP 단일 지점 에너지와 B3LYP-D3(BJ)/def2-QZVPP로부터 유도된 열적/엔트로피적 기여를 결합한 복합 스킴을 통해 도출되었다. 고리 변형 에너지($E_{strain}$)는 이소데스믹 수소 이동 반응(isodesmic hydrogen-transfer reactions)을 통해 계산되었다.
• 기술자 추출: 최적화된 구조로부터 기하학적 기술자($\Delta\theta^2$, $\Delta R^2$)를 도출하였다. 전자적 기술자(ADCH 전하, 다중극 모멘트)는 Multiwfn을 통해 계산되었다. 질소 고리를 위해 특정된 비공유 전자쌍 쿨롱 반발($E_{coul}$) 항이 공식화되었다.
• 통계 분석: CCSD(T) 벤치마크 에너지에 대해 전체 부분 최소제곱(Full-subset OLS) 회귀를 수행하였다. 모델 선택에는 $R^2$, Leave-One-Out 교차 검증(LOOCV) RMSE, 분산 팽창 계수(VIF), 그리고 AIC/BIC 기준을 활용하였다.
• 전자적 특성 규명: N5의 이상 거동을 CASSCF(6,5) 계산과 핵 독립 화학적 이동(NICS) 계산을 통해 조사하여 $\sigma$-호모아로마티시티( $\sigma$-homoaromaticity)를 탐구하였다.

주요 기여 및 결과

1. 2주기 고리의 안정성 이분법:
   계산적 스크리닝 결과, 뚜렷한 안정성 경계가 드러났다. 탄소(C)와 질소(N)만이 넓은 크기 범위(C3–C30+ 및 N3–N14+)에 걸쳐 안정적인 단일 결합 호모원자 고리를 지지한다.

   • 베릴륨 (Be): 전자 결핍으로 인해 전자 결핍 케이지(electron-deficient cages)로 붕괴된다.
   • 붕소 (B): B4보다 큰 고리는 다중 중심 결합 불안정성으로 인해 불안정하다.
   • 산소 (O): O4 이상의 고리는 압도적인 비공유 전자쌍 반발에 의해 유도된 퍼옥시 사슬 형성을 통해 분해된다.
     이 안정성은 대칭성 일치 원리로 설명된다. C와 N은 순환적 토폴로지와 호환 가능한 전자 배치(사면체/삼각뿔형 $sp^3$)를 갖는 반면, Be, B, O는 대안적인 구조를 선호한다.

2. 탄소 고리를 위한 통합 기하학적 모델:
   탄소 고리(C3–C14)의 경우, 결합각 편차($\Delta\theta^2$)와 결합 길이 왜곡($\Delta R^2$)에 기반한 순수 기하학적 이매개변수 모델이 변형 분산의 **98.2%**를 포착하였다 ($R^2 = 0.982$, RMSE = 0.36 kcal mol$^{-1}$ atom$^{-1}$). 전자적 기술자를 포함해도 개선 효과가 미미하였으며, 이는 탄소 고리 변형이 기하학적 요인에 의해 지배됨을 확인시켜 준다.

3. 질소 고리에서의 비공유 전자쌍 반발의 역할:
   질소 고리의 경우, 기하학적 골격만으로 분산의 92.5%를 설명할 수 있다. 그러나 물리적으로 유도된 비공유 전자쌍 쿨롱 반발 항($E_{coul}$)을 포함하면 모델의 성능이 $R^2 = 0.968$ (RMSE = 0.64 kcal mol$^{-1}$ atom$^{-1}$)로 크게 향상된다. 이는 탄소에 비해 높은 수준으로 유지되는 질소 고리의 잔류 변형이 질소의 비공유 전자쌍 간의 반발 작용에 의해 구동됨을 입증한다.

4. N5의 이상 현상과 $\sigma$-호모아로마티시티:
   N5는 $-4.0$ kcal mol$^{-1}$의 잔차를 보이는 통계적 이상치로 식별되었으며, 이는 예상치 못한 안정화를 나타낸다. 이는 $\sigma$-호모아로마티시티에 기인한다: 5개의 평면 내 질소 비공유 전자쌍이 $\sigma$-골격을 통해 비편재화되어 디아트롭적 고리 전류(diatropic ring current)를 생성한다 (NICS(0) = $-15.18$ ppm). 이 효과는 비공유 전자쌍 반발을 상쇄하며, 이는 벤젠의 $\pi$-아로마티시티와는 구별되는 현상으로, 분자 평면 위에서의 고리 전류의 급격한 감쇠(NICS(1))를 통해 증명된다.

의의 및 제안된 프레임워크
본 논문은 기하학적, 정전기적, 그리고 양자 아로마틱 효과를 단일 정량적 프레임워크로 통합하는 **층위적 고리 변형 이론(layered theory of ring strain)**을 제안한다:
$$E_{strain} = E_{geom}(\Delta\theta^2, \Delta R^2) + E_{LP}(E_{coul}) + E_{\sigma-aromatic}(NICS)$$

• $E_{geom}$: 모든 주족 원소 호모원자 고리에 공통적인 보편적 기하학적 골격으로, 바이어(Baeyer)의 각 변형(angle strain)을 회복하고 결합 길이 왜곡을 추가한다.
• $E_{LP}$: 비공유 전자쌍을 가진 원소(질소계 원소)에 필요한 비공유 전자쌍 쿨롱 보정 항이며, 질소에 대해 가장 지배적인 비기하학적 요인이 된다.
• $E_{\sigma-aromatic}$: 기하학적 및 전자적 조건이 결합을 통한 비공유 전자쌍 비편재화를 가능하게 하는 특정 사례(N5와 같은 경우)에 대한 비가산적 양자 보정 항이다.

본 연구는 탄소에 대해서는 고리 변형이 근본적으로 기하학적이지만, 질소에 대해서는 비공유 전자쌍 반발에 대한 명시적 처리가 필요하며, $\sigma$-아로마티시티가 작은 고리에서 중요한 보정 역할을 한다는 것을 확립한다. 이 연구는 경험적인 힘의 장과 양자 역학적 실재 사이의 간극을 메우며, CCSD(T) 에너지와 직접 연결된 기술자를 제공함으로써 투명하고 물리적 근거가 확실한 모델을 제시한다.