Molecular simulations using the Grand Canonical Monte Carlo (GCMC) method have been performed in order to obtain physical insights on how the interaction between ibuprofen (IBP) and activated carbons (ACs) in aqueous mixtures affects IBP removal from water by ACs. A nanoporous carbon model based on units of polyaromatic molecules with different number of rings, defects and polar-oxygenated sites is described. Individual effects of factors such as porous features and chemical heterogeneities in the adsorbents are investigated and quantified. Results are in good agreement with experimental adsorption data, highlightening the ability of GCMC simulation to describe the macroscopic adsorption performance in drug removal applications, while also providing additional insights into the IBP/water adsorption mechanism. The simulation results allow finding the optimal type of activated carbon material for separating this pollutant in water treatment. (C) 2017 Elsevier Inc. All rights reserved.