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.

Density functional theory study on the activation of molecular oxygen on a stepped gold surface in an aqueous environment: a new approach for simulating reactions in solution

The activation of oxygen molecules is an important issue in the gold-catalyzed partial oxidation of alcohols in aqueous solution. The complexity of the solution arising from a large number of solvent molecules makes it difficult to study the reaction in the system. In this work, O-2 activation on an Au catalyst is investigated using an effective approach to estimate the reaction barriers in the presence of solvent. Our calculations show that O-2 can be activated, undergoing OOH* in the presence of water molecules. The OOH* can readily be formed on Au(211) via four possible pathways with almost equivalent free energy barriers at the aqueous-solid interface: the direct or indirect activation of O-2 by surface hydrogen or the hydrolysis of O-2 following a Langmuir-Hinshelwood mechanism or an Eley-Rideal mechanism. Among them, the Eley-Rideal mechanism may be slightly more favorable due to the restriction of the low coverage of surface H on Au(211) in the other mechanisms. The results shed light on the importance of water molecules on the activation of oxygen in gold-catalyzed systems. Solvent is found to facilitate the oxygen activation process mainly by offering extra electrons and stabilizing the transition states. A correlation between the energy barrier and the negative charge of the reaction center is found. The activation barrier is substantially reduced by the aqueous environment, in which the first solvation shell plays the most important role in the barrier reduction. Our approach may be useful for estimating the reaction barriers in aqueous systems.

On the solvation of L-aspartic acid

We use molecular statics and dynamics to study the stability of L-aspartic acid both in vacuo and solvated by polar and non-polar molecules using density functional theory in the generalized gradient approximation. We find that structures stable in vacuo are unstable in aqueous solution and vice versa. From our simulations we are able to come to some conclusions about the mechanism of stabilisation of zwitterions by polar protic solvents, water and methanol.

Inner-shell photoexcitation of Fe XV and Fe XVI

The configuration-interaction method as implemented in the computer code CIV3 is used to determine energy levels, electric dipole radiative transition wavelengths, oscillator strengths and transition probabilities for inner-shell excitation of transitions in Fe XV and Fe XVI. Specifically, transitions are considered of the type 1s(2) 2s(2) 2p(6) 3s(2) -1s(2) 2s(2) 2p(5) 3l3l' 3l" (l, l' and l" = s,p or d) in FeXV and 1s(2) 2s(2) 2p(6) 3s- 1s(2) 2s(2) 2p(5) 3l3l' (l and l' = s,p or d) in FeXVI, using the relativistic Breit-Pauli approach. An assessment of the accuracy of the derived atomic data is performed.

Radiation damping effects on L-shell photoionization cross- sections of O-like Fe xix and N-like recombination rate coefficients for Fe xx

Photoionization cross-sections out of the fine-structure levels (2S(2)2p(4) P-3(2,0,1)) of the O-like Fe ion Fe XIX have been reinvestigated. Data for photoionization out of each of these finestructure levels have been obtained, where the calculations have been performed with and without the inclusion of radiation damping on the resonance structure in order to assess the importance of this process. Recombination rate coefficients are determined using the Milne relation, for the case of an electron recombining with N-like Fe ions (Fe XX) in the ground state to form O-like Fe (Fe XIX) existing in each of the fine- structure ground-state levels. Recombination rates are presented over a temperature range similar to 4.0 less than or equal to log T-e less than or equal to 7.0, of importance to the modelling of X-ray emission plasmas.

Ionization from the outer shell of Ar by proton impact

Recent results for proton-argon total ionization cross sections [Kirchner Phys. Rev. Lett. 79, 1658 (1997)] show large disagreement between theory and experiment for energies below 80 keV. To address this problem we have employed a recently developed theoretical method with a more pragmatic approach to the charge screening both in the initial and final channels. The target is considered as a one-electron atom and the interactions between this active electron and remaining target electrons are treated by a model potential including both short- and long-range effects. In the final channel the usual product of two continuum distorted wave functions each associated with a distinct electron-nucleus interaction is used. New results in the present calculation show good agreement in total cross sections for the energy range 10-300 keV with the measurement of Rudd [Rev. Mod. Phys. 57, 965 (1985)].

Fixed-node diffusion Monte Carlo computations for closed-shell jellium spheres

Fixed-node diffusion Monte Carlo computations are used to determine the ground state energy and electron density for jellium spheres with up to N = 106 electrons and background densities corresponding to the electron gas parameter 1 less than or equal to r(s)less than or equal to5.62. We analyze the density and size dependence of the surface energy, and we extrapolate our data to the thermodynamic limit. The results agree well with the predictions of density functional computations using the local density approximation. In the case of N = 20, we extend our computation to higher densities and identify a transition between atomic- and jelliumlike nodal structures occurring at the background density corresponding to r(s)=0.13. In this case the local density approximation is unable to reproduce the changes in the correlation energy due to the discontinuous transition in the ground state nodal structure. We discuss the relevance of our results for nonlocal approximations to density functional theory.

Solvation effects on equilibria: Triazoles and N-methyl piperidinol

We have performed calculations of the solvation effects on a number of equilibrium constants in water using a recently proposed hybrid quantum classical scheme in which the liquid environment is modelled using classical solvent molecules and the solute electronic structure is computed using modern quantum chemical methods. The liquid phase space is sampled from a fully classical simulation. We find that solvation effects on both triazole tautomeric equilibrium constants and piperidinol conformational equilibrium constants can be interpreted in terms of subtle differences in the local environment which can be seen in probability densities and radial distribution functions. Lower level calculations were performed for comparison and we conclude that the solvation thermodynamics can be predicted from a good classical model of solvent and solute molecules, but the implicit models that we tried are less successful.

Solvation Structure and Transport of Acidic Protons in Ionic Liquids: A First-principles Simulation Study

Ab initio simulations of a single molecule of HCl in liquid dimethyl imidazolium chloride [dmim][Cl] show that the acidic proton exists as a symmetric, linear ClHCl- species. Details of the solvation structure around this molecule are given. The proton-transfer process was investigated by applying a force along the antisymmetric stretch coordinate until the molecule broke. Changes in the free energy and local solvation structure during this process were investigated. In the reaction mechanism identified, a free chloride approaches the proton from the side. As the original ClHCl- distorts and the incoming chloride forms a new bond to the proton, one of the original chlorine atoms is expelled and a new linear molecule is formed.