13.3: Derivación del Lagrangiano
( \newcommand{\kernel}{\mathrm{null}\,}\)
El propósito de este capítulo es encontrar el voltajeV(r) y la densidad de cargaρch(r) alrededor de un átomo, y utilizaremos cálculo de variaciones para lograr esta tarea. Necesitamos hacer algunas suposiciones bastante severas para que este problema sea manejable. Considera un átomo neutro aislado con muchos electrones a su alrededor. SupongamosT≈0 K, por lo que todos los electrones ocupan los niveles de energía más bajos posibles. Supongamos que el átomo es esféricamente simétrico. Todas las cantidades que encontramos, como voltaje, densidad de carga y Lagrangian, varían conr pero no varían conθ oϕ. Utilizaremos coordenadas esféricas con el origen en el núcleo del átomo. Si bien las cantidades varían con la posición, supongamos que ninguna cantidad varía con el tiempo. La densidad de carga nosρch(r) indica dónde es más probable que se encuentren los electrones en promedio. Se relaciona con la función de onda mecánica cuántica,ψ, por
ρch=−q⋅|ψ|2
dondeq está la magnitud de la carga de un electrón. Supongamos que todos los electrones que rodean al átomo están distribuidos uniformemente y pueden tratarse como si fueran una nube electrónica uniforme de alguna densidad de carga.
Escoge uno de los electrones del átomo, y considera lo que sucede cuando el electrón se mueve radialmente hacia dentro y hacia afuera. La figura13.3.1 ilustra esta situación. A medida que el electrón se mueve, ocurre la conversión de energía El objetivo de esta sección es anotar lo hamiltoniano y lagrangiano para este proceso de conversión energética. Escribimos estas cantidades en las unidades de energía por unidad de volumen por electrón de valencia bajo consideración.

