# Synchronization-based computation through networks of coupled oscillators

^{1}Departament de Física i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Terrassa, Spain^{2}Neuroheuristic Research Group, HEC Lausanne, University of Lausanne, Lausanne, Switzerland^{3}Laboratory of Biological Networks, Center for Biomedical Technology, Universidad Politécnica de Madrid, Madrid, Spain^{4}Complex Systems Group and GISC, Universidad Rey Juan Carlos, Madrid, Spain^{5}Department of Experimental and Health Sciences, Universitat Pompeu Fabra, Barcelona Biomedical Research Park, Barcelona, Spain

The mesoscopic activity of the brain is strongly dynamical, while at the same time exhibits remarkable computational capabilities. In order to examine how these two features coexist, here we show that the patterns of synchronized oscillations displayed by networks of neural mass models, representing cortical columns, can be used as substrates for Boolean-like computations. Our results reveal that the same neural mass network may process different combinations of dynamical inputs as different logical operations or combinations of them. This dynamical feature of the network allows it to process complex inputs in a very sophisticated manner. The results are reproduced experimentally with electronic circuits of coupled Chua oscillators, showing the robustness of this kind of computation to the intrinsic noise and parameter mismatch of the coupled oscillators. We also show that the information-processing capabilities of coupled oscillations go beyond the simple juxtaposition of logic gates.

## 1. Introduction

It has been established that the healthy brain operates in a highly coordinated way that involves different neural oscillations spanning through multiple spatiotemporal scales (Freeman, 1975; Singer, 1993, 1999; Başar et al., 2000; Varela et al., 2001; Ward, 2003; Buzsáki and Draguhn, 2004). Even though these oscillatory rhythms may have different synchronization properties (Kopell et al., 2000), they have been explained as a result of the balance between excitatory and inhibitory neurons in a network (Börgers and Kopell, 2003; Börgers et al., 2005). In turn, this synchronous firing may subserve complex coordinated patterns of spiking activity which may be transmitted in large neural networks with high temporal accuracy over long distances (Abeles, 1991; Rodriguez et al., 1999; Abeles et al., 2004; Asai et al., 2008; Asai and Villa, 2012; Barardi et al., 2014b). It has also been accepted that the oscillatory activity exhibited by brain signals such as local field potentials (LFP), electroencephalograms (EEG) and magnetoencephalograms (MEG), arises from the synchronized activity of large neuronal assemblies. Such collective dynamics throughout the different scales in the brain is likely to determine the functional role of normal and aberrant synchronization mechanisms during adaptive and cognitive processes as well as brain diseases (Del Prete et al., 2004; Iglesias and Villa, 2010; Pons et al., 2010; Villa and Tetko, 2010). Synchronization's role in coordinating and processing information at different spatiotemporal scales has been also stressed (Lachaux et al., 1999; Stam and de Bruin, 2004; Busáki, 2006; Malagarriga et al., 2015). For instance, synchronization-based selectivity of visual response has been studied in the context of monkeys and cats (Castelo-Branco et al., 2000; Womelsdorf et al., 2006) or even in humans (Rodriguez et al., 1999). Besides, synchronization participates in the odor perception (Stopfer et al., 1997; Laurent et al., 2001; Blumhagen et al., 2011) and coherence of stimuli also affects the selective capability of oscillatory networks (Börgers and Kopell, 2008; Börgers et al., 2008). The processing and computation mechanisms based in all this diversity of synchronized elements has also been studied in detail (Engel et al., 2001; Fries, 2009; Maris et al., 2013; Nikolić et al., 2013; Womelsdorf et al., 2014). Thus, the interaction of different synchronized ensembles of neurons (Womelsdorf et al., 2007; Wulff et al., 2009) plays a role in tasks like learning item-context associations (Tort et al., 2009), selective attention (Fries et al., 2001, 2008; Womelsdorf and Fries, 2007; Bosman et al., 2012), or even conscious perception (Melloni et al., 2007; Levy et al., 2013). At the larger spatial scale in the brain, synchronization participates in the control of task-switching (Phillips et al., 2013) and is studied routinely in normal and abnormal EEG and MEG recordings (Stam, 2005).

Even though much progress in the understanding of these synchronization mechanisms has been gained during many years, it is not fully understood yet how these synchronization relations are established with the participation of different scales simultaneously, or how they operate at the same time without interfering with each other (Barardi et al., 2014a). So, for instance, the information processing capacity of the brain operating under multiple scales has been described very often in terms of logic calculus. At the most microscopic level, the idea of logic calculus based on neuronal activities was embedded in the seminal work of McCulloch and Pitts (1990). Neuronal circuitry performing logic operations was physically implemented in cell cultures of *in vitro* models of selected brain areas (Feinerman et al., 2008; Wolf and Geisel, 2008). This approach is mainly based on action potentials and on the connectivity within the network, rather than on a dynamical analysis of the ongoing activity. At the cellular level neurons have revealed that, in addition to behaving as a bistable system, they can be driven into a continuous oscillation by means of selected voltage-dependent inward currents controlled by intracellular calcium concentrations (Contreras and Steriade, 1995; Hughes et al., 2002; Crunelli et al., 2005). Besides, from the microscopic point of view, neurons may coordinate their firing in response to incoming stimuli, opening the way to a neurocomputing paradigm characterized by different synchronized states where the neurons oscillate with equal frequencies and specific phase relationships (Hoppensteadt and Izhikevich, 2000; Zanin et al., 2011). By associating logical states to the dynamics of coupled oscillators, all usual Boolean operations can be implemented and a full computational model can be obtained (Xu et al., 2004). Beyond the cellular level of neuronal oscillators, it was recently demonstrated that circuits of neurons embedded within a large-scale network of cortical cells were able to express logic functions that are dependent on complex spatiotemporal patterns (Vardi et al., 2013; Menon and Sinha, 2014). This type of analysis can also be made at the mesoscale. Large brain circuits are frequently described as networks of nodes associated with neuronal assemblies, evolving at the mesoscopic scale, in such a way that their dynamics can be considered as that of limit-cycle oscillators subjected to weak forcing and coupling. Phase-reduction theory has revealed synchronization to be among the most relevant features that determine the dynamical states of these systems (Pikovsky et al., 2003; Brown et al., 2004). Furthermore, coupled oscillator theory has established the conditions that allow all the nodes, or a subgroup of them, to operate in one of several synchronization regimes, including complete, lag, generalized, and phase synchronization (Boccaletti et al., 2002; Li and Chen, 2004). In fact, recent work shows that in networks of mesoscopic brain oscillators different forms of synchronization might coexist (Malagarriga et al., 2015). This phenomenon enlarges the processing capacity of neural oscillators, and endows the corresponding networks with stability, flexibility and robustness against perturbations (Zanette, 2004).

In this paper we present a combination of a theoretical approach and its experimental implementation that may shed light on possible mechanisms of brain computation based on synchronized oscillations. Specifically, we show that networks of neural mass oscillators may process inputs in a complex Boolean-like manner. So, the combination of all the fluctuating inputs received by each oscillator in the network determines the global network dynamical response. Every combination of inputs received by the nodes produces a synchronization pattern that relates the dynamics of all the nodes of the network. We interpret the different synchronization states (i.e., phase, generalized, lag, or complete synchronization) of two oscillators in the network as Boolean variables that allow to classify the response of inputs onto pairwise logic gates of different nature. By tuning the characteristics of the oscillatory input acting upon the neural masses, the pairwise coordinated activity changes accordingly. Each form of synchronization brings different information in terms of the synchronized motion in phase space, thus providing additional characterization of such incoming stimuli in the form of coordination evolution. We show that several distinct logical operations can be implemented, in this way, by the same type of neural mass network. However, far from suggesting that the processing capacity of the brain is the result of a more or less complex boolean circuitry, we postulate that it results from its very complex collective dynamics in which synchronization may play a very relevant role (Buzsáki and Draguhn, 2004). The capacity of synchronized oscillations to perform Boolean-like operations is demonstrated by a physical implementation in the form of an electronic system consisting of coupled Chua oscillators driven by oscillatory input signals. Finally, we show that the processing possibilities of larger networks composed by this type of systems go beyond the simple juxtaposition of logic gates.

