# Experiments with arbitrary networks in time-multiplexed delay systems

###### Abstract

We report a new experimental approach using an optoelectronic feedback loop to investigate the dynamics of oscillators coupled on large complex networks with arbitrary topology. Our implementation is based on a single optoelectronic feedback loop with time delays. We use the space-time interpretation of systems with time delay to create large networks of coupled maps. Others have performed similar experiments using high-pass filters to implement the coupling; this restricts the network topology to the coupling of only a few nearest neighbors. In our experiment, the time delays and coupling are implemented on a field-programmable gate array, allowing the creation of networks with arbitrary coupling topology. This system has many advantages: the network nodes are truly identical, the network is easily reconfigurable, and the network dynamics occur at high speeds. We use this system to study cluster synchronization and chimera states in both small and large networks of different topologies.

Complex networks of coupled oscillators have proven to be systems that can display incredibly rich dynamical behaviors. Despite great theoretical advances in our understanding of coupled oscillator networks, it has proven difficult to design experiments that permit the study of dynamics on large networks with arbitrary topology. Here we present a new experimental approach that allows for the investigation of large networks of truly identical nodes with arbitrary topology. Our approach relies upon the space-time interpretation of systems with time delay in order to construct a network of coupled maps using a single nonlinear, time-delayed feedback loop. Our experiment provides a significant improvement over previous network experiments in terms of speed, ease of implementation, and cost. This paper describes our experimental approach and its capabilities, including the ability to study cluster synchronization and chimera states in both small and large networks.

## I Introduction

Systems of coupled oscillators have been of great theoretical and practical interest since Huygens’s observation of the anti-phase synchronization of two pendulum clocks Huygens (1897). Over the last 30 years, our understanding of coupled oscillators has grown from the two coupled pendula of Huygens to the synchronization Pikovsky, Rosenblum, and Kurths (2003); Pecora and Carroll (2015) and control Sorrentino et al. (2007) of nonlinear oscillators coupled together in complex networks of interactions. Recently, there has been an explosion of interest in patterns of network synchronization, such as chimera states Abrams and Strogatz (2004); Panaggio and Abrams (2015) and cluster synchronization Zhou and Kurths (2006); Ji et al. (2013)(see, e.g., the Chaos special issue on Patterns of Network Synchronization Abrams, Pecora, and Motter (2016)).

However, experimental progress on complex networks has been slow compared to theoretical advances. For example, it took a decade from the initial prediction of chimera states in 2002 Kuramoto and Battogtokh (2002) to their first experimental observations in 2012 Hagerstrom et al. (2012); Tinsley, Nkomo, and Showalter (2012). Only recently has there been significant progress in the experimental study of large networks.

Many of the experiments on large networks of oscillators are limited in the types of coupling they can achieve. Electrochemical Schmidt et al. (2014) and passively mode locked laser Viktorov et al. (2014) systems with nonlinear global coupling have been used to study chimera states. Electronic and optoelectronic systems with time delay have used the virtual space-time Arecchi et al. (1992) interpretation to implement large networks of coupled oscillators for the study of both chimera states Larger, Penkovsky, and Maistrenko (2013, 2015) and reservoir computing Appeltant et al. (2011); Paquot et al. (2012); Larger et al. (2012); Antonik, Haelterman, and Massar (2017); however, these systems so far have been restricted to networks with cyclic symmetry. Coupled metronomes have been used to study chimera states in hierarchical networks Martens et al. (2013); Blaha et al. (2016), but this approach is impractical for the implementation of networks with complex topologies. The authors are aware of only a few experiments that can be used to study large networks with arbitrary coupling topology. Chemical oscillators Tinsley, Nkomo, and Showalter (2012); Totz et al. (2015); Wickramasinghe and Kiss (2013, 2014) and spatial light modulator feedback loops Hagerstrom et al. (2012); Pecora et al. (2014) have been used to study hundreds of coupled oscillators, but these experiments run on slow timescales, can be expensive and have inherent heterogeneities between nodes.

