ISRN Condensed Matter PhysicsVolume 2011 (2011), Article ID 980192, 6 pagesdoi:10.5402/2011/980192

Research Article

Heisenberg Model with Added Dzyaloshinskii-Moria Interaction

M. R. Soltani,¹ J. Vahedi,¹ M. R. Abolhassani,¹ and A. A. Masoudi²

¹Department of Physics, Science and Research Branch, Islamic Azad University, Tehran, Iran²Department of Physics, Alzahra University, P.O. Box 19585-466, Tehran, Iran

Received 25 October 2011; Accepted 30 November 2011

Academic Editors: C. Andreani and W. Selke

Copyright © 2011 M. R. Soltani et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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 is given by where is the exchange coupling and denotes the magnetic field. The exact solution for 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 . 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: where is the DM vector and chooses in the 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 [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 rotation around X, to convert :

Using operators [15],

The interacting Hamiltonian terms can be calculated as follows: Substituting these terms, we simply obtain the following: Finally one can find the effective Hamiltonian: in which is known as staggered field and . 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 () 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 direction (unrotated z direction). When , the magnetization in the direction is an uniform and increasing function of .

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. To solve this, we choose the singlet and triplet states for Hilbert space. We know that its eigenstates relate to as follows: Using these eigenstates, the singlet and triplet states are defined as follows: In this [base kets], the matrix representation of Hamiltonian is as follows: We diagonalize the Hamiltonian (11) and determine the GS of the systems as follows: In which , , and are as follows: 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 and as follows: The averaging is calculated on GS. The magnetization along the applied field is plotted via for different DM vectors in Figure 1. It can be seen that for small , the magnetization for small fields remains close to zero. Increasing up to , the magnetization becomes saturated. Due to quantum fluctuations in the presence of , there is not any sharp transient in saturation point. For greater values of , the magnetization increases by imposing magnetic field. In XXZ model, in magnetization has some fluctuations and for 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.

Figure 1: Magnetization along applied field,, as a function of external field , for exchange constant and different values of DM parameter .

We show that in previous section that by imposing magnetic field in the presence of DM interaction a staggered magnetization will be created normal to surface using the definition of staggered magnetization we can obtain the following: The staggered magnetization via for . is plotted in Figure 2. Figure 2 shows that applying homogeneous , creates Neel order in 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 . 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: where refers to -, -, -axes. There are two types of gapped and gapless energy spectrum chiral phase [19–21].

Figure 2: Staggered magnetization perpendicular to the plan consists of applied magnetic field and DM vector as a function of external magnetic field for exchange energy and different values of DM interaction .

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: In Figure 3, this chiral is plotted as a function of for different values of . 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 , 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.

Figure 3: Chiral order parameter in direction, as a function of for and different values of DM interaction parameter .

4. Three Particle Spin-(1/2) Systems

Consider a three-particle spin-(1/2) system whose Hamiltonian is as follows: In order to solve the problem, we choose the [base kets] of the Hilbert space as eigenstates of the operator as follows that are In these states, the Hamiltonian matrix is In which We diagonalzed the Hamiltonian and found and the ground state of the system. Then order parameters as follows: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 , 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.

Figure 4: Magnetization along applied field, for three particle, as a function of external field , for exchange constant and different values of DM parameter .

Figure 5: Staggered magnetization perpendicular to the plan consists of applied magnetic field and DM vector as a function of external magnetic field for exchange energy and different values of DM interaction .

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 -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 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 to saturated region , [but increasing DM interaction favor to produce a nonzero magnetization at 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.

References

C. N. Yang and C. P. Yang, “One-dimensional chain of anisotropic spin-spin interactions. I. Proof of Bethe's hypothesis for ground state in a finite system,” Physical Review B, vol. 150, no. 1, pp. 321–327, 1966.

V. L. Pokrovsky and A. L. Talapov, “theory of two-dimensional incommensurate crystals,” Soviet Physics, vol. 51, pp. 134–148, 1980.

R. B. Griffiths, “Magnetization curve at zero temperature for the antiferromagnetic heisenberg linear chain,” Physical Review, vol. 133, no. 3, pp. A768–A775, 1964.

D. C. Dender, P. R. Hammar, D. H. Reich, C. Broholm, and G. Aeppli, “Direct observation of field-induced incommensurate fluctuations in a one-dimensional S = 1/2 antiferromagnet,” Physical Review Letters, vol. 79, no. 9, pp. 1750–1753, 1997.

J. Sirker, A. Weisse, and O. P. Sushkov, “The field-induced magnetic ordering transition in TlCuCl₃,” Journal of the Physical Society of Japan, vol. 74, pp. 129–134, 2005.

T. Sakai and H. Shiba, “Numerical study of a model for NENP: one-dimensional S = 1 antiferromagnet in a staggered field,” Journal of the Physical Society of Japan, vol. 63, no. 3, pp. 867–871, 1994.

G. Chaboussant, M. H. Julien, Y. Fagot-Revurat et al., “Zero temperature phase transitions in spin-ladders: phase diagram and dynamical studies of Cu₂(C₅H₁₂N₂)₂Cl₄,” European Physical Journal B, vol. 6, no. 2, pp. 167–181, 1998.

H. Kageyama, K. Yoshimura, R. Stern et al., “Exact dimer ground state and quantized magnetization plateaus in the two-dimensional spin system SrCu₂(BO₃)₂,” Physical Review Letters, vol. 82, no. 15, pp. 3168–3171, 1999.

O. Cépas, K. Kakurai, L. P. Regnault et al., “Dzyaloshinski-moriya interaction in the 2D spin gap system SrCu₂(BO₃)₂,” Physical Review Letters, vol. 87, no. 16, Article ID 167205, pp. 1–4, 2001.

M. Jaime, V. F. Correa, N. Harrison et al., “Magnetic-field-induced condensation of triplons in Han Purple pigment BaCuSi₂O₆,” Physical Review Letters, vol. 93, no. 8, Article ID 087203, 2004.

I. Dzyaloshinsky, “A thermodynamic theory of "weak" ferromagnetism of antiferromagnetics,” Journal of Physics and Chemistry of Solids, vol. 4, no. 4, pp. 241–255, 1958.

T. Moriya, “New mechanism of anisotropic superexchange interaction,” Physical Review Letters, vol. 4, no. 5, pp. 228–230, 1960.

M. Oshikawa and I. Affleck, “Field-induced gap in S = 1/2 antiferromagnetic chains,” Physical Review Letters, vol. 79, no. 15, pp. 2883–2886, 1997.

M. Sato and M. Oshikawa, “Coupled S = 1/2 Heisenberg antiferromagnetic chains in an effective staggered field,” Physical Review B, vol. 69, no. 5, Article ID 54406, 12 pages, 2004.

Y. Onose, T. Ideue, H. Katsura, Y. Shiomi, N. Nagaosa, and Y. Tokura, “Observation of the magnon hall effect,” Science, vol. 329, no. 5989, pp. 297–299, 2010.

O. A. Tretiakov and A. Abanov, “Current driven mnagnetization dynamics in ferromagnetic nanowires with a dzyaloshinskii-moriya interaction,” Physical Review Letters, vol. 105, no. 15, Article ID 157201, 2010.

I. G. Bostrem, J. I. Kishine, and A. S. Ovchinnikov, “Theory of spin current in chiral helimagnets,” Physical Review B, vol. 78, no. 6, Article ID 064425, 2008.

A. A. Nersesyan, A. O. Gogolin, and F. H. L. Eler, “Incommensurate spin correlations in spin-1/2 frustrated two-Leg heisenberg ladders,” Physical Review Letters, vol. 81, no. 4, pp. 910–913, 1998.

M. Kaburagi, H. Kawamura, and T. Hikihara, “Spin and chiral orderings of frustrated quantum spin chains,” Journal of the Physical Society of Japan, vol. 68, no. 10, pp. 3185–3188, 1999.

T. Hikihara, M. Kaburagi, and H. Kawamura, “Ground-state phase diagrams of frustrated spin-S XXZ chains: chiral ordered phases,” Physical Review B, vol. 63, no. 17, pp. 1744301–17443011, 2001.

M. Takahashi, Thermodynamics of One Dimensional Solvable Models, Cambridge University Press, 1999.