CO2/离子液体体系热力学性质的分子动力学模拟

超临界CO2和离子液体(ILs)是两种绿色溶剂. 离子液体可以溶解超临界CO2, 而超临界CO2不能溶解离子液体. 由此设计构成的CO2/IL二元系统, 同时具备了超临界CO2和离子液体的许多优点: 既可以降低离子液体的粘度, 还便于相分离, 是新型的耦合绿色溶剂. 其物理化学性质对于设计反应、分离等过程非常重要. 因此, 本文以CO2/IL二元系统为研究对象, 通过选择合适的分子力场和系综, 运用分子动力学(MD)模拟方法研究了CO2/[bmim][PF6]、CO2/[bmim][NO3]等体系的热力学性质. 结果表明, CO2对ILs膨胀度的影响非常小, 当CO2摩尔分数为0.5时, ILs膨胀仅为15%. CO2/ILs的扩散系数远小于CO2膨胀甲醇、乙醇溶液的扩散系数. 随着CO2含量的增加, ILs的扩散系数提高, 粘度显著下降, 表明CO2能有效地改善ILs扩散性, 减小其粘度. 因此CO2可用以改善离子液体溶剂体系的传递特性, 增强反应分离过程在其中的进行.     

A new force field model of 1-butyl-3-methylimidazolium tetrafluoroborate ionic liquid and acetonitrile mixtures

Recently, we introduced a new force field (FF) to simulate transport properties of imidazolium-based room-temperature ionic liquids (RTILs) using a solid physical background. In the present work, we apply this FF to derive thermodynamic, structure, and transport properties of the mixtures of 1-butyl-3-methylimidazolium tetrafluoroborate, [BMIM][BF(4)], and acetonitrile (ACN) over the whole composition range. Three approaches to derive a force field are formulated based on different treatments of the ion-ion and ion-molecule Coulomb interactions: unit-charge, scaled-charge and floating-charge approaches. The simulation results are justified with the help of experimental data on specific density and shear viscosity for these mixtures. We find that a phenomenological account (particularly, a simple scaled-charge model) of electronic polarization leads to the best-performing model. Remarkably, its validity does not depend on the molar fraction of [BMIM][BF(4)] in the mixture. The derived FF is so far the first molecular model which is able to simulate all transport properties of the mixtures, comprising RTIL and ACN, fully realistically.     

Molecular dynamics simulations of gas separations using faujasite-type zeolite membranes

Gas separations with faujasite zeolite membranes have been examined using the method of molecular dynamics. Two binary mixtures are investigated, oxygen/nitrogen and nitrogen/carbon dioxide. These mixtures have been found experimentally to exhibit contrasting behavior. In O(2)/N(2) mixtures the ideal selectivity (pure systems) is higher than the mixture selectivity, while in N(2)/CO(2) the mixture selectivity is higher than the ideal selectivity. One of the key goals of this work was to seek a fundamental molecular level understanding of such divergent behavior. Our simulation results (using previously developed intermolecular models for both the gases and zeolites investigated) were found to replicate this experimental behavior. By examining the loading of the membranes and the diffusion rates inside the zeolites, we have been able to explain such contrasting behavior of O(2)/N(2) and N(2)/CO(2) mixtures. In the case of O(2)/N(2) mixtures, the adsorption and loading of both O(2) and N(2) in the membrane are quite competitive, and thus the drop in the selectivity in the mixture is primarily the result of oxygen slowing the diffusion of nitrogen and nitrogen somewhat increasing the diffusion of oxygen when they pass through the zeolite pores. In N(2)/CO(2) systems, CO(2) is rather selectively adsorbed and loaded in the zeolite, leaving very little room for N(2) adsorption. Thus although N(2) continues to have a higher diffusion rate than CO(2) even in the mixture, there are so few N(2) molecules in the zeolite in mixtures that the selectivity of the mixture increases significantly compared to the ideal (pure system) values. We have also compared simulation results with hydrodynamic theories that classify the permeance of membranes to be either due to surface diffusion, viscous flow, or Knudsen diffusion. Our results show surface diffusion to be the dominant mode, except in the case of N(2)/CO(2) binary mixtures where Knudsen diffusion also makes a contribution to N(2) transport.     

Modified free volume theory of self-diffusion and molecular theory of shear viscosity of liquid carbon dioxide

