arxiv:cond-mat/ v1 [cond-mat.str-el] 17 Jan 1997

Size: px

Start display at page:

Download "arxiv:cond-mat/ v1 [cond-mat.str-el] 17 Jan 1997"

Simon Austen Shaw
5 years ago
Views:

1 Field-induced Commensurate-Incommensurate phase transition in a Dzyaloshinskii-Moriya spiral antiferromagnet. A. Zheludev, S. Maslov, and G. Shirane Department of Physics, Brookhaven National Laboratory, Upton, New York arxiv:cond-mat/ v1 [cond-mat.str-el] 17 Jan 1997 Y. Sasago, N. Koide, and K. Uchinokura Department of Applied Physics, The University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113, Japan (February 24, 2018) We report an observation of a commensurate-incommensurate phase transition in a Dzyaloshinskii- Moriya spiral magnet Ba 2CuGe 2O 7. The transition is induced by applying a magnetic field in the plane of spin rotation. In this experiment we have direct control over the strength of the commensurate potential, while the preferred incommensurate period of the spin system remains unchanged. Experimental results for the period of the soliton lattice and bulk magnetization as a function of external magnetic field are in quantitative agreement with theory. Studies of Commensurate-Incommensurate (CI) phase transitions have a long history, dating back to pioneering works of Frenkel and Kontorova [1] and Frank and van der Merwe. [2]. Since then CI transitions were discovered and studied in a number of such seemingly unrelated systems as noble gas monolayers adsorbed on graphite surface [3,4], charge density wave materials [5,6], and rareearth magnets [7] (For comprehensive reviews see for example Ref. [8]). As a rule, CI transitions result from a competition between two distinct terms in the Hamiltonian that have different built-in spatial periodicities and are often referred to as potential and elastic energy, respectively. The potential energy by definition favors a structure commensurate with the crystal lattice. The elastic term is intrinsic to the system where the transition occurs, and has a different natural built-in period. In many known realizations of CI, such as adsorbed gas monolayers, it is the period set by the elastic term that can be varied in an experiment to drive the transition, whereas both the strength and the period of the potential remain constant. In other systems, such as rareearth magnets, both the elastic term (exchange coupling between spins) and the potential (magnetic anisotropy) can be changed, but only indirectly, by varying the temperature. Fromtheverystartitwasclearthatinitspurestforma CI transition may be driven by the change of the strength of the potential alone, with the two built-in periods remaining constant. An elegant experimental realization of this type of CI was first proposed by Dzyaloshinskii [9], who considered an incommensurate spiral magnetic structure in a magnetic field applied in the plane of rotation of spins. In this model incommensurability is intrinsic to the spin system and results from spin-spin interactions. The role of the potential is played solely by an external magnetic field H, that favors a commensurate spin-flop state. In a real magnetic material with this type of CI transition an experimentalist would have a convenient handle on the strength of the potential term, adjusting it by simply changing the magnetic field. Such systems are not easy to find. In most spiral magnets, e.g., cubic MnSi [10] and FeGe [11], even a small field is sufficient to realign the spin plane perpendicular to the field direction. While RbMnBr 3 [12] may actually be one material where this does not occur, to our knowledge, the most crucial quantities, namely the incommensurability parameter ζ and magnetization M, have not been measured as a function of the external field in any clean realization of Dzyaloshinskii s model to date. In the present paper we report the first direct experimental observation of a Dzyaloshinskii-type field-induced CI transition in Ba 2 CuGe 2 O 7, also presenting experimental data for ζ(h) and M(H). On the theoretical side we go beyond a qualitative analysis of the critical properties (close to the phase transition), as was previously done by Dzyaloshinskii. Dealing with this particular system, we construct an exactly solvable model and derive exact results that are in quantitative agreement with experiment throughout the entire phase diagram. The structural and magnetic properties of Ba 2 CuGe 2 O 7 are discussed in detail in Ref. [13]. In the layered tetragonal crystal structure the magnetic Cu 2+ sites form a square lattice with nearest-neighbor (nn) distances of 6 Å along the (1,1,0) and (1,-1,0) directions, respectively. The lowtemperature(t N = 3.26K)magneticphaseisaweakdistortion of the Néel spin arrangement, with all spins confined to the (1, 1,0) plane and the staggered magnetization slowly rotating upon translation along the (1, 1, 0) direction. The propagation vector is (1+ζ,ζ,0), where 1

