arxiv: v1 [math.ap] 1 Jun 2018

Transcription

1 REFRACTORS IN ANISOTROPIC MEDIA ASSOCIATED WITH NORMS arxiv: v [math.ap] Jun 208 CRISTIAN E. GUTIÉRREZ, QINGBO HUANG AND HENOK MAWI Abstract. We show existence of interfaces between two anisotropic materials so that light is refracted in accordance with a given pattern of energy. To do this we formulate a vector Snell law for anisotropic media when the wave fronts are given by norms for which the corresponding unit spheres are strictly convex. Contents. Introduction 2 2. Preliminaries on Norms and Convexity 4 3. A vector Snell s law for anisotropic media Physical constraints for refraction 8 4. Uniformly refracting surfaces 0 5. The refractor problem when κ <, κ in (4.2) 6. Propagation of light in anisotropic materials Case considered for the application of our results Connection with optimal mass transport 25 References 25 June 4, 208. The first author was partially supported by NSF grant DMS , and the third author partially supported by NSF grant HRD

2 2 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI. Introduction Anisotropic materials are those whose optical properties vary according to the direction of propagation of light. Typical examples are crystals, where the refractive index depends on the direction of the incident light, see [BW59, Chapter XV], [LL84, Chapter XI] and [Som54, Chapter IV]. Important research was done on this subject because of it multiple applications, see the fundamental work [KK65], and [YY84], [Sch07] for more recent applications and references. Mathematically, in these materials wave fronts satisfy the Fresnel partial differential equation which in the particular case of isotropic materials is the eikonal equation. A difficulty with anisotropic materials in the geometrical optics regime is that incident rays may be refracted into two rays: an ordinary ray and an extraordinary ray. This is the phenomenon of bi-refringence, observed experimentally in crystals, and is a consequence from the fact that the Fresnel equation splits in general as the product of two surfaces, see (6.8). The main purpose of this paper is to show existence of interfaces between two homogenous and anisotropic materials so that light is refracted in accordance with a given pattern of energy. As a main step to achieve this, we give a formulation of a vector Snell s law in anisotropic materials when the wave fronts are given by norms in R n which has independent interest. More precisely, suppose N i (x), i =, 2, are norms in R n, Σ i = {x : N i (x) = }, Ω i Σ i are domains, f > 0 is an integrable function on Ω, and µ is a Radon measure in Ω 2 with f (x) dx = µ(ω Ω 2 ). We have two anisotropic media I and II such that the wave fronts in I are given by N and the wave fronts in II given by N 2. Light rays emanate from the origin, located in medium I, with intensity f (x) for each x Ω. We seek a surface S separating media I and II so that all rays emanating from the origin and with directions in Ω are refracted into rays with directions in Ω 2 and the conservation There is no mathematical objection to look at this problem in n dimensions but the physical problem is three dimensional.

3 REFRACTORS IN ANISOTROPIC MEDIA 3 of energy condition τ(e) f (x) dx = µ(e) holds for each Borel set E Ω 2 where τ(e) = {x Ω : x is refracted into E}. This is called the refractor problem, and when media I and II are homogeneous and isotropic it is solved in [GH09] using optimal mass transport and in [Gut4] with a different method. A main difficulty to solve this problem is to lay down the mathematical formulation of the physical laws and constraints in anisotropic media to our setting with norms. To place our results in perspective we mention the following. A similar problem for norms but for reflection was studied in [CH09]. Once the Snell law and physical constraints for anisotropic media are formulated and proved in Section 3, our existence results use the abstract method developed in [GH4], where existence results for the near field refractor problem in homogenous and isotropic media are obtained. Further results on geometric optics problems for refraction in homogeneous and isotropic media have been studied in [GM3], [LGM7], and [Kar6]. On the other hand, the mathematical literature for these problems in anisotropic media is lacking. The paper is organized as follows. Section 2 recalls a few results on norms and convexity that will be used later. The Snell law in anisotropic media is proved in Section 3 as a consequence of Fermat s principle of least time. The discussion on the physical constraints for refraction in anisotropic media is in Section 3.. Section 4 introduces and analyzes the surfaces refracting all rays into a fixed direction which are used in Section 5 to show the existence Theorems 5.6 and 5.7. Section 6 introduces and analyzes the Fresnel pde for the wave fronts in general materials non homogenous and anisotropic. In Section 6. we apply the results from the previous sections to materials having permittivity and permeability coefficients ɛ and µ that are constant matrices with one a constant multiple of the other. Finally, in Section 7 we relate our problem with optimal mass transport.

4 4 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI 2. Preliminaries on Norms and Convexity Consider a norm N(x) = x in R n, and let Σ = {x R n : N(x) = } be the unit sphere in the norm such that Σ is a strictly convex surface. Given a vector ν S n, the support function of Σ is defined by ϕ(ν) = sup x ν. x Σ Clearly ϕ is strictly positive (ϕ(ν) / ν ). Since Σ is compact, there is x 0 Σ such that ϕ(ν) = x 0 ν, and since Σ is strictly convex x 0 is unique. The hyperplane Π ν with equation {x : l(x) = x ν ϕ(ν) = 0} is a supporting hyperplane to Σ at x 0. That is, l(x) 0 for all x Σ and l(x 0 ) = 0. We then have a map ν S n x 0 Σ Φ : S n Σ, Φ(ν) = x 0. This map assigns to each vector ν S n a unique point Φ(ν) Σ (uniqueness follows from the strict convexity of Σ) so that Π ν is a supporting hyperplane to Σ at Φ(ν). In addition, also from the strict convexity of Σ, and consequently the uniqueness of maximizer of ϕ, we have Φ( ν) = Φ(ν). Φ is called the support map. The dual norm of N, denoted by N (x) = x is defined as follows. Since R n is finite dimensional, each linear functional l : R n R can be represented as l(x) = y x for a unique y R n. Given y R n, the dual norm of N is then N (y) = sup x Σ x y, and writing y = λ ν with ν S n, we get that N (y) = λ ϕ(ν). Hence the dual norm sphere of Σ is Σ = {y R n : N (y) = } = {ν/ϕ(ν) : ν S n }. We recall the following. Lemma 2.. [CH09, Lemma 2.3] For each x Σ and ν = ν/ϕ(ν) with ν S n we have (a) x ν ; and (b) x ν = if and only if x Φ(ν). Since Σ is a convex surface, for each x Σ, there is a supporting hyperplane to Σ at x and let ν(x) be the outer unit normal to such a supporting hyperplane. If Σ

