In this thesis, novel Monte Carlo methods for precisely calculating the critical phenomena of the effectively frustrated quantum spin system are developed and applied to the critical phenomena of the spin-Peierls systems. Three significant methods are introduced for the first time: a new optimization algorithm of the Markov chain transition kernel based on the geometric weight-allocation approach, the extension of the worm (directed-loop) algorithm to nonconserved particles, and the combination with the level spectroscopy. Utilizing these methods, the phase diagram of the one-dimensional XXZ...