Pith. sign in

REVIEW 3 major objections 3 minor 20 references

The Weighted Mean Curvature Derivative of a Space-Filling Diagram

T0 review · 3 major / 3 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper derives an explicit gradient formula for the weighted mean curvature of a space-filling diagram, completing the derivative set needed for the morphometric solvation free energy.

desk verdict The stress-test is right: Lemma 9 and Eq. (63) have a missing factor r_ij, so Theorem 10's s_i is dimensionally wrong as written — but this looks fixable, and the paper otherwise deserves a rigorous referee. read the letter →

arxiv 1908.06779 v2 pith:KDDUARKJ submitted 2019-08-19 cs.CG physics.bio-ph

Arseniy Akopyan , Herbert Edelsbrunner This is my paper · ORCID
classification cs.CGphysics.bio-ph
keywords moleculardynamicsspace-fillingdiagramsintrinsicvolumeweightedmeancurvaturealphashapesinclusion-exclusionderivativesdiscontinuities
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper derives a closed-form formula for the derivative of the weighted mean curvature of a space-filling diagram, the union of balls used to model a molecule, with respect to the positions of the ball centers. That derivative is the missing piece of the morphometric approach to solvation free energy, which writes the nonpolar part of the solvation potential as a linear combination of weighted volume, area, mean curvature, and Gaussian curvature; gradients of the other three terms were already known. The main result, Theorem 10, expresses the derivative at state $x$ with momentum $t$ as $\langle m, t\rangle$, where each centered gradient block $m_i$ is $p_i + q_i + s_i$ and each piece is computed from boundary simplices of the $\alpha$ complex. Theorem 11 shows the gradient is continuous except on a $(3n-1)$-dimensional set $M_{II}$ of critical sphere incidences. If correct, the paper completes the information needed to evaluate the force from the morphometric solvation potential inside a molecular dynamics inner loop.

What carries the argument

The central object is the $\alpha$ complex of the space-filling diagram, the nerve of the Voronoi decomposition, which encodes the fractions of sphere patches, intersection circles, and corner points that belong to the boundary through quantities $\sigma_i$, $\sigma_{ij}$, $\sigma_{ijk}$, and $\nu_{ijk}$. The argument differentiates the weighted mean curvature by decomposing it into a sphere-patch term and an arc term, then splitting each center's motion into a direction-preserving stretch and a distance-preserving rotation. For the arc fractions, Lemma 8 changes the velocity field so that the center and plane of the intersection circle stay fixed, which makes the derivative of $\sigma_{ij}$ tractable; Archimedes' cap-area formula and Heron's formula supply the derivatives of $\sigma_i$, $r_{ij}$, and $\phi_{ij}$. The gradient is obtained by redistributing these local derivatives over the boundary edges and triangles of the $\alpha$ shape.

What would settle it

Evaluate equation (7) at a generic configuration, perturb a single center by $\pm\varepsilon$ along a random direction, and compare the central finite difference of mean with $\langle m, t\rangle$ from Theorem 10; any nonzero gap as $\varepsilon \to 0$ would refute the gradient formula.

Watch

Extended reading notes

Core claim

The paper's central claim is that the weighted mean curvature of a space-filling diagram, as defined in equation (7) with the contribution of each circular arc split equally between the two intersecting spheres, has a derivative that can be written in closed form as $D\mathrm{mean}_x(t) = \langle m, t\rangle$. The gradient $m$ is the sum of three components: $p$, which tracks changes in the exposed-area fractions $\sigma_i$ of the spheres; $q$, which tracks changes in the radii $r_{ij}$ of intersection circles and the dihedral angles $\phi_{ij}$ between outward normals; and $s$, which tracks changes in the fractions $\sigma_{ij}$ of the intersection circles that lie on the boundary. Each component is obtained by differentiating the relevant local quantity and redistributing the resulting derivative over the boundary simplices of the $\alpha$ complex, giving formulas (45), (50), and (64). The paper further proves in Theorem 11 that mean and its gradient are continuous at all states outside $M_{II}$, a $(3n-1)$-dimensional subset characterized by violations of the general-position condition on sphere incidences; across $M_{II}$ the gradient may jump or become infinite according to the cases in Table 1. Together with the previously known gradients of weighted volume, area, and Gaussian curvature, this yields the derivative of the morphometric expression of the solvation free energy.

