Научная статья на тему 'Simulation of the magnetic structure upon ordering eg orbitals to the quasi-oneand quasi-two-dimensional magnets by quantum Monte Carlo metod'

Simulation of the magnetic structure upon ordering eg orbitals to the quasi-oneand quasi-two-dimensional magnets by quantum Monte Carlo metod Текст научной статьи по специальности «Физика»

CC BY
172
43
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
АЛЬТЕРНИРОВАНИЕ ОБМЕНА / ОРБИТАЛЬНОЕ И МАГНИТНОЕ УПОРЯДОЧЕНИЕ / СТРАЙПЫ / ALTERNATION OF EXCHANGE / ORBITAL AND MAGNETIC ORDERING / STRIPES

Аннотация научной статьи по физике, автор научной работы — Aplesnin S. S., Moskvin A. I., Piskunova N. I.

The exchange mechanism influence on the electrons ordering on eg orbitals in a chain and in a two-dimensional Heisenberg model with exchange anisotropy for a S = 1/2 spin is determined by a quantum Monte Carlo method. The existence regions of a long-range quasi-oneand two-dimensional antiferromagnetic order with the special exchange topology are determined. The plateau existence region in the field dependence of magnetization, as well as the wave vector of the magnetic structure modulation with Q = p/2 in the (magnetic field exchange alternating) plane, is determined.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Simulation of the magnetic structure upon ordering eg orbitals to the quasi-oneand quasi-two-dimensional magnets by quantum Monte Carlo metod»

Библиографический список

1. Вологдин, Э. Как возникают и звучат ошибки квантования / Э. Вологдин // Звукорежиссер. 2006. С. 28-41 ; 2007. С. 32-40.

2. Оболонин, И. А. Оценка влияния аналого-цифрового преобразования на искажения при компрессии цифровых аудиоданных / И. А. Оболонин // Вестн. Сиб. гос. ун-та телекоммуникаций и информатики : сб. науч. тр. / Сиб. гос. ун-т телекоммуникаций и информатики. Новосибирск, 2008. №> 2. С. 67-71.

I. A. Obolonin, N. A. Rygovskaya

FILTRATION INFLUENCE ANALYSIS AT PRIMARY PULSED-CODE CONVERSION ON INPUT SIGNALS DISTORTION OF CODERS WITH AUDIO DATA COMPRESSION

The influence of selective circuits of digital-analogue converters at primary pulsed-code conversion on input signal distortion of coders with audio data compression is considered.

Keywords: audio data compression, signal distortion, digital-analogue converters, analogue-digital converters, group delay time.

УДК 539.21; 537.86

S. S. Aplesnin, A. I. Moskvin, N. I. Piskunova

SIMULATION OF THE MAGNETIC STRUCTURE UPON ORDERING EG ORBITALS TO THE QUASI-ONE- AND QUASI-TWO-DIMENSIONAL MAGNETS BY QUANTUM MONTE CARLO METOD1

The exchange mechanism influence on the electrons ordering on eg orbitals in a chain and in a two-dimensional Heisenberg model with exchange anisotropy for a S = 1/2 spin is determined by a quantum Monte Carlo method. The existence regions of a long-range quasi-one- and two-dimensional antiferromagnetic order with the special exchange topology are determined. The plateau existence region in the field dependence of magnetization, as well as the wave vector of the magnetic structure modulation with Q = p/2 in the (magnetic field - exchange alternating) plane, is determined.

Keywords: alternation of exchange, orbital and magnetic ordering, stripes.

Magnetic materials with double orbital quasidegeneration are characterized not only by spin-dependent interactions, but also by the dependence of the exchange integral on the mutual position of the orbitals. The interplay of spins with the charge-orbital ordering changes not only magnetic properties, but also transport characteristics, for example, gives rise to giant magnetoresistance in manganites. The orbital ordering forms a quasi-low-dimensional structure in KCuF3 [1] and two-dimensional antiferromagnetic in NH4CuCl3 [2] at low temperatures.

