Gelee – Wikipedia

Physikalisches Modell fester Metalle als Elektronengase

Gelee, auch bekannt als die einheitliches Elektronengas (UEG) oder homogenes Elektronengas (HEG), ist ein quantenmechanisches Modell wechselwirkender Elektronen in einem Festkörper, bei dem angenommen wird, dass die positiven Ladungen (dh Atomkerne) gleichmäßig im Raum verteilt sind; die Elektronendichte ist auch im Raum eine einheitliche Größe. Dieses Modell ermöglicht es, sich auf die Effekte in Festkörpern zu konzentrieren, die aufgrund der Quantennatur von Elektronen und ihrer gegenseitigen abstoßenden Wechselwirkungen (aufgrund gleicher Ladung) auftreten, ohne explizit das Atomgitter und die Struktur eines realen Materials einzuführen. Jellium wird in der Festkörperphysik häufig als einfaches Modell delokalisierter Elektronen in einem Metall verwendet, wo es Eigenschaften realer Metalle wie Abschirmung, Plasmonen, Wigner-Kristallisation und Friedel-Oszillationen qualitativ reproduzieren kann.

Bei Nulltemperatur hängen die Eigenschaften von Jellium ausschließlich von der konstanten Elektronendichte ab. Dies verleiht ihm eine Behandlung innerhalb der Dichtefunktionaltheorie; der Formalismus selbst liefert die Grundlage für die Näherung der lokalen Dichte an das Austauschkorrelations-Energiedichtefunktional.

Der Begriff Gelee wurde von Conyers Herring geprägt und spielt auf den “positiven Jelly”-Hintergrund und das typische metallische Verhalten an, das es zeigt.[1]

Hamiltonian[edit]

Das Jellium-Modell behandelt die Elektron-Elektron-Kopplung rigoros. Die künstliche und strukturlose Hintergrundladung wechselwirkt elektrostatisch mit sich selbst und den Elektronen. Der Jellium-Hamiltonian für n Elektronen, die in einem Raumvolumen eingeschlossen sind Ω, und mit elektronischer Dichte ρ(R) und (konstante) Hintergrundladungsdichte n(R) = n/Ω ist[2][3]

h^=h^el+h^BeinCk+h^el−BeinCk,{displaystyle {hat{H}}={hat{H}}_{mathrm {el}}+{hat{H}}_{mathrm {zurück} }+{hat{H}} _{textrm{el-back}},,}

wo

  • hel ist der elektronische Hamilton-Operator, der aus den kinetischen und den Elektron-Elektronen-Abstoßungstermen besteht:
h^el=Σich=1nPich22m+Σich<Jne2|Rich−RJ|{displaystyle {hat{H}}_{mathrm {el}}=sum_{i=1}^{N}{frac {p_{i}^{2}}{2m}}+ Summe _{i

  • hzurück ist der Hamiltonoperator der elektrostatisch mit sich selbst wechselwirkenden positiven Hintergrundladung:
h^BeinCk=e22∫ΩDR∫ΩDRIch n(R)n(RIch)|R−RIch|=e22(nΩ)2∫ΩDR∫ΩDRIch 1|R−RIch|{displaystyle {hat{H}}_{mathrm {zurück} }={frac {e^{2}}{2}}int _{Omega}mathrm {d} mathbf {R} int_{Omega}mathrm{d}mathbf{R} ‘{frac{n(mathbf{R})n(mathbf{R}’)}{|mathbf{R} – mathbf {R} ‘|}}={frac {e^{2}}{2}}left({frac {N}{Omega}}right)^{2}int _{Omega }mathrm{d}mathbf{R}int_{Omega}mathrm{d}mathbf{R} ‘{frac{1}{|mathbf{R} -mathbf{R}’ |}}}

  • hel-zurück ist der Hamilton-Operator der Elektron-Hintergrund-Wechselwirkung, wiederum eine elektrostatische Wechselwirkung:
h^el−BeinCk=∫ΩDR∫ΩDR ρ(R)n(R)|R−R|=−e2nΩΣich=1n∫ΩDR 1|Rich−R|{displaystyle {hat{H}}_{mathrm {el-back} }=int _{Omega}mathrm {d} mathbf {r} int _{Omega}mathrm {d} mathbf{R}{frac{rho(mathbf{r})n(mathbf{R})}{|mathbf{r} -mathbf{R} |}}=-e^{2 }{frac{N}{Omega}}sum_{i=1}^{N}int_{Omega}mathrm{d}mathbf{R}{frac{1}{| mathbf{r}_{i}-mathbf{R} |}}}

hzurück ist eine Konstante und im Grenzfall eines unendlichen Volumens divergent zusammen mit hel-zurück. Die Divergenz wird durch einen Term aus der Elektron-Elektron-Kopplung aufgehoben: Die Hintergrundwechselwirkungen heben sich auf und das System wird von der kinetischen Energie und der Kopplung der Elektronen dominiert. Eine solche Analyse wird im Fourier-Raum durchgeführt; die verbleibenden Wechselwirkungsterme des Hamilton-Operators entsprechen der Fourier-Entwicklung der Elektronenkopplung, für die Q0.

Beiträge zur Gesamtenergie[edit]

Der traditionelle Weg, das Elektronengas zu untersuchen, besteht darin, mit nicht wechselwirkenden Elektronen zu beginnen, die nur durch den kinetischen Energieteil des Hamilton-Operators bestimmt werden, auch Fermi-Gas genannt. Die kinetische Energie pro Elektron ist gegeben durch

KE=35EF=35ℏ2kF22me=35(9π4)231(RS/ein0)2Ry≈2.21(RS/ein0)2Ry{displaystyle KE={frac {3}{5}}E_{F}={frac {3}{5}}{frac {hbar^{2}k_{F}^{2}}{ 2m_{e}}}={frac {3}{5}}{biggl (}{frac {9pi }{4}}{biggr )}^{frac {2}{3}} {frac {1}{(r_{s}/a_{0})^{2}}}{textrm {Ry}}approx {frac {2.21}{(r_{s}/a_{0}) )^{2}}}{textrm {Ry}}}

wo

EF{displaystyle E_{F}}

ist die Fermi-Energie,

kF{displaystyle k_{F}}

ist der Fermi-Wellenvektor, und der letzte Ausdruck zeigt die Abhängigkeit vom Wigner-Seitz-Radius

RS{displaystyle r_{s}}

wo Energie in Rydbergs gemessen wird.

Ohne viel Arbeit zu leisten, kann man vermuten, dass die Elektron-Elektron-Wechselwirkungen wie die Umkehrung des durchschnittlichen Elektron-Elektron-Abstands skalieren und somit als

1/R12{displaystyle 1/r_{12}}

(da die Coulomb-Wechselwirkung wie eine über den Abstand zwischen Ladungen verläuft), wenn wir die Wechselwirkungen also als kleine Korrektur der kinetischen Energie betrachten, beschreiben wir den Grenzwert von small

RS{displaystyle r_{s}}

(dh

1/RS2{displaystyle 1/r_{s}^{2}}

größer sein als

1/RS{displaystyle 1/r_{s}}

) und damit eine hohe Elektronendichte. Leider haben echte Metalle typischerweise

RS{displaystyle r_{s}}

zwischen 2-5, was bedeutet, dass dieses Bild ernsthaft überarbeitet werden muss.

Die erste Korrektur des Freie-Elektronen-Modells für Jellium stammt vom Fock-Austausch-Beitrag zu Elektron-Elektron-Wechselwirkungen. Addiert man dies hinzu, hat man eine Gesamtenergie von

E=2.21RS2−0,916RS{displaystyle E={frac {2.21}{r_{s}^{2}}}-{frac {0.916}{r_{s}}}}

wobei der negative Term auf Austausch zurückzuführen ist: Austauschwechselwirkungen senken die Gesamtenergie. Korrekturen höherer Ordnung der Gesamtenergie sind auf die Elektronenkorrelation zurückzuführen, und wenn man sich entscheidet, in einer Reihe für kleine . zu arbeiten

RS{displaystyle r_{s}}

, findet man

E=2.21RS2−0,916RS+0,0622ln⁡(RS)−0,096+Ö(RS){displaystyle E={frac {2,21}{r_{s}^{2}}}-{frac {0,916}{r_{s}}}}+0,0622ln(r_{s})-0,096+O (r_{s})}

