Increasing Surface Plasmons Propagation via Photonic Nanojets with Periodically Spaced 3 D Dielectric Cuboids

A structure based on periodically arranged 3D dielectric cuboids connected by photonic nanojets (PNJs) is proposed with the aim of increasing the propagation distance of surface plasmon polaritons (SPPs) at the telecom wavelength of 1550 nm. The performance of the structure is evaluated and compared with the case without the cuboids demonstrating that the SPPs propagation length is enhanced by a factor greater than 2, reaching a value of approximately 19λ0, when the gap between the cuboids is 2.5λ0. Also, the dependence of the propagation length with the height of the cubes is evaluated, showing that this parameter is critical for a good performance of the chain. A subwavelength resolution is obtained for all the jets generated at the output of the cuboids.


Introduction
The intrinsic limit imposed by the diffraction of the electromagnetic waves has been intensively studied by the scientific community in order to improve the performance of microscopy applications [1].Different techniques have been reported to surpass this limit, using for instance metamaterials [2,3], lenses [4][5][6][7] and microspherical particles [8,9].
Within this context, photonic nanojets (PNJs) were reported and experimentally demonstrated in the last decade with the aim to achieve a resolution below the diffraction limit [10][11][12].A PNJ is a region with a strong electric field concentration just at the output surface of micro-scaled dielectric particles (typically spheres or cylinders) illuminated with a plane wave.It has been demonstrated that the ratio between the refractive index of the dielectric and the background medium should be <2 in order to excite a PNJs with a subwavelength focus [12].Based on this, different configurations have been reported to improve the resolution of PNJs such as non-spherical particles [13], Luneburg lenses [14,15] and multilayered elements [16].Demonstrating that these structures can be used in applications such as nanoscopes and optical switching [8,17].
Recently, 3D dielectric cuboids have been proposed and experimentally demonstrated as an alternative structure to produce PNJs at sub-THz frequencies [18,19].It has been demonstrated that the same refractive index ratio constraint applies to cuboids in order to produce jets at its output.Also 3D dielectric cuboids have been studied working in the reflection mode regime, demonstrating that the PNJs can also be excited with this configuration [20].Additionally, it has been reported that this structure can be used as a waveguide by simply implementing a periodic array of 3D dielectric cuboids [21].Interestingly, PNJs have also been proposed at optical wavelengths using surface plasmon polaritons (SPPs) instead of a plane wave as an excitation source [22,23], demonstrating that the location of the PNJ can be tailored by changing the height of the dielectric particles (cuboids or cylinders).
In this communication, we report a plasmonic waveguide consisting of periodically spaced 3D dielectric cuboids on top of a metal film excited by SPPs, able to almost double the propagation distance of SPPs at the telecommunication wavelength of λ 0 = 1550 nm (which is a standard wavelength for optical applications).The structure is evaluated by periodically placing a total number of 6 cuboids and changing the gap between them.A study of the behavior when the height of the cuboids is changed from 80 nm to 240 nm is carried out for different distances between the cuboids, showing that the best performance in terms of the propagation distance of SPPs is achieved when the height of the cuboids is 120 nm and the gap between them is 2.5λ 0 .For this optimum case, the propagation of SPPs can be extended up to approximately 10λ 0 compared with the case without cuboids.Moreover, it is shown that a subwavelength focus is achieved for all the cuboids, demonstrating that the propagation distance of SPPs can be extended without deteriorating the performance of PNJs.

