Black spots in a surfactant-rich Belousov Zhabotinsky reaction dispersed in a water-in-oil microemulsion system

Size: px
Start display at page:

Download "Black spots in a surfactant-rich Belousov Zhabotinsky reaction dispersed in a water-in-oil microemulsion system"

Transcription

1 THE JOURNAL OF CHEMICAL PHYSICS 122, Black spots in a surfactant-rich Belousov Zhabotinsky reaction dispersed in a water-in-oil microemulsion system Akiko Kaminaga, Vladimir K. Vanag, and Irving R. Epstein Department of Chemistry and Volen Center for Complex Systems, MS 015, Brandeis University, Waltham, Massachusetts Received 23 November 2004; accepted 16 February 2005; published online 3 May 2005 The Belousov Zhabotinsky BZ reaction dispersed in water-in-oil aerosol OT AOT microemulsion has been studied at small radius R d of water nanodroplets (R d nm 0.17, = H 2 O / AOT =9). Stationary spotlike and labyrinthine Turing patterns are found close to the fully oxidized state. These patterns, islands of high concentration of the reduced state of the Ru bpy 3 2+ catalyst, can coexist either with black reduction waves or, under other conditions, with the white oxidation waves usually observed in the BZ reaction. The experimental observations are analyzed with the aid of a new Oregonator-like model and qualitatively reproduced in computer simulations American Institute of Physics. DOI: / I. INTRODUCTION The Belousov Zhabotinsky BZ reaction, 1,2 the catalytic periodic oxidation of malonic acid MA by bromate in acidic aqueous solution, is the prototype nonlinear reaction for studying chemical oscillations in a stirred tank reactor and waves in spatially extended systems gels, 3 membranes, 4 mesoporous glasses, 5 and ion-exchange resins 6. The BZ reaction dispersed in water-in-oil aerosol OT microemulsion BZ-AOT system has been shown to be an excellent system for studying pattern formation in reaction-diffusion systems, including stationary Turing patterns, 7 oscillons, 8 and new wave patterns such as packet waves and antispirals, 9,10 dash waves, and segmented spirals. 11,12 In general, global Turing patterns 13 are of three main types: hexagons, stripes or labyrinthine patterns, or honeycombs or reverse hexagons. All of these patterns have been observed in the chlorite-iodide-malonic acid CIMA reaction in a gel reactor, 14,15 but until now only hexagons and labyrinthine patterns have been found in the BZ-AOT system. 7 Hexagons consist of white spots, since the patterns are typically recorded at the maximum of the absorption band of the reduced state of the catalyst ferroin or Ru bpy 3 2+, and the spots are islands of the oxidized catalyst, which has significantly less absorption and consequently transmits more light than its surroundings. Black spots, stationary islands of the reduced catalyst, patterns inverse to the white-spot hexagons, have not previously been observed in the BZ-AOT system. An AOT sodium bis 2-ethylhexyl sulfosuccinate microemulsion consists of nanometer-sized water droplets dispersed in a continuous oil octane phase. Each droplet is surrounded by a surfactant AOT monolayer. The radius of the water droplet cores is determined by the molar ratio = H 2 O / AOT and is roughly equal to R w =0.17 nm. 16,17 Most of our previous experiments were performed at = In the present set of experiments we explore smaller water droplets that form in surfactant-rich emulsions with =9. Water molecules in such droplets are adjacent to the polar SO 3 2 groups of the AOT surfactant molecules and therefore have special physical properties e.g., high viscosity, small dielectric constant. 18 Chemical reaction rates generally change in such an environment. In preliminary experiments with the BZ-AOT system in a stirred reactor CSTR, we found unusual oscillations at =9, with spikes in the concentration of oxidized catalyst ferriin during the autocatalytic stage of the oscillations. 19 A decrease in also implies an increase in the volume fraction of surfactant s at constant volume fraction of water w since s 21.6 w /. As a consequence, there will be a repartitioning of bromine which acts as a fast-diffusing inhibitor in the BZ-AOT system between the continuous oil phase and the nanodroplets. These observations suggest that new patterns may occur at small. In this paper we present the first experimental and theoretical results on black spots in the BZ-AOT system that can be observed either alone or accompanied by waves black or white. In Sec. II, we briefly describe our experimental and theoretical methods. In Sec. III, we report our experimental results. In Sec. IV, we suggest a new model of the BZ-AOT system that qualitatively captures the main experimental results found in this work. In Sec. V, we discuss plausible explanations for black spots and their interactions with waves and finally link these patterns to Turing patterns that can be altered by the excitability of the BZ-AOT system. II. MATERIALS AND METHODS The experimental arrangement and methods are described in detail elsewhere. 7,10,11 Two stock microemulsions ME, with the same molar ratio and volume droplet fraction d d = s + w, were prepared by mixing 1.5 M AOT in octane with aqueous solutions of malonic acid and H 2 SO 4 ME-1 or with NaBrO 3 and Ru bpy 3 2+ ME-2. The final reactive ME was obtained by mixing ME-1 with ME-2 volume ratio 1:1 and with enough additional octane to obtain the desired d. We worked mostly with microemulsions below the percolation level, as confirmed by measurements of /2005/ /174706/11/$ , American Institute of Physics

2 Kaminaga, Vanag, and Epstein J. Chem. Phys. 122, the electrical conductivity YSI The radius distribution of the water droplets was measured by quasielastic light scattering Protein Solutions, DynaPro. A small drop 0.1 ml of reactive ME was sandwiched between two flat glass windows separated by an annular teflon gasket Zefluor membrane of thickness 80 m, inner diameter 25 mm, and outer diameter 47 mm. The edges of the glass windows were sealed with teflon tape, and then the two windows were pressed together. This sandwich constitutes our reactor. Patterns in the reactor were observed through a microscope SZH10, Olympus in transmitted light passed through a 450 nm interference filter. The microscope is equipped with a charge-coupled device CCD camera TM-7CN, Pulnix connected to a personal computer. Snapshots and movies were recorded and analyzed with IMAGE- PRO PLUS software. The wavelength of a pattern was determined from the maximum of its fast fourier transform FFT. Numerical simulations of ordinary differential equations and reaction-diffusion equations were performed with the software FLEXPDE, 20 in which a Newton Raphson iteration process is used with a variable time step and mesh. FLEXPDE refines the triangular finite element mesh and/or time step until the estimated error in any variable is less than a specified tolerance, which we chose as 10 4 or in some cases as 10 5 at every cell of the mesh. The accuracy of FLEXPDE at a chosen tolerance was checked by comparing the parameters found for the onset of Hopf and/or Turing instability with those obtained analytically from linear stability analysis, as well as by varying the tolerance. III. EXPERIMENTAL RESULTS A. Turing patterns An AOT microemulsion loaded with the components of the BZ reaction sulfuric acid, sodium bromate, malonic acid, and Ru bpy 3 generates Turing patterns in a wide range of initial reagent concentrations at =15 and d 0.5 below the onset of percolation. 7,8 At lower 9 with similar or identical reactant concentrations, we observe a new pattern of black spots, which are actually islands of high concentration of the reduced catalyst Ru bpy 3 2+ in a continuous white background of oxidized Ru bpy Figure 1 demonstrates patterns at different initial concentrations of bromate. The left column shows patterns that emerge within 1 2 min after the first autocatalytic oxidation of Ru bpy 3 2+, which occurs after an induction period of min. The middle column demonstrates the relatively stationary patterns that develop from the transient patterns in the left column. White spot Turing patterns b are found at small NaBrO 3, while novel black spot patterns c f emerge at larger NaBrO 3. We also see in Fig. 1 d that black spots may be irregular in shape. Figure 1 f and its cloudy FFT show that these patterns need not have an intrinsic wavelength and that the spots may vary in size. To understand the origin of these black spots, we first measure the conductivity of the microemulsion to check that the BZ-AOT system is still below the percolation threshold at =9. This is indeed the case for d 0.45 Fig. 2 a. We then analyze the behavior of the BZ-AOT system in a CSTR FIG. 1. Turing patterns in the BZ-AOT system at different concentrations of NaBrO 3. Stationary patterns are in the middle column. Left column shows transient patterns that first emerge in the reactor. Parameters: =9, d =0.35, MA =0.25 M, H 2 SO 4 =0.32 M, Ru bpy 3 2+ =4.0 mm, NaBrO 3 /M= a,b 0.18, c, d 0.24, e, f 0.3. All concentrations are for the aqueous pseudophase not the total volume of microemulsion. Frame size= mm 2. Two-dimensional FFTs of the snapshots in the middle column are shown in the right column. Two images of the FFT correspond to two different scaling factors 1 and 2 of the FFT and allow one to see the FFT at different contrast levels. White dots correspond to FFT peaks in the allowable contrast range, 1 black 256 white. Figs. 2 b 2 d. After an induction period, which decreases with NaBrO 3, autocatalysis starts and then, depending on the bromate concentration, oscillations begin with a delay T ox Fig. 2 c. The larger NaBrO 3, the longer T ox, and the system stays for a longer time in the oxidized state Fig. 2 d. Comparing the behavior of the system in our spatially extended reactor and in a stirred tank reactor, we conclude FIG. 2. a Conductivity of an AOT microemulsion loaded with the reactants of the BZ reaction as a function of droplet fraction at two different. Concentrations: MA =0.25 M, H 2 SO 4 =0.32 M, NaBrO 3 =0.30 M. b, c Time series of bulk oscillations in a CSTR recorded as absorption at =650 nm arbitrary wavelength in the absorption band of the oxidized catalyst in a 1 cm square cell. Compositions of the BZ-AOT system as in Fig. 1; panels b and c correspond to a, b and e, f in Fig. 1, respectively. Period of oscillations 1 min. d Delay time T ox see c vs NaBrO 3. T ox is dependent on stirring rate and thus can vary slightly from experiment to experiment.

