The effect of temperature on the dynamics of a homogeneous oscillatory system operated in batch and under flow

The effect of temperature on a network of chemical reactions is not obvious, especially when compared to the effect exerted on elementary steps. There are few reports regarding the estimation of parameters such as activation energies for oscillating chemical systems. Still less investigated is the importance of the relative distance from thermodynamic equilibrium on the way in which temperature influences the oscillators' dynamics – a crucial aspect for the understanding of chemical and bio-chemical oscillating networks. In this paper we use the bromate–oxalic acid–acetone–cerium oscillatory system to study the influence of temperature under close and far-from-equilibrium regimes. The research was carried out under identical conditions for batch and flow (in a continuous flow stirred tank reactor, CSTR) regimes, and the main oscillation features were preserved, so that it was possible to isolate the effect of flow. Overall, increasing the flow results in an increase of the oscillatory frequency. The apparent oscillatory activation energy was found to decrease from 72 ± 6 kJ mol−1, for the system operated in batch, to 50 ± 2 kJ mol−1, under the flow regime. The role of the distance from the thermodynamic equilibrium on the temperature dependence is generalized and discussed in connection with other systems. Numerical simulations using the Brusselator model under batch and flow regimes further helped the discussion of the main experimental results.


Introduction
A network of chemical reactions under oscillatory regime might display a rather unusual response to temperature changes. Examples include non-Arrhenius behaviors such as temperature compensation and overcompensation. [1][2][3][4][5] These infrequent effects play a decisive role, for instance, in temporally organized patterns found in living systems which are commonly known to be rather insensitive to temperature and other parameters such as pH, etc. Temperature compensation is particularly important as part of the wide-ranging homeostatic mechanism, and it underlies many oscillatory mechanisms in living systems (cf. circadian, ultradian and some neuronal rhythms), cf. ref. 6 and references therein. In terms of individual reactions, temperature compensation results from the balance among the weighted activation energies of frequency-increasing and frequency-decreasing elementary steps in a reaction network. 1,2,7 Despite the comparable simplicity of chemical oscillators, results obtained in such designed in vitro systems can provide valuable information on some structural aspects of more complex and sometimes less tractable biochemical oscillators. [8][9][10] Indeed, the observation of temperature compensation in comparable simpler (electro)chemical oscillators 5,11,12 attests their importance as model systems. In most cases, the effect of temperature in these oscillators has been discussed in terms of the temperature dependence of the oscillatory frequency, as rstly suggested by Körös. 13 One of the most prominent features of living systems is that their internal organization is kept at the expense of exporting entropy to the environment. In order to stay alive, these open systems import (or are fed with) high free energy materials and reject the degraded ones. Already in 1950, von Bertalanffy 14 recognized that "from the physical point of view, the characteristic state of the living organism is that of an open system", and as such, living systems can thrive by "maintaining themselves in exchange of materials with environment, and in continuous building up and breaking down of their components". Therefore, the position with respect to the thermodynamic equilibrium is a key parameter to the understanding and description of the actual state of open systems. In chemical systems, the departure from the state of thermodynamic equilibrium can be tuned by the rate at which reactants are fed and products (and also unreacted species and intermediates) are removed from the reaction vessel. Tuning this parameter results on non-trivial changes in most reaction rates. Hence, the effect exerted by temperature is expected to depend on the distance from thermodynamic equilibrium.
A literature survey reveals a lack of reports dealing with the inuence of the relative position with respect to thermodynamic equilibrium on the way in which temperature inuences the oscillators' dynamics. In a seminal report, Rössler 15 conjectured on the importance of ow when designing a temperature-compensated homogenous chemical oscillator. The discussion was carried out in terms of a theoretical model for a relaxation oscillator. Insofar, however, it seems that no verication of this expectation has been provided. Ruoff and coworkers 16 studied the effect of temperature on the dynamics of some bromate oscillators catalyzed by cerium, and reported a considerable decrease in the activation energy when the system was operated in a continuous ow stirred tank reactor (CSTR), when compared with the batch regime. The effect of ow and/or temperature on the dynamics has been reported for some oscillatory systems such as the hydrogen peroxide-thiosulfate-sulte ow system, 11,17 pH oscillators, 12 and the Bray reaction. 18 We have recently studied 19 the time evolution of the apparent oscillatory activation energy for a BZ-like oscillator operated in batch as it approaches the thermodynamic equilibrium. As the main result, we observed that the apparent oscillatory activation energy increases as reactants are consumed and products accumulate inside the reactor. A detailed investigation of formic acid electro-oxidation on platinum 5 revealed the predominance of highly non-Arrhenius temperature dependence over a wide parameter window. The system was characterized in terms of the distance from the thermodynamic equilibrium and temperature compensation was found to prevail at considerably high applied currents or, equivalently, at the farthest distance from equilibrium.
The present contribution results from our interest to understand the combined effect of temperature and distance from the thermodynamic equilibrium. We decided to undergo such an investigation using the bromate-oxalic acid-acetonecerium system 20,21 as a model system. This system was chosen mainly because of its robustness and low sensitivity of its properties, such as the oscillation morphology, on temperature and ow rates. The constancy of such features might indicate that the underlying chemistry remains unchanged at distinct ow rates and temperatures. As far as the effect of temperature is concerned, this system has been studied under both batch and ow regimes. 22,23,24,25 In contrast to previous reports by other authors, however, the systematic experiments reported here allowed at isolating as much as possible the contribution of the ow on the temperature dependence. We initially optimized the system in order to compare the effect of temperature in batch and ow. In this way, the open and closed systems were studied under identical conditions of stirring rate, temperature, and concentrations. In order to go deeper in the understanding of the role played by the distance from thermodynamic equilibrium on the system's dynamics, we carried out numerical simulations using the Brusselator model under batch and ow regimes. Finally, the results are discussed in connection with other processes in electro-, bio-and chemical systems.

