Binary mixtures of supercritical carbon dioxide with methanol. A molecular dynamics simulation study

Molecular Dynamics Simulation on Thickening and Solubility Properties of Novel Thickener in Supercritical Carbon Dioxide

Supercritical CO2 has wide application in enhancing oil recovery, but the low viscosity of liquid CO2 can lead to issues such as poor proppant-carrying ability and high filtration loss. Therefore, the addition of thickening agents to CO2 is vital. Hydrocarbon polymers, as a class of green and sustainable materials, hold tremendous potential for acting as thickeners in supercritical CO2 systems, and PVAc is one of the best-performing hydrocarbon thickeners. To further improve the viscosity enhancement and solubility of PVAc, here we designed a novel polymer structure, PVAO, by introducing CO2-affine functional groups to PVAc. Molecular dynamics simulations were adopted to analyze viscosity and relevant solubility parameters systematically. We found that PVAO exhibits superior performance, with a viscosity enhancement of 1.5 times that of PVAc in supercritical CO2. While in the meantime, PVAO maintains better solubility characteristics than PVAc. Our findings offer insights for the future design of other high-performance polymers.

Thermodynamics of mixing methanol with supercritical CO2 as seen from computer simulations and thermodynamic integration

The changes in extensive thermodynamic quantities, such as volume, energy, Helmholtz free energy and entropy, occurring upon mixing liquid methanol with supercritical CO₂, are calculated using Monte Carlo simulations and thermodynamic integration for all eight combinations of four methanol and two CO₂ potential models in the entire composition range at 313 K. The obtained results are also compared with experimental data whenever possible. The transition of the system from liquid to a supercritical state is found to occur at this temperature around a CO₂ mole fraction value of 0.95 with all model combinations considered. This liquid to supercritical transition is always accompanied by positive Helmholtz free energy of mixing values and, consequently, by the non-miscibility of the two components. Furthermore, both this non-miscibility around the liquid to supercritical transition and also the miscibility of the two components below this transition, in the liquid regime, are found to be primarily of the energetic rather than entropic origin; the entropy of mixing turns out to be very close to zero, and around the liquid to supercritical transition even its qualitative behaviour is strongly model dependent. Finally, it is found that the methanol expansion coefficient is not sensitive to the details of the potential models, and it is always in excellent agreement with the experimental data. On the other hand, both the volume and the energy of mixing depend strongly on the molar volume of neat CO₂ in the model being used, and in this respect the TraPPE model of CO₂ [J. J. Potoff and J. I. Siepmann, AIChE J., 2001, 47, 1676] performs considerably better than that of Zhang and Duan [Z. Zhang and Z. Duan, J. Chem. Phys., 2005, 122, 214507].

Electron Donor-Acceptor Nature of the Ethanol-CO2 Dimer

Supercritical CO₂ is an appealing nontoxic, environmentally friendly solvent for the industrial extraction of many classes of compounds, from caffeine to natural product drug precursors to petrochemical impurities. Apolar in isolation, the ability of supercritical CO₂ to dissolve polar species has been empirically shown to be greatly enhanced by the addition of a small molar percentage of a polar cosolvent, often ethanol. Computational work predicts that the isolated ethanol-CO₂ complex can exist either in an electron-donor configuration or through a hydrogen-bonding one; yet, neither has been previously experimentally observed. Here, we demonstrate by rotational spectroscopy that the isolated, gas-phase ethanol-CO₂ dimer is an electron donor-acceptor complex.

Hydrogen bond, electron donor-acceptor dimer, and residence dynamics in supercritical CO(2)-ethanol mixtures and the effect of hydrogen bonding on single reorientational and translational dynamics: A molecular dynamics simulation study

