Theory for the three-dimensional Mercedes-Benz model of water

Isothermal-isobaric algorithm to study the effects of rotational degrees of freedom-Benz water model

We have developed isothermal-isobaric algorithm for non-equilibrium Monte Carlo simulations. As first we have shown that the new method correctly predict density by comparing it to the density determined in canonical Monte Carlo simulations through the virial pressure. The new method was then used to study the effect of translational and rotational degrees of freedom on the structural and thermodynamic properties of the simple Mercedes-Benz water model. By holding one of the temperatures constant and varying the other one, we investigated how the position of the density maxima changes. We have observed that upon increase of rotational temperature the fluid become more Lennard-Jones like and the density maxima disappears.

A new one-site coarse-grained model for water: Bottom-up many-body projected water (BUMPer). I. General theory and model

Water is undoubtedly one of the most important molecules for a variety of chemical and physical systems, and constructing precise yet effective coarse-grained (CG) water models has been a high priority for computer simulations. To recapitulate important local correlations in the CG water model, explicit higher-order interactions are often included. However, the advantages of coarse-graining may then be offset by the larger computational cost in the model parameterization and simulation execution. To leverage both the computational efficiency of the CG simulation and the inclusion of higher-order interactions, we propose a new statistical mechanical theory that effectively projects many-body interactions onto pairwise basis sets. The many-body projection theory presented in this work shares similar physics from liquid state theory, providing an efficient approach to account for higher-order interactions within the reduced model. We apply this theory to project the widely used Stillinger-Weber three-body interaction onto a pairwise (two-body) interaction for water. Based on the projected interaction with the correct long-range behavior, we denote the new CG water model as the Bottom-Up Many-Body Projected Water (BUMPer) model, where the resultant CG interaction corresponds to a prior model, the iteratively force-matched model. Unlike other pairwise CG models, BUMPer provides high-fidelity recapitulation of pair correlation functions and three-body distributions, as well as N-body correlation functions. BUMPer extensively improves upon the existing bottom-up CG water models by extending the accuracy and applicability of such models while maintaining a reduced computational cost.

Hierarchy of anomalies in the two-dimensional Mercedes-Benz model of water

We investigate by Monte Carlo simulations density, diffusion, and structural anomalies of the simple two-dimensional Mercedes-Benz (MB) model of water, which is a very simple toy model for explaining the origin of water properties. MB water molecules are modeled as two-dimensional Lennard-Jones disks, with three orientation-dependent hydrogen-bonding arms, arranged as in the MB logo. The model is in a way also a variance of silica-like models. Beside the known thermodynamic anomaly for the model we also found diffusion and structural anomalies and map out the cascade of density, structural, pair entropy, and diffusivity anomalies for MB model. The orientational order parameters with three and six-fold symmetry were determined and maximum for each one observed. The anomalies occur in hierarchy order, which is a slight variation of the hierarchy order in real water. The diffusion anomaly region is the innermost in the hierarchy while for water it is the density anomaly region.

Modelling water with simple Mercedes-Benz models

The structures and properties of biomolecules like proteins, nucleic acids, and membranes depend on water. Water is also very important in industry. Overall, water is unusual substance with more than 70 anomalous properties. The understanding of water is advancing significantly due to theoretical and computational modeling. There are different kind of models, models with fine-scale properties and increasing structural detail with increasing computational expense and simple models which focus on global properties of water like thermodynamics, phase diagram and are less computational expensive. Simplified models give a better understanding of water in ways that complement more complex models. Here, we review a simple model, the two dimensional Mercedes-Benz (MB) model of water. We present results by Monte Carlo simulations for anomalies and phase diagram and application of various theoretical methods.

Hierarchical lattice models of hydrogen-bond networks in water

