Skip to content

Unitary Time Evolution

Unitary time evolution is the reversible, inner-product-preserving dynamics of a closed quantum system. If the state is known at a reference time t0t_0, then

∣ψ(t)⟩=U(t,t0)∣ψ(t0)⟩,U†(t,t0)U(t,t0)=I.\begin{aligned} \lvert\psi(t)\rangle &= U(t,t_0)\lvert\psi(t_0)\rangle, \\ U^\dagger(t,t_0)U(t,t_0) &=I. \end{aligned}

The second relation is not merely a normalization trick. It says that evolution preserves the entire geometry of Hilbert space: norms, angles, orthogonality, and transition probabilities. For a well-posed Schrödinger equation, this structure follows from the self-adjointness of the Hamiltonian.

This page develops that core statement. Detailed construction and composition of U(t,t0)U(t,t_0) belong to Time-Evolution Operator, while the functional-analytic converse belongs to Stone Theorem.

A linear operator UU on a Hilbert space H\mathcal H is unitary when

U†U=UU†=I.U^\dagger U = UU^\dagger = I.

Equivalently, UU is a surjective linear map that preserves inner products. Its adjoint is its inverse:

U−1=U†.U^{-1}=U^\dagger.

Consequently, a unitary evolution can always be reversed at the level of the closed-system state. For a two-time evolution operator,

U(t,t0)−1=U(t0,t)=U†(t,t0).U(t,t_0)^{-1} = U(t_0,t) = U^\dagger(t,t_0).

Linearity also preserves superposition. If

∣χ(t0)⟩=a∣ψ(t0)⟩+b∣ϕ(t0)⟩,\lvert\chi(t_0)\rangle = a\lvert\psi(t_0)\rangle + b\lvert\phi(t_0)\rangle,

then

∣χ(t)⟩=a∣ψ(t)⟩+b∣ϕ(t)⟩.\lvert\chi(t)\rangle = a\lvert\psi(t)\rangle + b\lvert\phi(t)\rangle.

In finite dimensions, U†U=IU^\dagger U=I already implies UU†=IUU^\dagger=I. In an infinite-dimensional Hilbert space, the first identity alone says that UU is an isometry; surjectivity must also be established. The unilateral shift in Exercise 7 is the standard counterexample.

For the abstract operator theory behind these statements, see Unitary Operators.

Let two states evolve under the same U(t,t0)U(t,t_0):

∣ψ(t)⟩=U(t,t0)∣ψ0⟩,∣ϕ(t)⟩=U(t,t0)∣ϕ0⟩.\begin{aligned} \lvert\psi(t)\rangle &= U(t,t_0)\lvert\psi_0\rangle,\\ \lvert\phi(t)\rangle &= U(t,t_0)\lvert\phi_0\rangle. \end{aligned}

Their inner product is unchanged:

⟨ϕ(t)∣ψ(t)⟩=⟨ϕ0∣U†(t,t0)U(t,t0)∣ψ0⟩=⟨ϕ0∣ψ0⟩.\begin{aligned} \langle\phi(t)\vert\psi(t)\rangle &= \langle\phi_0\vert U^\dagger(t,t_0)U(t,t_0) \vert\psi_0\rangle\\ &= \langle\phi_0\vert\psi_0\rangle. \end{aligned}

Several consequences follow at once.

Setting ∣ϕ0⟩=∣ψ0⟩\lvert\phi_0\rangle=\lvert\psi_0\rangle gives

∥ψ(t)∥2=∥ψ0∥2.\lVert\psi(t)\rVert^2 = \lVert\psi_0\rVert^2.

A normalized state therefore remains normalized. In wave mechanics this is the global statement that the integral of the probability density remains one. Its local refinement is the continuity equation, derived in Unitarity and Conservation of Probability.

If ⟨ϕ0∣ψ0⟩=0\langle\phi_0\vert\psi_0\rangle=0, then the evolved states remain orthogonal. More generally,

∣⟨ϕ(t)∣ψ(t)⟩∣2=∣⟨ϕ0∣ψ0⟩∣2.\left| \langle\phi(t)\vert\psi(t)\rangle \right|^2 = \left| \langle\phi_0\vert\psi_0\rangle \right|^2.

Thus two closed-system states do not become more or less distinguishable merely because both undergo the same unitary evolution. A unitary may rotate their representation, but it cannot change their overlap.

Unitary maps preserve linear independence, orthonormal bases, and distances:

∥U∣ψ⟩−U∣ϕ⟩∥=∥∣ψ⟩−∣ϕ⟩∥.\left\lVert U\lvert\psi\rangle - U\lvert\phi\rangle \right\rVert = \left\lVert \lvert\psi\rangle - \lvert\phi\rangle \right\rVert.

They may nevertheless change relative phases between components of a single state. Those phases can alter later interference probabilities even though the state remains normalized.

Unitarity does not imply that the probability of every fixed measurement outcome is constant. For a fixed projector Πa\Pi_a,

pa(t)=⟨ψ0∣U†(t,t0)ΠaU(t,t0)∣ψ0⟩,p_a(t) = \langle\psi_0\vert U^\dagger(t,t_0)\Pi_a U(t,t_0) \vert\psi_0\rangle,

which generally depends on time. The probability is constant for all initial states only when the relevant projector is invariant under the evolution. For a time-independent Hamiltonian, this is closely related to

[H,Πa]=0.[H,\Pi_a]=0.

The general criterion, including explicit time dependence, is developed in Conservation Laws.

From the Schrödinger Equation to Unitarity

Section titled “From the Schrödinger Equation to Unitarity”

Suppose ∣ψ(t)⟩\lvert\psi(t)\rangle and ∣ϕ(t)⟩\lvert\phi(t)\rangle satisfy the same Schrödinger equation,

iℏ∣ψ˙⟩=H∣ψ⟩,iℏ∣ϕ˙⟩=H∣ϕ⟩.i\hbar\lvert\dot\psi\rangle = H\lvert\psi\rangle, \qquad i\hbar\lvert\dot\phi\rangle = H\lvert\phi\rangle.

Differentiating their inner product gives

ddt⟨ϕ∣ψ⟩=⟨ϕ˙∣ψ⟩+⟨ϕ∣ψ˙⟩=iℏ⟨ϕ∣(H†−H)∣ψ⟩.\begin{aligned} \frac{d}{dt} \langle\phi\vert\psi\rangle &= \langle\dot\phi\vert\psi\rangle + \langle\phi\vert\dot\psi\rangle\\ &= \frac{i}{\hbar} \langle\phi\vert \left(H^\dagger-H\right) \vert\psi\rangle. \end{aligned}

If H=H†H=H^\dagger, the derivative vanishes. The Schrödinger equation therefore preserves every inner product, not only the norm of one chosen state.

The same argument can be written directly for the evolution operator. If

iℏ∂U∂t=H(t)U,U(t0,t0)=I,i\hbar\frac{\partial U}{\partial t} = H(t)U, \qquad U(t_0,t_0)=I,

and H(t)H(t) is self-adjoint, then

∂∂t(U†U)=iℏU†HU−iℏU†HU=0.\begin{aligned} \frac{\partial}{\partial t} \left(U^\dagger U\right) &= \frac{i}{\hbar}U^\dagger H U - \frac{i}{\hbar}U^\dagger H U\\ &=0. \end{aligned}

Because U†U=IU^\dagger U=I at t=t0t=t_0, it remains II throughout the evolution.

In finite-dimensional systems, this calculation is ordinarily sufficient. For unbounded Hamiltonians, symbols such as H†=HH^\dagger=H conceal domain questions: the Schrödinger equation must generate a well-defined evolution on the Hilbert space. This is why the precise requirement is self-adjointness, not merely symmetry of matrix elements on an unspecified domain. See Hermitian vs. Self-Adjoint Operators for the distinction.

The implication also works locally in the other direction. Write an infinitesimal evolution as

U(t+δt,t)=I−iℏG(t) δt+O(δt2).U(t+\delta t,t) = I - \frac{i}{\hbar}G(t)\,\delta t + O(\delta t^2).

Then

U†U=I+iℏ(G†−G)δt+O(δt2).U^\dagger U = I + \frac{i}{\hbar} \left(G^\dagger-G\right)\delta t + O(\delta t^2).

Unitarity to first order requires

G†=G.G^\dagger=G.

The generator of differentiable unitary time evolution is therefore self-adjoint. In quantum dynamics that generator is the Hamiltonian.

For a strongly continuous one-parameter unitary group, Stone Theorem makes this correspondence exact: there is a unique self-adjoint generator HH such that

U(τ)=exp⁡ ⁣(−iℏHτ).U(\tau) = \exp\!\left( -\frac{i}{\hbar}H\tau \right).

