Chemical potential vs pressure and temperature; difficulty with Fermi gases

In summary, the chemical potential decreases with increasing temperature but increases with increasing pressure in 1D Fermi gas.
  • #1
jfizzix
Science Advisor
Insights Author
Gold Member
757
355
For any system where the thermodynamic limit exists, we know that the internal energy U, the entropy σ, the total particle number N and the total volume V are all extensive. Because of this, we know that the Euler relation holds true

U = -PV + [itex]\tau[/itex]σ + [itex]\mu[/itex]N

and that the chemical potential is just the Gibbs free energy per particle.

G = [itex]\mu[/itex] N

We can use the differential relation of the Gibbs free energy

dG = V dP -σ d[itex]\tau[/itex] + [itex]\mu[/itex] dN

to find expressions for the derivatives of the chemical potential.
d[itex]\mu[/itex]/d[itex]\tau[/itex] = 1/N * dG/d[itex]\tau[/itex] =-σ/N
d[itex]\mu[/itex]/dP = 1/N * dG/dP = V/N

Since N, [itex]\sigma[/itex], and V are all positive, it appears that the chemical potential must increase with pressure at constant temperature, and decrease with temperature at constant pressure.

My question is this:
the 1D fermi gas is a case where the chemical potential actually increases with temperature for small temperatures, though we would think in all cases the chemical potential should be decreasing with temperature.

How do we reconcile this result with the previous assumptions about the chemical potential?

-James
 
Science news on Phys.org
  • #2
jfizzix said:
My question is this:
the 1D fermi gas is a case where the chemical potential actually increases with temperature for small temperatures, though we would think in all cases the chemical potential should be decreasing with temperature.
I don't have this stuff in mind. Could you please for our convenience write down some formula that shows this?
I suspect you are considering a Fermi gas with fixed particle number instead of fixed pressure or equivalently U.
 
  • #3
Excuse me,but I don't very understand how it is done with this:
dμ/dτ = 1/N * dG/dτ =-σ/N
dμ/dP = 1/N * dG/dP = V/N

I am a little confused.
 
  • #4
He is using e.g.
[itex](\partial \mu/\partial T)_{P,N}=-(\partial S/\partial N)_{T,P}=-S/N[/itex]
The first equality is a Maxwell relation, the second one follows from S(N,p,T) being an extensive function whence it must be proportional to N.
 
  • #5
DrDu said:
I don't have this stuff in mind. Could you please for our convenience write down some formula that shows this?
I suspect you are considering a Fermi gas with fixed particle number instead of fixed pressure or equivalently U.


For the 1D (or 2D or 3D) Fermi gas we can calculate the chemical potential in the following way. In a nutshell, we say that the number of particles is constant whether the temperature is zero or not.

First, we can find the number of particles at zero temperature in terms of the Fermi energy. Starting at the ground state, we put in fermions at each successive orbital until all N fermions have been placed. The highest energy (that of the final fermion) is the Fermi energy.

Next we can calculate the number of particles at nonzero temperature by integrating over the density of states. For a nonzero temperature, the integral is harder to evaluate, but can be carried out for small temperatures with what is known as the Sommerfeld expansion/approximation.

Setting these two expressions for the number of particles equal to each other, we can calculate the chemical potential in terms of the Fermi energy.

for small temperatures the chemical potential of the 1D Fermi gas increases in the following way with temperature:
[itex]\frac{\mu}{E_{F}}[/itex] ≈ 1 + [itex]\frac{\pi^{2}}{12}[/itex]*([itex]\frac{k_{b}T}{E_{F}})^{2}[/itex]

In the 2D case, the chemical potential is just the Fermi energy at low temperatures, and in the 3D case, the chemical decreases with temperature in the same way that it increases for the 1D gas. Kittel and Kroemer has a nice derivation, as well as Schroeder's book.
 
Last edited:
  • Like
Likes sourav kundu
  • #6
DrDu said:
I don't have this stuff in mind. Could you please for our convenience write down some formula that shows this?
I suspect you are considering a Fermi gas with fixed particle number instead of fixed pressure or equivalently U.

The natural variables for the Gibbs free energy are G(P,T,N), just as for the internal energy
U=U(S,V,N)
H=H(S,P,N) (for the enthalpy)
F=F(T,V,N) (for the Helmholtz free energy)

Since P and T are intensive, and G is extensive, we can write G as N times some other function of P and T.
Taking the derivative of G with respect to N, we get that [itex]\mu[/itex] is that same function of P and T.
 
  • #7