Experimental section
The cylindrical glass reactor employed in all experiments is 90 mm high and 47 mm in diameter. The temperature was controlled by water circulation through its glass jacket, with the aid of a thermostatic bath (Microquimica, model  and monitored using a digital thermometer (Tec-Lab). The Teon reactor cap contains the holes to allow the use of a platinum electrode and a reversible hydrogen electrode, lled with aqueous sulfuric acid solution with concentration identical to that found in the reactor. The potential difference was measured with a multimeter (Minipa, ET2201). Further experimental details can be found elsewhere. 19,26 For open reactor experiments (CSTR) two peristaltic pumps (Milan, 640) were used. The system was fed with three aqueous H 2 SO 4 solutions: solution A (CH 3 COCH 3 and (COOH) 2 ), solution B (NaBrO 3 ), and solution C (Ce 2 (SO 4 ) 3 ). In this case, all solutions were placed in jacketed asks and kept inside a thermostatic bath to assure the same temperature as that inside the reaction vessel. The solution inside the reactor was mixed by the use of a magnetic stirrer (Marconi, MA089) and a magnetic Teon bar 15 mm long and 5.7 mm of diameter. All experiments were carried out with a stirring rate of 700 rpm.
The oscillating frequency, u, for batch experiments was taken as an average for the high amplitude set of oscillations, eventually discarding the rst or the last value, if it is discrepant. The frequency of low amplitude oscillations was not considered because this regime was shown to be irregular.

Experiments
Oscillations in batch. The effect of temperature was rstly investigated in the system operated in batch. Fig. 1 shows the time-series for the bromate-oxalic acid-acetone-cerium system operated in batch at different temperatures from 5 to 35 C. As a rst noteworthy aspect, it is seen that temperature affects only slightly the morphology of the main, high amplitude, oscillations, suggesting that, in the range investigated, there is no change in the mechanism underlying the oscillatory dynamics.
In contrast to what is observed when Ce 4+ ions are present, there is no induction period when Ce 3+ is used. 26 Moreover, the use of Ce 4+ results in considerable scattering in the Arrhenius plots. Therefore, all experiments presented here were carried out with Ce 3+ , performed in triplicate, and all time-series presented are representative of the system's dynamics and reproducible.
At 35 C the system did not oscillate under the given conditions. The addition of a small aliquot (0.37 mL) of acetone at the instant indicated by the vertical line, produced the emergence of transient oscillations, as shown in Fig. 1(d). To circumvent the absence of oscillations at 35 C, we performed experiments with slightly higher acetone concentration, namely 0.157 mol L À1 , and the results are presented in Fig. 2.
Overall, the high amplitude oscillations presented in Fig. 1  and 2 have comparable waveform and amplitude. Only a small change in the oscillatory frequency is observed for distinct acetone concentrations. The initial and uniform high amplitude oscillations were used to estimate the apparent oscillatory activation energy. Fig. 3 shows the Arrhenius plot for the data presented in Fig. 1 and 2, and illustrates the effect of temperature on the dynamics of the bromate-oxalic acid-acetonecerium oscillatory system. The, apparent oscillatory activation energies, E u , calculated in these plots 13 were found to amount to 72 AE 2 kJ mol À1 (for [CH 3 COCH 3 ] ¼ 0.115 mol L À1 ), and to 74 AE 1 kJ mol À1 (for [CH 3 COCH 3 ] ¼ 0.157 mol L À1 ). Experiments were done in triplicate and these results represent the averaged values. Pastapur and Kulkarni 22 studied the effect of temperature in the bromate-oxalic acid-acetone system in batch using Ce 3+ and Mn 2+ , as catalysts. The study was performed at 25, 30, 35, and 40 C and the activation energy obtained under oscillatory regime for the system with Ce 3+ was 73 or 82 kJ mol À1 , depending on the method used to estimate. The discrepancies between these values and the ones we found in the present work can be attributed to the differences in the concentrations and temperature range used in each case.
An important remark to be done at this point is that although no oscillations were found in batch at 35 C for the lower acetone concentration used, the apparent oscillatory activation energy depicted in Fig. 3 ($72-74 kJ mol À1 ) is representative of the system operated in batch in the temperature range investigated. The ow experiments were carried out with [CH 3 COCH 3 ] ¼ 0.115 mol L À1 and the main results are summarized in the following.
Experiments in a CSTR. As already mentioned, one of the aims of the present work consists in the investigation of the effect that temperature exerts on the systems' dynamics under different regimes, namely the distance from thermodynamic equilibrium. Therefore, the study in the open system started with the investigation of the ow rate effect on the dynamics. The oscillations' morphology and amplitude were barely affected by the ow rate in the range investigated. The effect of the ow rate on the oscillatory dynamics at 25 C is displayed in Fig. 4. The vertical lines indicate the instant at which the ow was interrupted.

Fig. 3
Arrhenius plots using the oscillatory frequency u, for the bromate-oxalic acid-acetone-cerium system operated in batch at two distinct initial concentrations of acetone: 0.115 mol L À1 (black squares) and 0.157 mol L À1 (red circles). Other concentrations: s À1 ), and 78 s (at k 0 ¼ 0.020 s À1 ). Fig. 5 summarizes the effect of ow rate on the oscillatory frequency for all experiments. As can be seen, the oscillatory frequency increases, almost linearly, with the ow rate. Data points in batch, i.e. k 0 ¼ 0, account for three different experiments.
Our results are in line with that reported by Pereira and Faria, 21 who observed a comparable behavior and reported a decrease of about 20% on the oscillation period for an increase in k 0 from 0.0119 to 0.0352 s À1 , for the oscillatory bromateoxalic acid-acetone-cerium system. General increase of the oscillatory frequency when increasing ow rates has been also found for the hydrogen peroxide-thiosulfate-sulte ow system. 11 The difficulty in comparing previously published data is due to the commonly observed problem that the system is studied under different conditions (mainly concentration of chemicals) in batch and ow. As already mentioned, this is not a problem in the present case.
In homogeneous systems, the decrease in oscillatory frequency, which follows the reduction of ow rate can be tentatively interpreted as a damping caused by a decrease in the concentration of reacting species and the corresponding reduction of reaction rates of all reactions. As the ow rate increases, the concentrations of reacting species inside the reaction vessel are maintained high because of the efficient replenishment of fresh reactants, and also the removal of products.
For most of the homogeneous systems operated in batch, the degree of decrease of the oscillatory frequency in time reects the strength of the damping; see below more on this aspect. We have explained this phenomenon before in somewhat different terms, 19 but the concept is equivalent in both cases. This is a known effect that calls for a special care on the denition of the procedure for estimating the activation energy under oscillatory regime. 27 Moreover, the damping in oscillatory frequency, observed under batch conditions, is also accompanied by a decrease in amplitude. This interpretation is valid for the same type of oscillations, and thus does not include sequential oscillations witnessed in batch. 19,[28][29][30][31] In fact, sequential oscillations are dynamic states located at different regions in the phase space that are visited as an uncontrollable parameter slowly varies. 32,33 As such, there is no particular trend to be expected for the oscillations' frequency and amplitude in sequential oscillations.
When compared to chemical systems, the relative position with respect to the thermodynamic equilibrium in electrochemical systems can be controlled either via potentiostatic or galvanostatic mode. 34 In both cases, the parameter, i.e. the potential or the current, respectively, readily informs on the distance with respect to the thermodynamic equilibrium. Nevertheless, the departure from the equilibrium is even more transparent since the ow of electrons through the interface is controlled. In agreement with the present ndings, the oscillatory frequency also increases with the applied current for the electro-oxidation of formic acid on platinum 5 and on platinumtin surfaces. 35 Again, these are examples of a situation where the same kind of oscillations are compared at different applied currents. As a nal observation on the comparison between chemical and electrochemical systems in the context of the present discussion, it is remarkable the peculiarity of chemical systems in the sense that the situation for k 0 ¼ 0 cannot be obtained in electrochemical systems. Consequently, the range of k 0 explored in Fig. 5 is in principle broader than that registered in electrochemical systems. Even considering the problem in properly normalizing the distance from thermodynamic equilibrium, this fact allows for further generalizations.
The effect of temperature was also investigated using a system with identical composition as in Fig. 4, but with k 0 ¼ 0.0098 s À1 . Fig. 6 shows the effect of temperature on the timeseries for the bromate-oxalic acid-acetone-cerium oscillatory system under ow regime. The vertical lines indicate when the ow was stopped. Aer interrupting the feeding, the system oscillates for a while, just as observed when operated in batch, vide supra. In agreement with the results at 35 C for the system operated in batch for this acetone concentration, the oscillations cease earlier when the temperature is increased, and as shown in Fig. 6(d). As it was observed in batch, a regular, Fig. 4 The influence of the flow rate on the dynamics of the bromateoxalic acid-acetone-cerium oscillatory system at 25 C.

