Adsorption Kinetics in Nanoscale Porous Coordination Polymers.

Nanoscale porous coordination polymers were synthesized using simple wet chemical method. The effect of various polymer surfactants on colloidal stability and shape selectivity was investigated. Our results suggest that the nanoparticles exhibited significantly improved adsorption kinetics compared to bulk crystals due to decreased diffusion path lengths and preferred crystal plane interaction.

Separation of polar compounds using a flexible metal-organic framework.

A flexible metal-organic framework constructed from a flexible linker is shown to possess the capability of separating mixtures of polar compounds (propanol isomers) by exploiting the differences in the saturation capacities of the constituents. Transient breakthrough simulations show that these sorption-based separations are in favor of the component with higher saturation capacity.

A Combined Experimental and Computational Study on the Stability of Nanofluids Containing Metal Organic Frameworks.

Computational studies on nanofluids composed of metal organic frameworks were performed using molecular modeling techniques. Grand Canonical Monte Carlo simulations were used to study the adsorption behavior of 1,1,1,3,3-pentafluoropropane (R-245fa) in a MIL-101 metal organic frameworks at various temperatures. To understand the stability of the nanofluid composed of MIL-101 particles, we performed molecular dynamics simulations to compute potentials of mean force between hypothetical MIL-101 fragments terminated with two different kinds of modulators in R-245fa and water. Our computed potentials of mean force results indicate that the metal organic frameworks particles tend to disperse better in water than in R-245fa. The reasons for this difference in dispersion were analyzed and are discussed in the paper. Our results agree with experimental results indicating that the potential models employed and modeling approaches provide good descriptions of molecular interactions and the reliabilities.

Fluorocarbon adsorption in hierarchical porous frameworks.

Metal-organic frameworks comprise an important class of solid-state materials and have potential for many emerging applications such as energy storage, separation, catalysis and bio-medical. Here we report the adsorption behaviour of a series of fluorocarbon derivatives on a set of microporous and hierarchical mesoporous frameworks. The microporous frameworks show a saturation uptake capacity for dichlorodifluoromethane of >4 mmol g(-1) at a very low relative saturation pressure (P/Po) of 0.02. In contrast, the mesoporous framework shows an exceptionally high uptake capacity reaching >14 mmol g(-1) at P/Po of 0.4. Adsorption affinity in terms of mass loading and isosteric heats of adsorption is found to generally correlate with the polarizability and boiling point of the refrigerant, with dichlorodifluoromethane > chlorodifluoromethane > chlorotrifluoromethane > tetrafluoromethane > methane. These results suggest the possibility of exploiting these sorbents for separation of azeotropic mixtures of fluorocarbons and use in eco-friendly fluorocarbon-based adsorption cooling.

Understanding the rates and molecular mechanism of water-exchange around aqueous ions using molecular simulations.

Solvation processes occurring around aqueous ions are of fundamental importance in physics, chemistry, and biology. Over the past few decades, several experimental and theoretical studies were devoted to understanding ion solvation and the processes involved in it. In this article, we present a summary of our recent efforts that, through computer simulations, focused on providing a comprehensive understanding of solvent-exchange processes around aqueous ions. To accomplish these activities, we have looked at the mechanistic properties associated with the water-exchange process, such as potentials of mean force, time-dependent transmission coefficients, and the corresponding rate constants using transition state theory, the reactive flux method, and Grote-Hynes treatments of the dynamic response of the solvent.

Water exchange rates and molecular mechanism around aqueous halide ions.

Molecular dynamics simulations were performed to systematically study the water-exchange mechanism around aqueous chloride, bromide, and iodide ions. Transition state theory, Grote-Hynes theory, and the reactive flux method were employed to compute water exchange rates. We computed the pressure dependence of rate constants and the corresponding activation volumes to investigate the mechanism of the solvent exchange event. The activation volumes obtained using the transition state theory rate constants are negative for all the three anions, thus indicating an associative mechanism. Contrary to the transition state theory results, activation volumes obtained using rate constants from Grote-Hynes theory and the reactive flux method are positive, thus indicating a dissociative mechanism.

