Finite temperature fermion Monte Carlo simulations of frustrated spin-Peierls systems

This approach provides a manifold of equivalent actions over which the negative sign problem can be optimised. As a result, we can reach temperature scales well below the magnetic scale.

Here, we show how to generalise this algorithm to spin-Peierls systems. In contrast to exact diagonalisation approaches, Monte Carlo methods are not Hilbert space bound such that the computational effort per sweep remains invariant when adding phonons. However, the computational effort required to generate Independent configurations increases in the presence of phonons.