No Arabic abstract
We have investigated temperature trends of the microscopic structure of the SPC/E and TIP4P-Ew water models in terms of the pair distribution functions, coordination numbers, the average number of hydrogen bonds, the distribution of bonding states of a single molecule as well as the angular distribution of molecules by using the constant pressure molecular dynamics simulations. The evolution of the structure is put in correspondence with the dependence of water density on high temperatures down to the region of temperatures where the system becomes supercooled. It is shown that the fraction of molecules with three and four bonds determine the maximum density for both models. Moreover, the temperature dependence of the dielectric constant is obtained and analyzed.
Isothermal-isobaric molecular dynamics simulations have been performed to examine an ample set of properties of the model water-N,N-dimethylformamide (DMF) mixture as a function of composition. The SPC-E and TIP4P-Ew water models together with two united atom models for DMF [Chalaris M., Samios J., J. Chem. Phys., 2000, 112, 8581; Cordeiro J., Int. J. Quantum Chem., 1997, 65, 709] were used. Our principal analyses concern the behaviour of structural properties in terms of radial distribution functions, and the number of hydrogen bonds between molecules of different species as well as thermodynamic properties. Namely, we explore the density, excess mixing molar volume and enthalpy, the heat capacity and excess mixing heat capacity. Finally, the self-diffusion coefficients of species and the dielectric constant of the system are discussed. In addition, surface tension of water-DMF mixtures has been calculated and analyzed.
We have investigated thermodynamic and dynamic properties as well as the dielectric constant of water-metha-nol model mixtures in the entire range of composition by using constant pressure molecular dynamics simulations at ambient conditions. The SPC/E and TIP4P/Ew water models are used in combination with the OPLS united atom modelling for methanol. Changes of the average number of hydrogen bonds between particles of different species and of the fractions of differently bonded molecules are put in correspondence with the behavior of excess mixing volume and enthalpy, of self-diffusion coefficients and rotational relaxation times. From the detailed analyses of the results obtained in this work, we conclude that an improvement of the description of an ample set of properties of water-methanol mixtures can possibly be reached, if a more sophisticated, carefully parameterized, e.g., all atom, model for methanol is used. Moreover, exploration of parametrization of the methanol force field, with simultaneous application of different combination rules for methanol-water cross interactions, is required.
Isothermal-isobaric molecular dynamics simulations are used to examine the microscopic structure and other properties of a model solution consisting of NaCl salt dissolved in water-methanol mixture. The SPC/E water model and the united atom model for methanol are combined with the force field for ions by Dang [J. Amer. Chem. Soc., 1995, 117, 6954] to describe the entire system. Our principal focus is to study the effects of two variables, namely, the solvent composition and ion concentrations on the solutions density, on the structural properties, self-diffusion coefficients of the species and the dielectric constant. Moreover, we performed a detailed analysis of the first coordination numbers of the species. Trends of the behaviour of the average number of hydrogen bonds between solvent molecules are evaluated.
The properties of model solutions consisting of a solute --- single curcumin molecule in water, methanol and dimethyl sulfoxide solvents have been studied using molecular dynamics (MD) computer simulations in the isobaric-isothermal ensemble. The united atom OPLS force field (OPLS-UA) model for curcumin molecule proposed by us recently [J. Mol. Liq., 2016, 223, 707] in combination with the SPC/E water, and the OPLS-UA type models for methanol and dimethyl sulfoxide have been applied. We have described changes of the internal structure of the solute molecule induced by different solvent media in very detail. The pair distribution functions between particular fragments of a solute molecule with solvent particles have been analyzed. Statistical features of the hydrogen bonding between different species were explored. Finally, we have obtained a self-diffusion coefficient of curcumin molecules in three model solvents.
The ultraviolet (UV) photodissociation of amorphous water ice at different ice temperatures is investigated using molecular dynamics (MD) simulations and analytical potentials. Previous MD calculations of UV photodissociation of amorphous and crystalline water ice at 10 K [S. Andersson et al., J. Chem. Phys. 124, 064715 (2006)] revealed -for both types of ice- that H atom, OH, and H2O desorption are the most important processes after photoexcitation in the uppermost layers of the ice. Water desorption takes place either by direct desorption of recombined water, or when, after dissociation, an H atom transfers part of its kinetic energy to one of the surrounding water molecules which is thereby kicked out from the ice. We present results of MD simulations of UV photodissociation of amorphous ice at 10, 20, 30, and 90 K in order to analyze the effect of ice temperature on UV photodissociation processes. Desorption and trapping probabilities are calculated for photoexcitation of H2O in the top four monolayers and the main conclusions are in agreement with the 10 K results: desorption dominates in the top layers, while trapping occurs deeper in the ice. The hydrogen atom photodesorption probability does not depend on ice temperature, but OH and H2O photodesorption probabilities tend to increase slightly (~30%) with ice temperature. We have compared the total photodesorption probability (OH+H2O) with the experimental total photodesorption yield, and in both cases the probabilities rise smoothly with ice temperature. The experimental yield is on average 3.8 times larger than our theoretical results, which can be explained by the different time scales studied and the approximations in our model.