Abstract:We introduce QICS (Quantum Information Conic Solver), an open-source primal-dual interior point solver fully implemented in Python, which is focused on solving optimization problems arising in quantum information theory. QICS has the ability to solve optimization problems involving the quantum relative entropy, noncommutative perspectives of operator convex functions, and related functions. It also includes an efficient semidefinite programming solver which exploits sparsity, as well as support for Hermitian matrices. QICS is also currently supported by the Python optimization modelling software PICOS. This paper aims to document the implementation details of the algorithm and cone oracles used in QICS, and serve as a reference guide for the software. Additionally, we showcase extensive numerical experiments which demonstrate that QICS outperforms state-of-the-art quantum relative entropy programming solvers, and has comparable performance to state-of-the-art semidefinite programming solvers.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop an efficient, open - source optimization solver specifically for solving optimization problems in quantum information theory. Specifically, the paper introduces QICS (Quantum Information Conic Solver), a primal - dual interior - point solver fully implemented in Python, which focuses on solving optimization problems involving quantum relative entropy, non - commutative operator convex function perspectives and their related functions.
### Main Problems and Goals
1. **Solve optimization problems in quantum information theory**:
- QICS can handle optimization problems involving quantum relative entropy, non - commutative operator convex function perspectives and their related functions.
- Quantum relative entropy is defined as:
\[
S(X \| Y) := \text{tr}[X(\log(X) - \log(Y))]
\]
where \( X, Y \in H_n^+ \) and satisfy \(\ker(Y) \subseteq \ker(X)\).
2. **Improve solution efficiency**:
- QICS solves asymmetric cone programming problems by using the Skajaa - Ye algorithm and the stepping algorithm in Hypatia.
- The self - consistency of the natural barrier functions of the quantum relative entropy cone and other related cones has been proven, enabling QICS to use the interior - point method to solve these problems.
3. **Support sparsity and Hermitian matrices**:
- QICS uses sparsity techniques to more efficiently construct Schur complement matrices, thus accelerating the solution process.
- It supports complex Hermitian matrix variables, not just real symmetric matrices, which is especially important in quantum information theory.
4. **Comparison with other solvers**:
- The paper presents extensive numerical experiments, showing that QICS outperforms existing state - of - the - art solvers in solving quantum relative entropy programming problems, and its performance in solving semidefinite programming problems is comparable to that of state - of - the - art solvers.
### Implementation Methods
- **Primal - dual interior - point algorithm**: Use the Nesterov - Todd algorithm to solve symmetric cone programming problems, including linear, semidefinite and second - order cone programming.
- **Support for quantum relative entropy cones**: Implement effective algorithms for quantum relative entropy cones and their related cones.
- **Efficient linear system solution**: Simplify the solution steps of the Hessian matrix through pre - image techniques and effective algorithms for specific slices.
### Application Areas
QICS is mainly applied to optimization problems in quantum information theory, such as security certification of quantum key distribution protocols, quantum conditional entropy calculation, etc. In addition, it also supports classical entropy cones and non - commutative perspective cones, making it potentially useful in a wider range of application scenarios.
In summary, this paper aims to provide an efficient, flexible and easy - to - use tool to solve complex optimization problems in quantum information theory.