1
|
Bulutoglu PS, Wang S, Boukerche M, Nere NK, Corti DS, Ramkrishna D. An investigation of the kinetics and thermodynamics of NaCl nucleation through composite clusters. PNAS NEXUS 2022; 1:pgac033. [PMID: 36713321 PMCID: PMC9802385 DOI: 10.1093/pnasnexus/pgac033] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/24/2021] [Revised: 03/04/2022] [Accepted: 03/28/2022] [Indexed: 02/01/2023]
Abstract
Having a good understanding of nucleation is critical for the control of many important processes, such as polymorph selection during crystallization. However, a complete picture of the molecular-level mechanisms of nucleation remains elusive. In this work, we take an in-depth look at the NaCl homogeneous nucleation mechanism through thermodynamics. Distinguished from the classical nucleation theory, we calculate the free energy of nucleation as a function of two nucleus size coordinates: crystalline and amorphous cluster sizes. The free energy surface reveals a thermodynamic preference for a nonclassical mechanism of nucleation through a composite cluster, where the crystalline nucleus is surrounded by an amorphous layer. The thickness of the amorphous layer increases with an increase in supersaturation. The computed free energy landscape agrees well with the composite cluster-free energy model, through which phase specific thermodynamic properties are evaluated. As the supersaturation increases, there is a change in stability of the amorphous phase relative to the solution phase, resulting in a change from one-step to two-step mechanism, seen clearly from the free energy profile along the minimum free energy path crossing the transition curve. By obtaining phase-specific diffusion coefficients, we construct the full mesoscopic model and present a clear roadmap for NaCl nucleation.
Collapse
Affiliation(s)
- Pelin S Bulutoglu
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN 47907-2100, USA
| | - Shiyan Wang
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN 47907-2100, USA
| | - Moussa Boukerche
- Process Research and Development , AbbVie Inc, North Chicago, IL 60064, USA
| | - Nandkishor K Nere
- Process Research and Development , AbbVie Inc, North Chicago, IL 60064, USA
| | - David S Corti
- Davidson School of Chemical Engineering, Purdue University, West Lafayette, IN 47907-2100, USA
| | | |
Collapse
|
2
|
Nie C, Geng J, Marlow WH. Formation free energy of an i-mer at spinodal. J Chem Phys 2021; 154:234108. [PMID: 34241258 DOI: 10.1063/5.0051879] [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
In statistical mechanics, the formation free energy of an i-mer can be understood as the Gibbs free energy change in a system consisting of pure monomers after and prior to the formation of the i-mer. For molecules interacting via Lennard-Jones potential, we have computed the formation free energy of a Stillinger i-mer [F. H. Stillinger, J. Chem. Phys. 38, 1486 (1963)] and a ten Wolde-Frenkel (tWF) [P. R. ten Wolde and D. Frenkel, J. Chem. Phys. 109, 9901 (1998)] i-mer at spinodal at reduced temperatures from 0.7 to 1.2. It turns out that the size of a critical Stillinger i-mer remains finite and its formation free energy is on the order of kBT, and the size of a critical tWF i-mer remains finite and its formation free energy is even higher. This can be explained by Binder's theory [K. Binder, Phys. Rev. A 29, 341 (1984)] that for a system, when approaching spinodal, if the Ginzburg criterion is not satisfied, a gradual transition will take place from nucleation to spinodal decomposition, where the free-energy barrier height is on the order of kBT.
Collapse
Affiliation(s)
- Chu Nie
- School of Power and Mechanical Engineering, Wuhan University, Wuhan 430072, China
| | - Jun Geng
- State Power Investment Cooperation Hydrogen Energy Development Co., Ltd., South Park, Bejing Future Science & Technology Park, Beijing 100029, China
| | - William H Marlow
- Nuclear Engineering Department, Texas A&M University, 3133 TAMU, College Station, Texas 77843-3133, USA
| |
Collapse
|
3
|
Extracting meaningful standard enthalpies and entropies of activation for surface reactions from kinetic rates. REACTION KINETICS MECHANISMS AND CATALYSIS 2020. [DOI: 10.1007/s11144-020-01747-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
4
|
Nie C, Geng J, Marlow WH. Formation free energies of clusters at high supersaturations. J Chem Phys 2019; 151:134111. [DOI: 10.1063/1.5111943] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Affiliation(s)
- Chu Nie
- School of Power and Mechanical Engineering, Wuhan University, Wuhan 430072, China
| | - Jun Geng
- State Power Investment Corporation Hydrogen Energy Development Co. Ltd., South Park, Bejing Future Science & Technology Park, Beijing 100029, China
| | - William H. Marlow
- Nuclear Engineering Department, Texas A&M University, 3133 TAMU, College Station, Texas 77843-3133, USA
| |
Collapse
|
5
|
Rehner P, Gross J. Surface tension of droplets and Tolman lengths of real substances and mixtures from density functional theory. J Chem Phys 2018; 148:164703. [DOI: 10.1063/1.5020421] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Affiliation(s)
- Philipp Rehner
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| | - Joachim Gross
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart, Pfaffenwaldring 9, 70569 Stuttgart, Germany
| |
Collapse
|
6
|
Xiao Q, Liu Y, Guo Z, Liu Z, Frenkel D, Dobnikar J, Zhang X. What experiments on pinned nanobubbles can tell about the critical nucleus for bubble nucleation. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2017; 40:114. [PMID: 29273950 DOI: 10.1140/epje/i2017-11604-7] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/15/2017] [Accepted: 12/04/2017] [Indexed: 05/19/2023]
Abstract
The process of homogeneous bubble nucleation is almost impossible to probe experimentally, except near the critical point or for liquids under large negative tension. Elsewhere in the phase diagram, the bubble nucleation barrier is so high as to be effectively insurmountable. Consequently, there is a severe lack of experimental studies of homogenous bubble nucleation under conditions of practical importance (e.g., cavitation). Here we use a simple geometric relation to show that we can obtain information about the homogeneous nucleation process from Molecular Dynamics studies of bubble formation in solvophobic nanopores on a solid surface. The free energy of pinned nanobubbles has two extrema as a function of volume: one state corresponds to a free-energy maximum ("the critical nucleus"), the other corresponds to a free-energy minimum (the metastable, pinned nanobubble). Provided that the surface tension does not depend on nanobubble curvature, the radius of the curvature of the metastable surface nanobubble is independent of the radius of the pore and is equal to the radius of the critical nucleus in homogenous bubble nucleation. This observation opens the way to probe the parameters that determine homogeneous bubble nucleation under experimentally accessible conditions, e.g. with AFM studies of metastable nanobubbles. Our theoretical analysis also indicates that a surface with pores of different sizes can be used to determine the curvature corrections to the surface tension. Our conclusions are not limited to bubble nucleation but suggest that a similar approach could be used to probe the structure of critical nuclei in crystal nucleation.
Collapse
Affiliation(s)
- Qianxiang Xiao
- State Key Laboratory of Organic-Inorganic Composites, Beijing University of Chemical Technology, 100029, Beijing, China
| | - Yawei Liu
- State Key Laboratory of Organic-Inorganic Composites, Beijing University of Chemical Technology, 100029, Beijing, China
| | - Zhenjiang Guo
- State Key Laboratory of Organic-Inorganic Composites, Beijing University of Chemical Technology, 100029, Beijing, China
| | - Zhiping Liu
- State Key Laboratory of Organic-Inorganic Composites, Beijing University of Chemical Technology, 100029, Beijing, China
| | - Daan Frenkel
- Department of Chemistry, University of Cambridge, Lensfield Road, CB21EW, Cambridge, UK
| | - Jure Dobnikar
- Institute of Physics, Chinese Academy of Sciences, 100190, Beijing, China.
| | - Xianren Zhang
- State Key Laboratory of Organic-Inorganic Composites, Beijing University of Chemical Technology, 100029, Beijing, China
| |
Collapse
|
7
|
Urrutia I, Paganini IE. Virial series for inhomogeneous fluids applied to the Lennard-Jones wall-fluid surface tension at planar and curved walls. J Chem Phys 2016; 144:174102. [DOI: 10.1063/1.4947587] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Ignacio Urrutia
- Departamento de Física de la Materia Condensada, Centro Atómico Constituyentes, CNEA, Av.Gral. Paz 1499, 1650 Pcia. de Buenos Aires, Argentina and CONICET, Avenida Rivadavia 1917, C1033AAJ Buenos Aires, Argentina
| | - Iván E. Paganini
- Departamento de Física de la Materia Condensada, Centro Atómico Constituyentes, CNEA, Av.Gral. Paz 1499, 1650 Pcia. de Buenos Aires, Argentina and CONICET, Avenida Rivadavia 1917, C1033AAJ Buenos Aires, Argentina
| |
Collapse
|
8
|
Bourasseau E, Malfreyt P, Ghoufi A. Surface tension and long range corrections of cylindrical interfaces. J Chem Phys 2015; 143:234708. [DOI: 10.1063/1.4937924] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
|
9
|
Wilhelmsen Ø, Bedeaux D, Kjelstrup S, Reguera D. Communication: superstabilization of fluids in nanocontainers. J Chem Phys 2014; 141:071103. [PMID: 25149768 DOI: 10.1063/1.4893701] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
One of the main challenges of thermodynamics is to predict and measure accurately the properties of metastable fluids. Investigation of these fluids is hindered by their spontaneous transformation by nucleation into a more stable phase. We show how small closed containers can be used to completely prevent nucleation, achieving infinitely long-lived metastable states. Using a general thermodynamic framework, we derive simple formulas to predict accurately the conditions (container sizes) at which this superstabilization takes place and it becomes impossible to form a new stable phase. This phenomenon opens the door to control nucleation of deeply metastable fluids at experimentally feasible conditions, having important implications in a wide variety of fields.
Collapse
Affiliation(s)
- Øivind Wilhelmsen
- Department of Chemistry, Norwegian University of Science and Technology, NO-7391 Trondheim, Norway
| | - Dick Bedeaux
- Department of Chemistry, Norwegian University of Science and Technology, NO-7391 Trondheim, Norway
| | - Signe Kjelstrup
- Department of Chemistry, Norwegian University of Science and Technology, NO-7391 Trondheim, Norway
| | - David Reguera
- Departament de Física Fonamental, Universitat de Barcelona, Martí i Franquès 1, Barcelona, Spain
| |
Collapse
|
10
|
Joswiak MN, Duff N, Doherty MF, Peters B. Size-Dependent Surface Free Energy and Tolman-Corrected Droplet Nucleation of TIP4P/2005 Water. J Phys Chem Lett 2013; 4:4267-72. [PMID: 26296177 DOI: 10.1021/jz402226p] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
Classical nucleation theory is notoriously inaccurate when using the macroscopic surface free energy for a planar interface. We examine the size dependence of the surface free energy for TIP4P/2005 water nanodroplets (radii ranging from 0.7 to 1.6 nm) at 300 K with the mitosis method, that is, by reversibly splitting the droplets into two subclusters. We calculate the Tolman length to be -0.56 ± 0.09 Å, which indicates that the surface free energy of water droplets that we investigated is 5-11 mJ/m(2) greater than the planar surface free energy. We incorporate the computed Tolman length into a modified classical nucleation theory (δ-CNT) and obtain modified expressions for the critical nucleus size and barrier height. δ-CNT leads to excellent agreement with independently measured nucleation kinetics.
Collapse
Affiliation(s)
- Mark N Joswiak
- †Department of Chemical Engineering and ‡Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| | - Nathan Duff
- †Department of Chemical Engineering and ‡Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| | - Michael F Doherty
- †Department of Chemical Engineering and ‡Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| | - Baron Peters
- †Department of Chemical Engineering and ‡Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
11
|
Zhukhovitskii DI. Molecular dynamics study of nanobubbles in the equilibrium Lennard-Jones fluid. J Chem Phys 2013; 139:164513. [DOI: 10.1063/1.4826648] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
12
|
Torabi K, Corti DS. Toward a molecular theory of homogeneous bubble nucleation: II. Calculation of the number density of critical nuclei and the rate of nucleation. J Phys Chem B 2013; 117:12491-504. [PMID: 24020901 DOI: 10.1021/jp404151h] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
In the present paper, we develop a method to calculate the rate of homogeneous bubble nucleation within a superheated L-J liquid based on the (n,v) equilibrium embryo free energy surface introduced in the first paper (DOI: 10.1021/jp404149n). We express the nucleation rate as the product of the concentration of critical nuclei within the metastable liquid phase and the relevant forward rate coefficient. We calculate the forward rate coefficient of the critical nuclei from their average lifetime as determined from MD simulations of a large number of embryo trajectories initiated from the transitional region of the metastable liquid configuration space. Therefore, the proposed rate coefficient does not rely on any predefined reaction coordinate. In our model, the critical nuclei belong to the region of the configuration space where the committor probability is about one-half, guaranteeing the dynamical relevance of the proposed embryos. One novel characteristic of our approach is that we define a limit for the configuration space of the equilibrium metastable phase and do not include the configurations that have zero committor probability in the nucleation free energy surface. Furthermore, in order to take into account the transitional degrees of freedom of the critical nuclei, we develop a simulation-based approach for rigorously mapping the free energy of the (n,v) equilibrium embryos to the concentration of the critical nuclei within the bulk metastable liquid phase.
Collapse
Affiliation(s)
- Korosh Torabi
- School of Chemical Engineering, Purdue University , West Lafayette, Indiana 47907-2100, United States
| | | |
Collapse
|
13
|
Kusaka I, Talreja M, Tomasko DL. Beyond classical theory: Predicting the free energy barrier of bubble nucleation in polymer foaming. AIChE J 2013. [DOI: 10.1002/aic.14062] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Isamu Kusaka
- William G. Lowrie Dept. of Chemical and Biomolecular Engineering; The Koffolt Laboratories, The Ohio State University; Columbus OH 43210
| | - Manish Talreja
- School of Chemical and Biomolecular Engineering; Georgia Institute of Technology; NW, Atlanta GA 30332
| | - David L. Tomasko
- William G. Lowrie Dept. of Chemical and Biomolecular Engineering; The Koffolt Laboratories, The Ohio State University; Columbus OH 43210
| |
Collapse
|
14
|
Malijevský A, Jackson G. A perspective on the interfacial properties of nanoscopic liquid drops. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2012; 24:464121. [PMID: 23114181 DOI: 10.1088/0953-8984/24/46/464121] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
The structural and interfacial properties of nanoscopic liquid drops are assessed by means of mechanical, thermodynamical, and statistical mechanical approaches that are discussed in detail, including original developments at both the macroscopic level and the microscopic level of density functional theory (DFT). With a novel analysis we show that a purely macroscopic (static) mechanical treatment can lead to a qualitatively reasonable description of the surface tension and the Tolman length of a liquid drop; the latter parameter, which characterizes the curvature dependence of the tension, is found to be negative and has a magnitude of about a half of the molecular dimension. A mechanical slant cannot, however, be considered satisfactory for small finite-size systems where fluctuation effects are significant. From the opposite perspective, a curvature expansion of the macroscopic thermodynamic properties (density and chemical potential) is then used to demonstrate that a purely thermodynamic approach of this type cannot in itself correctly account for the curvature correction of the surface tension of liquid drops. We emphasize that any approach, e.g., classical nucleation theory, which is based on a purely macroscopic viewpoint, does not lead to a reliable representation when the radius of the drop becomes microscopic. The description of the enhanced inhomogeneity exhibited by small drops (particularly in the dense interior) necessitates a treatment at the molecular level to account for finite-size and surface effects correctly. The so-called mechanical route, which corresponds to a molecular-level extension of the macroscopic theory of elasticity and is particularly popular in molecular dynamics simulation, also appears to be unreliable due to the inherent ambiguity in the definition of the microscopic pressure tensor, an observation which has been known for decades but is frequently ignored. The union of the theory of capillarity (developed in the nineteenth century by Gibbs and then promoted by Tolman) with a microscopic DFT treatment allows for a direct and unambiguous description of the interfacial properties of drops of arbitrary size; DFT provides all of the bulk and surface characteristics of the system that are required to uniquely define its thermodynamic properties. In this vein, we propose a non-local mean-field DFT for Lennard-Jones (LJ) fluids to examine drops of varying size. A comparison of the predictions of our DFT with recent simulation data based on a second-order fluctuation analysis (Sampayo et al 2010 J. Chem. Phys. 132 141101) reveals the consistency of the two treatments. This observation highlights the significance of fluctuation effects in small drops, which give rise to additional entropic (thermal non-mechanical) contributions, in contrast to what one observes in the case of planar interfaces which are governed by the laws of mechanical equilibrium. A small negative Tolman length (which is found to be about a tenth of the molecular diameter) and a non-monotonic behaviour of the surface tension with the drop radius are predicted for the LJ fluid. Finally, the limits of the validity of the Tolman approach, the effect of the range of the intermolecular potential, and the behaviour of bubbles are briefly discussed.
Collapse
Affiliation(s)
- Alexandr Malijevský
- E Hála Laboratory of Thermodynamics, Institute of Chemical Process Fundamentals of the ASCR, Prague 6, Czech Republic.
| | | |
Collapse
|