Continuous Fractional Component Monte Carlo (CFCMC) and molecular dynamics (MD) simulations are performed to calculate the solubilities and self-diffusion coefficients of four light n-alkanes (methane, ethane, propane, and n-butane) in aqueous NaCl solutions as well as the thermodynamic properties of their corresponding hydrate crystals. Correction factors k_ij to the Lorentz-Berthelot combining rules for alkane groups (CH₃) and water are optimized (k_ij = 1.04) by fitting excess chemical potentials to experimental data at 1 bar and 298.15 K. Using these values of k_ij, we calculate the solubilities of the four alkanes in aqueous NaCl solutions with different molalities (0-6) mol/kg at different temperatures (278.15-308.15) K and pressures (1, 100, 200, 300) bar. The diffusion coefficients of the four alkanes in NaCl solutions (0-6) mol/kg are calculated at different temperatures (278.15-308.15) K and 1 bar and corrected for the finite-size effects. The lattice parameters of the corresponding hydrates with different guest molecules are computed using MD simulations at different temperatures (150-290) K and pressures (5-700) MPa. Isothermal compressibilities at 287.15 K and thermal expansion coefficients at 14.5 MPa for the corresponding hydrates are calculated. We present an extensive collection of thermodynamic data related to gas hydrates that contribute to a fundamental understanding of natural gas hydrate science.

@en

Knowledge of the microscopic behavior of CO₂ hydrates in oceanic sediments is crucial to evaluate the efficiency and stability of hydrate-based CO₂ sequestration in oceans. Here, systematic molecular dynamics simulations are executed to investigate the growth and dissociation of CO₂ hydrates, and the mechanical instability of CO₂ hydrate-Illite interface in the brine-urea-Illite system. Simulation results show that the CO₂ hydrate growth is jointly affected by the confined space, Illite surface properties, and presence of urea. Specifically, the interfacial H₂O and the ion layer on the Illite surface hinder the growth of CO₂ hydrate crystals toward Illite surfaces. Urea molecules can bind salt ions and increase CO₂ concentrations in the water, thus kinetically promoting CO₂ hydrate growth. The dissociation of the CO₂ hydrate is affected by Illite surface properties and the CO₂ hydrate structure. CO₂ hydrate starts from the regions where hydrate particles are minimally in contact and extends on both sides. The mechanical tension and compression of the CO₂ hydrate-Illite interface exhibit nonlinear characteristics by changing the hydrogen bonds and the CO₂ hydrate structure. The molecular insight into the microscopic behavior of CO₂ hydrates in the brine-urea-Illite system contributes to a broader understanding of hydrate-based CO₂ sequestration.

@en

Knowledge on the kinetics of gas hydrate dissociation in clay pores at static and dynamic fluid conditions is a fundamental scientific issue for improving gas production efficiency from hydrate deposits using thermal stimulation and depressurization respectively. Here, molecular dynamics simulations were used to investigate poly- and mono-crystalline methane hydrates in Na-montmorillonite clay nanopores. Simulation results show that hydrate dissociation is highly sensitive to temperature and pressure gradients, but their effects differ. Temperature changes increase thermal instability of water and gas molecules, leading to layer-by-layer dissociation from the outer surface. Under flow conditions, laminar flow predominates in nano-pores, and non-Darcy flow occurs due to clay-fluid interactions. Viscous flow disrupts hydrogen bonding at the hydrate surface, enhancing kinetic instability of water. Grain boundaries of polycrystalline hydrates are less stable compared to bulk phases and preferentially decompose, forming new dissociation fronts. This accelerates dissociation compared to monocrystalline hydrates. Fracture occurs at the grain boundaries of polycrystalline hydrate in the fluid, resulting in separate hydrate crystal grains. This fracture process further accelerates hydrate dissociation. In flow systems, methane nanobubbles form in fluid and readily transport with fluid flow. Unlike surface nanobubbles at static conditions, these liquid nanobubbles exhibit mobility. The findings of this study can contribute to a better understanding of the complex phase transition behavior of hydrate in confined environment, and provide theoretical support for improving production control technology.

@en

