{"title":"On the Numerical Integration of the Multidimensional Kuramoto Model","authors":"Marcus A. M. de Aguiar","doi":"10.1007/s13538-024-01493-z","DOIUrl":null,"url":null,"abstract":"<p>The Kuramoto model, describing the synchronization dynamics of coupled oscillators, has been generalized in many ways over the past years. One recent extension of the model replaces the oscillators, originally characterized by a single phase, by particles with <span>\\(D-1\\)</span> internal phases, represented by a point on the surface of the unit D-sphere. Particles are then more easily represented by <i>D</i>-dimensional unit vectors than by <span>\\(D-1\\)</span> spherical angles. However, numerical integration of the state equations should ensure that the propagated vectors remain unit and that particles rotate on the sphere as predicted by the dynamical equations. As discussed in (Lee et al. in Journal of Statistical Mechanics: Theory and Experiment 2023(4):043403, 2023), integration of the three-dimensional Kuramoto model using Euler’s method with time step <span>\\(\\Delta t\\)</span> not only changes the norm of the vectors but produces a small rotation of the particles around the wrong axis. Importantly, the error in the axis’ direction does not vanish in the limit <span>\\(\\Delta t \\rightarrow 0\\)</span>. Therefore, instead of displacing the unit vectors in the direction of the velocity, one should perform a sequence of direct small rotations, as dictated by the equations of motion. This keeps the particles on the sphere at all times, ensuring exact norm preservation, and rotates the particles around the proper axis for small <span>\\(\\Delta t\\)</span> (Lee et al. in Journal of Statistical Mechanics: Theory and Experiment 2023(4):043403, 2023). Here, I propose an alternative way to do such integration by rotations in 3D that can be generalized to more dimensions using Cayley-Hamilton’s theorem. Explicit formulas are provided for 2, 3, and 4 dimensions. I also compare the results with the fourth-order Runge–Kutta method, which seems to provide accurate results even requiring renormalization of the vectors after each integration step.</p>","PeriodicalId":499,"journal":{"name":"Brazilian Journal of Physics","volume":null,"pages":null},"PeriodicalIF":1.5000,"publicationDate":"2024-05-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Brazilian Journal of Physics","FirstCategoryId":"101","ListUrlMain":"https://doi.org/10.1007/s13538-024-01493-z","RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"PHYSICS, MULTIDISCIPLINARY","Score":null,"Total":0}
引用次数: 0
Abstract
The Kuramoto model, describing the synchronization dynamics of coupled oscillators, has been generalized in many ways over the past years. One recent extension of the model replaces the oscillators, originally characterized by a single phase, by particles with \(D-1\) internal phases, represented by a point on the surface of the unit D-sphere. Particles are then more easily represented by D-dimensional unit vectors than by \(D-1\) spherical angles. However, numerical integration of the state equations should ensure that the propagated vectors remain unit and that particles rotate on the sphere as predicted by the dynamical equations. As discussed in (Lee et al. in Journal of Statistical Mechanics: Theory and Experiment 2023(4):043403, 2023), integration of the three-dimensional Kuramoto model using Euler’s method with time step \(\Delta t\) not only changes the norm of the vectors but produces a small rotation of the particles around the wrong axis. Importantly, the error in the axis’ direction does not vanish in the limit \(\Delta t \rightarrow 0\). Therefore, instead of displacing the unit vectors in the direction of the velocity, one should perform a sequence of direct small rotations, as dictated by the equations of motion. This keeps the particles on the sphere at all times, ensuring exact norm preservation, and rotates the particles around the proper axis for small \(\Delta t\) (Lee et al. in Journal of Statistical Mechanics: Theory and Experiment 2023(4):043403, 2023). Here, I propose an alternative way to do such integration by rotations in 3D that can be generalized to more dimensions using Cayley-Hamilton’s theorem. Explicit formulas are provided for 2, 3, and 4 dimensions. I also compare the results with the fourth-order Runge–Kutta method, which seems to provide accurate results even requiring renormalization of the vectors after each integration step.
期刊介绍:
The Brazilian Journal of Physics is a peer-reviewed international journal published by the Brazilian Physical Society (SBF). The journal publishes new and original research results from all areas of physics, obtained in Brazil and from anywhere else in the world. Contents include theoretical, practical and experimental papers as well as high-quality review papers. Submissions should follow the generally accepted structure for journal articles with basic elements: title, abstract, introduction, results, conclusions, and references.