Investigating Interfacial Effects on Surface Nanobubbles without Pinning Using Molecular Dynamics Simulation

Mechanism of Charged Graphene Substrate Effects on the Stability of Interfacial Nanobubbles: Molecular Dynamics Simulations

Charged solid substrates play a crucial role in influencing the behavior of interfacial nanobubbles, although the underlying mechanisms are not yet fully understood. To explore this process in greater depth, we employed molecular dynamics (MD) simulations to systematically examine the effects of charged graphene on the morphological evolution, solid interface structure, and stability of interfacial nanobubbles, thereby revealing the intrinsic mechanisms. Our findings indicate that as surface charge density increases, the gas-solid interactions gradually diminish while the liquid-solid interactions significantly intensify. This results in a progressive reduction in both the contact angle and radius of the nanobubbles, eventually causing their detachment from the substrate and transformation to bulk-phase nanobubbles. Moreover, the enhanced gas accumulation effect at the solid interface leads to a reduction in the internal pressure of the bubbles, thus improving the stability of the interfacial nanobubbles. Additionally, the increase in the surface charge density elevates the water molecule density at the solid interface, which in turn strengthens the hydrogen bond network of interfacial water molecules, further stabilizing the liquid-solid interface structure. In summary, this study highlights the critical role of surface charge in regulating interfacial nanobubble behavior, providing new theoretical guidance for optimizing electrode materials and controlling bubble behavior in electrochemical systems.

Theoretical model of dynamics and stability of nanobubbles on heterogeneous surfaces

Surface nanobubbles have revealed a new mechanism of gas-liquid-solid interaction at the nanoscale; however, the nanobubble evolution on real substrates is still veiled, because the experimental observation of contact line motions at the nanoscale is too difficult.

HYPOTHESIS

This study proposes a theoretical model to describe the dynamics and stability of nanobubbles on heterogeneous substrates. It simultaneously considers the diffusive equilibrium of the liquid-gas interface and the mechanical equilibrium at the contact line, and introduces a surface energy function to express the substrate's heterogeneity.

VALIDATION

The present model unifies the nanoscale stability and the microscale instability of surface bubbles. The theoretical predictions are highly consistent to the nanobubble morphology on heterogeneous surfaces observed in experiments. As the nanobubbles grow, a lower Laplace pressure leads to weaker gas adsorption, and the mechanical equilibrium can eventually revert to the classical Young-Laplace equation above microscale.

FINDINGS

The analysis results indicate that both the decrease in substrate surface energy and the increase in gas oversaturation are more conducive to the nucleation and growth of surface nanobubbles, leading to larger stable sizes. The larger surface energy barriers result in the stronger pinning, which is beneficial for achieving stability of the pinned bubbles. The present model is able to reproduce the continual behaviors of the three-phase contact line during the nanobubble evolution, e.g., "pinning, depinning, slipping and jumping" induced by the nanoscale defects.

Molecular Dynamics Study on the Storage of Carbon Dioxide in Single-Walled Carbon Nanotubes at Low Pressures

The storage of carbon dioxide (CO2) in single-walled carbon nanotubes was studied with molecular dynamics simulation. The influences of the temperature, system average density, and nanotube size on the CO2 pressure, density distribution, and intermolecular forces were investigated. Multilayer adsorption inside nanotubes was observed as average density increases at lower pressures, which is desirable in industry. Meanwhile, a nanobubble was gradually formed in the center of the nanotube, and the system with the nanobubble was stabilized by the balance between the positive Laplace pressure and the negative liquid pressure when the size of the nanobubble was higher than the critical size. The adsorption effect of the nanotube wall leads to high local condensed density near the wall and stronger intermolecular repulsion, while Laplace pressure results in a low local condensed density in the adsorbed CO2 near the bubble interface and stronger intermolecular attraction. The stretching effect that originates from the intermolecular force dominated by attraction in the condensed phase leads to low pressure. At the critical nanobubble size, a higher CO2 average density can be achieved by lowering the temperature and increasing the nanotube radius or length. When the adsorption impact of the nanotube wall on bubble destabilization becomes negligible as the adsorption layer thickens, further increasing the nanotube radius leads to limited increase of the average density at the critical nanobubble size. The simulation of a graphene-sealed nanotube confirmed the formation of a vapor nanobubble under more realistic conditions. This work provides insights into utilizing carbon nanotubes as a material for CO2 capture with multilayer adsorption at lower pressures.

Molecular Dynamics Study of the Microscopic Mechanical Balance at the Three-Phase Contact Line of Interfacial Nanobubble

