REVIEW 2 major objections 4 minor 23 references
Three-phase equilibria in density-functional theory: interfacial tensions
T0 review · 2 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read A two-density mean-field model reduces every non-wet interfacial tension to one Euclidean triangle formula, which also delivers the standard tricritical exponents.
desk verdict Clean geometric tension formula for a toy two-density DFT, with an honest gap: the general-case proof rests on a numerically checked but unproven integral identity. 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 central object is the product free-energy density $F=V_\alpha V_\beta V_\gamma$ with isotropic parabolic wells $V_\nu$ centred at the three bulk phase points. The argument's load-bearing mechanism is the equivalence between the full two-density interface problem and a reduced problem in which $V_\beta$ is a constant set to its midpoint value; this rests on the first integral $\frac{1}{2}(\rho_1')^2+\frac{1}{2}(\rho_2')^2=F$ and on two integral identities, equations (15)–(20), expressing the invariance of $\sigma_{\alpha\gamma}$ under rotations of the $\beta$ point about the $\alpha\gamma$ midpoint. Those identities reduce the variational problem to a straight-line chord integral and yield $\sigma_{\alpha\gamma}=(\sqrt{2}/6)p^3\ell$.
What would settle it
Solve the Euler–Lagrange equations for a non-wet $\alpha\gamma$ interface in a triangle with $\alpha=(1,0)$, $\gamma=(-1,0)$, and $\beta=(\ell\cos\theta,\ell\sin\theta)$ for $\ell=1$, $\theta=\pi/6$, compute $\sigma_{\alpha\gamma}$ numerically to high precision, and compare it with $(\sqrt{2}/6)p^3\ell$; a relative disagreement beyond the numerical error would break equation (22).
Extended reading notes
Core claim
Working with the square-gradient free-energy functional $\Psi = \frac{1}{2}(\rho_1')^2+\frac{1}{2}(\rho_2')^2 + F$, with $F=V_\alpha V_\beta V_\gamma$ and $V_\nu=(\rho_1-\rho^\nu_1)^2+(\rho_2-\rho^\nu_2)^2$, the authors find that the spectator well $V_\beta$ in the interface integral for $\sigma_{\alpha\gamma}$ can be replaced by its value at the midpoint of the $\alpha\gamma$ segment, provided the $\alpha\gamma$ trajectory does not pass through $\beta$. The replacement turns the curved interface trajectory into the straight chord between $\alpha$ and $\gamma$, and the integral evaluates to $(\sqrt{2}/6)p^3\ell$. The authors verify the required rotational-invariance identities numerically for arbitrary triangle shapes, recover earlier conjectures from the isosceles case, and then apply the formula near a tricritical point to obtain $\sigma_{\beta\gamma}\propto\epsilon^{3/2}$ at a critical endpoint and $\mu_t=2$ on approach to the tricritical point.
Load-bearing premise
Everything rests on the numerical claim that rotating the spectator phase point around the midpoint of the other two phase points leaves the interfacial tension exactly unchanged; the paper has no analytic proof of this identity for a general triangle.
Editorial extensions
If this is right
- Interfacial tensions for any non-wet interface in this model can be read off directly from the three bulk density points, eliminating the need to solve the Euler–Lagrange equations for the density profile.
- Near a critical endpoint, the formula gives the mean-field exponent $\mu_c=3/2$ for the critical interfacial tension and a linear approach of the non-critical tension, implying the generic power law $\sigma_{\beta\gamma}\propto(\sigma_{\alpha\beta}-\sigma_{\alpha,\beta\gamma})^{3/2}$.
- On approach to the tricritical point the formula gives the tricritical exponent $\mu_t=2$ for the vanishing interfacial tension, matching mean-field expectations.
- The result sharpens and generalizes earlier conjectures for isosceles three-phase triangles to arbitrary triangle geometry.
Reading between the lines
- If the midpoint-replacement identity persists for products of more than three wells, the same reasoning would predict an $n$-vertex polygon formula for interfacial tensions in $n$-phase equilibria, a testable extension of equation (22).
- The formula suggests that, at least in this class of mean-field models, three-phase interfacial tensions are determined by the metric geometry of the bulk phase points alone, which may make experimental estimates from measured coexistence densities possible without solving profile equations.
- A direct numerical check on a deliberately asymmetric triangle, such as $\theta=\pi/6$ with unequal side lengths, would probe the unproven rotational-invariance identity more stringently than the isosceles and colinear cases.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies a mean-field square-gradient density-functional model for three-phase coexistence, with two density fields and a local free-energy density F = V_alpha V_beta V_gamma, where each V_nu is an isotropic quadratic well centered on the bulk density of phase nu. For a non-wet alpha-gamma interface the authors propose the purely geometric formula sigma_alpha_gamma = (sqrt(2)/6) p^3 ell, with p the distance between the alpha and gamma density points and ell the distance from beta to the midpoint of the alpha-gamma segment. The derivation proceeds through a Hellmann-Feynman-type stationarity relation under bulk-density shifts, an exact collinear calculation, and two integral identities (15)-(20) that are verified numerically but not proved analytically. The formula is then applied near a tricritical point under the scaling rho_2 = -rho_1^2 and field variables s,t; the authors derive the mean-field exponents mu_t = 2 for the approach to the tricritical point, mu_c = 3/2 for the vanishing of the near-critical interfacial tension at a critical endpoint, and the relation sigma_beta_gamma proportional to (sigma_alpha_beta - sigma_alpha,beta_gamma)^{3/2}, all reported to agree with high-precision numerical computation and with earlier special-case conjectures in Ref. [9].
Significance. If Eq. (22) is exact, this is a striking and useful result: in this three-phase DFT the interfacial tension of a non-wet interface is obtained without solving for density profiles, and the formula has a transparent geometric interpretation. The paper also clarifies an apparent discrepancy with model T in Ref. [9] concerning the exponent relating the vanishing critical tension to the noncritical tension difference, and it does so through a derivation with no fitted parameters. The numerical checks reported are consistent with the formula for the geometries tested. The main reservation is that general-triangle exactness rests on unproven integral identities; until those identities are proved, or the claim is appropriately qualified, the central result should be regarded as a well-supported conjecture rather than an exact theorem. The paper is potentially publishable after major revision.
major comments (2)
- [§II, Eqs. (15)–(20)] The exactness of Eq. (22) for a general triangle rests entirely on the rotational-invariance step, which in turn rests on the integral identities (15)–(20). The paper states explicitly that no analytic argument is known for these identities and supports them only by numerical computation. This is load-bearing: if any of these identities fails for some triangle, the replacement of V_beta by its midpoint value is unjustified and Eq. (22) is at best approximate. No numerical method, grid spacing, residual magnitude, or convergence test is reported, so the phrase "hold numerically exactly for all theta" cannot be assessed. The manuscript should either provide a proof of (15)–(20), or state Eq. (22) as a conjecture with detailed quantitative numerical evidence, and correspondingly soften the exactness claims in the abstract and conclusion.
- [§II, around Eq. (22)] The invariance argument moves beta along a full circle about the midpoint of the alpha-gamma segment, but the non-wet character of the alpha-gamma interface is not shown to be preserved along the entire circle. For sufficiently small ell, some positions of beta on that circle would correspond to a wet alpha-gamma interface, in which case the non-wet interfacial tension is not the quantity under discussion and the reduction to the collinear non-wet configuration is invalid. The domain of validity of Eqs. (15)–(22) in the (p, ell, theta) parameter space should be stated explicitly and checked numerically.
minor comments (4)
- [Abstract and Conclusion] The phrase "proven to be numerically exact" is an overstatement; numerical computation cannot prove exactness, and the abstract should distinguish the exact analytic results (collinear case and critical-endpoint limits) from the numerically supported identities underlying the general formula.
- [Eqs. (15)–(20)] The z-dependence of the density profiles is suppressed in the notation V_alpha(rho_1,rho_2), etc.; please state explicitly that these are functionals of the equilibrium profiles and define the polar variables r(z) and phi(z) once in a single place for readability.
- [§III A, Eq. (29)] Eq. (29) is described as a refinement of the conjecture in Eq. (22) of Ref. [9]; the nature of the refinement (exact prefactor versus scaling form) should be stated in one sentence.
- [§III C] The remark that an exceptional circumstance can give an exponent 3/4, postponed to future work, is too cryptic for the present paper; either add a one-sentence indication of the mechanism or delete the remark.
Circularity Check
No significant circularity: the central tension formula is derived from stationarity arguments and checked against independent full-model numerics.
full rationale
The paper's central result, Eq. (22), is not obtained by fitting parameters or by importing a conclusion from prior work. The derivation proceeds from the model free energy (1)-(3), the first integral (5), and the stationarity condition (14). The key step is the rotational-invariance conjecture supported by integral identities (15)-(20), which the paper explicitly states are verified numerically rather than proven analytically. This is an evidentiary gap, not circularity: the numerical verification is independent of the formula it supports. The special collinear case is handled analytically, and the general formula is then checked against high-precision numerical solution of the full two-density model. Applications near the tricritical point use the general formula with scaling variables, and the resulting exponents are not imposed but derived from the analytic expressions, with numerical confirmation. Citations to prior work, especially [9], provide context and special-case conjectures but are not load-bearing: the central derivation does not reduce to those citations. No fit, no parameter, and no conclusion is introduced by definition of the quantity it is said to predict. Therefore no circular step is exhibited.
Assumptions & free parameters
assumptions (4)
- domain assumption The free energy is a square-gradient functional with F(rho1,rho2) = V_alpha V_beta V_gamma, each V_nu an isotropic quadratic well centered at a bulk phase density point.
- standard math Equilibrium density profiles satisfy the first integral (1/2)(rho1')^2 + (1/2)(rho2')^2 = F(rho1,rho2).
- ad hoc to paper The integral identities (15)-(20) hold for all triangle angles theta, expressing rotational invariance of sigma_alpha_gamma under motion of the beta point around the alpha-gamma midpoint.
- domain assumption Near a tricritical point, rho2 = -rho1^2 and the bulk densities are roots of phi(rho1) = rho1^3 - 3t rho1 + 2s.
Cite this review
Pith. "Pith review of Three-phase equilibria in density-functional theory: interfacial tensions." pith.science (2026). https://pith.science/paper/QZSTG523
@misc{pith2026190804721,
author = {Pith},
title = {Pith review of: Three-phase equilibria in density-functional theory: interfacial tensions},
year = {2026},
howpublished = {\url{https://pith.science/paper/QZSTG523}},
note = {Machine review of arXiv:1908.04721}
}
read the original abstract
A mean-field density-functional model for three-phase equilibria in fluids (or other soft condensed matter) with two spatially varying densities is analyzed analytically and numerically. The interfacial tension between any two out of three thermodynamically coexisting phases is found to be captured by a surprisingly simple analytic expression that has a geometric interpretation in the space of the two densities. The analytic expression is based on arguments involving symmetries and invariances. It is supported by numerical computations of high precision and it agrees with earlier conjectures obtained for special cases in the same model. An application is presented to three-phase equilibria in the vicinity of a tricritical point. Using the interfacial tension expression and employing the field variables compatible with tricritical point scaling, the expected mean-field critical exponent is derived for the vanishing of the critical interfacial tension as a function of the deviation of the noncritical interfacial tension from its limiting value, upon approach to a critical endpoint in the phase diagram. The analytic results are again confirmed by numerical computations of high precision.
Figures
Reference graph
Works this paper leans on
-
[9]
R. B. Griffiths, J. Chem. Phys. 60, 195 (1974)
work page 1974
-
[1]
In contrast, in our model the bulk densities are constrained by the physical requirement that they be solutions of (24) for given field values s and t. In closing this subsection we note that in exceptional circumstances within our model an exponent different from the generic value 3/2 may be obtained (notably the value 3/4). We postpone a discussion of thi...
-
[2]
in the (ρ1,ρ 2)-plane represent the densities of three coexisting bulk phases α, β and γ, respectively. Note that F , along with |∂F/∂ρ 1| and|∂F/∂ρ 2|, is zero at any of the densities of the three coexisting phases and otherwise positive so as to describe a state of thermodynamic three-phase equilibrium. The structure of the interface between, say, phase...
-
[3]
be the density point in the ( ρ1,ρ 2)- plane such that V∗ β≡Vβ(ρ∗ 1,ρ∗
-
[4]
is the effective constant value that the spectator well ought to take in order for the correct interfacial tension to be reproduced. With this V∗ β , from (8), we anticipate the following identity: σαγ = 2 ∫ actual dz Vα(ρ1,ρ 2)Vβ(ρ1,ρ 2)Vγ(ρ1,ρ 2) = 2V∗ β ∫ straight dz Vα(ρ1,ρ 2)Vγ(ρ1,ρ 2), (9) where the first integral involves the actual equilibrium densi...
-
[5]
so that this equality holds for an arbitrary location of the β phase point in the (ρ1,ρ 2)-plane. The point (ρ∗ 1,ρ∗
-
[6]
that satisfies (9) can be calculated analytically for the special case when β lies on the straight line through α and γ in the (ρ1,ρ 2)-plane, and outside the open line segment (α,γ ). (If β lay inside, the αγ interface would be necessarily wet by the β phase, and so we excluded the possibility.). Let τ denote the coordinate along this line in the (ρ1,ρ 2)...
-
[7]
satisfies ρ∗ 1 = (ρα 1 +ργ 1)/2 ρ∗ 2 = (ρα 2 +ργ 2)/2 (21) This insight leads to a simple analytic expression for the interfacial tension, applicable to a three-phase triangle of general geometry, and for a non-wet interface, which reads σαγ = √ 2 6 p3𝓁, (22) with p the Euclidean distance from α to γ and 𝓁 the Euclidean distance from β to the midpo...
Show all 23 references
-
[8]
Widom, J
B. Widom, J. Phys. Chem. 77, 2196 (1973)
1973
-
[10]
Lang and B
J.C. Lang and B. Widom, Physica 81A, 190 (1975)
1975
-
[11]
J. C. Lang, P. K. Lim, and B. Widom, J. Phys. Chem. 80, 1719 (1976)
1976
-
[12]
Widom, J
B. Widom, J. Chem. Phys. 67, 872 (1977)
1977
-
[13]
Kerins and B
J. Kerins and B. Widom, J. Chem. Phys. 77, 2016 (1982)
1982
-
[14]
M. M. Telo da Gama and R. Evans, Mol. Phys. 48, 229 (1983)
1983
-
[15]
M. M. Telo da Gama and R. Evans, Mol. Phys. 48, 251 (1983)
1983
-
[16]
Koga and B
K. Koga and B. Widom, J. Phys.: Condens. Matter 28, 244016 (2016)
2016
-
[17]
Rowlinson and B
J.S. Rowlinson and B. Widom, Molecular Theory of Capillarity (Oxford: Clarendon) (1982)
1982
-
[18]
For a review, see e.g., D. Bonn, J. Eggers, J.O. Indekeu, J. Meunier, and E. Rolley, Rev. Mod. Phys. 81, 739 (2009)
2009
-
[19]
Widom, Physica 95A, 1 (1979)
B. Widom, Physica 95A, 1 (1979)
1979
-
[20]
Widom, Chem
B. Widom, Chem. Soc. Rev. 14, 121 (1985)
1985
-
[21]
M. E. Fisher and P.J. Upton, Phys. Rev. Lett. 65, 2402 (1990)
1990
-
[22]
Fisher, Physica A 172, 77 (1991)
M.E. Fisher, Physica A 172, 77 (1991)
1991
-
[23]
Sundar and B
G. Sundar and B. Widom, J. Chem. Phys. 91, 4802 (1987). 14 Figure 1. Equivalence of two ways of obtaining the interfacial tension σαγ of a non-wet interface in a three-phase equilibrium. (a) Three-phase triangle in the plane of two densities ρ1 and ρ2. The bulk phase points ar...
1987
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.