Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Plasmonic array nanoantennas on layered substrates: modeling and radiation characteristics

Open Access Open Access

Abstract

In this paper, we theoretically characterize the performance of array of plasmonic core-shell nano-radiators located over layered substrates. Engineered substrates are investigated to manipulate the radiation performance of nanoantennas. A rigorous analytical approach for the problem in hand is developed by applying Green’s function analysis of dipoles located above layered materials. It is illustrated that around the electric scattering resonances of the subwavelength spherical particles, each particle can be viewed as an induced electric dipole which is related to the total electric field upon that particle by a polarizability factor. Utilizing this, we can effectively study the physical performance of such structures. The accuracy of our theoretical model is validated through using a full-wave finite difference time domain (FDTD) numerical technique. It is established that by novel arraying of nano-particles and tailoring their multilayer substrates, one can successfully engineer the radiation patterns and beam angles. Several optical nanoantennas designed on layered substrates are explored. Using the FDTD the effect of finite size substrate is also explored.

©2009 Optical Society of America

1. Introduction

Recently, the field of optical nano-metals and plasmonics has experienced an extensive growth in research activities, owing to the fascinating possibility of confining electromagnetic radiation to subwavelength spatial domains via surface plasmon polaritons (SPPs) [1]–[6]. These surface modes are established by the couplings of electromagnetic fields to conduction electrons at the boundary between a noble metal and a dielectric. This phenomenon offers unique opportunity for nanooptics physics enabling, areas such as plasmonic waveguides, nano-interconnects and plasmonic light sources [1]. Other newly proposed applications include plasmonic THz nano-antennas [7]–[14] and solar-cells energy-efficient guiding [15]–[19]. Patterning plasmonic particles can offer significant advantage for achieving desired optical properties.

One of the interesting applications of the plasmonic particles is to build antenna devices in optical frequencies featuring photonic communication. However, design and fabrication of nanoantennas is the most challenging part mainly due to the optical properties of metal. For instance, at optical frequencies metal is well described by a frequency dispersive negative permittivity material rather than a high-conductivity structure as what is observed in microwave spectrum. As a consequent, one can offer resonant nanoparticles using core-shell dielectric-plasmonic elements enabling nanoantenna radiation [14]. The goal of this article is to theoretically exploit the performance of arrays of plasmonic nanoparticles featuring enhanced radiation characteristics when placed over layered substrates.

At microwave frequencies, in order to achieve a high gain antenna performance, one needs to use the concept of array antennas where different radiators are tailored in unique arrangements to control the fields amplitudes and phases in the far-field spectrum enabling directive emission. This concept can be realized in optics with the use of plasmonic particles acting like small dipole radiators. Li et al. [7] used this concept for providing a working THz Yagi-Uda antenna in free-space. Similar performance is observed in [14]. From practical point of view, however, one needs to deposit the nano radiators on a substrate. In this case, the radiation performance changes drastically as for instance illustrated for an array of plasmonic rods located above a half-space medium through a numerical analysis [11].

It can be of great benefit if one establishes a theoretical formulation for performance analysis of nanoantennas composed of an array of core-shell spherical particles located above layered substrates. This will allow successful study of the complex structure via understanding the physics and obtaining the optimal radiation behavior and beam scanning for the application of interest.

To establish the theory of nanoantennas, array of plasmonic core-shells over a layered material is approximated with electric dipoles where their induced dipoles are calculated through the Green’s function analysis of dipoles over the layered material [20]–[22]. The obtained Sommerfeld integrals are efficiently calculated by defining a robust integral path [20]. We apply our theory-modeling to obtain the radiation performance of a dipole exciting an array of plasmonic nanoantennas located above a planar layered material. The induced dipole on each nano coreshell is the key parameter determining the antenna radiation characteristics. A full-wave finite difference time domain (FDTD) technique [23] is also applied to characterize the performance of array of plasmonic nanoantennas on layered substrate, validating our theoretical model. The behavior of nanoantennas array on finite-size substrate is also explored.

Optimal patterning of plasmonic particles along with engineering their layered substrate can successfully tailor the optical parameters for achieving a directive emission. It is illustrated that a focused-beam can be obtained for a 2D array of plasmonic particles-nanoantenna deposited on a layered materials and arranged in a Yagi-Uda type configuration. We show that by engineering the layered substrate one can control the radiation beam for focusing it in the desired direction, or suppressing the back radiation.

2. Theory and modeling

The configuration of a layered substrate plasmonic particle-antenna is shown in Fig. 1. The geometry consists of subwavelength nano core-shells implanted above an infinitely long (in x and y directions) semiconductor substrate where layer(s) of plasmonic or dielectric materials are coating the semiconductor surfaces. Here, the plasmonic particles are spherical core-shells with dielectric cores (kc,εrc) and plasmonic coatings (ks,εrs). The use of plasmonic material as the coating allows greater control on scattering characteristics of the nano core-shells [7].

At the plasmonic resonance the scattered field from a small spherical core-shell under the influence of an incident field is dominated by dipole terms. Thus, a subwavelength particle at the scattering resonance can be viewed as an induced dipole (p), which is related to the total electric field upon that particle by a polarizability factor (α).

 figure: Fig. 1.

Fig. 1. Array of plasmonic nanoantennas located above a layered substrate.

Download Full Size | PDF

p=α.Etotal=i6πε0k03Γ1(e).Etotal,

where k 0 is the wave number of the surrounding medium and, Γ(e) 1 is the first electric scattering coefficient of the plasmonic core-shell [7, 24]:

Γn(e)=Un(e)Un(e)+iVn(e),

U (e) n and V (e) n are given by:

Un(e)=jn(kcr1)jn(ksr1)yn(ksr1)0jn(k˜cr1)εcjn(k˜sr1)εsyn(k˜sr1)εs00jn(ksr2)yn(ksr2)jn(k0r2)0j˜n(ksr2)y˜n(ksrs)εsj˜n(k0r2)ε0,
Vn(e)=jn(kcr1)jn(ksr1)yn(ksr1)0jn(k˜cr1)εcjn(k˜sr1)εsyn(k˜sr1)εp00jn(ksr2)yn(ksr2)yn(k0r2)0j˜n(ksr2)y˜n(ksrs)εsy˜n(k0r2)ε0.

