Search PubMed⌕ Search

Biomedical subjects

Edward J Maginn

Publications and source records attributed to Edward J Maginn.

13 recordsLinked to original sources

Molecular simulation study of some thermophysical and transport properties of triazolium-based ionic liquids.

Results of a molecular dynamics study of several triazolium-based ionic liquids are reported. Triazolium cations include 1,2,4-triazolium, 1,2,3-triazolium, 4-amino-1,2,4-triazolium, and 1-methyl-4-amino-1,2,4-triazolium. Each cation was paired with a nitrate or perchlorate anion. These materials are part of a class of ionic compounds that have been synthesized recently but for which little physical property data are available. Properties of the more common ionic liquid, 1-n-butyl-3-methylimidazolium nitrate, are also computed and compared with the properties of the triazolium-based compounds. A molecular mechanics force field was developed for these materials using a mix of ab initio calculations and parameter fitting using the molecular compound 1H-1,2,4-triazole as a basis for the triazolium cations. Liquid-phase properties that were computed include heat capacities, cohesive energy densities, gravimetric densities/molar volumes as a function of temperature and pressure, self-diffusivities, rotational time constants, and various pair correlation functions. In the solid phase, heat capacities and lattice parameters were computed. Of all of these properties, only lattice parameters have been measured experimentally (and only for four of the triazolium compounds). The agreement with the experimental crystal structures was good. When compared with that of the imidazolium-based ionic liquid, the triazolium-based materials have much smaller molar volumes, higher cohesive energy densities, and larger specific heat capacities. They also tend to be less compressible, have a higher gravimetric density, and have faster rotational dynamics but similar translational dynamics.

Journal Article↗

Measurement of SO2 solubility in ionic liquids.

Measurements of the solubility of sulfur dioxide (SO(2)) in the ionic liquids 1-n-hexyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide ([hmim][Tf(2)N]) and 1-n-hexyl-3-methylpyridinium bis(trifluoromethylsulfonyl)imide ([hmpy][Tf(2)N]) at temperatures from 25 to 60 degrees C and pressures up to 4 bar indicate that large amounts (up to 85 mol %) of SO(2) dissolve in ionic liquids by simple physical absorption.

Imidazoles↗

Liquid phase behavior of ionic liquids with alcohols: experimental studies and modeling.

Ionic liquids (ILs) have been suggested as potential "green" solvents to replace volatile organic solvents in reaction and separation processes due to their negligible vapor pressure. To develop ILs for these applications, it is important to gain a fundamental understanding of the factors that control the phase behavior of ionic liquids with other liquids. In this work, we continue our study of the effect of chemical and structural factors on the phase behavior of ionic liquids with alcohols, focusing on pyridinium ILs for comparison to imidazolium ILs from our previous studies. The impact of different alcohol and IL characteristics, including alcohol chain length, cation alkyl chain length, anion, different substituent groups on the pyridinium cation, and type of cation (pyridinium vs imidazolium) will be discussed. In general, the same type of behavior is observed for pyridinium and imidazolium ILs, with all systems studied exhibiting upper critical solution temperature behavior. The impacts of alcohol chain length, cation chain length, and anion, are the same for pyridinium ILs as those observed previously for imidazolium ILs. However, the effect of cation type on the phase behavior is dependent on the strength of the cation-anion interaction. Additionally, all systems from this study and our previous work for imidazolium ILs were modeled using the nonrandom two-liquid (NRTL) equation using two different approaches for determining the adjustable parameters. For all systems, the NRTL equation with binary interaction parameters with a linear temperature dependence provided a good fit of the experimental data.

Journal Article↗

Atomistic simulation of solid-liquid coexistence for molecular systems: application to triazole and benzene.

A method recently developed to rigorously determine solid-liquid equilibrium using a free-energy-based analysis has been extended to analyze multiatom molecular systems. This method is based on using a pseudosupercritical transformation path to reversibly transform between solid and liquid phases. Integration along this path yields the free energy difference at a single state point, which can then be used to determine the free energy difference as a function of temperature and therefore locate the coexistence temperature at a fixed pressure. The primary extension reported here is the introduction of an external potential field capable of inducing center of mass order along with secondary orientational order for molecules. The method is used to calculate the melting point of 1-H-1,2,4-triazole and benzene. Despite the fact that the triazole model gives accurate bulk densities for the liquid and crystal phases, it is found to do a poor job of reproducing the experimental crystal structure and heat of fusion. Consequently, it yields a melting point that is 100 K lower than the experimental value. On the other hand, the benzene model has been parametrized extensively to match a wide range of properties and yields a melting point that is only 20 K lower than the experimental value. Previous work in which a simple "direct heating" method was used actually found that the melting point of the benzene model was 50 K higher than the experimental value. This demonstrates the importance of using proper free energy methods to compute phase behavior. It also shows that the melting point is a very sensitive measure of force field quality that should be considered in parametrization efforts. The method described here provides a relatively simple approach for computing melting points of molecular systems.

