The nuclear quadrupole interaction (NQI) of quadrupolar nuclei in molecules in the gas phase has recently been quantified experimentally by the technique of Perturbed Angular Correlation (PAC) spectroscopy. Interpretation of the NQI data must rely on high-level electronic structure methods for the calculation of the electric field gradient (EFG). These calculations commonly omit contributions from anharmonicity, rovibrational coupling and (if applicable) decay-induced nuclear recoil to the molecular structure and the response of the EFG. Here we choose an alternative approach, conducting Born-Oppenheimer molecular dynamics (BOMD) simulations, and account for the effects of molecular motion on the EFG. We apply this approach to the analysis of PAC spectroscopic data recorded for gas phase CdI2 (Haas et al., 2021). With this methodology, we obtain a value of the nuclear quadrupole moment Q(Cd-111, 5/2(+)) = 0.656(7) b.