Phase Diagram of spin 1/2 XXZ Model With DzyaloshinskiiMoriya Interaction
Abstract
We have studied the phase diagram of the one dimensional XXZ model with DzyaloshinskiiMoriya (DM) interaction. We have applied the quantum renormalization group (QRG) approach to get the stable fixed points, critical point and the scaling of coupling constants. This model has three phases, ferromagnetic, spinfluid and Néel phases which are separated by a critical line which depends on the DM coupling constant. We have shown that the staggered magnetization is the order parameter of the system and investigated the influence of DM interaction on the chiral ordering as a helical magnetic order.
pacs:
75.10.Pq,73.43.Nq,03.67.Mn,64.60.aeI Introduction
Quantum phase transition has been one of the most interesting topic in the area of strongly correlated systems in the last decade. It is a phase transition at zero temperature where the quantum fluctuations play the dominant role vojta (). Suppression of the thermal fluctuations at zero temperature introduces the ground state as the representative of the system. The properties of the ground state may be changed drastically shown as a nonanalytic behavior of a physical quantity by reaching the quantum critical point (QCP). This can be done by tuning a parameter in the Hamiltonian, for instance the magnetic field or the amount of disorder. The ground state of a typical quantum many body systems consist of a superposition of a huge number of product states. Understanding this structure is equivalent to establishing how subsystems are interrelated, which in turn is what determines many of the relevant properties of the system. In Mott insulators the Heisenberg interaction is in most cases the dominant source of coupling between local moments, and most theoretical investigations are based on modeling in which only this type of interaction is included. Recently some novel magnetic properties in antiferromagnetic (AF) systems were discovered in the variety of quasione dimensional materials that are known to belong to an antisymmetric interaction of the form which is known as the DzyaloshinskiiMoriya (DM)interaction. The relevance of antisymmetric superexchange interactions in spin Hamiltonian which leads to either a week ferromagnetic (F) or helical magnetic distortion in quantum AF systems, has been introduced phenomenologically by DzyaloshinskiiDzyaloshinskii (). A microscopic model of antisymmetric exchange interaction was first proposed by MoriyaMoriya () which showed that such interactions arise naturally in perturbation theory due to the spinorbit coupling in magnetic systems with low symmetry and is essentially an extension of the Anderson superexchange mechanismAnderson () that shows for spinflip hopping of electrons. Since it (DM interaction) breaks the fundamental SU(2) symmetry of the Heisenberg interactions, it is at the origin of many deviations from pure heisenberg behavior such as cantingCoffey () or small gapsDender1 (); Oshikawa1 (); Zhao (); Fouet (); Chernyshev (). A number of AF systems expected to be described by DM interaction, such as Dender1 (); Dender2 (), Kohgi (); Fulde (); Oshikawa2 (), Tsukada (), , Grande () and Greven (), which exhibit unusual and interesting magnetic properties in the presence of quantum fluctuations and/or applied magnetic fieldGrande (); Yildirim (); Katsumata (). Also belonging to the class of DM antiferromagnets is , which is a parent compound of hightemperature superconductorsKastner (). This has stimulated extensive investigation on the physical properties of the DM interaction. However, This interaction is rather difficult to Handel analytically, which has brought much uncertainty in the interpretation of experimental data and has limited our understanding of many interesting quantum phenomena of lowdimensional magnetic materials.
In the present paper, we have considered the one dimensional XXZ model with DM interaction by implementing the quantum renormalization group (QRG) method. In the next section the QRG approach will be explained and the renormalization of coupling constant are obtained. In section (III), we will obtain the phase diagram, fixed points, critical points and calculate the staggered magnetization as the order parameter of the underling quantum phase transition. We will also introduce the chiral order as an ordering which is produced by DM interaction. The exponent which shows the divergence of correlation function close to the critical point , the dynamical exponent and the exponent which shows the vanishing of staggered magnetization near the critical point will also be calculated. In Sec (IV), discussion concludes the paper.
Ii Quantum renormalization group
The main idea of the RG method is the elimination or the thinning of the degrees of freedom caried out step by step in an iteration procedure. Here, we used the well known Kadanoff block method as it is both well suited to perform analytical calculations in the lattice models and is conceptually straightforward to be extended to the higher dimensionsmiguel1 (); miguelbook (); Langari (); Jafari (). In the Kadanoff’s method, the lattice is divided into blocks in which the Hamiltonian can be exactly diagonalized. Selecting a number of lowlying eigenstates of the blocks the full Hamiltonian is projected onto these eigenstates and an the effective (renormalized) Hamiltonian is obtained.
The Hamiltonian of XXZ model with DM interaction in the direction on a periodic chain of sites is
(1)  
where the is exchange constant, is the strength of component of DM interaction and the easyaxis anisotropy defined by which can be positive and negative. The positive and negative corresponds to the antiferromagnetic and ferromagnetic (F) cases, respectively. refers to the component of the Pauli matrix at site . By implement rotation around axis on odd or even sites, the AF case of Hamiltonian () is mapped on the F case () with opposite sign of anisotropy,
(2)  
So we can restrict ourselves to AF case () with and arbitrary anisotropy ( and ) without loss of generality.
The effective Hamiltonian to the first order RG approximation is
We consider a threesite block procedure defined in Fig.(1). The block Hamiltonian , its eigenstates and eigenvalues are given in Appendix A. The threesite block Hamiltonian has four doubly degenerate eigenvalues (see Appendix A). is the projection operator of the ground state subspace defined by , Where and are the doubly degenerate ground states, and are the renamed base kets in the effective Hilbert space. For each block we keep two states ( and ) to define the effective (new) site. Thus, the effective site can be considered as having a spin 1/2. Due to the level crossing which occurs for the eigenstates of the block Hamiltonian, the projection operator () can be different depending on the coupling constants. Therefore, we must specify different regions with the corresponding ground states. As Fig.(5) shows there are two regions with different eigenstates which are separated by where a level crossing occurs. In region (A) the ground state is the doublydegenerate ferromagnetic state and while in region (B) and are the degenerate ground states. At the level crossing () the ground state is 4fold degenerate (, , , ). A summary of this information is given in Fig.(5) of appendix A.
In the following, we will classify the RG equation of the regions where each of this states represent the ground state.
ii.1 Region (A): is the ground state.
In this region the effective Hamiltonian in the first order correction is similar to the initial one, i.e,
(3) 
where , and are the renormalized coupling constants. The new renormalized coupling constants are found to be functions of the original ones given by the following equations,
(4)  
ii.2 Region (B): is the ground state.
In this region the effective Hamiltonian to the first order corrections leads to the ferromagnetic Ising model
where
Iii Phase Diagram
iii.1 Region (A)
For simplicity we have separated this region into positive anisotropy and negative anisotropy sectors.

