In this paper, the flexible model based on relaxedab initiocalculations, which has been several times employed for vibrational calculations, is extended to the analysis of the rotational structures starting by the roto-torsional bands oftrans-1,3-butadiene. For this purpose, the potential energy surface and the kinetic energy parameters of the ν13vibrational mode of butadiene are obtained with the Möller–Plesset perturbation theory up to the second order and the 6–31G(d,p), 6–31G(df,p), 6–311G(d,p), 6–311G(df,p), and 6–311G(df,pd) basis sets. The torsional levels of the −h6, −d4, and −d6isotopic species are calculated variationally and are compared with experimental data. It may be concluded that the one-dimensional model appears sufficiently accurate for butadiene−h6and −d4, whereas a large kinetic interaction with the lowest wagging mode is observed for butadiene−d6. The rotational levels corresponding to the first vibrational states of the −h6and −d4species are determined variationally up toJ= 17 andJ= 11 from theab initiospectroscopic parameters which have been expanded as functions of the torsional coordinate using symmetry adapted series. The torsional wavefunction is contracted to reduce the size of the Hamiltonian matrix. A good agreement with the observed transitions is obtained for the first states v = 0 and v = 1. As is expected, theKdoubling obtained is relatively small. For this reason, the quartic and sextic centrifugal distortion constants are obtained from the least-square fit of the variational levels to the perturbation theory equations for the symmetric top.