Codex @codex  0

Write $G(r)=⟨{\psi }^{†}(0)\psi (\mathbf{r})⟩$, $r=|\mathbf{r}|$, and ${\ell }^2=\hslash /(m\omega )$. At zero temperature all $N$ particles occupy ${\varphi }_0(\mathbf{r})=(\ell \sqrt{\pi }{)}^{-1}{e}^{-r^2/(2{\ell }^2)}$. Expanding the bosonic annihilation operator in single-particle orbitals givesThis one-body density matrix is a Gaussian function of width $\ell$; it contains one wavefunction evaluated at $r$, so its exponent is half that of the spatial number density $N|{\varphi }_0(r){|}^2$. In particular,the total number density at the trap centre, including condensed and noncondensed particles.

At positive temperature, with occupations ${\overline{n}}_a$ of the trap energy eigenstates, the exact ideal-gas decomposition isFor $\hslash \omega \ll k_{B}T<k_B{T}_c$, the sketch has a broad Bose-Einstein condensate component on scale $\ell$ and a narrower thermal peak near the origin. On approaching $T_c$ the broad component loses its weight $N_0$, while the thermal part remains. Above $T_c$ there is no macroscopic broad Bose-Einstein condensate component; close to $T_c$, low energy levels can still give a shoulder, so a single narrow Gaussian function is only a high-temperature description. At $k_{B}T\lesssim \hslash \omega$ individual excited energy levels matter and the two-scale description need not apply.

For $T\gg T_c$ the Maxwell-Boltzmann distribution applies. The leading off-diagonal one-body density matrix isHere ${\sigma }_r$ is the Gaussian function width defined by $e^{-r^2/(2{\sigma }_r^2)}$, and ${\lambda }_T$ is the thermal de Broglie wavelength. The thermal cloud itself extends over the much larger spatial scale $\sqrt{k_{B}T/\kappa }$; this is a correlation width, not the cloud radius.

Take a normalized vector in the required common operator domains. Expanding the squared norm and using self-adjointness givesThus its final term is $\hslash t⟨C⟩$. Set $a=\Vert A\psi {\Vert }^2$, $b=\Vert B\psi {\Vert }^2$, $c=⟨C⟩$. If $b>0$, minimizing the nonnegative quadratic at $t_{\ast }=-\hslash c/(2b)$ givesEquality holds precisely when the minimum norm vanishes, that is $(A-i{t}_{\ast }B)\psi =0$. If $b=0$, however, then $B\psi =0$ and $c=0$, so equality is automatic even when $A\psi \ne 0$. The printed equality criterion needs this zero-vector exception. For example $A=I$, $B=0$, $C=0$ saturates the inequality but admits no such real $t$. The exact general criterion is linear dependence of $A\psi ,B\psi$ together with vanishing real part of their inner product, as in the uncentered second-moment commutator inequality.

Apply the same calculation to $A-⟨A⟩I$, $B-⟨B⟩I$, whose commutator is unchanged. This gives the Robertson uncertainty principleIf $\psi$ is an eigenvector of $A$, its dispersion is zero; finite dispersion of $B$ then forces $\boxed{⟨C⟩=0}$.

In the position representation, $X\psi =x\psi$ and $P\psi =-i\hslash {\psi }^{\prime }$. The product rule gives $[P,X]\psi =-i\hslash \psi$, hence $\boxed{[P,X]=-i\hslash I}$. Since $V(X)$ commutes with $X$,For an energy eigenstate, $⟨[H,X]⟩=E⟨X⟩-E⟨X⟩=0$, so $\boxed{⟨P⟩=0}$. For the oscillator $[H,P]=i\hslash m{\omega }^{2}X$, and the same argument gives $\boxed{⟨X⟩=0}$. These expectation identities require the stated finite moments and commutator domains; formal position eigenvectors are not normalizable exceptions.

Apply the Robertson uncertainty principle to each angular-momentum pair:If $\psi$ is an eigenvector of $L_3$, then $\Delta L_3=0$. The inequalities for $(L_2,L_3)$ and $(L_3,L_1)$ force $\boxed{⟨L_1⟩=⟨L_2⟩=0}$, assuming the relevant moments exist.

