arxiv: v2 [physics.optics] 9 Aug 2017

Transcription

1 Single freeform surface design for prescribed input wavefront and target irradiance Christoph Bösel, and Herbert Gross, Friedrich-Schiller-Universität Jena, Institute of Applied Phsics, Albert-Einstein-Str. 5, 7745 Jena, German Fraunhofer-Institut für Angewandte Optik und Feinmechanik, Albert-Einstein-Str. 7, 7745 Jena, German arxiv:75.454v [phsics.optics] 9 Aug 7 In beam shaping applications, the minimization of the number of necessar optical elements for the beam shaping process can benefit the compactness of the optical sstem and reduce its cost. The single freeform surface design for input wavefronts, which are neither planar nor spherical, is therefore of interest. In this work, the design of single freeform surfaces for a given zero-étendue source and comple target irradiances is investigated. Hence, not onl collimated input beams or point sources are assumed. Instead, a predefined input ra direction vector field and irradiance distribution on a source plane, which has to be redistributed b a single freeform surface to give the predefined target irradiance, is considered. To solve this design problem, a partial differential equation (PDE) or PDE sstem, respectivel, for the unknown surface and its corresponding ra mapping is derived from energ conservation and the ra-tracing equations. In contrast to former PDE formulations of the single freeform design problem, the derived PDE of Monge-Ampère tpe is formulated for general zero-étendue sources in cartesian coordinates. The PDE sstem is discretized with finite differences and the resulting nonlinear equation sstem solved b a root-finding algorithm. The basis of the efficient solution of the PDE sstem builds the introduction of an initial iterate constuction approach for a given input direction vector field, which uses optimal mass transport with a quadratic cost function. After a detailed description of the numerical algorithm, the efficienc of the design method is demonstrated b appling it to several design eamples. This includes the redistribution of a collimated input beam beond the paraial approimation, the shaping of point source radiation and the shaping of an astigmatic input wavefront into a comple target irradiance distribution. PACS numbers: 4.5.-i, 4.5.Eq. I. INTRODUCTION The development of numerical design algorithms for the calculation of freeform surfaces for irradiance control in the geometrical optics limit is of great interest in current research. The task of these methods is thereb the calculation of freeform surfaces, which redistribute a given light source into a required target irradiance pattern. In particular the calculation of continuous freeform surfaces, which is considerd here, is of importance for manufacturabilit and the reduction of possible diffraction effects. Especiall the irradiance control with continuous single freeform surfaces for input beams with planar or nonspherical wavefronts lacks of detailed discussion in the literature. Such design algorithms can potentiall reduce the number of necessar optical elements or increase the efficienc in practical applications and are investigated in this work. In principle, the freeform design for shaping of ideal sources can be divided into three major groups. The first group consists of design approaches using Oliker s supporting quadrics method (SQM) [ 4], which is arguabl the theoreticall best justified, most reliable freeform design method and usable in comple optical sstems []. As it was shown b Oliker [], the geometric properties of quadric surfaces can be used effectivel for the freeform design process. The shaping of e.g. a point source with a christoph.boesel@uni-jena.de freeform mirror can be done b placing the point source in the common focal point of a number of ellipsoids. B changing the focal parameters of each ellipsoid appropriatel, a freeform mirror, which gives the required irradiance distribution on the target area, can be designed b a unification of the ellipsoids. Due to the high computational effort, necessar for high resolution target irradiances and large number of quadrics, the major difficult, considering the SQM, is the development of efficient algorithms. The second group of design approaches consists of methods, which are modelling the freeform design problem b a highl nonlinear partial differential equation (PDE) of Monge-Ampère tpe for point sources [5 7] and collimated beams [8, 9]. The main problem is thereb the development of numerical algorithms to solve the PDE. An advantage of the PDE methods over the SQM is the lower computational effort, enabling them to handle high resolution target irradiances. The major drawback of the PDE methods is the lack of theoretical justification and reliabilit as pointed out b Oliker [4]. The third group are ra-mapping or source-target mapping methods, respectivel [ 5]. These are based on the calculation of a ra mapping between the source and the target irradiance distribution an the subsequent construction of the freeform surface from the mapping. The major difficult of these methods is the finding of an integrable mapping, enabling the design of a freeform surface, which is on the one hand continuous and on the other hand maps the source onto the target distribution. Some of the ra-mapping methods [ ] eplore a mapping