3 Black spots Belousov Zhabotinsky reaction J. Chem. Phys. 122, FIG. 3. Dependence of wavelength of Turing patterns on a MA at d =0.39 and b d at MA =0.27 M. Other parameters: =9, H 2 SO 4 =0.32 M, Ru bpy 3 2+ =4.0 mm, NaBrO 3 = a 0.27 M, b 0.24 M open triangles and 0.30 M closed triangles. FIG. 4. Evolution of a Turing pattern in the BZ-AOT system. a 5min corresponds to 36 min in e, b 45 min, c 95 min after emergence of the first pattern. Frame size= mm 2. d Time dependence of the wavelength of Turing patterns. e Time series of bulk oscillations in a CSTR recorded as absorption at =650 nm. Parameters: =9, d =0.39, NaBrO 3 =0.27 M, MA =0.10 M, H 2 SO 4 =0.32 M, Ru bpy 2+ 3 =4.0 mm. that black spots emerge at the time when the system first exhibits autocatalysis and stays in the oxidized state. Thus, the emergence of black spots at larger bromate concentrations is related to the fully oxidized state of the catalyst, or in other words to the oxidized state of the system. Increasing the sulfuric acid concentration from 0.24 M to 0.32 M at NaBrO 3 =0.18 M, MA =0.25 M and Ru bpy 3 2+ =3.2 mm gives the same tendency: white spot Turing patterns at smaller H 2 SO 4 transform to black spot Turing patterns at larger H 2 SO 4. A decrease in the catalyst concentration from Ru bpy 3 2+ =4 mm to 3.2 mm at NaBrO 3 =0.18 M, H 2 SO 4 =0.32 M, and MA =0.25 M brings the system from white spot to black spot Turing patterns. These dependences are in agreement with our model 10 and 11 described below and with the response of the model to the parameter m see Eq. 17 below. In all the above-mentioned experiments, emergence of black spots in the spatially extended reactor is accompanied by an increase in the time T ox in the CSTR. Decreasing the malonic acid concentration from 0.4 M to 0.1 M results in an increase of the Turing wavelength by a factor of 2 Fig. 3 a but no significant change in T ox. The wavelengths are obtained from patterns, such as those in the middle columns of Fig. 1, observed within 5 min after the first patterns emerge. At high MA M waves also develop; we discuss wavepattern interaction later. At smaller MA M, Turing patterns are extremely stable against waves: no waves emerge in the layer, though bulk oscillations occur in a CSTR see Fig. 4. The Turing patterns survive for about 3 h, evolving somewhat with time in our batch reactor, as neighboring black spots merge and the pattern s wavelength increases Fig. 4 d, for example, from 0.5 mm Fig. 4 a to 1mm Fig. 4 c after ca. 1.5 h. The stability of black spot Turing patterns at small MA may result from the fact that the steady state of the system is close to the fully oxidized state during the entire period of observation 3 h. As is seen from the kinetics in a CSTR Fig. 4 e, oscillations cease at the oxidized state of the system final absorption 0.2, while at larger MA 0.15 M oscillations conclude at the reduced state final absorption 0. We also studied how the emergence of black spots depends on droplet fraction. An increase in d from 0.35 to 0.5 at =9, NaBrO 3 = M, H 2 SO 4 =0.32 M, MA = M, and Ru bpy 3 2+ =4 mm results in a transition from black spots to white spots or labyrinthine Turing patterns. In all cases, the wavelength of the Turing patterns black or white spots decreases with d Fig. 3 b. At even larger d , above the onset of percolation, no stationary patterns could be found, and only waves were observed. The dependences of the wavelength of Turing patterns on MA Fig. 3 a and d Fig. 3 b and the time dependence shown in Fig. 4 d note that MA decreases with time suggest that one of the most important parameters affecting is the bulk concentration of malonic acid, which is equal to d MA. As Fig. 1 shows, a moderate increase in NaBrO 3 or in H 2 SO 4 leads to black spot patterns with a characteristic wavelength about 0.35 mm Fig. 1 d, but further increase in NaBrO 3 yields black spot patterns without a characteristic wavelength Fig. 1 f, as seen also from the blurred ring of the corresponding FFT. In general, patterns without an intrinsic wavelength can be considered as a set of localized stationary patterns. The latter can emerge via a subcritical Turing instability or bistability. In the case of a subcritical Turing instability, the system should be very sensitive to local perturbations, for example, due to small dust particles. Our previous experiments with dash waves 11 suggest that, in order to be bistable in a closed reactor, the BZ-AOT system should possess two different pools of water nanodroplets with significantly different radii. To choose between these alternatives, we performed experiments with dust particles and with quasielastic light scattering. Light scattering experiments revealed that an AOT microemulsion loaded with the BZ reactants NaBrO 3 H 2 SO M, MA =0.25 M at d 0.5 and =9

4 Kaminaga, Vanag, and Epstein J. Chem. Phys. 122, FIG. 5. Black spots in the BZ-AOT system at large NaBrO 3 for a, c specially washed glass windows and filtered microemulsion; b, d microemulsion without filtration and cleaning of window surfaces. Parameters: =9, d =0.35, MA =0.25 M, H 2 SO 4 =0.32 M, NaBrO 3 =0.33 M, Ru bpy 3 2+ =4.0 mm. Frame sizes are mm 2. Time lapse between snapshots a and c, as well as between b and d, is2min. has only a single narrow peak in the intensity of scattered light as a function of the droplet radius. Thus we have a monodisperse microemulsion, which suggests that the BZ- AOT system has only one steady state. Small dust particles are inevitably present in the BZ- AOT system. To eliminate these particles we first carried out a simple filtration of the microemulsion through a 0.2 m teflon filter, as we usually do in our light scattering experiments. We then carefully cleaned the surfaces of the glass windows. Some dust particles, probably precipitate formed in a reaction of the catalyst with bromine, 21 are firmly adsorbed on the windows and require removal with several solvents e.g., ethanol, acetone. Even such thorough cleaning procedures did not allow us to eliminate dust particles completely, but did significantly reduce their number. In Fig. 5 we show black spots in two microemulsions with the same composition at high concentrations of bromate with Fig. 5 a and without Fig. 5 b the above cleaning/ filtration procedures. After the induction period, black spots emerge only around the dust particles which are about 40 m in diameter. The microemulsion with fewer dust particles has a smaller number of black spots. Thus at larger concentrations of bromate or sulfuric acid, the BZ-AOT system becomes more sensitive to dust particles local large perturbations. Note that the black spots here are almost circular and are approximately equal in size 0.2 mm diameter. This observation suggests that the irregular shape of the black spots in Fig. 1 d is connected with irregular prepatterns Fig. 1 c that emerge immediately after the first autocatalytic transformation of the catalyst from the reduced to the oxidized state. B. Turing patterns in combination with waves and oscillations FIG. 6. Coexistence of black spots with a black reduction waves and b white oxidation waves at different times. c Space-time dependence at the gray line in snapshots a and b. Arrows in plot c indicate times at which snapshots a and b were taken. Time lapse between snapshots a and b is 7 min. Parameters: =9, d =0.35, NaBrO 3 =0.30 M, MA =0.25 M, H 2 SO 4 =0.32 M, Ru bpy 3 2+ =3.2 mm. Frame sizes: a, b mm 2, c 3mm 30 min. Wave velocity 50 m/s for black waves first waves in c, 17 m/s for white waves in b, and 12 m/s for last set of waves in c. Arrows in a and b indicate direction of wave propagation. Another interesting phenomenon found in this system is the interaction between black spots and waves or bulk oscillations presumably, Turing Hopf and/or Turing-excitability interactions. We observed two types of waves in our system: reduction or black waves at early stages in the life of black spot patterns Fig. 6 a and white waves waves of oxidation typically observed in the aqueous BZ reaction at later stages Fig. 6 b or, at larger concentrations of bromate Figs. 5 c and 5 d, shortly after the appearance of black spots. In the last case, reduction waves do not appear at all. The emergence of black waves in a thin layer and the start of bulk oscillations in a CSTR occur at approximately the same time; thus it is most likely that these waves arise from a Hopf instability rather than excitability. Since the oscillations start and end abruptly at large amplitude Figs. 2 c and 4 e, respectively, the transition between the oxidized steady state and the oscillatory regime appears to be subcritical. Black reduction waves were found in the aqueous BZ reaction at low MA / NaBrO 3 ratios, 22,23 but have not been observed previously in the BZ-AOT system. The reduction waves we find here have an ill-defined width and a steep trailing front. An unprecedented phenomenon found in this system is the propagation of black waves, as well as white waves, through Turing patterns Fig. 6. To see more clearly what happens in the system, we extract a space-time plot from a recorded movie using pixels along a line segment parallel to the direction of wave propagation and crossing several black spots. In a typical space-time plot Fig. 6 c, stationary black spots are represented by vertical black bands. Only stationary black spots are seen in the reactor for the first 5 min after their emergence. Then, black waves appear and propagate at about 50 m/s, significantly faster than white trigger waves in the BZ-AOT system, which have speeds of 2 10 m/s. 7,12 As seen in Fig. 6 c, black spots maintain their positions during the passage of the broad

5 Black spots Belousov Zhabotinsky reaction J. Chem. Phys. 122, FIG. 7. Stationary a, b and oscillatory c, d Turing patterns in the BZ-AOT system at a 5 min, b 25 min, c 35 min, d 35.5 min after the emergence of the pattern. e Space-time dependence at the gray line in snapshot b. Arrows in plot e correspond to snapshots a d. Parameters: =9, d =0.39, NaBrO 3 =0.27 M, MA =0.20 M, H 2 SO 4 =0.32 M, Ru bpy 2+ 3 =4.0 mm. Frame sizes: a d mm 2, e 3mm 40 min. black waves despite the fact that we temporarily lose information about the location of the spots while the waves are passing through them. Figure 6 c also shows that black spots do not oscillate during the passage of black waves. The black spots are essentially within the accuracy of our experiment stationary during the propagation of black waves. Black waves can be observed only for a few minutes, when the system is very close to the fully oxidized state. Then, black waves gradually transform to white waves, and after ca. 5 7 min we clearly see white waves propagating through the still stationary black spots Fig. 6 b. The velocity of these white waves is about 17 m/s, which is still higher than that of trigger waves. The propagation of waves through stationary black spots resembles the behavior of trigger waves in the presence of small obstacles. If the linear size of the obstacle is small, as in our case, the wave proceeds unchanged. Large obstacles cause waves to break, giving rise to spirals. The difference is that our medium is homogeneous and our obstacles are islands of reduced catalyst. Theoretical questions about how and under what circumstances waves are able to propagate through Turing structures will be discussed below. After the appearance of white waves, the black spots slowly fade out as the vertical black stripes in Fig. 6 c vanish. Finally bottom of Fig. 6 c we observe only thin white spiral waves with velocities around 12 m/s. At larger bromate concentrations, we find another scenario for the emergence of white waves. As shown in Figs. 5 c and 5 d, after a stationary period black spots start to emit white circular waves. Figure 5 clearly demonstrates that the simultaneous presence of black spots and white waves is possible if the stationary steady state of the system gray background lies between the fully reduced and fully oxidized states, so that the catalyst in the black spots is between the fully reduced state black and the background gray, while in the areas occupied by white waves, the catalyst concentration is between the background and the fully oxidized white levels. Figure 2 c suggests that the steady state indeed lies between the fully reduced and fully oxidized states, since the concentration during the period T ox corresponds to the pseudosteady state. Note also that black spots remain clearly visible when white waves pass through them. In areas between two waves on the gray background black spots are seen only at the early stages of white waves for several minutes after the waves emergence. At later stages, black spots become nearly invisible on the gray almost black background, as in Fig. 6 b. We interpret this behavior as arising from the continuously changing steady state level, which approaches the reduced steady state of the catalyst see also Fig. 2 c. As shown in Fig. 1 a, black spots can be transformed into white spots we use the term spot here for any small isolated area regardless of its shape round, elongated, or even short stripes, by decreasing NaBrO 3. This transition is also favored by increases in Ru bpy 3, as well as in d. Making small corresponding changes in the concentrations of these reactants and d and holding H 2 SO 4 constant and MA above 0.15 M, we obtained first black spots labyrinthine on a white background, Fig. 7 a that quickly transformed to a white labyrinth with continuous black background Fig. 7 b. After about 30 min, the white spots begin to oscillate Figs. 7 c 7 e with a period T 1 min. The period of the corresponding CSTR oscillations is 1.3 min. Two antiphase snapshots separated in time by T/2 are shown in Figs. 7 c and 7 d. In the corresponding movie, these oscillations look like very fast phase or kinematic waves moving chaotically in different directions, since oscillations in different parts of the reactor are not synchronized. But, unlike the continuous waves in Fig. 6 that pass through or between black spots, the very fast phase waves in Fig. 7 are not continuous. The black background is stationary within the accuracy of our CCD camera and never becomes white. Therefore this behavior consists not of waves but of oscillatory Turing patterns. In experiments with another oil tetradecane at almost the same reactant concentrations as in Fig. 7 NaBrO 3 =0.24 M, Ru bpy 3 2+ =3.2 mm, =10, d =0.28, we observed Turing patterns that oscillated synchronously without phase waves in a large area, about 6 5 mm 2. Comparison of Figs. 6 and 7, particularly the time-space plots 6 c and 7 e, and the fact that concentrations of reactants and d are quite similar in the two experiments suggests that the oscillatory white-spot Turing patterns and the waves propagating through black-spot Turing patterns are related phenomena. IV. THEORETICAL CONSIDERATIONS In this section we present a qualitative description of the observed patterns on the basis of a simple model. The BZ-