In previous work on the density fluctuation theory of transport coefficients of liquids, it was necessary to use empirical self-diffusion coefficients to calculate the transport coefficients (e.g., shear viscosity of carbon dioxide). In this work, the necessity of empirical input of the self-diffusion coefficients in the calculation of shear viscosity is removed, and the theory is thus made a self-contained molecular theory of transport coefficients of liquids, albeit it contains an empirical parameter in the subcritical regime. The required self-diffusion coefficients of liquid carbon dioxide are calculated by using the modified free volume theory for which the generic van der Waals equation of state and Monte Carlo simulations are combined to accurately compute the mean free volume by means of statistical mechanics. They have been computed as a function of density along four different isotherms and isobars. A Lennard-Jones site-site interaction potential was used to model the molecular carbon dioxide interaction. The density and temperature dependence of the theoretical self-diffusion coefficients are shown to be in excellent agreement with experimental data when the minimum critical free volume is identified with the molecular volume. The self-diffusion coefficients thus computed are then used to compute the density and temperature dependence of the shear viscosity of liquid carbon dioxide by employing the density fluctuation theory formula for shear viscosity as reported in an earlier paper (J. Chem. Phys. 2000, 112, 7118). The theoretical shear viscosity is shown to be robust and yields excellent density and temperature dependence for carbon dioxide. The pair correlation function appearing in the theory has been computed by Monte Carlo simulations.     

A new gas chromatographic methodology for the estimation of the composition of binary gas mixtures

The relatively new technique of reversed-flow gas chromatography (RFGC) is used to determine the diffusion coefficients of pure gases into gas mixtures (D(mix)(exp)). The pure gases are CO and CO(2), and the mixtures consist of H(2) and He in various volume percentage compositions. A linear regression analysis of D(mix)(exp) of CO and CO(2) in various mixtures of H(2) and He against the percentage composition (X(H2) or X(He)) of the mixtures at different temperatures results in an empirical equation relating D(mix)(exp) to the corresponding theoretical values of the diffusion coefficients of CO and CO(2) in the pure gases H(2) and He, as they are calculated from the Fuller-Schettler-Giddings equation. The empirical equation shows that the diffusion coefficient of an analyte gas in a gas mixture is the partial sum of its diffusion coefficients in the component gases, therefore making possible the determination of the mole fractions of the components of the mixture. The found percentage volume compositions are very close to those determined independently by routine gas chromatography, indicating that the proposed RFGC methodology could be successfully applied to the accurate determination of the volume composition of binary gas mixtures.     

Quasi-elastic neutron scattering study of dimethyl-sulfoxide-water mixtures: probing molecular mobility in a nonideal solution

The translational and rotational motions of water and dimethyl sulfoxide, [DMSO, (CH(3))(2)SO] have been investigated using quasi-elastic neutron scattering. Water-DMSO mixtures at five DMSO mole fractions, chi(DMSO), ranging from 0 to 0.75, were measured. Hydrogen-deuterium substitution was used to extract independently the water proton dynamics (d-DMSO-H(2)O), the DMSO methyl proton dynamics (h-DMSO-D(2)O) and to obtain background corrections (d-DMSO-D(2)O). The translational diffusion of water slows down significantly compared to bulk water at all chi(DMSO)>0. The rotational time constant for water exhibits a maximum at chi(DMSO)=0.33 that corresponds to the observed maximum of the viscosity of the mixture. Data for DMSO can be analyzed in terms of a relatively slow tumbling of the molecule about its center-of-mass in conjunction with random translational diffusion. The rotational time constant for this motion exhibits some dependence on chi(DMSO), while the translational diffusion constant shows no clear variation for chi(DMSO)>0. The results presented reinforce the idea that due to the stronger associative nature of DMSO, DMSO-water aggregates are formed over the whole composition range, disturbing the tetrahedral natural arrangement of the water molecules. As a consequence adding DMSO to water causes a drastic slowing down of the dynamics of the water molecule, and vice versa.     

Effects of molecular association on polarizability relaxation in liquid mixtures of benzene and hexafluorobenzene

In this work we have studied the relaxation dynamics of the many-body polarizability anisotropy in liquid mixtures of benzene (Bz) and hexafluorobenzene (Hf) at room temperature by femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES) experiments and molecular dynamics (MD) simulations. The computed polarizability response arising from intermolecular interactions was included using the first-order dipole-induced-dipole model with the molecular polarizability distributed over the carbon sites of each molecule. We found good qualitative agreement between experiments and simulations in the features exhibited by the nuclear response function R(t) for pure liquids and mixtures. The long-time diffusive decay of R(t) was observed to vary substantially with composition, slowing down noticeably with dilution of each of the species as compared with that in the corresponding pure liquids. MD simulation shows that the effect on R(t) is due to the formation of strong and localized intermolecular association between Bz and Hf species that hinder the rotational diffusive dynamics. The formation of these Bz-Hf complexes in the liquid mixtures also modifies the rotational diffusive dynamics of the component species in such a way that cannot be explained solely in terms of a viscosity effect. Even though the computed orientational diffusive relaxation times associated with Bz and Hf are larger by a factor of approximately 2 than those from experiments, we found similar trends in experiments and simulations for these characteristic times as a function of composition. Namely, the collective and single-molecule orientational correlation times associated with Bz are observed to grow monotonically with the dilution of Bz, while those corresponding to Hf species exhibit a maximum at the equimolar composition. We attribute the quantitative discrepancy between experiments and simulations to the use of the Williams potential, which seems to overestimate the intermolecular interactions and thus predicts not only a slower translational dynamics but also a slower rotational diffusion dynamics than in real fluids.     