Load-bearing premise

The formula differentiates the weighted mean curvature defined in equation (7), which splits each arc's contribution equally between the two spheres; if the physical nonpolar solvation energy requires a different split, this gradient is the derivative of a geometric quantity that may not be the target energy.

Editorial extensions

If this is right

  • The morphometric solvation force can be assembled in closed form by combining the new mean-curvature gradient with the known gradients of weighted volume, area, and Gaussian curvature, weighted by the coefficients of the morphometric expansion.
  • Molecular dynamics integrators can treat the nonpolar solvation force as continuous except when the state crosses the exceptional set of critical sphere incidences, where the gradient jumps or diverges according to the scaling laws in Table 1.
  • Because all terms in Theorem 10 are expressed using boundary simplices of the alpha complex, the gradient can be evaluated locally and efficiently, making it suitable for the inner loop of a simulation.
  • The formulas adapt without difficulty to splitting arc mean curvature according to the Voronoi tessellation instead of the equal Apollonius split, as the paper notes in its discussion.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • This reader notes that the equal-split rule in equation (7) is a modeling choice; comparing forces computed from the Apollonius and Voronoi splits on small solvated molecules could show whether the choice matters for the physical free energy.
  • The discontinuity classification suggests that an event-driven integrator could pause at crossings of the exceptional set and use the jump scaling from Table 1 rather than smoothing the force; the paper does not implement this.
  • A natural numerical check, not reported in the paper, is to compare the closed-form gradient against central finite differences of equation (7) on random generic configurations; agreement within rounding error would be strong evidence for correctness.
  • The same decomposition into direction-preserving and distance-preserving motions might extend to derivatives of other weighted curvature functionals defined on the boundary of a union of balls, though the paper does not claim this.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 3 minor

Summary. The paper derives an explicit formula for the gradient of the weighted mean curvature of a union of balls (a space-filling diagram). Using the Voronoi decomposition and the alpha complex, it decomposes the weighted mean curvature into sphere-patch and arc contributions and computes their derivatives under concurrent motion of all centers. The main result, Theorem 10, expresses the gradient as m_i = p_i + q_i + s_i with the three terms given in Eqs. (45), (50), and (64). Theorem 11 states that this gradient is continuous outside M_II, a (3n-1)-dimensional set of non-generic configurations. The paper frames the result as the last missing piece, together with earlier volume, area, and Gaussian curvature derivatives, for computing the morphometric solvation free energy in molecular dynamics.

Significance. If correct, the formula is a useful and nontrivial contribution: it is explicit, parameter-free, and built from geometric quantities computable from the alpha complex, with no fitted data. The derivation is not circular; it follows from the explicit definition of weighted mean curvature in Eq. (7). The continuity analysis in Section 5, if fully proved, would be valuable for robust implementation. A caveat, acknowledged by the authors in Section 6, is that the gradient is for the equal-split definition; if a physically different split of the arc contribution is required, the formula would need to be re-derived. The paper reports no numerical validation, which would be especially helpful in view of the dimensional problem noted below.

