[论文解读] New results concerning the stability of equilibria of a delay differential equation modeling leukemia

Anca Veronica Ion|arXiv (Cornell University)|Jan 26, 2010

Chronic Myeloid Leukemia Treatments参考文献 6被引用 6

一句话总结

本文修正并完善了慢性骨髓性白血病延迟微分方程模型中平衡点稳定性的分析。通过李雅普诺夫方法与特征方程分析，提供了局部稳定性和霍普夫分支的准确条件，纠正了先前研究中关于参数依赖性和分支阈值的错误，特别是针对正平衡点的分析。

ABSTRACT

The paper is devoted to the study of stability of equilibrium solutions of a delay differential equation that models leukemia. The equation was previously studied in [5] and [6], where the emphasis is put on the numerical study of periodic solutions. Some stability results for the equilibria are also presented in these works, but they are incomplete and contain some errors. Our work aims to complete and to bring corrections to those results. Both Lyapunov first order approximation method and second Lyapunov method are used.

研究动机与目标

• 修正先前文献[5]和[6]中关于白血病延迟微分方程模型平衡点稳定性分析的不准确之处。
• 为该模型的零平衡点和正平衡点提供完整且严格的稳定性分析。
• 阐明延迟r、增殖率k和希尔系数n等参数在决定稳定性与分支行为中的作用。
• 利用李雅普诺夫方法与线性化系统特征值分析，建立霍普夫分支与局部稳定性的正确条件。

提出的方法

• 利用抽象泛函微分方程理论与区间[−r, 0]上连续函数空间，严格证明了在[−r, ∞)上有界解的存在性与唯一性。
• 推导出平衡解为x₁ = 0和x₂ = [(β₀/δ)(k−1)−1]^{1/n}，其生物学可行性取决于r < r_max。
• 通过特征方程λ + B + δ − kB e^{−λr} = 0分析每个平衡点附近的线性化系统，其中B依赖于平衡点与参数。
• 采用第一类与第二类李雅普诺夫方法评估稳定性，重点关注特征值实部的符号以及霍普夫分支的首阶李雅普诺夫系数。
• 通过函数g(r) = T^{−1}(−(δ+B)r) − arccos((δ+B)/(kB))分析临界延迟r*处稳定性变化的条件，以确定符号变化。
• 通过验证特征值关于延迟的导数非零以及首阶李雅普诺夫系数非零，确立霍普夫分支的条件。

实验结果

研究问题

• RQ1 在白血病延迟微分方程模型中，零平衡点与正平衡点的正确稳定性条件是什么？
• RQ2 为何[5]和[6]中的稳定性结果不完整且错误，特别是关于希尔系数n与延迟r的依赖性？
• RQ3 B与δ + B的符号如何影响正平衡点的稳定性，其对分支行为有何影响？
• RQ4 霍普夫分支发生的精确条件是什么？它与先前研究中提出的错误条件有何不同？
• RQ5 参数r如何影响稳定与不稳定动力学之间的转换？参数k = 2e^{−γr}在此过程中起什么作用？

主要发现

• 零平衡点x₁ = 0当且仅当k < 1时稳定，与希尔系数n无关，这与先前认为n ∈ [0,1]会影响其稳定性的观点相矛盾。
• 正平衡点x₂ > 0仅在r < r_max = −(1/γ)ln(½(1 + δ/β₀))时存在，此条件要求δ/β₀ < 1。
• 对于正平衡点，稳定性取决于B的符号：当B > 0时，稳定性成立的条件为β₀(k−1)/δ > 1，等价于x₂ > 0。
• 当B < 0且δ + B > 0时，稳定性条件为r < arccos((δ + B)/(kB)) / ω₀，其中ω₀ = √(k²B² − (δ + B)²)，此结果修正了先前研究中错误使用√(k²B² − (δ + B)²)的问题。
• 本文指出，[5]和[6]中的稳定性域是错误的，因其使用了错误的ω₀表达式，导致分支阈值不准确。
• 当g(r*) = 0且dμ/dα(α*) ≠ 0时，霍普夫分支在r*处发生，稳定性由中心流形上首阶李雅普诺夫系数的符号决定。