The ability to predict the thermodynamic
properties of amine species
in CO2-loaded aqueous solutions, including their deprotonation
(pK
a) and carbamate to bicarbonate reversion
(pK
c) equilibrium constants and their
corresponding standard reaction enthalpies, is of critical importance
for the design of improved carbon capture solvents. In this study,
we used isocoulombic forms of both reactions to determine these quantities
for a large set of aqueous alkanolamine solvent systems. Our hybrid
approach involves using classical molecular dynamics simulations with
the general amber force field (GAFF) and semi-empirical AM1–BCC
charges (GAFF/AM1–BCC) in the solution phase, combined with
high-level composite quantum chemical ideal-gas calculations. We first
determined a new force field (FF) for the hydronium ion (H3O+) by matching to the single experimental pK
a data point for the well-known monoethanolamine system
at 298.15 K. We then used this FF to predict the pK
a values for 76 other amines at 298.15 K and for all 77
amines at elevated temperatures. Additionally, we indirectly relate
the H3O+ hydration free energy to that of H+ and provide expressions for intrinsic hydration free energy
and enthalpy of the proton. Using the derived H3O+ FF, we predicted the pK
a values of a
diverse set of alkanolamines with an overall average absolute deviation
of less than 0.72 pK
a units. Furthermore,
the derived H3O+ FF is able to predict the protonation
enthalpy of these amines when used with the GAFF. We also predicted
the carbamate reversion constants of the primary and secondary amine
species in the data set and their corresponding standard heats of
reaction, which we compared with the scarcely available experimental
data, which are often subject to significant uncertainty. Finally,
we also described the influence of electronic and steric effects of
different molecular fragments/groups on the stabilities of the carbamates.