Is it possible to sum over indices below or above fermi surface in sympy.secondquant package

Question

I am trying to do a summation over fermion index in the variable "expr" using sympy.secondquant. Some of my code is shown below

from sympy import *