Pith. sign in

REVIEW 3 major objections 4 minor 25 references

Hydrodynamic Limit of the Symmetric Zero-Range Process with Slow Boundary

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

Pith's one-line read The zero-range process with slow boundary reservoirs has a rigorous hydrodynamic limit: its empirical density converges to the unique weak solution of a nonlinear heat equation, with boundary conditions that switch from Robin (at the critic

desk verdict Rigorous hydrodynamic limit for a standard boundary-driven zero-range model, but the proof leans on an unproved boundary replacement lemma that carries the θ-dependent boundary conditions. read the letter →

arxiv 2508.19447 v1 pith:6IWMKNDW submitted 2025-08-26 math.PR

classification math.PR MSC 60K3582C2235K55
keywords hydrodynamiclimitzero-rangeprocessslowboundaryRobinconditionnonlinearheatequationentropymethodreplacementlemmaspectralgap
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 proves that the large-scale density of a symmetric zero-range process on a finite interval, in contact with boundary reservoirs that inject and remove particles at rates slowed by a factor N^{-θ} with θ ≥ 1, converges in probability to the unique weak solution of a nonlinear heat equation. The limiting equation is ∂_t ρ = ΔΦ(ρ), where Φ is the flux function determined by the jump rate g, and the boundary conditions depend sharply on θ: for θ = 1 they are Robin-type, while for θ > 1 they degenerate to zero-flux Neumann conditions as the reservoirs become asymptotically invisible. This is the first rigorous derivation in this boundary-driven setting, closing a gap left open by an earlier heuristic result. The proof follows the entropy method but requires substantial new ingredients because the invariant measure with reservoirs is not translation-invariant; the authors replace it by the invariant measure of the dynamics restricted to the bulk and control the discrepancy via a spectral gap estimate.

What carries the argument

The central objects are the flux function Φ = R^{-1}, where R(φ) is the expected occupation under the product invariant measure with fugacity φ (so that Φ(α) is the expected jump rate at density α), and the entropy method with a translation-invariant reference measure: the invariant measure of the dynamics restricted to the bulk, denoted ¯ν_φ. The one-block estimate uses a spectral gap lower bound (Assumption (SG)) together with the Rayleigh estimate, while the two-block estimate requires a new spectral gap bound for a coupled two-box zero-range process. The indicator κ̃ = κ1_{θ=1} encodes the critical slowdown: only at θ = 1 do the reservoirs survive in the macroscopic boundary conditions.

What would settle it

Run microscopic simulations for a specific jump rate that satisfies (ND) and (2.1) but violates (SG) — for example, a rate with strong blocking that makes the spectral gap decay faster than ℓ^{-2} — and test whether the empirical density still converges to the nonlinear heat equation with the stated boundary conditions; a persistent mismatch as N grows would refute the sufficiency of the theorem's hypotheses, while convergence would indicate that (SG) is not necessary.

Watch

Extended reading notes

Core claim

Theorem 2.12: for θ ≥ 1, the sequence of empirical measures π_t^N converges in probability to the absolutely continuous measure ρ(t,u)du, where ρ is the unique weak solution of the nonlinear heat equation ∂_t ρ = ΔΦ(ρ) with boundary conditions ∂_uΦ(ρ(0)) = κ̃(λΦ(ρ(0)) − α) and ∂_uΦ(ρ(1)) = κ̃(β − δΦ(ρ(1))), with κ̃ = κ1_{θ=1}. The proof is built on the entropy method of Guo, Papanicolaou and Varadhan, adapted to the non-conservative, non-translation-invariant setting by using the bulk-restricted invariant measure as a reference, establishing one-block and two-block replacement lemmas via spectral gap estimates, and proving a spectral gap bound for a coupled two-box process. The restriction θ

Load-bearing premise

Assumption (SG): the spectral gap of the zero-range process on a box of size ℓ with j particles must be bounded below by C/ℓ², uniformly in j; if this rate fails, the denominators in the Rayleigh estimates inside the replacement lemmas can lose positivity and the proof collapses.

Editorial extensions