In this paper, we describe a new approach that allows for the experimental investigation of the dynamics of arbitrary networks of coupled maps. Our approach is similar to that of previous opto-electronic experiments Kouomou et al. (2005) that use the space-time interpretation to create coupled node networks with a single time-delayed system Appeltant et al. (2011); Larger et al. (2012); Larger, Penkovsky, and Maistrenko (2013, 2015); Antonik, Haelterman, and Massar (2017). However, these experiments are restricted to cyclically symmetric networks due to the time-invariant (in these cases, highpass) filtering in their electronics. We use optics and electronics designed to have no highpass filtering; this removes the few-nearest neighbor coupling that was always present in previous work. Additionally, we use a field-programmable gate array (FPGA) to perform time multiplexing in order to implement arbitrary coupling topologies between the nodes. This approach has the additional advantage that it makes it possible to obtain a network of truly identical nodes because all nodes utilize the same electronic and optical components. This eliminates the (impossible) requirement to build many individual identical systems, and replaces it with the much easier task of obtaining temporal stability. Our experimental system is low cost, high speed, and easily reconfigurable compared with other experiments on large networks with arbitrary topology.

Our system provides a versatile experimental platform to study dynamics on all sorts of networks. The coupling between nodes can be instantaneous and/or time-delayed. The network can be static, time-varying, or adaptive. The nodes are truly identical since they all use the same experimental apparatus; however, they can be made inhomogeneous by controlling node parameters in the FPGA for each node. The ability to control the inhomogeneity of the nodes should enable the experimental study of multi-layer networks Boccaletti et al. (2014) and Asymmetry-Induced Synchronization Nishikawa and Motter (2016). One can introduce targeted perturbations to individual nodes, which may be useful for testing basin stability Menck et al. (2013, 2014) and controllability Liu, Slotine, and Barabási (2011).

The paper is organized as follows. In Section II, we describe the experiment for a single nonlinear map and investigate its dynamics. In Section III, we discuss how we create two coupled nonlinear maps using the space-time interpretation of time-delayed systems. In Section IV, we show how this same system can be used to implement larger networks of coupled maps, and we use it to study cluster and chimera states in both small and large networks of different topologies.

## Ii A single nonlinear map

We first discuss the simplest version of the experiment: operation as a single nonlinear map. Figure 1 provides a schematic of the experiment. Light from an 850 nm continuous-wave, fiber-coupled distributed feedback laser passes through an integrated LiNO_{3} electro-optic intensity modulator ( = 1.85 V) and is converted into an electrical signal by a photoreceiver. A flipflop performs a sample-and-hold operation on the electrical signal at a rate of = 10 kHz. This is implemented on an FPGA (Altera Cyclone V GT), and the delay caused by the clocked sample-and-hold is much longer than the optical delay. This makes time discrete and decouples consecutive time steps. The electronic output of the flipflop is amplified then applied to the RF port of the electro-optic modulator. The modulator provides the nonlinearity. An independent power supply controls the DC bias of the electro-optic modulator.

This experiment can be described by the following nonlinear map:

(1) |

where is the normalized voltage applied to the modulator at discrete time , and is the normalized round-trip gain. The normalized intensity of the light passing through the electro-optic modulator can be modeled as , where is the DC bias point of the electro-optic modulator. The sine-squared nonlinearity is intrinsic to all wave-interference devices, including our intensity modulator. With this experimental system, one could obtain any desired nonlinearity by performing a nonlinear operation on in the FPGA. Alternatively, one could introduce a different nonlinearity by using different, sufficiently fast nonlinear optical or electronic components.

Figure 2a shows an experimentally measured bifurcation diagram of the single nonlinear map with . Our experimental map can exhibit fixed point, periodic, and chaotic behaviors. A bifurcation diagram numerically simulated from Eq. (1) and shown in Fig. 2b, agrees well with the experiment and suggests that Eq. (1) is an accurate model.

If is too large, the voltage applied to the modulator may exceed the allowed limit for the device (for our modulator, the maximum voltage range is about 8V). This would limit the permissible values of to about 2.5. However, since is periodic, we resolve this problem by calculating modulo in the FPGA, allowing our experiment to operate at very high values of .

