Molecular Transport Behavior of CO2 in Ionic Polyimides and Ionic Liquid Composite Membrane Materials.

Ionic polyimides (i-PI) are a new class of polymer materials that are very promising for CO2 capture membranes, and recent experimental studies have demonstrated their enhanced separation performance with the addition of imidazolium-based ionic liquids (ILs). However, there is very little known about the molecular-level interactions in these systems, which give rise to interesting gas adsorption and diffusion characteristics. In this study, we use a combination of Monte Carlo and molecular dynamics simulations to analyze the equilibrium and transport properties of CO2 molecules in the i-PI and i-PI + IL composite materials. The addition of several different common ILs are modeled, which have a plasticization effect on the i-PI, lowering the glass transition temperature (Tg). The solubility of CO2 strongly correlates with the Tg, but the diffusion demonstrates more unpredictable behavior. At low concentrations, the IL has a blocking effect, leading to reduced diffusion rates. However, as the IL surpasses a threshold value, the relationship is inverted and the IL has a facilitating effect on the gas transport. This behavior is attributed to the simultaneous contributions of the increased i-PI plasticization at higher IL concentrations (facilitating gas hopping rates from cavity-to-cavity) and the increased IL continuity throughout the system, enabling more favorable transport pathways for CO2 diffusion.

Polycyclic aromatic hydrocarbons as model solutes for carbon nanomaterials in ionic liquids.

The aim of this work is to understand the details of the interactions of ionic liquids with carbon nanomaterials (graphene and nanotubes) using polyaromatic compounds as model solutes. We have combined the measurements of thermodynamic quantities of solvation with molecular dynamics simulations to provide a microscopic view. The solubility of five polycyclic aromatic hydrocarbons (naphthalene, anthracene, phenanthrene, pyrene and coronene) was determined in seven ionic liquids ([C4C1im][C(CN)3], [C4C1pyrr][Ntf2], [C10C1im][Ntf2], [C2C1im][C(CN)3], [C2C1im][Ntf2], [C3C1pyrr][N(CN)2] and [C4C1im][N(CN)2]) at 298 K. The enthalpies of the dissolution of naphthalene, anthracene and pyrene were measured in four of the ionic liquids. Free energies were estimated from those measurements in order to analyse the entropic or enthalpic contributions to the dissolution process. Molecular dynamics simulations provided solvation free energies that were compared to experimental and structural information. Spatial distributions of solvent ions around the solutes when combined with IR measurements elucidate the structure of solvation environments. Interactions between the imidazolium rings of cations and the π system of the solutes have been identified. However, ionic liquids with pyrrolidinium cations appeared as better solvents due to favourable enthalpic contributions compared to imidazolium cations. Long alkyl side chains on cations lead to higher solubility and lower enthalpy of dissolution by creating a "softer" solvation environment. Considering the effect of anions, small and planar anions lead to higher solubilities and lower enthalpies of dissolution of polyaromatic hydrocarbons. These findings provide the design principles based on molecular interactions and the structure of solvation environments to choose or formulate ionic liquids in view of their affinity for carbon nanomaterials.

Exfoliation of graphene and fluorographene in molecular and ionic liquids.

We use molecular dynamics simulations to study the exfoliation of graphene and fluorographene in molecular and ionic liquids, by performing computer experiments in which one layer of the 2D nanomaterial is peeled from a stack, in vacuum and in the presence of solvent. The liquid media and the nanomaterials are represented by fully flexible, atomistic force fields. From these simulations we calculate the potential of mean force, or reversible work, required to exfoliate the materials. Calculations in water and organic liquids showed that small amides (NMP, DMF) are among the best solvents for exfoliation, in agreement with the experiment. We tested ionic liquids with different cation and anion structures, allowing us to learn about their solvent qualities for the exfoliation of the nanomaterials. First, a long alkyl side chain on the cation is favourable for exfoliation of both graphene and fluorographene. The presence of aromatic groups on the cation is also favourable for graphene. No beneficial effect was found between fluorine-containing anions and fluorographene. We also analysed the ordering of ions in the interfacial layers with the materials. Near graphene, nonpolar groups are found along with charged groups, whereas near fluorographene almost exclusively non-charged groups are found, with ionic moieties segregated to the second layer. Therefore, fluorographene appears to be the more hydrophobic surface, as expected.

The Structure of Ethylbenzene, Styrene and Phenylacetylene Determined by Total Neutron Scattering.

Organic solvents such as phenylacetylene, styrene and ethylbenzene are widely used in industrial processes, especially in the production of rubber or thermoplastics. Despite their important applications detailed knowledge about their structure is limited. In this paper the structures of these three aromatic solvents were investigated using neutron diffraction. The results show that many of their structural characteristics are similar, although the structure of phenylacetylene is more ordered and has a smaller solvation sphere than either ethylbenzene or styrene. Two regions within the first coordination sphere, in which the surrounding molecules show different preferable orientations with respect to the central molecule, were found for each liquid. Additionally, the localisation of the aliphatic chains reveals that they tend to favour closer interactions with each other than to the aromatic rings of the adjacent molecules.

Mechanism of OH radical hydration: a comparative computational study of liquid and supercritical solvent.

