An adaptive algorithm for n-body field expansions
Weinberg, Martin D.
1998-05-28
Confined spin waves reveal an assembly of nanosize domains in
ferromagnetic La(1-x)CaxMnO3 (x=0.17,0.2) 
Creator
Date
2004-08-27
Description
Comment: 4 pages, 5 figures
We report a study of spin-waves in ferromagnetic La$_{1-x}$Ca$_{x}$MnO$_3$,
at concentrations x=0.17 and x=0.2 very close to the metallic transition
(x=0.225). Below T$_C$, in the quasi-metallic state (T=150K), nearly
q-independent energy levels are observed. They are characteristic of standing
spin waves confined into finite-size ferromagnetic domains, defined in {\bf a,
b) plane for x=0.17 and in all q-directions for x=0.2. They allow an estimation
of the domain size, a few lattice spacings, and of the magnetic coupling
constants inside the domains. These constants, anisotropic, are typical of an
orbital-ordered state, allowing to characterize the domains as "hole-poor". The
precursor state of the CMR metallic phase appears, therefore, as an assembly of
small orbital-ordered domains.
at concentrations x=0.17 and x=0.2 very close to the metallic transition
(x=0.225). Below T$_C$, in the quasi-metallic state (T=150K), nearly
q-independent energy levels are observed. They are characteristic of standing
spin waves confined into finite-size ferromagnetic domains, defined in {\bf a,
b) plane for x=0.17 and in all q-directions for x=0.2. They allow an estimation
of the domain size, a few lattice spacings, and of the magnetic coupling
constants inside the domains. These constants, anisotropic, are typical of an
orbital-ordered state, allowing to characterize the domains as "hole-poor". The
precursor state of the CMR metallic phase appears, therefore, as an assembly of
small orbital-ordered domains.
Type
Database
Show preview
Hide preview
ar
X
iv :c
on d-
m at
/0 40
86 11
v2 [
co nd
-m at.
str -el
] 11
Ja n 2
00 5
Confined spin-waves reveal an assembly of nanosize domains in ferromagnetic La1−xCaxMnO3 (x=0.17, 0.2)
M. Hennion1, F. Moussa1, P. Lehouelleur1, F. Wang1,2, A. Ivanov3, Y. M. Mukovskii4 and D. Shulyatev4 1 Laboratoire Le´on Brillouin, CEA-CNRS, CE-Saclay, 91191 Gif sur Yvette, France.
2 Institute of Physics Chinese Academy of Sciences, 100080 Beijing, China 3 Institut Laue-Langevin, 156 X, 38042 Grenoble Cedex 9, France 4 Moscow State Steel and Alloys Institute, Moscow 119991 Russia
We report a study of spin-waves in ferromagnetic La1−xCaxMnO3, at concentrations x=0.17 and x=0.2 very close to the metallic transition (x=0.225). Below TC , in the quasi-metallic state (T=150K), nearly q-independent energy levels are observed. They are characteristic of standing spin waves confined into finite-size ferromagnetic domains, defined only in (a, b) plane for x=0.17, and in all directions for x=0.2. They allow an estimation of the domain’s size, a few lattice spacings, and of the magnetic coupling constants inside the domains. These constants, anisotropic, are typical of an orbital-ordered state, allowing to characterize the domains as ”hole-poor”. The precursor state of the CMR metallic phase appears, therefore, as an assembly of small orbital-ordered domains.
PACS numbers: PACS numbers: 74.25.Ha 74.72.Bk, 25.40.Fq
In metallic La1−xCaxMnO3, for x≈0.3, the origin of the colossal magneto-resistance (CMR) which occurs at the metal-insulator transition is still a challenge for the- ories. The role played by an inhomogeneous ground state has been outlined by several theoretical[1] and ex- perimental works[2, 3]. In the semi-conducting parent LaMnO3, the magnetic coupling is of super-exchange (SE) type, ferromagnetic (F) in the (a, b) plane (Ja,b> 0) and antiferromagnetic (AF) along c (Jc<0), consistent with the orbital ordering which occurs at the Jahn-Teller transition TJT . Doping with holes induces a canted AF state (CAF). There, neutron diffuse scattering in- dicates the existence of magnetic inhomogeneities at- tributed to hole-rich clusters embedded in a hole-poor matrix[4], and a new spin-wave branch appears at small- q around F Bragg peaks. This new branch coexists with the spin-wave branch close to that of LaMnO3 (SE type of coupling), attributed to the hole-poor matrix[5], which disappears at x=0.125 (F state). The ferromagnetic state which occurs in the range 0.125≤x<0.225 shows a quasi-metallic behavior below TC and an insulating one at lower temperature[6]. There, several works have also proposed an inhomogeneous picture using NMR[7], Mossbauer[8] or magnetization measurements[9].
In the present paper, we report a determination of spin-waves in ferromagnets La1−xCaxMnO3 at x=0.17 and 0.2 very close to the ”true” metallic state (x≥0.225). Below TC , in the quasi-metallic state, q-independent en- ergy levels of magnetic origin are observed in the large q-range. The q-dependence of their intensity, their evolu- tion with temperature and magnetic applied field suggest that they can be ascribed to standing spin waves charac- teristic of a confinement into small domains. The shape and the size (a few lattice spacings) of the domains can be estimated as well as the magnetic coupling inside the domains. All the observations are consistent with a pic- ture of hole-poor and orbital-ordered domains.
Inelastic neutron experiments have been carried out at the reactor Orphe´e (Laboratoire Le´on Brillouin) and at the Institut Laue-Langevin, on triple axis spectrom- eters, using fixed final wave vector of neutrons varying from 1.05 to 4.1A˚−1, and appropriate filters. The TC values are 175K and 180K for x=0.17 and x=0.2 respec- tively. Both samples have been characterized by resis- tivity measurements. In La0.83Ca0.17MnO3, the high- temperature pseudo-cubic structure, with lattice param- eter a0=3.89A˚, becomes orthorhombic at TJT=240K, whereas such an orthorhombicity is unobservable for La0.8Ca0.2MnO3. All Q values (Q=q+τ and q=ζ) are given in reduced lattice units with pseudo-cubic indexa- tion. The two single crystals are twinned, so that each di- rection is superimposed on related symmetry directions. Data are fitted by a lorentzian shape weighted by a Bose factor. This analysis yields the energy E(q), the damping Γ(q) with Γ(q)/E(q)≈1/5 and the intensity.
We first consider La0.83Ca0.17MnO3. The spin wave modes obtained at T=150K, 100K and 10K along [001]+[010]+[001] directions are reported in Fig. 1, left panel. At 150K and 100K, one may distinguish two dis- persive modes labelled (1) and (3). They are very close to previous measurements at x=0.125 reported in or- thorhombic indexation[5]. This comparison allows us to assign the curve (1), dispersing up to ≈18 meV, to [100] or [010] directions, which have been shown to be equiv- alent, and the curve (3), with a down-turn at q=0.25, to the [001] or c direction. We mention that this down- turn, also observed at x=0.125 and in the CAF state, is surprising for a F coupling. It denotes an underlying AF coupling reminiscent of the orbital ordering of LaMnO3. The very new feature is the gap-opening in the dispersive curve (1). Within this gap lies a nearly q-independent energy level, (2), at a mean value of 4.5 meV at 150K. With decreasing temperature, all the energies increase. At 100K and below, a single dispersion curve is observed
2 0.00 0.25 0.500.00 0.25 0.50 0
5
10
15
20
En er
gy (m
eV )
[100] [001]
0.0 0.1 0.2 0.3 0.4 0.5 0
5
10
15
20
½
½
0.0 0.1 0.2 0.3 0.4 0.5 0
5
10
15
20
E ne
rg y
(m eV
)
[001] B
x=0.2
0.00
(1)
(0)
E
T=150K
,
[001]
E [100] B
+related symmetry[1+ζ,0,0]+related symmetry[1+ ζ,0,0]
(0)
(2)
(1)
(3)
(1)
(2)
T=10K
(2)
(3) (2) E B E B
[100]
0.0 0.1 0.2 0.3 0.4 0.5 0
5
10
15
20 T=17KT=50K T=50K
x=0.17 T=100KT=150K
(3) (2)
(1) (1)
(3)
FIG. 1: Magnetic excitations measured along [100]+[010]+[001] for x=0.17 (left panel) and x=0.2 (right panel) at several temperatures. Full (empty) circles correspond to modes with main (weak) intensity. In left panel, the solid and broken lines are guides to the eyes, and the vertical lines locate 1/8, 1/4 and 3/8 q values. In right-upper panel, the dashed line indicates the spin-wave dispersion of a virtual large-size domain (see the text).
for q<0.125, larger gaps open at q=1/8, 1/4 and 3/8 and the level labelled (2) appears modulated with a maxi- mum close to 3/8. Measurements of spin waves along [110] and [111] at 150K also reveal q-independent energy levels lying at nearly the same values in the two direc- tions.
These observations call for the following comments. A q-independent energy level indicates localized excita- tions. The existence of a gap-opening in the branch (1) points out its interaction with level (2). This implies that level (2) belongs to the same macroscopic domain, namely to the (a, b) plane, but not to the c direction. The localized excitations or standing spin waves are con- fined within these planes, exhibiting a two-dimensionnal (2D) character. This is consistent with observations along [110] and [111], as indicated below. Since the dis- persionless level appears above qcr≈0.2, an approximate size ξ=a0/qcr≈ 20A˚ may be provided for the domain. A more quantitative determination of the size can be done by considering the energy value of the q-independent level. Following the arguments developped below for the
ζ=0.5 T=150K
ζ=0.275 x=0.2, [100]+[010]+[001]
(3) (2)
(1)
ph
ph
ph
(1) (2)
(3)
0 4 8 12 16 20 24 28 Energy (meV)
0
100
200
300
N eu
tr on
c ou
nt s/
12 0s
0 4 8 12 16 20 24 28 Energy (meV)
0
50
100
150
200
250
300
C ou
nt s/1
20 s
150K T=210K
FIG. 2: Examples of energy spectra along [100]+[010]+[001] at H=0, after application of a magnetic field H=2T. Left panel: ζ=0.275 at T=150K (circles) and 210K (triangles, shifted by 100 neutron counts). Right panel: ζ=0.5. The modes (1), (2) and (3) correspond to spin-wave branches in- dicated in Fig. 1, and ”ph” is a phonon mode.
x=0.2 compound, the ratio, ≈1/4, between the energies of the level (2) and of the zone boundary, determines a size of 4 lattice spacings (16A˚) along [100] (or [010]). In the insulating state and specially below 100K, this pic- ture is no longer valid. We may explain the origin of gaps and the modulation of the level (2) by assuming an un- derlying periodicity of 4a0 in the (a, b) plane. At the new zone boundaries of the super-cell, the dispersion is folded and new gaps occur. This description of the spin dynam- ics is likely related to observations in La7/8Sr1/8MnO3, where static superstructures are observed[11].
The study of La0.8Ca0.2MnO3 reveals a very interest- ing evolution with x. This sample has been first stud- ied in zero field. Preliminary observations at 17K along [001]+[010]+[001] have been reported[5]. Results are de- scribed with decreasing temperature. At 150K, unlike the x=0.17 case, a single dispersion branch, labelled (0), is observed in the small q range up to q≈0.3, indicat- ing that the [100] (a and b) and [001] ( c) directions are now equivalent and the dispersion isotropic (see Fig 1- right upper panel). For q≥ 0.3, the magnetic intensity is shared between several nearly q-independent levels de- fined by their mean values: a broad energy level, (1), at 18±1 meV, likely coupled to the lower optical phonon mode, a level (2) at 9±0.5 meV and an overdamped level, quasielastic. A magnetic field (H=2T) applied at 150K, suppresses the quasi-elastic intensity attributed to spin fluctuations, and reduces the energy linewidth of the dis- persionless magnetic levels. When removing the field at 150K, the spin fluctuations remain considerably reduced, and the levels keep the characteristics observed in applied field. A low-energy mode (3) is now defined at a mean value 4.5 meV (see Fig 1) and the mode (1) appears sep- arated from the phonon mode (24 meV), not reported
3 in Fig 1. Corresponding energy spectra are displayed in Fig. 2, showing the q-dependence of the intensity. As q increases from q=0.275, to 0.5, the intensity of the modes (2) and (3) is transferred towards the energy mode (1) (18 meV), so that the intensity follows the q-dependence of a F dispersion curve. Very interestingly, these modes persist above TC (180K). This is shown by the spectrum at 210K in Fig. 2, left panel, where a slight decrease of the two energy levels (2) and (3) is observed whereas a quasielastic intensity (spin fluctuations) increases above TC . In contrast, the small-q dispersive branch (0) renor- malizes to zero at TC . By lowering temperature in zero field, the lower-energy mode (3) nearly disappears and level (2) becomes more dispersed, which reduces the en- ergy gaps (Fig. 1, lower-right panel).
Spin wave measurements along [110]+[101]+[011] and [111] directions are reported in Fig 3 at T=150K. At small q, the quadratic dispersion ω=Dq2, (0), defines an isotropic stiffness constant (D=30±1meVA˚2 becom- ing 48meVA˚2 at 10K). At larger q values, the magnetic intensity is transferred to nearly q-independent energy levels. Along [111], a level (2) is determined at 23±1 meV and a slightly dispersive one (1) around 46±2 meV. Fig 4 displays the balance in intensity between the modes (0) and (2) as q varies. Along [110]+[101]+[011], the en- tanglement between magnon and phonon branches (see dashed lines in Fig. 3-left panel), led us to use polar- ized neutrons (PN) which allow to identify magnons and phonons. Several levels have been determined, more or less dispersive. Dispersion could have been induced by the applied magnetic field (H=2T) which aligns macro- scopic domains. Excitations at energy larger than 35 meV could not been detected because of a too weak inten- sity. By lowering temperature, the energy levels slightly increase.
These overall observations can be explained as it fol- lows. Unlike La0.83Ca0.17MnO3 where the confined spin waves have been assigned to the (a, b) plane only (2D character), for x=0.2, the isotropy of the small-q disper- sion suggests that all directions are concerned. The levels are attributed to standing spin-waves inside 3D clusters. Let us use the energy of the propagating wave for a 3D Heisenberg model (see dotted line in Fig. 1 upper-right panel):
E(qx,qy,qz)=8S{Ja,b[1-0.5(cos(2piqx)+cos(2piqy))]+ Jc/2[1-cos(2piqz)]}, where qx, qy, qz are related to a, b, c directions. This expression yields relations between Ja,b, Jc and the zone boundary energies
EB: E [100] B =E
[010] B =8SJa,b, E
[001] B =8SJc, E
[110] B =2E
[100] B ,
E [101] B =E
[011] B =E
[100] B +E
[001] B and E
[111] B =2E
[100] B +E
[001] B
These latter relations are also valid for finite-size clus- ters. According to predictions[10], the cluster’s size ξ is related to the level with the lowest-energy value EL. In the case of small clusters, it corresponds to a ”half-wave” or a half-period of wave. For the confined wave of energy
0.0 0.1 0.2 0.3 0.4 0.5 0
10
20
30
40
50
0.0 0.1 0.2 0.3 0.4 0.5 0
10
20
30
40
50 mg mg (PN)
�����
�����
�����
�����
�����
�����
�������������
�������������
�������������
�������������
�������������
�������������
�������������
�������������
�������������
�������������
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
2E +EB
(3)
(2)
B
(1)
+1/2E
Ca 20%
B [100]
EB
[100]
[001]
Ca 20%
(0)
[1+ζ,ζ,0] +related symmetry [1+ +related symmetryζ,ζ,ζ]
������
������
������
������
������
������
En er
gy (m
eV )
T=150K
(2)
(1)
(0)
T=150K [001]
(4)
FIG. 3: Magnetic (mg) excitations for q along [110]+[101]+[011] (left panel) and [111] (right panel) at T=150K, measured with unpolarised (circles) and half- polarized (triangles) neutrons. Full (empty) symbols refer to modes with main (weak) intensity. The hatched area correspond to calculated levels (see the text). Dashed lines are guide to the eyes for phonon branches measured at the same temperature.
0 5 10 15 20 25 30 Energy (meV)
0
100
200
300
400
500
600
N eu
tr on
c ou
nt s/3
0s
ζ=0.135ζ=0.15ζ=0.175ζ=0.2ζ=0.25
x=0.2 [111]
T=150K
(2)
(0)
FIG. 4: Energy spectra for several ζ along [111] at 150K, showing a dispersive mode (0) and a q-independent one (2).
EB, for which the nearest neighbour spins fluctuate in phase opposition, a half-period corresponds to one lattice spacing. Therefore, one expects that EB is a multiple of EL. The number of lattice spacings defining the size may be provided by the ratio EB/EL. For a size of two lattice spacings along [100], one expects 2 levels, with values at EL=EB/2 and EB, corresponding to confined half-wave and full-wave respectively. Along the other directions, if the cluster is isotropic or nearly cubic, one expects the same number of levels with the same ratio between the energies, the energy values being deduced from the above relations. We show now that the observations support this assumption. In Fig. 1, right
panel, level (1), E [100] B ≈18 meV, which defines the energy
of the first neighbor coupling along [100], and level (2) at half this value, are assigned respectively to the full-wave and half-wave, standing along [100] (or [010] equivalent)
4 �
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
J(m eV
)
c
0.00 0.05 0.10 0.15 0.20 0.25 x(Ca)
−1.0
−0.5
0.0
0.5
1.0
1.5
2.0
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
J
Ja,b
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
��
��
��
��
��
��
��
��
��
��
��
��
��
��
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
����������
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
a
b
a
ξ∼16 Αξ∼20 Α ξ∼30 Α
[001]
[010]
ξ∼8 Α x(Ca)
[100]
0.05<x<0.125 x=0.2x=0.17
CAF FMF
c
b
c
FIG. 5: Upper panel: Variation of the magnetic coupling values Ja,b (full symbols) and Jc (empty symbols) with x at 10K (circles) and 150K (squares). Lower panel: Schematic drawing showing hole-rich (empty) and hole-poor (hatched) media. Increasing x from left to right, it shows 2D F hole-rich spin clusters, 2D F hole-poor spin clusters and isotropic 3D hole-poor spin clusters. A half (dashed line) and a full (solid line) standing waves are shown inside a cluster.
confined in a cluster with 2 lattice spacings (ξ=8A˚). In the same way, level (2) at ≈9 meV and level (3) at half this value (≈4.5 meV) are respectively assigned to the full and half wave along c or [001]. Level (2) corresponds
therefore to 1/2E [100] B and E
[001] B values, superimposed
because of twinning. From the E [100] B and E
[001] B values,
Ja,b (1.12±0.1 meV) and Jc (0.56±0.05 meV) are determined. Within the assumption of isotropic or nearly-cubic shape, two levels at EB and EB/2 are also determined along the other symmetry directions. Due to the experimental uncertainty, they are shown by hatched area instead of narrow lines in Fig 3. In left panel, the sets (1), (2) and (3), (4), belong respectively to [110] and [101]+[011] directions, unequivalent, weighted by 1/3 and 2/3. The agreement with experimental values is good. It supports the proposed assignation of the levels and the assumption of isotropic shape for the clusters. In the same way, for x=0.17, Ja,b may be
determined from E [100] B (or E
[110] B =2E
[100] B ) and Jc, from
the difference between E [111] B and E
[110] B . These values
are reported in Fig. 5, upper panel, with those obtained at x=0.2. The small value of Jc is consistent with the 2D character found for the confined spin waves for this x=0.17 compound.
The existence of clusters with a typical size indicates that the magnetic characteristics differ inside from out- side the cluster. This implies a change of the charge density on the cluster’s scale. The nature, ”hole-poor”, suggested by the anisotropy of the coupling constants de- termined above, is supported by a comparison with obser-
vations at lower doping. In Fig. 5 upper panel, the values Ja,b and Jc, determined at 150K from the q-independent levels are reported with the J values obtained at 10K from the spin wave branch of SE type observed in the low dop- ing and undoped compounds[5]. Besides a small temper- ature effect (12% in the 10K-150K range), the variation of Ja,b(x) and Jc(x) is monotonous across the CAF/F boundary (Jc=0 at x=0.125), showing a tendency to be- come isotropic at x=0.225 (metallic state). This varia- tion implies a common origin for the coupling on the F and CAF sides, of SE type, characteristic of an orbital- ordered state with a roughly linear effect of the doping. This variation is close to theoretical predictions[12]. As in LaMnO3 or in the CAF state, this SE coupling is at- tributed to ”hole-poor” or without holes regions. There- fore, the cluster could consist of one Mn3+ with its first Mn3+ neighbourgs, the mobiles holes being confined at their boundaries. The small-q dispersion (0), which de- termines TC , could be attributed to a coupling induced by mobile-holes through the clusters. These clusters dif- fer from those found in the CAF state where the ”hole- poor” region corresponds to the matrix and the clusters, described as ”hole-rich” platelets[4]. The general evolu- tion with x in direct space is schematically shown in Fig. 5. Finally, a comparison with other techniques[7, 8] al- lows to estimate the lifetime τ of these clusters for x=0.2, 10−6s < τ < 10−9s, very large for the neutron probe. In conclusion, this study reports a quantitative descrip-
tion of the precursor state of the metallic phase which occurs for x≥0.225. An evolution of the shape and of the size of clusters with x is proposed, thanks to the exis- tence of confined spin waves which are observed for the first time in small F clusters. Very interestingly, these magnetic clusters have features (size, isotropy) very sim- ilar to those observed at larger x in CMR compounds around TC [2]. The present observations could be there- fore crucial for the understanding of the CMR properties. The authors are very indebted to L. P. Regnault and J.
Kulda for their help for experiments with field and polar- ized neutrons, and N. Shannon, T. Ziman, D. Khomskii and A. M. Oles´ for helpful discussions.
[1] A. Moreo et al. Science 283, 2034 (1999) [2] J. W. Lynn et al. Phys. Rev. Lett. 76, 4046 (1996) [3] J. M. De Teresa et al. Nature 386, 256 (1997) [4] M. Hennion et al. Phys. Rev. Lett. 81, 1957 (1998) and
Phys. Rev. B 61, 9513 (2000) [5] G. Biotteau et al. Phys. Rev. B 64 104421 (2001) [6] T. Okuda et al. Phys. Rev. B 61 8009 (2000) [7] G. Papavassiliou et al. Phys. Rev. Lett.91, 147205 (2003) [8] V. Chechersky et al. Phys. Rev. B 59 497 (1999) [9] V. Markovitch et al. Phys. Rev. B 66 094409 (2002) [10] P. V. Hendriksen et al. Phys. Rev. B 48, 7259 (1993) [11] Y. Yamada et al. Phys. Rev. B 62, 11600 (2000) [12] L. F. Feiner and A. M. Oles´ Physica B 259,796 (1999)
iv :c
on d-
m at
/0 40
86 11
v2 [
co nd
-m at.
str -el
] 11
Ja n 2
00 5
Confined spin-waves reveal an assembly of nanosize domains in ferromagnetic La1−xCaxMnO3 (x=0.17, 0.2)
M. Hennion1, F. Moussa1, P. Lehouelleur1, F. Wang1,2, A. Ivanov3, Y. M. Mukovskii4 and D. Shulyatev4 1 Laboratoire Le´on Brillouin, CEA-CNRS, CE-Saclay, 91191 Gif sur Yvette, France.
2 Institute of Physics Chinese Academy of Sciences, 100080 Beijing, China 3 Institut Laue-Langevin, 156 X, 38042 Grenoble Cedex 9, France 4 Moscow State Steel and Alloys Institute, Moscow 119991 Russia
We report a study of spin-waves in ferromagnetic La1−xCaxMnO3, at concentrations x=0.17 and x=0.2 very close to the metallic transition (x=0.225). Below TC , in the quasi-metallic state (T=150K), nearly q-independent energy levels are observed. They are characteristic of standing spin waves confined into finite-size ferromagnetic domains, defined only in (a, b) plane for x=0.17, and in all directions for x=0.2. They allow an estimation of the domain’s size, a few lattice spacings, and of the magnetic coupling constants inside the domains. These constants, anisotropic, are typical of an orbital-ordered state, allowing to characterize the domains as ”hole-poor”. The precursor state of the CMR metallic phase appears, therefore, as an assembly of small orbital-ordered domains.
PACS numbers: PACS numbers: 74.25.Ha 74.72.Bk, 25.40.Fq
In metallic La1−xCaxMnO3, for x≈0.3, the origin of the colossal magneto-resistance (CMR) which occurs at the metal-insulator transition is still a challenge for the- ories. The role played by an inhomogeneous ground state has been outlined by several theoretical[1] and ex- perimental works[2, 3]. In the semi-conducting parent LaMnO3, the magnetic coupling is of super-exchange (SE) type, ferromagnetic (F) in the (a, b) plane (Ja,b> 0) and antiferromagnetic (AF) along c (Jc<0), consistent with the orbital ordering which occurs at the Jahn-Teller transition TJT . Doping with holes induces a canted AF state (CAF). There, neutron diffuse scattering in- dicates the existence of magnetic inhomogeneities at- tributed to hole-rich clusters embedded in a hole-poor matrix[4], and a new spin-wave branch appears at small- q around F Bragg peaks. This new branch coexists with the spin-wave branch close to that of LaMnO3 (SE type of coupling), attributed to the hole-poor matrix[5], which disappears at x=0.125 (F state). The ferromagnetic state which occurs in the range 0.125≤x<0.225 shows a quasi-metallic behavior below TC and an insulating one at lower temperature[6]. There, several works have also proposed an inhomogeneous picture using NMR[7], Mossbauer[8] or magnetization measurements[9].
In the present paper, we report a determination of spin-waves in ferromagnets La1−xCaxMnO3 at x=0.17 and 0.2 very close to the ”true” metallic state (x≥0.225). Below TC , in the quasi-metallic state, q-independent en- ergy levels of magnetic origin are observed in the large q-range. The q-dependence of their intensity, their evolu- tion with temperature and magnetic applied field suggest that they can be ascribed to standing spin waves charac- teristic of a confinement into small domains. The shape and the size (a few lattice spacings) of the domains can be estimated as well as the magnetic coupling inside the domains. All the observations are consistent with a pic- ture of hole-poor and orbital-ordered domains.
Inelastic neutron experiments have been carried out at the reactor Orphe´e (Laboratoire Le´on Brillouin) and at the Institut Laue-Langevin, on triple axis spectrom- eters, using fixed final wave vector of neutrons varying from 1.05 to 4.1A˚−1, and appropriate filters. The TC values are 175K and 180K for x=0.17 and x=0.2 respec- tively. Both samples have been characterized by resis- tivity measurements. In La0.83Ca0.17MnO3, the high- temperature pseudo-cubic structure, with lattice param- eter a0=3.89A˚, becomes orthorhombic at TJT=240K, whereas such an orthorhombicity is unobservable for La0.8Ca0.2MnO3. All Q values (Q=q+τ and q=ζ) are given in reduced lattice units with pseudo-cubic indexa- tion. The two single crystals are twinned, so that each di- rection is superimposed on related symmetry directions. Data are fitted by a lorentzian shape weighted by a Bose factor. This analysis yields the energy E(q), the damping Γ(q) with Γ(q)/E(q)≈1/5 and the intensity.
We first consider La0.83Ca0.17MnO3. The spin wave modes obtained at T=150K, 100K and 10K along [001]+[010]+[001] directions are reported in Fig. 1, left panel. At 150K and 100K, one may distinguish two dis- persive modes labelled (1) and (3). They are very close to previous measurements at x=0.125 reported in or- thorhombic indexation[5]. This comparison allows us to assign the curve (1), dispersing up to ≈18 meV, to [100] or [010] directions, which have been shown to be equiv- alent, and the curve (3), with a down-turn at q=0.25, to the [001] or c direction. We mention that this down- turn, also observed at x=0.125 and in the CAF state, is surprising for a F coupling. It denotes an underlying AF coupling reminiscent of the orbital ordering of LaMnO3. The very new feature is the gap-opening in the dispersive curve (1). Within this gap lies a nearly q-independent energy level, (2), at a mean value of 4.5 meV at 150K. With decreasing temperature, all the energies increase. At 100K and below, a single dispersion curve is observed
2 0.00 0.25 0.500.00 0.25 0.50 0
5
10
15
20
En er
gy (m
eV )
[100] [001]
0.0 0.1 0.2 0.3 0.4 0.5 0
5
10
15
20
½
½
0.0 0.1 0.2 0.3 0.4 0.5 0
5
10
15
20
E ne
rg y
(m eV
)
[001] B
x=0.2
0.00
(1)
(0)
E
T=150K
,
[001]
E [100] B
+related symmetry[1+ζ,0,0]+related symmetry[1+ ζ,0,0]
(0)
(2)
(1)
(3)
(1)
(2)
T=10K
(2)
(3) (2) E B E B
[100]
0.0 0.1 0.2 0.3 0.4 0.5 0
5
10
15
20 T=17KT=50K T=50K
x=0.17 T=100KT=150K
(3) (2)
(1) (1)
(3)
FIG. 1: Magnetic excitations measured along [100]+[010]+[001] for x=0.17 (left panel) and x=0.2 (right panel) at several temperatures. Full (empty) circles correspond to modes with main (weak) intensity. In left panel, the solid and broken lines are guides to the eyes, and the vertical lines locate 1/8, 1/4 and 3/8 q values. In right-upper panel, the dashed line indicates the spin-wave dispersion of a virtual large-size domain (see the text).
for q<0.125, larger gaps open at q=1/8, 1/4 and 3/8 and the level labelled (2) appears modulated with a maxi- mum close to 3/8. Measurements of spin waves along [110] and [111] at 150K also reveal q-independent energy levels lying at nearly the same values in the two direc- tions.
These observations call for the following comments. A q-independent energy level indicates localized excita- tions. The existence of a gap-opening in the branch (1) points out its interaction with level (2). This implies that level (2) belongs to the same macroscopic domain, namely to the (a, b) plane, but not to the c direction. The localized excitations or standing spin waves are con- fined within these planes, exhibiting a two-dimensionnal (2D) character. This is consistent with observations along [110] and [111], as indicated below. Since the dis- persionless level appears above qcr≈0.2, an approximate size ξ=a0/qcr≈ 20A˚ may be provided for the domain. A more quantitative determination of the size can be done by considering the energy value of the q-independent level. Following the arguments developped below for the
ζ=0.5 T=150K
ζ=0.275 x=0.2, [100]+[010]+[001]
(3) (2)
(1)
ph
ph
ph
(1) (2)
(3)
0 4 8 12 16 20 24 28 Energy (meV)
0
100
200
300
N eu
tr on
c ou
nt s/
12 0s
0 4 8 12 16 20 24 28 Energy (meV)
0
50
100
150
200
250
300
C ou
nt s/1
20 s
150K T=210K
FIG. 2: Examples of energy spectra along [100]+[010]+[001] at H=0, after application of a magnetic field H=2T. Left panel: ζ=0.275 at T=150K (circles) and 210K (triangles, shifted by 100 neutron counts). Right panel: ζ=0.5. The modes (1), (2) and (3) correspond to spin-wave branches in- dicated in Fig. 1, and ”ph” is a phonon mode.
x=0.2 compound, the ratio, ≈1/4, between the energies of the level (2) and of the zone boundary, determines a size of 4 lattice spacings (16A˚) along [100] (or [010]). In the insulating state and specially below 100K, this pic- ture is no longer valid. We may explain the origin of gaps and the modulation of the level (2) by assuming an un- derlying periodicity of 4a0 in the (a, b) plane. At the new zone boundaries of the super-cell, the dispersion is folded and new gaps occur. This description of the spin dynam- ics is likely related to observations in La7/8Sr1/8MnO3, where static superstructures are observed[11].
The study of La0.8Ca0.2MnO3 reveals a very interest- ing evolution with x. This sample has been first stud- ied in zero field. Preliminary observations at 17K along [001]+[010]+[001] have been reported[5]. Results are de- scribed with decreasing temperature. At 150K, unlike the x=0.17 case, a single dispersion branch, labelled (0), is observed in the small q range up to q≈0.3, indicat- ing that the [100] (a and b) and [001] ( c) directions are now equivalent and the dispersion isotropic (see Fig 1- right upper panel). For q≥ 0.3, the magnetic intensity is shared between several nearly q-independent levels de- fined by their mean values: a broad energy level, (1), at 18±1 meV, likely coupled to the lower optical phonon mode, a level (2) at 9±0.5 meV and an overdamped level, quasielastic. A magnetic field (H=2T) applied at 150K, suppresses the quasi-elastic intensity attributed to spin fluctuations, and reduces the energy linewidth of the dis- persionless magnetic levels. When removing the field at 150K, the spin fluctuations remain considerably reduced, and the levels keep the characteristics observed in applied field. A low-energy mode (3) is now defined at a mean value 4.5 meV (see Fig 1) and the mode (1) appears sep- arated from the phonon mode (24 meV), not reported
3 in Fig 1. Corresponding energy spectra are displayed in Fig. 2, showing the q-dependence of the intensity. As q increases from q=0.275, to 0.5, the intensity of the modes (2) and (3) is transferred towards the energy mode (1) (18 meV), so that the intensity follows the q-dependence of a F dispersion curve. Very interestingly, these modes persist above TC (180K). This is shown by the spectrum at 210K in Fig. 2, left panel, where a slight decrease of the two energy levels (2) and (3) is observed whereas a quasielastic intensity (spin fluctuations) increases above TC . In contrast, the small-q dispersive branch (0) renor- malizes to zero at TC . By lowering temperature in zero field, the lower-energy mode (3) nearly disappears and level (2) becomes more dispersed, which reduces the en- ergy gaps (Fig. 1, lower-right panel).
Spin wave measurements along [110]+[101]+[011] and [111] directions are reported in Fig 3 at T=150K. At small q, the quadratic dispersion ω=Dq2, (0), defines an isotropic stiffness constant (D=30±1meVA˚2 becom- ing 48meVA˚2 at 10K). At larger q values, the magnetic intensity is transferred to nearly q-independent energy levels. Along [111], a level (2) is determined at 23±1 meV and a slightly dispersive one (1) around 46±2 meV. Fig 4 displays the balance in intensity between the modes (0) and (2) as q varies. Along [110]+[101]+[011], the en- tanglement between magnon and phonon branches (see dashed lines in Fig. 3-left panel), led us to use polar- ized neutrons (PN) which allow to identify magnons and phonons. Several levels have been determined, more or less dispersive. Dispersion could have been induced by the applied magnetic field (H=2T) which aligns macro- scopic domains. Excitations at energy larger than 35 meV could not been detected because of a too weak inten- sity. By lowering temperature, the energy levels slightly increase.
These overall observations can be explained as it fol- lows. Unlike La0.83Ca0.17MnO3 where the confined spin waves have been assigned to the (a, b) plane only (2D character), for x=0.2, the isotropy of the small-q disper- sion suggests that all directions are concerned. The levels are attributed to standing spin-waves inside 3D clusters. Let us use the energy of the propagating wave for a 3D Heisenberg model (see dotted line in Fig. 1 upper-right panel):
E(qx,qy,qz)=8S{Ja,b[1-0.5(cos(2piqx)+cos(2piqy))]+ Jc/2[1-cos(2piqz)]}, where qx, qy, qz are related to a, b, c directions. This expression yields relations between Ja,b, Jc and the zone boundary energies
EB: E [100] B =E
[010] B =8SJa,b, E
[001] B =8SJc, E
[110] B =2E
[100] B ,
E [101] B =E
[011] B =E
[100] B +E
[001] B and E
[111] B =2E
[100] B +E
[001] B
These latter relations are also valid for finite-size clus- ters. According to predictions[10], the cluster’s size ξ is related to the level with the lowest-energy value EL. In the case of small clusters, it corresponds to a ”half-wave” or a half-period of wave. For the confined wave of energy
0.0 0.1 0.2 0.3 0.4 0.5 0
10
20
30
40
50
0.0 0.1 0.2 0.3 0.4 0.5 0
10
20
30
40
50 mg mg (PN)
�����
�����
�����
�����
�����
�����
�������������
�������������
�������������
�������������
�������������
�������������
�������������
�������������
�������������
�������������
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
2E +EB
(3)
(2)
B
(1)
+1/2E
Ca 20%
B [100]
EB
[100]
[001]
Ca 20%
(0)
[1+ζ,ζ,0] +related symmetry [1+ +related symmetryζ,ζ,ζ]
������
������
������
������
������
������
En er
gy (m
eV )
T=150K
(2)
(1)
(0)
T=150K [001]
(4)
FIG. 3: Magnetic (mg) excitations for q along [110]+[101]+[011] (left panel) and [111] (right panel) at T=150K, measured with unpolarised (circles) and half- polarized (triangles) neutrons. Full (empty) symbols refer to modes with main (weak) intensity. The hatched area correspond to calculated levels (see the text). Dashed lines are guide to the eyes for phonon branches measured at the same temperature.
0 5 10 15 20 25 30 Energy (meV)
0
100
200
300
400
500
600
N eu
tr on
c ou
nt s/3
0s
ζ=0.135ζ=0.15ζ=0.175ζ=0.2ζ=0.25
x=0.2 [111]
T=150K
(2)
(0)
FIG. 4: Energy spectra for several ζ along [111] at 150K, showing a dispersive mode (0) and a q-independent one (2).
EB, for which the nearest neighbour spins fluctuate in phase opposition, a half-period corresponds to one lattice spacing. Therefore, one expects that EB is a multiple of EL. The number of lattice spacings defining the size may be provided by the ratio EB/EL. For a size of two lattice spacings along [100], one expects 2 levels, with values at EL=EB/2 and EB, corresponding to confined half-wave and full-wave respectively. Along the other directions, if the cluster is isotropic or nearly cubic, one expects the same number of levels with the same ratio between the energies, the energy values being deduced from the above relations. We show now that the observations support this assumption. In Fig. 1, right
panel, level (1), E [100] B ≈18 meV, which defines the energy
of the first neighbor coupling along [100], and level (2) at half this value, are assigned respectively to the full-wave and half-wave, standing along [100] (or [010] equivalent)
4 �
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
�
J(m eV
)
c
0.00 0.05 0.10 0.15 0.20 0.25 x(Ca)
−1.0
−0.5
0.0
0.5
1.0
1.5
2.0
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
J
Ja,b
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
�����������
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
��
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
������
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
��
��
��
��
��
��
��
��
��
��
��
��
��
��
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
����������
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
�����
a
b
a
ξ∼16 Αξ∼20 Α ξ∼30 Α
[001]
[010]
ξ∼8 Α x(Ca)
[100]
0.05<x<0.125 x=0.2x=0.17
CAF FMF
c
b
c
FIG. 5: Upper panel: Variation of the magnetic coupling values Ja,b (full symbols) and Jc (empty symbols) with x at 10K (circles) and 150K (squares). Lower panel: Schematic drawing showing hole-rich (empty) and hole-poor (hatched) media. Increasing x from left to right, it shows 2D F hole-rich spin clusters, 2D F hole-poor spin clusters and isotropic 3D hole-poor spin clusters. A half (dashed line) and a full (solid line) standing waves are shown inside a cluster.
confined in a cluster with 2 lattice spacings (ξ=8A˚). In the same way, level (2) at ≈9 meV and level (3) at half this value (≈4.5 meV) are respectively assigned to the full and half wave along c or [001]. Level (2) corresponds
therefore to 1/2E [100] B and E
[001] B values, superimposed
because of twinning. From the E [100] B and E
[001] B values,
Ja,b (1.12±0.1 meV) and Jc (0.56±0.05 meV) are determined. Within the assumption of isotropic or nearly-cubic shape, two levels at EB and EB/2 are also determined along the other symmetry directions. Due to the experimental uncertainty, they are shown by hatched area instead of narrow lines in Fig 3. In left panel, the sets (1), (2) and (3), (4), belong respectively to [110] and [101]+[011] directions, unequivalent, weighted by 1/3 and 2/3. The agreement with experimental values is good. It supports the proposed assignation of the levels and the assumption of isotropic shape for the clusters. In the same way, for x=0.17, Ja,b may be
determined from E [100] B (or E
[110] B =2E
[100] B ) and Jc, from
the difference between E [111] B and E
[110] B . These values
are reported in Fig. 5, upper panel, with those obtained at x=0.2. The small value of Jc is consistent with the 2D character found for the confined spin waves for this x=0.17 compound.
The existence of clusters with a typical size indicates that the magnetic characteristics differ inside from out- side the cluster. This implies a change of the charge density on the cluster’s scale. The nature, ”hole-poor”, suggested by the anisotropy of the coupling constants de- termined above, is supported by a comparison with obser-
vations at lower doping. In Fig. 5 upper panel, the values Ja,b and Jc, determined at 150K from the q-independent levels are reported with the J values obtained at 10K from the spin wave branch of SE type observed in the low dop- ing and undoped compounds[5]. Besides a small temper- ature effect (12% in the 10K-150K range), the variation of Ja,b(x) and Jc(x) is monotonous across the CAF/F boundary (Jc=0 at x=0.125), showing a tendency to be- come isotropic at x=0.225 (metallic state). This varia- tion implies a common origin for the coupling on the F and CAF sides, of SE type, characteristic of an orbital- ordered state with a roughly linear effect of the doping. This variation is close to theoretical predictions[12]. As in LaMnO3 or in the CAF state, this SE coupling is at- tributed to ”hole-poor” or without holes regions. There- fore, the cluster could consist of one Mn3+ with its first Mn3+ neighbourgs, the mobiles holes being confined at their boundaries. The small-q dispersion (0), which de- termines TC , could be attributed to a coupling induced by mobile-holes through the clusters. These clusters dif- fer from those found in the CAF state where the ”hole- poor” region corresponds to the matrix and the clusters, described as ”hole-rich” platelets[4]. The general evolu- tion with x in direct space is schematically shown in Fig. 5. Finally, a comparison with other techniques[7, 8] al- lows to estimate the lifetime τ of these clusters for x=0.2, 10−6s < τ < 10−9s, very large for the neutron probe. In conclusion, this study reports a quantitative descrip-
tion of the precursor state of the metallic phase which occurs for x≥0.225. An evolution of the shape and of the size of clusters with x is proposed, thanks to the exis- tence of confined spin waves which are observed for the first time in small F clusters. Very interestingly, these magnetic clusters have features (size, isotropy) very sim- ilar to those observed at larger x in CMR compounds around TC [2]. The present observations could be there- fore crucial for the understanding of the CMR properties. The authors are very indebted to L. P. Regnault and J.
Kulda for their help for experiments with field and polar- ized neutrons, and N. Shannon, T. Ziman, D. Khomskii and A. M. Oles´ for helpful discussions.
[1] A. Moreo et al. Science 283, 2034 (1999) [2] J. W. Lynn et al. Phys. Rev. Lett. 76, 4046 (1996) [3] J. M. De Teresa et al. Nature 386, 256 (1997) [4] M. Hennion et al. Phys. Rev. Lett. 81, 1957 (1998) and
Phys. Rev. B 61, 9513 (2000) [5] G. Biotteau et al. Phys. Rev. B 64 104421 (2001) [6] T. Okuda et al. Phys. Rev. B 61 8009 (2000) [7] G. Papavassiliou et al. Phys. Rev. Lett.91, 147205 (2003) [8] V. Chechersky et al. Phys. Rev. B 59 497 (1999) [9] V. Markovitch et al. Phys. Rev. B 66 094409 (2002) [10] P. V. Hendriksen et al. Phys. Rev. B 48, 7259 (1993) [11] Y. Yamada et al. Phys. Rev. B 62, 11600 (2000) [12] L. F. Feiner and A. M. Oles´ Physica B 259,796 (1999)
Comments
