pith. sign in

arxiv: 2606.30439 · v1 · pith:X6VU4A4Bnew · submitted 2026-06-29 · 🧮 math.DS

On a Generalized Compartment Model for Ethanol Metabolism in the Human Body

Manh Tuan Hoang , Thi Kim Quy Ngo , Benjamin Wacker This is my paper

Pith reviewed 2026-06-30 03:56 UTC · model grok-4.3

classification 🧮 math.DS
keywords ethanol metabolismcompartment modelglobal asymptotic stabilityLyapunov functionnonlinear rate functionsdiscrete-time modelpositivity and boundedness
0
0 comments X

The pith

Generalized compartment model for ethanol metabolism proves global asymptotic stability of its unique equilibrium using a quadratic Lyapunov function.

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

The paper extends a compartment model of ethanol in the human body by replacing the standard Michaelis-Menten liver metabolism rate with a general class of nonlinear functions. This change is shown to preserve positivity and boundedness of solutions while allowing a quadratic Lyapunov function to establish that the unique equilibrium is globally asymptotically stable. A discrete-time version of the model is constructed and proven to match the continuous version's qualitative behavior when the time step meets an explicit upper bound. Numerical experiments with several different nonlinear rate functions are used to illustrate the results. The work therefore supplies both greater modeling flexibility and rigorous guarantees on long-term dynamics.

Core claim

The generalized continuous-time compartment model with nonlinear ethanol metabolism rate functions has positive and bounded solutions and possesses a unique equilibrium that is globally asymptotically stable, as shown by constructing an appropriate quadratic Lyapunov function whose derivative is negative definite along trajectories. The corresponding discrete-time model reproduces positivity, boundedness, and global asymptotic stability whenever the time step size satisfies a derived restriction.

What carries the argument

The general class of nonlinear rate functions for hepatic ethanol metabolism, which replaces the Michaelis-Menten form and permits construction of a quadratic Lyapunov function proving global stability.

If this is right

  • Solutions remain positive and bounded for nonnegative initial conditions.
  • The unique equilibrium is globally asymptotically stable under the stated conditions on the rate functions.
  • The discrete-time model inherits positivity, boundedness, and global stability when the time step is sufficiently small.
  • Numerical simulations with varied nonlinear rate functions confirm convergence to the equilibrium.

Where Pith is reading between the lines

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

  • The same Lyapunov construction may apply to other compartment models whose elimination rates obey comparable growth bounds.
  • The time-step restriction supplies a practical rule for choosing simulation step sizes that preserve qualitative behavior.

Load-bearing premise

The nonlinear rate functions must satisfy positivity, continuity, and suitable monotonicity or growth conditions so that a unique equilibrium exists and a quadratic Lyapunov function can be built with negative definite derivative.

What would settle it

A concrete nonlinear rate function satisfying the paper's technical conditions for which numerical integration shows either multiple equilibria or trajectories that fail to converge to the predicted equilibrium.

Figures

Figures reproduced from arXiv: 2606.30439 by Benjamin Wacker, Manh Tuan Hoang, Thi Kim Quy Ngo.