2 calculated from optimal mass transport (OMT) theor with a quadratic cost function. This special mapping solves the collimated beam shaping problem in a paraial approimation [], but is in general not integrable. As we will demonstrate in this work, the quadratic cost function mapping is nevertheless of great usefulness for the freeform design process for more general optical configurations. In this work, we want to generalize the design approach for collimated beam shaping with continuous single freeform surfaces presented in Ref. [] to general input wavefronts of ideal sources, whereb, contrar to former PDE formulations, parallel input beams or point sources are not assumed. Instead, we derive a sstem of PDE s and a Monge-Ampère equation, respectivel, for a predefined input direction vector field and present a detailed numerical algorithm to solve the corresponding problem. As indicated in the previous paragraph, the basis of the numerical algorithm builds thereb a ra mapping, derived from OMT with a quadratic cost function, which will serve as an initial iterate for solving the derived PDE sstem. As we will demonstrate eplicitel, the design algorithm etends the work of Ref. [] firstl to collimated beam shaping in the nonparaial regime, secondl to the irradiance control of point sources, and thirdl, to other input wavefront shapes. To introduce the design algorithm, this work is structured as follows. In section II and III, a PDE sstem is derived from energ conservation and the ra-tracing equations for single freeform lens and mirror surfaces. For the numerical solution of the PDE sstem, a method for the calculation of an appropriate initial iterate is introduced in section IV. In section V, the numerical algorithm for determining the desired freeform surface is presented and in section VI several design eamples, illustrating the effectiveness and usefulness of the design algorithm, are discussed. II. GEOMETRY In the following, the problem of mapping an incoming beam, defined b the ra direction vector field ŝ = (ŝ ) (ŝ ), () (ŝ ) z and the irradiance distribution I S (, ) on the plane z = z, onto a predefined output irradiance distribution I T (, ) on the target plane z = z T is considered (see Fig. ). The goal is thereb the calculation of a continuous reflecting or refracting surface, which enables the mapping of the source onto the target distribution. In a geometrical optics framework, this problem is full governed b two equations. This is on the one hand the Jacobian equation I S (, ) I T (, ) Ω S Ω T FIG.. The input I S(, ) on the plane z = z, the output irradiance I T (, ) on the plane z = z T and the input directions ŝ are given. The input vector field is redirected b the freeform surface z(, ) to the target points u(, ) according to the law of refraction/reflection to give the desired ouput distribution. The coordinates ( S, S) represent the intersection points of the input vectors with z(, ). det( u(, ))I T (u(, )) = I S (, ), () describing the relation of the input and output irradiance b the ra mapping or coordinate transformation u(, ), respectivel, and on the other hand the ra tracing equations for refractive and reflective surfaces { ŝ = nŝ + n ˆn ŝ + } n [ (ˆn ŝ ) ] ˆn, ŝ = ŝ (ˆn ŝ )ˆn () with the surface unit normal vector field ˆn. Equation () relates the input coordinates (, ) and the target coordinates u(, ) b the unknown surface z(, ). Unit vector fields are hereb denoted b the hat. To determine the unknown surface, we want to derive a PDE sstem for the mapping u(, ) and the surface z(, ). This can be done b epressing the considered geometr, according to Fig., through the vector fields s = S S, s = u (, ) u (, ), z( S, S ) z z T z n = S z( S, S ) S z( S, S ), (4) and b calculating the intersection points u(, ) with the target plane z = z T from the relation

