The stability of charged interfaces varies as a function of the nature of their neutralizing counterions. Changing the counterions from monovalent to divalent cations enhances the interionic correlations and drastically reduces the strong interfacial repulsion, leading, in some cases, to a net attraction. We have used Monte Carlo simulations in the canonical ensemble in order to obtain a direct derivation of the free energy of two charged lamellae, including the contribution of the interionic correlations in the derivations of the interfacial energy and entropy. This interfacial free energy may be used as an effective potential describing the stability of such charged interfaces without requiring the knowledge of the counterion distribution at the solid/liquid interfaces. This procedure is validated by an independent derivation of the net interfacial pressure.