## Iii Two coupled nonlinear maps

In order to explore the dynamics of coupled nonlinear systems, one can envision coupling two nominally identical nonlinear maps. Such a system can be modeled as

(2) |

where is the coupling strength.

However, there are significant experimental difficulties in directly coupling two nonlinear maps. Two separate systems must be built and coupled; even then, the individual systems are only nominally identical.

One can, however, modify the single map described above to realize two coupled maps by adding time-delays as shown in Fig. 3. This is our approach. In this setup, the single oscillator with time-delays and a multiplexer (MUX) functions as two individual nonlinear maps coupled together. Every even time step of the full system is interpreted as an iteration of node 0, and every odd time step of the full system is interpreted as an iteration of node 1. This system can also be modeled by Eq. 2, but requires the assembly of only one experimental apparatus. Further, each of the two nodes are identical, since they utilize the same electronic and optical components. The individual maps now update at half the rate since the measurement of is performed serially.

In practice, the flip flop, time delays, time multiplexing, multiplication by and , and addition of the self-feedback and coupling components are all implemented on an FPGA. This allows for a compact and flexible implementation of the experiment: The experiment can be switched between the single nonlinear map described in section II and the coupled map system simply by reconfiguring the FPGA.

Figure 4 shows the measured and simulated root-mean-square (RMS) synchronization error of the coupled map system depicted in Fig. 3 and described by Eq. 2. We define the RMS synchronization error as

(3) |

where denotes an average over time. is zero in the case of a completely synchronized solution and approaches 1 in the limit that the and are uncorrelated.

One can determine the stability of the synchronized state by linearizing about the synchronized solution and calculating the Lyapunov exponent of the variational equation for the perturbations transverse to the synchronization manifold. The variational equation is

(4) |

where is a perturbation transverse to the synchronization manifold. We have calculated the Lyapunov exponent of Eq. 4 as a function of for fixed ; the results are shown in Fig. 4. The values of with negative Lyapunov exponent correspond exactly to the values where we observe synchronization in the simulations without noise. Further, all of the values where we observe synchronization in the experiment correspond to with negative Lyapunov exponent; however there are some narrow regions of with negative Lyapunov exponent where we do not observe synchronization in the experiment. This is due to noise in the experiment. Noise in the experiment comes from a variety of sources; some of these noise sources include discretization noise in the ADC and DAC, electronic noise in the DAC amplifier, and Johnson noise in the photoreceiver. Because the voltage applied to the modulator can have a DC component, significant electrical power can be dissipated in the modulator’s 50 input resistor. This can result in heating, and so there are also fluctuations in the intensity at the output of the modulator.

We model all of these sources of noise by applying additive Gaussian noise with standard deviation 0.02 to each normalized intensity at each time step:

(5) |

where are independent, identically distributed random variables taken from a Gaussian distribution with zero mean and unit variance. The synchronization error from the simulation with , shown in red in Fig. 4, shows good agreement with the experimentally measured result.

## Iv Complex networks of nonlinear maps

The experiment described in Sec. III can be adapted as shown in Fig. 5 to create an arbitrary network of coupled nonlinear maps. The only modification is in the FPGA configuration. The normalized voltage output of the FPGA is now a linear combination of the previous normalized intensity measurements. This can be interpreted as a network of coupled oscillators according to the space-time interpretation described in ref. Arecchi et al., 1992 and utilized in ref. Appeltant et al., 2011.

In the space-time interpretation, there are two “time steps”: the standard time and the network time . In particular, for a network of nodes, at standard time , the node is updated as follows: , where is the network time. is the normalized voltage applied to the modulator at standard time . The FPGA then measures and records , the normalized intensity that passes through the modulator. is stored in a shift register on the FPGA, and is used to update nodes in the future. The network updates at a rate of .

The system can be modeled in network time as

(6) | |||

(7) |

Eq. 6 clearly shows that this system can be understood to be a network of coupled maps.