2 ζ = Only the nearest-neighbor in-plane antiferromagnetic exchange constant is significant and is equal to J ab = 0.48 mev. The coupling between adjacent Cuplanes is ferromagnetic, with J c / J ab 1/37. We have previously suggestedthat the incommensurate structure is a result of Dzyaloshinskii-Moriya (DM) [14,15] antisymmetric exchange interactions. The corresponding term in the Hamiltonian can be written as D(S 1 S 2 ), where denotes the vector product and D is a vector associated with the oriented bond between the two interacting spins. From the symmetry properties of the lattice we deduce that for a bond between two nn Cu sites along the (1,1,0) direction, the only allowed components of D are those along (1, 1,0) and (0,0,1), respectively (Fig. 1). The (1, 1, 0)-component does not change sign from one bond to the next, while the (0,0,1) component is sign-alternating. It is the uniform component that is liable for the incommensurate distortion of the Néel structure, and for the rest of the paper we shall ignore the oscillating component, assuming that D is in the (a, b) plane. The interaction energy is minimized when all spins are perpendicular to D. The total exchange energy of the pair of nearest-neighbor spins is given by 2J ab cosφ + Dsinφ = 4Jab 2 +D2 cos(φ α), where φ is the angle between spins, and α = arctand/2j ab. The energy is a minimum at φ = π+α. The classical ground state is therefore a spin-spiral, with all spins in the (1,- 1,0) plane, and the angle between subsequent spins equal to π +α. The angle α is related to the propagation vector ζ by α = 2πζ 10. Precisely this spin configuration was found in the initial zero-field neutron diffraction experiments [13]. The central experimental result of this paper is the observation of a CI transition in Ba 2 CuGe 2 O 7, induced by a magnetic field applied along the (0,0,1) direction. Single-crystal magnetization measurements were performed using a conventional DC-squid magnetometer in the temperature range K. Experimental χ(h) dm dh for T = 2 K are shown in Fig. 2. For H a no anomalies are observed. In contrast, when the field is applied along the c axis, a distinct feature is seen around H = 2 T and indicates the presence of a magnetic phase transition. Neutron diffraction experiments were carried out on the H9 (cold beam) and H4M (thermal beam) 3-axis spectrometers at the High Flux Beam Reactor (HFBR) at Brookhaven National Laboratory on a mm 3 single-crystal sample, in the temperature range K and magnetic fields up to 6.5 T, applied along the c-axis of the crystal [16]. In zero field, at T = 2.4 K< T N, elastic scans along the (1 + ζ,ζ,0) direction show magnetic Bragg reflections centered at an incommensurate position ζ = ± (Ref. [13], Fig. 3b). As the magnetic field increases, the peak moves in closer to the antiferromagnetic zone-center at (1, 0, 0) [Fig. 3, insert (a)]. At H > H c 2.3 T the satellites at (1±ζ,±ζ,0)arenolongerobserved,but arereplacedbya singlepeakatthec-point(1,0,0)[fig.3, insert(b)]. The magnetic structure thus becomes commensurate. The experimental field-dependence of the propagation vector ζ is shown in the main panel of Fig. 3. To quantitatively describe the field-dependent behaviour, we follow the approach of Dzyaloshinskii[9]. We assume that the vector of local staggered magnetization at point r remains in the (1, 1,0) plane and forms angle θ(r) with respect to the c-axis. The free energy per Cu-plane in the continuous limit is then given by: ( ( θ(r) F = [ ρ s 2 x α ) 2 ( ) 2 ( ) ) 2 θ(r) θ(r) + +γ Λ y z (χ χ )H 2 2 sin 2 θ(r)]dxdy. (1) Here the axes x, y and z run along the (1,1,0), (1, 1,0) and (0, 0, 1) directions, respectively, and r = (x, y, z). The first term is the total elastic energy of exchange interactions, and favorsa spiral structure of period 2πΛ α. In Eq. (1)Λisthe in-planenncu-cudistance, andρ s isthe in-plane spin stiffness, that for classical spins at zero temperature is given by ρ s (0) = Jab 2 +D2 S 2. γ is the spin stiffness anisotropy defined by γ(λ c /Λ) 2 = J c / J ab 1/37 for Ba 2 CuGe 2 O 7. The second term represents the Zeeman energy. χ (T) and χ (T) are defined as magnetic susceptibilities with respect to fields that rotate along with the spiral structure and are parallel or perpendicular to the local staggered magnetization, respectively. In the paramegnetic phase χ (T) = χ (T), while at T = 0 the classicalresult is χ (0) = (gµ B ) 2 /(8J ab Λ 2 ), and χ (0) = 0. The equilibrium spin configuration should minimize the free energy (1), and therefore satisfy 2 θ x 2 = (χ χ )H 2 sinθcosθ = 1 ρ s 2Γ 2 sin2θ, (2) where Γ = (ρ s /H 2 (χ χ )) 1/2. Expression (2) has the form of the sine-gordon equation, that is central to describing CI transitions in many systems, and its soliton lattice solutions are wellknown: θ(x) = am(x/βγ,β), (3) where am(x, β) is the Jacobi elliptic function of modulus β. Analogues of Eqs. (1) and (2) were derived in Ref. [9]. To obtain exact results for ζ(h) and M(H), however, we make one additional crucial step. For each value of H of all valid solutions, labeled by β, one has to choose the one that indeed corresponds to the global minimum of the free energy. This is done by substituting Eq. (3) into Eq. (1) and minimizing the resulting expression with respect to β. After some algebra we find that F is minimized when 2

