Abstract:Quantum linear system algorithms (QLSAs) for gate-based quantum computing can provide exponential speedups for solving linear systems but face challenges when applied to finite element problems due to the growth of the condition number with problem size. Furthermore, QLSAs cannot use an approximate solution or initial guess to output an improved solution. Here, we present Quantum Relaxation for Linear System (qRLS), as an iterative approach for gate-based quantum computers by embedding linear stationary iterations into a larger block linear system. The condition number of the block linear system scales linearly with the number of iterations independent of the size and condition number of the original system. The well-conditioned system enables a practical iterative solution of finite element problems using the state-of-the-art Quantum Signal Processing (QSP) variant of QLSAs, for which we provide numerical results using a quantum computer simulator. The iteration complexity demonstrates favorable scaling relative to classical architectures, as the solution time is independent of system size and requires O(log(N)) qubits. This represents an exponential efficiency gain, offering a new approach for iterative finite element problem-solving on quantum hardware.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is: when using the finite element method (FEM) to solve linear systems, how to use gate - based quantum computing to overcome the problem in classical algorithms where the condition number grows with the problem size. Specifically, the paper proposes a new method called Quantum Relaxation for Linear System (qRLS). By embedding linear stable iterations into a larger block - linear system, the condition number is only related to the number of iterations, and has nothing to do with the size and condition number of the original system. This enables qRLS to provide a practical iterative solution while maintaining exponential - level efficiency, and is especially suitable for solving finite - element problems.
### Key Point Summary:
1. **Background Problems**:
- Classical algorithms (such as the conjugate gradient method) usually have a complexity that is linearly related to the number of unknowns \( N \) when solving large - scale linear systems.
- Quantum Linear System Algorithms (QLSAs) can theoretically provide exponential - level acceleration, but in practical applications, they face the challenge that the condition number grows with the problem size.
2. **Proposed Solutions**:
- **qRLS Method**: By embedding linear stable iterations into a larger block - linear system, the condition number of the new system is only related to the number of iterations, thus avoiding the influence of the condition number of the original system.
- **Condition Number Control**: The condition number of the new system has a linear relationship with the number of iterations, rather than being proportional to the size or condition number of the original system.
- **Quantum Signal Processing (QSP)**: Use QLSA variants of QSP to solve the new system, achieving efficient and accurate solutions.
3. **Advantages**:
- **Exponential - level Efficiency**: qRLS is independent of the system size in terms of solution time, and only requires \( O(\log(N)) \) qubits, achieving an exponential - level efficiency gain.
- **Scope of Application**: Applicable to positive - definite linear systems, especially for solving finite - element problems.
4. **Technical Details**:
- **Block Encoding**: Encode the iteration process into a larger block - diagonal linear system, which is implemented by a quantum circuit.
- **State Preparation**: Use a quantum oracle to efficiently prepare the initial state.
- **Measurement and Precision**: Extract the final iteration result by measuring the first register of the quantum state, and analyze the success probability and precision requirements of the measurement.
### Formula Summary:
- **Block Matrix**:
\[
M_{l,c} = \sum_{i = 0}^{l + c} |i\rangle\langle i| \otimes I - \sum_{i = 1}^{l} |i\rangle\langle i - 1| \otimes R - \sum_{i = l + 1}^{l + c} |i\rangle\langle i - 1| \otimes I
\]
where \( R = I - \tau A \), and \( \tau \) is selected to ensure \( \| \tau A \| \leq 1 \).
- **Condition Number**:
\[
\kappa(M) \leq 2l
\]
- **Success Probability**:
- For the case of zero initial guess:
\[
p \geq \frac{1}{2}
\]
- For the case of non - zero initial guess:
\[
p \geq \frac{1}{2} \left( \frac{\| x^* \| - \epsilon}{\| x^* \| + \epsilon} \right)^2
\]
- **Precision Requirements**:
- For the case of zero initial guess:
\[
\delta \leq \frac{\epsilon}{2\sqrt{2l}}
\]
- For non - zero initial guess: