arxiv: v1 [math-ph] 20 Jul 2007

Transcription

1 Many-body wave scattering by small bodies and applications arxiv: v [math-ph] 20 Jul 2007 A. G. Ramm (Mathematics epartment, Kansas St. University, Manhattan, KS 66506, USA and TU armstadt, Germany) PACS: Kf, tg, dt MSC: 35J05, 35J0, 70F0, 74J25, 8U40, 8V05 Key words: many-body problem, wave scattering by small bodies, small particles, smart materials. Abstract A rigorous reduction of the many-body wave scattering problem to solving a linear algebraic system is given bypassing solving the usual system of integral equation. The limiting case of infinitely many small particles embedded into a medium is considered and the limiting equation for the field in the medium is derived. The impedance boundary conditions are imposed on the boundaries of small bodies. The case of Neumann boundary conditions (acoustically hard particles) is also considered. Applications to creating materials with a desired refraction coefficient are given. It is proved that by embedding suitable number of small particles per unit volume of the original material with suitable boundary impedances one can create a new material with any desired refraction coefficient. The governing equation is a scalar Helmholtz equation, which one obtains by Fourier transforming the wave equation. Introduction This paper is a continuation of [4] and uses some of the results from [9], [7], [0], [5], [6]. Applications of our theory to creating materials with desired refraction coefficient, including negative refraction, are discussed in [], [2], [3], [7]. Wave scattering by small bodies is a classical branch of science: it was originated by Rayleigh in 87. In [2] one finds a discussion of wave scattering by a small particle. In [] there is a review of the low frequency scattering theory and

2 formulas for scattering by small balls and ellipsoids are given. In [9] the theory is developed for small bodies of arbitrary shapes. In [4] the many-body scattering problem was reduced to solving linear algebraic systems bypassing the usual study of a system of integral equations. In this paper we apply the approach proposed in [4] and study the limiting behavior of the scattering solution when the number of small bodies tends to infinity in such a way that the characteristic size a of the small particles is related to their number M so that M = O( ) in Theorem 2, and a M = O( ) in Theorem 3. Sufficient conditions for convergence of the scattering a 3 solution in this limiting process are given. We prove that these conditions are, in some sense, also necessary for convergence. The limit of the scattering solution is a function, which satisfies some differential or integral-differential equations. These equations describe the behavior of the wave field in the new medium, obtained in the limit. There is a large literature on the calculation of the effective dielectric permittivity and magnetic permeability of the composite materials (Maxwell-Garnett and Bruggeman recipes and their numerous versions, and newer theories [8], [4]). In the literature mostly a randomly uniform distribution of the inclusions is assumed and the resulting homogenized medium is described by effective constant dielectric permittivity and magnetic permeability, which can be tensors. In this work the propagation and scattering of scalar waves are discussed, and the homogenized medium is described not by a constant refraction coefficient, but by a refraction coefficient which is a function of spatial variables. Let us formulate the problem. Consider first a bounded domain R 3 filled with a material with a known refraction coefficient n 0 (x). The governing equation is: L 0 u 0 := ( 2 + k 2 n 0 (x) ) u 0 = 0 in R 3. (.) We assume that n 0 (x) = in = R 3 \, k = const > 0, and n 0 = max x n 0 (x) <. The operator L 0 can be written as a Schrödinger operator: L 0 = 2 + k 2 q 0 (x), q 0 (x) := k 2 [ n 0 (x)], (.2) and q 0 = 0 in. One has n 0 (x) = k 2 q 0 (x), so there is a one-to-one correspondence between n 0 (x) and q 0 (x). If n 0 (x) is known, then one knows the scattering solution: L 0 u 0 = 0 in R 3, u 0 (x) = e ikα x + A 0 (β, α) eikr r + o( ) x, r = x, β := r r. (.3) The coefficient A 0 (β, α) is called the scattering amplitude, the unit vector α S 2 is given, α is the direction of the incident plane wave e ikα x, S 2 is the unit sphere 2

3 in R 3, β S 2 is the direction of the scattered wave, k > 0 is a wave number, which we assume fixed throughout the paper. By this reason we do not show the k-dependence of A and u 0. Let G(x, y) be the resolvent kernel of L 0 satisfying the radiation condition (or the limiting absorption principle): L 0 G(x, y) = δ(x y) in R 3. (.4) This function G(x, y) is known because q 0 (x) is known. Consider now the scattering problem for many small bodies m embedded in, m M: L 0 u M = 0 in R 3 \ u M M m, (.5) = u 0 + A M (β, α) eikr r + o( r ), r = x, x r = β, (.6) u M N = ζ mu M on := m, m M, (.7) where u 0 is the solution ot the scattering problem (.3). Here N is the normal to pointing out of m, ζ m is a complex number, the boundary impedace, Im ζ m 0, is uniformly C,λ with respect to m, m M. By C,λ surface we mean the surface with local equation x 3 = f(x, x 2 ), where f C,λ, λ > 0. We assume throughout this paper that n 0 ka, d a, (.8) where a = 2 max m diam m, d = min m j dist( m, j ). (.9) By V m := m the volume of m is denoted, and by the surface area of is denoted. One can prove (see Section 3) that problem (.5) (.7) has at most one solution if Im ζ m 0, m M, and Im q 0 (x) 0. We look for the solution to problem (.5) (.7) of the form u M (x) = u 0 (x) + G(x, s)σ m (s)ds, (.0) where σ m should be found from the boundary conditions (.7). For any σ m the function (.0) solves equation (.5) and satisfies condition (.6): A M (β, α) = 4π 3 u 0 (s, β)σ m (s)ds. (.)

4 Formula (.) follows from (.6), (.0) and the Ramm s lemma ([8], formulas (5..3), (5..36)): ( G(x, y) = eik x ) 4π x u 0(y, α) + o, x, x x x = α, (.2) where u 0 (x, α) is the scattering solution. A similar formula was proved earlier in [6], p. 46, for the resolvent kernel of the Laplacian in the exterior of a bounded obstacle, (and even earlier, in [5], for some unbounded obstacles). The scattering amplitude for problem (.5) (.7) is A(β, α) = A 0 (β, α) + A M (β, α), (.3) where A 0 is defined in (.3) and A m is defined in (.6). If ka is sufficiently small, then k 2 is not a irichlet eigenvalue of the operator 2 q 0 (x) in m, m M. If Im ζ m 0, m M; Im q 0 (x) 0, (.4) then the unique solution to problem (.5) (.7) can be found in the form (.0). Theorem Assume (.8) and (.4). Then problem (.5) (.7) has a solution of the form (.0) and this is the unique solution of the problem (.5) (.7). Proof of Theorem is given in Section 3. In [3] there is a detailed study of boundary value problems in domains of the type := R 3 \ M m. In [3] the case of irichlet boundary condition on was studied, the case of Neumann boundary condition was mentioned as an open problem, and the case of impedance boundary condition was not studied. Let Note that We need two lemmas. Lemma If g(x, y) := eik x y 4π x y, g 0(x, y) := G(x, y) = g(x, y) 4π x y. (.5) g(x, z)q 0 (z)g(z, y)dz. (.6) t x a, x y d a, (.7) then ( a g(t, y) g(x, y) c d + ka ), (.8) 2 d where c > 0 stands for various positive constants independent of a and d. 4

