The associating lattice gas in the presence of interacting solutes.

We have investigated the phase diagram of a statistical model for hydrogen-bonding solutions for polar solutes. The structured solvent is represented by an associating lattice gas, which presents anomalous density and liquid-liquid coexistence. Polar solute particles and solvent particles interact attractively, while the solvent-solvent interaction is made directional through bonding arms, which mimic hydrogen bonds. The model behavior is obtained via Monte Carlo simulations in the grand-canonical ensemble, for different sets of parameters. For small solute chemical potential and weak attraction between solute and solvent particles, addition of solute yields a shift in the transition lines of the pure solvent. This is the scenario explored by different authors, in the pursuit of stabilizing the water liquid-liquid coexistence line. However, as we show, in the case of larger solute chemical potentials, or of stronger solute-solvent attractions, new phases may arise.

Structure and anomalous solubility for hard spheres in an associating lattice gas model.

In this paper we investigate the solubility of a hard-sphere gas in a solvent modeled as an associating lattice gas. The solution phase diagram for solute at 5% is compared with the phase diagram of the original solute free model. Model properties are investigated both through Monte Carlo simulations and a cluster approximation. The model solubility is computed via simulations and is shown to exhibit a minimum as a function of temperature. The line of minimum solubility (TmS) coincides with the line of maximum density (TMD) for different solvent chemical potentials, in accordance with the literature on continuous realistic models and on the "cavity" picture.

Lipid bilayer pre-transition as the beginning of the melting process.

We investigate the bilayer pre-transition exhibited by some lipids at temperatures below their main phase transition, and which is generally associated to the formation of periodic ripples in the membrane. Experimentally we focus on the anionic lipid dipalmytoylphosphatidylglycerol (DPPG) at different ionic strengths, and on the neutral lipid dipalmytoylphosphatidylcholine (DPPC). From the analysis of differential scanning calorimetry traces of the two lipids we find that both pre- and main transitions are part of the same melting process. Electron spin resonance of spin labels and excitation generalized polarization of Laurdan reveal the coexistence of gel and fluid domains at temperatures between the pre- and main transitions of both lipids, reinforcing the first finding. Also, the melting process of DPPG at low ionic strength is found to be less cooperative than that of DPPC. From the theoretical side, we introduce a statistical model in which a next-nearest-neighbor competing interaction is added to the usual two-state model. For the first time, modulated phases (ordered and disordered lipids periodically aligned) emerge between the gel and fluid phases as a natural consequence of the competition between lipid-lipid interactions.

Ionization and structural changes of the DMPG vesicle along its anomalous gel-fluid phase transition: a study with different lipid concentrations.

Dispersions of saturated anionic phospholipid dimyristoyl phosphatidylglycerol (DMPG) have been extensively studied regarding their peculiar thermostructural behavior. At low ionic strength, the gel-fluid transition is spread along nearly 17 degrees C, displaying several thermal events in the calorimetric profile that is quite different from the single sharp peak around 23 degrees C found for higher ionic strength DMPG dispersions. To investigate the role of charge in the bilayer transition, we carefully examine the temperature dependence of the electrical conductivity of DMPG dispersions at different concentrations, correlating the data with the corresponding differential scanning calorimetry (DSC) traces. Electrical conductivity together with electrophoretic mobility measurements allowed the calculation of the dependence of the degree of ionization of DMPG vesicles on lipid concentration and temperature. It was shown that there is a decrease in vesicle charge as the lipid concentration increases, which is probably correlated with the increase in the concentration of bulk Na(+). Apart from the known increase in the electrical conductivity along the DMPG temperature transition region, a sharp rise was observed at the bilayer pretransition for all lipid concentrations studied, possibly indicating that the beginning of the chain melting process is associated with an increase in bilayer ionization. It is confirmed here that the gel-fluid transition of DMPG at low ionic strength is accompanied by a huge increase in the dispersion viscosity. However, it is shown that this measured macroviscosity is distinct from the local viscosity felt by either charged ions or DMPG charged aggregates in measurements of electrical conductivity or electrophoretic mobility. Data presented here give support to the idea that DMPG vesicles, at low ionic strength, get more ionized along the temperature transition region and could be perforated and/or deformed vesicle structures.

Dynamic transitions in a three dimensional associating lattice gas model.

We investigate the thermodynamic and dynamic properties of a three dimensional associating lattice gas (ALG) model through Monte Carlo simulations. The ALG model combines a soft core potential and orientational degrees of freedom. The competition of directional attractive forces and the soft core potential results in two coexisting liquid phases which are also connected through order-disorder critical transitions. The model presents structural order, density, and diffusion anomalies. Our study suggests that the dynamic fragile-to-strong transitions are associated to changes in structural order.