The hydrogen bonding and dynamics in a supercritical mixture of carbon dioxide with ethanol as a cosolvent (X(ethanol) approximately 0.1) were investigated using molecular dynamics simulation techniques. The results obtained reveal that the hydrogen bonds formed between ethanol molecules are significantly more in comparison with those between ethanol-CO(2) molecules and also exhibit much larger lifetimes. Furthermore, the residence dynamics in the solvation shells of ethanol and CO(2) have been calculated, revealing much larger residence times for ethanol molecules in the ethanol solvation shell. These results support strongly the ethanol aggregation effects and the slow local environment reorganization inside the ethanol solvation shell, reported in a previous publication of the authors [Skarmoutsos et al., J. Chem. Phys. 126, 224503 (2007)]. The formation of electron donor-acceptor dimers between the ethanol and CO(2) molecules has been also investigated and the calculated lifetimes of these complexes have been found to be similar to those corresponding to ethanol-CO(2) hydrogen bonds, exhibiting a slightly higher intermittent lifetime. However, the average number of these dimers is larger than the number of ethanol-CO(2) hydrogen bonds in the system. Finally, the effect of the hydrogen bonds formed between the individual ethanol molecules on their reorientational and translational dynamics has been carefully explored showing that the characteristic hydrogen bonding microstructure obtained exhibits sufficiently strong influence upon the behavior of them.

Solvation and solvatochromism in CO2-expanded liquids. 1. Simulations of the solvent systems CO2 + cyclohexane, acetonitrile, and methanol

Molecular dynamics simulations of CO(2)-expanded cyclohexane, acetonitrile, and methanol are reported at various compositions along the experimental bubble-point curve at 298 K. Simulated properties include energies, local compositions, viscosities, diffusion coefficients, and dielectric constants and relaxation times. On the basis of the limited comparisons to experimental data currently available, the results indicate that simple intermolecular potential models previously developed for simulating the pure components provide reasonable representations of the energetics and dynamics of these gas-expanded liquids.

Hydrogen Bonding of Methanol in Supercritical CO2: Comparison between 1H NMR Spectroscopic Data and Molecular Simulation Results

Molecular dynamics simulation results on hydrogen bonding in mixtures of methanol with CO2 at supercritical, liquid-like conditions are compared to 1H NMR spectroscopic data that have recently become available. The molecular models are parametrized using vapor-liquid equilibrium data only, which they reliably describe. A new molecular model for methanol of Lennard-Jones plus point charge type is presented. This molecular methanol model is investigated in terms of its capability to yield hydrogen-bonding statistics. Simple assumptions are made regarding the assignment of NMR chemical shifts to the different types of hydrogen-bonded species. Only two state-independent parameters are fitted to the large NMR data set on the basis of hydrogen-bonding statistics from molecular simulations. Excellent agreement between the molecular simulation results and the NMR data is found. This shows that the molecular models of the simple type studied here cannot only describe thermodynamic properties but also structural effects of hydrogen bonding in solutions.

Investigation of the local composition enhancement and related dynamics in supercritical CO2-cosolvent mixtures via computer simulation: The case of ethanol in CO2

The supercritical mixture ethanol-carbon dioxide (EtOH-CO2) with mole fraction of ethanol X(EtOH) congruent with 0.1 was investigated at 348 K, by employing the molecular dynamics simulation technique in the canonical ensemble. The local intermolecular structure of the fluid was studied in terms of the calculated appropriate pair radial distribution functions. The estimated average local coordination numbers and mole fractions around the species in the mixture reveal the existence of local composition enhancement of ethanol around the ethanol molecules. This finding indicates the nonideal mixing behavior of the mixture due to the existence of aggregation between the ethanol molecules. Furthermore, the local environment redistribution dynamics have been explored by analyzing the time correlation functions (TCFs) of the total local coordination number (solvent, cosolvent) around the cosolvent molecules in appropriate parts. The analysis of these total TCFs in the auto-(solvent-solvent, cosolvent-cosolvent) and cross-(solvent-cosolvent, cosolvent-solvent) TCFs has shown that the time dependent redistribution process of the first solvation shell of ethanol is mainly determined by the redistribution of the CO2 solvent molecules. These results might be explained on the basis of the CO2-CO2 and EtOH-CO2 intermolecular forces, which are sufficiently weaker in comparison to the EtOH-EtOH hydrogen bonding interactions, creating in this way a significantly faster redistribution of the CO2 molecules in comparison with EtOH. Finally, the self-diffusion coefficients and the single reorientational dynamics of both the cosolvent and solvent species in the mixture have been predicted and discussed in relationship with the local environment around the species, which in the case of the EtOH molecules seem to be strongly affected.

Solvation and Solvatochromism in CO2-Expanded Liquids. 2. Experiment−Simulation Comparisons of Preferential Solvation in Three Prototypical Mixtures

