🔬 condensed matter

Statistical Mechanics of a Quantum Harmonic Oscillator with Folded Gaussian Frequency

This paper presents a statistical-mechanical analysis of quantum harmonic oscillators with frequencies drawn from a folded Gaussian distribution, deriving exact thermodynamic quantities and demonstrating that the distribution's kink at zero frequency induces a soft-mode infrared tail that leads to a $T^d$ power-law heat capacity in lattices, thereby clarifying that the resulting low-temperature behavior stems from single-site spectral properties rather than rare Lifshitz tails.

Original authors: Liu Zhao

Published 2026-08-18

📖 7 min read🧠 Deep dive

CC BY 4.0

Original authors: Liu Zhao

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the quiet world of quantum physics, the simplest building block is the harmonic oscillator. Imagine a tiny weight attached to a spring, bobbing back and forth. In the classical world, this motion can happen at any speed, but in the quantum realm, the energy comes in specific, fixed packets, and the speed of the bobbing is determined by a single number: the frequency. For decades, physicists have used this simple model to understand how matter behaves when it is cold, how heat moves through solids, and why the laws of thermodynamics hold true even at the coldest temperatures imaginable. Usually, scientists assume this frequency is a fixed, unchanging value. However, real materials are rarely perfect. They contain impurities, defects, and random variations that make the properties of their atoms fluctuate from place to place. When these random variations are introduced into the frequency of the oscillator, the system becomes a window into the messy, disordered reality of condensed matter. The question then becomes: how does this randomness change the way the material stores energy and responds to temperature?

A researcher at Nankai University has recently tackled this question by studying a specific kind of randomness. Instead of assuming the frequency is a fixed number or a completely chaotic variable, the study treats the frequency as a value drawn from a "folded Gaussian" distribution. To visualize this, imagine taking a standard bell curve of possible values and folding it in half so that all negative numbers are flipped to positive. This creates a distribution where the frequency is always positive, but it has a sharp, non-smooth point, or "kink," right at zero. This kink is crucial because it allows for a significant number of oscillators to have very low frequencies, creating a "soft mode" that behaves differently than the smooth, predictable tails found in other models of disorder. The researcher calculated exactly how a single such oscillator, and then a whole lattice of them, would behave thermodynamically, deriving precise mathematical descriptions for their energy, heat capacity, and entropy without relying on approximations that might fail in extreme conditions.

The findings reveal a clear distinction between a single isolated oscillator and a large collection of them. For a single oscillator with this specific type of randomness, the study shows that as the temperature drops toward absolute zero, the system does not freeze instantly. Instead, its ability to absorb heat, known as heat capacity, decreases in a straight line as the temperature falls. This linear drop means the system still has some freedom to move even at very low temperatures, but it eventually settles down completely. Crucially, the entropy, which measures the disorder of the system, also drops to zero as the temperature approaches absolute zero. This confirms that the third law of thermodynamics holds true even with this disorder; the system does not get stuck in a state of permanent confusion. The randomness simply slows down the freezing process, making it gradual rather than abrupt, but it does not break the fundamental laws of physics.

When the researcher expanded the model to include a large number of independent oscillators, the behavior remained consistent with the laws of thermodynamics, but the scale changed. The total heat capacity of the group became proportional to the number of oscillators, a property known as extensivity. Furthermore, the study demonstrated that the variations between different samples of these oscillators become negligible as the group gets larger. This "self-averaging" effect means that if you were to measure a large block of this disordered material, the result would be predictable and stable, regardless of the specific random arrangement of frequencies in that particular sample. The fluctuations in the measurement shrink as the square root of the number of particles, ensuring that the macroscopic properties are reliable and well-defined.

