Sybyl commands to animate vibrations

The vibrations calculated by MOPAC can be visualized by Sybyl. This option is to be carried out in command mode. We use the SN2 TS, discussed in the paragraphs 1A and 2A as an example, and repeat the MOPAC calculations first within Sybyl.
The .dat file for the last FORCE calculation
 AM1 CHARGE=-1 T=8000 FORCE NOINTER NOXYZ
 sn2 reaction, location of ts
 
 Br    0.00000000  0    0.0000000  0    0.0000000  0    0    0    0     -0.7504
  C    2.44997974  1    0.0000000  0    0.0000000  0    1    0    0      0.1841
  H    1.10095987  1   73.6994822  1    0.0000000  0    2    1    0      0.1649
  C    1.49462660  1   95.6968652  1  121.2687668  1    2    1    3     -0.2039
  C    1.49462171  1   95.7511836  1 -121.2379141  1    2    1    3     -0.2039
 Cl    2.19290408  1  156.5867701  1    0.0000000  0    2    1    3     -0.6776
  H    1.11862609  1  110.8750551  1   59.8135012  1    5    2    4      0.0846
  H    1.11650825  1  110.4026520  1 -179.7385318  1    5    2    4      0.0808
  H    1.11937240  1  110.6240975  1  -59.6531514  1    5    2    4      0.0780
  H    1.11937154  1  110.6527622  1   60.3571905  1    4    2    5      0.0784
  H    1.11654691  1  110.3952703  1 -179.5465968  1    4    2    5      0.0809
  H    1.11864270  1  110.8480230  1  -59.1139105  1    4    2    5      0.0842

contains the correct geometry; it can be transformed into a Sybyl .mol2 file using Babel (_C_ comp. chem. menu).
Since we use the internal MOPAC, an earlier version, we repeat the TS calculation. This can be done through the menu, Compute, Interfaces, MOPAC, select optimization, make the charge -1 and add the keyword TS. Also select the medium queue and submit.
After completion, check the new geometry for anomalies and, if correct, enter the following commands (in which snvib becomes the filename and the > represents the sybyl prompt. Some answers are presented as defaults):
>qcpe m1 mopac setup snvib
>charge
>-1
>method
>am1
>output
>force
>(yes)
and 'return's to leave setup mode and submit the job under the name snvib.
After completion, enter:
>qcpe m1 mopac retrieve snvib
>normal_modes
>animation
upon which Sybyl reads the vibration information from the .out file, which is also displayed on the screen. Mark the value for the first vibration (which should be negative), as Sybyl prompts for the frequency to be animated.
Nine geometries are calculated, and shown alternatingly.
Translated to an MPEG file it looks like this.
This is a side step from paragraph 2A.