5 Lemma 2 If (.7) holds, then ( a G(t, y) G(x, y) c d + ka ). (.9) 2 d These lemmas are proved in Section 3. Let us formulate our results under simplifying but physically reasonable assumptions. Theorem 2 Assume that lim x m m ζ m J m 4π = h(x m), where J m := dsdt s t, (.20) and for any subdomain the following relation holds = N(x)dx. (.2) a m Assume that = c a 2, and J m = c 2 a 3, where c, c 2 > 0 are constants independent of m. Finally we assume that M = O(a ) and d = O(a /3 ) as a 0. Under these assumptions there exists the limit: This u(x) solves the equations: u(x) = u 0 (x) and where the potential q is of the form: lim u M(x) = u(x) = u(x, α). (.22) M G(x, y) p(y) u(y)dy, Lu := [ 2 + k 2 q(x)]u = 0 in R 3, (.23) q(x) = q 0 (x) + p(x), and u satisfies the radiation condition: p(x) = 4πc2 N(x)h(x) c 2 [ + h(x)], (.24) ( u = e ikα x + A(β, α) eikr ) r + o, r = x, (.25) r where A(β, α) = A 0 (β, α) u 0 (y, β)p(y)u(y, α)dy, (.26) 4π and u 0 (y, β) is the scattering solution defined in (.3). 5

6 Theorem 3 Assume that ζ m = 0, m M, and the following limits exist: lim V m β (m) pj = β pj (y)ν(y)dy, (.27) m lim m V m = ν(y)dy, (.28) where ν(y) 0 and β pj (y) are continuous functions in, and β (m) pj is the magnetic polarizability tensor of the body m, defined in (2.38)-(2.39), see below. Then the function u M (x), defined in (.0), tends to the limit: and U(x) solves the equation: U(x) = u 0 (x) + G(x, y) U(y)ν(y)dy lim u M(x) = U(x) = U(x, α), (.29) M 3 p,j= G(x, y) y p U(y) y j β pj (y)ν(y)dy. If all the small particles are balls of radius a > 0, then V m = 4πa3 3, = 4πa 2, J m = 6π 2 a 3 dt, 4π s t = a, s. In this case dt s t = dtds 4π s t, (.30) that is, the mean value of the integral dt on the surface S s t m equals to this integral. If is not a sphere, this mean value is an approximate value of the above integral. Note that under the assumptions of Theorem 2 one has M = O(a ), while under the assumptions of Theorem 3 one has M = O(a 3 ) (see formula (2.50) below). Therefore, one needs many more particles to deal with the Neumann boundary condition, that is, with acoustically hard particles, than with the impedance boundary condition with large boundary impedance ζ = O(a ). We will discuss at the end of Section 4 in more detail the question concerning the compatibility of the assumption (.8), namely d a, and the existence of the limits (.27) and (.28). It will be shown that the assumption d a is compatible with the existence of the limit (.28) only if ν(y) is sufficiently small, and in this case the existence of the limit (.27) is also compatible with the assumption d a. In Section 2 Theorems 2 and 3 are proved. In Section 3 Theorem 3 and Lemmas, 2 are proved. In Section 4 some examples are given, the significance of the compatibility of the assumptions d a and (.2), (.27) (.28) is discussed, and a possible application of our results to creating materials with a desired refraction coefficient is described. 6

7 2 Proof of Theorem 2 Let us look for the solution to problem (.5) (.7) of the form: u M = u 0 (x) + G(x, s)σ m (s)ds, (2.) where G(x, y) is the resolvent kernel of L 0, see (.4), and σ m are arbitrary functions at the moment. For any σ m the function (2.) solves equation (.5) and satisfies the radiation condition (.6). Since problem (.5) (.7) has at most one solution, the function (2.) is the unique solution to (.5) (.7) provided that σ m are chosen so that the boundary conditions (.7) are satisfied. Since diam m, m M, are small, let us write (2.) as u M = u 0 (x) + G(x, x m )Q m + [G(x, s) G(x, x m )]σ m (s)ds, (2.2) where x m m is a point inside m and Q m := σ m (s)ds. (2.3) The choice of x m m is arbitrary because diam m 2a is small. We will prove that Q m 0, give an analytic formula for Q m (formula (2.20) below), and approximate the field u M in (2.2) by the expression: u M = u 0 (x) + G(x, x m )Q m. The error of this approximate formula is of order max( a, ka), see estimate (2.7) d below. Therefore this error tends to zero as a 0 since d = O(a /3 ). Let us estimate the term E m := [G(x, s) G(x, x m )]σ m (s)ds. (2.4) By the inequality (.9) one gets ( a E m c d + ka ) Q 2 m, x x m d a. (2.5) d We will prove below that Q m = O(a), see formula (2.20), and, since G(x, x m ) cd if x x m d > 0, one has: (a) G(x, x m )Q m = O. (2.6) d 7

8 Let us prove that under our assumptions the term E m is much smaller than O( a d ). Using again inequality (.9), one gets: Therefore, E m c(ad 2 + kad )O(a). E m O( a2 d + kaa 2 d ) O(a ), (2.7) d because ka and a d by assumption. So, our claim is verified. Moreover, E m G(x, x m )Q m if x x m d a, because M = O( ). a To find Q m, we use the boundary condition (.7). Let us write u(x) in a neighborhood of as u M (x) = u e (x) + G(x, s)σ j (s)ds, x x j 2a, (2.8) where u e is the effective field acting on the j th small particle from outside: u e (x) := u M (x) G(x, s)σ j (s)ds = G(x, x m )Q m + O( a ). (2.9) d m j We neglect the error term O( a ) in what follows. From (2.8) and (.7) one gets: d 0 = u en (s) ζ j u e (s) + A jσ j σ j 2 ζ j T j σ j, s, (2.0) where u en (s) is the normal derivative of u e at the point s. One can rewrite this equation as: σ j = A j σ j 2ζ j T j σ j 2ζ j u e (s) + 2u en (s). Here the operators A j and T j are defined as follows: σ j (t)dt T j σ j := G(s, t)σ j (t)dt 4π s t, (2.) G(s, t) A j σ j := 2 σ j (t)dt 2 N s N s 4π s t σ j(t)dt := Aσ j, (2.2) and we have used the following approximations: G(x, y) = g 0 (x, y)[ + O( x y )], x y 0; g 0 (x, y) := 4π x y, (2.3) G(x, y) y p = g 0 y p [ + O ( x y 2 ln x y )], x y 0. (2.4) 8

