Devlog #69 – Fala 49: transport błonowy, gaz Van der Waalsa & wzrost ziaren

Fala 49 w skrócie

Transport błonowy — równanie Goldmana-Hodgkina-Katza

Każdy spoczynkowy neuron znajduje się na poziomie mniej więcej −70 mV, ponieważ jego błona jest selektywnie przepuszczalna: K⁺ wycieka na zewnątrz szybciej, niż Na⁺ wnika do wnętrza. Ilościowe ujęcie tego zjawiska wymaga równania napięcia Goldmana-Hodgkina-Katza (GHK), wyprowadzonego w latach 1943–1949 poprzez niezależne rozwiązanie równania elektrodyfuzji Nernsta-Plancka dla każdego przepuszczalnego jonu, przy założeniu stałego pola elektrycznego w poprzek błony:

V_m = (RT/F) · ln[(P_Na[Na]_zewn + P_K[K]_zewn + P_Cl[Cl]_wewn) / (P_Na[Na]_wewn + P_K[K]_wewn + P_Cl[Cl]_zewn)]

Zwróć uwagę na asymetryczny znak dla Cl⁻: ponieważ niesie on ładunek ujemny, stężenia wewnętrzne i zewnętrzne są zamienione miejscami w liczniku/mianowniku względem kationów.

Stężenia jonów i potencjały Nernsta

Spoczynkowe stężenia fizjologiczne są kanoniczne:

• Na⁺: 145 mM zewnątrzkomórkowo, 12 mM wewnątrzkomórkowo → E_Na ≈ +67 mV (siła napędowa skierowana do wnętrza w spoczynku)
• K⁺: 4 mM zewnątrzkomórkowo, 140 mM wewnątrzkomórkowo → E_K ≈ −94 mV (siła napędowa skierowana na zewnątrz w spoczynku)
• Cl⁻: 120 mM zewnątrzkomórkowo, 4 mM wewnątrzkomórkowo → E_Cl ≈ −89 mV

Błona spoczynkowa znajduje się na poziomie −70 mV, ponieważ P_K ≫ P_Na w spoczynku — komórka zachowuje się niemal jak elektroda potasowa. Podczas potencjału czynnościowego bramkowane napięciowo kanały Na⁺ otwierają się przejściowo, a P_Na gwałtownie rośnie, przesuwając V_m w stronę +35 mV.

Animacja i fizyka cząstek

Symulacja przedstawia przekrój poprzeczny dwuwarstwy lipidowej z główkami fosfolipidów (ikony kółko-linia) tworzącymi dwie warstwy. Cząstki Na⁺ (niebieskie), K⁺ (złote) i Cl⁻ (zielone) poruszają się w przedziałach zewnątrz- i wewnątrzkomórkowym. Każda cząstka otrzymuje niewielki kierunkowy dryf proporcjonalny do jej elektrochemicznej siły napędowej (V_m − E_jon), co tworzy netto jednokierunkowy przepływ przez otwarte kanały za każdym razem, gdy błona jest przesunięta względem potencjału Nernsta danego jonu. Pompa Na⁺/K⁺-ATPaza jest przedstawiona jako fioletowe białko, które w jednym cyklu wymienia 3 Na⁺ na zewnątrz na 2 K⁺ do wnętrza, wnosząc niewielki prąd elektrogenny.

Gaz Van der Waalsa — konstrukcja Maxwella i punkt krytyczny

Prawo gazu doskonałego PV = nRT pomija oddziaływania międzycząsteczkowe oraz skończoną objętość cząsteczek. Van der Waals (1873) poprawił oba te założenia:

(P + a/V_m²)(V_m − b) = RT

Człon a/V_m² reprezentuje spadek ciśnienia wynikający z sił przyciągania; b to objętość wykluczona (czterokrotność objętości cząsteczki). W wysokich temperaturach izoterma maleje monotonicznie i zachowuje się gazowo. Poniżej temperatury krytycznej T_c izoterma rozwija niefizyczną „pętlę van der Waalsa” — obszar, w którym ∂P/∂V > 0 — który trzeba zastąpić poziomą linią wiązania (współistnienie dwóch faz) wyznaczaną regułą równych pól Maxwella.

Punkt krytyczny i zmienne zredukowane

Jednoczesne przyjęcie ∂P/∂V = 0 oraz ∂²P/∂V² = 0 daje:

T_c = 8a/(27Rb),  P_c = a/(27b²),  V_c = 3b