## 2. Materials and Methods

In this Section we describe how the networks of neural mass models and Chua oscillators are implemented. We also describe the numerical methods used and how the synchronization relations between the different nodes of the network are analyzed.

### 2.1. Networks of Neural Mass Oscillators

In order to describe theoretically the computation capabilities of a network of cortical columns acting as a cluster of interacting logic gates, we consider a network of neural mass models as described by Jansen et al. (Jansen et al., 1993; Jansen and Rit, 1995). In the Jansen model, all the neurons in each cortical column are classified into three different groups: pyramidal neurons (responsible of the signal measured in the scalp), excitatory interneurons and inhibitory interneurons. The three populations interact with each other through excitatory and inhibitory connections. Specifically, the population of pyramidal neurons receives both excitatory and inhibitory inputs, in the form of a feedback loop, from the interneuron populations. At the same time, the pyramidal population sends an excitatory input to the interneuron populations inside its cortical column in a feed-forward manner. The pyramidal population also receives an excitatory (inhibitory) input from the pyramidal (inhibitory interneurons) neighboring columns that form the network of neural mass oscillators.

The dynamics of each column in the network is described as follows. Each population converts the average pulse density of action potentials arriving to the population from different origins, ∑*p*_{i}(*t*), into an average excitatory postsynaptic membrane potential, *y*_{e}(*t*), through the following expression:

A similar expression with different constant terms, *B* and *b*, substituting *A* and *a*, respectively, works for average inhibitory postsynaptic membrane potentials, *y*_{i}. *A* (*B*) is related with the maximum height of the excitatory (inhibitory) postsynaptic potential, whereas *a* (*b*) represents the inverse of the membrane time constants and the dendritic delays. At the somas of the neurons forming the populations, the net average postsynaptic potential (PSP) input for the considered population, *m*(*t*) = *y*_{e}(*t*)−*y*_{i}(*t*), is transformed into an average density of spikes, *p*_{i}(*t*) = *C*_{i}*S*[*m*(*t*)], which will become part of the input of other populations. The transfer function *S*[*m*(*t*)] is given by

where *e*_{0} determines the maximum firing rate of the neural population, ν_{0} sets the net PSP for which a 50% firing rate is achieved and *r* is the steepness of this sigmoidal transformation.

The dynamical activity of the three populations in each cortical column *i* of the network follows from (Malagarriga et al., 2014, 2015):

where ${{y}}_{{0}}^{{i}}$ represents the PSP that feeds the interneurons populations, and ${{y}}_{{1}}^{{i}}$ (${{y}}_{{2}}^{{i}}$) represents the excitatory (inhibitory) PSP that feeds the pyramidal population. The intensity of the excitatory (inhibitory) coupling of columns with their neighboring columns is given by α_{ij} (β_{ij}). Moreover, each column may receive a time dependent input *I*^{i}(*t*) composed by a constant input, ${\overline{{p}}}^{{i}}$, and a periodic external stimulus coming from other brain structures or the sensory system. We represent this periodic input as a sinusoidal driving, i.e., δ^{i}(*t*)*sin*(2π*f*^{i}*t*+ϕ^{i}). Besides, each column may receive a random excitatory contribution onto the pyramidal cells, χ_{i}(*t*), which can be associated with a stochastic process occurring at a cellular level. The contribution of the column *i* to the EEG activity measured in the scalp is given by ${{y}}_{{1}}^{{i}}{-}{{y}}_{{2}}^{{i}}$. Thus, we will analyze the activity of each cortical column considering the evolution of ${{y}}^{{i}}{=}{{y}}_{{1}}^{{i}}{-}{{y}}_{{2}}^{{i}}$. The model presented here has an extensive repertoire of dynamical states, being able to produce periodic, quasi-periodic or chaotic (Skarda and Freeman, 1987; Spiegler et al., 2011; Malagarriga et al., 2015) behavior. It also exhibits excitatory/inhibitory segregation (Malagarriga et al., 2014, 2015) depending on the choice of α_{ij} and β_{ij} values in the network. These coupling parameters also fix a whole set of synchronized regimes that may coexist in the network (Malagarriga et al., 2015). By choosing carefully the input protocol in terms of ${\overline{{p}}}^{{i}}$, δ^{i}(*t*)*sin*(2π*f*^{i}*t*+ϕ^{i}), χ^{i}(*t*), α_{ij}, and β_{ij}, we will establish different synchronization-based logic gates.

### 2.2. Networks of Chua Oscillators

The feasibility and robustness of our theoretical results are proven experimentally by considering a network of electronic oscillators. In this case, the dynamics for every node *i* in the network is determined by the classical Chua circuit (Madan, 1993) (see Figure 1). This system is described by the following equations (Kennedy, 1992; Wagemakers et al., 2007):

where ${{{v}}_{{A}}}^{{i}}$ and ${{{v}}_{{B}}}^{{i}}$ are the voltages of the two capacitors (*C*_{1} and *C*_{2}) and ${{{i}}_{{L}}}^{{i}}$ is the intensity through the coil of the circuit (*L*_{1}). *s*(*t*) is an external oscillatory input signal that we can (“1”) or can not (“0”) activate, controlling the dynamical state at which the Chua operates. The strength of the input signal is proportional to the inverse of *R*_{in}. We will tune its value depending the type of logic gate we will be considering. ${f}({{{v}}_{{A}}}^{{i}})$ is a piece-wise (non-linear) function given by:

with *B*_{p} = 1.7 V being the breaking point of the piece-wise function. The values of *G*_{1} and *G*_{2} are obtained from different resistances of the electronic circuit:

and the rest of the parameters are *V*_{cc} = 15V, *V*_{EE} = −15V, *R*_{1} = 222 Ω, *R*_{2} = 22 kΩ, *R*_{3} = 2.2 kΩ, *R*_{4} = 3.3 kΩ, *C*_{1} = 10 nF, *C*_{2} = 100 nF, *L*_{1} = 20 mH, and *R*_{5} = 1.38 kΩ.

**Figure 1. Electronic implementation of the Chua circuit**. Two TL082 operational amplifiers are the core of the non-linear part of the circuit which follows the function given in Equation 10. The input signal (0/1) is introduced through the capacitor *C*_{1}, while the output of the circuit is the voltage (*v*_{B}) of both *C*_{2} and *L*_{1}.

The coupling circuit between the different Chua circuits (not shown in Figure 1), represented as the last term in Equation (chua_Secondeq), consists on a voltage follower placed at the output of ${{v}}_{{B}}^{{j}}$ combined with a coupling resistance *R*_{c}, whose value controls the amount of unidirectional coupling from ${{v}}_{{B}}^{{j}}$ to ${{v}}_{{B}}^{{i}}$. The coupling from unit *i* to unit *j* is introduced in a similar way, thus leading to an effective bidirectional coupling. Values of the coupling and input resistances of the different gates implemented in the experiments are summarized in Table 1.

The input signal *s*(*t*) used as the “1” input for the logic gates is introduced through the variable *v*_{A} by means of a voltage follower and a coupling resistance *R*_{in}. The signal *s*(*t*) was recorded previously from a Chua circuit with the same parameters as the ones used in the experiments. Both the input and the output signals are characterized by a chaotic behavior (see Madan, 1993 for a detailed analysis of the bifurcation diagram of the Chua circuit, and Wagemakers et al., 2007 for an implementation of a similar circuit in the chaotic regime). Finally, outputs of both gates are sent to an *integrator* system (*C*), by means of a bidirectional coupling, adjusted by an integrating resistance *R*_{int} (see Figure 4 for details).