5 REFRACTORS IN ANISOTROPIC MEDIA 5 is such that at each x one can pick a supporting hyperplane with normal ν(x) in such a way that ν(x) is continuous for all x Σ, that is, Σ has a continuous normal field, then from the proof of [CGH08, Theorem 6.2], Σ has a unique tangent plane at each point; that is, Σ is differentiable. With the notation from [CH09], the Minkowski functional of Σ defined by M Σ (x) = inf{r > 0 : x r B} satisfies M Σ (x) = N(x), where B = {x R n : N(x) }. It is proved in [CH09, Lemma 2.4] that Σ has a continuous normal field ν(x) if and only if N C (R n \ {0}); and p(x) := N(x) = ν(x)/ϕ(ν(x)). We also recall [CH09, Lemma 2.5] saying that if Σ is C, then Σ is strictly convex. Also, if Σ is strictly convex, then Σ C, [CH09, Lemma 2.7]. Therefore, if Σ is strictly convex and C, then p : Σ Σ is an homeomorphism and p p = Id; p = N. 3. A vector Snell s law for anisotropic media We have two homogenous 2 and anisotropic media I and II so that the surfaces for the wave fronts are given by a norm N in I, and given by a norm N 2 in II. 3 We are assuming that the norms N i C (R n \ {0}) and the corresponding unit spheres Σ i are strictly convex, i =, 2. Suppose I and II are separated by a plane having normal ν from medium I to medium II as in Figure. We formulate the Snell law in anisotropic media as follows: Each incident ray traveling in medium I with direction x Σ with x ν 0 and striking the plane Π at some point P 0 is refracted in medium II into a direction m Σ 2 if (3.) p 2 (m) p (x) ν, with p i = N i, i =, 2; see Figure. For each x Σ with x ν 0, there is at most one m Σ 2 satisfying (3.) with m ν 0. In fact, from (3.) there is λ R with 2 In homogenous media light travels in straight lines. This follows from [KK65, Equation (3.97)] because the function H there is independent of x, y, z. 3 A wave front is a surface in 3d space described by ψ(x) = c t where ψ is a function, c is the speed of light in vacuum and t is time. This means the points x on the wave front that are travelling for a time t are located on the surface ψ(x) = c t; see [KK65, Chapter II, Sec. 2].

6 6 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI p (x) I x Π ν II p 2 (m) m Figure. Snell s law p 2 (m) = p (x) + λ ν. That is, if l is the line passing through p (x) with direction ν, then p 2 (m) l Σ 2. So if there were m m 2 satisfying (3.), then p 2 (m i ) l Σ 2, i =, 2. The outer normal to Σ 2 at p ( 2(m i ) equals to N 2 p2 (m i ) ) ( = p 2 p2 (m i ) ) = m i because p 2 : Σ 2 Σ 2 is a homeomorphism since Σ 2 is C and strictly convex. And also p 2 (m ) p 2 (m 2 ). On the other hand, since the points p 2 (m i ) are also on the line l and Σ 2 is strictly convex, the normals at p 2(m i ) must satisfy that its dot products with ν must have different signs. Therefore, it cannot both happen m ν 0 and m 2 ν 0. Therefore, only one of the two m i can satisfy m i ν 0. Physically, the norm N(x) = represents the location of the points x after traveling for a time t, with c t =, from the origin into the given medium. For example, if we are in an homogenous and isotropic medium with refractive index n, then the wave propagates from the origin with velocity v = c/n. So if x satisfies N(x) =, then the Euclidean distance from O to x must satisfy x /t = v. Since t = /c, we obtain x = v /c = /n so N(x) = n x. Therefore, if medium I

7 REFRACTORS IN ANISOTROPIC MEDIA 7 has refractive index n and medium II has refractive index n 2, then N (x) = n x x and N 2 (x) = n 2 x. We then have p i (x) = N i (x) = n i, i =, 2, x 0, and so from x (3.) we recover the standard Snell law: the unit incident direction x is refracted into the unit direction m when n x n 2 m ν, see [GH09, Formula (2.)]. We shall prove that (3.) is equivalent to Fermat s principle of least time with respect to the norms N, N 2 where Σ and Σ 2 strictly convex. Suppose two anisotropic media I, II with norms N in I and N 2 in II, are separated by a plane Π as in Figure. Given X I and Y II, then Fermat s principle states that the (minimal) optical path from X to Y through the plane Π is the path XP 0 Y where P 0 Π is the unique point such that (3.2) min{n (P X) + N 2 (Y P) : P Π} = N (P 0 X) + N 2 (Y P 0 ). Equation (3.2) implies that (3.3) p 2 (Y P 0 ) p (P 0 X) ν where ν is the normal to Π. In fact, for P Π we can write P = P 0 + n i= t i e i where e,, e n is a basis for Π. From (3.2) n n t j N P 0 X + t i e i + N 2 P 0 Y + t i e i i= n n = p P 0 X + t i e i e j + p 2 P 0 Y + t i e i e j = 0 i= when t j = 0, for j =,, n. Since p i ( X) = p i (X), (3.3) follows. Also since p i are homogenous of degree zero, from (3.3) we then obtain (3.) with m = Y P 0 N 2 (Y P 0 ) and x = P 0 X N (P 0 X). Vice versa, from (3.) we deduce (3.2). In fact, let us fix X I and Y II, and consider a path XP 0 Y with P 0 Π. The ray from X to P 0 has direction x = P 0 X N (P 0 x), and the ray from P Y P 0 0 to Y has direction m =. If x is ( ) ( N) 2 (Y P 0 ) Y P 0 P 0 X refracted into m, then from (3.), p 2 p ν, and since p i N 2 (Y P 0 ) N (P 0 x) i= i=

