Atomistic molecular dynamics simulations of H2O diffusivity in liquid and supercritical CO2

Analysis of the Distribution and Influencing Factors of Diffusion Coefficient Model Parameters Based on Molecular Dynamics Simulations

The establishment of mathematical models to predict the diffusion coefficients of gas and liquid systems have important theoretical significance and practical value. In this work, based on the previously proposed diffusion coefficient model D_LV, the distribution and influencing factors of the model parameters characteristic length (L) and diffusion velocity (V) were further investigated using molecular dynamics simulations. The statistical analysis of L and V for 10 gas systems and 10 liquid systems was presented in the paper. New distribution functions were established to describe the probability distributions of molecular motion L and V. The mean values of the correlation coefficients were 0.98 and 0.99, respectively. Meanwhile, the effects of molecular molar mass and system temperature on the molecular diffusion coefficients were discussed. The results show that the effect of molecular molar mass on the diffusion coefficient mainly affects the molecular motion L, and the effect of system temperature on the diffusion coefficient mainly affects V. For the gas system, the average relative deviation of D_LV and D_MSD is 10.73% and that of D_LV and experimental value is 12.63%; for the solution system, the average relative deviation of D_LV and D_MSD is 12.93% and that of D_LV and experimental value is 18.86%, which indicates the accuracy of the new model results. The new model reveals the potential mechanism of molecular motion and provides a theoretical basis for further study of the diffusion process.

Strongly diluted dimethyl-imidazolium chloride–alcohol solutions: solvents are structurally different but dynamic heterogeneities are similar

Based on the analysis of dynamic properties of ionic liquid solutions, the descriptions of diffusion mechanisms are built for dimethylimidazolium chloride (dmim⁺/Cl⁻)–alcohol solute systems and the influence of the monohydric alcohols' molecular structure on their diffusion mechanisms in dmim⁺/Cl⁻–alcohol at T = 400 K by molecular dynamics simulations are studied. From the analysis of radial distribution functions, MSDs, velocity autocorrelation function, and autocorrelation functions of dispersion we found that the motion of all components in IL dmim⁺/Cl⁻–alcohol (ethanol, propanol) systems at T = 400 K occurs in the sub-diffuse regime and that the dynamics of the dmim⁺/Cl⁻–alcohol (ethanol, propanol) systems is heterogeneous. The increase of the alkyl chain length of the alcohol molecule does not affect the motion of the ionic liquid components; instead, it increases the characteristic times describing the model representation of alcohol molecule diffusion at short and medium times, without affecting diffusion mechanisms.

The increase of the alkyl chain length of the alcohol molecule does not affect the motion of the ionic liquid components.

Diffusivity and hydrophobic hydration of hydrocarbons in supercritical CO2 and aqueous brine

CO₂ injection (EOR and sequestration technique) creates the amalgamation of hydrocarbons, CO₂, and aqueous brine in the subsurface. In this study, molecular dynamics (MD) simulations were used to investigate the diffusivity of hydrocarbon molecules in a realistic scenario of supercritical CO₂ (SC-CO₂) injection in the subsurface over a wide range of pressures (50 < P < 300 bar) and aqueous brine concentrations (0, 2, and 5% brine). To overcome existing challenges in traditional diffusivity calculation approaches, we took advantage of fundamental molecular-based methods, along with further verification of results by previously published experimental data. In this regard, computational methods and MD simulations were employed to compute diffusion coefficients of hydrocarbons (benzene and pentane). It was found that the presence of water and salt affects the thermodynamic properties of molecules where the intermolecular interactions caused the hydrophobic hydration of hydrocarbons coupled with ionic hydration due to hydrogen bond and ion-dipole interactions. Based on these results, it is demonstrated that the formation of water clusters in the SC-CO₂ solvent is a major contributor to the diffusion of hydrophobic molecules. The outcome at different pressure conditions showed that hydrocarbons always would diffuse less in the presence of water. The slopes of linearly fitted MSD of benzene and pentane infinitely diluted in SC-CO₂ is around 13 to 20 times larger than the slope with water molecules (4 wt%). When pressure increases (100–300 bar), the diffusion coefficients (D) of benzene and pentane decreases (around 1.2 × 10⁻⁹ to 0.4 × 10⁻⁹ and 2 × 10⁻⁹ to 1 × 10⁻⁹ m² s⁻¹, respectively). Furthermore, brine concentration generally plays a negative role in reducing the diffusivity of hydrocarbons due to the formation of water clusters as a result of hydrophobic and ionic hydration. Under the SC-CO₂ rich (injection) system in the shale reservoir, the diffusion of hydrocarbon is correlated to the efficiency of hydrocarbon flow/recovery. Ultimately, this study will guide us to better understand the phenomena that would occur in nanopores of shale that undergo EOR or are becoming a target of CO₂ sequestration.

CO₂ injection (EOR and sequestration technique) creates the amalgamation of hydrocarbons, CO₂, and aqueous brine in the subsurface.

Finite-Size Effects of Binary Mutual Diffusion Coefficients from Molecular Dynamics

