비국소 지연 반응확산 종양 모델의 동역학 분석 및 치료 효과

초록

본 논문은 비국소 지연을 포함한 반응‑확산 종양 모델에 치료 항목을 도입하고, Lyapunov‑Schmidt 축소를 이용해 비자명 정상상태의 존재와 근사식을 도출한다. 이후 선형화된 연산자의 스펙트럼을 분석하여 안정성 조건을 제시하고, 지연에 의존하는 Hopf 분기 구간을 규명한다. 마지막으로 구체적인 예시와 수치 시뮬레이션을 통해 치료 파라미터 변화가 시스템의 안정성과 주기적 진동에 미치는 영향을 확인한다.

상세 분석

이 연구는 종양 성장과 치료 효과를 공간적으로 이질적인 환경에서 기술하기 위해, 확산 항(d Δu)과 비국소 지연 항 ∫ΩS(x,y)u(y,t‑τ)dy를 포함한 비선형 반응‑확산 방정식을 설정한다. 치료 항 R(x,t)=β q(u)+r(x) 로 모델링하여, β는 전반적인 치료 강도, q(u)는 세포 밀도에 대한 억제 효과, r(x)는 공간적 치료 효율 차이를 반영한다.
주요 수학적 기법은 다음과 같다. 첫째, β를 분기 매개변수로 잡고, β=β에서의 주특잇값 문제 d Δφ+F(0,0)φ−r(x)φ=βφ 를 통해 양의 주특잇값 β와 양의 고유함수 φ를 확보한다. β에서 연산자 Lβ는 비가역이므로 직접적인 암시적 함수 정리를 적용할 수 없으며, 대신 Lyapunov‑Schmidt 축소를 수행한다. X=Ker(Lβ*)⊕X₁, Y=Ker(Lβ*)⊕Y₁ 로 분해하고, 투영 연산자 P와 I‑P를 이용해 무한 차원 문제를 1차원 스칼라 방정식 f(ν,β)=0 로 환원한다. 여기서 ν는 φ에 대한 스칼라 계수이며, f는 적절히 전개된 형태로 ϱ(β)+κ(β)ν+o(ν) 로 표현된다. κ(β)≠0이면 암시적 함수 정리를 적용해 ν=ν(β) 를 얻고, uβ(x)=ν(β)φ*(x)+g(ν(β)φ*,β) 로 근사 정상상태를 구성한다. 특히, 2차 항까지 전개하면 uβ(x)≈θ₄φ*(x)κ(β*)(β‑β*)+m(x)(β‑β*)² 형태가 도출된다. 이 식을 통해 정상상태의 부호와 존재 구간을 β와 κ(β*)의 부호 관계로 명시한다.
두 번째 단계는 정상상태의 선형 안정성 분석이다. 정상상태 uβ(x) 를 기준으로 선형화하면 지연 항이 e^{‑λτ} 로 나타나는 특성 방정식 Π(τ,β,λ)φ=0 를 얻는다. 여기서 Π는 확산, 치료, 비국소 상호작용을 모두 포함한다. β>β* 일 때는 0 상태가 안정하지만, 양의 정상상태에 대해서는 λ의 실수부가 음이 되도록 하는 충분조건을 κ(β)와 지연 τ의 조합으로 제시한다. 특히, Lemma 3.1 에서는 λ가 실수부 비음이 되면 κ(β*)·(λβ−β*) 가 유계임을 보이며, 이는 Hopf 분기 발생 가능성을 제한한다.
세 번째 단계는 Hopf 분기 조건이다. 특성 방정식을 λ=iω (ω>0) 로 두고 실부와 허부를 분리하면, ω와 τ 사이에 파라미터 곡선 τ=τ₀(β) 가 존재함을 확인한다. 이 곡선 위에서 복소 고유값 쌍이 실축을 가로지며, 서브크리티컬 Hopf 분기가 일어난다. 논문은 κ(β*)·q′(0)·θ₃ 와 같은 계수들의 부호가 분기 방향을 결정한다는 점을 강조한다.
마지막으로, 저자들은 구체적인 1‑차원 구간 Ω=