Molecular geometric phase from the exact electron-nuclear factorization

Abstract

The Born-Oppenheimer electronic wave function ${\mathrm{\ensuremath{\Phi}}}_{R}^{\text{BO}}(r)$ picks up a topological phase factor $\ifmmode\pm\else\textpm\fi{}1$, a special case of Berry phase, when it is transported around a conical intersection of two adiabatic potential energy surfaces in $R$ space. We show that this topological quantity reverts to a geometric quantity ${e}^{i\ensuremath{\gamma}}$ if the geometric phase $\ensuremath{\gamma}=\ensuremath{\oint}\mathrm{Im}\ensuremath{\langle}{\mathrm{\ensuremath{\Phi}}}_{R}|{\mathbf{\ensuremath{\nabla}}}_{\ensuremath{\mu}}{\mathrm{\ensuremath{\Phi}}}_{R}\ensuremath{\rangle}\ifmmode\cdot\else\textperiodcentered\fi{}d{\mathbf{R}}_{\ensuremath{\mu}}$ is evaluated with the conditional electronic wave function ${\mathrm{\ensuremath{\Phi}}}_{R}(r)$ from the exact electron-nuclear factorization ${\mathrm{\ensuremath{\Phi}}}_{R}(r)\ensuremath{\chi}(R)$ instead of the adiabatic function ${\mathrm{\ensuremath{\Phi}}}_{R}^{\text{BO}}(r)$. A model of a pseudorotating triatomic molecule, also applicable to dynamical Jahn-Teller ions in bulk crystals, provides examples of nontrivial induced vector potentials and molecular geometric phase from the exact factorization. The induced vector potential gives a contribution to the circulating nuclear current that cannot be removed by a gauge transformation. The exact potential energy surface is calculated and found to contain a term depending on the Fubini-Study metric for the conditional electronic wave function.