3 s + λ ŝ = s λ = z T z( S, S ) (ŝ ) z. (5) Hereb S, S and z( S, S ) are the intersection coordinates of the input vector field s with the surface z(, ). From Eq. (5) we then directl get u (, ) = λ (ŝ ) + (s ), u (, ) = λ (ŝ ) + (s ), (6) which define a relation between the surface z(, ), the surface gradient and the target coordinates. For the case of collimated input beams, analogues equations were derived b Wu et al. [8] and Oliker [6]. Together with Eq. (), Eq. (6) builds a PDE sstem for the unknown functions u (, ), u (, ) and z(, ). According to Eq. (4), the remaining difficult is thereb the dependenc of Eq. (6) on the surface intersection points ( S, S ), which are determined b the surface z(, ) and therefore unknown. A similar problem arising in their work on the reconstruction of an optical surface from eperimental ra data [7], was solved b Rubinstein and Wolansk using the geometrical relation S z( S, S ) z = (ŝ ) (ŝ ) z S = (ŝ ) (ŝ ) z [z( S, S ) z ] +, S = (ŝ ) S = (ŝ ) z( S, S ) z (ŝ ) z (7) between the given ra directions defined in Eq. () and the vector field s in Eq. (4). B appling this coordinate transformation to the coordinates S, S and the gradient S z( S, S ) = ( S z( S, S ), S z( S, S )) occuring in Eq. (4), Eq. (6) can be epressed b the initial coordinates (, ) on the plane z = z. After using the chain rule it follows Monge-Ampère equation in Eq. (), we will therefore S z( S, S ) = ( S ) [ z( S, S )] ( S )[ z( S, S )] omit the discretization of second order derivatives. In, addition to that, the boundar conditions can be directl ( S ) ( S ) ( S ) ( S ) S z( S, S ) = ( controlled through the mapping u(, ) b requiring S )[ z( S, S )] + ( S )[ z( S, S )]. that boundar points of the support of I S (, ) are ( S ) ( S ) ( S ) ( S ) (8) The remaining derivatives of S and S can be evaluated with Eq. (7). With these transformations and the definition z S (, ) z( S, S ), (9) Eq. (6) is full epressed in terms of the -coordinates. Hence, the Jacobian equation () and Eq. (6) build a sstem of PDEs for the unknown functions u(, ) and z S (, ). After solving the equation sstem for z S (, ), the coordinates ( S, S ) can than be calculated from Eq. (7), which gives the surface z(, ) on the scattered S - S grid. III. PDE SYSTEM AND MONGE-AMPÈRE EQUATION In the previous section it was shown that the single freeform design problem for zero-étendue sources is governed b the PDE sstem [( u )( u ) ( u )( u )] I T (u, u ) = I S (, ), u (, ) = [z T z S (, )] (ŝ ) (ŝ ) z + [z S (, ) z ] (ŝ ) (ŝ ) z, u (, ) = [z T z S (, )] (ŝ ) (ŝ ) z + [z S (, ) z ] (ŝ ) (ŝ ) z, () whereb ŝ is the predefined input ra direction vector field and ŝ is defined according to section II in terms of z S (, ), z S (, ) and z S (, ). For a given source and target distribution I S (, ) and I T (, ), and a given input ra direction vector field ŝ the equation sstem in Eq. () has to be solved simultaneousl for the surface z S (, ) and the ra mapping u(, ). Equation () can be further simplified b plugging u (, ) and u (, ) into the Jacobian equation to obtain a nonlinear second order PDE of Monge-Ampère tpe for the function z S (, ): [z( S, S ) z ] + A, z S + B z S + C z S + [ ( z S )( z S ) ( z S ) ] = D. (ŝ ) z () Hereb, A, B, C, D are functions of,, z S (, ), z S (, ) and z S (, ), which can be calculated straightforeward. Compared to the Eq. sstem (), the Eq. () has numericall the advantage of a lesser number of design variables. Despite this advantage, the numerical method presented in the following will nevertheless be based on the description in Eq. (). In contrast to the mapped onto boundar points of the support of I T (, ). This is contrar to the Monge-Ampère equation for which the same goal has to be accomplished through the surface gradients z S (, ) at the boundar of the surface, which leads to the dependenc of the boundar conditions on the input field ŝ. Furthermore, the properties of the Jacobian equation and the mapping equations in Eq. () can be controlled separatel. IV. INITIAL ITERATE The construction of freeform surfaces for comple target irradiances from nonlinear PDEs, like Eqs. () or () is usuall done b discretizatizing the PDEs and

4 4 solving the resulting nonlinear equation sstem b an iterative process. A high resolution target irradiance distribution and therefore large number of unknowns requires in general an appropriate initial iterate to ensure (fast) convergence. Especiall for discretization schemes, which are not provabl convergent, those initial iterates are of great importance. For the published PDE methods for single freeform surface design an inital iterate producing a uniform illumation on the target plane or other simple surface shapes are usuall chosen. Here, we present a different method to construct an initial iterate for a given input direction vector field ŝ. It is based on the OMT with the quadratic cost function and inspired b Ref. [7]. Following the argumentation from Ref. [7], an unknown surface can be eperimentall reconstructed from the known input ra directions and measured target plane intersection points of the ras send out from (, ). Since the target points are not given for the freeform design problem considered here, the approach from Ref. [7] can not be directl applied. Nevertheless, as we will demonstrate in the following, the approach is of great usefullness to generate an initial iterate for the freeform surface design. () as an initial mapping for solving Eq. (). As in the previous publications [, 9], we will make use of the algorithm b Sulman et al. [] for the eamples in section VI for demonstrational purposes. B. Initial surface After its calculation, the mapping u (, ) can be used to construct the initial surface zs (, ) as demonstrated in the following. The resulting initial surface construction equations are thereb a generalization of the surface construction equations used in Ref. [] to noncollimated beams. The derivation can be done b using the law of refraction/reflection n( S, S ) (n( S, S )) z = n ŝ n ŝ (n ŝ n ŝ ) z = S (z z S (, )) () and appling again the coordinate transformation defined in Eq. (7) to invert the gradient b using A. Initial mapping For the single freeform surface design process a frequentl used mapping u(, ) is calculated from OMT b considering Eq. () together with the quadratic cost function condition d(i S, I T ) = inf u M u() I S ()d, () which we will refer to in the following b u (, ). M in Eq. () denotes the set of energ preserving mappings. This is arguabl the most studied OMT problem in numerical mathematics and provides us therefore with efficient algorithms for the mapping calculation. As it was demonstrated in several publications, the corresponding mapping u (, ) in general will not generate a continuous surface z(, ) fulfilling the required design goal of mapping the given input I S (, ) onto the desired output irradiance I T (, ). However, due to its properties, the mapping u (, ) appears to be a good candidate for creating an appropriate initial iterate. According to Eq. (), the quadratic cost function can in principle be thought of as describing a ra propagation in a homogeneous medium [8]. Since a continuous surface can be interpreted as a perturbation of this free propagation, it can be epected that the mapping u (, ) serves as a reasonable initial iterate for solving Eq. (). From a numerical point of view it is also helpful that u (, ) alread fulfills the Jacobian equation in Eq. (), which therefore basicall constrains all numerical corrections to u (, ) to the space of energ preserving mappings. Hence, we will use the mapping defined b Eqs. () and S z S (, ) = n (ŝ ) n (ŝ ) n (ŝ ) z n (ŝ ) z n (, ), S z S (, ) = n (ŝ ) n (ŝ ) n (ŝ ) z n (ŝ ) z n (, ). (4) After using the chain rule, this gives the differential equations [7] zs (, ) = [ ( ) ( )] n (, ) + n (, ) (ŝ) (ŝ) (ŝ ) z + n (, ) (ŝ ) z (zs (, ) z ) n (, ) (ŝ) (ŝ ) z n (, ) (ŝ) (ŝ ) z zs (, ) = [ ( ) ( )] n (, ) + n (, ) (ŝ) (ŝ) (ŝ ) z + n (, ) (ŝ ) z (zs (, ) z ) n (, ) (ŝ) (ŝ ) z n (, ) (ŝ) (ŝ ) z (5) for the initial surface zs (, ). For a given input direction vector field, defined b Eq. (), and given target points u (, ), Eq. (5) can be integrated along an arbitrar path on the support of I S (, ) to generate the initial surface zs (, ). Together with the initial mapping u (, ), the surface zs (, ) builds an initial iterate for solving the PDE sstem (). We want to remark that the initial surface construction approach, presented in this section can also be used for the Monge-Ampère equation in Eq. (). To implement the boundar conditions in that case, the surface boundar normal vectors defined b Eq. (5) have to be chosen in a wa that the boundar of I S (, ) is mapped onto the boundar of I T (, ).

