Single-Qubit Rotation Gates

Every operation a quantum computer performs on a single qubit is a 2×22\times2 unitary, and once the physically invisible global phase is stripped away it is an element of SU(2)SU(2) — the group of unit-determinant 2×22\times2 unitaries. That group is generated by exactly three matrices, the Pauli operators X,Y,ZX,Y,Z, so every single-qubit gate is the exponential of a real combination of Paulis. The rotation gates Rx(θ),Ry(θ),Rz(θ)R_x(\theta),R_y(\theta),R_z(\theta) are the cleanest slice of this fact: each is the exponential of a single Pauli about a single Cartesian axis, and together with rotations about tilted axes they sweep out all of SU(2)SU(2). This article is where the discrete Pauli gates of the previous nodes become a continuous, angle-tunable family, and where the Bloch sphere stops being a static picture and starts to spin. Two threads run through everything below. First, an algebraic gift: because n^⋅σ⃗\hat n\cdot\vec\sigma squares to the identity, its exponential resums into the finite closed form cos⁡(θ/2) I−isin⁡(θ/2) (n^⋅σ⃗)\cos(\theta/2)\,I-i\sin(\theta/2)\,(\hat n\cdot\vec\sigma) — the Pauli plays the imaginary unit, and Euler's formula does the rest. Second, a geometric one: conjugating a state by Rn^(θ)R_{\hat n}(\theta) rotates its Bloch vector rigidly by the full angle θ\theta about n^\hat n, though the gate itself carries only the half-angle θ/2\theta/2. That mismatch is the fingerprint of the two-to-one map SU(2)→SO(3)SU(2)\to SO(3) — the reason a qubit is a spinor, and the reason a 2π2\pi turn flips a state's sign while a 4π4\pi turn is needed to truly do nothing. From here the Pauli gates reappear as half-turns, the Hadamard and phase gates as short rotation products, and the Z–Y decomposition delivers single-qubit universality.

Key takeaways
  • A rotation gate is the exponential of a Pauli about an axis: Rn^(θ)=e−iθ n^⋅σ⃗/2R_{\hat n}(\theta)=e^{-i\theta\,\hat n\cdot\vec\sigma/2}, a continuous, angle-tunable single-qubit unitary.
  • Because (n^⋅σ⃗)2=I(\hat n\cdot\vec\sigma)^2=I, the exponential collapses to the closed form cos⁡(θ/2) I−isin⁡(θ/2) (n^⋅σ⃗)\cos(\theta/2)\,I-i\sin(\theta/2)\,(\hat n\cdot\vec\sigma) — no infinite series survives.
  • On the Bloch sphere the gate is an honest rigid rotation by angle θ\theta about n^\hat n; the gate carries θ/2\theta/2, the sphere turns θ\theta, and that factor of two is the SU(2)→SO(3)SU(2)\to SO(3) double cover.
  • Special angles recover the discrete gates: Rx(π)=−iXR_x(\pi)=-iX, Rz(π/2)=e−iπ/4SR_z(\pi/2)=e^{-i\pi/4}S, and Rn^(2π)=−IR_{\hat n}(2\pi)=-I (not II).
  • Any single-qubit unitary factors as U=eiαRz(β)Ry(γ)Rz(δ)U=e^{i\alpha}R_z(\beta)R_y(\gamma)R_z(\delta) (Z–Y decomposition), so three rotations plus a phase realise all of SU(2)SU(2) — the root of single-qubit universality.
  • Rotations are the native, calibrated primitives of real hardware: virtual-ZZ frame changes, resonant Rx(π/2)R_x(\pi/2) pulses, and the parameter-shift gradients that train variational circuits.
  • Key idea: each Pauli squares to identity, so its exponential resums to Rn^(θ)=cos⁡(θ/2) I−isin⁡(θ/2) (n^⋅σ⃗)R_{\hat n}(\theta)=\cos(\theta/2)\,I-i\sin(\theta/2)\,(\hat n\cdot\vec\sigma) — a half-angle SU(2)SU(2) gate that drives a full-angle SO(3)SO(3) rotation of the Bloch vector.
The core insight

Each Pauli squares to identity, so its exponential resums into a finite Euler formula, Rn^(θ)=cos⁡(θ/2) I−isin⁡(θ/2) (n^⋅σ⃗)R_{\hat n}(\theta)=\cos(\theta/2)\,I-i\sin(\theta/2)\,(\hat n\cdot\vec\sigma) — a half-angle rotation in SU(2)SU(2) that rotates the Bloch vector by the full angle θ\theta about n^\hat n in SO(3)SO(3).

Where this lives today

Rotation gates are not a pedagogical abstraction; they are the literal instruction set of every quantum processor now running. On superconducting hardware the ZZ-rotation is a virtual-ZZ gate implemented for free — the control software advances the phase reference of all subsequent microwave pulses, so Rz(θ)R_z(\theta) takes zero time and zero error (McKay and co-workers, 2017); the only driven single-qubit gate is then a calibrated Rx(π/2)R_x(\pi/2) (the "X\sqrt{X}"), from which any rotation compiles as RzRxRzR_z R_x R_z. Trapped-ion machines expose driven Rx(θ),Ry(θ)R_x(\theta),R_y(\theta) directly, the axis set by the drive phase and the angle by the pulse area Ωt\Omega t. Above the hardware, variational algorithms (VQE, QAOA) build ansätze from parametrised rotations and train them by the parameter-shift rule, an exact gradient that exploits the generator's two-eigenvalue structure. And at the fault-tolerant frontier the Solovay–Kitaev theorem and number-theoretic compilers (Ross–Selinger) approximate any Rz(θ)R_z(\theta) to accuracy ε\varepsilon with O(log⁡(1/ε))O(\log(1/\varepsilon)) Clifford+TT gates — the bridge from continuous angles to the discrete, error-corrected circuits a scalable machine will run.

