REVIEW 4 major objections 5 minor 37 references
Partial-Wetting Phenomena in Active Matter
T0 review · 4 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read This paper demonstrates that partial wetting occurs in a nonequilibrium ternary laning system and verifies that the contact angles of the C-rich phase obey Young's equation, with interfacial tensions scaling as γ_αβ ∝ (Δq_αβ E)^0.23.
desk verdict Lens-shaped partial wetting in a ternary laning system is genuinely new; the Young's-equation 'verification' is a self-consistent check, not an independent test. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The argument runs through two stages. First, the Onsager variational principle, applied to driven particles at a fixed external field, maps the ternary laning system onto an equilibrium-like binary mixture: with fast relaxation along the field direction and slow lateral dynamics, each species forms drifting pseudo-particle columns, and the longitudinal dissipation from their relative velocities acts as an effective interaction χ, yielding an effective free energy à = Σ φ_α ln φ_α + χ φ_α φ_β. Second, Landau mean-field theory for interfaces is borrowed from equilibrium: the A–B interfacial profile is fitted to a tanh form, giving the interfacial width λ ∝ (Δq_AB E)^−0.23, and the relations γ = 2√2 c φ_0²/(3λ) and λ = (c/(2−χ))^{1/2} convert this width into γ ∝ (Δq E)^0.23. This scaling is then inserted into Young's equation for all three interfaces, producing closed-form predictions for the contact angles that depend only on q_C.
What would settle it
Measure the interfacial tension of the A–C and B–C interfaces directly (for example from the capillary-wave spectrum of their fluctuation profiles in the ternary simulation) and check whether both collapse onto γ = C(Δq E)^0.23 with the same C and exponent inferred from the A–B interface. Alternatively, simulate a system with the same Δq but a deliberately different C for one pair, such as by rescaling the species diameters, and test whether the contact angles still follow Eq. (6); if they change with E, the Young's-equation verification fails.
Extended reading notes
Core claim
The central claim is that partial-wetting phenomena in an active ternary laning system are governed by the same mechanical balance of interfacial tensions at the contact line as equilibrium wetting, i.e., by Young's equation. For a C-rich lens between A and B bulk phases, the force balance reads γ_AC cos θ_AC + γ_BC cos θ_BC = γ_AB and γ_AC sin θ_AC = γ_BC sin θ_BC. The paper further claims that each interfacial tension scales as γ_αβ = C(Δq_αβ E)^0.23 with a common constant C, so that, since Δq_AB = 2, Δq_AC = 1−q_C, Δq_BC = 1+q_C, the contact angles become functions of q_C only and are independent of E. Simulation measurements of θ_AC and θ_BC for q_C = 0, 0.3, and 0.5 across the partial-wetting window are reported to match the theoretical predictions (θ = 54°, θ_AC ≈ 59° and θ_BC ≈ 48°, and so on), including the discontinuous jumps to 0° and 180° at the complete-wetting and drying boundaries. The authors take this as quantitative verification of Young's equation in a nonequilibrium active system.
Load-bearing premise
The claim rests on assuming that the A–C and B–C interfaces obey the same power law γ = C(Δq E)^0.23 with the same constant C as the A–B interface; this was not measured directly, and if the prefactor or exponent differs between interfaces, the predicted contact angles and the verification of Young's equation would not follow.
Editorial extensions
If this is right
- Partial wetting in this active system is genuinely partial: C-rich droplets sit at the A–B interface with finite, tunable contact angles over an intermediate range of drive strengths, with well-defined drying and complete-wetting boundaries.
- Contact angles are set solely by the coupling charge q_C of the wetting species; changing the overall drive E moves the system between wetting regimes but does not alter the shape of the droplet within the partial-wetting window.
- The effective free energy of the laning system takes an equilibrium Flory–Huggins-like form, so equilibrium tools (spinodal and binodal lines, mean-field tension predictions) transfer to this driven system.
- Tuning q_C (or the relative drift speeds) provides a dynamic handle on interfacial tension, which the authors propose for designing responsive active coatings, microfluidic devices, and structured robotic swarms.
Reading between the lines
- A stronger and more direct test would measure γ_AC and γ_BC independently, for example from capillary-fluctuation spectra of their interfaces; the paper's verification rests on the assumption that all three tensions share one common prefactor C and the same exponent 0.23.
- The columnar lens geometry (translationally invariant along the field, straight contact lines) is a new testing ground for line-tension effects in nonequilibrium interfaces, since line tension would show up as a size dependence of the contact angles for small droplets.
- If the scaling survives direct measurement, the ratio γ_αβ/γ_AB is predicted to be a universal function of (Δq_αβ/Δq_AB)^0.23; checking this ratio across different total densities and species fractions would extend the claim beyond the single parameter set reported.
- The same tunable-wetting mechanism should appear in other driven multi-species systems, such as oppositely charged colloids in a field or sheared binary mixtures, offering a testable prediction of a nonequilibrium wetting phase diagram.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper studies a three-dimensional A-B-C laning system under an external drive, in which A and B with opposite charges phase-separate and a third component C accumulates at the A-B interface. It identifies drying, partial wetting, and complete wetting states as functions of E and qC, marks phase boundaries with order-parameter fluctuations and unsupervised PCA, and fits the A-B interface profile to a tanh form to obtain lambda proportional to (Delta q E)^-0.23. Using equilibrium mean-field relations gamma proportional to 1/lambda and lambda = (c/(2-chi))^{1/2}, it converts this scaling into gamma_alpha,beta proportional to (Delta q_alpha,beta E)^0.23 and inserts it into Young's equation, giving contact angles that depend only on qC. The simulated theta_AC and theta_BC are reported to agree with these predictions, and the paper claims quantitative validation of Young's equation in active matter.
Significance. The observation of three wetting states in a driven laning system and the apparent E-independence of the contact angles is a substantive result. The use of unsupervised machine learning to locate phase boundaries and the explicit tuning of three interfacial interactions through qC are valuable innovations. However, because the predicted curves are generated by the same scaling relation used to fit the binary A-B data, the comparison in the paper is a consistency check rather than an independent verification of Young's equation. If the common-prefactor tension scaling is independently confirmed by direct measurements of gamma_AC and gamma_BC, and if the mean-field gamma-lambda relation is corrected, the result would be a notable extension of Young's equation to a nonequilibrium active system.
major comments (4)
- [Main text, Eq. (6); SM Sec. 7] The central verification of Young's equation is not independent of the scaling hypothesis being tested. The exponent k = 0.23 is obtained by fitting the binary A-B interfacial width (Fig. 1(c), SM Sec. 5), and the same mean-field relation gamma_alpha,beta proportional to (Delta q_alpha,beta E)^0.23 with a single common prefactor is then substituted into Young's equation to obtain Eq. (6). Since gamma_AC and gamma_BC are never measured directly, the agreement in Fig. 3 and SM Fig. S6 is a consistency check between the contact-angle data and an assumed tension ratio, not a measurement of the tension balance. The authors should either measure the three interfacial tensions independently (for example from pressure tensors or capillary shapes) or present an explicit comparison of the tension ratios gamma_AC/gamma_BC against the assumed power-law form.
- [SM Sec. 7; main text Fig. 3(d)] The common-prefactor power law is applied to interface pairs for which no phase-separated interface exists in a binary reference system. At qC = 0.5 and E = 120, Delta q_AC E = 60, which is below the binary phase-separation threshold (Delta q E)_c approximately 76, so the A-C pair does not form a bulk interface and the relation gamma_AC proportional to (Delta q_AC E)^0.23 cannot be defined for that state point. The same issue arises for qC values near the lower partial-wetting boundary. The paper should either justify the extrapolation of the interfacial-tension scaling into the subcritical regime for the A-C pair or restrict the comparison to state points where all three interfaces exist in their respective binary systems.
- [Main text, interfacial-width discussion; SM Eq. (S44) and SM Sec. 5] The mean-field relation between the interfacial width and the effective interaction parameter is internally inconsistent as written. For the Flory-Huggins-like free energy in Eq. (4), phase separation requires chi > chi_s = 2, in which case the expression lambda = (c/(2-chi))^{1/2} in SM Eq. (S44) gives an imaginary width. In addition, the statement in SM Sec. 5 that chi is proportional to 2 - (Delta q E)^0.46 implies that chi decreases as the drive increases, which contradicts the phase diagram and the physical picture that stronger driving strengthens phase separation. The authors should correct the lambda-chi relation and re-derive gamma(lambda), or state explicitly that lambda is used only as an empirical fit and that chi is not inferred from Eq. (S44).
- [Main text, contact-angle predictions; SM Eq. (S53)] The theoretical curves compared with simulation are not parameter-free predictions. Both the exponent k and the common prefactor in the tension scaling are taken from the binary A-B fit and are then inserted into Young's equation. The paper should be explicit that Eq. (6) is a derived consequence of the fitted scaling, and it should discuss how a failure of the common-prefactor assumption would change the predicted angles. Without such a discussion, the claimed quantitative verification of Young's equation is overstated.
minor comments (5)
- [Abstract] The phrase "which described by Young's equation" should read "which is described by Young's equation."
- [Fig. 3 and SM Fig. S6] The figures do not report error bars or the number of independent runs used for the contact-angle averages; this information should be added so the claimed quantitative agreement can be assessed.
- [SM Sec. 5] The expression "chi proportional to 2 - (Delta q E)^0.46" is ambiguous; it should be replaced with a properly normalized equation with a constant and an explicit statement of whether chi increases or decreases with Delta q E.
- [Main text, Eq. (6)] The exponent k is used in Eq. (6) before it is defined; the text should state that k = 0.23 is the fitted exponent from Fig. 1(c).
- [SM Sec. 7] For qC = 0.5 the numerical predicted values theta_AC approximately 63.34 degrees and theta_BC approximately 43.96 degrees should be given in the caption of Fig. S6 as well as in the text.
Circularity Check
No significant circularity: the contact angles are independently measured simulation outputs, not fit inputs; the fitted exponent k=0.23 is transferred from the binary A-B interface to predict A-C/B-C angles, and the unverified common-prefactor scaling is a robustness limitation, not a circular reduction.
full rationale
The central derivation (Eq. 6 and SM Eqs. S50-S53) combines Young's equation (Eq. 5) with the scaling γαβ ∝ (Δqαβ E)^0.23, whose exponent k=0.23 is fit to the width of the binary A-B interface (Fig. 1(c), SM Sec. 5). The predicted quantities are the ternary contact angles θAC and θBC. No contact-angle measurement, lens shape, or three-phase simulation data enter the determination of k or of the prefactor; the angles are therefore empirically free outputs of the proposed model. If the common-prefactor assumption for A-C and B-C interfaces failed, the theoretical curves would differ, so the reported agreement is a genuine, falsifiable test of the combination of Young's equation and the assumed scaling. The main weakness—that γAC and γBC are never measured directly and are assumed to obey the same power law with the same prefactor as γAB—is a model-assumption/robustness concern, not a circular reduction: the paper nowhere fits the contact angles or the Young-equation balance to itself. The mean-field γ–λ relation is an explicitly stated analogy adopted from standard textbooks, not a suppressed ansatz imported from the authors' own prior work. No load-bearing self-citation chain is present. Consequently, no circularity step can be exhibited with the required quote-and-reduction evidence.
Assumptions & free parameters
free parameters (3)
- k =
0.23
- c
- tau
assumptions (4)
- domain assumption Time-scale separation: longitudinal relaxation along E is fast, lateral dynamics slow; convective flow ignored.
- domain assumption Pseudo-particle picture: field-aligned clusters act as diffusing pseudo-particles with drift velocity v_alpha,z = q_alpha E / xi.
- ad hoc to paper Equilibrium mean-field theory applies to the nonequilibrium interfaces: tanh profile, lambda = (c/(2-chi))^1/2, gamma = 2 sqrt(2) c phi0^2/(3 lambda).
- ad hoc to paper Universal interfacial-tension scaling: gamma_alpha_beta proportional to (Delta q_alpha_beta E)^0.23 with the same prefactor for AB, AC, and BC.
invented entities (1)
-
Pseudo-particles (field-aligned laning clusters treated as coarse-grained particles)
Cite this review
Pith. "Pith review of Partial-Wetting Phenomena in Active Matter." pith.science (2026). https://pith.science/paper/JHRWF6LI
@misc{pith2026250414818,
author = {Pith},
title = {Pith review of: Partial-Wetting Phenomena in Active Matter},
year = {2026},
howpublished = {\url{https://pith.science/paper/JHRWF6LI}},
note = {Machine review of arXiv:2504.14818}
}
read the original abstract
Abundant interfacial phenomena in nature, such as water droplets on lotus leaves and water transport in plant vessels, originate from partial-wetting phenomena, which can be well described by Young's equation. It remains an intriguing question whether similar behaviors exist in active matter. In this letter, we present a clear demonstration of the partial-wetting phenomenon in a ternary laning system, which is a typical active system. A phase diagram is constructed in which the relative drift velocities of different components govern the transitions among drying, partial wetting, and complete wetting states. The mechanical balance on the contact lines of the partial-wetting phase described by Young's equation is verified. A theoretical picture is proposed to explain the analogy of partial wetting in the laning system to that in the equilibrium system.
Figures
Reference graph
Works this paper leans on
-
[1]
Doi, Soft Matter Physics (oxford university press, 2013)
M. Doi, Soft Matter Physics (oxford university press, 2013)
work page 2013
-
[2]
Safran, Statistical Thermodynamics Of Surfaces, In- terfaces, And Membranes(CRC Press, 2018)
S. Safran, Statistical Thermodynamics Of Surfaces, In- terfaces, And Membranes(CRC Press, 2018)
work page 2018
-
[3]
P. G. De Gennes, Rev. Mod. Phys.57, 827 (1985)
work page 1985
-
[4]
X. Yao, Y. Song, and L. Jiang, Advanced Materials23, 719 (2011)
work page 2011
- [5]
-
[6]
Stenhammar, R
J. Stenhammar, R. Wittkowski, D. Marenduzzo, and M. E. Cates, Phys. Rev. Lett.114, 018301 (2015)
2015
- [7]
- [8]
Show all 37 references
-
[9]
Virnau, Phys
J.T.Siebert, F.Dittrich, F.Schmid, K.Binder, T.Speck, and P. Virnau, Phys. Rev. E98, 030601 (2018)
2018
-
[10]
Chacón, F
E. Chacón, F. Alarcón, J. Ramírez, P. Tarazona, and C. Valeriani, Soft Matter18, 2646 (2022)
2022
-
[11]
J. U. Klamser, S. C. Kapfer, and W. Krauth, Nat Com- mun 9, 5045 (2018)
2018
-
[12]
M. E. Cates and J. Tailleur, Annu. Rev. Condens. Matter Phys. 6, 219 (2015)
2015
- [13]
-
[14]
Das and R
S. Das and R. Chelakkot, Soft Matter16, 7250 (2020)
2020
-
[15]
Rojas-Vega, P
M. Rojas-Vega, P. De Castro, and R. Soto, Phys. Rev. E 107, 014608 (2023)
2023
-
[16]
Tjhung, C
E. Tjhung, C. Nardini, and M. E. Cates, Phys. Rev. X 8, 031080 (2018)
2018
-
[17]
Caballero, A
F. Caballero, A. Maitra, and C. Nardini, Phys. Rev. Lett. 134, 087105 (2025)
2025
-
[18]
Hermann, D
S. Hermann, D. De Las Heras, and M. Schmidt, Phys. Rev. Lett.123, 268002 (2019)
2019
-
[19]
Fausti, E
G. Fausti, E. Tjhung, M. E. Cates, and C. Nardini, Phys. Rev. Lett.127, 068001 (2021)
2021
-
[20]
K. R. Sütterlin, A. Wysocki, A. V. Ivlev, C. Räth, H. M. Thomas, M. Rubin-Zuzic, W. J. Goedheer, V. E. Fortov, A. M. Lipaev, V. I. Molotkov, O. F. Petrov, G. E. Morfill, and H. Löwen, Phys. Rev. Lett.102, 085003 (2009)
2009
-
[21]
Reichhardt, Phys
C.Reichhardt, J.Thibault, S.Papanikolaou,andC.J.O. Reichhardt, Phys. Rev. E98, 022603 (2018)
2018
-
[22]
K. A. Bacik, B. S. Bacik, and T. Rogers, Science379, 923 (2023)
2023
-
[23]
Poncet, O
A. Poncet, O. Bénichou, V. Démery, and G. Oshanin, 6 Phys. Rev. Lett.118, 118002 (2017)
2017
-
[24]
Jose and L
R. Jose and L. Santen, Phys. Rev. Lett. 124, 198103 (2020)
2020
-
[25]
C.-R. Du, K. R. Sütterlin, K. Jiang, C. Räth, A. V. Ivlev, S. Khrapak, M. Schwabe, H. M. Thomas, V. E. Fortov, A. M. Lipaev, V. I. Molotkov, O. F. Petrov, Y. Ma- lentschenko, F. Yurtschichin, Y. Lonchakov, and G. E. Morfill, New J. Phys.14, 073058 (2012)
2012
-
[26]
Piel, Phys
C.Killer, T.Bockwoldt, S.Schütt, M.Himpel, A.Melzer, and A. Piel, Phys. Rev. Lett.116, 115002 (2016)
2016
-
[27]
Mullin, Phys
T. Mullin, Phys. Rev. Lett.84, 4741 (2000)
2000
-
[28]
Dzubiella, G
J. Dzubiella, G. P. Hoffmann, and H. Löwen, Phys. Rev. E 65, 021402 (2002)
2002
-
[29]
G. S. Redner, M. F. Hagan, and A. Baskaran, Phys. Rev. Lett. 110, 055701 (2013)
2013
-
[30]
G. Xu, T. Huang, Y. Han, and Y. Chen, Soft Matter17, 9607 (2021)
2021
-
[31]
See Supplemental Material at [URL will be inserted by publisher] for model and simulation details, theoretical derivations, and additional results on phase separation and partial wetting phenomena
-
[32]
X. Xu, Q. Wei, H. Li, Y. Wang, Y. Chen, and Y. Jiang, Phys. Rev. E99, 043307 (2019)
2019
-
[33]
P. M. Chaikin, T. C. Lubensky, and T. A. Witten,Prin- ciples of condensed matter physics, Vol. 10 (Cambridge university press Cambridge, 1995)
1995
-
[34]
Zhang, E
W. Zhang, E. D. Gomez, and S. T. Milner, Phys. Rev. Lett. 119, 017801 (2017)
2017
-
[35]
Q. Yuan, J. Yang, Y. Sui, and Y.-P. Zhao, Langmuir33, 6464 (2017)
2017
-
[36]
Semprebon, G
C. Semprebon, G. McHale, and H. Kusumaatmaja, Soft Matter 13, 101 (2017)
2017
-
[37]
PARTIAL-WETTING PHENOMENA IN ACTIVE MATTER
Y. Jonosono, S.-i. Tsuda, T. Tokumasu, and H. Na- gashima, Langmuir40, 8440 (2024). 7 SUPPLEMENTAL MATERIAL FOR“PARTIAL-WETTING PHENOMENA IN ACTIVE MATTER” 1.Model and method 1.1.Model To describe the wetting behavior of active particles in a nonequilibrium system, this study ...
2024
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.