Electronic absorption and emission spectra of 10-bis(phenylethynyl)anthracene (PEA) and coumarin 153 (C153) are measured as functions of composition along the bubble-point curve at 25 degrees C in CO2-expanded cyclohexane (c-C6H12), acetonitrile (CH3CN), and methanol (CH3OH). The nonlinear dependence of the spectral frequencies on composition suggests substantial preferential solvation of both solutes by the liquid components of these mixtures. Estimates of enrichment factors (local mole fraction of a component divided by its bulk value) based on the assumption that spectral shifts are linearly related to local composition are quite large (approximately 10) in the cases of the C153/CH3CN + CO2 and C153/CH3OH + CO2 systems at high xCO2. Computer simulations of anthracene, the chromophore of PEA, and C153 in these three CO2-expanded liquids are used to clarify the relationship between local composition and spectral shift. A semiempirical model consisting of additive electrostatic and dispersive interactions is able to capture the main features observed experimentally in all six solute/solvent combinations. The simulations show that the commonly used assumption of a linear relation between spectral shifts and local compositions grossly exaggerates the extent of preferential solvation in these mixtures. The collective nature of electrostatic solvation and the composition dependence of the solute's coordination number are shown to be responsible for the breakdown of this assumption.

Local density augmentation and dynamic properties of hydrogen-and non-hydrogen-bonded supercritical fluids: A molecular dynamics study

The local density inhomogeneities in neat supercritical fluids were investigated via canonical molecular dynamics simulations. The selected systems under investigation were the polar and hydrogen-bonded fluid methanol as well as the quadrupolar non-hydrogen-bonded carbon dioxide one. Effective local densities, local density augmentation, and enhancement factors were calculated at state points along an isotherm close to the critical temperature of each system (T(r)=1.03). The results obtained reveal strong influence of the polarity and hydrogen bonding upon the intensity of the local density augmentation. It is found that this effect is sufficiently larger in the case of the polar and associated methanol in comparison to those predicted for carbon dioxide. For both fluids the local density augmentation values are maximized in the bulk density region near 0.7rho(c), a result that is in agreement with experiment. In addition, the local density dynamics of each fluid were investigated in terms of the appropriate time correlation functions. The behavior of these functions reveals that the bulk density dependence of the local density reorganization times is very sensitive to the specific intermolecular interactions and to the size of the local region. Also, the estimated local density reorganization time as a function of bulk density of each fluid was further analyzed and successfully related to two different time-scale relaxation mechanisms. Finally, the results obtained indicate a possible relationship between the single-molecule reorientational dynamics and the local density reorganization ones.

Molecular Dynamics Simulation of the Cybotactic Region in Gas-Expanded Methanol−Carbon Dioxide and Acetone−Carbon Dioxide Mixtures

Local solvation and transport effects in gas-expanded liquids (GXLs) are reported based on molecular simulation. GXLs were found to exhibit local density enhancements similar to those seen in supercritical fluids, although less dramatic. This approach was used as an alternative to a multiphase atomistic model for these mixtures by utilizing experimental results to describe the necessary fixed conditions for a locally (quasi-) stable molecular dynamics model of the (single) GXL phase. The local anisotropic pair correlation function, orientational correlation functions, and diffusion rates are reported for two systems: CO2-expanded methanol and CO2-expanded acetone at 298 K and pressures up to 6 MPa.

Local Density Inhomogeneities and Dynamics in Supercritical Water: A Molecular Dynamics Simulation Approach

Molecular dynamics atomistic simulations in the canonical ensemble (NVT-MD) have been used to investigate the "Local Density Inhomogeneities and their Dynamics" in pure supercritical water. The simulations were carried out along a near-critical isotherm (Tr = T/Tc = 1.03) and for a wide range of densities below and above the critical one (0.2 rho(c) - 2.0 rho(c)). The results obtained reveal the existence of significant local density augmentation effects, which are found to be sufficiently larger in comparison to those reported for nonassociated fluids. The time evolution of the local density distribution around each molecule was studied in terms of the appropriate time correlation functions C(Delta)rhol(t). It is found that the shape of these functions changes significantly by increasing the density of the fluid. Finally, the local density reorganization times for the first and second coordination shell derived from these correlations exhibit a decreasing behavior by increasing the density of the system, signifying the density effect upon the dynamics of the local environment around each molecule.

Pressure and Temperature Dependence of the Hydrogen Bonding in Supercritical Ethanol: A Computer Simulation Study

