UNIFAC[1][2] (Abkürzung für Universal Quasichemical Functional Group Activity Coefficients) ist ein Verfahren zur Abschätzung von Aktivitätskoeffizienten, das in der Verfahrenstechnik und der Technischen Chemie häufig eingesetzt wird.
UNIFAC ist von UNIQUAC abgeleitet. UNIQUAC verwendet stoffspezifische Parameter, die aus experimentell ermittelten Daten bestimmt werden müssen. Diese Parameter werden in UNIFAC vorhergesagt, so dass keine gemessenen Daten mehr benötigt werden.
UNIFAC ist eine Gruppenbeitragsmethode, die auf dem Prinzip der Mischung aus Strukturgruppen basiert. Dies steht im Gegensatz zu der normalen Betrachtungsweise der Mischung aus Molekülen. Die Gruppen sind dabei üblicherweise funktionelle Gruppen, bspw. Alkohol- oder Carbonylgruppe, aber auch kleinere Molekülfragmente wie einzelne Atome, jedoch werden diese kleinen Fragmente stets unter Berücksichtigung ihrer chemischen Umgebung betrachtet. Einige wenige kleinere Moleküle (etwa Wasser) werden vollständig als eigenständige Gruppe definiert.
Das Prinzip der Gruppenbeitragsmethoden besteht im Wesentlichen darin, dass nur die Eigenschaften einiger Dutzend Strukturgruppen bekannt sein müssen, statt die Eigenschaften von einigen Millionen Stoffen für die Berechnung von Stoffgemischeigenschaften zu benötigen. Mit diesen wenigen Strukturgruppen als Bausteinen kann eine sehr große Anzahl von Molekülen konstruiert werden.
UNIFAC benutzt zur Vorhersage zwei Arten von Gruppenbeiträgen:
Die Gruppenvolumina- und oberflächen basieren auf angenäherten van der Waalsschen Oberflächen und Volumina[3] und sind somit Konstanten im Modell, die einen physikalisch-chemischen Hintergrund haben.
Die Wechselwirkungsparameter werden an experimentelle Aktivitätskoeffizienten sowie an Phasengleichgewichtsdaten, aus denen Aktivitätskoeffizienten abgeleitet werden können, mit nichtlinearen Optimierverfahren angepasst. Die Wechselwirkungsparameter werden somit ausschließlich empirisch ermittelt.
Das UNIFAC-Modell erlaubt die Vorhersage von Aktivitätskoeffizienten in Gemischen mit Alkanen, Alkenen, Alkinen, Alkoholen, Aromaten, Estern, Ethern, Aminen, Carbonsäuren, organischen Fluoriden, organischen Chloriden, organischen Bromiden, organischen Iodiden, Thiolen, Sulfonen, Wasser, Furfural, Thiophenen, Pyridinen, Morpholin, Isocyanaten, Silanen, Siloxanen, Amiden und organischen Nitraten.
Für diese funktionellen Gruppen sind Wechselwirkungsparameter mit mindestens einer der anderen Gruppen vorhanden. Es sind jedoch nicht alle Kombinationen vollständig parametrisiert, bspw. kann ein Aktivitätskoeffizient eines Gemisch aus 3-Methylthiophen und n-Hexan berechnet werden, eine Vorhersage eines Aktivitätskoeffizienten in einem Gemisch des Thiophens mit 3-Hexen scheitert jedoch an fehlenden Parametern.
UNIFAC kann nicht verwendet werden, um Aktivitätskoeffizienten in Gemischen zu berechnen, die Salze bzw., allgemeiner, Elektrolyte oder Polymere enthalten. Für diese beiden Stoffklassen gibt es aber Erweiterungen des hier beschriebenen UNIFAC-Modells.
Ausführliche Formulierungen können im Originalartikel gefunden werden oder auch in Lehrbüchern.[5][6]
Die Berechnung des Aktivitätskoeffizienten erfolgt additiv über zwei Terme:
$ \ln \gamma _{i}=\ln \gamma _{i}^{C}+\ln \gamma _{i}^{R} $
γC wird als kombinatorischer Anteil bezeichnet und wird aus van der Waalsschen Oberflächen (F) und Volumina (V) sowie Molenbrüchen (x) berechnet.
$ \ln \gamma _{i}^{C}=1-V_{i}+\ln V_{i}-5q_{i}\left(1-{{V_{i}} \over {F_{i}}}+\ln {{V_{i}} \over {F_{i}}}\right) $
mit
$ V_{i}={{r_{i}} \over {\sum _{j}{r_{j}x_{j}}}} $
$ F_{i}={{q_{i}} \over {\sum _{j}{q_{j}x_{j}}}} $
Neben den Molenbrüchen werden van der Waalsche Oberflächen qi und Volumina ri der Moleküle benötigt. Diese lassen sich aus tabellierten Werten für die Gruppen (Gruppenoberfläche Qk und Gruppenvolumen Rk) bestimmen.
$ q_{i}=\sum {\nu _{k}^{(i)}Q_{k}} $
$ r_{i}=\sum {\nu _{k}^{(i)}R_{k}} $
νk(i) ist die Häufigkeit der Gruppe k im Molekül i.
γR wird als Restanteil bezeichnet und wird schlussendlich aus angepassten Wechselwirkungsparametern berechnet.
Der Restanteil wird aus Gruppenaktivitätskoeffizienten Γk berechnet.
$ \ln \gamma _{i}^{R}=\sum _{k}{\nu _{k}^{(i)}\left(\ln \Gamma _{k}-\ln \Gamma _{k}^{(i)}\right)} $
Die Gruppenaktivitätskoeffizienten in der Mischung Γk und im Reinstoff Γk(i) werden über die Beziehung
$ \ln \Gamma _{k}=Q_{k}\left[1-\ln \left(\sum _{m}{\Theta _{m}\Psi _{mk}}\right)-\sum _{m}{{\Theta _{m}\Psi _{km}} \over {\sum _{n}{\Theta _{n}\Psi _{nm}}}}\right] $
berechnet.
Θ ist der Oberflächenanteil
$ \Theta _{m}={{Q_{m}X_{m}} \over {\sum _{m}{Q_{m}X_{m}}}} $
und X der Gruppenmolenbruch.
$ X_{m}={{\sum _{j}{\nu _{m}^{(j)}x_{j}}} \over {\sum _{j}\sum _{m}{\nu _{n}^{(j)}x_{j}}}} $
Der Parameter ψ enthält die anpassbaren Parameter des UNIFAC-Modells.
$ \Psi _{nm}=\exp {\frac {-a_{nm}}{T}} $ mit anm als Wechselwirkungsparameter zwischen den Gruppen.
UNIFAC wurde in den 1970er Jahren entwickelt, wobei der Schwerpunkt zu Beginn alleine auf die Vorhersage von Dampf-Flüssig-Gleichgewichten von Mischungen einfacher, im Wesentlichen organischer Substanzen und von Wassergemischen gelegt wurde. Da für die Modellentwicklung Aktivitätskoeffizienten für eine Vielzahl von Gemischen benötigt wurden, wurde mit dem Aufbau einer Faktendatenbank (der Dortmunder Datenbank) begonnen, die noch heute existiert, wenn auch in stark erweiterter und veränderter Form.
Da Aktivitätskoeffizienten es auch erlauben, mit einfachen thermodynamischen Beziehungen Fest-Flüssig-Gleichgewichte (Löslichkeiten von Feststoffen in Flüssigkeiten) und Flüssig-Flüssig-Gleichgewichte (siehe hierzu Mischungslücke) zu berechnen, wurde das UNIFAC-Modell in den achtziger Jahren auch zunehmend für die Berechnung dieser Phasengleichgewichte verwendet. Die dabei zutage tretenden Modellschwächen des Originalmodells führte zum einen zur Entwicklung von speziell parametrisierten Modellen für die Berechnung von Flüssig-Flüssig-Gleichgewichten (UNIFAC-LLE) und später zu einem Modell für die Abschätzung von Octanol-Wasser-Verteilungskoeffizienten.[7] Zum anderen wurde das Modell selbst erweitert, bspw. um temperaturabhängige Wechselwirkungsparameter:
$ \psi =\exp \left({{\frac {-a_{nm}}{T}}+{\frac {-b_{nm}\cdot T}{T}}+{\frac {-c_{nm}\cdot T^{2}}{T}}}\right) $
Diese erweiterten Modelle werden als modifizierte UNIFAC-Modelle[8] bezeichnet.
UNIFAC (insbesondere die neueren Varianten) ist heute das am häufigsten verwendete Verfahren zur Abschätzung von Stoffgemischdaten. Die Vorhersagen werden bspw. in der Prozesssimulation verwendet, dem heute vorherrschenden Verfahren zur Auslegung und Optimierung chemischer Prozesse, Anlagen und ganzer Fabriken. UNIFAC wird außerdem bei der Prozesssynthese verwendet, bei dem, sehr allgemein formuliert, Stoffe mit bestimmten Eigenschaften für bestimmte Aufgaben gesucht werden. Diese Aufgabe kann bspw. die eines Schleppmittels für die Azeotrop- oder Extraktivrektifikation sein.
UNIFAC wird in einer Reihe von Forschungsgruppen genutzt und weiterentwickelt. Aktuelle Entwicklungen (Auswahl) zielen auf