pith. sign in

arxiv: 2605.18279 · v1 · pith:UI3OWUERnew · submitted 2026-05-18 · ❄️ cond-mat.soft · cond-mat.mtrl-sci· cond-mat.stat-mech

Coalescence of Polymer Droplets Moving on a Surface with Stiffness Gradient

Pith reviewed 2026-05-20 00:11 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.mtrl-scicond-mat.stat-mech
keywords polymer dropletscoalescencestiffness gradientdurotaxispower-law growthbridge heightcapillarityviscoelasticity
0
0 comments X

The pith

Droplet coalescence on stiffness gradients shows two distinct power-law regimes in bridge growth.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper examines the merging of two polymer droplets that are driven to move in the same direction by a gradient in surface stiffness. As a reference, it also looks at coalescence when the droplets are stationary on uniform stiffness surfaces. The main observation is that the height of the connecting bridge grows according to a power law in time, but the exponent changes from a higher value to a lower one after some time. This change marks the shift from a regime where surface tension (capillarity) dominates the process to one where the viscoelastic properties of the material take over. The speed difference between the two droplets and the strength of attractions within and between droplets also influence how the coalescence proceeds, with potential relevance to printing technologies and biological systems.

Core claim

The temporal growth of the bridge height (h) follows a power law (h ∼ t^α), with a transition from a higher to a lower value of α as a function of time, pointing to the presence of two distinct power-law growth regimes, where the transition signifies the crossover from the capillarity-dominated regime to the viscoelasticity-dominated regime of coalescence for durotaxis-driven droplet motion on a stiffness gradient surface.

What carries the argument

The power-law growth of bridge height h ∼ t^α and the time-dependent transition in the exponent α that indicates regime crossover.

If this is right

  • The coalescence process depends on the velocity ratio of the leading and trailing droplets.
  • Both the rate and extent of coalescence are affected by the attractive interaction strengths between and within droplets.
  • These findings apply to practical settings like microfluidics and ink-jet printing on substrates with varying stiffness.
  • Results contribute to understanding interactions among multicellular aggregates moving on biological surfaces with stiffness variations.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The observed crossover might be tunable by changing the stiffness gradient to control merging times in applications.
  • Similar power-law transitions could appear in other gradient-driven soft matter phenomena, such as cell migration or emulsion dynamics.
  • Experimental validation with real polymer droplets on fabricated stiffness-gradient surfaces could confirm the simulation results.

Load-bearing premise

The simulation model with the selected interaction potentials and velocity ratios accurately reproduces the physical mechanisms of durotaxis-driven coalescence without dominant numerical artifacts or unphysical boundary effects.

What would settle it

Direct measurement of bridge height versus time in an experiment with polymer droplets on a stiffness-gradient surface showing no transition in the growth exponent would falsify the claim of two distinct regimes.

Figures

Figures reproduced from arXiv: 2605.18279 by Divyansh Tripathi, Panagiotis E. Theodorakis, Swarn Lata Singh, Vimal Kishore.