major comments (3)
  1. [§3.3, Lemma 9, Eq. (38); §4, Eqs. (63)–(64)] There is a dimensional inconsistency in the radius-change term of Lemma 9 that propagates into Theorem 10. In Eq. (38), the first summand inside the bracket, ⟨T_ijk,x_k−P⟩/⟨x_k−P,u^P_ij⟩, has units length/time, while the second summand, (dα_P/dr_ij)(dr_ij/d‖x_i−x_j‖)⟨u_ij,t_j−t_i⟩, has units 1/time, because dα_P/dr_ij has units 1/length and dr_ij/d‖x_i−x_j‖ is dimensionless. Since the entire bracket is multiplied by 1/(2π r_ij), the second term contributes units 1/(length·time) to σ′_ij, which must be 1/time. The same missing factor appears in Eq. (63): \bar{s}_ijk has units 1/length, so \bar{s}_ijk⟨u_ij,t_j−t_i⟩ in Eq. (60) has units 1/time, whereas s′ in Eq. (41) must have units length/time; correspondingly, the last two terms in Eq. (64) have units 1/length while all other terms in s_i are dimensionless. The likely correction is an additional factor r_ij in the radius-change term of Eq. (38) and in Eq. (63). As printed, the central gradient formula is dimensionally inconsistent and cannot be correct.
  2. [§3.1, Eq. (20) and Lemma 6] The derivation of the distance-preserving contribution to σ′_i is only a sketch. The text asserts that the net area loss or gain can be measured by projecting the relevant portions of the Voronoi edges B_ijk onto the NS axis and then states Eq. (20). This is not a routine chain-rule step: the signed contributions, the role of the fraction ν_ijk, and the conversion from a projection angle to the derivative with respect to ∠x_ix_j need to be spelled out. Since this term feeds directly into p_i in Eq. (45), the proof should be completed or a complete reference supplied.
  3. [§5, Table 1 and Theorem 11] The continuity claims underlying Theorem 11 are asserted rather than proved. In Cases C1–C3 and N01–N23, the statements such as Δarea = ε, Δlength = √ε, Angle = 1, and the corresponding order of Δmean and Δ‖∇mean‖ are introduced with 'by easy analysis' and no derivation. Since Theorem 11 is a stated theorem and the scaling exponents determine which singularities are excluded from the continuity set, a precise local calculation (or a reference containing it) is needed for each case.
minor comments (3)
  1. [Eqs. (30), (38)] Eq. (38) contains a stray closing parenthesis after 1/(2πr_ij), and the square-root symbol in Eq. (30) appears corrupted in the typeset text; the denominator should be the positive square root of the Heron-type expression.
  2. [§6] The equal-split rule for distributing the mean curvature of an arc is a modeling choice. The paper correctly notes that a Voronoi split is also feasible, but it gives no physical or numerical argument that the equal split is the appropriate one for solvation free energy; this should be stated more prominently as a limitation.
  3. [No numerical section] The paper would be strengthened by a small numerical verification of Theorem 10 against finite differences on random configurations; such a check is cheap and would have caught the dimensional inconsistency above.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: Theorem 10 is a direct term-by-term derivative of the explicitly defined weighted mean curvature, with no fitted inputs.

full rationale

The paper's central object, mean(x), is defined explicitly in Eq. (7), restated as Eq. (15), as a sum of sphere-patch terms 4π Σ w_i σ_i r_i minus arc terms (π/2) Σ (w_i+w_j) σ_ij φ_ij r_ij. Theorem 10 is obtained by differentiating this definition term by term: p' from Lemma 6 via Eqs. (39)-(45), q' from Lemma 7 via Eqs. (46)-(50), and s' from Lemma 9 via Eqs. (60)-(64). Every differentiated quantity (σ_i, r_ij, φ_ij, σ_ij) is a geometric function of the ball positions, not a fitted parameter, and no subset of data is used to train or calibrate the formula. The only modeling choice, splitting the arc contribution equally between the two spheres, is made at the definition stage in Section 2, stated openly ('we partition the mean curvature in equal parts to the two intersecting spheres'), and explicitly flagged in Section 6 as an open physical question; changing the split would change the function being differentiated, not make the differentiation circular. The companion reference [1] (same authors) supplies the weighted Gaussian curvature derivative and a corner-partition detail for Eq. (8), but neither supports the derivation of Theorem 10, which is self-contained from Lemmas 6-9. The self-citations [3,7] are published background for the analogous volume and area gradients and are not load-bearing here. No step reduces by construction to its own input; any alleged dimensional error in Lemma 9/Eq. (63) would be a technical correctness issue, not a circularity.

Assumptions & free parameters 0 free parameters · 7 assumptions · 0 invented entities

The paper derives the mean curvature gradient from an explicit definition of weighted mean curvature, using classical geometry (Steiner, Hadwiger, Crofton, Archimedes, Heron) and the morphometric modeling assumption from physics. There are no fitted free parameters and no new entities.

