The adsorption of Si on goethitc (alpha-FeOOH) has been studied in batch experiments that cover the natural range of Si concentrations as found in the environment. The results have been interpreted and quantified with the charge distribution (CD) and multi-site surface complexation (MUSIC) model in combination with an extended Stern (ES) layer model option. This new double layer approach (ES) accounts for ordering of interfacial water molecules leading to stepwise changes in the location of electrolyte ions near the surface [T. Hiemstra, W.H. Van Riemsdijk, J. Colloid Interface Sci. 301 (2006) 1]. The Si adsorption on goethite peaks at a pH of similar to 9 and decreases at lower and higher pH values. Thermodynamically, the pH-dependency of silicic acid adsorption is related to the value of the proton co-adsorption and can also be linked to the Si charge distribution in the interface as is discussed. Based on published EXAFS data, the adsorption of Si on goethite was modeled as the formation of a bidentate surface complex. The ionic charge distribution (CD) of this complex can be calculated from the geometry of this surface complex, optimized with molecular orbital/density functional theory (MO/DFT), and combined with a correction for water dipole orientation. The resulting CD value has been applied successfully in the description of the adsorption data. The use of a theoretical CD value has the practical advantage of a reduction of the number of adjustable parameters with a factor 2. To describe the adsorption at a high Si loading, formation of a Si polymer, e.g. a tetramer, is proposed. Such a species is only contributing to the overall adsorption at solution concentrations above about 10(-4) M, where super saturation with respect to quartz exists. The adsorbed silica polymer hydrolyzes at high pH. The reactive ligand of the polymer is quite acid (log K approximate to 6.5-7.1), which is typically for the SiO-1 surface groups of polymerized Si, like amorphous SiO2(s), and the SiO-1 ligand of the aqueous dimer Si2O(OH)(5)O-1(aq). The applied model correctly predicts the change of particles charge and the shift in IEP due to proton release upon Si adsorption. (c) 2007 Elsevier Inc. All rights reserved.