5 5 V. NUMERICS A. Discretization After the calculation of the initial iterate, the PDE sstem () has to be solved. Hence, we have to discretize () and solve the resulting nonlinear equation sstem b a root-finding algorithm. To do that, we follow Ref. [9] b introducing the corrections u(, ) and z S (, ) to u (, ) and z (, ): B. Design algorithm The design algorithm developed in the previous sections can be summarized as follows:. Define the input ra direction vector field ŝ and the input irradiance I S (, ) on z = z and the output irradiance I T (, ) on z = z T.. Calculate the initial map u (, ) between I S (, ) and I T (, ). u(, ) = u (, ) + u(, ), z S (, ) = zs (, ) + z S (, ) (6). Integrate Eq. (5) along an arbitrar path on the support of I S (, ) to get the initial surface zs (, ). and redefine the initial iterates for the solving process to be u(, ) = and z(, ) =. Accordingl, Eq. () becomes a PDE sstem for u(, ) and z S (, ). To solve the PDE sstem, we discretize the functions u (, ), u (, ), z S (, ) and their derivatives in Eq. () b central finite differences ( z S ) [( z S) i;j+ ( z S ) i;j ], ( z S ) [( z S) i+;j ( z S ) i ;j ] (7) for the inner points i; j =,..., N and second order finite differences ( z S ) [( z S) i;j+ 4( z S ) i;j+ + ( z S ) i;j ], ( z S ) [( z S) i+;j 4( z S ) i+;j + ( z S ) i;j ], (8) for the boundar points. Additionall, we appl boundar conditions b requiring that the boundar points of the source distribution I S (, ) are mapped onto the boundar points of the target distribution I T (, ). For the mapping of square supports, which is demonstrated in section 6, the boundar condition reads therefore: u (.5, ) = u (.5, ) =, [.5,.5], u (,.5) = u (,.5) =, [.5,.5]. (9) After the discretization of Eq. (), the resulting nonlinear Eq. sstem can then be solved b a root-find algorithm for the N unknown values ( u ) i;j, ( u ) i;j and ( z S ) i;j. The solution will lead to a discretized surface z S (, ) from which, according to Eq. (7), the scattered grid points ( S, S ) can be calculated. As in Ref. [9], we use MATLAB s fsolve() function with the trustregion-reflective solver for the solution of the nonlinear Eq. sstem (). 4. Discretize the PDE sstem in Eq. () to get a nonlinear equation sstem. 5. Appl the boundar conditions and solve the nonlinear equation sstem with u (, ) and zs (, ) as the initial iterate to get u(, ) and z S (, ). As mentioned above, for a general input vector field ŝ, the presented design algorithm gives us the freeform surface z(, ) on the scattered grid points ( S, S ). The qualit of the surface z S (, ) is thereb alread characterized b Eq. (), which is equivalent to a ra-tracing calculation. For manufacturing processes it might still be useful to evaluate the surface qualit in a ra-tracing simulation. We therefore interpolate the scattered data points z( S, S ) onto a regular grid and import it to a ra-tracing software. The interpolation can be done b scattered data interpolation with e.g. MATLAB 5b s griddata() function. However, in our eperience this function has two major drawbacks for the interpolation of freeform surfaces. This is its inabilit (for large data sets) to etrapolate outside of the conve hull of the grid points, which leads to inaccuracies at the boundar of the freeform surface in form of zig-zag patterns at the target irradiance boundar. Furthermore, the interpolation with griddata() leads to oscillations in the freeform surface, depending on the stepsize of the interpolation grid and the interpolation method. To overcome this limitation, we therefore use the RegularizeDataD() function from the MATLAB file echange website []. This function provides the possibilit to choose how smooth the interpolated data is fitted to the scattered data. In our eamples, we made good eperience with a small smootheness factor for low noise data. Additionall, RegularizeDataD() etrapolates outside the conve hull of the scattered data grid, which ma also lead to slight inaccuracies at the boundar of the freeform surface, which lower the energ efficienc, as we will point out in the eample section VI. Since the initial mapping u (, ) defined b Eqs. () and () is in general not integrable, the integration of Eq. (5) will be path dependent and therefore lead to different initial surfaces zs (, ) for different integration