where r 1 and r 2 are the inner and outer radii respectively. It is well known that, the scattering resonance of a plasmonic concentric core-shell structure can be obtained at desired frequencies by tuning the (r 1/r 2) radii ratio [7]. Figures. 2(a) and (b) present the magnitudes and phases of normalized polarizability versus the r 1/r 2 ratio for a concentric core-shell with silver coating and dielectric core made of SiO 2, wherein the total radius r 2/0.1λ 0 is taken to be constant and the core’s radius (r 1) varies. As it is shown, by tuning the r 1/r 2 ratio, we can tailor the resonance frequency of the plasmonic nano-particle, for instance, the nano-sphere goes to resonance at f=650 THz0=461 nm) if r 1/r 2=0.71. Another interesting feature in these figures, is that at the resonance the phase of polarizability is close to π/2 (phase of i), which means that at the resonance the induced dipole is almost in phase with the incident field. Notice that, to tune the resonant frequency for a simple metallic nano-particle, one should change the total size of the particle, while in many applications it is required to maintain a desired size. Nevertheless, with the use of core-shell structures, the total size can be remained the same, whereas at the same time the resonant frequency can be tuned by changing the radii ratio.

 figure: Fig. 2.

Fig. 2. Magnitude and phase of normalized polarizability of a plasmonic core-shell sphere versus r 1/r 2 ratio for different frequencies for a structure with r 2=0.1λ, εrc=2.2+0.01i and silver coating. (a) Polarization magnitude, and (b) Polarization phase, (c) Silver permittivity behavior.

Download Full Size | PDF

In order to analyze the frequency dependency of the scattering resonance, it is of value to consider the frequency dispersion of the plasmonic material. The permittivity of the silver in optical regime can be described with a Drude model, in which the numerical parameters are achieved by fitting the available experimental data of the Silver’s permittivity to the Drude model [25] (see also Fig. 2(c)).

εrs=11.513×1032ω(ω+i7.9×1013).

In addition, note that silver is the material of choice, due to its small loss in optical frequencies.

Now, let us investigate the case of Fig. 1. If a non-periodic array of plasmonic core-shells are placed over a layered material, with the presence of an incident field, each nano-particle can be viewed as an induced dipole around the scattering resonance of that nano core-shell. From (1) it is clear that, the induced dipole on each nano core-shell is proportional to the total field upon that particle. In this scenario, the total field upon each nano core-shell can be expressed as the summation of three terms. The first part is associated with the total incident field in the absence of the nano core-shells (E totalinc), the second part is the electric field due to the couplings between the nano core-shells in the absence of the layered material (G̿l dipole(r l,r q)p q). Since this term represents the couplings between the nano core-shells, for the lth nano core-shell we consider the fields of all other nano core-shell except itself (excluding the field of the lth particle). The last term, is associated with the reflected fields from the layered substrate (G̿lreflected(r l,r q)p q). Note that for the computing the last two terms we approximate each nano core-shell with an electric dipole. Hence for the second term we calculate the dipolar couplings and the Green’s function analysis of dipoles over layered material is applied for evaluation of the final part [20]–[22]. Also, it is worth mentioning that the fields associated with every nano core-shells (both couplings and reflected field) is directly proportional to the induce dipole moment, thus the induced dipole moment for each particle is derived by solving the following linear system of equations. For l,q ∊ {1,2; …,N} with N being the total number of particles, we obtain:

pl=αl(Einctotal(rl)+q,qlG̿dipolel(rl,rq)pq+qG̿reflectedl(rl,rq)pq),

where Etotalinc denotes the sum of the incident field and its reflection from the layered material in the absence of the nano-spheres. G̿l dipole is the dyadic Green’s function of the qth nano core-shell evaluated at the position of the lth particle. G̿l reflected is the reflected Greens function of the qth nano-sphere (from the layered material) computed at the location of the lth one, i.e.,

G̿dipolel(rl,rq)=[(k12+22x)2xy2xz2xy(k12+22y)2yz2xz2yz(k12+22z)]eik1rlrq4πε1rlrq,
G̿reflectedl(rlrq)=[Gerxlx(rlrq)Gerylx(rlrq)Gerzlx(rlrq)Gerxly(rlrq)Geryly(rlrq)Gerzly(rlrq)Gerxlz(rlrq)Gerylz(rlrq)Gerzlz(rlrq)].

By defining r=r l-r q, the elements of the matrix in equation (6b) are derived as [20, 21]:

Gerxlz=18πε1cosϕdkρkρ2H1(1)(kρρ)R˜1,2TMeik1z(z+2d1),
Gerylz=18πε1sinϕdkρkρ2H1(1)(kρρ)R˜1,2TMeik1z(z+2d1),
Gerzlz=i8πε1dkρkρ3k1zH0(1)(kρρ)R˜1,2TMeik1z(z+2d1),

where for example, Gelzrx is Green’s function for the z-directed electric field associated with a x-directed dipole. Let G̃ be the kernel of the integral, for w ∊ {x,y,z} and s ∊ {x,y} we have [20]:

G˜erwls(kρ)=1kρ2[zG˜erwls(kρ)zμẑ×sG˜hrwlz(kρ)],

where Ghlzrw is the Green’s function for the z-directed magnetic field due to a w-directed dipole located above the layered material [20]–[22]. TM/TE 1,2 is the generalized reflection coefficient (that incorporates subsurface reflections) at the interface between region 1 and region 2 [20]. The generalized reflection coefficient at the interface of region i and i+1 is defined as [20]:

R˜i,i+1TMTE=Ri,i+1TMTE+R˜i+1,i+2TMTEe2iki+1,z(di+1di)1+Ri,i+1TMTER˜i+1,i+2TMTEe2iki+1,z(di+1di),
Ri,jTMTE=(εjμj)kiz(εiμi)kjz(εjμj)kiz+(εiμi)kjz,

In the above equations, ki is the wave number in the layer i and:

kiz=ki2kρ2.

