Comparison of different water models from ambient to supercritical conditions: A Monte Carlo simulation and molecular Ornstein–Zernike study

Theoretical Study of Raman Intensities of p-Nitroaniline in Different Solvent Conditions by Using a Reference Interaction Site Model Self-Consistent Field Explicitly Including Constrained Spatial Electron Density Distribution

Raman spectroscopy is one of the most powerful tools to understand and characterize the states and structures of systems in several environments. To obtain highly accurate changes in Raman intensities of systems in solution, theoretical treatment, which can deal with not only the states and structures of systems but also the environment around molecules, proves to be significant. Hence, in this study, we developed the calculation of changes in Raman intensities of systems in different solvent conditions by using the reference interaction site model self-consistent field study explicitly including constrained spatial electron density distribution; this model is designed based on elements from both quantum mechanics and statistical mechanics. We showed that our calculation method could reproduce the changes in Raman intensities of p-nitroaniline (pNA) under different solvent conditions, including supercritical water, which has been observed in previous experimental studies. Based on the analysis of the calculation results, we observed that the ratio of the Raman intensity change of pNA in different solvent conditions is strongly correlated with the charge-transfer character of pNA.

Mechanisms of Efficient Desalination by a Two-Dimensional Porous Nanosheet Prepared via Bottom-Up Assembly of Cucurbit[6]urils

Many researchers have examined the desalination performance of various kinds of two-dimensional (2D) porous nanosheets prepared by top-down approaches such as forming pores on the plain based on molecular dynamics (MD) simulations. In contrast, it is rare to find MD simulations addressing the desalination performance of a 2D porous nanosheet prepared by bottom-up approaches. We investigated the desalination performance of a 2D porous nanosheet prepared by the assembly of cucurbit[6]uril (CB[6]) via MD simulation. The model 2D CB[6] nanosheet features CB[6] with the carbonyl-fringed portals of 3.9 Å and the interstitial space filled with hydrophobic linkers and dangling side chains. Our MD simulation demonstrated that the 2D porous CB[6] nanosheet possesses a 70 to 140 times higher water permeance than commercial reverse osmosis membranes while effectively preventing salt passage. The extremely high water permeance and perfect salt rejection stem from not only CB[6]’s nature (hydrophilicity, negative charge, and the right dimension for size exclusion) but also the hydrophobic and tightly filled interstitial space. We also double-checked that the extremely high water permeance was attributable to only CB[6]’s nature, not water leakage, by contrasting it with a 2D nanosheet comprising CB[6]-spermine complexes. Lastly, this paper provides a discussion on a better cucurbituril homologue to prepare a next-generation desalination membrane possessing great potential to such an extent to surpass the 2D porous CB[6] nanosheet based on quantum mechanics calculations.

Singular value decomposition analysis of the electron density changes occurring upon electrostatic polarization of water

In-depth elucidation of how molecules are electrically polarized would be one key factor for understanding the properties of those molecules under various thermodynamic and/or spatial conditions. Here this problem is tackled for the case of hydrogen-bonded water by conducting singular value decomposition of the electron density changes that occur upon electrostatic polarization. It is shown that all those electron density changes are approximately described as linear combinations of ten orthonormal basis “vectors”. One main component is the interatomic charge transfer through each OH bond, while some others are characterized as the atomic dipolar polarizations, meaning that both of these components are important for the electrostatic polarization of water. The interaction parameters that reasonably well reproduce the induced dipole moments are derived, which indicate the extent of mixing of the two components in electrostatic polarization.

The main features of the electron density changes that occur upon electrostatic polarization of water are elucidated by conducting singular value decomposition analysis of those changes.

Simulations of acid dissociation constants of polyprotic acids in near-critical and supercritical water

This paper reports a molecular dynamics study on the dissociation of sulfuric acid and phosphoric acid in near-critical and supercritical water. pK(a) is known to vary as the temperature and pressure vary, and this variation has important implications for corrosion in supercritical water reaction vessels. This work uses the SPC/E water model and solutes based upon DFT calculations to examine both structural and thermodynamic properties of the dissociation processes. An increase in solute-solvent ordering is observed for larger charges, and this also corresponds to a lower rate of diffusion. All dissociation reactions become less favored with increasing temperature except pK(a1) for sulfuric acid which becomes significantly more favored until 748 K.

TOWARDS A MOLECULAR THEORY FOR THE VAN DER WAALS–MAXWELL DESCRIPTION OF FLUID PHASE TRANSITIONS

