About the AIOMFAC model
A brief introduction
What is AIOMFAC?
The abbreviation AIOMFAC stands for "Aerosol Inorganic-Organic Mixtures Functional groups Activity Coefficients". As this name suggests, AIOMFAC is a thermodynamic model designed for the calculation of activity coefficients of different chemical species in inorganic-organic mixtures – with a focus on those inorganic and organic species known or suggested to be present in atmospheric aerosol particles. The AIOMFAC model is a general thermodynamic activity coefficient model also applicable to other scientific fields.
In the context of chemical thermodynamics activity coefficients are temperature-dependent, component- and mixture-specific parameters describing non-ideal mixing behavior in solutions. Interactions between organic and inorganic mixture components almost always lead to significant deviations from ideal mixing (i.e., deviations from Raoult's law).
The model version available on this website is termed AIOMFAC-web, emphasizing the web-interface and restrictions of this implementation compared to the more general AIOMFAC framework. Specific help and information regarding the use of AIOMFAC-web are given on the Hints & Examples page.
The AIOMFAC model framework and an initial model parameterization for aqueous electrolyte solutions and interactions between inorganic ions and alcohols/polyols at room temperature have been developed by Andreas Zuend (Zünd), Claudia Marcolli, Beiping Luo, and Thomas Peter at ETH Zurich, Switzerland (Zuend et al., 2008). An extension of AIOMFAC to many more organic functional groups and a comprehensive re-parameterization of interaction parameters between inorganic ions and organic functional groups has been realized in subsequent years as part of a collaboration between scientists from the California Institute of Technology (Andreas Zuend, John H. Seinfeld), ETH Zurich (Claudia Marcolli, Thomas Peter, Daniel M. Lienhard, Vacharaporn Soonsin, Ulrich K. Krieger), and the University of Manchester, UK (Alastair M. Booth, David O. Topping, Gordon McFiggans), published in 2011 (Zuend et al., 2011).
The AIOMFAC model has been developed and tested against hundreds of experimental datasets and the model has been shown to be useful, reasonably accurate, and practical for the calculation of activity coefficients in complex multicomponent mixtures as well as for thermodynamic equilibrium computations based on AIOMFAC (vapor-liquid equilibria, solid-liquid equilibria, liquid-liquid equilibria, deliquescence relative humidities, gas-particle partitioning) (Zuend et al. 2008; Zuend et al. 2010; Ciobanu et al. 2010; Zuend et al. 2011; Song et al., 2012; Zuend and Seinfeld, 2012; 2013).
In comparison to other thermodynamic activity coefficient models, AIOMFAC offers several features that most other models lack:
The AIOMFAC model details are given in the peer-reviewed articles describing the model and its parameterization (Zuend et al., 2008; Zuend et al. 2011). In these papers, a clear description of the thermodynamic basis is presented, including the definition of reference states for molality-based and mole fraction-based activity coefficients and their derivation from expressions of Gibbs excess energy. We encourage people interested in the mathematical expressions and physicochemical basis of AIOMFAC to read the papers. A basic introduction to the thermodynamic principles defining and governing expressions for thermodynamic equilibrium, chemical potentials, and activity coefficients, as well as related consequences such as the Gibbs-Duhem relation and Gibbs' phase rule can be found in advanced textbooks on chemical thermodynamics and are also described in Zuend (2007) and applied in Zuend et al. (2010). An overview of the experimental database used for the parameterization of AIOMFAC is given in the Supplement to the article of Zuend et al. (2011), available here.
Gibbs energy contributions: combination of different types of interactions
The excess Gibbs energy Gex characterizes for the overall non-ideality of a thermodynamic system. Different types of interactions involving charged and non-charged (neutral) solution species, interacting over different geometric ranges, can be identified. The total excess Gibbs energy of a multicomponent mixture in the AIOMFAC model is expressed in terms of three contributions from long-range (LR), middle-range (MR) and short-range (SR) interactions between molecules (or functional groups) and ions (Zuend et al., 2008), as shown in Fig. 1. The long-range part, an extended Debye-Hückel expression, and the semi-empirical middle-range part form together a Pitzer-like (group-contribution) model, enabling accurate descriptions of electrolyte solutions, from dilute to highly concentrated (potentially supersaturated) conditions. Short-range interactions between functional groups are calculated with the widely-used UNIFAC model (Fredenslund et al., 1975; Hansen et al. 1991) as modified by Zuend et al. (2011) to include improved specific interaction parameters from Marcolli and Peter (2005) and from Peng et al. (2001).
Since atmospheric aerosols likely contain tens, hundreds, or even thousands of different organic compounds, a practical way to consider this variety of species with a manageable degree of model complexity is by means of organic compound classes and functional groups instead of individual compounds. Therefore, AIOMFAC relies on a group-contribution concept that was introduced with the well-known UNIFAC model (Fredenslund et al. 1975; Hansen et al. 1991), which is also the basis of the LIFAC model (Yan et al. 1999; Kiepe et al. 2006). The LIFAC model and different versions of the UNIFAC model served as a basis for the development of AIOMFAC. In fact, AIOMFAC includes a modified UNIFAC model in addition to a Pitzer-like ion-interaction model and an organic-inorganic mixing term similar to the one in LIFAC, as shown in Fig. 1 above.
The group-contribution concept is based on the idea that physicochemical properties of organic molecules can be represented reasonably well by segmenting a molecule into different "functional" groups and consideration of specific properties of these functional groups. This leads to the description of a liquid mixture as a "solution of groups" (instead of a solution of molecules, see also Fig. 3 below). This concept has been proven useful; it is often a good approximation for the real behavior of organic mixtures and aqueous organic mixtures (see references to UNIFAC and mod. UNIFAC models). However, there are also limitations to a group-contribution approach as discussed in Zuend et al. (2011), e.g., when applied to complex, multifunctional organic compounds, effects of intramolecular interactions between functional groups in close proximity are usually not considered. Such intramolecular effects can be considered only by introducing additional types of functional groups or second order correction coefficients, which is often not a viable approach as it would further complicate the approach instead of keeping it as simple as possible.
Main groups and subgroups
AIOMFAC and UNIFAC distinguish between two levels of functional group classification: main groups and subgroups. Main groups comprise a less resolved classification of functional groups ignoring the number of hydrogen atoms bonded to carbon atoms, e.g., CHn (alkyl main group), or CHnCO (ketone main group). Subgroups further differentiate between groups of the same main group consisting of different numbers of H-atoms, e.g., between CH3, CH2, CH, and C of the same (standard) alkyl main group or between the subgroups CH3CO and CH2CO of the ketone main group. The AIOMFAC main groups and examples of their assignments to organic structures are shown in Fig. 2.
Main groups of organic compounds are used to describe middle-range interactions between dissolved inorganic ions and functional groups, as the number of H-atoms bonded to C-atoms plays a minor role compared to the polarity of carbon-oxygen bonds. That is, organic compounds are still treated as a combination of subgroups, but the activity coefficient contributions from different subgroups of a certain main group interacting with ions are assigned the same value (the main group interaction value). However, in the UNIFAC short-range part of AIOMFAC, interactions between organic functional groups and their interactions with water molecules (main group = subgroup = H2O) are resolved on the subgroup level. Therefore, in the input to AIOMFAC, the organic compounds need to be defined in terms of subgroups (main groups are then assigned automatically based on the stated subgroups). An example of the segmentation of an organic substance (ferulic acid) in terms of subgroups and their mixing with water and a salt (here NaCl) as a solution of groups is shown in Fig. 3.
The procedure for the definition of mixtures other than the one shown in Fig. 3, including systems containing multiple organic and inorganic substances, follows the same basic steps. While the number of different main groups and subgroups remain limited and manageable, the number of organic species which can be represented by AIOMFAC subgroups is, in principle, unlimited. Hence, adding hundreds of different species to a mixture for a specific AIOMFAC calculation only seems to make the system more complex, yet most of the organic components will share some (or all) of their subgroups with the pool of subgroups already present due to the definition of other system components. Adding compounds consisting of subgroups already present in the mixture will therefore change only the relative abundance (mole fractions) of subgroups, but not the complexity of the system as seen by AIOMFAC. Thus, as a limiting case, the solution-of-groups representation cannot become more complex than the maximum number of different subgroups available in AIOMFAC (currently ~40 organic subgroups (13 main groups), 12 ions, and water). This reduced complexity is one of the aspects that make the group-contribution concept a powerful and practical method. Another advantage of the group-contribution concept is its predictive capability: it can be used to predict properties of molecular structures that have not been measured or modeled before. This is important with regard to predictions of activity coefficients and pure compound vapor pressures of organic compounds in the atmosphere.
As described above, the list of chemical species available for calculations with AIOMFAC is enormous. Possible organic compounds are only limited by the types of organic functional groups currently available in AIOMFAC. AIOMFAC essentially comprises the whole list of organic functional groups available in standard UNIFAC (see Hansen et al., 1991), but since the focus of AIOMFAC is on organic-inorganic mixtures, only organic main groups for which also middle-range interaction parameters with inorganic ions have been determined are made available in AIOMFAC-web (with a few exceptions). Figure 4 shows the list of cations, anions, and organic main groups currently available. If a main group is available, all pertaining subgroups are available too – and these can be selected from the drop-down lists on the model input page (Run Model). Missing interaction parameters are indicated by a question mark in Fig. 4. Missing parameters could be the cause of an error message from AIOMFAC-web in case a mixture that would require undetermined parameters had been defined as model input. The reason for missing interaction parameters is the lack of experimental data required to determine corresponding AIOMFAC model parameters. In some cases, missing model parameters can be estimated based on an analogy approach as described in Zuend and Seinfeld (2012) using the functional group polarity series of Zuend et al. (2011) as constraints.
The valid temperature range covered by AIOMFAC depends on the type of a mixture activity coefficients are computed for. In case of aqueous electrolyte solutions and (aqueous) mixtures of organic components and electrolytes, the temperature range is limited to room temperature, in a strict sense to 298.15 K (25 °C). However, activity coefficients are rather weak functions of temperature (Gibbs-Helmholtz equation), e.g., in contrast to vapor pressures. Therefore, AIOMFAC calculations in a temperature range close to 298 K, say ±10 K, will be applicable. In many cases a larger deviation (±20 K, ±30 K) from 298 K will still provide a good estimate for activity coefficients in organic-inorganic mixtures, as the AIOMFAC model has been parameterized using also experimental data deviating from 298 K and the model contains several terms which at least to some extent account for a temperature effect (likely better constrained towards temperatures higher than 298 K). The rather weak temperature dependence of activity coefficients should not be confused with the temperature dependence of solubility limits of organic compounds and inorganic salts or the temperature dependence of dissociation constants, because equilibrium "constants" have also a temperature dependence on their own.
The temperature dependence of electrolyte-free mixtures (aqueous organic solutions) is described over a wider range of temperatures, since in this case AIOMFAC reduces to a modified UNIFAC model, including the relatively simple temperature dependent main group interaction parameterization of standard UNIFAC. In this case, the valid / applicable temperature range is ~280 K to ~400 K. Note that there are modified versions like UNIFAC Dortmund, which allow a physically better and more reliable description of temperature dependency. In all cases AIOMFAC will likely perform better at temperatures higher than room temperature as compared to similar deviations to lower temperatures, because the experimental database used to determine interaction parameters was biased towards more data at higher temperatures. An improved description of the temperature dependence of activity coefficients will be made available in a future version of AIOMFAC.
References↑ back to text
|© 2012–2022 Andreas Zuend, California Institute of Technology | McGill University last page update: 2022-05-06|
|Home | Disclaimer | Citation | About | Credits | Contact|