1
|
Baran Ł, Rżysko W, MacDowell LG. Self-diffusion and shear viscosity for the TIP4P/Ice water model. J Chem Phys 2023; 158:064503. [PMID: 36792509 DOI: 10.1063/5.0134932] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
With an ever-increasing interest in water properties, many intermolecular force fields have been proposed to describe the behavior of water. Unfortunately, good models for liquid water usually cannot provide simultaneously an accurate melting point for ice. For this reason, the TIP4P/Ice model was developed for targeting the melting point and has become the preferred choice for simulating ice at coexistence. Unfortunately, available data for its dynamic properties in the liquid state are scarce. Therefore, we demonstrate a series of simulations aimed at the calculation of transport coefficients for the TIP4P/Ice model over a large range of thermodynamic conditions, ranging from T = 245 K to T = 350 K, for the temperature, and from p = 0 to p = 500 MPa, for the pressure. We have found that the self-diffusion (shear viscosity) exhibits smaller (increased) values than TIP4P/2005 and experiments. However, rescaling the temperature with respect to the triple point temperature, as in a corresponding states plot, we find that TIP4P/Ice compares very well with TIP4P/2005 and experiment. Such observations allow us to infer that despite the different original purposes of these two models examined here, one can benefit from a vast number of reports regarding the behavior of transport coefficients for the TIP4P/2005 model and utilize them following the routine described in this paper.
Collapse
Affiliation(s)
- Łukasz Baran
- Department of Theoretical Chemistry, Institute of Chemical Sciences, Faculty of Chemistry, Maria-Curie-Sklodowska University in Lublin, Pl. M Curie-Sklodowskiej 3, 20-031 Lublin, Poland
| | - Wojciech Rżysko
- Department of Theoretical Chemistry, Institute of Chemical Sciences, Faculty of Chemistry, Maria-Curie-Sklodowska University in Lublin, Pl. M Curie-Sklodowskiej 3, 20-031 Lublin, Poland
| | - Luis G MacDowell
- Departamento de Química-Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
2
|
Jakse N, Sandberg J, Granz LF, Saliou A, Jarry P, Devijver E, Voigtmann T, Horbach J, Meyer A. Machine learning interatomic potentials for aluminium: application to solidification phenomena. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 51:035402. [PMID: 36301702 DOI: 10.1088/1361-648x/ac9d7d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/06/2022] [Accepted: 10/25/2022] [Indexed: 06/16/2023]
Abstract
In studying solidification process by simulations on the atomic scale, the modeling of crystal nucleation or amorphization requires the construction of interatomic interactions that are able to reproduce the properties of both the solid and the liquid states. Taking into account rare nucleation events or structural relaxation under deep undercooling conditions requires much larger length scales and longer time scales than those achievable byab initiomolecular dynamics (AIMD). This problem is addressed by means of classical molecular dynamics simulations using a well established high dimensional neural network potential trained on a set of configurations generated by AIMD relevant for solidification phenomena. Our dataset contains various crystalline structures and liquid states at different pressures, including their time fluctuations in a wide range of temperatures. Applied to elemental aluminium, the resulting potential is shown to be efficient to reproduce the basic structural, dynamics and thermodynamic quantities in the liquid and undercooled states. Early stages of crystallization are further investigated on a much larger scale with one million atoms, allowing us to unravel features of the homogeneous nucleation mechanisms in the fcc phase at ambient pressure as well as in the bcc phase at high pressure with unprecedented accuracy close to theab initioone. In both cases, a single step nucleation process is observed.
Collapse
Affiliation(s)
- Noel Jakse
- Université Grenoble Alpes, CNRS, Grenoble INP, SIMaP, F-38000 Grenoble, France
| | - Johannes Sandberg
- Université Grenoble Alpes, CNRS, Grenoble INP, SIMaP, F-38000 Grenoble, France
- Institut für Materialphysik im Weltraum, Deutsches Zentrum für Luft- und Raumfahrt (DLR), 51170 Köln, Germany
- Department of Physics, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany
| | - Leon F Granz
- Institut für Materialphysik im Weltraum, Deutsches Zentrum für Luft- und Raumfahrt (DLR), 51170 Köln, Germany
- Department of Physics, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany
| | - Anthony Saliou
- Université Grenoble Alpes, CNRS, Grenoble INP, SIMaP, F-38000 Grenoble, France
| | - Philippe Jarry
- C-TEC, Parc Economique Centr'alp, 725 rue Aristide Bergès, CS10027, Voreppe 38341 CEDEX, France
| | - Emilie Devijver
- Université Grenoble Alpes, CNRS, Grenoble INP, LIG, F-38000 Grenoble, France
| | - Thomas Voigtmann
- Institut für Materialphysik im Weltraum, Deutsches Zentrum für Luft- und Raumfahrt (DLR), 51170 Köln, Germany
- Department of Physics, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany
| | - Jürgen Horbach
- Institut für Theoretische Physik II, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, 40225 Düsseldorf, Germany
| | - Andreas Meyer
- Institut für Materialphysik im Weltraum, Deutsches Zentrum für Luft- und Raumfahrt (DLR), 51170 Köln, Germany
- Institut Laue-Langevin (ILL), 38042 Grenoble, France
| |
Collapse
|
3
|
Attia E, Dyre JC, Pedersen UR. Extreme case of density scaling: The Weeks-Chandler-Andersen system at low temperatures. Phys Rev E 2021; 103:062140. [PMID: 34271644 DOI: 10.1103/physreve.103.062140] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Accepted: 05/28/2021] [Indexed: 11/07/2022]
Abstract
This paper studies numerically the Weeks-Chandler-Andersen system, which is shown to obey hidden scale invariance with a density-scaling exponent that varies from below 5 to above 500. This unprecedented variation makes it advantageous to use the fourth-order Runge-Kutta algorithm for tracing out isomorphs. Good isomorph invariance of structure and dynamics is observed over more than three orders of magnitude temperature variation. For all state points studied, the virial potential-energy correlation coefficient and the density-scaling exponent are controlled mainly by the temperature. Based on the assumption of statistically independent pair interactions, a mean-field theory is developed that rationalizes this finding and provides an excellent fit to data at low temperatures.
Collapse
Affiliation(s)
- Eman Attia
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, 4000 Roskilde, Denmark
| | - Jeppe C Dyre
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, 4000 Roskilde, Denmark
| | - Ulf R Pedersen
- Glass and Time, IMFUFA, Department of Science and Environment, Roskilde University, P.O. Box 260, 4000 Roskilde, Denmark
| |
Collapse
|
4
|
Rozas RE, MacDowell LG, Toledo PG, Horbach J. Crystal growth of bcc titanium from the melt and interfacial properties: A molecular dynamics simulation study. J Chem Phys 2021; 154:184704. [PMID: 34241033 DOI: 10.1063/5.0049131] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
The crystal growth kinetics and interfacial properties of titanium (Ti) are studied using molecular dynamics computer simulation. The interactions between the Ti atoms are modeled via an embedded atom method potential. First, the free solidification method (FSM) is used to determine the melting temperature Tm at zero pressure where the transition from liquid to body-centered cubic crystal occurs. From the simulations with the FSM, the kinetic growth coefficients are also determined for different orientations of the crystal, analyzing how the coupling to the thermostat affects the estimates of the growth coefficients. At Tm, anisotropic interfacial stiffnesses and free energies as well as kinetic growth coefficients are determined from capillary wave fluctuations. The so-obtained growth coefficients from equilibrium fluctuations and without the coupling of the system to a thermostat agree well with those extracted from the FSM calculations.
Collapse
Affiliation(s)
- R E Rozas
- Department of Physics, University of Bío-Bío, Av. Collao 1202, P.O. Box 5-C, Concepción, Chile
| | - L G MacDowell
- Departamento de Química Física, Universidad Complutense de Madrid, Madrid 28040, Spain
| | - P G Toledo
- Chemical Engineering Department and Surface Analysis Laboratory (ASIF), University of Concepción, P.O. Box 160-C, Correo 3, Concepción, Chile
| | - J Horbach
- Institut für Theoretische Physik II: Weiche Materie, Heinrich-Heine-Universität Düsseldorf, Universitätsstraße 1, D-40225 Düsseldorf, Germany
| |
Collapse
|
5
|
Yurchenko SO, Komarov KA, Kryuchkov NP, Zaytsev KI, Brazhkin VV. Bizarre behavior of heat capacity in crystals due to interplay between two types of anharmonicities. J Chem Phys 2018; 148:134508. [DOI: 10.1063/1.5022969] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Stanislav O. Yurchenko
- Bauman Moscow State Technical University, 2nd Baumanskaya Street 5, Moscow 105005, Russia
| | - Kirill A. Komarov
- Bauman Moscow State Technical University, 2nd Baumanskaya Street 5, Moscow 105005, Russia
| | - Nikita P. Kryuchkov
- Bauman Moscow State Technical University, 2nd Baumanskaya Street 5, Moscow 105005, Russia
| | - Kirill I. Zaytsev
- Bauman Moscow State Technical University, 2nd Baumanskaya Street 5, Moscow 105005, Russia
| | - Vadim V. Brazhkin
- Institute for High Pressure Physics RAS, Kaluzhskoe Shosse, 14, Troitsk, Moscow 108840, Russia
| |
Collapse
|
6
|
Reinhart WF, Panagiotopoulos AZ. Crystal growth kinetics of triblock Janus colloids. J Chem Phys 2018; 148:124506. [DOI: 10.1063/1.5021347] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Wesley F. Reinhart
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | | |
Collapse
|
7
|
Şopu D, Rogal J, Drautz R. Thermodynamic and kinetic solid-liquid interface properties from transition path sampling. J Chem Phys 2016; 145:244703. [DOI: 10.1063/1.4972583] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
|
8
|
Glass-Forming Tendency of Molecular Liquids and the Strength of the Intermolecular Attractions. Sci Rep 2016; 6:36934. [PMID: 27883011 PMCID: PMC5121653 DOI: 10.1038/srep36934] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2016] [Accepted: 10/24/2016] [Indexed: 12/24/2022] Open
Abstract
When we cool down a liquid below the melting temperature, it can either crystallize or become supercooled, and then form a disordered solid called glass. Understanding what makes a liquid to crystallize readily in one case and form a stable glass in another is a fundamental problem in science and technology. Here we show that the crystallization/glass-forming tendencies of the molecular liquids might be correlated with the strength of the intermolecular attractions, as determined from the combined experimental and computer simulation studies. We use van der Waals bonded propylene carbonate and its less polar structural analog 3-methyl-cyclopentanone to show that the enhancement of the dipole-dipole forces brings about the better glass-forming ability of the sample when cooling from the melt. Our finding was rationalized by the mismatch between the optimal temperature range for the nucleation and crystal growth, as obtained for a modeled Lennard-Jones system with explicitly enhanced or weakened attractive part of the intermolecular 6–12 potential.
Collapse
|
9
|
Abstract
Although the freezing of liquids and melting of crystals are fundamental for many areas of the sciences, even simple properties like the temperature–pressure relation along the melting line cannot be predicted today. Here we present a theory in which properties of the coexisting crystal and liquid phases at a single thermodynamic state point provide the basis for calculating the pressure, density and entropy of fusion as functions of temperature along the melting line, as well as the variation along this line of the reduced crystalline vibrational mean-square displacement (the Lindemann ratio), and the liquid's diffusion constant and viscosity. The framework developed, which applies for the sizable class of systems characterized by hidden scale invariance, is validated by computer simulations of the standard 12-6 Lennard-Jones system. Melting is a classic first-order phase transition, but an accurate thermodynamic description is still lacking. Here, Pedersen et al. develop a theory, validated by simulations of the Lennard-Jones system, for the melting thermodynamics applicable to all systems characterized by hidden scale invariance.
Collapse
|