We give a power sum (boson) realization for the higher Macdonald operators and apply it to calculating the homological invariants of the Hopf link, which include Khovanov-Rozansky as a special case.