Як утворюються кристали — зародження, ґратки Браве та агрегація з обмеженою дифузією

Чи то досконалий діамант, похований мільярд років, чи сніжинка, що росте за секунди, — кристали постають з одного й того ж фундаментального процесу: атоми або іони самочинно організуються у періодичну ґратку, бо таке розташування мінімізує вільну енергію. Отримана геометрія кодує фізику атомних зв'язків у кожному куті й грані.

1. Пересичення та рушійна сила

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

Коефіцієнт пересичення: S = c / c* Відносне пересичення: σ = (c - c*) / c* = S - 1, де c = фактична концентрація розчиненої речовини, c* = рівноважна концентрація насичення. Різниця хімічних потенціалів: Δμ = kT \cdot ln(S) = kT \cdot ln(1 + σ). Для S = 1 (немає пересичення): Δμ = 0 \to немає рушійної сили. Для S > 1: Δμ > 0 \to кристалізація вигідна

Пересичення досягається охолодженням (для більшості солей розчинність спадає з T), випаровуванням (вода випаровується → концентрація зростає), зміною pH (для білків і біомінералів) або додаванням антирозчинника (наприклад, додаванням етанолу до водного розчину солі).

2. Класична теорія зародження

Перш ніж кристал зможе рости, має утворитися стабільний зародок. Це енергетично затратно, бо потрібно створити нову межу поділу кристал-рідина. Класична теорія зародження (CNT) врівноважує об'ємний виграш вільної енергії й вартість межі поділу:

Вільна енергія Гіббса сферичного зародка радіусом r: ΔG(r) = -(4/3)πr³ \cdot Δg_v + 4πr² \cdot γ, де: Δg_v = -(kT/Ω) ln(S) [об'ємний виграш вільної енергії, Ω = молекулярний об'єм], γ = міжфазна енергія кристал-рідина [Дж/м²]. Перший доданок: об'ємний виграш (росте як r³, сприяє великим кристалам). Другий доданок: вартість поверхні (росте як r², протидіє росту). Критичний радіус r*, де ΔG максимальна: r* = 2γΩ / (kT \cdot ln S). Бар'єр вільної енергії зародження: ΔG* = (16π/3) \cdot γ³Ω² / (kT ln S)². Класична швидкість зародження: J = A \cdot exp(-ΔG* / kT) [зародків на об'єм за секунду]

Гетерогенне зародження: на практиці зародження майже завжди починається на поверхні, частинці пилу чи іоні — де міжфазна енергія нижча, а ΔG* зменшується на множник f(θ) = (2+cosθ)(1−cosθ)²/4, де θ — крайовий кут. Саме тому прозорі розчини миттєво каламутніють, коли їх потурбувати.

CNT передбачає, що швидкість зародження надзвичайно чутливо залежить від пересичення: дворазове збільшення σ може підняти J на 10 порядків. Це пояснює, чому вище критичного порогу зародження виглядає майже миттєвим.

3. Ґратки Браве та кристалічні системи

Усі кристали — це періодичні розташування атомів. Всю повторювану структуру можна описати елементарною коміркою — найменшим повторюваним блоком, — що характеризується трьома довжинами ребер (a, b, c) і трьома кутами (α, β, γ). Існує рівно 7 кристалічних систем і 14 ґраток Браве, які вичерпують усі можливості 3D періодичної симетрії:

Система | Обмеження | Ґратки Браве | Приклад мінералу ─────────────┼────────────────────┼──────────────────────┼────────────────── Кубічна | a=b=c, α=β=γ=90° | P, I, F | NaCl, алмаз, пірит Тетрагональна | a=b\neqc, α=β=γ=90° | P, I | циркон, рутил Ромбічна | a\neqb\neqc, α=β=γ=90° | P, C, I, F | сірка, арагоніт Гексагональна | a=b\neqc, α=β=90°, γ=120° | P | кварц, лід, графіт Ромбоедрична | a=b=c, α=β=γ\neq90° | R | кальцит, турмалін Моноклінна | a\neqb\neqc, α=γ=90°\neqβ | P, C | гіпс, слюда Триклінна | a\neqb\neqc, α\neqβ\neqγ | P | польовий шпат, аксиніт P=примітивна, I=об'ємноцентрована, F=гранецентрована, C=базоцентрована, R=ромбоедрична

У поєднанні з 32 точковими групами (симетріями обертання й відбиття) 14 ґраток Браве породжують рівно 230 просторових груп — повний каталог 3D кристалічної симетрії. Рентгенівська дифракція (піки Брегга) виявляє, до якої просторової групи належить кристал.

Закон Брегга: коли рентгенівські промені відбиваються від паралельних кристалічних площин, розділених відстанню d, конструктивна інтерференція виникає при кутах, що задовольняють nλ = 2d sinθ. Вимірювання того, які кути дають яскраві плями, визначає міжплощинні відстані d, а отже й ґратку — основа рентгеноструктурного аналізу.

4. Режими росту й габітус кристала

Після зародження ріст кристала відбувається через приєднання нових атомів або іонів до наявних вузлів ґратки. Механізм залежить від пересичення:

Пошарово (Франка-ван дер Мерве): при низькому σ атоми приєднуються на сходинках (у зламах) і завершують один шар перед початком наступного. Дає плоскі, добре огранені кристали.
Спіральний ріст: гвинтова дислокація, що виходить на поверхню, створює самопідтримну сходинку, яка розкручується назовні — пояснює, чому ріст відбувається навіть при дуже низькому пересиченні.
3D-острівцевий (Фольмера-Вебера): при високому σ або поганому змочуванні зародки утворюються поверх незавершених шарів, даючи шорсткі, заокруглені грані кристала.
Дендритний ріст: при дуже високому σ дифузія стає лімітуючою стадією. Кути й вершини ростуть швидше за плоскі грані → розгалужені, фракталоподібні гілки.

Отриманий габітус кристала (зовнішня форма) відображає відносні швидкості росту різних граней. NaCl кубічного габітусу швидко росте на гранях {100}; додавання домішок сечовини вибірково адсорбується на цих гранях, уповільнюючи ріст {100} і утворюючи натомість октаедричні кристали.

5. Утворення сніжинок і дендритний ріст

Лід кристалізується в гексагональній системі (просторова група P6₃/mmc). Шестикратна симетрія снігових кристалів безпосередньо відображає гексагональну ґратку. Але чому така нескінченна різноманітність форм?

Падаючи крізь атмосферу, кристал льоду проходить крізь зони з різною температурою (−2°C до −22°C) і вологістю (пересичення 0–30%). Ріст надзвичайно чутливий до цих місцевих умов:

Режим росту залежно від температури (приблизно): 0 до -3°C: тонкі пластини, -3 до -5°C: голки / стовпчики, -5 до -10°C: порожнисті стовпчики / пластини, -10 до -22°C: секторні / зіркові пластини \to зіркові дендрити, нижче -22°C: знову пластини. Діаграма Накаї зіставляє температуру + пересичення \to 41 виявлену морфологію

Усі шість променів сніжинки переживають однакову історію температури й вологості в кожну мить (бо вони на відстані мікрометрів). Саме тому кожен промінь майже ідентичний — вони ростуть синхронно за однакових фізичних умов. Та все ж немає двох однакових сніжинок, бо кожна проходить унікальну мікротраєкторію крізь хмару.

Модель Ліббрехта: Кен Ліббрехт (Caltech) показав, що перехід між пластинчастим і стовпчастим ростом зумовлений відносною кінетикою приєднання молекул до базисної грані (вісь c) проти грані призми (вісь a). При температурах близько −5°C переважає кінетика призми → голки. Близько −15°C переважає базисна кінетика → зіркові пластини.

6. Агрегація з обмеженою дифузією (DLA)

DLA — це проста обчислювальна модель, що відтворює розгалужену морфологію дендритних кристалів, електроосаджень і мінеральних дендритів. Запропонована Віттеном і Сандером (1981):

Розмістіть частинку-зародок у центрі
Випустіть випадкового блукача здалеку
Блукач дифундує, доки не торкнеться кластера → прилипає назавжди
Повторіть — результат фрактальний агрегат

Фрактальна розмірність DLA: D \approx 1.71 (2D), \approx 2.50 (3D). Масштабування радіуса інерції: R_g ~ N^(1/D) (N = кількість частинок). Фрактал виникає тому, що: - вершини кластера «екранують» внутрішню частину від нових блукачів; - поле дифузії \nabla²φ = 0 еквівалентне рівнянню Лапласа; - швидкість межі v \propto -\nablaφ (лапласівський ріст); - та сама математика керує в'язким пальцеутворенням (комірка Хеле-Шоу)

// 2D DLA з позаґратковим випадковим блуканням
function runDLA(numParticles = 3000) {
  const cluster = [[0, 0]]; // частинка-зародок у початку координат
  const r_max = () => Math.max(...cluster.map(([x,y]) => Math.hypot(x,y)));

  for (let i = 0; i < numParticles; i++) {
    const R = r_max() + 5; // радіус випуску
    let angle = Math.random() * 2 * Math.PI;
    let x = R * Math.cos(angle);
    let y = R * Math.sin(angle);

    while (true) {
      x += (Math.random() - .5) * 2;
      y += (Math.random() - .5) * 2;
      if (Math.hypot(x, y) > R * 3) break; // втік — перезапуск

      const stuck = cluster.some(([cx, cy]) => Math.hypot(x-cx, y-cy) < 1.5);
      if (stuck) { cluster.push([x, y]); break; }
    }
  }
  return cluster;
}

Отриманий кластер візуально збігається з мінеральними дендритами (оксид марганцю на вапняку), каналами блискавки та електроосадженим цинком — усе це приклади лапласівського росту в природі.

7. Самоцвіти та мінеральні кристали

Геологічні кристали ростуть мільйони років за екстремальних тисків і температур. Умови визначають, які мінерали можуть утворитися та їхній розмір:

Алмаз: чистий вуглець, кубічна система. Утворюється на глибині 150–200 км (5–6 ГПа, ~1200°C) у літосферній мантії, потім виноситься на поверхню кімберлітовими виверженнями. Швидкість росту: ~1 мм за мільйон років.
Кварц (SiO₂): тригональна система. Росте гідротермально з багатих на кремнезем флюїдів у жилах і пегматитах. Кольори основних різновидів: аметист (Fe³⁺), цитрин (Fe⁴⁺), рожевий кварц (Ti).
Сніжинковий обсидіан: не кристал — вулканічне скло (аморфне). Кристаліти кристобаліту утворюються при відпалі, створюючи «сніжинковий» візерунок: девітрифікація.
Гігантські кристали селеніту (Печера кристалів, Мексика): гіпс (CaSO₄·2H₂O), моноклінний. Завдовжки 10–11 м, віком ~500 000 років. Виросли у гідротермальному розсолі при 58°C — трохи нижче температури, за якої ангідрит знову розчиняється в гіпс.

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

❄️ Спробуйте симуляцію росту сніжинки →