9 Note that (see [9], p. 96, formula (7.2)): Aσds = σ j ds (2.5) Indeed, ds The integral N S 2π s t σ(t)dt = dtσ j (t) ds N S 2π s t = σ j (t)dt. N S 2π s t ds =, t, is well known in potential theory for surfaces C,λ. Integrating (2.0) over, using formula (2.5), and the divergence theorem, one gets: Q j = ζ j u e (x j ) ζ j ds The function u e (y) is smooth, so σ j (t)dt 4π s t + u e dy. (2.6) j j u e (y)dy = V j u e (x j )[ + o()], a 0, (2.7) where V j = j is the volume of j and we have used the smallness of the diameter of j, that is, the smallness of a. Let us write ds = Q j σ j (t)dt 4π s t = dt σ j (t) S j dt ds 4π s t ds 4π s t = Q jj j 4π, J j := Here we approximated the continuous on function ds. s t dt If is a sphere of radius a, then ds s t s =a so in this case equation (2.8) is exact. From (2.6) (2.8) one finds a formula for Q j : ds dt s t. (2.8) ds by its mean value s t = 4πa, t = a, (2.9) Q j = ζ j + ζ jj j u e (x j ). (2.20) 4π 9

10 We neglected the term V j u e (x j ) = O(a 3 ) which is much smaller than ζ j = O(a) as a 0, because = O(a 2 ) and ζ j = O( ). The quantity J a j = O(a 3 ). Therefore ζ jj j 4π = O(). We choose ζ j = H(x j) a, (2.2) where H(x) is a continuous function in, which we can choose as we wish subject to the condition ImH 0, because Im ζ j 0. If the small particles are all of the same shape and size then = c a 2, and J j = c 2 a 3, where c, c 2 > 0 are constants independent of j, j M, and then and Formulas (2.2), (2.20), (2.22) and (2.23) imply: u M (x) = u 0 (x) ζ j J j 4π = H(x j)c 2 4πc := h(x j ), (2.22) ζ j = H(x j )c a. (2.23) G(x, x m ) 4πc2 c 2 h(x m)a u M (x m ), (2.24) + h(x m ) where x x m d a, and we replaced u e (x m ) by u M (x m ) because their difference (see (2.9)) is of order O( a ). Indeed d u M (x) u e (x) G(x, s) σ j (s) ds c d Q j c a d, x x j d a, where c, c > 0 are some constants independent of a. If the assumption (.2) holds, then lim G(x, x m ) 4πc2 c 2 h(x m ) u M (x m )a = + h(x m ) To pass to the limit in (2.26) one can use the following lemma. (2.25) G(x, y) 4πc2 c 2 h(y) u(y)n(y)dy. + h(y) (2.26) Lemma 3 Assume that x m m, diam m 2a, f is a continuous function in with a possible exception of a point y 0 in a neighborhood of which it is absolutely integrable, for example, it admits an estimate f(y) c y y 0, b < 3, and assume b that lim a = N(x)dx (2.27) m 0

11 for any subdomain, where N(x) is a continuous function. Then there exists the limit lim f(y m )a = f(y)n(y)dy. (2.28) Remark In our case f(y) = G(x, y) 4πc2 h(y) c 2 [+h(y)] u M(y) and (2.27) is the assumption (.2). Proof of Lemma 3 Let = P p= p, where p and q do not intersect each other, p is the closure of the domain p, and lim P max p diam p = 0. Choose any point y (p) p and note that Therefore lim sup f(y (p) ) f(y m ) < ε p 0 as diam p 0. (2.29) y m m, m p f(y m )a = lim = = P a P f(y m ) = p= m p p= P [f(y (p) ) + O(ε p )] N(y)dy p P [f(y (p) ) + O(ε p )] [N(y (p) ) + O(ε p)] p, (2.30) p= p= where lim P max p ε p = 0. Let P in (2.30). Then lim P P [f(y (p) ) + O(ε p )] [N(y (p) ) + O(εc p)] p = p= [f(y (p) ) + O(ε p )] lim a m p f(y)n(y)dy. (2.3) In the above argument we assumed that f is continuous in. If f has an integrable singularity at a point x 0, then we choose a ball B(x 0, δ ε ) centered at x 0 of radius δ ε such that sup 0<δ<δε B(x 0 f(y) dy < ε, where ε > 0 is an,δ) arbitrary small fixed number. Then f(y) N(y) dy < cε, sup 0<δ<δ ε B(x 0,δ) where c = max y N(y) > 0 is a constant independent of ε. Now we apply the above argument to the region \B(x 0, δ), where f is continuous and get: lim y m B(x 0,δ) f(y m )a = \B(x 0,δ) f(y)n(y)dy. (2.32)

12 The left side of (2.28) in the case of f having integrable singularity at the point x 0 and continuous in \x 0 is understood as the limit of the expression on the left of (2.32) as δ 0. This yields (2.28). Lemma 3 is proved. Passing to the limit M, or a 0, in equation (2.24) and using Lemma 3, one gets u(x) = u 0 (x) G(x, y) p(y) u(y)dy, p(x) = 4πc2 c 2 h(y)n(y) + h(y). Applying the operator L 0 = 2 + k 2 q 0 (x) to this equation and using (.3) and (.4), one obtains equation (.23). Formulas (.25) and (.26) follow from the above equation and from formula (.2). This concludes the proof of Theorem 2. Remark 2 It is possible (and not difficult) to generalize Theorem 2 to the case of particles with different shapes. Since this does not lead to an essentially new result, we do not go into detail. In [8] one can find analytical formulas for the S matrix for wave scattering by small bodies of arbitrary shapes. Proof of Theorem 3. Now we assume ζ m = 0, m M, which means that all the small particles are acoustically hard. In this case equation (2.0) takes the form σ j = A j σ j + 2u en (s), s, j M, (2.33) where u e (x) := u 0 (x) + m j G(x, s)σ m (s)ds. (2.34) We cannot use approximation (2.2) because the quantity Q m now is of the same order of magnitude as the integral [G(x, s) G(x, x m )] σ m (s)ds, or even smaller than this integral. This is established below. While under the assumptions of Theorem 2 we had Q m = O(a), now, under the assumptions of Theorem 3, we have Q m = O(ka 2 a 3 ), which is a much smaller quantity than O(a) because ka. To estimate the order of magnitude of Q m, we integrate (2.33) over and use (2.5). The result is: Q j = u en ds = u e dx u e (x j )V j, (2.35) j where V j is the volume of j, and we have used the assumption d a. This assumption allows one to claim that u e (x) is practically constant in the domain 2

