Self-oscillations and chemical waves in CO oxidation on Pt and Pd: Kinetic Monte Carlo models

Theoretical studies of the spatiotemporal dynamics of CO oxidation on Pt(100) and Pd(110) single crystal surfaces have been carried out by the kinetic Monte Carlo method. For both surfaces, Monte Carlo simulation has revealed oscillations of the CO₂ formation rate and of the concentrations of adsorbed species. The oscillations are accompanied by wave processes on the model surface. Simulations have demonstrated that there is a narrow reaction zone when an oxygen wave propagates over the surface. The existence of this zone has been confirmed by experimental studies. Taking into account the anisotropy of the Pd(110) crystal has no effect on the oscillation period and amplitude, but leads to the formation of elliptic oxygen patterns on the surface. It is possible to obtain a wide variety of chemical waves (cellular and turbulent structures, spirals, rings, and strips) by varying the parameters of the computational experiment.     

Control and coupling of spiral waves in excitable media

We report the complex dynamics of spiral waves observed in the ferroin-catalyzed BZ reaction. The reaction is run in an open unstirred reactor (CFUR) with the catalyst immobilized on a polysulfone membrane. The catalyst-loaded membrane is placed between two well stirred compartments which are fed with solutions of sulfuric acid/malonic acid/bromide and sulfuric acid/bromate, respectively. An electrical field perpendicular to the membrane can be applied via Pt-ring electrodes or, alternatively, via transparent electrodes made of ITO-coated glass. In the field-free case relatively simple target and spiral patterns are observed in the membrane. If an alternating electrical field is applied the spiral core drifts through the membrane. The actual trajectory of the spiral tip depends on the amplitude and frequency of the applied electrical field. If the perturbation parameters are chosen properly the wave fronts break up and new spiral cores emerge under the influence of the alternating field. Complex spatio-temporal patterns may be induced which are reminiscent of "spiral-chaos". After switching off the perturbation the system returns to its previous, "simple" behaviour. Our experimental observations are confirmed by model calculations based on the Barkley model of spiral waves. The technique of using modulated excitability to control the dynamics of spiral waves is further extended to the coupling of two spirals in two CFURs. We present numerical simulations based on two identical excitable reaction-diffusion (RD) systems which are mutually coupled. The coupling is based on the observation of an arbitrarily chosen point inside each of the RD systems: If a chemical wave passes the point of optical observation in system 1 an electric field is applied to system 2 and vice versa. Thus, a local observation made in one system is transformed in a global perturbation of the second CFUR. We report the observation of CFUR states where the two spiral waves are spatially and temporally coupled to each other.     

Chemical turbulence and line defects induced by gradient effects in a three-dimensional reaction-diffusion system

We report experimental observations of chemical turbulence and line defects in a three-dimensional (3-D) Belousov-Zhabotinsky (BZ) reaction-diffusion system. Transitions from spiral waves to 3-D chemical turbulence to line defects are observed. These transitions are caused by concentration gradients across the third dimension in the 3-D reaction medium, indicating the observed line defects have a 3-D structure. The line defects come out of the 3-D turbulent state, and become smooth with the increase of the control parameter. Simulation with the two-variable Oregonator model in the 3-D system reproduces similar line defects.     

11,12,13-三甲基-1,4,7,10-四氮杂十三环-10,13-二烯-镍(II)配合物催化的新型振荡化学反应

本文报道了一种新型振荡化学反应体系--BrOˉ~3-MA-H2SO4-[Ni(L-H)]X(X为配阴离子)体系, 其中MA为丙二酸, L代表题给大环11,12,13-三甲基-1,4,7,10-四氮杂十三 环-10, 13-二烯。测得了该体系发生振荡的组分浓度范围; 研究了各物浓度、自由基抑制剂和还原剂对振荡反应的影响; 还测得了振荡的轨迹为螺旋线而非极限环, 由于催化剂消耗振荡衰减; 初步提出了振荡反应发生的可能机理。     

Chemical applications of group theory and topology

