Skip to content

Bloch Sphere for Quantum Information

The Bloch ball is an exact three-real-parameter representation of every one-qubit density operator:

ρ=12(I+r⋅σ),r∈R3,∥r∥≤1.\begin{gathered} \rho = \frac12 \left( I+\mathbf r\cdot\boldsymbol\sigma \right), \\ \mathbf r\in\mathbb R^3, \qquad \lVert\mathbf r\rVert\leq1. \end{gathered}

The vector r\mathbf r is the Bloch vector. Pure states lie on the unit sphere ∥r∥=1\lVert\mathbf r\rVert=1, mixed states lie inside it, and the maximally mixed state is the center. Its components are experimentally meaningful Pauli expectation values:

r=(⟨X⟩, ⟨Y⟩, ⟨Z⟩).\mathbf r = \left( \langle X\rangle,\, \langle Y\rangle,\, \langle Z\rangle \right).

In quantum information, the picture is more than a state visualization. It is an operational dictionary:

  • a state corresponds to its vector r\mathbf r;
  • a projective measurement corresponds to an axis n^\hat{\mathbf n};
  • a one-qubit unitary corresponds to a proper rotation of r\mathbf r;
  • a one-qubit channel corresponds to an affine map r↦Tr+t\mathbf r\mapsto T\mathbf r+\mathbf t.

This page owns that operational QI dictionary. Bloch Sphere Geometry owns the coordinate geometry, Bloch Sphere for Spin owns the spinor and SU(2) interpretation, and Bloch Sphere for Density Operators derives the mixed-state ball from density-matrix positivity.

From a Density Matrix to Three Measurable Numbers

Section titled “From a Density Matrix to Three Measurable Numbers”

The identity and Pauli matrices form an orthogonal basis for 2×22\times2 operators under the Hilbert–Schmidt inner product:

Tr⁡(σiσj)=2δij,Tr⁡(σi)=0.\operatorname{Tr}(\sigma_i\sigma_j)=2\delta_{ij}, \qquad \operatorname{Tr}(\sigma_i)=0.

Consequently, every trace-one Hermitian matrix has a unique expansion

ρ=12(I+rxX+ryY+rzZ).\rho = \frac12 \left( I+r_xX+r_yY+r_zZ \right).

The coefficients are recovered by traces:

ri=Tr⁡(ρσi).r_i=\operatorname{Tr}(\rho\sigma_i).

In the computational basis,

ρ=12(1+rzrx−iryrx+iry1−rz).\rho = \frac12 \begin{pmatrix} 1+r_z&r_x-ir_y \\ r_x+ir_y&1-r_z \end{pmatrix}.

Thus rzr_z is the computational-basis population imbalance, while rxr_x and ryr_y encode the real and imaginary parts of the coherence. A basis label such as ZZ is a convention fixed by the chosen encoding and control frame. It need not be a literal spatial direction.

The eigenvalues of ρ\rho are

λ±=1±∥r∥2.\lambda_\pm = \frac{1\pm\lVert\mathbf r\rVert}{2}.

The condition ρ≥0\rho\geq0 is therefore exactly

∥r∥≤1.\lVert\mathbf r\rVert\leq1.

The purity is

Tr⁡(ρ2)=1+∥r∥22.\operatorname{Tr}(\rho^2) = \frac{1+\lVert\mathbf r\rVert^2}{2}.

It ranges from 1/21/2 at the center to 11 on the surface. The radius measures purity for a one-qubit state; it does not by itself identify a temperature, a noise source, or a preparation procedure.

After choosing the computational basis, a pure state can be parameterized as

∣ψ(θ,ϕ)⟩=cos⁡θ2∣0⟩+eiϕsin⁡θ2∣1⟩.|\psi(\theta,\phi)\rangle = \cos\frac{\theta}{2}|0\rangle + e^{i\phi} \sin\frac{\theta}{2}|1\rangle.

Its Bloch vector is

r=(sin⁡θcos⁡ϕ, sin⁡θsin⁡ϕ, cos⁡θ).\mathbf r = \left( \sin\theta\cos\phi,\, \sin\theta\sin\phi,\, \cos\theta \right).

