The Nernst Equation — Single-Ion Equilibrium
Start with the equilibrium potential for a single ion species crossing a membrane. At electrochemical equilibrium, the concentration gradient (chemical potential) exactly balances the electrical potential difference:
E_ion = (RT / zF) · ln([X]_out / [X]_in)
At 37 °C (310 K) and z = +1, the Nernst factor RT/F ≈ 26.7 mV × 2.303 ≈ 61.5 mV per decade of concentration ratio.
Open the Membrane Transport simulator. In the stats panel, read E_Na, E_K, and E_Cl at resting conditions (T = 310 K). Confirm that E_K ≈ −94 mV (K⁺ wants to leave the cell) and E_Na ≈ +67 mV (Na⁺ wants to enter).
The Goldman-Hodgkin-Katz Equation — Multiple Ions
When two or more ion species are present, the membrane potential is a permeability-weighted average of their Nernst potentials — not a simple arithmetic mean:
V_m = (RT/F) ln[(P_Na[Na]_out + P_K[K]_out + P_Cl[Cl]_in) / (P_Na[Na]_in + P_K[K]_in + P_Cl[Cl]_out)]
At rest, P_K dominates (P_K : P_Na : P_Cl ≈ 1 : 0.04 : 0.45), so V_m sits close to E_K at −70 mV.
In the simulator, slowly drag P_Na from 0.01 to 5.0. Watch V_m increase from −70 to near +35 mV as the P_Na/P_K ratio rises from 0.04 to 5. This is exactly what happens when voltage-gated Na⁺ channels open during an action potential.
Action Potential Presets — Hodgkin-Huxley in Practice
The Hodgkin-Huxley model (1952) describes how time-varying Na⁺ and K⁺ conductances produce action potentials. The GHK equation gives a static snapshot at each moment.
Use the four presets: Resting (−70 mV, K⁺-dominated), Action Potential Peak (+35 mV, Na⁺-dominated), Fully Depolarized (0 mV, equal permeabilities), and Hyperpolarized (−90 mV, very high P_K). For each, note which ions have inward vs outward driving forces.
The driving force on Na⁺ is (V_m − E_Na) : at −70 mV this is −137 mV (strongly inward); at +35 mV it is approximately 0 (near reversal potential). This is why Na⁺ current self-limits at the action potential peak even while Na⁺ channels remain open.
From Ions to Molecules — Ideal vs Real Gas
The Nernst equation contains a logarithm of a concentration ratio — exactly like the entropy term in an ideal free-energy change: ΔG = RT ln(Q/K). Ideal behaviours (ideal gas, ideal solution, Raoult’s law) all share this logarithmic form.
Now open the Van der Waals Gas simulator. At T/T_c > 1.5, the P-V isotherms look like ideal-gas hyperbolas. The van der Waals equation reduces to PV ≈ RT at high V_m when the correction terms a/V_m² and b/V_m become negligible.
This is the first connection: both the Nernst equation and ideal-gas law emerge from dilute, non-interacting limits. When interactions matter — either between ions across a membrane or between gas molecules near their boiling point — you need corrections.
The Van der Waals Equation — Cubic and Its Consequences
The VdW equation is cubic in volume at fixed T and P:
PV_m³ − (Pb + RT)V_m² + aV_m − ab = 0
Below T_c this has three positive real roots: V_liquid, and V_vapour (and an unphysical middle root). The critical conditions ∂P/∂V = 0 and ∂²P/∂V² = 0 give T_c = 8a/(27Rb), P_c = a/(27b²), V_c = 3b.
In the simulator, set the H₂O preset and move T/T_c from 1.2 to 0.7. Watch the isotherm develop a local maximum and minimum (the “van der Waals loop”). The region between the two extrema is mechanically unstable.
Maxwell Equal-Area Construction — Finding the Phase Boundary
The physical two-phase coexistence pressure P_eq is not the VdW minimum — it is the pressure at which the areas above and below the horizontal tie line are equal:
∫(V_liq to V_vap) P_VdW dV = P_eq · (V_vap − V_liq)
This equal-area condition is equivalent to equal Gibbs free energy (equal chemical potential) of the liquid and vapour phases. Thermodynamically, it is the condition that the two phases are in equilibrium — they can coexist without net interconversion.
Enable “Maxwell Construction” in the simulator. The flat horizontal segment replaces the loop. Move T/T_c from 0.5 to 1.0 and observe how both liquid and vapour volumes approach V_c as the critical point is reached. At T = T_c, the loop shrinks to the inflection point and the two phases become identical.
Arrhenius Kinetics — The Bridge to Materials
Both the van der Waals critical point and grain growth share a common mathematical ancestor: Arrhenius kinetics. Any thermally-activated process has a rate that depends exponentially on temperature:
k = A · e^(−E_a / RT)
For grain growth, the rate constant is K(T) = K₀ exp(−Q/RT) where Q is the activation energy for grain-boundary diffusion. Just as Na⁺ channels have an activation voltage V_half, grain boundaries have an activation energy Q that must be overcome by thermal fluctuations.
In the Grain Growth simulator, change Temperature from 400 °C to 1200 °C (steel preset) and observe the growth rate stat K(T). It changes by many orders of magnitude — at 1200 °C steel grain boundaries migrate thousands of times faster than at 400 °C.
Grain Growth Kinetics — The Parabolic Law
Normal grain growth during annealing follows a parabolic law (Beck 1952):
d̄² − d₀² = K(T) · t
This arises because the driving force for boundary migration is inversely proportional to the grain radius (boundary curvature), while the migration velocity is proportional to driving force. Combining these gives a differential equation d(d̄)/dt ∝ 1/d̄, which integrates to d̄² − d₀² = const.
Using the simulator, press “+1 h Anneal” five times at 900 °C, then five more times at 1100 °C. The first five hours doubled d̄ by a small amount; the second five hours (at higher T) produce a much larger jump because K(T) grows exponentially with temperature.
Hall-Petch Strengthening — Structure Controls Properties
The Hall-Petch relation connects microstructure to mechanical performance:
σ_y = σ₀ + k_HP / √d̄
The 1/√d̄ dependence arises from dislocation pile-up theory: grain boundaries act as barriers to dislocation motion. Finer grains mean more grain boundaries per unit volume, more pile-up stress, and higher yield strength.
Observe the Hall-Petch plot in the simulator (bottom-right panel). The orange dot moves along the curve as grain size evolves. Notice that the slope is steep at small d̄ (large Hall-Petch benefit for small further refinement) and nearly flat at large d̄ (diminishing returns from further coarsening).
Synthesis: You have now traced a line from the logarithmic ion equilibrium in a cell membrane, through polynomial real-gas equations and phase transitions, to exponential kinetics and power-law structural mechanics. All three areas — biophysics, physical chemistry, and materials science — use the same mathematical tools: logarithms of ratios, cubic roots, and exponential Boltzmann factors.