Experimental and numerical polarization analysis of the 3D transfer behavior in microsphere-assisted interferometry for 1D phase gratings

Lucie Hüser; Tobias Pahl; Peter Lehmann

doi:10.1051/jeos/2023029

Abstract

Enhancing the lateral resolution in optical microscopy and interferometry is of great interest in recent research. In order to laterally resolve structures including feature dimensions below the Abbe resolution limit, microspheres in the optical near-field of the specimen are shown to locally improve the resolution of the imaging system. Experimental and simulated results following this approach are obtained by a high NA Linnik interferometer and analyzed in this contribution. They show the reconstructed surface of a 1D phase grating below the resolution limit. For further understanding of the transfer characteristics, measured interference data are compared with FEM (finite element method) based simulations with respect to the polarization dependency of the relevant image information for 1D phase gratings. Therefore, the implemented Koehler illumination as well as the experimental setup utilize polarized light.

Keywords

CSI Interferometry Microsphere Polarization Resolution

1　Introduction

With the ongoing trend towards miniaturization overcoming Abbe’s diffraction limit in optical metrology is of particular interest regarding research and industrial applications. Microsphere assistance is known as a technique to achieve topographical interferometric [1–3] as well as microscopic [4, 5] measurements of structures below the resolution limit through near-field support. Photonic nanojets are frequently mentioned as the most decisive mechanism [6–8] in this context. Also the role of evanescent waves and whispering gallery modes [9, 10] is discussed. It is also found that the polarization of the illuminating light has an influence on the resolution capabilities when using microspheres [5]. More recently, a combination of the small field of view with the effectively enhanced numerical aperture of the system is found to be the most likely underlying physical mechanism [11, 12].

In optical imaging microscopy the fundamental resolution limit for periodic structures is according to Ernst Abbe given by,(1) $Λ_{\min} = \frac{λ}{2 NA},$ where Λ_min is the minimal period length of a structure that can be resolved, NA is the numerical aperture of the system and λ is the central illumination wavelength. Placing microspheres directly on the measurement object enables to reconstruct surfaces with periodic features below this limit. Several publications show also attempts to include microspheres in more practicable positioning systems, e.g. the attachment of microspheres on cantilevers or trapping them with optical tweezers [12].

In this contribution the effects of illumination polarization are examined with respect to interference microscopy. Thus, measured as well as rigorously simulated data sets of a grating with periodicity close to the resolution limit are compared with respect to polarization. Throughout this paper for TM polarization the electric field component is defined in the same plane as the 1D phase grating (x, z) similar to [13]. The results are analyzed in the 3D spatial frequency domain in order to find out the transfer characteristics and compare them for simulation and measurement results. Furthermore, a grating with period below the resolution limit was measured and the result of the interference phase analysis is shown.

2　Experimental setup and measurement results

Experimental results are obtained using a high NA Linnik interferometer (100×, NA = 0.9). The measurement process follows the principle of coherence scanning interferometry (CSI) and is described in further detail in [1]. For image acquisition a depth scan with a step size of 20 nm is performed by the piezo scanner. The resulting image stack is analyzed with respect phase and envelope of the interference signals resulting in a reconstructed surface topography of the measured specimen. Koehler illumination with an LED (royalblue, central wavelength: approximately λ = 440 nm) is used in order to provide spatially incoherent illumination. This leads to a resolution limit of Λ_min = 244 nm. A schematic representation of the experimental setup is shown in Figure 1.

Figure 1.Schematic representation of the microsphere assisted Linnik interferometer comprising two high NA microscope objective lenses (MO) and a beam splitter (BS). The illumination (blue, representing Koehler illumination) and imaging (red) beam paths are shown. The illumination beam path appearing as a plane wave below the objective lens is omitted for better visibility. The figure is not drawn to scale to enhance the comprehensibility.

When using microsphere enhancement, grating periods of 230 nm length were successfully measured. A result of a reconstructed surface topography is shown in Figure 2. In this result the magnification factor of 1.4 induced by the microsphere can be observed [11]. A phase evaluation algorithm based on a lock-in algorithm further elaborated in [14] is used to analyze the acquired interference signals. The measured specimen was a linewidth/pitch standard fabricated by supracon.

Figure 2.Measurement result of a 230 nm grating (linewidth/pitch standard). For phase analysis of the interference signals an evaluation wavelength of 650 nm was used. The dependency of the evaluation wavelength on the phase evaluation algorithm used for the surface reconstruction is elaborated in [1].

