1
|
Ibagon I, Furlan AP, Oliveira TJ, Dickman R. Phase diagram and critical properties of a two-dimensional associating lattice gas. Phys Rev E 2021; 104:064120. [PMID: 35030842 DOI: 10.1103/physreve.104.064120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Accepted: 12/01/2021] [Indexed: 06/14/2023]
Abstract
We revisit the associating lattice gas (ALG) introduced by Henriques et al. [Phys. Rev. E 71, 031504 (2005)PLEEE81539-375510.1103/PhysRevE.71.031504] in its symmetric version. In this model, defined on the triangular lattice, interaction between molecules occupying nearest-neighbor sites depends on their relative orientation, mimicking the formation of hydrogen bonds in network-forming fluids. Although all previous studies of this model agree that it has a disordered fluid (DF), a low-density liquid (LDL), and a high-density liquid (HDL) phase, quite different forms have been reported for its phase diagram. Here, we present a thorough investigation of its phase behavior using both transfer matrix calculations and Monte Carlo (MC) simulations, along with finite-size scaling extrapolations. Results in striking agreement are found using these methods. The critical point associated with the DF-HDL transition at full occupancy, identified by Furlan et al. [Phys. Rev. E 100, 022109 (2019)2470-004510.1103/PhysRevE.100.022109] is shown to be one terminus of a critical line separating these phases. In opposition to previous simulation studies, we find that the transition between the DF and LDL phases is always discontinuous, similar to the LDL-HDL transition. The associated coexistence lines meet at the point where the DF-HDL critical line ends, making it a critical-end-point. Overall, the form of the phase diagram observed in our simulations is very similar to that found in the exact solution of the model on a Husimi lattice. Our results confirm that, despite the existence of some waterlike anomalies in this model, it is unable to reproduce key features of the phase behavior of liquid water.
Collapse
Affiliation(s)
- I Ibagon
- Departamento de Física, ICEx, Universidade Federal de Minas Gerais, C.P. 702, 30123-970, Belo Horizonte, Minas Gerais-Brazil
| | - A P Furlan
- Departamento de Física, ICEx, Universidade Federal de Minas Gerais, C.P. 702, 30123-970, Belo Horizonte, Minas Gerais-Brazil
| | - T J Oliveira
- Departamento de Física, Universidade Federal de Viçosa, 36570-900, Viçosa, Minas Gerais, Brazil
| | - R Dickman
- Departamento de Física and National Institute of Science and Technology for Complex Systems, ICEx, Universidade Federal de Minas Gerais, C.P. 702, 30123-970 Belo Horizonte, Minas Gerais-Brazil
| |
Collapse
|
2
|
Moradi N, Greiner A, Melchionna S, Rao F, Succi S. A hydro-kinetic scheme for the dynamics of hydrogen bonds in water-like fluids. Phys Chem Chem Phys 2015; 16:15510-8. [PMID: 24953220 DOI: 10.1039/c4cp00921e] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A hydro-kinetic scheme for water-like fluids, based on a lattice version of the Boltzmann equation, is presented and applied to the popular TIP3P model of liquid water. By proceeding in much larger steps than molecular dynamics, the scheme proves to be very effective in attaining global minima of classical pair potentials with directional and radial interactions, as confirmed by further simulations using the three-dimensional Ben-Naim water potential. The scheme is shown to reproduce the propensity of water to form nearly four hydrogen bonds per molecule, as well as their statistical distribution in the presence of thermal fluctuations, at a linear cost of computational time with the system size.
Collapse
Affiliation(s)
- Nasrollah Moradi
- Freiburg Institute for Advanced Studies (FRIAS), University of Freiburg, Albertstrasse 19, 79104, Freiburg, Germany.
| | | | | | | | | |
Collapse
|
3
|
Roberts CJ, Blanco MA. Role of anisotropic interactions for proteins and patchy nanoparticles. J Phys Chem B 2014; 118:12599-611. [PMID: 25302767 PMCID: PMC4226310 DOI: 10.1021/jp507886r] [Citation(s) in RCA: 46] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
![]()
Protein–protein
interactions are inherently anisotropic
to some degree, with orientation-dependent interactions between repulsive
and attractive or complementary regions or “patches”
on adjacent proteins. In some cases it has been suggested that such
patch–patch interactions dominate the thermodynamics of dilute
protein solutions, as captured by the osmotic second virial coefficient
(B22), but delineating when this will
or will not be the case remains an open question. A series of simplified
but exactly solvable models are first used to illustrate that a delicate
balance exists between the strength of attractive patch–patch
interactions and the patch size, and that repulsive patch–patch
interactions contribute significantly to B22 for only those conditions where the repulsions are long-ranged.
Finally, B22 is reformulated, without
approximations, in terms of the density of states for a given interaction
energy and particle–particle distance. Doing so illustrates
the inherent balance of entropic and energetic contributions to B22. It highlights that simply having strong
patch–patch interactions will only cause anisotropic interactions
to dominate B22 solution properties if
the unavoidable entropic penalties are overcome, which cannot occur
if patches are too small. The results also indicate that the temperature
dependence of B22 may be a simple experimental
means to assess whether a small number of strongly attractive configurations
dominate the dilute solution behavior.
Collapse
Affiliation(s)
- Christopher J Roberts
- Department of Chemical and Biomolecular Engineering, and Center for Molecular and Engineering Thermodynamics, University of Delaware , Newark, Delaware 19716, United States
| | | |
Collapse
|
4
|
Gavazzoni C, Gonzatti GK, Pereira LF, Ramos LHC, Netz PA, Barbosa MC. The role of the anisotropy on the solid-fluid phase transition in core-softened shoulder-dumbbells systems. J Chem Phys 2014. [DOI: 10.1063/1.4871110] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
|
5
|
Pȩkalski J, Ciach A, Almarza NG. Periodic ordering of clusters in a one-dimensional lattice model. J Chem Phys 2013; 138:144903. [DOI: 10.1063/1.4799264] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
6
|
Cartwright JHE, Piro O, Sánchez PA, Sintes T. Ice polyamorphism in the minimal Mercedes-Benz model of water. J Chem Phys 2012; 137:244503. [PMID: 23277941 DOI: 10.1063/1.4772202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Julyan H E Cartwright
- Instituto Andaluz de Ciencias de la Tierra, CSIC-Universidad de Granada, E-18071 Granada, Spain.
| | | | | | | |
Collapse
|
7
|
Szortyka MM, Fiore CE, Barbosa MC, Henriques VB. Hydration and anomalous solubility of the Bell-Lavis model as solvent. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2012; 86:031503. [PMID: 23030919 DOI: 10.1103/physreve.86.031503] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/08/2012] [Indexed: 06/01/2023]
Abstract
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.
Collapse
Affiliation(s)
- Marcia M Szortyka
- Departamento de Física, Universidade Federal de Santa Catarina, Caixa Postal 476, 88010-970 Florianópolis, Santa Catarina, Brazil.
| | | | | | | |
Collapse
|
8
|
Oliveira TJ, Stilck JF, Barbosa MAA. Solution of an associating lattice-gas model with density anomaly on a Husimi lattice. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2010; 82:051131. [PMID: 21230461 DOI: 10.1103/physreve.82.051131] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/09/2010] [Indexed: 05/30/2023]
Abstract
We study a model of a lattice gas with orientational degrees of freedom which resemble the formation of hydrogen bonds between the molecules. In this model, which is the simplified version of the Henriques-Barbosa model, no distinction is made between donors and acceptors in the bonding arms. We solve the model in the grand-canonical ensemble on a Husimi lattice built with hexagonal plaquettes with a central site. The ground state of the model, which was originally defined on the triangular lattice, is exactly reproduced by the solution on this Husimi lattice. In the phase diagram, one gas and two liquid [high density liquid (HDL) and low density liquid (LDL)] phases are present. All phase transitions (GAS-LDL, GAS-HDL, and LDL-HDL) are discontinuous, and the three phases coexist at a triple point. A line of temperatures of maximum density in the isobars is found in the metastable GAS phase, as well as another line of temperatures of minimum density appears in the LDL phase, part of it in the stable region and another in the metastable region of this phase. These findings are at variance with simulational results for the same model on the triangular lattice, which suggested a phase diagram with two critical points. However, our results show very good quantitative agreement with the simulations, both for the coexistence loci and the densities of particles and of hydrogen bonds. We discuss the comparison of the simulations with our results.
Collapse
Affiliation(s)
- Tiago J Oliveira
- Instituto de Física, Universidade Federal Fluminense, Av. Litorânea s/n, 24210-340, Niterói, RJ, Brazil
| | | | | |
Collapse
|
9
|
Szortyka MM, Fiore CE, Henriques VB, Barbosa MC. Diffusion anomaly and dynamic transitions in the Bell–Lavis water model. J Chem Phys 2010; 133:104904. [DOI: 10.1063/1.3479001] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
10
|
Barros de Oliveira A, Neves EB, Gavazzoni C, Paukowski JZ, Netz PA, Barbosa MC. Liquid crystal phase and waterlike anomalies in a core-softened shoulder-dumbbells system. J Chem Phys 2010; 132:164505. [DOI: 10.1063/1.3386384] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
11
|
Szortyka MM, Girardi M, Henriques VB, Barbosa MC. Dynamic transitions in a three dimensional associating lattice gas model. J Chem Phys 2010; 132:134904. [DOI: 10.1063/1.3354112] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
12
|
Pretti M, Buzano C, De Stefanis E. Revisiting waterlike network-forming lattice models. J Chem Phys 2010; 131:224508. [PMID: 20001058 DOI: 10.1063/1.3270000] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We revisit different three-dimensional network-forming lattice models proposed in the literature to investigate water anomalies. We perform a semianalytical calculation based on a cluster-variation technique, showing a quite good agreement with independent Monte Carlo results. The method allows us to clarify the structure of the phase diagrams, which turn out to exhibit different kinds of orientationally ordered phases. We point out that certain "waterlike" thermodynamic anomalies, claimed by previous studies, are indeed artifacts of a homogeneity assumption made in the analytical treatment. We argue that such a difficulty is common to a whole class of lattice models for water and suggest a possible way to overcome the problem in terms of "equivalent" models defined on random lattices.
Collapse
Affiliation(s)
- M Pretti
- Center for Statistical Mechanics and Complexity, CNR-INFM Roma 1, Piazzale Aldo Moro 2, I-00185 Roma, Italy.
| | | | | |
Collapse
|
13
|
Buzano C, De Stefanis E, Pretti M. Cluster-variation approximation for a network-forming lattice-fluid model. J Chem Phys 2008; 129:024506. [DOI: 10.1063/1.2919126] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
14
|
Barbosa MAA, Henriques VB. Frustration and anomalous behavior in the Bell-Lavis model of liquid water. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:051204. [PMID: 18643057 DOI: 10.1103/physreve.77.051204] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/18/2008] [Indexed: 05/26/2023]
Abstract
We have reconsidered the Bell-Lavis model of liquid water and investigated its relation to its isotropic version, the antiferromagnetic Blume-Emery-Griffiths model on the triangular lattice. Our study was carried out by means of an exact solution on the sequential Husimi cactus. We show that the ground states of both models share the same topology and that fluid phases (gas and low- and high-density liquids) can be mapped onto magnetic phases (paramagnetic, antiferromagnetic, and dense paramagnetic, respectively). Both models present liquid-liquid coexistence and several thermodynamic anomalies. This result suggests that anisotropy introduced through orientational variables play no specific role in producing the density anomaly, in agreement with a similar conclusion discussed previously following results for continuous soft core models. We propose that the presence of liquid anomalies may be related to energetic frustration, a feature common to both models.
Collapse
|
15
|
Girardi M, Balladares AL, Henriques VB, Barbosa MC. Liquid polymorphism and density anomaly in a three-dimensional associating lattice gas. J Chem Phys 2007; 126:064503. [PMID: 17313225 DOI: 10.1063/1.2434974] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Mauricio Girardi
- Universidade Federal de Pelotas-UNIPAMPA/Bagé, Rua Carlos Barbosa SN, CEP 96400-970 Bagé, Rio Grande do Sul, Brazil.
| | | | | | | |
Collapse
|
16
|
Guisoni N, Henriques VB. Hydrophobic Hydration in an Orientational Lattice Model. J Phys Chem B 2006; 110:17188-94. [PMID: 16928016 DOI: 10.1021/jp060729f] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
To shed light on the microscopic mechanism of hydrophobic hydration, we study a simplified lattice model for water solutions in which the orientational nature of hydrogen bonding as well as the degeneracy related to proton distribution are taken into account. Miscibility properties of the model are looked at for both polar (hydrogen bonding) and nonpolar (non-hydrogen bonding) solutes. A quasichemical solution for the pure system is reviewed and extended to include the different kinds of solute. A Monte Carlo study of our model yields a novel feature for the local structure of the hydration layer: energy correlation relaxation times for solvation water are larger than the corresponding relaxation times for bulk water. This result suggests the presence of ordering of water particles in the first hydration shell. A nonassociating model solvent, represented by a lattice gas, presents opposite behavior, indicating that this effect is a result of the directionality of the interaction. In presence of polar solutes, we find an ordered mixed pseudophase at low temperatures, indicating the possibility of closed loops of immiscibility.
Collapse
Affiliation(s)
- Nara Guisoni
- Instituto de Pesquisa e Desenvolvimento, Universidade do Vale do Paraíba (UNIVAP), Av. Shishima Hifumi, 2911, Urbanova cep 12244-000, São José dos Campos, SP, Brazil.
| | | |
Collapse
|
17
|
Pretti M, Buzano C. Thermodynamic anomalies in a lattice model of water: Solvation properties. J Chem Phys 2005; 123:24506. [PMID: 16050758 DOI: 10.1063/1.1950628] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
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. The model is similar to the one proposed by Roberts and Debenedetti [J. Chem. Phys. 105, 658 (1996)], simplified by removing distinction between "donors" and "acceptors." We focus on the solvation properties, mainly as far as an ideally inert (hydrophobic) solute is concerned. As in our previous analysis, devoted to neat water [J. Chem. Phys. 121, 11856 (2004)], we make use of a generalized first-order approximation on a tetrahedral cluster. We show that the model exhibits quite a coherent picture of water thermodynamics, reproducing qualitatively several anomalous properties observed both in pure water and in solutions of hydrophobic solutes. As far as supercooled liquid water is concerned, the model is consistent with the second critical-point scenario.
Collapse
Affiliation(s)
- M Pretti
- Istituto Nazionale per la Fisica della Materia (INFM) and Dipartimento di Fisica, Politecnico di Torino, Corso Duca degli Abruzzi 24, I-10129 Torino, Italy.
| | | |
Collapse
|
18
|
Lishchuk SV, Lokotosh TV, Malomuzh NP. Properties of the H-bond network for two-dimensional lattice water model. J Chem Phys 2005; 122:244504. [PMID: 16035779 DOI: 10.1063/1.1940030] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A microscopic Hamiltonian of the hydrogen-bond network in two-dimensional lattice water is proposed, which describes the formation and disruption of the H bonds, their bending, and which satisfies the Bernal-Fowler rules [J. D. Bernal and R. H. Fowler, J. Chem. Phys. 1, 515 (1933)]. The thermodynamic properties of the H-bond network are studied using the method of many-particle irreducible distribution functions, which is a generalization of the Kikuchi cluster approach [R. Kikuchi, Phys. Rev. 81, 988 (1951)] and the Bethe-Peierls quasiactivities method [H. A. Bethe, Prog. R. Soc. A 150, 552 (1935)]. The temperature dependencies of the average number of H bonds per molecules, the contribution of the H bonds into the heat capacity of the system, and the parameters describing the correlations between the states of molecules on the neighboring sites are investigated. It is shown that depending on the magnitude of the interaction between the H bonds in the H-bond subsystem either smooth or sharp first-order phase transition can occur. The role of different factors in the formation of the properties of the H-bond network is discussed.
Collapse
|