Elastic behavior of spherical nanodroplets in head-on collision

Simulation results for head-on collisions of equal-sized spherical polymer nanodroplets using molecular dynamics are presented. Elastic behavior of an initial compressed phase for the colliding droplets is analyzed. Deformations and contact radii of the nanodroplets are compared with the Hertzian model of elastic solid balls. It is found that at least the initial phase of collision can be explained by this continuum model, except at the very moment of the beginning of collision.

💡 Research Summary

The paper presents a systematic molecular dynamics (MD) investigation of head‑on collisions between equal‑sized spherical polymer nanodroplets, focusing on the very early, compression‑dominated stage of the impact. Each droplet is constructed from thousands of coarse‑grained beads linked by finitely extensible nonlinear elastic (FENE) springs, with non‑bonded interactions described by a Lennard‑Jones (LJ) potential. After equilibrating the droplets at a reduced temperature of 0.8 ε/k_B, they are propelled toward one another with identical velocities (v₀ ≈ 0.5 σ/τ) so that the initial center‑to‑center distance is just larger than twice the droplet radius. The simulations are performed using a velocity‑Verlet integrator with a timestep of 0.005 τ, and the full stress tensor, particle positions, and contact geometry are recorded throughout the impact.

The central question is whether the classical Hertzian contact theory—originally derived for elastic solid spheres—can describe the deformation and contact radius of nanodroplets that are, in principle, liquid‑like. To answer this, the authors extract the normal force F(t) acting between the droplets, the indentation depth δ(t) (the overlap of the two spherical caps), and the instantaneous contact radius a(t). According to Hertz theory, for two identical spheres of radius R, the relationships are F ∝ δ^{3/2} and a = √(R δ). When plotted on log–log axes, the MD data for the early compression phase follow a straight line with a slope of ≈1.5, confirming the F ∝ δ^{3/2} scaling. Moreover, the measured a(t) matches the Hertz prediction within a few percent, indicating that the droplets behave as if they possess an effective elastic modulus comparable to that of a solid.

Nevertheless, the agreement breaks down at the very instant the droplets first touch. In this pre‑contact regime the normal force rises more sharply than the Hertz law predicts, and the apparent contact radius is smaller than √(R δ). The authors attribute these deviations to surface tension (γ) and the discrete nature of the molecular interactions, which dominate when the contact area is only a few molecular diameters across. Once the compression proceeds beyond this initial “adhesion‑controlled” stage, the continuum description becomes valid.

Beyond the elastic stage, the simulations reveal rapid energy dissipation. The kinetic energy transferred into internal modes, chain stretching, and surface wave propagation leads to a marked reduction in rebound velocity. Consequently, the force–displacement curve departs from the Hertzian power law, and viscous effects become significant. The authors argue that while Hertzian mechanics can be used to model the initial shock‑absorption capability of nanodroplets, any comprehensive description of the full collision must incorporate viscoelastic or plastic dissipation mechanisms.

The discussion places these findings in a broader context. The authors identify three conditions under which a continuum elastic model is appropriate for nanoscale droplets: (1) the contact radius must be large enough to average out atomistic fluctuations, (2) the impact velocity must be moderate so that inertial and viscous forces do not overwhelm elastic response, and (3) surface tension effects must be relatively weak compared with the elastic restoring forces. When these criteria are met, the Hertzian framework provides a surprisingly accurate, computationally inexpensive tool for predicting contact geometry and peak forces in nanodroplet collisions.

In conclusion, the study demonstrates that the early compression phase of head‑on nanodroplet collisions can be quantitatively described by the Hertzian contact theory, despite the liquid nature of the droplets. This insight is valuable for designing nanoscale impact‑absorbing materials, micro‑capsule delivery systems, and for interpreting experimental observations of droplet coalescence at the nanometer scale. At the same time, the work highlights the limitations of purely elastic models: the onset of contact and the later, dissipative stages require explicit molecular‑level treatment. Future work is suggested to explore the influence of droplet composition, temperature, off‑axis collisions, and multi‑droplet interactions, thereby extending the applicability of the combined continuum‑molecular approach.