Bulbul Chakraborty, PhD , Brandeis University
We have achieved a numberof new results on the behavior of repulsively coupled chemical oscillators, andmade significant progress in developing a framework for stochastic differentialequations that model reactiondiffusion systems.
Fig 1: (Top) Frequency order parameter and bifurcation diagram illustrating phase synchronization with increasing K. (Bottom) (a) The different helicity patterns, (b) Domains of different helicities (white and green). The red region shows the oscillators that are synchronized.

Patterns on a Triangular Lattice: On a triangular lattice, the optimallocal pattern of nearestneighbor oscillators locking into a p phase difference is incompatible with thegeometry. Results from ournumerical simulations show that the frequencies synchronize beyond a criticalcoupling strength (for a finite size system), however the phase pattern remainsdisordered with domains of different helicities frozen in (Fig. 1). Experimentsperformed in the Fraden group show clear patterns of domains with differenthelicities, and the structures are qualitatively similar to our numerical results. The experimental system, however, has acoupling that is quite different from the sine function in the Kuramotomodel. Specifically, the effectivecoupling, measured experimentally, is always positive and can be modeled as aquadratic function cutoff at 2p. We are now extending ournumerical studies to analyze domain formation in oscillators that are locallycoupled through this purely positive interaction.
Lattice Models of ReactionDiffusion Systems: Theoretical studies of chemical oscillations broadlyfall into two categories: (1) a rate equation approach in which the reactantsare assumed to be well mixed, and (2) a partial differential equation approachin which spatial dependence is taken into account by adding diffusion terms tothe rate equations. These two approaches qualitatively explain experimentalmost of the phenomena observed in the microfluidic droplets. It is, however, known fromtheoretical studies of individual level models of reaction diffusion systems(Butler and Goldenfeld, Phys. Rev. E 80, 030902 (R) (2009)) in the context ofpopulation dynamics, that the intrinsic noise in these systems can change thephase diagram, and enhance pattern formation. The parameter that controls the fluctuations is the volumeof a wellmixed patch. This volumecan be estimated as, where D is a typical diffusion constant, R is a typicalreaction rate and d is the spatial dimension of the system. This parameter can be controlled in theBZ microfluidic array tuning the concentration of the reactants. We have initiated a study of aonedimensional lattice in which each site has chemical reactions correspondingto a single Brusselator, and the activator and inhibitor can hop from one latticesite to another with specified diffusion rates.
Fig 2: (Left) Power Spectrum in the meanfield fixedpoint regime. The clear peak at finite frequency indicates coherent oscillations (Right) (a) Limit cycle corresponding to the same parameters as the left figure.

Fig 3: Phase diagram of the meanfield 1D Brusselator model. The vertical axis represents the ratio of the diffusion constants of the inhibitor and the activator species and the horizontal axis represents the creation rate of the activator. 
Results from a 1D lattice of stochastic Brusselators: Fig. 3 shows the meanfield phase diagram of aonedimensional Brusselator model with the inhibitory species represented by y and the activator species represented by x. Therate R_{0} in Fig. 3 represents the creation rate of theactivator. We studied thestochastic model over a range of parameters. Fig. 4 shows spacetime plots at R_{0 }= 3000, and D_{y}/D_{x}varying from 1 to 6. There is a clear signature of emergenceof Turing patterns in the meanfield turingoscillatory regime. Even at D_{y}/D_{x} =1, the spacetime plot shows clear deviations from a pure oscillatorybehavior. We are in the process ofcarefully analyzing the pattern formation in the onedimensional stochasticmodel. In addition, we haveinitiated studies of an inhomogeneous lattice that is meant to model the BZmicrofluidic array. In this lattersystem, the diffusion constants are different in different parts of the latticethat represent the oil and water regions. Preliminary results show an even richer phase diagram of the stochasticmodel in this spatially structured medium.
Field Theory of the Stochastic Brusslelator: Stochastic models can berepresented by field theories using operators that correspond to the processes(reaction and diffusion) occurring in these models (U. C. Tauber et al, J.Phys. A 38, R79 (2005)). We have constructed a field theory for the stochasticBrusselator model, and in collaboration with Rajesh Ravindran from theMathematical Sciences Institute in India, we are in the process of analyzingthis field theory. One of our aimsis to characterize the phase transitions in these nonequilibrium systems.