Voir la source de Unifac

[[Image:SattleAzeotropACMChemSep.gif|thumb|250px|« Selle » [[azéotrope]] calculée avec UNIFAC à {{unité|1|atm}}. Les lignes rouges sont des compositions de vapeur et les lignes bleues sont des compositions liquides. L'image est en rotation pour montrer plus clairement la forme en selle des équilibres vapeur-liquide.]]
La méthode '''''UNI''versal ''F''unctional ''A''ctivity ''C''oefficient''' (UNIFAC en anglais, en français ''C''oefficient d'''A''ctivité ''F''onctionnel ''UNI''versel) est un système [[Méthode quantique semi-empirique|semi-empirique]] de prévision d'une [[activité chimique]] non-électrolyte dans des mélanges non idéaux.

UNIFAC utilise le [[groupe fonctionnel]] présent sur les molécules qui constituent le mélange liquide pour calculer les coefficients d'activité. En utilisant les interactions de chacun des groupes fonctionnels présents sur les molécules, ainsi que des coefficients d'interaction binaires, l'activité de chacune des solutions peut être calculée. Ces informations peuvent être utilisées pour obtenir des renseignements sur les équilibres liquides, ce qui est utile dans de nombreux calculs thermodynamiques, tels que la conception d'un [[Réacteur (chimie)|réacteur chimique]], et les calculs de [[distillation]].

Le [[Modèle mathématique|modèle]] UNIFAC a d'abord été publié en 1975 par Fredenslund, Jones et Prausnitz, un groupe de chercheurs en [[génie chimique]] de l'[[université de Californie]]. Par la suite, ils ont publié, tout comme d'autres auteurs, une large gamme d'articles sur UNIFAC, étendant les capacités du modèle. Pour cela, ils ont dû développer les paramètres existants du modèle UNIFAC ou bien en créer de  nouveaux. UNIFAC est une tentative de ces chercheurs de fournir un modèle souple des équilibres liquides pour une utilisation plus large en [[chimie]], ainsi que dans des disciplines telles que le génie chimique ou encore le [[génie des procédés]].

== Introduction ==
Un problème particulier dans le domaine de l'état liquide [[thermodynamique]] est l'approvisionnement des constantes thermodynamiques fiables. Ces constantes sont nécessaires pour une prédiction réussie de l'état d'[[énergie libre]] du système. Sans cette information, il est impossible de modéliser les [[équilibre chimique|équilibres]] des phases du système.

L'obtention de ces données sur l'énergie libre n'est pas un problème commun et nécessite des expériences minutieuses, comme la [[calorimétrie]], pour réussir à mesurer l'énergie du système. Même lorsque ce travail est effectué, il est impossible de tenter de l'effectuer pour chaque classe de produits chimiques possible, en plus du binaire, ou même encore des mélanges de ceux-ci. Pour pallier ce problème, les modèles de prévision de l'énergie libre, tels qu'UNIFAC, sont utilisés pour prédire l'énergie du système reposant sur quelques constantes mesurées en amont.

Bien qu'il soit théoriquement possible de calculer certains de ces paramètres en utilisant les méthodes de simulation sur ordinateur [[Méthode ab initio de chimie quantique|''ab initio'']], plusieurs problèmes majeurs avec cette approche existent. D'une part et principalement, les ressources de calcul pour ces raisonnements sont immenses, et mis  à l'échelle, cela est extrêmement défavorable pour les systèmes avec plusieurs atomes. Deuxièmement, les énergies obtenues à partir de ces calculs provenant des simulations ''ab initio'', ont souvent besoin d'une vérification expérimentale pour confirmer leurs résultats. Enfin, ces calculs exigent un degré élevé d'expertise et une bonne compréhension de la [[chimie quantique]]. D'où la nécessité, pour les modèles simplifiés comme UNIFAC, de continuer de réussir à prédire l'état thermodynamique du système.

== La corrélation UNIFAC == 
La corrélation UNIFAC tente de décomposer le problème de prévision des interactions entre les molécules, en décrivant les interactions moléculaires basées sur les groupes fonctionnels attachés à la molécule. Cette démarche est faite dans le but de réduire le grand nombre d'interactions binaires qu'il faudrait mesurer pour prévoir l'état du système.

=== Activité chimique ===
{{article général|Activité chimique}}
Le [[Activité chimique|coefficient d'activité]] des composantes d'un système est un facteur de correction qui tient compte des écarts de systèmes réels par rapport à ceux d'une [[solution idéale]]. Il peut à la fois être mesuré par une expérience ou bien estimé à partir de modèles chimiques (tels que UNIFAC). En ajoutant un facteur de correction - l'activité de la fraction en phase liquide d'un mélange liquide (<math>a_i</math>, l'activité de la composante i<sup>th</sup>) - une partie des effets de la vraie solution ne peut être prise en compte. L'activité d'un produit chimique réel est une [[Fonction d'état|fonction de l'état]] thermodynamique du système, à savoir la température et la pression.

Il est possible de calculer des phénomènes tels que la [[Modèle d'Ising|séparation de phases]] et les [[UNIQUAC|équilibres liquide-vapeur]], en étant équipé des coefficients d'activité et en ayant une connaissance des composantes et de leurs quantités relatives. UNIFAC s'efforce d'être un modèle général pour la prédiction réussie des coefficients d'activité.

=== Les paramètres du modèle ===
Le modèle UNIFAC divise en deux composants le coefficient d'activité pour chaque espèce du système : un combinatoire <math>\gamma ^c</math> et un composant résiduel <math>\gamma^r </math>. Pour la <math>i</math>{{e}} molécule, les coefficients d'activité sont ventilés selon l'équation suivante :
	
<math>\ln \gamma_i =  \ln \gamma_i^c + \ln \gamma_i^r</math>.

Dans le modèle UNIFAC, il existe trois paramètres principaux, nécessaires pour déterminer l'activité de chaque molécule dans le système. Tout d'abord, il y a l'aire de la surface du groupe <math>R</math> et le volume des contributions <math>Q</math> obtenu à partir de l'aire de la surface et des volumes de [[rayon de van der Waals|van der Waals]]. Ces paramètres reposent uniquement sur les groupes fonctionnels individuels et sur les molécules hôtes. Enfin, il y a le paramètre d'interaction binaire <math>\tau_{ij}</math>, qui est lié à l'énergie d'interaction de paires moléculaires <math>U_i</math> (l'équation se trouve dans la section « résiduel »). Ces paramètres doivent être obtenus soit par le biais d'expériences, ''via'' le montage de données, soit par la simulation moléculaire.

==== Combinatoire ==== 
Le composant [[combinatoire]] de l'activité est composé par plusieurs termes de son équation (ci-dessous). Il est le même que pour le modèle [[UNIQUAC]].

<math>\ln \gamma_i^c = \ln \frac{\phi_i}{x_i} + \frac{z}{2} q_i \ln\frac{\theta_i}{\phi_i} + L_i - \frac{\phi_i}{x_i}\displaystyle\sum_{j=1}^{n} x_j L_j </math>

Ici, <math>\theta_i</math> et <math>\phi_i</math> sont le segment molaire pondérée et la zone ''fractionnaire'' composant la <math>i</math>{{e}} molécule dans l'ensemble du système. Ils sont définis par l'équation suivante, où <math>L_i</math> est un paramètre composé de <math>r</math>, <math>z</math> et <math>q</math>. <math>z</math> est le [[nombre de coordination effectif|numéro de la coordination]] du système, mais le modèle se trouve être relativement insensible à sa valeur, par conséquent, il est souvent cité comme une constante ayant la valeur de 10.

<math>\theta_i = \frac{x_i q_i}{\displaystyle\sum_{j=1}^{n} x_j q_j} \mathrm{\,\,;\,\,}  \phi_i = \frac{x_i r_i}{\displaystyle\sum_{j=1}^{n} x_j r_j} \mathrm{\,\,;\,\,} L_i = \frac{z}{2}(r_i - q_i)-(r_i-1)\mathrm{\,\,;\,\,} z=10</math>

Ici, <math>q_i</math> et <math>r_i</math> sont calculés à partir de la surface du groupe et des contributions volumiques <math>Q</math> et <math>R</math> (habituellement obtenues par valeurs tabulées) ainsi que par le nombre d'occurrences du groupe fonctionnel sur chaque molécule <math>\nu_k</math> de telle sorte que :

<math>r_i = \displaystyle\sum_{k=1}^{n} \nu_k R_k \mathrm{\,\,;\,\,} q_i = \displaystyle\sum_{k=1}^{n}\nu_k Q_k </math>.

==== Résiduel ====
Le composant résiduel de l'activité <math>\gamma^{r}</math> est dû à des interactions entre les groupes présents dans le système, et le document original se référant au concept d'une « solution-de-groupes ». Le composant résiduel de l'activité pour la molécule <math>i</math> contenant <math>n</math> groupes fonctionnels et uniques, peut être écrit comme suit :

<math> \ln \gamma_i^r = \displaystyle\sum_{k}^n \nu_k^{(i)} \left[ \ln \Gamma_k - \ln \Gamma_k^{(i)} \right] </math>

où, <math>\Gamma_k^{(i)}</math> est l'activité d'un groupe isolé dans une solution constituée uniquement de molécules de type <math>i</math>. La formulation de l'activité résiduelle assure que, dans le cas limité à une seule molécule dans une solution de composante pure, l'activité est égale à 1 ; que par la définition de <math>\Gamma_k^{(i)}</math>, on constate que <math> \ln \Gamma_k - \ln \Gamma_k^{(i)}</math> est égal à zéro. La formule suivante est utilisée à la fois pour <math>\Gamma_k</math> et pour <math>\Gamma_k^{(i)}</math> :

<math>\ln \Gamma_k = Q_k \left[1 - \ln \displaystyle\sum_m \Theta_m \Psi_{mk} - \displaystyle\sum_m \frac{\Theta_m \Psi_{km}}{\displaystyle\sum_n \Theta_n \Psi_{nm}}\right]</math>.

Dans cette formule, <math>\Theta_m</math> est la somme de la fraction de la surface du groupe <math>m</math>, sur l'ensemble des différents groupes.  <math>\Theta_m</math>, est aussi à peu près semblable dans la forme, mais il n'est pas le même que <math>\theta_i</math>. <math>\Psi_{mn}</math> est le paramètre d'interaction du groupe ainsi qu'une mesure de l'énergie d'interaction entre les groupes. Il est calculé en utilisant une [[loi d'Arrhenius|équation d'Arrhenius]] (mais avec une pseudo-constante de valeur 1). <math>X_n</math> est la [[fraction molaire]] du groupe, qui est le nombre de groupes <math>n</math> dans la solution, divisé par le nombre total de groupe.

<math>\Theta_m = \frac{Q_m X_m}{\displaystyle\sum_{n} Q_n X_n} </math>

<math>\Psi_{mn} = \mathrm{exp} -\left[ \frac{U_{mn} - U_{nm} }{RT} \right]\mathrm{\,;\,\,}</math>  <math>X_m = \frac{ \displaystyle\sum_j \nu^j_m x_j}{\displaystyle\sum_j \displaystyle\sum_n \nu_n^j x_j}</math>

<math>U_{mn}</math> est l'énergie d'interaction entre les groupes de <math>m</math> et <math>n</math>, exprimée en joules par mole (unités [[Système international d'unités|SI]]), et <math>R</math> est la [[constante universelle des gaz parfaits]]. Notons que ce n'est pas le cas de <math>U_{mn} = U_{nm}</math>, donnant lieu à un paramètre {{quoi|non-réflexive}}. L'équation pour le paramètre d'interaction du groupe peut être simplifiée comme suit :

<math>\Psi_{mn} = \mathrm{exp}\frac{-a_{mn}}{T}</math>.

Ainsi, <math>a_{mn}</math> représente toujours l'énergie nette de l'interaction entre les groupes <math>m</math> et <math>n</math>, mais il a des unités quelque peu inhabituelles de [[température  absolue]] (SI de [[kelvin]]). Ces valeurs d'énergie d'interaction sont obtenues à partir des données expérimentales, et sont généralement sous forme de tableaux.

== Références ==
* ''Group-Contribution Estimation of Activity Coefficients in Nonideal Liquid Mixtures'', Journal de l'institut américain de l'ingénierie chimique, {{vol.}}21, {{n°}}6, {{date-|novembre 1975}}, {{p.}}{{formatnum:1086}}-{{formatnum:1099}}.

== Voir aussi ==
=== Articles connexes ===
* [[Thermodynamique]]
* [[Activité chimique]]
* [[UNIQUAC]]
* [[Cosmospace|COSMOSPACE]]
* {{Lien|trad=Non-random two-liquid model|texte=NRTL}}

=== Lecture complémentaire ===
* Équilibre vapeur-liquide utilisant UNIFAC : Aage Fredenslund, {{Lien|fr=Jürgen Gmehling|lang=en|trad=Jürgen Gmehling|texte=Jürgen Gmehling}} et Peter Rasmussen, {{Lien|fr=Group contribution method|lang=en|trad=Group contribution method|texte=''Group contribution method''}}, [[Elsevier]] Scientific, New York, 1979

=== Liens externes ===
* {{Lien|fr=PSRK|lang=en|trad=PSRK|texte=PSRK}} : modèle d'estimation pour les calculs des phases d'équilibre des mélanges de composants chimiques. Appelé ''Predictive Soave-Redlich-Kwong''
*  {{Lien|fr=MOSCED|lang=en|trad=MOSCED|texte=MOSCED}} : modèle pour estimer la limitation des coefficients d'activité à dilution infinie.
* {{en}} [http://www.aim.env.uea.ac.uk/aim/info/UNIFACgroups.html UNIFAC groupes structurels et ses paramètres]
* {{en}} [http://www.aiomfac.caltech.edu/index.html AIOMFAC modèle en ligne] UNIFAC, pour le calcul des coefficients d'activité dans des mélanges organiques et inorganiques.


{{portail|chimie}}

{{DEFAULTSORT:UNIFAC}}
[[Catégorie:Thermodynamique]]
[[Catégorie:Thermodynamique chimique]]
[[Catégorie:Chimie générale]]
[[Catégorie:Thermochimie]]