Corrections to the Reflectance of Graphene by Light Emission

1. Introduction

The visual detection of graphene on Si substrates involves a complex scientific aspect related to the interference effects of light [1,2,3]. To enhance visibility, it is crucial for the contrast between the reflectances from the substrate and the graphene layer to be sufficiently large. Si substrates with a specific thickness of ${\mathrm{SiO}}_2$ ($d_{{\mathrm{SiO}}_2}$) are known to provide advantages in this context [4,5]. Specifically, when the optical path length difference of the ${\mathrm{SiO}}_2$ layer is one-half or three-half the wavelength of the incident light, the contrast reaches its maximum.

The enhanced visibility of graphene is primarily attributed to a substrate-induced enhancement of light absorption. This enhancement is valuable not only for graphene detection but also for exploring notable phenomena. For instance, the reflectance of specific graphene multilayers on SiO₂/Si substrates can be reduced to zero for normally incident visible light with a wavelength ($\lambda$) approximately equal to $2d_{{\mathrm{SiO}}_2}$. Zero reflectance is achieved through destructive interference caused by ${\mathrm{SiO}}_2$ and a significant absorption coefficient of graphene [6,7,8]. More importantly, the reflectance of multilayer graphene on SiO₂/Si substrates is intricately determined because the significant absorption may result in non-negligible luminescence. This luminescence also functions as a secondary light source, leading to more sophisticated interference effects of light [9,10,11,12,13,14,15].

Since graphene lacks an energy bandgap, unlike semiconductors, it is naturally expected that the impact of light emission on its optical properties is negligible. However, electrons that have absorbed light undergo various relaxation processes towards the ground state, [16,17] generally involving luminescence contributions. For instance, Lui et al. [12] and Suemoto et al. [14] have reported that graphene and graphite emit light under femtosecond laser pulse excitation, respectively. In this paper, we demonstrate that, despite the tiny branching ratio of light emission to absorption, light emission constitutes the primary correction to visible reflectance, thanks to the substantial absorption of graphene. We have found this notable feature for the first time by leveraging both the destructive interference effects of substrates and the coherence of the multilayer structure. Additionally, we can show that light emission significantly contributes to enhancing the visibility of graphene on a substrate.

In this study, we delineate coherent and incoherent corrections to the reflectance of multilayer graphene (on SiO₂/Si substrates) caused by the light it emits after absorbing incident light. The coherent corrections are dependent on only two parameters: the branching ratio ($\mathcal{B}$) of coherent light emission to absorbed light and the phase (+ or −) of the emitted light. Both parameters are assumed to be independent of the number of layers (N) and the wavelength of light. Contrasts (reflectivities) calculated using these phenomenologically determined parameters show reasonable agreement with measured values for various N. Despite the branching ratio, determining the strength of light emitted from a graphene monolayer, being less than one percent, the coherent components become the primary corrections to the reflectance of multilayer graphene. This is due to the coherence increasing the amplitude of reflection through constructive interference. Furthermore, based on the observation that the phase of reflection and transmission coefficients of each layer translates into the amplitude of the emitted light, we hypothesize that the origin of the coherent components is coherent light emission stimulated by the incident light. Meanwhile, the Raman effect is explained as incoherent corrections. These conclusions are independent of the specific values of the two parameters.

Recently, the optical constants of single-layer graphene were measured with high accuracy [18,19]. However, classical electromagnetic models were employed for fitting and interpreting the results. Our approach advances the physical understanding of the measured optical constants of single-layer graphene to a more fundamental level, incorporating principles of quantum electrodynamics that describe the creation of photons (light emission). Moreover, our formulation is versatile and generally applicable to any layered material and superlattice, promising an accurate description of their optical processes.

Multilayer graphene composed of tens to hundreds of layers is an interesting research subject. However, it has not been explored much, partly because the success of the exfoliation method has rapidly shifted the interest of many researchers from infinite layers of graphite to few-layer graphene. The obvious advantage of multilayer graphene is that it can increase the signal strength, [20] whereas the signal strength of few-layer graphene is low and difficult to measure. Besides that, multilayer graphene hosts various intriguing phenomena. For instance, studies have revealed that the optimal number of layers for absorbing infrared radiation is 87 [21]. Additionally, a notable nonlinear optical effect has been attributed to multilayer graphene [22,23,24]. Yang et al., for instance, demonstrated the maximum third-harmonic signal from 24-layer graphene on a quartz substrate [25].

This paper is structured as follows: In Section 2, we present fundamental insights into the reflectance of multilayer graphene on SiO₂/Si substrates. By utilizing the Fresnel equation, we can replicate measured contrasts and assert that the optical constants obtained by El-Sayed et al. [19] sufficiently and accurately describe the reflectance. However, it is noteworthy that these optical constants lead to an unexpectedly large effective coupling between light and graphene, a phenomenon inconsistent with existing experiments. Section 4 introduces our formulation of corrections to the reflectance arising from light emissions. Our model effectively describes the reflectance without introducing such inconsistencies. These corrections are categorized into two types: coherent emission, corresponding to stimulated emission with a common phase, and incoherent emission, featuring a random phase identified as the Raman effect. A detailed comparison between calculated and measured contrasts is presented in Section 5. Finally, Section 6 offers a discussion of the findings, and Section 7 provides the conclusions.

2. Reflectance of Graphene Multilayer on Si Substrate

In this section, we demonstrate that the measured reflectance of various graphene multilayers with different thickness aligns well with the Fresnel equation when adopting optical constants (n and k) obtained by El-Sayed et al. [19] through ellipsometric measurements of chemical vapor deposited (CVD) graphene monolayer. The key observation from this section is that electron-light coupling constant, extracted from the established n and k values, deviates significantly from the fine-structure constant $\alpha \simeq 1/137$. In fact, the difference is beyond the level of small corrections of order of ${\alpha }^2$ caused by such as a change in the band dispersion at high energy from linear dispersion (commonly known as trigonal warping effects) and Fermi velocity renormalization, [26,27] which have been discussed theoretically thus far. This observation partly motivates the introduction of a new reflectance model developed in the subsequent sections.

2.1. Fresnel Equation

The reflectance of N-layer graphene on SiO₂/Si substrate to normally incident light of wavelength $\lambda$ can be formulated using the reflection coefficient $r_N\left(\lambda \right)$ as $R_N\left(\lambda \right)={\left|r_N\left(\lambda \right)\right|}^2$ (Fresnel equation), where

