The Bloch Sphere

A single qubit lives in a two-dimensional complex vector space. Yet its physical state can be drawn as a single point on an ordinary sphere. This is not a metaphor — it is an exact, faithful geometric picture. Every pure state of a qubit corresponds to exactly one point on the surface of a unit sphere. This holds no matter how the state was prepared. Every point on that sphere corresponds to exactly one pure qubit state. This correspondence holds up to an unobservable global phase. This picture matters because quantum algorithms are usually described as sequences of gates acting on qubits. Gates are unitary matrices acting on complex amplitudes. That algebra is correct but visually opaque. The Bloch sphere turns matrix multiplication into rotation. It turns superposition into geometric position. It turns "how much does this gate change the state" into "how far does this point move." Hardware engineers use it daily to reason about pulse sequences on real qubits. Textbooks use it to make single-qubit gates memorable and provable at a glance. The big picture is a change of coordinates. A qubit state is normally written as a pair of complex numbers subject to a normalization constraint. This leaves two real degrees of freedom after removing an overall phase. Those two real numbers are exactly the latitude and longitude of a point on a sphere. This article makes that correspondence precise. It derives the polar and azimuthal angles from the amplitudes. It explains why global phase disappears and relative phase survives as the sphere's longitude. It also gives worked numerical examples you can check by hand. By the end you will be able to convert any single-qubit state vector into Bloch angles and back. You will be able to explain geometrically why ∣0⟩|0\rangle and ∣1⟩|1\rangle sit at opposite poles. You will also see why this picture is the natural stage for describing Single-Qubit Gates: X, Y, Z, H, Rotation Gates and Arbitrary Single-Qubit Unitaries as sphere rotations.

Key takeaways
  • Any pure single-qubit state can be written ∣ψ⟩=cos⁡(θ/2)∣0⟩+eiφsin⁡(θ/2)∣1⟩|\psi\rangle = \cos(\theta/2)|0\rangle + e^{i\varphi}\sin(\theta/2)|1\rangle.
  • (θ,φ)(\theta, \varphi) are spherical coordinates: θ\theta is the polar (latitude) angle, φ\varphi is the azimuthal (longitude) angle.
  • Global phase (an overall factor eiγe^{i\gamma}) is physically invisible and is discarded before mapping to the sphere.
  • Relative phase between the ∣0⟩|0\rangle and ∣1⟩|1\rangle amplitudes is physically meaningful and becomes the sphere's longitude φ\varphi.
  • The poles are ∣0⟩|0\rangle (θ=0\theta=0) and ∣1⟩|1\rangle (θ=π\theta=\pi); the equator holds equal-superposition states like ∣+⟩|+\rangle, ∣−⟩|-\rangle, ∣+i⟩|+i\rangle, ∣−i⟩|-i\rangle.
  • This representation only works for pure states of a single qubit; mixed states live inside the sphere, and multi-qubit states need a different (higher-dimensional) tool.
The core insight

A pure single-qubit state is fully and faithfully described by two angles — its Bloch sphere coordinates — because global phase is unobservable and normalization removes one more real degree of freedom.

The idea

A general single-qubit pure state is a unit vector in C2\mathbb{C}^2:

∣ψ⟩=α∣0⟩+β∣1⟩,∣α∣2+∣β∣2=1,|\psi\rangle = \alpha|0\rangle + \beta|1\rangle, \qquad |\alpha|^2 + |\beta|^2 = 1,

with α,β∈C\alpha,\beta \in \mathbb{C}. Naively this state is specified by four real numbers (two real and two imaginary parts). The normalization constraint removes one real degree of freedom, leaving three. But quantum mechanics postulates that two state vectors differing only by an overall complex phase, ∣ψ⟩|\psi\rangle and eiγ∣ψ⟩e^{i\gamma}|\psi\rangle for real γ\gamma, are physically indistinguishable: every measurement probability, computed via ∣⟨ϕ∣ψ⟩∣2|\langle\phi|\psi\rangle|^2, is unchanged by such a factor. This is the global phase, and it removes one further real degree of freedom, leaving exactly two real parameters. Two real parameters are precisely what is needed to label a point on the surface of a sphere: a polar angle θ∈[0,π]\theta \in [0,\pi] and an azimuthal angle φ∈[0,2π)\varphi \in [0,2\pi), exactly as for latitude and longitude on a globe. Writing the amplitudes in this parametrization,

