Magnet-Free Nonreciprocal Edge Plasmons in Optically Pumped Bilayer Graphene

TL;DR

This paper shows that bilayer graphene can support nonreciprocal edge plasmons without magnets, using light to control their behavior more efficiently than monolayer systems.

Contribution

The study introduces stronger nonreciprocal plasmons in bilayer graphene with lower pump requirements and tunable energy gaps.

Findings

Bilayer graphene shows stronger nonreciprocity in edge plasmons than monolayer systems.

Lower pump amplitude is needed in bilayer graphene to achieve plasmon splitting.

The gate-tunable energy gap in bilayer graphene offers enhanced control over plasmonic behavior.

Abstract

Recent theoretical studies have shown that gapped Dirac materials (such as gapped monolayer graphene) optically pumped with circularly polarized light can host edge-localized plasmon modes with nonreciprocal dispersions driven by valley population imbalance. Here, we extend this framework to Bernal-stacked bilayer graphene. Using the Wiener–Hopf method, we compute the exact edge plasmon dispersion, confinement length, and electric potential. Our results show that bilayer graphene exhibits stronger nonreciprocity in edge plasmons, requiring approximately one order of magnitude lower pump amplitude to achieve splitting compared with monolayer Dirac systems. Furthermore, the gate-tunable energy gap of bilayer graphene provides an additional degree of control, positioning optically pumped bilayer graphene as a versatile platform for magnet-free nonreciprocal plasmonics.

Figures

Full text

1. Introduction

Plasmons can confine light far below the diffraction limit, enabling the development of compact tunable elements for integrated photonics [1,2,3,4,5]. Breaking reciprocity so that plasmon modes preferentially propagate in one direction is particularly important for applications such as on-chip isolators and circulators [6,7,8,9,10,11]. Achieving nonreciprocal plasmon dispersion requires breaking time-reversal symmetry to lift the degeneracy between counterpropagating plasmons. The most common method of producing nonreciprocal plasmons uses static magnetic fields [12,13,14,15,16,17,18]. However, this method requires strong external magnetic fields and is incompatible with compact integrated photonic platforms [19,20,21].

Recently, alternative magnetic-free approaches to plasmonic nonreciprocity have been proposed [17,18]. In gapped Dirac materials, a nonzero Berry flux gives rise to collective edge modes with nonreciprocal plasmon dispersions in opposite propagation directions, even in the absence of an external magnetic field. This nonreciprocity occurs because the Berry curvature induces an anomalous velocity that assists or opposes charge motion along the edge. One way to induce a nonzero net Berry flux is to drive gapped Dirac materials with circularly polarized light, which selectively pumps carriers for different valleys. This valley-selective pumping leads to a population imbalance between valleys with opposite Berry flux, yielding a nonzero net flux [17]. These studies have established that optically pumped gapped monolayer graphene supports chiral edge plasmons; however, the potential benefits of multilayer systems beyond monolayer Dirac-like systems remain largely unexplored.

Our central question is whether bilayer graphene can deliver quantitatively stronger and practically tunable nonreciprocal edge plasmons. Unlike monolayer graphene, bilayer graphene has a gate-tunable energy gap that arises from a perpendicular displacement field breaking layer inversion symmetry [22,23,24,25,26,27]. This tunability provides additional control to optimize the system’s optical response. Moreover, when gapped, bilayer graphene exhibits a Mexican-hat dispersion, in which the density of states near the energy gap diverges as the inverse square root of energy—unlike the moderately linear increase in monolayer graphene. Under circular pumping, this feature implies that bilayer graphene can support a larger valley population imbalance.

