We consider a multi-phase energy functional whose ellipticity rate and growth radically change according to variable coefficients. The aim of the talk is to show the Hölder regularity for the gradients of local minimizers of the multi-phase energy functional under sharp assumptions relating the growth couples to the Hölder exponents of the modulating coefficients. The talk is based on the joint work with Cristiana De Filippis (Oxford).