3 β E(β) = H ; (4) H c H c = πα ρs, (5) 2Λ χ χ where E(β) is the elliptic integral of the second kind. H c is the critical field at which the CI transition occurs [9]. Indeed, for H > H c the spin structure is given by the soliton-free solution θ(r) π/2, which corresponds to a commensurate spin-flop phase, as visualized in Fig. 4(a). The staggered magnetization in this case is parallel to the x-axis and the spins are slightly tilted in the direction of the field. In the limit H = 0 one has β 0 x Λ and θ(x) = α 2π, which corresponds to an unperturbed sinusoidal spin-spiral [Fig. 4(c)]. Most interesting is the case 0 < H < H c, where the spin structure may be described as a soliton lattice: regions of the spin-flop phase are interrupted at regular intervals by magnetic domain walls, or solitons [Fig. 4(b)]. In each soliton the direction of staggered magnetization rotates by an angle π. At H H c, the density of solitons starts to decrease very rapidly, as 1/ ln(h c H) [9]. The transition is thus almost first order: as H H c, the two magnetic satellites at (1±ζ,±ζ,0) converge to eventually produce a single peak at (1,0,0). The exact expression for β(h) [Eq. (4)] enables us to derive parametric equations for the field-dependence of magnetization and incommensurability parameter, and directly compare these predictions to experimental resultsforba 2 CuGe 2 O 7. UsingtheformulaM = F/ H and the equalities for derivatives of elliptic functions [17] one gets: M = χ H +(χ χ )H 1 β 2 ( 1 E(β) ) K(β). (6) The magnetization curve is continuous at the critical field, while χ(h) dm dh diverges as H c is approached from below, and is constant and equal to χ at H > H c. In zero field χ is equal to (χ +χ )/2. We now use Eq. (6) together with the formula (4) to fit the experimental χ(h) for Ba 2 CuGe 2 O 7 measured at T = 2 K. With χ = emu/g, χ = emu/g and H c = 1.88 T a very good fit is obtained (Fig. 2, solid line for H c). In the refinement we have included a linear contribution to χ(h). This term is present in both χ c (H) and χ a (H) and is a resultof non-linearcorrections to the local susceptibility, neglected in our treatment. As H H c the soliton density decreases, i.e., the period of the magnetic structure increases with increasing field. In our case expressing ζ in terms of β is rather straightforward and yields: ζ(h) ζ(0) = H H c π 2 4βK(β) = π 2 4E(β)K(β). (7) Again, Eq. (7), when combined with Eq. (4), gives a parametric curve for ζ(h), that can be fit to the experimental data for Ba 2 CuGe 2 O 7 using H c as the only adjustable parameter [ζ(0) = is measured independently]. The result with H c = 2.13 T is shown in a solid line in Fig. 3. A remarkable agreement is obtained except very close to H c. Discrepancies close to the transition point are to be expected, since the transition is almost first order, the soliton lattice is very soft and easily pinned by any impurities. The same effects prevent us from experimentally observing a true divergence in χ(h) at the critical field (Fig. 2). Finally, we can estimate the actual value of the critical field using the previously measured exchange constants and ζ(0). With J ab = 0.48 mev [13] the exchange energy per bond is J = 2J ab S mev. For χ and χ we can use the values for a classical Heisenberg antiferromagnet at T = 0. ESR measurements [18] provide us with the gyromagnetic ratios of Cu 2+ : g a = and g c = Using Eq. (5) and χ = (gµ B ) 2 /(8J) we obtain H c = 3.3 T, that should be compared to the experimental value H c 2.1 T. Considering that in these estimates we have completely ignored quantum- and temperature corrections to χ and ρ s, a 30% consistency is indeed acceptable. In summary, we have observed a rare type of CI transition that is driven exclusively by the changing strength of the commensurate potential. The latter is directly controlled in an experiment by varying the magnetic field. A transition of this kind was envisioned over three decades ago by Dzyaloshinskii, and now we find that Ba 2 CuGe 2 O 7 exhibits it in its original form. Now that the underlying physics is rather well understood, Ba 2 CuGe 2 O 7 can be used as a very neat and simple model system for further studies of spiral magnetism. For example, it would be veryinterestingto lookcloserat the magnetic critical behavior. Another promising direction for future work is the study of the effect of an in-plane magnetic field. We wish to thank P. Bak, V. Emery and I. Zaliznyak for useful discussions. This study was supported in part by NEDO (New Energy and Industrial Technology Development Organization) International Joint Research Grant and the U.S. -Japan Cooperative Program on Neutron Scattering. Work at Brookhaven National Laboratory was supported by the U.S. Department of Energy Division of Material Science, under contract DE-AC02-76CH [1] Y. Frenkel and T. Kontorova, Zh. Eksp. Teor. Fiz 8, 1340 (1938). [2] F. Frank and J. van der Merwe, Proc. R. Soc. London A198, 205 (1949). 3

