{"id":"15c710b4-39f2-4f9f-a76c-61521e7caaa6","arxiv_id":"1908.04721","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"For the two-density square-gradient model, the tension of a non-wet interface equals (sqrt(2)/6) times the cube of one density-space side times the median distance from the third phase.","lead":"A mean-field model of fluids with two densities is shown to have a simple geometric formula for the tension of interfaces between coexisting phases. The formula is applied near tricritical points to recover known scaling exponents for how interfacial tensions vanish.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (22) rests on unproven integral identities (15)-(20); numerical confirmation for all theta is not a proof, so exactness for arbitrary triangles remains open.","rationale":"The reader's weakest_assumption identifies exactly the same load-bearing concern: the unproven integral identities (15)-(20). The paper's own text concedes the lack of an analytic argument and substitutes numerical confirmation. That is the correct point at which the central claim is least secure. I agree with the conditional verdict: the result is plausible, supported for special cases and by numerical evidence, and the critical-exponent applications are valuable, but the general formula (22) is not fully established analytically. No additional concern of comparable weight emerged. The collinear case is exact, the isosceles case is exact, and the numerical evidence is consistent, so the verdict should remain CONDITIONAL rather than being upgraded or downgraded. The concrete test proposed would settle whether the missing identity holds to numerical precision and, if it fails, would directly invalidate Eq. (22) as an exact statement.","tokens_in":9706,"tokens_out":3710,"duration_ms":40239,"concrete_test":"Solve the full Euler-Lagrange equations for a non-symmetric triangle, e.g. alpha=(1,0), gamma=(-1,0), beta=(cos(pi/4), sin(pi/4)), on a long interval with high-order shooting or relaxation and an adaptive grid. Compute I = integral V_alpha(rho) V_gamma(rho) [rho_1 sin(theta) - rho_2 cos(theta)] dz and the positive normalization N = integral V_alpha V_gamma (|rho_1 sin(theta)| + |rho_2 cos(theta)|) dz. If |I|/N is not below about 1e-10, identity (15), and hence Eq. (22), is not exact. Repeat for theta=0.3 and theta=1.0 to cover non-symmetric cases. In addition, directly test finite-rotation invariance by computing sigma_alpha_gamma at theta and theta+Delta(theta) and verifying the difference is zero to the same relative tolerance; first-order stationarity alone would not guarantee the finite result used in (22).","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is exactness of sigma_alpha_gamma = (sqrt(2)/6) p^3 ell for arbitrary triangle geometry. The derivation has two steps: (i) invariance of sigma_alpha_gamma under rotation of beta about the midpoint of the alpha-gamma segment, established from stationarity condition (14) and the integral identities (15)-(20); and (ii) reduction to the collinear case, where V_beta at the midpoint acts as the effective constant potential. Step (i) is the load-bearing point. The paper states explicitly, 'For the general case we know of no analytic argument for proving the integral identities,' and relies on numerical confirmation. This is a genuine gap: a finite set of computations cannot establish 'numerically exactly for all theta,' and no error tolerance, grid-convergence study, or residual magnitude is reported. If identity (15) or (18) fails for any triangle, rotational invariance fails, the replacement of V_beta by its midpoint value is unjustified, and Eq. (22) is at best approximate. The tricritical-exponent applications in Section III inherit this uncertainty, although the critical-endpoint limits (29) and (31) are separately exact and would be unaffected.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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].","tokens_in":9939,"tokens_out":9960,"duration_ms":97833,"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":[{"comment":"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.","section":"§II, Eqs. (15)–(20)"},{"comment":"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.","section":"§II, around Eq. (22)"}],"minor_comments":[{"comment":"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.","section":"Abstract and Conclusion"},{"comment":"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.","section":"Eqs. (15)–(20)"},{"comment":"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.","section":"§III A, Eq. (29)"},{"comment":"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.","section":"§III C"}],"recommendation":"major_revision","confidential_remarks":"For the editor: the central issue is the unproven integral identities (15)-(20). If the authors can supply a proof, acceptance is straightforward. If not, the manuscript should be revised to present Eq. (22) as a numerically supported conjecture and to provide a proper convergence study. I do not see grounds for rejection: the collinear and critical-endpoint results are exact, and the near-tricritical predictions are internally consistent."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this paper gives you a clean, geometric formula for interfacial tensions in a toy two-density DFT, and the formula is probably right, but the general proof hinges on an unproven integral identity the authors are upfront about.\n\nWhat's actually new: Eq. (22) — sigma_alpha_gamma = (sqrt2/6) p^3 ell — is a genuine generalization of earlier special-case conjectures from Koga-Widom, and it's nice because you can compute the tension without solving the density profile. The application to tricritical scaling is also substantive: using field variables compatible with the tricritical point, they recover the expected mean-field exponent 3/2 for the critical tension as a function of the noncritical tension difference, and explain why the earlier model T gave 3 instead. That correction is credible.\n\nThe derivation is clean in the collinear case: the spectator potential V_beta gets replaced by its midpoint value exactly, using parity. The extension to arbitrary triangles goes through rotational invariance of the tension with respect to moving beta around a circle about the midpoint of alpha-gamma. That invariance rests on integral identities (15)-(20), and the authors say plainly that no analytic argument is known for the general case. They rely on numerical confirmation. That is a real gap. A finite set of high-precision computations can't establish exactness for all theta, and no convergence or residual analysis is reported. If the identity fails for some triangle geometry, formula (22) is approximate rather than exact. The tricritical results in Section III inherit that uncertainty, though the CEP limit (29) is separately exact and the near-CEP expansion (31) follows from a direct straight-line calculation in the limit.\n\nIs this fatal? I don't think so. The model is deliberately simple, the collinear case and the isosceles case are exactly right, and the numerical evidence looks strong. But the paper does claim exactness for arbitrary triangles, and that claim is supported by a conjecture, not a proof. For a physics paper that's borderline acceptable; for anyone relying on (22) away from the numerically checked configurations, it's worth flagging.\n\nThe reliance on Koga-Widom [9] is not circular — the earlier conjectures are special cases and published, and the new exponent correction is argued from physical field constraints. Citation pattern is fine.\n\nWho this is for: people working on wetting and interface theory in mean-field DFT, especially those mapping wetting phase diagrams. It's a useful result in a controlled model, not a broad general theorem.\n\nVerdict: deserves serious peer review. A referee should ask the authors to either prove the identities or soften the exactness claim to a numerically supported conjecture. I'd suggest publishing after that revision.","headline":"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.","tokens_in":10435,"tokens_out":2187,"would_cite":true,"duration_ms":20879,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["82B26","82B27","82D15"],"pacs":[],"model":"deepseek-v4-flash","headline":"A two-density mean-field model reduces every non-wet interfacial tension to one Euclidean triangle formula, which also delivers the standard tricritical exponents.","keywords":["density-functional theory","interfacial tension","three-phase equilibria","tricritical point","critical exponents","square-gradient model","wetting","mean-field theory"],"falsifier":"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).","tokens_in":9495,"feed_emoji":"🔺","tokens_out":6791,"duration_ms":59790,"temperature":0.7,"pith_summary":"This paper establishes that, in a mean-field density-functional model where the local free-energy density is the product of three isotropic parabolic wells in the plane of two densities, the interfacial tension of any non-wet interface between two coexisting phases is given exactly by the simple Euclidean expression $\\sigma_{\\alpha\\gamma}=(\\sqrt{2}/6)p^{3}\\ell$. Here $p$ is the distance between the two phase points in the density plane and $\\ell$ is the distance from the third, spectator phase point to the midpoint of the segment joining them. Because the formula needs only the bulk phase densities and not the density profiles, it turns a variational problem into a piece of triangle geometry. The paper also uses the formula, together with tricritical scaling variables, to derive the mean-field critical exponents for interface tensions near a critical endpoint and near a tricritical point.","feed_headline":"One triangle formula gives every interface tension","feed_subtitle":"For a two-density mean-field fluid, the non-wet interface tension is exactly the triangle quantity (√2/6)p^3ℓ.","key_machinery":"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$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the model T free-energy density and the earlier conjectures for isosceles three-phase triangles that equation (22) generalizes.","marker":"[9]"},{"why":"Provides the Hellmann–Feynman-type variation identity used to derive the rotational invariance conditions (15)–(20).","marker":"[12]"},{"why":"Gives the tricritical-point scaling and the shape of the three-phase region used in Section III.","marker":"[2]"},{"why":"States the triangle inequality for three-phase interfacial tensions that defines the non-wet regime.","marker":"[10]"},{"why":"Prediction of the linear field dependence of the non-critical interfacial tension at a critical endpoint, which equation (33) confirms.","marker":"[13]"},{"why":"Makes the scaling prediction for the singular contribution to the non-critical tension with exponent $\\mu_c$, compared with equation (35).","marker":"[14]"},{"why":"Provides the scaling theory of the critical-endpoint singularities used to interpret the exponent 3/2.","marker":"[15]"}],"fun_headline_variants":["Triangle geometry gives all three interface tensions","Simple triangle formula nails three-phase tensions","Three-phase interfaces from one geometric quantity","Geometric triangle formula fixes interface tensions","All interface tensions from a single triangle expression"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Triangle geometry gives all three interface tensions","Simple triangle formula nails three-phase tensions","Three-phase interfaces from one geometric quantity","Geometric triangle formula fixes interface tensions","All interface tensions from a single triangle expression"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00071,"raw_usage":{"total_tokens":3206,"prompt_tokens":965,"completion_tokens":2241,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":581,"completion_tokens_details":{"reasoning_tokens":2189}},"tokens_in":581,"tokens_out":2241,"duration_ms":16089,"temperature":1.0,"reasoning_tokens":2189,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:33:30.499838+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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).","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the model T free-energy density and the earlier conjectures for isosceles three-phase triangles that equation (22) generalizes."},{"cited_title":"Widom, J","cited_arxiv_id":null,"evidence_quote":"Provides the Hellmann–Feynman-type variation identity used to derive the rotational invariance conditions (15)–(20)."},{"cited_title":"constant of the motion","cited_arxiv_id":null,"evidence_quote":"Gives the tricritical-point scaling and the shape of the three-phase region used in Section III."},{"cited_title":"Lang and B","cited_arxiv_id":null,"evidence_quote":"States the triangle inequality for three-phase interfacial tensions that defines the non-wet regime."},{"cited_title":"Kerins and B","cited_arxiv_id":null,"evidence_quote":"Prediction of the linear field dependence of the non-critical interfacial tension at a critical endpoint, which equation (33) confirms."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Makes the scaling prediction for the singular contribution to the non-critical tension with exponent $\\mu_c$, compared with equation (35)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the scaling theory of the critical-endpoint singularities used to interpret the exponent 3/2."}],"review_version":1}