The integrals in (7) and (8) are calculated by defining a robust numerical path whereby the discretization of the integration steps is chosen to adapt to the integrand’s oscillation [20]. For the verification purposes, the integrals in (7) and (8) are computed in a homogenous medium (ε=ε 0,µ=µ 0) where their exact values are given by the response of a dipole in the same homogenous regime. Figure .3 shows the maximum percentage of error between all the computed integrals (nine integrals) in (7) and (8) for each observation point marked as (ρ,z 0). As it is exhibited, the maximum error is less than 0:1% for different observation points, validating our method of integration.

 figure: Fig. 3.

Fig. 3. Maximum percentage of error in obtaining the integrals in (7) and (8) using our theoretical model compared to their exact values for different observation points. ρ and z0 denote the location of the observation point.

Download Full Size | PDF

It is worth noticing that, if instead of the array only one nano core-shell is implanted above a layered material then the induced dipole is achieved utilizing the following equation:

p=(I̿αG̿reflected)1αEtotal,

where G̿reflected in this particular case simplifies to a diagonal matrix, meaning that there is no coupling between the components of the induced dipole. It is worth highlighting again, that in this case the induced dipole is only related to the reflected field of that nano-particle from the layered material and the incident field.

As mentioned earlier, regardless of the structure, we assume that each nano core-shell can be modeled as an electrically small dipole with the induced dipole p. Yet, the induced dipole for each nano core-shell is derived by using (5), in which the parameters of the G̿reflected are determined by the geometry and materials of the layered structure. Moreover, the far-field radiation pattern is also affected by the presence of the substrate. Hence, to understand the behavior of the antenna located above a layered material, let us first analyze the radiation pattern of such antennas. Since each nano core-shell behaves like an electric dipole, the far-field radiation pattern of such a structure depicted in Fig. 1, consists of summation of the radiation patterns of electric dipoles located above a layered material.

The far zone electric field for a dipole deposited over a layered-substrate can be evaluated using the conventional steepest decent contour (SDC) technique with the transformation kρ=k sinθ, where the θ is the spherical angle from the z axis. So for each dipole (px,py,pz) located at (x 0,y 0,z 0), the upper half-space and lower half-space far-field radiation pattern can be represented as [26]:

E=EθEϕ=kJ24πεJeikJrr(pxcosϕ+pysinϕ)cosθΦJ2pzsinθΦJ1(pxsinϕpysinϕ)ΦJ3.

Where the index J ∊ [1;N] is to distinguish between the upper-half (ε 1,µ 1) and lower-half (εNN). The potential parameters are defined as [26]:

Φ11=eik1r[1+R˜1,2TM(θ)eik1z0],
Φ12=eik1r[1R˜1,2TM(θ)eik1z0],
Φ13=eik1r[1+R˜1,2TM(θ)eik1z0],
ΦN1=ΦN2=nNn1cosθs˜zTTM(θ)eikN[z0s˜z+(dNd1)cosθ],
ΦN3=cosθs˜zTTM(θ)eikN[z0s˜z+(dNd1)cosθ],

with n 1 and nN being the refractive indices for the upper and lower half-spaces respectively. The TTM/TE is the generalized TM/TE transmission coefficient for a layered material [20]. And kJr and kJz are specified as:

kJr=kJ[x0cosϕsinθ+y0sinϕsinθ+z0cosθ],
kJz=2kJ(z0+d1)cosθ,

also:

s˜z=k1zkN=(n1nN)2sin2θ.

Notice that, from (13) a remarkable result is achieved. The first term in the potentials (Φ1 J3 J) can be interpreted as the far zone radiation pattern of a dipole, where the second term can be identified as the radiation from a dipole located at the image plane (at the distance z0 beneath the d 1) weighted by generalized reflection coefficients. It is also worth mentioning that, at θ=π/2, k 1z=0 and hence RTM/TE 1,2 is zero, which leads to R̃TM/TE 1,2=-1 and consequently, the electric field is zero at θ=π/2 direction. So clearly the radiation pattern of a nanoantenna located above a planar layered material, is much different from that in the absence of layered-substrate. By utilizing the dipole approximation for the nano-particles, in the next section we discuss the required steps for designing a narrow beamwidth antenna.

3. Antenna applications

In practice, it is often required to design an antenna system that will yield the prescribed far-field radiation pattern. For example, a very common request is to design an antenna whose pattern has a narrow beamwidth and small side lobes. At radio frequencies, an approach to achieve this, is based on the concept of Yagi-Uda antennas [27]. The high directivity is achieved by placing a reflector and several directors arranged in proper locations around a resonant feed element. The directors are designed to be more capacitive while the reflector is more inductive, resulting in a narrow radiation beam toward the direction of directors and a minimum toward the reflector. To take advantage of this principle in optical regime, it has been suggested to use plasmonic nano-particles as the antenna element and to put an optical dipole source as the feed element [7]–[14]. In reality the antenna must be placed on a substrate where surface wave couplings-interactions with the plasmonic particles antenna elements can drastically change the radiation performance. This phenomenon is addressed carefully in our work by applying the developed formulations. A composite of indium gallium arsenide (InGaAs)-silver (Ag)-silicon dioxide (SiO2) substrate is engineered to first suppress the surface waves propagation and second tailor the radiation beam in a desired direction enabling an on-chip nanoantenna.

