Thermodynamic anomalies in a lattice model of water

We investigate a lattice-fluid model of water, defined on a three-dimensional body centered cubic lattice. Model molecules possess a tetrahedral symmetry, with four equivalent bonding arms, aiming to mimic the formation of hydrogen bonds. The model is similar to the one proposed by Roberts and Debenedetti [J. Chem. Phys. 105, 658 (1996)], simplified in that no distinction between bond "donors" and "acceptors" is imposed. Bond formation depends both on orientation and local density. In the ground state, we show that two different ordered (ice) phases are allowed. At finite temperature, we analyze homogeneous phases only, working out phase diagram, response functions, the temperature of maximum density locus, and the Kauzmann line. We make use of a generalized first-order approximation on a tetrahedral cluster. In the liquid phase, the model exhibits several anomalous properties observed in real water. In the low temperature region (supercooled liquid), there are evidences of a second critical point and, for some range of parameter values, this scenario is compatible with the existence of a reentrant spinodal.     

Thermodynamic and dynamic anomalies in a one-dimensional lattice model of liquid water

We investigate the occurrence of waterlike thermodynamic and dynamic anomalous behavior in a one dimensional lattice gas model. The system thermodynamics is obtained using the transfer matrix technique and anomalies on density and thermodynamic response functions are found. When the hydrogen bond (molecules separated by holes) is more attractive than the van der Waals interaction (molecules in contact) a transition between two fluid structures is found at null temperature and high pressure. This transition is analogous to a 'critical point' and intimately connects the anomalies in density and in thermodynamic response functions. Monte Carlo simulations were performed in the neighborhood of this transition and used to calculate the self diffusion constant, which increases with density as in liquid water.     

Thermodynamic properties of lattice hard-sphere models

Thermodynamic properties of several lattice hard-sphere models were obtained from grand canonical histogram- reweighting Monte Carlo simulations. Sphere centers occupy positions on a simple cubic lattice of unit spacing and exclude neighboring sites up to a distance sigma. The nearestneighbor exclusion model, sigma = radical2, was previously found to have a second-order transition. Models with integer values of sigma = 1 or 2 do not have any transitions. Models with sigma = radical3 and sigma = 3 have weak first-order fluid-solid transitions while those with sigma = 2 radical2, 2 radical3, and 3 radical2 have strong fluid-solid transitions. Pressure, chemical potential, and density are reported for all models and compared to the results for the continuum, theoretical predictions, and prior simulations when available.     

Thermodynamic properties of the Cu-Au system using a face-centered-cubic lattice model with a renormalized potential

A Monte Carlo simulation is carried out to study thermodynamic properties of Cu-Au alloys using a face-centered-cubic (fcc) lattice-gas model. To obtain quantitatively accurate results, a Finnis-Sinclair-type potential, which has been widely used for molecular dynamics (MD) simulations, is employed. To overcome some shortcomings of lattice-gas models such as neglecting vibrational entropy, the potential is mapped onto the fcc lattice using the renormalization technique. The renormalized potential gives an improved Cu-Au phase diagram compared to the original MD potential applied directly on the lattice.     

Thermodynamic stability and kinetic foldability of a lattice protein model

By using serial mutations, i.e., a residue replaced by 19 kinds of naturally occurring residues, the stability of native conformation and folding behavior of mutated sequences are studied. The 3 x 3 x 3 lattice protein model with two kinds of interaction potentials between the residues, namely the original Miyazawa and Jernigan (MJ) potentials and the modified MJ potentials (MMJ), is used. Effects of various sites in the mutated sequences on the stability and foldability are characterized through the Z-score and the folding time. It is found that the sites can be divided into three types, namely the hydrophobic-type (H-type), the hydrophilic-type (P-type) and the neutral-type (N-type). These three types of sites relate to the hydrophobic core, the hydrophilic surface and the parts between them. The stability of the native conformation for the serial mutated sequences increases (or decreases) as the increasing in the hydrophobicity of the mutated residues for the H-type sites (or the P-type sites), while varies randomly for the N-type sites. However, the foldability of the mutated sequences is not always consistent with the thermodynamic stability, and their relationship depends on the site types. Since the hydrophobic tendency of the MJ potentials is strong, the ratio between the number of the H-type sites and the number of the P-type sites is found to be 1:2. Differently, for the MJJ potentials it is found that such a ratio is about 1:1 which is relevant to that of real proteins. This suggests that the modification of the MJ potentials is rational in the aspect of thermodynamic stability. The folding of model proteins with the MMJ potentials is fast. However, the relationship between the foldability and the thermodynamic stability of the mutated sequences is complex.     

