An
algorithm for the calculation of hyperfine structure and spectra
of diatomic molecules based on the variational nuclear motion is presented.
The hyperfine coupling terms considered are Fermi-contact, nuclear
spin-electron spin dipole–dipole, nuclear spin–orbit,
nuclear spin-rotation, and nuclear electric quadrupole interactions.
Initial hyperfine-unresolved wave functions are obtained for a given
set of potential energy curves and associated couplings by a variation
solution of the nuclear-motion Schrödinger equation. Fully
hyperfine-resolved parity-conserved rovibronic Hamiltonian matrices
for a given final angular momentum,
F
, are constructed and then diagonalized to give hyperfine-resolved
energies and wave functions. Electric transition dipole moment curves
can then be used to generate a hyperfine-resolved line list by applying
rigorous selection rules. The algorithm is implemented in
Duo
, which is a general program for calculating spectra of diatomic
molecules. This approach is tested for NO and MgH, and the results
are compared to experiment and shown to be consistent with those given
by the well-used effective Hamiltonian code PGOPHER.