Optimization of basis sets in valence bond calculations

Using the symmetry of the molecular point group, the Hamiltonian matrix is factorized. At the same time, the chemically inequivalent configurations are used to calculate the entries of the Hamiltonian matrix directly. This methodology saves a large amount of computer memory and computational time, especially for π-systems with more than 20 π-electrons.