Para entender lo que sucede cuando se mueve el electrón, considera la energía del átomo con más detalle. La ley de Coulomb, introducida en la Ecuación 1.6.2, nos dice que los objetos cargados ejercen fuerzas sobre otros objetos cargados. Más específicamente, la intensidad del campo eléctrico→E debida a una carga puntual deQ culombios a unar distancia rodeada por un material con permitividadϵ viene dada por
→E=Qˆar4πϵr2.
El átomo está compuesto por protones cargadosN positivamente. El electrón bajo consideración siente una fuerza de Coulomb atractiva debido a estos protones. Adicionalmente, el átomo tieneN electrones, yN−1 de estos ejercen una fuerza repulsiva de Coulomb sobre el electrón bajo consideración. Dado que existe una separación de carga y un campo eléctrico, se almacena la energía. Llamar al componente de la energía del átomo debido a la interacción Coulomb entre los protones del núcleo y el electrón bajo consideraciónECoulombenucl. Llamar a la interacción de Coulomb entre el electrón bajo consideración y todos los demás electronesEeeinteract. El átomo también tiene energía cinética. Llama a la energía cinética del núcleoEkineticnucl y a la energía cinética de todos los electronesEkinetice. La energía del átomo es la suma de todos estos términos.
Eatom=ECoulombenucl.+Ekineticnucl+Eeeinteract+Ekinetice
La energía debida al giro de los electrones y protones se ignora al igual que la energía debida a la interacción con cualquier otro objeto cargado cercano. AT≈0 K, la energía cinética del núcleo será cercana a cero, por lo que podemos ignorar el término,Ekineticnucl≈0. La cantidadEkinetice no puede ser exactamente cero. En el capítulo 6 trazamos diagramas de nivel de energía para electrones alrededor de un átomo. Incluso enT=0 K, los electrones tienen algo de energía interna, y esta energía se denota por el nivel de energía ocupado.
Si tenemos un átomo grande con muchos electrones a su alrededor, la interacción Coulomb entre cualquier electrón y el núcleo está protegida por la interacción Coulomb de todos los demás electrones. Más específicamente, supongamos que tenemos un átomo aislado conN protones en el núcleo yN electrones a su alrededor. Si escogemos uno de los electrones,ECoulombenucl para ese electrón se describe la energía almacenada en el campo eléctrico debido a la separación de carga entre el núcleo de carga positivaNq y ese electrón. Sin embargo, también hayN−1 otros electrones que tienen una carga negativa. El términoEeeinteract describe la energía almacenada en el campo eléctrico debido a la separación de carga entre losN−1 otros electrones y el electrón bajo consideración. Estos términos se cancelan un poco entre sí porque el electrón bajo consideración interactúa conN protones cada uno de carga positiva q yN−1 electrones cada uno de carga negativa−q. No obstante, los términos no desaparecen por completo. Calculando
ECoulombenucl+Eeeinteract
es complicado porque los electrones están en movimiento, y no sabemos realmente dónde están o incluso dónde es más probable que se encuentren. De hecho, estamos tratando de resolver por dónde es probable que se encuentren.
A medida que movemos el electrón bajo consideración hacia adentro y hacia afuera radialmente, la energía se transfiere entre (ECoulombenucl+Eeeinteract) yEkinetice. El hamiltoniano es la suma de estas dos formas de energía por unidad de volumen, y el lagrangiano es la diferencia de estas dos formas de energía por unidad de volumen. Ambas cantidades tienen las unidadesJm3. Elija voltajeV(r) como ruta generalizada y densidad de cargaρch(r) como potencial generalizado. La variable independiente de estas cantidades es la posición radialr, no el tiempo. Ya podemos escribir el hamiltoniano y lagrangiano.
H(r,V,dVdr)=(ECoulombenuclV+EeeinteractV)+EkineticeV
L(r,V,dVdr)=(ECoulombenuclV+EeeinteractV)−EkineticeV
El siguiente paso es escribir
ECoulombenuclV+EeeinteractV
en cuanto al caminoV. Como se detalla en la Tabla 12.2.3, la densidad de energía debida a un campo eléctrico→E viene dada por
EV=12ϵ|→E|2.
Recuerda queE representa la energía mientras que→E representa el campo eléctrico. El campo eléctrico es el gradiente negativo de la tensiónV(r).
→E=−→∇V.
Podemos combinar estas expresiones y la Ecuación 13.2.6 para escribir el primer término del hamiltoniano y el lagrangiano en términos del camino generalizado.
ECoulombenuclV+EeeinteractV=12ϵ|→∇V|2
H(r,V,dVdr)=(12ϵ|→∇V|2)+EkineticeV
L(r,V,dVdr)=(12ϵ|→∇V|2)−EkineticeV
La siguiente tarea es describir el término restanteEkineticeV como una función del camino generalizado también. Esta tarea es un poco más desafiante. Seguimos tomando el enfoque de hacer aproximaciones severas hasta que sea manejable. Tenemos queρch(r) expresarnos en función deV(r). Entonces con algo de álgebra, seEkineticeV puede escribir puramente en función deV(r).
Queremos generalizar sobre la energía cinética de los electrones. Sin embargo, cada electrón tiene su propia velocidad→v e impulso→M. Estas cantidades dependen de la posición
→r=rˆar+θˆaθ+ϕˆaϕ
de alguna manera desconocida. Además, el cálculo deEkineticeV depende de la densidad de cargaρch(r), que es la cantidad desconocida que estamos tratando de encontrar. Tenemos más suerte al describir estas cantidades en el espacio recíproco, introducido en la Sec. 6.3. La posición se denota en el espacio recíproco por un vector de onda
→k=˜rˆar+˜θˆaθ+˜ϕˆaϕ
Podemos describir las propiedades de un material describiendo cómo varían con la posición en el espacio real. Por ejemplo,ρch(r) representa la densidad de carga de los electrones en función de la distancia r desde el centro del átomo. Nos puede interesar cómo otras cantidades, como la energía requerida para arrancar un electrón o la energía cinética interna a un electrón, varían también con la posición en el espacio real. En lugar de describir cómo varían las cantidades con la posición en el espacio real, podemos describir cómo varían las cantidades con la frecuencia espacial de los electrones. Esta es la idea detrás de representar cantidades en el espacio recíproco. Nos puede interesar cómo varía la densidad de carga de los electrones en función de la frecuencia espacial de las cargas en un cristal u otro material, y esta es la idea representada por funciones de vector de onda comoρch(→k). Estamos tratando de resolver por densidad de cargaρch(r). Esperamos que los electrones sean más propensos a encontrarse a ciertas distanciasr del centro del átomo que a otras distancias. Sin embargo, no hay patrón a la densidad de carga en función del vector de onda,ρch(→k). Supongamos queρch es aproximadamente constante con respecto a|→k| hasta algún nivel. Con algo más de trabajo, esta suposición nos permitirá resolver por densidad de cargaρch(r).
La energía cinética de un solo electrón viene dada por
Ekineticee−=12m|→v|2
dondem está la masa del electrón. Podemos escribir esta energía en términos de ímpetu,→M=m→v. (Tenga en cuenta que el impulso→M y el impulso generalizadoM son diferentes y tienen diferentes unidades.)
Ekineticee−=|→M|22m
No sabemos cómo varía la energía en función de la posición r. En cambio, podemos escribir la energía en función del momento cristalino→Mcrystal o del vector de onda→k, y sabemos algo sobre la variación de estas cantidades. El impulso del cristal es igual al vector de onda escalado por la constante de Planck.
→Mcrystal=ℏ→k
Tiene las unidades de impulsokg⋅ms, y se introdujo en la Sec. 6.3.2. La energía cinética de un electrón en función del momento cristalino viene dada por
Ekineticee−=(→Mcrystal)22m=(ℏ|→k|)22m.
Un vector en el espacio recíproco se representa Ecuación\ ref {13.3.14}, y la Ecuación\ ref {13.3.18} puede simplificarse porque estamos asumiendo simetría esférica˜θ=˜ϕ=0. La magnitud del vector de onda se convierte|→k|=˜r, y podemos escribir la energía como
Ekineticee−=ℏ2˜r22m.
Así como cada electrón tiene su propio impulsom|→v|, cada electrón tiene su propio impulso cristalinoℏ|→k|. Sin embargo, conocemos alguna información sobre el vector|→k| de onda de los electrones en el átomo. AT=0 K, los electrones ocupan los estados de energía más bajos permitidos. Los estados energéticos están ocupados hasta algún estado ocupado más alto llamado la energía FermiEf. Mientras que los ingenieros eléctricos usan el término energía Fermi, los químicos a veces usan el término potencial químicoμchem. Los estados energéticos más bajos, están ocupados mientras que los superiores están vacíos. De manera similar, los vectores de onda están ocupados hasta algún vector de onda ocupado más alto llamado vector de onda Fermikf.
|→k|={filled state˜r<kfempty state˜r>kf
La energía de Fermi y el vector de onda Fermi están relacionados por
Ef=ℏ2k2f2m.
Utilizamos la idea de espacio recíproco para escribir una expresión de la energía cinética de los electrones por unidad de volumen [136, p. 49]. La energía cinética debida a cualquier electrón en función de la posición en el espacio recíproco viene dada por la Ecuación\ ref {13.3.19}. Tenga en cuenta que a cada valor de|→k|=˜r, el electrón tiene una energía cinética diferente. Para encontrar la energía cinética por unidad de volumen debida a todos los electrones, nos integramos sobre todo|→k|=˜r en coordenadas esféricas que están ocupadas por electrones, y luego dividimos por el volumen ocupado en el→k espacio.
EkineticeV=1vol. occupied in k space⋅∫filled k levels(Ekineticee−)(e−volume)d(vol. all k space)
El número de electrones por unidad de volumen viene dado por
(e−volume)=−ρchq.
El volumen ocupado en el espacio recíproco es43πk3f, el volumen de una esfera de radiokf.
EkineticeV=143πk3f⋅∫filled k levels(ℏ2˜r22m)(−ρchq)d(vol. all k space)
Un elemento diferencial del volumen se expresa como
d3|→k|=˜r2sin˜θd˜rd˜θd˜ϕ.
EkineticeV=143πk3f⋅∫filled k levels(ℏ2|→k|22m)(−ρchq)(˜r2sin˜θd˜rd˜θd˜ϕ)
Como se describió anteriormente, los electrones ocupan estados en el espacio recíproco solo con0≤˜r≤kf.
EkineticeV=143πk3f⋅kf∫˜r=0π∫˜θ=02π∫˜ϕ=0(ℏ2˜r22m)(−ρchq)(˜r2sin˜θd˜rd˜θd˜ϕ)
La integral anterior puede ser evaluada directamente. El primer paso para evaluarlo es sacar constantes al exterior. Como se describió anteriormente,ρch varía conr pero no˜r, por lo que también se puede tirar fuera de la integral.
EkineticeV=−143πk3f⋅ℏ2ρch2mqkf∫˜r=0π∫˜θ=02π∫˜ϕ=0˜r4sin˜θd˜rd˜θd˜ϕ
La integral se separa y puede ser evaluada.
EkineticeV=−143πk3f⋅ℏ2ρch2mq(π∫˜θ=02π∫˜ϕ=0sin˜θd˜rd˜θd˜ϕ)(kf∫˜r=0˜r4d˜r)
EkineticeV=−143πk3f⋅ℏ2ρch2mq4π(kf55)
EkineticeV=−3ρchk2fℏ210mq
La densidad de carga es una función de la posición en el espacio realr, y estamos en el proceso de resolver para esta función,ρch(r). Sin embargo, también depende de la energía de FermiEf, y de ahí del vector de onda Fermikf, para el átomo. A continuación, encontramos la relación entreρch ykf. Se permiten dos electrones por nivel de energía (spin up y spin down), de ahí pork estado lleno. El número de estados llenos por átomo en el espacio recíproco está relacionado con la densidad de carga.
ρch=−2q(no. filled k statesunit vol. in k space)
En la Sec. 6.3.1, vimos que una celda primitiva en el espacio recíproco era(2π)3 veces la celda primitiva en el espacio real, entonces
(unit vol. k space)=(2π)3⋅(unit vol. real space)=(2π)3.
Sabemos algo sobre los vectores de onda de estados rellenos en el espacio recíproco. AT=0 K, los estados más bajos se llenan, y todos los demás están vacíos, y se llenan hasta un radio dekf. El volumen de una esfera de radiokf viene dado por43πk3f, y esto representa el número dek estados llenos por volumen de espacio recíproco. Por lo tanto, podemos simplificar la expresión anterior.
ρch=−2q⋅43πk3f⋅1(2π)3
ρch=−q3π2k3f
kf=(−3π2qρch)1/3
Queremos escribirEkineticeV en función del camino generalizadoV. Ahora podemos lograr esta tarea combinando Ecuaciones\ ref {13.3.31} y\ ref {13.3.36}.
EkineticeV=−3ℏ210mqρch(−3π2qρch)2/3
EkineticeV=−3ℏ210mq(−3π2q)2/3ρ5/3ch
La energía eléctrica es producto de la carga y el voltaje. Más específicamente, a partir de la Ecuación 2.2.7, viene dada por
E=12QV.
La densidad de energía eléctrica viene dada por
EV=12ρchV.
Utilice la ecuación\ ref {13.3.40} para relacionarρch yV.
EkineticeV=12ρchV=−3ℏ210mq(−3π2q)2/3ρ5/3ch
Ahora hemos relacionado el camino generalizado y el potencial generalizado.
V=−3ℏ25mq(−3π2q)2/3ρ2/3ch
ρch=(−5mq3ℏ2⋅(−3π2q)−2/3)3/2V3/2
ρch=[(−5mq3ℏ2)3/2(−q3π2)]⋅V3/2
Por último, podemos escribirEkineticeV en función deV.
EkineticeV=[(−5mq3ℏ2)3/2(−q3π2)]V5/2
Observe que la cantidad entre paréntesis arriba es constante. El coeficientec0 se define a partir del término entre paréntesis.
c0=(−5mq3ℏ2)3/2(−q3π2)
EkineticeV=c0V5/2
Ahora podemos describir todos los términos del lagrangiano en términos de nuestro camino generalizado.
ECoulombenuclV+EeeinteractV=12ϵ|→∇V|2
EkineticeV=c0V5/2
El hamiltoniano representa la densidad total de energía, y el lagrangiano representa la diferencia de densidad energética de estas formas de energía. El hamiltoniano y lagrangiano tienen la formaH=H(r,V,dVdr) yL=L(r,V,dVdr) dónder está la posición en coordenadas esféricas. No hayθ oϕ dependencia deH oL. Todo es esféricamente simétrico.
H=12ϵ|→∇V|2+c0V5/2
L=12ϵ|→∇V|2−c0V5/2
Como un aparte, consideremos unaEf=μchem vez más la energía Fermi. Con algo de álgebra, podemos escribirlo en función del voltaje. Usa Ecuaciones\ ref {13.3.21},\ ref {13.3.36}, y\ ref {13.3.46}.
Ef=ℏ2k2f2m=ℏ22m(−3π2ρchq)2/3
Ef=ℏ22m(−3π2q)2/3[(−5mq3ℏ2⋅(−3π2q)−2/3)3/2V3/2]2/3
Ef=−5q6V
Observe que la energía Fermi es solo una versión escalada del voltajeV con respecto a un nivel del suelo enr=∞. Los ingenieros eléctricos suelen utilizar la palabra voltaje como sinónimo de potencial. Cuando los químicos usan el término potencial químico, se refieren a la misma cantidad simplemente escalada por una constante. Así como el voltaje es una cantidad fundamental de ingeniería eléctrica que representa lo difícil que es mover electrones, el potencial químico es una cantidad fundamental de química que representa lo difícil que es mover electrones alrededor.