Simulation of SOA Formation from the Photooxidation of Monoalkylbenzenes in the Presence of Aqueous Aerosols Containing Electrolytes under Various NO x Levels

The formation of secondary organic aerosols (SOAs) from the photooxidation of three monoalkylbenzenes (toluene, ethylbenzene, and n-propylbenzene) in the presence of inorganic seeds (SO4-NH4-H2O system) under varying NOx levels has been simulated using the Unified Partitioning-Aerosol Phase Reaction (UNIPAR) model. The evolution of the volatility10 reactivity distribution (mass-base stoichiometric coefficient, αi) of oxygenated products, which were created by the nearexplicit gas kinetic mechanism, was integrated with the model using the parameters linked to the concentrations of HO2 and RO2 radicals. This dynamic distribution was applied to estimate the model parameters related to the thermodynamic constants of the products in multiple phases (e.g., the gas phase, organic phase, and inorganic phase) and the reaction rate constants in the aerosol phase. The SOA mass was predicted through the partitioning and aerosol chemistry processes of the oxygenated 15 products in both the organic phase and aqueous solution containing electrolytes, with the assumption of organic-inorganic phase separation. The prediction of the time series SOA mass (12-hr), against the aerosol data obtained from an outdoor photochemical smog chamber, was improved by the dynamic αi set compared to the prediction using the fixed αi set. Overall, the effect of an aqueous phase containing electrolytes on SOA yields was more important than that of the NOx level under our simulated conditions or the utilization of the age-driven αi set. Regardless of the NOx conditions, the SOA yields for the three 20 aromatics were significantly higher in the presence of wet electrolytic seeds than those obtained with dry seeds or no seed. When increasing the NOx level, the fraction of organic matter (OM) produced by aqueous reactions to the total OM increased due to the increased formation of relatively volatile organic nitrates and peroxyacyl nitrate like products. The predicted partitioning mass fraction increased as the alkyl chain length increases but the organic mass produced via aerosol phase reactions decreased due to the increased activity coefficient of the organic compounds containing longer alkyl chains. Overall, 25 the lower mass-base SOA yield was seen in the longer alkyl-substituted benzene in both the presence and absence of inorganic seeded aerosols. However, the difference of mole-base SOA yields of three monoalkylbenzenes becomes small because the highly reactive organic species (i.e., glyoxal) mainly originates from ring opening products without alkyl side chain. UNIPAR predicted the conversion of hydrophilic, acidic sulfur species to non-electrolytic dialkyl-organosulfate (diOS) in the aerosol. Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2018-963 Manuscript under review for journal Atmos. Chem. Phys. Discussion started: 26 October 2018 c © Author(s) 2018. CC BY 4.0 License.


Introduction
Anthropogenic volatile organic compounds (VOCs) have significant impacts on urban and regional atmospheric chemistry, despite fewer global emissions compared with biogenic VOCs (McDonald et al., 2018).As an important group of anthropogenic VOCs, aromatic hydrocarbons (HCs) are emitted from automobile exhaust (Zhang et al., 2018) and solvent use (Cheng et al., 2018) and are known to be precursors for secondary organic aerosols (SOAs), which are formed during Published by Copernicus Publications on behalf of the European Geosciences Union.
SOA formation has attracted substantial interest from scholars because of its vital role in affecting climate change (IPCC, 2015;Seinfeld and Pandis, 2016), urban visibility (Chen et al., 2012;Ren et al., 2018), and health (Requia et al., 2018).The prediction of SOA formation was first fulfilled by a gas-particle partitioning model.The partitioningbased SOA model uses two surrogate products (Odum et al., 1996) or several semivolatile surrogates (e.g., volatility basis set, VBS) (Donahue et al., 2006), with semiempirical parameters (e.g., the product stoichiometric coefficient, α; and gasparticle partitioning coefficient, K p ) for each HC system under a given NO x condition.Due to its simplicity and high efficiency, the partitioning-based model has been widely used in regional and global models.Nonetheless, the models and their predecessors are limited to predict SOAs formed from in-particle chemistry due to the loss of product structures, which govern the reactivity of organic species in the aerosol phase.Overall, regional air quality models have historically underestimated fine particulate matter in summertime (Appel et al., 2017;Huang et al., 2017) due to the lack of in-particle chemistry, particularly in the presence of an aqueous phase containing electrolytes (Ervens et al., 2011;Tsigaridis et al., 2014;Kelly et al., 2018).
A few models have attempted to implement in-particle chemistry into SOA models.For instance, Johnson et al. (2004Johnson et al. ( , 2005) ) simulated aromatic SOA chamber data, with a modified K p , to obtain experimentally comparable results, while the delayed simulated SOA mass indicated the occurrence of chemical reactions in the aerosol phase.McNeill et al. (2012) developed the Gas-Aerosol Model for Mechanism Analysis (GAMMA) to predict the formation of SOAs via aqueous-phase chemistry, which was further applied to the production of isoprene SOAs.Im et al. (2014) advanced the Unified Partitioning Aerosol Phase Reaction (UNIPAR) model, which predicted the SOA mass from partitioning processes and aerosol-phase reactions (reactions in both organic and inorganic phases and organosulfate (OS) formation).In that study, toluene and 1,3,5-trimethylbenzene SOAs were modeled using near-explicit products with the organicinorganic phase separation mode.Beardsley and Jang (2016) extended UNIPAR to simulate isoprene SOAs in the single homogeneously mixed phase (organic-inorganic mixture).Despite the reasonable prediction of SOA masses, UNIPAR faced inaccuracies in predicting time series SOA data due to the use of a fixed (nonage-driven) mass-based stoichiometric coefficient (α i ) set.
Age-driven functionalization and fragmentation alter the volatility and reactivity of products and their molecular struc-tures (Donahue et al., 2006;Rudich et al., 2007;Shilling et al., 2007;Hartikainen et al., 2018), which, in turn, varies the in-particle chemistry.Cappa and Wilson (2012) employed tunable parameters to kinetically demonstrate the evolution of SOA mass and the bulk oxygen-to-carbon atomic ratio (O:C ratio) during photochemical aging.However, oligomerization reactions in the aerosol phase were excluded.Donahue et al. (2011) developed a 2D-VBS method, which represented product aging by remapping the volatility and polarity (O : C ratio) of the products in 2-D space.Zhao et al. (2015) reported a discrepancy in the simulated toluene SOAs and α-pinene SOAs within the same 2D-VBS configuration, which may result from the different reactivities of the oxidation products of the precursors in aerosol-phase reactions.In this study, we have attempted to improve the UNIPAR model by using dynamic (age-driven) α i and applying the resulting model to predict the SOA formation of three monoalkylbenzenes (i.e., toluene, ethylbenzene, and npropylbenzene) under a wide range of environmental conditions (i.e., NO x , temperature, humidity, sunlight, and aerosol acidity).To consider the effect of the aging process on SOA formation, model parameters related to the organic molecular structures (i.e., the molecular weight, MW; and O : C ratio) and the α i set are calculated as the system ages, allowing for the internally dynamic estimation of the activity coefficient of the products (lumping species) in the aqueous phase containing electrolytes.Hence, the model is able to dynamically compute the partitioning coefficient of organics in the inorganic phase (K in ) by reflecting the photochemical evolution of the products in the gas phase and, consequently, improving SOA prediction.Organosulfate (OS), which has been identified in both laboratory and field studies (Hettiyadura et al., 2015;J. Li et al., 2016;Estillore et al., 2016;Chen et al., 2018), is an important chemical species due to its low volatility and ability to modulate the hygroscopicity of sulfate constituents.In the presence of acidic sulfate constituents, UNI-PAR also predicts the production of non-electrolytic sulfates (i.e., dialkyl-substituted OS, diOS) and the ensuing modification of aqueous-phase reactions.The feasibility of unified rate constants for aerosol-phase reactions was evaluated by extending the preexisting rate constants, which has been employed for toluene and 1,3,5-trimethylbenzene (Im et al., 2014) and isoprene (Beardsley and Jang, 2016), to the three monoalkylbenzenes in this study.