6 6 paths. Nevertheless, in our eperience the path dependenc has no significant influence on the suitabilit for using zs (, ) as an initial iterate. b a stepwise reduction of the target plane distance z T. However, in our eperience, this stepwise reduction is in a practical calculation not necessar to ensure numerical convergence. VI. EXAMPLES In this section, we want to demonstrate basicall three things. Firstl, the design algorithm presented in section V etends the work of Ref. [] on collimated beam shaping with single freeform surfaces to the nonparaial regime. Secondl, the design algorithm can be applied to the design of single freeform surfaces for shaping of point light sources, which is frequentl discussed in literature. And thirdl, we want to demonstrate the applicabilit to a nonstandard design eample b dealing with an input wavefront, which is neither planar nor spherical. This is done in the following three subsections. All the eamples in this sections were calculated on an Intel Core i at.4ghz with 6GB RAM. To make the ra-tracing results comparable, we used for each of the eamples the target distribution lake (see Fig. (b)) with a resolution of 5 5 piels. For the calculation of the initial surface b Eq. (5), we used MATLAB s ode45 solver with tolerances of, which is in our eperience necessar for a legitimate comparison of the initial and final surface in a ra-tracing simulation. For the root finding of Eq. () b fsolve() tolerances of were used. For the eamples with lenses refractive indices of n =.5 were used. The ra-tracing simulations for the following eamples are done with 6 ras and detectors with 5 5 piels b using a self-programmed MATLAB ra-tracing toolbo. Since the design eamples are calculated in a geometrical optics framework, we will use arbitrar units as length units and omit the eplicit declaration of lengths b a.u.. A. Collimated beam shaping beond the paraial approimation In Ref. [] it was shown that the single freeform surface design problem for collimated input beams ŝ = () can be solved b the quadratic cost function ra mapping, defined b Eqs. () and (), in the paraial regime. Here, we demonstrate that the presented algorithm from section V directl etends the results from Ref. [] to nonparaial calculations. Since the quadratic cost function mapping solves the single freeform surface design problem eactl for infinite distances between the freeform surface and the target plane[], it can be argued analogousl to Ref. [9] that the numerical convergence of the algorithm can theoreticall be ensured FIG.. (a) Sstem laout. The incoming collimated Gaussian beam is reflected b the freeform mirror to give the required target distribution lake at the predefined position. (b) Freeform mirror. From an application point of view, the etension to the nonparaial regime is especiall important for optical configurations in which the distance between the target plane and the freeform surface is small, the target area is strongl shifted or scaled, or if there are caustics between the target plane and the freeform surface. In the following, we therefore demonstrate the utilit of the presented algorithm b designing a freeform mirror, which shifts and scales the target area compared to the position and size of the input beam. The corresponding sstem laout, considered in this section, can be seen in Fig., showing a target area with a side length, which is times as large as the input beams side length of and shifted b 5. The distance between the freeform mirror