We briefly review developments of theories for phase transitions of molecular fluids and mixtures, from semi-phenomenological approaches providing equations of state with adjustable parameters to first-principles microscopic methods qualitatively correct for a variety of molecular models with realistic interaction potentials. We further present the generalization of the van der Waals–Maxwell description of fluid phase diagrams to account for chemical specificities of polar molecular fluids, such as hydrogen bonding. Our theory uses the reference interaction site model (RISM) integral equation approach complemented with the new closure we have proposed (KH approximation), successful over a wide range of density from gas to liquid. The RISM/KH theory is applied to the known three-site models of water, methanol, and hydrogen fluoride. It qualitatively reproduces their vapor-liquid phase diagrams and the structure in the gas as well as liquid phases, including hydrogen bonding. Furthermore, phase transitions of water and methanol sorbed in nanoporous carbon aerogel are described by means of the replica generalization of the RISM approach we have developed. The changes as compared to the bulk fluids are in agreement with simulations and experiment. The RISM/KH theory is promising for description of phase transitions in various associating fluids, in particular, electrolyte as well as non-electrolyte solutions and ionic liquids.

Influence of Salts During Hydrothermal Biomass Gasification: The Role of the Catalysed Water-Gas Shift Reaction

Near-critical and supercritical water is an unusual reaction medium with extraordinary properties, varying with temperature and density. This opens the opportunity for new reactions and new technical processes. In the studies of synthesis reactions, total oxidation and biomass gasification interesting effects of salts are observed. These effects are assumed to be caused by complex formation, acidity/basicity or the presence of active hydrogen due to the water-gas shift reaction. The water-gas shift reaction is catalysed by salts and lead to the formation of hydrogen, which might react with other compounds.

This article focuses on the salt effects during hydrogen production by biomass gasification in supercritical water. This is a very interesting process to make use of “wet biomass”, which is up to now not applied in a technical process. Here salts influenced the main reactions pathways. Possible reasons are discussed.

Molecular modeling and dynamics studies with explicit inclusion of electronic polarizability. Theory and applications

A current emphasis in empirical force fields is on the development of potential functions that explicitly treat electronic polarizability. In the present article, the commonly used methodologies for modelling electronic polarization are presented along with an overview of selected application studies. Models presented include induced point-dipoles, classical Drude oscillators, and fluctuating charge methods. The theoretical background of each method is followed by an introduction to extended Langrangian integrators required for computationally tractable molecular dynamics simulations using polarizable force fields. The remainder of the review focuses on application studies using these methods. Emphasis is placed on water models, for which numerous examples exist, with a more thorough discussion presented on the recently published models associated with the Drude-based CHARMM and the AMOEBA force fields. The utility of polarizable models for the study of ion solvation is then presented followed by an overview of studies of small molecules (e.g. CCl(4), alkanes, etc) and macromolecule (proteins, nucleic acids and lipid bilayers) application studies. The review is written with the goal of providing a general overview of the current status of the field and to facilitate future application and developments.

Kirkwood-Buff integrals of aqueous alcohol binary mixtures

The Kirkwood-Buff integrals of some binary aqueous alcohol mixtures are computed from the available vapor pressure measurements and compared with previous results as well as small angle neutron scattering experiments. The emphasis of the present report is on accuracy of the results that can be achieved by these two different types of measurements. This seems to be needed, mainly in view of the discrepancies between the various published results, as shown herein. It is argued that agreement in peak positions is more important than that in magnitude. In general, very good agreement is obtained by both methods, and sources of disagreements are discussed. The issue of the computer simulations of aqueous systems and the problematics related to correlations, microheterogeneity, and consequently the Kirkwood-Buff integrals are equally discussed herein.

Percolation Transition in Supercritical Water: A Monte Carlo Simulation Study

Computer simulations of water have been performed on the canonical ensemble at 15 different molecular number densities, ranging from 0.006 to 0.018 A-3, along the supercritical isotherm of 700 K, in order to characterize the percolation transition in the system. It is found that the percolation transition occurs at a somewhat higher density than what is corresponding to the supercritical extension of the boiling line. We have shown that the fractal dimension of the largest cluster and the probability of finding a spanning cluster are the most appropriate properties for the location of the true percolation threshold. Thus, percolation transition occurs when the fractal dimension of the largest cluster reaches 2.53, and the probability of finding a cluster that spans the system in at least one dimension and in all the three dimensions reaches 0.97 and 0.65, respectively. On the other hand, the percolation threshold cannot be accurately located through the cluster size distribution, as it is distorted by appearance of clusters crossing the finite simulated system even far below the percolation threshold. The structure of the largest water cluster is dominated by a linear, chainlike arrangement, which does not change noticeably until the largest cluster becomes infinite.

