A method for computing second-order multiconfigurational perturbation theory (CASPT2) energy gradients numerically has been implemented and applied to a range of elementary organic chromophores, including 1,3 butadiene, acrolein, and two protonated Schiff bases. Geometries of ground and excited states-as well as conical intersections-are compared with the corresponding CASSCF structures, illustrating the effect of including the correction for dynamical electron correlation. It is shown that the differences between the two methods are not readily categorized, but that, while individual changes in bond lengths can be quite large ( approximately 0.01-0.02 A), the natures and CASPT2 energetics of the structures remain similar. Exceptions to this tend to be systems that have a strong ionic character and that are not well described at the CASSCF level. Basis set effects (double- vs. triple-zeta) were examined for a limited number of examples, and found to be quite dramatic at both levels of theory.