-
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.