Analysis of Nanodipoles in Optical Nanocircuits Fed by Gaussian Beam

We analyze impedance matching and excitation properties of a plasmonic optical nanocircuit composed by a receiving and an emitting dipole connected by a two-wire optical transmission line. The circuit is fed by a circular Gaussian beam focused on the receiving dipole. The numerical analysis is performed by linear method of moments with a given surface impedance of gold cylindrical conductors. With this model, we analyze the variation of standing-wave response along the circuit in function of some geometrical parameters. We present some conclusions concerning impedance matching between the transmission line and the emitting dipole and coupling between the receiving dipole and the incident Gaussian beam.


Introduction
Optical antennas are metal nanostructures used to transmit or receive optical fields [1][2][3][4].This definition is similar to that of conventional antennas in RF-microwave regimes.However, optical antennas can also be used to enhance and confine optical fields in nanometer regions smaller than the operating wavelength, that is, beyond the diffraction limit of light.This is possible due to the plasmonic effect of metals in optical frequencies which produces strong electrons oscillations named plasmonic resonances [5,6].
The concept of optical antennas was firstly used in design of efficient near-field optical probes.Classical papers about this application are [7], where a microwave version of optical probe based on nanoantenna was experimentally analyzed, and [8], where several configurations and new ideas of optical probes based on optical antennas were theoretically presented.More recent experimental and theoretical works on optical probes based on nanoantennas can be found in [9,10].Other examples of applications are enhancement of single molecule fluorescence [11][12][13][14], ultrahigh-density optical data storage devices [15], chemical enhancement for surface-enhanced Raman scattering [16], nanobioimaging [17], enhancement in the directional excitation and emission of single emitters [18,19], plasmonic laser antennas [20], and plasmonic photovoltaic cells [21,22].Potential applications of optical antennas and plasmonic waveguides are in nanophotonics.These elements can be used, for example, as fundamental components for the design of highly integrated photonic signal-processing systems, because they can manipulate optical fields in subwavelength nanoscales, beyond the diffraction limit [23].In this case, plasmonic waveguides can be used to interconnect different devices in an optical nanocircuit, and optical antennas can act as terminal elements that transform far field radiation into guided waves and vice versa.Some examples of plasmonic waveguide are two-wire optical transmission line (OTL) [24] and array of nanoparticles, acting as lumped nanocircuit elements [25].
Optical nanocircuits based on two-wire OTL and nanoantennas have been investigated in [26,27].In [26], a theoretical analysis of an optical nanocircuit composed by a receiving and an emitting nanoantenna connected by a two-wire OTL is presented.In this paper, the authors performed the input impedance matching varying the length of the nanodipoles for a fixed frequency, radius, and gap of the OTL.In [27], analysis of an optical nanocircuit fed by an aperture probe is fulfilled.The circuit is composed by a receiving and an emitting nanoantenna connected by a two-wire OTL and the coupling between the aperture probe and the receiving dipole is modeled by an equivalent voltage source.

International Journal of Antennas and Propagation
In this paper, we present an impedance matching and excitation analysis of a plasmonic optical nanocircuit composed by a receiving and an emitting nanodipoles connected by a two-wire OTL.The feeding of this circuit is made by an aperture probe that illuminates the receiving dipole.This circuit is similar to that presented in [27], but in the present case we model the radiated field from the aperture probe by a focused Gaussian beam [28], with  = 830 nm, beam waist of 340 nm, and electric field with polarization matched to the axis of the receiving dipole.The numerical analysis of this circuit is made by a simple and efficient computational method based on method of moments (MoM) [27].We study the variation of the induced current, voltage reflection coefficient and near-field distribution for different length and radius of the receiving and emitting dipole, and different radius and gap of the OTL.We present some conclusions about the conditions where one can obtain a better impedance matching between the emitting dipole and the OTL.We also discuss the conditions which allow one to obtain an efficient electromagnetic coupling between the Gaussian beam and the receiving dipole.The next sections present the theoretical development, impedance matching analysis, excitation analysis, and conclusions.