assumptions (7)
  • standard math Steiner's formula and Hadwiger's characterization theorem: every rigid-motion-invariant, continuous, additive valuation on convex bodies is a linear combination of intrinsic volumes.
    Used in Section 1 to justify the morphometric expression Wnp as a linear combination of intrinsic volumes.
  • standard math Crofton's integral formula extends intrinsic volumes to non-convex bodies such as space-filling diagrams.
    Invoked in Section 1 to define volume, area, mean curvature, and Gaussian curvature for the union of balls.
  • standard math Archimedes' theorem on sphere cap area linear in height.
    Used in Section 3.1 (Formula 1) to compute derivatives of exposed sphere fractions.
  • standard math Heron's formula for triangle area.
    Used in Section 3.2 to differentiate the radius of the intersection circle.
  • domain assumption The morphometric approach: the non-polar solvation free energy is a linear combination of weighted volume, area, mean curvature, and Gaussian curvature.
    Adopted from physics literature [12,17] in Section 1; the paper does not derive it, it differentiates this expression.
  • domain assumption The weighted mean curvature is defined via rolling-ball smoothing and the equal split of arc contributions between the two intersecting spheres, as in (7).
    The derivative formula depends on this choice; Section 6 notes an alternative split.
  • domain assumption Generic position assumptions (Conditions I and II in Section 5) and the alpha complex bookkeeping.
    The gradient formula is expressed through fractions computed from the alpha complex, which is well-defined only under generic conditions; the paper characterizes exceptions.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The Weighted Mean Curvature Derivative of a Space-Filling Diagram." pith.science (2026). https://pith.science/paper/KDDUARKJ

@misc{pith2026190806779,
  author       = {Pith},
  title        = {Pith review of: The Weighted Mean Curvature Derivative of a Space-Filling Diagram},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/KDDUARKJ}},
  note         = {Machine review of arXiv:1908.06779}
}
abstract

Representing an atom by a solid sphere in $3$-dimensional Euclidean space, we get the space-filling diagram of a molecule by taking the union. Molecular dynamics simulates its motion subject to bonds and other forces, including the solvation free energy. The morphometric approach [HRC13,RHK06] writes the latter as a linear combination of weighted versions of the volume, area, mean curvature, and Gaussian curvature of the space-filling diagram. We give a formula for the derivative of the weighted mean curvature. Together with the derivatives of the weighted volume in [EdKo03], the weighted area in [BEKL04], and the weighted Gaussian curvature [AkEd19], this yields the derivative of the morphometric expression of the solvation free energy.

Figures

Figures reproduced from arXiv: 1908.06779 by the authors.

Figure 1
Figure 1. After this modification, the point xj moves with speed v along uji, in which v is the length of the new velocity vector. Accordingly, we change the velocity vector of xk to Vijk = tk − ti − ωj × (xk − xi) for every k, and note that Viji = 0 and Vijj = Vij . We finally fix the plane of the circle. To this end, we recall that ξi and ξj given in (17) are the signed distances of xi and xj from the plane of Sij . We have… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

20 extracted references · 20 canonical work pages

  1. [1]

    The weighted Gaussian curvature derivative of a space-filling diagram

    A.\ Akopyan and H.\ Edelsbrunner. The weighted Gaussian curvature derivative of a space-filling diagram. Manuscript, IST Austria, Klosterneuburg, Austria, 2019

    work page 2019

  2. [2]

    The Morse theory of Cech and Delaunay complexes

    U.\ Bauer and H.\ Edelsbrunner. The Morse theory of Cech and Delaunay complexes. Trans.\ Amer.\ Math.\ Soc. 369 (2017), 3741--3762

    work page 2017

  3. [3]

    The area derivative of a space-filling diagram

    R.\ Bryant, H.\ Edelsbrunner, P.\ Koehl and M.\ Levitt. The area derivative of a space-filling diagram. Discrete Comput.\ Geom. 32 (2004), 293--308

    work page 2004

  4. [4]

    M.W.\ Crofton. On the theory of local probability, applied to straight lines drawn at random in a plane; the methods used being also extended to the proof of certain new theorems in the integral calculus. Phil.\ Trans.\ Royal Soc.\ London 158 (1868), 181--199

    work page

  5. [5]

    The union of balls and its dual shape

    H.\ Edelsbrunner. The union of balls and its dual shape. Discrete Comput.\ Geom. 13 (1995), 415--440

    work page 1995

  6. [6]

    Geometry and Topology for Mesh Generation

    H.\ Edelsbrunner. Geometry and Topology for Mesh Generation. Cambridge Univ.\ Press, Cambridge, England, 2001

    work page 2001

  7. [7]

    The weighted-volume derivative of a space-filling diagram

    H.\ Edelsbrunner and P.\ Koehl. The weighted-volume derivative of a space-filling diagram. Proc.\ Natl.\ Acad.\ Sci. 100 (2003), 2203--2208

    work page 2003

  8. [8]

    Three-dimensional alpha shapes

    H.\ Edelsbrunner and E.P.\ M\" u cke. Three-dimensional alpha shapes. ACM Trans.\ Graphics 13 (1994), 43--72

    work page 1994

