X-Git-Url: https://bilbo.iut-bm.univ-fcomte.fr/and/gitweb/LiCO.git/blobdiff_plain/e0ade9eff447e924d284a39ac11e50707aa6529e..refs/heads/master:/PeCO-EO/articleeo.tex?ds=inline diff --git a/PeCO-EO/articleeo.tex b/PeCO-EO/articleeo.tex index f9a5f6d..3516881 100644 --- a/PeCO-EO/articleeo.tex +++ b/PeCO-EO/articleeo.tex @@ -24,9 +24,10 @@ The most important problem in a Wireless Sensor Network (WSN) is to optimize the use of its limited energy provision, so that it can fulfill its monitoring task as long as possible. Among known available approaches that can be used to improve power management, lifetime coverage optimization provides activity -scheduling which ensures sensing coverage while minimizing the energy cost. An approach called Perimeter-based Coverage Optimization protocol -(PeCO) is proposed. It is a hybrid of centralized and distributed methods: the region of -interest is first subdivided into subregions and the protocol is then +scheduling which ensures sensing coverage while minimizing the energy cost. In +this paper an approach called Perimeter-based Coverage Optimization protocol +(PeCO) is proposed. It is a hybrid of centralized and distributed methods: the +region of interest is first subdivided into subregions and the protocol is then distributed among sensor nodes in each subregion. The novelty of the approach lies essentially in the formulation of a new mathematical optimization model based on the perimeter coverage level to schedule sensors' activities. @@ -83,16 +84,17 @@ This paper makes the following contributions : architecture. \item A new mathematical optimization model is proposed. Instead of trying to cover a set of specified points/targets as in most of the methods proposed in - the literature, a mixed-integer program based on the perimeter - coverage of each sensor is formulated. The model involves integer variables to capture the - deviations between the actual level of coverage and the required level. + the literature, a mixed-integer program based on the perimeter coverage of + each sensor is formulated. The model involves integer variables to capture + the deviations between the actual level of coverage and the required level. Hence, an optimal schedule will be obtained by minimizing a weighted sum of these deviations. \item Extensive simulation experiments are conducted using the discrete event - simulator OMNeT++, to demonstrate the efficiency of the PeCO protocol. The PeCO protocol has been compared to two approaches found in the literature: - DESK~\citep{ChinhVu} and GAF~\citep{xu2001geography}, and also to the - protocol DiLCO published in~\citep{Idrees2}. DiLCO uses the same framework as - PeCO but is based on another optimization model for sensor scheduling. + simulator OMNeT++, to demonstrate the efficiency of the PeCO protocol. The + PeCO protocol has been compared to two approaches found in the literature: + DESK~\citep{ChinhVu} and GAF~\citep{xu2001geography}, and also to the protocol + DiLCO published in~\citep{Idrees2}. DiLCO uses the same framework as PeCO but + is based on another optimization model for sensor scheduling. \end{enumerate} The rest of the paper is organized as follows. In the next section some related @@ -205,8 +207,8 @@ decomposed into 4 phases: information exchange, leader election, decision, and sensing. The simulations show that DiLCO is able to increase the WSN lifetime and provides improved coverage performance. {\it In the PeCO protocol, a new mathematical optimization model is proposed. Instead of trying to cover a set - of specified points/targets as in the DiLCO protocol, an integer - program based on the perimeter coverage of each sensor is formulated. The model involves + of specified points/targets as in the DiLCO protocol, an integer program based + on the perimeter coverage of each sensor is formulated. The model involves integer variables to capture the deviations between the actual level of coverage and the required level. The idea is that an optimal scheduling will be obtained by minimizing a weighted sum of these deviations.} @@ -214,7 +216,6 @@ and provides improved coverage performance. {\it In the PeCO protocol, a new \section{ The P{\scshape e}CO Protocol Description} \label{sec:The PeCO Protocol Description} - \subsection{Assumptions and Models} \label{CI} @@ -243,11 +244,11 @@ $k$~sensors) if and only if each sensor in the network is $k$-perimeter-covered (perimeter covered by at least $k$ sensors). Figure~\ref{figure1}(a) shows the coverage of sensor node~$0$. On this figure, -sensor~$0$ has nine neighbors. For each neighbor the two points -resulting from the intersection of the two sensing areas have been reported on its perimeter (the -perimeter of the disk covered by the sensor~$0$). These points are -denoted for neighbor~$i$ by $iL$ and $iR$, respectively for left and right from -a neighboring point of view. The resulting couples of intersection points +sensor~$0$ has nine neighbors. For each neighbor the two points resulting from +the intersection of the two sensing areas have been reported on its perimeter +(the perimeter of the disk covered by the sensor~$0$). These points are denoted +for neighbor~$i$ by $iL$ and $iR$, respectively for left and right from a +neighboring point of view. The resulting couples of intersection points subdivide the perimeter of sensor~$0$ into portions called arcs. \begin{figure}[ht!] @@ -288,7 +289,7 @@ from the first intersection point after point~zero, and the maximum level coverage is determined for each interval defined by two successive points. The maximum level of coverage is equal to the number of overlapping arcs. For example, between~$5L$ and~$6L$ the maximum level of coverage is equal to $3$ -(the value is highlighted in yellow at the bottom of Figure~\ref{figure2}), which +(the value is given at the bottom of Figure~\ref{figure2}), which means that at most 2~neighbors can cover the perimeter in addition to node $0$. Table~\ref{my-label} summarizes for each coverage interval the maximum level of coverage and the sensor nodes covering the perimeter. The example discussed @@ -409,7 +410,6 @@ The pseudocode implementing the protocol on a node is given below. More precisely, Algorithm~\ref{alg:PeCO} gives a brief description of the protocol applied by a sensor node $s_k$ where $k$ is the node index in the WSN. - \begin{algorithm2e} \label{alg:PeCO} \caption{PeCO pseudocode} @@ -454,17 +454,21 @@ in the current period. Each sensor node determines its position and its subregion using an embedded GPS or a location discovery algorithm. After that, all the sensors collect position coordinates, remaining energy, sensor node ID, and the number of their one-hop live neighbors during the information exchange. -Both INFO packet and ActiveSleep packet contain two parts: header and data payload. The sensor ID is included in the header, where the header size is 8 bits. The data part includes position coordinates (64 bits), remaining energy (32 bits), and the number of one-hop live neighbors (8 bits). Therefore the size of the INFO packet is 112 bits. The ActiveSleep packet is 16 bits size, 8 bits for the header and 8 bits for data part that includes only sensor status (0 or 1). -The sensors inside a same region cooperate to elect a leader. The selection -criteria for the leader are (in order of priority): +Both INFO packet and ActiveSleep packet contain two parts: header and data +payload. The sensor ID is included in the header, where the header size is 8 +bits. The data part includes position coordinates (64 bits), remaining energy +(32 bits), and the number of one-hop live neighbors (8 bits). Therefore the size +of the INFO packet is 112 bits. The ActiveSleep packet is 16 bits size, 8 bits +for the header and 8 bits for data part that includes only sensor status (0 or +1). The sensors inside a same region cooperate to elect a leader. The +selection criteria for the leader are (in order of priority): \begin{enumerate} \item larger number of neighbors; \item larger remaining energy; \item and then, in case of equality, larger indexes. \end{enumerate} Once chosen, the leader collects information to formulate and solve the integer -program which allows to build the set of active sensors in the sensing -stage. +program which allows to build the set of active sensors in the sensing stage. \section{Perimeter-based Coverage Problem Formulation} \label{cp} @@ -604,12 +608,12 @@ coverage task. This value corresponds to the energy needed by the sensing phase, obtained by multiplying the energy consumed in the active state (9.72 mW) with the time in seconds for one period (3600 seconds), and adding the energy for the pre-sensing phases. According to the interval of initial energy, a sensor may -be active during at most 20 periods. the information exchange to update the coverage -is executed every hour, but the length of the sensing period could be reduced -and adapted dynamically. On the one hand a small sensing period would allow the network to -be more reliable but would have resulted in higher communication costs. On the -other hand the choice of a long duration may cause problems in case of nodes -failure during the sensing period. +be active during at most 20 periods. the information exchange to update the +coverage is executed every hour, but the length of the sensing period could be +reduced and adapted dynamically. On the one hand a small sensing period would +allow the network to be more reliable but would have resulted in higher +communication costs. On the other hand the choice of a long duration may cause +problems in case of nodes failure during the sensing period. The values of $\alpha^j_i$ and $\beta^j_i$ have been chosen to ensure a good network coverage and a longer WSN lifetime. Higher priority is given to the @@ -642,19 +646,20 @@ approach. subregions during the current sensing phase and $N$ is total number of grid points in the sensing field. A layout of $N~=~51~\times~26~=~1326$~grid points is considered in the simulations. -\item {\bf Active Sensors Ratio (ASR)}: a major objective of the proposed protocol is to - activate as few nodes as possible, in order to minimize the communication - overhead and maximize the WSN lifetime. The active sensors ratio is defined as - follows: +\item {\bf Active Sensors Ratio (ASR)}: a major objective of the proposed + protocol is to activate as few nodes as possible, in order to minimize the + communication overhead and maximize the WSN lifetime. The active sensors ratio + is defined as follows: \begin{equation*} \scriptsize - \mbox{ASR}(\%) = \frac{\sum\limits_{r=1}^R \mbox{$|A_r^p|$}}{\mbox{$|J|$}} \times 100 + \mbox{ASR}(\%) = \frac{\sum\limits_{r=1}^R \mbox{$|A_r^p|$}}{\mbox{$|S|$}} \times 100 \end{equation*} where $|A_r^p|$ is the number of active sensors in the subregion $r$ in the sensing period~$p$, $R$ is the number of subregions, and $|J|$ is the number of sensors in the network. -\item {\bf Energy Saving Ratio (ESR)}:this metric, which shows the ability of a protocol to save energy, is defined by: +\item {\bf Energy Saving Ratio (ESR)}: this metric, which shows the ability of a + protocol to save energy, is defined by: \begin{equation*} \scriptsize \mbox{ESR}(\%) = \frac{\mbox{Number of alive sensors during this round}} @@ -684,16 +689,16 @@ approach. \subsection{Simulation Results} - -The PeCO protocol has been implemented in OMNeT++~\citep{varga} simulator in order to assess and analyze its performance. -The simulations were run on a DELL laptop with an Intel Core~i3~2370~M (1.8~GHz) -processor (2 cores) whose MIPS (Million Instructions Per Second) rate is equal -to 35330. To be consistent with the use of a sensor node based on Atmels AVR -ATmega103L microcontroller (6~MHz) having a MIPS rate equal to 6, the original -execution time on the laptop is multiplied by 2944.2 $\left(\frac{35330}{2} -\times \frac{1}{6} \right)$. Energy consumption is calculated according to the -power consumption values, in milliWatt per second, given in Table~\ref{tab:EC}, -based on the energy model proposed in \citep{ChinhVu}. +The PeCO protocol has been implemented in OMNeT++~\citep{varga} simulator in +order to assess and analyze its performance. The simulations were run on a DELL +laptop with an Intel Core~i3~2370~M (1.8~GHz) processor (2 cores) whose MIPS +(Million Instructions Per Second) rate is equal to 35330. To be consistent with +the use of a sensor node based on Atmels AVR ATmega103L microcontroller (6~MHz) +having a MIPS rate equal to 6, the original execution time on the laptop is +multiplied by 2944.2 $\left(\frac{35330}{2} \times \frac{1}{6} \right)$. Energy +consumption is calculated according to the power consumption values, in +milliWatt per second, given in Table~\ref{tab:EC}, based on the energy model +proposed in \citep{ChinhVu}. \begin{table}[h] \centering @@ -715,16 +720,16 @@ based on the energy model proposed in \citep{ChinhVu}. The modeling language for Mathematical Programming (AMPL)~\citep{AMPL} is used to generate the integer program instance in a standard format, which is then -read and solved by the optimization solver GLPK (GNU linear Programming Kit +read and solved by the optimization solver GLPK (GNU Linear Programming Kit available in the public domain) \citep{glpk} through a Branch-and-Bound method. In practice, executing GLPK on a sensor node is obviously intractable due to the -huge memory use. Fortunately, to solve the optimization problem, the use of +huge memory use. Fortunately, to solve the optimization problem, the use of commercial solvers like CPLEX \citep{iamigo:cplex} which are less memory -consuming and more efficient is possible, or a lightweight heuristic may be implemented. For example, -for a WSN of 200 sensor nodes, a leader node has to deal with constraints -induced by about 12 sensor nodes. In that case, to solve the optimization -problem a memory consumption of more than 1~MB can be observed with GLPK, -whereas less than 300~KB would be needed with CPLEX. +consuming and more efficient is possible, or a lightweight heuristic may be +implemented. For example, for a WSN of 200 sensor nodes, a leader node has to +deal with constraints induced by about 12 sensor nodes. In that case, to solve +the optimization problem a memory consumption of more than 1~MB can be observed +with GLPK, whereas less than 300~KB would be needed with CPLEX. Besides PeCO, three other protocols will be evaluated for comparison purposes. The first one, called DESK, is a fully distributed coverage algorithm @@ -770,13 +775,13 @@ allows later a substantial increase of the coverage performance. \subsubsection{Active Sensors Ratio} -Minimizing the number of active sensor nodes in each period is essential to minimize the -energy consumption and thus to maximize the network lifetime. +Minimizing the number of active sensor nodes in each period is essential to +minimize the energy consumption and thus to maximize the network lifetime. Figure~\ref{figure6} shows the average active nodes ratio for 200 deployed -nodes. DESK and GAF have 30.36~\% and 34.96~\% active nodes for -the first fourteen rounds, and the DiLCO and PeCO protocols compete perfectly with -only 17.92~\% and 20.16~\% active nodes during the same time interval. As the -number of periods increases, the PeCO protocol has a lower number of active nodes in +nodes. DESK and GAF have 30.36~\% and 34.96~\% active nodes for the first +fourteen rounds, and the DiLCO and PeCO protocols compete perfectly with only +17.92~\% and 20.16~\% active nodes during the same time interval. As the number +of periods increases, the PeCO protocol has a lower number of active nodes in comparison with the three other approaches and exhibits a slow decrease, while keeping a greater coverage ratio as shown in Figure \ref{figure5}. @@ -789,18 +794,17 @@ keeping a greater coverage ratio as shown in Figure \ref{figure5}. \subsubsection{Energy Saving Ratio} - -The simulation results show that the protocol PeCO saves - efficiently energy by turning off some sensors during the sensing phase. As - shown in Figure~\ref{figure7}, GAF provides better energy saving than PeCO for - the first fifty rounds. Indeed GAF balances the energy consumption among - sensor nodes inside each small fixed grid and thus permits to extend the life - of sensors in each grid fairly. However, at the same time it turns on a large - number of sensors and that leads later to quickly deplete sensor's batteries. - DESK algorithm shows less energy saving compared with other approaches. In - comparison with PeCO, DiLCO protocol usually provides lower energy saving - ratios. Moreover, it can be noticed that after round fifty, PeCO protocol - exhibits the slowest decrease among all the considered protocols. +The simulation results show that the protocol PeCO saves efficiently energy by +turning off some sensors during the sensing phase. As shown in +Figure~\ref{figure7}, GAF provides better energy saving than PeCO for the first +fifty rounds. Indeed GAF balances the energy consumption among sensor nodes +inside each small fixed grid and thus permits to extend the life of sensors in +each grid fairly. However, at the same time it turns on a large number of +sensors and that leads later to quickly deplete sensor's batteries. DESK +algorithm shows less energy saving compared with other approaches. In +comparison with PeCO, DiLCO protocol usually provides lower energy saving +ratios. Moreover, it can be noticed that after round fifty, PeCO protocol +exhibits the slowest decrease among all the considered protocols. \begin{figure}[h!] %\centering @@ -817,15 +821,15 @@ The effect of the energy consumed by the WSN during the communication, computation, listening, active, and sleep status is studied for different network densities and the four approaches compared. Figures~\ref{figure8}(a) and (b) illustrate the energy consumption for different network sizes and for -$Lifetime_{95}$ and $Lifetime_{50}$. The results show that the PeCO protocol is the most -competitive from the energy consumption point of view. As shown by both figures, -PeCO consumes much less energy than the other methods. One might think that the -resolution of the integer program is too costly in energy, but the results show -that it is very beneficial to lose a bit of time in the selection of sensors to -activate. Indeed the optimization program allows to reduce significantly the -number of active sensors and also the energy consumption while keeping a good -coverage level. The energy overhead when increasing network -size is the lowest with PeCO. +$Lifetime_{95}$ and $Lifetime_{50}$. The results show that the PeCO protocol is +the most competitive from the energy consumption point of view. As shown by both +figures, PeCO consumes much less energy than the other methods. One might think +that the resolution of the integer program is too costly in energy, but the +results show that it is very beneficial to lose a bit of time in the selection +of sensors to activate. Indeed the optimization program allows to reduce +significantly the number of active sensors and also the energy consumption while +keeping a good coverage level. The energy overhead when increasing network size +is the lowest with PeCO. \begin{figure}[h!] \centering @@ -839,16 +843,17 @@ size is the lowest with PeCO. \subsubsection{Network Lifetime} -In comparison with the two other approaches, PeCO and DiLCO protocols are better for prolonging the network lifetime. In -Figures~\ref{figure9}(a) and (b), $Lifetime_{95}$ and $Lifetime_{50}$ are shown for -different network sizes. As can be seen in these figures, the lifetime -increases with the size of the network, and it is clearly larger for the DiLCO and -PeCO protocols. For instance, for a network of 300~sensors and coverage ratio -greater than 50\%, it can be observed on Figure~\ref{figure9}(b) that the lifetime is -about twice longer with PeCO compared to the DESK protocol. The performance -difference is more obvious in Figure~\ref{figure9}(b) than in -Figure~\ref{figure9}(a) because the gain induced by protocols (PeCO and DiLCO) increases with -time, and the lifetime with a coverage over 50\% is far longer than with 95\%. +In comparison with the two other approaches, PeCO and DiLCO protocols are better +for prolonging the network lifetime. In Figures~\ref{figure9}(a) and (b), +$Lifetime_{95}$ and $Lifetime_{50}$ are shown for different network sizes. As +can be seen in these figures, the lifetime increases with the size of the +network, and it is clearly larger for the DiLCO and PeCO protocols. For +instance, for a network of 300~sensors and coverage ratio greater than 50\%, it +can be observed on Figure~\ref{figure9}(b) that the lifetime is about twice +longer with PeCO compared to the DESK protocol. The performance difference is +more obvious in Figure~\ref{figure9}(b) than in Figure~\ref{figure9}(a) because +the gain induced by protocols (PeCO and DiLCO) increases with time, and the +lifetime with a coverage over 50\% is far longer than with 95\%. \begin{figure}[h!] \centering @@ -860,17 +865,21 @@ time, and the lifetime with a coverage over 50\% is far longer than with 95\%. \label{figure9} \end{figure} -Figure~\ref{figure10} compares the lifetime coverage of the DiLCO and PeCO protocols -for different coverage ratios. Protocol/70, Protocol/80, -Protocol/85, Protocol/90, and Protocol/95 correspond to the amount of time during which the +Figure~\ref{figure10} compares the lifetime coverage of the DiLCO and PeCO +protocols for different coverage ratios. Protocol/70, Protocol/80, Protocol/85, +Protocol/90, and Protocol/95 correspond to the amount of time during which the network can satisfy an area coverage greater than $70\%$, $80\%$, $85\%$, $90\%$, and $95\%$ respectively, where the term Protocol refers to DiLCO or -PeCO. Indeed there are applications that do not require a 100\% coverage of the -area to be monitored. For example, forest -fire application might require complete coverage -in summer seasons while only require 80$\%$ of the area to be covered in rainy seasons~\citep{li2011transforming}. As another example, birds habit study requires only 70$\%$-coverage at nighttime when the birds are sleeping while requires 100$\%$-coverage at daytime when the birds are active~\citep{1279193}. - PeCO always outperforms DiLCO for the three lower coverage ratios, moreover the -improvements grow with the network size. DiLCO outperforms PeCO when the coverage ratio is required to be $>90\%$, but PeCO extends the network lifetime significantly when coverage ratio can be relaxed. +PeCO. Indeed there are applications that do not require a 100\% coverage of the +area to be monitored. For example, forest fire application might require +complete coverage in summer seasons while only require 80$\%$ of the area to be +covered in rainy seasons~\citep{li2011transforming}. As another example, birds +habit study requires only 70$\%$-coverage at nighttime when the birds are +sleeping while requires 100$\%$-coverage at daytime when the birds are +active~\citep{1279193}. PeCO always outperforms DiLCO for the three lower +coverage ratios, moreover the improvements grow with the network size. DiLCO +outperforms PeCO when the coverage ratio is required to be $>90\%$, but PeCO +extends the network lifetime significantly when coverage ratio can be relaxed. \begin{figure}[h!] \centering \includegraphics[scale=0.55]{figure10.eps} @@ -887,16 +896,14 @@ hand, the choice of $\beta \gg \alpha$ prevents the overcoverage, and also limits the activation of a large number of sensors, but as $\alpha$ is low, some areas may be poorly covered. This explains the results obtained for $Lifetime_{50}$ with $\beta \gg \alpha$: a large number of periods with low -coverage ratio. On the other hand, when $\alpha \gg \beta$ is chosen, -the coverage is favored even if some areas may be overcovered, so a high coverage ratio is -reached, but a large number of sensors are activated to achieve this goal. -Therefore the network lifetime is reduced. The choice $\alpha=0.6$ and +coverage ratio. On the other hand, when $\alpha \gg \beta$ is chosen, the +coverage is favored even if some areas may be overcovered, so a high coverage +ratio is reached, but a large number of sensors are activated to achieve this +goal. Therefore the network lifetime is reduced. The choice $\alpha=0.6$ and $\beta=0.4$ seems to achieve the best compromise between lifetime and coverage -ratio. That explains why this setting has been chosen for the experiments +ratio. That explains why this setting has been chosen for the experiments presented in the previous subsections. - - \begin{table}[h] \centering \caption{The impact of $\alpha$ and $\beta$ on PeCO's performance} @@ -922,33 +929,31 @@ $\alpha$ & $\beta$ & $Lifetime_{50}$ & $Lifetime_{95}$ \\ \hline \section{Conclusion and Future Works} \label{sec:Conclusion and Future Works} -In this paper the problem of perimeter coverage optimization in -WSNs has been studied. A new protocol called Perimeter-based Coverage -Optimization is designed. This protocol schedules nodes' activities (wake up and sleep stages) with -the objective of maintaining a good coverage ratio while maximizing the network -lifetime. This protocol is applied in a distributed way in regular subregions -obtained after partitioning the area of interest in a preliminary step. It works -in periods and is based on the resolution of an integer program to select the -subset of sensors operating in active status for each period. This work is -original in so far as it proposes for the first time an integer program -scheduling the activation of sensors based on their perimeter coverage level, -instead of using a set of targets/points to be covered. Several simulations have -been carried out to evaluate the proposed protocol. The simulation results show -that PeCO is more energy-efficient than other approaches, with respect to -lifetime, coverage ratio, active sensors ratio, and energy consumption. - -This framework will be extented so that the schedules are planned for multiple -sensing periods. The integer program would be improved to take into -account heterogeneous sensors from both energy and node characteristics point of -views. Finally, it would be interesting to implement the PeCO protocol using a +In this paper the problem of perimeter coverage optimization in WSNs has been +studied. A new protocol called Perimeter-based Coverage Optimization is +designed. This protocol schedules nodes' activities (wake up and sleep stages) +with the objective of maintaining a good coverage ratio while maximizing the +network lifetime. This protocol is applied in a distributed way in regular +subregions obtained after partitioning the area of interest in a preliminary +step. It works in periods and is based on the resolution of an integer program +to select the subset of sensors operating in active status for each period. +This work is original in so far as it proposes for the first time an integer +program scheduling the activation of sensors based on their perimeter coverage +level, instead of using a set of targets/points to be covered. Several +simulations have been carried out to evaluate the proposed protocol. The +simulation results show that PeCO is more energy-efficient than other +approaches, with respect to lifetime, coverage ratio, active sensors ratio, and +energy consumption. + +This framework will be extented so that the schedules are planned for multiple +sensing periods. The integer program would be improved to take into account +heterogeneous sensors from both energy and node characteristics point of views. +Finally, it would be interesting to implement the PeCO protocol using a sensor-testbed to evaluate it in real world applications. \subsection*{Acknowledgments} -The authors are deeply grateful to the anonymous reviewers for their -constructive advice, which improved the technical quality of the paper. As a -Ph.D. student, Ali Kadhum Idrees would like to gratefully acknowledge the -University of Babylon - Iraq for financial support and Campus France for the -received support. This work is also partially funded by the Labex ACTION program +Ali Kadhum Idrees' PhD thesis is financially supported in part by University of Babylon (Iraq). +This work is also partially funded by the Labex ACTION program (contract ANR-11-LABX-01-01). \bibliographystyle{gENO}