Skip to Main Content

Freeform Optics: Optimal Transport, Minkowski Method, and Monge–Ampère-Type Equations

Henok Mawi

Communicated by Notices Associate Editor Reza Malek-Madani

Article cover

1. Description and Background

A freeform optical surface, simply stated, refers to a surface whose shape lacks rotational symmetry. The use of such surfaces allows generation of complex, compact, and highly efficient imaging systems. Ever since lenses without symmetries were used in World War I in periscopes, the engineering and design of freeform optical surfaces have gone through remarkable evolution, with applications in a wide range of areas, including medical devices, clean energy technology, military surveillance equipments, mobile displays, remote sensing, and several other areas of imaging and nonimaging optics that can benefit from distributing light from a source to a target in a controlled fashion using spatially and energy-efficient systems. See 17 and the references therein.

Mathematically, the design of freeform optical surfaces is an inverse problem related to optimal transportation theory and leads to a class of nonlinear partial differential equations (PDE) called generated Jacobian equations for which the Monge–Ampère equation is a prototype. The modeling of the problem is based on the systematic application of the laws of reflection and refraction in geometric optics along with energy conservation principles.

For completeness we recap the laws of geometric optics. Consider a reflecting or refracting surface and suppose a unit vector in the direction of is incident on the surface at a point where the normal to is given by the unit vector . Reflection and refraction happen in tandem. (Fig. 1). However, for simplicity of the models we treat them separately. If is a unit vector in the direction of the reflected ray, the law of reflection (angle of incidence is equal to angle of reflection) in vector form says:

If medium and medium are two homogeneous isoptropic media with respective refractive indices and and is the direction of the refracted ray into medium , the law of refraction (Snell’s law), which states where is the angle between and (angle of incidence) and is the angle between and (angle of refraction) in vector form, states that:

where is the relative refractive index, and points into medium .

Figure 1.

Laws of reflection and refraction.

Graphic without alt text

When medium I is optically denser than medium II or equivalently , the refracted ray bends away from the normal. As a result, there is a critical value of the angle of incidence for which there is no value of unless . Thus for we impose the condition that

Similar physical restriction applies for .

By using equation 1.1 or 1.2 we can define a map that sends points in (incident rays) to points in (reflected rays if we use 1.1 or refracted rays if we use 1.2) according to laws of geometric optics.

Generally, in an optical surface design problem the data consists of two bounded domains , contained in and two probability measures defined on and respectively, satisfying the energy conservation statement

and correspond to the directions of incident light and the directions of light after reflection or refraction respectively. Likewise and are the intensity of light from the source and a desired illumination intensity on the target , respectively. We regard as a source pair and as a target pair.

The objective is then to find an optical surface (lens/mirror) such that

for all .

In practice, is a reflecting mirror or refracting lens and for , represents the set of all directions from the source that will, after reflection or refraction off of the lens , propagate in the direction of . We use the notation for . The map is called the ray-tracing map of .

The physical meaning of equation 1.5 is that the prescribed intensity of reflected or refracted light on the set , which is , is equal to the intensity of light that comes from the source, which is . Depending on whether is a reflector or a refractor, the design problem is called a reflector/refractor problem.

The models used to obtain numerical and analytical solutions for reflector/refractor problems are developed by using one of the following three approaches. One of the approaches is to use variational method where the problem is cast as an optimal mass transportation problem from to for an appropriate cost function and use Kantorovich duality, 20. Another method is to derive the second order nonlinear PDE of Monge–Ampère-type satisfied by the surface defining function from the energy conservation relation 1.5 between the light energy emitted by the source and the light energy received by the target and analyze the PDE, 11. A third approach is to exploit mainly the inherent geometric features of the problem and use a classical approach of Minkowski and Aleksandrov which was used in solving the Minkowski problem 16 of finding a closed convex surface whose Gaussian curvature is a given function of the exterior unit normal.

In the discussion below, we will focus on refractor problems and briefly describe the aforementioned approaches by using prototype refractor problems. We will state the problems explicitly and show how the available methods are used to solve the corresponding design problems. We will also mention problems that are worth studying in the future. We point out to the reader that unless mentioned otherwise, we consider .

2. Optimal Transport: The Far Field Refractor Problem

Several freeform lens design problems have been successfully modeled by using optimal transport framework on the sphere. See 41020. The models have also been adopted by the optics community. Among these problems is the far field refractor problem.

2.1. Statement of the far field refractor problem with point source

In this problem, we are given two domains contained in the unit sphere of (in practice ), ; a punctual source of light located at the origin , surrounded by media I, from which a monochromatic ray of light is issued in each direction with intensity density function given by for a.e. on and a prescribed intensity density distribution which satisfies

The objective is to find a refracting surface parametrized as , between medium and medium such that all rays emitted from the point with directions are transported via refraction by the surface into media with directions in in such a way that the prescribed illumination pattern is achieved on .

In particular, if we let and where is the surface measure on , this means and satisfy the mass balance condition and a weak solution to the far field refractor problem is defined as an interface so that the associated map satisfies

We recall for that the measure is the push forward measure. Because of 1.3, we assume for and .

By using Snell’s law 1.2 it can be shown that the semi-ellipsoid of revolution (ellipse in 2D case) with one focus at , given in polar representation as (10)

for some and , has a uniform refracting property. This means that if is an interface between media and , all rays issued from in a direction with will be refracted in the direction after hitting . Motivated by this, a solution for the far field refractor problem is sought among surfaces which are supported from outside by semi-ellipsoids at each point. That is, for each point there exists a semi-ellipsoid with such that with equality at . We now discuss how optimal transport technique is used to find such a solution. First, a brief review of optimal transport.

2.2. Review of optimal transport problems