The polar angle controls the population imbalance, and the azimuthal angle records relative phase. Global phase is absent because the density operator

∣ψ⟩⟨ψ∣|\psi\rangle\langle\psi|

is unchanged by ∣ψ⟩↦eiχ∣ψ⟩|\psi\rangle\mapsto e^{i\chi}|\psi\rangle. The derivation from normalized rays is given in Bloch Sphere for Spin and Projective Hilbert Space.

The eigenstates of XX, YY, and ZZ define six useful reference points.

DirectionStatePauli eigenvalue
+z+z∣0⟩\lvert0\rangleZ=+1Z=+1
−z-z∣1⟩\lvert1\rangleZ=−1Z=-1
+x+x∣+⟩=(∣0⟩+∣1⟩)/2\lvert+\rangle=(\lvert0\rangle+\lvert1\rangle)/\sqrt2X=+1X=+1
−x-x∣−⟩=(∣0⟩−∣1⟩)/2\lvert-\rangle=(\lvert0\rangle-\lvert1\rangle)/\sqrt2X=−1X=-1
+y+y∣+i⟩=(∣0⟩+i∣1⟩)/2\lvert+i\rangle=(\lvert0\rangle+i\lvert1\rangle)/\sqrt2Y=+1Y=+1
−y-y∣−i⟩=(∣0⟩−i∣1⟩)/2\lvert-i\rangle=(\lvert0\rangle-i\lvert1\rangle)/\sqrt2Y=−1Y=-1

These states are common calibration targets. The names xx, yy, and zz refer to Pauli coordinates in the logical frame. A hardware implementation must establish which pulses and readout settings realize those coordinates.

If preparations ρa\rho_a occur with classical probabilities pap_a, the averaged state is

ρ=∑apaρa.\rho=\sum_a p_a\rho_a.

Bloch vectors average in the same way:

r=∑apara.\mathbf r=\sum_a p_a\mathbf r_a.

The Bloch ball is therefore convex. Mixing surface points produces an interior point unless all weight lies on the same pure state.

The decomposition is not unique. For example,

I2=12(∣0⟩⟨0∣+∣1⟩⟨1∣)\frac I2 = \frac12 \left( |0\rangle\langle0| + |1\rangle\langle1| \right)

and

I2=12(∣+⟩⟨+∣+∣−⟩⟨−∣).\frac I2 = \frac12 \left( |+\rangle\langle+| + |-\rangle\langle-| \right).

Both ensembles give r=0\mathbf r=\mathbf0. A point in the ball specifies the density operator, not a unique story about which pure state “really occurred.” It can also be the reduced state of a qubit entangled with another system, with no ignorance ensemble preferred by the formalism.

Operational Bloch-sphere dictionary showing a state and measurement axis, a unitary rotation, and an affine channel deforming the ball

The Bloch representation organizes three one-qubit tasks. A projective measurement compares the state vector r\mathbf r with an axis n^\hat{\mathbf n}. A unitary rotates the ball rigidly, preserving radii and angles. A noisy channel acts affinely as r↦Tr+t\mathbf r\mapsto T\mathbf r+\mathbf t and can contract, distort, or translate the accessible region. Not every affine map taking the ball into itself is automatically a completely positive quantum channel.

For a unit vector n^\hat{\mathbf n}, the observable

σn^=n^⋅σ\sigma_{\hat{\mathbf n}} = \hat{\mathbf n}\cdot\boldsymbol\sigma

has eigenvalues ±1\pm1. Its projectors are

Π±(n^)=12(I±n^⋅σ).\Pi_\pm^{(\hat{\mathbf n})} = \frac12 \left( I\pm\hat{\mathbf n}\cdot\boldsymbol\sigma \right).

The Born rule gives

p±=Tr⁡(ρΠ±(n^))=12(1±n^⋅r).\begin{aligned} p_\pm &= \operatorname{Tr} \left( \rho\Pi_\pm^{(\hat{\mathbf n})} \right) \\ &= \frac12 \left( 1\pm\hat{\mathbf n}\cdot\mathbf r \right). \end{aligned}

The expectation value is the signed projection

