arxiv: v1 [math.ap] 2 May 2018

Transcription

1 . REFLECTION AND REFRACTION PROBLEMS FOR METASURFACES RELATED TO MONGE-AMPÈRE EQUATIONS CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI arxiv: v [math.ap] May 08 Abstract. A metasurface is a surface, tipically a plane, on which a function called phase discontinuity is chosen so that the metasurface produces a desired reflection or refraction job. We derive the equations that the phase discontinuity function must satisfy, they are Monge-Ampère partial differential equations, and we prove eistence of solutions. Contents. Introduction.. Background 3. Reflection 4.. Collimated case 4.. Point Source Reflection 7 3. Refraction 3.. Collimated case 3.. Point Source Refraction 3 4. Conclusion 7 References 8. Introduction The subject of metalenses is a flourishing area of research having multiple potential applications and is one of the nine runners-up for Science s Breakthrough of the Year 06 []. A typical problem is given a surface, in most applications a plane, find a function on the surface a phase discontinuity so that the pair, surface together with the phase discontinuity the metalens refracts or reflects light in a desired manner. This leads to ultra thin optical components that produce abrupt changes over the scale of the free-space wavelength in the phase. This is in contrast with classical lens design, where the question is to engineer the gradual Research supported by NSF grant DMS

2 CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI accumulation of phase delay as the wave propagates in the device, reshaping the scattered wave front and beam profile at will. In particular, in lenses light propagates over distances much larger than the wavelength to shape wavefronts. A vast literature appeared in recent years on this subject, see for eample [], [3], [4], and the comprehensive review articles [5] and [6]. For more recent work in the area and an etensive up to date bibliography, we refer to [7]; see also [8]-[8] and [9]. In our previous paper [0] we gave a mathematically rigorous foundation to deal with general metasurfaces and to determine the relationships between the curvature of the surface and the phase discontinuity. The present paper builds on that work and shows that phase discontinuities theoretically eist so that there are metalenses that can refract or reflect energy with prescribed energy patterns. More precisely, a question considered and answered in this paper is the following. Let Ω be a set on the, y-plane, Ω a set of unit vectors in R 3, and Γ = {z = a} a horizontal plane above the, y-plane. We are given two intensities, i.e, two non negative functions, f defined in Ω and g defined in Ω satisfying the energy conservation condition f, y ddy = gz dσz, Ω Ω where dσz denotes as usual the element of area in the unit sphere of R 3. A collimated beam is emanating from Ω. That is, for each, y Ω a ray is emitted in the vertical direction e 3 = 0, 0, with intensity f, y and strikes the plane Γ at the point, y, a = P. According to the generalized law of reflection.4 this ray is reflected by the meta surface Γ, ψ into a ray having direction T, y = e 3 λ e 3 ψ, y, a, since νp = e 3. The question is then to find a function ψ defined in Γ, called a phase discontinuity, such that the meta lens, i.e., the pair Γ, ψ reflects all rays from Ω into Ω, that is, TΩ = Ω, and the energy conservation balance. E f, y ddy = TE gz dσz holds for each subset E of Ω, see Figure. We show in Section. that this problem is mathematically solvable. We also consider a similar question when the incident rays emanate from a point source into a set of unit directions Ω, see Figure. Such a problem is solved in Section..

