1. Introduction
Graphene is the first and most mature of bidimensional materials that has been isolated [
1,
2], and it has attracted a lot of attention from the scientific community due to its unique physical and optoelectronic properties [
3]. In particular, in addition to its saturable absorption [
4] and tunable electro-absorption [
5], several studies presented very promising results concerning its Kerr nonlinear optical response at telecom wavelengths (i.e., around 1.5 μm) both theoretically [
6,
7,
8,
9] and experimentally [
10,
11,
12]. These findings potentially make graphene a good candidate to improve the performance of nonlinear photonic devices for datacom applications. The nonlinear efficiency of a standard dielectric waveguide is typically quantified by the nonlinear parameter
, which roughly gives the nonlinear effect induced per unit of waveguide length and power. By using various types of hybrid graphene/dielectric waveguides, very high effective values of
have been measured in the literature, by means of self-phase modulation (SPM) [
12,
13] or four-wave mixing (FWM) [
10,
14,
15,
16] experiments. Typical values for
that range from a few hundred up to a few thousand per watts per meter have been reported around
μm, which is more than one order of magnitude the value achieved with tightly confining silicon waveguides (for instance, [
17]).
However, a great disparity exists in the literature regarding the nonlinear response of graphene [
18], and the abovementioned results have not led to nonlinear graphene-based hybrid devices with outstanding performance so far [
19]. This is partly due to the typically high linear absorption of graphene that tends to mitigate the net nonlinear response of the device, especially at telecom wavelengths. In addition, the origin of the nonlinear effects in graphene, which are mediated by photogenerated carriers at near-IR wavelengths, as was convincingly highlighted in recent papers [
20,
21,
22], seems to restrict the use of the graphene-effective nonlinearity to relatively low power levels, thereby restricting the absolute magnitude of these effects. Practically, if one is to use graphene (or any other 2D material) for locally enhancing the nonlinear response of an otherwise passive photonic integrated circuit, one has to better assess and understand the tradeoff associated with the loss penalty and nonlinear contribution from the resulting hybrid graphene/dielectric section integrated within such circuits.
Although Si is a mature photonic platform, which has been used with graphene in initial demonstrations [
11,
12,
13], it suffers, at telecom wavelengths, from relatively high two-photon absorption, and an associated free carrier penalty which severely limits the resulting nonlinear device performance under increasing powers. Si
N
on insulator represents another mature platform for creating low loss waveguides with a lower nonlinear response and no two-photon absorption, thereby making it particularly attractive for use with 2D material patches that can locally enhance the relatively modest nonlinear response of the Si
N
waveguides (
γ∼1 m
−1W
−1).
In this paper, we investigate the use of short graphene/Si
N
hybrid waveguide sections, and aim to assess the potential of this approach as a way to control and locally enhance the nonlinear response of a mature and low loss Si
N
waveguide circuit. More specifically, we conduct degenerate four-wave mixing measurements at
on Si
N
waveguides partially covered with millimeter-scale-long graphene patches to probe the nonlinear response of these waveguides at telecom wavelengths. The specific use of pulsed pump and probe signals with 7
peak power levels allows us to boost, a priori, the nonlinear effects with respect to prior CW four-wave mixing measurements that have been performed on waveguides covered by graphene under tens of milliwatt powers [
10,
16]. Our four-wave mixing measurements, using different graphene patch lengths and different pump powers, are compared with simulations that show that graphene locally changes the sign of the nonlinear
parameter of the waveguide, and enhances its magnitude from 1 m
−1W
−1 to −10 m
−1W
−1. Our simulations allow us to take into account the global response of the waveguide as well as to differentiate between the linear loss penalty and the nonlinear benefit induced by graphene on the whole waveguide structure. Surprisingly, the outstanding nonlinear response (
10
−13 m
2W
−1 and up to
10
−12 m
2W
−1) of graphene that was reported in a few papers [
12,
16,
23] does not translate, in our measurements, into a tremendously high local nonlinear response of our hybrid graphene/Si
N
waveguide. The modest nonlinear enhancement afforded by our hybrid graphene/Si
N
waveguides does not fully compensate for linear loss penalty induced by graphene absorption along the hybrid waveguide patch. Most critically, our work contributes to assessing the nonlinear performance and limits of waveguides locally covered with graphene patches within passive photonic circuits at telecom wavelengths, and their potential in applications.
3. Nonlinear Measurements on the Hybrid Waveguides
In order to characterize the nonlinear response of the graphene-covered section of the waveguide, we probe the idler generation as a function of both the peak pump power and the length of graphene covering the waveguide.
Figure 3 shows the spectrum measured for the waveguides covered by different graphene lengths, and for the case where the coupled peak power
.
The idler signal is clearly detected at approximately 1541
. As the graphene length in contact with the Si
N
waveguide increases, we observe an overall reduction of the whole signal (probe, idler, and pump). This signature is directly correlated to the high linear losses induced by the graphene-covered section of the waveguide as was measured in the
Section 2.1. The four-wave mixing conversion efficiency (
CE’), after propagation along a waveguide of length
L, is usually defined as [
28]
where the input probe power at the entrance of the waveguide is typically considered. In our case, each waveguide has a different drop in transmission, depending on the graphene length covering it (between ∼2 dB for
and ∼9 dB for
), which is relatively constant with wavelength across the C-band and therefore equally affects the pump, probe, and idler signals. Therefore, to somewhat separate this linear loss penalty from the FWM conversion efficiency and quantitatively analyse the impact of the nonlinear response of graphene on the idler generation, we use instead the following ratio, referred to as the FWM conversion efficiency in the rest of the paper:
This expression gives higher values than Equation (
1) for the
CE because the probe power at the output of the waveguides is decreased with respect to that at the entrance of the waveguides by the propagation loss. Equation (
2) thus allows us to leave aside the graphene-induced optical power drop equally affecting the output signals from the
CE estimation. Note that this ratio can also be directly extracted from the measured FWM spectra.
Figure 4 shows the conversion efficiency given by Equation (
2) as a function of the graphene length interacting with the waveguides, and for different coupled pump powers. We first observe that the
CE increases with the coupled pump power for each waveguide (
Figure 4a), as expected, and that this increase can be well fitted with a quadratic behavior (dotted line on
Figure 4a). Regarding the impact of graphene, we observe on
Figure 4b that for each pump power, the
CE decreases with the greater length of graphene. Furthermore, the maximum drop of
CE, measured by comparing
with the reference waveguide, increases in amplitude with power and varies from
down to
, for
and
pump power, respectively.
At first sight, various phenomena might explain the unexpected and apparent
CE reduction caused by graphene in these experiments. The first and simpler explanation is the high linear propagation loss caused by graphene, which might hide a potential increase of the nonlinear Kerr response along the hybrid graphene/Si
N
section compared to the bare Si
N
. Another explanation might be found in the sign of the nonlinear contribution of graphene, which was measured to be negative [
12,
23], i.e., of opposite sign to the Si
N
waveguide nonlinear response before and after the hybrid graphene/Si
N
section. The fact that the whole chip consists of three subsequent and distinct waveguide sections, respectively without/with/without graphene, indeed makes it more difficult to directly account for both the linear and nonlinear local contribution of graphene to the cumulative four-wave mixing response measured across the entire chip. Therefore, to tell these different effects apart, we carry out some simulations in the next section, which take into account the response of the hybrid graphene/Si
N
waveguides, and that of the bare Si
N
waveguides before and after the graphene-covered waveguide section. These simulations will thus clarify the impact of the linear loss and nonlinear response of graphene as induced locally within the hybrid graphene/Si
N
waveguide section on the response of the whole structure.
4. Comparison with Simulations and Discussion
Our waveguides are composed of three consecutive sections of waveguides (clad Si
N
; unclad graphene covered Si
N
; and clad Si
N
). We model the degenerate four-wave mixing response of each waveguide section by a coupled system of nonlinear Schrödinger equations involving the interacting pump, probe, and idler signals. The system of Equation (
3) describes the evolution of the electric field envelope,
, with
for the pump, probe, and idler, respectively [
28]. We have
and
represent the temporal and spatial derivative, respectively.
i is the imaginary unit,
represents the dispersion coefficient of the second order,
the linear propagation loss (kept constant here for all three signals), and
the nonlinear parameter of the considered waveguide section.
is associated with the probe or idler walk-off, in the frame of the pump. Graphene saturable absorption was not included in Equation (
3), because we did not observe any signature of this effect under our experimental conditions.
We use the split-step Fourier method (SSFM) to numerically solve this system of equations. It allows us to change the parameters along the propagation direction according to the specific response of the local waveguide structure, while feeding the output signal of one given section as the input for the simulation of the subsequent one.
Table 1 contains the parameters used in the simulations to model each waveguide section (with and without graphene). The second-order dispersion, effective area, and nonlinear parameter were computed by using mode profile simulation with the software Lumerical
®.
The loss at the etched interface corresponds to the estimated loss when the guided mode crosses the boundary between the clad section of the Si
N
waveguide and the unclad one (i.e., covered with graphene). By using these parameters, we run simulations in order to identify which value for the nonlinear coefficient
, associated with the hybrid graphene/Si
N
section of the waveguide, best reproduces the four-wave mixing measurements for all graphene lengths and pump power values. The simulation results are shown in
Figure 5 along with the measured
CE as a function of the graphene patch length, and for different coupled pump powers.
From the simulated curves, the case (yellow curves) is equivalent to considering that graphene has no particular nonlinear contribution to the whole waveguide response, but only affects it negatively via adding some linear loss along the hybrid graphene/SiN section. We obtain from these curves the direct loss penalty caused by graphene on the four-wave mixing response of the whole waveguide. Accordingly, graphene-induced losses give rise to a drop in the conversion efficiency of approximately 10 between the 2 -long SiN structure without graphene and the SiN waveguide locally covered by 1.1 of graphene. This CE decrease primarily reflects the missing nonlinear contribution from the 1.7 -long SiN waveguide following the hybrid graphene/SiN section, whenever the latter strongly absorbs the signal, i.e., for increasing graphene lengths. However, the measurements (in black) suggest an even stronger reduction of the conversion efficiency induced by the presence of graphene, which reaches up to −5 with respect to the yellow curve, for the maximum pump power (7 W) and graphene length (1.1 mm).
As observed on
Figure 5, the four-wave mixing
CE measurements are relatively well reproduced considering −10 m
−1W
−1 <
< −7 m
−1W
−1 i.e., an absolute value almost one order of magnitude larger than that of the clad Si
N
waveguide. This single value for
consistently reproduces the measurements for the whole range of graphene lengths and pump powers shown in
Figure 5, apart from the sole
measurement, which is a bit less reliable due to lower S/N and the lack of sensitivity of our setup. This estimated
value is also found to be negative, because positive values would further deviate from the measured trend as compared with the case
=
(yellow curve). Considering that the nonlinear response of Si
N
yields a positive nonlinear parameter, a negative effective
implies that the nonlinearity induced by graphene itself is negative, as was suggested by earlier reports [
12,
23], and is strong enough to overcompensate the nonlinearity of the underlying Si
N
waveguide on which it is deposited. Qualitatively, this
CE-enhanced reduction can be understood by the opposite contributions to the idler signal generation arising, respectively, from the hybrid graphene/Si
N
section (with a negative
) and the clad Si
N
sections (with a positive
) before and after the section covered by graphene. Although our measurements and simulations allow us to quantify and demonstrate some nonlinear enhancement of the waveguide locally provided by graphene, this suggests that our particular Si
N
chip is here not ideal to practically exploit the nonlinear effects of graphene, as the contributions of consecutive sections with opposite signs are undoing each other.
A comparison of our results with the literature on graphene-covered nonlinear waveguides is shown in
Table 2. The references used here correspond to similar SPM and FWM experiments at telecom wavelengths, yet under slightly different experimental conditions (either CW or pulsed signals and various power levels), as indicated in the table. Regarding first the sign of the nonlinear response of graphene that was extracted from the different experiments, our negative value is consistent with the reports from Vermeulen’s group [
12], but it is the opposite of what was found by others. We argue, in particular, that the extraction of graphene nonlinearity is quite tricky from measurements on hybrid graphene/Si waveguides, in which the positive nonlinear contribution from the underlying Si waveguides (almost two orders of magnitude larger than for Si
N
waveguides) cannot be ignored. We also highlight that it is not straightforward to extract the sign of the nonlinearity from FWM measurements. In our case, the opposite nonlinear contribution from the Si
N
waveguide sections (for which
) combined with our simulations enabled us to reliably access the negative value for the effective nonlinear parameter of our hybrid graphene/Si
N
section.
Focusing on the four-wave mixing measurements presented in the
Table 2, our work is conducted in a pulsed regime rather than by using CW signals, allowing us to test the waveguides under several watts instead of tens [
10] or hundreds of milliwatt [
16]) powers, which should boost, in principle, the nonlinear effects observed. According to the usual perturbative description of the nonlinear response of dielectric materials, the power variation should yet not affect the hybrid waveguide nonlinear response. We observe, however, a striking difference between the effective nonlinear
for the graphene-covered waveguides, which is, in our case, of approximately 10 m
−1W
−1, i.e., several orders of magnitude less than the other values indicated in the Table, where up to a few thousand were reported. Considering the 2D nature of graphene covering the waveguide and assuming that the nonlinear response of graphene is a few orders of magnitude greater than the one of the underlying waveguide (which is particularly true for Si
N
), one could rightly argue that the net nonlinear effect of the hybrid section is constrained by the interaction between the guided mode and the graphene covering the waveguide. By using the effective thickness approach to describe graphene as a thin (typically ∼0.3-nm thick) but standard material, an equivalent nonlinear effective index
can be inferred from the mode field distribution and its overlap with graphene [
11]. Although this approach is questionable for 2D materials, and a more relevant approach could be considered [
22], it allows us to more simply compare the values inferred for the graphene nonlinear response between different hybrid waveguides, after factoring out the variations in the light–graphene interaction between the different underlying waveguide geometries used in
Table 2. By using this approximated approach, our results suggest an equivalent nonlinear index for graphene of
m
2W
−1. This value remains one order of magnitude lower than the one extracted from the references in
Table 2, which lead to
m
2W
−1 for most of them.
From a practical point of view, to evaluate the efficiency of graphene-covered waveguides within photonic circuits, we should take into account both the additional loss penalty and the nonlinear contribution provided by the local addition of graphene. The graphene-induced linear propagation loss of the hybrid mode (denoted as
) reasonably reflects the degree of interaction between the guided mode and graphene. This loss varies quite significantly between the different waveguide geometries of
Table 2, and ranges between 86 dB cm
−1 for our work up to several hundred and over 1000 dB cm
−1 for others. To some extent, the nonlinear contribution of graphene to the overall nonlinear response of the hybrid guided mode should similarly increase with this interaction. It is thus relevant to compare the ratio
/
for the different reports so as to quantify the amount of nonlinear effects that a particular hybrid graphene/dielectric structure can produce per unit of power, normalised with respect to the associated graphene-induced loss penalty. The last column of
Table 2 shows this quantity for the different graphene/dielectric waveguide geometries. This ratio remains much weaker, in our case—between one and two orders of magnitude lower with respect to the other references. Therefore, despite the relative lower loss of our structures, reflecting a weaker interaction of graphene with the guided mode, much lower nonlinear effects can be comparatively achieved. In our geometry, the gain in the nonlinear parameter thus remains marginal relative to the loss penalty locally induced by graphene.
Having thus factored out the effect of the loss and the light–graphene interaction, we are left to explain the striking difference in the graphene nonlinear response from that of the literature. Although the perturbative approach for describing nonlinearities traditionally leads to a power-independent
response for dielectric materials, this description might fail for nontransparent graphene, in which nonlinear effects tend to be mediated by significant power-dependent carrier dynamics. This might explain that nonlinear effects in graphene might not increase under larger powers as for traditional dielectric media. Instead, the
parameter of graphene-covered waveguides might provide an effective nonlinear response that is only valid for a limited range of powers, and the absolute value of this effective response would appear to decrease with increasing powers according to our experiments with respect to prior works conducted at lower powers. This hypothesis is supported by the recent theoretical description of the nonlinear response of graphene [
21,
22,
31], in which a strong link between the carrier dynamics in graphene and its nonlinear response has been established.
5. Design Rules for Optimizing Hybrid Nonlinear Waveguides Locally Coated with 2D Materials
One last question arising from this investigation is a practical one: are there some optimum conditions to leverage the graphene nonlinearity when exploiting graphene patches, or more generally 2D material patches, integrated within a passive photonic circuit? Let us suppose that a patch of 2D material (either graphene or any other 2D material) of length
is locally positioned on top of a longer passive and low loss waveguide, which exhibits a positive (and relative low) nonlinear parameter value
, associated to the bare waveguide. Taking four-wave mixing as an example, two scenarios are considered, one for each sign associated with the nonlinear coefficient
along the 2D material-covered section of the waveguide. In order to derive simpler analytical expressions (see
supplementary information), we only consider a 2D material-covered section of waveguide followed by a bare section of waveguide without 2D material. Compared to the case used in our experiments, it is equivalent to neglect the contribution from the first short bare waveguide section, which just adds an offset to the nonlinear effect measured from the whole waveguide.
Our analysis (see
supplementary information) shows that, under these conditions, the trend observed on the idler generation as a function of the 2D material length critically depends on whether the nonlinear parameter ratio,
/
, between the hybrid 2D material/dielectric waveguide and the bare waveguide is greater or lower than the product
, with
the effective length of the bare section of waveguide after the 2D material. Note that the quantity
eventually equals the
/
loss ratio of the hybrid waveguide with respect to bare waveguide, for long passive circuits.
More quantitatively, we show that the 2D material patch increases the net idler generation of the overall waveguide as long as one of the following conditions are met, depending on the sign of
:
In both cases, there is an optimum 2D material length that optimizes the nonlinear idler generation, which is found to be
Eventually, if the 2D material nonlinear response is strong enough to satisfy Equation (
4), then the 2D material nonlinearity can benefit the net nonlinear effect accumulated across the entire waveguide. Otherwise, the net nonlinear effect decreases for increasing 2D material length, because the loss penalty is larger than the nonlinear enhancement locally imparted by the 2D material. Both scenarios are illustrated in
Figure 6, for either sign of
.
We note that the constraint set by Equation (
4), on
is harder to meet when the latter has an opposite sign to that of the bare waveguide (
Figure 6b), with respect to the situation where both the 2D material nonlinearity and that of the underlying waveguide material jointly contribute to the nonlinearity (
Figure 6a). Consistently, the optimum 2D material length given by Equation (
5) also tends to be larger when the sign of the nonlinear parameters are different. In the case in which
, we also observe the existence of a minimum idler power, which is almost canceled out, for a length of 2D material that compensates for the opposite nonlinear contribution from the bare waveguide.
In our experiments, . Therefore, a net increase of idler intensity with graphene could appear with values greater than m−1W−1 (if ) or than m−1W−1 (if ). Our experimental results, and the comparison to more refined simulations (which yield m−1W−1), indicate that the benefit of graphene, in our chip, is not only overcompensated by its loss penalty, but is also not enough to compensate for the nonlinear effect induced by the surrounding sections of SiN.
One simple solution to experimentally observe a benefit of graphene on photonic circuits would be to shorten the surrounding bare waveguides or, more practically, to decrease its nonlinearity. For instance, other waveguide platforms might be used, such as SiO
glass waveguides with 10 times lower nonlinearity that might boost the ratio
via a reduction of
for the surrounding waveguide. Without changing the material platform, one could alternatively decrease the
parameter of the surrounding waveguide by enlarging its cross-section while keeping it locally small to increase the graphene–light interaction only where graphene is to be integrated. Finally, other 2D materials with a better loss/nonlinearity trade-off (i.e., a higher
/
), like graphene oxide [
32,
33], might be better suited for nonlinear applications.