W zmiennych zredukowanych (T_r = T/T_c, P_r = P/P_c, V_r = V/V_c) równanie van der Waalsa staje się uniwersalne: wszystkie substancje spełniające założenia VdW podlegają temu samemu równaniu stanu — jest to prawo stanów odpowiadających sobie. Dlatego CO₂ i H₂O przy T/T_c = 0.9 wyglądają jakościowo podobnie na diagramie P-V.

Implementacja konstrukcji Maxwella

Znalezienie ciśnienia równowagi P_eq dla danej podkrytycznej temperatury T wymaga rozwiązania równania sześciennego. Izoterma VdW jako funkcja V_m jest obniżonym równaniem sześciennym i poniżej T_c ma trzy rzeczywiste pierwiastki (V_ciecz, V_spinodalny1, V_spinodalny2, V_para). Ciśnienie Maxwella spełnia:

∫(V_ciecz do V_para) P_VdW dV = P_eq · (V_para − V_ciecz)

Symulator wykorzystuje wyszukiwanie binarne P_eq (pomiędzy P_min a P_max pętli), a dla każdej wartości próbnej P_eq rozwiązuje równanie sześcienne (P + a/V²)(V − b) = RT trygonometryczną metodą Cardana, aby znaleźć wszystkie trzy rzeczywiste pierwiastki, a następnie ocenia bilans pól. Zbieżność do względnej tolerancji 10⁻⁶ osiągana jest w mniej niż 60 iteracjach.

Wzrost ziaren — kinetyka mikrostruktury i Hall-Petch

Gdy metal odkształcony na zimno jest wyżarzany, zgromadzona energia odkształcenia napędza migrację granic ziaren. Granice przemieszczają się w stronę swojego środka krzywizny, eliminując małe ziarna na korzyść dużych. Makroskopowe prawo (Burke i Turnbull, 1952) ma charakter paraboliczny:

d̄² − d₀² = K(T) · t,  K(T) = K₀ exp(−Q/RT)

gdzie Q to energia aktywacji migracji granic ziaren (zależna od dyfuzyjności sieciowej gatunku ograniczającego szybkość procesu), K₀ to czynnik przedwykładniczy, a R = 8.314 J/mol·K. Stal ma Q ≈ 280 kJ/mol; aluminium jedynie 142 kJ/mol, ponieważ Al ma wyższą energię błędu ułożenia, a jego granice migrują łatwiej.

Wzmocnienie Halla-Petcha

Granice ziaren blokują poślizg dyslokacji. W 1951 roku Hall i Petch niezależnie wykazali, że granica plastyczności podlega zależności:

σ_y = σ₀ + k_HP / √d̄

gdzie σ₀ to naprężenie tarcia sieciowego materiału objętościowego (niezależne od rozmiaru ziarna), a k_HP to stała materiałowa mierząca, jak trudno jest przenieść poślizg przez granicę ziarna. Dla stali k_HP ≈ 0.74 MPa·m^½; dla aluminium jedynie 0.07, ponieważ granice ziaren Al są słabszymi barierami.

Oznacza to, że za każdym razem, gdy wyżarzasz stal w 1000 °C i podwajasz średnią średnicę ziarna, zmniejszasz o połowę wkład Halla-Petcha do granicy plastyczności — jest to bezpośredni kompromis między przetwarzalnością (grube ziarna łatwiej się formuje) a właściwościami mechanicznymi (drobniejsze ziarna są mocniejsze).

Wizualna mikrostruktura

Lewe płótno wykorzystuje teselację Voronoi zasianą N₀ losowymi punktami, aby zbudować początkową strukturę ziaren. Każde ziarno otrzymuje odrębny pastelowy kolor; granice ziaren są renderowane w niemal czarnym kolorze. Podczas animacji najmniejsze ziarna są stopniowo scalane z najbliższymi sąsiadami poprzez przypisanie ich pikseli na nowo, co jakościowo odtwarza grubienie napędzane krzywizną. Prawe płótno pokazuje paraboliczną krzywą wzrostu d̄(t) oraz logarytmiczną krzywą σ_y(d) Halla-Petcha — pomarańczowa kropka na każdej z nich oznacza bieżący stan symulacji.

Co dalej

Cztery najsłabsze kategorie — biologia komórkowa, chemia fizyczna, inżynieria materiałowa i immunologia — zostały w tej fali wzmocnione. Kolejne fale będą nadal wypełniać podobne luki, z planowanymi symulacjami dotyczącymi: osmozy i ciśnienia turgorowego (biologia komórkowa), kinetyki reakcji i wykresów Arrheniusa (chemia fizyczna), propagacji pęknięć zmęczeniowych (inżynieria materiałowa) oraz aktywacji kaskady dopełniacza (immunologia).