[논문 리뷰] Turing and wave instabilities in hyperbolic reaction-diffusion systems: The role of second-order time derivatives and cross-diffusion terms on pattern formation

Joshua Ritchie, Andrew L. Krause외 1명|arXiv (Cornell University)|2022. 04. 28.

Nonlinear Dynamics and Pattern Formation참고 문헌 54인용 수 26

한 줄 요약

이 논문은 2차 시간 도함수와 교차확산 항을 포함한 쌍곡 반응-확산 시스템에서 튜링 및 파동 불안정성을 조사한다. 두 불안정성에 대한 必要 및 充분 조건을 유도하여, 시간 도함수가 튜링 패턴의 형태를 변화시키지 않고 그 존재 가능한 매개변수 영역만 제약한다는 것을 보여주며, 이는 고전적 시스템에 존재하지 않는 시공간 패턴인 파동 불안정성을 가능하게 한다. 핵심적으로, 파동 불안정성은 튜링 패턴과 상호 배타적인 조건에서 발생하므로, 활성화제가 억제제보다 빠르게 확산되는 경우에도 대칭성 붕괴가 가능하다.

ABSTRACT

Hyperbolic reaction-diffusion equations have recently attracted attention both for their application to a variety of biological and chemical phenomena, and for their distinct features in terms of propagation speed and novel instabilities not present in classical two-species reaction-diffusion systems. We explore the onset of diffusive instabilities and resulting pattern formation for such systems. Starting with a rather general formulation of the problem, we obtain necessary and sufficient conditions for the Turing and wave instabilities in such systems, thereby classifying parameter spaces for which these diffusive instabilities occur. We find that the additional temporal terms do not strongly modify the Turing patterns which form or parameters which admit them, but only their regions of existence. This is in contrast to the case of additional space derivatives, where past work has shown that resulting patterned structures are sensitive to second-order cross-diffusion and first-order advection. We also show that additional temporal terms are necessary for the emergence of spatiotemporal patterns under the wave instability. We find that such wave instabilities exist for parameters which are mutually exclusive to those parameters leading to stationary Turing patterns. This implies that wave instabilities may occur in cases where the activator diffuses faster than the inhibitor, leading to routes to spatial symmetry breaking in reaction-diffusion systems which are distinct from the well studied Turing case.

연구 동기 및 목표

• 일반적인 쌍곡 반응-확산 시스템에서 튜링 및 파동 불안정성이 발생하는 조건을 분류하는 것.
• 2차 시간 도함수와 교차확산 항이 유래 패턴의 실현 가능성과 구조에 미치는 영향을 규명하는 것.
• 패턴 형성 역학을 결정하는 시간적 고차 도함수와 공간적 고차 도함수의 역할을 대비하는 것.
• 파동 불안정성이 발생하는 매개변수 영역을 특정하여, 고전적 폴리노미얼 시스템에서는 불가능한 경우를 밝혀내는 것.
• 이러한 불안정성이 생물학적 및 물리적 패턴 형성에 미치는 영향을 분석하며, 특히 고전적 튜링 가정을 위반하는 경우(예: 활성화제의 빠른 확산)에 대한 의미를 탐색하는 것.

제안 방법

• 2차 시간 도함수와 교차확산 항을 포함한 일반적인 이종 쌍곡 반응-확산 시스템을 수립: τ₁∂²u/∂t² = ∇·(d₁₁∇u + d₁₂∇v) + f(∂u/∂t, ∂v/∂t, u, v), τ₂∂²v/∂t² = ∇·(d₂₁∇u + d₂₂∇v) + g(∂u/∂t, ∂v/∂t, u, v).
• 균일한 평형 상태에 대한 선형 안정성 분석을 적용하여 분산 관계 및 불안정성에 대한 고유값 조건을 도출한다.
• 튜링 불안정성(공간적으로 비균일한 정적 패턴)과 파동 불안정성(시공간적 진동 패턴)에 대한 必要 및 充분 조건을 유도한다.
• 카타네오-크리스토프 플럭 법칙과 반응-텔레그래프 방정식 체계를 활용하여 쌍곡적 구조가 물리적 원리로부터 타당함을 정당화한다.
• 완전 비선형 수치 시뮬레이션을 수행하여 이론적 예측을 검증하고 패턴 동역학을 시각화한다.
• 고전적 폴리노미얼 시스템, 교차확산이 있는/없는 쌍곡 시스템, 기저 상태(예: 스칼라 방정식) 등 다양한 경우를 비교한다.

실험 결과

연구 질문

• RQ1 2차 시간 도함수는 쌍곡 반응-확산 시스템에서 튜링 불안정성의 발생과 구조에 어떤 영향을 미치는가?
• RQ2 파동 불안정성이 쌍곡 시스템에서 어떤 조건에서 발생하는가? 그리고 고전적 튜링 또는 튜링-홉프 불안정성과 어떻게 다를까?
• RQ3 활성화제가 억제제보다 빠르게 확산되는 경우에도 파동 불안정성이 발생할 수 있는가? 이는 고전적 튜링 메커니즘을 넘어서 패턴 형성의 범위를 어떻게 확장하는가?
• RQ4 교차확산 항은 튜링 및 파동 패턴의 실현 가능성과 형태에 어떤 영향을 미치며, 이러한 항이 없는 시스템과 비교해 어떻게 다를까?
• RQ5 파동 불안정성을 통해 시공간 패턴을 생성하기 위한 최소한의 시스템 요구 조건(시간 도함수의 차수 및 종 수)은 무엇인가?

주요 결과

• 2차 시간 도함수는 표준 튜링 불안정성 조건 (3.15)–(3.16)를 변경하지 않으며, 오직 존재 가능한 매개변수 영역만 수정할 뿐이므로, 이는 때로는 이를 축소시킬 수 있다.
• 2차 시간 도함수를 추가하더라도 최종 정적 튜링 패턴의 형태는 본질적으로 변화하지 않으며, 이는 튜링 패턴 형태의 강건성을 시사한다.
• 파동 불안정성은 오직 2차 시간 도함수 덕분에 가능하며, 고전적 폴리노미얼 반응-확산 시스템에서는 발생하지 않는다.
• 파동 불안정성이 발생하는 매개변수 영역는 정적 튜링 패턴을 지지하는 영역와 상호 배타적이므로, 패턴 형성에 대해 서로 다른 경로를 제공할 수 있다.
• 파동 불안정성은 주로 한 가지 공간 파장에 의해 지배되는 시공간 진동 패턴을 유도하며, 이는 혼돈적인 튜링-홉프 행동과 뚜렷이 다른 동역학을 보인다.
• 교차확산 항은 튜링 및 파동 불안정성의 실현 가능성과 불안정성 조건을 크게 변화시키며, 이는 유래 패턴의 형태까지 변화시킬 수 있다. 이는 시간 도함수만으로는 달성할 수 없는 효과이다.