% Template article for preprint document class `elsart'
% SP 2006/04/26

\documentclass{elsart}

% Use the option doublespacing or reviewcopy to obtain double line spacing
 %\documentclass[doublespacing]{elsart}

% if you use PostScript figures in your article
% use the graphics package for simple commands
% \usepackage{graphics}
% or use the graphicx package for more complicated commands
% \usepackage{graphicx}
% or use the epsfig package if you prefer to use the old commands
% \usepackage{epsfig}

% The amssymb package provides various useful mathematical symbols
\usepackage{amssymb}
\usepackage{amsmath}
\usepackage{epsfig}
% The lineno packages adds line numbers. Start line numbering with
% \begin{linenumbers}, end it with \end{linenumbers}. Or switch it on
% for the whole article with \linenumbers.
% \usepackage{lineno}

% \linenumbers

\usepackage{stmaryrd}
%my definitions
\def\norm#1{\|#1\|} 
\newcommand{\pop}{\mathcal{I}}
\newcommand{\stat}{\mathcal{S}}
%\newcommand{\traj}[1]{\mathcal{T}^{ #1} \left(I ,\mathcal{E}_b \right)}
\newcommand{\traj}{\mathcal{T}}
\newcommand{\eps}{\varepsilon}
\newcommand{\glo}[1]{#1}
\newcommand{\ind}{\mathds{1}}
\newcommand{\To}{\longrightarrow}
\newcommand{\h}{\mathcal{H}}
\newcommand{\A}{\mathcal{A}}
\newcommand{\F}{\mathcal{F}}
\newcommand{\J}{\mathcal{J}}
\newcommand{\M}{\mathcal{M}}
\newcommand{\W}{\mathcal{W}}
\newcommand{\X}{\mathcal{X}}
\newcommand{\BOP}{\mathbf{B}}
\newcommand{\BH}{\mathbf{B}(\mathcal{H})}
\newcommand{\KH}{\mathcal{K}(\mathcal{H})}
\newcommand{\R}{\mathbb{R}}
\newcommand{\pr}{\mathbb{P}}
\newcommand{\N}{\mathbb{N}}
\newcommand{\Z}{\mathbb{Z}}
\newcommand{\Q}{\mathbb{Q}}
%\newcommand{\E}{\mathbb{E}}
\newcommand{\T}{\mathbb{T}}
\newcommand{\Complex}{\mathbb{C}}
\newcommand{\Field}{\mathbb{F}}
\newcommand{\RPlus}{\Real^{+}}
\newcommand{\Poly}{\mathcal{P}(E)}
\newcommand{\EssD}{\mathcal{D}}
\newcommand{\Lom}{\mathcal{L}}
\newcommand{\abs}[1]{\left\vert#1\right\vert}
\newcommand{\set}[1]{\left\{#1\right\}}
\newcommand{\seq}[1]{\left<#1\right>}
%\newcommand{\norm}[1]{\left\Vert#1\right\Vert}
\newcommand{\essnorm}[1]{\norm{#1}_{\ess}}
\newcommand{\leftexp}[2]{{\, \vphantom{#2}}^{#1}{#2}}
\newcommand{\leps}[1]{\leftexp{\eps}{#1}}
\newcommand{\ec}[2]{\E \left[ \left. #1 \right| #2 \right]}
%\renewcommand{\proofname}{D\'emonstration}
%\newcommand{\eps}{\varepsilon}
%\newcommand{\leftexp}[2]{{\, \vphantom{#2}}^{#1}{#2}}
%\newcommand{\leps}[1]{\leftexp{\eps}{#1}}
\newcommand{\ieg}{ \left[\kern-0.15em\left[ }
\newcommand{\ied}{ \right]\kern-0.15em\right] }
% \nbOne : 1I : symbol one
\def\nbOne{{\mathchoice {\rm 1\mskip-4mu l} {\rm 1\mskip-4mu l}
{\rm 1\mskip-4.5mu l} {\rm 1\mskip-5mu l}}}
\newcommand{\app}[5]{
\begin{array}{rccl}
 #1 : 	& #2  	& \longrightarrow & #3 \\
 			&	#4 	& \longmapsto 		& #5
\end{array} 
}

\begin{document}

\begin{frontmatter}

% Title, authors and addresses

% use the thanksref command within \title, \author or \address for footnotes;
% use the corauthref command within \author for corresponding author footnotes;
% use the ead command for the email address,
% and the form \ead[url] for the home page:
% \title{Title\thanksref{label1}}
% \thanks[label1]{}
% \author{Name\corauthref{cor1}\thanksref{label2}}
% \ead{email address}
% \ead[url]{home page}
% \thanks[label2]{}
% \corauth[cor1]{}
% \address{Address\thanksref{label3}}
% \thanks[label3]{}

\title{Frequencies locking in countable cellular systems, localization of quasi-periodic solutions of autonomous differential systems}

% use optional labels to link authors explicitly to addresses:
% \author[label1,label2]{}
% \address[label1]{}
% \address[label2]{}

\author{Laurent Gaubert}

\address{Centre Europ\'een de R\'ealit\'e Virtuelle, 
LISYC EA3883 UBO/ENIB,
25 rue Claude Chappe, 29280 Plouzan\'e, France}

\begin{abstract}
We address the question of frequencies locking in coupled differential systems, related to the existence of quasi-periodic solutions of differential systems. Our tool is what we call ``cellular systems'', quite general as it deals with countable number of coupled systems in some general Banach spaces. Moreover, the inner dynamics of each subsystem does not have to be specified. We reach some general results about how the frequencies locking phenomenon is related to the structure of the coupling map, and therefore about the localization of quasi-periodic solutions of some differential systems that may be seen as cellular systems. This paper gives some explanations about how and why synchronized behaviors naturally occur it a wide variety of complex systems.
\end{abstract}

\begin{keyword}
% keywords here, in the form: keyword \sep keyword
coupled differential systems \sep synchronization \sep frequencies locking
% PACS codes here, in the form: \PACS code \sep code
\PACS 
\end{keyword}
\end{frontmatter}
%\tableofcontents
\section{Introduction}
\label{sec:introduction}
% Always give a unique label
% and use \ref{<label>} for cross-references
% and \cite{<label>} for bibliographic references
% use \sectionmark{}
% to alter or adjust the section heading in the running head

Synchronization is an extremely important and interesting emergent property of complex systems. The first example found in literature goes back to the 17th century with Christiaan Huygens' work \cite{Huge73,Benn02}. This kind of emergent behavior can be found in artificial systems as well as in natural ones and at many scales (from cell to whole ecological systems). Biology abounds with periodic and synchronized phenomena and the work of Ilya Prigogine shows that such behaviors arise within specific conditions: a dissipative structure generally associated to a nonlinear dynamics \cite{Prig67}. Biological systems are open, they evolve far from thermodynamic equilibrium and are subject to numerous regulating processes, leading to highly nonlinear dynamics. Therefore periodic behaviors appear (with or without synchronization) at any scale \cite{Prig69}. More generally, life itself is governed by circadian rhythms \cite{Gonz04}. Those phenomena are as much attractive as they are often spectacular: from cicada populations that appear spontaneously every ten or thirteen years \cite{Hopp76} or networks of heart cells that beat together \cite{Mich87} to huge swarms in which fireflies, gathered in a same tree, flash simultaneously \cite{Buck88}. This synchronization phenomenon occupies a privileged position among emergent collective phenomena because of its various applications in neuroscience, ecology, earth Science, for instance \cite{Winf90,Stro03,Manr04}, as well as in the field of coupled dynamical systems, especially through the notion of chaotic systems' synchronization \cite{Peco97,Goed03} and the study of coupled-oscillators \cite{Kope00}. This wide source of examples leads the field of research to be highly interdisciplinary, from pure theory to concrete applications and experimentations.

The classical concept of synchronization is related to the locking of the basic frequencies and instantaneous phases of regular oscillations. One of the most successful attempts to explore this emergent property is due to Kuramoto \cite{Kura75,Kura84}. As in Kuramoto's work, those questions are usually addressed by studying specific kinds of coupled systems (see for instance \cite{Erme86,Ren00,Goel02}). Using all the classical methods available in the field of dynamical systems, researchers study specific trajectories of those systems in order to get information on possible attracting synchronized state \cite{Wu94,Kope00,Ren00,Piko01,Goel02,Izhi03}. 

The starting point of this work was the following question : ``Why synchronization is such a widely present phenomena ?''
In order to give some mathematical answer to this question, the first step is to build a model of coupled systems that is biologically inspired. This is what is done in the second section, after having described some basic material, we define what we call cellular systems and cellular coupler. If one would summarize the specificities of cellular system, one could say that each cell (subsystem) of a cellular system receives information from the whole population (the coupled system) according to some constraints:
\begin{itemize}
	\item a cell has access to linear transformations of all the others cell's states
	\item the way this information is gathered depends (not linearly) on the cell's state itself  
\end{itemize}
In other words, a cell interprets its own environment via the states of the whole population and according to its own state.

It's a bit surprising that despite this model arising very naturally, it gives a good framework to address the main question. Indeed, in the third section we expose a localization result concerning periodic trajectories of cellular systems, according to some sub-periods dependencies. It exhibits some links between the coupler's properties and the structure of periodic trajectories. 

The fourth section gives some example of general results that may be proved using the localization lemma. Moreover, it goes out of the scope of coupled systems as, synchronization is strongly related to the more abstract field of dynamical systems. If one thinks about presence of regular attractors (in opposition with strange attractor) in a differential system, one may for example classify those as:
\begin{itemize}
	\item point attractor
	\item limit cycle
	\item limit torus
\end{itemize}
Those attractors can be related to coupled systems in an obvious way: roughly speaking, a point attractor may be seen as a solution of coupled systems for which each of the subsystems has a constant behavior. Similarly, a limit cycle may be thought as the situation where every subsystem oscillate, all frequencies among the whole system being locked. A limit torus is a similar situation which differs from the previous one by the fact that the frequencies are not locked (non commensurable periods of a quasi-periodic solution of the whole coupled system). Hence, the three previous cases may be translated into the coupled dynamical systems context:
\begin{itemize}
	\item point attractor $\leftrightarrow$ constants trajectories of each subsystem
	\item limit cycle $\leftrightarrow$ periodic sub-trajectories with locked frequencies
	\item limit torus $\leftrightarrow$ periodic sub-trajectories with unlocked frequencies
\end{itemize}

Therefore, we deduce some results about the localization of solutions of the third type, quasi-periodic solutions, using the point of view of coupled dynamical systems. The results of this third section may help to understand why the second case is the most observed in natural systems, which may be seen as coupled dynamical systems, at many levels. Indeed, the section ends with a sketch of how the cellular systems point of view may be applied to a wide class of differential systems in order to systematically address those questions with algebraic tools.

\section{Basic material and notations}
\label{sec:basics}
As our model is inspired by cellular tissues, some terms clearly come from the vocabulary used to describe those kinds of complex systems. 

\subsection{Model of population behavior}
\label{sec:ModelOfPopulation}
Here are the basic compounds and notations of our model: 

A \emph{population} $\pop$ is a countable set, so we may consider it as a subset $\pop \subset \N$. Moreover, as it's only the cardinality of $\pop$ that's important, $\pop$ may be chosen as an interval of integer. Elements of $\pop$ are called \emph{cells}. \\
	
We suppose that the systems we want to study are valued in some Banach spaces. Thus, for any $i \in \pop$, $(E_i,\|.\|_i)$ is a Banach space, and the \emph{state space} of $\pop$ is the vector space $\displaystyle \mathcal{S}=\prod_{i \in I}E_i$.\\

We will sometimes identify $E_i$ with 
$$ \prod_{j < i}\{0\} \times E_i \times \prod_{j > i}\{0\} \subset \stat$$
and then consider it as a subspace of $\stat$.\\

We denote $\mathcal{S}_b$ the space of uniformly bounded states:
$$
\mathcal{S}_b = \left\{x \in \mathcal{S}, \; \sup_{i \in \pop} \| x_i \|_i < \infty \right\}
$$
This subspace will sometimes be useful as, embodied with the norm $\|x\|_{\infty}=\sup_{i \in \pop} \| x_i \|_i$, it's a Banach space, allowing the classic Picard�Lindel�f theorem to be valid.\\	
	 
Given an interval $\Omega \subset \R$, a \emph{trajectory} $x$ of $\pop$ is an element of $\mathcal{C}^{\infty}\left(\Omega , \stat \right)$. Such $x$ is then described by a family of $\mathcal{C}^{\infty}$ applications $(x_i)_{i \in \pop}$ such that $\forall i \in \pop$:
	$$
	\app{x_i}{\Omega }{E_i}{t}{x_i(t)}
	$$
	The space of trajectories on $I$ is denoted $\traj$.\\

A \emph{period} on $\pop$ is a map $\tau : \pop \to \R_+^*$. A trajectory $x \in \mathcal{T}$ is said to be \emph{$\tau$-periodic} if for any $i \in \pop$, $x_i$ is $\tau(i)$-periodic and non constant. $\tau(i)$ is then said to be \emph{a period of the cell} $i$. The space of such trajectories is written $\mathcal{T}_{\tau}$.


Each cell $i$ is supposed to behave according to an autonomous differential system given by a vector field $F_i: E_i \to E_i$. Thus, given a family of functions $\{F_i\}_{i \in \pop}$ we define the vector field $F_{\pop}$ on $\mathcal{S}$:
	$$
	\app{F_{\pop}}{\mathcal{S}}{\mathcal{S}}{x}{F_{\pop}(x)}
	$$ 
	with, for any $i \in \pop$:
	$$
	\left[ F_{\pop}(x) \right]_i=F_i(x_i)
	$$
	
	%$F_{\pop}$ is said to be $\mathcal{C}^r$ if every $F_p$ is $\mathcal{C}^r$ and if for any bounded subset $B$ of $\R^n$ we have:
	%$$
	%\sup_{x \in B \times p \in \pop}\norm{dF_p(x)} < \infty
	%$$
	\begin{rem}
	The definition of periodic trajectory handle both classical concepts of periodic and quasi-periodic solutions of a differential system. From the point of view of coupled systems, it describes the situation in which each subsystem of the whole system oscillates.	We stress the point that a period of a periodic trajectory needs not to be a minimal period ($\tau(i)$ isn't necessarily a generator of the group of periods of $x_i$). Nevertheless, our definition of $\mathcal{T}_{\tau}$ avoid any trajectory which contains some constant component (none of the $x_i$ can be a constant map) as they may be seen as degenerate (localized into an hyperplane of $\mathcal{S}$).
	\end{rem}  
 
We recall that a (finite) subset $\left\{ \tau_1, \ldots , \tau_k \right\}$ of $\R$ is said to be \textbf{rationally dependent} if there exists some integers $l_1, \ldots , l_k$ non all zero and such that:
$$
l_1\tau_1 + \ldots + l_k \tau_k = 0
$$
Then there exists a unique lowest common multiple (\emph{lcm}) $\tau_0$ for which there exits $n_1, \ldots ,n_k$ such that:
$$
n_1 \tau_1= \ldots = n_k \tau_k=\tau_0
$$
An infinite set of real numbers is said to be rationally dependent if any finite subset is rationally dependent.
 
Now, any period $\tau$ on $\pop$ (or, equivalently, any periodic trajectory) defines a equivalence relation on $\pop$ as:
$$
i \sim j \; \Leftrightarrow \; \{\tau(i),\tau(j) \} \textnormal{ is a dependent set}
$$ 
Hence we may consider the (countable) partition $\pop(\tau)$ of $\pop$ into equivalence classes ($K$ countable):
$$
\pop(\tau)=\left\{\pop_k \right\}_{k \in K} 
$$

	
\subsection{Cellular coupler and cellular systems}
\label{sec:CouplerOnPop}
In this section we build what we call \emph{cellular systems} by way of \emph{cellular coupler}. Most of the works always deal with a specific way of coupling dynamical systems: one adds a quantity (that models interactions between subsystems) to the derivative of the systems. This leads to equations with the following typical shape (here, there are only two coupled systems):

\begin{align*}
x_1'(t)=F \big(x_1(t) \big)+  G_1 \big(x_1(t),x_2(t)\big) \\
x_2'(t)=F \big(x_2(t) \big)+  G_2 \big(x_1(t),x_2(t)\big) 
\end{align*}
The functions $G_1$ and $G_2$ are the coupling functions. The problem is then restated in terms of phase-shift variables and efforts are made to detect stable states and to prove their stability.

Our approach to coupling is different. We study exclusively a way of coupling where the exchanges are made on the current state of the system. This means that the coupling quantity applies inside the map $F$, which leads us to the following type of equation:

\begin{align}
\label{eq:coupling_source} x_1'(t)=F \big(x_1(t)+  H_1(x_1(t),x_2(t)) \big) \\
\nonumber x_2'(t)=F \big(x_2(t)+  H_2(x_1(t),x_2(t)) \big) 
\end{align}
\begin{rem}
We stress the point that those two different ways of handling coupled systems are quite equivalent in most cases. Indeed, starting with the first two equations, as soon as $G_1$ and $G_2$ stay in the range of $F$ (which is likely if the coupling functions are small), we can rewrite them in the second shape involving $H_1$ and $H_2$.
\end{rem}
The last kind of coupled systems is sometimes studied (for instance in \cite{Izhi03}) but never broadly (indeed, if one wants some quantitative results about convergence of trajectories, one must work with specific equations and dynamical systems). Even in a few papers that are quite general (as the very interesting \cite{Stew02}) some strong assumptions are made (in \cite{Stew02} authors deal with symmetric periodic solutions). The kind of coupled systems we handle are a generalization of the one describe in equation (\ref{eq:coupling_source}). Its general shape is:
$$
x_i'(t)=F_i \left( \sum_{j \in \pop} c_{ij}(x_i(t))  x_j(t) \right)
$$
Each cell $i \in \pop$ holds it's own differential system represented by a map $F_i$ (hence, all the dynamical systems are not forcibly identical nor have the same shape, nor that they are weakly coupled (as in the classical paper of Art Winfree \cite{Winf67}). A cell $i$ ``interprets'' it's own environment by mean of the functions $c_{ij}$.
 
Before giving the exact definition of a cellular coupler, we recall that, as $\mathcal{S}$ may be seen as a module on the ring $\displaystyle \prod_{i \in \pop} \mathcal{L}(E_i)$, $\mathcal{L}(\mathcal{S})$ has to be understood as the space of linear operators on $\mathcal{S}$ with coefficients in the $\mathcal{L}(E_i,E_j)$. Any $M \in \mathcal{L}(\mathcal{S})$ may then be written as an infinite (if $\pop$ isn't finite) matrix:
$$
M=\left[m_{ij} \right]_{(i,j) \in I^2}, \; m_{ij} \in \mathcal{L}(E_j,E_i)
$$
In this context, here is the definition of a cellular coupler on $\pop$: 
\begin{defn}
\label{def:coupler}
A \emph{cellular coupling map} on $\pop$ is a map $c$:
$$
\app{c}{\mathcal{S}}{\mathcal{L}(\mathcal{S})}{x}{c(x)}
$$
such that the matrix $\left[c_{ij} \right]_{(i,j) \in \pop^2}$ satisfies:
\begin{enumerate}

	\item $\displaystyle \forall (i,j) \in \pop^2, \;  \forall x \in \mathcal{S}, \; c_{ij}(x) \textnormal{ depends only on } x_i, \textnormal{ so that we may consider it as a map } \\ c_{ij} : E_i \to \mathcal{L}(E_j,E_i)$\\
	
	\item $ \displaystyle \forall (i,j) \in \pop^2, \; c_{ij}$ has continuous values: $c_{ij} : E_i \to \mathcal{L}_c(E_j,E_i)$\\ \textbf{a garder ?}
	
	\item $ \displaystyle \forall i \in \pop, \; \forall x_i \in E_i,\; \sum_{j \in \pop} \|c_{ij}(x_i) \|_i < + \infty$\\ 
\end{enumerate}
Then, $c$ defines a \emph{cellular coupler} $\widetilde{c}$ on $\pop$ in the following way:
$$
\app{\widetilde{c}}{\mathcal{S}}{\mathcal{S}}{x}{c(x).x}
$$
\end{defn}


In other words (for the sake of simplicity, we only take examples with a finite population), for any $x \in \mathcal{S}$, the matrix $c(x)$ has the following shape:
$$
c(x)=\left[\begin{array}{ccc} c_{11}(x_1) &  \cdots & c_{1k}(x_1) \\
												\vdots & \ddots & \vdots \\
												c_{k1}(x_k)  & \cdots & c_{kk}(x_k)\end{array} \right]
												\in \mathcal{L}(\mathcal{S})
$$
And then :
$$
\widetilde{c}(x)=c(x).x=
c(x)=\left[\begin{array}{c} c_{11}(x_1).x_1+ \ldots + c_{1k}(x_1).x_k\\
												\vdots \\
												c_{k1}(x_k).x_1+ \ldots + c_{kk}(x_k).x_k \end{array} \right]
												\in \stat
$$
Now we can define a cellular system:
\begin{defn}
Let $F_{\pop}$ be a vector field on $\stat$ given by a family $\{F_i\}_{i \in \pop}$ of vector fields on the $E_i$. Let $\widetilde{c}$ a cellular coupler on $\pop$. $(\pop,F_{\pop},\widetilde{c})$ is called a \emph{cellular system}. A trajectory of this system is a trajectory $x \in \mathcal{T}$ that satisfies:
$$
x'=F_{\pop} \circ \widetilde{c}(x)=F_{\pop} \big(c(x).x \big)
$$
in other words:
$$
\forall i \in \pop, \; \forall t \in \Omega, \; x_i'(t)=F_i\left(\sum_{j \in \pop} c_{ij}(x_i(t)).x_j(t)\right)
$$
\end{defn}
This equation may be naturally interpreted in biological terms: the cell $i$ behaves according to a mean of the states of all other cells $x_j$, but only its state defines how this mean is computed (the cell interprets its own environment), and this link \emph{state$\leftrightarrow$ interpreting function} has no reason to be linear in $x_i$.\\

\begin{rem}
In order to avoid any confusion, we stress the differences between trajectory and solution regarding periodic behaviors. In this paper, periodic trajectory has a specific meaning related to the cells. A periodic trajectory of a cellular system is a trajectory for which each cell has a periodic behavior. From the classic point of view of differential equations, periodic and quasi-periodic solutions of the cellular system are periodic trajectories. In the case of a periodic solution, $\tau(\pop)$ admit a \emph{lcm}, which false in the case of a quasi-periodic solution.   
\end{rem}
In the next section we start of by exposing algebraic links between a cellular coupler and a periodic trajectory, and then we turn to our localization lemma.

\section{Localisation lemma}
\label{sec:PropertiesOfACouplerOnPop}
 
If $M$ is a matrix indexed on $\pop^2$, and if $J \subset \pop$, we write $\overline{J} =\pop-J $ and we define $M^J$ as the matrix:
$$
M = (m_{ij})_{(i,j) \in J \times \overline{J}}
$$
For $x \in \stat$ (resp. $x \in \traj$) we denote $x^J$ the vector (resp. the map) $[x_i]_{i \in J}$ (see figure \ref{fig:matrice_extraite}).\\
 
\begin{figure}[ht]
	\centering
		\includegraphics[width=200pt]{matrice_extraite.eps}
	\caption{Matrix and vectors associated to a subset $J$ of $\pop$.}
	\label{fig:matrice_extraite}
\end{figure}
If $\pop(\tau)=\{I_1,...,I_K\}$ is a partition of $\pop$, we define the matrix $M^{\pop(\tau)}$ as (see figure \ref{fig:matrice_extraite_partition}):
$$
m^{\pop(\tau)}_{ij}=\left\{ \begin{array}{ll} 0 & \textnormal{ if } (i,j) \in I_1 \times I_1 \cup \ldots \cup I_K \times I_K \\
																						m_{ij} & \textnormal{ if not } \end{array} \right.
$$
\begin{figure}[ht]
	\centering
		\includegraphics[width=200pt]{matrice_extraite_partition.eps}
	\caption{Matrix associated to a partition of $\pop$.}
	\label{fig:matrice_extraite_partition}
\end{figure}

We can then go further and begin to work on the heart of our matter. The forthcoming result that can be used in many ways and generalized as, for the sake of simplicity, we did not use the weakest assumptions under which it holds (for example, the series convergence in the proof can be insured in many other contexts).
\begin{lem}
\label{lem:localization}
Let $\left(\pop,F_{\pop},\widetilde{c} \right)$ be a cellular system and $\tau$ a period on $\pop$. Let $U \subset \mathcal{S}$ on which $F_{\pop}$ is injective.
If $x \in \mathcal{T}^{\tau}$ is a periodic trajectory of cellular system that satisfies:
\begin{enumerate}
	\item $x(\Omega) \subset \stat_b $ ;
	\item $\widetilde{c}(x)(\Omega) \subset U$
\end{enumerate}
then there exists $b \in \mathcal{S}_b$ such that for any $t \in \Omega$:
$$
 x \left(t \right) - b \in   \ker \left( c(x(t))^{\pop(\tau)}\right)
$$
\end{lem}
\begin{rem}Note that the first condition on $x$ is useless if $\pop$ is finite. \end{rem}

The previous result is not very practical as it involves the trajectory $x$ itself, which is unknown. As there is no ambiguity, we define the kernel of $c^{\pop(\tau)}$ as:
$$
\ker \left(c^{\pop(\tau)} \right) = \bigcup_{x \in \mathcal{S}}\ker\left( c(x)^{\pop(\tau)}\right)
$$
Hence we may give a weaker version of the previous lemma 
\begin{cor}
Under the conditions of lemma \ref{lem:localization} there exists $b \in \stat$ such that:
$$
 x \left(\Omega \right) - b \in   \ker \left(c^{\pop(\tau)} \right)
$$
\end{cor}
\begin{pf}[of lemma \ref{lem:localization}]
First of all, let's check that $\widetilde{c}(x)$ is $\tau$-periodic.\\
For any $i \in \pop$, $x'_i$ is $\tau(i)$-periodic and non constant for $x_i$ is so. Let's $U_i=U \cap E_i$, $F_i$ has to be injective on $U_i$. Hence, as $x$ is a trajectory of the cellular system, $F_i \left( \widetilde{c}(x)_i \right)$ must be periodic and then $\widetilde{c}(x)_i$ is $\tau(i)$-periodic. Therefore, $c(x)$ is $\tau$-periodic.\\
 
Now, according to the partition $\pop(\tau)=\left\{\pop_k \right\}_{k \in K}$ defined by $\tau$ (see section \ref{sec:ModelOfPopulation}), let $k \in K$ and $i \in \pop_k$. For any $M \in \N$ we define the following set:
$$
\pop_k^M=\pop_k \cap \llbracket 0, M \rrbracket
$$
The set $\tau \left( \pop^M_k \right)$ is now a finite dependent set, so that we can consider its \emph{lcm} $\tau_k^M$. Now, for any $j \in \pop^M_k $, $x_j$ and $\widetilde{c}(x)_j$ are $\tau_j^M$-periodic, so that, for any integer $N$:
\begin{align*}
\widetilde{c}(x)_i(t)
 & =  \frac{1}{N+1} \sum_{l=0}^{N} \widetilde{c}(x)_i\left(t+l \tau_k^M \right) \\
 & =  \frac{1}{N+1} \sum_{l=0}^{N} \sum_{j \in \pop} c_{ij} \big(x_i\left(t+l \tau_k^M \right)\big).x_j\left(t+l \tau_k^M \right) \\
 & =  \frac{1}{N+1} \sum_{l=0}^{N} \sum_{j \in \pop} c_{ij} \big(x_i(t)\big).x_j\left(t+l \tau_k^M \right) \\
 & =  \frac{1}{N+1} \sum_{l=0}^{N} \sum_{j \in \pop_k^M} c_{ij} \big(x_i(t)\big).x_j(t) \\
 & \quad +  \frac{1}{N+1}\sum_{l=0}^{N} \sum_{j \in \pop_k-\pop_k^M} c_{ij} \big(x_i(t)\big).x_j\left(t+l \tau_k^M \right) \\
 & \quad +  \frac{1}{N+1}\sum_{l=0}^{N} \sum_{j \in \overline{\pop_k}} c_{ij} \big(x_i(t)\big).x_j\left(t+l \tau_k^M \right) \\
 & =  \sum_{j \in \pop_k^M} c_{ij} \big(x_i(t)\big).x_j(t) \\
 & \quad +  \frac{1}{N+1}\sum_{l=0}^{N} \sum_{j \in \pop_k-\pop_k^M} c_{ij} \big(x_i(t)\big).x_j\left(t+l \tau_k^M \right) \\
 & \quad +  \frac{1}{N+1}\sum_{l=0}^{N} \sum_{j \in \overline{\pop_k}} c_{ij} \big(x_i(t)\big).x_j\left(t+l \tau_k^M \right) \\
\end{align*}
As $x$ is uniformly bounded, using the 3rd property of a coupling map (def. \ref{def:coupler}), we may invert the summation orders in the previous equation and compute the limits when $M \to +\infty$ and $N \to + \infty$ in any order. Thus we have:
\begin{align*}
 \widetilde{c}(x)_i(t) & =  \sum_{j \in \pop_k^M} c_{ij} \big(x_i(t)\big).x_j(t) \\
 & +   \sum_{j \in \pop_k-\pop_k^M} c_{ij} \big(x_i(t)\big) \left[ \frac{1}{N+1}\sum_{l=0}^{N}x_j\left(t+l \tau_k^M \right) \right] \\
 &  +   \sum_{j \in \overline{\pop_k}} c_{ij} \big(x_i(t)\big) \left[ \frac{1}{N+1}\sum_{l=0}^{N}x_j\left(t+l \tau_k^M \right) \right] \\
\end{align*}
for the same reasons, it's easy to show that:
and
$$
\lim_{M \to + \infty } \lim_{N \to + \infty } \sum_{j \in \pop_k-\pop_k^M} c_{ij} \big(x_i(t)\big) \left[ \frac{1}{N+1}\sum_{l=0}^{N}x_j\left(t+l \tau_k^M \right) \right] = 0 
$$
Now, as for all $j \in \overline{I_k}$, $\tau_k^M$ and $\tau(j)$ are non commensurable, if we denote $tau_j'$ the generator of $x_j$ group of period, as $\tau(j)=n_j \tau_j'$ for a a certain integer $n_j$, $\tau_k^M$ and $\tau_j'$ as well are non commensurable. Therefore, the sequence $\displaystyle \left(\frac{t+l \tau_k^M}{\tau_j'}\right)_{l \in \N}$ is equidistributed mod $1$, and we may apply some classic ergodic theorem (see for instance \cite{Sina76,Buni00}) and write:
$$
\lim_{N \to + \infty } \frac{1}{N+1}\sum_{l=0}^{N}x_j\left(t+l \tau_k^M \right)  = \frac{1}{\tau_j'} \int_0^{\tau(j)}x_j(s)ds=\frac{n_j}{\tau(j)} \int_0^{\tau(j)}x_j(s)ds
$$
We can now define the state $b$ as:
$$
b=[b_j]_{j \in \pop}  , \; b_j=\frac{n_j}{\tau(j)} \int_0^{\tau(j)}x_j(s)ds %\textnormal{ if } i \in I_j
$$
so that:
$$
\lim_{N \to + \infty } \sum_{j \in \overline{\pop_k}} c_{ij} \big(x_i(t)\big) \left[ \frac{1}{N+1}\sum_{l=0}^{N}x_j\left(t+l \tau_k^M \right) \right]=\sum_{j \in \overline{\pop_k}} c_{ij} \big(x_i(t)\big)b_j
$$ 
hence, we have shown that:
$$
\widetilde{c}(x)_i(t) =\sum_{j \in \pop_k} c_{ij} \big(x_i(t)\big).x_j(t) + \sum_{j \in \overline{\pop_k}} c_{ij} \big(x_i(t)\big).b_j 
$$
But, obviously, from the beginning we had:
$$
\widetilde{c}(x)_i(t) =\sum_{j \in \pop_k} c_{ij} \big(x_i(t)\big).x_j(t) + \sum_{j \in \overline{\pop_k}} c_{ij} \big(x_i(t)\big).x_j(t) 
$$
So that:
$$
\sum_{j \in \overline{\pop_k}}  c_{ij} \big(x_i(t)\big).x_j(t) = \sum_{j \in \overline{\pop_k}} c_{ij} \big(x_i(t)\big).b_j
$$
The previous work can be done for any $i$ which belongs to $\pop_k$, thus we can summarize in the following way (see previously defined notations):
$$
c \big(x(t) \big)^{\pop_k}.\left( x(t)- b \right)^{\overline{\pop_k}} = 0
$$
Again, the previous conclusion still holds for any $k \in K$, hence we may conclude using our notations:
$$
c \big(x(t) \big)^{\pop(\tau)}.\left( x(t)- b \right) = 0
$$
\end{pf}

In the next section we give some examples of results based upon this lemma. We will mainly show how lemma \ref{lem:localization} may be applied to turn synchronization issues (and existence of quasi-periodic solutions to a differential system) into algebraic problems. One of the main argument is that one wants to avoid periodic trajectories for which one cell is inert (a constant map), as it may be discarded from the population (in the case of an infinite population, this may lead to recurrence reasoning). 
\section{Applications}
\label{sec:Applications}

\subsection{Weakly injective coupler}
\label{sec:WeaklyInjectiveCoupler}
In this example we just write down an elementary property of $\widetilde{c}$ which ensure that a periodic trajectory must have an inert cell.\\

\subsection{Localization results with bounded states and finite population}
\label{sec:boundedstatesFinitePopulation}

\begin{defn}
Let $\widetilde{c}$ be a cellular coupler on $\pop$. $\widetilde{c}$ is said to be weakly injective if for any non trivial partition $\pop(\tau)$ of $\pop$ there exist $i \in \pop$ such that:
$$
\forall x \in \stat, \; \ker \left(c(x)^{\pop(\tau)} \right) \cap E_i = \{0\}
$$ 
\end{defn}
Now we can state a simple result:   
\begin{prop}
\label{prop:weaklyinjective}
Under the conditions of lemma \ref{lem:localization}, if $\widetilde{c}$ is weakly injective and if $x$ is a $\tau$-periodic trajectory of the cellular system, then $\tau(\pop)$ is a dependent set.
\end{prop}
\begin{pf}
Let suppose that $\pop(\tau)$ is not trivial, applying lemma \ref{lem:localization} we know that:
$$
c \big(x(t) \big)^{\pop(\tau)}.\left( x(t)- b \right) = 0
$$
As $\widetilde{c}$ is weakly injective, there exists $i \in \pop$ such that:
$$
\forall t \in \Omega, \; x(t)_i=b_i
$$
which contradicts the definition of a periodic trajectory.
\end{pf}
This result may be restated in terms of quasi-periodic solution of the cellular system:
\begin{prop}
\label{prop:weaklyinjectivebis}
Under the conditions of lemma \ref{lem:localization}, if $\widetilde{c}$ is weakly injective and if $\tau$ is bounded, the cellular system has no quasi-periodic solution.
\end{prop}

The next example deals with some topological properties of a coupler (how it connects cells together).
\subsection{Chained cellular system}

We first study the case of differential systems for which the spaces $E_i$ have same finite dimension and are coupled with $k$-nearest neighbors (the finite dimension condition isn't necessary, but it makes the exposure simpler). 
This case is formally described by a cellular system $(\pop,F_{\pop},\widetilde{c})$ where $\pop$ is countable, all $\dim(E_i)=n$ and $\widetilde{c}$ satisfies:
$$
\forall i,j \in \pop, |j-i| > k \Rightarrow c_{ij}=0
$$
This is what we call a \emph{chained cellular system}. 
Adding the following condition on the coupler, we may reach a general result:
\begin{defn}
A cellular coupler $\widetilde{c}$ is said to have \emph{full rank} if for any $i,j \in \pop$ and $x \in \mathcal{S}$ the map $c_{ij}(x)$ has full rank
\end{defn}
\begin{prop}
\label{theorem:dependency1}
Let $(\pop,F_{\pop},\widetilde{c})$ be a chained cellular system coupled with $k$-nearest neighbors (all $E_i$ having same finite dimension). Let $F_{\pop}$ be injective on $U \subset \mathcal{S}$ and $x$ a $\tau$-periodic trajectory that stays in $U$. If $\widetilde{c}$ has maximal rank and if there exists $I \in \pop(\tau)$ which contains $2k$ consecutive cells, i.e. there exists $i \in \pop$ such that:
$$
\llbracket i, i+2k-1 \rrbracket \subset I 
$$ 
Then $\pop(\tau)=\{\pop\}$ (equivalently, $\tau(\pop)$ is a dependent set).
\end{prop}  
\begin{pf}
Let suppose that $I \ne \pop$. There must exist $\llbracket i, i+2k \rrbracket \subset I$, such that $i-1 \notin I$. Then, line $i+k-1$ of the matrix $c(x(t))^{\pop(\tau)}$ contains only one non zero element $c_{i+k-1,i-1}$. As this linear map is injective for any $t\in \Omega$, we know that:
$$
\ker(c(x(t))^{\pop(\tau)}) \bigcap E_{i-1} = \{0\}
$$
Applying lemma \ref{lem:localization} we know that there exists $b_{i-1} \in E_{i-1}$ such that for any $t \in \Omega$:
$$
x_{i-1}(t)- b_{i-1}  \in \ker \left(c(x(t))^{\pop(\tau)} \right) \bigcap E_{i-1}
$$
\emph{i.e.} $x_{i-1}(t)=b_{i-1}$ is a constant map, which contradicts the definition of a periodic trajectory. So we can conclude that $I= \pop$. 
\end{pf}
If we assume that $\tau$ is bounded, this result may be restated as: ``as soon as $k$ consecutive cells are synchronized (locked frequencies), then all the population is syncrohnized''.\\

Moreover, we may drop the assumptions made on the $E_i$ and reach an interesting connecting result concerning the case when $k=1$. 
\begin{prop}
\label{theorem:dependency2}
Let $(\pop,F_{\pop},\widetilde{c})$ be a chained cellular system coupled with $1$-nearest neighbor. Let $F_{\pop}$ be injective on $U \subset \mathcal{S}$ and $x$ a $\tau$-periodic trajectory that stays in $U$.  If $\widetilde{c}$ has maximal rank and if there exists two sets $I_1$ and $I_2$ in $\pop(\tau)$ such that:
$$
\llbracket i, i+1 \rrbracket \subset I_1 \quad \llbracket i+2, i+3 \rrbracket \subset I_2
$$ 
Then $I_1 = I_2$.
\end{prop}
\begin{pf}
Let suppose that the cells $i+1$ have non commensurable periods with those of the cells $i+2$ (\emph{i.e.} $I_1 \ne I_2$). Following the previous proof, we know that the lines $i+1$ and $i+2$ of the matrix $c(x(t))^{\pop(\tau)}$ contains only one non zero element, respectively $c_{i+1,i+2}$ and $c_{i+2,i+1}$. But, we recall that for any $t \in \Omega$: 
$$
c_{i+1,i+2}(x_{i+1}(t)) : E_{i+2} \to E_{i+1} \quad ,\quad c_{i+2,i+1}(x_{i+2}(t)) : E_{i+1} \to E_{i+2}
$$
As the coupler has maximal rank, one of the previous map must be injective for all $t \in \Omega$ (here we use the Axiom of choice to ensure that there must be a linear injection between the two Banach spaces). Using the same argument we may conclude that either $x_{i+1}$ is a constant map, either it's $x_{i+2}$, leading to a contradiction.
\end{pf}
Moreover, one could restate those results in terms of quasi periodic solutions of differential systems, but it may sound less intuitive. We will do it in the next sections.\\

For the next example, we add some regularity conditions on the cellular system which lead to a interesting description of $\stat$.  
\subsection{Localization results with bounded states}
\label{sec:boundedstates}
As $\left(\mathcal{S}_b,\|.\|_{\infty} \right)$ is a Banach space, the classic Picard�Lindel�f theorem is valid and we can give a version adapted to cellular systems (we stress the point that in this section, any vector field $F_{\pop}$ has to be a vector field on $\mathcal{S}_b$, as well for any cellular coupler $\widetilde{c}$, which brings some constraint on the families $(F_i)_{i \in \pop}$ and $(c_{ij})_{(i,j) \in \pop^2}$). 
\begin{prop}
\label{prop:existence}
If $F_{\pop}: \mathcal{S}_b \to \mathcal{S}_b$ and $\widetilde{c}$ are locally lipshcitz, which is the case if for any $x \in \mathcal{S}_b$ there exists a neighborhood $V=\displaystyle \prod_{i \in \pop} V_i$, a positive number $k$ and a sequence $(k_j)_{j \in \pop}$ of positive numbers such that:
\begin{enumerate}
	\item $ \displaystyle \forall y,z \in V, \forall i \in \pop, \; \|F_i(y_i)-F_i(z_i)\|_i \le k \|y_i-z_i\|_i$
	\item $ \displaystyle \forall y,z \in V, \forall i \in \pop, \; \|c_{ij}(y_i)-c_{ij}(z_i)\|_{(E_j,E_i)} \le k_j \|y_i-z_i\|_i$
	\item $ \displaystyle \sum_{j \in \pop} k_j < +\infty$ 
\end{enumerate}
then, given any initial condition $(t^0,x^0)$ in $\R \times \mathcal{S}_b$, the cellular coupling admits a unique maximal solution $x$ that satisfies $x(t^0)=x^0$.
\end{prop}

Before stating our localization result, we need to define the sets that any non synchronized periodic trajectory of the cellular system must avoid (or, with the classical point of view, any quasi periodic solution).

\begin{defn}
Let $\widetilde{c}$ be a cellular coupler on $\pop$. The set of regular points for $\widetilde{c}$ is defined as:
$$
R(\widetilde{c})=\left\{x \in \mathcal{S}, \; \forall \pop(\tau) \textnormal{ partition of }\pop, \; c(x)^{\pop(\tau)} \textnormal{ is injective }\right\}
$$ 
We say that $\widetilde{c}$ is regular if $R(\widetilde{c})=\mathcal{S}$.
\end{defn}
Now we can state a localization result:   
\begin{prop}
\label{prop:localization2}
Under the conditions of lemma \ref{lem:localization} and proposition \ref{prop:existence}, if there exists a infinite compact subset $V \subset \Omega$ such that:   
$$
\forall t \in V, \; x(t) \in R(\widetilde{c})
$$ 
then $\tau(\pop)$ is a dependent set.
\end{prop}

One can rewrite this result in terms of differential systems:
\begin{prop}
\label{prop:dependency2bis}
Under the conditions of lemma \ref{lem:localization} and proposition \ref{prop:existence}, and if $\tau$ is bounded, a quasi-periodic trajectory must ``avoid'' $R(\widetilde{c})$ (it can't cross this set on an infinite compact subset of $\Omega$).
\end{prop}
\begin{pf}[of proposition \ref{prop:localization2}]
Let suppose that $\pop(\tau)$ is not trivial, applying lemma \ref{lem:localization} we know that:
$$
c \big(x(t) \big)^{\pop(\tau)}.\left( x(t)- b \right) = 0
$$
the assumptions made on $\widetilde{c}$ ensure that:
$$
\forall t \in V, \; x(t)= b 
$$
As $V$ has an accumulation point, we may conclude that there exists $t_0 \in V$ such that:
$$
x'(t_0)=0
$$
Proposition \ref{prop:existence} may be applied, hence we know that $t \mapsto x(t)$ is a constant map, which contradicts the definition of a periodic trajectory.
\end{pf}

The next example gives a more precise result in the case where the maps $c_{ij}$ don't depend on the state of the system (homogeneous coupler).
\subsection{Exact frequencies locking with homogeneous cellular coupler}
\label{sec:ExactFrequenciesLockingInFinitePopulationAndHomogeneousLinearCoupler}

If $x \in \mathcal{T}_{\tau}$, for any $i \in \pop$ the map $x_i$ equals its Fourier's series. We write:
$$
e^k_{\tau(i)}(t)=\exp\left(\frac{2 i \pi k t }{\tau(i)}\right) 
$$
and we define :
$$
\widehat{x}_i(k) = \frac{1}{\tau(p)} \int_{0}^{\tau(p)} x_i(t) e_{\tau(i)}^k(t) dt
$$
so that we have :
$$
x= \sum_{k \in \Z}  \widehat{x}(k) e^k   
$$
\textit{i.e.} $\forall i \in \pop$:
$$
x_i(t)= \sum_{k \in \Z} \widehat{x}_i(k) e_{\tau(i)}^k(t)
$$
with normal convergence (note that $\widehat{x}_i(k)$ is $E_i$-valued).\\

\begin{thm}
\label{theorem:k/2}
Under the conditions of lemma \ref{lem:localization}, let $\widetilde{c}$ be homogeneous and regular. If $\tau$ is a bounded period on $\pop$ and $x$ a $\tau$-periodic trajectory of the cellular system then $\tau$ is constant on $\pop$.
\end{thm}
\begin{rem}
As this result is true as soon as $\tau$ is a period of $x$, it may be applied to the minimal periods of each $x_i$, then its conclusion is that all cells have exactly the same minimal period.
\end{rem}
\begin{pf}
As $\widetilde{c}$ is homogeneous, we may identify it with $c$.\\
Moreover, applying lemma \ref{lem:localization} we know that $\tau(\pop)$ is a dependent set (unless at least one of the $x_i$ would be a constant map). We now have to prove that $\tau$ is constant on $\pop$.\\ 

Let's write a partition of $\pop$ according to $\tau$'s values on $\pop$ (we must recall that $\tau$ is supposed bounded):
$$
\left\{\pop_1,\pop_2,\ldots,\pop_K\right\}
$$
such that 
$$
\forall 1 \le k \le K, \; \tau(\pop_k)=\tau_k 
$$
and $\tau_l \ne \tau_k$ if $l \ne k$.\\

We now suppose that $K>1$.\\

As $\tau(\pop)$ is a finite dependent set, there exists $n_1, \ldots , n_K$ integers and $\tau_0$ (the \emph{lcm}) 
such that:
$$
\tau_0 = n_1 \tau_1 = n_2 \tau_2 = \ldots = n_K \tau_K
$$

The trajectory $x$ is $\tau_0$-periodic. We may therefore write its Fourier's series:
$$
x(t) = \sum_{l \in \Z} \widehat{x}(l) e_{\tau_0}^l(t)
$$
and as well for $c.x$ :
$$
(c.x)(t)=\sum_{l \in \Z} \widehat{c.x}(l) e_{\tau_0}^l(t)
$$
uniqueness of Fourier coefficients forces them to satisfy:
$$
\widehat{c.x}(l) = c \, \widehat{x}(l)  
$$
So that, for any $i \in \pop$:
$$
\widehat{c.x}_i(l)  = \sum_{j=1}^k c_{ij} \widehat{x}_j(l) 
$$
Now, let $i \in \pop_k$, the properties of Fourier decomposition ensure that $\widehat{x}_i(l)$ and $c.\widehat{x}_i(l)$ are zero as soon as $ n_k$ does not divide $l$ (as $(c.x)_i$ and $x_i$ are $\tau_k$-periodic and $\tau_0 = n_k \tau_k$).\\
So, if $l \in \Z$, let's define $I(l)$ as:
$$
I(l) =\left\{ k \in \{1, \ldots ,K \}  , \;  n_k \not| \; l \right\}
$$
For any integer $l$, if $k \in I(l)$ and $i \in \pop_k$, then $\widehat{x}_i(l)=c.\widehat{x}_i(l)=0$, so that (with similar convergence arguments that in the proof of lemma \ref{lem:localization}):
\begin{align*}
c.\widehat{x}_i(l) & = \sum_{j=1}^k c_{ij} \widehat{x}_j(l) \\
 0 & = \sum_{j \in I(l) }^k c_{ij} \widehat{x}_j(l) + \sum_{j \notin I(l) }^k c_{ij} \widehat{x}_j(l) \\
 0 & =  \sum_{j \notin I(l) }^k c_{ij} \widehat{x}_j(l) \\
\end{align*}
This last property, (observable on figure \ref{fig:matrice_extraite_produit}), can be writen as:
$$
\forall l \in \Z \quad c^{I(l)} \, \widehat{x}(l)^{\overline{I(l)}} = 0
$$
\begin{figure}[ht]
	\centering
		\includegraphics[width=200pt]{matrice_extraite_produit.eps}
	\caption{Constraints on the Fourier's coefficients $\displaystyle \hat{x}(l)$.}
	\label{fig:matrice_extraite_produit}
\end{figure}

This property holds for any integer $l$, and is empty when $l$ is a multiple of all the $n_i$. So that, if $\pop(l)$ is the partition of $\pop$ defined as:
$$
\pop(l)=\left\{ I(l), \overline{I(l)}\right\}
$$
we can re-write it as:
$$
\forall l \in \Z \quad c^{\pop(l)} \, \widehat{x}(l) = 0
$$

Let's now consider $\pop_1\ne \pop_2$ (this is possible as $K>1$). As those two classes are distinct, there exists $l$ such that $n_1$ does not divide $l$ and $n_2$ divides $l$. As $c$ is regular, $c^{\pop(l)}$ is thereby injective. We deduce that:
$$
\widehat{x}(l)^{\overline{I(l)}} = 0
$$
This proves that for any $l$ divisible by $n_2$ and not by $n_1$, $\widehat{x}(l)^{\overline{I(l)}}$ is zero. Thus, for any coefficient of $\widehat{x}(l)^{\overline{I(l)}}$ to be non zero, $n_1$ must divide $l$, and consequently (as none of the $x_i$ is a constant map) for all $i \in \pop_2$, $x_i(t)$ is $n_1 \tau_0$ periodic. This is incompatible with the partition of $\pop$. Thus, $K=1$ and thereby $\tau$ is a constant map (in other words, $\pop$ is synchronized). 
\end{pf}

\subsection{Perspectives of application to classical differential systems}
\label{sec:ApplicationDuClassicalDifferentialSystems}
In this last section, we show how the cellular systems point of view may be applied to classic differential systems and how dealing with different Banach spaces $E_i$ may be useful. This discussion will be enlightened with a really simple example (finite population).\\

Let $E$ be a Banach space and $F$ a vector field on $E$. We want to see how this differential equation may be seen as a cellular system. For instance, one could consider a simple conservative system on $E=\R^4$ with an Hamilton's equation given by (see \cite{Luo06})
$$
\begin{array}{rcl}
x'_1&=& y_1 \\
y'_1&=& \alpha x_1 - \beta x_1^3 + \eps x_2 \\
x'_2&=& y_2 \\
y'_2&=& - \gamma x_2 + \eps x_1
\end{array}
$$  
The first step is to identify the different cells of $\pop$. The first step is to factorize each term in the equations according to the different variables. For example, the second equation may be seen as:
$$
y'_1=(\alpha - \beta x_1^2) x_1 + \eps x_2
$$
So that the term $(\alpha - \beta x_1^2)$ has to be a part of the coupler we are building. Moreover, as it's the equation giving $y'_1$, and as the way a cell computes how it interprets the population's state depends only on its own state, $x_1$ and $y_1$ have to belong to the same cell. In this simple example it's the only case where to variables have to be gathered in the same cell. In the end, this leads to the following structure of cellular system:
$$
\pop=\{1,2,3\}
$$
with the Banach spaces:
$$
E_1= \R^2, E_2=E_3=\R
$$
As it should often be the case, the associated vector fields are just identity maps on $E_i$, and the coupler is then:
$$
\left[ \begin{array}{ccc} c_{11} & c_{12} & c_{13} \\
													c_{21} & c_{22} & c_{23} \\
													c_{31} & c_{32} & c_{33} \\
													\end{array}\right]
$$
with
$$
\app{c_{11}}{E_1}{\mathcal{L}(E_1)}{(x_1,y_1)}{\left[ \begin{array}{cc} 0 & \; 1  \\
													\alpha - \beta x_1^2 & \; 0\\
													\end{array}\right]}
\quad 
\app{c_{12}}{E_2}{\mathcal{L}(E_2,E_1)}{x_2}{\left[ \begin{array}{c} 0   \\
													 \eps \\
													\end{array}\right]}
\quad
\app{c_{13}}{E_3}{\mathcal{L}(E_3,E_1)}{y_2}{\left[ \begin{array}{c} 0   \\
													 0 \\
													\end{array}\right]}
$$
$$
\app{c_{21}}{E_1}{\mathcal{L}(E_1,E_2)}{(x_1,y_1)}{\left[ \begin{array}{cc} 0 \; & \; 0   \\
													\end{array}\right]}
\quad
\app{c_{22}}{E_2}{\mathcal{L}(E_2)}{x_2}{\left[ \begin{array}{c} 0   \\
													\end{array}\right]}
\quad
\app{c_{23}}{E_3}{\mathcal{L}(E_3,E_2)}{y_2}{\left[ \begin{array}{c} 1   \\
													\end{array}\right]}
$$
$$
\app{c_{31}}{E_1}{\mathcal{L}(E_1,E_3)}{(x_1,y_1)}{\left[ \begin{array}{cc} \eps \; &  \; 0   \\
													\end{array}\right]}
\quad
\app{c_{32}}{E_2}{\mathcal{L}(E_2,E_3)}{x_2}{\left[ \begin{array}{c} - \gamma   \\
													\end{array}\right]}
\quad
\app{c_{33}}{E_3}{\mathcal{L}(E_3)}{y_2}{\left[ \begin{array}{c} 0   \\
													\end{array}\right]}
$$
Now, before applying some of the previous techniques, we may compute the different decomposition of $c$ upon different non trivial partitions of $\pop$. Those partitions are:
$$
P_1=\big\{\{ 1\},\{ 2\},\{ 3\big\} \} ,\;  P_2=\big\{\{ 1,2\},\{ 3\} \big\} ,\; P_3=\big\{\{ 1,3\},\{ 2\} \big\} ,\; P_4=\big\{\{ 1\},\{ 2,3\} \big\}
$$
which give:
$$
c^{P_1}=\left[ \begin{array}{ccc} 0 & c_{12} & c_{13} \\
													c_{21} & 0 & c_{23} \\
													c_{31} & c_{32} & 0 \\
													\end{array}\right]
\quad 
c^{P_2}=\left[ \begin{array}{ccc} 0 & 0 & c_{13} \\
													0 & 0 & c_{23} \\
													c_{31} & c_{32} & 0 \\
													\end{array}\right]
\quad 
c^{P_3}=\left[ \begin{array}{ccc} 0 & c_{12} & 0 \\
													c_{21} & 0 & c_{23} \\
													0 & c_{32} & 0 \\
													\end{array}\right]
\quad
c^{P_4}=\left[ \begin{array}{ccc} 0 & c_{12} & c_{13} \\
													c_{21} & 0 & 0 \\
													c_{31} & 0 & 0 \\
													\end{array}\right]
$$
Now, in order to simplify, we replace the $c_{ij}$ that are identically zero by $0$, we obtain the following different matrices:
$$
c^{P_1}=\left[ \begin{array}{ccc} 0 & c_{12} & 0 \\
													0 & 0 & c_{23} \\
													c_{31} & c_{32} & 0 \\
													\end{array}\right]
\quad 
c^{P_2}=\left[ \begin{array}{ccc} 0 & 0 & 0 \\
													0 & 0 & c_{23} \\
													c_{31} & c_{32} & 0 \\
													\end{array}\right]
\quad 
c^{P_3}=\left[ \begin{array}{ccc} 0 & c_{12} & 0 \\
													0 & 0 & c_{23} \\
													0 & c_{32} & 0 \\
													\end{array}\right]
\quad
c^{P_4}=\left[ \begin{array}{ccc} 0 & c_{12} & 0 \\
													0 & 0 & 0 \\
													c_{31} & 0 & 0 \\
													\end{array}\right]
$$
In the end, writing the coupler as an application from $\mathcal{S}$ to $\mathcal{L}(\mathcal{S})$, one finds those four matrices:
$$
\left[ \begin{array}{cccc} 0 & 0 & 0 & 0\\
													0 &0& \eps & 0 \\
													0&0&0&1 \\
													\eps&0& -\gamma &0 \\
													\end{array}\right]
\quad 
\left[ \begin{array}{cccc} 0 & 0 & 0 & 0\\
													0 &0& 0 & 0 \\
													0&0&0&1 \\
													\eps & 0& -\gamma &0 \\
													\end{array}\right]
\quad 
\left[ \begin{array}{cccc} 0 & 0 & 0 & 0\\
													0 &0& \eps & 0 \\
													0&0&0&1 \\
													0&0& -\gamma &0 \\
													\end{array}\right]
\quad
\left[ \begin{array}{cccc} 0 & 0 & 0 & 0\\
													0 &0& \eps & 0 \\
													0&0&0&0 \\
													\eps&0& 0 &0 \\
													\end{array}\right]$$
													
At this point, we just have to check that the coupler is weakly injective:
$$
\ker\left( c^{P_1}\right) \cap E_2 = \ker\left( c^{P_4}\right) \cap E_2 = \{0\}
$$
$$
 \ker\left( c^{P_2}\right) \cap E_3 = \ker\left( c^{P_3}\right) \cap E_3 = \{0\} 
$$
So, we may apply the proposition \ref{prop:weaklyinjective} and without any analytic calculus, state that this differential system may not admit any quasi-periodic solution. In other words, in case there exists periodic trajectories (which is well known to be true) they must be synchronized.\\

Moreover, those conclusions may hold in a more general case were the $c_{ij}$ are less simple, and we can easily produce a result without any effort:
\begin{prop}[Generalized coupled pendulum]
Let's consider a differential system which is driven by the following equations:
$$
\begin{array}{rcl}
x'_1&=& a_{1}(x_1,y_1)x_1 + a_{2}(x_1,y_1)y_1+a_{3}(x_1,y_1)x_2+a_{4}(x_1,y_1)y_2 \\
y'_1&=& a_{5}(x_1,y_1)x_1 + a_{6}(x_1,y_1)y_1 \\
x'_2&=& a_7(x_2)x_2+u(x_2) y_2 \\
y'_2&=&  \eps(y_2) x_1 + a_8(y_2)y_1 - \gamma(y_2) x_2 +a_9(y_2)y_2
\end{array}
$$  
If the maps $u$ and $\eps$ never vanish, then the systems has no continuous quasi-periodic solution.
\end{prop}
Indeed, in this really simple example, the continuity of a solution is the only property needed to ensure the validity of lemma \ref{lem:localization}. Moreover, this result doesn't have to be deep in itself, neither has it to be the most general one we could have deduced from the previous discussion. It's just a sketch of how one can handle some structure properties of a differential system, applying \ref{lem:localization}, without going into deep and specific calculus.
   
\section{Conclusion and discution}
\label{sec:Conclusion}

In this work we have built a general framework of cellular systems in order to handle a wide variety of coupled systems, and therefore a wide class of complex systems. We focused on an emergent property of those dynamical systems: the frequencies locking phenomenon. Usually one observes solutions of particular coupled systems and shows that within suitable conditions synchronization must occur. Those results are qualitatively dependent on the systems of interest and do not stand in the general cases. We tried to change our point of view and to bring out completing results. As we choose not to address the problem of the convergence to a periodic solution, we don't prove that synchronization ultimately happens. Instead, we consider the problem at its end: if one supposes that some coupled systems ``truly'' oscillates, then they must be synchronized, regardless to the individual dynamical systems as soon as the maps which define each of them are injective nearby the trajectories. The loss in time evolution information is compensated by very general results, almost independent from the individual differential systems to be coupled. For example, we proved results concerning countable coupled systems, each of them needing not to be finite dimensional. In most papers (see for instance \cite{Kope00}) this population of coupled systems is implicitly defined and has only two cells (sometimes a finite number $N$, and more rarely an infinity). Moreover, on the contrary of what most studies about synchronization issues state, we do not assume anything concerning the cells dynamics (especially, we don't assume that they are oscillators). We only assume that they exhibit periodic behaviors under the coupling effects (the first assumption implies the second, but the opposite is clearly false). 

We believe that this way of reaching general results about cellular systems gives some explanations about why the frequencies locking phenomenon emerges naturally in a large variety of coupled dynamical systems. Our results show that the following alternative is natural in many cases: either the whole population is synchronized, or its cells can't all have periodic behaviors. 

Another interesting perspective is to apply this strategy to differential systems, as we outlined in the end of the fourth section. We think that it could prove useful to understand the stability or instability of limit torus when one adds some perturbations to a differential system. For example, if one already knows that there's no periodic solution, even with perturbations, and if on is able to prove, using our strategy, that a quasi-periodic solutions disappears under the effects of those perturbation, some conclusions about strange attractor emergence may be reachable.  

Moreover, we have achieved some similar work on a natural generalization of this strategy to non countable population (in order to model natural systems, it's often necessary to handle continuous populations). We truly think that all those results are only a part of what can be done using cellular systems and that this work enlarges the possibilities of studying synchronization issues. But the scope of those kind of cellular systems may be beyond synchronization questions, as it's quite general and allows some theoretical studies. It could be a promising theoretical tool to model complex systems by the way of coupled differential systems.\\



% Bibliographic references with the natbib package:
% Parenthetical: \citep{Bai92} produces (Bailyn 1992).
% Textual: \citet{Bai95} produces Bailyn et al. (1995).
% An affix and part of a reference:
%   \citep[e.g.][Ch. 2]{Bar76}
%   produces (e.g. Barnes et al. 1976, Ch. 2).

%\begin{thebibliography}{}

\bibliographystyle{elsart-num-sort}
%\bibliographystyle{alpha}
 %\bibliography{bibliographiejabref}
 \begin{thebibliography}{10}
\expandafter\ifx\csname url\endcsname\relax
  \def\url#1{\texttt{#1}}\fi
\expandafter\ifx\csname urlprefix\endcsname\relax\def\urlprefix{URL }\fi

\bibitem{Luo06}
A.~C.~J.~Luo, Predictions of quasi-periodic and chaotic motions in nonlinear Hamiltonian systems, Chaos, Solitons \& Fractals 28~(3) (2006) 627--649.

\bibitem{Benn02}
M.~Bennett, M.~F. Schatz, H.~Rockwood, K.~Wiesenfeld, Huygens's clocks, Proc.
  R. Soc. Lond. 458 (2002) 563--579.

\bibitem{Buck88}
J.~Buck, Synchronous rhythmic flashing of fireflies. ii, Quarterly Review of
  Biology 63~(3) (1988) 265--289.

\bibitem{Buni00}
L.~Bunimovich, Dynamical Systems, Ergodic Theory and Applications,
  Encyclopaedia of Mathematical Sciences, 2nd ed., Springer, 2000.

\bibitem{Erme86}
G.~B. Ermentrout, W.~C. Troy, Phaselocking in a reaction-diffusion system with
  a linear frequency gradient., SIAM-J.-Appl.-Math. 46~(3) (1986) 359--367.

\bibitem{Gaub07}
L.~Gaubert, Auto-organisation et \'emergence dans les syst\`emes coupl\'es,
  individuation de donn\'ees issues de syst\`emes biologiques coupl\'es., Ph.D.
  thesis, Universit\'e de Bretagne Occidentale (2007).

\bibitem{Goed03}
J.-P. Goedgebuer, P.~Levy, L.~Larger, Laser cryptography by optical chaos, in:
  J.-P. Goedgebuer, N.~Rozanov, S.~T. andA.S. Akhmanov, V.~Panchenko (eds.),
  Optical Information, Data Processing and Storage, and Laser Communication
  Technologies, vol. 5135, 2003.

\bibitem{Goel02}
P.~Goel, B.~Ermentrout, Synchrony, stability, and firing patterns in
  pulse-coupled oscillators, PHYSICA D 163~(3-4) (2002) 191--216.

\bibitem{Gonz04}
D.~Gonze, J.~Halloy, A.~Goldbeter, Stochastic models for cicardian
  oscillations: Emergence of a biological rythm, International Journal of
  Quantum Chemistry 98 (2004) 228--238.

\bibitem{Hopp76}
F.~C. Hoppensteadt, J.~B. Keller, Synchronization of periodical cicada
  emergences, Science 194~(4262) (1976) 335--337.

\bibitem{Huge73}
C.~Huygens, Christiani Hugenii Zulichemii, Const. F. Horologium oscillatorium
  sive De motu pendulorum ad horologia aptato demonstrationes geometricae,
  Parisiis : Apud F. Muguet, Paris, France, 1673.

\bibitem{Izhi03}
E.~M. Izhikevich, F.~C. Hoppensteadt, Slowly coupled oscillators: Phase
  dynamics and synchronization, SIAM J. APPL. MATH. 63~(6) (2003) 1935--1953.

\bibitem{Kope00}
N.~Kopell, G.~B. Ermentrout, Handbook of Dynamical Systems II: Toward
  Applications, chap. Mechanisms of phase-locking and frequency control in
  pairs of coupled neural oscillators, Elsevier, Amsterdam, 2002, pp. 3--54.

\bibitem{Kura75}
Y.~Kuramoto, Self-entrainment of a population of coupled non-linear
  oscillators, in: International Symposium on Mathematical Problems in
  Theoretical Physics, vol.~39, Springer Berlin / Heidelberg, 1975, pp.
  420--422.

\bibitem{Kura84}
Y.~Kuramoto, Chemical oscillations, waves, and turbulence, Springer-Verlag, New
  York, 1984.

\bibitem{Manr04}
S.~C. Manrubia, A.~S. Mikhailov, D.~H. Zanette, Emergence of Dynamical Order.
  Synchronization Phenomena in Complex Systems, World Scientific, Singapore,
  2004.

\bibitem{Mich87}
D.~C. Michaels, E.~P. Matyas, J.~Jalife, Mechanisms of sinoatrial pacemaker
  synchronization: A new hypothesis, Circ. Res. 61 (1987) 704--714.

\bibitem{Peco97}
L.~M. Pecora, T.~L. Carroll, G.~A. Johnson, D.~J. Mar, J.~F. Heagy,
  Fundamentals of synchronization in chaotic systems, concepts, and
  applications, Chaos 7~(4) (1997) 520--543.

\bibitem{Piko01}
A.~Pikovsky, M.~Rosenblumand, J.~Kurths, Synchronization: A Universal Concept
  in Nonlinear Sciences, Cambridge University Press, 2001.

\bibitem{Prig67}
I.~Prigogine, Introduction to Thermodynamics of Irreversible Processes, Wiley,
  New York, 1967.

\bibitem{Prig69}
I.~Prigogine, Theoretical Physics and Biology, M Marois, North-Holland:
  Amesterdam, 1969.

\bibitem{Ren00}
L.~Ren, B.~G. Ermentrout, Phase locking in chains of multiple-coupled
  oscillators, PHYSICA D 143~(1-4) (2000) 56--73.

\bibitem{Sina76}
Y.~G. Sinai, Introduction to ergodic theory, Princeton University Press, 1976.

\bibitem{Stew02}
I.~Stewart, M.~Golubitsky, Patterns of oscillation in coupled cell systems, in:
  P.~P.Holmes, A.Weinstein (eds.), Geometry, Dynamics and Mechanics: 60th
  Birthday Volume for J.E. Marsden, Springer-Verlag, New York, 2002, pp.
  243--286.

\bibitem{Stro03}
S.~Strogatz, Sync: The Emergeing Science of Spontaneous Order, New York:
  Hyperion, 2003.

\bibitem{Winf67}
A.~Winfree, Biological rhythms and the behavior of populations of coupled
  oscillators, Journal of Theoretical Biology 16~(1) (1967) 15--42.

\bibitem{Winf90}
A.~Winfree, The geometry of biological time, Springer- Verlag, New York, 1990.

\bibitem{Wu94}
C.~Wu, L.~O. Chua, A unified framework for synchronization and control of
  dynamical systems, Tech. Rep. UCB/ERL M94/28, EECS Department, University of
  California, Berkeley (1994).

\end{thebibliography}


% \bibitem[Names(Year)]{label} or \bibitem[Names(Year)Long names]{label}.
% (\harvarditem{Name}{Year}{label} is also supported.)
% Text of bibliographic item

%\bibitem[]{}

%\end{thebibliography}



\end{document}













 



 