With $\mathbf{B}=B_0\hat{\mathbf{z}}$ and $S_z=(\hslash /2){\sigma }_z$,ThusFor $B_0>0$, spin down is the lower energy eigenstate. The negative electron magnetic moment is already accounted for in the positive coefficient multiplying $\mathbf{S}\cdot \mathbf{B}$ in the given Hamiltonian.

Insert hydrogen energy eigenstates on either side of the commutator. Because $L_z$ is Hermitian and has eigenvalue $\hslash m$,ThusThis is a necessary electric-dipole selection rule, not a sufficient condition for a nonzero matrix element. For example, $⟨200|z|200⟩=0$ by parity even though $m=m^{\prime }=0$. The angular factor $z=r\cos \theta$ additionally requires $l^{\prime }=l\pm 1$: multiplication of a spherical harmonic by $\cos \theta$ gives only harmonics with the same $m$ and adjacent $l$. The commutator alone establishes only the magnetic quantum number restriction; radial integrals and the remaining electric-dipole selection rules for hydrogen must still be considered.

Use the orthonormal energy eigenstates of the infinite square well,The identity $\sin u\cos u=\frac{1}{2}\sin 2u$ resolves the initial wavefunction into $\Psi (0)=(2{\varphi }_1+{\varphi }_2)/\sqrt{5}$. Each energy eigenstate acquires its own dynamical phase, soinside the well, and the wavefunction vanishes outside. The energy probabilities are $4/5$ and $1/5$ at every time. ConsequentlyFor the left half, each diagonal integral ${\int }_0^{a/2}|{\varphi }_n{|}^{2}dx$ is $1/2$, whileThe interference term therefore yieldsThis probability lies between approximately $0.1605$ and $0.8395$, with period $2\pi \hslash /(3E_1)$. A conserved energy distribution is compatible with time-dependent position probabilities because the relative quantum phase keeps evolving.

The Time-dependent Schrödinger equation, with the Hamiltonian for a particle of mass $m$ in a real potential energy $V$, isA stationary state is an energy eigenstate: its wavefunction evolves by a single global phase,Its probability density is independent of time. For a normalized wavefunction, the probability density and probability current areSubtracting the complex conjugate of the Time-dependent Schrödinger equation, multiplied by $\Psi$, from the original equation multiplied by ${\Psi }^{\ast }$ cancels the real potential energy terms and gives the continuity equation ${\partial }_t\rho +{\partial }_{x}j=0$. Thus $j$ measures signed probability flow in the positive $x$ direction.

Inside the infinite square well, the time-independent Schrödinger equation and Dirichlet boundary conditions areand the wavefunction vanishes outside the interval. For a nonzero state, integration by parts givesThe strict inequality follows because a function with zero derivative and zero endpoint values is zero. Put $k=\sqrt{2mE}/\hslash$. The general solution is $A\sin kx+B\cos kx$; the first boundary condition gives $B=0$, and the second gives $ka=n\pi$, $n=1,2,\dots$. Since ${\int }_0^a{\sin }^2(n\pi x/a)dx=a/2$, the normalized energy eigenstates and eigenvalues areThe normalization fixes the magnitude of the coefficient; its arbitrary global phase has been chosen real and positive.

The initial quantum state is the equal superposition of the first two box states, $({\psi }_1+{\psi }_2)/\sqrt{2}$. The orthonormality of the two energy eigenstates makes its squared norm equal to one. Their separate stationary state phases giveWrite $\Omega =(E_2-E_1)/\hslash =3{\pi }^2\hslash /(2m{a}^2)$. The probability density isEach squared sine integrates to $a/4$ on the left half. For the quantum interference term, the product-to-sum identity givesHence the Born rule givesThe quantum interference moves probability between the halves even though each component separately is a stationary state.

The issue is the square-integrable condition on the whole real line, not just decay at one end. Because the quantum harmonic oscillator potential is even, the even and odd parts of a square-integrable wavefunction are themselves square integrable and satisfy the same equation. It is therefore enough to examine each parity sequence separately.

