Quantum control of molecular motion including electronic polarization effects with a two-stage toolkit

A method for incorporating strong electric field polarization effects into optimal control calculations is presented. A Born-Oppenheimer-type separation, referred to as the electric-nuclear Born-Oppenheimer (ENBO) approximation, is introduced in which variations of both the nuclear geometry and the external electric field are assumed to be slow compared with the speed at which the electronic degrees of freedom respond to these changes. This assumption permits the generation of a potential energy surface that depends not only on the relative geometry of the nuclei but also on the electric field strength and on the orientation of the molecule with respect to the electric field. The range of validity of the ENBO approximation is discussed in the paper. A two-stage toolkit implementation is presented to incorporate the polarization effects and reduce the cost of the optimal control dynamics calculations. As an illustration of the method, it is applied to optimal control of vibrational excitation in a hydrogen molecule aligned along the field direction. Ab initio configuration interaction calculations with a large orbital basis set are used to compute the H-H interaction potential in the presence of the electric field. The significant computational cost reduction afforded by the toolkit implementation is demonstrated.