Molecular dynamics simulation of nanobubble nucleation on rough surfaces

Molecular Mechanisms of the Generation and Accumulation of Gas at the Interface

Gas-evolving reactions are widespread in chemical and energy fields. However, the generated gas will accumulate at the interface, which reduces the rate of gas generation. Understanding the microscopic processes of the generation and accumulation of gas at the interface is crucial for improving the efficiency of gas generation. Here, we develop an algorithm to reproduce the process of catalytic gas generation at the molecular scale based on the all-atom molecular dynamics simulations and obtain the quantitative evolution of the gas generation, which agrees well with the experimental results. In addition, we demonstrate that under an external electric field, the generated gas molecules do not accumulate at the electrode surface, which implies that the electric field can significantly increase the rate of the gas generation. The results suggest that the external electric field changes the structure of the water molecules near the electrode surface, making it difficult for gas molecules to accumulate on the electrode surface. Furthermore, it is found that gas desorption from the electrode surface is an entropy-driven process, and its accumulation at the electrode surface depends mainly on the competition between the entropy and the enthalpy of the water molecules under the influence of the electric field. These results provide deep insight into gas generation and inhibition of gas accumulation.

Surface Morphology Enriching the Energy Barrier Leads to the Adsorption Characteristic of Nanobubbles

With the emergence of nanobubble research, nanobubble distribution morphology at the interface and its stability control become the bottlenecks of nanobubble resistance reduction applications. In this paper, the evolutionary behavior of nanobubbles on smooth and step HOPG surfaces was compared through molecular dynamics studies. The results show that the surface energy barrier provided by the step HOPG surface restricts diffusion of gas molecules. Then, a method of multisolvent evaporation for preparing hydrophobic nanoindent surfaces was proposed, which can achieve phase separation through different evaporation rates of multisolvents, thus realizing the preparation of surface structures with uniform distribution of nanoindents. In this paper, the nucleation processes of nanobubbles on PS nanoindent hydrophobic surface, HOPG flat hydrophobic surface, and HOPG nanostep hydrophobic surface were compared by using atomic force microscopy in liquid experiment. The evolution of the volume and distribution morphology of nanobubbles on the three nanostructures was observed by 24 h in situ tests, revealing that the energy barrier effect arising from the uneven surface structure can effectively prevent adjacent nanobubbles from merging in close proximity to each other. It is also pointed out that the hydrophobic nanoindents prepared by using the multivariate solvent evaporation method in this paper can cover most of nanobubbles for stable adsorption. It can be seen from the results that the volume drop of the nanobubbles on the HOPG flat hydrophobic surface is 27% and that on the HOPG nanostep and the PS nanoindent hydrophobic surface it is reduced to 19% and 3% under the effect of structural energy barriers, respectively. The density of the nanostructures determines whether the existence of nanobubbles is stable. The coverage of nanobubbles on the HOPG flat hydrophobic surface was 3.313% when the existence of nanobubbles was mostly stable. The HOPG nanostep and PS nanoindent sizes were positively correlated with the morphological size of the nanobubbles, which increased the coverage of the nanobubbles on the hydrophobic surface of the HOPG nanostep and PS nanoindent to 5.229% and 4.437%, respectively, when the existence of nanobubbles was mostly stable.

Second-Level Microgroove Convexity is Critical for Air Plastron Restoration on Immersed Hierarchical Superhydrophobic Surfaces

Air plastrons trapped on the surfaces of underwater superhydrophobic surfaces are critical for their function. Fibrillar morphologies offer a natural pathway, yet they are limited to a narrow range of liquid-surface systems and are vulnerable to pressure fluctuations that irreversibly destroy the air layer plastron. Inspired by the convexly grooved bases of water fern (Salvinia) leaves that support their fibrous outgrowths, we focus on the effect of such second-level grooved structures or microgrooves on the plastron restoration on immersed three-dimensional (3D)-printed hierarchical surfaces. Elliptical, interconnected microgrooves are fabricated with varying surface curvatures to study the effect of their morphology. Immersion experiments reveal that the convex groove curvature stabilizes a seed gas layer (SGL) that facilitates plastron restoration for all immersed hydrophobic surfaces. Theoretical calculations and atomic-scale computations reveal that the SGL storage capacity that sets the SGL robustness follows from the liquid menisci adaption to the groove geometry and pressure, from micro- to nanoscales, and it can be further tuned using separated grooves. Our study highlights groove convexity as a key morphological feature for the design of second-level architectures for underwater air plastron restoration on hierarchical superhydrophobic surfaces.