$$\begin{array}{cc}\hfill & r_N\left(\lambda \right)=\hfill \\ \hfill & \frac{\left[(1-n_{\mathrm{Si}})cos\phi -i\left(\frac{n_{\mathrm{Si}}}{n_{{\mathrm{SiO}}_2}}-n_{{\mathrm{SiO}}_2}\right)sin\phi \right]cos\varphi -\left[\left(\frac{n_{{\mathrm{SiO}}_2}}{\sqrt{{\epsilon }_{\mathrm{g}}}}-\frac{n_{\mathrm{Si}}\sqrt{{\epsilon }_{\mathrm{g}}}}{n_{{\mathrm{SiO}}_2}}\right)sin\phi +i\left(\frac{n_{\mathrm{Si}}}{\sqrt{{\epsilon }_{\mathrm{g}}}}-\sqrt{{\epsilon }_{\mathrm{g}}}\right)cos\phi \right]sin\varphi }{\left[(1+n_{\mathrm{Si}})cos\phi -i\left(\frac{n_{\mathrm{Si}}}{n_{{\mathrm{SiO}}_2}}+n_{{\mathrm{SiO}}_2}\right)sin\phi \right]cos\varphi -\left[\left(\frac{n_{{\mathrm{SiO}}_2}}{\sqrt{{\epsilon }_{\mathrm{g}}}}+\frac{n_{\mathrm{Si}}\sqrt{{\epsilon }_{\mathrm{g}}}}{n_{{\mathrm{SiO}}_2}}\right)sin\phi +i\left(\frac{n_{\mathrm{Si}}}{\sqrt{{\epsilon }_{\mathrm{g}}}}+\sqrt{{\epsilon }_{\mathrm{g}}}\right)cos\phi \right]sin\varphi }.\hfill \end{array}$$

(1)

Here, $\phi \equiv n_{{\mathrm{SiO}}_2}{d}_{{\mathrm{SiO}}_2}\frac{2\pi }{\lambda }$ is the phase acquired by light after it propagates through a distance $d_{{\mathrm{SiO}}_2}$ in SiO₂, and $\varphi \equiv \sqrt{{\epsilon }_{\mathrm{g}}}\frac{2\pi }{\lambda }{d}_N$ is the complex phase acquired when light passes through N-layer graphene of thickness $d_N\equiv Nd$. Multilayer graphene is treated as an effective medium whose unit length is the interlayer spacing d ($=0.335$ nm) and its dielectric constant is given by ${\epsilon }_g={(n+ik)}^2$ [19]. $n_{\mathrm{Si}}$ and $n_{{\mathrm{SiO}}_2}$ are the refractive indexes of Si and SiO₂, respectively. Si is treated as an absorbing substrate having a semi-infinite thickness ($n_{\mathrm{Si}}$ is a complex number) whose dispersion is taken into account, [28] while SiO₂ is treated as an absorption-free film [29]. $R_N\left(\lambda \right)$ depends sensitively on the two phases, $\phi$ and $\varphi$. When $N=0$ or $\varphi =0$ in Equation (1), $R_0\left(\lambda \right)$ corresponds to the reflectance of the substrate without graphene. It can be minimized for a specific $\lambda$ by destructive interference; namely, $R_0\left(\lambda \right)$ is at a minimum when $cos\phi =0$ as ${\left|\frac{n_{\mathrm{Si}}-n_{{\mathrm{SiO}}_2}^2}{n_{\mathrm{Si}}+n_{{\mathrm{SiO}}_2}^2}\right|}^2$ [30]. Monolayer graphene is most easily detectable on SiO₂/Si substrates when destructive interference occurs, because $|R_0\left(\lambda \right)-R_1\left(\lambda \right)|$ takes a maximum when $cos\phi =0$ (i.e., when $\lambda \simeq 2d_{\mathrm{Si}{0}_2}$ because $n_{{\mathrm{SiO}}_2}\simeq 1.46$) [4].

2.2. Comparison of Measured and Calculated Contrasts

Multilayer graphene was prepared by exfoliating highly oriented pyrolytic graphite (HOPG) on the same SiO₂/Si substrate. The reflectance of the multilayer graphene was measured with a spectroscopic reflectometer (TohoSpec3100, Toho Technology) using a $\times 50$ objective lens. First, we determined that $d_{{\mathrm{SiO}}_2}=268$ nm from the reflectance of the substrate (Appendix A). This value is used consistently in all the calculations reported in this paper, and it results in that destructive interference occurs for $R_0$ at $\lambda \simeq 520$ nm. The thickness of the graphene flakes was determined by atomic force microscopy (Dimension XR, Bruker).

The representative measured spectral contrasts are depicted as black dots (circles) in Figure 1 (The error bars for the data are within each circle). It is important to note that we present contrasts ($C_N\equiv R_N/R_0$) instead of the reflectivities ($R_N$) to prevent any artificial shifts in the reflectivities (see Appendix A for more details). The general feature of the spectral shapes can be elucidated as follows: for thin samples with fewer than 40 layers, the contrast is subdued due to destructive interference from the substrate, resulting in a concave structure near $\lambda =520$ nm. For thick samples with over 60 layers, $R_N$ is predominantly influenced by contributions from the N-layer graphene and is minimally impacted by the substrate. Consequently, given that $R_0$, suppressed by destructive interference, is in the denominator of the contrast, a convex structure appears near $\lambda =520$ nm.

The contrasts, calculated using Equation (1), are represented by the red solid curves in Figure 1. A satisfactory fit with no significant deviations is achieved for various multilayer graphene samples, except for $C_1$ and $C_3$. The observed discrepancy in these thin samples likely arises from differences in the interface between graphene and the substrate compared to that between adjacent graphene layers. This discrepancy is specific to thin samples and becomes negligible at 6 layers. The notable agreement between measured and calculated reflectance for visible light clearly indicates the following two facts. First, multilayer graphene can be treated as a collection of independent single layers, which is plausible since the effect of interlayer stacking does not manifest in the dynamical conductivity within the visible light range. Second, reflectance measurements provide a reliable value of layer number N, averaged within a light spot. This is true even when the thickness of the graphene flakes determined by atomic force microscopy shows positional fluctuations to a certain extent. These facts will be used to validate the underlying assumption of our theoretical model of reflection presented in the subsequent sections.

$C_{14}$ and $C_{15}$ vanish at $\lambda \simeq 2d_{\mathrm{Si}{0}_2}$, which has been referred to as zero reflection [5,30]. Zero reflection is useful in knowing the values of basic parameters later (at the end of Section 3). To capture the essential role of graphene in achieving zero reflection, let us consider Equation (1) when $cos\phi =0$. For the numerator to vanish, $i(n_{\mathrm{Si}}-n_{{\mathrm{SiO}}_2}^2)cos\varphi -\left(\frac{n_{{\mathrm{SiO}}_2}^2}{\sqrt{{\epsilon }_{\mathrm{g}}}}-n_{\mathrm{Si}}\sqrt{{\epsilon }_{\mathrm{g}}}\right)sin\varphi =0$ has to be satisfied. Since $\varphi$ is small, this equation can be simplified as $i(n_{\mathrm{Si}}-n_{{\mathrm{SiO}}_2}^2)-\left(n_{{\mathrm{SiO}}_2}^2-n_{\mathrm{Si}}{\epsilon }_{\mathrm{g}}\right)\frac{d}{\lambda }2\pi N=0$, which shows that $N\sim \frac{i{d}_{{\mathrm{SiO}}_2}}{\pi {\epsilon }_{\mathrm{g}}d}(1-\frac{n_{{\mathrm{SiO}}_2}^2}{n_{\mathrm{Si}}})$ is an approximate layer number that gives zero reflection. This argument makes it easy to understand that the dominant imaginary part of ${\epsilon }_{\mathrm{g}}$ is essential for zero reflection to occur.