Algorithms↗

Molecular modeling and experimental studies of the thermodynamic and transport properties of pyridinium-based ionic liquids.

A combined experimental and molecular dynamics study has been performed on the following pyridinium-based ionic liquids: 1-n-hexyl-3-methylpyridinium bis(trifluoromethanesulfonyl)imide ([hmpy][Tf(2)N]), 1-n-octyl-3-methylpyridinium bis(trifluoromethanesulfonyl)imide ([ompy][Tf(2)N]), and 1-n-hexyl-3,5-dimethylpyridinium bis(trifluoromethanesulfonyl)imide ([hdmpy][Tf(2)N]). Pulsed field gradient nuclear magnetic resonance spectroscopy was used to determine the self-diffusivities of the individual cations and anions as a function of temperature. Experimental self-diffusivities range from 10(-11) to 10(-10) m(2)/s. Activation energies for diffusion are 44-49 kJ/mol. A classical force field was developed for these compounds, and molecular dynamics simulations were performed to compute dynamic as well as thermodynamic properties. Evidence of glassy dynamics was found, preventing accurate determination of self-diffusivities over molecular dynamics time scales. Volumetric properties such as density, isothermal compressibility, and volumetric expansivity agree well with experiment. Simulated heat capacities are within 2% of experimental values.

Computer Simulation↗

Rapid shear viscosity calculation by momentum impulse relaxation molecular dynamics.

Recently, Arya et al. [J. Chem. Phys. 113, 2079 (2000)] introduced a new molecular dynamics method to rapidly compute the viscosity of fluids. The technique, termed momentum impulse relaxation (MIR), involves the imposition of a Gaussian velocity profile on an equilibrated system, after which the decay in the profile is monitored as a function of time. The shear viscosity is computed by matching the rate of decay of the velocity profile to the corresponding solution of the Navier-Stokes equation. The method was originally applied to simple systems (argon and n-butane) and found to give a comparable accuracy to conventional equilibrium and nonequilibrium methods with more than an order of magnitude reduction in computing time. In this work, we extend and generalize the method to examine larger molecules with higher viscosities than have been examined previously. A detailed analysis of the method is given, including the effect the velocity boundary conditions have on the viscosity, the sensitivity of the results to the velocity profile fitting procedure, the effect of preequilibration of the Gaussian profile, and the effect the system size and box shape have on the accuracy and speed of the method. It is shown that the MIR method can be extended to treat multiatom systems without loss of accuracy or computational efficiency.

Journal Article↗

Monte Carlo simulations of gas solubility in the ionic liquid 1-n-butyl-3-methylimidazolium hexafluorophosphate.

The Henry's constants of water, carbon dioxide, ethane, ethene, methane, oxygen, and nitrogen are computed in the ionic liquid 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF(6)]) using test particle insertion and expanded ensemble Monte Carlo methods. The partial molar enthalpy and partial molar entropy of solvation are also computed for water, carbon dioxide, and oxygen. The results from the simulations are compared against experimental data from the literature. In addition, the accuracy and precision of the two methods in determining the Henry's constant are examined. Local organization of the ionic liquid around a solute molecule is analyzed, and the interactions responsible for the experimentally observed solubility trends are identified.

Journal Article↗

Anion effects on gas solubility in ionic liquids.

This work presents the results of solubility measurements for a series of gases in 1-n-butyl-3-methyl imidazolium tetrafluoroborate and 1-n-butyl-3-methyl imidazolium bis(trifluoromethylsulfonyl) imide. The gases considered include benzene, carbon dioxide, nitrous oxide, ethylene, ethane, oxygen, and carbon monoxide. Carbon dioxide and oxygen solubilities are also reported in methyl-tributylammonium bis(trifluoromethylsulfonyl) imide, butyl-methyl pyrrolidinium bis(trifluoromethylsulfonyl) imide, and tri-isobutyl-methyl phosphonium p-toluenesulfonate. We report the associated Henry's constants and enthalpies and entropies of absorption. In general, benzene, followed by carbon dioxide and nitrous oxide, have the highest solubilities and strongest interactions with the ionic liquids, followed by ethylene and ethane. Oxygen had very low solubilities and weak interactions. Carbon monoxide had a solubility below the detection limit of our apparatus. Ionic liquids with the bis(trifluoromethylsulfonyl) imide anion had the largest affinity for CO(2), regardless of whether the cation was imidazolium, pyrrolidinium, or tetraalkylammonium. These results suggest that the nature of the anion has the most significant influence on the gas solubilities.

Journal Article↗

Synthesis, structural characterization, and molecular modeling of dodecaniobate keggin chain materials.