An in-depth review of the available experimental and molecular simulation studies of CO₂ diffusion in H₂O, which is a central property in important industrial and environmental processes, such as carbon capture and storage, enhanced oil recovery, and in the food industry is presented. The cases of both bulk and confined systems are covered. The experimental and molecular simulation data gathered are analyzed, and simple and computationally efficient correlations are devised. These correlations are applicable to conditions from 273 K and 0.1 MPa up to 473 K and 45 MPa. The available experimental data for diffusion coefficients of CO₂ in brines are also collected, and their dependency on temperature, pressure, and salinity is examined in detail. Other engineering models and correlations reported in literature are also presented. The review of the simulation studies focuses on the force field combinations, the data for diffusivities at low and high pressures, finite-size effects, and the correlations developed based on the Molecular Dynamics data. Regarding the confined systems, we review the main methods to measure and compute the diffusivity of confined CO₂ and discuss the main natural and artificial confining media (i.e., smectites, calcites, silica, MOFs, and carbon materials). Detailed discussion is provided regarding the driving force for diffusion of CO₂ and H₂O under confinement, and on the role of effects such as H₂O adsorption on hydrophilic confining media on the diffusivity of CO₂. Finally, an outlook of future research paths for advancing the field of CO₂ diffusivity in H₂O at the bulk phase and in confinement is laid out.

@en

Hydrogen is a clean-burning fuel that can be converted to other forms. of energy without generating any greenhouse gases. Currently, hydrogen is stored either by compression to high pressure (>700 bar) or cryogenic cooling to liquid form (<23 K). Therefore, it is essential to develop safe, reliable, and energy-efficient storage technology that can store hydrogen at lower pressures and temperatures. In this work, we systematically designed 2902 Mg-alkoxide-functionalized covalent-organic frameworks (COFs) and performed high-throughput (HT) computational screening for hydrogen storage applications at 111, 231, and 296 K. To accurately model the interaction between Mg-alkoxide sites and molecular hydrogen, we performed MP2 calculations to compute the hydrogen binding energy for different types of functionalized models, and the data were subsequently used to fit modified-Morse force field (FF) parameters. Using the developed FF models, we conducted HT grand canonical Monte Carlo (GCMC) simulations to compute hydrogen uptakes for both original and functionalized COFs. The generated data were subsequently used to evaluate the materials’ gravimetric and volumetric storage performance at various temperatures (111, 231, and 296 K). Finally, we developed machine learning (ML) models to predict the hydrogen storage performance of functionalized structures based on the features of the original structures. The developed model showed excellent performance with a mean absolute error (MAE) of 0.061 wt % and 0.456 g/L for predicting the gravimetric and volumetric deliverable capacities, enabling a quick evaluation of structures in a hypothetical COF database. The screening results demonstrated that the Mg-alkoxide functionalization yields greater improvements in volumetric H2 storage capacities for COFs with smaller pores compared to those with larger (mesoporous) pores.@en

H₂-CO₂ mixtures find wide-ranging applications, including their growing significance as synthetic fuels in the transportation industry, relevance in capture technologies for carbon capture and storage, occurrence in subsurface storage of hydrogen, and hydrogenation of carbon dioxide to form hydrocarbons and alcohols. Here, we focus on the thermodynamic properties of H₂-CO₂ mixtures pertinent to underground hydrogen storage in depleted gas reservoirs. Molecular dynamics simulations are used to compute mutual (Fick) diffusivities for a wide range of pressures (5 to 50 MPa), temperatures (323.15 to 423.15 K), and mixture compositions (hydrogen mole fraction from 0 to 1). At 5 MPa, the computed mutual diffusivities agree within 5% with the kinetic theory of Chapman and Enskog at 423.15 K, albeit exhibiting deviations of up to 25% between 323.15 and 373.15 K. Even at 50 MPa, kinetic theory predictions match computed diffusivities within 15% for mixtures comprising over 80% H₂ due to the ideal-gas-like behavior. In mixtures with higher concentrations of CO₂, the Moggridge correlation emerges as a dependable substitute for the kinetic theory. Specifically, when the CO₂ content reaches 50%, the Moggridge correlation achieves predictions within 10% of the computed Fick diffusivities. Phase equilibria of ternary mixtures involving CO₂-H₂-NaCl were explored using Gibbs Ensemble (GE) simulations with the Continuous Fractional Component Monte Carlo (CFCMC) technique. The computed solubilities of CO₂ and H₂ in NaCl brine increased with the fugacity of the respective component but decreased with NaCl concentration (salting out effect). While the solubility of CO₂ in NaCl brine decreased in the ternary system compared to the binary CO₂-NaCl brine system, the solubility of H₂ in NaCl brine increased less in the ternary system compared to the binary H₂-NaCl brine system. The cooperative effect of H₂-CO₂ enhances the H₂ solubility while suppressing the CO₂ solubility. The water content in the gas phase was found to be intermediate between H₂-NaCl brine and CO₂-NaCl brine systems. Our findings have implications for hydrogen storage and chemical technologies dealing with CO₂-H₂ mixtures, particularly where experimental data are lacking, emphasizing the need for reliable thermodynamic data on H₂-CO₂ mixtures.

