Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2608.10671

Quantum Computing for Industrial Electromagnetics: Applicability and Case Studies in Solving Maxwell's Equations

Francesco Turro, Marco Maronese, Daniele Dragoni

quant-phSubmitted Aug 11, 2026• Announced Aug 12, 2026• 0 citations
arXivHTMLPDF

Abstract

Computational electromagnetics plays a central role in many industrial applications but often requires substantial computational resources, particularly when fine spatial discretizations are needed. While classical approaches remain the standard, quantum computing offers the potential to accelerate large-scale simulations by encoding them with a limited number of qubits. Here, we investigate the performance and resource scaling of the Harrow-Hassidim-Lloyd (HHL) and Quantum Singular Value Transformation (QSVT) algorithms for solving linear systems generated by the finite-difference time-domain (FDTD) method, a widely adopted numerical scheme for discretizing Maxwell's equations. We benchmark their performance across representative industrial use cases, including radar propagation, lens simulations, and beamforming processes. Our results demonstrate the validity of the approaches, achieving state infidelities smaller than $2\cdot 10^{-3}$ with success probabilities greater than $10^{-3}$, compatible with practical quantum state sampling. Overall, we observe that the QSVT method consistently delivers higher accuracy. We further observe that the condition number of the linear matrix, a key factor governing the performance of quantum solvers, saturates as the number of spatial lattice points increases. This implies that the spatial grid can be scaled to realistic industrial dimensions without increasing the HHL or QSVT circuit depth due to ill-conditioned matrices.

Comments
19 pages, 13 figures

Related papers

Quantum-Based Methods for Solving Partial Differential Equations in Electromagnetics Engineering

Emanuel Colella, Luca Bastianelli, Valter Mariani Primiani +2 more

Quantum methods for solving partial differential equations (PDEs) generally fall into a few general categories, each leveraging different quantum mechanisms to achieve potential speedups. One well-known approach is the quantum linear system solver, epitomized by the Harrow-Hassidim-Lloyd (HHL) algorithm, which uses phase estimation and eigenvalue inversion to solve large, sparse linear systems arising from discretized linear PDEs [1]. This method can yield exponential improvements over classical solvers under certain conditions (mostly when the matrix is well-conditioned and sparse). More recent developments, such as Quantum Singular Value Transformation (QSVT), provide a framework for matrix inversion in a more flexible way, often improving on the cost and error dependence of earlier solvers [2]. Time-dependent PDEs can be tackled via Hamiltonian simulations. The Trotter-Suzuki method decomposes the exponential of a Hamiltonian into products of exponentials of simpler sub-Hamiltonians, allowing a piecewise approximation of continuous dynamics [3]. Continuous-time quantum walks, which treat a PDE operator as a Hamiltonian that governs the walk, similarly simulate processes like diffusion or wave propagation on discretized grids or graphs [4]. Beyond direct simulation, variational quantum algorithms have emerged as a promising route in the noisy intermediate-scale quantum era. These approaches encode the PDE solution as a parameterized quantum state and optimize the parameters via classical feedback, thereby minimizing a cost function tied to the PDE residual or energy functional [5]. Quantum annealing or adiabatic methods, meanwhile, can solve optimization problems derived from PDE constraints by adiabatically evolving a simple initial Hamiltonian to the final problem Hamiltonian and staying in the system's ground state [3]. In our study, we implement a new parameterized quantum model, specifically designed to solve large linear systems on quantum hardware. This model is less deep and computationally faster than the proposed one in [5], reducing the complexity and speeding up convergence times. Such improvements are particularly valuable in cases where classical algorithms may fail due to the problem's size or complexity. Ongoing research is focusing on generalized quantum models that can be adapted to practical PDE problems for solving large electromagnetic problems where classical computers can fail.

USNC-URSI Radio Science Meeting•Jul 13, 2025•0 citations

Practical Application of the Quantum Carleman Lattice Boltzmann Method in Industrial CFD Simulations

Francesco Turro, Alessandra Lignarolo, Dragoni, Daniele

Computational Fluid Dynamics simulations are crucial in industrial applications but require extensive computational resources, particularly for extreme turbulent regimes. While classical digital approaches remain the standard, quantum computing promises a breakthrough by enabling a more efficient encoding of large-scale simulations with a limited number of qubits. This work presents a practical numerical assessment of a hybrid quantum-classical approach to CFD based on the Lattice Boltzmann Method (LBM). The inherently non-linear LBM equations are linearized via a Carleman expansion and solved using the quantum Harrow Hassidim Lloyd algorithm (HHL). We evaluate this method on three benchmark cases featuring different boundary conditions, periodic, bounceback, and moving wall, using statevector emulation on high-performance computing resources. Our results confirm the validity of the approach, achieving median error fidelities on the order of $10^{-3}$ and success probabilities sufficient for practical quantum state sampling. Notably, the spectral properties of small lattice systems closely approximate those of larger ones, suggesting a pathway to mitigate one of HHL's bottlenecks: eigenvalue pre-evaluation.