Four new isostructural one-dimensional dodecaniobate Keggin materials, Na12[Ti2O2][TNb12O40] x xH2O and Na10[Nb2O2][TNb12O40] x xH2O with T = (Si or Ge), have been synthesized hydrothermally using a Lindqvist-ion salt, Na7[Nb6O19H] x 15H2O, as the precursor. Their structure, consisting of chains of Keggin ions [TNb12O40]16- linked by [Ti2O2]4+ or [Nb2O2]6+ bridges, was solved ab initio from powder diffraction data. The location of the charge-balancing sodium atoms and the water molecules was further investigated by molecular simulations. These compounds were also characterized by IR and solid-state 1H, 29Si, and 23Na MAS NMR spectroscopies. The structural relationships between these and related phases based on similar Keggin ion building units are discussed.

Journal Article↗

Isomolar semigrand ensemble molecular dynamics: development and application to liquid-liquid equilibria.

An extended system molecular dynamics method for the isomolar semigrand ensemble (fixed number of particles, pressure, temperature, and fugacity fraction) is developed and applied to the calculation of liquid-liquid equilibria (LLE) for two Lennard-Jones mixtures. The method utilizes an extended system variable to dynamically control the fugacity fraction xi of the mixture by gradually transforming the identity of particles in the system. Two approaches are used to compute coexistence points. The first approach uses multiple-histogram reweighting techniques to determine the coexistence xi and compositions of each phase at temperatures near the upper critical solution temperature. The second approach, useful for cases in which there is no critical solution temperature, is based on principles of small system thermodynamics. In this case a coexistence point is found by running N-P-T-xi simulations at a common temperature and pressure and varying the fugacity fraction to map out the difference in chemical potential between the two species A and B (mu(A)-mu(B)) as a function of composition. Once this curve is known the equal-distance/equal-area criterion is used to determine the coexistence point. Both approaches give results that are comparable to those of previous Monte Carlo (MC) simulations. By formulating this approach in a molecular dynamics framework, it should be easier to compute the LLE of complex molecules whose intramolecular degrees of freedom are often difficult to properly sample with MC techniques.

Journal Article↗

Toward a robust and general molecular simulation method for computing solid-liquid coexistence.

A rigorous and generally applicable method for computing solid-liquid coexistence is presented. The method overcomes some of the technical difficulties associated with other solid-liquid simulation procedures and can be implemented within either a molecular dynamics or Monte Carlo framework. The method consists of three steps: First, relative Gibbs free energy curves are created for the solid and liquid phases using histogram reweighting. Next, the free energy difference between the solid and liquid phases is evaluated at a single state point by integrating along a pseudosupercritical transformation path that connects the two phases. Using this result, the solid and liquid free energy curves are referenced to a common point, allowing a single coexistence point to be determined. Finally, Gibbs-Duhem integration is used to determine the full coexistence curve. To evaluate its utility, this method is applied to the Lennard-Jones and NaCl systems. Results for solid-liquid coexistence agree with previous calculations for these systems. In addition, it is shown that the NaCl model does not correctly describe solid-liquid coexistence at high pressures. An analysis of the accuracy of the method indicates that the results are most sensitive to the transformation free energy calculation.

Computer Simulation↗

Why Is CO2 so soluble in imidazolium-based ionic liquids?

Experimental and molecular modeling studies are conducted to investigate the underlying mechanisms for the high solubility of CO2 in imidazolium-based ionic liquids. CO2 absorption isotherms at 10, 25, and 50 degrees C are reported for six different ionic liquids formed by pairing three different anions with two cations that differ only in the nature of the "acidic" site at the 2-position on the imidazolium ring. Molecular dynamics simulations of these two cations paired with hexafluorophosphate in the pure state and mixed with CO2 are also described. Both the experimental and the simulation results indicate that the anion has the greatest impact on the solubility of CO2. Experimentally, it is found that the bis(trifluoromethylsulfonyl)imide anion has the greatest affinity for CO2, while there is little difference in CO2 solubility between ionic liquids having the tetrafluoroborate or hexafluorophosphate anion. The simulations show strong organization of CO2 about hexafluorophosphate anions, but only small differences in CO2 structure about the different cations. This is consistent with the experimental finding that, for a given anion, there are only small differences in CO2 solubility for the two cations. Computed and measured densities, partial molar volumes, and thermal expansion coefficients are also reported.

Journal Article↗

Knudsen diffusivity of a hard sphere in a rough slit pore.

An analytic theory for the Knudsen self-diffusivity D(s) of hard spheres in an atomically rough slit-shaped pore is presented which quantitatively matches simulation results. The theory assumes that, due to chaotic molecular trajectories caused by surface morphology, collisions of gas molecules with the wall are partly diffuse and partly specular, the relative magnitude of each depending upon the magnitude of the tangential momentum accommodation coefficient f. The theory thus represents a universal Knudsen fluctuation-dissipation correlation between longitudinal momentum loss and diffusivity that can simplify efforts to estimate D(s). It is also found that D(s) computed using Maxwell's theory of slip, in which collisions with the walls are assumed to be purely diffuse or specular, overpredicts the simulated D(s) by a large margin.

Journal Article↗