• Register
PhysicsOverflow is a next-generation academic platform for physicists and astronomers, including a community peer review system and a postgraduate-level discussion forum analogous to MathOverflow.

Welcome to PhysicsOverflow! PhysicsOverflow is an open platform for community peer review and graduate-level Physics discussion.

Please help promote PhysicsOverflow ads elsewhere if you like it.


PO is now at the Physics Department of Bielefeld University!

New printer friendly PO pages!

Migration to Bielefeld University was successful!

Please vote for this year's PhysicsOverflow ads!

Please do help out in categorising submissions. Submit a paper to PhysicsOverflow!

... see more

Tools for paper authors

Submit paper
Claim Paper Authorship

Tools for SE users

Search User
Reclaim SE Account
Request Account Merger
Nativise imported posts
Claim post (deleted users)
Import SE post

Users whose questions have been imported from Physics Stack Exchange, Theoretical Physics Stack Exchange, or any other Stack Exchange site are kindly requested to reclaim their account and not to register as a new user.

Public \(\beta\) tools

Report a bug with a feature
Request a new functionality
404 page design
Send feedback


(propose a free ad)

Site Statistics

205 submissions , 163 unreviewed
5,047 questions , 2,200 unanswered
5,345 answers , 22,709 comments
1,470 users with positive rep
816 active unimported users
More ...

  Multiqubit state tomography by performing measurement in the same basis

+ 13 like - 0 dislike

For a $n$-qubit state $\rho$ we perform all projective measurement consisting of one-particle measurements in the same basis, that is,

$$p_{i_1i_2\ldots i_n}(\theta,\varphi) = \text{Tr}\left \{ \rho [P_{i_1}(\theta,\varphi)\otimes P_{i_2}(\theta,\varphi) \otimes \ldots\otimes P_{i_n}(\theta,\varphi)] \right\},$$ where $P_0 = |\psi(\theta,\varphi)\rangle\langle \psi(\theta,\varphi) |$ and $P_1=\mathbb{I}-P_0$.

Does it suffice to recover the state?

If so, how can one do it?

This post has been migrated from (A51.SE)
asked Mar 15, 2012 in Theoretical Physics by Piotr Migdal (1,260 points) [ no revision ]
Have you checked whether it works for $n=2$?

This post has been migrated from (A51.SE)
@PeterShor Yes, it does work for $n=2$ (checked straightforwardly by looking at the Fourier components).

This post has been migrated from (A51.SE)
I made a sign error, making two equations linearly independent while they weren't. Sorry.

This post has been migrated from (A51.SE)

6 Answers

+ 6 like - 0 dislike

Working out an explicit case of Alex's construction, you get the two states $$ \rho_1 = \frac{1}{2} \left(|0\rangle\langle 0| \otimes |+ \rangle \langle + | + |1\rangle\langle 1| \otimes |- \rangle \langle - |\right) = \frac{1}{4}\left(\begin{array}{cccc}1 & 1 && \\ 1 & 1 & & \\ & & \phantom{-}1 & -1 \\ & & -1 & \phantom{-}1\end{array}\right)$$ and its image under the swap operator $$ \rho_2 = \frac{1}{2}\left(|+\rangle\langle +| \otimes |0 \rangle \langle 0 | + |-\rangle\langle -| \otimes |1 \rangle \langle 1 |\right) = \frac{1}{4}\left(\begin{array}{cccc}1 & &1& \\ & \phantom{-}1 & & -1\\ 1 & & 1 & \\ & -1& & \phantom{-}1\end{array}\right). $$

It's fairly easy to check that you cannot tell the difference between these. Let $|\psi\rangle = \alpha |0\rangle + \beta |1 \rangle$. Then $\langle \psi | = \alpha^* \langle 0 | + \beta^* \langle 1 |$ and $\langle \psi^\perp | = \beta \langle 0 | - \alpha \langle 1 |$.

Taking inner products, you find the probability of all the outcomes are equal for both $\rho_1$ and $\rho_2$. The outcomes are given by the following table: $$ \begin{array}{c c} \mathrm{outcome} &\mathrm{probability}\\ \langle \psi | \langle \psi | &\textstyle{\frac{1}{4}}\left(|\alpha|^2 |\alpha+\beta|^2 + |\beta|^2 |\alpha-\beta|^2\right)\\ \langle \psi | \langle \psi^\perp | &\textstyle{\frac{1}{4}}\left(|\alpha|^2 |\alpha-\beta|^2 + |\beta|^2 |\alpha+\beta|^2\right)\\ \langle \psi^\perp | \langle \psi | &\textstyle{\frac{1}{4}}\left(|\alpha|^2 |\alpha-\beta|^2 + |\beta|^2 |\alpha+\beta|^2\right)\\ \langle \psi^\perp | \langle \psi^\perp | &\textstyle{\frac{1}{4}}\left(|\alpha|^2 |\alpha+\beta|^2 + |\beta|^2 |\alpha-\beta|^2\right) \end{array} $$