In this study, we demonstrate that these unique features of bilayer graphene lead to significantly larger and more resolvable edge-mode splitting. We modeled the nonequilibrium carrier distribution under continuous circular pumping using a density-matrix approach and computed the corresponding nonequilibrium optical conductivity. Applying the Wiener–Hopf technique, we provide an exact analysis of edge plasmons, presenting the nonreciprocal dispersion of edge plasmons, the electric potential profile, and the confinement length. While our monolayer results reproduce the previously reported chiral-plasmon phenomenology, bilayer graphene exhibits much stronger nonreciprocal splitting: for comparable gap sizes, the required pump intensities are nearly an order of magnitude lower to resolve the frequency difference between counterpropagating modes. Furthermore, the gate-controllable energy gap of bilayer graphene enhances or suppresses nonreciprocity depending on pump strength. In addition, in previous theoretical work, the demonstration of nonreciprocal edge plasmons is shown under a large energy gap (∼0.5 $[eqn]$ ) and a strong pumping field with the electric field amplitude reaching $[eqn]$ [17]. Our work, by contrast, not only extends the scope to bilayer graphene but also considers smaller energy gaps and much weaker pump intensities (e.g., $[eqn]$ ) that are more readily achievable for on-chip photonic devices. We show that, in this parameter range relevant to on-chip devices, the nonreciprocal splitting is still sizable. Our findings identify optically pumped bilayer graphene as a versatile platform for tunable, magnet-free, nonreciprocal plasmonics and motivate experimental exploration of chiral edge plasmons in gate-tunable bilayer graphene devices.

2. Theoretical Framework

2.1. Edge Plasmons

We consider a semi-infinite two-dimensional graphene sheet strictly confined to the $[eqn]$ plane and occupying the $[eqn]$ region, with edge plasmons propagating along the y-axis (see Figure 1). The two-dimensional conductivity tensor of the system, which exhibits an abrupt jump at $[eqn]$ at the boundary of the graphene sheet, is expressed as follows:

$$[eqn]$$

where $[eqn]$ is the Heaviside step function and $[eqn]$ is the Dirac delta function. Assuming a plane wave dependence, the edge plasmon has an electric potential of the following form:

$$[eqn]$$

The associated oscillating charge density, strictly confined to the $[eqn]$ plane, is expressed as follows: $[eqn]$ , where $[eqn]$ for $[eqn]$ . The electrical potential and charge density are related through the continuity equation $[eqn]$ , where $[eqn]$ is the current density. Using Ohm’s law $[eqn]$ and $[eqn]$ , where $[eqn]$ is the electric field, the continuity equation becomes $[eqn]$ . By substituting Equations (1) and (2) into the continuity equation, we obtain $[eqn]$ , where $[eqn]$ , Given that $[eqn]$ and $[eqn]$ , owing to the system’s isotropy and broken time-reversal symmetry under circularly polarized light, the equation for $[eqn]$ is simplified to the following:

$$[eqn]$$

The term $[eqn]$ arises from the discontinuity of the system at the boundary $[eqn]$ , where the conductivity tensor changes abruptly. In the non-retarded regime $[eqn]$ , the electric potential is determined by the Poisson equation $[eqn]$ , where $[eqn]$ is the effective background permittivity. At $[eqn]$ , this reduces to the following integral form:

$$[eqn]$$

where

$$[eqn]$$

is the kernel of the integral equation and $[eqn]$ is the modified Bessel function of the second kind [12,14].

The exact electric potential and dispersion of edge plasmons can be obtained by solving Equations (3) and (4) using the Wiener–Hopf method, which is well-suited for problems with boundary conditions in a semi-infinite domain, as is the case here. The basic idea is that the real-space integral and differential equations are Fourier transformed, and the kernel of the transformed equation is factorized into two new functions, one of which is analytic in the upper half of the complex plane and the other in the lower half. This factorization separates the equation into two parts, each of which is valid on a different half-plane. This is solved using the arguments of analytic continuation. For technical details of the derivation, interested readers are referred to the existing literature [28,29,30,31,32,33,34]. In this paper, we present only the final expressions for the electrostatic potential and dispersion relation obtained using the Wiener–Hopf technique. The edge plasmon dispersion $[eqn]$ is determined from the following nonlinear equation [34]:

$$[eqn]$$

where $[eqn]$ and $[eqn]$ are dimensionless longitudinal and Hall conductivities, respectively. Here $[eqn]$ and $[eqn]$ , where $[eqn]$ is the polylogarithm function of order 2. The electric potential function of an edge plasmon in the transverse plane $[eqn]$ , where $[eqn]$ corresponds to the region $[eqn]$ and $[eqn]$ to $[eqn]$ , is given by the following expressions [35]:

$$[eqn]$$ $$[eqn]$$

where

$$[eqn]$$ $$[eqn]$$

