I am somewhat confused by the use of a Kronecker Product in Madsen, Georg K. H. et al. “BoltzTraP2, a program for interpolating band structures and calculating semi-classical transport coefficients.” Comput. Phys. Commun. 231 (2018): 140-145. (Link to paper on arXiv). The equation I'm having trouble with is equation 8 which reads
$$ \sigma(\epsilon,T)=\int \sum_b {\bf v}_{b,{\bf k}}\otimes {\bf v}_{b,{\bf k}}\tau_{b,{\bf k}} \delta(\epsilon-\epsilon_{b,{\bf k}}){{d{\bf k}}\over{8\pi^3}} $$
What is the Kronecker product supposed to do here? My understanding is that the result of such an operation should be a matrix but that makes no mathematical sense to me in this context, so what is going on?