View Article Online
Arrhenius-like behavior is also found for the ow regime and both oscillations' waveform and amplitude remain very similar at all temperatures. The effect of temperature on the oscillatory frequency under both batch and open regimes is summarized in Fig. 7(a). The missing point in the data for the batch system accounts for the fact that no oscillations are observed at 35 C for the lower acetone concentration, cf. Fig. 3. The oscillatory frequencies at different temperatures taken from Fig. 6 were used to estimate the activation energy for the CSTR regime. 19 The Arrhenius plots for the CSTR regime are presented in Fig. 7(b). An excellent linearity is observed and an apparent oscillatory activation energy of 50 AE 2 kJ mol À1 results.
As already anticipated in Fig. 5, it is clear from Fig. 7(a) that the open system oscillates faster than the one operated in batch, for the whole range of temperature studied. The apparent oscillatory activation energy decreases from 72 AE 2 kJ mol À1 , in batch, to 50 AE 2 kJ mol À1 under ow. As the batch regime is virtually at k 0 ¼ 0, if other k 0 values higher than 0.0098 s À1 were employed an oscillatory activation energy still lower should be nd. These two aspects are central in the present study and will be further discussed below.
Experimental results for other systems are in line with our ndings, as already mentioned in the Introduction. 16,19 Nogueira et al. 19 in particular, clearly show the increase of the apparent oscillatory activation energy as the system approaches the thermodynamic equilibrium. As already stressed, in order to isolate as much as possible the effects of interest, namely the operation mode and temperature, it is mandatory to keep, if possible, all the other parameters constant, as in the present case. In this respect, we are not aware of a comparable report. Taking specically the work by Kulkarni and co-workers, [22][23][24][25] who investigated the effect of temperature under comparable conditions for the same system studied here, besides the different concentrations used in some experiments in batch and under ow, the authors also used Ce 4+ as the catalyst in some experiments, which may not be appropriate. In fact, we have shown that Ce 3+ is the needed catalyst to observe oscillations in the bromate-oxalic acid-acetone-cerium system. 21,26 In batch, no induction period is observed when Ce 3+ is used, and in a CSTR, oscillations were found with Ce 3+ and with aged solution of Ce 4+ , aer partial conversion of Ce 4+ to Ce 3+ by reaction with acetone. The use of Ce 4+ results in long induction times and leads to imprecision. In terms of the inuence of temperature, we have observed that the Arrhenius plots are less linear when Ce 4+ is used instead of Ce 3+ .
Back to the main results, in summary, we observed that, as the system is moved out from the state of thermodynamic equilibrium, its oscillatory frequency generally increases and its apparent oscillatory activation energy decreases. The rst aspect, the frequency increasing with the distance from equilibrium was already discussed. We focus now on the temperature dependence, as discussed in terms of apparent activation energy for the oscillatory system. The activation energy of an elementary reaction step accounts for the energy barrier that molecules have to overcome in order to form products. For more complex processes, however, the presence of different steps and intermediates results in the existence of several energy barriers. In the case of an oscillating network consisting of many reactions and chemical species, the apparent  Fig. 7 (a) The effect of temperature on the oscillatory frequency u, for the bromate-oxalic acid-acetone-cerium system operated in batch and also in a CSTR. (b) Arrhenius plots using the oscillatory frequency, u, for the system operated in a CSTR. Data for experiments in batch as in Fig. 1 and for the open system as in Fig. 6. oscillatory activation energy, as measured by the oscillatory frequency, informs on the temperature dependence of the oscillator, as it was based on a single reaction. As presented by Körös, 13 the system can be regarded as 'a series of autocatalytic reaction bursts occurring with a certain frequency'. In our case, the apparent activation energy decreases as the system is moved away from equilibrium, which is equivalent to saying that it becomes less sensitive to temperature changes. As a consequence, the relative independence of the physiological rhythms in living systems on the environmental temperature, known as temperature compensation, might indeed be favored under signicantly far from equilibrium conditions. It is well-known that dynamic self-organization, 36 such as temporally organized rhythms, that characterizes all living structures, occurs only in open and far from equilibrium systems. What our results reinforce is that the distance from equilibrium favors the decrease of temperature sensitivity.
Curiously enough, but also in line with our observations, our previous study on the electro-oxidation of formic acid on platinum and in acidic media, 5 showed that, even for a system having a highly non-Arrhenius temperature dependence, temperature compensation is found under considerably far from equilibrium regime.
The role of ow in the temperature compensation has been explored by Otto Rössler in an oen cited proceeding paper. 15 The author discusses the temperature-sensitivity as resulting of a combination of kinetic and ow parameters and draw some particular conditions for achieving temperature compensation for the case studied, namely a relaxation oscillator. Importantly, the author mentions the role eventually played by membranes. In biological organisms, membranes are key components that control the inow and outow of chemicals and preserve the internal organization by keeping the system under far from equilibrium regime. However, besides the magnitude of ow itself, the selectivity, or its effect on particular species, is also a key aspect when discussing properties such as the temperature dependence on the oscillatory frequency. This aspect will be further discussed below.
In the following we present some numerical results that explore some aspects raised by our experiments. Given the lack of a well-established model for the studied system, we used a general and very simple model to investigate the effect of ow and temperature.