The most striking results appear when the disordered material is arranged in a lattice, like the atoms in a crystal. In a perfect crystal, the way heat capacity drops at low temperatures is governed by the geometry of the lattice and the speed of sound waves traveling through it. In this disordered lattice, the sharp kink at zero frequency in the distribution combines with the geometry of the lattice to produce a new power-law behavior. Instead of the heat capacity dropping exponentially fast as it does in perfect crystals with a gap in energy, it drops much more slowly, following a power law that depends on the dimension of the space. For a three-dimensional lattice, the heat capacity drops with the cube of the temperature. This creates a "soft-mode tail" where low-energy vibrations persist much longer than expected. The study clarifies that this behavior is not a violation of the third law, nor is it the result of rare, large-scale fluctuations often called Lifshitz tails in other contexts. Instead, it is a direct consequence of the simple, local kink in the frequency distribution. The system freezes, but it does so in a way that is dictated by the sheer number of available low-frequency states, rather than by a lack of energy gaps.

This work serves as a minimal benchmark for understanding how disorder affects quantum systems. By isolating the effect of a simple, folded frequency distribution, the researcher has shown that the "soft" behavior seen in some disordered materials does not require complex, large-scale spatial fluctuations to occur. It can arise from the simple statistical properties of the frequencies themselves. The study provides exact formulas for these behaviors, allowing physicists to predict how such systems will respond to temperature changes without needing to run complex simulations. It confirms that while disorder can change the rate at which a system freezes, it does not necessarily prevent it from reaching a state of zero entropy. The results offer a clear, mathematically rigorous picture of how randomness interacts with the fundamental laws of heat and motion, distinguishing between the behavior of a single particle and the collective behavior of a vast array of them.

The implications of these findings extend to how we understand the thermodynamics of real-world materials that are never perfectly ordered. The study suggests that the "third-law violation" sometimes discussed in literature regarding disordered phonons is actually a matter of terminology and spectral properties rather than a true failure of thermodynamics. The system does not retain a residual entropy; it simply freezes more slowly and follows a different mathematical path to get there. This distinction is vital for interpreting experimental data on disordered solids, where the presence of low-frequency modes can lead to unexpected heat capacity measurements. The research establishes that these measurements are consistent with standard thermodynamic principles, provided one accounts for the specific way the disorder is distributed. By providing a clean, solvable model, the work offers a reference point for more complex theories and helps separate the effects of simple frequency disorder from more exotic phenomena like Anderson localization or true Lifshitz tails.

Ultimately, the paper demonstrates that the statistical mechanics of a quantum harmonic oscillator with folded Gaussian disorder is a tractable and well-behaved system. It bridges the gap between the idealized world of perfect crystals and the messy reality of disordered matter. The key takeaway is that the specific shape of the frequency distribution, particularly the presence of a kink at zero, dictates the low-temperature behavior. This leads to a linear drop in heat capacity for a single oscillator and a power-law drop for a lattice, both of which satisfy the third law of thermodynamics. The study confirms that disorder can be a benign, predictable force that modifies the rate of freezing without destroying the fundamental order of the universe. Through exact calculations and careful analysis of limits, the researcher has provided a clear map of how these quantum systems behave, offering a solid foundation for future explorations into the thermodynamics of disordered matter.

Technical Summary: Statistical Mechanics of a Quantum Harmonic Oscillator with Folded Gaussian Frequency

Problem Statement
This work addresses the statistical mechanics of a quantum harmonic oscillator subject to frequency disorder, specifically where the frequency $\omega$ is drawn from a "folded Gaussian" distribution ( $\omega = |\xi|$ , with $\xi \sim \mathcal{N}(\mu, \sigma^2)$ ). While static disorder in oscillator systems is a known phenomenon linked to Anderson localization and Lifshitz tails, this model isolates a specific, minimal source of disorder: the non-analytic "kink" at $\omega=0$ introduced by the absolute value operation on a Gaussian variable. The paper seeks to determine how this specific spectral feature—a soft-mode tail rather than a genuine Lifshitz tail driven by rare spatial fluctuations—affects thermodynamic quantities (partition function, internal energy, free energy, heat capacity, and entropy) in both single-oscillator and lattice contexts.

Methodology
The study employs a rigorous statistical-mechanical framework combining exact integral representations with asymptotic analysis:

