AbstractWe study the role of disorder in the vibration spectra of molecules and atoms in solids. This disorder may be described phenomenologically by a fractional generalization of ordinary quantum-mechanical oscillator problem. To be specific, this is accomplished by the introduction of a so-called fractional Laplacian (Riesz fractional derivative) to the Scrödinger equation with three-dimensional (3D) quadratic potential. To solve the obtained 3D spectral problem, we pass to the momentum space, where the problem simplifies greatly as fractional Laplacian becomes simply $k^μ $ k μ , k is a modulus of the momentum vector and $μ $ μ is Lévy index, characterizing the degree of disorder. In this case, $μ → 0$ μ → 0 corresponds to the strongest disorder, while $μ → 2$ μ → 2 to the weakest so that the case $μ =2$ μ = 2 corresponds to “ordinary” (i.e. that without fractional derivatives) 3D quantum harmonic oscillator. Combining analytical (variational) and numerical methods, we have shown that in the fractional (disordered) 3D oscillator problem, the famous orbital momentum degeneracy is lifted so that its energy starts to depend on orbital quantum number l. These features can have a strong impact on the physical properties of many solids, ranging from multiferroics to oxide heterostructures, which, in turn, are usable in modern microelectronic devices.