### 2.3. Numerical Methods

The equations of the Jansen model were numerically solved by means of Heun's method (García-Ojalvo and Sancho, 1999; Toral and Colet, 2014), as performed in Malagarriga et al. (2015). We generated random numbers using standard GSL routines to set different initial conditions and to introduce noise in the dynamics. We implement numerically the noise term χ_{i}(*t*) using:

where ξ_{i} is the noise amplitude and Δ*t* is the integration time step, whereas η(*t*) is a number resulting from a white noise Gaussian distribution with zero mean and variance equal to 1 (García-Ojalvo and Sancho, 1999). Each simulation of the Jansen model had a time step of 1 ms and ran over a total time of *t* = 25 s in Figures 2, 3, and *t* = 50 s in Figure 5. An initial interval of 10 s was omitted to avoid transients in Figures 2, 3, whereas no transient dynamics were found for the simulations presented in Figure 5.

**Figure 2. Implementation of binary logic gates formed by pairwise connected cortical columns**. (**A**) Cartoon depicting two connected cortical columns. Both columns are bidirectionally coupled through both excitatory–excitatory (black solid arrows) and inhibitory–excitatory (blue solid lines) contacts. The intensity of such connections is weighted by α_{ij}=α_{ji} and β_{ij}=β_{ij} coupling strengths, respectively. The pyramidal population within each cortical column additionally receives a constant background pulse density $\overline{{{p}}^{{i}}}$, an oscillatory input defined by an amplitude δ^{i} and a frequency *f*^{i} and a stochastic input contribution (see Equation 5). The characteristics of these inputs define states which feed this binary logic gate. The out-coming signals, ${{y}}^{{i}}{=}{{y}}_{{1}}^{{i}}{-}{{y}}_{{2}}^{{i}}$, of the two columns may (*C*(*y*^{1}, *y*^{2}) ≈ 1) or may not (*C*(*y*^{1}, *y*^{2}) ≈ 0) be correlated. The correlation value of the timetraces of the two nodes defines the output state of the gate for the given input stimulation (see Table 3). (**B**) Several binary logic gates can be obtained from the system shown in panel **(A)** if the appropriate combination of parameters and input protocols are selected (see Table 2). Correlated signals, in the form of PS, LS, or CS, result in a “1” state of synchronization, however, not correlated signals are interpreted as a “0” state of synchronization. In the panel **(B)**, both, the time traces of the oscillators and the state resulting from their synchronization are shown for all the binary combinations of input protocols (00, 01, 10, and 11). Several truth tables defined in this way for the logic gates AND, NAND, NOR, OR, XNOR, and XOR are shown.

**Figure 3. Implementation of two logic gates embedded in a network of coupled cortical units**. (**A**) Cartoon depicting the network of coupled units. Nodes *A*_{1} and *A*_{2} implement a NOR gate (in a CS regime) while nodes *B*_{1} and *B*_{2} implement an AND gate (in a PS regime). (**B**) The implementation of such logic gates is not altered by the state of synchronization of neighboring pairs, i.e., nodes *B*_{1} and *B*_{2} are capable of implementing a AND gate for any synchronized state of nodes *A*_{1} and *A*_{2} (shown only for a single input state combination for nodes *A*_{1} and *A*_{2}). However, the dynamics of *B*_{1} and *B*_{2} are indeed altered, as time traces are no longer the same for each of the four realizations. In fact, such changes in the dynamics may produce long transients before falling into the corresponding synchronized regime (see third pair of time traces between 30 and 35 s). (**C**) In some cases, when tuning the input, the dynamics of the oscillators no longer allows to implement the *desired* logic gates. For an *A*_{1} and *A*_{2} (constant) input configuration of 00, the implementation of an AND gate in nodes *B*_{1} and *B*_{2} fails (see first time traces, where no synchronized output configuration is obtained for the 11 input for *B*_{1} and *B*_{2} as required by the AND truth table). Another type of failure, observed in the second example of timetraces in the panel **(C)**, is the impossibility to return to the initial state of desynchronization. This results in a history-dependent logic gate.

### 2.4. Analysis

In order to determine the state of our logic gates, we need to quantify the synchronization between pairs of nodes. Each synchronized state is defined by a characteristic functional relationship between the dynamics of the interacting elements (Boccaletti et al., 2002). Phase synchronization (PS) entails a constant phase difference in time between the coupled oscillators, whereas amplitudes remain uncorrelated (Rosenblum et al., 1996). Generalized synchronization (GS) is characterized by a complex functional relationship between the dynamics of the oscillators that can only be unveiled by the auxiliary system approach (Abarbanel et al., 1996) or the nearest neighbor method (Moskalenko et al., 2012). In turn, lag synchronization (LS) imply a constant time shift between the signals of the two oscillators, with amplitudes being completely correlated, whereas in complete synchronization (CS) no time shift is present. Although CS is not likely to be achieved in a neural context, here we want to emphasize the possibility of the system to reach several synchronized states and, thus, we include here this type of coordinated dynamical evolution as a possible coding state. Accordingly, we computed in the neural mass model the cross-correlation and a phase synchronization index between the output signals of the cortical columns, *y*^{m}(*t*) (Lachaux et al., 1999). Cross-correlations allow us to distinguish between zero-lag (complete) synchronization (CS) and lag synchronization (LS), whereas the phase synchronization index provides evidence of phase synchronization (PS). Generalized synchronization entails the implementation of a three input logic gate, which has been excluded for the sake of clarity. The phase ϕ^{m}(*t*) of the output signal is obtained from the Hilbert transform of *y*^{m}(*t*) (Rosenblum et al., 1996; Mormann et al., 2000). From the phase, the phase synchronization index γ of two oscillators, 1 and 2, is calculated from Δϕ^{12}(*t*)=ϕ^{1}(*t*)−ϕ^{2}(*t*) as (Mormann et al., 2000; Quian Quiroga et al., 2002):

For the Chua circuits, we evaluate the synchronization error as the average of the difference between the outputs of two systems (e.g., units *A*_{1} and *A*_{2} in Figure 4A):

with *T* being the total number of time steps. We use a similar expression in terms of *y*^{1} and *y*^{2} for the neural mass oscillators. We consider two oscillators to have complete synchronization (or lag synchronization if there is a time shift of the signals) when the values of the synchronization error are lower than a certain threshold ϵ_{th} (ϵ_{th} = 0.10 V in the case of Chua oscillators, ϵ_{th} = 0.01 mV for Jansen oscillators). On the other hand, phase synchronization arises for high values of the phase synchronization index (with a threshold of γ_{th} = 0.85 for both cases) and, at the same time, high values of the synchronization error (ϵ > ϵ_{th} = 0.10 V, ϵ_{th} = 0.01 mV, for the Chua and Jansen oscillators, respectively).

**Figure 4. Experimental implementation of integrated logic gates**. **(A)** Qualitative description of the experimental setup: Nodes *A*_{1} and *A*_{2} form a XNOR gate (in a CS regime) while nodes *B*_{1} and *B*_{2} implement an AND gate (in a PS regime). The output of the two sync-based gates is integrated through node C. **(B)** Time series of the XNOR gate. Functioning of the gate relies on the complete synchronization of units *A*_{1} and *A*_{2}. The upper time trace, obtained for low values of the coupling with node C (*R*_{int} = 100 kΩ), shows the different outputs of the truth table of the XNOR gate (see Table 5). In the bottom signal, coupling with node C is increased (*R*_{int} = 25 kΩ), and the gate begins to fail. **(C)** Time series of the AND gate. In this case, the functioning of the gate relies on the phase synchronization of units *B*_{1} and *B*_{2}. The upper time trace, obtained for low values of the coupling (*R*_{int} = 100 kΩ) with node C, corresponds to the different outputs of the truth table of the AND gate (see Table 5). In the bottom time trace, coupling with node C is increased (*R*_{int} = 25 kΩ), and the AND gate also begins to fail.

