Abstract
Graphene absorbs light in accordance with the fine-structure constant α ≃ 1/137. The α value is universal in the sense that it is unrelated to material parameters but solely related to the elementary charge, which is a fundamental constant governing the coupling of light and matter. A new universality governed only by α has yet to be discovered in graphene optics. However, since α is a dimensionless quantity, it can be anticipated that the reciprocal of α appears as a characteristic number of layers in the light absorptance of an N-layer graphene. Here, we report that this number is 2/πα. We show simultaneously that for light in the infrared to visible range, an N-layer graphene with N above 1500 may be regarded as graphite. This also enables us to obtain the simplest expression for the optical constants of graphite, which is written in terms of α and interlayer distance only.
Similar content being viewed by others
Introduction
By a mechanical cleavage method1, the number of layers in graphite is decreased from a huge number to unity, and graphene can be simply defined as a single layer of graphite. But, above how many layers of graphene can we regard it as if “graphite”? This is a natural question for many readers. We answer the question by calculating the optical properties of an N-layer graphene as a function of the number of layers N. In the process of solving this problem, we noticed that various interesting aspects lie behind it. They include not only scientific problems but also important problems related to applications. For example, what is the optimum N for absorbing light most efficiently? Interestingly, the answer is related to a fundamental physics constant, namely the fine-structure constant α ≃ 1/137.
The exfoliation of a single-layer graphene from graphite provided a great opportunity to explore the behavior of massless Dirac fermions2,3. It is widely recognized that because of the conical band structure unique to massless Dirac fermions, graphene absorbs light in accordance with the fine-structure constant2,4,5. Considering the fact that graphite generally possesses massive Dirac fermions6, there should be interesting physics relevant to a change in N from ∞ to 1. Furthermore, the “mass” is closely related to stacking order and interlayer distance which are changed by thermal expansion caused by absorption of light. Therefore, it is meaningful to investigate the optical properties of an N-layer graphene as a function of N.
In this paper, we show the universal layer number in graphite to be 2/πα. By employing the transfer matrix method, we show simultaneously that for light in the infrared to visible range, an N-layer graphene with N above 1500 may be regarded as graphite, on the basis of a convergence criterion. This also enables us to obtain the simplest expression for the optical constants of graphite, which is written in terms of α and interlayer distance only5. The validity of the method is proved by the fact that the calculated reflectance reproduces experimentally observed reflectivity of graphite7. When similar considerations are applied to the general van der Waals heterostructures, many insights besides those into graphite will be obtained.
Results
Universal layer number and optical definition of graphite
Figure 1a shows the N-dependence of the absorptance AN. The solid curves are obtained by the transfer matrix method for different wavelengths, while the dashed curve at the top is from a simpler model that neglects the reflection. There are three noticeable points in this figure. First, for a small N value below 20, it is hard to distinguish between the results of the two models. This means that the reflection is suppressed and negligible for such N values, which are the systems that were studied extensively after the discovery of graphene. The simplified model is justified for a small N although it fails to explain the fact that graphite gleams. Second, there is a lump-like structure at around 90 that is almost independent of the wavelength. This suggests the presence of a universal number of layers. Figure 1b, which shows the N-dependence of the reflectance RN and transmittance TN for a wavelength at 1550 nm, reveals that the lump-like structure corresponds to the layer number at which RN and TN cross each other. In other words, because RN increases and TN decreases with increasing N, RN + TN takes an extremal value at the crossing point where the condition RN = TN is satisfied. To confirm that the crossing point is exactly given by a universal layer number N = 2/πα, we take the ω → 0 limit in Eq. (4) which immediately gives \({R}_{N}={(\frac{\pi \alpha N}{2+\pi \alpha N})}^{2}\) and \({T}_{N}={(\frac{2}{2+\pi \alpha N})}^{2}\). The emergence of the universal layer number in light absorption is intriguing, because the lump-like structure exhibits a maximum efficiency of 50% for a wide range of wavelengths (>1500 nm) at this universal layer number (2/πα ≈ 87). Physically, maximum efficiency is achieved by the fact that the electric field is not decaying and almost constant (due to boundary reflection). As a result, each layer absorbs an almost equal amount of light most efficiently. Third, AN is almost saturated above about 500 and 670 for 532 and 1550 nm, respectively. The critical layer number, defined more precisely by the tolerance ∣AN − A∞∣ ≤ 10−3, increases with decreasing photon energy, as shown in Fig. 1c. The critical layer number is closely related to the light frequency of interest, and the energy dependence must be relevant to that of the field localization, because the convergence of AN means that TN approaches 0. For light in a wide range of frequencies from infrared to visible, Fig. 1c clearly shows that for N > 1500, AN meets the convergence criterion. Therefore, the condition N > 1500 serves as the definition of “graphite” or more exactly, the lower limit of the number of layers that can be considered graphite from the optics point of view.
The results in Fig. 1a, b, c can be qualitatively understood with the basic knowledge about the absorption of light in thin films8. For a very thin film, absorption indeed scales linearly with the real part of the dynamical conductivity and thus also with the number of layers, shown by the simplified model. If the dynamical conductivity does not depend on the frequency, then also the absorption is necessarily dispersionless. Above a certain thickness, this linear approximation fails, while the wave nature of light starts to play a role (thickness is still small, but no longer fully negligible as compared with the wavelength). The absorption then becomes dispersive and it does not go linearly (not even monotonically) with the thickness anymore. Such behavior is well described in classical works on the optical properties of solids8. The only difference here is the non-dispersive conductivity of graphene due to interband transitions.
Simple formula for optical constant of graphite
We now derive a simple formula for the refractive index n and extinction coefficient κ of graphite. The proof consists of two logical steps. First, the optical constants are determined by the dynamical conductivity of graphite σgraphite as (n + iκ)2 = 1 + iσgraphite/ϵ0ω. As a result, we can obtain a theoretical Rgraphite by the standard formula; Rgraphite = [(n − 1)2 + κ2]/[(n + 1)2 + κ2], if σgraphite is known. Second, by the ansatz σgraphite = σgraphene/d5,9, it can be shown that the resultant Rgraphite exactly reproduces the infinitely large N of the transfer matrix model (A∞ = 1 − Rgraphite). Therefore, with the transfer matrix model we establish the following simple analytical expression of the optical constants of graphite,
These can be used to capture the changes in d caused, for example, by a high pressure10, temperature increase (thermal expansion) due to light absorption, and the propagation of a strain pulse11, through the reflectivity of graphite. We note that the frequency dependence of the critical layer number shown in Fig. 1c is explained by the field localization characterized by the skin-depth δ(ω) ≡ c/ωκ(ω) = (2/πα)n(ω)d, since n(ω) increases as ω decreases. Note also that it contains 2/πα as the characteristic layer number.
We note that the α value cannot change with pressure, so any pressure-induced change to the linear-band structure results in a change of the reflectance. But it is not due to a pressure-induced changes of α but due to a breakdown of the simple description of the reflectance in terms of α and d only.
Discussion
Here, we compare the calculation with an experiment for graphite (highly oriented pyrolithic graphite)7. By using the above-mentioned formula with the optical constants measured for an ordinary ray, we plot the reflectivity as a function of photon energy in Fig. 1d (dashed curve). We compare it with the result calculated with N = 1500, for which there is little discrepancy between theory and experiment at least for photon energy below 0.5 eV. The difference seen from 0.5 to 1.5 eV means extra absorption, which may be the effects of thermal expansion and diffusion12. This speculation is consistent with the fact that Rgraphite decreases with increasing d. The discrepancy between theory and experiment seen at high energy above 1.5 eV is suppressed by the inclusion of the correction πα(1 + 0.05(ℏω/t)2), where t = 3 eV is the hopping integral. The need to correct the conductivity arises partly because of the breakdown of the Dirac cone approximation at high energy13. Our model of graphite, based on the assumption that electrons at different layers interact only through electromagnetic fields14, satisfactory reproduces the experimental result for the reflectance of graphite, despite its simplicity. Our aim was to describe the optical properties of an N-layer graphene in a simplified situation, namely, in the absence of a substrate. In particular, for a small N with non-negligible transmittance, there is a possibility that the transmittance (and therefore the absorptance) is modified by a special substrate. The transfer matrix method can be modified to incorporate the effect of a substrate, which might be necessary if we are to understand the results for a complicated situation.
Finally, we show the correction to AN caused by an electronic coupling between layers that is expected to exist for natural graphite with AB stacking dominance10,15. The energy dispersion relation changes from massless Dirac fermions to massive Dirac fermions with the “mass” \({m}_{r}={\gamma }_{1}\cos \left(\frac{r\pi }{N+1}\right)\), where r (=1, ⋯ , N) is the wavenumber of the standing wave formed along the c-axis by an interlayer hopping of γ1 (= 0.4 eV)6,16. We calculated interband dynamical conductivity of an N-layer graphene with AB stacking and found that it is expressed in terms of the mass variables as
where Θ(x) denotes the step function satisfying Θ(x) = 1 for x ≥ 0 and 0 otherwise. This formula reproduces the results obtained previously for small N values17,18,19 and infinite N20. By defining the relative permittivity as \({\varepsilon }_{N}(\omega )\equiv 1+i\frac{{\sigma }_{N}(\omega )}{{\epsilon }_{0}\omega }\), we can obtain the reflection amplitude as \({c}_{0}^{r}=\frac{i \big(\sqrt{{\varepsilon }_{N}}-\frac{1}{\sqrt{{\varepsilon }_{N}}}\big)\sin \left(\sqrt{{\varepsilon }_{N}}\frac{\omega ({z}_{N}-{z}_{1})}{c}\right)}{2\cos \left(\sqrt{{\varepsilon }_{N}}\frac{\omega ({z}_{N}-{z}_{1})}{c}\right)-i\big(\sqrt{{\varepsilon }_{N}}+\frac{1}{\sqrt{{\varepsilon }_{N}}}\big)\sin \left(\sqrt{{\varepsilon }_{N}}\frac{\omega ({z}_{N}-{z}_{1})}{c}\right)}\). Figure 2 shows the resultant AN for three photon energies ℏω = 0.8, 0.4, and 0.2 eV. The correction to the lump-like structure is seen. Interestingly, the correction is most pronounced when photon energy is near 2γ1 and is suppressed by decreasing photon energy. The former originates from the fact that σN(ω) is enhanced by direct interband transitions at ℏω = 2γ110,15. The latter is relevant to the fact that the effect of the γ1 is removed from σN(ω) in the ω → 0 limit: indeed, \({\mathrm{lim}\,}_{\omega \to 0}{\sigma }_{N}(\omega )=\pi \alpha {\epsilon }_{0}c/d\) provided that the system is undoped. Our analysis presented here is limited to AB stacking, and there is a possibility that σN(ω) undergoes further correction because the low-energy band structure is perturbed by an interlayer interaction in a stacking dependent manner6,21.
In summary, we have proposed the concept of a universal number of layers 2/πα ≈ 87. When the optical properties of an N-layer graphene are viewed as a function of N, 2/πα is much larger than the layer numbers that have been studying extensively after the discovery of the graphene monolayer. Our approaches based on the simple approximation of separated, electrically uncoupled graphene N layers is applicable to other graphene-like materials, multilayer epitaxial graphene and turbostratic graphite22, besides highly oriented pyrolithic graphite. We extended the study to natural graphite, by including interlayer electronic coupling, and found that the N-dependence of absorptance would exhibit a correction to the universal layer number for specific photon energy. When similar considerations are applied to the general van der Waals heterostructures, many insights besides those into graphite will be obtained.
Methods
The transfer matrix method is useful for determining the propagation of waves such as electrons, light, and phonons23 in a superlattice and has also been used for a superlattice containing graphene24,25,26. Our superlattice is shown in Fig. 3, where the previously unreported key idea relates to the interlayer space between adjacent graphene layers as a constituent of the superlattice. In other words, graphite is a superlattice consisting of a vacuum and an atomic layer. The arrows along the z-axis indicate the propagation direction of light; the first graphene layer transmits and reflects the incident light in the forward and backward directions with the amplitudes \({c}_{1}^{t}\) and \({c}_{0}^{r}\), which are to be determined. Such transmission and reflection are repeated at each layer. The absolute square of \({c}_{N}^{t}\) and \({c}_{0}^{r}\) corresponds to the transmittance TN and reflectance RN, respectively. The total absorptance is given by AN = 1 − RN − TN, which is equivalent to the sum of the energies absorbed by each graphene layer.
The electromagnetic fields of adjacent layers are related by the transfer matrix T given by
The T matrix is obtained from Maxwell’s equations and is actually the product of two matrices. The first matrix expresses the propagation of electromagnetic fields in the interlayer vacuum space, which is fixed by the interlayer distance d, angular frequency of light ω, and speed of light c. The second matrix represents the boundary condition of graphene; the electric field is continuous (\({E}_{i+1}^{\prime}={E}_{i}\)), while the magnetic field is discontinuous (\({B}_{i+1}^{\prime}={B}_{i}-(\pi \alpha /c){E}_{i}\)). The discontinuity is given by the current in graphene, Ji = σgrapheneEi, where σgraphene = παϵ0c is the graphene’s dynamical conductivity, divided by −ϵ0c2 according to Ampére’s circuital law. By multiplying the transfer matrix with the field at an infinitesimal distance above the top layer N − 1 times, we can determine the field at an infinitesimal distance below the N-th layer,
Because \({E}_{0}({z}_{1})={e}^{i\omega {z}_{1}/c}+{c}_{0}^{r}{e}^{-i\omega {z}_{1}/c}\) and \(c{B}_{0}({z}_{1})={e}^{i\omega {z}_{1}/c}-{c}_{0}^{r}{e}^{-i\omega {z}_{1}/c}\) for the light at the entrance, and \({E}_{N}^{\prime}({z}_{N})={c}_{N}^{t}{e}^{i\omega {z}_{N}/c}\) and \(c{B}_{N}^{\prime}({z}_{N})={c}_{N}^{t}{e}^{i\omega {z}_{N}/c}\) for the field at the exit, Eq. (4) provides two equations for determining the amplitudes \({c}_{0}^{r}\) and \({c}_{N}^{t}\). Once \({c}_{0}^{r}\) is known, \({c}_{i}^{t}\) and \({c}_{i}^{r}\) can be calculated by using Eq. (3), and the field configuration is locally determined. Please note that in the ω → 0 limit in Eq. (4), we obtain
which immediately gives \({R}_{N}={\left(\frac{\pi \alpha N}{2+\pi \alpha N}\right)}^{2}\) and \({T}_{N}={\left(\frac{2}{2+\pi \alpha N}\right)}^{2}\).
For comparison with the transfer matrix model, we use a simpler model that neglects the reflection. The first layer absorbs πα of the incident light and transmits the remaining 1 − πα. Because the second layer absorbs πα of the transmitted light, the amount of transmitted light after the second layer becomes (1 − πα)2. By repeating the calculation, we find that the transmittance after N layers is given by (1 − πα)N and that the absorptance becomes 1 − (1 − πα)N. This result is independent of ω in contrast to AN.
Data availability
The data that support the findings of this study are available from the corresponding author upon reasonable request.
References
Novoselov, K. S. et al. Two-dimensional atomic crystals. Proc. Natl Acad. Sci. USA 102, 10451–3 (2005).
Novoselov, K. S. et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature 438, 197–200 (2005).
Zhang, Y., Tan, Y.-W., Stormer, H. L. & Kim, P. Experimental observation of the quantum Hall effect and Berryas phase in graphene. Nature 438, 201–204 (2005).
Nair, R. R. et al. Fine structure constant defines visual transparency of graphene. Science 320, 1308–1308 (2008).
Kuzmenko, A. B., van Heumen, E., Carbone, F. & van der Marel, D. Universal optical conductance of graphite. Phys. Rev. Lett. 100, 117401 (2008).
Min, H. & MacDonald, A. H. Electronic structure of multilayer graphene. Prog. Theor. Phys. Suppl. 176, 227–252 (2008).
Djurišić, A. B. & Li, E. H. Optical properties of graphite. J. Appl. Phys. 85, 7404–7410 (1999).
Palik, E. D. & Furdyna, J. K. Infrared and microwave magnetoplasma effects in semiconductors. Rep. Prog. Phys. 33, 307 (1970).
Blinowski, J. et al. Band structure model and dynamical dielectric function in lowest stages of graphite acceptor compounds. J. de. Phys. 41, 47–58 (1980).
Hanfland, M., Syassen, K. & Sonnenschein, R. Optical reflectivity of graphite under pressure. Phys. Rev. B 40, 1951–1954 (1989).
Thomsen, C., Grahn, H. T., Maris, H. J. & Tauc, J. Surface generation and detection of phonons by picosecond light pulses. Phys. Rev. B 34, 4129–4138 (1986).
Marsden, B., Mummery, A. & Mummery, P. Modelling the coefficient of thermal expansion in graphite crystals: implications of lattice strain due to irradiation and pressure. Proc. R. Soc. A: Math., Phys. Eng. Sci. 474, 20180075 (2018).
Stauber, T., Peres, N. M. R. & Geim, A. K. Optical conductivity of graphene in the visible region of the spectrum. Phys. Rev. B 78, 085432 (2008).
Painter, G. S. & Ellis, D. E. Electronic band structure and optical properties of graphite from a variational approach. Phys. Rev. B 1, 4747–4752 (1970).
Taft, E. A. & Philipp, H. R. Optical properties of graphite. Phys. Rev. 138, A197–A202 (1965).
Partoens, B. & Peeters, F. M. Normal and Dirac fermions in graphene multilayers: Tight-binding description of the electronic structure. Phys. Rev. B 75, 193402 (2007).
Abergel, D. S. L. & Fal’ko, V. I. Optical and magneto-optical far-infrared properties of bilayer graphene. Phys. Rev. B 75, 155430 (2007).
Koshino, M. & Ando, T. Magneto-optical properties of multilayer graphene. Phys. Rev. B 77, 115313 (2008).
Orlita, M. & Potemski, M. Dirac electronic states in graphene systems: optical spectroscopy studies. Semiconductor Sci. Technol. 25, 063001 (2010).
Ichikawa, Y. H. & Kobayashi, K. Optical properties of graphite in the infrared region. Carbon 3, 401–406 (1966).
Johnson, L. G. & Dresselhaus, G. Optical properies of graphite. Phys. Rev. B 7, 2275–2285 (1973).
Hass, J. et al. Why multilayer graphene on 4H–SiC (000\(\overline1\)) behaves like a single sheet of graphene. Phys. Rev. Lett. 100, 125504 (2008).
Mizuno, S., Ito, M. & Tamura, S.-i Time advance and delay of phonon packets scattered off superlattice systems. Jpn. J. Appl. Phys. 33, 2880–2885 (1994).
Kaipa, C. S. R. et al. Enhanced transmission with a graphene-dielectric microstructure at low-terahertz frequencies. Phys. Rev. B 85, 245407 (2012).
Zhan, T., Shi, X., Dai, Y., Liu, X. & Zi, J. Transfer matrix method for optics in graphene layers. J. Phys.: Condens. Matter 25, 215301 (2013).
Cheng, Y. H., Chen, C., Yu, K. Y. & Hsueh, W. J. Extraordinary light absorptance in graphene superlattices. Opt. Express 23, 28755 (2015).
Acknowledgements
We thank Dr. Keiko Kato for helpful discussions.
Author information
Authors and Affiliations
Contributions
K.S. and K.H. conceived of the presented idea. K.S. developed the theory and performed the computations. All authors discussed the results and contributed to the final paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Sasaki, K., Hitachi, K. Universal layer number in graphite. Commun Phys 3, 90 (2020). https://doi.org/10.1038/s42005-020-0354-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s42005-020-0354-y
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.