8 8 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI are homogeneous of degree zero (3.4) p 2 (Y P 0 ) p (P 0 X) ν. Consider the functional F(P) = N (P X) + N 2 (Y P) for P Π. Write P Π as P = P 0 + n i= t i e i where e,, e n is a basis for Π. We can write F(P) = F(t,, t n ). As before and from (3.4) n n F(t,, t n ) = p t j P 0 X + t i e i e j + p 2 P 0 Y + t i e i e j = 0 i= for j =,, n and t = = t n = 0. On the other hand, since the functional F on Π is strictly convex, there are unique t 0,, t0 such that the minimum of n F is attained at t 0,, t0 n. Therefore, t0 = = t0 n = 0 and the point P 0 is the unique point in Π minimizing F, obtaining Fermat s principle. 3.. Physical constraints for refraction. Since the incident ray x is in media I we must have x ν 0, where ν is the normal to the hyperplane separating I and II, ν having direction from media I to media II. Similarly, since the refracted ray m is in media II, we also have m ν 0; Figure. We analyze here the meaning of these two physical constraints x ν 0 and m ν 0 in the following two cases. Case : Let us assume first that Σ is contained inside Σ 2, i.e., N 2 (x) for all N (x) =. If x Σ, then p (x) Σ and from Lemma 2.(b) x p (x) =. Hence, if x Σ and m Σ 2, then x (p (x) p 2 (m) ) = x p 2 (m) = N 2 (x) x N 2 (x) p 2(m) N 2 (x) i= from Lemma 2.(a) 0 from the assumption Σ contained inside Σ 2. Thus, if x ν 0 and x is refracted into m, then from (3.) p (x) p 2 (m) = λ ν and so λ > 0. Hence, m ν 0 and Lemma 2.(b) imply that m p (x) = m (p (x) p 2 (m) ) = λ m ν 0.

9 REFRACTORS IN ANISOTROPIC MEDIA 9 Therefore, when Σ is contained in the interior of Σ 2, we obtain the physical constraint (3.5) m p (x), for refraction of x Σ into m Σ 2. Case 2: Let us assume now that Σ 2 is contained inside Σ, i.e., N (x) for all N 2 (x) =. Reasoning as in the first case, we have for x Σ and m Σ 2, that m (p 2 (m) p (x) ) = m p (x) = N (m) m N (m) p (x) N (m) from Lemma 2.(a) 0 from the assumption Σ 2 contained inside Σ. Thus, if m ν 0 and x is refracted into m, then from (3.) p 2 (m) p (x) = λ ν and so λ > 0. Hence, x ν 0 and Lemma 2.(b) imply that x p 2 (m) = x (p 2 (m) p (x) ) = λ x ν 0. Therefore, when Σ 2 is contained in the interior of Σ, we obtain the physical constraint (3.6) x p 2 (m), for refraction of x Σ into m Σ 2. Notice that, as explained before, if medium I has refractive index n and medium x II has refractive index n 2, then N (x) = n x and N 2 (x) = n 2 x ; and p i (x) = n i x, i =, 2. Hence, if n > n 2 we are in Case above, and so (3.5) reads x m n 2 /n for x, m unit vectors. If n < n 2 then we are then in Case 2, and so (3.6) reads x m n /n 2 for x, m unit vectors. Therefore, when the media I and II are homogenous and isotropic we recover the physical constraints showed in [GH09, Lemma 2.].

10 0 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI 4. Uniformly refracting surfaces In this section, we shall describe the surfaces separating two anisotropic materials I and II, like in Section 3, so that rays emanating from a point source, the origin, located in medium I are refracted in medium II into a fixed direction m Σ 2. These surfaces will have the form (4.) { X R n : N (X) = p 2 (m) X + b }, where b R. If we write X = ρ(x) x for N (x) =, then the polar radius ρ(x) = b x p 2 (m) for x Σ. To show that these surfaces do the desired refraction job, as in Section 3. we distinguish two cases. Case I: Σ is strictly contained in the interior of Σ 2, that is, (4.2) κ = sup N 2 (x) <. N (x)= In this case, given m Σ 2, the desired surface is { } b (4.3) S I (m, b) = ρ(x)x : ρ(x) = x p 2 (m), x Σ with m p (x) with b > 0. In fact, to verify that each ray with direction x Σ such that m p (x) is refracted by S I (m, b) into m, we need to verify that (3.) holds, and the physical constraints x ν 0 and m ν 0 are met with ν the normal from medium I to II. From (4.), the outward normal at a point X is ν = p (x) p 2 (m) with x = X/N (X) and so (3.) holds. From Lemma 2. x ν = x (p (x) p 2 (m) ) = x p 2 (m) (4.4) = N 2 (x) x N 2 (x) p 2(m) N 2 (x) > κ > 0 from (4.2). Also m ν = m (p (x) p 2 (m) ) = m p (x) m p 2 (m) = m p (x) 0 by the definition of S I.