DrDu said:
He is using e.g.
[itex](\partial \mu/\partial T)_{P,N}=-(\partial S/\partial N)_{T,P}=-S/N[/itex]
The first equality is a Maxwell relation, the second one follows from S(N,p,T) being an extensive function whence it must be proportional to N.

Now that I think about it, the maxwell relation here says it all.

in 1-D, the density of states decreases with increasing energy as [itex]1/\sqrt{\epsilon}[/itex]. This means the number of accessible states decreases with increasing particle number (and hence increasing Fermi energy).

As such, the entropy must decrease with the addition of more fermions, which means that the chemical potential increases with increasing temperature.

I'm satisfied as to why the chemical potential decreases using the maxwell relations, but why then can we not say G = [itex]\mu[/itex]*N.

Thanks to all who helped,

-James
 
  • #8
jfizzix said:
, but why then can we not say G = [itex]\mu[/itex]*N.

But this equation holds true!
 
  • #9
jfizzix said:
For the 1D (or 2D or 3D)
[itex]\frac{\mu}{E_{F}}[/itex] ≈ 1 + [itex]\frac{\pi^{2}}{12}[/itex]*([itex]\frac{k_{b}T}{E_{F}})^{2}[/itex]

What you need is mu(P,T), so you should eliminate E_F in terms of P and T.
 
  • #10
Ok,
I think I understand it now:
To lowest order in the Sommerfeld expansion independent of dimension [itex]\mu(T,P)=E_F(T,P)[/itex]
so
[itex] \partial \mu /\partial T|_P=\partial E_F /\partial T|_P=- (\partial P/\partial T)_{E_F}(\partial E_F/\partial P)_T[/itex]
As the Fermi energy is a monotonously increasing function of the density, only, in the first derivative on the right we can alternatively keep the density constant. But P always increases with increasing temperature at constant density, so this term is positive. Considering the second derivative, the density (and E_F) increases with p at constant T, so this term is also positive. So to lowest order, the change of mu with T is always negative as it should be. The small higher order corrections will not change this.

Or stated differently, when T is increasing, density or E_F has to disminish for P to remain constant.
 
Last edited:
  • #11
DrDu said:
Ok,
I think I understand it now:
To lowest order in the Sommerfeld expansion independent of dimension [itex]\mu(T,P)=E_F(T,P)[/itex]
so
[itex] \partial \mu /\partial T|_P=\partial E_F /\partial T|_P=- (\partial P/\partial T)_{E_F}(\partial E_F/\partial P)_T[/itex]
As the Fermi energy is a monotonously increasing function of the density, only, in the first derivative on the right we can alternatively keep the density constant. But P always increases with increasing temperature at constant density, so this term is positive. Considering the second derivative, the density (and E_F) increases with p at constant T, so this term is also positive. So to lowest order, the change of mu with T is always negative as it should be. The small higher order corrections will not change this.

Or stated differently, when T is increasing, density or E_F has to diminish for P to remain constant.




The Fermi energy [itex]E_{f}[/itex] does not depend explicitly on temperature; only on particle number [itex]N[/itex] and the space in which the particles are contained (i.e. in 1-D, the length; in 2D, the area; in 3D, the volume).

One can express the Fermi energy in terms of N and length and get the chemical potential in terms of N, length, and temperature.
Integrating this with respect to N gives you the Helmholtz free energy [itex]F(T,V,N)[/itex] since [itex]μ(T,V,N)=(∂F/∂N)_{T,V}[/itex]. we can then take the derivative of F with respect to volume to get the pressure [itex]P(T,V,N)=-(∂F/∂V)_{T,N}[/itex]. We can also take the derivative of F with respect to temperature to get the entropy [itex]S(T,V,N)=-(∂F/∂T)_{V,N}[/itex]. With the state equation of the pressure, we can substitute this into our expression for the entropy to get the entropy in terms of N,P, and T.
If you take the derivative of this entropy with respect to N, you find it is negative.. which means that where
[itex](∂μ/∂T)_{P,N}=−(∂S/∂N)_{T,P}[/itex]

the chemical potential must increase with temperature for the 1D degenerate Fermi gas.

Now, in 2D and 3D (taking higher order terms) the chemical potential decreases with temperature, but the question persistently remains...

How do we say that G is not μ*N? Is G no longer extensive for a 1D Fermi gas?

Thanks again,