3 REFLECTION AND REFRACTION PROBLEMS FOR METASURFACES 3 In addition, we consider and solve similar problems for refraction using the generalized law. both in the collimated and point source cases, see Figures 3 and 4, Sections 3. and 3.. It is our purpose to show that each of these problems has a theoretical solution. To do this we derive the partial differential equation, for each problem, satisfied by the phase discontinuity ψ and show it is a Monge-Ampère equation. Net we show that the resulting equations have solutions by application of a result by Urbas []. The equations corresponding to the four problems considered are.9,.0, 3. and 3.3, and they can be regarded as particular cases of 3.4. Monge-Ampère equations appear naturally in optics for freeform lens design that have been the subject of recent research, see for eample []-[3]. Therefore it is natural that these type of equations appear also for metasurfaces. Monge- Ampère equations have been recently the subject of important mathematical research due to their connections with various topics such us optimal mass transport. We refer the reader to [3] and [33] for details and references therein. We hope our contributions in this paper may be useful to understand theoretically what kind of phase discontinuities are possible to design in the applications of metasurfaces. We mention that recent work using the ideas from [0] for reflection is done in [34] to design graphene-based metasurfaces that can be actively tuned between different regimes of operation, such as anomalous beam steering and focusing, cloaking, and illusion optics, by applying electrostatic gating without modifying the geometry of the metasurface. Finally, if the surface Γ is not necessarily a plane, then is possible to derive the corresponding partial differential equation that the phase discontinuity ψ satisfies, in both the reflection and refraction cases. These are equations of Monge-Ampère type that require a more complicated derivation and analysis beyond the scope of this article and will appear as a part of the forthcoming PH.D. thesis of the second author... Background. We begin describing the formulation of the generalized laws of reflection and refraction with phase discontinuities from [0]. For refraction consider two homogenous and isotropic media I and II with refractive indices n and n respectively. Suppose we have an interface surface Γ separating media I and II, and a function ψ, called the phase discontinuity, defined in a small neighborhood of Γ. If a ray with unit direction emanating from medium I strikes the surface Γ at some point P, then it is refracted into medium II into a ray

4 4 CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI having unit direction m such that. n n m = λ νp + ψp, see [0, Formula 6], where νp is the outer unit normal to the surface Γ at P, ψ denotes the gradient, and λ is a constant depending on, νp, ψp and the refractive indices. In fact, it is proved in [0, Formula ] that λ can be calculated:.3 λ = n ψ ν n n ψ + n ψ ν. This law is derived in [0, Section 3] using wave fronts. The job of the function ψ, concentrated around Γ, is to change the direction of the incoming rays. As a difference from the standard Snell law of refraction, the rays here are mainly bent due to the function ψ and not by the change in the refractive indices of the surrounding media. In particular, the generalized law of refraction. makes sense when n = n. When ψ is constant and n n, we obtain the standard Snell law of refraction. The case of reflection is when n = n and since now the reflected vector must be on the same side of Γ, that is, m ν 0, the generalized reflection law has the form.4 m = λ/n νp + ψ/n P, with the unit incident ray, m the unit reflected ray, νp the normal to Γ at P, and.5 λ = n ψ ν + n n ψ + n ψ ν, see [0, Section 3, Remark ]. standard reflection law. Once again when ψ is constant this yields the. Reflection.. Collimated case. Here we solve the first problem stated in the introduction. From.4 and.5 with n =, the vertical ray emanating from the point, y Ω is reflected by the metasurface Γ, ψ into the unit direction.6 T, y = i, y λν, y ψ, y, where i, y = 0, 0,, ν, y is the normal to Γ = {z = }, and λ = i ψ ν + i ψ [i ψ ν] = + ψ ψ y.

5 REFLECTION AND REFRACTION PROBLEMS FOR METASURFACES 5 We remark that in the last identity we have used that ψ z = 0 since we seek a phase discontinuity ψ tangential to the surface Γ. Γ,ψ m Ω Figure. Reflection from an etended source rays are monochromatic; colors are used only for visual purposes. Therefore T, y = T, T, T 3 = ψ, y, ψ y, y, ψ, y ψ y, y. From the conservation of energy condition. and the formula of change of variables for surface integrals.7 f d = gy dσy = gtz J T z dz, E TE for each open set E Ω, and where J T = T, y T y, y. From.7 we obtain.8 f = gt, y J T, y for, y Ω. To calculate J T, y, since T, y =, differentiating with respect to and y yields the equations T T = T T y = 0. Hence, assuming T 3, y 0 and solving these equations in T 3 and T 3 y we get T 3 = T T + T T T 3 and T 3 y = T T y + T T y T 3. E

