Active Focusing of Light in Plasmonic Lens via Kerr Effect
 Author: Nasari Hadiseh, Abrishamian Mohammad Sadegh
 Organization: Nasari Hadiseh; Abrishamian Mohammad Sadegh
 Publish: Journal of the Optical Society of Korea Volume 16, Issue3, p305~312, 25 Sep 2012

ABSTRACT
We numerically demonstrate the performance of a plasmonic lens composed of an array of nanoslits perforated on thin metallic film with slanted cuts on the output surface. Embedding Kerr nonlinear material in nanoslits is employed to modulate the output beam. A two dimensional nonlineardispersive finitedifference timedomain (2D NDFDTD) method is utilized. The performance parameters of the proposed lens such as focal length, fullwidth halfmaximum, depth of focus and the efficiency of focusing are investigated. The structure is illuminated by a TMpolarized plane wave and a Gaussian beam. The effect of the beam waist of the Gaussian beam and the incident light intensity on the focusing effect is explored. An exact formula is proposed to derive electric field
from electric flux densityE in a KerrDispersive medium. Surface plasmon (SPs) modes and FabryPerot (FP) resonances are used to explain the physical origin of the light focusing phenomenon. Focused ion beam milling can be implemented to fabricate the proposed lens. It can find valuable potential applications in integrated optics and for tuning purposes.D

KEYWORD
Plasmonics , MIM waveguide , NDFDTD

