Symmetry Analysis of Multiferroics A. B. Harris (NIST, 2007) The objective of this set of lectures is to describe in detail how to apply symmetry considerations to the magnetic and ferroelectric states of multiferroics which have a combined magnetic and ferroelectric phase transition. The aim is to do this is as transparent a way as possible. 1
2 OUTLINE I. II. III. IV. V. Setting the Stage Symmetry of the Quadratic Free Energy A. Physics of this Free Energy B. Eigenvector Condition (Tb. Mn. O 3) C. Group of the Wavevector (Tb. Mn 2 O 5 , YMn 2 O 5) D. Impose Inversion Symmetry E. Introduce and Analyze Order Parameters Consider Quartic Terms in the Free Energy Construct the Magnetoelectric Coupling using the Symmetry of the Order Parameters Discussion
HISTORICAL REMARKS Prior to our 2005 Physical Review Letter 1 it was not recognized that, in addition to the symmetry operations which leave the waveector invariant, one could exploit inversion symmetry. (This may have been recognized in some quarters, but certainly none of the multiferroic community knew it!) This is documented in Ref. 2, on which these notes are based and to which citation may be made. A reader who understands representation theory can skip to slide #27!! 3
INVERSE SUSCEPTIBILITY The susceptibility is a key quantity in continuous transitions like those we are about to study. One is used to saying that the susceptibility diverges at a continuous transition. Alternatively, one can plot the free energy, F, (see the appendix) of a ferromagnet versus the magnetization, M. By time reversal only even powers of M can appear. The actual value of M is the one which minimizes F. Free energy F = (T-Tc)M 2 + u. M 4 for a sequence of temperatures T< < Tc < T 1 < T 2. Note that the free energy is unstable relative to the formation of long-range order for T < Tc 4
INVERSE SUSCEPTIBILITY For small H the magnetization M is 5 Ferromagnet T 2 > T 1 > Tc The coefficient of M 2 is the inverse susceptibility: In the Appendix we allow for all Fourier components: When the first Fourier component condenses, the others are inhibited from condensing because there is a finite amount of spin to go around. For a ferromagnet the instability first occurs at q=0 as one reduces the temperature.
INVERSE SUSCEPTIBILITY For the antiferromagnet the instability leads to up-down-up-down spins, i. e. aq = p. Observe that if you have a plot of the susceptibility versus wavevector, you don’t have to know whether the system is ferromagnetic or antiferromagnetic: the location of the minimum in inverse chi tells you which it is. Antiferromagnet The instability occurs at wavevector q=p/a 6
WAVEVECTOR SELECTION: INCOMMENSURATE SPINS For competing interactions, the minimum in inverse chi can be anywhere in the zone, as at right. The minimum in inverse chi locates the wavevector of the ordered state. This is often referred to as “wavevector selection. ” So far we have treated M as a scalar, as it is for the Ising model. What happens if M is a vector? 7 Incommensurate Magnet T = Tc T = T 1 > Tc T = T 2 > T 1
8 INVERSE SUSCEPTIBILITY Suppose we have vector spins with an easy axis along z. Then Where K>0 is the anisotropy energy. Here as T approaches Tc from above, only Mz becomes unstable relative to long-range order. As another example suppose the free energy is (Continued next slide)
INVERSE SUSCEPTIBILITY How do we deal with Now we have an anisotropy tensor, which we can diagonalize to find the anisotropies associated with each principal axis. Keeping only the order parameter associated with the easiest axis, , we have This example teaches us that if the free energy is a quadratic form, ordering takes place in the channel defined by the eigenvector associated with the eigenvalue that is vanishing. 9
INCOMMENSURATE CHI Now we consider the free energy of a system which has n>1 vector spins in the unit cell: Let Sa(R, t) denote the a-component of the spin at position t within the unit cell at R. Then where sa(q, t) is complex with sa(-q, t)=sa(q, t)*. The quadratic free energy is the Hermitian form This quadratic form has complex-valued eigenvectors which we can normalize, but they have an arbitrary phase factor (which I will discuss later). Obviously we focus on the critical eigenvalue, i. e. the one that first becomes unstable. 10
INTERPRETATION For each value of the wavevector we have 3 n eigenectors and 3 n eigenvalues. As the temperature is lowered we have wavevector selection: at some q, one of these eigenvalues will become unstable. This determines the wavevector of the ordered state (at least for temperature just below Tc). Now consider the eigenvector. Its normalization is not fixed by the quadratic free energy. The quartic terms fix the normalization and give the square root temperature dependence characteristic of mean field theory. Reminder: all the coefficients of the critical eigenvector (which represents the Fourier transform of the spin distribution) are complex-valued. Note: we reject the possibility of accidental degeneracy. Degeneracy can only occur if symmetry requires it. 11
12 INTERPRETATION So we write the eigenvector Y as where the components inside the square bracket are normalized so that the sum of the squares of their magnitudes is unity. I will discuss the phase f later. The amplitude X varies as (Tc-T)1/2. sa, t is the complex amplitude of the wave of the a-component of the t sublattice. So the eigenvector tells us the pattern into which the spins order, as we will see in a moment.
APPLICATION 13 One might attempt to write down and anlyze the matrix of coefficients F. That is a hopeless task because our microscopic models are not reliable. Instead, we use the symmetry of the crystal to limit the possible eigenvectors which are determined by fitting the diffraction data. The standard approach is to invoke group theory. I will avoid that as much as possible. Instead I will invoke the principle that first-year graduate students learn: if a set G 1, G 2, …of mutually commuting operators commute with the Hamiltonian, then the eigenvectors of the Hamiltonian can be classified according to their eigenvalues g 1, g 2, … In the language of group theory, this case is the case when all the irreducible representations are one dimensional.
14 CALCULATION The essential input we require is the symmetry of the incommensurate wavevector (usually as determined by neutron diffraction) i. e. the symmetry operations of the crystal which also leave this wavevector invariant. Since we are concerned with phase transitions in which ferroelectricity appears in coincidence with magnetic ordering, spatial inversion I is a significant symmetry operation, although it does not usually leave the wavevector invariant. So I’m afraid we have to actually confront the crystal structure of the paramagnetic phase.
15 EXAMPLE A: Tb. Mn. O 3 The space group operations 3 are defined in Table I and the ion positions are listed in Table II (x=0. 9836, y=0. 0810)4. The figure below shows the mirror plane (z=1/4) for mz and the glide plane (x=3/4) for mx with associated translation uy = b/2. From neutron diffraction 5, 6 the wavevector lies along b (i. e. y) and in rlu is q= 0. 28. The wavevector is invariant under mx and mz. Table I: General positions for Pbnm (TMO). 3, 4 Here x means –x. Table II: Mn and Tb positions. 4 = Mn = Tb
Eigenfunctions for Tb. Mn. O 3 (TMO) 16 In Table III I give the eigenfunctions for TMO. First I explain how to read this table. In each column headed by the symmetry label (Gn) I list the components of eigenfunctions on the sites 1 – 8 as labeled in Table II. The three numbers in each box are for the x, y, and z components. These eigenfunctions are eigenfunctions of mx and mz with respective eigenvalues l(mx) and l(mz). Here L= exp(ipq), where q is the incommensurate wavevector in reciprocal lattice units (rlu), so Table III: Eigenin ``real” units, it is 2 pq/b. vectors. L = exp(ipq).
EIGENFUNCTIONS 17 The eigenfunction for G 1 has five parameters x. M, y. M, z. T 1, and z. T 2 indicating that there are five independent eigenfunctions of this symmetry. So far, there is no restriction on these five complex parameters. As we will see below, inversion symmetry places important constraints on these parameters. On the next slides we will show the eigenvectors are easily constructed. Table III (repeated)
EIGENFUNCTIONS To construct eigenfunctions of mx and mz it is convenient to write Y in terms of subvectors as where each symbol is a vector with four components. X(M) has the x-components of the four Mn sublattices, Y(M) and Z(M) the y and z components, respectively, of the four Mn sublattices. Similarly, X(T), Y(T) , and Z(M) have the x, y, and z components, respectively, of the four Tb sublattices. The representation of Eq. (2) is convenient because the symmetry operations only interrelate the four components within a given symbol. Consequently, each eigenfunction has components only with in a single four-component subvector. 18
EIGENFUNCTIONS The main ideas involved in determining how spin functions transform are: 1) The transformed spin at r is obtained by OS acting on the spin which was at the original location [OR]-1 r, 2) exp(iqr) is a number and O does not act on it, and 3) the subscripts ``i” and ``f” mean ``before” and ``after” transformation. In Table IV we list the initial and final sites under the operations we will invoke. Table IV (at right): initial site, ti, sites, tf, after mx, mz, and I. The sites are numbered according to Table II. 19
EIGENFUNCTIONS To construct the eigenfunctions of mx and mz, we need to discuss how the symmetry operations act on Fourier transforms of spin functions like sx(q, 1). First of all, spin is a pseudovector, so the mirror ma does’t reverse the sign of the a-component of spin but does reverse the sign of the other components. Let O denote an operator acting on spin and space, OR (OS) an operator acting only on space (spin), and lower case denote a Fourier transform. On the next slide we indicate the effect of a symmetry operator O=OR OS on a Fourier transform. 20
EIGENFUNCTIONS OR = a point group operation GO + a translation u. O: so that (since GO leaves q invariant) Thus, the final result is 21
EIGENFUNCTIONS The result of the previous slide was For mz, u 0=0, and for mx, q. u 0=pq (recall rlu!). Also, OS gives a plus or minus sign. Thus where the sign depends on O and a and t (t’) is the initial (final) sublattice. Thus for the pseudovector spins, with L = exp(ipq), where xab = -1, except that xaa = +1, so that mz 2 = 1 and mx 2 = L 2. Note: ( t’)’=t. 22
EIGENFUNCTIONS Now we show the eigenvectors are constructed. For instance, to construct the eigenvector scaled by x. M, we assign m 1 x the value 1 and then deduce the other components according to the symmetry (l ‘s). So we introduce projection operators To check the effect of Px, use mx 2 = L 2 = l(mx)2: 23
EIGENFUNCTIONS Thus P Px. Pz projects an eigenfunction of both mx and mz with respective eigenvalues l(mx) and l(mz). We extend the notation of Eq. (2) so that (a, b, c, d)s. T represents a subvector S(T), whose four components are a, b, c, and d, and where s=x, y, or z (to label components) and T=M or T (for Mn or Tb). We now illustrate the use of these projection operators to construct eigenfunctions of mx and mz. To get the eigenfunction Y scaled by x. M , we let the projection operator act on X(M) = [1, 0, 0, 0]x. M in the notation of Eq. (2): 24
EIGENFUNCTIONS In the first line we used Table IV which tells us that mz takes site #1 into site #3 and Eq. (5) with xzx=-1. In the third line Table IV tells us that mx takes site #1 into site #2 and #3 into #4 and we used Eq. (4) with xxx=1. The above eigenvector regulates the x-components on the Mn sites. To get the amplitudes for G 1 set l(mz)=1 and l(mx)=L. This leads to Y=[1, 1, -1]x. M in agreement with the eigenvector of Table IV scaled by x. M. 25
26 EIGENFUNCTIONS To further illustrate the construction of the eigenvectors we apply the projection operator to the Tb wavefunction [1, 0, 0, 0]x. T : So for l(mz)=+1, there are no eigenvectors of the x-component. Otherwise the eigenfunction covers only sites #5 and #8. One can verify that the above result agrees with Table III. Enough checks!
EIGENFUNCTIONS: SUMMARY For convenience I repeat Table III at right. For those rejoining here, each column contains the most general eigenvector having the symmetry specified either by the eigenvalues of the symmetry operators, or by the irreducible representation G. At this stage of the development the 5 (or 7) free parameters may assume arbitrary complex values. The guiding principle is that we do not admit the possibility of accidental degeneracy. When ordering takes place only one symmetry eigen-vector can condense. The Fourier transform of the spin distribution is proportional to this ``critical” eigenvector. The symmetry (G) and the amplitudes (x, y, z) are determined to best fit to diffraction data. 27
28 INVERSION SYMMETRY Up to now our results are absolutely standard. It is often said (in reviews 7) that the only symmetries that can be exploited are those that leave the wavevector invariant. This is not so because when one has inversion symmetry, the free energy clearly must be invariant with respect to inversion. 1, 2, 8 We now analyze the consequences of that. We assume that only a single representation is critical and write the free energy as an Hermitian form in in the amplitudes (which we here denote Qi) of one of the columns (G) of Table III: So Q 1(G) = x. M, Q 2(G) = y. M, etc.
INVERSION Clearly, we now have to determine how inversion I transforms the coordinates Q. For that purpose we modify the derivation of Eq. (3) so that it applies to I. Replace ``leaves q invariant” by ``changes the sign of q. ” Thereby we get What does this mean? The initial and final sites are given in Table IV. For the Mn sites, they are at centers of symmetry so tf = ti. So Isa(q, n)=sa(q, n)* for Mn sites (n=1, 2, 3, 4). For the Tb sites the pairs of initial and final sites are 5 and 7 and 6 and 8, so Isa(q, 5)=sa(q, 7)* and Isa(q, 6)=sa(q, 8)*. 29
30 INVERSION Now we deduce the effect of inversion on the symmetry adapted coordinates, denoted Qn(G), which are the amplitudes in terms of the eigenfunctions of Table III. As an illustration we consider the coordinates of G 1 from column 1 of Table III:
INVERSION From the discussion of Eq. (6) one sees that Now write the free energy (for G = G 1) in terms of the Q’s and require invariance under inversion: Here and below we invoke the convention that Roman letters are real and Greek ones complex. Thus 31
INVERSION To repeat: We now need to study how inversion affects the matrix. For example, IQ 1* Q 4 = Q 5* Q 1. After inversion the 5, 1 element of the matrix is the old 1, 4 element. A simpler example: IQ 1* Q 2 = Q 2* Q 1, so after inversion the 2, 1 element is the old 1, 2 element. 32
INVERSION Since these two forms have to be equal no matter how we choose the Q’s, the matrices have to be equal. Thus a, b, and e must be real and Thus the free energy is of the form (with Roman letters real and Greek ones complex). The critical eigenvector is an eigenvector of the above matrix. 33
INVERSION To repeat: we found the free energy to be Because of the special form of this matrix (due to invoking inversion symmetry) the eigenvector is shown in the Appendix to be of the form 6 where the subscript Q says that the components are Q 1, Q 2, Q 3, … and we normalize by a 2+b 2+c 2+2|r|2=1. Since the eigenvalue problem takes place in a complex vector space, we can not omit the arbitrary phase f. We discuss this result below. 34
INVERSION We now discuss the result of the last slide: Previously (in Table III) we showed that symmetry led to an eigenvector of symmetry G 1 which had five complex coefficients. Now we see that apart from an overall phase factor, three of these are real, and the other two (for Tb), rather than being independent complex numbers are the complex conjugates of one another: So instead of having 5 complex amplitudes to fit the diffraction data, we have three real and one complex parameter. The result for G = G 3 is similar. The result for G 2 (G 4) is (see the Appendix) 35
INVERSION So at right we have updated Table III to take account of inversion symmetry. The eigenfunctions of each symmetry Gn should be assigned an independent phase factor fn and amplitude An. If only one symmetry is present, the phase is usually irrelevant. But if two symmetries are simultaneously present, the phases matter because the relative phase has significance, as we shall see. Reminder: Roman letters are real and Greek ones complex. 36 As Table III, but with Inversion symmetry.
37 ORDER PARAMETERS We now discuss how the magnetic structure can be described by order parameters. 2, 8 Suppose that the structure has Gn symmetry and the spin wavefunction is put into the form of Eq. (8). How do we think about an order parameter? Here the order parameter has to be complex because any phase is allowed in Eq. (8). As a function of temperature the normalized eigenvector will be essentially constant, and only the magnitude A will vary. So we identify the complex order parameter as where An is real.
38 ORDER PARAMETERS The order parameter inherits the symmetry of the representation Gn. Thus At right is the phase diagram 9 of TMO. As one cools at H=0 through. T> = 40 K, order in s 3 appears 6 in the HTI (high-temperature incommensurate) phase. 8 Then at T=28 K additional order in s 2 appears 6 in the low temperature (LTI) phase. Cooling at H>10 T leads to a different phase. The polarization is nonzero only in the LTI phases. TMO Para P=0
MAGNETIC STRUCTURE Here I discuss qualitatively the magnetic structure. Neutron diffraction 6 indicates that in the HTI phase the ordering with s 3 (i. e. G 3) involves the spins pointing along the b (i. e. y) axis. Since the a, b, and c directions are inequivalent, one axis must be easiest, and obviously this is the b axis. Note that the pattern of sublattice spins is (++--). The quartic terms in the free energy favor fixed spin length 8 which is not consistent with a collinear sinusoidal state. Thus, it is not surprising 8, 10 that a phase transition occurs at lower temperature where transverse order appears along the c axis scaled by the order parameter s 2, for G 2. The pattern of these transverse moments is again (++--) indicating nearly isotropic interactions. 39
PHENOMENOLOGY This behavior is consistent with the following simple phenomenological Landau expansion: and is consistent with the famous review of incommensurate systems by Nagamiya. 10 It is not obvious that the above free energy ensures that the wavevectors of the two order parameters are equal. If the interactions are nearly isotropic, one would expect this to be the case. Within Landau theory the term that locks to two order parameters is of the form 8 where u<0 and this term only exists if q 2 = q 3. 40
MAGNETOELECTRIC COUPLING Now we discuss the Landau theory of magnetoelectric (ME) coupling V. 1, 2, 8 Order B is induced by order A if there is a term in the free energy linear in B times 2 (or more) powers of A. Time reversal invariance ensures that such a term linear in the magnetization can not occur, so the polarization P can not induce magnetic order. But magnetization can induce P. The interaction we seek is therefore of the form V = ss P and the ME free energy is then where the electric susceptibility c. E is finite since we assume that without magnetism there is little tendency to have a spontaneous polarization. 41
42 ME COUPLING The ME coupling is schematically of the form Where c is a constant. Now minimize with respect to Pa: so this interaction does just what we want: it induces a polarization only when ss is nonzero.
43 ME COUPLING As we argued, the magnetoelectric interaction V is linear in P and quadratic in the magnetic order parameters. To conserve wavevector V is where G is a reciprocal lattice vector and delta is unity if q+q’=G and is zero otherwise. We drop the Umklapp term with G nonzero (for which see Ref. 11). Then we have
ME COUPLING Warning: the discussion that follows applies to TMO. For different symmetry systems the ME Interaction need not assume the form found below. 2, 12 In TMO, V must be invariant with respect to mx, mz, and I. To handle the cases that will be of interest to us, let us assume that there are two irreducible representations (irreps) present. Gr and Gs. In Eq. (10) the term with n=m is zero because it changes sign under I. Physically this reflects the fact that a single wave has a center of symmetry which, for an incommensurate q, can be as close to an inversion center of the lattice as you like. Thus 44
45 ME COUPLING To repeat: This has to be invariant under mz. Using Eq. (9) we have This equals. V only if mz. Pa=-Pa, i. e. for a=z (i. e. c). So One can show, using Eq. (9) that V is invariant under mx.
46 ME COUPLING We have V as This has to be invariant under inversion. We have To have IV = V, az must be imaginary: az=ir, with r real. We set sn(q) = |sn| exp (ifn). Then The polarization in TMO (at low H) is only found 8 in the LTI phase where it lies along c (i. e. z) as in (11).
ME COUPLING The appearance of the relative phase is easily understood. The phase of the irrep regulates the location of origin of the cosine wave. The wave has inversion symmetry about this point. So, if the two irreps have their origin at the same point, i. e. if f 2 = f 3, then the entire structure has inversion symmetry about a lattice site. The fact that P is proportional to |s 2 s 3| also agrees with experiment 8 in that the temperature dependence of Pz looks like that of an order parameter and not (order parameter)2. This happens because the HTI order parameter s 3 is more or less saturated by the time one gets into the LTI phase. 47
HIGH FIELD PHASE In the high field phase (see the figure on slide #38) the polarization lies alog a (i. e. x) in the analog of the LTI phase. Thus the magnetic structure must be s 3 + sn such that their product transforms like Pa. For the product to be odd under mx, we need n=1 or 4. For it to be even under mz, we need n=1 or 3. So this analysis indicates that n=3. (This was first predicted in Ref. 13. ) It should be noted that having a different irrep does not necessarily mean a totally unrelated structure. It may simply mean that the structure has been rotated. Because different spin components transform different under symmetry operations, this causes a change in the irrep. 48
49 Tb. Mn 2 O 5 (T 25) T 25’s space group is Pbam, 14, 15 whose symmetry operations are in Table VI and whose ion positions are in Table VII. Below T=43 K, T 25 has a commensurate phase with q=(1/2, 0, 1/4) and an incommensurate phase with q=(0. 48, 0, 0. 32). 16 Table VI. Pbam space group operations. 3 Since the irreps for T 25 are two dimensional, the simple method used above is clumsy. 2 An Table VII. Positions of alternate way, used here, is to supplement the standard analysis spins in T 25: x=0. 09, y=-0. 15, z=0. 25, 14 by including inversion symmetry. X=0. 14 and Y=0. 17. 15
50 T 25 The standard analysis (without inversion symmetry) for the wavevector (0. 50, 0, qz) leads to Table VIII, which is read just as Table III, except that the two columns do not represent different symmetry irreps, but rather transform like two columns of the same irrep. Since the two columns are degenerate in free energy, the actual spin function is a linear combination (determined by a fit to diffraction data) of the two columns. The complex coefficients s 1 and s 2 are the order parameters. Table VIII for T 25.
51 T 25 TRANSFORMATIONS Before discussing inversion symmetry let us consider the transformation properties of the eigenfunctions of Table VIII. To do that we need to know how the spin functions transform under the symmetry operations. Here it is convenient to redefine the spin Fourier transforms so that the exponential is exp (iq. R) rather than exp [iq. (R+t)]. With this change, a transformation 1) changes the orientation of the spin, 2) takes the spin from an initial sublattice to a final sublattice, and 3) introduces a phase factor (given In columns 4, 6, and 8 in the Table IX, below) depending on the position change. (See Ref. 2. )
T 25 TRANSFORMATIONS 52 In Table IX we give the way the sublattices transform. The spin transformation factor xab (for the symmetry operation ma) are the same as before: xab =1 if a=b and is -1 otherwise. Thus In the first equation, xxx=1, and in the second, xyz=-1. One needs the phase factors in Table IX, taken from Ref. 2. I does not carry a x factor. Table IX. As Table IV. For exp(if) and exp(if’) see the text. Inversion takes q into – q. xab is not included.
INVERSION SYMMETRY Now we discuss the effect of inversion symmetry. To do this we will study the free energy when it is written in terms of the symmetry adapted coordinates of Table VIII. If we take the 12 symmetry coordinates of column #1, then the free energy is a quadratic form with a 12 by 12 matrix Fnm. In terms of Fourier components the free energy is written as For convenience I will absorb the ½ into c-1. Now we introduce symmetry adapted coordinates for the first column of Table IX as we did in slide #30. These are listed on the next two slides. 53
SYMMETRY COORDINATES Below we give the Q’s which transform like the first column of the two-dimensional irrep. 54
SYMMETRY COORDINATES Below we give the Q’s which transform like the second column of the two-dimensional irrep. 55
SYMMETRY OPERATIONS Now we explicitly display the effect of the symmetry operations using Table IX and xab = 1 if a=b but is -1 otherwise: Using Table IX one can determine the effect of inversion on the Q’s and R’s. We find that where the necessary data is in Table X. Transformation parameters for inversion. 56
DOUBLE DEGENERACY The meaning of a two-dimensional irrep is that the free energy matrix F (shown on the next slide) for the R’s is exactly the same as that for the Q’s: This guarantees that all eigenvalues are doubly degenerate. What this means is that symmetry dictates that these two modes (or any linear combination of them) become unstable (relative to ordering) simultaneously. Thus the free energy in terms of their order parameters is 57
58 Here Sn is a 4 x 4 real symmetric matrix, Rn a 4 x 4 real matrix. The other matrices assume the forms shown below. F= Below we will verify this form. Meanwhile, see that F acting on Y, given below, yields the a vector of the same form. Roman = real, i 2=-1. Y= [a, b, c, d; z, z*; e, f, g, h; x, x*; im, in, io, ip; r, -r*] where a is real.
59 VERIFYING THE FORM OF F Now we use invariance under inversion I: (d) (a) (e) (b) (f) (c) (g) NOTES: When we have a result of the form Rn. XRm, then we can conclude (in the second line of these equations) that X=Hnm. In (a) H 12 is real. In (b) H 15=H 16*. In (d) H 1, 13 is imaginary. In (e) H 1, 17=-H 1, 18*. No info in (g).
T 25 SPIN STRUCTURE Note that the spin Fourier transform comes from a linear combination of the Q’s and the R’s. So we parametrize the eigenvector as (see bottom of slide #57) Y= s 1 [a, b, c, d; z, z*; e, f, g, h; x, x*; im, in, io, ip; r, -r*]Q Here the subscript Q indicates that these are the values of Q 1, Q 2, … Q 18 and we will normalize the vector so that the sum of the magnitudes squared of its components is unity. The matrix for the R’s is identical to that for the Q’s. This guarantees that the Q mode and the R mode become unstable simultaneously and that the critical R eigenvector has the same form as Y. 60
61 T 25 SPIN STRUCTURE Since the variables and matrices are degenerate, the eigenvector for the R’s is the same as that for the Q’s: F = s 2 [a, b, c, d; z, z*; e, f, g, h; x, x*; im, in, io, ip; r, -r*]R , where the subscript R indicates that these are the components of R. We want the spin Fourier components, which we get by writing When s 2=0, then all the R’s on slide #55 are zero, which gives relations between pairs of sa(q, t)’s. Likewise, when s 1=0, all the Q’s are zero and we have other relations between pairs of sa(q, t)’s.
T 25 SPIN STRUCTURE To illustrate: when s 2=0, then so that These results give the first column in Table X, below. The second column can be obtained similarly, or maybe more easily by applying my to the left column. In any event, the final result is shown in Table XI. 62
63 T 25 SPIN STRUCTURE In any event, the result for the spin structure of T 25 after taking account of inversion is given in Table XI, where the Roman letters are real and the Greek ones complex. Note The effect of inversion: the two Tb orbits (sites 5 -6 and 7 -8) are no longer independent, but now are complex conjugates of one another. Also the others are now restricted to be real. Of course, we have an overall undetermined phase which we put into the order parameters. Table XI. As Table III.
SYMMETRY OF THE ORDER PARAMETERS Up to now we have studied how the spin components transform under symmetry operations. We can summarize these results by instead attributing the transformation to the transformation of the order parameters. This is very convenient because this enables us to analyze the symmetry of Landau expansions in terms of these order parameters. Applying the transformation law of Table IX to the eigenvectors of Table XI, we see that (15) This information is crucial to constructing the ME interaction, as we shall see. 64
YMn 2 O 5 (YMO) The results for YMO are included in those for T 25 because the only difference is that in YMO there is no analog of the magnetic Tb sites. It is interesting to express the structure determined by Chapon et al. , 17 in terms of order parameters sn. The simplest way to identify the order parameters is to note the symmetry of the spin structure: mx takes it into the negative of itself and my takes it into an orthogonal state. Thus only s 2 is nonzero. 65
ME INTERACTION Here I discuss the ME interaction. Of course the Landau expansion will be the same for both T 25 and YMO. We start with the general trilinear form: Eq. (15) says that I[s 1 s 2*] = s 1 s 2*. So only terms with n=m in V can be inversion invariant. So Eq. (15) says that the factor in square brackets is Odd under my. Therefore only a=y contributes: 66
67 ME INTERACTION To repeat: the final result for the ME interaction is which gives the spontaneous polarization (along y, as observed 18, 19) as This relation is also obtained for Rb. Fe(Mo. O 4)2, 12 where its prediction for the temperature dependence was confirmed by comparing P(T) with the neutron diffraction measurement of |s|2.
68 CONCLUSION The invocation of inversion symmetry usually cuts down to half the number of fitting parameters needed for magnetic structure determinations in incommensurate systems. The characterization of the magnetic ordering by the introduction of order parameters has many advantages. The most important one is that the symmetry of the unit cell (important in these magnetically frustrated systems) is properly taken into account, which is not the case 20 with continuum theories. 21
APPENDIX A: FREE ENERGY To help understand the free energy, F, I now give a simple approximation for it. F is the minimum of where H is the Hamiltonian and r is the trial density matrix (normalized to have unit trace) and which is varied to minimize F. Take H to be -JSi. Sj for nearest neighbors on a cubic lattice and let each S assume the values +1 and -1 (Ising spins). Take the density matrix to be such that the value +1 occurs with probability (1+ s )/2 and -1 with probability (1 -s)/2. Then Tr [Sir] = s = M is the order parameter and Where z=6 is the coordination number of the cubic lattice of N sites and -TS is the second term in Eq. (A 1) for F. (Continued next slide. ) 69
70 We now evaluate –TS: So up to quartic order the free energy is For T -> Tc-, s ~ (Tc – T)1/2.
71 I now describe the result when we allow the ordering to be according to a wavevector q, so that instead of the probability p of a spin being up being uniform, now it is that of a spin-density wave given by Of course, this kind of state for incommensurate q is not favored by a nearest neighbor interaction which gives either a ferromagnet or an antiferromagnet, depending on the sign of J. Therefore we allow for competing interactions and we introduce (Continued next slide. )
72 For this system the free energy is This formula holds for a ferrr- or antiferro-magnet. It is useful to identify the wavevector-dependent susceptibility c(q) as For nearest neighbor interactions on a cubic lattice this gives For a ferromagnet J is negative and the inverse susceptibility is minimal at q=0, as expected.
73 FORM OF EIGENVECTOR The matrix under consideration for irreps G 3 and G 1 is G = One can show that when G operates on a vector of the form [a, b, c, d, x+iy, x-iy], where all the letters represent real quantities, it gives back a vector of this same form. That is enough to say that any eigenvector is of this form. (continued on the next slide)
74 The matrix we are analyzing is G = Now we form U-1 GU, where Since U-1 GU is a real symmetric matrix, its eigenvectors fn are real valued. Then the eigenvector of G, namely Ufn has the form asserted.
75 The analysis for the G 2 and G 4 eigenvectors is similar. The free energy is now a function of 7 coordiates: = H(G 2) By equating H and IH we get relation, such as a, b, and e are real, g=d*, g’ = d’*, etc. Thus, inversion symmetry leads to the matrix H assuming the form Its eigenvectors are of the form [a, b, c, r, r*, t, t*].
76 A FEW REFERENCES 1. 2. 3. G. Lawes, et al. , Phys. Rev. Lett. 95, 087205 (2005). A. B. Harris, cond-mat/0610241 (submitted to PRB). A. J. Wilson, International Tables for Crystallography (Kluwer Academic, Dordrecht, 1995) Vol. A. 4. E. E. Sauerbrei et al. , Acta Crystallogr. B 29, 2304 (1975). 5. R. Kajimoto et al. , Phys. Rev. B 70, 012401 (2004). 6. M. Kenzelmann et al. , Phys. Rev. Lett. 96, 087206 (2005). 7. E. F. Bertaut, J. de Physique, Colloque C 1, 32, 462 (1971). 8. M. Kenzelmann et al. , Phys. Rev. B 74, 014429 (2006). 9. T. Kimura et al. , Nature 426, 55 (2003). 10. T. Nagamiya, in Solid State Physics, ed. F. Seitz and D. Turnbull (Academic, New York, 1967), Vol. 20, p 346. Continued on the next slide
77 11. 12. 13. 14. 15. 16 17. 18. 19. 20. 21. J. J. Betouras et al. , Phys. Rev. Lett. 98, 257602 (2007). M. Kenzelmann et al. , Phys. Rev. Lett. 98, 267205 (2007). A. B. Harris and G. Lawes, “Ferroelectricity in Incommmensurate Magnets, ” in The Handbook of Magnetism and Advanced Magnetic Materials, ” Ed. H. Kronmuller and S. Parkin, (Wiley, 2008); cond-mat/0508617. G. Buisson, Phys. Stat. Sol. 16, 533 (1973). G. Buisson, Phys. Stat. Sol. 17, 191 (1973). L. Chapon et al. , Phys. Rev. Lett. 93, 177402 (2004). L. Chapon et al. , Phys. Rev. Lett. 96, 097601 (2006). A. Inomata and K. Kohn, J. Phys. Cond. Mat. 8, 2673 (1996). Data for T 25 and YMO. I. Kagomiya et al. , Ferroelectrics 286, 167 (2003). Data for T 25. M. Kenzelmann and A. B. Harris, cond-mat/0610471 accepted as a comment (re Ref. 21) by Phys. Rev. Lett. M. Mostovoy, Phys. Rev. Lett. 96, 067601 (2006).
Скачать презентацию Symmetry Analysis of Multiferroics A B Harris NIST
e32763a5c4ae85089e6d2288dee9d5ca.ppt