11 REFRACTORS IN ANISOTROPIC MEDIA Case II: Σ 2 is strictly contained in the interior of Σ, that is, (4.5) κ = inf N (x)= N 2(x) >. In this case, given m Σ 2, the desired surface is { } b (4.6) S II (m, b) = ρ(x)x : ρ(x) = x p 2 (m), x Σ with x p 2 (m) > with b > 0. In fact and once again, to verify that each ray with direction x Σ such that x p 2 (m) > is refracted by S II (m, b) into m, we need to verify that (3.) holds, and the physical constraints x ν 0 and m ν 0 are met with ν the normal towards medium II. From the definition of S II, the outward normal at a point X is ν = p 2 (m) p (x) with x = X/N (X) and so (3.) holds. From Lemma 2. m ν = m (p 2 (m) p (x) ) = m p (x) = N (m) m N (m) p (x) N (m) κ N 2(m) = κ > 0 from (4.5). Also x ν = x (p 2 (m) p (x) ) = x p 2 (m) > 0 by the definition of S II. Remark 4.. If medium I is homogeneous and isotropic with refractive index n, then N (x) = n x. Also, if II is also similar with refractive index n 2, then N 2 (x) = n 2 x. In this case, condition (4.2) is equivalent to n > n 2, and the surface S I (m, b) is a half ellipsoid of revolution with axis m, recovering the surfaces from [GH09, Formula (2.8)]. Similarly, condition (4.5) is equivalent to n < n 2, and S II (m, b) is one of the branches of a hyperboloid of two sheets as in [GH09, Formula (2.9)]. 5. The refractor problem when κ <, κ in (4.2) Using the uniformly refracting surfaces introduced in Section 4, we state and solve here the refraction problem we are interested in. We are given two closed domains Ω Σ, Ω 2 Σ 2, a non negative function f L (Ω ), and a Radon measure µ in Ω 2 satisfying the following conditions: (a) the surface measure of the boundary of Ω is zero;

12 2 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI (b) inf x Ω,m Ω 2 m p (x) ; (c) Σ and Σ 2 are C and strictly convex. Refractors are then defined as follows. Definition 5.. The surface S = {ρ(x) x : x Ω }, with ρ C(Ω ), ρ > 0, is a refractor from Ω to Ω 2, if for each x 0 Ω there exist m Ω 2 and b > 0 such that the surface S I (m, b) supports S at x 0, that is, ρ(x) b x p 2 (m) for all x Ω with equality at x = x 0. The refractor mapping associated with the refractor S is the set valued function (5.) R S (x 0 ) = {m Ω 2 : there exists S I (m, b) supporting S at x 0 }. We have the following lemma. Lemma 5.2. If a refractor S is parametrized by ρ(x) x, then ρ is Lipschitz continuous in Ω. Proof. Let x 0, x Ω and S I (m, b) supporting S at x 0. Then ρ(x) ρ(x 0 ) b x p 2 (m) b x 0 p 2 (m) = b (x x 0 ) p 2 (m) ( x p2 (m) ) ( x 0 p 2 (m) ) = ρ(x 0 ) (x x 0) p 2 (m) x p 2 (m) C x x 0 max ρ x x 0 p 2 (m) Ω κ from (4.4) since N i and are all equivalent norms, i =, 2. Reversing the roles of x and x 0 we obtain the lemma. Following the notation from [GH4, Section 2], we denote by C S (Ω, Ω 2 ) the class of set-valued maps Φ : Ω Ω 2 that are single valued for a.e. x Ω, with respect to f dx, that are continuous in Ω, and Φ(Ω ) = Ω 2. Continuity of Φ at x 0 Ω means that if x k x 0 and y k Φ(x k ), then there is a subsequence y kj and y 0 Φ(x 0 ) such that y kj y 0. Lemma 5.3. If S is a refractor from Ω to Ω 2, then the refractor map R S C S (Ω, Ω 2 ).

