Homogeneous bubble nucleation in water at negative pressure: A Voronoi polyhedra analysis

Interfacial Thermal Fluctuations Stabilize Bulk Nanobubbles

Consensus on bulk nanobubble stability remains elusive, despite accepted indirect evidence for longevity. We develop a nanobubble evolution model by incorporating thermal capillary wave theory that reveals that dense nanobubbles generated by acoustic cavitation tend to shrink and intensify interfacial thermal fluctuations; this significantly reduces surface tension to neutralize enhanced Laplace pressure, and secures their stabilization at a finite size. A stability criterion emerges: thermal fluctuation intensity scales superlinearly with curvature: sqrt[⟨h^{2}⟩]∝(1/R)^{n}, n>1. The model prolongs the time frame for nanobubble contraction to 2 orders of magnitude beyond classical theory estimates, and captures the equilibrium radius (90-215 nm) within the experimental range.

Stochastic computer experiments of the thermodynamic irreversibility of bulk nanobubbles in supersaturated and weak gas-liquid solutions

Spontaneous gas-bubble nucleation in weak gas-liquid solutions has been a challenging topic in theory, experimentation, and computer simulations. In analogy with recent advances in crystallization and droplet formation studies, the diffusive-shielding stabilization and thermodynamic irreversibility of bulk nanobubble (bNB) mechanisms are revisited and deployed to characterize nucleation processes in a stochastic framework of computer experiments using the large-scale atomic/molecular massively parallel simulator code. Theoretical bases, assumptions, and limitations underlying the irreversibility hypothesis of bNBs, and their computational counterparts, are extensively described and illustrated. In essence, it is established that the irreversibility hypothesis can be numerically investigated by converging the system volume (due to the finiteness of interatomic forces) and the initial dissolved-gas concentration in the solution (due to the single-bNB limitation). Helium nucleation in liquid Pb17Li alloy is selected as a representative case study, where it exhibits typical characteristics of noble-gas/liquid-metal systems. The proposed framework lays down the bases on which the stability of gas-bNBs in weak and supersaturated gas-liquid solutions can be inferred and explained from a novel perspective. In essence, it stochastically marches toward a unique irreversible state along out-of-equilibrium nucleation/growth trajectories. Moreover, it does not attempt to characterize the interface or any interface-related properties, neither theoretically nor computationally. It was concluded that bNBs of a few tens of He-atoms are irreversible when dissolved-He concentrations in the weak gas-liquid solution are at least ∼50 and ∼105 mol m-3 at 600 and 1000 K (and ∼80 MPa), respectively, whereas classical molecular dynamics -estimated solubilities are at least two orders of magnitude smaller.

A Deep Potential model for liquid-vapor equilibrium and cavitation rates of water

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.

Anomalous Vapor and Ice Nucleation in Water at Negative Pressures: A Classical Density Functional Theory Study

In contrast to the abundance of work on the anomalous behavior of water, the relationship between the water's thermodynamic anomalies and kinetics of phase transition from metastable water is relatively unexplored. In this work, we have employed classical density functional theory to provide a unified and coherent picture of nucleation (both vapor and ice) from metastable water at negative pressure conditions. Our results suggest a peculiar nonmonotonic temperature dependence of vapor-liquid surface tension at temperatures where vapor-liquid coexistence is metastable with respect to the ice phase. The vapor nucleation barrier on isochoric cooling also shows a nonmonotonic temperature dependence. We further report that, for low density isochores, the temperature of the minimum vapor nucleation barrier (TΔΩv/min*) does not coincide with the temperature of maximum density (TMD) where metastability is maximum. The difference between the TΔΩv/min* and the TMD, however, decreases with increasing the density of the isochore. The vapor nucleation barrier along isobars shows an interesting crossover behavior in the vicinity of the Widom line on lowering the temperature. Our results on the ice nucleation suggest an anomalous retracing behavior of the nucleation barrier along isotherms at negative pressures and theoretically validate the recent findings that the reentrant ice(Ih)-liquid coexistence line can induce a drastic change in the kinetics of ice nucleation. Thus, this study establishes a direct connection between the metastable water's thermodynamic anomalies and the (vapor and ice) nucleation kinetics. In addition, this study provides deeper insights into the origin of the isothermal compressibility maximum on isochoric cooling.

Dynamic Surface Tension Enhances the Stability of Nanobubbles in Xylem Sap

