The initial distribution {R, P} is then obtained by transforming { R, P} from the normal mode representation to the primitive coordinates using the unitary transformation that diagonalizes the Hessian matrix. The positions and momenta for all molecules except for fulvene were sampled from the above Wigner distribution. For fulvene, only the coordinates were sampled and the momenta were set to zero, such that the molecule encounters only the slanted conical intersection.