Spotlight #37 – Materiałoznawstwo: dyslokacje, diagramy fazowe, pasma półprzewodnikowe i nanomateriały

Od atomowych defektów krystalicznych po nanostruktury z ograniczeniem kwantowym — materiałoznawstwo łączy fizykę materii skondensowanej z inżynierią. Ten przegląd śledzi mechaniczne zachowanie metali poprzez mechanikę dyslokacji, mapuje przestrzeń składu stopów za pomocą diagramów fazowych, wyjaśnia, dlaczego krzem jest półprzewodnikiem o skośnej przerwie energetycznej, i kończy się na pograniczu grafenu oraz superparamagnetycznych nanocząstek.

Każde ciało stałe jest termodynamicznym kompromisem między idealnym kryształem periodycznym a strukturą losową o maksymalnej entropii. Rzeczywiste materiały są pełne defektów — wakansji, dyslokacji, granic ziaren — które nie są niedoskonałościami do wyeliminowania, lecz dźwigniami, którymi inżynierowie posługują się, aby regulować twardość, przewodność i magnetyzm. Jednocześnie zmniejszenie materiału do rozmiarów nanometrowych odsłania zjawiska kwantowe niewidoczne w skali makroskopowej: zależną od rozmiaru absorpcję optyczną w kropkach kwantowych, zależną od chiralności przewodność w nanorurkach węglowych oraz bezszczelinową dyspersję Diraca w grafenie.

1. Defekty krystaliczne i dyslokacje

Idealny kryształ znajduje się w minimum energii swobodnej tylko w 0 K; skończona temperatura gwarantuje obecność defektów punktowych. Defekty liniowe (dyslokacje) są głównym mechanizmem, dzięki któremu metale odkształcają się plastycznie: dyslokacja pozwala atomom przesuwać się pojedynczo, zamiast wymagać jednoczesnego przesunięcia całej płaszczyzny, co obniża teoretyczne naprężenie ścinające o czynnik ~1000.

Defekty punktowe, dyslokacje i wektor Burgersa

Defekty punktowe:
  Wakans: brakujący atom;  E_f ≈ 0.5–3 eV (Cu: 1.28 eV)
  Międzywęzłowy: dodatkowy atom poza węzłem sieci (własny lub obcy)
  Para Frenkla: wakans + atom międzywęzłowy w tej samej podsieci (kryształy jonowe)
  Para Schottky'ego: para wakansów kationowego i anionowego (neutralność elektryczna)
  Równowagowy udział wakansów: n_v/N = exp(−E_f / k_B T)

Typy dyslokacji:
  Krawędziowa: dodatkowa półpłaszczyzna; wektor Burgersa b ⊥ linia dyslokacji
  Śrubowa: spiralna ścieżka wokół linii; b ∥ linia dyslokacji
  Mieszana: kombinacja obu charakterów (większość rzeczywistych dyslokacji)

Wektor Burgersa b:
  Zdefiniowany przez obwód Burgersa otaczający rdzeń dyslokacji
  Błąd zamknięcia = b  (dla kryształu idealnego: b = wektor sieci)
  |b| ≈ a/√2 ≈ 0.25 nm (metale RSC, kierunek gęsto upakowany ⟨110⟩)

Systemy poślizgu (metale RSC: Al, Cu, Ni, Au):
  Płaszczyzny poślizgu: {111} (płaszczyzny gęsto upakowane)
  Kierunki poślizgu: ⟨110⟩ (rzędy gęsto upakowane) → 12 systemów poślizgu

Naprężenie Peierlsa-Nabarro: τ_PN ~ G exp(−2πd / b(1−ν))
  d = odstęp międzypłaszczyznowy, ν = współczynnik Poissona, G = moduł ścinania
  Wyjaśnia, dlaczego płaszczyzny gęsto upakowane ślizgają się preferencyjnie

Gęstość dyslokacji ρ (m−²):
  Metal wyżarzony: ρ ≈ 10¹&sup0; m−²
  Silnie odkształcony: ρ ≈ 10¹&sup5; m−²
          

