How Multivalency controls Ionic Criticality
  To understand how multivalency influences the reduced critical temperatures, Tce (z), and densities, roce (z), of z : 1 ionic fluids, we study equisized hard-sphere models with z = 1-3. Following Debye, Hueckel and Bjerrum, association into ion clusters is treated with, also, ionic solvation and excluded volume. In good accord with simulations but contradicting integral-equation and field theories, Tce falls when z increases while roce rises steeply: that 80-90% of the ions are bound in clusters near T_c serves to explain these trends. For z \neq 1 interphase Galvani potentials arise and are evaluated. 
Michael E. Fisher, Jean-No\"el Aqua and Shubho Banerjee