Caling nwchem with a mpi communicator won't work for the parts of nwchem written on top of global arrays. My suggestion is to develop a simple interface using a scripting language like python. This isn't the fastest, but it's probably not realistic to do direct dynamics efficiently except for really small systems or the plane-wave code, so the interface time probably won't matter.