위상장 연구: 피트 패턴 Si(001) 기판의 표면 확산 구동 형태 진화

요약

이종 에피택셜 섬의 측면 정렬은 증착 전 기판의 피트 패터닝으로 일상적으로 달성됩니다. 그러나 이러한 피트의 모양, 방향, 깊이를 제어하는 것은 비자명합니다. 왜냐하면 준안정 기하학이 어닐링이나 후속 성장 중에 진화하기 때문입니다. 본 연구에서는 초기 깊이와 모양에 따른 Si(001) 상의 전형적인 피트 형태를 예측하는 연속체 위상장 모델을 제시합니다. 이 모델은 표면 확산과 실리콘 울프 형태를 재현하는 강한 이방성 표면 에너지 밀도를 포함합니다. 실험적으로 관찰된 준안정 구성을 성공적으로 재현하고, 초기 깊이에 관계없이 피트의 전체적 평탄화를 향한 공통 운동론적 경로를 밝혀냈습니다. 두 개의 서로 다른 평형 면을 가진 두 물질의 이종 에피택시로 모델을 확장하여, 예를 들어 Ge 오버그로스 중 피트 윤곽의 회전을 보여주며, 이는 최근 실험적 관찰과 일치합니다.

배경

Ge/Si 또는 InGaAs/GaAs와 같은 격자 부적합 이종 에피택시는 Stranski-Krastanow(SK) 성장 모드를 통해 3차원 섬 형성을 자주 유도합니다. 자기 조립으로 양자점을 생성할 수 있지만, 이 과정에 내재된 무작위 핵형성은 넓은 크기 및 모양 분포를 초래하여 소자 성능을 제한합니다.

지난 수십 년 동안 이종 에피택셜 성장을 정렬된 구조로 유도하기 위한 다양한 기술이 개발되었습니다. 그중 피트 패턴 기판은 섬의 높은 측면 정렬과 정밀한 크기 제어를 달성하는 가장 다재다능한 접근법 중 하나로 부상했습니다.

피트 패턴은 나노임프린트 리소그래피, 전자빔 리소그래피와 반응성 이온 식각 또는 습식 화학 식각 결합, 나노인덴테이션 등 하향식 방법으로 제작됩니다. 이 기술들은 나노미터 정밀도로 정렬된 피트 어레이 설계를 가능하게 하며, 최적화된 성장 조건에서 섬의 거의 완벽한 측면 정렬을 유도합니다.

피트 모양은 총 표면 에너지에 영향을 미치며, 결과적으로 섬 핵형성에 영향을 주므로 피트 형태의 정확한 제어가 필수적입니다. 고온에서 모세관력은 피트를 에너지적으로 유리한 구성으로 진화시키며, 때로는 완전한 치유로 이어집니다. 따라서 어닐링이나 추가 물질 증착을 일상적으로 사용하여 준안정 피트 모양을 안정화합니다. 안정화 후에도 이종 에피택셜 성장 단계에서 추가 진화가 발생할 수 있습니다.

본 연구에서는 표면 에너지 감소를 통한 표면 확산에 의해 구동되는 피트 패턴 기판의 진화를 설명합니다. 확산 제한 동역학을 이종 에피택셜 시스템에 적용한 위상장 공식을 채택했습니다. 이 모델은 현실적인 이방성 표면 에너지로 평형 모양을 재현할 수 있으며 실험 데이터로 검증되었습니다.

우리는 광범위하게 연구된 피트 패턴 Si(001) 표면에 초점을 맞춥니다. 논문 구성은 다음과 같습니다: 다음 섹션에서 위상장 모델 소개; 3절에서 Si(001) 피트 평탄화와 평형을 향한 운동론적 경로 논의; 4절에서 피트 회전을 유도하는 Ge 오버그로스 적용 제시; 결론에서 발견의 함의 논의.

방법

위상장 모델

위상장 모델은 φ=1(고체)에서 φ=0(진공)으로 매끄럽게 변하는 연속 질서 매개변수 φ를 도입합니다. φ의 진화는 표면 확산을 포착하면서 표면 에너지의 강한 이방성을 존중하는 퇴화된 Cahn-Hilliard 방정식으로 지배됩니다. 총 자유 에너지 함수형은 다음과 같습니다.

F = ∫Ω γ(n̂) [ ε/2 |∇φ|² + 1/ε B(φ) ] d𝑟 + ∫Ω β/(2ε) [ -ε∇²φ + 1/ε B′(φ) ]² d𝑟,

여기서 Ω은 시뮬레이션 영역, ε은 계면 두께, B(φ)=18φ²(1-φ)²는 이중 우물 포텐셜, β는 강한 이방성 영역에서 날카로운 모서리를 둥글게 하는 Willmore 정규화 매개변수입니다. 화학 포텐셜 μ=δF/δφ가 다음을 통해 진화를 구동합니다.

∂φ/∂t = D ∇·[ M(φ) ∇μ ],

여기서 D는 확산 계수, M(φ)=36/ε φ²(1-φ)²는 표면 이동도입니다. μ의 표현에는 이방성 표면 에너지 밀도 γ(n̂)와 Willmore 정규화 항이 포함되어 면 진화의 물리적으로 현실적인 묘사를 보장합니다.

이방성 표면 에너지

표면 에너지 밀도를 다음과 같이 모델링합니다.

