Molecular orbitals

Molecuar or crystal Kohn-Sham orbitals can be output in the Gaussian cube format, and thereby visualized by many software such as VESTA [103] and XCrySDen [105]. The relevant keywords are given by

     MO.fileout         on    # on|off, default=off
     num.HOMOs           1    # default=2
     num.LUMOs           1    # default=2
     MO.Nkpoint          2    # default=1

     <MO.kpoint
       0.0  0.0  0.0
       0.5  0.0  0.0
     MO.kpoint>
When you want to generate the cube files, please swtich on the keyword 'MO.fileout'. The numbers of the highest occupied molecular (crystal) orbitals (HOMOs) and the lowest occupied molecular (crystal) orbitals (LUMOs) to be output can be specified by the keywords 'num.HOMOs' and 'num.LUMOs', respectively. In case of a band calculation, the k-points at which the HOMOs and LUMOs are calculated are specified by the keywords: 'MO.Nkpoint' and 'MO.kpoint'. The keyword 'MO.Nkpoint' gives the number of k-points at which the HOMOs and LUMOs are calculated, and by the keyword 'MO.kpoint' you can specifiy the k-points explicitly as shown above, where the specification is made based on the reciprocal vectors for the unit cell vectors given by 'Atoms.UnitVectors'. The output cube files are summarized as below:

Cluster cases:

If 'MO.fileout=ON' and 'scf.EigenvalueSolver=Cluster', the following files are also generated:

Bulk cases:

If 'MO.fileout=ON' and 'scf.EigenvalueSolver=Band', the following files are also generated:

As an example, Fig. 31 show the HOMO and LUMO of a valorphin molecule.

Figure 31: (a) HOMO of a valorphin molecule, and (b) LUMO of a valorphin molecule. The figures were visualized by VESTA [103]. The input file is 'Valorphin_MO.dat' in the directory 'work'.
\includegraphics[width=15.0cm]{Val_MOs.eps}