2.3. Effective Coupling Constant

As the difference between the calculated and experimental values of contrasts proves to be sufficiently small for various multilayer graphene samples when utilizing the experimental ${\epsilon }_g={(n+ik)}^2$ values in Equation (1), the optical constants [19] are the results that theory should ultimately elucidate. It can be inferred that nearly all the optical information of graphene multilayers is encompassed in the optical constant of monolayer graphene. This inference is partially attributed to the fact that the effect of the stacking order on the reflectance does not manifest in the visible regime. Consequently, we must inquire to what extent the theory of graphene optics can account for the n and k values.

Because of a conical energy-band structure of graphene known as the Dirac cone, the dynamical conductivity is well approximated by $\pi \alpha$ for visible light [31]. As a result, (suspended) monolayer graphene absorbs $\sim 2.3$ percent ($=\pi \alpha$) of the incident visible light [32]. A straightforward calculation of the Kubo formula shows that the dynamical conductivity of graphite is given by that of graphene divided by the interlayer spacing d ($=0.335$ nm): ${\sigma }_{\mathrm{graphite}}=\pi \alpha /d$ [21,33,34]. The reflectance in the visible regime is free from the effects of the stacking order, [35,36,37] Fermi energy position, and temperature at room temperature ranges [38]. Thus, the relative permittivity of graphite for visible light wavelengths $\lambda$ is written as

$$\begin{array}{c}\hfill {\epsilon }_g={\epsilon }_r+i\frac{\alpha \lambda }{2d},\end{array}$$

(2)

where ${\epsilon }_r$ is the dielectric constant of the interlayer space.

By equating the right-hand side of Equation (2) with ${(n+ik)}^2$, we define an effective coupling constant ${\alpha }_{\mathrm{eff}}=4nkd/\lambda$, which is compared with $\alpha$ in Figure 2a. The difference between ${\alpha }_{\mathrm{eff}}$ and $\alpha$ is actually larger than the order of 10 percent of $\alpha$ which is beyond the level of small corrections of order of ${\alpha }^2$ considered theoretically so far. Similarly, we define ${\epsilon }_{\mathrm{eff}}=n^2-k^2$ and plot it in Figure 2b. If the interlayer space is a vacuum, an appropriate choice of ${\epsilon }_r$ would be 1. However, ${\epsilon }_{\mathrm{eff}}$ is very different from unity because the electronic wave function of the $\pi$-orbital spreads into the interlayer space, light propagating in it is subjected to the spread of the wave function [39,40]. Fang et al. [39] calculate ${\epsilon }_r=6.9$ using a microscopic Poisson equation which has been tested by experiment. [40] The calculated effective dielectric thickness of graphene is found to be 0.22 nm, and the microscopic dielectric permittivity decays from 6.9 in the carbon-atom plane to the vacuum permittivity within approximately 0.1 nm. Though the calculation is for the electric field pointing in the direction perpendicular to the graphene sheet, similar (but slightly smaller) value is expected for the direction parallel to the sheet (let us assume it is $4\sim 6$ here).

3. Corrections to Reflectance by Light Emission

A substantial correction to the dynamical conductivity is deemed impermissible because, if allowed, it would result in an inconsistency with the experimental observation that a (suspended) monolayer graphene absorbs approximately 2.3 percent ($=\pi \alpha$) of the incident visible light [32]. An accurate theory capable of describing the reflectance does not necessitate corrections to dynamical conductivity but rather requires a mechanism that explains reflectance without altering $\alpha$. The model of light emission from graphene introduced in this section serves as an illustration of such a mechanism. Our model inherently incorporates the crucial concept of wave interference, specifically coherence or incoherence, as the initial phase of light emitted from each graphene layer.

3.1. Basic Idea

Figure 5a illustrates our model of reflection, where horizontal lines on the substrate represent N-layer graphene, and the vertical lines depict light rays with arrows indicating the directions of light propagation. The light rays on the left side (black in color) of Figure 5a show the primary processes of reflection (excluding contributions from light emission). In this process, incident light from a light source is transmitted and reflected by graphene, while some energy of light being absorbed by each layer. The reflection coefficient, $r_N$, is calculated from a primary model which is defined in Section 4.2. The light rays on the right side (red in color) correspond to the light emission. Suppose that the jth layer emits light. The emitted light is transmitted and reflected by graphene until the light escapes the system, and it contributes to the reflectance of the system. Thus, there is another “reflection coefficient” when N-layer graphene emits light which is defined in Section 4.3. Let $z_N$ denotes the sum over such amplitudes from all layers. Once we know what $r_N$ and $z_N$ are, then the reflectance is given by $R_N={|r_N+z_N|}^2$.

We use a transfer matrix method to calculate $r_N$ and $z_N$. [21,34] Transfer matrix method is useful in calculating reflection (up arrow) and transmission (down) coefficients at each layer [$r_N\left(j\right)$ and $t_N\left(j\right)$ in Figure 5a] in addition to the electric field $E_j^N$ ($j=1,\dots ,N$) that determines the absorption of the jth layer as $A_j^N\equiv \pi \alpha {\left|E_j^N\right|}^2$. The total absorption of N-layer graphene is ${\sum }_{j=1}^N{A}_j^N$. [21]

3.2. Primary model

Basically, primary model means the Fresnel equation of Equation (1) where ${\epsilon }_g$ is given by Equation (2). ${\epsilon }_g$ has a single unknown constant ${\epsilon }_r$ which can be estimated from the following observations. Ultimately, we conclude that ${\epsilon }_r=4.5$.

