Set $u=A\psi$, $v=B\psi$ and $z=⟨u,v⟩$. The self-adjointness of $A,B$ impliesThe Cauchy-Schwarz inequality followed by $|z{|}^2\ge (\mathrm{Im}z{)}^2$ givesBoth inequalities must be equalities. The exact conditions are $A\psi$ and $B\psi$ are linearly dependent, and $\mathrm{Re}⟨A\psi ,B\psi ⟩=0$. Equivalently, equality holds if either vector is zero; when both are nonzero it holds precisely whenThe zero-vector cases matter: they saturate the inequality but cannot all be expressed by solving for $B\psi$ as a multiple of $A\psi$. The real-part condition can also be written ${\mathbb{E}}_{\psi }(AB+BA)=0$, using the anticommutator. This is the uncentered second-moment commutator inequality.

For a normalized wavefunction, the expectation value of a function of position isprovided the integral exists. For an unnormalized state, divide this expression by $\int |\psi {|}^{2}dx$. Put $b=m\omega /\hslash >0$. The ground state has probability density $\sqrt{b/\pi }{e}^{-b{x}^2}$, which is even. ThusThe zero mean position ${\mathbb{E}}_{{\psi }_0}(x)=0$ expresses symmetry, whereas this positive expectation value measures the mean distance from the centre.

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.

Substitute the power series and collect each power of $\xi$:Matching the power series coefficients gives the recurrence relationThe even and odd coefficients form separate sequences, determined by $a_0$ and $a_1$, respectively. For example, $a_2=(1-ϵ)a_0/2$ and $a_3=(3-ϵ)a_1/6$. These are the two independent power series solutions of the Hermite differential equation.

The dimensionless coordinate for the quantum harmonic oscillator gives $d/dx=\sqrt{m\omega /\hslash }d/d\xi$. The time-independent Schrödinger equation therefore becomesFor $\psi =e^{-{\xi }^2/2}f$, direct differentiation yieldsCancelling the nonzero exponential reduces the equation to the Hermite differential equation

At the midpoint the second sine vanishes, whereas the derivative of the first sine vanishes. Thus the wavefunction and its spatial derivative simplify toSubstitution in the probability current givesThe continuity equation, integrated over the left half, gives $P_{\mathrm{L}}^{\prime }(t)=j(0,t)-j(a/2,t)$. The Dirichlet boundary condition makes $j(0,t)=0$, and the preceding probability formula givesPositive midpoint current removes probability from the left half. This is precisely the midpoint current of a two-mode box state.

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.

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 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.

A density operator is a positive self-adjoint trace-class operator with trace one. In finite dimensions this is a positive semidefinite Hermitian matrix with $\mathrm{tr}\rho =1$. Its expectation of an observable is $\boxed{⟨A⟩=\mathrm{tr}(A\rho )}$ whenever the expectation is defined.

In the source's square-root notation, take the unique positive square root $W=\sqrt{\rho }$. It has Hilbert-Schmidt norm one and satisfies $W{W}^{\ast }=\rho$, soThe standard vector purification is the associated vector $|W⟩⟩={\sum }_{r,s}{W}_{rs}|{\xi }_r⟩|{\eta }_s⟩$ on a system and a chosen copy. Tracing out the copy gives $\rho$; measuring $A\otimes I$ gives the same trace formula. Thus the operator square root and the purified vector are two related objects, distinguished by this vectorization.

For a general coefficient matrix $C$, expansion of the expectation givesThus $\boxed{{\rho }_1=C{C}^{\ast }}$. Positivity is immediate from $v^{\ast }C{C}^{\ast }v=\Vert C^{\ast }v{\Vert }^2$, and its trace is one by normalization of $\Psi$.

If $C$ is nonsingular, it is square and ${\rho }_1$ is positive definite. DefineThen $U{U}^{\ast }=C^{-1}{\rho }_1(C^{\ast }{)}^{-1}=I$, which for a square matrix also implies $U^{\ast }U=I$. Thus $U$ is unitary, $\sqrt{{\rho }_1}=CU$, and $C=\sqrt{{\rho }_1}{U}^{\ast }$.