@en

Thermodynamic factors for diffusion connect the Fick and Maxwell-Stefan diffusion coefficients used to quantify mass transfer. Activity coefficient models or equations of state can be fitted to experimental or simulation data, from which thermodynamic factors can be obtained by differentiation. The accuracy of thermodynamic factors determined using indirect routes is dictated by the specific choice of an activity coefficient model or an equation of state. The Permuted Widom’s Test Particle Insertion (PWTPI) method developed by Balaji et al. enables direct determination of thermodynamic factors in binary and multicomponent systems. For highly dense systems, for example, typical liquids, it is well known that molecular test insertion methods fail. In this article, we use the Continuous Fractional Component Monte Carlo (CFCMC) method to directly calculate thermodynamic factors by adopting the PWTPI method. The CFCMC method uses fractional molecules whose interactions with their surrounding molecules are modulated by a coupling parameter. Even in highly dense systems, the CFCMC method efficiently handles molecule insertions and removals, overcoming the limitations of the PWTPI method. We show excellent agreement between the results of the PWTPI and CFCMC methods for the calculation of thermodynamic factors in binary systems of Lennard-Jones molecules and ternary systems of Weeks-Chandler-Andersen molecules. The CFCMC method applied to calculate the thermodynamic factors of realistic molecular systems consisting of binary mixtures of carbon dioxide and hydrogen agrees well with the NIST REFPROP database. Our study highlights the effectiveness of the CFCMC method in determining thermodynamic factors for diffusion, even in densely packed systems, using relatively small numbers of molecules.

@en

We explore the impact of force field parameters and reaction equilibrium on the scaling behavior towards the critical point in reactive binary systems, focusing on NO ₂/N ₂O ₄. This system can be considered as a special single-component system since NO ₂ and N ₂O ₄ are in chemical equilibrium via the chemical reaction 2NO ₂⇌N ₂O ₄. We simplify the system by representing both components as single LJ particles, achieving excellent agreement with densities computed using molecular simulations in which all-atom force fields were used. We investigate the effect of force field parameters (ɛ and σ) on phase behavior and show that the critical exponent β remains constant, which means that intermolecular interactions do not affect the scaling to the critical point when the chemical reaction takes place. We also investigate the sensitivity of the reaction equilibrium constant and show that even small changes in isolated molecule partition functions lead to large differences in chemical equilibria. We show that the critical exponent β is different for systems with different reaction equilibrium constants, so a careful parameterization of β is needed for an accurate computation of critical temperatures of reactive mixtures. We perform a screening of reactive binary mixtures for a wide range of ideal gas reaction equilibrium constants, revealing key insights into the thermodynamic behavior and critical properties. Thereby we facilitate the efficient screening of reactive binary mixtures for various applications. Our results emphasize the importance of accurately parameterizing β and provide valuable insights into the critical scaling behavior of complex reactive systems.

@en

Focused ultrasound has experimentally been found to enhance the diffusion of nanoparticles; our aim with this work is to study this effect closer using both experiments and non-equilibrium molecular dynamics. Measurements from single particle tracking of 40 nm polystyrene nanoparticles in an agarose hydrogel with and without focused ultrasound are presented and compared with a previous experimental study using 100 nm polystyrene nanoparticles. In both cases, we observed an increase in the mean square displacement during focused ultrasound treatment. We developed a coarse-grained non-equilibrium molecular dynamics model with an implicit solvent to investigate the increase in the mean square displacement and its frequency and amplitude dependencies. This model consists of polymer fibers and two sizes of nanoparticles, and the effect of the focused ultrasound was modeled as an external oscillating force field. A comparison between the simulation and experimental results shows similar mean square displacement trends, suggesting that the particle velocity is a significant contributor to the observed ultrasound-enhanced mean square displacement. The resulting diffusion coefficients from the model are compared to the diffusion equation for a two-time continuous time random walk. The model is found to have the same frequency dependency. At lower particle velocity amplitude values, the model has a quadratic relation with the particle velocity amplitude as described by the two-time continuous time random walk derived diffusion equation, but at higher amplitudes, the model deviates, and its diffusion coefficient reaches the non-hindered diffusion coefficient. This observation suggests that at higher ultrasound intensities in hydrogels, the non-hindered diffusion coefficient can be used.