It does seem possible that you can distinguish between any two pure states using the OP's measurements, but I haven't checked this carefully.

This post has been migrated from (A51.SE)
answered Mar 25, 2012 by Peter Shor (790 points) [ no revision ]
Indeed, I see - $\rho_1 = (1 + \sigma_z \otimes \sigma_x)/4$, $\rho_2 = (1 + \sigma_x \otimes \sigma_z)/4$.

This post has been migrated from (A51.SE)
+ 5 like - 0 dislike

It does not suffice for $n>1$. The idea is that density matrix for qubit may be expressed as $\frac{1}{2}\bigl(\hat{1}+ \sum_{k=1}^3 v_k \hat\sigma_k\bigr)$, there $v$ is a 3D real unit vector. Now we should consider span of tensor products of such terms. For any fixed $v$ we have $2^n$ tensor products of $\hat{P}_i = \bigl(\hat{1} + (-1)^i\sum_{k=1}^3 v_k \hat\sigma_k\bigr)/2$, $i =0,1$. Linear span of given products is equal with span generated by $2^n$ different combinations of tensor products of $2 \times 2$ unit matrix and (fixed for given $v$) matrix $\hat v = \sum_{k=1}^3 v_k \hat\sigma_k$. We want to know span of all such products for arbitrary $v$. This span is less than maximal $4^n$ due to linear dependence between coefficients (see example below).

Already consideration of the case $n=2$ shows, why it does not work. Four projectors relevant for this case may be written as $$\frac{1}{4}(\hat{1} \pm \sum_{k=1}^3 v_k \hat\sigma_k) \otimes (\hat{1} \pm \sum_{k=1}^3 v_k \hat\sigma_k)$$ so coefficients for $\hat\sigma_k \otimes \hat\sigma_j$ and $\hat\sigma_j \otimes \hat\sigma_k$ are equal for $k,j = 1,2,3$ in all four cases ($\pm v_k v_j$). They are equal for any $v_k$ if bases for both qubits are equal, so consideration of many different measurements does not change the property. So decomposition with Pauli matrices based on such terms always has three equal coefficients and may not represent decomposition of general state.

This post has been migrated from (A51.SE)
answered Mar 22, 2012 by Alex V (300 points) [ no revision ]
You can easily distinguish between the states $|01\rangle$ and $|10\rangle$ with OP's measurements. How is this compatible with your statement that coefficients for $\sigma_k \otimes \sigma_j$ and $\sigma_j \otimes \sigma_k$ are equal?

This post has been migrated from (A51.SE)
@Peter Shor: they are equal only for $k,j \neq 0$ but it is not so for terms with $\sigma_0 =1$ (so my note about "a swap operator" in earlier answer was misleading).

This post has been migrated from (A51.SE)
@Peter Shor: I already wrote, I do not require that - we may distinuish the swap of *all* components due to assymetry between $\sigma_0 \otimes \sigma_k$ and $\sigma_k \otimes \sigma_0$.

This post has been migrated from (A51.SE)
I know. The two states I gave were indeed an example of your construction, but it turns out I made a miscalculation ... you cannot distinguish between them, so your answer is correct. See my answer.

This post has been migrated from (A51.SE)
For me it is intuitive but not obvious: why out of _nonorthogonal_ measurements $P_i$ defined as in the question one, doing the reconstruction, need to constrain oneself to states $P_i$?

This post has been migrated from (A51.SE)
@Piotr Migdal: Using other words I am simply showing that the set of that nonorthogonal measurements is not "informationally complete".

This post has been migrated from (A51.SE)
+ 1 like - 0 dislike

If it can be done, then you will need at least $4^n$ angles $(\theta_k, \phi_k)$. However, the code below suggests that this method will generate at most $O\left(2^{\frac32n}\right)$ independent probabilities. Thus, it only works for $n=1$.

enter image description here

Suppose for the moment that it can be done; here is how to reconstruct the state from a set of angles. Arrange the angles into one index set $\Lambda = \{(\theta_k,\phi_k)| k\in\{0,1,2,\ldots,4^n\}\}$. So now you have an overcomplete (and necessarily non-orthogonal) basis of operators $\{P(\lambda) | \lambda\in \Lambda\}$. In other branches of mathematics, this thing is called a frame. To recover the state from the set of probabilities $\{p(\lambda)|\lambda\in\Lambda\}$ you need to compute a dual frame, call it $\{Q(\lambda)|\lambda\in\Lambda\}$. Then the state is given by $$ \rho = \sum_\lambda p(\lambda) Q(\lambda) \stackrel{\text{or}}{=} \sum_\lambda \text{Tr}(\rho Q(\lambda)) P(\lambda). $$ (Just like the $P$ and $Q$ function in quantum optics -- except in reverse notation.)