We develop a graph-based model of the hydrogen-bond network in water, with a view toward quantitatively modeling the molecular-level correlational structure of the network. The networks formed are studied by the constructing the model on two infinite-dimensional lattices. Our models are built bottom up, based on microscopic information coming from atomistic simulations, and we show that the predictions of the model are consistent with known results from ab initio simulations of liquid water. We show that simple entropic models can predict the correlations and clustering of local-coordination defects around tetrahedral waters observed in the atomistic simulations. We also find that orientational correlations between bonds are longer ranged than density correlations, determine the directional correlations within closed loops, and show that the patterns of water wires within these structures are also consistent with previous atomistic simulations. Our models show the existence of density and compressibility anomalies, as seen in the real liquid, and the phase diagram of these models is consistent with the singularity-free scenario previously proposed by Sastry and coworkers [Phys. Rev. E 53, 6144 (1996)1063-651X10.1103/PhysRevE.53.6144].

Liquid part of the phase diagram and percolation line for two-dimensional Mercedes-Benz water

Monte Carlo simulations and Wertheim's thermodynamic perturbation theory (TPT) are used to predict the phase diagram and percolation curve for the simple two-dimensional Mercedes-Benz (MB) model of water. The MB model of water is quite popular for explaining water properties, but the phase diagram has not been reported till now. In the MB model, water molecules are modeled as two-dimensional Lennard-Jones disks, with three orientation-dependent hydrogen-bonding arms, arranged as in the MB logo. The liquid part of the phase space is explored using grand canonical Monte Carlo simulations and two versions of Wertheim's TPT for associative fluids, which have been used before to predict the properties of the simple MB model. We find that the theory reproduces well the physical properties of hot water but is less successful at capturing the more structured hydrogen bonding that occurs in cold water. In addition to reporting the phase diagram and percolation curve of the model, it is shown that the improved TPT predicts the phase diagram rather well, while the standard one predicts a phase transition at lower temperatures. For the percolation line, both versions have problems predicting the correct position of the line at high temperatures.

Analytical theory of the hydrophobic effect of solutes in water

We develop an analytical statistical-mechanical model for hydrophobic solvation in water. In this three-dimensional Mercedes-Benz-like model, two neighboring waters have three possible interaction states: a radial van der Waals interaction, a tetrahedral orientation-dependent hydrogen-bonding interaction, or no interaction. Nonpolar solutes are modeled as van der Waals particles of different radii. The model is sufficiently simple that we can calculate the partition function and thermal and volumetric properties of solvation versus temperature, pressure, and solute radius. Predictions are in good agreement with results of Monte Carlo simulations. And their trends agree with experiments on hydrophobic solute insertion. The theory shows that first-shell waters are more highly structured than bulk waters, because of hydrogen bonding, and that that structure melts out faster with temperature than it does in bulk waters. Because the theory is analytical, it can explore a broad range of solvation properties and anomalies of water, at minimal computational expense.

Integral equation and thermodynamic perturbation theory for a two-dimensional model of chain-forming fluid

In this paper we applied analytical theories for the two dimensional chain-forming fluid. Wertheims thermodynamic perturbation theory (TPT) and integral equation theory (IET) for associative liquids were used to study thermodynamical and structural properties of the chain-forming model. The model has polymerizing points at arbitrary position from center of the particles. Calculated analytical results were tested against corresponding results obtained by Monte Carlo computer simulations to check the accuracy of the theories. The theories are accurate for the different positions of patches of the model at all values of the temperature and density studied. The IET's pair correlation functions of the model agree well with computer simulations. Both TPT and IET are in good agreement with the Monte Carlo values of the energy, chemical potential and ratios of free, once and twice bonded particles.

Explicit-water theory for the salt-specific effects and Hofmeister series in protein solutions

