Skip to content

Two-Level and Finite-Dimensional Canonical Systems

Two-level systems are the smallest quantum models that support nontrivial superposition, interference, basis dependence, level repulsion, and coherent state transfer. Their state space is only two-dimensional, yet the same algebra describes localized states in a double well, two site orbitals, a selected atomic transition, spin-1/21/2, and many engineered quantum devices.

The small dimension is both the strength and the danger of the model. A 2×22\times2 Hamiltonian can often be solved exactly, but a physical system is genuinely two-level only when the neglected states remain dynamically irrelevant on the energy and time scales of interest. This chapter develops the reusable two-state structure and keeps that validity question visible throughout.

Choose an orthonormal basis {∣1⟩,∣2⟩}\{\lvert1\rangle,\lvert2\rangle\} for a two-dimensional Hilbert space. A normalized pure state is

∣ψ⟩=c1∣1⟩+c2∣2⟩,∣c1∣2+∣c2∣2=1.\lvert\psi\rangle = c_1\lvert1\rangle+c_2\lvert2\rangle, \qquad \lvert c_1\rvert^2+\lvert c_2\rvert^2=1.

Normalization removes one real parameter, and the physically irrelevant global phase removes another. A pure two-level ray therefore has two real degrees of freedom. This is why it admits a Bloch-sphere representation.

In the chosen basis, the most general Hermitian Hamiltonian is

H=(E1ΔΔ∗E2),E1,E2∈R.H = \begin{pmatrix} E_1 & \Delta\\ \Delta^* & E_2 \end{pmatrix}, \qquad E_1,E_2\in\mathbb R.

The diagonal entries are basis-state energies or detunings. The off-diagonal element Δ\Delta mixes the basis states. Its magnitude controls coherent coupling, while its phase depends partly on basis conventions.

The detailed physical setup and the conditions for using this truncation begin in Two-Level Systems. The complete static diagonalization and transition formula are developed in Two-State Hamiltonians.

Every Hermitian 2×22\times2 matrix has a unique expansion

H=c0I+b⋅σ,H = c_0I+\mathbf b\cdot\boldsymbol\sigma,

where c0∈Rc_0\in\mathbb R, b=(bx,by,bz)∈R3\mathbf b=(b_x,b_y,b_z)\in\mathbb R^3, and σ=(σx,σy,σz)\boldsymbol\sigma=(\sigma_x,\sigma_y,\sigma_z). For the matrix above,

c0=E1+E22,bz=E1−E22,bx=Re⁡Δ,by=−Im⁡Δ.\begin{aligned} c_0&=\frac{E_1+E_2}{2}, & b_z&=\frac{E_1-E_2}{2},\\ b_x&=\operatorname{Re}\Delta, & b_y&=-\operatorname{Im}\Delta. \end{aligned}

The scalar term c0Ic_0I shifts both energy levels equally. In closed-system dynamics it contributes only a global phase, so transition probabilities are governed by b⋅σ\mathbf b\cdot\boldsymbol\sigma.

Using

(b⋅σ)2=∣b∣2I,(\mathbf b\cdot\boldsymbol\sigma)^2 = \lvert\mathbf b\rvert^2I,

one obtains the two eigenvalues and spectral projectors without solving a quadratic equation component by component:

E±=c0±∣b∣,E_\pm = c_0\pm\lvert\mathbf b\rvert, P±=12(I±b^⋅σ),b^=b∣b∣.P_\pm = \frac12 \left( I\pm\hat{\mathbf b}\cdot\boldsymbol\sigma \right), \qquad \hat{\mathbf b}=\frac{\mathbf b}{\lvert\mathbf b\rvert}.

The level splitting is therefore

ΔE=E+−E−=2∣b∣.\Delta E = E_+-E_- = 2\lvert\mathbf b\rvert.

Pauli-Matrix Hamiltonians is the canonical page for this dictionary, its basis dependence, and the resulting time-evolution operator.

