pith. sign in

arxiv: 2501.08247 · v3 · submitted 2025-01-14 · 🧮 math.NA · cs.NA

A Convergent Geometry-Aware Reduction for Diffusion in Branched Tubular Networks

Zachary M. Miksis , Gillian Queisser This is my paper

Pith reviewed 2026-05-23 05:09 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Fick-Jacobs equationdiffusion reductiontubular networksgeometry-aware expansionone-dimensional approximationnumerical stabilityconvergencebranched networks
0
0 comments X

The pith

Treating the Fick-Jacobs derivation as a local Taylor expansion removes geometry-dependent instability and produces the first convergent one-dimensional reduction.

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

The paper shows that instability in the Fick-Jacobs equation arises from a structural inconsistency created by truncation in its original derivation. Recasting the same steps as a locally defined Taylor expansion yields a geometry-aware model whose truncation error stays independent of tube radius variation. A corresponding discretization is then provably stable and converges to the correct reduced solution, while all prior corrections from the literature remain unable to reach that solution no matter how fine the mesh. The same cost-efficient schemes extend directly to branched networks and recover full three-dimensional behavior on neurobiological test cases where the classical reduction fails.

Core claim

The central claim is that a geometry-aware expansion of the Fick-Jacobs model, obtained by treating the derivation as a locally defined Taylor expansion, produces geometry-independent error; its finite-difference discretization is stable and converges to the true reduced solution, whereas standard literature corrections cannot converge to this solution under any spatial refinement.

What carries the argument

The locally defined Taylor expansion of the concentration field, which supplies the geometry-aware correction terms while keeping truncation error independent of tube geometry.

Load-bearing premise

The root cause of instability is the structural inconsistency from truncation in the classical derivation, and recasting the derivation as a local Taylor expansion removes that inconsistency without introducing new uncontrolled errors.

What would settle it

A sequence of successively refined one-dimensional meshes on a fixed test geometry with a known three-dimensional reference solution, showing whether the new method's error decreases while the error from each standard correction stays bounded away from the geometry-aware limit.

Figures

Figures reproduced from arXiv: 2501.08247 by Gillian Queisser, Zachary M. Miksis.

