Universal tight binding model for chemical reactions in solution and at surfaces. II. Water

A revised water model intended for use in condensed phase simulations in the framework of the self consistent polarizable ion tight binding theory is constructed. The model is applied to water monomer, dimer, hexamers, ice, and liquid, where it demonstrates good agreement with theoretical results obtained by more accurate methods, such as DFT and CCSD(T), and with experiment. In particular, the temperature dependence of the self diffusion coefficient in liquid water predicted by the model, closely reproduces experimental curves in the temperature interval between 230 K and 350 K. In addition, and in contrast to standard DFT, the model properly orders the relative densities of liquid water and ice. A notable, but inevitable, shortcoming of the model is underestimation of the static dielectric constant by a factor of two. We demonstrate that the description of inter and intramolecular forces embodied in the tight binding approximation in quantum mechanics leads to a number of valuable insights which can be missing from ab initio quantum chemistry and classical force fields. These include a discussion of the origin of the enhanced molecular electric dipole moment in the condensed phases, and a detailed explanation for the increase of coordination number in liquid water as a function of temperature and compared with ice-leading to insights into the anomalous expansion on freezing. The theory holds out the prospect of an understanding of the currently unexplained density maximum of water near the freezing point. 

Universal tight binding model for chemical reactions in solution and at surfaces. III. Stoichiometric and reduced surfaces of titania and the adsorption of water

We demonstrate a model for stoichiometric and reduced titanium dioxide intended for use in molecular dynamics and other atomistic simulations and based in the polarizable ion tight binding theory. This extends the model introduced in two previous papers from molecular and liquid applications into the solid state, thus completing the task of providing a comprehensive and unified scheme for studying chemical reactions, particularly aimed at problems in catalysis and electrochemistry. As before, experimental results are given priority over theoretical ones in selecting targets for model fitting, for which we used crystal parameters and band gaps of titania bulk polymorphs, rutile and anatase. The model is applied to six low index titania surfaces, with and without oxygen vacancies and adsorbed water molecules, both in dissociated and non-dissociated states. Finally, we present the results of molecular dynamics simulation of an anatase cluster with a number of adsorbed water molecules and discuss the role of edge and corner atoms of the cluster. (C) 2014 AIP Publishing LLC.

Universal tight binding model for chemical reactions in solution and at surfaces. I. Organic molecules

As is now well established, a first order expansion of the Hohenberg-Kohn total energy density functional about a trial input density, namely, the Harris-Foulkes functional, can be used to rationalize a non self consistent tight binding model. If the expansion is taken to second order then the energy and electron density matrix need to be calculated self consistently and from this functional one can derive a charge self consistent tight binding theory. In this paper we have used this to describe a polarizable ion tight binding model which has the benefit of treating charge transfer in point multipoles. This admits a ready description of ionic polarizability and crystal field splitting. It is necessary in constructing such a model to find a number of parameters that mimic their more exact counterparts in the density functional theory. We describe in detail how this is done using a combination of intuition, exact analytical fitting, and a genetic optimization algorithm. Having obtained model parameters we show that this constitutes a transferable scheme that can be applied rather universally to small and medium sized organic molecules. We have shown that the model gives a good account of static structural and dynamic vibrational properties of a library of molecules, and finally we demonstrate the model's capability by showing a real time simulation of an enolization reaction in aqueous solution. In two subsequent papers, we show that the model is a great deal more general in that it will describe solvents and solid substrates and that therefore we have created a self consistent quantum mechanical scheme that may be applied to simulations in heterogeneous catalysis.

Dissolution and processing of cellulose using ionic liquids, cellulose solution, and regenerating cellulose.

Cellulose is dissolved in an ionic liq. without derivatization, and is regenerated in a range of structural forms without requiring the use of harmful or volatile org. solvents. Cellulose soly. and the soln. properties can be controlled by the selection of the ionic liq. constituents, with small cations and halide or pseudohalide anions favoring soln.; dissoln. can be aided by irradn. An ionic liq., [C4mim]Cl, proved to be the best for dissolving cellulose. [on SciFinder(R)]

Arsenic-glutathione complexes - Their stability in solution and during separation by different HPLC modes

Complexes of arsenic compounds and glutathione are believed to play an essential part in the metabolism and transport of inorganic arsenic and its methylated species. Up to now, the evidence of their presence is mostly indirect. We studied the stability and Chromatographic behaviour of glutathione complexes with trivalent arsenic: i.e. As^III(GS)₃, MA ^III(GS)₂ and DMA^III(GS) under different conditions. Standard ion chromatography using PRP X-100 and carbonate or formic acid buffer disintegrated the complexes, while all three complexes are stable and separable by reversed phase chromatography (0.1% formic acid/acetonitrile gradient). As^III(GS)₃ and MA^III(GS)₂ were more stable than DMA^III(GS), which even under optimal conditions tended to degrade on the column at 25 °C. Chromatography at 6 °C can retain the integrity of the samples. These results shed more light on the interpretation of a vast number of previously published arsenic speciation studies, which have used Chromatographic separation techniques with the assumption that the integrity of the arsenic species is guaranteed. © The Royal Society of Chemistry 2004.