## 3. Results

### 3.1. Small Networks of Neural Mass Oscillators

The input received by two Jansen oscillators within a network may determine their state of synchronization (Malagarriga et al., 2015). Thus, labeling the input signals as “0” and “1” (arbitrarily) and the synchronization of the two oscillators as “0” (“1”) when the oscillators are not (are) synchronized, we can interpret the dynamical response of these two nodes to their inputs in terms of binary logic gates. By changing the characteristics of the inputs and the excitatory/inhibitory coupling strengths between the two oscillators, α_{ij}∕β_{ij}, several types of binary logic gates may be created.

In order to understand how these clusters of gates operate in networks of oscillators we first study a simpler situation formed only by two bidirectionally coupled Jansen oscillators (See Figure 2, Table 2). In this configuration, each oscillator (*A*_{1} and *A*_{2} from now on) receives an external oscillatory and stochastic input, i.e., Input *A*_{1} (δ_{A1}, *f*_{A1}, ξ_{A1}) and Input *A*_{2} (δ_{A2}, *f*_{A2}, ξ_{A2}), that adds to the other pulse density contributions. The output signal of each oscillator is the difference between the Excitatory Postsynaptic Potential (EPSP) and the Inhibitory Postsynaptic Potential acting upon the pyramidal population (${{y}}^{{i}}{=}{{y}}_{{1}}^{{i}}{-}{{y}}_{{2}}^{{i}}$) (see triangles in Figure 2 and the Section 2). The evaluation of the synchronization of the output of the two elements that form the gate, *y*^{i} and *y*^{j}, determines the dynamical response of the system. Other characteristics of the response may also inform about the nature of the inputs received by the oscillators. For instance, it has been shown (Huang et al., 2011; Spiegler et al., 2011) that the driving of Jansen oscillators by a periodic input may result in chaotic, quasi-periodic or periodic dynamical evolutions. So, taken as a whole, both the synchronization state of two oscillators and their dynamical state (e.g., chaotic or oscillatory), inform about the characteristics of the input received by the gate. As mentioned above, we have assigned arbitrarily the “0” and “1” values to the absence of synchronization and its presence, respectively. However, our approach can be generalized in several ways. We could, as well, have labeled each synchronization type (PS, LS, or CS) with a different label increasing the possible number of output states. We could also consider the correlation value to define the output state allowing for fuzzy logic computation. We will not study such complex scenarios, in order to keep our approach as simple as possible.

**Table 2. Values of the coupling parameters and input stimulation protocols for inputs labeled as “1” in the implementation of the binary logic gates shown in Figure 2**.

**Table 3. Synchronization errors and phase synchronization indices for each logic gate implemented in Figure 2**.

With all these ingredients we have constructed truth tables based on synchronized states as shown in Figure 2B. The time traces show the online implementation of AND, NAND, NOR, OR, XNOR, and XOR gates that operate in different synchronization regimes. To obtain them, different input protocols and coupling strength relations between the two nodes are needed (see Table 2). However, a single pair of connected oscillators with fixed coupling strengths α_{ij}∕β_{ij} may organize its response in different ways when the input characteristics (frequency, amplitude, etc.) change (e.g., see AND, OR, XOR gates in Table 2). So the same brain circuit represented by a single pair of oscillators may behave as different logic gates depending of the inputs they receive. This ability to classify the response of the system in different ways depending of the inputs received results in a higher flexibility for the information processing capacity of the network. The network is not a static circuit which computes passively the response to the inputs but a whole collection of circuits (based in the logic gates described here) which reconfigures itself depending of the input received. This type of complex networks of logic gates which reconfigure themselves dynamically to adapt its response to specific inputs has been found in other natural systems such as the signaling networks in eukaryotic cells (Domedel-Puig et al., 2011; Rué et al., 2012). As said above, we do not claim, however, that the brain is a circuit of boolean gates but a highly complex dynamical system able to process information in a very sophisticated manner at the mesoscale.

Another trait of these logic gate implementations is related with the complex dynamical evolutions that the neural mass oscillators may show. For instance, for chaotic states, time evolutions depend strongly on the initial conditions. This feature introduces a very strong link between the dynamical evolution of one oscillator and its input. Resulting from this, we argue that our implementation may possess *reversible logic* characteristics (Bennett, 1973), as it allows to recognize from the out coming signals which of the two elements in the logic gate is receiving an input. The details of such paradigm are however out of the scope of this work.

The extension of the previous results to a network of interacting cortical columns leads to spatially distributed computation. In this case, each node in the network receives inputs that determine their synchronization state with other network nodes, following one of the implemented logic gates shown in Figure 2B. In that way, different regions of the network may act as gates that process inputs from the same or different origins, enriching the processing capabilities of the network even further. In this paper, we present a first step toward analyzing these capabilities by considering only relatively simple network motifs. Even if these motifs behave as well-established logic gates in isolation, their behavior when operating within larger networks might be more complex or even unexpected, as shown in Figure 3. In those cases, the whole network can process the whole set of its inputs following complex multidimensional logic rules.

Figure 3A shows a network of five interacting cortical columns that implement two logic gates: NOR (with CS synchronization) and AND (with PS synchronization). The peripheral pairs of columns are capable of working independently, even though their dynamics are influenced by the behavior of the whole network. Figure 3B shows the truth table and the corresponding time traces for the response of the AND gate (columns *B*_{1} and *B*_{2}), while columns *A*_{1} and *A*_{2}, which implement the NOR gate, receive different inputs. In this case it is worth mentioning that all four time traces for the AND implementation, which are symmetric realizations in terms of initial conditions, display distinctive dynamics depending on the *A*_{1} and *A*_{2} inputs. However, they do not lose the capacity of remaining synchronized/non-synchronized. Nevertheless, some realizations display states of long transient dynamics (see the “0” return state in all pairs of time traces) which are indeed not beneficial for fast brain computation. One possible application of transient dynamics in this system can be related with the ability to discern between the different input scenarios thanks to the different length of the transient dynamical evolutions. Such transient dynamics coding may be related with brain functions as proposed in Rabinovich and Varona (2011). Thus, overall, such dynamical behavior shows that the network implementation of logic gates is stable in terms of synchronization and arises from interdependent dynamics.

Despite being a robust dynamical feature, logic gate implementation strongly depends on the type of input and the coupling strength ranges between columns. In this regard, Figure 3C shows two situations in which there is a wrong output configuration of an AND gate depending on the input applied to nodes *A*_{1} and *A*_{2} (here δ_{A1} = δ_{A2} = 300 Hz, while all other input characteristics are the same as those in Table 4). *B*_{1} and *B*_{2} both receive the standard oscillatory input protocol. In the first pair of time traces shown, the output displayed by nodes *B*_{1} and *B*_{2} does not show correlated dynamics (when *B*_{1} and *B*_{2} receive “1” inputs) as expected for a correct AND implementation. In turn, the second pair of time traces shows how, after displaying a correct output for the initial four pairs of input states, the subsequent return to the first output configuration is no longer possible. Such history dependent behavior entails the impossibility, for these conditions, of a forward implementation of the previous logic gate (AND) but results in the implementation of a state- (or history-) dependent one. These situations result only for specific input patterns. However, the response to other inputs may produce well-behaved logic gates.

**Table 4. Values of the parameters used in the implementation of the NOR and AND logic gates appearing in Figure 3**.

### 3.2. Small Networks of Chua Oscillators

