Ergodicity of the Hybridization-Expansion Monte Carlo Algorithm for   Broken-Symmetry States
  With the success of dynamical mean field theories, solvers for quantum-impurity problems have become an important tool for the numerical study of strongly correlated systems. Continuous-time Quantum Monte Carlo sampling of the expansion in powers of the hybridization between the ``impurity'' and the bath provides a powerful solver when interactions are strong. Here we show that the usual updates that add or remove a pair of creation-annihilation operators are rigorously not ergodic for several classes of broken-symmetries that involve spatial components. We show that updates with larger numbers of simultaneous updates of pairs of creation-annihilation operators remedy this problem. As an example, we apply the four operator updates that are necessary for ergodicity to the case of d-wave superconductivity in plaquette dynamical mean-field theory for the one-band Hubbard model. While the results are qualitatively similar to those previously published, they are quantitatively better that previous ones, being closer to those obtained by other approaches. 
P. S\'emon, G. Sordi, A.-M.S. Tremblay
cond-mat.str-el cond-mat.supr-con