Abstract Hypothesis The interactions of organic molecules with mineral surfaces are influenced by several factors such as adsorbate speciation, surface atomic and electronic structure, and environmental conditions. When coupled with thermodynamic techniques, energetics from atomistic modeling can provide a molecular-level picture of which factors determine reactivity. This is paramount for evaluating the chemical processes which control the fate of these species in the environment. Experiments Inner-sphere adsorption of oxalate and pyrocatechol on (001), (110), and (012) α-Fe2O3 surfaces was modeled using Density Functional Theory (DFT). Unique bidentate binding modes were sampled along each facet to study how different adsorbate and surface factors govern site preference. Adsorption energetics were then calculated using a DFT + thermodynamics approach which combines DFT energies with tabulated data and Nernst-based corrective terms to incorporate different experimental parameters. Findings Instead of a universal trend, each facet displays a unique factor that dominates site preference based on either strain (001), functional groups (110), or topography (012). Adsorption energies predict favorable inner-sphere adsorption for both molecules but opposite energetic trends with varying pH. Additionally, vibrational analysis was conducted for each system and compared to experimental IR data. The work presented here provides an effective, computational methodology to study numerous adsorption processes occurring at the surface-aqueous interface.