Experimental techniques
The SOA formation from the photooxidation of monoalkylbenzenes was conducted in the University of Florida Atmospheric PHotochemical Outdoor Reactor (UF APHOR) (Table 1).The concentrations of HCs, trace gasses (O x , SO 2 , and O 3 ), inorganic ions, aerosol acidity, and organic carbon (OC) of particles were monitored, as were the meteorological factors (i.e., relative humidity; temperature; and ul-  S7k a "Tol", "EB", and "PB" represent toluene, ethylbenzene, and n-propylbenzene oxidation experiments, respectively.b "E" or "W" that follows the experiment date represents the east or west chamber for the UF APHOR, respectively.c SOA data obtained from Im et al. (2014).d "SA", "wAS", and "dAS" denote directly injected sulfuric-acid-seeded aerosol, wet ammonium-sulfate-seeded aerosol, and dry ammonium-sulfate-seeded aerosol, respectively (dry: RH < ERH; wet: RH > ERH).For those without indication, SO 2 (in the unit of ppb) was injected into the chamber to generate sulfuric acidic seeds under the sun light.e SOA yield is estimated using Y SOA = OM/ HC, where OM is formed organic matter and HC is consumed HC.Yield in the table was estimated where SOA mass reached to the maximum over the course of the experiments.f This column denotes in which figures the corresponding data were used.The accuracy of relative humidity (RH) is 5 %.The accuracy of temperature is 0.5 K.
traviolet, UV, radiation).The configurations of the chamber and instrumentations were described by Im et al. (2014), J. Li et al. (2016), Beardsley and Jang (2016), Yu et al. (2017), andJiang et al. (2017).Aerosol acidity ([H + ], mol L −1 of aerosol) is monitored using colorimetry integrated with the reflectance UV-visible spectrometer (C-RUV) technique (Li et al., 2015) (Sect.S1 in the Supplement).The diOS concentration (µmol m −3 ) in an aerosol is estimated by the difference [H + ] obtained from ion chromatography (IC) interfaced with a particle-into-liquid sampler (PILS) (Li et al., 2015) and C-RUV method.Each HC was studied under at least two NO x levels (high NO x : HC/NO x < 5.5; low NO x : HC/NO x > 5.5) with or without inorganic-seeded aerosols (i.e., sulfuric acid, SA; or ammonium sulfate, AS).HONO was added into the system as a reaction initiator.To inves-tigate the effect of the liquid water content (LWC) on ASseeded SOA, two RH conditions were applied: (1) dry: RH < efflorescence RH (ERH) of the AS seed; (2) wet: RH > 50 % to prevent crystallization of AS seed.The ratio of organic matter (OM) to OC was experimentally determined to be 1.9 (Table 1, EB4), which was similar to the reported value of 2.0 for a series of toluene-NO x oxidation study (Kleindienst et al., 2007).