Theoretical Development
This section presents the geometry of the problem, the linear method of moments (MoM) model of the circuit, the Gaussian beam used to feed the circuit, and a numerical example of simulation.

Description of the Problem.
The analyzed optical circuit is shown in Figure 1 where an aperture probe is electromagnetically coupled to a nanocircuit.The nanocircuit is composed of a receiving dipole, two-wire optical transmission line, and an emitting dipole.The material of the circuit is gold and the medium is vacuum with no substrate.We model the coupling between the aperture probe and the receiving antenna by an equivalent focused Gaussian beam radiated from the aperture.The polarization of this field is parallel to the axis of the receiving dipole (axis ) to maximize the electromagnetic coupling and efficiency.
We calculate the impedance matching between the OTL and the emitting dipole and optimize the electromagnetic coupling between the receiving dipole and the Gaussian beam.We solve numerically the scattered problem by the method of moments and calculate the induced current, voltage reflection coefficient, and near-field distribution.

Method of Moments Model.
To analyze the equivalent nanocircuit of the system presented in Figure 1, we use the linear method of moments (MoM) approximation with sinusoidal basis functions and equivalent surface impedance [27].Figure 2 shows the geometry of the original problem and the MoM equivalent model.In Figure 2(a), we depicted receiving dipole 1 (left), an OTL, and an emitting dipole 2 (right).
The dimensions of this circuit are ℎ 1 ,  ℎ1 and ℎ 2 ,  ℎ2 the arm length and radius of the nanodipoles 1 and 2, respectively,  and   the length and radius of the OTL wires, respectively,  the distance between the axes of the OTL, and  =  − 2  the gap distance between the surfaces of the OTL.The total length of nanodipoles 1 and 2 is  1 = 2ℎ 1 + and  2 = 2ℎ 2 +, respectively.
The losses in metal are described by surface impedance   .This surface impedance can be obtained approximately by considering cylindrical waveguide with the mode TM 01 [29]: where  =  0 √ 1 and  0 =  √  0  0 .The boundary condition for the electric field satisfied at the surface's conductor is (  +   ) ⋅   =   , where   is a unitary vector tangential to the surface of the metal,   is the scattered electric field due to the induced linear current  on the conductor,   is the incident electric field of the Gaussian beam source (Figure 1), and  is the induced longitudinal current in a given point of the nanodipole.The integral equation of the scattered field along the length  of the nanodipole is given by [27] where () =   conventional testing procedure.The following linear system of equations is obtained: where   is the mutual impedance between sinusoidal current elements  and , Δ  = [Δ  + Δ +1 ]/2.The solution of (4) produces the current along the nanocircuit.With these results, it is possible to calculate the near and far field distributions and other parameters.

Gaussian Beam Source.
The electromagnetic wave radiated from the aperture probe of Figure 1 illuminates the receiving dipole.It can be approximated by a Gaussian beam.This kind of wave is obtained by solving the scalar Helmholtz wave equation with the paraxial approximation [28].The magnetic vector potential of a Gaussian beam polarized in the axis  and traveling in the + direction is given by where  0 =(2P 0 /k) 1/2 , P is the power of the beam,  is the beam radius (Figure 3),  is the curvature radius of the phase front, and  is the phase of the beam.The three principal parameters that define the beam are the power , the radius of the beam waist  0 (Figure 3), and the operating wavelength .
With the vector potential (5), the  component of electric field   of the Gaussian beam can be obtained from the Maxwell equations [28].The beam of Figure 1 is focused on the receiving dipole 1 with polarization along the dipole axis (axis ), the direction of propagation is +, the beam axis is along the axis , and the minimum waist ( 0 ) is localized at  = 0, which is the plane of the nanocircuit.In all the analysis presented in this work, we consider a fixed Gaussian beam with power  = 10 −12 W, wavelength  = 830 nm, and beam waist  0 = Distance along the circuit (nm) 340 nm.We use these values for comparisons with the results of paper [26], but the analysis presented here can be applied for beams with other parameters.The field distribution of this beam is shown in Figure 3, where the electric field amplitude abs(  ) is presented at the planes  and  and the phase distribution angle of   is presented at the plane .The phase at plane  is constant.

