Abstract
We developed an open-source chemical reaction equilibrium solver in Python (CASpy, https://github.com/omoultosEthTuDelft/CASpy) to compute the concentration of species in any reactive liquid-phase absorption system. We derived an expression for a mole fraction-based equilibrium constant as a function of excess chemical potential, standard ideal gas chemical potential, temperature, and volume. As a case study, we computed the CO(2) absorption isotherm and speciation in a 23 wt % N-methyldiethanolamine (MDEA)/water solution at 313.15 K, and compared the results with available data from the literature. The results show that the computed CO(2) isotherms and speciations are in excellent agreement with experimental data, demonstrating the accuracy and the precision of our solver. The binary absorptions of CO(2) and H(2)S in 50 wt % MDEA/water solutions at 323.15 K were computed and compared with available data from the literature. The computed CO(2) isotherms showed good agreement with other modeling studies from the literature while the computed H(2)S isotherms did not agree well with experimental data. The experimental equilibrium constants used as an input were not adjusted for H(2)S/CO(2)/MDEA/water systems and need to be adjusted for this system. Using free energy calculations with two different force fields (GAFF and OPLS-AA) and quantum chemistry calculations, we computed the equilibrium constant (K) of the protonated MDEA dissociation reaction. Despite the good agreement of the OPLS-AA force field (ln[K] = -24.91) with the experiments (ln[K] = -23.04), the computed CO(2) pressures were significantly underestimated. We systematically investigated the limitations of computing CO(2) absorption isotherms using free energy and quantum chemistry calculations and showed that the computed values of μ(i)(ex) are very sensitive to the point charges used in the simulations, which limits the predictive power of this method.