1. Introduction

Graphene plasmons (GPs), the intrinsic collective oscillations of electrons coupled with electromagnetic waves in doped graphene, have attracted enormous interests for their unique properties, including such as inherently highly controllable, long-lived and extremely electromagnetic field confinement and enhancement in mid-infrared and terahertz spectral regimes [1–3]. Since theoretically proposed by Hanson [4] and Jablan et al [5], GPs have been widely studied for electro-optical modulation [6,7], quantum plasmonics [8], light harvesting [9], transformation optics [10] and infrared biosensors [11]. Due to large wavevector mismatch between GPs and free light [12,13], propagating GPs are usually excited by deep sub-wavelength antennas such as a linear antenna [14], cross shape antennas [15], a column of nanoantenna pairs [16] and carbon nanotubes [17]. However, the propagation direction of excited GPs is usually isotropic in the graphene plane as a result of the symmetry of structures and excitation configurations. And the unidirectional launching of GPs is still unsolved problem although it is important in ultracompact plasmonic devices. On the other hand, to achieve highly directional launching of surface waves in metal, attempts have been made to break the symmetry by introducing oblique incidence [18, 19], double slits [20], metal nanoantenna arrays [21], magneto-optical substrate [22] and circularly polarized incident waves [23,24], etc. Among the numerous reported methods, the use of near field interference of a circularly polarized wave is cornerstone for active switching due to the spin-orbit coupling of light [25,26], along with very high extinction ratio between different directions. In these inhomogeneous plasmonic fields, the right- and left-handed spin states and orbital properties become strongly coupled with each other [27]. However, the current experimental effort to mimic a two-dimensional rotating dipole by oblique incidence is controversial due to the inevitable magnetic induction currents [28]. It is well-known that two orthogonally oriented dipoles can be induced by two orthogonally polarized incident plane waves. Aware of that a circularly polarized dipole can be efficiently mimicked by a nanowire illuminated by two orthogonal plane waves due to the extremely localization of GPs, it is natural for us to consider realizing unidirectional generation of GPs by a combined circularly polarized dipole in unstructured graphene. Furthermore, the combined dipoles can be separated in space compared to a single circularly polarized dipole, which provides us another degree of freedom to control the excitation of GPs.

In this study, two orthogonally oriented dipoles are employed to efficiently excite symmetric and anti-symmetric charge ordering modes in flat and arc graphene planes. As long as the constructive and destructive interferences of near-fields take place in different propagation directions, the unidirectional launching of GPs occurs. Due to the inherent phase difference between symmetric and anti-symmetric evanescent modes induced by x-polarized and z-polarized dipoles, the extra phase difference of π/2, e.g., circularly polarized dipoles, should be introduced. Moreover, when the circularly polarized dipole is decomposed into two linear polarized ones and separated on different sides of the graphene, the behaviors of induced charge distribution and handedness of the background field are opposite. In further, the circularly polarized dipole can be efficiently mimicked by a semiconductor nanowire illuminated by two orthogonally polarized plane waves in experiments, and the problem about magnetic induction currents induced by oblique incidence can be solved in our system. And our findings should be found applications in compact plasmonic circuits in mid-infrared and terahertz regimes.

2. Theoretical background

The phenomena of unidirectional excitation of GPs can be understood by considering a dipole placed at a subwavelength distance close to a free standing graphene sheet. Figure 1(a) illustrates the scheme employed in our design. A two-dimensional (2D) dipole with momentum p_2D = [p_x, p_z] is placed above a graphene sheet. A Cartesian coordinate system is chosen with the graphene sheet laying in z = 0 and the position of the dipole is (0, z_dip). Without loss of generality, the result can be extended to three-dimensional (3D) treatment directly [23]. The vector potential A induced by the dipole without graphene can be expressed as A(r) = −iωµ₀G(r′, r)p, where $G(\mathbf{r},{\mathbf{r}}^{\prime })=\frac{i}{4}{H}_0^{(1)}\left(k_0\left|\mathbf{r}-{\mathbf{r}}^{\prime }\right|\right)$ is the 2D Green’s function in free space [29, 30], and k₀ = ω/c is wavenumber in vacuum. The angular spectrum decomposition of the vector potential can be written as $\mathbf{A}(\mathbf{r})=\frac{\omega {\mu }_0\mathbf{p}}{4\pi }{\int }_{-\infty }^{\infty }\frac{\text{exp}\left(i{k}_z\left|z-z_{\text{dip}}\right|\right)}{k_z}\text{exp}(i{k}_{x}x)d{k}_x$, where $k_z={\left(k_0^2-k_x^2\right)}^{1/2}$ is the longitude wavenumber. Thus the magnetic field can be deduced as

 