⟨σn^⟩=p+−p−=n^⋅r.\langle\sigma_{\hat{\mathbf n}}\rangle = p_+-p_- = \hat{\mathbf n}\cdot\mathbf r.

This is the operational meaning of the angle between the state and measurement vectors. Parallel pure-state and measurement axes give a certain +1+1 outcome; antiparallel axes give a certain −1-1 outcome; perpendicular axes give equal probabilities.

State axes and measurement axes are different objects

Section titled “State axes and measurement axes are different objects”

The same three-dimensional coordinates describe a state vector and a measurement direction, but they play different roles. The state belongs to the preparation model. The axis belongs to the measurement effect. If both are unknown, observed probabilities alone do not determine which one is miscalibrated.

A real detector may be biased or unsharp. A binary effect can be expanded as

E+=αI+m⋅σ,E_+ = \alpha I+\mathbf m\cdot\boldsymbol\sigma,

with positivity constraints on α\alpha and m\mathbf m. The resulting probability

p(+∣ρ)=α+m⋅rp(+|\rho) = \alpha+\mathbf m\cdot\mathbf r

need not correspond to an ideal projective axis. Measurement Tomography owns reconstruction of unknown detector effects.

If ideal XX, YY, and ZZ measurements are calibrated, repeated outcomes estimate

rx=pX(+)−pX(−),ry=pY(+)−pY(−),rz=pZ(+)−pZ(−).\begin{aligned} r_x&=p_X(+)-p_X(-), \\ r_y&=p_Y(+)-p_Y(-), \\ r_z&=p_Z(+)-p_Z(-). \end{aligned}

For counts Ni,+N_{i,+} and Ni,−N_{i,-} in basis ii,

r^i=Ni,+−Ni,−Ni,++Ni,−.\widehat r_i = \frac{N_{i,+}-N_{i,-}} {N_{i,+}+N_{i,-}}.

These estimates fluctuate. Independent linear inversion can produce

∥r^∥>1,\lVert\widehat{\mathbf r}\rVert>1,

which corresponds to no positive density operator. This is not evidence for a state outside quantum mechanics; it is a sign that statistical noise and model assumptions must be handled. Constrained least squares, maximum likelihood, or Bayesian estimation can enforce physicality, but the estimator then has its own bias and uncertainty properties.

Tomography also inherits calibration error. If state preparation and measurement are both imperfect, a reconstructed vector can absorb errors from either side. A complete report states:

  • the preparation and measurement model;
  • shots per setting and uncertainty intervals;
  • drift monitoring and data-selection rules;
  • the estimator and positivity constraint;
  • readout correction or calibration procedure;
  • whether quoted fidelity is raw, corrected, or model dependent.

Three Pauli settings are informationally complete for one-qubit state tomography, but tomography scales poorly for many qubits. An arbitrary nn-qubit density matrix has 4n−14^n-1 real parameters. The single-qubit sphere does not remove that scaling.

Ignoring an overall phase, every one-qubit unitary can be written

U(n^,ϑ)=exp⁡(−iϑ2n^⋅σ).U(\hat{\mathbf n},\vartheta) = \exp \left( -\frac{i\vartheta}{2} \hat{\mathbf n}\cdot\boldsymbol\sigma \right).

Conjugation by UU maps

ρ′=UρU†\rho' = U\rho U^\dagger

to a real-space rotation of the Bloch vector:

r′=Rn^(ϑ)r.\mathbf r' = R_{\hat{\mathbf n}}(\vartheta)\mathbf r.

The rotation preserves

∥r′∥=∥r∥,\lVert\mathbf r'\rVert = \lVert\mathbf r\rVert,

so a unitary cannot purify or mix an isolated qubit. It moves pure states along the surface and mixed states on a sphere of fixed radius.

The factor of 1/21/2 in the exponential matters. A Bloch-vector rotation by angle ϑ\vartheta is generated by ϑn^⋅σ/2\vartheta\hat{\mathbf n}\cdot\boldsymbol\sigma/2. The underlying spinor changes sign under a 2π2\pi rotation, while its density operator and Bloch vector return after 2π2\pi. Spin Rotations owns the SU(2)-to-SO(3) derivation.