13 j in the absence of j-th particle. ifferentiation with respect to x brings a factor k. Since we assume that k > 0 is fixed, this factor is not important for our argument, but to make the dimensionality of the term V j u e clear, we may write V j u e = O(k 2 a 3 ). This quantity has dimensionality of length since ka is dimensionless. We now prove that the term E m := [G(x, s) G(x, x m )]σ m (s)ds, which was neglected under the assumptions of Theorem 2, because it was much smaller than G(x, x m )Q m, is now, under the assumption ζ m = 0, m M, of the same order of magnitude as G(x, x m )Q m, namely O(k 2 a 3 d ), or even larger. We have [G(x, s) G(x, x m )]σ m (s)ds = y G(x, a) y=xm (s x m ) σ m (s)ds, x x m d a, (2.36) where we have used the assumption x x m a and kept the main term in the Taylor s expansion of the function G(x, s) G(x, x m ). Recall, that (s x m ) p σ m (s)ds = V m β (m) pj u e (y) y j y=xm, (2.37) where β (m) pj is the magnetic polarizability tensor defined in [9], (p.55, formulas (5.3)-(5.5) and p.62, formula (5.62)), and (s x m ) p is the p-th component of the vector s x m. Namely, if then S σ = Aσ 2N j, (2.38) s p σ(s)ds = V β pj, (2.39) where V is the volume of the body with boundary S, N j is the j-th component of the exterior unit normal N to S, the role of the point x m from equation (2.37) is played by the origin, which is located inside S, and the role of is played by S. Equation (2.33) with j = m can be written as ( σ m = A m σ m 2N j u e(x m ) ). (2.40) y j Compare (2.40) and (2.38) and get (2.37). Formulas for the tensor β pj = α pj (γ) for bodies of arbitrary shapes are γ= derived in [9], p.55, formula (5.5), so one may consider the tensor β pj known for bodies of arbitrary shapes. The parameter γ = ǫ i ǫ e ǫ i +ǫ e, where ǫ i is the dielectric 3

14 permittivity of the body and ǫ e is the dielectric permittivity of the surrounding medium. The case γ = occurs when ǫ i = 0. This is the case, for example, in the problem of calculation the magnetic dipole moment of a superconductor placed in a homogeneous magnetic field: in the superconductor the magnetic induction vector B = 0, which means that the magnetic permeability µ i of such body is zero, µ i = 0, see [2]. That is why the tensor β pj is called magnetic polarizability tensor in [9]. From (2.36) and (2.37) it follows that [ G(x, s) G(x, xm ) ] G(x, y) y=xm u e (y) y=xm σ m (s)ds = V m β (m) pj, y p y j (2.4) where one sums up over the repeated indices p, j, but nor over m. The quantity on the right of (2.4) is of the order O(k 2 a 3 d ) if kd, that is, of the same order as G(x, x m )Q m, provided that x x m d a, and it is of the order O(ka 3 d 2 ) if kd <. Indeed, β (m) pj = O(), V m = O(a 3 ), and y G(x, y) ( c max k, ) d d. 2 Let us prove the estimate y G(x, y) c max (k d, d 2 ) for x y d a, where c > 0 is a constant independent of d. We have G(x, y) = g(x, y) g(x, z)q 0 (z)g(z, y)dz := g TG, where T is compact as an operator in L p (), p under our assumptions, namely, R 3 is a bounded domain, q 0 (x) is a bounded piecewise-continuous function. From this equation we get so Clearly, y G(x, y) = y g(x, y) T y G. y g = g(ik x y ) y x x y, y g(x, y) 2 max ( k 2πd, 4πd 2 ) = 2π max( k d, d 2 ), x y d > 0. Thus y G(x, y) y g(x, y) [ +c x z x y 2 ] yg(z, y) dz := y g (+ci), ik x y 4

15 where I := One has and dz x y 2 y G(z, x) x z + k2 x y c 2 dz z y 2 x z c ln x y, sup ln x y x y 2 c, x,y where c = c() is a constant. Therefore c I + k2 d c, 2 dz x y 2 z y 2 x z + k2 d 2. and y G(x, y) c max ( k d, ) d 2 + k2 d c max( k 2 d, ), d 2 as claimed. If k, i.e. kd, then d d 2 y G(x, y) c k, x y d > 0. d If k <, i.e. kd <, then d d 2 y G(x, y) c, x y d > 0. d 2 Therefore, the right side of (2.4) is O ( ) k 2 a 3 d if kd, in which case it is of the same order as the term G(x, x m )Q m. If kd <, then the right side of (2.4) is O ( ) ka 3 d, in which case it may become larger than the term G(x, 2 xm )Q m because the ratio ka3 / k2 a 3 = > provided that kd <. d 2 d kd Writing the field (2.) in the form [ u M (x) = u 0 (x)+ G(x, x m )Q m + G(x, s) G(x, xm ) ] σ m (s)ds (2.42) and using formulas (2.35) and (2.4), one gets: u M (x) = u 0 (x)+ G(x, x m ) u e (x m ) V m G(x, x m ) u e (x m ) V m β (m) pj (x m ), y p y j (2.43) and over the repeated indices p, j one sums up. Let a 0, M. We want to give sufficient conditions for passing to this limit in (2.43). Lemma 4 Assume that for any subdomain the following limits exist: lim V m β (m) pj (x m ) = β pj (y) ν(y)dy, (2.44) m lim m V m = 5 ν(y)dy. (2.45)

16 Then the limiting form of equation (2.43) is: U(x) = u 0 (x) + G(x, y) U(y) ν(y)dy G(x, y) U(y) β pj (y) ν(y)dy, (2.46) y p y j where one sums up over the repeated indices p, j. Remark 3 If one assumes that ν(y) vanishes near the boundary S of and integrates the last integral in (2.46) by parts, one gets U(x) = u 0 (x) + { G(x, y) U(y)ν(y) + 3 p,j= ( U(y) )} β pj (y) ν(y) dy. y p y j (2.47) Applying the operator L 0 to both sides of (2.47) and using (.4) one gets: L 0 U + ν(y) U(x) + 3 p,j= ( U(y) ) β pj (y) ν(y) = 0. (2.47 ) y p y j Remark 4 If all the small particles are identical, then V m = c 3 a 3, where the positive constant c 3 does not depend on m, and β (m) pj = β pj. Then [ lim V m = lim c a ] 3 = lim[c a 3 N( )], (2.48) m m where N( ) is the number of small particles in the domain. For the limit (2.48) to exist it is sufficient that N( ) = ν(y)dy, (2.49) c a 3 where ν(y) 0 is a continuous function, and the limit in (2.48) is equal to ν(y)dy. One can write (2.49) as N(y)dy = ν(y) dy. (2.50) c a3 In contrast to Theorem 2, where M = O( a ), we now have M = O( a 3 ). Similarly, lim m V m β (m) pj (x j ) = lim[c a 3 β pj N( )] = β pj ν(y)dy. (2.5) 6

