Using a VMD script written by me (MOisoall.tcl) I try to visualize the behavior of HOMO in the course of an inversion of ammonia (Fig. 1), but the VMD script does not work correctly (Fig. 2). Could someone check the MOisoall.tcl script and correct it ?. I attach the .cub files of the three-point HOMO of said reaction (MO_SP000 * .cub).
NOTE: to run the script it is necessary to move the .cub and MOisoall.tcl files to the VMD directory and write:
source MOisoall.tcl
in the VMD console (Fig. 3) and press the Enter button.