기저 상태, 엔벨로프 함수 및 다중대역 해밀토니안

벌크 k·p 해밀토니안을 양자우물 문제로 변환하려면 두 가지 요소가 필요하다. 하나는 벌크 행렬이 공간적으로 변화하는 구조를 기술할 수 있게 해주는 엔벨로프 함수 근사법이며, 다른 하나는 해당 행렬을 표현하기 위한 일관된 기저이다. 본 페이지에서는 이 두 요소를 종합하여 OghmaNano에서 사용되는 정확한 기저 순서를 기록함으로써 각 모델별 페이지들이 단일 정의를 기준으로 참조할 수 있도록 한다.

1. 세 종류의 함수

혼동하기 쉬운 세 가지 개체를 구분해두는 것이 중요합니다. 블로흐 함수\(\psi_{n\mathbf{k}} = e^{i\mathbf{k}\cdot\mathbf{r}} u_{n\mathbf{k}}\)는 결정체의 정확한 고유상태입니다. 브릴루앙 영역 중심에서의 셀주기 함수\(u_{n0}\)는 고정된 기저함수 중 하나로서 전도대에서는 s형태를, 가전대에서는 p형태(\(X,Y,Z\))를 띱니다. 엔벨로프 함수\(F_n(z)\)는 각 기저함수에 곱해지는 완만하게 변하는 진폭으로서 양자구속 효과를 나타냅니다. 이종구조에서는 전체 상태가 여러 기저함수들의 합으로 표현됩니다.

\[ \Psi(\mathbf{r}) = e^{i\mathbf{k}_\parallel\cdot\boldsymbol{\rho}} \sum_n F_n(z)\, u_{n0}(\mathbf{r}), \]

해석기가 계산하는 것은 바로 이러한 엔벨로프 함수들의 집합 \(\{F_n(z)\}\)입니다.

2. 엔벨로프 함수 근사법

엔벨로프 함수 근사는 두 가지 가정을 전제로 한다. 첫째, 모든 층에서 동일한 영역중심 기저 함수 \(u_{n0}\)를 사용할 수 있다는 점이며 둘째, 재료의 물성값(밴드갭, 질량, 케인 및 루팅거 파라미터 등)이 격자 상수보다 훨씬 큰 규모에서만 위치에 따라 변화한다는 점이다. 이에 따라 벌크 영역의 k·p 행렬은 엔벨로프 함수에 작용하는 미분 연산자로 변환되며 그 계수 역시 위치에 따라 달라진다. 핵심적인 변환 과정은 성장 방향의 파동벡터가 미분 연산자로 대체되는 것이다.

\[ k_z \;\longrightarrow\; -i\,\frac{d}{dz}, \]

면내 방향의 파동벡터 \(k_x, k_y\)들은 구조가 면 내에서 병진 대칭성을 가지므로 여전히 양호한 양자수로 작용합니다. 벌크 해밀토니안 내의 모든 \(k_z\) 항은 이 연산자로 대체되며, 모든 \(k_z^2\) 항은 이전 페이지에 설명된 대로 대칭적인 2차 미분 연산자로 치환됩니다. \(k_x k_z\)와 같은 곱셈 항들은 \(k_x(-i\,d/dz)\) 형태로 변환되는데, 이때 계수가 변하더라도 연산자가 헤르미티안 성질을 유지하도록 대칭화 처리가 이루어집니다. 그 결과로 엔벨로프 함수들에 대한 행렬 미분방정식이 도출되며, 이는 유한차분 페이지에서 설명된 방식대로 이산화됩니다.

3. 블록 형태의 해밀토니안 읽는 방법

다중대역 해밀토니안은 개별 성분들을 검토하기에 앞서 소수의 블록들로 나누어 해석하는 것이 가장 적절합니다. 가장 기본적인 분할 방식은 전도대와 유사한 상태들과 가전대와 유사한 상태들을 구분하는 것입니다.

\[ H = \begin{pmatrix} H_{cc} & H_{cv} \\[2pt] H_{vc} & H_{vv} \end{pmatrix}, \]

여기서 \(H_{cc}\)은 s형 전도대 상태들을, \(H_{vv}\)는 p형 가전대 상태들을, 그리고 \(H_{cv}=H_{vc}^\dagger\)은 이들 사이의 케인 운동량 결합을 나타낸다. 전도대 블록은 본질적으로 비포물선형 유효질량 연산자이며, 가전대 블록에는 무거운정공, 가벼운정공 및 스플릿오프 현상에 관한 정보가 포함된다. 결합 블록은 전도대의 비포물선적 특성을 가전대 구조와 연결함으로써 이 모델을 진정한 의미에서 ‘다중대역’ 모델로 만드는 핵심 요소이다. 응력 및 스핀-궤도 효과 항들도 각각 적절한 블록에 추가된다.

4. OghmaNano에서 사용되는 기저 순서

해밀토니안 항목의 수치값은 기저 상태들이 나열되는 순서에 따라 결정됩니다. OghmaNano는 아연석 및 석회석 구조 모델에 대해서는 고정된 직교좌표계 기저를 사용하며, Luttinger–Kohn 모델에는 총각동량 기저를 활용합니다. 기저 순서는 다음과 같습니다:

• 8-밴드 아연석 구조 및 8-밴드 석회석 구조: \(\{\,|S\!\uparrow\rangle,\ |S\!\downarrow\rangle,\ |X\!\uparrow\rangle,\ |Y\!\uparrow\rangle,\ |Z\!\uparrow\rangle,\ |X\!\downarrow\rangle,\ |Y\!\downarrow\rangle,\ |Z\!\downarrow\rangle\,\}\)
• 10-밴드 아연석 구조: 위에 언급된 8개의 상태에 두 개의 국소화된 질소 상태가 추가된 형태입니다 \(\{\,|N\!\uparrow\rangle,\ |N\!\downarrow\rangle\,\}\).
• 6-밴드 루팅거-콘 모델: \(\{\,\text{HH}\!\uparrow,\ \text{LH}\!\uparrow,\ \text{LH}\!\downarrow,\ \text{HH}\!\downarrow,\ \text{SO}\!\uparrow,\ \text{SO}\!\downarrow\,\}\)
• 실리콘 Δ-밸리: 각 밸리마다 단일(스칼라) 밴드가 존재하며, 6개의 밸리 방향 각각에 대해 한 번씩 해석됩니다.

직교좌표계 기반의 \(|X\rangle,|Y\rangle,|Z\rangle\) 오비탈을 사용하는 것은 총각동량 기반의 \(|J,m_j\rangle\) 조합을 사용하는 것보다 운동량 및 변형 결합을 실수 형태로 간단히 구성할 수 있게 해주지만, 그 대가로 스핀-궤도 결합은 명시적인 비대각 항으로 남게 됩니다. 반면 Luttinger–Kohn 모델에서는 \(|J,m_j\rangle\) 기저를 채택하여 스핀-궤도 결합을 대각 항으로 만들어주며, 운동량 관련 항들은 전형적인 \(P,Q,R,S\) 형태를 갖게 됩니다. 이 두 기저 선택 방식은 물리적으로 동등하며 각각 특정 모델에 적합하므로 각각의 설명 페이지에서 상세히 기술되어 있습니다.