Dear Kwant core team, Dirac or Majorana Fermions might be used in kwant with Discretizing continuous Hamiltonian on a square lattice. However the way that kwant descritizes the Dirac Hamiltonian gives rise to what is called doubling Fermions (even the desicritization preserves the hermicity ’symmetric derivative') . In this context, several approaches have been used to get rid of the extra Dirac cones in the Brillouin zone. Without going in deeper details, each of those proposed methods has some issue related to preserving either chiral or symplectic symmetry.To avoid this issue recent methods goes beyond nearest neigbour interaction and then my question come to raise. *How can we include second or third nearest interaction in kwant.continuum.discretize *. Best regards Adel