I would like to calculate both the rotational constant B and the centrifugal distortion constant D of HBr (for start) using Gaussian 03. Does anyone know what the proper command for this task is?
Try first reoptimizing the structure with Opt(tight) SCF(Tight) NoSymm keywords, since your system is apparently not in the minimum. The run the VibRot job as a separate one.