Correction: Camiola, V.D., et al. Equilibrium Wigner Function for Fermions and Bosons in the Case of a General Energy Dispersion Relation. Entropy 2020, 22, 1023

In Section 5 of Equilibrium Wigner Function for Fermions and Bosons in the Case of a General Energy Dispersion Relation [...].

In Section 5 of Equilibrium Wigner Function for Fermions and Bosons in the Case of a General Energy Dispersion Relation (Entropy 2020, 22, 1023 [1]), the Laplace transform was used in the wrong way to find the solution of Equation (13); therefore, the derived Equation (14) is not correct. Here, we rewrite Section 5, the determination of examples of the solution by adopting suitable asymptotic expansions. This correction does not modify the main conclusion of the article, except the fact that we find additional solutions that do not represent only a modulation of the solutions in the non-degenerate case.
The authors apologize for the oversight and hope that it has not had affected the work of other colleagues.
The new section is as follows.

Particular Cases
Let us consider the homogeneous case. We have Φ(x) = Φ = constant, and the Wigner function does not depend on x, that is w = w(k, β). Moreover, we assume that the first Brillouin zone is finite, that is, w eq (k, β) is zero outside a compact set B, which is symmetric with respect to the origin. Under these assumption, Equation (10) reads Solving this equation is a very difficult task, and here, we propose two possible ways to find an approximate solution that satisfies the condition lim β→0 + w eq (k, β) = 1 (see [24]). We assume the ansatz with n ∈ N and λ : B → R.
Substituting in (13) and considering the lower order terms, one obtains where µ(B) is the measure of the first Brillouin zone. The equilibrium Wigner function approximated at the first order reads The other terms λ m (k) with m ≥ 2 can be found with an iterative procedure. Another analytical result can be found around the origin of the first Brillouin zone, that is, for |k| 1. About the dispersion relation, the natural assumption ε(k) = ak 2 + o(|k| 2 ), with a ∈ R + , is made for |k| 1. We look for solutions according to the ansatz w eq (k, β) = γ(k, β)e −βε(k) with γ(k, β), a slowly varying function with respect to k: By substituting in (13), one obtains At the lowest order, we have Due to the symmetry of B B e −2aβk | 2 k · k dk = 0 and, therefore, γ 0 (β) satisfies the equation for which the solution is Therefore, the equilibrium Wigner function around the center of the first Brillouin zone is In order to present some numerical results, let us consider a tridimensional gas of non-interacting particles with a quadratic dispersion relation ε(k) = ak 2 , a ∈ R + . To analytically evaluate the term α(β), B is extended to R 3 , obtaining which, inserted into (17), gives In Figure 1, the Wigner function (18) is plotted versus energy for several values of the parameters β and φ − φ F . For comparison, the Maxwell-Boltzmann distribution w * = e −β(ε+φ−φ F ) is also shown. It is important to observe that, at high temperatures, that is, low β's, the equilibrium Wigner functions for Bosons and Fermions are both close to the Maxwell-Boltzmann distribution; therefore, Bosons and Fermions tend to have the same behaviour for β → 0 + . Moreover, in the plotted cases, the equilibrium Wigner function is positive. At low temperatures (high values of β), the behaviour strongly depends on φ − φ F . When φ − φ F = 0.5, the equilibrium Wigner functions for Bosons and Fermions are still very close and positive. If φ − φ F = −0.5, the equilibrium Wigner functions of Bosons and Fermions and the Maxwell-Boltzmann distribution have relevant differences: the most remarkable one is that the Fermions Wigner function is negative while the others are positive.
As a last remark, although apparently simple, the case of constant potential has physically relevant applications. If one considers the transport of phonons in a crystal lattice without any mechanical deformation, they do not undergo any external field but have a dispersion relation that is not usually quadratic. For example, acoustic phonons have a linear dispersion relation near the center of the first Brillouin zone (the Debye approximation), that is ε(p) = c sh |p| where c s is the sound speed.