SO2 Saturation of the Room Temperature Ionic Liquid [C(2)mim][NTf2] Much Reduces the Activation Energy for Diffusion

The physical effect of high concentrations of reversibly dissolved SO2 on [C(2)mim][NTf2] was examined using cyclic voltammetry, chronoamperometry, and ESR spectroscopy. Cyclic voltammetry of the oxidation of solutions of ferrocene, N,N,N',N'-tetramethyl-p-phenylenediamine (TMPD), and chloride in the room temperature ionic liquid (RTIL) 1-ethyl-3-methylimidazolium bis(trifluoromethanesufonyl)imide ([C(2)mim][NTf2]) reveals an increase in limiting current of each species corresponding to the addition of increasing concentrations of sulfur dioxide. Quantitative chronoamperometry reveals an increase in each species' diffusion coefficient with SO2 concentration. When chronoamperometric data were obtained for ferrocene in [C(2)mim][NTf2] at a range of temperatures, the translational diffusion activation energy (29.0 +/- 0.5 kJ mol(-1)) was found to be in good agreement with previous studies. Adding SO2 results in apparent near-activationless translational diffusion. A significant decrease in the activation energy of rotational diffusion with the SO2 saturation of a 2,2,6,6-tetramethyl-1-piperidinyloxyl (TEMPO) solution in [C(2)mim][NTf2] (29.9 +/- 2.0 to 7.7 +/- 5.3 kJ mol(-1)) was observed using electron spin resonance (ESR) spectroscopy. The reversible physical absorption Of SO2 by [C(2)mim][NTf2] should have no adverse effect on the ability of that ionic liquid to be employed as a solvent in an electrochemical gas sensor, and it is possible that the SO2-mediated reduction of RTIL viscosity could have intrinsic utility.

High energy resolution fluorescence detection XANES - an in situ method to study the interaction of adsorbed molecules with metal catalysts in the liquid phase

The change in the Pt electronic structure following the adsorption of an a,ß-unsaturated aldehyde and ketone was followed by in situ HERFD-XANES in the liquid phase. The resulting shift in the Pt Fermi energy is in good agreement with the molecule adsorption energy trends calculated by DFT and provides insight into the reaction selectivity.

Solid-liquid interfacial free energy out of equilibrium

The properties of the interface between solid and melt are key to solidification and melting, as the interfacial free energy introduces a kinetic barrier to phase transitions. This makes solidification happen below the melting temperature, in out-of-equilibrium conditions at which the interfacial free energy is ill defined. Here we draw a connection between the atomistic description of a diffuse solid-liquid interface and its thermodynamic characterization. This framework resolves the ambiguities in defining the solid-liquid interfacial free energy above and below the melting temperature. In addition, we introduce a simulation protocol that allows solid-liquid interfaces to be reversibly created and destroyed at conditions relevant for experiments. We directly evaluate the value of the interfacial free energy away from the melting point for a simple but realistic atomic potential, and find a more complex temperature dependence than the constant positive slope that has been generally assumed based on phenomenological considerations and that has been used to interpret experiments. This methodology could be easily extended to the study of other phase transitions, from condensation to precipitation. Our analysis can help reconcile the textbook picture of classical nucleation theory with the growing body of atomistic studies and mesoscale models of solidification.

A note on the prediction of liquid hold-up with the stratified roll wave regime for gas/liquidco-current flow in horizontal pipes.

