Spotlight #26 — Матеріалознавство

Матеріалознавство перебуває на перетині фізики, хімії та інженерії. Будь-яку властивість матеріалу — твердість, провідність, пластичність, температуру плавлення — можна простежити до атомарної структури та дефектів. Цей огляд проведе вас крізь вісім симуляцій, які роблять ці зв'язки відчутними: ви можете обертати кристалічні ґратки, спостерігати, як атоми дифундують уздовж меж зерен, згинати віртуальний сталевий стрижень за межу плинності та досліджувати, чому кремній має заборонену зону, а мідь — ні.

1. Кристалічні структури — BCC, FCC, HCP та алмаз

Кожен твердий метал упорядковує свої атоми в одну з чотирнадцяти ґраток Браве. Три найпоширеніші — об'ємноцентрована кубічна (BCC: залізо при кімнатній температурі, вольфрам, хром), гранецентрована кубічна (FCC: алюміній, мідь, золото, нікель) та гексагональна щільно упакована (HCP: титан, магній, цинк). Структура алмазної кубічної ґратки (кремній, германій, алмаз) лежить в основі всієї сучасної електроніки.

Структура    Атомів/комірку   Коорд.число   ЩЗ      Приклад
BCC          2                8             0.680   Fe-α, W, Cr
FCC          4                12            0.740   Al, Cu, Au
HCP          6 (ефект. 2)     12            0.740   Ti, Mg, Zn
Алмаз        8                4             0.340   Si, Ge, C(алмаз)

ЩЗ = коефіцієнт щільності упаковки = об'єм атомів / об'єм елементарної комірки
FCC та HCP обидві щільно упаковані (ЩЗ 0.740, різняться лише порядком шарів ABCABC vs ABABAB)

Параметр ґратки a (FCC з атомного радіуса r): a = 2√2 · r
Міжплощинна відстань (індекси Міллера hkl, кубічна ґратка): d_hkl = a / √(h²+k²+l²)
          

3D-переглядач кристалічних структур дозволяє обертати всі шість поширених структур — включно з NaCl (іонна) та алмазом — з рендерингом за глибиною (алгоритм художника), щоб внутрішні атоми залишалися видимими. Увімкніть суперкомірку 2×2×2, щоб порахувати відстані координаційних оболонок, і спостерігайте, як панель координаційного числа оновлюється при перемиканні між структурами.

2. Атомна дифузія та закони Фіка

Атоми в твердому тілі не нерухомі — вони перескакують між вузлами ґратки через вакансії, міжвузлові положення або межі зерен. Макроскопічний опис — другий закон Фіка: ∂C/∂t = D·∇²C. Але сам коефіцієнт дифузії D залежить від температури через співвідношення Ареніуса — експоненційний енергетичний бар'єр, що виникає через енергію, необхідну для того, щоб протиснутися повз сусідні атоми.

Перший закон Фіка (стаціонарний потік):
  J = −D · ∂C/∂x        [атомів м⁻² с⁻¹]

Другий закон Фіка (залежний від часу):
  ∂C/∂t = D · ∂²C/∂x²   (1D, стала D)

Дифузійність за Ареніусом:
  D(T) = D₀ · exp(−Qd / kT)    або еквівалентно exp(−Qd / RT)
  D₀ = передекспоненційний множник [м² с⁻¹]
  Qd = енергія активації дифузії [еВ або Дж/моль]

Розв'язок для навуглецювання (напівнескінченне тіло, стала поверхнева концентрація Cs):
  (Cs − C(x,t)) / (Cs − C₀) = erf( x / (2√(Dt)) )
  де C₀ = початкова об'ємна концентрація

Типові енергії активації (самодифузія):
  Fe у Fe-α (BCC):  0.84 еВ
  Cu у Cu (FCC):    2.04 еВ
  B у Si:           3.46 еВ  (повільно — важливо для виготовлення транзисторів)
          

Симуляція дифузії в кристалі реалізує 2D-модель вакансійного механізму на квадратній ґратці, зі швидкостями перескоку, заданими законом Ареніуса. Підніміть повзунок температури і спостерігайте, як фронти дифузії поширюються крізь спочатку сегреговану структуру — профіль Фіка виникає спонтанно з мільйонів окремих атомних стрибків.

3. Механічна поведінка — напруження, деформація та руйнування

Тягніть металевий стрижень із наростаючою силою, і він спочатку пружно розтягується (закон Гука, оборотно), потім переходить у пластичну деформацію (незворотну, спричинену рухом дислокацій), деформаційно зміцнюється по мірі заплутування дислокацій, досягає межі міцності на розрив, утворює шийку і зрештою руйнується. Уся крива кодує мікроструктуру сплаву: дрібніший розмір зерна підвищує межу плинності (закон Холла-Петча), а частинки преципітату закріплюють дислокації (дисперсійне зміцнення — основа авіакосмічних алюмінієвих сплавів).