A two-level Hamiltonian is not meaningful until its basis has been identified. Several common bases answer different physical questions:

BasisTypical meaningWhat off-diagonal coupling does
Localized basisLeft/right well or site 1/21/2Transfers amplitude between locations
Bare internal-state basisTwo uncoupled atomic or molecular levelsDrive or interaction mixes the levels
Spin-component basisEigenstates of a selected SiS_iA transverse field mixes the outcomes
Instantaneous energy basisEigenstates of H(t)H(t) at fixed ttTime dependence can couple the branches
Diabatic basisStates continued through an uncoupled crossingFixed coupling opens an avoided crossing

A basis transformation changes the matrix entries and rotates the Pauli vector. It does not change invariant quantities such as the eigenvalues, trace, determinant, or level splitting. Statements such as “the system stays in the same state” are incomplete unless the basis or observable is specified.

This distinction is especially important near avoided crossings. Following one instantaneous energy branch can require changing from one localized or diabatic basis state to the other.

After removing global phase, any normalized pure state can be written as

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

Its Bloch vector is

r=⟨ψ∣σ∣ψ⟩=(sin⁡θcos⁡ϕ,sin⁡θsin⁡ϕ,cos⁡θ),\mathbf r = \langle\psi\vert\boldsymbol\sigma\vert\psi\rangle = (\sin\theta\cos\phi, \sin\theta\sin\phi, \cos\theta),

and its rank-one projector is

∣ψ⟩⟨ψ∣=12(I+r⋅σ).\lvert\psi\rangle\langle\psi\rvert = \frac12 \left( I+\mathbf r\cdot\boldsymbol\sigma \right).

The polar coordinate records the population imbalance in the chosen basis; the azimuthal coordinate records relative phase. The same sphere can describe a spin, a double-well state, or two internal levels, but the physical meaning of its axes changes with the realization.

For a time-independent Hamiltonian H=c0I+b⋅σH=c_0I+\mathbf b\cdot\boldsymbol\sigma, the Bloch vector obeys

drdt=2ℏ b×r.\frac{d\mathbf r}{dt} = \frac{2}{\hbar}\, \mathbf b\times\mathbf r.

Thus unitary two-level evolution is a rotation around the Hamiltonian axis at angular frequency 2∣b∣/ℏ2\lvert\mathbf b\rvert/\hbar. Bloch Sphere: Wave-Mechanics Perspective develops the state geometry, overlap formula, and measurement probabilities. Mixed states and the Bloch ball belong in Bloch Sphere for Density Operators.

Consider the real symmetric form

H(λ)=c0(λ)I+ϵ(λ)σz+Kσx.H(\lambda) = c_0(\lambda)I + \epsilon(\lambda)\sigma_z + K\sigma_x.

Its energies are

E±(λ)=c0(λ)±ϵ(λ)2+K2.E_\pm(\lambda) = c_0(\lambda) \pm \sqrt{\epsilon(\lambda)^2+K^2}.

If K=0K=0, the levels cross where ϵ=0\epsilon=0. If K≠0K\ne0, the minimum gap is 2∣K∣2\lvert K\rvert. Far from resonance, ∣ϵ∣≫∣K∣\lvert\epsilon\rvert\gg\lvert K\rvert, the eigenstates are close to the displayed basis states. Near resonance they are strongly mixed.

This simple spectrum carries several lessons:

  • off-diagonal coupling produces level repulsion in a two-state subspace;
  • the minimum gap measures the coupling in this convention;
  • the eigenvectors can change character rapidly even though the eigenvalues remain smooth;
  • an exact crossing can persist only when the relevant coupling vanishes, often because of symmetry or parameter tuning;
  • a small gap does not by itself justify discarding nearby third and higher levels.

Coupled Wells and Avoided Crossings derives this model from localized wavefunctions and tunneling. Tight-Binding Dimer gives the identical matrix in the language of onsite energies, hopping, and bonding or antibonding orbitals.

