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. These systems are modelled 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.