Effects of addition of salts on stability of aqueous protein solutions are studied theoretically and the results are compared with experimental data. In our approach, all the interacting species, proteins, ions, and water molecules, are accounted for explicitly. Water molecules are modeled as hard spheres with four off-center attractive square-well sites. These sites serve to bind either another water or to solvate the ions or protein charges. The ions are represented as charged hard spheres, and decorated by attractive sites to allow solvation. Spherical proteins simultaneously possess positive and negative groups, represented by charged hard spheres, attached to the surface of the protein. The attractive square-well sites, mimicking the protein-protein van der Waals interaction, are located on the surface of the protein. To obtain numerical results, we utilized the energy route of Wertheim's associative mean spherical approximation. From measurable properties, we choose to calculate the second virial coefficient B2, which is closely related to the tendency of proteins to aggregate and eventually crystalize. Calculations are in agreement with experimental trends: (i) For low concentration of added salt, the alkali halide salts follow the inverse Hofmeister series. (ii) At higher concentration of added salt, the trend is reversed. (iii) When cations are varied, the salts follow the direct Hofmeister series. (iv) In contrast to the colloidal theories, our approach correctly predicts the non-monotonic behavior of B2 upon addition of salts. (v) With respect to anions, the theory predicts for the B2 values to follow different sequences below and above the iso-ionic point, as also confirmed experimentally. (vi) A semi-quantitative agreement between measured and calculated values for the second virial coefficient, as functions of pH of solution and added salt type and concentration, is obtained.

Liquid-liquid critical point in a simple analytical model of water

A statistical model for a simple three-dimensional Mercedes-Benz model of water was used to study phase diagrams. This model on a simple level describes the thermal and volumetric properties of waterlike molecules. A molecule is presented as a soft sphere with four directions in which hydrogen bonds can be formed. Two neighboring waters can interact through a van der Waals interaction or an orientation-dependent hydrogen-bonding interaction. For pure water, we explored properties such as molar volume, density, heat capacity, thermal expansion coefficient, and isothermal compressibility and found that the volumetric and thermal properties follow the same trends with temperature as in real water and are in good general agreement with Monte Carlo simulations. The model exhibits also two critical points for liquid-gas transition and transition between low-density and high-density fluid. Coexistence curves and a Widom line for the maximum and minimum in thermal expansion coefficient divides the phase space of the model into three parts: in one part we have gas region, in the second a high-density liquid, and the third region contains low-density liquid.

Protein Aggregation and Molecular Crowding: Perspectives From Multiscale Simulations

Cells are extremely crowded environments, thus the use of diluted salted aqueous solutions containing a single protein is too simplistic to mimic the real situation. Macromolecular crowding might affect protein structure, folding, shape, conformational stability, binding of small molecules, enzymatic activity, interactions with cognate biomolecules, and pathological aggregation. The latter phenomenon typically leads to the formation of amyloid fibrils that are linked to several lethal neurodegenerative diseases, but that can also play a functional role in certain organisms. The majority of molecular simulations performed before the last few years were conducted in diluted solutions and were restricted both in the timescales and in the system dimensions by the available computational resources. In recent years, several computational solutions were developed to get close to physiological conditions. In this review we summarize the main computational techniques used to tackle the issue of protein aggregation both in a diluted and in a crowded environment.

Recent developments in the theoretical, simulational, and experimental studies of the role of water hydrogen bonding in hydrophobic phenomena

