Introducing Corrections to the Reflectance of Graphene by Light Emission

Monolayer graphene absorbs 2.3 percent of the incident visible light. This 'small' absorption has been used to emphasize the visual transparency of graphene, but it in fact means that multilayer graphene absorbs a sizable fraction of incident light, which causes non-negligible fluorescence. In this paper, we formulate the light emission properties of multilayer graphene composed of tens to hundreds of layers using a transfer matrix method and confirm the method's validity experimentally. We could quantitatively explain the measured contrasts of multilayer graphene on SiO$_2$/Si substrates and found sizable corrections, which cannot be classified as incoherent light emissions, to the reflectance of visible light. The new component originates from coherent emission caused by absorption at each graphene layer. Multilayer graphene thus functions as a partial coherent light source of various wavelengths, and it may have surface-emitting laser applications.


I. 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 SiO 2 (d SiO2 ) are known to provide advantages in this context.[4,5] Specifically, when the optical path length difference of the SiO 2 layer is one-half or three-halves of 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 2 /Si substrates can be reduced to zero for normally incident visible light with a wavelength (λ) approximately equal to 2d SiO2 .Zero reflectance is achieved through destructive interference caused by SiO 2 and a significant absorption coefficient of graphene.[6][7][8] More importantly, the reflectance of multilayer graphene on SiO 2 /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 2 /Si substrates) caused by the light it emits after absorbing incident light.The coherent corrections are dependent on only two parameters: the branching ratio (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 interfer-ence.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 Sec.II, we present fundamental insights into the reflectance of multilayer graphene on SiO 2 /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 III 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 Sec.IV.Finally, Section V offers a discussion of the findings, and Section VI provides the conclusions.

II. 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 electronlight coupling constant, extracted from the established n and k values, deviates significantly from the finestructure constant α ≃ 1/137.In fact, the difference is beyond the level of small corrections of order of α 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.

A. Fresnel equation
The reflectance of N -layer graphene on SiO 2 /Si substrate to normally incident light of wavelength λ can be formulated using the reflection coefficient r N (λ) as Here, φ ≡ n SiO2 d SiO2 2π λ is the phase acquired by light after it propagates through a distance d SiO2 in SiO 2 , and ϕ ≡ √ ε g 2π λ d N is the complex phase acquired when light passes through N -layer graphene of thickness d N ≡ N d.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 ε g = (n + ik) 2 .[19] n Si and n SiO2 are the refractive indexes of Si and SiO 2 , respectively.Si is treated as an absorbing substrate having a semi-infinite thickness (n Si is a complex number) whose dispersion is taken into account, [28] while SiO 2 is treated as an absorption-free film.[29] R N (λ) depends sensitively on the two phases, φ and ϕ.When N = 0 or ϕ = 0 in Eq. ( 1), R 0 (λ) corresponds to the reflectance of the substrate without graphene.It can be minimized for a specific λ by destructive interference; namely, R 0 (λ) is at a minimum when cos φ = 0 as . [30] Monolayer graphene is most easily detectable on SiO 2 /Si substrates when destructive interference occurs, because |R 0 (λ) − R 1 (λ)| takes a maximum when cos φ = 0 (i.e. when λ ≃ 2d Si02 because n SiO2 ≃ 1.46).[4] B. Comparison of measured and calculated contrasts Multilayer graphene was prepared by exfoliating highly oriented pyrolytic graphite (HOPG) on the same SiO 2 /Si substrate.The reflectance of the multilayer graphene was measured with a spectroscopic reflectometer (TohoS-pec3100, Toho Technology) using a ×50 objective lens.First, we determined that d SiO2 = 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 λ ≃ 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 Fig. 1 (The error bars for the data are within each circle).It is important to note that we present contrasts (C N ≡ 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 λ = 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 λ = 520 nm.
The contrasts, calculated using Eq. ( 1), are represented by the red solid curves in Fig. 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, re- Measured (black circles) and calculated (red solid) spectral contrasts of N -layer graphene on the same SiO2/Si substrate.The measurments were conducted using a white-light source at room temperature.The red solid curves were obtained from Eq. ( 1) with εg = (n+ik) 2 using optical constants obtained by El-Sayed et al. [19] The horizontal axis is λ (nm), and the reliable range of our spectrometer is 450 to 800 nm.
flectance 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 λ ≃ 2d Si02 , 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 Sec.III).To capture the essential role of graphene in achieving zero reflection, let us consider Eq. ( 1) when cos φ = 0.For the numerator to vanish, nSi ) is an approximate layer number that gives zero reflection.This argument makes it easy to understand that the dominant imaginary part of ε g is essential for zero reflection to occur.

