One major source of complexity in the implementation of nanoparticles in aqueous electrolytes arises from the strong influence that biological environments has on their physicochemical properties. A key parameter for understanding the molecular mechanisms governing the physicochemical properties of nanoparticles is the formation of the surface charge density. In this article, we present an efficient and accurate approach that combines a recently introduced classical solvation density functional theory for spherical electrical double layers with a surface complexation model to account for ion-ion correlation and excluded volume effects on the surface titration of spherical nanoparticles. We apply the proposed computational approach to account for the charge-regulated mechanisms on the surface chemistry of spherical silica (SiO2) nanoparticles. We analyze theeffects of the nanoparticle size, as well as pH level and electrolyte concentration of the aqueous solution on the nanoparticle’s surface charge density and Zeta potential. We validate our predictions for 580 Å and 200 Å nanoparticles immersed in acid, neutral and alkaline mono-valent aqueous electrolyte solutions against experimental data. Our results on mono-valent electrolyte show that the excluded volume and ion-ioncorrelations contribute significantly to the surface charge density and Zeta potential of the nanoparticle at high electrolyte concentration and pH levels, where the solvent crowding effects and electrostatic screening have shown a profound influence on the protonation/deprotonation reactions at the liquid/solute interface. The success of this approach in describing physicochemical properties of silica nanoparticles supports its broader application to study other spherical metal oxide nanoparticles.