The adsorption equilibrium constants of monovalent and divalent cations to material surfaces in aqueous media are central to many technological, natural, and geochemical processes. Cation adsorption-desorption is often proposed to occur in concert with proton transfer on hydroxyl-covered mineral surfaces, but to date this cooperative effect has been inferred indirectly. This work applies density functional theory-based molecular dynamics simulations of explicit liquid water/mineral interfaces to calculate metal ion desorption free energies. Monodentate adsorption of Na, Mg, and Cu on partially deprotonated silica surfaces are considered. Na is predicted to be unbound, while Cu exhibits binding free energies to surface SiO groups that are larger than those of Mg. The predicted trends agree with competitive adsorption measurements on fumed silica surfaces. As desorption proceeds, Cu dissociates one of the HO molecules in its first solvation shell, turning into Cu(OH)(HO), while Mg remains Mg(HO). The protonation state of the SiO group at the initial binding site does not vary monotonically with cation desorption.