Fig. 1 Schematic of directional excitation of graphene plasmons by a dipole source. The dipole source can be mimicked by a semiconductor nanowire.

Download Full Size | PDF

Next we turn to the case where the dipole is laid on top of graphene. Due to the extremely large wavenumber of GPs, i.e., k_x = k_spp ≫ k₀, $k_z={\left(k_0^2-k_x^2\right)}^{1/2}\approx i{k}_{\text{spp}}$ is satisfied. Then the induced magnetic field has a relation of $H_y^0(k_{\text{spp}},z)\propto p_z/p_x\mp i$. One can conclude that the complete interferences take place as long as p_z and p_x have equally modulus with phase difference of ±π/2. The contribution from the graphene can be included by introducing reflected and transmitted fields, which can be calculated via simply multiplying the individual angular spectrums with corresponding Fresnel coefficients r^p and t^p, respectively. When the dimensionless conductivity α = 2πσ/c in units of the fine-structure constant α₀ ≃ 1/137 was adopted, the Fresnel coefficients can be written as [12,13]:

Noting that r^p and t^p depend on the modulus of k_x only, naturally, one can divide the contributions of a circularly polarized dipole into two parts, i.e., ${\rho }_s^{\mathit{ind}}(k_x)={\rho }^{p_z}(k_x)+{\rho }^{p_x}(k_x)$. The former induced by z-polarized dipole (abbreviated as p_z for convenience) satisfies ${\rho }^{p_z}(k_x)\propto p_z{k}_x^2/k_z$, while the latter induced by x-polarized dipole satisfies ρ^p_x(k_x) ∝ p_xk_x. The fundamental mechanism for directional generation is the charge density induced by p_z has an even parity both in angular spectrum and real space [31], whereas the opposite holds true for a p_x dipole. The superposition of ρ^p_z and ρ^p_x (not the H_y with opposite parities) leads to the constructive and destructive interferences in different directions.

When the p_z or p_x dipole is moved to its image position (0, −z_dip), noting that H̃_y(k_x, 0⁻) − H̃_y(k_x, 0⁺) = H₀(k_x, 0) (1 + r^p − t^p) = H_y(k_x, 0⁺) − H_y(k_x, 0⁻) and the minus sign should be adopted before p_x in the expressions of H_y thanks to z− (−z_dip) > 0, therefore the charge density can be written as

There are a lot of parameters to quantitative describe the asymmetrical excitation. Among them, the angular spectrum ratios satisfy $R_k[F]\equiv \left|F(k_x)/F(-k_x)\right|=\left|\left(\frac{p_z}{p_x}\frac{k_x}{k_z}+1\right)/\left(-\frac{p_z}{p_x}\frac{k_x}{k_z}+1\right)\right|$, (F ∈ {H_y, E_z, E_x, ρ}) which depend on (p_z/p_x)(k_x/k_z) only. Moreover, the spatial dependent near field ratios defined as R_x[F] ≡ |F(x)/F(−x)|, (F ∈ {H_y, ρ, P_x}) are also very important for GPs which can be obtained via full-field simulations and verified by near field experiments directly. In engineering, another important parameter to quantify the asymmetrical transmission is extinction ratio, which is defined as the logarithm of energy flux ratio in opposite directions η = 10log(P_r/P_l). The right and left energy flux along the graphene can be obtained by integrating the relative Poynting vector along z direction far from the dipole source,

3. Directional Generation of GPs by a single dipole

