Vladimir Vlasov, Maxim Komarov, Arkadij Pikovskij

• We describe synchronization transitions in an ensemble of globally coupled phase oscillators with a bi-harmonic coupling function, and two sources of disorder-diversity of the intrinsic oscillators' frequencies, and external independent noise forces. Based on the self-consistent formulation, we derive analytic solutions for different synchronous states. We report on various non-trivial transitions from incoherence to synchrony, with the following possible scenarios: simple supercritical transition (similar to classical Kuramoto model); subcritical transition with large area of bistability of incoherent and synchronous solutions; appearance of a symmetric two-cluster solution which can coexist with the regular synchronous state. We show that the interplay between relatively small white noise and finite-size fluctuations can lead to metastability of the asynchronous solution.