When studying the ground state and low-temperature effects in low-dimensional systems, one should remember that the quantum fluctuations in a spin system are of significant importance and, in the case of the strong interaction of spins with orbitals through the exchange interaction, it must be taken into account that the hopping integrals between the neighboring 3d ions depends on both the orbital type and on the mutual position of the sites, because the electron-density distribution is not spherically symmetric.

The interrelation between the spin order and orbital ordering is clearly illustrated in a Kugel’-Khomsky model [1] for the Hamiltonian of perovskites (eg ions at the sites of a simple cubic lattice), which is obtained from the multielectron Hamiltonian with the same exchange parameters J = 4t2/U, where t is the hopping integral and U is the Coulomb-repulsion parameter at the site between the d „ - d „ and d _ - d „ „ orbitals.

z2 z2 *2 - y2 *2 - y2

The physical origin of this peculiar situation is the strong spatial anisotropy of the eg orbital wave functions. This type of the orbital ordering provides the largest energy gain due to the quantum spin fluctuations. An orbital flip modulates strength of the neighboring exchange bond and causes the existence of the strong antiferromagnetic exchange interaction directed along the overlap of the orbitals with strong space anisotropy of the exchange interactions that are typical of the quasi-one-dimensional system under the condition of the presence of a gap in an orbital excitation spectrum.

1 This work was supported by the Russian Basic Research Foundation (RFFI-BRFFI, project No. 08-02-90031 Bela), RFFI (project

No. 09-02-00554-1).

In addition to the ferro- and antiferromagnetic ordering of the d3z2 r2 orbitals, the ordering of the d3z2 _ r2 and dx2 _ y2 orbitals is also possible with account for the fact that the hopping integral between these orbitals may substantially differ due to the strong electron correlations. Strength of the antiferromagnetic interaction is determined by a region of the overlap of the 2p-3d orbitals and depends on the orbital eg state of an electron. The overlap integral between the dx2 _ y2 andpxorbitals is Exa(l, m, n) = -^3/2 l(l2 - m2)(pds), where (pds) is the overlap integral between the d and p orbitals and (l, m, n) is the unit vector along the direction from a cation to an anion. The overlap integral between d2x2 _ z2 _ y2 and px orbitals is expressed as Ex (l, m, n) = l[n2 - (l2 +m2)/2](pds) [3]. Then, the hopping amplitude between adjacent copper ions via the p orbital along the x axis is evaluated as = Exa(1, 0, 0) Ex(-1, 0, 0)/(ep - e) [3], where ep and edare the energy levels for the d and p orbitals, respectively, and Exa and Ex p are the overlap integrals of the dx2 _ y2 and px orbitals and px and d2x2 _ z2 _y2 orbitals, which will be denoted as d^.

The account for the Coulomb interaction between electrons located both at a site and between different orbitals will also lead to modification of exchange parameters and alteration of the Jaa /Jpp relation between them. Below we consider a model with one electron (hole) on the eg orbital with an exchange interaction between them specified by a parameter. This model is applicable to perovskites containing Mn3+, Ni3+, Cu2+ and Fe2+ ions in octahedral environment.

This study is aimed to determination of the effect of the exchange interaction between eg electrons on the ordering of d and d orbitals in a Mott insulator and

2x2 - z2 - y2 x2 - y2

estimation of a region of the parameters, at which the magnetic order dimensionality changes from quasi-two-dimensional to quasi-one-dimensional antiferromagnetic. It is known that low-dimensional systems with alternating exchange imitating the spin Peierls transition are in a more stable state and have an energy gain of E/J ~ 8, where 8 = J..+1 - Jj+lj+l- The interaction of the electrons on the e orbitals with various

g

octahedron oscillation modes forms orbital ordering. The most widespread type of the orbital ordering is associated with the alternating of the d3x2 _ r2 and d3y2 _ r2 orbitals whose schematic arrangement has the form «. The exchange altenating in the spin Peierls quasi-two-dimensional antiferromagnetic gives rise to dissappearing long range order and to the formation of a singlet state. Possible the alteration of exchange along of one direction of lattice as a result of orbital ordering causes the instability AF order.

Let us determine change of the magnetic properties and the magnetic order type at the formation of orbital ordering.

1D Model. Let us consider the quasi-one-dimensional models with the symmetric alternating exchange (fig. 1, b) described by the Hamiltonian:

H = -£[( J + 5) StSt+i +JS1+1S1 + 2 +

. _ (1)

(J - 5)

H = -£[( J + 25) S,SI+l +JSI+lSI + 2 + i

+ (J -5)Si + 2Si+з +JSi+3§i + 4] -XHi SZ.

(2)

where J ~ ia, p is the exchange interactions between electrons located at the nearest d3 2 - d3 2, d 2 2 - d 2 2 d3 2 - d 2 2

3z2 3z^ x2 - y2 x2 - y2 3z2 x2 - y2

orbitals, there are J(1 ± 8), J(1 + 28), J(J < 0), 8 - a value of the alternating exchange; H is the external magnetic field. When the ratio of the hopping integrals is written as taa /t = 3/4 and t /tpp = ^/3/2 , the exchange interactions between electrons on the d and d orbitals differ by a factor of

x2 - y2 3z2 J

almost: J /JRR = 0.56, Jb /JRR = 0.75.

aa pp ’ ap pp

Л+2 Si+3 +JSi+3Si+4 ]-X HiSZ

and with the asymmetric alternating exchange(fig. 1, c):

Fig. 1. Illustration of the orbital ordering d3z2 and dx2 _ y2 for Neel AFM with homogeneous distribution of an exchange interaction (a), with the symmetric alternating exchange (eq. (1)) (b) and with the asymmetric alternating exchange (eq. (2)) (c) on a chain

As a calculation method, we take the quantum Monte Carlo method unifying two algorithms, worldlines and continuous time [4]. The continuous time world-line Monte Carlo approach is based on expansion of a statistical evolution operator exp(-H/T) by exchange interaction strength. The world-line configuration of spins flips under the action of operators of the creation and annihilation of a spin state in imaginary time and real space [4; 5]. For the spins S = 1/2 located at the sites of the chain L = 400 with the nonuniform exchange distribution and periodic boundary conditions in the Trotter direction and on the chain are used. The calculation method was described in detail in [6]. The following quantities are calculated in the framework of this method: the magnetization m = 2< Sz >, spin-spin correlation function < S z(0)S z(r) >, correlation radius, static

susceptibility c = mlH in the external magnetic field directed along the quantization axis, staggered magnetization ms, energy E, a magnetic structure factor

S= N X Г] exp(-iqrj )Sz (0)Sz (rj ).

(3)

The exchange alternating in the spin Peierls chain induces a gap in the spectrum of the triplet excitations and gives rise to a finite correlation radius in the ground state. Using the calculated dependences of the magnetization and spin correlation functions on the external magnetic field, we determine the critical fields Hc corresponding to the formation of the long-range magnetic order and the gap in the triplet-excitation spectrum. The calculations provide a linear dependence of D » 28 the gap on the exchange alternating and are in good agreement with the limiting case for one dimer, Hc = gSzJ(1 + 8) = 2J. The thermodynamic characteristics - the specific heat and susceptibility -decrease according to the exponential law below this temperature.

The ordered arrangement of the orbital pairs described by Hamiltonian (eq. (2)) gives rise to the appearance of the plateau m = 1/2 on the magnetization curve m(H) at a critical external field Hc1 as shown in fig. 2, a. The magnetic structure factor determined from the spin-spin correlation function has the main maximum and a number of satellites. The magnetic field aligns the spins of the kinks along the field and stabilizes the long-range ferromagnetic order, which does not coincide with the shortrange order calculated by chain-averaging of the spin-spin correlation function in the first coordination sphere (fig. 2, b). As the magnetic field increases, the modulation of the magnetic structure with the wave vector q = p/2 is observed near the magnetic field H ~ Hc1. The arrangement of the spins in the chain in the plateau region Hcl < H < Hc2 can be represented as . The linear

dependence m(H) for low external fields (see fig. 2, a) indicates the absence of the gap in the spectrum of triplet spin excitations and the magnetic state is a quantum spin liquid with the finite correlation radius.