@en

Vapor-Liquid Equilibria (VLE) of hydrogen (H₂) and aqueous electrolyte (KOH and NaCl) solutions are central to numerous industrial applications such as alkaline electrolysis and underground hydrogen storage. Continuous fractional component Monte Carlo simulations are performed to compute the VLE of H₂ and aqueous electrolyte solutions at 298-423 K, 10-400 bar, 0-8 mol KOH/kg water, and 0-6 mol NaCl/kg water. The densities and activities of water in aqueous KOH and NaCl solutions are accurately modeled (within 2% deviation from experiments) using the non-polarizable Madrid-2019 Na⁺/Cl⁻ ion force fields for NaCl and the Madrid-Transport K⁺ and Delft Force Field of OH⁻ for KOH, combined with the TIP4P/2005 water force field. A free energy correction (independent of pressure, salt type, and salt molality) is applied to the computed infinite dilution excess chemical potentials of H₂ and water, resulting in accurate predictions (within 5% of experiments) for the solubilities of H₂ in water and the saturated vapor pressures of water for a temperature range of 298-363 K. The compositions of water and H₂ are computed using an iterative scheme from the liquid phase excess chemical potentials and densities, in which the gas phase fugacities are computed using the GERG-2008 equation of state. For the first time, the VLE of H₂ and aqueous KOH/NaCl systems are accurately captured with respect to experiments (i.e., for both the liquid and gas phase compositions) without compromising the liquid phase properties or performing any refitting of force fields.

@en

The use of reactive working fluids in thermodynamic cycles is currently being considered as an alternative to inert working fluids, because of the preliminarily attested higher energy-efficiency potential. The current needs to simulate their use in thermodynamic cycles, which may operate in liquid, vapour or vapour-liquid state, are an accurate real-fluid equation of state and ideal gas thermochemical properties of each molecule constituting the mixture, to calculate the equilibrium constant. To this end, the appeal to a multi-scale theoretical methodology is paramount and its definition represents the objective of the present work. This methodology is applied and validated on the system N₂O₄ ⇌ 2NO₂. Firstly, the equations solved for simultaneous two-phase and reaction equilibrium are presented. Secondly, ideal gas thermochemical properties of N₂O₄ and NO₂ are computed at atomic scale by quantum mechanics simulations. Then, to apply the selected cubic equation of state, pure-component properties of the species forming the reactive mixture (critical point coordinates and acentric factor) are required as input. However, these properties are not measurable, since NO₂ and N₂O₄ do not exist in nature as pure components. To get around this difficulty, the methodology relies on molecular Monte Carlo simulations of the pure N₂O₄ and NO₂, as well as on the reactive N₂O₄ ⇌ 2NO₂, enabling the determination of those missing pure-component properties and thus the calculation, on a macroscopic scale, of the reactive mixture properties. Finally, the comparison of calculated mixture properties with available experimental data leads to validate the accuracy of the proposed methodology.

@en

Non-polarizable force fields fail to accurately predict free energies of aqueous electrolytes without compromising the predictive ability for densities and transport properties. A new approach is presented in which (1) TIP4P/2005 water and scaled charge force fields are used to describe the interactions in the liquid phase and (2) an additional Effective Charge Surface (ECS) is used to compute free energies at zero additional computational expense. The ECS is obtained using a single temperature-independent charge scaling parameter per species. Thereby, the chemical potential of water and the free energies of hydration of various aqueous salts (e.g., NaCl and LiCl) are accurately described (deviations less than 5% from experiments), in sharp contrast to calculations where the ECS is omitted (deviations larger than 20%). This approach enables accurate predictions of free energies of aqueous electrolyte solutions using non-polarizable force fields, without compromising liquid-phase properties.@en

