pith. sign in

arxiv: 2601.06863 · v2 · pith:BNNPBT6Enew · submitted 2026-01-11 · 🧮 math.PR · cs.NA· math.NA

Surface Dean--Kawasaki equations

John Bell , Ana Djurdjevac , Nicolas Perkowski This is my paper

Pith reviewed 2026-05-21 16:36 UTC · model grok-4.3

classification 🧮 math.PR cs.NAmath.NA
keywords Dean-Kawasaki equationLangevin dynamicshypersurfacestochastic partial differential equationmartingale formulationfinite-volume discretizationsurface geometry
0
0 comments X

The pith

Stochastic particle density on a hypersurface obeys a Dean-Kawasaki equation whose operators encode the surface metric.

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

The paper starts from microscopic Langevin dynamics of particles constrained to a hypersurface and derives the corresponding surface Dean-Kawasaki equation as a stochastic PDE. The equation is posed in the martingale sense so that the geometry enters directly through the induced metric and the associated differential operators. A reader cares because this supplies a continuum model that automatically respects curvature and can be extended to interacting particles or to a surface that itself evolves according to an SDE coupled to the particles.

Core claim

Starting from the underlying Langevin system, we derive the surface Dean-Kawasaki equation and formulate it in the martingale sense. The resulting SPDE explicitly reflects the geometry of the hypersurface through the induced metric and its differential operators. The framework accommodates pairwise interactions and environmental potentials, and extends to evolving hypersurfaces driven by an SDE that interacts with the particles.

What carries the argument

The surface Dean-Kawasaki equation, the stochastic PDE obtained by passing from the Langevin particle system to the empirical density while keeping the induced metric and its operators explicit in the drift and noise terms.

If this is right

  • The same derivation applies when particles interact through pairwise potentials or feel external environmental potentials.
  • Weak uniqueness holds for the surface DK equation in the non-interacting case.
  • A finite-volume scheme can be constructed that preserves the fluctuation-dissipation balance of the continuous equation.
  • The framework produces a corresponding surface DK equation when the hypersurface itself evolves according to an SDE coupled to the particles.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The geometric form may simplify analysis of curvature-driven segregation or patterning in systems such as proteins on cell membranes.
  • The martingale formulation could be used to obtain large-deviation principles or hydrodynamic limits directly on manifolds.
  • The discretization method offers a practical route to test whether observed surface effects in experiments arise from the metric terms derived here.

Load-bearing premise

The particles obey Langevin dynamics on a hypersurface that admits a Monge gauge parametrization.

What would settle it

Direct numerical simulation of many non-interacting particles obeying the original Langevin dynamics on a fixed surface (for example a paraboloid) and comparison of the observed density fluctuations against solutions of the derived surface DK equation would falsify the claim if the statistics disagree.

Figures

Figures reproduced from arXiv: 2601.06863 by Ana Djurdjevac, John Bell, Nicolas Perkowski.