Design
The chain of 3D dielectric cuboids considered in this work is schematically shown in Figure 1.It consists of an array of 6 cuboids with dimensions d x ˆdy ˆdz separated by a gap p placed on top of a metal layer of gold.A dispersive model is used for gold following the Johnson and Christy experimental data [24].The height of the metal layer is 100 nm, i.e., well above the skin depth at 1550 nm which is below 30 nm, so that it can be considered semi-infinite at the design wavelength [25].The material of the cuboids is Silicon Nitride Si 3 N 4 with has a refractive index n d = 1.97.Note that a fixed value has been selected instead of a dispersive model to simplify the computation (this can be done because the variation of n d is very small-from 1.97 to 1.99-within the spectral band of 750-2000 nm [23]).The whole structure is immersed in air with a refractive index n 0 = 1.With this configuration, the SPPs are excited from the back of the chain and propagate along the z-axis.The numerical simulations along this manuscript are done using the transient solver of the commercial software CST Microwave Studio™.As it is known, there are several techniques to excite SPPs, such as dipoles with different polarizations [26,27].In our case we use a waveguide placed at the back of the structure with dimensions 5 µm ˆ0.15 µm along the xand zaxis, respectively [23].The total dimensions of the simulation space are 5 µm ˆ41 µm along the xand zaxis, respectively.An extra air region of 1 µm was added at the top of the structure and open boundary conditions were used for the entire simulation domain.With this condition, the air domain at the top of the structure is considered as semi-infinite, enabling to avoid undesirable reflections.
cuboids [21].Interestingly, PNJs have also been proposed at optical wavelengths using surface plasmon polaritons (SPPs) instead of a plane wave as an excitation source [22,23], demonstrating that the location of the PNJ can be tailored by changing the height of the dielectric particles (cuboids or cylinders).
In this communication, we report a plasmonic waveguide consisting of periodically spaced 3D dielectric cuboids on top of a metal film excited by SPPs, able to almost double the propagation distance of SPPs at the telecommunication wavelength of λ0 = 1550 nm (which is a standard wavelength for optical applications).The structure is evaluated by periodically placing a total number of 6 cuboids and changing the gap between them.A study of the behavior when the height of the cuboids is changed from 80 nm to 240 nm is carried out for different distances between the cuboids, showing that the best performance in terms of the propagation distance of SPPs is achieved when the height of the cuboids is 120 nm and the gap between them is 2.5λ0.For this optimum case, the propagation of SPPs can be extended up to approximately 10λ0 compared with the case without cuboids.Moreover, it is shown that a subwavelength focus is achieved for all the cuboids, demonstrating that the propagation distance of SPPs can be extended without deteriorating the performance of PNJs.

Design
The chain of 3D dielectric cuboids considered in this work is schematically shown in Figure 1.It consists of an array of 6 cuboids with dimensions dx × dy × dz separated by a gap p placed on top of a metal layer of gold.A dispersive model is used for gold following the Johnson and Christy experimental data [24].The height of the metal layer is 100 nm, i.e., well above the skin depth at 1550 nm which is below 30 nm, so that it can be considered semi-infinite at the design wavelength [25].The material of the cuboids is Silicon Nitride Si3N4 with has a refractive index nd = 1.97.Note that a fixed value has been selected instead of a dispersive model to simplify the computation (this can be done because the variation of nd is very small-from 1.97 to 1.99-within the spectral band of 750-2000 nm [23]).The whole structure is immersed in air with a refractive index n0 = 1.With this configuration, the SPPs are excited from the back of the chain and propagate along the z-axis.The numerical simulations along this manuscript are done using the transient solver of the commercial software CST Microwave Studio™.As it is known, there are several techniques to excite SPPs, such as dipoles with different polarizations [26,27].In our case we use a waveguide placed at the back of the structure with dimensions 5 μm × 0.15 μm along the x-and z-axis, respectively [23].The total dimensions of the simulation space are 5 μm × 41 μm along the x-and z-axis, respectively.An extra air region of 1 μm was added at the top of the structure and open boundary conditions were used for the entire simulation domain.With this condition, the air domain at the top of the structure is considered as semi-infinite, enabling to avoid undesirable reflections.As it can be observed in Figure 1b, two different regions are recognized in the chain, namely, with and without cuboids, which correspond to an air-dielectric-metal and an air-metal multilayer As it can be observed in Figure 1b, two different regions are recognized in the chain, namely, with and without cuboids, which correspond to an air-dielectric-metal and an air-metal multilayer structures, respectively.In both of them, the original multilayer problem can be reduced to an effective refractive index applying an homogenization technique, as described in [22,23].
Photonics 2016, 3, 10 3 of 7 Following this procedure, we find that the effective index of the region without cuboids is very similar to the index of air (~1.0043) at 1550 nm, indicating that SPPs are not strongly coupled at the design wavelength.This is corroborated in the simulation results of Figure 2 where we have represented the E y -field distribution on the metal surface.It can be observed that SPPs are weakly coupled to the surface of the metal and hence decay rapidly at 1550 nm, reaching a short propagation length (i.e., distance at which the magnitude of the electric field decays by a factor 1/e) of ~9λ 0 , see green curve in Figure 3h.
Our aim in this work is to increase the propagation length by means of a chain of PNJ cuboids.In the next section we will perform a parametrical study to find the distance between cuboids as well as the cuboid height that optimizes the propagation length and the full-width at half-maximum along the transversal x-axis (FWHM x ).
Photonics 2016, 3, 10 3 of 7 structures, respectively.In both of them, the original multilayer problem can be reduced to an effective refractive index applying an homogenization technique, as described in [22,23].Following this procedure, we find that the effective index of the region without cuboids is very similar to the index of air (~1.0043) at 1550 nm, indicating that SPPs are not strongly coupled at the design wavelength.This is corroborated in the simulation results of Figure 2 where we have represented the Ey-field distribution on the metal surface.It can be observed that SPPs are weakly coupled to the surface of the metal and hence decay rapidly at 1550 nm, reaching a short propagation length (i.e., distance at which the magnitude of the electric field decays by a factor 1/e) of ~9λ0, see green curve in Figure 3h.
Our aim in this work is to increase the propagation length by means of a chain of PNJ cuboids.In the next section we will perform a parametrical study to find the distance between cuboids as well as the cuboid height that optimizes the propagation length and the full-width at half-maximum along the transversal x-axis (FWHMx).

