A simple theoretical approach is developed that can be used to predict the preference of ion adduct formation (with alkali Li+, Na+, K+ and alkaline earth Ca2+, Mg2+ metals) in electrospray ionization mass spectrometry (ESI-MS) of succinic acid, associated with several protonation/deprotonation equilibria. The applied strategy consists of using a vacuum environment as well as both implicit and explicit solvation of reactive sites and density functional theory as the method of choice. These distinct levels of theory mimic the smooth transition between the condensed environment and free ion in the gas phase. Good correlation between the Gibbs free energies for protonation/adduct formation processes with peak observation in the obtained mass spectra provide insight into the physical basis behind adduct preference and selectivity. This signifies the relationship between microscopic interactions, ionization efficiency, and types of ions that reach the detector.