imported>Mabschaaf |
imported>Π π π |
| Zeile 1: |
Zeile 1: |
| Die klassische '''Dichtefunktionaltheorie''' (DFT) ('''Klassische Dichtefunktionaltheorie''') ist in der [[Statistische Physik|statistischen Physik]] eine Methode, das Verhalten eines Vielteilchensystems (etwa eines [[Gas|Gases]] in einem Behälter) zu beschreiben. Die DFT ist heutzutage eine Standardtechnik in der Flüssigkeitstheorie. Im Gegensatz zur älteren [[Dichtefunktionaltheorie (Quantenphysik)|quantenmechanischen Dichtefunktionaltheorie]] wird sie auf Vielteilchensysteme angewandt, die mit der klassischen Physik beschrieben werden.
| | #WEITERLEITUNG [[Dichtefunktionaltheorie (statistische Physik)]] |
| | |
| Die klassische DFT ermöglicht es für gegebene Parameter (u. a. [[Temperatur]] und von außen vorgegebene [[fundamentale Wechselwirkung|Wechselwirkungen]]) die ortsabhängige [[Teilchendichte|Dichte]] dieses Systems, Korrelationsfunktionen (u. a. [[Radiale Verteilungsfunktion]]) und thermodynamische Eigenschaften ([[Freie Energie]], [[Zustandsgleichung]], [[Phasenübergang|Phasenübergänge]]) zu berechnen. Die Stärke der DFT liegt darin, dass sie für ''inhomogene Systeme'' (ortsabhängige Teilchendichte) anwendbar ist.
| |
| | |
| Die Theorie wurde durch [[Robert Evans (Physiker)|Robert Evans]] (Universität Bristol) 1979 begründet, der das zugrundeliegende Variationsprinzip bewies, wobei er auf ein entsprechendes Prinzip bei quantenmechanischen Vielteilchensystemen bei endlicher Temperatur von [[N. David Mermin]] (1965, eine Verallgemeinerung des [[Hohenberg-Kohn-Theorem]]s) zurückgriff. Die Theorie hat historische Vorläufer in klassischen Untersuchungen von [[Johannes Diderik van der Waals]] über die Flüssig-Gas-Grenzfläche (1893) und von [[Lars Onsager]] über Phasenübergänge in Flüssigkristallen (1949).
| |
| | |
| == Anwendungsgebiete ==
| |
| | |
| Die klassische DFT ermöglicht es thermodynamische Eigenschaften und Korrelationsfunktionen für Systeme zu berechnen, deren Translationsinvarianz und/oder Rotationsinvarianz gebrochen ist. In inhomogenen Situationen können Effekte auftreten, die es in der homogenen Phase nicht gibt. Beispiele hierfür sind:
| |
| | |
| * Einschränkende Geometrien erzeugen ein externes Potential und verändern die Struktur des Fluids und sein Phasenverhalten: Aufgrund von Packungseffekten oszilliert die Dichte in der Nähe einer Wand, während sich in größerer Entfernung die Effekte herausmitteln. An festen Substraten und Wänden kann Adsorption zu geordneten Strukturen in der Nähe der Wand führen. Dabei können makroskopisch messbare Effekte wie Benetzung, Sedimentation, oder Oberflächen-Phasenübergänge auftreten. Bei Mikrokanälen oder porösen Medien können Effekte wie Kapillarkondensation oder capillary filling auftreten.
| |
| | |
| * Grenzflächen zwischen verschiedenen koexistierenden Phasen (flüssig-gasförmig, kristallin-flüssig, flüssig-flüssig-[[Mischungslücke|Entmischung]] bei Flüssigkeitsmischungen, ...). Die Teilchendichte variiert stark im Gebiet der Phasengrenzfläche. Diese Variation führt zum makroskopischen Effekt der Oberflächenspannung.
| |
| | |
| * Äußere Felder, z. B. Sedimentation von suspendierten Kolloiden im Gravitationsfeld, Ausrichtung von Dipolen im elektrischen Feld (Moleküle eines [[Flüssigkristall]]s).
| |
| | |
| Auch die kristalline Phase lässt sich grundsätzlich im Rahmen der DFT behandeln, da man diese als periodische Dichte darstellen kann. Man kann also den Phasenübergang des Gefrierens berechnen.
| |
| | |
| Im Vergleich zu Simulationen ([[Monte-Carlo-Simulation|MC]] oder [[Moleküldynamik|MD]]) ist die DFT-Lösung meist deutlich schneller zu berechnen.
| |
| | |
| Die klassische DFT ist nur im klassischen Limit anwendbar, also nicht dort, wo quantenmechanische Effekte dominieren. Ein Kriterium hierfür lautet, dass die [[thermische Wellenlänge]] viel kleiner als die mittlere nächste-Nachbar-Entfernung sein muss.
| |
| | |
| == Variationsprinzip ==
| |
| | |
| Das Prinzip der Dichtefunktionaltheorie beruht darauf, dass das thermodynamische Potenzial (z. B. die [[freie Energie]] <math>F</math> oder das [[großkanonisches Potenzial|großkanonische Potenzial]] <math>\Omega</math>) eines [[Ensemble (Physik)|Ensembles]] sich als [[Funktional]] der mikroskopischen [[Teilchendichte|Dichte]] <math>\rho</math> schreiben lässt, also <math>F[\rho]</math> bzw. <math>\Omega[\rho]</math>. Dieses Funktional wird minimal bei der physikalisch realisierten Dichte <math>\rho_0</math>, der Gleichgewichtsdichte. Für das großkanonische Potential gilt also
| |
| | |
| :<math>\Omega[\rho]>\Omega[\rho_{0}]\ \text{für}\ \rho\neq\rho_{0}</math>
| |
| | |
| wobei <math>\Omega[\rho_{0}]=\Omega_0</math>. Ist das Dichtefunktional für ein gegebenes System bekannt, kann die Gleichgewichtsdichte also durch Minimieren des Funktionals gefunden werden. Die [[Funktionalableitung]] ist bei der Gleichgewichtsdichte gleich null:
| |
| | |
| :<math>\left.\frac{\delta\Omega[\rho]}{\delta\rho(\boldsymbol{r})}\right|_{\rho_{0}}=0</math>
| |
| | |
| === Herleitung ===
| |
| | |
| Im Folgenden verwendete Symbole: <math>H</math> die [[Hamiltonfunktion]], <math>\mu</math> das [[chemisches Potential|chemische Potential]], <math>N</math> die [[Teilchenzahl]], <math>k_\mathrm{B}</math> die [[Boltzmannkonstante]], <math>T</math> die Temperatur, <math>\beta=(k_\mathrm{B}T)^{-1}</math> die inverse [[thermische Energie]], <math>Z_\text{gk}</math> die [[großkanonische Zustandssumme]], <math>\operatorname{Tr}</math> die klassische Spurbildung im großkanonischen Ensemble.
| |
| | |
| Aus der Statistik des großkanonischen Ensembles ist die Gleichgewichtsdichte und das großkanonische Gleichgewichtspotential bekannt.
| |
| :<math>
| |
| f_{0}=Z_\mathrm{gk}^{-1}\exp\!\left[-\beta(H-\mu N)\right]
| |
| \quad\implies\quad
| |
| \Omega_{0}=-k_{\mathrm{B}}T\ \ln Z_{\text{gk}}
| |
| </math>
| |
| | |
| Das großkanonische Potential kann man allgemein (auch für Nichtgleichgewichtszustände) als Funktional einer beliebigen Wahrscheinlichkeitsdichte <math>f</math> im Phasenraum schreiben:
| |
| | |
| :<math>\Omega_{\text{gk}}[f]=\operatorname{Tr}\left\{ f\left(H-\mu N+k_\mathrm{B}T\ln f\right)\right\}</math>
| |
| | |
| Man kann mit einer [[Gibbs-Ungleichung]] beweisen, dass die Gleichgewichtsdichte <math>f_0</math> das Funktional minimiert:
| |
| | |
| :<math>\Omega_0=\min_{f}\Omega_{\text{gk}}[f] </math>
| |
| | |
| Der entscheidende Schritt ist der Übergang von einem Funktional von <math>f</math> zu einem Funktional von <math>\rho</math>. Dabei hängt <math>f</math> von <math>6N</math> Koordinaten ab (für <math>N</math> Teilchen jeweils <math>3</math> Impulskoordinaten <math>\boldsymbol{p}_{i}</math> und <math>3</math> Ortskoordinaten <math>\boldsymbol{r}_{i}</math>) dagegen <math>\rho</math> nur von <math>3</math> Ortskoordinaten <math>\boldsymbol{r}</math>. Der Zusammenhang ist gegeben durch: <math>\textstyle \rho(\boldsymbol{r})=\operatorname{Tr}\{ f \sum_{i=1}^{N}\delta(\boldsymbol{r}-\boldsymbol{r}_{i}) \}</math>.
| |
| | |
| Die Minimierung nach <math>f</math> kann in eine doppelte Minimierung umgeschrieben werden (Levy-Methode):
| |
| | |
| :<math>\Omega_{0}=\min_{\rho}\Omega[\rho] \quad\text{wobei}\quad \Omega[\rho]=\min_{f\to\rho}\Omega_{\text{gk}}[f]</math>
| |
| | |
| Die innere Minimierung bedeutet, dass <math>\Omega_{\text{gk}}[f]</math> minimiert wird unter der Bedingung, dass ein bestimmtes <math>\rho</math> durch <math>f</math> generiert wird. Aus dieser Minimierung lässt sich jedoch im Allgemeinen kein analytischer Ausdruck für (den Exzessteil von) <math>\Omega[\rho]</math> ableiten.
| |
| | |
| == Intrinsische Freie Energie ==
| |
| | |
| Die Hamiltonfunktion kann man aufspalten <math>H=K_{N}+\Phi_{N}+V_{N}</math> in kinetische Energie, interne Wechselwirkung und externe Wechselwirkung. Diese lauten:
| |
| :<math>\textstyle K_{N}=\sum_{i=1}^{N}\frac{\boldsymbol{p}_{i}^{2}}{2m}\,,\quad\Phi_{N}=\sum_{i<j}^{N}\phi(\boldsymbol{r}_{i},\boldsymbol{r}_{j})\,,\quad V_{N}=\sum_{i=1}^{N}V(\boldsymbol{r}_{i})</math>
| |
| | |
| Hier ist <math>\phi(\boldsymbol{r}_{i},\boldsymbol{r}_{j})</math> das effektive interne Potential zwischen den Flüssigkeitsteilchen und <math>V(\boldsymbol{r}_i)</math> das externe Potential.
| |
| | |
| Einsetzen und Aufspalten:
| |
| :<math>\Omega[\rho]=\underbrace{\underbrace{\overbrace{\min_{f\to\rho}\operatorname{Tr}\left\{ f\left(K_{N}+k_{\mathrm{B}}T\ln f\right)\right\} }^{\mathcal{F}_{\text{id}}[\rho]}+\overbrace{\min_{f\to\rho}\operatorname{Tr}\left\{ f\Phi_{N}\right\} }^{\mathcal{F}_{\text{ex}}[\rho]}}_{\mathcal{F}[\rho]}+\overbrace{\min_{f\to\rho}\operatorname{Tr}\left\{ fV_{N}\right\} }^{\int\mathrm{d}\boldsymbol{r}\,\rho(\boldsymbol{r})\, V(\boldsymbol{r})}}_{F[\rho]}+\overbrace{\min_{f\to\rho}\operatorname{Tr}\left\{ f\left(-\mu N\right)\right\} }^{-\int\mathrm{d}\boldsymbol{r}\,\rho(\boldsymbol{r})\,\mu}</math>
| |
| | |
| Hierbei wurden die Funktionale freie Energie <math>F[\rho]</math> und intrinsische freie Energie <math>\mathcal{F}[\rho]</math> definiert. Man führt das intrinsische Freie-Energie-Funktional derart ein, dass es nur von internen Wechselwirkungen <math>\phi</math>, nicht aber von externen Wechselwirkungen <math>V</math> abhängt. Somit ist die analytische Form von <math>\mathcal{F}[\rho]</math> auf alle Inhomogenitäten anwendbar, da diese von <math>V</math> hervorgerufen werden.
| |
| | |
| Der Zusammenhang zwischen <math>\Omega[\rho]</math> und <math>\mathcal{F}[\rho]</math> ist also gegeben durch
| |
| | |
| :<math>\Omega[\rho]=\mathcal{F}[\rho]-\int\mathrm{d}\boldsymbol{r}\,\rho(\boldsymbol{r})\mu_{\text{in}}(\boldsymbol{r})</math>
| |
| | |
| wobei <math>\mu_{\text{in}}(\boldsymbol{r}):=\mu-V(\boldsymbol{r})</math> das intrinsische chemische Potential definiert. Dies entspricht einer [[Legendre-Transformation]] zwischen den thermodynamischen Potentialen.
| |
| | |
| === Ideal und Exzess-Funktional ===
| |
| | |
| Man teilt die intrinsische freie Energie in einen idealen und einen exzess Teil auf. Ersterer beschreibt den wechselwirkungsfreien Anteil (siehe [[ideales Gas]]), letzterer beschreibt die Wechselwirkungen innerhalb der Flüssigkeit.
| |
| | |
| :<math>\mathcal{F}[\rho] =\mathcal{F}_{\text{id}}[\rho]+\mathcal{F}_{\text{ex}}[\rho]</math>
| |
| | |
| Der ideale Teil lässt sich analytisch exakt berechnen (<math>\Lambda</math> ist die [[thermische Wellenlänge]]):
| |
| | |
| :<math>\mathcal{F}_{\text{id}}[\rho]=k_\mathrm{B}T\int\mathrm{d}\boldsymbol{r}\,\rho(\boldsymbol{r})\left\{ \ln\!\left[\Lambda^{3}\rho(\boldsymbol{r})\right]-1\right\} </math>
| |
| | |
| Das Exzess-Funktional hängt von dem jeweiligen internen Wechselwirkungspotential <math>\phi</math> ab und ist im Allgemeinen unbekannt. Man definiert die direkten Korrelationsfunktionen als Funktionalableitung:
| |
| | |
| :<math>c^{(n)}(\boldsymbol{r}_{1},\ldots,\boldsymbol{r}_{n},[\rho]):=-\beta\left.\frac{\delta^{n}\mathcal{F}_{\text{ex}}[\rho]}{\delta\rho(\boldsymbol{r}_{1})\ldots\delta\rho(\boldsymbol{r}_{n})}\right|_{\rho}</math>
| |
| | |
| Für das ideale Gas ist <math>\Phi=0</math> und somit <math>\mathcal{F}_{\text{ex}}=0</math> und <math>c^{(n)}=0</math>
| |
| | |
| == Anwendung des Variationsprinzips ==
| |
| | |
| Für obiges Variationsprinzip gilt:
| |
| | |
| :<math>0=\left.\frac{\delta\Omega[\rho]}{\delta\rho(\boldsymbol{r})}\right|_{\rho=\rho_{0}}=
| |
| k_\mathrm{B}T\ \ln\!\left[\Lambda^{3}\rho_{0}(\boldsymbol{r})\right]-k_\mathrm{B}T\ c^{(1)}(\boldsymbol{r},[\rho_{0}])-\mu_{\text{in}}(\boldsymbol{r})</math>
| |
| | |
| Durch Auflösen nach der Dichte erhält man die verallgemeinerte [[barometrische Höhenformel]]
| |
| | |
| :<math>\rho_{0}(\boldsymbol{r})=\Lambda^{-3}\exp\!\left[\beta\mu_{\text{in}}(\boldsymbol{r})+c^{(1)}(\boldsymbol{r},[\rho_{0}])\right]</math>
| |
| | |
| Zur barometrische Höhenformel des idealen Gases kommt die Funktion <math>c^{(1)}(\boldsymbol{r})</math> hinzu, die also den Einfluss der Teilchenwechselwirkung auf das Dichteprofil beinhaltet. Schreibt man den Exponenten als <math>\beta\mu-\beta\left\{ V(\boldsymbol{r})-k_{\mathrm{B}}T\, c^{(1)}(\boldsymbol{r})\right\}</math>, sieht man, dass das externe Potential um <math>-k_{\mathrm{B}}T\, c^{(1)}(\boldsymbol{r})</math> moduliert wird.
| |
| | |
| Da <math>c^{(1)}(\boldsymbol{r})</math> ein Funktional der Dichte ist, kann man die verallgemeinerte barometrische Höhenformel verwenden die Gleichung mittels [[Fixpunktiteration]] selbstkonsistent zu lösen.
| |
| | |
| == Näherungen für das Exzess-Funktional ==
| |
| | |
| Unterschiedliche Wechselwirkungen zwischen den betrachteten Teilchen (z. B. [[Lennard-Jones-Potential]], harte Kugeln, weiche Repulsion zwischen Polymer-Knäuel) erfordern verschiedene Exzess-Funktionale. Ist jedoch für eine bestimmte (interne) Wechselwirkung ein Funktional bekannt, lassen sich damit sämtliche inhomogene Situationen (für sämtliche externe Potentiale) berechnen.
| |
| | |
| Nur für harte Stäbchen in einer Dimension lässt sich das Exzess-Funktional exakt konstruieren, für alle anderen wechselwirkenden Systeme müssen geeignete Approximationen verwendet werden. Somit liegt das zentrale Problem der DFT in der Beschaffung einer geeigneten Näherung für dieses Funktional. Die Funktionalentwicklung ausgehend von mikroskopischen Eigenschaften (von einer effektiven Hamiltonfunktion) erfordert viel Erfahrung. Es gibt jedoch einige Standardfunktionale, die sehr vielseitig einsetzbar sind.
| |
| | |
| Häufig verwendete Näherungen sind:
| |
| * [[Lokale Dichtenäherung]] (Local Density Approximation, LDA)
| |
| * [[Molekularfeldnäherung]] (Mean Field Approximation, MFA bzw. Random Phase Approximation, RPA)
| |
| * [[Gewichtete Dichtenäherung]] (Weighted Density Approximation, WDA)
| |
| * [[Rosenfeld-Funktionale]] (Fundamental measure theory, FMT)
| |
| | |
| == Dynamische Dichtefunktionaltheorie ==
| |
| Neben der DFT, die Gleichgewichtszustände betrachtet, gibt es für Nichtgleichgewichtszustände auch die DDFT (dynamische DFT), mit der man die zeitliche Entwicklung eines Systems berechnen kann (z. B. [[Kolloid|kolloidale Suspensionen]], die der [[Brownsche Bewegung|Brownschen Bewegung]] unterliegen).
| |
| | |
| == Literatur ==
| |
| * J.-P. Hansen, I. R. McDonald: ''Theory of Simple Liquids''. Academic Press/Elsevier, 4. Auflage, 2013, ISBN 978-0-12-387032-2
| |
| * D. Henderson: ''Fundamentals of Inhomogeneous Fluids''. Dekker, 1992, ISBN 978-0-82-478711-0, Chap. 3 by R. Evans.
| |
| * R. Evans: ''The nature of the liquid-vapour interface and other topics in the statistical mechanics of non-uniform, classical fluids''. Adv. Phys. '''28''', 143-200 (1979). {{doi|10.1080/00018737900101365}}
| |
| * Y. Rosenfeld: ''Free-energy model for the inhomogeneous hard-sphere fluid mixture and density-functional theory of freezing''. Phys. Rev. Lett. '''63''', 980–983 (1989). {{doi|10.1103/PhysRevLett.63.980}}
| |
| *M. Schmidt, M. Burgis, W. S. B. Dwandaru, G. Leithall, P. Hopkins: Recent developments in classical density functional theory: Internal energy functional and diagrammatic structure of fundamental measure theory, Condensed Matter Physics, Band 15, 2002, 43603, 1-15, [http://www.icmp.lviv.ua/journal/zbirnyk.72/43603/art43603.pdf pdf]
| |
| | |
| == Siehe auch ==
| |
| * [[Dichtefunktionaltheorie (Quantenphysik)]] und [[Hohenberg-Kohn-Theorem]].
| |
| | |
| ==Weblinks==
| |
| *[http://agenda.albanova.se/getFile.py/access?contribId=260&resId=251&materialId=250&confId=2509 R. Evans, Density Functional Theory for Inhomogeneous Liquids I, 2009, pdf]
| |
| | |
| [[Kategorie:Statistische Physik]]
| |
| [[Kategorie:Variationsrechnung]]
| |