First, there must not be a large discrepancy between ${\epsilon }_g$ and the experimentally determined optical constants [19,41]. Experimental values for n and k [19] are shown as • and ∘ in Figure 6a. The lines depict (bare) optical constants plotted using Equation (2) with the refractive index $n_g=\mathrm{Re}\left[\sqrt{{\epsilon }_g}\right]$ and absorption coefficient $k_g=\mathrm{Im}\left[\sqrt{{\epsilon }_g}\right]$ for ${\epsilon }_r=4.5$ and 5.5. When ${\epsilon }_r=4.5$, $k\simeq k_g$ but n has a certain difference from $n_g$. When ${\epsilon }_r=5.5$, $n\simeq n_g$ but k has a certain difference from $k_g$. Because $n_g$ increases with increasing ${\epsilon }_r$ while $k_g$ decreases, there is no ${\epsilon }_r$ value that can reproduce n and k simultaneously. This suggests that there should be such a discrepancy between them which is attributed to the corrections by light emission. Second, ${\epsilon }_g$ has to roughly reproduce the behavior of the reflectance of multilayer graphene. It exhibits a minimum at a certain wavelength, primarily due to destructive interference caused by SiO₂. As shown in Figure 6b, the position is red-shifted by increasing N, indicating that even thin graphite samples significantly impact the light interference effect. When ${\epsilon }_r=1$, the position changes little, and a sizable artificial shift in wavelength is needed to ensure consistency between theory and experiment, which cannot be explained as a correction. When ${\epsilon }_r=4.5$, a small difference between theory and experiment still remains. However, as we show later, the corrections provide better agreement not only for the wavelengths giving minimum reflectivity but also for the minimum reflectivity values, thus accounting for the difference. Third, ${\epsilon }_r=4.5$ roughly reproduces the reflectivity of graphite in the infrared region [42]. Similar ${\epsilon }_r$ values have been used to reproduce the observed reflectivities of graphite and graphene [18,19,41]. We note that the value of ${\epsilon }_r$ is less than the magnitude of the imaginary part of ${\epsilon }_g$, since visible light has a much longer $\lambda$ (400∼800 nm) than d, although $\alpha$ is certainly a small quantity. The optical properties of multilayer graphene are thus characterized by the large imaginary part of ${\epsilon }_g$.

Figure 4. ${\epsilon }_r$ of the primary model (a) Experimental n and k values are represented by dots, taken from Ref. [19]. The lines depict bare optical constants without corrections. (b) Dots indicate measured wavelengths corresponding to the minimum reflectance. A comparison between the measured and calculated results suggests that ${\epsilon }_r=4.5$ is a reasonable value.

Figure 4. ${\epsilon }_r$ of the primary model (a) Experimental n and k values are represented by dots, taken from Ref. [19]. The lines depict bare optical constants without corrections. (b) Dots indicate measured wavelengths corresponding to the minimum reflectance. A comparison between the measured and calculated results suggests that ${\epsilon }_r=4.5$ is a reasonable value.

4. Corrections to Reflectance by Light Emission

A substantial correction to the dynamical conductivity is deemed impermissible because, if allowed, it would result in an inconsistency with the experimental observation that a (suspended) monolayer graphene absorbs approximately 2.3 percent ($=\pi \alpha$) of the incident visible light [32]. An accurate theory capable of describing the reflectance does not necessitate corrections to dynamical conductivity but rather requires a mechanism that explains reflectance without altering $\alpha$. The model of light emission from graphene introduced in this section serves as an illustration of such a mechanism. Our model inherently incorporates the crucial concept of wave interference, specifically coherence or incoherence, as the initial phase of light emitted from each graphene layer.

4.1. Basic Idea

Figure 5a illustrates our model of reflection, where horizontal lines on the substrate represent N-layer graphene, and the vertical lines depict light rays with arrows indicating the directions of light propagation. The light rays on the left side (black in color) of Figure 5a show the primary processes of reflection (excluding contributions from light emission). In this process, incident light from a light source is transmitted and reflected by graphene, while some energy of light being absorbed by each layer. The reflection coefficient, $r_N$, is calculated from a primary model which is defined in Section 4.2. The light rays on the right side (red in color) correspond to the light emission. Suppose that the jth layer emits light. The emitted light is transmitted and reflected by graphene until the light escapes the system, and it contributes to the reflectance of the system. Thus, there is another “reflection coefficient” when N-layer graphene emits light which is defined in Section 4.3. Let $z_N$ denotes the sum over such amplitudes from all layers. Once we know what $r_N$ and $z_N$ are, then the reflectance is given by $R_N={|r_N+z_N|}^2$.

Figure 5. Model description (a) The primary (left side) and secondary (right side) processes involved in the reflection are physically interconnected through light absorption at each graphene layer. (b) The self-consistent calculation of $Z_j$ is explained in detail in the text.

Figure 5. Model description (a) The primary (left side) and secondary (right side) processes involved in the reflection are physically interconnected through light absorption at each graphene layer. (b) The self-consistent calculation of $Z_j$ is explained in detail in the text.

We use a transfer matrix method to calculate $r_N$ and $z_N$ [21,34]. Transfer matrix method is useful in calculating reflection (up arrow) and transmission (down) coefficients at each layer [$r_N\left(j\right)$ and $t_N\left(j\right)$ in Figure 5a] in addition to the electric field $E_j^N$ ($j=1,\dots ,N$) that determines the absorption of the jth layer as $A_j^N\equiv \pi \alpha {\left|E_j^N\right|}^2$. The total absorption of N-layer graphene is ${\sum }_{j=1}^N{A}_j^N$. [21]

4.2. Primary Model

Basically, primary model means the Fresnel equation of Equation (1) where ${\epsilon }_g$ is given by Equation (2). ${\epsilon }_g$ has a single unknown constant ${\epsilon }_r$ which can be estimated from the following observations. Ultimately, we conclude that ${\epsilon }_r=4.5$.

First, there must not be a large discrepancy between ${\epsilon }_g$ and the experimentally determined optical constants [19,41]. Experimental values for n and k [19] are shown as • and ∘ in Figure 4a. The lines depict (bare) optical constants plotted using Equation (2) with the refractive index $n_g=\mathrm{Re}\left[\sqrt{{\epsilon }_g}\right]$ and absorption coefficient $k_g=\mathrm{Im}\left[\sqrt{{\epsilon }_g}\right]$ for ${\epsilon }_r=4.5$ and 5.5. When ${\epsilon }_r=4.5$, $k\simeq k_g$ but n has a certain difference from $n_g$. When ${\epsilon }_r=5.5$, $n\simeq n_g$ but k has a certain difference from $k_g$. Because $n_g$ increases with increasing ${\epsilon }_r$ while $k_g$ decreases, there is no ${\epsilon }_r$ value that can reproduce n and k simultaneously. This suggests that there should be such a discrepancy between them which is attributed to the corrections by light emission. Second, ${\epsilon }_g$ has to roughly reproduce the behavior of the reflectance of multilayer graphene. It exhibits a minimum at a certain wavelength, primarily due to destructive interference caused by SiO₂. As shown in Figure 4b, the position is red-shifted by increasing N, indicating that even thin graphite samples significantly impact the light interference effect. When ${\epsilon }_r=1$, the position changes little, and a sizable artificial shift in wavelength is needed to ensure consistency between theory and experiment, which cannot be explained as a correction. When ${\epsilon }_r=4.5$, a small difference between theory and experiment still remains. However, as we show later, the corrections provide better agreement not only for the wavelengths giving minimum reflectivity but also for the minimum reflectivity values, thus accounting for the difference. Third, ${\epsilon }_r=4.5$ roughly reproduces the reflectivity of graphite in the infrared region [42]. Similar ${\epsilon }_r$ values have been used to reproduce the observed reflectivities of graphite and graphene [18,19,41]. We note that the value of ${\epsilon }_r$ is less than the magnitude of the imaginary part of ${\epsilon }_g$, since visible light has a much longer $\lambda$ (400∼800 nm) than d, although $\alpha$ is certainly a small quantity. The optical properties of multilayer graphene are thus characterized by the large imaginary part of ${\epsilon }_g$.