This note presents a simple model for prediction of liquid hold-up in two-phase horizontal pipe flow for the stratified roll wave (St+RW) flow regime. Liquid hold-up data for horizontal two-phase pipe flow [1, 2, 3, 4, 5 and 6] exhibit a steady increase with liquid velocity and a more dramatic fall with increasing gas rate as shown by Hand et al. [7 and 8] for example. In addition the liquid hold-up is reported to show an additional variation with pipe diameter. Generally, if the initial liquid rate for the no-gas flow condition gives a liquid height below the pipe centre line, the flow patterns pass successively through the stratified (St), stratified ripple (St+R), stratified roll wave, film plus droplet (F+D) and finally the annular (A+D, A+RW, A+BTS) regimes as the gas rate is increased. Hand et al. [7 and 8] have given a detailed description of this progression in flow regime development and definitions of the patterns involved. Despite the fact that there are over one hundred models which have been developed to predict liquid hold-up, none have been shown to be universally useful, while only a handful have proven to be applicable to specific flow regimes [9, 10, 11 and 12]. One of the most intractable regimes to predict has been the stratified roll wave pattern where the liquid hold-up shows the most dramatic change with gas flow rate. It has been suggested that the momentum balance-type models, which give both hold-up and pressure drop prediction, can predict universally for all flow regimes but particularly in the case of the difficult stratified roll wave pattern. Donnelly [1] recently demonstrated that the momentum balance models experienced some difficulties in the prediction of this regime. Without going into lengthy details, these models differ in the assumed friction factor or shear stress on the surfaces within the pipe particularly at the liquid–gas interface. The Baker–Jardine model [13] when tested against the 0.0454 m i.d. data of Nguyen [2] exhibited a wide scatter for both liquid hold-up and pressure drop as shown in Fig. 1. The Andritsos–Hanratty model [14] gave better prediction of pressure drop but a wide scatter for liquid hold-up estimation (cf. Fig. 2) when tested against the 0.0935 m i.d. data of Hand [5]. The Spedding–Hand model [15], shown in Fig. 3 against the data of Hand [5], gave improved performance but was still unsatisfactory with the prediction of hold-up for stratified-type flows. The MARS model of Grolman [6] gave better prediction of hold-up (cf. Fig. 4) but deterioration in the estimation of pressure drop when tested against the data of Nguyen [2]. Thus no method is available that will accurately predict liquid hold-up across the whole range of flow patterns but particularly for the stratified plus roll wavy regime. The position is particularly unfortunate since the stratified-type regimes are perhaps the most predominant pattern found in multiphase lines.

A reactive force field simulation of liquid-liquid phase transitions in phosphorus

A force field model of phosphorus has been developed based on density functional (DF) computations and experimental results, covering low energy forms of local tetrahedral symmetry and more compact (simple cubic) structures that arise with increasing pressure. Rules tailored to DF data for the addition, deletion, and exchange of covalent bonds allow the system to adapt the bonding configuration to the thermodynamic state. Monte Carlo simulations in the N-P-T ensemble show that the molecular (P-4) liquid phase, stable at low pressure P and relatively low temperature T, transforms to a polymeric (gel) state on increasing either P or T. These phase changes are observed in recent experiments at similar thermodynamic conditions, as shown by the close agreement of computed and measured structure factors in the molecular and polymer phases. The polymeric phase obtained by increasing pressure has a dominant simple cubic character, while the polymer obtained by raising T at moderate pressure is tetrahedral. Comparison with DF results suggests that the latter is a semiconductor, while the cubic form is metallic. The simulations show that the T-induced polymerization is due to the entropy of the configuration of covalent bonds, as in the polymerization transition in sulfur. The transition observed with increasing P is the continuation at high T of the black P to arsenic (A17) structure observed in the solid state, and also corresponds to a semiconductor to metal transition. (C) 2004 American Institute of Physics.

Density functional and Monte Carlo studies of sulfur. II. Equilibrium polymerization of the liquid phase

The equilibrium polymerization of sulfur is investigated by Monte Carlo simulations. The potential energy model is based on density functional results for the cohesive energy, structural, and vibrational properties as well as reactivity of sulfur rings and chains [Part I, J. Chem. Phys. 118, 9257 (2003)]. Liquid samples of 2048 atoms are simulated at temperatures 450less than or equal toTless than or equal to850 K and P=0 starting from monodisperse S-8 molecular compositions. Thermally activated bond breaking processes lead to an equilibrium population of unsaturated atoms that can change the local pattern of covalent bonds and allow the system to approach equilibrium. The concentration of unsaturated atoms and the kinetics of bond interchanges is determined by the energy DeltaE(b) required to break a covalent bond. Equilibrium with respect to the bond distribution is achieved for 15less than or equal toDeltaE(b)less than or equal to21 kcal/mol over a wide temperature range (Tgreater than or equal to450 K), within which polymerization occurs readily, with entropy from the bond distribution overcompensating the increase in enthalpy. There is a maximum in the polymerized fraction at temperature T-max that depends on DeltaE(b). This fraction decreases at higher temperature because broken bonds and short chains proliferate and, for Tless than or equal toT(max), because entropy is less important than enthalpy. The molecular size distribution is described well by a Zimm-Schulz function, plus an isolated peak for S-8. Large molecules are almost exclusively open chains. Rings tend to have fewer than 24 atoms, and only S-8 is present in significant concentrations at all T. The T dependence of the density and the dependence of polymerization fraction and degree on DeltaE(b) give estimates of the polymerization temperature T-f=450+/-20 K. (C) 2003 American Institute of Physics.

