Quantum State Tomography Via Compressed Sensing
Quantum State Tomography Via Compressed Sensing
net/publication/49751450
CITATIONS READS
665 205
5 authors, including:
Jens Eisert
Freie Universität Berlin
366 PUBLICATIONS 20,437 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Jens Eisert on 31 May 2014.
The tasks of reconstructing the quantum states and minimum-rank matrix subject to linear constraints is NP-
processes produced by physical systems—known respec- hard in general [8].
tively as quantum state and process tomography [1]—are In addition to a reduction in experimental complexity,
of increasing importance in physics and especially in one might hope that a postprocessing algorithm taking as
quantum information science. Tomography has been used input only OðrdÞ d2 numbers could be tuned to run
to characterize the quantum state of trapped ions [2] and an considerably faster than standard methods. Since the out-
optical entangling gate [3] among many other implemen- put of the procedure is a low-rank approximation to the
tations. But a fundamental difficulty in performing tomo- density operator and only requires OðrdÞ numbers be
graphy on many-body systems is the exponential growth in specified, it becomes conceivable that the run time scales
the state space dimension. For example, to get a maximum- better than Oðd2 Þ, clearly impossible for naive approaches
likelihood estimate of a quantum state of 8 ions, Ref. [2] using dense matrices.
required hundreds of thousands of measurements and In this Letter, we introduce a method to achieve such
weeks of postprocessing. drastic reductions in measurement complexity, together
Still, one might hope to overcome this obstacle because with efficient algorithms for postprocessing. The approach
the vast majority of quantum states are not of physical further develops ideas that have recently been studied under
interest. Rather, one is often interested in states with special the label of ‘‘compressed sensing.’’ Compressed sensing [9]
properties: pure states, states with particular symmetries, provides techniques for recovering a sparse vector from a
ground states of local Hamiltonians, etc., and tomography small number of measurements [10]. Here, sparsity means
might be more efficient in such special cases [4]. that this vector contains only a few nonzero entries in a
In particular, consider pure or nearly pure quantum specified basis, and the measurements are linear functions
states, i.e., states with low entropy. More precisely, con- of its entries. When the measurements are chosen at random
sider a quantum state that is essentially supported on an (in a certain precise sense), then with high probability two
r-dimensional space, meaning the density matrix is close surprising things happen: the vector is uniquely determined
(in a given norm) to a matrix of rank r, where r is small. by a small number of measurements, and it can be recovered
Such states arise in very common physical settings, e.g., a by an efficient convex optimization algorithm [9].
pure state subject to a local noise process [5]. Matrix completion [11–13] is a generalization of com-
A standard implementation of tomography [6,7] would pressed sensing from vectors to matrices. Here, one recov-
use d2 or more measurement settings, where d ¼ 2n for an ers certain ‘‘incoherent’’ low-rank matrices X from a small
n-qubit system. But a simple parameter counting argument number of matrix elements Xi;j . The problem of low-rank
suggests that OðrdÞ settings could possibly suffice—a quantum state tomography bears a strong resemblance to
significant improvement. However, it is not clear how to matrix completion. However, there are important differ-
achieve this performance in practice, i.e., how to choose ences. We wish to use measurements that can be more
these measurements, or how to efficiently reconstruct the easily implemented in an experiment than obtaining ele-
density matrix. For instance, the problem of finding a ments i;j of density matrices. Previous results [11–13]
cannot be applied to this more general situation. We would the solution is unique if for all deviations :¼
also like to avoid any unnatural incoherence assumptions away from either R Þ 0 or k þ ktr > kktr .
crucial in prior work [11]. We will ascertain this by using a basic idea from convex
Our first result is a protocol for tomography that over- optimization: constructing a strict subgradient Y for the
comes both of these difficulties: it uses Pauli measurements norm. A matrix Y is a strict subgradient if k þ ktr >
only, and it works for arbitrary density matrices. We prove kktr þ trY for all Þ 0. The main contribution below
that only Oðrdlog2 dÞ measurement settings suffice. What is a method for constructing such a Y which is also in the
is more, our proof introduces some new techniques, which range of R. For then R ¼ 0 implies that is orthogonal
both generalize and vastly simplify the previous work on to the range of R, thus trY ¼ 0 and the subgradient
matrix completion. condition reads k þ ktr > kktr . This implies unique-
In a real experiment, the measurements are noisy, and ness. (In fact, it suffices to approximate the condition.)
the true state is only approximately low rank. We show that Let E be the projection onto the range of , let T be the
our method is robust to these sources of error. We also space spanned by those operators whose row or column
describe ways to certify that a state is nearly pure without space is contained in range . Let P T be the projection
any a priori assumptions. onto T, P ?T onto the orthogonal complement. Decompose
Finally, we present fast algorithms for reconstructing the ¼ T þ ? T , the parts of that lie in the subspaces T
density matrix from the measurement statistics based on and T ? . We distinguish two cases: (i) kT k2 > d2 k? T k2 ,
semidefinite programming—a feature not present in earlier and (ii) kT k2 d2 k? k
T 2 [16].
methods for pure-state tomography [4,6,7]. Reconstructing Case (i) is easier. In this case, is well approximated by
a low-rank density matrix for 8 qubits takes about 1 min on T and essentially we only have to show that the restriction
an ordinary laptop computer. A :¼ P T RP T of R to T is invertible. Using a noncom-
While our methods do not overcome the exponential mutative large-deviation bound (see Refs. [18,19]),
growth in measurement complexity (which is provably 2
Pr½kA 1T jj > t < 4dret =8 (2)
impossible for any protocol capable of handling generic
pure states), they do significantly push the boundary of where ¼ m=ðdrÞ [16]. Hence the probability that kA
what can be done in a realistic setting [14]. 1T k > 12 is smaller than 4dre=32 ¼: p1 . If that is not the
Matrix recovery using Pauli measurements.—We con- case, one easily sees that kRk2 > 0, concluding the
sider the case of n spin-1=2 systems in an unknown Nstate proof for this case.
[15]. An n-qubit Pauli matrix is of the form w ¼ ni¼1 wi , Case (ii) is more involved. A matrix Y 2
where wi 2 f1; x ; y ; z g. There are d2 such matrices, spanðwðA1 Þ; . . . ; wðAm ÞÞ is an almost subgradient [20] if
labeled wðaÞ, a 2 ½1; d2 . The protocol proceeds as fol-
lows: choose m integers A1 ; . . . ; Am 2 ½1; d2 at random kP T Y Ek2 1=ð2d2 Þ; kP ?
T Yk < 1=2: (3)
and measure the expectation values trwðAi Þ. One then First, suppose such a Y exists. Then a simple calculation
solves a convex optimization problem: minimize kktr (see Ref. [18]) using the condition (ii) shows that R ¼ 0
[16] subject to indeed implies k þ ktr > kktr as hinted at above. This
tr ¼ 1; trwðAi Þ ¼ trwðAi Þ: (1) proves uniqueness in case (ii). The difficult part consists in
showing that an almost subgradient exists.
Theorem 1 (low-rank tomography)—Let be an arbi-
To this end, we design a recursive process (the ‘‘golfing
trary state of rank r. If m ¼ cdrlog2 d randomly chosen
scheme’’ [17]) which converges to a subgradient exponen-
Pauli expectations are known, then can be uniquely
tially fast. Assume we draw l batches of 0 rd Pauli
reconstructed by solving the convex optimization problem
observables independently at random (0 will be chosen
(1) with probability of failure exponentially small in c.
later). Define recursively X0 ¼ E,
The proof is inspired by, but technically very different
from, earlier work on matrix completion [11]. Our methods Xi
Yi ¼ Rj Xj1 ; Xi ¼ E P T Yi ; (4)
are more general, can be tuned to give tighter bounds, and j¼1
are much more compact, allowing us to present a fairly
complete argument in this Letter. A more detailed presen- Y ¼ Yl . Let Ri be the sampling operator associated with
tation will be published elsewhere [17]. the ith batch, and Ai its restriction to T. Assume that in
Proof.—Here we sketch the argument and explain the each run kAi 1T k2 < 1=2. Denote the probability of
main ideas; detailed calculations are in the supplementary this event not occurring by p2 . Then
material [18]. kXi k2 ¼ kXi1 P T Ri Xi1 k2 ¼ kð1T Ai ÞXi1 k2
Note that the linear constraints (1) depend only on the
projection of onto the span of the measured observables 1=2kXi1 k2 ;
pffiffiffi
wðA1 Þ; . . . ; wðAm Þ. This is preciselyPthe range of the so that kXi k2 2i kX0 k ¼ 2i r. Hence,pYffiffiffi ¼ Yl fulfills
‘‘sampling operator’’ R: ° md m i¼1 wðAi ÞtrwðAi Þ. the first part of (3), as soon as l log2 ð2d2 rÞ. We turn to
(Note that E½RðÞ ¼ .) Indeed, the convex program the second part. Again using large-deviation
pffiffiffi techniques
can be written as min kktr s.t. R ¼ R. Evidently, [18] we find kP ? T Ri Xi1 k 1=ð4 rÞkXi1 k2 with some
150401-2
week ending
PRL 105, 150401 (2010) PHYSICAL REVIEW LETTERS 8 OCTOBER 2010
Oversampling Ratio: m 2dr r2 [1] In Quantum State Estimation, edited by M. Paris and
3 4 5 6 7 8 9 10 J. Řeháček, Lect. Notes Phys. Vol. 649 (Springer,
1. 1.
0.9 0.9 Heidelberg, 2004).
0.8 0.8 [2] H. Häffner et al., Nature (London) 438, 643 (2005).
Hybrid Random
0.7 0.7 [3] J. L. O’Brien et al., Nature (London) 426, 264 (2003).
Trace Distance
Fidelity [4] M. S. Kaznady and D. F. V. James, Phys. Rev. A 79,
0.6 0.6
Fidelity
150401-4