REVIEW 2 major objections 5 minor 39 references
This paper provides numerical evidence that floating drops in a capillary tube can have two distinct equilibrium shapes with exactly equal potential energy, making the energy minimizer non-unique, and that in two dimensions an asymmetric wa
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · deepseek-v4-flash
2026-08-03 08:21 UTC pith:4DBX3JS7
load-bearing objection A solid, reproducible numerical study with genuinely new equal-energy observations, but the global non-uniqueness-of-minimizer claim outruns the evidence because only a restricted class of configurations is compared. the 2 major comments →
Comparison of the potential energy for different equilibrium configurations of symmetric and asymmetric floating drops
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For the model of a floating drop in a laterally bounded container, the paper's central numerical finding is that the Euler-Lagrange equations never have unique solutions: for every choice of the nine physical parameters considered, both a centrally located drop and a wall-bound drop exist as equilibria, and in the two-dimensional model also wall-bound and split-wall configurations. Comparing the full potential energy (surface, gravitational, and wetting terms) across these classes, the authors find many parameter values where the centrally located and wall-bound drops have the same energy, so the presumed global minimizer is not unique. They also find a two-dimensional case where the asymmet
What carries the argument
The argument is carried by a high-accuracy numerical solver for the free-boundary Young-Laplace equations, which uses Newton's method with Chebyshev spectral collocation to compute interfaces that meet at the triple junction with the prescribed contact angles, plus nested root-finding loops that match the free-boundary parameter and the drop volume. The energy functional, consisting of surface-tension, gravitational, and wetting terms, is then evaluated on each computed configuration via Clenshaw-Curtis quadrature. This allows a systematic sweep of the nine-dimensional parameter space and direct comparison of energies across configuration classes, which is the mechanism behind every non-uniq
Load-bearing premise
The conclusion that the energy minimizer is sometimes non-unique assumes that the only relevant candidate configurations are the centered, wall-bound, and split-wall drops computed here; if an off-center or fully three-dimensional asymmetric drop with lower energy exists at those parameters, the non-uniqueness claim fails.
What would settle it
At a parameter set where the computed central and wall-bound drops have equal energy, solve for an off-center drop in two dimensions (or a fully three-dimensional asymmetric drop in three dimensions) and evaluate its potential energy: a strictly lower value would disprove the claimed non-unique minimizer. Alternatively, rerun the equal-energy cases with a finer Chebyshev grid or tighter Newton tolerance and check whether the energy difference shrinks to zero or resolves to a nonzero gap.
If this is right
- For open regions of parameter space, the floating-drop variational problem has at least two global minimizers with equal energy, so the ground state is not unique.
- In the two-dimensional model, an asymmetric shape attached to a single wall can be the global minimizer, so symmetry breaking occurs even though the data and domain are symmetric.
- The widely used heuristic that the concavity of the drop-free interface decides between centered and wall-bound minimizers is not reliable; explicit counterexamples exist.
- Along a volume-increase path, the energies of the two symmetric configurations can cross twice, giving two distinct volumes where the minimizer switches or ties.
- The parameter space contains curves along which equal-energy degeneracy persists, so degenerate minimizers are not isolated accidents.
Where Pith is reading between the lines
- Beyond the paper: if the equal-energy sets are codimension-one surfaces in parameter space, then generic small perturbations of parameters will select one configuration, which could make the non-uniqueness hard to observe experimentally unless parameters are tuned.
- Beyond the paper: the two-dimensional asymmetry result suggests that in full three dimensions, off-center or tilted drops might also compete with the symmetric ones; the paper leaves this open, and a numerical method for those PDEs could test it.
- Beyond the paper: the equal-energy crossings resemble exchange-of-stability or pitchfork bifurcations; treating volume or density as a continuation parameter could reveal a connecting bifurcation structure not visible in the present energy profiles.
- Beyond the paper: the claim that off-center drops are unlikely minimizers rests on reflection arguments, but a direct numerical search for such configurations at the equal-energy parameters would be a cheap and decisive check.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a numerical method, based on Chebyshev spectral collocation and Newton iteration, for computing equilibrium configurations of floating drops in laterally bounded containers, under the assumption that configurations are described by generating curves. The method is applied to centrally located and wall-bound drops in R^3, and to central, single-wall, and split-wall drops in R^2. The paper explores a nine-dimensional physical parameter space, computes potential energies of the different configurations, and reports non-uniqueness of solutions to the Euler–Lagrange equations. Its headline claim is that there are many parameter values where a centrally located drop and a wall-bound drop have equal potential energy, giving 'strong evidence of non-uniqueness of energy minimizers' and, in the authors' words, 'no uniqueness at all' for solutions of the floating drop problem.
Significance. If the equal-energy observations are correct, they constitute a striking numerical discovery in a classical capillary free-boundary problem, and the symmetry-breaking examples in R^2 are also valuable. The paper provides a reproducible computational framework (code is on GitHub), explicit quadrature error bounds, and a systematic parameter-space exploration. The main weakness is that the global interpretation of the results rests on an unproven assumption that the computed configuration classes exhaust all relevant minimizers; the authors themselves use hedged language ('presumed', 'likely') but the abstract and conclusions state the stronger claim as established. The numerical evidence for equal-energy crossings also lacks a convergence study.
major comments (2)
- [Abstract; §1; §4.2; §5] The central claim of non-uniqueness of energy minimizers is stated more strongly than the evidence supports. The paper only compares generating-curve families: central and wall-bound drops in R^3, and central, single-wall, and split wall-bound drops in R^2. Non-computed configurations are dismissed heuristically: 'Reflection arguments make the off-center drops unlikely to be energy minimizing' (§4.2) and R^3 asymmetric solutions are 'outside the scope' (§5). Equal energies between the computed families do not establish that the energy minimizer is non-unique in the full problem. Even within the computed families, the paper compares energies of stationary solutions but does not prove they are local or global minimizers under arbitrary perturbations. The hedged phrase '(presumed) energy minimizers' in §5 is appropriate, but the Abstract's 'strong evidence of non-uniqueness of energy minimi
- [§3.1, §4.3, Eqs. (31)–(35)] The equal-energy crossings in Figures 9, 23, and 24 are the quantitative basis for the paper's most striking claim, yet no convergence study is reported for the computed energy values. The paper states 14-digit Newton tolerance and 10-digit BVP tolerance and cites the Clenshaw–Curtis quadrature error bound in Eq. (32), but it does not report actual errors, mesh-refinement checks, or residual norms for the energy computations. Near a crossing, a small systematic error in either energy curve could create or remove an intersection. The authors should provide a convergence test (e.g., n+1 = 14 versus 28 or 56 collocation points) and error estimates for the energy differences at the reported crossings, or at least for one representative case from each of Figures 9, 23, and 24.
minor comments (5)
- [Figure 18 caption] The caption reads 'On the left is a drop adjacent to the left wall, and this is the energy minimizer. On the left is a drop evenly split...' The second 'On the left' should be 'On the right'.
- [§2.2] The text says Bagley and Treinen [2] present 'a united approach'; this should likely be 'a unified approach'.
- [§2.3–2.4] The two fzero loops (for F(ψ̄) and for volume matching) are described with initial guesses and observed sign changes, but no proof of continuity or uniqueness of the zero is given. This is acceptable for a computational paper, but a sentence noting that convergence is empirical and pointing to the code would improve reproducibility.
- [§2.6] The paper says the three surface tensions are 'arbitrarily normalized to add up to a value of 16'; it would help to state explicitly that this is a scaling convention with no loss of generality for the parameter studies reported.
- [§3.2] The heuristic for when a centrally located drop or wall-bound drop is the energy minimizer is stated concisely. The counterexample in §4.1 is a nice addition, but the caption of Figure 13 should say 'tube radius X=2' rather than 'R=2' for consistency with the R^2 setup.
Circularity Check
No significant circularity: the equal-energy and non-uniqueness results emerge from solving boundary-value problems and independent energy quadrature, not from fitting parameters to the target conclusion.
full rationale
The paper's central claims are (i) non-uniqueness of Euler-Lagrange solutions and (ii) common crossings where central and wall-bound configurations have equal potential energy. Both are obtained by solving the underlying boundary-value problems for fixed physical parameters, then evaluating the energy functional (5) using Clenshaw-Curtis quadrature. No parameter is fitted to make the energies equal: the equal-energy points arise as intersections of independently computed energy curves in parameter sweeps (Figures 9, 11, 23, 24). The paper explicitly states it makes no attempt to fit solutions to experimental data (Section 2.6), and the volume-matching and free-boundary matching use fzero on equations (16)-(19) derived from the model, not on energy targets. The self-citations are to a spectral solver [36], existence results [33], and symmetry results [35]; these are tools and supporting theory, not the target non-uniqueness/equal-energy claim, and the solver is accompanied by public code. The restriction to generating-curve configurations and the treatment of off-center and fully 3D asymmetric drops as outside scope is an acknowledged completeness limitation, not a circular derivation. Thus no step reduces by construction to its inputs.
Axiom & Free-Parameter Ledger
free parameters (7)
- Drop volume Vol =
0.4, 0.8, 1, 1.5, 2 in examples
- Tube radius R (R3) / half-width X (R2) =
2
- Drop density ρ1 =
0.01 to 14.9 (ρ0=0, ρ2=15)
- Surface tensions σ01, σ02, σ12 =
triples normalized to sum 16, e.g. (3,7,6), (7.9999,2.0001,6)
- Plate angles γ2_0p, γ1_0p =
e.g. π/2, 2, 0.01; γ2_1p derived from equation (4)
- Chebyshev collocation points n+1 =
14
- Initial-guess factor for \bar r volume loop =
0.2 (adjustable)
axioms (6)
- domain assumption The three-fluid energy functional (5) and Young-Laplace free boundary conditions (6)-(9) from [8] correctly describe floating-drop equilibrium.
- standard math Force balance and the Neumann triangle relation (1)-(2) hold at the triple junction.
- domain assumption Wetting energies satisfy Finn's contact-angle relation (4), giving the third plate angle from two.
- domain assumption Configurations can be reduced to generating curves: radial symmetry in R3 and planar curves in R2.
- standard math Existence and symmetry of the considered configurations follow from [33], [35], [37].
- domain assumption The spectral solver [36] returns numerical solutions within the stated 14-digit/10-digit relative tolerances.
read the original abstract
We provide a numerical method for computing solutions to a free boundary problem arising from the equilibrium state of a floating drop. This numerical method is based on a Newton's method for the underlying nonlinear boundary value problems, and at each iterative step a Chebyshev spectral collocation method is employed. The problems considered here are those that can be described by using generating curves, and include problems in $\mathbb{R}^2$ and $\mathbb{R}^3$. The resulting nine-dimensional space of physical parameters is explored, and examples are given that highlight the potential energy of centrally located drops, wall-bound drops, and asymmetrical configurations in $\mathbb{R}^2$. Non-uniqueness of solutions to the corresponding Euler-Lagrange equations is displayed, and also strong evidence of non-uniqueness of energy minimizers is given.
Figures
Reference graph
Works this paper leans on
-
[1]
Jared L. Aurentz and Lloyd N. Trefethen,Block operators and spectral discretizations, SIAM Rev.59(2017), no. 2, 423–446, DOI 10.1137/16M1065975. MR3646500
-
[2]
Zachary Bagley and Ray Treinen,On the classification and asymptotic behavior of the symmetric capillary surfaces, Exp. Math.27(2018), no. 2, 215–229, DOI 10.1080/10586458.2016.1245641. MR3798195
arXiv 2018
-
[3]
J.C.,An attempt to test the theories of capillary action by comparing the theoretical and measured forms of drops of fluid, Cambridge, 1883
Francis Bashforth and Adams. J.C.,An attempt to test the theories of capillary action by comparing the theoretical and measured forms of drops of fluid, Cambridge, 1883
-
[4]
Ivan Blank, Alan Elcrat, and Raymond Treinen,Geometry of the triple junction between three fluids in equilibrium, Electron. J. Differential Equations (2019), Paper No. 101, 35. MR4012582
2019
-
[5]
Driscoll and Nicholas Hale,Rectangular spectral collocation, IMA J
Tobin A. Driscoll and Nicholas Hale,Rectangular spectral collocation, IMA J. Numer. Anal.36(2016), no. 1, 108–132, DOI 10.1093/imanum/dru062. MR3463435
-
[6]
T. A. Driscoll, N. Hale, and L. N. Trefethen (eds.),Chebfun Guide, Pafnuty Publications, Oxford, 2014
2014
-
[7]
Alan Elcrat, Tae-Eun Kim, and Ray Treinen,Annular capillary surfaces, Arch. Math. (Basel)82(2004), no. 5, 449–467, DOI 10.1007/s00013-003-0101-0. MR2061451
-
[8]
Alan Elcrat, Robert Neel, and David Siegel,Equilibrium configurations for a floating drop, J. Math. Fluid Mech.6 (2004), no. 4, 405–429
2004
-
[9]
Alan Elcrat and Ray Treinen,Numerical results for floating drops, Discrete Contin. Dyn. Syst.suppl.(2005), 241–249. 21
2005
-
[10]
Alan Elcrat and Ray Treinen,Floating drops and functions of bounded variation, Complex Anal. Oper. Theory5 (2011), no. 1, 299–311, DOI 10.1007/s11785-009-0032-2. MR2773068
-
[11]
284, Springer-Verlag, New York, 1986
Robert Finn,Equilibrium capillary surfaces, Grundlehren der Mathematischen Wissenschaften [Fundamental Princi- ples of Mathematical Sciences], vol. 284, Springer-Verlag, New York, 1986
1986
-
[12]
,The contact angle in capillarity, Phys. Fluids18(2006), no. 4, 047102, 7, DOI 10.1063/1.2185655. MR2259293
-
[13]
Gibbs,Ph.D
S.T. Gibbs,Ph.D. Thesis Research Proposal, University of Waterloo (1989)
1989
-
[14]
James Gordon and David Siegel,Properties of annular capillary surfaces with equal contact angles, Pacific J. Math. 247(2010), no. 2, 353–370, DOI 10.2140/pjm.2010.247.353. MR2734153
-
[15]
,Approximating annular capillary surfaces with equal contact angles, Pacific J. Math.247(2010), no. 2, 371–387, DOI 10.2140/pjm.2010.247.371. MR2734154
-
[16]
Jonas Haug and Ray Treinen,Multi-scale spectral methods for bounded radially symmetric capillary surfaces, Electron. Trans. Numer. Anal.60(2024), 20–39, DOI 10.1553/etna vol60s20. MR4695961
doi:10.1553/etna 2024
-
[17]
Henry Ickes and Ray Treinen,The existence of an energy minimizing configuration for multiple solid objects floating in a bath of three liquids, Ann. Mat. Pura Appl. (4)199(2020), no. 2, 821–831, DOI 10.1007/s10231-019-00902-4. MR4079662
-
[18]
Marquis de La Place,Celestial mechanics. Vols. I–IV, Translated from the French, with a commentary, by Nathaniel Bowditch, Chelsea Publishing Co., Inc., Bronx, N.Y., 1966
1966
-
[19]
Math.166(1994), no
Gary Lawlor and Frank Morgan,Paired calibrations applied to soap films, immiscible fluids, and surfaces or networks minimizing other norms, Pacific J. Math.166(1994), no. 1, 55–83
1994
-
[20]
Gian Paolo Leonardi,Infiltrations in immiscible fluids systems, Proc. Roy. Soc. Edinburgh Sect. A131(2001), no. 2, 425–436, DOI 10.1017/S0308210500000937. MR1830418
-
[21]
135, Cambridge University Press, Cambridge, 2012
Francesco Maggi,Sets of finite perimeter and geometric variational problems, Cambridge Studies in Advanced Math- ematics, vol. 135, Cambridge University Press, Cambridge, 2012. An introduction to geometric measure theory. MR2976521
2012
-
[22]
Massari,The parametric problem of capillarity: the case of two and three fluids, Ast´ erisque118(1984), 197–203 (English, with French summary)
U. Massari,The parametric problem of capillarity: the case of two and three fluids, Ast´ erisque118(1984), 197–203 (English, with French summary)
1984
-
[23]
John McCuan,New geometric estimates for Euler elastica, J. Elliptic Parabol. Equ.1(2015), 387–402, DOI 10.1007/BF03377387. MR3487343
-
[24]
,Self-intersection of nod(oid)al curves, Geom. Dedicata216(2022), no. 6, Paper No. 70, 20, DOI 10.1007/s10711-022-00715-5. MR4487570
-
[25]
John McCuan and Ray Treinen,Capillarity and Archimedes’ principle of flotation, Pacific J. Math.265(2013), no. 1, 123–150, DOI 10.2140/pjm.2013.265.123. MR3095116
-
[26]
,On floating equilibria in a laterally finite container, SIAM J. Appl. Math.78(2018), no. 1, 551–570, DOI 10.1137/16M1088818. MR3765926
-
[27]
Frank Morgan,Immiscible fluid clusters inR 2 andR 3, Michigan Math. J.45(1998), no. 3, 441–450, DOI 10.1307/mmj/1030132292
arXiv 1998
-
[28]
A beginner’s guide; Illus- trated by James F
,Geometric measure theory, 5th ed., Elsevier/Academic Press, Amsterdam, 2016. A beginner’s guide; Illus- trated by James F. Bredt. MR3497381
2016
-
[29]
David Siegel,Approximating symmetric capillary surfaces, Pacific J. Math.224(2006), no. 2, 355–365, DOI 10.2140/pjm.2006.224.355. MR2231935
-
[30]
Slobozhanin,Equilibrium and stability of three capillary surfaces with a common line of contact, Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza176(1986), no
L.A. Slobozhanin,Equilibrium and stability of three capillary surfaces with a common line of contact, Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza176(1986), no. 3, 170–173
1986
-
[31]
Trefethen,Spectral methods in MATLAB, Software, Environments, and Tools, vol
Lloyd N. Trefethen,Spectral methods in MATLAB, Software, Environments, and Tools, vol. 10, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2000. MR1776072
2000
-
[32]
MR3012510
,Approximation theory and approximation practice, Society for Industrial and Applied Mathematics (SIAM), Philadelphia, PA, 2013. MR3012510
2013
-
[33]
Ray Treinen,A general existence theorem for symmetric floating drops, Arch. Math. (Basel)94(2010), no. 5, 477–488
2010
-
[34]
,Extended annular capillary surfaces, J. Math. Fluid Mech.14(2012), no. 4, 619–632, DOI 10.1007/s00021- 012-0100-7. MR2992032
doi:10.1007/s00021- 2012
-
[35]
,On the symmetry of solutions to some floating drop problems, SIAM J. Math. Anal.44(2012), no. 6, 3834– 3847, DOI 10.1137/110855569. MR3023432
-
[36]
,Spectral methods for capillary surfaces described by bounded generating curves, Appl. Math. Comput.450 (2023), Paper No. 127886, 17, DOI 10.1016/j.amc.2023.127886. MR4566044 22
arXiv 2023
-
[37]
equilibrium configurations for a floating drop
Raymond Treinen,Discussion of a uniqueness result in “equilibrium configurations for a floating drop”, Electron. J. Differential Equations, posted on 2023, Paper No. 32, 11, DOI 10.58997/ejde.2023.32. MR4574300
-
[38]
Wente,New exotic containers, Pacific J
Henry C. Wente,New exotic containers, Pacific J. Math.224(2006), no. 2, 379–398, DOI 10.2140/pjm.2006.224.379. MR2231937
-
[39]
Brian White,Existence of least-energy configurations of immiscible fluids, J. Geom. Anal.6(1996), no. 1, 151–161. 23
1996
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.