My statistical mech. book's got some info about it. If you take the quantum mechanical ideal gas where the wavefunctions are subject to perdodic boundary conditions (period L) the single-particle eigenstates are plane-waves of momentum [itex](p_x,p_y,p_z)=\frac{h}{L}(n_x,n_y,n_z)[/itex] where [itex]n_x,n_y,n_z[/itex] are integers the average occupation of one of the eigenstates is
[tex]\bar N_p = \frac{1}{\exp(\alpha+p^2/2mkT) \pm 1}[/tex]
where the + is taken for fermions and the minus for bosons. The thermal wavelength comes into play in the paragraph of the classical limit.
If [itex]e^\alpha gg 1[/itex] then the [itex]\pm 1[/itex] in the above expression is neglible, giving:
[tex]\bar N_p \approx e^{-\alpha}e^{-p^2/2mkT}[/tex]
which is the classical Maxwell distribution.
One way to interpret [itex]e^\alpha \gg 1[/itex] is to note that it implies that the average occupation of every momentum state is much less than one. In order to express [itex]e^\alpha[/itex] in terms of the macroscopic parameters of the system, let's assume that it is valid and evaluate the normalization sum for the qm-distribution.
[tex]\sum_{\vec p}\bar N_p = \sum_{\vec p} \approx e^{-\alpha}e^{-p^2/2mkT}=N[/tex]
where N is the # of particles.The sum is taken over all momentum eigenstate values. The density of momentum e.v.'s in momentum space is [itex](L/h)^3=V/h^3[/itex], where V is the volume of the system, the sum over the discrete value of p may be converted to an integral over a continuous vector variable.
[tex]\frac{Ve^{-\alpha}}{h^3}\int \exp(-p^2/2mkT) d^3\vec p = N[/tex]
This is a standard Guassian integral (the product of three actually, one for each component):
[tex]\int \exp(-p^2/2mkT) d^3\vec p=(2\pi mkT)^{3/2}[/tex]
Now the book defines a \emph{thermal deBroglie wavelength} [itex]\lambda[/itex] by:
[tex]\lambda \equiv \frac{h}{\sqrt{2\pi mkT}}[/tex]
then the normalization gives:
[tex]e^{-\alpha}=(N/V)\lambda^3[/tex]
Therefore, [itex]\e^{-\alpha}[/itex] is equal to the average number of particles in volume [itex]\lambda^3[/itex]. If that number is much less than one, the classical approximation is justified. [itex]\lambda[/itex] is the deBroglie wavelength of a particle with energy [itex]\pi kT[/itex], which is about twice the average energy of the particles in a gas at temperature T.
Well, I hope that was at least a little bit enlightening. I guess the reason for the definition should be gotten from:
[tex]\lambda = \frac{h}{p}[/tex]
where
[tex]p=\int \exp(-p^2/2mkT) d p[/tex]