Basis Subroutine and System Specific Parameters: Collision of an Atom with a CH₂(X ³B₁) (0,v₂,0) Bender Vibrational Level

System subroutine: sych2x

Basis subroutine: hibach2x.f

NTERM: the number of potential surfaces involved. This parameter should equal the number of MU terms in the angular expansion of the potential.. This parameter can not be changed

NUMPOT: an index representing the particular potential used. This variable is passed to the POT subroutine

IPOTSY: cylindrical symmetry of potential. The variables (theta,phi) describing the angular expansion of the potential should be defined with the a inertial axis defined as the body-frame z axis and, if possible, the xz plane as a plane of symmetry of the molecule (in this case, POTSY = 2). If the flag IHOMO = .true., only terms with LAMBDA + MU equal to an integer multiple of IPOTSY can be included in the potential. Example: for H₂O, IPOTSY = 2

IOP: ortho/para label for molecular states of the asymmetric top. If IHOMO =.true. then

if IOP = 1: only para states included in channel expansion
if IOP = −1: only ortho states included in channel expansion

IVBEND: the v₂ bending vibrational quantum number [0 <= IVBEND <= 3]

JMAX: the maximum rotational angular momentum included in the channel expansion

EMAX: the maximum energy of a state to be included in the rotational state basis

The CH₂(³B₁) state is well described by Hund's case (b) coupling, in which the electron spin [S=1] is weakly coupled to the molecular frame. As a result, the rotational rotational angular momentum n is coupled to the electron to yield states of total angular momentum j = n − 1, n, n + 1 [if n >= 1]. This subroutine sets up a calculation of spin-free cross sections, for transitions between states of differing n.

Because of the large effective A rotational constant in the bender levels, the rotational wave functions are well approximated by symmetric top wave functions. The latter can be expressed in a symmetrized basis as [S. Green, J. Chem. Phys. 64, 3463 (1976)]

| n k m s > = [2(1+δ_k0)]^−1/2 [ | n k m > + s | n −k m >]

where s is the symmetry index (+1 or −1). By setting BASTST=.TRUE., you can output the values of n, s, the values of the asymmetric top prolate and oblate projection quantum numbers [k_p and k_o], and the internal energies.

---


Return to Hibridon Help