Computational studies of water exchange around aqueous Li+ with polarizable potential models.

To enhance our understanding of the mechanism of water-exchange around aqueous Li(+), we carried out a systematic study on this system using molecular dynamics simulations with polarizable potential models. The mechanistic properties associated with the water-exchange process, such as potentials of mean force, time dependent transmission coefficients, and the corresponding rate constants, were examined using transition rate theory, the reactive flux method, and Grote-Hynes treatments of the dynamic response of the solvent. We compared the computed rate theory results with results from previous corresponding studies in which classical non-polarizable force fields were used. Our computed barrier heights for water exchange are significantly larger than those obtained using classical non-polarizable force fields. We also studied the effect of pressure on water-exchange rates and the corresponding activation volume. Our computed rate results for water exchange increase with pressure; therefore, a small negative activation volume is observed.

Understanding ion-ion interactions in bulk and aqueous interfaces using molecular simulations.

In addition to its scientific significance, the distribution of ions in the bulk and at aqueous interfaces is also very important for practical reasons. Providing a quantitative description of the ionic distribution, and describing interactions between ions in different environments, remains a challenge, and is the subject of current debate. In this study, we found that interionic potentials of mean force (PMFs) and interfacial properties are very sensitive to the ion-ion interaction potential models. Our study predicted a Sr(2+)--CI- PMF with no contact ion-pair state and a shallow solvent-separated ion-pair state. In addition, we were able to quantitatively capture the experimental X-ray reflectivity results of the aqueous salt interface of the Sr(2+)--Cl- ion-pair, and provided a detailed physical description of the interfacial structure for this system. We also predicted the Xray reflectivity results for SrBr2 and SrI2 systems.

Pairing mechanism among ionic liquid ions in aqueous solutions: a molecular dynamics study.

In this study, we carried out molecular dynamics simulations to examine the molecular mechanism for ionic liquid pair association in aqueous solutions. We chose the commonly studied imidazolium-based ionic liquid pairs. We computed potentials of mean force (PMF) for four systems: 1,3-dimethlylimidazoliumchloride, 1,3-dimethlylimidazolium iodide, 1-methly-3-octylimidazolium chloride, and 1-methly-3-octylimidazolium iodide. Our PMF studies show a stronger interaction for the ion pairs of systems involving dimethlylimidazolium as the cation species compared with that of the systems containing octylimidazolium. This result indicates a decrease in ion-pair association as the cation alkyl tail length increases. We also studied the kinetics of ion-pair dissociation using different rate theories such as the Grote-Hynes and Kramer's theories. As expected, the computed rate results significantly deviated from results obtained from transition state theory because it does not account for dynamical solvent effects. Dissociative barrier curvatures are found to be very small for the systems investigated because the transmission coefficients computed using Grote-Hynes theory and Kramer's theory are approximately equal. Our analysis of the rotational dynamics of cations revealed that the time scales for molecular reorientation are longer for cations with longer alkyl tails.

Anions, the Reporters of Structure in Ionic Liquids.

In this work we compare the role that different anions play in the structure function S(q) for a set of liquids with the same cation. It is well established that because of their amphiphilic nature and their often larger size, cations play a fundamental role in the structural landscape of ionic liquids. On the other hand, it is often atoms in the anions that display the largest X-ray form factors and therefore play a very significant role as reporters of structure in small- and wide-angle X-ray scattering (SAXS/WAXS)-type experiments. For a set of liquids with similar topological landscape, how does S(q) change when the anionic scattering is deemphasized? Also, how do we computationally recover the typical length scale of important and perhaps universal ionic liquid structural features such as charge alternation when these are experimentally inaccessible from S(q) because of interference cancellations? We answer these questions by studying three different tetrapentylammonium-based liquids with the I(-), PF6(-) and N(CN)2(-) anions.

