Structural, density, entropy, and diffusivity anomalies of the TIP4P/2005 model of water are mapped out over a wide range of densities and temperatures. The locus of temperatures of maximum density (TMD) for this model is very close to the experimental TMD locus for temperatures between 250 and 275 K. Four different water models (mTIP3P, TIP4P, TIP5P, and SPC/E) are compared with the TIP4P/2005 model in terms of their anomalous behavior. For all the water models, the density regimes for anomalous behavior are bounded by a low-density limit at around 0.85-0.90 g cm(-3) and a high-density limit at about 1.10-1.15 g cm(-3). The onset temperatures of the density anomaly in the various models show a much greater variation, ranging from 202 K for mTIP3P to 289 K for TIP5P. The order maps for the various water models are qualitatively very similar with the structurally anomalous regions almost superimposable in the q(tet)-τ plane. Comparison of the phase diagrams of water models with the region of liquid-state anomalies shows that the crystalline phases are much more sensitive to the choice of water models than the liquid state anomalies; for example, SPC/E and TIP4P/2005 show qualitatively similar liquid state anomalies but very different phase diagrams. The anomalies in the liquid in all the models occur at much lower pressures than those at which the melting line changes from negative to positive slope. The results in this study demonstrate several aspects of structure-entropy-diffusivity relationships of water models that can be compared with experiment and used to develop better atomistic and coarse-grained models for water.

Anomalous behavior of the excess entropy (S(e)) and the associated scaling relationship with diffusivity are compared in liquids with very different underlying interactions but similar water-like anomalies: water (SPC/E and TIP3P models), tetrahedral ionic melts (SiO(2) and BeF(2)), and a fluid with core-softened, two-scale ramp (2SRP) interactions. We demonstrate the presence of an excess entropy anomaly in the two water models. Using length and energy scales appropriate for onset of anomalous behavior, we show the density range of the excess entropy anomaly to be much narrower in water than in ionic melts or the 2SRP fluid. While the reduced diffusivities (D*) conform to the excess-entropy-scaling relation, D* = A exp(alphaS(e)) for all the systems (Rosenfeld, Y. Phys. Rev. A 1977, 15, 2545), the exponential scaling parameter, alpha, shows a small isochore dependence in the case of water. Replacing S(e) by pair correlation-based approximants accentuates the isochore dependence of the diffusivity scaling. Isochores with similar diffusivity-scaling parameters are shown to have the temperature dependence of the corresponding entropic contribution. The relationship between diffusivity, excess entropy, and pair correlation approximants to the excess entropy are very similar in all the tetrahedral liquids.

Molecular dynamics simulations of the Oeffner-Elliot model of germania (GeO(2)) are performed to identify nested regions of anomalous behavior in structural order, diffusivity, and pair entropy in the density-temperature plane, analogous to that seen in BeF(2), SiO(2), and H(2)O. The decreasing constraint of local tetrahedrality in GeO(2), compared to SiO(2) and BeF(2), substantially lowers the onset temperatures for anomalous behavior relative to the experimental melting temperatures (T(m)). Germania resembles water, more strongly than the ionic melts, in terms of temperatures for onset of anomalous behavior as well as in the order maps; for example, the structural anomaly sets in at 3.42T(m) in BeF(2), 3.09T(m) in SiO(2), 1.43T(m) in GeO(2), and 1.21T(m) in H(2)O. The detailed shapes of the anomalous regimes vary for different systems but the relative temperatures of onset for different anomalies are very similar in the different systems. The pair correlation entropy is shown to be a crucial and experimentally accessible quantity for relating structure, entropy, and diffusivity that could be potentially useful for a large class of inorganic ionic liquids.

Molecular dynamics simulations of water, liquid beryllium fluoride and silica melt are used to study the accuracy with which the entropy of ionic and molecular liquids can be estimated from atom-atom radial distribution function data. The pair correlation entropy is demonstrated to be sufficiently accurate that the density-temperature regime of anomalous behaviour as well as the strength of the entropy anomaly can be predicted reliably for both ionic melts as well as different rigid-body pair potentials for water. Errors in the total thermodynamic entropy for ionic melts due to the pair correlation approximation are of the order of 10% or less for most state points but can be significantly larger in the anomalous regime at very low temperatures. In the case of water, the rigid-body constraints result in larger errors in the pair correlation approximation, between 20 and 30%, for most state points. Comparison of the excess entropy, S e , of ionic melts with the pair correlation entropy, S 2 , shows that the temperature dependence of S e is well described by T −2/5 scaling across both the normal and anomalous regimes, unlike in the case of S 2 . The residual multiparticle entropy, ∆S = S e − S 2 , shows a strong negative correlation with tetrahedral order in the anomalous regime.