Model descriptions
The structure of the UNIPAR model is illustrated in Fig. 1.The simulation of aromatic SOA formation in the aqueous phase containing electrolytes was performed under the assumption of complete organic-inorganic phase separation.C and K denote the concentration and the partitioning coefficient of organic compounds, respectively, in gas phase (g), organic phase (or), and inorganic phase (in).k or,i denotes the reaction rate constant of oligomerization of organic compounds in the organic phase.k AC,i denotes the reaction rate constant of acid-catalyzed oligomerization of organic compounds in the inorganic phase and is determined as a function of aerosol acidity ([H + ]) and ambient humidity (RH)."OM" represents the concentration of organic matter.Subscripts "AR", "P ", and "T " indicate OM formed from aerosol-phase reactions, OM formed from the partitioning process, and total OM, respectively.Subscript i represents each lumping species.diOS represents the concentration of organosulfate -dialkyl sulfate (diOS) in this study.Bertram et al. (2011) modeled the separation RH (SRH) in the liquid-liquid phase of the mixture of organic and AS using the bulk O : C ratio.When ambient RH < SRH, the system undergoes organic-inorganic phase separation.The reported O : C ratios of the toluene, ethylbenzene, and npropylbenzene SOAs were 0.62 (Sato et al., 2012), 0.55 (Sato et al., 2012), and 0.45 (L.J. Li et al., 2016), respectively, which caused the corresponding SRH values to be 65 %, 80 %, and 93 %, respectively.Most RH for active photooxidation of HCs under ambient sunlight were under 65 %, which supported the assumption of organic-inorganic phase separation.In addition, as less soluble oligomers formed in the aerosol phase, an SRH higher than 65 % was more likely to be yielded.
To simulate age-dependent SOA formation, α i is reconstructed over time by a weighted average method using a pair of gas-phase oxidation compositions with different aging statuses: fresh composition and highly oxidized composition.The weighting factor at time = t is related to an aging scale factor (f A (t)), which is defined as where [RO 2 ] and [HO 2 ] represent the concentrations (ppb) of RO 2 and HO 2 radicals, respectively, and [HC] 0 represents the initial HC concentration (ppbC).The lower boundary of f A (t) (t = fresh) to determine the fresh α i set is equal to −7.2 at HC/NO x = 2 (high NO x levels) and −3.7 at HC/NO x = 14 (low NO x levels) for all three HCs.The upper boundary of f A (t) (t = highly aged) to determine the highly aged α i set is equal to −5.2 and −2.9 under the same high and low NO x levels, respectively.Both the fresh α i and highly aged α i are functions of HC/NO x .f A (t) is further converted into a fractional aging scale (f A (t) ) ranging from 0 (fresh composition) to 1 (highly aged composition) using a weight average method at each NO x level.Then, α i is dynamically reconstructed based on f A (t) under varying NO x conditions.
The molecular structures, including O : C i , MW (MW i ), and hydrogen bonding (HB i ) parameters, of each species (i) are also dynamically represented by a similar method, as shown in Sects.S3 and S4.

SOA formation: partitioning
The partitioning coefficient (K P ) from the gas (g) phase to the organic (or) phase (K or,i , m 3 µg −1 ) and from the g phase to the inorganic (in) phase (K in,i , m 3 µg −1 ) of each species is estimated using the following gas-particle absorption model (Pankow, 1994). and where R represents the gas constant (8.314J mol −1 K −1 ).
T represents the ambient temperature (K).MW or and MW in represent the average MW (g mol −1 ) of organic and inorganic aerosols, respectively.p o l,i represents the subcooled liquid vapor pressure (mm Hg) of a species, i.In the organic phase, we assume that the activity coefficient (γ or,i ) of a species (i) is unity (Jang and Kamens, 1998).In the inorganic phase, γ in,i is semi-empirically predicted by a regression equation, which was fit the theoretical activity coefficients of various organic compounds to RH, fractional sulfate (FS), and molecular structures (i.e., MW i , O : C i , and HB i ).FS is a numerical indicator for inorganic compositions related to aerosol acidity FS = , where [SO 2− 4 ] and [NH + 4 ] are the concentration of the total sulfate and the total ammonium, respectively).The theoretical activity coefficients were estimated at a given humidity and an aerosol composition through a thermodynamic model (Aerosol Inorganic-Organic Mixtures Functional Groups Activity Coefficients, AIOMFAC) (Zuend et al., 2011).
The statistical information for Eq. ( 4) is shown in Sect.S4 and Fig. S1 in the Supplement.The resulting K or,i and K in,i are employed to calculate the concentration (µg m −3 ) of the lumping species in multiple phases (C g,i , C or,i , C in,i , and Schell et al. ( 2001) developed a partitioning model to predict SOA formation.This model was reconstructed by Cao and Jang (2010) to include OM formed via aerosol-phase reactions (OM AR,i ) for a species (i), which is estimated in Sect.3.3.OM formed during the partitioning process (OM P ) is estimated by utilizing the mass balance shown in the following equation.
) is the effective saturation concentration and OM 0 represents the concentration (mol m −3 ) of the preexisting OM.MW oli,i represents the average MW of oligomeric products.Equation ( 5) is solved via iterations using the globally converging Newton-Raphson method (Press et al., 1992).

