Adaptive Kink Filtration: Achieving Asymptotic Size-Independence of Path Integral Simulations Utilizing the Locality of Interactions

Amartya Bose

2024-09-23

Abstract:Recent method developments involving path integral simulations have come a long way in making these techniques practical for studying condensed phase non-equilibrium phenomena. The main difficulty that still needs to be surmounted is that the scaling of the algorithms still depends upon the system dimensionality. The majority of recent techniques have only changed the order of this scaling and not eased the dependence on the system size. In this current work, we introduce an adaptive kink filtration technique for path generation approach that leverages the locality of the interactions present in the system and the consequent sparsity of the propagator matrix to remove the asymptotic size dependence of the simulations for the propagation of reduced density matrices. This enables the simulation of larger systems at a significantly reduced cost. The only place where a system size-dependence still remains is the calculation of the dynamical maps or propagators. The cost of calculating this solvent-renormalized propagator is the same as the cost of propagating all the elements of the reduced density matrix, which scales as the square of the size. The adaptive kink filtration technique introduced here is compatible with previously introduced approaches to making path integral simulations more efficient, and can additionally be used for calculation of correlation functions.

Chemical Physics