Air seeded nanobubbles have recently been observed within tree sap under negative pressure. They are stabilized by an as yet unidentified process, although some embolize their vessels in extreme circumstances. Current literature suggests that a varying surface tension helps bubbles survive, but few direct measurements of this quantity have been made. Here, we present calculations of dynamic surface tension for two biologically relevant lipids using molecular dynamics simulations. We find that glycolipid monolayers resist expansion proportionally to the rate of expansion. Their surface tension increases with the tension applied, in a similar way to the viscosity of a non-Newtonian fluid. In contrast, a prototypical phospholipid was equally resistant to all applied tensions, suggesting that the fate of a given nanobubble is dependent on its surface composition. By incorporating our results into a Classical Nucleation Theory (CNT) framework, we predict nanobubble stability with respect to embolism. We find that the metastable radius of glycolipid coated nanobubbles is approximately 35 nm, and that embolism is in this case unlikely when the external pressure is less negative than -1.5 MPa.

Cavitation Mean Expectation Time in a Stretched Lennard-Jones Fluid under Confinement

We investigate the nucleation of cavitation bubbles in a confined Lennard-Jones fluid subjected to negative pressures in a cubic enclosure. We perform molecular dynamics (MD) simulations with tunable interatomic potentials that enable us to control the wettability of solid walls by the liquid, that is, its contact angle. For a given temperature and pressure, as the solid is taken more hydrophobic, we put in evidence, an increase in nucleation probability. A Voronoi tessellation method is used to accurately detect the bubble appearance and its nucleation rate as a function of the contact angle. We adapt classical nucleation theory (CNT) proposed for the heterogeneous case on a flat surface to our situation where bubbles may appear on flat walls, edges, or corners of the confined box. We finally calculate a theoretical mean expectation time in these three cases. The ratio of these calculated values over the homogeneous case is computed and compared successfully against MD simulations. Beyond the infinite liquid case, this work explores the heterogeneous nucleation of cavitation bubbles, not only in the flat surface case but for more complex confining geometries.

fromage: A library for the study of molecular crystal excited states at the aggregate scale

The study of photoexcitations in molecular aggregates faces the twofold problem of the increased computational cost associated with excited states and the complexity of the interactions among the constituent monomers. A mechanistic investigation of these processes requires the analysis of the intermolecular interactions, the effect of the environment, and 3D arrangements or crystal packing on the excited states. A considerable number of techniques have been tailored to navigate these obstacles; however, they are usually restricted to in‐house codes and thus require a disproportionate effort to adopt by researchers approaching the field. Herein, we present the FRamewOrk for Molecular AGgregate Excitations (fromage), which implements a collection of such techniques in a Python library complemented with ready‐to‐use scripts. The program structure is presented and the principal features available to the user are described: geometrical analysis, exciton characterization, and a variety of ONIOM schemes. Each is illustrated by examples of diverse organic molecules in condensed phase settings. The program is available at https://github.com/Crespo-Otero-group/fromage.

Cavitation inception of water with solid nanoparticles: A molecular dynamics study

Cavitation in liquid with impurities is important in heterogeneous nucleation applications. One of the most widely existing kinds of impurities is solid particles, which can be found in natural water from rivers and specially prepared water such as nanofluids. Understanding the effects caused by the existence of nanoparticles on cavitation in water is vital to the rapidly developed nanotechnologies and medical researches. In this study, cavitation in water with nanoparticles is investigated through molecular dynamics simulations. The effects by nanoparticle materials and sizes on cavitation are discussed by using SiO₂ and polyethylene spherical nanoparticles with different diameters. The nucleation rate and the formation of critical bubbles in cavitation are studied via the Voronoi tessellation and the mean first passage time methods. The hydrogen bond network in water is also analyzed. Results reveal that SiO⁠₂ and polyethylene nanoparticles may destabilize the hydrogen bond network in water. With the same particle size, cavitation in water with polyethylene nanoparticles is promoted to a greater extent than that with SiO₂ nanoparticles. With the same nanoparticle material, cavitation is promoted with the increase in particle size in a range spanning half to ten times the critical bubble radius. Beyond this range, particle size has little influence on cavitation. Reasons for those effects on cavitation due to the presence of solid nanoparticles are discussed by analysing the changes of hydrogen bonds network in water.

Molecular mechanism for cavitation in water under tension

Despite its relevance in biology and engineering, the molecular mechanism driving cavitation in water remains unknown. Using computer simulations, we investigate the structure and dynamics of vapor bubbles emerging from metastable water at negative pressures. We find that in the early stages of cavitation, bubbles are irregularly shaped and become more spherical as they grow. Nevertheless, the free energy of bubble formation can be perfectly reproduced in the framework of classical nucleation theory (CNT) if the curvature dependence of the surface tension is taken into account. Comparison of the observed bubble dynamics to the predictions of the macroscopic Rayleigh-Plesset (RP) equation, augmented with thermal fluctuations, demonstrates that the growth of nanoscale bubbles is governed by viscous forces. Combining the dynamical prefactor determined from the RP equation with CNT based on the Kramers formalism yields an analytical expression for the cavitation rate that reproduces the simulation results very well over a wide range of pressures. Furthermore, our theoretical predictions are in excellent agreement with cavitation rates obtained from inclusion experiments. This suggests that homogeneous nucleation is observed in inclusions, whereas only heterogeneous nucleation on impurities or defects occurs in other experiments.