6 Kaminaga, Vanag, and Epstein J. Chem. Phys. 122, AOT system is so complex that no single model captures all the patterns found to date. The Field Kőrös Noyes model 24 describes the chemistry of the BZ system quite well, but contains a large number of variables. The Oregonator model 7,25 for the BZ reaction is simple, but does not explicitly take into account the catalyst concentration and cannot explain phenomena occurring at the oxidized state of the catalyst, like reduction waves. 23 This is a severe limitation for studying the phenomena of interest here, since black spots emerge at large concentrations of bromate and/or sulfuric acid, when the system is in the fully oxidized state. A model that nicely describes oscillations close to the oxidized state 26 has been applied successfully to explain collision-stable reduction waves, 27 which we also observed in our experiments. This many-variable model, however, is still quite complex and time consuming for simulations of the spatially extended BZ-AOT system. We have therefore developed a new Oregonator-like model with just one additional phenomenological term that shuts down the autocatalytic reaction when the concentration of oxidized catalyst approaches its total concentration. This model is based on the following chemical equations: 2H + + A + Y X k 1, 1 H + + X + Y 0 k 2, 2 H + + A + X + 2red 2X +2Z k 3, 3 2X 0 k 4, 4 B + Z hy + red k 5, 5 where A=BrO 3, B is a mixture of malonic and bromomalonic acids, X=HBrO 2, Y =Br,red is the reduced state of the catalyst, and Z is its oxidized state red + Z C 0, h is a stoichiometric coefficient, k 3 and 3 the rate of reaction 3 depend on red as k 3 = k 3 red / red + c, 6 3 = k 3 X A H + red/ red + c, where c C 0, the total catalyst concentration. If red c, then k 3 k 3 and 3 =k 3 X A H + =k 3 X, as in the classical Oregonator model H + and A=BrO 3 are present in such high initial concentrations that they remain essentially unchanged during the reaction. If Z is close to C 0 and consequently red c, then k 3 tends to zero. Equations 6 and 7 reflect the competition between the forward reaction BrO * 2 + red HBrO 2 + Z 8 and the back reaction of radical recombination BrO * 2 + BrO * 2 BrO 3 + HBrO 2. At large red, reaction 8 dominates and k 3 k 3. At small red, however, reaction 9 prevails and the autocatalytic reaction 3 stops. On introducing dimensionless variables X =k 3 x/ 2k 4, Y =k 3 y/k 2, Z =k 3 2 z/ k 4 k 5, t= /k 5, and parameters f =2h, q=2k 4 k 1 / k 2 k 3, m=k 3 2 / C 0 k 4 k 5, 1 =c/c 0, =k 5 /k 3, where 7 9 FIG. 8. a Oscillations in 0D system 10 and 11 with f =1.2, q=0.001, m=182, 1 =0.01, and =0.01 curve 1, curve 2 ; max z =1/m = The same shape of oscillations as in curve 2 is obtained at m=181 and =0.01. b Curves 1 and 2 are nullclines for Eqs. 10 and 11, respectively, with parameters: f =1.2, q=0.001, m=182, 1 =0.01, and = curve 3, trajectory of large amplitude oscillations, 0.01 curve 4, trajectory of small amplitude oscillations. Steady state cross point of curves 1 and 2 is inside the loops 3 and 4. c Phase diagram of system 10 and 11 in logarithmic coordinates. Parameters f =0.98, q=0.001, 1 =0.01, D z /D x =2.75. SS is steady state with fully oxidized catalyst, T+H means Turing+Hopf instabilities. Curve 1 white squares corresponds to the onset of Turing instability, curve 2 to the onset of Hopf instability. Points of curve 2 are fitted by a linear trend line with slope d Dispersion curves for system 10 and 11 with f =1.2, q=0.001, m=182, 1 =0.01, =0.01, D X =1. Curves 1 correspond to Fig. 10 a D Z =3.2, and curves 2 to Fig. 10 b D Z =2.72. Curves 1r and 2r are real part of eigenvalue, curves 1i and 2i are corresponding imaginary parts. Two horizontal dashes indicate half of the imaginary part at k=0. Two vertical dashes indicate half of k T Re has a maximum at k=k T. k 1 =k 1 H + 2 A, k 2 =k 2 H +, k 3 =k 3 H + A, k 5 =k 5 B, the rate equations for 1 5 with diffusion terms take the form x/ = 1/ q x fz/ q + x + x 1 mz / 1 +1 mz x 2 + D x x, 10 z/ = x 1 mz / 1 +1 mz z + D z z, 11 where the quasistationary state approximation dy/d =0 has been made. The diffusion coefficients D x and D z are also dimensionless, and only their ratio is important. The Laplacian 2 / r 2 in 1D and 2 / x / y 2 in 2D, where r, x, and y here represent spatial coordinates. Equations 10 and 11 have only a single steady state. They correctly reproduce the oscillatory behavior near the fully oxidized state of the catalyst Fig. 8 a found in previous experiments. 26 For modeling the BZ-AOT system, Eqs. 10 and 11 can be extended by introducing a fast-diffusing activator s, BrO * 2 radical and a fast-diffusing inhibitor u,br 2 in the oil phase: 7,10 x/ = fz q x / q + x + x 1 mz / 1 +1 mz x 2 x + s / + D x x, 12

7 Black spots Belousov Zhabotinsky reaction J. Chem. Phys. 122, z/ = x 1 mz / 1 +1 mz z z + u + D z z, 13 s/ = x s / 2 + D s s, u/ = z u / 3 + D u u This extension increases flexibility and allows us to simulate black spots with D z =D x. However, if we take D z D x, the two-variable model 10 and 11 is sufficient to obtain black spots. Since these two equations are simpler than the fourvariable model 12 15, we employ them. The system 10 and 11 can have both supercritical and subcritical Hopf bifurcations between the fully oxidized steady state and the oscillatory state, can exhibit excitability, and can show canard behavior see Figs. 8 a and 8 b, in which a small parameter change near a supercritical Hopf instability transforms small amplitude oscillations to large amplitude oscillations. The combination of these bulk temporal 0D behaviors with supercritical and subcritical Turing instabilities gives rise to a variety of spatiotemporal patterns: stationary Turing patterns black spots, localized stationary Turing patterns, several types of oscillatory Turing patterns, including subharmonic Turing and subharmonic Hopf modes, and reduction waves. Two parameters of model 10 and 11, namely, and m, depend on the concentrations of bromate, sulfuric acid, and malonic acid, which allows us to make a comparison with the experimental data; = k 5 /k 3 B / H + A, m = k 3 2 / k 4 k 5 H + 2 A 2 / C 0 B The parameter f depends on the ratio between the bromomalonic and malonic acid concentrations and should increase with time in our batch experiment, since initially there is no bromomalonic acid, so f =0. Thus we assume that f is small when we are looking for black spots in the model. As a first step in investigating the dynamical behavior of system 10 and 11, we carried out a linear stability analysis and constructed dispersion curves, like those shown in Fig. 8 d. An eigenvalue of the Jacobian matrix of Eqs. 10 and 11 with Re =0, Im 0 atk=0 implies the supercritical onset of a Hopf instability. If Re =0 and Im =0 at k 0, we have the onset of a Turing instability also supercritical. Typical dynamical behavior of system 10 and 11 in the -m parameter plane at small f is shown in Fig. 8 c. Depending on D z /D x, the Turing line onset of Turing instability can lie above or below the Hopf line, but the two lines are always nearly parallel. The slope of these lines in logarithmic coordinates is about 0.3. If we increase, e.g., bromate A, as in the experiments shown in Fig. 1, decreases as A 1 and m increases as A 2, in accord with Eqs. 16 and 17. Since the slope of ln A 2 vs ln A 1 is 2, i.e., significantly less than 0.3, an increase in A pushes the system toward the steady state domain, as we find in experiment. Thus this simple test of model 10 and 11 is consistent with our experimental results. FIG. 9. Turing patterns in model 10 and 11. a Parameters: f =1.2, q =0.001, m=190, 1 =0.01, =0.01, D X =1, D Z =3.25, size=5 T 5 T, T =1.76. At these parameters, Re,k=0 2.1 and Re,k=k T 0.29, thus the system exhibits pure Turing instability; T =2 /k T. b An increase in Re,k=k T by increasing D Z to 100, transforms black spots in a into stripes. c Localized Turing spot. Parameters: f =0.8, q=0.001, m=23, 1 =0.01, =0.01, D X =1, D Z =10; size=8 T 8 T, T =1.92. Re,k= , Re,k=k T All eigenvalues have negative real part at all wave numbers k. Stationary single spot emerges only after a large local perturbation of the steady state, x ini =0.1 x SS at the center of the square area, where x SS is the steady state value of x. d Increasing Re,k=k T to 3 by decreasing m to 22 leads to the slow spreading of Turing spots, which finally occupy the entire area; all other parameters as in c. Examples of hexagonal black spots are shown in Figs. 9 a and 9 d. Atf 1.2 Figs. 9 a and 9 b, the Turing instability is supercritical. If Re decreases below zero at k=k T k T is the wave number at which the Turing peak reaches its maximum, i.e., at which Re is the most positive, k T 0 because of, e.g., an increase in m, we observe a slow disappearance of Turing patterns. Increasing a positive Re, e.g., by increasing D z or decreasing m, leads to transformation of black spots into stripes Fig. 9 b. At smaller f, around 0.8, the Turing instability is subcritical, and Turing patterns that emerge at positive Re survive at negative Re, or can be induced locally from the homogeneous steady state by a sufficiently large initial perturbation Fig. 9 c. These localized Turing patterns are analogous to those found theoretically by Jensen et al. 27,28 in the Lengyel Epstein model 29 of the CIMA reaction or those found in a model of the BZ-AOT system when the steady state is close to the fully reduced state of the catalyst. 8 Note that in the latter case the localized patterns are white spots. If we decrease m equivalent to a decrease in BrO 3 in the experiment, black spots start to occupy the entire area Fig. 9 d. The Turing wavelength T T =2 /k T calculated with model 10 and 11 increases slightly with MA, in contradiction to the experimental data shown in Fig. 3 a. We attribute this discrepancy to the fact that the model neglects the

