I was running computations with background charges using
Chrgfield = QMMM()
Chrgfield.extern.addCharge( 5.05881, 6.24500, 0.28003, -0.70717)
but I wonder if there is any way to add point dipoles to a QM calculation, probably in a QM/MM scheme, similar to the background charges shown above? I mean, add a set of dipoles at certain xyz positions, similar to what some force fields do.
If yes, how to do it in PsiAPI and/or Psithon mode?