Controlled excitations of the Belousov-Zhabotinsky reaction: Experimental procedures

The purpose of this research was to explore the unstirred, ferroin-catalyzed Belousov-Zhabotinsky (BZ) reaction as an experimental model for the response of excitable media to small perturbations (slightly larger than the threshold for excitations). Following Showalter et al. (Showalter, K.; Noyes, R. M.; Turner, H. J.Am. Chem. Soc. 1979, 101, 7463-69), we used a positively biased silver electrode to release silver ions into a BZ reaction mixture, removing bromide ions and causing an excitation if sufficient bromide was removed. We found (1) a scaling region in which the delay before activation increased linearly as the size of the perturbation decreased, qualitatively consistent with but not fully explained by the Oregonator of Field et al. (Field, R. J.; K?r?s, E.; Noyes, R. M. J. Am. Chem. Soc. 1972, 94, 8649-64); (2) evidence for a 10 s oligomerization time scale; and (3) that activations were always delayed until after the end of a pulse of current, with the delay essentially constant for sufficiently long pulses, an effect not seen in simple ODE models but consistent with the anomalously large current apparently required for activation (Showalter, K.; Noyes, R. M. J. Am. Chem. Soc. 1976, 98, 3730-31) and explainable by bromide transport. Overall, the BZ system appeared to be well-suited as an experimental prototype, despite its complexity.     

The role of radicals in the Belousov-Zhabotinsky reaction

A new theory of the Belousov-Zhabotinsky (BZ) reaction, the Radicalator model, is presented. This model is based on a negative feedback loop involving a fast reaction between malonyl and bromine dioxide radicals. Experimental evidence for the validity of the model is given for BZ systems in 3 M and 1 M sulfuric acid solution.     

Paramagnetic nuclear relaxation study of the structure and dynamics of lyotropic lamellar phases

Lamellar phases of double tailed sodium alkyl-phosphates/water systems have been investigated by ³¹P and ¹³C relaxations enhanced by paramagnetic divalent ions exchanging rapidly among the polar heads of the surfactant molecules. In order to extract both structural and dynamic information from these experiments, the paramagnetic probes have been chosen such that their electron spin-lattice relaxation times are significantly larger (Mn²⁺, VO²⁺) or smaller (Ni²⁺) than the time scales of molecular motions generally found in the 0.1-2.0 ns range. The paramagnetic relaxation of surfactant nuclei results from intra- and intermolecular contributions, the latter corresponding to dipolar interactions in a two dimensional lattice. The sum of these contributions has been computed as an average over all accessible conformers resulting from the trans-gauche isomerizations of alkyl chains and fitted to the experimental relaxation rates at several magnetic field strengths. This procedure allows the determination of the populations of the main conformers, of the reorientational correlation times and of the lateral diffusion coefficient of the surfactant, which are sometimes difficult to obtain by other N.M.R. methods.     

Paramagnetic nuclear relaxation study of the structure and dynamics of lyotropic lamellar phases

Abstract

Lamellar phases of double tailed sodium alkyl-phosphates/water systems have been investigated by ³¹ P and ¹³ C relaxations enhanced by paramagnetic divalent ions exchanging rapidly among the polar heads of the surfactant molecules. In order to extract both structural and dynamic information from these experiments, the paramagnetic probes have been chosen such that their electron spin-lattice relaxation times are significantly larger (Mn²⁺, VO²⁺) or smaller (Ni²⁺) than the time scales of molecular motions generally found in the 0.1-2.0 ns range. The paramagnetic relaxation of surfactant nuclei results from intra- and intermolecular contributions, the latter corresponding to dipolar interactions in a two dimensional lattice. The sum of these contributions has been computed as an average over all accessible conformers resulting from the trans-gauche isomerizations of alkyl chains and fitted to the experimental relaxation rates at several magnetic field strengths. This procedure allows the determination of the populations of the main conformers, of the reorientational correlation times and of the lateral diffusion coefficient of the surfactant, which are sometimes difficult to obtain by other N.M.R. methods.     

Noise-induced order in the chaos of the Belousov-Zhabotinsky reaction

Noise can stabilize a metastable state in such a way that the system remains in this state for a longer time than in the absence of noise. When this phenomenon is observed in chaos, it is called "noise-induced order." We have experimentally detected noise-induced order in the Belousov-Zhabotinsky reaction. That is, when noise is added to the chaos with the flow rate near the period-three oscillation, a decrease of the maximum Lyapunov exponent and a convergence of the Fourier spectrum are observed. Moreover, the analysis on the one-dimensional return map reveals that noise-induced order is caused by the convergence of the chaotic trajectory into the laminar region.     

