Statistical Mechanics of a Quantum Harmonic Oscillator with Folded Gaussian Frequency

A self-contained statistical-mechanics treatment of a single quantum harmonic oscillator is presented, whose frequency $\omega$ is drawn from a folded Gaussian distribution: $\omega =\mathrm{\mid }\xi \mathrm{\mid }$ with $\xi \sim \mathcal{N}(\mu ,{\sigma }^2)$. The exact integral representations for the partition function, internal energy, free energy, heat capacity, and entropy are derived, and analytic approximations are given in two complementary limits—small variance ($\sigma \ll \mu$) via a cumulant expansion, and the zero-center case ($\mu =0$) via low-frequency asymptotic analysis. The model is extended to $N$ independent oscillators, where the heat capacity is shown to be extensive with self-averaging fluctuations $\propto N^{-1\mathrm{/}2}$, and finally to a disordered oscillator lattice, where the folded-Gaussian kink at $\omega =0$ produces a soft-mode infrared tail. For a single isolated oscillator with $\mu =0$, both $C$ and $S$ vanish linearly at low $T$. In the lattice case, the van Hove factor converts this to a $T^d$ power law. The oft-quoted ``third-law violation’’ for disordered phonons is here shown to be a spectral property—the absence of an energy gap and a power-law freeze-out—driven by the single-site distribution kink rather than by a genuine Lifshitz tail (which requires rare large-scale spatial fluctuations). The folded Gaussian thus serves as a minimal benchmark for soft-mode disorder thermodynamics.