We study the mid-infrared plasmonic response in Bernal-stacked bilayer graphene. Unlike its monolayer counterpart, bilayer graphene accommodates optically active phonon modes and a resonant interband transition at infrared frequencies. They strongly modifies the plasmonic properties of bilayer graphene, leading to Fano-type resonances, giant plasmonic enhancement of infrared phonon absorption, narrow window of optical transparency, and a new plasmonic mode at higher energy than the classical plasmon.Plasmonics[1] is an important subfield of photonics that deals with the excitation, manipulation, and utilization of plasmons-polaritons [2]. It is a key element of nanophotonics [3], metamaterials with novel electromagnetic phenomena [4,5] and also has potential applications in biosensing [6].Recently, graphene has emerged as a promising platform for plasmonics [7]. It has many desirable properties such as gate-tunability, extreme light confinement, long plasmon lifetime, and plasmonic resonances in the terahertz to mid-infrared (IR) regime [8][9][10][11][12][13]. Spatially resolved propagating plasmons has been observed with scanning near-field optical microscope [14,15] In this paper, we discuss why Bernal AB-stacked bilayer graphene is important and interesting in its own right as a plasmonic material. Apart from a few theoretical studies of plasmons in bilayer graphene [23][24][25][26][27][28], there is still no experimental studies of bilayer graphene plasmonics. First indication that the plasmonic response in bilayer graphene might be very different than that of monolayer is its two prominent IR structures in its optical conductivity. IR optical measurements of bilayer graphene reveal a phonon peak at ω ≈ 0.2 eV, with a strong dependence of peak intensity and Fano-type lineshape on the applied gate voltage [29,30]. The interlayer coupling in bilayer graphene also results in two nested bands, which presents a set of doping dependent IR features [31][32][33]. This interband transitions between the two nested bands produced a conductivity peak at ω ≈ 0.4 eV in optical IR measurements [34][35][36]. The impact of these IR structures on the bilayer plasmonic response has not been studied. We found several novel plasmonic effects in bilayer graphene: (i) giant plasmonic enhancement of infrared phonon absorption, (ii) an extremely narrow optical transparency window, and (iii) a new plasmonic mode at higher energy than the classical plasmon.Bilayer graphene arranged in the Bernal AB stacking order is considered, with basis atoms A 1 , B 1 and A 2 , B 2 in the top and bottom layers respectively. The intralayer coupling is γ 0 ≈ 3 eV and the interlayer coupling between A 2 and B 1 is γ 1 ≈ 0.39 eV, an average of values reported in optical IR and photoemission measurements [34][35][36][37][38]. We work within the 4 × 4 atomic p z orbitals basis, i.e., where a † i and b † i are creation operators for the i th layer on the A/B sublattices. Within this basis, the Hamiltonian near the K point can be written as:, where σ i and...