A Quantitative Analysis of Electrochemical CO2 Reduction on Copper in Organic Amide and Nitrile-Based Electrolytes.

Aqueous electrolytes used in CO2 electroreduction typically have a CO2 solubility of around 34 mM under ambient conditions, contributing to mass transfer limitations in the system. Non-aqueous electrolytes exhibit higher CO2 solubility (by 5-8-fold) and also provide possibilities to suppress the undesired hydrogen evolution reaction (HER). On the other hand, a proton donor is needed to produce many of the products commonly obtained with aqueous electrolytes. This work investigates the electrochemical CO2 reduction performance of copper in non-aqueous electrolytes based on dimethylformamide (DMF), n-methyl-2-pyrrolidone (NMP), and acetonitrile (ACN). The main objective is to analyze whether non-aqueous electrolytes are a viable alternative to aqueous electrolytes for hydrocarbon production. Additionally, the effects of aqueous/non-aqueous anolytes, membrane, and the selection of a potential window on the electrochemical CO2 reduction performance are addressed in this study. Experiments with pure DMF and NMP mainly produced oxalate with a faradaic efficiency (FE) reaching >80%; however, pure ACN mainly produced hydrogen and formate due to the presence of more residual water in the system. Addition of 5% (v/v) water to the non-aqueous electrolytes resulted in increased HER and formate production with negligible hydrocarbon production. Hence, we conclude that aqueous electrolytes remain a better choice for the production of hydrocarbons and alcohols on a copper electrode, while organic electrolytes based on DMF and NMP can be used to obtain a high selectivity toward oxalate and formate.

Isobaric Vapor-Liquid Equilibrium Data for Tetrahydrofuran + Acetic Acid and Tetrahydrofuran + Trichloroethylene Mixtures.

Vapor-liquid equilibrium (VLE) data for the binary systems tetrahydrofuran (THF) + acetic acid (AA) and THF + trichloroethylene (TCE) were measured under isobaric conditions using an ebulliometer. The boiling temperatures for the systems (THF + AA/THF + TCE) are reported for 13/15 compositions and five/six different pressures ranging from 50.2/60.0 to 101.1/101.3 kPa, respectively. The THF + AA system shows simple phase behavior with no azeotrope formation. The THF + TCE system does not exhibit azeotrope formation but seems to have a pinch point close to the pure end of TCE. The nonrandom two-liquid (NRTL) and universal quasichemical (UNIQUAC) activity coefficient models were used to accurately fit the binary (PTx) data. Both models were able to fit the binary VLE data satisfactorily. However, the NRTL model was found to be slightly better than UNIQUAC model in fitting the VLE data for both systems. The results can be used for designing liquid-liquid extraction and distillation processes involving mixtures of THF, AA, and TCE.

Solubility of CO2 in Aqueous Formic Acid Solutions and the Effect of NaCl Addition: A Molecular Simulation Study.

There is a growing interest in the development of routes to produce formic acid from CO2, such as the electrochemical reduction of CO2 to formic acid. The solubility of CO2 in the electrolyte influences the production rate of formic acid. Here, the dependence of the CO2 solubility in aqueous HCOOH solutions with electrolytes on the composition and the NaCl concentration was studied by Continuous Fractional Component Monte Carlo simulations at 298.15 K and 1 bar. The chemical potentials of CO2, H2O, and HCOOH were obtained directly from single simulations, enabling the calculation of Henry coefficients and subsequently considering salting in or salting out effects. As the force fields for HCOOH and H2O may not be compatible due to the presence of strong hydrogen bonds, the Gibbs-Duhem integration test was used to test this compatibility. The combination of the OPLS/AA force field with a new set of parameters, in combination with the SPC/E force field for water, was selected. It was found that the solubility of CO2 decreases with increasing NaCl concentration in the solution and increases with the increase of HCOOH concentration. This continues up to a certain concentration of HCOOH in the solution, after which the CO2 solubility is high and the NaCl concentration has no significant effect.

Electrochemical Reduction of CO2 to Oxalic Acid: Experiments, Process Modeling, and Economics.

