arxiv:cond-mat/ v1 [cond-mat.dis-nn] 28 Dec 2004

Transcription

1 arxiv:cond-mat/ v1 [cond-mat.dis-nn] 28 Dec 2004 EFFECTIVE CONDUCTIVITY OF 2D ISOTROPIC TWO-PHASE SYSTEMS IN MAGNETIC FIELD S.A.Bulgadaev 1, F.V.Kusmartsev 2 Landau Institute for Theoretical Physics, Chernogolovka, Moscow Region, Russia, Department of Physics, Loughborough University, Loughborough, LE11 3TU, UK Using the linear fractional transformation, connecting effective conductivities ˆσ e of isotropic two-phase systems with and without magnetic field, explicit approximate expressions for ˆσ e in a magnetic field are obtained. They allow to describe ˆσ e of various inhomogeneous media at arbitrary phase concentrations and magnetic fields. The obtained results are applicable for different two-phase systems regular and nonregular as well as random), satisfying the symmetry and self-duality conditions. All these results admit a direct experimental checking. 1. Introduction PACS: Ak, Ng, Tm, r Last time, under investigation of magneto-resistive properties of new materials, which are connected with the high-temperature superconductivity such as oxide materials with the perovskite type structure), it was established that they often have unusual transport properties. For example, the magnetoresistance becomes very large the so called colossal magnetoresistance in such materials as manganites) [1] or grows approximately linearly with magnetic field up to very high fields in silver chalcogenides) [2]. There is an opinion that these properties take place due to phase inhomogeneities of these materials [2]. For this reason a calculation of the effective conductivity σ e of inhomogeneous heterophase systems without and with magnetic field at arbitrary partial conductivities and phase concentrations is very important problem. Unfortunately, the existing effective medium approximations EMA) cannot give an explicit simple formulas for magneto-resistivity convenient for description of the experimental results in a wide range of partial parameters even for two-phase random systems [3]. Such formulas can be obtained for some systems only in high magnetic field [3, 4]. The situation is much better for 2D inhomogeneous systems in the perpendicular magnetic field. Here, due to their exact duality properties, not only a few exact results have been obtained for the effective conductivity of random inhomogeneous systems [5-10], but also the transformations, connecting effective conductivities of isotropic two-phase inhomogeneous systems with and without magnetic field, have been constructed [9, 11]. This transformation let us call it a magnetic one), in principle, permits to obtain the explicit expressions for σ e in a magnetic field, if the corresponding expressions are known for σ e without magnetic field. Such explicit approximate expressions for σ e have been obtained in 1 bulgad@itp.ac.ru 2 F.Kusmartsev@lboro.ac.uk 1

2 some limiting cases, for example, at small concentration of one phase, for weakly inhomogeneous media [11], for inclusions of super- or non-conducting phases [6-8,11]. In this letter we will present the explicit approximate expressions for σ e in a wide region of concentrations and partial conductivities, using the full dual transformation recently constructed in [10] and the corresponding expressions for σ e at H = 0 from [12]. 2. Dual transformation between systems with H 0 and H = 0 The effective conductivity of two-phase isotropic systems in a magnetic field has the following form ˆσ = σ ik = σ d δ ik + σ t ǫ ik, σ d H) = σ d H), σ t H) = σ t H), 1) here δ ik is the Kronecker symbol, ǫ ik is the unit antisymmetric tensor. The effective conductivity ˆσ e of 2-phase random or regular symmetric systems with the partial conductivities σ id, σ it i = 1, 2) we assume that σ id 0) and concentrations x i must be a symmetric function of pairs of arguments ˆσ i, x i ) and a homogeneous a degree 1) function of σ di,ti. For this reason it is invariant under permutation of pairs of partial parameters ˆσ e ˆσ 1, x 1 ˆ, x 2 ) = ˆσ e ˆ, x 2 ˆσ 1, x 1 ). 2) The effective conductivity of 2-phase systems must reduce to some partial ˆσ i, if x i = 1 i = 1, 2). Futher it will be more convenient to use the complex representation for coordinates and vector fields [13] z = x + iy, j = j x + ij y, e = e x + ie y, σ = σ d + iσ t. The complex conductivity transforms as σ = Tσ) = cσ ib idσ + a, 3) which generalizes the inversion transformations of systems without magnetic field [5]. Thus, the dual transformations DT) in systems with magnetic field have a more richer structure due to the fact that they are connected with some subgroup of group of linear fractional transformations, conserving the imaginary axis in the complex conductivity plane [5,7,9]. The transformation of σ e depends on 3 real parameters since one of 4 parameters can be factored due to the fractional structure of T). There are various ways to choose 3 parameters. In our treatment it will be convenient to factor d. This gives 3 parameters ā = a/d, b = b/d, c = c/d, determining a transformation T. The transformation 3) has the following form in terms of conductivity components σ d and σ t σ d = σ ac + bd d dσ d ) 2 + a + dσ t ) 2 = cσ ā + b/ c d σ d ) 2 + ā + σ t ) 2, 2

