I have run an NVT AIMD spin-polarised calculation on an isolated O2 molecule in a box. I would like to perform motional averaging of the O2 spin density. My understanding was that one could extract the spin desnity frame by frame form the AIMD trajectory, convert the spin density to the cube file format and perform the averaging of the spin density over the entire trajectory.