Development of complex classical force fields through force matching to ab initio data: Application to a room-temperature ionic liquid

Recent experimental neutron diffraction data and ab initio molecular dynamics simulation of the ionic liquid dimethylimidazolium chloride ([dmim]Cl) have provided a structural description of the system at the molecular level. However, partial radial distribution functions calculated from the latter, when compared to previous classical simulation results, highlight some limitations in the structural description offered by force fieldbased simulations. With the availability of ab initio data it is possible to improve the classical description of [dmim]Cl by using the force matching approach, and the strategy for fitting complex force fields in their original functional form is discussed. A self-consistent optimization method for the generation of classical potentials of general functional form is presented and applied, and a force field that better reproduces the observed first principles forces is obtained. When used in simulation, it predicts structural data which reproduces more faithfully that observed in the ab initio studies. Some possible refinements to the technique, its application, and the general suitability of common potential energy functions used within many ionic liquid force fields are discussed.

Application of static charge transfer within an ionic-liquid force field and its effect on structure and dynamics

The effects of linear scaling of the atomic charges of a reference potential on the structure, dynamics, and energetics of the ionic liquid 1,3-dimethylimidazolium chloride are investigated. Diffusion coefficients that span over four orders of magnitude are observed between the original model and a scaled model in which the ionic charges are +/- 0.5 e. While the three-dimensional structure of the liquid is less affected, the partial radial distribution functions change markedly-with the positive result that for ionic charges of +/- 0.7 e, an excellent agreement is observed with ab initio molecular dynamics data. Cohesive energy densities calculated from these partial-charge models are also in better agreement with those calculated from the ab initio data. We postulate that ionic-liquid models in which the ionic charges are assumed to be +/- 1 e overestimate the intermolecular attractions between ions, which results in overstructuring, slow dynamics, and increased cohesive energy densities. The use of scaled-charge sets may be of benefit in the simulation of these systems-especially when looking at properties beyond liquid structure-thus providing on alternative to computationally expensive polarisable force fields.

Electrode kinetics and mechanism of iodine reduction in the room-temperature ionic liquid [C(4)mim][NTf2]

The fast electrochemical reduction of iodine in the RTIL 1-butyl-3-methylimidazolium bis(trifluoromethyl-sulfonyl)imide, [C(4)mim][NTf2], is reported and the kinetics and mechanism of the process elucidated. Two reduction peaks were observed. The first reduction peak is assigned to the process

Electrooxidation of the iodides [C(4)mim]I, LiI, NaI, KI, RbI, and CsI in the room temperature ionic liquid [C(4)mim][NTf2]

The electrochemical oxidation of 1-butyl-3-methylimidazolium iodide, [C(4)mim]I, has been investigated by cyclic voltammetry at a platinum microelectrode at varying concentrations in the RTIL 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide, [C(4)mim][NTf2]. Two oxidation peaks were observed. The first peak is assigned to the oxidation of iodide to triiodide, in an overall two-electron process: 3I(-)- 2e(-) -> I-3(-). At higher potentials, the electrogenerated triiodide oxidizes to iodine, in an overall one-electron process: I-3(-) - e(-) -> 3/2I(2). An average diffusion coefficient, D, for I- of 1.55 x 10(-11) m(2) s(-1) was obtained. A digital simulation program was used to simulate the voltammetric response, and kinetic parameters were successfully extracted. The parameters deduced from the simulation include D for I-, I-3(-), and I-2 and K-eq,K-2, the equilibrium constant for the reaction of iodide and iodine to form triiodide. Values for these parameters are of the same order as those previously published for the oxidation of Br- in the same RTIL [Allen et al. J. Electroanal. Chem. 2005, 575, 311]. Next, the cyclic voltammetry of five different inorganic iodide salts was studied by dissolving small amounts of the solid in [C(4)mim][NTf2]. Similar oxidation peaks were observed, revealing diffusion coefficients of ca. 0.55, 1.14, 1.23, 1.44, and 1.33 x 10(-11) m(2) s(-1) and solubilities of 714, 246, 54, 83, and 36 mM for LiI, NaI, KI, RbI, and CsI, respectively. The slightly smaller diffusion coefficients for the XI salts (compared to [C(4)mim]I) may indicate that I- is ion-paired with Li+, Na+, K+, Rb+, and Cs+ in the RTIL medium.

