Distribution of proton dissociation constants for model humic and fulvic acid molecules.

The intrinsic proton binding constants of 10 model humic acid and six model fulvic acid molecules are calculated using SPARC Performs Automated Reasoning in Chemistry (SPARC). The accuracy of the SPARC calculations is examined using estimated microscopic binding constants of various small organic acids. An equimolar mixture of the appropriate hypothetical molecules is used as a representation of soil and aqueous humic acid and fulvic acid. The probability distributions of the mixture microscopic proton binding constants and the intrinsic proton binding constants in the metal speciation models WHAM V and WHAM VI (Windermere humic aqueous models) are compared. The idea is to assess the predictive value of the molecular mixture models as representations of heterogeneous natural organic matter. For aqueous humic and fulvic acids, the results are comparable to the WHAM distribution. For soil humic acid, the WHAM probability distribution is less acidic for the carboxylic sites but similar to that of the phenolic sites. Computations made using the WHAM molecular distributions and WHAM VI are comparable to titration data for Suwannee River fulvic acid. These results suggest that mixture molecular models can be used to investigate and predict the binding of metal cations to humic and fulvic acids.