Articulo de referencia

cribado de Thomas-Fermi

El apantallamiento de Thomas-Fermi es un enfoque teórico para calcular los efectos del apantallamiento del campo eléctrico por electrones en un sólido. [ 1 ] Es un caso especial...

El apantallamiento de Thomas-Fermi es un enfoque teórico para calcular los efectos del apantallamiento del campo eléctrico por electrones en un sólido. [ 1 ] Es un caso especial de la teoría de Lindhard más general ; en particular, el apantallamiento de Thomas-Fermi es el límite de la fórmula de Lindhard cuando el vector de onda (el recíproco de la escala de longitud de interés) es mucho menor que el vector de onda de Fermi, es decir, el límite de larga distancia. [ 1 ] Recibe su nombre de Llewellyn Thomas y Enrico Fermi .

El vector de onda de Thomas-Fermi (en unidades Gaussian-cgs ) es [ 1 ].k02=4πmi2norteμ,{\displaystyle k_{0}^{2}=4\pi e^{2}{\frac {\partial n}{\partial \mu }},} donde μ es el potencial químico ( nivel de Fermi ), n es la concentración de electrones y e es la carga elemental .

Por ejemplo, en el caso de semiconductores que no están demasiado dopados, la densidad de carga ne μ / k B T , donde k B es la constante de Boltzmann y T es la temperatura. En este caso, k02=4πmi2nortekBT,{\displaystyle k_{0}^{2}={\frac {4\pi e^{2}n}{k_{\rm {B}}T}},}

es decir, 1 / k0 viene dado por la fórmula familiar para la longitud de Debye . En el extremo opuesto, en el límite de baja temperatura T = 0 , los electrones se comportan como partículas cuánticas ( fermiones ). Tal aproximación es válida para metales a temperatura ambiente, y el vector de onda de apantallamiento de Thomas-Fermi kTF dado en unidades atómicas es kTF2=4(3norteπ)1/3.{\displaystyle k_{\rm {TF}}^{2}=4\left({\frac {3n}{\pi }}\right)^{1/3}.}

Si restauramos la masa del electrónmetromi{\displaystyle m_{e}}y la constante de Planck{\displaystyle \hbar }, el vector de onda de apantallamiento en unidades gaussianas esk02=kTF2(metromi/2){\displaystyle k_{0}^{2}=k_{\rm {TF}}^{2}(m_{e}/\hbar ^{2})}.

Para obtener más detalles y un análisis más profundo, incluyendo los casos unidimensionales y bidimensionales, consulte el artículo sobre la teoría de Lindhard .

Derivación

Relación entre la densidad electrónica y el potencial químico interno

El potencial químico interno (estrechamente relacionado con el nivel de Fermi , véase más adelante) de un sistema de electrones describe la energía necesaria para añadir un electrón al sistema, sin tener en cuenta la energía potencial eléctrica. A medida que aumenta el número de electrones en el sistema (con temperatura y volumen constantes), el potencial químico interno también aumenta. Esto se debe principalmente a que los electrones cumplen el principio de exclusión de Pauli : solo un electrón puede ocupar un nivel de energía, y los niveles de energía más bajos ya están ocupados, por lo que los nuevos electrones deben ocupar niveles de energía cada vez más altos.

Dado un gas de Fermi de densidadnorte{\displaystyle n}, el estado de momento ocupado más alto (a temperatura cero) se conoce como momento de Fermi,kF{\displaystyle k_{\rm {F}}}.

Entonces, la relación requerida se describe mediante la densidad numérica de electrones.norte(μ){\displaystyle n(\mu )}en función de μ , el potencial químico interno. La forma funcional exacta depende del sistema. Por ejemplo, para un gas de Fermi tridimensional , un gas de electrones no interactuantes, a temperatura de cero absoluto, la relación esnorte(μ)μ3/2{\displaystyle n(\mu )\propto \mu ^{3/2}}.

Prueba: Incluyendo la degeneración de espín,

