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
|
Furlan AP, Oliveira TJ, Stilck JF, Dickman R. Order-disorder transition in a two-dimensional associating lattice gas. Phys Rev E 2019; 100:022109. [PMID: 31574678 DOI: 10.1103/physreve.100.022109] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/26/2019] [Indexed: 11/07/2022]
Abstract
We study an associating lattice gas (ALG) using Monte Carlo simulation on the triangular lattice and semianalytical solutions on Husimi lattices. In this model, the molecules have an orientational degree of freedom and the interactions depend on the relative orientations of nearest-neighbor molecules, mimicking the formation of hydrogen bonds. We focus on the transition between the high-density liquid (HDL) phase and the isotropic phase in the limit of full occupancy, corresponding to chemical potential μ→∞, which has not yet been studied systematically. Simulations yield a continuous phase transition at τ_{c}=k_{B}T_{c}/γ=0.4763(1) (where -γ is the bond energy) between the low-temperature HDL phase, with a nonvanishing mean orientation of the molecules, and the high-temperature isotropic phase. Results for critical exponents and the Binder cumulant indicate that the transition belongs to the three-state Potts model universality class, even though the ALG Hamiltonian does not have the full permutation symmetry of the Potts model. In contrast with simulation, the Husimi lattice analyses furnish a discontinuous phase transition, characterized by a discontinuity of the nematic order parameter. The transition temperatures (τ_{c}=0.51403 and 0.51207 for trees built with triangles and hexagons, respectively) are slightly higher than that found via simulation. Since the Husimi lattice studies show that the ALG phase diagram features a discontinuous isotropic-HDL line for finite μ, three possible scenarios arise for the triangular lattice. The first is that in the limit μ→∞ the first-order line ends in a critical point; the second is a change in the nature of the transition at some finite chemical potential; the third is that the entire line is one of continuous phase transitions. Results from other ALG models and the fact that mean-field approximations show a discontinuous phase transition for the three-state Potts model (known to possess a continuous transition) lends some weight to the third alternative.
Collapse
Affiliation(s)
- A P Furlan
- Departamento de Física, ICEx, Universidade Federal de Minas Gerais, C. P. 702, 30123-970 Belo Horizonte, Minas Gerais, Brazil
| | - Tiago J Oliveira
- Departamento de Física, Universidade Federal de Viçosa, 36570-900, Viçosa, Minas Gerais, Brazil
| | - Jürgen F Stilck
- Instituto de Física and National Institute of Science and Technology for Complex Systems, Universidade Federal Fluminense, Niterói, Rio de Janeiro, Brazil
| | - Ronald 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
|
3
|
Noya EG, Almarza NG, Lomba E. Assembly of trivalent particles under confinement: from an exotic solid phase to a liquid phase at low temperature. SOFT MATTER 2017; 13:3221-3229. [PMID: 28398440 DOI: 10.1039/c7sm00217c] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/25/2023]
Abstract
Using computer simulations, we study the phase diagram of a two-dimensional system of disk particles with three patches distributed symmetrically along the particle equator. The geometry of the particles is compatible with a honey-comb lattice at moderately low temperature and pressure, whereas it is expected that the system forms a close-packed triangular lattice at high temperature and pressure. The effect of patch size within the single bond per patch regime was investigated, and it was found that the topology of the phase diagram changes drastically with patch size. Interestingly, in particles with small patches (with a half opening angle of 10°), the fluid transforms upon increasing the pressure into a rather exotic phase that can be understood as a honey-comb lattice whose voids are filled continuously with additional particles that remain, on average, unbound. Eventually, all the voids are occupied so that particles are located at the positions of a triangular lattice, but only two thirds of the particles are orientationally ordered whereas the remaining one third can rotate almost freely as in a plastic crystal. At moderately low temperature, the fluid transforms into a nearly empty honey-comb lattice, whereas at high temperature it transforms directly into the almost filled lattice. Interestingly, for particles with big patches (with a half opening angle of 20°), the honey-comb and triangular lattices are separated by a liquid phase that remains stable down to fairly low temperatures. Less surprisingly, only particles with big patches exhibit an equilibrium gas-liquid separation.
Collapse
Affiliation(s)
- Eva G Noya
- Instituto de Química Física Rocasolano, Consejo Superior de Investigaciones Cientficas, CSIC, Calle Serrano 119, 28026 Madrid, Spain.
| | | | | |
Collapse
|
4
|
Gorbunov VA, Akimenko SS, Myshlyavtsev AV. Cross-impact of surface and interaction anisotropy in the self-assembly of organic adsorption monolayers: a Monte Carlo and transfer-matrix study. Phys Chem Chem Phys 2017. [DOI: 10.1039/c7cp01863k] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We have theoretically studied the features of self-assembly in organic adsorption layers where both “molecule–surface” and “molecule–molecule” interactions are anisotropic.
Collapse
Affiliation(s)
| | | | - A. V. Myshlyavtsev
- Omsk State Technical University
- Omsk
- Russian Federation
- Institute of Hydrocarbon Processing SB RAS
- Omsk
| |
Collapse
|
5
|
Šimėnas M, Ibenskas A, Tornau EE. Phase transition properties of the Bell-Lavis model. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 90:042124. [PMID: 25375455 DOI: 10.1103/physreve.90.042124] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/20/2014] [Indexed: 06/04/2023]
Abstract
Using Monte Carlo calculations we analyze the order and the universality class of phase transitions into a low-density (honeycomb) phase of a triangular antiferromagnetic three-state Bell-Lavis model. The results are obtained in a whole interval of chemical potential μ corresponding to the honeycomb phase. Our results demonstrate that the phase transitions might be attributed to the three-state Potts universality class for all μ values except for the edges of the honeycomb phase existence. At the honeycomb phase and the low-density gas phase boundary the transitions become of the first order. At another, honeycomb-to-frustrated phase boundary, we observe the approach to the crossover from the three-state Potts to the Ising model universality class. We also obtain the Schottky anomaly in the specific heat close to this edge. We show that the intermediate planar phase, found in a very similar antiferromagnetic triangular Blume-Capel model, does not occur in the Bell-Lavis model.
Collapse
Affiliation(s)
- M Šimėnas
- Semiconductor Physics Institute, Center for Physical Sciences and Technology, Goštauto 11, LT-01108 Vilnius, Lithuania
| | - A Ibenskas
- Semiconductor Physics Institute, Center for Physical Sciences and Technology, Goštauto 11, LT-01108 Vilnius, Lithuania
| | - E E Tornau
- Semiconductor Physics Institute, Center for Physical Sciences and Technology, Goštauto 11, LT-01108 Vilnius, Lithuania
| |
Collapse
|
6
|
Ibenskas A, Šimėnas M, Tornau EE. Antiferromagnetic triangular Blume-Capel model with hard-core exclusions. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:052144. [PMID: 25353776 DOI: 10.1103/physreve.89.052144] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/06/2014] [Indexed: 06/04/2023]
Abstract
Using Monte Carlo simulation, we analyze phase transitions of two antiferromagnetic (AFM) triangular Blume-Capel (BC) models with AFM interactions between third-nearest neighbors. One model has hard-core exclusions between the nearest-neighbor (1NN) particles (3NN1 model) and the other has them between the nearest-neighbor and next-nearest-neighbor particles (3NN12 model). Finite-size scaling analysis reveals that in these models, the transition from the paramagnetic to long-range order (LRO) AFM phase is either of the first order or goes through an intermediate phase which might be attributed to the Berezinskii-Kosterlitz-Thouless (BKT) type. The properties of the low-temperature phase transition to the AFM phase of the 1NN, 3NN1, and 3NN12 models are found to be very similar for almost all values of a normalized single-ion anisotropy parameter, 0 < δ < 1.5. Higher temperature behavior of the 3NN12 and 3NN1 models is rather different from that of the 1NN model. Three phase transitions are observed for the 3NN12 model: from the paramagnetic phase to the phase with domains of the LRO AFM phase at T(c), from this structure to the diluted frustrated BKT-type phase at T(2), and from the frustrated phase to the AFM LRO phase at T(1). For the 3NN12 model, T(c) > T(2) > T(1) at 0 < δ < 1.15 (range I), T(c) ≈ T(2) > T(1) at 1.15 < δ < 1.3 (range II), and T(c) = T(2) = T(1) at 1.3 < δ < 1.5 (range III). For the 3NN1 model, T(c) ≈ T(2) > T(1) at 0 < δ < 1.2 (range II) and T(c) = T(2) = T(1) at 1.2 < δ < 1.5 (range III). There is only one first-order phase transition in range III. The transition at T(c) is of the first order in range II and either of a weak first order or a second order in range I.
Collapse
Affiliation(s)
- A Ibenskas
- Semiconductor Physics Institute, Center for Physical Sciences and Technology, Goštauto 11, LT-01108 Vilnius, Lithuania
| | - M Šimėnas
- Semiconductor Physics Institute, Center for Physical Sciences and Technology, Goštauto 11, LT-01108 Vilnius, Lithuania
| | - E E Tornau
- Semiconductor Physics Institute, Center for Physical Sciences and Technology, Goštauto 11, LT-01108 Vilnius, Lithuania
| |
Collapse
|
7
|
Badasyan A, Tonoyan SA, Giacometti A, Podgornik R, Parsegian VA, Mamasakhlisov YS, Morozov VF. Unified description of solvent effects in the helix-coil transition. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:022723. [PMID: 25353524 DOI: 10.1103/physreve.89.022723] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/02/2014] [Indexed: 06/04/2023]
Abstract
We analyze the problem of the helix-coil transition in explicit solvents analytically by using spin-based models incorporating two different mechanisms of solvent action: explicit solvent action through the formation of solvent-polymer hydrogen bonds that can compete with the intrinsic intra-polymer hydrogen bonded configurations (competing interactions) and implicit solvent action, where the solvent-polymer interactions tune biopolymer configurations by changing the activity of the solvent (non-competing interactions). The overall spin Hamiltonian is comprised of three terms: the background in vacuo Hamiltonian of the "Generalized Model of Polypeptide Chain" type and two additive terms that account for the two above mechanisms of solvent action. We show that on this level the solvent degrees of freedom can be explicitly and exactly traced over, the ensuing effective partition function combining all the solvent effects in a unified framework. In this way we are able to address helix-coil transitions for polypeptides, proteins, and DNA, with different buffers and different external constraints. Our spin-based effective Hamiltonian is applicable for treatment of such diverse phenomena as cold denaturation, effects of osmotic pressure on the cold and warm denaturation, complicated temperature dependence of the hydrophobic effect as well as providing a conceptual base for understanding the behavior of intrinsically disordered proteins and their analogues.
Collapse
Affiliation(s)
- Artem Badasyan
- Materials Research Laboratory, University of Nova Gorica, Vipavska 13, SI-5000 Nova Gorica, Slovenia
| | - Shushanik A Tonoyan
- Department of Molecular Physics, Yerevan State University, A. Manougian Str. 1, 375025, Yerevan, Armenia
| | - Achille Giacometti
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca' Foscari Venezia, Calle Larga S. Marta DD2137, I-30123 Venezia, Italy
| | - Rudolf Podgornik
- Department of Theoretical Physics, J. Stefan Institute and Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana - SI-1000 Ljubljana, Slovenia and Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003-9337, USA
| | - V Adrian Parsegian
- Department of Physics, University of Massachusetts, Amherst, Massachusetts 01003-9337, USA
| | - Yevgeni Sh Mamasakhlisov
- Department of Molecular Physics, Yerevan State University, A. Manougian Str. 1, 375025, Yerevan, Armenia
| | - Vladimir F Morozov
- Department of Molecular Physics, Yerevan State University, A. Manougian Str. 1, 375025, Yerevan, Armenia
| |
Collapse
|
8
|
Fiore CE, da Luz MGE. Exploiting a semi-analytic approach to study first order phase transitions. J Chem Phys 2013; 138:014105. [DOI: 10.1063/1.4772809] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
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
|
10
|
Fiore CE, da Luz MGE. General approach for studying first-order phase transitions at low temperatures. PHYSICAL REVIEW LETTERS 2011; 107:230601. [PMID: 22182076 DOI: 10.1103/physrevlett.107.230601] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/11/2011] [Indexed: 05/31/2023]
Abstract
By combining different ideas, a general and efficient protocol to deal with discontinuous phase transitions at low temperatures is proposed. For small T's, it is possible to derive a generic analytic expression for appropriate order parameters, whose coefficients are obtained from simple simulations. Once in such regimes simulations by standard algorithms are not reliable; an enhanced tempering method, the parallel tempering-accurate for small and intermediate system sizes with rather low computational cost-is used. Finally, from finite size analysis, one can obtain the thermodynamic limit. The procedure is illustrated for four distinct models, demonstrating its power, e.g., to locate coexistence lines and the phase density at the coexistence.
Collapse
Affiliation(s)
- C E Fiore
- Departamento de Física, Universidade Federal do Paraná, 81531-980, Curitiba-PR, Brazil.
| | | |
Collapse
|
11
|
Fiore CE. Comparing different protocols of temperature selection in the parallel tempering method. J Chem Phys 2011; 135:114107. [DOI: 10.1063/1.3636382] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Badasyan AV, Tonoyan SA, Mamasakhlisov YS, Giacometti A, Benight AS, Morozov VF. Competition for hydrogen-bond formation in the helix-coil transition and protein folding. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2011; 83:051903. [PMID: 21728568 DOI: 10.1103/physreve.83.051903] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/15/2011] [Indexed: 05/31/2023]
Abstract
The problem of the helix-coil transition of biopolymers in explicit solvents, such as water, with the ability for hydrogen bonding with a solvent is addressed analytically using a suitably modified version of the Generalized Model of Polypeptide Chains. Besides the regular helix-coil transition, an additional coil-helix or reentrant transition is also found at lower temperatures. The reentrant transition arises due to competition between polymer-polymer and polymer-water hydrogen bonds. The balance between the two types of hydrogen bonding can be shifted to either direction through changes not only in temperature, but also by pressure, mechanical force, osmotic stress, or other external influences. Both polypeptides and polynucleotides are considered within a unified formalism. Our approach provides an explanation of the experimental difficulty of observing the reentrant transition with pressure and underscores the advantage of pulling experiments for studies of DNA. Results are discussed and compared with those reported in a number of recent publications with which a significant level of agreement is obtained.
Collapse
Affiliation(s)
- A V Badasyan
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca' Foscari Venezia, Venezia, Italy.
| | | | | | | | | | | |
Collapse
|
13
|
Almarza N, Noya E. Phase transitions of a lattice model for patchy particles with tetrahedral symmetry. Mol Phys 2011. [DOI: 10.1080/00268976.2010.523521] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
14
|
Fiore CE, da Luz MGE. A simple protocol for the probability weights of the simulated tempering algorithm: applications to first-order phase transitions. J Chem Phys 2010; 133:244102. [PMID: 21197971 DOI: 10.1063/1.3519813] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023] Open
Abstract
The simulated tempering (ST) is an important method to deal with systems whose phase spaces are hard to sample ergodically. However, it uses accepting probabilities weights, which often demand involving and time consuming calculations. Here it is shown that such weights are quite accurately obtained from the largest eigenvalue of the transfer matrix--a quantity straightforward to compute from direct Monte Carlo simulations--thus simplifying the algorithm implementation. As tests, different systems are considered, namely, Ising, Blume-Capel, Blume-Emery-Griffiths, and Bell-Lavis liquid water models. In particular, we address first-order phase transition at low temperatures, a regime notoriously difficulty to simulate because the large free-energy barriers. The good results found (when compared with other well established approaches) suggest that the ST can be a valuable tool to address strong first-order phase transitions, a possibility still not well explored in the literature.
Collapse
Affiliation(s)
- Carlos E Fiore
- Departamento de Física, Universidade Federal do Paraná, CP 19044, 81531-980 Curitiba-PR, Brazil.
| | | |
Collapse
|
15
|
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
|
16
|
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
|
17
|
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
|