Figure 6. ${\epsilon }_r$ of the primary model (a) Experimental n and k values are represented by dots, taken from Ref. [19]. The lines depict bare optical constants without corrections. (b) Dots indicate measured wavelengths corresponding to the minimum reflectance. A comparison between the measured and calculated results suggests that ${\epsilon }_r=4.5$ is a reasonable value.

Figure 6. ${\epsilon }_r$ of the primary model (a) Experimental n and k values are represented by dots, taken from Ref. [19]. The lines depict bare optical constants without corrections. (b) Dots indicate measured wavelengths corresponding to the minimum reflectance. A comparison between the measured and calculated results suggests that ${\epsilon }_r=4.5$ is a reasonable value.

4.3. Model of Light Emission

The corrections to the reflectance are the main subjects of this paper [9,10]. Specifically, we consider corrections where some fraction of the energy absorbed by the jth layer (of N-layer graphene) is transferred to light emitted from that layer [see the right side of Figure 3(a)]. The amplitude of the emitted light is assumed to be the square root of the layer absorption $A_j^N\equiv \pi \alpha {\left|E_j^N\right|}^2$,34] multiplied by the branching ratio, $\mathcal{B}$, i.e., $\sqrt{\left(\mathcal{B}/2\right)A_j^N}$, where $1/\sqrt{2}$ means that the light emission is direction-independent along the c-axis. Note that $A_j^N$ depends not only on j and N but also on $\lambda$ and $d_{{\mathrm{SiO}}_2}$.

To examine how light emitted from the jth layer affects the reflectance, we define two subsystems, as shown in Figure 7b: one is an isolated $(j-1)$-layer graphene in the air; the other is $(N-j)$-layer graphene on SiO₂/Si substrate. Using the transfer matrix method, we can obtain the transmission and reflection coefficients of an isolated $(j-1)$-layer graphene in the air [denoted as $t_{j-1}^g$ and $r_{j-1}^g$] and the reflection coefficient of $(N-j)$-layer graphene on the SiO₂/Si substrate [denoted as $r_{N-j}$] [21,34]. Let the reflection coefficients be $X_j$ and $Y_j$ and transmission coefficient be $Z_j$ for the combined subsystems [see Figure 7b]. These can be obtained by a self-consistent manner as follows. After calculating $X_j^{\left(n\right)}$, we add it to $\sqrt{\left(\mathcal{B}/2\right)A_j^N}$ of the incident light to the $(j-1)$-layer graphene (in the air) as $\sqrt{\left(\mathcal{B}/2\right)A_j^N}+X_j^{\left(n\right)}$ and recalculate $Y_j^{(n+1)}=r_{j-1}^g\left(\sqrt{\left(\mathcal{B}/2\right)A_j^N}+X_j^{\left(n\right)}\right)$ and $Z_j^{(n+1)}=t_{j-1}^g\left(\sqrt{\left(\mathcal{B}/2\right)A_j^N}+X_j^{\left(n\right)}\right)$. Then, we add a new $Y_j^{(n+1)}$ to $\sqrt{\left(\mathcal{B}/2\right)A_j^N}$ of the light incident to the $(N-j)$-layer graphene on the SiO₂/Si substrate as $\sqrt{\left(\mathcal{B}/2\right)A_j^N}+Y_j^{(n+1)}$ and recalculate $X_j^{(n+1)}=r_{N-j}\left(\sqrt{\left(\mathcal{B}/2\right)A_j^N}+Y_j^{(n+1)}\right)$. These computations are repeated until $X_j$ and $Y_j$ converge. In this way, we can obtain analytical expressions for the converged $X_j$, $Y_j$, and $Z_j$ for a given $\mathcal{B}$:

$$\begin{array}{cc}\hfill & X_j\left(\mathcal{B}\right)=r_{N-j}\left\{\frac{1+r_{j-1}^g}{1-r_{N-j}{r}_{j-1}^g}\right\}\sqrt{\left(\frac{\mathcal{B}}{2}\right)A_j^N},\hfill \\ \hfill & Y_j\left(\mathcal{B}\right)=r_{j-1}^g\left\{\frac{1+r_{N-j}}{1-r_{N-j}{r}_{j-1}^g}\right\}\sqrt{\left(\frac{\mathcal{B}}{2}\right)A_j^N},\hfill \\ \hfill & Z_j\left(\mathcal{B}\right)=t_{j-1}^g\left\{\frac{1+r_{N-j}}{1-r_{N-j}{r}_{j-1}^g}\right\}\sqrt{\left(\frac{\mathcal{B}}{2}\right)A_j^N}.\hfill \end{array}$$

(3)

The “corrected” electric fields at an infinitesimal distance above and below the jth layer become $\sqrt{\left(\mathcal{B}/2\right)A_j^N}+X_j^{\left(n\right)}+Y_j^{(n+1)}$ and $\sqrt{\left(\mathcal{B}/2\right)A_j^N}+Y_j^{\left(n\right)}+X_j^{(n+1)}$, respectively. Self-consistency, whereby ${lim}_{n\to \infty }{X}_j^{\left(n\right)}=X_j$ and ${lim}_{n\to \infty }{Y}_j^{\left(n\right)}=Y_j$, is therefore essential to ensuring that the corrected electric field is continuous at the jth layer, which is a requirement of Maxwell equations. The corrected amplitude of the emitted light is written as

$$\begin{array}{c}\hfill \sqrt{\left(\frac{\mathcal{B}}{2}\right)A_j^N}+X_j+Y_j=(1+r_{j-1}^g)\left\{\frac{1+r_{N-j}}{1-r_{N-j}{r}_{j-1}^g}\right\}\sqrt{\left(\frac{\mathcal{B}}{2}\right)A_j^N}.\end{array}$$

(4)

By comparing this with $Z_j\left(\mathcal{B}\right)$, we see that more accurate value of the amplitude of the emitted light is given by multiplying $\{\cdots \}\sqrt{\left(\mathcal{B}/2\right)A_j^N}$ with $1+r_{j-1}^g$ as the renormalization constant, and $(1+r_{j-1}^g)\{\cdots \}\sqrt{\left(\mathcal{B}/2\right)A_j^N}$ is what $b_N\left(j\right)$ in Figure 3(a) represents. Therefore, we redefine $Z_j$ as

$$\begin{array}{cc}\hfill & Z_j\left(\mathcal{B}\right)\equiv t_{j-1}^g{b}_N\left(j\right),\hfill \end{array}$$

(5)

$$\begin{array}{cc}\hfill & b_N\left(j\right)\equiv (1+r_{j-1}^g)\left\{\frac{1+r_{N-j}}{1-r_{N-j}{r}_{j-1}^g}\right\}\sqrt{\left(\frac{\mathcal{B}}{2}\right)A_j^N}.\hfill \end{array}$$