8 Kaminaga, Vanag, and Epstein J. Chem. Phys. 122, FIG. 10. Space-time plot for oscillatory Turing patterns in 1D. Model parameters: f =1.2, q=0.001, m=182, 1 =0.01, =0.01, D x =1, space=6 T, T = a 1.78, b 1.7; D z = a 3.2, b White corresponds to maximum of activator x and black to minimum. Total time= a 0.6, b reaction of MA with Br 2.IfBr 2 is the fast-diffusing inhibitor, then its diffusion length l d, the average distance a molecule of Br 2 diffuses before undergoing a chemical reaction, should be strongly dependent on MA and d. The Turing wavelength T should increase with l d. The latter quantity would decrease with MA if the MA-Br 2 reaction were taken into account. Model 10 and 11 possesses a solution for oscillatory Turing patterns due to Turing Hopf interaction and can therefore mimic the interaction between black or white Turing spots and waves or oscillations. If the Turing peak in the dispersion curve like those shown in Fig. 8 d is too large, only stationary Turing patterns are found. Oscillatory Turing patterns emerge at the smaller Re shown in Fig. 8 d. Figure 10 shows two different types of oscillatory Turing patterns in 1D found at system parameters corresponding to the dispersion curves in Fig. 8 d. Figure 10 a demonstrates patterns in which all peaks oscillate synchronously and are fixed in space. In accord with the terminology introduced by DeWit et al., 30 these patterns are referred to as mixed mode. In-phase oscillations of all Turing peaks start supercritically as we gradually decrease Re at k=k T in the dispersion curves e.g., by decreasing D z Ref. 31. Figure 10 b demonstrates patterns in which odd and even peaks oscillate antiphase, a subharmonic mode. 30 This mode emerges at still smaller Re, as shown in Fig. 8 d, curve 2r. The transition of the Turing peaks from in-phase oscillations to out-of-phase oscillations occurs through a behavior that we characterize as chaoticlike, since chaos is not welldefined in spatially extended systems. Here, chaoticlike means that in a range of the parameter D z , all other parameters as in Fig. 10 the pattern oscillates aperiodically in time with no long-range spatial correlation. We measured long time sequences at many spatial points. As shown in Fig. 11, for D z =3.2 and D z =2.72, the FFTs of these sequences have well-defined peaks at multiples of the fundamental frequencies and Hz, respectively if the dimensionless time is measured in seconds at all spatial points examined. For D z =2.80, 2.85, 2.9, and 3.0, the FFT spectra appear noisy, lose their simple harmonic structure, and differ at different spatial points. More significantly, we calculated correlation functions for time sequences measured at pairs of spatial points separated by a distance d. For D z =3.2 and D z =2.72, these correlation functions are periodic with period T and 2 T, respectively. For D z =2.8 3, the FIG. 11. FFT spectra for time series of 1D oscillatory Turing patterns shown in Fig. 10 for D z = a 3.2, b 2.9, c 2.8, d Bold solid and thin dotted lines in b and c correspond to different spatial points. correlation functions are aperiodic, and can be approximated by the function cos 2 d/ T exp d/l, where the correlation length L characterizes the decay of the correlation. The smaller the correlation length, the more chaotic is the behavior. The smallest L values, 1 2 T, are found near D z =2.85. A theoretical interpretation of the mixed and subharmonic modes remains elusive. The subharmonic mode is expected to emerge as the result of a resonance between a subharmonic Hopf mode at k=k 1/2 defined as Im,k =k 1/2 = 1/2 Im,k=0, see horizontal dashes in Fig. 8 d and a Turing mode at k T. Thus we should have 2k 1/2 k T. For the subharmonic mode, however, 2k 1/2 is significantly larger than k T, and Re at k=2k 1/2 is negative. On the other hand, for the mixed mode, the relation 2k 1/2 k T is approximately true. Note also that in both cases k T /2 is in the range of wave numbers in which Im 0 see vertical dashes in Fig. 8 d. It is worth noting also that out-of-phase oscillations of Turing peaks can occur even when Re 0 atk=0 but Re 0 atk=k T due to either excitability or subcriticality of the Hopf instability at some f. Another explanation of mixed and subharmonic modes 32 involves Floquet multipliers that can be calculated for the homogeneous limit cycle as a function of the wave number k. Subharmonic modes can give rise to a variety of patterns in 2D. Twinkling hexagons, first found theoretically as a resonant interaction between wave and Turing instabilities, 33 and oscillatory labyrinthine patterns are shown in Fig. 12. At small f 1.2, there is a large area in the parameter space m, f, ; q=0.001, 1 =0.01, where the Hopf bifurcation is subcritical and the steady state is close to the fully oxidized state. In this region, model 10 and 11 produces reduction waves when the catalyst diffusion coefficient D Z is smaller than D X and significantly smaller than the critical value at which Turing instability emerges. Figure 13 a demonstrates such waves. They are characterized by broad leading spatial coordinate and steep trailing spatial coordinate fronts. The shape of a newly emerged reduction wave spatial coordinate changes with

9 Black spots Belousov Zhabotinsky reaction J. Chem. Phys. 122, FIG. 14. White Turing spots at large f. f =2.18, q=0.01, m=10, 1 =0.01, =0.1, D x =1, D z =100. a, b space=16 T, T =11.3, 1D case, z SS = For c g, t=5, 17, 20, 45, and 300, respectively, size=8 T 2 T, h t=21, size=11 T 11 T. In all cases, vertical thin line in the center is perturbed initially. FIG. 12. Oscillatory Turing patterns in 2D for model 10 and 11. Parameters: f =1.2, q=0.001, m=182, 1 =0.01, =0.01, D x =1. For a l, D z =2.8 and t=0.054 between snapshots ; one full period 0.6, size=5 T 6.9 T, T =1.76. For m, n, D z =2.68 and t=0.345 half a period, size=5 T 5 T, T =1.68. time, since the trailing front propagates faster than the leading front until the wave reaches a stationary shape space coordinate 70 90, which depends on m, f,, and D Z /D X. The shape of the corresponding temporal oscillations for the 0D system can be seen in Fig. 13 b. Comparing Figs. 13 b and 13 a, we observe that the time profile for the oscillations and the spatial profile of the newly emerged wave are identical. Simulations show that the reduction waves require D Z D X, while Turing patterns can be found in model 10 and 11 only at D Z D X. Reduction waves found earlier in other models of the BZ reaction also require the relation FIG. 13. a One-dimensional Reduction waves in model 10 and 11. Bold line is catalyst z; thin dotted line is activator x. Arrow shows direction of wave propagation. At initial moment of time, left edge of the segment with zero flux boundary conditions was perturbed. Parameters: f =0.55, q=0.001, m=2, 1 =0.01, =0.072, D X =1, D Z =0 similar waves are obtained at D Z =0.2. The system is close to the subcritical Hopf bifurcation. At these parameters, Re = at k=0, while at =0.071, Re = atk=0. b Homogeneous large-amplitude oscillations on a segment of length 80 after homogeneous large-amplitude perturbation x 0 =x SS ; x SS = , z SS = , D Z =30. Oscillations are the same as in 0D case. D Z D X. 27 Thus, a more complex model, perhaps model 12 15, will be required to explain propagation of waves reduction and/or oxidation through Turing patterns. The last phenomenon we consider is the interaction of Turing instability with excitability or subcritical Hopf instability. Turing patterns can develop from the steady state close to the fully reduced state white spots at large f or from the fully oxidized state black spots at small f. In model 10 and 11, at large f and consequently at small product z SS m, the system is nearly independent of m and behaves like the classical Oregonator model. Due to excitability or subcritical Hopf instability of the system in 0D, Turing instability can lead to unusual patterns with several characteristic wavelengths in the BZ-AOT system Fig. 14. The Turing instability causes a small initial deviation from the steady state to grow. If this deviation crosses a threshold, excitability, rather than the Turing mode, amplifies this deviation, leading to a large, broad peak of catalyst that serves as an inhibitor in the two-variable Oregonator model. Catalyst concentrations exceed the steady state z SS = at all spatial points. In normal Turing patterns, the maxima and minima of all variables lie above and below the steady state levels, respectively. The broadly distributed inhibitor in this system suppresses patterns in the neighborhood of the initial disturbance, and finally we observe in 1D stationary patterns with a very large wavelength Figs. 14 a and 14 b, about three Turing wavelengths. We stress that, in analyzing dynamical patterns, it is important not only to determine the type of instability Turing, Hopf or finite wavelength Hopf that destabilizes the homogeneous steady state, but also to understand what processes stabilize the resulting patterns. In 2D, we observe patterns with two different wavelengths Figs. 14 e 14 h. If we perturb an initial narrow vertical stripe, Turing spots first emerge with the characteristic Turing wavelength Fig. 14 c. Then, these spots split into two Fig. 14 d before distant spots emerge Fig. 14 e with the same wavelength as in 1D. Finally, after a long period of time, a stationary structure is established Fig. 14 g. In larger systems Fig. 14 h, we find patterns with