∣ψ⟩=cos⁡ ⁣(θ2)∣0⟩+eiφsin⁡ ⁣(θ2)∣1⟩,|\psi\rangle = \cos\!\left(\frac{\theta}{2}\right)|0\rangle + e^{i\varphi}\sin\!\left(\frac{\theta}{2}\right)|1\rangle,

is called the Bloch sphere representation. The half-angle θ/2\theta/2 is a deliberate convention: it makes θ\theta range over the full [0,π][0,\pi] while ∣0⟩|0\rangle and ∣1⟩|1\rangle, which are physically orthogonal (perfectly distinguishable), end up diametrically opposite on the sphere (θ=0\theta = 0 and θ=π\theta = \pi) rather than at right angles. This matches physical intuition better than the raw amplitude angle would. The key conceptual split is between global phase and relative phase. Global phase multiplies the entire state vector, eiγ∣ψ⟩e^{i\gamma}|\psi\rangle, and is unobservable — it never appears in the Bloch angles at all. Relative phase is the phase difference between the ∣1⟩|1\rangle amplitude and the ∣0⟩|0\rangle amplitude; after factoring out global phase to make the ∣0⟩|0\rangle coefficient real and non-negative, this relative phase survives as exactly the angle φ\varphi. It is directly observable, for instance by measuring in the XX or YY basis, and it is what distinguishes states like ∣+⟩=12(∣0⟩+∣1⟩)|+\rangle = \frac{1}{\sqrt2}(|0\rangle+|1\rangle) from ∣−⟩=12(∣0⟩−∣1⟩)|-\rangle = \frac{1}{\sqrt2}(|0\rangle - |1\rangle), which have the same θ\theta but different φ\varphi.

How it works

Step 1 — start from amplitudes. Given ∣ψ⟩=α∣0⟩+β∣1⟩|\psi\rangle = \alpha|0\rangle + \beta|1\rangle with ∣α∣2+∣β∣2=1|\alpha|^2+|\beta|^2=1, write α=∣α∣eiγα\alpha = |\alpha|e^{i\gamma_\alpha} and β=∣β∣eiγβ\beta = |\beta|e^{i\gamma_\beta} in polar form.

Step 2 — strip global phase. Multiply the whole state by e−iγαe^{-i\gamma_\alpha}. This does not change the physical state, but now the ∣0⟩|0\rangle coefficient is real and non-negative: α′=∣α∣\alpha' = |\alpha|, and β′=∣β∣ei(γβ−γα)\beta' = |\beta|e^{i(\gamma_\beta - \gamma_\alpha)}. Define φ=γβ−γα(mod2π)\varphi = \gamma_\beta - \gamma_\alpha \pmod {2\pi} — this is the relative phase.

Step 3 — read off angles. Since ∣α∣2+∣β∣2=1|\alpha|^2 + |\beta|^2 = 1, there is a unique θ∈[0,π]\theta \in [0,\pi] with

∣α∣=cos⁡(θ/2),∣β∣=sin⁡(θ/2).|\alpha| = \cos(\theta/2), \qquad |\beta| = \sin(\theta/2).

Concretely, θ=2arccos⁡(∣α∣)\theta = 2\arccos(|\alpha|) (equivalently θ=2arcsin⁡(∣β∣)\theta = 2\arcsin(|\beta|), with the branch fixed by θ∈[0,π]\theta\in[0,\pi]).

Step 4 — plot the point. (θ,φ)(\theta,\varphi) are spherical coordinates on a unit sphere. The corresponding Cartesian (Bloch vector) coordinates are:

x=sin⁡θcos⁡φ,y=sin⁡θsin⁡φ,z=cos⁡θ.x = \sin\theta\cos\varphi, \qquad y = \sin\theta\sin\varphi, \qquad z = \cos\theta.

These are exactly the expectation values of the Pauli operators: x=⟨ψ∣X∣ψ⟩x = \langle\psi|X|\psi\rangle, y=⟨ψ∣Y∣ψ⟩y=\langle\psi|Y|\psi\rangle, z=⟨ψ∣Z∣ψ⟩z=\langle\psi|Z|\psi\rangle, which is why (x,y,z)(x,y,z) is called the Bloch vector.

Step 5 (inverse) — angles back to amplitudes. Given (θ,φ)(\theta,\varphi), reconstruct

