On the efficiency of pairwise Hamiltonian control to desynchronize the higher-order Kuramoto model

기술 요약: 고차 쿠라모토 모델의 탈동기화를 위한 쌍체 해밀토니안 제어의 효율성에 관하여

문제 정의
결합된 진동자 시스템에서의 동기화는 편무적인 현상으로, 전력망과 같은 맥락에서는 유익하지만, 간질 발작과 같은 다른 맥락에서는 해로울 수 있는 이중적인 함의를 갖는다. 전통적인 제어 전략은 흔히 쌍체 상호작용(pairwise interactions)에 의존하지만, 최근의 증거들은 많은 자연 및 공학 시스템이 고차(다체) 상호작용을 나타낸다는 것을 보여준다. 이러한 상호작용은 시스템 역학을 근본적으로 변화시키며, 특히 동기화된 상태의 선형 안정성을 증가시키는 동시에 그 인력 영역(basin of attraction)을 축소시킨다. 이는 "더 깊지만 더 작은" 영역이라는 현상으로 설명된다.

본 연구가 다루는 핵심 문제는 이러한 고차 상호작용의 상충하는 효과들이, 고차 쿠라모토 모델을 대상으로 설계된 최소 침습적 쌍체 제어 전략의 효율성에 어떻게 영향을 미치는가이다. 구체적으로, 저자들은 쌍체 네트워크를 위해 개발된 제어 기법을 고차 쿠라모토 모델에 적용할 때, 증가된 선형 안정성(탈동기화를 방해함)과 축소된 인력 영역(탈동기화를 용이하게 함) 중 어느 것이 지배적인지를 조사한다.

방법론
본 연구는 쌍체 상호작용($K_1$)과 3체 상호작용($K_2$)을 모두 포함하는 $N$개의 위상 진동자로 구성된 고차 쿠라모토 모델을 사용한다. 시스템은 다음과 같이 제어된다:
$$\dot{\theta}_i = \omega_i + \frac{K_1}{\langle k^{(1)} \rangle} \sum_{j} A^{(1)}_{ij} \sin(\theta_j - \theta_i) + \frac{K_2}{\langle k^{(2)} \rangle} \sum_{j,k} A^{(2)}_{ijk} \sin(\theta_j + \theta_k - 2\theta_i) + \mu \delta_i p_i$$
여기서 $\mu$는 제어 강도이며, $\delta_i$는 노드 $i$가 제어되는지 여부를 나타낸다.

제어 메커니즘은 Asllani 등 [35]으로부터 변형된 최소 침습적 쌍체 해밀토니안 제어이다. 제어항 $p_i$는 다음과 같이 정의된다:
$$p_i = -\frac{1}{2K_1^2} \frac{M^2}{\langle k^{(1)} \rangle^2} R_M \tilde{R}_{M,i} \cos(\Psi_M - \tilde{\Psi}_{M,i})$$
이 제어는 고정된 $M$개의 노드 부분집합에만 작용한다. 이는 제어되는 노드들의 동기화 수준을 측정하고, 동기화 크기에 비례하는 힘을 가하며, 시스템이 완전히 탈동기화되면 0이 된다. 이 접근 방식은 고차 제어 항을 도출하는 계산 복잡성을 피한다.

저자들은 두 가지 질서 매개변수를 사용하여 제어 효율성을 평가한다:

1. 위상 동기화 ($\hat{R}$): 시간 평균 쿠라모토 질서 매개변수.
2. 주파수 동기화 ($\hat{R}_{\dot{\theta}}$): 위상 잠금(주파수 동기화)의 척도.

수치 실험은 다음 두 가지 토폴로지에서 수행되었다:

• 하이퍼링(Hyperrings): 꼬임 상태(twisted states)를 지원하는 규칙적인 구조.
• 무작위 에르되시-레니 하이퍼그래프(Random Erdős-Rényi hypergraphs): 회전 불변성이 없는 불규칙한 구조.

연구는 고차 결합 강도($K_2$), 제어 강도($\mu$), 제어되는 노드의 수($M$), 그리고 시스템 상태와 동기화 영역 중심 사이의 초기 거리($\epsilon$)를 체계적으로 변화시켰다.

주요 결과

1. 동기화된 초기 조건:
   시스템이 동기화 인력 영역의 중심(완전 동기화 상태) 근처에서 시작될 때, 고차 결합 강도 $K_2$를 증가시키면 탈동기화가 더 어려워진다. 임계 제어 강도($\mu_c$)는 $K_2$에 따라 단조 증가한다. 이는 이 영역에서 고차 상호작용에 의해 유도된 선형 안정성 증가가 인력 영역의 축소보다 지배적임을 나타낸다. 쌍체 제어는 $\mu$가 충분히 크다면 여전히 효과적이지만, 제어 비용은 $K_2$와 함께 상승한다.

2. 섭동된 초기 조건:
   시스템이 동기화 상태에서 더 멀리 떨어져 있을 때(더 큰 $\epsilon$), 비단조적 행동이 나타난다.

   • 중간 단계의 $K_2$: 탈동기화가 가장 어려우며, 제어 효율이 가장 낮다.
   • 큰 $K_2$: 탈동기화가 더 쉬워지며, 제어 효율이 향상된다.
     이러한 행동은 큰 $K_2$의 경우, 인력 영역의 축소가 지배적인 요인이 되어, 높은 선형 안정성에도 불구하고 제어가 궤적을 영역 밖으로 밀어내는 것을 더 쉽게 만든다는 것을 시사한다.

3. 부분 제어:
   본 연구는 오직 일부 노드($M < N$)만을 제어하더라도 탈동기화가 달성될 수 있음을 확인하였다. 최소 필요 노드 수($M_c$)는 동기화된 초기 조건에서는 일반적으로 $K_2$에 따라 증가하지만, 섭동된 초기 조건에서는 제어 강도에서 관찰된 경향과 유사하게 비단조적 행동을 보인다.

4. 주파수 vs 위상 탈동기화:
   제어는 위상 동기화를 효과적으로 억제한다. 그러나 주파수 동기화(예: 하이퍼링에서의 꼬임 또는 클러스터 상태)는 구조적 대칭성 및 제어가 주파수 잠금을 직접 겨냥하지 않는다는 점 때문에 특히 중간 정도의 $K_2$에서 더 견고하게 유지된다.

의의 및 주장
본 논문은 고차 상호작용과 제어 효율성 사이의 상호작용에 대한 명확한 그림을 제공한다고 주장한다. 주요 기여는 고차 상호작용의 "더 깊지만 더 작은" 효과가 초기 상태에 따라 서로 다른 제어 영역을 생성한다는 것을 입증한 것이다:

• 동기화 상태 근처에서는 고차 상호작용이 안정성 강화로 인해 제어의 장애물로 작용한다.
• 동기화 상태에서 멀어질 경우, 충분히 강한 고차 상호작용은 인력 영역을 축소함으로써 오히려 제어를 용이하게 할 수 있다.

저자들은 단순하고 최소 침설적인 쌍체 제어 전략이, 제어 강도나 제어 노드 수가 적절하다면 고차 쿠라모토 시스템을 탈동기화하기에 충분하다는 것을 주장한다. 이 발견은 고차 제어 항(Ref. [100]에서 제안된 것과 같은)을 사용하는 것이 반드시 필수적이지는 않다는 것을 시사하며, 이는 더 실용적인 구현 경로를 제공한다. 결과적으로 본 연구는 고차 상호작용을 가진 복잡한 네트워크를 위한 제어 전략을 설계할 때 선형 안정성과 인력 영역의 기하학적 구조를 모두 고려하는 것이 중요함을 강조한다.