C. 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 ε g = (n + ik) 2 values in Eq. ( 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 πα for visible light.[31] As a result, (suspended) monolayer graphene absorbs ∼ 2.3 percent (= πα) 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): σ graphite = πα/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 λ is written as where ε r is the dielectric constant of the interlayer space.By equating the right-hand side of Eq. ( 2) with (n + ik) 2 , we define an effective coupling constant α eff = 4nkd/λ, which is compared with α in Fig. 2(a).The difference between α eff and α is actually larger than the order of 10 percent of α which is beyond the level of small corrections of order of α 2 considered theoretically so far.Similarly, we define ε eff = n 2 − k 2 and plot it in Fig. 2(b).If the interlayer space is a vacuum, an appropriate choice of ε r would be 1.However, ε eff is very different from unity because the electronic wave function of the π-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 ε 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 ∼ 6 here).

III. 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 (= πα) 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 α.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.
A. Basic idea Figure 3(a) 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 Fig. 3(a) 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 Sec.III B. 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 Sec.III C. 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 We use a transfer matrix method to calculate r N and z N .[21,34] Transfer matrix method is useful in calculat- ing reflection (up arrow) and transmission (down) coefficients at each layer [r N (j) and t N (j) in Fig. 3(a)] in addition to the electric field E N j (j = 1, . . ., N ) that determines the absorption of the jth layer as The total absorption of N -layer graphene is

B. Primary model
Basically, primary model means the Fresnel equation of Eq. ( 1) where ε g is given by Eq. (2).ε g has a single unknown constant ε r which can be estimated from the following observations.Ultimately, we conclude that ε r = 4.5.
First, there must not be a large discrepancy between ε g and the experimentally determined optical constants.[19,41] Experimental values for n and k [19] are shown as • and • in Fig. 4(a).The lines depict (bare) optical constants plotted using Eq. ( 2) with the refractive index n g = Re[ √ ε g ] and absorption coefficient k ≃ k g but n has a certain difference from n g .When ε r = 5.5, n ≃ n g but k has a certain difference from k g .Because n g increases with increasing ε r while k g decreases, there is no ε 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, ε 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 2 .As shown in Fig. 4(b), the position is red-shifted by increasing N , indicating that even thin graphite samples significantly impact the light interference effect.When ε 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 ε 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, ε r = 4.5 roughly reproduces the reflectivity of graphite in the infrared region.[42] Similar ε r values have been used to reproduce the observed reflectivities of graphite and graphene.[18,19,41] We note that the value of ε r is less than the magnitude of the imaginary part of ε g , since visible light has a much longer λ (400∼800 nm) than d, although α is certainly a small quantity.The optical properties of multilayer graphene are thus characterized by the large imaginary part of ε g .