In the previous section we have shown theoretically that information processing in terms of synchronization-based logic gates for small networks of neural mass oscillators is possible. In this section, we show that these results can also be obtained experimentally by means of electronic circuits. Even though we will use another type of system, the Chua oscillator (Madan, 1993), that can show several dynamic regimes, we are going to consider only situations in which the circuits operate in the chaotic regime (see Section 2.2).

Such an experimental implementation will be a proof of the fact that systems as those considered in the previous section are robust against the parameter mismatch of electronic circuits and the intrinsic noise of real systems. Chaotic systems, such as Chua circuits, may display different kinds of synchronization that can be tuned through the coupling strength between the networked elements and the particular dynamics of an input signal (Arenas et al., 2008). Even though the nature of the oscillators is different from those considered theoretically, these systems have the ability to process information in the same way as those described in the previous section. We have designed and implemented experimentally several logic gates as those described in Figure 2B. Nevertheless, in what follows we are going to focus on the description and integration of a XNOR and AND gates in a network, the former based on complete synchronization and the latter on phase synchronization.

Figure 4A shows schematically the system integrating the outputs of a XNOR and AND gate. This circuit represents a network similar to that studied in Figure 3. It consists of an integrated circuit formed by five Chua oscillators, two of them implementing a XNOR gate (nodes *A*_{1} and *A*_{2}), other two forming an AND gate (nodes *B*_{1} and *B*_{2}) and the fifth Chua circuit (nodes *C*) integrating the output of both gates. The input of each dynamical unit can be either “1” (when a complex signal is injected into the node) and “0” (in the absence of an input signal). In order to assess the emergence of the synchronized time evolutions we have computed the synchronization errors of the two pairs of nodes that implement the logic gates (see Equation 15). Such errors allow us the determine a proper threshold to define each synchronized state, and thus each output state, in the truth tables of the logic gates. Time traces and the corresponding truth table for the XNOR gate are shown in Figure 4B, while synchronization errors determining the gates' outputs are summarized in Table 5.

**Table 5. Synchronization errors and truth tables of the XNOR gate (with complete synchronization) and the AND gate (with phase synchronization) embedded in a network of coupled Chua oscillators. The thresholds to establish synchronized states are shown in Materials and Methods section**.

**Table 6. Values of the parameters used in the implementation of the flip-flop operational gate shown in Figure 5**.

The AND gate is based on phase synchronization, and it is obtained by a fine tuning of the input and coupling resistances. In this case, we must reach a high phase synchronization but preventing complete synchronization, the latter leading to unavoidable (i.e., trivial) matching of the phases of the oscillators. To guarantee that we have phase, and not complete, synchronization we checked at the same time the synchronization error ϵ and the phase synchronization index γ (see Equation 14). A combination of a high ϵ and a high γ is the signature of phase synchronization between the two units forming the gate. Figure 4C shows the time series and the corresponding truth table for the AND gate. Table 5 summarizes the values of ϵ and γ that lead to a successful performance of both gates.

Finally, it is worth mentioning that a tuning of the integrator resistances (see Table 1) allows the correct functioning of the ensemble of gates. When the integrator resistances are decreased (i.e., the coupling with the central node *C* is increased) to values close to *R*_{int}≤25 kΩ the functioning of the gates begins to fail due to interferences between them (see bottom panels of Figures 4B,C). Such feature was also present in the case of the Jansen oscillators, in which increasing the coupling strength between the peripheral cortical columns and the central node led to an unstable implementation of logic gates.

### 3.3. Neural Mass Implementation of a Flip Flop Circuit

Other extended systems able to perform more complex logical operations may be implemented by means of networks of coupled oscillators. We have constructed a Set-Reset Flip Flop circuit, which is capable of storing a bit of information at each forward step. To do so, the system has two inputs (Set, *S*, and Reset, *R*) and two outputs *C* and $\overline{{C}}$, which follow the truth table shown in Figure 5. When neither *S* nor *R* receive an input the output states are *C* = 0 and $\overline{{C}}{=}{1}$. This output state is preserved for the following input step, which is *S* = 0 and *R* = 1. Such feature keeps memory of the previous output state. The following step is *S* = 1 and *R* = 0 which “flips” the output states, being *C* = 1 and $\overline{{C}}{=}{0}$. An undesired state is achieved when both *S* and *R* are 1, which leads to *C* = 0 and $\overline{{C}}{=}{0}$ state. Figure 5A shows the network used to implement such memory: two inputs (*S* and *R*) feed two cortical columns which in turn project unidirectionally toward neural masses *B*_{1} and *A*_{2}, via intermediate cortical columns *I*_{1} and *I*_{2}, respectively. The state of synchronization of nodes *A*_{1} and *A*_{2} will give the output *C*, whereas the state of synchronization of nodes *B*_{1} and *B*_{2} gives $\overline{{C}}$. Figure 5B shows the online implementation of the flip-flop memory, with “flipped” output states and an undesired ${C}{=}\overline{{C}}{=}{0}$ state in which none of the pairs *A*_{1}–*A*_{2} or *B*_{1}–*B*_{2} are synchronized. Note that transients affect the performance of the flip-flop implementation (see third pair of time traces below Figure 5B) but such obstacle may be dependent on initial conditions. Moreover, as in previous implementations, each output state is characterized by distinctive time evolutions, which, on top of memory storage, also give information about which terminal (*S* or *R*) receives the input.

**Figure 5. Implementation of a flip-flop memory**. (**A**) Cartoon depicting the network of cortical columns that fulfill the flip-flop operational gate. *S* and *R* are cortical columns that receive oscillatory inputs, respectively. These columns excite upstream connected columns (*A*_{1} and *B*_{2}) which are bidirectionally coupled to columns *A*_{2} and *B*_{1}, accordingly. The latter receive inputs from two columns that, in turn, receive the output from *A*_{1} and *A*_{2} (output *C*) and *B*_{1} and *B*_{2} (output $\overline{{C}}$.) (**B**) Online implementation of a flip-flop. The state of synchronization is preserved when inputs *S* = 0 and *R* = 0, or *S* = 1 and *R* = 0, giving *C* = 0 and $\overline{{C}}{=}{1}$, but “flips” when *S* = 0 and *R* = 1, giving *C* = 1 and $\overline{{C}}{=}{0}$. A not desired situation is the one in which *C* = 0 and $\overline{{C}}{=}{0}$, i.e., neither *A*_{1} and *A*_{2} or *B*_{1} and *B*_{2} are synchronized. This happens when *S* = 1 and *R* = 1, which fulfills a flip-flop truth table. (**C**) Magnification of the time scale of the output signals for the sequence of five *S*−*R* bit pairs of panel **(B)**. The full scale for each sub panel with a pair of signals is corresponding to 5 s.

## 4. Discussion

In order to address how information can be processed, from the perspective of synchrony, at the mesoscopic scale we have analyzed theoretically a network of coupled neural mass oscillators. This network uses synchronization as the essential ingredient to process the information arriving to/from each of its nodes. We have seen that by interpreting the inputs arriving to the oscillators of the network as “0” or “1” and defining the output of the binary operation in terms of the synchronized state of the two oscillators also as “0” or “1”, several binary logic gates can be constructed. This dependence of the synchronization level of two columns on their stimulation has been observed experimentally, for instance, in the cat visual cortex (Gray et al., 1989). Interestingly, different binary logic gates constructed using the same physiological circuitry result only from changes in the input signals received by the oscillators (e.g., AND, OR, and XOR gates in Figure 2). This rich behavior shown by only two coupled cortical columns may be very fruitful when many other columns are considered. In this sense, the ability to analyze input signals with very different characteristics (average density of spikes, amplitude and frequency of oscillations or noise level) is multiplied by the simple addition of this type of binary logic gates in a network. Nevertheless, this simplistic view may be even more sophisticated when putting the binary motifs together in a larger network. As shown in Figure 3, simply by connecting two different logic gates through a hub may result in a system where the two gates operate in parallel independently of each other (Figure 3A) or operate in a different way (Figure 3B). In this case, outputs may depend on the input of both gates at the same time or on the history of the input states driving the nodes (see examples of both behaviors in Figure 3B). This type of dynamics, in larger networks, makes selectivity of the state in terms of the input protocol even richer than just the repetition of simple binary logic gates in a network. In order to show the generality of this type of networks (and also its robustness in terms of the dynamical oscillators used to build the network) we have constructed several binary logic gates with electronic circuits operating in a chaotic regime. We have shown experimentally that a network built by coupling two of these gates through a hub (using the same simple motif as before) is able to process information as expected (see Figure 4). Finally, we have shown theoretically that by using a network of oscillators we can implement a Set-Reset Flip Flop circuit (Hahnloser et al., 2000), which is an example of another stimulus selector, in this case, that is able to store information.