4 [3] R. Clarke, in Ordering in two dimensions, edited by S. Sinha (North Holland, Amsterdam, 1980), p. 53. [4] H. Zabel, ibid., p. 61. [5] J. A. Wilson, F. J. Disalvo, and S. Mahajan, Adv. Phys. 24, 117 (1975). [6] R. M. Fleming, D. E. Moncton, and D. McWhan, Phys. Rev. B 18, 5560 (1978). [7] W. C. Koehler, in Magnetic Properties of Rare Earth Metals, edited by R. J. Eliott (Plenum, New York, 1972), p. 81. [8] P. Bak, Rep. Prog. Phys. 45, 587 (1982); V. L. Pokrovsky, A. L. Talapov, and P. Bak, in Solitons, edited by S. E. Trullinger, V. E. Zakharov, and V. L. Pokrovsky (Elsevier Science, Amsterdam, 1986), p. 71. [9] I. E. Dzyaloshinskii, Sov. Phys.- JETP 20, 665 (1965). [10] Y.Ishikawa, K.Tajima, D.Bloch, andm.roth, SolidState Comm. 19, 525 (1976). [11] B. Lebech, J. Bernhard, and T. Flertoft, J. Phys: Condens. Matter 1, 6105 (1989). [12] M. Zhitomirsky, O.Petrenko, and L. Prozorova, Phys. Rev. B 52, 3511, (1995). [13] A. Zheludev et al., Phys. Rev. B 54, (1996). [14] I. Dzyaloshinskii, Soviet Physics-JETP 5, 1259 (1957). [15] T. Moriya, Phys. Rev. 120, 91 (1960). [16] A complete experimental report will be published elsewhere. [17] S. Gradstein and I. Ryzhik, Tables of Integrals, Series, and Products (Academic Press, New York, 1980). [18] Y. Sasago, unpublished data. FIG. 1. A schematic view of a Cu-Ge-O layer in Ba 2CuGe 2O 7. The arrows indicate the components of Dzyaloshinskii vector D, allowed by the symmetry. FIG. 2. Field dependence of the magnetic susceptibility measured in Ba 2CuGe 2O 7 at T = 2 K for the magnetic field applied along the c (circles) and a (crosses) axes of the crystal, respectively. The solid line is a theoretical fit to the data, as described in the text. FIG. 3. Field dependence of the magnetic propagation vector in Ba 2CuGe 2O 7 measured at T = 2.4 K. The solid line is a theoretical fit given by Eqs. (4,7). Insert: Elastic scans across the antiferromagnetic zone center measured in Ba 2CuGe 2O 7 at T = 2.4K for two different values of magnetic field applied along the (0,0,1) direction. FIG. 4. Spin configurations for the spin-flop phase at H > H c (a), the soliton lattice at H = 0.997H c (b) and the circular spin spiral at H = 0 (c). 4

5 Cu a D z GeO 4 D y D z c D y b Zheludev et al. Fig. 1

6 60 Ba 2 CuGe 2 O 7 T=2K H c χ (10-6 emu/g) H a H (T) Zheludev et al. Fig. 2

7 ζ (r. l. u.) I (counts/40 s) 400 (a) Q=(1+ζ,ζ,0) H=1.9 T (b) Q=(1+ζ,ζ,0) H=2.4 T Ba 2 CuGe 2 O 7 T=2.4 K H c =2.13 T ζ (r.l.u) H (T) Zheludev et al Fig. 3

8 (a) (b) H/H c =1.10 H/H c = (c) H=0 Zheludev et al. Fig. 4

Luigi Paolasini

Luigi Paolasini paolasini@esrf.fr LECTURE 5: MAGNETIC STRUCTURES - Mean field theory and magnetic order - Classification of magnetic structures - Collinear and non-collinear magnetic structures. - Magnetic