Here, $[eqn]$ , $[eqn]$ and $[eqn]$ are defined as the complex roots of the two-dimensional dielectric function $[eqn]$ with $[eqn]$ ( $[eqn]$ ) having a positive (negative) imaginary part, and $[eqn]$ . Here, we denote by $[eqn]$ the upper/lower half-plane factors of the kernel with $[eqn]$ . $[eqn]$ are explicitly written as $[eqn]$ and $[eqn]$ with the following:

$$[eqn]$$

Throughout we use the principal branches for the square roots.

Alternatively, we can use the Fetter approximation to obtain the edge plasmon dispersion, where the integral kernel $[eqn]$ [Equation (5)] is replaced by $[eqn]$ with the same area and second momentum [12]. This yields the simple equation for the dispersion relation of the edge plasmon, given by $[eqn]$ where the +(−) sign in the last term corresponds to the $[eqn]$ ( $[eqn]$ ) plasmon edge mode propagating in the left (right) direction [see Figure 1]. In the low-frequency and clean limit, the longitudinal optical conductivity exhibits typical Drude behavior $[eqn]$ , where D is the Drude weight, and the Hall conductivity $[eqn]$ is approximated as a constant independent of $[eqn]$ . Using these, it is straightforward to obtain the edge plasmon dispersion in the long-wavelength limit ( $[eqn]$ ), which is expressed as follows:

$$[eqn]$$

where $[eqn]$ is the bulk plasmon dispersion.

Notably, Equations (6) and (12) highlight the central role of optical conductivity in determining the nonreciprocity of the edge plasmon. In particular, the Hall component of the conductivity $[eqn]$ breaks the symmetry between modes traveling in opposite directions, which is directly responsible for the nonreciprocal propagation. Accordingly, in the next section, we present a theoretical model for the conductivity of graphene under circularly polarized light.

2.2. Graphene Model

Our analysis is based on a low-energy effective Hamiltonian for gapped monolayer and bilayer graphene, which captures the essential features of low-energy physics near the $[eqn]$ and $[eqn]$ valleys. The Hamiltonians for gapped monolayer and Bernal-stacked bilayer graphene with an energy gap U are given by $[eqn]$ and $[eqn]$ , respectively, where $[eqn]$ is the momentum vector for valley index $[eqn]$ , $[eqn]$ is the Fermi velocity of graphene, U is the interlayer potential bias, and t is the interlayer hopping energy [36,37]. The matrices $[eqn]$ and $[eqn]$ are Pauli matrices for sublattice and layer pseudospin, respectively, while $[eqn]$ and $[eqn]$ denote $[eqn]$ identity matrices.

2.3. Steady-State Charge Distribution and Optical Conductivity

Under continuous optical pumping, the carrier distribution reaches a time-independent steady-state nonequilibrium. The steady-state distribution function is obtained by solving the von Neumann equation, which describes the evolution of the density matrix under the influence of an electromagnetic field: $[eqn]$ where $[eqn]$ is the bare Hamiltonian without interaction with the electromagnetic field, $[eqn]$ is the density matrix, and $[eqn]$ is the interaction term that considers coupling with the electromagnetic field with $[eqn]$ being the electric field of the circularly polarized monochromatic pumping light. Using the rotating-wave and relaxation-time approximations, the diagonal components of the density matrix are obtained by solving the following system of linear equations [17]:

$$[eqn]$$

where $[eqn]$ , $[eqn]$ , $[eqn]$ , and $[eqn]$ with $[eqn]$ and $[eqn]$ denoting the electron-hole and valley relaxation times, respectively. $[eqn]$ is the dephasing time required for decoherence between the conduction and valence bands and $[eqn]$ is the valley-dependent dipole matrix element that determines the optical selection rule, where $[eqn]$ is the electronic wavefunction and $[eqn]$ is the lowest conduction (highest valence) energy band. Unless otherwise specified, all numerical results were obtained with realistic parameter values: $[eqn]$ and $[eqn]$ . Without loss of generality, we assume right circularly polarized light (i.e., $[eqn]$ ) where $[eqn]$ denotes the pumping field amplitude.

