The low-temperature thermodynamic properties of a spin-\textonehalf{} one-dimensional random anisotropic (Heisenberg-Ising) antiferromagnet described by the Hamiltonian $H=\ensuremath{\Sigma}{i}^{}{J}_{i}({\ensuremath{\sigma}}_{x}^{i}{\ensuremath{\sigma}}_{x}^{i+1}+{\ensuremath{\sigma}}_{y}^{i}{\ensuremath{\sigma}}_{y}^{i+1}+\ensuremath{\gamma}{\ensuremath{\sigma}}_{z}^{i}{\ensuremath{\sigma}}_{z}^{i+1})$ are studied as a function of disorder and anisotropy. The
Ji> 0 are independent random variables obeying a probability distribution
P(J), and $0<~\ensuremath{\gamma}<~\ensuremath{\infty}$. The approach used is a numerical implementation of a real-space renormalization-group (RG) method previously introduced. The isotropic Heisenberg case ($\ensuremath{\gamma}=1$), the
XY case ($\ensuremath{\gamma}=0$). and the Ising case ($\ensuremath{\gamma}=\ensuremath{\infty}$) are fixed points of the RG transformation. It is found that in the
XY region $\ensuremath{\gamma}<~1$, including the Heisenberg point, the system exhibits singular behavior in the thermodynamic properties for arbitrary probability distributions for the couplings. For the
XY limit ($\ensuremath{\gamma}=0$) this is in agreement with known exact results. The functional form for the low-temperature susceptibility is found to be $\ensuremath{\chi}\ensuremath{\sim}\frac{1}{(T{\mathrm{ln}}^{m}(\frac{T}{{T}_{0}}))}$ in the entire region $0<~\ensuremath{\gamma}<~1$ for arbitrary probability distributions. In the Ising region ($\ensuremath{\gamma}>1$) the susceptibility shows an approximate power-law divergence for small anisotropy but goes eventually to zero as $T\ensuremath{\rightarrow}0$. Possible relevance of these results to recent experiments on Qn
(TCNQ)2 is discussed.
The low-temperature thermodynamic properties of a random antiferromagnetic Heisenberg chain in zero magnetic field are studied approximately by a renormalisation group method. The couplings J obey a probability distribution function P(J) approximately 1/J alpha , with alpha a constant. It is found that the susceptibility chi and the specific heat C behave like chi approximately T - alpha chi and C approximately T 1- alpha c at low temperatures alpha chi and alpha c are weakly varying functions of temperature and slowly approach unity as T goes to zero. Comparison of these results with other theoretical analysis as well as experimental results on organic charge transfer compounds is also discussed.
The magnetic susceptibility ${\ensuremath{\chi}}^{*}(t)$ and specific heat
C(t) versus temperature t of the spin
S=1/2 antiferromagnetic (AF) alternating-exchange
(J1 and
J2) Heisenberg chain are studied for the entire range $0<~\ensuremath{\alpha}<~1$ of the alternation parameter $\ensuremath{\alpha}\ensuremath{\equiv}{J}_{2}{/J}_{1}{(J}_{1},$
J2> 0, J2< J1, t=kBT/J1, ${\ensuremath{\chi}}^{*}=\ensuremath{\chi}{J}_{1}{/Ng}^{2}{\ensuremath{\mu}}_{B}^{2}).$ For the uniform chain $(\ensuremath{\alpha}=1),$ the high-accuracy ${\ensuremath{\chi}}^{*}(t)$ and
C(t) Bethe ansatz data of Kl\"umper and Johnston (unpublished) are shown to agree very well at low t with the respective exact theoretical low-
t logarithmic correction predictions of Lukyanov [Nucl. Phys. B
522, 533 (1998)]. Accurate $(\ensuremath{\sim}{10}^{\ensuremath{-}7})$ independent empirical fits to the respective data are obtained over t ranges spanning 25 orders of magnitude, $5\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}25}<~t<~5,$ which contain extrapolations to the respective exact
t=0 limits. The infinite temperature entropy calculated using our
C(t) fit function is within 8 parts in
108 of the exact value
ln2. Quantum Monte Carlo (QMC) simulations and transfer-matrix density-matrix renormalization group (TMRG) calculations of ${\ensuremath{\chi}}^{*}(\ensuremath{\alpha},t)$ are presented for
0.002< t< 10 and $0.05<~\ensuremath{\alpha}<~1,$ and an accurate $(2\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}4})$ two-dimensional $(\ensuremath{\alpha},t)$ fit to the combined data is obtained for
0.01< t< 10 and $0<~\ensuremath{\alpha}<~1.$ From the low-
t TMRG data, the spin gap $\ensuremath{\Delta}(\ensuremath{\alpha})$ is extracted for $0.8<~\ensuremath{\alpha}<~0.995$ and compared with previous results, and a fit function is formulated for $0<~\ensuremath{\alpha}<~1$ by combining these data with literature data. We infer from our data that the asymptotic critical regime near the uniform chain limit is only entered for $\ensuremath{\alpha}\ensuremath{\gtrsim}0.99.$ We examine in detail the theoretical predictions of Bulaevskii [Sov. Phys. Solid State
11, 921 (1969)], for ${\ensuremath{\chi}}^{*}(\ensuremath{\alpha},t)$ and compare them with our results. To illustrate the application and utility of our theoretical results, we model our experimental $\ensuremath{\chi}(T)$ and specific heat
Cp(T) data for
NaV2O5 single crystals in detail. The $\ensuremath{\chi}(T)$ data above the spin dimerization temperature ${T}_{\mathrm{c}}\ensuremath{\approx}34\mathrm{K}$ are not in quantitative agreement with the prediction for the
S=1/2 uniform Heisenberg chain, but can be explained if there is a moderate ferromagnetic interchain coupling and/or if J changes with T. Fitting the $\ensuremath{\chi}(T)$ data using our ${\ensuremath{\chi}}^{*}(\ensuremath{\alpha},t)$ fit function, we obtain the sample-dependent spin gap and range $\ensuremath{\Delta}{(T=0)/k}_{\mathrm{B}}=103(2)\mathrm{K},$ alternation parameter $\ensuremath{\delta}(0)\ensuremath{\equiv}(1\ensuremath{-}\ensuremath{\alpha})/(1+\ensuremath{\alpha})=0.034(6)$ and average exchange constant
J(0)/kB=640(80)K. The $\ensuremath{\delta}(T)$ and $\ensuremath{\Delta}(T)$ are derived from the data. A spin pseudogap with magnitude $\ensuremath{\approx}0.4\ensuremath{\Delta}(0)$ is consistently found just above
Tc, which decreases with increasing temperature. From our
Cp(T) measurements on two crystals, we infer that the magnetic specific heat at low temperatures $T\ensuremath{\lesssim}15\mathrm{K}$ is too small to be resolved experimentally, and that the spin entropy at
Tc is too small to account for the entropy of the transition. A quantitative analysis indicates that at
Tc, at least 77% of the entropy change due to the transition at
Tc and associated order parameter fluctuations arise from the lattice and/or charge degrees of freedom and less than 23% from the spin degrees of freedom.
The quantum renormalization-group technique recently developed by Chakravarty and the present author is used to study weak random couplings in quantum-mechanical and classical nearest-neighbor Heisenberg ferromagnets in D=2+\ensuremath{\epsilon} dimensions and at low temperatures T. Momentum-shell recursion relations are obtained by considering the renormalization-group transformations of the cumulants of the probability distribution of the random exchange couplings. From a one-loop calculation it is found that, if the second cumulants ${\ensuremath{\sigma}}^{2}$ of the random couplings are smaller than a critical value ${\ensuremath{\sigma}}_{c}^{2}$, the only modification due to randomness in the low-temperature behavior of the order-parameter correlation length \ensuremath{\xi} in two dimensions is a renormalization of the spin stiffness ${\ensuremath{\rho}}_{s}$ according to ${\ensuremath{\rho}}_{s}$(X)=(1-X${)}^{\ensuremath{\varphi}}$${\ensuremath{\rho}}_{s}$(0), where X=${\ensuremath{\sigma}}^{2}$/${\ensuremath{\sigma}}_{c}^{2}$ and \ensuremath{\varphi}=(1/3 within the one-loop approximation. At zero temperature \ensuremath{\xi} diverges for with ${\ensuremath{\nu}}_{x}$=1/D. Precisely at X=1 the critical behavior of the system is characterized by new exponents; the dynamical critical exponent z is given by z=2+\ensuremath{\varphi}/${\ensuremath{\nu}}_{x}$. For D(1-\ensuremath{\varphi})>2 a new phase appears in the T-X phase diagram, and it is argued that this phase is characterized by spin-glass order.
The one-dimensional quantum spin-\textonehalf{} Heisenberg antiferromagnetic model with randomly distributed interaction strengths is solved approximately for several different distributions. Ground-state energy and low-temperature properties are evaluated. Universal qualitative features are found in the specific heat and the magnetic susceptibility, which display a power-law dependence on temperature. Such features hold for nonsingular distributions as well as for distributions with power-law divergence at the origin. The approximate method of solution is based on successive eliminations of spins coupled by the maximum coupling constant.