We present several scenarios in which the proper choices of dipoles close to graphene sheet provide possibilities for directional excitation of GPs. First of all, we consider the basic model described in Sec 2, and compare the simulated results to the analytical results calculated from angular spectra. In the simulation, the frequency of electromagnetic field emitted by the dipole is 30 THz (corresponding to λ ≈ 10 μm). The dipole is situated at a distance of d = 0.01λ on top of the graphene plane and has a momentum of [1, p_z/p_x]p_x, where the unit length momentum ṗ_x = −iωp_x = 1A · m is adopted for convenience, and the ratio p_z/p_x is discussed later. Chemical potential of the doped graphene is set as μ=0.4 eV, and ambient temperature is set as T=300 K, the in-plane conductivity of the graphene is computed within the local-random phase approximation (RPA) [32,33] with an intrinsic relaxation time τ = 120 fs (indicating the mobility of μ = 3000 cm²/Vs), which is a typical parameter derived from experiments [11, 34, 35]. Commercial software COMSOL Multiphysics based on FEM method is adopted to solve the Maxwell equations. From the dispersion relations of GPs, the wavenumber of GPs for this free-standing graphene is $k_{\text{spp}}=\sqrt{1-{\alpha }^2}\hspace{0.17em}\hspace{0.17em}{k}_0\approx 21.78\hspace{0.17em}\hspace{0.17em}{k}_0$, indicating the plasmon wavelength λ_spp is 459 nm, and the extinction parameter is a^* ≡ k_z/k_x|_kx=Re{kspp} = 0.9989i. Therefor p_z/p_x = ±0.9989i will lead to completely destructive interference of the excited GPs in a certain direction.

The simulated distribution of H_y field |Re{H_y}| is depicted in Fig. 2(a). Clearly, the plasmon mode is unidirectionally excited by the circularly polarized dipole, which has a much larger amplitude along +x than −x. Besides, one can find that the background field is anticlockwise rotational due to the individual rotational dipole source. Thus the angular momentum density of the background field is $\mathbf{L}={\epsilon }_0\mathbf{r}\times (\mathbf{E}\times \mathbf{B})=\frac{1}{c^2}\mathbf{r}\times \mathbf{S}$, which is along −y direction, where S is the Poynting vector [36,37]. Moreover, the angular momentum direction of preferred exited GPs is consistent with the angular momentum direction of background fields due to the conservation of angular momentum. When the phase difference between p_x and p_z changes from π/2 to −π/2 or placing the rotational polarized dipole on bottom of the graphene plane, the directions of angular momentum as well as the propagating GPs inverse. Therefore one can determine the preferred excitation directions simply by the direction of angular momentum. On the other hand, the directional excitation of GPs in real space can be understood by the asymmetry of the angular spectra in different directions. The angular spectra of initial, reflected and total magnetic fields are shown in Fig. 2(b). One can find the angular spectra have constructive and destructive interferences at +k_x and −k_x, respectively. To understand this effect, the Fresnel coefficient r^p is considered within plasmon pole approximation (k_x ≃ ±k_spp). The coefficient $r^p\propto {\left(k_x^2-k_{\text{spp}}^2\right)}^{-1}$ has peaks when k_x = ±Re{k_spp}. Specifically, when p_z/p_x equals to a^*, $H_y^0(k_{\text{spp}},0^{+})\propto \frac{p_z}{p_x}{\frac{k_x}{k_z}|}_{k_x=\text{Re}\left\{k_{\text{spp}}\right\}}+1=2$, this is a peak due to constructive interference and the intensity is twice as large as the magnitude of magnetic field induced by p_x or p_z individually. In the same manner, the valley exists due to constructive interference occurs at $\left|H_y^0\left(-k_{\text{spp}},0^{+}\right)\right|=0$. As a result, the angular spectra of initial, reflected and total fields have peaks at +k_spp as well as valleys at −k_spp for p_z/p_x = a^*. As to an ideal rotational polarized dipole, e.g. p_z/p_x equals to 1i, and $H_y^0\left(k_{\text{spp}},0^{+}\right)\propto k_x/{k_z|}_{k_x=\text{Re}\left\{k_{\text{spp}}\right\}}+1$ is about 2 owing to a^* ≈ 1i. However, there is a remarkable difference near k_x = −k_spp, where it is a peak instead of valley for the ideal circularly polarized dipole. This comes from that |r^p| is maxima at −k_spp and $\left|H_y^0(-k_{\text{spp}},0^{+})\right|\propto \left|-i/a^{*}+1\right|$ is a slowly varying quantity. In further, the spatial distribution of charge density is a vital physical quantity to describe the collective oscillations, such as plasmons. In Fig. 2(c), the simulated charge density distribution in the graphene plane is compared to the analytical result from Eq. (5). They are in perfect agreement and have apparently constructive and destructive interferences in x > 0 and x < 0, respectively. The spatial dependent near field ratios of exited GPs are shown in Fig. 2(d). We can see that the near field ratios of reflected magnetic, charge density and energy flux are over 100 for x < 3λ_spp.

 

