Devlog #67 – Fala 47: kinetyka elektrodowa, łańcuch polimerowy i domeny magnetyczne

Fala 47 w skrócie

Kinetyka elektrodowa — Butler-Volmer w akcji

Szybkość, z jaką przebiega reakcja elektrochemiczna na powierzchni elektrody, jest opisana równaniem Butlera-Volmera, które ujmuje dwa konkurujące procesy wykładnicze — utlenianie i redukcję — z których każdy jest przyspieszany przez przyłożone nadpotencjał η:

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

gdzie j₀ to gęstość prądu wymiany (bazowy prąd płynący w obu kierunkach w stanie równowagi), α to współczynnik przeniesienia (asymetria między energiami aktywacji reakcji w przód i w tył), a η = E − E_rów to nadpotencjał.

Trzy tryby wizualizacji

Prawy canvas przełącza się między trzema wykresami diagnostycznymi, rutynowo stosowanymi przez elektrochemików:

• Krzywa polaryzacji (j w funkcji η): dobrze znana, esowata odpowiedź prądowo-napięciowa, wysycająca się w kierunku prądu granicznego przy dużych wartościach η.
• Wykres Tafla (log|j| w funkcji η): linearyzuje równanie Butlera-Volmera w reżimie wysokich wartości η. Nachylenie Tafla b = RT/(αF) dla gałęzi anodowej; parametry kinetyczne można odczytać graficznie.
• Woltamperometria cykliczna: przykładane jest liniowe wymuszenie potencjału, przemiatane w przód, a następnie w tył, kreśląc charakterystyczną krzywą CV w kształcie „dzioba kaczki". Przemiatania anodowe (do przodu) rysowane są na czerwono, katodowe (wsteczne) na niebiesko.

Lewa animacja elektrody pokazuje jony dryfujące w stronę i od powierzchni elektrody, gdzie zdarzenia transferu elektronu pojawiają się z częstością proporcjonalną do chwilowej wartości |j|. Przy wysokim nadpotencjale anodowym animacja roi się od transferów elektronów skierowanych na zewnątrz; przy potencjałach katodowych kierunek się odwraca.

Łańcuch polimerowy — od błądzenia losowego do skalowania Flory'ego

Najprostszym modelem polimeru jest swobodnie przegubowy łańcuch (FJC): N segmentów o jednakowej długości b, z których każdy jest zorientowany izotropowo. Przy braku oddziaływań objętości wykluczonej (granica „rozpuszczalnika theta”) statystyka łańcucha jest identyczna jak w błądzeniu losowym, co daje średniokwadratową odległość koniec-koniec R_ee = b√N oraz promień żyracji R_g = b√(N/6).

Po włączeniu oddziaływań objętości wykluczonej (dobry rozpuszczalnik) łańcuch pęcznieje. Teoria Flory'ego przewiduje skalowanie potęgowe R_g ∝ N^ν z ν ≈ 0,588 w trzech wymiarach. W złym rozpuszczalniku łańcuch się zapada, a ν = 1/3.

Pivotowe Monte Carlo

Czysta dyfuzja termiczna jest niezwykle powolna w próbkowaniu przestrzeni konfiguracyjnej polimeru. Symulator wykorzystuje ruchy pivotowe: wybierane jest losowe wiązanie wzdłuż łańcucha, a cały podłańcuch po jednej stronie jest obracany o losowy kąt wokół tego punktu obrotu. Kryterium Metropolisa akceptuje nową konfigurację na podstawie wagi boltzmannowskiej zmiany energii. Dla jakości rozpuszczalnika model energii faworyzuje zmiany Rg — penalizując ekspansję w złym rozpuszczalniku, a kompresję w dobrym.

Prawy histogram akumuluje próbki R_g z tysięcy zaakceptowanych ruchów, nakładając na nie żółte dopasowanie gaussowskie. Przybliżenie gaussowskie jest dokładne w warunkach theta; odchyla się w sposób wyraźny w reżimach dobrego i złego rozpuszczalnika, gdy rozkład nabiera niegaussowskiego ogona.