If this is right

  • The hydrostatic limit follows as a corollary: the stationary measure ¯ν_N converges to the stationary profile ¯ρ_θ, the solution of the elliptic problem with the corresponding Robin or Neumann boundary conditions.
  • For θ > 1 the reservoirs have no macroscopic influence: the limiting equation is the nonlinear heat equation with zero-flux Neumann boundary conditions, so the system behaves as isolated in the limit.
  • The uniqueness of the weak solution, proved in Lemma 2.10, ensures that the identification of the limit is complete: any limit point must equal the unique solution.
  • The replacement lemmas are stated for the jump rate g but extend to any Lipschitz cylinder function Ψ, yielding a Boltzmann-Gibbs-type principle for a whole class of observables, not just the jump rate.
  • The spectral gap assumption (SG) is only used in the replacement lemmas; the paper notes that a positive lower bound depending on g, ℓ and j would suffice.

Reading between the lines

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

  • The method likely extends to long-range jumps, where the same entropy approach with a translation-invariant reference measure and replacement lemmas would produce convergence in probability to a fractional nonlinear heat equation—this is suggested by the paper's discussion of the exclusion analogue, but is not proved here.
  • For θ < 1 the paper's strategy would require the stationary measure as reference, at the cost of losing translation invariance; a plausible outcome is that Dirichlet boundary conditions emerge, but the necessary equivalence-of-ensembles and spectral gap estimates are not yet available.
  • The stochastic domination hypothesis µ_N ≤ ¯ν_N may be replaceable by a weaker uniform moment bound combined with a truncation argument, since its only role is to bound high-density configurations and moments of g; a testable variation is to verify the hydrodynamic limit for product initial measures with slowly varying parameter without imposing domination.
  • One can quantitatively test the critical slowdown: for θ = 1, the boundary flux should equal λΦ(ρ(0)) − α in the limit, while for θ > 1 it should vanish; finite-size simulations measuring the boundary flux for increasing N would directly probe the sharp transition.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper studies the symmetric nearest-neighbour zero-range process on the finite interval {1,...,N-1} in contact with reservoirs at the boundary sites 1 and N-1, with injection/removal rates scaled by N^{-theta}, theta >= 1. Under assumptions that the jump rate g is nondecreasing with bounded increments, the initial measures have relative entropy O(N) with respect to the NESS and are stochastically dominated by it, and that a spectral-gap lower bound (SG) holds, the authors claim that the empirical density converges in probability to the unique weak solution of the nonlinear heat equation partial_t rho = Delta Phi(rho) with Robin boundary conditions for theta=1 and Neumann conditions for theta>1 (Eqs. (2.23)-(2.24)). The proof follows the entropy method: Dynkin martingales, tightness, absolute continuity of limit points, bulk one-block and two-block replacement lemmas, a boundary replacement lemma, and an energy estimate for the regularity of Phi(rho).

Significance. If the proof is completed, this would be the first rigorous hydrodynamic limit for the boundary-driven zero-range process in the nearest-neighbour case with slow reservoirs, closing a gap left open in [FMN21] and extending the slow-boundary exclusion results to a model with unbounded occupation numbers. The theta=1 Robin-type transition is a genuinely new feature for zero-range. The paper has clear strengths: the explicit product NESS is used intelligently, the translation-invariant reference measure is a good choice, Lemma 2.7 on the coupled two-box spectral gap is a useful tool, and the uniqueness proof for the nonlinear Robin problem is careful. The bulk replacement estimates are detailed. However, the boundary treatment is incomplete in a load-bearing way: the boundary replacement lemma is stated without proof, and the remainder R5 in Section 3.4.1 is asserted to vanish by an argument that does not apply to it. The energy estimate is also delegated to [KL99] with little adaptation. These gaps must be addressed before the main theorem can be considered established.

