The user documentation for Hibridon
BASISTYPE = 9
System subroutine: systpln
Basis subroutine: hiba09_stpln.F90
Ref. C. Rist, M. H. Alexander, and P. Valiron, J. Chem. Phys. 98, 4662 (1993).
NTERM: (the number of potential surfaces involved: this should be 2 in this case). This parameter can not be changedNUMPOT: an index representing the particular potential used. This variable is passed to the POT subroutineIPOTSY: cylindrical symmetry of potential. If the flag IHOMO = .True., only terms with mu1 equal to an integer multiple of IPOTSY can be included in the potential. Example: for NH3, IPOTSY = 3IOP: ortho/para label for rotational states of diatomic. If IHOMO = .True. then
IOP = 1: only para states are included in channel expansionIOP = -1: only ortho states are included in channel expansionNINV: number of inversion doublets included:
NINV = +1, only + inversion levels includedNINV = -1, only - inversion levels includedNINV = 2, both inversion levels includedKMAX: the maximum projection quantum number includedJMAX0: the maximum rotational angular momenta or the k = 0 stackJMAX1: the maximum rotational angular momenta or the k = 1 stackJMAX2: the maximum rotational angular momenta or the k = 2 stackJMAX3: the maximum rotational angular momenta or the k = 3 stackJMAX4: the maximum rotational angular momenta or the k = 4 stackJMAX5: the maximum rotational angular momenta or the k = 5 stackJMAX6: the maximum rotational angular momenta or the k = 6 stackIPOTSY2: symmetry of potential. if the linear molecule is homonuclear then IPOTSY = 2 and only terms with lambda2 even can be included in the potential, else IPOTSY = 1J2MAX: the maximum rotational angular momentum for the linear moleculeJ2MIN: the minimum rotational angular momentum for the linear moleculeBROT, CROT: rotational constants of symmetric topDROT: rotational constant of diatomic molecule
E(jkj2) = BROT j (j + 1) + (CROT - BROT) k2 + DROT j2 (j2 + 1)