Effects of SO2 on optical properties of secondary organic aerosol generated from photooxidation of toluene under different relative humidity conditions

Secondary organic aerosol (SOA) has great impacts on air quality, climate change and human health. The composition and physicochemical properties of SOA differ greatly because they form under different atmospheric conditions and from various precursors as well as differing oxidation. In this work, photooxidation experiments of toluene were performed under four conditions (dry, dry with SO2, wet and wet with SO2) to investigate the effect of SO2 under different relative humidities on the composition and optical properties of SOA at wavelengths of 375 and 532 nm. According to our results, the increase in humidity enhances not only light absorption but also the scattering property of the SOA. Oligomers formed through multiphase reactions might be the reason for this phenomenon. Adding SO2 slightly lowers the real part of the complex refractive index, RI(n), of toluene-derived SOA (RI(n)dry,SO2 < RI(n)dry, RI(n)wet,SO2 < RI(n)wet), which might be a result of the partitioning of low-oxidation-state products. The imaginary part of the complex refractive index, RI(k), is enhanced under dry conditions with SO2 compared to that of only dry conditions, which might be due to acid-catalyzed aldol condensation reactions. Wet conditions with SO2 shows the combined effect of SO2 and humidity. The extinction properties of toluene-derived SOA under wet conditions with SO2 increased by approximately 30 % compared to that of toluene-derived SOA formed under dry conditions. Our results suggest that various atmospheric conditions will affect the composition and optical proprieties of SOA, which has significant implications for evaluating the impacts of SOA on the rapid formation of regional haze, global radiative balance and climate change.


