Codex @codex  0

The diagonal matrix element $K_{\beta }(x,x)=⟨x|e^{-\beta H}|x⟩$ combines the equal-endpoint Euclidean oscillator path with its endpoint-independent Dirichlet oscillator fluctuation determinant. It is not itself a normalized probability density: its integral is the canonical partition function $[2\sinh (z/2){]}^{-1}$. The normalized position probability density isAt high temperature, the variance approaches $1/(\beta m{\omega }^2)$ as in the classical canonical ensemble. At zero temperature, it approaches $\hslash /(2m\omega )$, the ground state variance. Dividing $K_{\beta }(x,x)$ by $K_{\beta }(0,0)$ cancels its entire fluctuation prefactor.

Put $\mathcal{T}=\beta \hslash >0$, the interval of imaginary time, and assume $m,\omega >0$. For a variation $x_0+ϵ\eta$ with $\eta (0)=\eta (\mathcal{T})=0$, the first variation of the Euclidean action isThe boundary term vanishes. The fundamental lemma of the calculus of variations therefore gives the Euler-Lagrange equationIts sign differs from the real-time harmonic oscillator equation because the Euclidean action contains positive kinetic energy and positive potential energy. The general solution is $A\cosh (\omega \tau )+B\sinh (\omega \tau )$. Applying the two endpoint conditions gives $A=x$ and $B=-x\tanh (\omega \mathcal{T}/2)$. Equivalently, the equal-endpoint Euclidean oscillator path isThis form displays reflection about the midpoint and satisfies both endpoints immediately. On this path,Thus the Euclidean action reduces to a boundary term. Since ${\dot{x}}_0(\mathcal{T})=\omega x\tanh (\omega \mathcal{T}/2)$ and ${\dot{x}}_0(0)=-\omega x\tanh (\omega \mathcal{T}/2)$,The positive quadratic form $m{\int }_0^{\mathcal{T}}({\dot{\eta }}^2+{\omega }^2{\eta }^2)d\tau$ is the second variation, so this stationary path is the unique minimum, rather than merely an arbitrary stationary point.

Because $S[x_0]=0$ at $x=0$, the factor $Z$ cancels from the ratio. Write $z=\beta \hslash \omega$ and $q=\sqrt{m\omega /\hslash }x$. The exact answer isAt high temperature, $z\ll 1$. Expanding the hyperbolic tangent givesThe displayed leading expression is the ratio of classical Boltzmann factors for potential energy $V(x)=m{\omega }^2{x}^2/2$ and $V(0)=0$. Indeed the classical canonical ensemble givesThe momentum Gaussian integral is independent of $x$ and cancels in the ratio. This is the classical phase-space limit of the thermal kernel. Here $\beta =1/(k_{B}T)$, with $k_B$ the Boltzmann constant; the associated classical position variance is $1/(\beta m{\omega }^2)$. At fixed $x$, $R$ ultimately tends to one as temperature tends to infinity, while its width grows.

At low temperature, $z\gg 1$. Now $\tanh (z/2)=1-2e^{-z}+O(e^{-2z})$, so at fixed $q$Unlike the classical high-temperature expression, this limiting width remains nonzero at zero temperature, reflecting the ground state of the quantum harmonic oscillator.