In an eigenbasis of ${\rho }_1$, its square root is $\mathrm{diag}(\sqrt{p_r})$. Hence $c_{rs}=\sqrt{p_r}\overline{U_{sr}}$. DefineUnitarity gives $⟨{\zeta }_r|{\zeta }_t⟩={\delta }_{rt}$, and substitution gives the Schmidt decompositionFor singular or rectangular $C$, the same conclusion follows on its support by a singular value decomposition; the stated nonsingular case needs only the unitary calculation above.

For the supplied Bell state, the requested matrices areThe Schmidt vectors can be chosen as ${\zeta }_1={\eta }_2$, ${\zeta }_2=-{\eta }_1$, and the two coefficients are $1/\sqrt{2}$. The square root's canonical vector purification is $(|{\xi }_1{\eta }_1⟩+|{\xi }_2{\eta }_2⟩)/\sqrt{2}$; the original Bell vector is a different purification of the same reduced state, related by a unitary on the second factor. This illustrates the unitary freedom of purification.

Let $E=f_H(\psi )$. For any admissible variation $\eta$, differentiation with respect to a real parameter givesStationarity for all variations, including their imaginary multiples, forces $(H-E)\psi =0$; conversely this eigenvector equation makes the derivative vanish. Thus stationary Rayleigh quotients are exactly eigenvectors, with their eigenvalues as the quotient values.

Put $\kappa =e^2/(4\pi {ϵ}_0)$ and assume $\alpha >0$. The squared norm of the radial trial function is $4\pi {\int }_0^{\infty }{r}^2{e}^{-2\alpha r}dr=\pi /{\alpha }^3$. Its kinetic expectation, computed through $\int |\nabla \psi {|}^2$, is ${\hslash }^2{\alpha }^2/(2m)=\kappa a{\alpha }^2/2$. Splitting the potential integral at $R$ givesFor example, the first integral is $[1-e^{-2s}(1+2s+2s^2)]/(4{\alpha }^3)$ and the second is $e^{-2s}(R/(2\alpha )+1/(4{\alpha }^2))$, where $s=\alpha R$. ThereforeIts derivative is $\kappa [a\alpha -(1+2\alpha R)e^{-2\alpha R}]$. The actual stationary condition is consequentlyThe extra Coulomb factor in the printed stationary equation is an error: the common prefactor cancels on differentiation. The corrected equation is dimensionless on both sides. Its left side increases from zero, while its right side decreases from one, so it has one positive solution and gives the minimum in this trial family.

Let $\epsilon =R/a$ and $y=a\alpha$. Expansion of $(1+2s)e^{-2s}$ gives $1-2s^2+O(s^3)$. ThusThere is no linear correction in $R/a$. By the quantum variational principle, any normalized trial state gives an upper bound on the true ground energy. Choosing $\alpha =1/a$ already gives the exact boundIts small-radius expansion isUsing the optimized $\alpha$ improves this bound starting at order ${\epsilon }^4$, and hence gives the same leading finite-size correction. The true energy is also at least the point-nucleus value $-\kappa /(2a)$, since the finite-radius potential is nowhere more attractive.

A permutation $\pi$ of identical particles multiplies a fermionic state by $\mathrm{sgn}(\pi )$; a bosonic state is invariant. Particle exchange keeps $y_{+}$ fixed and sends $y_{-}$ to $-y_{-}$. Consequently the spatial wave function is even in $y_{-}$ for the stated bosonic model and odd for the stated fermionic model. This assumes no additional spin variables in the supplied model; exchange of full particle states is the underlying rule.

