Quantum Algorithm for Solving the Wave Equation on a Hexagonal Grid
Numerical methods based on finite differences for solving the wave equation are prone to various errors, a challenge that also applies to the semi-discrete case of a recently proposed quantum algorithm that preserves continuous time while discretizing space into a quantum state vector. The inherent...
Saved in:
| Main Author: | |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
IEEE
2024-01-01
|
| Series: | IEEE Access |
| Subjects: | |
| Online Access: | https://ieeexplore.ieee.org/document/10679139/ |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1850263526516457472 |
|---|---|
| author | Roman Novak |
| author_facet | Roman Novak |
| author_sort | Roman Novak |
| collection | DOAJ |
| description | Numerical methods based on finite differences for solving the wave equation are prone to various errors, a challenge that also applies to the semi-discrete case of a recently proposed quantum algorithm that preserves continuous time while discretizing space into a quantum state vector. The inherent limitations of quantum processing, which affect the generality and extensibility of the quantum algorithm compared to classical numerical solvers, are offset by the logarithmic scaling of the number of qubits with the number of grid points. The higher grid density in combination with the higher-order approximation of the continuous Laplacian operator improves the numerical accuracy. In contrast to the Dirichlet boundary condition, however, the Neumann boundary condition in the quantum algorithm cannot exceed a second-order Laplacian approximation. Moreover, the anisotropy error is inherent to the lattice geometry. In this work, we propose a quantum wave equation solver on a hexagonal grid, which exhibits intrinsic grid dispersion comparable to that of a fourth-order accurate method on a rectangular grid and provides a quadratic improvement in the anisotropy measures. We show that the graph Laplacian on a hexagonal grid embedding the Neumann boundary condition is a bidirected graph. To make it suitable for use as a system Hamiltonian, a symmetrization and factorization based on the asymmetric adjacency matrix is proposed. Dispersion relations for quantum algorithms solving the wave equation on different lattices are derived and compared. The algorithms are analyzed as a limiting case of Yee’s finite-difference time-domain method in the context of electromagnetic wave propagation. |
| format | Article |
| id | doaj-art-41d902520dee4ee4a4ef0b4e45619423 |
| institution | OA Journals |
| issn | 2169-3536 |
| language | English |
| publishDate | 2024-01-01 |
| publisher | IEEE |
| record_format | Article |
| series | IEEE Access |
| spelling | doaj-art-41d902520dee4ee4a4ef0b4e456194232025-08-20T01:54:57ZengIEEEIEEE Access2169-35362024-01-011213100913102210.1109/ACCESS.2024.345947110679139Quantum Algorithm for Solving the Wave Equation on a Hexagonal GridRoman Novak0https://orcid.org/0000-0001-5757-6703Department of Communication Systems, Jožef Stefan Institute, Ljubljana, SloveniaNumerical methods based on finite differences for solving the wave equation are prone to various errors, a challenge that also applies to the semi-discrete case of a recently proposed quantum algorithm that preserves continuous time while discretizing space into a quantum state vector. The inherent limitations of quantum processing, which affect the generality and extensibility of the quantum algorithm compared to classical numerical solvers, are offset by the logarithmic scaling of the number of qubits with the number of grid points. The higher grid density in combination with the higher-order approximation of the continuous Laplacian operator improves the numerical accuracy. In contrast to the Dirichlet boundary condition, however, the Neumann boundary condition in the quantum algorithm cannot exceed a second-order Laplacian approximation. Moreover, the anisotropy error is inherent to the lattice geometry. In this work, we propose a quantum wave equation solver on a hexagonal grid, which exhibits intrinsic grid dispersion comparable to that of a fourth-order accurate method on a rectangular grid and provides a quadratic improvement in the anisotropy measures. We show that the graph Laplacian on a hexagonal grid embedding the Neumann boundary condition is a bidirected graph. To make it suitable for use as a system Hamiltonian, a symmetrization and factorization based on the asymmetric adjacency matrix is proposed. Dispersion relations for quantum algorithms solving the wave equation on different lattices are derived and compared. The algorithms are analyzed as a limiting case of Yee’s finite-difference time-domain method in the context of electromagnetic wave propagation.https://ieeexplore.ieee.org/document/10679139/Electromagneticsfinite difference methodshexagonal latticesnumerical dispersionnumerical modelingradio wave propagation |
| spellingShingle | Roman Novak Quantum Algorithm for Solving the Wave Equation on a Hexagonal Grid IEEE Access Electromagnetics finite difference methods hexagonal lattices numerical dispersion numerical modeling radio wave propagation |
| title | Quantum Algorithm for Solving the Wave Equation on a Hexagonal Grid |
| title_full | Quantum Algorithm for Solving the Wave Equation on a Hexagonal Grid |
| title_fullStr | Quantum Algorithm for Solving the Wave Equation on a Hexagonal Grid |
| title_full_unstemmed | Quantum Algorithm for Solving the Wave Equation on a Hexagonal Grid |
| title_short | Quantum Algorithm for Solving the Wave Equation on a Hexagonal Grid |
| title_sort | quantum algorithm for solving the wave equation on a hexagonal grid |
| topic | Electromagnetics finite difference methods hexagonal lattices numerical dispersion numerical modeling radio wave propagation |
| url | https://ieeexplore.ieee.org/document/10679139/ |
| work_keys_str_mv | AT romannovak quantumalgorithmforsolvingthewaveequationonahexagonalgrid |