norte=21(2π)343πkF3,μ=2kF22metro.{\displaystyle n=2{\frac {1}{(2\pi )^{3}}}{\frac {4}{3}}\pi k_{\rm {F}}^{3}\quad ,\quad \mu ={\frac {\hbar ^{2}k_{\rm {F}}^{2}}{2m}}.}

(en este contexto, es decir, el cero absoluto, el potencial químico interno se denomina más comúnmente energía de Fermi ).

Como otro ejemplo, para un semiconductor de tipo n con una concentración de electrones baja a moderada,norte(μ)miμ/kBT{\displaystyle n(\mu )\propto e^{\mu /k_{\rm {B}}T}}.

aproximación local

La principal suposición del modelo de Thomas-Fermi es que existe un potencial químico interno en cada punto r que depende únicamente de la concentración de electrones en ese mismo punto . Este comportamiento no puede ser del todo cierto debido al principio de incertidumbre de Heisenberg . Ningún electrón puede existir en un solo punto; cada uno se distribuye en un paquete de ondas de tamaño ≈ 1 / k F , donde k F es el número de onda de Fermi , es decir, un número de onda típico para los estados en la superficie de Fermi . Por lo tanto, no es posible definir un potencial químico en un solo punto, independientemente de la densidad de electrones en puntos cercanos.

Sin embargo, es probable que el modelo de Thomas-Fermi sea una aproximación razonablemente precisa siempre que el potencial no varíe mucho en longitudes comparables o menores que 1 / k F. Esta longitud suele corresponder a unos pocos átomos en los metales.

Electrones en equilibrio, ecuación no lineal

Finalmente, el modelo de Thomas-Fermi supone que los electrones están en equilibrio, lo que significa que el potencial químico total es el mismo en todos los puntos. (En terminología electroquímica, "el potencial electroquímico de los electrones es el mismo en todos los puntos". En terminología de física de semiconductores, "el nivel de Fermi es plano"). Este equilibrio requiere que las variaciones en el potencial químico interno se correspondan con variaciones iguales y opuestas en la energía potencial eléctrica . Esto da lugar a la "ecuación básica de la teoría no lineal de Thomas-Fermi": [ 1 ]ρinducido(r)=mi[norte(μ0+miϕ(r))norte(μ0)]{\displaystyle \rho ^{\text{induced}}(\mathbf {r} )=-e[n(\mu _{0}+e\phi (\mathbf {r} ))-n(\mu _{0})]} donde n ( μ ) es la función discutida anteriormente (densidad electrónica en función del potencial químico interno), e es la carga elemental , r es la posición yρinducido(r){\displaystyle \rho ^{\text{induced}}(\mathbf {r} )}es la carga inducida en r . El potencial eléctricoϕ{\displaystyle \phi }se define de tal manera queϕ(r)=0{\displaystyle \phi (\mathbf {r} )=0}en los puntos donde el material es eléctricamente neutro (el número de electrones es exactamente igual al número de iones), y de manera similar μ 0 se define como el potencial químico interno en los puntos donde el material es eléctricamente neutro.

Linealización, función dieléctrica

Si el potencial químico no varía demasiado, la ecuación anterior se puede linealizar: ρinducido(r)mi2norteμϕ(r){\displaystyle \rho ^{\text{induced}}(\mathbf {r} )\approx -e^{2}{\frac {\partial n}{\partial \mu }}\phi (\mathbf {r} )} dóndenorte/μ{\displaystyle \partial n/\partial \mu }se evalúa en μ 0 y se trata como una constante.

Esta relación se puede convertir en una función dieléctrica dependiente del vector de onda : [ 1 ] (en unidades cgs-gaussianas )ε(q)=1+k02q2{\displaystyle \varepsilon (\mathbf {q} )=1+{\frac {k_{0}^{2}}{q^{2}}}} dónde k0=4πmi2norteμ.{\displaystyle k_{0}={\sqrt {4\pi e^{2}{\frac {\partial n}{\partial \mu }}}}.} A grandes distancias ( q → 0 ), la constante dieléctrica tiende al infinito, lo que refleja el hecho de que las cargas se aproximan cada vez más a un apantallamiento perfecto a medida que se observan desde más lejos.