Die Serie ist ziemlich genau für kleine

RS{displaystyle r_{s}}

aber von zweifelhaftem Wert für

RS{displaystyle r_{s}}

Werte in tatsächlichen Metallen gefunden.

Für die gesamte Palette von

RS{displaystyle r_{s}}

, kann die Korrelationsenergiedichte von Chachiyo als Korrektur höherer Ordnung verwendet werden. In diesem Fall,

E=2.21RS2−0,916RS+einln⁡(1+BRS+BRS2){displaystyle E={frac {2,21}{r_{s}^{2}}}-{frac {0,916}{r_{s}}}+aln left(1+{frac {b }{r_{s}}}+{frac {b}{r_{s}^{2}}}right)}

, [4] was ziemlich gut (in der Größenordnung von Milli-Hartree) mit der Quanten-Monte-Carlo-Simulation übereinstimmt.

Null-Temperatur-Phasendiagramm von Jellium in drei und zwei Dimensionen[edit]

Die Physik des Null-Temperatur-Phasenverhaltens von Jellium wird durch die Konkurrenz zwischen der kinetischen Energie der Elektronen und der Elektron-Elektron-Wechselwirkungsenergie angetrieben. Der kinetische Energieoperator in den Hamiltonschen Skalen als

1/RS2{displaystyle 1/r_{s}^{2}}

, wo

RS{displaystyle r_{s}}

ist der Wigner-Seitz-Radius, während der Wechselwirkungsenergieoperator skaliert als

1/RS{displaystyle 1/r_{s}}

. Daher dominiert die kinetische Energie bei hoher Dichte (klein

RS{displaystyle r_{s}}

), während bei geringer Dichte die Wechselwirkungsenergie dominiert (groß

RS{displaystyle r_{s}}

).

Die Grenze der hohen Dichte ist dort, wo Jellium einem nicht wechselwirkenden freien Elektronengas am ähnlichsten ist. Um die kinetische Energie zu minimieren, werden die Einzelelektronenzustände delokalisiert, in einem Zustand, der sehr nahe an der Slater-Determinante (nicht wechselwirkender Zustand) aus ebenen Wellen liegt. Hier werden die Zustände der ebenen Welle mit dem niedrigsten Impuls doppelt von Spin-up- und Spin-down-Elektronen besetzt, was eine paramagnetische Fermi-Flüssigkeit ergibt.

Bei niedrigeren Dichten, wo die Wechselwirkungsenergie wichtiger ist, ist es energetisch vorteilhaft, dass das Elektronengas spinpolarisiert (dh ein Ungleichgewicht in der Anzahl der Spin-up- und Spin-down-Elektronen aufweist), was zu einem ferromagnetischen Fermi . führt Flüssigkeit. Dieses Phänomen ist bekannt als wandernder Ferromagnetismus. Bei ausreichend niedriger Dichte wird der kinetische Energieverlust, der sich aus der Notwendigkeit ergibt, ebene Wellenzustände mit höherem Impuls zu besetzen, durch die Verringerung der Wechselwirkungsenergie aufgrund der Tatsache, dass Austauscheffekte nicht unterscheidbare Elektronen voneinander fernhalten, mehr als ausgeglichen.

Eine weitere Reduzierung der Wechselwirkungsenergie (auf Kosten der kinetischen Energie) kann durch die Lokalisierung der Elektronenorbitale erreicht werden. Als Ergebnis bildet Jellium bei einer Temperatur von Null und einer ausreichend niedrigen Dichte einen sogenannten Wigner-Kristall, in dem die Einteilchenorbitale eine ungefähr Gaußsche Form haben, die auf Kristallgitterplätzen zentriert ist. Sobald sich ein Wigner-Kristall gebildet hat, kann es bei den Wigner-Kristallen mit abnehmender Dichte grundsätzlich zu weiteren Phasenübergängen zwischen verschiedenen Kristallstrukturen und zwischen verschiedenen magnetischen Zuständen (zB antiferromagnetische zu ferromagnetische Spinkonfigurationen) kommen. Wenn eine Wigner-Kristallisation auftritt, erhält Jellium eine Bandlücke.

