|
Tags: Blanking Manual revert |
| (34 intermediate revisions by the same user not shown) |
| Line 1: |
Line 1: |
| == Introduction ==
| |
|
| |
|
| Conventional optical imaging systems are bulky and complex, requiring multiple elements to correct aberrations. Optical metasurfaces, planar structures that are capable of manipulating light at subwavelength scales, are compact alternatives to conventional refractive optical elements. Their miniature volume is suitable for technologies like AR/VR displays and wearables. However, existing metalenses face significant challenges from monochromatic (e.g., coma) and chromatic aberrations, limiting their applicability. Here we present an end-to-end AI-based computational method that parametrizes the profile of metalenses and optimizes it based on customized loss functions. This innovation enables wide-angle imaging with corrected aberrations while retaining a single-layer form factor, overcoming the key limitations of existing metalenses and advancing their potential for miniaturized imaging systems.
| |
|
| |
| == Background ==
| |
|
| |
| === What Is Metasurface? ===
| |
|
| |
| A metasurface is an two-dimensional artificially engineered material composed of subwavelength-scaled patterns. Metasurfaces manipulate electromagnetic waves, e.g. light waves, through specific boundary conditions imposed at their interfaces, unlike conventional materials which manipulate EM waves through its bulk properties, such as refractive index. Metasurface's unique approach enables precise control over wavefronts, allowing for innovative applications such as planar lenses (metalenses) and holograms. Their thin profile and design flexibility are important in integrated photonics and advanced imaging systems.
| |
|
| |
| [[File:Fig-metasurface.png|thumb|center|600 px|Figure 1: Metasurfaces with various nanostructures. (Source: Meinzer 2014)]]
| |
|
| |
| === Metasurface Optic ===
| |
|
| |
| Transfer function of conventional thin refractive lenses has amplitude and phase. Phase change, of an incident plane wave propagating through the lens, will leads to a change in the wave vector. That is, the incident plane wave will change its propagating direction. In a metalens, the phase is induced via the response of nanostructures. By applying local, gradient phase shifts to incoming waves, metasurfaces generalize the conventional laws of reflection and refraction, called generalized Snell's law. In the equation below, <math> \theta_t </math>, <math> \theta_i </math> are the reflection/refraction and incident angle, respectively; <math> n_t </math>, <math> n_i </math> are the refractive index of the corresponding material; <math> \lambda </math> is the wavelength of the light; and <math> \Phi </math> is the phase profile of the metalens.
| |
|
| |
| <math display="block"> n_t \sin \theta_t - n_i \sin \theta_i = \frac{\lambda}{2\pi} \frac{\text{d}\Phi}{\text{d}x} </math>
| |
|
| |
| For conventional laws of reflection and refraction, the term on the right hand side of the equation is 0. The generalized Snell's law implies that we can control the refraction angle by designing the phase profile of the metalens.
| |
|
| |
| [[File:Gsl.png|thumb|center|350 px|Figure 2: Generalized Snell's law. (Source: Yu 2011)]]
| |
|
| |
| === Hyperbolic Phase Profile ===
| |
|
| |
| Hyperbolic phase profile is commonly used for metalens, given by the equation below, where <math> \rho = \sqrt{x^2 + y^2} </math> and <math> (x,y) </math> represents the coordinate at the aperture plane; <math> f </math> is the focal length of the metalens; and <math> \lambda </math> is the wavelength of light that the metalens designed for.
| |
|
| |
| <math display="block"> \Phi (\rho) = - \frac{2\pi}{\lambda} (\sqrt{\rho ^ 2 + f ^ 2} - f) </math>
| |
|
| |
| In Fig. 3, the amplitude, which is fixed to 1, and the hyperbolic phase is shown. This metalens has the following specs:
| |
| * Radius <math> R = 100 \text{ } \mu \text{m} </math>,
| |
| * Focal length <math> f = 200 \text{ } \mu \text{m} </math>,
| |
| * Numerical aperture <math> NA = 0.89 </math>,
| |
| * Pitch size <math> 330 \text{ nm} </math>.
| |
|
| |
| [[File:Hyperbolic_phase.png|thumb|center|500 px|Figure 3: The constant amplitude profile and the hyperbolic phase profile of metalens.]]
| |
|
| |
| The point spread function (PSF) at the focal plane of the metalens is shown in Fig. 4. Strel ratio (SR) is a quantitative metric that describes the quality of the PSF. After normalizing the PSF and a perfect diffration-limited airy disk by the total energy, SR is the ratio of the peak intensity of the PSF to the peak intensity of the airy disk. At normal incident angle, the PSF is spherical-aberration-free and outperform the conventional lenses. As the incident angle increasing, the metalens shows an severe coma aberration, which limits the field of view (FOV) of metalens to great extent.
| |
|
| |
| [[File:Psf_before.png|thumb|center|500 px|Figure 4: PSFs of the hyperbolic phase metalens at different incident angles and their Strel ratio SR.]]
| |
|
| |
| == Methods ==
| |
|
| |
| === Angular Spectrum Method ===
| |
|
| |
| The Angular Spectrum Method (ASM) is a computational technique used in wave propagation analysis, particularly in optics. It models the propagation of waves by decomposing a complex wavefront into a spectrum of plane waves with spatial frequencies and angles. Each plane wave propagates independently according to its wave vector (direction of propagation), allowing the reconstruction of the wavefront at any distance from the source. ASM is especially useful for simulating wave behavior in free space or through homogeneous media.
| |
|
| |
| In the project, the optical system is straight forward. There is only one metalens with given amplitude and phase profile, and the sensor is located at the focal plane of the metalens. We assume that the object is in the far field, so its spatial structure can be represented as a distribution of spatial frequencies in the optical wave's angular spectrum. The metalens, with its specific amplitude and phase profile, acts as a spatial filter, modifying the incoming wavefront accordingly.
| |
|
| |
| By applying ASM, we can compute how the wavefront evolves from the metalens to the sensor located at its focal plane. This involves transforming the wavefront into its angular spectrum using a Fourier transform, applying the appropriate transfer function that accounts for metalens, and then inverse Fourier transforming the result to reconstruct the wavefront at the sensor plane. This approach captures diffraction and interference effects, allowing us to analyze your system's optical performance.
| |
|
| |
| === Parametrization of the Phase ===
| |
|
| |
| To model the phase profile of the metalens, we use its radial symmetry by parametrizing the phase distribution using polynomials in the radial coordinate. This approach simplifies the representation while preserving the essential features of the lens’s optical behavior. Specifically, the phase profile is expressed as a polynomial expansion in terms of the radial distance from the lens center, with coefficients determining the lens's focusing characteristics.
| |
|
| |
| <math display="block"> \Phi(\rho) = \sum_n cn(\frac{\rho}{R})^n </math>
| |
|
| |
|
| |
|
| |
| Commonly used polynomials include Zernike polynomials or simple power series, which efficiently capture aberrations and design features. This parametrization not only reduces the computational complexity of simulating the metalens but also provides flexibility in optimizing its optical performance for specific applications.
| |
|
| |
| === Optimization Algorithm ===
| |
|
| |
| == Results ==
| |
|
| |
| === Optimized Phase Profile ===
| |
|
| |
| The phase profile after optimization is shown in Figure 5. On the right, a slice of the phase at <math> y = 0 \text{ } \mu \text{m}</math> is compared to the same slice of the hyperbolic phase.
| |
|
| |
| [[File:Optimized_phase.png|thumb|center|550 px|Figure 5: Optimized phase profile and the slice of optimized phase profile.]]
| |
|
| |
| === PSF Comparison ===
| |
|
| |
| The PSFs at the focal plane of the metalens with optmized phase is shown in Fig. 6. We observed that the coma aberration is corrected to a great extent. With <math> SR > 0.1 </math> for incident angles at <math> 10^{\circ} </math>, <math> 20^{\circ} </math>, <math> 25^{\circ} </math>. On the other hand, The <math> SR </math> at the normal incident angle decreases, and we observe scattering rings which indicates spherical aberration. This implies that the optimzation algorithm is looking for a balance of the tradeoff between coma aberration and spherical aberration.
| |
|
| |
| [[File:Psf_after.png|thumb|center|500 px|Figure 6: PSFs at different incident angles of the metalens with optimized phase.]]
| |
|
| |
| To better compare the performance of the metalens, Fig. 7 shows the modulation transfer function (MTF) of the metalens imaging system. For the metalens with the hyperbolic phase, we observe that the contrast drops steeply for incident angles at <math> 10^{\circ} </math>, <math> 20^{\circ} </math>, <math> 25^{\circ} </math>. While after optimization, the MTFs across all incident angles are balanced with MTF at the normal incident angle being compromised. This is in accord with the observation from the PSF.
| |
|
| |
| [[File:MTF.png|thumb|center|500 px|Figure 7: MTFs along x-axis and y-axis of the metalens of hyperbolic phase and optimized phase.]]
| |
|
| |
| === Image Simulations ===
| |
|
| |
|
| |
|
| |
| [[File:S_sim.png|thumb|center|500 px|Figure 8: ]]
| |
|
| |
| The following images are simulated using space-varying 2D convolution algorithm. In space-varying convolution, we divide the whole FOV into <math> n </math>-by-<math> n </math>
| |
| isoplanatic patches, where we assume PSF is invariant. We implement space-invariant 2D convolution on isoplanatic patches and then merge the matches together by apply a Gaussian window. The periphery of the simulated images with hyperbolic phase metalens are a lot more blurry compared to the ones with optimized phase, which agrees with the expectation. In this simulation, we assume that the PSF across the RGB channel are the same, so color in this case is only for visual purpose.
| |
|
| |
| [[File:Simulation_img.png|thumb|center|400 px|Figure 9: Simulated images from space-varying convolution of the target image with PSFs of two metalenses.]]
| |
|
| |
| == Conclusions ==
| |
|
| |
| == Appendix ==
| |
|
| |
| You can write math equations as follows:
| |
| <math>y = x + 5 </math>
| |
|
| |
| You can include images as follows (you will need to upload the image first using the toolbox on the left bar, using the "Upload file" link).
| |
|
| |
| [[File:Snip 20210106183207.png|200px]]
| |