The coordinate change is orthogonal, with unit Jacobian determinant, so the kinetic energy becomeson the square $|y_{+}|,|y_{-}|<a$ with zero boundary values. Product states ${\psi }_r(y_{+}){\psi }_s(y_{-})$ have energyThe one-dimensional eigenfunction is even for odd $s$ and odd for even $s$. Bosons therefore allow odd $s$, whereas fermions require even $s$. The lowest bosonic pair is $(r,s)=(1,1)$; its first excited pair is $(2,1)$. The lowest fermionic pair is $(1,2)$. ThusThe printed bosonic ground-energy claim omits the factor ${\pi }^2/4$; the value above follows directly from the supplied one-dimensional spectrum and kinetic energy.

The separation is $|x_1-x_2|=\sqrt{2}|y_{-}|$. The normalized $y_{+}$ factor integrates out. In the fermionic ground state, the relative-coordinate density is $a^{-1}{\sin }^2(\pi y_{-}/a)$, givingFor the bosonic ground state it is $a^{-1}{\cos }^2(\pi y_{-}/(2a))$, givingBoth integrals follow by replacing squared trigonometric functions by their constant and cosine parts and integrating once by parts. Henceand the difference is $\boxed{2\sqrt{2}a/{\pi }^2}$. The odd relative state vanishes at coincidence, which increases its expected separation.

For a normalized state and a self-adjoint operator, the expectation value and dispersion areUse the physics convention that the inner product is conjugate-linear in its first argument, and assume the displayed products and second moments are defined. Self-adjointness gives $\overline{⟨\psi |A\psi ⟩}=⟨A\psi |\psi ⟩=⟨\psi |A\psi ⟩$, so the expectation is real. MoreoverIt vanishes exactly when $\psi$ is an eigenstate, with eigenvalue its expectation.

The product rule gives $XPf=-i\hslash x{f}^{\prime }$ and $PXf=-i\hslash (f+x{f}^{\prime })$, hence $\boxed{[X,P]f=i\hslash f}$ on their common differentiable domain. For general operators, put $z=⟨A\psi |B\psi ⟩$. ThenThe Cauchy-Schwarz inequality yieldsApply this to centered position and momentum operators, whose commutator is unchanged. Their squared norms are the variances, givingEquality requires both proportionality in Cauchy-Schwarz and a purely imaginary inner product. Thus for real means $x_0,p_0$,with real $\gamma$. Solving this first-order equation givesNormalizability on the real line forces $\gamma >0$. Writing ${\sigma }^2=\hslash /(2\gamma )$, all normalized equality states are thereforeThey have $\Delta X=\sigma$ and $\Delta P=\hslash /(2\sigma )$. A nonconstant quadratic phase would increase the product, so the full family here has a real positive Gaussian width and only a linear spatial phase.

For bound hydrogen atom states, the radial wave function is regular and finite at the origin, and decays at infinity so that $4\pi {\int }_0^{\infty }{r}^2|\Psi (r){|}^{2}dr<\infty$. Equivalently the reduced radial function $w(r)=r\Psi (r)$ has $w(0)=0$ and is square integrable on the half-line. A singular $1/r$ branch is not an admissible regular Coulomb eigenfunction; its derivative energy diverges and it introduces an origin singularity absent from the Hamiltonian domain.

Write $\kappa =e^2/(4\pi {ϵ}_0)$, so $a={\hslash }^2/(m\kappa )$. Multiplication by $r$ turns the radial equation intoSince $w=a{e}^{-\lambda R}f(R)$ and $r=aR$, substitution givesChoose $\boxed{\lambda =\mu a>0}$ to remove the constant coefficient and factor out the decaying exponential. For $f={\sum }_{k\ge 0}{a}_k{R}^{k+c}$, the lowest power has coefficient $c(c-1)a_0$, so $\boxed{c(c-1)=0}$. The $c=0$ behavior would give $\Psi \sim a_0/R$, violating regularity; choose $c=1$.

Matching successive powers then yieldsA nonterminating tail eventually has one sign, and its coefficient ratio is asymptotic to $2\lambda /k$. More precisely, for any $0<\delta <\lambda$, the ratio is eventually at least $(2\lambda -\delta )/(k+1)$. Its same-sign tail consequently dominates a positive constant times the exponential series of rate $2\lambda -\delta$, up to finitely many terms. The factor $e^{-\lambda R}$ cannot then make $\Psi$ decay. Thus square integrability forces termination.