The ordering of orbital pairs in the magnetic field provides three types of the magnetic order: short-range antiferromagnetic order, ferromagnetic order, and modulated ferromagnetic order with the wave vector is equal to Q = p/2, depending on the relation of the magnetic field and exchange alternating. Figure 3 shows the regions of these phases in the (magnetic field, exchange alternating) plane.

The decrease in the energy of the antiferromagnetic chain with the nonuniform periodic distribution of exchange is caused by the decrease in the effective chain length owing to a change in the correlation radius. In particular, in the limiting cases, the energies per spin in the infinite antiferromagnetic chain and dimer are E.. = -2ln 2 + 0.5 »

° chain

» - 0.443 3J and Ed = -0.75J, respectively. As a result, the energy of the alternating chain depends on the wave vector of the exchange modulation; i. e., the larger the number of the spins coupled by the strong exchange, the lower the change in the energy of the chain after the exchange alternating. The unit cell contains two and four spins in a spin Peierls magnet and in a magnet with orbital ordering, respectively. Figure 4 shows the Monte Carlo calculations of the relative energy change as a function of the exchange-alternating magnitude.

Fig. 2. Field dependences of the magnetization m (a) and spinspin correlation functions <S z(0)S z(r)> ?(b) at the distance r = 1 (1, 3, 5) and r = 31 (2, 4, 6) for 8 = 0.2 (1), 0.4 (2), and 0.6 (3). The symbol sizes correspond to the errors

Fig. 3. Phase diagram of the magnet with the orbital ordering of pairs containing the regions with the short-range antiferromagnetic order (AF), ferromagnetic order (FM), and modulated ferromagnetic order (FM mod) with the plateau m(H) = 1/2 in the field region H (1) - Hc2(2) in the (magnetic field, exchange alternating) plane. For the model with the symmetric alternating exchange (eq. (1)) are the critical fields Hc1(1, 2), and asymmetric alternating exchange (eq. (2)) are H (3, 4, 5, 6). The critical fields H (5, 6) correspond to the singlet-triplet transition in a four-spin cluster and are calculated with 8 ® 1 by exact diagonalization method in this cluster

Using the typical parameters of the hopping integrals from the cation to anion t. _ = 2 eV, t.. = 0.5 eV, and t,.

dz2 - p 7 dx2 - p 7 d(x2

- y2) - p = 1.5 eV, charge gap ep - ed = 3 eV, and U = 6 eV, we estimate the gain in the exchange energy as AEex ~ 0.05 and

0.017 eV for the bb«« ordering ofpairs of orbitals and for the dz2 - dz2 - dx2 - y2 - dx2 - y2 ordering, respectively. Owing to the competition between the Coulomb and exchange interactions between the electrons on the neighboring orbitals, a certain orbital order with the structure wave vector Q = p/2 can be induced. In magnets with narrow optical bands and weak dispersion of the optical oscillation mode and the electron-phonon coupling constant ~ g^©-1 near the band edge, the exchange mechanism possibly gives rise to change in the magnetic and orbital orders with decreasing the temperature.

5

Fig. 4. The dependence of the change normalized energy state AE(8)/£(8 = 0) = (E(8) - £(8 = 0))/E(S = 0) from the alternating exchange for the ordering of the orbitals (1) (eq. (2)) and (2) (eq. (1)), the spin Peierls 1D model (3), and power law approximation functions AE(8)/£(g = 0) = 48“ with he parameters (3) 4 = 0.67 (2) and a = 1.50(4) and (1) 4 = 0.15(1) and a = 1 (orbital ordering).

The errors correspond to the symbol sizes

Thus, in a quasi-low-dimensional magnet with one electron (hole) on the eg orbital and competing Coulomb and exchange interactions, the appearance of the orbital order with the structure wave vector Q = p/2 or the softening of the elastic oscillation mode near this vector is possible. The regions of existence of a plateau and of the modulated ferromagnet with the special exchange topology in the (magnetic field, exchange alternating) plane are determined.

2D Model. Let us consider a model with antiferromagnetic exchange alternation and a stripe structure for a spin S = 1 /2. The Hamiltonian has the following form:

H = - X [r 1+5)JaaS"+1, jSj+2, j +

i, j, a

_l_z^tttt/ott ntt _|_ ntt ott \ _|_