7 7 and the target plane z T is approimatel (a) (b) Initial surface (a) Final surface (b) FIG.. (a) Normalized Gaussian input irradiance distribution I S(, ) and (b) normalized output irradiance distribution I T (, ) lake. The output irradiance in (b) has the side length of, but will be scaled and shifted for comparison with the ra tracing simulations according to the laouts of each eample. I S(, ) and the input ra direction vector field in Eq. () at z = z and I T (, ) at z = z T build the input of the design algorithm presented in section V. For the irradiance distribution, we choose the Gaussian distribution, which is represented in Fig. (a) and the test image lake as the required target distribution (Fig. (b)). According to section V, we first calculate the initial mapping between the input and output irradiances (sec), then the initial surface from Eq. (5) (68sec) and finall solve the PDE sstem in Eq. () (8sec). After the construction of the initial and final freeform mirror, we evaluate their qualit b a ra-tracing simulation. The corresponding results and comparisons with the algorithm presented in Ref. [], which is equal to the initial iterate of the root finding, are shown in Fig. 4. Additionall, the rms and the peakto-valle of I T (, ) = I T (, ) I T,RT (, ), and the correlation coefficient corr IT between the predefined irradiance I T (, ) and the irradiance from the ra tracing I T,RT (, ), and the energ efficienc η are given in Table I. The results thereb show clearl the improvement of the presented design algorithm compared to Ref. [] and the mapping from the quadratic cost function (), respectivel. TABLE I. Comparison of I T for eample lake with a collimated input beam. z (, ) z(, ) rms IT pv IT corr IT η 85.66% 99.98% B. Point source Besides the collimated beam shaping with a single freeform surface, we want to appl the design algorithm (c) (d) FIG. 4. Target irradiance distribution I T (, ) from the ra tracing simulation produced b (a) the initial surface and (b) the optimized surface. Difference between the predefined distribution I T (, ) and the ra tracing simulation I T,RT (, ) for (c) the initial surface and (d) the optimized surface. According to (c) the illumation pattern from the initial surface is strongl distorted, whereas the optimized surface produces onl slight deviations from the predefined pattern. to the other standard configuration considered in literature, which is the mapping of a point source onto a required irradiance distribution on a target plane z T. In the case of single freeform lens design, the point source is assumed to be in the medium or the ras emitted b the point source first hit a spherical surface, respectivel. The ra directions are therefore first unaltered and then redirected b the freeform surface to give the predefined target distribution I T (, ). Since we want to appl the algorithm from section V, we first have to project the point source onto the input plane z = z, giving the distribution I S (, ). Subsequentl, the mapping u (, ) between the projected distribution I S (, ) on z = z and the target distribution I T (, ) on z = z T can be calculated. For the point source, we choose a lambertian distribution with a maimum opening angle of degrees, which refers to the angle between the z-ais and the corner points of the quadratic aperture and leads to the irradiance distribution I S (, ) presented in Fig. 5(a). The input vector field ŝ on z = z is given b ŝ = + + z z -. () The placement of the plane z = z relative to the point source does not matter. Depending on the value of z onl the mapping u (, ) has to be scaled b a constant

8 (a ).5 -. (b ). FIG. 5. (a) Projection of the lambertian source onto the input plane z = z, which serves as the irradiance distribution IS (, ) for the calculation of u (, ) b OMT and together with the input ra direction vector field in Eq. () and IT (, ) as the input of the design algorithm. (b) After solving the PDE sstem in Eq. () for the surface zs (, ), the scattered grid points (S, S ) can be calculated from Eq. (7), which gives z(s, S ). This function is then interpolated and etrapolated onto a regular grid. factor to create the predefined target distribution area size, which we choose to have a side length of 5. For the distance between the point source and the freeform surface we choose z(, ) =, which enters the design process as an integration constant in Eq. (5), and for the distance between the point source and the target plane, we choose zt =. The corresponding sstem laout can be seen in Fig. 6. The calculation of the initial mapping u (, ) took approimatel 78sec, the initial surface zs (, ) 8sec and the final surface zs (, ) about sec. After the calculation of the freeform surface zs (, ) b the presented design algorithm an additional complication, compared to the collimated beam shaping, arises since the freeform surface is calculated on the nonuniform S -S -grid. Fig. 5(b) shows ever 4th point of the scattered grid points. The outermost points define the aperture of the freeform lens. As described in section V, the function zs (, ) = z(s, S ) is interpolated and etrapolated onto a square, uniform grid b RegularizeDataD() to give z(, ). For a fine boundar interpolation, we choose an interpolation grid of 5 5 points, whereb the minimum and maimum and values correspond to the minimum and maimum S and S values (see Fig. 5(b)). The function values are then imported to a ra-tracing software. In Fig. 7 and Table II the results of the ra tracing are shown. From Fig. 7(a) and 7(c) it can be seen that the quadratic cost function in Eq. () provides a reasonable initial iterate for the root-finding process, as argued in section IV. Figure 7(b) and 7(c) and Table II show results, which are comparable to the eample from the previous subsection, ecept for the slightl lowered efficienc η due to the interpolation of zs (, ) at the surface boundar. FIG. 6. (a) Sstem Laout. The point source at z = with a lambertian intensit profile is redistributed b the freeform lens to give the required target irradiance at z =. (b) Freeform lens surface. Since RegularizeDataD() etrapolates the surface to a square area, the corner areas are outside of the maimum opening angle of degrees and therefore irrelevant and not considered in the ra-tracing simulation. TABLE II. Comparison of IT for eample lake with a point source. z (, ) z(, ) rms IT pv IT corrit η 99.6% 99.55% From the viewpoint of manufacturing, it is usuall preferred to have a plane instead of a spherical surface as the entrance surface of the lens. In the following section it is demonstrated that the design algorithm can be applied directl to these configurations.