Up to global phase:

Gate familyBloch action
XXRotation by π\pi about xx
YYRotation by π\pi about yy
ZZRotation by π\pi about zz
Rx(ϑ)R_x(\vartheta)Rotation by ϑ\vartheta about xx
Ry(ϑ)R_y(\vartheta)Rotation by ϑ\vartheta about yy
Rz(ϑ)R_z(\vartheta)Rotation by ϑ\vartheta about zz
HHRotation by π\pi about (x^+z^)/2(\hat{\mathbf x}+\hat{\mathbf z})/\sqrt2, up to phase

Single-Qubit Gates owns gate matrices, phase conventions, Euler decompositions, native controls, and compilation distinctions. Here the sphere provides the geometric action.

One may rotate the state while holding the measurement axis fixed:

r↦Rr,\mathbf r\mapsto R\mathbf r,

or hold the state fixed and rotate the observable oppositely:

n^↦R−1n^.\hat{\mathbf n}\mapsto R^{-1}\hat{\mathbf n}.

The probability is unchanged because

n^⋅Rr=(R−1n^)⋅r.\hat{\mathbf n}\cdot R\mathbf r = \left( R^{-1}\hat{\mathbf n} \right)\cdot\mathbf r.

This distinction is practical. A laboratory may implement an XX measurement by applying a basis-change pulse and then using a fixed ZZ readout. Whether the diagram shows a rotated state or rotated measurement axis depends on the chosen description.

Rotations generally do not commute:

Rx(α)Ry(β)≠Ry(β)Rx(α).R_x(\alpha)R_y(\beta) \neq R_y(\beta)R_x(\alpha).

On a state column, the rightmost gate acts first. Bloch diagrams can hide this convention when arrows are drawn without an explicit time order. For pulse sequences, state the active rotation convention, axis orientation, and multiplication order.

Every trace-preserving linear map on one-qubit operators acts affinely on the Bloch vector:

r′=Tr+t,\mathbf r' = T\mathbf r+\mathbf t,

where TT is a real 3×33\times3 matrix and t∈R3\mathbf t\in\mathbb R^3. A physical quantum channel must be completely positive as well as trace preserving; this imposes constraints stronger than merely mapping the Bloch ball into itself.

A channel is unital when

E(I)=I.\mathcal E(I)=I.

In Bloch coordinates this means

t=0.\mathbf t=\mathbf0.

Unital channels leave the center fixed, though they may contract or anisotropically distort the ball. Nonunital channels translate the center.

Define the dephasing family

Dλ(ρ)=1+λ2ρ+1−λ2ZρZ,−1≤λ≤1.\begin{aligned} \mathcal D_\lambda(\rho) &= \frac{1+\lambda}{2}\rho + \frac{1-\lambda}{2}Z\rho Z, \\ -1\leq\lambda\leq1. \end{aligned}

Its Bloch action is

(rx,ry,rz)⟼(λrx,λry,rz).(r_x,r_y,r_z) \longmapsto (\lambda r_x,\lambda r_y,r_z).

At λ=0\lambda=0, transverse coherence vanishes and the sphere collapses to the zz-axis segment. Populations are unchanged. Negative λ\lambda includes an additional phase inversion; in many dynamical dephasing models λ\lambda decays from 11 toward 00. Dephasing Channel owns the Kraus, master-equation, and coherence-time treatments.

For the convention

Pq(ρ)=(1−q)ρ+qI2,0≤q≤1,\mathcal P_q(\rho) = (1-q)\rho+q\frac I2, \qquad 0\leq q\leq1,

the Bloch vector transforms isotropically:

r⟼(1−q)r.\mathbf r\longmapsto(1-q)\mathbf r.

The ball contracts toward the center. Other sources use different depolarizing parameters, so the channel definition should accompany any quoted error probability.

With ∣0⟩|0\rangle chosen as the ground state and damping probability γ\gamma,

rx′=1−γ rx,ry′=1−γ ry,rz′=(1−γ)rz+γ.\begin{aligned} r_x'&=\sqrt{1-\gamma}\,r_x, \\ r_y'&=\sqrt{1-\gamma}\,r_y, \\ r_z'&=(1-\gamma)r_z+\gamma. \end{aligned}

