Revisiting Theoretical Modeling of Seebeck Coefficient of Semiconductors

기술 요약: 반도체의 제베크 계수(Seebeck Coefficient)에 대한 이론적 모델링의 재고

문제 제기
본 논문은 반도체의 제베크 계수($\alpha_S$)를 모델링하는 데 사용되는 기존 방법론, 특히 볼츠만 수송 방정식(BTE)과 비열 유속(specific heat flux) 정의에 의존하는 방식에서 발견되는 근본적인 열역학적 결함을 식별한다. 저자는 기존 접근 방식의 세 가지 주요 문제를 다음과 같이 강조한다:

1. 잠재적 혼동: 전위(electrochemical potential, $E_F$)와 정전 전위(electrostatic potential, $\varphi$)의 정의를 빈번하게 혼용하여 잘못된 유도를 초래함.
2. 임의적 무시: BTE 유도 과정에서, 개방 회로(open-circuit, OC) 조건에서도 샘플 전체에 걸쳐 변화하는 변수인 전도대 하단($E_C$)과 페르미 준위($E_F$)의 공간적 구배(gradient)를 분포 함수의 구배를 계산할 때 임의로 무시함.
3. 부적절한 열 유속: 국소 열 유속($J_Q = T J_S = J_E - \mu J_N$)의 정의는 엔트로피 생성(entropy production)이 없는 국소 준정적(가역적) 과정이라는 강력한 가정에 의존함. 이는 근본적인 비평형 열역학 원리에 위배되나, BTE, 랜드아우어-뷔티커(Landauer-Büttiker), 쿠보(Kubo) 공식 적용 시 널리 사용됨.

방법론
이러한 결함을 해결하기 위해, 저자는 개방 회로 조건에서의 **소레 효과(Soret effect)가 포함된 드리프트-확산 방정식(drift-diffusion equation)**에 기반한 대안적 방법론을 제 제안한다. 이 접근 방식은 명시적인 열 유속 제시와 전기화 // 정전 전위 사이의 모호성을 피한다.

유도는 세 단계로 진행된다:

1. 근평형 상태 (밴드 항): 국소 평형(Local Equilibrium, LEA)과 국소 전기적 중성을 가정하여, $\alpha_S$의 밴드 항에 대한 폐쇄형 모델을 유도한다. 순 전하 유속을 zero로 설정하고 일반화된 아인슈타인 관계식을 활용함으로써, "가역 모델(reversible model)"을 회복한다. 이 모델은 전자 속도 및 완화 시간 항을 명시적으로 제거하고, 대신 페르미 준위의 공간적 구배에 의존한다.
2. 포논 드리프트(Phonon-drag) 항: 비평형 포논과의 운동량 교환에서 발생하는 전하 유속($J_d$) 항을 추가하여 포논 드리프트 효과를 포함하도록 방법론을 확장한다. 그 결과 도출된 포논 드리프트 항의 식은 기존 BTE 기반 모델과 동일한 형태임을 보여준다.
3. 탄성(Ballistic) 사례: 탄성 영역(LEA가 유효하지 않은 경우)에서는 랜드아우어 공식을 사용하여 핫 터미널(hot terminal)과 콜드 터미널(cold terminal)로부터의 전하 유속을 계산하여, 기존의 탄성 제베크 계수 식을 유도한다.

이 논문은 이 모델들(BTE, 가역 모델, 탄성 모델)을 에너지 지수 $r$에 의존하는 단일 적분 형태로 통합하며, $r$은 산란 메커니즘과 물질의 차원에 따라 달라진다.

주요 기여

• 열역학적 교정: 본 논문은 표준 BTE 유도 및 열 유속 정의에 존재하는 열역학적 불일치를 엄격히 비판한다.
• 가역 모델의 유도: 저자는 이전 연구에서 제시된 가역 모델이 속도 및 완화 항을 생략했음에도 불구하고, 근평형 사례에서 제베크 계수의 밴드 항으로서 타당함을 입증한다.
• 이론적 경계 설정: 본 연구는 가역 모델이 제베크 계수 밴드 항의 이론적 **하한선(lower bound)**을 제공함을 확립한다(3D 물질의 경우 $r=1/2$에 해당). 반면, 약하게 스크리닝된 이온화 불순물 산란을 사용하는 BTE 모델은 상한선($r=3$)을 제공한다.
• 통합 프레임워크: 상태 밀도(density of states), 속도 및 완화 시간의 에너지 의존성에 기초하여 가역, 탄성, BTE 모델을 비교할 수 있는 통합된 수학적 프레임워크를 제시한다.

결과
고농도 n형 실리콘($n_d > 1 \times 10^{19} \text{ cm}^{-3}$)에 대해 사례 연구를 수행하여, 세 가지 모델을 Geballe 등과 Sadhu 등의 실험 데이터와 비교하였다:

• 가역 모델: 낮은 계산 비용으로 "상당히 양호한" 예측을 제공했다. 이 모델은 실험 데이터의 하한선 역할을 지속적으로 수행하였으며, 도핑 농도가 $2.7 \times 10^{19}$, $6.0 \times 10^{19}$, $7.0 \times 10^{19} \text{ cm}^{-3}$일 때 각각 13%, 18%, 5%의 평균 상대 편차를 보였다.
• 탄성 모델: 혼합된 성능을 보였으며, 편차는 3%, 6%, 16%였다. 나노와이어 케이스($7.0 \times 10^{19} \text{ cm}^{-3}$)에서 데이터를 과대평가하였다.
• BTE 모델: 약하게 스크리닝된 이온화 불순물 산란을 사용할 경우, 모델은 실험값을 크게 과대평가하였으며 평균 상대 편차는 60%에서 90%에 달했다. 이는 BTE 결과가 완화 시간 공식의 선택에 매우 민감함을 강조한다.

의의 및 주장
본 논문은 제안된 가역 모델이 열전 특성, 특히 근평형 및 가역 사례를 평가하는 데 있어 간결하고, 이론적으로 견고하며, 계산 효율적인 대안을 제공한다고 주장한다. 복잡한 산란 시간 계산을 제거하고 열역학적 역설을 피함으로써, 이 모델은 고처리량 계산 및 재료 스크리닝을 위한 유효한 도구로 제시된다. 저자는 이 모델이 비가역적 메커니즘(산란)을 무시함에도 불구하고 제베크 계수의 하한선을 정확하게 포착하므로, 열전 특성 평가를 위한 강력한 기준점(baseline)이 된다고 강조한다.