We note that can be a Laplacian matrix, which indicates diffusive coupling and can lead to negative values of . This is not a problem because is determined from the intensities using signed arithmetic in the FPGA. We take advantage of the -periodicity of the modulator to allow us to output only positive voltages from our unipolar digital-to-analog converter.

Similar optoelectronic experiments have used the space-time interpretation to realize networks for reservoir computing Appeltant et al. (2011); Larger et al. (2012); Antonik, Haelterman, and Massar (2017) and to observe chimeras Larger, Penkovsky, and Maistrenko (2015). These systems use a time-invariant high-pass filter, which creates few nearest-neighbor coupling. In contrast, we are able to create arbitrary networks of coupled maps in a table-top experimental network by removing the high-pass filter and using a flipflop to sample-and-hold in the FPGA to obtain discrete time. Our nodes are maps; discrete time is a natural consequence of the space-time interpretation of systems with time delays Arecchi et al. (1992). Further, since is controlled by the FPGA in our experiment, we are able to explore networks with either identical (by keeping constant in time) or non-identical nodes (by varying periodically in time).

### iv.1 Chimeras in large ring network

In order to demonstrate that our experimental system can be used to study large networks, we now investigate chimera states in a ring network of 128 nodes. As shown in Fig. 6a, each node is coupled to the nodes on either side of it with equal strength ; this is the same coupling scheme used in ref. Hagerstrom et al., 2012 to observe chimeras. The nodes are coupled in a Laplacian fashion. The network can be modeled by Eq. 6 with the adjacency matrix given by

For this experiment, we set 2.67, 0, =128, and =52, which is a scaled version of the parameters used in ref. Hagerstrom et al., 2012 At these values of and , a single uncoupled map behaves chaotically. However, when coupled together with 0.0107, the chimera state displayed in Fig. 6b is observed. In this state, all nodes have period-4 dynamics, but there are two regions of spatial coherence separated by two regions of spatial incoherence. The space-time network interpretation is shown in Fig. 6c. The system starts with the oscillators uncoupled. They oscillate chaotically and independently, as shown. There is a long transient after the coupling is turned on (not shown). At the end of the transient, the system settles into the chimera state (shown starting at time ).

At a size of 128 nodes, the network updates at 10 kHz/128 78 Hz, which is an order of magnitude faster than previous coupled map experiments implemented with a SLMHagerstrom et al. (2012). Another advantage of our experiment, as noted earlier, is that all the nodes are identical. This removes any question of potential inhomogeneities among oscillators leading to the chimera states.

Since this network exhibits cyclic symmetry, it might be possible to implement in the previous versions of experiments using time delays to implement virtual nodes Appeltant et al. (2011); Larger et al. (2012); Larger, Penkovsky, and Maistrenko (2013, 2015); Antonik, Haelterman, and Massar (2017) as discussed above; however, this has never been done. The chimeras that have been observed in these systems are qualitatively different than those observed here.

### iv.2 Cluster synchronization in a 5 node network with symmetries

In the previous section, we demonstrated that our experimental system can be used to create realizations of networks with large numbers of nodes. In order to illustrate that our approach can be used to realize networks with arbitrary topology, we implement the network shown in Fig. 7a. Previous experiments using time delays and time-invariant filters to implement cyclically symmetric networks Appeltant et al. (2011); Larger et al. (2012); Larger, Penkovsky, and Maistrenko (2013, 2015); Antonik, Haelterman, and Massar (2017) cannot be used to study this network. We reproduce a result from ref. Sorrentino et al., 2016; this work explored the stability of cluster synchronization in a 5 node network with Laplacian coupling implemented on a SLM feedback system. We show that our time-delayed feedback network can also display cluster synchronization, but runs at rates two orders of magnitude faster than the SLM network. The network can be modeled by Eq. 6 with the adjacency matrix shown in Fig. 7a.

