Noora
Hyttinen
Department of Chemistry, Nanoscience Center, University of Jyväskylä, FI-40014 Jyväskylä, Finland. E-mail: noora.x.hyttinen@jyu.fi
First published on 5th April 2023
Liquid–liquid phase separation (LLPS) affects the water uptake of aerosol particles in the atmosphere through Kelvin and Raoult effects. This study investigates LLPS in ternary mixtures containing water and two organic compounds, using a conductor-like screening model for real solvents (COSMO-RS). COSMO-RS found LLPS in all of the studied mixtures containing water and proxies for primary and secondary organic aerosol (POA and SOA, respectively), due to the limited solubility of the hydrophobic POA proxies in water. The computations predict additional three-phase states in some of the SOA–POA–water mixtures at relative humidity (RH) close to 100%, which was not observed in experiments, likely due to the relatively low RH (90%) used in the experiments. A computational method, such as COSMO-RS, allows for the estimation of new information on mixing states and mixtures that cannot be accessed experimentally. Comparison with experiments can also provide insight into which types of compounds may be present in SOA. Additionally, the possibility of LLPS can be assessed faster with rough estimates rather than by computing the whole phase diagram.
Water activity is generally lower in concentrated solutions than in dilute solutions, due to differences in water concentrations. However, phase separated organic–water mixtures strongly deviate from ideality. Water activity can even be close to unity in the equilibrium of binary two-phase states, while in one-phase mixed states with equal water content, water activity can be significantly lower.4 If the particles contain multiple organic species with different polarities, phase separation can occur at relatively low relative humidities (water activities).5
LLPS in aerosol particles is highly dependent on the size of the particle.6–11 Experiments have shown that smaller aerosol particles (specific size depending on the mixture) are homogeneous, while larger particles with the same composition have phase separated.12 This indicates that a higher relative humidity (RH) is needed for smaller particles to phase separate compared to larger particles or bulk solutions, in cases where a high RH leads to phase separation.
In recent years, LLPS has been observed in particles containing one or more different organic compounds.5,13–20 For example, Mahrt et al. investigated particles containing proxies of primary and secondary organic aerosol (POA and SOA, respectively) under varying RH conditions. They found that the particles phase separate at all RH if the oxygen-to-carbon ratio of the POA proxy is low enough (O:
C < 0.11).5 Additionally, Mahrt et al. found that SOA–POA particles phase separate if the difference in the oxygen-to-carbon ratios (O
:
C) of the POA (single compound) and SOA (mixture of volatile organic compound oxidation products) is higher than 0.265.20 Recently, three liquid phases were observed in aerosol particles consisting of water, inorganic salt and two organic species with different polarities.21
Experimentally, it is not possible to investigate the phase separation behavior of real atmospheric aerosols with known compositions. Models are therefore necessary for estimating the phase behavior of atmospheric aerosols. Thermodynamic models, such as conductor-like screening model for real solvents (COSMO-RS22–24) and Aerosol Inorganic–Organic Mixtures Functional groups Activity Coefficients (AIOMFAC25–27), can be used to estimate the equilibrium liquid-phase compositions of mixtures. Both models are used for bulk liquid solutions, which makes them optimal for modeling LLPS in larger aerosol particles. COSMO-RS, implemented in the COSMOtherm program,28 has been shown to accurately estimate condensed-phase properties, such as solubilities and activity coefficients, of atmospherically relevant multifunctional compounds and water, if suitable conformer distributions are used in the calculation.29–31
Here, COSMOtherm was used to calculate the phase diagrams of some of the ternary SOA–POA–water mixtures, investigated experimentally by Mahrt et al.5 Additionally, the phase separation behavior was predicted for a large set of SOA–POA–water mixtures with possible SOA constituents (psat < 10−5 Pa), as well as POA proxies with relatively high O:
C ratios (O
:
C ≥ 0.125). The SOA constituents were selected from previous studies that have identified possible oxidation products of α-pinene, Δ3-carene, cyclohexene and isoprene.4,32–34
![]() | (1) |
![]() | (2) |
The pure compound (xi = 1) is used as the reference state composition x°. The use of pseudo-chemical potential allows for the calculation of activity coefficient in infinitely diluted states (xi → 0).
Phase equilibrium calculations are based on activity coefficients. At LLE, the LLE condition applies for all compounds i in the mixture:
aαi(xa) = aβi(xβ) | (3) |
In COSMOtherm LLE calculations, the accuracy threshold is by default 10−5. The number of iteration points in the LLE calculations was increased to 100000, because the LLE was not found for all mixtures in the default 1000 iteration points.
All conformers were found using the systematic conformer sampling and Merck molecular force field (MMFF9438) of Spartan20.39 For compounds that have more than 200 conformers, a subset of conformers were selected for further calculations based on their pseudo-chemical potentials in infinite dilution in water and water-insoluble organic matter (WIOM, CC(O)C1OC2OC(OC2(C)O1)C1
CC(C)
CC(C)
C140). The pseudo-chemical potentials were predicted using a machine learning model.37,41 200 lowest pseudo-chemical potential conformers in each of the two solvents (at most 400 conformers) were selected for density functional theory calculations using the COSMOconf program.42,43
The COSMOconf calculations included single-point calculations at the BP/SV(P) level of theory, followed by removal of conformers with similar pseudo-chemical potentials in a set of commonly used solvents. The geometries of all conformers were then optimized at the BP/SV(P) and BP/def-TZVP levels of theory with duplicate removal after both optimization steps based on both geometries and pseudo-chemical potentials. The input files for COSMOtherm calculations were obtained through BP/def2-TZVPD-FINE level of theory single-point calculations. The input files of cyclohexene and isoprene oxidation products were taken directly from Hyttinen et al.44,45
In order to find the 10 most stable condensed-phase conformers (5 in both water and WIOM) for the thermodynamic property calculations, the pseudo-chemical potential of each conformer was calculated using the BP_TZVPD_FINE_21 parametrization of COSMOtherm. Additionally, only conformers within 2 kcal mol−1 from the lowest pseudo-chemical potential were selected. This ensures that unfavorable conformers are not included when the compound has very few conformers, as COSMOtherm tends to favor low COSMO-energy conformers even if their pseudo-chemical potential in the solution is high.30
Mahrt et al.5 reported LLPS in ternary mixtures with known organic:
organic mass ratios (1
:
1) under three different RH conditions: dry (<5%), 45% and 90%. There were three distinct cases in the ternary SOA–POA–water mixtures:
1. Phase separation occurring at all RH, e.g., glycerol–DES.
2. Phase separation occurring at relatively low RH (<90%), e.g., maleic acid–DES and cis-pinonic acid–DOP.
3. No LLPS at RH below 90%, e.g., cis-pinonic acid–DES.
Here, COSMOtherm calculations are compared with the phase separation behavior of the three example cases using the COSMOtherm-computed water activities (and LLE) in mixing states with varying water contents.
For the first case, Mahrt et al.5 observed LLPS in glycerol–DES particles in both dry and humid conditions. COSMOtherm calculations predict phase-separation in all water contents, see Fig. 1. All of the experimental mixtures (1:
1 organic mass ratio, dashed black line) fall into the unstable region, indicating LLPS into one phase containing water and glycerol, and a DES-rich phase that contains only small amounts of water and glycerol.
![]() | ||
Fig. 1 COSMOtherm-derived phase diagram of glycerol, diethyl sebacate and water at 296 K. The solid lines represent tie lines, where the color is water activity in both liquid phases (solid circles) at phase separation. The dashed line represents the 1![]() ![]() |
For the second case, Fig. 2 shows the calculated ternary phase diagram of the maleic acid–DES–water mixture. The calculations agree with the experimental observations of one phase at 45% RH and two phases at 90% RH. COSMOtherm-predicted phase separation RH is around 75% for the experimental 1:
1 mass ratio of the two organics. An additional crystal phase was seen in the experiments under dry conditions (RH <5%), which is likely caused by maleic acid being a solid at 296 K. However, the solid solubility of maleic acid was not investigated further in this study.
![]() | ||
Fig. 2 COSMOtherm-derived phase diagram of the ternary maleic acid–diethyl sebacate–water mixture at 296 K. Solid lines represent tie lines. The dashed contour lines represent water activity in a one-phase state and are derived from individual points calculated in different mixing ratios. The dashed black line represents the 1![]() ![]() |
For the third case, Mahrt et al.5 observed no phase separation in the cis-pinonic acid–DES particles. On the other hand, LLPS was observed at high RH (90%) in the cis-pinonic acid–DOP particles,5 corresponding to the second case. Fig. 3(a) and (b) show COSMOtherm-derived liquid-phase diagrams of the ternary aqueous mixtures containing DES and DOP, respectively.
![]() | ||
Fig. 3 COSMOtherm-derived phase diagram of cis-pinonic acid, water and (a) diethyl sebacate and (b) dioctyl phthalate at 296 K. The dashed black line represents the 1![]() ![]() |
In the cis-pinonic acid–DES–water mixture, an almost pure aqueous phase is formed at close to 99% RH (Fig. 3(a)). The cis-pinonic acid–DES particle would therefore not take up enough water in the experimental RH of 90% to result in LLPS. This is in agreement with the experiments.5
In the cis-pinonic acid–DOP–water mixture, COSMOtherm predicts an additional phase separation into two phases, both containing significant amounts of all three components. The water activity at the tie line with the lowest water mole fraction found by COSMOtherm is 0.41. This is in agreement with the experiments,5 where two phases were observed at 90% RH in particles containing DOP but not DES. Even though Mahrt et al.5 did not observe LLPS at 45% RH, the small size of the particles and uncertainties of COSMOtherm calculations may explain why a higher RH was needed for LLPS in the experiments than what is predicted by COSMOtherm.
COSMOtherm found additional LLE (dotted lines in Fig. 3(b)), where the binodal compositions fulfill the LLE conditions, but the binodal compositions are inside the other phase separation region (empty markers in Fig. 3(b)). These compositions are likely in the metastable region, which may become relevant in smaller particles. Additionally, the solid lines form a triangle, marking a 3-phase region. This was not seen in the experiments with up to 90% RH. A higher RH (close to 100%) would likely be required to reach high enough water content in the particles to form an additional water-rich phase predicted by COSMOtherm.
With the exception of the cis-pinonic acid–DOP–water mixture, there is a clear correlation between phase separation in the three binary mixtures of the ternary mixture and the RH required for phase separation in the ternary mixtures. For example, when LLE is found between the two organic compounds (in addition to one of the organic compounds and water), the mixture phase separates at all RH (Fig. 1). On the other hand, if both organic phases separate with water, water activity is likely to be very high at phase separation and RH close to 100% is needed for LLPS (Fig. 3(a)). If one of the organic compounds is miscible with both water and the other organic, the LLPS RH is more likely to be below 95% (Fig. 2).
Fig. 4 shows how the infinite dilution activity coefficients of the studied SOA constituents correlate with the phase separation behavior of the ternary mixtures. The cases similar to the cis-pinonic acid–DOP–water mixture (i.e., the existence of a second phase separation region) were identified by a single LLE calculation at SOA:
POA:water mixing ratio of 0.3
:
0.4
:
0.3.
The green markers represent mixtures that phase separate at all RH. In general, the binary SOA–POA mixtures phase separate and most of the water will partition to the SOA-rich phase. On the other hand, the mixtures that require a high RH to phase separate (blue markers in Fig. 4) contain SOA that phase separates with water. In these solutions, water forms its own phase, while the POA and SOA form one organic-rich phase that also contains water. The intermediate cases (magenta markers in Fig. 4) that phase separate at RH between 1 and 98% contain SOA that does not phase separate with water or POA. Both phases will contain significant amounts of all 3 components. It should be noted that the SOA:
POA mixing ratio affects the LLPS behavior of the particle. For example, “LLPS at all RH” means that there exists some SOA
:
POA mixing ratio region, in which the binary SOA–POA mixture phase separates into two liquid phases.
Table S1 of the ESI† summarizes the types of LLPS behaviors of all of the studied ternary mixtures. Several mixtures that are predicted to form three liquid phases at RH close to 100% were also identified. Three phases can form, when all of the three constituents are partially miscible with the other two constituents. Additionally, there exist mixtures similar to cis-pinonic acid–DOP–water, where three phases were predicted, even though full miscibility was predicted between cis-pinonic acid and DOP (Fig. 3(b)).
Octanol, WIOM and DES (Fig. 4(a)) present similar phase behavior in relation to the SOA activity coefficients in infinite dilution in water and POA, while the mixtures containing DOP are likely to phase separate with relatively low γPOASOA values (Fig. 4(b)). Chemically DES and DOP are quite similar compounds, as they have identical oxygen containing functional groups. However, DOP contains a benzene ring and longer alkyl groups than DES. Still, COSMOtherm predicts LLPS in the SOA–DOP mixture with much lower SOA activity coefficients (γDOPSOA ∼ 1) compared to the corresponding SOA–DES mixtures (γDESSOA ∼ 10). These results explain why two phases were observed in more mixtures containing DOP than other POA proxies with similar O:
C ratios.5 The main difference between DOP and the other studied POA proxies is the higher tendency of DOP to phase separate even when the infinite dilution activity coefficients are relatively low (high solubility).
Generally, if LLPS exists in the binary SOA–water and POA–water mixtures, the ternary mixture phase separates at RH close to 100%. Additional phase separations are mainly predicted for DOP as the POA proxy as well as some mixtures containing WIOM and DES, when the SOA constituent is large (molar mass > 290 g mol−1). Three-phase mixtures are predicted for all four POA proxies, when all three binary pairs of the ternary mixture phase separate.
Fig. 5 shows phase diagrams of caronic acid–POA–water mixtures with the four studied POA proxies. Caronic acid is an example SOA constituent, that exhibits the same LLPS behavior (0% < LLPS RH < 99%) with the four POA proxies. However, the predicted phase separation RH is much higher for mixtures containing DES, octanol and WIOM (above 90%; Fig. 5(a), (c) and (d), respectively) than DOP (RH < 35%; Fig. 5(b)). The lower LLPS RH of the caronic acid–DOP–water mixture is caused by caronic acid partitioning more to the aqueous phase than the DOP-rich phase, especially at high caronic acid mole fractions. Additionally, the aqueous phase in the caronic acid–WIOM–water mixture contains more caronic acid than the equivalent aqueous phases in the mixtures containing DES and octanol. Unfortunately, the differences in LLPS RH cannot be easily determined with binary LLE calculations, and the full ternary phase diagrams need to be computed to obtain this information.
For example, previous COSMOtherm calculations showed that the interaction between water and the different isomers (trans-β, cis-β, δ1, δ4) of isoprene epoxydiol (IEPOX; C5H10O3) is similar.4 However, the addition of a POA proxy (WIOM or DES) to the mixture leads to different phase behavior between the different IEPOX isomers. Cis- and trans-β-IEPOX have limited solubilities in WIOM and DES, while δ1- and δ4-IEPOX are miscible with both WIOM and DES. Fig. 6(a) and (b) show phase diagrams of WIOM, water and δ1- and cis-β-IEPOX, respectively. The LLPS RH of the δ1-IEPOX–WIOM–water mixture is around 30%, while cis-β-IEPOX–WIOM–water particles phase separate at all RH.
![]() | ||
Fig. 6 COSMOtherm-derived phase diagrams of (a) δ1-IEPOX and (b) cis-β-IEPOX, WIOM and water at 296 K. The lines represent tie lines (two-phase states) of the mixtures. |
Compared to the glycerol–DES–water mixture, a significant amount of WIOM mixes in the aqueous SOA-rich phase of the IEPOX–WIOM–water mixtures. Additionally, the WIOM-rich phase in the δ1-IEPOX–WIOM–water mixture (Fig. 6(a)) contains 0.07–0.12 mole fraction of water. In this case, adding δ1-IEPOX to the mixture enhances the solubility of water in the WIOM-rich phase compared to the binary WIOM–water mixture. Cis-β-IEPOX has the opposite effect, worsening the solubility of water in the WIOM-rich phase.
In general, compounds containing more than 11 carbon atoms are more likely to phase separate with either water or POA. This indicates that large SOA constituents form particles that have LLPS behavior independent of the RH (see Table S1 of the ESI†). Instead, the large SOA constituents form one-phase mixtures under RH < 100% more often than the smaller SOA constituents.
Mahrt et al.47 found that a SOA–SOA particle phase separates, when the difference in the average O:
C ratios, Δ(O
:
C), of the two SOA is equal to or higher than 0.47, while two SOA with similar O
:
C (Δ(O
:
C) < 0.47) form one-phase particles. No RH dependence was observed in these experiments, corresponding to the glycerol–DES–water mixture. Two phases were observed in particles containing oxidation products of one sesquiterpene with a low O
:
C ratio (valencene, β-caryophyllene, farnesene) and one aromatic compound with a high O
:
C ratio (toluene, catechol). All of the combinations containing α-pinene ozonolysis products formed one-phase particles.47
The O:
C ratio of the studied SOA constituents ranges from 0.26 to 1.2. COSMOtherm predicts miscibility between most pairs of the studied SOA constituents, regardless of their O
:
C ratios. An exception is glycerol, which phase separates with some isomers of the larger studied SOA constituents, C12 and C18 compounds with O
:
C ratios between 0.33 and 1.17. The Δ(O
:
C) values of these phase separating binary mixtures range from 0 to 0.67. Fig. 7 shows the phase diagram of glycerol–C18H26O6–water mixture (Δ(O
:
C) = 0.67). Here, water is more soluble in the SOA constituent than any of the POA proxies. This leads to a smaller two-phase region compared to the glycerol–DES–water mixture (Fig. 1). The occurrence of LLPS is therefore more dependent on the SOA
:
SOA mixing ratio compared to SOA
:
POA mixing ratios.
![]() | ||
Fig. 7 COSMOtherm-derived phase diagram of glycerol, C18H26O6 (CC1(C)C(CC![]() ![]() ![]() ![]() |
As was hypothesized by Mahrt et al.47 the actual range of O:
C ratios of each SOA likely plays a role in the phase state of SOA–SOA particles. Unfortunately, only the average O
:
C ratios of SOA are available from experiments.47 Perhaps the SOA from some aromatic precursors (high average O
:
C) contains compounds that are similar to glycerol (and water). These compounds could separate from other highly oxygenated SOA constituents even under dry conditions and form phase separated particles under all RH conditions. Or possibly some oxidation products of sesquiterpenes are similar to the POA proxies (low O
:
C and hydrophobic), while having low enough volatilities to remain in the particle phase (high molar mass).
At LLPS, water activity may be below ideality (enhanced water uptake), if the SOA-rich phase contains only small amounts of POA, as is the case in the glycerol–DES–water mixture (Fig. 1). Even in the cis-β-IEPOX–WIOM–water mixture, the activity coefficients of water at phase separation (aqueous phase) are slightly above unity for all of the computed tie-lines. Based on calculations presented here, POA lowers water uptake into aerosol particles. Water uptake is only enhanced (in terms of decreased water activity) by LLPS in SOA–POA aerosol particles, when the SOA-rich phase contains very little POA (xPOA < 0.001).
Aerosol process models that consider LLPS through lowered surface tension of a hydrophobic organic shell phase and the non-ideality of water have been able to explain experimental hygroscopic growth better than models that assume one-phase particles.3 Additionally, when the water content of the particle is high enough at equilibrium under high RH conditions, the particle becomes mixed, leading to a surface tension equal to that of pure water.3 In these aqueous particles, the enhancing effect of the lower surface tension on water uptake is lost. However, some of the low surface tension organics may still act as surfactants and partition to the particle–air interface of one-phase particles.
Based on comparisons between COSMOtherm calculations and experiments, COSMOtherm is able to predict experimental phase equilibria in binary and ternary supermicron particles. The calculated results complement experimental data on LLPS to improve the description of different SOA and POA types in atmospheric models. In order to further validate the accuracy of COSMOtherm in LLPS calculations, more measurements are needed under different RH conditions using a range of organic mixing ratios.
Footnote |
† Electronic supplementary information (ESI) available: Comparing COSMOtherm calculations with experiments; and table detailing the LLPS behavior of all studied mixtures. See DOI: https://doi.org/10.1039/d3cp00691c |
This journal is © the Owner Societies 2023 |