2. Mechanizmy umacniania

Wszystkie mechanizmy umacniania działają poprzez utrudnianie ruchu dyslokacji. Wkłady w granicę plastyczności pochodzące od różnych mechanizmów są w przybliżeniu addytywne (reguła dodawania pierwiastków kwadratowych obowiązuje, gdy kilka mechanizmów działa jednocześnie).

Halla-Petcha, umocnienie odkształceniowe Taylora i inne mechanizmy

Umocnienie odkształceniowe (Taylor, 1934):
  σ = σ_0 + M α G b √ρ   (umocnienie Taylora)
  M ≈ 3.06 (czynnik Taylora dla polikryształu RSC)
  α ≈ 0.2–0.4 (stała oddziaływania)
  Wykładnik umocnienia odkształceniowego n w równaniu Hollomona: σ = K ε^n

Kryterium Considère'a dla powstawania przewężenia (n = ε_u):
  dσ/dε = σ   →  równomierne wydłużenie ε_u = n

Umacnianie granicami ziaren (Hall-Petch, 1951):
  σ_y = σ_0 + k_y / √d
  d = średnica ziarna, k_y ≈ 0.7 MPa·m^(1/2) dla stali miękkiej
  Metale nanokrystaliczne (d < 20 nm): możliwy odwrotny efekt Halla-Petcha

Umacnianie roztworowe:
  Δτ ∝ c^(1/2) · ε^(3/2) (model Fleischera)
  c = udział pierwiastka rozpuszczonego, ε = odkształcenie niedopasowania wynikające z różnicy rozmiarów

Umacnianie wydzieleniowe (Orowan, 1948):
  Ominięcie Orowana: Δτ_or ≈ G b / (L − 2r)  gdzie L = odstęp między cząstkami, r = promień
  Przecinanie cząstek (wydzielenia koherentne): Δτ_sh = F_max / (b L)
  Promień krytyczny cząstki r*, przy którym ominięcie staje się łatwiejsze niż przecinanie
          

3. Dwuskładnikowe diagramy fazowe

Dwuskładnikowy diagram fazowy przedstawia termodynamicznie stabilne fazy w funkcji składu i temperatury przy stałym ciśnieniu. Zawiera on ogromną ilość eksperymentalnych danych termochemicznych w formie bezpośrednio przydatnej do projektowania stopów i obróbki cieplnej.

Cechy diagramu fazowego, reguła dźwigni i przemiany

Reguła faz Gibbsa (dla stopu dwuskładnikowego przy 1 atm):
  F = C − P + 1   (C=2 składniki, ciśnienie stałe)
  F = 3 − P   → maks. 2 fazy w obszarze dwufazowym (linia wiążąca), 1 w punkcie eutektycznym

Reguła dźwigni (skład w obszarze dwufazowym α+L):
  x_0 = skład stopu,  x_α = solvus α,  x_L = likwidus
  Udział α: f_α = (x_L − x_0) / (x_L − x_α)
  Udział cieczy: f_L = (x_0 − x_α) / (x_L − x_α)

Reakcja eutektyczna (np. Pb-Sn przy 61.9% wag. Sn, 183°C):
  L → α + β   (trójfazowy punkt niezmienniczy: F=0)
  Mikrostruktura eutektyczna: odstęp lamelarny α-β λ ≈ 1–10 μm

Reakcja perytektyczna:
  L + α → β   (rzadsza; Fe-C ma perytektykę przy 0.09% wag. C)

Przemiany w stanie stałym (diagram TTT stali):
  Perlit: γ → α + Fe_3C  (dyfuzyjna; "nos" przy ~550°C)
  Bainit: szybsze hartowanie; drobniejsza dyspersja węglików w ferrycie
  Martenzyt: bezdyfuzyjna (ścinaniowa); poniżej M_s ≈ 200–400°C
    Sieć tetragonalna centrowana przestrzennie; twardość do 900 HV; bardzo krucha

Udział martenzytu (Koistinen-Marburger):
  f_M = 1 − exp[−0.011(M_s − T)]   dla T < M_s
          

