1
|
Sanchez-Burgos I, Muniz MC, Espinosa JR, Panagiotopoulos AZ. A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water. J Chem Phys 2023; 158:2889532. [PMID: 37158636 DOI: 10.1063/5.0144500] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Accepted: 03/30/2023] [Indexed: 05/10/2023] Open
Abstract
Computational studies of liquid water and its phase transition into vapor have traditionally been performed using classical water models. Here, we utilize the Deep Potential methodology-a machine learning approach-to study this ubiquitous phase transition, starting from the phase diagram in the liquid-vapor coexistence regime. The machine learning model is trained on ab initio energies and forces based on the SCAN density functional, which has been previously shown to reproduce solid phases and other properties of water. Here, we compute the surface tension, saturation pressure, and enthalpy of vaporization for a range of temperatures spanning from 300 to 600 K and evaluate the Deep Potential model performance against experimental results and the semiempirical TIP4P/2005 classical model. Moreover, by employing the seeding technique, we evaluate the free energy barrier and nucleation rate at negative pressures for the isotherm of 296.4 K. We find that the nucleation rates obtained from the Deep Potential model deviate from those computed for the TIP4P/2005 water model due to an underestimation in the surface tension from the Deep Potential model. From analysis of the seeding simulations, we also evaluate the Tolman length for the Deep Potential water model, which is (0.091 ± 0.008) nm at 296.4 K. Finally, we identify that water molecules display a preferential orientation in the liquid-vapor interface, in which H atoms tend to point toward the vapor phase to maximize the enthalpic gain of interfacial molecules. We find that this behavior is more pronounced for planar interfaces than for the curved interfaces in bubbles. This work represents the first application of Deep Potential models to the study of liquid-vapor coexistence and water cavitation.
Collapse
Affiliation(s)
- Ignacio Sanchez-Burgos
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue,Cambridge CB3 0HE, United Kingdom
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Jorge R Espinosa
- Maxwell Centre, Cavendish Laboratory, Department of Physics, University of Cambridge, J J Thomson Avenue,Cambridge CB3 0HE, United Kingdom
- Departamento de Química Fisica, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | | |
Collapse
|
2
|
Hasan F, Al Mahmud KAH, Khan MI, Kang W, Adnan A. Effect of random fiber networks on bubble growth in gelatin hydrogels. SOFT MATTER 2021; 17:9293-9314. [PMID: 34647568 DOI: 10.1039/d1sm00587a] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
In hydrodynamics, the event of dynamic bubble growth in a pure liquid under tensile pressure is known as cavitation. The same event can also be observed in soft materials (e.g., elastomers and hydrogels). However, for soft materials, bubble/cavity growth is either defined as cavitation if the bubble growth is elastic and reversible or as fracture if the cavity growth is by material failure and irreversible. In any way, bubble growth can cause damage to soft materials (e.g., tissue) by inducing high strain and strain-rate deformation. Additionally, a high-strength pressure wave is generated upon the collapse of the bubble. Therefore, it is crucial to identify the critical condition of spontaneous bubble growth in soft materials. Experimental and theoretical observations have agreed that the onset of bubble growth in soft materials requires higher tensile pressure than pure water. The extra tensile pressure is required since the cavitating bubble needs to overcome the elastic and surface energy in soft materials. In this manuscript, we developed two models to study and quantify the extra tensile pressure for different gelatin concentrations. Both the models are then compared with the existing cavitation onset criteria of rubber-like materials. Validation is done with the experimental results of threshold tensile pressure for different gelatin concentrations. Both models can moderately predict the extra tensile pressure within the intermediate range of gelatin concentrations (3-7% [w/v]). For low concentration (∼1%), the network's non-affinity plays a significant role and must be incorporated. On the other hand, for higher concentrations (∼10%), the entropic deformation dominates, and the strain energy formulation is not adequate.
Collapse
Affiliation(s)
- Fuad Hasan
- Department of Mechanical and Aerospace Engineering, The University of Texas at Arlington, USA.
- Woolf Hall, Room 315C, Arlington, TX 76019, USA
| | - K A H Al Mahmud
- Department of Mechanical and Aerospace Engineering, The University of Texas at Arlington, USA.
- Woolf Hall, Room 315C, Arlington, TX 76019, USA
| | - Md Ishak Khan
- Department of Mechanical and Aerospace Engineering, The University of Texas at Arlington, USA.
- Woolf Hall, Room 315C, Arlington, TX 76019, USA
| | - Wonmo Kang
- School for Engineering of Matter, Transport and Energy, Arizona State University, USA
| | - Ashfaq Adnan
- Department of Mechanical and Aerospace Engineering, The University of Texas at Arlington, USA.
- Woolf Hall, Room 315C, Arlington, TX 76019, USA
| |
Collapse
|
3
|
Berkowitz ML. Stretch-Induced Cavitation: How Critical Cavity Radius and Barrier Energy, Radius, and Energy of a Stable Cavity Depend on the Stretching Factor. J Phys Chem B 2021; 125:4409-4414. [PMID: 33896169 DOI: 10.1021/acs.jpcb.1c00096] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Characteristics of tension-induced cavitation, such as free energy at the barrier for cavitation, the size of the critical (barrier) cavity, the stable cavity size, and the free energy of the stable cavity, depend on the amount of tension (stretch) and the initial size of the sample. In this work, we study how the characteristics of the cavitation mentioned above scale with the amount of applied tension. We consider two models characterizing the properties of cavitating liquid: (a) a simple model with a linear tension-strain relation and neglect of curvature dependence of cavity surface tension and (b) a more realistic model with a nonlinear tension-strain relation and curvature-dependent surface tension. For both models, we find the relevant scaling relations when we stretch the initial volume of the liquid sample in the interval between 1% and 20% of the initial volume. Specific numerical tests are performed for the case of liquid water when the initial volume of the sample is a sphere with a radius of 100 nm.
Collapse
Affiliation(s)
- Max L Berkowitz
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
4
|
Hewage SA, Kewalramani J, Meegoda JN. Stability of nanobubbles in different salts solutions. Colloids Surf A Physicochem Eng Asp 2021. [DOI: 10.1016/j.colsurfa.2020.125669] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
|
5
|
Montero de Hijes P, Shi K, Noya EG, Santiso EE, Gubbins KE, Sanz E, Vega C. The Young–Laplace equation for a solid–liquid interface. J Chem Phys 2020; 153:191102. [DOI: 10.1063/5.0032602] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- P. Montero de Hijes
- Faculty of Chemistry, Chemical Physics Department, Universidad Complutense de Madrid, Plaza de las Ciencias, Ciudad Universitaria, Madrid 28040, Spain
| | - K. Shi
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, USA
| | - E. G. Noya
- Instituto de Química Física Rocasolano, Consejo Superior de Investigaciones Científicas, CSIC, Calle Serrano 119, 28006 Madrid, Spain
| | - E. E. Santiso
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, USA
| | - K. E. Gubbins
- Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, North Carolina 27606, USA
| | - E. Sanz
- Faculty of Chemistry, Chemical Physics Department, Universidad Complutense de Madrid, Plaza de las Ciencias, Ciudad Universitaria, Madrid 28040, Spain
| | - C. Vega
- Faculty of Chemistry, Chemical Physics Department, Universidad Complutense de Madrid, Plaza de las Ciencias, Ciudad Universitaria, Madrid 28040, Spain
| |
Collapse
|
6
|
Manning GS. On the thermodynamic stability of bubbles, immiscible droplets, and cavities. Phys Chem Chem Phys 2020; 22:17523-17531. [PMID: 32728677 DOI: 10.1039/d0cp02517h] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Nanobubbles filled with air or a variety of pure gases are observed to persist in bulk water for weeks and months. Nanoemulsions consisting of oil droplets in water are also remarkably stable against coagulation, with lifetimes up to weeks even if not coated with surfactants. The inverse system of nanodroplets of water in oil is also accessible for study and application. Voids on the nanoscale are formed in simulations of water under strong tension and are stable during the time of the simulation. The stability of these nano-entities is ultimately determined by the molecular-level structure of their interfaces. However, a thermodynamic theory might also be capable of providing some insight. We therefore consider spherical gas nanobubbles, immiscible liquid nanodroplets, and nanocavities formed in water under negative pressure on the same footing, and give a unified thermodynamic analysis of these systems. In all cases, mechanical equilibrium (local free energy maximum or minimum) is expressed by the Laplace equation, and thermodynamic stability (local free energy minimum) follows from the radius dependence of surface tension. All of them would be unstable if their surface tensions were constant. Data from the literature allow construction of numerical examples for cavities and gas nanobubbles. Spectroscopic data are cited in support of an interfacial water structure in gas nanobubbles and water droplets in oil that differ from their flat surface counterparts. The observed longevity of nanobubbles in particular has been thought to violate fundamental principles of diffusion and solubility. A close look at the Laplace equation and its derivation shows why this widespread belief is incorrect.
Collapse
Affiliation(s)
- Gerald S Manning
- Department of Chemistry and Chemical Biology, Rutgers University, 123 Bevier Road, Piscataway, NJ 08854-8087, USA.
| |
Collapse
|
7
|
Min SH, Wijesinghe S, Berkowitz ML. Enhanced Cavitation and Hydration Crossover of Stretched Water in the Presence of C 60. J Phys Chem Lett 2019; 10:6621-6625. [PMID: 31609628 DOI: 10.1021/acs.jpclett.9b02511] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
We performed molecular dynamics simulations on systems containing stretched water and a C60 buckyball molecule. Our goals were to understand how the presence of the hydrophobic impurity influences the rate of cavitation in stretched water and how the change in pressure (an increase in the value of negative pressure) affects the nature of hydrophobic hydration. Our simulations show that the presence of a buckyball increases the rate of cavitation in water under negative pressure. When studying the influence of the degree of stretching on hydration, we observed that at pressures above -100 MPa the mechanism of hydrophobic hydration is the one that characterizes hydration of a small particle. At some pressure below -100 MPa, there is a crossover in the mechanism of hydration where dewetting occurs by forming cavities next to the surface of the buckyball, and this is characteristic of hydrophobic hydration of large particles.
Collapse
Affiliation(s)
- Sa Hoon Min
- Department of Chemistry , The University of North Carolina at Chapel Hill , Chapel Hill , North Carolina 27599 , United States
| | - Sidath Wijesinghe
- Department of Chemistry , The University of North Carolina at Chapel Hill , Chapel Hill , North Carolina 27599 , United States
| | - Max L Berkowitz
- Department of Chemistry , The University of North Carolina at Chapel Hill , Chapel Hill , North Carolina 27599 , United States
| |
Collapse
|
8
|
Malek SMA, Poole PH, Saika-Voivod I. Surface tension of supercooled water nanodroplets from computer simulations. J Chem Phys 2019; 150:234507. [DOI: 10.1063/1.5096990] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Affiliation(s)
- Shahrazad M. A. Malek
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John’s, Newfoundland A1B 3X7, Canada
| | - Peter H. Poole
- Department of Physics, St. Francis Xavier University, Antigonish, Nova Scotia B2G 2W5, Canada
| | - Ivan Saika-Voivod
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John’s, Newfoundland A1B 3X7, Canada
| |
Collapse
|