Molecular mechanism of the adsorption process of an iodide anion into liquid-vapor interfaces of water-methanol mixtures.

To enhance our understanding of the molecular mechanism of ion adsorption to the interface of mixtures, we systematically carried out a free energy calculations study involving the transport of an iodide anion across the interface of a water-methanol mixture. Many body affects are taken into account to describe the interactions among the species. The surface propensities of I(-) at interfaces of pure water and methanol are well understood. In contrast, detailed knowledge of the molecular level adsorption process of I(-) at aqueous mixture interfaces has not been reported. In this paper, we explore how this phenomenon will be affected for mixed solvents with varying compositions of water and methanol. Our potential of mean force study as function of varying compositions indicated that I(-) adsorption free energies decrease from pure water to pure methanol but not linearly with the concentration of methanol. We analyze the computed density profiles and hydration numbers as a function of concentrations and ion positions with respect to the interface to further explain the observed phenomenon.

Molecular mechanism of specific ion interactions between alkali cations and acetate anion in aqueous solution: a molecular dynamics study.

Specific ion interactions between alkali cations (i.e., Li(+), Na(+), and K(+)) and an acetate anion in aqueous solution were studied using molecular dynamics simulation techniques and polarizable potential models. The ions-acetate systems were used as a model for understanding the interactions between ions and protein surfaces. We computed free energy profiles for different ion pairs using constrained mean force methods. Upon analyzing the computed free energy profiles for the Na(+)/K(+)-acetate ion-pairs, we observed a deeper contact ion minimum and also a larger association constant for the Na(+)-acetate pair as compared to the corresponding K(+)-acetate pair. These observations help to demonstrate the preferential binding of Na(+) over K(+) to protein surfaces. We also applied various rate theories to study the kinetics of ion pair interconversion.

SAXS anti-peaks reveal the length-scales of dual positive-negative and polar-apolar ordering in room-temperature ionic liquids.

Structural patterns that have the same spatial periodicity but a phase offset give rise to peaks and anti-peaks (negative-going peaks) at the same q value in the SAXS structure function S(q). As an example, in ionic liquids we often find charge alternation, and at the distance where one finds a density enhancement of charges of the same type one also finds a depletion of charges of opposite sign. Another such situation arises with polar-apolar densities. At distances where there is enhancement of same-type (polar-polar or apolar-apolar) densities there is also a depletion of opposite-type (polar-apolar) density. This gives rise to prepeaks and what we call same spatial periodicity anti-prepeaks.

Temperature-dependent structure of ionic liquids: X-ray scattering and simulations.

In this article we determine the temperature-dependent structure of the tetradecyltrihexylphosphonium bis(trifluoromethylsulfonyl)amide ionic liquid using a combination of X-ray scattering and molecular dynamics simulations. As in many other room-temperature ionic liquids three characteristic intermolecular peaks can be detected in the structure function S(q). A prepeak or first sharp diffraction peak is observed at about q = 0.42 A(-1). Long range anion-anion correlations are the most important contributors to this peak. In all systems we have studied to date, this prepeak is a signature of solvation asymmetry. The peak in S(q) near q = 0.75 A(-1) is the signature of ionic alternation and arises from the charge ordered separation of ions of the same charge. The most intense diffraction peak near q = 1.37 A(-1) arises from short-range separation between ions of opposite charge combined with a significant contribution from cationic carbon-carbon interactions, indicating that cationic hydrophobic tails have significant contacts.

Dry excess electrons in room-temperature ionic liquids.