This study reveals the microscopic mechanical balance at the three-phase contact line (TPCL) of an interfacial nanobubble on a substrate with a wettability pattern using molecular dynamics simulations. The apparent contact angle was compared to that evaluated using Young's equation, in which the interfacial tensions were computed using a mechanical route. The comparison was conducted by changing the wettability of the substrate from hydrophilic to neutral while maintaining a hydrophobic region in the center of the substrate. When the wettability pattern pins the TPCL at the wettability boundary, the contact angle computed by Young's equation is larger than the apparent contact angle because a pinning force exists in the inward direction of the nanobubble. Conversely, on the surfaces where the wettability pattern does not pin the TPCL, the contact angle computed by Young's equation agrees with the apparent contact angle because the pinning force disappears. The distribution of principal stresses around the TPCL, which was visualized for the first time in this study, indicates that large compressive principal stresses exist between the liquid phase and the solid substrate interface, which pin the TPCL at the surface wettability boundary, and that the maximum principal stress occurs in the inward direction of the nanobubbles at the TPCL. The normalized pinning force estimated from the maximum principal stress is equivalent to that measured experimentally.

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.

Manipulating Trapped Nanobubbles Moving and Coalescing with Surface Nanobubbles

Trapped nanobubbles are observed nucleating at nanopits on a pitted substrate, while surface nanobubbles are usually formed on the smooth solid surface in water. In this work, trapped nanobubbles and surface nanobubbles were captured by a tapping-mode atomic force microscope (AFM) on a nanopitted substrate based on the temperature difference method. A single trapped nanobubble was manipulated to change into a surface nanobubble, then to change into the trapped nanobubble again. At the same time, surface nanobubbles can be moved to merge into a trapped nanobubble. Our results show that the scan load and the size of the scan area were the main factors that significantly affect the mobility of surface/trapped nanobubbles. The coalescence and mutual transformation of the two kinds of nanobubbles indicate that trapped nanobubbles and surface nanobubbles have the same chemical nature, which also provides vital experimental proof of the existence of nanobubbles in the course of contact line depinning. Our results are of great significance for understanding nanobubble stability and providing guidelines in some industrial applications.

Quantifying interfacial tensions of surface nanobubbles: How far can Young's equation explain?

Nanobubbles at solid-liquid interfaces play a key role in various physicochemical phenomena and it is crucial to understand their unique properties. However, little is known about their interfacial tensions due to the lack of reliable calculation methods. Based on mechanical and thermodynamic insights, we quantified for the first time the liquid-gas, solid-liquid, and solid-gas interfacial tensions of submicron-sized nitrogen bubbles at graphite-water interfaces using molecular dynamics (MD) analysis. It was revealed that Young's equation holds even for nanobubbles with different radii. We found that the liquid-gas and solid-liquid interfacial tensions were not largely affected by the gas density inside the nanobubbles. In contrast, the size effect on the solid-gas interfacial tension was observed, namely, the value dramatically decreased upon an increase in the gas density due to gas adsorption on the solid surface. However, our quantitative evaluation also revealed that the gas density effect on the contact angles is negligible when the footprint radius is larger than 50 nm, which is a typical range observed in experiments, and thus the flat shape and stabilization of submicron-sized surface bubbles observed in experiments cannot be explained only by the changes in interfacial tensions due to the van der Waals interaction-induced gas adsorption, namely by Young's equation without introducing the pinning effect. Based on our analysis, it was clarified that additional factors such as the differences in the studied systems are needed to explain the unresolved open issues - a satisfactory explanation for the nanobubbles in MD simulations being ultradense, non-flat, and stable without pinning.

Shock-induced collapse of surface nanobubbles

The collapse of cavitation bubbles often releases high-speed liquid jets capable of surface damage, with applications in drug delivery, cancer treatment, and surface cleaning. Spherical cap-shaped surface nanobubbles have previously been found to exist on immersed substrates. Despite being known nucleation sites for cavitation, their collapsing dynamics are currently unexplored. Here, we use molecular dynamics simulations to model the shock-induced collapse of different surface nanobubble sizes and contact angles. Comparisons are made with additional collapsing spherical nanobubble simulations near a substrate, to investigate the differences in their jet formation and resulting substrate pitting damage. Our main finding is that the pitting damage in the surface nanobubble simulations is greatly reduced, when compared to the spherical nanobubbles, which is primarily caused by the weaker jets formed during their collapse. Furthermore, the pit depths for surface nanobubble collapse do not depend on bubble size, unlike in the spherical nanobubble cases, but instead depend only on their contact angle. We also find a linear scaling relationship for all bubble cases between the final substrate damage and the peak pressure impulse at the impact centre, which can now be exploited to assess the relative damage in other computational studies of collapsing bubbles. We anticipate the more controlled surface-damage features produced by surface nanobubble cavitation jets will open up new applications in advanced manufacturing, medicine, and precision cleaning.

Recent advances for understanding the role of nanobubbles in particles flotation

Traditional froth flotation is the primary method for the separation and upgrading of fine mineral particles. However, it is still difficult for micro-fine and low-quality minerals to effectively separate. It is generally believed that bubble miniaturization is of great significance to improve flotation efficiency. Due to their unique physical and chemical properties, the application of nanobubbles (NBs) in ore flotation and other fields has been widely investigated as an important means to solve the problems of fine particle separation. Therefore, a fundamental understanding of the effect of NBs on flotation is a prerequisite to adapt it for the treatment of fine and low-quality minerals for separation. In this paper, recent advances in the field of nanobubble (NB) formation, preparation and stability are reviewed. In particular, we highlight the latest progress in the role of NBs on particles flotation and focus in particular on the particle-particle and particle-bubble interaction. A discussion of the current knowledge gap and future directions is provided.