10 Kaminaga, Vanag, and Epstein J. Chem. Phys. 122, V. DISCUSSION AND CONCLUSION FIG. 15. Stationary Turing patterns in model 10 and 11 in the oxidized steady state of the system. Parameters: f =0.55, q=0.001, m=2, 1 =0.01, =0.072, D X =1 as in Fig. 8, D Z =30. a c small random initial perturbation x 0 =x SS RANDOM 1 ; x SS = , z SS = ; a, b x light line and z dark line profiles at t=75 and t=370, respectively; c z profiles at t= 1 75, 2 79, 3 90, 4 135, 5 375, z is shifted down by 1 0, , , 4 0.1, and Small amplitude oscillations at each spatial point with a period of 2.14 gradually die and at t=50 75 first Turing patterns shown in a start to emerge. d Large initial perturbation at one edge of the segment. Horizontal lines in d indicate steady states. Wavelength = a 8.64, b 12.3; d significantly different distances between neighboring spots. This last case is extremely time consuming to simulate, and we were able to calculate only a relatively short period of time up to t=21. Similar behavior can occur with Turing patterns in a BZ- AOT system close to the oxidized state and possessing a subcritical Hopf instability. Figure 15 shows Turing patterns black spots in system 10 and 11 when the steady state is close to the fully oxidized state. When we apply a smallamplitude random initial perturbation to the steady state, Turing patterns with the theoretically predicted wavelength T =8.64 start to grow Fig. 15 a. Some peaks grow faster than others Fig. 15 c, curve 2 and significantly affect the wavelength of neighboring peaks Fig. 15 c, curve 3. All Turing peaks emerge in a relatively short time t 15 between curves 1 and 3, and we can observe Turing patterns with different wavelengths for an extended interval t 200. The peaks of these patterns slowly drift and rearrange, ultimately forming stationary Turing patterns Fig. 15 b with a wavelength of about If we make a large perturbation of the same system at one edge, Turing patterns with wavelength 14.5 emerge Fig. 15 d. All values of z in these patterns are below the steady state as in Fig. 14 b, except that there z lies above the steady state. We see that the wavelength of the final patterns depends on the initial perturbation and can deviate significantly from the theoretically predicted Turing wavelength T, depending on the parameters of the system. This theoretical finding suggests an explanation for the fact that the Turing patterns in some of our experiments do not form regular hexagons and that their wavelength can vary from spot to spot, giving a cloudy FFT. We have found a series of stationary patterns in surfactant-rich microemulsions with small nanodroplet radii at 9. At low bromate or sulfuric acids concentrations, the usual white spots appear. At higher concentrations, a new phenomenon, black spots, emerges. At still larger concentrations of bromate these black spots become irregular, without an intrinsic wavelength. In general, such patterns can arise in several ways. Three possible mechanisms are 1 subcritical Turing instability with negative Re at k=k T, 2 Turing instability with a broad Turing peak spanning a large enough range of wave numbers to induce interacting subharmonic modes, and 3 bistability between steady states. Usually, the BZ reaction has only one steady state, but in freshly prepared microemulsions, two pools of nanodroplets with different radii can emerge. 11 Interaction between these pools can lead to bistability. Our light scattering experiments reveal, however, that for the small used here, even freshly prepared microemulsions have a monodisperse distribution of droplet radii at small droplet fraction below the percolation threshold. Thus we do not have evidence for bistability. Choosing among mechanisms 1 and 2, we found strong evidence in support of option 1. First, our experimental results reveal that the BZ-AOT system at large bromate concentrations irregular black spots is very sensitive to dust particles. Special filtration and cleaning the surface of glass windows significantly decreases the number of black spots. Theoretical analysis with model 10 and 11 shows that localized Turing spots due to subcritical Turing instability are possible. In this case, spots can emerge anywhere and do not necessarily have a characteristic wavelength. Localized black spots are analogous to the white localized spots found recently at the opposite side of the oscillatory domain, close to the boundary with the reduced catalyst steady state. 8 We have found that malonic acid strongly affects the Turing wavelength. Our experiments suggest that MA 1/2. This dependence can be attributed to the diffusion wavelength of Br 2, l d D/k 1/2, where k is the characteristic reaction rate constant for disappearance of Br 2, for example, by reaction with MA; thus k=k MA k =29 M 1 s 1 in the aqueous phase 34. We have also found oscillatory Turing patterns and black spot Turing patterns through which black reduction and white oxidation waves can propagate. Though these patterns resemble one another in some respects, oscillatory Turing patterns can be explained by the simple two-variable model 10 and 11, while waves propagating through Turing patterns require a more complex model. We believe that model 10 and 11 is a good building block for a more detailed model of the BZ-AOT system. Although Turing Hopf interaction more precisely, interaction between instabilities leading to spatial or temporal symmetry breaking is well documented in physical systems e.g., fluid convection, there are only a few experiments in chemical systems where this interaction has been studied. 39,40 All of these experiments were performed on the CIMA reaction in a CFUR, in which a gradient of concentrations in a direction perpendicular to the gel layer is inevi-

Front velocity in models with quadratic autocatalysis

Front velocity in models with quadratic autocatalysis JOURNAL OF CHEMICAL PHYSICS VOLUME 117, NUMBER 18 8 NOVEMBER 2002 Front velocity in models with quadratic autocatalysis Vladimir K. Vanag a) and Irving R. Epstein Department of Chemistry and Volen Center

More information

Excitation of Waves in a Belousov-Zhabotinsky System in Emulsion Media

Excitation of Waves in a Belousov-Zhabotinsky System in Emulsion Media Excitation of Waves in a Belousov-Zhabotinsky System in Emulsion Media Bull. Korean Chem. Soc. 2008, Vol. 29, No. 11 2241 Excitation of Waves in a Belousov-Zhabotinsky System in Emulsion Media C. Basavaraja,

More information

Temperature control of pattern formation in the Ru(bpy) catalyzed BZ-AOT system

Temperature control of pattern formation in the Ru(bpy) catalyzed BZ-AOT system PAPER www.rsc.org/pccp Physical Chemistry Chemical Physics Temperature control of pattern formation in the Ru(bpy) 3 2+ -catalyzed BZ-AOT system Rachel McIlwaine, Vladimir K. Vanagw and Irving R. Epstein*

More information

Time-periodic forcing of Turing patterns in the Brusselator model

Time-periodic forcing of Turing patterns in the Brusselator model Time-periodic forcing of Turing patterns in the Brusselator model B. Peña and C. Pérez García Instituto de Física. Universidad de Navarra, Irunlarrea, 1. 31008-Pamplona, Spain Abstract Experiments on temporal

More information

7. Well-Stirred Reactors IV

7. Well-Stirred Reactors IV 7. Well-Stirred Reactors IV The Belousov-Zhabotinsky Reaction: Models and Experiments Oregonator [based on the FKN mechanism; Field, R. J. & Noyes, R. M.: Oscillations in chemical systems. IV. Limit cycle

More information

Resonant suppression of Turing patterns by periodic illumination

Resonant suppression of Turing patterns by periodic illumination PHYSICAL REVIEW E, VOLUME 63, 026101 Resonant suppression of Turing patterns by periodic illumination Milos Dolnik,* Anatol M. Zhabotinsky, and Irving R. Epstein Department of Chemistry and Volen Center

More information

Dynamic Mechanism of Photochemical Induction of Turing Superlattices in the Chlorine Dioxide-Iodine-Malonic Acid Reaction-Diffusion System

Dynamic Mechanism of Photochemical Induction of Turing Superlattices in the Chlorine Dioxide-Iodine-Malonic Acid Reaction-Diffusion System 5382 J. Phys. Chem. A 2005, 109, 5382-5387 Dynamic Mechanism of Photochemical Induction of Turing Superlattices in the Chlorine Dioxide-Iodine-Malonic Acid Reaction-Diffusion System Igal Berenstein, Lingfa

More information

Diffusion-Induced Instabilities near a Canard

Diffusion-Induced Instabilities near a Canard J. Phys. Chem. 1995,99, 1593-1511 1593 Diffusion-nduced nstabilities near a Canard Frank Buchholtz" nstitute of Physical Chemistry, University of Wurzburg, Marcusstr. 9-11, 977 Wurzburg, Germany Milos

More information

arxiv:chao-dyn/ v1 12 Feb 1996

arxiv:chao-dyn/ v1 12 Feb 1996 Spiral Waves in Chaotic Systems Andrei Goryachev and Raymond Kapral Chemical Physics Theory Group, Department of Chemistry, University of Toronto, Toronto, ON M5S 1A1, Canada arxiv:chao-dyn/96014v1 12

More information

Diffusive instabilities in heterogeneous systems

Diffusive instabilities in heterogeneous systems JOURNAL OF CHEMICAL PHYSICS VOLUME 119, NUMBER 14 8 OCTOBER 2003 Diffusive instabilities in heterogeneous systems Vladimir K. Vanag and Irving R. Epstein Department of Chemistry and Volen Center for Complex

More information

Resonance in periodically inhibited reaction diffusion systems

Resonance in periodically inhibited reaction diffusion systems Physica D 168 169 (2002) 1 9 Resonance in periodically inhibited reaction diffusion systems Karl Martinez, Anna L. Lin, Reza Kharrazian, Xaver Sailer, Harry L. Swinney Center for Nonlinear Dynamics and

More information

Pattern Formation and Wave Propagation in the Belousov-Zhabotinskii Reaction

Pattern Formation and Wave Propagation in the Belousov-Zhabotinskii Reaction Pattern Formation and Wave Propagation in the Belousov-Zhabotinskii Reaction Dan Hill and Timothy Morgan University of California, San Diego Physics Department INTRODUCTION The Belousov-Zhabotinskii (BZ)

More information

Complex Patterns in a Simple System

Complex Patterns in a Simple System Complex Patterns in a Simple System arxiv:patt-sol/9304003v1 17 Apr 1993 John E. Pearson Center for Nonlinear Studies Los Alamos National Laboratory February 4, 2008 Abstract Numerical simulations of a

More information

Patterns in the Belousov Zhabotinsky reaction in water-in-oil microemulsion induced by a temperature gradient

Patterns in the Belousov Zhabotinsky reaction in water-in-oil microemulsion induced by a temperature gradient PAPER www.rsc.org/pccp Physical Chemistry Chemical Physics Patterns in the Belousov Zhabotinsky reaction in water-in-oil microemulsion induced by a temperature gradient Jorge Carballido-Landeira, ab Vladimir

More information

Semi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical

Semi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical ANZIAM J. 53 (EMAC211) pp.c511 C524, 212 C511 Semi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical M. R. Alharthi 1 T. R. Marchant 2 M. I.

More information

Spatial Periodic Perturbation of Turing Pattern Development Using a Striped Mask

Spatial Periodic Perturbation of Turing Pattern Development Using a Striped Mask 4428 J. Phys. Chem. A 2003, 107, 4428-4435 Spatial Periodic Perturbation of Turing Pattern Development Using a Striped Mask Igal Berenstein, Milos Dolnik,* Anatol M. Zhabotinsky, and Irving R. Epstein

More information

Theoretical Models for Chemical Kinetics

Theoretical Models for Chemical Kinetics Theoretical Models for Chemical Kinetics Thus far we have calculated rate laws, rate constants, reaction orders, etc. based on observations of macroscopic properties, but what is happening at the molecular

More information

CHEM 515: Chemical Kinetics and Dynamics

CHEM 515: Chemical Kinetics and Dynamics Alejandro J. Garza S01163018 Department of Chemistry, Rice University, Houston, TX email: ajg7@rice.edu, ext. 2657 Submitted December 12, 2011 Abstract Spontaneous antispiral wave formation was observed

More information

Chemical Wave Packet Propagation, Reflection, and Spreading

Chemical Wave Packet Propagation, Reflection, and Spreading 11676 J. Phys. Chem. A 2002, 106, 11676-11682 Chemical Wave Packet Propagation, Reflection, and Spreading Lingfa Yang and Irving R. Epstein* Department of Chemistry and Center for Complex Systems, MS 015,

More information

Nature-inspired Analog Computing on Silicon

Nature-inspired Analog Computing on Silicon Nature-inspired Analog Computing on Silicon Tetsuya ASAI and Yoshihito AMEMIYA Division of Electronics and Information Engineering Hokkaido University Abstract We propose CMOS analog circuits that emulate

More information

Georgia Institute of Technology. Nonlinear Dynamics & Chaos Physics 4267/6268. Faraday Waves. One Dimensional Study

Georgia Institute of Technology. Nonlinear Dynamics & Chaos Physics 4267/6268. Faraday Waves. One Dimensional Study Georgia Institute of Technology Nonlinear Dynamics & Chaos Physics 4267/6268 Faraday Waves One Dimensional Study Juan Orphee, Paul Cardenas, Michael Lane, Dec 8, 2011 Presentation Outline 1) Introduction