Figure 4(a) illustrates a Yagi-Uda type optical nanoantenna constructed from a reflector core-shell sphere located at r ref=-0.2λ 0y and a director located at r dir=0.3λ 0y. Both particles have the core made of SiO 2 (εrc=2.2+0.01i) and the shell made of silver (εrs=-8.06+0.175i). For the working frequency of f=650 THz, the outer radius of each particle is 0.1λ 0 whereas the radii ratios for reflector and director are r 1r/r 2r=0.75 and r 1d/r2d=0.65 respectively. The ratios of radii r 1/r 2 are chosen to set the nano spheres operate around their resonances. Further, under this design, the phase of the scattering coefficient for the reflector compared to the incident field is greater than π/2 (inductive), and that for the director is less than π/2 (capacitive). These conditions ensure successful Yagi-Uda antenna realization. Note that, unless otherwise mentioned, the optical dipole source throughout this work is perpendicular to the top surface of the multilayer substrate (p=z). The radiation patterns for the 2-element Yagi-Uda depicted in Fig. 4(a), are obtained using our theoretical model and are demonstrated in Figs. 4(b) and 4(c) in xy and yz planes, respectively. In this scenario, since the optical dipole is polarized in the z direction, the induced dipoles on the nano core-shells are also polarized in the z direction. The polarizabilities for the reflector and the director associated with the radii ratio are αr=0.059e i0.77π ε oλ30 and αd=0.052e i0.21π ε o λ 3 0 respectively. From the polarizabilities one can notice that the phases of the reflector and the director are shifted from π/2 indicating the inductive and capacitive behavior for the reflector and the director respectively. It is worth noticing that, all the radiation patterns plotted here are the radiated power normalized to the maximum value of the radiated power of the optical dipole source in the absence of the nano core-shells (and the layered substrate).

To validate the accuracy of our dipole-mode model, a full wave numerical analysis based on an in-house developed FDTD method [23] is applied for characterizing the plasmonic core-shell structures and determining the radiation patterns (as shown in Figs. 4(b) and (c) with the red dashed lines). The near-fields based on theory and numerical techniques are also compared in Figs. 5(a) and (b). Good comparisons are observed. It is worth highlighting that, FDTD characterizes the complete core-shell particles while in our theoretical model we approximate them with dipole elements located at the core-shells’ centers. Therefore, the near field characteristics of the antenna obtained by applying our theoretical approach show the field behavior of dipoles. This implies that, the electric field has local maxima at the location of the dipoles (centers of the core-shells). While the FDTD provides the dipole behavior outside the core-shells. Namely one does not observe the local maxima in the center of the core-shells. Nevertheless, considering the small size of particles and the fact that we are interested in the field performance outside the particles, this is a very good approximation as observed from the comparisons.

It must also be mentioned that due to the spherical complexity of the particles and our FDTD cubical meshing, a fine discretization with stair-casing around the boundary is used. The dispersive property of the plasmonic shells adds another difficulty for the FDTD analysis. The FDTD obtained resonant frequencies for the core-shell nanoantenna particles have almost 8% error in compared to the theoretical prediction. Using smaller cell-size can reduce the error to about 5%. Here, we scale properly the FDTD operating frequency to be around the resonance in accordance with our theory implementation for the particles. Notice that, unless otherwise mentioned, the near-field distributions are normalized to the maximum value of Ez in the yz plane.

 figure: Fig. 4.

Fig. 4. (a) A Yagi-Uda type antenna constructed from a dipole source exciting two nano core-shells, operating at frequency f=650 THz. The reflector’s radii ratio is r 1r/r 2r=0.75 and the radii ratio of the director is r 1d/r 2d=0.65. Both core-shells have core made of SiO 2(εrc=2.2+0.01i) and shell made of silver (εrs=-8.06+0.175i). The normalized radiated power obtained by theory and FDTD, (b) in xy plane, and (c) in yz plane. A good comparison is observed. Note that FDTD characterizes the actual plasmonic core-shell structure whereas in theoretical model we approximate them with dipole modes.

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. Magnitude of Ez(dB) for the nanoantenna in Fig. 4(a) in xy (z=0) plane and yz (x=0) plane. (a) Theory, (b) FDTD. It is worth highlighting again, that in our theory each nano core-shell is modeled with an electric dipole, where FDTD calculates the field for the actual structure. Outside the plasmonic core-shells one can expect the similar performance.

Download Full Size | PDF

The next step is to investigate the performance of the array antenna when it is located above an InGaAs substrate. This offers potential advantage for making on-chip nanoantennas. The 2-particle nanoantenna of the previous example is now located above an InGaAs slab of 0.35λ 0, see Fig. 6(a) (in our theoretical modeling a refractive index of about n ~3.4 is assumed). The radiation pattern and the electric fields in xy and yz planes are shown in Figs. 6(b) and 7. As highlighted earlier, due to the surface wave propagation along the substrate and its strong interactions with the plasmonic antennas the radiation beams are degraded and the radiated power is split in different directions. Similar observation has been reported in [11, 26]. Far field behaviors reveals the antenna beam has a maximum in the xz plane. Although the source is polarized along the z, but because of the plasmonic particles locations along the y and their slab interactions, there are some y-polarized induced dipoles for the plasmonic particles, cf. Fig. 7(a). This can be seen from (8) where Gelyrz is not zero, and hence the induced dipole have both y and z components. In addition, the near field in the xz plane in Fig. 7(b) shows that the electric field above and below the substrate are comparable verifying the field propagation through the substrate. The FDTD results for the near-fields are also illustrated in Fig. 8 validating our Green’s function base formulations and concluding discussions.

 figure: Fig. 6.

Fig. 6. (a) Nanoantenna configuration in Fig. 4(a) located above a slab of InGaAs with the thickness of 0.35λ 0. (b) The 3D radiation pattern of the radiated power. The substrate considerably degrades the antenna radiation pattern.

Download Full Size | PDF

For the sake of completion, we also investigate the behavior of a nanoantenna, when the source feed is in parallel direction with the slab surface, cf. Fig. 9(a). The radiation pattern and the near-field (normalized to the maximum value of Ex) in yz plane are demonstrated in Figs. 9 (b) and (c), respectively. In this scenario, the feed dipole is polarized along x (p=x) and hence the induced dipoles on the nano core-shells are also x-polarized. Fig. 9 (b) exhibits that the surface wave generation has a dominant effect on the radiation performance.

 figure: Fig. 7.

Fig. 7. Near field distribution for the antenna in Fig. 6(a) obtained by using the theoretical approach of this paper, (a) The magnitude of Ey(dB) in xy plane, (b) The magnitude of Ez(dB) in yz plane. Field penetration and propagation through the substrate is illustrated.

Download Full Size | PDF