The physical picture

Hold the operator n^⋅σ⃗\hat n\cdot\vec\sigma in your hand as if it were the number ii. It is Hermitian, traceless, and — the crucial fact — it squares to the identity, exactly the property i2=−1i^2=-1 upgraded to a matrix. Feed it into an exponential and you do not get an unruly infinite series; you get sines and cosines, because the even powers pile up on II and the odd powers pile up on n^⋅σ⃗\hat n\cdot\vec\sigma, just as eiϕe^{i\phi} splits into cos⁡ϕ\cos\phi and isin⁡ϕi\sin\phi. A rotation gate is therefore a smooth, sinusoidal dial between doing nothing (the II term at θ=0\theta=0) and applying the full Pauli reflection (the n^⋅σ⃗\hat n\cdot\vec\sigma term at θ=π\theta=\pi). Geometrically that dial is a turntable for the Bloch sphere: choose an axis, choose an angle, and spin the whole sphere. The first figure makes the half-angle-in, full-angle-out relationship concrete for one specific gate.

Figure
Bloch-sphere turntable: the gate R_x(π/2) = cos(π/4)I − i sin(π/4)X = 0.7071(I − iX) rotates the state's Bloch vector rigidly by the full angle θ = 90° about x̂, carrying the north pole (the 0 state) to the −ŷ point, while the gate itself carries only the half-angle θ/2 = 45° — the fingerprint of the SU(2)→SO(3) double cover
Gate carries θ/2 = 45°, the Bloch sphere turns by the full θ = 90° z, |0⟩ |1⟩ +ŷ −ŷ x̂ = rotation axis n̂ θ = 90° r₀ = +ẑ r₁ = −ŷ R_x(π/2) = cos(π/4)I − i sin(π/4)X = 0.7071 (I − iX) gate phase θ/2 = 45° sphere angle θ = 90°

Push the dial further and the double cover announces itself. As θ\theta runs through π\pi, 2π2\pi, and beyond, the matrix entries cos⁡(θ/2)\cos(\theta/2) and sin⁡(θ/2)\sin(\theta/2) trace sinusoids of period 4π4\pi, not 2π2\pi. A full 2π2\pi turn brings the Bloch vector exactly home, yet the state vector has quietly picked up a global −1-1: the gate is −I-I, not II. Only after a second full turn, at θ=4π\theta=4\pi, does it return to the true identity. This holds for every axis n^\hat n, and it is the same −1-1 a neutron beam accumulates when its spin is precessed through 2π2\pi. The second figure plots both entries across two full Bloch turns so the 4π4\pi periodicity is impossible to miss.

Figure
The double cover made visible: the matrix entries of R_n(θ) are cos(θ/2) and sin(θ/2), so the gate has period 4π rather than 2π — at θ = 2π the gate equals −I, a full Bloch turn that flips the state's global sign, and only at θ = 4π does it return to the identity I
The punchline: R_n(θ) has period 4π — a 2π turn gives −I, not I +1 0 −1 0 π 2π 3π 4π rotation angle θ cos(θ/2) sin(θ/2) θ = 2π: R = −I θ = 4π: R = I θ = π: R = −i n·σ

The core equations

The defining object is the exponential of a Pauli about an axis. A real unit vector n^\hat n picks the axis, the angle θ\theta sets the amount of turn, and the generator is −i/2-i/2 times the corresponding component of the Pauli vector — the −i-i makes the result unitary, and the 12\tfrac12 is the half-angle that the double cover forces on us.

Rn^(θ)=exp⁡ ⁣(−i θ2 n^⋅σ⃗)(1)R_{\hat n}(\theta)=\exp\!\left(-i\,\frac{\theta}{2}\,\hat n\cdot\vec\sigma\right) \qquad (1)

where: θ\theta — rotation angle on the Bloch sphere (radians); n^=(nx,ny,nz)\hat n=(n_x,n_y,n_z) — real unit rotation axis; σ⃗=(X,Y,Z)\vec\sigma=(X,Y,Z) — the Pauli vector; n^⋅σ⃗=nxX+nyY+nzZ\hat n\cdot\vec\sigma=n_xX+n_yY+n_zZ — a traceless Hermitian involution. Reading: the generator n^⋅σ⃗\hat n\cdot\vec\sigma lives in the Lie algebra su(2)\mathfrak{su}(2); exponentiating it sweeps out the group SU(2)SU(2). The single most useful identity in the subject is that this exponential is finite. Because the generator squares to the identity, the Taylor series folds up exactly into a cosine on II and a sine on the generator — Euler's formula with a matrix in the seat of ii.

Rn^(θ)=cos⁡ ⁣θ2 I−isin⁡ ⁣θ2 (n^⋅σ⃗)(2)R_{\hat n}(\theta)=\cos\!\frac{\theta}{2}\,I-i\sin\!\frac{\theta}{2}\,(\hat n\cdot\vec\sigma) \qquad (2)

