Adsorption of organic molecules on carbon surfaces: Experimental data and molecular dynamics simulation considering multiple protonation states

• Adsorption model of organics with multiple protonation and zwitterionic forms. • Loadings of maleic acid and phenylalanine over a wide pH and concentration range. • Prediction of the model parameters using molecular dynamics (MD). • Model parameters fitted via genetic algorithm compared to MD model parameters. • Explanation of pH impact on adsorption and vice versa. Owing to their high specific surface and low production cost, carbon materials are among the most important adsorption materials. Novel usages, for instance in pharmaceutical applications, challenge existing methods because charged and strongly polar substances need to be adsorbed. Here, we systematically investigate the highly complex adsorption equilibria of organic molecules having multiple protonation states as a function of pH. The adsorption behavior depends on intermolecular interactions within the solution (dissociation equilibria) and between adsorbed molecules on the carbon surface (electrostatic forces). For the model substances maleic acid and phenylalanine, we demonstrate that a custom-made genetic algorithm is able to extract up to nine parameters of a multispecies isotherm from experimental data covering a broad pH-range. The parameters, including adsorption affinities, interaction energies, and maximum loadings were also predicted by molecular dynamics simulations. Both approaches obtained a good qualitative and mostly also quantitative description of the adsorption behavior within a pH-range of 2–12. By combining the determined isotherms with mass balances, the final concentrations and pH-shifts of batch adsorption experiments can be predicted. The developed modeling tools can be easily adapted to other types of pH-dependent, multispecies adsorbates and therefore will help to optimize adsorption-based processes in different fields.