Figure 2 shows the numerically calculated steady-charge distribution functions at $[eqn]$ and $[eqn]$ for both monolayer and bilayer graphene with a finite energy gap of $[eqn]$ . The pump frequency was chosen to match the energy gap to maximize the optical pumping effects, and the temperature was set to $[eqn]$ to clearly reveal valley-polarized optical excitation. Unless otherwise stated, the remainder of this study uses a room temperature of $[eqn]$ . The selection rule is captured in Figure 2g,h which plot the magnitude of the field-projected dipole matrix element, $[eqn]$ . For monolayer graphene, the selection rule forbids interband optical pumping at the $[eqn]$ Dirac point, whereas the selection rule in bilayer graphene forbids optical pumping at both Dirac points (i.e., $[eqn]$ and $[eqn]$ valleys), i.e., the rule is completely valley independent [38,39]. The valley-dependent selection rule reemerges at finite momentum $[eqn]$ around each valley, where $[eqn]$ . Importantly, Bernal-stacked bilayer graphene exhibits a Mexican-hat energy dispersion, resulting in a diverging density of states near the energy gap, in contrast to gapped monolayer graphene, where the density of states increases linearly without peaks near the gap [see Figure 2c,f]. This enhanced density of states leads to more efficient optical absorption in the bilayer graphene. Consequently, more charges are pumped in the bilayer graphene for the same pumping field strength, as shown in Figure 3a, where we plot the valley imbalance $[eqn]$ , with $[eqn]$ and $[eqn]$ denoting the photoexcited carrier densities at $[eqn]$ and $[eqn]$ valleys, respectively.

We take an approximate approach to calculate the optical conductivity of graphene driven by continuous circularly polarized light using the standard Kubo formula with the equilibrium Fermi–Dirac distribution replaced by the nonequilibrium distribution function evaluated using Equation (13) [40,41,42]:

$$[eqn]$$

where g is the spin degeneracy, $[eqn]$ is the band energy, $[eqn]$ , $[eqn]$ are band indices, $[eqn]$ , and $[eqn]$ . where $[eqn]$ denotes the velocity operator. The parameter $[eqn]$ is a phenomenological broadening parameter that incorporates multi-body and disorder effects. This approximation neglects coherent dressing of the electronic bands, which is justified in the weak-driving regime quantified by the dimensionless quantity $[eqn]$ [43,44,45,46,47,48,49,50,51], where $[eqn]$ is the pump frequency. Typical experimental parameters relevant to our work (e.g., $[eqn]$ , $[eqn]$ , $[eqn]$ ) ensure that we remain within the weak-driving limit. In all numerical calculations, we use a realistic broadening parameter $[eqn]$ [52,53,54], which is appropriate for moderately clean graphene. The momentum integral is carried out over $[eqn]$ and $[eqn]$ using an adaptive quadrature method with the cutoff momentum $[eqn]$ chosen so that the conductivity barely changes upon doubling $[eqn]$ . Because the integrand contains weak and integrable singularities, we enforce a minimum of $[eqn]$ subintervals in the adaptive quadrature to ensure robust convergence and accuracy.

As Hall conductivity breaks reciprocity between counterpropagating plasmon edge modes, it is instructive to first examine the photoinduced Hall conductivity before presenting our numerical results for edge plasmons. Although our full numerical calculations of the edge plasmons used the dynamic photoinduced Hall conductivity, we present the static photoinduced Hall conductivity contribution at $[eqn]$ and $[eqn]$ valleys in Figure 3. The $[eqn]$ dependence of the Hall conductivity is weak unless the plasmon frequency is comparable to the energy gap, so the static Hall conductivity captures the relevant trend. Note that the Hall conductivity has opposite signs owing to the opposite chirality of the $[eqn]$ and $[eqn]$ valleys. For a weak electric field strength, only a small fraction of the charges is pumped, and the Hall conductivity for both monolayer and bilayer is negligibly small. However, above $[eqn]$ , $[eqn]$ in bilayer graphene grows very rapidly, reaching $[eqn]$ , whereas in monolayer graphene, $[eqn]$ reaches the same magnitude only at much higher field strengths $[eqn]$ . Thus, bilayer graphene exhibits significantly stronger nonreciprocity than monolayer graphene.

After obtaining the optical conductivity, the next section presents the numerical results for plasmon dispersion and localization length, obtained using Equations (6)–(8) as presented in the previous section.

3. Results and Discussion

