Water density, computer simulation

Adsorption energy, effect on density, computer simulation of structure and dynamics of water near metal surfaces, 34-36... 

The visualization of volumetric properties is more important in other scientific disciplines (e.g., computer tomography in medicine, or convection streams in geology). However, there are also some applications in chemistry (Figure 2-125d), among which only the distribution of water density in molecular dynamics simulations will be mentioned here. Computer visualization of this property is usually realized with two or three dimensional textures [203]. 

Tarek et al. [388] studied a system with some similarities to the work of Bocker et al. described earlier—a monolayer of n-tetradecyltrimethylammonium bromide. They also used explicit representations of the water molecules in a slab orientation, with the mono-layer on either side, in a molecular dynamics simulation. Their goal was to model more disordered, liquid states, so they chose two larger molecular areas, 0.45 and 0.67 nm molecule Density profiles normal to the interface were calculated and compared to neutron reflectivity data, with good agreement reported. The hydrocarbon chains were seen as highly disordered, and the diffusion was seen at both areas, with a factor of about 2.5 increase from the smaller molecular area to the larger area. They report no evidence of a tendency for the chains to aggregate into ordered islands, so perhaps this work can be seen as a realistic computer simulation depiction of a monolayer in an LE state. 

Data obtained on the electrode-oxygen distances and the reorientation of water in the inner layer are consistent with the results of computer simulations. The relative density profile for oxygen at the water/platinum interface shows two distinct maxima at 2.5 and 5.4 A, respectively. The first maximum corresponds to water directly bound to the surface whereas the less pronounced second maximum can be ascribed to the second layer of water linked through hydrogen bonds to the first layer. 

The calculations of the stmcture of water between charged flat walls show that the density profile becomes asymmetric and that there is enhanced structuring. This enhanced structuring is intimately connected with the possibility of a continuous phase transition in quasi two-dimensional systems, a subject of recent intense interest. ° Most of the molecular dynamics computer simulations on the effects of an external field have been carried out in an attempt to clarify the field-induced restructuring of water molecules at the metal surface, for which recent experimental data have become available. ... 

The molecular dynamic technique has been validated for water structures through comparison of calculated properties with experimental thermodynamic water data, such as the density maximum, the high heat capacity, and diffraction patterns (Stillinger and Rahman, 1974) as well as the hydrate infrared (vibrational) spectral data by Bertie and Jacobs (1977, 1982). With acceptable comparisons of many computed and experimental properties of water structures, there is little doubt that a substance similar to water has been simulated. 

The minimum information covers chemical formula, molecular weight, normal boiling point, freezing point, liquid density, water solubility and critical properties. Additional properties are enthalpies of phase transitions, heat capacity of ideal gas, heat capacity of liquid, viscosity and thermal conductivity of liquid. Computer simulation can estimate missing values. The use of graphs and tables of properties offers a wider view and is strongly recommended. 

Tunon, I., Martins-Costa M. T. C, Millot C., Ruiz-Lopez M. F. and Rivail J. L., A Coupled Density Functional-Molecular Mechanics Monte Carlo Simulation Method The Water Molecule in Liquid Water. J.Comput.Chem. (1996) 17 19-29. 

In this chapter, we have reviewed some of our own work on solvation properties in supercritical fluids using molecular dynamics computer simulations. We have presented the main aspects associated with the solvation structures of purine alkaloids in CO2 under different supercritical conditions and in the presence of ethanol as co-solvent, highlighting the phenomena of solvent density augmentation in the immediate neighborhood of the solute and the effects from the strong preferential solvation by the polar co-solvent. We have also presented a summary of our results for the structure and dynamics of supercritical water and ammonia, focusing on the dielectric behavior of supercritical water as functions of density and temperature and the behavior of excess solvated electrons in aqueous and non-aqueous associative environments. 

Our findings for rs and th may be compared with results of computer simulations for water. Values between 1 and 2 ps are stated for the average lifetime of a hydrogen bond by different authors (121-123), in satisfactory agreement with our experimental values. It is also interesting to compare with the frequency shift correlation function of the vibrational modes of water obtained from MD computations (124). Recently a slower component of this function with an exponential time constant of 0.8 ps was predicted for HDO in D20 at 300 K and a density of 1.1 g/cm3 (pressure %2 kbar). The existence of the slow component is a necessary prerequisite for the observation of spectral holes and the spectral relaxation time rs reported here. The faster component of the frequency shift correlation function with rc = 50 fs (124) represents rapid fluctuations that contribute to the spectral bandwidths of the spectral species and of the spectral holes. 

Can the Density Maximum of Water Be Found by Computer Simulation ... 

The results summarised here are for pure water at the temperature 25°C and the density 1.000 g cm , and are obtained by solving numerictdly the Ornstein-Zemike (OZ) equation for the pair correlation functions, using a closure that supplements the hypernet-ted chain (HNC) approximation with a bridge function. The bridge function is determined from computer simulations as described below. The numerical method for solving the OZ equation is described by Ichiye and Haymet and by Duh and Haymet. ... 

The complex cation ND4+ has been studied hy X-ray diffraction (49) and by the neutron isotopic difference method (31a), and results show that the coordination is weak, as might be expected for the low charge density. An MD computer simulation study using modified L-J potentials for the ion-water interactions gave good agreement with the ND results (75). 

Interaction potential models of water developed for computer simulations are typically fitted to the properties of the liquid phase. The most frequently used experimental data to be matched are the heat of vaporization (or the configurational internal energy), the structure at the level of pair correlations and the density. In the case of the most popular models tests have been carried out for further properties to check their performance Thanks to their classical, nonpolarizable character rigid planar models of water like the SPC/E and TIP4P are inexpensive to implement in computer simulations. In the following table we present some alternative parametrizations suggested in the literature recently. 

The Monte Carlo computer simulations by Kim and Landau [88] prove that also in the second and higher adsorbed layers the critical surface coverage is equal to one half. This means that in water adsorption on oxides, the density of the slab-like multilayer phase is below one half at a statistical coverage of three adsorbed layers. 

---