The back radiation can be suppressed by coating the silicon surface with a layer of silver material with ε=-8.06+0.175i and thickness of 0.2λ 0. The two plasmonic radiators are now implanted at the top of the silver, refer to Fig. 10(a). Figures. 10(b) and 11 represent the far field radiation and the near field behavior of the antenna. As it is demonstrated, the back-radiation is stopped and the nanoantenna array radiates towards the up where the side lobe level gets smaller. Also the the near field distribution of the antenna proves that the magnitude of the field below the substrate is considerably smaller (around 50dB) than that in the forward region (upper half-space). That is to say, by adding a plasmonic layer we prevent the wave to transmit on the other side of the substrate.

 figure: Fig. 8.

Fig. 8. Near field distribution for the antenna in Fig. 6(a) obtained utilizing FDTD numerical technique, (a) The magnitude of Ey(dB) in xy plane, (b) The magnitude of Ez(dB) in yz plane.

Download Full Size | PDF

 figure: Fig. 9.

Fig. 9. (a) A nanoantenna array of a horizontal dipole, px, and two nano core-shells located above a InGaAs slab of 0.35λ0 thickness. (b) The 3D radiation pattern. (c) The magnitude of Ex(dB) in yz plane.

Download Full Size | PDF

From the above, and considering the strong effect of the layered materials on the plasmonic nanoantenna performance, one may engineer the substrate of an array antenna properly to achieve a desired radiation pattern characteristic. Basically, the substrate and the reflected wave from that can control the dipolar polarization of each plasmonic-particle array elements, shaping the radiation beam. To provide some physical understanding, it is valuable to study the performance of a single optical dipole located at the top surface of a planar layered material, refer to Fig. 12(a). Fig. 12(b) shows the angle of the maximum beam (in the yz plane) for an optical dipole located above an InGaAs slab (of 0.35λ 0) which its top surface is coated with 0.2λ 0 of silver. The dashed line exhibits the angle of maximum radiation for the antenna, versus the distance from the layered material (ε r2=1). As observed, with changing the distance between the dipole and the multi layer substrate the angle of radiation changes from 30± to around 75°. Note that, as we increase the distance between the dipole and the layered material the number of side lobes increases (see also [26]) so the beam of the main lobe changes and a jump in the maximum radiation is exhibited. Adding another layer to this configuration will provide better opportunity for controlling the beam radiation angle. We consider a case where a third layer of SiO 2 medium is grown over the silver and the dielectric slab (InGaAs). The solid line in Fig. 12(b) demonstrates that the beam angle (for maximum radiation) varies between 30° to around 90°. Evidently, by depositing the third dielectric layer, we have a better control on the angle of radiation.

 figure: Fig. 10.

Fig. 10. (a) Configuration of the nanoantenna in Fig. 4(a) where the nano core-shells are located above a InGaAs slab of 0.35λ 0 which its top surface is coated with 0.2λ 0 of silver. (b) The 3D radiation pattern. The silver layer suppresses the back radiation.

Download Full Size | PDF

 figure: Fig. 11.

Fig. 11. Near field distribution for the antenna in Fig. 10(a). (a) The magnitude of Ey(dB) in xy plane, (b) The magnitude of Ez(dB) in yz plane.

Download Full Size | PDF

A three-layers substrate can be engineered to manipulate the plasmonic particles dipolar modes and tilt the beam of a two-particles free-space Yagi-Uda antenna from endfire-type radiation to close broadside performance. The layered-substrate includes an InGaAs slab with thickness of 0.35λ 0, silver with thickness of 0.2λ 0 and a coating dielectric with εr=2.2+0.01i and thickness of 0.1λ 0 (refer to Fig. 13(a)). The radiation performance is obtained in Fig. 13(b). It shows the main lobe around 16° with suppressed back radian (due to negative permittivity of the silver layer). This is an interesting result, because we started with a z-directed optical dipole source which has a maximum radiation around θ=90° and we end up with a maximum radiation around 16°, by tailoring the plasmonic particles array and the substrate. The near-field performance is plotted in Fig. 14(a) and (b).

As observed from Fig. 13, the two nano core-shells parasitic radiators array and its engineered substrate provides a narrow beam only in the yz plane. To control the 3D radiation performance in both planes (xz and yz) and achieve a pencil-beam broadside radiation one may change the distances of the nano core-shells and/or place four plasmonic nano core-shells around the dipole excitation. Figure 15(a) shows the configuration for an optimized design where four plasmonic nano core-shells are placed such that the antenna has a maximum radiation around θ=0. In this case by using the data from Fig. 2(a), the radii ratios for the nano particles are chosen to be closer to the 0.71 (which is the radii ratio of the resonance for f=650 THz). Hence the polarizability of the particles are larger and the antenna can have higher power transmission. Figure 15(b) demonstrates the 3D radiation pattern of the radiated power. As expected a pencil-beam around the z-axis is successfully achieved. The four plasmonic nanoantennas manipulate the induced dipole polarizations, and hence control the radiation characteristic to the application of interest.

 figure: Fig. 12.

Fig. 12. (a) A vertical dipole located at the top of a planar layered material. The layered substrate includes a slab of InGaAs with thickness of 0.35λ 0 where its top surface is coated with silver of thickness 0.2λ 0. A third layer of dielectric ε r2 is deposited at top of the silver. (b) The angle of maximum radiation in the yz plane vs. the thickness of the dielectric layer. The jump in angle is because the maximum radiation occurs in another direction.

Download Full Size | PDF

 figure: Fig. 13.

Fig. 13. (a) Nanoantenna configuration in Fig. 4(a) located above a 3 layered substrate including InGaAs slab of 0.35λ 0, silver with thickness of 0.2λ 0 and a 0.1λ 0 layer of SiO 2. (b) The 3D radiation pattern of the radiated power. Adding a third layer controls better the radiation pattern.

Download Full Size | PDF

From practical points of view, it is challenging to precisely fabricate the core-shell antenna. Specially one might be concerned about the exact thickness of the plasmonic shell that goes beyond one nanometer. Therefore, it is important to investigate the effect of changing the radii ratio. Figure 16 shows the change in the far field behavior, if the thickness of the shell varies. As it is evident the direction of the maximum beam is not so sensitive to the plasmonic shell thickness. Nevertheless the power transmission of the antenna changes.

