Confluent supersymmetric algorithm for bilayer graphene

Abstract

External magnetic field profiles leading to equidistant and partially equidistant bilayer graphene spectra within the tight-binding model are obtained. This is achieved by implementing the integral and differential versions of the second-order confluent algorithm to the harmonic oscillator for arbitrary real factorization energies. Additionally, new Barut–Girardello and Gilmore–Perelomov coherent states for bilayer graphene are derived, for both diagonal and non-diagonal ladder operators. Their time evolution is analyzed, finding temporal stability and cyclic evolution in some cases. This fact is contrasted with the non-cyclic evolution of bilayer graphene coherent states obtained when using two different factorization energies. Likewise, the geometric phase and uncertainty product of the quadratures for the previously obtained coherent states are studied.