Cassie State Stability and Gas Restoration Capability of Superhydrophobic Surfaces with Truncated Cone-Shaped Pillars

The gas layer stability on superhydrophobic surfaces and gas restoration on the immersed superhydrophobic surfaces have been great challenges for their practical applications in recent years. Inspired by the naturally existing mushroom-like super-repellent superhydrophobic patterns, we choose superhydrophobic surfaces with truncated cone-shaped pillars as our research objects to tackle such challenges by tuning their geometrical parameters. We perform molecular dynamics simulations to investigate the Cassie-Wenzel transition under external pressure and the Wenzel-Cassie transition due to underwater spreading of compressed bubbles. Theories based on the Young-Laplace equation and total free-energy variation are developed to explore the influence of geometrical parameters of pillars on the pressure resistance and underwater gas restoration, which is in good agreement with simulation results. These simulation results and theoretical analysis suggest that cork-shaped pillars, analogous to the surface structures of natural organisms like springtails and Salvinia leaves, can be super-repellent to the liquid and favorable for the gas spreading process. Our study provides theoretical guidance for the design of superhydrophobic surfaces with both Cassie state stability and gas restoration capability.

Molecular Dynamics Study of Bubble Nucleation on an Ideally Smooth Substrate

Questions regarding bubble nucleation on an ideally smooth surface are seemingly endless, but it can not be adequately verified yet because of the scale limitation (microscopic scale). Hence, in this study, bubble nucleation on an ideally smooth substrate is explored using the molecular dynamics simulation method. An ideally smooth hydrophilic platinum substrate at 145 K is conducted to heat the simple L-J liquid argon. Results show that a visible bubble nucleus successfully forms on the ideally smooth substrate without any additional disturbance, which is common in boiling studies using the traditional numerical simulation methods. However, the nucleation position is unpredictable. At the atomic level, the thermal energy transfer from an ideally smooth substrate to liquid atoms is inhomogeneous due to atomic inhomogeneous distribution and irregular movement, which are the key influencing factors for achieving bubble nucleation. The inhomogeneity will be highlighted with the heating process. As a result, some local liquid atoms near the ideally smooth surface absorb more thermal energy to overcome their potential barrier at a specific moment, causing the emergence of a distinct nucleus there. Furthermore, nanostructure substrates are introduced to make a comparison with the smooth substrate in bubble nucleation. There is no significant difference in the inception temperature of nucleation between the ideally smooth and nanostructure substrates, but the latter has better performance in improving the bubble nucleation rate.

Stability of pinned surface nanobubbles against expansion: Insights from theory and simulation

While growth and dissolution of surface nanobubbles have been widely studied in recent years, their stability under pressure changes or a temperature increase has not received the same level of scrutiny. Here, we present theoretical predictions based on classical theory for pressure and temperature thresholds (p_c and T_c) at which unstable growth occurs for the case of air nanobubbles on a solid surface in water. We show that bubbles subjected to pinning have much lower p_c and higher T_c compared to both unpinned and bulk bubbles of similar size, indicating that pinned bubbles can withstand a larger tensile stress (negative pressure) and higher temperatures. The values of p_c and T_c obtained from many-body dissipative particle dynamics simulations of quasi-two-dimensional (quasi-2D) surface nanobubbles are consistent with the theoretical predictions, provided that the lateral expansion during growth is taken into account. This suggests that the modified classical thermodynamic description is valid for pinned bubbles as small as several nanometers. While some discrepancies still exist between our theoretical results and previous experiments, further experimental data are needed before a comprehensive understanding of the stability of surface nanobubbles can be achieved.