Fig. 2 Directional excitation of GPs by a circularly polarized dipole. (a) Magnetic field distributions for GPs excited by a 2D circularly polarized dipole p_2D = [1, a^*]p₀, where −iωp₀ = 1A · m. (b) Angular momentum spectra of initial ( $H_y^0$), reflective ( $H_y^{\text{ref}}$) and total (H_y) magnetic field magnitude of the polarized dipole (solid lines) and an ideal circularly polarized dipole (dashed lines). The dotted lines indicate the wavenumbers of GPs. (c) The simulated (solid line) and analytically calculated (marked by red circle) spatial dependent charge density in graphene. (d) Spatial dependent near field ratio $R_x\left[H_y^{\text{ref}}\right]=\left|H_y^{\text{ref}}(x)/H_y^{\text{ref}}(-x)\right|$(colored in black), R_x[ρ] = |ρ(x)/ρ(−x)|(colored in red) and R_x[P_x] = |P_x(x)/P_x(−x)|(colored in blue).

Download Full Size | PDF

4. Directional Generation of GPs by two dipoles

A circularly polarized dipole can be decomposed into two orthogonal polarized dipoles, thus it is interesting to see what will happen if these dipoles are placed on the different sides of graphene. Specifically, two orthogonal polarized dipoles with momenta p_x = [1, 0]p₀ and p_z = [0, a^*]p₀ are considered, where they are located at (0, 0.01λ) and (0, −0.01λ), respectively. The simulated results are shown in Fig. 3. One can find that the rotational direction of background field in Fig. 3(a) is the same as that in Fig. 2(a), however, the propagation direction of excited GPs inverses. This result means that one cannot distinguish these two cases from the far fields excepted for the preferred propagation direction of excited GPs. Considering that the H_y is discontinuous across the graphene plane, the $H_y^{\text{ref}}$ changes from upside to downside of graphene when the dipole is placed to its image position. Meanwhile, $H_y^0$ and the induced charge ρ remain unique. Similar to the results in Fig. 2(b), the angular spectra of $H_y^0$ and ρ are also shown in Fig. 3(b) to demonstrate the mechanism of directional excitation in this scenario. Remarkably, the preferred directions of background field and induced charge are opposite, which are along +x and −x, respectively. This result is rather counterintuitive and means that the angular momentum is ’non-conservation’ at first sight. Actually, one can find this puzzling result comes from the magnetic field discontinuity at upper and lower sides of graphene (r ≠ t − 1), and the angular momentum is conservation. In our proposed scheme, the angular momentum is zero at the original point with both positive and negative signs in z = 0 plane simultaneously. Thus the use of the conservation of angular momentum cannot determine the preferred direction of GPs directly. To demonstrate the physical factor to determine the preferred directions, we turn to see the dependence of directional generation on phase difference Δϕ = ϕ^p_z − ϕ^p_x = arg{p_z/p_x}. Similar to the superposition of polarizations, these two opposite sense of rotations lead to a classification of vibration ellipses according to their handedness, which is decided by the phase difference of two vibration vectors. If the phase difference satisfies Δϕ = mπ, m = 0, 1, 2···, the superpositions are linear polarized dipoles, and their angular momenta is zero due to r//S, thus the excited GPs should be isotropic without any other asymmetry to fulfill the conservation of angular momentum. When the phase difference satisfies Δϕ = π/2±2mπ, the near fields rotates in the anticlockwise sense, it is said to be left-handed. If extra π phase is introduce to the Δϕ, the handedness and preferred direction of excitation will change. In the proposed system, there are three important factors to determine the handedness and preferred direction. The first one is the initial phase difference Δϕ which is from the dipoles themselves, e.g., if the initial arg{p_z/p_x} changes from π/2 to π/2 ± π, the handedness, i.e., rotational direction of background field and preferred direction will change. The second factor is the dipole position relative to the graphene. From the relation $H_y=H_y^{p_z}+H_y^{p_x}=i/\mu \left[k_x{A}_z+\mathbf{sgn}\left(z_{\text{dip}}-z\right)k_z{A}_x\right]$, one can known that moving p_x dipole to its image position will introduced a minus sign due to the sign function, which is equivalent to introduce extra π phase difference when talking about the handedness and the initial electromagnetic field, while there is no extra phase difference when moving the p_z dipoles to their image position. The total extra phase difference from aforementioned two factors will determine the handedness and the preferred direction of initial field. However, they are insufficient to determine the preferred direction of exited GPs. Noting that the scattering field of upper and lower sides of graphene satisfied −r = t − 1 for free standing graphene, which introduced a minus sign compared to the continuous boundary condition r = t − 1. This is the last vital factor to determine the preferred direction of induced field. Actually, the aforementioned counterintuitive result is originated from the minus sign, which can not be treated as extra phase difference as before because it only acts on induced field and do not affect the handedness and the distribution of initial magnetic field. In a word, there are three factors for p_x to affect the preferred direction of excited GPs, while only two factors for p_z to affect the preferred direction of excited GPs in our considered system. When the p_x and p_z locate in the same side of graphene, $H_y^{\text{ref}}(z=0)$ always denotes the magnetic field in the dipole side, thus the angular momentums of $H_y^0$ and ρ have the same sign owing to the conservation of angular momentum. When they are in the different side, $H_y^{\text{ref}}$ in z = 0 induced by the two dipoles denotes different sides of graphene, and the preferred direction of GPs should be decided by ρ rather than $H_y^0$. Due to the extra minus sign, the preferred direction of ρ and $H_y^0$ is always opposite in this condition. The spatial dependent induced charge density is plotted in Fig. 3(c). One can find that the charge oscillates only in the −x direction which is in good agreement with the analytical result. The comparison of spatial dependent extinction ratio of a single circularly polarized dipole (named after case I) and two orthogonal polarized dipoles placed at both sides of the graphene (named after case II) is shown in Fig. 3(d). One can see the extinction ratio is over 20 for x < 4λ_spp, and the ratio in case II is less than that in case I, this originates from the opposite preferred directions of the initial and induced fields in case II. The difference on extinction ratio between these two cases can be ignored when x ≥ 5λ_spp.

 

