Lifshitz–Kosevich Quantum Oscillation Calculator

JJ Ben-Joseph headshot JJ Ben-Joseph

Enter parameters to compute.

Quantum Oscillations in Metals and the Lifshitz–Kosevich Formula

Quantum oscillations provide a powerful window into the Fermi surface of metals and semimetals. When a clean metallic system is placed in a strong magnetic field, the energy levels of conduction electrons collapse into discrete Landau levels separated by the cyclotron energy. As the field varies, these levels sweep past the Fermi energy and cause oscillations in thermodynamic and transport properties, a phenomenon known as the de Haas–van Alphen effect (for magnetization) and the Shubnikov–de Haas effect (for resistivity). The Lifshitz–Kosevich (LK) formula encapsulates how the amplitude of these oscillations depends on temperature, impurity scattering, magnetic field strength, and the geometry of the Fermi surface. Our calculator implements the LK expression for the amplitude reduction factors, allowing researchers to estimate how observable the oscillations will be under given conditions.

At the heart of the LK theory lies the assumption of a quasi-two-dimensional slice of the Fermi surface perpendicular to the applied field. The extremal cross-sectional area A_F determines the basic oscillation frequency via F = \frac{\hbar}{2\pi e} A_F, measured in tesla. The oscillatory part of the magnetization or resistivity typically includes a factor \cos\left(2\pi\left(\frac{F}{B} - \gamma\right)\right), where the phase \gamma depends on Berry phase and dimensionality. However, finite temperature and disorder damp the amplitude, and the LK formula expresses this suppression as a product of separate factors. The thermal factor is R_T = \frac{X}{\sinh X} with X = \frac{2\pi^2 k_B T m^*}{\hbar e B}. The Dingle factor accounting for impurity scattering is R_D = \exp\left(-\frac{2\pi^2 k_B T_D m^*}{\hbar e B}\right). Multiplying these together yields the overall damping. The calculator evaluates these expressions numerically using user-specified values of B, T, effective mass m^*, and Dingle temperature T_D. It also computes the oscillation period in inverse field units, \Delta(1/B) = 1/F, revealing how the oscillations space out in reciprocal magnetic field.

The thermal reduction factor arises from the smearing of the Fermi–Dirac distribution at finite temperature. In a magnetic field, each Landau level is broadened by the derivative of the Fermi function, and as temperature increases, the sharpness of the occupancy step at the Fermi energy diminishes. The factor R_T = X/\sinh X captures this effect with remarkable precision. Here, X = 14.69 \frac{m^* T}{B} when temperature is measured in kelvin, field in tesla, and effective mass in units of the bare electron mass. As T \to 0, X \to 0 and the factor tends to unity; at high temperature, the hyperbolic sine dominates, causing the amplitude to fall exponentially. Users can exploit the calculator to determine what temperatures are necessary to resolve oscillations for a given field and effective mass: heavier quasiparticles (larger m^*) suffer stronger suppression and thus require lower temperatures or higher fields.

The Dingle factor reflects impurity scattering and finite quasiparticle lifetimes. Introducing a lifetime broadens the Landau levels by an energy γ = ħ/(2τ), equivalently characterized by the Dingle temperature T_D = \frac{\hbar}{2\pi k_B \tau}. The exponential factor R_D = e^{-X_D} with X_D = 14.69 \frac{m^* T_D}{B} conveys that cleaner samples (smaller T_D) exhibit stronger oscillations. By adjusting the Dingle temperature in the input, experimentalists can assess how much improvement in sample quality is needed for a detectable signal. Conversely, fitting observed amplitudes as a function of field can yield T_D, providing a sensitive probe of disorder.

The frequency F itself encodes the extremal Fermi-surface area. In quasi-two-dimensional compounds such as the cuprate superconductors or organic conductors, the presence of multiple orbits leads to several frequencies, each with its own amplitude and phase. While our calculator focuses on a single frequency, users can apply it iteratively to different F values to analyze multi-orbit systems. The period \Delta(1/B) = 1/F indicates how widely spaced the oscillations are in inverse field: for F = 100\,\text{T}, the period is 0.01 T⁻¹, meaning that sweeping 1/B by this amount moves through one oscillation cycle. Experiments often plot magnetization versus 1/B to reveal the periodic structure clearly.