Introduction
Secondary organic aerosol (SOA) accounts for a major fraction of the atmospheric fine particulate matters (PM 2.5 ) and has significant impacts on air quality, climate change and human health (Seinfeld and Pandis, 2006;Jimenez et al., 2009; could influence the radiative balance directly by scattering and absorbing solar and terrestrial radiation and indirectly by impacting cloud formation and their lifetimes (Ciais et al., 2013;Andreae and Gelencser, 2006;Moise et al., 2015). On the urban and regional scales, SOA contributes to the degradation of visibility and adversely influences human health (Watson, 2002;Pope et al., 2002). Light-absorbing aerosols (including black carbon, BC, mineral dust and brown carbon, BrC) are recognized as playing important roles in climate radiative forcing because of the strong dependence of their optical properties on the aerosol composition, the complexity of their production and the poor constraints on their contribution to radiative forcing Wang et al., 2013). Quantifying the optical properties of SOA is one of the key problems in accessing anthropogenic pollution data for visibility, air quality and climate change, as well as one of the most urgent issues in atmospheric sciences Moise et al., 2015;Andreae and Gelencser, 2006).
Aromatic compounds, accounting for 20 %-40 % v/v of gasoline fuel, are part of the representative anthropogenic volatile organic compounds (VOCs) and play an important role in the formation of tropospheric ozone and SOA Ng et al., 2007;Zhang et al., 2015;Odum et al., 1997). Among these aromatic hydrocarbons, toluene is thought to be one of the most important SOA precursors because of its strong emissions and SOA formation potential (Hildebrandt et al., 2009;Odum et al., 1997Odum et al., , 1996Tao et al., 2017). The optical properties of aromatic-derived SOA have received significant attention over the last decade, most of which focused on the effect of nitrogen oxides (NO x ) on the optical properties of aromatic SOA (Li et al., 2014;Kim and Paulson, 2013;Lin et al., 2015;Nakayama et al., 2013). In particular, K.  studied the effects of multiphase processes on the optical properties of m-xylene SOA (K. . However, the optical properties of toluene-derived SOA under various atmospheric relevant conditions remain unclear, which limits our understanding of severe and complex air pollution.
Recent studies have shown that SOA is a dominant fraction of the total atmospheric aerosol burden, particularly in heavily polluted areas (Guo et al., 2014;Tao et al., 2017). Reported field experiments have shown that the atmospheric concentration of sulfur dioxide (SO 2 ) could be up to nearly 200 ppb in heavy haze pollution episodes in China, while the formation and growth rates of SOA and sulfate, which significantly contribute to severe haze pollution, were much faster than that during clean periods (H. . According to An's study, the concentrations of NO 2 and nitrate were also quite high in China, especially in the North China Plain (NCP) . Previous studies have mostly focused on the enhancement of SOA yields introduced by the presence of SO 2 for isoprene, α-pinene and anthropogenic precursors (Santiago et al., 2012;Kleindienst et al., 2006;Liggio and Li, 2013). They also revealed that an enhancing effect of SOA yields is because of the acidic aerosol prod-ucts of SO 2 and the formation of new particles. Surratt et al. (2010) found that reactions of isoprene epoxydiols in the presence of acidic aerosols produced from the oxidation of SO 2 could be a substantial source of "missing urban SOA" (Surratt et al., 2010). Jaoui et al. (2008) showed that the addition of SO 2 would make the color of the SOA generated via the photooxidation of α-pinene and toluene darker and browner, and oligomers and nitrogenous organic compounds were detected in the SOA extracts (Jaoui et al., 2008). Nakayama et al. (2018) investigated the effect of SO 2 on the optical properties of isoprene under various NO x concentrations and oxidation pathways (Nakayama et al., 2015(Nakayama et al., , 2018. However, more research is needed on the role of SO 2 in the subsequent optical properties of SOA formed from various VOCs. Field studies have shown that haze events were often accompanied by high relative humidity (RH) (Sun et al., 2016a, b), and organic aerosols were mostly liquid (Shiraiwa et al., 2017;Liu et al., 2017). Multiphase reactions also play an important role in SOA formation, as shown in previous studies. For example, Jia and Xu (2018) have shown that the yield of toluene SOA almost doubled at a relative humidity of 85 % compared to dry conditions (Jia and Xu, 2018). K.  found that multiphase reactions would enhance light scattering and radiative forcing of the m-xylene SOA (K. . Understanding the impact of complicated air conditions on SOA generation, especially the effect of polluting gases (e.g., SO 2 ) and phase state, is urgently needed for forecasting, assessing and controlling air pollution.
In this work, we investigated the effect of SO 2 on the optical properties and chemical composition of the SOA derived from toluene under different humidity. The results will greatly help the evaluation of the toluene-derived SOA on atmospheric visibility and climate change under complicated pollution conditions. More importantly, the data in the current study will be highly useful for the simulation of models and field observations performed under various pollution conditions.

Smog chamber experiments and online measurements
Experiments were performed in a 5 m 3 dual-reactor smog chamber, the details of which were described elsewhere . Briefly, the reactors were placed in a thermally isolated enclosure in which the temperature and RH could be well controlled by blowers and air conditioners. The variation range of temperature was 0.5 • C, and the limit for variability in RH was below 5 %. Experiments were performed using multiple UV (ultraviolet) light sources, which resulted in a sim-ilar spectrum as that of solar radiation .
During the experiments, the concentrations of NO, NO 2 and SO 2 were measured by corresponding gas analyzers (Teledyne API T200UP and Thermo Scientific 48i). The concentration of toluene was measured by proton transfer reaction quadrupole mass spectrometry in H 3 O + mode (PTR-QMS; Ionicon). The size distributions, number concentrations and mass concentrations of the SOA were determined by a scanning mobility particle sizer (SMPS), which consisted of an electrostatic classifier (EC; TSI 3080), a differential mobility analyzer (DMA; TSI 3081) and a condensation particle counter (CPC; TSI 3776). A density of 1.4 g cm −3 was used for the toluene SOA according to previous studies (Ng et al., 2007;Li et al., 2014).
The extinction coefficients at 532 nm were detected with custom-built cavity ring-down spectrometers (CRDSs), the details of which were given previously (Wang et al., 2012;K. Li et al., 2014. The scattering, absorption and extinction coefficients at 375 nm were measured by a photoacoustic extinctiometer (PAX-375, Droplet Measurement Technologies).

Experimental conditions and offline measurements
As shown in Table 1, four sets of experiments were conducted as follows: D, dry conditions (RH < 5 %; SO 2 < 1 ppb); DS, dry conditions with SO 2 (RH < 5 %; SO 2 = 30-50 ppb); W, wet conditions (RH > 80 %; SO 2 < 1 ppb); and WS, wet conditions with SO 2 (RH > 80 %; SO 2 = 30-50 ppb). To ensure that the results were reproducible, each experiment was performed twice. Before each experiment, the chamber was flushed at least three times using zero air (AADCO 737-15) to achieve clean conditions in which the number concentration of particles lower than 50 cm −3 and concentrations of NO x and SO 2 lower than 1 ppb. After preparation, a known volume of organic precursor was added to the chamber by a glass U tube through zero air. HONO was then added into the chamber by bubbling a small flow of zero air through the immediately mixed solution of 1 mL of 1 wt % NaNO 2 and 2 mL of 10 wt % H 2 SO 4 . The byproducts NO and NO 2 were also introduced into the chamber. The concentrations of toluene in the chamber were approximately 180 ppb, while the total NO x levels were approximately 500 ppb. Parenthetically, unlike NO x analyzers equipped with molybdenum converters, HONO was not detected in the T200UP NO x analyzer with a patented highefficiency photolytic converter. SO 2 was introduced into the chamber from a 29.6 ppm gas cylinder. For experiments under wet conditions, zero air was passed through a fuel cell humidifier (FC-200-780-7MP, Perma Pure), and the RH was above 80 % before the experiments started. It is important to highlight that dew drops did not appear in the chamber during the experiments due to accurate control of the chamber temperature (±0.5 • C). The contents in the chamber were mixed for approximately 10 min, and then the UV lamps were turned on to initiate the photoreaction. When the SOA mass concentration stopped rising, the lights were switched off to stop the reaction. All experiments were conducted at 25.0 ± 0.5 • C.
After each experiment, the SOA was collected by PTFE filters (0.2 µm, 47 mm, MerckMillipore, TYPE FGLP) and then extracted in 5 mL methanol in an ultrasonic bath for 30 min for mass spectrometry analysis and in 10 mL for UV-Vis absorption analysis. The methanol solutions were analyzed by electrospray ionization time-of-flight mass spectrometry (ESI-TOF-MS) in positive and negative mode. A typical mass-resolving power of > 30 000 was achieved at m/z 200, with an absolute mass error below 3 ppm. Each sample extract was measured in triplicate with an injection volume of 10 µL. The ions detected in the filter blank were subtracted, and molecular formulas were assigned to all ions in the samples having signal-to-noise ratios greater than 10. The settings for ESI positive mode were as follows: capillary voltage, 3 kV; nebulizer gas pressure, 0.4 bar; dry gas flow rate, 3 L per minute; dry gas temperature, 180 • C; and mass detection range, m/z 100-1000. For negative mode, the capillary voltage was set to 2.5 kV, the nebulizer gas pressure was 0.3 bar, the dry gas flow rate was 4 L per minute, the dry gas temperature was 150 • C, and the mass detection ranged from m/z 200 to 1000.

Calculation of the complex refractive indexes
The complex refractive index (RI; RI = n + ki) is the only intrinsic optical property of particles, including the real part of complex refractive index -RI(n), standing for scatteringand the imaginary part of complex refractive index -RI(k), standing for absorption (Bohren and Huffman, 1983;Bond and Bergstrom, 2006). The RI calculation method used in this study has been applied and approved in previous studies (K. Li et al., 2014Peng et al., 2018). Briefly, the extinction coefficients (α ext ) at 532 nm measured by CRDS and the scattering, absorption and extinction coefficients at 375 nm measured by PAX-375 were calculated on the basis of Mie theory. For monodispersed spherical particles, α ext can be represented by where σ ext is the extinction cross section, N is the number concentration of particles, D is the diameter of particles, and Q ext is the extinction efficiency, which is the ratio of the extinction cross section to the geometric area of the particles. For polydispersed particles with lognormal size distribution and a geometric standard deviation (σ g ) smaller than 1.5, α ext can be expressed as follows: where S tot is the total surface concentration. A hypothesis was made for simplifying data processing, as follows: the Q ext value of the polydispersed particles in the whole size distribution range was the same as the Q ext value of particles with the surface mean diameter (D sm ). The S tot and D sm values were given by SMPS directly, of which the uncertainties were ±1 % and ±5 %, respectively. The value of Q ext of particles with a given D sm was subsequently calculated through the Mie program. The best-fit complex index was determined by minimizing the following reduced merit function (χ r ):

Calculation of RI of products
The RI(n) of the products was calculated by the quantitative structure-property relationship, developed by Redmond and Thomson (Redmond and Thompson, 2011). In brief, the RI(n) values were estimated by polarizability (α), degree of unsaturation (µ), mass density (ρ m ) and molecular weight (M) via the application of Eq. (4).
The mass density of each molecule was estimated by the E-AIM model (Extended Aerosol Inorganics Model; http:// www.aim.env.uea.ac.uk/aim/density/density.php, last access: 22 October 2019), which had high accuracy for 166 organic compounds (Girolami, 1994). It should be noted that the model should be used at 550 nm. It would cause only small, acceptable errors for its use at 532 nm. This method was employed in our previous studies as well Peng et al., 2018). Similar to previous studies, negligible light absorption was found for toluene SOA at wavelength > 500 nm, as shown in Fig. S3 in the Supplement as well, so only the real part of RI is considered at 532 nm (J. K. Li et al., , 2014. For brevity, at 532 nm, "RI" represents the real part of the complex refractive index, RI(n), in the following text.
3 Results and discussion

General results of the experiments
As shown in Fig. 1, particles were generated only several minutes after the lamps turned on. Although the concentrations of NO differed under different conditions, the maximum concentrations of NO x were approximately the same, ∼ 500 ppb. The consumption of the organic precursor was similar in all experiments (∼ 100 ppb), indicating similar OH exposures and oxidation levels. As shown in Fig. S1, the maximum total number concentrations of the toluene SOA were almost doubled under conditions with SO 2 , which implied that the addition of SO 2 could promote new particle formation. At the end of the experiments, the S tot of particles were approximately 2.9 × 10 9 , 4.3 × 10 9 , 3.8 × 10 9 and 6.9 × 10 9 nm 2 cm −3 under the D, DS, W and WS conditions, respectively. The addition of SO 2 caused an enhancement of the total surface concentrations of the particles, probably because of the formation of new particles and/or acidiccatalyzed reactions, as implied in previous studies (Chu et al., 2016;Deng et al., 2017;Hallquist et al., 2009). Moreover, the S tot under wet conditions was higher than that under dry conditions as well, probably due to aqueous reactions of water soluble products, as indicated previously (Jia and Xu, 2018;Kamens et al., 2011;Hinks et al., 2018). The extinction efficiencies of the SOA under four different conditions at 375 and 532 nm are illustrated in Fig. 2. The fitting curves match well with the experiment data points under each condition. At 532 nm, the average values of RI(n) were 1.412, 1.348, 1.504 and 1.468 under the D, DS, W and WS conditions, respectively. The complex refractive indexes at 375 nm under the D, DS, W and WS conditions were 1.45, 1.37 + 0.014i, 1.4566 + 0.022i and 1.51 + 0.012i, re-  Table 1 for details. The concentrations of gas species are shown on the left axes, while the concentrations of particles are shown on the right axes.
spectively. The complex refractive indexes of the SOA as a function of surface mean diameter (D sm ) under various conditions are shown in Fig. S2. As the reaction proceeded, the surface mean diameter of particles increased, while RI(n) decreased. This phenomenon was observed in all experiments in this study and agreed well with previous studies (K. Li et al., 2014;Peng et al., 2018).
The complex refractive indexes under dry conditions (Figs. 2a and 3a) were comparable with those reported in previous studies (Nakayama et al., 2013;Kim and Paulson, 2013;Li et al., 2014;Moise et al., 2015). For example, Kim and Paulson (2013) discovered that the RI(n) values of toluene-derived SOA at 532 nm were in the range of 1.35 to 1.61 (experimental conditions: 3.2 ppm toluene and 3.1 ppm NO x ) (Kim et al., 2010;Kim and Paulson, 2013). Li et al. (2014) found RI(n) of toluene-derived SOA was approximately 1.5 under HONO conditions, and the concentration of NO x was approximately 1 ppm (Li et al., 2014). The slightly lower RI(n) in the current study (1.412) might be a result of the lower NO x level, as Nakayama et al. (2013) found that higher NO x levels would lead to higher RI(n) values (Nakayama et al., 2013).
As shown in Fig. S3, the absorption of toluene-derived SOA under the D condition at 375 nm was approximately 0.002 a.u., from which we could calculate that RI(k) under the D condition as approximately 0.001 (Sun et al., 2007). PAX-375 showed no absorption of the toluene-derived SOA under dry conditions. In other words, the absorption of SOA under dry conditions was below the detection limit, and the calculated RI(k) was approximately 0.0009. Under this circumstance, RI(k) was set to zero (±0.001), which was lower than other studies. Table S1 in the Supplement lists the complex refractive indexes of toluene-derived SOA. For instance, Nakayama et al. (2013) reported that RI(k) of toluene-derived SOA was 0.05 at 355 nm and 0.002-0.007 at 405 nm (experimental conditions: 4.0 ppm toluene and 0.1-0.6 ppm NO x ) (Nakayama et al., 2010(Nakayama et al., , 2013, and  suggested that the RI(k) of toluene-derived SOA was 0.01-0.03 at 320 nm and 0.002-0.02 at 405 nm (experimental conditions: 5.0 ppm toluene and 0-10 ppm NO x ) . This phenomenon may have occurred because the ratios of NO to NO 2 during the formation of SOA were high, and light-absorbing substances, such as nitrophenols, would be hardly produced in our system in that case.    As shown in Fig. 3, RI(n) and RI(k) at 532 and 375 nm under the four different conditions differ greatly. At 375 nm, RI(k) under the W condition was the highest, while RI(k) under the DS and WS conditions were the second highest, approximately 0.014. For RI(n), the highest value appeared in the W condition, while the lowest value showed in the DS condition. All values of RI(n) at 375 nm were larger than those at 532 nm with difference of 0.022-0.062. The real and imaginary parts of RI increased as the wavelength decreased, which was in good accordance with previous studies (Moise et al., 2015;Flores et al., 2014a, b;Lin et al., 2015;Liu et al., 2015;Nakayama et al., 2013Nakayama et al., , 2018Nakayama et al., , 2015. For example, Liu et al. (2015) found the values of RI(n) of toluene-derived secondary organic material decreased by approximately 0.1 when wavelength changed from 240 to 800 nm . In this situation, when concerning the optical properties of SOA, the wavelength should be considered as well.
To explain the RI variations in toluene-derived SOA under the different conditions, aerosols were collected and analyzed by ESI-TOF-MS both in positive and negative mode. Mass spectrometry has enabled the molecular characterization of organic oligomers in SOA because of its high resolution and mass accuracy (Lin et al., , 2012Hinks et al., 2018). The products were inferred based on mass spectra (MS) and previous studies (Kelly et al., 2010;Hinks et al., 2018;T. Liu et al., 2016;Nakayama et al., 2015). Most products were concentrated around m/z 200-400 (Figs. 4, S4, S5 and S6; Tables 2, S2 and S3), and some of them contained organonitrogen groups, which could be seen more clearly in negative mode. However, previous studies mainly focused on a lower mass range (< 200 Da) (Hinks et al., 2018;Forstner et al., 1997;Staudt et al., 2014;Wang et al., 2018;Birdsall et al., 2010); e.g., Birdsall et al. (2010) detected benzaldehyde, cresol, phenol and butenedial in the photooxidation process of toluene under both low-NO x and high-NO x conditions (Birdsall et al., 2010). In this work, we analyzed the products around m/z 200-400 and discussed the effect of SO 2 under different humidity. The representative identified molecular weight, molecular formula and calculated RI(n)  Tables S2 and S3.

Effects of high humidity
High relative humidity causes a substantial increase in the RI(n) values for the toluene-derived SOA (Figs. 2 and 3). The real part of complex refractive indexes increased from 1.45 to 1.566 at 375 nm and 1.412 to 1.504 at 532 nm in absence of SO 2 . Previous studies reported similar findings; e.g., K.  found that high RH would enhance the light scattering of SOA by multiphase reactions (K. ; J.  discovered that mass absorption coefficient values at 365 nm of SOA produced by photooxidation of trimethylbenzene and toluene were enhanced when reaction RH increased from < 5 % to 80 % (J. .
Wet conditions led to higher RI(n) values compared with those of dry conditions and induced the formation of oligomers, as shown in Figs. 5 and S4. Figure S4 and previous studies (K.  have shown that oligomers above 500 Da appear to have high relative intensities under wet conditions, and these products possess high RI(n) values, as shown in Tables 2, S2 and S3. Multiphase reactions in the presence of liquid water led to the formation of oligomers from intermediate products such as glyoxal and methylglyoxal, resulting in large enhancements in RI(n), which was proven in previous studies (K. . As was proven in the study of Molteni et al. (2018), highly oxygenated molecules (HOMs) were also important oxidation products of toluene. HOMs that contained a phenyl ring structure could not absorb light above 300 nm; however, the subsequent HOMs of nitration and ring-opening might absorb UV-Vis light. HOMs were detected in our experiments, especially under high-humidity conditions. Oligomers pro- duced from multiphase reactions, such as acid-catalyzed aldol condensations, were found in mass spectra from products derived under the W condition, especially in positive mode, which could extend light absorption to longer wavelengths (Lambe et al., 2013). Lambe et al. (2013) found that the conjugated double bonds could enhance the absorption in the UV-Vis light irradiation as well (Lambe et al., 2013). This phenomenon was also proved in previous reports (Jia and Xu, 2018;Faust et al., 2017;Liu et al., 2019).

Effects of SO 2 under dry conditions
Figures 2 and 3 show that adding SO 2 resulted in a decrease in the RI(n) values of toluene-derived SOA from 1.45 to 1.37 at 375 nm and from 1.412 to 1.348 at 532 nm under dry conditions. Compared with that under the D condition, RI(k) under the DS condition increased by 0.014. A similar phenomenon was also found in other systems; e.g., Nakayama et al. (2015) discovered that adding SO 2 caused the RI(n) values of isoprene-derived SOA to decrease, while causing the RI(k) values to increase (Nakayama et al., 2015). Nakayama et al. (2018) further found that adding SO 2 would cause significant light absorption at short visible and ultraviolet wavelengths of isoprene-derived SOA under OH oxidation progresses (Nakayama et al., 2018), which was quite similar to our results.
To clearly see the change in chemical compositions under the D and DS conditions, the subtraction plots of mass spectra are shown in Fig. S6. As we can see in Fig. S6a, adding SO 2 caused the relative intensities of products with lower molecular weight (< 200 Da) to increase. These products are mainly alcohols and esters with smaller molecular weight and RI(n) values of which are lower than 1.4. The reason for this is most likely that the addition of SO 2 caused large amounts of new particle formation and high particle number concentrations in the system (Chu et al., 2015(Chu et al., , 2016Deng et al., 2017;Liu et al., 2019), which could provide a larger mass concentration and adsorb highly volatile small molecules into the particle phase (K. . These small molecule products usually possess high volatility and a low oxidation state, which would also reduce the oxidation state of aerosols generated by toluene. Previous studies have found similar phenomena, for example, Zhao et al. (2018) found that high SO 2 concentrations decreased the ratio O/C in α-pinene or limonene systems (Zhao et al., 2018), while T.  discovered that the oxidation state of carbon was −0.51 ± 0.06 for SOA formed from light-duty gasoline vehicle exhaust with SO 2 and −0.19 ± 0.08 without SO 2 (T. , which all implied that adding SO 2 reduced the oxidation state of the resulting SOA. It should be noted that our offline analytical method caused the loss of a large part of alcohols and esters with small molecular weights, which might overestimate the values of their RI(n).
Light absorption properties of SOA are related to its composition, the contribution of each product to light absorption and so on Moise et al., 2015). Ji et al. (2017) reported that toluene oxidation proceeded dominantly via the cresol pathway and formed highly functionalized products such as polyhydroxytoluenes (Ji et al., 2017), which would undergo subsequent reactions with OH to form the precursors including α-carbonyl compounds, organic acids and other highly oxygenated low-volatility products. Although the peak of 3-methylcatechol, one of the major polyhydroxytoluenes of toluene oxidation products, appeared around 275 nm in the UV-Vis spectrum (Ferris et al., 1971), the subsequent products like nitrocatechol and low-volatility oligomers formed by small α-carbonyl compounds might absorb near-ultraviolet and visible light. Parenthetically, the reaction rate of 3-methylcatechol with OH radicals was 7.44 × 10 −11 cm 3 molecule −1 s −1 , and the estimation of atmospheric lifetime was 22 min approximately (Coeur-Tourneur et al., 2010). Small α-carbonyl compounds and organic acids might undergo particle phase reactions, e.g., acid-catalyzed aldol condensation reactions under SO 2 conditions, plausibly contributing to the observed light absorption (Fu et al., 2009(Fu et al., , 2008. Nakayama et al. (2018) reported the same phenomenon on isoprene SOA (Nakayama et al., 2018(Nakayama et al., , 2015, while Marrero-Ortiz et al. (2019) found BrC particles were formed from small α-dicardonyls and amines (Marrero-Ortiz et al., 2019). Formation of charge transfer (CT) complexes might be another reason for light absorption enhancement, which could lead to optical transitions through a transfer of charge from a donor group, such as hydroxyls, to an acceptor group, such as a ketone or aldehyde (Phillips and Smith, 2014). For organosulfate, another kind of BrC, we did not detect them under the DS condition, which is in accordance with previous studies (Staudt et al., 2014).

Effects of SO 2 under wet conditions
Under wet conditions with SO 2 , the average RI(n) values of toluene-derived SOA were 1.51 at 375 nm and 1.468 at 532 nm, which were higher than those of the DS condition and lower than those of the W condition (Figs. 2 and 3). As for absorption, RI(k) under the WS condition was lower than the W condition and similar to that under the DS condition. Figures 4 and S5 show the results of mass spectra difference in toluene-derived SOA under the DS or W condition minus the WS condition. Low-oxidation state organic matters and oligomers were both found in the mass spectra. Under the WS condition, relative intensities of products above 400 Da (oligomers) were higher than those under the DS condition (negative values), suggesting more types of oligomers were produced. More types of low organic matters were observed compared to those under wet conditions without SO 2 . Under these combined effects, the values of RI(n) of toluenederived SOA of the WS condition were lower than those under the W condition and higher than those under the DS condition.
The RI(k) under the WS condition is almost equal to that under the DS condition, while lower than the values under the W condition. Products with lower oxidation state and less conjugated oligomers caused by addition of SO 2 might be the reason for this phenomenon, as proved in previous studies (Nakayama et al., 2015;. The concentrations of the donor group under the WS condition were lower than those under the DS condition, resulting in lower concentrations of CT complexes, which might be another reason for reducing the values of RI(k). Organosulfate compounds were not found under the WS condition. The combined effect of SO 2 and wet conditions on optical properties of toluene-derived SOA is the first described, to the best of our knowledge, and has a significant influence on light absorption, extinction, visibility and direct radiative forcing of regional air, especially in complex polluted areas. Our results provide some explanations for the observed variation, and further research is needed to quantify this synergistic effect.

Atmospheric and climate implication
The values of mass cross section (MAC) were calculated for the DS, W and WS conditions at 375 nm, and the method was described in the Supplement. The SOA under the D condition was not calculated because no absorption was observed. The average values of MAC were 0.2749, 0.3082 and 0.2131 g m −2 under the DS, W and WS conditions, respectively. These results are similar to Liu's work, which noted that wet conditions would cause higher values of MAC than those observed under dry conditions (J. . However, the MAC values in Liu's work were approximately 0.4 g m −2 under wet conditions and 0.01 g m −2 under dry conditions at 380 nm, higher than our results, which is likely due to the high ratios of NO to NO 2 in their study.
The impacts of the atmosphere and climate were assessed by comparing the ratio of light extinction efficiency and simple forcing efficiency (SFE) under the four different conditions, and the results are shown in Figs. 6 and S7. Aerosol sizes between 100 and 250 nm were circled because these are atmospherically relevant aerosols sizes Tao et al., 2017). As shown in Fig. 6, adding SO 2 caused a light extinction efficiency reduction of approximately 16 %-35 %, with an average of 25 % at 532 nm, and humid conditions enhanced the light extinction efficiency by approximately 36 %-64 %, with an average of 50 %. For the comprehensive impact made by SO 2 and high humidity, the light extinction efficiency increased about 16 %-47 % with an average of 30 %. These results confirm that SOA generated under synergistic pollution conditions -conditions that contained high emissions of primary particles and gaseous precursors from multiple sources, efficient secondary matter formation, as well as adverse meteorological and climate conditions and regional transport -might have a greater impact on the visibility reduction, atmospheric photochemical reactions and secondary species formation (Dickerson et al., 1997;. We estimate the clear-sky direct radiative forcing per unit optical depth with the help of the SFE concept, as mentioned in the Supplement. Figure S7 shows the SFE of toluenederived SOA under the D, DS, W and WS conditions at 375 and 532 nm, which could also reflect the change in direct radiative forcing (DRF). The forcing efficiency crossed over from negative (warming) to positive (cooling) values at diameter ≈ 200 nm under the DS, W and WS conditions. The lowest SFE values, approximately −35 W g −1 , were observed under the W condition, while at diameter approximately 200 nm, the SFE values were −10 and −20 W g −1 under the WS and DS conditions, respectively. All values of SFE at 532 nm were above zero, which indicates that toluenederived SOA has cooling effects on climate change. The SFE values under the DS condition decreased by approximately 25 % compared to that under the D condition, and the SFE values under the W and WS conditions increased by approximately 50 % and 30 % compared to that under the D condition, respectively, which is a similar trend to that of the extinction efficiency. The SFE values in our system coincided with those of organic aerosols from wood combustion and the burning of boreal peatlands (Chen and Bond, 2010;Chakrabarty et al., 2016). Our model does not include hygroscopicity and other factors, which would increase particle size and negative forcing. In this situation, our forcing is much lower than those in global climate models (Schulz et al., 2006). The combined effects of SO 2 and humidity should be considered in the modified climate model.

Conclusion
The effect of SO 2 under different humidities on the optical properties of SOA photooxidized by toluene was investigated in this study. The results show that for the experimental system with RH greater than 80 %, as expected, the increase in humidity greatly enhanced the real part of RI, from 1.412 to 1.504 at 532 nm, and from 1.45 to 1.566 at 375 nm, the imaginary part of RI was enhanced as well, which is probably because of the oligomers formation from multiphase reactions. Adding SO 2 can reduce the RI(n) values of toluenederived SOA at 375 and 532 nm, whether under low or high humidity. The RI values of toluene-derived SOA produced under dry conditions with SO 2 are 1.37 + 0.014i at 375 nm and 1.348 at 532 nm, while the RI(n) values under dry and SO 2 -free conditions are 1.412 at 532 nm and 1.45 at 375 nm.
The reason for this phenomenon might be that adding SO 2 caused large amounts of new particle formation and high particle mass concentrations in the system, which could have absorbed highly volatile molecules into the particle phase. Highly volatile molecules produce lower oxidation state and lower RI(n) values, resulting in the decrease in RI(n) for toluene-derived SOA. The increase in RI(k) is probably related to acid-catalyzed reactions on acidic particles. For the experimental system under high-humidity conditions with sulfur dioxide, the RI(n) was higher than SOA derived from dry conditions without SO 2 . The RI(k) under this condition are lower than those under wet conditions without SO 2 because fewer oligomers formed. The extinction properties under the WS condition are approximately 30 % higher than under the D condition. The results here highlighted that the combined effect of SO 2 and high humidity could greatly enhance the refractive index, light scattering, and direct radiative forcing of toluene-derived SOA and, potentially, others. These results will improve our understanding of SOA optical properties, especially under complex atmospheric conditions.
Data availability. All data in this work are available by contacting the corresponding author Maofa Ge (gemaofa@iccas.ac.cn) or Weigang Wang (wangwg@iccas.ac.cn).
Author contributions. WW, MG and WZ conceived and led the studies. WZ, JL and CP performed chamber simulation and data analysis. KL, LZ, BS, YC and ML discussed the results and commented on paper. WZ prepared the paper with contributions from all coauthors.
Competing interests. The authors declare that they have no conflict of interest.