The following procedure is described for investigating the qualitative dynamics of simple chemical systems: 1) A so-called influence diagram is generated representing the relationships between the reference reactants (phase-determining intermediates); 2) This influence diagram is used to generate a truth table indicating possible transitions between state vectors representing the signs of the time derivatives of of the reference reactant concentrations; 3) The truth table is used to determine a state transition diagram representing the flow topology around unstable equilibrium points; 4) The characteristic equation of the adjacency matrix of the influence diagram is solved in order to determine the presence of such unstable equilibrium points. The two types of qualitative dynamics possible for chemical systems containing two reference reactants and one feedback circuit are bifurcation between two attracting regions (bistability) and limit cycle oscillation. However, in two reference reactant systems oscillation requires an additional self-activating loop to generate the unstable equilibrium point required for its realization. Bistability and limit cycle oscillation are also two of the possible types of qualitative dynamics for chemical systems containing three reference reactants. However, chemical systems with three reference reactants and two or more feedback circuits can also contain interlocking limit cycles, which can lead to toroidal oscillations or chaos. The influence diagrams are given for the systems exhibiting these various types of dynamic behavior along with a summary of the important properties of all 729 possible influences for simple chemical systems containing three reference reactants.     

硫化物在铂电极上电氧化时空动力学的模型模拟和分析(英文)

根据铂电极上硫化物电催化氧化的反应机理,本文提取动力学模型并利用数值模拟研究了N型负微分阻抗(N-NDR)振荡区域的电极表面时空反应动力学.在均相体系模拟中观察到电流简单振荡和复杂振荡,其来源于双电层电势自催化与传质限制和毒化物种吸附负反馈的相互耦合.为了更接近于真实体系,在模型中考虑了平行和垂直于电极表面两个方向的传质过程.模拟结果发现了与实验现象具有相同演化行为的复杂斑图,如行波和闪烁波;同时在传质耦合体系模拟中观察到双电层电势双臂螺旋波.本研究工作促进对电化学体系时空斑图的理解和预测.     

Chemical wave studies in the bromate–pyrocatechol beads system

We studied the chemical wave activity of the pyrocatechol‐acidic bromate system in the presence of ferroin‐loaded beads. The wave activity lasted for more than 24 h while meandering spirals continued for up to 10 h. Rigid and meandering spiral waves were investigated. We have analyzed the wave propagation speed and spiral tip trajectory versus the initial concentrations of all reagents as well as the age of the solution. Wave velocity depends on [H⁺] and [BrO] concentrations by the relationship v=k[H⁺]^1/2[BrO]^1/2, which is in agreement with other studies. This system is ideal to study wave activity and spiral waves as it does not produce precipitates under the studied conditions. © 2011 Wiley Periodicals, Inc. Int J Chem Kinet 43: 198–203, 2011     

Electronic coherences and vibrational wave-packets in single molecules studied with femtosecond phase-controlled spectroscopy

Employing femtosecond pulse-shaping techniques we investigate ultrafast, coherent and incoherent dynamics in single molecules at room temperature. In first experiments single molecules are excited into their purely electronic 0-0 transition by phase-locked double-pulse sequences with pulse durations of 75 fs and 20 nm spectral band width. Their femtosecond kinetics can then be understood in terms of a 2-level system and modelled with the optical Bloch equations. We find that we observe the coherence decay in single molecules, and the purely electronic dephasing times can be retrieved directly in the time domain. In addition, the Rabi-frequencies and thus the transition dipole moments of single molecules are determined from these data. Upon excitation of single molecules into a vibrational level of the electronically excited state also incoherent intra-molecular vibrational relaxation is recorded. Increasing the spectral band width of the excitation pulses to up to 120 nm (resulting in a transform-limited pulse width of 15 fs) coherent superpositions of excited state vibrational modes, i.e. vibrational wave packets, are excited. The wave-packet oscillations in the excited state potential energy surface are followed in time by a phase-controlled pump-probe scheme, which permits to record wave packet interference, and to determine the energies of vibrational modes and their coupling strengths to the electronic transition.     

Effects of delayed linear electrical perturbation of the Belousov-Zhabotinsky reaction: a case of complex mixed mode oscillations in a batch reactor