major comments (3)
  1. [Section 3.4.1, Eq. (3.27), R5] The proof asserts that E_mu[|R5|] tends to 0 'by the same argument given for R3'. This is not justified and is load-bearing. R3 is a bulk remainder with an explicit 1/N factor and a smooth factor Delta G_s; after summation by parts it is controlled by the continuity of Delta G. R5 has no such factor: it contains g(eta_s(1)) - (1/epsilon N) sum_{y=1}^{epsilon N} g(eta_s(y)) (and the analogous right-boundary term), multiplied by partial_u G_s(0), plus for theta=1 additional terms with G_s(0). A single-site rate does not converge in L1 to its spatial block average under product-type measures with slowly varying profile; the L1 difference is of order one. Vanishing of R5 is a boundary one-block replacement, not a corollary of the R3 estimate. Without a proof of this step, the boundary terms in (3.27) cannot be replaced by Phi of directed density averages, and the boundary conditions in (2.24
  2. [Section 3.6, Lemma 3.14] The boundary replacement lemma is stated without proof, with only the remark that the proof is 'almost identical' to the bulk argument. This lemma is exactly what sends R6 to zero in Section 3.4.1 and hence produces the Robin (theta=1) and flux (theta>1) boundary terms in (2.24). The directed averages in (3.28) are one-sided, the boxes touch the boundary, and the reference measure is not translation-invariant near the boundary; the bulk one-block and two-block estimates (Lemmas 3.12-3.13) are formulated for symmetric bulk boxes. Surface terms in the Feynman-Kac/Rayleigh estimates (3.43)-(3.48) could survive after summation. A complete proof, or a precise reduction to Lemmas 3.12-3.13, is required.
  3. [Appendix A.1, Lemma A.1 / Eq. (A.3)] The proof of the energy estimate stops at (A.3) with 'the exact same argument given in the proof of [KL99, Lemma 7.3]'. This estimate is used to prove Proposition 3.8, i.e., item ii) of Definition 2.9, namely Phi(rho) in L^2([0,T],H^1). The adaptation is not routine: the process has boundary reservoirs, the Dirichlet form D_N^0 in (3.44) is only the bulk part, and the auxiliary Lemmas B.1-B.3 involve boundary corrections. The reader cannot verify that the expression W_N is controlled without seeing the argument. Please include the details or a precise statement of the modified estimate.
minor comments (4)
  1. [Appendix A.1, first paragraph] The text says 'By the bulk replacement lemma, namely Lemma 3.14'; Lemma 3.14 is the boundary replacement lemma. The reference should be to Lemma 3.9.
  2. [Eq. (3.27)] In the last term there is a misplaced parenthesis: 'delta { g(eta_s(N-1) - 1/(epsilon N) ...' should read 'delta ( g(eta_s(N-1)) - 1/(epsilon N) ...'.
  3. [Lemma 3.13] The statement uses 'sup_{|y|<=epsilon N}' inside the expectation while y subsequently appears in the sums; clarify whether y is a fixed index over which the sup is taken. Also the order of limits in the lemma statement should match the order used at the end of the proof (N -> infinity, then epsilon -> 0, then ell -> infinity).
  4. [Remark 3.2] The notation 'Q_N' with 'N >= 0' appears to be a typo for 'N >= 1'.

Circularity Check

0 steps flagged · score 2.0 of 10

No circular derivation: the hydrodynamic equation and its boundary conditions are obtained from the generator via martingale and replacement arguments; the omitted boundary replacement proof is a rigor gap, not a circular step.

full rationale