Ejemplo: Una carga puntual

Si se coloca una carga puntual Q en r = 0 en un sólido, ¿qué campo producirá, teniendo en cuenta el apantallamiento electrónico?

Se busca una solución coherente a dos ecuaciones:

  • La fórmula de apantallamiento de Thomas-Fermi proporciona la densidad de carga en cada punto r en función del potencial.ϕ(r){\displaystyle \phi (\mathbf {r} )}en ese momento.
  • La ecuación de Poisson (derivada de la ley de Gauss ) relaciona la segunda derivada del potencial con la densidad de carga.

Para la fórmula no lineal de Thomas-Fermi, resolver estas ecuaciones simultáneamente puede ser difícil, y generalmente no existe una solución analítica. Sin embargo, la fórmula linealizada tiene una solución simple (en unidades cgs-gaussianas ):ϕ(r)=Qrmik0r{\displaystyle \phi (\mathbf {r} )={\frac {Q}{r}}e^{-k_{0}r}} Con k 0 = 0 (sin apantallamiento), esto se convierte en la conocida ley de Coulomb .

Cabe señalar que, además del apantallamiento aquí descrito, puede existir permitividad dieléctrica ; por ejemplo, debido a la polarización de los electrones internos inmóviles. En ese caso, sustituya Q por Q / ε , donde ε es la permitividad relativa debida a estas otras contribuciones.

Gas de Fermi a temperatura arbitraria

Temperatura efectiva para el apantallamiento de Thomas-Fermi. La forma aproximada se explica en el artículo y utiliza la potencia p=1,8.

Para un gas de Fermi tridimensional (gas de electrones no interactuantes), el vector de onda de apantallamientok0{\displaystyle k_{0}}puede expresarse como una función tanto de la temperatura como de la energía de Fermi.miF{\displaystyle E_{\rm {F}}}El primer paso es calcular el potencial químico interno.μ{\displaystyle \mu }, que implica la inversa de una integral de Fermi-Dirac , μkBT=F1/21[23Γ(3/2)(miFT)3/2].{\displaystyle {\frac {\mu }{k_{\rm {B}}T}}=F_{1/2}^{-1}\left[{2 \over {3\Gamma (3/2)}}\left({E_{\rm {F}} \over T}\right)^{3/2}\right].}

Podemos expresark0{\displaystyle k_{0}}en términos de una temperatura efectivaTmiFF{\displaystyle T_{\rm {eff}}}:k02=4πmi2norte/kBTmiFF{\displaystyle k_{0}^{2}=4\pi e^{2}n/k_{\rm {B}}T_{\rm {eff}}}, okBTmiFF=norteμ/norte{\displaystyle k_{\rm {B}}T_{\rm {eff}}=n\partial \mu /\partial n}. El resultado general paraTmiFF{\displaystyle T_{\rm {eff}}}esTmiFFT=43Γ(1/2)(miF/kBT)3/2F1/2(μ/kBT).{\displaystyle {T_{\rm {eff}} \over T}={4 \over 3\Gamma (1/2)}{(E_{\rm {F}}/k_{\rm {B}}T)^{3/2} \over F_{-1/2}(\mu /k_{\rm {B}}T)}.}En el límite clásicokBTmiF{\displaystyle k_{\rm {B}}T\gg E_{\rm {F}}}, encontramosTmiFF=T{\displaystyle T_{\rm {eff}}=T}, mientras que en el límite degeneradokBTmiF{\displaystyle k_{\rm {B}}T\ll E_{\rm {F}}}encontramos kBTmiFF=(2/3)miF.{\displaystyle k_{\rm {B}}T_{\rm {eff}}=(2/3)E_{\rm {F}}.} Una forma aproximada simple que recupera ambos límites correctamente es kBTmiFF=[(kBT)pag+(2miF/3)pag]1/pag,{\displaystyle k_{\rm {B}}T_{\rm {eff}}=\left[(k_{\rm {B}}T)^{p}+(2E_{\rm {F}}/3)^{p}\right]^{1/p},} para cualquier potenciapag{\displaystyle p}. Un valor que ofrece una concordancia decente con el resultado exacto para todoskBT/miF{\displaystyle k_{\rm {B}}T/E_{\rm {F}}}espag=1.8{\displaystyle p=1.8}, [ 2 ] que tiene un error relativo máximo de < 2,3%.