No abstract

Density functional theory (DFT) calculations were performed to study the mechanism of carbon dioxide (CO 2 ) reduction to carbon monoxide (CO) and methanol (CH 3 OH) on CeO 2 (110) surface. CO 2 dissociates to CO on interacting with the oxygen vacancy on reduced ceria surface.The oxygen atom heals the vacancy site and regenerates the stoichiometric surface via a redox mechanism with intrinsic activation and reaction energies of 259.2 kJ/mole and 238.6 kJ/mole respectively. Lateral interaction of oxygen vacancies were studied by the generation of two oxygen vacancies per unit of CeO 2 surface. Compared to a single isolated vacancy, the activation and reaction energies of CO 2 dissociation on a di-vacancy were approximately reduced to half of its value. Hydrogen atom co-adsorbed on the surface was observed to assist CO 2 dissociation by forming a carboxyl intermediate, CO 2 +H→COOH (∆E act = 39.0 kJ/mole, ∆H = -69.2 kJ/mole) which on subsequent dissociation produces CO via the redox mechanism. On hydrogenation, CO is likely to produce methanol. The energetics of CO hydrogenation to produce methanol showed exothermic steps with activation barriers comparable to the DFT calculations reported for Cu (111) and CeO 2-x /Cu(111) interface. While on the stoichiometric surface, COOH dissociation COOH→CO+OH (∆E act = 55.6 kJ/mole, ∆H =5.7 kJ/mole) is likely to be difficult as compared to rest of the elementary steps, whereas on the reduced surface the energetics of the same step were significantly lowered (∆E act = 47.4 kJ/mole, ∆H = -80.4 kJ/mole). In comparison, hydrogenation of methoxy, H 3 CO+H→H 3 COH, appears to be relatively difficult (∆E act = 58.7 kJ/mole) on the reduced surface.

The extent to which the presence of a biomolecular solute modifies the local energetics of water molecules, as measured by the tagged molecule potential energy (TPE), is examined using molecular dynamics simulations of the beta-hairpin of 2GB1 and the alpha-helix of deca-alanine in water. The CHARMM22 force field, in conjunction with the TIP3P solvent water model, is used for the peptides, with simulations of TIP3P and SPC/E water used as benchmarks for the behavior of bulk solvent. TIP3P water is shown to have significantly lower local tetrahedral order and higher binding energy than SPC/E at the same state point. The TIP3P and SPC/E water models show very similar dynamical correlations in the TPE fluctuations on frequency scales greater than 0.1 cm(-1). In addition, the two models show the same linear correlation between mean tetrahedral order and binding energy, suggesting that the relationship between choice of water models and simulated hydration behavior may involve a complex interplay of static and dynamic factors. The introduction of a peptide in water modifies the local TPE of water molecules as a function of distance from the biomolecular interface. There is an oscillatory variation in the TPE with distance from the peptide for water molecules lying outside a 3 A radius and extending to at least 10 A. These variations are of the order of 2-5% of the bulk TPE value and are anticorrelated with variations in local tetrahedral order in terms of locations of maxima and minima, which may be understood in terms of the relative contribution of van der Waals and Coulombic contributions to the TPE. The distance-dependent variations in local order and energetics are essentially the same for the beta-hairpin of 2GB1 as well as deca-alanine. Within a radius of 3 A, the perturbation of the solvent structure is very significant with local TPEs that are 10-15% lower than the bulk value. The chemical identity of side-chain residues and the secondary structure play an important role in determining residue-dependent variations in the TPEs. The variation in the residue-dependent tagged molecule potential energies is of the order of 3-5%, while the local residence times vary by a factor of approximately 5. The correlation of the local residence times with the local energetics within the innermost hydration layer is weak, though charged residues typically have low binding energies and large residence times.

Diffusivity, ionic conductivity, and viscosity of network-forming ionic melts are examined using molecular dynamics simulations of BeF2 and SiO2 melts. These tetrahedral, network-forming ionic melts are shown to possess a conductivity anomaly, in addition to waterlike viscosity and diffusivity anomalies, corresponding to a striking breakdown of the Nernst-Einstein relation. The contrasting scaling behavior of the different mobility measures with different structural contributions to the excess entropy is demonstrated.