3　Rigorous simulation of the imaging process

Several attempts have been discussed in recent publications on theoretical and simulation studies of microsphere-assisted interferometry and microscopy. One of the most common approaches, which is also directly related to the analysis of photonic nanojets, is a rigorous simulation of a microsphere in free space illuminated by a plane wave source [7]. While this illustrates in detail the formation of nanojets, several aspects relevant to the overall imaging process in a microscope-based measurement setup are neglected [15].

Pahl et al. [13] developed a rigorous CSI simulation model based on FEM calculation of the electric field distribution near the sample and Fourier optics modeling of the imaging process. This model accounts for 3D conical illumination in a microscope and has recently been extended to include a microcylinder to approximate microsphere-assisted interferometry [16]. It represents a rigorous approach to the theoretical representation of interferometric imaging processes. Conical illumination is implemented with a discretization of angles of incidence homogeneously filling the numerical aperture of the system. This approach enables to consider Koehler illumination. Using the calculated near-field data sets, a far-field expansion is performed for each angle of incidence separately, leading to a Fourier optical representation of the imaging system considering the reference plane. The resulting intensity distributions are incoherently superimposed, resulting in [13]:(2) $\begin{matrix} I_{k} (x, z) = \int_{0}^{2 π} d φ_{in} \int_{0}^{θ_{in, \max}} d θ_{in} \\ P^{2} (θ_{in}) I_{k, θ_{in}, φ_{in}} (x, z) k^{2} \sin (θ_{in}) \cos (θ_{in}), \end{matrix}$ with φ_in representing the conical incident angle and θ_in,max = arcsin(NA) the maximum opening angle of the aperture. With P(θ_in) the pupil function for both, the reference and the object arm is considered. The intensity $I_{k, θ_{in}, φ_{in}}$ is the intensity distribution for each angle φ_in and θ_in,max calculated from the simulated electric field distributions. Taking the spectral distribution S(k) of the light source with the wave number k = 2π/λ into account the total intensity $I (x, z) = \int d k S (k) I_{k} (x, z)$ is calculated.

Finally, signal processing algorithms (e.g. phase and envelope evaluation algorithms) are applied to the a stack of simulated interference images. This corresponds to the signal processing performed on experimentally obtained image stacks, thus providing direct comparability.

4　Analysis of the results in the 3D spatial frequency domain

The transfer characteristics of an interferometer can be obtained from the 3D spatial frequency domain representation of a measured data set [17], which enables a profound understanding of the 3D measurement process. This method is applied to microsphere-assisted interferometry to gain further insight into the mechanisms underlying the lateral resolution enhancement.

For further analysis the data is transferred to the 3D spatial frequency domain also called q–space, which is based on the incident wave vectors k_in and the scattered wave vectors k_s defined as, $k_{in} = k (\begin{array}{l} \sin θ_{in} \cos φ_{in} \\ \sin θ_{in} \sin φ_{in} \\ - \cos θ_{in} \end{array}), k_{s} = k (\begin{array}{l} \sin θ_{s} \cos φ_{s} \\ \sin θ_{s} \sin φ_{s} \\ \cos θ_{s} \end{array}),$ and therefore q = k_s − k_in holds.

The resulting 3D spatial frequency distribution is limited by the Ewald limiting sphere, which is shown in Figure 3 and explained in further detail in [17]. Thus, the diffraction orders q_x,n = 2πn/Λ with n ∈ {…, −1, 0, 1,…} for a grating of period length Λ in x-direction appear as distinct vertical lines within the Ewald limiting sphere, as exemplarily shown in Figure 3.

$Representation of the Ewald limiting sphere for monochromatic light of wave number k0 in q – space. The q – values belonging to specular reflection (qx = 0) and backscattering are marked. The most outlying edges at qx = ±2k0NA correspond to the highest resolvable spatial frequencies for the maximum incident angle θin,max. kr is representing the reflected wave vector. As an example, the first diffraction orders for a periodic structure with a period length Λ were inserted.$

Figure 3.Representation of the Ewald limiting sphere for monochromatic light of wave number k₀ in q – space. The q – values belonging to specular reflection (q_x = 0) and backscattering are marked. The most outlying edges at q_x = ±2k₀NA correspond to the highest resolvable spatial frequencies for the maximum incident angle θ_in,max. k_r is representing the reflected wave vector. As an example, the first diffraction orders for a periodic structure with a period length Λ were inserted.

