1. Introduction
Silicon photonics based on the mature silicon-on-insulator (SOI) platform has aroused tremendous interest for building compact, high performance, and energy-efficient photonic integrated circuits (PICs), which has also promoted the rapid development of optical interconnects, especially for data center and telecom applications [
1,
2,
3,
4]. To keep up with the ever-increasing demand of high capacity on-chip optical interconnects, much attention has been focused on the development of advanced multiplexing technologies (e.g., wavelength-division-multiplexing [
5], polarization-division-multiplexing [
6], mode-division-multiplexing [
7]) and their hybrid multiplexing types [
8]. However, these multiplexing technologies will inevitably impose higher performance requirements on optical modulators, e.g., higher modulation efficiency (ME), larger optical bandwidth, ultra-fast modulation speed, ultracompact size, lower energy consumption, and complementary metal-oxide-semiconductor (CMOS) processing compatibility [
1]. These requirements cannot be well satisfied by silicon modulators which are based on weak plasma dispersion effects [
9]. As a result, new materials, structures, and mechanisms have been introduced to the SOI platform to improve the modulation performance of silicon modulators. For instance, III–V epitaxial materials were transferred to an SOI wafer forming a hybrid silicon electro-absorption modulator (EAM) [
10], which achieves a high ME of 0.25 dB/µm and over 67 GHz modulation bandwidth at the wavelength of 1.3 µm, while at the cost of incompatible CMOS processes and complex fabrication procedures. The silicon-organic hybrid (SOH) platform—filling nano-slot waveguides with organic cladding materials—was developed, and the corresponding SOH-based electro-optic modulator consumes only 0.7 fJ/bit for the device length of 1 mm with 12.5 GHz modulation bandwidth [
11]. The effort to decrease the device footprint is likely to be in progress. Other waveguide structures could also be applied onto this platform, e.g., slot photonic crystal waveguide [
12]. Moreover, by means of the large thermo-optic coefficient of silicon [
13], the thermal modulation method has also been explored for silicon modulators, but the modulation speed was limited and the energy consumption was relatively high. Therefore, new hybrid silicon modulators with better performance and smaller size are highly desirable to fulfill the modulation requirements of high capacity on-chip optical interconnects [
14].
Recently, graphene, a single-atom-layer of graphite with its atoms arranged in a hexagonal lattice [
15], has become a promising two-dimensional (2D) material in the field of photonics and electronics with its exceptional optical/electrical properties such as ultra-high carrier mobility (e.g., >200,000 cm
2·V
−1·s
−1 at room temperature) [
16], tunable conductivity, broad bandwidth, zero bandgap, and CMOS compatible processes [
17,
18,
19]. These remarkable properties match up well with the performance requirements and developing directions of silicon modulators if graphene can be effectively integrated on the SOI platform [
20,
21]. In 2011, the first single-layer graphene-on-silicon EAM was reported with an ME of 0.1 dB/µm and modulation bandwidth of 1.2 GHz [
22]. In 2012, a double-layer version of graphene-on-silicon EAM was also developed with an ME of ~0.16 dB/µm and modulation bandwidth of 1 GHz [
23]. Notwithstanding, the reported ME are relatively low because of the ultra-thin thickness of graphene (~0.34 nm) and the weak evanescent wave in the cladding regions of silicon waveguide interaction with the graphene layers coated, resulting in reduced light–graphene interaction (LGI). In order to effectively improve LGI, special waveguide structures with enhanced light intensities have been exploited to combine with graphene flake, such as slot waveguides [
24], surface plasmon polariton waveguides [
25], and hybrid plasmonic waveguides [
26], instead of the commonly used rid/ridge waveguides. The device ME is significantly increased (e.g., 0.84 dB/µm) [
27], thus the device length can be shortened and the energy consumption reduced. However, these devices are intrinsically polarization-dependent, particularly for some strong light-confinement structures, operating under only a single-polarization state. As a result, these devices would not satisfy the requirements of dual-polarization operation for the polarization-division-multiplexing [
6] and multi-dimensional hybrid multiplexing [
8] systems. It is therefore necessary to develop polarization-independent or polarization-insensitive high performance graphene-based optical modulators. By transferring graphene layers on an isosceles trapezoid silicon nanowire and capping with another silicon nanowire, a polarization-independent EAM has been proposed, with an extinction ratio (ER) higher than 20 dB in a length of 30 µm [
28]. However, its ER discrepancy for the two polarizations is high (~2.3 dB) and the precise control of the slanted angle of the silicon nanowire might pose challenges for practical fabrication processes. Furthermore, based on a typical double-layer graphene-on-silicon structure [
23], polarization-insensitive modulation can also be realized by optimizing the dimensions of the bottom silicon nanowire, while its obtained performance is limited with an ME of only 0.29 dB/µm and a transmission loss of ~2 dB [
29]. Therefore, in order to effectively enhance the performance of graphene-based polarization-insensitive optical modulators, new structures and operating mechanisms remain to be developed.
In this paper, we propose a highly-efficient, ultra-broadband, and polarization-insensitive EAM based on a multi-deposited and multi-transferred hybrid graphene-silicon waveguide structure. The input and output waveguides are conventional silicon nanowires. In order to enhance the ME or LGI for both polarizations, the proposed hybrid graphene-silicon waveguide was designed such that the light interaction with graphene layers happens in the intense field distribution region rather than in the cladding region with weak evanescent waves as in previously reported work. Meanwhile, hexagonal boron nitride (hBN) was used as the spacer or insulating layer to encapsulate graphene layers and maintain the high carrier mobility of graphene [
30]. By optimizing the dimensions of all hybrid graphene-silicon waveguide layers, polarization-insensitive modulation was realized with an ME of ~1.11 dB/µm and an ME discrepancy of lower than 0.006 dB/µm for both polarizations, which are better than previous reports [
28,
29]. We then propose a hybrid graphene-silicon-based polarization-insensitive EAM with a modulation length of 20 µm. The modulation depth (MD) and insertion loss (IL) were higher than 22 dB and lower than 0.23 dB, respectively, at 1.55 µm. The corresponding available bandwidth reached over 300 nm (from 1.367 to 1.668 µm, covering the whole S to L bands and most of the E and U bands) with a MD greater than 20 dB and MD discrepancy less than 2 dB, simultaneously. We also study the electrical properties such as 3 dB modulation bandwidth and energy consumption. Finally, we discuss some feasible ways to further enhance the device performance. The rest of this paper is organized as follows:
Section 2 presents the device design and operating principle of the proposed hybrid graphene-silicon based EAM;
Section 3 analyzes the optical modulation properties of the hybrid graphene-silicon waveguide structure;
Section 4 gives the numerical results and discusses the proposed device; and the conclusions are drawn in
Section 5.
3. Optical Properties of Hybrid Graphene-Silicon Waveguide
The discovery of graphene has accelerated the development of on-chip photonic devices to be employed as high-performance building blocks for optical interconnects. Graphene can be analyzed theoretically using either the isotropic [
32,
33] or anisotropic [
34,
35] models. In the isotropic model, ultra-high ME can be easily obtained (e.g., a record high ME of ~4.5 dB/µm [
33]) by insetting graphene layers into the slot region of the horizontal slot waveguide where its light field is greatly confined in the graphene layer region by the epsilon-near-zero effect [
33]. However, such theoretical predictions cannot be achieved in experiments. More research revealed that we ought to treat graphene as an anisotropic material since it has only one single-atom layer; the carriers are strongly confined to move in this layer [
36,
37]. Therefore, only the in-plane permittivity and conductivity of graphene can be adjusted by changing its Fermi level (or chemical potential
µc). The corresponding out-of-plane permittivity remains unchanged at around 1 or 2.5 [
37]. For the in-plane conductivity of graphene, an analytic expression can be deduced from the Kubo formula, which comprises intra- and inter-band contributions as follows [
38]:
and
where
e is the electron charge,
ħ is the reduced Planck’s constant,
T is the temperature,
kB is the Boltzmann constant,
τ is the momentum relaxation time,
ω is the angular frequency,
ε is the energy, and
is the Fermi-Dirac distribution. If
is satisfied, the inter-band conductivity can be estimated as [
38]:
Here, we have set
τ = 0.5 ps (based on a conservatively estimated carrier mobility of 10,000 cm
2·V
−1·s
−1) and
T = 300 K. Moreover, thin 3D volumes and 2D sheets are the commonly assumed models in simulation because considering the graphene flake as one-atom layer thick and modeling the graphene layer as a 2D sheet gives better consistency by comparing the experimental and simulation results [
39]. Therefore, in the following analyses, we will model the graphene layers as typical anisotropic 2D sheets and combine them with the multilayer silicon nanowires to form the hybrid graphene-silicon waveguide.
Figure 2a,b shows the effective modal index (real [Re(
neff)], imaginary [Im(
neff)]) of the hybrid graphene-silicon waveguide as a function of the graphene chemical potential (
µc) for both transverse electric (TE) and transverse magnetic (TM) modes at
λ = 1.55 µm where the modal calculation is based on the above-mentioned surface conductivity model of graphene. The calculation engine is available from a commercial software package [
40]. The widths and thicknesses of five inverted U-shaped graphene-silicon layers are
w1 = 80 nm,
w2 = 210 nm,
w3 = 340 nm,
w4 = 470 nm,
w5 = 600 nm, and
h1 = 40 nm,
h2 = 95 nm,
h3 = 150 nm,
h4 = 205 nm,
h5 = 260 nm, respectively. We note that when
µc increases, Re(
neff) increases slightly first and then decreases quickly. The maximum value is located at
µc = 0.4 eV, while Im(
neff) reveals a large change around
µc = 0.4 eV. The reason is that if
µc is under the half of the photon energy (
hv/2 = 0.4 eV), inter-band transitions of graphene layers will dominate and generate a strong optical absorption, resulting in the high imaginary part Im(
neff). If
µc >
hv/2, inter-band transitions will be blocked owing to the Pauli blocking mechanism, and the corresponding graphene layers will become transparent [
41,
42]. To quantify this absorption manner, we used the parameter of MPA, which is derived from the modal effective index [
27]:
Figure 2a,b shows that MPA presents a huge variation at around
µc = 0.4 eV and gradually remains constant when
µc deviates from 0.4 eV. To render the proposed EAM to have a relatively large tolerance for
µc, we set
µc = 0.1 eV and
µc = 0.7 eV as the “OFF” and “ON” states, corresponding to the MPA of ~1.1223 dB/µm (1.1167 dB/µm) and ~0.0033 dB/µm (0.0032 dB/µm) for the TE (TM) mode, respectively.
Figure 2a,b also depicts the electric field profiles at the “ON” state for both polarizations, while those at the “OFF” state are not shown, since their field profiles are very similar. Considering the anisotropic property of the ultra-thin graphene layer, only the graphene layers located in the horizontal plane can contribute to the EA modulation under the TE polarization, while the slot effect generated by different refractive indices between adjacent materials (hBN and Si) is on the vertical plane. Thus, the vertical slot effect could hardly affect the horizontal in-plane modulation performance due to their orthogonal polarization states. Similar behaviors can also be observed for the TM polarization. By applying suitable voltages on the graphene layers through metal contacts on both sides of the central hybrid graphene-silicon waveguide, the
µc of the graphene will vary correspondingly, and the relation between them can be calculated via the parallel plate capacitor model as follows [
43]:
where
VF is the Fermi velocity (≈10
6 m/s),
ε0 is the permittivity of vacuum,
εd and
hd are the relative dielectric constant and thickness of the spacer layer between two graphene layers, respectively, and
V0 denotes the bias voltage generated by the natural doping. Here, we treat |
Vg −
V0| as the applied voltage for simplicity, and |
Vg −
V0| versus
µc for different spacer and insulating materials (Al
2O
3, Si
3N
4, hBN) is shown in
Figure 2c where the spacer or insulating thickness
hd is 6 nm. The inset calculates the applied voltages at the “OFF” and “ON” states for different spacer materials. Among these typical spacer and insulating materials, hBN is the best choice, since hBN encapsulation of graphene layers can maintain graphene’s optical and electrical properties, e.g., a high carrier mobility (~140,000 cm
2·V
−1·s
−1) [
30], which is important to realize high-speed EAM. Besides, with the relatively low dielectric constant of hBN (~3.9), the electrical resistance-capacitance (
RC) constant is reduced, which increases the 3 dB modulation bandwidth of the hybrid graphene-silicon-based polarization-insensitive EAM and reduces its energy consumption to some extent [
44].
4. Results and Discussion
To realize polarization-insensitive modulation with high ME and ultra-broad bandwidth, we designed a hybrid graphene-silicon waveguide using multi-deposited and multi-transferred methods, which can make the light interaction with graphene layers in its waveguide interior instead of in the cladding region with weak evanescent waves. Therefore, the dimensions of each of the inverted U-shaped graphene-silicon waveguide layer and the number of layers (
N) are the key parameters. Here, the width-thickness ratio of every inverted U-shaped graphene-silicon waveguide layer was kept the same with that of the input/output silicon nanowire (
W = 600 nm,
H = 300 nm), and the width difference (Δ
w) between adjacent layers was also kept the same for simplicity.
Figure 3 shows the calculated MEs and their absolute difference for both polarizations as a function of the innermost waveguide width
w1 of the hybrid graphene-silicon waveguide with one to eight layers, where a width relation of Δ
w·(
N − 1) +
w1 =
W applies based on the assumptions above, and ME = MPA (
µc = 0.1 eV) − MPA (
µc = 0.7 eV) [
29]. We note that when the layer number
N increases, ME gradually increases for both polarizations. On the other hand, when
N increases, the minimum value of the absolute ME difference between the TE and TM modes decreases first and then increases, with the lowest value attained at
N = 5 and
w1 = 80 nm, which corresponds to the lowest polarization-dependence of the ME (<0.006 dB/µm). The obtained ME is relatively high for both polarizations (~1.11 dB/µm). When
N > 5, the required fabrication process will become increasingly complex and the polarization-dependence of ME obviously increases when compared with that of
N < 5. Therefore,
N = 5 is our optimal choice, and the corresponding
w1 and Δ
w are 80 nm and 130 nm, respectively. Other waveguide widths (
w2~
w5) and thicknesses (
h1~
h5) were obtained using the relationship mentioned above.
In the following, by leveraging the zero bandgap (or broadband absorption) property of graphene [
18,
19], we improve the optical bandwidth of the proposed EAM to cover the main optical communication bandwidth.
Figure 4a shows the variation of MPA at both the “ON” and “OFF” states for the TE and TM modes in the wavelength range from 1.2 to 1.9 µm based on the modal analysis with material dispersion considered [
45]. At the “ON” state, both polarization modes show similar wavelength dependence. At the “OFF” state, the wavelength dependences of the two polarization states behave differently. When the wavelength increases, the MPA of TE mode gradually increases, but that of TM mode first increases and then decreases. To better illustrate this wavelength-dependence of the modulation feature,
Figure 4b shows the MEs and their absolute difference for both polarizations. The wavelength ranges from O to U band are also marked. From
Figure 4b, the MEs for both polarizations are almost equal in the S and C bands, corresponding to polarization-independent operation. If the ME discrepancy for both polarizations is set at <0.1 dB/µm, the available bandwidth is from 1.2 to 1.68 µm, which covers the full band of optical communications, giving rise to broadband EA modulation. We then studied the modulation performance of the proposed hybrid graphene-silicon waveguide as a function of the spacer layer (hBN) thickness
h for both polarizations.
Figure 5 shows the calculated MEs and their absolute difference. From
Figure 5, the ME curves of the TE and TM modes exhibit an opposite variation trend. The two curves intercept near
h = 10 nm, which corresponds to the polarization-insensitive modulation.
Based on the obtained high ME, ultra-broadband, and polarization-insensitive properties, we further propose a hybrid graphene-silicon-based polarization-insensitive EAM with a modulation length of only
L = 20 µm. We analyzed its transmission properties by using a three dimensional-finite difference time domain (3D-FDTD) method [
40].
Figure 6a shows the transmission spectra of the proposed EAM at the “ON” and “OFF” states for both polarizations, where the simulation wavelength spectrum is from 1.2 to 1.9 µm and the material dispersion is included based on Palik’s book [
45]. It is noteworthy that the difference of optical transmittances between the “ON” and “OFF” states is very close for the TE and TM modes within the full band of optical communications, especially for the S and C bands. To better show this phenomenon, we also plotted MDs and their absolute difference for both polarizations as a function of the same wavelength, as shown in
Figure 6b, where MD = transmittance (
µc = 0.7 eV) − transmittance (
µc = 0.1 eV) [
29]. The wavelength-dependence of the MD curves is very similar to that of the ME curves for both polarizations shown in
Figure 4b based on the modal analysis, which also demonstrates the reliability of the 3D-FDTD method. If we use the criterion of the MD discrepancy being lower than 2 dB, the allowable bandwidth is from 1.221 to 1.668 µm, which is slightly smaller than that based on the result of modal analysis, since the input/output silicon nanowire will incur weak coupling loss in the transmission. Moreover, to make the proposed hybrid graphene-silicon-based EAM operate at high MD, we set another criterion, i.e., MD > 20 dB, and the corresponding allowable bandwidth is from 1.367 to 1.771 µm. By considering these two criteria, the final available bandwidth is 301 nm (from 1.367 to 1.668 µm), covering the whole S, C, L bands and the most of E and U bands, where ILs and back reflection losses of both polarizations are lower than 0.4 dB and −27.8 dB, respectively, within this range.
To fabricate the proposed EAM, an SOI wafer with a 40-nm thick top silicon layer and a 2 µm thick buried oxide layer was employed. The detailed fabrication process can be found in
Figure 7. For the first hybrid graphene-silicon waveguide layer, the bottom silicon nanowire with a width of
w1 = 80 nm and a length of
L = 20 µm was patterned with an E-beam lithography and etched to the buried oxide layer with an inductively coupled plasma reactive ion etching process (step 1 and step 2). Then, two graphene layers grown by CVD were transferred onto the surface of the bottom silicon nanowire, including the spacer layer of hBN [
46,
47], where the hBN encapsulation of two graphene layers (hBN-graphene-hBN-graphene-hBN heterostructure) was fabricated by employing multi-level stacking of the 2D material with a van der Waals assembly process (step 3) [
30,
48]. The thickness of the spacer layer between the two graphene layers was 6 nm. In this process, the spacer layer including its encapsulated graphene layers was extended on both sides of the central hybrid graphene-silicon waveguide to ensure effective metal contacts (step 4) [
49,
50]. The lateral distance between the metal contact and the central waveguide was set as 1.5 μm on each side to prevent metal absorption. For the second hybrid graphene-silicon waveguide layer, the second silicon nanowire was grown on the first hybrid graphene-silicon waveguide layer by epitaxy technology and then patterned and etched to form the required waveguide dimension (
w2 = 210 nm,
h2 = 95 nm, and
L = 20 µm) (step 5). The procedures of fabrication and transferring of graphene layers and spacer layers are similar with those described in the first hybrid graphene-silicon waveguide layer, including the metal contacts on both sides (step 6). We then repeated this process three more times to fabricate the third to the fifth hybrid graphene-silicon waveguide layers consecutively (step 7 to step 9). The corresponding width and thickness of the silicon nanowire were
w3 = 340 nm,
w4 = 470 nm,
w5 = 600 nm, and
h3 = 150 nm,
h4 = 205 nm,
h5 = 260 nm, respectively. The length was
L = 20 µm. Finally, the input and output silicon nanowires (
W = 600 nm,
H = 300 nm, length = 5 µm) were also fabricated on the same wafer to be butt-coupled with the hybrid graphene-silicon based EAM, allowing for easy and direct connection with other integrated components. Finally, by applying switching voltages on the graphene layers using fabricated metal electrodes, the input light could be modulated by this highly-efficient, ultra-broadband, and polarization-insensitive EAM.
Next, to investigate the performance deterioration brought by fabrication imperfections in practice, we analyzed the fabrication tolerances of the width
W, thickness
H, and length
L of the hybrid graphene-silicon-based EAM. The deviations of
W and
H were assumed to be evenly distributed in every hybrid graphene-silicon waveguide layer based on the same fabrication processing flow.
Figure 8 shows the transmittance and MD of the devices as its
W,
H, and
L deviate from their designed values (
W = 600 nm,
H = 300 nm,
L = 20 µm) at the “ON” and “OFF” states for both polarizations. When
W and
H increase, MD gradually decreases for both polarizations owing to the increase in mode field area and weakened LGI, except for the TM mode in the range of Δ
H = [−50 nm, −10 nm] where the MD of the TM mode has a peak value at Δ
H = −20 nm. Moreover, the polarization-insensitive property of the proposed EAM remains unaffected if ΔW and ΔH are kept within the ranges of [0, 100 nm] and [0, 50 nm], respectively, which gives a relatively large tolerance for fabrication. For the device length
L, the variations of transmittance and MD of the TE and TM modes are nearly the same, showing excellent polarization-insensitive properties. In addition, a higher MD for this device can be achieved only by increasing the device length because its insertion loss is very low and increases slowly owing to the small absorption loss at the “ON” state for both polarizations. Finally,
Figure 9 plots the evolution of the optical field intensity for both polarizations along the optimized hybrid graphene-silicon-based polarization-insensitive EAM at both the “ON” and “OFF” states. The structural dimensions are based on the optimized values above, and the operating wavelength is 1.55 µm. When the EAM operates at the “ON” state, the input TE and TM modes will pass through the modulation region with very low loss. In contrast, when the EAM operates at the “OFF” state, the input TE and TM modes will attenuate quickly as they enter into the modulator, which is owed to the high absorption loss of graphene layers.
Besides the above-analyzed optical properties of the proposed hybrid graphene-silicon-based polarization-insensitive EAM, its corresponding electrical characteristics are also critical, such as the 3 dB modulation bandwidth (
f3dB) and the energy consumption (
E), which are associated with the total capacitance (
C) and resistance (
R) of the optical device,
f3dB = 1/(2
πRC),
E =
CV2/4 (
V stands for the switching voltage between “OFF” and “ON” states) [
40,
42]. Here, we estimated the electrical characteristics of the proposed EAM using the parallel plate capacitor and the equivalent electrical circuit model for simplicity [
27,
28,
29]. By considering the multilayer hybrid graphene-silicon waveguide structure and the overlapping region between graphene layers, the total resistance of the device is estimated to be ~80 Ω (including the graphene sheet resistance and the metal-graphene contact resistance) and the total capacitance (
C) is estimated to be ~326 fF. Thus,
f3dB and
E are derived to be ~6.1 GHz and ~7.8 pJ/bit, respectively. To further enhance the electrical properties of the proposed EAM, its total resistance and capacitance should be effectively reduced, especially for the capacitance, which can be realized by cutting down the overlapping area between two corresponding graphene layers, choosing the spacer/insulating layer with the lower dielectric constant, and increasing the thickness of the spacer/insulating layer. In addition, traveling-wave electrode designs could also be leveraged to enhance the modulation performance of the proposed EAM [
51].