Reply to "Comment on 'Simple improvements to classical bubble nucleation models' "

We reply to the Comment by Schmelzer and Baidakov [Phys. Rev. E 94, 026801 (2016)].10.1103/PhysRevE.94.026801 They suggest that a more modern approach than the classic description by Tolman is necessary to model the surface tension of curved interfaces. Therefore we now consider the higher-order Helfrich correction, rather than the simpler first-order Tolman correction. Using a recent parametrization of the Helfrich correction provided by Wilhelmsen et al. [J. Chem. Phys. 142, 064706 (2015)]JCPSA60021-960610.1063/1.4907588, we test this description against measurements from our simulations, and find an agreement stronger than what the pure Tolman description offers. Our analyses suggest a necessary correction of order higher than the second for small bubbles with radius ≲1 nm. In addition, we respond to other minor criticism about our results.

Effects of ensembles on methane hydrate nucleation kinetics

By performing molecular dynamics simulations to form a hydrate with a methane nano-bubble in liquid water at 250 K and 50 MPa, we report how different ensembles, such as the NPT, NVT, and NVE ensembles, affect the nucleation kinetics of the methane hydrate. The nucleation trajectories are monitored using the face-saturated incomplete cage analysis (FSICA) and the mutually coordinated guest (MCG) order parameter (OP). The nucleation rate and the critical nucleus are obtained using the mean first-passage time (MFPT) method based on the FS cages and the MCG-1 OPs, respectively. The fitting results of MFPT show that hydrate nucleation and growth are coupled together, consistent with the cage adsorption hypothesis which emphasizes that the cage adsorption of methane is a mechanism for both hydrate nucleation and growth. For the three different ensembles, the hydrate nucleation rate is quantitatively ordered as follows: NPT > NVT > NVE, while the sequence of hydrate crystallinity is exactly reversed. However, the largest size of the critical nucleus appears in the NVT ensemble, rather than in the NVE ensemble. These results are helpful for choosing a suitable ensemble when to study hydrate formation via computer simulations, and emphasize the importance of the order degree of the critical nucleus.

Simple improvements to classical bubble nucleation models

We revisit classical nucleation theory (CNT) for the homogeneous bubble nucleation rate and improve the classical formula using a correct prefactor in the nucleation rate. Most of the previous theoretical studies have used the constant prefactor determined by the bubble growth due to the evaporation process from the bubble surface. However, the growth of bubbles is also regulated by the thermal conduction, the viscosity, and the inertia of liquid motion. These effects can decrease the prefactor significantly, especially when the liquid pressure is much smaller than the equilibrium one. The deviation in the nucleation rate between the improved formula and the CNT can be as large as several orders of magnitude. Our improved, accurate prefactor and recent advances in molecular dynamics simulations and laboratory experiments for argon bubble nucleation enable us to precisely constrain the free energy barrier for bubble nucleation. Assuming the correction to the CNT free energy is of the functional form suggested by Tolman, the precise evaluations of the free energy barriers suggest the Tolman length is ≃0.3σ independently of the temperature for argon bubble nucleation, where σ is the unit length of the Lennard-Jones potential. With this Tolman correction and our prefactor one gets accurate bubble nucleation rate predictions in the parameter range probed by current experiments and molecular dynamics simulations.

Bubble nucleation in simple and molecular liquids via the largest spherical cavity method

In this work, we propose a methodology to compute bubble nucleation free energy barriers using trajectories generated via molecular dynamics simulations. We follow the bubble nucleation process by means of a local order parameter, defined by the volume of the largest spherical cavity (LSC) formed in the nucleating trajectories. This order parameter simplifies considerably the monitoring of the nucleation events, as compared with the previous approaches which require ad hoc criteria to classify the atoms and molecules as liquid or vapor. The combination of the LSC and the mean first passage time technique can then be used to obtain the free energy curves. Upon computation of the cavity distribution function the nucleation rate and free-energy barrier can then be computed. We test our method against recent computations of bubble nucleation in simple liquids and water at negative pressures. We obtain free-energy barriers in good agreement with the previous works. The LSC method provides a versatile and computationally efficient route to estimate the volume of critical bubbles the nucleation rate and to compute bubble nucleation free-energies in both simple and molecular liquids.

Evaluation of finite-size effects in cavitation and droplet formation

