Impact of dynamic electrostatic disorder on hole mobility in rubrene: a nonadiabatic molecular dynamics investigation

기술 요약: 루브렌(Rubrene) 내 정공 이동도에 미치는 동적 정전기적 무질서의 영향

문제 제기
루브렌과 같은 고이동도 유기 분자 결정은 유기 전자 공학에 매우 중요하지만, 전하 수송을 정량적으로 예측 가능한 방식으로 기술하는 것은 여전히 어려운 과제로 남아 있다. 이러한 물질에서의 전하 수송은 소폴라론 호핑(small-polaron hopping)과 밴드 수송(band transport) 사이의 중간 영역에서 발생하며, 이는 전하 운반체의 파동 함수가 열적 요동에 의해 변조되는 과도적 비국소화(transient delocalization) 특성을 갖는다. 혼합 양자-고전 비단열 분자 역학(특히 Fragment Orbital 기반 Surface Hopping, 즉 FOB-SH)은 특정 수송 메커니즘을 가정하지 않고 이 영역을 시뮬레이션할 수 있는 유망한 경로를 제공하지만, 루브렌과 같은 무극성 결정에 대한 기존의 대규모 시뮬레이션들은 일반적으로 동적 정전기적 무질서를 간과해 왔다. 이러한 누락은 장거리 쿨롱 상호작용(일반적으로 Ewald 합산법 필요)을 평가하는 데 드는 계산 비용 때문이었으며, 이러한 상호작용이 무극성 시스템에서는 미미할 것이라는 가정에 기인했다. 그러나 이러한 상호작력을 무시하면 재배열 에너지(reorganization energy)에 대한 정전기적 기여가 제거되고 사이트 에너지 요동이 억제되어, 운반체의 비국소화 및 이동도 예측에 있어 부정확한 결과를 초래할 수 있다.

방법론
본 연구는 루브렌의 실온 정공 수송에 대한 FOB-SH 시뮬레이션에 동적 정전기적 무질서를 통합한다. 저자들은 여러 전하 국부 상태(charge-localized states)에 대한 정전기적 기여를 평가할 때 발생하는 계산 병목 현상을 극복하기 위해, 효율적인 가감법(addition-subtraction scheme)과 결합된 Damped Shifted-Force (DSF) 실공간 정전기 합산법을 사용한다.

주요 방법론적 구성 요소는 다음과 같다:

• 정전기 처리: DSF 방법은 유한한 컷오프($R_{cut} = 16$ Å)에서 0으로 이동된 감쇠된 실공간 상호작용을 사용하여 장거리 쿨롱 상호작용을 근사한다. 가감법은 전하 국부 상태의 정전기 에너지를 중성 시스템의 에너지와 보정 항의 합으로 분해하여, 계산 스케일링을 단일 계산 대비 $O(N_{mol})$에서 거의 $O(1)$ 수준으로 줄여준다.
• 편극 효과: 명시적인 편극력장(polarizable force fields)의 비용 없이 전자 편극을 반영하기 위해, 원자 부분 전하에 전하 스케일링(계수 $\gamma = 0.80$)을 적용한 암시적 기법을 적용하였다.
• 시뮬레이션 프로토콜: FOB-SH 시뮬레이션은 300 K에서 $50 \times 13 \times 1$ 단위 격자(1300개 분자)의 슈퍼셀을 사용하여 수행되었다. 시뮬레이션은 핵을 고전적으로 진화시키는 동안 정공 파동 함수를 양자 역학적으로 전파하였다. 결과는 "제로 전하(zero-charge)" 모델(부분 전하를 0으로 설정한 모델)과 비교되었으며, 재배열 에너지에 대해 Smooth Particle-Mesh Ewald(SPME) 계산을 통해 검증되었다.
• 분석: 본 연구는 재배열 에너지($\lambda$), 사이트 에너지 무질서($\sigma_E$), 전자 결합 요동, 파동 함수 비국소화 정도를 측정하기 위한 역참여율(Inverse Participation Ratio, IPR), 그리고 평균 제곱 변위(MSD)로부터 유도된 전하 운반체 이동도($\mu$)를 정량화하였다.

주요 결과

1. 재배열 에너지: 정전기적 상호작용(DSF 또는 SPME를 통해)을 포함하면 제로 전하 모델에 비해 정공 전달을 위한 재배열 에너지가 약 20–40 meV 증가한다(기초값 152 meV에서 ~175–191 meV로). 스펙트럼 분석 결과, 이러한 증가는 저주파 격자 운동(<200 cm$^{-1}$)과 고주파 분자 내 모드 모두에 걸친 결합 강화에서 기인함을 보여준다. DSF 방법은 이러한 정전기적 기여의 대부분을 포착하며, SPME 벤치마크와 잘 일치하는 모습을 보인다.
2. 사이트 에너지 무질서: 정전기 포함은 사이트 에너지의 제곱평균제곱근(RMS) 요동($\sigma_E$)을 제로 전하 모델의 63.5 meV에서 DSF의 69.2 meV로 유의미하게 증가시킨다. 이러한 증가는 증가된 재배열 에너지와 정량적으로 일치하며, 이는 정전기가 오프-다이아고날(off-diagonal) 전자 결합에는 거의 영향을 주지 않으면서 주로 다이아고날(diagonal) 에너지 무질서를 수정한다는 것을 확인시켜 준다.
3. 파동 함수 비국소화: 강화된 사이트 에너지 무질서는 정공 파동 함수의 공간적 범위를 감소시킨다. 평균 역참여율(IPR)은 제로 전하 모델의 13에서 DSF의 9로 감소하며, 이는 전하 운반체가 더 국부화됨을 나타낸다.
4. 전하 이동도: 비국소화의 감소는 예측된 정공 이동도의 상당한 감소로 이어진다. 고이동도 $a$-방향을 따라, 이동도는 제로 전하 모델의 34.8 cm$^2$V$^{-1}$s$^{-1}$에서 DSF의 21.1 cm$^2$V$^{-1}$s$^{-1}$로 떨어진다. 이 값은 실험적 참조값인 20.2 cm$^2$V$^{-1}$s$^{-1}$와 밀접하게 일치한다. 또한, 정전기 포함은 예측된 이방성 비율($\mu_a/\mu_b$)을 개선하여 실험 관측치에 더 가깝게 만든다.

의의 및 주장
본 논문은 동적 정전기적 무질서가 루브렌과 같이 무극성 분자로 구성된 고이동도 유기 반도체에서도 정량적으로 예측 가능한 전하 수송 시뮬레이션을 위한 필수적인 요소임을 입증한다. 저자들은 효율적인 DSF 체계를 통해 이러한 상호작용을 명시적으로 포함함으로써, 정전력을 무시한 기존 모델들이 운반체의 비국소화와 이동도를 과대평가했음을 보여준다. 본 연구는 정전기적 사이트 에너지 요동이 고이동도 유기 결정에서 관찰되는 동적 국부화의 주요 메커니즘임을 확립하며, 이론적 시뮬레이션과 실험적 현실 사이의 간극을 메우기 위해 정전력의 포함이 필수적임을 밝혀냈다. 또한, 본 연구는 주기적 시스템에서의 대규모 비단열 역학을 위한 확장 가능한 대안으로서 DSF 방법의 타당성을 검증하였다.