The simplest population transfer occurs for

H=Kσx,∣ψ(0)⟩=∣1⟩.H=K\sigma_x, \qquad \lvert\psi(0)\rangle=\lvert1\rangle.

Then

P2(t)=sin⁡2(Ktℏ).P_2(t) = \sin^2\left(\frac{Kt}{\hbar}\right).

This is not an irreversible jump. Probability amplitude flows coherently between the two basis states and later returns. The same equation can describe tunneling between wells, hopping across a dimer, or precession between spin-component states after translating the symbols into the relevant physical observables.

With detuning,

H=ϵσz+Kσx,H=\epsilon\sigma_z+K\sigma_x,

the transition probability becomes

P1→2(t)=K2ϵ2+K2sin⁡2(ϵ2+K2ℏt).P_{1\to2}(t) = \frac{K^2}{\epsilon^2+K^2} \sin^2\left( \frac{\sqrt{\epsilon^2+K^2}}{\hbar}t \right).

Detuning increases the oscillation frequency but suppresses the maximum transfer. This single formula is the static ancestor of the rotating-frame Rabi result.

Two canonical protocols add explicit time dependence in different ways.

The Landau–Zener model uses

H(t)=vt2σz+Δσx.H(t) = \frac{vt}{2}\sigma_z+\Delta\sigma_x.

It asks whether a state follows an instantaneous energy branch as the detuning is swept through the minimum gap. In the standard infinite-time convention, the nonadiabatic transition probability is

PLZ=exp⁡(−2πΔ2ℏv).P_{\mathrm{LZ}} = \exp\left( -\frac{2\pi\Delta^2}{\hbar v} \right).

Landau–Zener Problem: First Encounter fixes the convention and interprets the formula. Its asymptotic derivation and relation to adiabatic theory belong in Landau–Zener Transition.

For a sinusoidally driven transition, a rotating frame and the rotating-wave approximation can produce

Heff=ℏ2(δσz+Ωσx).H_{\mathrm{eff}} = \frac{\hbar}{2} \left( \delta\sigma_z+\Omega\sigma_x \right).

Starting in the lower basis state, the excited-state probability is

Pe(t)=Ω2Ω2+δ2sin⁡2(t2Ω2+δ2).P_e(t) = \frac{\Omega^2}{\Omega^2+\delta^2} \sin^2\left( \frac{t}{2}\sqrt{\Omega^2+\delta^2} \right).

Rabi Oscillations: First Encounter develops this coherent rotating-frame picture. The Rotating-Wave Approximation is an approximation with a domain of validity, not a change of notation.

The common algebra should not erase the differences among physical systems.

RealizationBasis statesCoupling or controlImportant omitted structure
Double wellLeft/right localized statesBarrier tunneling and biasHigher intrawell levels
Tight-binding dimerTwo site orbitalsHopping and onsite detuningOther sites, particles, interactions
Selected internal transitionTwo atomic or molecular levelsElectromagnetic driveSelection rules, other levels, spontaneous emission
Spin-1/21/2Two spin-component eigenstatesMagnetic fieldSpatial motion and full rotation theory
Effective qubitDevice-dependent logical statesCalibrated control fieldsLeakage, noise, measurement apparatus

Spin-1/2 as a Canonical System: First Encounter supplies the minimal dictionary Si=(ℏ/2)σiS_i=(\hbar/2)\sigma_i. The canonical treatment of spin, spinors, rotations, and magnetic-field dynamics is in Symmetry, Angular Momentum, and Spin. Effective Hamiltonians in Quantum Information explains how physical multilevel systems reduce to leakage-aware logical Hamiltonians; complete protocols belong in the planned Quantum Information and Computation volume.

When a two-level truncation is trustworthy

Section titled “When a two-level truncation is trustworthy”

Let PP project onto the proposed two-state subspace and Q=I−PQ=I-P onto all omitted states. A useful model requires more than two nearby eigenvalues. One must also control the couplings QHPQHP and any drive frequencies or bandwidths that can populate the omitted sector.