To illustrate the connection between the q–space representation of the interferometric measurement data stack ${\tilde{I}}_{IF} (q)$ and the 3D transfer function H_IF(q) the following equation can be used for the interference part of the signal [17]:(3) ${\tilde{I}}_{IF} (q) = [{\tilde{U}}_{0} (q) \otimes {\tilde{U}}_{ref} (q)]_{q_{x}, q_{y}} H (q),$ (4) $= {\tilde{U}}_{0} (q) H_{IF} (q) = O (q) H_{IF} (q),$ where ${\tilde{U}}_{0} (q)$ denotes the scattered light field of the phase object in q–space. Since the reference plane is a mirror, all contributions to the signal formation lie on the q_z-axis, and the contribution of ${\tilde{U}}_{ref} (q)$ to the correlation (denoted with ⊗) can be considered by H_IF(q). O(q) represents the 3D object spectrum if the surface is represented by a thin foil [18]. Using the inverse Fourier transform, the interference measurement data stack in the spatial domain is given by the convolution of the surface foil o(x, y, z) and the real part h_IF(x, y, z) of the inverse Fourier transformed function H_IF(q), therefore (with * denoting the convolution),(5) $I (x, y, z) = 2 o (x, y, z) * Re {h_{IF} (x, y, z)} .$

This shows that by analyzing the transfer behavior of the system including the microspheres, which strongly affect the imaging capabilities, important insight into relevant mechanisms can be gained.

Figures 4a–4d show measured and simulated interferometric data sets for a rectangular silicon grating (RS-N, SiMETRICS) with Λ = 300 nm and a nominal height of 140 nm. The phase grating leads to a phase modulation of the interference fringes along the x-axis, which is centered around the virtual image plane of the microsphere [1].

Figure 4.Interferometric measurement data obtained from a rectangular grating (SiMETRICS RS-N, Λ = 300 nm) using royalblue (a) TM and (b) TE polarized light through a microsphere (SiO₂, 5–9 μm diam.), simulated data sets of a similar grating with (c) TM and (d) TE polarized light assuming a microcylinder of 5 μm diameter. The corresponding 3D spatial frequency representations are depicted in (e) – (h). The results shown in (a) and (e) are also analyzed in [11].

Comparing the different polarization configurations it can be observed that for TM polarized light the phase modulation induced by the grating is much stronger in both, measurement and simulation results. Besides, a qualitative correspondence of measured and simulated results can be seen. Therefore, the simulation approach in principle reproduces the transfer characteristics of microsphere-assisted interferometry.

In Figures 4e–4h the 3D spatial frequency domain representations in the q_xq_z-plane of the different data sets are depicted. The first order diffracted light components appear at q_x ≈ ±15 μm⁻¹ (considering the additional magnification through the microsphere). Instead of distinct vertical lines the diffraction maxima are blurred. This effect can be explained by the limited field of view while imaging through a microsphere. In Fourier optics the field of view corresponds to the area of integration for the Fourier integral, which can be considered by multiplication with an appropriate function describing the limited extension of the illuminated area in the object space under the microsphere. This leads to a convolution with the Fourier transformed illumination function in the Fourier domain and thus explains the broadening of the diffraction orders.

Comparing the frequency representations for the TM and TE case, the intensity maxima related to first order diffraction are lower for the TE case, which is in agreement with the weaker phase modulation visible in Figure 4b. In comparison with the simulated results in the 3D spatial frequency domain, some effects occur differently. The magnification introduced by the microcylinder is higher (compare the extension in x-direction of the phase modulation in Figs. 4a and 4c), leading to slightly different positions of the diffraction orders. Furthermore, due to the use of a 2D microcylinder instead of a sphere in the simulation further effects appear [11]. Also slight differences in the diameter of the utilized microspheres compared to the microcylinder affect the results.

For a better differentiation of the results shown in Figure 4 of simulated interferometric data sets without microspheres are shown in Figure 5. For the TM polarization in Figure 5a, a stronger modulation of the phase laterally can be seen here. This corresponds to the results recorded with microsphere assistance.

Figure 5.Simulated interferometric data set with a rectangular grating structure similar to the SiMETRICS RS-N grating as a specimen (Λ = 300 nm) for the (a) TM and (b) TE polarization case.