If the last nonzero term is $a_{n-1}{R}^n$, termination requires $\lambda n=1$. ThereforeFor $n=1$, $\lambda =1$, $f$ is proportional to $R$, andThe radial volume element is $4\pi r^{2}dr$, so the mean radius is the ratio of the radius-weighted norm integral to the norm integral. With $r=aR$, the common angular and normalization factors cancel:In the ground state $\boxed{⟨r⟩=3a/2}$.

In the position representation the Time-dependent Schrödinger equation isFor a normalized state define probability density and probability current byMultiplying the equation by $\overline{\psi }$ and its complex conjugate by $\psi$, then subtracting, cancels the real potential. The resulting identity isThus $\boxed{{\partial }_t\rho +{\partial }_{x}j=0}$. Integrating over an interval gives $d{\int }_b^c\rho dx/dt=j(b)-j(c)$: probability changes only through flux across the endpoints. This is local conservation of quantum probability.

Continuity with the zero exterior imposes zero endpoint values on the infinite square well. Solving the stationary equation with these boundary conditions giveswith zero extension outside. The integral of ${\sin }^2(n\pi x/a)$ is $a/2$, which fixes the normalization. The initial coefficients are $3/5$ and $4/5$, whose squared moduli sum to one. Each energy component acquires its stationary phase:At the midpoint, ${\Psi }_1=\sqrt{2/a}$, ${\Psi }_1^{\prime }=0$, ${\Psi }_2=0$, ${\Psi }_2^{\prime }=-(2\pi /a)\sqrt{2/a}$. HenceTaking its imaginary part proves

For classical events, $A\setminus C\subset (A\setminus B)\cup (B\setminus C)$: an outcome in $A$ but not in $C$ either fails $B$ or belongs to $B\setminus C$. Monotonicity and the union bound prove the Fréchet-Bell inequalityMalus' law gives transmission ${\cos }^2\gamma$ for a photon polarized by the first filter, so its conditional rejection probability at the second is $\boxed{{\sin }^2\gamma }$. Conditional and joint probabilities must be distinguished: $P(A\setminus B)=P(A)P(B^c\mid A)$.

For equal unbiased marginals $P(A)=P(B)=1/2$, the quantum two-setting rejection probabilities are $\frac{1}{2}{\sin }^2({\theta }_A-{\theta }_B)$. A classical common assignment of all three outcomes would therefore require ${\sin }^2({\theta }_A-{\theta }_B)+{\sin }^2({\theta }_B-{\theta }_C)\ge {\sin }^2({\theta }_A-{\theta }_C)$. Choose directions $0,\gamma ,2\gamma$. ThenFor example $\gamma =\pi /6$ gives $1/2<3/4$. More generally ${\sin }^{2}a+{\sin }^{2}b-{\sin }^2(a+b)=-2\sin a\sin b\cos (a+b)$. In a triangle whose side directions are these filters, the internal angle supplementary to $a+b$ is obtuse exactly in this violation range. This specifies the angle convention behind the triangle description: with polarization axes regarded as unoriented lines, a claim that every sufficiently large physical axis separation violates the inequality would not be correct. The contradiction concerns a single classical joint model for all settings, not the validity of ordinary probability or a comparison of arbitrarily renormalized sequential experiments.

Use $|0⟩=|{\alpha }_0⟩$ and $|1⟩=|{\alpha }_1⟩$. Projecting the first photon with $|0⟩⟨0|$ gives the unnormalized product vector $c_0|00⟩$, with success probability $|c_0{|}^2$. Projection of the second onto $|\theta ⟩=\cos \theta |0⟩+\sin \theta |1⟩$ givesConditioned on success, normalize these vectors; a zero amplitude means that outcome never occurs. For the maximally entangled equal-amplitude pair, rotation of both real polarization bases gives the same correlations for arbitrary settings. Comparing the appropriately normalized coincidence/rejection frequencies at three setting pairs tests the inequality above. A general unequal-amplitude state need not reproduce the unbiased rotationally invariant probabilities, so the maximally entangled choice supplies the clean explicit example.