where: II — the 2×22\times2 identity. Reading: every single-qubit rotation is one cosine, one sine, and the axis — nothing else to compute, and the derivation below shows why it ends. Specialising the axis to x^\hat x, y^\hat y, z^\hat z gives the three canonical rotation matrices. Written out, RxR_x and RzR_z are complex while RyR_y is the lone all-real member — a fact that makes RyR_y the workhorse of amplitude-only rotations.

Rx(θ)=(cos⁡θ2−isin⁡θ2−isin⁡θ2cos⁡θ2)(3)R_x(\theta)=\begin{pmatrix}\cos\frac{\theta}{2}&-i\sin\frac{\theta}{2}\\[2pt]-i\sin\frac{\theta}{2}&\cos\frac{\theta}{2}\end{pmatrix} \qquad (3)
Ry(θ)=(cos⁡θ2−sin⁡θ2sin⁡θ2cos⁡θ2)(4)R_y(\theta)=\begin{pmatrix}\cos\frac{\theta}{2}&-\sin\frac{\theta}{2}\\[2pt]\sin\frac{\theta}{2}&\cos\frac{\theta}{2}\end{pmatrix} \qquad (4)
Rz(θ)=(e−iθ/200eiθ/2)=e−iθ/2 diag ⁣(1,eiθ)(5)R_z(\theta)=\begin{pmatrix}e^{-i\theta/2}&0\\0&e^{i\theta/2}\end{pmatrix}=e^{-i\theta/2}\,\mathrm{diag}\!\left(1,e^{i\theta}\right) \qquad (5)

where: each θ\theta is the turn about the named axis. Reading: RzR_z is diagonal because ZZ is; the rightmost form of (5) exposes that Rz(θ)R_z(\theta) equals the phase gate diag(1,eiθ)\mathrm{diag}(1,e^{i\theta}) only up to the global factor e−iθ/2e^{-i\theta/2} — a distinction the Common-pitfall box below insists on. The geometric payoff is that conjugating the Bloch-vector operator by a rotation gate produces a genuine 3×33\times3 real rotation of the Bloch vector. The gate acts on 22-component amplitudes; its shadow on R3\mathbb R^3 is an ordinary rotation matrix.

Rn^(θ) (r⃗⋅σ⃗) Rn^(θ)†=(Rn^(θ) r⃗)⋅σ⃗(6)R_{\hat n}(\theta)\,(\vec r\cdot\vec\sigma)\,R_{\hat n}(\theta)^\dagger=\big(\mathcal R_{\hat n}(\theta)\,\vec r\big)\cdot\vec\sigma \qquad (6)

where: r⃗\vec r — Bloch vector, ρ=12(I+r⃗⋅σ⃗)\rho=\tfrac12(I+\vec r\cdot\vec\sigma); Rn^(θ)∈SO(3)\mathcal R_{\hat n}(\theta)\in SO(3) — the 3×33\times3 rotation by angle θ\theta about n^\hat n. Reading: the map Rn^(θ)↦Rn^(θ)R_{\hat n}(\theta)\mapsto\mathcal R_{\hat n}(\theta) is the two-to-one homomorphism SU(2)→SO(3)SU(2)\to SO(3); both ±Rn^(θ)\pm R_{\hat n}(\theta) map to the same Rn^(θ)\mathcal R_{\hat n}(\theta). Two rotations about different axes compose into a single rotation about a third — the quaternion multiplication rule in disguise, and the algebraic core of gate synthesis.

Ra^(θ)Rb^(ϕ)=[cθcϕ−sθsϕ (a^ ⁣⋅ ⁣b^)]I−i[cθsϕ b^+sθcϕ a^+sθsϕ (a^ ⁣× ⁣b^)] ⁣⋅ ⁣σ⃗(7)R_{\hat a}(\theta)R_{\hat b}(\phi)=\big[c_\theta c_\phi-s_\theta s_\phi\,(\hat a\!\cdot\!\hat b)\big]I-i\big[c_\theta s_\phi\,\hat b+s_\theta c_\phi\,\hat a+s_\theta s_\phi\,(\hat a\!\times\!\hat b)\big]\!\cdot\!\vec\sigma \qquad (7)

where: cθ=cos⁡(θ/2), sθ=sin⁡(θ/2)c_\theta=\cos(\theta/2),\ s_\theta=\sin(\theta/2) (likewise for ϕ\phi); a^,b^\hat a,\hat b — the two axes. Reading: the scalar part gives the resultant half-angle via cos⁡(Θ/2)=cθcϕ−sθsϕ(a^⋅b^)\cos(\Theta/2)=c_\theta c_\phi-s_\theta s_\phi(\hat a\cdot\hat b); the vector part gives the resultant axis. Finally, three rotations exhaust the group: any single-qubit unitary is a global phase times a ZZ–YY–ZZ sandwich. This is the Z–Y decomposition, and it is the reason the rotation family alone is universal for one qubit.

U=eiα Rz(β) Ry(γ) Rz(δ)(8)U=e^{i\alpha}\,R_z(\beta)\,R_y(\gamma)\,R_z(\delta) \qquad (8)

where: α\alpha — global phase in U(1)U(1); β,γ,δ\beta,\gamma,\delta — three real Euler angles fixing the SU(2)SU(2) part. Reading: U(2)U(2) has four real parameters; one is the phase α\alpha and three are the Euler angles, matching the count exactly.

