First introduced by M. Pollak,^[1] the Coulomb gap is a soft gap in the single-particle density of states (DOS) of a system of interacting localized electrons. Due to the long-range Coulomb interactions, the single-particle DOS vanishes at the chemical potential, at low enough temperatures, such that thermal excitations do not wash out the gap.

At zero temperature, a classical treatment of a system gives an upper bound for the DOS near the Fermi energy, first suggested by Efros and Shklovskii.^[2] The argument is as follows: Let us look at the ground state configuration of the system. Defining ${\displaystyle E_{i}}$ as the energy of an electron at site ${\displaystyle i}$, due to the disorder and the Coulomb interaction with all other electrons (we define this both for occupied and unoccupied sites), it is easy to see that the energy needed to move an electron from an occupied site ${\displaystyle i}$ to an unoccupied site ${\displaystyle j}$ is given by the expression:

${\displaystyle \Delta E=E_{j}-E_{i}-e^{2}/r_{ij}}$.

The subtraction of the last term accounts for the fact that ${\displaystyle E_{j}}$ contains a term due to the interaction with the electron present at site ${\displaystyle i}$, but after moving the electron this term should not be considered. It is easy to see from this that there exists an energy ${\displaystyle E_{f}}$ such that all sites with energies above it are empty, and below it are full (this is the Fermi energy, but since we are dealing with a system with interactions it is not obvious a-priori that it is still well-defined). Assume we have a finite single-particle DOS at the Fermi energy, ${\displaystyle g(E_{f})}$. For every possible transfer of an electron from an occupied site ${\displaystyle i}$ to an unoccupied site ${\displaystyle j}$, the energy invested should be positive, since we are assuming we are in the ground state of the system, i.e., ${\displaystyle \Delta E>=0}$. Assuming we have a large system, consider all the sites with energies in the interval ${\displaystyle [E_{f}-\epsilon ,E_{f}+\epsilon ].}$ The number of these, by assumption, is ${\displaystyle N=2\epsilon g(E_{f}).}$ As explained, ${\displaystyle N/2}$ of these would be occupied, and the others unoccupied. Of all pairs of occupied and unoccupied sites, let us choose the one where the two are closest to each other. If we assume the sites are randomly distributed in space, we find that the distance between these two sites is of order: ${\displaystyle R\sim (N/V)^{-1/d}}$, where ${\displaystyle d}$ is the dimension of space. Plugging the expression for ${\displaystyle N}$ into the previous equation, we obtain the inequality: ${\displaystyle E_{j}-E_{i}-Ce^{2}(\epsilon g(E_{f})/V)^{1/d}>0}$ where ${\displaystyle C}$ is a coefficient of order unity. Since ${\displaystyle E_{j}-E_{i}<2\epsilon }$, this inequality will necessarily be violated for small enough ${\displaystyle \epsilon }$. Hence, assuming a finite DOS at ${\displaystyle E_{f}}$ led to a contradiction. Repeating the above calculation under the assumption that the DOS near ${\displaystyle E_{f}}$ is proportional to ${\displaystyle (E-E_{f})^{\alpha }}$ shows that ${\displaystyle \alpha >=d-1}$. This is an upper bound for the Coulomb gap. Efros^[3] considered single electron excitations, and obtained an integro-differential equation for the DOS, showing the Coulomb gap in fact follows the above equation (i.e., the upper bound is a tight bound).

Other treatments of the problem include a mean-field numerical approach,^[4] as well as more recent treatments such as,^[5] also verifying the upper bound suggested above is a tight bound. Many Monte Carlo simulations were also performed,^[6][7] some of them in disagreement with the result quoted above. Few works deal with the quantum aspect of the problem.^[8] Classical Coulomb gap in clean system without disorder is well captured within Extended Dynamical Mean Field Theory (EDMFT) supported by Metropolis Monte Carlo simulations. ^[9]

Direct experimental confirmation of the gap has been done via tunneling experiments, which probed the single-particle DOS in two and three dimensions.^[10][11] The experiments clearly showed a linear gap in two dimensions, and a parabolic gap in three dimensions. Another experimental consequence of the Coulomb gap is found in the conductivity of samples in the localized regime. The existence of a gap in the spectrum of excitations would result in a lowered conductivity than that predicted by Mott variable-range hopping. If one uses the analytical expression of the single-particle DOS in the Mott derivation, a universal ${\displaystyle e^{-1/T^{1/2}}}$ dependence is obtained, for any dimension.^[12] The observation of this is expected to occur below a certain temperature, such that the optimal energy of hopping would be smaller than the width of the Coulomb gap. The transition from Mott to so-called Efros–Shklovskii variable-range hopping has been observed experimentally for various systems.^[13] Nevertheless, no rigorous derivation of the Efros–Shklovskii conductivity formula has been put forth, and in some experiments ${\displaystyle e^{-1/T^{\alpha }}}$ behavior is observed, with a value of ${\displaystyle \alpha }$ that fits neither the Mott nor the Efros–Shklovskii theories.

• Coulomb's law