To conclude, it is worth mentioning that, in this work, we have considered only the simplest interpretation of input and output states (leading to Boolean logic). However, our results may be analyzed in a wider view, for instance, if we explicitly consider the degree of synchronization of the different elements (resulting in fuzzy logic) or if we consider as possible output states all types of synchronization (phase, generalized, lag, complete, …) between the different elements which form the network of oscillators. The fact that we consider only one of the different dynamical characteristics of the system, in our case its degree of synchronization, is a coarse simplification. The dynamical response of the network is not determined only by its degree of synchronization. For instance, the frequencies involved in the dynamics, or the degree of excitation/inhibition segregation, may also inform about the input stimulus characteristics, enlarging in this way the computational capabilities of the system.

## Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

Work supported by the Ministerio de Economia y Competividad (Spain) and FEDER through projects FIS2012-37655-C02-01 and FIS2013-41057-P. JG acknowledges support from the ICREA Academia programme and from the Generalitat de Catalunya (project 2014SGR0974).

## References

Abarbanel, H., Rulkov, N., and Sushchik, M. (1996). Generalized synchronization of chaos: the auxiliary system approach. *Phys. Rev. E. Stat. Phys. Plasmas. Fluids. Relat. Interdiscip. Topics* 53, 4528–4535. doi: 10.1103/PhysRevE.53.4528

Abeles, M. (1991). *Corticonics: Neural Circuits of the Cerebral Cortex, 1st Edn.* Cambridge: Cambridge University Press.

Abeles, M., Hayon, G., and Lehmann, D. (2004). Modeling compositionality by dynamic binding of synfire chains. *J. Comput. Neurosci.* 17, 179–201. doi: 10.1023/B:JCNS.0000037682.18051.5f

Arenas, A., Díaz-Guilera, A., Kurths, J., Moreno, Y., and Zhou, C. (2008). Synchronization in complex networks. *Phys. Rep.* 469, 93–153. doi: 10.1016/j.physrep.2008.09.002

Asai, Y., Guha, A., and Villa, A. E. P. (2008). Deterministic neural dynamics transmitted through neural networks. *Neural Netw.* 21, 799–809. doi: 10.1016/j.neunet.2008.06.014

Asai, Y., and Villa, A. E. P. (2012). Integration and transmission of distributed deterministic neural activity in feed-forward networks. *Brain Res.* 1434, 17–33. doi: 10.1016/j.brainres.2011.10.012

Barardi, A., Malagarriga, D., Sancristobal, B., García-Ojalvo, J., and Pons, A. J. (2014a). Probing scale interaction in brain dynamics through synchronization. *Philos. Trans. R. Soc. B Biol. Sci.* 369:20130533. doi: 10.1098/rstb.2013.0533

Barardi, A., Sancristóbal, B., and Garcia-Ojalvo, J. (2014b). Phase-coherence transitions and communication in the gamma range between delay-coupled neuronal populations. *PLoS Comput. Biol.* 10:e1003723. doi: 10.1371/journal.pcbi.1003723

Başar, E., Başar-Eroğlu, C., Karakaş, S., and Schürmann, M. (2000). Brain oscillations in perception and memory. *Int. J. Psychophysiol.* 35, 95–124. doi: 10.1016/S0167-8760(99)00047-1

Bennett, C. (1973). Logical reversibility of computation. *IBM J. Res. Dev.* 17, 525–532. doi: 10.1147/rd.176.0525

Blumhagen, F., Zhu, P., Shum, J., Schärer, Y., Yaksi, E., Deisseroth, K., et al. (2011). Neuronal filtering of multiplexed odour representations. *Nature* 479, 493–498. doi: 10.1038/nature10633

Boccaletti, S., Kurths, J., Osipov, G., Valladares, D., and Zhou, C. (2002). The synchronization of chaotic systems. *Phys. Rep.* 366, 1–101. doi: 10.1016/S0370-1573(02)00137-0

Börgers, C., Epstein, S., and Kopell, N. J. (2005). Background gamma rhythmicity and attention in cortical local circuits: a computational study. *Proc. Natl. Acad. Sci. U.S.A.* 102, 7002–7007. doi: 10.1073/pnas.0502366102

Börgers, C., Epstein, S., and Kopell, N. J. (2008). Gamma oscillations mediate stimulus competition and attentional selection in a cortical network model. *Proc. Natl. Acad. Sci. U.S.A.* 105, 18023–18028. doi: 10.1073/pnas.0809511105

Börgers, C., and Kopell, N. (2003). Synchronization in networks of excitatory and inhibitory neurons with sparse, random connectivity. *Neural Comput.* 15, 509–538. doi: 10.1162/089976603321192059

Börgers, C., and Kopell, N. (2008). Gamma osillations and stimulus selection. *Neural Comput.* 20, 383–414. doi: 10.1162/neco.2007.07-06-289

Bosman, C. A., Schoffelen, J. M., Brunet, N., Oostenveld, R., Bastos, A. M., Womelsdorf, T., et al. (2012). Attentional stimulus selection through selective synchronization between monkey visual areas. *Neuron* 75, 875–888. doi: 10.1016/j.neuron.2012.06.037

Brown, E., Moehlis, J., and Holmes, P. (2004). On the phase reduction and response dynamics of neural oscillator populations. *Neural Comput.* 16, 673–715. doi: 10.1162/089976604322860668

Buzsáki, G., and Draguhn, A. (2004). Neuronal oscillations in cortical networks. *Science* 304, 1926–1929. doi: 10.1126/science.1099745

Castelo-Branco, M., Goebel, R., Neuenschwander, S., and Singer, W. (2000). Neural synchrony correlates with surface segregation rules. *Nature* 405, 685–689. doi: 10.1038/35015079

Contreras, D., and Steriade, M. (1995). Cellular basis of EEG slow rhythms: a study of dynamic corticothalamic relationships. *J. Neurosci*. 15(1 Pt 2), 604–622.

Crunelli, V., Tóth, T. I., Cope, D. W., Blethyn, K., and Hughes, S. W. (2005). The ‘window’t-type calcium current in brain dynamics of different behavioural states. *J. Physiol.* 562, 121–129. doi: 10.1113/jphysiol.2004.076273

Del Prete, V., Martignon, L., and Villa, A. E. (2004). Detection of syntonies between multiple spike trains using a coarse-grain binarization of spike count distributions. *Network* 15, 13–28. doi: 10.1088/0954-898X/15/1/002

Domedel-Puig, N., Rué, P., Pons, A. J., and García-Ojalvo, J. (2011). Information routing driven by background chatter in a signaling network. *PLoS Comput. Biol.* 7:e1002297. doi: 10.1371/journal.pcbi.1002297

Engel, A. K., Fries, P., and Singer, W. (2001). Dynamic predictions: oscillations and synchrony in top–down processing. *Nat. Rev. Neurosci.* 2, 704–716. doi: 10.1038/35094565

