  Topological phases exhibit unconventional order that cannot be detected by
any local order parameter. In the framework of Projected Entangled Pair
States(PEPS), topological order is characterized by an entanglement symmetry of
the local tensor which describes the model. This symmetry can take the form of
a tensor product of group representations, or in the more general case a
correlated symmetry action in the form of a Matrix Product Operator(MPO), which
encompasses all string-net models. Among other things, these entanglement
symmetries allow for the description of ground states and anyon excitations.
Recently, the idea has been put forward to use those symmetries and the anyonic
objects they describe as order parameters for probing topological phase
transitions, and the applicability of this idea has been demonstrated for
Abelian groups. In this paper, we extend this construction to the domain of
non-Abelian models with MPO symmetries, and use it to study the breakdown of
topological order in the double Fibonacci (DFib) string-net and its Galois
conjugate, the non-hermitian double Yang-Lee (DYL) string-net. We start by
showing how to construct topological order parameters for condensation and
deconfinement of anyons using the MPO symmetries. Subsequently, we set up
interpolations from the DFib and the DYL model to the trivial phase, and show
that these can be mapped to certain restricted solid on solid(RSOS) models,
which are equivalent to the $((5\pm\sqrt{5})/2)$-state Potts model,
respectively. The known exact solutions of the statistical models allow us to
locate the critical points, and to predict the critical exponents for the order
parameters. We complement this by numerical study of the phase transitions,
which fully confirms our theoretical predictions; remarkably, we find that both
models exhibit a duality between the order parameters for condensation and
deconfinement.
