Molecular simulation of osmometry in aqueous solutions of the BMIMCl ionic liquid: a potential route to force field
1Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA. ddhabal@seas.upenn.edu debiit09@gmail.com.
Abstract:
Despite widespread development and use of ionic liquids (ILs) in both academic and industrial research, computational force fields (FFs) for most of those are not available for a precise description of inter-species interactions in aqueous environments. In the scope of this study, by means of molecular simulations, the osmotic coefficient of an aqueous solution of an IL is calculated and used as a basis to reparameterize popular IL-FFs existing in the literature. We first calculate the osmotic coefficients (at 298.15 K and 1 atm pressure) of aqueous solutions of 1-butyl-3-methylimidazolium chloride (BMIMCl), a generic IL, popularly used in biomass processing and the subsequent conversion to value-added intermediates. The performance of two popular atomic, nonpolarizable FFs developed for BMIMCl, one by Lopes, Pádua, and coworkers (FF-LP) and the other by Sambasivarao, Acevedo, and coworkers (FF-SA), when mixed with the SPC/E water model, is tested with respect to their ability to reproduce the experimental osmotic coefficient data. Interestingly, the osmotic coefficient is found to be increasing with a gradual increase in IL molality within the concentration range of our investigation, which is contrary to the experimental trend reported in the literature for the same IL-water mixture. Henceforth, necessary corrections to the nonbonded ion-ion and ion-water interactions are made to match the experimental osmotic coefficient. To further assess the reliability of the new FF, we extensively explore the thermodynamic (density, isothermal compressibility, and thermal expansion coefficient), dynamic (diffusivity and viscosity), and association/dissociation properties (rationalized with the help of radial distribution functions) with both the original and reparameterized FF for a wider range of concentrations up to a molality of 18.50 mol kg-1. The calculated quantities are compared against experimental data wherever available. The modified FF parameters exhibit significant improvements in terms of its ability to match experimental solution properties, such as density, viscosity, association/dissociation, etc. We report that excessive dissociation of BMIMCl in water is responsible for the shortcomings observed in the original FFs and improved prediction of physicochemical properties could be achieved using the modified FFs.
More Related Videos
08:54Vibrational Spectra of a N719-Chromophore/Titania Interface from Empirical-Potential Molecular-Dynamics Simulation, Solvated by a Room Temperature Ionic Liquid
Published on: January 25, 2020
10:52Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics
Published on: April 12, 2019
Related Concept Videos
03:13Intermolecular Forces
02:56Intermolecular Forces and Physical Properties
02:28Intermolecular Forces in Solutions
When the strengths of the intermolecular forces of attraction between solute and solvent species in a solution are no different than those present in the separated components, the solution is formed with no accompanying energy change. Such a solution is called an ideal solution. A mixture of ideal gases (or gases such as helium and argon,...
02:34Comparing Intermolecular Forces: Melting Point, Boiling Point, and Miscibility
Temporary attractive forces like dispersion are present in all molecules, whether they are polar or nonpolar. They...
03:12Freezing Point Depression and Boiling Point Elevation
The boiling point of a liquid is the temperature at which its vapor pressure is equal to ambient atmospheric pressure. Since the vapor pressure of a solution is lowered due to the presence of nonvolatile solutes, it stands to reason that the solution’s boiling point will subsequently be increased. Vapor pressure increases with temperature, and so a solution will require a higher temperature than will pure solvent to achieve any given vapor pressure, including one...
01:34The Fluid Mosaic Model