Thermodynamic equations of state from molecular solvation

A general expression of the canonical partition function for mixture fluids in terms of solvation free energy is presented. Following the same approach as set forth in generalized van der Waals theory, we show that the physical assumptions made in existing thermodynamic models from the perspective of molecular solvation can be prevailed. For example, the temperature dependence on the coordination number, i.e. the number of solvent molecules surrounding the solute, has an impact not only on the temperature dependence of the solvation free energy but on the non-linear solvent reaction field response as well. More importantly, the new formulation provides a unified scheme for deriving the commonly used equations of state (EOS) and newly developed COSMO-type liquid activity coefficient models. We show that from this new formulation it is possible to develop a new class of thermodynamic models that behave like existing EOS and liquid models in the low and high-density limits, respectively.     

Individual degrees of freedom and the solvation properties of water

Using molecular dynamics simulations in conjunction with home-developed Split Integration Symplectic Method we effectively decouple individual degrees of freedom of water molecules and connect them to corresponding thermostats. In this way, we facilitate elucidation of structural, dynamical, spectral, and hydration properties of bulk water at any given combination of rotational, translational, and vibrational temperatures. Elevated rotational temperature of the water medium is found to severely hinder hydration of polar molecules, to affect hydration of ionic species in a nonmonotonous way and to somewhat improve hydration of nonpolar species. As proteins consist of charged, polar, and nonpolar amino-acid residues, the developed methodology is also applied to critically evaluate the hypothesis that the overall decrease in protein hydration and the change in the subtle balance between hydration of various types of amino-acid residues provide a plausible physical mechanism through which microwaves enhance aberrant protein folding and aggregation.     

Thermodynamic properties and conductivities of some dodecylsurfactants in water

Densities, heat capacities, enthalpies of dilution, osmotic coefficients and conductivities are reported for dodecylamine hydrochloride, dodecyldimethylammonium and dodecyltrimethylammonium chloride in water over a wide range of concentration. The last two properties were also measured for dodecyltrimethylammonium bromide. From the thermodynamic data partial molar volumes, heat capacities and relative enthalpies and nonideal free energies and entropies were derived as a function of the surfactant concentration. The cmc's and degree of counterion dissociation were also calculated from the transport properties. It is shown that the trends of volumes, enthalpies, free energies and entropies are quite regular whereas heat capacities present maxima and minima at concentrations which depend on the nature of surfactants. Corresponding changes were observed in the osmotic coefficients and specific conductivities. The thermodynamic functions of micellization were evaluated on the basis of the pseudo-phase transition model. Finally, the effects of the introduction of methyl groups in the hydrophilic moiety of the surfactant and of the nature of the counterion on the thermodynamic properties of monomers and micelles are examined.     

Thermodynamic properties of N-octyl- and N-dodecylnicotinamide chlorides in water

Densities, heat capacities and enthalpies of dilution at 25°C and osmotic coefficients at 37°C were measured for N-octyl- and N-dodecylnicotinamide chlorides in water over an extended concentration region. Partial molar volumes, heat capacities, relative enthalpies and nonideal free energies and entropies at 25°C were derived as a function of the surfactant concentration. For both surfactants, plots of volumes, enthalpies and free energies vs. concentration are regular whereas those of heat capacities and entropies present anomalies at about 0.8 and 0.1m for the octyl and dodecyl compounds, respectively. Changes in the slope of a plot of osmotic coefficients times molality vs. molality were also observed at these same concentrations. These peculiarities are ascribed to micelle structural transitions. The nonideal free energies do not seem to depend on the alkyl chain length when they are plotted vs. m/C _cmc . Also, a plot of the nonideal free energy vs. logm/C _cmc is roughly independent of the nature of the surfactant because of the constant activity of surfactants in micellar solutions. Nonideal free energies, enthalpies and entropies have been calculated at 15 and 35°C. At each concentration the nonideal free energy is temperature independent as a result of a compensatory effect between enthalpy and entropy. The thermodynamic functions of micellization were graphically evaluated on the basis of the pseudo-phase transition model. These data suggest that the nicotinamide group possesses less hydrophilic character than the ammonium group.     