Fig. 3 Directional excitation of GPs by two separated orthogonal polarized dipoles. (a) Magnetic field distributions for GPs excited by two orthogonal polarized dipoles p_x = [1, 0]p₀ and p_z = [0, a^*]p₀ located at (0, 0.01λ) and (0, −0.01λ), respectively. The insert figure shows the enlarged excitation region. (b) Angular momentum spectra of initial magnetic field magnitude $H_y^0$ and charge density ρ (multiplied by c) in the graphene plane, the initial and induced quantities have opposite preferred direction for excitation of GPs. (c) The simulated (solid line) and analytically calculated (marked by red circle) spatial dependent charge density in graphene. (d) The spatial dependent extinction ratio for a single circularly polarized dipole (case I, colored in black), and two separated orthogonal polarized dipoles (case II, colored in red).

Download Full Size | PDF

5. Directional generation of GPs on a curved system

The curvature breaks the mirror symmetry relative to graphene plane, then the induced radiative loss will affect the excitation and propagation efficiency of GPs. In infinite long graphene cylinder system, the critical curvature radius which permits existence of confined wave can be calculated by $r_c\approx \frac{k_0}{\text{Im}\left\{k_z\right\}\left(k_{\text{spp}}-k_0\right)}\approx 0.048{\lambda }_{\text{spp}}$ [38], thus a circular radius of the graphene cylinder as r = 1.4 μm ≈ 3λ_spp is chosen, which is a typical value in flexible transformation plasmonics [38]. A circularly polarized dipole with momentum as the same as in flat graphene is placed above the infinite long graphene cylinder at a distance of 100 nm. The configuration and simulated H_y field amplitude are shown in Fig. 4(a). Remarkably, the mode propagates mainly along clockwise direction. which is in coincide with the result in flat graphene (case I). We can describe the induced charge density by

 

Fig. 4 Directional excitation of GPs on a free-standing infinite long graphene cylinder. The distribution of magnetic field |Re{H_z}| for GPs excited by the configuration of case I(a) and two separated orthogonal polarized dipoles (c). Simulated and analytically calculated spatial dependent charge density for the case I(b) and case II(d). The phase in (c) is set as π/4 in order to show the two individual sources.

Download Full Size | PDF

6. Experimental proposals