(6)

We can interpret $Z_j\left(\mathcal{B}\right)$ as follows. The transmission coefficient $t_{j-1}^g$ is the direct propagation of the renormalized light emitted from jth layer to the air, and $|t_{j-1}^g{|}^2$ monotonously decreases with increasing j [21]. The effects of scattering and absorption of the emitted light caused by surrounding layers is included by the part in the brace $\{\cdots \}$. It tends to suppress the magnitude of $Z_j$, but sometimes enhance. For example, when $N=1$ (i.e., monolayer on a substrate), the part becomes $1+r_0$ which is larger than unity when $r_0$ is positive.

$Z_j\left(\mathcal{B}\right)$ is the value at zero initial phase, so the transmission coefficient can be given a phase degree of freedom expressing the coherence or incoherence of the light emission from the different layers:

$$\begin{array}{c}\hfill z_N\equiv \sum _{j=1}^N{e}^{i{\theta }_j}{Z}_j\left(\mathcal{B}\right).\end{array}$$

(7)

Accordingly, the corrected reflectance is uniquely determined by $R_N\equiv {\left|r_N+z_N\right|}^2={\left|r_N\right|}^2+2\mathrm{Re}\left[r_N{z}_N^{*}\right]+{\left|z_N\right|}^2$. The value of $z_N$ depends on these phases ${\theta }_j$ [43].

We consider a case in which the phase is given by a coherent phase. The exact derivation of the phase will be shown elsewhere because it is beyond the scope of the present paper. Here, we concisely explain the basic logic leading to the coherent phase in terms of quantum electrodynamics. First, we can define a quantum mechanical state of light ($|{\Psi }_a\rangle$) that the primary model describes (see left side of Figure 5). All the information of light is expressed by the coefficients $r_N\left(j\right)$ and $t_N\left(j\right)$ ($j=1,\dots ,N$). Second, we can also define another quantum state of light ($|{\Psi }_b\rangle$) for the emitted light (see right side of Figure 5). All the information of emitted light is expressed by the coefficients $b_N\left(j\right)$. These two states have an overlap $b_N^{*}\left(j\right)t_N\left(j\right)+b_N^{*}\left(j\right)r_N\left(j\right)$ caused by jth layer graphene. Thus, if we consider a linear superposition of these states as $|{\Psi }_a\rangle +e^{i\theta }|{\Psi }_b\rangle$ to form energy eigenstates, the phase $e^{i{\theta }_j}$ must be chosen so that $e^{-i{\theta }_j}{b}_N^{*}\left(j\right)(t_N\left(j\right)+r_N\left(j\right))$ becomes a real number, namely

$$\begin{array}{c}\hfill e^{i{\theta }_j}=\pm \frac{t_N\left(j\right)+r_N\left(j\right)}{|t_N\left(j\right)+r_N\left(j\right)|}\frac{b_N{\left(j\right)}^{*}}{|b_N\left(j\right)|}.\end{array}$$

(8)

The factor ± is a global phase ($\theta$) in the sense that it is independent of the value of j. Because the scattered light ($r_N$) and the emitted light ($z_N$) form a two-level state, there are two possible linear superpositions of their energy eigenstates, $-1$ ($\theta =\pi$) or $+1$ ($\theta =0$). The minus sign ($e^{i\pi }$) is assigned to the lower energy state. From Equations (5) and (8), we obtain $e^{i{\theta }_j}{Z}_j\left(\mathcal{B}\right)=\pm \frac{t_N\left(j\right)+r_N\left(j\right)}{|t_N\left(j\right)+r_N\left(j\right)|}{t}_{j-1}^g\left|b_N\left(j\right)\right|$.

Including the correction due to coherent light emission leads to

$$\begin{array}{c}\hfill R_N={\left|r_N+\sum _{j=1}^N{e}^{i{\theta }_j}{Z}_j\left({\mathcal{B}}_{coh}\right)\right|}^2,\end{array}$$

(9)

where ${\mathcal{B}}_{coh}$ is the branching ratio of the energy of the emitted coherent photons to that of the absorbed photons. Since coherent photon emission is related to the electron-photon coupling strength of the annihilated photo-excited electron-hole pairs, ${\mathcal{B}}_{coh}$ should be on the order of ${\left(\pi \alpha \right)}^2$ and insensitive to changes in N.

Next, we apply $R_N={\left|r_N\right|}^2+2\mathrm{Re}\left[r_N{z}_N^{*}\right]+{\left|z_N\right|}^2$ to the case that ${\theta }_j$ in $z_N$ is a random variable. Here, the definition of randomness is that if we take the time average regarding ${\theta }_j$, we have $\langle \mathrm{Re}\left[r_N{z}_N^{*}\right]\rangle =0$ and $\langle {\left|z_N\right|}^2\rangle ={\sum }_{j=1}^N{\left|Z_j\left(\mathcal{B}\right)\right|}^2$. We will refer to this case as incoherent corrections, which also include the cases that the global phase takes 0 and $\pi$ if there is a perturbation that can mix the two energy levels. An interference term is now included in ${\left|z_N\right|}^2$ as the last term of

$$\begin{array}{c}\hfill {\left|z_N\right|}^2=\sum _{i=1}^N{\left|Z_i\left(\mathcal{B}\right)\right|}^2+\sum _{i\ne j}{e}^{i({\theta }_i-{\theta }_j)}{Z}_i\left(\mathcal{B}\right)Z_j^{*}\left(\mathcal{B}\right),\end{array}$$

(10)

but it vanishes when taking the time average and only the first term of the incoherent corrections remains. [44] Inelastic scattering of light such as Raman scattering is usually considered to give rise to incoherent photons. Let ${\mathcal{B}}_{inc}$ be the branching ratio of the energy of the emitted incoherent photons to that of absorbed photons. Since $Z_j\left({\mathcal{B}}_{inc}\right)$ is proportional to $\sqrt{{\mathcal{B}}_{inc}{A}_j^N}$ [Equation (3)], the incoherent corrections are proportional to ${\mathcal{B}}_{inc}{A}_j^N$. For Raman scattering, the parameter ${\mathcal{B}}_{inc}$ is fundamentally determined by the electron-photon and electron-phonon coupling strengths, and it should not be so sensitive to the change in N. Indeed, the incoherent corrections with a constant ${\mathcal{B}}_{inc}$ follow the measured N dependence of the G band Raman intensity [Section 5.3]. The G band consists of optical phonons at the $\Gamma$ point, whose lattice vibrations are in-plane.

A generalized reflection formula covering the above two cases (coherent and incoherent corrections) can be written as