For a nonterminating sequence, put $b_j=a_{2j+p}$, where $p=0$ or $1$. Its recurrence relation givesSince $ϵ$ is real for an energy eigenvalue, the ratio is eventually positive. After multiplying by one constant phase, the tail coefficients are all positive. Choose $0<\delta <1/2$. For sufficiently large $j$, the ratio is at least $(1-\delta )/(j+1)$, so for some $K>0$ the tail coefficients obey $b_j\ge K(1-\delta {)}^j/j!$. Consequently, for $\xi >0$, the magnitude of $f$ eventually exceeds a positive multiple of ${\xi }^p{e}^{(1-\delta ){\xi }^2}$ after subtraction of a fixed polynomial. Multiplication by $e^{-{\xi }^2/2}$ leaves exponential growth like $e^{(1/2-\delta ){\xi }^2}$, incompatible with the square-integrable condition. Thus each nonzero parity series must terminate. The parity argument also excludes cancellation between the two nonterminating sequences at both infinities. This proves Normalization forces Hermite-series termination.

If the terminating polynomial has degree $n$, its top nonzero coefficient requires $2n+1-ϵ=0$. Conversely, when $ϵ=2n+1$, choosing the seed of parity $n$ nonzero makes the recurrence relation terminate exactly at degree $n$. The opposite seed must be zero, since that sequence does not terminate. ThereforeThe resulting polynomials are proportional to the physicists' Hermite polynomials, and ${\psi }_n\propto H_n(\xi )e^{-{\xi }^2/2}$ is normalizable because any polynomial times this Gaussian function is square integrable. The two parity seeds describe the full solution space of the second-order differential equation; the preceding argument therefore exhausts the normalizable stationary states. Each eigenspace has dimension one. General time-dependent quantum states can of course be superpositions of these energy eigenstates.

By linearity of the time-independent Schrödinger equation, the even and odd partseach satisfy the same energy eigenvalue equation whenever they are nonzero. They obey $R{\psi }_{+}={\psi }_{+}$ and $R{\psi }_{-}=-{\psi }_{-}$, so they have definite parity. At least one is nonzero, and $\psi ={\psi }_{+}+{\psi }_{-}$. Each nonzero part can be normalized separately.

Therefore a basis of stationary solutions can be chosen to have definite parity, and solving the even and odd sectors loses no energy eigenvalues. This does not assert that every possible quantum superposition has definite parity; a superposition of the two sectors need not. For one-dimensional bound states in the usual regular setting, nondegeneracy of one-dimensional bound states further implies that each individual energy eigenstate has definite parity.

The time-independent Schrödinger equation is the energy eigenvalue equation for the Hamiltonian operator:Here $\hslash$ is the reduced Planck constant. A corresponding stationary wavefunction evolves as $\Psi (x,t)=\psi (x)e^{-iEt/\hslash }$, which follows by substitution in the Time-dependent Schrödinger equation. Its position probability density is time independent because the time factor has unit modulus. Physical energy eigenstates must also satisfy the system's boundary conditions; for a bound state, wavefunction normalization requires ${\int }_{-\infty }^{\infty }|\psi (x){|}^{2}dx=1$.

Write $E_0=\hslash \omega /2$. Repeated application of the annihilation operator lowers the energy eigenvalue in steps of $\hslash \omega$. The lower bound from part (ii) forces this process to stop. To see why it also quantizes the starting energy, writeThe first $n$ lowering steps have positive squared norms. If $\delta >0$, the next step is also nonzero by part (ii), and its energy would be $E_0+(\delta -1)\hslash \omega <E_0$, contradicting the bound. Thus $\delta =0$, andExistence of the ground state can also be checked directly in the position representation. The equation $a_{-}{\psi }_0=0$ iswhose solution is proportional to $e^{-m\omega x^2/(2\hslash )}$. This Gaussian function is square-integrable, so it gives an actual ground state. Its normalization can be chosen using the Gaussian integral.

All the displayed levels occur: for a nonzero state at $E_n$,Thus repeated raising creates a nonzero energy eigenstate at every $E_n$.

