Twisted traces of modular functions on hyperbolic $3$-space

We compute analogues of twisted traces of CM values of harmonic modular functions on hyperbolic $3$-space and show that they are essentially given by Fourier coefficients of the j-invariant. From this we deduce that the twisted traces of these harmonic modular functions are integers. Additionally, we compute the twisted traces of Eisenstein series on hyperbolic $3$-space in terms of Dirichlet L-functions and divisor sums.