This channel contracts and translates the ball toward +z+z. It is not unital because

I2⟼12(I+γZ).\frac I2 \longmapsto \frac12 \left( I+\gamma Z \right).

Amplitude-Damping Channel gives the canonical channel derivation and physical interpretation.

A picture is not a complete process certification

Section titled “A picture is not a complete process certification”

Drawing an ellipsoid inside the ball is useful but insufficient to certify a channel. Complete positivity constrains the allowed contractions and translations, and correlated initial states can invalidate a state-independent reduced channel model. Process tomography additionally depends on trusted preparations and measurements. The Quantum Channels and Noise volume section owns Kraus, Choi, and Stinespring representations and their consistency conditions.

Several one-qubit quantities become elementary functions of Bloch vectors.

For qubit states ρ\rho and σ\sigma with vectors r\mathbf r and s\mathbf s,

D(ρ,σ)=12∥r−s∥.D(\rho,\sigma) = \frac12 \lVert\mathbf r-\mathbf s\rVert.

Thus optimal one-shot distinguishability is related to Euclidean separation in the ball. Antipodal pure states have distance 11 and are orthogonal. The canonical definition and operational theorem are in Trace Distance.

Using the squared Uhlmann-fidelity convention,

F(ρ,σ)=1+r⋅s2+1−∥r∥21−∥s∥22.\begin{aligned} F(\rho,\sigma) &= \frac{ 1+\mathbf r\cdot\mathbf s }{2} \\ &\quad + \frac{ \sqrt{1-\lVert\mathbf r\rVert^2} \sqrt{1-\lVert\mathbf s\rVert^2} }{2}. \end{aligned}

If ρ=∣ψ⟩⟨ψ∣\rho=|\psi\rangle\langle\psi| is pure, this reduces to

F(ρ,σ)=⟨ψ∣σ∣ψ⟩=1+r⋅s2.F(\rho,\sigma) = \langle\psi|\sigma|\psi\rangle = \frac{1+\mathbf r\cdot\mathbf s}{2}.

Some sources call F\sqrt F the fidelity, so conventions must be stated. Fidelity fixes the site’s formula convention.

Because the eigenvalues are (1±r)/2(1\pm r)/2 with r=∥r∥r=\lVert\mathbf r\rVert, the von Neumann entropy is

S(ρ)=h2(1+r2),S(\rho) = h_2 \left( \frac{1+r}{2} \right),

where h2h_2 is binary entropy. Entropy decreases monotonically from one bit at the center to zero on the surface. This simple radial dependence is special to one qubit.

Prepare ∣0⟩|0\rangle, apply a π/2\pi/2 pulse about yy, and obtain the +x+x state. A relative phase accumulation

Uϕ=exp⁡(−iϕZ2)U_\phi = \exp \left( -\frac{i\phi Z}{2} \right)

rotates the vector to

rϕ=(cos⁡ϕ,sin⁡ϕ,0).\mathbf r_\phi = (\cos\phi,\sin\phi,0).

An XX measurement gives

pX(+∣ϕ)=1+cos⁡ϕ2,p_X(+|\phi) = \frac{1+\cos\phi}{2},

while a YY measurement gives

pY(+∣ϕ)=1+sin⁡ϕ2.p_Y(+|\phi) = \frac{1+\sin\phi}{2}.

The phase is invisible to a direct ZZ measurement because rz=0r_z=0 throughout. This is why a final analysis pulse is needed to convert phase into a population difference.

If the same equatorial state undergoes Dλ\mathcal D_\lambda,

rϕ⟼λ(cos⁡ϕ,sin⁡ϕ,0).\mathbf r_\phi \longmapsto \lambda (\cos\phi,\sin\phi,0).

Ramsey fringe contrast is multiplied by λ\lambda. The radial contraction summarizes the loss of coherence for this state family, but it does not by itself identify whether the cause was stochastic detuning, entanglement with an environment, averaging over drift, or another mechanism.

Suppose a prepared state has