$$\begin{array}{c}\hfill R_N(\lambda ,\theta ,{\mathcal{B}}_{coh},{\mathcal{B}}_{inc})\equiv {\left|r_N+e^{i\theta }\sum _{j=1}^N\frac{t_N\left(j\right)+r_N\left(j\right)}{|t_N\left(j\right)+r_N\left(j\right)|}\frac{b_N{\left(j\right)}^{*}}{|b_N\left(j\right)|}{Z}_j\left({\mathcal{B}}_{coh}\right)\right|}^2+\sum _{j=1}^N{\left|Z_j\left({\mathcal{B}}_{inc}\right)\right|}^2.\end{array}$$

(11)

When ${\mathcal{B}}_{coh}={\mathcal{B}}_{inc}=0$, $R_N(\lambda ,\theta ,{\mathcal{B}}_{coh},{\mathcal{B}}_{inc})$ reduces to Equation (1) with Equation (2). Figure 7 shows the N dependence of $R_N(\lambda ,\theta ,0,0)$ (dashed), $R_N(\lambda ,\pi ,{\mathcal{B}}_{coh}=0.0007,0)$ (red), $R_N(\lambda ,0,{\mathcal{B}}_{coh}=0.0007,0)$ (blue), and ${\sum }_{j=1}^N{\left|Z_j({\mathcal{B}}_{inc}=0.1)\right|}^2$ (black), for a fixed $\lambda =540$ nm. Note that the incoherent corrections always increase the reflectance and preclude zero reflections at $N\sim 15$, which is in contrast to the coherent corrections. Moreover, the G band Raman intensity is enhanced when zero reflection occurs [10,45,46,47]. This situation is called interference-enhanced Raman scattering, [9] and it is reasonably reproduced by Equation (11). The experimental fact that zero reflection is observed at $N\sim 15$ [see 14 and 15-layer in Figure 1] shows that ${\mathcal{B}}_{coh}\sim 0.0007$, ${\mathcal{B}}_{inc}$ is much smaller than 0.1, and $\theta =\pi$ (red curve in Figure 7).

Figure 7. Reflectance of N-layer graphene on SiO₂/Si substrate. Dependence of the corrections to $R_N\left(540\mathrm{nm}\right)$ [Equation (11)] on N, where $d_{{\mathrm{SiO}}_2}=268$ nm, ${\mathcal{B}}_{inc}=0.1$, ${\mathcal{B}}_{coh}=0.0007$ and the global phase is − or +. The incoherent components are enhanced (which is called interference-enhanced Raman scattering) at around 20 layers.

Figure 7. Reflectance of N-layer graphene on SiO₂/Si substrate. Dependence of the corrections to $R_N\left(540\mathrm{nm}\right)$ [Equation (11)] on N, where $d_{{\mathrm{SiO}}_2}=268$ nm, ${\mathcal{B}}_{inc}=0.1$, ${\mathcal{B}}_{coh}=0.0007$ and the global phase is − or +. The incoherent components are enhanced (which is called interference-enhanced Raman scattering) at around 20 layers.

5. Applications of Model

In this section, we show that the discrepancy between the measured reflectance of multilayer graphene [Figure 1] and the prediction of the model [Equation (11)] is sufficiently small for the present purpose. Our model is, therefore, nearly equivalent to the Fresnel equation with the experimental optical constants (n and k), [19] while our model can describe the interesting aspects of reflection. Using monolayer graphene, we provide a detailed explanation of the mechanism modifying the reflectance without introducing any artificial change in the dynamical conductivity. To showcase the versatility of our model, we also explore Raman scattering as incoherent corrections.

5.1. Comparison of Theory and Experiment

Calculated contrasts are shown in Figure 8 as red solid curves, which include coherent corrections only (i.e., ${\mathcal{B}}_{inc}=0$). Green dashed curves represent the primary model of Equation (1), which does not include any corrections due to light emission. All calculations were performed with ${\mathcal{B}}_{coh}=0.0007$, where this value was chosen so that we could obtain good agreement between the calculations and observations for all layers. Note that this value is consistent with zero reflection being observed for $N\sim 15$.

From the consistency between the calculated and measured contrasts shown in Figure 8 (except $C_1$ and $C_3$), we can draw two main conclusions. First, the $\pi$ phase ($\theta =\pi$) of the coherent corrections is essential. If we adopt 0 phase ($\theta =0$), a serious discrepancy arises, as can be readily imagined from the relative location of the red solid curves with respect to the green dashed ones. Second, the incoherent corrections are rather small. In fact, for most of the layers examined (not shown in Figure 8), the incoherent corrections did not improve the fitting. Our estimated reasonable range of ${\mathcal{B}}_{inc}$ is less than 0.01.

Only for the 226-layer, there is a slight but non-negligible deviation of the red solid curve from the measured contrast. A relatively small difference between the red solid and green dashed curves shows that the strength of $|r_N+z_N{|}^2-{\left|r_N\right|}^2$ is suppressed and that $z_N$ is under some special phase balance by interference for $N\sim 226$. Thus, a slight shift in ${\theta }_j$ might improve the fitting. For example, second order corrections which arise due to a self-consistent calculation of $E_j^N$ (and $A_j^N$) might be relevant to this.

5.2. Monolayer

Unfortunately, the Fresnel equation is inconsistent with the measured contrast of monolayer graphene (see $C_1$ in Figure 1), probably because reflectance depends on the condition of the interface between graphene and substrates. However, as we have carefully confirmed that almost all the information of the corrections from the emitted light is included in the reflectance of monolayer graphene, we believe that the contrast $C_1$ calculated from the Fresnel equation with the experimental n and k values is the result that we should compare with the model.

In Figure 9a, we present simulated (black dotted) and calculated (green dashed and red solid) spectral contrasts of monolayer graphene on the SiO₂/Si substrate. Clearly, the corrections are of physical significance; black dots and red solid curve almost perfectly match. In $R_1={|r_1+z_1|}^2$, $r_1$ becomes a positive number only near $\lambda =2d_{{\mathrm{SiO}}_2}$ as shown in Figure 9c. This is due to the destructive interference caused by SiO₂, which also increases the absorption because $A_1^1=\pi \alpha {|1+r_1|}^2$. This enhanced absorption leads to the main difference between the reflectances from the substrate ($R_0$) and from the graphene on it ($R_1$), increasing the visibility of graphene. The correction due to light emission $z_1$ is a negative number due to the negative global phase of an energetically stable configuration of light. Thus, the corrections increase $|R_1-R_0|$. Namely, the increase in the visibility of graphene is mainly due to the substrate-induced enhancement of light absorption and is partly due to light emission.

It is important to note that the primary model can explain the measured contrasts if $\alpha$ is more than 20 percent larger than $1/137$. However, this immediately leads to an inconsistency with the experimental fact that (suspended) monolayer graphene absorbs $\sim 2.3$ percent ($=\pi \alpha$) of the incident visible light [32]. Additionally, the primary model with such corrections to $\alpha$ does not reproduce the measured contrasts of many samples with different thickness. Furthermore, the primary model is not applicable to the Raman effect, while our model can include it in a natural way, as shown below.