I found no circular step in which a prediction is equivalent to an input by construction. The main result, Theorem 2.12, is not obtained by fitting parameters to data or by renaming a known pattern: the limiting PDE (2.23) and its θ-dependent Robin boundary terms are derived from the Dynkin martingale decomposition of the empirical measure, after bulk and boundary replacement lemmas. The boundary conditions are read off from the generator terms after passing to limits, not imposed by hand. Two passages assert missing support and should be weighed as rigor gaps rather than circularity. First, Lemma 3.14 (Boundary Replacement Lemma) is stated and then its proof is omitted: 'The proof of this lemma follows from an argument which is almost identical to the one given in the previous section, so we omit it. The reader is invited to retrace all the steps of the proof of Lemma 3.9, replacing the centred averages with the directed averages, and will see that no additional argument is required.' This lemma directly controls the remainder R6_{N,ε} in Section 3.4.1, and hence the θ=1 Robin boundary terms. The omission is a completeness gap in the proof of the boundary condition, but it is not circular: Lemma 3.14 is not a restatement of Theorem 2.12, and the analogous bulk lemma is proved in detail in Section 3.5. Second, the energy estimate in Appendix A.1 is delegated to [KL99]: 'From here, the conclusion follows from the exact same argument given in the proof of [KL99, Lemma 7.3].' This is another proof gap, but again not a circular reduction. The paper does lean on several results coauthored by the present authors: [FMN21] (with A. Neumann), [BDGN20] (with P. Gonçalves and A. Neumann), and [BMNS17] (with A. Neumann). These are used for auxiliary ingredients: a Feynman-Kac formula, spectral gap examples, and a uniqueness framework. These citations are not load-bearing in the sense that the central convergence theorem reduces to them; the paper provides proofs or adaptations for the claims it directly uses, and the cited items are standard tools with independent content. There is no invocation of a self-authored uniqueness theorem to forbid alternatives; uniqueness of weak solutions is proved in Appendix A.2 using a generalization of [BDGN20, Lemma 7.3] with proof included. Assumption (SG), the spectral gap lower bound, is an explicit hypothesis of the theorem rather than a derived consequence, and the examples satisfying it are cited from the external literature

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

No new physical entities are postulated. The 'coupled' two-box zero-range process in Lemma 2.7 is a proof device for the spectral gap bound, not a new physical object. The only auxiliary free choice is the reference fugacity phi, which does not enter the final theorem statement.

free parameters (1)
  • reference fugacity phi = phi >= max{sup_{N,x} phi_bar_N(x), alpha/lambda, beta/delta}
    Introduced in the replacement lemmas (Section 3.5) to define a translation-invariant reference measure nu_bar_phi; chosen by hand to dominate the boundary rates, not fitted to data.
assumptions (5)
  • domain assumption Assumption (ND): g is nondecreasing (2.16)
    Used to prove attractiveness (Lemma 2.6), essential for transferring stochastic domination and controlling unbounded occupations.
  • domain assumption Assumption (SG): spectral gap gap(ell,j) >= C/ell^2 for box size ell with j particles (2.18)
    Load-bearing for one-block and two-block estimates (Lemmas 3.12, 3.13) via the Rayleigh estimate; examples cited from [LSV96], [Mor06], [Nag10].
  • domain assumption lim_{phi uparrow phi*} Z(phi) = infinity (after Remark 2.2)
    Ensures the range of R is [0, infinity) and Phi = R^{-1} is well defined; needed to write the hydrodynamic PDE.
  • domain assumption Condition (2.6) for the unique invariant product measure nu_bar_N
    Taken from [LMS05, Equation 20]; guarantees the NESS exists with explicit fugacity profile (2.7).
  • domain assumption Initial measures satisfy H(mu_N | nu_bar_N) <~ N and mu_N <= nu_bar_N (Theorem 2.12 hypotheses)
    Entropy bound permits the change of measure to the reference measure; stochastic domination controls the unbounded particle number via attractiveness.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Hydrodynamic Limit of the Symmetric Zero-Range Process with Slow Boundary." pith.science (2026). https://pith.science/paper/6IWMKNDW

@misc{pith2026250819447,
  author       = {Pith},
  title        = {Pith review of: Hydrodynamic Limit of the Symmetric Zero-Range Process with Slow Boundary},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/6IWMKNDW}},
  note         = {Machine review of arXiv:2508.19447}
}
abstract

We study the hydrodynamic behaviour of the symmetric zero-range process on the finite interval $\{1, \ldots, N-1\}$ in contact with slow reservoirs at the boundary. Particles are injected and removed at sites $1$ and $N-1$ at rates that scale like $N^{-\theta}$ with $\theta\ge1$. Under mild assumptions on the jump rate and the sequence of initial measures, we show that the empirical density evolves on the diffusive scale according to a nonlinear heat equation, with boundary conditions reflecting the strength of the reservoirs.