More information

Chaos in the Belousov-Zhabotinsky Reaction

Chaos in the Belousov-Zhabotinsky Reaction Chaos in the Belousov-Zhabotinsky Reaction David Connolly Brad Nelson December 2, 2011 Abstract In this project, we investigate two different models of the Belousov-Zhabotinsky Reaction, the Oregonator

More information

Competition of Spatial and Temporal Instabilities under Time Delay near Codimension-Two Turing Hopf Bifurcations

Competition of Spatial and Temporal Instabilities under Time Delay near Codimension-Two Turing Hopf Bifurcations Commun. Theor. Phys. 56 (2011) 339 344 Vol. 56, No. 2, August 15, 2011 Competition of Spatial and Temporal Instabilities under Time Delay near Codimension-Two Turing Hopf Bifurcations WANG Hui-Juan ( )

More information

Flow- and Diffusion Distributed Structures with noise at the inlet

Flow- and Diffusion Distributed Structures with noise at the inlet Flow- and Diffusion Distributed Structures with noise at the inlet Pavel V. Kuptsov a,, Razvan A. Satnoianu b arxiv:0710.3812v2 [nlin.ps] 19 Nov 2007 a Department of Informatics, Saratov State Law Academy,

More information

On the mystery of differential negative resistance

On the mystery of differential negative resistance On the mystery of differential negative resistance Sebastian Popescu, Erzilia Lozneanu and Mircea Sanduloviciu Department of Plasma Physics Complexity Science Group Al. I. Cuza University 6600 Iasi, Romania

More information

Festival of the Mind Chaotic Chemical Waves: Oscillations and waves in chemical systems. Dr. Jonathan Howse Dr.

Festival of the Mind Chaotic Chemical Waves: Oscillations and waves in chemical systems. Dr. Jonathan Howse Dr. Festival of the Mind 2014 Chaotic Chemical Waves: Oscillations and waves in chemical systems Dr. Jonathan Howse Dr. Annette Taylor Mark Fell www.markfell.com Some history The BZ reaction In 1951, Belousov,

More information

Oregonator model of the Belousov-Zhabotinsky reaction Richard K. Herz,

Oregonator model of the Belousov-Zhabotinsky reaction Richard K. Herz, Oregonator model of the Belousov-Zhabotinsky reaction Richard K. Herz, rherz@ucsd.edu Boris Belousov in the 1950's discovered that a mixture of malonic acid, potassium bromate, and cerium sulfate an acidic

More information

Comparative Analysis of Packet and Trigger Waves Originating from a Finite Wavelength Instability

Comparative Analysis of Packet and Trigger Waves Originating from a Finite Wavelength Instability 11394 J. Phys. Chem. A 2002, 106, 11394-11399 Comparative Analysis of Packet and Trigger Waves Originating from a Finite Wavelength Instability Vladimir K. Vanag*, and Irving R. Epstein Department of Chemistry

More information

Heterogeneous Sources of Target Patterns in Reaction-Diffusion Systems

Heterogeneous Sources of Target Patterns in Reaction-Diffusion Systems J. Phys. Chem. 1996, 100, 19017-19022 19017 Heterogeneous Sources of Target Patterns in Reaction-Diffusion Systems Andrej E. Bugrim, Milos Dolnik, Anatol M. Zhabotinsky,* and Irving R. Epstein Department

More information

6. Well-Stirred Reactors III

6. Well-Stirred Reactors III 6. Well-Stirred Reactors III Reactors reaction rate or reaction velocity defined for a closed system of uniform pressure, temperature, and composition situation in a real reactor is usually quite different

More information

Locking of Turing patterns in the chlorine dioxide iodine malonic acid reaction with one-dimensional spatial periodic forcing

Locking of Turing patterns in the chlorine dioxide iodine malonic acid reaction with one-dimensional spatial periodic forcing PCCP Dynamic Article Links View Online Cite this: Phys. Chem. Chem. Phys., 2011, 13, 12578 12583 www.rsc.org/pccp PAPER Locking of Turing patterns in the chlorine dioxide iodine malonic acid reaction with

More information

Computational models: Reaction Diffusion. Matthias Vigelius Week 6

Computational models: Reaction Diffusion. Matthias Vigelius Week 6 Computational models: Reaction Diffusion Matthias Vigelius Week 6 Reactions (Chemical) reactions Reaction is a transformation from reactants into products: A + 2B D + 3E + k k [ ] Reactions occur with

More information

Pattern formation arising from wave instability in a simple reaction-diffusion system

Pattern formation arising from wave instability in a simple reaction-diffusion system Pattern formation arising from wave instability in a simple reaction-diffusion system Anatol M. Zhabotinsky, Milos Dolnik, and Irving R. Epstein Department of Chemistry and Volen Center for Complex Systems,

More information

TWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations

TWO DIMENSIONAL FLOWS. Lecture 5: Limit Cycles and Bifurcations TWO DIMENSIONAL FLOWS Lecture 5: Limit Cycles and Bifurcations 5. Limit cycles A limit cycle is an isolated closed trajectory [ isolated means that neighbouring trajectories are not closed] Fig. 5.1.1

More information

Dependence of Wave Speed on Acidity and Initial Bromate Concentration in the Belousov-Zhabotinsky Reaction-Diffusion System

Dependence of Wave Speed on Acidity and Initial Bromate Concentration in the Belousov-Zhabotinsky Reaction-Diffusion System Subscriber access provided by BRANDEIS UNIV Dependence of Wave Speed on Acidity and Initial Bromate Concentration in the Belousov-Zhabotinsky Reaction-Diffusion System Matthew D. Eager, Marianne Santos,

More information

Oscillatory Turing Patterns in a Simple Reaction-Diffusion System

Oscillatory Turing Patterns in a Simple Reaction-Diffusion System Journal of the Korean Physical Society, Vol. 50, No. 1, January 2007, pp. 234 238 Oscillatory Turing Patterns in a Simple Reaction-Diffusion System Ruey-Tarng Liu and Sy-Sang Liaw Department of Physics,

More information

MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT

MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT MAGNETIC NOZZLE PLASMA EXHAUST SIMULATION FOR THE VASIMR ADVANCED PROPULSION CONCEPT ABSTRACT A. G. Tarditi and J. V. Shebalin Advanced Space Propulsion Laboratory NASA Johnson Space Center Houston, TX

More information

Jumping solitary waves in an autonomous reaction diffusion system with subcritical wave instability

Jumping solitary waves in an autonomous reaction diffusion system with subcritical wave instability COMMUNICATION www.rsc.org/pccp Physical Chemistry Chemical Physics Jumping solitary waves in an autonomous reaction diffusion system with subcritical wave instability Lingfa Yang, Anatol M. Zhabotinsky

More information

Controlling chemical oscillations in heterogeneous Belousov-Zhabotinsky gels via mechanical strain

Controlling chemical oscillations in heterogeneous Belousov-Zhabotinsky gels via mechanical strain PHYSICAL REVIEW E 79, 046214 2009 Controlling chemical oscillations in heterogeneous Belousov-Zhabotinsky gels via mechanical strain Victor V. Yashin, 1 Krystyn J. Van Vliet, 2 and Anna C. Balazs 1, *

More information

Dynamics and Mechanism of Bromate Oscillators with 1,4-Cyclohexanedione

Dynamics and Mechanism of Bromate Oscillators with 1,4-Cyclohexanedione 10074 J. Phys. Chem. A 2003, 107, 10074-10081 Dynamics and Mechanism of Bromate Oscillators with 1,4-Cyclohexanedione István Szalai, Krisztina Kurin-Cso1rgei, Irving R. Epstein,*, and Miklós Orbán*, Department

More information

EXCITATION OF GÖRTLER-INSTABILITY MODES IN CONCAVE-WALL BOUNDARY LAYER BY LONGITUDINAL FREESTREAM VORTICES

EXCITATION OF GÖRTLER-INSTABILITY MODES IN CONCAVE-WALL BOUNDARY LAYER BY LONGITUDINAL FREESTREAM VORTICES ICMAR 2014 EXCITATION OF GÖRTLER-INSTABILITY MODES IN CONCAVE-WALL BOUNDARY LAYER BY LONGITUDINAL FREESTREAM VORTICES Introduction A.V. Ivanov, Y.S. Kachanov, D.A. Mischenko Khristianovich Institute of

More information

Controlling Coupled Chemical Oscillators: Toward Synchronization Engineering and. Chemical Computation

Controlling Coupled Chemical Oscillators: Toward Synchronization Engineering and. Chemical Computation Controlling Coupled Chemical Oscillators: Toward Synchronization Engineering and Chemical Computation Senior Thesis Presented to The Faculty of the School of Arts and Sciences Brandeis University Undergraduate

More information

Chemically Driven Convection in the Belousov-Zhabotinsky Reaction Evolutionary Pattern Dynamics