Results and Discussion
First, we evaluate the performance of the dielectric chain as the height of the 3D cuboids and the period p between them are varied.From our previous works [18,23] we know that the optimal lateral dimensions of the cuboids should be equal to the operation wavelength, dx = dz = λ0 = 1550 nm in order to produce the PNJ at the output face of the dielectric cuboids.With these dimensions, we evaluate the performance of the chain for p ranging from 0 to 3λ0 with a step of 0.5λ0 and a height of cuboids dy varying from 80 nm to 240 nm with a step of 20 nm.As it has been demonstrated in previous works, this dimension is critical in the performance of the PNJ, as it directly allows tuning the effective index of the region where the cuboid is present [22,23].
Contour plots showing the magnitude of the Ey field at the design wavelength as a function of the propagation z-axis and height of the dielectric cuboids are shown in Figure 3a-g for each value of p under study.All the results have been normalized to the incident electric field amplitude at the back of the structure.This figure gives us all the information to evaluate the performance of the chain.Let us first study the influence of the period between the cuboids.As shown in the figure, the worst case is observed when p = 0 where the decay of the SPPs is stronger compared with the other values of p.This is an expected result, as this period corresponds to the case when the cuboids are touching and hence only one PNJ is excited at the end of the structure (z ~ 7λ0).As the cuboids are taken far apart, several peaks of high electric field concentration appear between cuboids.These peaks are the PNJs excited by each cuboid.The propagation length increases as p is increased due to the PNJ excited by each cuboid produces an enhancement of the electric field at the focal length [23].Note that this propagation length increases until it reaches a maximum at p = 2.5λ0, and from there on it decreases again.The best case is when p = 2.5λ0 and dy = 120 nm (highlighted with a dashed white line in Figure 3f).We have plotted in Figure 3h the magnitude of the Ey field along the propagation axis for this optimum case compared to a flat metal film.The SPP propagation length is 19λ0, i.e., it is enhanced by a factor greater than two compared to the case without the cuboids (which reached only ~9λ0).

