Fill Kwant system with position dependent parameters from Qsymm model
It would be great if one could create a Kwant system by filling a Qsymm model with position dependent parameters. The motivation for this is to study defects in systems with Hamiltonians that break particular symmetries.
Things to consider:
The workhorse that converts onsites and hoppings to Kwant value functions is Model.lambdify . That's over here.
Unlike kwant.continuum, which knows how to treat position dependence, Model.lambdify does not know that
So a minimal change that would allow to create models with position dependence would be to:
- Subclass a model and implement a different lambdify method
- Replace model parameters so that they acquire position dependece