Resonance frequencies for the first four eigenmodes (), nm, eV, meV, , , , and .
This chapter discusses various approaches for calculating the modified Mie-Lorenz coefficients of the graphene-based multilayered cylindrical and spherical geometries. Initially, the Kubo model of graphene surface conductivity is discussed. Then, according to it, the formulations of scattering from graphene-based conformal particles are extracted. So, we have considered a graphene-wrapped cylinder and obtained its scattering coefficients by considering graphene surface currents on the shell. Later, a layered nanotube with multiple stacked graphene-dielectric interfaces is introduced, and for analyzing the plane wave scattering, graphene surface conductivity is incorporated in the transfer matrix method (TMM). Unlike the previous section, the dielectric model of graphene material is utilized, and the boundary conditions are applied on an arbitrary graphene interface, and a matrix-based formulation is concluded. Then, various examples ranging from super-scattering to super-cloaking are considered. For the scattering analysis of the multilayered spherical geometries, recurrence relations are introduced for the corresponding modified Mie-Lorenz coefficients by applying the boundary conditions at the interface of two adjacent layers. Later, for a sub-wavelength nanoparticle with spherical morphology, the full electrodynamics response is simplified in the electrostatic regime, and an equivalent circuit is proposed. Various practical examples are included to clarify the importance of scattering analysis for graphene-based layered spheres in order to prove their importance for developing novel optoelectronic devices.
- modified Mie-Lorenz theory
- 2D material
Cylindrically layered structures have various exotic applications. For instance, a metal-core dielectric-shell nano-wire has been proposed for the cloaking applications in the visible spectrum. The functionality of this structure is based on the induction of antiparallel currents in the core and shell regions, and the design procedure is the so-called scattering cancelation technique . Experimental realization of a hybrid gold/silicon nanowire photodetector proves the practicality of these structures . As an alternative approach for achieving an invisible cloak, cylindrically wrapped impedance surfaces are designed by a periodic arrangement of metallic patches, and the approach is denominated as mantle cloaking . Conversely, cylindrically layered structures can be designed in a way that they exhibit a scattering cross-section far exceeding the single-channel limit. This phenomenon is known as super-scattering and has various applications in sensing, energy harvesting, bio-imaging, communication, and optical devices [4, 5]. Moreover, a cylindrical stack of alternating metals and dielectrics behaves as an anisotropic cavity and exhibits a dramatic drop of the scattering cross-section in the transition from hyperbolic to elliptic dispersion regimes [6, 7]. The Mie-Lorenz theory is a powerful, an exact, and a simple approach for designing and analyzing the aforementioned structures.
Multilayered spherical structures have also attracted lots of interests in the field of optical devices. A dielectric sphere made of a high index material supports electric and magnetic dipole resonances which results in peaks in the extinction cross-section . Moreover, by covering the dielectric sphere with a plasmonic metal shell, an invisible cloak is realizable, which is useful for sensors and optical memories . By stacking multiple metal-dielectric shells, an anisotropic medium for scattering shaping can be achieved .
From the above discussions, it can be deduced that tailoring the Mie-Lorenz resonances in the curved particles results in developing novel optical devices. In this chapter, we are going to extend the realization of various optical applications based on the excitations of localized surface plasmons (LSP) in graphene-wrapped cylindrical and spherical particles. To this end, initially we introduce a brief discussion of modeling graphene material based on corresponding surface conductivity or dielectric model. Later, we extract the modified Mie-Lorenz coefficients for some curved structures with graphene interfaces. The importance of developed formulas has been proven by providing various design examples. It is worth noting that graphene-wrapped particles with a different number of layers have been proposed previously as refractive index sensors, waveguides, super-scatterers, invisible cloaks, and absorbers [11, 12, 13, 14, 15]. Our formulation provides a unified approach for the plane wave and eigenmode analysis of graphene-based optical devices.
Graphene is a 2D carbon material in a honeycomb lattice that exhibits extraordinary electrical and mechanical properties. In order to solve Maxwell’s equations in the presence of graphene, two approaches are applied by various authors, and we will review them in the following paragraphs. It should be noted that although we are discussing the graphene planar model, we will use the same formulas for the curved geometries when the number of carbon atoms exceeds 104, letting us neglect the effect of defects. Therefore, the radii of all cylinders and spheres are considered to be greater than 5 nm . Moreover, bending the graphene does not have any considerable impact on the properties of its surface plasmons, except for a small downshift of the frequency. Figure 1 shows the propagation of the graphene surface plasmons on the S-shaped and G-shaped curves .
Since graphene material is atomically thin, in order to consider its impact on the electromagnetic response of a given structure, boundary conditions at the interface can be simply altered. To this end, graphene surface currents that are proportional to its surface conductivity should be accounted for ensuring the discontinuity of tangential magnetic fields. In the infrared range and below, we can describe the graphene layer with a complex-valued surface conductivity which may be modeled using the Kubo formulas [18, 19]. The intraband and interband contributions of graphene surface conductivity under local random phase approximations read as :
The parameters , , , , and are reduced Plank’s constant, electron charge, chemical potential, charge carriers scattering rate, and temperature, respectively. The above equations are valid for the positive valued chemical potentials. Moreover, graphene-based structures can be analyzed by assigning a small thickness of to the graphene interface and later approaching it to the zero. In this method, by defining the volumetric conductivity as and using it in Maxwell’s curl equations, the equivalent complex permittivity of the layer can be obtained as :
where subscripts and represent the real and imaginary parts of the surface conductivity, respectively. Both models will be used in the following sections.
Figure 2(a) and (b) shows the real and imaginary parts of graphene surface conductivity at the temperature of T = 300°K. The real part of the conductivity accounts for the losses, while the positive valued imaginary parts represent the plasmonic properties . Moreover, the real and imaginary parts of the graphene equivalent bulk permittivity are shown in Figure 2(c) and (d). The negative valued real relative permittivity represents the plasmonic excitation, and the imaginary part of the permittivity represents the losses . It should be noted that all of the formulas of this chapter are adapted with time-harmonic dependency.
2. Graphene-coated cylindrical tubes
In this section, the modified Mie-Lorenz coefficients of a single-layered graphene-coated cylindrical tube will be extracted. The formulation is expanded into the multilayered graphene-based tubes through exploiting the TMM method, and later, various applications of the analyzed structures, including emission and radiation properties, complex frequencies, super-scattering, and super-cloaking, will be explained.
2.1 Scattering from graphene-coated wires
Let us consider a graphene-wrapped infinitely long cylindrical tube. The structure is shown in Figure 3(a), and it is considered that a TEz-polarized plane wave illuminates the cylinder. In general, TE and TM waves are coupled in the cylindrical geometries. For the normally incident plane waves, they become decoupled, and they can be treated separately. For simplicity, we consider the normal incidence of plane waves where the wave vector is perpendicular to the cylinder axis.
In order to obtain the modified Mie-Lorenz coefficients, the incident, scattered, and internal electromagnetic fields are expanded in terms of cylindrical coordinates special functions which are, respectively, the Bessel functions and exponentials in the radial and azimuthal directions. In order to exploit a terse mathematical notation, the vector wave functions are introduced as :
The complete explanation of the above vector wave functions and their self and mutual orthogonally relations can be found in the classic electromagnetic books . In the above equation, is the solution of the Bessel differential equation, and is its order. It is clear that in the environment, the radial field contains the Hankel function of the first kind in order to account for the radiation condition at infinity, while in the medium region, the Bessel function is utilized to satisfy the finiteness condition at the origin of the structure.
In the graphene-based cylindrical structures, the plasmonic state is achieved via illuminating a TEz wave to the structure. Therefore, for the normal illumination, the incident, scattered, and dielectric electromagnetic fields are shown with the superscripts , respectively, and they read as :
where and are the magnitudes of the incident electric and magnetic fields, respectively, and they are related via the intrinsic impedance of the free space. The wavenumber in the region is denoted by . The coefficients are, respectively, , , and for the incident, scattered, and core regions. Moreover, and bn are the well-known Mie-Lorenz coefficients, which are called the modified Mie-Lorenz coefficients for the scattering analysis of graphene-based structures.
The boundary conditions at the graphene interface at are the continuity of the tangential electric fields along with the discontinuity of tangential magnetic fields. Therefore:
By applying the boundary conditions in the expanded fields, the linear system of equations for extracting the unknowns can be readily obtained. The solution of the extracted equations for the scattering coefficients leads to:
The same procedure can be repeated for the TMz illumination. The normalized scattering cross-section (NSCS) reads as:
where the normalization factor is the single-channel scattering limit of the cylindrical structures. In order to have some insight into the scattering performance of graphene-wrapped wires, the scattering efficiency for ε1 = 3.9 and μc = 0.5 eV is plotted in Figure 3(b) by varying the radius of the wire. As the figure illustrates, a peak valley line shape occurs in each wavelength. They correspond to invisibility and scattering states and will be further manipulated in the next sections to develop some novel devices. The excitation frequency of the plasmons is the complex poles of the extracted coefficients  which will be discussed in the next subsection. Interestingly, the scattering states of graphene-coated dielectric cores are polarization-dependent. By using a left-handed metamaterial as a core, this limitation can be obviated .
2.1.1 Eigenmode problem and complex frequencies
As in any resonant problem, additional information can be obtained by studying the solutions to the boundary value problem in the absence of external sources (eigenmode approach). Although, from a formal point of view, this approach has many similar aspects with those developed in previous sections, the eigenmode problem presents an additional difficulty related to the analytic continuation in the complex plane of certain physical quantities. Due to the fact that the electromagnetic energy is thus leaving the LSP (either by ohmic losses or by radiation towards environment medium), the LSP should be described by a complex frequency where the imaginary part takes into account the finite lifetime of such LSP. The eigenmode approach is not new in physics, but its appearance is associated to any resonance process (at an elementary level could be an RLC circuit), where the complex frequency is a pole of the analytical continuation to the complex plane of the response function of the system (e.g., the current on the circuit). Similarly, in the eigenmode approach presented here, the complex frequencies correspond to poles of the analytical continuation of the multipole terms (Mie-Lorenz coefficients) in the electromagnetic field expansion.
In order to derive complex frequencies of LSP modes in terms of the geometrical and constitutive parameters of the structure, we use an accurate electrodynamic formalism which closely follows the usual separation of variable approach developed in Section 2.1. We can obtain a set of two homogeneous equations for the –th LSP mode :
where the prime denotes the first derivative with respect to the argument of the function and (). For this system to have a nontrivial solution, its determinant must be equal to zero, a condition which can be written as:
where , . This condition is the dispersion relation of LSPs, and it determines the complex frequencies in terms of all the parameters of the wire cylinder.
2.1.2 Non-retarded dispersion relation
When the size of the cylinder is small compared to the eigenmode wavelength, i.e., , where c is the speed of light in free space, Eq. (13) can be approximated by using the quasistatic approximation as follows. Using the small argument asymptotic expansions for Bessel and Hankel functions, the functions and . Thus, the dispersion relation (14) adopts the form:
Taking into account that in the non-retarded regime the propagation constant of the plasmon propagating along perfectly flat graphene sheet can be approximated by:
it follows that the dispersion relation (14) for LSPs in dielectric cylinders wrapped with a graphene sheet can be written as:
where is the LSP multipole order. The dispersion Eq. (17), known as Bohr condition, states that the –th LSP mode of a graphene-coated cylinder accommodates along the cylinder perimeter exactly oscillation periods of the propagating surface plasmon corresponding to the flat graphene sheet.
For large doping () and relatively low frequencies (), the intraband contribution to the surface conductivity (1), the Drude term, plays the leading role. In this case, the non-retarded dispersion equation Eq. (14) is written as:
which can be analytically solved for the plasmon eigenfrequencies,
where is the effective plasma frequency of the graphene coating. It is worth noting that the real part of is proportional to , and as a consequence, the net effect of the chemical potential increment is to increase the spectral position of the resonance peaks when the structure is excited with a plane wave or a dipole emitter [28, 29].
In the following example, we consider a graphene-coated wire with a core radius nm, made of a non-magnetic dielectric material of permittivity immersed in the vacuum. The graphene parameters are eV, meV, and K. Table 1 shows the first four eigenfrequencies calculated by solving the full retarded (FR) dispersion Eq. (14) (second column) and by using the analytical approximation (AA) given by Eq. (19) (third column). Since the radius of the wire is small compared with the eigenmode wavelengths, good agreement is obtained between the complex FR and AA values, even when the AA assigns to all multipolar plasmon modes.
2.2 Multilayered graphene-based cylindrical structures
In this section, multilayered cylindrical tubes with multiple graphene interfaces are of interest. In order to ease the derivation of the unknown expansion coefficients, matrix-based TMM formulation is generalized to the tubes with several graphene interfaces. Initially, consider a layered cylinder constructed by the staked ordinary materials under TEz plane wave illumination, as shown in Figure 4. The total magnetic field at the environment can be expressed as the superposition of incident and scattered waves as in Section 2.1. The unknown expansion coefficients of the scattered wave can be determined by means of the matrix defined as :
where C represents the core layer. In the above equation, the dynamical matrix of each region is constructed based on its constitutive and geometrical parameters distinguished through the subscript q (q = 1, 2, …, N). We have:
The argument of the above special functions is , and the TEz wave impedance equals . After generating matrix for the structure, coefficients can be calculated as:
In order to incorporate the graphene surface conductivity in the above formulas, let us consider each graphene interface as a thin dielectric with the equivalent complex permittivity defined in Eq. (3) and utilize the TMM formulation in the limiting case of a small radius at the graphene interface with the wave number of kg, i.e., . At each boundary, using the Taylor expansion as in the matrix, the graphene interface can be represented by the following matrices. We have:
where the free-space impedance equals 377 ohms. Once matrix is generated, the modified Mie-Lorenz coefficients and thus scattering cross-section are readily attainable. In the following subsections, the above equations will be used to design some novel optoelectronic devices.
2.2.1 Application in mantle cloaking
Widely tunable scattering cancelation is feasible by using patterned graphene-based patch meta-surface around the dielectric cylinder as shown in Figure 5. The surface impedance of the graphene patches can be simply and accurately calculated by closed-form formulas, to be inserted in the modified Mie-Lorenz theory .
2.2.2 Application in super-scattering
Let us consider a triple shell graphene-based nanotube under plane wave illumination, as shown in Figure 6(a). This structure is used to design a dual-band super-scatterer in the infrared frequencies. To this end, modified Mie-Lorenz coefficients of various scattering channels should have coincided with the proper choice of geometrical and optical parameters. In order to construct the Tn matrix for this geometry, one needs to multiply nine 2 × 2 dynamical matrices, which is mathematically complex for analytical scattering manipulation. Therefore, the associated planar structure, shown in Figure 6(b), is used to develop the dispersion engineering method as a quantitative design procedure of the super-scatter. The separations of the free-standing graphene layers are d1 = d2 = 45 nm in the planar structure, and the transmission line model is used to analyze it. Moreover, the chemical potential of lossless graphene material is μc = 0.2 eV in all layers. The dispersion diagram of the planar structure is illustrated in Figure 7(a), which predicts the presence of three plasmonic resonances in each scattering channel of the tube at around the frequencies that fulfill , where is the mean of the radii of all layers and β is the propagation constant of the plasmons in the planar structure. This condition is known as Bohr’s quantization formula , and its validity for our specific structure is proven by means of the previously developed formulas in Figure 7(b). Eqs. (12), (22), and (23) are used to obtain this figure.
In order to design a dual-band super-scatterer, the plasmonic resonances of two scattering channels have coincided by fine-tuning the results of the Bohr’s model. The optimized geometrical and constitutive parameters are Rc = 45.45 nm, d1 = 45.05 nm, d2 = 43.23 nm, ε1 = 3.2, ε2 = 2.1, ε3 = 2.2, and ε4 = 1. Figure 8 shows the NSCS and magnetic field distribution for the dual operating bands of the structure. It is clear that NSCS exceeds the single-channel limit by the factor of 4, and in the corresponding magnetic field, there is a large shadow around the nanometer-sized cylinder at each operating frequency. Other designs are also feasible by altering optical and geometrical parameters. Furthermore, the far-field radiation pattern is a hybrid dipole-quadrupole due to simultaneous excitation of the first two channels. It should be noted that an inherent characteristic of the super-scatterer design using plasmonic graphene material is extreme sensitivity to the parameters. Moreover, in the presence of losses, the scattering amplitudes do not reach the single-channel limit anymore, and this restricts the practical applicability of the concepts to low-frequency windows.
2.2.3 Application in simultaneous super-scattering and super-cloaking
As another example, the dispersion diagram of Figure 7(a) along with Foster’s theorem has been used to conclude that each scattering channel of the triple shell tube contains two zeros which are lying between the plasmonic resonances, predicted by the Bohr’s model. Later, we have coincided the zeros and poles of the first two scattering channels in order to observe super-scattering and super-cloaking simultaneously . The optimized material and geometrical parameters are εc = 3.2, ε1 = ε2 = 2.1, Rc = 45.45 nm, d1 = 46.25 nm, and d2 = 46.049 nm. The NSCS curves corresponding to the super-cloaking and super-scattering regimes are illustrated in Figure 9(a) and (b), as well as the expected phenomenon, is clearly observed. The corresponding magnetic field distributions, shown in Figure 9(c) and (d), also manifest the reduced and enhanced scatterings in the corresponding operating bands, respectively. Similar to the dual-band super-scatterer of the previous section, the performance of this structure is very sensitive to the optical, material, and geometrical parameters. By further increasing the number of graphene shells, other plasmonic resonances and zeros can be achieved for the manipulation of the optical response.
3. Graphene-coated spherical structures
In this section, multilayered graphene-coated particles with spherical morphology are investigated, and corresponding modified Mie-Lorenz coefficients are extracted by expanding the incident, scattered, and transmitted electromagnetic fields in terms of spherical harmonics. It is clear that by increasing the number of graphene layers, further degrees of freedom for manipulating the optical response can be achieved. For the simplicity of the performance optimization, an equivalent RLC circuit is proposed in the quasistatic regime for the sub-wavelength plasmons, and various practical examples are presented.
3.1 Multilayered graphene-based spherical structures
In this section, the most general graphene-based structure with N dielectric layers, as shown in Figure 10, is considered, and plane wave scattering is analyzed through extracting recurrence relations for modified Mie-Lorenz coefficients. It should be noted that since, in the TMM method, multiple matrix inversions are necessary, unlike the cylindrically layered structures of the previous section, the spherical geometries are analyzed through recurrence relations. Also, scattering from a single graphene-coated sphere has been formulated elsewhere , and it can be simply attained as the special case of our formulation.
The scattering analysis is very similar to that of the single-shell sphere , unless the Kronecker delta function is used in the expansions in order to find the electromagnetic fields of any desired layer with terse expansions. Therefore :
By considering as either or , which stand for the spherical Bessel and Hankel functions of the first kind with order n, respectively, and as the associated Legendre function of order (n, m), the vector wave functions are defined as follows:
where super-indices (1) in the vector wave functions show that the Hankel functions are used in the field expansions. The boundary conditions at the interface of adjacent layers read as:
Therefore, the linear system of equations resulting from the above conditions is:
where , , , and (d is defined as a symbol for the derivative with respect to the radial component). By rearranging the above equations, the coefficients of the layer (p + 1) can be written in terms of the coefficients of the layer p as:
where the sub/superscripts H and V represent the TE and TM waves, respectively. The directions of propagation of these waves are realized thought the subscripts F (outgoing waves) and P (incoming waves). The effective reflection coefficients are extracted as:
Moreover, it can be readily shown that the transmission coefficients read as:
where . By using , the recurrence relations can be started, and the field expansion coefficients in any desired layer can be obtained. The extinction efficiency is related to the external modified Mie-Lorenz coefficients via:
where symbol represents the real part of the summation. In order to verify the extracted coefficients, the extinction efficiencies of three graphene-coated structures is provided in Figure 11. In the graphical representation of the structures, the dashed lines illustrate graphene interfaces, while the solid line shows a PEC core. The optical and geometrical parameters are R1 = 200 nm, R2 = 100 nm, R3 = 50 nm, μc = 0.3 eV, T = 300° K, and τ = 0.02 ps. The analytical results are compared with the numerical results of CST 2017 commercial software, and good agreement is achieved. Moreover, the analytical formulation provides a fast and accurate tool for the scattering shaping of various spherical geometries.
In order to realize the priority of the closed-form analytical formulation with respect to the numerical analysis, the simulation times of both methods are included in Table 2. Considerable time reduction using the exact solution is evident. Moreover, since 3D meshing and perfectly matched layers are not required in this method, it is efficient in terms of memory as well.
3.1.1 Quasistatic approximation and RLC model
Based on the results of Section 3.1, the modified Mie-Lorenz coefficients of the graphene-based spherical particles form infinite summations in terms of spherical Bessel and Hankel functions. In general, graphene plasmons are excited in the sub-wavelength regime, and only the leading order term of the summation is sufficient for achieving the results with acceptable precision. In this regime, the polynomial expansion of the special functions can also be truncated in the first few terms . Later, the extracted modified Mie-Lorenz coefficients can be rewritten in the form of the polynomials. To further simplify the real-time monitoring and performance optimization of the graphene-coated nanoparticles, an equivalent RLC circuit can be proposed by representing the rational functions in the continued fraction form as :
The equivalent circuit corresponding to the above representation is shown in Figure 12.
The continued fraction representation for the TM coefficients is:
where , , , , , .
The elements of the equivalent circuit for the TM coefficients read as:
3.1.2 Application in emission
In order to illustrate the application of Mie analysis for the graphene-wrapped structures, let us consider vertical and horizontal dipoles in the proximity of a graphene-coated sphere, as shown in Figure 13. Although in the Mie analysis, the excitation is considered to be a plane wave, by using the scattering coefficients, the total decay rates can be calculated for the dipole emitters, and it can be proven that the localized surface plasmons of the graphene-wrapped spheres can enhance the total decay rate, which is connected to the Purcell factor [16, 37]. The amount of electric field enhancement for the radial-oriented and tangential oscillating dipoles with the distance of xd, respectively, read as:
Figure 13(b) shows the local field enhancement for the average orientation of the dipole emitter in the vicinity of the sphere with R1 = 20 nm, coated by a graphene material with the chemical potential of μc = 0.1 eV. As the figure shows, an enhanced electric field in the order of ∼104 is obtained for the dipole distance of 1 nm with averaged orientation, and it decreases as the dipole moves away from the sphere.
3.1.3 Application in super-scattering
The possibility of a super-scatterer design using graphene-coated spherical particles is illustrated in Figure 14. The design parameters are ε1 = 1.44, R1 = 0.24 μm, and μc = 0.3 eV. The structure can be simply analyzed by the modified Mie-Lorenz coefficients. The general design concepts are similar to their cylindrical counterparts, namely, dispersion engineering using the associated planar structure, as shown in the inset of the figure. Due to the excitation of TM surface plasmons, the normalized extinction cross-section is five times greater than the bare dielectric sphere. Moreover, similar to the cylindrical super-scatterers, by considering a small amount of loss for the graphene coating by assigning , the performance is considerably degraded .
3.1.4 Application in wide-band cloaking
By pattering graphene-based disks with various radii around a dielectric sphere, it is feasible to design a wide-band electromagnetic cloak at infrared frequencies. The geometry of this structure is illustrated in Figure 15. In order to analyze the proposed cloak by the modified Mie-Lorenz theory, the polarizability of the disks can be inserted in the equivalent conductivity method. The extracted equivalent surface conductivity can be used to tune the surface reactance of the sphere for the purpose of cloaking .
3.1.5 Application in multi-frequency cloaking
The other application that can be adapted to our proposed formulation of multilayered spherical structures is multi-frequency cloaking. As Figure 16 shows, by proper design, a single graphene coating can eliminate the dipole resonace in a single reconfigurable frequency. The radius of the sphere is R1 = 100 nm and its core permittivity is ε1 = 3. It can be concluded that double graphene shells can suppress the scattering in the dual frequencies since each graphene shell with different geometrical and optical properties can support localized surface plasmon resonances in a specific frequency. By further increase of the graphene shells, other frequency bands can be generated. Figure 16(b) shows the cloaking performance of a spherical particle with multiple graphene shells. The radii of the spheres are 107.5, 131.5, and 140 nm, and the corresponding chemical potentials are 900, 500, and 700 meV, respectively. The permittivity of the dielectric filler is 2.1 .
3.1.6 Application in electromagnetic absorption
As another example, a dielectric-metal core-shell spherical resonator (DMCSR) with the resonance frequency lying in the near-infrared spectrum is considered. In order to increase the optical absorption, the outer layer of the structure is covered with graphene. The localized surface plasmons of graphene are mainly excited in the far-infrared frequencies and in the near-infrared and visible range; it behaves like a dielectric. By hybridizing the graphene with a resonator, its optical absorption can be greatly enhanced. Figure 17 shows the performance of the structure for various core radii .
The provided examples are just a few instances for scattering analysis of graphene-based structures. Based on the derived formulas, other novel optoelectronic devices based on graphene plasmons can be proposed. Moreover, since assemblies of polarizable particles fabricated by graphene exhibit interesting properties such as enhanced absorption, negative permittivity, giant near-field enhancement, and large enhancements in the emission and the radiation of the dipole emitters [40, 41, 42, 43], the research can be extended to the multiple scattering theory.