Determining whether diabolical singularities limit the accuracy of molecular property based diabatic representations: The 1,21A states of methylamine.

An efficient, easily implemented method for locating singularities attributable to the failure of the defining equations in a molecular property based diabatization, termed diabolical singular points, is reported. For two state diabatizations, the singular points form a seam of dimension N int - 2, where N int is the number of internal degrees of freedom. The dynamical outcomes of nuclear trajectories that reach the region of this seam are flawed. The algorithm easily identifies these otherwise hard to anticipate regions of fallaciously large derivative coupling. The fact that the algorithm is easily incorporated into a two state diabatization code based on molecular properties makes it a practical tool for determining whether the existence of diabolical singularities is relevant to the problem being considered. The algorithm is illustrated using a multireference single and double excitation configuration interaction description of the 1,21A states of CH3NH2.