Numerical solutions of the Brusselator model
The chosen model to test the ideas just discussed was the Brusselator, 37,38 a very simple, classical example of a limit cycle oscillator. The model consists of the following hypothetical steps: where A and B are the initial reactants whose concentrations are assumed to be constant, D and E are nal products, X and Y are intermediate species that oscillate under a suitable choice of the rate of reactions and initial concentrations of A and B. This mechanism can be translated to a set of differential equations: These equations describe how the concentrations of the two intermediate species X and Y change as a function of their concentration and that of the initial reactants A and B. In the original formulation of the Brusselator model, 37,38 [A] and [B] are held constant which implies continuous inow for those species, and the system represents, thus, a semi-batch oscillator. In our case, the original equations were modied to directly incorporate the effect of reactant inow and outow uxes, as we shall describe in the following. Specic routines were built using the Mathematica® soware to solve the resulting ODEs as well as analyze the resulting time series as a function of selected parameter changes. Runge-Kutta with variable step size was the preferred method for numerical integration.
The effect of the ow. In order to explicitly incorporate the effect of the ow of reactants in a continuously fed well-stirred tank reactor (CSTR) it is necessary to add two other equations accounting for the addition and removal of the initial reactants A and B, where k f stands for the ow rate: The nal set of equations describing the CSTR system can be achieved with a minor intervention on the original ODEs to account for the removal of the intermediate species. The nal set of equations for the CSTR oscillator then becomes: Oscillations may be found for this system of equations aer numerical integration and a sensible choice of the parameters. In what follows, the values of A 0 and B 0 were xed at 20, and the oscillatory frequency analyzed as a function of the ow rate, k f , for some selected sets of the constants {k 1 , k 2 , k 3 , k 4 } which were representative of the overall model behavior. Numerical results of the effect of the ux of reactants through the reactor on the oscillatory frequency are displayed in Fig. 8.
The majority of curves presented in Fig. 8 reveal a general increase in the values of oscillatory frequency as the ow rate increases, in line with our experimental ndings, cf. Fig. 5. This behavior was found for six combinations of the rate constants {k 1 , k 2 , k 3 , k 4 } and for the whole range of ow in which oscillations exist. One exception was found for the set {1, 2, 1, 3}, but the deviation occurs only for a very high k f , namely k f higher than 0.5 s À1 , i.e. in an extreme parameter region. In terms of the mechanism, the decrease of oscillatory frequency when increasing the ow rate was found only for a combination of the smallest value for the rate of step 1 (A / X) with the highest value for the rate of step 4 (X / E), and only for very high ow rates. Under these conditions the decrease of the oscillatory frequency can be attributed to the small concentrations of X (low formation and fast removal) and also to its low residence time within the reactor. Finally, we observed that the range of ow rate in which oscillations can be found generally increases with the decrease of k 3 and with the increase of k 4 . It was mentioned in the experimental section that, in batch, a possible explanation for the slower dynamics was the general decrease in the concentration of active species caused by their consumption in a closed system. This effect is absent in our simulations since we used constant concentrations for A and B in the simulations for the original Brusselator, which prevents any damping and produces only stable oscillations as for the system under explicit ow.
As seen in the set ODE (3), the ow rate parameter k f is found to impact several terms of the differential equations' set. It is possible to evaluate how each one of these terms affects the overall oscillatory frequency. To do so, we can discriminate each k f,i which appear in the set ODE (4): Altogether, there are six ow rate terms, four of them which increase (k f1 and k f3 ) and decrease (k f2 and k f4 ) the concentrations of A and B, and two related with the withdraw of the intermediate species X and Y (k f5 and k f6 ). We can vary each of those k f,i individually while keeping the remaining values constant at a control value. It is important to note, at this point, that this is merely a mathematical procedure, aiming to explore how these parameters, individually, inuence the overall oscillatory dynamics for the model, and that there is no direct physical counterpart associated to a variation on a specic value of k f,i , since it would violate volume conservation. Still, one can see the variation in those parameters as a mean of increasing or decreasing the relative abundances of the reactants A, B, X, and Y. With this thought in mind, and the tool in hands, important insights can be achieved as of the role of availability of each of the reactant and its effect on oscillatory properties. Fig. 9 shows the effect of the variation of each k f,i value on the oscillatory frequency for the CSTR version of the Brusselator. On panel (a) the results of such procedure are shown for values of k f,i starting from a low ow rate value (k 0 ¼ 40), while (b) displays the results for a high ow rate value (k 0 ¼ 200). The set of four k i values used were {1, 2, 1, 2}, but the general trend is similar for other sets of k i . Fig. 9 shows that each different constant is associated with different effects on the oscillatory frequency. A common feature, however, is that a given rate constant plays the same role either at low or high ow rate. It is possible to quantify the magnitude of individual ow rate parameters k f,i on the overall behavior of frequency by introducing a sensitivity coefficient. Such coefficients have been used in biological contexts to evaluate the response of specic living organelles' functions for variations of important chemical species or physical parameters such as temperature. 39 In our specic case, it is intuitive to introduce the sensitivity coefficient for the oscillatory frequency: here, a specic sensitivity coefficient S k f,i u describes how the oscillatory frequency u changes as function of a variation on a specic ow rate parameter k f,i . The partial derivative is evaluated in the vicinities of the control ow rate, rendering thus a special physical meaning to the sensitivity coefficients, since the change in volume of the reactor is innitesimal for that limit. Table 1 presents the overall dependency of the oscillatory frequency on the increase in each of the k f,i constants as well as the specic control coefficients. It can be seen that parameters k f1 and k f2 related with the addition and removal of A, respectively, increase and decrease the oscillatory frequency. On the other hand, parameters k f3 and k f4 , associated with the addition and removal of B, display the opposite effect. This result is in accordance with the simulated original Brusselator system where integration of the original ODE (1) set reveals that an increase in the concentration of A leads to higher oscillatory frequencies, while an increase for B promotes the opposite trend. Finally, Table 1 shows that an increase in the values of k f5 and k f6 promotes a decrease in oscillatory frequency, which means that depletion on the concentration of the intermediates X and Y has a negative effect on this property.
We observe that the individual ow rate parameter operates on several terms of the set of differential equations and that each one of them has a particular effect on the oscillatory frequency, since it is specically associated to the increase/ withdraw of a particular reactant. Back to our specic results, a closer look in Fig. 9 reveals that the constants k f1 and k f3 related to the addition of reactants A and B, respectively, display the most pronounced effect on frequency, with k f1 promoting its increase while k f3 producing its decrease. Thus we can explain the overall increase of frequency as the ow rate is increased (curve in black in Fig. 9), as a function of a higher weight of the term related with k f1 . In summary, the results shown in Fig. 9 can be rationalized as the net increase in the oscillatory frequency is a result of an increasingly surplus of the initial reactant A.
As represented by the black curve in Fig. 9, simple experiments in a CSTR would result in overall increase of the oscillatory frequency with the increase in the ow rate. The selective removal or addition of a given species in the reaction network would result in the increase or decrease of the oscillatory frequency. An interesting aspect that emerges from these results consists of the opposite effects exerted by the removal or addition of a given species. This feature opens the possibility of exploring the selective role of a given species to achieve properties such as temperature compensation. Indeed, the high selectivity of cellular membrane might play a role in keeping biological rhythms unaltered despite of changes in the environment. The importance of the interactions between membrane ion transport and ion concentration gradients in circadian rhythms has been discussed in the Njus-Sulzman-Hastings or membrane model. 40 Despite the criticisms to some aspects of this earlier model, Nitabach et al. 41 stressed the essential role that the ion uxes through the membrane play in the mechanism of the core oscillator. Therefore, this discussion suggests that, instead of the robustness of a given biochemical network, it would be thus more adequate to consider the conned system (chemistry and membrane) as whole. The use of methods to investigate the reaction mechanism in complex systems 42,43 is certainly of help in this direction.
The effect of temperature. Temperature variation will affect directly the overall oscillatory frequency through its inuence on each specic reaction rates k i . This effect can be mathematically treated through the use of control coefficients: 1-4 The overall effect of temperature on frequency can be understood as the sum: Through reference to Arrhenius' law, the following equation can be derived: here it is possible to see that each specic reaction will inuence the overall frequency change with respect to its specic control coefficient and activation energy. The overall activation energy can be calculated as: For the Brusselator model, this expression actually translates into: which indicates how the apparent oscillatory activation energy relates with the activation energies of individual steps. Fig. 10 portraits the inuence of the reaction rates on the oscillatory frequency for the CSTR version of the Brusselator model. Results in Fig. 10 describe the impact of the rate constants on the oscillatory frequency. Coefficients C ki u were estimated and indicate that: steps 1 (A / X) and 3 (Y + 2X / 3X) contribute to increase the oscillatory frequency as the temperature increases, whereas steps 2 (B + X / Y + D) and 4 (X / E) belong to the frequency decreasing set. In short, steps that produce X increase the oscillatory frequency and the ones that consume X decrease the oscillatory frequency. In general, this trend was also found for almost all combinations of constants investigated. According to eqn (e6), the apparent oscillatory activation energy for the Brusselator model would thus result of positive contributions of the activation energies of steps 1 and 3, and negative ones of steps 2 and 4. The impact of individual rate constants on the oscillatory frequency can be readily compared to the role of selective feeding of different species discussed above. For instance, saying that step 1 (A / X) contributes to increase the oscillatory frequency as the temperature increases, as evidenced by its positive control coefficient, is equivalent to the observation that the addition/removal of species A contributes to the increase/decrease of the oscillatory frequency, as discussed above.

