[논문 리뷰] New results concerning the stability of equilibria of a delay differential equation modeling leukemia

Anca Veronica Ion|arXiv (Cornell University)|2010. 01. 26.

Chronic Myeloid Leukemia Treatments참고 문헌 6인용 수 6

한 줄 요약

이 논문은 만성골수성백혈병을 모델링하는 지연미분방정식에서 평형점의 안정성 분석을 수정하고 보완한다. 르아프노프 방법과 특성방정식 분석을 통해 국소 안정성과 허프 분기 조건을 정확히 제시하며, 이전 연구에서 파arameter 의존성과 분기 임계값에 관해 오류가 있었던 점을 수정한다. 특히 양의 평형점에 대해 그러하다.

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는 평형점과 매개변수에 따라 달라진다.
• 안정성은 고유값의 실수부의 부호와 허프 분기의 첫 번째 르아프노프 계수를 중심으로 첫 번째 및 두 번째 르아프노프 방법을 사용하여 평가한다.
• 안정성 변화가 일어나는 임계 지연 r*를 구하기 위해 함수 g(r) = T^{−1}(−(δ+B)r) − arccos((δ+B)/(kB))를 유도하고 부호 변화를 분석한다.
• 고유값이 지연에 대해 미분계수가 0이 아니며 첫 번째 르아프노프 계수가 0이 아니라는 조건을 확인함으로써 허프 분기 조건을 설정한다.

실험 결과

연구 질문

• 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*에서 허프 분기가 발생하며, 안정성은 중심다양체 위에서 첫 번째 르아프노프 계수의 부호에 의해 결정된다.