Hydrate dissociation is often accompanied by the formation of nanobubbles. Knowledge of the effects of nanobubbles on hydrate dissociation is essential for understanding the dynamic behavior of the hydrate phase change and improving the gas production efficiency. Here, molecular dynamics simulations were performed to study the methane hydrate dissociation kinetics with and without a pre-existing methane nanobubble. The results show that the hydrate cluster in the liquid phase dissociates layer-by-layer. This process is shown to be independent of the temperature and nanobubble presence at the simulation conditions. Hydrate dissociation does not always lead to nanobubble formation because the supersaturated methane solution can be stable for a long time. A steep methane concentration gradient was observed between the hydrate cluster surface and the methane nanobubble, which can enhance the directional migration of methane and effectively minimize the methane concentration in the liquid phase, thereby increasing the driving force for the hydrate dissociation. Our findings indicate that the presence of a nanobubble near the hydrate surface does not decrease the activation energy of hydrate dissociation, but it can increase the intrinsic decomposition rate. The average hydrate dissociation rate is linearly correlated with the mass flow rate towards the nanobubble. The mass flow rate is determined by the nanobubble size and hydrate-nanobubble distance. Our findings contribute to the fundamental understanding of the dissociation mechanism of gas hydrates in the liquid phase, which is crucial for the design and optimization of efficient gas hydrate production techniques.@en

The chemical feedback between building blocks in templated polymerization of diblock copolymers and their consecutive micellization was studied for the first time by means of coarse-grained molecular dynamics simulations. Using a stochastic polymerization model, we were able to reproduce the experimental findings on the effect of chemical feedback on the polymerization rates at low and high solution concentrations. The size and shape of micelles were computed using a newly developed software in Python conjugated with graph theory. In full agreement with the experiments, our simulations revealed that micelles formed by the templated micellization are more spherical and have a lower radius of gyration than those formed by the traditional two-step micellization method. The advantage of molecular simulation over the traditional kinetic models is that with the simulation, one studies in detail the heterogeneous polymerization in the presence of the oppositely charged template while also accounting for the incompatibility between reacted species, which significantly influences the reaction process.

@en

Due to the intermittency of renewable energy sources, alkaline water electrolyzers are typically operated at partial load compared to the nominal design value. It is well-known that gas crossover is dominant at low current densities leading to higher anodic hydrogen content and higher cathodic oxygen content in the separator tanks. High anodic hydrogen content is tantamount to loss of product hydrogen which results in an explosive atmosphere in the gas phase if the volumetric hydrogen content in oxygen exceeds 4%. We have developed a transient model of a multi-cell stack which can describe the operation of the electrolyzer with mixed electrolyte flows (anolyte and catholyte), separated flows, or a combination thereof (dynamic switching). This is a major extension of the steady-state model developed by Haug et al. (International Journal of Hydrogen Energy, 2017, 42, 15,689–15707). In sharp contrast to the steady-state model by Haug et al., the transient model can calculate the gas crossover as the operating conditions (e.g. electrolyte flow cycles) dynamically change in time. Depending on the size of the stack and the separator tanks, the model estimates different rates for impurities to build up. The transient model is validated using independent experimental results by Haug et al. and Brauns et al. (Electrochimica Acta, 2022, 404, 139,715) The results show that the dynamic model can follow experimental results for fluctuating current densities for a period of several days. We found that the dynamic response and transition time to steady state depend significantly on the geometrical volume of the gas separators with respect to the single-cell stack. For a multi-cell stack, we find that the impurities build-up faster when increasing the number of cells in the stack. This model serves as a tool for sizing and process management of the electrolyzer system and the separator tanks especially with respect to explosion safety.

@en

Data for several key thermodynamic and transport properties needed for technologies using hydrogen (H₂), such as underground H₂ storage and H₂O electrolysis are scarce or completely missing. Force field-based Molecular Dynamics (MD) and Continuous Fractional Component Monte Carlo (CFCMC) simulations are carried out in this work to cover this gap. Extensive new data sets are provided for (a) interfacial tensions of H₂ gas in contact with aqueous NaCl solutions for temperatures of (298 to 523) K, pressures of (1 to 600) bar, and molalities of (0 to 6) mol NaCl/kg H₂O, (b) self-diffusivities of infinitely diluted H₂ in aqueous NaCl solutions for temperatures of (298 to 723) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H₂O, and (c) solubilities of H₂ in aqueous NaCl solutions for temperatures of (298 to 363) K, pressures of (1 to 1000) bar, and molalities of (0 to 6) mol NaCl/kg H₂O. The force fields used are the TIP4P/2005 for H₂O, the Madrid-2019 and the Madrid-Transport for NaCl, and the Vrabec and Marx for H₂. Excellent agreement between the simulation results and available experimental data is found with average deviations lower than 10%.

