Abstract
This essay examines the coarse-grained simulation framework of Foley and Johnson, which couples over-damped Brownian dynamics of protein assemblies to Fourier-space propagation of Helfrich elastic membranes, for potential structural analogies to the Riemann Hypothesis via dynamical zeta functions.
Download Full Article
This article is available as a downloadable PDF with complete code listings and syntax highlighting.
Source Paper and Its Mathematics
The source paper by Foley and Johnson presents a hybrid simulation framework for membrane bending by flexible protein lattices. The mathematical core consists of two coupled systems: (1) rigid-body Brownian dynamics for protein triskelia governed by over-damped Langevin equations with stochastic forces, and (2) a Helfrich elastic membrane propagated in Fourier space via Fourier-Space Brownian Dynamics (FSBD). The membrane height field evolves under the Helfrich bending energy functional, quadratic in mean curvature, and the protein-membrane coupling is mediated by harmonic springs. The dynamics are entirely stochastic, living in the diffusive, over-damped regime with no inertial terms.
Attempted Bridge to Riemann
We approached this paper through the suggested angle of dynamical systems and transfer operators, seeking structures analogous to Ruelle zeta functions or Selberg zeta functions. The superficial resonance came from three observations: the use of Fourier-space methods (spectral analysis), the term "dynamical simulation," and the presence of a quadratic energy functional reminiscent of actions in path integrals. One might hope to map the Fourier modes of the membrane to eigenvalues of a Laplacian, the protein assembly configurations to periodic orbits, and the stochastic evolution to a transfer operator whose spectrum encodes zeta zeros.
Assessment: Analogy Strength
The correspondence fails to achieve even the level of formal analogy. The source dynamics are stochastic differential equations (Brownian motion), not deterministic hyperbolic maps or flows. The Fourier decomposition is a spatial spectral method for solving linear PDEs, not the quantum chaos spectrum of the Laplacian on a hyperbolic manifold. There are no periodic orbits to form a dynamical zeta function, no expanding/contracting foliations, and no Ruelle transfer operator. The "spectral gap" in this context refers to the relaxation rate of the Fokker-Planck operator (real, negative eigenvalues), not the complex eigenvalues associated with resonances of chaotic flows.
What Would Be Required
To construct a viable bridge, the source paper would need to study deterministic Anosov flows or expanding maps—for example, geodesic flow on a surface of negative curvature, or billiard dynamics—where the Ruelle zeta function is defined by periodic orbits and converges via hyperbolicity. The current stochastic membrane model, while dynamically rich in the biological sense, lacks the structural ingredients for the dynamical zeta approach to the Riemann Hypothesis.
This essay was produced by an automated research pipeline and has not been peer reviewed; conjectures herein are unproven.