Electrochemical Generation of Individual Nanobubbles Comprising H2, D2, and HD

The electrochemical reduction of deuterons (2D⁺ + 2e^- → D₂) at Pt nanodisk electrodes (radius = 15-100 nm) in D₂O solutions containing deuterium chloride (DCl) results in the formation of a single gas nanobubble at the electrode surface. Analogous to that previously observed for the electrochemical generation of H₂ nanobubbles, the nucleation and growth of a stable D₂ nanobubble is characterized in voltammetric experiments by a highly reproducible and well-resolved sudden drop in the faradaic current, a consequence of restricted mass transport of D⁺ to the electrode surface following the liquid-to-gas phase transition. D₂ nanobubbles are stable under potential control due to a dynamic equilibrium existing between D₂ gas dissolution and electrochemical generation of D₂ at the circumference of the Pt nanodisk electrode. Remarkably, within the error of the experimental measurement (<6%), the electrochemical current required to nucleate a D₂ gas phase in a D₂O solution is identical to that for the H₂ gas phase in a H₂O solutions, indicating that the concentration required for nucleating a D₂ nanobubble in D₂O (0.29 M) is ∼1.25 times larger than that for a H₂ nanobubble (0.23 M), while the supersaturation is ∼300 in each case. We further demonstrate that individual nanobubbles can be electrogenerated in mixed D₂O/H₂O solutions containing both D⁺ and H⁺ at respective individual concentrations well below those required to nucleate a gas phase containing either pure D₂ or H₂. This latter finding indicates that the resulting nanobubbles comprise a mixture of D₂, H₂, and HD molecules with the chemical composition of a nanobubble determined by the concentrations and diffusivities of D⁺ and H⁺ in the mixed D₂O/H₂O solutions.

Nucleation and Growth of a Nanobubble on Rough Surfaces

We study the nucleation and growth of a nanobubble on rough surfaces using molecular dynamics simulations. A nanobubble nucleates and grows by virtue of a heterogeneous surface reaction which results in the production of gas molecules near the surface. We study the role of surface roughness in the nucleation and growth behavior of a nanobubble. We perform simulations at various reaction rates and surface morphology and quantified the growth dynamics of a nanobubble. Our simulations show that after the onset of nucleation, the nanobubble grows rapidly with radius following t^1/3 behavior followed by a diffusive growth regime which is marked by t^1/2 growth behavior. This growth behavior remains independent of surface roughness and reaction rates over the range considered in this study. We also quantified the oversaturation of gas required for nucleation of a nanobubble and demonstrated its dependence on the surface morphology.

Measuring temperature effects on nanobubble nucleation via a solid-state nanopore

In this study, we designed SiN_X solid-state nanopores to detect the temperature effect on the hydrogen nanobubble formation. Here, we integrated a temperature controller with the highly sensitive nanopore. As the temperature decreases from 25 °C to 5 °C, the occurrence of the nanobubble nucleation inside a 12.3 nm SiN_X nanopore confined space decreased from 102 s^-1 to 23 s^-1, and the life-time of nanobubbles increased from 1.16 ms to 4.78 ms. The results further gave the activation energy for nanobubble nucleation which was 8.1 × 10^-20 J with a 12.3 nm SiN_X nanopore. Our method provides an efficient analytical tool for revealing the temperature-dependent nanobubble nucleation, which further benefits the fundamental understanding of nanobubble nucleation.

Enhanced Cavitation and Hydration Crossover of Stretched Water in the Presence of C60

We performed molecular dynamics simulations on systems containing stretched water and a C₆₀ 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.

Monitoring nanobubble nucleation at early-stage within a sub-9 nm solid-state nanopore