Delayed feedbacks are quite common in many physical and biological systems and in particular many physiological systems. Delay can cause a stable system to become unstable and vice versa. One of the well-studied non-biological chemical oscillators is the Belousov-Zhabotinskii (BZ) reaction. This gives relaxation oscillations for a considerable period of time under batch conditions. This paper deals with the effect of perturbing the limit cycle oscillation of BZ reaction by employing a delayed electrical feedback to the system under batch conditions. The parameters chosen to study are external resistance and delay. For various resistances and delays the system was electrically perturbed and found to exhibit various complex mixed mode oscillations. The dynamic features are accounted for by the Oregonator model, with time delay incorporated in one of the variables. This revised version was published online in June 2006 with corrections to the Cover Date.     

Spiralling optical morphologies in spherulites of poly(hydroxybutyrate)

Some rather unusual optical morphologies in thin films of poly(hydroxybutyrate) in which a temperature gradient is imposed across the thickness of the film are reported. Spherulites in which the Maltese cross degenerates into a continuous spiral grow under these conditions, and the direction of the rotation of the spirals depends on the direction of heating. These morphologies are explained with the existing understanding of spherulite optics. The helicoidally twisting crystallites are modeled as twisting around an axis at a fixed angle to the radius of the spherulite (and the plane of the film). The possible implications for future, inclusive models of banding in spherulites are discussed. Further observations on the temperature dependence of the optical banding pattern in poly(hydroxybutyrate) are also reported, and an unexpected minimum in band spacing and fine optically visible fibrillar texture is discussed. © 2000 John Wiley & Sons, Inc. J Polym Sci B: Polym Phys 38: 1575–1583, 2000     

Response analysis in biochemical chain reactions with negative feedforward and feedbackward loops

Finite difference equations can be used to study the responses of biochemical chain reactions at any step of the chain to an external stimulus. In this study, we developed mathematical models for two hypothetical chain reactions involving loops to study the responses in the chain as the length of the chain gets longer, so called transient and steady state responses. The first model is for a chain with a negative feedforward loop, and the second one is for a chain that has a negative feedback loop. Although both of the models have the same steady state equations and values, we showed that the chain with negative feedforward and negative feedback loops can produce significantly different behaviors. The former can bring the chain into oscillations with various periods and eventually chaos when the feedback is strong enough as the length of the reaction chain increases, whereas the latter is not capable of producing oscillations and more complicated dynamics.     

Iodine-starch clathrate complexes under the impact of a low-frequency acoustic field

The results from research on the kinetics of destruction of clathrate complexes (amyloidine, amylopectoiodine and iodinol) in low-frequency acoustic waves from 5 to 25 Hz are presented. Features of the transfer of energy from low-frequency acoustic oscillations are determined. The ability of systems containing biologically active clathrate structures to interact in resonance with external acoustic stimuli according to the energy state of the active part of iodine-containing compexes is confirmed. In the first approximation, such compounds can be regarded as model biochemical systems.     

Photodynamic Action of Methylene Blue on the Paramecium Membrane

An electrophysiological study of photodynamic action on the Paramecium membrane was carried out. In the presence of methylene blue (MB), light-spot stimulation of an anterior and a posterior part induced a depolarization and a hyperpolarization of the membrane, respectively. Under voltage-clamping, the anterior stimulation induced an inward current, while the posterior stimulation induced an outward current. The amplitudes of these currents were dependent on the membrane potential. When K⁺ channels were blocked with Cs⁺ and tetraethylammonium (TEA⁺), the posterior outward current was inhibited, while the anterior inward current was not inhibited. Intracellular application of the Ca²⁺ chelator, 1,2 -bis (2-aminophenoxy) ethane- N,N,N',N' -tetraacetic acid (BAPTA) also inhibited the posterior outward current, but the anterior inward current was unaffected. These results suggest that photodynamic action on the Paramecium membrane primarily opens the Ca²⁺ channels and the following influx of Ca²⁺ activates the Ca²⁺-dependent K⁺ channels localized mainly on the posterior part of the membrane.     

On the mechanism initiating bursting oscillatory patterns during the pitting corrosion of a passive rotating iron-disc electrode in halide-containing sulphuric acid solutions