Glucose solvation by the ionic liquid 1,3-dimethylimidazolium chloride: A simulation study

Simulations of beta-glucose in the ionic liquid 1,3-dimethylimidazoliurn chloride have been performed in order to examine the solvation environment of the carbohydrate. Both single molecule and 1:5 glucose:ionic liquid (16.7 wt %) solutions are Studied, and the hydrogen bonding between sugar and solvent is examined. The primary solvation shell around the perimeter of the glucose ring consists predominantly of chloride anions which hydrogen bond to the hydroxyl groups. A small presence of the cation is also found, with the association Occurring through the weakly acidic hydrogen at the 2-position of the imidazolium ring interacting with the oxygen atoms of the sugar secondary hydroxyls. An average chloride coordination number of 4 is found around the glucose for both the single molecule and high concentration Simulations, despite the reduced chloride:glucose ratio in the latter case. In relation to the cation, the glucose molecules occupy positions above and below the plane of the imidazolium ring. Importantly, even at high glucose concentrations, no significant change in the anion-cation interactions and overall liquid structure of the ionic liquid is found, indicating that the glucose is readily accommodated by the solvent at this concentration. Dominant contributions to the sugar-ionic liquid interaction energy come from favorable hydrogen bonding (electrostatic) interactions between hydroxyls and chlorides, although a small favorable van der Waals energy contribution is also seen between the sugar and cations suggesting that the cation could be tailored in order to further improve the dissolution of glucose/cellulose in ionic liquid systems.

Electrochemical kinetics of Ag vertical bar Ag+ and TMPD vertical bar TMPD+center dot in the room-temperature ionic liquid [C(4)mpyrr][NTf2]; toward optimizing reference electrodes for voltammetry in RTILs

The voltammetry and kinetics of the Ag vertical bar Ag+ system (commonly used as a reference electrode material in both protic/aprotic and RTIL solvents) was studied in the room-temperature ionic liquid N-butyl-N-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide, [C(4)mpyrr][NTf2] on a 10 mu m diameter Pt electrode. For the three silver salts investigated (AgOTf, AgNTf2, and AgNO3, where OTf- = trifluoromethanesulfonate, NTf2- = bis(trifluoromethylsulfonyl)imide, and NO3- = nitrate), the voltammetry gave rise to a redox couple characteristic of a

A comparative electrochemical study of diffusion in room temperature ionic liquid solvents versus acetonitrile

Measurements on the diffusion coefficient of the neutral molecule N,N,N',N'-tetramethyl-para-phenylenediamine and the radical cation and dication generated by its one- and two-electron oxidation, respectively, are reported over the range 298-348 K in both acetonitrile and four room temperature ionic liquids (RTILs). Data were collected using single and double potential step chronamperometry at a gold disk electrode of micrometer dimension, and analysed via fitting to the appropriate analytical expression or, where necessary, to simulation. The variation of diffusion coefficient with temperature was found to occur in an Arrhenius-type manner for all combinations of solute and solvent. For a given ionic liquid, the diffusional activation energies of each species were not only closely equivalent to each other, but also to the RTIL's activation energy of viscous flow. In acetonitrile supported with 0.1 M tetrabutylammonium perchlorate, the ratio in diffusion coefficients of the radial cation and dication tot he neutral molecule were calculated as 0.89 +/- 0.05 and 0.51 +/- 0.03, respectively. In contrast, amongst the ionic liquids the same ratios were determined to be on average 0.53 +/- 0.04 and 0.33 +/- 0.03. The consequences of this dissimilarity are considered in terms of the modelling of voltammetric data gathered within ionic liquid solvents.

Measurements of the XUV mass absorption coefficient of an overdense liquid metal

The increase in the XUV mass absorption coefficient of liquid aluminium, produced by high-power-laser shock-compression, is measured using XUV laser radiography. At a photon energy of 63 eV a change in the mass absorption coefficient by up to a factor of similar to2.2 is determined at densities close to twice that of solid and electron temperatures of the order of 1 eV. Comparison with hydrodynamic simulations indicate that the absorption coefficient scales with density as rho (1.3 +/-0.2).