Search

This paper presents a detailed study of a water adlayer adsorbed on Pt(111) and Rh(111) surfaces using periodic density functional theory methods. The interaction between the metal surface and the water molecules is assessed from molecular dynamics simulation data and single point electronic structure calculations of selected configurations. It is argued that the electron bands around the Fermi level of the metal substrate extend over the water adlayer. As a consequence in the presence of the water layer the surface as a whole still maintains its metallic conductivity-a result of a crucial importance for understanding the process of electron transfer through the water/metal interface and electrochemical reactions in particular. Our results also indicate that there exists a weak bond between the hydrogen of the water and the Rh metal atoms as opposed to the widespread (classical) models based on purely repulsive interaction. This suggests that the commonly used classical interactions potentials adopted for large scale molecular dynamics simulations of water/metal interfaces may need revision. Two adsorption models of water on transition metals with the OH bonds pointing towards or away of the surface are also examined. It is shown that due to the very close values of their adsorption energies one should consider the real structure of water on the surface as a mixture of these simple "up" and "down" models. A model for the structure of the adsorbed water layer on Rh(111) is proposed in terms of statistical averages from molecular dynamics simulations.

We report results from molecular dynamics simulations of the freezing transition of TIP5P water molecules confined between two parallel plates under the influence of a homogeneous external electric field, with magnitude of 5 V/nm, along the lateral direction. For water confined to a thickness of a trilayer we find two different phases of ice at a temperature of T=280 K. The transformation between the two, proton-ordered, ice phases is found to be a strong first-order transition. The low-density ice phase is built from hexagonal rings parallel to the confining walls and corresponds to the structure of cubic ice. The high-density ice phase has an in-plane rhombic symmetry of the oxygen atoms and larger distortion of hydrogen bond angles. The short-range order of the two ice phases is the same as the local structure of the two bilayer phases of liquid water found recently in the absence of an electric field [J. Chem. Phys. 119, 1694 (2003)]. These high- and low-density phases of water differ in local ordering at the level of the second shell of nearest neighbors. The results reported in this paper, show a close similarity between the local structure of the liquid phase and the short-range order of the corresponding solid phase. This similarity might be enhanced in water due to the deep attractive well characterizing hydrogen bond interactions. We also investigate the low-density ice phase confined to a thickness of 4, 5, and 8 molecular layers under the influence of an electric field at T=300 K. In general, we find that the degree of ordering decreases as the distance between the two confining walls increases. (C) 2004 American Institute of Physics.

The fifth-order two-dimensional (2D) Raman signals have been calculated from the equilibrium and nonequilibrium (finite field) molecular dynamics simulations. The equilibrium method evaluates response functions with equilibrium trajectories, while the nonequilibrium method calculates a molecular polarizability from nonequilibrium trajectories for different pulse configurations and sequences. In this paper, we introduce an efficient algorithm which hybridizes the existing two methods to avoid the time-consuming calculations of the stability matrices which are inherent in the equilibrium method. Using nonequilibrium trajectories for a single laser excitation, we are able to dramatically simplify the sampling process. With this approach, the 2D Raman signals for liquid xenon, carbon disulfide, water, acetonitrile, and formamide are calculated and discussed. Intensities of 2D Raman signals are also estimated and the peak strength of formamide is found to be only five times smaller than that of carbon disulfide.

An examination of the finite system size effects on the calculated interfacial properties of the molten salt KI is described, as an extension of the work described in the preceding paper [J. Chem. Phys. 115, 8603 (2001)]. It is shown that, for system sizes beyond ∼1000 ions, the finite size effects are accounted for by the predictions of capillary wave (CW) theory, as in van der Waals liquids. Hence, there do not appear to be additional long-range effects arising from the Coulombic interactions. Techniques used to avoid truncation effects in the Coulombic and dispersion interactions appear to have eliminated system size dependent artefacts from these sources. The scaling behavior predicted from CW theory is used to predict the bulk surface tension for the chosen interaction potential, which is then compared with experimental values, and very good agreement found. Examination of the interfacial layers show that the local structure remains strongly charge ordered, but the surface is highly roughened by clus...

In this paper we introduce a three-dimensional version of the Mercedes-Benz model to describe water molecules. In this model van der Waals interactions and hydrogen bonds are given explicitly through a Lennard-Jones potential and a Gaussian orientation-dependent terms, respectively. At low temperature the model freezes forming Ice-I and it reproduces the main peaks of the experimental radial distribution function of water. In addition to these structural properties, the model also captures the thermodynamical anomalies of water: The anomalous density profile, the negative thermal expansivity, the large heat capacity, and the minimum in the isothermal compressibility.

Molecular dynamics simulation has been performed on a liquid film that is sheared in between solid surfaces. As a shear is given to the liquid film, a Couette-like flow is generated in the liquid and energy conversion occurs from the macroscopic flow to the thermal energy, which is discharged back to the solid walls. In such a way, momentum and thermal energy fluxes are present simultaneously. And all these thermal and fluid phenomena take place in highly nonequilibrium state where thermal energy is not distributed equally to each degree of freedom of molecular motion in the vicinities of the solid-liquid interface. In the present paper, platinum and water are employed as solid and liquid, respectively. First, the structure and orientation of water molecules in the vicinities of the solid surfaces are analyzed and how these structure and orientation are influenced by the shear is considered. Based on this result, momentum and thermal energy transfer in the vicinities of and at the solid-liquid interfaces are investigated in detail. Results are compared with those of our previous study, in which monatomic and diatomic molecules are employed as liquid.