Hydrophobic effects (hydrophobic hydration and hydrophobic interaction) constitute an important element of a wide variety of phenomena relevant to biological, physical, chemical, environmental, engineering, and pharmaceutical sciences, such as the immiscibility of oil and water, self-assembly of amphiphiles leading to micelle and membrane formation, folding and stability and unfolding of the native structure of a biologically active protein, gating of ion channels, wetting, froth floatation, and adhesion. On the other hand, the hydrogen bonding ability of water plays a major (if not crucial) role in hydrophobic phenomena. We present a review of most important and relatively recent experimental, simulational, and theoretical research on hydrophobic phenomena in various systems. With a particular interest we survey investigations clarifying the role of water hydrogen bonding therein, because it has been the main object of our own recent research. We have developed a probabilistic hydrogen bond (PHB) model that allows one to obtain an analytic expression for the number of bonds per water molecule as a function of its distance to a hydrophobe, hydrophobe radius, and temperature. Knowing that function, one can explicitly identify a water hydrogen bond contribution to the external potential whereto a water molecule is subjected near a hydrophobe. Combining the PHB model with the classical density functional theory (DFT), one can examine the contribution of water hydrogen bonding to the temperature and lengthscale effects on the hydration of particles and on their solvent-mediated interactions over the entire low-to-high temperature and small-to-large lengthscale ranges. We applied the combined DFT/PHB model to study a variety of hydrophobic phenomena such as (liquid) water in contact with a hydrophobic plate, solvation of spherical solutes of various radii in associated and non-associated liquids at various temperatures, the solvent-mediated interaction of spherical solutes and its temperature dependence, interaction of C60 fullerenes in water, temperature effect on the evaporation lengthscale of water confined between two hydrophobes, temperature dependence of the effective width of the solute-solvent transition layer and average density therein. These applications demonstrated that the DFT/PHB model can serve as a valuable tool in studying hydrophobic phenomena because it constitutes a balanced combination of simplicity, accuracy, and detail. The predictions of the combined DFT/PHB approach for the solvent density profiles and thermodynamic aspects of hydrophobic phenomena are generally in good agreement with experiments and simulations. For example, it predicts the small-to-large crossover lengthscale of its mechanism to be approximately in the range from 1nm to 4nm, and decreasing with increasing temperature. It also suggests that, in terms of the average fluid density in the solute-solvent transition layer, the transition layer for small hydrophobes (of radii ≲2 nm) becomes enriched with rather than depleted of fluid when both the solvent-solute affinity and hb-energy alteration ratio become large enough. The boundary values of these parameters, needed for the depletion-to-enrichment crossover, are predicted to decrease with increasing temperature.

Protein aggregation in salt solutions

Protein aggregation is broadly important in diseases and in formulations of biological drugs. Here, we develop a theoretical model for reversible protein-protein aggregation in salt solutions. We treat proteins as hard spheres having square-well-energy binding sites, using Wertheim's thermodynamic perturbation theory. The necessary condition required for such modeling to be realistic is that proteins in solution during the experiment remain in their compact form. Within this limitation our model gives accurate liquid-liquid coexistence curves for lysozyme and γ IIIa-crystallin solutions in respective buffers. It provides good fits to the cloud-point curves of lysozyme in buffer-salt mixtures as a function of the type and concentration of salt. It than predicts full coexistence curves, osmotic compressibilities, and second virial coefficients under such conditions. This treatment may also be relevant to protein crystallization.

A molecular-based approach to the thermodynamics of aqueous solutions: binary mixture of water and carbon dioxide

A simple model and theory of molecular fluids is applied to a binary mixture of water and carbon dioxide. An approach based on the perturbation theory is followed using a reference system of so-called pseudo-hard bodies for water and hard triatomics for carbon dioxide. Pseudo-hard bodies bear the traits of the non-additive nature of association supplementing the common excluded volume effect. The reference term is parametrized using Monte Carlo simulation data on the compressibility factor. After adding a simple mean-field term to the reference equation, fluid phase equilibria are qualitatively reproduced.

Existence of a liquid-liquid phase transition in methanol

A simple model is constructed to study the phase diagram and thermodynamic properties of methanol, which is described as a dimer of an apolar sphere mimicking the methyl group and a sphere with core-softened potential as the hydroxyl group. Performing classical Monte Carlo simulations, we obtained the phase diagram, showing a second critical point between two different liquid phases. Evaluating systems with a different number of particles, we extrapolate to infinite size in accordance with Ising universality class to obtain bulk values for critical temperature, pressure, and density. Strong evidence that the structure of the liquid changes upon transition from high- to low-density phase was provided. From the experimentally determined hydrogen bond strength and length in methanol and water, we propose where the second critical point of methanol should be.

Thermodynamics and the hydrophobic effect in a core-softened model and comparison with experiments