Formation of mesoscopic water networks in aqueous systems

Formation of the macroscopically-infinite hydrogen-bonded water network in various aqueous systems occurs via 3D percolation transition when the probability of finding a spanning water cluster exceeds 95%. As a result, in a wide interval of water content below the percolation threshold, rarefied quasi-2D water networks span over the mesoscopic length scale. Formation and topology of spanning water networks, which affect various properties of aqueous systems, can be described within the framework of the percolation theory.

Spinodal of supercooled polarizable water

We develop a series of molecular dynamics computer simulations of liquid water, performed with a polarizable potential model, to calculate the spinodal line and the curve of maximum density inside the metastable supercooled region. After analyzing the structural properties, the liquid spinodal line is followed down to T=210K . A monotonic decrease is found in the explored region. The curve of maximum density bends on approaching the spinodal line. These results, in agreement with similar studies on nonpolarizable models of water, are consistent with the existence of a second critical point for water.

Solvation in supercritical water

Solvation in supercritical water under equilibrium and nonequilibrium conditions is studied via molecular dynamics simulations. The influence of solute charge distributions and solvent density on the solvation structures and dynamics is examined with a diatomic probe solute molecule. It is found that the solvation structure varies dramatically with the solute dipole moment, especially in low-density water, in accord with many previous studies on ion solvation. This electrostrictive effect has important consequences for solvation dynamics. In the case of a nonequilibrium solvent relaxation, if there are sufficiently many water molecules close to the solute at the outset of the relaxation, the solvent response measured as a dynamic Stokes shift is almost completely governed by inertial rotations of these water molecules. By contrast, in the opposite case of a low local solvent density near the solute, not only rotations but also translations of water molecules play an important role in solvent relaxation dynamics. The applicability of a linear response is found to be significantly restricted at low water densities.

Molecular dynamics simulations of liquid methanol and methanol–water mixtures with polarizable models

A polarizable model for simulation of liquid methanol, compatible with the COS/G2 water model, has been developed using the Charge-on-Spring (COS) technique. The model consists of three point charges, with one polarizable center on the oxygen atom. The Lennard-Jones parameters on the oxygen atom together with the molecular polarizability were varied to reproduce the experimental heat of vaporization and density of liquid methanol at ambient conditions. We examined the energies of various methanol dimers in the gas phase and compared them with values obtained from ab initio calculations. The model was then used to study the thermodynamic, dynamic, structural, and dielectric properties of liquid methanol as well as of a methanol-water mixture. A microscopic picture of the structure of pure liquid methanol and of the methanol-water mixture is provided. Good agreement was found between the results from our model simulations and available experimental and ab initio calculation data. In particular, the experimental dielectric permittivity of 32 could be reproduced, which had been shown to be difficult when using nonpolarizable models.

Structure of the Liquid−Vapor Interface of Water−Methanol Mixtures as Seen from Monte Carlo Simulations

Monte Carlo simulation of the vapor-liquid interface of water-methanol mixtures of different compositions, ranging from pure water to pure methanol, have been performed on the canonical (N, V, T) ensemble at 298 K. The analysis of the systems simulated has revealed that the interface is characterized by a double layer structure: methanol is strongly adsorbed at the vapor side of the interface, whereas this adsorption layer is followed at its liquid side by a depletion layer of methanol of lower concentration than in the bulk liquid phase of the system. The dominant feature of the interface has been found to be the adsorption layer in systems of methanol mole fractions below 0.2, and the depletion layer in systems of methanol mole fractions between 0.25 and 0.5. The orientation of the molecules located at the depletion layer is found to be already uncorrelated with the interface, whereas the methanol molecules of the adsorption layer prefer to align perpendicular to the interface, pointing straight toward the vapor phase by their methyl group. Although both the preference of the molecular plane for a perpendicular alignment with the interface and the preference of the methyl group for pointing straight to the vapor phase are found to be rather weak, the preference of the methyl group for pointing as straight toward the vapor phase as possible within the constraint imposed by the orientation of the molecular plane is found to be fairly strong. One of the two preferred orientations of the interfacial water molecules present in the neat system is found to disappear in the presence of methanol, because methanol molecules aligned in their preferred orientation can replace these water molecules in the hydrogen-bonding pattern of the interface.

Line of percolation in supercritical water

