Discretizing Helmholtz problems via finite elements yields linear systems whose efficient solution remains a major challenge for classical computation. In this paper, we investigate how variational quantum algorithms could address this challenge. We first show that, for regular meshes, a block encoding of the operators A and A†A arising from the high-order finite element discretisation of Helmholtz problems can be designed, resulting in a quantum circuit of depth O(p3polylog(Np)) with N the number of elements and p the order of the finite elements. Then we apply our algorithm to a one-dimensional Helmholtz problem with Dirichlet and Neumann boundary conditions for various wavenumbers.
@article{arxiv.2512.22665,
title = {Variational quantum algorithm for solving Helmholtz problems with high order finite elements},
author = {Arnaud Rémi and François Damanet and Christophe Geuzaine},
journal= {arXiv preprint arXiv:2512.22665},
year = {2025}
}