3 σ t = cdσ2 d + a + dσ t)cσ t b) dσ d ) 2 + a + dσ t ) 2 = c σ2 d + ā + σ t)σ t b/c) σ d ) 2 + ā + σ t ) 2. 4) These DT allow to construct the transformation, connecting effective conductivities of two-phase systems with magnetic field and without it. Analogous connections have been found firstly on a basis of solutions of the corresponding Laplace and boundary equations [11], but it has some ambiguity problems with a determination of the parameters of the artificial system. Later, this transformation has been constructed directly from the DT under two simplifying conjectures in [9]. Recently, it was constructed in the full form, using all 3 parameter and explicitly reproducing the known exact values for σ e, in [10]. The parameters of such transformation let us call it T h, we will also omit the bars over its parameters) a, b = b/c, c depend on the partial conductivities and have the following form a ± = c = 2 σ 1 2 ± B 2σ 1t t ), b ± = σ ± B, 2σ 1t t ) B = [σ 1t t ) 2 + σ 1d d ) 2 ][σ 1t t ) 2 + σ 1d + d ) 2 ], 5) where σ i 2 = σid 2 + σ2 it, and, evidently, B 0. The diagonal or real) parts of σ i transform under T h as σ id = σ id ca + b ) σ id ) 2 + a + σ it ) 2 = cσ id σ ai, σ ai = a + σ it. 6) The parameters a, b, c satisfy also the additional relations [10] [ ) ] 2 1/2 σ1t t A = = a + b ) ) 1/2, 7) σ 1d + d σ 1d d + σ 1d t + d σ 1t σ 1d + d = c σ 2 e ab σ e 2 + a = cσ 1dd b /a), 8) 2 σ 1d d + which ensure the reproduction of the exact results for σ e at the equal phase concentrations x 1 = x 2 = 1/2. The relations 7)and 8) give us a highly nontrivial check of a self-consistency of the transformation. The direct check of them is a rather complicated task. Fortunately, as it was noted in [9] and shown in [10], T h transforms a circumference see fig.1.) with a radius R, centered at the imaginary axis at z 0 and passing through the two points σ 1 and, where R = a + b 2 = B 2σ 1t t ), z a + b 0 = = σ σ 1t t ), 9) into the real axis. All equalities necessary for fulfillment of 7),8) correspond to the known equalities between lengths of various chords, their projections on a diameter and the radius of the circles. For example, the equality 7) takes the form 2Rσ a A = σa 2 + σ = 1/1 δ/2rσ a ), σa 2 1dσ =, 10) 2d 3

4 ImΣ b Σ 2 Σ a Σ 1 ReΣ z 0 Σ x a Fig.1. A schematic picture of the chord, defining a geometrical sense of the factor A from the expression for the real part of the exact σ e. where a is a squared geometrical average of σ ai, i = 1, 2) which always have the same sign), and δ defines a difference between x = 2Rσ a a, a squared real projection of the chord, having an imaginary projection σ a, and a squared geometric average of σ id δ = x σ 1dd. 11) 3. Inhomogeneous systems with compact inclusions Now, having all formulas for the transformation, one can construct the explicit expressions for σ e of inhomogeneous system in a magnetic field. Below, in order to simplify a view of the subsequent formulas, we will omit the subindex d in the partial diagonal parts σ id i = 1, 2). Firstly we consider the case of compact inclusions of one phase into another. The corresponding effective conductivity σ e has the next form x 1 = x, x 2 = 1 x) [12] σ e {σ}, {x}) = σ x 1σ 1 x 2 12) Substituting into 12) the primed conductivities one obtains the primed effective conductivity ) x ) 1 x σ ed{σ}, {x}) = c 13) Then, substituting this into 4), one obtains for the diagonal part σ ed of the effective conductivity of these systems in a magnetic field the following expres- 4