Derivation AResumming the exponential into a closed-form rotation

The goal is to prove that the matrix exponential in equation (1) is not an infinite object at all but the two-term closed form of equation (2), and to see precisely which property of the generator makes the collapse happen. The lever is a single algebraic fact — that n^⋅σ⃗\hat n\cdot\vec\sigma is an involution, squaring to the identity — which lets us sort the power series into even and odd terms and recognise the cosine and sine series hiding inside. This is the exact analogue of deriving Euler's formula eiϕ=cos⁡ϕ+isin⁡ϕe^{i\phi}=\cos\phi+i\sin\phi from i2=−1i^2=-1, with the Pauli combination playing the role of ii.

Step 1. Abbreviate A=n^⋅σ⃗A=\hat n\cdot\vec\sigma and expand the exponential as a power series, e−i(θ/2)A=∑k≥01k!(−iθ2)kAke^{-i(\theta/2)A}=\sum_{k\ge0}\frac{1}{k!}\left(-i\tfrac{\theta}{2}\right)^k A^k. Why: nothing is assumed yet — this is just the definition of a matrix exponential, valid because the series converges absolutely for any finite matrix.

Step 2. Show A2=IA^2=I. Squaring, A2=∑j,knjnk σjσkA^2=\sum_{j,k}n_jn_k\,\sigma_j\sigma_k; the symmetric part uses σjσj=I\sigma_j\sigma_j=I and the antisymmetric part cancels because distinct Paulis anticommute, σjσk+σkσj=0\sigma_j\sigma_k+\sigma_k\sigma_j=0. What survives is (nx2+ny2+nz2)I=I(n_x^2+n_y^2+n_z^2)I=I since n^\hat n is a unit vector. Why: this is the one fact that makes the whole series terminate — every power of AA is now either II or AA.

Step 3. Split by parity. Even powers give A2m=IA^{2m}=I and odd powers give A2m+1=AA^{2m+1}=A, so the sum separates into (∑m(−1)m(θ/2)2m(2m)!)I−i(∑m(−1)m(θ/2)2m+1(2m+1)!)A\big(\sum_m\frac{(-1)^m(\theta/2)^{2m}}{(2m)!}\big)I-i\big(\sum_m\frac{(-1)^m(\theta/2)^{2m+1}}{(2m+1)!}\big)A. Why: gathering the II terms and the AA terms is legitimate because they are linearly independent, so no cross-contamination occurs.

Step 4. Recognise the two bracketed series as cos⁡(θ/2)\cos(\theta/2) and sin⁡(θ/2)\sin(\theta/2), giving Rn^(θ)=cos⁡(θ/2) I−isin⁡(θ/2) AR_{\hat n}(\theta)=\cos(\theta/2)\,I-i\sin(\theta/2)\,A. Why: this is equation (2); the exponential of any axis-Pauli is a finite, hand-computable 2×22\times2 matrix — the practical foundation for everything that follows. *As a concrete check: Rx(π/2)=0.7071(I−iX)R_x(\pi/2)=0.7071(I-iX) and Rz(π)=−iZR_z(\pi)=-iZ reproduce the closed form exactly, and the involution A2=IA^2=I holds for a tilted axis such as n^=(1,1,−1)/3\hat n=(1,1,-1)/\sqrt3.*

Go deeper · Derivation B · 1 minConjugation turns the gate into a rigid Bloch-sphere rotation
Derivation BConjugation turns the gate into a rigid Bloch-sphere rotation

The goal is to demonstrate equation (6) concretely for Rz(θ)R_z(\theta) and, in doing so, to expose the factor of two between the gate's argument and the geometric turn. We compute how Rz(θ)R_z(\theta) conjugates the operator XX; the result is a linear combination of XX and YY with coefficients cos⁡θ\cos\theta and sin⁡θ\sin\theta — the full angle, not the half-angle — which is exactly a rotation of the xx–yy plane of the Bloch sphere about the zz-axis. This is where the half-angle in the gate and the full angle on the sphere are reconciled, and it corrects a sign that is easy to drop.

Step 1. Write the conjugation in generator form, Rz(θ)XRz(θ)†=e−iθZ/2 X eiθZ/2R_z(\theta)XR_z(\theta)^\dagger=e^{-i\theta Z/2}\,X\,e^{i\theta Z/2}, and apply the Hadamard lemma eGXe−G=X+[G,X]+12[G,[G,X]]+⋯e^{G}Xe^{-G}=X+[G,X]+\tfrac12[G,[G,X]]+\cdots with G=−iθZ/2G=-i\theta Z/2. Why: the nested-commutator expansion turns operator conjugation into a series in θ\theta that we can resum, mirroring Step 1 of the previous derivation.

Step 2. Evaluate the first commutator using [Z,X]=2iY[Z,X]=2iY, so [G,X]=−iθ2(2iY)=+θ Y[G,X]=-\tfrac{i\theta}{2}(2iY)=+\theta\,Y. The double commutator uses [Z,Y]=−2iX[Z,Y]=-2iX, giving [G,[G,X]]=−θ2X[G,[G,X]]=-\theta^2X, and the pattern then cycles X→Y→−X→−YX\to Y\to -X\to -Y. Why: the commutators reproduce the derivative structure of sine and cosine, which is the algebraic reason a rotation appears.