Deep eutectic solvents formed between choline chloride and carboxylic acids: versatile alternatives to ionic liquids

Deep Eutectic Solvents (DES) can be formed between a variety of quaternary ammonium salts and carboxylic acids. The physical properties are significantly affected by the structure of the carboxylic acid but the phase behavior of the mixtures can be simply modeled by taking account of the mole fraction of carboxylic acid in the mixture. The physical properties such as viscosity, conductivity, and surface tension of these DES are similar to ambient temperature ionic liquids and insight into the cause of these properties is gained using hole-theory. It is shown that the conductivity and viscosity of these liquids is controlled by ion mobility and the availability of voids of suitable dimensions, and this is consistent with the fluidity of other ionic liquids and molten salts. The DES are also shown to be good solvents for metal oxides, which could have potential application for metal extraction.     

Low pressure solubility and thermodynamics of solvation of oxygen,carbon dioxide,and carbon monoxide in fluorinated liquids

The solubility of oxygen, carbon dioxide, and carbon monoxide in three fluorinated liquids – perfluorohexylethane, perfluorooctane and bromoperfluorooctane – is presented. Mole fraction solubilities were calculated from new experimental Ostwald coefficient data for CO₂ and CO, and from previously published values for O₂, associated with original values of density and vapour pressure for the pure solvents. Carbon dioxide is the most soluble gas with mole fraction solubilities of the order of 10⁻². Oxygen and carbon monoxide are one order of magnitude less soluble. The measurements were done as a function of temperature between (288 and 313) K and from the variation of the calculated Henry’s law constants with temperature, the thermodynamic properties of solvation such as the Gibbs free energy, the enthalpy and the entropy were calculated. The precision of the experimental data, considered as the average absolute deviation of the Henry’s law constants from appropriate smoothing equations is of 1% for carbon dioxide and oxygen and of 3% for carbon monoxide. The data obtained here are judged accurate to within ±5%.     

Mutual diffusion in binary mixtures of ionic liquids and molecular liquids by dynamic light scattering (DLS)

The present study shows that dynamic light scattering (DLS) is capable of measuring mutual diffusion coefficients for binary mixtures of ionic liquids (ILs) with different molecular liquids over the complete composition range. Evidence is given that the light scattering signals are related to true molecular binary diffusion. The method stands out due to its ability to work non-invasively in macroscopic thermodynamic equilibrium with reasonable accuracy and within convenient measurement periods. Compared with other techniques, mixtures with distinctly higher viscosities can be probed. For exemplary binary mixtures of 1-ethyl-3-methylimidazolium ethylsulfate ([EMIM][EtSO(4)]) with acetone, acetonitrile, dichloromethane, ethanol, or water as well as of 1-ethyl-3-methylimidazolium methanesulfonate ([EMIM][MeSO(3)]) with acetone, water, or methanol, mutual diffusivity data were measured over a wide range of composition at a temperature of 293.15 K. In general, the mutual diffusivity increases with increasing mole fraction of the molecular liquid and similarities to aqueous solutions of classical inorganic salts can be found. The characteristic behavior of the mutual diffusion coefficients is influenced by the nature of the chosen molecular liquid. For IL water mixtures, low light scattering intensities were observed despite the large refractive index difference of the pure components. The reason for this behavior may be the existence of water clusters in the mixtures. Additional measurements for IL acetone mixtures at temperatures ranging from 278.15 K to 323.15 K showed that the temperature dependence of the mutual diffusivity can be represented by Arrhenius functions and is increasing for decreasing mole fractions of acetone.     

Relaxation of the optical Kerr effect of anisotropic molecules in mixed liquids

The relaxation time of the Kerr effect of nitrobenzene and m-nitrotoluene in various mixtures with carbon tetrachloride and various alcohols was determined by measuring the kinetics of the Kerr effect using picosecond laser techniques. These measurements yield information on the rotational motion of molecules in liquids. The relaxation time data are interpreted in terms of an effective local viscosity effect, pair correlation, and coupling of rotational motion with shear modes.     