4. Teoria pasmowa półprzewodników

Elektronowa struktura pasmowa kryształu wynika z periodyczności sieci. Twierdzenie Blocha gwarantuje, że stany własne w potencjale periodycznym przyjmują postać fali płaskiej modulowanej funkcją periodyczną względem sieci, a relacja dyspersji E(k) w pierwszej strefie Brillouina określa, czy materiał jest metalem, półprzewodnikiem czy izolatorem.

Twierdzenie Blocha, masa efektywna i złącze p-n

Twierdzenie Blocha:
  ψ_k(r) = u_k(r) e^(ik·r)
  u_k(r) = funkcja periodyczna względem sieci;  k = pęd krystaliczny w pierwszej SB

Masa efektywna m* (krzywizna pasma):
  1/m* = (1/ℏ²) d²E/dk²
  Dziury lekkie, dziury ciężkie w paśmie walencyjnym Si/Ge (anizotropowa m*)
  Pasmo przewodnictwa GaAs: m*_e = 0.067 m_e  (bardzo lekka → wysoka ruchliwość)

Typy przerwy energetycznej:
  Przerwa prosta (GaAs, InP): maksimum pasma walencyjnego i minimum pasma przewodnictwa przy tym samym k
    Efektywne przejścia optyczne → diody LED, lasery
  Przerwa skośna (Si, Ge): do zachowania pędu potrzebny fonon
    E_g(Si) = 1.12 eV przy 300 K;  E_g(Ge) = 0.66 eV;  E_g(GaAs) = 1.42 eV

Równowaga złącza p-n:
  Potencjał wbudowany: V_bi = (k_BT/e) ln(N_a N_d / n_i²)
  Szerokość warstwy zubożonej: W = √(2εε_0 V_bi (1/N_a + 1/N_d) / e)

Równanie diody Shockleya:
  I = I_0 (e^(eV/k_BT) − 1)
  I_0 = A e n_i² (D_n/(N_a L_n) + D_p/(N_d L_p))  (prąd nasycenia wsteczny)
  Współczynnik idealności η ≈ 1 przy ograniczeniu dyfuzyjnym; η ≈ 2 przy ograniczeniu rekombinacyjnym
          

5. Fizyka polimerów

Polimery to długołańcuchowe cząsteczki, których właściwości mechaniczne i termodynamiczne wynikają ze statystycznych konformacji szkieletu łańcucha. Model łańcucha swobodnie przegubowego oraz teoria pola średniego Flory'ego wspólnie wyjaśniają, dlaczego guma jest sprężysta, dlaczego polimery ulegają rozdzieleniu fazowemu i jak temperatura zeszklenia zależy od sztywności łańcucha.

Statystyka łańcucha, sprężystość gumy i Flory-Huggins

Łańcuch swobodnie przegubowy (N segmentów, długość b):
  Odległość między końcami: ⟨R²⟩ = Nb²
  Długość persystencji l_p: korelacja orientacji ⟨t(0)·t(s)⟩ = exp(−s/l_p)
  Łańcuch robakowaty: ⟨R²⟩ = 2 l_p L (1 − l_p/L (1 − e^(−L/l_p)))

Sprężystość gumy (sieć afiniczna) (Kuhn, 1946):
  Moduł ścinania: G = n k_B T   (n = gęstość usieciowania na jednostkę objętości)
  Naprężenie σ = G(λ − λ−²)   (λ = współczynnik rozciągnięcia, nieściśliwy)

Temperatura zeszklenia T_g:
  Równanie WLF: log(a_T/a_T_g) = −C_1(T−T_g)/(C_2+(T−T_g))
  C_1 ≈ 17.4, C_2 ≈ 51.6 K (uniwersalne stałe WLF w pobliżu T_g)
  Teoria wolnej objętości: T_g rośnie wraz ze sztywniejszym szkieletem, krótszymi łańcuchami bocznymi

