From the Formalism of the Delayed Mackey-Glass Equation and Its Application to Leukemia Growth
Research Article  ·  Published: 20 September 2026
Issue cover
Journal of Numerical Simulations in Physics and Mathematics
Volume 2, Issue 2, 2026: 130-140
Research Article Open Access

From the Formalism of the Delayed Mackey-Glass Equation and Its Application to Leukemia Growth

1 Institute for Energy and Nuclear Research (IPEN), University of São Paulo, São Paulo 05508-220, Brazil
* Corresponding Author: Matheus dos Santos Farias, [email protected]
Volume 2, Issue 2
You have full access to this open access article · CC BY 4.0 License

Article Information

Abstract

This study develops a mathematically explicit and computationally reproducible framework based on the delayed Mackey--Glass equation for leukemia-related population dynamics. The classical model is formulated from a population balance law, delayed maturation, and saturating nonlinear feedback, and is then extended by a stage-dependent treatment-control term. Parameters that are not identifiable from the available clinical information are taken from established mathematical and hematopoietic literature, while the treatment-control coefficients are constrained phenomenologically by the modeled burden states. A brief systematic literature search was used to characterize previous Mackey--Glass applications in hematopoiesis and leukemia and to define the methodological gap addressed here. A published case of acute myeloid leukemia-associated myeloid sarcoma treated with radiotherapy and azacitidine provides the clinical treatment schedule and volumetric constraint. The model output is mapped explicitly to a three-dimensional point-cloud representation through a deterministic population-to-point transformation, while the observed reduction from 560~cm$^3$ to 157~cm$^3$ after 19.8~Gy provides an independent constraint for a population-to-volume mapping. Sensitivity to the delay and Hill exponent is examined, and the delay-adapted fourth-order Runge--Kutta implementation is compared with an independent adaptive method-of-steps reference. The resulting framework integrates formal derivation, structured literature evidence, treatment modeling, numerical verification, and clinically anchored visualization without presenting the point cloud as a patient-specific radiological reconstruction.

Graphical Abstract

From the Formalism of the Delayed Mackey-Glass Equation and Its Application to Leukemia Growth

Keywords

Mackey-Glass equation leukemia growth delay differential equations mathematical oncology radiotherapy numerical simulation

Data Availability Statement

Not applicable.

Funding

This work was supported without any funding.

Conflicts of Interest

The author declares no conflicts of interest.

AI Use Statement

The author declares that no generative AI was used in the preparation of this manuscript.

Ethical Approval and Consent to Participate

This study did not involve recruitment, intervention, or new data collection from human or animal participants. The clinical values used to define the treatment scenario were obtained from a previously published case report and were used only as literature-based clinical constraints for a mathematical simulation. No identifiable patient data were collected or analyzed; therefore, no new ethical approval or informed consent was required for the present simulation study.