Chemically Driven Convection in the Belousov-Zhabotinsky Reaction Evolutionary Pattern Dynamics Review Forma, 30, S33 S53, 2015 Chemically Driven Convection in the Belousov-Zhabotinsky Reaction Evolutionary Pattern Dynamics Hidetoshi Miike 1,Tatsunari Sakurai 2, Atsushi Nomura 3 and Stefan C. Müller

More information

Krisztina Kurin-Cso1rgei, Anatol M. Zhabotinsky,*, Miklós Orbán,*, and Irving R. Epstein

Krisztina Kurin-Cso1rgei, Anatol M. Zhabotinsky,*, Miklós Orbán,*, and Irving R. Epstein J. Phys. Chem. 1996, 100, 5393-5397 5393 Bromate-1,4-Cyclohexanedione-Ferroin Gas-Free Oscillating Reaction. 1. Basic Features and Crossing Wave Patterns in a Reaction-Diffusion System without Gel Krisztina

More information

Numerical simulations of the edge tone

Numerical simulations of the edge tone Numerical simulations of the edge tone I. Vaik, G. Paál Department of Hydrodynamic Systems, Budapest University of Technology and Economics, P.O. Box 91., 1521 Budapest, Hungary, {vaik, paal}@vizgep.bme.hu

More information

PhD Thesis. Investigation of the oxidation of iron(iii) complexes of porphyrin derivatives by bromate. Dénesné Rácz Krisztina

PhD Thesis. Investigation of the oxidation of iron(iii) complexes of porphyrin derivatives by bromate. Dénesné Rácz Krisztina PhD Thesis Investigation of the oxidation of iron(iii) complexes of porphyrin derivatives by bromate Dénesné Rácz Krisztina University L. Eötvös, Institute of Chemistry 2008 I. Introduction Nonlinear phenomena

More information

Control of waves, patterns and turbulence in chemical systems

Control of waves, patterns and turbulence in chemical systems Physics Reports 425 (2006) 79 194 www.elsevier.com/locate/physrep Control of waves, patterns and turbulence in chemical systems Alexander S. Mikhailov a,, Kenneth Showalter b a Abteilung Physikalische

More information

The new physiology of vision-chapter The colours of interference. SIR C V RAMAN Received October 14, 1965

The new physiology of vision-chapter The colours of interference. SIR C V RAMAN Received October 14, 1965 Proc. Indian Acad. Sci. A62 243-248 (1965) The new physiology of vision-chapter The colours of interference XXVII. SIR C V RAMAN Received October 14, 1965 The characteristic features and properties of

More information

Control of Neo-classical tearing mode (NTM) in advanced scenarios

Control of Neo-classical tearing mode (NTM) in advanced scenarios FIRST CHENGDU THEORY FESTIVAL Control of Neo-classical tearing mode (NTM) in advanced scenarios Zheng-Xiong Wang Dalian University of Technology (DLUT) Dalian, China Chengdu, China, 28 Aug, 2018 Outline

More information

Nonlinear Dynamical Response of Ferrofluid to Magnetic Field. Kipp Schoenwald, Daniel Potter, Amir Hamid and Amin Agha

Nonlinear Dynamical Response of Ferrofluid to Magnetic Field. Kipp Schoenwald, Daniel Potter, Amir Hamid and Amin Agha Nonlinear Dynamical Response of Ferrofluid to Magnetic Field Kipp Schoenwald, Daniel Potter, Amir Hamid and Amin Agha Previous Work and Objective Mahr and Rehberg, "Nonlinear Dynamics of a Single Ferro-fluid

More information

A SPIRAL STRUCTURE FOR ELEMENTARY PARTICLES

A SPIRAL STRUCTURE FOR ELEMENTARY PARTICLES A SPIRAL STRUCTURE FOR ELEMENTARY PARTICLES Suraj Kumar Government Holkar Science College, A.B. Road, Indore (M.P.) 452017 INDIA Abstract In this paper we have tried to deduce the possible origin of particle

More information

CHEM-UA 652: Thermodynamics and Kinetics

CHEM-UA 652: Thermodynamics and Kinetics CHEM-UA 65: Thermodynamics and Kinetics Notes for Lecture I. THE COMPLEXITY OF MULTI-STEP CHEMICAL REACTIONS It should be clear by now that chemical kinetics is governed by the mathematics of systems of

More information

Two-stage Rydberg charge exchange in a strong magnetic field

Two-stage Rydberg charge exchange in a strong magnetic field Two-stage Rydberg charge exchange in a strong magnetic field M. L. Wall, C. S. Norton, and F. Robicheaux Department of Physics, Auburn University, Auburn, Alabama 36849-5311, USA Received 21 June 2005;

More information

Well Stirred Reactor Stabilization of flames

Well Stirred Reactor Stabilization of flames Well Stirred Reactor Stabilization of flames Well Stirred Reactor (see books on Combustion ) Stabilization of flames in high speed flows (see books on Combustion ) Stabilization of flames Although the

More information

Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore

Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Lecture No. # 26 Problem solving : Heterogeneous reactions Friends, in last few

More information

Cubic autocatalytic reaction diffusion equations: semi-analytical solutions

Cubic autocatalytic reaction diffusion equations: semi-analytical solutions T&T Proof 01PA003 2 January 2002 10.1098/rspa.2001.0899 Cubic autocatalytic reaction diffusion equations: semi-analytical solutions By T. R. Marchant School of Mathematics and Applied Statistics, The University

More information

Reactions. John Vincent Department of Chemistry University of Alabama

Reactions. John Vincent Department of Chemistry University of Alabama Oscillating Chemical Reactions John Vincent Department of Chemistry University of Alabama Kinetics In kinetics we study the rate at which a chemical process occurs. Besides information about the speed

More information

Effects of self-steepening and self-frequency shifting on short-pulse splitting in dispersive nonlinear media

Effects of self-steepening and self-frequency shifting on short-pulse splitting in dispersive nonlinear media PHYSICAL REVIEW A VOLUME 57, NUMBER 6 JUNE 1998 Effects of self-steepening and self-frequency shifting on short-pulse splitting in dispersive nonlinear media Marek Trippenbach and Y. B. Band Departments

More information

Spatiotemporal pattern formation in a prey-predator model under environmental driving forces

Spatiotemporal pattern formation in a prey-predator model under environmental driving forces Home Search Collections Journals About Contact us My IOPscience Spatiotemporal pattern formation in a prey-predator model under environmental driving forces This content has been downloaded from IOPscience.

More information

Experiment 7: SIMULTANEOUS EQUILIBRIA

Experiment 7: SIMULTANEOUS EQUILIBRIA Experiment 7: SIMULTANEOUS EQUILIBRIA Purpose: A qualitative view of chemical equilibrium is explored based on the reaction of iron(iii) ion and thiocyanate ion to form the iron(iii) thiocyanate complex

More information

Turing mechanism for patterning

Turing mechanism for patterning Turing mechanism for patterning Hyun Youk Delft University of Technology (Dated: July 23, 2018). I. MAIN IDEA OF THE TURING MECHANISM The Turing mechanism for spatial-pattern formation describes a uniformly

More information

Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores

Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores PHYSICAL REVIEW E VOLUME 62, NUMBER 1 JULY 2000 Velocity cross-correlations and atomic momentum transfer in simple liquids with different potential cores A. Verdaguer and J. A. Padró Departament de Física

More information

Period Doubling Cascade in Diffusion Flames

Period Doubling Cascade in Diffusion Flames Period Doubling Cascade in Diffusion Flames Milan Miklavčič Department of Mathematics, Michigan State University, East Lansing, MI 48824, USA Combustion Theory and Modelling 11 No 1 (2007), 103-112 Abstract

More information

Stochastic resonance in the absence and presence of external signals for a chemical reaction

Stochastic resonance in the absence and presence of external signals for a chemical reaction JOURNAL OF CHEMICAL PHYSICS VOLUME 110, NUMBER 7 15 FEBRUARY 1999 Stochastic resonance in the absence and presence of external signals for a chemical reaction Lingfa Yang, Zhonghuai Hou, and Houwen Xin

More information

PATTERN FORMATION IN NOISY SELF-REPLICATING SPOTS

PATTERN FORMATION IN NOISY SELF-REPLICATING SPOTS PTTERN FORMTION IN NOISY SELF-REPLICTING SPOTS ndreea Munteanu 1, and Ricard V. Solé 1, 2 1 ICRE-Complex Systems Lab, Universitat Pompeu Fabra (GRIB), Dr iguader 8, 83 Barcelona, Spain 2 Santa Fe Institute,

More information

3/9/2011. Outline Chapter 7 Waves Water Waves Water Waves. Water waves are really circular. They are an example of Mechanical waves.

3/9/2011. Outline Chapter 7 Waves Water Waves Water Waves. Water waves are really circular. They are an example of Mechanical waves. Outline Chapter 7 Waves 7-1. Water Waves 7-2. Transverse and Longitudinal Waves 7-3. Describing Waves 7-4. Standing Waves 7-5. Sound 7-6. Doppler Effect 7-7. Musical Sounds 7-8. Electromagnetic Waves 7-9.

More information

Two-Dimensional simulation of thermal blooming effects in ring pattern laser beam propagating into absorbing CO2 gas

Two-Dimensional simulation of thermal blooming effects in ring pattern laser beam propagating into absorbing CO2 gas Two-Dimensional simulation of thermal blooming effects in ring pattern laser beam propagating into absorbing CO gas M. H. Mahdieh 1, and B. Lotfi Department of Physics, Iran University of Science and Technology,

More information

The Reynolds experiment

The Reynolds experiment Chapter 13 The Reynolds experiment 13.1 Laminar and turbulent flows Let us consider a horizontal pipe of circular section of infinite extension subject to a constant pressure gradient (see section [10.4]).

More information

Plasma heating in stellarators at the fundamental ion cyclotron frequency

Plasma heating in stellarators at the fundamental ion cyclotron frequency PHYSICS OF PLASMAS VOLUME 7, NUMBER FEBRUARY 000 Plasma heating in stellarators at the fundamental ion cyclotron frequency V. A. Svidzinski and D. G. Swanson Department of Physics, Auburn University, Auburn,

More information

TORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR

TORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR TORSION PENDULUM: THE MECHANICAL NONLINEAR OSCILLATOR Samo Lasič, Gorazd Planinšič,, Faculty of Mathematics and Physics University of Ljubljana, Slovenija Giacomo Torzo, Department of Physics, University

More information

Dynamics of plumes driven by localized heating in a stably stratified ambient

Dynamics of plumes driven by localized heating in a stably stratified ambient Dynamics of plumes driven by localized heating in a stably stratified ambient Abstract Juan M. Lopez 1 and Francisco Marques 2 1 School of Mathematical and Statistical Sciences, Arizona State University,

