Solving the Nonlinear Vlasov Equation on a Quantum Computer

We present a mapping of the nonlinear, electrostatic Vlasov equation with Krook-type collision operators, discretized on a (1+1) dimensional grid, onto a recent Carleman linearization-based quantum algorithm for solving ordinary differential equations (ODEs) with quadratic nonlinearities. We derive upper bounds for the query- and gate complexities of the quantum algorithm in the limit of large grid sizes. We conclude that these are polynomially larger than the time complexity of the corresponding classical algorithms. We find that this is mostly due to the dimension, sparsity and norm of the Carleman linearized evolution matrix. We show that the convergence criteria of the quantum algorithm places severe restrictions on potential applications. This is due to the high level of dissipation required for convergence, that far exceeds the physical dissipation effect provided by the Krook operator for typical plasma physics applications.