Step 3. Resum the two interleaved series to Rz(θ)XRz(θ)†=cos⁡θ X+sin⁡θ YR_z(\theta)XR_z(\theta)^\dagger=\cos\theta\,X+\sin\theta\,Y. A direct matrix multiplication confirms the sign: Rz(θ)XRz(θ)†R_z(\theta)XR_z(\theta)^\dagger has off-diagonal entries e∓iθe^{\mp i\theta}, which equal cos⁡θ X+sin⁡θ Y\cos\theta\,X+\sin\theta\,Y, not the minus-sign variant. Why: keeping this sign correct is what makes the operator identity consistent with the state-vector rotation in the next step — a subtlety the thin migrated draft got backwards.

Step 4. Translate to the Bloch vector. Since ri′=Tr(ρ R†σiR)r_i'=\mathrm{Tr}(\rho\,R^\dagger\sigma_iR), the transformed components are rx′=rxcos⁡θ−rysin⁡θr_x'=r_x\cos\theta-r_y\sin\theta, ry′=rxsin⁡θ+rycos⁡θr_y'=r_x\sin\theta+r_y\cos\theta, rz′=rzr_z'=r_z — the standard SO(3)SO(3) rotation by the full angle θ\theta about z^\hat z. Why: the gate argument is θ/2\theta/2 but the sphere turns by θ\theta, precisely the two-to-one SU(2)→SO(3)SU(2)\to SO(3) relationship, verified against the state Rz(θ)∣+⟩R_z(\theta)|+\rangle whose Bloch vector is (cos⁡θ,sin⁡θ,0)(\cos\theta,\sin\theta,0). *As a concrete check: Rz(0.7)XRz(0.7)†=cos⁡(0.7)X+sin⁡(0.7)YR_z(0.7)XR_z(0.7)^\dagger=\cos(0.7)X+\sin(0.7)Y; the alternative −sin⁡θ Y-\sin\theta\,Y form gives the wrong sign and fails this check.*

Go deeper · Derivation C · 1 minEvery single-qubit gate is three rotations
Derivation CEvery single-qubit gate is three rotations

The goal is to prove the Z–Y decomposition of equation (8): that an arbitrary U∈U(2)U\in U(2) can be written as a global phase times Rz(β)Ry(γ)Rz(δ)R_z(\beta)R_y(\gamma)R_z(\delta), so the rotation gates alone generate every single-qubit operation. The strategy is to strip the global phase to reduce to SU(2)SU(2), note that any SU(2)SU(2) matrix has a rigid two-parameter shape, and then show that the explicit product Rz(β)Ry(γ)Rz(δ)R_z(\beta)R_y(\gamma)R_z(\delta) realises that shape for suitable angles. This is the concrete engine behind single-qubit universality and behind every compiler that turns an abstract gate into hardware rotations.

Step 1. Factor the phase. Any U∈U(2)U\in U(2) has ∣det⁡U∣=1|\det U|=1, so det⁡U=e2iα\det U=e^{2i\alpha} for some real α\alpha; then e−iαU∈SU(2)e^{-i\alpha}U\in SU(2). Every special-unitary 2×22\times2 matrix has the form (a−bˉbaˉ)\begin{pmatrix}a&-\bar b\\ b&\bar a\end{pmatrix} with ∣a∣2+∣b∣2=1|a|^2+|b|^2=1. Why: this two-complex-parameter shape, constrained by one real equation, carries exactly the three real degrees of freedom of SU(2)SU(2).

Step 2. Multiply out the ansatz. Using equations (3)–(5), Rz(β)Ry(γ)Rz(δ)R_z(\beta)R_y(\gamma)R_z(\delta) has top-left entry e−i(β+δ)/2cos⁡(γ/2)e^{-i(\beta+\delta)/2}\cos(\gamma/2) and bottom-left entry ei(β−δ)/2sin⁡(γ/2)e^{i(\beta-\delta)/2}\sin(\gamma/2), and the other two entries are fixed by the (a−bˉbaˉ)\begin{pmatrix}a&-\bar b\\ b&\bar a\end{pmatrix} pattern. Why: the product already lands in SU(2)SU(2), so we only need to hit the correct aa and bb.

Step 3. Match parameters. Set a=e−i(β+δ)/2cos⁡(γ/2)a=e^{-i(\beta+\delta)/2}\cos(\gamma/2) and b=ei(β−δ)/2sin⁡(γ/2)b=e^{i(\beta-\delta)/2}\sin(\gamma/2); then ∣a∣=∣cos⁡(γ/2)∣|a|=|\cos(\gamma/2)| and ∣b∣=∣sin⁡(γ/2)∣|b|=|\sin(\gamma/2)| fix γ\gamma, while the two phases fix β+δ\beta+\delta and β−δ\beta-\delta, hence β\beta and δ\delta separately. Why: three real equations for three angles always have a solution, so no SU(2)SU(2) target is unreachable.

Step 4. Reassemble. Restoring the phase gives U=eiαRz(β)Ry(γ)Rz(δ)U=e^{i\alpha}R_z(\beta)R_y(\gamma)R_z(\delta), the claimed decomposition. Why: the same argument with RyR_y replaced by RxR_x, or RzR_z by RxR_x, gives the Z–X and X–Y variants — the freedom compilers exploit to match a device's native axes. *As a concrete check: for (β,γ,δ)=(0.9,1.3,0.4)(\beta,\gamma,\delta)=(0.9,1.3,0.4) the assembled product's entries match aa and bb exactly.*