In the positive anisotropy sector the RG equations show that the coupling, representing the energy scale, approaching zero by iterating RG procedure. Thus, at the zero temperature,the quantum phase transition is the result of competition between the anisotropy () and the DM coupling constant (). In the region of planar anisotropy , the symmetric interactions () is known not to support any kind of long range order and the ground state is the so called spinfluid (SF) state. Increasing the amount of anisotropy is necessary to stabilize the spin alignment. For the ground state is the Néel ordered state. In the case of , the anisotropy constant () and antisymmetric (DM) coupling are in competition with each other. The latter thus destroys the ordering tendency of the former and defers creating of Néel order. Our RG equations show that the phase boundary between the SF and Néel phases which depends on the DM coupling is (see Fig.(2)) which agrees with the phase boundary reported in Ref.Alcaraz, . This critical line coincide with boundary line which obtained by classical approximation (see appendix C). The RG equations (Eq.(4)) express the DM coupling dose not flow under RG transformations, and the anisotropy coupling goes to zero () in SF phase while it scales to infinity () in the Néel phase.
We have linearized the RG flow at the critical line and found one relevant and one marginal directions. The eigenvalues of the matrix of linearized flow are , . The corresponding eigenvectors in the coordinates are , . The marginal direction corresponds to the tangent line of the critical line and the relevant direction shows the direction of anisotropy’s flow (Fig.(2)).
However we have found the boundary of the SFNéel transition by calculating the staggered magnetization (see appendix B) in the direction as an order parameter (Fig.(3) and Fig.(4)),
(5) is zero in the SF phase and has a nonzero value in the Néel phase. Thus the staggered magnetization is the proper order parameter to represent the SFNéel transition. We have plotted versus and versus in Fig.(3) and Fig.(4), respectively. In Fig.(3) it is obvious that the staggered magnetization goes to zero continuously at the critical value of which shows the destruction of Néel order. The critical value where the staggered magnetization vanishes above it, increases by increasing of anisotropy, this means SFNéel transition point () depends on the DM interaction. It is seen in Fig.(4) that the staggered magnetization is zero for (SF phase) and has a nonzero value for (Néel phase), while enhancing of DM coupling defer the creation of Néel order.
Moreover, to study the influence of DM coupling we have calculated the chiral order Kawamura (); Kaburagi ()() in the direction as the helical magnetization in one dimensionJafari2 () which is created by DM interaction.
Unfortunately the chiral order is not a self similar operator under RG transformations. In fact an term of Hamiltonian () shows up in the renormalized chiral order under RG. Thus, we have to calculate the sum of chiral and term which we have represented it by in the following equation,
(6) In Fig.(3) and Fig.(4) the DC order has plotted versus and respectively, for different values of and . The figures manifest that the order enhances by increasing of DM coupling and reduces with increasing of anisotropy parameter.
We have also calculated the critical exponents at the critical line . In this respect, we have obtained the dynamical exponent, the exponent of order parameter and the diverging exponent of the correlation length. This corresponds to reach the critical point from the Néel phase by approaching . The dynamical exponent is , the staggered magnetization goes to zero like with . The correlation length diverges with exponent . The remarkable result of these exponents is the independence of their values on the D value and equality of them with the correspondening ones in the model. The detail of this calculation is similar to what is presented in Ref.Langari ().

