I'm defining the phosphorene unit cell perfectly. However, when the leads are attached to the scattering region, the kwant increases the unit cell, and thus "duplicate" energy levels appear. I'm using the model with 10 interlayer hoppings and 5 interlayer hoppings. Apparently, in the model of 5 interlayer hoppings it doesn't happen! Thank you for your help now.