We investigate the low-temperature properties of the spin-1 antiferromagnetic Heisenberg chain with bond alternation by the quantum Monte Carlo method (loop algorithm). The strength of bond alternation at the gapless point is estimated as ${\ensuremath{\delta}}_{c}=0.2595\ifmmode\pm\else\textpm\fi{}0.0005.$ We confirm numerically that the low-temperature properties at the gapless point are consistent with field-theoretical predictions. The numerical results are compared with those of the spin-
1/2 antiferromagnetic Heisenberg chain and recent experimental results for $[{{\mathrm{Ni}(333\ensuremath{-}\mathrm{t}\mathrm{e}\mathrm{t})(\mathrm{\ensuremath{\mu}}\ensuremath{-}\mathrm{N}}_{3}){}}_{n}]({\mathrm{ClO}}_{4}{)}_{n}$ [$333\ensuremath{-}\mathrm{t}\mathrm{e}\mathrm{t}=\mathrm{t}\mathrm{e}\mathrm{t}\mathrm{r}\mathrm{a}\mathrm{a}\mathrm{m}\mathrm{i}\mathrm{n}\mathrm{e}$ ${N,N}^{\ensuremath{'}}$-bis(3-aminopropyl)-1,3-propanediamine].
We study a class of integrable alternating ( S 1 , S 2 ) quantum spin chains with critical ground state properties. Our main result is the description of the thermal Drude weight of the one-dimensional alternating spin chain as a function of temperature. We have identified the thermal current of the model with alternating spins as one of the conserved currents underlying the integrability. This allows for the derivation of a finite set of non-linear integral equations for the thermal conductivity. Numerical solutions to the integral equations are presented for specific cases of the spins S 1 and S 2 . In the low-temperature limit a universal picture evolves where the thermal Drude weight is proportional to temperature T and central charge c .
Affleck has made predictions for the critical exponents of a class of integrable spin chains of general spin
s. The expressions are
s dependent. Our interest here is a finite-size scaling analysis of spin-1 systems to determine the validity of the prediction for the alternation critical exponent. We discuss these results in the light of our numerical results for the conformal anomaly
c, the exponent $\ensuremath{\eta}$, and the mass-gap exponent associated with biquadratic exchange. The numerical results are not generally in good agreement with theory, and hence the possible presence of logarithmic corrections is discussed.