Show all 20 references
  1. [9]

    Solvation energy in protein folding and binding

    D.\ Eisenberg and A.D.\ McLachlan. Solvation energy in protein folding and binding. Nature (London) 319 (1986), 199--203

    1986

  2. [10]

    u r konvexe K\

    H.\ Hadwiger. Beweis eines Funktionalsatzes f\" u r konvexe K\" o rper. Abh.\ Math.\ Sem.\ Univ.\ Hamburg 17 (1951), 11--23

    1951

  3. [11]

    Solvation of proteins: linking thermodynamics to geometry

    H.\ Hansen-Goos, R.\ Roth, K.R.\ Mecke and S.\ Dietrich. Solvation of proteins: linking thermodynamics to geometry. Phys.\ Rev.\ Lett. 99 (2007), 128101:1--4

    2007

  4. [12]

    A morphometric approach for the accurate solvation thermodynamics of proteins and ligands

    Y.\ Harano, R.\ Roth and S.\ Chiba. A morphometric approach for the accurate solvation thermodynamics of proteins and ligands. J.\ Comput.\ Chem. 34 (2013), 1969--1974

    2013

  5. [13]

    Morphological thermodynamics of fluids: shape dependence of free energies

    P.-M.\ K\" o nig, R.\ Roth and K.R.\ Mecke. Morphological thermodynamics of fluids: shape dependence of free energies. Phys.\ Rev.\ Lett. 93 (2004), 160601:1--4

    2004

  6. [14]

    The interpretation of protein structure: estimation of static accessibility

    B.\ Lee and F.\ Richards. The interpretation of protein structure: estimation of static accessibility. J.\ Mol.\ Biol. 55 (1971), 379--400

    1971

  7. [15]

    Hydrophobicity at small and large length scales

    K.\ Lum, D.\ Chandler and J.D.\ Weeks. Hydrophobicity at small and large length scales. J.\ Phys.\ Chem.\ B 103 (1999), 4570--4577

    1999

  8. [16]

    A morphological model for complex fluids

    K.R.\ Mecke. A morphological model for complex fluids. J.\ Phys.: Condens.\ Matter 8 (1996), 9663--9667

    1996

  9. [17]

    Morphometric approach to the solvation free energy of complex molecules

    R.\ Roth, Y.\ Harano and M.\ Kinoshita. Morphometric approach to the solvation free energy of complex molecules. Phys.\ Rev.\ Lett. 97 (2006), 078101

    2006

  10. [18]

    Electrostatics and dynamics of proteins

    T.\ Simonson. Electrostatics and dynamics of proteins. Rep.\ Prog.\ Phys. 66 (2003), 737--787

    2003

  11. [19]

    Solvation free-energy estimated from macroscopic continuum theory: an accuracy assessment

    T.\ Simonson and A.T.\ Br\" u nger. Solvation free-energy estimated from macroscopic continuum theory: an accuracy assessment. J.\ Phys.\ Chem. 98 (1994), 4683--4694

    1994

  12. [20]

    U ber parallele Fl\

    J.\ Steiner. \" U ber parallele Fl\" a chen. Monatsber.\ Preuss.\ Akad.\ Wiss. (1840), 114--118, [Ges.\ Werke, Vol II (Reimer, Berlin, 1882), 245--308]

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.