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
Struktury krystaliczne
Wizualizator 3D sieci BCC, FCC, HCP i innych sieci Bravais z oznaczeniami płaszczyzn i kierunków Millera.
Diagram fazowy stopu
Dwuskładnikowe układy eutektyczne i perytektyczne: przeciągnij punkt skład/temperatura, aby odczytać udziały fazowe wg reguły dźwigni na żywo.
Wzrost kryształu
Kinetyczny model wzrostu powierzchni metodą Monte Carlo z regulowaną temperaturą i szybkością osadzania, przejście od wzrostu warstwowego do wyspowego.
Płyn Lennarda-Jonesa
Symulacja MD atomów gazu szlachetnego oddziałujących poprzez potencjał 6-12; przejście fazowe, funkcja korelacji parowej i wyświetlanie ciśnienia.