The cluster state we investigate is depicted in Fig. 7a: nodes 0 and 2 are synchronized, nodes 1 and 3 are synchronized, and node 4 is a singleton. The system parameters are , , and , in agreement with the parameters used to observe the same cluster state in the SLM experiment in ref. Sorrentino et al., 2016. An experimental time series of this cluster state is shown in Fig. 7b. The background color of each time step indicates the corresponding node (Fig. 7a) in the space-time network interpretation. The dotted black lines separate the network time steps . The dynamics of the system in the space-time network interpretation is shown in Fig. 7c. In Fig. 7b and 7c, it is clear that the desired cluster synchronization state is displayed by the experiment.

### iv.3 Chimeras in a 5 node, globally-coupled network

The network we discussed in the last section has previously been studied in great detail in ref. Sorrentino et al., 2016, and was used to demonstrate the flexibility and speed of our experimental system. We now investigate chimera states in a new network topology.

Experimental observations and stability analysis of chimera states in small networks of four globally-coupled continuous time oscillators have been studied previously Hart et al. (2016). It was predicted that similar chimera states should exist for larger globally coupled networks, and that it should be possible to determine their stability using the analysis techniques described in ref. Pecora et al., 2014; Hart et al., 2016; Sorrentino et al., 2016. However, due to experimental limitations on the number of nodes and links, chimeras in larger networks of globally coupled oscillators have not been experimentally observed in opto-electronic networks.

We now investigate chimera states in a network of five globally-coupled nodes in our experimental system. The chimera state we consider here is a dynamical state in which three of the nodes are synchronized, and the other two nodes are desynchronized both from the cluster of three and from each other. Group theoretical arguments suggest that such chimera states should exist in globally coupled networksHart et al. (2016). However, their stability determines whether they should be observable in experiments.

As in ref. Hart et al., 2016, we added a coupling time delay to observe chimera states in our globally-coupled network. This is straightforward to do with the FPGA hardware via shift registers. The network can be modeled as

(8) |

where as before and is the coupling delay. For this experiment, we take = 1 iteration.

Figure 8b shows a typical time series from a chimera in the experiment. The parameter values are , , and . As one can see from Fig. 2, an uncoupled oscillator is chaotic with these parameters. In this chimera, nodes 0, 2, and 3 are in the synchronous (coherent) region, and nodes 1 and 4 are in the desynchronized (incoherent) region. When we start from different initial conditions, different nodes end up in the coherent and incoherent regions. This is expected, since all of our nodes are identical.

In order to investigate the stability of these chimera states, we performed a stability analysis according to the methods described in ref. Pecora et al., 2014; Hart et al., 2016. We linearize about the chimera state to obtain the variational equations for the network, then use group theory-based techniques to pick out the perturbation directions transverse to the synchronization manifold. The stability of the chimera state is determined by the largest Lyapunov exponent (LLE) of these transverse variational equations: if the LLE is negative, the chimera state is stable. The result of this calculation for and is shown in Fig. 8c. We see that the chimera state is stable in the region from to . We were not able to calculate the Lyapunov exponent in the transverse direction for some values of . In regions where the chimera state is unstable and a more symmetric state (such as the globally synchronized state or the triplet-doublet state) is stable, one cannot calculate the stability of the chimera state. In this case the trajectory of the chimera state cannot be determined so one cannot linearize about the chimera state and calculate its stability in the usual way.

We also calculate the stability of the triplet-doublet cluster synchronous state, in order to show how the chimera state forms. The results of this calculation are also shown in Fig. 8c. For lower values of , the triplet-doublet cluster state is stable. As increases, the doublet cluster undergoes isolated desynchronization Pecora et al. (2014) and becomes unstable; however, the triplet cluster state remains stable. This results in a triplet-singlet-singlet state, which we call a chimera state.

## V Conclusions

We have developed a new experimental approach to study dynamics on large networks with arbitrary topology. This system provides a significant advance in network experiments in terms of speed, cost, size, and flexibility of coupling topologies. The experiment relies on the space-time interpretation of systems with time delayed feedback in order to create a network of coupled maps out of a single time delayed feedback loop. Because we are able to use a single loop, all of our nodes are truly identical We use optical and electronic components with no highpass filtering in order to remove the nearest-neighbor coupling found in previous systems that rely on the space-time interpretation to create networksAppeltant et al. (2011); Larger et al. (2012); Larger, Penkovsky, and Maistrenko (2013, 2015); Antonik, Haelterman, and Massar (2017). This, combined with an FPGA in the feedback loop to implement time delays, allows us to study networks with arbitrary topology, with adaptive capabilities, and at high speed. The FPGA also allows us to intentionally introduce heterogeneities in the nodes or links as well as additional sources of noise into our experiments, if desired.

