Abstract
Gaseous domains formed on solid–liquid interface have attracted scientists’ attentions in recent 2 decades, and the existence of interfacial nanobubble (INB) has been basically confirmed. However, an overall understanding on INB is still lacking. This research studied the influence of surface chemical heterogeneity on the morphology of INB by molecular dynamics simulations technique. The results showed that the gaseous domains could not nucleate on the hydrophilic substrate, while only dense gas layer (DGL) could be observed from the time-averaged density map for homogeneously hydrophobic substrate due to the random moving of INB. If there was a hydrophobic patch on the hydrophilic surface, INB could form on the hydrophobic patch with contact line being pinned at the boundary of the patch. In this case, the contact angle (gas-side) increased with the gas oversaturation degree and decreased with surface hydrophobicity of the patch. For the case that there existed a more hydrophobic patch/site on the hydrophobic surface, the INB could have moved randomly along the hydrophobic surface, but its receding contact line was pinned by the more hydrophobic patch/site. Hence, the INB could only move in the vicinity of this pinning patch/site, so that an INB profile instead of a DGL formed due to the pinning effect, and the apparent contact angle of the INB is significantly lower than the actual one. Throughout this study, the apparent INB we observed from experiments may be different from its instantaneous state and is significantly affected by surface heterogeneity.
Introduction
Interfacial nanobubble (INB) was first proposed by Parker et al. (1994) in 1994 due to the steps and the long-range attractive force in the force curve between two approaching hydrophobic surfaces. Until 2000, the first image of solid-water interfacial nanobubble was independently obtained through tapping mode of the atomic force microscopy by Lou et al. (2000) and Ishida et al. (2000). Since then, research on INBs have been extensively focused on due to their potential applications in many fields including froth flotation, protein adsorption, drag reduction, catalysis, and electrolysis (Lohse and Zhang, 2015).
INBs are usually referred to as gas aggregates that are confined in a spherical cap with height less than 100 nm and contact line diameter less than 1 μm (Alheshibri et al., 2016). According to Young–Laplace equation, the internal pressure of INB could be described aswhere P is the internal pressure of INB, is the gas–liquid interfacial tension, R is the curvature radius of the INB, is the ambient pressure, L is the width of INB, and is the contact angle of the INB. For an INB with L = 100 nm and (from gas side) in pure water, the calculated internal pressure of INB is p = 1.47 MPa = 14.7 atm, where . The huge internal pressure of INB would lead to the diffusive outflux of gas molecules, and the INB would dissolve at a time scale of according to the diffusion equation proposed by Epstein and Plesset (1950). Therefore, the research on INB had always been companied with the query about its existence or not (Ljunggren and Eriksson, 1997; Mao et al., 2004; McKee and Ducker, 2005; Takata et al., 2006). However, extensive studies in the last 2 decades demonstrated that INBs with long life-time indeed exist through various methods including total internal reflection fluorescence microscopy (Chan and Ohl, 2012; Tan et al., 2017), neutron reflectometry (Steitz et al., 2003), quartz crystal microbalance (Seo et al., 2007), attenuated total reflectance infrared spectroscopy (Zhang et al., 2008), and surface plasmon resonance (Martinez and Stroeve, 2007). In addition, INBs could exist on various substrates, including hydrophobic surfaces, such as HOPG (Fang et al., 2016; Lu et al., 2014), polystyrene (Li et al., 2016), molybdenite (Wang et al., 2019), OTS-coated Si (Zhang et al., 2006), and decanethiol-coated gold (Zhang, 2008), as well as hydrophilic surfaces, such as mica (Lou et al., 2000; Wang et al., 2019), gold (Holmberg et al., 2003), and metals (Cavicchi and Avedisian, 2007). Throughout these studies, INBs are found to have abnormal contact angle, which is much lower than the contact angle derived from Young’s equation. INBs usually form on hydrophobic surface, but it is beyond understanding why they could also form on the hydrophilic substrates, as well as why they have various contact angles even on the same surface (Lohse and Zhang, 2015). To data, the nature of INB has not been fully understood yet.
Molecular dynamics (MD) simulation, as a state-of-the-art technology to simulate the nanoscale phenomenon, is an ideal method to research the nucleation and dynamics of INBs from the molecular-level insights (Pu et al., 2020). By using MD simulation, many discoveries that are hardly acquired from experimental techniques were made. Dammer and Lohse (2006) discovered the gas enrichment phenomenon at liquid–solid interface. Xiao et al. (2017) found that the MD method could simulate the process of solvent exchange that have been widely used to produce INBs in experiments. Weijs et al. (2012) researched the formation of INBs on the homogeneous surface and found that the INB could form under a low gas solubility. Maheshwari et al. (2016), Maheshwari et al. (2018), Maheshwari et al. (2020) studied the stability of INBs and found that the INBs could be pinned by the more hydrophobic part, and the INB nucleation process is affected by the surface morphology. Recently, Yen et al. (2021) demonstrated that the gas molecules might migrate from gas enrichment layer to the inside of INB through the study in potential of mean force. Wang et al. (2022) researched the interfacial phenomenon in cement fluidity and found the interfacial structure of water is influenced by the shearing fluid, which is significant in many interfacial systems such as boundary slip. However, seldom researches have been conducted to study the formation and morphology of INBs on different types of substrates.
In this study, we employed MD simulations to study how INBs are influenced by the surface heterogeneity through the formation of INB on various types of substrates including homogeneous surface with different hydrophobicity, hydrophilic surface patterned with hydrophobic patch, and medium hydrophobic surface patterned with more hydrophobic patch/site. The current work facilitates understanding the interfacial phenomenon of surface nanobubble or nanodroplet and how it affects the morphology we observed by analytical instruments.
Simulation Methodology
MD simulations were performed with the open source code GROMACS 2019.6 package (Abraham et al., 2015) and were visualized by VMD software (Humphrey et al., 1996). Like many previous studies (Lu, 2019; Maheshwari et al., 2016; Maheshwari et al., 2020), this work employed Lennard–Jones beads to represent the liquid (L), solid (S), and gas (G) particles. Solid particles were fixed in a fcc lattice, while liquid and gas particles could move freely. The expression of LJ potential used to describe the interaction between different particles is shown as follows:where , , and is the interaction strength, characteristic size, and separation distance between particles i and j. The LJ parameters of various particles, shown in Table 1, were determined by referring to previous reports (Lu, 2019; Maheshwari et al., 2020), where SL, SM, and SH represent low, medium, and high hydrophobic surfaces, respectively. A cutoff radius of 1.2 nm was used to compute the LJ potential. The leap-frog algorithm was used to integrate the equation of motion, with a time step of 2 fs. Periodic boundary conditions were used in all three directions. The simulations were conducted in NVT ensemble for 1 ns where the system volume was unchanged and the temperature was kept at 300 K, and then were conducted in the NPzT ensemble where the pressure maintained at 1 bar. Semi-isotropic pressure coupling was used, and the simulation box could scale in Z direction to maintain constant pressure. The substrate was divided into two parts, as shown in red/orange and silver colors in Figure 1. The upper part of the solid was the surface we studied, while the bottom part was always hydrophilic particles (SL). Supplementary Figure S1 shows the typical initial configurations of Figures 2, 8. In the initial configuration, the liquid and gas particles were randomly dispersed above the solid substrate, and the gas particles were mainly arranged near the surface to aid the bubble nucleation for decreasing the equilibrium time. Supplementary Figure S2 shows the evolution of potential, temperature, and pressure of the system in Figure 2 with time, which clearly indicates that the simulation system already reaches equilibrium before 10 ns. In this study, coordinates were output every 1,000 steps, and the typical simulation time was 120 ns, unless otherwise specified.
TABLE 1
| Molecules | (nm) | (kJ/mol) |
|---|---|---|
| L-L | 0.34 | 3.00 |
| G-G | 0.50 | 1.00 |
| L-G | 0.41 | 1.73 |
| SL-L | 0.34 | 2.10 |
| SL-G | 0.40 | 1.56 |
| SM-L | 0.34 | 1.73 |
| SM-G | 0.40 | 1.73 |
| SH-L | 0.34 | 1.56 |
| SH-G | 0.40 | 2.10 |
LJ parameters used in simulations.
FIGURE 1
Typical snapshots of the equilibrated box for various substrates. (A) SL, (B) SM, (C) SH. In this case, 300 gas particles and 4,000 liquid particles were randomly filled into the simulation box with a typical size of 10.20 × 2.75 × 11.00 nm3 in x-, y-, and z-direction. The blue and green beads represent the liquid and gas particles, respectively, and the silver, orange, and red beads represent the fixed solid particles of SL, SM and SH.
FIGURE 2
Snapshots of simulation for SM surface. (A) 20.1 ns, (B) 21.0 ns, (C) 24.5 ns. In this case, 300 gas particles and 4,000 liquid particles were randomly filled into the simulation box with a typical size of 10.20 × 2.75 × 11.00 nm3 in x-, y-, and z-direction. The blue and green beads represent the liquid and gas particles, respectively, and the silver and orange beads represent the fixed solid particles of SL and SM. A movie for the process is available as Supplementary Video S1.
After the simulation, both the average densities of gas particles and of liquid particles in each bin of 0.1 × 2.75 × 0.1 nm3 within a certain simulation time range were calculated. The obtained density distribution of gas particles was presented as a color map such as Figure 3. The density distribution of liquid particles was analyzed to calculate the contact angle. The Z coordinate of mass center of the top layer of solid plus was considered as the solid–liquid interface. The gas–liquid interface was defined as the position where the liquid density equals half of the bulk liquid density. The height, radius, and contact angle of the INBs were obtained by fitting the gas–liquid interface to a circle function.
FIGURE 3
Time-averaged gas density distribution color maps (20–120 ns). (A) SM surface, (B) SH surface. The DGL profile were seen on both SM and SH substrates.
Results and Discussion
Gaseous Domains on Homogeneous Surfaces
Figure 1 shows the simulation snapshots of SL, SM, and SH. We could see an ordered layer of liquid particles formed adjacent to the SL surface. Bulk nanobubble is nucleated due to the gas oversaturation, but they could not adhere to the substrate because the hydrated layer on the hydrophilic surface is difficult to rupture. For SH surface, a dense gas layer (DGL) rather than an INB forms at the solid–liquid interface, which could also be seen from the averaged density map of gas molecules shown in Figure 3B. Only on a medium hydrophobic surface (SM substrate) could INB be formed.
When we inspected the trajectory of simulation for SM surface, we found that the INB moved along the solid–liquid interface randomly all the time. Figures 2A–C show the snapshots of INB formed on SM surface at three different simulation moments. The movie for the process is also available as Supplementary Video S1. The INB located at different positions at different moments. A similar phenomenon was also reported by Wu et al. (2017) who proposed that this random movement of INB was caused by the thermal fluctuation, and the movement frequency increased with the surface hydrophobicity. We believe Brownian movement might be also responsible for this phenomenon. When we averaged the density distribution of gas molecules over a time range of 100 ns, as shown in Figure 3A, we can see only a DGL profile instead of an INB on SM surface due to the rapid movement of INB, which is a visual effect from time-average rather than the instantaneous situation shown in Figure 2. Considering that the movement of INB is much faster than the temporal resolution of experimental instruments as that the former is at nanosecond magnitude while the latter is usually larger than millisecond for imaging INBs, this movement of INB is expected to influence the morphology and the calculated contact angle of the INB in experiments.
Gaseous Domains on Heterogeneous Surfaces
Hydrophilic Substrate Patterned by Hydrophobic Patch
The ideal homogeneous surface without pinning site rarely exists, and almost all the substrates coexist with heterogeneous pinning sites more or less. For hydrophilic substrate, it is usually hard to totally avoid the existence of hydrophobic patches on the surface. Figures 4, 5 show the gas density distribution maps for the hydrophilic substrates patterned by medium hydrophobic and strong hydrophobic patches, respectively, under different degrees of gas oversaturation. These results indicate that INBs could form on the hydrophobic patches patterned on the hydrophilic substrates. As shown in Figures 4, 5, the height and contact angle of INB increase with increasing gas particle number, while the width is unchanged due to contact line pinning, indicating that the height and contact angle of INBs are positively correlated with the degrees of gas oversaturation for INBs formed on the hydrophobic patches patterned on the hydrophilic substrate. However, if the degree of gas oversaturation is not high enough, such as 150 gas particles in Figure 4A or 200 gas particles in Figure 5A, a DGL instead of an INB would form, which is due to the rapid movement of small gas domain along hydrophobic patch just like what Figure 2 shows. In addition, the minimal gas particle numbers required to form an INB for SH patch is higher than that for SM patch, indicating that the minimal degree of gas oversaturation required to form INBs increases with the increase of hydrophobicity of the surface below INB.
FIGURE 4
Gas density distribution map averaged from 20 to 120 ns for SM patterned hydrophilic substrate. (A) 150, (B) 200, (C) 300, (D) 400 gas molecules.
FIGURE 5
Gas density distribution map averaged from 20 to 120 ns for SH patterned hydrophilic substrate. (A) 200, (B) 250, (C) 300, (D) 400 gas molecules.
Figures 6A,B show the density distribution curves along the central axis of the INBs. For medium hydrophobic patch as shown in Figure 6A, the intensities of the peaks in these curves increase with increasing gas particle number, while the peak intensities are almost unchanged with the increase of gas number for high hydrophobic pattern as shown in Figure 6B. It indicates that the gas density distribution next to the surface is affected by the degree of gas oversaturation for medium hydrophobic surface, but it is not affected by the gas oversaturation for strong hydrophobic surface. Figure 6C shows the contact angles of INBs formed on SM and SH patches when using different gas number. We can see that the contact angle increases with the increase of gas number, and the contact angles for SH are lower than those for SM, indicating that the contact angle of INBs increases with the gas oversaturation degree for a given surface hydrophobicity and decreases with the surface hydrophobicity for a given gas oversaturation. Figure 6D shows the number density of gas particles in the bulk increases with the increase of INB curvature, which might be due to the higher Laplace pressure for INB with a larger curvature.
FIGURE 6
Density distribution curves of gas particles along Z-axis and contact angles of INBs for hydrophilic substrates patterned by hydrophobic patch: (A) SM patch, (B) SH patch, (C) contact angles of INBs for SM and SH patches, (D) gas densities in bulk as a function of curvature of INBs. Note: the number of 400, 300, 250, 200, and 150 represent the number of gas molecules used in the simulations, where the liquid molecule number was 4,000 for all of these cases; the density distribution curves were obtained from the gas density of the center axis of INBs (averaged from 4.5 to 5.5 nm of X-axis).
Hydrophobic Substrate Patterned by Stronger Hydrophobic Patch
For a hydrophobic substrate, there might exist some patches with stronger hydrophobicity. Figure 7 shows the INBs formed on SM substrate patterned by SH patch. We could see that a cap-shaped INB forms on SH patterned SM surface, but the contact line stands on the SM surface rather than being pinned at boundary between SH and SM, which is not like the case of Figure 5. Although the contact line locates on the same surface (SM), contact angle is found increasing with the increase of gas oversaturation rather than keeping constant. The contact angles of Figures 7B–D are respectively 39.02°, 55.98°, and 66.76°, and their heights are 2.02, 2.93, and, 4.12 nm, respectively, namely the contact angle of INB increases with the increase of height, which is consistent with previous experiments (Zhang et al., 2012). The gas domain of Figure 7A is more like a DGL, which is because that the gas oversaturation degree is too low to form an INB. The gas enrichment layer is also found existing outside of the INB, which is also reported by previous researches (Yen, 2020) and is considered to account for the stability of INB (Tortora et al., 2020).
FIGURE 7
Time-averaged gas density distribution map for SH patterned SM substrate. (A) 400, (B) 500, (C) 600, (D) 800 gas molecules. There are 9,000 liquid molecules in these simulations, and the initial box size is 21.19 × 2.75 × 13.00 nm3. The yellow beads represent SM particles, the red beads represent SH particles, and the silver beads represent SL particles.
When we examined the simulation trajectory, we found that the advancing contact line of INB was not pinned by SH patch, but the receding contact line was always pinned. Hence, the INB could only move around the SH pattern, producing an apparent INB with contact line on SM surface and center on SH pattern. This phenomenon is more obvious when the SH pattern is diminished to a site, just as shown in Figure 8 where an INB forms on the SM surface patterned by SH site. It is apparent that the receding contact line of INB was always pinned by SH site, but the advancing contact line could not be pinned, so that the INB moved in a round area centered on SH site. These moving trails formed an apparent INB shown in Figure 9A. In this case, the apparent INB we observed is just a shadow rather than the actual nanobubble, and the time-averaged contact angle is also not the actual contact angle. Figure 9B shows a sketch for the influence of pinning site on contact angle of the INB. In this case, the white spherical cap represents the actual INB, while the red solid line represents the apparent INB. Obviously, the apparent contact angle we could see is much lower than the true value due to the pinning effect.
FIGURE 8
Snapshots of INBs formed on SH patterned SM substrate. (A) Snapshot of 20.1 ns represents the receding contact line was pinned when INB was moving toward the left; (B) 22.6 ns. (C) Snapshot of 24.8 ns represents the receding contact line was pinned when INB was moving toward the right. In this case, 600 gas particles and 15,000 liquid particles were added in the simulation boxes with the typical size of 21.19 × 2.75 × 17.00 nm3. A movie for the process is available as Supplementary Video S2.
FIGURE 9
(A) Time-averaged density distribution map of gas particles. (B) Sketch of INB pinned by the more hydrophobic site on the hydrophobic substrate. The yellow substrate represents the hydrophobic solid, and the red site represents the more hydrophobic solid. The white spherical cap represents the actual INB, and the red solid line represents the apparent INB we usually used to calculate contact angle.
For the case with several patches, there is similar phenomenon. Figures 10A,B show the INB formed on the SM substrate with two SH patches where gas particle number are 600 and 1,000, respectively. For the case of 600 gas particles, the INB was not large enough to occupy two patches. Hence, two INBs were formed on the two patches, respectively (Figure 10A). When the gas particle number increased to 1,000, a large INB crossing two patches was formed (Figure 10B). This INB kept moving randomly around the two patches, with receding contact line always being pinned by the patch (Supplementary Video S3). In this case, the role of the two small patches is similar to a large patch just like what Figure 7 shows, and the influence of the two patches on morphology of the INB is also like what Figure 9B illustrates.
FIGURE 10
Time-averaged gas density distribution map for substrate with two patches under gas density number of (A) 600 and (B) 1,000. Both the liquid particle number is 15,000. A movie for the process of (B) is available as Supplementary Video S3.
Conclusion
In this study, the influence of surface heterogeneity on interfacial nanobubble (INB) was studied by using MD method. The following conclusions were obtained:
(1) For smoothly homogeneous ideal-surface, the gaseous domain could not nucleate on the hydrophilic surface, and only dense gas layer (DGL) forms on strong hydrophobic surface. For medium hydrophobic surface, although INB could nucleate on the substrate, it moves freely along the solid–liquid interface, so that the time-averaged density distribution map of gas particles exhibits a DGL shape.
(2) For hydrophilic substrate patterned by hydrophobic patch, the INB could form on the patch with the contact line being pinned at the boundary. The contact angle of INB increases with gas oversaturation degree, and decreases with the surface hydrophobicity of the patch. In addition, the minimal gas oversaturation degree required to form an INB increases with the surface hydrophobicity of the patch. If the degree of gas oversaturation is lower than this required minimal value, a DGL instead of an INB would form.
(3) For hydrophobic substrate patterned by stronger hydrophobic patch/site, the INB could also form. We found that the receding contact line rather than advancing contact line was always pinned by the boundary of the patch, so that the INB randomly moved around the patch. Hence, a cap-shaped INB centered on the patch was seen from the averaged gas density map. For this case, the width of the apparent INB was much larger than that of instantaneous state, and thus the contact angle of the apparent INB was lower than the actual one. If the INB is large enough to occupy several patches, these patches behave like a whole large patch.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
HY conceived and preformed experiments, analyzed the data, and wrote the article. FZ analyzed the data. YX conceived the experiments. XG and YC wrote and edited the article.
Funding
This work was supported by the National Nature Science Foundation of China (Grant no. 51920105007), the Jiangsu Natural Science Fund–Youth Fund (BK20190639), and National Key R&D Program of China (2020YFC1908803).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmats.2021.824125/full#supplementary-material
References
1
AbrahamM. J.MurtolaT.SchulzR.PállS.SmithJ. C.HessB.et al (2: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. Softwarex1-2 (C), 19–25. 10.1016/j.softx.2015.06.001
2
AlheshibriM.QianJ.JehanninM.CraigV. S. J. (2016). A History of Nanobubbles. Langmuir32 (43), 11086–11100. 10.1021/acs.langmuir.6b02489
3
CavicchiR. E.AvedisianC. T. (2007). Bubble Nucleation and Growth Anomaly for a Hydrophilic Microheater Attributed to Metastable Nanobubbles. Phys. Rev. Lett.98 (12), 124501. 10.1103/PhysRevLett.98.124501
4
ChanC. U.OhlC.-D. (2012). Total-internal-reflection-fluorescence Microscopy for the Study of Nanobubble Dynamics. Phys. Rev. Lett.109 (17), 174501. 10.1103/PhysRevLett.109.174501
5
DammerS. M.LohseD. (2006). Gas Enrichment at Liquid-Wall Interfaces. Phys. Rev. Lett.96 (20), 206101. 10.1103/PhysRevLett.96.206101
6
EpsteinP. S.PlessetM. S. (1950). On the Stability of Gas Bubbles in Liquid‐Gas Solutions. J. Chem. Phys.18, 1505–1509. 10.1063/1.174818210.1063/1.1747520
7
FangC.-K.KoH.-C.YangC.-W.LuY.-H.HwangI.-S. (2016). Nucleation Processes of Nanobubbles at a Solid/water Interface. Sci. Rep.6, 24651. 10.1038/srep24651
8
HolmbergM.KühleA.Jørgen GarnæsJ.MørchK. A.BoisenA. (2003). Nanobubble Trouble on Gold Surfaces. Langmuir19 (25), 10510–10513. 10.1021/la0352669
9
HumphreyW.DalkeA.SchultenK. (1996). VMD: Visual Molecular Dynamics. J. Mol. Graphics14 (1), 33–38. 10.1016/0263-7855(96)00018-5
10
IshidaN.InoueT.MiyaharaM.HigashitaniK. (2000). Nano Bubbles on a Hydrophobic Surface in Water Observed by Tapping-Mode Atomic Force Microscopy. Langmuir16 (16), 6377–6380. 10.1021/la000219r
11
LiD.JingD.PanY.BhushanB.ZhaoX. (2016). Study of the Relationship between Boundary Slip and Nanobubbles on a Smooth Hydrophobic Surface. Langmuir32 (43), 11287–11294. 10.1021/acs.langmuir.6b02877
12
LjunggrenS.ErikssonJ. C. (1997). The Lifetime of a Colloid-Sized Gas Bubble in Water and the Cause of the Hydrophobic Attraction. Colloids Surf. A: Physicochemical Eng. Aspects129-130 (30), 151–155. 10.1016/S0927-7757(97)00033-2
13
LohseD.ZhangX. (2015). Surface Nanobubbles and Nanodroplets. Rev. Mod. Phys.87 (3), 981–1035. 10.1103/RevModPhys.87.981
14
LouS.-T.OuyangZ.-Q.ZhangY.LiX.-J.HuJ.LiM.-Q.et al (2000). Nanobubbles on Solid Surface Imaged by Atomic Force Microscopy. J. Vac. Sci. Technol. B18 (5), 2573–2575. 10.1116/1.1289925
15
LuY.-H.YangC.-W.FangC.-K.KoH.-C.HwangI.-S. (2014). Interface-induced Ordering of Gas Molecules Confined in a Small Space. Sci. Rep.4, 7189. 10.1038/srep07189
16
LuY. (2019). Drag Reduction by Nanobubble Clusters as Affected by Surface Wettability and Flow Velocity: Molecular Dynamics Simulation. Tribology Int.137, 267–273. 10.1016/j.triboint.2019.05.010
17
MaheshwariS.van der HoefM.Rodrı́guez Rodrı́guezJ.LohseD. (2018). Leakiness of Pinned Neighboring Surface Nanobubbles Induced by Strong Gas-Surface Interaction. ACS Nano12 (3), 2603–2609. 10.1021/acsnano.7b08614
18
MaheshwariS.van der HoefM.ZhangX.LohseD. (2016). Stability of Surface Nanobubbles: A Molecular Dynamics Study. Langmuir32 (43), 11116–11122. 10.1021/acs.langmuir.6b00963
19
MaheshwariS.van KruijsdijkC.SanyalS.HarveyA. D. (2020). Nucleation and Growth of a Nanobubble on Rough Surfaces. Langmuir36 (15), 4108–4115. 10.1021/acs.langmuir.0c00635
20
MaoM.ZhangJ.YoonR.-H.DuckerW. A. (2004). Is There a Thin Film of Air at the Interface between Water and Smooth Hydrophobic Solids? Langmuir20 (10), 4310. 10.1021/la040053k
21
MartinezJ.StroeveP. (2007). Transient Behavior of the Hydrophobic Surface/Water Interface: From Nanobubbles to Organic Layer. J. Phys. Chem. B111 (51), 14069–14072. 10.1021/jp077110d
22
McKeeC. T.DuckerW. A. (2005). Refractive Index of Thin, Aqueous Films between Hydrophobic Surfaces Studied Using Evanescent Wave Atomic Force Microscopy. Langmuir21 (26), 12153–12159. 10.1021/la051008v
23
ParkerJ. L.ClaessonP. M.AttardP. (1994). Bubbles, Cavities, and the Long-Ranged Attraction between Hydrophobic Surfaces. J. Phys. Chem.98 (34), 8468–8480. 10.1021/j100085a029
24
PuJ. H.SunJ.WangW.WangH. S. (2020). Generation and Evolution of Nanobubbles on Heated Nanoparticles: A Molecular Dynamics Study. Langmuir36 (9), 2375–2382. 10.1021/acs.langmuir.9b03715
25
SeoH.YooM.JeonS. (2007). Influence of Nanobubbles on the Adsorption of Nanoparticles. Langmuir23 (4), 1623–1625. 10.1021/la062763r
26
SteitzR.GutberletT.HaussT.KlösgenB.KrastevR.SchemmelS.et al (2003). Nanobubbles and Their Precursor Layer at the Interface of Water against a Hydrophobic Substrate. Langmuir19 (6), 2409–2418. 10.1021/la026731p
27
TakataY.ChoJ.-H. J.LawB. M.AratonoM. (2006). Ellipsometric Search for Vapor Layers at Liquid-Hydrophobic Solid Surfaces. Langmuir22 (4), 1715–1721. 10.1021/la052599s
28
TanB. H.AnH.OhlC.-D. (2017). Resolving the Pinning Force of Nanobubbles with Optical Microscopy. Phys. Rev. Lett.118 (5), 054501. 10.1103/PhysRevLett.118.054501
29
TortoraM.MeloniS.TanB. H.GiacomelloA.OhlC.-D.CasciolaC. M. (2020). The Interplay Among Gas, Liquid and Solid Interactions Determines the Stability of Surface Nanobubbles. Nanoscale12, 22698–22709. 10.1039/d0nr05859a
30
WangM.ZhangK.JiX.WangP.MaH.ZhangJ.et al (2022). Molecular Insight into the Fluidity of Cement Pastes: Nano-Boundary Lubrication of Cementitious Materials. Construction Building Mater.316, 125800. 10.1016/j.conbuildmat.2021.125800
31
WangY.LuoX.QinW.JiaoF. (2019). New Insights into the Contact Angle and Formation Process of Nanobubbles Based on Line Tension and Pinning. Appl. Surf. Sci.481, 1585–1594. 10.1016/j.apsusc.2019.01.292
32
WeijsJ. H.SnoeijerJ. H.LohseD. (2012). Formation of Surface Nanobubbles and the Universality of Their Contact Angles: a Molecular Dynamics Approach. Phys. Rev. Lett.108 (10), 104501. 10.1103/PhysRevLett.108.104501
33
WuC.-J.ChuK.-C.ShengY.-J.TsaoH.-K. (2017). Sliding Dynamic Behavior of a Nanobubble on a Surface. J. Phys. Chem. C121 (33), 17932–17940. 10.1021/acs.jpcc.7b04924
34
XiaoQ.LiuY.GuoZ.LiuZ.LohseD.ZhangX. (2017). Solvent Exchange Leading to Nanobubble Nucleation: A Molecular Dynamics Study. Langmuir33 (32), 8090–8096. 10.1021/acs.langmuir.7b01231
35
YenT.-H. (2020). Influence of Gas Aggregation on Water-Solid Interface: Molecular Simulation. Mol. Simulation46 (17), 1373–1382. 10.1080/08927022.2020.1828881
36
YenT.-H.LinC.-H.ChenY.-L. (2021). Effects of Gas Adsorption and Surface Conditions on Interfacial Nanobubbles. Langmuir37 (8), 2759–2770. 10.1021/acs.langmuir.0c03511
37
ZhangX. H.MaedaN.CraigV. S. J. (2006). Physical Properties of Nanobubbles on Hydrophobic Surfaces in Water and Aqueous Solutions. Langmuir22 (11), 5025–5035. 10.1021/la0601814
38
ZhangX. H. (2008). Quartz crystal Microbalance Study of the Interfacial Nanobubbles. Phys. Chem. Chem. Phys.10 (45), 6842–6848. 10.1039/b810587a
39
ZhangX. H.QuinnA.DuckerW. A. (2008). Nanobubbles at the Interface between Water and a Hydrophobic Solid. Langmuir24 (9), 4756–4764. 10.1021/la703475q
40
ZhangX.UddinM. H.YangH.ToikkaG.DuckerW.MaedaN. (2012). Effects of Surfactants on the Formation and the Stability of Interfacial Nanobubbles. Langmuir28 (28), 10471–10477. 10.1021/la301851g
Summary
Keywords
interfacial nanobubble, pinning, molecular dynamics, contact angle, surface heterogeneity
Citation
Yang H, Zhang F, Xing Y, Gui X and Cao Y (2022) Influence of Surface Heterogeneity on Morphology of Interfacial Nanobubble. Front. Mater. 8:824125. doi: 10.3389/fmats.2021.824125
Received
29 November 2021
Accepted
28 December 2021
Published
01 March 2022
Volume
8 - 2021
Edited by
Lijuan Zhang, Shanghai Advanced Research Institute (CAS), China
Reviewed by
Youguo Yan, China University of Petroleum, Huadong, China
Muhan Wang, Qingdao University of Technology, China
Updates
Copyright
© 2022 Yang, Zhang, Xing, Gui and Cao.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yaowen Xing, cumtxyw@126.com; Yijun Cao, yijuncao@126.com
This article was submitted to Colloidal Materials and Interfaces, a section of the journal Frontiers in Materials
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.