Edward J. Maginn

Assessing the factors responsible for ionic liquid toxicity to aquatic organisms via quantitative structure–property relationship modeling

Using previously published toxicity data as well as a small set of heretofore-unpublished results, quantitative structure–property relationship models are developed to assess the factors that govern the toxicity of a range of different ionic liquids to two aquatic organisms (Vibrio fischeri and Daphnia magna). With at most four molecular descriptors, log10 EC50 and log10 LC50 data are reproduced with an R2 of 0.78–0.88. Besides the well-established link between toxicity and alkyl chain length on imidazolium, pyridinium and quaternary ammonium-based ionic liquids, the models predict that toxicity increases slightly with the number of nitrogen atoms in an aromatic cation ring. All other things being equal, toxicity is expected to show the trend with cation type of ammonium < pyridinium < imidazolium < triazolium < tetrazolium. In addition, toxicity is expected to decrease with ring methylation as well as with an increase in the number of negatively charged atoms on the cation. The anion plays a secondary role in toxicity for the compounds studied here, although the presence of positively charged atoms on the anion are predicted to slightly increase toxicity.

Amine-Functionalized Task-Specific Ionic Liquids: A Mechanistic Explanation for the Dramatic Increase in Viscosity upon Complexation with CO2 from Molecular Simulation

The capture of CO2 from fossil fuel combustion, particularly in coal-fired power plants, represents a critical component of efforts aimed at stabilizing greenhouse gas levels in the atmosphere. Alkanolamines have traditionally been used to this end; however, drawbacks such as volatility, degradation, and regeneration costs have been drivers for the development of new, superior technologies. Recently, several seminal studies with ionic liquids (ILs), both experimental and computational, have demonstrated their potential as CO2 capture agents. In traditional ILs, experimental studies with CO2 have revealed its unusually high physical solubility in these media. Complementary simulation studies have provided evidence that this is attributable to CO2 occupying void space within the liquid and favorably interacting with the anion. Recently, a series of second-generation task-specific ionic liquids (TSILs) containing amine functional groups have been synthesized and demonstrated to have much higher capacities for CO2 due to their reactivity with CO2, as well unusually high viscosities in both the neat and complexed states. The current work extends the seminal studies of CO2 capture with ILs by providing insight from simulations into the mechanism responsible for the dramatic increase in viscosity upon complexation. Simulations conclusively demonstrate that the slow translational and rotational dynamics, which are manifest in the high viscosity, may be attributable to the formation of a strong, pervasive hydrogen-bonded network. Semiquantitative estimates of the cation and anion self-diffusion coefficients and rotational time constants, as well as detailed hydrogen bond analysis, are consistent with the experimentally observed formation of glassy or gel-like materials upon contact with CO2. This has significant implications for the design of new approaches or materials involving ILs that take advantage of these preconceived limitations, in the synthesis or manipulation of new TSIL frameworks for CO2 capture, and in novel experimental studies of chemistries and dynamics in persistent heterogeneous environments.

Molecular simulation of ionic liquids: current status and future opportunities

Ionic liquids are salts that are liquid near ambient conditions. Interest in these unusual compounds has exploded in the last decade, both at the academic and commercial level. Molecular simulations based on classical potentials have played an important role in helping researchers understand how condensed phase properties of these materials are linked to chemical structure and composition. Simulations have also predicted many properties and unexpected phenomena that have subsequently been confirmed experimentally. The beneficial impact molecular simulations have had on this field is due in large part to excellent timing. Just when computing power and simulation methods matured to the point where complex fluids could be studied in great detail, a new class of materials virtually unknown to experimentalists came on the scene and demanded attention. This topical review explores some of the history of ionic liquid molecular simulations, and then gives examples of the recent use of molecular dynamics and Monte Carlo simulation in understanding the structure of ionic liquids, the sorption of small molecules in ionic liquids, the nature of ionic liquids in the vapor phase and the dynamics of ionic liquids. This review concludes with a discussion of some of the outstanding problems facing the ionic liquid modeling community and how condensed phase molecular simulation experts not presently working on ionic liquids might help advance the field.

