4 votos

Potencial químico para el gas fermi ideal con energías de partículas individuales$\epsilon_k^{\pm}=\pm \hbar c |\mathbf k|$

Me encontré con la siguiente pregunta en algunos problemas de ejemplo para mi de posgrado de la mecánica estadística curso examen final:

Considere un gas ideal de fermiones de la densidad de $n$ en tres dimensiones con la sola partícula eigenstate energías dadas por $\epsilon_k^{\pm}=\pm \hbar c |\mathbf k|$. Suponga que el potencial químico $\mu = 0$$T=0$. Demostrar que si el fermión de la densidad es constante, $\mu(T)=0$ todos los $T$.

Mi primer acercamiento fue el uso de la de Fermi-Dirac distribución (el uso de unidades de energía para la temperatura): $$N=\sum_k\langle n_k\rangle = \sum_k \frac{1}{e^{\frac {\epsilon_k-\mu}T}+1}=\sum_k \frac{1}{e^{\frac {\hbar c |\mathbf k|-\mu}T}+1}$$ La conversión de este en un integrante en el impulso de espacio: $$N=V\int_{\Bbb R^3}\frac {d^3k}{(2 \pi)^3}\frac{1}{e^{\frac {\hbar c |\mathbf k|-\mu}T}+1}=\frac {V}{2 \pi^2}\int_0^\infty dk \ \frac{k^2}{z^{-1}e^{\frac {\hbar c k}T}+1}$$ Donde $z=e^{\frac {\mu}{T}}$ es la fugacidad. Mediante el cambio de variables, obtenemos: $$N=\frac {V}{2 \pi^2} \left(\frac{T}{\hbar c}\right)^3\int_0^\infty dx \ \frac{x^2}{z^{-1}e^{x}+1}$$ Que en términos de la de Fermi-Dirac funciones, es : $$n=\frac {1}{ \pi^2} \left(\frac{T}{\hbar c}\right)^3f_3(z)$$ Donde $n$ es la densidad. A bajas temperaturas, el uso de la Sommerfeld expansión $f_3\left(e^{\frac{\mu}{T}}\right) \simeq \frac 16\left(\frac {\mu}{T}\right)^3\left(1+\pi^2 \left(\frac T{\mu}\right)^2\right)$, obtenemos: $$n\simeq\frac {1}{ \pi^2} \left(\frac{T}{\hbar c}\right)^3\frac 16\left(\frac {\mu}{T}\right)^3\left(1+\pi^2 \left(\frac T{\mu}\right)^2\right)\simeq \frac {V}{ 6\pi^2} \left(\frac{\mu}{\hbar c}\right)^3 \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ T\rightarrow0^+ $$ Lo que da : $$\mu \simeq \hbar c\left(6 \pi^2 n\right)^{\frac{1}{3}} \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ T\rightarrow0^+ $$ Yo no puedo ver cómo podemos asumir que el potencial químico a ser cero en el cero de temperatura. Claramente, la solución anterior muestra que eso es imposible. ¿Cómo podría el gas cero de la energía de fermi? Puedo obtener el mismo resultado a la hora de calcular el grand canónica de la función de partición desde cero y utilizando el gran potencial. Un cero de potencial químico no hace sentido intuitivo, porque implica que la adición o sustracción de las partículas del gas no tendría coste de energía de cualquier tipo; que no puede ser verdad, porque cuando se añade un fermión no puede compartir un nivel de energía con otra partícula (principio de Pauli), y por lo tanto tiene que ir a un nivel de energía más alto, lo que significa que hay un valor distinto de cero el costo de la energía.

Otro problema es que yo no puede entender lo que la pregunta se entiende por la $\pm$ $\epsilon_k^{\pm}=\pm \hbar c |\mathbf k|$ relación. Qué significa que cada eigenstate de impulso corresponde a dos autoestados de energía, uno negativo y uno positivo? Por ejemplo, considere el caso de $\epsilon_k=- \hbar c |\mathbf k|$. El mismo razonamiento anterior nos da: $$n=\frac {1}{2 \pi^2} \left(\frac{T}{\hbar c}\right)^3\int_0^\infty dx \ \frac{x^2}{z^{-1}e^{-x}+1} \ \rightarrow +\infty$$ Y esta integral diverge claramente para todos los valores de $z$.

Me siento como que estoy interpretando mal la pregunta. Cualquier ayuda se agradece.

5voto

Bemmu Puntos 108

Antes de profundizar en las matemáticas, vamos a tratar de conseguir algo de intuición acerca de la función de la energía que se está tratando. Se le da

$$\epsilon_{\boldsymbol{k}}^{\pm}=\pm\hbar c\left|\boldsymbol{k}\right|$$

$$\epsilon_{\rm F}=\mu\left(T=0\right)=0$$

Esto significa que en $T=0$ de toda la energía negativa de los estados están ocupados y toda la energía positiva que los estados están vacías. Esto es típico de los sistemas en la física de estado sólido. Básicamente tenemos aquí una Estructura de Banda: las energías negativas forman lo que se conoce como la Cenefa de la Banda, y las energías positivas de lo que se conoce como la Banda de Conducción. La cenefa de la banda está totalmente ocupada con los electrones, lo que la hace difícil de sobrellevar. Por lo tanto, los físicos crean el concepto de Agujeros, o falta de electrones. Mientras que la cenefa de la banda está llena de electrones, es equivalente a decir que está vacía de agujeros - fácil de tratar. En resumen, los agujeros de vivir en la cenefa de la banda y de los electrones en vivo en la banda de conducción.

