We study the thermodynamic properties of a microscopic model of coupled oscillators that exhibits a dynamical phase transition from a desynchronized to a synchronized phase. We consider two different configurations for the thermodynamic forces applied on the oscillators, one resembling the macroscopic power grids, and one resembling autonomous molecular motors. We characterize the input and the output power as well as the efficiency at maximum power, providing analytic expressions for such quantities near the critical coupling strength. We discuss the role of the quenched disorder in the thermodynamic force distributions and show that such a disorder may lead to an enhancement of the efficiency at maximum power.OPEN ACCESS RECEIVED exclusion process on a single lattice [8,25], or on a network [26]. In these studies, we found an increase of the EMP in a many-motor system with respect to the single motor case, for a suitable choice of the model parameters. Remarkably,in [8,25] we found that the enhancement of the EMP occurs in a range of parameter values compatible with the biological estimates for the molecular motor Kinesin. From those studies we concluded that after a dynamical phase transition the dynamical response of the system to an external drive can change, leading in turn to a change in the thermodynamic properties. Specifically the dependence of the delivered power on the driving thermodynamic forces may vary.One of the main limitation that one faces when studying the thermodynamic properties of exclusion processes on a lattice, is that the intensity of interaction between the motors can only be indirectly tuned by changing the kinetic parameters and thus the density of motors on the lattice [8,25,26]. Furthermore, if one wants to study the effect of force disorder on the motor particles, one has to resort to numerical simulations, as no exact result exists for the exclusion process with heterogeneous particles.Instead, here we consider a model of N interacting microscopic particles, where the particle-particle interaction is an explicit parameter, that can be tuned in order to drive a dynamical phase transition, from a weakly interacting-incoherent system to a strongly interacting-coherent system. This model was originally introduced by Sakaguchi in [27] as an extension of the Kuramoto model (KM) [28], to study the synchronization of a group of interacting oscillators in contact with a reservoir at constant temperature [29,30]. Furthermore, the effect of quenched disorder in the thermodynamic force distribution can be taken into account within the present model. The model is introduced and discussed in section 2. The N interacting particles can be viewed as a network of energy producers and users or as a system of interacting autonomous motors under the effect of thermodynamic forces. Since the dynamical phase diagram can be obtained in terms of the particle interaction strength, temperature and force distribution, in section 3 we will discuss how to calculate the relevant thermodynamic quantities, name...