6 6 CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI From an elementary calculation with the determinant we obtain Hence T T y = T T y det T. T 3 T T y J T = T 3, y detd ψ, where D ψ is the matri of the second derivatives in and y. Therefore from.8 the phase discontinuity ψ satisfies the following Monge-Ampère equation.9 detd ψ = ψ, y ψ y, y f, y g T, y. To show that.9 has solutions we invoke [, Theorem ], which in simpler terms says the following: if D, D are uniformly conve smooth domains in R n, f > 0 is a smooth function in D, f > 0 is a smooth function in D, satisfying.0 f d = f p dp, D then the boundary value problem D detd u = f f u in D, ud = D, has a conve C solution, and any two such solutions differ by a constant. In fact, to apply this result to show eistence of solutions to.9, set n =, let f, y = f, y for, y D = Ω, g p, p, p p f p, p =, p p for p, p D = ΠΩ, where Π is the orthogonal projection from a set on the unit sphere onto the, y-plane. In particular, Ω is a subset of the lower unit hemisphere z 0. We need to verify.0. From the conservation of energy assumption f d = gy dσy, Ω Ω

7 REFLECTION AND REFRACTION PROBLEMS FOR METASURFACES 7 and using the parametrization q = q, q q, q we can write gy dσy = Ω = ΠΩ ΠΩ g q, q dq q g p, p p dp = D f p, p dp. Therefore.0 holds and hence the eistence of solutions to.9 follows. Γ,ψ m O Figure. Reflection from a point source.. Point Source Reflection. We now have a domain Ω of the unit sphere in R 3 and rays emanate from the origin with intensity f 0 for each Ω. Ω is as in the previous section, i.e., a domain of the unit sphere and g > 0 is a function in Ω such that the following energy conservation condition holds. f dσ = gy dσy. Ω Ω Again Γ is the plane z =. Of course, we assume that rays from the origin with unit direction in Ω reach the plane Γ. The question is then to find a phase discontinuity ψ on Γ such that all rays emitted from the origin with direction

8 8 CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI Ω and intensity f, are reflected by the meta surface Γ, ψ into Ω such that. f dσ = gy dσy E for each subset E of Ω and TΩ = Ω, where T is the reflection map. In order to find the equation ψ satisfies, we parametrize the domains in the sphere using spherical coordinates: su, v = cos u sin v, sin u sin v, cos v, 0 u π, 0 v π/. Parametrizing Ω in these coordinates we obtain Ω = so, for some domain O [0, π] [0, π/]. Re writing the integrals in. in spherical coordinates, and letting su = E, we have f su, v s u s v dudv = f dσ U E = gy dσy = gtsu, v T s u T s v dudv. TE U Since this equation must hold for all open sets U O, it follows that T satisfies the equation TE.3 T s u T s v s u s v = f su, v gtsu, v. The plane Γ is described in spherical coordinates by the polar radius.4 ru, v = su, v = cos u tan v, sin u tan v,. cos v From.4 with n =, if the incident ray has direction i = su, v, then the reflected ray that has unit direction Tsu, v = su, v λν ψru, v, where ν = 0, 0, is the normal to Γ at the incident point. Since we seek, as before, for a phase ψ tangential to Γ, we have ψ, y, = ψ, y,, ψ y, y,, 0. In addition, from.5 λ = i ν + i ψ i ν = cos v + cos u sin v ψ ru, v sin u sin v ψ y ru, v = cos v +,

9 REFLECTION AND REFRACTION PROBLEMS FOR METASURFACES 9 where = cos u sin v ψ ru, v sin u sin v ψ y ru, v. writing T in components Therefore Tsu, v = T su, v, T su, v, T 3 su, v.5 = cos u sin v ψ ru, v, sin u sin v ψ y ru, v,. Since Tsu, v =, it follows as in Section. that.6 T s u T s v = T 3 s det T s u T s v. T s u T s v From.5 T s u = sin u sin v ψ ru, v sin u tan v ψ y ru, vcos u tan v, cos u sin u T s v = cos u cos v ψ ru, v ψ cos y ru, v, v cos v T s u = cos u sin v ψ y ru, v sin u tan v ψ yy ru, vcos u tan v, cos u sin u T s v = sin u cos v ψ y ru, v ψ cos yy ru, v. v cos v Inserting these in.6 yields.7 T s u T s v = T 3 s det Au, v D,y ψru, vbu, v where sin u sin v cos u cos v Au, v =, cos u sin v sin u cos v sin u tan v cos u cos Bu, v = v. cos u tan v sin u cos v We can re write the above quantities in rectangular coordinates noticing that = cos u tan v, y = sin u tan v, ru, v =, y,, + y + = cos v and + y =