Alignment of lamellar diblock copolymer phases under shear: insight from dissipative particle dynamics simulations

Sheared self-assembled lamellar phases formed by symmetrical diblock copolymers are investigated through dissipative particle dynamics simulations. Our intent is to provide insight into the experimental observations that the lamellar phases adopt parallel alignment at low shear rates and perpendicular alignment at high shear rates and that it is possible to use shear to induce a transition from the parallel to perpendicular alignment. Simulations are initiated either from lamellar structures prepared under zero shear where lamellae are aligned into parallel, perpendicular, or transverse orientations with respect to the shear direction or from a disordered melt obtained by energy minimization of a random structure. We first consider the relative stability of the parallel and perpendicular phases by applying shear to lamellar structures initially aligned parallel and perpendicular to the shear direction, respectively. The perpendicular lamellar phase persists for all shear rates investigated, whereas the parallel lamellar phase is only stable at low shear rates, and it becomes unstable at high shear rates. At the high shear rates, the parallel lamellar phase first transforms into an unstable diagonal lamellar phase; and upon further increase of the shear rate, the parallel lamellar phase reorients into a perpendicular alignment. We further determine the preferential alignment of the lamellar phases at low shear rate by performing the simulations starting from either the initial transverse lamellar structure or the disordered melt. Since the low shear-rate simulations are plagued by the unstable diagonal lamellar phases, we vary the system size to achieve the natural spacing of the lamellae in the simulation box. In such cases, the unstable diagonal lamellar phases disappear and lamellar phases adopt the preferential alignment, either parallel or perpendicular. In agreement with the experimental observations, the simulations show that the lamellar phase preferentially adopts the parallel orientation at low shear rates and the perpendicular orientation at high shear rates. The simulations further reveal that the perpendicular lamellar phase has lower internal energy than the parallel lamellar phase, whereas the entropy production of the perpendicular lamellar phase is higher with respect to the parallel lamellar phase. Values of the internal energy and entropy production for the unstable diagonal lamellar phases lie between the corresponding values for the parallel and perpendicular lamellar phases. These simulation results suggest that the relative stability of the parallel and perpendicular lamellar phases at low shear rates is a result of the interplay between competing driving forces in the system: (a) the system's drive to adopt a structure with the lowest internal energy and (b) the system's drive to stay in a stationary nonequilibrium state with the lowest entropy production.     

Mechanism of the transition between lamellar and gyroid phases formed by a diblock copolymer in aqueous solution

The mechanism of the transition from a lamellar phase to a gyroid phase in an aqueous solution of a diblock copolymer has been studied by time-resolved synchrotron small-angle X-ray scattering. The transition occurs via a metastable perforated lamellar structure. The perforations initially have liquidlike ordering before developing hexagonal packing. The transient phase of irregularly perforated layers is revealed by the development of diffuse scattering peaks, just below the Bragg peaks of the lamellar structure. The diffuse scattering is modeled by Monte Carlo simulations of perforated layers. Following the formation of perforations, Bragg peaks characteristic of a hexagonal structure signal an ordering into a hexagonal lattice (with the concomitant loss of diffuse scattering). Computer simulations based on a dynamic density functional model reproduce these features. The hexagonal perforated lamellar phase is rapidly replaced by the gyroid phase. The domain spacing of the gyroid phase is larger than that of the perforated lamellar structure. The perforated lamellar and gyroid phases coexist for a defined period. The reverse transition from gyroid to lamellae occurs directly, with no transient or metastable intermediates.     

The source of the carbon monoxide in the classical Belousov-Zhabotinsky reaction

CO and CO2 evolution was measured in a cerium and in a ferroin-catalyzed Belousov-Zhabotinsky (BZ) reaction. These gases were stripped from the reaction mixture by a N2 carrier gas, mixed with H2, converted to methane on a Ni catalyst, and then measured by a flame ionization detector (FID). CO could be detected separately by absorbing CO2 on a soda lime column. In separate experiments it was proven that CO is produced in a reaction of BrO2* radicals with bromomalonic acid (BrMA). To this end BrO2(.-) radicals were generated in two different ways: (i) in the reaction HBrO2 + HBrO3 <--> 2 BrO2(.-) + H2O and (ii) by reducing HBrO3 to BrO2(.-) by Fe(2+). It was found that (.-)OH radicals--produced by Fenton's reagent--can also generate CO from BrMA. We propose that CO can be formed when an inorganic radical (like BrO2(.-) or (.-)OH) reacts with the enol form of BrMA producing an acyl radical which decarbonylates in the next step. Malonic acid (MA)-BrMA mixtures were prepared by a new method modifying Zaikin and Zhabotinsky's original recipe to minimize the production of dibromomalonic acid (Br2MA).     