17 We gave in this Remark some practically realizable sufficient conditions for the existence of the limits (2.44) and (2.45). Let us verify that if the limits (2.44) (2.45) exist, then the limit of the right side of equation (2.43) exists, and, denoting this limit by U(x) = lim u M (x), one obtains the limiting form of equation (2.43): U(x) = u 0 (x) + G(x, y) U(y) ν(y)dy which is equation (2.46). We took into account that 3 p,j= G(x, y) y p lim u e(x) = U(x). U(y) y j β pj (y) ν(y)dy, (2.52) This is so because, as a 0, the input of a single particle into the field U(x) tends to zero. Let us verify the existence of the limit of the right side of equation (2.43). We use, as in the proof of Lemma 3, a representation of of the form = P p= p, and assume that Then lim lim max diam p = 0. (2.53) P p P G(x, y m ) u e (x m )V m = = = P lim G(x, x m ) u e (x m )V m m p P G(x, y (p) ) u e (y (p) )( + ε p ) lim p= p= V m m p P G(x, y (p) ) u e (y (p) )( + ε p ) ν 2 (y (p) )( + ε p ) p, (2.54) p= where lim max ( ε p + ε p ) = 0. (2.55) P p 7

18 Let P in (2.54) and use (2.55) to get lim P P G(x, y (p) ) u e (y (p) ) ν(y (p) ) p ( + ε p + ε p + ε p ε p) p= = G(x, y) U(y) ν(y)dy. (2.56) We have replaced u e (y (p) ) in the limit P by U(y), because U(y) u e (y) = G(s, t) σ(t)dt = o() as a 0. (2.57) From (2.56) and (2.54) one gets: lim G(x, x m ) u e (x m )V m = G(x, y) U(y) ν(y)dy. (2.58) The singular points x = y of G(x, y) are treated as in the proof of Theorem 2. The function G(x, y) c x y as x y 0, so G(x, y) L () as a function of y for any x. Similar arguments, applied to the last sum in (2.43), lead to the formula lim G(x, x m ) u e (x m ) V m β (m) G(x, y) pj (x m ) = y p y j y p where one sums up over the repeated indices p, j. Theorem 3 is proved. U y j β pj (y) ν(y)dy, (2.59) In Section 4 we discuss the compatibility of the condition d a and the existence of the limit (.28). 3 Auxiliary results In this Section we prove Theorem and Lemmas, 2. Proof of Lemma Let us start with the following observations: t y x y t y (x y) = t x a, (3.) sup e is a, (3.2) a s a e ik t y e ik x y = e ik( t y x y ) ka, (3.3) 8

19 where the last inequality follows from (3.2). One has g(t, y) g(x, y) = x y e ik t y t y e ik x y x y t y x y t y x y t y + + t y e ik t y e ik x y (3.4) x y t y a x y t y + ka x y a d 2 ( a) + ka d O( a ) ka + d d 2 d. Lemma is proved. Proof of Lemma 2 Let us start with the equation: G(x, y) = g(x, y) g(x, z)q 0 (z)g(z, y)dz, (3.5) where q 0 is defined in (.2). From (3.5) one gets: G(t, y) G(x, y) g(t, y) g(x, y) + [g(t, z) g(x, z)]q 0 (z)g(z, y)dz O ( a ka d2) + d + c dz g(t, z) g(x, z) z y. (3.6) Here we have used Lemma and the estimates sup q 0 (z) c 4, G(z, y) c 5 z y, (3.7) z where c 4, c 5 > 0 are some constants. Let us estimate the integral dz I := g(t, z) g(x, z) z y g(t, z) g(x, z) dz g(t, z) g(x, z) dz = + x z d 4,z z y x z d 4,z z y := I + I 2. (3.8) By Lemma, which is applied to I with d replaced by d, one gets 4 I c ( a d + ka ) 2 d d x z 4 dz z y c ( a d + ka ). (3.9) 2 d Here and below we do not write z under the integration sign to simplify the notations. 9

20 Let us estimate I 2 : I 2 dz e ik t z x z e ik x z t z. (3.0) 4π x z d z y t z x z 4 One has e ik t z x z t z e ik x z x z t z + t z e ik t z e ik x z x t + t z k t z x z (3.) x t + k t z t x. Thus, with x y d a and t x a, one has: I 2 dz( x t + k t z t x ) 4π z y t z x z t x 4π x z d 4 ( x z d 4 dz z y t z x z + k x z d 4 dz ) z y x z ca ( d + k ). (3.2) 2 d From (3.9) and (3.2) the estimate (.9) follows. Lemma 2 is proved. Proof of Theorem Let us first prove that if conditions (.4) hold, then problem (.5) (.7) has at most one solution. It is sufficient to prove that the homogeneous problem ( 2 + k 2 q 0 )u = 0 in R 3 \ M m, (3.3) u r iku = o( r ), u = O( ), r = x, (3.4) r u N = ζ m u on, m M, (3.5) has only the trivial solution if conditions (.4) hold. Taking complex conjugate of (3.3) (3.5) one gets: ( 2 + k 2 q 0 (x) ) u = 0 in R 3 \ M m, (3.6) u r + iku = o( ) (), u = O, r = x, (3.7) r r u N = ζ m u on, m M. (3.8) Multiply (3.3) by u, (3.6) by u, subtract from the first equation the second one, and integrate over the region (R 3 \ M m) B R := R, where B R is the 20