I. INTRODUCTION
While onchip information processing can efficiently speed up by employing photonic integrated circuits (PICs), there was a significant challenge to realize them by a density comparable to the electronic counterparts due to the presence of the optical diffraction limit. Surface plasmonbased photonics or plasmonics, a recently emerged device technology, can propose a solution to this dilemma, because plasmonics has both the large data carrying capacity of photonics and the miniaturization of electronics [13]. Although surface plasmon polaritons (SPPs), surface electromagnetic waves coupled to collective oscillations of free electrons in a metal [4], were first theoretically proposed in 1957 [5], the significant growth and interest in plasmonics is formed since Ebbesen first reported the extraordinary optical transmission through a 2D metallic hole array in 1998 [6]. Nano fabrication techniques, such as electronbeam lithography and focused ion beam milling can be implemented to pattern metallic structures at the nanoscale and make it possible to study the interaction of these nanostructures with light [7, 8].
The strong localization of SPPs at the metaldielectric interfaces has attracted tremendous interests of researchers to SPPbased subwavelength guiding components [9] such as arrayed nanoparticles [10], Vgrooves [11], wedges [12], nanowires [13], dielectricloaded metal films [14] and plasmonic slot waveguides [15]. Among them, MetalInsulatorMetal (MIM) waveguides are promising in design of a variety of nanoscale plasmonic devices because of possessing high group velocity over a wide frequency range from DC to visible, strong confinement of light waves and also simple fabrication and integration into optical circuits [16, 17]. Although InsulatorMetalInsulator (IMI) waveguides have the advantage of less loss for the longer propagation distance, but their light confinement capability into subwavelength scales is poor [18]. Numerous functional plasmonic MIM structures like directional couplers [19], filters [20], Yshaped combiners [21], sensors [22], Ushaped waveguides [23], Bragg reflectors [24], etc, have been numerically and/or experimentally investigated.
Refractive lenses are one of the most ubiquitous optical components with applications ranging from imaging to concentrating light, but their light confinement capability deteriorates as their size approaches the wavelength of light due to the diffraction limit. Surfaceplasmon based lenses or plasmonic lenses as an alternative to the ordinary refractive lenses are capable of super focusing beyond the diffraction limit and have great applications in optical data storage, nanofabrication, single molecular biosensing, circular polarizer analyzer and etc [7, 2527]. Single subwavelength slit surrounded by surface corrugation [28] or by chirped dielectric surface gratings [29], chirped circular slits corrugated on metallic film [30], quasiperiodic array of nanoholes [31] and nanometric crossshaped aperture arrays [32] in a metal screen are some of the reported design principles to implement the focusing capability of plasmonic lenses. Nanoslits perforated on thin metallic films are recently employed to manipulate the phase front profile [3335]. The phase shift experienced by light passing through nanoslits is sensitive to its width and depth and also the refractive index of incorporated material. The required phase front profile for focusing action can be achieved by appropriate adjustment of the properties of each slit.
Active control of plasmons which is one of the greatest challenges in recent years is needed to achieve tunable plasmonic devices, besides that, more specifically; the ability to control light with light is required in alloptical signal processing in PICs and optical computing. This can be realized by employing active materials whose properties can be altered by some form of stimulation [3, 36]. Switches and modulators on the basis of the electrooptic effect [37, 38], the magnetooptic effect [39], the thermooptic effect [40], the plasma dispersion effect [41], the plasmonic excitation of quantum dots [42] and the Kerr nonlinearity [43] are investigated. In the plasmonic lenses formed by nanoslits perforated on thin films, modulating the output beam can be accomplished by incorporating nanoslits with a type of active material like Kerr nonlinear material [44] or anisotropic nematic liquid crystal [45].
In this paper we numerically explore the focusing effect of a plasmonic lens formed by nanoslits perforated on a thin metallic film with slanted cuts on the output surface and incorporated by Kerr nonlinear material. The performance parameters of the proposed lens including focal length (FL), depth of focus (DOF) defined at the half maximum intensity, fullwidth halfmaximum (FWHM) and the efficiency of focusing are investigated. The lens is illuminated by a Gaussian beam and plane wave and the possibility of tuning the focal length by the intensity of incident light and the beam waist of the Gaussian beam is investigated.
This paper is organized as follows; In Section 2 the structure of proposed lens and the analysis method are presented. Section 3 elucidates results of the detailed simulations of the lens and the conclusion will be in Section 4.
II. DEVICE STRUCTURE AND ANALYSIS METHOD
The simulated structure is displayed in Fig. 1(a). The slit widths in consequence from top to bottom are 120, 100, 80, 60, 80, 100 and 120 nm with 400 nm center to center spacing between any two adjacent slits. The underlying physics involved in design of our proposed lens can be discussed as follows. A narrow slit surrounded by metallic walls is the basic building block of the lens. If the space between every two adjacent slits is much more than the metal’s skin depth at operating wavelength, these slits can be assumed as isolated MIMs besides that the slits widths in our structure are much smaller than the operating wavelength, so just the propagation of fundamental SPP mode in MIM waveguide is considered. The complex propagation constant
β of the transverse magnetic (TM) mode in this MIM waveguide is determined by the following dispersion relation [15]:Where
and
are the propagation constants of the dielectric and metal respectively.
k _{0}=2π/λ_{0},ε_{m} ,ε_{d} andω are wave number in free space, the dielectric constants of metal cladding and dielectric core and the width of the MIM waveguide, respectively. The phase delay introduced by a slit is calculated from [34]:Where
Δ？ _{1} andΔ？ _{2} are accompanied phase changes at the entrance and exit interfaces and cancel each other if the medium on the illuminated and unilluminated sides of the lens be the same,d is the slit depth,θ originates from multiple reflections of light between the entrance and exit surfaces. But the phase delay is dominated by theRe (βd ) [34]. So if the slits be the same in width and incorporated material, the phase delay will be proportional to slit depth, on the other hand narrower slits bring more phase delays, in the result of more residence of mode in metal cladding which causes the light to travel slower [7]. In our proposed lens the slit depth and width increases and decreases respectively toward the center of the structure and hence creates the required phase front profile for light focusing which is shown in Fig. 1(b).The performance of the proposed lens has been simulated by a twodimensional nonlineardispersive finitedifference timedomain (2D NDFDTD) numerical method, assuming the slit lengths (in z direction) to be infinite. This assumption is acceptable for slit lengths larger than 15 μm [46] and by using the fact that the error in 2D simulations in comparison with 3D simulations will be negligible (Less than 10 percent). Obviously more accurate results by investigating the effect of finite slit length can be obtained by 3D NDFDTD but at a cost of large amounts of computer memory and increased computing time. The outer boundary of the computation lattice is terminated to the convolutional perfectly matched layer (CPML) to dissipate outgoing waves [47]. The secondorder Lorentz dispersion model is employed to characterize the frequencydependent relative permittivity of silver, when time dependency is taken
exp (jωt ) [48]:Where ω_{0} and Γ are resonant frequency and damping coefficient.
ε _{∞}=ε_{0} (1+χ_{∞}) determines the response of the medium for frequencies far above resonant frequency andχ _{0} is the DC response of the polarization to the electric field. Knowingε_{R} at operating frequency which can be easily obtained from hand books like [49] and setting reasonable values forε _{∞} andχ _{0} likeχ _{0}=10 andε _{∞}=2.89 ε_{0}, we can calculateω _{0} and Γ from [48]:A TMpolarized plane wave of 850 nm wavelength consisting of
E_{x} ,E_{y} , andH_{z} field components illuminates the structure. The relative permittivity of silver at this wavelength isε_{m} = 33.22 j 1.17. The grid sizes are chosen Δx = Δy = 5 nm and the time step is setachieved by Courant stability condition, where
c is the speed of light in free space.The nanoslits are filled with Kerr nonlinear material. For a Kerr nonlinear medium the dielectric constant
ε_{d} depends on the intensity of incident light by [50]:The linear dielectric constant
ε_{l} is set as 2.25. The thirdorder nonlinear susceptibility χ^{(3)} = 1.4×10^{10} esu is chosen as a typical value of nonlinear optical materials such as InGaAsP [50].To simulate the performance of proposed lens, at first step, The electric flux density
D is obtained by solution to Maxwell’s equations; ∇× = μ∂E /∂H t , and ∇× =∂H /∂D t where andE are the electric and magnetic fields, respectively and μ is the magnetic permeability, according to Yee algorithm. The electric field can be derived from following equation [47]:H Where ε_{0} is the free space permittivity,
P ^{ L} andP ^{ NL} are linear and nonlinear polarization vectors associated with the Lorentz model of silver and Kerr nonlinearity, respectively [47].Where
A = (2ω^{2} _{0}Δt ^{2})/((Γ/2)Δt + 1),B = (ΓΔt  2)/(ΓΔt + 2),C = (ε_{0}χ_{0}ω_{0} ^{2}Δt ^{2})/((Γ/2)Δt + 1) and Δt is the time step. By substituting Eq. (8) and Eq. (9) in Eq. (7), the electric field can be obtained by considering two methods, by using Newton iteration [47]:Where
m = 0 ,1 ,2 ,？ , ^{<m>} is the approximation ofE ^{n+1} at the mth iteration andE ^{<0>}=E E ^{n} . This iteration process continues until sufficient accuracy is obtained. Another way is by considering Eq. (7) as a cubic equation:Where
a =ε _{0}χ^{(3)},b =ε _{0}ε _{∞}, andthat has a solution for
^{n+1}:E This result is used in following simulations and provides more accuracy than the first method for few iterations.
III. SIMULATION RESULTS AND DISCUSSION
As mentioned previously, change in the dielectric constant of incorporated material in nanoslits can be implemented to modulate the output beam. For different values of
ε_{d} the dispersion equation is solved forβ by usingMathematica . The phase delay of transmitted light through nanoslits in the situation of illuminating the structure by plane wave is calculated according to Eq. (2) and the results are presented in Fig. 2(a). The results of change in the propagation length defined asL = (2Im (β ))^{？1} [51] by varyingε_{d} are illustrated in Fig. 2(b).The nanoslits in our structure are filled with Kerr nonlinear material. By increasing the intensity of incident light, in the situation that FabryPerot (FP) resonances
occur in the slits in the result of impedance mismatch at the entrance and exit surfaces of each slit, the dielectric constant of incorporated nonlinear material increase. By increasing
ε_{d} , the relative phase difference and hence the phase front profile changes so tuning the focal length by incident light intensity becomes possible. When the incident light intensity is 13.2 MW/cm^{2}, all of the slits are out of FP resonances and the calculated electric field intensity on the xaxis of the lens and on the focal plane which represent the FL, DOF and FWHM respectively are shown in Fig. 3(a), (b).The results depicted in Fig. 4(a) illustrate the movement of focal point toward lens by increasing the intensity of incident light, which is in the result of formation FP resonances in three central slits, while the outer slits remain out of FP resonances. Figure 4(b) represents the dependence of the dielectric constant of incorporated Kerr nonlinear material in nanoslits located at y=0 and 400 nm to the incident light intensity.
Table 1 lists the focusing properties of the lens for different values of incident light intensity. FWHM is approximately unchanged and remains equal to 590 nm for all four cases but the DOF parameter which determines the working distance for beaming application, decreases.
Certainly, rather than the absorption of a portion of light in the structure, there is reflection of light from the surface of the lens in the illuminated side, so the proposed lens cannot transmit 100% of the incident light. Also, it should be pointed out that in these kinds of lenses which are formed on the basis of nanoslits perforated on thin metallic films, because of the low decoupling rate of surface plasmon modes to propagating electromagnetic wave modes, most of the energy remains at the output surface of the lens and hence the efficiency of energy passing through the lens will be much more than the efficiency of focusing [25]. The efficiency of focusing which is obtained by dividing
the intensity of focal point to the incident intensity increases by increasing the incident light intensity; attributed to the formation of FP resonances in central slits and the results are also presented in Table 1.
In comparison with the situation of illuminating the lens by a plane wave, when the lens is illuminated by a Gaussian beam, the central slits provide more phase retardation and hence the focal length decreases. Also the amount of decrease in focal length depends on the beam waist of the Gaussian beam which is shown in Fig. 5.
The performance parameters of the lens for the cases of illuminating the structure by plane wave, Gaussian beam with beam waist of 2 μm and 1.5 μm are listed in Table 2.
We are limited in how much we can increase the incident
light intensity to reduce the focal length, but we can increase the maximum change in focal length by employing a Gaussian beam. Our simulations reveal that when the structure is illuminated by 5.3 GW/cm^{2} Gaussian beam with 1.5 μm beam waist and plane wave the focal length will be 515 and 660 nm respectively, so if illuminating the structure by 13.2
MW/cm^{2} plane wave be considered as the initial situation, the maximum change in focal length increases 64% in the case of employing a Gaussian beam instead of a plane wave as is shown in Fig. 6.
Time？average electric field intensity distribution E^{2} for two cases depicted in Fig. 6, are presented in Fig. 7.
The maximum intensity of incident light in our simulations is 5.3 GW/cm^{2} which is below the threshold that can result in material damage, but temperature increase as the result of optical absorption, should be considered. Electron thermal dynamics induced by the interaction of a high intensity light with a metal can be precisely described using the TwoTemperature Model. Two temperatures refer to that of the electrons (
T_{e} ) and the lattice(T_{l} ), and are functions of time and depth into the sample. The lattice thermalizes with the electron but because of the large value of lattice heat capacity, ~2.4 J/cm^{3}/C° for silver, changes of hundreds of degree toT_{e} lead to changes of only tens of degrees toT_{l} . In experimental cases, by heat transfer to the substrate and cooling the lattice, this temperature rise will be lower [52, 53].IV. CONCLUSION
In conclusion, the performance of a plasmonic lens composed of an array of nanoslits introduced in thin metallic film with slanted cuts on the output surface and filled with Kerr nonlinear material has been explored. The performance parameters of the lens for the cases of illuminating by plane wave and Gaussian beam with different intensity and beam waist have been investigated. An exact formula for obtaining electric field
from electric flux densityE in Kerrdispersive medium which provide more accuracy than conventional ittirative method has been proposed. The maximum focusing efficiency of 53.5% has been achieved for the situation that the structure is illuminated by a plane wave with the intensity of 5.3 GW/cm^{2}. The focal length can be tuned on the xaxis of the lens by amount of 370 nm. According to our knowledge there are not any structural damages for the values of incident light intensity up to 5.3 GW/cm^{2}. Our proposed plasmonic lens offers great applications at many regions like near field imaging, sensing, subwavelength optics and etc.D