Theory of electron solvation in polar liquids: a continuum model

The solvation of electrons in polar liquids is analyzed on the basis of an extended continuum model. In addition to the long-range electron-dipole interaction two short-range interactions are introduced. Among others one accounts for interactions with groups capable of forming hydrogen bonds and the second for quadrupolar characteristics of the liquid molecules. Both are induced by the orientation of the molecular dipole. Applying the scaling method a proper reaction coordinate is introduced and the solvation dynamics are discussed for the electron in the electronic ground state and after excitation to the p-type excited state. The observed spectral evolution of the transient absorption spectra, after two photon excitations for electrons in water and in methanol, is well described by this theory. An analytic estimate for the nonradiative deactivation from the electronically excited solvated electron is found to be consistent with an observed lifetime of 50 fs for the electron in water. The theory predicts an about three times slower internal conversion in methanol as solvent in comparison with water.     

Ion solvation in water clusters

We describe here molecular dynamics computer simulations performed to study the solvation of ions (Cl^? and Na⁺) in water clusters. Our simulations show that the calculated structure and dynamics of the clusters is very sensitive to the potential model which is used to describe the interactions. From the comparison with thermodynamic data and data from the photoelectron spectra we conclude that in Cl^?(H₂O)_n (n≤20) clusters the ion is located on the surface of the cluster.     

Thermodynamic properties and structure of water films in flotation processes

The thermodynamic properties of water layers on barite and marble samples covered with various amounts of tetradecylammonium chloride (TDAHCl) after flotation processes were investigated by thermal analysis. The relation between the activation energy, the enthalpy and the entropy of the bonded water and the amount of TDAHCl previously deposited on the bare mineral surfaces was determined. An interpretation of the changes in these parameters, the water structure and the flotability is described.

---

Zusammenfassung Mittels Thermoanalyse wurden nach Schwimmaufbereitungsverfahren die thermodynamischen Eigenschaften der Wasserschichten über Baryt- und Marmorproben mit verschiedenem Gehalt an Tetradecylammoniumchlorid (TDACl) untersucht. Es wurde die Beziehung zwischen der Aktivierungsenergie, der Enthalpie und der Entropie des gebundenen Wassers und der Menge an zuvor an der Oberfläche des bloßen Minerals festgesetzen TDACl festgestellt. Die Veränderlichkeit dieser Parameter sowie die bindung des Wassers und di Flotierbarkeit werden interpretiert.

     

A molecular model of water based on the lattice gas model

A lattice gas model is used to describe the vapor-liquid state of water molecules. The orientationally directed interaction of the water molecules via their tetrahedral structure and dipole-dipole interaction are considered in the theory, along with the Lennard-Jones contributions to the potential of molecular interaction, which stabilize the system with dipole interaction. We studied how the radius of the molecular interaction potential affects the equilibrium characteristics of the system (the phase separation curves of the vapor-liquid system, and the relationship between the fluid density and the chemical potential value).     

VBSM: a solvation model based on valence bond theory

A new solvation model, called VBSM, is presented. The model combines valence bond (VB) theory with parameters determined for the SM6 solvation model (Kelly, C. P.; Cramer, C. J.; Truhlar, D. G. J. Chem. Theo. Comp. 2005, 1, 1133-1152). VBSM, like SM6, is based on the generalized Born (GB) approximation for bulk electrostatics and atomic surface tensions to account for cavitation, dispersion, and solvent structure (CDS). The solvation free energy of VBSM includes (i) a self-consistent polarization term obtained by using VB atomic charges in a GB reaction field with a VB self-consistent field procedure that minimizes the total energy of the system with respect to the valence bond orbitals and (ii) a geometry-dependent CDS term to account for deviations from bulk-electrostatic solvation. Test calculations for a few systems show that the liquid-phase partial atomic charges obtained by VBSM are in good agreement with liquid-phase charges obtained by charge model CM4 (Kelly, C. P.; Cramer, C. J.; Truhlar, D. G. J. Chem. Theo. Comp. 2005, 1, 1133-1152). Free energies of solvation are calculated for two prototype test cases, namely, for the degenerate S(N)2 reaction of Cl(-) with CH(3)Cl in water and for a Menshutkin reaction in water. These calculations show that the VBSM method provides a practical alternative to single-configuration self-consistent field theory for solvent effects in molecules and chemical reactions.     

