The dynamical systems found in nature are rarely isolated. Instead they interact and influence each other. The coupling functions that connect them contain detailed information about the functional mechanisms underlying the interactions and prescribe the physical rule specifying how an interaction occurs. A coherent and comprehensive review is presented encompassing the rapid progress made recently in the analysis, understanding, and applications of coupling functions. The basic concepts and characteristics of coupling functions are presented through demonstrative examples of different domains, revealing the mechanisms and emphasizing their multivariate nature. The theory of coupling functions is discussed through gradually increasing complexity from strong and weak interactions to globally coupled systems and networks. A variety of methods that have been developed for the detection and reconstruction of coupling functions from measured data is described. These methods are based on different statistical techniques for dynamical inference. Stemming from physics, such methods are being applied in diverse areas of science and technology, including chemistry, biology, physiology, neuroscience, social sciences, mechanics, and secure communications. This breadth of application illustrates the universality of coupling functions for studying the interaction mechanisms of coupled dynamical systems.

Spontaneous explosive emergent behavior takes place in heterogeneous networks when the frequencies of the nodes are positively correlated to the node degree. A central feature of such explosive transitions is a hysteretic behavior at the transition to synchronization. We unravel the underlying mechanisms and show that the dynamical origin of the hysteresis is a change of basin of attraction of the synchronization state. Our findings hold for heterogeneous networks with star graph motifs such as scale-free networks, and hence, reveal how microscopic network parameters such as node degree and frequency affect the global network properties and can be used for network design and control.

Heterogeneity in the degree distribution is known to suppress global synchronization in complex networks of symmetrically coupled oscillators. Scale-free networks display a great deal of heterogeneity, containing a few nodes, termed hubs, that are highly connected, while most nodes receive only a few connections. Here, we show that a group of synchronized nodes may appear in scale-free networks: hubs undergo a transition to synchronization while the other nodes remain unsynchronized. This general phenomenon can occur even in the absence of global synchronization. Our results suggest that scale-free networks may have evolved to complement various levels of synchronization.

For centuries isolation has been the main control strategy of unforeseen epidemic outbreaks. When implemented in full and without delay, isolation is very effective. However, flawless implementation is seldom feasible in practice. We present an epidemic model called SIQ with an isolation protocol, focusing on the consequences of delays and incomplete identification of infected hosts. The continuum limit of this model is a system of Delay Differential Equations, the analysis of which reveals clearly the dependence of epidemic evolution on model parameters including disease reproductive number, isolation probability, speed of identification of infected hosts and recovery rates. Our model offers estimates on minimum response capabilities needed to curb outbreaks, and predictions of endemic states when containment fails. Critical response capability is expressed explicitly in terms of parameters that are easy to obtain, to assist in the evaluation of funding priorities involving preparedness and epidemics management.