Skip to content

Sudden Approximation

The sudden approximation describes a Hamiltonian change completed so quickly that the quantum state has negligible time to evolve during the change. Immediately after an ideal sudden quench, the ket is therefore the same vector as immediately before it. What changes is the Hamiltonian, its eigenbasis, and the energy distribution assigned to that unchanged state.

This page is the canonical home for abrupt Hamiltonian changes as an approximation method. Quantum Quenches owns the subsequent many-body initial-value protocol, correlation fronts, entanglement growth, and return amplitudes. The exact framework for arbitrary time dependence lives in Time-Dependent Hamiltonians. Work Distributions owns the thermodynamic measurement protocol, while Squeezed States: First Encounter owns the general squeeze-operator formalism.

Let a control parameter change from λi\lambda_i to λf\lambda_f over a quench interval

ti≤t≤tf,τq≡tf−ti,t_i\le t\le t_f, \qquad \tau_q\equiv t_f-t_i,

and let

H(t)=H(λ(t)).H(t)=H\big(\lambda(t)\big).

The exact evolution across the interval is

Uq≡U(tf,ti)=Texp⁡[−iℏ∫titfH(t) dt].\begin{aligned} U_q &\equiv U(t_f,t_i) \\ &= \mathcal T \exp\left[ -\frac{i}{\hbar} \int_{t_i}^{t_f}H(t)\,dt \right]. \end{aligned}

The sudden approximation is

Uq≃eiϕqIU_q \simeq e^{i\phi_q}I

on the states or subspace relevant to the calculation. The scalar phase ϕq\phi_q has no effect on probabilities. Thus the physical statement is not that the Hamiltonian changes by a small amount. The change Hi→HfH_i\to H_f may be large. The statement is that the integrated nontrivial action accumulated during the short switching interval is small.

A common time-scale estimate is

τqΔEcharℏ≪1,\frac{\tau_q\Delta E_{\mathrm{char}}}{\hbar} \ll 1,

or equivalently

τq≪1ωchar,ωchar=ΔEcharℏ.\tau_q \ll \frac{1}{\omega_{\mathrm{char}}}, \qquad \omega_{\mathrm{char}} = \frac{\Delta E_{\mathrm{char}}}{\hbar}.

The phrase “characteristic energy” must be tied to the prepared state and to the matrix elements activated by the protocol. It is not automatically the smallest gap, the largest eigenvalue of an unbounded Hamiltonian, or the final level spacing.

For a bounded Hamiltonian on a finite-dimensional space or controlled subspace, separate an arbitrary scalar part:

H(t)=Eref(t)I+K(t).H(t) = E_{\mathrm{ref}}(t)I + K(t).

Define the integrated action

Aq≡1ℏ∫titf∥K(t)∥ dt.\mathcal A_q \equiv \frac{1}{\hbar} \int_{t_i}^{t_f} \lVert K(t)\rVert\,dt.

The scalar term contributes

Φq=1ℏ∫titfEref(t) dt.\Phi_q = \frac{1}{\hbar} \int_{t_i}^{t_f} E_{\mathrm{ref}}(t)\,dt.

Duhamel’s formula gives the useful bound

∥Uq−e−iΦqI∥≤Aq.\left\lVert U_q-e^{-i\Phi_q}I \right\rVert \le \mathcal A_q.

Therefore

Aq≪1\mathcal A_q\ll1

is a sufficient suddenness condition on that space. Subtracting ErefIE_{\mathrm{ref}}I matters: a large common energy offset generates only a global phase and should not make an otherwise sudden protocol appear inaccurate.

For an unbounded Hamiltonian, the full operator norm may be infinite. One then needs a state-specific estimate, an energy-restricted subspace, a quadratic-form argument, or a theorem adapted to the model. A useful first-order diagnostic for a fixed initial ket is

ϵi≡1ℏ∫titfΔiH(t) dt,\epsilon_i \equiv \frac{1}{\hbar} \int_{t_i}^{t_f} \Delta_i H(t)\,dt,