Figure 4 compares the edge plasmon dispersion (solid lines) with the bulk graphene plasmon (black dashed line) for monolayer [panels (a–c)] and bilayer [panels (d–f)] graphene as the pump amplitude increases from $[eqn]$ , 5, to $[eqn]$ . For the weakest pump amplitude, $[eqn]$ [panels (a) and (d)], both edge plasmon modes are displaced relative to the bulk plasmon mode by nearly the same amount, indicating weak nonreciprocity. As $[eqn]$ increases, the nonreciprocity becomes stronger: the left (right-)-propagating mode approaches (moves away from) the bulk mode. This trend follows the increase in the photoinduced Hall conductivity, as shown in Figure 3. At $[eqn]$ , the left-propagating branch edge mode merges with the bulk mode at a critical momentum $[eqn]$ ; for $[eqn]$ , only the right-propagating mode remains, rendering the edge plasmon effectively unidirectional. We emphasize that the edge plasmon nonreciprocity is much more pronounced in bilayer graphene, consistent with its higher photoinduced Hall conductivity discussed in the previous section. In particular, achieving a comparable level of nonreciprocity in bilayer graphene requires approximately an order of magnitude smaller $[eqn]$ than in monolayer graphene.

The stronger nonreciprocity in bilayer graphene is more pronounced in Figure 5, where we present the dimensionless ratio $[eqn]$ for the three pump amplitudes used in Figure 4. Here $[eqn]$ is the plasmon frequency splitting, $[eqn]$ denotes the right/left–propagating branches, and $[eqn]$ is the plasmon decay rate. The dimensionless ratio $[eqn]$ quantifies the spectral resolvability of the two counterpropagating edge modes with $[eqn]$ , implying that their plasmon frequencies are easily resolved. For the weak pump, $[eqn]$ , $[eqn]$ in both mono- and bilayer graphene. Thus, the edge plasmon nonreciprocity is effectively unobservable. Increasing $[eqn]$ enhances $[eqn]$ . In bilayer graphene, $[eqn]$ is well above unity over a broad momentum range of $[eqn]$ . By contrast, in monolayer graphene, a larger pump amplitude $[eqn]$ is required to reach $[eqn]$ . At the same pump amplitude ( $[eqn]$ ), $[eqn]$ reaches $[eqn]$ in bilayer graphene, significantly larger than that for monolayer graphene. The bilayer curves terminate at the critical momentum $[eqn]$ , where the left–propagating edge mode merges with the bulk plasmon and ceases to exist, consistent with the onset of unidirectionality.

In Figure 6, we plot the spatial profile of the normalized electrostatic potential $[eqn]$ of the edge plasmons propagating in monolayer (a,c) and bilayer (b,d) graphene in the transverse plane ( $[eqn]$ ) at fixed $[eqn]$ . The black curve corresponds to the $[eqn]$ contour and serves as a measure of transverse confinement. A comparison of panels (a) and (c) with panels (b) and (d) shows that the potential profile depends on the propagation direction, indicating that nonreciprocity manifests not only in the dispersion but also in the transverse confinement length. For $[eqn]$ , both monolayer and bilayer edge modes are tightly localized with nearly identical confinement lengths, and the bilayer is only slightly more localized. This observation is consistent with the monolayer and bilayer edge-plasmon dispersions being well-separated from the bulk plasmon dispersion for $[eqn]$ . In contrast, for $[eqn]$ , which is close to the critical momentum $[eqn]$ where the bilayer edge and bulk plasmon modes merge, the bilayer mode substantially broadens into the interior, exhibiting a much longer confinement length. The monolayer mode is well-localized, similar to the $[eqn]$ case, because the plasmon dispersion is well-displaced from the bulk mode. We emphasize that the overlap between the counterpropagating edge fields is much smaller in the bilayer, implying that elastic backscattering of edge plasmons by imperfections is more strongly suppressed in the bilayer graphene [16,55,56,57].