Feinerman, O., Rotem, A., and Moses, E. (2008). Reliable neuronal logic devices from patterned hippocampal cultures. *Nat. Phys.* 4, 967–973. doi: 10.1038/nphys1099

Fries, P. (2009). Neuronal gamma-band synchronization as a fundamental process in cortical computation. *Annu. Rev. Neurosci.* 32, 209–224. doi: 10.1146/annurev.neuro.051508.135603

Fries, P., Reynolds, J. H., Rorie, A. E., and Desimone, R. (2001). Modulation of oscillatory neuronal synchronization by selective visual attention. *Science* 291, 1560–1563. doi: 10.1126/science.1055465

Fries, P., Womelsdorf, T., Oostenveld, R., and Desimone, R. (2008). The effects of visual stimulation and selective visual attention on rhythmic neuronal synchronization in macaque area v4. *J. Neurosci.* 28, 4823–4835. doi: 10.1523/JNEUROSCI.4499-07.2008

García-Ojalvo, J., and Sancho, J. (1999). *Noise in Spatially Extended Systems*. New York, NY: Springer Verlag.

Gray, C. M., Konig, P., Engel, A. K., and Singer, W. (1989). Oscillatory responses in cat visual cortex exhibit inter-columnar synchronization which reflect global stimulus properties. *Nature* 338, 334–337. doi: 10.1038/338334a0

Hahnloser, R. H., Sarpeshkar, R., Mahowald, M. A., Douglas, R. J., and Seung, H. S. (2000). Digital selection and analogue amplification coexist in a cortex-inspired silicon circuit. *Nature* 405, 947–951. doi: 10.1038/35016072

Hoppensteadt, F. C., and Izhikevich, E. M. (2000). Synchronization of laser oscillators, associative memory, and optical neurocomputing. *Phys. Rev. E* 62, 4010. doi: 10.1103/physreve.62.4010

Huang, G., Zhang, D., Meng, J., and Zhu, X. (2011). Interactions between two neural populations: a mechanism of chaos and oscillation in neural mass model. *Neurocomputing* 74, 1026–1034. doi: 10.1016/j.neucom.2010.11.019

Hughes, S. W., Cope, D. W., Blethyn, K. L., and Crunelli, V. (2002). Cellular mechanisms of the slow (<1 Hz) oscillation in thalamocortical neurons *in vitro*. *Neuron* 33, 947–958. doi: 10.1016/S0896-6273(02)00623-2

Iglesias, J., and Villa, A. E. P. (2010). Recurrent spatiotemporal firing patterns in large spiking neural networks with ontogenetic and epigenetic processes. *J. Physiol. Paris* 104, 137–146. doi: 10.1016/j.jphysparis.2009.11.016

Jansen, B. H., and Rit, V. G. (1995). Electroencephalogram and visual evoked potential generation in a mathematical model of coupled cortical columns. *Biol. Cybern.* 73, 357–366. doi: 10.1007/BF00199471

Jansen, B. H., Zouridakis, G., and Brandt, M. E. (1993). A neurophysiologically-based mathematical model of flash visual evoked potentials. *Biol. Cybern.* 68, 275–283. doi: 10.1007/BF00224863

Kennedy, M. P. (1992). Robust OP AMP realization of Chua's circuit. *Frequenz* 46, 66–80. doi: 10.1515/FREQ.1992.46.3-4.66

Kopell, N., Ermentrout, G. B., Whittington, M. A., and Traub, R. D. (2000). Gamma rhythms and beta rhythms have different synchronization properties. *Proc. Natl. Acad. Sci. U.S.A.* 97, 1867–1872. doi: 10.1073/pnas.97.4.1867

Lachaux, J.-P., Rodriguez, E., Martinerie, J., and Varela, F. J. (1999). Measuring phase synchrony in brain signals. *Hum. Brain Mapp.* 8, 194–208. doi: 10.1002/(SICI)1097-0193(1999)8:4<194::AID-HBM4>3.0.CO;2-C

Laurent, G., Stopfer, M., Friedrich, R. W., Rabinovich, M. I., Volkovskii, A., and Abarbanel, H. D. (2001). Odor encoding as an active, dynamical process: experiments, computation, and theory. *Ann. Rev. Neurosci.* 24, 263–297. doi: 10.1146/annurev.neuro.24.1.263

Levy, J., Vidal, J. R., Oostenveld, R., FitzPatrick, I., Démonet, J.-F., and Fries, P. (2013). Alpha-band suppression in the visual word form area as a functional bottleneck to consciousness. *NeuroImage* 78, 33–45. doi: 10.1016/j.neuroimage.2013.04.020

Li, C., and Chen, G. (2004). Synchronization in general complex dynamical networks with coupling delays. *Phys. A Stat. Mech. Appl.* 343, 263–278. doi: 10.1016/j.physa.2004.05.058

Madan, R. (1993). *Chua's Circuit: A Paradigm for Chaos*. River Edge, NJ: World Scientific Publishing Company.

Malagarriga, D., Villa, A. E. P., García-Ojalvo, J., and Pons, A. J. (2014). “Excitation/inhibition patterns in a system of coupled cortical columns,” in *Artificial Neural Networks and Machine Learning–ICANN 2014*, eds S. Wermter, C. Weber, W. Duch, T. Honkela, P. Koprinkova-Hristova, S. Magg, G. Palm, and A. E. P. Villa (Hamburg: Springer), 651–658.

Malagarriga, D., Villa, A. E. P., Garcia-Ojalvo, J., and Pons, A. J. (2015). Mesoscopic segregation of excitation and inhibition in a brain network model. *PLoS Comput. Biol.* 11:e1004007. doi: 10.1371/journal.pcbi.1004007

Maris, E., Womelsdorf, T., Desimone, R., and Fries, P. (2013). Rhythmic neuronal synchronization in visual cortex entails spatial phase relation diversity that is modulated by stimulation and attention. *NeuroImage* 74, 99–116. doi: 10.1016/j.neuroimage.2013.02.007

McCulloch, W. S., and Pitts, W. (1990). A logical calculus of the ideas immanent in nervous activity. 1943. *Bull. Math. Biol.* 52, 99–115. discussion 73–97. doi: 10.1007/BF02459570

Melloni, L., Molina, C., Pena, M., Torres, D., Singer, W., and Rodriguez, E. (2007). Synchronization of neural activity across cortical areas correlates with conscious perception. *J. Neurosci.* 27, 2858–2865. doi: 10.1523/JNEUROSCI.4623-06.2007

Menon, S., and Sinha, S. (2014). “‘Defective logic: using spatiotemporal patterns in coupled relaxation oscillator arrays for computation,” in *International Conference on Signal Processing and Communications (SPCOM), 2014* (Bangalore: IEEE), 1–6.

Mormann, F., Lehnertz, K., David, P., and Elger, C. E. (2000). Mean phase coherence as a measure for phase synchronization and its application to the {EEG} of epilepsy patients. *Phys. D* 144, 358–369. doi: 10.1016/S0167-2789(00)00087-7

Moskalenko, O. I., Koronovskii, A. A., Hramov, A. E., and Boccaletti, S. (2012). Generalized synchronization in mutually coupled oscillators and complex networks. *Phys. Rev. E* 86:036216. doi: 10.1103/PhysRevE.86.036216

Nikolić, D., Fries, P., and Singer, W. (2013). Gamma oscillations: precise temporal coordination without a metronome. *Trends Cogn. Sci.* 17, 54–55. doi: 10.1016/j.tics.2012.12.003

Phillips, J. M., Vinck, M., Everling, S., and Womelsdorf, T. (2013). A long-range fronto-parietal 5- to 10-hz network predicts “top-down” controlled guidance in a task-switch paradigm. *Cerebral Cortex* 24, 1996–2008. doi: 10.1093/cercor/bht050