r=(0.6,0,0.8),\mathbf r=(0.6,0,0.8),

which is pure because ∥r∥=1\lVert\mathbf r\rVert=1. Measuring along

n^=x^+z^2\hat{\mathbf n} = \frac{ \hat{\mathbf x}+\hat{\mathbf z} }{\sqrt2}

gives

p+=12[1+1.42]≈0.995.p_+ = \frac12 \left[ 1+\frac{1.4}{\sqrt2} \right] \approx0.995.

The near-certain outcome reflects geometric alignment, not simultaneous pre-existing XX and ZZ values.

An nn-qubit density matrix requires 4n−14^n-1 real parameters. One may draw a reduced Bloch vector for each qubit, but those vectors omit correlations and entanglement. For the Bell state

∣Φ+⟩=∣00⟩+∣11⟩2,|\Phi^+\rangle = \frac{|00\rangle+|11\rangle}{\sqrt2},

each qubit has r=0\mathbf r=\mathbf0, yet the joint state is pure and maximally entangled. Two centered local balls do not describe the global state.

Generalized Bloch vectors are not ordinary balls

Section titled “Generalized Bloch vectors are not ordinary balls”

For a qudit, one can expand a state in d2−1d^2-1 traceless generators. Positivity does not fill a simple Euclidean ball when d>2d>2. The elegant equivalence

ρ≥0⟺∥r∥≤1\rho\geq0 \quad\Longleftrightarrow\quad \lVert\mathbf r\rVert\leq1

is special to d=2d=2 under the standard normalization.

The zz axis usually names the computational basis, but whether ∣0⟩|0\rangle is north or south, ground or excited, horizontal or vertical polarization is conventional. Rotating-frame phases determine laboratory xx and yy. A diagram without those conventions may reverse signs or rotation directions.

A trajectory does not identify a mechanism

Section titled “A trajectory does not identify a mechanism”

The same visible path can arise from a Hamiltonian, a time-dependent control frame, postselection, or a noisy channel. A shrinking vector signals reduced one-qubit purity, but not a unique microscopic environment. Mechanistic claims require a declared dynamical model and independent tests.

Interior points do not reveal which classical mixture produced them. They may also be reduced states of entangled systems. The ball captures all one-qubit measurement statistics, not a hidden decomposition.

The sphere does not make noncommuting observables classical

Section titled “The sphere does not make noncommuting observables classical”

The coordinates rxr_x, ryr_y, and rzr_z are expectation values measurable on separate, similarly prepared systems. They are not three jointly readable values carried by one specimen. Uncertainty and measurement disturbance remain.

  • Calling every point in the ball a pure state; only the surface is pure.
  • Treating the Bloch vector as the state vector in C2\mathbb C^2.
  • Forgetting that global phase is absent but relative phase sets the azimuth.
  • Interpreting logical xx, yy, and zz as literal laboratory directions without an encoding map.
  • Using r↦Tr+t\mathbf r\mapsto T\mathbf r+\mathbf t without checking complete positivity.
  • Inferring a calibrated state from data when the measurement axes are themselves unknown.
  • Comparing depolarizing or dephasing parameters without checking conventions.
  • Using local Bloch vectors as a complete description of a multi-qubit state.
  • Treating linear-inversion estimates with ∥r^∥>1\lVert\widehat{\mathbf r}\rVert>1 as physical states.
  • Assuming a shrinking vector identifies one unique decoherence mechanism.

Ideal measurements give

pX(+)=0.65,pY(+)=0.30,pZ(+)=0.75.\begin{aligned} p_X(+)&=0.65, \\ p_Y(+)&=0.30, \\ p_Z(+)&=0.75. \end{aligned}

Find the Bloch vector and density matrix. Is the estimate physical?

Solution

For a Pauli measurement,

ri=2pi(+)−1.r_i=2p_i(+)-1.

Therefore

r=(0.30,−0.40,0.50).\mathbf r=(0.30,-0.40,0.50).

Its length is

∥r∥=0.09+0.16+0.25=0.50<1,\begin{aligned} \lVert\mathbf r\rVert &= \sqrt{0.09+0.16+0.25} \\ &= \sqrt{0.50} <1, \end{aligned}