For swapping, label the four systems in order ${\xi }_1,{\eta }_1,{\xi }_2,{\eta }_2$ and let $|\Phi {⟩}_{\eta }=2^{-1/2}{\sum }_{j=0}^1|{\eta }_j{\eta }_j⟩$. Its overlap with $|{\eta }_r{\eta }_s⟩$ is ${\delta }_{rs}/\sqrt{2}$. Applying its rank-one projector to Alice's two systems givesAfter regrouping the remote and local systems, this factorizes into $[{\sum }_r{b}_r^2|{\xi }_r{\xi }_r⟩/\sqrt{2}]\otimes |\Phi {⟩}_{\eta }$. The selected outcome has probability $p=\frac{1}{2}{\sum }_r|b_r{|}^4$. Therefore the normalized conditional remote state isIt has Schmidt rank at least two whenever the original pair was entangled, since squaring a nonzero coefficient does not make it zero. Thus the two photons that never interacted become entangled conditional on Alice's joint measurement outcome. This is entanglement swapping. Alice is projecting onto a specified Bell state, not performing a universal test of whether an arbitrary unknown state is entangled. Averaging all her outcomes leaves the remote reduced density operator unchanged, so no faster-than-light signal is available without her classical outcome label.

For a nonzero trial vector use the Rayleigh quotient $\mathcal{E}(\psi )=⟨\psi ,H\psi ⟩/⟨\psi ,\psi ⟩$. Its first variation in any allowed direction $\eta$ isVanishing for all complex directions implies $(H-\mathcal{E})\psi =0$; an eigenvector conversely makes it vanish. The ground energy is the infimum of this quotient for a lower-bounded Hamiltonian, so minimizing over a trial family provides an upper bound to it, exact when that family contains a ground eigenvector.

Set $a=⟨\varphi ,H\varphi ⟩$, $b=⟨\varphi ,I\varphi ⟩$ and $c=⟨\varphi ,HI\varphi ⟩$. All three are real: $H,I$ are self-adjoint and commute. Also $I^{†}I=I^2=1$. Expansion of the numerator and norm givesDifferentiate the quotient, retaining only nonzero trial vectors:Thus the stationarity equation is $\boxed{(t^2-1)(c-ab)=0}$. If $c=ab$, the quotient is identically $a$ wherever defined. Otherwise $t=\pm 1$ are the two extrema, withIn this nonconstant case $|b|<1$, so both denominators are positive. If $c-ab>0$, the plus value is the maximum and the minus value the minimum; the assignments reverse when $c-ab<0$. One can also see global extremality by resolving the trial vector into the orthogonal $I=\pm 1$ sectors: its energy is a weighted mean of their two fixed projected energies.

The null-vector exception matters. If $|b|=1$, then $I\varphi =\pm \varphi$, the covariance $c-ab$ vanishes and one of $\varphi \pm I\varphi$ is zero. That value has no Rayleigh quotient and is not an extra stationary physical state.

For particle exchange, change of integration variables shows that $I\varphi (x_1,x_2)=\varphi (x_2,x_1)$ is unitary and $I^2=1$, hence self-adjoint. Exchanging coordinates interchanges the two kinetic derivatives and leaves the symmetric potential unchanged, so $[I,H]=0$. The vectors $\varphi \pm I\varphi$ lie in its symmetric and antisymmetric sectors. Their normalized nonzero versions are respectively bosonic and fermionic spatial trial wavefunctions.

Because the energy of a general trial state is the weighted average of its sector energies, at least one nonzero symmetry projection has energy no greater than the original trial. Thus minimizing can be performed in a definite exchange-symmetry sector. This is a variational approximation within the trial subspace; it does not assert that projecting an arbitrary trial vector produces an exact eigenvector of $H$. For particles with a prescribed statistics sector one minimizes within that physical sector.