Molecule Types Asymmetric Tops <Prev Next>

Asymmetric Top Hamiltonian

Closed Shell Molecules

The asymmetric top Hamiltonian used is either of the A or S standard reduced forms proposed by Watson (Watson, 1977). The SReduction flag at the Molecule level determines which one, with a setting of false implying the A reduction is used. When reporting these constants it is important to specify both the reduction and the axis system used. Note that, for open shell systems with S > 0, J should be replaced with N in the discussion below, though not in the parameter names.

The rigid rotor part can be specified in two ways, independent of the reduction:

or the equivalent form (which can give better determined parameters):

where:

and:

Versions before 8.0.226 required the a and z axes to be same, i.e. for a Ir or Il representation to be used.

The different forms of the centrifugal distortion terms do not have a simple relationship:

 

A reduction

Quartic

Sextic

Octic


 

S reduction

Quartic

Sextic

Octic

Note that various symbols are in use for many of the centrifugal distortion constants; in particular different symbols are typically used for the A and S reduction constants rather than as here. Some of these are listed below, together with the name PGOPHER uses.

PGOPHER

A reduction

S reduction

A

A

B

B

C

C

BBar

Bdelta

DJ

ΔJ

DJ

DJK  

ΔJK

DJK

DK

ΔK

DK

deltaJ

δJ

d1

deltaK

δK

d2

HJ

ΦJ

HJ

HJK

ΦJK

HJK

HKJ

ΦKJ

HKJ

HK

ΦK

HK

phiJ

φJ

h1

phiJK

φJK

h2

phiK

φK

h3

LJ

LJ

LJJK

LJJK

LJK

LJK

LKKJ

LKKJ

LK

LK

llJ

lJ

llJK

lJK

llKJ

lKJ

llK

lK

Open Shell Molecules

For open shell molecules (S > 0) a Hund’s case (b) basis is used. The rotational part is as above, though for J and its components read N throughout. The following additional terms are present:

Spin Rotation Interaction

εaaNaSa +  εbbNbSb +  εccNcSc
+ ½ε‾ab(NaSb + SbNa + NbSa + SaNb)
+ ½ε‾ac(NaSc + ScNa + NcSa + SaNc)
+ ½ε‾bc(NbSc + ScNb + NcSb + SaNb)

where ε‾ab = eabbar = ½(εab + εba), ε‾ac = eacbar = ½(εac + εca) and ε‾bc = ebcbar = ½(εbc + εcb). The slightly long winded form arises because the operators Np and Sq only commute if p = q. Following Brown and Sears (1979) the centrifugal terms come in two slightly different forms, depending on the representation:

A reduction: ΔsNN2N.S + ½ΔsNK(N2NzSz+NzSzN2) + ΔsKNN.SNz2 + ΔsKNz3Sz
+ δsNN.S(N+2 + N-2) + ½δsK{(N+2 + N-2)NzSz + NzSz(N+2 + N-2)}
S reduction: DsNN2N.S + ½DsNK(N2NzSz+NzSzN2) + DsKNN.SNz2 + DsKNz3Sz
+ ds1N.S(N+2 + N-2) + ds2(N+3S+ + N-3S-)

Only that last term actually has a different operator form. Note that the S reduction names are used to label the parameters. 

Spin Spin Interaction

α(3Sz2-S2) + β(Sx2-Sy2) = α(3Sz2-S2) + ½β(S+2+S-2)

Also in use are D = 3α and E = β.

References