En la temperatura efectiva dada anteriormente, la temperatura se utiliza para construir un modelo clásico efectivo. Sin embargo, esta forma de la temperatura efectiva no recupera correctamente el calor específico y la mayoría de las demás propiedades del sistema finito.T{\displaystyle T}fluido electrónico incluso para el gas de electrones no interactuantes. Por supuesto, no intenta incluir los efectos de interacción electrón-electrón. Una forma simple para una temperatura efectiva que recupera correctamente todas las propiedades funcionales de la densidad incluso del gas de electrones interactuantes , incluidas las funciones de distribución de pares en un valor finito.T{\displaystyle T}, se ha dado utilizando el modelo clásico de cadena hiperredada ( CHNC ) del fluido de electrones. Es decir, TmiFFmiF=(T2miF2+Tq2miF2)1/2{\displaystyle {\frac {T_{\rm {eff}}}{E_{\rm {F}}}}=\left({\frac {T^{2}}{E_{\rm {F}}^{2}}}+{\frac {T_{q}^{2}}{E_{\rm {F}}^{2}}}\right)^{1/2}} donde la temperatura cuánticaTq{\displaystyle T_{q}}se define como: TqmiF=1a+brs+dors{\displaystyle {\frac {T_{q}}{E_{\rm {F}}}}={\frac {1}{a+b{\sqrt {r_{\rm {s}}}}+cr_{\rm {s}}}}} donde a = 1,594 , b = −0,3160 , c = 0,0240 . Aquírs{\displaystyle r_{\rm {s}}}es el radio de Wigner-Seitz correspondiente a una esfera en unidades atómicas que contiene un electrón. Es decir, sinorte{\displaystyle n}es el número de electrones en un volumen unitario utilizando unidades atómicas donde la unidad de longitud es el Bohr, es decir,5,291 77 × 10 −9  cm , entonces rs=(34πnorte)1/3.{\displaystyle r_{\rm {s}}=\left({\frac {3}{4\pi n}}\right)^{1/3}.} Para un gas de electrones denso, por ejemplo, conrs1{\displaystyle r_{\rm {s}}\approx 1}o menos, las interacciones electrón-electrón se vuelven insignificantes en comparación con la energía de Fermi, entonces, usando un valor ders{\displaystyle r_{\rm {s}}}cerca de la unidad, vemos que la temperatura efectiva CHNC enT=0{\displaystyle T=0}se aproxima a la forma2miF/3{\displaystyle 2E_{\rm {F}}/3}. También se han dado otros mapeos para el caso 3D, [ 3 ] y fórmulas similares para la temperatura efectiva para el mapeo clásico del gas de electrones bidimensional. [ 4 ]

Véase también

Referencias

  1. 1 2 3 4 5 N. W. Ashcroft y ND Mermin, Física del estado sólido (Thomson Learning, Toronto, 1976)
  2. Stanton, Liam G.; Murillo, Michael S. (2016-04-08). "Transporte iónico en materia de alta densidad energética" . Physical Review E. 93 ( 4) 043203. American Physical Society (APS). Bibcode : 2016PhRvE..93d3203S . doi : 10.1103/physreve.93.043203 . ISSN 2470-0045 . PMID 27176414 .  
  3. ^ Yu Liu y Jianzhong Wu, J. Chem. Física. 141 064115 (2014)
  4. ^ François Perrot y MWC Dharma-wardana, Phys. Rev. Lett. 87 , 206404 (2001)