The percolation transition of the hydrogen-bonded clusters of molecules is investigated in supercritical water by Monte Carlo computer simulations. Simulations have been performed at four thermodynamic state points located above the supercritical extension of the vapor-liquid coexistence curve on the p-T phase diagram and at four state points located below this curve. It is found in a temperature range of a few hundred Kelvin that the extension of the vapor-liquid coexistence curve separates the supercritical thermodynamic states in which the water molecules form infinite hydrogen-bonded clusters from those in which the hydrogen-bonded clusters are isolated oligomers. However, the difference between the size of the hydrogen-bonded clusters at thermodynamic states located at the two sides of the extension of the coexistence curve is found to decrease with increasing temperature, and the present results suggest that this difference is likely to vanish at high enough temperatures.

Simulations of vapor water clusters at vapor–liquid equilibrium

The Gibbs-ensemble Monte Carlo methods based on the extended single point charge [H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, J. Phys. Chem. 91, 6269 (1987)] potential-energy surface have been used to study the clustering of vapor phase water under vapor-liquid equilibrium conditions between 300 and 600 K. It is seen that the number of clusters, as well as the cluster size, increase with temperature. This is primarily due to the increase in vapor density that accompanies the temperature increase at equilibrium. In addition, due to entropic effects, the percentage of clusters that have linear (or open) topologies increases with temperature and dominates over the minimum-energy cyclic topologies at the temperatures studied here. These results are insensitive to the number of molecules used in the simulations and the criterion used to define a water cluster.

Liquid–vapor and liquid–liquid phase equilibria of the Brodholt–Sampoli–Vallauri polarizable water model

Liquid-vapor and liquid-liquid phase equilibria of the polarizable Brodholt-Sampoli-Vallauri water model have been investigated by Gibbs ensemble Monte Carlo computer simulations. The coexisting liquid and vapor densities and energy of vaporization of the model is found to be in a reasonable agreement with experimental data in the entire temperature range of liquid-vapor coexistence. The critical temperature and density of the model are found to be 615 K and 0.278 gcm(3), respectively, close to the experimental values of 647.1 K and 0.322 gcm(3). In the supercooled state two distinct liquid-liquid coexistence regions are observed. The existence of liquid-liquid phase separation of a polarizable water model is demonstrated for the first time.

Charge-on-spring polarizable water models revisited: From water clusters to liquid water to ice

The properties of two improved versions of charge-on-spring (COS) polarizable water models (COS/G2 and COS/G3) that explicitly include nonadditive polarization effects are reported. In COS models, the polarization is represented via a self-consistently induced dipole moment consisting of a pair of separated charges. A previous polarizable water model (COS/B2), upon which the improved versions are based, was developed by Yu, Hansson, and van Gunsteren. To improve the COS/B2 model, which overestimated the dielectric permittivity, one additional virtual atomic site was used to reproduce the water monomer quadrupole moments besides the water monomer dipole moment in the gas phase. The molecular polarizability, residing on the virtual atomic site, and Lennard-Jones parameters for oxygen-oxygen interactions were varied to reproduce the experimental values for the heat of vaporization and the density of liquid water at room temperature and pressure. The improved models were used to study the properties of liquid water at various thermodynamic states as well as gaseous water clusters and ice. Overall, good agreement is obtained between simulated properties and those derived from experiments and ab initio calculations. The COS/G2 and COS/G3 models may serve as simple, classical, rigid, polarizable water models for the study of organic solutes and biopolymers. Due to its simplicity, COS type of polarization can straightforwardly be used to introduce explicit polarization into (bio)molecular force fields.

Structural and electronic evolution of the As(OH)3 molecule in high temperature aqueous solutions: An x-ray absorption investigation

The geometrical and electronic structure of the arsenious acid molecule As(OH)(3) in aqueous solutions has been investigated by x-ray absorption spectroscopy (XAS) within extended x-ray absorption spectroscopy (EXAFS) and x-ray absorption near edge structure (XANES), using realistic first-principle calculations in the latter case. This investigation was performed on aqueous solutions of arsenious acid from ambient to supercritical conditions (P = 250 and 600 bars, T <or= 500 degrees C) using a new optical cell. The analysis of the XAS spectra is consistent with (1) a constant As-O distance, (2) an opening of the O-As-O angles within the C(3V) pyramidal structure in the range 30-200 degrees C. This structural evolution comes along with a small decrease of the partial charges of the atoms in the As(OH)(3) molecule. The explanation invoked for both structural and electronic modifications observed is the weakening of the interactions, through hydrogen bonds, between the As(OH)(3) complex and water molecules. This is a fingerprint of the similar weakening of hydrogen bonding interactions in the solvent itself.