As a step toward deeper insight on the "hydrogen bonding" in supercritical ethanol (scEtOH), we carried out NVT molecular dynamics simulations of the fluid over a wide range of temperatures and pressures. The fluid was studied at SC conditions for which thermodynamic and spectroscopic (NMR, infrared, Raman, dielectric) data are available. The various site-site pair distribution functions (pdf's) were calculated, and their temperature and pressure dependence was obtained. It was found that over the thermodynamic conditions investigated here, scEtOH remains highly structured. Moreover, the characteristic behavior of the first peaks in H-H, O-O, and H-O pdf's reveals that hydrogen bonds still exist in scEtOH. The analysis focuses also on the reorientational dynamics of the bond unit vectors O-H, C-O, and of the permanent dipole moment of the molecules as well as the total dipole moment of the sample. The corresponding Legendre time correlation functions were discussed in connection to the "hydrogen bonding" in the fluid and in the context of experimental results. Specifically, the behavior of the O-H dynamics exhibits the well-known associative nature of the molecules in the system. A further analysis of the hydrogen bonds was carried out, and the degree of aggregation (average number of H-bonds per molecule) was obtained and compared with results from NMR chemical shift studies. Also the estimated monomer and free O-H groups in the fluid were compared with results from IR and Raman vibrational spectroscopy. The percentage analysis fi of the liquid and scEtOH molecules, with i = 0, 1, 2, 3, ... hydrogen bonds per molecule, has been obtained. The results show the existence of small, linear-chain oligomers formed mainly by two molecules, whereas the number of the three body oligomers, and specifically that of four body oligomers in the sample, is relatively small.

Excited-State Proton Transfer in Gas-Expanded Liquids: The Roles of Pressure and Composition in Supercritical CO2/Methanol Mixtures

Excited-state proton transfer from 5,8-dicyano-2-naphthol to methanol takes place in CO2/methanol mixtures, in the pressure and temperature ranges of supercritical CO2. The efficiency of the proton-transfer step decreases with the pressure. This is assigned to the perturbation of the methanol clusters solvating the naphthol.

Dielectric Constant and Density Dependence of the Structure of Supercritical Carbon Dioxide Using a New Modified Empirical Potential Model: A Monte Carlo Simulation Study

Two modified versions of the Elementary Physical Model (EPM) [J. Phys. Chem. 1995, 99, 12021] for supercritical carbon dioxide have been proposed in this work and their validities are affirmed by computing the thermodynamic properties and dielectric constant up to 910 kg/m3 with use of canonical ensemble Monte Carlo simulation. Simulations performed for 500 molecules with the EPM2-M model reproduce the experimental data accurately at all thermodynamic states. The structural analyses demonstrate that the aggregation is strong at low density while the coordination number is large at high density. In addition, a detailed study on the radial and angular correlation functions reveals that the T-shaped geometry is dominate while a variety of other structures still appear in the first coordination shell. Furthermore, the angular correlation functions show that the probability of a molecule being oriented toward the convex side of another molecule is equal to that pointing toward the concave side since the molecular nonlinearity of carbon dioxide is only marginal. As the distance between two molecules increases, the preferred orientations disappear quickly and all the results are in good agreement with the prior ab initio calculation [J. Chem. Phys. 2004, 120, 9694].

Binary phase behavior and aggregation of dilute methanol in supercritical carbon dioxide: A Monte Carlo simulation study

Configurational-bias Monte Carlo simulations in the Gibbs and isobaric-isothermal ensembles using the transferable potentials for phase equilibria force field were carried out to investigate the thermophysical properties of mixtures containing supercritical carbon dioxide and methanol. The binary vapor-liquid coexistence curves were calculated at 333.15 and 353.15 K and are in excellent agreement with experimental measurements. The self-association of methanol in supercritical carbon dioxide was investigated over a range of temperatures and pressures near the mixture critical point. The temperature dependence of the equilibrium constants for the formation of hydrogen-bonded aggregates (from dimer to heptamer) allowed for the determination of the enthalpy of hydrogen bonding, DeltaHHB, in supercritical carbon dioxide with values for DeltaHHB of about 15 kJ mol(-1) falling within the range of previously proposed values. No strong pressure dependence was observed for the formation of aggregates. Apparently the decrease of the entropic penalty and of the enthalpic benefit upon increasing pressure or solvent density mostly cancel each other's effect on aggregate formation.