That group form applies directly to time-translation-invariant dynamics. A driven Hamiltonian H(t)H(t) instead supplies a self-adjoint instantaneous generator; its two-time propagator need not depend only on t−t0t-t_0.

Time dependence does not by itself destroy unitarity. If H(t)H(t) is self-adjoint and the time-dependent Schrödinger problem is well posed, then U(t,t0)U(t,t_0) remains unitary.

What changes is the formula used to construct it. In general,

U(t,t0)≠exp⁡ ⁣[−iℏ∫t0tH(s) ds]U(t,t_0) \ne \exp\!\left[ -\frac{i}{\hbar} \int_{t_0}^{t}H(s)\,ds \right]

when Hamiltonians at different times fail to commute:

[H(t1),H(t2)]≠0.[H(t_1),H(t_2)]\ne 0.

Time ordering is then required. The canonical derivations are in Time-Evolution Operator and Time-Dependent Hamiltonians.

Unitary evolution acts on a density operator by conjugation:

ρ(t)=U(t,t0)ρ(t0)U†(t,t0).\rho(t) = U(t,t_0)\rho(t_0)U^\dagger(t,t_0).

This transformation preserves the defining properties of a density operator.

Trace. Cyclicity gives

Tr⁡ρ(t)=Tr⁡(U†Uρ(t0))=Tr⁡ρ(t0).\operatorname{Tr}\rho(t) = \operatorname{Tr} \left( U^\dagger U\rho(t_0) \right) = \operatorname{Tr}\rho(t_0).

Positivity. For every ∣χ⟩\lvert\chi\rangle,

⟨χ∣ρ(t)∣χ⟩=⟨U†χ∣ρ(t0)∣U†χ⟩≥0.\langle\chi\vert\rho(t)\vert\chi\rangle = \langle U^\dagger\chi\vert \rho(t_0) \vert U^\dagger\chi\rangle \ge 0.

Spectrum. If

ρ(t0)∣rj⟩=rj∣rj⟩,\rho(t_0)\lvert r_j\rangle = r_j\lvert r_j\rangle,

then

ρ(t)U∣rj⟩=rjU∣rj⟩.\rho(t)U\lvert r_j\rangle = r_jU\lvert r_j\rangle.

The eigenvalues {rj}\{r_j\} are unchanged. It follows that unitary evolution preserves purity and von Neumann entropy:

Tr⁡ρ(t)2=Tr⁡ρ(t0)2,S ⁣(ρ(t))=S ⁣(ρ(t0)).\begin{aligned} \operatorname{Tr}\rho(t)^2 &= \operatorname{Tr}\rho(t_0)^2,\\ S\!\left(\rho(t)\right) &= S\!\left(\rho(t_0)\right). \end{aligned}

A pure state therefore cannot become a mixed state under a unitary acting on that system alone. A subsystem can become mixed, however, when a larger closed system evolves unitarily and develops entanglement. Tracing out the environment discards correlations and is not a unitary operation on the subsystem.

Consider

H=ℏΩ2σz,τ=t−t0.H = \frac{\hbar\Omega}{2}\sigma_z, \qquad \tau=t-t_0.

The unitary evolution is

U(τ)=(e−iΩτ/200eiΩτ/2).U(\tau) = \begin{pmatrix} e^{-i\Omega\tau/2} & 0\\ 0 & e^{i\Omega\tau/2} \end{pmatrix}.

Start from the σx\sigma_x eigenstate

∣+⟩=∣0⟩+∣1⟩2.\lvert+\rangle = \frac{ \lvert0\rangle+\lvert1\rangle }{\sqrt2}.

At time tt,

∣ψ(t)⟩=12(e−iΩτ/2∣0⟩+eiΩτ/2∣1⟩).\lvert\psi(t)\rangle = \frac{1}{\sqrt2} \left( e^{-i\Omega\tau/2}\lvert0\rangle + e^{i\Omega\tau/2}\lvert1\rangle \right).

The probabilities in the zz basis remain

p0(t)=p1(t)=12.p_0(t)=p_1(t)=\frac12.

In the xx basis, however,

p+(t)=cos⁡2 ⁣(Ωτ2),p−(t)=sin⁡2 ⁣(Ωτ2).\begin{aligned} p_+(t) &= \cos^2\!\left(\frac{\Omega\tau}{2}\right),\\ p_-(t) &= \sin^2\!\left(\frac{\Omega\tau}{2}\right). \end{aligned}

