Using lattice QCD we are calculating the hadronic mixing matrix elements relevant to a large class of BSM neutral B meson mixing scenarios. I'll discuss recent motivations for improving the precision of these calculations, the relevant class of mixing operators (effective four-quark, \Delta flavor =2 interactions), and outline preliminary results.