Quantized tensor networks for solving the Vlasov-Maxwell equations

Erika Ye,Nuno Loureiro
2024-06-04
Abstract:The Vlasov-Maxwell equations provide an \textit{ab-initio} description of collisionless plasmas, but solving them is often impractical because of the wide range of spatial and temporal scales that must be resolved and the high dimensionality of the problem. In this work, we present a quantum-inspired semi-implicit Vlasov-Maxwell solver that utilizes the quantized tensor network (QTN) framework. With this QTN solver, the cost of grid-based numerical simulation of size $N$ is reduced from $\mathcal{O}(N)$ to $\mathcal{O}(\text{poly}(D))$, where $D$ is the ``rank'' or ``bond dimension'' of the QTN and is typically set to be much smaller than $N$. We find that for the five-dimensional test problems considered here, a modest $D=64$ appears to be sufficient for capturing the expected physics despite the simulations using a total of $N=2^{36}$ grid points, \edit{which would require $D=2^{18}$ for full-rank calculations}. Additionally, we observe that a QTN time evolution scheme based on the Dirac-Frenkel variational principle allows one to use somewhat larger time steps than prescribed by the Courant-Friedrichs-Lewy (CFL) constraint. As such, this work demonstrates that the QTN format is a promising means of approximately solving the Vlasov-Maxwell equations with significantly reduced cost.
Computational Physics,Plasma Physics,Quantum Physics
What problem does this paper attempt to address?
This paper attempts to solve the computational complexity and high - dimensional problems in solving the Vlasov - Maxwell equations. Specifically, the Vlasov - Maxwell equations provide an ab - initio description of the behavior of collisionless plasmas, but due to the wide range of spatio - temporal scales to be resolved and the highly multi - dimensional nature of the problem, direct solution of these equations is usually computationally impractical. ### Main problems and solutions in the paper 1. **Computational complexity**: - Solving the Vlasov - Maxwell equations involves a wide range of spatio - temporal scales and high dimensions (for example, 6 + 1 - dimensional PDE), which makes traditional numerical simulation methods very costly in terms of computational cost. - The author proposes a quantum - inspired semi - implicit Vlasov - Maxwell solver, using the Quantized Tensor Networks (QTN) framework. This method reduces the cost of grid - based numerical simulation from \(O(N)\) to \(O(\text{poly}(D))\), where \(D\) is the "rank" or "bond dimension" of QTN, which is usually much smaller than \(N\). 2. **High - dimensional problem**: - To deal with the high - dimensional problem, the author uses the QTN format to represent and manipulate data. In this way, even for a five - dimensional test problem, using a relatively small \(D = 64\) can capture the expected physical phenomena, without the need to use a huge number of grid points (such as \(N = 2^{36}\)) as in full - rank calculations. 3. **Time - step limitation**: - Traditional grid - based methods are limited by the Courant - Friedrichs - Lewy (CFL) condition and must use very small time steps. The QTN time - evolution scheme based on the Dirac - Frenkel variational principle allows the use of larger time steps, thereby improving computational efficiency. ### Key contributions - **Improvement in computational efficiency**: By introducing the QTN framework, the computational resources required to solve the Vlasov - Maxwell equations are significantly reduced. - **Feasibility of large - scale simulations**: It makes plasma simulations with higher dimensions and finer resolutions possible, which is crucial for understanding astrophysical phenomena and fusion energy systems. - **Flexibility in time steps**: A new time - evolution scheme is proposed, breaking through the time - step limitations of traditional methods and further improving computational efficiency. ### Formula explanations The formulas involved in the paper are accurately presented in Markdown format, for example: - Vlasov equation: \[ \frac{\partial f_s}{\partial t}+v_s\cdot\nabla_x f_s+\frac{1}{m_s}F_s\cdot\nabla_v f_s = 0 \] - Lorentz force: \[ F_s=q_s(v_s\times B + E) \] - Maxwell equations: \[ \frac{\partial B}{\partial t}=-\nabla\times E \] \[ \frac{\partial E}{\partial t}=c^2\nabla\times B-\frac{1}{\varepsilon_0}J \] Through these improvements, the paper shows that the QTN format is a promising method for approximately solving the Vlasov - Maxwell equations with a significant reduction in computational cost.