Thermodynamic properties of coupled mixed-valence molecules using a cooperative PKS theory: A second-order localized—delocalized transition