Time Evolution Operators

In quantum mechanics, the state of a closed system evolves in time by a unitary operator. If the system's energy is described by a Hermitian operator called the Hamiltonian $H$ , the Schrodinger equation gives the time-evolution operator

$U(t) = e^{-iHt}$

(setting $\hbar = 1$ throughout). Starting from state $|\psi_0\rangle$ , after time $t$ the system is in the state $|\psi(t)\rangle = U(t)|\psi_0\rangle$ .

Why the matrix exponential?

For any Hermitian matrix $H$ with eigenvalues $E_k$ and eigenstates $|E_k\rangle$ ,

$e^{-iHt} = \sum_k e^{-iE_k t}|E_k\rangle\langle E_k|.$

Each energy eigenstate acquires a phase $e^{-iE_k t}$ ; superpositions of different energies therefore oscillate relative to each other — that oscillation is quantum dynamics.

Evolution under a single Pauli

The simplest nontrivial Hamiltonians are proportional to the Pauli matrices $X$ , $Y$ , $Z$ . Because $Z^2 = I$ , its matrix exponential closes into trigonometric functions:

$e^{-i\theta Z} = \cos\theta\, I - i\sin\theta\, Z = \begin{pmatrix} e^{-i\theta} & 0 \\ 0 & e^{i\theta} \end{pmatrix}.$

Comparing with the standard definition of the $R_z$ rotation gate,

$R_z(\phi) = \begin{pmatrix} e^{-i\phi/2} & 0 \\ 0 & e^{i\phi/2} \end{pmatrix} = e^{-i(\phi/2)Z},$

we get the identity

$e^{-i\theta Z} = R_z(2\theta).$

This is the key bridge between physics (Hamiltonian evolution) and gate-level quantum programming: evolving for time $t$ under the Hamiltonian $H = Z$ with coupling strength $\theta/t$ is exactly the same as applying an $R_z$ gate with angle $2\theta$ .

A concrete example

Take $\theta = \pi/4$ . Then

$e^{-i(\pi/4)Z} = \begin{pmatrix} e^{-i\pi/4} & 0 \\ 0 & e^{i\pi/4} \end{pmatrix},$

which equals $R_z(\pi/2)$ . This is sometimes called the $S$ gate (up to a global phase), and it rotates the Bloch sphere by $90°$ around the $z$ -axis.

Why this matters for algorithms

Simulating the dynamics of a physical system (molecules, materials, quantum magnets) on a quantum computer reduces to implementing $e^{-iHt}$ as a sequence of gates. For multi-term Hamiltonians such as $H = \sum_j h_j P_j$ (a sum of Pauli products), the key technique is Trotterization:

$e^{-i(A+B)t} \approx \left(e^{-iAt/n}\, e^{-iBt/n}\right)^n,$

which becomes exact as $n \to \infty$ . Each factor $e^{-iP_j \Delta t}$ can then be compiled into one or a few primitive gates. Upcoming lessons build on this foundation.

Try it

Implement the unitary $e^{-i(\pi/4)Z}$ on a single qubit. Use the rotation-gate identity above to find the right call. The grader checks the full unitary matrix, so a do-nothing circuit will fail.

Run your code to see the quantum state.