We need to highlight again that the designed structures presented in this paper are achieved for infinite-size substrates (along the x and y directions), where a finite substrate size antenna can have different performance due to the type of surface wave generation, the fields propagation in and along the substrate, and the edges diffractions. It will be instructive to investigate the performance of an array of four plasmonic core-shells located over a finite-size substrate with dimensions of 1.3λ 0×1.3λ 0, as depicted in Fig. 17(a). The FDTD results for the radiation patterns and near-fields are shown in Fig. 17(b) and (c). The finite size substrate has less effect on the radiation pattern of the 2D Yagi-Uda antenna compared to its counterpart studies in Fig.15 (infinite layers). This is because of the small size of the substrate where a surface wave cannot be supported much and instead the waves penetrate through the layers from the edges. Since the infinite case is the limit of the finite one, when the size of the structure approaches to infinity the field at the back can be rather small leading to a strong directional beam above, as the one we obtain from the theoretical computation. Although, in this case the edges diffractions can have some undesirable effects on the radiation pattern increasing the side-lobes. The optimum design for a plasmonic array nanoantenna located above a finite substrate is currently under investigation.

 figure: Fig. 14.

Fig. 14. Near field distribution for the antenna in Fig. 13(a), (a) The magnitude of Ey(dB) in xy plane, (b) The magnitude of Ez(dB) in yz plane.

Download Full Size | PDF

 figure: Fig. 15.

Fig. 15. (a) 4-particle nanoantenna array with an engineered substrate demonstrating a broadside radiation characteristic. The radii ratio for the particles are r 1 r/r 2r=0.72 and r 1d/r 2d=0.69, and dref=0.25λ 0 while ddir=0.65λ 0. The layered material includes 0.5λ 0 of InGaAs, 0.2λ 0 of silver coating and a dielectric layer (SiO 2) with the thickness of 0.1λ 0. (b) The 3D radiation pattern.

Download Full Size | PDF

4. Conclusion

In this paper, the radiation performance of an array of plasmonic core-shell nanoantennas located above engineered substrates are explored. We develop a theoretical-modeling framework base on the Green’s function analysis of dipole excitation above layered structures, to successfully characterize the antenna performance. The plasmonic core-shells are approximated with electric dipoles where their induced dipoles are calculated through the Greens function modeling of dipoles over the layered materials. The theoretical approach and developed formulations are verified through applying a full wave numerical analysis utilizing FDTD method. Both the radiation patterns and near-field characteristics are determined.

The effect of the dielectric substrate layers on the radiation performance is highlighted. It is obtained that a dielectric substrate can drastically change the radiation pattern of a plasmonic nanoantenna from what is obtained for that when is located in free space. Integrating a silver layer can reduce the antenna interaction with the dielectric substrate and suppress the back radiation. A composite substrate for an optimized 2D array of four-plasmonic nano core-shell radiators is engineered to tailor the dipolar modes of the nanoantennas and accomplish a pencil beam radiation.

 figure: Fig. 16.

Fig. 16. Radiation pattern sensitivity to the change of r 1r for the configuration shown in Fig. 13(a). The radiation performance in yz plane, where r 2r=r 2d=46 nm and r 1d=30 nm.

Download Full Size | PDF

 figure: Fig. 17.

Fig. 17. (a) 4-particle nanoantenna array located above a finite-size substrate. In this case the radii ratios are r 1r/r 2r=0.75 and r 1d/r 2d=0.65. (a) The radiation pattern in xy and yz planes, (b) The near-field distribution of Ez(dB) in yz plane. The FDTD technique demonstrates the effect of finite-size substrate on the field penetration through the substrate and the edge diffractions (as this will establish different radiation characteristic compared to what is obtained for the infinite-size substrate model).

Download Full Size | PDF

The FDTD is applied to manipulate the performance of an array of plasmonic nanoantennas located above a finite-size substrate. A small size substrate may not have much effect on the radiation performance, where increasing the size can tailor the plasmonic particles polarizations and control the radiation pattern of interest.

The results of this paper can be of great value for creating on-chip nanoantennas with desired high-performance radiation characteristics.

Acknowledgement

This work was supported in part by the U.S. Air Force Office of Scientific Research (AFOSR) and National Science Foundation (NSF).

References and links

1. E. Ozbay, Plasmonics: Merging photonics and electronics at nanoscale dimensions, Science, 311, 189–193, (2006). [CrossRef]   [PubMed]  

2. S. I. Bozhevolnyi and V. M. Shalaev, Nanophotonics with surface plasmonsPart I, Science, 40, 58–66, (2006).

3. S. I. Bozhevolnyi and V. M. Shalaev, Nanophotonics with surface plasmonsPart II, Science, 40, 66–72, (2006).

4. S. A. Maier and H. A. Atwater, Plasmonics: Localization and guiding of electromagnetic energy in metal/dielectric structures, J. Appl. Phys., 98, 011101, (2005). [CrossRef]  

5. L. A. Sweatlock, S. A. Maier, H. A. Atwater, J. J. Penninkhof, and A. Polman, Highly confined electromagnetic fields in arrays of strongly coupled ag nanoparticles, Phys. Rev. B., 71, 235408, (2005). [CrossRef]  

6. V. Podolskiy, A. Sarychev, and V. Shalaev, Plasmon modes in metal nanowires, J. Nonlinear Opt. Phys. Mater., 11, 65–74, (2002). [CrossRef]  

7. J. Li, A. Salandrino, and N. Engheta, Shaping light beams in the nanometer scale: A yagi-uda nanoantenna in the optical domain, Phys. Rev. B., 76, 245403, (2007). [CrossRef]  

8. J. Li and N. Engheta, Optical leaky-wave nano-antennas using plasmonic nanowires with periodical variation of permittivity, (presented at the 2005 Annual Meeting of the OSA., Tucson, Arizona, 1620, 2005).

9. J. Li, A. Salandrino, and N. Engheta, Radiation characteristics and beam forming of multi-particles nanoantennas at optical frequencies, (iWAT 06: Small Antennas and Novel Metamaterials, White Plains, NY, 432–433, 2006).

10. P. Muhlschlegel, H. J. Eisler, O. J. F. Martin, B. Hecht, and D.W. Pohl, Resonant optical antennas, Science, 308, 16071609, (2005). [CrossRef]  

