Open Access Article
Benjamin L.
Deming
ab and
Paul J.
Ziemann
*ab
aDepartment of Chemistry, University of Colorado, Boulder, Colorado 80309, USA. E-mail: paul.ziemann@colorado.edu; Fax: +1-303-492-1149
bCooperative Institute for Research in Environmental Sciences (CIRES), Boulder, Colorado 80309, USA
First published on 20th January 2021
Partitioning of nitric acid and sulfuric acid between aqueous and organic phases may play a role in the chemistry of aqueous/organic phase-separated atmospheric aerosol particles, as well as other environmental systems. To assess the extent to which these acids partition to organic phases, organic compounds with O
:
C ratios ranging from 0 to 0.75 were mixed with acidified aqueous phases. The organic phase was removed and its acid content quantified with a mass-based measurement. Up to 60% of the acid was observed to partition to the organic phase, present either entirely in the dissociated form, the undissociated form, or both depending on the organic phase. Values of partitioning coefficients, Ddiss and Dun, corresponding to partitioning of the dissociated and undissociated acid, were determined by fitting a simple partitioning model to the data. Values of Dun correlated well with the O
:
C ratio of the organic component and to a lesser extent with the water content of the organic phase. The addition of ammonium sulfate led to an increase or decrease in partitioning of sulfuric acid to the organic phase, apparently depending on the balance between the resulting enhancement in the HSO4− concentration and the reduction in the water content of the organic phase. These results demonstrate that significant fractions of sulfuric and nitric acid can partition to the organic phase in phase-separated systems in undissociated and dissociated forms, and thus potentially participate in specific or general acid catalysis of organic reactions.
Environmental significanceAqueous/organic phase-separated systems are common to the atmosphere in the form of aerosols, but are also found in other environments such as soils and the surface microlayer of oceans and lakes. In many of these systems inorganic acids are presumed to reside predominantly in the aqueous phase, but data on this partitioning are relatively rare. Measurements of the partitioning of nitric acid and sulfuric acid between an aqueous phase and a number of model organic phases covering a range of polarity demonstrate that significant quantities of acid can partition to the organic phase, to an extent that depends on the O : C ratio of the organic phase, as well as other factors. These acids can then participate in organic reactions.
|
:
C ratio).13 This value serves as a surrogate for polarity or oxidation state, as it is readily determined by mass spectrometry in laboratory and field measurements.14
In addition to having important physical properties, atmospheric aerosols also serve as sites for condensed-phase chemical reactions, essentially acting as suspended micro-reactors. Particle-phase reactions of organic compounds can be initiated by OH radicals, NO3 radicals, and O3, similar to atmospheric gas-phase oxidation, and they can occur via the coupling of functional groups (such as carbonyl, hydroxyl, hydroperoxy, and carboxyl) to form oligomers.15,16 These oligomerization reactions can modify particle chemical and physical properties and also contribute to particle growth by reducing effective compound vapor pressures, thus drawing more mass into the particle through subsequent gas-to-particle partitioning. Like many condensed-phase reactions, a number of those that can occur in aerosol particles are catalyzed by strong acids. For example, the dehydration of alcohols, the formation of hemiacetals and hemiketals, the formation of acetals and ketals, the conversion of hydroperoxides to carbonyls, and ester formation can be significantly faster in the presence of strong acid.16 Notably, many of these reactions produce water as a product, so that in an aqueous environment there will be a strong driving force against their occurrence. But in a phase-separated system, one in which organic reactants are concentrated in an organic phase where the water content is much reduced, such reactions may be highly favorable as long as sufficient acid partitions to the organic phase to drive catalysis.
The most common inorganic acids in the atmosphere are sulfuric acid, predominantly formed from oxidation of SO2, and nitric acid, formed from various reactions of nitrogen oxides.17 Though strong acids in aqueous solution, they may exist in either their dissociated or undissociated form in predominantly organic phases. This distinction is especially important for the catalysis of particle-phase reactions. When the solvent mixture is sufficiently polar to allow significant dissociation, the rate is likely controlled by the H+ concentration (referred to as specific acid catalysis). On the other hand, if the undissociated form of the acid dominates in the solution, the rate can instead be determined by the molecular form (general acid catalysis).18 Depending on the reaction and conditions, one pathway may be significantly faster than the other. It is therefore important to have some understanding of the extent to which sulfuric and nitric acid partition to organic phases, as well as the forms in which they reside.
The purpose of this study was to investigate the partitioning of sulfuric and nitric acid in model organic/aqueous phase-separated systems that were chosen to act as proxies for phase-separated aerosol. Understanding the partitioning behavior of these acids in these simple model systems can help determine the extent to which the processes discussed above are relevant to real environmental systems.
:
C ratios and therefore polarities. To be compatible with the partitioning measurements they also needed to be non-reactive towards both strong acid and ammonia for the concentrations and exposure times used here. Hexane, 2-hexanone, 2-butanone, ethyl acetate, and propylene carbonate were thus used as organic phases. Other relevant organic compounds tested, notably carboxylic acids, alcohols, and aldehydes, were found to react with the ammonia and/or strong acid, and therefore could not be used. The phases and relevant properties are listed in Table 1, and the structures of the organic components are given in ESI.†
| Organic or aqueous phase | O : C ratio |
Molecular weight (g mol−1) | Organic phase water (vol%) | Density (g mL−1) |
|---|---|---|---|---|
| Hexane | 0.00 | 86.2 | 0.01 | 0.66 |
| 2-Hexanone | 0.17 | 100.2 | 2 | 0.81 |
| 2-Butanone | 0.25 | 72.1 | 1 | 0.81 |
| Ethyl acetate | 0.50 | 88.1 | 3 | 0.90 |
| Propylene carbonate | 0.75 | 102.1 | 8 | 1.21 |
| Water | — | 18.0 | — | 1.00 |
| Aqueous ammonium sulfate (30%, w/w) | — | — | — | 1.17 |
The aqueous phases were acidified to create solutions between approximately 0 and 1 M, not considering any reaction with ammonium sulfate. Aliquots (typically 1.5 mL each) of organic phase and acidified aqueous phase were added to test tubes and then mixed by vortexing for approximately 30 s, a time that is shown in the ESI† to be sufficient to establish partitioning equilibrium. The test tubes were then centrifuged and 1 mL of organic phase was transferred to a vial pre-weighed on a Mettler Toledo XS3DU Microbalance. Because the addition of acid or ammonium sulfate could change the overall volume of each phase, the total aqueous and organic volumes after mixing were measured using a syringe that is accurate to approximately 0.05 mL. The sampled organic solution was exposed to a slow stream of ammonia (∼600 cm3 min−1) for ∼15 s, which is shown in ESI† to be sufficient to convert any nitric or sulfuric acid to the corresponding ammonium salt. This was often accompanied by the visible crystallization of these salts. Since strong acids were used, experiments were performed wearing gloves, a lab coat, and goggles. Samples were exposed to ammonia in a fume hood and remained within one for at least 10 min to allow for any residual vapor to dissipate.
With the exception of propylene carbonate, the organic phases were all sufficiently volatile to be removed via evaporation in a hood under a steady stream of nitrogen. Since ammonium nitrate might also evaporate during this process (it is much more volatile than ammonium sulfate), experiments described in ESI† were conducted to verify that under the conditions used here no ammonium nitrate was lost by evaporation. The vials, then containing only ammonium salt, assumed to be either ammonium sulfate or ammonium nitrate, were again weighed on the microbalance. The balance is accurate to ±0.5 μg, corresponding to minimum detectable organic-phase concentrations of 6.3 μM for nitric acid and 3.8 μM for sulfuric acid. Blanks were measured for all systems by conducting identical experiments but without adding acid, and in all cases yielded negligible mass after the organic phase was dried. Samples containing propylene carbonate were dried to a constant mass with nitrogen in a hood in an attempt to remove any organic-phase water, and then weighed. The mass of propylene carbonate was corrected for by treating the samples containing no added acid as a blank. This approach worked well for nitric acid but resulted in nonphysical results (blank mass higher than sample mass) for sulfuric acid, apparently because some residual water could not be fully removed by drying when ammonium sulfate was present. The results for experiments with propylene carbonate have larger uncertainties due to these additional complications.
We note that although the assumption that the mass of salt measured in all these experiments corresponded to either ammonium nitrate (NH4NO3) or ammonium sulfate ((NH4)2SO4) is valid for the nitric acid experiments, since ammonium nitrate is the only possible salt, in the sulfuric acid experiments ammonium bisulfate (NH4HSO4), letovicite ((NH4)3H(SO4)2), or ammonium sulfate hydrate ((NH4)2SO4·H2O) might also have formed. If the salt was one of these three instead of ammonium sulfate, then the reported concentration of sulfuric acid would be in error by −15%, −7%, and 12%, respectively. Considering the other sources of uncertainty in these measurements, these alternate forms could be present without significantly affecting the results presented here. Nonetheless, in an attempt to verify that the salt formed was ammonium sulfate, the behavior of a phase-separated salt sample and ammonium sulfate and ammonium bisulfate standards were compared when heated in an oven at 230 °C. Although the thermal decomposition of ammonium sulfate salts is complicated,22 it is expected that ammonium sulfate would be stable at this temperature while ammonium bisulfate, letovicite, and ammonium sulfate hydrate would not. Based on the results of these experiments, which are described in ESI,† it appears that the phase-separated salt sample was ammonium sulfate.
The acidity of the organic phases in these systems was also measured over a larger concentration range using pH strips (Sigma Chemical Company). Aqueous phases containing between 10−6 and 10 M of either HNO3 or H2SO4 were prepared and then 2 mL of each aqueous phase and 2 mL of either hexane, 2-butanone, ethyl acetate, and propylene carbonate were added to a test tube and mixed by vortexing for 30 s. The organic layer was then removed using a syringe and a pH strip was submerged in the solvent. The four dye squares were compared visually to the provided key and matched to the closest pH value.
The dissociation of these acids in the aqueous phase is well understood and the corresponding values of the acid dissociation constants for the aqueous and organic phases, Ka,aq and Ka,org, are defined using the species concentrations as
![]() | (1) |
![]() | (2) |
Concentrations are designated in units of moles of species per liter of the noted phase (either organic or aqueous). The partitioning of undissociated and dissociated acid between the organic and aqueous phases is quantified by the distribution coefficients, Dun and Ddiss, defined as
![]() | (3) |
![]() | (4) |
Small distribution coefficients indicate that acid (in either the dissociated or undissociated form) prefers the aqueous phase over the organic phase. Eqn (1)–(4) and the equations Haq+ = Aaq− and Horg+ = Aorg− can also be combined to obtain
![]() | (5) |
| Torg = HAorg + Aorg− | (6) |
![]() | (7) |
. Lastly, we define the ratio of total acid in the organic phase to the total in the aqueous phase, R, as![]() | (8) |
Because nitric and sulfuric acid are strong acids, they are almost entirely dissociated in the aqueous phase. Using values of Ka,aq of 40 M for nitric acid and 104 M (a lower limit) for sulfuric acid,23 at most 3% and 0.01% of each should be in the undissociated form in the aqueous phase in these experiments. Assuming that Taq ≈ Aaq− = Haq+ and that
, eqn (8) becomes
![]() | (9) |
Substituting in the equations defining Dun and Ddiss, this becomes
![]() | (10) |
Plotting values of R as a function of Taq should therefore result in straight lines. The slope of the line corresponds to
and the y-intercept is equal to Ddiss.
It is worth considering some limiting cases of this relationship in order to orient one to the figures that follow. As noted above, the values of Ka,aq for nitric and sulfuric acid are 40 and >104, so that these acids are predominantly dissociated in water. If Dun is on the order of 1 or less, indicating that undissociated nitric and sulfuric favor the aqueous phase over the organic phase, then
. Eqn (10) will therefore reduce to R = Ddiss and the plots will appear as flat, straight lines with y-intercepts equal to Ddiss. In this case, the fraction of total acid in the organic phase is constant and does not depend on how much acid is present. Furthermore, since
![]() | (11) |
and a y-intercept ≈ 0. In cases when significant fractions of the acid in the organic phase are present in the undissociated and dissociated forms, then the plot will have a non-zero slope and y-intercept. We note here that the inclusion of the second dissociation of sulfuric acid has a slight effect on this framework, as shown by derivations in the ESI† that take this dissociation into account. Under the assumption that SO42− does not significantly partition to the organic phase (discussed in more detail in Section 3.4 and the ESI†), these considerations result in correction factors of 1.11 and 1.46 applied to the derived values of Dun and Ddiss, respectively.
AIOMFAC is available for download (https://github.com/andizuend/AIOMFAC), and here the liquid–liquid equilibrium model was built on top of it using Zuend and Seinfeld (2013)24 as a guide. To ensure that this version of the model performs as expected, Fig. 5d in Zuend and Seinfeld (2013)24 was captured using the “IgorThief” tool in Igor Pro and compared with the corresponding model outputs. This figure was chosen because it most closely resembles the kinds of systems and results the model was used for in this work. The comparison, shown in Fig. S5,† indicates that this version of the model adequately reproduces those results. In addition to using AIOMFAC for the partitioning modeling, it also was used to calculate the equilibrium concentrations of species in aqueous sulfuric acid and sulfuric acid/ammonium sulfate solutions.
:
C ratios of these organic components are 0.0, 0.17, 0.25, 0.50, and 0.75, respectively, indicating that the polarity increases in the order shown. This trend is also roughly reflected in the vol% water in these components of 0.01, 2, 1, 3, and 8, respectively. One can expect that the addition of water to the organic phase will increase the polarity of the resulting solution, thus enhancing the extent to which both the undissociated and dissociated forms of the acid can be solvated and increasing partitioning of both species to the organic phase. In the figures, values of R(Torg/Taq) are plotted against Taq, as suggested by the discussion above and eqn (10). Least-squares lines were first fit to the data, and then because of the significant scatter, in cases where the uncertainty in the slope was greater than the value the slope was set to zero. And in one case where the y-intercept was slightly negative, and thus nonphysical, that value was also set to zero. Visual inspection of the resulting lines indicates that this approach is reasonable, and it allows the profiles to be classified according to zero or positive slope and zero or positive y-intercept, thus simplifying the data interpretation and comparisons between systems.
The ratio Torg/Taq is approximately 0.001 for hexane; 0.01–0.1 for 2-hexanone, 2-butanone, and ethyl acetate; and 0.1–1 for propylene carbonate. These values correlate better with the vol% water for each group of 0.01, 1–3, and 8 than with the O
:
C ratios of 0.0, 0.17–0.5, and 0.75, although in either case there is clear evidence that partitioning increases with increasing polarity of the organic phase. The results also indicate that nitric acid partitions more strongly to the organic phase than sulfuric acid. As discussed in Section 2.3, the presence of three types of lines: (1) zero slope and positive y-intercept, (2) positive slope and zero y-intercept, and (3) positive slope and positive y-intercept are indicative of acids being fully dissociated, undissociated, and partially dissociated, respectively, in the organic phase. In addition, the y-intercept gives the value of the Ddiss, the distribution coefficient for NO3− or HSO4−, and the slope gives the value of
, the ratio of the distribution coefficient for HNO3 or H2SO4 divided by its acid dissociation constant in water. Values of Ddiss and Dun determined from the lines are given in Table 2, although for lines with zero slope the value of Dun cannot be determined. The uncertainties were taken from the fits to the data as output by Igor Pro. Values of Ddiss and Dun range from approximately 0.001–0.1 and 2–60, respectively, indicating that larger fractions of undissociated acid partition to the organic phase than dissociated acid. It should be remembered, however, that for values of Ka,aq of 40 M for nitric acid and 104 M for sulfuric acid (Levanov et al., 2018)23 at most about 3% and 0.01% of these acids should be undissociated in the aqueous phase. As a result, except in the case of ethyl acetate, a significant fraction of the acid present in the organic phase is always dissociated. When both Ddiss and Dun are known, Ka,org can be calculated using eqn (5). This is only possible with this data set for nitric acid in 2-hexanone and propylene carbonate, resulting in values of 0.11 M and 0.022 M, respectively. Unsurprisingly, the dissociation of nitric acid is therefore significantly reduced in predominantly organic phases as compared to the aqueous phase. Although the values of Ddiss and Dun indicate that in some systems the undissociated form of the acid dominates, in others the dissociated form does, and in still others both forms are significantly present in the organic phase, one could expect that for the complex mixture of oxidized organic compounds present in SOA, that significant amounts of both undissociated and dissociated forms of nitric and sulfuric acid should be present. This should allow for the occurrence of both specific and general acid catalysis in the organic phase, even though reactions in the aqueous phase are probably limited to specific acid catalysis.
| Organic phase | HNO3 | H2SO4 | ||
|---|---|---|---|---|
| D diss | D un | D diss | D un | |
| Hexane | 0.0014 ± 0.0004 | — | 0.0026 ± 0.0004 | — |
| 2-Hexanone | 0.072 ± 0.01 | 1.9 ± 0.6 | 0.016 ± 0.003 | — |
| 2-Butanone | 0.077 ± 0.003 | — | 0.019 ± 0.001 | — |
| Ethyl acetate | 0 | 3.3 ± 0.2 | 0 | 59 ± 10 |
| Propylene carbonate | 0.13 ± 0.06 | 31 ± 8 | — | — |
The pH of the organic phases measured after equilibration with a series of aqueous HNO3 or H2SO4 solutions ranging in concentration from 10−6 M to 10 M are shown in Fig. 3. It was assumed that the H+ activity measured by the pH strips was equal to the H+ concentration. This assumption ignores any activity effects, which are likely negligible compared to the uncertainty of the pH measurement (roughly one or two pH units), but may be important in extreme cases. These organic phase H+ concentrations are compared with those estimated using Ddiss values derived above and shown by the lines in Fig. 3. Note that not all systems had measurable Ddiss values and so there is not a line for each solvent. Although there is a fair amount of scatter, the H+ concentrations obtained from pH measurements largely line up with those obtained using Ddiss values. These results suggest that the mass-based measurements discussed above extend over a wider concentration range than measured. The notable exception is hexane, which is either overestimated by the mass-based measurements (probably because of the small amount of acid that partitioned to the organic phase) or underestimated by the pH measurements (perhaps due to activity effects or a relative lack of water in the organic phase).
![]() | ||
| Fig. 3 Measured and modeled organic phase H+ concentrations. Markers represent pH strip measurements and lines correspond to values derived from measured Ddiss values. | ||
:
C ratio of the organic phase. To supplement the measurements made here, relevant data sets from Stephen and Stephen (1963)27 were also fit to eqn (10) (Fig. S6†). The values of Dun (including measured and literature-derived values) were found to be approximately linear with respect to the O
:
C ratio of the organic phase, following the equationDun = 42 × (O : C) − 6 | (12) |
:
C value. Also included in the plot is a value of 8.2 calculated from the Henry's law coefficients of 1.6 × 105 M atm−1 for partitioning of gaseous HNO3 to water28 and 3.7 × 104 M atm−1 for partitioning to SOA formed from the oxidation of n-pentadecane by OH radicals.29 This SOA contained carbonyl, hydroxyl, carboxyl, ester, hydroperoxy, and peroxide groups and had an average molecular weight of 260 g mol−1 and an O
:
C ratio of 0.19.30 This data point, shown as a green diamond, lies slightly above the fitted line in Fig. 5, but still suggests that the values derived in this work are relevant to real aerosol systems despite the differences in functionality and molecular weight. Unlike the case for Dun, values of Ddiss do not correlate with O
:
C ratios, and although there is a general increase with increasing vol% water, there is too much scatter for this value to be a reliable predictor of this parameter (plots not shown).
![]() | ||
Fig. 5 Relationship between Dun and the O : C ratio of the organic phase for partitioning of HNO3. Values correspond to the following organic phase and O : C ratios: dibutyl ether (0.13),27 2-hexanone (0.17), SOA (0.19),29,30 ethyl ether (0.25),27 diethylene glycol dibutyl ether (0.25),27 ethyl acetate (0.50), and propylene carbonate (0.75). The error bars correspond to ± one standard deviation as determined from the linear fits used to calculate Dun, and in most cases fall within the symbol. The line is a linear least-squares fit to the data. | ||
| Organic phase | Water | Ammonium sulfate solution | ||
|---|---|---|---|---|
| D diss | D un | D diss | D un | |
| Hexane | 0.0026 ± 0.0004 | — | 0.0034 ± 0.0004 | — |
| 2-Hexanone | 0.016 ± 0.003 | — | 0.0439 ± 0.0044 | — |
| 2-Butanone | 0.019 ± 0.001 | — | 0.0067 ± 0.0007 | — |
| Ethyl acetate | 0 | 59 ± 10 | — | 68 ± 11 |
:
C ratio and vol% water of the organic phase, and to be influenced positively or negatively by the presence of ammonium sulfate, a dominant salt in atmospheric aerosol. Since studies were conducted with organic components having O
:
C ratios ranging from 0 to 0.75, the results should be applicable to the atmosphere, where the O
:
C ratio of organic aerosol typically ranges from 0.2 to 0.8.14 By fitting the data to a simple partitioning model, values of distribution coefficients, Ddiss and Dun, for partitioning of dissociated and undissociated acid to the organic phase were obtained, providing quantitative information on the amounts and form of partitioned acids. The results of this study suggest that partitioning of inorganic acids to predominantly organic phases can be significant and warrants further study. It would therefore be useful if future studies could expand on the limited selection of single organic phases employed here to include more complex mixtures and compounds with other functional groups, which might be possible by using more dilute acid solutions (which would be less likely to catalyze reactions of functional groups) and more sensitive acid detection methods such as ion chromatography or electrochemical measurements. Furthermore, it is an open question how well investigations of acidity and partitioning of the bulk solutions studied here apply to aerosol particles. For example, studies by Wei et al.31 suggest that pH gradients exist within particles, with protons accumulating at the surface and correspondingly more basic conditions towards the center of the particle. But on the other hand, studies by Dallemagne et al.32 showed that the pH of an organic shell in a phase-separated particle essentially matches that in the corresponding bulk solution. Comparison of the measurements made here with predictions of a liquid–liquid equilibrium model showed that the model seriously underpredicted the extent of acid partitioning to the organic phase, which was due at least in part because the model only considers the dissociated acid. The inclusion of undissociated acid should improve predictions, though more partitioning measurements are likely needed before this can be implemented. Depending on the acid and the particular model organic phase, these acids can exist primarily in either the dissociated or undissociated form, or both. One could then expect that for the complex mixture of oxidized organic compounds present in atmospheric aerosol that significant amounts of both forms of nitric and sulfuric acid should be present, allowing for the occurrence of both specific (H+) and general (HA) acid catalysis in the organic phase, even though reactions in the aqueous phase are probably limited to specific acid catalysis.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/d0ea00003e |
| This journal is © The Royal Society of Chemistry 2021 |