Incorporating Memory into Propagation of 1-Electron Reduced Density Matrices
<|reference_start|>Incorporating Memory into Propagation of 1-Electron Reduced Density Matrices: For any linear system where the unreduced dynamics are governed by unitary propagators, we derive a closed, time-delayed, linear system for a reduced-dimensional quantity of interest. We apply this method to understand 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 self-contained, symmetry/constraint-preserving method to propagate reduced TDCI electron density matrices. Our method preserves properties of density matrices such as Hermitian symmetry and constant trace. In numerical tests on two model systems ($\text{H}_2$ and $\text{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. This derivation enables a proof that the trace of the reduced TDCI density matrix is constant and equals the number of electrons.<|reference_end|>
@article{bhat2024incorporating,
    title={Incorporating Memory into Propagation of 1-Electron Reduced Density
  Matrices},
    author={Harish S. Bhat and Hardeep Bassi and Karnamohit Ranka and Christine M.
  Isborn},
    journal={arXiv preprint arXiv:2403.15596},
    year={2024},
    archivePrefix={arXiv},
    eprint={2403.15596},
    primaryClass={math.DS cs.CE math-ph math.MP physics.chem-ph quant-ph}
}