Entropy is a measure of the randomness of a system. Estimating the entropy of a quantum state is a basic problem in quantum information. In this paper, we introduce a time-efficient quantum approach to estimating the von Neumann entropy
S(ρ) and Rényi entropy
Sα(ρ) of an
N-dimensional quantum state
ρ, given access to independent samples of
ρ. Specifically, we provide the following: 1. A quantum estimator for
S(ρ) with time complexity
O~(N2), improving the prior best time complexity
O~(N6) by Acharya, Issa, Shende, and Wagner (2020) and Bavarian, Mehraba, and Wright (2016). 2. A quantum estimator for
Sα(ρ) with time complexity
O~(N4/α−2) for
0<α<1 and
O~(N4−2/α) for
α>1, improving the prior best time complexity
O~(N6/α) for
0<α<1 and
O~(N6) for
α>1 by Acharya, Issa, Shende, and Wagner (2020), though at a cost of a slightly larger sample complexity. Moreover, these estimators are naturally extensible to the low-rank case. We also provide a sample lower bound for estimating
Sα(ρ). Technically, our method is quite different from the previous ones that are based on weak Schur sampling and Young diagrams. At the heart of our construction, is a novel tool called samplizer, which can "samplize" a quantum query algorithm to a quantum algorithm with similar behavior using only samples of quantum states; this suggests a unified framework for estimating quantum entropies. Specifically, when a quantum oracle
U block-encodes a mixed quantum state
ρ, any quantum query algorithm using
Q queries to
U can be samplized to a
δ-close (in the diamond norm) quantum algorithm using
Θ~(Q2/δ) samples of
ρ. Moreover, this samplization is proven to be optimal, up to a polylogarithmic factor.
The von Neumann and quantum R\'enyi entropies characterize fundamental properties of quantum systems and lead to many theoretical and practical applications. Quantum algorithms using a purified quantum query model can speed up quantum entropy estimation, while little is known about the complexity of using identical copies of the quantum state. This paper presents quantum entropy estimation algorithms with a cost of copies scaling polynomially in the rank of the state. In contrast to current methods that depend on the dimension of the system, our methods could provide exponential resource savings in the scenario of low-rank states. Furthermore, we show how to construct quantum circuits using primitive single-qubit or two-qubit gates efficiently and thus provide practical methods for estimating quantum entropies of quantum systems. We also conduct simulation experiments to show the effectiveness and noise robustness of our algorithms.
We describe a quantum algorithm to estimate the $\ensuremath{\alpha}$-Renyi entropy of an unknown density matrix $\ensuremath{\rho}\ensuremath{\in}{\mathbb{C}}^{d\ifmmode\times\else\texttimes\fi{}d}$ for $\ensuremath{\alpha}\ensuremath{\ne}1$ by combining the recent technique of quantum singular value transformations with the method of estimating normalized traces in the one clean qubit model. We consider an oracular input model where the input state is prepared via a quantum oracle that outputs a purified version of the state, assumed to be nonsingular. Our method outputs an estimate of the $\ensuremath{\alpha}$-Renyi entropy to additive precision $\ensuremath{\epsilon}$, using an expected total number $O(\frac{1}{{(x\ensuremath{\epsilon})}^{2}})$ of independent applications of a quantum circuit which coherently queries the input unitary $O(\frac{1}{\ensuremath{\delta}}log\frac{d}{\ensuremath{\epsilon}})$ times, in each case measuring a single output qubit. Here $\ensuremath{\delta}$ is a lower cutoff on the smallest eigenvalue of $\ensuremath{\rho}$ and $x=\frac{1}{d}\mathrm{Tr}({\ensuremath{\rho}}^{\ensuremath{\alpha}})$. The expected number of measurements made in this method can be compared to results in the sample complexity model that generally require $\mathrm{\ensuremath{\Theta}}(\frac{{d}^{2}}{{\ensuremath{\epsilon}}^{2}})$ samples. Furthermore, we also show that multiplicative approximations can be obtained by iteratively using additive approximations, with an overhead logarithmic in the dimension
d.
We investigate the computational hardness of estimating the quantum
α-Rényi entropy
SαR(ρ)=1−αlnTr(ρα) and the quantum
q-Tsallis entropy
SqT(ρ)=q−11−Tr(ρq), both of which converge to the von Neumann entropy as the order approaches
1. The promise problems Quantum
α-Rényi Entropy Approximation (RényiQEA
α) and Quantum
q-Tsallis Entropy Approximation (TsallisQEA
q) ask whether $ {\rm S}^ {\tt R}_α(ρ)$ or
SqT(ρ), is at least
τ1 or at most
τ2, where
τ1−τ2 is typically a positive constant. Previous hardness results cover only the von Neumann entropy (order
1) and some cases of the quantum
q-Tsallis entropy, while existing approaches do not readily extend to other orders.
We establish that for all positive real
α and
q, and also for
α=∞, the rank-
2 variants Rank2RényiQEA
α and Rank2TsallisQEA
q are BQP-hard. Combined with prior (rank-dependent) quantum query algorithms in Wang, Guan, Liu, Zhang, and Ying (TIT 2024), Wang, Zhang, and Li (TIT 2024), and Liu and Wang (SODA 2025), as well as the one derived from O'Donnell and Wright (STOC 2016), our results imply:
- For all real orders
α>0 or
α=∞, and for all real orders
0<q≤1, LowRankRényiQEA
α and LowRankTsallisQEA
q are BQP-complete, where both are restricted versions of RényiQEA
α and TsallisQEA
q with
ρ of polynomial rank.
- For all real order
q>1, TsallisQEA
q is BQP-complete.
Our hardness results stem from reductions based on new inequalities relating the
α-Rényi or
q-Tsallis binary entropies of different orders. These reductions differ substantially from previous approaches, and the inequalities are of independent interest.
The entropy of a quantum system is a measure of its randomness and is useful in quantifying entanglement. We study the problem of measuring the von Neumann and Rényi entropies of an unknown mixed quantum state given access to independent copies of the state. For Rényi entropy of integral order exceeding one, we determine the order-optimal copy complexity and show that it is strictly lower than the number of copies required to learn the underlying state. The main technical innovation is a concentration result for certain polynomials that arise in the Kerov algebra of Young diagrams, which is proven using the cycle structure of compositions of certain types of permutations. For von Neumann entropy and Rényi entropy of non-integral orders, we provide upper and lower bounds on the sample complexity of the Empirical Young Diagram (EYD) algorithm, which is the analogue of the empirical plug-in estimator in classical estimation.
Entropy is a measure of the randomness of a system. Estimating the entropy of a quantum state is a basic problem in quantum information. In this paper, we introduce a time-efficient quantum approach to estimating the von Neumann entropyS(ρ) and Rényi entropySα(ρ) of anN-dimensional quantum state ρ, given access to independent samples of ρ. Specifically, we provide the following quantum estimators. • A quantum estimator forS(ρ) with time complexity Õ (N2),1improving the prior best time complexity Õ(N6) by Acharya, Issa, Shende, and Wagner (2020) and Bavarian, Mehraba, and Wright (2016). • A quantum estimator forSα(ρ) with time complexity Õ (N4/α−2) for 0N4−2/α) for α > 1, improving the prior best time complexity Õ(N6/α) for 0N6) for α > 1 by Acharya, Issa, Shende, and Wagner (2020), though at a cost of a slightly larger sample complexity. Moreover, these estimators are naturally extensible to the low-rank case. We also provide a sample lower bound Ω(max{N/ε,N1/α−1/ε1/α}) for estimatingSα(ρ). Technically, our method is quite different from the previous ones that are based on weak Schur sampling and Young diagrams. At the heart of our construction, is a novel tool calledsamplizer, which can “samplize” a quantum query algorithm to a quantum algorithm with similar behavior using only samples of quantum states; this suggests a new framework for estimating quantum entropies. Specifically, when a quantum oracleUblockencodes a mixed quantum state ρ, any quantum query algorithm usingQqueries toUcan be samplized to a δ-close (in the diamond norm) quantum algorithm using Θ(Q2/δ) samples of ρ. Moreover, this samplization is proven to be optimal, up to a polylogarithmic factor.