Incorporating memory into propagation of 1-electron reduced density matrices

For any linear system with unreduced dynamics governed by invertible propagators, we derive a closed, time-delayed, linear system for a reduced-dimensional quantity of interest. This method does not target dimensionality reduction: rather, this method helps shed light on the memory-dependence of 1-electron reduced density matrices in time-dependent configuration interaction (TDCI), a scheme to solve for the correlated dynamics of electrons in molecules. Though time-dependent density functional theory has established that the 1-electron reduced density possesses memory-dependence, the precise nature of this memory-dependence has not been understood. We derive a symmetry/constraint-preserving method to propagate reduced TDCI electron density matrices. In numerical tests on two model systems (H2 and HeH+), we show that with sufficiently large time-delay (or memory-dependence), our method propagates reduced TDCI density matrices with high quantitative accuracy. We study the dependence of our results on time step and basis set. To implement our method, we derive the 4-index tensor that relates reduced and full TDCI density matrices. Our derivation applies to any TDCI system, regardless of basis set, number of electrons, or choice of Slater determinants in the wave function.