Domeny magnetyczne — model Isinga i histereza

Materiały ferromagnetyczne poniżej temperatury Curie T_c samoczynnie się magnesują, jednak cała próbka może mieć zerową magnetyzację netto, ponieważ dzieli się na domeny Weissa — obszary uporządkowanych spinów oddzielone ścianami domenowymi. Dwuwymiarowy model Isinga, dokładnie rozwiązany przez Onsagera w 1944 roku, uchwytuje istotę tej fizyki:

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

Każdy węzeł sieci niesie spin s = ±1. Najbliżsi sąsiedzi oddziałują ze sprzężeniem wymiennym J > 0 (ferromagnetycznym). Temperatura Curie wynosi k_BT_c = 2J/ln(1+√2) ≈ 2,269J.

Renderowanie sieci spinów w czasie rzeczywistym

Sieć spinów 80×60 jest renderowana w każdej klatce za pomocą ImageData na niewielkim canvasie poza ekranem (80×60 pikseli), a następnie skalowana do pełnego canvasu wyświetlania 480×300 przy pomocy drawImage z ustawieniem imageSmoothingEnabled = false. Daje to ostre, pikselowe granice domen przy znikomym koszcie GPU. Spiny w górę są indygo (#3f51b5); spiny w dół są czerwone (#e53935). Ściany domenowe pojawiają się jako ostre granice kolorystyczne, gdy T jest znacznie poniżej T_c; rozmywają się i fluktuują, gdy T zbliża się do T_c, a cała sieć staje się nieuporządkowanym wzorem „sól i pieprz" powyżej T_c.

Animowana pętla B-H

W trybie automatycznego przemiatania przyłożone pole H narasta liniowo od −1 do +1 i z powrotem, rejestrując pary (H, M) na każdym kroku. Zgromadzony ślad jest rysowany na prawym canvasie: przemiatanie w przód na czerwono, wsteczne na niebiesko. Powstaje charakterystyczna pętla histerezy — jej szerokość odzwierciedla pole koercji H_c, jej wysokość magnetyzację nasycenia M_s, a jej pole powierzchni energię rozpraszaną w jednym cyklu (prawo Steinmetza ∝ f B_max^1,6).

Pięć gotowych ustawień materiałowych obejmuje zakres od magnesów twardych (kobalt, J=2,0, szeroka pętla) po magnesy miękkie (permaloj, J=0,8, oraz miękki ferryt, J=0,6, z wąskimi, niemal liniowymi pętlami optymalnymi dla transformatorów mocy).

Materiały towarzyszące

Fali 47 towarzyszą dwa wpisy na blogu poświęcone pokrewnym tematom:

• Spotlight #48 — nauka o materiałach i inżynieria materiałowa: dogłębne omówienie struktur krystalicznych, diagramów fazowych stopów, krzywych naprężenie-odkształcenie, stopów z pamięcią kształtu, struktury pasmowej, nadprzewodnictwa i adsorpcji powierzchniowej — wszystko z odnośnikami do symulatorów dostępnych na platformie.
• Learning #36 — mechanika płynów i hydrogeologia: dziewięcioetapowy przewodnik obejmujący równania Naviera-Stokesa, liczbę Reynoldsa, prawo Darcy'ego, równanie przepływu wód gruntowych oraz transport adwekcyjno-dyfuzyjny.

Co dalej

Platforma nadal uzupełnia luki w pokryciu tematycznym. Wśród priorytetowych dodatków rozważanych obecnie znajdują się trybologia (fizyka tarcia, zużycia i smarowania), propagacja fal sejsmicznych (fale P i S w ośrodkach warstwowych) oraz spektroskopia kropek kwantowych (ograniczenie kwantowe i fluorescencja zależna od rozmiaru). Po stronie bloga zaplanowane są zarówno wpis Spotlight poświęcony symulacjom nauk o Ziemi, jak i wpis Learning o mechanice statystycznej i przejściach fazowych.