Malonic acid concentration as a control parameter in the kinetic analysis of the Belousov-Zhabotinsky reaction under batch conditions

The influence of the initial malonic acid concentration [MA]0 (8.00 x 10(-3) < or = [MA]0 < or = 4.30 x 10(-2) mol dm(-3)) in the presence of bromate (6.20 x 10(-2) mol dm(-3)), bromide (1.50 x 10(-5) mol dm(-3)), sulfuric acid (1.00 mol dm(-3)) and cerium sulfate (2.50 x 10(-3) mol dm(-3)) on the dynamics and the kinetics of the Belousov-Zhabotinsky (BZ) reactions was examined under batch conditions at 30.0 degrees C. The kinetics of the BZ reaction was analyzed by the earlier proposed method convenient for the examinations of the oscillatory reactions. In the defined region of parameters where oscillograms with only large-amplitude relaxation oscillations appeared, the pseudo-first order of the overall malonic acid decomposition with a corresponding rate constant of 2.14 x 10(-2) min(-1) was established. The numerical results on the dynamics and kinetics of the BZ reaction, carried out by the known skeleton model including the Br2O species, were in good agreement with the experimental ones. The already found saddle node infinite period (SNIPER) bifurcation point in transition from a stable quasi-steady state to periodic orbits and vice versa is confirmed by both experimental and numerical investigations of the system under consideration. Namely, the large-amplitude relaxation oscillations with increasing periods between oscillations in approaching the bifurcation points at the beginning and the end of the oscillatory domain, together with excitability of the stable quasi-steady states in their vicinity are obtained.     

On the problem of “mysterious” adducts in the Belousov-Zhabotinsky reaction

The formation of an adduct betweeno-phenanthroline and bromine, a feasible intermediate in the Belousov-Zhabotinsky reaction (BZ) catalyzed by ferroin, was studied. The stability constant of the adduct was determined. One of the assumptions on the composition of the “red precipitate,” the substance to which several differences between the behavior of the ferroin and cerium BZ reactions are due, was confirmed.     

Spatio-temporal perturbation of the dynamics of the ferroin catalyzed Belousov-Zhabotinsky reaction in a batch reactor caused by sodium dodecyl sulfate micelles

The effects of the anionic surfactant sodium dodecyl sulfate (SDS) on the spatio-temporal and temporal dynamics of the ferroin-catalyzed Belousov-Zhabotinsky (BZ) reaction have been studied over a wide surfactant concentration range. For the first time, investigations were performed also for unstirred systems. The presence of SDS in the reaction mixture influences the oscillatory parameters to an extent that significantly depends on the surfactant concentration. The trend of the wave speed v upon the increasing amount of SDS was found to have a maximum at [SDS] = 0.075 mol dm (-3) ( v = 0.071 mm s (-1)), after which the speed decreased to 0.043 mm s (-1) at [SDS] = 0.5 mol dm (-3), which is below the value found in the absence of the surfactant ( v = 0.055 mm s (-1)). The response of the oscillatory BZ system to the addition of SDS has been ascribed to two different causes: (a) the peculiar capability of the organized surfactant assemblies to affect the reactivity by selectively sequestering some key reacting species and (b) the modifications induced by SDS on the physical properties of the medium. These hypotheses have been corroborated by performing spectrophotometric investigations on the stirred BZ system. Complementary viscosity measurements gave useful hints for the clarification of the surfactant role.     

Pentanary cross-diffusion in water-in-oil microemulsions loaded with two components of the Belousov-Zhabotinsky reaction

