Metal cations are ubiquitous components in biological environments and play an important role in regulating cellular functioning and membrane properties. By applying metadynamics simulations, we have performed systematic free-energy calculations of Na + , K + , Ca 2+ , and Mg 2+ bound to phospholipid membrane surfaces for the first time. The free-energy landscapes unveil specific binding behaviors of metal cations on phospholipid membranes. Na + and K + are more likely to stay in the aqueous solution, and can easily bind to a few lipid oxygens by overcoming low free-energy barriers. Ca 2+ is most stable when bound to four lipid oxygens of the membranes, rather than being * To whom correspondence should be addressed † Technical University of Catalonia-Barcelona Tech ‡ Boston University ¶ University of Cambridge 1 hydrated in the aqueous solution. Mg 2+ is tightly hydrated, and can hardly lose a hydration water and bind directly to the membranes. When bound to the membranes, the cations' most favorable total coordination numbers with water and lipid oxygens are the same as their corresponding hydration numbers in aqueous solution, indicating a competition between ion binding to water and lipids. The binding specificity of metal cations on membranes is then highly correlated with the hydration free-energy and the size of the hydration shell.