Nucleation of bubbles and droplets is of fundamental interest in science and technology and has been widely investigated through experiments, theory, and simulations. Giving the rare event nature of these phenomena, nucleation simulations are computationally costly and require the use of a limited number of particles. Moreover, they are often performed in the canonical ensemble, i.e., by fixing the total volume and number of particles, to avoid the additional complexities of implementing a barostat. However, cavitation and droplet formation take place differently depending on the ensemble. Here, we analyze the importance of finite-size effects in cavitation and droplet formation. We present simple formulas which predict the finite-size corrections to the critical size, the nucleation barrier, and the nucleation rates in the canonical ensemble very accurately. These results can be used to select an appropriate system-size for simulations and to get a more precise evaluation of nucleation in complex substances, by using a small number of molecules and correcting for finite-size effects.

Bubble evolution and properties in homogeneous nucleation simulations

We analyze the properties of naturally formed nanobubbles in Lennard-Jones molecular dynamics simulations of liquid-to-vapor nucleation in the boiling and the cavitation regimes. The large computational volumes provide a realistic environment at unchanging average temperature and liquid pressure, which allows us to accurately measure properties of bubbles from their inception as stable, critically sized bubbles, to their continued growth into the constant speed regime. Bubble gas densities are up to 50% lower than the equilibrium vapor densities at the liquid temperature, yet quite close to the gas equilibrium density at the lower gas temperatures measured in the simulations: The latent heat of transformation results in bubble gas temperatures up to 25% below those of the surrounding bulk liquid. In the case of rapid bubble growth-typical for the cavitation regime-compression of the liquid outside the bubble leads to local temperature increases of up to 5%, likely significant enough to alter the surface tension as well as the local viscosity. The liquid-vapor bubble interface is thinner than expected from planar coexistence simulations by up to 50%. Bubbles near the critical size are extremely nonspherical, yet they quickly become spherical as they grow. The Rayleigh-Plesset description of bubble-growth gives good agreement in the cavitation regime.

Direct simulations of homogeneous bubble nucleation: Agreement with classical nucleation theory and no local hot spots

We present results from direct, large-scale molecular dynamics simulations of homogeneous bubble (liquid-to-vapor) nucleation. The simulations contain half a billion Lennard-Jones atoms and cover up to 56 million time steps. The unprecedented size of the simulated volumes allows us to resolve the nucleation and growth of many bubbles per run in simple direct micro-canonical simulations while the ambient pressure and temperature remain almost perfectly constant. We find bubble nucleation rates which are lower than in most of the previous, smaller simulations. It is widely believed that classical nucleation theory (CNT) generally underestimates bubble nucleation rates by very large factors. However, our measured rates are within two orders of magnitude of CNT predictions; only at very low temperatures does CNT underestimate the nucleation rate significantly. Introducing a small, positive Tolman length leads to very good agreement at all temperatures, as found in our recent vapor-to-liquid nucleation simulations. The critical bubbles sizes derived with the nucleation theorem agree well with the CNT predictions at all temperatures. Local hot spots reported in the literature are not seen: Regions where a bubble nucleation event will occur are not above the average temperature, and no correlation of temperature fluctuations with subsequent bubble formation is seen.

Size-Dependent Surface Free Energy and Tolman-Corrected Droplet Nucleation of TIP4P/2005 Water

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.

Instantaneous normal mode analysis for intermolecular and intramolecular vibrations of water from atomic point of view

By exploiting the instantaneous normal mode (INM) analysis for models of flexible molecules, we investigate intermolecular and intramolecular vibrations of water from the atomic point of view. With two flexible SPC/E models, our investigations include three aspects about their INM spectra, which are separated into the unstable, intermolecular, bending, and stretching bands. First, the O- and H-atom contributions in the four INM bands are calculated and their stable INM spectra are compared with the power spectra of the atomic velocity autocorrelation functions. The unstable and intermolecular bands of the flexible models are also compared with those of the SPC/E model of rigid molecules. Second, we formulate the inverse participation ratio (IPR) of the INMs, respectively, for the O- and H-atom and molecule. With the IPRs, the numbers of the three species participated in the INMs are estimated so that the localization characters of the INMs in each band are studied. Further, by the ratio of the IPR of the H atom to that of the O atom, we explore the number of involved OH bond per molecule participated in the INMs. Third, by classifying simulated molecules into subensembles according to the geometry of their local environments or their H-bond configurations, we examine the local-structure effects on the bending and stretching INM bands. All of our results are verified to be insensible to the definition of H-bond. Our conclusions about the intermolecular and intramolecular vibrations in water are given.

Toward a molecular theory of homogeneous bubble nucleation: II. Calculation of the number density of critical nuclei and the rate of nucleation

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.