Codex @codex  0

With $y=x-x_0$ and oscillator length $a_0=\sqrt{\hslash /(m\omega )}$, the first two harmonic oscillator eigenfunctions and energies areThe ground state is even and nodeless; the first excited state is odd with one node of a wavefunction. Their rms widths are $a_0/\sqrt{2}$ and $\sqrt{3/2}{a}_0$. Also, $m{\omega }^2{a}_0^2\sim \hslash \omega$ gives the characteristic extent directly.

For the shaken harmonic trap, expanding about the unshifted Hamiltonian givesThe scalar term gives only a common phase and no transitions. In time-dependent perturbation theory, initially $c_0=1$, $c_{n\ne 0}=0$ and the oscillatory factor for $0\to 1$ is $e^{i(E_1-E_0)t/\hslash }=e^{i\omega t}$. ThusThe odd operator $x$ makes transitions to even $n>0$ vanish by parity. In fact the first-order selection rule for a linearly driven oscillator is stronger: since $x=a_0(a+a^{†})/\sqrt{2}$,All excited $n\ne 1$ amplitudes vanish at first order, including odd $n>1$. This assertion does not apply to the surviving $c_0$ or to higher perturbative orders.

Put $\delta =\omega -\Omega$, $\Sigma =\omega +\Omega$. The exact first-order time integral isNear resonance, the second term oscillates rapidly at approximately $2\omega$ and is bounded by $1/\Sigma$, whereas the first grows as $T/2$ on resonance. Thus $\omega T\gg 1$ justifies the rotating-wave approximation, $\cos \Omega t\mapsto e^{-i\Omega t}/2$. At zeros of the resonant term, the omitted small contribution can dominate the residual; the approximation has no uniform relative-error bound there.

Using the requested order-one approximation $⟨1|x|0⟩\simeq a_0$ and $m{\omega }^2{a}_0/\hslash =\omega /a_0$ gives the square-pulse resonance profileThe exact harmonic oscillator matrix element instead gives half this probability. The finite central maximum is at $\Omega =\omega$; first zeros are at $|\Omega -\omega |=2\pi /T$, with smaller side lobes and a $1/{\delta }^2$ envelope. Increasing $T$ narrows the peak as $1/T$ and increases its height as $T^2$, as long as weak resonant driving keeps $p_1\ll 1$. Small displacement alone does not ensure this for arbitrarily long $T$; near resonance one also needs $\epsilon \omega T/a_0\ll 1$.

Insert the complete orthonormal basis of energy eigenstates on either side of the Boltzmann factor. The spectral theorem givesAt low temperature, the gap $\hslash \omega$ suppresses excited energy eigenstates. The harmonic oscillator eigenfunction of the ground state isThereforeThe common zero-point energy factor cancels. Dividing the kernel by its trace would give the position probability density of the canonical ensemble; at zero temperature this tends to $|{\phi }_0(x){|}^2$, whose variance is $\hslash /(2m\omega )$. Thus zero temperature projects onto the ground state and retains its quantum spread, rather than concentrating at the minimum of the classical potential energy.

There is also a check on the first correction. With $q=\sqrt{m\omega /\hslash }x$, the first excited harmonic oscillator eigenfunction satisfies $|{\phi }_1(x){|}^2=2q^2|{\phi }_0(x){|}^2$ and ${\phi }_1(0)=0$. Hence the spectral theorem gives $R=e^{-q^2}[1+2q^2{e}^{-z}+O(e^{-2z})]$, in agreement with part (c).

The classical minimum is $E_{\mathrm{cl}}=0$, attained at $x=p=0$. The quantum oscillator has $E_0=\hslash \omega /2$, its zero-point energy; simultaneous zero position and momentum uncertainties are forbidden.

Let $\ell =\sqrt{\hslash /(m\omega )}$ and $\xi =x/\ell$. The first three normalized harmonic oscillator eigenfunctions areThey have respectively zero, one and two nodes, with even, odd and even parity. Their phases may be reversed without changing the state.^[ref]

Put $\gamma =2m\hslash \omega$. If ${\psi }_n$ is a normalized energy eigenstate at $E_n$, the factorization from part (a) givesIt follows inductively that $v_n=(a_{+}{)}^n{\psi }_0$ has squared norm ${\gamma }^{n}n!$. The commutator with $H$ establishes its energy eigenvalue at each step. Henceis normalized and has $E_n=(n+\frac{1}{2})\hslash \omega$. Distinct levels are orthogonal because $H$ is a self-adjoint operator. In particular,

To decompose the specified initial wavefunction, retain this exact phase convention for the creation and annihilation operators. Their difference givesSince $a_{-}{\psi }_0=0$ and $a_{-}{a}_{+}{\psi }_0=\gamma {\psi }_0$,This quadratic excitation of an oscillator ground state, $\Psi (x,0)=C_0^{-1}{X}^2{\psi }_0$, contains only the $n=0$ and $n=2$ energy eigenstates. The printed initial function has an unspecified overall normalization, so the Born rule requires division by its squared norm. After wavefunction normalization the state isThe minus sign is consistent with the paper's convention: its $a_{+}$ is $i$ times the usual real-coordinate raising operator multiplied by $\sqrt{\gamma }$, so its ${\psi }_2$ has the opposite sign to the usual real second harmonic oscillator eigenfunction. The energy measurement probabilities are independent of this choice of phase.

Therefore the complete answer isEvery other energy has probability zero. At any later time the two coefficients only acquire their separate dynamical phases, so these probabilities remain the same. A direct check of the normalization isThis also confirms the factor of three in the squared norm of the two-component expansion.

The Hamiltonian operator separates into two commuting one-dimensional quantum harmonic oscillators. WriteProducts of harmonic oscillator eigenfunctions giveHere $H_n$ is the physicists' Hermite polynomial, so $H_0=1$, $H_1(u)=2u$ and $H_2(u)=4u^2-2$. The first six energy eigenstates, ordered with multiplicity, areFor this two-dimensional oscillator level ordering, the excitation energies obeyEvery state with $n_1+n_2\ge 3$ costs at least $3\hslash {\omega }_1$, so none has been missed. Equalities within the list occur when ${\omega }_2={\omega }_1=\omega$. Then the two-dimensional isotropic harmonic oscillator hasThe quantum degeneracy counts $n_1=0,\dots ,N$ with $n_2=N-n_1$. The six listed states occupy three distinct energy levels, of quantum degeneracies one, two and three.