11. T. H. Taminiau, F. D. Stefani, and N. F. van Hulst, Enhanced directional excitation and emission of single emitters by a nano-optical yagi-uda antenna, Optics express, 16, 10858–10866, (2008). [CrossRef]   [PubMed]  

12. D. R. Jackson, T. Zhao, J. T. Williams, and A. A. Oliner, Leaky surface-plasmon theory for dramatically enhanced transmission through a sub-wavelength aperture, Part II: Leaky-wave antenna model, (IEEE International Symposium on Antenna and Propagations, 2, 1095–1098, 2003).

13. N. C. Panoiu and R. M. Osgood, Optical antenna arrays in the visible range, Opt. Lett., 32, 2825, (2007). [CrossRef]   [PubMed]  

14. H. F. Hofmann, Terukazu, and Y. Kadoya, Design parameters for a nano-optical yagi-uda antenna, New J. Phys., 9, 217, (2007). [CrossRef]  

15. K. Nakayama, K. Tanabe, and H. A. Atwater, Plasmonic nanoparticle enhanced light absorption in gaas solar cells, Applied Physics Letters, 93, 121904, (2008). [CrossRef]  

16. D. R. Matthews, H. D. Summers, K. Njoh, S. Chappell, R. Errington, and P. Smith, Optical antenna arrays in the visible range, Opt. Express, 15, 3479–3487, (2007). [CrossRef]  

17. D. M. Schaadt, B. Feng, and T. F. Yu, Optical antenna arrays in the visible range, Appl. Phys. Lett., 86, 063106, (2005). [CrossRef]  

18. V. E. Ferry, L. A. Sweatlock, D. Pacifici, and H. A. Atwater, Plasmonic nanostructure design for efficient light coupling into solar cells, Nane letters, 8, 12, (2008).

19. S. Pillai, K. R. Catchpole, T. Trupke, and M. A. Green, Surface plasmon enhanced silicon solar cells, Journal of Applied Physics, 101, 093105, (2007). [CrossRef]  

20. W. C. Chew, Waves and Fields in Inhomogeneous Media, IEEE Press, New York, (1995).

21. S. Ghadarghadr and H. Mosallaei, Array of plasmonic nanoparticles enabling energy coupling-guiding in solar systems: A theoretical analysis, (Optics and Photonics for Advanced Energy Technology, Cambridge, MA, 2009).

22. S. Ghadarghadr and H. Mosallaei, Nanoantennas Array Enabling Optical Communication, (Frontiers in Optics 2009/Laser Science XXV, San Jose, CA, 2009).

23. H. Mosallaei, FDTD-PLRC technique for modeling of anisotropic-dispersive media and metamaterial devices, IEEE Trans. Electromagn. Compat., 49, 649–660, (2007). [CrossRef]  

24. J. A. Stratton, Electromagnetic Theory, McGraw Hill, New York, (1941).

25. P. B. Johnson and R. W. Christy., Optical constants of the nobel metals, Phys. Rev. B, 6, 43704379, (1972).

26. L. Novotny and B. Hecht, Principles of Nano-Optics, Cambrdige University Press, United Kingdom, (2006).

27. C. A. Balanis, Antenna Theory: Analysis and Design, John Wiley & Sons, (2005).

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (17)

Fig. 1.
Fig. 1. Array of plasmonic nanoantennas located above a layered substrate.
Fig. 2.
Fig. 2. Magnitude and phase of normalized polarizability of a plasmonic core-shell sphere versus r 1/r 2 ratio for different frequencies for a structure with r 2=0.1λ, εrc =2.2+0.01i and silver coating. (a) Polarization magnitude, and (b) Polarization phase, (c) Silver permittivity behavior.
Fig. 3.
Fig. 3. Maximum percentage of error in obtaining the integrals in (7) and (8) using our theoretical model compared to their exact values for different observation points. ρ and z0 denote the location of the observation point.
Fig. 4.
Fig. 4. (a) A Yagi-Uda type antenna constructed from a dipole source exciting two nano core-shells, operating at frequency f=650 THz. The reflector’s radii ratio is r 1r /r 2r =0.75 and the radii ratio of the director is r 1d /r 2d =0.65. Both core-shells have core made of SiO 2(εrc =2.2+0.01i) and shell made of silver (εrs =-8.06+0.175i). The normalized radiated power obtained by theory and FDTD, (b) in xy plane, and (c) in yz plane. A good comparison is observed. Note that FDTD characterizes the actual plasmonic core-shell structure whereas in theoretical model we approximate them with dipole modes.
Fig. 5.
Fig. 5. Magnitude of Ez (dB) for the nanoantenna in Fig. 4(a) in xy (z=0) plane and yz (x=0) plane. (a) Theory, (b) FDTD. It is worth highlighting again, that in our theory each nano core-shell is modeled with an electric dipole, where FDTD calculates the field for the actual structure. Outside the plasmonic core-shells one can expect the similar performance.
Fig. 6.
Fig. 6. (a) Nanoantenna configuration in Fig. 4(a) located above a slab of InGaAs with the thickness of 0.35λ 0. (b) The 3D radiation pattern of the radiated power. The substrate considerably degrades the antenna radiation pattern.
Fig. 7.
Fig. 7. Near field distribution for the antenna in Fig. 6(a) obtained by using the theoretical approach of this paper, (a) The magnitude of Ey (dB) in xy plane, (b) The magnitude of Ez (dB) in yz plane. Field penetration and propagation through the substrate is illustrated.
Fig. 8.
Fig. 8. Near field distribution for the antenna in Fig. 6(a) obtained utilizing FDTD numerical technique, (a) The magnitude of Ey (dB) in xy plane, (b) The magnitude of Ez (dB) in yz plane.
Fig. 9.
Fig. 9. (a) A nanoantenna array of a horizontal dipole, px, and two nano core-shells located above a InGaAs slab of 0.35λ0 thickness. (b) The 3D radiation pattern. (c) The magnitude of Ex (dB) in yz plane.
Fig. 10.
Fig. 10. (a) Configuration of the nanoantenna in Fig. 4(a) where the nano core-shells are located above a InGaAs slab of 0.35λ 0 which its top surface is coated with 0.2λ 0 of silver. (b) The 3D radiation pattern. The silver layer suppresses the back radiation.
Fig. 11.
Fig. 11. Near field distribution for the antenna in Fig. 10(a). (a) The magnitude of Ey (dB) in xy plane, (b) The magnitude of Ez (dB) in yz plane.
Fig. 12.
Fig. 12. (a) A vertical dipole located at the top of a planar layered material. The layered substrate includes a slab of InGaAs with thickness of 0.35λ 0 where its top surface is coated with silver of thickness 0.2λ 0. A third layer of dielectric ε r2 is deposited at top of the silver. (b) The angle of maximum radiation in the yz plane vs. the thickness of the dielectric layer. The jump in angle is because the maximum radiation occurs in another direction.
Fig. 13.
Fig. 13. (a) Nanoantenna configuration in Fig. 4(a) located above a 3 layered substrate including InGaAs slab of 0.35λ 0, silver with thickness of 0.2λ 0 and a 0.1λ 0 layer of SiO 2. (b) The 3D radiation pattern of the radiated power. Adding a third layer controls better the radiation pattern.
Fig. 14.
Fig. 14. Near field distribution for the antenna in Fig. 13(a), (a) The magnitude of Ey (dB) in xy plane, (b) The magnitude of Ez (dB) in yz plane.
Fig. 15.
Fig. 15. (a) 4-particle nanoantenna array with an engineered substrate demonstrating a broadside radiation characteristic. The radii ratio for the particles are r 1 r/r 2r =0.72 and r 1d /r 2d =0.69, and dref=0.25λ 0 while ddir =0.65λ 0. The layered material includes 0.5λ 0 of InGaAs, 0.2λ 0 of silver coating and a dielectric layer (SiO 2) with the thickness of 0.1λ 0. (b) The 3D radiation pattern.
Fig. 16.
Fig. 16. Radiation pattern sensitivity to the change of r 1r for the configuration shown in Fig. 13(a). The radiation performance in yz plane, where r 2r =r 2d =46 nm and r 1d =30 nm.
Fig. 17.
Fig. 17. (a) 4-particle nanoantenna array located above a finite-size substrate. In this case the radii ratios are r 1r /r 2r =0.75 and r 1d /r 2d =0.65. (a) The radiation pattern in xy and yz planes, (b) The near-field distribution of Ez (dB) in yz plane. The FDTD technique demonstrates the effect of finite-size substrate on the field penetration through the substrate and the edge diffractions (as this will establish different radiation characteristic compared to what is obtained for the infinite-size substrate model).

