1
|
Lascaris E, Marchese F, Gaspar N. Crystallization and the liquid-liquid critical point in nonbonded modified-WAC models. J Chem Phys 2024; 161:044503. [PMID: 39037140 DOI: 10.1063/5.0215601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2024] [Accepted: 07/03/2024] [Indexed: 07/23/2024] Open
Abstract
For decades, it has been known that Liquid-Liquid Critical Points (LLCPs) can exist in one-component liquids, yet a comprehensive understanding of the conditions under which they arise remains elusive. To better comprehend the possible interplay between the LLCP and the crystalline phase, we conduct molecular dynamics simulations using the nonbonded family of modified-WAC (mWAC) models, which are known to exhibit a LLCP for certain parameter values. By comparing different versions of the mWAC model-those featuring a LLCP and those lacking one-we identify several key differences between the models relating to crystallization. Those models that do have a LLCP are found to have multiple stable crystalline phases, one of them being a solid-state ionic conductor similar to superionic ice. Moreover, we find that for models that do not have a LLCP, the liquid becomes a glass at a larger range of temperatures, possibly preventing the occurrence of a LLCP. Further studies are required to determine if these results are general or model-specific.
Collapse
Affiliation(s)
- Erik Lascaris
- Department of Chemistry & Physical Sciences, Pace University, New York, New York 10038, USA
| | - Francesca Marchese
- Department of Chemistry & Physical Sciences, Pace University, New York, New York 10038, USA
| | - Nicole Gaspar
- Department of Chemistry & Physical Sciences, Pace University, New York, New York 10038, USA
| |
Collapse
|
2
|
Dubey V, Dueby S, Daschakraborty S. Breakdown of the Stokes-Einstein relation in supercooled water: the jump-diffusion perspective. Phys Chem Chem Phys 2021; 23:19964-19986. [PMID: 34515269 DOI: 10.1039/d1cp02202d] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Although water is the most ubiquitous liquid it shows many thermodynamic and dynamic anomalies. Some of the anomalies further intensify in the supercooled regime. While many experimental and theoretical studies have focused on the thermodynamic anomalies of supercooled water, fewer studies explored the dynamical anomalies very extensively. This is due to the intricacy of the experimental measurement of the dynamical properties of supercooled water. Violation of the Stokes-Einstein relation (SER), an important relation connecting the diffusion of particles with the viscosity of the medium, is one of the major dynamical anomalies. In absence of experimentally measured viscosity, researchers used to check the validity of SER indirectly using average translational relaxation time or α-relaxation time. Very recently, the viscosity of supercooled water was accurately measured at a wide range of temperatures and pressures. This allowed direct verification of the SER at different temperature-pressure thermodynamic state points. An increasing breakdown of the SER was observed with decreasing temperature. Increasing pressure reduces the extent of breakdown. Although some well-known theories explained the above breakdown, a detailed molecular mechanism was still elusive. Recently, a translational jump-diffusion (TJD) approach has been able to quantitatively explain the breakdown of the SER in pure supercooled water and an aqueous solution of methanol. The objective of this article is to present a detailed and state-of-the-art analysis of the past and present works on the breakdown of SER in supercooled water with a specific focus on the new TJD approach for explaining the breakdown of the SER.
Collapse
Affiliation(s)
- Vikas Dubey
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | - Shivam Dueby
- Department of Chemistry, Indian Institute of Technology Patna, Bihar 801106, India.
| | | |
Collapse
|
3
|
Yang M, Karmakar T, Parrinello M. Liquid-Liquid Critical Point in Phosphorus. PHYSICAL REVIEW LETTERS 2021; 127:080603. [PMID: 34477397 DOI: 10.1103/physrevlett.127.080603] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/04/2021] [Revised: 07/07/2021] [Accepted: 07/19/2021] [Indexed: 06/13/2023]
Abstract
The study of liquid-liquid phase transitions has attracted considerable attention. One interesting example of this phenomenon is phosphorus, for which the existence of a first-order phase transition between a low density insulating molecular phase and a conducting polymeric phase has been experimentally established. In this Letter, we model this transition by an ab initio quality molecular dynamics simulation and explore a large portion of the liquid section of the phase diagram. We draw the liquid-liquid coexistence curve and determine that it terminates into a second-order critical point. Close to the critical point, large coupled structure and electronic structure fluctuations are observed.
Collapse
Affiliation(s)
- Manyi Yang
- Italian Institute of Technology, Via Melen 83, 16152 Genova, Italy
| | - Tarak Karmakar
- Italian Institute of Technology, Via Melen 83, 16152 Genova, Italy
| | | |
Collapse
|
4
|
Affiliation(s)
- Hajime Tanaka
- Department of Fundamental Engineering, Institute of Industrial Science, University of Tokyo, Meguro-ku, Tokyo 153-8505, Japan
| |
Collapse
|
5
|
Pal T, Vogel M. On the relevance of electrostatic interactions for the structural relaxation of ionic liquids: A molecular dynamics simulation study. J Chem Phys 2019; 150:124501. [DOI: 10.1063/1.5085508] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Affiliation(s)
- Tamisra Pal
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstraße 6, 64289 Darmstadt, Germany
| | - Michael Vogel
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstraße 6, 64289 Darmstadt, Germany
| |
Collapse
|
6
|
Hestand NJ, Skinner JL. Perspective: Crossing the Widom line in no man’s land: Experiments, simulations, and the location of the liquid-liquid critical point in supercooled water. J Chem Phys 2018; 149:140901. [DOI: 10.1063/1.5046687] [Citation(s) in RCA: 54] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Nicholas J. Hestand
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| | - J. L. Skinner
- Institute for Molecular Engineering, University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
7
|
Palmer JC, Poole PH, Sciortino F, Debenedetti PG. Advances in Computational Studies of the Liquid–Liquid Transition in Water and Water-Like Models. Chem Rev 2018; 118:9129-9151. [DOI: 10.1021/acs.chemrev.8b00228] [Citation(s) in RCA: 116] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jeremy C. Palmer
- Department of Chemical and Biomolecular Engineering, University of Houston, Houston, Texas 77204, United States
| | - Peter H. Poole
- Department of Physics, St. Francis Xavier University, Antigonish, NS B2G 2W5, Canada
| | - Francesco Sciortino
- Dipartimento di Fisica and CNR-ISC, Sapienza Universita’ di Roma, Piazzale A. Moro 5, 00185 Rome, Italy
| | - Pablo G. Debenedetti
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, United States
| |
Collapse
|
8
|
Pafong Sanjon E, Drossel B, Vogel M. Effects of the bond polarity on the structural and dynamical properties of silica-like liquids. J Chem Phys 2018; 148:104506. [PMID: 29544292 DOI: 10.1063/1.5017681] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Silica is a network-forming liquid that shares many properties with water due to its tetrahedral structure. It undergoes a transition from a fragile to a strong liquid as the temperature is decreased, which is accompanied by a structural change to lower density and higher tetrahedral order. In order to disentangle the effects of Coulomb and van der Waals interactions on the structure and dynamics of liquid silica, we modify the bond polarity by changing the partial charges assigned to each atom. Using molecular dynamics simulations, we show that density, tetrahedral order, and structural relaxation times decrease when reducing bond polarity. Moreover, we find that the density maximum and the fragile-to-strong transition move to lower temperatures until they eventually vanish when the partial charges are decreased below approximately 75% of their regular value. Irrespective of whether strong or fragile behavior exists, structural relaxation is governed by hopping motion at sufficiently low temperatures. As long as there is a strong regime, the energy barrier associated with strong dynamics decreases with decreasing partial charges, but the dependence on the bond polarity differs from that of the activation energy in the Arrhenius regime at high temperatures. We show that the fragile-to-strong transition is associated with structural changes occurring between the first and second coordination shells that lead to a decrease in density and an increase in tetrahedral order. In particular, independent of the value of the partial charges, the distribution of the local structures is the same at this dynamic crossover, but we find no evidence that the effect occurs upon crossing the Widom line. In the fragile regime at intermediate temperatures, the relaxation times are well described by a previously proposed model which decomposes the apparent activation energy into a constant single-particle contribution and a temperature-dependent collective contribution. However, our results for silica-like melts do not obey several common relations of the model parameters reported for molecular glass formers.
Collapse
Affiliation(s)
- E Pafong Sanjon
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| | - B Drossel
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| | - M Vogel
- Institut für Festkörperphysik, Technische Universität Darmstadt, Hochschulstr. 6, 64289 Darmstadt, Germany
| |
Collapse
|
9
|
Zhao G, Wang H, Hu DM, Ding MC, Zhao XG, Yan JL. Anomalous phase behavior of first-order fluid-liquid phase transition in phosphorus. J Chem Phys 2018; 147:204501. [PMID: 29195280 DOI: 10.1063/1.4999009] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Although the existence of liquid-liquid phase transition has become more and more convincing, whether it will terminate at a critical point and what is the order parameter are still open. To explore these questions, we revisit the fluid-liquid phase transition (FLPT) in phosphorus (P) and study its phase behavior by performing extensive first-principles molecular dynamics simulations. The FLPT observed in experiments is well reproduced, and a fluid-liquid critical point (FLCP) at T = 3000 ∼ 3500 K, P = 1.5-2.0 Kbar is found. With decreasing temperature from the FLCP along the transition line, the density difference (Δρ) between two coexisting phases first increases from zero and then anomalously decreases; however, the entropy difference (ΔS) continuously increases from zero. These features suggest that an order parameter containing contributions from both the density and the entropy is needed to describe the FLPT in P, and at least at low temperatures, the entropy, instead of the density, governs the FLPT.
Collapse
Affiliation(s)
- G Zhao
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| | - H Wang
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| | - D M Hu
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| | - M C Ding
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| | - X G Zhao
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| | - J L Yan
- School of Physics and Optoelectronic Engineering, Ludong University, Yantai 264025, People's Republic of China
| |
Collapse
|
10
|
Affiliation(s)
- Jingxiang Guo
- Department of Chemical and Biomolecular Engineering, University of Houston , Houston, TX, USA
| | - Rakesh S. Singh
- Department of Chemistry, Johns Hopkins University , Baltimore, MD, USA
| | - Jeremy C. Palmer
- Department of Chemical and Biomolecular Engineering, University of Houston , Houston, TX, USA
| |
Collapse
|
11
|
Impact of local symmetry breaking on the physical properties of tetrahedral liquids. Proc Natl Acad Sci U S A 2018; 115:1980-1985. [PMID: 29439201 DOI: 10.1073/pnas.1717233115] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Water and silica are the most important materials with local tetrahedral symmetry. They have similar crystalline polymorphs and exhibit anomalous density maximum in the liquid state. However, water and silica also show very different characteristics. For instance, the density of water varies much more sharply than that of liquid silica near the maximum as temperature changes. More notably, silica is a very good glass-former, but water is an extremely poor one. The physical origins of these similarities and differences still remain elusive, due to the lack of a microscopic understanding of the structural ordering in these two important liquids. Here, by accessing microscopic structural information by computer simulations, we reveal that local translational symmetry breaking is responsible for the density anomalies. On the other hand, the difference in the degree of local orientational symmetry breaking between water and silica, which originates from the difference in their bonding nature, causes not only the difference in the sharpness of density anomalies, but also their distinct glass-forming abilities. Our work not only shows the crucial roles of local translational and orientational symmetry breaking in the physical properties of the two extremely important materials, water and silica, but also provides a unified scenario applicable for other tetrahedral liquids such as Si, Ge, C, BeF2, and GeO2.
Collapse
|
12
|
Guo J, Palmer JC. Fluctuations near the liquid–liquid transition in a model of silica. Phys Chem Chem Phys 2018; 20:25195-25202. [DOI: 10.1039/c8cp04237c] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Molecular dynamics simulations reveal anomalous small-angle scattering and liquid–liquid phase separation in an ionic model of silica.
Collapse
Affiliation(s)
- Jingxiang Guo
- Department of Chemical and Biomolecular Engineering
- University of Houston
- Houston
- USA
| | - Jeremy C. Palmer
- Department of Chemical and Biomolecular Engineering
- University of Houston
- Houston
- USA
| |
Collapse
|
13
|
Chen R, Lascaris E, Palmer JC. Liquid–liquid phase transition in an ionic model of silica. J Chem Phys 2017. [DOI: 10.1063/1.4984335] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
14
|
Gallo P, Amann-Winkel K, Angell CA, Anisimov MA, Caupin F, Chakravarty C, Lascaris E, Loerting T, Panagiotopoulos AZ, Russo J, Sellberg JA, Stanley HE, Tanaka H, Vega C, Xu L, Pettersson LGM. Water: A Tale of Two Liquids. Chem Rev 2016; 116:7463-500. [PMID: 27380438 PMCID: PMC5424717 DOI: 10.1021/acs.chemrev.5b00750] [Citation(s) in RCA: 447] [Impact Index Per Article: 55.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
Water is the most abundant liquid
on earth and also the substance
with the largest number of anomalies in its properties. It is a prerequisite
for life and as such a most important subject of current research
in chemical physics and physical chemistry. In spite of its simplicity
as a liquid, it has an enormously rich phase diagram where different
types of ices, amorphous phases, and anomalies disclose a path that
points to unique thermodynamics of its supercooled liquid state that
still hides many unraveled secrets. In this review we describe the
behavior of water in the regime from ambient conditions to the deeply
supercooled region. The review describes simulations and experiments
on this anomalous liquid. Several scenarios have been proposed to
explain the anomalous properties that become strongly enhanced in
the supercooled region. Among those, the second critical-point scenario
has been investigated extensively, and at present most experimental
evidence point to this scenario. Starting from very low temperatures,
a coexistence line between a high-density amorphous phase and a low-density
amorphous phase would continue in a coexistence line between a high-density
and a low-density liquid phase terminating in a liquid–liquid
critical point, LLCP. On approaching this LLCP from the one-phase
region, a crossover in thermodynamics and dynamics can be found. This
is discussed based on a picture of a temperature-dependent balance
between a high-density liquid and a low-density liquid favored by,
respectively, entropy and enthalpy, leading to a consistent picture
of the thermodynamics of bulk water. Ice nucleation is also discussed,
since this is what severely impedes experimental investigation of
the vicinity of the proposed LLCP. Experimental investigation of stretched
water, i.e., water at negative pressure, gives access to a different
regime of the complex water diagram. Different ways to inhibit crystallization
through confinement and aqueous solutions are discussed through results
from experiments and simulations using the most sophisticated and
advanced techniques. These findings represent tiles of a global picture
that still needs to be completed. Some of the possible experimental
lines of research that are essential to complete this picture are
explored.
Collapse
Affiliation(s)
- Paola Gallo
- Dipartimento di Matematica e Fisica, Università Roma Tre , Via della Vasca Navale 84, 00146 Rome, Italy
| | - Katrin Amann-Winkel
- Department of Physics, AlbaNova University Center, Stockholm University , SE-106 91 Stockholm, Sweden
| | - Charles Austen Angell
- Department of Chemistry and Biochemistry, Arizona State University , Tempe, Arizona 85287, United States
| | - Mikhail Alexeevich Anisimov
- Institute for Physical Science and Technology and Department of Chemical and Biomolecular Engineering, University of Maryland , College Park, Maryland 20742, United States
| | - Frédéric Caupin
- Institut Lumière Matière, UMR5306 Université Claude Bernard Lyon 1-CNRS, Université de Lyon, Institut Universitaire de France , 69622 Villeurbanne, France
| | - Charusita Chakravarty
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas , New Delhi 110016, India
| | - Erik Lascaris
- Center for Polymer Studies and Department of Physics, Boston University , Boston, Massachusetts 02215, United States
| | - Thomas Loerting
- Institute of Physical Chemistry, University of Innsbruck , 6020 Innsbruck, Austria
| | | | - John Russo
- Institute of Industrial Science, University of Tokyo , 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan.,School of Mathematics, University of Bristol , Bristol BS8 1TW, United Kingdom
| | - Jonas Alexander Sellberg
- Biomedical and X-ray Physics, Department of Applied Physics, AlbaNova University Center, KTH Royal Institute of Technology , SE-106 91 Stockholm, Sweden
| | - Harry Eugene Stanley
- Center for Polymer Studies and Department of Physics, Boston University , Boston, Massachusetts 02215, United States
| | - Hajime Tanaka
- Institute of Industrial Science, University of Tokyo , 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan
| | - Carlos Vega
- Departamento de Quimica Fisica, Facultad de Ciencias Quimicas, Universidad Complutense de Madrid , 28040 Madrid, Spain
| | - Limei Xu
- International Centre for Quantum Materials and School of Physics, Peking University , Beijing 100871, China.,Collaborative Innovation Center of Quantum Matter, Beijing 100871, China
| | | |
Collapse
|