Девлог #67 – Хвиля 47: кінетика електродів, полімерний ланцюг & магнітні домени

Хвиля 47 одним поглядом

Кінетика електродів — рівняння Батлера-Волмера в дії

Швидкість, з якою електрохімічна реакція протікає на поверхні електрода, визначається рівнянням Батлера-Волмера, яке кодує два конкуруючих експоненційних процеси — окислення та відновлення — кожен з яких прискорюється прикладеною перенапругою η:

j = j₀[exp(αFη/RT) − exp(−(1−α)Fη/RT)]

де j₀ — густина обмінного струму (базовий струм, що протікає в обох напрямках у рівновазі), α — коефіцієнт переносу (асиметрія між прямою та зворотною енергіями активації), а η = E − E_рівн — перенапруга.

Три режими візуалізації

Праве полотно перемикається між трьома діагностичними графіками, якими постійно користуються електрохіміки:

• Крива поляризації (j від η): звичний S-подібний відгук струму на напругу, що насичується до граничного струму за великих η.
• Графік Тафеля (log|j| від η): лінеаризує рівняння Батлера-Волмера в режимі великих η. Нахил Тафеля b = RT/(αF) для анодної гілки; кінетичні параметри можна отримати графічно.
• Циклічна вольтамперометрія: застосовується лінійний розгортка потенціалу, вперед і назад, вимальовуючи характерну качкодзьобу криву ЦВ. Анодні (прямі) розгортки намальовані червоним, катодні (зворотні) — синім.

Ліва анімація електрода показує іонні частинки, що дрейфують до і від поверхні електрода, з подіями переносу електронів, що з'являються зі швидкістю, пропорційною миттєвому |j|. За високої анодної перенапруги анімація рясніє вихідними переносами електронів; за катодних потенціалів напрямок змінюється на протилежний.

Полімерний ланцюг — від випадкових блукань до скейлінгу Флорі

Найпростіша модель полімеру — це вільно-зчленований ланцюг (FJC): N сегментів однакової довжини b, кожен орієнтований ізотропно. За відсутності взаємодій виключеного об'єму (границя «тета-розчинника») статистика ланцюга ідентична випадковому блуканню, що дає середньоквадратичну відстань між кінцями R_ee = b√N та радіус інерції R_g = b√(N/6).

Коли взаємодії виключеного об'єму увімкнені (хороший розчинник), ланцюг набухає. Теорія Флорі передбачає степеневий скейлінг R_g ∝ N^ν з ν ≈ 0.588 у трьох вимірах. У поганому розчиннику ланцюг колапсує, і ν = 1/3.

Поворотне Монте-Карло

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

Права гістограма накопичує вибірки R_g протягом тисяч прийнятих рухів, накладаючи гаусову апроксимацію жовтим кольором. Гаусова апроксимація точна за тета-умов; вона помітно відхиляється в режимах хорошого та поганого розчинника, оскільки розподіл набуває негаусового хвоста.

Магнітні домени — модель Ізінга та гістерезис

Феромагнітні матеріали нижче температури Кюрі T_c спонтанно намагнічуються, але об'ємний зразок може мати нульову чисту намагніченість, тому що він поділений на домени Вейсса — області вирівняних спінів, розділені доменними стінками. 2D-модель Ізінга, точно розв'язана Онсагером у 1944 році, охоплює основну фізику:

H = −J ∑<ij> sisj − μH ∑i si

Кожен вузол ґратки несе спін s = ±1. Найближчі сусіди взаємодіють з обмінним зв'язком J > 0 (феромагнітний). Температура Кюрі дорівнює k_BT_c = 2J/ln(1+√2) ≈ 2.269J.

Рендеринг спінової ґратки в реальному часі

Спінова ґратка 80×60 рендериться на кожному кадрі за допомогою ImageData на невеликому позаекранному полотні (80×60 пікселів), потім масштабується до повного дисплейного полотна 480×300 через drawImage з imageSmoothingEnabled = false. Це дає чіткі межі доменів у стилі піксель-арт при незначних витратах GPU. Спіни «вгору» — індиго (#3f51b5); спіни «вниз» — червоні (#e53935). Доменні стінки з'являються як різкі кольорові межі, коли T значно нижче T_c; вони розмиваються та коливаються, коли T наближається до T_c, а вся ґратка стає невпорядкованим «сіль-і-перець» візерунком вище T_c.

Анімована петля B-H

У режимі автоматичного розгортання прикладене поле H лінійно змінюється від −1 до +1 і назад, записуючи пари (H, M) на кожному кроці. Накопичений слід малюється на правому полотні: пряме розгортання червоним, зворотне синім. Виникає характерна петля гістерезису — її ширина відображає коерцитивне поле H_c, її висота — намагніченість насичення M_s, а її площа — енергію, розсіяну за цикл (закон Штейнмеца ∝ f B_max^1.6).

П'ять пресетів матеріалів охоплюють діапазон від твердих магнітів (кобальт, J=2.0, широка петля) до м'яких магнітів (пермалой, J=0.8, та м'який ферит, J=0.6, з вузькими, майже лінійними петлями, оптимальними для силових трансформаторів).

Супровідний контент

Хвилю 47 супроводжують два пости блогу на суміжні теми:

• Спотлайт #48 — Матеріалознавство та інженерія: глибоке занурення в кристалічні структури, фазові діаграми сплавів, криві напруження-деформації, сплави з пам'яттю форми, зонну структуру, надпровідність та адсорбцію на поверхнях — усе з посиланнями на симулятори платформи.
• Навчання #36 — Механіка рідин та гідрогеологія: дев'ятикроковий керований виклад, що охоплює рівняння Нав'є-Стокса, число Рейнольдса, закон Дарсі, рівняння підземного потоку та адвективно-дифузійне перенесення.

Що далі

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