4 Induced Current Density Programs

4.10 TIPOMO - orbital selection

The package provides the program (TIPOMO) to select the π/σ orbitals to be used for the calculation of the contributions to the various kind of density functions discussed before. The current version of the program works well in case of molecules containing second-row atoms and hydrogen atoms. Future releases will fix eventual limitations. When the program is run without argument, a fully conjugated π-system is assumed. Alternatively, one can issue the number of double bonds as an argument, i.e., TIPOMO nπ. The search of π-orbitals goes as described in the following.

  • Count how many Hydrogen atoms are in the molecules (nH).

  • Calculate the number of core orbitals ncore=NATOnH, where NATO is the total number of atoms.

  • When run without argument, set nπ=ncore/2, which is the correct value in case of fully conjugated π-systems.

  • For each occupied molecular orbitals, calculate the maximum value of the electron density on the nuclei and sorts these densities by descending order. The first ncore molecular orbitals in the ordered list are designed to be core orbitals.

  • For each of the remaining orbitals, calculate the maximum value of the electron density for each pair of bonded atoms in the bond center position and sorts these densities by ascending order. The first nπ are designed to be π-orbitals.

The program can be used in four different ways:

TIPOMO    : write PIG file with pi-orbital indices
TIPOMO -  : write SIG file with sigma-orbital indices
TIPOMO  n : write PIG file with n indices
TIPOMO -n : write SIG file with n indices

Both PIG and SIG files can be used by the IC programs by means of the -qf option. Note that these files contain n+1 lines, the first with the value of n and the others with the molecular orbital indices (one per line) designed to be π or σ orbitals.