9 9 Initial surface Final surface (a) - (b) (c) - - (d) FIG. 7. Target irradiance distribution IT (, ) from the ra tracing simulation produced b the (a) initial surface and (b) optimized surface. Difference between the predefined distribution IT (, ) and the ra tracing simulation IT,RT (, ) for (c) the initial surface and (d) the optimized surface. Astigmatic wavefront FIG. 8. (a) Sstem Laout. The incoming collimated gaussian beam is refracted b a predefined surface to give an astigmatic input wavefront. After another refraction b the plane entrance surface the ras are refracted b the calculated freeform surface to give the required target distribution. (b) Freeform lens surface In this section, we want to eplore an eample with a nonstandard input vector field s. We therefore predefine a surface zpre (, ) = 9, deforming a collimated beam with a gaussian distribution (Fig. 9(a)) into a beam with an astigmatic wavefront. The corresponding sstem laout can be seen in Fig. 8. As in the previous eamples, we need the input vector field s and the irradiance distribution on an intermediate plane z = z. This plane is chosen to be z = 5 and therefore after the refraction of the astigmatic beam b the plane entrance surface (z = in Fig. 8(a)), and the aperture at z = z is a square with a side length of. The vector field s and the source distribution IS (, ) on z = z can then be calculated numericall, e.g. from the law of refraction and the Jacobian equation (). IS (, ) is shown in Fig. 9(b) and serves as the input for the calculation of u (, ) (47sec). For the calculation of the intial surface with Eq. (5), we have to fi one point of the freeform surface relative to the plane z = z, which we choose to be z(, ) = z + and the side length of the target area is chosen to be. The calculation of zs (, ) took 45sec and the root finding of Eq. () about 5sec. Since the input vector field is noncollimated, the surface z(, ) is again calculated on a scattered grid and has to be interpolated for the ra-tracing simulation. The results of the ra tracing are shown in Fig. and Table III. Also for this eample, the results show that the initial map is appropriate for the root finding of Eq. () and C (a) (b) FIG. 9. (a) Input gaussian beam before the refraction b the predefined astigmatic surface. (b) Source irradiance distribution IS (, ) at z = 5 (Fig. 8(a)) after refraction b the plane surface. IS (, ) is used as the input distribution for the calculation of u (, ) with optimal mass transport and together with the input ra direction vector field s and IT (, ) as the input of the design algorithm. that b the presented design algorithm, single freeform surfaces for comple target irradiances and nonplanar or

10 nonspherical input wavefronts can be calculated. The quantitative differences between the ra-tracing results for the eamples are thereb mainl due to interpolation and etrapolation effects. It is also a priori not clear, if for the particular predefined sstem parameters of each eample an eact solution to Eq. () eists. This might also lead to a degradation of the qualiti of the ra-tracing results, since an accurate root-finding of the respective nonlinear equation sstems might not be possible. In contrast to the first two eamples, additional errors are introduced in the third design eample b the numerical calculation of I S (, ) and the input direction vector field. - Initial surface - (a) - - (c) Final surface - (b) - - (d) -5-5 FIG.. Target irradiance distribution I T (, ) from the ra tracing simulation produced b the (a) initial surface and (b) optimized surface. Difference between the predefined distribution I T (, ) and the ra tracing simulation I T,RT (, ) for (c) the initial surface and (d) the optimized surface. TABLE III. Comparison of I T for eample lake with an astigmatic input wavefront z (, ) z(, ) rms IT pv IT corr IT η 99.% 99.% VII. CONCLUSION A single freeform surface design algorithm for prescribed input wavefronts and comple irradiance distributions was introduced. This was done b deriving the PDE sstem () and Monge-Ampère equation (), which governs the single freeform design process for ideal - sources in the geometrical optics limit, from the energ conservation and ra tracing equations in cartesian coordinates. After introducing a possible wa to calculate an appropriate initial iterate in section IV, the PDE sstem is then solved for the unknown functions u(, ) and z(, ) b a finite difference discretization and a root finding of the resulting nonlinear equation sstem. The design eamples in section VI show that the presented design algorithm etends the results of Ref. [] for collimated beam shaping to the nonparaial regime and that the algorithm is applicable to point sources and more general input wavefronts. This might be helpful in applications in which one has to deal with the shaping of e.g. astigmatic wavefronts produced b laser diodes or simplif manufacturing processes. Also the relaation of phsical situations, e.g. for point sources with large opening angles, might be reached b a predefined entrance surface of a freeform lens. Despite the presented results, there are still interesting problems left for future research and improvements, respectivel. For eample, a mathematical rigorous justification of the initial iterate construction method presented in section IV is of interest, since there is no guarantee that the proposed method will alwas give a suitable intial iterate, despite the arguments given in section IV. Also the numerical etension of the presented algorithm to more complicated boundaries of the irradiance distributions would help to increase the number of possible applications or energ efficienc, respectivel. This requires the numerical implementation of the transport boundar conditions, which are realized for irradiance distributions with square boundar shapes b Eq. (9), for more comple boundar shapes and will be investigated in the future. As mentioned in section III, a numerical investigation of the Monge-Ampère equation in Eq. () is of interest. Despite the listed additional complications, the reduced number of design variables can benefit the speed of the root-finding process. Another interesting future goal is the generalization of the presented design algorithm to wavefront and irradiance control with double freeform surfaces b combining it with the results presented in Ref. [9] and the investigation of the manufacturabilit of the calculated freeform surfaces. FUNDING We acknowledge support from the German Federal Ministr of Education and Research (BMBF) under contract number PT69X within the project KoSimO. ACKNOWLEDGMENTS The authors want to thank Norman Worku for providing and helping them with the MatLightTracer and Johannes Stock for comments on the manuscript.