We performed H-cell and flow cell experiments to study the electrochemical reduction of CO2 to oxalic acid (OA) on a lead (Pb) cathode in various nonaqueous solvents. The effects of anolyte, catholyte, supporting electrolyte, temperature, water content, and cathode potential on the Faraday efficiency (FE), current density (CD), and product concentration were investigated. We show that a high FE for OA can be achieved (up to 90%) at a cathode potential of -2.5 V vs Ag/AgCl but at relatively low CDs (10-20 mA/cm2). The FE of OA decreases significantly with increasing water content of the catholyte, which causes byproduct formation (e.g., formate, glycolic acid, and glyoxylic acid). A process design and techno-economic evaluation of the electrochemical conversion of CO2 to OA is presented. The results show that the electrochemical route for OA production can compete with the fossil-fuel based route for the base case scenario (CD of 100 mA/cm2, OA FE of 80%, cell voltage of 4 V, electrolyzer CAPEX of $20000/m2, electricity price of $30/MWh, and OA price of $1000/ton). A sensitivity analysis shows that the market price of OA has a huge influence on the economics. A market price of at least $700/ton is required to have a positive net present value and a payback time of less than 10 years. The performance and economics of the process can be further improved by increasing the CD and FE of OA by using gas diffusion electrodes and eliminating water from the cathode, lowering the cell voltage by increasing the conductivity of the electrolyte solutions, and developing better OA separation methods.

Surface Coverage as an Important Parameter for Predicting Selectivity Trends in Electrochemical CO2 Reduction.

The electrochemical CO2 reduction reaction (CO2RR) is important for a sustainable future. Key insights into the reaction pathways have been obtained by density functional theory (DFT) analysis, but so far, DFT has been unable to give an overall understanding of selectivity trends without important caveats. We show that an unconsidered parameter in DFT models of electrocatalysts-the surface coverage of reacting species-is crucial for understanding the CO2RR selectivities for different surfaces. Surface coverage is a parameter that must be assumed in most DFT studies of CO2RR electrocatalysts, but so far, only the coverage of nonreacting adsorbates has been treated. Explicitly treating the surface coverage of reacting adsorbates allows for an investigation that can more closely mimic operating conditions. Furthermore, and of more immediate importance, the use of surface coverage-dependent adsorption energies allows for the extraction of ratios of adsorption energies of CO2RR intermediates (COOHads and HCOOads) that are shown to be predictive of selectivity and are not susceptible to systematic errors. This approach allows for categorization of the selectivity of several monometallic catalysts (Pt, Pd, Au, Ag, Zn, Cu, Rh, W, Pb, Sn, In, Cd, and Tl), even problematic ones such as Ag or Zn, and does so by only considering the adsorption energies of known intermediates. The selectivity of the further reduction of COOHads can now be explained by a preference for Tafel or Heyrovsky reactions, recontextualizing the nature of selectivity of some catalysts. In summary, this work resolves differences between DFT and experimental studies of the CO2RR and underlines the importance of surface coverage.

Electro-osmotic Drag and Thermodynamic Properties of Water in Hydrated Nafion Membranes from Molecular Dynamics.

One of the important parameters in water management of proton exchange membranes is the electro-osmotic drag (EOD) coefficient of water. The value of the EOD coefficient is difficult to justify, and available literature data on this for Nafion membranes show scattering from in experiments and simulations. Here, we use a classical all-atom model to compute the EOD coefficient and thermodynamic properties of water from molecular dynamics simulations for temperatures between 330 and 420 K, and for different water contents between λ = 5 and λ = 20. λ is the ratio between the moles of water molecules to the moles of sulfonic acid sites. This classical model does not capture the Grotthuss mechanism; however, it is shown that it can predict the EOD coefficient within the range of experimental values for λ = 5 where the vehicular mechanism dominates proton transfer. For λ > 5, the Grotthuss mechanism becomes dominant. To obtain the EOD coefficient, average velocities of water and ions are computed by imposing different electric fields to the system. Our results show that the velocities of water and hydronium scale linearly with the electric field, resulting in a constant ratio of ca. 0.4 within the error bars. We find that the EOD coefficient of water linearly increases from 2 at λ = 5 to 8 at λ = 20 and the results are not sensitive to temperature. The EOD coefficient at λ = 5 is within the range of experimental values, confirming that the model can capture the vehicular transport of protons well. At λ = 20, due to the absence of proton hopping in the model, the EOD coefficient is overestimated by a factor of 3 compared to experimental values. To analyze the interactions between water and Nafion, the partial molar enthalpies and partial molar volumes of water are computed from molecular dynamics simulations. At different water uptakes, multiple linear regression is used on raw simulation data within a narrow composition range of water inside the Nafion membrane. The partial molar volumes and partial molar excess enthalpies of water asymptotically approach the molar volumes and molar excess enthalpies of pure water for water uptakes above 5. This confirms the model can capture the bulklike behavior of water in the Nafion at high water uptakes.

