Solvent-mediated modification of thermodynamics and kinetics of monoethanolamine regeneration reaction in amine-stripping carbon capture: Computational chemistry study.

A major limitation of amine-based post-combustion carbon capture technology is the necessity to regenerate amines at high temperatures, which dramatically increases operating costs. This paper concludes the effect of solvent choice as a possible route to modify the thermodynamics and kinetics characterizing the involved amine regeneration reactions and discusses whether these modifications can be economically beneficial. We report experimentally benchmarked computational chemistry calculations of monoethanolamine regeneration reactions employing aqueous and non-aqueous solvents with a wide range of dielectric constants. Unlike previous studies, our improved computational chemistry framework could accurately reproduce the right experimental activation energy of zwitterion formation. From the thermodynamics and kinetics of the predicted reactions, the use of non-aqueous solvents with small dielectric constants led to reductions in regeneration Gibbs free energies, activation barriers, and enthalpy changes. This can reduce energy consumption and give an opportunity to run desorption columns at relatively lower temperatures, thus offering the possibility of relying on low-grade waste heat as an energy input.

Preferential heating of aqueous amine solutions using infrared radiation at selected vibrational frequencies: A molecular dynamics study.

A recent CO2 capture experiment suggests that microwaves might be beneficial for regeneration of aqueous amine solutions due to both thermal and nonthermal effects [S. J. McGurk et al., Appl. Energy 192, 126 (2017)]. We use classical molecular dynamics to simulate heating of aqueous amine solutions using electromagnetic radiation with different frequencies in both microwave and infrared regions. The infrared frequencies were selected based on the partial vibrational density of states of water and amine. Unlike the microwave case, we found that preferential heating of water or amine can be achieved using their relevant vibrational frequencies in the infrared region, suggesting that microwave heating is not an optimal choice for an efficient amine regeneration reported in a recent carbon capture experiment. Interestingly, only flexible water models augmented with an anharmonic O-H bond stretching potential were able to quantitatively predict the expected differential heating profiles of systems involving water.

Molecular dynamics simulation of microwave heating of liquid monoethanolamine (MEA): An evaluation of existing force fields.

We present a complete classical molecular dynamics (MD) study of the dielectric heating of liquid monoethanolamine (MEA) at microwave (MW) frequencies ranging from 1.0 to 10.0 GHz. The detailed dielectric properties predicted by a series of existing empirical force fields of MEA were carefully compared to experimental results. We find that all the evaluated force fields were unable to accurately predict experimental static dielectric constant, frequency-dependent dielectric spectra, and MW heating profiles of liquid MEA, although GROMOS-aa (all-atom GROningen molecular simulation) is the most accurate of those tested. With an isotropic scaling of partial atomic charges, the modified GROMOS-aa and OPLS-aa (all-atom optimized potentials for liquid simulations) force fields could accurately reproduce the experimental static dielectric constant and frequency-dependent dielectric spectra, but they failed to predict MW heating rates directly from MD heating simulations. Thus, the recently presented approach [F. J. Salas et al., J. Chem. Theory Comput. 11, 683 (2015); A. P. de la Luz et al., ibid. 11, 2792 (2015)] to tune existing force fields is not an ideal approach to produce force fields suitable for accurate dielectric heating studies.

Classical molecular dynamics simulation of microwave heating of liquids: The case of water.

We perform a complete classical molecular dynamics study of the dielectric heating of water in the microwave (MW) region. MW frequencies ranging from 1.0 to 15.0 GHz are used together with a series of well-known empirical force fields. We show that the ability of an empirical force field to correctly predict the dielectric response of liquids to MW radiation should be evaluated on the basis of a joint comparison of the predicted and experimental static dielectric constant, frequency-dependent dielectric spectra, and heating profiles. We argue that this is essential when multicomponent liquids are studied. We find that both the three-site OPC3 and four-site TIP4P-ϵ empirical force fields of water are equally superior for reproducing dielectric properties at a range of MW frequencies. Despite its poor prediction of the static dielectric constant, the well-known SPCE force field can be used to accurately describe dielectric heating of water at low MW frequencies.

Survey of classical density functionals for modelling hydrogen physisorption at 77K.

This work surveys techniques based on classical density functionals for modeling the quantum dispersion of physisorbed hydrogen at 77K. Two such techniques are examined in detail. The first is based on the "open ring approximation" (ORA) of Broukhno et al., and it is compared with a technique based on the semiclassical approximation of Feynman and Hibbs (FH). For both techniques, a standard classical density functional is used to model hydrogen molecule-hydrogen molecule (i.e., excess) interactions. The three-dimensional (3D) quantum harmonic oscillator (QHO) system and a model of molecular hydrogen adsorption into a graphitic slit pore at 77K are used as benchmarks. Density functional results are compared with path-integral Monte Carlo simulations and with exact solutions for the 3D QHO system. It is found that neither of the density functional treatments are entirely satisfactory. However, for hydrogen physisorption studies at 77K the ORA based technique is generally superior to the FH based technique due to a fortunate cancellation of errors in the density functionals used. But, if more accurate excess functionals are used, the FH technique would be superior.