11 [] D. Michaelis, D. Schreiber, and A. Bräuer, Cartesian oval representation of freeform optics in illumination sstems, Opt. Lett. 6(6), 98 9 (). [] V.I. Oliker, Mathematical aspects of design of beam shaping surfaces in geometrical optics, in Trends in Nonlinear Analsis, M. Kirkilionis, S. Kromker, R. Rannacher, and F. Tomi, eds. (Springer-Verlag, ), pp. 9-. [] F.R. Fournier, W.J. Cassarl, and J.P. Rolland, Fast freeform reflector generation using source-target maps, Opt. Epress 8(5), (). [4] V.I. Oliker, Controlling light with freeform multifocal lens designed with supporting quadric method(sqm), Opt. Epress 5(4), A58-A7 (7) [5] H. Ries and J. Muschaweck, Tailored freeform optical surfaces, J. Opt. Soc. Am. A 9(), (). [6] R. Wu, L. Xu, P. Liu, Y. Zhang, Z. Zheng, H. Li, and X. Liu, Freeform illuminaton design: a nonlinear boundar problem for the elliptic Monge-Ampère equation, Opt. Lett. 8(), 9 (). [7] K. Bri, Y. Hafizogullari, and A. Platen, Designing illumination lenses and mirrors b the numerical solution of MongeAmpère equations, J. Opt. Soc. Am. A (), 7-6 (5) [8] R. Wu, P. Liu, Y. Zhang, Z. Zheng, H. Li, and X. Liu, A mathematical model of the single freeform surface design for collimated beam shaping, Opt. Epress (8), (). [9] C.R. Prins. J.H.M. ten Thije Boonkkamp, J. Van Roosmalen, W.L. IJzerman, and T.W. Tukker, A Monge- Ampre-Solver for free-form reflector design, SIAM J. Sci. Comput. 6(), B64 B66 (4). [] V. Oliker, J. Rubinstein, and G. Wolansk, Ra mapping and illumination control, J. Photon. Energ (), 5599 (). [] A. Bruneton, A. Bäuerle, P. Loosen, and R. Wester, High resolution irradiance tailoring using multiple freeform surfaces, Opt. Epress (9), (). [] Z. Feng, L. Huang, G. Jin, and M. Gong, Designing double freeform optical surfaces for controlling both irradiance and wavefront, Opt. Epress (), (). [] C. Bösel and H. Gross, Ra mapping approach for the efficient design of continuous freeform surfaces, Opt. Epress 4(), (6). [4] D. Ma, Z. Feng, and R. Liang, Tailoring freeform illumination optics in a double-pole coordinate sstem, Appl. Opt. 54(9), (5) [5] C. Gannon and R. Liang, Ra mapping with surface information for freeform illumination design, Opt. Epress 5(8), (7). [6] V. Oliker, Differential equations for design of a freeform single lens with prescribed irradiance properties, Opt. Eng. 5(), (4). [7] J. Rubinstein, G. Wolansk, Reconstruction of Optical Surfaces from Ra Data, Optical Review 8(4), 8-8 (). [8] J. Rubinstein, G. Wolansk, A variational principle in optics, J. Opt. Soc. Am. A (), 64-7 (4) [9] C. Bösel and H. Gross, Ra mapping approach in double freeform surface design for collimated beam shaping beond the paraial approimation, Appl. Opt. 56(), (7). [] M.M. Sulman, J.F. Williams, and R.D. Russel, An efficient approach for the numerical solution of Monge- Ampère equation, Appl. Numer. Math. 6(), 98-7 (). [] regularizedatad