Diffusion anomaly and dynamic transitions in the Bell-Lavis water model.

In this paper we investigate the dynamic properties of the minimal Bell-Lavis (BL) water model and their relation to the thermodynamic anomalies. The BL model is defined on a triangular lattice in which water molecules are represented by particles with three symmetric bonding arms interacting through van der Waals and hydrogen bonds. We have studied the model diffusivity in different regions of the phase diagram through Monte Carlo simulations. Our results show that the model displays a region of anomalous diffusion which lies inside the region of anomalous density, englobed by the line of temperatures of maximum density. Further, we have found that the diffusivity undergoes a dynamic transition which may be classified as fragile-to-strong transition at the critical line only at low pressures. At higher densities, no dynamic transition is seen on crossing the critical line. Thus evidence from this study is that relation of dynamic transitions to criticality may be discarded.

Dynamic transitions in a two dimensional associating lattice gas model.

Using Monte Carlo simulations we investigate some new aspects of the phase diagram and the behavior of the diffusion coefficient in an associating lattice gas (ALG) model on different regions of the phase diagram. The ALG model combines a two dimensional lattice gas where particles interact through a soft core potential and orientational degrees of freedom. The competition between soft core potential and directional attractive forces results in a high density liquid phase, a low density liquid phase, and a gas phase. Besides anomalies in the behavior of the density with the temperature at constant pressure and of the diffusion coefficient with density at constant temperature are also found. The two liquid phases are separated by a coexistence line that ends in a bicritical point. The low density liquid phase is separated from the gas phase by a coexistence line that ends in tricritical point. The bicritical and tricritical points are linked by a critical lambda-line. The high density liquid phase and the fluid phases are separated by a second critical tau-line. We then investigate how the diffusion coefficient behaves on different regions of the chemical potential-temperature phase diagram. We find that diffusivity undergoes two types of dynamic transitions: a fragile-to-strong transition when the critical lambda-line is crossed by decreasing the temperature at a constant chemical potential; and a strong-to-strong transition when the critical tau-line is crossed by decreasing the temperature at a constant chemical potential.

Liquid polymorphism, order-disorder transitions and anomalous behavior: A Monte Carlo study of the Bell-Lavis model for water.

The Bell-Lavis model for liquid water is investigated through numerical simulations. The lattice-gas model on a triangular lattice presents orientational states and is known to present a highly bonded low density phase and a loosely bonded high density phase. We show that the model liquid-liquid transition is continuous, in contradiction with mean-field results on the Husimi cactus and from the cluster variational method. We define an order parameter which allows interpretation of the transition as an order-disorder transition of the bond network. Our results indicate that the order-disorder transition is in the Ising universality class. Previous proposal of an Ehrenfest second order transition is discarded. A detailed investigation of anomalous properties has also been undertaken. The line of density maxima in the HDL phase is stabilized by fluctuations, absent in the mean-field solution.

Liquid polymorphism and density anomaly in a lattice gas model.

We present a simple model for an associating liquid in which polymorphism and density anomaly are connected. Our model combines a two dimensional lattice gas with particles interacting through a soft core potential and orientational degrees of freedom represented through thermal "ice variables." The competition between the directional attractive forces and the soft core potential leads to a phase diagram in which two liquid phases and a density anomaly are present. The coexistence line between the low density liquid and the high density liquid has a positive slope contradicting the surmise that the presence of a density anomaly implies that the high density liquid is more entropic than the low density liquid.

Microcanonical Monte Carlo simulation of lattice gas models.

We introduce a true microcanonical ensemble appropriate for Monte Carlo simulations of lattice gas models in which both the energy and the number of particles are held fixed. We also deduce formulas that allow us to do the calculation of the intensive quantities: temperature and chemical potential. The numerical results coming from Monte Carlo simulations defined according to this microcanonical ensemble compare well with results coming from the ordinary canonical ensemble.

Lattice solution model for order-disorder transitions in membranes and Langmuir monolayers.

