The role of asymmetric time delay and its structure in 1D swarmalators

This paper investigates a one-dimensional swarmalator model with asymmetric time delay, revealing that the delay's internal structure fundamentally reshapes the collective phase diagram by systematically expanding the active $\pi$ state and establishing that the delay's form, rather than just its magnitude, is a decisive factor in emergent swarmalator behavior.

The Gist

Imagine a massive crowd of tiny robots, or perhaps a school of fish, all moving around and trying to coordinate their actions. In the world of physics, these are called swarmalators. They are special because they do two things at once: they move together in space (like a flock of birds) and they sync up their internal "beats" or rhythms (like a group of people clapping in unison).

Usually, when these groups interact, they assume everyone hears and reacts to each other instantly. But in the real world, nothing is instant. There is always a tiny delay—like the time it takes for a sound to travel or for a brain to process a signal. This paper asks a simple but crucial question: Does it matter how that delay happens?

The "Asymmetric" Twist

Most previous studies assumed the delay was "symmetric," meaning every part of the interaction was slowed down equally. It's like if everyone in a choir had to wait the same amount of time before singing their next note.

The authors of this paper decided to test an asymmetric delay. Imagine a scenario where:

• The robots wait to react to where their neighbors are (spatial delay).
• But they react instantly to what rhythm their neighbors are keeping (phase instant).
• Or vice versa.

They found that where you put the delay changes everything. It's not just about how long the wait is, but which part of the conversation is delayed.

The Five Ways the Crowd Behaves

By running computer simulations, the researchers discovered that these delayed swarmalators settle into five distinct "moods" or states:

1. The Async State (The Chaos): Everyone is doing their own thing. They are scattered randomly in space and out of sync in rhythm. It's a noisy, disorganized crowd.
2. The Static $\pi$ State (The Frozen Pairs): The crowd splits into two perfect groups. One group is at one spot, the other is exactly opposite. They are frozen in place, perfectly synchronized, but separated by a "half-turn" (180 degrees).
3. The Active $\pi$ State (The Dancing Pairs): This is the most interesting new discovery. Like the frozen pairs, they split into two groups. But instead of sitting still, they start marching in a circle together. The delay actually creates this motion. Without the delay, they would just sit still.
4. The Phase Wave (The Wave): Imagine a stadium wave. The crowd moves in a line where position and rhythm are perfectly linked. If you are at a certain spot, you are at a certain point in your rhythm.
5. The Unsteady State (The Jitter): The crowd can't decide. They oscillate back and forth between order and chaos, never settling down.

The Big Discovery: Delay as a "Volume Knob"

The most surprising finding is how the delay acts like a control knob for the crowd's behavior:

• When the delay is in the "movement" part (the sine terms): Increasing the delay acts like a magnet for the Active $\pi$ state. The longer the delay, the more likely the crowd is to split into two groups and start marching in circles. The delay stabilizes this dancing motion.
• When the delay is in the "rhythm" part (the cosine terms): The crowd tends to become unsteady and jittery. They can't find a stable rhythm and start shaking or oscillating wildly.

The "Symmetric" vs. "Asymmetric" Showdown

The authors compared their new "asymmetric" model to the old "symmetric" model (where delay hits everything equally).

• Symmetric Delay: Tends to make the crowd jittery and unsteady. It's hard for them to find a stable rhythm.
• Asymmetric Delay: Can actually help the crowd find a stable, organized way to move (the Active $\pi$ state).

The Takeaway

Think of a group of dancers trying to learn a routine.

• If they all have a slow, identical reaction time, they might just stumble and jitter.
• But if the delay is structured specifically—say, they wait a moment to see where their partner is, but react instantly to the music—the delay actually helps them lock into a beautiful, synchronized dance where they spin in pairs.

The paper concludes that the structure of the delay is just as important as the delay itself. It's not just about how slow the signal is; it's about which signal is slow. This changes the entire "map" of how these groups behave, showing that a little bit of the right kind of lag can turn chaos into a coordinated dance.

Technical Summary

Technical Summary: The Role of Asymmetric Time Delay and its Structure in 1D Swarmalators