Nothing nonunitary has occurred. The norm is fixed, but the Hamiltonian creates a changing relative phase, so probabilities in a basis not aligned with HH oscillate. On the Bloch sphere, this is a rotation about the zz axis.

Unitary evolution is the dynamical rule for a closed system over the interval being modeled. It is not the most general rule for every state change encountered in practice.

Suppose a system SS interacts with an environment EE. The joint state may evolve unitarily under VSEV_{SE}, while the reduced state follows

ρS(t)=Tr⁡E[VSE(ρS(t0)⊗ρE)VSE†].\rho_S(t) = \operatorname{Tr}_E \left[ V_{SE} \left( \rho_S(t_0)\otimes\rho_E \right) V_{SE}^\dagger \right].

This reduced map is generally a quantum channel, not conjugation by a unitary on SS. It can decrease purity, increase entropy, and make initially distinct subsystem states less distinguishable. The missing information resides in system–environment correlations.

Three other cases must also be separated from closed-system unitary evolution:

  • Selective measurement update conditions on a recorded outcome and is generally nonunitary and nonlinear after normalization.
  • Nonselective measurement and noise are described by completely positive trace-preserving maps, of which unitary channels are a special case.
  • Effective non-self-adjoint Hamiltonians can describe loss, decay, postselection, or a restricted sector. Their changing norm signals that some part of the physical description has been omitted or conditioned upon.

The canonical entry point for these general maps is Quantum Channels and Noise. Measurement-specific state changes are introduced in Generalized Measurements Overview.

StatementCorrect conclusion
UU is unitaryInner products and norms are preserved.
The state evolves unitarilyEvolution of the closed-system state is reversible.
A fixed observable is measured laterIts outcome probabilities may change with time.
[A,H]=0[A,H]=0 for time-independent AA and HHThe statistics of AA are conserved.
A global state evolves unitarilyA reduced subsystem need not evolve unitarily.
A state acquires phasesRelative phases may change interference; a common global phase does not.
The Hamiltonian is time dependentEvolution can still be unitary, but time ordering may be needed.
  • Checking only norm preservation for one state instead of inner-product preservation for all states.
  • Assuming U†U=IU^\dagger U=I automatically implies UU†=IUU^\dagger=I in every infinite-dimensional setting.
  • Calling a symmetric differential operator self-adjoint without specifying its domain and boundary conditions.
  • Interpreting unitary evolution as conservation of every observable.
  • Treating the ordinary exponential of ∫H(t) dt\int H(t)\,dt as valid when Hamiltonians at different times do not commute.
  • Expecting a reduced subsystem to evolve unitarily merely because the composite system does.
  • Describing selective measurement update as ordinary Hamiltonian evolution.
  • Confusing a physically irrelevant global phase with a dynamically relevant relative phase.
  • J. von Neumann, Mathematical Foundations of Quantum Mechanics, Princeton University Press, 1955.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020, ch. 2.
  • L. E. Ballentine, Quantum Mechanics: A Modern Development, World Scientific, 1998, chs. 3–4.
  • B. C. Hall, Quantum Theory for Mathematicians, Springer, 2013, chs. 5–6.
  • M. Reed and B. Simon, Methods of Modern Mathematical Physics, Vol. I: Functional Analysis, revised ed., Academic Press, 1980.
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, 10th anniversary ed., Cambridge University Press, 2010, ch. 8.
  1. Let UU act on a finite-dimensional Hilbert space and satisfy U†U=IU^\dagger U=I. Prove that UU is unitary and preserves every inner product.
Solution

For any ∣ψ⟩\lvert\psi\rangle,

∥U∣ψ⟩∥2=⟨ψ∣U†U∣ψ⟩=∥ψ∥2.\lVert U\lvert\psi\rangle\rVert^2 = \langle\psi\vert U^\dagger U\vert\psi\rangle = \lVert\psi\rVert^2.

Thus the kernel of UU is trivial, so UU is injective. In a finite-dimensional vector space, an injective linear map from the space to itself is also surjective. Hence U−1U^{-1} exists. Multiplying U†U=IU^\dagger U=I on the right by U−1U^{-1} gives U†=U−1U^\dagger=U^{-1} and therefore UU†=IUU^\dagger=I.

For arbitrary ∣ϕ⟩\lvert\phi\rangle and ∣ψ⟩\lvert\psi\rangle,