arXiv:2504.13033•Apr 17, 2025•0 citations

Quantum Method for Finite Element Simulation of Electromagnetic Problems

Jianan Zhang, Feng Feng, Qi‐Jun Zhang

Quantum simulation of electromagnetic (EM) structures is still in its infancy. In this paper, we investigate the possibility of applying quantum computing to solve the finite element equations for EM problems. Specifically, we propose to leverage a milestone in quantum computing, i.e., the Harrow/Hassidim/Lloyd (HHL) algorithm, to solve the finite element equations in the EM domain. To do this, we first present a systematic method to reformulate the finite element equations into quantum computation format and apply the HHL algorithm to solve the new equations. Taking advantage of the special properties of finite element matrices in EM problems, we then devise a stepwise algorithm to efficiently determine the suitable values for the hyperparameters of HHL. The proposed method can in theory solve the EM finite element equations in logarithmic time relative to the size of the finite element matrix. This is exponentially faster than that achievable by classical computation as the matrix size increases. A two-dimensional EM example is used to illustrate how the proposed method can be used to find the solutions to EM problems through quantum computation.

IEEE MTT-S International Microwave Symposium digest•Jun 7, 2021•13 citations

Toward Practical Application of the Quantum Carleman Lattice Boltzmann Method in Industrial CFD Simulations

Francesco Turro, Alessandra Lignarolo, Daniele Dragoni

Computational Fluid Dynamics simulations are crucial in industrial applications but require extensive computational resources, particularly for extreme turbulent regimes. While classical digital approaches remain the standard, quantum computing promises a breakthrough by enabling a more efficient encoding of large-scale simulations with a limited number of qubits. This work presents a practical numerical assessment of a hybrid quantum-classical approach to CFD based on the Lattice Boltzmann Method (LBM). The inherently non-linear LBM equations are linearized via a Carleman expansion and solved using the quantum Harrow-Hassidim-Lloyd algorithm (HHL). We evaluate this method on three benchmark cases featuring different boundary conditions—periodic, bounce-back, and moving wall—using state-vector emulation on high-performance computing resources. Our results confirm the validity of the approach, achieving median error fidelities on the order of$10^{-3}$and success probabilities sufficient for practical quantum state sampling. Notably, the spectral properties of small lattice systems closely approximate those of larger ones, suggesting a pathway to mitigate one of HHL's bottlenecks: eigenvalue pre-evaluation.

IEEE Transactions on Quantum Engineering•Jan 1, 2025•2 citations

Quantum Method of Calculating Multiport S-Parameter Magnitudes Based on FEM for Microwave Components: An HHL Framework

Xiaolong Li, Feng Feng, Jiali Zhang +5 more

Quantum computing in the electromagnetic (EM) field represents a transformative frontier with significant potential for advancement. Recent studies have demonstrated the potential of quantum computing for solving waveguide mode problems, solving finite element method (FEM) equations for EM problems, and scaling the FEM solution. Among these applications, the Harrow–Hassidim–Lloyd (HHL) algorithm is an algorithm with large quantum advantages and has attracted attention. However, existing quantum methods remain constrained by the problem that the quantum-based FEM solution only obtains electric field magnitude and lacks phase information. The lack of phase information results in the inability to solve an important characteristic of EM problems: scattering parameters (S-parameters). To address the above issue, this work proposes a novel quantum method of calculating multiport S-parameter magnitudes based on FEM and HHL algorithm for microwave components. The proposed method introduces novel mathematical formulations optimized for quantum architectures. To streamline the quantum computation process and avoid the need for multiple quantum circuit simulations for each S-parameter, we design an integrated quantum circuit capable of computing magnitudes of all multiport S-parameters within a single execution. Furthermore, a specific matrix-scaling method is also developed to simplify the proposed circuits. The specific matrix-scaling method is used to improve the condition of the FEM matrix and reduce the complexity of the designed quantum circuits. The proposed method is validated using two microwave structures, demonstrating how our proposed quantum method computes S-parameter magnitudes in practical EM applications.

IEEE Transactions on Microwave Theory and Techniques•Nov 5, 2025•4 citations

Quantum linear system solvers: A survey of algorithms and applications

Mauro E. S. Morales, Lirandë Pira, Philipp Schleich +7 more