Figure 1
Figure 1. Figure 1: The flow of ethanol through the human body. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Monotone and nonmonotone rate functions. [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The graphs of the rate functions in Table [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The solutions in the phase spaces with the parameter Set [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: The solutions in the phase spaces with the parameter Set [PITH_FULL_IMAGE:figures/full_fig_p012_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: The solutions in the phase spaces with the parameter Set [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: The solutions in the phase spaces with the parameter Set [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The B components corresponding to the functions f1, f2, f3, f4 in [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The C components corresponding to the functions f1, f2, f3, f4 in [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: The solutions of the discrete-time model with the initial data [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: The solutions of the discrete-time model with the initial data [PITH_FULL_IMAGE:figures/full_fig_p015_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: The solutions of the discrete-time model with the initial data [PITH_FULL_IMAGE:figures/full_fig_p016_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: The solutions of the discrete-time model with the initial data [PITH_FULL_IMAGE:figures/full_fig_p016_13.png] view at source ↗
read the original abstract

We introduce a generalized continuous-time compartment model of ethanol metabolism in the human body that extends a recently developed framework. In the proposed model, we replace the Michaelis-Menten mechanism of the liver's ethanol metabolism rate with a general class of nonlinear rate functions. This modification provides greater modeling flexibility and enables the model to capture a wider range of hepatic ethanol metabolism dynamics. The qualitative behavior of the proposed ethanol metabolism model is analyzed rigorously. More specifically, we investigate the positivity and boundedness of solutions, as well as the global asymptotic stability (GAS) of the unique equilibrium point using an appropriate quadratic Lyapunov function. Second, we formulate a discrete-time counterpart of the proposed continuous-time model and investigate its dynamical properties. We show that, under an appropriate condition on the time step size, the discrete-time model faithfully reproduces the qualitative dynamical behavior of the corresponding continuous-time system. Lastly, we conduct a series of numerical experiments employing several ethanol metabolism rate functions to support the theoretical results.

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

3 major / 2 minor

Summary. The paper introduces a generalized continuous-time compartment model for ethanol metabolism, replacing the Michaelis-Menten liver rate with a general class of nonlinear rate functions to increase modeling flexibility. It claims to prove positivity and boundedness of solutions, existence of a unique equilibrium, and its global asymptotic stability via a quadratic Lyapunov function. A discrete-time counterpart is formulated, shown to reproduce the continuous dynamics under a suitable time-step restriction, and supported by numerical experiments with several rate functions.

Significance. If the technical hypotheses on the rate class are made explicit and the Lyapunov analysis is fully detailed and verified, the work would provide a flexible yet rigorously analyzed framework for compartment models in mathematical biology. The combination of continuous/discrete analysis and numerics is a positive feature, but the current presentation leaves the scope and verifiability of the central claims unclear.

major comments (3)
  1. [Model formulation] Model formulation section: The general class of nonlinear ethanol metabolism rate functions is introduced without an explicit list of the technical conditions (positivity, continuity, monotonicity or growth bounds) required for existence of a unique equilibrium and for the subsequent qualitative analysis. This premise is load-bearing for all claims of positivity, boundedness, and GAS.
  2. [Global asymptotic stability analysis] Global asymptotic stability section: The abstract and text refer to an 'appropriate quadratic Lyapunov function' but provide neither its explicit form nor the derivative calculation showing negative-definiteness along trajectories. Without this, the GAS claim cannot be verified.
  3. [Discrete-time model] Discrete-time model section: The 'appropriate condition on the time step size' that ensures the discrete system reproduces the continuous qualitative behavior is stated but not given explicitly or derived; this condition is central to the discrete-time fidelity claim.
minor comments (2)
  1. [Abstract] Abstract: The summary of results should reference the specific theorems or conditions rather than using vague phrases such as 'an appropriate condition on the time step size'.
  2. [Numerical experiments] Numerical experiments: The choice of parameter values and initial conditions in the simulations should be justified with reference to physiological ranges or prior literature.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments, which help clarify the presentation of our results. We address each major comment below and will revise the manuscript to improve explicitness and verifiability while preserving the core contributions.

read point-by-point responses

  1. Referee: [Model formulation] Model formulation section: The general class of nonlinear ethanol metabolism rate functions is introduced without an explicit list of the technical conditions (positivity, continuity, monotonicity or growth bounds) required for existence of a unique equilibrium and for the subsequent qualitative analysis. This premise is load-bearing for all claims of positivity, boundedness, and GAS.

    Authors: We agree that the technical conditions on the rate functions should be stated explicitly rather than left implicit. In the revised manuscript we will add a dedicated paragraph or subsection listing the precise assumptions (positivity, continuity, monotonicity, and growth bounds) that are used to guarantee existence of a unique equilibrium and to carry out the positivity, boundedness, and stability analysis. revision: yes

  2. Referee: [Global asymptotic stability analysis] Global asymptotic stability section: The abstract and text refer to an 'appropriate quadratic Lyapunov function' but provide neither its explicit form nor the derivative calculation showing negative-definiteness along trajectories. Without this, the GAS claim cannot be verified.

    Authors: The referee correctly notes that the explicit quadratic Lyapunov function and the derivative calculation are not supplied in the current text. This omission prevents independent verification. In the revision we will state the precise form of the quadratic Lyapunov function and provide the full computation showing that its derivative is negative definite along trajectories under the model assumptions. revision: yes

  3. Referee: [Discrete-time model] Discrete-time model section: The 'appropriate condition on the time step size' that ensures the discrete system reproduces the continuous qualitative behavior is stated but not given explicitly or derived; this condition is central to the discrete-time fidelity claim.

    Authors: We acknowledge that the time-step restriction is mentioned only qualitatively. In the revised version we will derive and state the explicit upper bound on the time step h, together with the inequalities that guarantee preservation of positivity, boundedness, and global asymptotic stability for the discrete-time system. revision: yes

Circularity Check

0 steps flagged

No significant circularity

full rationale

The derivation consists of standard qualitative analysis (positivity, boundedness, GAS via quadratic Lyapunov) for a stated class of nonlinear rate functions satisfying explicit technical hypotheses (positivity, continuity, monotonicity/growth bounds). The discrete-time counterpart is shown to preserve the same properties under a time-step restriction. No step reduces by construction to a fitted parameter, self-definition, or load-bearing self-citation chain; the central claims remain conditional on the model-class assumptions and are independent of any internal renaming or ansatz smuggling. This is the normal case of a self-contained mathematical analysis.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claims rest on the assumption that the nonlinear rate functions belong to a sufficiently regular class permitting a unique equilibrium and a quadratic Lyapunov function with the required sign properties; no free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption The ethanol metabolism rate functions are positive, continuous, and satisfy the growth and monotonicity conditions required for existence of a unique equilibrium and for the derivative of the quadratic Lyapunov function to be negative definite.
    Invoked to establish positivity, boundedness, and global asymptotic stability of the continuous-time model.

pith-pipeline@v0.9.1-grok · 5701 in / 1369 out tokens · 68611 ms · 2026-06-30T03:56:43.448058+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

31 extracted references · 5 canonical work pages

  1. [1]

    L. J. S. Allen, Some discrete-time SI, SIR, and SIS epidemic models, Math. Biosci. 124 (1994) 83-125

    1994

  2. [2]

    L. J. S. Allen, An introduction to mathematical biology, Pearson/Prentice Hall, 2007

    2007

  3. [3]

    U. M. Ascher, L. R. Petzold, Computer Methods for Ordinary Differential Equations and Differential-Algebraic Equations, SIAM, Philadelphia, 1998

    1998

  4. [4]

    Capasso, G

    V. Capasso, G. Serio, A generalization of the Kermack-McKendrick deterministic epi- demic model, Mathematical Biosciences 42(1978) 43-61

    1978

  5. [5]

    J. Cui, X. Mu, H. Wan, Saturation recovery leads to multiple endemic equilibria and backward bifurcation, Journal of Theoretical Biology 254 (2008) 275-283. 17

    2008

  6. [6]

    Dawes, M.O

    J.H.P. Dawes, M.O. Souza, A derivation of Holling’s type I, II and III functional responses in predator-prey systems, Journal of Theoretical Biology 327:11–22, 2013, https://doi.org/10.1016/j.jtbi.2013.02.017

    work page doi:10.1016/j.jtbi.2013.02.017 2013

  7. [7]

    Edenberg, J.N

    H.J. Edenberg, J.N. McClintick, Alcohol Dehydrogenases, Aldehyde Dehydrogenases, and Alcohol Use Disorders: A Critical Review, Alcoholism: Clinical and Experimental Research 42(12):2281–2297, 2018, https://doi.org/10.1111/acer.13904

    work page doi:10.1111/acer.13904 2018

  8. [8]

    Elaydi, An introduction to difference equations, Springer-Verlag, New York, 2005

    S. Elaydi, An introduction to difference equations, Springer-Verlag, New York, 2005

    2005

  9. [9]

    Kalinkat, B.C

    G. Kalinkat, B.C. Rall, S.F. Uiterwall, W. Uszko, Empirical evidence of type III func- tional responses and why it remains rare, Frontiers in Ecology and Evolution 11:1033818, 2023, https://doi.org/10.3389/fevo.2023.1033818

    work page doi:10.3389/fevo.2023.1033818 2023

  10. [10]

    H. K. Khalil, Nonlinear systems, Third Edition, Prentice Hall, 2002

    2002

  11. [11]

    L. M. Ladino, C. Mammana, E. Michetti, J. C. Valverde, Discrete time population dynamics of a two-stage species with recruitment and capture, Chaos, Solitons and Fractals 85 (2016) 143-150

    2016

  12. [12]

    M. D. Levitt, D. G. Levitt, Use of a Two-Compartmental Model to Assess Pharma- cokinetics of Human Ethanol Metabolism, Alcoholism: Clinical and Experimental Re- search 22, no. 8 (1998): 1680-1688,https://doi.org/10.1111/j.1530-0277.1998. tb03966.x

    work page doi:10.1111/j.1530-0277.1998 1998

  13. [13]

    J. Li, Z. Ma , F. Brauer, Global analysis of discrete-time SI and SIS epidemic models, Mathematical Biosciences and Engineering 2007(4) 699-710

    2007

  14. [14]

    Ludwin, Blood Alcohol Content, Undergraduate Journal of Mathematical Modeling: One + two 3, no

    C. Ludwin, Blood Alcohol Content, Undergraduate Journal of Mathematical Modeling: One + two 3, no. 2 (2011): 1, https://doi.org/10.5038/2326-3652.3.2.1

    work page doi:10.5038/2326-3652.3.2.1 2011

  15. [15]

    R. M. May, Simple mathematical models with very complicated dynamics, Nature 261(1976) 459-467

    1976

  16. [16]

    R. E. Mickens, Nonstandard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994

    1994

  17. [17]

    R. E. Mickens, Applications of Nonstandard Finite Difference Schemes, World Scientific, Singapore, 2000

    2000

  18. [18]

    R. E. Mickens, Dynamic consistency: a fundamental principle for constructing nonstan- dard finite difference schemes for differential equations, Journal of Difference Equations and Applications 11 (2005) 645-653

    2005

  19. [19]

    R. E. Mickens, Advances in the Applications of Nonstandard Finite Difference Schemes, World Scientific, Singapore, 2005

    2005

  20. [20]

    R. E. Mickens, Nonstandard Finite Difference Schemes: Methodology and Applications, World Scientific, 2020. 18

    2020

  21. [21]

    Ruana, W

    S. Ruana, W. Wang, Dynamical behavior of an epidemic model with a nonlinear inci- dence rate, Journal of Differential Equations 188(2003) 135-163

    2003

  22. [22]

    Seibert, R

    P. Seibert, R. Suarez, Global stabilization of nonlinear cascade systems, Syst. Contr. Lett. 14 (1990), 347–352

    1990

  23. [23]

    H. L. Smith and P. Waltman, The Theory of the Chemostat: Dynamics of Microbial Competition, Cambridge University Press, 2009

    2009

  24. [24]

    A. M. Stuart, A. R. Humphries, Dynamical systems and numerical analysis, Cambridge university press, 1998

    1998

  25. [25]

    B. Wacker, Analysis of a Finite-Difference Method Based on Nonlocal Approximations for a Nonlinear, Extended Three-Compartmental Model of Ethanol Metabolism in the Human Body, Mathematical Methods in the Applied Sciences 48 (2025), 9975-9992

    2025

  26. [26]

    Wacker, Global asymptotic stability of one time-continuous and two time-discrete models for ethanol metabolism in the human body, Aequat

    B. Wacker, Global asymptotic stability of one time-continuous and two time-discrete models for ethanol metabolism in the human body, Aequat. Math. 99 (2025), 2401- 2422

    2025

  27. [27]

    X. Wang, Z. Wang, H. Shen, Dynamical analysis of a discrete-time SIS epidemic model on complex networks, Applied Mathematics Letters 94 (2019) 292-299

    2019

  28. [28]

    X. Wu, Z. Zhou, F. Xie, Multi-scale dynamics of a piecewise-smooth Bazykin’s prey- predator system, Nonlinear Dynamics 113:1969–1981, 2025, https://doi.org/10. 1007/s11071-024-10292-9

    1969

  29. [29]

    D. Xiao, S. Ruan, Global analysis of an epidemic model with nonmonotone incidence rate, Mathematical Biosciences 208 (2007) 419-429

    2007

  30. [30]

    R. Xu, Z. Ma, Global stability of a delayed SEIRS epidemic model with saturation incidence rate, Nonlinear Dyn (2010) 61: 229-239

    2010

  31. [31]

    Zhang, X

    X. Zhang, X. Liu, Backward bifurcation of an epidemic model with saturated treatment function, J. Math. Anal. Appl. 348 (2008) 433-443. 19

    2008