A generalization to the quantum case of a recently introduced algorithm (Y. Tomita and Y. Okabe, Phys. Rev. Lett. {\bf 86}, 572 (2001)) for the determination of the critical temperature of classical spin models is proposed. We describe a simple method to automatically locate critical points in (Quantum) Monte Carlo simulations. The algorithm assumes the existence of a finite correlation length in at least one of the two phases surrounding the quantum critical point. We illustrate these ideas on the example of the critical inter-chain coupling for which coupled antiferromagnetic S=1 spin chains order at T=0. Finite-size scaling relations are used to determine the exponents, $\nu=0.72(2)$ and $\eta=0.038(3)$ in agreement with previous estimates.

Open segments of S=1 antiferromagnetic spin chains are studied at finite temperatures and fields using continuous time Quantum Monte Carlo techniques. By calculating the resulting magnetization profiles for a large range of chain lengths with fixed field and temperature we reconstruct the experimentally measured NMR spectrum of impurity doped Y$_2$BaNi$_{1-x}$Mg$_x$O$_5$. For temperatures above the gap the calculated NMR spectra are in excellent agreement with the experimental results, confirming the existence of $S=1/2$ excitations at the end of open S=1 chain segments. At temperatures below the gap, neglecting inter chain couplings, we still find well defined peaks in the calculated NMR spectra corresponding to the $S=1/2$ chain end excitations. At low temperatures, inter chain couplings could be important, resulting in a more complicated phase.

We investigate the occurrence of singlet bound states in the Raman spectrum of dimerized spin 1/2 chains by Exact Diagonalization and Density Matrix Renormalization Group techniques. We predict that several bulk $\delta$-peaks could be observed in pure systems. Furthermore, we show that new low energy lines arise from non-magnetic impurity doping. These features are interpreted in terms of soliton-antisoliton and soliton-impurity bound states respectively. Energies and spectral weights associated with these bound states are sensitive to lattice relaxation effects. Our results are discussed in the context of the inorganic spin-Peierls compound CuGeO$_3$ and quantitatively compared to recent Raman experiments.

The role of the spin-phonon coupling in spin-Peierls chains doped with spin-0 or spin-1 impurities is investigated by various numerical methods such as exact diagonalization, quantum Monte Carlo simulations and Density Matrix Renormalization Group. Various treatments of the lattice, in a fully quantum mechanical way, classically in the adiabatic approximation or using a fixed three-dimensional dimerization pattern are compared. For an isolated chain, strong bonds form between the two spin-1/2 sites next to the impurity site, leading to the appearance of magneto-elastic solitons. We also show that these excitations do not bind to spin-0 impurities but are weakly attracted by spin-1 impurities. However, the interchain elastic coupling generates an effective confining potential at the non-magnetic impurity site which can lead to the formation of soliton-impurity bound states. We also predict that a soliton and an antisoliton bound to two impurities on the same chain can annihilate each other when the separation between the impurities is smaller than a critical value depending on the interchain elastic constant.

A simple intuitive picture of spin-Peierls antiferromagnets arises from regarding the elementary excitations as S=1/2 solitons. In a strictly one-dimensional system these excitations are assumed not to form bound-states and to be repelled by impurities. Couplings to the three-dimensional lattice are assumed to produce an effective confining potential which binds solitons to antisolitons and to impurities, with the number of bound-states increasing as the interchain coupling goes to 0. We investigate these various assumptions numerically in a phononless model where spontaneous dimerization arises from frustration and the interchain coupling is treated in mean field theory.

The role of dynamical magneto-elastic coupling in spin-Peierls chains is investigated by numerical and analytical techniques. We show that a Heisenberg spin chain coupled to dynamical optical phonons exhibits a transition towards a spontaneously dimerized state in a wide range of parameter space. The low energy excitations are characterized as solitons. No binding between solitons occurs in the isolated spin-phonon chain and the dynamical spin structure factor shows a broad magnon dispersion. However, elastic interchain coupling can lead to the formation of bound states.

The $S=1$ antiferromagnetic Heisenberg chain both with and without single ion anisotropy is studied. Using the recently proposed density matrix renormalization group technique we calculate the energy gaps as well as several different correlation functions. The two gaps, $\Delta_{||}, \Delta_\perp$, along with associated correlation lengths and velocities are determined. The numerical results are shown to be in good agreement with theoretical predictions derived from the nonlinear sigma model and a free boson model. We also study the $S=1/2$ excitations that occur at the ends of open chains; in particular we study the behavior associated with open boundary conditions, using a model of $S=1/2$ spins coupled to the free bosons.

The structure function, S(k), for the s=1, Haldane gap antiferromagnetic chain, is measured accurately using the recent density matrix renormalization group method, with chain-length 100. Excellent agreement with the nonlinear $\sigma$ model prediction is obtained, both at $k\approx \pi$ where a single magnon process dominates and at $k\approx 0$ where a two magnon process dominates. We repeat our calculation with crystal field anisotropy chosen to model NENP, obtaining good agreement with both field theory predictions and recent experiments. Correlation lengths, gaps and velocities are determined for both polarizations.

Universal properties of the zero temperature superconductor-insulator transition in two-dimensional amorphous films are studied by extensive Monte Carlo simulations of bosons in a disordered medium. We report results for both short-range and long-range Coulomb interactions for several different points in parameter space. In all cases we observe a transition from a superconducting phase to an insulating Bose glass phase. {}From finite-size scaling of our Monte Carlo data we determine the universal conductivity $\sigma^*$ and the critical exponents at the transition. The result $\sigma^* = (0.55 \pm 0.06) (2e)^2/h$ for bosons with long-range Coulomb interaction is roughly consistent with experiments reported so far. We also find $\sigma^* = (0.14 \pm 0.03) (2e)^2/h$ for bosons with short-range interactions.