Simulating quantum dynamics in two-dimensional lattices with tensor network influence functional belief propagation

Describing nonequilibrium quantum dynamics remains a significant computational challenge due to the growth of spatial entanglement. The tensor network influence functional (TN-IF) approach mitigates this problem for computing the time evolution of local observables by encoding the subsystem's influence functional path integral as a matrix product state, thereby shifting the resource governing computational cost from spatial entanglement to temporal entanglement. We extend the applicability of the TN-IF method to two-dimensional lattices by demonstrating its construction on tree lattices and proposing a belief propagation (BP) algorithm for the TN-IF, termed influence functional BP (IF-BP), to simulate local observable dynamics on arbitrary graphs. Even though the BP algorithm introduces uncontrolled approximation errors on arbitrary graphs, it provides an accurate description for locally treelike lattices. Numerical simulations of the kicked Ising model on a heavy-hex lattice, motivated by a recent quantum experiment, highlight the effectiveness of the IF-BP method, which demonstrates superior performance in capturing long-time dynamics where traditional tensor network state-based methods struggle. Our results further reveal that the temporal entanglement entropy only grows logarithmically with time for this model, resulting in a polynomial computational cost for the whole method. We further construct a cluster expansion of IF-BP to introduce loop correlations beyond the BP approximation, providing a systematic correction to the IF-BP estimate. We demonstrate the power of the cluster expansion of the IF-BP in simulating the quantum quench dynamics of the two-dimensional transverse field Ising model, obtaining numerical results that improve on the state-of-the-art.