Molecular dynamics study of nanobubbles in the equilibrium Lennard-Jones fluid

We employ a model, in which the density fluctuations in a bulk liquid are represented as presence of the clusters of molecules with the lowered number of nearest neighbors (number of bonds). The nanobubble size distribution is calculated on the basis of a close analogy between the surface part of the work of formation for a cluster and for a nanobubble. The pre-exponential factor for this distribution is related to the fluid compressibility. Estimates made for different liquids show that it can be noticeably different from that adopted in the classical nucleation theory (CNT). Molecular dynamics (MD) simulation is performed for a liquid inside a macroscopic droplet of molecules interacting via the Lennard-Jones potential plus a long-range tail. The nanobubbles are identified by clusters of bond-deficient particles with the optimum number of bonds that provide the maximum nanobubble number density and maximum resolvable nanobubble equimolar size. The results of MD simulation are in qualitatively better agreement with proposed theory than with CNT.