[Kwant] Spin orbit interaction for honeycomb

Dear all, I note that we can easily consider the next-nearest-neighbor hopping in honeycomb, however, what is the model (Hamiltonian) for this code? In graphene, the intrinsic SOI is included in the reference:PRL 95, 226801 (2005) and we also have that in silicene: PRL 110, 026603 (2013). Is the model for Kwant the same as the two reference papers? Moreover,can we consider the Rashba SOI using the Tinyarray package in honeycomb? the hoppings along x and y direction are more complicated in honeycomb. Regards, Qingtian Zhang