Thermodynamic properties of some N-alkyl-N-methylpiperidinium chlorides and N-alkylpiperidine hydrochlorides in water

Densities and heat capacities at 25°C were measured for N-octyl-, N-decyl- and N-dodecyl-N-methylpiperidinium chlorides and for N-octyl- and N-dodecylpiperidine hydrochlorides in water as functions of concentration. Enthalpies of dilution at 25°C and osmotic coefficients at 37°C of the N-methyl-N-alkylpiperidinium chlorides were also measured as functions of concentration. The partial molar volumes, heat capacities, relative enthalpies, nonideal Gibbs energies and entropies at 25°C were derived as functions of the surfactant concentration. By increasing the alkyl chain length of the surfactant, both the apparent molar volume vs. concentration curves are shifted toward greater values while the corresponding ones for the heat capacity are moved toward more negative values. These results are consistent with the higher hydrophobicity the longer the alkyl chain of the surfactant is. In the micellar region, the entropy and enthalpy vs. log m/m _cmc curves increase in a parallel manner by decreasing the alkyl chain length of the surfactant. Consequently, the negligible effect of the hydrophobicity of the surfactant on the Gibbs energy vs. log m/m _cmc trends is due to the enthalpy-entropy compensative effect. The thermodynamic functions of micellization were graphically evaluated on the basis of the pseudo-phase transition model. The absolute values of both the volume and heat capacity of micellization increase with an increasing number of carbon atoms in the alkyl chain (n _c ). The enthalpy and entropy of micellization vs. n _c are convex curves. Comparisons are also made between the present data and those of some alkylpyridinium chlorides reported elsewhere.     

Thermodynamic properties of N-octyl-, N-decyl- and N-dodecylpyridinium chlorides in water

Densities, heat capacities and enthalpies of dilution at 25°C and osmotic coefficients at 37°C were measured for N-octyl-, N-decyl- and N-dodecyl-pyridinium chlorides in water over a wide concentration region. Conductivity measurements were performed in order to evaluate the cmc and the degree of counterion dissociation. Partial molar volumes, heat capacities, relative enthalpies and nonideal free energies and entropies at 25°C were derived from the experimental data as functions of the surfactant concentration. The changes with concentration of these properties are quite regular with the exception of the heat capacities which display anomalies at about 0.9, 0.25 and 0.12 mol-kg^–1 for the octyl, decyl and dodecyl compounds, respectively. At these concentrations there were also changes in the slopes of the specific conductivity and of the product of the osmotic coefficients and the molality vs. concentration. These peculiarities can be ascribed to micelle structural transitions. The thermodynamic functions of micellization were graphically evaluated on the basis of the pseudo-phase transition model. These data have been compared to those for alkyltrimethylammonium bromides and alkylnicotinamide chlorides. It is shown that the introduction of the hydrophilic CONH₂ group lowers the hydrophilic character of the pyridinium ring.     

Nucleation in hydrophobic cylindrical pores: a lattice model

We consider the nucleation process associated with capillary condensation of a vapor in a hydrophobic cylindrical pore (capillary evaporation). The liquid-vapor transition is described within the framework of a simple lattice model. The phase properties are characterized both at the mean-field level and with Monte Carlo simulations. The nucleation process for the liquid to vapor transition is then specifically considered. Using umbrella sampling techniques, we show that nucleation occurs through the condensation of an asymmetric vapor bubble at the pore surface. Even for highly confined systems, good agreement is found with macroscopic considerations based on classical nucleation theory. The results are discussed in the context of recent experimental work on the extrusion of water in hydrophobic pores.