The dipole employed in the paper can be mimicked by a semiconductor nanowire illuminated by two orthogonally polarized plane waves. Due to the wave interference, the background standing wave satisfies $E_x\propto \sqrt{2}\text{cos}(k_{x}x+\pi /4)$, $E_z\propto \sqrt{2}i\text{sin}(k_{x}x+\pi /4)$ and p_z/p_x = E_z/E_x = itan(k_xx+ π/4), thus the amplitude condition of an ideal circular dipole requires k_xD ≪ π/4, where D is the dimension of nanowire. In the configuration of directional excitation of metallic plasmons, $k_{x}D~1/\sqrt{2}$, the induced dipole moment is too nonuniform to be treated as an ideal circular dipole source. This means that the method is not suitable for directional excitation of metallic plasmons. However, there is no limitation any more in excitation of GPs due to the deep subwavelength of the nanowire size in infrared spectrum, i.e., $k_{x}D\ll 1/\sqrt{2}$. In our simulation, an In_0.53Ga_0.47As nanowire with diameter D of 100 nm (0.01λ) is used to mimic the rotational polarized dipole. The Drude model was adopted to model the dielectric constant of In_0.53Ga_0.47As. In this model, the dielectric function is given by $\epsilon (\omega )={\epsilon }_{\infty }-{\omega }_p^2/\left[\omega (\omega +i\gamma )\right]$, where ${\omega }_p={\left(n{q}^2/m_e^{*}{\epsilon }_0\right)}^{1/2}$ is the plasma frequency, ε_∞ is the high frequency dielectric constant, and $\gamma =q/{\mu }_e{m}_e^{*}$ is damping rate. Extracted from the reference in [39], the parameters of In_0.53Ga_0.47As are ε_∞ = 12.15, τ = γ⁻¹ = 0.1 ps, and m^* = 0.523m_e. Moreover, n = 6.3 × 10¹⁸/cm³ is used to realize the resonance of the nanowire at 30 THz. The absorption cross length normalized to geometry cross length for different diameters of the nanowire is shown in the insert figure of Fig. 5(a), the absorption peak lays at 30 THz and is independent on the diameter of the nanowire because the electrostatic approximation is satisfied. Two time harmonic orthogonal incident plane waves with amplitude of 1 V/m and phase difference of π/2 are taken to illuminate the wire, the schematic and simulated electric field distributions are depicted in Fig. 5(a), where the incident field has been subtracted from the total field. The expression of the incident fields adopted in the simulation is expressed as

 

Fig. 5 Directional excitation of GPs using In_0.53Ga_0.47As nanowires illuminated by two orthogonally polarized plane waves. (a) Schematics representation and electric field distributions |E_z| for excited GPs, the insert figure is normalized absorption cross section of the In_0.53Ga_0.47As nanowire with radii of 20 nm, 50 nm and 100 nm. (b) The dependence of the angular spectrum ratios R_k(k_spp) on extinction parameter |p_z/p_x|. The blue circle indicates the parameter of the considered nanowire. The insert figure is the polarized circle of the nanowire with (solid line) and without (broken line) graphene sheet. (c) The simulated and theoretically calculated spatial dependent charge density in graphene. The thick line indicates the charge distribution induced by ideal circular polarized dipole with dipole momentum as −iωp_x = 76.72 pA · m, the thin line indicate the case of dipole momentum as −iωp_x = 61.08 pA · m and p_z/p_x = −0.0624 + 0.9953i, respectively. The dot marked line indicates the simulated result of semiconductor nanowire. (d) The energy flux ratio of ideal dipole (p_z/p_x = i, solid line), actual dipole(p_z/p_x = −0.0624 + 0.9953i, broken line) and simulated result of semiconductor nanowire (thick line).

Download Full Size | PDF

7. Conclusion

We have demonstrated here that a circularly polarized dipole and two mirror image symmetric dipoles with orthogonally polarizations can directional generate propagating GPs due to the near field interference. The viewpoint of angular momentum conservation is very efficient to determine the preferred propagation direction of exited GPs. When the dipoles are laid in different sides of graphene, the spatial charge density rather than the magnetic field should be adopted to analysis the excited GPs due to the extra minus sign from the discontinuous of magnetic field. On this occasion, the magnetic field of dipole and induced charge distribution have opposite preferred directions and the properties of excited GPs should be described by the behavior of induced charge. Moreover, the direction generation of GPs can be extended into arc surface directly. Furthermore, a semiconductor nanowire can be regarded as a localized source to mimic the polarized dipoles, which can be realized in real experiments.