In this presentation we delineate a Gibbs free energy technique based on the extension of the mass and energy exchange (stoichiometric) matrix of the Lagrange method. As the Lagrange multipliers represent the chemical potential of the Gibbs'ian system components, the exchange matrix is developed to include other contributions than those directly involved in the (elemental) mass balance of the system. For example, interfacial energy between immiscible equilibrium phases is a feature explained by classical thermodynamics, yet it can not be included in a multi-phase calculation with mere mass balance constraints. By introducing an extended matrix with a surface area constraint, the ‘surface phase' can be consistently incorporated. Similarly, the electrochemical potential difference in an aqueous multi-phase system with two solution phases which are separated by a membrane can be taken into account by utilizing the matrix extension. In an analogous fashion, a further extension of the matrix brings about the possibility to include reaction kinetic restrictions to control the extents of chemical reactions in terms of their affinity.