Molecular dynamics studies of melting and some liquid-state properties of 1-ethyl-3-methylimidazolium hexafluorophosphate [emim][PF6

Molecular dynamics simulations are used to study the liquid-state properties and melting of 1-ethyl-3-methylimidazolium hexafluorosphosphate [emim][PF6] using the force field of Canongia Lopes et al. [J. Phys. Chem. B 108, 2038 (2004)] and geometric constants from crystallographic data. The structures of the solid and liquid states are characterized by carbon-carbon, carbon-phosphorous, and phosphorous-phosphorous radial distribution functions. Spatial correlations among the ions are strong in the liquid state. The cohesive energy density and the temperature dependences of the molar volume and density of the liquid have been computed. The melting point is determined by equilibrating the solid-state supercells in which void defects have been introduced to eliminate the free-energy barrier for the formation of a solid-liquid interface. The computed melting point is 375+/-10 K, which is approximately 10% higher than the experimental value of 333 K.     

Molecular dynamics studies of melting and solid-state transitions of ammonium nitrate

Molecular dynamics simulations are used to calculate the melting point and some aspects of high-temperature solid-state phase transitions of ammonium nitrate (AN). The force field used in the simulations is that developed by Sorescu and Thompson [J. Phys. Chem. A 105, 720 (2001)] to describe the solid-state properties of the low-temperature phase-V AN. Simulations at various temperatures were performed with this force field for a 4 x 4 x 5 supercell of phase-II AN. The melting point of AN was determined from calculations on this supercell with voids introduced in the solid structure to eliminate superheating effects. The melting temperature was determined by calculating the density and the nitrogen-nitrogen radial distribution functions as functions of temperature. The melting point was predicted to be in the range 445 +/- 10 K, in excellent agreement with the experimental value of 442 K. The computed temperature dependences of the density, diffusion, and viscosity coefficient for the liquid are in good agreement with experiment. Structural changes in the perfect crystal at various temperatures were also investigated. The ammonium ions in the phase-II structure are rotationally disordered at 400 K. At higher temperatures, beginning at 530 K, the nitrate ions are essentially rotationally unhindered. The density and radial distribution functions in this temperature range show that the AN solid is superheated. The rotational disorder is qualitatively similar to that observed in the experimental phase-II to phase-I solid-state transition.     

Thermodynamic studies of ionic interactions in aqueous solutions of imidazolium-based ionic liquids [Emim][Br] and [Bmim][Cl

Experimental measurements of density at different temperatures ranging from 293.15 to 313.15 K, the speed of sound and osmotic coefficients at 298.15 K for aqueous solution of 1-ethyl-3-methylimidazolium bromide ([Emim][Br]), and osmotic coefficients at 298.15 K for aqueous solutions of 1-butyl-3-methylimidazolium chloride ([Bmim][Cl]) in the dilute concentration region are taken. The data are used to obtain compressibilities, expansivity, apparent and limiting molar properties, internal pressure, activity, and activity coefficients for [Emim][Br] in aqueous solutions. Experimental activity coefficient data are compared with that obtained from Debye-Hückel and Pitzer models. The activity data are further used to obtain the hydration number and the osmotic second virial coefficients of ionic liquids. Partial molar entropies of [Bmim][Cl] are also obtained using the free-energy and enthalpy data. The distance of the closest approach of ions is estimated using the activity data for ILs in aqueous solutions and is compared with that of X-ray data analysis in the solid phase. The measured data show that the concentration dependence for aqueous solutions of [Emim][Br] can be accounted for in terms of the hydrophobic hydration of ions and that this IL exhibits Coulombic interactions as well as hydrophobic hydration for both the cations and anions. The small hydration numbers for the studied ILs indicate that the low charge density of cations and their hydrophobic nature is responsible for the formation of the water-structure-enforced ion pairs.     

Extended Veytsman statistics for the solid–fluid transition in the framework of lattice fluid

Lattice fluid can describe a vapor–liquid transition but not a solid–fluid transition. In this work, we propose a simple and analytic term which yields a solid–fluid transition when coupled with a lattice based equation of state (EOS). The proposed term is derived based on the two assumptions that (1) solid can be considered as highly associated phase affected by strong attractive force and (2) this force is distinct from the conventional attractive forces yielding a vapor–liquid transition. To formulate these assumptions, we extend Veytsman statistics by modifying its density dependency. The derived term was combined with a quasi-chemical nonrandom lattice fluid theory (QLF) developed by the authors. The combined model was found to require only two parameters besides 3 QLF parameters for physical properties calculation of three phases. When tested against equilibrium properties of 8 components, the combined model was found to closely reproduce melting pressure, sublimation pressure, and vapor pressure, but underestimate solid density as well as heat of melting at the triple point temperature. It was found that the present approach can yield a solid–liquid transition at all temperatures.     

1‐Octanol/Water Partition Coefficients of 1Alkyl‐3‐methylimidazolium Chloride

The solubilities of 1alkyl‐3‐methylimidazolium chloride, [C_nmim][Cl], where n=4, 8, 10, and 12, in 1octanol and water have been measured by a dynamic method in the temperature range from 270 to 370 K. The solubility data was used to calculate the 1octanol/water partition coefficients as a function of temperature and alkyl substituent. The melting point, enthalpies of fusion, and enthalpies of solid–solid phase transitions were determined by differential scanning calorimetry, DSC. The solubility of [C_nmim][Cl], where n=10 or 12 in 1octanol is comparable and higher than that of [C₄mim][Cl] in 1octanol. Liquid 1n‐octyl‐3‐methylimidazolium chloride, [C₈mim][Cl], is not miscible with 1octanol and water, consequently, the liquid–liquid equilibrium, LLE was measured in this system. The differences between the solubilities in water for n=4 and 12 are shown only in α₁ and γ₁ solid crystalline phases. Additionally, the immiscibility region was observed for the higher concentration of [C₁₀mim][Cl] in water. The intermolecular solute–solvent interaction of 1butyl‐3‐methylimidazolium chloride with water is higher than for other 1alkyl‐3‐methylimidazolium chlorides. The data was correlated by means of the UNIQUAC ASM and two modified NRTL equations utilizing parameters derived from the solid–liquid equilibrium, SLE. The root‐mean‐square deviations of the solubility temperatures for all calculated data are from 1.8 to 7 K and depend on the particular equation used. In the calculations, the existence of two solid–solid first‐order phase transitions in [C₁₂mim][Cl] has also been taken into consideration. Experimental partition coefficients (log P) are negative at three temperatures; this is evidence for the possible use of these ionic liquids as green solvents.     

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.     

Molecular dynamics simulations of p-tert-butylcalix[4]arene with small guest molecules

Classical molecular dynamics simulations were used to study p-tert-butylcalix[4]arene inclusion compounds with xenon, nitrogen, hydrogen, methane, and sulfur dioxide guest molecules. The calixarene units were taken to be rigid and the intermolecular molecular interactions were modeled as a sum of the van der Waals interactions with parameters from the AMBER force field and electrostatic interactions. Simulations of the high-density alpha phase and low-density beta0 phase of p-tert-butylcalix[4]arene were used to test the force field. The predicted densities of the two phases were found to agree with experimental measurements at 173 K to within 5 %. Simulations were performed with guests placed inside the calixarene cages of the beta0 phase. Guest-host ratios of 1:1 to 1:4 were considered. Changes in the unit-cell volume and density of the phases with the addition of guest molecules and the inclusion energies for the guests were determined. Finally, the dynamics of the guest motion inside the cages were characterized by determining the root-mean-square displacements and velocity autocorrelation functions of the xenon and nitrogen guests.     

Molecular dynamics simulations of the melting of aluminum nanoparticles

Molecular dynamics simulations are performed to determine the melting points of aluminum nanoparticles of 55-1000 atoms with the Streitz-Mintmire [Phys. Rev. B 1994, 50, 11996] variable-charge electrostatic plus potential. The melting of the nanoparticles is characterized by studying the temperature dependence of the potential energy and Lindemann index. Nanoparticles with less than 850 atoms show bistability between the solid and liquid phases over temperature ranges below the point of complete melting. The potential energy of a nanoparticle in the bistable region alternates between values corresponding to the solid and liquid phases. This bistability is characteristic of dynamic coexistence melting. At higher temperatures, only the liquid state is stable. Nanoparticles with more than 850 atoms undergo a sharp solid-liquid-phase transition characteristic of the bulk solid phase. The variation of the melting point with the effective nanoparticle radius is also determined.     

A general purpose model for the condensed phases of water: TIP4P/2005

A potential model intended to be a general purpose model for the condensed phases of water is presented. TIP4P/2005 is a rigid four site model which consists of three fixed point charges and one Lennard-Jones center. The parametrization has been based on a fit of the temperature of maximum density (indirectly estimated from the melting point of hexagonal ice), the stability of several ice polymorphs and other commonly used target quantities. The calculated properties include a variety of thermodynamic properties of the liquid and solid phases, the phase diagram involving condensed phases, properties at melting and vaporization, dielectric constant, pair distribution function, and self-diffusion coefficient. These properties cover a temperature range from 123 to 573 K and pressures up to 40,000 bar. The model gives an impressive performance for this variety of properties and thermodynamic conditions. For example, it gives excellent predictions for the densities at 1 bar with a maximum density at 278 K and an averaged difference with experiment of 7 x 10(-4) g/cm3.     

Nanosolids, slushes, and nanoliquids: characterization of nanophases in metal clusters and nanoparticles

One of the keys to understanding the emergent behavior of complex materials and nanoparticles is understanding their phases. Understanding the phases of nanomaterials involves new concepts not present in bulk materials; for example, the phases of nanoparticles are quantum mechanical even when no hydrogen or helium is present. To understand these phases better, molecular dynamics (MD) simulations on size-selected particles employing a realistic analytic many-body potential based on quantum mechanical nanoparticle calculations have been performed to study the temperature-dependent properties and melting transitions of free Al n clusters and nanoparticles with n = 10-300 from 200 to 1700 K. By analyzing properties of the particles such as specific heat capacity (c), radius of gyration, volume, coefficient of thermal expansion (beta), and isothermal compressibility (kappa), we developed operational definitions of the solid, slush, and liquid states of metal clusters and nanoparticles. Applying the definitions, which are based on the temperature dependences of c, beta, and ln kappa, we determined the temperature domains of the solid, slush, and liquid states of the Al n particles. The results show that Al n clusters ( n or= 19, diameter of more than 1 nm) do have a melting transition and are in the liquid state above 900-1000 K. However, all aluminum nanoparticles have a wide temperature interval corresponding to the slush state in which the solid and liquid states coexist in equilibrium, unlike a bulk material where coexistence occurs only at a single temperature (for a given pressure). The commonly accepted operational marker of the melting temperature, namely, the peak position of c, is not unambiguous and not appropriate for characterizing the melting transition for aluminum particles with the exception of a few particle sizes that have a single sharp peak (as a function of temperature) in each of the three properties, c, beta, and ln kappa.     

Constrained fluid lambda-integration: constructing a reversible thermodynamic path between the solid and liquid state

A novel lambda-integration path is proposed for calculating the Gibbs free energy difference between any arbitrary solid and liquid state needed for the location of melting lines. This technique involves reversibly forcing a liquid state to a solid state across the phase transition along a nonphysical path, thermodynamically coupling the two states directly. The process eliminates the need for coupling to idealized reference states as is presently performed and hence simplifies the location of phase transitions for computer simulation systems. More specifically the path involves a three stage process, whereby, initially a liquid state is transformed to a weakly attractive fluid using linear lambda-integration scaling of the intermolecular potential. In the second stage, the resulting fluid is then constrained to the required solid configurational phase space via the insertion of a periodic lattice of 3D Gaussian wells. The final stage involves reversing to full strength the main intermolecular potential while gradually turning off the constraining 3D Gaussian lattice finally resulting in a stable (or metastable) solid state. Each stage was found to be completely reversible and the resulting change in free energy was thermodynamically integrable. The methodology is demonstrated and validated by calculating solid-liquid coexistence points using the new technique and comparing to those in present literature for the truncated and shifted Lennard-Jones system. The results are found to be in good agreement. The new method is not limited to melting phase transitions and is readily applicable to any simulation methodology, simulation cell size and/or intermolecular potential including ab initio methods.     

The salting-out effect and phase separation in aqueous solutions of tri-sodium citrate and 1-butyl-3-methylimidazolium bromide

The aim of this work is to obtain further evidence about the salting-out effect produced by the addition of tri-sodium citrate to aqueous solutions of water miscible ionic liquid 1-butyl-3-methylimidazolium bromide ([C₄mim][Br]) by evaluating the effect of tri-sodium citrate on the thermodynamic properties of aqueous solutions of this ionic liquid. Experimental measurements of density and sound velocity at different temperatures ranging from (288.15 to 308.15) K, the refractive index at 308.15 K and the liquid–liquid phase diagram at different temperatures ranging from (288.15 to 338.15) K for aqueous solutions containing 1-butyl-3-methylimidazolium bromide ([C₄mim][Br]) and tri-sodium citrate (Na₃Cit) are taken. The apparent molar volume of transfer of [C₄mim][Br] from water to aqueous solutions of Na₃Cit have positive values and it increases by increasing salt molality. Although at high IL molality, the apparent molar isentropic compressibility shows similar behaviour with that of the apparent molar volume. However at low concentrations of IL, the apparent molar isentropic compressibility of transfer of [C₄mim][Br] from water to aqueous solutions of Na₃Cit have negative values. The effects of temperature and the addition of Na₃Cit and [C₄mim][Br] on the liquid–liquid phase diagram of the investigated system have been studied. It was found that an increase in temperature caused the expansion of the one-phase region. The presence of Na₃Cit triggers a salting-out effect, leading to significant upward shifts of the liquid–liquid de-mixing temperatures of the system. The effect of temperature on the phase-forming ability in the system investigated has been studied based on a salting-out coefficient obtained from fitting the binodal values to a Setschenow-type equation for each temperature. Based on cloud point values, the energetics of the clouding process have been estimated and it was found that both of entropy and enthalpy are the driving forces for biphasic formation.     

Solubility of phosphonium ionic liquid in alcohols, benzene, and alkylbenzenes

The (solid + liquid) phase equilibria and (liquid + liquid) phase equilibria of binary mixtures containing quaternary phosphonium salt-tetrabutylphosphonium methanesulfonate and alcohols or alkylbenzenes were investigated. The systems {[(CH(3)CH(2)CH(2)CH(2))4P][CH(3)SO(3)] + 1-butanol, or 1-hexanol, 1-octanol, 1-decanol, or 1-dodecanol} and {[(CH(3)CH(2)CH(2)CH(2))4P][CH(3)SO(3)] + benzene, or toluene, ethylbenzene, or propylbenzene} have been measured by a dynamic method at a wide range of temperatures from 220 to 386 K. Solid-liquid equilibria with immiscibility in the liquid phase were detected with the aromatic hydrocarbons ethylbenzene and propylbenzene. The basic thermodynamic properties of pure ionic liquid--the melting point, enthalpy of fusion, enthalpy of solid-solid-phase transition, and glass transition--have been determined by differential scanning calorimetry. The experimental data of systems with alcohols were correlated by means of the UNIQUAC ASM and NRTL1 equations and of systems with alkylbenzenes with Wilson and NRTL equations utilizing parameters derived from the (solid + liquid) equilibrium. The root-mean-square deviations of the solubility temperatures for all calculated data are dependent upon the particular system and the particular equation used.     

A molecular dynamics study of the crystalline and liquid phases of pyridine

A molecular dynamics simulation study of the solid and liquid phases of pyridine (C₅H₅N) has been carried out based on an intermolecular potential parameterized to solid state data. Structural and thermodynamic properties, as well as selected time correlation functions have been calculated for both phases. The simulations show evidence of a strong coupling between translational and rotational motion, particularly in the liquid phase. Comparison with the experimental data for both phases reveals inadequacies in the potential model employed.     

The structure,thermodynamic properties and infrared spectra of liquid water and ice

Monte Carlo simulations are used, together with models of the intramolecular and intermolecular potential surfaces, to model liquid water and several phases of ice. Intramolecular relaxation makes important contributions to both thermodynamic and structural properties. A quantum local mode analysis of the Monte Carlo configurations is used to predict the density of states and infrared absorption intensities for the intramolecular bending and stretching vibrations. The large shifts from the gas phase OH stretch frequencies observed experimentally in the liquid and solid phases are due to anharmonic terms in the intramolecular surface rather than to harmonic intermolecular coupling. A significant contribution to observed changes in IR intensity on condensation arises from the large molecular polarisability.     

Molecular dynamics simulation of the energetic room-temperature ionic liquid, 1-hydroxyethyl-4-amino-1,2,4-triazolium nitrate (HEATN)

Molecular dynamics (MD) simulations have been performed to investigate the structure and dynamics of an energetic ionic liquid, 1-hydroxyethyl-4-amino-1,2,4-triazolium nitrate (HEATN). The generalized amber force field (GAFF) was used, and an electronically polarizable model was further developed in the spirit of our previous work (Yan, T.; Burnham, C. J.; Del Popolo, M. G.; Voth, G. A. J. Phys. Chem. B 2004, 108, 11877). In the process of simulated annealing from a liquid state at 475 K down to a glassy state at 175 K, the MD simulations identify a glass-transition temperature region at around 250-275 K, in agreement with experiment. The self-intermediate scattering functions show vanishing boson peaks in the supercooled region, indicating that HEATN may be a fragile glass former. The coupling/decoupling of translational and reorientational ion motion is also discussed, and various other physical properties of the liquid state are intensively studied at 400 K. A complex hydrogen bond network was revealed with the calculation of partial radial distribution functions. When compared to the similarly sized 1-ethyl-4-methyl-1,4-imidazolium nitrate ionic liquid, EMIM+/NO3-, a hydrogen bond network directly resulting in the poorer packing efficiency of ions is observed, which is responsible for the lower melting/glass-transition point. The structural properties of the liquid/vacuum interface shows that there is vanishing layering at the interface, in accordance with the poor ion packing. The effects of electronic polarization on the self-diffusion, viscosity, and surface tension of HEATN are found to be significant, in agreement with an earlier study on EMIM+/NO3- (Yan, T.; Burnham, C. J.; Del Popolo, M. G.; Voth, G. A. J. Phys. Chem. B 2004, 108, 11877).     

Bindungsverhältnisse in den kristallinen Phasen von Caesiumhydroxid und -deuterohydroxid,CsOH und CsOD

Bonding Conditions in Crystalline Phases of CsOH and CsOD The crystalline phases of CsOH and CsOD have been investigated between 23 K and the melting point by x-ray and neutron diffraction experiments as well as by calorimetric, IR, and Ramanspectroscopic measurements. For crystal data see ?Inhaltsübersicht”?. In the modification stable at low temperatures which resembles a TlI-type like structure the charge-asymmetric hydroxide ions are onedimensional hydrogen bonded. Certain forms of movements of the anions occur by warming up. In the high temperature NnC1-type structure the anions rotate quasi freely. The results of the various methods of investigation of the alkali metal hydroxides in the solid state are discussed with respect to the chemical bonding in these compounds.     

Far infrared spectra of liquid and solid benzene

The far infrared spectra of liquid and solid benzene are compared. The benzene was maintained at a temperature of 330 K and the applied pressure varied to alter the state. Liquid benzene is observed to possess a structure resembling that of the solid, close to the melting point.     

Coexistence of liquid and solid phases of Bmim-PF6 ionic liquid on mica surfaces at room temperature

In this communication, we report the coexistence of liquid and solid phases of room temperature ionic liquid (IL) [Bmim][PF6] on mica surfaces, observed by tapping mode atomic force microscopy (AFM) in air. Multilayers as well as drop-on-the-layer structures of the ionic liquid are revealed. The solid layers were very stable, and their orientations were affected by the mica lattice, indicating that the ionic liquid undergoes a template-induced ordering process. These results are helpful for advancing the understanding of interfacial structures of ILs on solid surfaces, the analogous structural patterns in both of its solid and liquid phases, and its heterogeneity.     

Solid–liquid equilibria for systems containing 1-butyl-3-methylimidazolium chloride

The solubility of 1-butyl-3-methylimidazolium chloride [C₄mim][Cl] in alcohols {ethanol, 1-butanol, 1-hexanol, 1-octanol, 1-decanol, 1-dodecanol, 2-butanol, 2-methyl-2-propanol (tert-butanol)} has been measured by a dynamic method from 270 K to the melting point of the ionic liquid or to the boiling point of the solvent. The melting point, enthalpy of fusion, and the temperature of the glass phase transition were determined by differential scanning calorimetry.The solubility data were correlated by means of the Wilson, UNIQUAC ASM and modified NRTL1 equations utilizing parameters derived from the solid–liquid equilibrium data. The root-mean-square deviations of the solubility temperatures for all calculated data were higher than 0.9 K and depended on the particular equation used.