In this article, we describe general trends to be expected at short times when an excess electron is generated or injected in different room-temperature ionic liquids (RTILs). Perhaps surprisingly, the excess electron does not localize systematically on the positively charged cations. Rather, the excess charge localization pattern is determined by the cation and anion HOMO/LUMO gaps and, more importantly, by their relative LUMO alignments. As revealed by experiments, the short-time (ps/ns) transient UV spectrum of excess electrons in RTILs is often characterized by two bands, a broad band at low energies (above 1000 nm) and another weaker band at higher energies (around 400 nm). Our calculations show that the dry or presolvated electron spectrum (fs) also has two similar features. The broad band at low energies is due to transitions between electronic states with similar character on ions of the same class but in different locations of the liquid. The lower-intensity band at higher energies is due to transitions in which the electron is promoted to electronic states of different character, in some cases on counterions. Depending on the chemical nature of the RTIL, and especially on the anions, excess electrons can localize on cations or anions. Our findings hint at possible design strategies for controlling electron localization, where electron transfer or transport across species can be facilitated or blocked depending on the alignment of the electronic levels of the individual species.

How is charge transport different in ionic liquids and electrolyte solutions?

In this article we show that, analyzed in a barycentric reference frame, the deviation in conductivity measured directly from impedance experiments with respect to that estimated indirectly from NMR diffusion experiments has different origins in electrolyte solutions and pure salts. In the case of electrolyte solutions, the momentum conservation law is satisfied by solvent + ions. Instead, in a molten salt or ionic liquid momentum conservation must be satisfied solely by the ions. This has significant implications. While positively correlated motion of ions of opposite charge is a well justified explanation for the reduction in impedance conductivity in the case of electrolyte solutions, it is not so in the case of ionic liquids and molten salts. This work presents a set of equations that in the case of ionic liquids and molten salts can be used to obtain from direct measurements of impedance and NMR the distinct part of the diffusion coefficient matrix in the barycentric reference frame. In other words, by using experimentally measurable quantities, these equations allow us to access the motional coupling between ions for which there is no single direct experimental measurement technique. While equations of this type have been proposed before, the ones presented here can be easily derived from the momentum conservation law and linear response theory. Our results indicate that the decrease in the impedance conductivity with respect to NMR conductivity in ionic liquids and molten salts is due to anticorrelated motion of ions of same charge. This scenario is different in electrolyte solutions, where the positively correlated motion of ions of opposite charge makes a significant contribution to the decrease in the impedance conductivity. In contrast, in a system comprising a single binary salt (a room temperature ionic liquid or a molten salt), the cation-anion distinct diffusion coefficient is negative definite and opposes the contribution from the cation-cation and anion-anion distinct diffusion coefficients. This property of the cation-anion distinct diffusion coefficient in systems comprising just two ion-constituents holds true not just in the barycentric reference frame but also in any of the internal reference frames of nonequilibrium thermodynamics.

Temperature-dependent structure of methyltributylammonium bis(trifluoromethylsulfonyl)amide: X ray scattering and simulations.

We report the combined results of computational and x ray scattering studies of amorphous methyltributylammonium bis(trifluoromethylsulfonyl)amide as a function of temperature. These studies included the temperature range for the normal isotropic liquid, a deeply supercooled liquid and the glass. The low q peaks in the range from 0.3 to 1.5 A(-1) in the structure function of this liquid can be properly accounted for by correlations between first and second nearest neighbors. The lowest q peak can be assigned to real space correlations between ions of the same charge, while the second peak arises mostly from nearest neighbors of opposite charge. Peaks at larger q values are mostly intramolecular in nature. While our simulated structure functions provide an excellent match to our experimental results and our experimental findings agree with previous studies reported for this liquid, the prior interpretation of the experimental data in terms of an interdigitated smectic A phase is not supported by our simulations. In this work, we introduce a set of general theoretical partitions of real and reciprocal space correlations that allow for unambiguous analysis of all intra- and interionic contributions to the structure function and coherent scattering intensity. We find that the intermolecular contributions to the x ray scattering intensity are dominated by the anions and cross terms between cations and anions for this ionic liquid.

What is the origin of the prepeak in the X-ray scattering of imidazolium-based room-temperature ionic liquids?