Limiting cases

Sending the angle or axis to special values recovers the discrete gates and the double-cover signature — each an entry the rotation family "contains" up to a global phase.

CaseResultMeaning
θ=π\theta=\pi about x^\hat xRx(π)=−iXR_x(\pi)=-iXa half-turn about x^\hat x is the Pauli XX gate up to the unobservable phase −i-i
θ=2π\theta=2\pi, any axisRn^(2π)=−IR_{\hat n}(2\pi)=-Ia full Bloch turn returns the vector but flips the state's sign — the spinor double cover
θ=π/2\theta=\pi/2 about z^\hat zRz(π/2)=e−iπ/4SR_z(\pi/2)=e^{-i\pi/4}Srecovers the phase gate SS up to a phase; likewise Rz(π/4)=e−iπ/8TR_z(\pi/4)=e^{-i\pi/8}T
θ→0\theta\to0Rn^(θ)≈I−iθ2 n^⋅σ⃗R_{\hat n}(\theta)\approx I-\tfrac{i\theta}{2}\,\hat n\cdot\vec\sigmainfinitesimal rotations are generated by the Pauli; the generator is −i2n^⋅σ⃗-\tfrac{i}{2}\hat n\cdot\vec\sigma

Worked examples

1 · Composing two π/2 rotations into a single axis-angle turn (fully worked)

A circuit applies Rx(π/2)R_x(\pi/2) and then Ry(π/2)R_y(\pi/2) to a qubit. Every product of rotations is again a rotation — find the single axis n^\hat n and angle Θ\Theta of the combined gate U=Ry(π/2)Rx(π/2)U=R_y(\pi/2)R_x(\pi/2).

Step 1. Write each factor in closed form. With cos⁡(π/4)=sin⁡(π/4)=1/2\cos(\pi/4)=\sin(\pi/4)=1/\sqrt2, both are Rx(π/2)=12(I−iX)R_x(\pi/2)=\tfrac{1}{\sqrt2}(I-iX) and Ry(π/2)=12(I−iY)R_y(\pi/2)=\tfrac{1}{\sqrt2}(I-iY).

Step 2. Multiply. U=12(I−iY)(I−iX)=12(I−iX−iY+i2 YX)U=\tfrac12(I-iY)(I-iX)=\tfrac12\big(I-iX-iY+i^2\,YX\big). Using YX=−iZYX=-iZ, the last term is i2(−iZ)=+iZi^2(-iZ)=+iZ, so U=12I−i2(X+Y−Z)U=\tfrac12 I-\tfrac{i}{2}(X+Y-Z).

Step 3. Read off the axis-angle by matching to cos⁡(Θ/2)I−isin⁡(Θ/2)(n^⋅σ⃗)\cos(\Theta/2)I-i\sin(\Theta/2)(\hat n\cdot\vec\sigma). The identity coefficient gives cos⁡(Θ/2)=12\cos(\Theta/2)=\tfrac12, so Θ/2=60∘\Theta/2=60^\circ and \boldsymbol{\Theta=120^\circ}. The Pauli coefficient is 12(X+Y−Z)\tfrac12(X+Y-Z); dividing by sin⁡(60∘)=3/2\sin(60^\circ)=\sqrt3/2 gives n^⋅σ⃗=13(X+Y−Z)\hat n\cdot\vec\sigma=\tfrac{1}{\sqrt3}(X+Y-Z), hence \boldsymbol{\hat n=(1,1,-1)/\sqrt3}.

Step 4. Cross-check with the scalar composition law (7): cos⁡(Θ/2)=cθcϕ−sθsϕ(x^⋅y^)=12⋅12−12⋅0=12\cos(\Theta/2)=c_\theta c_\phi-s_\theta s_\phi(\hat x\cdot\hat y)=\tfrac{1}{\sqrt2}\cdot\tfrac{1}{\sqrt2}-\tfrac12\cdot0=\tfrac12, agreeing with Θ=120∘\Theta=120^\circ; and n^\hat n is a unit vector since (1+1+1)/3=1(1+1+1)/3=1. Two orthogonal 90∘90^\circ turns compose into one 120∘120^\circ turn about the body-diagonal (1,1,−1)/3(1,1,-1)/\sqrt3 — a fact any circuit optimiser uses to fuse adjacent rotations.

2 · Preparing an equal superposition with a π/2 pulse

A qubit starts in ∣0⟩|0\rangle. Apply Rx(π/2)R_x(\pi/2). Find the output state, the measurement probabilities in the computational basis, and the Bloch vector.

Step 1. Apply the gate. Rx(π/2)∣0⟩=12(I−iX)∣0⟩=12(∣0⟩−i∣1⟩)R_x(\pi/2)|0\rangle=\tfrac{1}{\sqrt2}(I-iX)|0\rangle=\tfrac{1}{\sqrt2}\big(|0\rangle-i|1\rangle\big).

Step 2. Read the probabilities. P(0)=∣⟨0∣ψ⟩∣2=12P(0)=|\langle0|\psi\rangle|^2=\tfrac12 and P(1)=∣−i/2∣2=12P(1)=|{-i}/\sqrt2|^2=\tfrac12 — an equal superposition, P0=P1=12P_0=P_1=\tfrac12, just like a Hadamard but with a relative phase of −i-i instead of +1+1.

