Ilya Vilkoviskiy,Dmitry A. Abanin
Abstract:Spin-boson (SB) model plays a central role in studies of dissipative quantum dynamics, both due its conceptual importance and relevance to a number of physical systems. Here we provide rigorous bounds of the computational complexity of the SB model for the physically relevant case of a zero temperature Ohmic bath. We start with the description of the bosonic bath via its Feynman-Vernon influence functional (IF), which is a tensor on the space of spin's trajectories. By expanding the kernel of the IF functional via a sum of decaying exponentials, we obtain an analytical approximation of the continuous bath by a finite number of damped bosonic modes. We bound the error induced by restricting bosonic Hilbert spaces to a finite-dimensional subspace with small boson numbers, which yields an analytical form of a matrix-product state (MPS) representation of the IF. We show that the MPS bond dimension $D$ scales polynomially in the error on physical observables $\epsilon$, as well as in the evolution time $T$, $D\propto T^4/\epsilon^2$. This bound indicates that the spin-boson model can be efficiently simulated using polynomial in time computational resources.
Quantum Physics,Mesoscale and Nanoscale Physics,Strongly Correlated Electrons,Mathematical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to efficiently simulate the computational complexity of the Spin - Boson (SB) model in non - Markovian dynamics. Specifically, the paper focuses on the computational complexity of the SB model under the zero - temperature Ohmic bath condition and provides a strict boundary estimate. By introducing a finite number of damped bosonic modes to approximate the continuous bosonic bath, the author obtains an Influence Functional (IF) represented by a Matrix Product State (MPS). The main result of the paper shows that the bond dimension \(D\) of the MPS has the following relationship with time \(T\) and error \(\epsilon\):
\[ D \sim \left(\frac{\omega_c T}{\epsilon}\right)^4 \]
where \(\omega_c\) is the cutoff frequency, \(T\) is the evolution time, and \(\epsilon\) is the error of the physical observable. This result indicates that the SB model can be efficiently simulated using polynomial - time computational resources.
### Background of the Paper
The Spin - Boson model plays a central role in the study of dissipative quantum dynamics, not only because of its conceptual importance but also because it is related to many physical systems. Although rich analytical results have been obtained in some limit cases, it is still a challenge to accurately describe the dynamics of the system in the strong - coupling limit. To meet this challenge, researchers have developed a variety of numerical methods, such as chain mapping combined with tDMRG technology, truncating the Feynman - Vernon Influence Functional (IF) or enhancing the reduced - density tensor (ADT), etc.
### Main Contributions
1. **Strict Complexity Boundary**: The paper provides a strict boundary estimate of the computational complexity of the SB model under the zero - temperature Ohmic bath condition.
2. **MPS Representation**: By introducing a finite number of damped bosonic modes, the paper gives an MPS representation of the IF and proves the polynomial relationship between the bond dimension of the MPS and time and error.
3. **Error Analysis**: The paper analyzes in detail the errors introduced by truncating bosonic modes and the approximate function \(\eta_{i - j}\) and proves that these errors can be controlled.
### Method Overview
- **Analytical Approximation of the Influence Function**: By expanding the IF kernel as a sum of decaying exponentials, the author obtains an analytical approximation of the continuous bath.
- **MPS Representation**: Using auxiliary bosonic modes, the IF is represented in the form of MPS, and the bond dimension of the MPS is analyzed.
- **Error Estimation**: By analyzing the truncation of bosonic modes and the approximation error of the function \(\eta_{i - j}\), it is proved that the total error can be controlled within a given accuracy range.
### Conclusion
By introducing a finite number of damped bosonic modes, the paper successfully represents the Influence Functional of the SB model in the form of MPS and provides a strict complexity boundary estimate. This result indicates that the SB model can be efficiently simulated using polynomial - time computational resources, providing a theoretical basis for numerical calculations. In addition, the method in the paper can be extended to more general non - Gaussian Influence Functionals, which has important theoretical and practical application values.