∣ψ⟩=cos⁡(θ/2)∣0⟩+eiφsin⁡(θ/2)∣1⟩,|\psi\rangle = \cos(\theta/2)|0\rangle + e^{i\varphi}\sin(\theta/2)|1\rangle,

which is unique up to the (physically irrelevant) global phase you discarded in Step 2. ``` FUNCTION StateToBlochAngles(alpha, beta):

# alpha, beta are complex numbers, |alpha|^2 + |beta|^2 = 1 ASSERT |alpha|^2 + |beta|^2 ≈ 1 # normalization check r_alpha = |alpha| r_beta = |beta| theta = 2 arccos(clamp(r_alpha, -1, 1)) # theta in [0, pi] IF r_alpha < EPS: # alpha ≈ 0 (pure |1>) phi = arg(beta) # any global phase convention ELSE IF r_beta < EPS: # beta ≈ 0 (pure |0>) phi = 0 # azimuth undefined at pole; pick 0 ELSE: phi = arg(beta) - arg(alpha) # relative phase phi = phi mod 2pi RETURN (theta, phi) FUNCTION BlochAnglesToState(theta, phi): alpha = cos(theta / 2) beta = exp(i phi) sin(theta / 2) RETURN (alpha, beta) # canonical representative,

# alpha real and ≥ 0 ``` The pseudocode's edge-case handling matters: at the poles (θ=0\theta=0 or θ=π\theta=\pi) one of α,β\alpha,\beta is exactly zero, so the relative phase φ\varphi is geometrically meaningless (longitude is undefined at the poles of a globe too) — the function returns a conventional value rather than dividing by zero or taking `arg` of `0`.

Figure
The Bloch sphere: axes, poles, and equatorial states
|0⟩ (z=+1) |1⟩ (z=−1) x −x y |+⟩ (φ=0) |−⟩ (φ=π) |+i⟩ (φ=π/2) |−i⟩ (φ=3π/2) |ψ⟩ θ Poles: computational basis. Equator: equal superpositions.

Complexity & correctness

"Complexity" here is representational, not computational. Converting between amplitudes (α,β)(\alpha,\beta) and angles (θ,φ)(\theta,\varphi) costs a constant number of arithmetic operations: one `arccos`, one or two `arg` (complex argument / `atan2`) calls, and a subtraction. There are no loops, no data structures, and no dependence on any input size — it is O(1)O(1) time and O(1)O(1) space per qubit, because a single qubit's pure state is described by a fixed, small number of real parameters regardless of context.

Correctness — this is a bijection (up to global phase). The argument has three parts. Existence. Given any valid amplitudes α,β\alpha,\beta with ∣α∣2+∣β∣2=1|\alpha|^2+|\beta|^2=1, Step 2–3 always produce a valid θ∈[0,π]\theta \in [0,\pi]: this follows because ∣α∣∈[0,1]|\alpha| \in [0,1] by normalization, and arccos⁡\arccos is defined and single-valued on [−1,1][-1,1], mapping onto [0,π][0,\pi]. The value φ\varphi always exists modulo 2π2\pi as the difference of two arguments (except at the measure-zero poles, handled by convention). Uniqueness up to global phase. Suppose two amplitude pairs (α,β)(\alpha,\beta) and (α′,β′)(\alpha',\beta') give the same (θ,φ)(\theta,\varphi) via the canonical form (real non-negative α\alpha). Then α=α′=cos⁡(θ/2)\alpha=\alpha'=\cos(\theta/2) and β=β′=eiφsin⁡(θ/2)\beta=\beta'=e^{i\varphi}\sin(\theta/2) exactly, since the canonicalization removes all phase freedom. Conversely, if (α,β)(\alpha,\beta) and (α′,β′)(\alpha',\beta') describe the same physical state, they differ by a global phase eiγe^{i\gamma}; Step 2 strips exactly this factor by construction (it multiplies by e−iγαe^{-i\gamma_\alpha}, whatever that phase is), so both canonicalize to the identical (θ,φ)(\theta,\varphi). Hence the map from physical states (equivalence classes under global phase) to points (θ,φ)(\theta,\varphi) is a genuine bijection. Surjectivity onto the sphere. Every (θ,φ)∈[0,π]×[0,2π)(\theta,\varphi) \in [0,\pi]\times[0,2\pi) arises from some valid normalized state, namely the one produced by `BlochAnglesToState`, and one checks directly that ∣cos⁡(θ/2)∣2+∣sin⁡(θ/2)∣2=1|\cos(\theta/2)|^2 + |\sin(\theta/2)|^2 = 1 for all θ\theta — so the map covers the whole sphere, not just part of it.