Flexible models of the radical and water molecules including short-range interaction of hydrogen atoms have been employed in molecular dynamic simulation to understand mechanism of (●)OH hydration in aqueous systems of technological importance. A key role of H-bond connectivity patterns of water molecules has been identified. The behavior of (●)OH(aq) strongly depends on water density and correlates with topological changes in the hydrogen-bonded structure of water driven by thermodynamic conditions. Liquid and supercritical water above the critical density exhibit the radical localization in cavities existing in the solvent structure. A change of mechanism has been found at supercritical conditions below the critical density. Instead of cavity localization, we have identified accumulation of water molecules around (●)OH associated with the formation of a strong H-donor bond and diminution of non-homogeneity in the solvent structure. For all the systems investigated, the computed hydration number and the internal energy of hydration Δ(h)U showed approximately linear decrease with decreasing density of the solvent but a degree of radical-water hydrogen bonding exhibited non-monotonic dependence on density. The increase in the number of radical-water H-acceptor bonds is associated with diminution of extended nets of four-bonded water molecules in compressed solution at ~473 K. Up to 473 K, the isobaric heat of hydration in compressed liquid water remains constant and equal to -40 ± 1 kJ mol(-1).

Molecular dynamics study of the hydration of the hydroxyl radical at body temperature.

Classical molecular dynamics (MD) simulation of ˙OH in liquid water at 37 °C has been performed using flexible models of the solute and solvent molecules. We derived the Morse function describing the bond stretching of the radical and the potential for ˙OH-H(2)O interactions, including short-range interactions of hydrogen atoms. Scans of the potential energy surface of the ˙OH-H(2)O complex have been performed using the DFT method with the B3LYP functional and the 6-311G(d,p) basis set. The DFT-derived partial charges, ±0.375e, and the equilibrium bond-length, 0.975 Å, of ˙OH resulted in the dipole moment of 1.76 D. The radical-water radial distribution functions revealed that ˙OH is not built into the solvent structure but it rather occupies distortions or cavities in the hydrogen-bonded network. The solvent structure at 37 °C has been found to be the same as that of pure water. The hydration cage of the radical comprises 13-14 water molecules. The estimated hydration enthalpy -42 ± 5 kJ mol(-1) is comparable with the experimental value -39 ± 6 kJ mol(-1) for 25 °C. Inspection of hydrogen bonds showed the importance of short-range interaction of hydrogen atoms and indicated that neglect of the angular condition greatly overestimates the number of the H-acceptor radical-water bonds. The mean number ̅n = 0.85 of radical-water H-bonds has been calculated using geometric definition of H-bond and ̅n = 0.62 has been obtained when the energetic condition, E(da)≤-8 kJ mol(-1), was additionally considered. The continuous lifetimes of 0.033 ps and 0.024 ps have been estimated for the radical H-donor and the H-acceptor bonds, respectively. Within statistical uncertainty the radical self-diffusion coefficient, (2.9 ± 0.6) × 10(-9) m(2) s(-1), is the same as (3.1 ± 0.5) × 10(-9) m(2) s(-1) calculated for water in solution and in pure solvent. To the best of our knowledge, this is the first study of the ˙OH(aq) properties at a biologically relevant body temperature.

Transition from patchlike to clusterlike inhomogeneity arising from hydrogen bonding in water.

Assembling of water molecules via hydrogen bonding has been studied by molecular dynamics simulations using flexible potential model. The relationship between the number of H-bonds per molecule, n(HB), the size of H-bonded nets, k, and the size of patches of four-bonded molecules, k(4), has been examined for several thermodynamic states of water ranging from ambient to supercritical conditions. Two kinds of structural inhomogeneity have been found: the patchlike associated with the mean n(HB)> 2.0 and the clusterlike for n(HB)< 1.9. In compressed water up to ~473 K patches coexist with less ordered nets, both constituting the gel-like H-bonded network. The size of patches steeply decreases with the increasing temperature and the decreasing density of water. The inhomogeneity resulting from the presence of patches disappears above 473 K. This feature is associated with the rapid increase in the fraction of unbound molecules and with the breakage of the gel-like network into a variety of H-bonded clusters leading to the clusterlike structural inhomogeneity. In contrast to the patchlike inhomogeneity an increase in temperature and a decrease in density make this kind of inhomogeneity more pronounced. A degree of connectivity of H-bonds has been characterized by a parameter P(g) defined as the total fraction of molecules belonging to the H-bonded clusters of size k ≥ 5. The simulation-derived values of P(g) agree well with the predictions of the random bond theory giving the explicit expression for P(g) as a function of the mean n(HB). Going from ambient to supercritical conditions, we have found that the patchlike inhomogeneity is connected with the very slight reduction in P(g), whereas the clusterlike inhomogeneity generates a steep linear decrease of P(g) with the decreasing mean n(HB). The self-diffusion coefficient calculated for the thermodynamic states of water showing the clusterlike inhomogeneity has occurred to be inversely proportional to the density. We have also found that the clusterlike inhomogeneity is associated with the linear correlation between P(g) and the macroscopic properties of water: the static dielectric constant, the viscosity, and the density. The provided relationships allow one to estimate the degree of connectivity of hydrogen bonds from the measured macroscopic quantities.