Ahora viene la parte clave. Al excitar un electrón de la cenefa de la banda para la banda de conducción que son esencialmente de la creación de un par de electrón-hueco: un agujero (de electrones de la vacante) aparece en la cenefa de la banda y un electrón aparece en la banda de conducción. Es muy importante señalar que, debido a esta simetría, el número de electrones es igual al número de agujeros. Esto es cierto sólo porque el nivel de Fermi se encuentra exactamente en el medio entre las dos bandas.

Ahora entendemos la intuición detrás de este problema. El siguiente paso es calcular el número de electrones y el número de agujeros, exigir a los dos a ser el mismo y obtenga $\mu=\mu\left(T\right)$. Como usted dijo, la probabilidad de ocupación de electrones en una energía $\epsilon$ está dada por la de Fermi-Dirac Distribución de la función

$$f_{\rm e}\left(\epsilon\right)=\frac{1}{e^{\beta\left(\epsilon-\mu\right)}+1}$$

con $\beta\equiv \frac{1}{K_{\rm B}T}$. Por otro lado, la probabilidad de ocupación de un agujero en la cenefa de la banda es igual a la probabilidad de que la falta de electrones

$$f_{\rm h}\left(\epsilon\right)=1-f_{\rm e}\left(\epsilon\right)=1-\frac{1}{e^{\beta\left(\epsilon-\mu\right)}+1}=\frac{1}{e^{-\beta\left(\epsilon-\mu\right)}+1}$$

Por lo tanto el número de electrones y los huecos, respectivamente, está dada por

$$n_{\rm e}\left(T\right)=\int_{0}^{\infty}{\rm d}\epsilon D\left(\epsilon\right)f_{\rm e}\left(\epsilon\right)$$

$$n_{\rm h}\left(T\right)=\int_{-\infty}^{0}{\rm d}\epsilon D\left(\epsilon\right)f_{\rm h}\left(\epsilon\right)=\int_{0}^{\infty}{\rm d}\epsilon D\left(-\epsilon\right)f_{\rm h}\left(-\epsilon\right)$$

en todas las temperaturas. Aquí $D\left(\epsilon\right)$ es la densidad de estados, y se puede argumentar que a partir de la simetría $D\left(\epsilon\right)=D\left(-\epsilon\right)$. El taponamiento de las distribuciones, obtenemos

$$n_{\rm e}\left(T\right)=\int_{0}^{\infty}{\rm d}\epsilon D\left(\epsilon\right)\frac{1}{e^{\beta\left(\epsilon-\mu\right)}+1}$$

$$n_{\rm h}\left(T\right)=\int_{0}^{\infty}{\rm d}\epsilon D\left(\epsilon\right)\frac{1}{e^{\beta\left(\epsilon+\mu\right)}+1}$$

Ahora usted puede ver inmediatamente que si $n_{\rm e}\left(T\right)=n_{\rm h}\left(T\right)$ es cierto para todas las temperaturas, entonces usted debe tener

$$\mu\left(T\right)=0$$

para todas las temperaturas, como quería.

2voto

Kite.Y Puntos 141

@eranreches la respuesta es completa y suficiente. Sólo quiero agregar algunos pequeños comentarios acerca de la solución en tu pregunta original:

En primer lugar, como ya se ha visto, hay dos ramas, mientras que en su cálculo se utiliza sólo uno de ellos: eso significa que usted está tratando de llenar los electrones en único positivo de energía de los estados, entonces cualquier número finito de llenado resultaría en un positivo potencial químico $\mu$$T=0$, que es el nivel más alto de energía alcanzado después de rellenar todos los electrones.

En segundo lugar, si usted también toma la energía negativa en cuenta, pero aún así tratar de utilizar un estándar de Fermi-Dirac integral para el cálculo, entonces te darías cuenta de que hay un número infinito de energía negativa de los niveles puede ser llenado -- no hay límite inferior para los niveles de energía, lo que resulta en un desastre ya que no sabes ni por donde empezar el llenado de los electrones. Y eso es parte de la razón @eranreches utiliza la densidad de estados de $D(\epsilon)$ en su cálculo, lo que podría evitar el problema en un cálculo explícito. Mientras que en el mundo real, por ejemplo, en estado sólido, sistemas, cualquier banda de energía de curso está restringido en un número finito de ventana de energía, y también el impulso $k$ tendría un valor finito; luego, más específicamente, usted puede tener $D(|\epsilon| > \epsilon_0) = 0$.

En tercer lugar, la dispersión lineal es bastante común en estado sólido, se relaciona con muchos temas interesantes: por ejemplo, Luttinger líquido, conos de Dirac.... Por lo general, una aproximación para sólo un número finito de estados diseñados por $k$ en el Brilioun de la zona, y se ve sobre todo para la banda de cruce (o tocar) caso: cuando dos bandas de música de la cruz el uno al otro, la región alrededor del punto de cruce tendría una dispersión lineal $\epsilon \sim |k|$. Y el significado de $\mu(T=0)=0$ en este caso es simplemente que los electrones llenan todos los estados por debajo de este punto de cruce.

En general no serían otras que las dispersiones para la banda tocando, por ejemplo,$\epsilon\sim k^2$. (Usted podría encontrar algo relacionado aquí: arXiv 1603.03093.)

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X