Complete Optical Absorption in Periodically Patterned Graphene

Phys. Rev. Lett. 108, 047401
Optical response of doped graphene nanodisks. (a) Characteristic extinction cross section of a nanodisk (60 nm in diameter and Fermi energy EF=0.4eV), dominated by a pronounced dipolar plasmon and strong near-field enhancement. Solid curve: Full electrodynamic calculation. Dashed curve: Fit to a Lorentzian polarizability. Inset: Near-field intensity normalized to incident intensity for plane-wave illumination with polarization parallel to the disk (the disk is shown in perspective relative to the plane of the plot). (b),(c) Size dependence of the maximum extinction cross section and the polarizability parameters [see Eq. (1)].Optical response of doped graphene nanodisks. (a) Characteristic extinction cross section of a nanodisk (60 nm in diameter and Fermi energy EF=0.4eV), dominated by a pronounced dipolar plasmon and strong near-field enhancement. Solid curve: Full el... Show more

Light absorption plays a central role in optical detectors and photovoltaics. Inspired by nature, two different routes have been investigated to achieve perfect absorption. (i) A first one consists in relying on diffusion in disordered lossy surfaces (e.g., black silver and carbon). Engineered materials have been synthesized following this solution to produce extraordinary broadband light absorption (e.g., dense arrays of carbon nanotubes [1]). (ii) A second approach consists in using ordered periodic structures, as found in some nocturnal insects, where they produce the moth eye effect [2]. This alternative has been pioneered by experimental and theoretical work showing total light absorption (TLA) in the visible using metallic gratings [3–5]. In this context, the Salisbury screen [6,7], consisting of a thin absorbing layer placed above a reflecting surface, has been known to produce TLA, and it can be integrated in thin structures using magnetic-mirror metamaterials [8]. Similar phenomena have been reported at infrared (IR) [9–11] and microwave [12,13] frequencies, including omnidirectional TLA [14], which has been realized by using periodic surfaces supporting localized plasmon excitations.

The availability of high-quality graphene as a stable material with extraordinary (opto)electronic properties [15–17] makes a compelling case for exploring its ability to harvest light for potential application to optoelectronics, with the advantage of being optically tunable via electrostatic doping [18]. However, a single sheet of homogeneous graphene is poorly absorbing [19] (about 2.3% absorption), so the challenge is to transform it into a perfect absorber, for which we can rely on its power to host extremely confined plasmons [20,21].

In this Letter, we show that a single sheet of doped graphene, patterned into a periodic array of nanodisks, exhibits 100% light absorption. We first discuss the extinction cross section of graphene disks (i.e., the sum of light absorption and elastic scattering), which can exceed by over an order of magnitude their geometrical area [21]. These disks therefore belong to the class of absorbing particles that can be arranged in periodic planar arrays such that their cross section exceeds the area of the unit cell. We further assume here that all diffracted beams are evanescent. Under these conditions, we derive a universal maximum absorption for any thin layer separating two different media. In particular, TLA is predicted when the transmission channel is suppressed, for instance, under total internal reflection (TIR) by illumination from a prism substrate, and also when the particles are supported on a dielectric film on top of a thick metal [6,7]. We show that a patterned graphene sheet displays these effects. Full absorption by an atomically thin carbon film is thus possible by using currently available technology.

Optical response of graphene nanodisks.—Doped graphene has been predicted to support strongly confined, long-lived plasmons [20], which produce sharp IR resonances in nanodisks tens of times smaller than the light wavelength [21]. Figure 1(a) shows an example simulated with the boundary element method [22]. We find that the polarizability of these disks exhibits a resonance that is well described by a Lorentzian line shape as [23]

Absorption by a layer of doped graphene nanodisks. (a) Scheme of the geometry under consideration. (b),(c) Universal contours showing a maximum limit to the absorbance of a thin material layer sandwiched between two different media for s and p polarization as a function of incidence angle and dielectric contrast between both media. Light is incident from the ϵ1 medium. (d),(e) Realization of maximum absorption with arrays of graphene disks [see (a)] of diameter D=60nm, EF=0.4eV, and different lattice spacings and angles of incidence (see labels). The lower-index medium has ϵ=1 in all cases. The numerical labels in (d) and (e) correspond to the positions of the symbols in (b) and (c).Absorption by a layer of doped graphene nanodisks. (a) Scheme of the geometry under consideration. (b),(c) Universal contours showing a maximum limit to the absorbance of a thin material layer sandwiched between two different media for s and p polari... Show more