Problem Statement
Swarmalators are coupled oscillators that synchronize simultaneously in space and phase, modeling systems ranging from biological microswimmers to robotic swarms. While time delay is ubiquitous in such systems due to finite signal propagation and processing lags, previous research has largely focused on symmetric delays, where the same delay $\tau$ affects all coupling channels equally. The structural form of the delay—specifically whether it is symmetric or asymmetric (affecting only specific terms in the interaction)—has received little attention. This paper addresses the gap in understanding how the structure of asymmetric time delay influences the emergent collective behavior of swarmalator populations, distinct from the mere magnitude of the delay.

Methodology
The authors investigate a one-dimensional (1D) swarmalator model where agents possess a spatial position $x_i$ and an internal phase $\theta_i$, both defined on a circle $S^1$. The study focuses on a specific variant of asymmetric delay where the time lag $\tau$ enters only the self-interaction (sine) terms of the spatial and phase dynamics, while the cross-modal (cosine) coupling factors remain instantaneous. The governing equations are delay differential equations:
$$\dot{x}_i = \frac{J}{N} \sum_{j=1}^N \sin(x_j(t-\tau) - x_i(t)) \cos(\theta_j(t) - \theta_i(t))$$
$$\dot{\theta}_i = \frac{K}{N} \sum_{j=1}^N \sin(\theta_j(t-\tau) - \theta_i(t)) \cos(x_j(t) - x_i(t))$$
where $J$ and $K$ are coupling constants.

The methodology combines:

1. Numerical Simulations: Direct integration of the governing equations for $N=1000$ agents to identify collective states and construct bifurcation diagrams in the $(J, K)$ parameter plane for varying delays ($\tau$).
2. Analytical Stability Analysis: Linear stability analysis of key collective states (async, phase wave, and $\pi$ states) in the continuum limit. The authors derive characteristic equations using Fourier mode expansions and the Lambert $W$ function to determine closed-form stability boundaries.
3. Comparative Analysis: The asymmetric delay model is compared against a complementary asymmetric case (delay in cosine terms only) and a symmetric delay model (delay in all terms) to isolate the effects of delay structure.

Key Results
The study identifies five distinct collective states:

• Static $\pi$ state: Two clusters separated by spatial and phase distance $\pi$, with zero velocity.
• Active $\pi$ state: Similar two-cluster structure but with non-zero collective velocity.
• Async state: Uniform distribution in space and phase with no order.
• Phase wave state: Strong correlation between position and phase ($x = \pm \theta + c$).
• Unsteady state: Persistent oscillations of order parameters.

The primary findings regarding the impact of asymmetric delay structure are:

• Expansion of the Active $\pi$ State: In the model where delay enters the sine terms, increasing the delay $\tau$ systematically expands the stability region of the active $\pi$ state at the expense of other ordered states. This contrasts with the symmetric delay model, which more strongly promotes unsteady, non-stationary dynamics.
• Stability Boundaries: The authors derive analytical conditions for the stability of the async state, phase wave, and $\pi$ states. For the async state, the stability boundary is determined by the condition $2J + 2K + \tau JK = 0$. For the phase wave, stability depends on the interplay between coupling constants and delay, with Hopf bifurcations occurring under specific conditions.
• Bistability: The analysis reveals regions of bistability where the active $\pi$ state and the async state coexist as stable attractors. The asymptotic state selected by the system depends on the initial conditions, with the basins of attraction shifting as coupling parameters vary.
• Structural Sensitivity: The paper demonstrates that the internal structure of the delay is a decisive factor. When delay is introduced into the "sine" terms, it promotes ordered active states (active $\pi$). Conversely, when delay is introduced into the "cosine" terms (complementary case), it favors the emergence of unsteady or non-stationary states.

Significance and Claims
The paper claims that the internal structure of time delay, not merely its magnitude, fundamentally reshapes the collective phase diagram of swarmalator populations. By providing closed-form stability conditions validated against numerical simulations, the work establishes that asymmetric delay can act as a control parameter to promote specific ordered states (like the active $\pi$ state) that are not observed or are less stable in symmetric delay models or delay-free systems. The authors suggest these findings are relevant for understanding biological systems (e.g., predator-prey interactions, motor systems) and engineered swarms where signal processing and spatial alignment may operate on different timescales. The study serves as a foundational step in isolating the essential effects of delay structure in lower-dimensional variants before tackling more complex higher-dimensional settings.