The Peierls-Boltzmann transport equation (PBE), which governs nonequilibrium phonon transport, suffers from the curse of dimensionality because of its high-dimensional phase space including both real and modal spaces. We explore the use of matrix product states (MPS) for numerical simulation of the PBE. We show that an MPS configuration based on scattering events combined with a dimensionless form of the solution can drastically increase the locality of correlations between tensors in the MPS representation, enhancing its effectiveness in dimension reduction. We further examine the effects of index ordering in an MPS and find that the highest locality is achieved when tensor chains associated with both real and modal spaces are connected at the coarsest grid to each other in the center of the MPS. Using this optimal configuration and a solver inspired by the density matrix renormalization group, we solve the PBE discretized by a finite volume method (FVM). The solution is obtained for crystalline silicon across ballistic, quasiballistic, and diffusive transport regimes. An MPS truncated to a compression ratio of ${10}^{-3}$ suffices to reproduce reference solutions with high fidelity. The computational cost scales sublinearly with the number of grid points in both real and modal spaces, achieving roughly an order of magnitude reduction in computational time compared to the FVM with sparse matrix operation.

Close

PrevNext

View figure in article

• Phonons
• Thermal transport
• Transport phenomena
• Matrix product states