so it is physical. The density matrix is

ρ=12(1.500.30+0.40i0.30−0.40i0.50).\rho = \frac12 \begin{pmatrix} 1.50&0.30+0.40i \\ 0.30-0.40i&0.50 \end{pmatrix}.

Equivalently,

ρ=(0.750.15+0.20i0.15−0.20i0.25).\rho = \begin{pmatrix} 0.75&0.15+0.20i \\ 0.15-0.20i&0.25 \end{pmatrix}.

Tomographic linear inversion gives

r^=(0.8,0.6,0.3).\widehat{\mathbf r}=(0.8,0.6,0.3).

Why is this not a valid qubit state, and what should an analyst do?

Solution

The squared length is

∥r^∥2=0.82+0.62+0.32=1.09>1.\lVert\widehat{\mathbf r}\rVert^2 = 0.8^2+0.6^2+0.3^2 = 1.09>1.

The smaller eigenvalue

λ^−=1−1.092\widehat\lambda_- = \frac{ 1-\sqrt{1.09} }{2}

is negative, so the reconstructed matrix is not positive semidefinite.

The analyst should not simply interpret the vector as “more than pure.” They should retain the raw counts, quantify statistical and calibration uncertainty, and use a documented physical estimator such as constrained least squares, maximum likelihood, or a Bayesian method. The choice of estimator and correction procedure must be reported.

A qubit is in the +x+x state, r=(1,0,0)\mathbf r=(1,0,0). It is measured along

n^=x^+z^2.\hat{\mathbf n} = \frac{ \hat{\mathbf x}+\hat{\mathbf z} }{\sqrt2}.

Find p+p_+ and p−p_-.

Solution

The dot product is

n^⋅r=12.\hat{\mathbf n}\cdot\mathbf r = \frac1{\sqrt2}.

Hence

p+=12(1+12),p−=12(1−12).\begin{aligned} p_+ &= \frac12 \left( 1+\frac1{\sqrt2} \right), \\ p_- &= \frac12 \left( 1-\frac1{\sqrt2} \right). \end{aligned}

Numerically, p+≈0.854p_+\approx0.854 and p−≈0.146p_-\approx0.146.

The initial Bloch vector is r=(1,0,0)\mathbf r=(1,0,0). Apply

Rz(π2)=exp⁡(−iπZ4).R_z\left(\frac\pi2\right) = \exp\left(-\frac{i\pi Z}{4}\right).

Find the final state and the probabilities of XX, YY, and ZZ outcomes.

Solution

A positive π/2\pi/2 rotation about zz sends +x+x to +y+y:

r′=(0,1,0).\mathbf r'=(0,1,0).

The final state is

∣+i⟩=∣0⟩+i∣1⟩2|+i\rangle = \frac{|0\rangle+i|1\rangle}{\sqrt2}

up to global phase. Therefore

pX(+)=1/2,pY(+)=1,pZ(+)=1/2.\begin{aligned} p_X(+)&=1/2, \\ p_Y(+)&=1, \\ p_Z(+)&=1/2. \end{aligned}

The +x+x state passes through Dλ\mathcal D_\lambda with 0≤λ≤10\leq\lambda\leq1. Find its output Bloch vector and purity. Evaluate the result at λ=0\lambda=0 and λ=1\lambda=1.

Solution

The input vector is (1,0,0)(1,0,0), so

r′=(λ,0,0).\mathbf r'=(\lambda,0,0).

Its purity is

