Numerical study of the incommensurate phase in spin-Peierls systems

We analyze several properties of the lattice solitons in the incommensurate phase of spin-Peierls systems using exact diagonalization and quantum Monte Carlo simulations. These systems are modeled by an antiferromagnetic Heisenberg chain with nearest-and next-nearest-neighbor interactions coupled to the lattice in the adiabatic approximation. Several relations among features of the solitons and magnetic properties of the system have been determined and compared with analytical predictions. We have studied in particular the relation between the soliton width and the spin-Peierls gap. Although this relation has the form predicted by bosonized field theories, we have found some important quantitative differences which could be relevant to describe experimental studies of spin-Peierls materials.