In dual-gated bilayer graphene, a perpendicular displacement field breaks the inversion symmetry and induces an interlayer potential asymmetry U, allowing continuous gate tuning of the energy band gap [22,23,24,25,26,27]. Accordingly, we examine how the tunable interlayer bias U affects edge-plasmon nonreciprocity. In Figure 7a, we plot the photoinduced Hall conductivity as a function of pump amplitude $[eqn]$ for $[eqn]$ and $[eqn]$ . As established above, edge-mode nonreciprocity increases with photoinduced Hall conductivity $[eqn]$ ; thus, Figure 7a provides a practical guide to how nonreciprocity evolves with U. For a weak field, $[eqn]$ , $[eqn]$ is larger for the smaller gap. However, with increasing $[eqn]$ , $[eqn]$ increases more rapidly for the larger U, and the ordering reverses at $[eqn]$ . This behavior carries over to the edge-mode splitting shown in Figure 7b, where $[eqn]$ is plotted as a function of $[eqn]$ for fixed values of $[eqn]$ , 10, and $[eqn]$ . Note that, for $[eqn]$ and $[eqn]$ , $[eqn]$ follows the same order as $[eqn]$ . For a strong pump amplitude of $[eqn]$ , the ordering reverses with $[eqn]$ , which is almost twice as large as that for $[eqn]$ . In summary, our analysis shows that large nonreciprocity can be achieved with a large energy gap if a large $[eqn]$ is available. However, for weak $[eqn]$ , a small energy gap is more favorable for producing nonreciprocal edge plasmons. Our results show that, for sufficiently large pump amplitudes, increasing the gate-induced gap U enhances the nonreciprocal splitting $[eqn]$ , whereas in the low-field regime, a smaller gap yields a larger $[eqn]$ .

Our analysis so far assumes that the pump photon energy is resonant with the band gap U, which we vary between $[eqn]$ and $[eqn]$ . This mid-infrared band is well covered by established sources such as quantum cascade lasers (QCLs) [58,59,60,61,62], which are already mature for room-temperature continuous-wave operation across this range. Moreover, on-chip integration of QCLs has been shown via transfer-printing and heterogeneous bonding onto Si/Ge-on-Si platforms, and QCLs have been used to excite graphene on a chip-scale mid-IR device [58,59]. It should be noted, however, that operation beyond $[eqn]$ can be problematic due to heating effects such as rapid optical-phonon emission [63,64], resulting in shorter valley relaxation times. While operation up to $[eqn]$ still remains feasible with pulsed or low-duty excitation and proper thermal management, a practical upper is approximately < $[eqn]$ so that coupling to optical phonons is suppressed.

We end this section with a discussion on the role of stacking order on the edge plasmon nonreciprocity. For AA stacked bilayer graphene, the electronic band structure is gapless, with two Dirac cones shifted in energy by the interlayer hopping [65]. A perpendicular electric field does not open an energy gap, only pushing the two Dirac cones further apart in energy. In such metallic cases, where the Berry curvature (and thus the Hall conductivity) is negligible [66], the nonreciprocity of edge plasmons is expected to be extremely small, in contrast to the AB-stacked bilayer graphene studied in our work. For BA stacking, reversing the stacking from AB to BA at a fixed energy gap should flip the sign of the Hall conductivity and therefore reverse the direction of nonreciprocity. Furthermore, this result for BA stacking indicates that an AB-BA stacking domain wall under a perpendicular electric field flips the sign of the Hall conductivity across the interface. In our electrodynamic formulation, this acts as a 1D boundary with a discontinuous conductivity tensor [c.f., Equation (1)], yielding a guided domain-wall plasmon. The details of this plasmon mode can be obtained at the same approximation level as our edge-mode analysis and are left for future work.

4. Conclusions

We investigated edge plasmons in gapped monolayer and bilayer graphene under continuous circularly polarized light. Optical pumping induces a population imbalance between valleys, and the resulting photoinduced Hall conductivity breaks the reciprocity of counterpropagating plasmon edge modes without requiring a magnetic field. Using the exact Wiener–Hopf method, we calculated the edge plasmon dispersion, confinement length, and spatial profile of the plasmonic edge-mode electric field, demonstrating that circular pumping in gapped bilayer graphene also induces strong nonreciprocity in edge plasmons predicted for gapped Dirac materials. We found that bilayer graphene, owing to its Mexican-hat electronic band structure, supports significantly stronger nonreciprocal splitting, achieving resolvable chiral splitting ( $[eqn]$ ) at pump fields nearly an order of magnitude weaker than those required for gapped monolayer graphene. Our results further indicate that increasing the interlayer bias U of bilayer graphene enhances the plasmon edge-mode nonreciprocity in strong fields but reduces it in weak fields. Because U is gate-tunable, this additional control provides a practical means of optimizing nonreciprocity in bilayer-graphene plasmonic devices for a given pump. Overall, these findings suggest that optically pumped bilayer graphene is a promising candidate for compact, tunable, and nonreciprocal plasmonic devices operating in the mid-infrared regime.