Why Pauli expectation values match spherical coordinates. The identity x=⟨X⟩,y=⟨Y⟩,z=⟨Z⟩x=\langle X\rangle, y=\langle Y\rangle, z=\langle Z\rangle is not a coincidence; it is a direct computation. For ∣ψ⟩=cos⁡(θ/2)∣0⟩+eiφsin⁡(θ/2)∣1⟩|\psi\rangle=\cos(\theta/2)|0\rangle+e^{i\varphi}\sin(\theta/2)|1\rangle,

⟨Z⟩=cos⁡2(θ/2)−sin⁡2(θ/2)=cos⁡θ,\langle Z\rangle = \cos^2(\theta/2) - \sin^2(\theta/2) = \cos\theta,

using Z∣0⟩=∣0⟩,Z∣1⟩=−∣1⟩Z|0\rangle=|0\rangle, Z|1\rangle=-|1\rangle, and similarly expanding X=(0110)X=\begin{pmatrix}0&1\\1&0\end{pmatrix} and Y=(0−ii0)Y=\begin{pmatrix}0&-i\\i&0\end{pmatrix} against the same state gives ⟨X⟩=sin⁡θcos⁡φ\langle X\rangle=\sin\theta\cos\varphi and ⟨Y⟩=sin⁡θsin⁡φ\langle Y\rangle=\sin\theta\sin\varphi. This is exactly why a unitary gate, which conjugates the density matrix and therefore linearly transforms (⟨X⟩,⟨Y⟩,⟨Z⟩)(\langle X\rangle,\langle Y\rangle,\langle Z\rangle), acts as a rotation of the Bloch vector — the geometric intuition used throughout single-qubit gate theory is not a pedagogical shortcut, it is a provable consequence of the Pauli algebra.

Honest limitations. The Bloch sphere is a surface, and it represents only pure states — states with certainty of preparation, no classical uncertainty mixed in. A statistical mixture of pure states (a mixed state, described by a density matrix rather than a state vector) corresponds to a point strictly inside the sphere, with the fully mixed state (maximal uncertainty) at the exact center. The representation also does not extend simply to multiple qubits: an nn-qubit pure state has 2n+1−22^{n+1}-2 real parameters (after normalization and global phase), which for n≥2n\ge2 already exceeds what any single low-dimensional sphere can encode, and entangled states have no faithful "per-qubit Bloch vector" picture at all — tracing out the rest of the system to get a single-qubit Bloch vector for one qubit of an entangled pair necessarily loses information and lands strictly inside the sphere.

Worked example

Take the state produced by applying a Hadamard-like superposition with a phase, a common intermediate state in phase-estimation-style circuits:

∣ψ⟩=12∣0⟩+12i ∣1⟩=12∣0⟩+12eiπ/2∣1⟩.|\psi\rangle = \frac{1}{\sqrt2}|0\rangle + \frac{1}{\sqrt2} i\,|1\rangle = \frac{1}{\sqrt2}|0\rangle + \frac{1}{\sqrt2}e^{i\pi/2}|1\rangle.
StepComputationResult
Amplitudesα=12\alpha = \tfrac{1}{\sqrt2}, β=i2\beta = \tfrac{i}{\sqrt2}already normalized: ∣α∣2+∣β∣2=12+12=1|\alpha|^2+|\beta|^2 = \tfrac12+\tfrac12=1
Strip global phasearg⁡(α)=0\arg(\alpha)=0, so no rotation neededα′=12\alpha'=\tfrac{1}{\sqrt2} (already real, non-negative)
Polar angleθ=2arccos⁡(∣α∣)=2arccos⁡(1/2)\theta = 2\arccos(|\alpha|) = 2\arccos(1/\sqrt2)θ=2⋅π/4=π/2\theta = 2\cdot \pi/4 = \pi/2
Azimuthal angleφ=arg⁡(β)−arg⁡(α)=π/2−0\varphi = \arg(\beta)-\arg(\alpha) = \pi/2 - 0φ=π/2\varphi = \pi/2
Bloch vectorx=sin⁡θcos⁡φ, y=sin⁡θsin⁡φ, z=cos⁡θx=\sin\theta\cos\varphi,\ y=\sin\theta\sin\varphi,\ z=\cos\thetax=1⋅0=0x = 1\cdot 0 = 0,  y=1⋅1=1\ y = 1\cdot 1 = 1,  z=0\ z=0
InterpretationPoint on the equator, on the +y+y axisThis is exactly the state ∣+i⟩|+i\rangle