Numerical Example.
Based on the theory presented above, we developed a MoM code in Matlab to analyze the nanocircuit shown in Figure 2. In this section, we present an example of simulation of the nanocircuit shown in Figure 2(a) fed by the Gaussian beam depicted in Figure 3. Figure 4(a) shows the geometry and discretization parameters used in this simulation and Figure 4(b) presents the result of the current distribution along the circuit.The near-field distribution for this example is given in Figure 5.We observe in these results the stationary behavior in the OTL, which is due to the mismatching in the impedances of nanodipole 2 and OTL.To make a quantitative measure of the impedance matching degree, we calculate approximately the voltage stationary wave ratio (VSWR) near nanodipole 2 as VSWR =  max / min , where  max and  min are, respectively, the maximum and minimum current magnitude nearest to dipole 2. With this parameter, we calculate the voltage reflection coefficient as |Γ  | = (VSWR − 1)/(VSWR + 1).In this numerical example we obtained |Γ  | = 0.4.
The next sections present two different analyses of the nanocircuit.In the first analysis, we fix the dimensions of the receiving dipole and focus on the impedance matching problem between the optical transmission line and the emitting dipole; that is, we are interested in the search for the parameters of the nanocircuit which provide maximum energy transferred to the emitting dipole and minimum reflected energy.In the second analysis, we fix the dimensions of the OTL and focus on the coupling problem between the Gaussian beam and the receiving antenna; that is, we are interested in the search for the better dimensions of the receiving dipole which provide maximum receiving energy.