Nach der Hartree-Fock-Theorie wird das ferromagnetische Fluid bei einem Dichteparameter von . schlagartig stabiler als das paramagnetische Fluid

RS=5.45{displaystyle r_{s}=5.45}

in drei Dimensionen (3D) und

2.01{displaystyle 2.01}

in zwei Dimensionen (2D).[5] Nach der Hartree-Fock-Theorie tritt die Wigner-Kristallisation jedoch bei

RS=4.5{displaystyle r_{s}=4,5}

in 3D und

1.44{displaystyle 1.44}

in 2D, so dass Jellium kristallisieren würde, bevor wandernder Ferromagnetismus auftritt.[6] Darüber hinaus sagt die Hartree-Fock-Theorie ein exotisches magnetisches Verhalten voraus, wobei die paramagnetische Flüssigkeit gegenüber der Bildung einer spiralförmigen Spindichtewelle instabil ist.[7][8] Leider enthält die Hartree-Fock-Theorie keine Beschreibung von Korrelationseffekten, die energetisch überhaupt wichtig sind, außer den allerhöchsten Dichten, und so ist eine genauere Theorieebene erforderlich, um quantitative Aussagen über das Phasendiagramm von Jellium zu treffen.

Quantum Monte Carlo (QMC)-Methoden, die eine explizite Behandlung von Elektronenkorrelationseffekten ermöglichen, gelten allgemein als der genaueste quantitative Ansatz zur Bestimmung des Null-Temperatur-Phasendiagramms von Jellium. Die erste Anwendung der Diffusions-Monte-Carlo-Methode war die berühmte Berechnung des Null-Temperatur-Phasendiagramms von 3D-Jellium von Ceperley und Alder aus dem Jahr 1980.[9] Sie berechneten, dass der Übergang von paramagnetisch-ferromagnetischer Flüssigkeit bei auftritt

RS=75(5){displaystyle r_{s}=75(5)}

und Wigner-Kristallisation (zu einem kubisch raumzentrierten Kristall) bei

RS=100(20){displaystyle r_{s}=100(20)}

. Nachfolgende QMC-Berechnungen[10][11] haben ihr Phasendiagramm verfeinert: Es gibt einen Übergang zweiter Ordnung von einem paramagnetischen Fluidzustand zu einem partiell spinpolarisierten Fluid von

RS=50(2){displaystyle r_{s}=50(2)}

bis etwa

100{displaystyle 100}

; und Wigner-Kristallisation erfolgt bei

RS=106(1){displaystyle r_{s}=106(1)}

.

In 2D zeigen QMC-Rechnungen, dass der Übergang von paramagnetischer Flüssigkeit zu ferromagnetischer Flüssigkeit und die Wigner-Kristallisation bei ähnlichen Dichteparametern im Bereich