Tr⁡(ρ′2)=1+λ22.\operatorname{Tr}(\rho'^2) = \frac{1+\lambda^2}{2}.

At λ=1\lambda=1, the state remains pure with purity 11. At λ=0\lambda=0, the output is the center of the equatorial chord,

ρ′=I2,\rho'=\frac I2,

with purity 1/21/2. Complete ZZ dephasing of the equal superposition produces an equal incoherent mixture of ∣0⟩|0\rangle and ∣1⟩|1\rangle.

6. Show that amplitude damping is nonunital

Section titled “6. Show that amplitude damping is nonunital”

Apply amplitude damping with probability γ\gamma to the maximally mixed state. Find the output vector and explain why the result proves the channel is nonunital.

Solution

The maximally mixed input has

r=0.\mathbf r=\mathbf0.

Using the affine map,

r′=(0,0,γ).\mathbf r' = (0,0,\gamma).

Thus

Aγ(I2)=12(I+γZ).\mathcal A_\gamma \left( \frac I2 \right) = \frac12 \left( I+\gamma Z \right).

For γ>0\gamma>0, this is not I/2I/2. Equivalently,

Aγ(I)≠I.\mathcal A_\gamma(I)\neq I.

The channel translates the center toward the ground-state pole and is nonunital.

  • F. Bloch, “Nuclear induction,” Physical Review 70, 460–474, 1946, doi:10.1103/PhysRev.70.460.
  • U. Fano, “Description of states in quantum mechanics by density matrix and operator techniques,” Reviews of Modern Physics 29, 74–93, 1957, doi:10.1103/RevModPhys.29.74.
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, 10th anniversary ed., Cambridge University Press, 2010, doi:10.1017/CBO9780511976667.
  • J. Watrous, The Theory of Quantum Information, Cambridge University Press, 2018, doi:10.1017/9781316848142.
  • I. Bengtsson and K. Życzkowski, Geometry of Quantum States, 2nd ed., Cambridge University Press, 2017, doi:10.1017/9781139207010.
  • J. Preskill, Lecture Notes for Physics 219/Computer Science 219: Quantum Computation, California Institute of Technology, course materials.
  • D. F. V. James, P. G. Kwiat, W. J. Munro, and A. G. White, “Measurement of qubits,” Physical Review A 64, 052312, 2001, doi:10.1103/PhysRevA.64.052312.
  • M. G. A. Paris and J. Řeháček, eds., Quantum State Estimation, Springer, 2004, doi:10.1007/b98673.
  • M. B. Ruskai, S. Szarek, and E. Werner, “An analysis of completely-positive trace-preserving maps on M2M_2,” Linear Algebra and its Applications 347, 159–187, 2002, doi:10.1016/S0024-3795(01)00547-X.
  • A. Fujiwara and P. Algoet, “One-to-one parametrization of quantum channels,” Physical Review A 59, 3290–3294, 1999, doi:10.1103/PhysRevA.59.3290.
  • A. Gilchrist, N. K. Langford, and M. A. Nielsen, “Distance measures to compare real and ideal quantum processes,” Physical Review A 71, 062310, 2005, doi:10.1103/PhysRevA.71.062310.
  • R. Blume-Kohout, “Optimal, reliable estimation of quantum states,” New Journal of Physics 12, 043034, 2010, doi:10.1088/1367-2630/12/4/043034.
  • S. T. Merkel, J. M. Gambetta, J. A. Smolin, S. Poletto, A. D. Córcoles, B. R. Johnson, C. A. Ryan, and M. Steffen, “Self-consistent quantum process tomography,” Physical Review A 87, 062119, 2013, doi:10.1103/PhysRevA.87.062119.

The Bloch ball gives an exact operational coordinate system for one qubit. A state is represented by r=(⟨X⟩,⟨Y⟩,⟨Z⟩)\mathbf r=(\langle X\rangle,\langle Y\rangle,\langle Z\rangle) with ∥r∥≤1\lVert\mathbf r\rVert\leq1; pure states lie on the surface, and mixed states lie inside. An ideal projective measurement along n^\hat{\mathbf n} has probabilities (1±n^⋅r)/2(1\pm\hat{\mathbf n}\cdot\mathbf r)/2. Three calibrated Pauli measurements reconstruct the state, subject to finite-sample and calibration uncertainty.

One-qubit unitaries rotate the ball rigidly. General channels act affinely as r↦Tr+t\mathbf r\mapsto T\mathbf r+\mathbf t, with complete positivity constraining the allowed deformation. Dephasing contracts transverse components, depolarization contracts isotropically, and amplitude damping contracts and translates toward the ground-state pole. The picture is exact for one qubit but does not encode multi-qubit correlations, choose an ensemble decomposition, certify a channel, or replace calibration and uncertainty analysis.