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)