+K (Si+1, jSi+1, j+1 +Si+2, jSi+2, j+1) +

i Taa ca ca ,

+J Si, jSi+i, j+ (4)

+ (1 - 5)Jaa S“ jS* j+1 + Jaa S+2, jS+3, j +

+ (1 - 5)JaaS+3,jS+4, j + S+3, jS+3, j+1)] -

-X H,Sf,

where J(1 + 8), J(1 — 8), JJ < 0) are the exchange interactions between electrons located at the nearest d , - d o, ^ , -

2x2 2x2* x2 - y2

- dx2 - y2, dx2 - ^ - d2x2 orbitals; H is the external magnetic field, h = J zz — J xx(yy)/J zz is the exchange anisotropy value; K < 0 is the interchain exchange value fixed in our calculations as K/J = 1 /16. The ratio of the hopping integrals are t66 /tbb = 3 /4 and t /tpp =1/4. The exchange interactions between electrons on the d _ . and ^ . . orbitals are J /J„„ =

x2 - y2 2x2 - z2 - y2 aa pp

= 0.56 and Ja /Jpp = 0.063. This problem was solved by the described quantum Monte Carlo method. In the calculations a square lattice with L = 40; 48; 60; 72 and periodic boundary conditions in the Trotter direction and on the lattice are used. From 40 000 to 60 000 Monte Carlo steps (MCS) per site are spent to reach equilibrium and another 80 000.. .100 000 MCS are used for the averaging. The root mean square errors of the computed quantities lie in the range of 0 : 1 % to 0 : 6 %.

We consider two types of the orbital ordering that illustrated in fig. 5. The d2x2 - d2x2 orbitals overlap gives rise to the strong antiferromagnetic exchange and space anisotropy of the exchange interactions in a lattice, which differ in more than an order of magnitude. The magnetic properties of such a system are similar to those of a quasi-one-dimensional antiferromagnet. Alternation of pairs of the d2x2 - d2x2, dx2 - y2 - dx2 - y2 orbitals induces stripes along the [01] direction and alteration of the exchange interactions along the [10] direction, as one can see in fig. 5, c, d. The magnetic properties of systems with the two types of orbital ordering have been analyzed on the basis of a spin-spin correlation function, staggered magnetization m, a magnetic structure factor (eq. (3)), energy and specfic heat CkB /N = dE /dT.

Alternation of exchange along one of the lattice directions enhances the quantum fluctuations, which results in reducing spin on site in the limit of 5...9 % in the range of parameters 0 < 8 < 0 : 4 and decreasing Neel temperature. Normalized Neel temperature is well fitted by the linear function Tn(8)/Tn(0) = 1 — 0.68 for a series of anisotropy parameters (insert in fig. 6). Monte Carlo simulation of thermodynamic characteristic at larger 8 and interpolation of linear function TN(8)/TN(0) indicate to stability of AF order as compared to disordered spin state. Similar effect of reduction of spin arises from the exchange anisotropy. Staggered magnetization and Neel temperature of AF having the stripe structure rises at increasing anisotropy as shown in fig. 6. The dependence of TN(h) versus exchange anisotropy is better interpolated by a power function TN(h)/J = 1/4h16 than a logarithmic law. It differ from behavior anisotropic Heisenberg antiferromagnetic with S = 1/2 on a square lattice that reveals the logarithmic dependence TN(h)/J = 2 / ln(11/h).

We have determined the basic magnetic properties of a two-dimensional magnet with the orbital ordering presented in fig. 5, b. As a consequence of low dimensionality of a system, the quantum spin fluctuations may cause orbital correlations and the ordering of the d2x2 — d2x2 orbitals. In this case, spin interactions are quasi-one-dimensional along a chain and favor the rise of exchange energy, which is proportional to E/J = 0.44(1 + 8). Ising-like anisotropy suppresses the quantum spin fluctuations and decreases an absolute energy value from E/J = 0. 44 to 0.25 with an increase in the exchange anisotropy. With the growth of the exchange alternation along the [10] direction, energy of a 2D system