21 ball centered at the origin of radius R. Using Green s formula, one gets: 0 = [u 2 u u 2 u (q 0 q 0 ) u 2 ]dx R ( = 2i Im q 0 (x) u 2 u dx + u R x =R r u u ) ds r ( u u N u u ) ds. (3.9) N Using (3.7) and (3.8) one rewrites (3.9) as follows: 0 = 2i Im q 0 (x) u 2 dx + 2ik u 2 ds + o() 2i Im ξ m u 2 ds. R x =R (3.20) Letting R, taking into account that q 0 (x) = 0 in = R 3 \, and one gets: 0 \ S M m Im q 0 (x) u 2 dx + Im ζ m u 2 ds k lim sup R x =R u 2 ds. (3.2) Since all the terms on the right side of this relation are non-positive by the assumptions (.4), it follows that lim sup R x =R u 2 ds = 0. This implies that u = 0 (see, [8], p. 23). Thus, uniqueness of the solution to problem (.5).7 is proved. Let us prove the existence of the solution to (.5) (.7) of the form (.0). The existence of the solution of the form (.0) will be established if one proves the existence of σ m, m M, such that boundary condition (.7) is satisfied: u en ζ j u e + A jσ j σ j 2 ζ j T j σ j = 0, j M. (3.22) Here u e, which depends on j, is defined by the formula: u e := u G(x, s)σ j (s)ds = u 0 + G(x, s)σ m (s)ds. (3.23) m j Under our assumptions C,λ uniformly with respect to m. Therefore equation (3.22) is of Fredholm type in the space L 2 ( M ). The corresponding homogeneous equation, i.e., the equation with u 0 = 0, cannot have a nontrivial solution because such a solution would generate by formula (.0) with u 0 = 0 a 2

22 function u M (x) = M G(x, s)σ m (s)ds, which would solve the homogeneous problem (.5) (.7). We have already proved that such a function has to be zero in R 3 \ M m. Thus, u M = 0, m M, and u M solves the problem: Sm L 0 u M = 0 in m, u M Sm = 0, m M. (3.24) If diam m 2a is sufficiently small, then problem (3.24) has only the trivial solution for every m, m M. Therefore u M = 0 in m and in R 3 \ M m. Therefore, by the formula for the jumps of the normal derivatives of the single layer potential, u + M N u M N = σ m, we conclude that σ m = 0, m M. This implies the existence of the solution to problem (.5) (.7) of the form (.0). Theorem is proved. Let us return to the assumptions of Theorem 2, namely, ζ m = O( a ), an( b(y)) = N(y) b (y) ( + o()), where b (y) is a cube, centered at the point y with the side b > 0, and o() is related to the limiting process b 0. Under these assumptions let us establish an estimate for the function v M := u M u 0, which is uniform with respect to M, or a 0. From this estimate it follows that v M converges, as a 0, in L 2 (R 3, ( + x ) γ ), where γ > 0 is an arbitrary fixed constant. The function v M satisfies the radiation condition at infinity. The function u 0 Hloc 2 (R3 ) solves the equation L 0 u 0 = 0 in R 3. Let e := R 3 \ M m and S := M. Let ( ) /2, M ( ) /2. v := v(x) 2 ( + x ) γ dx v = ( v N 2 + v 2 )ds e The estimate we wish to prove is: v M c u 0. (3.25) Here and below c > 0 stand for various constants independent of a. Let us outline the proof of inequality (3.25). Step. If M = O( ), then the right side of (3.25) is bounded as a 0. a Indeed, the number of small particles is M = O( ) and u a 0 is Hloc 2 (R3 ), so that u 0 and u 0N are bounded in L 2 ( ) uniformly with respect to m, m M. Thus, u 0 c ( ) /2 a max ( u 0N 2 + u 0 2 c )ds m M a /2 c, 22

23 where c > 0 stand for various constants independent of a. Step 2. If the inequality (3.25) is false, then there is a sequence u (n) 0, u(n) 0 =, such that v (n) M := v(n) n. efine w (n) := v(n). Then v (n) w (n) =. (3.26) From the weak compactness of bounded sets in L 2, it follows, that one may select a subsequence, denoted again w (n), such that w (n) converges weakly in L 2 loc ( ) to a function w. The function w (n) solves the problem: w (n) N L 0 w (n) = 0 in e, ζ mw (n) = (ζ m u (n) 0 u (n) 0N )/ v(n) on, m M, (3.27) and w (n) satisfies the radiation condition. It follows from (3.27) that 2 w (n) < c, so w (n) H 2 loc ( e) < c, where Hloc 2 ( e) is the Sobolev space. Thus, one may assume, using the compactness of the embedding from Hloc 2 into L2 loc, that w(n) converges to w strongly in L 2 loc ( e). This and equation (3.27) imply that w (n) converges to w strongly in Hloc 2 ( e), so that w solves equation (3.27), satisfies the radiation condition and the homogeneous boundary condition (3.27), that is, w N ζ m w = 0 on. Therefore, by already proved uniqueness theorem (see the proof of Theorem ), we conclude that w = 0. The terms u (n) 0 / v (n) and u (n) 0N / v(n) tend to zero as n, because v (n) > n. Therefore, the limiting function w satisfies the homogeneous boundary condition w N = ζ m w on, m M. Let us prove that w (n) (x) < c, x > R, where R > 0 is sufficiently large x and c > 0 does not depend on n. For w (n) one has a representation by the Green formula in the region x > R, where R > 0 is large enough, so that the ball B R := {x : x < R} contains. Namely w (n) (x) = (w r (n) g(x, s) g r (x, s)w (n) )ds, x > R, (3.28) s =R where the derivatives with respect to r are the derivatives along the normal to the sphere S R := {s : s = R}, and g is defined in (.5). It follows from (3.28) that w (n) (x) < c for x > R, where c > 0 is a constant independent of n, x because local converegence in H 2 implies that the L 2 (S R )-norms of w (n) and of w r (n) are bounded uniformly with respect to n. Therefore lim n w(n) w = 0, (3.29) 23

24 because on compact sets lim n w (n) w H 2 loc ( e) = 0, and near infinity the inequality w (n) (x) < c implies that x w (n) (x) 2 ( + x ) γ = O(R γ ) 0, R, {x: x >R} so that (3.29) holds. Because of the uniqueness of the limit, not only a subsequence of w (n) but the sequence itself converges to w as n. This leads to a contradiction, because w = 0 and (3.26) together with (3.29) imply w =. This contradiction proves inequality (3.25). From inequality (3.25) and Step one concludes that that u M contains a weakly convergent in L 2 loc ( e) subsequence. By the arguments, similar to the given above, this subsequence converges in L 2 (R 3, ( + x ) γ ). Its limit solves equation (.23). The relation M = O( ) plays an important role in our proof of Theorem 2 a and in Step in the above argument. 4 Application to creating smart materials Let us ask the following question: can one make a material with a desired refraction coefficient n(x) in a bounded domain R 3, filled by a material with a known refraction coefficient n 0 (x), for example n 0 (x) = n 0 = const in, by embedding into a number of small particles, each of which is defined by its shape and its boundary impedance? Consider first the particles satisfying the assumptions of Theorem 2. More specifically, suppose that all the particles are balls of the same radius a. In this case so and formula (.24) yields = 4πa 2, J m = s =a c = 4π, c 2 = 6π 2, t =a ds dt s t = 6π2 a 3, 4πc 2 c 2 = 4π, p(x) = 4πN(x) h(x) + h(x), (4.) where h(x) is defined by the choice of the boundary impedances by formula (.20): and N(x) is defined by formula (.2). ζ(x) = h(x) a, (4.2) 24

