2025-11-23T07:58:16.677890

Periodic solutions in a tumor-immune competition system with time-delay and chemotherapy effects

Amster, Rivera, Arredondo
The main purpose of this paper is to analyze the dynamics of the system of time-delay differential equations (DDEs) \begin{equation*} \begin{split} \dot{T}(t)&=T(t) f(t,T(t))-γE(t)T(t),\\ \dot{E}(t)&=σ+ \frac{pE(t)T(t-τ_1)}{g+a T(t-τ_1)}-\frac{mE(t)T(t-τ_2)}{g+a T(t-τ_2)}-ηE(t), \end{split} \end{equation*} where $T=T(t)$ and $E=E(t)$ represent the concentrations of tumor and effector cells at the time $t$. The coefficients $σ$, $μ$, $γ$, and $η$ are all positive, and $f(t, T)$ represents the relative growth rate of tumor cells, corresponding to a generalized logistic growth function that describes periodic time chemotherapeutic effects. The parameter $τ_1 \in \mathbb{R}_{\ge 0}$ is the response time delay of the immune system (mediated by effector cells) to an invasion of tumor cells, while $τ_2 \in \mathbb{R}_{\ge 0}$ represents the time delay of tumor cells in response to the appearance of effector cells.
academic

Periodic solutions in a tumor-immune competition system with time-delay and chemotherapy effects

Basic Information

  • Paper ID: 2510.11135
  • Title: Periodic solutions in a tumor-immune competition system with time-delay and chemotherapy effects
  • Authors: Pablo Amster, Andrés Rivera, John A. Arredondo
  • Classification: math.DS (Dynamical Systems)
  • Publication Date: October 13, 2025
  • Paper Link: https://arxiv.org/abs/2510.11135

Abstract

This paper primarily analyzes the dynamical properties of a system of differential equations with time delays:

\dot{T}(t) &= T(t)f(t,T(t)) - \gamma E(t)T(t) \\ \dot{E}(t) &= \sigma + \frac{pE(t)T(t-\tau_1)}{g+aT(t-\tau_1)} - \frac{mE(t)T(t-\tau_2)}{g+aT(t-\tau_2)} - \eta E(t) \end{align}$$ where $T(t)$ and $E(t)$ represent the concentrations of tumor cells and effector cells at time $t$, respectively; $f(t,T)$ denotes the relative growth rate of tumor cells, corresponding to a generalized logistic growth function describing periodic chemotherapy effects; $\tau_1$ is the response delay of the immune system to tumor cell invasion, and $\tau_2$ is the response delay of tumor cells to the appearance of effector cells. ## Research Background and Motivation ### Problem Background 1. **Tumor-immune system interactions**: Tumors are abnormal growths of body tissues that are normally controlled by the immune system; however, tumor formation is often associated with immune system deficiencies 2. **Importance of time-delay effects**: Time delays are ubiquitous in biological systems; the immune system requires time to recognize and respond to tumors, and tumors also exhibit delayed responses to immune attacks 3. **Periodic effects of chemotherapy**: Chemotherapy typically employs periodic dosing schedules, necessitating consideration of its time-varying effects on tumor growth ### Research Motivation 1. **Limitations of existing models**: Classical tumor-immune models typically employ logistic growth, but experimental data show that Gompertz models provide better fitting 2. **Impact of delay parameters**: A deeper analysis of how different delay parameters affect system dynamics is needed 3. **Existence of periodic solutions**: Exploring the conditions for the existence of periodic solutions under periodic chemotherapy effects ## Core Contributions 1. **Established a tumor-immune time-delay model with generalized Richards-type growth**, incorporating two distinct delay parameters $\tau_1$ and $\tau_2$ 2. **Provided comprehensive analysis of equilibrium point existence and stability**, offering classification of equilibrium points under different parameter conditions 3. **Determined critical delay values for stability switching**, providing explicit stability intervals 4. **Proved the existence of periodic solutions**, establishing periodic solution branches from autonomous systems to non-autonomous systems via continuation methods 5. **Provided conditions for Hopf bifurcation**, analyzing the system's transition from stable equilibrium points to periodic solutions ## Methodology Details ### Problem Definition Investigation of the dynamical properties of a dual-delay tumor-immune competition system, including: - Existence and stability analysis of equilibrium points - Impact of delay parameters on system stability - Proof of periodic solution existence - Hopf bifurcation analysis ### Model Architecture #### Basic Model Framework System (2) describes the interaction between tumor cells T(t) and effector cells E(t): **Tumor cell equation**: $$\dot{T}(t) = T(t)f(t,T(t)) - \gamma E(t)T(t)$$ - First term: natural growth of tumor cells - Second term: killing of tumor cells by effector cells **Effector cell equation**: $$\dot{E}(t) = \sigma + \frac{pE(t)T(t-\tau_1)}{g+aT(t-\tau_1)} - \frac{mE(t)T(t-\tau_2)}{g+aT(t-\tau_2)} - \eta E(t)$$ - $\sigma$: constant input of effector cells - Second term: activation of effector cells by tumor cells (Holling II functional response) - Third term: exhaustion of effector cells - Fourth term: natural death of effector cells #### Growth Function Design Richards-type growth function is employed: $$f(t,T) = r\left[\frac{(1-T^\beta)}{\beta} - b(t)\right]$$ where: - $\beta \in (0,1]$: shape parameter - $b(t)$: periodic chemotherapy function - When $\beta=1$, it reduces to logistic growth - As $\beta \to 0$, it approaches Gompertz growth ### Technical Innovations 1. **Dual-delay mechanism**: Distinguishes between immune activation delay $\tau_1$ and immune exhaustion delay $\tau_2$, more accurately reflecting biological processes 2. **Generalized growth model**: Unifies logistic and Gompertz growth modes through parameter $\beta$ 3. **Parameter scaling technique**: Simplifies parameter analysis through transformation (14) 4. **Stability switching curves**: Provides stability boundaries in the $(\tau_1,\tau_2)$ plane ## Theoretical Analysis ### Equilibrium Point Analysis #### Tumor-free Equilibrium Point For all parameter values, the system always possesses a tumor-free equilibrium point $(0,\sigma/\eta)$. **Stability conditions** (Theorem 10): Define $\Delta := \gamma\sigma - rb\eta$, then: - If $\Delta < 0$: tumor-free equilibrium point is unstable - If $\Delta = 0$: tumor-free equilibrium point is stable - If $\Delta > 0$: tumor-free equilibrium point is asymptotically stable #### Coexistence Equilibrium Point When $a=0$, the coexistence equilibrium point $(T^*,E^*)$ satisfies equation (17): $$\mu(T^{*\beta} - b)T^* + T^{*\beta} = b - \sigma$$ **Existence conditions** (Proposition 8): - Condition †: $\sigma < b$ and $\mu \leq 0$ - Condition ‡: $\sigma \leq b$ and $0 < \mu \leq \mu_c$ where $\mu_c := \left(\frac{1}{b}\left(\frac{1-\beta}{1+\beta}\right)^{\beta-1}\right)^{1/\beta}$ ### Stability Analysis #### Small Delay Case **Theorem 12**: For sufficiently small $a,\tau_1,\tau_2 \geq 0$: 1. Under condition †, $(T^*,E^*)$ is locally asymptotically stable 2. Under condition ‡, if additional conditions are satisfied, $(T^*,E^*)$ is stable #### General Delay Case **Theorem 14**: When $\tau_1 = \tau_2 = \tau$, there exists a critical value $\tau_c$ such that: - $\tau < \tau_c$: equilibrium point is stable - $\tau > \tau_c$: equilibrium point loses stability and undergoes Hopf bifurcation The critical value is determined by the following equation: $$G(y) = y^4 + \left[(\lambda_1^*)^2 + (\lambda_2^*)^2\right]y^2 + (\lambda_1^*\lambda_2^*)^2 - N^2 = 0$$ ### Periodic Solution Existence **Theorem 16** (Continuation Theorem): Let $(T^*,E^*)$ be a feasible equilibrium point when $b \equiv b_0$. Then there exist constants $\phi^*,\tau_j^* > 0$ and $r^* > 0$ such that when $\|\phi\|_\infty < \phi^*$, $\tau_1 < \tau_1^*$, $\tau_2 < \tau_2^*$, system (2) possesses at least one $\omega$-periodic solution $(T,E)$ satisfying $\|(T-T^*,E-E^*)\|_\infty < r^*$. **Proof Strategy**: 1. Verify that the linearized system has no non-trivial periodic solutions (Lemma 4) 2. Apply the implicit function theorem and Fredholm alternative 3. Establish smooth branches from autonomous delay-free systems to non-autonomous delay systems ## Numerical Verification ### Parameter Settings The paper employs typical parameter values from the literature: - $g \approx 10^6$ (relatively large) - After scaling, $a \to a/g = O(10^{-6})$ (relatively small) ### Stability Switching Curves **Theorem 15** provides stability switching curves in the $(\tau_1,\tau_2)$ plane: $$C_{s,k} = \{(\tau_1^{\pm}(y), \tau_2^{\pm}(y)) \in \mathbb{R}_+^2 : y \in I, s,k \in \mathbb{Z}\}$$ where $\tau_1^{\pm}$ and $\tau_2^{\pm}$ are given by equations (35) and (36). ## Related Work ### Tumor Growth Models 1. **Classical models**: Exponential growth, logistic growth, Gompertz growth 2. **Experimental validation**: Multiple studies demonstrate that Gompertz models provide the best fit to experimental data 3. **Parameter correlations**: The two parameters in Gompertz models exhibit significant correlation ### Tumor-Immune Interaction Models 1. **Kuznetsov et al. (1994)**: Established the foundational framework for tumor-immune models 2. **Li et al. (2021, 2022)**: Investigated single and dual delay effects 3. **This paper's contribution**: Extended to generalized Richards growth and periodic chemotherapy effects ### Time-Delay Differential Equation Theory 1. **Stability analysis**: Characteristic equation method, Lyapunov function method 2. **Bifurcation theory**: Hopf bifurcation, stability switching 3. **Periodic solutions**: Continuation methods, fixed-point theory ## Conclusions and Discussion ### Main Conclusions 1. **Complete equilibrium point classification**: Provided existence conditions and classification for all possible equilibrium points 2. **Clear stability conditions**: Offered explicit stability criteria and critical delay values 3. **Periodic solution existence**: Proved the existence of periodic solutions under periodic chemotherapy effects 4. **Clear bifurcation phenomena**: Revealed the transition mechanism from stability to oscillation ### Limitations 1. **Parameter simplification**: The assumption $a=0$, while biologically motivated, limits generality 2. **Insufficient numerical verification**: Theoretical results lack adequate numerical simulation verification 3. **Biological validation**: Model predictions require more experimental data support ### Future Directions 1. **Three-dimensional extension**: Consider additional cell types (e.g., regulatory T cells) 2. **Spatial heterogeneity**: Introduce diffusion terms to study spatial patterns 3. **Stochastic effects**: Consider noise effects on system dynamics 4. **Optimal control**: Design optimal chemotherapy schemes based on the model ## In-Depth Evaluation ### Strengths 1. **Theoretical rigor**: Complete mathematical analysis with clear proof procedures 2. **Reasonable model**: Clear biological background and well-defined parameter meanings 3. **Methodological innovation**: Novel combination of dual-delay mechanisms and generalized growth functions 4. **Practical results**: Provides theoretical guidance for chemotherapy scheme design ### Weaknesses 1. **Weak numerical section**: Section 6 "Numerical Verification" is marked as "under construction" 2. **Missing parameter estimation**: No methods provided for actual parameter estimation 3. **Insufficient model validation**: Lacks comparison with clinical data ### Impact 1. **Theoretical contribution**: Enriches applications of time-delay differential equations in biomathematics 2. **Application value**: Provides new perspectives for mathematical modeling of tumor treatment 3. **Methodological significance**: The application of continuation methods in periodic solution analysis has general applicability ### Applicable Scenarios 1. **Theoretical research**: Time-delay differential equations, dynamical systems theory 2. **Biomathematics**: Tumor growth models, immune system modeling 3. **Medical applications**: Chemotherapy scheme optimization, treatment efficacy prediction ## References The paper cites 27 related references, primarily including: - Classical literature on tumor growth models - Stability theory of time-delay differential equations - Mathematical models of tumor-immune interactions - Bifurcation theory and periodic solution existence theorems --- **Overall Assessment**: This is a highly theoretical mathematical biology paper that makes important contributions to time-delay differential equations and tumor-immune system modeling. The theoretical analysis is thorough and the methodology is innovative; however, the numerical verification section requires further development. It has significant implications for understanding the complex dynamics of tumor-immune systems.