decreases due to a decrease in the interchain interaction. Figure 7 illustrates energy of a 2D magnet calculated for several lattice sizes. The intersection of energies of the 1D and 2D magnets with the arrangement of bonds in a lattice as in fig. 5, c, d determines a region of stability of the 2D magnetic ordering with a stripe structure. A phase boundary separating the region of existence of a magnets with the d2x2 - d2x2 and d2x2 - dx2 - y2 orbital ordering is given in fig. 8.

The quasi-one-dimensional AF ordering becomes stable at 8 > 0.2; h = 0, 8 >0.34; h ® 1. The ratio of exchanges in

KCuF3 is J(dx2 - y2 - dx2 - y2)/J(d2x2 - dx2 - ,2) » 0 : 78 that

corresponds to 8 ~ 0.22. Exchange anisotropy is small for the cubic crystal and according to our results the ordering of the d2x2 - d2x2 orbitals is more preferably. Orthorhombic distortion of a structure in NH4CuCl3 induces distortion of a crystal

field and the intra-atomic Coulomb interaction U between

dd

electrons on the orbital. The change in distance and angle of a bond between the nearest cations modifies the overlap of integral of the wave functions of the d2x2 and dx2 - y2 orbitals. As a result, the ratio of amplitudes of the hopping integrals may vary in the wide range 0.1 < taa /t < 1, because the hopping integral is t = E2/(ep — ed + Udd). The exchange interaction between electrons on the d2x2 - d2x2 orbitals decreases and the ratio of exchanges J(dx2 - y2 - dx2 - y2)/J(d2x2 - d2x2)

grows, that corresponds to the decreasing of an alternation parameter. As a consequence, the exchange interaction causes rearrangement of the d2x2 and dx2 - y2 orbitals and lead to strong decreasing Neel temperature in NH4CuCl3.

Thus, we study the exchange mechanism effect on the ordering of electrons on eg orbitals in a square lattice by the Heisenberg model with exchange anisotropy. The ordering of pairs of the d2x2 - d2x2, dx2 - y2 - dx2- y2 orbitals provides the two-dimensional antiferromagnetic state with a stripe structure. In the quasi-two-dimensional model with the alternating exchange the regions of the stable antiferromagnetic order in (alternating exchange-exchange anisotropy) plane are determined. The critical parameters of the alternating exchange and exchange anisotropy, at which the quasi-one-dimensional state becomes unstable are calculated.

We have analyzed the exchange mechanism of the ordering of electrons on the eg orbitals in the quasi-low-dimensional magnets using the Heisenberg model with the special exchange topology. The gain in the exchange energy for the orthogonal ordering of orbitals pairs bb«« and for the dz2 - dz2 - dx2 - y2 - dx2 - y2 ordering has estimated. We have found the exchange parameters associated with the relation of the exchange interaction of electrons on the

1 1 !K

p ' J

J(l+6)

J

I i i i J(l-o)

:ca

Fig. 5. dxl- dxl (a) and d.

be d

d10 (b) orbital ordering. Arrangement of the strong intrachain J{ 1 +8) and week interchain

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

v2 - y2

K exchanges in a square lattice; the J - exchange is related to the overlap of d2x2 A2 orbitales; J{ 1 + 8)

^v2-.v2 orbitales overlap

is related to the d2 2- d2 0 orbitals, J{ 1 - 8) is the d

,v2 - v2

0.0

0.5

П

10

n

а 6

Fig. 6. Neel temperature TN/J of the quasi-two-dimensional AF at 8 = 0.3 (1) and fitting function TN/J = h1/16/4 (solid line) as a function of exchange anisotropy (a). Insert: normalized Neel temperature TN (8) = TN (0) versus alternating exchange at h = 0.25 (1); 0.5 (2); 0.75 (3). Magnetization on the site s calculated by MC at T/TN = 0.2; L = 60 versus exchange anisotropy (b)

a

d2x2 - d2x2 and dx2 - y2 - d - y2 orbitals and exchange anisotropy, at which the ordering of pairs of the d2x2 - d2x2, d 0 „ - d , orbitals forms the two-dimensional

x2 - y2 x2 - y2