⟨Uϕ∣Uψ⟩=⟨ϕ∣U†U∣ψ⟩=⟨ϕ∣ψ⟩.\langle U\phi\vert U\psi\rangle = \langle\phi\vert U^\dagger U\vert\psi\rangle = \langle\phi\vert\psi\rangle.
  1. Two normalized states have transition probability
Pϕ←ψ=∣⟨ϕ∣ψ⟩∣2.P_{\phi\leftarrow\psi} = \left| \langle\phi\vert\psi\rangle \right|^2.

Show that applying the same unitary UU to both states leaves this probability unchanged. Why does this not imply that the probability of a fixed detector outcome is time independent?

Solution

The evolved overlap is

⟨Uϕ∣Uψ⟩=⟨ϕ∣U†U∣ψ⟩=⟨ϕ∣ψ⟩,\langle U\phi\vert U\psi\rangle = \langle\phi\vert U^\dagger U\vert\psi\rangle = \langle\phi\vert\psi\rangle,

so its absolute square is unchanged.

A fixed detector outcome corresponds to a fixed projector Π\Pi. Its probability is

p(t)=⟨ψ0∣U†(t,t0)ΠU(t,t0)∣ψ0⟩,p(t) = \langle\psi_0\vert U^\dagger(t,t_0)\Pi U(t,t_0) \vert\psi_0\rangle,

which can vary because the detector state is not being evolved alongside ∣ψ(t)⟩\lvert\psi(t)\rangle. Only the overlap of two states subjected to the same unitary is guaranteed to be constant.

  1. Suppose an infinitesimal evolution has the form
U(δt)=I−iℏG δt+O(δt2).U(\delta t) = I-\frac{i}{\hbar}G\,\delta t +O(\delta t^2).

Use unitarity to show that GG must be self-adjoint.

Solution

Taking the adjoint gives

U†(δt)=I+iℏG†δt+O(δt2).U^\dagger(\delta t) = I+\frac{i}{\hbar}G^\dagger\delta t +O(\delta t^2).

Therefore

U†U=I+iℏ(G†−G)δt+O(δt2).U^\dagger U = I + \frac{i}{\hbar} \left(G^\dagger-G\right)\delta t + O(\delta t^2).

The coefficient of δt\delta t must vanish if U†U=IU^\dagger U=I, so G†=GG^\dagger=G.

  1. For the qubit Hamiltonian
H=ℏΩ2σz,H=\frac{\hbar\Omega}{2}\sigma_z,

start from ∣0⟩\lvert0\rangle. Show that all σz\sigma_z probabilities are constant. Then start from ∣+⟩\lvert+\rangle and find the probability of obtaining +1+1 in a σx\sigma_x measurement after time τ\tau.

Solution

Because ∣0⟩\lvert0\rangle is an energy eigenstate,

U(τ)∣0⟩=e−iΩτ/2∣0⟩.U(\tau)\lvert0\rangle = e^{-i\Omega\tau/2}\lvert0\rangle.

Only a global phase changes, so p0=1p_0=1 and p1=0p_1=0 at every time.

For ∣+⟩\lvert+\rangle,

⟨+∣U(τ)∣+⟩=12(e−iΩτ/2+eiΩτ/2)=cos⁡ ⁣(Ωτ2).\begin{aligned} \langle+\vert U(\tau)\vert+\rangle &= \frac12 \left( e^{-i\Omega\tau/2} + e^{i\Omega\tau/2} \right)\\ &= \cos\!\left(\frac{\Omega\tau}{2}\right). \end{aligned}

Hence

p+(τ)=cos⁡2 ⁣(Ωτ2).p_+(\tau) = \cos^2\!\left(\frac{\Omega\tau}{2}\right).

The evolution is unitary even though this fixed-basis probability oscillates.

  1. Let ρ′=UρU†\rho'=U\rho U^\dagger. Prove that ρ′\rho' has the same eigenvalues as ρ\rho. Deduce that purity and von Neumann entropy are invariant.
Solution

If ρ∣rj⟩=rj∣rj⟩\rho\lvert r_j\rangle=r_j\lvert r_j\rangle, then

ρ′U∣rj⟩=UρU†U∣rj⟩=rjU∣rj⟩.\begin{aligned} \rho' U\lvert r_j\rangle &= U\rho U^\dagger U\lvert r_j\rangle\\ &= r_jU\lvert r_j\rangle. \end{aligned}