References

  1. American Cancer Society. (2025). Cancer facts & figures 2025. American Cancer Society. Retrieved from https://www.cancer.org/research/cancer-facts-statistics/all-cancer-facts-figures/2025-cancer-facts-figures.html
    [Google Scholar]
  2. National Cancer Institute. (n.d.). Cancer stat facts: Leukemia. Surveillance, Epidemiology, and End Results Program. Retrieved from https://seer.cancer.gov/statfacts/html/leuks.html
    [Google Scholar]
  3. Mackey, M. C., & Glass, L. (1977). Oscillation and chaos in physiological control systems. Science, 197(4300), 287-289.
    [CrossRef] [Google Scholar]
  4. Berezansky, L., Braverman, E., & Idels, L. (2013). Mackey-Glass model of hematopoiesis with monotone feedback revisited. Applied Mathematics and Computation, 219(9), 4892--4907.
    [CrossRef] [Google Scholar]
  5. Berezansky, L., Braverman, E., & Idels, L. (2013). Mackey--Glass model of hematopoiesis with non-monotone feedback: Stability, oscillation and control. Applied Mathematics and Computation, 219(11), 6268--6283.
    [CrossRef] [Google Scholar]
  6. Adimy, M., Crauste, F., & Ruan, S. (2005). A mathematical study of the hematopoiesis process with applications to chronic myelogenous leukemia. SIAM journal on applied mathematics, 65(4), 1328-1352.
    [CrossRef] [Google Scholar]
  7. Pujo-Menjouet, L. (2016). Blood cell dynamics: half of a century of modelling. Mathematical Modelling of Natural Phenomena, 11(1), 92-115.
    [CrossRef] [Google Scholar]
  8. Majdoul, S., Colson-Durand, L., To, N. H., & Belkacemi, Y. (2017). Adaptive radiotherapy for an uncommon chloroma. Case Reports in Oncology, 9(3), 593-598.
    [CrossRef] [Google Scholar]
  9. Shampine, L. F., & Thompson, S. (2001). Solving ddes in MATLAB. Applied Numerical Mathematics, 37(4), 441-458.
    [CrossRef] [Google Scholar]
  10. Crauste, F. (2009). Delay model of hematopoietic stem cell dynamics: asymptotic stability and stability switch. Mathematical Modelling of Natural Phenomena, 4(2), 28-47.
    [CrossRef] [Google Scholar]
  11. Halanay, A., Cândea, D., & Rădulescu, I. R. (2014). Existence and Stability of Limit Cycles in a Two-delays Model of Hematopoiesis Including Asymmetric Division. Mathematical Modelling of Natural Phenomena, 9(1), 58-78.
    [CrossRef] [Google Scholar]
  12. Badralexi, I., & Halanay, A. (2017, January). Stability and oscillations in a CML model. In AIP Conference Proceedings (Vol. 1798, No. 1, p. 020011). AIP Publishing LLC.
    [CrossRef] [Google Scholar]
  13. Cândea, D., Halanay, A., Rădulescu, R., & Tălmaci, R. (2017, January). Parameter estimation and sensitivity analysis for a mathematical model with time delays of leukemia. In AIP Conference Proceedings (Vol. 1798, No. 1, p. 020034). AIP Publishing LLC.
    [CrossRef] [Google Scholar]
  14. Parajdi, L. G., Precup, R., Bonci, E. A., & Tomuleasa, C. (2020). A mathematical model of the transition from normal hematopoiesis to the chronic and accelerated-acute stages in myeloid leukemia. Mathematics, 8(3), 376.
    [CrossRef] [Google Scholar]
  15. Bakst, R. L., Dabaja, B. S., Specht, L. K., & Yahalom, J. (2018). Use of radiation in extramedullary leukemia/chloroma: guidelines from the International Lymphoma Radiation Oncology Group. International Journal of Radiation Oncology* Biology* Physics, 102(2), 314-319.
    [CrossRef] [Google Scholar]
  16. Loscocco, G. G., & Vannucchi, A. M. (2023). Myeloid sarcoma: more and less than a distinct entity. Annals of Hematology, 102(8), 1973-1984.
    [CrossRef] [Google Scholar]
  17. Adimy, M., Crauste, F., & Ruan, S. (2006). Periodic oscillations in leukopoiesis models with two delays. Journal of Theoretical Biology, 242(2), 288-299.
    [CrossRef] [Google Scholar]
  18. Colijn, C., & Mackey, M. C. (2005). A mathematical model of hematopoiesis—I. Periodic chronic myelogenous leukemia. Journal of Theoretical Biology, 237(2), 117-132.
    [CrossRef] [Google Scholar]
  19. Safarishahrbijari, A., & Gaffari, A. (2013). Parameter identification of hematopoiesis mathematical model–periodic chronic myelogenous leukemia. Contemporary Oncology/Współczesna Onkologia, 17(1), 73-77.
    [CrossRef] [Google Scholar]
  20. Foley, C., & Mackey, M. C. (2009). Dynamic hematological disease: a review. Journal of mathematical biology, 58(1), 285-322.
    [CrossRef] [Google Scholar]

Cite This Article

