▎ 摘　　要

Here, we benefit from the strong nonlinear response of graphene and the rich variety of resonances provided by a graphene/dielectric Fibonacci multilayer to realize bistability and multistability in the terahertz (THz) frequency range. Toward this pursuit, we employ the nonlinear transfer matrix method. We examine the suitability of resonances in the Fibonacci multilayer for the bi/multistability purposes and determine the proper working point. We report various switching up/down manners via single or stepwise jumps between states of the same or different contrasts upon increasing followed by decreasing the intensity of the incident wave. We show that graphene samples of high quality are preferred for bi/multistable switching in terms of reducing the switch-up/-down thresholds and widening the multistable region. We also explore the possibility of tuning the bi/multistable behavior via the frequency and angle of the incident wave as well as the graphene Fermi level. We envision precious applications in THz switching, realizing logic gates, and so on for this system. (C) 2017 Optical Society of America