SOA formation: aerosol-phase reactions
The formation of OM AR,i is processed in both the organic and inorganic phases: oligomerization in the organic phase to form OM AR,or,i and oligomerization in the inorganic phase to form OM AR,in,i based on the assumption of a selfdimerization reaction (i.e., second-order reaction) (Odian, 2004) for organic compounds in media.Oligomerization in an aqueous phase can be accelerated under acidic environment (Jang et al., 2002).The oligomerization rate constants (L mol −1 s −1 ) in the organic phase and inorganic phase are k o,i and k AC,i , respectively, and the kinetic equations for oligomerizations are written as follows.
The bracketed terms in the equations indicate the conversion factors from aerosol-based concentrations (C or,i and C in,i : mol L −1 ) into air-based concentrations (µg m −3 ).The detailed derivations are shown in Sect.S5 and are illustrated in Fig. S2.ρ or and ρ in represent the density of the aerosol of organic and inorganic aerosol.ρ or was experimentally determined (EB4 in Table 1) to be 1.38 g cm −3 , which was similar to the reported value of 1.4 g cm −3 for aromatic SOA (Nakao et al., 2011;Chen et al., 2017;Ng et al., 2007).ρ in is obtained from a regression equation through the extended aerosol inorganic model (E-AIM) (Clegg et al., 1998).Due to atmospheric diurnal patterns (high RH at nighttime to low humidity during daytime), it is likely that the RH changes would be based on inorganic aerosol ERH.UNIPAR internally predicts the ERH using the equation derived by Colberg et al. (2003).k AC,i in Eq. ( 7) is estimated based on a semiempirical model developed by Jang et al. (2005) as a function of species reactivity (R i ), protonation equilibrium constant (pK BH + i ), excess acidity (X), water activity (a w ), and proton concentration ([H + ]), which are estimated by the E-AIM.
In the organic phase, k o,i is estimated by excluding the X and a w [H + ] terms.The formed OM AR can be calculated as a sum of OM AR,or,i and OM AR,or,i for each species assuming that OM AR is irreversibly formed and nonvolatile (Kleindienst et al., 2006;Cao and Jang, 2010).