25 If the original refraction coefficient is n 0 (x), then the corresponding potential is q 0 (x) = k 2 [ n 0 (x)] by formula (.2). If the desired refraction coefficient in is n(x), then the corresponding potential is q(x) = k 2 [ n(x)], so p(x) = q(x) q 0 (x) = k 2 [n 0 (x) n(x)]. (4.3) To create a material with the desired refraction coefficient n(x) it is sufficient to choose N(x) and h(x) so that (4.) holds with p(x) defined in (4.3). If the new material with the refraction coefficient n(x) has some absorption, that is, Im n(x) 0, and Im n 0 = 0, then Im p(x) 0. Let us prove that any function p(x) in with Im p 0, can be obtained (in many ways, non-uniquely) by formula (4.) with some choices of a nonnegative function N(x) and a function h(x) with Im h 0. Let p(x) = p (x) + ip 2 (x), p 2 (x) 0, and h(x) = h (x) + ih 2 (x), h 2 (x) 0. Assume that p(x) is given. Then (4.) implies Thus p + ip 2 = 4π (h + ih 2 )( + h ih 2 ) ( + h ) 2 + h 2 2 N(x). (4.4) p = 4π h + h 2 + h2 2 N(x), p ( + h ) 2 + h 2 2 = 4π N(x). (4.5) 2 ( + h ) 2 + h 2 2 There are many choices of the three functions: N(x) 0, h 2 (x) 0 and a realvalued function h (x) such that relations (4.5) hold. For example, if p > 0 and p 2 0, then one can choose h 2 h (x) = 0, h 2 (x) = p (x) p 2 (x), N(x) = p2 (x) + p2 2 (x) 4π p (x). (4.6) It is a simple matter to check that relations (4.5) hold with the choice (4.6). Since one has three functions h (x), h 2 (x) 0 and N(x) 0 to satisfy two equations (4.5) with p 2 (x) 0, there are many ways to do this. A particular choice of h(x) = h (x) + ih 2 (x) and N(x) 0 yields the surface impedance ζ(x) of the particles to be embedded around each point x, ζ(x) = h(x) by a formula (4.2), and the number of particles per unit volume around the point x, namely, by formula (.2) this number is N(x), so that the number of particles to a be embedded in the volume dx around point x is equal to N(x) dx. The smallest a distance d between the embedded particles should satisfy the inequality d a. One may try to take practically d > 0a. Example Suppose that the elementary subdomain p, used in the proof of Lemma 3, is a cube with the side b d, x p. Let, for example b = 0 2 cm, d = 0 3 cm, a = 0 5 cm. Then there are ( b 3 d) = 0 3 small particles in p around a point x, the center of p. The function N(x) in p in this example 25