10 0 CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI tan v. We obtain T = T = T 3 = = A = + y + ψ, y, = + y + ψ, y,, y + y + ψ y, y, = + y + ψ, y, + y + ψ, y, y + y + ψ y, y, + y + ψ, y, + y + ψ, y,, y y + y + + y + y + + y + y B = y + + y. + y + y + + y y + y + + y, y, Also,.8 s u s v = cos v sin u cos v sin v cos u cos v sin v + y = sin v = + y +. Now notice that deta D ψb = detb detab D ψ, with.9 and y + + y + + y B = + y + + y + y + y y y AB = b + c y b + y c y b + y c = D + y +, b + y c

11 REFLECTION AND REFRACTION PROBLEMS FOR METASURFACES where b = + y + y + / and c = + y + y + 3/. Therefore deta D ψb = detb det D + y + ψ. Letting φ, y = + y + ψ, y, using the last equation in.7, and using.8, we obtain from.3 that φ satisfies the following equation.0 + y + 3/ detd φ, y φ, y φ y, y = f, y, +y +. g φ, y, φ y, y, φ, y φ y, y The above equation holds for, y D where D is obtained as follows: for each direction e Ω, the ray with this direction intersect the plane z = at a unique point, y, this collection of and y is D. We now proceed as in the previous section to show eistence of solutions to.0. To this end we need to identity the functions f, f in.0. Parametrizing Ω by q : D Ω with q, y =, y,, we let f f q, y, y = +y + for, y D = D. Also let f p, p = g p, p, ΠΩ. With these choices and observing that f z dσz = Ω D f p p p p, y, +y + ddy, + y + 3/ + y + 3/, for p, p D = a similar calculation as at the end of last section, we obtain that. is equivalent to.0 and therefore eistence of solutions to.0 follows as before invoking [, Theorem ]. 3. Refraction Here we solve two problems similar to the ones considered in the previous section but for refraction.

12 CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI m n Γ,ψ n Ω Figure 3. Refraction from an etended source 3.. Collimated case. Incident rays are emitted from an open set Ω of the -y plane with direction i, y = e 3 = 0, 0,, and Γ is the plane z =. From the generalized law of refraction. and.3, the meta surface Γ, ψ refracts the incident ray i, y into a ray r, y with direction satisfying n i, y n r, y = λν, y + ψ, y, where n and n are the refractive indices of the two homogeneous and isotropic media separated by the plane Γ, ν is the unit normal to the plane Γ. Also λ = n i ψ ν n n i ψ + [n i ψ ν] = n n ψ + ψ y, since we seek ψ tangential to Γ; i.e., ψ z = 0. We then let T : Ω Ω to be T, y := r, y = ψ, y, ψ y, y, ψ n n n, y + ψ y, y. We seek ψ defined on Γ with TΩ = Ω and satisfying the conservation of energy balance E f d = gy dσy = gtz J T dz for each E Ω, TE E