A predictive approach based on the Simha–Somcynsky free‐volume theory for the effect of dissolved gas on viscosity and glass transition temperature of polymeric mixtures

The gas concentration and pressure effects on the shear viscosity of molten polymers were modeled by using a unified approach based on a free volume theory. A concentration and pressure dependent “shift factor,” which accounts for free volume changes associated with polymer‐gas mixing and with variation of absolute pressure as well as for dilution effects, has been herein used to scale the pure polymer viscosity, as evaluated at the same temperature and atmospheric pressure. The expression of the free volume of the polymer/gas mixture was obtained by using the Simha and Somcynsky equation of state for multicomponent fluids. Experimental shear viscosity data, obtained for poly(ε‐caprolactone) with nitrogen and carbon dioxide were successfully predicted by using this approach. Good agreement with predictions was also found in the case of viscosity data reported in the literature for polystyrene and poly(dimethylsiloxane) with carbon dioxide. Free volume arguments have also been used to predict the T_g depression for polystyrene/carbon dioxide and for poly(methyl methacrylate)/carbon dioxide mixtures, based on calculations performed, again, with the Simha and Somcynsky theory. © 2006 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 44: 1863–1873, 2006     

Correlation between the fluorescent response of microfluidity probes and the water content and viscosity of ionic liquid and water mixtures

Accurate data on transport properties such as viscosity are essential in plant and process design involving ionic liquids. In this study, we determined the absolute viscosity of the ionic liquid + water system at water mole fractions from 0 to 0.25 for three 1-alkyl-3-methylimidazolium ionic liquids: 1-butyl-3-methylimidazolium hexafluorophosphate, 1-butyl-3-methylimidazolium bis(trifluoromethane sulfonyl)imide and 1-ethyl-3-methylimidazolium bis(trifluoromethane sulfonyl)imide. In each case, the excimer to monomer ratio for 1,m-bis(1-pyrenyl)alkanes (m= 3 or 10) was found to increase linearly with the mole fraction of water. Of the probes studied only PRODAN and rhodamine 6G, both of which have the ability to participate in hydrogen bonding, exhibited Perrin hydrodynamic behavior in the lower viscosity bis(trifluoromethane sulfonyl)imides. As a result, these probes allow for the extrapolation of the absolute viscosity of the ionic liquid mixture from the experimental fluorescence steady-state polarization values.     

Ion association in [bmim][PF6]/naphthalene mixtures: an experimental and computational study

Mixtures of room temperature ionic liquids (IL) with neutral organic molecules provide a valuable testing ground to investigate the interplay of the ionic and molecular-dipolar state in dense Coulomb systems at near ambient conditions. In the present study, the viscosity eta and the ionic conductivity sigma of 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6])/naphthalene mixtures at T = 80 degrees C have been measured at 10 stoichiometries spanning the composition range from pure naphthalene to pure [bmim][PF6]. The viscosity grows nearly monotonically with increasing IL mole fraction ( x), whereas the conductivity per ion displays a clear peak at x approximately 15%. The origin of this maximum has been investigated using molecular dynamics simulations based on a classical force field. Snapshots of the simulated samples show that the conductivity maximum is due to the gradual transition in the IL component from an ionic state at high x to a dipolar fluid made of neutral ion pairs at low x. At concentrations x < 0.20 the ion pairs condense into molecular-thin filaments bound by dipolar forces and extending in between nanometric droplets of IL. These results are confirmed and complemented by the computation of dynamic and transport properties in [bmim][PF6]/naphthalene mixtures at low IL concentration.     

Coupling of translational and rotational motion in chiral liquids in electromagnetic and circularly polarised electric fields