Equations (25)

Equations on this page are rendered with MathJax. Learn more.

p = α . Etotal = i 6πε0k03 Γ1(e) . Etotal ,
Γn(e) = Un(e)Un(e)+iVn(e) ,
Un(e) = jn(kcr1)jn(ksr1)yn(ksr1)0jn(k˜cr1)εcjn(k˜sr1)εsyn(k˜sr1)εs00jn(ksr2)yn(ksr2)jn(k0r2)0j˜n(ksr2)y˜n(ksrs)εsj˜n(k0r2)ε0,
Vn(e) = jn(kcr1)jn(ksr1)yn(ksr1)0jn(k˜cr1)εcjn(k˜sr1)εsyn(k˜sr1)εp00jn(ksr2)yn(ksr2)yn(k0r2)0j˜n(ksr2)y˜n(ksrs)εsy˜n(k0r2)ε0.
εrs = 1 1.513×1032ω(ω+i7.9×1013) .
pl = αl (Einctotal(rl)+q,qlG̿dipolel(rl,rq)pq+qG̿reflectedl(rl,rq)pq) ,
G̿dipolel (rl,rq)=[(k12+22x)2xy2xz2xy(k12+22y)2yz2xz2yz(k12+22z)] eik1rlrq4πε1rlrq ,
G̿reflectedl (rlrq)=[Gerxlx(rlrq)Gerylx(rlrq)Gerzlx(rlrq)Gerxly(rlrq)Geryly(rlrq)Gerzly(rlrq)Gerxlz(rlrq)Gerylz(rlrq)Gerzlz(rlrq)].
G erxlz = 18πε1 cos ϕ dkρkρ2H1(1)(kρρ)R˜1,2TMeik1z(z+2d1),
G erylz = 18πε1 sin ϕ dkρkρ2H1(1)(kρρ)R˜1,2TMeik1z(z+2d1),
G erzlz = i8πε1 dkρkρ3k1zH0(1)(kρρ)R˜1,2TMeik1z(z+2d1),
G˜ erwls (kρ)=1kρ2[zG˜erwls(kρ)zμẑ×sG˜hrwlz(kρ)],
R˜i,i+1TMTE = Ri,i+1TMTE+R˜i+1,i+2TMTEe2iki+1,z(di+1di)1+Ri,i+1TMTER˜i+1,i+2TMTEe2iki+1,z(di+1di) ,
Ri,jTMTE = (εjμj)kiz(εiμi)kjz(εjμj)kiz+(εiμi)kjz ,
kiz = ki2kρ2 .
p = (I̿αG̿reflected)1 α Etotal ,
E = EθEϕ=kJ24πεJeikJrr(pxcosϕ+pysinϕ)cosθΦJ2pzsinθΦJ1(pxsinϕpysinϕ)ΦJ3.
Φ11 = eik1r [1+R˜1,2TM(θ)eik1z0],
Φ12 = eik1r [1R˜1,2TM(θ)eik1z0],
Φ13 = eik1r [1+R˜1,2TM(θ)eik1z0],
ΦN1 = ΦN2 = nNn1 cosθs˜z TTM (θ)eikN[z0s˜z+(dNd1)cosθ],
ΦN3 = cosθs˜z TTM (θ)eikN[z0s˜z+(dNd1)cosθ],
kJ r = kJ [x0cosϕsinθ+y0sinϕsinθ+z0cosθ] ,
kJ z = 2 kJ (z0+d1)cosθ,
s˜z = k1zkN = (n1nN)2sin2θ .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.