γ(n̂) = γ₀ [ 1 - Σᵢ αᵢ (n̂·mᵢ)^{wᵢ} Θ(n̂·mᵢ) ],

여기서 벡터 mᵢ는 최소 에너지 결정학 방향(실리콘의 경우 〈001〉, 〈113〉, 〈110〉, 〈111〉)에 해당합니다. 각 최소의 깊이 αᵢ는 문헌에 보고된 상대 표면 에너지 값 γᵢ에서 유도되며, 폭 매개변수 wᵢ는 면의 각도 범위를 제어합니다. 실리콘의 경우 γ₀=1, α_{〈001〉}=0.15, w 값은 50(단, w_{〈113〉}=100)으로 설정했습니다. Willmore 정규화 매개변수 β는 계산적 실현성을 유지하면서 둥근 모서리를 해상하기 위해 0.005로 선택했습니다.

초기 형태 및 시뮬레이션 설정

반경 L, 깊이 H의 원형 함몰로 (001) 표면의 초기 피트 기하학을 구성하며, 폭 σ=L/2의 가우시안 꼬리로 주변 평면과 매끄럽게 연결합니다. 인터페이스로부터의 부호 있는 거리 d(r)를 쌍곡 탄젠트 프로필 φ(r)=½[1-tanh(3d(r)/ε)]에 삽입하여 초기 조건을 생성합니다. 일반적인 영역 크기 L=1(임의 단위), 계면 두께 ε=0.2를 사용했습니다. 유한 요소 툴박스 AMDiS를 반암시적 시간 적분 스킴과 인터페이스에서의 메쉬 세분화로 사용했습니다. 측면 방향에 주기적 경계 조건, 수직 축에 무유속 조건을 적용했습니다. 확산 계수 D=1로 설정하여 시뮬레이션 시간을 임의 단위로 보고합니다.

결과 및 논의

Si(001) 피트 평탄화

Figure 2는 종횡비 R=0.25와 R=0.5인 피트의 초기 시간 진화를 보여줍니다. 처음에 특징 없는 프로파일은 실리콘 표면 에너지 밀도의 최소에 해당하는 면을 자발적으로 발달시킵니다. 얕은 피트(R=0.25)의 경우, 정사각형 〈001〉 바닥이 〈113〉 에지와 좁은 〈110〉 삼각형으로 둘러싸입니다. 더 깊은 피트(R=0.5)에서는 바닥에 추가 〈111〉 면이 나타나고 〈110〉 면이 넓어집니다. 이 예측은 피트 모양 Si(001) 표면의 실험적 관찰과 정량적으로 일치합니다.

Figure 3은 더 깊은 피트(R=0.5)의 장시간 진화를 추적합니다. 원근도와 평면도의 시퀀스는 가파른 〈111〉 면의 점진적 소멸, 인접 〈113〉 면의 조대화, 〈110〉 면의 최종 소멸을 보여주며, 이들이 함께 평탄한 표면을 생성합니다. 총 표면 에너지는 Figure 4에 보인 것처럼 단조롭게 감소하며, 동역학은 모든 R≤0.5에 대해 공통 경로를 따릅니다. 1100-1200°C에서 수십 나노미터 피트의 특성 시간 척도는 수 시간 규모로, 실리콘 표면 확산에 대한 문헌 값과 일치합니다.

Ge 오버그로스로 인한 모양 변화 모방

서로 다른 평형 면 세트를 가진 물질 증착의 영향을 평가하기 위해, R≈0.1인 피트 패턴 Si(001) 기판上的 Ge 오버그로스를 시뮬레이션했습니다. t≈5.0에서의 형태(Fig. 3)에서 시작하여, Ge/Si(001) 시스템의 특징인 〈105〉 방향을 따라 깊은 최소를 포함하도록 실리콘 표면 에너지 밀도를 대체했습니다. 폭 매개변수 w_{〈105〉}를 500으로 증가시켜 얕은 〈105〉 면을 포착했습니다.

Figure 5는 표면이 급속히 재구성됨을 보여줍니다: 〈105〉 면이 남은 〈113〉 면 사이에서 핵형성되어 후자를 희생시키며 확장합니다. 최종 피트는 오직 〈105〉 면으로만 경계가 지어지며, 평면도는 초기 방향에 비해 45° 회전된 피트 윤곽을 보여줍니다. 이 회전은 Ge 증착 중 실험적으로 보고되었으며, 이제 표면 에너지 최소화만으로 발생함이 밝혀졌습니다.

결론

실리콘 표면 에너지의 강한 이방성을 명시적으로 포함하는 연속체 위상장 프레임워크를 개발하여 피트 패턴 Si(001) 기판의 표면 확산 구동 진화를 포착했습니다. 이 모델은 준안정 면 형성, 피트의 최종 평탄화, 이러한 상태로 이끄는 운동론적 경로를 재현합니다. 〈105〉 면을 포함하도록 표면 에너지 설명을 확장함으로써, Ge 오버그로스 중 관찰된 피트 윤곽 회전이 표면 에너지 감소만으로 설명될 수 있음을 보였습니다. 이 모델은 일반적이며 이방성 재매개변수화로 다른 기판에 쉽게 적용 가능하여, 나노미터 정밀도로 이종 에피택셜 섬 배치를 안내하는 피트 기하학 설계를 위한 예측 도구를 제공합니다.