one can use the Bosonic operators of sympy.physics.secondquant, e.g $ b(k), b^{\\dagger}(k) $:
sympy.physics.secondquant
import sympy as sy import sympy.physics.secondquan