시간에 따라 변하는 고전자 궤도 분해를 통한 고조파 발생 스펙트럼 해석

초록

본 연구는 실시간 TD‑CIS(또는 TD‑DFT‑TDA) 방법을 이용해 복잡한 분자계의 고조파 발생(HHG) 신호를 개별 분자오비탈(MO) 별로 정확히 분해한다. 복소 에너지를 포함한 TDSE 전파와 고에너지·연속 상태를 위한 맞춤형 가우시안 기저함수를 사용해 CO₂와 H₂O의 정렬된 구조에서 강장 레이저와의 상호작용을 시뮬레이션하였다. MO 이온화 에너지와 레이저 편광 대칭성의 결합이 HHG 스펙트럼의 특정 구간을 지배한다는 것을 확인하고, 특히 H₂O에서 HOMO, HOMO‑1, HOMO‑2의 기여도가 레이저 편광 방향에 따라 달라짐을 밝혀냈다.

상세 분석

이 논문은 고조파 발생(HHG) 현상을 분자오비탈 단위로 정량화하기 위해 실시간 시간‑의존성 전자구조 이론을 확장하였다. 핵심은 TD‑CIS(또는 KS‑기반 TD‑DFT‑TDA) 파동함수를 사용해 전자 파동함수를 기저 상태와 단일 전자 여기 상태들의 선형 결합으로 표현하고, 복소 에너지를 도입해 이온화 손실을 자연스럽게 포함시킨 점이다. 이렇게 얻어진 계의 시간‑의존성 계수 cₙ(t)를 기반으로 전이 쌍극자 모멘트 μ(t)를 세 부분(ground‑ground, ground‑excited, excited‑excited)으로 정확히 분해하고, 각 오비탈 i에 대한 μ_i(t)를 도출한다. 이후 푸리에 변환을 통해 HHG 스펙트럼 P(ω)를 각 오비탈 기여 P_i(ω)와 오비탈 간 간섭 항으로 나누어 분석한다.

이론적 구현에서는 H₀를 대각화 가능한 형태로 두고, 레이저와의 상호작용 V(t)를 비대각 행렬로 처리한다. 분할 전파법(split‑operator)으로 시간 전파를 수행하며, V(t)의 지수 연산은 고정된 유니터리 변환 U를 이용해 대각화 후 계산한다. 전자 기저는 5‑fold augmented Dunning 계열에 연속 상태를 모사하는 최적화된 Gaussian K‑함수를 추가한 5aug‑cc‑pVDZ+5K(CO₂)와 5aug‑cc‑pVTZ+5K(H₂O)를 사용해 고에너지 및 연속 영역을 충분히 커버한다. 800개의 전자 여기 상태를 포함해 파동함수를 전개함으로써 실질적인 이온화 채널을 모두 고려한다.

실험적 검증을 위해 CO₂와 H₂O를 각각 800 nm, 8.5×10¹³ W cm⁻², 23 optical‑cycle 펄스로 조사하였다. CO₂에서는 HOMO와 HOMO‑1이 레이저 편광이 분자축에 평행할 때 강하게 기여하고, HOMO‑2는 거의 무시된다. 이는 기존 실험에서 보고된 최소(minimum) 현상이 다중‑MO 간 간섭에 기인함을 재현한다. H₂O에서는 레이저 편광이 분자 평면에 수직일 때 HOMO가 주된 기여를, 평면에 평행할 때는 HOMO‑1과 HOMO‑2가 각각 고조파 영역(특히 15‑30번째 고조파)에서 두드러진 기여를 보인다. 이러한 편광 의존성은 각 MO의 이온화 전위와 전이 쌍극자 행렬 요소가 레이저 전기장 대칭성과 어떻게 매칭되는가에 따라 달라진다.

결과적으로, 이 방법은 (1) MO 별 기여와 (2) 기여 간 위상·진폭 간섭을 정량적으로 제공함으로써 HHG 스펙트럼을 “분자 오비탈 망원경”으로 활용할 수 있음을 입증한다. 또한, ground‑excited와 excited‑excited 전이 항을 포함함으로써 기존 근사(ground‑excited만 고려)보다 더 정확한 스펙트럼 재현이 가능함을 보여준다.