APA Style
dos Santos Farias, M. (2026). From the Formalism of the Delayed Mackey-Glass Equation and Its Application to Leukemia Growth. Journal of Numerical Simulations in Physics and Mathematics, 2(2), 130-140. https://doi.org/10.62762/JNSPM.2026.604346
Export Citation
RIS Format
Compatible with EndNote, Zotero, Mendeley, and other reference managers
TY  - JOUR
AU  - Farias, Matheus dos Santos
PY  - 2026
DA  - 2026/09/20
TI  - From the Formalism of the Delayed Mackey-Glass Equation and Its Application to Leukemia Growth
JO  - Journal of Numerical Simulations in Physics and Mathematics
T2  - Journal of Numerical Simulations in Physics and Mathematics
JF  - Journal of Numerical Simulations in Physics and Mathematics
VL  - 2
IS  - 2
SP  - 130
EP  - 140
DO  - 10.62762/JNSPM.2026.604346
UR  - https://www.icck.org/article/abs/JNSPM.2026.604346
KW  - Mackey-Glass equation
KW  - leukemia growth
KW  - delay differential equations
KW  - mathematical oncology
KW  - radiotherapy
KW  - numerical simulation
AB  - This study develops a mathematically explicit and computationally reproducible framework based on the delayed Mackey--Glass equation for leukemia-related population dynamics. The classical model is formulated from a population balance law, delayed maturation, and saturating nonlinear feedback, and is then extended by a stage-dependent treatment-control term. Parameters that are not identifiable from the available clinical information are taken from established mathematical and hematopoietic literature, while the treatment-control coefficients are constrained phenomenologically by the modeled burden states. A brief systematic literature search was used to characterize previous Mackey--Glass applications in hematopoiesis and leukemia and to define the methodological gap addressed here. A published case of acute myeloid leukemia-associated myeloid sarcoma treated with radiotherapy and azacitidine provides the clinical treatment schedule and volumetric constraint. The model output is mapped explicitly to a three-dimensional point-cloud representation through a deterministic population-to-point transformation, while the observed reduction from 560~cm$^3$ to 157~cm$^3$ after 19.8~Gy provides an independent constraint for a population-to-volume mapping. Sensitivity to the delay and Hill exponent is examined, and the delay-adapted fourth-order Runge--Kutta implementation is compared with an independent adaptive method-of-steps reference. The resulting framework integrates formal derivation, structured literature evidence, treatment modeling, numerical verification, and clinically anchored visualization without presenting the point cloud as a patient-specific radiological reconstruction.
SN  - 3068-9082
PB  - Institute of Central Computation and Knowledge
LA  - English
ER  - 
BibTeX Format
Compatible with LaTeX, BibTeX, and other reference managers
@article{Farias2026From,
  author = {Matheus dos Santos Farias},
  title = {From the Formalism of the Delayed Mackey-Glass Equation and Its Application to Leukemia Growth},
  journal = {Journal of Numerical Simulations in Physics and Mathematics},
  year = {2026},
  volume = {2},
  number = {2},
  pages = {130-140},
  doi = {10.62762/JNSPM.2026.604346},
  url = {https://www.icck.org/article/abs/JNSPM.2026.604346},
  abstract = {This study develops a mathematically explicit and computationally reproducible framework based on the delayed Mackey--Glass equation for leukemia-related population dynamics. The classical model is formulated from a population balance law, delayed maturation, and saturating nonlinear feedback, and is then extended by a stage-dependent treatment-control term. Parameters that are not identifiable from the available clinical information are taken from established mathematical and hematopoietic literature, while the treatment-control coefficients are constrained phenomenologically by the modeled burden states. A brief systematic literature search was used to characterize previous Mackey--Glass applications in hematopoiesis and leukemia and to define the methodological gap addressed here. A published case of acute myeloid leukemia-associated myeloid sarcoma treated with radiotherapy and azacitidine provides the clinical treatment schedule and volumetric constraint. The model output is mapped explicitly to a three-dimensional point-cloud representation through a deterministic population-to-point transformation, while the observed reduction from 560~cm\$^3\$ to 157~cm\$^3\$ after 19.8~Gy provides an independent constraint for a population-to-volume mapping. Sensitivity to the delay and Hill exponent is examined, and the delay-adapted fourth-order Runge--Kutta implementation is compared with an independent adaptive method-of-steps reference. The resulting framework integrates formal derivation, structured literature evidence, treatment modeling, numerical verification, and clinically anchored visualization without presenting the point cloud as a patient-specific radiological reconstruction.},
  keywords = {Mackey-Glass equation, leukemia growth, delay differential equations, mathematical oncology, radiotherapy, numerical simulation},
  issn = {3068-9082},
  publisher = {Institute of Central Computation and Knowledge}
}

Article Metrics

Citations
Crossref
0
Scopus
0
Views
15
PDF Downloads
3

Publisher's Note

ICCK stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and Permissions

CC BY Copyright © 2026 by the Author(s). Published by Institute of Central Computation and Knowledge. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/), which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.
Journal of Numerical Simulations in Physics and Mathematics
Journal of Numerical Simulations in Physics and Mathematics
ISSN: 3068-9082 (Online) | ISSN: 3068-9074 (Print)
Portico
Preserved at
Portico