Effects of Gas Adsorption and Surface Conditions on Interfacial Nanobubbles

Gas aggregation and formation of interfacial nanobubbles (INBs) provide challenges and opportunities in the operation of micro-/nanofluidic devices. In the current study, we used molecular dynamics(MD) simulations to investigate the effects of hydrophobicity and various homogeneous surface conditions on gas aggregation and INB stability with a series of 3D argon-water-solid and water-solid systems. Among various signatures of surface hydrophobicity, the potential of mean force (PMF) minima exhibited the strongest correlation with the water molecular orientation at the liquid-solid interface, compared to the depletion layer width and the droplet contact angle. Our results indicated that argon aggregation on the substrate was a function of hydrophobicity as well as competition between gas-solid and water-solid PMFs. Thus, one precondition for gas aggregation on a surface is that the free energy minima of gas induced by the surface be much lower than that induced by water. We found that although the presence of gas molecules had little effect on the measures of wettability, it enhanced density fluctuations near liquid-solid interfaces. The PMF of gas along the surface tangential plane exhibited a small energy barrier between the epitaxial gas layer (EGL) in the bubble and the gas enrichment layer (GEL) in the liquid, which may benefit nanobubble stability. Much lower PMF in the EGL compared to that in the GEL indicated that gas molecules could migrate from the GEL to the nanobubble basement. However, density fluctuations enhanced by the GEL could reduce the energy barrier, thus reducing the stability of INBs.

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.

Surface nanobubbles: Theory, simulation, and experiment. A review

Surface nanobubbles (NBs) are stable gaseous phases in liquids that form at the interface with solid substrates. They have been particularly intriguing for their high stability that contradicts theoretical expectations and their potential relevance for many technological applications. Here, we present the current state of the art in this research area by discussing and contrasting main results obtained from theory, simulation and experiment, and presenting their limitations. We also provide future perspectives anticipating that this review will stimulate further studies in the research area of surface NBs.

Study on the Formation and Properties of Trapped Nanobubbles and Surface Nanobubbles by Spontaneous and Temperature Difference Methods

Trapped nanobubbles are gas domains trapped at nanopits on the solid-liquid interface. This is different from surface nanobubbles that usually form at the smooth surface. Herein, both trapped nanobubbles and surface nanobubbles formed on the nanopitted polystyrene film were studied by a spontaneous formation method and a temperature difference method. Trapped nanobubbles behave more flexibly than surface nanobubbles under different scanning loads. The nanopits under trapped nanobubbles appear after being subjected to large force scanning, and both trapped nanobubbles and surface nanobubbles can recover after reducing the scan load. The contact angles of the two kinds of nanobubbles were calculated and were found to be approximately constant. Configurations of trapped nanobubbles including under the pit mouth, protruding out but pinning at the pit mouth, and protruding out and extending around the pit mouth were experimentally observed. Gas oversaturation in the liquid after replacing the low-temperature water with high-temperature water was evaluated and was found to be a key factor for nanobubble formation and led to trapped nanobubbles protruding out and extending. Our study should be helpful in understanding the formation mechanism and properties of trapped nanobubbles and surface nanobubbles, and it will also be useful for further research on the control of nanobubble distribution.

Stability of Surface Nanobubbles without Contact Line Pinning

Although the stability of most surface nanobubbles observed can be well interpreted by contact line pinning and supersaturation theory, there is increasing evidence that at least for certain situations, contact line pinning is not required for nanobubble stability. This raises a significant question of what is the stability mechanism for those sessile nanobubbles. Through theoretical analysis and molecular dynamics simulations, in this work, we report two mechanisms for stabilizing surface nanobubbles on flat and homogeneous substrates. One is attributed to constant adsorption of trace impurities on the nanobubble gas?liquid interface, through which nanobubble growing or shrinking causes the increase and decrease of interfacial tension, acting as a restoring force to bring the nanobubble to its equilibrium size. The other is attributed to the deformation of a soft substrate induced by the formed nanobubble, which in turn stabilizes the nanobubble via impeding the contact line motion, similar to self-pinning of microdroplets on soft substrates. Both mechanisms can interpret, depending on the specified conditions, how surface nanobubbles can remain stable in the absence of contact line pinning.

Stability, Dynamics, and Tolerance to Undersaturation of Surface Nanobubbles

The theoretical understanding of surface nanobubbles-nanoscale gaseous domains on immersed substrates-revolves around two contrasting perspectives. One perspective, which considers gas transport in the nanobubbles' vicinity, explains numerous stability-related properties but systematically underestimates the dynamical response timescale by orders of magnitude. The other perspective, which considers gas transport as the bulk liquid equilibrates with the external environment, recovers the experimentally observed dynamical timescale but incorrectly predicts that nanobubbles progressively shrink until dissolution. We propose a model that couples both perspectives, which is capable of explaining the stability, dynamics, and unexpected tolerance of surface nanobubbles to undersaturated environments.