Nanobubble nucleation study is important for understanding the dynamic behavior of nanobubble growth, which is instructive for the nanobubble applications. Benefiting from nanopore fabrication, herein, we fabricated a sub-9 nm SiN_X nanopore with the comparable size to nanobubbles at early-stage. The confined nanopore interface serves as a generator for producing nanobubbles by the chemical reaction between NaBH₄ and H₂ O and as an ultra-sensitive sensor for monitoring the H₂ nanobubble nucleation process. By carrying out the NaBH₄ concentration-dependent experiments, we found the life-time of nanobubbles decreased 250 times and the frequency of nanobubble generation increased 38 times with the NaBH₄ concentration increasing from 6 to 100 mM. The long-time equilibrium between gas molecules inward flux and outward flux could prolong the life-time of nanobubbles to hundreds of milliseconds at low NaBH₄ concentration. The raw current trace depicted that the transient accumulation and dissolution of cavity occurred during all the life-time of nanobubbles. Therefore, the sub-9 nm SiN_X nanopore shows a strong ability for real-time monitoring the nanobubble nucleation at early-stage with high temporal and spatial resolution. This work provides a guide to study the dynamic and stochastic characteristics of nanobubbles.

Nanoscale Study of Bubble Nucleation on a Cavity Substrate Using Molecular Dynamics Simulation

In this paper, the molecular dynamics simulation method is utilized to investigate the phase transition behavior of an argon film placed on cavity substrates with different wettability conditions. A simple Lennard-Jones liquid is heated by a metal platinum substrate at different temperatures, and a complete process of bubble nucleation is successfully visualized on the cavity substrate at temperatures of 150 and 160 K. Moreover, the bubble nucleation behavior shows dependence on cavity wettability. A layer of liquid atom is attracted to the strongly hydrophilic cavity and obtains more energy to nucleate first. In contrast, the liquid atom suffers a large repulsive force from the metal atom in the hydrophobic cavity, thus an original small bubble nucleus stably stays inside before the incipient boiling time. With an increase in the heating time, the original bubble nucleus grows up from the hydrophobic cavity. This bubble nucleation behavior on a hydrophobic cavity is in agreement with macro theory, which states that a cavity provides an original nucleus for bubble formation and growth. Besides, cavity wettability plays a crucial role in the incipient boiling temperature of an argon film. The incipient boiling temperature increases with the weakening of the cavity hydrophobicity, and this trend is in accordance with macro experiments, which show that liquid is easier to boil on a more hydrophobic substrate.

Leakiness of Pinned Neighboring Surface Nanobubbles Induced by Strong Gas-Surface Interaction

The stability of two neighboring surface nanobubbles on a chemically heterogeneous surface is studied by molecular dynamics (MD) simulations of binary mixtures consisting of Lennard-Jones (LJ) particles. A diffusion equation-based stability analysis suggests that two nanobubbles sitting next to each other remain stable, provided the contact line is pinned, and that their radii of curvature are equal. However, many experimental observations seem to suggest some long-term kind of ripening or shrinking of the surface nanobubbles. In our MD simulations we find that the growth/dissolution of the nanobubbles can occur due to the transfer of gas particles from one nanobubble to another along the solid substrate. That is, if the interaction between the gas and the solid is strong enough, the solid-liquid interface can allow for the existence of a "tunnel" which connects the liquid-gas interfaces of the two nanobubbles to destabilize the system. The crucial role of the gas-solid interaction energy is a nanoscopic element that hitherto has not been considered in any macroscopic theory of surface nanobubbles and may help to explain experimental observations of the long-term ripening.

Theory of nanobubble formation and induced force in nanochannels

This paper presents a fundamental theory of nanobubble formation and induced force in confined nanochannels. It is shown that nanobubble formation between hydrophobic plates can be predicted from their surface tension and geometry, with estimated values for the surface free energy and the force acting on the plates in good agreement with the results of molecular dynamics simulation and experimentation. When a bubble is formed between two plates, vertical attractive force and horizontal retract force due to the shifted plates are applied to the plates. The net force exerted on the plates is not dependent on the distance between them. The short-range force between hydrophobic surfaces due to hydrophobic interaction appears to correspond to the force estimated by our theory. We compared between experimental and theoretical values for the binding energy of a molecular motor system to validate our theory. The tendency that the binding energy increases as the size of the protein increases is consistent with the theory.