Linked Product Approximation to the Global Electronic Overlap Matrix

The global many-electron wave function overlap matrix accounts for all effects beyond the Born-Oppenheimer approximation in the discrete variable local diabatic representation, a numerically exact framework for modeling nonadiabatic conical intersection wave packet dynamics. Nevertheless, calculating the electronic overlap matrix from the electronic structure is computationally expensive. Here, we introduce an approximation for constructing the electronic overlap matrix between any two long-range geometries by the product of nearest-neighbor overlap matrices (links) along a path connecting these two geometries. This approximation significantly reduces the computational effort by requiring electronic structure calculations only for the nearest-neighbor overlap matrices. The accuracy of this approximation is demonstrated through an exact simulation of a proton-coupled electron transfer model. Our results show that although the approximate overlap matrix can exhibit noticeable differences from the exact ones, the conical intersection dynamics is in almost exact agreement with those using the exact overlap matrix.