In this paper, we present the numerical evidence of multiple bifurcations in thermoacoustic instability as we decrease the equivalence ratio (ϕ) from fuel-rich (ϕ = 1.5) to fuel-lean (ϕ =0.59) within a single-element coaxial liquid rocket engine burning hydrogen peroxide and methane. An axisymmetric compressible Large Eddy Simulation solver in an open-source platform OpenFOAM combined with a global methane chemical kinetical mechanism was used to model the turbulent mixing and combustion process. After the validation of grid independence and simulation results, the synchronization phenomenon between hydrodynamic instability, acoustics and combustion is the main focus of this paper. The results here first revealed that in a direction of decreasing equivalence ratios, the first transition to limit cycle state from combustion noise state was observed via intermittency. A second bifurcation was also observed, leading to quasi-period thermoacoustic instability from the limit cycle state. Intrinsic vorticity modes and natural acoustic modes could be regarded as oscillators. The competition and cooperation between the two modes led to a phenomenon beyond limit cycle states. The vortex dynamics of thermoacoustic systems are dominated by the interaction of intrinsic vorticity modes and natural acoustic modes and exhibit an adjustment of vortex shedding frequencies, resulting in different stability regions. Herein, the flame index was performed in more detail to distinguish between the premixed flame and diffusion flame. Results indicated that the burning mode was an important driver in promoting instability whilst it also served as a medium for competition and cooperation between the natural acoustic modes and intrinsic vorticity modes.