Fizyka, technologia oraz modelowanie wzrostu kryształów
Zbigniew R. Żytkiewicz
Instytut Fizyki PAN02-668 Warszawa, Al. Lotników 32/46 tel: 22 843 66 01 ext. 3363
E-mail: zytkie@ifpan.edu.pl
Stanisław Krukowski i Michał Leszczyński Instytut Wysokich Ciśnień PAN
01-142 Warszawa, ul Sokołowska 29/37 tel: 22 88 80 244
e-mail: stach@unipress.waw.pl, mike@unipress.waw.pl
Wykład – 2 godz./tydzień – wtorek 9.15 – 11.00 Interdyscyplinarne Centrum Modelowania UW
Budynek Wydziału Geologii UW – sala 3075
http://www.icm.edu.pl/web/guest/edukacja
http://www.unipress.waw.pl/~stach/wyklad_ptwk_2009
II. semestr – Wstęp
Wzrost kryształu – dwa etapy
(4) (5) (6) (2) (1) (3)transport objętościowy substancji
procesy powierzchniowe
wolniejszy z 2 etapów transportu masy
decyduje o szybkości wzrostu kryształu
dyfuzja, konwekcja (swobodna i wymuszona)
wykłady 5 i 7
wykłady 4 i 6
transport ciepła
3
4 5
2
1
Ø – energia wiązania sąsiednich ścian zysk energetyczny: <1> Ø <2> 2 × Ø <3> 3 × Ø <4> 4 × Ø <5> 5 × Ø
Atom na powierzchni
koncentracja <5> i <4> maleje z czasem pozycja <3> sama się odtwarza
wzrost na załamaniach (kinkach) stopni
Vst Vgr
⇑
→
⇑
V
st 0λ
pr ze sy ce niemod płynących stopni wzrostowych (step flow mode)
Skąd wziąć stopnie ?
dyslokacje (wykład 9)
dyslokacja śrubowa – permanentne źródło stopni
Skąd wziąć stopnie ?
dezorientacja (miscut) powierzchni
⇑
⇑→
→
⇑
→
⇓
λ
0V
stopniaV
grα
αtaras
stopień
α
λ
tg
c
2
0=
c/2 – wysokość stopnia
STM Si(100)
Swartzentruber PRB 1993V
grnie zależy od α, bo:
t
c
V
gr2
=
const
V
t
stopnia=
=
λ
0 0λ
∝
stopniaV
Skąd wziąć stopnie ?
dwuwymiarowe zarodkowanie
przesycenie większe niż krytyczne - zbyt duży flux cząstek
- zbyt mała ruchliwość powierzchniowa - ….
to widać w RHEED
wykład II-8oscylacje RHEED brak oscylacji RHEED w step flow mode 0
λ
>
dyfL
0λ
<
dyfL
birth & spread model
cr
r
r
>
Ohara and Reid, Modelling Crystal Growth Rates from Solution (1973)
Dyfuzja
F C(t)t
C
F
div
∂
∂
−
=
C
grad
D
F
=
−
Ft
C
C
D
∂
∂
=
∇
2 równanie dyfuzji wykład 7 1. Prawo Ficka+ warunek początkowy
+ warunki brzegowe
)
,
,
,
(
x
y
z
t
C
Transport objętościowy – dyfuzja
Przykład: wzrost A z roztworu w B (B – rozpuszczalnik)
np. Si w GaZałożenie: rozpuszczalność A w B rośnie z T + równowaga fazowa
(∞ szybka kinetyka powierzchniowa)
C – koncent racja A T T1 T2 kryształ A źródło A t = 0, T(z) = T1 C z C(T 1 ) roztwór A w B
nuda; nic się nie dzieje
stan początkowy
kryształ A źródło A t > 0 C z C(T 1 ) roztwór A w B T1 T2 C(T 2 )dyfuzyjny transport cząstek
wzrost kryształu
stan końcowy
Uwagi: - pełne rozwiązanie C(z,t) z równania dyfuzji
- orientacja układu nieważna
-stan stacjonarny 0 = ∂ ∂ ∝ z gr z C D Vj
dyfKonwekcja naturalna
)
,
( C
T
ρ
ρ
=
konwekcja
termiczna
stężeniowa
konwekcja naturalna –
przepływ w cieczy wywołany lokalnymi zmianami gęstości i obecnością pola grawitacyjnego0 < ∂ ∂ T ρ 0 < ∂ ∂ C ρ
typowe roztwory III-V ρsolvent> ρsolute (Ga, In) (As, P)
wykład 7 kryształ A źródło A roztwór A w B T1 T2
j
dyf T2 > T1 C(T 2 )C z C(T 1 ) L 0 C(z)T2 > T1 → ρ(L) < ρ(0) stabilizujący termiczny rozkład ρ(T) C2 > C1 → ρ(L) < ρ(0) stabilizujący stężeniowy rozkład ρ(C)
pole grawitacyjne ρ(T 1 ) ρ(T) z ρ(T 2 ) ρ(C 1 ) ρ(C) z ρ(C 2 )
roztwór lekki na górze
roztwór ciężki na dole ρT(z)
ρ
C(z)
Konwekcja naturalna
wykład 7)
,
( C
T
ρ
ρ
=
konwekcja
termiczna
stężeniowa
0 < ∂ ∂ T ρ 0 < ∂ ∂ C ρ źródło A T1 T2j
dyf C(T 2 ) T2 > T1 C z C(T 1 ) L 0 C(z)T2 > T1 → ρ(L) > ρ(0) destabilizujący termiczny rozkład ρ(T) C2 > C1 → ρ(L) > ρ(0) destabilizujący stężeniowy rozkład ρ(C)
pole grawitacyjne ρ(T 1 ) ρ(T) z ρ(T 2 ) ρ(C 1 ) ρ(C) z ρ(C 2 )
roztwór ciężki na górze
roztwór lekki na dole ρT(z)
ρ
C(z)
układ niestabilny – przepływ konwekcyjny w objętości cieczy
kryształ AKonwekcja naturalna
wykład 7)
,
( C
T
ρ
ρ
=
konwekcja
termiczna
stężeniowa
0 < ∂ ∂ T ρ 0 < ∂ ∂ C ρ C(T2) C z C(T1) L 0 pole grawitacyjne T1 źród ło A T2j
dyf T2 > T1 kry sz ta łA C(z) ρT(z)ρ
C(z)
ρ(T1) ρ(T) z ρ(T2) ρ(C1) ρ(C) z ρ(C2) źród ło A T2j
dyf T2 > T1 kry sz ta łAukład niestabilny – przepływ
konwekcyjny w objętości cieczy
PYTANIE 1:
Znajdź rozkład gęstości roztworu
wywołany rozkładem temperatury i
koncentracji
PYTANIE 2:
Czy w roztworze pojawi się przepływ
konwekcyjny, a jeśli tak to jaki?
Konwekcja wymuszona
dyfuzyjny transport cząstek
kryształ A źródło A C z C(T 1 ) roztwór A w B T1 T2 C(T 2 )
j
dyf T2 > T1t = ∞
↑
⇔
↑
∂
∂
grV
z
C
Rzeczywistość:
nie można uniknąć grad T w poziomie
↓
niekontrolowany przepływ
kryształ A źródło A C z C(T 1 ) T1 T2 C(T 2 ) T2 > T1t = ∞
dyfuzyjny + konwekcyjny transport cząstek
strefa wymieszana (jednorodna) strefa dyfuzyjnaω
2 1 6 1 3 1ν
ω
δ
=
const
⋅
D
δ - grubość warstwy dyfuzyjnej ν - lepkość
ω - prędkość kątowa
Burton, Prim, Schlichter, J. Chem. Phys. 21 (1953) 1987. Zalety mieszania mechanicznego:
• zwiększenie prędkości wzrostu
• kontrola ewentualnych naturalnych przepływów w cieczy • większe ryzyko przesycenia stężeniowego
Koncepcja strefy dyfuzyjnej:
roztwór = strefa doskonale wymieszana + strefa transportu dyfuzyjnego
Następne wykłady – pakiet I
I. Wzrost kryształów objętościowych
1. Wstęp - mechanizmy transportu masy i ciepła
Z.R. Żytkiewicz (16 luty 2010)
2. Wzrost kryształów objętościowych z roztopu
T. Słupiński (23 luty 2010)
3. Metody wzrostu objętościowego z roztworu
T. Słupiński (2 marzec 2010)
4. Wzrost objętościowy z fazy gazowej
K. Grasza (9 marzec 2010)
I. Wzrost kryształów objętościowych
Cz Si – Ø 450 mm (1000 kg kryształu !!!)
wzrost z roztopu (melt)
główne cechy:
- skład kryształu = skład cieczy
- wzrost kontrolowany przechłodzeniem
- najważniejszy transport ciepła; transport masy tylko o tyle o ile wpływa na rozkład T
- duże prędkości wzrostu ~cm/h (ważna emisja ciepła krystalizacji)
I. Wzrost kryształów objętościowych
wzrost z roztworu (solution)
główne cechy:
- temperatura wzrostu poniżej punktu topnienia materiału
- bardzo szeroki zakres materiałów jeśli istnieje dobry rozpuszczalnik (najlepiej jeśli to składnik kryształu) - kluczowy dobór rozpuszczalnika
GaN z roztworu Ga-N
2 mm
- wzrost kontrolowany przesyceniem
- najważniejszy transport masy w roztworze; transport ciepła znacznie szybszy niż transport masy - średnie prędkości wzrostu < mm/h
GaN z roztworu w nadkrytycznym NH3
T ~500 – 600oC; p ~3000 atm
roztwór zasadowy (mineralizatory KNH2; NaNH2, …)
1 cm
Ammono Unipress
I. Wzrost kryształów objętościowych
wzrost z fazy gazowej
główne cechy:
- dla materiałów, które posiadają „rozsądną” prężność par składników - możliwość wspomagania transportu reakcją chemiczną
Physical vapor transport
albo
Chemical Vapor Transport
T2 > T1
Następne wykłady – pakiet II
II. Techniki wzrostu epitaksjalnego:
5. Epitaksja – wstęp
Z.R. Żytkiewicz (16 marzec 2010)
6. Wzrost epitaksjalny z fazy ciekłej (LPE)
Z.R. Żytkiewicz (23 marzec 2010)
7. Epitaksja z fazy gazowej (VPE)
M. Leszczyński (30 marzec 2010)
8. Epitaksja z wiązek molekularnych (MBE) Z.R. Żytkiewicz (13 kwiecień 2010)
kryształy objętościowe
Tranzystor polowy FET
(field effect transistor)
implantacja lub dyfuzja
domieszki
kryształy objętościowe
ogniwo słoneczne
dioda laserowa
struktury epitaksjalne
Potrzebne metody wytwarzania cienkich (~10 nm) warstw krystalicznych
n-AlGaAs metal (e.g. aluminum)
ohmic ohmic
source gate drain
Insulating substrate i-GaAs i-AlGaAs tb δ 2DEG
HEMT
2 DE G EFdefinicje
epitaksja = nakładanie warstw monokrystalicznych na monokrystaliczne podłoże wymuszające strukturę krystaliczną warstwy
II. Techniki wzrostu epitaksjalnego
epitaksja
zarodki nowej warstwy
podłoże niedopasowanie sieciowe - wzrost pseudomorficzny - relaksacja niedopasowania - dyslokacje niedopasowania podłoże as warstwa aerelax s s e
a
a
a
f
=
(
−
)
/
niedopasowanie sieci ( lattice misfit) e elf
h
E
∝
2⋅
warstwa naprężonametody redukcji defektów w niedopasowanych sieciowo warstwach epitaksjalnych naprężenia termiczne Tgr GaAs Si RT naprężenia rozciągające w GaAs/Si GaAs Si
α
α
<
II. Techniki wzrostu epitaksjalnego
Vgr ∼ μm/min (17 nm/s)
ruch suwaka grafit
podłoże suwak grafitowy ciekłe
roztwory
epitaksja z fazy ciekłej LPE
wzrost z roztworu
T
epitaksja z fazy gazowej VPE
NH3
H2
GaCl3
GaCl3 + NH3 → GaN + 3HCl
epitaksja z wiązek molekularnych MBE
Vgr∼ μm/h (0.3 nm/s = 1 ML/s)
duża precyzja wzrostu b. cienkich warstw