A simple and computationally inexpensive core-softened model, originally proposed by Franzese [G. Franzese, J. Mol. Liq. 136, 267 (2007)], was adopted to show that it exhibits properties of waterlike fluid and hydrophobic effect. The potential used between particles is spherically symmetric with two characteristic lengths. Thermodynamics of nonpolar solvation were modeled as an insertion of a modified Lennard-Jones particle. It was investigated how the anomalous predictions of the model as well as the nonpolar solvation compare with the experimental data for water anomalies and the temperature dependence of noble gases hydration. It was shown that the model qualitatively follows the same trends as water. The model is able to reproduce waterlike anomalous properties (density maximum, heat capacity minimum, isothermal compressibility, etc.) and hydrophobic effect (minimum solubility for nonpolar solutes near ambient conditions, increased solubility of larger noble gases, etc.). It is argued that the model yields similar results as more complex and computationally expensive models.

The hydrophobic effect in a simple isotropic water-like model: Monte Carlo study

Using Monte Carlo computer simulations, we show that a simple isotropic water-like model with two characteristic lengths can reproduce the hydrophobic effect and the solvation properties of small and large non-polar solutes. Influence of temperature, pressure, and solute size on the thermodynamic properties of apolar solute solvation in a water model was systematically studied, showing two different solvation regimes. Small particles can fit into the cavities around the solvent particles, inducing additional order in the system and lowering the overall entropy. Large particles force the solvent to disrupt their network, increasing the entropy of the system. At low temperatures, the ordering effect of small solutes is very pronounced. Above the cross-over temperature, which strongly depends on the solute size, the entropy change becomes strictly positive. Pressure dependence was also investigated, showing a "cross-over pressure" where the entropy and enthalpy of solvation are the lowest. These results suggest two fundamentally different solvation mechanisms, as observed experimentally in water and computationally in various water-like models.

Core-softened fluids as a model for water and the hydrophobic effect

An interaction model with core-softened potential in three dimensions was studied by Monte Carlo computer simulations and integral equation theory. We investigated the possibility that a fluid with a core-softened potential can reproduce anomalies found experimentally in liquid water, such as the density anomaly, the minimum in the isothermal compressibility as a function of temperature, and others. Critical points of the fluid were also determined. We provided additional arguments that the old notion, postulating that only angular-dependent interactions result in density anomaly, is incorrect. We showed that potential with two characteristic distances is sufficient for the system to exhibit water-like behavior and anomalies, including the famous density maximum. We also found that this model can properly describe the hydrophobic effect.

The application of the thermodynamic perturbation theory to study the hydrophobic hydration

The thermodynamic perturbation theory was tested against newly obtained Monte Carlo computer simulations to describe the major features of the hydrophobic effect in a simple 3D-Mercedes-Benz water model: the temperature and hydrophobe size dependence on entropy, enthalpy, and free energy of transfer of a simple hydrophobic solute into water. An excellent agreement was obtained between the theoretical and simulation results. Further, the thermodynamic perturbation theory qualitatively correctly (with respect to the experimental data) describes the solvation thermodynamics under conditions where the simulation results are difficult to obtain with good enough accuracy, e.g., at high pressures.

Ice polyamorphism in the minimal Mercedes-Benz model of water

We investigate ice polyamorphism in the context of the two-dimensional Mercedes-Benz model of water. We find a first-order phase transition between a crystalline phase and a high-density amorphous phase. Furthermore, we find a reversible transformation between two amorphous structures of high and low density; however, we find this to be a continuous and not an abrupt transition, as the low-density amorphous phase does not show structural stability. We discuss the origin of this behavior and its implications with regard to the minimal generic modeling of polyamorphism.

Coarse-Grained Molecular Models of Water: A Review

Coarse-grained (CG) models have proven to be very effective tools in the study of phenomena or systems that involve large time- and length-scales. By decreasing the degrees of freedom in the system and using softer interactions than seen in atomistic models, larger timesteps can be used and much longer simulation times can be studied. CG simulations are widely used to study systems of biological importance that are beyond the reach of atomistic simulation, necessitating a computationally efficient and accurate CG model for water. In this review, we discuss the methods used for developing CG water models and the relative advantages and disadvantages of the resulting models. In general, CG water models differ with regards to how many waters each CG group or bead represents, whether analytical or tabular potentials have been used to describe the interactions, and how the model incorporates electrostatic interactions. Finally, how the models are parameterized depends on their application, so, while some are fitted to experimental properties such as surface tension and density, others are fitted to radial distribution functions extracted from atomistic simulations.