antiferromagnetic state with a stripe structure. In a quasi-low-dimensional magnet with one electron (hole) on the eg orbital and competing Coulomb and exchange interactions, the appearance of the orbital order with the structure wave vector Q = p/2 or the softening of the elastic oscillation mode near this vector is possible. The phase diagram of the modulated quase-one-dimensional ferromagnet is determined. This diagram has a plateau m(H) = 1/2 on the magnetization curve in the (magnetic field, exchange alternating) plane. The region of the existence of an quasi-two-dimensional antiferromagnetic in the plane of alternation

- anisotropy of exchange have been calculated. Quantum reduction of spin on site for quasi-two-dimensional antiferromagnetic with stripe structure versus exchange anisotropy has been estimated.

0.00 0.15 0,30 0,45 0.-50

Fig. 7. Energy of a quasi-two-dimensional antiferromagnetic E2D (see fig. 5, d) and a quasi-one-dimensional antiferromagnetic with the strong intrachain J(1 + 8) exchange normalized to the energy of the E1D antiferromagnetic with the exchange ratio K/J = 1/16 for h = 0.25; L = 40(1); 48(2); 60(3); 72(4) calculated at T/TN = 0.2 as a function of the exchange alternation. Maximum error bar denoted by square

Bibliography

1. Kugel’, K. I. The Jahn-Teller effect and magnetism: transition metal compounds / K. Kugel’, D. I. Khomsky // Sov. Phys. Usp. 1982. Vol. 25. P. 231.

2. Shiramura, W. Magnetization Plateaus in NH4CuCl3 / W. Shiramura, K. J. Takatsu, B. Kurniawan // Phys. Soc. Jpn. 1998. Vol. 67. P. 1548.

3. Slater, J. C. Simplified LCAO Method for the Periodic Potential Problem / J. C. Slater, G. F. Koster // Phys. Rev. 1954. Vol. 94. P. 1458.

4. Prokof’ev, N. V. Polaron Problem by Diagrammatic Quantum Monte Carlo / N. V. Prokof’ev, B. V. Svistunov // Phys. Rev. Lett. 1998. Vol. 81. P. 2514.

5. Аплеснин, C. C. Моделирование димерного состояния в CuGeO3 в двумерной анизотропной модели Гейзенберга с альтернированным обменом / С. С. Аплеснин // Журн. эксперим. и теорет. физики. 1997. Т. 112, № 12. С. 2184-2197.

6. Аплеснин, C. C. Неадиабатическое взаимодействие акустических фононов со спинами S = 1/2 в двумерной модели Гейзенберга / С. С. Аплеснин // Журн. эксперим. и теорет. физики. 2003. Т. 124, №9 5. С. 1080-1089.

0,50 ------------'----------1---------------------1---

0.25

0,0 0,5 1.0

Fig. 8. Phase diagram of a magnetic containing the regions of quasi-two-dimensional AF with the orbital ordering d2x2;

dx2 - y2 (under line) and quasi-one-dimensional AF with d2x2- d2x2 orbital ordering (upper line) in a plane alternation -anisotropy exchange

С. С. Аплеснин, А. И. Москвин, Н. И. Пискунова

ВЫЧИСЛЕНИЕ МАГНИТНОЙ СТРУКТУРЫ ПРИ УПОРЯДОЧЕНИИ Бс ОРБИТАЛЕЙ В КВАЗИОДНОМЕРНЫХ И КВАЗИДВУМЕРНЫХ МАГНЕТИКАХ КВАНТОВЫМ МЕТОДОМ МОНТЕ-КАРЛО

Исследовано влияние обменного механизма на упорядочение электронов на е^ орбиталях в цепочках и в двумерной модели Гейзенберга с обменной анизотропией для спина Б = 1/2 квантовым методом Монте-Карло. Определены области существования дальнего квазиодномерного и квазидвумерного антиферромагнитного порядка со специальной топологий обмена. Найдено плато намагниченности в кривой намагниченности от поля и волновой вектор модуляции магнитной структуры с Q = р /2 на плоскости «магнитное поле - альтернирование обмена».

Ключевые слова: альтернирование обмена, орбитальное и магнитное упорядочение, страйпы.

i Надоели баннеры? Вы всегда можете отключить рекламу.