Let and be compact subsets of and . We refer to as the cost function. Denote by and the set of probability measures on and respectively. Given and a measurable map satisfying is called a transport map from to . Monge’s formulation of the optimal transport problem is to minimize

among all transport maps from to . A minimizer of 2.3 is called an optimal transport map. A generalization due to Kantorovich 18 of Monge’s problem is stated as

among all transport plans where

Here and are projection maps. If is lower semicontinuous and bounded from below, problem 2.4 has a minimizer in . A minimizer of 2.4 is called an optimal plan.

The -transform of a given function is defined as where

Also , the -transform of , is defined by

A function is called -concave if for some . We also define for any the -superdifferential at as

Problem 2.4 is a linear optimization problem with convex constraints and thus, it admits a dual problem of maximizing

over and 18.

One of the basic results in optimal transport theory is that if is continuous and and are compact, then the dual problem 2.5 has a maximizer of the form where the map , which is referred to as Kantorovich potential, is a -concave function. Furthermore, if is single valued for , then induces the optimal plan for 2.4. That is, minimizes 2.4. Consequently, is an optimal transport map from to .

2.3. The far field refractor problem as optimal transport

To transform the far field refractor problem 2.2 to an optimal transport problem, we introduce the logarithmic cost function given by

is continuous and bounded from below since . Thus, if we consider the dual problem 2.5 corresponding to given as

there exists 10 a -concave such that is a maximizer of 2.7. Let be -concave for the cost function . If we set

then for each there exists and such that

That is, at each point the surface is supported from outside by a semi-ellipsoid for some and .

If is smooth, then and have the same normals at . By the refraction property of a ray emanating from in a direction will be refracted off in the direction by the surface . Therefore, . Thus, satisfies, . Also, from 2.9 we get

and hence

Furthermore if is given by 2.6, for any -concave the -superdifferential mapping is single valued for . Thus if is a maximizer of 2.7, is a measure preserving map from to . Therefore -a.e., and more importantly proving that will be a solution of the far field refractor problem. From this we deduce that solving the far field refractor problem corresponds to finding a maximizer of the dual problem to the optimal transport problem 2.7. It is further shown in 10 that if and are two solutions of the far field refractor problem, then for some constant C.

The optimal transport approach can also be used in the case of anisotropic materials where optical properties vary according to the direction of propagation of light 6. In this type of media modeling the design problem is more complicated due to issues such as birefringence, a situation where the incident rays may be refracted into two rays. It is also used to develop numerical methods to approximate solutions 4 to freeform design problems. Recently, a more efficient method called entropic regularization is applied to the the reflector problem 2 and it is likely that this approach can be applied to refractor problems.

Finally, it is worth noting that the far field problem has not been studied in its full generality using optimal transport method. In practice, if a ray of light hits an interface between two media which have different optical properties, the ray is partly reflected back and partly transmitted (Fig. 1). It is an open problem whether or not this energy imbalance between the source and target can be modeled as a variational problem by using the optimal transport approach.

3. Minkowski Approach: The Near Field Refractor Problem

The Minkowski problem involves finding a closed convex surface whose Gaussian curvature at a point with exterior normal is given by a continuous positive function 16. To solve this, the approach was first to define a set function on by where is a surface element on the unit sphere, and consider the general problem of finding a convex hypersurface for which is a surface area function. That is, to find a convex hypersurface such that for a given subset of , the area of is equal to where . Then partition the sphere into smaller domains with area and obtain an approximation to . Replacing by gives a discrete problem of existence of convex polyhedron with faces of areas and outward normal . Let be a solution of the discrete problem. As (or the diameter of to ), it is proved that the sequence of surfaces converges to a surface which solves the Minkowski problem. This classical approach of Minkowski is used to obtain iterative methods to solve not only various geometric optics problems related to both reflection 3 and refraction, but also to develop numerical methods to solve PDEs 15. The approach is also generalized to approximate solutions to semi-discrete optimal mass transport problems and generated Jacobian equations. See 1 and the reference therein. It is also one of the methods used by the optics community to design freeform optical surfaces 14. Below we exhibit how this approach is used in the near field refractor problem.

3.1. Statement of the near field refractor problem with point source

For this problem we are given , , a hyperplane in with and the origin . For each , a light ray is issued in the direction of from a point source of light located at . The illumination intensity of the source is given by a density function , a.e. A prescribed irradiance distribution on the target set is given by a nonnegative density function . and satisfy the mass balance 2.1. We assume that is surrounded by medium I and is surrounded by medium II. Additional constraints are imposed on and to make the problem physically feasible 9.

The near field refractor problem involves finding an interface between media and , parametrized radially as , that redirects each ray with direction by refraction into so that a prescribed irradiance distribution is obtained on . More precisely we will require

holds for all Borel sets .

Similar to the semi-ellipsoids in the far field problem, there are surfaces that have uniform refracting property in the near field case. The Descartes ovoid refracts all rays to a single point (Fig. 2). Let and . The Descartes ovoid with one focus at the origin is given as

where

If the relative refractive index of the material inside the oval to that of outside is , then by using 1.2 it can be shown that any ray emanating from the origin and having direction with will pass through the point after refracting off of the ovoid .

Figure 2.

Refraction property of oval; , .

Graphic without alt text

The Descartes ovoids will be used as the building blocks of the near field refractors and we look for a solution for the near field refractor problem among surfaces which are supported from above by an ovoid at each point.

It is known that the near field refractor problem can’t be cast as an optimal mass transport problem. However, Minkowski’s approach can be used to solve the problem.

3.2. The near field refractor problem using Minkowski’s approach

Once again let and . The road map is as follows. First partition and approximate the measure on the target by a sequence of discrete measures concentrated at finite points in and that converge to as . Then solve the problem when the target distribution is . Finally let