Results and Discussion
First, we evaluate the performance of the dielectric chain as the height of the 3D cuboids and the period p between them are varied.From our previous works [18,23] we know that the optimal lateral dimensions of the cuboids should be equal to the operation wavelength, d x = d z = λ 0 = 1550 nm in order to produce the PNJ at the output face of the dielectric cuboids.With these dimensions, we evaluate the performance of the chain for p ranging from 0 to 3λ 0 with a step of 0.5λ 0 and a height of cuboids d y varying from 80 nm to 240 nm with a step of 20 nm.As it has been demonstrated in previous works, this dimension is critical in the performance of the PNJ, as it directly allows tuning the effective index of the region where the cuboid is present [22,23].
Contour plots showing the magnitude of the E y field at the design wavelength as a function of the propagation z-axis and height of the dielectric cuboids are shown in Figure 3a-g for each value of p under study.All the results have been normalized to the incident electric field amplitude at the back of the structure.This figure gives us all the information to evaluate the performance of the chain.Let us first study the influence of the period between the cuboids.As shown in the figure, the worst case is observed when p = 0 where the decay of the SPPs is stronger compared with the other values of p.This is an expected result, as this period corresponds to the case when the cuboids are touching and hence only one PNJ is excited at the end of the structure (z ~7λ 0 ).As the cuboids are taken far apart, several peaks of high electric field concentration appear between cuboids.These peaks are the PNJs excited by each cuboid.The propagation length increases as p is increased due to the PNJ excited by each cuboid produces an enhancement of the electric field at the focal length [23].Note that this propagation length increases until it reaches a maximum at p = 2.5λ 0 , and from there on it decreases again.The best case is when p = 2.5λ 0 and d y = 120 nm (highlighted with a dashed white line in Figure 3f).We have plotted in Figure 3h the magnitude of the E y field along the propagation axis for this optimum case compared to a flat metal film.The SPP propagation length is ~19λ 0 , i.e., it is enhanced by a factor greater than two compared to the case without the cuboids (which reached only ~9λ 0 ).Regarding the height of the cuboids, it is interesting to observe that in all cases the longest SPPs propagation lengths are obtained for values of dy between approximately 100 and 140 nm.This performance can be easily explained in terms of the effective index, neff.For this analysis we concentrate on p = 2.5λ0, which can be considered representative of all the studied cases.We have calculated the power distribution on the surface of the metal film for different values of dy, see Figure 4a-e.Then we have extracted the propagation length for all the dy values considered and have plotted the results in Figure 4f along with the analytical values of the effective index contrast between the regions with and without cuboids [23].For small height values (dy = 80; 100 nm) the effective index contrast is close to 1 and hence, the PNJs are produced far away from the output surface of the cuboids [18,23], giving as a result a poor illumination of the neighboring cuboid in the chain.On the other hand, for large height values (dy > 120 nm) the index contrast increases as dy is increased.This implies that the PNJs are produced inside of the cuboids, as observed in Figure 4d,e and also discussed in [18,23].Therefore, the excitation of consecutive cuboids is not efficient and the propagation length is reduced.(e) p = 2λ 0 ; (f) p = 2.5λ 0 and (g) p = 3λ 0 .(h) Numerical results of the normalized magnitude of the electric field along the propagation axis at the surface of the metal for the optimal case, p = 2.5λ 0 and d y = 120 nm (extracted from the dashed line shown in panel (f)) (blue curve) and a flat metal film (green curve).The red curve corresponds to a decay of 1/e with respect to the amplitude of the incident surface plasmon polaritons (SPP) at the input.The 3D cuboids are also shown as grey blocks.
Regarding the height of the cuboids, it is interesting to observe that in all cases the longest SPPs propagation lengths are obtained for values of d y between approximately 100 and 140 nm.This performance can be easily explained in terms of the effective index, n e f f .For this analysis we concentrate on p = 2.5λ 0 , which can be considered representative of all the studied cases.We have calculated the power distribution on the surface of the metal film for different values of d y , see Figure 4a-e.Then we have extracted the propagation length for all the d y values considered and have plotted the results in Figure 4f along with the analytical values of the effective index contrast between the regions with and without cuboids [23].For small height values (d y = 80; 100 nm) the effective index contrast is close to 1 and hence, the PNJs are produced far away from the output surface of the cuboids [18,23], giving as a result a poor illumination of the neighboring cuboid in the chain.On the other hand, for large height values (d y > 120 nm) the index contrast increases as d y is increased.This implies that the PNJs are produced inside of the cuboids, as observed in Figure 4d,e and also discussed in [18,23].Therefore, the excitation of consecutive cuboids is not efficient and the propagation length is reduced.
The longest propagation length is when d y = 120 nm, in which the effective index contrast is ~1.2.Note that with this value the PNJs are produced very near but slightly away from output surface of the cuboids (approximately at 0.16λ 0 ) as can be observed in Figure 4b.Since the PNJs are not just at the output surface of the cuboids, it allows obtaining the optimal illumination of the next cuboids of the chain.Note that this value of d y is different from the optimal case we found in our previous study [23].Recall that in that paper the aim was to produce the PNJ just at the output surface of the cuboids and not to extend the SPP propagation.Anyway, it is interesting to highlight that the SPP propagation distance is enhanced for a wide range of values of d y .This is an important feature in terms of the allowed tolerances for fabrication purposes.If we take 120 nm as a reference, the height of the cuboids may be fabricated with a tolerance of ˘40 nm without strongly altering the SPP propagation distance.
Photonics 2016, 3, 10 5 of 7 The longest propagation length is when dy = 120 nm, in which the effective index contrast is ~1.2.Note that with this value the PNJs are produced very near but slightly away from output surface of the cuboids (approximately at 0.16λ0) as can be observed in Figure 4b.Since the PNJs are not just at the output surface of the cuboids, it allows obtaining the optimal illumination of the next cuboids of the chain.Note that this value of dy is different from the optimal case we found in our previous study [23].Recall that in that paper the aim was to produce the PNJ just at the output surface of the cuboids and not to extend the SPP propagation.Anyway, it is interesting to highlight that the SPP propagation distance is enhanced for a wide range of values of dy.This is an important feature in terms of the allowed tolerances for fabrication purposes.If we take 120 nm as a reference, the height of the cuboids may be fabricated with a tolerance of ±40 nm without strongly altering the SPP propagation distance.For the sake of completeness, we also evaluate the PNJ performance in terms of the transversal resolution.We have represented in Figure 5 the full-width at half-maximum along the x-axis (FWHMx, defined at the distance along the x-axis at which the power decays to half of its maximum) of the PNJs produced by all the cuboids in the array, for the optimal case under study.It can be observed that a subwavelength focus is obtained in all cases with a resolution below 0.75λ0.