30<RS<40{displaystyle 30

.[12][13] Die neuesten QMC-Rechnungen zeigen, dass es für ein ferromagnetisches Fluid keinen Stabilitätsbereich gibt.[14] Stattdessen findet bei ein Übergang von einer paramagnetischen Flüssigkeit zu einem hexagonalen Wigner-Kristall statt

RS=31(1){displaystyle r_{s}=31(1)}

. Möglicherweise gibt es einen kleinen Stabilitätsbereich für einen (frustrierten) antiferromagnetischen Wigner-Kristall vor einem weiteren Übergang zu einem ferromagnetischen Kristall. Der Kristallisationsübergang in 2D ist nicht erster Ordnung, daher muss es eine kontinuierliche Reihe von Übergängen von Fluid zu Kristall geben, möglicherweise mit gestreiften Kristall/Fluid-Phasen.[15] Experimentelle Ergebnisse für ein 2D-Lochgas in einer GaAs/AlGaAs-Heterostruktur (die, obwohl sie sauber ist, möglicherweise nicht genau dem idealisierten Jellium-Modell entspricht) zeigen eine Wigner-Kristallisationsdichte von

RS=35,1(9){displaystyle r_{s}=35.1(9)}

.[16]

Anwendungen[edit]

Jellium ist das einfachste Modell wechselwirkender Elektronen. Es wird bei der Berechnung von Eigenschaften von Metallen verwendet, wo die Kernelektronen und die Kerne als einheitlicher positiver Hintergrund modelliert werden und die Valenzelektronen mit voller Strenge behandelt werden. Semi-infinite Jellium-Platten werden verwendet, um Oberflächeneigenschaften wie Austrittsarbeit und Oberflächeneffekte wie Adsorption zu untersuchen; in der Nähe von Oberflächen variiert die Elektronendichte oszillierend und fällt im Volumen auf einen konstanten Wert ab.[17][18][19]

In der Dichtefunktionaltheorie wird Jellium bei der Konstruktion der Näherung der lokalen Dichte verwendet, die wiederum eine Komponente komplexerer Austauschkorrelations-Energiefunktionale ist. Aus Quanten-Monte-Carlo-Berechnungen von Jellium wurden genaue Werte der Korrelationsenergiedichte für mehrere Werte der Elektronendichte erhalten,[9] die verwendet wurden, um semiempirische Korrelationsfunktionale zu konstruieren.[20]

Das Jellium-Modell wurde auf Superatome angewendet und in der Kernphysik verwendet.

Siehe auch[edit]

  • Freies Elektronenmodell – ein Modellelektronengas, bei dem die Elektronen mit nichts interagieren.
  • Fast freies Elektronenmodell – ein Modellelektronengas, bei dem die Elektronen nicht miteinander wechselwirken, aber ein (schwaches) Potenzial des Atomgitters spüren.

Verweise[edit]

  1. ^ Hughes, RIG (2006). “Theoretische Praxis: das Bohm-Pines Quartett” (PDF). Perspektiven der Wissenschaft. 14 (4): 457–524. mach:10.1162/posc.2006.14.4.457. S2CID 57569991.
  2. ^ Brutto, EKU; Runge, E.; Heinonen, O. (1991). Vielteilchentheorie. Bristol: Verlag Adam Hilger. S. 79–80. ISBN 978-0-7503-0155-8.
  3. ^ Giuliani, Gabriele; Vignale; Giovanni (2005). Quantentheorie der Elektronenflüssigkeit. Cambridge University Press. pp. 13–16. ISBN 978-0-521-82112-4.
  4. ^ Teepanis Chachiyo (2016). “Einfache und genaue gleichförmige Elektronengas-Korrelationsenergie für den gesamten Dichtebereich”. J.Chem. Physik. 145 (2): 021101. Bibcode:2016JChPh.145b1101C. mach:10.1063/1.4958669. PMID 27421388.
  5. ^ Giuliani, Gabriele; Vignale; Giovanni (2005). Quantentheorie der Elektronenflüssigkeit. Cambridge University Press. ISBN 978-0-521-82112-4.
  6. ^ JR-Weg; MD Towler; RJ-Bedürfnisse (2003). „Uneingeschränkte Hartree-Fock-Theorie der Wigner-Kristalle“. Phys. Rev. B. 68 (4): 045107. ARXIV:0909.5498. Bibcode:2003PhRvB..68d5107T. mach:10.1103/PhysRevB.68.045107. S2CID 8932393.
  7. ^ AW Overhauser (1960). “Riesige Spindichtewellen”. Phys. Rev. Lett. 4 (9): 462–465. Bibcode:1960PhRvL…4..462O. mach:10.1103/PhysRevLett.4.462.
  8. ^ AW Overhauser (1962). „Spindichtewellen in einem Elektronengas“. Phys. Rev. 128 (3): 1437-1452. Bibcode:1962PhRv..128.1437O. mach:10.1103/PhysRev.128.1437.
  9. ^ ein B DM Ceperley; BJ Erle (1980). “Grundzustand des Elektronengases nach einer stochastischen Methode”. Phys. Rev. Lett. (Eingereichtes Manuskript). 45 (7): 566–569. Bibcode:1980PhRvL..45..566C. mach:10.1103/PhysRevLett.45.566.
  10. ^ FH Zong; C. Lin; DM Ceperley (2002). „Spinpolarisation des dreidimensionalen Elektronengases niedriger Dichte“. Phys. Rev. E. 66 (3): 1–7. arXiv:cond-mat/0205339. Bibcode:2002PhRvE..66c6703Z. mach:10.1103/PhysRevE.66.036703. PMID 12366294. S2CID 11606173.
  11. ^ ND Drummond; Z. Radnai; JR-Weg; MD Towler; RJ-Bedürfnisse (2004). „Diffusionsquanten-Monte-Carlo-Studie von dreidimensionalen Wigner-Kristallen“. Phys. Rev. B. 69 (8): 085116. arXiv:0801.0377. Bibcode:2004PhRvB..69h5116D. mach:10.1103/PhysRevB.69.085116. S2CID 18176116.
  12. ^ B. Tanatar; DM Ceperley (1989). „Grundzustand des zweidimensionalen Elektronengases“. Phys. Rev. B. 39 (8): 5005–5016. Bibcode:1989PhRvB..39.5005T. mach:10.1103/PhysRevB.39.5005. PMID 9948889.
  13. ^ F. Rapisarda; G. Senator (1996). “Diffusions-Monte-Carlo-Studie von Elektronen in zweidimensionalen Schichten”. Aust. J. Physik. 49: 161. Bibcode:1996AuJPh..49..161R. mach:10.1071/PH960161.
  14. ^ ND Drummond; RJ-Bedürfnisse (2009). „Phasendiagramm des zweidimensionalen homogenen Elektronengases mit niedriger Dichte“. Phys. Rev. Lett. 102 (12): 126402. ARXIV:1002.2101. Bibcode:2009PhRvL.102l6402D. mach:10.1103/PhysRevLett.102.126402. PMID 19392300. S2CID 35125378.
  15. ^ B. Spivak; SA Kivelson (2004). „Phasen zwischen einer zweidimensionalen Elektronenflüssigkeit und Wigner-Kristall“. Phys. Rev. B. 70 (15): 155114. Bibcode:2004PhRvB..70o5114S. mach:10.1103/PhysRevB.70.155114.
  16. ^ J. Yoon; CCLi; D. Shahar; DC-Tsui; M. Shayegan (1999). “Wigner-Kristallisation und Metall-Isolator-Übergang von zweidimensionalen Löchern in GaAs at
    B=0{displaystyle B=0}

    “. Phys. Rev. Lett. 82 (8): 1744. ARXIV:cond-mat/9807235. Bibcode:1999PhRvL..82.1744Y. mach:10.1103/PhysRevLett.82.1744. S2CID 119371913.
  17. ^ Lang, ND (1969). „Selbstkonsistente Eigenschaften der Elektronenverteilung an einer Metalloberfläche“. Solid-State-Community. 7 (15): 1047–1050. Bibcode:1969SSCom…7.1047L. mach:10.1016/0038-1098(69)90467-0.
  18. ^ Lang, ND; Kohn, W. (1970). „Theorie der Metalloberflächen: Arbeitsfunktion“. Phys. Rev. B. 3 (4): 1215–223. Bibcode:1971PhRvB…3.1215L. mach:10.1103/PhysRevB.3.1215.
  19. ^ Lang, ND; Kohn, W. (1973). „Oberflächen-Dipol-Barrieren in einfachen Metallen“. Phys. Rev. B. 8 (12): 6010–6012. Bibcode:1973PhRvB…8.6010L. mach:10.1103/PhysRevB.8.6010.
  20. ^ Perdew, JP; McMullen, ER; Zunger, Alex (1981). „Dichtefunktionaltheorie der Korrelationsenergie in Atomen und Ionen: Ein einfaches analytisches Modell und eine Herausforderung“. Phys. Rev. A. 23 (6): 2785–2789. Bibcode:1981PhRvA..23.2785P. mach:10.1103/PhysRevA.23.2785.