Organosulfate formation
In the presence of aqueous acidic sulfate, UNIPAR predicts the formation of diOS ([diOS] model ) to compute the change in aerosol hygroscopicity and acidity.At each time step, free electrolytic sulfate ([SO 2− 4 free where f diOS represents the diOS conversion factor introduced by Im et al. (2014), which was semi-empirically determined to be 0.071 in this study.N diOS represents the numeric parameter for scaling lumping groups based on the effectiveness of the chemical species to form diOS.For example, the diOS scale factor is 1 for each alcohol and aldehyde group and 2 for each epoxide group (see Tables S1-S3 for functional groups).Then, N diOS is summed at each time step and applied to Eq. ( 9).

Prediction of SOA mass under the evolution of oxygenated products
As reported in former studies, the kinetic mechanism tends to underestimate the decay of aromatic HCs because of the low prediction of OH radicals (Johnson et al., 2005;Bloss et al., 2005).In this study, the addition of artificial OH radicals varies with the HC/NO x ratio by fitting the predicted decay of HCs using the kinetic mechanism in the experimental measurements.The time profiles of the decays of the three HCs are shown in Fig. S3 (Sect.S6).
When the NO x level is very low, the maximum additional OH radical production rate for monoalkylbenzenes is 2 × 10 8 molecules cm 3 s −1 , which is less than 4 × 10 8 (Bloss et al., 2005) but similar to the value reported by Im et al. (2014).When HC/NO x < 3, no addition of artificial OH radicals is needed for the chamber simulation of the decay of monoalkylbenzenes.For the make-up OH production rate constants of all three HCs under varying NO x , the mathematical weighting equation is written below: dynamic makeup OH rate = e 0.6×(HC/NO x ) e 0.6×(HC/NO x ) + 50 × 2.0 × 10 8 molecules cm 3 s −1 .(10) In our model, we assume that the oxidation of products progresses in the gas phase.Lambe et al. (2012) reported that the transition point of n-C 10 SOAs from a functionalizationdominant regime to a fragmentation-dominant regime is approximately 3 d (photochemical equivalent age under an atmospheric OH exposure of 1.5 × 10 6 molecules cm −3 ).Under this criterion, we exclude the aging of nonvolatile aerosol products (OM AR ).However, the oxidation of aerosol products for longer periods of time may decrease the volatility (George and Abbatt, 2010;Jimenez et al., 2009).
Figure 2 illustrates the evolution of the volatilityreactivity-based distribution of the mass-based stoichiometric coefficient (α i ) of toluene at the two different NO x levels (HC/NO x = 2.9 and 10.5).Collectively, most α i values at both NO x levels tend to decline as the reaction time lapses (Fig. 2a vs. Fig.2b; Fig. 2c vs. Fig.2d) since the evolution of some semivolatile organic compounds (SVOCs) forms very volatile molecules (i.e., CO 2 , formic acid, and formaldehyde).For example, the α i values of highly reactive carbonyls with high volatility (GLY and MGLY in Table S1 of Sect.S7) are high under the fresh condition and significantly decline as the system ages, because they undergo fast photolysis under sunlight (George et al., 2015;Henry and Donahue, 2012).Consequently, the decay of these highly reactive species leads to the decrease in the production of OM AR .
The high NO x level delays the oxidation of gas-phase products.Similar trends in the α i set can be found for ethylbenzene and n-propylbenzene, as shown in Sect.S7 (Tables.S2 and S3, Figs.S4 and S5).The α i of highly reactive species (e.g., GLY, 8VF, 3M, and 5S) decreases by increasing the NO x level due to the suppression of the HO x cycle via the reaction of NO 2 with OH radicals.As seen in Fig. 2d, some medium reactivity species -i.e., 2-methyl-4-oxo-3-nitro-2butenoic acid (3M), 2-methyl-4-oxo-2-butenoic acid (6M), and acetyl-3-oxopropanoate (7M) -start to form as NO x decreases.
In Fig. 3, the comparison between the simulations of SOA formation from toluene oxidation is based on two different α i -reconstruction strategies: dynamic α i and fixed α i .A clear improvement in the prediction of SOA formation is demonstrated when comparing the SOA mass using dynamic α i to that using fixed α i .The aged SOA growth from the three systems (i.e., low NO x level, Fig. 3a and d; moderate NO x level, Fig. 3b and e; and high NO x level with an inorganic seed, Fig. 3c and f) is even smaller than that predicted with the less-aged α i set, which is obtained when precursors being consumed make up half of the total simulated consumption.Our model simulation against the chamber data suggests that while aging may alter aerosol compositions (Fig. 2), it does not always increase SOA yields.Traditionally, the SOA mass has been predicted using fixed thermodynamic parameters (i.e., K p and α i ), which is inadequate when reflecting upon practical scenarios, where oxygenated product distributions vary dynamically with oxidation.

Effects of aerosol acidity and LWC on SOA formation
In the model, aerosol acidity was expressed using a fractional free sulfate (FFS), which is defined as . Humidity can influence both aerosol acidity and LWC, which are the model parameters in UNIPAR.Thus, UNIPAR has the capability to decouple the effect of aerosol acidity and humidity, as shown in Fig. 4 for toluene SOA.The impact of aerosol acidity and humidity on the yields of SOAs derived from ethylbenzene and n-propylbenzene is illustrated in Fig. S6 (Sect.S8).The dramatic difference in SOA yields appears between the RH above ERH and the RH below ERH.The LWC disappears below ERH, and there are no aqueous reactions.For example, the observed SOA yield of ethylbenzene with effloresced AS was significantly smaller than that with wet AS: 11 % (EB8 in Table 1) vs. 30 % (EB9 in Table 1).Kamens et al. (2011) and Liu et al. (2018) reported a significantly lower yield of toluene SOA for dry AS-seeded aerosols compared with its wet counterpart.The partitioning of polar carbonaceous products increases with increasing LWC and, thus, aqueous reactions.In the presence of wet aerosols, SOA yields gradually increase with increasing FFS (increasing acidity) at a given RH due to acid-catalyzed oligomerization.The oxygenated products of toluene are relatively more polar than those of ethylbenzene or propylbenzene and positively attributed to the increase in OM AR .
Compared to isoprene SOAs reported by Beardsley and Jang (2016), the impacts of humidity and acidity on the SOA yields of monoalkylbenzenes in this study are relatively weaker above the ERH (Fig. 4), except for the highly acidic conditions under high humidity.In this study, aromatic SOA mass is attributed to a few highly reactive species, such as GLY.Other aromatic oxidation products partitioned in the aerosol phase have moderate solubility and they are slow to react in the aqueous phase.Isoprene products are more hygroscopic than aromatic products and even mixable with an aqueous phase containing electrolytes.The reactions of medium reactivity polar products that formed during isoprene oxidation can be accelerated by an acid catalyst with higher sensitivities to acidity and humidity.

Figure 5 illustrates the time profiles of the predicted concentrations of diOS ([diOS] model ) and protons ([H + ]) with the measured concentrations of diOS ([diOS] exp ), [NH +
4 ], and [SO 2− 4 ] for different aromatic HCs under given experimental conditions (Fig. 5a-c).For the three SAseeded SOA experiments, the fractions of diOS to the total sulfate amount are 0.09, 0.15, and 0.06 for toluene (Exp.Tol8, HC/NO x = 2.9, FS changing from 0.64 to 0.39), ethylbenzene (Exp.EB7, HC/NO x = 12.3, FS changing from 0.82 to 0.46), and n-propylbenzene (Exp.PB5, HC/NO x = 14.4,FS changing from 0.76 to 0.38), respectively.The [diOS]   aerosols, was higher than those in Exp.Tol8 and Exp.PB5, indicating that the acidic condition was favorable for the formation of diOS (Surratt et al., 2010;Lin et al., 2013).Beardsley and Jang (2016) reported that the diOS fraction for isoprene SOAs was 0.26 (HC/NO x = 32.5,FS changing from 0.69 to 0.47), which was more than that for toluene SOAs, indicating that the oxidation products of isoprene may contain more reactive species to form diOS.For example, IEPOX products in isoprene SOAs are known to be reactive to SA (Budisulistiorini et al., 2017).Additionally, isoprene aerosol products are mixable with aqueous solutions containing electrolytes, and they can more effectively form diOS compared to the aromatic products in liquid-liquid phase separation.
To estimate the potential upper boundary of the concentration of diOS ([diOS] max ) in monoalkylbenzene SOA, the aerosol composition was predicted by the model in the presence of SA aerosols (without neutralization with ammonia) under the given experimental conditions shown in Fig. 5.The resulting diOS fractions were 0.29 (OM-to-sulfate mass ratio (OM : sulf) = 1.4), 0.25 (OM : sulf = 1.4), and 0.12 (OM : sulf = 0.7) for toluene, ethylbenzene, and npropylbenzene, respectively.The aerosol acidity of the ambient aerosol is generally lower than ammonium hydrogen sulfate (AHS), and, thus, the diOS fraction in ambient air would be much lower than the estimated upper boundary.Figure 5 suggests that the change in both aerosol acidity and hygroscopicity by the formation of non-electrolytic sulfate is important to predict SOA mass.

Effect of NO
x on SOA formation in the presence of an aqueous phase containing electrolytes Figure 6 shows the impact of NO x on the three aromatic SOAs in this study by producing SOAs at two different NO x levels in the presence and absence of SO 2 .Overall, regardless of the inorganic seed conditions, both the chamber observation and model simulation suggest that increasing the NO x level leads to the decreased formation of SOAs.This trend in the absence of inorganic seed aerosols has also been observed multiple times (Li et al., 2015;Ng et al., 2007;Song Figure 3.Comparison between simulated SOA mass using the fixed α i (α i is obtained when precursor being consumed make up half of the total simulated consumption) and dynamic α i (α i evolving as photooxidation) under (a) low-NO x condition (Exp.Tol6, HC/NO x = 10.5),(b) moderate-NO x condition (Exp.Tol7, HC/NO x = 6.4), and (c) high-NO x condition (Exp.Tol8, HC/NO x = 2.9 with sulfuric acid (SA)seeded aerosol).(d, e, f) represent the time-dependent SOA growth curve (SOA mass concentration against consumed HC in the unit of µg m −3 ) corresponding to the experimental conditions of (a, b, c), respectively.The solid circle represents the experimental measurements.The SOA mass is corrected for particle loss to the chamber wall.The experimental conditions are available in Table 1.et al., 2005).By increasing [NO x ], the path of an RO 2 radical progresses to the formation of organonitrate and peroxyacyl nitrate (PAN) products, which are less reactive to aerosolphase reactions.They are relatively volatile and, thus, insignificantly attributed to partitioning SOA mass.For example, the SOA yields under the low NO x level (HC(ppbC)/NO x (ppb) = 9.1-14.8,Table 1) in the presence of SO 2 , with a similar degree of ammonia titration (i.e., similar FS values by the end of the experiments), were higher than those without seeded aerosols: 42 % for toluene (Exp.Tol1, FS = 0.44), 26 % for ethylbenzene (Exp.EB1, FS = 0.37), and 66 % for propylbenzene (Exp.PB1, FS = 0.43).The impact of aerosol acidity was even greater for SOAs produced under a high NO x level (HC(ppbC)/NO x (ppb) = 2.8-5.0):65 % for toluene (Exp.Tol3, FS = 0.43), 146 % for ethylbenzene (Exp.EB5, FS = 0.39), and 77 % for propylbenzene (Exp.PB3, FS = 0.40).SOA formation under high-NO x conditions is generally more sensitive to aerosol acidity compared to that at low NO x levels (Fig. 6a-c vs. Fig.6d-f).The fractions of medium reactivity products are relatively high in high NO x levels, and their reactions in aerosol phase can be accelerated by an acid catalyst.Figure S7 (Sect.S9) also illustrates the simulation of SOA mass against chamber data under varying humidity, NO x levels, and aerosol acidity (Table 1).for three monoalkylbenzenes to important variables -i.e., humidity (Ai vs. Bi for AHS and Ci vs. Di for AS), temperature (Ai vs. Ei for AHS and Fi vs. Gi without inorganic aerosol), aerosol acidity (Ai vs. Ci at RH = 45 % and Bi vs. Di at RH = 65 %), and NO x levels (Ai vs. Hi with AHS-seeded aerosols and Fi vs. Ii without inorganic-seeded aerosols).The most drastic change appears by changing the temperature from 298 K (Ai) to 273 K (Ei).The SOA yield is known to increase by 20 %-150 %, which results from a 10 K decrease in temperature (Sheehan and Bowman, 2001).For all SOAs, noticeable changes are shown between the absence (Fi) and presence (Ai) of wet inorganic seeds, while a minor change appears between wet AHS (Ai) and wet AS (Ci).Within the wet acidic aerosols (Ai vs. Bi and Ci vs. Di), the effect of RH is insignificant in our simulation, as discussed in Sect.4.2.Although the impact of NO x (Ai vs. Hi and Fi vs. Ii) is less than that of temperature and inorganic seeds, SOA yields are still significantly altered, as discussed in Sect.4.4.The panel ii series in Fig. 7 illustrates SOA growth curves under various conditions shown in panel i.Overall, the simulated SOA yields (slopes) increase with a www.atmos-chem-phys.net/19/5719/2019/decreased alkyl chain length (toluene > ethylbenzene > npropylbenzene), which is consistent with our chamber observations (Table 1).Although the decrease in the vapor pressure of products benefits the increases in OM P as the alkyl chain length increases, the increase in the activity coefficient of the organic products containing longer alkyl chains in aqueous phase is unfavorable to the formation of OM AR via aqueous reactions.However, the difference of mole-based SOA yields of three monoalkylbenzenes becomes small because the highly reactive organic species (i.e., glyoxal), which are produced through ring opening reactions without an alkyl side chain, are significantly attributed to OM AR .Fig. 7ii confirms that the effect of an aqueous phase containing electrolytes on SOA yields is more critical than that of the NO x level under our simulated conditions.

Sensitivity of model prediction to major variables and model uncertainty
To determine the model sensitivity to these parameters, simulations (Exp.Tol9 in Table 1) were performed by increasing/decreasing vapor pressure (V p ), the enthalpy of vaporization (H vap ), γ in,i , and k AC,i by factors of 1.5, 1.1, 2, and 2, respectively.The corresponding changes in the SOA mass are −7.2 %/6.8 %, −1.9 %/1.7 %, −8.8 %/4.7 %, and 2.5 %/−3.5 %, respectively.The change in SOA mass from the reference for each simulation is shown in Fig. S8 (Sect.S10).
The uncertainty associated with the group contribution method used for V p estimation is a factor of 1.45 (Zhao et al., 1999).H vap has a reported error of 2.6 % (Kolska et al., 2005).γ in,i is estimated as a function of O : C, MW, RH, and FS (Eq. ).k AC,i is semi-empirically calculated based on [H + ], LWC, and species reactivity (Eq.8).The E-AIM is performed to estimate the LWC, which is reliable and based on a broadly used water activity dataset (Zhang et al., 2000).Yet, the inorganic thermodynamic models including E-AIM performed inadequately in the prediction of [H + ] under low-RH and ammonia-rich conditions (FS < 0.55) (Li and Jang, 2012).
In the model, non-electrolytic diOS was predicted and applied to the prediction of LWC and [H + ], which subsequently affect aerosol growth via aqueous reactions.Typically, the monoalkyl sulfate is identified as a product of the esterification of SA with reactive species (Hettiyadura et al., 2015;J. Li et al., 2016;Estillore et al., 2016;Chen et al., 2018).It is possible that monoalkyl sulfates can influence LWC and aerosol acidity differently from sulfuric acid, although they are strongly acidic and hygroscopic.Although Noziere et al. (2010) reported that OS could be produced by the reactions of GLY and sulfate radicals in the presence of aqueous AS under UV light, the amounts of formed monoalkyl OS and their influence on aerosol hygroscopicity are still not clear.Some other factors in recent investigations, such as organic vapor wall loss and aerosol viscosity, have not been accounted for by the UNIPAR model.The loss of organic vapor to the Teflon chamber wall can compete with the gasparticle partitioning process and the reactions in both the gas phase and aerosol phase to initiate a negative bias in the experimental measurements (Zhang et al., 2014;Mcvay et al., 2014).The modeling of the gas-wall process of semivolatile organic compounds can improve the prediction of SOA mass in regional scales.In addition, an increased aerosol viscosity via aging could modify the diffusivity of the partitioned The degree of neutralization is indicated by FS, ranging from 1 (for sulfuric acid) to 0.33 (for ammonium sulfate)."SA" stands for experiment with direct-injection sulfuric-acid-seeded aerosols.The ions and diOS concentrations were corrected for the particle loss to the chamber wall.The experimental conditions are available in Table 1.d-f) conditions in the presence (red color)/absence (green color) of SO 2 -derived sulfuric-acid-seeded aerosol.Solid, dashed, and dotted lines denote the total organic matter (OM T ), the OM from partitioning only (OM P ), and the OM from the aerosol-phase reactions (OM AR ), respectively.The degree of ammonia neutralization with sulfuric acid is indicated by the FS end , which is the FS at the end of the experimental run.The FS end ranges from 1 (for sulfuric acid) to 0.33 (for ammonium sulfate).The uncertainty associated with experimentally measured OM is about 9 %.The SOA mass was corrected for the particle loss to the chamber wall.The experimental conditions are available in Table 1.organic molecules (Abramson et al., 2013) and the reaction rate constant for oligomerization in the aerosol phase.

Conclusions and implications
Despite numerous studies in SOA characterization and formation mechanisms, substantial biases between the simulated and field-measured SOA mass were still found (Hodzic et al., 2016) due to the inadequacy of handling the dynamic multigenerational aging (Jathar et al., 2016) and aqueous reactions of the oxygenated products in the presence of an aqueous phase containing electrolytes (Ervens et al., 2011).In this study, the UNIPAR model addressed those issues using a dynamic age-driven α i set, multiphase partitioning of organic compounds, and in-particle chemistry.Although the utilization of the age-driven α i set improves the time series  .The initial concentrations of monoalkylbenzenes, preexisting OM (OM 0 ), NH 4 HSO 4 (AHS)seeded aerosol, and (NH 4 ) 2 SO 4 (AS)-seeded aerosol are 20 ppb, 2, 20, and 20 µg m −3 , respectively.The gas-phase simulation used the sunlight on 14 June 2018 (Exp.Tol9 in Table 1).OM P (diagonal stripes fill) and OM AR (solid fill) represent the organic matter from the partitioning process and aerosol-phase reactions.(ii) shows the time-dependent SOA growth curve for three monoalkylbenzenes under corresponding simulation conditions (top table).The concentrations that follow the legends refer to the mass concentrations of the consumed monoalkyl-substituted benzenes in each simulation under the high or low-NO x conditions.prediction of SOA mass, as shown in Fig. 3, the photochemical evolution of the gas-phase products via monoalkylbenzene oxidation (Figs. 2, S4, and S5) does not increase the SOA mass, as is commonly suggested.Overall, the effect of an aqueous phase containing electrolytes on SOA formation was more critical than that of the NO x level under our simulated conditions.By adding a wet inorganic seed to the nonseed SOA system, the mass-based SOA yields under high NO x levels increase more than those under low-NO x conditions (Fig. 6 in Sect. 4.4).The vapor pressure of volatile organonitrate and PAN-like species, which are formed at high NO x levels, is not low enough to increase partitioning SOA mass (Fig. 7Aii).Thus, SOA yields decreased by increasing NO x levels.Overall, both simulation and chamber data show that monoalkylbenzene SOA yields increase with a decreased alkyl chain length: toluene > ethylbenzene > npropylbenzene.This difference is most noticeable in the presence of an inorganic seed at high NO x levels (Fig. 7Aii and Hii) (Colberg et al., 2003).
Due to the pervasiveness and relatively high concentration of toluene in the urban situation, where HC/NO x < 5.5 and wet inorganic seeds typically exist, the importance of toluene SOAs to the urban SOA burden can increase.The oxidation products from aromatic HCs can also involve cloud condensation nuclei activity due to their high reactivity via heterogeneous chemistry (Molteni et al., 2018), resulting in a change in the properties of clouds and fog and the urban radiation balance (Gordon et al., 2016).The unified aerosol-phase reaction rate constants for three monoalkylbenzenes represent the feasibility of applying the UNIPAR model to more aromatic systems (dialkyl benzenes and trialkyl benzenes) and the complex urban mixture.
Data availability.The chamber data and simulation results as part of this publication are available upon request.
Author contributions.MJ designed the experiments and CZ and ZY carried them out.MJ developed the model, and CZ performed the calculation of model parameters and the simulations.CZ and MJ prepared the manuscript with contributions from ZY.
Competing interests.The authors declare that they have no conflict of interest.

Figure 1 .
Figure 1.Simplified scheme of the UNIPAR model.[HC] 0 represent the initial hydrocarbon (HC) concentration.The dynamic mass-based stoichiometric coefficient (dynamic α i ), the consumption of HCs ( HCs), the concentration of hydroperoxide radical ([HO 2 ]), and the concentration of organic peroxyl radical ([RO 2 ]) are simulated from the gas kinetic mechanism (MCM v3.3.1).The aging scale factor (f A ) is represented as a function of [HO 2 ], [RO 2 ], and [HC] 0 , which is detailed in Sect.3.1.C and K denote the concentration and the partitioning coefficient of organic compounds, respectively, in gas phase (g), organic phase (or), and inorganic phase (in).k or,i denotes the reaction rate constant of oligomerization of organic compounds in the organic phase.k AC,i denotes the reaction rate constant of acid-catalyzed oligomerization of organic compounds in the inorganic phase and is determined as a function of aerosol acidity ([H + ]) and ambient humidity (RH)."OM" represents the concentration of organic matter.Subscripts "AR", "P ", and "T " indicate OM formed from aerosol-phase reactions, OM formed from the partitioning process, and total OM, respectively.Subscript i represents each lumping species.diOS represents the concentration of organosulfate -dialkyl sulfate (diOS) in this study.

4. 5
Figure7illustrates the sensitivity of the SOA mass simulated at relatively low concentration of HCs (20 ppb) (panel i) for three monoalkylbenzenes to important variables -i.e., humidity (Ai vs. Bi for AHS and Ci vs. Di for AS), temperature (Ai vs. Ei for AHS and Fi vs. Gi without inorganic aerosol), aerosol acidity (Ai vs. Ci at RH = 45 % and Bi vs. Di at RH = 65 %), and NO x levels (Ai vs. Hi with AHS-seeded aerosols and Fi vs. Ii without inorganic-seeded aerosols).The most drastic change appears by changing the temperature from 298 K (Ai) to 273 K (Ei).The SOA yield is known to increase by 20 %-150 %, which results from a 10 K decrease in temperature(Sheehan and Bowman, 2001).For all SOAs, noticeable changes are shown between the absence (Fi) and presence (Ai) of wet inorganic seeds, while a minor change appears between wet AHS (Ai) and wet AS (Ci).Within the wet acidic aerosols (Ai vs. Bi and Ci vs. Di), the effect of RH is insignificant in our simulation, as discussed in Sect.4.2.Although the impact of NO x (Ai vs. Hi and Fi vs. Ii) is less than that of temperature and inorganic seeds, SOA yields are still significantly altered, as discussed in Sect.4.4.The panel ii series in Fig.7illustrates SOA growth curves under various conditions shown in panel i.Overall, the simulated SOA yields (slopes) increase with a

Figure 4 .
Figure 4. Simulated toluene SOA yields (Y SOA = OM/ HC at the end of the simulation, where OM is formed organic matter and HC is consumed HC) as a function of relative humidity (RH: 0.1-0.9)and fractional free sulfate (FFS: 0-1), where FFS = ([SO 2− 4 ] − 0.5[NH + 4 ])/[SO 2− 4 ], which is another numerical indicator that is used to estimate aerosol acidity ([H + ]) in the inorganic thermodynamic model.The RH and FFS are fixed in the simulations.The gas-phase simulations are based on the experimental condition of 14 June 2018 (Exp.Tol9 in Table 1) (initial HC concentration = 5 ppb, HC/NO x = 2, preexisting OM (OM 0 ) mass concentration = 2 µg m −3 , sulfate mass concentration = 20 µg m −3 , and the mass ratio of the consumed HC to sulfate ( HC : sulf) = 1).

Figure 6 .
Figure 6.Time profiles of measured and modeled SOA mass concentrations (µg m −3 ) for toluene, ethylbenzene, and n-propylbenzene SOA under low-NO x (a-c)/high-NO x (d-f) conditions in the presence (red color)/absence (green color) of SO 2 -derived sulfuric-acid-seeded aerosol.Solid, dashed, and dotted lines denote the total organic matter (OM T ), the OM from partitioning only (OM P ), and the OM from the aerosol-phase reactions (OM AR ), respectively.The degree of ammonia neutralization with sulfuric acid is indicated by the FS end , which is the FS at the end of the experimental run.The FS end ranges from 1 (for sulfuric acid) to 0.33 (for ammonium sulfate).The uncertainty associated with experimentally measured OM is about 9 %.The SOA mass was corrected for the particle loss to the chamber wall.The experimental conditions are available in Table1.

Figure 7 .
Figure7.The simulated SOA mass (i) for toluene (Tol), ethylbenzene (EB), and n-propylbenzene (PB) under different conditions.The simulation conditions are listed in the top table.The initial concentrations of monoalkylbenzenes, preexisting OM (OM 0 ), NH 4 HSO 4 (AHS)seeded aerosol, and (NH 4 ) 2 SO 4 (AS)-seeded aerosol are 20 ppb, 2, 20, and 20 µg m −3 , respectively.The gas-phase simulation used the sunlight on 14 June 2018 (Exp.Tol9 in Table1).OM P (diagonal stripes fill) and OM AR (solid fill) represent the organic matter from the partitioning process and aerosol-phase reactions.(ii) shows the time-dependent SOA growth curve for three monoalkylbenzenes under corresponding simulation conditions (top table).The concentrations that follow the legends refer to the mass concentrations of the consumed monoalkyl-substituted benzenes in each simulation under the high or low-NO x conditions.

Table 1 .
Experimental conditions and resulting SOA chamber data from the monoalkylbenzenes photooxidation experiments performed under various NO x conditions with/without inorganic-seeded aerosol in the dual outdoor UF APHOR chambers.
model reasonably agrees with[diOS]exp .The aerosols in Exp.Tol8 and Exp.PB5 show the cessation in diOS formation at approximately 10:00 EST since they became effloresced due to the neutralization of SA with ammonia under the reduction in humidity during the daytime.The diOS fraction in Exp.EB7, which contained wet acidic www.atmos-chem-phys.net/19/5719/2019/Atmos.Chem.Phys., 19, 5719-5735, 2019