We measure cross-diffusion coefficients in a five-component system, an aerosol OT (AOT) water-in-oil microemulsion loaded with two constituents of the Belousov-Zhabotinsky (BZ) reaction (H(2) O/AOT/BZ1/BZ2/octane). The species BZ1 is either NaBr, an inhibitor of the BZ reaction, or ferroin, a catalyst for the reaction. As species BZ2, we choose Br(2) , an intermediate in the reaction. The cross-diffusion coefficients between BZ1 and BZ2 are found to be negative, which can be understood in terms of complexation between these species. Using a four-variable model for the BZ reaction, we find that the cross-diffusion coefficients measured here can lead to a noticeable shift in the onset of Turing instability in the BZ-AOT system.     

Formation of GeO2 nanosheets using water thin layers in lamellar phase as a confined reaction field--in situ measurement of SAXS by synchrotron radiation

Highly crystallized GeO2 nanosheets were synthesized by hydrolysis and condensation reactions of germanium alkoxide using a 2-dimensional flat thin lamellar phase water layer containing surfactant molecules at the liquid-liquid interface as a confined reaction field.     

Semi-quantitative determination of cationic surfactants in aqueous solutions using gold nanoparticles as reporter probes

Concentrations of cationic surfactants in aqueous solutions have been estimated on the basis of changes in the color of gold nanoparticles, used as reporter probes. We have shown that the colors of gold nanoparticles with anionic protective groups on their surfaces shift from red to indigo/purple and then back to red in a range of cationic surfactant solutions in which concentrations vary from very low to above the theoretical CMCs. The color changes occur near the theoretical CMCs, presumably because the presence of surfactant micelles in the solution prevents the gold nanoparticles from aggregating. We have used gold nanoparticles as reporter probes to determine the concentrations of cationic surfactants in products such as hair conditioners, which often contain large amounts of alkyltrimethylammonium halides. Although this approach can only provide an estimate, it can be performed simply by addition of a given amount of gold nanoparticles to a series of diluted solutions, without the need for instruments or labor-intensive procedures. Figure Photographs of a series of diluted hair conditioner solutions with added gold nanoparticles

---

     

Fluorescent probe as reporter on the local structure and dynamics in hydrolysis-condensation process of organotrialkoxysilanes

To get some information on the aggregation behaviors of the products derived from different organotrialkoxysilanes, the hydrolysis-condensation processes of some organotrialkoxysilanes have been examined by means of pyrene as fluorescent probe. The organotrialkoxysilanes used in the research were n-octadecyltri-methoxysilane (ODTMS), n-octyltrimethoxysilane (OTMS), 3-glycidoxypropyltrimethoxysilane (GTMS), 3-methacryloxypropyltrimethoxysilane (MAPTMS), and propyltrimethoxy-silane (PTMS). The results show that pyrene as fluorescence probe can respond sensitively not only to the organization state of the hydrolysates but also to the change in the organization state during the condensation process. The organization states during the hydrolysis and condensation can be explained in terms of structures of the products. In the initial stage, the silanols with long organic chains are amphiphilic molecules, and such nature of the silanols can be compared to that of a surfactant. Therefore, the excimer emission of pyrene is extremely obvious because of such silanols being prone to form aggregates. In the case of silanols having short alkyl groups or epoxy groups, these silanols homogenously disperse in solution, which results in the appearance of an only monomer emission of pyrene. In the late stage, the fluorescence behavior of pyrene is also sensitive to structural evolution of the silicates. The fluorescence spectra of pyrene during the condensation of the silanols with short alkyl groups or epoxy groups are almost in silence, indicating that the condensation products, with a low condensation degree, homogeneously disperse in solution. For the silanols with long hydrophobic substituents in different lengths, the changes in fluorescence spectra of pyrene during the condensation are varied. Commonly, the excimer emission is noticeable, implying that the condensation products with high condensation degree inhomogenously disperse in solution. However, the relative excimer/monomer fluorescence intensity is alkyl chain-length dependent. The longer alkyl chains in the condensation products result in the appearance of the obvious excimer emission. These phenomena imply that the condensation degree of the products increases with the length of the alkyl chains. Additionally, the distorted spectrum of pyrene appears in the case of the organotrialkoxysilanes with side chain substituent, illustrating that the steric hindrance between the substituents can be monitored by fluorescence of pyrene. All these results are verified by the fluorescence-quenching measurements. The approach in the present study gives new insights into the local structure and dynamics in hydrolysis-condensation process of organotrialkoxysilanes and emphasizes the influence of the self-assembling behavior.     

Electrolyte-induced mesoscopic aggregation of thiacarbocyanine dye in aqueous solution: counterion size specificity

