Magnetoelastic effects in low-dimensional magnetic systems

We consider a new realization of magnetoelastic interactions in low-dimensional magnetic systems. We show that low-dimensional spin systems are unstable with respect to the spontaneous appearance of alternating distortions of the positions of the three-dimensional nonmagnetic atoms (ligands), that s...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Datum:2001
Hauptverfasser: Zvyagin, A.A., Zvyagina, G.A.
Format: Artikel
Sprache:English
Veröffentlicht: Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України 2001
Schriftenreihe:Физика низких температур
Schlagworte:
Online Zugang:http://dspace.nbuv.gov.ua/handle/123456789/130019
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Magnetoelastic effects in low-dimensional magnetic systems / A.A. Zvyagin , G.A. Zvyagina // Физика низких температур. — 2001. — Т. 27, № 4. — С. 406-411. — Бібліогр.: 27 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id irk-123456789-130019
record_format dspace
spelling irk-123456789-1300192018-02-05T03:02:48Z Magnetoelastic effects in low-dimensional magnetic systems Zvyagin, A.A. Zvyagina, G.A. Низкотемпеpатуpная магнитостpикция магнетиков и свеpхпpоводников We consider a new realization of magnetoelastic interactions in low-dimensional magnetic systems. We show that low-dimensional spin systems are unstable with respect to the spontaneous appearance of alternating distortions of the positions of the three-dimensional nonmagnetic atoms (ligands), that surround the magnetic ions. Those distortions are supplemented by the spontaneous onset of alternating effective g factors of the magnetic ions in the phase with short-range interactions. We discuss the possibility of observing the effect in an uniform external magnetic field, which in the situation considered produces both magnetization and staggered magnetization of the magnetic subsystem. The connection of the proposed theory with recent experiments on effectively low-dimensional magnetic systems (organic spin chains, heavy-fermion compounds, rare-earth molybdates) is discussed. 2001 Article Magnetoelastic effects in low-dimensional magnetic systems / A.A. Zvyagin , G.A. Zvyagina // Физика низких температур. — 2001. — Т. 27, № 4. — С. 406-411. — Бібліогр.: 27 назв. — англ. 0132-6414 PACS: 75.10.Jm, 71.70.Fk, 71.70.Ej http://dspace.nbuv.gov.ua/handle/123456789/130019 en Физика низких температур Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
language English
topic Низкотемпеpатуpная магнитостpикция магнетиков и свеpхпpоводников
Низкотемпеpатуpная магнитостpикция магнетиков и свеpхпpоводников
spellingShingle Низкотемпеpатуpная магнитостpикция магнетиков и свеpхпpоводников
Низкотемпеpатуpная магнитостpикция магнетиков и свеpхпpоводников
Zvyagin, A.A.
Zvyagina, G.A.
Magnetoelastic effects in low-dimensional magnetic systems
Физика низких температур
description We consider a new realization of magnetoelastic interactions in low-dimensional magnetic systems. We show that low-dimensional spin systems are unstable with respect to the spontaneous appearance of alternating distortions of the positions of the three-dimensional nonmagnetic atoms (ligands), that surround the magnetic ions. Those distortions are supplemented by the spontaneous onset of alternating effective g factors of the magnetic ions in the phase with short-range interactions. We discuss the possibility of observing the effect in an uniform external magnetic field, which in the situation considered produces both magnetization and staggered magnetization of the magnetic subsystem. The connection of the proposed theory with recent experiments on effectively low-dimensional magnetic systems (organic spin chains, heavy-fermion compounds, rare-earth molybdates) is discussed.
format Article
author Zvyagin, A.A.
Zvyagina, G.A.
author_facet Zvyagin, A.A.
Zvyagina, G.A.
author_sort Zvyagin, A.A.
title Magnetoelastic effects in low-dimensional magnetic systems
title_short Magnetoelastic effects in low-dimensional magnetic systems
title_full Magnetoelastic effects in low-dimensional magnetic systems
title_fullStr Magnetoelastic effects in low-dimensional magnetic systems
title_full_unstemmed Magnetoelastic effects in low-dimensional magnetic systems
title_sort magnetoelastic effects in low-dimensional magnetic systems
publisher Фізико-технічний інститут низьких температур ім. Б.І. Вєркіна НАН України
publishDate 2001
topic_facet Низкотемпеpатуpная магнитостpикция магнетиков и свеpхпpоводников
url http://dspace.nbuv.gov.ua/handle/123456789/130019
citation_txt Magnetoelastic effects in low-dimensional magnetic systems / A.A. Zvyagin , G.A. Zvyagina // Физика низких температур. — 2001. — Т. 27, № 4. — С. 406-411. — Бібліогр.: 27 назв. — англ.
series Физика низких температур
work_keys_str_mv AT zvyaginaa magnetoelasticeffectsinlowdimensionalmagneticsystems
AT zvyaginaga magnetoelasticeffectsinlowdimensionalmagneticsystems
first_indexed 2025-07-09T12:42:29Z
last_indexed 2025-07-09T12:42:29Z
_version_ 1837173259774197760
fulltext Fizika Nizkikh Temperatur, 2001, v. 27, No. 4, p. 406–411Z vya gin A. A . and Z vya gina G. A .Ma gne toelast ic e ffec ts in low-dim ensiona l m agn etic system sZvyag in A . A. a nd Zvyag in a G. A.Mag neto elastic eff ects in low-dimen sio nal mag net ic syst ems Magnetoelastic effects in low-dimensional magnetic systems A. A. Zvyagin1,2 and G. A. Zvyagina2 1 Max Planck Institut fur Physik Komplexer Systeme, No..thnitzer Str, 38, D-01187, Dresden, Germany 2 B. Verkin Institute for Low Temperature Physics and Engineering of the National Academy of Sciences of Ukraine 47 Lenin Ave., Kharkov 61164, Ukraine Received July 25, 2000 We consider a new realization of magnetoelastic interactions in low-dimensional magnetic systems. We show that low-dimensional spin systems are unstable with respect to the spontaneous appearance of alternating distortions of the positions of the three-dimensional nonmagnetic atoms (ligands), that surround the magnetic ions. Those distortions are supplemented by the spontaneous onset of alternating effective g factors of the magnetic ions in the phase with short-range interactions. We discuss the possibility of observing the effect in an uniform external magnetic field, which in the situation considered produces both magnetization and staggered magnetization of the magnetic subsystem. The connection of the proposed theory with recent experiments on effectively low-dimensional magnetic systems (organic spin chains, heavy-fermion compounds, rare-earth molybdates) is discussed. PACS: 75.10.Jm, 71.70.Fk, 71.70.Ej Interest in electron systems with substantial coup- ling between the charge, spin, and orbital degrees of freedom of the electronic and elastic subsystems of a crystal has grown considerably during the last decade. The prime examples of the manifestation of such cooperative effects are the phenomena of col- losal magnetoresistance of manganites, non-Fermi- liquid behavior of some heavy-fermion compounds, spin-Peierls and charge ordering behaviors in inor- ganic systems and unconventional superconducti- vity [1]. The Jahn-Teller effect [2] is probably the oldest known manifestation of such a coupling. Here the degeneracy of the orbital states of a molecule is removed by the deformation of the latter. The cooperative Jahn–Teller effect reveals itself in a structural phase transition. It has been observed in a number of compounds [3]. In the spin-Peierls transition also the degeneracy of the electronic (spin, not orbital) subsystem of a one- dimensional (1D) spin chain is removed due to the coupling with a longitudinal phonon of the 3D crystal lattice. A spin gap is opened for the low- lying spin excitation. On the other hand, the corre- sponding phonon mode possesses softening (Kohn anomaly). Some of magnetic compounds with the essential coupling between spin, orbital, and elastic subsystems manifest paramagnetic spin behavior with two essentially inequivalent magnetic centers at low temperatures (which are higher, though, than the temperature of the phase transition to magnetically ordered 3D state). Often the latter has not been observed [4–9]. For higher tempera- tures the inequivalence between the two magnetic centers smears out. The presence of inequivalent magnetic centers in low-dimensional quantum spin systems is usually connected with slightly different local surrounding of two types of magnetic ions and involves either staggered g factors of the magnetic ions [4] (another case pertains to two anisotropic g tensors canted with respect to the principal axis [8] or the Dzyaloshinskii–Moriya coupling in crystals without mirror magnetic symmetry [10]. Low-lying spin excitations are gapless for these low-dimen- sional systems in the absence of an external mag- netic field. However,in an external field the low- temperature specific heat of some of those systems reveals the emergence of a spin gap induced by an external magnetic field [8,11,12]. To explain that gapped behavior a 1D spin subsystem with staggered effective g factors has been sta- died [5,13,14]. The low-temperature, low-field electron spin resonance (ESR) [4,6–9] exhibited two inequivalent magnetic centers. On the other hand, the higher-field ESR treatment reported only © A. A. Zvyagin and G. A. Zvyagina, 2001 a single magnetic center [15]. Optical (Raman and infrared) measurements [16] have revealed the tendency for softening of some lattice modes. Ultra- sonic measurements also revealed anomalous beha- viors of the sound velocities for some phonon modes in such crystals [17,18]. In this work we propose a mechanism for expla- nation of the observed anomalous behavior of the spin, orbital, and elastic subsystems of that group of magnetic low-dimensional systems. Namely, we point out that the dimerization of the effective g factors in a quantum antiferromagnetic spin chain can be caused by the coupling of the spin and orbital subsystems to the elastic subsystem, fol- lowed by a distortion of the local surrounding of the magnetic ions. In paramagnets the value of the g factor differs from 2. The difference stems from the effect of the crystalline electric field of nonmag- netic ligands due to the presence of a spin-orbit interaction. The essence of our study is following. Let the configuration of ligands, which surround two neighboring magnetic ions along the chain direction, possesses small shifts (distortions) of op- posite signs (antiferrodistortions [8]). Such distor- tions will immediately produce a change of the crystalline electric fields of the ligands. Hence, the orbital moment of the magnetic ions will be alternately affected. Then the spin–orbit interaction yields, first, two different values of the effective g factors of the magnetic ions, and second, a stag- gered Dzyaloshinskii—Moriya interaction (because of the odd magnetic symmetry with respect to the principal axis). The effect of the latter can be also transferred to the effective staggered g factor by means of an alternating rotation of spins. Hence, the staggered distortion of the nonmagnetic sur- roundings (ligands) of the magnetic ions produces two inequivalent magnetic centers for the spin sub- system. The elastic subsystem loses energy due to the alternating distortions of ligands, while the electron subsystem realizes an energy gain. The steady-state configuration of the total crystal is determined by the competition between those two processes. From the magnetic standpoint that effect is collective because the spins interact with each other. On the other hand, the inequivalence of magnetic centers can be seen only in nonzero mag- netic field, so the effect is field-induced. To fix the stage, consider the low-temperature behavior of a quasi-low-dimensional magnetic crys- tal. The spin–spin coupling along one direction of this crystal is larger than the interactions along the other crystallographic directions. Suppose this quasi-1D spin subsystem can be described by a Hamiltonian of a gapless spin-1/2 antiferromag- netic (AF) chain. The fact that the low-lying exci- tation of the 1D AF spin chain is gapless implies that the system is critical, i.e., the ground-state correlation functions decay in a power-law manner. Suppose also that for the reasons explained above, some components of the effective g factors of the spins are alternating, i.e., g1,2 = g(1 ± δ), where δ << 1 is proportional to the small alternating distortions of the local environment. The parameter δ is not fixed, but has to be determined from the conditions of the steady state. It is equivalent to the mean field approximation for phonons. It is justi- fied for the 3D elastic subsystem. The 1D spin–spin interactions will be taken into account nonper- turbatively. In the case of alternating Dzyaloshin- skii—Moriya interactions present in the spin chain the alternating effective g factors are g1,2 = g × × cos ϑ(1 ± tan ϑ), where the angle ϑ originates from, e.g., the alternating components of the crys- talline electric field of the ligands [14]. For rare- earth molybdates, where the g tensors for two cen- ters are canted with respect to the main axis, this means that g1,2 = g cos α cos θ (1 ± tan α tan θ) [8], where α and θ are the angles determining the canting of the g tensors and field direction, respec- tively. The application of a field parallel to the crystallographic axes does not change the distor- tion. For rare-earth ions (with a Kramers- or non- Kramers-doublet crystalline field ground state) for sufficiently low temperatures T compared with the crystalline-field splitting D, one can use the two lowest levels of the ion as an effective spin-1/2 with magnetically anisotropic behavior. It turns out that no change of the g factors can be produced by distortions of ligands for the ground-state doublets. They do not have orbital degrees of freedom and are therefore unaffected by distortions. However, the excited crystalline electric field states are effec- tively included. This is why a distortion of the ligands can produce a change of the g factors, through the off-diagonal matrix elements of the relevant electric-multipole operators. This means that the conditions for the applicability of our effective spin-1/2 description for rare-earth com- pounds are: (i) T << D, i.e., there is no thermal population of the crystalline electric field-excited levels; (ii) off-diagonal matrix elements (between the ground-state doublet and excited levels) of the Zeeman term must be small compared with D; (iii) the matrix elements between the ground-state doublet and crystalline electric field-exited levels of the operator describing the distortion must not be too small compared to D (otherwise distortions Magnetoelastic effects in low-dimensional magnetic systems Fizika Nizkikh Temperatur, 2001, v. 27, No. 4 407 would produce a negligible change of the ground- state g factors); (iv) distortions of the ligands must not affect the exchange constant J; see below (with J << D). The Hamiltonian can be written as � sp = J ∑ j (S j S j+1 + (∆ − 1)Sj zSj+1 z ) − − µB H ∑ j (g1 S2j z + g2 S2j+1 z ) , (1) where 0 ≤ ∆ ≤ 1 is the parameter of the «easy- plane» magnetic anisotropy, H is an external mag- netic field, and µB is the Bohr magneton. It is easy to show that because of nonzero δ the application of an external field to this system yields two effects: homogeneous magnetization (mz) and staggered magnetization. The former changes the ground-state filling of the Dirac sea (H plays the role of the Fermi energy). For large enough values of the field H the spin subsystem possesses trimerization (for the Dirac sea of spinons being one-third filled), quadrimerization (for one-fourth filling), etc. For those values of the field 2nkF = π, with n being an integer and 2kF = π(1 + 2mz), i.e., we expect the series of transformations towards an inhomogeneous (incommensurate) magnetic structure due to the umklapp processes rather than the dimeriza- tion [19]. We restrict ourselves to small enough values of H here. The staggered magnetization is the relevant perturbation from the renormalization group (RG) viewpoint. It produces the gap for low-lying magnetic excitations of the system (spi- nons). We can calculate the response of our (criti- cal) spin chain to the relevant perturbation (stag- gered field). We perform such a calculation in a RG framework. Our study shows that the exponents are non-integer in general, in contrast to simple pertur- bation or mean-field theories [8,16], which of course are not legitimate for the low-dimensional quantum spin systems. An application of scaling relations provides a simple tool to understand some essential aspects of the behavior of a critical chain under a relevant perturbation. Recall that response of the free energy fcl and the correlation function ξcl of a classical critical d-dimensional system per- turbed by a relevant operator δ� ′ with RG eigen- value ν−1 > 0 is ∆fcl ∝ δdν , ξcl ∝ δ−ν . (2) A quantum critical d-dimensional system (which in our case is the spin 1D subsystem of the crystal) formally behaves in the scaling regime equivalently to a (d + z)-dimensional classical system, where z is the dynamical critical exponent. Hence, the ground- state energy and the gap of the low-lying spin excitations of the d-dimensional quantum critical system are formally proportional to the free ener- gy and to the inverse correlation function of the (d + z)-dimensional classical critical system, respec- tively. The RG eigenvalue ν−1 is related to the scaling dimension x of the particular operator by x + ν−1 = d + z. For the (conformally invariant) AF spin chain we have d = z = 1, i.e., ν = (2 − x)−1. Hence the renormalization of the ground-state ener- gy per site of the quantum critical chain and the low-lying spin excitation gap (which is equal to zero at the unperturbed point) to the staggered magnetic field are ∆Eq ∝ − (gδµ B H)2/(2−x e ), ∆ ∝ (gδµB H)1/(2−x e ) , (3) respectively, where xe is the minimal scaling expo- nent for energy—energy correlations. Here we have ignored logarithmic corrections. They are present due to the marginal operators in the RG sense and essential for the SU(2)-symmetric case. To find the scaling dimension for our critical spin chain we use the results of the conformal field theory (CFT). According to the CFT approach [20] the asymptotics of the correlation functions of pri- mary fields in the ground state are known to be 〈φ∆±(r, t)φ∆±(0, 0)〉 = exp (2iDPF r)(r − ivF t)−2∆+ × × (r + ivF t)−2∆− , where vF and PF are the Fermi velocity and the Fermi momentum, respectively [a (half-)integer D measures the momentum of the primary field in units of the Fermi momentum]. The scaling dimension and spin for each primary field are determined by xφ = ∆+ + ∆− and sφ = ∆+ − ∆−. The parameters ∆± can be calculated according to the finite-size analysis of the low-energy physics of the critical spin chain. Combining all the effects, we can write for the ground-state energy of the spin subsystem with nonzero δ Esp = − vF    gδµ B H vF    2/(2−x e ) , (4) where for our AF spin chain vF = πJ/2 × × √1 − ∆2 cos−1 ∆ is the Fermi velocity of spinons for H = 0. We point out that the scaling approach is only valid in the vicinity of the critical point, i.e., in principle the values of the magnetic field and δ are small. Note that the Fermi velocity monotonically decreases with the growth of H and becomes zero at the spin-saturation point (for A. A. Zvyagin and G. A. Zvyagina 408 Fizika Nizkikh Temperatur, 2001, v. 27, No. 4 H ≥ Hc ∼ J∆/4gµB). The exponent is equal to xe = π/2(π − cos−1 ∆) for H = 0 (it increases monot- onically with the increase of H and becomes 1 at H = Hc) [21]. Obviously Esp < 0, and hence one has an energy gain due to nonzero δ. On the other hand, in the lowest order in δ the elastic subsystem loses an energy proportional to Cδ2/2, where C is the elastic constant. Whether the ground-state steady-state configuration corresponds to zero or nonzero δ depends on the scaling exponent xe . For xe < 1 the ground-state steady-state configuration of the spin and elastic subsystems corresponds to nonzero δ, and hence, to two inequivalent spin centers (with two different g factors) and to the nonzero alternating distortions (antiferrodistor- tions) of the local environment of the magnetic ions (ligands). (Here we take into account that |δ| < 1.) The steady-state δ is equal to δ0 =    vF C    2−x e 2(1−x e )    gµ B H vF    1 1−x e . (5) Clearly δ0 = 0 for H = 0 and for H ≥ Hc . This is the main result of our report. For high temperatures (much higher than J), the steady-state δ is zero, naturally. Hence, there has to be a phase transition between the low-temperature phase with nonzero staggered magnetization, gapped low-lying spin ex- citations, and nonzero alternating distortions of the ligands, surrounding the magnetic ions, and the high-temperature phase in which the staggered mag- netization is zero, low-lying spin excitations are gapless and there are no distortions of the ligands. In fact we have considered the effect under a tacit assumption that the initial frequencies ω0 of the phonons which are coupled to the spin chain are small ω0 << J (adiabatic approximation). It is in principle possible to calculate the effect more pre- cisely for any ω0 [22,23]: there the conditions for the steady-state configuration of the total system correspond to the cross section of the respective phonon modes ω(q) (at q = 2kF) connected to the alternating distortions (antiferrodistortions) of li- gands with the magnetic modes [22,23]. Let us consider the interesting (RG-marginal with xe = 1) limiting case, namely ∆ = 0. That case corresponds to the XX spin-1/2 chain with alter- nating g factors. The Hamiltonian of the latter can be exactly mapped by means of the nonlocal Jor- dan—Wigner transformation on the Hamiltonian of the lattice noninteracting fermions [24]. The ground-state energy of the spin subsystem and the elastic subsystem can be written in the form Egs = Cδ2 2 − gµ B H    1 − λ c π    − − 1 π ∫ 0 λ c dλ √(δgµ B H)2 + [J cos (λ/2)]2 , (6) where λc = cos−1 [2(1 − δ2)(gµB H/J)2 − 1]. Mini- mizing the ground-state energy with respect to δ we obtain two possible solutions. The first one, δ = 0, corresponds to zero distortions. As to the second one, it satisfies the equation πC √J2 + (δgµ B H)2 (2gµ B H)2 = F    λ c 2 , k    , (7) where F(λc /2, k) is the incomplete elliptic integral of the first kind with k−1 = √1 + (δgµB H/J)2 . No- tice that for H ≥ Hc = J/gµB √1 − δ2 one has λc = π, and the integral is zero (hence for H ≥ Hc there is only one δ0 = 0 solution in the ground state, as should be the case). We emphasize again that we are studying the effect of sufficiently weak mag- netic fields, for which, on the one hand, the scaling approach is valid, and, on the other hand, we restrict consideration to the dimerization only, i.e., the field is not so strong as to produce a one-third filling of the Dirac sea, etc., or a series of trans- formations into inhomogeneous magnetic phases. Equation (7) can be re-written in the form √1 − (1 − δ2)(gµ B H/J)2 = sn (u, k) , (8) where u = πC √J2 + (δgµB H)2/(2gµB H)2 and sn (u, k) is the Jacobi elliptic sine function. This equation can be solved numerically. For k ∼ 1 (which corresponds to small δ), the asymptotic behavior of the elliptic function is known to be sn (u, k) ≈ tanh u + (1 − k2)(sinh u cosh u − u)/4 cosh2 u. (For large k one has sn (u, k) ≈ sin u − k2 × × cos u(u − sin u cos u)/4.) It is possible to write an analytical asymptotical expression for the gro- und-state steady-state nonzero δ0 as δ0 ≈ (2J/gµ B H) × ×      A tanh u0 − A2 2 − A tanh 2u0 − (Au0/cosh2 u0)      1/2 , (9) where A = √1 − (gµB H/J)2 and u0 = πCJ × × (2gµB H)−2. It is clear that nonzero δ0 can appear Magnetoelastic effects in low-dimensional magnetic systems Fizika Nizkikh Temperatur, 2001, v. 27, No. 4 409 only for tanh u > A. For nonzero temperatures we calculate the free energy of our system (assuming that the Debye energy of the phonons is large, and, therefore the elastic subsystem effectively remains in the ground state) F = Cδ2 2 − T π ∑ ± ∫ 0 π dλ ln    2 cosh    ε±(λ) 2T       , (10) where ε±(λ) = gµB H ± √(δgµB H)2 + [J cos (λ/2)]2 . Then the minimization of the free energy with respect to δ yields Cδ = 1 2π ∑ ± ∫ 0 π dλ tanh ε±(λ) 2T ∂ε±(λ) ∂δ . (11) For T >> J we can replace the hyperbolic tangent by its argument, and one clearly sees that there exists only one solution to Eq. (11): δ0 = 0. The numerical analysis shows that δ0 ≠ 0 exists at low temperatures. The critical temperature Tc of the transition between the low-temperature (effectively two-center) phase and the uniform high-tempera- ture phase is obtained by setting δ0 = 0 for the second order phase transition, or the free energies of the phases with δ0 = 0 and δ0 ≠ 0 for the first-order one. The critical temperature can be estimated (with the main input to the integral given by the divergent denominator, i.e., at the van Hove singu- larities) as Tc ∼ (J − gµ B H)/ln (πCJ/2g2µ B 2 H2) . (12) The analysis shows that δ0 = 0 pertains to the mini- mum of the free energy for T > Tc and to the maximum for T < Tc . Hence the phase transition is of the second order. We again emphasize that the phase transformations studied here are magnetic field-induced; they are absent without the latter. We point out also that the order of the phase transformation in the general case is connected with the orientation of the direction of the magnetic field with respect to the crystal axes and on the mutual ratios of the exchange constant, anisotropy, mag- netic field, and elastic constant. It turns out that (alternating) distortions of the ligands, which result in the inequivalence of the effective g factors of neighboring spins in quasi-1D chains, can be supplemented in the real compounds by alternating distortions of the magnetic ions themselves (which leads to the spin-Peierls insta- bilities), as has been observed, e.g., in CuGeO3 [25]. In such a situation the two effects are additive in the formation of a spin gap for low-lying spin excitations. The values of the effects are naturally different. Although the effect studied in this work manifests itself only in some domain of values of the magnetic field, the spin-Peierls period doubling is not caused by the magnetic field. The critical temperature is determined from the condi- tion of which of the instabilities becomes mani- fested first. We suppose that the cooperative effects similar to the one studied in this paper have been observed in some rare-earth molybdates. There at low-tem- perature ESR studies that have observed two ine- quivalent magnetic centers, (see, e.g., Refs. 6–9), while for higher temperatures [7,8] or higher mag- netic fields (higher frequencies) [15] only one mag- netic center was seen. The low-temperature specific heat exhibits features at temperatures higher than the transition to the magnetically ordered phase (which were hidden, though, due to the presence of a Schottky anomaly) [8]. Ultrasonic measurements have observed anomalies in the low-temperature dependences of the velocities of sound for acous- tic [17] and optical branches [18] of the phonon spectra. Optical and magnetic investigations have observed a cooperative effect (in nonzero magnetic field), that was classified as being of the Jahn—Tel- ler type [6,16,26]. Note that the transition to the magnetically ordered state was observed at much lower temperatures (see, e.g., Ref. 8, 18, 26), hence the cooperative effect we have discussed is observed in the paramagnetic (from the magnetic viewpoint) phase. It turns out that magnetic and elastic charac- teristics of some of those systems manifested the mentioned features only for nonzero external mag- netic fields, with the order of the transition (which of the derivatives of the thermodynamic potential possesses singular behavior) being determined by the direction of the field, while for zero field for some compounds there were no observations of the Jahn—Teller-like cooperative effect [27]. Summarizing, we have studied a new realization of magnetoelastic interactions in low-dimensional magnetic systems. We have shown that low dimen- sional spin systems are unstable with respect to the spontaneous appearance of alternating distortions of the positions of the three-dimensional nonmagnetic atoms (ligands) that surround the magnetic ions. Those distortions are supplemented by the sponta- neous onset of alternating effective g factors of magnetic ions in the phase with short-range interac- tions. We have discussed the possibility of obser- A. A. Zvyagin and G. A. Zvyagina 410 Fizika Nizkikh Temperatur, 2001, v. 27, No. 4 ving the effect in a uniform external magnetic field, which in the situation considered produces both magnetization and staggered magnetization of the magnetic subsystem. We suppose that the effects, investigated theoretically in this work, have pro- bably been observed in the low-temperature experi- ments in some quasi-low-dimensional magnetic com- pounds, in which the spin, charge, and orbital characteristics of the electron subsystem are sub- stantially coupled to the elastic subsystem. 1. Colossal Magnetoresistance, Charge Ordering and Related Properties of Manganese Oxides, C. N. R. Rao and B. Raveau (eds.), World Scientific, Singapore (1998); Co- lossal Magnetoresistance Oxides, Y. Tokura (ed.), Gordon and Breach, Monographs in Cond. Matt. Science (1999); H. V. Lo..hneysen, J. Magn. Magn. Mater. 200, 532 (1999); M. Hase, I. Terasaki, and K. Uchinokura, Phys. Rev. Lett. 70, 3651 (1993); M. Isobe and Y. Ueda, J. Phys. Soc. Jpn. 65, 1178 (1996); Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, T. Fujita, J. G. Bednorz, and F. Lichtenberg, Nature (London) 372, 532 (1994); T. Ishigiro, K. Yamaji, and G. Saito, Organic Supercon- ductors, Springer-Verlag, Berlin (1998). 2. H. A. Jahn and E. Teller, Proc. Roy. Soc. A161, 220 (1937). 3. K. I. Kugel’ and D. I. Khomskii, Usp. Fiz. Nauk 136, 621 (1982) [Sov. Phys. Usp. 25, 231 (1982)]. 4. M. Date, H. Yamazaki, M. Motokawa, and S. Tazawa, Prog. Theor. Phys. Suppl. 46, 194 (1978); K. Oshima, K. Okuda, and M. Date, J. Phys. Soc. Jpn. 41, 475 (1976); ibid. 44, 757 (1978); K. Takeda, Y. Yoshino, K. Matsumoto, and T. Haseda, ibid. 49, 162 (1980). 5. A. Ochiai, T. Suzuki, and T. Kasuya, J. Phys. Soc. Jpn. 59, 4129 (1990); P. Bonville, A. Ochiai, T. Suzuki, and E. Vincent, J. Phys. (Paris) 4, 595 (1994); M. Kohgi, K. Iwasa, J.-M. Mignot, A. Ochiai, and T. Suzuki, Phys. Rev. B56, R11388 (1997); M. Kohgi, K. Iwasa, J.-M. Mignot, N. Pyka, A. Ochiai, H. Aoki, and T. Suzuki, Physica B230–232, 638 (1997); Physica B259–261, 269 (1999); M. Oshikawa, K. Ueda, H. Aoki, A. Ochiai, and M. Kohgi, J. Phys. Soc. Jpn. 68, 3181 (1999). 6. T. S. Stetsenko, A. M. Pshisukha, S. D. Elchaninova, and A. I. Zvyagin, Opt. Spectrosc. 34, 227 (1973); A. I. Zvy- agin, T. S. Stetsenko, V. G. Yurko, and R. A. Vaishnoras, JETP Lett. 17, 135 (1973). 7. A. M. Pshisukha, A. S. Chernyi, and A. I. Zvyagin, Fiz. Nizk. Temp. 1, 473 (1975) [Sov. J. Low Temp. Phys. 1, 233 (1975)]; A. M. Pshisukha, A. I. Zvyagin, and A. S. Chernyi, ibid. 2, 339 (1976) [ibid. 2, 18 (1976)]. 8. A. H. Cooke, M. M. Davidson, N. J. England, M. J. M. Leask, J. B. Lowry, A. C. Tropper, and M. R. Wells, J. Phys. C9, L573 (1976); M. J. M. Leask, A. C. Tropper, and M. R. Wells, ibid. 14, 3481 (1981). 9. M. I. Kobets, V. A. Pashchenko, and E. N. Khatsko, Fiz. Nizk. Temp. 26, 370 (2000) [Low Temp. Phys. 26, 270 (2000)]. 10. A. A. Zvyagin, Fiz. Nizk. Temp. 21, 825 (1995) [Low Temp. Phys. 21, 636 (1995)]. 11. D. C. Dender, P. R. Hammar, D. H. Reich, C. Broholm, and G. Aeppli, Phys. Rev. Lett. 79, 1750 (1997). 12. M. Ko..ppen, M. Lang, R. Helfrich, F. Steglich, P. Thal- meier, B. Schmidt, B. Wand, D. Pankert, H. Benner, H. Aoki, and A. Ochiai, Phys. Rev. Lett. 82, 4548 (1999); F. Steglich, P. Gegenwart, R. Helfrich, C. Langhammer, P. Hellmann, L. Donnevert, C. Geibel, M. Lang, G. Sparn, W. Assmus, G. R. Stewart, and A. Ochiai, Z. Phys. B103, 235 (1997); F. Steglich, M. Ko..ppen, P. Gegenwart, T. Cichorek, B. Wand, M. Lang, P. Thalmeier, B. Schmidt, H. Aoki, and A. Ochiai, Acta Phys. Polon. A97, 1 (2000). 13. M. Oshikawa and I. Affleck, Phys. Rev. Lett. 79, 2883 (1997). 14. H. Shiba, K. Ueda, and O. Sakai, Preprint ISSP (1999). 15. A. G. Anders, S. V. Volotskii, and O. E. Zubkov, Fiz. Nizk. Temp. 19, 731 (1993) [Low Temp. Phys. 19, 526 (1993)]. 16. D. Mihailovic′, J. F. Ryan, and M. C. K. Wiltshire, J. Phys. C20, 3047 (1987); D. Mihailovic′ and J. F. Ryan, ibid. 20, 3063 (1987). 17. I. M. Vitebskii, S. V. Zherlitsyn, A. I. Zvyagin, A. A. Stepanov, and V. D. Fil, Fiz. Nizk. Temp. 12, 1108 (1986) [Sov. J. Low Temp. Phys. 12, 626 (1986)]. 18. G. A. Zvyagina and A. A. Zvyagin, Fiz. Nizk. Temp. 26, 482 (2000) [Low Temp. Phys. 26, 354 (2000)]. 19. L. N. Bulaevskii, Zh. Eksp. Teor. Fiz. 43, 968 (1962) [Sov. Phys. JETP 16, 685 (1963)]; M. C. Cross, Phys. Rev. B20, 4606 (1979). 20. A. A. Belavin, A. M. Polyakov, and A. B. Zamolodchikov, Nucl. Phys. B241, 333 (1984). 21. See, e.g. V. E. Korepin, N. M. Bogoliubov, and A. G. Izergin, Quantum Inverse Scattering Method and Correla- tion Functions, Cambridge University Press (1993) and references therein. 22. M. C. Cross and D. S. Fisher, Phys. Rev. B19, 402 (1979). 23. See, e.g., R. J. Bursill, R. H. McKenzie, and C. J. Hamer, Phys. Rev. Lett. 83, 408 (1999). 24. V. M. Kontorovich and V. M. Tsukernik, Zh. Eksp. Teor. Fiz. 53, 1167 (1967) [Sov. Phys. JETP 26, 687 (1968)]. 25. A. A. Gippius, E. N. Morozova, D. F. Khozeev, A. N. Vasil’ev, M. Baenitz, G. Dhalenne, and A. Revcolevschi, J. Phys. Condensed Matt. 9, L71 (2000). 26. V. I. Kut’ko, V. A. Pashchenko, and M. I. Kobets, Fiz. Nizk. Temp. 19, 1354 (1993) [Low Temp. Phys. 19, 962 (1993)]. 27. V. I. Kut’ko, Yu. M. Kharchenko, A. A. Stepanov, and N. M. Nesterenko, Fiz. Nizk. Temp. 20, 361 (1994) [Low Temp. Phys. 20, 288 (1994)]. Magnetoelastic effects in low-dimensional magnetic systems Fizika Nizkikh Temperatur, 2001, v. 27, No. 4 411