Figure 1
Figure 1. Figure 1: Jacobs flux through a truncated cone. truncated cone and subtracting the total flux through the right face, we get the total flux through the domain. This is equal to the change in concentration in time over the whole volume of the domain. Setting these equal to each other and ignoring infinitesimal ∆x terms, Jacobs derived what is now referred to as the Fick-Jacobs equation, ∂c ∂t = D0  ∂ 2 c ∂x2 + 2 R d… view at source ↗
Figure 2
Figure 2. Figure 2: 2.4.1 Derivation of the Expanded Flux Fick-Jacobs equation Using this reconsideration of the model, the flux from the left side of the truncated cone is −D0 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 2
Figure 2. Figure 2: Expanded Jacobs flux through a truncated cone. [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: A discretized branch point of a network. [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Possible branching 2-paths Pi with concentrations ci . For the first order term, we use second order upwinding. Here we must be careful, as this requires two nodes to the right (or left, depending on the ”wind” direction), and branching can occur at any node in the network. To account for this, we consider all 2-paths, Pi , originating from our origin node [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Numerical approximations by different models with exact solution at [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: L 1 error of different models over time in a truncated cone. 5.1.3 Convergence In [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Convergence errors of different models at [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The (a) 3D rendering of the branched domain, and (b) numerical solution at the large [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Numerical solution on a branched domain at [PITH_FULL_IMAGE:figures/full_fig_p015_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Converging numerical solutions at the left endpoint of a branched domain with different [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Relative L 1 error at different refinement levels at t = 10 for our expanded flux model and the classical Fick-Jacobs model. Additionally, [PITH_FULL_IMAGE:figures/full_fig_p016_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Numerical solutions on a branched domain with lateral influx and outflux at various [PITH_FULL_IMAGE:figures/full_fig_p018_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Converging numerical solutions of a branched domain with lateral flux at different levels [PITH_FULL_IMAGE:figures/full_fig_p018_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Nodes (S, A, B) and nearby nodes (S1, A1, B1). (a) The complete neuronal model used with the radius indicated, and regions of each measuring node highlighted. (b)-(d) Zoom in of each highlighted region, detailing the region local to each measuring node and labels of node locations. calcium influx through post-synaptic density, JSY N ; and the leak calcium flux across the plasma membrane, Jl,p. Details of … view at source ↗
Figure 15
Figure 15. Figure 15: Concentration at each node (S, A, B) and nearby nodes (S [PITH_FULL_IMAGE:figures/full_fig_p020_15.png] view at source ↗
read the original abstract

Diffusion through tubular networks with variable radius arises in a wide range of biological, engineering, and physical applications. The Fick-Jacobs equation is the standard one-dimensional reduction of this problem, briefly derived nearly a century ago in a classical textbook, but was shown to be unstable and inaccurate when the radial gradient is large by Zwanzig in 1992. Three decades of subsequent modifications have failed to resolve this instability because they all inherit a common structural inconsistency introduced by truncation in the original derivation - one that becomes immediately apparent from novel elementary analysis. In this work, we return to the foundations of the Fick-Jacobs derivation and treat it as a locally defined Taylor expansion, recovering a model with geometry-independent error that contrasts directly with the geometry-dependent instability of past corrections. The result is a new geometry-aware expansion of the Fick-Jacobs model, with a numerical discretization that is provably stable and convergent, and the first method known to the authors to converge spatially to the correct geometry-aware solution. Analysis shows that standard corrections from the literature cannot converge to this solution regardless of spatial refinement. We derive efficient numerical schemes for branched networks at equivalent computational cost, and demonstrate that a geometry-aware one-dimensional reduction can faithfully reproduce full three-dimensional results of a neurobiologically relevant problem that the standard reduction cannot achieve.

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

0 major / 3 minor

Summary. The manuscript derives a geometry-aware one-dimensional reduction for diffusion through branched tubular networks of variable radius by recasting the classical Fick-Jacobs approximation as a locally defined Taylor expansion. This produces a model whose truncation error is independent of geometry. The authors supply a discretization that is provably stable and convergent to the correct 3D-averaged solution, prove that standard literature corrections cannot converge to this solution under spatial refinement, develop equivalent-cost schemes for networks, and demonstrate faithful reproduction of full 3D results on a neurobiologically relevant branched geometry where prior reductions fail.

Significance. If the stability, convergence, and non-convergence claims hold, the work resolves a long-standing limitation of reduced-order models for diffusion in complex tubular geometries that has persisted since Zwanzig's 1992 analysis. The first-principles re-derivation, explicit identification of the structural inconsistency in prior truncations, and provision of the first spatially convergent 1D scheme constitute a substantive advance for applications in biology and engineering. The numerical validation on a realistic branched network strengthens the practical relevance.

minor comments (3)
  1. [Abstract] Abstract: the strong claim that 'standard corrections from the literature cannot converge to this solution regardless of spatial refinement' would benefit from an immediate parenthetical pointer to the specific theorem or section containing the supporting analysis.
  2. [Introduction] The introduction would be clearer if the precise assumptions on the scale of radius variation (relative to axial length) were stated explicitly before the derivation begins.
  3. Figure captions for the neurobiological example should include the specific mesh resolution or number of 3D degrees of freedom used in the reference solution for direct comparison.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive and accurate summary of the manuscript, for recognizing its potential significance in resolving long-standing issues with reduced-order models, and for recommending minor revision. No specific major comments were provided in the report.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The derivation chain begins with a standard Taylor expansion applied locally to the diffusion problem, which is an independent first-principles technique with no dependence on fitted parameters, prior self-citations for uniqueness, or renaming of known results. The geometry-aware reduction, stability proof, and convergence claims follow directly from this expansion and subsequent analysis without reducing to the inputs by construction. No load-bearing self-citations or self-definitional steps are present in the provided abstract or described claims.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Review based on abstract only; the derivation appears to rest on standard mathematical assumptions about local expansions and diffusion in cylindrical coordinates with no new free parameters or invented entities described.

axioms (1)
  • domain assumption The radial dependence of the concentration field admits a locally valid Taylor expansion whose truncation error is geometry-independent when properly formulated.
    Central modeling choice that replaces the original truncation.

pith-pipeline@v0.9.0 · 5764 in / 1253 out tokens · 74857 ms · 2026-05-23T05:09:10.590308+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages

  1. [1]

    Alvarez-Ramirez, L

    J. Alvarez-Ramirez, L. Dagdug, L. Inzunza, and E. Rodriguez , Asymmetrical dif- fusion across a porous medium-homogeneous fluid interface , Physica A: Statistical Mechanics and its Applications, 407 (2014), pp. 24–32

    work page 2014

  2. [2]

    Berezhkovskii and G

    A. Berezhkovskii and G. Hummer, Single-file transport of water molecules through a carbon nanotube, Phys. Rev. Lett., 89 (2002), p. 064503

    work page 2002

  3. [3]

    A. M. Berezhkovskii, L. Dagdug, and S. M. Bezrukov, Range of applicability of modified Fick-Jacobs equation in two dimensions , J. Chem. Phys., 143 (2015), p. 164102

    work page 2015

  4. [4]

    A. M. Berezhkovskii, M. A. Pustovoit, and S. M. Bezrukov , Diffusion in a tube of varying cross section: Numerical study of reduction to effective one-dimensional description , J. of Chem. Phys., 126 (2007), p. 134706

    work page 2007

  5. [5]

    Borole, CalciumSim: Simulator for calcium dynamics on neuron graphs using dimension- ally reduced model, Master’s thesis, Temple University, Philadelphia, PA, 2022

    P. Borole, CalciumSim: Simulator for calcium dynamics on neuron graphs using dimension- ally reduced model, Master’s thesis, Temple University, Philadelphia, PA, 2022

    work page 2022

  6. [6]

    Borole, J

    P. Borole, J. M. Rosado, M. R. Neal, and G. Queisser , Neuronal resilience and calcium signaling pathways in the context of synapse loss and calcium leaks: A computational modeling study and implications for Alzheimer’s disease , SIAM J. on Appl. Math., 83 (2023), pp. 2418–2442

    work page 2023

  7. [7]

    Breit and G

    M. Breit and G. Queisser , What is required for neuronal calcium waves? a numerical parameter study, The Journal of Mathematical Neuroscience, 8 (2018), pp. 1–22

    work page 2018

  8. [8]

    Breit, M

    M. Breit, M. Stepniewski, S. Grein, P. Gottmann, L. Reinhardt, and G. Queisser, Anatomically detailed and large-scale simulations studying synapse loss and synchrony using neurobox, Frontiers in Neuroanatomy, 10 (2016), p. 8

    work page 2016

  9. [9]

    Chac ´on-Acosta, A

    G. Chac ´on-Acosta, A. Le ´on-Ram´ırez, and O. Gonz ´alez-Gaxiola, Biharmonic Fick–Jacobs diffusion in narrow channels , Physica A, 628 (2023), p. 129155. 21

    work page 2023

  10. [10]

    Grein, M

    S. Grein, M. Stepniewski, S. Reiter, M. Knodel, and G. Queisser, 1d-3d hybrid mod- eling - from multi-compartment models to full resolution models in space and time , Frontiers in Neuroinformatics, 8 (2014), p. 68

    work page 2014

  11. [11]

    Guan and G

    Q. Guan and G. Queisser, Modeling calcium dynamics in neurons with endoplasmic reticu- lum: Existence, uniqueness and an implicit–explicit finite element scheme , Comm. in Nonlinear Sci. and Num. Sim., 109 (2022), p. 106354

    work page 2022

  12. [12]

    Gurung and Q

    A. Gurung and Q. Guan, Hybrid PDE-deep neural network model for calcium dynamics in neurons, 2024

    work page 2024

  13. [13]

    Hille, Ion Channels of Excitable Membranes, 3rd ed., Sinauer Associates, Inc., Sunderland, Massachusetts, 2001

    B. Hille, Ion Channels of Excitable Membranes, 3rd ed., Sinauer Associates, Inc., Sunderland, Massachusetts, 2001

    work page 2001

  14. [14]

    Hines, Efficient computation of branched nerve equations, Int

    M. Hines, Efficient computation of branched nerve equations, Int. J. Bio-Medical Computing, 15 (1984), pp. 69–76

    work page 1984

  15. [15]

    M. H. Jacobs, Diffusion Processes, Springer-Verlag, New York, 1967

    work page 1967

  16. [16]

    Kalinay, When is the next extending of Fick-Jacobs equation necessary? , J

    P. Kalinay, When is the next extending of Fick-Jacobs equation necessary? , J. Chem. Phys., 139 (2013), p. 054116

    work page 2013

  17. [17]

    Kalinay and J

    P. Kalinay and J. K. Percus , Corrections to the Fick-Jacobs equation , Phys. Rev. E, 74 (2006), p. 041203

    work page 2006

  18. [18]

    U. F. Keyser, B. N. Koeleman, S. van Dorp, R. M. M. S. D. Krapf, S. Lemay, N. Dekker, and C. Dekker , Direct force measurements on DNA in a solid-state nanopore , Nat. Phys., 2 (2006), pp. 473–477

    work page 2006

  19. [19]

    D. P. Lucero , Performance characteristics of permeation tubes , Anal. Chem., 43 (1971), pp. 1744–1749

    work page 1971

  20. [20]

    Z. M. Miksis and G. Queisser , A physics-informed neural network for coupled calcium dynamics in a cable neuron , in ICLR 2024 Workshop on AI4DifferentialEquations in Science, Vienna, AT, May 2024

    work page 2024

  21. [21]

    S. K. Mitusch, S. W. Funke, and M. Kuchta , Hybrid FEM-NN models: Combining artificial neural networks with the finite element method, J. Comp. Phys., 446 (2021), p. 110651

    work page 2021

  22. [22]

    Reguera and R

    D. Reguera and R. Rub´ı, Kinetic equations for diffusion in the presence of entropic barriers, Phys. Rev. E, 64 (2001), p. 061106

    work page 2001

  23. [23]

    J. M. Romero, O. Gonz ´alez-Gaxiola, and G. Chac ´on-Acosta, Exact solutions to Fick-Jacobs equation, Int. J. of Pure and Appl. Math., 82 (2013), pp. 41 – 52

    work page 2013

  24. [24]

    all- or-nothing

    J. Rosado, V. D. Bui, C. A. Haas, J. Beck, G. Queisser, and A. Vlachos , Cal- cium modeling of spine apparatus-containing human dendritic spines demonstrates an “all- or-nothing” communication switch between the spine head and dendrite , PLoS computational biology, 18 (2022), p. e1010069

    work page 2022

  25. [25]

    Seibold, D

    B. Seibold, D. Shirokoff, and D. Zhou , Unconditional stability for multistep ImEx schemes: Practice, J. of Comput. Phys., 376 (2019), pp. 295–391. 22

    work page 2019

  26. [26]

    Y. Shi, Y. T. Lee, and A. S. Kim , Knudsen diffusion through cylindrical tubes of varying radii: Theory and Monte Carlo simulations , Transport in Porous Media, 93 (2012), pp. 517– 541

    work page 2012

  27. [27]

    Shirinpour, N

    S. Shirinpour, N. Hananeia, J. Rosado, H. Tran, C. Galanis, A. Vlachos, P. Jedlicka, G. Queisser, and A. Opitz , Multi-scale modeling toolbox for single neu- ron and subcellular activity under transcranial magnetic stimulation , Brain Stimulation, 14 (2021), pp. 1470–1482

    work page 2021

  28. [28]

    J. W. Thomas , Numerical Partial Differential Equations: Finite Difference Methods , no. 22 in Texts in Applied Mathematics, Springer, 1995

    work page 1995

  29. [29]

    W. M. Yamada, C. Koch, and P. R. Adams , Multiple Channels and Calcium Dynamics , Bradford Books, Toronto, ON, 1998, ch. 2, pp. 137–170

    work page 1998

  30. [30]

    Zwanzig, Diffusion past an entropy barrier , J

    R. Zwanzig, Diffusion past an entropy barrier , J. Phys. Chem., 96 (1992), pp. 3296 – 3930. 23

    work page 1992