We consider S=1/2 antiferromagnetic Heisenberg chains with alternating bonds and quenched disorder, which represents a theoretical model of the compound CuCl_{2x}Br_{2(1-x)}(γ-{pic})_2. Using a numerical implementation of the strong disorder renormalization group method we study the low-energy properties of the system as a function of the concentration, x, and the type of correlations in the disorder. For perfect correlation of disorder the system is in the random dimer (Griffiths) phase having a concentration dependent dynamical exponent. For weak or vanishing disorder correlations the system is in the random singlet phase, in which the dynamical exponent is formally infinity. We discuss consequences of our results for the experimentally measured low-temperature susceptibility of CuCl_{2x}Br_{2(1-x)}(γ-{pic})_2.
The S = 1/2 Heisenberg chains with bond alternation and randomness on the strong bonds are studied by the density matrix renormalization group method. It is assumed that the odd-th bond is antiferromagnetic with strength J and the even-th bond can take the values JA and JF (JA > J > 0 > JF) randomly. The ground state of this model interpolates between the Haldane and dimer phases via a randomness dominated intermediate phase. On the basis of the scaling of the low energy spectrum and mean field treatment of the interchain coupling, we find that the magnetic long range order is induced by randomness in the intermediate regime. In the magnetization curves, there appears a plateau at the fractional value of the saturated magnetization. The fine structures of the magnetization curves and low energy spectrum are understood using the cluster picture. The relation to the recent experiment for (CH3)2CHNH3Cu(ClxBr1-x)3 is discussed.
The S = 1/2 Heisenberg chains with bond alternation and randomness on the strong bonds are studied by the density matrix renormalization group method. It is assumed that the oddth bond is antiferromagnetic with strength J and even-th bond can take the values JA and JF (JA > J > 0 > JF) randomly. The ground state of this model interpolates between the Haldane and dimer phases via a randomness dominated intermediate phase. Based on the scaling of the low energy spectrum and mean field treatment of the interchain coupling, it is found that the magnetic long range order is induced by randomness in the intermediate regime. In the magnetization curves, there appears a plateau at the fractional value of thesaturated magnetization. The fine structures of the magnetization curves and low energy spectrum are understood based on the cluster picture. The relation with the recent experiment for (CH3)2CHNH3Cu(ClxBr1−x)3 is discussed.
Thermodynamic properties of a tetrameric bond-alternating Heisenberg spin chain with ferromagnetic-ferromagnetic-antiferromagnetic-antiferromagnetic exchange interactions are studied using the transfer-matrix renormalization group and compared to experimental measurements. The temperature dependence of the uniform susceptibility exhibits typical ferrimagnetic features. Both the uniform and staggered magnetic susceptibilities diverge in the limit $T\ensuremath{\rightarrow}0$, indicating that the ground state has both ferromagnetic and antiferromagnetic long-range orders. A double-peak structure appears in the temperature dependence of the specific heat. Our numerical calculation gives a good account for the temperature and field dependence of the susceptibility, the magnetization, and the specific heat for $\mathrm{Cu}{(3\text{\ensuremath{-}}\mathrm{Clpy})}_{2}{({\mathrm{N}}_{3})}_{2}$ $(3\text{\ensuremath{-}}\mathrm{Clpy}=3\text{\ensuremath{-}}\text{Chloroyridine})$.
The low energy properties of the spin-1/2 random Heisenberg chain with ferromagnetic and antiferromagnetic interactions are studied by means of the density matrix renormalization group (DMRG) and real space renormalization group (RSRG) method for finite chains. The results of the two methods are consistent with each other. The deviation of the gap distribution from that of the random singlet phase and the formation of the large-spin state is observed even for relatively small systems. For a small fraction of the ferromagnetic bond, the effect of the crossover to the random singlet phase on the low temperature susceptibility and specific heat is discussed. The crossover concentration of the ferromagnetic bond is estimated from the numerical data.
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.