Abstract
Closed formulas are presented for the electromagnetic field of given power in the lens pupil, which maximizes the longitudinal electric field when focusing through an interface at arbitrary depth along the optical axis. The optimum pupil field is found to be a continuous, monotonously increasing function of the radial pupil coordinate, which differs considerably from the commonly used annular illumination. Several cases of pupil fields and focused fields are shown for different materials, NA, and focusing depth. Also, the effect of absorbing media is considered.
© 2016 Optical Society of America
1. INTRODUCTION
When focusing a polarized beam with a diffraction-limited lens of high numerical aperture, refraction of the vector electromagnetic field must be considered, using vectorial diffraction theory proposed first by Ignatovky and later by Richards and Wolf to obtain the vectorial Airy pattern. With a radially polarized incident beam, the strong bending of the light rays causes the longitudinal electric field component in the focal point to become dominant. When focusing in air or in an immersion system, the maximum of the longitudinal field component is attained precisely in the focal point. For the same numerical aperture, the symmetrical squared amplitude distribution of the longitudinal field component is more confined than the elliptically shaped squared modulus of the total electric field of the classical Airy pattern obtained by focusing a linear polarized plane wave.
A small and symmetric focused spot has applications in optical lithography [1], microscopy [2–4], laser writing [5], and optical trapping [6]. In many of these applications, the light is focused through an interface into a medium of a different refractive index, which can have considerable consequences for the distribution of the different components of the focused vectorial field. Annular illumination [7] and shaded-ring filters [3] are commonly used in order to increase the relative strength of the longitudinal field, although this can lead to undesired strong sidelobes and loss of power. A new optimization of the longitudinal field in focus is required to take into account the influence of the susbtrate while preserving the power transmission. This can be considered an extension of the work of Urbach in [8] for the case of a layered interface.
Here, we present a closed form expression for the electric field in the lens pupil, which maximizes the longitudinal electric field component in the focal point inside a single-layer medium for a fixed power. The optimum incident field in the pupil is still radially polarized, but it depends on the position of the distance of the focal point from the interface. The optimum amplitude distribution is a continuous, monotonically increasing function of the radial coordinate in the pupil of which the precise shape depends on the numerical aperture. Furthermore, for focusing behind an interface, there is an additional phase distribution of the pupil field, which depends on the depth of the focal point inside the material.
2. OPTIMIZATION PROBLEM
We proceed by maximizing the longitudinal field component on the optical axis at an arbitrary depth behind the interface and define the pupil field, which attains this goal as the optimum field. We will first derive the angular spectrum of the optimum electric field in the focal region as a generally valid result that does not depend on the diffraction model used to simulate the lens. Then, we determine the optimum pupil field using vectorial diffraction theory proposed first by Ignatovsky and later by Richards and Wolf [9,10].
With respect to Cartesian coordinates , we consider the complex amplitudes of a time–harmonic electromagnetic field and . The incident field is propagating in the positive direction in the half-space with real permittivity . It is incident on the substrate, which fills the half-space and has a permittivity with non-negative imaginary part. If the imaginary part is positive, the material the substrate consists of is absorbing. The time dependence of the fields is given by the factor , with , which is omitted throughout the rest of the paper.
The incident field in can be expanded in plane waves. It is assumed that the sources that generate the field are at a distance that is far enough so that the waves that are evanescent in the positive direction can be neglected and, hence, only plane waves that propagate in the positive direction are included. We furthermore assume that only plane waves of which the wave vectors are inside the cone are contributing to the field, where NA is the numerical aperture of the lens in vacuum. Hence, the plane wave expansion of the incident field in is
with being the wave vector with because the waves are propagating in the positive direction, and the wavenumber in vacuum, , where is the maximum angle between the positive direction and the wave vectors in the first medium. We will define the vectors of a basis in space for medium 1 and 2 labeled by the superscript j, respectively: where is defined as the transverse wavenumber. Every pair with and corresponds to a wave vector. Because the electric field amplitudes of the plane waves satisfy , the vector has only a and component:Note that when the substrate is absorbing, is complex and, hence, a -polarized wave is elliptically polarized. An -polarized wave, however, is always linearly polarized.
From Eq. (3), it becomes clear that the components of the plane waves do not contribute to the longitudinal field in the focal spot; hence, the optimum field consists only of the component.
The time-averaged incident power flow in the positive direction through any plane between the lens and the interface can be computed with the Poynting vector and Plancherel’s theorem:
The presence of the interface affects the power flow due to reflection in the negative direction, as well as possible absorption in the second medium.
For the field transmitted in , a similar plane wave expansion holds as Eqs. (1) and (2), except that the permittivity must be replaced by . We write instead of in :
We have
where is the Fresnel coefficient for -polarization, given following the convention of Born and Wolf as [11] where the subscripts 1 and 2 denote the medium in which the quantity is measured, and the dependence of on and is not explicitly written. Therefore, for , and, in particular, for ,In this paper, we always place the geometrical focus of the lens at the origin , i.e., in the interface. At a fixed depth , we will consider as a functional of . One may assume that is positive real because this assumption amounts to stating that the phase of the component of the electric field is zero at time , and this can always be achieved by a shift in time. Then, we will find the plane wave amplitudes , which maximize the axial longitudinal field amplitude in the point under the constraint that the mean power flow through the exit pupil be held constant, i.e.,
Solving Eq. (11) can be done with the Lagrange multiplier rule,
where denotes the Lagrange multiplier. From Eq. (12) follows a system of linear equations for the real and imaginary parts of the plane wave amplitudes through which they can be expressed in the Lagrange multiplier:The Lagrange multiplier is obtained by substituting Eq. (13) into the power constraint:
The magnitude of the plane wave amplitudes is independent on the location when medium 2 is nonabsorbing. There is, however, an additional phase modulation that is nonzero for all . The integral in Eq. (14) cannot be computed analytically due to the Fresnel coefficients and must be evaluated numerically.
After substituting Eq. (13) into Eq. (9), the incident electromagnetic field in a point with cylindrical coordinates and can be expressed as
and where , , and are unit vectors of the cylindrical coordinate system and where the function is defined byIn Eq. (17), the special functions represent Bessel functions of the first kind and order .
3. OPTIMUM PUPIL FIELD DISTRIBUTION
Now we will address the problem of how to realize the optimum field as computed in Section 2 by focusing with a high numerical aperture lens. The optimum pupil field is based on a model for high NA focusing derived by Ignatovski and Richards and Wolf [9,10]. Spatial light modulators can then be used to shape the amplitude, phase, and polarization in the pupil. In the model of Ignatovski and Richards and Wolf, the electric field at a pupil point with polar coordinates is related to the field of the plane wave with transverse wave vector such that
where denotes the focal distance. The complex pupil field is given by [12]The optimum pupil field is radially polarized and has rotational symmetry in amplitude and phase. Its amplitude is a monotonously increasing function of the radial coordinate . This amplitude distribution depends on the dielectric permittivities on both sides of the interface and on NA but not on the focal plane location if and only if the media are nonabsorbing (). The phase does, however, depend on the assumed position of the focal point.
When focusing behind the interface in the second medium , a -dependent phase modulation is introduced. In absorbing media, the exponential factor also contains additional amplitude modulation through the imaginary part of , which will make the optimum pupil amplitude become dependent. Figure 1 shows the pupil field amplitudes for several NA and the same total power with , and , and . As a function of radius, these amplitude distributions appear to be similar, although the angular distribution of these incident fields is significantly different for different values of NA. Nevertheless, the dependency of the optimum pupil field on NA is much weaker than when focusing in air [8]. Particularly for far off-axis rays (close to ), the Fresnel transmission coefficients in Eq. (19) have low values, thus slightly flattening the field amplitude in this regime. The choice of was made throughout this paper in order to demonstrate the largest possible effect of the interface. When the difference in permittivity before and after the interface becomes smaller, the influence of the interface is reduced and the system tends toward an immersion system. Other parameters are chosen in a way that they are representative for the applications mentioned in the introduction, but, since we derived closed expressions, they can easily be changed using Eqs. (17) and (19).
Figure 2 shows the pupil field amplitudes for several and the same total power with , and , with and . The total incident power is held constant. The pupil amplitudes in Fig. 1 and 2 are independent of because is real.
When absorption is included, it simultaneously affects the amplitude and phase of the pupil field. To demonstrate the effect on amplitude, Fig. 3 displays the optimum pupil amplitude for a medium with (, ) and without absorption (). The field is maximized at a depth of . Furthermore, , . Introducing absorption in medium 2 further reduces the optimum pupil field amplitude at large pupil radii. This effect can be understood, as light rays passing through the edge of the lens travel a greater distance to reach the focal spot, as compared with rays closer to the optical axis, and would be more strongly attenuated along the way.
Finally, the phase front of the optimum pupil field is generally nonuniform unless . Figure 4 shows some phase distributions of the optimum pupil fields, which maximize the longitudinal field in focus at different depths . Other conditions chosen for this example are , , and . Note that is measured in wavelengths units in medium 2.
To account for values of , a defocus term arises in the optimum pupil phase whereas absorption leads to a nonspherical wavefront.
4. VISUALIZATION OF FOCUSED FIELDS
The optimum electromagnetic field in focus is rotationally symmetric around the axis; hence, both and are independent of . For a given optical system with and an air/medium interface with and (note that common organic photoresists {Shipley-series [13], MAN-series [14]} have similar dielectric constants in the near-UV range), the focused field amplitudes resulting, which follow from Eq. (15), are plotted for the optimized as well as a uniform-amplitude (top hat) pupil field in Fig. 5 (note that the uniform amplitude field is zero in the center to avoid a polarization singularity). The interface is located such that the longitudinal fields peak at a depth . The longitudinal, transversal, and total focused field amplitudes are denoted by , , and , respectively, and the total field peaks at unity.
For the case of an incident field with a top hat amplitude profile, the transverse electric field is the dominant field component. After optimization, the longitudinal field dominates over the transverse field, leading to a narrower total field distribution. This effect becomes stronger when considering the field intensities (squared magnitudes). Interestingly, the transverse field component of the optimized pupil field also becomes narrower, thus further reducing the width of the total field distribution.
Figure 1 seems to suggest that the optimum pupil field amplitude as a function of pupil radius depends only weakly on the numerical aperture. To demonstrate the dependency of the focused fields on NA, Fig. 6 shows the maximum longitudinal field , which is positive real by definition, for a focused field maximized on the optical axis at . Taking different values of demonstrates that, as expected, the maximum longitudinal field decreases when either the real or imaginary part of the second dielectric permittivity increases.
5. CONCLUSION
To conclude, we present a closed-form expression for the pupil field, which, when focused through an interface, maximizes the longitudinal electric field component at a selected depth in the substrate for a fixed total power. The optimum pupil field is radially polarized, has rotational symmetry in amplitude and phase, and the field amplitude increases continuously with distance from the optical axis. When focusing into another medium, the phase front shows an extra depth-dependent curvature. In media with absorption, both the phase and amplitude of the optimum pupil field change due to the imaginary part of the respective dielectric constant.
Because the longitudinal field has spatial dimensions below the diffraction limit, this framework can be used to obtain a small and symmetric focused spot inside a possibly absorbing medium at arbitrary depth, which has been implemented in the field of optical lithography [15], and has further applications in microscopy and laser machining. It is also useful for applications that make use of a strong longitudinal field, such as particle trapping [16].
REFERENCES
1. B. Jia, H. Kang, J. Li, and M. Gu, “Use of radially polarized beams in three-dimensional photonic crystal fabrication with the two-photon polymerization method,” Opt. Lett. 34, 1918–1920 (2009). [CrossRef]
2. D. Biss and T. Brown, “Cylindrical vector beam focusing through a dielectric interface,” Opt. Express 9, 490–497 (2001). [CrossRef]
3. H. Lin, B. Jia, and M. Gu, “Generation of an axially super-resolved quasi-spherical focal spot using an amplitude-modulated radially polarized beam,” Opt. Lett. 36, 2471–2473 (2011). [CrossRef]
4. L. E. Helseth, “Roles of polarization, phase and amplitude in solid immersion lens systems,” Opt. Commun. 191, 161–172 (2001). [CrossRef]
5. A. Kumar and M. C. Gupta, “Laser machining of micro-notches for fatigue life,” Opt. Lasers Eng. 48, 690–697 (2010). [CrossRef]
6. O. M. Marago, P. H. Jones, P. G. Gucciardi, G. Volpe, and A. C. Ferrari, “Optical trapping and manipulation of nanostructures,” Nat. Nanotechnol. 8, 807–819 (2013). [CrossRef]
7. R. Dorn, S. Quabis, and G. Leuchs, “Sharper focus for a radially polarized light beam,” Phys. Rev. Lett. 91, 233901 (2003). [CrossRef]
8. H. P. Urbach and S. F. Pereira, “Field in focus with a maximum longitudinal electric component,” Phys. Rev. Lett. 100, 123904 (2008). [CrossRef]
9. V. S. Ignatovski, “Diffraction by a lens having arbitrary opening,” in Trans. Opt. Inst. Petrograd I, paper IV (1919, in Russian).
10. B. Richards and E. Wolf, “Electromagnetic diffraction in optical systems,” Proc. R. Soc. London 253, 358–379 (1959). [CrossRef]
11. M. Born and E. Wolf, Principles of Optics, 7th ed. (Cambridge University, 1999).
12. H. P. Urbach and S. F. Pereira, “Focused fields of given power with maximum electric field components,” Phys. Rev. A 79, 013825 (2009). [CrossRef]
13. “Shipley co. S1800 series data sheet,” cmnst.ncku.edu.tw/ezfiles/23/1023/img/127/s1800seriesDataSheet.pdf, accessed 6 October 2015.
14. “ma-n co. 2400 series data sheet,” www.microchem.com/PDFs_MRT/ma-N%202400%20overview.pdf, accessed 8 October 2015.
15. K. Ushakova, Q. Y. van den Berg, S. F. Pereira, and H. P. Urbach, “Demonstration of spot size reduction by focussing amplitude modulated radially polarized light on a photoresist,” J. Opt. 17, 125615 (2015). [CrossRef]
16. Q. Zhan, “Trapping metallic Rayleigh particles with radial polarization,” Opt. Express 12, 3377–3382 (2004). [CrossRef]