Lipid monolayers and bilayers have been used as experimental models for the investigation of membrane thermal transitions. The main transition takes place near ambient temperatures for several lipids and reflects the order-disorder transition of lipid hydrocarbonic chains, which is accompanied by a surface density gap. Equivalence between the transitions in the two systems has been argued by several authors. The two-state statistical model adopted by numerous authors for different properties of the membrane, such as permeability, diffusion, and mixture or insertion of cholesterol or protein, is inadequate for the description of charged membranes, since it lacks a proper description of surface density. We propose a lattice solution model which adds interactions with water molecules to lipid-lipid interactions and obtain its thermal properties under a mean-field approach. Density variations, although concomitant with chain order variations, are independent of the latter. The model presents both chain order and gas-liquid transitions, and extends the range of applicability of previous models, yielding Langmuir isotherms in the full range of pressures and areas.

Light scattering on the structural characterization of DMPG vesicles along the bilayer anomalous phase transition.

Highly charged vesicles of the saturated anionic lipid dimyristoyl phosphatidylglycerol (DMPG) in low ionic strength medium exhibit a very peculiar thermo-structural behavior. Along a wide gel-fluid transition region, DMPG dispersions display several anomalous characteristics, like low turbidity, high electrical conductivity and viscosity. Here, static and dynamic light scattering (SLS and DLS) were used to characterize DMPG vesicles at different temperatures. Similar experiments were performed with the largely studied zwitterionic lipid dimyristoyl phosphatidylcholine (DMPC). SLS and DLS data yielded similar dimensions for DMPC vesicles at all studied temperatures. However, for DMPG, along the gel-fluid transition region, SLS indicated a threefold increase in the vesicle radius of gyration, whereas the hydrodynamic radius, as obtained from DLS, increased 30% only. Despite the anomalous increase in the radius of gyration, DMPG lipid vesicles maintain isotropy, since no light depolarization was detected. Hence, SLS data are interpreted regarding the presence of isotropic vesicles within the DMPG anomalous transition, but highly perforated vesicles, with large holes. DLS/SLS discrepancy along the DMPG transition region is discussed in terms of the interpretation of the Einstein-Stokes relation for porous vesicles. Therefore, SLS data are shown to be much more appropriate for measuring porous vesicle dimensions than the vesicle diffusion coefficient. The underlying nanoscopic process which leads to the opening of pores in charged DMPG bilayer is very intriguing and deserves further investigation. One could envisage biotechnological applications, with vesicles being produced to enlarge and perforate in a chosen temperature and/or pH value.

Hydration and anomalous solubility of the Bell-Lavis model as solvent.

We address the investigation of the solvation properties of the minimal orientational model for water originally proposed by [Bell and Lavis, J. Phys. A 3, 568 (1970)]. The model presents two liquid phases separated by a critical line. The difference between the two phases is the presence of structure in the liquid of lower density, described through the orientational order of particles. We have considered the effect of a small concentration of inert solute on the solvent thermodynamic phases. Solute stabilizes the structure of solvent by the organization of solvent particles around solute particles at low temperatures. Thus, even at very high densities, the solution presents clusters of structured water particles surrounding solute inert particles, in a region in which pure solvent would be free of structure. Solute intercalates with solvent, a feature which has been suggested by experimental and atomistic simulation data. Examination of solute solubility has yielded a minimum in that property, which may be associated with the minimum found for noble gases. We have obtained a line of minimum solubility (TmS) across the phase diagram, accompanying the line of maximum density. This coincidence is easily explained for noninteracting solute and it is in agreement with earlier results in the literature. We give a simple argument which suggests that interacting solute would dislocate TmS to higher temperatures.

Liquid polymorphism and density anomaly in a three-dimensional associating lattice gas.

The authors investigate the phase diagram of a three-dimensional associating lattice gas (ALG) model. This model combines orientational icelike interactions and "van der Waals" that might be repulsive, representing, in this case, a penalty for distortion of hydrogen bonds. These interactions can be interpreted as two competing distances, making the connection between this model and continuous isotropic soft-core potentials. The authors present Monte Carlo studies of the ALG model showing the presence of two liquid phases, two critical points, and density anomaly.

Canonical and microcanonical Monte Carlo simulations of lattice-gas mixtures.

We propose strict canonical and microcanonical Monte Carlo algorithms for an arbitrary lattice-gas binary mixture. We deduce formulas that allow us to obtain field quantities over the ensembles in which their conjugate extensive quantities are conserved. As an example, we have considered a lattice-gas mixture that is equivalent to the spin-1 Blume-Emery-Griffiths model [Phys. Rev. A 4, 1071 (1971)]. For a finite system and near a phase coexistence, the field as a function of its extensive conjugate shows a loop that disappears in the thermodynamic limit giving rise to the usual tie line. The first-order phase transition was determined by the use of three criteria.