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:2601.06492 · IEEE International Symposium on Information Theory (ISIT)

Algorithms for Computing the Petz-Augustin Capacity

Chun-Neng Chu, Wei-Fu Tseng, Yen-Huan Li

cs.ITmath.OCquant-phSubmitted Jan 10, 2026• Announced Jan 13, 2026• 1 citation
arXivHTMLPDF

Abstract

We propose the first algorithms with non-asymptotic convergence guarantees for computing the Petz-Augustin capacity, which generalizes the channel capacity and characterizes the optimal error exponent in classical-quantum channel coding. This capacity can be equivalently expressed as the maximization of two generalizations of mutual information: the Petz-Rényi information and the Petz-Augustin information. To maximize the Petz-Rényi information, we show that it corresponds to a convex Hölder-smooth optimization problem, and hence the universal fast gradient method of Nesterov (2015), along with its convergence guarantees, readily applies. Regarding the maximization of the Petz-Augustin information, we adopt a two-layered approach: we show that the objective function is smooth relative to the negative Shannon entropy and can be efficiently optimized by entropic mirror descent; each iteration of entropic mirror descent requires computing the Petz-Augustin information, for which we propose a novel fixed-point algorithm and establish its contractivity with respect to the Thompson metric. Notably, this two-layered approach can be viewed as a generalization of the mirror-descent interpretation of the Blahut-Arimoto algorithm due to He et al. (2024).

Related papers

A Linearly Convergent Algorithm for Petz-Augustin Means with Applications to Computing Petz-Augustin Capacity and Equilibrium Prices in Fisher Markets

Chun-Neng Chu, Tseng, Wei-Fu, Yen-Huan Li

We study the computation of the Petz-Augustin mean of order $α\in (0,1) \cup (1,\infty)$, defined as the minimizer of a weighted sum of $n$ Petz-Rényi divergences of order $α$ over the set of $d$-by-$d$ quantum states, where the Petz-Rényi divergence is a quantum generalization of the classical Rényi divergence. We propose the first algorithm with a non-asymptotic convergence guarantee for solving this optimization problem. The iterates are guaranteed to converge to the Petz-Augustin mean at a linear rate of $O\left( \lvert 1 - 1/α\rvert^T \right)$ with respect to the Thompson metric for $α\in(1/2,1)\cup(1,\infty)$, where $T$ denotes the number of iterations. The algorithm has an initialization time complexity of $O\left(nd^3\right)$ and a per-iteration time complexity of $O\left(nd^2 + d^3\right)$. Two applications follow. The first application concerns the Petz-Augustin capacity of order $α\in(1/2,1)$, which generalizes the quantum channel capacity and characterizes the optimal error exponent in classical-quantum channel coding. To compute this quantity, we propose the first iterative method with a non-asymptotic convergence guarantee, which uses our algorithm for computing the Petz-Augustin mean as a subroutine. The second application concerns the equilibrium prices in Fisher markets with constant elasticity of substitution (CES) utilities. We establish that when all buyers have CES utilities of common elasticity, the equilibrium prices are equivalent to the Petz-Augustin mean when all quantum states commute. In this case, our proposed algorithm is a tâtonnement process à la Walras. We further extend the algorithm to inhomogeneous Fisher markets, where buyers have CES utilities with different elasticities, and prove a faster convergence rate than existing tâtonnement-type algorithms.

arXiv:2502.06399•Feb 10, 2025•0 citations

A Mirror-Descent Algorithm for Computing the Petz-Rényi Capacity of Classical-Quantum Channels

Yu-Hong Lai, Hao-Chung Cheng

We study the computation of the $α$-Rényi capacity of a classical-quantum (c-q) channel for $α\in(0,1)$. We propose an exponentiated-gradient (mirror descent) iteration that generalizes the Blahut-Arimoto algorithm. Our analysis establishes relative smoothness with respect to the entropy geometry, guaranteeing a global sublinear convergence of the objective values. Furthermore, under a natural tangent-space nondegeneracy condition (and a mild spectral lower bound in one regime), we prove local linear (geometric) convergence in Kullback-Leibler divergence on a truncated probability simplex, with an explicit contraction factor once the local curvature constants are bounded.

quant-pharXiv:2601.10558•IEEE International Symposium on Information Theory (ISIT)•Jan 16, 2026•1 citation

Computing Quantum Channel Capacities

Navneeth Ramakrishnan, Raban Iten, Volkher B. Scholz +1 more

The capacity of noisy quantum channels characterizes the highest rate at which information can be reliably transmitted and it is therefore of practical as well as fundamental importance. Capacities of classical channels are computed using alternating optimization schemes, called Blahut-Arimoto algorithms. In this work, we generalize classical Blahut-Arimoto algorithms to the quantum setting. In particular, we give efficient iterative schemes to compute the capacity of channels with classical input and quantum output, the quantum capacity of less noisy channels, the thermodynamic capacity of quantum channels, as well as the entanglement-assisted capacity of quantum channels. We give rigorousa priorianda posterioribounds on the estimation error by employing quantum entropy inequalities and demonstrate fast convergence of our algorithms in numerical experiments.

