# Kinetic theory of coupled oscillators

Last week, I gave a physics colloquium at the Catholic University of America about recent work on using kinetic theory and field theory approaches to analyze finite-size corrections to networks of coupled oscillators.  My slides are here although they are converted from Keynote so the movies don’t work.   Coupled oscillators arise in  contexts as diverse as the brain, synchronized flashing of fireflies, coupled Josephson junctions, or unstable modes of the Millennium bridge in London.  Steve Strogatz’s book Sync gives a popular account of the field.  My talk considers the Kuramoto model

$\dot{\theta}_i = \omega_i+\frac{K}{N}\sum_j\sin(\theta_j-\theta_i)$   (1)

where the frequencies $\omega$ are drawn from a fixed distribution $g(\omega)$. The model describes the dynamics of the phases $\theta$ of an all-to-all connected network of oscillators.  It can be considered to be the weak coupling limit of a set of nonlinear oscillators with different natural frequencies and a synchronizing phase response curve.

Kuramoto in the 1970’s showed in the mean field limit $N\rightarrow \infty$, that when the coupling strength $K$ is below a critical value $K_c$, the osciilators will be in an incoherent state but for $K>K_c$ a group of oscillators will phase lock. This transition could be observed with the order parameter

$Z=\frac{1}{N}\sum_j e^{i\theta_j}= r e^{i\psi}$

The bifurcation is reminiscent of a magnetization phase transition in equilibrium statistical mechanics.  In 1991, Steve Strogatz and Rennie Mirollo showed that the incoherent state loses stability at the critical point as expected but it is marginally stable below criticality.  This was surprising because simulations of the Kuramoto model showed that the magnitude of the order parameter decayed to zero below the critical point quite quickly.  This was partially resolved a year later by Strogatz, Mirollo and Matthews, who showed that the order parameter decays due to dephasing, an effect similar to Landau damping in plasma physics.

However, open questions about how the system behaved for a large but finite number of oscillators remained.  Did finite size corrections stabilize the marginal incoherent state and not just the order parameter?  Also for finite size, the phase transition is no longer sharp but smoothed over by fluctuations in the order parameter.  Hence, a calculation of the variance of the order parameter was necessary.  How could this be calculated? In two papers with Eric Hildebrand and Michael Buice (Hildebrand et al. PRL 98 , 054101 (2007) and Buice and Chow, PRE 76 031118 (2007), both available here) we showed that an application of methods from the kinetic theory of plasmas and field theory can answer these questions.

Kinetic theory was developed by Boltzmann and Maxwell to show how microscopic Hamiltonian dynamics of particles could account for the thermodynamic properties of gases.  We adapted a formalism previously applied to plasma physics.  The idea starts by encoding the microscopic dynamics of the oscillators into a single probability density function $\eta(\theta,\omega)=\sum_{ij}\delta(\theta-\theta_i)\delta(\omega-\omega_j)$.  The density will then obey a continuity equation in oscillator phase space of the form

$\frac{\partial \eta}{\partial t}+\frac{\partial}{\partial\theta}[v(\theta,\omega)\eta]=0$