5.3. Raman Scattering as Incoherent Corrections

We measured the G band Raman peak intensity as a function of layer number in order to verify the interpretation of Raman scattering being an incoherent light emission [10]. The incoherent correction is defined by the last term of Equation (11) as $I_N\equiv {\sum }_{j=1}^N{\left|Z_j\left({\mathcal{B}}_{inc}\right)\right|}^2$. $I_N$ is proportional to ${\mathcal{B}}_{inc}$ and does not depend on the value of ${\mathcal{B}}_{inc}$ when scaled. As shown in Figure 10, there is a reasonable similarity between the measured Raman intensity (dots) and calculated incoherent component (dashed curve). Also plotted is an approximation of $I_N$ (dot-dashed curve) defined by

$$\begin{array}{c}\hfill I_N^{\mathrm{direct}}\equiv \frac{{\mathcal{B}}_{inc}}{2}\sum _{j=1}^N{\left|t_{j-1}^g\right|}^2{A}_j^N,\end{array}$$

(12)

to show the effect of multiple scattering of incoherent light. The approximation overestimates the intensity for thick samples, as readily imagined.

On the other hand, there is a noticeable discrepancy between them for samples with fewer than 30 layers, where there is a dip in the reflectance that is similar to the observation by No et al. [47]. The assumption of a random phase for ${\theta }_j$ in Equation (7) is a possible reason for the discrepancy, because random phases can undergo synchronization or entrainment [43]. An intermediate state of the phase ${\theta }_j$ that is neither random nor perfectly coherent may account for the behavior.

6. Discussion

There is a possibility that substrates play a decisive role in determining the selection of the two states specified by the global phase $\theta$ (0 and $\pi$). To see this, let us consider monolayer graphene suspended in the air. From Equation (1), the reflection and transmission coefficients (without substrates) are $r_1\sim -\pi \alpha /2$ and $t_1\sim 1-\pi \alpha /2$, and the absorption is given by $A_1^1=1-r_1^2-t_1^2\sim \pi \alpha$. Light emission modifies the reflection coefficient as follows:

$$\begin{array}{c}\hfill r_1+e^{i\theta }{Z}_1\left(\mathcal{B}\right)=-\frac{\pi \alpha }{2}+e^{i\theta }\sqrt{\frac{\mathcal{B}}{2}\pi \alpha }.\end{array}$$

(13)

The magnitude of the second term is 0.0028 when $\mathcal{B}=0.0007$ which is about 25 percent the magnitude of the first term (0.011). The reflection increases or decreases depending on $\theta$. When $\theta =\pi$ (0), the correction term is negative (positive) in sign; therefore, the light emission increases (decreases) the reflectance. Mathematically speaking, the change in the reflectance is equivalent to a replacement of $\alpha$ as $\alpha \to {\alpha }_{\mathrm{eff}}=\alpha -e^{i\theta }\sqrt{\pi \mathcal{B}/2\alpha }$. This immediately leads to an inconsistency with the experimental fact that suspended monolayer graphene absorbs $\sim 2.3$ percent ($=\pi \alpha$, the measured uncertainty is within 5 percent of $\pi \alpha$) of the incident visible light [32]. This inconsistency is resolved by considering that for graphene suspended in the air the states of $\theta =\pi$ and 0 are degenerate and the effect of light emission effectively disappear by interference. This contrasts with monolayer graphene on SiO₂/Si substrates for which $\theta =\pi$ is selected.

Since coherent and incoherent emissions are two extreme cases (uniform and random) of the phase ${\theta }_j$ in Equation (7), a sharp distinction between the coherent and incoherent emissions is not always possible. The proper way to calculate the reflectance is to derive a dynamical model of ${\theta }_j$ at a microscopic level, [43,48] and use it to calculate $R_N={\left|r_N+{\sum }_{j=1}^N{e}^{i{\theta }_j}{Z}_j\left(\mathcal{B}\right)\right|}^2$. Especially in the case of monolayer, they are inseparable as

$$\begin{array}{c}\hfill R_1={\left|r_1+e^{i{\theta }_1}{Z}_1\left({\mathcal{B}}_{coh}\right)+e^{i{\varphi }_1}{Z}_1\left({\mathcal{B}}_{inc}\right)\right|}^2,\end{array}$$

(14)

when ${\varphi }_1$ and ${\theta }_1$ have some correlation. Then, the reflectance is always subject to fluctuations created by the last term (through electron-phonon interactions).

Our model postulates that the effect of the emitted light appears as a correction to the reflection (and transmission) coefficient, not as a correction to the dynamical conductivity. Meanwhile, spontaneous emission is generally treated as a loss, and it is often included as a phenomenological relaxation constant in the dynamical conductivity. For example, in Ref. [19], the authors introduce relaxation constants for Drude-Lorentz oscillator model to interpret the measured optical constants. However, an excessively large relaxation constant (or very short lifetime) for the Drude term (0.6 fs) already raises concerns about the naturalness of including such a relaxation parameter [49]. The justification of our postulate needs a theoretical clarification at a more fundamental level of quantum electrodynamics, which is capable of describing photon creation and annihilation, while the excellent agreement between the model and experiments clearly shows that this postulate works well. Our model is also consistent with a theoretical result that the dynamical conductivity is free from such a correction when graphene is undoped (i.e., charge neutrality condition is satisfied) [31]. Moreover, the model explains the N-dependence of the Raman intensity as the incoherent corrections (Section 5.3), besides the contrast, in a unified manner.

7. Conclusions

In summary, we succeeded in explaining the measured visible contrasts of multilayer graphene samples on an SiO₂/Si substrate by including coherent light emissions that come from some fraction (${\mathcal{B}}_{coh}=0.0007$) of the absorbed photon energy. The coherent corrections are essential for assuring the consistency between theory and experiment, while the incoherent corrections can be neglected for the contrast.

Photo-excited electrons contribute insignificantly to light emission when they are distant from the bottom of the conduction band. Namely, the value of $\mathcal{B}$ for those electrons would be suppressed. In fact, graphene lacks a bandgap and the branching ratio of coherent light emission to absorbed light is very small (${\mathcal{B}}_{coh}=0.0007$). What we have argued for in this paper is corrections (to the main effect) that have small branching ratios. However, whether light emission from those electrons can be entirely neglected depends on various factors. Graphene serves as an interesting example where corrections are greatly enhanced by its large absorption. The destructive interference effect from the substrate and the multilayer-induced coherence are the means by which it is observable in the reflectance. Similar emission-based corrections could be anticipated for other layered materials without band gaps, and the method developed here may prove useful in accurately understanding their optical properties.

Our formulation of the reflectance using the transfer matrix method has a descriptive ability for layered materials having defects and irregularities. For this reason, and considering the success it has had in describing the reflectance of relatively simple systems (graphene multilayer at visible range), we believe that some future form of this theory may be useful in describing the optical properties of any layered material with or without a band gap.