Multi-band pairing of ultrarelativistic electrons and holes in graphene bilayer

Multi-band pairing of effectively ultrarelativistic electrons and holes in asymmetrically biased graphene bilayer in strong coupling regime is considered. In this regime, the pairing affects both conduction and valence bands of the both graphene layers, and the order parameter is a matrix, which indices correspond to the bands. For band-diagonal s-wave pairing, we derive the system of multi-band gap equations for the gaps in the valence and conduction bands and solve it in the approximation of constant gaps and in the approximation of separable pairing potential. For a characteristic width of the pairing region of order of magnitude of the chemical potential, the gap values are not much different from single-band BCS estimations. However, if the pairing region is wider, then the gaps can be much larger and depend exponentially on its energy width. We also predict gapped and soliton-like oscillations of a relative phase of the gaps and unpairing of quarter-vortices at Kosterlitz-Thouless transition.