Pikovsky, A., Rosenblum, M., and Kurths, J. (2003). *Synchronization: A Universal Concept in Nonlinear Sciences, Vol. 12.* Cambridge: Cambridge University Press.

Pons, A. J., Cantero, J. L., Atienza, M., and Garcia-Ojalvo, J. (2010). Relating structural and functional anomalous connectivity in the aging brain via neural mass modeling. *NeuroImage* 52, 848–861. doi: 10.1016/j.neuroimage.2009.12.105

Quian Quiroga, R., Kraskov, A., Kreuz, T., and Grassberger, P. (2002). Performance of different synchronization measures in real data: a case study on electroencephalographic signals. *Phys. Rev. E* 65, 041903. doi: 10.1103/physreve.65.041903

Rabinovich, M. I., and Varona, P. (2011). Robust transient dynamics and brain functions. *Front. Comput. Neurosci.* 5:24. doi: 10.3389/fncom.2011.00024

Rodriguez, E., George, N., Lachaux, J., Martinerie, J. P., Renault, B., and Varela, F. (1999). Perception's shadow: long-distance synchronization of human brain activity. *Nature* 397, 430–433. doi: 10.1038/17120

Rosenblum, M. G., Pikovsky, A. S., and Kurths, J. (1996). Phase synchronization of chaotic oscillators. *Phys. Rev. Lett.* 76, 1804–1807. doi: 10.1103/PhysRevLett.76.1804

Rué, P., Domedel-Puig, N., Garcia-Ojalvo, J., and Pons, A. J. (2012). Integration of cellular signals in chattering environments. *Prog. Biophys. Mol. Biol.* 110, 106–112. doi: 10.1016/j.pbiomolbio.2012.05.003

Singer, W. (1993). Synchronization of cortical activity and its putative role in information processing and learning. *Ann. Rev. Physiol.* 55, 349–374. doi: 10.1146/annurev.ph.55.030193.002025

Singer, W. (1999). Neuronal synchrony: a versatile code for the definition of relations. *Neuron* 24, 49–65. doi: 10.1016/S0896-6273(00)80821-1

Skarda, C. A., and Freeman, W. J. (1987). How brains make chaos in order to make sense of the world. *Behav. Brain Sci.* 10, 161–195. doi: 10.1017/S0140525X00047336

Spiegler, A., Knösche, T. R., Schwab, K., Haueisen, J., and Atay, F. M. (2011). Modeling brain resonance phenomena using a neural mass model. *PLoS Comput. Biol.* 7:e1002298. doi: 10.1371/journal.pcbi.1002298

Stam, C. (2005). Nonlinear dynamical analysis of eeg and meg: review of an emerging field. *Clin. Neurophysiol.* 116, 2266–2301. doi: 10.1016/j.clinph.2005.06.011

Stam, C. J., and de Bruin, E. A. (2004). Scale-free dynamics of global functional connectivity in the human brain. *Hum. Brain Mapp.* 22, 97–109. doi: 10.1002/hbm.20016

Stopfer, M., Bhagavan, S., Smith, B., and Laurent, G. (1997). Impaired odour discrimination on desynchronization of odour-encoding neural assemblies. *Nature* 390, 70–74.

Toral, R. and Colet, P. (2014). *Stochastic Numerical Methods: An Introduction for Students and Scientists*. Wenheim: Wiley-VCH.

Tort, A. B. L., Komorowski, R. W., Manns, J. R., Kopell, N. J., and Eichenbaum, H. (2009). Theta–gamma coupling increases during the learning of item–context associations. *Proc. Natl. Acad. Sci. U.S.A.* 106, 20942–20947. doi: 10.1073/pnas.0911331106

Vardi, R., Guberman, S., Goldental, A., and Kanter, I. (2013). An experimental evidence-based computational paradigm for new logic-gates in neuronal activity. *Europhys. Lett.* 103, 66001. doi: 10.1209/0295-5075/103/66001

Varela, F., Lachaux, J. P., Rodriguez, E., and Martinerie, J. (2001). The brainweb: phase synchronization and large-scale integration. *Nat. Rev. Neurosci.* 2, 229–239. doi: 10.1038/35067550

Villa, A. E. P., and Tetko, I. V. (2010). Cross-frequency coupling in mesiotemporal eeg recordings of epileptic patients. *J. Physiol. Paris* 104, 197–202. doi: 10.1016/j.jphysparis.2009.11.024

Wagemakers, A., Buldú, J. M., and Sanjuán, M. A. F. (2007). Isochronous synchronization in mutually coupled chaotic circuits. *Chaos* 17:023128. doi: 10.1063/1.2737820

Ward, L. M. (2003). Synchronous neural oscillations and cognitive processes. *Trends Cogn. Sci.* 7, 553 – 559. doi: 10.1016/j.tics.2003.10.012

Wolf, F., and Geisel, T. (2008). Neurophysics: logic gates come to life. *Nat. Phys.* 4, 905–906. doi: 10.1038/nphys1142

Womelsdorf, T., and Fries, P. (2007). The role of neuronal synchronization in selective attention. *Curr. Opin. Neurobiol.* 17, 154–160. doi: 10.1016/j.conb.2007.02.002

Womelsdorf, T., Fries, P., Mitra, P. P., and Desimone, R. (2006). Gamma-band synchronization in visual cortex predicts speed of change detection. *Nature* 439, 733–736. doi: 10.1038/nature04258

Womelsdorf, T., Schoffelen, J.-M., Oostenveld, R., Singer, W., Desimone, R., Engel, A. K., et al. (2007). Modulation of neuronal interactions through neuronal synchronization. *Science* 316, 1609–1612. doi: 10.1126/science.1139597

Womelsdorf, T., Valiante, T. A., Sahin, N. T., Miller, K. J., and Tiesinga, P. (2014). Dynamic circuit motifs underlying rhythmic gain control, gating and integration. *Nat. Neurosci.* 17, 1031–1039. doi: 10.1038/nn.3764

Wulff, P., Ponomarenko, A. A., Bartos, M., Korotkova, T. M., Fuchs, E. C., Bähner, F., et al. (2009). Hippocampal theta rhythm and its coupling with gamma oscillations require fast inhibition onto parvalbumin-positive interneurons. *Proc. Natl. Acad. Sci. U.S.A.* 106, 3561–3566. doi: 10.1073/pnas.0813176106

Xu, D., Principe, J. C., and Harris, J. G. (2004). “Logic computation using coupled neural oscillators,” in *International Symposium on Circuits and Systems, 2004. ISCAS'04. Proceedings of the 2004*, Vol. 5 (Vancouver: IEEE), V788–V791.

Zanette, D. H. (2004). Propagation of small perturbations in synchronized oscillator networks. *Europhys. Lett.* 68, 356. doi: 10.1209/epl/i2004-10218-2

Keywords: synchronization, neural mass, Chua oscillators, complex networks, information processing, logic gate

Citation: Malagarriga D, García-Vellisca MA, Villa AEP, Buldú JM, García-Ojalvo J and Pons AJ (2015) Synchronization-based computation through networks of coupled oscillators. *Front. Comput. Neurosci*. 9:97. doi: 10.3389/fncom.2015.00097

Received: 08 April 2015; Accepted: 14 July 2015;

Published: 04 August 2015.

Edited by:

Andre Longtin, University of Ottawa, CanadaReviewed by:

Carmen Castro Canavier, LSU Health Sciences Center New Orleans, USAGianluigi Mongillo, Paris Descartes University, France

Copyright © 2015 Malagarriga, García-Vellisca, Villa, Buldú, García-Ojalvo and Pons. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Daniel Malagarriga, Departament de Física i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Edifici Gaia, Rambla Sant Nebridi 22, 08222 Terrassa, Spain, daniel.malagarriga@upc.edu