[논문 리뷰] An efficient method to solve the mathematical model of HIV infection for CD8+ T-cells

Samad Noeiaghdam, Emran Khoshrouye Ghiasi|arXiv (Cornell University)|2019. 06. 28.

Fractional Differential Equations Solutions참고 문헌 52인용 수 7

한 줄 요약

이 논문은 CD8+ T세포를 포함하는 비선형 수학적 모델인 HIV 감염 모델을 해결하기 위해 호모토피 분석 변환 방법(HATM)을 제안한다. 호모토피 분석 방법과 라플라스 변환을 결합함으로써, 수렴 분석, $\bar{h}$-곡선 및 잔차 오차도를 통해 $N=5$ 및 $N=10$에 대해 고정밀도 해를 도출하여 복잡한 생물학적 시스템에 대한 이 방법의 정확성과 강건성을 확인한다.

ABSTRACT

In this paper, the mathematical model of HIV infection for CD8+ T-cells is illustrated. The homotopy analysis method and the Laplace transformations are combined for solving this model. Also, the convergence theorem is proved to demonstrate the abilities of presented method for solving non-linear mathematical models. The numerical results for N = 5, 10 are presented. Several h-curves are plotted to show the convergence regions of solutions. The plots of residual error functions indicate the precision of presented method.

연구 동기 및 목표

• CD8+ T세포를 포함하는 HIV 감염 동역학을 모델링하는 비선형 상미분 방정식계를 효율적으로 해결하기 위한 반해석적 방법을 개발하는 것.
• 호모토피 분석 방법(HAM)과 라플라스 변환을 융합하여 수렴성과 정확도를 향상시킨 새로운 HATM 프레임워크를 구축하는 것.
• 주어진 HIV-CD8+ T세포 모델에 대해 HATM 해의 수렴성을 엄밀한 수렴 정리로 증명하는 것.
• 다양한 근사 차수($N=5$, $N=10$)에 대해 잔차 오차 함수 및 $ar{h}$-곡선 분 析를 통해 방법의 정밀도를 검증하는 것.

제안 방법

• HATM는 HIV-CD8+ T세포 모델(1)의 각 방정식에 라플라스 변환을 적용함으로써 적용되며, 이는 시스템을 대수방정식 집합으로 변환한다.
• 보조 매개변수 $ \bar{h} $를 사용하여 호모토피를 구성하고, 해는 임bedding 매개변수 $ q $의 거듭제곱에 대한 급수 전개로 표현되며, $ n $차 수준의 변형 방정식들이 유도된다.
• 역 라플라스 변환을 $ n $차 수준의 변형 방정식에 적용하여 $ T_n(t), I_n(t), V_n(t), Z_n(t), Z_{a,n}(t) $에 대한 재귀적 표현을 도출한다.
• 근사 해는 $ T_N(t) = \sum_{n=0}^{N} T_n(t) $로 잘라낸 급수 형태로 얻어지며, 다른 변수들에 대해서도 유사한 형태를 가진다.
• 수렴 제어 매개변수 $ \bar{h} $는 $ \bar{h} $-곡선을 그려서 수렴 영역을 식별함으로써 최적화된다.
• 잔차 오차 함수 $ E_{N,T}(t), E_{N,I}(t), \dots $가 계산되고 플롯되어 해의 정확도를 정량화한다.

실험 결과

연구 질문

• RQ1 HATM는 CD8+ T세포를 포함하는 HIV 감염 동역학을 모델링하는 비선형 상미분 방정식계를 효과적으로 해결할 수 있는가?
• RQ2 HIV-CD8+ T세포 모델에 대한 HATM 해의 수렴 영역은 무엇이며, 근사 차수 $ N $에 따라 어떻게 변화하는가?
• RQ3 잔차 오차 함수로 측정했을 때, HATM 해는 정확한 해와 비교해 얼마나 정밀한가?
• RQ4 수렴 제어 매개변수 $ \bar{h} $의 다양한 값은 HATM 해의 수렴성과 안정성에 어떻게 영향을 미치는가?
• RQ5 HATM는 HIV-CD8+ T세포 시스템과 같은 복잡한 비선형 생물수학적 모델을 해결하기 위한 신뢰할 수 있고 융통성 있는 프레임워크를 제공하는가?

주요 결과

• N=5에 대해 HATM 해는 $ t $, $ \bar{h} $ 및 그 거듭제곱에 대한 급수 형태로 명시적으로 유도되었으며, 계수는 모델의 매개변수를 포함한다.
• N=5일 때 해의 수렴 영역은 $ -0.9 \leq \bar{h}_T \leq -0.2 $, $ -0.8 \leq \bar{h}_I \leq -0.2 $, $ -0.8 \leq \bar{h}_V \leq -0.2 $, $ -1.2 \leq \bar{h}_Z \leq -0.6 $, 및 $ -0.8 \leq \bar{h}_{Z_a} \leq -0.4 $로 확인되었다.
• N=10일 때 수렴 영역은 $ -0.9 \leq \bar{h}_T \leq -0.4 $, $ -1 \leq \bar{h}_I \leq -0.2 $, $ -1 \leq \bar{h}_V \leq -0.3 $, $ -1.2 \leq \bar{h}_Z \leq -0.4 $, 및 $ -1 \leq \bar{h}_{Z_a} \leq -0.4 $로 확장되었으며, 이는 고차수 근사에서 안정성이 향상됨을 시사한다.
• $ \bar{h} = -0.8 $에서 $ N=5 $ 및 $ N=10 $에 대해 그린 잔차 오차 함수는 시간 영역 전반에 걸쳐 거의 영에 가까운 크기를 보이며, 높은 해의 정확성을 확인한다.
• 수렴 정리는 엄밀하게 증명되었으며, 이는 HATM가 HIV-CD8+ T세포 모델과 같은 비선형 시스템을 해결하는 데 신뢰할 수 있는 방법임을 입증한다.
• 이 방법은 강력한 유연성과 정밀도를 보이며, $ N=10 $일 때 잔차 오차가 항상 $ 10^{-3} $ 이하로 유지되어 고도로 신뢰할 수 있는 계산 성능을 나타낸다.