Figures

Figures reproduced from arXiv: 2508.19447 by the authors.

Figure 1
Figure 1. Microscopic dynamics of the symmetric zero-range process with open boundaries. [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

25 extracted references · 17 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type url doi eprint volume year label extra.label sort.label INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION not #0 #2 if INTEGERS l FUNCTION string.length #1 'l := ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...

  3. [3]

    D. Aldous . Stopping Times and Tightness . Ann. Probab. 6, no. 2, (1978), 335--340. doi:10.1214/aop/1176995579 https://doi.org/10.1214/aop/1176995579

  4. [4]

    Bernardin , P

    C. Bernardin , P. Cardoso , P. Gonçalves , and S. Scotta . Hydrodynamic Limit for a Boundary Driven Super-Diffusive Symmetric Exclusion . Stochastic Process. Appl. 165, (2023), 43--95. doi:10.1016/j.spa.2023.08.002 https://doi.org/10.1016/j.spa.2023.08.002

  5. [5]

    Bonorino , R

    L. Bonorino , R. De Paula , P. Gonçalves , and A. Neumann . Hydrodynamics of Porous Medium Model with Slow Reservoirs . J. Stat. Phys. 179, no. 3, (2020), 748--788. doi:10.1007/s10955-020-02550-y https://doi.org/10.1007/s10955-020-02550-y

  6. [6]

    Bernardin , P

    C. Bernardin , P. Gonçalves , and B. Jim\'enez-Oviedo . Slow to Fast Infinitely Extended Reservoirs for the Symmetric Exclusion Process with Long Jumps . Markov Process. Related Fields 25, no. 2, (2019), 217--274

  7. [7]

    Bernardin , P

    C. Bernardin , P. Gonçalves , and B. Jim\'enez-Oviedo . A Microscopic Model for a One Parameter Class of Fractional Laplacians with Dirichlet Boundary Conditions . Arch. Ration. Mech. Anal. 239, no. 1, (2021), 1--48. doi:10.1007/s00205-020-01549-9 https://doi.org/10.1007/s00205-020-01549-9

  8. [8]

    Bernardin , P

    C. Bernardin , P. Gonçalves , B. Jim\'enez-Oviedo , and S. Scotta . Non-Equilibrium Stationary Properties of the Boundary Driven Zero-Range process with Long Jumps . J. Stat. Phys. 189, no. 3, (2022), Paper No. 32, 32. doi:10.1007/s10955-022-02987-3 https://doi.org/10.1007/s10955-022-02987-3

Show all 25 references
  1. [9]

    Baldasso , O

    R. Baldasso , O. Menezes , A. Neumann , and R. R. Souza . Exclusion Process with Slow Boundary . J. Stat. Phys. 167, no. 5, (2017), 1112--1142. doi:10.1007/s10955-017-1763-5 https://doi.org/10.1007/s10955-017-1763-5

  2. [10]

    Floreani and A

    S. Floreani and A. G. Casanova . Non-Equilibrium Steady State of the Symmetric Exclusion Process with Reservoirs . ArXiv e-prints (2023). doi:10.48550/arXiv.2307.02481 https://doi.org/10.48550/arXiv.2307.02481

  3. [11]

    Franco , P

    T. Franco , P. Gonçalves , and A. Neumann . Equilibrium Fluctuations for the Slow Boundary Exclusion Process . In From Particle Systems to Partial Differential Equations , vol. 209 of Springer Proc. Math. Stat., 177--197. Springer, Cham, 2017. doi:10.1007/978-3-319-66839-0_9 h...

  4. [12]

    Franco , P

    T. Franco , P. Gonçalves , and A. Neumann . Non-Equilibrium and Stationary Fluctuations of a Slowed Boundary Symmetric Exclusion . Stochastic Process. Appl. 129, no. 4, (2019), 1413--1442. doi:10.1016/j.spa.2018.05.005 https://doi.org/10.1016/j.spa.2018.05.005

  5. [13]

    Fr\'ometa , R

    S. Fr\'ometa , R. Misturini , and A. Neumann . The Boundary Driven Zero-Range Process . In From Particle Systems to Partial Differential Equations , vol. 352 of Springer Proc. Math. Stat., 253--281. Springer, Cham, 2021. doi:10.1007/978-3-030-69784-6_12 https://doi.org/10.1007...

  6. [14]

    Funaki , P

    T. Funaki , P. Van Meurs , S. Sethuraman , and K. Tsunoda . Motion by Mean Curvature from Glauber-Kawasaki Dynamics with Speed Change . J. Stat. Phys. 190, no. 3, (2023), Paper No. 45, 30. doi:10.1007/s10955-022-03044-9 https://doi.org/10.1007/s10955-022-03044-9

  7. [15]

    Gonçalves

    P. Gonçalves . Hydrodynamics for Symmetric Exclusion in Contact with Reservoirs . In Stochastic Dynamics Out of Equilibrium , vol. 282 of Springer Proc. Math. Stat., 137--205. Springer, Cham, 2019. doi:10.1007/978-3-030-15096-9\\_4 https://doi.org/10.1007/978-3-030-15096-9\

  8. [16]

    M. Z. Guo , G. C. Papanicolaou , and S. R. S. Varadhan . Nonlinear Diffusion Limit for a System with Nearest Neighbor Interactions . Comm. Math. Phys. 118, no. 1, (1988), 31--59. http://projecteuclid.org/euclid.cmp/1104161907 http://projecteuclid.org/euclid.cmp/1104161907

  9. [17]

    Gonçalves and S

    P. Gonçalves and S. Scotta . Diffusive to Super-Diffusive Behavior in Boundary Driven Exclusion . Markov Process. Related Fields 28, no. 1, (2022), 149--178

  10. [18]

    Kipnis and C

    C. Kipnis and C. Landim . Scaling Limits of Interacting Particle Systems. Springer-Verlag, Berlin, 1999. doi:10.1007/978-3-662-03752-2 https://doi.org/10.1007/978-3-662-03752-2

  11. [19]

    T. M. Liggett . Interacting Particle Systems . Classics in Mathematics. Springer-Verlag, Berlin, 2005. doi:10.1007/b138374 https://doi.org/10.1007/b138374

  12. [20]

    Levine , D

    E. Levine , D. Mukamel , and G. M. Sch\"utz . Zero-Range Process with Open Boundaries . J. Stat. Phys. 120, no. 5-6, (2005), 759--778. doi:10.1007/s10955-005-7000-7 https://doi.org/10.1007/s10955-005-7000-7

  13. [21]

    Landim , S

    C. Landim , S. Sethuraman , and S. Varadhan . Spectral Gap for Zero-Range Dynamics . Ann. Probab. 24, no. 4, (1996), 1871--1902. doi:10.1214/aop/1041903209 https://doi.org/10.1214/aop/1041903209

  14. [22]

    Marahrens , A

    D. Marahrens , A. Menegaki , and C. Mouhot . A Consistency-Stability Approach to Scaling Limits of Zero-Range Processes . ArXiv e-prints (2024). doi:10.48550/arXiv.2412.16714 https://doi.org/10.48550/arXiv.2412.16714

  15. [23]

    B. Morris . Spectral Gap for the Zero Range Process with Constant Rate . Ann. Probab. 34, no. 5, (2006), 1645--1664. doi:10.1214/009117906000000304 https://doi.org/10.1214/009117906000000304

  16. [24]

    Nagahata

    Y. Nagahata . Spectral Gap for Zero-Range Processes with Jump Rate g(x)=x^ . Stochastic Process. Appl. 120, no. 6, (2010), 949--958. doi:10.1016/j.spa.2010.01.019 https://doi.org/10.1016/j.spa.2010.01.019

  17. [25]

    F. Spitzer . Interaction of Markov Processes . Advances in Math. 5, (1970), 246--290. doi:10.1016/0001-8708(70)90034-4 https://doi.org/10.1016/0001-8708(70)90034-4

Pith tools

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