Bursting current oscillations (trains of current spikes interrupted by a steady-state passive current) were found to occur during the iron anodic polarization in halide-containing sulphuric acid solutions. These oscillations occur at potentials located within a partially passive state of the iron. The partially passive state appears at potentials more positive than the unstable passive-active transition region of the current-potential curve of the Fe|2 M H₂SO₄ system for halide concentrations approximately ranged between 20 and 50 mM. The bursting oscillations appear after a certain induction period which ranged from minutes to several hours. Experiments were conducted in which an external resistance was connected in series between the ground and the working electrode for the Fe|2 M H₂SO₄ + 30 mM Br⁻ system. The results demonstrate that the appearance of oscillations is closely related to an IR-drop mechanism operating during the pitting corrosion of Fe induced by halides. When an appropriate value of the external resistance was imposed in the system the induction period was entirely eliminated. Bursting oscillations appear at once after setting the potential at the desired value. These results provide strong evidence that during the anodic electrodissolution of iron in halide-containing sulphuric acid solutions bursting current oscillations will be induced even though the electrode seems to be in the partially passive state. This is expected to occur after the passage of an induction period during which the system acquires a critical value of the ohmic potential drop. Received: 20 October 1997 / Accepted: 2 January 1998     

Front waves in the NO + NH3 reaction on Pt{100}

In the present work, we spatially extended a brand new kinetic mechanism of the NO + NH3 reaction on Pt{100} to simulate the experimentally observed spatiotemporal traveling waves. The kinetic mechanism developed by Irurzun, Mola, and Imbihl (IMI model) improves the former model developed by Lombardo, Fink, and Imbihl (LFI model) by replacing several elementary steps to take into account experimental evidence published since the LFI model appeared. The IMI model achieves a better agreement with the experimentally observed dependence of the oscillation period on temperature. In the present work, the IMI model is extended by considering Fickean diffusion and coupling via the gas phase. Traveling waves propagating across the surface are obtained at realistic values of temperature and partial pressure. A transition from amplitude to phase waves is observed, induced either by temperature or by the gas global coupling strength. The traveling waves simulated in the present work are not associated with fixed defects, in agreement with experimental evidence of spiral centers capable of moving on the surface. Also, the IMI model adequately predicts the presence of macroscopic oscillations in the partial pressures of the reactants coexisting with front wave patterns on the surface.     

Transition from traveling to standing waves as a function of frequency in a reaction-diffusion system

The addition of polyethylene glycol to the Belousov-Zhabotinsky reaction increases the frequency of oscillations, which in an extended system causes a transition from traveling to standing waves. A further increase in frequency causes another transition to bulk oscillations. The standing waves are composed of two domains, which oscillate out of phase with a small delay between them, the delay being smaller as the frequency of oscillations is increased.     

Outward photocurrent component in chloroplasts of Peperomia metallica and its assignment to the 'closed thylakoid' recording configuration

The photoinduced electrical events at energy-conserving chloroplast membranes can be studied in whole plastids using suction electrodes. In chloroplasts of Peperomia metallica the kinetic profile of photocurrent contains a minor outward component that occurs prior to and differs in polarity from the main component. The origin of this outward current was analyzed using single-turnover flashes in combination with prolonged light exposures and differential physicochemical treatments of tip-located (internal) and the exposed parts of a chloroplast. The outward current signal was higher after 10- to 20-s preillumination and gradually reduced in darkness. The relative amplitude of the outward peak current was enhanced when photosystem II (PS II) was excited by flashes given in the presence of far-red background light (lambda = 712 nm). The outward current was small or absent under conditions promoting activity of photosystem I (cyclic electron transport supported by artificial redox mediators in the presence of diuron) and was particularly high in the presence of PS II electron acceptors (e.g., p-phenylenediamine). This indicates the predominant association of the outward current with activity of PS II. The external application of diuron strongly inhibited the inward current, giving rise to a temporal increase in the outward current. On the contrary, when diuron was added into the suction pipette, the outward current was inhibited soon after sealing. The data suggest that the outward current originated in the tip-located portions of the thylakoid membrane that have orientation opposite to the exposed part of 'whole thylakoid'. These tip-located membrane portions are least accessible for inhibitors added into the outer medium and are highly sensitive to inhibitors (diuron), ionophores (gramicidin D), and detergents (Triton X-100) added into the pipette. Differential involvement of two photosystems in generation of the outward current may be caused by uneven structural distribution of photosystems I and II between appressed (granal) and nonappressed (stromal) thylakoids and by different recording configurations for these thylakoids.     