The observation of a first sharp diffraction peak (FSDP) at low frequency in the X-ray and neutron scattering spectra of different imidazolium-based room-temperature ionic liquids (RTILs) (the so-called prepeak) has often been experimentally interpreted as indicative of mesoscopic organization leading to nanoscale segregation and the formation of domains of different morphologies. This interpretation that has permeated the analysis of many recently published articles deserves an in depth theoretical analysis. In this article, we use several different computational techniques to thoroughly dissect the atomistic components giving rise to the low-frequency FSDP as well as other features in the structure function (S(q)). By understanding how S(q) changes as imidazolium-based ionic systems undergo solid-liquid phase transition, and by artificially perturbing the liquid structure in a way that directly couples to the intensity of the FSDP, we are able to identify in a rigorous way its geometric origin. Similar to the solid phase, the liquid phase is characterized by two typical length scales between polar groups. The shorter length scale gives rise to a shoulder peak in S(q) at about 0.9 Å(-1) whereas the longer one gives rise to the prepeak.

Controlling the outcome of electron transfer reactions in ionic liquids.

In this article, we investigate the excited state intramolecular electron transfer (ET) reaction of crystal violet lactone (CVL) in the room temperature ionic liquid (RTIL) N-propyl-N-methylpyrrolidinium bis(trifluoromethylsulfonyl) imide [Pr(31)(+)][Tf(2)N(-)]. This system was chosen in light of recent experimental observations by Maroncelli and co-workers (J. Phys. Chem. B 2007, 111, 13473), in which the kinetics of electron transfer between S(1) (commonly referred as LE) and S(2) (commonly referred as CT) emission states and, therefore, the ratio of emitting populations were shown to be absorption-wavelength-dependent. Our computational studies indicate that the kinetics of the intramolecular ET between S(1) and S(2) states of CVL in [Pr(31)(+)][Tf(2)N(-)] is local solvent-environment-dependent. Because emission time scales are smaller than solvent relaxation time scales, this behavior is characteristic of RTILs but uncommon in conventional solvents. Therefore, RTILs open a window of opportunity for manipulating the outcome of chemical reactions simply by tunning the initial excitation wavelength. Our studies show that when acetonitrile is used as a solvent instead of [Pr(31)(+)][Tf(2)N(-)] the ratio of populations of emission states is independent of excitation wavelength, eliminating the opportunity for influencing the outcome of reactions.

Molecular dynamics study of the temperature-dependent Optical Kerr effect spectra and intermolecular dynamics of room temperature ionic liquid 1-methoxyethylpyridinium dicyanoamide.

We have performed classical molecular dynamics simulations to calculate the Optical Kerr effect (OKE) spectra of 1-methoxyethylpyridinium dicyanoamide, a room-temperature ionic liquid (IL) which has been recently studied by Shirota and Castner (Shirota, H. ; Castner, E. J. Phys. Chem. A 2005, 109, 9388-9392) in comparison to its neutral isoelectronic solvent mixture. Our theoretical and computational studies show that the decay of the collective polarizability anisotropy correlation exhibits several different time scales originating from inter- and intramolecular dynamics, in good agreement with experiments. What's more, we find that the portion of the collective anisotropic polarizability relaxation due to "interaction-induced" phenomena is important at times much longer than those observed in normal solvents when these are far from their glass transition temperature. From our long (60 ns) molecular dynamics simulations, we are able to determine the appropriate time scales for orientational relaxation and interaction-induced processes occurring in the liquid. We find that the cationic contribution to the OKE signal is predominant. Because of the slow nature of relaxation processes in ILs, these calculations are very time, memory, and storage intensive. In the context of this research, we have developed a polarizable force field for this system and also theoretical methodology to generate molecular polarizabilities for arbitrarily shaped molecules and ions from corresponding atomic polarizabilities. We expect this methodology to have an important impact on the speed of molecular dynamics simulations of polarizable systems in the future.