arXiv:1905.01286•IEEE Transactions on Information Theory•Oct 28, 2020•33 citations

Minimizing Quantum Renyi Divergences via Mirror Descent with Polyak Step Size

Jun-Kai You, Hao–Chung Cheng, Yen-Huan Li

Quantum information quantities play a substantial role in characterizing operational quantities in various quantum information-theoretic problems. We consider numerical computation of four quantum information quantities: Petz-Augustin information, sandwiched Augustin information, conditional sandwiched Renyi entropy and sandwiched Renyi information. To compute these quantities requires minimizing some order-$α$ quantum Renyi divergences over the set of quantum states. Whereas the optimization problems are obviously convex, they violate standard bounded gradient/Hessian conditions in literature, so existing convex optimization methods and their convergence guarantees do not directly apply. In this paper, we propose a new class of convex optimization methods called mirror descent with the Polyak step size. We prove their convergence under a weak condition, showing that they provably converge for minimizing quantum Renyi divergences. Numerical experiment results show that entropic mirror descent with the Polyak step size converges fast in minimizing quantum Renyi divergences.

arXiv:2109.06054•Sep 13, 2021•1 citation

Minimizing Quantum Rényi Divergences via Mirror Descent with Polyak Step Size

Jun-Kai You, Hao–Chung Cheng, Yen-Huan Li

Quantum information quantities play a substantial role in characterizing operational quantities in various quantum information-theoretic problems. We consider the numerical computation of four quantum information quantities: Petz-Augustin information, sandwiched Augustin information, conditional sandwiched Rényi entropy, and sandwiched Rényi information. To compute these quantities requires minimizing some order-α quantum Rényi divergences over the set of quantum states. Whereas the optimization problems are obviously convex, they violate standard bounded gradient/Hessian conditions in the literature, so existing convex optimization methods and their convergence guarantees do not directly apply. In this paper, we propose a new class of convex optimization methods called mirror descent with the Polyak step size. We prove their convergence under a weak condition, showing that they provably converge for minimizing quantum Rényi divergences. Numerical experiment results show that entropic mirror descent with the Polyak step size converges fast in minimizing quantum Rényi divergences.

IEEE International Symposium on Information Theory (ISIT)•Jun 26, 2022•15 citations

A Bregman Proximal Perspective on Classical and Quantum Blahut-Arimoto Algorithms

Kerry He, James Saunderson, Hamza Fawzi

The Blahut-Arimoto algorithm is a well-known method to compute classical channel capacities and rate-distortion functions. Recent works have extended this algorithm to compute various quantum analogs of these quantities. In this paper, we show how these Blahut-Arimoto algorithms are special instances of mirror descent, which is a type of Bregman proximal method, and a well-studied generalization of gradient descent for constrained convex optimization. Using recently developed convex analysis tools, we show how analysis based on relative smoothness and strong convexity recovers known sublinear and linear convergence rates for Blahut-Arimoto algorithms. This Bregman proximal viewpoint allows us to derive related algorithms with similar convergence guarantees to solve problems in information theory for which Blahut-Arimoto-type algorithms are not directly applicable. We apply this framework to compute energy-constrained classical and quantum channel capacities, classical and quantum rate-distortion functions, and approximations of the relative entropy of entanglement, all with provable convergence guarantees.