Sanity check going backward: plug (θ,φ)=(π/2,π/2)(\theta,\varphi)=(\pi/2,\pi/2) into `BlochAnglesToState`: α=cos⁡(π/4)=1/2\alpha=\cos(\pi/4)=1/\sqrt2, β=eiπ/2sin⁡(π/4)=i/2\beta=e^{i\pi/2}\sin(\pi/4) = i/\sqrt2 — matches the original state exactly, confirming the round trip. Now contrast with a state that has both a nontrivial global phase and a nontrivial relative phase, to show the two effects don't get confused:

∣ϕ⟩=eiπ/3(32∣0⟩+12eiπ/4∣1⟩).|\phi\rangle = e^{i\pi/3}\left(\tfrac{\sqrt3}{2}|0\rangle + \tfrac12 e^{i\pi/4}|1\rangle\right).

Here α=32eiπ/3\alpha = \tfrac{\sqrt3}{2}e^{i\pi/3} and β=12eiπ/3eiπ/4=12ei7π/12\beta = \tfrac12 e^{i\pi/3}e^{i\pi/4} = \tfrac12 e^{i7\pi/12}.

  • ∣α∣=3/2≈0.866|\alpha| = \sqrt3/2 \approx 0.866, so θ=2arccos⁡(0.866)=2⋅(π/6)=π/3\theta = 2\arccos(0.866) = 2\cdot(\pi/6) = \pi/3.
  • arg⁡(α)=π/3\arg(\alpha) = \pi/3, arg⁡(β)=7π/12\arg(\beta)=7\pi/12, so φ=7π/12−π/3=7π/12−4π/12=π/4\varphi = 7\pi/12 - \pi/3 = 7\pi/12 - 4\pi/12 = \pi/4.

Notice: the leading factor eiπ/3e^{i\pi/3} (global phase) never appears in θ\theta or φ\varphi — it cancels exactly in the subtraction arg⁡(β)−arg⁡(α)\arg(\beta)-\arg(\alpha). Only the relative phase π/4\pi/4 between the two terms survives, exactly as the theory predicts. Both ∣ϕ⟩|\phi\rangle and 32∣0⟩+12eiπ/4∣1⟩\tfrac{\sqrt3}{2}|0\rangle+\tfrac12 e^{i\pi/4}|1\rangle map to the identical point (θ,φ)=(π/3,π/4)(\theta,\varphi)=(\pi/3,\pi/4) on the sphere, confirming they are the same physical state.

Common pitfalls & traps

PitfallWhy it happensFix
Using θ\theta (not θ/2\theta/2) in the amplitudesConfusing the sphere angle with the raw trigonometric angle between amplitudesAlways use half-angles: cos⁡(θ/2),sin⁡(θ/2)\cos(\theta/2), \sin(\theta/2), so orthogonal states land at antipodal points
Treating global phase as observableForgetting that eiγ∣ψ⟩e^{i\gamma}|\psi\rangle gives identical measurement statisticsAlways strip global phase (make α\alpha real, non-negative) before reading off φ\varphi
Confusing relative phase with global phaseBoth look like "an ei⋅e^{i\cdot} somewhere in the formula"Relative phase is a difference of two arguments; global phase is a common factor on both amplitudes
Assuming mixed states have a Bloch sphere pointThe Bloch ball (interior) is often glossed over in intro treatmentsRemember: pure states are on the surface (∣r⃗∣=1|\vec r|=1); mixed states are strictly inside (∣r⃗∣<1|\vec r|<1)
Trying to extend the picture directly to 2+ qubitsThe sphere trick feels like it should generalizeIt does not: entangled multi-qubit states have no single per-qubit Bloch vector that captures all information
Azimuth undefined at poles causing NaN/crashes in code`arg(0)` is undefined; dividing zero vectorsSpecial-case θ=0\theta=0 or θ=π\theta=\pi and return a fixed convention for φ\varphi

Pioneers & history