Solubilities and Transport Properties of CO2, Oxalic Acid, and Formic Acid in Mixed Solvents Composed of Deep Eutectic Solvents, Methanol, and Propylene Carbonate.

Recently, deep eutectic solvents (DES) have been considered as possible electrolytes for the electrochemical reduction of CO2 to value-added products such as formic and oxalic acids. The applicability of pure DES as electrolytes is hindered by high viscosities. Mixtures of DES with organic solvents can be a promising way of designing superior electrolytes by exploiting the advantages of each solvent type. In this study, densities, viscosities, diffusivities, and ionic conductivities of mixed solvents comprising DES (i.e., reline and ethaline), methanol, and propylene carbonate were computed using molecular simulations. To provide a quantitative assessment of the affinity and mass transport of CO2 and oxalic and formic acids in the mixed solvents, the solubilities and self-diffusivities of these solutes were also computed. Our results show that the addition of DES to the organic solvents enhances the solubilities of oxalic and formic acids, while the solubility of CO2 in the ethaline-containing mixtures are in the same order of magnitude with the respective pure organic components. A monotonic increase in the densities and viscosities of the mixed solvents is observed as the mole fraction of DES in the mixture increases, with the exception of the density of ethaline-propylene carbonate which shows the opposite behavior due to the high viscosity of the pure organic component. The self-diffusivities of all species in the mixtures significantly decrease as the mole fraction of DES approaches unity. Similarly, the self-diffusivities of the dissolved CO2 and the oxalic and formic acids also decrease by at least 1 order of magnitude as the composition of the mixture shifts from the pure organic component to pure DES. The computed ionic conductivities of all mixed solvents show a maximum value for mole fractions of DES in the range from 0.2 to 0.6 and decrease as more DES is added to the mixtures. Since for most mixtures studied here no prior experimental measurements exist, our findings can serve as a first data set based on which further investigation of DES-containing electrolyte solutions can be performed for the electrochemical reduction of CO2 to useful chemicals.

Electroreduction of CO2/CO to C2 Products: Process Modeling, Downstream Separation, System Integration, and Economic Analysis.

Direct electrochemical reduction of CO2 to C2 products such as ethylene is more efficient in alkaline media, but it suffers from parasitic loss of reactants due to (bi)carbonate formation. A two-step process where the CO2 is first electrochemically reduced to CO and subsequently converted to desired C2 products has the potential to overcome the limitations posed by direct CO2 electroreduction. In this study, we investigated the technical and economic feasibility of the direct and indirect CO2 conversion routes to C2 products. For the indirect route, CO2 to CO conversion in a high temperature solid oxide electrolysis cell (SOEC) or a low temperature electrolyzer has been considered. The product distribution, conversion, selectivities, current densities, and cell potentials are different for both CO2 conversion routes, which affects the downstream processing and the economics. A detailed process design and techno-economic analysis of both CO2 conversion pathways are presented, which includes CO2 capture, CO2 (and CO) conversion, CO2 (and CO) recycling, and product separation. Our economic analysis shows that both conversion routes are not profitable under the base case scenario, but the economics can be improved significantly by reducing the cell voltage, the capital cost of the electrolyzers, and the electricity price. For both routes, a cell voltage of 2.5 V, a capital cost of $10,000/m2, and an electricity price of <$20/MWh will yield a positive net present value and payback times of less than 15 years. Overall, the high temperature (SOEC-based) two-step conversion process has a greater potential for scale-up than the direct electrochemical conversion route. Strategies for integrating the electrochemical CO2/CO conversion process into the existing gas and oil infrastructure are outlined. Current barriers for industrialization of CO2 electrolyzers and possible solutions are discussed as well.

Effect of Water Content on Thermodynamic Properties of Compressed Hydrogen.

Force field-based molecular simulations were used to calculate thermal expansivities, heat capacities, and Joule-Thomson coefficients of binary (standard) hydrogen-water mixtures for temperatures between 366.15 and 423.15 K and pressures between 50 and 1000 bar. The mole fraction of water in saturated hydrogen-water mixtures in the gas phase ranges from 0.004 to 0.138. The same properties were calculated for pure hydrogen at 323.15 K and pressures between 100 and 1000 bar. Simulations were performed using the TIP3P and a modified TIP4P force field for water and the Marx, Vrabec, Cracknell, Buch, and Hirschfelder force fields for hydrogen. The vapor-liquid equilibria of hydrogen-water mixtures were calculated along the melting line of ice Ih, corresponding to temperatures between 264.21 and 272.4 K, using the TIP3P force field for water and the Marx force field for hydrogen. In this temperature range, the solubilities and the chemical potentials of hydrogen and water were obtained. Based on the computed solubility data of hydrogen in water, the freezing-point depression of water was computed ranging from 264.21 to 272.4 K. The modified TIP4P and Marx force fields were used to improve the solubility calculations of hydrogen-water mixtures reported in our previous study [Rahbari A.;J. Chem. Eng. Data2019, 64, 4103-4115] for temperatures between 323 and 423 K and pressures ranging from 100 to 1000 bar. The chemical potentials of ice Ih were calculated as a function of pressure between 100 and 1000 bar, along the melting line for temperatures between 264.21 and 272.4 K, using the IAPWS equation of state for ice Ih. We show that at low pressures, the presence of water has a large effect on the thermodynamic properties of compressed hydrogen. Our conclusions may have consequences for the energetics of a hydrogen refueling station using electrochemical hydrogen compressors.

Liquid-Liquid Extraction of Formic Acid with 2-Methyltetrahydrofuran: Experiments, Process Modeling, and Economics.

Formic acid (FA) is an interesting hydrogen (H2) and carbon monoxide (CO) carrier that can be produced by the electrochemical reduction of carbon dioxide (CO2) using renewable energy. The separation of FA from water is challenging due to the strong (cross)association of the components and the presence of a high boiling azeotrope. For the separation of dilute FA solutions, liquid-liquid extraction is preferred over conventional distillation because distilling large amounts of water is very energy-intensive. In this study, we use 2-methyltetrahydrofuran (2-MTHF) to extract FA from the CO2 electrolysis process, which typically contains <20 wt % of FA. Vapor-liquid equilibrium (VLE) data of the binary system 2-MTHF-FA and liquid-liquid equilibrium (LLE) data of the ternary system 2-MTHF-FA-water are obtained. Continuous extraction and distillation experiments are performed to test the extraction power and recovery of 2-MTHF from the extract. The VLE and LLE data are used to design a hybrid extraction and distillation process to produce a commercial grade product (85 wt % of FA). A detailed economic analysis of this hybrid extraction-distillation process is presented and compared with the existing FA separation methods. It is shown that 2-MTHF is a cost-effective solvent for FA extraction from dilute streams (<20 wt % FA).

Solubility of Carbon Dioxide, Hydrogen Sulfide, Methane, and Nitrogen in Monoethylene Glycol; Experiments and Molecular Simulation.

Knowledge on the solubility of gases, especially carbon dioxide (CO2), in monoethylene glycol (MEG) is relevant for a number of industrial applications such as separation processes and gas hydrate prevention. In this study, the solubility of CO2 in MEG was measured experimentally at temperatures of 333.15, 353.15, and 373.15 K. Experimental data were used to validate Monte Carlo (MC) simulations. Continuous fractional component MC simulations in the osmotic ensemble were performed to compute the solubility of CO2 in MEG at the same temperatures and at pressures up to 10 bar. MC simulations were also used to study the solubility of methane (CH4), hydrogen sulfide (H2S), and nitrogen (N2) in MEG at 373.15 K. Solubilities from experiments and simulations are in good agreement at low pressures, but deviations were observed at high pressures. Henry coefficients were also computed using MC simulations and compared to experimental values. The order of solubilities of the gases in MEG at 373.15 K was computed as H2S > CO2 > CH4 > N2. Force field modifications may be required to improve the prediction of solubilities of gases in MEG at high pressures and low temperatures.

Artificial intelligence and thermodynamics help solving arson cases.

In arson cases, evidence such as DNA or fingerprints is often destroyed. One of the most important evidence modalities left is relating fire accelerants to a suspect. When gasoline is used as accelerant, the aim is to find a strong indication that a gasoline sample from a fire scene is related to a sample of a suspect. Gasoline samples from a fire scene are weathered, which prohibits a straightforward comparison. We combine machine learning, thermodynamic modeling, and quantum mechanics to predict the composition of unweathered gasoline samples starting from weathered ones. Our approach predicts the initial (unweathered) composition of the sixty main components in a weathered gasoline sample, with error bars of ca. 4% when weathered up to 80% w/w. This shows that machine learning is a valuable tool for predicting the initial composition of a weathered gasoline, and thereby relating samples to suspects.

Brick-CFCMC: Open Source Software for Monte Carlo Simulations of Phase and Reaction Equilibria Using the Continuous Fractional Component Method.