Countercation size specificity is presented for the electrolyte-induced aggregation of 3,3'-disulfopropyl-5,5'-dichloro-9-methyl thiacarbocyanine (TCC) dye in aqueous solution. Addition of electrolytes having a small monovalent cation (Na+, NH4+, or Cs+) induced pure H aggregates of TCC, whereas J aggregates were preferentially promoted by electrolytes with a large monovalent cation ([N(CH3)4]+ or [N(C2H5)4]+). The electrolyte-induced H aggregate (HS aggregate) differed spectroscopically from that spontaneously self-assembled in aqueous solution. Mesoscopic structure of the HS aggregates was revealed via polarized-light microscopy and atomic force microscopy; a rodlike morphology of 50-70 nm wide and tens to hundreds of micrometers long with very strong negative birefringence. A simple structural model based on semiempirical molecular orbital calculations can explain the aggregation behaviors: The anionic TCC monomer shows a considerable planar geometry between two benzothiazole end groups when it involves a sodium cation, which favors the H-type molecular arrangements in a face-to-face orientation. On the other hand, the TCC dye has a twisted conformation when it implicates a large tetramethylammonium cation, resulting in the formation of the J aggregates.     

Characterization of the properties of stationary phases for liquid chromatography in aqueous mobile phases using aromatic sulphonic acids as the test compounds

We investigated the effects of the concentration of naphthalene sulphonic acids (NSAs) as anionic test compounds in the injected sample and of the salt additives to the mobile phase on ion-exclusion. The retention behaviour of NSAs sensitively reflects even minor changes in the ionic and hydrophobic interactions and can be useful for predicting the effects of the stationary phases in reversed-phase chromatography of polar and ionic compounds, both small ones and biopolymers, e.g., oligonucleotides. We studied chromatographic properties of several stationary phases intended for separations in aqueous mobile phases: a C18 column end-capped with polar hydrophilic groups, a densely bonded C8 column doubly end-capped with short alkyl groups, a short alkyl stationary phase designed to keep full pore accessibility in highly-aqueous mobile phases and a Bidentate column with “bridged” C18 groups attached to the silica hydride support. The chemistry and pore structure of various types of column packing materials and of the salt additives to the mobile phase affect the proportion of the pore volume non-accessible to anions due to ion-exclusion and consequently the peak asymmetry and hydrophobic selectivity in reversed-phase chromatography of organic acids. We also addressed the problems connected with the determination of column hold-up volume in aqueous mobile phases. The accessibility of the stationary phase for anionic compounds in contact with the sample zone is affected by ion-exclusion due to repulsive interactions with the negatively charged surface in the pores of the stationary phase. The accessible part of the stationary phase increases and consequently the migration velocity along the column decreases with increasing concentration of the sample in the zone moving along the column. Because of a limited access to the stationary phase, its capacity can be easily overloaded. The combination of the column overload and ion-exclusion effects may result in fronting or tailing peak asymmetry. To explain this behaviour, we proposed a modified Langmuir model, respecting the variation of the column capacity due to the effects of sample concentration on ion-exclusion.     

Estimation of the activation energy in the Belousov-Zhabotinsky reaction by temperature effect on excitable waves

We report the temperature effect on the propagation of excitable traveling waves in a quasi-two-dimensional Belousov-Zhabotinsky reaction-diffusion system. The onset of excitable waves as a function of the sulfuric acid concentration and temperature is identified, on which the sulfuric acid concentration exhibits an Arrhenius dependence on temperature. On the basis of this experimental data, the activation energy of the self-catalyzed reaction in the Oregonator model is estimated to be 83-113 kJ/mol, which is further supported by our numerical simulations. The estimation proceeds without analyzing detailed reaction steps but rather through observing the global dynamic behaviors in the BZ reaction. For a supplement, the wave propagation velocities are calculated based on our results and compared with the experimental observations.     

Coexistence of wave propagation and oscillation in the photosensitive Belousov-Zhabotinsky reaction on a circular route

The photosensitive Belousov-Zhabotinsky (BZ) reaction was investigated on a circular ring, which was drawn using computer software and then projected on a film soaked with BZ solution using a liquid-crystal projector. Under the initial conditions, a chemical wave propagated with a constant velocity on the black ring under a bright background. When the background was rapidly changed to dark, coexistence of the oscillation on part of the ring and propagation of the chemical wave on the other part was observed. These experimental results are discussed in relation to the nature of the photosensitive BZ reaction and theoretically reproduced based on a reaction-diffusion system using the modified Oregonator model.