비직교 국소 궤도 기반 2차 광학 비선형성 계산법

초록

본 논문은 비직교 의사원자 궤도(PAO)를 이용해 속도 게이지 섭동 이론으로 2차 광학 비선형성 χ^(2)(−2ω;ω,ω)를 계산하는 새로운 프레임워크를 제시한다. 슬레이터‑코스터 형태의 두 중심 모멘텀 적분을 대칭으로 식별하고, 이를 SIESTA 기반 비직교 기저에 적용한다. 3C‑SiC와 GaAs에 대한 테스트 결과는 기존 평면파 및 Wannier 기반 방법과 일치함을 보여준다.

상세 분석

이 연구는 비직교 의사원자 궤도(PAO)를 기반으로 한 속도 게이지 섭동 이론을 이용해 반도체의 두 번째 차수 광학 감수도 χ^(2)_{ijk}(−2ω;ω,ω)를 계산하는 체계적인 방법을 개발하였다. 기존의 대부분의 비선형 광학 계산은 평면파 기반의 직교 기저 혹은 Wannier 함수로 후처리하는 방식에 의존했으며, 이는 진공이 큰 시스템(분자, 나노튜브, 2D 물질 등)에서 계산 비용이 급증하거나 Wannier화 과정에서 손실이 발생하는 단점을 가지고 있었다. 저자들은 SIESTA에서 제공하는 비직교 PAO를 그대로 활용함으로써 이러한 문제를 회피하고, 특히 비직교 기저의 중첩 행렬 S(k)와 해밀토니안 H(k)를 일반화된 고유값 문제 H(k)C_n(k)=ε_n(k)S(k)C_n(k) 형태로 풀어 Bloch 상태를 구성한다.

핵심 기술은 두 중심 모멘텀 연산자 ⟨φ_α|p̂_i|φ_β⟩를 슬레이터‑코스터(Slater‑Koster) 방식으로 전개하고, 결정 대칭성을 이용해 필요한 적분을 최소화하는 것이다. 구체적으로, 각 PAO는 구형조화함수와 방사형 함수를 곱한 형태이며, 구형조화함수는 Condon‑Shortley 위상을 포함한 실구형조화함수로 정의된다. 이러한 정의는 각 궤도에 대한 회전 행렬을 명시적으로 제공함으로써, 서로 다른 원자 위치에 있는 PAO 간의 두 중심 적분을 좌표 변환을 통해 효율적으로 계산한다.

속도 게이지 섭동 이론에서 χ^(2)는 두 개의 항 A_{ijk}(2ω)와 B_{ijk}(ω)로 분해된다. A 항은 두 광자 흡수·방출 과정을, B 항은 하나의 광자와 전자·정공 사이의 비탄성 전이를 기술한다. 저자들은 이 두 항을 각각 Im