Non-equilibrium molecular dynamics simulations of R and S enantiomers of 1,1-chlorofluoroethane, both for pure liquids and racemic mixtures, have been performed at 298 K in the absence and presence of both electromagnetic (e/m) and circularly polarised electric (CP) fields of varying frequency (100-2200 GHz) and intensity (0.025-0.2 V ?(-1) (rms)). Significant non-thermal field effects were noted in the coupling of rotational and translational motion; for instance, in microwave and far-infrared (MW/IR) e/m fields, marked increases in rotational and translational diffusion vis-à-vis the zero-field case took place at 0.025-0.1 V ?(-1) (rms), with a reduction in translational diffusion vis-à-vis the zero-field case above 0.1 V ?(-1) (rms) above 100 GHz. This was due to enhanced direct coupling of rotational motion with the more intense e/m field at the ideal intrinsic rotational coupling frequency (approximately 700 GHz) leading to such rapidly oscillating rotational motion that extent of translational motion was effectively reduced. In the case of CP fields, rotational and translational diffusion was also enhanced for all intensities, particularly at approximately 700 GHz. For both MW/IR and CP fields, non-linear field effects were evident above around 0.1 V ?(-1) (rms) intensity, in terms of enhancements in translational and rotational motion. Simulation of 90-10 mol. % liquid mixtures of a Lennard-Jones solvent with R and S enantiomer-solutes in MW/IR and CP fields led to more limited promotion of rotational and translational diffusion, due primarily to increased frictional effects. For both e/m and CP fields, examination of the laboratory- and inertial-frame auto- and cross-correlation functions of velocity and angular velocity demonstrated the development of explicit coupling with the external fields at the applied frequencies, especially so in the more intense fields where nonlinear effects come into play. For racemic mixtures, elements of the laboratory- and inertial-frame velocity and angular velocity were found to couple with each other to a lesser extent.     

Ideal gas solubilities and solubility selectivities in a binary mixture of room-temperature ionic liquids

This study focuses on the solubility behaviors of CO2, CH4, and N2 gases in binary mixtures of imidazolium-based room-temperature ionic liquids (RTILs) using 1-ethyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([C2mim][Tf2N]) and 1-ethyl-3-methylimidazolium tetrafluoroborate ([C2mim][BF4]) at 40 degrees C and low pressures (approximately 1 atm). The mixtures tested were 0, 25, 50, 75, 90, 95, and 100 mol % [C2mim][BF4] in [C2mim][Tf2N]. Results show that regular solution theory (RST) can be used to describe the gas solubility and selectivity behaviors in RTIL mixtures using an average mixture solubility parameter or an average measured mixture molar volume. Interestingly, the solubility selectivity, defined as the ratio of gas mole fractions in the RTIL mixture, of CO2 with N2 or CH4 in pure [C2mim][BF4] can be enhanced by adding 5 mol % [C2mim][Tf2N].     

Phase equilibria in carbon dioxide expanded solvents: Experiments and molecular simulations

We present complementary molecular simulations and experimental results of phase equilibria for carbon dioxide expanded acetonitrile, methanol, ethanol, acetone, acetic acid, toluene, and 1-octene. The volume expansion measurements were done using a high-pressure Jerguson view cell. Molecular simulations were performed using the Gibbs ensemble Monte Carlo method. Calculations in the canonical ensemble (NVT) were performed to determine the coexistence curve of the pure solvent systems. Binary mixtures were simulated in the isobaric-isothermal distribution (NPT). Predictions of vapor-liquid equilibria of the pure components agree well with experimental data. The simulations accurately reproduced experimental data on saturated liquid and vapor densities for carbon dioxide, methanol, ethanol, acetone, acetic acid, toluene, and 1-octene. In all carbon dioxide expanded liquids (CXL's) studied, the molecular simulation results for the volume expansion of these binary mixtures were found to be as good as, and in many cases superior to, predictions based on the Peng-Robinson equation of state, demonstrating the utility of molecular simulation in the prediction of CXL phase equilibria.     

Picosecond time-resolved fluorescence study on solute-solvent interaction of 2-aminoquinoline in room-temperature ionic liquids: aromaticity of imidazolium-based ionic liquids

Time-resolved fluorescence spectra and fluorescence anisotropy decay of 2-aminoquinoline (2AQ) have been measured in eight room-temperature ionic liquids, including five imidazolium-based aromatic ionic liquids and three nonaromatic ionic liquids. The same experiments have also been carried out in several ordinary molecular liquids for comparison. The observed time-resolved fluorescence spectra indicate the formation of pi-pi aromatic complexes of 2AQ in some of the aromatic ionic liquids but not in the nonaromatic ionic liquids. The fluorescence anisotropy decay data show unusually slow rotational diffusion of 2AQ in the aromatic ionic liquids, suggesting the formation of solute-solvent complexes. The probe 2AQ molecule is likely to be incorporated in the possible local structure of ionic liquids, and hence the anisotropy decays only through the rotation of the whole local structure, making the apparent rotational diffusion of 2AQ slow. The rotational diffusion time decreases rapidly by adding a small amount of acetonitrile to the solution. This observation is interpreted in terms of the local structure formation in the aromatic ionic liquids and its destruction by acetonitrile. No unusual behavior upon addition of acetonitrile has been found for the nonaromatic ionic liquids. It is argued that the aromaticity of the imidazolium cation plays a key role in the local structure formation in imidazolium-based ionic liquids.