5　Conclusion

Considering the transfer behavior of an optical system enhanced with microspheres placed in the near-field of a specimen, gives insight into the relevant mechanisms. Here, the influence of the polarization of the illuminating light is examined. It is pointed out how the transfer behavior of the optical imaging process including the microsphere contributes to a deeper understanding of the relevant mechanisms for imaging capabilities and resolution enhancement.

Analyzing the polarization dependency, TM polarization should be preferred to obtain the relevant phase information of 2D topographies, which are shift-invariant in y-direction. This is in agreement with previous results obtained for grating structures using an interference microscope without microsphere assistance and could be confirmed by rigorous simulations based on an FEM model [13].

Further, the reconstructed topography of a grating structure with a period length below Abbe’s resolution limit is shown to demonstrate the improved resolution enabled by the microsphere introduced into the optical system.

References

[1] L. Hüser, P. Lehmann. Microsphere-assisted interferometry with high numerical apertures for 3D topography measurements. Appl. Opt., 59, 1695(2020).

[2] S. Perrin, Y.J. Donie, P. Montgomery, G. Gomard, S. Lecler. Compensated microsphere-assisted interference microscopy. Phys. Rev. Appl., 13, 1(2020).

[3] F. Wang, L. Liu, P. Yu, Z. Liu, H. Yu, Y. Wang, W.J. Li. Three-dimensional super-resolution morphology by near-field assisted white-light interferometry. Sci. Rep., 6, 24703(2016).

[4] Z. Wang, W. Guo, L. Li, B. Luk’yanchuk, A. Khan, Z. Liu, Z. Chen, M. Hong. Optical virtual imaging at 50 nm lateral resolution with a white-light nanoscope. Nat. Commun., 2, 218(2011).

[5] A. Darafsheh, G.F. Walsh, L. Dal Negro, V.N. Astratov. Optical super-resolution by high-index liquid-immersed microspheres. Appl. Phys. Lett., 101, 141128(2012).

[6] H. Yang, R. Trouillon, G. Huszka, M.A. Gijs. Super-resolution imaging of a dielectric microsphere is governed by the waist of its photonic nanojet. Nano Lett., 16, 4862(2016).

[7] A. Darafsheh. Photonic nanojets and their applications. J. Phys.: Photon., 3, 022001(2021).

[8] S. Lecler, S. Perrin, A. Leong-Hoi, P. Montgomery. Photonic jet lens. Sci. Rep., 9, 1(2019).

[9] Y. Ben-Aryeh. Increase of resolution by use of microspheres related to complex Snell’s law. J. Opt. Soc. Am. A, 33, 2284(2016).

[10] S. Zhou, Y. Deng, W. Zhou, M. Yu, H.P. Urbach, Y. Wu. Effects of whispering gallery mode in microsphere super-resolution imaging. Appl. Phys. B: Lasers Opt., 123, 1(2017).

[11] L. Hüser, T. Pahl, M. Künne, P. Lehmann. Microsphere assistance in interference microscopy with high numerical aperture objective lenses. J. Opt. Microsyst., 2, 044501(2022).

[12] A. Darafsheh. Microsphere-assisted microscopy. J. Appl. Phys., 131, 031102(2022).

[13] T. Pahl, S. Hagemeier, M. Künne, D. Yang, P. Lehmann. 3D modeling of coherence scanning interferometry on 2D surfaces using FEM. Opt. Express, 28, 39807(2020).

[14] P. Lehmann, S. Tereschenko, W. Xie. Fundamental aspects of resolution and precision in vertical scanning white-light interferometry. Surf. Topogr.: Metrol. Prop., 4, 024004(2016).

[15] Z. Wang, B. Luk’yanchuk. Super-resolution imaging and microscopy by dielectric particle-lenses, 371-406(2019).

[16] T. Pahl, L. Hüser, S. Hagemeier, P. Lehmann. FEM-based modeling of microsphere-enhanced interferometry. Light: Adv. Manuf., 3, 1(2022).

[17] P. Lehmann, M. Künne, T. Pahl. Analysis of interference microscopy in the spatial frequency domain. J. Phys.: Photon., 3, 1(2021).

[18] J. Coupland, R. Mandal, K. Palodhi, R. Leach. Coherence scanning interferometry: Linear theory of surface measurement. Appl. Opt., 52, 3662(2013).