Quantization rule and harmonics of matter waves.

When modern scientist is recalling the quantization rule for photons, usually the famous E=h*ν is recalled. However, in the original Planck's derivation the more general rule of quantization was assumed: E=n*h*ν [1]. A similar rule of quantization was assumed by Nield Bohr concerning the orbital moment. The values n=2,3,4 .. are responsible for the excited states of the quantum system.
In quantum electrodynamics the quantization rule for electromagnetic field is similar to oscillator:
E=n*h*ν +0.5*h*ν 
At the same time the wavelength of the de Broglie wave has only one value: λ=h/p=h/(m*v)  not λ=n*h/p, where n is a number 1,2,3..., p is the pulse of the particle (applicable for any particle), m - rest mass and v - velocity of the particle (applicable only for non-relativistic case). 
Since the matter waves are not really easy to investigate, it may happened that the de Broglie wavelength also follows the most general rule with the presence of the harmonics:
λ=n*h/p,
But they were simply overlooked and careful experiment would be necessary to discover them.
For de Broglie wave it is would not be easy to find presence of such harmonics, because during the interference experiment they would generate the maxima and minima, which coincide with maxima and minima of the main matter wave. If harmonic is present in the miniscule amount it will lead to some hardly observable effect. 
Let's consider for example the case of only one added harmonic. Let the first harmonic has maximum of 1 and decays as exp(-0.1*m), where m is the interference band number, m=0,1,2,3,4... In this case the amplitudes of the interference bands would be: Io=1, exp(-0.1), exp(-0.2), exp(-0.3), exp(-0.4)=1, 0.905, 0.819, 0.741, 0.670...
Second harmonic of de Broglie wave would have the amplitude of 0.01 and decays with distance away from the center according to the same law exp(-0.1*m), here m is the interference band for the second harmonic, which would coincide with a certain band for first harmonic (because wavelength is exactly 2 times larger). The amplitude would be I=0.01, 0.01*exp(-0.1), 0.01*exp(-0.2)… =0.01, 0.00905, 0.00819 ...
The sum of the amplitudes would be:
1.01; exp(-0.1); exp(-0.2)+0.01*exp(-0.1); exp(-0.3); exp(-0.4)+0.01*exp(-0.2); …..=
1.01; 0.905; 0.828; 0.741; 0.679; …..
In order to distinguish the case of the one and multiple harmonics the ratio of the amplitudes of the consecutive bands may be calculated: bamd1/band0; band2/band1, band3/band2.....
For exactly one de Broglie wavelength that would be monotonic function:
0.905, 0.905, 0.905 …. (constant in this example, because the chosen decay function was exponential)
For the sum of harmonics it would be: ratios are:
0.896; 0.915; 0.895; 0.916 …. - non-monotonic function, the superposition of monotonic function and ao*Cos(pi*m), where m is the interference band number.
The third and higher harmonics will add more "waviness" to the smooth function.
How to estimate the amplitude of the harmonics in matter wave? The idea of evaluation is inferred from the reciprocity principle: the particle is both matter and wave [2]. Assuming the matter wave is something real (similar to a photon, but permanently "attached" to the particle), the probability of the excitation of the second energy level would be similar to the idea proposed by Plank [1]:
population of each next level would follow Boltzmann rule [3]:
-log(N_i/N)~E_i/kT
But what is the expected energy of the initial de Broglie wave? Hypothesizing that  the particle is both matter and wave it is possible to speculate about this value.
Since de Broglie wave is "attached" to the particle, the only velocity it may have is equals to the velocity of particle v. From the general rule connecting velocity, wavelength and frequency of the wave it follows:

v=λ*ν or λ=v/ν
here v is the velocity of the particle, λ is the wavelength, ν is the frequency of the wave.
Substituting λ  into the formula for non-relativistic de Broglie wave:
λ=h/(m*v)   and v/ν=h/(m*v)   
which may be transformed as follows:
mv²=hν  or mv²/2=hν/2
For non-relativistic particle the energy of de Broglie wave can not be larger than the full kinetic energy of the particles and quantized as oscillator (quite reasonable idea, because the zero energy of electromagnetic field has the same value). Assuming the next harmonic will have the energy according to the Planks rule (or quantum electrodynamic rule, similar to oscillator), the difference in energy between two levels for de Broglie wave would be double the kinetic energy of the non-relativistic particle (for relativistic particle the quantization rule is simpler and coincides with the quantization rule for photons).
Than the ratio of the amplitude of the second harmonic of de Broglie wave to the initial amplitude would be equal (from Boltzmann rule): 

I/I_o=exp(-2E_k/kT)
where Ek is the kinetic energy of the non-relativistic particle. For example for electron with possible to reach energy of 0.1 eV, observed at the room temperature (300 K, at this temperature de Broglie wave is still well resolved since 0.1 eV > kT), the ratio would be:

I/I_o=exp(-2E_k/kT)=4.4*10exp(-4)
Which is small, but at numerous averaging is possible to reach and to discover.


References.
1. http://web.phys.ntnu.no/~stovneng/TFY4165_2013/BlackbodyRadiation.pdf
2. https://tipikin.blogspot.com/2019/09/the-possible-way-to-search-for-new.html
3. https://en.wikipedia.org/wiki/Maxwell–Boltzmann_distribution