Exact Formulation: The canonical partition function for a fixed frequency is integrated against the folded Gaussian probability density function (PDF). This yields exact integral representations for all thermodynamic observables.
Analytic Approximations: Two complementary limits are analyzed:
- Small Variance ( $\sigma \ll \mu$ ): A cumulant expansion is performed around the mean frequency $\mu$ . The free energy, internal energy, and heat capacity are expanded to second order in the variance.
- Zero-Center Limit ( $\mu = 0$ ): The distribution reduces to a half-normal distribution. Low-temperature asymptotic analysis is applied, focusing on the behavior of the Bose-Einstein kernel near $\omega = 0$ where the PDF is finite and non-zero.
Scaling and Lattice Extension: The model is extended to $N$ independent oscillators to test extensivity and self-averaging. Finally, the system is embedded in a $d$ -dimensional disordered lattice, where the single-site frequency distribution is convolved with the van Hove density of states (DOS) of the clean lattice ( $\rho(\omega) \sim \omega^{d-1}$ ).

Key Contributions and Results

Single Oscillator Thermodynamics ( $\mu = 0$ ):
- For a single oscillator with a zero-mean folded Gaussian frequency, the heat capacity $C$ and entropy $S$ do not vanish exponentially at low temperatures. Instead, they exhibit linear behavior: $C \propto T$ and $S \propto T$ .
- This linear scaling arises because the folded Gaussian has a finite, non-zero density at $\omega=0$ ( $p_{FG}(0) > 0$ ), creating a "soft-mode" tail.
- Crucially, the third law of thermodynamics is satisfied: the entropy vanishes as $T \to 0$ , albeit linearly rather than exponentially.
Cumulant Expansion Validity:
- In the regime $\mu \gg \sigma$ , the cumulant expansion provides accurate approximations for the internal energy.
- However, the expansion fails to accurately capture the heat capacity even at moderate $\sigma/\mu$ ratios due to the strong nonlinearity of the heat capacity function with respect to frequency. The relative error for $C$ can reach $\sim 80\%$ for $\sigma/\mu \approx 0.25$ .
Extensivity and Self-Averaging:
- For $N$ independent oscillators, the total heat capacity is extensive.
- Sample-to-sample fluctuations scale as $N^{-1/2}$ , confirming that the system is self-averaging in the thermodynamic limit. The distribution of the heat capacity per oscillator approaches a Gaussian.
Disordered Lattice Behavior:
- When embedded in a $d$ -dimensional lattice, the single-particle linear law ( $C \propto T$ ) is modified by the van Hove factor $\omega^{d-1}$ .
- For $\mu = 0$ , the lattice heat capacity follows a power law $C \propto N T^d$ .
- For $\mu > 0$ , the system exhibits standard Debye-like exponential freeze-out ( $C \propto T^d e^{-\hbar\mu/k_B T}$ ) modified by a disorder prefactor.
Clarification of "Third-Law Violation":
- The paper challenges the notion of a "third-law violation" in disordered phonon systems. It argues that the power-law freeze-out observed in the $\mu=0$ lattice is a spectral property resulting from the absence of an energy gap and the specific $\omega^{d-1}$ phase-space volume, driven by the single-site distribution kink.
- This is distinct from a genuine Lifshitz tail, which requires rare, large-scale spatial fluctuations. The folded Gaussian model demonstrates that a simple kink in the frequency distribution is sufficient to produce soft-mode thermodynamics without invoking large-deviation mechanisms.

Significance and Claims
The author positions the folded Gaussian oscillator as a "minimal benchmark" for studying soft-mode disorder thermodynamics. Its significance lies in its ability to:

Guarantee positivity of frequency naturally (unlike truncated Gaussians).
Isolate the thermodynamic consequences of a non-analytic kink at the origin.
Serve as a "pedagogical null test" for disorder effects, distinguishing between spectral kinks and true Lifshitz tails.

The paper concludes that while the lattice with $\mu=0$ freezes slower (power law) than a clean lattice (exponential), it strictly satisfies the third law ( $S \to 0$ ). The "violation" often cited in literature is reinterpreted as a difference in the rate of freezing (spectral scaling) rather than a residual entropy floor. The work provides exact integral solutions and controlled analytic limits that clarify the interplay between disorder-induced spectral features and thermodynamic phase space.

Technical Summary: Statistical Mechanics of a Quantum Harmonic Oscillator with Folded Gaussian Frequency

More like this