Unfortunately, you cannot compute a dual frame for an object that is not a frame (i.e. it doesn't span the operator space). The MATLAB code suggests that this is the case. It picks a bunch of random angles and computes the measurement operators for them. Then it checks if they are linearly independent and tries to compute a dual for them. Running it thousands of times produced no variation in the result that this procedure will produce roughly $O\left(2^{\frac32n}\right)$ linearly independent measurement operators.

% Number of qubits
n = 3;

% ===> Generate random qubit states
num_states = 2^n+2;

% Pick random Bloch sphere directions
phi = 2*pi*rand(1,num_states);
theta = acos(2*rand(1,num_states)-1);
state = [cos(theta/2) ; exp(1i*phi).*sin(theta/2)];
orth_state = [cos((pi-theta)/2) ; exp(1i*(phi+pi)).*sin((pi-theta)/2)];

% Now construct the tensor product states (there should be 2^n for each
% single qubit state since there are this many distinct outcomes) 
big_state = zeros(2^n,2^n*num_states);
for state_idx = 1:num_states
    for perm_idx = 1:2^n
        temp_state = 1;
        binrep = repmat(de2bi(perm_idx-1,n),2,1);
        to_tensor = repmat(state(:,state_idx),1,n).^binrep.*...
        for sys_idx = 1:n
            temp_state = kron(temp_state,to_tensor(:,sys_idx));
        big_state(:,(state_idx-1)*(2^n)+perm_idx) = temp_state;

% ===> Calculate the frame operator

% Vectorize the projectors onto the states
P = zeros(4^n,2^n*num_states);
for state_idx = 1:2^n*num_states
    P(:,state_idx) = kron(big_state(:,state_idx)...

% The frame operator
S = P*P';

% ===> Check for informational completeness and compute dual

% Compute the rank of S
fprintf('Rank of S is %d.\n', rank(S));

% The dual frame
Q = S\P;

% ===> Verify the reconstruction formula

% pick a random state
rho = eye(2^n)/2^n;
rho = rho(:);

% Compute reconstruction error
err = sum(abs(rho - ...
fprintf('Error in reconstruction was %d.\n', err);
This post has been migrated from (A51.SE)
answered Mar 15, 2012 by Chris Ferrie (660 points) [ no revision ]
Most voted comments show all comments
@ChrisFerrie I believe it to be true as well. I guess $\propto 2^n$ directions should suffice.

This post has been migrated from (A51.SE)
@PiotrMigdal - I wrote some code that suggests this method will produce at most $(n+1)^2$ independent parameters. Could you expand on your proof that it works for $n=2$?

This post has been migrated from (A51.SE)
@PiotrMigdal - I'm convinced that it doesn't work so I changed my answer. Maybe someone will come up with a more elegant proof.

This post has been migrated from (A51.SE)
If you just choose one basis, you get $2^n-1$ independent real parameters. Either your code is wrong, or I'm misinterpreting what you're saying.

This post has been migrated from (A51.SE)
@PeterShor - Yeah, the problem was that only the projector on to one outcome of the measurement was considered. The code is fixed to consider all $2^n$ outcomes -- but it still doesn't seem to work.

This post has been migrated from (A51.SE)
Most recent comments show all comments
Dummy comment: we are assuming that we can change the $(\theta,\phi)$ each time we measure, right?

This post has been migrated from (A51.SE)
@PiotrMigdal - OK, if I understand correctly, I'll claim you can choose $4^n-1$ random directions $(\theta_k,\phi_k)$ (or any non-colinear directions) and it should work.

This post has been migrated from (A51.SE)
+ 1 like - 0 dislike

It seems this cannot be done in general for the following reason: We can regard the process as a global rotation $G = R_X^{\otimes n}(\phi) R_Z^{\otimes n}(\theta)$ (i.e. a local rotation through the same angles on each qubit) followed by a $Z$ basis measurement. As $G$ is symmetric, it can be decomposed in terms of symmetric sums of Pauli operators. Since there is exactly one symmetric sum for each choice of $n_I$, $n_X$, $n_Y$ and $n_Z$, the numbers of $I$, $X$, $Y$ and $Z$ local operators in the Pauli terms. Since $n = n_I + n_X + n_Y + n_Z$, there are $\binom{n+3}{3}$ independent terms in including the identity. Since each of these can have a complex coefficient, there are $2\binom{n+3}{3}$ parameters. We have some constraints on these to enforce unitarity, and to avoid counting the global phase. I'll ignore these constraints here, so I will be over counting. Since we make a $Z$ measurement of each qubit and we can consider correlations from arbitrary combinations of these ($2^n - 1$ independent operators), at most we are determining $2(2^n - 1)\binom{n+3}{3}$ parameters.

Note that $$2(2^n - 1)\binom{n+3}{3} = \frac{2(2^n -1) (n+3)!}{n!3!} = \frac{1}{3}(2^n-1)(n+3)(n+2)(n+1)$$ however for an $n$ qubit state there are $4^n - 1$ free parameters, so the number of free parameters in the quantum state rapidly overtakes the maximum number of independent parameters you can measure with such measurements, since if you divide the former by the latter you get $\frac{(n+1)(n+2)(n+3)}{3(2^n+1)}$, which rapidly drops below $1$ (between $n=8$ and $n=9$).

This post has been migrated from (A51.SE)
answered Mar 26, 2012 by Joe Fitzsimons (3,575 points) [ no revision ]
Even measurement in single basis produces $2^n-1$ parameters. My preliminary estimation give equation $\sum_k C^n_k C^{k+2}_k$, but seems Chris' curve corresponds a bit different values.

This post has been migrated from (A51.SE)
@Alex: Ah I guess you mean in the correlations. Sorry I forgot to account for that. I'll fix it now.

This post has been migrated from (A51.SE)
+ 1 like - 0 dislike

I want to give you a general argument for the negative result: A general $n$-qubit density matrix cannot be contructed from this type of measurements.

The following observation will be needed:

Let $G$ be a compact Lie group acting on a vector space $V$ via a representation $\pi$. Let $\rho$ be a density matrix on $V$.

Consider the quantum characteristic function:

$\phi(g) = \mathrm{Tr}(\rho \pi(g)), \ \ g\in G$,

When $\pi$ is irreducible, then the density matrix can be reconstructed by virtue of the Peter-Weyl theorem:

$\rho_{ij} = \int \phi(g)\pi_{ij} (g) d\mu(g)\ \ i,j = 1,.,.,.,\mathrm{dim}(V)$

However, in the reducible case, there exists a basis in which all matrix representatives are block diagonal.Thus, in this basis, the matrix elements of the density matrix between basis vectors of different irreducible factors cannot be reconstructed, even given the quantum characteristic function at all points of the group manifold.

To apply that to our case consider the group $G = \bigotimes_{i=1}^{n} U(1) \bigotimes SU(2)$.

Where every Abelian factor acts on a different qubit by phase multiplication and the $SU(2)$ factor acts diagonally on the qubits.

Observing that an SU(2) group element acting on a qubit can be parameterized in terms of the Euler angles:

$g = e^{i\frac{\psi}{2}}P_1(\theta, \phi) +e^{-i\frac{\psi}{2}}P_0(\theta, \phi)$.

Then the quantum characteristic function can be obtained from the series of measurement expectations by a Fourier transform with respect to the Abelian coordinates and $\psi$ and vice versa.

But the group $G$ acts reducibly on the $n$-qubit, since every Abelian factor acts reducibly and $SU(2)$ acts irreducibly only on symmetric powers of qubits (thus reducibly on the full tensor power), which completes the argument.

This post has been migrated from (A51.SE)
answered Mar 26, 2012 by David Bar Moshe (4,355 points) [ no revision ]
+ 0 like - 0 dislike

I think that a concise way to understand why one-particle measurements cannot always recover any mixed state is to make use of the following fact: If a POVM does not span the space of operators on the Hilbert space in question, then there always exist two distinct states that give rise to exactly the same measurement probabilities.

(See http://www.cs.uwaterloo.ca/~watrous/CS766/ProblemSets/solutions1.pdf for a proof).

For a system of $n$ qubits, the space of operators has dimension $4^n$ but a one-particle measurement has only $2^n$ operators, so it cannot possibly span the space and hence not all mixed states can be distinguished unambiguously.

This post has been migrated from (A51.SE)
answered Apr 4, 2012 by Juan Miguel Arrazola (45 points) [ no revision ]

Your answer

Please use answers only to (at least partly) answer questions. To comment, discuss, or ask for clarification, leave a comment instead.
To mask links under text, please type your text, highlight it, and click the "link" button. You can then enter your link URL.
Please consult the FAQ for as to how to format your post.
This is the answer box; if you want to write a comment instead, please use the 'add comment' button.
Live preview (may slow down editor)   Preview
Your name to display (optional):
Privacy: Your email address will only be used for sending these notifications.
Anti-spam verification:
If you are a human please identify the position of the character covered by the symbol $\varnothing$ in the following word:
Then drag the red bullet below over the corresponding character of our banner. When you drop it there, the bullet changes to green (on slow internet connections after a few seconds).
Please complete the anti-spam verification

user contributions licensed under cc by-sa 3.0 with attribution required

Your rights