-James
 
  • #12
jfizzix said:
The Fermi energy [itex]E_{f}[/itex] does not depend explicitly on temperature; only on particle number [itex]N[/itex] and the space in which the particles are contained (i.e. in 1-D, the length; in 2D, the area; in 3D, the volume).
That is not generally true but depends on the parametrization. E_F can be written as a function of only N/L but this is not what is needed in your question as you are taking the derivative at constant P and N and not at constant L and N.
Think of the ideal gas equation which is also the limit of the Fermi equation of state at low densities or high temperature: N/V=P/RT.
 
  • #13
DrDu said:
That is not generally true but depends on the parametrization. E_F can be written as a function of only N/L but this is not what is needed in your question as you are taking the derivative at constant P and N and not at constant L and N.
Think of the ideal gas equation which is also the limit of the Fermi equation of state at low densities or high temperature: N/V=P/RT.

Point taken..

We start with the Fermi energy of a 1D Fermi gas.

[itex]E_{f} = (\frac{\hbar^2 \pi^2}{8m})(\frac{N}{L})^2 = \alpha_{1} (\frac{N}{L})^2[/itex]

and the formula for the chemical potential at low temperature where the Fermi energy has already been substituted

[itex]\mu(\tau , L, N) ≈ \alpha_{1} (\frac{N}{L})^2 + \frac{\pi^2 \tau^2}{12 \alpha_{1}}(\frac{L}{N})^2[/itex]

We integrate the chemical potential with respect to N from 0 to N to obtain the Helmholtz free energy F. Recall that in 1D...
[itex]dF = - P dL - \sigma d\tau + \mu dN[/itex]
so that

[itex]F( \tau, L, N) = \int_{0}^{N}{\mu (\tau, L, N) dN} = \frac{\alpha_{1}}{3}\frac{N^3}{L^2} - \frac{\pi^2 \tau^2}{12 \alpha_{1}} \frac{L^2}{N}[/itex]

Now we take the derivative of the Helmholtz free energy with respect to length to find the pressure in 1D

[itex]P = -(\frac{\partial F}{\partial L})_{\tau,N} = \frac{2 \alpha_{1}}{3}\frac{N^3}{L^3} + \frac{\pi^2 \tau^2}{6 \alpha_{1}} \frac{L}{N} ≈ \frac{2 \alpha_{1}}{3}\frac{N^3}{L^3} [/itex]

Substituting this expression for the pressure into our original expression for the chemical potential, we find

[itex] \mu(\tau , P, N) ≈ \alpha_{1} (\frac{3 P}{2 \alpha_{1}})^{2/3}+ \frac{\pi^2 \tau^2}{12 \alpha_{1}}(\frac{2 \alpha_{1}}{3 P})^{2/3}[/itex]

In short, for small temperatures, and in one dimension, the chemical potential does increase with temperature at constant pressure and particle number. I could have made a mistake with the math, but I think it's all there.

-James
 
  • #14
jfizzix said:
[itex]F( \tau, L, N) = \int_{0}^{N}{\mu (\tau, L, N) dN} = \frac{\alpha_{1}}{3}\frac{N^3}{L^2} - \frac{\pi^2 \tau^2}{12 \alpha_{1}} \frac{L^2}{N}[/itex]

Methinks you deliberately dropped an infinity from the second term at N=0 on integration!
Obviously what is happening here is that at small N and fixed L and T the Sommerfeld expansion breaks down.
 
Last edited:
  • #15
DrDu said:
Methinks you deliberately dropped an infinity from the second term at N=0 on integration!
Obviously what is happening here is that at small N and fixed L and T the Sommerfeld expansion breaks down.

Never attribute to malice, that which can be attributed to stupidity.

I actually forgot that infinity there.. I'll have to get back to you on how to do it right.

-James
 
  • #16
You really should calculate the Sommerfeld expansion of p in terms of E_F and T.
 
  • #17
DrDu said:
You really should calculate the Sommerfeld expansion of p in terms of E_F and T.

I think you're right.

I just substitute the Fermi pressure along with factors of N and [itex]\tau[/itex] into the expression, and that should do the trick..

more to follow...

-James
 
  • #18
So.. in 1D,

The fermi energy [itex]E_{f}[/itex] can be expressed in terms of the Fermi pressure [itex]P_{f}[/itex] using the expression for the internal energy [itex]U[/itex] at absolute zero. If we differentiate the internal energy [itex]U[/itex] with respect to the length of the 1D box [itex]L[/itex] and add a negative sign, we have an expression for the pressure.