Step 3. Find the Bloch vector. A state 12(∣0⟩+eiφ∣1⟩)\tfrac{1}{\sqrt2}(|0\rangle+e^{i\varphi}|1\rangle) has Bloch vector (cos⁡φ,sin⁡φ,0)(\cos\varphi,\sin\varphi,0); here eiφ=−ie^{i\varphi}=-i so φ=−π2\varphi=-\tfrac{\pi}{2} and r⃗=(0,−1,0)=\vec r=(0,-1,0)=\boldsymbol{-\hat y}. The north pole ∣0⟩|0\rangle has rotated 90∘90^\circ about x^\hat x down to the −y^-\hat y point — exactly the trajectory drawn in the first figure, and the physical content of a π/2\pi/2 drive pulse in NMR or a superconducting qubit.

Physics of the portrait

The rotation gate is a piece of nineteenth-century mathematics wearing a twentieth-century physical costume. Hamilton's quaternions already encoded the double cover of rotations in 1843; what was missing was a physical system whose states actually transform under SU(2)SU(2) rather than SO(3)SO(3). That system arrived with the electron's spin, and the person who wrote down the generators we now exponentiate — the matrices X,Y,ZX,Y,Z — was Wolfgang Pauli.

Origins · 1 minWolfgang Pauli and how this was discovered
WP
Wolfgang Pauli
1900–1958 · the spin matrices and the algebra behind single-qubit rotations

Wolfgang Pauli was a physicist of ferocious rigour. His decisive contribution here is the 1927 paper Zur Quantenmechanik des magnetischen Elektrons, where, to describe the electron's newly discovered two-valued spin, he introduced the three 2×22\times2 matrices σx,σy,σz\sigma_x,\sigma_y,\sigma_z that now bear his name. These are precisely the generators whose exponentials are the rotation gates, Rn^(θ)=exp⁡(−iθ n^⋅σ⃗/2)R_{\hat n}(\theta)=\exp(-i\theta\,\hat n\cdot\vec\sigma/2). Pauli was not thinking of computing; he was capturing how a spin-12\tfrac12 magnetic moment responds to a field, and the half-angle that puzzles students of gates was for him the plain statement that a spin lives in the double cover of the rotation group — a 2π2\pi turn returns the direction but negates the wavefunction. His 1925 exclusion principle had already forced the electron to carry a hidden two-valued label; the spin matrices gave that label its algebra. A century on, every calibrated Rx(π/2)R_x(\pi/2) pulse on a real qubit is an experimental realisation of Pauli's matrices, and every compiler that decomposes a gate into RzRyRzR_z R_y R_z is doing group theory Pauli would have recognised at once.

History
1843
— Hamilton discovers the quaternions, the algebra of SU(2) rotations
1925
— Pauli's exclusion principle forces a hidden two-valued electron degree of freedom
1927
— Pauli introduces the 2×2 spin matrices X, Y, Z that generate the rotation gates
1931
— Wigner brings SU(2) and its two-to-one cover of SO(3) into quantum theory
1938
— Rabi drives resonant spin flips, realizing R_n(θ) as a timed pulse
1995
— Barenco and co-workers prove single-qubit rotations plus CNOT are universal

Going further

Texts. Nielsen & Chuang, Quantum Computation and Quantum Information, §4.2 (single-qubit operations) and §4.3 · Preskill, Lecture Notes on Quantum Computation, Ch. 2–3 · Watrous, The Theory of Quantum Information, Ch. 2 · Kitaev, Shen & Vyalyi, Classical and Quantum Computation, §8 · Mermin, Quantum Computer Science, Ch. 1 · Barenco et al., "Elementary gates for quantum computation", Phys. Rev. A 52, 3457 (1995). Key results to keep in mind: (1) the resummation from A2=IA^2=I yields the closed form; (2) the explicit matrices satisfy Rz(θ)=e−iθ/2diag(1,eiθ)R_z(\theta)=e^{-i\theta/2}\mathrm{diag}(1,e^{i\theta}); (3) Rn^(π)=−in^⋅σ⃗R_{\hat n}(\pi)=-i\hat n\cdot\vec\sigma, recovering the Pauli gates; (4) the conjugation Rz(θ)XRz(θ)†=cos⁡θ X+sin⁡θ YR_z(\theta)XR_z(\theta)^\dagger=\cos\theta\,X+\sin\theta\,Y via the Hadamard lemma with [Z,X]=2iY[Z,X]=2iY — note the ++ sign, confirmed by direct matrix multiplication; (5) the composition Ry(π/2)Rx(π/2)=12I−i2(X+Y−Z)R_y(\pi/2)R_x(\pi/2)=\tfrac12 I-\tfrac{i}{2}(X+Y-Z), a 120∘120^\circ turn about (1,1,−1)/3(1,1,-1)/\sqrt3; (6) the Z–Y entries are a=e−i(β+δ)/2cos⁡(γ/2)a=e^{-i(\beta+\delta)/2}\cos(\gamma/2) and b=ei(β−δ)/2sin⁡(γ/2)b=e^{i(\beta-\delta)/2}\sin(\gamma/2).

