Prediction of experimental excited-state absorption spectra by using vibronic transition calculations: Its practical application to a {\pi}-conjugated molecule

문제 정의
들뜬 상태 흡수(ESA) 스펙트럼의 정확한 이론적 예측은 계산 화학 분야에서 여전히 중요한 과제로 남아 있다. 과도 흡수(TA) 및 광유기 흡수(PIA) 스펙트럼은 빛에 의해 구동되는 분자 과정(예: 광촉매 작용, 에너지 전달)을 연구하는 데 필수적이지만, 기태 상태 블리칭(GSB), 유도 방출(SE), 그리고 ESA 기여도가 중첩되어 분석이 복받다. 밀도 범함수 이론(DFT)과 시간 의존 밀도 범함수 이론(TDDFT)은 기태 상태 흡수 및 방출을 신뢰성 있게 예측하지만, ESA 스펙트럼을 계산하는 것은 어렵다. 이는 표준 선형 응답 TDDFT가 두 개의 전자 들뜬 상태 사이의 전이 특성을 제공할 수 없기 때문이다. 이를 해결하기 위해 제2 선형 응답(SLR), 이차 응답(QR), 실시간(RT) TDDFT와 같은 고급 방법들이 개발되었으나, 이들은 시간 의존 교환-상관 포텐셜로 인한 정확도 저하, 예측 불가능한 에너지 이동, 또는 최종 전자 상태를 정확하게 할당하지 못하는 문제 등의 한계를 겪는 경우가 많다. 결과적으로, 실험적인 TA 및 PIA 스펙트럼을 해석하기 위한 이러한 특화된 방법들의 실제적인 적용은 제한적이다.

방법론
저자들은 표준 양자 화학 패키지(구체적으로 Gaussian 16)에서 사용 가능한 일상적인 DFT 및 TDDFT 계산과 진동 준위 전이 계산을 결합하여 분자 ESA 스펙트럼을 계산하는 실용적인 접근법을 제안한다. 이 방법은 다음의 프레임워크에 기반한다:

1. 포텐셜 에너지 표면(PESs): DFT 및 TDDFT(구체적으로 APFD 범함수와 6-311++G(d,p) 기저 집합)를 사용하여 바닥 상태($S_0$)와 관련 들뜬 상태($S_n$)에 대한 최적화된 기하 구조, 전자 에너지, 그리고 진동 법선 모드를 계산한다.
2. 진동 준위 전이 계산: 298 K에서 단열 헤시안(adiabatic Hessian) 프레임워크 내의 프랭크-콘돈-허츠베르크(Franck–Condon–Herzberg–Teller) 근사를 사용하여 전자 상태 간의 진동 준위 전이를 계산한다.
3. 스펙트럼 구축: 계산된 진동 준위 선들을 가우시안 엔벨로프(Gaussian envelopes)와 컨볼루션하여 밴드 형태를 시뮬레이션한다. 반치폭(FWHM) 및 보정 계수는 계산된 피크를 실험 데이터와 비교하여 결정된다.
4. 강도 근사: 이 방법의 핵심적인 단순화는 선형 응답 TDDFT에서 구할 수 없는 두 들뜬 상태 사이의 특정 전이 모멘트($S_n \to S_{n'}$)를 직접 계산하지 않는다는 점이다. 대신, $S_n \to S_{n'}$ 전이에 대한 전자 전이 모멘트를 대응하는 $S_0 \to S_{n'}$ 전이의 모멘트를 사용하여 근사한다. 결과적으로, 이 방법은 절대적 또는 상대적 강도보다는 진동 밴드의 형태와 피크 위치를 재현하는 데 집중한다.

주요 기여 및 결과
본 연구는 비정상적인 $S_2$ 상태로부터의 형광과 복잡한 들뜬 상태 역학으로 알려진 분자인 아줄렌(azulene)을 모델 시스템으로 사용하여 이 접근법을 검증한다.

• 기태 및 방출 스펙트럼: 이 방법은 아줄렌의 실험적 흡수 및 방출 스펙트럼을 성공적으로 재현하였으며, 이를 통해 계산된 전이 에너지, PES, 그리고 진동 법선 모드가 상당히 정확함을 확인하였다.
• ESA 스펙트럼 할당: 저자들은 $S_1$ 및 $S_2$ 상태로부터 높은 들뜬 상태로의 전이에 대한 ESA 스펙트럼을 계산하였다.
  • $S_1$ 상태의 경우, ~450 nm 및 ~375 nm에서의 계산된 진동 특징을 각각 $S_1 \to S_3$ 및 $S_1 \to S_4$ 전이로 할당하였다. 이는 해당 피크를 $S_1 \to S_6$ 및 $S_1 \to S_7$로 제안했던 기존의 CASSCF 계산 기반 할당과 대조된다.
  • $S_2$ 상태의 경우, ~575 nm에서의 피크를 $S_2 \to S_8$ 및 $S_2 \to S_9$ 전이로 할당하였으며, 이는 기존에 할당되었던 $S_2 \to S_{11}$과 차이가 있다.
• 진동 모드 식별: 계산은 전자 전이에 결합된 특정 진동 법선 모드를 성공적으로 식별하였다. 예를 들어, 면외 굽힘 모드($\nu_2$)가 $S_0 \to S_4$ 및 $S_1 \to S_4$ 전이 동안 발생하는 구조적 변화(평면에서 굽은 구조로)의 주요 동인임을 확인하였다.

의의 및 주장
본 논문은 이 접근법이 특수한 응답 이론 구현(예: QR-TDDFT)이나 명시적인 들뜬 상태 간 전이 모멘트 계산 없이도 낮은 에너지의 들뜬 상태 사이의 전이로부터 발생하는 ESA 밴드를 할당할 수 있는 "실용적인 프레임워크"를 제공한다고 주장한다.

• 접근성: 이 방법은 표준 DFT 및 TDDFT 계산을 활용하므로, 고급 응답 이론 코드를 사용할 수 없는 연구자들도 접근 가능하다.
• 형태 및 위치에 초점: 저자들은 이 방법이 ESA 밴드를 직접 할당하는 데 필요한 진동 밴드 형태와 피크 위치를 제공한다고 명시적으로 밝힌다. 또한 전이 모멘트의 근사로 인해 정확한 절대적 또는 상대적 ESA 강도를 제공하지 못한다는 점을 인정한다.
• 한계: 저자들은 조화 근사(highly anharmonic PES의 경우 실패할 수 있음)에 의존하며, 전이 에너지의 정확도가 선택된 DFT 범함수에 달려 있다는 점을 언급하며 방법론의 범위에 대해 겸허한 태도를 유지한다. 그들은 이 방법이 낮은 에너지의 들뜬 상태를 포함하는 TA 및 PIA 스펙트럼을 해석하는 데 가장 적합하다고 강조한다.

결론적으로, 본 연구는 진동 준위 전이 계산이 표준 TDDFT와 결합될 때 실험적 ESA 스펙트럼의 진동 특징과 밴드 위치를 효과적으로 재현할 수 있음을 보여주며, 분자 시스템의 과도 흡수 밴드를 할당하는 데 있어 실행 가능한 경로를 제시한다.