Self-referential Monte Carlo method for calculating the free energy of crystalline solids.

A self-referential Monte Carlo method is described for calculating the free energy of crystalline solids. All Monte Carlo methods for the free energy of classical crystalline solids calculate the free-energy difference between a state whose free energy can be calculated relatively easily and the state of interest. Previously published methods employ either a simple model crystal, such as the Einstein crystal, or a fluid as the reference state. The self-referential method employs a radically different reference state; it is the crystalline solid of interest but with a different number of unit cells. So it calculates the free-energy difference between two crystals, differing only in their size. The aim of this work is to demonstrate this approach by application to some simple systems, namely, the face centered cubic hard sphere and Lennard-Jones crystals. However, it can potentially be applied to arbitrary crystals in both bulk and confined environments, and ultimately it could also be very efficient.

Gas adsorption in active carbons and the slit-pore model 1: Pure gas adsorption.

We describe procedures based on the polydisperse independent ideal slit-pore model, Monte Carlo simulation and density functional theory (a 'slab-DFT') for predicting gas adsorption and adsorption heats in active carbons. A novel feature of this work is the calibration of gas-surface interactions to a high surface area carbon, rather than to a low surface area carbon as in all previous work. Our models are used to predict the adsorption of carbon dioxide, methane, nitrogen, and hydrogen up to 50 bar in several active carbons at a range of near-ambient temperatures based on an analysis of a single 293 K carbon dioxide adsorption isotherm. The results demonstrate that these models are useful for relatively simple gases at near-critical or supercritical temperatures.

Gas adsorption in active carbons and the slit-pore model 2: Mixture adsorption prediction with DFT and IAST.

We use a fast density functional theory (a "slab-DFT") and the polydisperse independent ideal slit-pore model to predict gas mixture adsorption in active carbons. The DFT is parametrized by fitting to pure gas isotherms generated by Monte Carlo simulation of adsorption in model graphitic slit-pores. Accurate gas molecular models are used in our Monte Carlo simulations with gas-surface interactions calibrated to a high surface area carbon, rather than a low surface area carbon as in all previous work of this type, as described in part 1 of this work. We predict the adsorption of binary mixtures of carbon dioxide, methane, and nitrogen on two active carbons up to about 30 bar at near-ambient temperatures. We compare two sets of results; one set obtained using only the pure carbon dioxide adsorption isotherm as input to our pore characterization process, and the other obtained using both pure gas isotherms as input. We also compare these results with ideal adsorbed solution theory (IAST). We find that our methods are at least as accurate as IAST for these relatively simple gas mixtures and have the advantage of much greater versatility. We expect similar results for other active carbons and further performance gains for less ideal mixtures.

Weighted density-functional theory for simple fluids: prewetting of a Lennard-Jones fluid.

The prewetting of a Lennard-Jones fluid is studied using weighted density-functional theory. The intrinsic Helmholtz free-energy functional is separated into repulsive and attractive contributions. An accurate functional for hard spheres is used for the repulsive functional and a weighted density-functional method is used for the attractive part. The results for this theory are compared against mean-field density-functional theory, the theory of Velasco and Tarazona [E. Velasco and P. Tarazona, J. Chem. Phys. 91, 7916 (1989)] and grand canonical ensemble simulation results. The results demonstrate that the weighted density functional for attractive forces may offer a significant increase in accuracy over the other theories. The density-functional and simulation results also indicate that a previous estimate of the wetting temperature for a model of the interaction of argon with solid carbon dioxide, obtained from simulations [J. E. Finn and P. A. Monson, Phys. Rev. A, 39, 6402 (1989)], is incorrect. The weighted density-functional method indicates that triple-point prewetting is observed for this model potential.

Weighted density functional theory for simple fluids: supercritical adsorption of a Lennard-Jones fluid in an ideal slit pore.

The adsorption of a Lennard-Jones fluid in an ideal slit pore is studied using weighted density functional theory. The intrinsic Helmholtz free-energy functional is separated into repulsive and attractive contributions. Rosenfeld's accurate fundamental measure functional is employed for the repulsive functional while another weighted density functional method is employed for the attractive functional. This other method requires an accurate equation of state for the bulk fluid and an accurate pair-direct correlation function for a uniform fluid, determined analytically or numerically. The results for this theory are compared against mean-field density functional theory and grand canonical ensemble simulation results, modeling the adsorption of ethane in a graphite slit. The results indicate that the weighted density functional method applied to the attractive functional can offer a significant increase in accuracy over the mean-field theory.