One additional thing, in case it was not clear: I also wish to exclude the calculation of all charge-charge / coulomb interactions between 1,4 bonded neighbors, for the entire molecule. Does "nrexcl = 3" (after the molecular name) take care of this?