where ωp is the plasmon frequency, κ is its decay rate, and κr is the radiative contribution to κ. These fitting parameters emerge from our full electromagnetic simulations of the graphene, and they bear an intrinsic dependence on the geometrical and electronic properties of the carbon structure. The extinction cross section σext=(4πω/c)Im{α} obtained from α [Fig. 1(a), broken curve] is in good agreement with rigorous numerical solutions of Maxwell’s equations (solid curve). We represent the parameters ωp, κ, and κr in Figs. 1(b) and 1(c) for different doping levels (quantified by the Fermi energy EF) as a function of disk size. Here, we compute the conductivity of graphene in the local-RPA limit [24] with an intrinsic relaxation time τ=μEF/evF2, where vFc/300 is the Fermi velocity and μ=10000cm2/Vs is the measured dc mobility [15] (e.g., τ10-13s for EF=0.1eV). Interestingly, the maximum extinction cross section σmaxext(3λ2/2π)κr/κ at the resonant wavelength λ can exceed by over an order of magnitude the geometrical area, although the disks have κκr, so that they are far from perfect two-level systems (such as an atom without inelastic decay channels), for which σmaxext0.48λ2.

Universal limit to absorption by a thin layer.—In a symmetric environment, the reflection coefficient r of a thin material layer (much thinner than λ) determines the transmission coefficient 1±r (upper and lower signs are for s- and p-polarized light, respectively, as required to enforce the symmetry of the scattered field with respect to the layer plane [25]; we use this convention from now on). Therefore, the maximum absorption ( 1-|r|2-|1±r|2) is limited to 50% (corresponding to r=1/2). This value can be increased by considering asymmetric environments, such as the patterned graphene sheet depicted in Fig. 2(a).

A simple analysis leads to the maximum possible absorption produced by a thin, periodically structured material layer that is infinitesimally close to the interface between two media of different permittivities, ϵ1 and ϵ2. The reflection and transmission coefficients of this structure must read

where r0 and t0 are the Fresnel coefficients of the bare ϵ1|ϵ2 interface, while η is the self-consistent amplitude of the wave scattered by the layer in response to both the external field and the multiple reflections within the layer-interface cavity. The scattered field is obviously symmetric for a layer of thickness λ, and it contributes to Eq. (2) with direct ( η) and interface-reflected ( ±r0η) components, as well as the corresponding transmission ( ±t0η). The absorbance A=1-|R|2-Re{f}|T|2, obtained from the integral of the Poynting vector over planes far from the structure, is thus an analytical function of the complex variable η, and its maximum admits the closed-form expressions [25]

Angular dependence of the peak absorbance for s (a) and p (b) polarization in arrays of graphene nanodisks [see Fig. 2a]. Each curve represents a fixed sample with parameters as shown by labels.Angular dependence of the peak absorbance for s (a) and p (b) polarization in arrays of graphene nanodisks [see Fig. 2a]. Each curve represents a fixed sample with parameters as shown by labels. Show more


for s and p polarizations, respectively, where f=(ϵ2/ϵ1-sin2θ)1/2/cosθ and θ is the angle of incidence from medium 1. The dependence of the universal maximum absorption on θ and ϵ2/ϵ1 predicted by Eqs. (3) and (4) is shown in Figs. 2(b) and 2(c). Interestingly, it becomes 100% under TIR from the high-index medium (i.e., for ϵ1sin2θ>ϵ2, leading to Re{f}=0), indicating that it is possible to simply suppress the reflection when transmission is already prevented by TIR.

Absorption in particle and disk arrays.—The universal maximum is just an upper limit to the absorption, but its derivation does not provide any clue on how to structure the material to reach such a limit, or whether it is at all possible to obtain it. However, a systematic inspection of graphene nanodisk arrays reveals that the universal limit is indeed reached by engineering the disk size, doping, and spacing. We show several examples in Figs. 2(d) and 2(e), computed from a layer-multiple-scattering method [26], in which the disks are described through their multipolar polarizability obtained from the boundary element method (solid curves; see Ref. [25]).