Molecular dynamics simulations were performed for the prediction of the finite-size effects of Maxwell-Stefan diffusion coefficients of molecular mixtures and a wide variety of binary Lennard-Jones systems. A strong dependency of computed diffusivities on the system size was observed. Computed diffusivities were found to increase with the number of molecules. We propose a correction for the extrapolation of Maxwell-Stefan diffusion coefficients to the thermodynamic limit, based on the study by Yeh and Hummer ( J. Phys. Chem. B , 2004 , 108 , 15873 - 15879 ). The proposed correction is a function of the viscosity of the system, the size of the simulation box, and the thermodynamic factor, which is a measure for the nonideality of the mixture. Verification is carried out for more than 200 distinct binary Lennard-Jones systems, as well as 9 binary systems of methanol, water, ethanol, acetone, methylamine, and carbon tetrachloride. Significant deviations between finite-size Maxwell-Stefan diffusivities and the corresponding diffusivities at the thermodynamic limit were found for mixtures close to demixing. In these cases, the finite-size correction can be even larger than the simulated (finite-size) Maxwell-Stefan diffusivity. Our results show that considering these finite-size effects is crucial and that the suggested correction allows for reliable computations.

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 CO₂, CH₄, and N₂ 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.

Molecular Modeling of Thermodynamic and Transport Properties for CO2 and Aqueous Brines

Molecular simulation techniques using classical force-fields occupy the space between ab initio quantum mechanical methods and phenomenological correlations. In particular, Monte Carlo and molecular dynamics algorithms can be used to provide quantitative predictions of thermodynamic and transport properties of fluids relevant for geologic carbon sequestration at conditions for which experimental data are uncertain or not available. These methods can cover time and length scales far exceeding those of quantum chemical methods, while maintaining transferability and predictive power lacking from phenomenological correlations. The accuracy of predictions depends sensitively on the quality of the molecular models used. Many existing fixed-point-charge models for water and aqueous mixtures fail to represent accurately these fluid properties, especially when descriptions covering broad ranges of thermodynamic conditions are needed. Recent work on development of accurate models for water, CO₂, and dissolved salts, as well as their mixtures, is summarized in this Account. Polarizable models that can respond to the different dielectric environments in aqueous versus nonaqueous phases are necessary for predictions of properties over extended ranges of temperatures and pressures. Phase compositions and densities, activity coefficients of the dissolved salts, interfacial tensions, viscosities and diffusivities can be obtained in near-quantitative agreement to available experimental data, using relatively modest computational resources. In some cases, for example, for the composition of the CO₂-rich phase in coexistence with an aqueous phase, recent results from molecular simulations have helped discriminate among conflicting experimental data sets. The sensitivity of properties on the quality of the intermolecular interaction model varies significantly. Properties such as the phase compositions or electrolyte activity coefficients are much more sensitive than phase densities, viscosities, or component diffusivities. Strong confinement effects on physical properties in nanoscale media can also be directly obtained from molecular simulations. Future work on molecular modeling for CO₂ and aqueous brines is likely to be focused on more systematic generation of interaction models by utilizing quantum chemical as well as direct experimental measurements. New ion models need to be developed for use with the current generation of polarizable water models, including ion-ion interactions that will allow for accurate description of dense, mixed brines. Methods will need to be devised that go beyond the use of effective potentials for incorporation of quantum effects known to be important for water, and reactive force fields developed that can handle bond creation and breaking in systems with carbonate and silicate minerals. Another area of potential future work is the integration of molecular simulation methods in multiscale models for the chemical reactions leading to mineral dissolution and flow within the porous media in underground formations.

Physical Absorption of Green House Gases in Amines: The Influence of Functionality, Structure, and Cross-Interactions

Monte Carlo simulations were performed in the isothermal-isobaric ensemble (NPT) to calculate the Henry constants of methane (CH₄), nitrous oxide (N₂O), and carbon dioxide (CO₂) in pure H₂O, amines, and alkanolamines using the classical Lorentz-Berthelot combining rules (L-B). The Henry constants of N₂O and CO₂ in water are highly overestimated and motivated us to propose a new set of unlike interactions. Contrarily, the Henry constant of N₂O in MEA is underestimated by around 40%, and again, a new reoptimized cross unlike parameter is able to reproduce the constant to within 10%. An analysis is given of the relationship between the physical absorption of these gases and the chemical structure or functionality of 12 molecules including amines and alkanolamines using the anisotropic united atom intermolecular potential (AUA4). Finally, the solubility of N₂O in an aqueous solution of monoethanolamine (MEA) at 30% (wt) was also studied. A Henry constant within 7% of the experimental value was found by using the reoptimized parameters along with L-B to account for the MEA + H₂O unlike interactions. This very good agreement without additional adjustments for the MEA + H₂O system may be attributed to the good excess properties predictions found in previous works for the binary mixture (MEA + H₂O). However, further work, including additional alkanolamines in aqueous solutions at several concentrations, is required to verify this particular point.