Teoria Flory'ego-Hugginsa (mieszanie polimer-rozpuszczalnik):
  ΔG_mix / (nRT) = φ ln φ + (1−φ)ln(1−φ)/N + χ φ(1−φ)
  χ = parametr Flory'ego-Hugginsa (entalpia mieszania)
  χ_c = (1/2)(1 + 1/√N)² ≈ 1/2 dla dużego N (punkt krytyczny)
  Rozpad spinodalny (rozdzielenie fazowe bez zarodkowania), gdy:
    d²ΔG_mix/dφ² < 0
          

6. Nanomateriały i ograniczenie kwantowe

Gdy materiał zostaje zredukowany do rozmiarów nanometrowych, ograniczenie kwantowe przesuwa poziomy energetyczne, gwałtownie zwiększa stosunek powierzchni do objętości oraz zmienia temperatury topnienia, właściwości magnetyczne i odpowiedź optyczną. Nanostruktury umożliwiły powstanie wyświetlaczy LED z kropkami kwantowymi, tranzystorów z jednościennych nanorurek węglowych oraz czujników opartych na grafenie.

Kropki kwantowe, nanorurki węglowe, grafen i superparamagnetyzm

Ograniczenie kropki kwantowej (cząstka w sferze):
  Równanie Brusa: E_gap = E_bulk + (ℏ²π² / 2R²)(1/m*_e + 1/m*_h) − 1.8e²/(4πεε_0 R)
  R = promień kropki;  przesunięcie ku niebieskiemu przy zmniejszaniu R (kolor zależny od rozmiaru)
  Kropki CdSe: 2 nm → niebieski,  7 nm → czerwony  (obejmuje widmo widzialne)

Chiralność jednościennej nanorurki węglowej (SWCNT):
  Wektor chiralny: C = n a_1 + m a_2   (n,m liczby całkowite, a = 0.246 nm stała sieci)
  "Fotel" (n=m): metaliczna (zerowa przerwa energetyczna)
  "Zygzak" (n,0) i chiralna: metaliczna jeśli (n−m) mod 3 = 0, w przeciwnym razie półprzewodnikowa
  Przerwa energetyczna półprzewodnikowej SWCNT: E_g = 0.9 eV / d(nm)
  Wysoka ruchliwość nośników: μ ≈ 10&sup5; cm²/(V·s)

Grafen:
  Sieć plastra miodu; dwuatomowa komórka elementarna (podsieci A, B) → dwa pasma
  Dyspersja w pobliżu punktów Diraca K,K': E(k) = ± ℏ v_F |k|
  Prędkość Fermiego v_F ≈ 10&sup6; m/s (≈ c/300)
  Bezszczelinowy półprzewodnik (półmetal); przerwa otwiera się przez domieszkowanie, naprężenie, podłoże
  Ruchliwość nośników (grafen zawieszony): ≈ 2 × 10&sup5; cm²/(V·s) w temperaturze pokojowej

Superparamagnetyzm:
  Cząstki jednodomenowe poniżej promienia krytycznego r_SD ≈ 10–30 nm
  Energia termiczna randomizuje kierunek namagnesowania: temperatura blokowania T_B
  τ = τ_0 exp(K V / k_B T)   (prawo Néela-Arrheniusa)
  K = stała anizotropii, V = objętość cząstki
  Powyżej T_B: brak remanencji (superparamagnetyk);  poniżej T_B: ferromagnetyk
  Zastosowanie: środki kontrastowe do MRI (nanocząstki tlenku żelaza, d ≈ 10 nm)
          

Metoda CALPHAD: Obliczanie diagramów fazowych (CALPHAD) wykorzystuje ocenione bazy danych termodynamicznych, takie jak SGTE i TCFE, do obliczania wieloskładnikowych (do 25 pierwiastków) diagramów fazowych, których nie można wyznaczyć w pełni doświadczalnie. Energie swobodne Gibbsa są parametryzowane za pomocą wielomianów Redlicha-Kistera i optymalizowane względem danych eksperymentalnych. CALPHAD jest obecnie standardem w projektowaniu stopów turbinowych, stali i stopów wysokoentropowych; wśród komercyjnych programów znajdują się Thermo-Calc i FactSage.

Wypróbuj te symulacje