5 sion ) x1 ) x2 σ ed {σ}, {x}) = σ ed ac + b) a b/a) σ ed )2 + a 2 = ) 2x1 ) 2x2. 14) One can check, using the relation 7), that 14) correctly reduces to σ i, when x i = 1 and to the exact formula for σ e at x 1 = x 2 = 1/2. For example, at x 1 = 1, x 2 = 0 14) reduces to σ ed {σ}, x 1 = 1) = σ 1 a b/a) σ1 ) 2 = σ 1, where at the last step we have used the relation 7) at σ 1 =. For the transverse part σ et one obtains ) 2x1 ) 2x2 σ et {σ}, {x}) = c σ ed )2 ab b/a + a σ ed )2 + a 2 = ) 2x1 ) 2x2. 15) Again, using now the relation 8), one can check that 15) correctly reproduces boundary values at x i = 1 i = 1, 2) as well as the exact value at equal concentrations x 1 = x 2 = 1/2. 4. Systems with a random parquet structure In this section we present analogous formulas for 2D isotropic inhomogeneous systems with the structure of the inhomogeneities of the random parquet type. The effective conductivity of such systems without a magnetic field is [12] σ e {σ}, {x}) = σ / σ 1 = ) 1/2 x1 σ x 2 σ 1 16) x 1 + x 2 σ 1 Then the diagonal and transverse parts of the primed effective conductivity σ e has the form ) σ1 ) 1/2 σ ed{σ}, σ1 {x}) = c x x σ2 2 ) ) x σ2 1 + x σ1 2 Its substitution into 4) gives for the diagonal part σ ed σ ed {σ}, {x}) = a b/a) σ1 ) 1/2 σ1 x1 σ1 σ1 x1 σ a ) 1/2 ) ) σ2 )1/2 σ +x 2 a1 σ ) a2 ) σ2 σ1 x 1 σ +x 2 a2 ) ) σ2 σ +x 2 a1 σ ) a2 ) σ2 σ1 x 1 σ +x 2 a2 17) ) 18) 5

6 One can check, using the relation 7) that 18) reduces to the right values at x i = 1, i = 1, 2) as well as for x 1 = x 2 = 1/2. For the transverse part σ et one obtains ) ) ) ) b/a 2 + σ1 σ x 1 2 σ +x 2 a1 σ ) a2 ) σ2 σ1 x 1 σ +x 2 σ et {σ}, {x}) = a a2 ). 19) ) x σ1 σ 1 2 ) ) σ +x 2 a1 σ ) a2 ) σ2 σ1 x 1 σ +x 2 a2 Again, using now the relation 8), one can check that 19) correctly reproduces boundary values as well as the exact value at equal concentrations x 1 = x 2 = 1/2. 5. Effective medium approximation in a magnetic field In this section we find out for a completness a magnetic transformation for the traditional effective medium approximation EMA) for the effective conductivity. The EMA for σ e in inhomogeneous two-phase self-dual systems without a magnetic field has the form [14] σ e {σ}, {x}) = x 1 2 )σ + x 1 2 )2 + σ 1 ), 20) where σ = σ 1 ). Then the primed effective conductivity will be σ e{σ}, {x}) = c x 1 2 )σ a + x 1 2 )2 a + σ 1 ), 21) here σ a = σ1 σ2. Substituting 21) into 4), one obtains for the diagonal part σ ed the next expression ) a b/a) x 1 2 )σ a + x 1 2 )2 σa 2 + σ1 σ ed {σ}, {x}) = ) 2. 22) x 1 2 )σ a + x 1 2 )2 σa 2 + σ1 One can check, using the relation 7), that 22) correctly reduces to σ i, when x i = 1 and to the exact formula for σ e at x 1 = x 2 = 1/2. For example, at x 1 = 1, x 2 = 0 22) reduces to σ ed {σ}, x 1 = 1) = a b/a)σ 1 a σ2 1 = σ 1. Analogously, substituting 21) into 4), one obtains for the transverse part of the effective conductivity σ et {σ}, {x}) = a b/a 2 + x 1 2 )σ a + x 1 2 )σ a + x 1 2 )2 a + σ1 ) 2 x 1 2 )2 a + σ1 ) 2. 23) 6