A practical audit asks:

  1. Spectral isolation. Is the separation to omitted states large compared with the retained splitting, coupling, drive amplitude, thermal scale, and inverse protocol time?
  2. Weak leakage. Are matrix elements from the retained subspace to omitted states small, forbidden by symmetry, or strongly off resonance?
  3. Controlled preparation. Does the initial state lie predominantly in the retained subspace?
  4. Controlled observation time. Can small leakage or phase errors accumulate during the experiment or calculation?
  5. Environment. Are decoherence, relaxation, and measurement backaction negligible, or have they been included in an open-system model?
  6. Parameter range. Does the isolation persist along the entire sweep, pulse, or control path rather than only at its starting point?

If an omitted state ∣n⟩\lvert n\rangle is separated by an energy Δn\Delta_n and coupled with matrix element VnV_n, the ratio

∣Vn∣∣Δn∣\frac{\lvert V_n\rvert}{\lvert\Delta_n\rvert}

is a first perturbative warning indicator, not a universal error bound. Near resonance, during strong driving, or over long times, the effective model must be rederived or enlarged.

The two-level approximation is therefore a scale-dependent statement. A system can be accurately two-level for one protocol and badly non-two-level for another.

For a first pass through the chapter:

  1. Start with Two-Level Systems for the state space and physical examples.
  2. Continue to Two-State Hamiltonians for eigenvalues, mixing, and coherent oscillation.
  3. Learn the compact invariant language in Pauli-Matrix Hamiltonians.
  4. Use Bloch Sphere: Wave-Mechanics Perspective to connect amplitudes, relative phase, measurement, and rotation.
  5. Choose Coupled Wells and Avoided Crossings or Tight-Binding Dimer for a spatial realization.
  6. Add explicit time dependence through Landau–Zener Problem: First Encounter and Rabi Oscillations: First Encounter.
  7. Use Spin-1/2 as a Canonical System: First Encounter as the bridge to the full spin volume.
PageCanonical role
Two-Level SystemsPhysical definition, state space, examples, and truncation cautions
Two-State HamiltoniansGeneral diagonalization, mixing angle, spectrum, and population dynamics
Pauli-Matrix HamiltoniansPauli decomposition, projectors, invariants, and exact propagator
Bloch Sphere: Wave-Mechanics PerspectivePure-state geometry, measurement axes, overlaps, and precession
Coupled Wells and Avoided CrossingsLocalized wavefunctions, tunneling matrix elements, and level repulsion
Tight-Binding DimerTwo-site hopping, bonding states, and the lattice-model bridge
Landau–Zener Problem: First EncounterLinear sweep, diabatic and adiabatic labels, and transition probability
Rabi Oscillations: First EncounterNear-resonant drive, rotating-frame dynamics, and pulse times
Spin-1/2 as a Canonical System: First EncounterMinimal spin dictionary and magnetic-field realization
  • Treating “two-level” as a permanent property rather than a controlled approximation for specified energies, drives, and times.
  • Calling the displayed basis the energy basis before diagonalizing the Hamiltonian.
  • Dropping c0Ic_0I in a context where absolute phases relative to another sector matter, even though it is irrelevant within an isolated two-level system.
  • Confusing the state vector, the ray, and the Bloch vector.
  • Forgetting that the phase of an isolated off-diagonal coupling can often be changed by rephasing the basis states.
  • Interpreting coherent hopping or Rabi cycling as irreversible transition rates.
  • Quoting a Landau–Zener probability without stating whether it describes adiabatic or diabatic outcomes and which slope convention defines vv.
  • Assuming every two-level system is physically a spin-1/21/2 system.
  • Applying a closed-system model when relaxation, dephasing, or measurement backaction controls the observed dynamics.
  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • L. Allen and J. H. Eberly, Optical Resonance and Two-Level Atoms, Dover, 1987.
  • L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, 3rd ed., Pergamon, 1977.
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, 2010.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013.
  1. Write