The Bloch sphere is named for Felix Bloch, who introduced the geometric picture in the context of nuclear magnetic resonance, where a spin-1/2 particle's state precesses under a magnetic field in a way that is naturally visualized as a vector on a sphere. The same geometry was later recognized as exactly the right tool for describing any two-level quantum system, including the qubit of quantum computing.

Origins · 1 minFelix Bloch and how this was discovered
FB
Felix Bloch
1905–1983 · Developed the mathematical description of nuclear induction and spin precession, for which he shared the 1952 Nobel Prize in Physics; the sphere used to visualize spin states bears his name.

history

  • 1946 — Felix Bloch publishes the equations describing nuclear magnetic resonance, introducing the vector picture of spin precession that becomes the Bloch sphere.
  • 1957 — The Bloch sphere is adopted broadly in the physics of two-level systems (masers, lasers) beyond its original NMR context.
  • 1980s — As quantum information theory develops, the two-level system is formalized as the abstract "qubit," and the Bloch sphere becomes its standard visualization.
  • 1990s — Textbooks and courses on quantum computing (following the rise of Dirac Notation and the State Vector as a teaching foundation) adopt the Bloch sphere as the canonical picture for single-qubit gates.
  • 2000s–present — The Bloch sphere becomes a standard diagnostic tool in experimental quantum hardware labs, used to visualize measured qubit states and gate errors on real superconducting and trapped-ion devices.

Where this matters

Every major single-qubit gate — XX, YY, ZZ, Hadamard, phase gates, and general rotation gates — is most easily understood and remembered as a specific rotation of the Bloch sphere: XX is a π\pi rotation about the xx-axis, ZZ a π\pi rotation about the zz-axis, and a general rotation gate Rn^(θ)R_{\hat n}(\theta) rotates the Bloch vector by angle θ\theta about axis n^\hat n. This geometric fact is exactly what makes the transition to Single-Qubit Gates: X, Y, Z, H, Rotation Gates and Arbitrary Single-Qubit Unitaries intuitive rather than purely algebraic. In experimental physics, quantum hardware teams (superconducting circuits, trapped ions, photonic qubits) routinely plot measured qubit states on Bloch spheres to visualize calibration errors, decoherence (the Bloch vector shrinking toward the center as the state becomes mixed over time), and the fidelity of applied pulses. Quantum optimal control — designing microwave pulses that implement a target gate — is frequently framed as steering a point along a desired path on the sphere. Software toolkits for quantum computing (simulators and visualizers used alongside frameworks like Qiskit, Cirq, and similar platforms) include Bloch sphere plotting as a core diagnostic, because it is the fastest way for a human to sanity-check a single-qubit state or debug an unexpected phase. The concept also underlies how textbooks introduce the geometric meaning of The Qubit and Superposition, giving beginners a concrete picture before the abstract linear-algebra machinery of multi-qubit systems and entanglement is introduced.

Ask the tutor: How do the X, Y, and Z gates each correspond to a specific rotation axis on the Bloch sphere, and why does that fully explain their action on any input state?

Why does normalization plus removing global phase leave exactly two real parameters?

A general amplitude pair (α,β)∈C2(\alpha,\beta) \in \mathbb{C}^2 has four real degrees of freedom (real and imaginary parts of each). The normalization constraint ∣α∣2+∣β∣2=1|\alpha|^2+|\beta|^2=1 is one real equation, cutting this to three — geometrically, three real parameters describe a point on the unit sphere in four-dimensional real space (S3S^3). Global phase equivalence, ∣ψ⟩∼eiγ∣ψ⟩|\psi\rangle \sim e^{i\gamma}|\psi\rangle, identifies an entire circle's worth of points as physically identical, removing one more real parameter. Three minus one leaves two, matching the two angles θ\theta and φ\varphi of ordinary spherical coordinates. This dimension-counting argument is why the Bloch sphere, an ordinary 2-dimensional surface, is exactly the right shape — not a coincidence of the specific parametrization chosen.

Can I recover the original amplitudes exactly from Bloch angles, including phase?

You can recover the amplitudes up to global phase, which is all that is physically meaningful. Given (θ,φ)(\theta,\varphi), the formula ∣ψ⟩=cos⁡(θ/2)∣0⟩+eiφsin⁡(θ/2)∣1⟩|\psi\rangle=\cos(\theta/2)|0\rangle+e^{i\varphi}\sin(\theta/2)|1\rangle produces one specific representative state, with the ∣0⟩|0\rangle coefficient real and non-negative. Any other state that differs from this by an overall phase factor eiγe^{i\gamma} is physically the same state and maps to the same (θ,φ)(\theta,\varphi). So the reconstruction is exact for the equivalence class of physically indistinguishable states, but it deliberately discards information (the specific global phase) that no experiment could ever detect anyway.

