In this work, we introduced a brand-new interface between PySpawn, software for running nonadiabatic molecular dynamics with AIMS algorithm, and OpenMolcas, a package for running multireference electronic structure calculations. This interface allows the user to run AIMS/CASSCF and AIMS/CASPT2 dynamics at different flavors: MS-CASPT2, XMS-CASPT2, XDW-CASPT, and RMS-CASPT2. It exploits the best of the two packages: a Python, modular implementation of the AIMS algorithm and extensive analysis capability in PySpawn, and advanced implementation of electronic structure calculations in OpenMolcas, including the availability of analytical gradients at the CASPT2 level, crucial for an affordable dynamics at such an accurate, but costly, level of theory. We applied our interface to fulvene, a challenging and interesting system. After being excited to the first singlet electronic state, fulvene can decay to the ground state through two different channels, promoted by two distinct crossings between the S 0 and S 1 surfaces. The perturbational correction to the energy, changes importantly the potential energy surfaces, closing one of the deactivation channels and inducing a slower relaxation with less degree of the characteristic reflection from S 0 to S 1 . Additionally, the four CASPT2 flavors do not induce the same dynamics, but slight differences can be observed according to the different implementations. XMS-CASPT2 and RMS-CASPT2 are the two flavors that perform the best in nonadiabatic dynamics with the energies calculated on-the-fly, ensuring smooth potentials and energy conservation. The applications and comparison of CASPT2 flavors' performances in nonadiabatic dynamics is an area yet to be fully explored, and in this work, we gave an idea of its importance by pointing out the dependence of the simulated behaviors on the level of structure employed.