When operated as a single nonlinear map, our experiment has displayed fixed point, periodic, and chaotic behavior for different parameters. We have also experimentally investigated two bidirectionally-coupled maps and networks of coupled maps, both large and small. We have observed cluster synchronization and two different types of chimera states in our system: chimeras in a small globally-coupled network and chimeras in a large ring network. We calculated the stability of the chimera states and cluster states in the 5 node, globally-coupled network, and we showed that the chimera state emerges when the doublet cluster of a triplet-doublet state undergoes isolated desynchronization.

The experiment currently updates at a rate of 10 kHz, which means that a network will update at a rate of 10 kHz/, where is the number of nodes in the network. Therefore the network updates slowly for extremely large networks (>1000 nodes). Some network experiments with fixed, non-reconfigurable topology operate at much faster time scales; for example, the globally-coupled mode-locked laser system in Ref. Viktorov et al., 2014 has dynamics in the range of tens of GHz. However, for all the networks that we consider here, our system provides a significant improvement in speed over previous experiments that allow easily reconfigurable, arbitrary network topology Hagerstrom et al. (2012); Pecora et al. (2014); Tinsley, Nkomo, and Showalter (2012). With suitable design it should be possible to increase the clock rate of our experiment many orders of magnitude, allowing experiments on arbitrary networks to occur at unprecedented speeds.

## Acknowledgements

The authors thank Francesco Sorrentino and Lou Pecora for helpful discussions. JDH and RR acknowledge support from ONR Grant No. N000141612481.

## References