Nanoscale design to enable the revolution in renewable energy

The creation of a sustainable energy generation, storage, and distribution infrastructure represents a global grand challenge that requires massive transnational investments in the research and development of energy technologies that will provide the amount of energy needed on a sufficient scale and timeframe with minimal impact on the environment and have limited economic and societal disruption during implementation. In this opinion paper, we focus on an important set of solar, thermal, and electrochemical energy conversion, storage, and conservation technologies specifically related to recent and prospective advances in nanoscale science and technology that offer high potential in addressing the energy challenge. We approach this task from a two-fold perspective: analyzing the fundamental physicochemical principles and engineering aspects of these energy technologies and identifying unique opportunities enabled by nanoscale design of materials, processes, and systems in order to improve performance and reduce costs. Our principal goal is to establish a roadmap for research and development activities in nanoscale science and technology that would significantly advance and accelerate the implementation of renewable energy technologies. In all cases we make specific recommendations for research needs in the near-term (2–5 years), mid-term (5–10 years) and long-term (>10 years), as well as projecting a timeline for maturation of each technological solution. We also identify a number of priority themes in basic energy science that cut across the entire spectrum of energy conversion, storage, and conservation technologies. We anticipate that the conclusions and recommendations herein will be of use not only to the technical community, but also to policy makers and the broader public, occasionally with an admitted emphasis on the US perspective.

Thermodynamic properties of the ionic liquid 1-n-butyl-3-methylimidazolium hexafluorophosphate from monte Carlo simulations

We report results from the first molecular simulation study of 1-n-butyl-3-methylimidazolium hexafluorophosphate [bmim][PF6], a widely studied ionic liquid. Monte Carlo simulations are carried out in the isothermal–isobaric (NPT) ensemble to calculate the molar volume, cohesive energy density, isothermal compressibility, cubic expansion coefficient and liquid structure as a function of temperature and pressure. A united atom forcefield is developed using a combination of ab initio calculations and literature parameter values. Calculated molar volumes are within 5% of experimental values, and a reasonable agreement is obtained between calculated and experimental values of the isothermal compressibility and cubic expansion coefficient. PF6 anions are found to preferentially cluster in two favorable regions near the cation.

