We have considered the 1D spin-(1/2) Heisenberg model with added Dzyaloshinskii-Moriya interaction. The effect of a uniform magnetic field on the ground state phase diagram of the model is studied. We have mapped the model to an effective model which is known as the 1D XXZ in both uniform and staggered magnetic fields. By selecting a block of two or three spins, we have solved the Hamiltonian exactly. Our results show that the quantum phase transitions can be obtained from the block of pair or three spins.
1. Introduction
Study of the magnetic field effect on the ground state characteristics of chain model of antiferromagnetic (AF) spin-(1/2) has attracted much interest in recent years. The Hamiltonian of this model in a homogenous magnetic field h is given by
H=J∑j=1N(SjxSj+1x+SjySj+1y+SjzSj+1z)+h∑j=1NSjx,
where J>0 is the exchange coupling and h denotes the magnetic field. The exact solution for h=0 is obtained by Bethe ansatz [1]. The energy spectrum is gapless and the system is in the Luttinger liquid phase. In this phase, the spin correlation functions are in power form. When the system imposes by homogeneous magnetic field, its spectrum remains gapless up to the critical field hc=2J. Here the Pokrovsky-Talapov phase transition takes place and the ground state of the system becomes saturated ferromagnetic [2, 3]. In many experimental data, the resulting data are different by theoretical predications [4–10]. These differences are due to DM interaction [11, 12], its Hamiltonian can be written as follows:
H=J∑j=1N(SjxSj+1x+SjySj+1y+SjzSj+1z)+∑j=1ND⃗⋅(S⃗j×S⃗j+1)+h∑j=1NSjx,
where D⃗ is the DM vector and chooses in the z direction. [As DM interaction broken the fundamental SU(2) symmetry which is related to Heisenberg isotropic interaction]. It is known as the origin of many declination and creates many different qualitative effects. It specially creates an energy gap by the scale of gap≈(Dh)2/3 [13, 14]. DM interaction can act as a vector potential on the spin wave in the magnon spin Hall effect [15]. In ferromagnetic nanowires DM interaction has profound effect on the motion of domain walls [16]. It can also give rise to spin current and soliton in spin chains [17]. We emphasis that the studies of Heisenberg model commonly without DM interaction and the role of DM interaction on GS of Heisenberg model of AF by spin-(1/2) are considered less than the other models. As the exact solution of this model with DM interaction cannot be done, we need to do a great value of theoretical works.
2. Mapping to the XXZ Chain
In this section we try to map the Heisenberg chain in the presence of external magnetic field and DM interaction to the well-known XXZ chain, analytically. At first we have done a π/2 rotation around X, to convert z→y,y→z:H=J∑j=1N(SjxSj+1x+SjySj+1y+SjzSj+1z)+D∑j=1N(SjxSj+1y-SjySj+1x)+h∑j=1NSjx.
Using S± operators [15],
Sj±=Sjx±Sjy⟹Sjx=12(Sj++Sj-),Sjy=12i(Sj+-Sj-).
The interacting Hamiltonian terms can be calculated as follows:SjxSj+1y=14i(Sj++Sj-)(Sj+1+-Sj+1-)=14i(Sj+Sj+1+-Sj+Sj+1-+Sj-Sj+1+-Sj-Sj+1-),SjySj+1x=14i(Sj+-Sj-)(Sj+1++Sj+1-)=14i(Sj+Sj+1++Sj+Sj+1--Sj-Sj+1+-Sj-Sj+1-),SjxSj+1x=14(Sj++Sj-)(Sj+1++Sj+1-)=14(Sj+Sj+1++Sj+Sj+1-+Sj-Sj+1++Sj-Sj+1-),SjySj+1y=14i(Sj+-Sj-)(Sj+1+-Sj+1-)=14(-Sj+Sj+1++Sj+Sj+1-+Sj-Sj+1+-Sj-Sj+1-),SjxSj+1x+SjySj+1y=12(Sj+Sj+1-+Sj-Sj+1+).
Substituting these terms, we simply obtain the following:
HDM=D∑j=1N(SjxSj+1y-SjySj+1x)=D2i∑j=1N(Sj-Sj+1+-Sj+Sj+1-),HXXx=J2∑j=1N(Sj+Sj+1-+Sj-Sj+1++2SjzSj+1z).
Finally one can find the effective Hamiltonian:H1=Jcosφ∑j=1N((SjxSj+1x+SjySj+1y)+cosφSjzSj+1z)-h∑j=1NSjx-h′∑J=1N(-1)jSjy,
in which h′∝D⃗×h⃗ is known as staggered field and D=Jtanφ. This Hamiltonian is considered as XXZ chain Hamiltonian in homogenous and staggered field. It can be seen that AF can be created by DM interaction. It is shown theoretically that applying a staggered field on the AF chain of Heisenberg spin-(1/2), causes an energy gap in the spectrum of the system. In the absences of homogenous and staggered magnetic field (h=h′=0) the GS of the system is in the spin Luttinger liquid phase. In the presence of homogenous field, the exact solution is obtained by the Bethe ansatz methods. In the presence of staggered magnetic field, the exact solution is impossible. AF causes Neel ordered in y direction (unrotated z direction). When h′=0,h≠0, the magnetization in the h direction is an uniform and increasing function of h.
In what follows, we will consider two- and three-particle systems because they can be exactly calculated, and this phenomena takes place in some alternative spin chain that gives us a qualitative view of real interacting model which cannot be solved.
3. Two Particles Spin-(1/2) Systems
In this section we consider a special model which composed of two particle systems. We choose AF interacting for spins. H=J(S1⃗⋅S2⃗)+D(S1⃗×S2⃗)-h(S1x+S2x).
To solve this, we choose the singlet and triplet states for Hilbert space. We know that its eigenstates relate to S1zS2z as follows:
S1zS2z|↑↑〉=14|↑↑〉,S1zS2z|↑↓〉=-14|↑↓〉,S1zS2z|↓↑〉=-14|↓↑〉,S1zS2z|↓↓〉=14|↓↓〉.
Using these eigenstates, the singlet and triplet states are defined as follows:|S〉=12(|↑↓〉-|↓↑〉),|t1〉=|↑↑〉,|t0〉=12(|↑↓〉+|↓↑〉),|t-1〉=|↓↓〉.
In this [base kets], the matrix representation of Hamiltonian is as follows:H=(-34J-D220-D22-D22J4-22h00-22hJ4-22h-D220-22hJ4).
We diagonalize the Hamiltonian (11) and determine the GS of the systems as follows:
|GS〉=A|S〉+B(|t1〉+|t-1〉)+C|t0〉.
In which A, B, and C are as follows:
A=J4,B=D22,C=2h2.
In order to define the characteristics of the GS of the model in different subspaces of the GS phase diagram, we first calculated the magnetic order parameter Mx,y,z and Mstx,y,zas follows:Mx,y,z=〈GS|1N∑nSnx,y,z|GS〉,Mstx,y,z=〈GS|1N∑n(-1)nSnx,y,z|GS〉.
The averaging is calculated on GS. The magnetization along the applied field is plotted via h for different DM vectors in Figure 1. It can be seen that for small D=0.1, the magnetization for small fields remains close to zero. Increasing h up to h>hc, the magnetization becomes saturated. Due to quantum fluctuations in the presence of D, there is not any sharp transient in saturation point. For greater values of D, the magnetization increases by imposing magnetic field. In XXZ model, in h<hc magnetization has some fluctuations and for h>hc will be saturated. Imposing transverse external magnetic field, causes removing quantum fluctuation and the spins completely direct along the field, but the magnetization saturation cannot be seen because of broken symmetry in finite systems.
Magnetization along applied field,Mx, as a function of external field h, for exchange constant J=1 and different values of DM parameter D=0.1,0.5,1.0.
We show that in previous section that by imposing magnetic field h in the presence of DM interaction a staggered magnetization will be created normal to D⃗-h⃗ surface using the definition of staggered magnetization we can obtain the following:
Mstz=12〈GS|S1z-S2z|GS〉.
The staggered magnetization via h for D=0.1,0.5,1.0. is plotted in Figure 2. Figure 2 shows that applying homogeneous h, creates Neel order in z direction. That shows spin-flop ordering. There is a maximum for AM about 0.3. This is the same maximum of an independent dimer which weakly depends on D. Nersesyan et al. [18] shows that in one-dimensional spiral model by the interaction of second nearest neighbored which has a simple surface, a new phase is created due to breaking parity symmetry. This calls chiral phase and its order parameter is defined as follows:
χα=〈GS|(S1⃗×S2⃗)α|GS〉,
where α refers to x-, y-, z-axes. There are two types of gapped and gapless energy spectrum chiral phase [19–21].
Staggered magnetization perpendicular to the plan consists of applied magnetic field and DM vector as a function of external magnetic field for exchange energy J=1 and different values of DM interaction D=0.1,0.5,1.0.
For studying the GS phase diagram of the spin chain in the external magnetic field and DM interaction, we calculate the chiral order parameter as follows:
χy=〈GS|(S⃗1×S⃗2)y|GS〉=〈GS|S1zS2x-S1xS2z|GS〉.
In Figure 3, this chiral is plotted as a function of h for different values of D=0.1,0.5,1.0. As it can be seen, in the absence of DM interaction, the chiral order parameter does not exist in the ground state of the system. When external magnetic field applies, the order parameter changes and starts to increase by increasing the magnetic field. By continuing increase the magnetic field h, chiral order parameter will decrease and for sufficiently values of it will be disappeared. This is a state at which the system is in the saturated magnetization phase and the results are in agreement with analytical results of Section 2.
Chiral order parameter in y direction, χy as a function of h for J=1 and different values of DM interaction parameter D=0.1,0.5,1.0.
4. Three Particle Spin-(1/2) Systems
Consider a three-particle spin-(1/2) system whose Hamiltonian is as follows:H=J(S⃗1⋅S⃗2+S⃗2⋅S⃗3)+D(S⃗1×S⃗2+S⃗2×S⃗3)-h(S1x+S2x+S3x).
In order to solve the problem, we choose the [base kets] of the Hilbert space as eigenstates of the operator S1zS2zS3z as follows that areS1zS2zS3z|↑↑↑〉=18|↑↑↑〉,S1zS2zS3z|↑↑↓〉=-18|↑↑↓〉,S1zS2zS3z|↑↓↑〉=-18|↑↓↑〉,S1zS2zS3z|↓↑↑〉=-18|↓↑↑〉,S1zS2zS3z|↑↓↓〉=18|↑↓↓〉,S1zS2zS3z|↓↓↑〉=18|↓↓↑〉,S1zS2zS3z|↓↑↓〉=18|↓↑↓〉,S1zS2zS3z|↓↓↓〉=-18|↓↓↓〉.
In these states, the Hamiltonian matrix isH=(ABC0B000B0AF0B00C0-ACA0C00F000A0CB0A0-ABF00B0AB0AB00C0F00C000C0BCA).
In which
A=J4⟶B=D4-h2⟶C=-D2-h2⟶F=-h2.
We diagonalzed the Hamiltonian and found and the ground state of the system. Then order parameters as follows:Mstz=12〈GS|S1x+S2x+S3x|GS〉,Mstz=12〈GS|S1z-S2z+S3z|GS〉,χy=〈GS|(S1⃗×S2⃗)y+(S2⃗×S3⃗)y|GS〉,are calculated for this system, numerically that are shown in Figures 4 and 5. In the three spin case, we find the same behavior which we do not show. Specially three-particle spin-(1/2) system has exactly the same results of two particle system due to defining the order parameters in the particle number unit. In the XXZ model h<hc, the magnetization has some of the quantum fluctuation and the system will be in the magnetic field direction but the magnetization cannot be seen due to the symmetry breakdown in the finite systems.
Magnetization along applied field, Mx for three particle, as a function of external field h, for exchange constant J=1 and different values of DM parameter D=0.1,0.5,1.0.
Staggered magnetization perpendicular to the plan consists of applied magnetic field and DM vector as a function of external magnetic field for exchange energy J=1 and different values of DM interaction D=0.1,0.5,1.0.
5. Conclusion
In summary, we have studied Heisenberg model with added Dzyaloshinskii-Moriya interaction in the presence of uniform transverse magnetic field. To this end, we have first implemented a uniform rotation about x-axis in order to gauge away the added Dzyaloshinskii-Moriya interaction and produce an XXZ model with uniform magnetic and a staggered field normal to D⃗-h⃗ surface. We have chosen blocks of two and three spins to be able to solve the model exactly and have also shown how this scheme has capability to capture the right physics. We plotted magnetization, staggered magnetization, and chirality versus magnetic for different DM interactions strength. In the absence of DM interactions there is a sharp change from zero magnetization at h<hc to saturated region h>hc, [but increasing DM interaction favor to produce a nonzero magnetization at h<hc region] and there is no sharp transient any more. Staggered magnetization perpendicular to the plan consists of applied magnetic field and DM vector shows spin-flop ordering. The physical systems which prefer spiral chiral configuration usually have competing interactions. For example, when the nearest neighbor (NN) exchange interaction is weak, the next nearest neighbor (NNN) interaction or even Dzyaloshinskii-Moriya (DM) interaction becomes relatively significant. So, we measure chirality in our model versus magnetic field. Our calculations show that in the absence of DM interaction, the chiral order parameter is not in the ground state of the system. But by turning external magnetic field on, the order parameter will be changed and will be increased by increasing the magnetic field. Continuing the increase of magnetic decreases the chiral order and for sufficient values of it will disappear.
YangC. N.YangC. P.One-dimensional chain of anisotropic spin-spin interactions. I. Proof of Bethe's hypothesis for ground state in a finite system196615013213272-s2.0-3604905649810.1103/PhysRev.150.321PokrovskyV. L.TalapovA. L.theory of two-dimensional incommensurate crystals198051134148GriffithsR. B.Magnetization curve at zero temperature for the antiferromagnetic heisenberg linear chain19641333A768A7752-s2.0-3614901543010.1103/PhysRev.133.A768DenderD. C.HammarP. R.ReichD. H.BroholmC.AeppliG.Direct observation of field-induced incommensurate fluctuations in a one-dimensional S = 1/2 antiferromagnet1997799175017532-s2.0-0001514970SirkerJ.WeisseA.SushkovO. P.The field-induced magnetic ordering transition in TlCuCl32005741291342-s2.0-29544449791SakaiT.ShibaH.Numerical study of a model for NENP: one-dimensional S = 1 antiferromagnet in a staggered field19946338678712-s2.0-21344495712ChaboussantG.JulienM. H.Fagot-RevuratY.HansonM.LévyL. P.BerthierC.HorvatićM.PiovesanaO.Zero temperature phase transitions in spin-ladders: phase diagram and dynamical studies of Cu2(C5H12N2)2Cl41998621671812-s2.0-0001679422KageyamaH.YoshimuraK.SternR.MushnikovN. V.OnizukaK.KatoM.KosugeK.SlichterC. P.GotoT.UedaY.Exact dimer ground state and quantized magnetization plateaus in the two-dimensional spin system SrCu2(BO3)219998215316831712-s2.0-4243717268CépasO.KakuraiK.RegnaultL. P.ZimanT.BoucherJ. P.AsoN.NishiM.KageyamaH.UedaY.Dzyaloshinski-moriya interaction in the 2D spin gap system SrCu2(BO3)220018716142-s2.0-42749098364167205JaimeM.mjaime@lanl.govCorreaV. F.HarrisonN.BatistaC. D.KawashimaN.KazumaY.JorgeG. A.SternR.HeinmaaI.ZvyaginS. A.SasagoY.UchinokuraK.Magnetic-field-induced condensation of triplons in Han Purple pigment BaCuSi2O6200493810.1103/PhysRevLett.93.087203087203DzyaloshinskyI.A thermodynamic theory of "weak" ferromagnetism of antiferromagnetics195844241255MoriyaT.New mechanism of anisotropic superexchange interaction19604522823010.1103/PhysRevLett.4.228OshikawaM.AffleckI.Field-induced gap in S = 1/2 antiferromagnetic chains19977915288328862-s2.0-0001233904SatoM.OshikawaM.Coupled S = 1/2 Heisenberg antiferromagnetic chains in an effective staggered field20046951254406OnoseY.IdeueT.KatsuraH.ShiomiY.NagaosaN.TokuraY.Observation of the magnon hall effect201032959892972992-s2.0-7795479262010.1126/science.1188260TretiakovO. A.AbanovA.Current driven mnagnetization dynamics in ferromagnetic nanowires with a dzyaloshinskii-moriya interaction2010105152-s2.0-7795758410410.1103/PhysRevLett.105.157201157201BostremI. G.KishineJ. I.OvchinnikovA. S.Theory of spin current in chiral helimagnets20087862-s2.0-5084913991910.1103/PhysRevB.78.064425064425NersesyanA. A.GogolinA. O.ElerF. H. L.Incommensurate spin correlations in spin-1/2 frustrated two-Leg heisenberg ladders19988149109132-s2.0-0000248349KaburagiM.KawamuraH.HikiharaT.Spin and chiral orderings of frustrated quantum spin chains19996810318531882-s2.0-0033260420HikiharaT.KaburagiM.KawamuraH.Ground-state phase diagrams of frustrated spin-S XXZ chains: chiral ordered phases20016317174430117443011TakahashiM.1999Cambridge University Press