- Huygens (1897) C. Huygens, Huygens Christiaan: oeuvres complètes, Vol. 1 (M. Nijhoff, 1897).
- Pikovsky, Rosenblum, and Kurths (2003) A. Pikovsky, M. Rosenblum, and J. Kurths, Synchronization: a universal concept in nonlinear sciences, Vol. 12 (Cambridge university press, 2003).
- Pecora and Carroll (2015) L. M. Pecora and T. L. Carroll, Chaos 25, 097611 (2015).
- Sorrentino et al. (2007) F. Sorrentino, M. di Bernardo, F. Garofalo, and G. Chen, Physical Review E 75, 046103 (2007).
- Abrams and Strogatz (2004) D. M. Abrams and S. H. Strogatz, Phys. Rev. Lett. 93, 174102 (2004).
- Panaggio and Abrams (2015) M. J. Panaggio and D. M. Abrams, Nonlinearity 28, R67 (2015).
- Zhou and Kurths (2006) C. Zhou and J. Kurths, Chaos 16, 015104 (2006).
- Ji et al. (2013) P. Ji, T. K. D. Peron, P. J. Menck, F. A. Rodrigues, and J. Kurths, Phys. Rev. Lett. 110, 218701 (2013).
- Abrams, Pecora, and Motter (2016) D. M. Abrams, L. M. Pecora, and A. E. Motter, Chaos 26, 094601 (2016).
- Kuramoto and Battogtokh (2002) Y. Kuramoto and D. Battogtokh, Nonlinear Phenom. Complex Syst. 5, 380 (2002).
- Hagerstrom et al. (2012) A. M. Hagerstrom, T. E. Murphy, R. Roy, P. Hövel, I. Omelchenko, and E. Schöll, Nat.Phys. 8, 658 (2012).
- Tinsley, Nkomo, and Showalter (2012) M. R. Tinsley, S. Nkomo, and K. Showalter, Nat. Phys. 8, 662 (2012).
- Schmidt et al. (2014) L. Schmidt, K. Schönleber, K. Krischer, and V. García-Morales, Chaos 24, 013102 (2014).
- Viktorov et al. (2014) E. A. Viktorov, T. Habruseva, S. P. Hegarty, G. Huyet, and B. Kelleher, Phys. Rev. Lett. 112, 224101 (2014).
- Arecchi et al. (1992) F. Arecchi, G. Giacomelli, A. Lapucci, and R. Meucci, Phys. Rev. A 45, R4225 (1992).
- Larger, Penkovsky, and Maistrenko (2013) L. Larger, B. Penkovsky, and Y. Maistrenko, Phys. Rev. Lett. 111, 054103 (2013).
- Larger, Penkovsky, and Maistrenko (2015) L. Larger, B. Penkovsky, and Y. Maistrenko, Nat. Commun. 6 (2015).
- Appeltant et al. (2011) L. Appeltant, M. C. Soriano, G. Van der Sande, J. Danckaert, S. Massar, J. Dambre, B. Schrauwen, C. R. Mirasso, and I. Fischer, Nat. Commun. 2, 468 (2011).
- Paquot et al. (2012) Y. Paquot, F. Duport, A. Smerieri, J. Dambre, B. Schrauwen, M. Haelterman, and S. Massar, Sci. Rep. 2 (2012).
- Larger et al. (2012) L. Larger, M. C. Soriano, D. Brunner, L. Appeltant, J. M. Gutiérrez, L. Pesquera, C. R. Mirasso, and I. Fischer, Opt. Express 20, 3241 (2012).
- Antonik, Haelterman, and Massar (2017) P. Antonik, M. Haelterman, and S. Massar, Phys. Rev. Applied 7, 054014 (2017).
- Martens et al. (2013) E. A. Martens, S. Thutupalli, A. Fourrière, and O. Hallatschek, Proceedings of the National Academy of Sciences 110, 10563 (2013).
- Blaha et al. (2016) K. Blaha, R. J. Burrus, J. L. Orozco-Mora, E. Ruiz-Beltrán, A. B. Siddique, V. Hatamipour, and F. Sorrentino, Chaos 26, 116307 (2016).
- Totz et al. (2015) J. F. Totz, R. Snari, D. Yengi, M. R. Tinsley, H. Engel, and K. Showalter, Phys. Rev. E 92, 022819 (2015).
- Wickramasinghe and Kiss (2013) M. Wickramasinghe and I. Z. Kiss, PloS one 8, e80586 (2013).
- Wickramasinghe and Kiss (2014) M. Wickramasinghe and I. Z. Kiss, Phys. Chem. Chem. Phys. 16, 18360 (2014).
- Pecora et al. (2014) L. M. Pecora, F. Sorrentino, A. M. Hagerstrom, T. E. Murphy, and R. Roy, Nat. Commun. 5, 4079 (2014).
- Kouomou et al. (2005) Y. C. Kouomou, P. Colet, L. Larger, and N. Gastaud, Physical review letters 95, 203903 (2005).
- Boccaletti et al. (2014) S. Boccaletti, G. Bianconi, R. Criado, C. I. Del Genio, J. Gómez-Gardenes, M. Romance, I. Sendina-Nadal, Z. Wang, and M. Zanin, Phys. Rep. 544, 1 (2014).
- Nishikawa and Motter (2016) T. Nishikawa and A. E. Motter, Phys. Rev. Lett. 117, 114101 (2016).
- Menck et al. (2013) P. J. Menck, J. Heitzig, N. Marwan, and J. Kurths, Nat. Phys. 9, 89 (2013).
- Menck et al. (2014) P. J. Menck, J. Heitzig, J. Kurths, and H. J. Schellnhuber, Nat. Commun. 5 (2014).
- Liu, Slotine, and Barabási (2011) Y.-Y. Liu, J.-J. Slotine, and A.-L. Barabási, Nature 473, 167 (2011).
- Sorrentino et al. (2016) F. Sorrentino, L. M. Pecora, A. M. Hagerstrom, T. E. Murphy, and R. Roy, Sci. Adv. 2, e1501737 (2016).
- Hart et al. (2016) J. D. Hart, K. Bansal, T. E. Murphy, and R. Roy, Chaos 26, 094801 (2016).