Flux-Force Formalism

Building upon the local equilibrium hypothesis established previously, the dissipation function, $\mathrm{\Psi }=T\sigma$, dictates that the volumetric entropy production rate manifests as a bilinear sum of conjugate thermodynamic fluxes and forces. In the immediate vicinity of global equilibrium, the system response resides within the linear phenomenological regime. Here, the thermodynamic fluxes, $J_i$, are mathematically formalized as a Taylor expansion truncated at the first order relative to the thermodynamic forces, $X_k$. This yields the foundational mapping $J_i={\sum }_k{L}_{ik}{X}_k$, where $L_{ik}$ denotes the phenomenological conductivity matrix.

The Linear Phenomenological Regime

The phenomenological coefficients encapsulate both direct transport phenomena (where $i=k$) and interference or cross-coupling effects (where $i\mathrm{\ne }k$). In multicomponent chemical manifolds, these cross-coupling terms govern complex transport behaviors such as the Soret effect (thermal diffusion) and the Dufour effect (diffusion thermoeffect).

Transport Phenomenon	Thermodynamic Flux ($J_i$)	Conjugate Force ($X_i$)	Direct Coefficient ($L_{ii}$)
Chemical Reaction	Reaction Velocity ($v$)	Affinity ($A\mathrm{/}T$)	Reaction Rate Constant
Heat Conduction	Heat Flux (${\mathbf{J}}_q$)	Temperature Gradient ($\mathrm{\nabla }(1\mathrm{/}T)$)	Thermal Conductivity
Mass Diffusion	Mass Flux (${\mathbf{J}}_m$)	Chemical Potential Gradient ($-\mathrm{\nabla }(\mu \mathrm{/}T)$)	Diffusion Coefficient

The cross-coupling architecture ensures that a gradient in one thermodynamic domain can drive a flux in an entirely distinct domain, provided the tensorial ranks of the fluxes and forces permit coupling under Curie's principle.

Onsager Reciprocity and Broken Time-Reversal Symmetry

This linear mapping is rigorously constrained by the Onsager reciprocal relations, which mandate the strict symmetry of the cross-coupling coefficients ($L_{ik}=L_{ki}$). This symmetry is not a macroscopic coincidence but a direct consequence of microscopic reversibility. Fluctuation theory provides the statistical mechanics architecture for this equivalence, explicitly demonstrating how macroscopic transport coefficients governing coupled mass and heat transfer in binary ideal-gas mixtures emerge from the temporal regression of microscopic thermal fluctuations 1.

Crucially, the topological robustness of these reciprocal relations extends to regimes where time-reversal symmetry is explicitly broken. Recent theoretical frameworks confirm that in the presence of generic magnetic fields, the fundamental thermodynamic constraints dictated by Onsager reciprocity persist 2. Consequently, breaking time-reversal symmetry does not liberate a system from these constraints; it structurally precludes the theoretical possibility of an engine achieving Carnot efficiency while maintaining finite power output, thereby anchoring the absolute limits of heat-to-work conversion in non-equilibrium states 2.

Active Dissipation and Chemical Affinities

Within chemical kinetics, the flux-force formalism maps the reaction velocity, $v$, to the de Donder chemical affinity, $A=-\sum {\nu }_j{\mu }_j$. In the strictly linear regime, where the affinity remains significantly smaller than the thermal energy ($A\ll RT$), the reaction flux scales proportionally with the affinity. However, when applied to complex active matter, this formalism reveals highly non-trivial dissipation topologies.

Consider the active, semi-flexible polymer networks of the actin cytoskeleton. Here, myosin motors continuously drive the system out of equilibrium, generating active stresses that modulate the network between stable and contractile states. Empirical analyses of the characteristic spectrum of actin filament fluctuations permit the precise calculation of total mechanical work and entropy production 3. Contrary to naive assumptions of monotonic scaling with kinetic activity, the entropy production rate is maximized in the non-contractile, stable state of actomyosin 3.

This non-trivial maximization indicates that activity-dependent dissipation is fundamentally governed by the structural disorder inherent in the molecular interactions between actin and myosin, rather than mere mechanical throughput 3. The system dissipates maximal energy precisely when microscopic interactions are highly disordered yet macroscopically stable.

Ultimately, as chemical manifolds are driven further from equilibrium, the linear assumptions embedded in the phenomenological matrix degrade. The ensuing non-linearities precipitate symmetry breaking and pattern formation, necessitating a transition from the flux-force formalism to the rigorous evaluation of excess entropy production. This theoretical pivot directly informs the stability criteria for non-equilibrium steady states, setting the stage for the Glansdorff-Prigogine criterion.

Key Terms

• Phenomenological Coefficients — The matrix of proportionality constants ($L_{ik}$) that linearly map thermodynamic forces to their conjugate and non-conjugate thermodynamic fluxes.
• Onsager Reciprocal Relations 1 — The principle dictating the symmetry of the phenomenological cross-coupling matrix ($L_{ik}=L_{ki}$), mathematically derived from the temporal regression of microscopic thermal fluctuations.
• de Donder Affinity — The specific thermodynamic force driving a chemical reaction, defined as the negative sum of the chemical potentials of the reactants and products.