[FIG. 1.] (a) Schematic and (b) phase front profile of the proposed plasmonic lens.

[FIG. 2.] (a) Relative phase retardation caused by varying and (b) propagation length of SPPs as a function of the dielectric constant of incorporated material in nanoslits.

[FIG. 3.] Calculated normalized electric field intensity on the (a) xaxis (b) focal plane of the lens for the 13.2 MW/cm2 intensity of incident plane wave.

[FIG. 4.] (a) Decrease in focal length by increasing the intensity of incident plane wave and (b) the dependence of the dielectric constant of incorporated Kerr nonlinear material in nanoslits located at y=0 and 400nm to the incident light intensity.

[TABLE 1.] Performance parameters of the proposed lens for different values of incident plane wave intensity

[FIG. 5.] Calculated normalized electric field intensity on the xaxis of the lens for the cases of illuminating the lens with plane wave, Gaussian beam with wa = 2 μm and wa = 1.5 μm.

[FIG. 6.] 370 nm shift in focal point for the cases of illuminating the lens with 13.2 MW/cm2 plane wave and 5.3 GW/cm2 Gaussian beam with wa=1.5 μm.

[TABLE 2.] Performance parameters of the proposed lens for the cases of illuminating the lens with plane wave, Gaussian beam with wa = 2 μm and wa = 1.5 μm

[FIG. 7.] Time average electric field intensity distribution for illuming the lens by (a) 13.2 MW/cm2 plane wave and (b) 5.3 MW/cm2 Gaussian beam with wa=1.5 μm.