Atomistic Simulation of the Absorption of Carbon Dioxide and Water in the Ionic Liquid 1-n-Hexyl-3-methylimidazolium Bis(trifluoromethylsulfonyl)imide ([hmim][Tf2N]

The solubility of water and carbon dioxide in the ionic liquid 1-n-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([hmim][Tf2N]) is computed using atomistic Monte Carlo simulations. A newly developed biasing algorithm is used to enable complete isotherms to be computed. In addition, a recently developed pairwise damped electrostatic potential calculation procedure is used to speed the calculations. The computed isotherms, Henrys Law constants, and partial molar enthalpies of absorption are all in quantitative agreement with available experimental data. The simulations predict that the excess molar volume of CO2/ionic liquid mixtures is large and negative. Analysis of ionic liquid conformations shows that the CO2 does not perturb the underlying liquid structure until very high CO2 concentrations are reached. At the highest CO2 concentrations, the alkyl chain on the cation stretches out slightly, and the distance between cation and anion centers of mass increases by about 1 angstroms. Water/ionic liquid mixtures have excess molar volumes that are also negative but much smaller in magnitude than those for the case of CO2.

A biased grand canonical Monte Carlo method for simulating adsorption using all-atom and branched united atom models

Configurational-bias Monte Carlo sampling techniques have been developed which overcome the difficulties of sampling configuration space efficiently for all-atom molecular models and for branched species represented with united atom models. Implementation details of this sampling scheme are discussed. The accuracy of a united atom forcefield with non-bond parameters optimized for zeolite adsorption and a widely used all-atom forcefield are evaluated by comparison with experimental sorption isotherms of linear and branched hydrocarbons.

A critical comparison of equilibrium, non-equilibrium and boundary-driven molecular dynamics techniques for studying transport in microporous materials

Transport in an idealized model with variable pore diameter as well as an AlPO4-5 zeolite is examined using three different molecular dynamics techniques: (1) equilibrium molecular dynamics (EMD); (2) external field nonequilibrium molecular dynamics (EF–NEMD); and (3) dual control volume grand canonical molecular dynamics (DCV–GCMD). The EMD and EF–NEMD methods yield identical transport coefficients for all the systems studied. The transport coefficients calculated using the DCV–GCMD method, however, tend to be lower than those obtained from the EMD and EF–NEMD methods unless a large ratio of stochastic to dynamic moves is used for each control volume, and a streaming velocity is added to all inserted molecules. Through development and application of a combined reaction–diffusion–convection model, this discrepancy is shown to be due to spurious mass and momentum transfers caused by the control volume equilibration procedure. This shortcoming can be remedied with a proper choice of streaming velocity in co...

A molecular dynamics investigation of the structural and dynamic properties of the ionic liquid 1-n-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide

Molecular dynamics simulations have been performed to investigate the structure and dynamics of the ionic liquid, 1-n-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide ([C(4)mim][Tf(2)N]) in the temperature range of 283-460 K. Extensive analysis was carried out to characterize a number of structural and dynamic features. Transport properties were computed using a variety of equilibrium methods that employed the Green-Kubo and Einstein formulations. Nonequilibrium techniques were also used. In general, different methods mostly yielded consistent results, although some differences were observed. Computed self-diffusivities and ionic conductivities tended to be slightly lower than experimental values, while computed viscosities were significantly higher than experiment. Computed thermal conductivities agreed reasonably well with experimental data. Despite these discrepancies, the simulations capture the experimental temperature-dependent trends for all these transport properties. Single ion dynamics were studied by examining diffusional anisotropy, the self-part of the van Hove function, non-Gaussian parameters, and incoherent intermediate scattering functions. It is found that cations diffuse faster than anions and are more dynamically heterogeneous. A clear anisotropy is revealed in cation displacement, with the motion normal to the imidazolium ring plane being the most hindered and the motion along the alkyl chain in the plane of the ring being the most facile. Cations structurally relax faster than anions but they rotationally relax slower than anions. There is a pronounced temperature dependence to the rotational anisotropy of the cations, but only a weak temperature dependence for the anions. The ionic conductivity deviates from the Nernst-Einstein relation due to the correlated motion of cations and anions. The results suggest that the dynamical behavior of this and related ionic liquids is extremely complex and consists of many different modes with widely varying timescales, making the prediction of dynamical trends extremely difficult.

Making Sense of Enthalpy of Vaporization Trends for Ionic Liquids: New Experimental and Simulation Data Show a Simple Linear Relationship and Help Reconcile Previous Data

Vaporization enthalpy of an ionic liquid (IL) is a key physical property for applications of ILs as thermofluids and also is useful in developing liquid state theories and validating intermolecular potential functions used in molecular modeling of these liquids. Compilation of the data for a homologous series of 1-alkyl-3-methylimidazolium bis(trifluoromethane-sulfonyl)imide ([C(n)mim][NTf2]) ILs has revealed an embarrassing disarray of literature results. New experimental data, based on the concurring results from quartz crystal microbalance, thermogravimetric analyses, and molecular dynamics simulation have revealed a clear linear dependence of IL vaporization enthalpies on the chain length of the alkyl group on the cation. Ambiguity of the procedure for extrapolation of vaporization enthalpies to the reference temperature 298 K was found to be a major source of the discrepancies among previous data sets. Two simple methods for temperature adjustment of vaporization enthalpies have been suggested. Resulting vaporization enthalpies obey group additivity, although the values of the additivity parameters for ILs are different from those for molecular compounds.