where the velocity $v=\dot\theta_i$ is given by the oscillator equation (1).  The continuity equation is called the Klimontovich equation in plasma physics and is only defined in the weak sense of distributions.  It is entirely equivalent to the original oscillator system and is only made useful for calculations by taking averages.  However, the Klimontovich equation is nonlinear in $\eta$ so when you take the average you’ll find that the equation for the first moment $\langle \eta(\theta,\omega,t) \rangle$ will depend on the second moment $\langle \eta(\theta,\omega,t)\eta(\theta',\omega',t')\rangle$.  The second moment will in turn depend on the third and so on.  This is known as the BBGKY hierarchy.

The hierarchy is made tractable by truncating at some order with the idea that the higher  order moments (actually cumulants) will matter less and less as the number of oscillators increases.  The reason is that the cumulants represent direct oscillator interactions and the effect of a single oscillator goes to zero as the number of oscillators goes to infinity.  If two oscillator correlations are ignored then we end up with mean field theory, where there are no correlations or fluctuations.  We showed in Hildebrand et al. that the correlations could be computed by going to the next level in the hierarchy.  One thing that is important to note is that the mean field theory has the exact same form as the Klimontovich equation.  The only difference being that the mean field equation must have smooth solutions.  Almost everyone goes directly to mean field theory when analyzing coupled oscillators using this formalism and this distinction is glossed over.

Now, calculations using the moment hierarchy involve solving coupled integro-differential equations, which is still pretty unwieldy.  This can be simplified by noting that to compute the cumulants of $\eta$ we can write down a generating functional.  In this case it will be an infinite dimensional path integral.  However, it can then be analyzed using the methods of field theory.  The generating functional, like any other finite dimensional integral, can be analyzed perturbatively using the methods of steepest descents and field theory is basically asymptotic analysis in infinite dimensions.  In our case, 1/N is the small parameter that allows us to use steepest descents, called a loop expansion in field theory.  We can keep track of the terms in the expansion (which multiply quite rapidly) by using Feynman diagrams.  The number of loops in the diagrams indicate the order in 1/N of the expansion.  In Buice and Chow, we show that the field theory formalism and the moment hierarchy are equivalent and a natural truncation or closure of the moment hierarchy is given by the loop expansion.  We also show that the incoherent state is strictly stable for finite N.

There is one sad story that is related to this project.  Eric was actually a grad student of J.D. Crawford at the University of Pittsburgh where I was on the math faculty from 1998 to 2004.  I had known about J.D.’s work in dynamical systems for several years prior to arriving at Pitt and was looking forward to finally meeting him.  However, he contracted cancer right around that time so I would miss him when I visited Pitt during my interview phase.  He then died the day I arrived to start my job so I never got to meet him.  He was only 44 and in his prime.  Eric eventually finished his PhD with me working on this project.

## 12 thoughts on “Kinetic theory of coupled oscillators”

1. […] By Carson Chow I was visiting the University of Utah this past week.  I gave talks on the Kinetic Theory of Coupled Oscillators and on Deriving Moment Equations for Neural Networks. On my way to the airport I wondered what […]

Like

2. […] Here is the longer and more informative version of the talk.  I summarized the papers on this talk here. Possibly related posts: (automatically generated)Revised version of paperFreiburg! im Breisgau!Off […]

Like

3. […] and I also applied the density functional idea to studying finite size effects in the Kuramoto model.  Eric Hildebrand and I had derived a moment hierarchy to describe fluctuations and correlations […]

Like

4. […] really get across why a kinetic theory of coupled oscillators is interesting and useful.   Here is the post giving more of the backstory on the topic, which has a link to an older version of the […]

Like

5. […] I gave a talk on finite system size effects in neural networks.  I mostly covered the work on the Kuramoto model with a little bit on synaptically coupled phase neurons at the end.  My slides are […]

Like

6. […] and Applied Mechanics to give a talk in one of the sessions.  I gave a condensed version of the kinetic theory of coupled oscillators talk I gave in Warwick last month.  The theme of the session was on recent advances in nonlinear […]

Like

7. […] adapted the formalism we used to analyze the finite size effects of the Kuramoto system (see here for summary) to a system of synaptically coupled phase […]

Like

8. […] we perform the calculation explicitly for the Kuramoto model of coupled oscillators (e.g. see here) but it can be generalized to any network of coupled elements.  The calculation relies on the path […]

Like

9. […] theory and is similar to what we used in our previous work on the Kuramoto model (see here and  here) and the “Spike model” (see here).  Our heuristic paper on path integral methods is […]

Like

10. […] spiking neurons. I’ve written about this work multiple times (e.g. here,  here, and here). In this series of papers, we looked at how you can compute fluctuations around the infinite […]

Like

11. […] size effects in a network of coupled deterministic oscillators (e.g. see here, here, here and here).  The second was to derive a set of generalized Wilson-Cowan equations that includes correlation […]

Like