13 REFLECTION AND REFRACTION PROBLEMS FOR METASURFACES 3 where J T = T, y T y, y. Since T, y = and similarly as for reflection, we have that J T = T 3, y D det ψ. n Therefore proceeding as in the reflection case, the phase discontinuity ψ must satisfy the following Monge-Ampére equation 3. n ψ, y + ψ y, y det D n ψ = g n ψ, n ψ y, f, y n ψ + ψ y ; notice that this equation is independent of the value of n. Similar to the reflection case, TΩ = Ω implies that n ψ, ψ y ΠΩ where Π is once again the orthogonal projection onto the -y plane. We claim, also in this case, that [, Theorem ] can be applied to obtain a solutions ψ to 3.. Indeed, letting f, y = f, y for, y D = Ω, g n p, n p, n p + p f p, p = n p +, p for p, p D = n ΠΩ, and proceeding as before we obtain that f d = gy dσy. Ω Ω is equivalent to.0 and so eistence of solutions follows as before. 3.. Point Source Refraction. We now analyze a problem similar to the one in Section. for refraction. That is, rays emanate for a point source and we seek a phase discontinuity ψ defined on the plane Γ = {z = } so that the refraction map T to be calculated in a moment satisfies the conservation of energy condition.. As in Section., this implies.3, i.e., T s u T s v s u s v = f su, v gtsu, v, and TΩ = Ω. Let us calculate the refraction map T. As in Section., the plane Γ is described by the polar radius.4. Then from., the refracted ray has unit direction Tsu, v = n n su, v n λν n ψru, v,

14 4 CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI m n Γ,ψ n O Figure 4. Refraction from a point source where ν = 0, 0, is the normal to Γ at the incident point, su, v are spherical coordinates, ψ, y, = ψ, y,, ψ y, y,, 0 since we seek a phase discontinuity ψ tangential to Γ, and λ = n i ν + n n i ψ n i ν = n cos v + n n cos u sin v ψ ru, v n sin u sin v ψ y ru, v = n cos v +, where = n n cos u sin v ψ ru, v n sin u sin v ψ y ru, v. Therefore Tsu, v = T su, v, T su, v, T 3 su, v n = cos u sin v ψ ru, v, n sin u sin v ψ y ru, v,. n n n n n Since Tsu, v = we have as in.6 that 3. T s u T s v = T 3 s det T s u T s v T s u T s v.

15 REFLECTION AND REFRACTION PROBLEMS FOR METASURFACES 5 On the other hand, T s u = n sin u sin v ψ ru, v sin u tan v ψ y ru, vcos u tan v, n n n T s v = n cos u cos v cos u ψ ru, v sin u ψ n n cos y ru, v, v n cos v T s u = n cos u sin v ψ y ru, v sin u tan v ψ yy ru, vcos u tan v, n n n T s v = n sin u cos v ψ y ru, v cos u sin u ψ n n cos yy ru, v. v n cos v Inserting these into 3. yields T s u T s v = T 3 s det Au, v D,y ψru, vbu, v n where Au, v = n sin u sin v cos u cos v, n cos u sin v sin u cos v sin u tan v cos u cos Bu, v = v. cos u tan v sin u cos v As in the point source reflection case, Section., we can re write the above quantities in rectangular coordinates noticing that = cos u tan v, y = sin u tan v,

16 6 CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI ru, v =, y,, + y + = cos v and + y = tan v. We obtain T = n n T = n T 3 = = n A = n n + y + n ψ, y, = n n + y + ψ, y,, n y + y + n ψ y, y, = n n + y + ψ, y, n n n + y + ψ, y, n n n y + y + ψ y, y, n n n + y + n ψ, y, n n + y + ψ, y, n y + y + + y + y + + y + y B = y + + y. + y Also from.8 s u s v = + y + + y y + y + + y, +y. Now notice that +y + det A D ψb = detb det AB D ψ, n n y, y, with B as in.9 and AB = n y b + c n y b + y c y b + y c b + y = n D n + y +, c where b = + y + y + / and c = + y + y + 3/. Therefore det A D ψb = detb det D n n n + y + ψ. n

17 REFLECTION AND REFRACTION PROBLEMS FOR METASURFACES 7 Letting φ, y = n n + y + n ψ, y, we obtain that φ satisfies the following equation y + 3/ detd φ, y φ, y φ y, y = f, y, +y +. g φ, y, φ y, y, φ, y φ y, y The above equation holds for, y D where D is obtained as at the end of Section.. Eistence of solutions to this equation follows as before letting f, y = f, y, g p +y +, p, p p for, y D + y + 3/ = D, and f p, p = for p p p, p D = ΠΩ, where Π is once again the orthogonal projection. Remark 3.. If a ray is emitted from a point Q and strikes the plane Γ = {z = } at the point P =, y,, let d Q, y be the distance from Q to P. In the collimated case, since all rays are vertical d Q, y =. And when the point source Q is the origin, d Q, y = + y +. Then writing φ, y = n n d Q, y n ψ, y, and noticing that n = n = in the reflection cases, the equations.9,.0, 3. and 3.3 can be written as 3.4 d 3/, y Q detd φ, y = φ, y φ y, y f, y g T, y, where f, y = f, y in the collimated case, and f, y = f, y, +y + the point source case. in 4. Conclusion We have derived the equations that a phase discontinuity defined on a plane must satisfy in order that the resulting metasurface reflects or refracts light emanating with certain given variable intensity into a set of directions having also another variable intensity. The cases considered are when light emanates in a

18 8 CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI collimated beam and when light emanates from one point source. The resulting equations for the phase discontinuity are Monge-Ampère partial differential equations which we show have solutions. References [] The runners-up. Science, :58 53, 06. [] Nanfang Yu, Patrice Genevet, Mikhail A. Kats, Francesco Aieta, Jean-Philippe Tetienne, Federico Capasso, and Zeno Gaburro. Light propagation with phase discontinuities: generalized laws of reflection and refraction. Science, : , 0. [3] Francesco Aieta, Patrice Genevet, Nanfang Yu, Mikhail A. Kats, Zeno Gaburro, and Federico Capasso. Out-of-plane reflection and refraction of light by anisotropic optical antenna metasurfaces with phase discontinuities. Nano letters, 3:70 706, 0. [4] Francesco Aieta, Ali Kabiri, Patrice Genevet, Nanfang Yu, Mikhail A Kats, Zeno Gaburro, and Federico Capasso. Reflection and refraction of light from metasurfaces with phase discontinuities. Journal of Nanophotonics, 6: , 0. [5] Nanfang Yu and Federico Capasso. Flat optics with designer metasurfaces. Nature Materials, 3:39 50, 04. [6] Hou-Tong Chen, Antoinette J Taylor, and Nanfang Yu. A review of metasurfaces: physics and applications. Reports on Progress in Physics, 797:07640, 06. [7] Patrice Genevet, Federico Capasso, Francesco Aieta, Mohammad Rezakhorasaninejad, and Robert Devlin. Recent advances in planar optics: from plasmonic to dielectric metasurfaces. Optica, 4: , January 07 [8] M. Khorasaninejad, A. Y. Zhu, C. Roques-Carmes, W. T. Chen, J. Oh, I. Mishra, R. C. Devlin, and F. Capasso. Polarization-insensitive metalenses at visible wavelengths. Nano Lett., 6:79 734, 06. [9] Mohammadreza Khorasaninejad, Wei Ting Chen, Robert C. Devlin, Jaewon Oh, Aleander Y. Zhu, and Federico Capasso. Metalenses at visible wavelengths: Diffraction-limited focusing and sub-wavelength resolution imaging. Science, 35690:90 94, June 06. [0] E. F. Kuester, M. A. Mohamed, M. Piket-May, and C. L. Holloway,Averaged transition conditions for electromagnetic fields at a metafilm. IEEE Transactions on Antennas and Propagation : [] K. Achouri, M. A. Salem, and C. Caloz. General metasurface synthesis based on susceptibility tensors. IEEE Transactions on Antennas and Propagation : [] Habib Ammari, Matias Ruiz, Wei Wu, Sanghyeon Yu, and Hai Zhang. Mathematical and numerical framework for metasurfaces using thin layers of periodically distributed plasmonic nanoparticles. In Proc. R. Soc. A, volume 47. The Royal Society, 06. [3] S. M. Kamali, A. Arbabi, E. Arbabi, Y. Horie, and A. Faraon, Decoupling optical function and geometrical form using conformal fleible dielectric metasurfaces, Nature Communications, [4] J. Cheng, S. Jafar-Zanjani, and H. Mosallaei, All-dielectric ultrathin conformal metasurfaces:lensing and cloaking applications at 53 nm wavelength, Scientific Reports, 6 06.