@en

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₂ 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₂ isotherms and speciations are in excellent agreement with experimental data, demonstrating the accuracy and the precision of our solver. The binary absorptions of CO₂ and H₂S in 50 wt % MDEA/water solutions at 323.15 K were computed and compared with available data from the literature. The computed CO₂ isotherms showed good agreement with other modeling studies from the literature while the computed H₂S isotherms did not agree well with experimental data. The experimental equilibrium constants used as an input were not adjusted for H₂S/CO₂/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₂ pressures were significantly underestimated. We systematically investigated the limitations of computing CO₂ 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.

@en

Water transport is one of the major factors responsible for moisture damage in asphalt pavements. To study the thermodynamics and kinetics of water transport in bitumen and to uncover microscale mechanisms of moisture-induced damage, molecular dynamics simulations were performed for up to 600 ns for water–bitumen systems with realistic water contents that varied from 0 to 1.76 wt%. Hydrogen bonding interactions and clustering of water molecules at various combinations of temperature and water content were investigated, and their effects on the self-diffusion coefficient of water and bitumen properties are computed and discussed. It is shown that the saturated water concentration in bitumen is small, especially at low temperatures, and additional water molecules tend to form large water clusters via hydrogen bonding, indicating micro-phase separation of the water and bitumen phases inside the simulation box. Hydrogen bonding and water clustering play a crucial role on the magnitude of the self-diffusion coefficient of water. Physical properties of bitumen that include viscosity and cohesive energy are affected by water. The presence of large water clusters is indicative of how degradation in cohesion is observed on the microscale.

@en

Experimentally measuring the diffusivities of CO₂ and H₂S in aqueous alkanolamine solutions presents an extremely challenging task. To overcome this challenge, we performed Molecular Dynamics (MD) simulations to study the effects of temperature and N-methyldiethanolamine (MDEA) concentration on self-diffusivities of CO₂ (D_CO2) and H₂S (D_H2S) in aqueous MDEA solutions. We compute the densities and viscosities of aqueous MDEA solutions for an MDEA concentration range of 10–50 wt% and a temperature range of 288–333 K showing an excellent agreement with experimental data from literature. We compute the self-diffusivity of MDEA (D_MDEA) in aqueous MDEA solutions and our findings show that the computed values of D_MDEA are in excellent agreement with experimental and simulation results from literature. The self-diffusivities D_CO2 and D_H2S in aqueous MDEA solutions are computed for a wide range of temperatures and MDEA concentrations and our results show that both D_CO2 and D_H2S depend significantly on temperature and MDEA concentration. We also show that both CO₂ and H₂S diffuse slower in aqueous MDEA solutions than in aqueous MEA solutions. By comparing the radial distribution functions of CO₂, H₂S, water, and MDEA, we show that H₂S has stronger interactions with the surrounding molecules than CO₂, which makes H₂S diffuse slower in aqueous MDEA solutions. We also investigate the densities and viscosities of acid gas loaded aqueous MDEA solutions and self-diffusivities of the reaction products of CO₂ and H₂S with aqueous MDEA solutions. We show that the self-diffusivities of CO₂-loaded solutions significantly decrease with increasing CO₂ loading while the self-diffusivities of H₂S-loaded solutions do not change with changing H₂S loading. Our results will be helpful in the design and optimization of acid gas removal units.

@en

Electrochemical reduction of carbon dioxide (CO₂) to useful products is an emerging power-to-X concept, which aims to produce chemicals and fuels with renewable electricity instead of fossil fuels. Depending on the catalyst, a range of chemicals can be produced from CO₂ electrolysis at industrial-scale current densities, high Faraday efficiencies, and relatively low cell voltages. One of the main challenges for up-scaling the process is related to (bi)carbonate formation (carbonation), which is a consequence of performing the reaction in alkaline media to suppress the competing hydrogen evolution reaction. The parasitic reactions of CO₂ with the alkaline electrolytes result in (bi)carbonate precipitation and flooding in gas diffusion electrodes, CO₂ crossover to the anode, low carbon utilization efficiencies, electrolyte carbonation, pH-drift in time, and additional cost for CO₂ and electrolyte recycling. We present a critical review of the causes, consequences, and possible solutions for the carbonation effect in CO₂ electrolyzers. The mechanism of (bi)carbonate crossover in different cell configurations, its effect on the overall process design, and the economics of CO₂ and electrolyte recovery are presented. The aim is to provide a better understanding of the (bi)carbonate problem and guide research directions to overcome the challenges related to low-temperature CO₂ electrolysis in alkaline media.

@en