  We consider the performance of Glauber dynamics for the random cluster model
with real parameter $q>1$ and temperature $\beta>0$. Recent work by Helmuth,
Jenssen and Perkins detailed the ordered/disordered transition of the model on
random $\Delta$-regular graphs for all sufficiently large $q$ and obtained an
efficient sampling algorithm for all temperatures $\beta$ using cluster
expansion methods. Despite this major progress, the performance of natural
Markov chains, including Glauber dynamics, is not yet well understood on the
random regular graph, partly because of the non-local nature of the model
(especially at low temperatures) and partly because of severe bottleneck
phenomena that emerge in a window around the ordered/disordered transition.
  Nevertheless, it is widely conjectured that the bottleneck phenomena that
impede mixing from worst-case starting configurations can be avoided by
initialising the chain more judiciously. Our main result establishes this
conjecture for all sufficiently large $q$ (with respect to $\Delta$).
Specifically, we consider the mixing time of Glauber dynamics initialised from
the two extreme configurations, the all-in and all-out, and obtain a pair of
fast mixing bounds which cover all temperatures $\beta$, including in
particular the bottleneck window. Our result is inspired by the recent approach
of Gheissari and Sinclair for the Ising model who obtained a similar-flavoured
mixing-time bound on the random regular graph for sufficiently low
temperatures. To cover all temperatures in the RC model, we refine
appropriately the structural results of Helmuth, Jenssen and Perkins about the
ordered/disordered transition and show spatial mixing properties "within the
phase", which are then related to the evolution of the chain.
