First, the accuracy of the implemented gradient and the smoothness of the potential energy surface (PES) are examined using acrolein in water. The reference wavefunctions were obtained at the three-state-averaged (SA3)-CASSCF(6e,5o) level, where the active space consisted of six electrons in five orbitals. CASPT2 calculations were carried out with an imaginary level shift of 0.2i.