Conclusions
In this work, a dielectric chain consisting of periodically spaced 3D dielectric cuboids has been thoroughly studied, with the aim to increase the propagation distance of SPPs at the telecom For the sake of completeness, we also evaluate the PNJ performance in terms of the transversal resolution.We have represented in Figure 5 the full-width at half-maximum along the x-axis (FWHM x , defined at the distance along the x-axis at which the power decays to half of its maximum) of the PNJs produced by all the cuboids in the array, for the optimal case under study.It can be observed that a subwavelength focus is obtained in all cases with a resolution below 0.75λ 0 .
Photonics 2016, 3, 10 5 of 7 The longest propagation length is when dy = 120 nm, in which the effective index contrast is ~1.2.Note that with this value the PNJs are produced very near but slightly away from output surface of the cuboids (approximately at 0.16λ0) as can be observed in Figure 4b.Since the PNJs are not just at the output surface of the cuboids, it allows obtaining the optimal illumination of the next cuboids of the chain.Note that this value of dy is different from the optimal case we found in our previous study [23].Recall that in that paper the aim was to produce the PNJ just at the output surface of the cuboids and not to extend the SPP propagation.Anyway, it is interesting to highlight that the SPP propagation distance is enhanced for a wide range of values of dy.This is an important feature in terms of the allowed tolerances for fabrication purposes.If we take 120 nm as a reference, the height of the cuboids may be fabricated with a tolerance of ±40 nm without strongly altering the SPP propagation distance.For the sake of completeness, we also evaluate the PNJ performance in terms of the transversal resolution.We have represented in Figure 5 the full-width at half-maximum along the x-axis (FWHMx, defined at the distance along the x-axis at which the power decays to half of its maximum) of the PNJs produced by all the cuboids in the array, for the optimal case under study.It can be observed that a subwavelength focus is obtained in all cases with a resolution below 0.75λ0.

Conclusions
In this work, a dielectric chain consisting of periodically spaced 3D dielectric cuboids has been thoroughly studied, with the aim to increase the propagation distance of SPPs at the telecom

Conclusions
In this work, a dielectric chain consisting of periodically spaced 3D dielectric cuboids has been thoroughly studied, with the aim to increase the propagation distance of SPPs at the telecom wavelength of 1550 nm.It has been shown that the propagation distance can be increased by a factor greater than two (up to ~19λ 0 ) for the optimal case when the gap between them is 2.5λ 0 and the height of the cuboids is 140 nm.It has been shown that for heights >180 nm, the propagation distance of the SPPs is strongly reduced which is a direct consequence of the ratio between the effective refractive indexes of the air-metal-dielectric and air-metal regions of the structure.Moreover, it has been demonstrated that a subwavelength transversal resolution of the PNJs can be achieved at the output of each cuboid of the array.The results here presented may be useful in the design of different photonic applications such as beam steerers and SPP optical tweezers, which are now under study.