Figure 1
Figure 1. Figure 1: Mean of ρ. The left panel is the finite volume scheme, the middle panel is a Brownian particles simulation and the right panel is the theoretical result. Note that we have set the range to correspond to the theoretical mean ±10%. Next we consider the variance of the equilibrium distributions computed from the simu￾lations. Although the mean value of ρ is constant over the domain, the variance of ρ is not. … view at source ↗
Figure 2
Figure 2. Figure 2: Mean number of particles in each cell. The left panel is the finite volume scheme, [PITH_FULL_IMAGE:figures/full_fig_p018_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Variance of ρ. The left panel is the finite volume scheme, the middle panel is a Brownian particles simulation and the right panel is the theoretical result. The variance of Ni,j is shown in [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The variance of the number of particles in each cell. The left panel is the finite [PITH_FULL_IMAGE:figures/full_fig_p019_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Time evolution of ρ. The range of ρ in the image is restricted to [0, 0.4] so the values are clipped particularly at early times. The peak values of ρ, left to right, are 0.166, 0.0786, 0.0542 and 0.0457. The height of the peaks on the surface have been scaled by 0.3 so that ρ is not obscured. In [PITH_FULL_IMAGE:figures/full_fig_p020_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Time evolution of Ni,j . The range of Ni,j in the image is restricted to [0, 80] so the values are clipped particularly at early times. The peak values of Ni,j , left to right are 177, 135, 103 and 81. As in [PITH_FULL_IMAGE:figures/full_fig_p021_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Time evolution of ρ with an external potential. The range of ρ in the image is restricted to [0, 0.4] so the values are clipped particularly at early times. The peak values of ρ, left to right are 0.203, 0.125, 0.0957 and 0.0856. As before, the height of the peaks on the surface have been scaled by 0.3 so that ρ is not obscured. Appendix A Sinusoidal perturbation of a square In order to avoid leaving the p… view at source ↗
Figure 8
Figure 8. Figure 8: Time evolution of Ni,j with an external potential. The range of Ni,j in the image is restricted to [0, 80] so the values are clipped particularly at early times. The peak values of Ni,j , left to right are 196, 144, 123 and 109. As before, the height of the peaks on the surface have been scaled by 0.3 so that the solution is not obscured. and G −1 =   1 + q 2 s − pq s − pq s 1 + p 2 s   and G −1/2 … view at source ↗
read the original abstract

We consider stochastic particle dynamics on hypersurfaces represented in Monge gauge parametrization. Starting from the underlying Langevin system, we derive the surface Dean-Kawasaki (DK) equation and formulate it in the martingale sense. The resulting SPDE explicitly reflects the geometry of the hypersurface through the induced metric and its differential operators. Our framework accommodates both pairwise interactions and environmental potentials, and we extend the analysis to evolving hypersurfaces driven by an SDE that interacts with the particles, yielding the corresponding surface DK equation for the coupled surface-particle system. We establish a weak uniqueness result in the non-interacting case, and we develop a finite-volume discretization preserving the fluctuation-dissipation relation. Numerical experiments illustrate equilibrium properties and dynamical behavior influenced by surface geometry and external potentials.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The paper starts from Langevin dynamics for particles on a hypersurface in Monge-gauge parametrization, derives the corresponding surface Dean-Kawasaki SPDE, and formulates it as a martingale problem. The SPDE incorporates the induced metric and surface differential operators. The framework is extended to pairwise interactions, external potentials, and coupled evolving hypersurfaces driven by an interacting SDE. Weak uniqueness is proved for the non-interacting case, a finite-volume scheme that preserves the fluctuation-dissipation relation is constructed, and numerical experiments illustrate equilibrium and dynamical behavior.

Significance. If the derivation and uniqueness result hold, the work supplies a geometrically consistent SPDE description of fluctuating particle systems on curved surfaces, together with a structure-preserving discretization. These elements would be useful for modeling stochastic dynamics on membranes or interfaces where both geometry and noise matter.

major comments (1)
  1. [evolving hypersurfaces section] § on evolving hypersurfaces (the paragraph beginning 'We extend the analysis to evolving hypersurfaces...'): the claim that the derivation proceeds 'without additional geometric constraints' is not yet secured. The Monge-gauge representation assumes the surface remains a graph over a fixed domain; the driving SDE for the surface can produce tilting or folding that violates this graph property, at which point the induced metric and the surface operators used in the SPDE are no longer defined. A concrete mechanism for reparametrization or a restriction on the admissible surface motions is needed to close the argument.
minor comments (2)
  1. [martingale formulation] The statement of the martingale problem (presumably around the definition of the surface DK equation) should explicitly list the test functions and the precise form of the quadratic variation term that encodes the fluctuation-dissipation balance.
  2. [numerical experiments] In the numerical section, the reported equilibrium statistics would be strengthened by a direct comparison against the known invariant measure for the non-interacting case on a flat torus or sphere.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive overall assessment of the manuscript and for the detailed comment on the evolving hypersurfaces section. We address this point below and will revise the manuscript to incorporate the suggested clarification.

read point-by-point responses

  1. Referee: [evolving hypersurfaces section] § on evolving hypersurfaces (the paragraph beginning 'We extend the analysis to evolving hypersurfaces...'): the claim that the derivation proceeds 'without additional geometric constraints' is not yet secured. The Monge-gauge representation assumes the surface remains a graph over a fixed domain; the driving SDE for the surface can produce tilting or folding that violates this graph property, at which point the induced metric and the surface operators used in the SPDE are no longer defined. A concrete mechanism for reparametrization or a restriction on the admissible surface motions is needed to close the argument.

    Authors: We appreciate the referee's careful identification of this subtlety in the Monge-gauge setting for evolving surfaces. The derivation in the paper is carried out under the assumption that the hypersurface admits a global Monge parametrization over a fixed base domain, which by definition requires the surface to remain a graph. We agree that the coupled SDE for the surface height can, in principle, produce evolutions that violate this property through excessive tilting or folding. The current text implicitly relies on local-in-time well-posedness of the surface SDE to stay within the graph regime, but does not spell out explicit safeguards. We will therefore revise the relevant section to state explicitly that the analysis holds under the additional (but standard for Monge gauge) proviso that the surface gradient remains bounded on the time interval of interest, and we will supply sufficient conditions on the drift and diffusion coefficients of the surface SDE that guarantee this bound. A brief remark on local reparametrization to restore the gauge when the bound is approached will also be added. These changes secure the derivation without introducing new geometric constraints beyond those already inherent to the Monge representation. revision: yes

Circularity Check

0 steps flagged

Derivation from external Langevin dynamics is self-contained with no reductions to inputs by construction

full rationale

The paper starts from the standard Langevin particle system on a hypersurface in Monge gauge parametrization and derives the surface Dean-Kawasaki SPDE in the martingale sense, incorporating the induced metric and differential operators. It extends this to interactions, potentials, and coupled evolving surfaces while establishing weak uniqueness in the non-interacting case and a fluctuation-dissipation-preserving discretization. No steps reduce a claimed prediction or first-principles result to a fitted parameter, self-definition, or load-bearing self-citation chain; the central claims remain independent of the outputs they produce.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The derivation rests on the domain assumption that particles follow Langevin dynamics on a Monge-gauge hypersurface; no free parameters or invented entities are indicated in the abstract.

axioms (2)
  • domain assumption Particles obey Langevin dynamics on the hypersurface
    Explicitly stated as the starting point for deriving the surface DK equation.
  • domain assumption Monge gauge parametrization is valid for the hypersurface
    Used to represent the surface geometry in the derivation.

pith-pipeline@v0.9.0 · 5656 in / 1239 out tokens · 38048 ms · 2026-05-21T16:36:13.283687+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Frustrated Fields: Statistical Field Theory for Frustrated Brownian Particles on 2D Manifolds

    cond-mat.stat-mech 2026-05 unverdicted novelty 7.0

    Frustrated Brownian particles on 2D manifolds form a precessing ring whose orientation dynamics are governed by the RP^2 nonlinear sigma model, with the rotational diffusion coefficient and ring profile fixed by simul...

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages · cited by 1 Pith paper

  1. [1]

    Stochastic models of intracellular transport.Reviews of Modern Physics, 85(1):135–196, 2013

    Paul C Bressloff and Jay M Newby. Stochastic models of intracellular transport.Reviews of Modern Physics, 85(1):135–196, 2013. 25

    work page 2013

  2. [2]

    Brownian motion meets riemann curvature.Journal of Statistical Mechanics: Theory and Experiment, 2010(08):P08006, 2010

    Pavel Castro-Villarreal. Brownian motion meets riemann curvature.Journal of Statistical Mechanics: Theory and Experiment, 2010(08):P08006, 2010

    work page 2010

  3. [3]

    Multi-agent system for target tracking on a sphere and its asymptotic behavior.Communications in Nonlinear Science and Numerical Simulation, 117:106967, 2023

    Sun-Ho Choi, Dohyun Kwon, and Hyowon Seo. Multi-agent system for target tracking on a sphere and its asymptotic behavior.Communications in Nonlinear Science and Numerical Simulation, 117:106967, 2023

    work page 2023

  4. [4]

    Giovanni Ciccotti, Tony Lelievre, and Eric Vanden-Eijnden. Projection of diffusions on submanifolds: Application to mean force computation.Communications on Pure and Applied Mathematics: A Journal Issued by the Courant Institute of Mathematical Sciences, 61(3):371–408, 2008

    work page 2008

  5. [5]

    The Dean-Kawasaki equation and the structure of density fluctuations in systems of diffusing particles.Arch

    Federico Cornalba and Julian Fischer. The Dean-Kawasaki equation and the structure of density fluctuations in systems of diffusing particles.Arch. Ration. Mech. Anal., 247(5):Paper No. 76, 59, 2023

    work page 2023

  6. [6]

    Density fluc- tuations in weakly interacting particle systems via the Dean-Kawasaki equation.arXiv preprint arXiv:2303.00429, 2023

    Federico Cornalba, Julian Fischer, Jonas Ingmanns, and Claudia Raithel. Density fluc- tuations in weakly interacting particle systems via the Dean-Kawasaki equation.arXiv preprint arXiv:2303.00429, 2023

    work page arXiv 2023

  7. [7]

    Cambridge University Press, Cambridge, second edition, 2014

    Giuseppe Da Prato and Jerzy Zabczyk.Stochastic equations in infinite dimensions, volume 152 ofEncyclopedia of Mathematics and its Applications. Cambridge University Press, Cambridge, second edition, 2014

    work page 2014

  8. [8]

    Langevin equation for the density of a system of interacting Langevin processes.Journal of Physics A: Mathematical and General, 29(24):L613, 1996

    David S Dean. Langevin equation for the density of a system of interacting Langevin processes.Journal of Physics A: Mathematical and General, 29(24):L613, 1996

    work page 1996

  9. [9]

    Weak error of Dean-Kawasaki equation with smooth mean-field interactions.arXiv preprint arXiv:2502.20929, 2025

    Ana Djurdjevac, Xiaohao Ji, and Nicolas Perkowski. Weak error of Dean-Kawasaki equation with smooth mean-field interactions.arXiv preprint arXiv:2502.20929, 2025

    work page arXiv 2025

  10. [10]

    Weak error analysis for a non- linear SPDE approximation of the Dean-Kawasaki equation.Stoch

    Ana Djurdjevac, Helena Kremp, and Nicolas Perkowski. Weak error analysis for a non- linear SPDE approximation of the Dean-Kawasaki equation.Stoch. Partial Differ. Equ. Anal. Comput., 12(4):2330–2355, 2024

    work page 2024

  11. [11]

    On the accuracy of finite-volume schemes for fluctuating hydrodynamics.Communications in Applied Mathematics and Computational Science, 5(2):149–197, 2010

    Aleksandar Donev, Eric Vanden-Eijnden, Alejandro Garcia, and John Bell. On the accuracy of finite-volume schemes for fluctuating hydrodynamics.Communications in Applied Mathematics and Computational Science, 5(2):149–197, 2010

    work page 2010

  12. [12]

    Duncan, Charles M

    Andrew B. Duncan, Charles M. Elliott, Grigorios A. Pavliotis, and Andrew M. Stuart. A multiscale analysis of diffusions on rapidly varying surfaces.J. Nonlinear Sci., 25(2):389– 449, 2015

    work page 2015

  13. [13]

    Weinan E, Tiejun Li, and Eric Vanden-Eijnden.Applied stochastic analysis, volume 199 ofGrad. Stud. Math.Providence, RI: American Mathematical Society (AMS), 2019

    work page 2019

  14. [14]

    Non-equilibrium large deviations and parabolic- hyperbolic PDE with irregular drift.Invent

    Benjamin Fehrman and Benjamin Gess. Non-equilibrium large deviations and parabolic- hyperbolic PDE with irregular drift.Invent. Math., 234(2):573–636, 2023

    work page 2023

  15. [15]

    Well-posedness of the Dean-Kawasaki and the nonlinear Dawson-Watanabe equation with correlated noise.Arch

    Benjamin Fehrman and Benjamin Gess. Well-posedness of the Dean-Kawasaki and the nonlinear Dawson-Watanabe equation with correlated noise.Arch. Ration. Mech. Anal., 248(2):Paper No. 20, 60, 2024. 26

    work page 2024

  16. [16]

    Consensus- based optimization on the sphere: convergence to global minimizers and machine learn- ing.J

    Massimo Fornasier, Hui Huang, Lorenzo Pareschi, and Philippe S¨ unnen. Consensus- based optimization on the sphere: convergence to global minimizers and machine learn- ing.J. Mach. Learn. Res., 22:Paper No. 237, 55, 2021

    work page 2021

  17. [17]

    Membrane undulations driven by force fluctuations of active proteins.Physical review letters, 93(26):268104, 2004

    Nir Gov. Membrane undulations driven by force fluctuations of active proteins.Physical review letters, 93(26):268104, 2004

    work page 2004

  18. [18]

    From interacting agents to density-based modeling with stochastic PDEs.Commun

    Luzie Helfmann, Nataˇ sa Djurdjevac Conrad, Ana Djurdjevac, Stefanie Winkelmann, and Christof Sch¨ utte. From interacting agents to density-based modeling with stochastic PDEs.Commun. Appl. Math. Comput. Sci., 16(1):1–32, 2021

    work page 2021

  19. [19]

    Hsu.Stochastic analysis on manifolds, volume 38 ofGraduate Studies in Math- ematics

    Elton P. Hsu.Stochastic analysis on manifolds, volume 38 ofGraduate Studies in Math- ematics. American Mathematical Society, Providence, RI, 2002

    work page 2002

  20. [20]

    Stochastic model of slow dynamics in supercooled liquids and dense colloidal suspensions.Physica A: Statistical Mechanics and its Applications, 208(1):35– 64, 1994

    Kyozi Kawasaki. Stochastic model of slow dynamics in supercooled liquids and dense colloidal suspensions.Physica A: Statistical Mechanics and its Applications, 208(1):35– 64, 1994

    work page 1994

  21. [21]

    On Dean–Kawasaki dy- namics with smooth drift potential.J

    Vitalii Konarovskyi, Tobias Lehmann, and Max von Renesse. On Dean–Kawasaki dy- namics with smooth drift potential.J. Stat. Phys., pages 1–16, 2020

    work page 2020

  22. [22]

    von Renesse

    Vitalii Konarovskyi, Tobias Lehmann, and Max-K. von Renesse. Dean-Kawasaki dy- namics: ill-posedness vs. triviality.Electron. Commun. Probab., 24:Paper No. 8, 9, 2019

    work page 2019

  23. [23]

    A math- ematical perspective

    Tony Leli` evre, Mathias Rousset, and Gabriel Stoltz.Free energy computations. A math- ematical perspective. London: Imperial College Press, 2010

    work page 2010

  24. [24]

    Nonequilibrium membrane fluctua- tions driven by active proteins.The Journal of chemical physics, 124(7), 2006

    Lawrence C-L Lin, Nir Gov, and Frank LH Brown. Nonequilibrium membrane fluctua- tions driven by active proteins.The Journal of chemical physics, 124(7), 2006

    work page 2006

  25. [25]

    Membrane proteins: molecular dynamics simula- tions.Current opinion in structural biology, 18(4):425–431, 2008

    Erik Lindahl and Mark SP Sansom. Membrane proteins: molecular dynamics simula- tions.Current opinion in structural biology, 18(4):425–431, 2008

    work page 2008

  26. [26]

    from cells to vesicles/II

    Reinhard Lipowsky and Erich Sackmann.Structure and dynamics of membranes: I. from cells to vesicles/II. generic and specific interactions. Elsevier, 1995

    work page 1995

  27. [27]

    Diffusion on ruffled membrane surfaces.The Journal of chemical physics, 126(23), 2007

    Ali Naji and Frank LH Brown. Diffusion on ruffled membrane surfaces.The Journal of chemical physics, 126(23), 2007

    work page 2007

  28. [28]

    Lateral diffusion of a protein on a fluctuating membrane

    Ellen Reister and Udo Seifert. Lateral diffusion of a protein on a fluctuating membrane. Europhysics Letters, 71(5):859, 2005

    work page 2005

  29. [29]

    Hybrid simulations of lateral diffusion in fluctuating membranes.Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 75(1):011908, 2007

    Ellen Reister-Gottfried, Stefan M Leitenberger, and Udo Seifert. Hybrid simulations of lateral diffusion in fluctuating membranes.Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 75(1):011908, 2007

    work page 2007

  30. [30]

    Diffusing proteins on a fluctuating membrane: Analytical theory and simulations.Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 81(3):031903, 2010

    Ellen Reister-Gottfried, Stefan M Leitenberger, and Udo Seifert. Diffusing proteins on a fluctuating membrane: Analytical theory and simulations.Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 81(3):031903, 2010

    work page 2010

  31. [31]

    Single-particle tracking: applications to membrane dynamics.Annual review of biophysics and biomolecular structure, 26(1):373–399, 1997

    Michael J Saxton and Ken Jacobson. Single-particle tracking: applications to membrane dynamics.Annual review of biophysics and biomolecular structure, 26(1):373–399, 1997. 27

    work page 1997

  32. [32]

    Veretennikov

    Alexander Yu. Veretennikov. On strong solution and explicit formulas for solutions of stochastic integral equations.Math. USSR Sb., 39:387–403, 1981. 28

    work page 1981