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

Resumo Técnico: Sobre a eficiência do controle Hamiltoniano de pares para dessincronizar o modelo de Kuramoto de ordem superior

Enunciado do Problema
A sincronização em sistemas de osciladores acoplados é um fenômeno ubíquo com implicações duais: é benéfica em contextos como redes elétricas, mas prejudicial em outros, como crises epilépticas. Enquanto estratégias de controle tradicionais frequentemente dependem de interações de pares, evidências recentes indicam que muitos sistemas naturais e de engenharia exibem interações de ordem superior (de muitos corpos). Essas interações alteram fundamentalmente a dinâmica do sistema, notadamente ao aumentar a estabilidade linear dos estados sincronizados enquanto, simultaneamente, reduzem seus bacias de atração — um fenômeno descrito como bacias "mais profundas, porém menores".

O problema central abordado é como esses efeitos conflitantes das interações de ordem superior influenciam a eficiência de uma estratégia de controle de pares minimamente invasiva projetada para dessincronizar tais sistemas. Especificamente, os autores investigam se o aumento da estabilidade linear (que dificulta a dessincronização) ou a redução do tamanho da bacia (que facilita a dessincronização) domina quando um esquema de controle originalmente desenvolvido para redes de pares é aplicado a modelos de Kuramoto de ordem superior.

Metodologia
O estudo emprega um modelo de Kuramoto de ordem superior composto por $N$ osciladores de fase com interações de pares ($K_1$) e de três corpos ($K_2$). O sistema é governado por:
$$\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$$
onde $\mu$ é a força de controle e $\delta_i$ indica se o nó $i$ é controlado.

O mecanismo de controle é um controle Hamiltoniano de pares minimamente invasivo, adaptado de Asllani et al. [35]. O termo de controle $p_i$ é definido como:
$$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})$$
Este controle atua apenas em um subconjunto de $M$ nós fixos. Ele mede o nível de sincronização dos nós controlados e aplica uma força proporcional à magnitude da sincronização, que desaparece quando o sistema está totalmente dessincronizado. Esta abordagem evita a complexidade computacional de derivar termos de controle de ordem superior.

Os autores avaliam a eficiência do controle utilizando dois parâmetros de ordem:

1. Sincronização de fase ($\hat{R}$): O parâço de ordem de Kuramoto com média temporal.
2. Sincronização de frequência ($\hat{R}_{\dot{\theta}}$): Uma medida de travamento de fase (sincronização de frequência).

Experimentos numéricos foram conduzidos em duas topologias:

• Hiperanéis (Hyperrings): Estruturas regulares que suportam estados torcidos (twisted states).
• Hipergrafos aleatórios de Erdős-Rényi: Estruturas irregulares sem invariância rotacional.

O estudo variou sistematicamente a força de acoplamento de ordem superior ($K_2$), a força de controle ($\mu$), o número de nós controlados ($M$) e a distância inicial ($\epsilon$) do estado do sistema em relação ao centro da bacia de sincronização.

Resultados Principais

1. Condições Iniciais Sincronizadas:
   Quando o sistema parte próximo ao estado totalmente sincronizado (centro da bacia de atração), o aumento da força de acoplamento de ordem superior $K_2$ torna a dessincronização mais difícil. A força de controle crítica ($\mu_c$) necessária para atingir a dessincronização aumenta monotonicamente com $K_2$. Isso indica que, neste regime, o aumento da estabilidade linear induzido pelas interações de ordem superior domina sobre a redução da bacia de atração. O controle de pares permanece eficaz se $\mu$ for suficientemente grande, mas o custo do controle aumenta com $K_2$.

2. Condições Iniciais Perturbadas:
   Quando o sistema parte de um ponto mais distante do estado sincronizado (maior $\epsilon$), emerge um comportamento não monotônico:

   • $K_2$ Intermediário: A dessincronização é mais difícil; o controle é menos eficiente.
   • $K_2$ Grande: A dessincronização torna-se mais fácil; a eficiência do controle melhora.
     Este comportamento sugere que, para $K_2$ grande, o encolhimento da bacia de atração torna-se o fator dominante, permitindo que o controle empurre as trajetórias para fora da bacia mais facilmente, apesar da maior estabilidade linear.

3. Controle Parcial:
   O estudo confirma que a dessincronização pode ser alcançada mesmo quando apenas uma fração de nós ($M < N$) é controlada. O número mínimo de nós necessários ($M_c$) geralmente aumenta com $K_2$ para condições iniciais sincronizadas, mas exibe comportamento não monotônico para condições iniciais perturbadas, espelhando as tendências observadas com a força de controle.

4. Dessincronização de Frequência vs. Fase:
   O controle suprime efetivamente a sincronização de fase. No entanto, a sincronização de frequência (por exemplo, em estados torcidos ou de agrupamento) é mais robusta, particularmente para $K_2$ moderado em hiperanéis, devido às simetrias estruturais e ao fato de o controle não visar diretamente o travamento de frequência.

Significância e Alegações
O artigo afirma fornecer uma visão mais clara da interação entre interações de ordem superior e eficiência de controle. A principal contribuição é demonstrar que o efeito "mais profundo, porém menor" das interações de ordem superior leva a regimes de controle distintos dependendo do estado inicial:

• Próximo ao estado sincronizado, as interações de ordem superior atuam como um obstáculo ao controle devido à estabilidade aumentada.
• Longe do estado sincronizado, interações de ordem superior suficientemente fortes podem, de fato, facilitar o controle ao encolher a bacia de atração.

Os autores afirmam que uma estratégia de controle de pares, minimamente invasiva e simples, é suficiente para dessincronizar sistemas de Kuramoto de ordem superior, desde que a intensidade do controle ou o número de nós controlados seja adequado. Este achado é significativo porque sugere que termos de controle de ordem superior complexos e computacionalmente caros (como proposto na Ref. [100]) podem não ser estritamente necessários para a dessincronização, oferecendo um caminho de implementação mais prático. Os resultados destacam a importância de considerar tanto a estabilidade linear quanto a geometria da bacia ao projetar estratégias de controle para redes complexas com interações de ordem superior.