Effects of magnetic anisotropy on spin dynamics of ferromagnetic frustrated chain
By exploiting density-matrix renormalization group techniques, we investigate the spin dynamics of a spin-1/2 one-dimensional - XXZ model with competing ferromagnetic and antiferromagnetic exchange couplings under applied magnetic fields. Numerical results of spin excitation spectra show that in the field-induced spin quadrupole regime, the longitudinal component has a gapless mode and the transverse component has a gapped mode irrespective of the exchange anisotropy. The excitation gap of the transverse spin excitation increases as the exchange anisotropy increases over the XY-like and Ising-like regions, demonstrating that two-magnon bound states are stabilized due to the easy-axis anisotropy.
Frustrated quantum spin systems in low dimensions have provided a fascinating playground to explore strange new phases that have no classical analogue, and rich physics of quantum phase transitions between them. One of the simplest models with frustration and low-dimensionality is a spin-1/2 one-dimensional - Heisenberg model, in which the nearest-neighbor exchange interaction and the next-nearest-neighbor exchange interaction compete with each other. Note that when is antiferromagnetic, frustration occurs in either case of antiferromagnetic or ferromagnetic . Recent studies on the ferromagnetic case in a magnetic field have revealed that the ground state is a vector chiral (VC) state in low magnetic fields, while it turns to an exotic Tomonaga-Luttinger-liquid state in high magnetic fields [1, 2, 3, 4]. In the high-field phase, the formation of two-magnon bound states leads to a spin quadrupole state. The longitudinal spin and quadrupole correlations are quasi-long-ranged, while the transverse spin correlation is short-ranged. The so-called spin nematic (SN) state is realized in the vicinity of the saturation, while the ground state shows a crossover to a spin-density-wave (SDW) state with lowering the magnetic field. The ground-state phase diagram in the presence of the exchange anisotropy has also been analyzed numerically . As for the experimental realization of the quadrupole state, LiCuVO has been studied as a possible candidate material. Neutron diffraction experiments have shown that the field dependence of an incommensurate magnetic peak appearing in high magnetic fields is well explained by theoretical predictions of the quadrupole state . Inelastic neutron scattering measurements have been carried out at zero magnetic field, indicating the existence of a two-spinon continuum , whereas the excitation spectrum in the high-field phase has not been reported yet.
To clarify the property of the field-induced quadrupole state from the viewpoint of the spin dynamics, we have studied the dynamical spin structure factor of the spin-1/2 one-dimensional - Heisenberg model by exploiting numerical methods such as a dynamical density-matrix renormalization group . In this paper, to gain further insight into the effect of the exchange anisotropy, we analyze the - XXZ model. We present numerical results on the anisotropic behavior of longitudinal and transverse spin excitation spectra, and discuss the dependence on the exchange anisotropy.
2 Model and method
We consider a spin-1/2 - XXZ model on a one-dimensional chain with sites, described by
where is the component of the spin-1/2 operator at site , and denote the ferromagnetic nearest-neighbor and the antiferromagnetic next-nearest-neighbor exchange interactions, respectively, is the anisotropy of the exchange interaction, and is the magnetic field along the direction. In the present study, we fix and , and examine the dependence on .
We investigate dynamical properties by using density-matrix renormalization group (DMRG) methods in open boundary conditions [9, 10]. To clarify the characteristics of anisotropic magnetic excitations, we calculate the and components of the dynamical spin structure factor, given by
where is the ground state and is its eigenenergy. We set an infinitesimal value to . The momentum representation of the spin operator in open boundary conditions is given by
For the calculation of the dynamical spin structure factor at a given magnetization with , we set the magnetic field to be the midpoint of the magnetization plateau of . We employ a dynamical DMRG method , in which we need three target states, i.e., , , and . We perform dynamical DMRG calculations with 32 sites at and for typical values of .
3 Numerical results
First, let us look at spin excitation spectra at zero magnetic field. In Fig. 1, we present intensity plots of and at . Figure 1(a) shows at for the isotropic Heisenberg case, where the system is in a Haldane dimer phase . Here, due to the SU(2) rotational invariance, as shown in Figs. 1(a) and 1(b). There is a sharp peak at , and a large amount of spectral weight is concentrated on this peak. We note that seems to be a finite value because of the finite-size effect. In fact, becomes small as the system size increases, and it is extrapolated to almost zero in the thermodynamic limit. This is consistent with the renormalization group study suggesting an exponentially small gap . The lower boundary of a continuum is found to form a sinusoidal dispersion, which corresponds to the des Cloizeaux-Pearson mode of spinon excitations in the limit of , i.e., decoupled antiferromagnetic chains . The sinusoidal dispersion is clearly visible due to the large intensity for , while that for is less distinct, since the spectral weight is distributed to a high energy region. Thus the distribution of the spectral weight is asymmetric with respect to [7, 14]. Note that the intensity extends to high energies beyond the upper boundary of the two-spinon continuum for , which is indicative of a four-spinon continuum suggested in LiCuVO .
In Figs. 1(c) and 1(d), we show and , respectively, at for the XY-like case, where the system is in a gapless chiral phase . Both and resemble those in the isotropic Heisenberg case in their overall spectral shapes, while they move to lower energy. We observe that the intensity at the lowest-energy peak is suppressed for , while it is enhanced for . Note that and have the same total intensity, and the difference of the distribution of the spectral weight in the space is essential. The spin excitation spectra are supposed to be gapless in the gapless chiral phase, but we do not clearly see the gapless nature in the present results because of the finite-size effect. In Figs. 1(e) and 1(f), we show and , respectively, at for the Ising-like case, where the system is in the so-called uudd phase [15, 16]. The overall spectral shapes are again similar to those of the isotropic Heisenberg case, while they are pushed up to higher energy in contrast to the XY-like case. In addition, the intensity at the lowest-energy peak is enhanced for , while it is suppressed for , which is the reverse of what we find in the XY-like case. In the uudd phase, the longitudinal spin excitation spectrum is gapless, while the transverse one is gapped, but these features are not clear because of the finite-size effect. The detailed finite-size scaling analysis is an important future problem, which will be reported elsewhere.
Now, we move on to the investigation in the magnetic field. In Fig. 2, we present intensity plots of and at , corresponding to the field-induced SDW regime. Figure 2(a) shows at for the isotropic Heisenberg case. We find a lowest-energy peak at . moves toward small momentum from the position at zero magnetic field. By analyzing for various values of , we have confirmed that changes with following a relation in the field-induced SDW/SN phase , which supports the bosonization result . On the other hand, is nearly zero, indicating a gapless mode of the longitudinal spin excitation. Note that we also find a ferromagnetic peak at the origin due to the finite magnetization. In Fig. 2(b), we show at . We find a lowest-energy peak at . In particular, we have verified that is extrapolated to a finite excitation energy in the thermodynamic limit, indicating a gapped mode of the transverse spin excitation. We notice that , i.e., it represents the spin excitation changing the magnetization by one. Since a finite energy is required to break a bound magnon pair, we have an excitation gap corresponding to the binding energy.
All these features found in the isotropic Heisenberg case are similarly observed for both the XY-like [Figs. 2(c) and 2(d)] and Ising-like [Figs. 2(e) and 2(f)] cases. That is, for , changes with as in the field-induced SDW/SN phase, and is nearly zero, indicating a gapless mode. For , is found at a finite energy, while it is extrapolated to a finite excitation energy in the thermodynamic limit, indicating a gapped mode. Here, we clearly see the difference of of between the XY-like and Ising-like cases. Compared to the isotropic Heisenberg case [Fig. 2(b)], decreases in the XY-like case [Fig. 2(d)], while it increases in the Ising-like case [Fig. 2(f)]. Thus the excitation gap is suppressed in the XY-like case, and it is enhanced in the Ising-like case.
To clarify the dependence of the excitation gap of the transverse spin excitation, we calculate a spin excitation energy, defined by
where is the lowest energy of the -site system in the subspace of at . Here, we use an ordinary DMRG method, not dynamical one, to obtain . In Fig. 3, we plot for several values of . The results in the field-induced SDW/SN regime are shown. Note that we find reentrant behavior in the XY-like case, i.e., the ground state changes from the VC state to the SDW/SN state and turns to the VC state again as the magnetization increases from zero to the saturation . We find that increases as increases over the XY-like and Ising-like regions. This indicates that the easy-plane anisotropy disfavors the formation of two-magnon bound states, and the easy-axis anisotropy stabilizes the two-magnon bound state, which is consistent with the behavior of correlation functions .
We have investigated the spin excitation dynamics of the spin-1/2 one-dimensional - XXZ model in the magnetic field by numerical methods. We have observed that in the field-induced SDW/SN phase, the longitudinal spin excitation is gapless, while the transverse spin excitation is gapped, irrespective of the exchange anisotropy. The excitation gap of the transverse spin excitation increases as the exchange anisotropy increases, which is indicative that the easy-axis anisotropy stabilizes the two-magnon bound state.
The author thanks T. Masuda, S. Maekawa, M. Mori, and T. Sugimoto for useful discussions. Numerical calculations have been supported by the supercomputer centers at the Japan Atomic Energy Agency and the Institute for Solid State Physics, the University of Tokyo.
-  Vekua T, Honecker A, Mikeska H J and Heidrich-Meisner F 2007 Phys. Rev. B 76 174420
-  Hikihara T, Kecke L, Momoi T and Furusaki A 2008 Phys. Rev. B 78 144404
-  Sudan J, Lüscher A and Läuchli M 2009 Phys. Rev. B 80 140402
-  Sato M, Momoi T and Furusaki A 2009 Phys. Rev. B 79 060406
-  Heidrich-Meisner F, McCulloch I P and Kolezhuk A K 2009 Phys. Rev. B 80 144417
-  Masuda T, Hagihala M, Kondoh Y, Kaneko K and Metoki N 2011 J. Phys. Soc. Jpn. 80 113705
-  Enderle M, Fåk B, Mikeska H J, Kremer R K, Prokofiev A and Assmus W 2010 Phys. Rev. Lett. 104 237207
-  Onishi H (Unpublished)
-  White S R 1992 Phys. Rev. Lett. 69 2863
-  Jeckelmann E 2002 Phys. Rev. B 66 045114
-  Furukawa S, Sato M, Onoda S and Furusaki A 2012 Phys. Rev. B 86 094417
-  Itoi C and Qin S 2001 Phys. Rev. B 63 224423
-  des Cloizeaux J and Pearson J J 1962 Phys. Rev. 128 2131
-  Ren J and Sirker J 2012 Phys. Rev. B 85 140410
-  Igarashi J 1989 J. Phys. Soc. Jpn. 58 4600
-  Tonegawa T, Harada I and Igarashi J 1990 Prog. Theor. Phys. Suppl. 101 513