Пружна область:     σ = E · ε            (модуль Юнга E, зазвичай 70–210 ГПа для металів)
Коефіцієнт Пуассона: ν = −ε_поперечна / ε_осьова   (ν ≈ 0.3 для більшості металів)
Критерій плинності (фон Мізес): σ_ефект = √½·[(σ₁−σ₂)²+(σ₂−σ₃)²+(σ₃−σ₁)²] = σ_y

Степеневе зміцнення (пластика): σ = σ_y + K · ε_p^n
  K = коефіцієнт зміцнення, n = показник деформаційного зміцнення (0.1–0.5)

Холл-Петч (зміцнення розміром зерна): σ_y = σ₀ + k / √d
  d = середній діаметр зерна, k ≈ 0.7 МПа·м^½ (сталь)

Критерій руйнування Гриффітса:
  σ_f = √(2Eγ / πa)    для наскрізної тріщини напівдовжиною a
  В'язкість руйнування:   K_Ic = σ_f · Y · √(πa)   [МПа·м^½]
          

Симуляція «напруження-деформація» анімує зразок на розтяг з утворенням шийки та руйнуванням у реальному часі. Оберіть зі сталі (E = 200 ГПа), алюмінію (E = 69 ГПа), каучуку (гіперпружний), кістки (анізотропна) або полімеру (в'язкопружна повзучість). Площа в'язкості — інтеграл σdε — оновлюється наживо, тож можна порівняти поглинання енергії до руйнування. Симуляція механіки руйнування додає концентрацію напружень біля вершини тріщини: спостерігайте, як тріщина поширюється в міру накопичення втомного пошкодження при циклічному навантаженні.

4. Дислокації — чому метали значно слабші, ніж передбачає теорія

Теоретичне зсувне напруження, необхідне для зсуву однієї кристалічної площини відносно іншої, становить τ_th = G/(2π) ≈ 10–30 ГПа. Але реальні метали плинуть при 50–500 МПа — у 100 разів слабше. Пояснення, знайдене незалежно Тейлором, Ораваном і Полані у 1934 році, — це дислокація: лінійний дефект, у якому ковзання вже почалося, але ще не завершилося. У будь-який момент рухаються лише атоми поблизу ядра дислокації; дислокація ковзає крізь кристал, як зморшка крізь килим.

Вектор Бюргерса b: напрямок і величина спотворення ґратки
  |b| = a/2 · √(h²+k²+l²) для FCC системи ковзання ⟨110⟩{111}

Поле напружень крайової дислокації:
  σ_xx = − Gb/(2π(1−ν)) · y(3x²+y²)/(x²+y²)²
  σ_yy =   Gb/(2π(1−ν)) · y(x²−y²)/(x²+y²)²
  τ_xy =   Gb/(2π(1−ν)) · x(x²−y²)/(x²+y²)²

Напруження Пайєрлса-Набарро (ґраткового тертя):
  τ_PN ≈ 2G · exp(−2πw/b)  де w = ширина дислокації

Дислокаційне зміцнення (зміцнення Тейлора):
  Δσ = M · α · G · b · √ρ
  ρ = густина дислокацій [м⁻²], α ≈ 0.3, M = фактор Тейлора ≈ 3.06

Розмноження джерелом Франка-Ріда: дислокації вигинаються між двома точками закріплення
  і утворюють повні петлі, коли τ ≥ Gb/(2L),  L = відстань між точками закріплення
          

Симуляція ковзання дислокацій показує крайові та гвинтові дислокації, що ковзають крізь 2D-кристал під дією прикладеного зсувного напруження. Підвищуйте напруження і спостерігайте, як дислокації зароджуються, розмножуються через джерела Франка-Ріда, накопичуються біля меж зерен і врешті-решт спричиняють плинність. Кольоровий шар показує локальну концентрацію зсувного напруження навколо кожного ядра дислокації.

5. Фазові діаграми — двокомпонентні сплави

Бінарна фазова діаграма відображає, які фази (тверда, рідка, змішана) термодинамічно стабільні залежно від температури та складу. Евтектична точка — це пара склад-температура, при якій рідина замерзає безпосередньо у дві тверді фази одночасно — найнижча температура плавлення в системі. Правило важеля дає частку кожної фази у двофазній області без розв'язання жодних рівнянь.

Правило фаз Гіббса:   F = C − P + 2
  F = число ступенів вільності (інтенсивні змінні, які можна вільно змінювати)
  C = число компонентів, P = число фаз

При сталому тиску (ізобарично): F = C − P + 1

Бінарна евтектика (C=2):
  Однофазна область (рідина або тверде тіло): F = 2 (температура + склад)
  Двофазна область (ліквідус + солідус): F = 1 (зафіксуй T → обидва склади зафіксовані)
  Евтектична точка (3 фази, рідина + 2 тверді): F = 0 (інваріантна)

Правило важеля (частка фази α при складі X_0 між Xα і Xβ):
  f_α = (X_β − X_0) / (X_β − X_α)
  f_β = (X_0 − X_α) / (X_β − X_α)

Мартенситне перетворення: бездифузійне зсувне (аустеніт Fe-C → мартенсит)
  Температура початку: Ms ≈ 539 − 423(%C) − 30.4(%Mn) − 17.7(%Ni) − 12.1(%Cr) [°C]
          

Симуляція фазової діаграми сплаву будує повну бінарну евтектичну фазову діаграму в стилі Pb-Sn. Клацніть у будь-якому місці діаграми, і розрахунок за правилом важеля з'явиться автоматично — відносні довжини сегментів лінії рівноваги масштабують частки фаз. Ви також можете симулювати шлях охолодження: слідкуйте за складом вертикально вниз і спостерігайте, як мікроструктура формується на анімованій вставці.

6. Зонна структура напівпровідників

Чому кремній проводить електрику лише при легуванні, тоді як мідь проводить вільно, а алмаз не проводить взагалі? Відповідь криється в електронній зонній структурі — дозволених енергетичних рівнях для електронів у періодичній ґратці. У металах рівень Фермі перебуває всередині частково заповненої зони. У напівпровідниках заборонена зона відокремлює заповнену валентну зону від порожньої зони провідності; теплова енергія, фотони (або атоми домішки) можуть подолати цей зазор. В ізоляторах зазор занадто широкий.

Заборонена зона майже вільних електронів на межі зони при k = π/a:
  E_gap = 2|V_G|   (V_G = компонента Фур'є періодичного потенціалу)

Розподіл Фермі-Дірака (ймовірність заповнення стану E):
  f(E) = 1 / (exp((E − E_F) / kT) + 1)

Власна концентрація носіїв:
  n_i = √(N_c · N_v) · exp(−E_g / 2kT)
  N_c = 2(2πm_e*kT/h²)^(3/2),  N_v аналогічно з m_h*

Легування n-типу (донорні атоми з рівнем E_d поблизу зони провідності):
  Виморожування → насичення → власний режим при зростанні T
  n ≈ N_d (при помірній T, повна іонізація)

Коефіцієнт Холла (визначення знаку носіїв):
  R_H = −1/(ne)  (n-тип)    R_H = +1/(pe) (p-тип)

Прямий проти непрямого зазору:
  Прямий (GaAs, GaN): оптичні переходи зі збереженням імпульсу → світлодіод, лазер
  Непрямий (Si, Ge): для переходу потрібен фонон → слабкий випромінювач
          

Симуляція напівпровідникових зон будує дисперсію E–k для 1D-моделі Кроніга-Пенні та дозволяє налаштовувати глибину потенціальної ями, відстань і температуру. Хвіст Фермі-Дірака заповнюється при підвищенні температури — можна побачити перехід від зовнішнього режиму (густина носіїв задається домішками) до власного режиму (густина носіїв задається забороненою зоною та kT). Окремо панель легування показує рівні донорів n-типу та акцепторів p-типу, що розташовані одразу всередині забороненої зони.

Чому матеріалознавство заслуговує на більше симуляцій

Матеріалознавство — одна з найбільш зручних для симуляцій галузей у фізичних науках. Молекулярна динаміка, аналіз методом скінченних елементів, теорія функціонала густини, моделі фазового поля — усе це регулярно застосовується в дослідженнях матеріалів, і їхні візуальні результати одразу інтуїтивно зрозумілі. Дислокація, що рухається крізь кристал, виглядає саме так, як і має виглядати. Межа фаз, що мігрує через зерно, виглядає саме так, як і має виглядати.

Проте матеріалознавство часто відсутнє в онлайн-освіті з фізики. Більшість вступних курсів переходять безпосередньо від атомних орбіталей до термодинаміки, оминаючи кристалічне тверде тіло. Цей огляд намагається заповнити цю прогалину симуляціями, прив'язаними до реальних параметрів матеріалів і реальних рівнянь — а не іграшками.

У майбутніх хвилях: еволюція густини дислокацій під час холодної обробки, парний потенціал Леннарда-Джонса, що пов'язує міжатомні сили з пружними константами, ріст зерен при відпалі (нормальний і аномальний), а також повна симуляція p-n переходу з областю збіднення та ВАХ.