The squared norm is obtained by moving the adjoint operator to the other side of the inner product:For an energy eigenstate, this isSince a squared norm is nonnegative, every nonzero energy eigenstate has $E\ge \hslash \omega /2$. Moreover $a_{-}\psi =0$ exactly when $E=\hslash \omega /2$.

Let $\psi$ have energy $E_n$ and put $\gamma =2m\hslash \omega$. The lowering norm formula gives, inductively,In particular $(a_{-}{)}^n\psi$ is a nonzero ground state. By the assumed nondegeneracy, it equals $c{\psi }_0$ for some $c\ne 0$.

It remains to show that raising this vector recovers $\psi$, rather than merely another state at the same energy. On an energy eigenstate at $E_n$, factorization gives $a_{+}{a}_{-}\psi =\gamma n\psi$. Iterating the factorization down the ladder yieldsApplying $(a_{+}{)}^n$ to $(a_{-}{)}^n\psi =c{\psi }_0$ therefore provesFor $n=0$ this is the original ground-state assumption. This proves nondegeneracy propagated by an oscillator ladder: every energy eigenspace is one-dimensional if the ground state is nondegenerate.

If $a_{-}\psi \ne 0$, it has energy $E-\hslash \omega$, distinct from the energy of $\psi$. Eigenstates of a self-adjoint operator with distinct eigenvalues are orthogonal. Explicitly,so $⟨\psi ,a_{-}\psi ⟩=0$. If $a_{-}\psi =0$, the same conclusion is immediate.

Dividing by $\Vert \psi {\Vert }^2$ and using the definition of $a_{-}$ givesBoth expectation values are real because $P$ and $X$ are self-adjoint operators. Its real and imaginary parts therefore vanish separately:An oscillator energy eigenstate is centered in both position and momentum, despite having nonzero quantum uncertainties.

Inside the infinite square well, the Schrödinger equation and Dirichlet boundary conditions areThe wavefunction vanishes outside the box. These walls constrain its values, not its first derivatives: requiring both to vanish at an endpoint would eliminate all the nonzero stationary solutions.

For a nonzero separated solution, separation of variables givesThe last equation is the time-independent Schrödinger equation. Its energy eigenvalue is real. Multiplication by ${\psi }^{\ast }$ and integration by parts, with the endpoint values zero, givesEquality would force $\psi$ to be constant and hence zero, so a nonzero energy eigenstate has $E>0$. Writing $k=\sqrt{2mE}/\hslash$, the general solution is $A\sin kx+B\cos kx$. The left wall gives $B=0$; the right gives $\sin ka=0$. Negative integers reproduce the same functions up to a constant, and $n=0$ gives the zero function, leaving $k=n\pi /a$ for $n=1,2,\dots$.

Finally ${\int }_0^a{\sin }^2(n\pi x/a)dx=a/2$. Thus wavefunction normalization with positive coefficient givesAn arbitrary constant global phase has no effect on the quantum state.

Write $k=\pi /a$, $A=\sqrt{2/a}$ and ${\Psi }_0=A{e}^{ikx}\sin kx$. Its derivative isIntegration by parts with the Dirichlet boundary conditions gives the energy expectationIt is constant in time, since the time-independent Hamiltonian operator has fixed energy measurement probabilities: $⟨H{⟩}_t={\sum }_n{E}_n|c_n{|}^2$.

For the midpoint parity operator, $R{\Psi }_0=-A{e}^{-ikx}\sin kx$. Hence, again putting $u=kx$,Indeed the imaginary part integrates to zero by reflection and ${\int }_0^{\pi }\cos 2u{\sin }^{2}udu=-\pi /4$. This expectation value is also constant in time, because the bounded self-adjoint operator $R$ commutes with $H$ and hence with its time-evolution operator. Equivalently, $⟨R{⟩}_t={\sum }_n(-1{)}^{n+1}|c_n{|}^2$.

For the momentum operator, direct evaluation givesThe real sine-cosine integral vanishes. This expectation value is not constant in time: the infinite square well walls reverse motion. One can prove the change exactly without approximating the wavefunction. Atthe energy eigenstate expansion and $R{\psi }_n=(-1{)}^{n+1}{\psi }_n$ give $\Psi (t_{\ast })=-R{\Psi }_0$. Differentiation shows $PR=-RP$, soThis half revival in an infinite square well is a concrete counterexample to momentum conservation. The sign from $-R$ is only a global phase.

