[논문 리뷰] Simple remedy for spurious states in discrete k(dot)p models of semiconductor structures

William R. Frensley, Raja N. Mir|arXiv (Cornell University)|2014. 12. 22.

Spectral Theory in Mathematical Physics인용 수 4

한 줄 요약

이 논문은 이산 k·p 모델에서 표준 중심차분 이산화를 일阶 상류 유한차분 형식으로 대체하여 임의의 상태(이차형식에서 비단조적 분산 관계로 인한 산물)를 제거한다. 이 방법은 물리적 대칭성을 유지하고 연속체 밴드 구조에 더 가까운 정확도를 보이며 절대적 안정성을 확보하며 고파동수 진동 해를 손상 없이 제거한다.

ABSTRACT

The spurious states found in numerical implementations of envelope function models for semiconductor heterostructures and nanostructures are artifacts of the use of the centered-difference formula. They are readily removed by employing a first-order difference scheme. The technique produces none of the loss of accuracy that is commonly feared from lower-order formulations; the fidelity to the continuum band structure is actually improved. Moreover the stability is absolute, thus no adjustments to the model parameters are required. These properties are demonstrated in an explicit calculation of a 3-band model of a semiconductor superlattice.

연구 동기 및 목표

• 반도체 이종구조 및 나노구조의 수치적 k·p 모델에서 지속적인 문제로 남아 있는 임의의 상태를 해결한다.
• k·p 이산화에서 정확도와 에르미트성 확보를 위해 중심차분이 필수적이라는 기존의 믿음을 도전한다.
• 일阶 비대칭 차분이 물리적 대칭성과 안정성을 유지하면서도 임의의 해를 제거할 수 있음을 보여준다.
• 일阶 차분 형식이 비단조적 분산 관계를 가지며 이차형식 중심차분의 역행적 행동을 피하는 단조적 밴드 분산을 생성함을 보여준다.
• 세부 조정이 필요하거나 복잡한 수정이 요구되는 기존 방법의 대안으로서, 매개변수 조정이 불필요한 강력한 대체 방법을 제공한다.

제안 방법

• 표준 중심차분 기울기 연산자 $\hat{K}_C$ 를 $\hat{K}_L$ 과 $\hat{K}_R$ 로 구성된 일阶 상류 연산자로 대체한다. 여기서 $\hat{K}_L f_n = -i/\Delta (f_n - f_{n-1})$ 이고 $\hat{K}_R f_n = -i/\Delta (f_{n+1} - f_n)$ 이다.
• 항등식 $\hat{K}_L \hat{K}_R = \hat{K}_R \hat{K}_L = \hat{L}$ 을 사용하며, 여기서 $\hat{L}$ 은 이산 라플라시안으로서 이차형식 운동에너지 항의 일致성을 보장한다.
• 표준 $\hat{K}_C$ 기반 형식을 대체하기 위해 이산 해밀토니안을 $\hat{H}_1 = \begin{bmatrix} E_C & \hat{K}_R P \\ P \hat{K}_L & E_V \end{bmatrix}$ 로 구성한다.
• 결과적으로 유도된 분산 관계 $K_{L,R}(k)$ 는 복소수이며, 중심차분의 실수값 $K_C(k) = \frac{1}{\Delta} \sin(k\Delta)$ 와 달리 $k = \pm \pi/\Delta$ 에서 임의의 영점이 발생하지 않는다.
• 2밴드 및 3밴드 k·p 모델을 사용하여 일계 및 이차형식의 밴드 구조를 비교함으로써, 일계 형식에서는 단조적 성질과 임의의 상태의 부재를 확인한다.
• 3밴드 슈퍼격자 모델을 이용해 방법을 검증하였으며, 파동함수는 고파동수 진동을 보이지 않으며 에너지 밴드는 여전히 단조적임을 확인하였다.

실험 결과

연구 질문

• RQ1 왜 중심차분 근사에서 사용하는 이산 k·p 모델에서 임의의 상태가 나타나는가?
• RQ2 일계 차분 형식이 정확도나 대칭성에 영향을 주지 않으면서도 임의의 상태를 제거할 수 있는가?
• RQ3 일계와 중심차분 기울기 연산자의 분산 관계는 밴드 구조에 어떤 영향을 미치는가?
• RQ4 일계 형식은 다양한 반도체 재료와 밴드 모델에서 안정적이고 정확한가?
• RQ5 일계 형식은 비대칭 이산화의 산물들을 피하면서도 물리적 대칭성을 유지할 수 있는가?

주요 결과

• 일계 상류 차분 형식은 중심차분 형식에서 기인하는 비단조적 분산 관계로 인해 발생하는 고파동수 임의의 상태를 제거한다.
• 일계 형식의 분산 관계인 $|K_{L,R}| = \frac{1}{\Delta} \sqrt{2 - 2\cos(k\Delta)}$ 는 단조적이며 중심차분의 $K_C(k) = \frac{1}{\Delta} \sin(k\Delta)$ 보다 진짜 $k$-의존성과 더 가깝다.
• 일계 형식은 2밴드 및 3밴드 모델 모두에서 단조적인 에너지 밴드를 생성하는 반면, 중심차분 형식은 비단조적(역행적) 행동을 보인다.
• 일계 형식으로 계산된 파동함수는 고파동수 진동을 보이지 않으며, 중심차분 형식으로 계산된 파동함수는 에너지 갭을 제외한 영역에서 이러한 산물이 나타난다.
• 이 방법은 조건부 안정성이며, 매개변수 조정이 필요로 하는 다른 방법들과 달리 모델 매개변수 조정이 불필요하다.
• 일계 형식은 물리적 대칭성과 연속체 밴드 구조에 대한 충실도를 유지하며, 메조스코픽 반도체 모델의 맥락에서 이차형식 대비 더 높은 정확도를 보인다.