TeNPy Forum

Hello,
Is there an easy way to access the matrix representation of the dN operator for fermion sites? Specifically if the filling is 0.5, and so dN = n - 0.5, how can this be represented as a matrix?

Thank you