Coalescence of Polymer Droplets Moving on a Surface with Stiffness Gradient
Pith reviewed 2026-05-20 00:11 UTC · model grok-4.3
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.
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
- 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
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.
Referee Report
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)
- [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.
- [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)
- [Abstract] Abstract contains a typographical error: 'viscoelasticty' should read 'viscoelasticity'.
- [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
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
-
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
-
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
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
free parameters (2)
- velocity ratio of leading to trailing droplet
- attractive strengths for droplet-droplet and intra-droplet interactions
axioms (2)
- domain assumption Droplet motion is induced by the surface stiffness gradient
- domain assumption Bridge height growth obeys a power-law form before saturating
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation 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.leanJ_uniquely_calibrated_via_higher_derivative unclear?
unclearRelation 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
-
[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]
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)
work page 2025
-
[3]
C. L. Tucker, III and P. Moldenaers, Annual Review of Fluid Mechanics34, 177 (2002)
work page 2002
-
[4]
J. Kamp, J. Villwock, and M. Kraume, Reviews in Chemical Engineering33, 1–47 (2016)
work page 2016
- [5]
-
[6]
T. K. Pradhan and P. K. Panigrahi, Colloids and Surfaces A: Physicochemical and Engineering Aspects608, 125555 (2021)
work page 2021
-
[7]
B. Hafskjold, T. B. Morrow, H. K. B. Celius, and D. R. Johnson, inSPE Annual Technical Conference and Exhibition, 94SPE (SPE, 1994)
work page 1994
-
[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]
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)
work page 2020
-
[10]
W. P¨ onisch, C. A. Weber, G. Juckeland, N. Biais, and V. Zaburdaev, New Journal of Physics 19, 015003 (2017)
work page 2017
-
[11]
R. A. Foty and M. S. Steinberg, Developmental Biology278, 255–263 (2005)
work page 2005
- [12]
-
[13]
S. Sart, A.-C. Tsai, Y. Li, and T. Ma, Tissue Engineering Part B: Reviews20, 365–380 (2014). 21
work page 2014
-
[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)
work page 2009
-
[15]
N. Jain, J. Moeller, and V. Vogel, Annual Review of Biomedical Engineering21, 267–297 (2019)
work page 2019
-
[16]
W. Li, X. Tang, and L. Wang, Science Advances6, 10.1126/sciadv.abc1693 (2020)
-
[17]
G. M. Whitesides, E. Ostuni, S. Takayama, X. Jiang, and D. E. Ingber, Annual Review of Biomedical Engineering3, 335–373 (2001)
work page 2001
-
[18]
H. Geng, J. Feng, L. M. Stabryla, and S. K. Cho, Lab on a Chip17, 1060–1068 (2017)
work page 2017
-
[19]
H. Zhao, D. Orejon, K. Sefiane, and M. E. R. Shanahan, Langmuir41, 3986–3994 (2025)
work page 2025
- [20]
- [21]
- [22]
-
[23]
H. Yue, J. C. Burton, and D. M. Sussman, Physical Review Research6, 10.1103/physrevre- search.6.023115 (2024)
-
[24]
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]
S. Ryu, H. Zhang, and U. J. Anuta, Micromachines14, 2046 (2023)
work page 2046
-
[26]
H. P. Kavehpour, Annual Review of Fluid Mechanics47, 245–268 (2015)
work page 2015
-
[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]
V. Arumugam and A. Vasudevan, in4TH INTERNATIONAL CONFERENCE ON ROBOTICS, INTELLIGENT AUTOMATION AND CONTROL TECHNOLOGIES (RI- ACT2023), Vol. 3216 (AIP Publishing, 2024) p. 030001
work page 2024
-
[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]
R. W. Hopper, Journal of Fluid Mechanics213, 349–375 (1990)
work page 1990
-
[31]
L. DUCHEMIN, J. EGGERS, and C. JOSSERAND, Journal of Fluid Mechanics487, 167–178 (2003). 22
work page 2003
-
[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)
work page 2012
-
[33]
N. K. Chandra and A. Kumar, Soft Matter21, 3168–3183 (2025)
work page 2025
-
[34]
J. D. Paulsen, J. C. Burton, and S. R. Nagel, Physical Review Letters106, 10.1103/phys- revlett.106.114501 (2011)
-
[35]
J. D. Paulsen, Physical Review E88, 10.1103/physreve.88.063010 (2013)
-
[36]
A. U. Joshi, O. A. Basaran, and D. S. Corti, Physical Review E111, 10.1103/fpb9-pbcc (2025)
-
[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]
S. C. Varma, A. Saha, S. Mukherjee, A. Bandopadhyay, A. Kumar, and S. Chakraborty, Soft Matter16, 10921–10927 (2020)
work page 2020
-
[39]
M. R. Hassan, J. Zhang, and C. Wang, Langmuir37, 5823–5837 (2021)
work page 2021
-
[40]
Lohse, Annual Review of Fluid Mechanics54, 349–382 (2022)
D. Lohse, Annual Review of Fluid Mechanics54, 349–382 (2022)
work page 2022
- [41]
-
[42]
S. C. Varma, A. Saha, and A. Kumar, Physics of Fluids33, 10.1063/5.0073936 (2021)
-
[43]
S. Arbabi and P. E. Theodorakis, Macromolecular Theory and Simulations32, 10.1002/mats.202300017 (2023)
-
[44]
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]
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]
R. Roy, R. L. Seiler, J. A. Weibel, and S. V. Garimella, Advanced Materials Interfaces7, 10.1002/admi.202000731 (2020)
-
[47]
M. Sokuler, G. K. Auernhammer, M. Roth, C. Liu, E. Bonacurrso, and H.-J. Butt, Langmuir 26, 1544–1547 (2009)
work page 2009
-
[48]
R. Kajouri, P. E. Theodorakis, and A. Milchev, Langmuir40, 17779–17785 (2024)
work page 2024
-
[49]
P. E. Theodorakis, S. A. Egorov, and A. Milchev, The Journal of Chemical Physics146, 10.1063/1.4990436 (2017)
-
[50]
M. Palmiero, I. Cantarosso, L. di Blasio, V. Monica, B. Peracino, L. Primo, and A. Puliafito, Molecular Oncology17, 1699–1725 (2023). 23
work page 2023
- [51]
-
[52]
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]
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]
P. E. Theodorakis and A. Milchev, Advances in Colloid and Interface Science353, 103882 (2026)
work page 2026
-
[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]
-
[57]
S. K. Cho, H. Moon, and C.-J. Kim, Journal of Microelectromechanical Systems12, 70–80 (2003)
work page 2003
-
[58]
Knoche, Crop Protection13, 163–178 (1994)
M. Knoche, Crop Protection13, 163–178 (1994)
work page 1994
-
[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)
work page 2022
-
[60]
D. R. Heine, G. S. Grest, and E. B. Webb, Physical Review E70, 10.1103/physreve.70.011606 (2004)
-
[61]
J. S. Allen, Journal of Colloid and Interface Science261, 481–489 (2003)
work page 2003
-
[62]
V. A. Lubarda and K. A. Talke, Langmuir27, 10705–10713 (2011)
work page 2011
- [63]
-
[64]
P. E. Theodorakis, E. A. M¨ uller, R. V. Craster, and O. K. Matar, Soft Matter11, 9254 (2015)
work page 2015
-
[65]
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
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.