Constrained nuclear-electronic orbital density functional theory: Energy surfaces with nuclear quantum effects

The nuclear-electronic orbital (NEO) framework enables the incorporation of nuclear quantum effects by treating both electrons and specific key nuclei quantum-mechanically. The conventional NEO method predicates on the controversial Born-Oppenheimer separation between classical and quantum nuclei, and its potential energy surface only depends on the coordinates of classical nuclei. In this paper, based on the fact that quantum nuclei are relatively localized, we develop the constrained nuclear-electronic orbital density functional theory (cNEO-DFT) by imposing a constraint on the expectation value of the quantum nuclear position. In this way, an extended NEO energy surface is obtained, which also depends on the quantum nuclear position. Compared to the potential energy surface obtained from conventional DFT, the extended NEO energy surface incorporates the nuclear quantum effects, which have notable impacts on the energy profile. Furthermore, cNEO-DFT can facilitate the location of NEO stationary states. It potentially can be used in geometry optimization, transition states search, and the calculation of reaction dynamics.