Impedance Matching Analysis
This section presents a parametric analysis of the impedance matching of the nanocircuit for different values of ℎ 2 ,  ℎ2 ,   , and .In this analysis, we fixed the dimensions of the receiving dipole ℎ 1 = 173 nm,  ℎ1 = 10 nm and the length of the OTL  = 1200 nm. Figure 6  Analyzing these curves we come to some conclusions.We note that the nanocircuits possess in general smaller degree of input impedance matching (higher |Γ  |) when the gap of the OTL  is increased.The exceptions are the cases of Figures 6(a) and 6(b) in the range  2 < 400 nm, where we have a better matching for higher .
We also observe that in general the impedance matching is better when  ℎ2 >   .This means that the values of |Γ  | in Figure 6(d depends mainly on the attenuation of the current along the OTL, that is, depends on the loss constant  of the OTL.This parameter is constant for the principal mode that propagates on the OTL and can be obtained approximately by the average inclination of the current versus distance along the OTL.In equation form we have  = Δ/Δ, where Δ is the variation of the average amplitude of the current in decibels (dB) along a given distance Δ (nm) in the OTL.With this definition the unit of this parameter is dB/nm.In the numerical example presented in Figure 4(b) we have  ≈ 0.007393 dB/nm.This result is very close to that of  = 0.007296 dB/nm obtained in [26] where the OTL is similar to our case shown in Figure 4(a) (  = 15 nm and  = 10 nm).
To better understand the behavior of the impedance matching and efficiency characteristic of the results presented in Figure 6, we plot in Figures 7 and 8 the current distributions for different geometric parameters.Figure 7 shows the currents for two cases with the same   = 10 nm but different voltage reflection coefficient of |Γ  | = 0.31 and 0.67.In these results we observe a higher stationary wave for the case |Γ  | = 0.67 than for the case |Γ  = 0.31.This shows that our approximate method to calculate |Γ  | provides a good measure of the degree of impedance matching.We also observe that the two cases present the same attenuation along the circuit (i.e., in both cases one has approximately the same loss constant  = 0.0111 dB/nm) because the OTL are constructed with the same radius of wires   = 10 nm.
Figure 8 presents the current distribution for two cases with good impedance matching |Γ  | = 0.26 and 0.36 but with different loss constant of  = 0.0119 dB/nm and 0.0084 dB/nm, respectively.This difference is mainly due to the difference of the radius   of the OTL.For lower values of   the attenuation is higher in OTL and this result is similar to that observed in RF-microwave regimes.This can be explained by the surface impedance model of ( 2), where smaller radius produces higher   and, consequently, higher loss in the conductors.

Analysis of Receiving Dipole.
In this section, we investigate the response of the isolated receiving dipole illuminated by a Gaussian beam polarized along nanodipole's axis (left side of Figure 3).We search the values of  1 and  ℎ1 which maximize the coupling between the beam and the antenna, that is, maximize the resonant response of the dipole.Figure 9 shows the variation of the current amplitude in the middle of the nanodipole abs(  ) in function of the total length  1 and radius  ℎ1 .For comparison, this figure shows also the results calculated with the software Comsol [30].These International Journal of Antennas and Propagation currents are approximately the same input currents when a given nanodipole is connected to the input terminal of the nanocircuit in Figure 1.Thus when this current has the maximum value the power entering the OTL is also maximal.
The results of Figure 9 show that the maximum induced current   occurs for the resonant lengths  1 corresponding to the resonances /2 and 3/2.The values of these resonant lengths are given in Table 1.One example of the characteristic field distributions in these resonances is presented in Figure 10 where  1 = 111.2nm and  ℎ1 = 10 nm.Also, the results of Figure 9 show that for higher values of  ℎ1 we have larger   and the resonant lengths  1 are increased.We observed that the resonance /2 produces higher current in the nanodipole than in the resonance 3/2.We can conclude that we have a better coupling when the nanodipole operates in the first resonance and possesses a large radius.The latter condition implies high effective area of the nanodipole [31].

Voltage Reflection Coefficient.
In this section we analyze the voltage reflection coefficient of the nanocircuit in Figure 1 where we use the resonant receiving nanodipoles with dimensions given in Table 1. Figure 11 presents the variation of |Γ  | in function of  2 for the first resonance /2 and the second resonance 3/2, with different values of  ℎ1 .In these simulations we fixed the following parameters of the OTL and  emitting dipole:   = 15 nm,  ℎ2 = 20 nm,  = 1200 nm, and   = 40.We choose these parameters because they presented a better impedance matching and lower attenuation.The results of Figure 11 show that the reflection coefficient does not vary significantly with the radius  ℎ1 and length  1 for the first resonance (Figure 11(a)).However, for the second resonance (Figure 11(b)) we have a better input impedance matching for larger values of  1 and  ℎ1 .The best impedance matching is obtained for the case  ℎ1 = 20nm and  1 = 637.7 nm, where we have |Γ  | = 0.21.We also note that the results of impedance matching for the second resonance (Figure 11(b)) are better than those of the first resonance (Figure 11(a)).These results prove that the degree of impedance matching along the circuit not only is a function of the OTL and the emitting dipole but also depends on the geometry of the receiving dipole.

Current and Field Distribution.
Figure 12 shows the current distribution along the circuit for the best cases of |Γ  | presented in Figure 11.In Figure 12(a) we plot the current for  2 = 470 nm which corresponds approximately to the minimum voltage reflection coefficient of the curves in Figure 11(a).In Figure 12(b) we plot the current for  2 = 500 nm, which corresponds to the minimum voltage reflection coefficient of the curves of Figure 11(b).
We observe that the amplitude of the input current in the nanocircuit is in general in accordance with the results presented in Figure 9 for the isolated receiving dipole, where for larger  1 and  ℎ1 we have higher current along the circuit for both resonances.However, for fixed  ℎ1 = 10 nm the first  resonance (Figure 12(a)) has the smaller input current amplitude as compared to the second resonance (Figure 12(b)).For  ℎ1 = 15 nm, the first resonance (Figure 12(a)) presents the same level of the input current amplitude as that of the second resonance (Figure 12(b)).And for  ℎ1 = 20 nm the first resonance (Figure 12(a)) presents higher input current amplitude than that of the second resonance (Figure 12(b)).However, the cases of Figure 12(b) give a better impedance matching (lower values of |Γ  |) than the cases of Figure 12(a).
We conclude with this analysis that a good impedance matching does not necessarily mean a good efficiency in the receiving dipole, that is, higher input current amplitude.For example, the case of  ℎ1 = 20nm and  1 = 176.6 nm (Figure 12(a)) presents a higher intensity of the input current but a smaller voltage reflection coefficient than the case of  ℎ1 = 20 nm and  1 = 637.7 nm (Figure 12(b)).The total electric field distribution at plane  = 30 nm of these two opposite situations is shown in Figure 13.We can see in this figure the differences in the field intensity, showing higher values in the emitting dipole for case of the first resonance (Figure 13(a)).

Conclusions
In this paper, we presented a theoretical analysis of the impedance matching and excitation characteristics of a plasmonic optical nanocircuit fed by a Gaussian beam, which is composed of a receiving and an emitting dipole connected by a two-wire optical transmission line.The method of moments with equivalent surface impedance of the conductors was used for numerical calculations.
In the impedance matching analysis, we presented a parametric study of the reference circuit for different geometrical dimensions such as line and emitting dipole radius, gap distance of the line, and length of emitting dipole.We observed that better impedance matching is obtained when the radius of the emitting dipole is larger than the radius of the OTL.We also noted that the impedance matching can  be optimized varying the gap distance of the line and length of emitting dipole, and there are many possible solutions for good impedance matching.With respect to the attenuation losses in the line, we observed that a good impedance matching does not necessarily mean a good transmission efficiency along the OTL.In other words, we can have circuits with the same degree of impedance matching but with different transmission efficiency, that is, with different attenuation losses.This occurs because the efficiency is mainly a function of the radius of the OTL, where for a smaller radius the conductors possess a larger equivalent surface impedance and consequently higher conduction loss.
In the excitation analysis, our results show that receiving dipoles with larger radius produce a better impedance matching and a higher input current in the circuit.However, good reception efficiency of the receiving dipole does not necessarily mean good impedance matching in the nanocircuit.In general, receiving dipole in the first resonance presents better reception efficiency but poor impedance matching in comparison with that of the second resonance.This means that to design an efficient circuit one has to make a compromise between reception efficiency, transmission efficiency, and impedance matching.However, in terms of efficiency, the best circuit that we have obtained is that presented in International Journal of Antennas and Propagation Figure 13(a), where the level of the signal that arrives in the emitting dipole is the maximum one obtained.This situation is only accomplished when the receiving dipole operates in the first resonance and the emitting dipole operates in the second resonance.
The analysis presented here can be useful as guidelines to design efficient plasmonic optical nanocircuits for applications in integrated nanophotonics for optical signal data processing.In future works, we will investigate the overall efficiency of the nanocircuit and consider some other

Figure 5 :
Figure 5: Normalized distribution of the electric field amplitude at plane  = 30 nm: (a) incident field of the Gaussian beam abs(   ), (b) scattered field abs(   ), and (c) total field abs(  ).

Figure 9 :
Figure 9: Variation of current amplitude in middle of nanorod (  ) in function of total length  1 and radius  ℎ1 .

Figure 12 :Figure 13 :
Figure 12: Current distribution along nanocircuit for different receiving dipole in different resonances.(a) Results for the first resonance (/2), where the parameters used in these curves are as follows:   =15 nm,  ℎ2 = 20 nm, and  2 = 470 nm.(b) Results for the second resonance (3/2), where the parameters used in these curves are as follows:   = 15 nm,  ℎ2 = 20 nm, and  2 = 500 nm.
[27]2 ℎ1 , and Δ  > 2  .Then the current in each segment is approximated by sinusoidal basis functions[27].The expansion constants   are shown in Figure2(b)where each constant defines one triangular sinusoidal current.To determine these constants, we use  =   − 2 rectangular pulse test functions with unitary amplitude and perform the