Figure 1
Figure 1. Figure 1: FIG. 1: System setup consisting of two polymer droplets marked as [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Time evolution of droplet coalescence on a uniform substrate for different stiffness [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Temporal evolution of bridge height for stationary coalescence on a surface with [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Temporal evolution of bridge height for the coalescence of two droplets moving on [PITH_FULL_IMAGE:figures/full_fig_p013_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Time-sequenced images of coalescence of two moving droplets for [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: Difference of centers of mass of the coalescing droplets (∆ [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7 [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: Distance between centers of mass of two coalescing droplets on a gradient surface [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: Time sequenced snapshot for coalescence of moving droplets for [PITH_FULL_IMAGE:figures/full_fig_p018_9.png] view at source ↗
read the original abstract

Here, we study the coalescence of two droplets that are moving in the same direction on a soft surface; the motion of the droplets is caused by a gradient in the surface stiffness. As reference, stationary coalescence of the same droplets is also studied on the corresponding uniform surfaces for different stiffness values. To describe the coalescence phenomenon on a surface with stiffness gradient, a relevant range of velocity ratios of the leading and the trailing droplet was considered to elucidate the effect of this parameter on coalescence. Moreover, to analyze the dynamics of the process, the temporal growth of the bridge height $(h)$ was investigated, which follows a power law $(h \sim t^{\alpha})$, before eventually attaining a constant value. The obtained values of $\alpha$ show a transition from a higher to a lower value as a function of time, pointing to the presence of two distinct power-law growth regimes, where the transition signifies the crossover from the capillarity-dominated regime to the viscoelasticty-dominated regime of coalescence. In addition, varying attractive strengths for droplet--droplet and intra-droplet interactions were considered. The results indicate that both the dynamics and the degree of the coalescence strongly depend on these interaction parameters. Thus, we anticipate that our results will shed more light on the durotaxis-driven coalescence of polymeric droplets for various relevant system parameters, which will have practical implications for applications ranging from microfluidics to ink-jet printing, where substrate properties may vary. In addition, results may add to the fundamental understanding of the interactions among multicellular aggregates moving on biological surfaces.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript reports molecular dynamics simulations of two polymer droplets undergoing durotaxis-driven coalescence on a soft substrate with a stiffness gradient. As a reference, stationary coalescence is also examined on uniform-stiffness surfaces. The central observation is that the bridge height h grows temporally as a power law h ∼ t^α, with α exhibiting a transition from a higher early-time value to a lower late-time value; this transition is interpreted as the crossover from a capillarity-dominated regime to a viscoelasticity-dominated regime. The study further explores the dependence on the velocity ratio of the leading to trailing droplet and on the strengths of droplet–droplet and intra-droplet attractive interactions, reporting that both the coalescence dynamics and the final degree of merging are sensitive to these parameters.

Significance. If the reported power-law transition and its physical attribution are confirmed by scaling analysis and robustness checks, the work would add a useful parameter study to the literature on durotaxis and droplet coalescence on heterogeneous soft substrates. The explicit variation of velocity ratios and interaction strengths provides concrete, falsifiable trends that could guide experiments in microfluidics or ink-jet printing on gradient surfaces.

major comments (2)
  1. [Results section on temporal growth of bridge height] The central claim that the observed decrease in α with time marks the capillarity-to-viscoelasticity crossover is load-bearing for the interpretation. No scaling relation is shown between the measured crossover time and the ratio of capillary time (ηR/γ) to the viscoelastic relaxation time, nor is the late-time α compared quantitatively to the expected exponent for viscoelastic coalescence on a gradient substrate. Without these checks the attribution remains observational.
  2. [Methods and Results on power-law fitting] The simulation employs specific interaction potentials and an imposed velocity ratio. It is not demonstrated that the extracted α values and the transition are insensitive to the precise definition of bridge height h, to finite-size effects, or to the details of how the power-law fits are performed (e.g., fitting windows, error estimation). These controls are necessary to rule out that the transition is an artifact of the model or analysis choices.
minor comments (2)
  1. [Abstract] Abstract contains a typographical error: 'viscoelasticty' should read 'viscoelasticity'.
  2. [Methods] The manuscript would benefit from a brief statement of how the stiffness gradient is implemented numerically and from a table or figure summarizing the range of velocity ratios and interaction strengths explored.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments and positive assessment of the significance of our work. We address each major comment below and outline the revisions that will be incorporated to strengthen the scaling support and robustness analysis.

read point-by-point responses
  1. Referee: [Results section on temporal growth of bridge height] The central claim that the observed decrease in α with time marks the capillarity-to-viscoelasticity crossover is load-bearing for the interpretation. No scaling relation is shown between the measured crossover time and the ratio of capillary time (ηR/γ) to the viscoelastic relaxation time, nor is the late-time α compared quantitatively to the expected exponent for viscoelastic coalescence on a gradient substrate. Without these checks the attribution remains observational.

    Authors: We agree that a more explicit connection to the relevant time scales would strengthen the interpretation. In the revised manuscript we will add estimates of the capillary time ηR/γ using the effective viscosity and surface tension extracted from the simulations, and we will compare these estimates directly to the observed crossover time. For the late-time regime we will discuss the measured α in the context of known viscoelastic coalescence exponents, noting that the stiffness gradient and imposed velocity ratio are expected to yield a modified effective exponent relative to uniform-substrate cases. These additions will be placed in the Results section together with a brief scaling argument. revision: yes

  2. Referee: [Methods and Results on power-law fitting] The simulation employs specific interaction potentials and an imposed velocity ratio. It is not demonstrated that the extracted α values and the transition are insensitive to the precise definition of bridge height h, to finite-size effects, or to the details of how the power-law fits are performed (e.g., fitting windows, error estimation). These controls are necessary to rule out that the transition is an artifact of the model or analysis choices.

    Authors: We acknowledge the need for explicit robustness checks. In the revised manuscript and supplementary material we will (i) test alternative operational definitions of bridge height h (different density cut-offs and geometric neck measures) and demonstrate that the two-stage power-law behavior and the transition time remain consistent; (ii) report results from additional runs with larger lateral system sizes that confirm the transition is preserved; and (iii) describe the fitting protocol in the Methods, including the time windows employed, the criterion for selecting the early- and late-time regimes, and the uncertainty estimation obtained from multiple independent trajectories. These controls will be presented to show that the reported transition is not an artifact of the chosen analysis procedure. revision: yes

Circularity Check

0 steps flagged

Simulation observations of power-law regimes are self-contained with no circular reduction

full rationale

The paper reports results from molecular simulations of droplet coalescence driven by stiffness gradients, including direct measurement of bridge height h(t) and extraction of power-law exponents α in early and late regimes. These are outputs from varying input parameters (velocity ratios, interaction strengths) rather than any analytical derivation or fitted model that presupposes the capillarity-to-viscoelasticity crossover. No equations or claims reduce by construction to the simulation inputs, and the central interpretation follows from the observed trajectories without load-bearing self-citations or ansatzes imported from prior author work. The study is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central observations rest on simulation setup choices and standard assumptions about droplet physics rather than new derivations.

free parameters (2)
  • velocity ratio of leading to trailing droplet
    A relevant range was considered to study its effect on coalescence dynamics.
  • attractive strengths for droplet-droplet and intra-droplet interactions
    Varied to determine dependence of coalescence dynamics and degree.
axioms (2)
  • domain assumption Droplet motion is induced by the surface stiffness gradient
    Core setup of the moving-droplet case.
  • domain assumption Bridge height growth obeys a power-law form before saturating
    Basis for identifying the two regimes and crossover.

pith-pipeline@v0.9.0 · 5840 in / 1292 out tokens · 59015 ms · 2026-05-20T00:11:00.480160+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    the temporal growth of the bridge height (h) was investigated, which follows a power law (h∼t^α), before eventually attaining a constant value. The obtained values of α show a transition from a higher to a lower value as a function of time, pointing to the presence of two distinct power-law growth regimes, where the transition signifies the crossover from the capillarity-dominated regime to the viscoelasticty-dominated regime of coalescence.

  • IndisputableMonolith/Foundation/AlphaCoordinateFixation.lean J_uniquely_calibrated_via_higher_derivative unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    All non-bonded interactions ... modeled using the truncated and shifted Lennard-Jones (LJ) potential ... FENE potential ... harmonic interaction potential U_harmonic(r) = −½Kr²

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

65 extracted references · 65 canonical work pages

  1. [1]

    Y. Li, A. A. Pahlavan, Y. Chen, S. Liu, Y. Li, H. A. Stone, and S. Granick, Proceedings of the National Academy of Sciences120, 10.1073/pnas.2214657120 (2023)

  2. [2]

    T. Xiao, X. Ma, H. Hu, F. Xiang, X. Zhang, Y. Zheng, H. Dong, B. Adhikari, Q. Wang, and A. Shi, Food Chemistry: X29, 102792 (2025)

  3. [3]

    C. L. Tucker, III and P. Moldenaers, Annual Review of Fluid Mechanics34, 177 (2002)

  4. [4]

    J. Kamp, J. Villwock, and M. Kraume, Reviews in Chemical Engineering33, 1–47 (2016)

  5. [5]

    Sarig, Y

    I. Sarig, Y. Starosvetsky, and A. D. Gat, Journal of Fluid Mechanics800, 264–277 (2016)

  6. [6]

    T. K. Pradhan and P. K. Panigrahi, Colloids and Surfaces A: Physicochemical and Engineering Aspects608, 125555 (2021)

  7. [7]

    Hafskjold, T

    B. Hafskjold, T. B. Morrow, H. K. B. Celius, and D. R. Johnson, inSPE Annual Technical Conference and Exhibition, 94SPE (SPE, 1994)

  8. [8]

    A. D. Pizarro, C. L. A. Berli, G. J. A. A. Soler-Illia, and M. G. Bellino, Nature Communications 13, 10.1038/s41467-022-30834-2 (2022)

  9. [9]

    C. H. Meredith, P. G. Moerman, J. Groenewold, Y.-J. Chiu, W. K. Kegel, A. van Blaaderen, and L. D. Zarzar, Nature Chemistry12, 1136–1142 (2020)

  10. [10]

    P¨ onisch, C

    W. P¨ onisch, C. A. Weber, G. Juckeland, N. Biais, and V. Zaburdaev, New Journal of Physics 19, 015003 (2017)

  11. [11]

    R. A. Foty and M. S. Steinberg, Developmental Biology278, 255–263 (2005)

  12. [12]

    Oriola, M

    D. Oriola, M. Marin-Riera, K. Anla¸ s, N. Gritti, M. Sanaki-Matsumiya, G. Aalderink, M. Ebisuya, J. Sharpe, and V. Trivedi, Soft Matter18, 3771–3780 (2022)

  13. [13]

    Sart, A.-C

    S. Sart, A.-C. Tsai, Y. Li, and T. Ma, Tissue Engineering Part B: Reviews20, 365–380 (2014). 21

  14. [14]

    P. W. Oakes, D. C. Patel, N. A. Morin, D. P. Zitterbart, B. Fabry, J. S. Reichner, and J. X. Tang, Blood114, 1387–1395 (2009)

  15. [15]

    N. Jain, J. Moeller, and V. Vogel, Annual Review of Biomedical Engineering21, 267–297 (2019)

  16. [16]

    W. Li, X. Tang, and L. Wang, Science Advances6, 10.1126/sciadv.abc1693 (2020)

  17. [17]

    G. M. Whitesides, E. Ostuni, S. Takayama, X. Jiang, and D. E. Ingber, Annual Review of Biomedical Engineering3, 335–373 (2001)

  18. [18]

    H. Geng, J. Feng, L. M. Stabryla, and S. K. Cho, Lab on a Chip17, 1060–1068 (2017)

  19. [19]

    H. Zhao, D. Orejon, K. Sefiane, and M. E. R. Shanahan, Langmuir41, 3986–3994 (2025)

  20. [20]

    Eggers, J

    J. Eggers, J. E. Sprittles, and J. H. Snoeijer, Annual Review of Fluid Mechanics57, 61 (2025)

  21. [21]

    EGGERS, J

    J. EGGERS, J. R. LISTER, and H. A. STONE, Journal of Fluid Mechanics401, 293–310 (1999)

  22. [22]

    Eggers, J

    J. Eggers, J. E. Sprittles, and J. H. Snoeijer, Annual Review of Fluid Mechanics57, 61–87 (2025)

  23. [23]

    H. Yue, J. C. Burton, and D. M. Sussman, Physical Review Research6, 10.1103/physrevre- search.6.023115 (2024)

  24. [24]

    Menchaca-Rocha, A

    A. Menchaca-Rocha, A. Mart´ ınez-D´ avalos, R. N´ u˜ nez, S. Popinet, and S. Zaleski, Physical Review E63, 10.1103/physreve.63.046309 (2001)

  25. [25]

    S. Ryu, H. Zhang, and U. J. Anuta, Micromachines14, 2046 (2023)

  26. [26]

    H. P. Kavehpour, Annual Review of Fluid Mechanics47, 245–268 (2015)

  27. [27]

    D. G. A. L. Aarts, H. N. W. Lekkerkerker, H. Guo, G. H. Wegdam, and D. Bonn, Physical Review Letters95, 10.1103/physrevlett.95.164503 (2005)

  28. [28]

    Arumugam and A

    V. Arumugam and A. Vasudevan, in4TH INTERNATIONAL CONFERENCE ON ROBOTICS, INTELLIGENT AUTOMATION AND CONTROL TECHNOLOGIES (RI- ACT2023), Vol. 3216 (AIP Publishing, 2024) p. 030001

  29. [29]

    R. W. Hopper, Journal of the American Ceramic Society67, C-262 (1984), https://ceramics.onlinelibrary.wiley.com/doi/pdf/10.1111/j.1151-2916.1984.tb19692.x

  30. [30]

    R. W. Hopper, Journal of Fluid Mechanics213, 349–375 (1990)

  31. [31]

    DUCHEMIN, J

    L. DUCHEMIN, J. EGGERS, and C. JOSSERAND, Journal of Fluid Mechanics487, 167–178 (2003). 22

  32. [32]

    J. D. Paulsen, J. C. Burton, S. R. Nagel, S. Appathurai, M. T. Harris, and O. A. Basaran, Proceedings of the National Academy of Sciences109, 6857–6861 (2012)

  33. [33]

    N. K. Chandra and A. Kumar, Soft Matter21, 3168–3183 (2025)

  34. [34]

    J. D. Paulsen, J. C. Burton, and S. R. Nagel, Physical Review Letters106, 10.1103/phys- revlett.106.114501 (2011)

  35. [35]

    J. D. Paulsen, Physical Review E88, 10.1103/physreve.88.063010 (2013)

  36. [36]

    A. U. Joshi, O. A. Basaran, and D. S. Corti, Physical Review E111, 10.1103/fpb9-pbcc (2025)

  37. [37]

    P. J. Dekker, M. A. Hack, W. Tewes, C. Datt, A. Bouillant, and J. H. Snoeijer, Physical Review Letters128, 10.1103/physrevlett.128.028004 (2022)

  38. [38]

    S. C. Varma, A. Saha, S. Mukherjee, A. Bandopadhyay, A. Kumar, and S. Chakraborty, Soft Matter16, 10921–10927 (2020)

  39. [39]

    M. R. Hassan, J. Zhang, and C. Wang, Langmuir37, 5823–5837 (2021)

  40. [40]

    Lohse, Annual Review of Fluid Mechanics54, 349–382 (2022)

    D. Lohse, Annual Review of Fluid Mechanics54, 349–382 (2022)

  41. [41]

    Klestova, E

    A. Klestova, E. Sergeeva, and A. V. Vinogradov, Coatings9, 275 (2019)

  42. [42]

    S. C. Varma, A. Saha, and A. Kumar, Physics of Fluids33, 10.1063/5.0073936 (2021)

  43. [43]

    Arbabi and P

    S. Arbabi and P. E. Theodorakis, Macromolecular Theory and Simulations32, 10.1002/mats.202300017 (2023)

  44. [44]

    Rostami, A

    P. Rostami, A. Erb, R. Azizmalayeri, J. Steinmann, R. W. Stark, and G. K. Auernhammer, Physical Review Fluids10, 10.1103/physrevfluids.10.063603 (2025)

  45. [45]

    P. R. Kaneelil, K. Tojo, P. K. Farsoiya, L. Deike, and H. A. Stone, Journal of Fluid Mechanics 1026, 10.1017/jfm.2025.10878 (2026)

  46. [46]

    R. Roy, R. L. Seiler, J. A. Weibel, and S. V. Garimella, Advanced Materials Interfaces7, 10.1002/admi.202000731 (2020)

  47. [47]

    Sokuler, G

    M. Sokuler, G. K. Auernhammer, M. Roth, C. Liu, E. Bonacurrso, and H.-J. Butt, Langmuir 26, 1544–1547 (2009)

  48. [48]

    Kajouri, P

    R. Kajouri, P. E. Theodorakis, and A. Milchev, Langmuir40, 17779–17785 (2024)

  49. [49]

    P. E. Theodorakis, S. A. Egorov, and A. Milchev, The Journal of Chemical Physics146, 10.1063/1.4990436 (2017)

  50. [50]

    Palmiero, I

    M. Palmiero, I. Cantarosso, L. di Blasio, V. Monica, B. Peracino, L. Primo, and A. Puliafito, Molecular Oncology17, 1699–1725 (2023). 23

  51. [51]

    Kim and T

    J. Kim and T. Ma, Biotechnology Progress29, 441–451 (2013)

  52. [52]

    Porter, C

    R. Porter, C. Trenado-Yuste, A. Martinez-Calvo, M. Su, N. S. Wingreen, S. S. Datta, and K. C. Huang, Nature Reviews Physics 10.1038/s42254-025-00849-x (2025)

  53. [53]

    Welker, T

    A. Welker, T. Cronenberg, R. Z¨ ollner, C. Meel, K. Siewering, N. Bender, M. Hennes, E. R. Old- ewurtel, and B. Maier, Physical Review Letters121, 10.1103/physrevlett.121.118102 (2018)

  54. [54]

    P. E. Theodorakis and A. Milchev, Advances in Colloid and Interface Science353, 103882 (2026)

  55. [55]

    H. Zhu, Q. Deng, J. Li, L. Yang, H. Li, Z. Zhao, Z. Wang, C. Pang, Y. Zhang, V. C.-H. Lui, W. Li, X. Yin, and L. Wang, Science Advances11, 10.1126/sciadv.adv6314 (2025)

  56. [56]

    Alistar and U

    M. Alistar and U. Gaudenz, Bioengineering4, 45 (2017)

  57. [57]

    S. K. Cho, H. Moon, and C.-J. Kim, Journal of Microelectromechanical Systems12, 70–80 (2003)

  58. [58]

    Knoche, Crop Protection13, 163–178 (1994)

    M. Knoche, Crop Protection13, 163–178 (1994)

  59. [59]

    A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, Comp. Phys. Comm.271, 108171 (2022)

  60. [60]

    D. R. Heine, G. S. Grest, and E. B. Webb, Physical Review E70, 10.1103/physreve.70.011606 (2004)

  61. [61]

    J. S. Allen, Journal of Colloid and Interface Science261, 481–489 (2003)

  62. [62]

    V. A. Lubarda and K. A. Talke, Langmuir27, 10705–10713 (2011)

  63. [63]

    Iwamoto, H

    M. Iwamoto, H. Ino, T. Kimura, and Y. Kishimoto, Nano Studies8, 29 (2013)

  64. [64]

    P. E. Theodorakis, E. A. M¨ uller, R. V. Craster, and O. K. Matar, Soft Matter11, 9254 (2015)

  65. [65]

    Katre, B

    P. Katre, B. Mahapatra, M. Karmakar, S. Jagdish, N. K. Chandra, and A. Kumar, Proceedings of the National Academy of Sciences123, 10.1073/pnas.2528750123 (2026). 24