[논문 리뷰] On global stability of the Lotka reactions with generalized mass-action kinetics

Balázs Boros, Josef Hofbauer외 1명|arXiv (Cornell University)|2016. 11. 17.

Gene Regulatory Network Analysis참고 문헌 18인용 수 15

한 줄 요약

이 논문은 임의의 거듭제곱법칙 동역학을 갖는 일반화된 Lotka-Volterra 시스템에서 유일한 양의 평형점의 전역 점근적 안정성을 연구한다. 정교하게 구성된 Dulac 함수를 사용한 Bendixson-Dulac 테스트를 적용하여, 주기적 궤도가 제거됨을 증명하고, 유계성 및 경계로의 수렴 방지를 함께 고려함으로써, 동역학 계수들이 야생행렬의 행렬식과 트레이스를 포함한 특정 부등식을 만족할 경우 전역 안정성이 확립된다.

ABSTRACT

Chemical reaction networks with generalized mass-action kinetics lead to power-law dynamical systems. As a simple example, we consider the Lotka reactions with two chemical species and arbitrary power-law kinetics. We study existence, uniqueness, and stability of the positive equilibrium, in particular, we characterize its global asymptotic stability in terms of the kinetic orders.

연구 동기 및 목표

• 일반화된 Lotka-Volterra 시스템에서 임의의 거듭제곱법칙 동역학을 갖는 유일한 양의 평형점의 전역 점근적 안정성을 분석하는 것.
• 비정수 동역학 계수를 갖는 시스템에 대해 이전의 국소 안정성 및 분기 분석을 전역 프레임워크로 확장하는 것.
• 비율 상수에 관계없이 양의 평형점이 전역 점근적으로 안정해지는 조건을 특징짓는 것.
• 주기적 해, 무한대로 발산하는 궤도, 양의 제1사분면 경계로의 수렴을 배제하는 것.

제안 방법

• 거듭제곱법칙 동역학에서 임의의 실수 지수를 갖는 일반화된 질량작용법칙 시스템으로 시스템을 수립한다.
• 분석을 단순화하기 위해 두 개의 지수를 줄이는 궤도 동치 시스템을 얻기 위해 변수를 변경한다.
• 형태 h(x,y) = x^{-p} y^{-q}인 Dulac 함수를 사용하여 Bendixson-Dulac 기준을 적용하여 주기적 궤도를 배제한다.
• 양의 제1사분면에서 그 발산이 엄격히 음수(또는 0)가 되도록 하는 Dulac 함수를 구성함으로써, 닫힌 궤도가 존재하지 않음을 보장한다.
• nullcline과 벡터장 방향 분석을 통한 비교를 통해 모든 해의 유계성을 증명한다.
• 궤도가 유한 시간 내에 nullcline을 가로질러 경계로의 수렴을 방지함으로써, 원점으로의 수렴 방지를 확립한다.

실험 결과

연구 질문

• RQ1 거듭제곱법칙 동역학을 갖는 일반화된 Lotka-Volterra 시스템에서, 어떤 조건에서 유일한 양의 평형점이 전역 점근적으로 안정해지는가?
• RQ2 임의의 실수 지수를 갖는 시스템에 대해, Bendixson-Dulac 기준을 적용하여 주기적 해를 배제할 수 있는가?
• RQ3 이러한 시스템 클래스에서 평형점에서의 야생행렬의 행렬식과 트레이스는 전역 안정성과 어떻게 관련되는가?
• RQ4 해가 양의 제1사분면의 경계로 수렴하거나 무한대가 되지 않도록 보장하기 위한 조건은 무엇인가?

주요 결과

• 행렬 C가 지수 차이의 행렬일 때, det C < 0 이고 야생행렬의 트레이스가 음수일 경우, 유일한 양의 평형점은 전역 점근적으로 안정해진다.
• 특화된 Dulac 함수 h(x,y) = x^{-p} y^{-q}를 사용한 Bendixson-Dulac 테스트는 (hf, hg)의 발산이 R²₊에서 엄격히 음수일 경우 주기적 해가 존재하지 않음을 보여준다.
• 벡터장의 방향이 nullcline에 대해 유의미하게 작용함으로써, 양의 제1사분면 내에서 해는 유계가 된다. 이는 무한대 성장 방지를 보장한다.
• 어느 해도 원점으로 수렴하지 않는다. x-nullcline 또는 y-nullcline를 가로질러 경로가 유한 시간 내에 원점 주변을 빠져나오게 된다.
• 트레이스가 0인 경우(Andronov-Hopf 분기), 안정성은 유도된 양수 d₁의 부호에 따라 달라진다: d₁ < 0 이면 평형점은 점근적으로 안정하다.
• 특수한 경우 α = 3/2, β = 1/2 에서는 평형점이 전역 점근적으로 안정하며, 이 경우 Dulac 함수는 평형점에서만 발산이 0이 된다.