We can gain further insight into thin-layer absorption by describing each graphene nanodisk (or any other particle) through the polarizability [Eq. (1)]. For arrays of small period a compared to the wavelength, the reflection coefficient of the layer reduces to [25,27]

Total absorption with graphene on a coated metallic substrate. (a) Normal-incidence absorbance by graphene disk arrays supported on a dielectric-coated gold surface (see the inset) for various values of the dielectric film thickness (ϵdielectric=2.1, ϵAu from Ref. 32). (b) Peak absorbance as a function of dielectric thickness d and lattice spacing a for the same disks as in (a). (c),(d) Incidence-angle and photon-energy dependence of the absorbance under the conditions of (a) for d=1.12μm.Total absorption with graphene on a coated metallic substrate. (a) Normal-incidence absorbance by graphene disk arrays supported on a dielectric-coated gold surface (see the inset) for various values of the dielectric film thickness (ϵdielectric=2.1,... Show more

where G is a lattice sum whose imaginary part is exactly given by Im{G}=S-2(ω/c)3/3, whereas its real part reduces to Re{G}g/a3 for λa, with g=5.52 ( g=4.52) for hexagonal (square) arrays. Here, S=2πω/cAcosθ ( S=2πωcosθ/cA) for s-polarized ( p-polarized) light, and A is the unit-cell area. Equation (5), together with Eq. (1), permits obtaining the transmission and reflection coefficients of the combined array-interface system analytically [they are simply given by Eq. (2) with η=r(1±r0)/(1-r0r) [25]], leading to the dashed curves of Figs. 2(d) and 2(e), in excellent agreement with full simulations.

In the symmetric configuration ( ϵ1=ϵ2), by combining Eqs. (1) and (5) and noticing that κκr, the maximum absorption (i.e., 50%, corresponding to r=1/2) is obtained at a frequency ωωp-3gκr/4(ωpa/c)3 (slightly redshifted with respect to the single-disk resonance ωp), under the condition ζ=1/2, where we define the parameter

This condition can be fulfilled by adjusting the spacing between disks, provided σmaxext is sufficiently large compared to the disk area. For asymmetric environments ( ϵ1ϵ2), the analysis becomes more involved, but similar redshifts are obtained for large enough σmaxext.

The robustness of TLA under TIR is explored in Fig. 3, where the peak absorbance is represented as a function of incidence angle for various values of the lattice spacing a, using the same disks in all cases [see parameters in Fig. 3(b)]. Nearly TLA is observed within the 20°–45° range of incidence angles under critical-coupling conditions [14].

Omnidirectional total light absorption.—The above prediction of TLA under TIR can be generalized to systems in which the transmission channel is suppressed [6,7]. In particular, we consider a graphene array placed above a metallic substrate. The graphene plasmons disappear close to a metal, but this problem is solved by introducing an intermediate dielectric coating layer [see the inset in Fig. 4(a)]. In this configuration, the metal can naturally act as a back-gate element to dope the graphene. It is easy to see that reflection is canceled when the substrate (i.e., metal plus coating) reflection coefficient is r0=-r/(1±2r) [25]. The coating-layer thickness d can be adjusted to match the required phase in this expression, and the condition is simply stated as 1/|r0|=|2±1/r|1, where the equality applies to nonabsorbing substrates (e.g., a dielectric under TIR or noble metals in the IR). Figure 4(b) shows TLA under this scheme for properly adjusted values of a and d. Obviously, the phase-matching condition is periodically satisfied along the d axis. Interestingly, the incidence-angle and polarization dependence of the absorption is very weak [Figs. 4(c) and 4(d)], pointing at the possibility of omnidirectional absorption. To understand this, we combine Eqs. (1) and (5) together with the relation |2±1/r|=1 (i.e., the condition of no reflection for |r0|=1), from which we find two frequencies for TLA: ω=ωp-3gκr/4(ωpa/c)3±(κ/2)1/ζ-1, under the condition ζ<1; only the right-most term depends on polarization and angle of incidence (through ζ), so total absorption can be made omnidirectional when the second term in the resonant ω dominates [i.e., for (ωa/c)3κr/κ], provided σmaxext is of the order of A.

Concluding remarks.—Planar textured materials can produce large light absorption, which reaches 100% when the transmission channel is suppressed. We predict that this effect takes place by using a single sheet of doped patterned graphene either under TIR or when the carbon layer is deposited on dielectric-coated metal. The absorbed energy is generally dissipated into phonons and electron-hole ( e-h) pairs, which eventually thermalize the structure to be either reemitted as thermal radiation or dissipated through heat diffusion. An interesting question is posed by the energy balance between absorption and emission in these resonant structures if they are made thermally isolated, with a view to converting incident light of a certain frequency into emission at another resonant frequency. In a different direction, thermally heated samples can be useful as narrow-band IR sources in virtue of Kirchhoff’s law, whereby emission at selected wavelengths determined via electrostatic doping takes place by heating the structure.

We also note that full absorption in graphene opens interesting possibilities to boost the efficiency of terahertz and IR detection (e.g., via e-h separation [28], thermoelectric reading [29,30], etc.). This type of detector can inherit the electrical tunability of graphene, to be used for direct, efficient spectral analysis of IR light with relative spectral resolution κ. Specifically, we think of a battery of perfectly absorbing graphene arrays operating at different wavelengths and covering a wide spectral region, so that optical or thermoelectric reading of each array provides the intensity of each spectral channel. On a more speculative direction, e-h separation in graphene [29] may potentially be employed to convert light into electrical power. The tunability of graphene can be used to absorb the full IR spectrum for application to solar cells by first spectrally separating the incident light and subsequently absorbing it at different resonant graphene structures.

It should be mentioned that the configuration of doped nanodisks presents a challenge. Besides the possibility of having narrow graphene bridges electrically connecting them to achieve electrostatic doping, which might damage the optical performance, there are other alternatives, such as decorating the structure with bridges made of a transparent, conductive material. Chemical doping is another option, which can be useful for sensing changes in the surrounding fluid. Besides, further work on ribbon arrays with similar structure is in progress, also leading to TLA, and allowing conductive contacts to be placed far from the bulk of the structure for long ribbons. Finally, we predict perfect absorption in planar particle arrays under the condition that the absorption cross section of an individual particle is comparable to the area of the unit cell. This condition is actually fulfilled by many nonideal two-level atoms or molecules. Our results could be tested with optically trapped atomic arrays, although graphene nanodisks are more convenient because they provide a practical implementation that is stable under ambient conditions. Noble-metal nanoparticle arrays can also exhibit similar effects in the near IR, with their spacing controlled through dielectric coating [31].


  1. K. Mizuno, J. Ishii, H. Kishida, Y. Hayamizu, S. Yasuda, D. N. Futaba, M. Yumura, and K. Hata, Proc. Natl. Acad. Sci. U.S.A. 106, 6044 (2009).
  2. P. B. Clapham and M. C. Hutley, Nature (London) 244, 281 (1973).
  3. M. C. Hutley and D. Maystre, Opt. Commun. 19, 431 (1976).
  4. M. Nevière, D. Maystre, R. MacPhedran, G. Derrick, and M. Hutley, in Proceedings of the ICO-11 Conference, Madrid, Spain, 1978 (unpublished), pp. 609–612.
  5. E. Popov, D. Maystre, R. C. McPhedran, M. Nevière, M. C. Hutley, and G. H. Derrick, Opt. Express 16, 6146 (2008).
  6. W. W. Salisbury, U.S. Patent No. 2 599 944 (1952).
  7. R. L. Fante and M. T. McCormack, IEEE Trans. Antennas Propag. 36, 1443 (1988).
  8. N. Engheta, in Digest of the 2002 IEEE AP-S International Symposium, San Antonio, TX (IEEE, New York, 2002), Vol. 2, pp. 392–395.
  9. E. Popov, L. Tsonev, and D. Maystre, Appl. Opt. 33, 5214 (1994).
  10. N. Liu, M. Mesch, T. Weiss, M. Hentschel, and H. Giessen, Nano Lett. 10, 2342 (2010).
  11. C. Hägglund, S. P. Apell, and B. Kasemo, Nano Lett. 10, 3135 (2010); ; C. Hägglund and S. P. Apell, Opt. Express 18, A343 (2010).
  12. Y. P. Bliokh, J. Felsteiner, and Y. Z. Slutsker, Phys. Rev. Lett. 95, 165003 (2005).
  13. N. I. Landy, S. Sajuyigbe, J. J. Mock, D. R. Smith, and W. J. Padilla, Phys. Rev. Lett. 100, 207402 (2008).
  14. T. V. Teperik, F. J. García de Abajo, A. G. Borisov, M. Abdelsalam, P. N. Bartlett, Y. Sugawara, and J. J. Baumberg, Nature Photon. 2, 299 (2008).
  15. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, Science 306, 666 (2004).
  16. K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and A. A. Firsov, Nature (London) 438, 197 (2005).
  17. F. Bonaccorso, Z. Sun, T. Hasan, and A. C. Ferrari, Nature Photon. 4, 611 (2010).
  18. C. F. Chen et al., Nature (London) 471, 617 (2011).
  19. K. F. Mak, M. Y. Sfeir, Y. Wu, C. H. Lui, J. A. Misewich, and T. F. Heinz, Phys. Rev. Lett. 101, 196405 (2008).
  20. M. Jablan, H. Buljan, and M. Soljačić, Phys. Rev. B 80, 245435 (2009).
  21. F. H. L. Koppens, D. E. Chang, and F. J. García de Abajo, Nano Lett. 11, 3370 (2011).
  22. F. J. García de Abajo and A. Howie, Phys. Rev. B 65, 115418 (2002).
  23. D. V. van Coevorden, R. Sprik, A. Tip, and A. Lagendijk, Phys. Rev. Lett. 77, 2412 (1996).
  24. B. Wunsch, T. Stauber, F. Sols, and F. Guinea, New J. Phys. 8, 318 (2006).
  25. See Supplemental Material at for more details on the theory.
  26. N. Stefanou, V. Yannopapas, and A. Modinos, Comput. Phys. Commun. 113, 49 (1998); ; 132, 189 (2000).
  27. F. J. García de Abajo, Rev. Mod. Phys. 79, 1267 (2007).
  28. F. N. Xia, T. Mueller, Y. M. Lin, A. Valdes-Garcia, and P. Avouris, Nature Nanotech. 4, 839 (2009).
  29. M. C. Lemme et al., Nano Lett. 11, 4134 (2011).
  30. J. C. W. Song et al., Nano Lett. 11, 4688 (2011).
  31. T. Ung, L. M. Liz-Marzán, and P. Mulvaney, J. Phys. Chem. B 105, 3441 (2001).
  32. P. B. Johnson and R. W. Christy, Phys. Rev. B 6, 4370 (1972).

About the Authors

Image of Sukosin Thongrattanasiri
Image of Frank H. L. Koppens
Image of F. Javier García de Abajo

Javier García de Abajo obtained his Ph.D. in condensed matter theory from the University of the Basque Country (Spain) in 1993. He is currently a Research Professor at the Instituto de Óptica (CSIC) in Madrid. He is a Fellow of the APS and the OSA. His main interests are plasmons in nanoparticles, the optical response of patterned surfaces, and the interaction of swift electrons with nanostructured materials.

Related Articles

Q&amp;A: Bringing Space Down to Earth

Q&A: Bringing Space Down to Earth

Chiara La Tessa hopes to contribute to cancer treatments and human-missioned trips to Mars by exploring the effects of high-energy particles on living things. Read More »

Viewpoint: Synchronized Cell Motion without Fluid Interactions
Biological Physics

Viewpoint: Synchronized Cell Motion without Fluid Interactions

A study of the breaststroke-motion in swimming algae shows that hydrodynamic forces cannot explain the synchronization between the microorganism’s two “arms.” Read More »

Focus: Water Molecules Need Help to Evaporate
Chemical Physics

Focus: Water Molecules Need Help to Evaporate

Each time a liquid water molecule enters the vapor phase, a coordinated dance of several molecules is involved, according to simulations. Read More »

More Articles