Solving linear systems of equations plays a fundamental role in numerous computational problems from different fields of science. The widespread use of numerical methods to solve these systems motivates investigating the feasibility of solving linear systems problems using quantum computers. In this work, we provide a survey of the main advances in quantum linear systems algorithms, together with some applications. We summarize and analyze the main ideas behind some of the algorithms for the quantum linear systems problem in the literature. The analysis begins by examining the Harrow-Hassidim-Lloyd (HHL) solver. We note its limitations and reliance on computationally expensive quantum methods, then highlight subsequent research efforts which aimed to address these limitations and optimize runtime efficiency and precision via various paradigms. We focus in particular on the post-HHL enhancements which have paved the way towards optimal lower bounds with respect to error tolerance and condition number. By doing so, we propose a taxonomy that categorizes these studies. Furthermore, by contextualizing these developments within the broader landscape of quantum computing, we explore the foundational work that have inspired and informed their development, as well as subsequent refinements. Finally, we discuss the potential applications of these algorithms in differential equations, quantum machine learning, and many-body physics.

arXiv:2411.02522•Reviews of Modern Physics•Jan 8, 2026•4 citations
Computational electromagnetics plays a central role in many industrial applications but often requires substantial computational resources, particularly when fine spatial discretizations are needed. While classical approaches remain the standard, quantum computing offers the potential to accelerate large-scale simulations by encoding them with a limited number of qubits.
Here, we investigate the performance and resource scaling of the Harrow-Hassidim-Lloyd (HHL) and Quantum Singular Value Transformation (QSVT) algorithms for solving linear systems generated by the finite-difference time-domain (FDTD) method, a widely adopted numerical scheme for discretizing Maxwell's equations. We benchmark their performance across representative industrial use cases, including radar propagation, lens simulations, and beamforming processes. Our results demonstrate the validity of the approaches, achieving state infidelities smaller than 2⋅10−32\cdot 10^{-3}2⋅10−3 with success probabilities greater than 10−310^{-3}10−3, compatible with practical quantum state sampling. Overall, we observe that the QSVT method consistently delivers higher accuracy.
We further observe that the condition number of the linear matrix, a key factor governing the performance of quantum solvers, saturates as the number of spatial lattice points increases. This implies that the spatial grid can be scaled to realistic industrial dimensions without increasing the HHL or QSVT circuit depth due to ill-conditioned matrices.
Quantum methods for solving partial differential equations (PDEs) generally fall into a few general categories, each leveraging different quantum mechanisms to achieve potential speedups. One well-known approach is the quantum linear system solver, epitomized by the Harrow-Hassidim-Lloyd (HHL) algorithm, which uses phase estimation and eigenvalue inversion to solve large, sparse linear systems arising from discretized linear PDEs [1]. This method can yield exponential improvements over classical solvers under certain conditions (mostly when the matrix is well-conditioned and sparse). More recent developments, such as Quantum Singular Value Transformation (QSVT), provide a framework for matrix inversion in a more flexible way, often improving on the cost and error dependence of earlier solvers [2]. Time-dependent PDEs can be tackled via Hamiltonian simulations. The Trotter-Suzuki method decomposes the exponential of a Hamiltonian into products of exponentials of simpler sub-Hamiltonians, allowing a piecewise approximation of continuous dynamics [3]. Continuous-time quantum walks, which treat a PDE operator as a Hamiltonian that governs the walk, similarly simulate processes like diffusion or wave propagation on discretized grids or graphs [4]. Beyond direct simulation, variational quantum algorithms have emerged as a promising route in the noisy intermediate-scale quantum era. These approaches encode the PDE solution as a parameterized quantum state and optimize the parameters via classical feedback, thereby minimizing a cost function tied to the PDE residual or energy functional [5]. Quantum annealing or adiabatic methods, meanwhile, can solve optimization problems derived from PDE constraints by adiabatically evolving a simple initial Hamiltonian to the final problem Hamiltonian and staying in the system's ground state [3]. In our study, we implement a new parameterized quantum model, specifically designed to solve large linear systems on quantum hardware. This model is less deep and computationally faster than the proposed one in [5], reducing the complexity and speeding up convergence times. Such improvements are particularly valuable in cases where classical algorithms may fail due to the problem's size or complexity. Ongoing research is focusing on generalized quantum models that can be adapted to practical PDE problems for solving large electromagnetic problems where classical computers can fail.
Computational Fluid Dynamics simulations are crucial in industrial applications but require extensive computational resources, particularly for extreme turbulent regimes. While classical digital approaches remain the standard, quantum computing promises a breakthrough by enabling a more efficient encoding of large-scale simulations with a limited number of qubits. This work presents a practical numerical assessment of a hybrid quantum-classical approach to CFD based on the Lattice Boltzmann Method (LBM). The inherently non-linear LBM equations are linearized via a Carleman expansion and solved using the quantum Harrow Hassidim Lloyd algorithm (HHL). We evaluate this method on three benchmark cases featuring different boundary conditions, periodic, bounceback, and moving wall, using statevector emulation on high-performance computing resources. Our results confirm the validity of the approach, achieving median error fidelities on the order of 10−310^{-3}10−3 and success probabilities sufficient for practical quantum state sampling. Notably, the spectral properties of small lattice systems closely approximate those of larger ones, suggesting a pathway to mitigate one of HHL's bottlenecks: eigenvalue pre-evaluation.
Quantum simulation of electromagnetic (EM) structures is still in its infancy. In this paper, we investigate the possibility of applying quantum computing to solve the finite element equations for EM problems. Specifically, we propose to leverage a milestone in quantum computing, i.e., the Harrow/Hassidim/Lloyd (HHL) algorithm, to solve the finite element equations in the EM domain. To do this, we first present a systematic method to reformulate the finite element equations into quantum computation format and apply the HHL algorithm to solve the new equations. Taking advantage of the special properties of finite element matrices in EM problems, we then devise a stepwise algorithm to efficiently determine the suitable values for the hyperparameters of HHL. The proposed method can in theory solve the EM finite element equations in logarithmic time relative to the size of the finite element matrix. This is exponentially faster than that achievable by classical computation as the matrix size increases. A two-dimensional EM example is used to illustrate how the proposed method can be used to find the solutions to EM problems through quantum computation.
Computational Fluid Dynamics simulations are crucial in industrial applications but require extensive computational resources, particularly for extreme turbulent regimes. While classical digital approaches remain the standard, quantum computing promises a breakthrough by enabling a more efficient encoding of large-scale simulations with a limited number of qubits. This work presents a practical numerical assessment of a hybrid quantum-classical approach to CFD based on the Lattice Boltzmann Method (LBM). The inherently non-linear LBM equations are linearized via a Carleman expansion and solved using the quantum Harrow-Hassidim-Lloyd algorithm (HHL). We evaluate this method on three benchmark cases featuring different boundary conditions—periodic, bounce-back, and moving wall—using state-vector emulation on high-performance computing resources. Our results confirm the validity of the approach, achieving median error fidelities on the order of10−310^{-3}10−3and success probabilities sufficient for practical quantum state sampling. Notably, the spectral properties of small lattice systems closely approximate those of larger ones, suggesting a pathway to mitigate one of HHL's bottlenecks: eigenvalue pre-evaluation.
Quantum computing in the electromagnetic (EM) field represents a transformative frontier with significant potential for advancement. Recent studies have demonstrated the potential of quantum computing for solving waveguide mode problems, solving finite element method (FEM) equations for EM problems, and scaling the FEM solution. Among these applications, the Harrow–Hassidim–Lloyd (HHL) algorithm is an algorithm with large quantum advantages and has attracted attention. However, existing quantum methods remain constrained by the problem that the quantum-based FEM solution only obtains electric field magnitude and lacks phase information. The lack of phase information results in the inability to solve an important characteristic of EM problems: scattering parameters (S-parameters). To address the above issue, this work proposes a novel quantum method of calculating multiport S-parameter magnitudes based on FEM and HHL algorithm for microwave components. The proposed method introduces novel mathematical formulations optimized for quantum architectures. To streamline the quantum computation process and avoid the need for multiple quantum circuit simulations for each S-parameter, we design an integrated quantum circuit capable of computing magnitudes of all multiport S-parameters within a single execution. Furthermore, a specific matrix-scaling method is also developed to simplify the proposed circuits. The specific matrix-scaling method is used to improve the condition of the FEM matrix and reduce the complexity of the designed quantum circuits. The proposed method is validated using two microwave structures, demonstrating how our proposed quantum method computes S-parameter magnitudes in practical EM applications.
Solving linear systems of equations plays a fundamental role in numerous computational problems from different fields of science. The widespread use of numerical methods to solve these systems motivates investigating the feasibility of solving linear systems problems using quantum computers. In this work, we provide a survey of the main advances in quantum linear systems algorithms, together with some applications. We summarize and analyze the main ideas behind some of the algorithms for the quantum linear systems problem in the literature. The analysis begins by examining the Harrow-Hassidim-Lloyd (HHL) solver. We note its limitations and reliance on computationally expensive quantum methods, then highlight subsequent research efforts which aimed to address these limitations and optimize runtime efficiency and precision via various paradigms. We focus in particular on the post-HHL enhancements which have paved the way towards optimal lower bounds with respect to error tolerance and condition number. By doing so, we propose a taxonomy that categorizes these studies. Furthermore, by contextualizing these developments within the broader landscape of quantum computing, we explore the foundational work that have inspired and informed their development, as well as subsequent refinements. Finally, we discuss the potential applications of these algorithms in differential equations, quantum machine learning, and many-body physics.