7 Again, using now the relation 8), one can check that 23) correctly reproduces boundary values as well as the exact value at equal concentrations x 1 = x 2 = 1/2. It follows from the property of the transformation T h that σ e belongs to the circumference 9) for all phase concentrations x [0, 1] [9]. Thus the effective conductivity maps the concentration segment [0, 1] into the corresponding arc of this circumference, connecting the points σ 1 and. The different expressions for σ e correspond to the different mappings. As an example, we present on the fig.2. the plots of the normalized diagonal part σ ed /σ 1d for three explicit expressions obtained in the sections 3,4,5 let us denote these expressions in an order of the obtaining 1,2,3), when σ it have the different signs. In this case the difference in the behaviour of σ ed with H = 0 [12] and H 0 is the most demonstrative. In particular, σ ed increases more rapidly with an increase of x and shows a saturation behaviour at x near 1. Thus, we have found the desired formulas. Substituting into 14),15), 17),18),22) and 23) the known functions σ id H), σ it H), i = 1, 2), one obtains the explicit expressions for the diagonal σ ed H) and transverse σ et H) parts of the effective conductivity of inhomogeneous systems in a magnetic field. It follows from these formulas that the corresponding σ e satisfy all necessary properties, enumerated above. More detailed analysis of the obtained formulas will be done in subsequent paper. 6. Conclusion Using the exact duality transformation we have found three explicit approximate expressions for the effective conductivity of 2D isotropic two-phase systems in a magnetic field. The three plots of the dependence of σ ed on the phase concentration at different values of the inhomogeneity and signs of the partial transverse conductivities are constructed. The obtained results can be applied for a description of σ e of various two-phase systems regular and nonregular as well as random), satisfying the symmetry and self-duality conditions, in a wide range of partial conductivities and at arbitrary concentrations and magnetic fields. All these results admit a direct experimental checking. Acknowledgments The authors are thankful to Prof.A.P.Veselov for very useful discussions of some mathematical questions. This work was supported by the RFBR grants , , and by the Royal Society UK) grant 2004/R4-EF. References [1] G.Allodi et al., Phys.Rev.B ) 6036; M.Hennion et al., Phys.Rev.Lett ) 1957; Y.Moritomo et al., Phys.Rev. B )

8 [2] R.Xu et al., Nature ) 57. [3] D.J.Bergmann, D.Stroud, Phys.Rev.B ) [4] Yu.A.Dreizin, A.M.Dykhne ZhETF ) 242 Sov.Phys. JETP ) 127); I.M.Kaganova, M.I.Kaganov, cond-mat/ ). [5] A.M.Dykhne, ZhETF ) 641, JETP ) 348). [6] A.L.Efros, B.I.Shklovskii, Phys.Stat.Sol. b), ) 475. [7] B.I.Shklovskii, ZhETF ) 288. [8] D.G.Stroud, D.J.Bergmann, Phys.Rev.B ) 447. [9] G.W.Milton, Phys.Rev. B ) [10] S.A.Bulgadaev, F.V.Kusmartsev, Phys.Lett. A 2005), in press. [11] B.Ya.Balagurov, ZhETF, ) [12] S.A.Bulgadaev, Pis ma v ZhETF, ) 615; Europhys.Lett ) 482, cond-mat/ , to be published. [13] L.D.Landau, E.M.Lifshitz, Electrodynamics of condensed media, Moscow, 1982 in Russian). [14] S.Kirkpatrick, Rev.Mod.Phys )

9 Σ ed Σ 1 d x a) Σ ed Σ 1 d x b) Σ ed Σ 1 d x c) Fig.2. The plots of the normalized diagonal part σ ed /σ 1d for three explicit expressions obtained above from the left to the right, respectively, 3,1,2) at σ 1t/σ 1d = 0.3, t/σ 1d = 0.7 and at the three different ratios of d /σ 1d a) d /σ 1d = 0.1, b) d /σ 1d = 0.01, c) d /σ 1d =