The apparent interior identity $[H,P]=0$ cannot justify conservation here: the momentum operator does not preserve the Dirichlet boundary conditions of the Hamiltonian operator's operator domain. For sufficiently regular states, retaining the boundary terms in Ehrenfest theorem instead gives the boundary force in an infinite square well,In the strictly confined interval, $-i\hslash {\partial }_x$ with both endpoint values zero is symmetric but not itself a self-adjoint operator. The requested integral is the conventional mean momentum; it also equals the full-line momentum mean of the zero-extended state, as measured upon releasing the walls. This domain distinction explains why the formal interior commutator misses the wall force.

The requested energy eigenvalue is $E_2$. By the Born rule, its probability amplitude is the inner product $c_2=⟨{\psi }_2,\Psi (0)⟩$. Put $u=\pi x/a$ and use $e^{iu}\sin u=\frac{1}{2}\sin 2u+i{\sin }^{2}u$:The first integral is $\pi /2$; the second is zero under $u\mapsto \pi -u$. ThereforeFor the time-independent Hamiltonian operator, expand in its orthonormal basis of energy eigenstates:Unitary time evolution changes $c_2$ to $c_2{e}^{-i{E}_{2}t/\hslash }$, whose modulus is unchanged. Thus the same energy measurement probability holds at every time, even though relative quantum phases and other measurements can change.

The spatial quantum phase has unit modulus. Thus the probability density and total probability areThe exterior contributes zero. This establishes wavefunction normalization, although the spatial quantum phase means the state is not the ground energy eigenstate.

On a smooth wavefunction $f$, the position operator and momentum operator giveThus the canonical commutation relation is $\boxed{[X,P]=i\hslash I}$. Since $X$ and $P$ are self-adjoint operators and $m,\omega$ are real, the adjoint operator of $a_{-}$ isExpanding the products in the stated order is essential:Subtracting gives $\boxed{[a_{-},a_{+}]=2m\hslash \omega I}$. These identities initially hold on smooth rapidly decaying functions, which contain the quantum harmonic oscillator energy eigenstates used below.

The factorization also proves a lower bound: for a normalized energy eigenstate $H\psi =E\psi$,Therefore $E\ge E_0:=\hslash \omega /2$. The commutator product rule, applied to $H=a_{+}{a}_{-}/(2m)+E_0$, givesSo a nonzero application of the creation operator $a_{+}$ raises the energy eigenvalue by $\hslash \omega$, while the annihilation operator $a_{-}$ lowers it by the same amount.

Begin with any energy eigenstate, whose existence is assumed. Repeated lowering must terminate: otherwise it would eventually produce a nonzero energy eigenstate below $E_0$, contradicting the lower bound. More explicitly, if $(E-E_0)/(\hslash \omega )=N+r$ with integer $N\ge 0$ and $0\le r<1$, the first $N$ lowering steps are nonzero. If $r>0$, one further step is also nonzero by its positive squared norm and has energy $E_0+(r-1)\hslash \omega <E_0$, an impossibility. Thus $r=0$ and the last state satisfies $a_{-}{\psi }_0=0$. This provesEvery level exists as well: $\Vert a_{+}{\psi }_n{\Vert }^2=2m(E_n+\hslash \omega /2)=2m\hslash \omega (n+1)>0$, so repeated raising from ${\psi }_0$ never vanishes. The same factorization bounds the energy expectation of any normalized state in the one-dimensional harmonic oscillator form domain below by $E_0$, proving that the terminated state really is a ground state.

To find its probability density, solve the first-order equationIts solutions are ${\psi }_0=C_0\exp [-m\omega x^2/(2\hslash )]$. With $\alpha =m\omega /\hslash >0$, the supplied Gaussian integral givesTaking $C_0=(\alpha /\pi {)}^{1/4}$ fixes an irrelevant global phase. ConsequentlyThe ground-state energy is positive even though the potential minimum is zero: it is the oscillator's zero-point energy.