Builds on. Pauli Gates (X, Y, Z) · Bloch Sphere Representation of Gates · Quantum Circuit Model. Leads to. Z–Y Decomposition of Single-Qubit Gates · Phase Gates (S and T) · Controlled-U Gates · Universal Gate Sets.

Ask the tutor: Ask me to walk you through the Z–Y decomposition of any specific gate you name, or to show how a resonant pulse of a given area realises it on real hardware.

Why does the gate carry half the angle that the Bloch sphere rotates through?

Because a qubit is a spinor, not a direction in space. Its state lives in a complex 22-dimensional space on which SU(2)SU(2) acts, and SU(2)SU(2) is the double cover of SO(3)SO(3): the map from gates to physical rotations is two-to-one. The gate Rn^(θ)R_{\hat n}(\theta) and its negative −Rn^(θ)-R_{\hat n}(\theta) give the identical SO(3)SO(3) rotation of the Bloch vector, so a single geometric turn corresponds to two gate values differing by θ→θ+2π\theta\to\theta+2\pi. The half-angle in the exponent is exactly what makes one full geometric turn (θ=2π\theta=2\pi) send the gate to −I-I rather than back to II — you are seeing the spinor nature of the qubit directly in the arithmetic.

Is R_z(θ) the same as the phase gate diag(1, e^{iθ})?

Not as matrices, only as operations on a lone qubit. Writing them out, Rz(θ)=diag(e−iθ/2,eiθ/2)R_z(\theta)=\mathrm{diag}(e^{-i\theta/2},e^{i\theta/2}) while P(θ)=diag(1,eiθ)P(\theta)=\mathrm{diag}(1,e^{i\theta}); they differ by the overall factor e−iθ/2e^{-i\theta/2}. A global phase is unobservable, so acting alone the two give identical statistics and Bloch rotations, and you may swap them freely. The catch is control: in a controlled version the phase that was global for the target becomes a relative phase between the "control off" and "control on" branches, now fully measurable. So controlled-Rz(θ)R_z(\theta) and controlled-P(θ)P(\theta) are genuinely different gates — always check which convention a paper or compiler means by "RzR_z".

How do I get an arbitrary single-qubit gate if hardware only gives me rotations?

You use the Z–Y decomposition: every U∈U(2)U\in U(2) equals eiαRz(β)Ry(γ)Rz(δ)e^{i\alpha}R_z(\beta)R_y(\gamma)R_z(\delta) for some phase and three Euler angles. Finding them is short — read γ\gamma off the entry magnitudes and β,δ\beta,\delta off their phases. If your device's native axes are not yy and zz, an equivalent Z–X or X–Y form does the same job, since the derivation needs only two non-parallel axes. This is exactly what a compiler does: it takes the abstract unitary and emits a short list of native rotations — which is why "single-qubit universality" and "three rotations" are the same statement.

What is the physical meaning of R_n(θ) on real qubit hardware?

It is a driven oscillation. Apply a resonant pulse — microwave for superconducting qubits, laser or microwave for ions — of fixed amplitude for a time tt, and the qubit undergoes Rn^(Ωt)R_{\hat n}(\Omega t), where Ω\Omega is the Rabi frequency and the equatorial axis n^\hat n is set by the phase of the drive. A pulse area Ωt=π\Omega t=\pi gives a π\pi-pulse, a bit flip; Ωt=π/2\Omega t=\pi/2 gives the π/2\pi/2-pulse that creates superpositions and opens and closes Ramsey interferometry. A rotation gate is thus no metaphor: it is the closed-form solution of the driven two-level Schrödinger equation, and the angle is literally how long you drive.

Why is R_y the only rotation with real matrix entries, and does that matter?

It comes down to which Pauli is imaginary. Of X,Y,ZX,Y,Z, only Y=(0−ii0)Y=\begin{pmatrix}0&-i\\ i&0\end{pmatrix} is itself imaginary, so −isin⁡(θ/2) Y-i\sin(\theta/2)\,Y comes out real, leaving Ry(θ)R_y(\theta) with entirely real entries cos⁡(θ/2)\cos(\theta/2) and ±sin⁡(θ/2)\pm\sin(\theta/2); RxR_x has imaginary off-diagonals and RzR_z complex diagonal phases. This is more than a curiosity: because RyR_y keeps amplitudes real, it is the natural gate for adjusting the populations of ∣0⟩|0\rangle and ∣1⟩|1\rangle without adding relative phase — which is why the Z–Y decomposition puts RyR_y in the middle, with the two flanking RzR_z's handling phase and the single RyR_y handling the "tilt".

How do rotation gates make variational quantum algorithms trainable?

Because their angle dependence is a clean sinusoid, their gradients are exact and cheap. A rotation generator has just two eigenvalues, ±12\pm\tfrac12, so any circuit expectation ⟨H⟩(θ)\langle H\rangle(\theta) is a pure sinusoid in each angle. That gives the parameter-shift rule: the exact derivative is ∂θ⟨H⟩=12[⟨H⟩θ+π/2−⟨H⟩θ−π/2]\partial_\theta\langle H\rangle=\tfrac12[\langle H\rangle_{\theta+\pi/2}-\langle H\rangle_{\theta-\pi/2}], obtained by running the same circuit at two shifted angles — no finite-difference error, no ancillas. This turns the circuit into a differentiable function of its angles that classical optimisers can train, which is why variational ansätze are parametrised in rotation gates specifically rather than in arbitrary unitaries.

Before this

Where this leads