arXiv:2306.04492•IEEE Transactions on Information Theory•Jun 11, 2024•7 citations
We propose the first algorithms with non-asymptotic convergence guarantees for computing the Petz-Augustin capacity, which generalizes the channel capacity and characterizes the optimal error exponent in classical-quantum channel coding. This capacity can be equivalently expressed as the maximization of two generalizations of mutual information: the Petz-Rényi information and the Petz-Augustin information. To maximize the Petz-Rényi information, we show that it corresponds to a convex Hölder-smooth optimization problem, and hence the universal fast gradient method of Nesterov (2015), along with its convergence guarantees, readily applies. Regarding the maximization of the Petz-Augustin information, we adopt a two-layered approach: we show that the objective function is smooth relative to the negative Shannon entropy and can be efficiently optimized by entropic mirror descent; each iteration of entropic mirror descent requires computing the Petz-Augustin information, for which we propose a novel fixed-point algorithm and establish its contractivity with respect to the Thompson metric. Notably, this two-layered approach can be viewed as a generalization of the mirror-descent interpretation of the Blahut-Arimoto algorithm due to He et al. (2024).
We study the computation of the Petz-Augustin mean of order α∈(0,1)∪(1,∞)α\in (0,1) \cup (1,\infty)α∈(0,1)∪(1,∞), defined as the minimizer of a weighted sum of nnn Petz-Rényi divergences of order ααα over the set of ddd-by-ddd quantum states, where the Petz-Rényi divergence is a quantum generalization of the classical Rényi divergence. We propose the first algorithm with a non-asymptotic convergence guarantee for solving this optimization problem. The iterates are guaranteed to converge to the Petz-Augustin mean at a linear rate of O(∣1−1/α∣T)O\left( \lvert 1 - 1/α\rvert^T \right)O(∣1−1/α∣T) with respect to the Thompson metric for α∈(1/2,1)∪(1,∞)α\in(1/2,1)\cup(1,\infty)α∈(1/2,1)∪(1,∞), where TTT denotes the number of iterations. The algorithm has an initialization time complexity of O(nd3)O\left(nd^3\right)O(nd3) and a per-iteration time complexity of O(nd2+d3)O\left(nd^2 + d^3\right)O(nd2+d3). Two applications follow. The first application concerns the Petz-Augustin capacity of order α∈(1/2,1)α\in(1/2,1)α∈(1/2,1), which generalizes the quantum channel capacity and characterizes the optimal error exponent in classical-quantum channel coding. To compute this quantity, we propose the first iterative method with a non-asymptotic convergence guarantee, which uses our algorithm for computing the Petz-Augustin mean as a subroutine. The second application concerns the equilibrium prices in Fisher markets with constant elasticity of substitution (CES) utilities. We establish that when all buyers have CES utilities of common elasticity, the equilibrium prices are equivalent to the Petz-Augustin mean when all quantum states commute. In this case, our proposed algorithm is a tâtonnement process à la Walras. We further extend the algorithm to inhomogeneous Fisher markets, where buyers have CES utilities with different elasticities, and prove a faster convergence rate than existing tâtonnement-type algorithms.
We study the computation of the ααα-Rényi capacity of a classical-quantum (c-q) channel for α∈(0,1)α\in(0,1)α∈(0,1). We propose an exponentiated-gradient (mirror descent) iteration that generalizes the Blahut-Arimoto algorithm. Our analysis establishes relative smoothness with respect to the entropy geometry, guaranteeing a global sublinear convergence of the objective values. Furthermore, under a natural tangent-space nondegeneracy condition (and a mild spectral lower bound in one regime), we prove local linear (geometric) convergence in Kullback-Leibler divergence on a truncated probability simplex, with an explicit contraction factor once the local curvature constants are bounded.
The capacity of noisy quantum channels characterizes the highest rate at which information can be reliably transmitted and it is therefore of practical as well as fundamental importance. Capacities of classical channels are computed using alternating optimization schemes, called Blahut-Arimoto algorithms. In this work, we generalize classical Blahut-Arimoto algorithms to the quantum setting. In particular, we give efficient iterative schemes to compute the capacity of channels with classical input and quantum output, the quantum capacity of less noisy channels, the thermodynamic capacity of quantum channels, as well as the entanglement-assisted capacity of quantum channels. We give rigorousa priorianda posterioribounds on the estimation error by employing quantum entropy inequalities and demonstrate fast convergence of our algorithms in numerical experiments.
Quantum information quantities play a substantial role in characterizing operational quantities in various quantum information-theoretic problems. We consider numerical computation of four quantum information quantities: Petz-Augustin information, sandwiched Augustin information, conditional sandwiched Renyi entropy and sandwiched Renyi information. To compute these quantities requires minimizing some order-ααα quantum Renyi divergences over the set of quantum states. Whereas the optimization problems are obviously convex, they violate standard bounded gradient/Hessian conditions in literature, so existing convex optimization methods and their convergence guarantees do not directly apply. In this paper, we propose a new class of convex optimization methods called mirror descent with the Polyak step size. We prove their convergence under a weak condition, showing that they provably converge for minimizing quantum Renyi divergences. Numerical experiment results show that entropic mirror descent with the Polyak step size converges fast in minimizing quantum Renyi divergences.
Quantum information quantities play a substantial role in characterizing operational quantities in various quantum information-theoretic problems. We consider the numerical computation of four quantum information quantities: Petz-Augustin information, sandwiched Augustin information, conditional sandwiched Rényi entropy, and sandwiched Rényi information. To compute these quantities requires minimizing some order-α quantum Rényi divergences over the set of quantum states. Whereas the optimization problems are obviously convex, they violate standard bounded gradient/Hessian conditions in the literature, so existing convex optimization methods and their convergence guarantees do not directly apply. In this paper, we propose a new class of convex optimization methods called mirror descent with the Polyak step size. We prove their convergence under a weak condition, showing that they provably converge for minimizing quantum Rényi divergences. Numerical experiment results show that entropic mirror descent with the Polyak step size converges fast in minimizing quantum Rényi divergences.
The Blahut-Arimoto algorithm is a well-known method to compute classical channel capacities and rate-distortion functions. Recent works have extended this algorithm to compute various quantum analogs of these quantities. In this paper, we show how these Blahut-Arimoto algorithms are special instances of mirror descent, which is a type of Bregman proximal method, and a well-studied generalization of gradient descent for constrained convex optimization. Using recently developed convex analysis tools, we show how analysis based on relative smoothness and strong convexity recovers known sublinear and linear convergence rates for Blahut-Arimoto algorithms. This Bregman proximal viewpoint allows us to derive related algorithms with similar convergence guarantees to solve problems in information theory for which Blahut-Arimoto-type algorithms are not directly applicable. We apply this framework to compute energy-constrained classical and quantum channel capacities, classical and quantum rate-distortion functions, and approximations of the relative entropy of entanglement, all with provable convergence guarantees.