We present a new molecular simulation code, Brick-CFCMC, for performing Monte Carlo simulations using state-of-the-art simulation techniques. The Continuous Fractional Component (CFC) method is implemented for simulations in the NVT/NPT ensembles, the Gibbs Ensemble, the Grand-Canonical Ensemble, and the Reaction Ensemble. Molecule transfers are facilitated by the use of fractional molecules which significantly improve the efficiency of the simulations. With the CFC method, one can obtain phase equilibria and properties such as chemical potentials and partial molar enthalpies/volumes directly from a single simulation. It is possible to combine trial moves from different ensembles. This enables simulations of phase equilibria in a system where also a chemical reaction takes place. We demonstrate the applicability of our software by investigating the esterification of methanol with acetic acid in a two-phase system.

High Pressure Electrochemical Reduction of CO2 to Formic Acid/Formate: A Comparison between Bipolar Membranes and Cation Exchange Membranes.

A high pressure semicontinuous batch electrolyzer is used to convert CO2 to formic acid/formate on a tin-based cathode using bipolar membranes (BPMs) and cation exchange membranes (CEMs). The effects of CO2 pressure up to 50 bar, electrolyte concentration, flow rate, cell potential, and the two types of membranes on the current density (CD) and Faraday efficiency (FE) for formic acid/formate are investigated. Increasing the CO2 pressure yields a high FE up to 90% at a cell potential of 3.5 V and a CD of ∼30 mA/cm2. The FE decreases significantly at higher cell potentials and current densities, and lower pressures. Up to 2 wt % formate was produced at a cell potential of 4 V, a CD of ∼100 mA/cm2, and a FE of 65%. The advantages and disadvantages of using BPMs and CEMs in electrochemical cells for CO2 conversion to formic acid/formate are discussed.

OCTP: A Tool for On-the-Fly Calculation of Transport Properties of Fluids with the Order- n Algorithm in LAMMPS.

We present a new plugin for LAMMPS for on-the-fly computation of transport properties (OCTP) in equilibrium molecular dynamics. OCTP computes the self- and Maxwell-Stefan diffusivities, bulk and shear viscosities, and thermal conductivities of pure fluids and mixtures in a single simulation. OCTP is the first implementation in LAMMPS that uses the Einstein relations combined with the order- n algorithm for the efficient sampling of dynamic variables. OCTP has low computational requirements and is easy to use because it follows the native input file format of LAMMPS. A tool for calculating the radial distribution function (RDF) of the fluid beyond the cutoff radius, while taking into account the system size effects, is also part of the new plugin. The RDFs computed from OCTP are needed to obtain the thermodynamic factor, which relates Maxwell-Stefan and Fick diffusivities. To demonstrate the efficiency of the new plugin, the transport properties of an equimolar mixture of water-methanol were computed at 298 K and 1 bar.

Combined Steam Reforming of Methane and Formic Acid To Produce Syngas with an Adjustable H2:CO Ratio.

Syngas is an important intermediate in the chemical process industry. It is used for the production of hydrocarbons, acetic acid, oxo-alcohols, and other chemicals. Depending on the target product and stoichiometry of the reaction, an optimum (molar) ratio between hydrogen and carbon monoxide (H2:CO) in the syngas is required. Different technologies are available to control the H2:CO molar ratio in the syngas. The combination of steam reforming of methane (SRM) and the water-gas shift (WGS) reaction is the most established approach for syngas production. In this work, to adjust the H2:CO ratio, we have considered formic acid (FA) as a source for both hydrogen and carbon monoxide. Using thermochemical equilibrium calculations, we show that the syngas composition can be controlled by cofeeding formic acid into the SRM process. The H2:CO molar ratio can be adjusted to a value between one and three by adjusting the concentration of FA in the reaction feed. At steam reforming conditions, typically above 900 K, FA can decompose to water and carbon monoxide and/or to hydrogen and carbon dioxide. Our results show that cofeeding FA into the SRM process can adjust the H2:CO molar ratio in a single step. This can potentially be an alternative to the WGS process.

Absorption Refrigeration Cycles with Ammonia-Ionic Liquid Working Pairs Studied by Molecular Simulation.

For absorption refrigeration, it has been shown that ionic liquids have the potential to replace conventional working pairs. Due to the huge number of possibilities, conducting lab experiments to find the optimal ionic liquid is infeasible. Here, we provide a proof-of-principle study of an alternative computational approach. The required thermodynamic properties, i.e., solubility, heat capacity, and heat of absorption, are determined via molecular simulations. These properties are used in a model of the absorption refrigeration cycle to estimate the circulation ratio and the coefficient of performance. We selected two ionic liquids as absorbents: [emim][Tf2N], and [emim][SCN]. As refrigerant NH3 was chosen due to its favorable operating range. The results are compared to the traditional approach in which parameters of a thermodynamic model are fitted to reproduce experimental data. The work shows that simulations can be used to predict the required thermodynamic properties to estimate the performance of absorption refrigeration cycles. However, high-quality force fields are required to accurately predict the cycle performance.

Thermodynamic and Transport Properties of Crown-Ethers: Force Field Development and Molecular Simulations.

Crown-ethers have recently been used to assemble porous liquids (PLs), which are liquids with permanent porosity formed by mixing bulky solvent molecules (e.g., 15-crown-5 ether) with solvent-inaccessible organic cages. PLs and crown-ethers belong to a novel class of materials, which can potentially be used for gas separation and storage, but their performance for this purpose needs to be assessed thoroughly. Here, we use molecular simulations to study the gas separation performance of crown-ethers as the solvent of porous liquids. The TraPPE force field for linear ether molecules has been adjusted by fitting a new set of torsional potentials to accurately describe cyclic crown-ether molecules. Molecular dynamics (MD) simulations have been used to compute densities, shear viscosities, and self-diffusion coefficients of 12-crown-4, 15-crown-5, and 18-crown-6 ethers. In addition, Monte Carlo (MC) simulations have been used to compute the solubility of the gases CO2, CH4, and N2 in 12-crown-4 and 15-crown-5 ether. The computed properties are compared with available experimental data of crown-ethers and their linear counterparts, i.e., polyethylene glycol dimethyl ethers.

Efficient Application of Continuous Fractional Component Monte Carlo in the Reaction Ensemble.

A new formulation of the Reaction Ensemble Monte Carlo technique (RxMC) combined with the Continuous Fractional Component Monte Carlo method is presented. This method is denoted by serial Rx/CFC. The key ingredient is that fractional molecules of either reactants or reaction products are present and that chemical reactions always involve fractional molecules. Serial Rx/CFC has the following advantages compared to other approaches: (1) One directly obtains chemical potentials of all reactants and reaction products. Obtained chemical potentials can be used directly as an independent check to ensure that chemical equilibrium is achieved. (2) Independent biasing is applied to the fractional molecules of reactants and reaction products. Therefore, the efficiency of the algorithm is significantly increased, compared to the other approaches. (3) Changes in the maximum scaling parameter of intermolecular interactions can be chosen differently for reactants and reaction products. (4) The number of fractional molecules is reduced. As a proof of principle, our method is tested for Lennard-Jones systems at various pressures and for various chemical reactions. Excellent agreement was found both for average densities and equilibrium mixture compositions computed using serial Rx/CFC, RxMC/CFCMC previously introduced by Rosch and Maginn (Journal of Chemical Theory and Computation, 2011, 7, 269-279), and the conventional RxMC approach. The serial Rx/CFC approach is also tested for the reaction of ammonia synthesis at various temperatures and pressures. Excellent agreement was found between results obtained from serial Rx/CFC, experimental results from literature, and thermodynamic modeling using the Peng-Robinson equation of state. The efficiency of reaction trial moves is improved by a factor of 2 to 3 (depending on the system) compared to the RxMC/CFCMC formulation by Rosch and Maginn.

Simulating the reactions of CO2 in aqueous monoethanolamine solution by reaction ensemble Monte Carlo using the continuous fractional component method.

Molecular simulations were used to compute the equilibrium concentrations of the different species in CO2/monoethanolamine solutions for different CO2 loadings. Simulations were performed in the Reaction Ensemble using the continuous fractional component Monte Carlo method at temperatures of 293, 333, and 353 K. The resulting computed equilibrium concentrations are in excellent agreement with experimental data. The effect of different reaction pathways was investigated. For a complete understanding of the equilibrium speciation, it is essential to take all elementary reactions into account because considering only the overall reaction of CO2 with MEA is insufficient. The effects of electrostatics and intermolecular van der Waals interactions were also studied, clearly showing that solvation of reactants and products is essential for the reaction. The Reaction Ensemble Monte Carlo using the continuous fractional component method opens the possibility of investigating the effects of the solvent on CO2 chemisorption by eliminating the need to study different reaction pathways and concentrate only on the thermodynamics of the system.