Interesting question:
One way can be to set the attractive term = 0*(rmin/r)^6, which means that one can use the USERE-hookup in CHARMM for a user-energy defined potential. Have just begun to learn more about this facility, but see usage.doc