Correlation between the photocatalysed oxidation of methylene blue in solution and the reduction of resazurin in a photocatalyst activity indicator ink (Rz Paii)

Determining the electronic structure and chemical potentials of molecules in solution

A simulation scheme is proposed for determining the excess chemical potential of a substance in solution. First, a Monte Carlo simulation is performed with classical models for solute and solvent molecules. A representative sample of these configurations is then used in a hybrid quantum/classical (QM/MM) calculation, where the solute is treated quantum-mechanically, and the average electronic structure is used to construct an improved classical model. This procedure is iterated to self-consistency in the classical model, which in practice is attained in one or two steps, depending on the quality of the initial guess. The excess free energy of the molecule within the QM/MM approach is determined relative to the classical model using thermodynamic perturbation theory with a cumulant expansion. The procedure provides a method of constructing classical point charge models appropriate for the solution and gives a measure of the importance of solvent fluctuations.

Boron in copper: A perfect misfit in the bulk and cohesion enhancer at a grain boundary

Using first principles electronic structure methods, we calculate the effects of boron impurities in bulk copper and at surfaces and grain boundaries. We find that boron segregation to the Sigma5(310)[001] grain boundary should strengthen the boundary up to 1.5 ML coverage (15.24 at./nm2). The maximal effect is observed at 0.5 ML and corresponds to boron atoms filling exclusively grain boundary interstices. In copper bulk, B causes significant distortion both in interstitial and regular lattice sites, for which boron atoms are either too big or too small. The distortion is compensated to a large extent when the interstitial and substitutional boron combine together to form a strongly bound dumbbell. Our prediction is that bound boron impurities should appear in a sizable proportion if not dominate in most experimental conditions. A large discrepancy between calculated heats of solution and experimental terminal solubility of B in Cu is found, indicating either a significant failure of the density functional approach or, more likely, strongly overestimated solubility limits in the existing B-Cu phase diagram.

Mutual Solubilities of Water and the [Cnmim][Tf2N] Hydrophobic Ionic Liquids

Ionic liquids (ILs) have recently garnered increased attention because of their potential environmental benefits as "green" replacements over conventional volatile organic solvents. While ILs cannot significantly volatilize and contribute to air pollution, even the most hydrophobic ones present some miscibility with water posing environmental risks to the aquatic ecosystems. Thus, the knowledge of ILs toxicity and their water solubility must be assessed before an accurate judgment of their environmental benefits and prior to their industrial applications. In this work, the mutual solubilities for [C2-C8mim][Tf2N] (n-alkyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide) and water between 288.15 and 318.15 K at atmospheric pressure were measured. Although these are among the most hydrophobic ionic liquids known, the solubility of water in these compounds is surprisingly large, ranging from 0.17 to 0.36 in mole fraction, while the solubility of these ILs in water is much lower ranging from 3.2 × 10-5 to 1.1 × 10-3 in mole fraction, in the temperature and pressure conditions studied. From the experimental data, the molar thermodynamic functions of solution and solvation such as Gibbs energy, enthalpy, and entropy at infinite dilution were estimated, showing that the solubility of these ILs in water is entropically driven. The predictive capability of COSMO-RS, a model based on unimolecular quantum chemistry calculations, was evaluated for the description of the binary systems investigated providing an acceptable agreement between the model predictions and the experimental data both with the temperature dependence and with the ILs structural variations.

Assessment of the solubility and bioaccessibility of arsenic in realgar wine using a simulated gastrointestinal system

Consumption of arsenic (As) wine is a traditional activity during the classic Chinese festival of Duanwu, colloquially known worldwide as the Dragon Boat Day. Arsenic wine is drunk on the morning of the fifth day of the fifth lunar calendar month to commemorate the death of Qu Yuan, a famed Chinese poet who drowned himself in protest of a corrupt government, and to protect against ill fortune. Although realgar minerals are characteristically composed of sparingly soluble tetra-arsenic tetra-sulfides (As(4)S(4)), purity does vary with up to 10% of As being present as non-sulfur bound species, such as arsenate (As(v)) and arsenite (As(III)). Despite, the renewed interest in As speciation and the bioaccessibility of the active As components in realgar based Chinese medicines, little is known about the safety surrounding the cultural practice of drinking As wine. In a series of experiments the speciation and solubility of As in a range of wines were investigated. Furthermore, a simulated gastrointestinal system was employed to predict the impact of digestive processes on As bioavailability. The predominant soluble As species found in all the wines were As(III) and As(v). Based on typical As wine recipes employing 0.1 g realgar mL(-1) wine, the concentration of dissolved As ranged from ca. 100 to 400 mg L(-1) depending on the ethanol content of the preparation: with the As solubility found to be higher in wines with a lower proportion of ethanol. Based on a common 100 mL measure of wine with a concentration of 400 mg As L(-1), the amount of soluble As would equate to around half of the acute minimal lethal dose for adults. This is likely an underestimate of the bioaccessible concentration, as a three-fold increase in bioaccessibility could be observed in the intestinal phase based on the results from the stimulated gastrointestinal system. (C) 2011 Elsevier B.V. All rights reserved.