Multiplicity of Singular Synchronous States in the Kuramoto Model of Coupled Oscillators

We study the Kuramoto model of globally coupled oscillators with a biharmonic coupling function. We develop an analytic self-consistency approach to find stationary synchronous states in the thermodynamic limit and demonstrate that there is a huge multiplicity of such states, which differ microscopically in the distributions of locked phases. These synchronous regimes already exist prior to the linear instability transition of the fully asynchronous state. In the presence of white Gaussian noise, the multiplicity is lifted, but the dependence of the order parameters on coupling constants remains nontrivial.