H=(31+i1−i−1)H = \begin{pmatrix} 3 & 1+i\\ 1-i & -1 \end{pmatrix}

in the form c0I+b⋅σc_0I+\mathbf b\cdot\boldsymbol\sigma, and find its eigenvalues.

Solution

Using

(c0+bzbx−ibybx+ibyc0−bz),\begin{pmatrix} c_0+b_z & b_x-ib_y\\ b_x+ib_y & c_0-b_z \end{pmatrix},

the diagonal entries give

c0=1,bz=2.c_0=1, \qquad b_z=2.

The off-diagonal entry bx−iby=1+ib_x-ib_y=1+i gives

bx=1,by=−1.b_x=1, \qquad b_y=-1.

Therefore

H=I+σx−σy+2σz.H = I+\sigma_x-\sigma_y+2\sigma_z.

Since ∣b∣=6\lvert\mathbf b\rvert=\sqrt6, the eigenvalues are

E±=1±6.E_\pm=1\pm\sqrt6.
  1. A Hamiltonian is H=ϵσz+KσxH=\epsilon\sigma_z+K\sigma_x. Show that its minimum spectral gap as ϵ\epsilon varies is 2∣K∣2\lvert K\rvert, and identify the corresponding eigenstates when K>0K\gt0.
Solution

The eigenvalues are

E±=±ϵ2+K2,E_\pm = \pm\sqrt{\epsilon^2+K^2},

so the gap is

ΔE=2ϵ2+K2.\Delta E = 2\sqrt{\epsilon^2+K^2}.

It is minimized at ϵ=0\epsilon=0, where ΔEmin⁡=2∣K∣\Delta E_{\min}=2\lvert K\rvert. At that point H=KσxH=K\sigma_x. For K>0K\gt0, the upper and lower eigenstates are respectively

∣+x⟩=∣1⟩+∣2⟩2,∣−x⟩=∣1⟩−∣2⟩2.\lvert+_x\rangle = \frac{\lvert1\rangle+\lvert2\rangle}{\sqrt2}, \qquad \lvert-_x\rangle = \frac{\lvert1\rangle-\lvert2\rangle}{\sqrt2}.
  1. A state has Bloch vector r=(0,3/5,4/5)\mathbf r=(0,3/5,4/5). Find the probabilities for measuring σz=±1\sigma_z=\pm1, and explain whether the state is pure.
Solution

The probabilities are

p±=12(1±rz).p_\pm = \frac12(1\pm r_z).

Thus

p+=910,p−=110.p_+=\frac{9}{10}, \qquad p_-=\frac{1}{10}.

The vector length is

∣r∣=925+1625=1.\lvert\mathbf r\rvert = \sqrt{\frac{9}{25}+\frac{16}{25}} = 1.

It therefore represents a pure state. A Bloch vector of length less than one would represent a mixed two-level density operator.

  1. A proposed two-level model retains states ∣1⟩\lvert1\rangle and ∣2⟩\lvert2\rangle but omits ∣3⟩\lvert3\rangle. The drive is resonant with the 1↔21\leftrightarrow2 transition, while the 2↔32\leftrightarrow3 detuning is Δ3\Delta_3 and its drive matrix element is V3V_3. Give a qualitative validity condition and name two additional checks needed before trusting the truncation.
Solution

A first off-resonant condition is

∣V3∣≪∣Δ3∣.\lvert V_3\rvert \ll \lvert\Delta_3\rvert.

This suppresses direct population of ∣3⟩\lvert3\rangle in a weak-drive perturbative regime. It is not sufficient by itself. One should also check, for example, that the pulse bandwidth does not overlap the omitted transition, that the evolution time is not long enough for small leakage to accumulate, that the initial state has negligible ∣3⟩\lvert3\rangle component, and that other omitted states or environmental processes are not more important. Strong driving can invalidate the condition even when the undriven spectrum appears well isolated.