Analytical model for three-dimensional Mercedes-Benz water molecules

We developed a statistical model which describes the thermal and volumetric properties of water-like molecules. A molecule is presented as a three-dimensional sphere with four hydrogen-bonding arms. Each water molecule interacts with its neighboring waters through a van der Waals interaction and an orientation-dependent hydrogen-bonding interaction. This model, which is largely analytical, is a variant of a model developed before for a two-dimensional Mercedes-Benz model of water. We explored properties such as molar volume, density, heat capacity, thermal expansion coefficient, and isothermal compressibility as a function of temperature and pressure. We found that the volumetric and thermal properties follow the same trends with temperature as in real water and are in good general agreement with Monte Carlo simulations, including the density anomaly, the minimum in the isothermal compressibility, and the decreased number of hydrogen bonds upon increasing the temperature.

Simple model of hydrophobic hydration

Water is an unusual liquid in its solvation properties. Here, we model the process of transferring a nonpolar solute into water. Our goal was to capture the physical balance between water's hydrogen bonding and van der Waals interactions in a model that is simple enough to be nearly analytical and not heavily computational. We develop a 2-dimensional Mercedes-Benz-like model of water with which we compute the free energy, enthalpy, entropy, and the heat capacity of transfer as a function of temperature, pressure, and solute size. As validation, we find that this model gives the same trends as Monte Carlo simulations of the underlying 2D model and gives qualitative agreement with experiments. The advantages of this model are that it gives simple insights and that computational time is negligible. It may provide a useful starting point for developing more efficient and more realistic 3D models of aqueous solvation.

Water and aqueous solutions: simple non-speculative model approach

Different ways of molecular modeling of water are analyzed and their similarities and differences identified. An up-to-date summary of achievements of a general approach to common rigid site-site interaction models of molecular fluids applied to water and aqueous solutions is then presented and discussed. The method is based on considering only a short-range part of a total realistic potential (such as SPC/E or TIPxP) which determines the structure of water (and fluids in general). A simplification of the interactions at short intermolecular separations leads then to simple models, called primitive models. Quite accurate results in an analytic form for the thermodynamic properties of the models are obtained using the thermodynamic perturbation theory. It is shown that the properly constructed primitive models reproduce, qualitatively, anomalies of pure water and basic characteristics of hydrophobic hydration. The concept of an extended excluded volume, based on pseudo-hard bodies, is introduced and exemplified by considering the partial molar volume of apolar solutes. Potential future development towards a theory of water based on the primitive models as a reference with the long-range contributions added as a perturbation is discussed.

A statistical mechanical theory for a two-dimensional model of water

We develop a statistical mechanical model for the thermal and volumetric properties of waterlike fluids. Each water molecule is a two-dimensional disk with three hydrogen-bonding arms. Each water interacts with neighboring waters through a van der Waals interaction and an orientation-dependent hydrogen-bonding interaction. This model, which is largely analytical, is a variant of the Truskett and Dill (TD) treatment of the "Mercedes-Benz" (MB) model. The present model gives better predictions than TD for hydrogen-bond populations in liquid water by distinguishing strong cooperative hydrogen bonds from weaker ones. We explore properties versus temperature T and pressure p. We find that the volumetric and thermal properties follow the same trends with T as real water and are in good general agreement with Monte Carlo simulations of MB water, including the density anomaly, the minimum in the isothermal compressibility, and the decreased number of hydrogen bonds for increasing temperature. The model reproduces that pressure squeezes out water's heat capacity and leads to a negative thermal expansion coefficient at low temperatures. In terms of water structuring, the variance in hydrogen-bonding angles increases with both T and p, while the variance in water density increases with T but decreases with p. Hydrogen bonding is an energy storage mechanism that leads to water's large heat capacity (for its size) and to the fragility in its cagelike structures, which are easily melted by temperature and pressure to a more van der Waals-like liquid state.