Spatiotemporal chaos arising from standing waves in a reaction-diffusion system with cross-diffusion

We show that quasi-standing wave patterns appear in the two-variable Oregonator model of the Belousov-Zhabotinsky reaction when a cross-diffusion term is added, no wave instability is required in this case. These standing waves have a frequency that is half the frequency of bulk oscillations displayed in the absence of diffusive coupling. The standing wave patterns show a dependence on the systems size. Regular standing waves can be observed for small systems, when the system size is an integer multiple of half the wavelength. For intermediate sizes, irregular patterns are observed. For large sizes, the system shows an irregular state of spatiotemporal chaos, where standing waves drift, merge, and split, and also phase slips may occur.     

Understanding the electronic reorganization in the thermal isomerization reaction of trans‐3,4‐dimethylcyclobutene. Origins of outward Pseudodiradical {2n + 2π} torquoselectivity

The thermal isomerization reaction of trans‐3,4‐dimethylcyclobutene (1,2,3,4‐DMC) to produce the isomer (2E, 4E)‐hexadiene have been studied using density functional theory at the B3LYP/6‐31+G level. For this reaction, two different channels of the conrotatory torquoselectivity allowing the formation of the two isomeric (E,E) and (Z,Z) have been characterized. The isomer (Z,Z) occurs through the inward conrotatory mechanism, whereas the isomer (E,E) occurs through the outward conrotatory mechanism. The outward conrotatory is favored by 11.3 kcal/mol with respect to inward conrotatory one. This behavior is consistent with the thermodynamic properties: enthalpy, free energy, and entropy calculated in both reaction pathways. The topological analysis of the electron localization function at the outward conrotatory transition state explicated the electronic reorganization through a pseudodiradical {2n + 2π} process and not a pericyclic reorganization. © 2012 Wiley Periodicals, Inc.     

Settling and deformation of a thin elastic shell on a thin fluid layer lying on a solid surface

placement of a soft contact lens onto the cornea, the upper eyelid deforms and settles the lens by squeezing fluid out of the post-lens tear film or POLTF (i.e., the tear fluid layer sandwiched between the lens and the cornea). This paper studies the physical mechanisms that control the dynamic state of the contact lens during blinking, i.e., its shape and its distance from the cornea, especially a long time after insertion. We model the lens as a deformable elastic shell and the cornea as a flat nondeformable body. The tear fluid is assumed to be Newtonian, and the lens is characterized by an elastic modulus and a Poisson ratio. Lubrication equations under creeping flow are used to solve the fluid problem, while the thin-shell approximation is applied to the solid lens. The solid and fluid mechanics problems are coupled by maintaining continuity of stress and velocity at the solid/liquid interface. Lid applied pressure causes the lens to approach the cornea by squeezing tear fluid out and also leads to the deformation of the lens. Subsequently, in the interblink period, since there is no applied force, the elastic energy stored in the lens due to its deformation is released causing it to move away from the cornea by imbibing tear fluid into the POLTF. If the POLTF thickness is large, the inward motion of the lens in the blink is more than the outward motion during interblink, and this causes the lens to settle closer to the cornea. Eventually, there may be a balance of the inward motion during the blink and the outward motion during the interblink. If so, the lens subsequently exhibits periodic steady-state motion. However, it is also possible that a balance of inward and outward motion is never achieved, and the lens continues to settle endlessly. If this happens, then the thinfilm interactions between the mucin-covered corneal surface and the lens material determine whether the lens actually touches the cornea and possibly adheres. Our elastohydrodynamic analysis serves as a useful tool to elucidate the effects of various lens parameters on the final settled state of the lens. In particular, we are concerned about eventual adherence and/or mechanical abrasion to the cornea, which is very important to the ocular health of soft contact lens wearers.