26 is found from the formula N(x) b 3 = 0 3 (use (.2) with = a p ), so N(x) = = 0 4. The number of small particles, embedded in the cube p around point x, the center of this cube, is 0 3 in this example. The relative volume of these particles in p is π = , so it is quite small, which is in full agreement with our theory. The assumption (.8), specifically, d a, is compatible with the requirement (.2). Indeed, if one denotes by N( ) the left side of (.2), then N( ) = O ( ) a for any. Let us assume that is a unit cube, and denote by N( ) the left side of (.2). The assumption d a implies that the number N( ) of particles in is O ( ) ( d. These relations are compatible if and only if O ( a) = O ) 3 d, i.e., 3 d = O(a /3 ). Therefore, it is possible to have a 0, a 0 and equation (.2) d satisfied. Let us discuss the new material properties, specifically, anisotropy, when acoustically hard particles are embedded in the domain, and the assumptions of Theorem 3 are valid. The physical situation is now quite different from the one in Theorem 2. From the physical point of view one can anticipate the drastic difference because the wave scattering by one small acoustically soft particle of the characteristic size a is isotropic and the scattering amplitude is of order a, while the wave scattering by a small acoustically hard particle is anisotropic and the corresponding scattering amplitude is of order k 2 a 3, (see [9], chapter 7). We assume that ka, say ka < 0., so that the quantity k 2 a 3 = (ka) 2 a is 00 times less than a. Example 2 Let us assume again that the small particles are all balls of the same radius a. Then V m = 4 3 πa3, ν(y) p = 4 3 πa3 N( p ), where N( p ) is the number of small particles in a small cube p centered at the point y. If b is the size of the edge of the cube p, then ν(y) = 4.8 a3 N( b 3 p ), where 4.8 is an approximate value of 4π. The magnetic polarizability tensor β 3 pj of a ball of radius a is β pj = 3 δ 2 pj, while the electric polarizability tensor of a perfectly conducting ball is 3δ ij, where {, p = j, δ pj = 0, p j. These values differ by the factor 4π from the values in [2] because we use the formula ϕ = for the potential of a point charge, while in [2] this potential is 4π x. In our example β x pj does not depend on m. Therefore the limit (.27) exists if the limit (.28) exists. The limit (.28) exists if and only if the following limit exists: 4π 3 lim a3 = ν(y)dy, (4.7) m 26

27 where ν(y) is the function defined in (.28). Thus, in contrast to Example, where N( p ) = O ( ) a, now we have N( p ) = O ( ) a. The relative volume of the 3 small particles in Example 2 is not negligible and does not go to zero as a 0, in contrast to Example. Let us discuss the compatibility of the condition d a and the existence of the limits (.27) and (.28). If the condition d a is compatible with the existence of the limit (.28), then it is compatible with the existence of the limit (.27). If the limit (.28) exists, then a 3 N( ) = O(), so N( ) = O(a 3 ). On the other hand, N( ) = O(d 3 ). These relations, in general, are not compatible because d a. Let us argue more precisely. Let = p, where p is a cube with the edge of size b. Let us assume that the small particles in p are identical and their characteristic size is a. If (.28) holds, where ν(y) is continuous, and if b is small, then the right side of (.28) equals to ν(y)b 3, y p =. The left side of (.28) equals to c 3 a 3 N( p ). Thus N( p ) = c 3 ν(y) b3. On the other hand, a 3 N( p ) = b3, prvided that one assumes that the centers of the small particles d 3 are at the uniform grid, so that there are b centers on the segment of length b. d If c 3 ν(y) b3 = b3, then a = ( ) ν(y) /3 a 3 d 3 d c 3. Therefore the condition d a is satisfied only if ( ν(y) ) /3 c 3, say ( ν(y) ) /3 c The number c 3 depends on the shape of the particle. If the particles are balls of radius a, then c 3 = 4.8. Therefore ν(y) The conclusion is: The condition d a is compatible with the existence of the limit (.28) only if the function ν(y) in (.28) is sufficiently small. In general, equation (.30) cannot be reduced to a local differential equation for U(x). However, if ν(y) is small, one may use perturbation theory to study equation (.30). However, under an additional assumption, reasonable from the physical point of view, one can reduce integral-differential equation (.30) to a differential equation. Namely, let us assume that ν(y) is a continuously differentiable function in which vanishes near the boundary S. Under this assumption one can integrate by parts the last integral in (.30) and get: U(x) = u 0 (x) + [ G(x, y) U(y)ν(y) + 3 p,j= ( U(y) )] β pj (y)ν(y). (4.8) y p y j Let us apply the operator L 0 = 2 + k 2 q 0 (x) to (4.8) and use (.4) to get: [ 2 + k 2 q 0 (y)]u + ν(y) 2 U(x) + 3 p,j= ( U(x) ) β pj (y) ν(y) = 0, (4.9) y p y 27

28 where U(x) satisfies the radiation condition of the type (.6). This is an elliptic equation and the perturbation P of the operator L 0 is: PU := ν(x) 2 U(x) + 3 p,j= ( U ) β pj (y) ν(x). (4.0) y p y j This perturbation is the sum of the terms with positive small coefficient ν(y) in front of the second derivatives of U and a term with the first order derivatives of U: [ PU = ν(x) 2 U(x) + 3 p,j= ( U(x) )] β pj (x) + x p x j 3 p,j= U(x) x j β pj (x) ν(x) x p. (4.) If both ν(x) and ν(x) are small, this equation can be studied by perturbation methods. The physical effect on the properties of the new material, created by embedding into small acoustically hard particles, consists in appearing of anisotropy in the new material: the propagation of waves is described by the integral-differential equation (4.8) or (under the additional assumption on ν(y), namely: ν(y) vanishes near the boundary S of ) by the differential equation (4.9) with variable coefficients in front of the senior (second order) derivatives and the terms with the first order derivatives. The role of the compatibility of the assumption d a and of the assumption (.28) is quite important. Although passing to the limit a 0, justified in the proof of Theorem 3, is based on the assumptions (.27) and (.28), but without the assumption d a one cannot expect, in general, that the effective field u e (x), acting on any single particle, is practically constant on the distances of the order 2a. This physical assumption is important for our theory. From the mathematical point of view, if ν(x) is not sufficiently small, then the existence of the unique solution to equation (4.8) or of the solution to equation (4.9), satisfying the radiation condition, is not guaranteed. If, on the other hand, the quantity ( ) sup ν(x) + ν(x), x R 3 that is, this quantity is sufficiently small, then one can argue that the norm of the integral operator in (4.8) in L 2 () is small, so that equation (4.8) has a unique solution in L 2 (). This solution admits a natural extension to the whole space R 3 by the right side of (4.8) because ν(y) vanishes outside. Since G(x, y) satisfies the radiation condition, the solution to (4.8) also satisfies this condition. Without the assumption that ν(x) + ν(x) is sufficiently small, one cannot use the above argument. With this assumption one may solve equation (4.8) by iterations and find in this way an approximate solution to this equation. The first iteration yields the 28

29 following approximate solution to equation (4.8): U(x) = u 0 (x) + [ G(x, y) u 0 (y)ν(y) + 3 p,j= ( u0 (y) )] β pj (y) ν(y) dy. y p y j (4.2) Formula (4.2) gives the correction to the solution u 0 (x) of the unperturbed scattering problem, i.e., the scattering problem in the absence of small bodies. Since one has u 0 = k 2 n 0 (x)u 0, (4.2) can be rewritten as: U(x) = u 0 (x) k 2 + G(x, y)n 0 (y)u 0 (y)ν(y)dy G(x, y) 3 p,j= ( u0 (y) ) β pj (y)ν(y) dy. (4.3) y p y j In [3], Chapter 3, Section 3, the Neumann problem for the Helmholtz equation with n 0 (x) = was studied in the domain, similar to the one in equation (.5) and it was proved under the assumptions used in [3], that the main term of the asymptotics of the solution, as the relative volume of the particles tends to zero, is the incident field, while the next term is proportional to this relative volume. References [] assios, G., Kleinman, R., Low frequency scattering, Oxford Univ. Press, Oxford, [2] Landau, L.., Lifshitz, E.M., Electrodynamics of continuous media, Pergamon Press, Oxford, 960. [3] Marchenko, V., Khruslov, E., Boundary-value problems in domains with fine-grained boundary, Naukova umka, Kiev, 974. (in Russian). [4] Milton, G., The theory of composites, Cambr. Univ. ress, Cambridge, [5] Ramm, A.G., Investigation of the scattering problem in some domains with infinite boundaries I, II, Vestnik 7, (963), 45-66; 9, (963), [6] Ramm, A.G., Scattering by obstacles,.reidel, ordrecht, 986. [7] Ramm, A.G., Equations for the self-consistent field in random medium, Phys.Lett. A, 32, N3-4, (2003), [8] Ramm, A.G., Inverse problems, Springer, New York,

30 [9] Ramm, A.G., Wave scattering by small bodies of arbitrary shapes, World Sci. Publ., Singapore, [0] Ramm, A.G., Electromagnetic wave scattering by many small particles, Phys. Lett. A, 360, N6, (2007), [] Ramm, A.G., istribution of particles which produces a desired radiation pattern, Communic. in Nonlinear Sci. and Numer. Simulation, 2, N7, (2007), 5-9. [2] Ramm, A.G., istribution of particles which produces a smart material, Jour. Stat. Phys., 27, N5, (2007), [3] Ramm, A.G., istribution of particles which produces a desired radiation pattern, Physica B, 394, N2, (2007), [4] Ramm, A.G., Many-body wave scattering by small bodies, J. Math. Phys., 48, N2, 02352, (2007). [5] Ramm, A.G., Wave scattering by small particles in a medium, Phys. Lett. A 367, (2007), [6] Ramm, A.G., Wave scattering by small impedance particles in a medium, Phys. Lett A, (2007) (to appear) [7] Ramm, A.G., Materials with a desired refraction coefficient can be made by embedding small particles, Phys. Lett A, (2007) (to appear) [8] Shivola, A., Electromagnetic mixing formulas and applications, IEE,