Finite-Field Ground State of the S=1 Antiferromagnetic-Ferromagnetic Bond-Alternating Chain

We investigate the finite-field ground state of the S=1 antiferromagnetic-ferromagnetic bond-alternating chain described by the Hamiltonian {\calH}=\sum\nolimits_{\ell}\bigl\{\vecS_{2\ell-1}\cdot\vecS_{2\ell} +J\vecS_{2\ell}\cdot\vecS_{2\ell+1}\bigr\} +D\sum\nolimits_{\ell} \bigl(S_{\ell}^z)^2 -H\textstyle\sum\nolimits_\ell S_\ell^z, where \hbox{$J\leq0$} and \hbox{$-\infty<D<\infty$}. We find that two kinds of magnetization plateaux at a half of the saturation magnetization, the 1/2-plateaux, appear in the ground-state magnetization curve; one of them is of the Haldane type and the other is of the large-$D$-type. We determine the 1/2-plateau phase diagram on the $D$ versus $J$ plane, applying the twisted-boundary-condition level spectroscopy methods developed by Kitazawa and Nomura. We also calculate the ground-state magnetization curves and the magnetization phase diagrams by means of the density-matrix renormalization-group method

Spin and chiral orderings of frustrated quantum spin chains

Ordering of frustrated S=1/2 and 1 XY and Heisenberg spin chains with the competing nearest- and next-nearest-neighbor antiferromagnetic couplings is studied by exact diagonalization and density-matrix renormalization-group methods. It is found that the S=1 XY chain exhibits both gapless and gapped `chiral' phases characterized by the spontaneous breaking of parity, in which the long-range order parameter is a chirality, $\kappa_i = S_i^xS_{i+1}^y-S_i^yS_{i+1}^x$, whereas the spin correlation decays either algebraically or exponentially. Such chiral phases are not realized in the S=1/2 XY chain nor in the Heisenberg chains.Comment: 4 pages, 5 EPS-figures, LaTeX(RevTeX),to appear in J.Phys.Soc.Japa

How to distinguish the Haldane/Large-D state and the intermediate-D state in an S=2 quantum spin chain with the XXZ and on-site anisotropies

We numerically investigate the ground-state phase diagram of an S=2 quantum spin chain with the $XXZ$ and on-site anisotropies described by ${\mathcal H}=\sum_j (S_j^x S_{j+1}^x+S_j^y S_{j+1}^y+\Delta S_j^z S_{j+1}^z) + D \sum_j (S_j^z)^2$, where $\Delta$ denotes the XXZ anisotropy parameter of the nearest-neighbor interactions and $D$ the on-site anisotropy parameter. We restrict ourselves to the $\Delta>0$ and $D>0$ case for simplicity. Our main purpose is to obtain the definite conclusion whether there exists or not the intermediate-$D$ (ID) phase, which was proposed by Oshikawa in 1992 and has been believed to be absent since the DMRG studies in the latter half of 1990's. In the phase diagram with $\Delta>0$ and $D>0$ there appear the XY state, the Haldane state, the ID state, the large-$D$ (LD) state and the N\'eel state. In the analysis of the numerical data it is important to distinguish three gapped states; the Haldane state, the ID state and the LD state. We give a physical and intuitive explanation for our level spectroscopy method how to distinguish these three phases.Comment: Proceedings of "International Conference on Frustration in Condensed Matter (ICFCM)" (Jan. 11-14, 2011, Sendai, Japan

Second order quantum renormalisation group of XXZ chain with next nearest neighbour interactions

We have extended the application of quantum renormalisation group (QRG) to the anisotropic Heisenberg model with next-nearest neighbour (n-n-n) interaction. The second order correction has to be taken into account to get a self similar renormalized Hamiltonian in the presence of n-n-n-interaction. We have obtained the phase diagram of this model which consists of three different phases, i.e, spin-fluid, dimerised and Ne'el types which merge at the tri-critical point. The anisotropy of the n-n-n-term changes the phase diagram significantly. It has a dominant role in the Ne'el-dimer phase boundary. The staggered magnetisation as an order parameter defines the border between fluid-Ne'el and Ne'el-dimer phases. The improvement of the second order RG corrections on the ground state energy of the Heisenberg model is presented. Moreover, the application of second order QRG on the spin lattice model has been discussed generally. Our scheme shows that higher order corrections lead to an effective Hamiltonian with infinite range of interactions.Comment: 10 pages, 4 figures and 1 tabl

Universal emergence of the one-third plateau in the magnetization process of frustrated quantum spin chains

We present a numerical study of the magnetization process of frustrated quantum spin-S chains with S=1, 3/2, 2 as well as the classical limit. Using the exact diagonalization and density-matrix renormalization techniques, we provide evidence that a plateau at one third of the saturation magnetization exists in the magnetization curve of frustrated spin-S chains with S>1/2. Similar to the case of S=1/2, this plateau state breaks the translational symmetry of the Hamiltonian and realizes an up-up-down pattern in the spin component parallel to the external field. Our study further shows that this plateau exists both in the cases of an isotropic exchange and in the easy-axis regime for spin-S=1, 3/2, and 2, but is absent in classical frustrated spin chains with isotropic interactions. We discuss the magnetic phase diagram of frustrated spin-1 and spin-3/2 chains as well as other emergent features of the magnetization process such as kink singularities, jumps, and even-odd effects. A quantitative comparison of the one-third plateau in the easy-axis regime between spin-1 and spin-3/2 chains on the one hand and the classical frustrated chain on the other hand indicates that the critical frustration and the phase boundaries of this state rapidly approach the classical result as the spin S increases.Comment: 15 pages RevTex4, 13 figure

Excitations with fractional spin less than 1/2 in frustrated magnetoelastic chains

We study the magnetic excitations on top of the plateaux states recently discovered in spin-Peierls systems in a magnetic field. We show by means of extensive density matrix renormalization group (DMRG) computations and an analytic approach that one single spin-flip on top of $M=1-\frac2N$ ($N=3,4,...$) plateau decays into $N$ elementary excitations each carrying a fraction $\frac1N$ of the spin. This fractionalization goes beyond the well-known decay of one magnon into two spinons taking place on top of the M=0 plateau. Concentrating on the $\frac13$ plateau (N=3) we unravel the microscopic structure of the domain walls which carry fractional spin-$\frac13$, both from theory and numerics. These excitations are shown to be noninteracting and should be observable in x-ray and nuclear magnetic resonance experiments.Comment: 6 pages, 5 figures. Accepted to be published in Phys. Rev.

Magnetic properties of the S=1/2S=1/2 distorted diamond chain at T=0

We explore, at T=0, the magnetic properties of the $S=1/2$ antiferromagnetic distorted diamond chain described by the Hamiltonian {\cal H} = \sum_{j=1}^{N/3}{J_1 ({\bi S}_{3j-1} \cdot {\bi S}_{3j} + {\bi S}_{3j} \cdot {\bi S}_{3j+1}) + J_2 {\bi S}_{3j+1} \cdot {\bi S}_{3j+2} + J_3 ({\bi S}_{3j-2} \cdot {\bi S}_{3j} + {\bi S}_{3j} \cdot {\bi S}_{3j+2})} \allowbreak - H \sum_{l=1}^{N} S_l^z with $J_1, J_2, J_3\ge0$, which well models ${\rm A_3 Cu_3 (PO_4)_4}$ with ${\rm A = Ca, Sr}$, ${\rm Bi_4 Cu_3 V_2 O_{14}}$ and azurite $\rm Cu_3(OH)_2(CO_3)_2$. We employ the physical consideration, the degenerate perturbation theory, the level spectroscopy analysis of the numerical diagonalization data obtained by the Lanczos method and also the density matrix renormalization group (DMRG) method. We investigate the mechanisms of the magnetization plateaux at $M=M_s/3$ and $M=(2/3)M_s$, and also show the precise phase diagrams on the $(J_2/J_1, J_3/J_1)$ plane concerning with these magnetization plateaux, where $M=\sum_{l=1}^{N} S_l^z$ and $M_s$ is the saturation magnetization. We also calculate the magnetization curves and the magnetization phase diagrams by means of the DMRG method.Comment: 21 pages, 29 figure

Onset of incommensurability in quantum spin chains

In quantum spin chains, it has been observed that the incommensurability occurs near valence-bond-solid (VBS)-type solvable points, and the correlation length becomes shortest at VBS-type points. Besides, the correlation function decays purely exponentially at VBS-type points, in contrast with the two-dimensional (2D) Ornstein-Zernicke type behavior in the other region with an excitation gap. We propose a mechanism to explain the onset of the incommensurability and the shortest correlation length at VBS-like points. This theory can be applicable for more general cases.Comment: 9 pages, 2 figure