A variational quantum algorithm for tackling multi-dimensional Poisson equations with inhomogeneous boundary conditions

Minjin Choi,Hoon Ryu

2024-11-05

Abstract:We design a variational quantum algorithm to solve multi-dimensional Poisson equations with mixed boundary conditions that are typically required in various fields of computational science. Employing an objective function that is formulated with the concept of the minimal potential energy, we not only present in-depth discussion on the cost-efficient & noise-robust design of quantum circuits that are essential for evaluation of the objective function, but, more remarkably, employ the proposed algorithm to calculate bias-dependent spatial distributions of electric fields in semiconductor systems that are described with a two-dimensional domain and up to 10-qubit circuits. Extending the application scope to multi-dimensional problems with mixed boundary conditions for the first time, fairly solid computational results of this work clearly demonstrate the potential of variational quantum algorithms to tackle Poisson equations derived from physically meaningful problems.

Quantum Physics,Computational Physics

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to design a Variational Quantum Algorithm (VQA) for solving multi - dimensional Poisson equations with mixed boundary conditions. Specifically, the goals of the paper include:

1. **Solving multi - dimensional Poisson equations**:
- The Poisson equation is a second - order partial differential equation and plays an important role in many computational science fields such as fluid dynamics, nanoelectronics, and computational electromagnetics.
- Traditional numerical methods have high computational costs when dealing with large - scale problems, while quantum computing is expected to provide more efficient solutions by exploiting quantum phenomena.

2. **Handling mixed boundary conditions**:
- Multi - dimensional Poisson equations usually need to handle different types of boundary conditions, such as Dirichlet boundary conditions and Neumann boundary conditions.
- The method proposed in the paper can handle these mixed boundary conditions, expanding the application range of the variational quantum algorithm.

3. **Improving the cost - efficiency and noise - robustness of the algorithm**:
- By adopting the concept of minimum potential energy to construct the objective function, the paper designs a quantum circuit that is cost - efficient and robust to noise.
- In particular, the paper optimizes the design of the quantum circuit so that it only contains one - qubit and two - qubit gates, thus achieving better performance on NISQ (Noisy Intermediate - Scale Quantum) devices.

4. **Verification in practical applications**:
- The paper verifies the effectiveness of the algorithm by calculating the electric field distribution in semiconductor systems, especially in two - dimensional domains, using circuits with up to 10 qubits.
- The results show that the algorithm can solve physically meaningful multi - dimensional Poisson equation problems with high precision.

### Main contributions

- **Algorithm design**: Proposed a variational quantum algorithm that can effectively solve multi - dimensional Poisson equations on NISQ devices.
- **Cost and noise - robustness**: By optimizing the quantum circuit design, improved the cost - efficiency and noise - robustness of the algorithm.
- **Practical application**: Verified the effectiveness of the algorithm in practical problems by calculating the electric field distribution in semiconductor systems.

### Formula summary

1. **Poisson equation**:
\[
-\Delta u(x)=f(x), \quad x\in(0, 1)^d
\]
where \(\Delta\) is the Laplace operator and \(f(x)\) is a given smooth function.

2. **Linear system problem**:
\[
A\mathbf{v}=\mathbf{b}
\]
where \(A\) is the system matrix, \(\mathbf{b}\) is the perturbation vector given by the function \(f\), and \(\mathbf{v}\) is the approximate solution of the Poisson equation.

3. **Objective function**:
\[
E(\theta)=-\frac{1}{2}\frac{\left(\text{Re}\langle\mathbf{b}|\psi(\theta)\rangle\right)^2}{\langle\psi(\theta)|A|\psi(\theta)\rangle}
\]
where \(\mathbf{b}=\frac{\mathbf{b}}{\|\mathbf{b}\|}\), \(\psi(\theta)\) is the parameterized trial state, and \(\text{Re}\langle\mathbf{b}|\psi(\theta)\rangle\) is the real part of \(\langle\mathbf{b}|\psi(\theta)\rangle\).

4. **Partial derivative**:
\[
\frac{\partial E(\theta)}{\partial}