19 REFLECTION AND REFRACTION PROBLEMS FOR METASURFACES 9 [5] S. Jafar-Zanjani, J. Cheng, and H. Mosallaei, Light manipulation with flat and conformal inhomogeneous dispersive impedance sheets: an efficient FDTD modeling, Applied optics, 55, [6] M. Khorasaninejad, Z. Shi, A. Y. Zhu, W. T. Chen, V. Sanjeev, A. Zaidi, and F. Capasso. Achromatic Metalens over 60 nm Bandwidth in the Visible and Metalens with Reverse Chromatic Dispersion. Nano Lett., 7 3, pp [7] A. Y. Zhu, Wei-Ting Chen, M. Khorasaninejad, J. Oh, A. Zaidi, I. Mishra, R. C. Devlin, and F. Capasso. Ultra-compact visible chiral spectrometer with meta-lenses. APL Photonics, [8] A. She, S. Zhang, S. Shian, D. R. Clarke and F. Capasso. Adaptive metalenses with simultaneous electrical control of focal length, astigmatism, and shift. Science Advances, Vol. 4, no. 08. [9] W. T. Chen, A. Y. Zhu, V. Sanjeev, M. Khorasaninejad, Z. Shi, E. Lee and F. Capasso. A broadband achromatic metalens for focusing and imaging in the visible Nature Nanotechnology vol. 3, pa [0] C. E. Gutiérrez, L. Pallucchini, and E. Stachura, General refraction problems with phase discontinuities on nonflat metasurfaces, JOSA A 347, [] J. Urbas, On the second boundary value problem for equations of Monge-Ampére type, J. Reine Angew. Math , 5-4. [] X. J. Wang, On the design of a reflector antenna, Inverse Problems 996, [3] C. E. Gutiérrez and Qingbo Huang, The refractor problem in reshaping light beams, Arch. Rational Mech. Anal , no., [4] C. E. Gutiérrez. Refraction problems in geometric optics. In Lecture Notes in Mathematics, vol. 087, pages Springer-Verlag, 04. [5] C. E. Gutiérrez and A. Sabra. Aspherical lens design and imaging. SIAM J. Imaging Sci., 9:386 4, 06. [6] C. E. Gutiérrez and Henok Mawi, The far field refractor with loss of energy, Nonlinear Analysis: Theory, Methods & Applications 8 03, 46. [7] C. E. Gutiérrez and Ahmad Sabra, The reflector problem and the inverse square law, Nonlinear Analysis: Theory, Methods & Applications 96 04, [8] C. E. Gutiérrez and Federico Tournier, The parallel refractor, Development in Mathematics 8 03, [9] Roberto De Leo, C. E. Gutiérrez and Henok Mawi, On the Numerical Solution of the Far Field Refractor Problem, Nonlinear Analysis: Theory, Methods & Applications 57 07, [30] Cristian E. Gutiérrez. Refraction problems in geometric optics. In Lecture Notes in Mathematics, volume 087, pages Springer, 04. [3] C. E. Gutiérrez and A. Sabra. Freeform lens design for scattering data with general radiant fields. Archive for Rational Mechanics and Analysis, pp , May 08, Volume 8. [3] C. E. Gutiérrez. The Monge-Ampère equation. Birhäuser, Boston, 06, nd Edition. [33] A. Figalli. The Monge-Ampère Equation and Its Applications. EMS Zurich Lectures in Advanced Mathematics, European Mathematical Society, vol., 07.

20 0 CRISTIAN E. GUTIÉRREZ AND LUCA PALLUCCHINI [34] S. R. Biswas, C. E. Gutiérrez, A. Nemilentsau, In-Ho Lee, Sang-Hyun Oh, P. Avouris, and Tony Low. Tunable Graphene Metasurface Reflectarray for Cloaking, Illusion, and Focusing. Physical Review Applied, 0340, volume 9, 08. [35] Rudolf Karl Luneburg. Mathematical theory of optics. Univ of California Press, 964. Department of Mathematics, Temple University, Philadelphia, PA, 9 address: gutierre@temple.edu, luca.pallucchini@temple.edu