Graphene coated subwavelength wires: A theorical investigation of the emission and radiation properties
This work analyzes the emission and radiation properties of a single optical emitter embedded in a graphene–coated subwavelength wire. We discuss the modifications of the spontaneous emission rate and the radiation efficiency as a function of the position and orientation of the dipole inside the wire. Our results show that these quantities can be enhanced by several orders of magnitude when the emission frequency coincides with one of the resonance frequencies of the graphene–coated wire. In particular, high–order plasmon resonances are excited when the emitter is moved from the wire center. The modifications by varying the orientation of the dipole in the near field distribution and in the far field intensities are shown.
Keywords:graphene, surface plasmons, quantum electrodynamics, plasmonics
Surface plasmon polariton results from the coherent coupling of photons to surface charge density oscillations [1, 2]. In natural plasmonic materials such as metals, constitutive parameters such as the conductivity and charge density are fixed, whereas they can be tuned in graphene electrically or by chemical doping , which has a dramatic effect on its optical properties. As the linear band structure of graphene causes the plasmon mass to depend on the Fermi–level position, electrically tunable surface plasmons can be supported by graphene from microwaves to the mid–infrared regimes . This has aroused new interest in studying graphene in the context of optical and plasmonic applications, including graphene quantum dots as a new generation of light–emitting devices , sensing , solar cells , to mention just a few.
Surface plasmons polaritons can be roughly divided into two categories: surface plasmon (SPs) propagating along waveguiding structures, such as an infinite flat graphene monolayer, and localized surface plasmons (LSPs) supported by spatially limited structures, such as scattering particles. Both kind of plasmon modes can be excited when a single optical emitter (such as quantum dots and single molecules) is placed next to one of the aforementioned structures. As a result of the high light confinement, an enhanced decay rate of the emitter into the plasmonic mode via the Purcell effect take place . In fact, it was recently showed that the interplay between an optical emitter and a graphene–coated sphere leads to an enhancement of the spontaneous emission into confined plasmonic modes . A variety of structures such as infinite graphene monolayers [10, 11, 12, 13], ribbons or nanometer sized disks  and double graphene waveguides [15, 16] have been the object of intensive research over the last few years due to the possibility to engineer surface plasmon mode density of states to control emission properties.
In this paper we consider a cylindrical dielectric core coated with a graphene layer and we investigate the role of the LSPs in modifying the emission and radiation rates of an emitter placed inside the graphene–coated cylinder. In this context, several works focused on the influence that the eigenmodes play on circular dielectric and metallic waveguides [17, 18, 19] or on single wall carbon nanotubes . Such systems offer a new tool for the interplay between light and matter at the nanometer scale [21, 22]. Other works have presented two dimensional calculations that extend such studies to examining the role of the shape of the waveguide [23, 24]. Two dimensional setting has important properties that allow to understand better the three dimensional case [17, 28], since qualitatively similar trends will hold in three dimensions . Experimental advances in the fabrication of micro–structures based on insulating materials such as SiO doped with molecular ions , the possibility of encapsulate single atom, molecules and compounds into graphene wires  and the fact that thanks to the van der Waals force, a graphene sheet can be tightly coated on a fiber surface , encourages an investigation about the emission and radiation characteristics of such graphene based systems.
The paper is organized as follows. First, in Section 2 we sketch an analytical method based on the separation of variables approach and obtain a solution for the electromagnetic field scattered by a graphene–coated wire when an oscillating line dipole source is located at an arbitrary position inside the wire cylinder. We derive analytical expressions for the power emitted and radiated by the source. In section 3 we present examples of emission and radiation decay rates corresponding to wires tightly coated with a graphene layer and compared the results with those obtained in the absence of the graphene coating. Finally, concluding remarks are provided in Section 4. The Gaussian system of units is used and an time–dependence is implicit throughout the paper, with as the angular frequency, as the time, and . The symbols Re and Im are used for denoting the real and imaginary parts of a complex quantity, respectively.
We consider a graphene coated cylinder with circular cross section (radius ) centered at , (figure 1). The wire substrate is characterized by the electric permittivity and the magnetic permeability . The coated wire is embedded in a transparent medium with electric permittivity and magnetic permeability . The graphene layer is considered as an infinitesimally thin, local and isotropic two–sided layer with frequency–dependent surface conductivity given by the Kubo formula [29, 30], which can be read as , with the intraband and interband contributions being
where is the chemical potential (controlled with the help of a gate voltage), the carriers scattering rate, the electron charge, the Boltzmann constant and the reduced Planck constant. The intraband contribution dominates for large doping and is a generalization of the Drude model for the case of arbitrary band structure, whereas the interband contribution dominates for large frequencies . A line dipole source (whose axis lie along the axis) with a dipole moment is placed inside the cylinder, at position (). The dipole is aligned at an angle with respect to the axis, as indicated in figure 1. The current density of the electric dipole is
In an unbounded medium, the dipole fields are obtained from the vector potential (refer to appendix for its derivation),
where is the non–zero component of the total magnetic field along the axis of the wire ( axis), is the transverse part of the operator, is the modulus of the photon wave vector in vacuum, is the angular frequency, is the vacuum speed of light, () is the smaller (larger) of and , and are the projection of vector on the and axis, respectively, and and are the nth Bessel and Hankel functions of the first kind, respectively. From Eqs. (4) and (5), we obtain the primary magnetic field emitted by the dipole,
for , and
When the dipole is located inside the coated wire cylinder, the scattered magnetic field along the axis of the wire, i.e, the component, denoted by (), are expanded as a series of cylindrical harmonics, one for the internal region (, superscript ) and another one for the external region (, superscript ),
where and are unknown complex coefficients. Taking into account that because of the graphene coating the tangential components of the magnetic field are no longer continuous across the boundary –as they were in the case of uncoated cylinders– the boundary conditions for our case can be expressed as
and , . Here, the prime denotes the derivative with respect to the argument. The denominator in Eqs. (13) and (14) are the same as that of the Mie scattering coefficients and may become very small in magnitude when the emission frequency coincides with one of the resonance frequencies [31, 32].
In the limit of , , and, as a consequence the amplitudes and are zeros except and , as expected from symmetry arguments. Once the amplitudes are determined, the scattered field, given by Eqs. (9) and (10) can be calculated at every point in the interior region (medium 1) and in the exterior region (medium 2). The total electromagnetic field allows us to calculate optical characteristics such as the power emitted and the power radiated by the dipole. The time average emitted power can be calculated from the integral of the normal component of the complex Poynting vector flux through an imaginary cylinder of length and radius that encloses the dipole (see Figure 1)
After some algebraic manipulation, the following results are found (see appendix)
Similarly, the time average radiative power emitted can be evaluated by calculating the complex Poynting vector flux through an imaginary cylinder of length and radius that encloses the graphene–coated cylinder (see Figure 1)
In the far–field region the calculation of the scattered fields given by Eq. (10) can be greatily simplified using the asymptotic expansion of the Hankel function for large argument. After some algebraic manipulation, we obtain
To characterize quantitatively the effect of the graphene coating, we define the normalized spontaneous emission rate as the ratio between the power emitted by the dipole, given by Eq. (25), and the power emitted by the same dipole embedded in an unbounded medium 1. In a similar way, the radiative efficiency is defined as the ratio between the power radiated by the dipole, given by Eq. (27), and the power emitted by the dipole in the unbounded medium 1.
We first briefly illustrate the emission and radiation characteristics in a dielectric wire without graphene coating (bare wire). Next, we consider a dielectric wire wrapped with a graphene coating. In all the examples the core (radius m) is made of a transparent material (, ) and is embedded in vacuum (). Rotational symmetry allows the emitter to be positioned on the axis ().
3.1 Uncoated dielectric wire
In Figure 2a we plot the frequency dependence of the normalized spontaneous emission (per unit length) for an emitter localized at , on the wire axis, and for m. The dipole moment is oriented in the radial direction (). Due to the dielectric wire is lossless, the curve in this figure also represents the radiation decay rate efficiency of the emitter. The peaks on the curve for correspond to the first two Mie dipolar resonances and are associated with the complex poles at the zeros of the denominators of the coefficients and in Eqs. (13) and (14) with .
The condition generates a series of complex roots indexed , whose real parts determine the resonant frequencies and the imaginary parts determine the quality factor of the resonances (). By using a Newton Raphson method, we have obtained m and m, whose real part is in very good agreement with the spectral position of the peaks. Consistent with causality, is strictly negative. The emission curve for m shows the peaks corresponding to resonances whose eigenfrequencies are zeroes of the denominator in Eqs. (13) and (14) with and with (see table 1). The values of the damping rates decrease with , resulting in high resonances as frequency is increased. This fact is confirmed by the curve for m in Figure 2a, since the width at half maximum of resonance peaks are decreasing with increasing frequency.
On the other hand, similar values between real parts of the dipolar and of the octupolar eigenfrequencies, suggest that these two resonances can be excited at the same emission frequency. This fact is showed in Figure 2a, where two peaks near m, one on the curve for and the other one on the curve for m, are marked with a dot. The insets show the spatial distribution of the magnetic field calculated at the emission frequency m for the case in which the emitter is localized at the wire center (left) and in case where the emitter is localized at m (right). This result, is a clear confirmation that different resonances (dipolar and octupolar in this case) can be excited by moving the emitter from the wire center.
In Figure 2b the normalized spontaneous emission curves are depicted for different orientation angles when the source is placed at m. To illustrate the effects of varying the orientation angle in the field near to the dielectric wire, in Figure 2c we have plotted the spatial distribution of the magnetic field for the octupolar resonance and for . We observe that these shapes are rotated and highly distorted with respect the shape calculated for showed in the inset of Figure 2a.
3.2 Graphene–coated wire
Having studied bare dielectric wires, we next explore the effects that a graphene coating has on the emission and the radiation spectrum of a dipole emitter localized inside the wire. In particular, the coating is expected to introduce LSPs mechanisms, which were absent in the bare wire, which are able to enhance the emission properties at resonance frequencies. The dashed and solid curves in Figure 3 display the frequency dependence of the spontaneous emission and the radiation efficiencies for the same wire considered in the previous example but now it is wrapped with a graphene coating. We used Kubo parameters eV, meV, K, emission frequencies in the range between m ( THz or wavelength m) and m ( THz or wavelength m), and for three different locations of the dipole, , m, and m.
Figure 3a shows that both and efficiencies are enhanced at a frequency near m (m) corresponding to the dipolar plasmon resonance of the graphene–coated wire. The correspondence between the spectral position of these peaks and the dipolar plasmon resonance can be clearly seen by calculating the complex root of the common denominator in equations (13) and (14). We have obtained m whose real part is in very good agreement with the spectral position of the peaks. In contrast to bare wires, where the emission decay rate reaches values close to unity at multipolar resonances, we obtain enhancement factors in both the emission and radiation decay rate efficiencies of around five orders of magnitude when the wire is wrapped with a graphene coating. Furthermore, we find that high–order plasmon resonances are excited when the emitter is moved from the wire center, as can be seen in Figure 3b where the curves of and have been plotted for m. In addition to the maximum enhancement in both the emission and the radiation spectra at a frequency near m ( or dipolar resonance), a local maximum enhancement near m corresponding to the quadrupolar plasmon resonance () appear. Other local maxima, particularly noticeable in the emission spectra, occur at frequencies near m ( or hexapolar resonance) and m ( or octupolar resonance). By calculating the complex root of the common denominator in equations (13) and (14), we have obtained m, m and m for which their real parts are in very good agreement with the spectral position of the peaks. When is increased, the enhancement factors in the spontaneous emission curve at frequencies corresponding to high–order resonances becomes more and more pronounced. This fact can be viewed in figure 3c in which the emission efficiencies for reach similar values ().
It is worth noting that the width at half maximum of the resonance peaks in Figure 3c, i.e., the quality factor of the resonance , almost not shows any significant dependence on the position of the peak. This agree well with the fact that the damping rate take the value m for large enough values . Although the enhancement factors in the radiation curve are less noticeable than those in the emission curve, these peaks are increased when is increased, as can be seen by comparing the radiation curves in Figures 3b and 3c. Moreover, we observe that a peak at a frequency near m (in which hexapolar order resonance occur), absent in the curve corresponding to m, becomes just visible.
In Figure 4 we plot the spatial distribution of the near magnetic field for the wire considered in Figure 3 and for the first four plasmon resonances. We have considered the particular case of the emitter localized at m as in Figure 3c.
The emission frequencies are m (Figure 4a), m (Figure 4b), m (Figure 4c), and m (Figure 4d), i.e., the values corresponding to the real part of the complex poles of the multipole coefficients and () and for which the strongest maxima in the emission efficiency occur. At these frequencies, the near field distributions follow the typical dipolar, quadrupolar, hexapolar and octupolar patterns. Contrary to what happens in bare wires, the presence of a surface current density induced on the graphene coating and the boundary condition in Eq. (12) lead to a discontinuity at m as can be observed in Fig. 4. For example, the map of the magnetic field showed in Figure 4c is characterized by six zones, three of them are red inside the cylinder and are yellow outside the cylinder, pointing out that the magnetic field changes from negative values to positive values when passing from inside to outside of the cylinder. In the other three zones, yellow inside the cylinder and red outside the cylinder, the magnetic field changes from positive to negative values when passing from inside to outside of the cylinder.
As well as the near field maps at resonance frequencies are closely related with the amplitude and location of the maxima in the emission spectra, the far field maps are related with the peak characteristics in the radiation spectra. The inset in Figure 4a shows that, for the emission frequency m (Re), the far field intensities, namely which is the fraction of the emitted power that is scattered into the angular region about the scattering direction , are dominated by an dipole pattern. The integral of over all angles, i.e., the radiative decay rate, is found to be times larger than that of the same dipole emitter in the absence of the graphene–coated wire, as can be seen in Figure 3c. When the emission frequency is m (Re), the radiation efficiency showed in Figure 3c exhibits a strong enhancement (an enhancement factor near three orders of magnitude greater than that of the same emitter in the absence of the graphene–coated wire) and the far field intensities are dominated by an electromagnetic quadrupole pattern, as indicated in the inset in figure 4b. Similarly, when the emission frequency is m (Re), the radiation efficiency plotted in Figure 3c exhibits a very weak peak, and as a result, a small hexapolar pattern overlaps the dipole pattern in the far–field intensities, as indicated in the inset in figure 4c. On the other hand, no peak is observed in the radiation efficiency curve when the emission frequency is m (Re). At this frequency, the far–field intensities are dominated by an electromagnetic dipole, as indicated in the inset in Figure 4d.
Contrary to the case of bare wires, where both the emission and radiation efficiencies are highly modified when the dipole is rotated, in the graphene coated dielectric wire we have observed that neither nor efficiencies show a noticeable dependence on the orientation angle . However, the shape of the near field distribution and of the far field intensities strongly depend on . This fact can be seen in Figure 5, where we have plotted the near field and the power patterns for four dipole orientations: , at a specifically frequency corresponding to the quadrupolar plasmon resonance.
When the dipole is rotated, both the near field map and the far field intensities are rotated with respect the case. Unlike the case of bare wires, where the near field distortion due to the rotation of the source is remarkable (see Figure 2c), Figure 5 shows that when the dielectric wire is wrapped with a graphene coating only a rotation of the near field is appreciable, being its distortion negligible. On the other hand, as is increased, the shape of the far field intensities are more and more distorted until arriving at the shape shown in Figure 5d corresponding to . Same behavior has been observed [not shown in Figure 5] at the other plasmon resonance frequencies.
To examine the effects of the emitter localization on plasmon resonances, Figure 6 shows the emission and the radiation decay rates as a function of the distance to the wire center. The emission frequency is chosen such that the emission is resonant at the orders: dipolar (a), quadrupolar (b), hexapolar (c) and octupolar (d). At distance m, the emission and the radiation decay rates for the dipole resonance are both times larger than in the absence of the graphene–coated cylinder. When is increased from this value, both efficiency curves decrease to reach their lowest value at m as its shown in Figure 6a. In contrast, for the quadrupole resonance, Figure 6b shows that both efficiency curves take a value close to unity at m and that these curves increase with the increasing of . The same behavior is observed in Figure 6c for the hexapole resonance, although the maximum value in the radiative decay efficiency is much lower than the maximum value in the quadrupolar resonance curve. For the octupole resonance, the value of the radiative decay efficiency does not shows any significant dependence on the position of the emitter, as can be seen in Figure 6d.
In conclusion, we have presented an analytical classical method based on the separation of variables approach to find the emission and radiation characteristics of an optical emitter inside a circular cross–section wire coated with a graphene layer. The method has been used to investigate how the location of the emitter, as well as the orientation of its dipole moment, affect the power emitted by the source. We found that the interplay between the optical emitter and the LSPs in the graphene coating strongly influences the spontaneous decay rate as well as the radiation characteristics. We have compared the results with those obtained for the same dielectric wire but without a graphene layer (bare wire), in which the decay rates are enhanced at resonances associated with cavity modes.
As might be expected from symmetry arguments, for the optical emitter placed on the wire center, we found that both the emission and radiation decay rates are enhanced just at dipolar resonance frequency. Instead, an enhancement at frequency of high–order resonances appear when the emitter is moved from the axis of the wire. The correspondence between the position of the spectral peaks and the multipolar resonances (cavity or plasmonic) has been clearly shown by calculating the complex poles of the coefficients of the multipole expansion of the electromagnetic field. In case of bare wires, cavity effects are intensified as frequency increases, i.e., where high resonances take place. On the contrary, for graphene coated wire, the resonance damping rate take the quasistatic value for large enough values. As a consequence, the resonance factor shows any dependence on the position of the spectral peaks. Moreover, whereas the emission and radiation decay rates are highly modified with the orientation of the dipole emitter in case of bare wires, we have found that neither emission nor radiation decay rates show a noticeable dependence on the orientation angle in case of graphene coated wires.
The spatial distribution of the electromagnetic field near the wire for different resonance frequencies and positions of the emitter have been investigated. The multipolar order revealed by the topology of the near field agrees well with the multipolar order revealed by the spectral position of the emission decay rate peak. We have shown the modifications in the near field distribution by varying the orientation of the dipole.
The author acknowledge the financial support of Consejo Nacional de Investigaciones Científicas y Técnicas, (CONICET, PIP 451).
5 Appendix I
In an unbounded medium (constitutve parameters , ) the vector potential satisfy the inhomogeneous Helmholtz equation 
where and denotes the position of the line dipole source, with , and . Since
we can write
Taking into account the finiteness at the origin and the radially outgoing property, we obtain
where () is the smaller (larger) of and . The coefficient is determined by the condition
By using where denotes the Wronskian, we obtain . Finally,
6 Appendix II
Here, we derive an analytical expression for the first term in Eq. (24), i.e., the contribution of the primary field emitted by the dipole. It is given by
where we have used , and .
It is worth noting that the same result can be obtained by using the Poynting theorem ,
where the field is the primary electric dipole field evaluated at the dipole position .
-  H. Raether 1988 Surface Plasmons On Smooth And Rough Surfaces And On Gratings (Berlin: Springer–Verlag)
-  Maier S A 2007 Plasmonics: Fundamentals and Applications (New York: Springer)
-  Fengnian Xia, Comment on Graphene versus metal plasmons, Nature Photonics 7, 420 (2013).
-  Jablan J., Soljacic M., Buljan H. 2013 Plasmons in graphene: fundamental properties and potential applications Proc. IEEE 101 1689–1704
-  S.H. Song, S.H. Jang, J. Chung, S.H. Jin, B.H. Kim, S.H. Hur, S. Yoo, Y.-H. Cho, S. Jeon Highly efficient light–emitting diode of graphene quantum dots fabricated from graphite intercalation compounds Adv. Opt. Mater., 2 (2014), pp. 1016–1023
-  María C. Dalfovo, Gabriela I. Lacconi, Moónica Moreno, Marta C. Yappert, Gamini U. Sumanasekera, Roberto C. Salvarezza, and Francisco J. Ibanez, Synergy between Graphene and Au Nanoparticles (Heterojunction) towards Quenching, Improving Raman Signal, and UV Light Sensing, Acs Appl. Mater. Interfaces, 6 (2014), pp 6384–8391
-  Francesco Bonaccorso, Luigi Colombo, Guihua Yu, Meryl Stoller, Valentina Tozzini, Andrea C. Ferrari, Rodney S. Ruoff, Vittorio Pellegrini, Graphene, related two–dimensional crystals, and hybrid systems for energy conversion and storage, Science 347, (6217) (2015).
-  Tame M S, McEnery K R, Ozdemir K, Lee J, Maier S A and Kim M S 2013 Quantum plasmonics Nature Physics 9 329–40
-  Christensen T, Jauho A-P, Wubs M and Mortensen N 2015 Localized plasmons in graphene–coated nanospheres Phys. Rev. B 91 125414
-  Hanson G W, Forati E, Linz W, Yakovlev A B 2012 Excitation of terahertz surface plasmons on graphene surfaces by an elementary dipole and quantum emitter: strong electrodynamic effect of dielectric support Phys. Rev. B 86 235440
-  Huidobro P A, Nikitin A Y, González–Ballestero C. , Martín–Moreno L, and García–Vidal F J 2012 Superradiance mediated by graphene surface plasmons Phys. Rev. B 85 155438
-  Kort–Kamp W J M, Amorim B, Bastos G, Pinheiro F A, Rosa F S S, Peres N M R, and Farina C 2015 Active magneto–optical control of spontaneous emission in graphene Phys. Rev. B 92 205415
-  Karanikolas V D, Marocico C A, and Bradley A L 2015 Effect of surface-plasmon polaritons on spontaneous emission and intermolecular energy–transfer rates in multilayered geometries Phys. Rev. B 91 125422
-  Koppens F H L, Chang D E and García de Abajo F J 2011 Graphene plasmonics: A platform for strong light–matter interactions Nano Lett. 11 3370–77
-  Zhang L, Fu X, Zhang M, and Yang J 2013 Spontaneous emission in paired graphene plasmonic waveguide structures Optic Epress 21 7897–907
-  Mauro Cuevas, Surface plasmon enhancement of spontaneous emission in graphene waveguides, Journal of Optics 18, 105003 (2016)
-  D. P. Fussell, R. C. McPhedran, and C. Martijn de Sterke, Decay rate and level shift in a circular dielectric waveguide, Phys. Rev. A 71, 013815 (2005).
-  Julien Barthes, Alexandre Bouhelier, Alain Dereux and Gerard Colas des Francs, Coupling of a dipolar emitter into one–dimensional surface plasmon, Scientific Reports 3, 2734, (2013).
-  Vasilios Karanikolas, Cristian A. Marocico, and A. Louise Bradley, Spontaneous emission and energy transfer rates near a coated metallic cylinder, Phys. Rev. A 89, 063817 (2014).
-  Luis Martín–Moreno, F. Javier García de Abajo, and Francisco J. García–Vidal, Ultraefficient Coupling of a Quantum Emitter to the Tunable Guided Plasmons of a Carbon Nanotube, Phys. Rev. Lett. 115, 173601 (2015).
-  David Dzsotjan, Anders S. Sorensen, and Michael Fleischhauer, Quantum emitters coupled to surface plasmons of a nanowire: A Green’s function approach, Phys. Rev. B 82, 075427 (2010).
-  A. V. Akimov, A. Mukherjee1, C. L. Yu, D. E. Chang, A. S. Zibrov, P. R. Hemmer, H. Park and M. D. Lukin, Generation of single optical plasmons in metallic nanowires coupled to quantum dots, Nature 450, 402–406 (2007).
-  L. Rogobete, C. Henkel, Spontaneous emission in a subwavelength environment characterized by boundary integral equations, Phys. Rev. 70, 063815 (2004).
-  Lavinia Rogobete, Hannes Schniepp, Vahid Sandoghdar, Carsten Henkel, Spontaneous emission in nanoscopic dielectric particles, Optics Letters 28, (2003).
-  S P Feofilov, A B Kulinkin, D A Eurov, D A Kurdyukov and V G Golubev, Fluorescence spectroscopy study of mesoporous SiO2 particles containing Gd2O3:Eu3Ã, Materials Research Express 1 (2014) 025019
-  G. H. Jeong, A. A. Farajian, R. Hatakeyama, T. Hirata, T. Yaguchi, K. Tohji, H. Mizuseki, and Y. Kawazoe, Cesium encapsulation in single-walled carbon nanotubes via plasma ion irradiation: Application to junction formation and ab initio investigation, Phys. Rev. B 68, 075410 (2003).
-  He X, Zhang X, Zhang H and Xu M 2014 Graphene covered on microfiber exhibiting polarization and polarizationdependent saturable absorption IEEE J. Sel. Top. Quantum Electron. 20 4500107
-  Vahid Siahpoush, Thomas Sondergaard, and Jesper Jung, Green’s function approach to investigate the excitation of surface plasmon polaritons in a nanometer-thin metal film, Phys. Rev. B 85, 075305 (2012).
-  Falkovsky F A 2008 Optical properties of graphene and IV–VI semiconductors. Phys Usp 51 887–97
-  Milkhailov S A and Siegler K 2007 New electromagnetic mode in graphene Phys Rev Lett 99 016803
-  Riso M, Cuevas M and Depine R A 2015 Tunable plasmonic enhancement of light scattering and absorption in graphenecoated subwavelength wires Journal of Optics 17 075001
-  M Cuevas , M Riso, and R A Depine 2016 Complex frequencies and field distributions of localized surface plasmon modes in graphene–coated subwavelength wires Journal of Quantitative Spectroscopy and Radiative Transfer 173 26–33
-  Novotny L, and Hecht B Principles of Nano–Optics; Cambridge University Press: New York, 2006.
-  Abramowitz M. and Stegun I. A., Handbook of Mathematical Functions (New York: Dover) (1965)