13 REFRACTORS IN ANISOTROPIC MEDIA 3 Proof. Let m Ω 2 and define b = max x Ω ( ρ(x) ( x p2 (m) ). From (4.4), b max x Ω ( ρ(x) ( κ) ) > 0. Also, there is x0 Ω with b = ρ(x 0 ) ( x 0 p 2 (m) ) and so m R S (x 0 ), showing that R S (Ω ) = Ω 2. Next, let us show that R S (x) is single valued for a.e. x Ω. In fact, if at x 0 Ω there exist m m 2 Ω 2 with m i R S (x 0 ), i =, 2, then x 0 is a singular point to the surface S. Otherwise, since S I (m i, b i ), i =, 2 support S at x 0, they would have the same tangent plane at x 0 ρ(x 0 ). Therefore, by the Snell law and since there is at most one m satisfying (3.), we obtain m = m 2. From Lemma 5.2, S is Lipschitz, and since Ω = 0, we obtain that R S (x) is single valued a.e. in Ω. It remains to show that R S is continuous. Let x i x 0 Ω and let m i R S (x i ). Hence ρ(x) b i x p 2 (m i ) for all x Ω with equality at x = x i. As before, b i = ( ) ρ(x i ) ( x i p 2 (m i ) ) ( κ) min Ω ρ, and b i = ρ(x i ) N 2 (x i ) N 2 (x i ) p 2(m i ) max Ω ρ (+κ) from Lemma 2.(a) and (4.2). We have m i Ω 2 Σ 2 and p 2 C (Σ 2 ). By compactness there are subsequences m ik m 0 Ω 2 and b ik b 0 > 0 so that b 0 ρ(x) x p 2 (m 0 ) for all x Ω with equality at x = x 0. This completes the proof to the lemma. Using [GH4, Lemma 2.], we obtain from Lemma 5.3 that if S is a refractor from Ω to Ω 2, then the set function (5.2) M S, f (E) = f (x) dx R S (E) is a Borel measure in Ω 2, that is called the refractor measure. Continuing using the set up from [GH4, Section 2], we recall [GH4, Definition 2.2]: given F C(Ω ) and a map T : F C S (Ω, Ω 2 ), we say T is continuous at φ F if whenever φ j F, φ j φ uniformly in Ω, x 0 Ω and y j T (φ j )(x 0 ), then there exists a subsequence y jl y 0 with y 0 T (φ)(x 0 ). If we let (5.3) F = {ρ C(Ω ) : S ρ is a refractor from Ω to Ω 2 parametrized by ρ(x)x} then we have the following lemma. x i

14 4 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI Lemma 5.4. The mapping T : F C S (Ω, Ω 2 ) defined by T (ρ) = R Sρ at each ρ F. is continuous Proof. Let ρ j, ρ F with ρ j ρ uniformly in Ω, x 0 Ω and m j R Sρj (x 0 ). Hence b j ρ j (x) x p 2 (m j ) for all x Ω with equality at x = x 0. As in the last part of the proof of Lemma 5.3, b j are bounded away from 0 and. Therefore there exist b subsequences b jl b and m jl m with ρ(x) x p 2 (m) for all x Ω with equality at x = x 0. Thus m R Sρ (x 0 ) and we are done. As a consequence of Lemma 5.4 we obtain from [GH4, Lemma 2.3] that if ρ j ρ uniformly in Ω, then M Sρj, f M Sρ, f weakly. In addition, properties (A)-(A3) from [GH4, Section 2.] translate to the present case as follows: (A) if S ρ and S ρ2 are refractors from Ω to Ω 2, then S ρ ρ 2 is a refractor from Ω to Ω 2 with ρ ρ 2 = min{ρ, ρ 2 }; (A2) if ρ (x 0 ) ρ 2 (x 0 ), then R Sρ (x 0 ) R Sρ ρ 2 (x 0 ); b (A3) we let h b,m (x) =, and we have x p 2 (m) {h b,m (x) : m Ω 2, 0 < b < } F, with F defined by (5.3). In addition we have the following (a) m R Shb,m (x) for all x Ω, from Section 4, Case I; (b) h b,m h b2,m for b b 2 ; (c) h b,m 0 uniformly in Ω as b 0; (d) h b,m h b0,m uniformly in Ω as b b 0. We then introduce the following definition. Definition 5.5. Let f L (Ω ) and let µ be a Radon measure in Ω 2 with Ω f dx = µ(ω 2 ). The refractor S from Ω to Ω 2 is a weak solution of the refractor problem if M S, f (E) = µ(e)

15 REFRACTORS IN ANISOTROPIC MEDIA 5 for each Borel set E Ω 2, where M S, f is the refractor measure defined by (5.2). Using the above set up and the existence results from [GH4, Section 2] we obtain the following theorems showing solvability of the refractor problem for anisotropic media when κ <. We first show solvability when the measure µ is discrete. Theorem 5.6. Let f L (Ω ) with f > 0 a.e., m,, m N Ω 2 be distinct points, and g,, g N positive numbers satisfying f dx = N Ω i= g i. Then for each 0 < b < there exist unique positive b 2,, b N such that S = {ρ(x)x : x Ω } with ρ(x) = min i N h b i,m i (x), is a weak solution to the refractor problem. In addition, M S, f ({m i }) = g i for i N. Proof. To prove this theorem, we use [GH4, Theorem 2.5] with the set up from above, for which we need to verify that the assumptions of that theorem are met. In fact, we need to show that we can choose positive numbers b 0,, 2 b0 such that N ρ 0 (x) = min i N h b 0 i,m i (x) such that M S ρ0, f (m i ) g i for 2 i N, where b 0 = b. We have h b 0,m (x) = b 0 x p 2 (m ) b κ from (4.4). Also h b 0 i,m (x) = b 0 i i x p 2 (m i ) b 0 i + κ for 2 i N from Lemma 2.(b) and (4.2). Therefore choosing b0,, 2 b0 N suitable so that ρ 0 (x) = h b 0,m (x), the assumptions of [GH4, Theorem 2.5] are met and the existence follows. The uniqueness follows from [GH4, Theorem 2.7] since f > 0 a.e. We are now ready to prove the following existence theorem for a general Radon measure µ. Theorem 5.7. Let f L (Ω ) with f > 0 a.e, and let µ be a Radon measure in Ω 2 such that f (x) dx = µ (Ω Ω 2 ). Then for each x 0 Ω and R 0 > 0, there exists S weak solution to the refractor problem passing through the point X 0 = R 0 x 0. Proof. Let µ l = N i= g i δ mi be a sequence of discrete measures with µ l µ weakly and µ l (Ω 2 ) = µ (Ω 2 ) for l =, 2,. From Theorem 5.6 and for the measure µ l,

16 6 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI there exists a refractor S ρ parametrized by ρ. Notice that S l l C l ρ is also a solution l to the same refractor problem since R Cl ρ = R ρ for each positive constant C l. Then l l pick C l so that C l ρ l (x 0) = R 0. Now we use the existence result [GH4, Theorem 2.8], and in order to do that we need to verify that the hypotheses (i) and (ii) of that theorem hold in the present case. To verify (i) we show that if R Range ( h b,m ), then κ R + κ h + κ b,m R κ. b In fact, there exists x Ω with R = h b,m (x ) = x p 2 (m), so h b,m(x) = b x p 2 (m) = R x p 2 (m), and the desired inequalities follow from Lemma x p 2 (m) 2.(b) and (4.2). The verification of (ii), that is, the family {ρ F : C 0 ρ C } is compact in C(Ω ), follows from Lemma 5.2 and the proof of Lemma 5.4. Remark 5.8. In the same way we can state and solve the refractor problem when κ >, i.e., (4.5) holds, using instead the uniformly refracting surfaces S II (m, b) b defined by (4.6). Now the functions h b,m are defined by h b,m (x) = x p 2 (m) and the properties (A)-(A3) defined after Lemma 5.4 must be changed in accordance with properties (A )-(A3 ) in [GH4, Section 2.2]. All lemmas in this section then hold true with obvious changes. For the existence of solutions we now need to use [GH4, Theorems 2.9 and 2.]. In Theorem 5.7 uniqueness follows from optimal mass transport, Section Propagation of light in anisotropic materials We begin this section with some background on the propagation of light in anisotropic materials. Let us assume we have a material whose permittivity and permeability are given by positive definite and symmetric matrices ɛ(x, y, z) and µ(x, y, z), respectively. Assuming we are in the geometric optics regime, i.e., the wave length of the radiation is very small compared with the objects considered, it is known [KK65, Chap. III, Sect. 4] that the function ψ = ψ(x, y, z) defining the wave fronts ψ(x, y, z)=constant, satisfies the following first order pde, the Fresnel

17 REFRACTORS IN ANISOTROPIC MEDIA 7 differential equation: ɛ R (6.) det R µ = 0, where R is the 3 3 skew-symmetric matrix 0 ψ z ψ y R = ψ z 0 ψ x. ψ y ψ x 0 We can re write Fresnel s equation in a simpler form using the following Schur s determinant identity: if A is an n n invertible matrix, B is n m, C is m n and D is m m, then A B det C D = det A det ( D C A B ) = det D det ( A B D C ), of course for the last identity D is invertible. We then get ɛ R det R µ = det ɛ det ( µ + R ɛ R ) = det µ det ( ɛ + R µ R ) and since ɛ, µ are positive definite, (6.) is equivalent to either (6.2) det ( µ + R ɛ R ) = 0, or (6.3) det ( ɛ + R µ R ) = 0. Letting (6.4) τ = µ /2 ɛµ /2, τ is symmetric and positive definite, so there is an orthogonal matrix O and a diagonal matrix D such that τ = ODO t.

18 8 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI For a column vector v = v v 2 v 3 define 0 v 3 v 2 Skew(v) = v 3 0 v. v 2 v 0 Given a 3 3 matrix B we have the formula B t Skew(B v) B = det B Skew(v). We then re write (6.3) as follows: det ( ɛ + R µ R ) = det µ det ( µ /2 ɛµ /2 + ( µ /2 Rµ /2) ( µ /2 Rµ /2)). Also since µ /2 is symmetric, we have µ /2 Rµ /2 = µ /2 Skew( ψ)µ /2 = µ /2 Skew ( µ /2 µ /2 ψ ) µ /2 = det ( µ /2) Skew ( µ /2 ψ ). Also µ /2 ɛµ /2 + ( µ /2 Rµ /2) ( µ /2 Rµ /2) = ODO t + OO ( t µ /2 Rµ /2) OO ( t µ /2 Rµ /2) OO t = O ( D + O ( t µ /2 Rµ /2) OO ( t µ /2 Rµ /2) O ) O t. We have R := O t ( µ /2 Rµ /2) O = det ( µ /2) O t Skew ( µ /2 ψ ) O = det ( µ /2) O t Skew ( O O t µ /2 ψ ) O = det ( µ /2) det O Skew ( O t µ /2 ψ ),

19 REFRACTORS IN ANISOTROPIC MEDIA 9 so D + R R = D + ( det ( µ /2)) 2 (det O) 2 Skew ( O t µ /2 ψ ) Skew ( O t µ /2 ψ ) = D + ( det ( µ /2)) 2 ( Skew O t µ /2 ψ ) Skew ( O t µ /2 ψ ) ) ) µ = D + Skew (O /2 t det ( µ µ /2) ψ Skew (O /2 t det ( µ /2) ψ. Therefore det ( ɛ + R µ R ) ( ) )) µ = det µ det D + Skew (O /2 t det ( µ µ /2) S Skew (O /2 t det ( µ /2) S (6.5) = det µ det ( D + Skew(p) Skew(p) ) = 0, where p = p p 2 p 3 := O t µ /2 det ( µ /2) S. Notice that this calculation is done at a fixed point (x, y, z) since the matrices ɛ and µ depend on the point (x, y, z); therefore the matrices D and O depend also on (x, y, z). Next we have Skew(p) Skew(p) = p p t p t p Id = p p (p p) Id, so by (6.5) the Fresnel equation for the wave fronts (6.3) is then 0 = det ( D + Skew(p) Skew(p) ) = det ( D + p p (p p) Id ). To write this equation in a more convenient form, set τ 0 0 D = 0 τ 2 0, 0 0 τ 3

20 20 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI (a matrix depending on (x, y, z)), so 0 = det ( D + p p (p p) Id ) τ + p 2 p 2 p p 2 p p 3 = det p 2 p τ 2 + p 2 2 p 2 p 2 p 3 p 3 p p 3 p 2 τ 3 + p 2 3 p 2 τ p 2 2 p2 p 3 p 2 p p 3 = det p 2 p τ 2 p 2 p2 3 p 2 p 3 p 3 p p 3 p 2 τ 3 p 2 p2 2 Let us now define for an arbitrary vector (p, p 2, p 3 ) the following functions, which depend on the point (x, y, z) since τ i depend on (x, y, z) Φ(p, p 2, p 3 ) = ( + ) p 2 2 τ2 τ + ( + ) p τ τ + ( + ) p τ τ, 2 and ( Ψ(p, p 2, p 3 ) = (p 2 + p2 2 + p2 3 ) p 2 τ 2 τ + p τ τ + ) p τ τ 2 It is easy to check that (6.6) τ p 2 2 p2 p 3 p 2 p p 3 det p τ τ 2 τ 2 p τ 2 p 2 p2 p 3 2 p 3 = 2Φ(p, p 2, p 3 ) + Ψ(p, p 2, p 3 ). 3 p 3 p p 3 p 2 τ 3 p 2 p2 2 Now write Next notice that 2Φ + Ψ = 2Φ + Φ 2 Φ 2 + Ψ = ( Φ) 2 (Φ 2 Ψ). (6.7) Φ 2 Ψ, which follows using Legrange multipliers since Φ 2 Ψ is homogenous of degree four. So we can write 2Φ + Ψ = ( Φ Φ 2 Ψ ) ( Φ + Φ 2 Ψ )..

21 REFRACTORS IN ANISOTROPIC MEDIA 2 We then obtain that the Fresnel equation of wave fronts (6.3) can be split as the following two equations (6.8) Φ Φ 2 Ψ = 0 or Φ + Φ 2 Ψ = 0. Each of these equations describes a three dimensional surface that depends of the point (x, y, z) chosen at the beginning; see Figure 2. That is, in this way each point (x, y, z) in the space has associated a pair of surfaces, one enclosing the other. The inner surface is convex and the outer surface is neither convex nor concave. We have shown that the vector (6.9) p = p p 2 p 3 := O t µ /2 det ( µ /2) ψ, belongs to one of the surfaces, with all quantities calculated at (x, y, z), and the matrix O is orthogonal and diagonalizes the matrix τ. In other words, we have shown that the gradient ψ(x, y, z) of the wave front ψ=constant, when multiplied µ /2 by the matrix det ( µ /2) and conveniently rotated by Ot, belongs to one of the surfaces described by the equations (6.8). Notice that when the permittivity matrix is ɛ Id and the permeability matrix is µ Id, where ɛ and µ are scalar functions depending only on position, then we recover the eikonal equation ψ 2 = ɛ µ. In fact, in this case the matrix τ = (ɛ/µ) Id, so τ i = ɛ/µ for i =, 2, 3, O = Id, Φ(p, p 2, p 3 ) = µ ɛ (p2 + p2 2 + p2 3 ), and ( µ Ψ(p, p 2, p 3 ) = ɛ ) 2 (p 2 + p2 2 + p2 3 )2. So Φ 2 = Ψ and both surfaces in (6.8) are identically equal to Φ(p, p 2, p 3 ) = µ ɛ (p2 + p2 2 + p2 3 ) = 0.

22 22 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI Figure 2. Fresnel surfaces when τ =, τ 2 = 2, τ 3 = 3; µ = Id So the vector p in (6.9) satisfies the last equation and p = µ ψ, therefore ψ 2 = ɛ µ. 6.. Case considered for the application of our results. For the application of our results from Sections 3 5 we consider materials having permittivity and permeability tensors ɛ and µ that are positive definite symmetric constant matrices with µ = a ɛ, where a is a positive number. These are homogeneous materials that when ɛ is not the identity matrix are anisotropic. We will associate with such a material a norm as follows. From the calculations above, the Fresnel equation in this case is as follows. From (6.4) we get τ = Id, so a Φ(p) = a p 2, Ψ(p) = a 2 p 4. Obviously, Φ 2 = Ψ and so the Fresnel equation is Φ = 0, i.e., p 2 = /a and therefore it has only one sheet. Then from (6.9) the vector ψ satisfies the det µ /2 µ /2

23 REFRACTORS IN ANISOTROPIC MEDIA 23 equation a µ /2 det µ /2 ψ 2 =. The last expression induces the following dual norm N (p) = a µ /2 det µ p /2. The norm N is the dual to the norm given by N(x) = sup x p = N (p)= det µ/2 a µ /2 x, which is the norm we associate to the material. Notice that if µ is the identity matrix, then ɛ = Id, the material is isotropic and has index of refraction n = a ɛ µ = / a. The norm obtained this way is then N(x) = n x, in agreement with the physical explanation for isotropic media given after (3.). Now, if N(x) = A x with A a constant matrix, then N(x) = N(x) At Ax. Therefore, having two materials I and II so that the wave fronts are given by norms N (x) = A x and N 2 (x) = A 2 x, respectively, the Snell law (3.) takes the following form: Each incident ray traveling in medium I with direction x Σ, i.e., N (x) =, with x ν 0 and striking the plane Π at some point P 0 is refracted in medium II into a direction m Σ 2, i.e., N 2 (m) =, if A t 2 A 2m A t A x ν, where ν is the unit normal at P 0 from medium I to medium II. In our application we have materials I and II having constant tensors (ɛ, a ɛ ) and (ɛ 2, a 2 ɛ 2 ), respectively, and therefore the associated norms to I and II are N (x) = det µ/2 a µ /2 x, N2 (m) = det µ/2 2 a2 µ /2 2 m, respectively. If we let A = det µ/2 a µ /2, A 2 = det µ/2 2 a2 µ /2 2,

24 24 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI then N (x) = A x and N 2 (m) = A 2 m. To apply the results of the previous sections to this case, from the definition of κ in (4.2) we have κ = sup A 2 x = sup A 2 A z A2 = A A x = z = the norm of the matrix A 2 A induced by the standard Euclidean norm. Hence when A2 A < the results from Section 5 are applicable to this case. On the other hand, when κ is defined by (4.5), setting A = A 2 A we get minimum eigenvalue of A t A = inf A 2 A z >, and the results from Remark 5.8 are applicable in this case. We can also write A 2 A = det µ/2 2 a2 µ /2 2 a det µ /2 µ /2 = z = a det µ /2 2 µ /2 µ /2. a 2 det µ /2 2 Once again notice that if µ i is the identity matrix, then ɛ i = a i Id, the materials are isotropic and have index of refraction n i = ɛ i µ i = / a i. The norms are then N i (x) = n i x, i =, 2, and κ = a /a 2 = n 2 /n in agreement with the physical explanation for isotropic media given after (3.). Finally, we remark that for the materials considered light rays travel in straight lines and they do not exhibit bi refringence, that is, each incident ray is refracted into only one ray. one sheet. The last property is because the Fresnel equation has only That rays travel in straight lines follows from Fermat s principle of least time explained in Section 3. Indeed, let X, Y be two points in space, γ(θ) = ( θ) X + θ Y, 0 θ, and let φ(θ) be any curve from X to Y. Then the optical length T for each curve satisfies T(γ) = 0 γ (θ) dθ = X Y, and T(φ) = 0 φ (θ) dθ 0 φ (θ) dθ = X Y. For general anisotropic materials when ɛ is not a multiple of µ, the Fresnel equation has two sheets, see Figure 2, and as mentioned before bi-refringence. 4 That is, κ is the spectral norm of the matrix A := A 2 A, i.e., κ = maximum eigenvalue of (A t A) 4

25 REFRACTORS IN ANISOTROPIC MEDIA 25 occurs. This is the case for crystals, that is, when ɛ is a diagonal constant matrix and µ = Id. 7. Connection with optimal mass transport The setting up, analysis, and results from the previous sections allow us to cast the refraction problem in optimal transport terms. However, the method used in Section 5 to prove existence of solutions relies more on a deeper insight of the physical and geometric features of the refractor problem. To apply the optimal mass transport approach, we use the abstract set up in [GH09, Section 3.2 and 3.3] and from Definition 5. introduce the cost function ( ) c(x, m) = log x p 2 (m) for x Ω Σ, m Ω 2 Σ 2 and κ < with κ in (4.2). With [GH09, Definition 3.9] of c-concavity, we have that S = {ρ(x)x : x Ω } is a refractor in the sense of Definition 5. above if and only if log ρ is c-concave. From the definition of c-normal mapping N c,φ given in [GH09, Definition 3.0], and the definition of refractor mapping R S given by (5.), we have that R S = N c,log ρ. One can easily check that S is a weak solution of the refractor problem if and only if log ρ is c- concave and N c,log ρ is a measure preserving map in the sense of [GH09, Equation (3.9)] from f (x) dx to µ. Hence existence and uniqueness up to dilations of the refractor problem follows as in [GH09, Theorem 3.5]. From Remark 5.8 and using the cost function c(x, m) = log(x p 2 (m) ) we obtain similar results for κ >, defined by (4.5). References [BW59] M. Born and E. Wolf. Principles of Optics, Electromagnetic theory, propagation, interference and diffraction of light. Cambridge University Press, seventh (expanded), 2006 edition, 959. [CGH08] L. A. Caffarelli, C. E. Gutiérrez, and Qingbo Huang. On the regularity of reflector antennas. Ann. of Math., 67: , 2008.

26 26 C. E. GUTIÉRREZ, QINGBO HUANG AND H. MAWI [CH09] L. A. Caffarelli and Qingbo Huang. Reflector problem in R n endowed with non- Euclidean norm. Arch. Rational Mech. Anal., 93(2): , [GH09] C. E. Gutiérrez and Qingbo Huang. The refractor problem in reshaping light beams. Arch. Rational Mech. Anal., 93(2): , [GH4] C. E. Gutiérrez and Qingbo Huang. The near field refractor. Annales de l Institut Henri Poincaré (C) Analyse Non Linéaire, 3(4): , July-August temple.edu/ gutierre/papers/nearfield.final.version.pdf. [GM3] C. E. Gutiérrez and H. Mawi. The far field refractor with loss of energy. Nonlinear Analysis: Theory, Methods & Applications, 82:2 46, 203. [Gut4] C. E. Gutiérrez. Refraction problems in geometric optics. In Lecture Notes in Mathematics, vol. 2087, pages Springer-Verlag, 204. [Kar6] A. Karakhanyan. An inverse problem for the refractive surfaces with parallel lighting. SIAM J. Math. Anal., 48(): , 206. [KK65] M. Kline and I. W. Kay. Electromagnetic theory and geometrical optics, volume XII of Pure and Applied Mathematics. Wiley, 965. [LGM7] R. De Leo, C. E. Gutiérrez, and H. Mawi. On the numerical solution of the far field refractor problem. Nonlinear Analysis: Theory, Methods & Applications, 57:23 45, 207. [LL84] L. D. Landau and E. M. Lifshitz. Electrodynamics of Continuous Media, volume 8 of Course of Theoretical Physics. Pergamon Press, 2nd revised and enlarged edition, 984. [Sch07] Toralf Scharf. Polarized light in liquid crystals and polymers. Wiley, [Som54] Arnold Sommerfeld. Optics, volume IV of Lectures on theoretical physics. Academic Press, 954. [YY84] A. Yariv and P. Yeh. Optical waves in crystals. John Wiley & Sons, 984. Department of Mathematics, Temple University, Philadelphia, PA 922 address: gutierre@temple.edu Department of Mathematics & Statistics, Wright State University, Dayton, OH address: qingbo.huang@wright.edu Department of Mathematics, Howard University, Washington, D.C address: henok.mawi@howard.edu