where

(ΔiH(t))2=⟨ψi∣H2(t)∣ψi⟩−⟨ψi∣H(t)∣ψi⟩2.\begin{aligned} \big(\Delta_iH(t)\big)^2 &= \langle\psi_i\rvert H^2(t)\lvert\psi_i\rangle \\ &\quad- \langle\psi_i\rvert H(t)\lvert\psi_i\rangle^2. \end{aligned}

If ϵi≪1\epsilon_i\ll1, the state changes only weakly beyond an overall phase at leading order. This is a diagnostic, not a universal error theorem for every unbounded or many-body problem.

The idealized protocol is

H(t)={Hi,t<0,Hf,t>0.H(t) = \begin{cases} H_i, & t\lt0,\\ H_f, & t\gt0. \end{cases}

Integrating the Schrödinger equation across a shrinking interval around t=0t=0 gives

∣ψ(0+)⟩−∣ψ(0−)⟩=−iℏlim⁡ε→0+×∫−εεH(t)∣ψ(t)⟩ dt.\begin{aligned} \lvert\psi(0^+)\rangle - \lvert\psi(0^-)\rangle &= -\frac{i}{\hbar} \lim_{\varepsilon\to0^+} \\ &\quad\times \int_{-\varepsilon}^{\varepsilon} H(t)\lvert\psi(t)\rangle\,dt. \end{aligned}

If the Hamiltonian has no singular impulse and the integral vanishes with ε\varepsilon, then

∣ψ(0+)⟩=∣ψ(0−)⟩.\lvert\psi(0^+)\rangle = \lvert\psi(0^-)\rangle.

Likewise, for a density operator,

ρ(0+)=ρ(0−).\rho(0^+)=\rho(0^-).

The time derivative generally changes:

iℏddt∣ψ(t)⟩∣0±=Hf/i∣ψ(0)⟩.i\hbar \left. \frac{d}{dt} \lvert\psi(t)\rangle \right|_{0^\pm} = H_{f/i} \lvert\psi(0)\rangle.

Thus the state is continuous while its subsequent trajectory changes abruptly.

If the quench changes a domain or even the Hilbert-space realization, “the same state” requires a specified embedding or matching map. Sudden changes of ideal hard-wall boundaries are useful models, but their domain assumptions should not be hidden.

Let the final Hamiltonian have spectral resolution

Hf=∑aEafPaf,H_f = \sum_a E_a^f P_a^f,

where PafP_a^f projects onto the full eigenspace at energy EafE_a^f. Immediately after an ideal sudden quench, the probability of finding that final energy is

paf=⟨ψi∣Paf∣ψi⟩.p_a^f = \langle\psi_i\rvert P_a^f \lvert\psi_i\rangle.

For a nondegenerate discrete spectrum,

Pmf=∣mf⟩⟨mf∣,P_m^f = \lvert m_f\rangle \langle m_f\rvert,

so

pmf=∣⟨mf∣ψi⟩∣2.p_m^f = \left\lvert \langle m_f\vert\psi_i\rangle \right\rvert^2.

If the initial state is an eigenstate ∣ni⟩\lvert n_i\rangle of HiH_i, then

pm←nsud=∣⟨mf∣ni⟩∣2.p_{m\leftarrow n}^{\mathrm{sud}} = \left\lvert \langle m_f\vert n_i\rangle \right\rvert^2.

These are overlap probabilities. No weak-coupling expansion is required, and no transition occurs during the zero-duration idealization. The word “transition” refers to the outcome obtained when the unchanged vector is analyzed in the new energy basis.

For a mixed initial state,

paf=Tr⁡(Pafρi).p_a^f = \operatorname{Tr} \left( P_a^f\rho_i \right).

Using projectors avoids basis-dependent statements inside a degenerate final eigenspace. For a continuum, the sum becomes a spectral measure or a probability density whose normalization must match the continuum eigenstates.

Control change, frozen state, and projection onto the final energy basis in the sudden approximation

An ideal sudden quench changes HiH_i to HfH_f over a time τq\tau_q too short for appreciable state evolution. The unchanged state is then resolved into final eigenstates, and each component acquires its own phase under the later evolution generated by HfH_f.

Write the initial vector in the final basis:

∣ψi⟩=∑mcm∣mf⟩,cm=⟨mf∣ψi⟩.\lvert\psi_i\rangle = \sum_m c_m\lvert m_f\rangle, \qquad c_m = \langle m_f\vert\psi_i\rangle.

For t>0t\gt0 and time-independent HfH_f,

∣ψ(t)⟩=∑mcme−iEmft/ℏ∣mf⟩.\lvert\psi(t)\rangle = \sum_m c_m e^{-iE_m^ft/\hbar} \lvert m_f\rangle.

The final energy probabilities ∣cm∣2\lvert c_m\rvert^2 are constant, but relative phases evolve. An observable AA therefore has expectation value

⟨A(t)⟩=∑m,ℓcm∗cℓei(Emf−Eℓf)t/ℏ×⟨mf∣A∣ℓf⟩.\begin{aligned} \langle A(t)\rangle &= \sum_{m,\ell} c_m^*c_\ell e^{i(E_m^f-E_\ell^f)t/\hbar} \\ &\qquad\times \langle m_f\rvert A\lvert\ell_f\rangle. \end{aligned}

Post-quench oscillations are not evidence that the sudden approximation failed. They are the exact evolution of the coherent superposition created by expressing the frozen state in the new basis.

For an ideal quench, the mean energy changes because the observable called the Hamiltonian changes:

ΔE‾=Tr⁡[ρi(Hf−Hi)].\Delta\overline E = \operatorname{Tr} \left[ \rho_i(H_f-H_i) \right].

For an initial eigenstate,

ΔE‾n=⟨ni∣Hf∣ni⟩−Eni.\Delta\overline E_n = \langle n_i\rvert H_f\lvert n_i\rangle - E_n^i.

The final energy variance is

(ΔEf)2=⟨Hf2⟩i−⟨Hf⟩i2.\big(\Delta E_f\big)^2 = \langle H_f^2\rangle_i - \langle H_f\rangle_i^2.

The mean change need not be positive for an arbitrary prepared state and protocol. Calling it mean work also requires a stated work convention. For an initial state diagonal in HiH_i, the two-point measurement mean agrees with this energy change, while the full distribution contains the probabilities for individual energy differences. That measurement framework belongs in Work Distributions.

Warm-Up: A Sudden Rotation of a Two-Level Hamiltonian

Section titled “Warm-Up: A Sudden Rotation of a Two-Level Hamiltonian”

Consider

Hi=−Δ2σzH_i = -\frac{\Delta}{2}\sigma_z

and

Hf=−Δ2cos⁡θ σz−Δ2sin⁡θ σx,Δ>0.\begin{aligned} H_f &= -\frac{\Delta}{2} \cos\theta\,\sigma_z \\ &\quad- \frac{\Delta}{2} \sin\theta\,\sigma_x , \\ \Delta &\gt 0. \end{aligned}

Prepare the ground state ∣gi⟩\lvert g_i\rangle of HiH_i. The final ground-state Bloch direction is rotated by angle θ\theta, so

∣⟨gf∣gi⟩∣2=cos⁡2θ2,\left\lvert \langle g_f\vert g_i\rangle \right\rvert^2 = \cos^2\frac{\theta}{2},

and

∣⟨ef∣gi⟩∣2=sin⁡2θ2.\left\lvert \langle e_f\vert g_i\rangle \right\rvert^2 = \sin^2\frac{\theta}{2}.

The mean injected energy is

ΔE‾=Δsin⁡2θ2.\Delta\overline E = \Delta \sin^2\frac{\theta}{2}.

For θ=π\theta=\pi, the old ground state is the new excited state. For θ=0\theta=0, the Hamiltonian basis is unchanged and the quench creates no excitation even if a common energy offset changes.

Example: Sudden Expansion of an Infinite Well

Section titled “Example: Sudden Expansion of an Infinite Well”

Take an initial infinite well on

0<x<Li0\lt x\lt L_i

and suddenly move the right wall to

Lf=αLi,α>1.L_f = \alpha L_i, \qquad \alpha\gt1.

The initial eigenfunction is

ϕni(x)=2Lisin⁡(nπxLi)\phi_n^i(x) = \sqrt{\frac{2}{L_i}} \sin\left( \frac{n\pi x}{L_i} \right)

for 0<x<Li0\lt x\lt L_i. Immediately after the expansion, the wavefunction has the same values on the old interval and is zero on the newly opened region:

ψ(x,0+)={ϕni(x),0<x<Li,0,Li<x<Lf.\psi(x,0^+) = \begin{cases} \phi_n^i(x), & 0\lt x\lt L_i,\\ 0, & L_i\lt x\lt L_f. \end{cases}

The final eigenfunctions are

ϕmf(x)=2Lfsin⁡(mπxLf).\phi_m^f(x) = \sqrt{\frac{2}{L_f}} \sin\left( \frac{m\pi x}{L_f} \right).

Their overlap amplitudes are

cm=⟨mf∣ni⟩=2α∫01dy ×sin⁡(mπyα)sin⁡(nπy)=1αsin⁡ ⁣[π(m/α−n)]π(m/α−n)−1αsin⁡ ⁣[π(m/α+n)]π(m/α+n).\begin{aligned} c_m &= \langle m_f\vert n_i\rangle \\ &= \frac{2}{\sqrt{\alpha}} \int_0^1dy\, \\ &\quad\times \sin\left( \frac{m\pi y}{\alpha} \right) \sin(n\pi y) \\ &= \frac{1}{\sqrt{\alpha}} \frac{ \sin\!\left[\pi(m/\alpha-n)\right] }{ \pi(m/\alpha-n) } \\ &\quad- \frac{1}{\sqrt{\alpha}} \frac{ \sin\!\left[\pi(m/\alpha+n)\right] }{ \pi(m/\alpha+n) } . \end{aligned}

At a removable point m=αnm=\alpha n, the limiting value is

cm=1α.c_m = \frac{1}{\sqrt{\alpha}}.

The final probabilities are

Pm=∣cm∣2,∑m=1∞Pm=1.P_m = \lvert c_m\rvert^2, \qquad \sum_{m=1}^{\infty}P_m=1.

For a doubling of the well, α=2\alpha=2, from the initial ground state n=1n=1,

P1=329π2≃0.360,P2=12,P3=3225π2≃0.130.\begin{aligned} P_1 &= \frac{32}{9\pi^2} \simeq 0.360, \\ P_2 &= \frac12, \\ P_3 &= \frac{32}{25\pi^2} \simeq 0.130. \end{aligned}

All other even-mm amplitudes vanish, while higher odd levels carry the remaining probability. The final energies are

Emf=π2ℏ2m22Mα2Li2.E_m^f = \frac{\pi^2\hbar^2m^2}{ 2M\alpha^2L_i^2 }.

For this ideal expansion, the quadratic-form expectation of the kinetic energy is unchanged at the instant of the quench:

∑mPmEmf=Eni.\sum_mP_mE_m^f = E_n^i.

The energy distribution broadens even though its mean is unchanged. Subsequent interference among the final eigenstates drives nonstationary motion inside the larger well. The stationary spectrum and eigenfunctions themselves belong in Infinite Square Well.

Example: Sudden Change of Oscillator Frequency

Section titled “Example: Sudden Change of Oscillator Frequency”

Let

H(ω)=p22M+12Mω2x2,H(\omega) = \frac{p^2}{2M} + \frac12M\omega^2x^2,

and quench

ωi⟶ωf.\omega_i \longrightarrow \omega_f.

Prepare the initial ground state ∣0i⟩\lvert0_i\rangle. It remains the same Gaussian immediately after the switch, but its width does not generally match the final ground-state width. Define

r≡12ln⁡(ωfωi).r \equiv \frac12 \ln\left( \frac{\omega_f}{\omega_i} \right).

In the final oscillator basis, the initial ground state is a squeezed vacuum up to the sign convention used for the squeeze operator. Only even final number states occur:

P2k+1=0,P_{2k+1}=0,

and

P2k=(2k)!22k(k!)2tanh⁡2krcosh⁡r.P_{2k} = \frac{(2k)!}{ 2^{2k}(k!)^2 } \frac{\tanh^{2k}r}{\cosh r}.

The final ground-state survival probability is

P0=1cosh⁡r=2ωiωfωi+ωf.\begin{aligned} P_0 &= \frac{1}{\cosh r} \\ &= \frac{ 2\sqrt{\omega_i\omega_f} }{ \omega_i+\omega_f }. \end{aligned}

The mean final occupation is

⟨Nf⟩=sinh⁡2r=14(ωfωi+ωiωf−2).\langle N_f\rangle = \sinh^2r = \frac14 \left( \frac{\omega_f}{\omega_i} + \frac{\omega_i}{\omega_f} - 2 \right).

Consequently,

⟨Hf⟩=ℏωf(sinh⁡2r+12)=ℏ4(ωi+ωf2ωi).\begin{aligned} \langle H_f\rangle &= \hbar\omega_f \left( \sinh^2r+\frac12 \right) \\ &= \frac{\hbar}{4} \left( \omega_i + \frac{\omega_f^2}{\omega_i} \right). \end{aligned}

Parity, not energy level number, is preserved by the basis overlap because both oscillator Hamiltonians are even under x↦−xx\mapsto-x. The general geometry and number-state expansion of squeezed states belong in Squeezed States: First Encounter.

For a short but nonzero quench, remove the scalar phase and expand the residual evolution:

eiΦqUq=I−iℏ∫titfK(t) dt+O ⁣(Aq2).\begin{aligned} e^{i\Phi_q}U_q &= I - \frac{i}{\hbar} \int_{t_i}^{t_f}K(t)\,dt \\ &\quad+ O\!\left( \mathcal A_q^2 \right). \end{aligned}

For an initial eigenstate ∣ni⟩\lvert n_i\rangle, the final-basis amplitude is therefore

eiΦqcm(τq)=⟨mf∣ni⟩−iℏ∫titf⟨mf∣K(t)∣ni⟩ dt+O ⁣(Aq2).\begin{aligned} e^{i\Phi_q}c_m(\tau_q) &= \langle m_f\vert n_i\rangle \\ &\quad- \frac{i}{\hbar} \int_{t_i}^{t_f} \langle m_f\rvert K(t) \lvert n_i\rangle\,dt \\ &\quad+ O\!\left( \mathcal A_q^2 \right). \end{aligned}

The first term is the ideal sudden overlap. The second is actual evolution during the ramp. These amplitudes interfere; one should not add the corresponding probabilities as if they described mutually exclusive processes.

If the first correction is too large, solve the short protocol with the exact time-evolution operator, a Magnus expansion, or direct numerical propagation. The Dyson Expansion for Transition Amplitudes gives the systematic ordered expansion.

A Short Pulse Is Not Always a Frozen Quench

Section titled “A Short Pulse Is Not Always a Frozen Quench”

Suppose the Hamiltonian contains a singular pulse,

H(t)=Hreg(t)+G δ(t−t0),H(t) = H_{\mathrm{reg}}(t) + G\,\delta(t-t_0),

with self-adjoint impulse operator GG. Even though the pulse duration is idealized as zero, its integrated action is finite:

∫dt G δ(t−t0)=G.\int dt\, G\,\delta(t-t_0) = G.

The state changes discontinuously according to

∣ψ(t0+)⟩=e−iG/ℏ∣ψ(t0−)⟩.\lvert\psi(t_0^+)\rangle = e^{-iG/\hbar} \lvert\psi(t_0^-)\rangle.

This is an impulsive unitary kick, not the frozen-state sudden approximation. The relevant small parameter is pulse area divided by ℏ\hbar, not duration alone. A family of ever shorter pulses whose amplitude grows as 1/τq1/\tau_q can retain a nontrivial unitary limit.

The sudden and adiabatic approximations answer opposite limiting questions.

FeatureSudden limitAdiabatic limit
Protocol scaleshort compared with relevant dynamical timeslong compared with inverse-gap transition scales
What approximately followsthe fixed state vectorthe corresponding instantaneous eigenspace
Final probabilitiesoverlaps with the final spectral projectorspopulation remains in the followed eigenspace, up to error
Phase during protocolusually negligible except for a common phasedynamical and geometric phases matter
Main dangerfinite integrated action or unresolved high-energy modessmall gaps, crossings, and accumulated nonadiabatic coupling

For a large rotation of an eigenbasis, the two limits make sharply different predictions. In the sudden limit the old eigenvector is projected onto many final eigenvectors. In the adiabatic limit an initially isolated eigenstate follows its corresponding instantaneous eigenspace. Runtime diagnostics and leakage estimates are developed in Adiabatic Approximation as a Method, while the theorem-level statement is summarized in the Adiabatic Theorem.

Intermediate protocols are not generally described by interpolating probabilities between the two limits. They require solving the actual dynamics. Landau–Zener Transition is the canonical two-level example.

If EafE_a^f is degenerate, only

paf=Tr⁡(Pafρi)p_a^f = \operatorname{Tr} \left( P_a^f\rho_i \right)

is basis independent. Probabilities assigned to arbitrary basis vectors inside the eigenspace are meaningful only when an additional commuting observable or later measurement resolves them.

For generalized final eigenstates ∣E,α⟩\lvert E,\alpha\rangle,

pα(E)=∣⟨E,α∣ψi⟩∣2p_\alpha(E) = \left\lvert \langle E,\alpha\vert\psi_i\rangle \right\rvert^2

is a density with respect to the completeness measure. Its dimensions depend on whether the states are normalized in energy, momentum, or a finite box. The normalization dictionary in Density of States in Transition Rates applies here as well, although a sudden projection is not itself a golden-rule rate.

For an extensive Hamiltonian, a full-system operator norm and global fidelity can scale poorly with volume. A protocol may fail to be globally sudden even when local observables barely change during the switch. High-energy sectors, ultraviolet regularization, locality, and the order of thermodynamic and zero-duration limits can matter. The elementary sudden approximation should not be promoted to a uniform many-body theorem without specifying those controls.

  1. Specify the protocol. State HiH_i, HfH_f, the switching function, and the duration τq\tau_q.
  2. Identify the relevant state space. Name the prepared state, occupied energy window, or controlled subspace.
  3. Estimate integrated action. Remove scalar phases and test Aq\mathcal A_q, a state-specific energy spread, or model-specific transition amplitudes.
  4. Check for an impulse. A finite pulse area can produce a nontrivial unitary even as the duration vanishes.
  5. Match the Hilbert spaces. If boundaries or domains change, state the embedding of the pre-quench wavefunction into the post-quench problem.
  6. Project with spectral projectors. Use pa=Tr⁡(Pafρi)p_a=\operatorname{Tr}(P_a^f\rho_i), especially in degenerate sectors.
  7. Separate the switch from later evolution. The overlap distribution is set at the quench; later phases generate observable dynamics.
  8. Check normalization and energy moments. Verify ∑apa=1\sum_ap_a=1 or the corresponding continuum integral.
  9. Estimate corrections. Compare the finite-ramp amplitude with the ideal overlap before trusting the sudden result quantitatively.
  • Saying only that the change is “fast.” Fast must be compared with matrix elements, gaps, and the integrated nontrivial action.
  • Assuming a large Hamiltonian change invalidates suddenness. A large endpoint difference can still be sudden if the switching interval is sufficiently short on the relevant subspace.
  • Assuming every zero-duration protocol leaves the state fixed. A delta-like impulse with finite area applies a unitary kick.
  • Projecting before matching Hilbert spaces. Boundary quenches require a clear embedding of the initial state into the final domain.
  • Using basis-vector probabilities in a degenerate eigenspace. The invariant object is the probability of the whole spectral projector.
  • Calling post-quench oscillations transitions during the quench. They come from phase evolution under HfH_f after the overlap distribution has been created.
  • Confusing overlap probability with a golden-rule rate. A sudden quench can be large and nonperturbative; no long-time continuum limit is involved.
  • Ignoring normalization conventions in a continuum. Probability densities change when the generalized eigenstate normalization changes.
  • Treating an unbounded-spectrum time-scale slogan as a proof. State domains and high-energy tails require model-specific control.
  • Assuming sudden and adiabatic probabilities interpolate linearly. Intermediate dynamics retains phases and interference.

A two-level Hamiltonian is suddenly rotated through an angle θ\theta as in the warm-up example. The system starts in the initial ground state. Find the final energy probabilities and the final energy variance.

Solution

The probabilities are

pg=cos⁡2θ2,pe=sin⁡2θ2.p_g = \cos^2\frac{\theta}{2}, \qquad p_e = \sin^2\frac{\theta}{2}.

The final energies are −Δ/2-\Delta/2 and +Δ/2+\Delta/2. Hence

⟨Hf⟩=−Δ2cos⁡θ.\langle H_f\rangle = -\frac{\Delta}{2}\cos\theta.

Since Hf2=(Δ2/4)IH_f^2=(\Delta^2/4)I,

(ΔEf)2=Δ24−Δ24cos⁡2θ=Δ24sin⁡2θ.\begin{aligned} \big(\Delta E_f\big)^2 &= \frac{\Delta^2}{4} - \frac{\Delta^2}{4}\cos^2\theta \\ &= \frac{\Delta^2}{4} \sin^2\theta. \end{aligned}

The variance vanishes for θ=0\theta=0 and θ=π\theta=\pi because the old ground state is then a definite final energy state.

For a sudden expansion from Li=LL_i=L to Lf=2LL_f=2L, starting from n=1n=1, calculate P1P_1, P2P_2, and P3P_3. Why is P2P_2 obtained from a limit?

Solution

Substituting α=2\alpha=2 and n=1n=1 into the overlap formula gives

c1=83π2,c3=85π2.c_1 = \frac{8}{3\pi\sqrt2}, \qquad c_3 = \frac{8}{5\pi\sqrt2}.

Therefore

P1=329π2,P3=3225π2.P_1 = \frac{32}{9\pi^2}, \qquad P_3 = \frac{32}{25\pi^2}.

For m=2m=2, the closed form contains a removable 0/00/0 term because m/α=nm/\alpha=n. Returning to the integral gives

c2=22∫01sin⁡2(πy) dy=12,c_2 = \frac{2}{\sqrt2} \int_0^1\sin^2(\pi y)\,dy = \frac{1}{\sqrt2},

so

P2=12.P_2=\frac12.

The three probabilities sum to approximately 0.9900.990; higher odd final levels contain the remainder.

An oscillator frequency is suddenly changed from ωi\omega_i to 4ωi4\omega_i. Find the final ground-state probability and the mean final excitation number.

Solution

Here

r=12ln⁡4=ln⁡2.r = \frac12\ln4 = \ln2.

The ground-state probability is

P0=2ωi(4ωi)ωi+4ωi=45.\begin{aligned} P_0 &= \frac{ 2\sqrt{\omega_i(4\omega_i)} }{ \omega_i+4\omega_i } \\ &= \frac45. \end{aligned}

The mean occupation is

⟨Nf⟩=14(4+14−2)=916.\begin{aligned} \langle N_f\rangle &= \frac14 \left( 4+\frac14-2 \right) \\ &= \frac{9}{16}. \end{aligned}

Although the final ground-state probability is large, the remaining probability occupies an infinite set of even levels and gives a nonzero mean excitation.

On a controlled two-level subspace, suppose

∥K(t)∥≤ℏΩ2\lVert K(t)\rVert \le \frac{\hbar\Omega}{2}

throughout a quench of duration τq\tau_q. Give a sufficient suddenness condition and bound the deviation from identity up to a scalar phase.

Solution

The integrated action obeys

Aq=1ℏ∫titf∥K(t)∥ dt≤Ωτq2.\begin{aligned} \mathcal A_q &= \frac{1}{\hbar} \int_{t_i}^{t_f} \lVert K(t)\rVert\,dt \\ &\le \frac{\Omega\tau_q}{2}. \end{aligned}

Therefore

∥Uq−e−iΦqI∥≤Ωτq2.\left\lVert U_q-e^{-i\Phi_q}I \right\rVert \le \frac{\Omega\tau_q}{2}.

A sufficient suddenness condition is

Ωτq≪1.\Omega\tau_q\ll1.

This bound is conservative; a particular state can experience smaller corrections because of symmetry or cancellations.

Let

Hτ(t)=1τg(t/τ)G,H_\tau(t) = \frac{1}{\tau} g(t/\tau)G,

where g(s)g(s) has support on 0≤s≤10\le s\le1 and

∫01g(s) ds=1.\int_0^1g(s)\,ds=1.

Find the τ→0\tau\to0 evolution across the pulse when GG is time independent.

Solution

The integrated Hamiltonian is independent of τ\tau:

∫0τHτ(t) dt=∫0τ1τg(t/τ)G dt=G.\begin{aligned} \int_0^\tau H_\tau(t)\,dt &= \int_0^\tau \frac{1}{\tau} g(t/\tau)G\,dt \\ &= G. \end{aligned}

All pulse Hamiltonians commute because they are proportional to the same GG, so

Uτ=e−iG/ℏ.U_\tau = e^{-iG/\hbar}.

The limit is nontrivial rather than the identity. Shrinking duration does not imply a frozen state when the amplitude grows so that the pulse area remains finite.

The final Hamiltonian has a two-dimensional eigenspace with orthonormal basis ∣a⟩,∣b⟩\lvert a\rangle,\lvert b\rangle. Explain which sudden-quench probability is basis independent.

Solution

The spectral projector is

P=∣a⟩⟨a∣+∣b⟩⟨b∣.P = \lvert a\rangle\langle a\rvert + \lvert b\rangle\langle b\rvert.

For an initial density operator ρi\rho_i, the probability of the degenerate energy is

p=Tr⁡(Pρi).p = \operatorname{Tr}(P\rho_i).

Rotating ∣a⟩,∣b⟩\lvert a\rangle,\lvert b\rangle within the eigenspace changes the two individual basis probabilities but leaves their sum pp invariant. Individual outcomes require an additional observable that resolves a preferred basis inside the degenerate space.

  • L. D. Landau and E. M. Lifshitz, Quantum Mechanics: Non-Relativistic Theory, 3rd ed., Pergamon Press, 1977.
  • A. Messiah, Quantum Mechanics, Vol. II, North-Holland, 1962.
  • 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.
  • D. J. Griffiths and D. F. Schroeter, Introduction to Quantum Mechanics, 3rd ed., Cambridge University Press, 2018.
  • A. Polkovnikov, K. Sengupta, A. Silva, and M. Vengalattore, “Colloquium: Nonequilibrium dynamics of closed interacting quantum systems,” Reviews of Modern Physics 83, 863–883 (2011).