In this sector, the effective Hamiltonian is similar to the positive anisotropy case with the same coupling constants. For the ground state is the spinfluid phase and decreasing the anisotropy causes the ground state of the three site Hamiltonian changes by level crossing at where the RG equations should be reconstructed. However, the remarkable result is that the level crossing line which got by three site Hamiltonian coincides the critical line of this model in the thermodynamic systemAlcaraz (). The RG equations express the DM coupling dose not flow under RG transformations, and the anisotropy coupling goes to zero (). Thus the line is the position of stable fixed points where the RG flow is freezed.
iii.2 Region (B)
As we pointed out in sec.IIB the original Hamiltonian is mapped to the ferromagnetic Ising model. Ising model remains unchanged under RG as the stable fixed point and its properties are well known. We call this region as the ferromagnetic Ising phase.
Iv Summery and conclusions
We have applied the RG transformation to obtain the phase diagram, staggered magnetization and helical magnetization of model with DM interaction. In the positive anisotropy region, tuning the anisotropy coupling makes the system to fall into different phases, i.e Néel phase with nonzero order parameter and spinfluid one with vanishing order parameter as characterized by the staggered magnetization. The RG equations state that the system has fixed points at , and . The fixed points and are attractive and correspond to the SpinFluid and Néel phases, respectively. The fixed points at are repulsive and correspond to the critical points of this model, in the other word it is the critical line of this Hamiltonian. However, in the negative anisotropy region, the level crossing line which is obtained by three site block Hamiltonian eigenvalues, is the critical line of infinite size system and separates the ferromagnetic and spinfluid phases. Unfortunately we can not calculate the chiral order explicitly by RG method. To survey the influence of DM interaction and helical magnetization, the numerical Lanczos computation is in progress.
Acknowledgements.
The authors would like to thank J. Abouie, M. Kargarian and M. Siahatgar for fruitful discussions. This work was supported in part by the Center of Excellence in Complex Systems and Condensed Matter (www.cscm.ir).V Appendix
v.1 The block Hamiltonian of three sites, its eigenvectors and eigenvalues
We have considered the threesite block (Fig.(1)) with the following Hamiltonian
The interblock ( and intrablock () Hamiltonian for the three sites decomposition are
where refers to the component of the Pauli matrix at site of the block labeled by . The exact treatment of this Hamiltonian leads to four distinct eigenvalues which are doubly degenerate. The ground, first, second and third excited state energies have the following expressions in terms of the coupling constants.
(7)  
(8)  
(9)  
where is .
and are the eigenstates of . In Fig.(5) we have presented the different regions where the specified state is the ground state of the block Hamiltonian. In the region (A) the projection operator is
The Pauli matrices in the effective Hilbert space have the following transformations
In the region (B) () the projection operator is
and the Pauli matrices in the effective Hilbert space have the following transformations
v.2 Order Parameter and Chiral Order
v.2.1 Staggered magnetization
Generally, any correlation function can be calculated in the QRG scheme. In this approach, the correlation function at each iteration of RG is connected to its value after an RG iteration. This will be continued to reach a controllable fixed point where we can obtain the value of the correlation function. The staggered magnetization in direction can be written
(10) 
where is the Pauli matrix in the th site and is the ground state of chain. The ground state of the renormalized chain is related to the ground state of the original one by the transformation, .
This leads to the staggered configuration in the renormalized chain. The staggered magnetization in direction is obtained
(11)  
where is the staggered magnetization at the th step of QRG and is defined by .
v.2.2 Chiral Order
The chiral order which is the proper function to detect the helical magnetization in the systems can be written
(12) 
As we mentioned in the section III, the term of the Hamiltonian shows up to the chiral order under RG. The term order is written
(13) 
In this case the calculating of the chiral order being elaborate, because of the unknown effect of the term on the ground state of system at fixed point (). To simplification of the calculation, we transform the with DM interaction Hamiltonian (Eq.(1)) to the Ising model with DM interactionJafari2 () (IDM) by implement a nonlocal transformation which shows a rotation about the axis at site where . We have calculated the chiral order (Eq.(12)) of IDM in Ref.[Jafari2 ()]. By implement the inverse of transformation, the chiral order in IDM model transforms to the DC order where introduced in Eq.(6). The DC order has been plotted in Fig.(3) and Fig.(4) versus and .
v.3 Classical Approximation
In the classical approximation the spins are considered as classical vectors which form the spiral structure with a pitch angle between neighboring spins and canted angle
The classical energy per site for the with DM interaction Hamiltonian (Eq.(1)) is
The minimization of classical energy with respect to the angles and shows that there are two different regions. (I) , the minimum of energy is obtained by arbitrary and which show the spins projection on axis is nonzero and spins have the helical structure (see Fig.6) in the plain. In this region the minimum classical energy is
(14) 
(II) , the energy is minimized by and arbitrary or arbitrary and , which correspond respectively to the configurations with nonzero value of spins projection on axis with helical structure of spins projection in the plain and disorder configuration. In this region the minimum classical energy is
(15) 
Vi Canonical Transformation
The Hamiltonian of model with DM interaction (Eq.(1)) has the global symmetry. This Hamiltonian is mapped to the well known chain via a canonical transformationAristov (); Alcaraz (),
(16) 
which gives
(17) 
The symmetry of initial Hamiltonian survive in the transformed Hamiltonian too, but at the symmetry breaks to the local symmetry.
References
References
 (1) M. Vojta, Rep. Prof. Phys. 66, 2069 (2003) and references therein.
 (2) I. Dzyaloshinskii, J. Phys. Chem. Solids 4, 241 (1958).
 (3) T. Moriya, Phys. Rev 120, 91 (1960).
 (4) P. W. Anderson, Phys. Rev 115, 2 (1959).
 (5) D. Coffey, K. S. Trugman, Phys. Rev. B. 42, 6509 (1990).
 (6) D. C. Dender, P. R. Hammar,D. H. Reich, C. Broholm , and G. Aeppli Phys. Rev. Lett 79, 1750 (1997).
 (7) M. Oshikawa, I. Affleck, Phys. Rev. Lett 79, 2883 (1997).
 (8) J. Z. Zhao, X. Q. Wang, T. Xiang, Z. B. Su, L. Yu, Phys. Rev. Lett 90, 20204 (2003).
 (9) J.B. Fouet, O. Tchernyshyov, F. Mila, Phys. Rev. B. 70, 174427 (2004); J.B. Fouet, F. Mila, D. Clarke, H. Youk, O. Tchernyshyov. P.Fendley, R. M. Noack, Phys. Rev. B. 73, 214405 (2006).
 (10) A. L. Chernyshev, Phys. Rev. B. 72, 174414 (2005).
 (11) D. C. Dender, D. Dvidovic, D. H. Reich, C. Broholm , and G. Aeppli Phys. Rev. B 53, 2583 (1996).
 (12) M. Kohgi, K. Iwasa, J. Mignot, B. Fak, P. Gegenwart, M. Lang, A. Ochiai, H. Aoki, and T. Suzuki, Phys. Rev. Lett 86, 2439 (2000).
 (13) P. Fulde B. Schmidt, and P. Thalmeier, Europhys. Lett 31, 323 (1995).
 (14) M. Oshikawa, K. Ueda, H. Aoki, A. Ochiai and M. Kohgi, J. Phys. Soc. Jpn 68, 3181 (1999); H. Shiba, K. Udea, and O. Sakai, J. Phys. Soc. Jpn 69, 1493 (2000)
 (15) I. Tsukada, J. T. Takeya, T. Masuda and K. Uchinokura, Phys. Rev. Lett 87, 127203 (2001)
 (16) b. Grande and Hk. MllerBuschbaum, Z. Anorg. Allg. Chem 417, 68 (1975).
 (17) M. Greven, R. J. Birgeneau, Y. Endoh, M. A. Kastner, M. Matsuda, and G. Shirane, Z. Phys. B 96, 465 (1995).
 (18) A. B. Harris, A. Aharony and O. E. Wohlman, Phys. Rev. B. 52, 10239 (1995).
 (19) K. Katsumata, M. Hagiwara,Z. Honda, J. Satooka, A. Aharoy, R. J. Birgeneau, F. C. Chou, O. E. Wohlman, A. B. Harris,M. A. Kastner, Y. J. Kim, and Y. S. Lee, Europhys. Rev. Lett 54, 508(2001).
 (20) M. A. Kastner, R. J. Birgeneau, G. Shirane and Y. Endoh, Rev. Mod. Phys 70, 897 (1998).
 (21) M. A. MartinDelgado and G. Sierra, Int. J. Mod, Phys. A 11, 3145 (1996).
 (22) G. Sierra and M. A. Martin Delgado, in Strongly Correlated Magnetic and Superconducting Systems, Lecture Notes in Physics Vo1. 478 (springer, Berlin, 1997).
 (23) A. Langari, Phys. Rev. B 58, 14467 (1998); it ibid 69, 100402(R) (2004).
 (24) R. Jafari, A. Langari, Phys. Rev. B 76, 014412 (2007); R. Jafari, A. Langari, Physica A 364, 213 (2006).
 (25) F. C. Alcaraz and W. F. Wreszinski, J. Stat. Phys. 58, 45 (1990).
 (26) H. Kawamura, Phys. Rev. B 38, 4916 (1988).
 (27) M. Kaburagi, and H. Kawamura, T. Hikihara, J. Phys. Soc. Jpn. 68, 3185 (1999);T. Hikihara, M. Kaburagi, and H. Kawamura, arXiv:condmat/0007095v2; T. Hikihara, M. Kaburagi, and H. Kawamura, arXiv:condmat/0010283v1.
 (28) R. Jafari, M. Kargarian, A. Langari, and M. Siahatgar, Phys. Rev. B , (2008).
 (29) D. N. Aristov, S. V. Maleyev, Phys. Rev. B 62, R751