Why is the half-angle θ/2 used instead of just using the raw angle between the amplitudes?

Two orthogonal states, meaning states that are perfectly distinguishable by measurement, should intuitively sit as far apart as possible on the sphere — at opposite poles. If the amplitudes were parametrized with the full angle θ\theta directly (no halving), then ∣0⟩|0\rangle and ∣1⟩|1\rangle would end up 90 degrees apart on the sphere rather than 180 degrees apart, which would clash with the physical intuition that orthogonal states are "maximally different." The halving in cos⁡(θ/2),sin⁡(θ/2)\cos(\theta/2), \sin(\theta/2) is precisely the convention that fixes this, making perfectly distinguishable states antipodal and, more generally, making the geometric angle between two Bloch vectors relate cleanly to their quantum mechanical overlap.

Does every point inside the sphere correspond to some real physical situation, or is it just decoration?

It is fully physical. Points strictly inside the sphere correspond to mixed states, described by a density matrix ρ=12(I+r⃗⋅σ⃗)\rho = \tfrac12(I + \vec r\cdot\vec\sigma) where r⃗\vec r is the Bloch vector with ∣r⃗∣<1|\vec r|<1 and σ⃗=(X,Y,Z)\vec\sigma=(X,Y,Z). Such states arise naturally from decoherence (a qubit interacting with an uncontrolled environment), from classical uncertainty about which pure state was prepared, or from looking at one qubit of an entangled pair while ignoring the other. The exact center, r⃗=0\vec r = 0, is the maximally mixed state, representing complete ignorance — equivalent to a fair coin flip between ∣0⟩|0\rangle and ∣1⟩|1\rangle with no phase information at all.

Why can't this same trick be used to visualize two-qubit or n-qubit states?

Dimension counting again explains the limitation. An nn-qubit pure state lives in C2n\mathbb{C}^{2^n}, so after normalization and removing global phase it has 2(2n)−22(2^n)-2 real parameters — for n=2n=2 that is already 6 real parameters, far more than the 2 that a single sphere's surface can hold. Worse, entanglement means the two qubits' individual states are not even independently well-defined pure states; tracing out one qubit to look at the other generally produces a mixed state (a point inside its own local Bloch ball), and no assignment of two separate Bloch vectors recovers the correlations between the qubits. Genuinely faithful pictures for multi-qubit entangled states require higher-dimensional or more abstract tools.

How does the Bloch sphere make it easy to understand what a quantum gate does?

Because single-qubit unitary gates act on the Bloch vector exactly as three-dimensional rotations (and, if applied to a mixed state, contractions are absent for unitaries — only rotations occur, since unitary evolution preserves ∣r⃗∣|\vec r| for pure states and in fact preserves it for all states). This follows from the fact that conjugating the density matrix by a unitary is a linear, norm-preserving, orientation-preserving map on the real 3-vector (⟨X⟩,⟨Y⟩,⟨Z⟩)(\langle X\rangle,\langle Y\rangle,\langle Z\rangle), which any linear algebra course shows must be a rotation matrix. So instead of multiplying 2×22\times2 complex matrices to see what a gate does, you can just picture spinning the sphere about an axis by some angle — the intuition later used to describe XX, YY, ZZ, Hadamard, and general rotation gates geometrically.

What happens to the Bloch angles if I measure the qubit?

Measurement in the computational basis projects the state onto either ∣0⟩|0\rangle (the north pole) or ∣1⟩|1\rangle (the south pole), with probabilities cos⁡2(θ/2)\cos^2(\theta/2) and sin⁡2(θ/2)\sin^2(\theta/2) respectively — read directly off the polar angle. This is the Born rule expressed geometrically: the closer the Bloch point is to the north pole, the more likely a ∣0⟩|0\rangle outcome. After the measurement, the qubit's state vector collapses to whichever pole was observed, discarding all information about φ\varphi; the azimuthal angle, which encoded the relative phase, becomes physically irrelevant once the state collapses to a basis state, since a pole has no well-defined longitude.

Before this

Where this leads