Figure 1 .
Figure 1.Schematic representation of the cuboid chain: (a) perspective and (b) lateral view.The dielectric cuboids have dimensions dx × dy × dz and are periodically placed along the z-axis with an air gap p between them.

Figure 1 .
Figure 1.Schematic representation of the cuboid chain: (a) perspective and (b) lateral view.The dielectric cuboids have dimensions d x ˆdy ˆdz and are periodically placed along the z-axis with an air gap p between them.

Figure 2 .
Figure 2. Ey-field distribution on the xz-plane (a) and (yz-plane) (b) for a flat gold film at 1550 nm.

Figure 2 .
Figure 2. E y -field distribution on the xz-plane (a) and (yz-plane) (b) for a flat gold film at 1550 nm.

Figure 3 .
Figure 3. Normalized magnitude of the Ey field as a function of the propagation z-axis and height dy of the cuboids when the distance between cuboids is: (a) p = 0; (b) p = 0.5λ0; (c) p = λ0; (d) p = 1.5λ0;(e) p = 2λ0; (f) p = 2.5λ0 and (g) p = 3λ0.(h) Numerical results of the normalized magnitude of the electric field along the propagation axis at the surface of the metal for the optimal case, p = 2.5λ0 and dy = 120 nm (extracted from the dashed line shown in panel (f)) (blue curve) and a flat metal film (green curve).The red curve corresponds to a decay of 1/e with respect to the amplitude of the incident surface plasmon polaritons (SPP) at the input.The 3D cuboids are also shown as grey blocks.

Figure 3 .
Figure 3. Normalized magnitude of the E y field as a function of the propagation z-axis and height d y of the cuboids when the distance between cuboids is: (a) p = 0; (b) p = 0.5λ 0 ; (c) p = λ 0 ; (d) p = 1.5λ 0 ;(e) p = 2λ 0 ; (f) p = 2.5λ 0 and (g) p = 3λ 0 .(h) Numerical results of the normalized magnitude of the electric field along the propagation axis at the surface of the metal for the optimal case, p = 2.5λ 0 and d y = 120 nm (extracted from the dashed line shown in panel (f)) (blue curve) and a flat metal film (green curve).The red curve corresponds to a decay of 1/e with respect to the amplitude of the incident surface plasmon polaritons (SPP) at the input.The 3D cuboids are also shown as grey blocks.

Figure 4 .
Figure 4. Power distribution on the xz-plane for the case of 3D dielectric cuboids placed periodically with an air gap of p = 2.5λ0 and height: (a) dy = 80 nm; (b) dy = 120 nm; (c) dy = 160 nm; (d) dy = 200 nm and (e) dy = 240 nm.(f) Numerical results of SPP propagation distance (black curve) and index contrast (blue curve) as a function of the cuboid height.

Figure 5 .
Figure 5. Numerical results of the full-width at half-maximum along the transversal x-axis of the photonic nanojets (PNJ) produced by each cuboid for the optimal case when p = 2.5λ0 and dy = 120 nm.

Figure 4 .
Figure 4. Power distribution on the xz-plane for the case of 3D dielectric cuboids placed periodically with an air gap of p = 2.5λ 0 and height: (a) d y = 80 nm; (b) d y = 120 nm; (c) d y = 160 nm; (d) d y = 200 nm and (e) d y = 240 nm.(f) Numerical results of SPP propagation distance (black curve) and index contrast (blue curve) as a function of the cuboid height.

Figure 4 .
Figure 4. Power distribution on the xz-plane for the case of 3D dielectric cuboids placed periodically with an air gap of p = 2.5λ0 and height: (a) dy = 80 nm; (b) dy = 120 nm; (c) dy = 160 nm; (d) dy = 200 nm and (e) dy = 240 nm.(f) Numerical results of SPP propagation distance (black curve) and index contrast (blue curve) as a function of the cuboid height.

Figure 5 .
Figure 5. Numerical results of the full-width at half-maximum along the transversal x-axis of the photonic nanojets (PNJ) produced by each cuboid for the optimal case when p = 2.5λ0 and dy = 120 nm.

Figure 5 .
Figure 5. Numerical results of the full-width at half-maximum along the transversal x-axis of the photonic nanojets (PNJ) produced by each cuboid for the optimal case when p = 2.5λ 0 and d y = 120 nm.