a |ψ β ⟩ is the second-derivative coupling, which were defined in Eqs. 4 and 5, respectively. N R is the number of nuclear DOFs in the system. We refer the reader to Ref. 114 and Ref 115 for evaluating quantities (e.g., the nuclear kinetic energy) using the grid or spectral basis for the nuclear DOFs.