More information

Systems Biology: A Personal View XX. Biological oscillators: Hopf Bifurcation & Glycolysis. Sitabhra Sinha IMSc Chennai

Systems Biology: A Personal View XX. Biological oscillators: Hopf Bifurcation & Glycolysis. Sitabhra Sinha IMSc Chennai Systems Biology: A Personal View XX. Biological oscillators: Hopf Bifurcation & Glycolysis Sitabhra Sinha IMSc Chennai www.bio.brandeis.edu Oscillations in the Chlorine dioxide-iodine- Malonic Acid (ClO

More information

Course Name: AP Physics. Team Names: Jon Collins. Velocity Acceleration Displacement

Course Name: AP Physics. Team Names: Jon Collins. Velocity Acceleration Displacement Course Name: AP Physics Team Names: Jon Collins 1 st 9 weeks Objectives Vocabulary 1. NEWTONIAN MECHANICS and lab skills: Kinematics (including vectors, vector algebra, components of vectors, coordinate

More information

arxiv:chao-dyn/ v1 5 Mar 1996

arxiv:chao-dyn/ v1 5 Mar 1996 Turbulence in Globally Coupled Maps M. G. Cosenza and A. Parravano Centro de Astrofísica Teórica, Facultad de Ciencias, Universidad de Los Andes, A. Postal 26 La Hechicera, Mérida 5251, Venezuela (To appear,

More information

Self-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model

Self-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model Letter Forma, 15, 281 289, 2000 Self-Replication, Self-Destruction, and Spatio-Temporal Chaos in the Gray-Scott Model Yasumasa NISHIURA 1 * and Daishin UEYAMA 2 1 Laboratory of Nonlinear Studies and Computations,

More information

Point Vortex Dynamics in Two Dimensions

Point Vortex Dynamics in Two Dimensions Spring School on Fluid Mechanics and Geophysics of Environmental Hazards 9 April to May, 9 Point Vortex Dynamics in Two Dimensions Ruth Musgrave, Mostafa Moghaddami, Victor Avsarkisov, Ruoqian Wang, Wei

More information

The Study of Chemical Reactions. Mechanism: The complete, step by step description of exactly which bonds are broken, formed, and in which order.

The Study of Chemical Reactions. Mechanism: The complete, step by step description of exactly which bonds are broken, formed, and in which order. The Study of Chemical Reactions Mechanism: The complete, step by step description of exactly which bonds are broken, formed, and in which order. Thermodynamics: The study of the energy changes that accompany

More information

8.1 Bifurcations of Equilibria

8.1 Bifurcations of Equilibria 1 81 Bifurcations of Equilibria Bifurcation theory studies qualitative changes in solutions as a parameter varies In general one could study the bifurcation theory of ODEs PDEs integro-differential equations

More information

Turing patterns in a simple gel reactor

Turing patterns in a simple gel reactor Physica A 188 (1992) 17-25 North-Holland Turing patterns in a simple gel reactor R. Dennis Vigil, Q. Ouyang and Harry L. Swinney Center for Nonlinear Dynamics, Department of Physics, The University of

More information

Electromagnetic fields and waves

Electromagnetic fields and waves Electromagnetic fields and waves Maxwell s rainbow Outline Maxwell s equations Plane waves Pulses and group velocity Polarization of light Transmission and reflection at an interface Macroscopic Maxwell

More information

Coupled oscillations in a 1D emulsion of Belousov Zhabotinsky droplets

Coupled oscillations in a 1D emulsion of Belousov Zhabotinsky droplets Soft Matter View Online Dynamic Article Links C < Cite this: Soft Matter, 2011, 7, 3155 www.rsc.org/softmatter Coupled oscillations in a 1D emulsion of Belousov Zhabotinsky droplets Jorge Delgado, a Ning

More information

Abstract: Complex responses observed in an experimental, nonlinear, moored structural

Abstract: Complex responses observed in an experimental, nonlinear, moored structural AN INDEPENDENT-FLOW-FIELD MODEL FOR A SDOF NONLINEAR STRUCTURAL SYSTEM, PART II: ANALYSIS OF COMPLEX RESPONSES Huan Lin e-mail: linh@engr.orst.edu Solomon C.S. Yim e-mail: solomon.yim@oregonstate.edu Ocean

More information

Chemistry 11. Unit 3 The Physical Properties and Physical Changes of Substances

Chemistry 11. Unit 3 The Physical Properties and Physical Changes of Substances Chemistry 11 1 Unit 3 The Physical Properties and Physical Changes of Substances 2 1. Definitions in science Science is the observation, identification, description, experimental investigation, and theoretical

More information

36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson

36. TURBULENCE. Patriotism is the last refuge of a scoundrel. - Samuel Johnson 36. TURBULENCE Patriotism is the last refuge of a scoundrel. - Samuel Johnson Suppose you set up an experiment in which you can control all the mean parameters. An example might be steady flow through

More information

EXPERIMENT 2-6. e/m OF THE ELECTRON GENERAL DISCUSSION

EXPERIMENT 2-6. e/m OF THE ELECTRON GENERAL DISCUSSION Columbia Physics: Lab -6 (ver. 10) 1 EXPERMENT -6 e/m OF THE ELECTRON GENERAL DSCUSSON The "discovery" of the electron by J. J. Thomson in 1897 refers to the experiment in which it was shown that "cathode

More information

Clearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e.

Clearly the passage of an eigenvalue through to the positive real half plane leads to a qualitative change in the phase portrait, i.e. Bifurcations We have already seen how the loss of stiffness in a linear oscillator leads to instability. In a practical situation the stiffness may not degrade in a linear fashion, and instability may

More information

Physics 142 Wave Optics 1 Page 1. Wave Optics 1. For every complex problem there is one solution that is simple, neat, and wrong. H.L.

Physics 142 Wave Optics 1 Page 1. Wave Optics 1. For every complex problem there is one solution that is simple, neat, and wrong. H.L. Physics 142 Wave Optics 1 Page 1 Wave Optics 1 For every complex problem there is one solution that is simple, neat, and wrong. H.L. Mencken Interference and diffraction of waves The essential characteristic

More information

THE problem of phase noise and its influence on oscillators

THE problem of phase noise and its influence on oscillators IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: EXPRESS BRIEFS, VOL. 54, NO. 5, MAY 2007 435 Phase Diffusion Coefficient for Oscillators Perturbed by Colored Noise Fergal O Doherty and James P. Gleeson Abstract

More information

Car-Following Models as Dynamical Systems and the Mechanisms for Macroscopic Pattern Formation

Car-Following Models as Dynamical Systems and the Mechanisms for Macroscopic Pattern Formation Car-Following Models as Dynamical Systems and the Mechanisms for Macroscopic Pattern Formation R. Eddie Wilson, University of Bristol EPSRC Advanced Research Fellowship EP/E055567/1 http://www.enm.bris.ac.uk/staff/rew

More information

Threshold of singularity formation in the semilinear wave equation

Threshold of singularity formation in the semilinear wave equation PHYSICAL REVIEW D 71, 044019 (2005) Threshold of singularity formation in the semilinear wave equation Steven L. Liebling Department of Physics, Long Island University-C.W. Post Campus, Brookville, New

More information

Self-Excited Vibration in Hydraulic Ball Check Valve

Self-Excited Vibration in Hydraulic Ball Check Valve Self-Excited Vibration in Hydraulic Ball Check Valve L. Grinis, V. Haslavsky, U. Tzadka Abstract This paper describes an experimental, theoretical model and numerical study of concentrated vortex flow

More information

(1) The time t required for N generations to elapse is merely:

(1) The time t required for N generations to elapse is merely: 19 Changes In Reactor Power With Time The two preceding modules discussed how reactivity changes increase or decrease neutron flux and hence, change the thermal power output from the fuel. We saw how the

More information

Time Resolved Pulsed Laser Photolysis Study of Pyrene Fluorescence Quenching by I - Anion

Time Resolved Pulsed Laser Photolysis Study of Pyrene Fluorescence Quenching by I - Anion 1 Time Resolved Pulsed Laser Photolysis Study of Pyrene Fluorescence Quenching by I - Anion Cameron Incognito, Ryan Bella, Cassandra Smith, Brandon Alexander Department of Chemistry, The Pennsylvania State

More information

The Bernoulli theorem relating velocities and pressures along a streamline comes from the steady momentum equation for a constant density fluid,

The Bernoulli theorem relating velocities and pressures along a streamline comes from the steady momentum equation for a constant density fluid, Flow Science Report 0-14 ADDING FLOW LOSSES TO FAVOR METHODOLOGY C.W. Hirt Flow Science, Inc. June 014 Background From time to time it has been observed that large velocities may appear in the vicinity

More information

Behavior of sink and source defects in a one-dimensional traveling finger pattern

Behavior of sink and source defects in a one-dimensional traveling finger pattern PHYSICAL REVIEW E, VOLUME 63, 066305 Behavior of sink and source defects in a one-dimensional traveling finger pattern Piotr Habdas, Matthew J. Case, and John R. de Bruyn Department of Physics and Physical

More information

Long Range Correlations in Chemical Oscillations

Long Range Correlations in Chemical Oscillations Journal of Control Science and Engineering 4 (2016) 95-101 doi: 10.17265/2328-2231/2016.02.006 D DAVID PUBLISHING Dola Roy Chowdhury 1, Sudeshna Lahiri 2, A. N. S. Iyengar 3 and M. S. Janaki 3 1. Techno

More information

196 7 atmospheric oscillations:

196 7 atmospheric oscillations: 196 7 atmospheric oscillations: 7.4 INTERNAL GRAVITY (BUOYANCY) WAVES We now consider the nature of gravity wave propagation in the atmosphere. Atmospheric gravity waves can only exist when the atmosphere

More information

Introduction LECTURE 1

Introduction LECTURE 1 LECTURE 1 Introduction The source of all great mathematics is the special case, the concrete example. It is frequent in mathematics that every instance of a concept of seemingly great generality is in

More information

Chapter 37 Early Quantum Theory and Models of the Atom. Copyright 2009 Pearson Education, Inc.

Chapter 37 Early Quantum Theory and Models of the Atom. Copyright 2009 Pearson Education, Inc. Chapter 37 Early Quantum Theory and Models of the Atom Planck s Quantum Hypothesis; Blackbody Radiation Photon Theory of Light and the Photoelectric Effect Energy, Mass, and Momentum of a Photon Compton

More information

5.1 2D example 59 Figure 5.1: Parabolic velocity field in a straight two-dimensional pipe. Figure 5.2: Concentration on the input boundary of the pipe. The vertical axis corresponds to r 2 -coordinate,

More information