Generalized Gibbs Phase Rule for Nonequilibrium Steady States

Extending the fundamental phase rule F = C - P + 2 to systems maintained out of equilibrium by external driving forces

F_neq = C + N_A - P + 2
C = components, N_A = affinities, P = phases
764
Coexistence Points
2.15
NESS Manifold Dim.
1.30
EQ Manifold Dim.
17
Critical Line Points

Key Findings

1

Each affinity adds one degree of freedom

The coexistence manifold dimension grows linearly with the number of independent nonequilibrium driving parameters N_A. For C=1, P=2: F=1 (EQ) versus F=2 (NESS with N_A=1).

2

Critical line replaces critical point

In equilibrium, the critical point is isolated (F=0). Under driving, T_c(A) traces a one-dimensional curve in the (T, A) plane, consistent with F_crit = 1.

3

Current matching is a new constraint

At steady-state interfaces, particle currents must be continuous. This extra coexistence condition is absent in equilibrium, partially offsetting the enlarged state space.

4

Manifold dimensionality confirmed

Correlation dimension analysis yields dim=2.15 for the NESS manifold (predicted: 2) and dim=1.30 for the equilibrium slice (predicted: 1).

Interactive Phase Rule Calculator

Compute Degrees of Freedom

1
1
2
Degrees of Freedom
2
F = 1 + 1 - 2 + 2 = 2

Degrees of Freedom vs. N_A

Mean-Field Coexistence Surface

Critical Line T_c(A)

Free Energy Landscape

0.0

KLS Monte Carlo Simulation

Steady-State Current vs. Field

Density Profiles at T=1.2

Manifold Dimensionality

Predicted vs. Measured

Constraint Counting Summary

Quantity Equilibrium NESS (N_A=1) NESS (N_A=2)
Components C111
Affinities N_A012
Phases P222
Unknowns/phase D234
Constraints/interface K345
Degrees of Freedom F123

Critical Temperature Data

Affinity A T_c(A) Shift from EQ Status
0.002.0030.0%Confirmed
0.122.0030.0%Confirmed
0.241.996-0.4%Confirmed
0.371.988-0.8%Confirmed
0.491.974-1.5%Confirmed
0.611.959-2.2%Confirmed
0.731.937-3.3%Confirmed
0.861.908-4.8%Confirmed
0.981.878-6.3%Confirmed
1.101.841-8.1%Confirmed
1.221.797-10.3%Confirmed
1.351.746-12.9%Confirmed
1.471.680-16.1%Confirmed
1.591.606-19.8%Confirmed
1.711.518-24.2%Confirmed
1.841.401-30.1%Large shift
1.961.202-39.9%Large shift

KLS Phase Diagram Data

Order Parameter (density-profile variance) across (T, E) space