C. 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 Fig. 3(a)].The amplitude of the emitted light is assumed to be the square root of the layer absorption A N j ≡ πα|E N j | 2 , [34] multiplied by the branching ratio, B, i.e., (B/2) A N j , where 1/ √ 2 means that the light emission is directionindependent along the c-axis.Note that A N j depends not only on j and N but also on λ and d SiO2 .
To examine how light emitted from the jth layer affects the reflectance, we define two subsystems, as shown in Fig. 3(b): one is an isolated (j − 1)-layer graphene in the air; the other is (N − j)-layer graphene on SiO 2 /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 g j−1 and r g j−1 ] and the reflection coefficient of (N − j)layer graphene on the SiO 2 /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 Fig. 3(b)].These can be obtained by a self-consistent manner as follows.After calculating X (n) j , we add it to (B/2) A N j of the incident light to the (j − 1)-layer graphene (in the air) as (B/2) . Then, we add a new Y  ) strate as (B/2) . 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 B: The "corrected" electric fields at an infinitesimal distance above and below the jth layer become (B/2) , respectively.Self-consistency, whereby lim n→∞ X (n) j = X j and lim n→∞ Y (n) j = 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 By comparing this with Z j (B), we see that more accurate value of the amplitude of the emitted light is given by multiplying {• • • } (B/2) A N j with 1 + r g j−1 as the renormalization constant, and 3(a) represents.Therefore, we rede-fine Z j as b N (j) ≡ (1 + r g j−1 ) We can interpret Z j (B) as follows.The transmission coefficient t g j−1 is the direct propagation of the renormalized light emitted from jth layer to the air, and |t g j−1 | 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 {• • • }.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 (B) 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: Accordingly, the corrected reflectance is uniquely determined by The value of z N depends on these phases θ 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 (|Ψ a ⟩) that the primary model describes (see left side of Fig. 3(a)).All the information of light is expressed by the coefficients r N (j) and t N (j) (j = 1, . . ., N ).Second, we can also define another quantum state of light (|Ψ b ⟩) for the emitted light (see right side of Fig. 3(a)).All the information of emitted light is expressed by the coefficients b N (j).These two states have an overlap b * N (j)t N (j) + b * N (j)r N (j) caused by jth layer graphene.Thus, if we consider a linear superposition of these states as |Ψ a ⟩ + e iθ |Ψ b ⟩ to form energy eigenstates, the phase e iθj must be chosen so that e −iθj b * N (j)(t N (j) + r N (j)) becomes a real number, namely The factor ± is a global phase (θ) 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 twolevel state, there are two possible linear superpositions of their energy eigenstates, −1 (θ = π) or +1 (θ = 0).The minus sign (e iπ ) is assigned to the lower energy state.From Eqs. ( 5) and ( 8), we obtain e iθj Z j (B) = ± t N (j)+r N (j) |t N (j)+r N (j)| t g j−1 |b N (j)|.Including the correction due to coherent light emission leads to where 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, B coh should be on the order of (πα) 2 and insensitive to changes in N .Next, we apply 2 to the case that θ j in z N is a random variable.Here, the definition of randomness is that if we take the time average regarding θ j , we have ⟨Re We will refer to this case as incoherent corrections, which also include the cases that the global phase takes 0 and π if there is a perturbation that can mix the two energy levels.An interference term is now included in |z N | 2 as the last term of i̸ =j e i(θi−θj ) Z i (B)Z * j (B), (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 B inc be the branching ratio of the energy of the emitted incoherent photons to that of absorbed photons.Since Z j (B inc ) is proportional to B inc A N j [Eq.( 3)], the incoherent corrections are proportional to B inc A N j .For Raman scattering, the parameter 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 B inc follow the measured N dependence of the G band Raman intensity [Sec.IV C].The G band consists of optical phonons at the Γ point, whose lattice vibrations are in-plane.

IV. APPLICATIONS OF MODEL
In this section, we show that the discrepancy between the measured reflectance of multilayer graphene [Fig.1] and the prediction of the model [Eq.(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.

A. Comparison of theory and experiment
Calculated contrasts are shown in Fig. 6 as red solid curves, which include coherent corrections only (i.e., B inc = 0).Green dashed curves represent the primary model of Eq. ( 1), which does not include any corrections due to light emission.All calculations were performed with 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 ∼ 15.
From the consistency between the calculated and measured contrasts shown in Fig. 6 (except C 1 and C 3 ), we can draw two main conclusions.First, the π phase (θ = π) of the coherent corrections is essential.If we adopt 0 phase (θ = 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 Fig. 6), the incoherent corrections did not improve the fitting.Our estimated reasonable range of B inc is less than 0.01.
Only for the 226-layer, there is a slight but nonnegligible 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 − |r N | 2 is suppressed and that z N is under some special phase balance by interference for N ∼ 226.Thus, a slight shift in θ j might improve the fitting.For example, second order corrections which arise due to a self-consistent calculation of E N j (and A N j ) might be relevant to this.

B. Monolayer
Unfortunately, the Fresnel equation is inconsistent with the measured contrast of monolayer graphene (see C 1 in Fig. 1), probably because reflectance depends on the condition of the interface between graphene and sub-strates.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 Fig. 7(a), we present simulated (black dotted) and calculated (green dashed and red solid) spectral contrasts of monolayer graphene on the SiO 2 /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 λ = 2d SiO2 as shown in Fig. 7(c).This is due to the destructive interference caused by SiO 2 , which also increases the absorption because . 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 α 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 ∼ 2.3 percent (= πα) of the incident visible light.[32] Additionally, the primary model with such corrections to α 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.

C. 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 Eq. ( 11) as proportional to B inc and does not depend on the value of B inc when scaled.As shown in Fig. 8, 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 (dotdashed curve) defined by to show the effect of multiple scattering of incoherent light.The approximation overestimates the intensity for thick samples, as readily imagined.1) with Eq. ( 2)] which does not include corrections.The red solid curves include coherent corrections.The horizontal axis is λ (nm), and the reliable range of our spectrometer is 450 to 800 nm.
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 θ j in Eq. ( 7) is a possible reason for the discrepancy, because random phases can undergo synchronization or entrainment.[43] An intermediate state of the phase θ j that is neither random nor perfectly coherent may account for the behavior.

V. DISCUSSION
There is a possibility that substrates play a decisive role in determining the selection of the two states specified by the global phase θ (0 and π).To see this, let us consider monolayer graphene suspended in the air.From Eq. ( 1), the reflection and transmission coefficients (without substrates) are r 1 ∼ −πα/2 and t 1 ∼ 1 − πα/2, and the absorption is given by A Light emission modifies the reflection coefficient as follows: The magnitude of the second term is 0.0028 when B = 0.0007 which is about 25 percent the magnitude of the first term (0.011).The reflection increases or decreases depending on θ.When θ = π (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 α as α → α eff = α − e iθ πB/2α.The dots represent the measured peak intensity of the Raman G band when using a light source with a wavelength of 532 nm.The dashed curve corresponds to the incoherent correction, while the dot-dashed curve is an approximation that ignores multiple scattering.
immediately leads to an inconsistency with the experimental fact that suspended monolayer graphene absorbs ∼ 2.3 percent (= πα, the measured uncertainty is within 5 percent of πα) of the incident visible light.[32] This inconsistency is resolved by considering that for graphene suspended in the air the states of θ = π and 0 are degenerate and the effect of light emission effectively disappear by interference.This contrasts with monolayer graphene on SiO 2 /Si substrates for which θ = π is selected.Since coherent and incoherent emissions are two extreme cases (uniform and random) of the phase θ j in Eq. ( 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 θ j at a microscopic level, [43,48] and use it to Especially in the case of monolayer, they are inseparable as when ϕ 1 and θ 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 (Sec IV C), besides the contrast, in a unified manner.

VI. CONCLUSION
In summary, we succeeded in explaining the measured visible contrasts of multilayer graphene samples on an SiO 2 /Si substrate by including coherent light emissions that come from some fraction (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 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 (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 multilayerinduced 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.] and gray [run 2 (run 3)] dots) and calculated (green and blue dashed curves) spectral reflectivities of N -layer graphene on the same SiO2/Si substrate.A white-light source was used, and the reflectivity was measured at room temperature.The horizontal axis is λ (nm), and the reliable range of our spectrometer is 450 to 800 nm.
FIG.1.Measured (black circles) and calculated (red solid) spectral contrasts of N -layer graphene on the same SiO2/Si substrate.The measurments were conducted using a white-light source at room temperature.The red solid curves were obtained from Eq. (1) with εg = (n+ik) 2 using optical constants obtained by El-Sayed et al.[19]The horizontal axis is λ (nm), and the reliable range of our spectrometer is 450 to 800 nm.

FIG. 3 .
FIG. 3. 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 Zj is explained in detail in the text.

5 and 5 . 5 .
When ε r = 4.5, 2) A N j of the light incident to the (N − j)-layer graphene on the SiO 2 /Si sub-

FIG. 4 .
FIG. 4. ε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 εr = 4.5 is a reasonable value.

FIG. 6 .
FIG.6.Measured (black circles) and calculated (green dashed and red solid) spectral contrasts of N -layer graphene on the same SiO2/Si substrate.The green dashed curves represent the primary model [Eq.(1) with Eq. (2)] which does not include corrections.The red solid curves include coherent corrections.The horizontal axis is λ (nm), and the reliable range of our spectrometer is 450 to 800 nm.

FIG. 8 .
FIG. 8. Raman intensity as incoherent corrections.The dots represent the measured peak intensity of the Raman G band when using a light source with a wavelength of 532 nm.The dashed curve corresponds to the incoherent correction, while the dot-dashed curve is an approximation that ignores multiple scattering.

FIG. 9 .
FIG.9.Determination of SiO2 thickness.Measured (black [run 1] and gray [run 2 (run 3)] dots) and calculated (green and blue dashed curves) spectral reflectivities of N -layer graphene on the same SiO2/Si substrate.A white-light source was used, and the reflectivity was measured at room temperature.The horizontal axis is λ (nm), and the reliable range of our spectrometer is 450 to 800 nm.
Contrasts of monolayer graphene on SiO2/Si substrates.(a)Onecontrast is calculated using Eq.(1) with the experimental εg (black dotted), while the other two curves are obtained using the model with corrections (red solid curve) and without corrections (green dashed curve).The increased visibility of graphene is primarily attributed to substrate-induced enhancement of light absorption and partly to light emission.(b) Layer absorption of monolayer graphene on SiO2/Si substrates A 1 1 = πα|1+r1| 2 .The destructive interference effect of the substrate enhances absorption near λ ≃ 2dSiO 2 , where r1 is a positive number (c).