# Fermi energy of electron gas with electrostatic interaction

I have been given the following exam question and am unsure how I would go about solving it:

Consider the case of a one-dimensional metal, consisting of a chain of $$N$$ positive charges $$+q$$ separated by a distance $$2R$$ and immersed in a neutralizing background of electrons with density per unit length $$n_e$$. The electrostatic energy due to the interaction of the electrons with the ion cores and among themselves is: $$\mathcal{E}_{el}=-\frac{\alpha q^{2}}{4\pi \epsilon_{0}}\frac{N}{R}$$

Assuming that the electrons form a non-interacting Fermi gas, calculate the Fermi energy. Write down an expression for the total kinetic energy of the electrons.

Now, ordinarily I would compute the Fermi energy as follows: First write down the density of states for a gas of Fermions:

$$g(\epsilon)=\frac{V m^{3/2}}{\sqrt{2}\pi^{2}\hbar^3}\epsilon^{1/2}$$

Then we have:

$$N = \int_{0}^{\infty}g(\epsilon)n_{F}(\epsilon)\:\mathrm{d}\epsilon = \int_{0}^{\infty}g(\epsilon)\Theta(E_{F}-\epsilon)\:\mathrm{d}\epsilon$$

Where $$\Theta(\epsilon)$$ is the Heaviside step function. We can thus calculate:

$$N=\int_{0}^{E_{F}}g(\epsilon)\:\mathrm{d}\epsilon = \frac{Vm^{3/2}}{\sqrt{2}\pi^{2}\hbar^{3}}\frac{2}{3}\epsilon_{F}^{3/2}$$

So we have the Fermi energy:

$$\epsilon_{F} = \left(\frac{N}{V}\frac{2\sqrt{2}\pi^{2}\hbar^{3}}{3m^{3/2}}\right)^{2/3}=\frac{2\hbar^{2}}{m}\left(\frac{\pi^{2}}{3}\right)^{2/3}n_{e}^{2/3}$$

To calculate the total energy, we have:

$$E = \int_{0}^{E_{F}}\epsilon g(\epsilon)\:\mathrm{d}\epsilon$$

But none of this takes into account the electrostatic energy, so I fear that I have terribly misunderstood something.

Thanks!

Kinetic energy, method 2. Number of filled states \begin{equation} N=\frac{kL}{\pi} \end{equation} where $L=\frac{N}{n_e}$ is the chain length. This gives the wave vector \begin{equation} k=\frac{\pi N}{L}=\pi n_e \end{equation} and thus the kinetic energy \begin{equation} E_{kin}=\frac{\hbar^2}{2m}k^2=\frac{\hbar^2\pi^2n_e^2}{2m} \end{equation}