Quantitative interpretation of oscillation amplitudes also requires accounting for spin splitting. The LK formula can be extended with a spin-splitting factor R_S = \cos\left(\pi g m^*/2 m_e\right), where g is the effective g-factor. When the cosine vanishes, spin-up and spin-down contributions cancel, leading to “spin zeros.” While our calculator omits spin for simplicity, the amplitude output can be interpreted as the baseline before spin modulation. Researchers interested in g-factor measurements can combine the damping factors computed here with additional analysis to predict spin-zero conditions.

Historically, the LK theory marked a triumph of quantum many-body physics. Developed in the 1950s by I. M. Lifshitz and A. M. Kosevich, it unified semiclassical quantization with statistical mechanics, providing a robust tool to map Fermi surfaces long before angle-resolved photoemission became feasible. The method has enjoyed a renaissance in the study of topological materials, where Berry phases extracted from oscillation phases reveal nontrivial electronic structure. The thermal and Dingle factors remain central: only when R_T R_D is sufficiently large can phase analysis proceed accurately.

To demonstrate the influence of various parameters, the table below presents sample damping factors computed by the calculator for typical values encountered in experiments:

B (T)T (K)m*/mₑT_D (K)R_TR_DR_T R_D
101150.860.480.41
302230.550.760.42

These examples illustrate that high magnetic fields and low temperatures enhance the amplitude, but heavy effective masses and large Dingle temperatures suppress it. The interplay often guides experimental design: if the calculated R_T R_D is too small, one must either increase the field, reduce the temperature, or improve sample quality. Conversely, observing sizable oscillations can serve as a benchmark for the success of material synthesis and purification efforts.

Beyond simple metals, quantum oscillations have been observed in Dirac and Weyl semimetals, high-temperature superconductors, and two-dimensional electron gases. In these systems, the LK formula sometimes requires modifications to account for non-parabolic dispersion or strong correlations, yet the central role of thermal and disorder damping persists. The calculator thus provides a first approximation that can be refined with system-specific factors. For instance, in Dirac materials, the effective mass entering R_T is energy-dependent, reflecting the linear dispersion. Still, the expressions implemented here offer a baseline for understanding trends.

Finally, it is worth mentioning the practical challenge of measuring oscillations: high fields often necessitate pulsed magnets, where the field changes rapidly, complicating data acquisition. In such environments, knowing the expected amplitude ahead of time is invaluable for optimizing measurement protocols. The Lifshitz–Kosevich calculator empowers researchers to perform these estimates quickly, fostering efficient experiment planning. Whether exploring the Fermi surface of a newly synthesized material or teaching students about quantum oscillation phenomenology, this tool bridges theory and practice.

In summary, the Lifshitz–Kosevich formula elegantly captures how temperature, disorder, magnetic field, and quasiparticle mass determine the visibility of quantum oscillations. The calculator implemented here computes the thermal factor R_T, the Dingle factor R_D, and their product, while also reporting the oscillation period in inverse field units. By experimenting with different inputs, users can gain intuition about the delicate balance of parameters required for successful de Haas–van Alphen measurements and related phenomena.

Related Calculators

Van de Graaff Generator Calculator

Estimate dome capacitance, maximum voltage, charging current, and charge time for a Van de Graaff electrostatic generator.

Van de Graaff calculator electrostatic generator high voltage capacitance charging time

Camper Van Conversion Budget Calculator - Plan Your Adventure Rig

Estimate the total cost of converting a van into a camper including materials, labor, and unexpected expenses.

van conversion calculator camper budget estimator travel van cost

De Broglie Wavelength Calculator - Matter Waves

Compute the de Broglie wavelength of any particle by entering its mass and velocity. Explore wave-particle duality and the quantum realm.

de Broglie wavelength calculator matter wave quantum mechanics