Thus each eigenvalue rjr_j is preserved, including multiplicity. Both

Tr⁡ρ2=∑jrj2\operatorname{Tr}\rho^2 = \sum_j r_j^2

and

S(ρ)=−∑jrjln⁡rjS(\rho) = -\sum_j r_j\ln r_j

depend only on these eigenvalues, so they are unchanged.

  1. Let U(t,t0)U(t,t_0) satisfy
iℏ∂U∂t=H(t)Ui\hbar\frac{\partial U}{\partial t} = H(t)U

with H(t)=H†(t)H(t)=H^\dagger(t). Derive directly that U†(t,t0)U(t,t0)=IU^\dagger(t,t_0)U(t,t_0)=I.

Solution

The evolution equation and its adjoint imply

∂U∂t=−iℏH(t)U,∂U†∂t=iℏU†H(t).\begin{aligned} \frac{\partial U}{\partial t} &= -\frac{i}{\hbar}H(t)U, \\ \frac{\partial U^\dagger}{\partial t} &= \frac{i}{\hbar}U^\dagger H(t). \end{aligned}

Therefore

∂∂t(U†U)=∂U†∂tU+U†∂U∂t=iℏU†HU−iℏU†HU=0.\begin{aligned} \frac{\partial}{\partial t} \left(U^\dagger U\right) &= \frac{\partial U^\dagger}{\partial t}U + U^\dagger\frac{\partial U}{\partial t}\\ &= \frac{i}{\hbar}U^\dagger H U - \frac{i}{\hbar}U^\dagger H U\\ &=0. \end{aligned}

Since U(t0,t0)=IU(t_0,t_0)=I, the constant product is U†U=IU^\dagger U=I.

  1. On the Hilbert space ℓ2(N0)\ell^2(\mathbb N_0) with orthonormal basis {∣n⟩}n=0∞\{\lvert n\rangle\}_{n=0}^{\infty}, define the unilateral shift by
S∣n⟩=∣n+1⟩.S\lvert n\rangle = \lvert n+1\rangle.

Show that S†S=IS^\dagger S=I but SS†≠ISS^\dagger\ne I. What lesson does this give about infinite-dimensional isometries?

Solution

The adjoint acts as

S†∣0⟩=0,S†∣n⟩=∣n−1⟩(n≥1).\begin{aligned} S^\dagger\lvert0\rangle &=0,\\ S^\dagger\lvert n\rangle &= \lvert n-1\rangle \quad(n\ge1). \end{aligned}

Hence

S†S∣n⟩=∣n⟩S^\dagger S\lvert n\rangle = \lvert n\rangle

for every nn, so S†S=IS^\dagger S=I. But

SS†∣0⟩=0≠∣0⟩.SS^\dagger\lvert0\rangle=0 \ne \lvert0\rangle.

In fact,

SS†=I−∣0⟩⟨0∣.SS^\dagger = I-\lvert0\rangle\langle0\rvert.

The shift preserves all norms but is not surjective, so it is an isometry rather than a unitary operator. In infinite dimensions, the two identities in the unitary definition are not automatically equivalent.

  1. A joint unitary on two qubits maps ∣00⟩\lvert00\rangle to the Bell state
∣Φ+⟩=∣00⟩+∣11⟩2.\lvert\Phi^+\rangle = \frac{ \lvert00\rangle+\lvert11\rangle }{\sqrt2}.

Find the reduced state of the first qubit before and after the evolution. Explain why its reduced evolution cannot be unitary.

Solution

Initially,

ρA(t0)=∣0⟩⟨0∣,\rho_A(t_0) = \lvert0\rangle\langle0\rvert,

which is pure. After the joint evolution,

ρA(t)=Tr⁡B(∣Φ+⟩⟨Φ+∣)=12(∣0⟩⟨0∣+∣1⟩⟨1∣)=I2.\begin{aligned} \rho_A(t) &= \operatorname{Tr}_B \left( \lvert\Phi^+\rangle\langle\Phi^+\rvert \right)\\ &= \frac12 \left( \lvert0\rangle\langle0\rvert + \lvert1\rangle\langle1\rvert \right)\\ &= \frac{I}{2}. \end{aligned}

The purity changes from 11 to 1/21/2. Conjugation by a unitary on the first qubit would preserve the eigenvalues and purity of its density operator, so no unitary acting on that qubit alone can produce this change. The composite evolution is unitary; the reduced evolution is not.