Summary and conclusions
We investigate in this article the coupled effects of temperature and of the distance from thermodynamic equilibrium on the dynamics of a homogeneous oscillatory system. Experiments were carried out with the bromate-oxalic acid and acetonecerium reaction under different temperatures and ow rates, by means of a continuous ow stirred tank reactor (CSTR). The conclusions on the effect of the distance from thermodynamics equilibrium were validated by the nearly identical oscillation's features under both batch and ow regimes. The experimental results were discussed in connection with numerical simulations using the Brusselator model, in its conventional formulation and also in a modied one to explicitly include the effect of ow.
The experiments reveal a general increase of the oscillatory frequency when the ow rate is increased, in line with our previously published experiments in some homogeneous (chemical) and also heterogeneous (electrochemical) oscillators. Numerical simulations with the CSTR version of the Brusselator conrm this tendency but also disclose a tiny parameter region (for one specic set of rate constants and at extremely high ow rates) in which the increase of the ow rate results in a decrease in the oscillatory frequency. The effect of ow rate was also explored for each species of the Brusselator. Remarkably, we observed that, although the increase of the ow rate typically increases the oscillatory frequency, the feed concentration of each species impacts in a particular way the overall oscillatory frequency. This nding further stresses the role of ow rate in reacting systems, an important aspect for engineering articial oscillators for specic applications.
The apparent oscillatory activation energy obtained via Arrhenius plots using the experimental oscillatory frequency, was found to amount to 72 AE 6 kJ mol À1 , in batch, and to 50 AE 2 kJ mol À1 , under ow. Unlike previous reports, these results were obtained under identical concentrations of all species and thus allowed to separate the effects of ow and temperature dependence. The decrease of the apparent oscillatory activation energy from batch to ow was interpreted as a result of the maintenance of the concentration of reactants, in contrast to their decrease caused by consumption in batch, which is in line with our previous results. 19 Numerical simulations with the CSTR version of the Brusselator were explored in terms of the control coefficients, vln u/vln k i , 1-4 for each reaction step i of the model. The results uncover the existence of two steps that are associated to the frequency-increasing set and two to the frequency-decreasing one, whose role remains unchanged despite the operational conditions. Furthermore, and as already mentioned, the particular manner at which the feeding rate of each species affects the oscillatory frequency further reinforces the role played by the selectivity of biological membranes in this respect.