[itex]P = -(\frac{\partial U}{\partial L})_{\sigma,N}[/itex]

In 1D at absolute zero...

[itex]U = \frac{1}{3} N E_{f}[/itex] where [itex]E_{f} = \frac{\hbar^{2}}{8 m}(\frac{N \pi}{L})^{2}[/itex]

From these expressions, we can express the Fermi energy in terms of the pressure at absolute zero.

[itex]E_{f} = (\frac{3 \pi}{2})^{\frac{2}{3}}(\frac{\hbar^{2}}{8 m})^{\frac{1}{3}} P^{\frac{2}{3}} = C_{1} P^{\frac{2}{3}}[/itex]

Then the expression for the chemical potential becomes

[itex]\mu ≈ C_{1} P^{\frac{2}{3}} + \frac{\pi^{2}}{12} \frac{\tau^{2}}{C_{1} P^{\frac{2}{3}}}[/itex]

With this as our form for the chemical potential, we see that it does decrease with temperature at constant pressure and particle number.

I'm pretty sure I've got it right this time, but feel free to correct me. Does this make sense?

-James
 
  • #19
No, you need the expansion of U or P in T in the same order as that for mu.
Come on, it's not so difficult...
 
  • #20
[itex] E_F \approx C_1 P^{2/3}(1-\frac{\pi^2 k^2T^2}{9 (C_1P^{2/3})^2})[/itex]
Sufficient to render the dependence of mu on T negative.
 
Last edited:
  • #21
DrDu said:
[itex] E_F \approx C_1 P^{2/3}(1-\frac{\pi^2 k^2T^2}{9 (C_1P^{2/3})^2})[/itex]
Sufficient to render the dependence of mu on T negative.


could you show how you derived this?

-James
 
  • #22
jfizzix said:
could you show how you derived this?

-James
[itex] Lp=-\Omega=\frac{2}{\beta}\int_0^\infty d\omega D(\omega) \ln (1 +\exp(-\beta(\omega-\mu)) [/itex]
Here, [itex] \Omega [/itex] is the grand potential, D is the density of states [itex] L\frac{\sqrt{2m}}{\pi \hbar}\omega^{-1/2}=AL \omega^{-1/2}[/itex] with [itex]A=\alpha_1^{-1/2}/2[/itex].
Partial integration:
[itex]
p=-\frac{2A}{\beta} \int_0^\infty d\omega \left(\int_0^\omega d\omega' \omega'^{-1/2} \right) \frac{-\beta}{\exp(\beta(\omega-\mu))+1} [/itex]
Sommerfeld expansion...
[itex]p\approx A\left[ \frac{3}{2}E_F^{3/2}+\frac{\pi^2}{6\beta^2}E_F^{-1/2}\right] [/itex]
Solve for E_F at same level of accuracy.
 
Last edited:

What is chemical potential?

Chemical potential is a measure of the potential energy of a substance's molecules to undergo a chemical reaction. It is related to the concentration of a substance and its temperature and pressure.

How does chemical potential vary with pressure?

As pressure increases, the molecules of a substance are forced closer together, leading to a higher concentration and therefore a higher chemical potential. This is why gases tend to become more reactive at higher pressures.

How does chemical potential vary with temperature?

Temperature affects chemical potential by changing the average kinetic energy of molecules. As temperature increases, molecules move faster and collide more frequently, leading to a higher chemical potential.

What are Fermi gases?

Fermi gases are systems of particles that obey Fermi-Dirac statistics, meaning they follow the Pauli exclusion principle and cannot occupy the same quantum state. Examples include electrons in a metal or neutrons in a neutron star.

Why is it difficult to study Fermi gases?

Fermi gases are difficult to study because they can exhibit complex behavior due to their interactions and quantum nature. Additionally, experimental conditions must be carefully controlled to observe their properties accurately.

Similar threads

Replies
6
Views
858
  • Introductory Physics Homework Help
Replies
4
Views
750
Replies
1
Views
718
Replies
3
Views
958
  • Advanced Physics Homework Help
Replies
1
Views
631
Replies
6
Views
3K
  • Introductory Physics Homework Help
Replies
3
Views
1K
  • Materials and Chemical Engineering
Replies
4
Views
2K
  • Introductory Physics Homework Help
Replies
1
Views
997
Replies
4
Views
1K
Back
Top