You asked about the function of Bq atoms in a BSSE calculation. Basically, in order to do a Counterpoise correction you need to define two "pieces" of the molecule. NWChem chooses to do this by defining "normal" atoms and "ghost (bq)" atoms. Gaussian, for example, defines "layer 1" and "layer 2".