Photochemical degradation of iron(III)-citrate/citric acid aerosol quantified with the combination of three complementary experimental techniques and a kinetic process model

Abstract. Iron(III) carboxylate photochemistry plays an important role in aerosol aging, especially in the lower troposphere. These complexes can absorb light over a broad wavelength range, inducing the reduction of iron(III) and the oxidation of carboxylate ligands. In the presence of O2, ensuing radical chemistry leads to further decarboxylation, and the production of ·OH, HO2·, peroxides, and oxygenated volatile organic compounds, contributing to particle mass loss. The ·OH, HO2·, and peroxides in turn re-oxidize iron(II) back to iron(III), closing a photocatalytic cycle. This cycle is repeated resulting in continual mass loss due to the release of CO2 and other volatile compounds. In a cold and/or dry atmosphere, organic aerosol particles tend to attain highly viscous states. While the impact of reduced mobility of aerosol constituents on dark chemical reactions has received substantial attention, studies on the effect of high viscosity on photochemical processes are scarce. Here, we choose iron(III)-citrate (FeIII(Cit)) as a model light absorbing iron carboxylate complex that induces citric acid (CA) degradation to investigate how transport limitations influence photochemical processes. Three complementary experimental approaches were used to investigate kinetic transport limitations. The mass loss of single, levitated particles was measured with an electrodynamic balance, the oxidation state of deposited particles was measured with X-ray spectromicroscopy, and HO2· radical production and release into the gas phase was observed in coated wall flow tube experiments. To quantitatively compare these experiments and determine important physical and chemical parameters, a numerical multi-layered photochemical reaction and diffusion (PRAD) model that treats chemical reactions and transport of various species was developed. We observed significant photochemical degradation, with up to 80 % mass loss within 24 hours of light exposure. Interestingly, we also observed that mass loss always accelerated during irradiation, resulting in an increase of the mass loss rate by about a factor of 10. When we increased relative humidity, the observed particle mass loss rate also increased. This is consistent with strong kinetic transport limitations for highly viscous particles. The PRAD model was tuned to reproduce all experimental results and captured the essential chemistry and transport during irradiation. In particular, the photolysis rate of FeIII, the re-oxidation rate of FeII, HO2· production, and the diffusivity of O2 in aqueous FeIII(Cit)/CA system as function of relative humidity and FeIII(Cit) / CA molar ratio could be constrained. Photochemical degradation under atmospheric conditions predicted by the PRAD model shows that release of CO2 and re-partitioning of organic compounds to the gas phase may be very significant to accurately predict organic aerosol aging processes.


Abstract. Iron(III) carboxylate photochemistry plays an important role in aerosol aging, especially in the lower troposphere. These complexes can absorb light over a broad wavelength range, inducing the reduction of iron(III) and the oxidation of carboxylate ligands. In the presence of O 2 , the ensuing radical chemistry leads to further decarboxylation, and the production of . OH, HO .
2 , peroxides, and oxygenated volatile organic compounds, contributing to particle mass loss. The . OH, HO .
2 , and peroxides in turn reoxidize iron(II) back to iron(III), closing a photocatalytic cycle. This cycle is repeated, resulting in continual mass loss due to the release of CO 2 and other volatile compounds. In a cold and/or dry atmosphere, organic aerosol particles tend to attain highly viscous states. While the impact of reduced mobility of aerosol constituents on dark chemical reactions has received substantial attention, studies on the effect of high viscosity on photochemical processes are scarce. Here, we choose iron(III) citrate (Fe III (Cit)) as a model light-absorbing iron carboxylate complex that induces citric acid (CA) degradation to investigate how transport limitations influence photochemical processes. Three complementary experimental approaches were used to investigate kinetic transport limitations. The mass loss of single, levitated particles was measured with an elec-trodynamic balance, the oxidation state of deposited particles was measured with X-ray spectromicroscopy, and HO .
2 radical production and release into the gas phase was observed in coated-wall flow-tube experiments. We observed significant photochemical degradation with up to 80 % mass loss within 24 h of light exposure. Interestingly, we also observed that mass loss always accelerated during irradiation, resulting in an increase of the mass loss rate by about a factor of 10. When we increased relative humidity (RH), the observed particle mass loss rate also increased. This is consistent with strong kinetic transport limitations for highly viscous particles. To quantitatively compare these experiments and determine important physical and chemical parameters, a numerical multilayered photochemical reaction and diffusion (PRAD) model was developed that treats chemical reactions and the transport of various species. The PRAD model was tuned to simultaneously reproduce all experimental results as closely as possible and captured the essential chemistry and transport during irradiation. In particular, the photolysis rate of Fe III , the reoxidation rate of Fe II , HO .
2 production, and the diffusivity of O 2 in aqueous Fe III (Cit) / CA system as function of RH and Fe III (Cit) / CA molar ratio could be constrained. This led to satisfactory agreement within model

Introduction
Photochemistry in the atmosphere (either in the gas phase or in the particle phase) plays an important role in aerosol aging processes. Photochemically produced free radicals in the gas phase (mainly . OH) can be taken up by aerosol particles, inducing multiphase chemistry. However, uptake is limited by the collision rate and condensed phase molecular transport when diffusion coefficients are sufficiently low, which restricts chemical reactions to the near-surface region of the particle. In contrast, photochemically generated radicals in aerosol particles can be present throughout their bulk due to light penetrating their whole volume. Radical production induced by a direct photochemical reaction occurs when the energy of light quanta is high enough (mostly the UV part of the solar spectrum) to cause bond cleavage or rearrangement in a molecule. In the lower troposphere where UV light intensity is low, indirect photochemistry initiated by near UV and visible radiation may become significant. Important indirect photochemical processes are transition metal complex photochemistry and photosensitized processes (Corral Arroyo et al., 2018;George et al., 2015). This work focuses on iron carboxylate-catalyzed photochemistry due to its abundance and reactivity in the atmosphere.
Iron is the most abundant transition metal in the earth's crust. Wind erosion is the main source of iron in the atmosphere, but anthropogenic activities such as industrial processes, traffic, and combustion processes can also potentially release iron in particulate form (Deguillaume et al., 2005). Depending on parameters such as temperature, pH value, ionic strength, and the concentration of involved substances, iron can combine with inorganic or organic ligands to form complexes (Deguillaume et al., 2005;Faust and Hoigné, 1990;Kieber et al., 2005). Iron can be found complexed with low molecular weight inorganic species such as the hydroxide anion (OH − ), sulfate (SO 2− 4 ), and sulfite (SO 2− 3 ) (Brandt and van Eldik, 1995;Hofmann et al., 1991;Weschler et al., 1986). Quantifying iron atmospheric processing and solubility is of global importance, especially for nutrient input into the world's oceans (Hamilton et al., 2019;Kanakidou et al., 2018). Heterogeneous chemistry involving particulate iron and SO 2 can result in sulfate formation and increase aerosol loading (Grgić et al., 1998(Grgić et al., , 1999Grgić, 2009). Additionally, iron photochemical processing in aerosol particles, fog droplets, and cloud water is an important radical source (Bianco et al., 2020;Abida et al., 2012) and sink for or-ganic compounds (Weller et al., 2014(Weller et al., , 2013Herrmann et al., 2015). Organic compounds are a major component in atmospheric aerosol particles and have received more and more attention as potential ligands for iron(III) complexation. For instance, humic-like substances derived from water-soluble organic compounds have been reported to be strong chelating ligands with iron(III) (Dou et al., 2015;Kieber et al., 2003;Okochi and Brimblecombe, 2002;Willey et al., 2000). Oxalate and other carboxylates have been identified to be important ligands for iron(III) because they are available in sufficient amounts (Chebbi and Carlier, 1996;Kahnt et al., 2014;Kawamura et al., 1985) and the carboxylate groups are acidic enough to dissociate and chelate with iron(III) at atmospheric pH values (Okochi and Brimblecombe, 2002). Field studies have confirmed that soluble iron is mostly in complexes with carboxylate functions (Tapparo et al., 2020;Tao and Murphy, 2019).
In atmospheric aqueous phases, iron normally exists in oxidation states (II) and (III) and they can convert into each other via redox cycling. The ratio between iron(II) and iron(III) in aerosol particles is quite variable and depends on several factors such as the presence of light, oxidizing compounds, and ligands. For example, Grgić et al. (1999) reported that the concentration ratio of iron(II) to iron(III) varied between 0.9 and 3.1 in urban aerosol particles with a size range of 0.4-1.6 µm. Iron(III) carboxylate complexes [Fe III (OOC−R)] 2+ are well-known photoactive compounds Weller et al., 2013Weller et al., , 2014. They can easily get excited by light in the UV-VIS range, inducing ligand-tometal charge transfer (LMCT) (Cieśla et al., 2004), which is an inner-sphere electron transfer (i.e., the electron transfer occurs via a covalently bound bridging ligand) from the carboxylate group to the iron. Investigations using time-resolved transient spectroscopy reported the formation of long-lived radical complexes, [Fe II ( . OOC−R)] 2+ , with lifetimes on the order of a millisecond, followed by the dissociation to the organic radical R−COO . and an Fe II aqua complex (Feng et al., 2007;Glebov et al., 2011;Pozdnyakov et al., 2009;Zhang et al., 2009): R−COO . will decarboxylate almost instantaneously (k R3 ≈ 10 9 -10 12 s −1 ) (Abel et al., 2003;Bockman et al., 1997;Hilborn and Pincock, 1991): The alkyl radical R . will react rapidly with dissolved O 2 , producing a peroxy radical with k R4 ≈ 2 × 10 9 M −1 s −1 (von Sonntag and Schuchmann, 1991): Subsequent reactions of R . and RO . 2 are specific depending on the type of ligand and its substitution.
In this work we investigated iron(III) citrate ([Fe III (OOCCH 2 ) 2 C(OH)(COO)], in short Fe III (Cit)), as a model species to better understand iron carboxylate photochemistry in atmospheric aerosol particles. Our Fe III (Cit) system undergoes LMCT reactions in the same way as countless other iron(III) carboxylate compounds (Cieśla et al., 2004;Weller et al., 2013Weller et al., , 2014. Its photochemical reaction scheme is well established in both solution (Abida et al., 2012;Faust and Zepp, 1993;Pozdnyakov et al., 2012) and solid states (Abrahamson et al., 1994). Citric acid (CA) is an established proxy for oxygenated atmospheric organic matter with its thermodynamic properties, water diffusivity, and viscosity having been well studied (Lienhard et al., 2012(Lienhard et al., , 2014Song et al., 2016). For these reasons, it is a valid and reliable proxy for atmospheric iron carboxylate photochemical processes.
As schematically described in Fig. 1, Fe III (Cit) absorbs light up to 500 nm, inducing LMCT, followed by immediate decarboxylation of the central carboxyl group since the hydroxyl group adjacent to a carboxyl group facilitates decarboxylation (Weller et al., 2013). In the presence of O 2 , oxidants such as HO .
2 and H 2 O 2 will be produced, which can oxidize Fe II back to Fe III via Fenton reactions (Fenton, 1894) with additional oxidant production. Fe III then combines with another CA molecule in this aqueous system, closing the photocatalytic cycle, in which iron acts as a catalyst for CA degradation. In addition, the generation of reactive oxygen species (ROS) and peroxy radicals leads to further decarboxylation and more production of oxygenated volatile organic compounds (OVOCs) (e.g., acetone) (Pozdnyakov et al., 2008;Wang et al., 2012). Therefore, this photodegradation process is potentially an important sink of carboxylate groups in the troposphere.
We expect that the chemistry subsequent to the initial photochemical reaction steps in the aerosol phase will be significantly altered by diffusion limitations when Fe III (Cit) particles mixed with CA attain a high viscosity. As viscosity increases, molecular diffusion coefficients tend to decrease  and therefore photochemical cycling will also be slow. Increasing water content is expected when the relative humidity (RH) increases and will effectively plasticize particles  leading to better mixing of all compounds within the particles and faster photochemical cycling when compared with lower RH. However, these effects have been investigated in dark systems but not in photochemical systems (Berkemeier et al., 2016;Shiraiwa et al., 2011;Shiraiwa and Seinfeld, 2012;Steimer et al., 2015a). In order to better understand this system and how it reacts to RH, we used a triad of photochemical experiments including electrodynamic balance (EDB) measurements, scanning transmission X-ray microscopy coupled with near-edge X-ray absorption fine structure (STXM/NEXAFS) spectroscopy, and a coated-wall flow-tube (CWFT) experiments Figure 1. Photocatalytic cycle of Fe III (Cit) complex in an aqueous particle containing CA with explicit charge balance given in Table 2. R corresponds to the carboxylate side chain −CH 2 COO − . to investigate how particle size, mass, and indicators of chemical composition change during photochemical processes. In this work we mostly focus on the humidity dependence of this photochemical degradation, while Alpert et al. (2021) focus on the impacts on ROS and the fate of free radicals during this photochemical degradation.
To perform a quantitative comparison of these experiments and determine their relevant properties, a numerical multilayered photochemical reaction and diffusion (PRAD) model was developed that treats chemical reactions and the transport of various species. In addition, we will use the PRAD model to simulate photochemical aging processes under atmospheric conditions. In the following we briefly discuss our experimental approaches in Sect. 2 and include a detailed explanation of the PRAD model in Sect. 2.5. A comparison between the experimental results and the PRAD model is presented in Sect. 3. Finally, we discuss the impact and atmospheric importance of kinetic limitations to photochemical degradation in Sect. 4.

Solution preparation
CA (≥ 99.5 %) and iron(III) citrate tribasic monohydrate (18 %-20 % Fe basis) were purchased from Sigma-Aldrich. Iron(II) citrate (Fe II (HCit)) was purchased from Dr. Paul Lohmann GmbH KG. Dilute aqueous solutions of Fe III (Cit) / CA and Fe II (HCit) / CA were made in ultrapure water (18 M cm −1 , Milli-Q). Since Fe III (Cit) only dissolves slowly in water, CA solutions with suspended Fe III (Cit) crystals were sonicated for at least 24 h. The same dissolving procedure was also applied to the Fe II (HCit) powders. Note that all procedures were done under red light illumination because Fe III (Cit) is light sensitive. The molar ratio between Fe III (Cit) and CA was different for each experimental method used in this study. For EDB, STXM/NEXAFS, and CWFT experiments, stock solutions were prepared with molar ratios of 0.05, 1.0, and 0.07, respectively.

Bulk property measurements by EDB
We used an EDB to measure the mass loss in single, levitated particles under irradiation. The experimental setup has been described previously (Steimer et al., 2015b). In short, an electrically charged aqueous particle (radius ∼ 10 µm) is injected into an EDB. The balance is of the double ring design (Davis et al., 1990) with a high AC voltage applied to the two parallel electrode rings and a DC voltage across the hyperbolic endcaps. The DC field compensates the gravitational force of the particle and is used as a measure for the mass of the particle. The EDB is placed in a three-wall glass chamber with a cooling liquid (ethanol) pumped through the two inner walls and an insulation vacuum between two outer walls to control the temperature (T ) at the location where the particle levitates. The RH within the chamber is regulated by adjusting the ratio of a dry and a humidified gas flow through the chamber. In the experiments described in this work, we used a typical total flow of 40 sccm and set the total pressure inside the cell at 8 × 10 4 Pa.
The spherical particles were characterized by two methods based on Mie-resonance spectroscopy: (i) a narrow bandwidth tunable diode laser (TDL, tuning range 765-781 nm) was used to determine the refractive index and radius simultaneously with high-precision (Steimer et al., 2015b); and (ii) simultaneously, a broad-band LED centered around 640 nm was used to illuminate the particle. The backscatter signal from the LED was recorded using a spectrograph with a slow-scan back-illuminated charge-coupled device (CCD) array detector to follow the resonance wavelength shift of the particle (Zardini et al., 2006). If we assume the refractive index stays constant during the experiment, the radius change of a particle is easily calculated from this resonance wavelength shift: We illuminated particles to induce a photochemical reaction with either a CW diode laser emitting at 375 nm (LuxX 375-20, Omicron Laserage) or a frequency-doubled diode laser emitting at 473 nm (gem 473, Laser Quantum). At wavelengths of 375 and 473 nm, Fe III (Cit) is reported to have a molar absorptivity of 796 and 60.7 M −1 cm −1 , respectively (Pozdnyakov et al., 2008). with the initial value of the particle radius determined using highresolution Mie-resonance spectroscopy (not shown), and the change in radius determined from low-resolution Mie-resonance data given in (b). (d) Normalized mass remaining ratio deduced from DC voltage (black), and from Mie-resonance shift (red) -assuming constant density. Yellow shaded region in panels (a), (c), and (d) indicate the time of laser irradiation with 0.25 W cm −2 at 375 nm.
In a typical EDB experiment, we let the particle equilibrate to RH and T in a pure O 2 gas phase for up to 10 h in the dark before irradiation. Exemplary raw data of an experiment at 46 % RH and 293.5 K are shown in Fig. 2. The measured DC voltage compensating for the gravitational force, as well as the radius of the particle deduced from Mie-resonance spectroscopy, decreased dramatically during illumination in the first 18 h, with more than half of the initial mass lost to the gas phase. Note that the radius and mass loss rates increased as seen in Fig. 2c. We assumed the refractive index and density of the particle did not change upon photochemistry and therefore the mass loss calculated from the DC voltage could be directly compared with size change by calculating the par-ticle mass remaining ratio, where m 0 is the particle mass prior to irradiation. Mass loss derived from both ways independently is shown in Fig. 2d and reveal that there is a little difference between the mass loss up to t ≈ 65 000 s, corresponding to when m(t) m 0 < 0.4. Therefore, the refractive index and density are mostly governed by those of aqueous CA up until half the particle mass is lost. The total mass loss over 24 h irradiation is more significant and drops by 80 % for the particular experiment shown in Fig. 2. In addition, we observed the mass loss rate was initially ∼ 1.3 % h −1 and increased to ∼ 14 % h −1 when 40 % to 60 % of the initial mass was lost. This mass loss acceleration is discussed in further detail with the help of the PRAD model simulations in the Sect. 3.1. At t ≈ 65 000 s, the mass loss slowed down considerably when we observed a distortion in the Mie-resonance pattern ( Fig. 2b and video in the Supplement). The distortion may be attributable to partial crystallization of iron citrate in the particle, which would explain the slowing photochemical degradation.

Chemical characterization by STXM/NEXAFS
STXM/NEXAFS measurements were performed at the Pol-Lux endstation located at the Swiss Light Source (SLS) to obtain the Fe oxidation state of particles between 0.2 and 2 µm in diameter (Flechsig et al., 2007;Frommherz et al., 2010;Raabe et al., 2008). Particles containing Fe III (Cit) / CA were nebulized from aqueous solution with a mole ratio between Fe III (Cit) : CA of 1 : 1. They were dried in air at RH < 30 % and impacted onto silicon nitride membranes mounted in portable sample holders. The sample holders were transported to the endstation in an evacuated container and shielded against ambient light. Once there, they were mounted in the PolLux environmental microreactor (Huthwelker et al., 2010) and kept under a total pressure of 1.5 × 10 4 Pa, T = 293.5 K, RH = 40, 50 or 60 %, with a controlled gas flow. Further details of sample preparation have been previously reported Huthwelker et al., 2010;Steimer et al., 2014). The microreactor was mounted in a vacuum chamber for in situ STXM/NEXAFS analysis. When desired, the microreactor could operate in vacuum conditions without a gas flow. The transmission of X-ray photons through the particles was measured and converted to optical density, OD = − ln(I /I 0 ), where I and I 0 are the transmitted and incident photon flux as a function of X-ray energy. The Fe L-edge absorption was probed over the X-ray energy range of 700-735 eV. Figure 3 shows an example of NEXAFS spectra of Fe III (Cit) / CA particles before (orange) and after (red) irradiation with UV light. X-ray energy calibration was consistently performed using FeCl 2 and compared with the existing literature for FeCl 2 and FeCl 3 salts  and a mixture of Figure 3. Iron L-edge NEXAFS spectra of Fe III (Cit) / CA particles before and after irradiation with UV light shown as orange and red, respectively. The previously recorded spectrum from mixed xanthan gum (XG) and FeCl 2 particles exposed to ozone is shown as the purple line, and a spectrum from FeCl 2 particles is shown as the blue line . FeCl 2 and FeCl 3 spectra from Moffet et al. (2012) are shown as the blue and green shading, respectively. The vertical dashed lines indicate that peak X-ray absorption at 707.9 eV shifted to 708.3 eV for Fe II and 709.6 eV shifted to 710.0 eV for Fe III . xanthan gum and FeCl 2 oxidized by O 3 . The peak absorption for iron(II) and iron(III) are at X-ray energies of 708.3 and 710.0 eV. We were capable of resolving peaks separated by 0.4 eV at the Fe L-edge. Ferrous and ferric iron peaks are separated by 1.7 eV and thus clearly distinguishable. Following a previous procedure , we imaged particles at these two energies to determine the OD ratio between them. The fraction of Fe III out of total Fe, β, was then determined using the parameterization from Moffet et al. (2012). It is important to note that the Xray energy absorption peaks observed for FeCl 2 and FeCl 3 were identical for FeCl 2 mixed with xanthan gum either unexposed or exposed to O 3 . However, we have found that these peaks shifted by about +0.4 eV, possibly due to the strong complexation with CA. Small shifts in energy can occur depending on the chemical environment surrounding Fe atoms (Garvie et al., 1994;Moffet et al., 2012). In agreement with Alpert et al. (2019), the peak absorption energies for our particles were independent of RH from 0 % to 60 % within ±0.2 eV. When calculating β, we always imaged particles at 708.3 and 710.0 eV. 2 production determined by CWFT The HO . 2 release upon irradiation of a Fe III (Cit) / CA thin film was measured by scavenging HO .
2 with an excess of NO in a CWFT reactor (Duran glass, inside diameter 1.2 cm, length 50 cm). The film was composed of Fe III (Cit) / CA, and deposited inside the tubular glass flow tube with a thickness between 0.15 and 0.2 µm and an error of about 20 %. Details of the film preparation have been described previously (Corral Arroyo et al., 2018;González Palacios et al., 2016). Seven UV lamps (UV-A range, Philips Cleo Effect) were mounted surrounding the glass reactor held at 298.15 K. The total light output between 300 and 590 nm was 210 W m −2 . The flows of N 2 , O 2 , and NO passing through the reactor were controlled. The NO concentration during CWFT photochemical experiments was in excess (> 10 13 molecules cm −3 ) to efficiently scavenge 99 % of HO .
2 produced by the film. The concentration of NO was tracked by a chemiluminescence detector (Ecophysics CLD 77 AM). In an example CWFT experiment at RH = 29.3 %, a clear NO loss was observed when the UV lights were switched on as shown in Fig. 4, which was due to the release of HO .
2 radicals into the gas phase and reaction of NO with HO .
2 forming NO 2 and OH . . OH . is then scavenged by NO producing HONO. The production of HONO was routinely checked as described in González Palacios et al. (2016). Therefore, the production rate of HO .
2 , P HO . 2 , was calculated from the loss rate of NO assuming a 2 : 1 ratio to HO .
2 conversion: where [NO] is the loss of gas-phase concentration of NO in molecules cm −3 , flow is the volumetric gas flow in the CWFT in cm 3 s −1 , and S film is the surface area of the film in cm 2 .

Development of the photochemical reaction and diffusion (PRAD) model
We developed a PRAD model to interpret our experiments and to understand any feedback between transport limitations and photochemistry, especially under low RH conditions, corresponding to high viscosity of the particle phase. The PRAD model consists of two modules: a detailed chemical process module, treating equilibria and chemical reactions, and a transport module handling the physical transport of all species (including diffusion in the aqueous phase as well as gas-particle phase partitioning). Conceptually, the PRAD model relies on the kinetic model framework for aerosol surface chemistry and gas-particle interactions (Pöschl et al., 2007), which is similar, for example, to the KM-GAP model . Numerically, the PRAD model uses a Euler forward step method as explained in detail below, while KM-GAP solves coupled differential equations. In passing, there are alternative approaches; for example, Kinetiscope (Houle et al., 2015) does not integrate sets of coupled differential equations to predict the time history of a chemical system. Instead, it uses a general stochastic algorithm to propagate a reaction. As illustrated in Fig. 5, the PRAD model divides a spherical droplet into a number of shells, n, which exchange molecules after each chemical time step. Shell thickness and the number of shells were adjusted to enable the resolution of steep concentration gradients within a reasonable computation time. The volume of each shell was constant instead of the thickness and the shells become thinner and thinner from the center to the surface of the particle. For each shell and at each time step, we first calculated the composition using the thermodynamic equilibria of the Fe III (Cit) / CA system, as listed in Table 2. The Newton-Raphson method (Burden and Faires, 2011) was then used to calculate the turnover and the concentration of products and reactants over time for the chemical reactions also listed in Table 2 with a fixed time step of 0.2 s. After each time step, diffusion of all species between the shells and the evaporation of products (or condensation of the reactant O 2 ) were computed. The time step, t, for physical transport processes was determined dynamically to ensure both numerical stability and computational efficiency.
For each species, the molar flux from shell i to the next shell i + 1 was calculated as where D l is the liquid phase diffusion coefficient of the corresponding species. Shell i extends from r i−1 to r i , while shell i + 1 extends from r i to r i+1 , with r being the distance from the particle center. In Eq. (4), c is the molar concentration in Figure 5. Schematic illustration of the PRAD model showing the shells, transport fluxes (green arrows), and chemical processes (red arrows) of each species. The particle is radially symmetric with the surface of the particle marked as r n .
each shell of the aqueous particle, defined as where N i is moles of a particular species in shell i, and V i is the total volume of shell i. At the outermost shell n, the gas-particle phase partitioning of each species was determined by the modified Raoult's law. The flux from shell n into the gas phase, f n , was calculated to be where D g is the gas phase diffusion coefficient of the corresponding species, and R is the gas constant. p partial is the partial pressure of the species in the gas phase at infinite distance from the particle. Based on Henry's law, vapor pressure of the species p vapor is defined as where c n is the molar concentration in the outermost shell n, and H is the Henry's law coefficient. Substituting Eq. (7) into Eq. (6) leads to We calculated the partial pressure of H 2 O from RH, and took the partial pressure of O 2 based on the total pressure and the fraction of O 2 gas flow used in experiments. For other species, p partial is negligible and assumed to be zero. So, for all species other than H 2 O and O 2 , f n simplifies to Based on the calculated fluxes, the change in moles, N i , of each species in t was given by from which the concentration and corresponding shell thickness using the molar volume of each species were recalculated for the next time step. The chemical loss rate of O 2 was treated in the liquid phase diffusion module instead of in the chemical module, since the loss can be very fast and its lifetime can be smaller than 10 −6 s. If O 2 loss due to reaction was determined in the chemical module, a chemical time step of 10 −6 s would result in extensive computational time beyond what is reasonable for this study. Therefore, within each liquid phase diffusion time step, in addition to the physical transport, the chemical loss of O 2 molecules was calculated in each shell: Integration yields where k O 2 is the total chemical reaction rate of O 2 . In the outermost shell n, O 2 molar concentration is assumed to be at steady state all times, meaning that the chemical loss of O 2 is always compensated by the condensation of O 2 from the gas phase: Substituting Eq. (5) in Eq. (13), we calculate the moles of O 2 in this shell, N n,O 2 , as where H O 2 is Henry's law coefficient of O 2 . All iron-containing compounds and all free ions are assumed not to evaporate from the particle phase. All other species have Henry's law coefficients given in Table 1. If the evaporation rate of some species is fast enough, then their concentration in the outermost shell n can always be at steady state, which means that is, From Eq. (16), c n can be deduced to be c n = c n−1 1 1 + 0.5 (r n − r n−2 ) r n D g / H RT r 2 n−1 D l .
And the number of moles of the volatile species in the outermost shell N n as N n = c n V n .
We have parameterized aqueous and gas phase diffusion coefficients, D j l and D j g , respectively, for all species j given in Appendix A1 and A2. There are several other assumptions and approximations made to adapt this model to the aqueous Fe III (Cit) / CA system: 1. We always set water activity in the particle phase in equilibrium with the gas phase, since the RH of the gas phase did not change during each experiment.
2. The bulk accommodation coefficients of all species were assumed to be 1.
3. Neglecting the influence of Fe III (Cit), the water activity was taken from that of CA (a CA w ), which was determined from the mass fraction of CA, W CA , by Lienhard et al. (2012): and r = 6.59108 − 0.05294T + 1.06028 × 10 −4 T 2 . (21) The water activity of citrate (a Cit w ) was calculated using the same equation where the mass fraction of citrate W Cit was treated in two fractions: citrate without Fe and citrate containing Fe. The water activity of an aqueous 1 M Fe III (Cit) solution was determined using a water activity meter (AquaLab water, Model 3B, Decagon Devices, USA) at room temperature. We found this water activity corresponds to that of a 0.81 M aqueous CA solution. Therefore, for calculating water activity, all iron-containing citrate complexes (either Fe II or Fe III ) were treated with a reduced concentration: namely with a factor of 0.81 of the corresponding citrate molarity. Hence, the overall amount of citrate was calculated as accordingly and W Cit was determined using where M Cit and M H 2 O are the molar mass of CA and water, respectively. For other species j , the contribution to the water activity is proportional to their molar volume (MV j ), so that in total, 4. Fe III (Cit) photolysis, decarboxylation, and oxidation of the alcohol group in presence of O 2 yields the compounds O=C(CH 2 COO) 2− 2 or O=C(CH 2 COOH) 2 , which are C 5 species. We assumed that half of the C 5 species undergo photochemical reactions to produce CO 2 and compounds with 2-4 carbon atoms, C 4 , C 3 and C 2 (see Reactions R10-R14 shown in Table 2), all of which are capable of being released to the gas phase depending on their solubility. 5. We estimated the quantum yield in Reactions (R1) and (R2) in Table 2 In total, the PRAD model includes 13 equilibria and 17 chemical reactions among 32 species, as well as their condensed phase diffusivities and Henry's law coefficients. Some of these parameters are known from previous studies (see Tables 1 and 2 for references), while others are not known and difficult to estimate. For instance, even though absorption spectra of Fe III (Cit) have been measured in aqueous solution (Pozdnyakov et al., 2012), the corresponding quantum yield has not, which leaves the photolysis rate of Fe III (Cit), j , unknown. Also, there are no data reported on the diffusivity of O 2 in aqueous CA solutions, and the chemical reaction rate of the oxidation of the Fe II citrate complex by O 2 is quite uncertain (Gonzalez et al., 2017). In order to find the optimal parameter set, we compared the experimental data of the three setups taken under well-controlled conditions with model predictions and tuned the unknown parameters manually.
We restricted our tuning of the parameters to reach satisfactory agreement with all experimental data simultaneously. The equilibrium constants and rate coefficients that were tuned are indicated in Table 2 (the sensitivity of the PRAD model results to a few of its parameters is shown in Appendix A5). The parameters were adjusted in a wide and acceptable range until a good representation of our data could be obtained. For example, the fraction of iron(III) in a photoactive complex (equilibrium E5 in Table 2) must have been high enough to reproduce STXM/NEXAFS observations that iron could be reduced to low levels as seen in Fig. 7 described below. In comparison, equilibrium constant of E7 must have been much lower than equilibrium constant of E5 so that the amount of iron(III) in a nonphotoactive complex was small compared to being in complex with citrate. As another example, oxidation of Fe 2+ (Reactions R5-R8 in Table 2) is fairly well referenced, and therefore we adjusted the rate of Reaction (R9) until the model reoxidation rates matched those observed. Tuning of individual bulk diffusion coefficients for all species was not attempted. Instead, we Table 1. Liquid phase diffusivity factors (normalized to water) and Henry's law coefficients (Sander, 2015) of major species in the Fe III (Cit) photochemistry system.

Number Name
Formula unkn. prod C 5 C 5 2.69 × 10 −6 1 × 10 7 8000 32 unkn. prod C 5 C 5 _stable 2.69 × 10 −6 5 × 10 8 8000 a l j f is a factor of the diffusion coefficient of each species j normalized to that of water. b Henry's law is described as a function of temperature T : simplified the representation of diffusion coefficients using a parameterization as a function of molar mass, as described in Appendix A1. The two constants in Eq. (A8) and the two constants in Eq. (A3) were tuned, resulting in the absolute diffusion coefficients shown in Fig. A1. Henry's law coefficients for gasses were tuned, however purposefully set at values higher than expected for pure water or highly dilute aqueous solution. This was inspired by previous studies regularly reporting solubility of, e.g., O 2 and CO 2 higher in a variety of organic liquids than water (Fogg, 1992;Battino et al., 1983). It is important to note that the result of this tuning does not mean that we found the global minimum in the parameter space (see, e.g., (Berkemeier et al., 2017)). A thorough search for a global minimum for our model with 16 tuning parameters for chemistry, four tuning parameters (and our parameterization) for diffusion, and nine tuning parameters for solubility is computationally very expensive and beyond the scope of this paper. However, for our purpose here, namely modeling typical timescales of photochemical degradation of organic aerosol under atmospheric conditions (see Sect. 3.5), the PRAD model framework should allow sufficiently accurate predictions. In other words, we expect similar mass degradation in atmospheric particles due to the fact that many other relevant iron carboxylate compounds undergo LMCT similarly as to our model system (Weller et al., 2013(Weller et al., , 2014. Additionally, if a particular system requires parameter values that significantly differ from ours, the PRAD model framework itself should still be valid. Note that careful evaluation is needed when picking a single parameter of the PRAD model for use in another context. Comparisons of the refined model with our experimental data are shown in the next section.

The effect of RH on photocatalytic degradation efficiency
We performed experiments with single, levitated particles under continuous UV irradiation (375 nm) in pure O 2 at different RH to assess the effects of RH on the photocatalytic cycle shown in Fig. 1 and tested the model performance under these conditions. Qualitatively, a continuous decrease of particle mass and size is expected to occur due to evaporation of volatile products, as shown in Fig. 2. Figure 6 shows the fraction of particle mass remaining with the irradiation time at three different RH calculated from resonance wavelength shifts (Eqs. 1 and 2). Clearly, particle mass was lost to the gas phase with time due to the evaporation of photochemical products and, similarly to Fig. 2, all data show a very significant acceleration of mass loss with time. After tuning some of the parameters of the model as further discussed below, the PRAD model simulations reproduces our data with a very similar trend and magnitude over all, which gives us confidence that the PRAD model captures the essential chemistry and transport during irradiation. (How particle mass evolves subsequently until 80 % mass loss in both experiments and models is shown in Fig. A4.) However, the model is not able to capture the full degree of acceleration of the degradation rate, as it does not attempt to include the complete multigenerational oxidation chemistry at the level of individual components after initial radical production. The degradation processes were faster at higher RH. At lower RH, the particle was expected to be more viscous, diffusion coefficients were expected to be lower, products were generated at a lower rate, and volatile products moved more slowly to the surface to evaporate. More importantly, O 2 taken up by the particle from the gas phase diffused more slowly into the bulk of the particle at lower RH. Thus, less HO .
2 and H 2 O 2 formed and less Fe II could be reoxidized from the surface to the center of the particle. The observed gradient in the Fe III fraction, β, and the modeled gradients in O 2 and ROS in the particle have been shown with radial profiles in Alpert et al. (2021). This resulted in fewer photochemically active Fe III complexes available for photocatalytic degradation. The characteristic degradation time shortened by a factor of 5.5 when RH increased from 46 % and 61 %, which demonstrates that photochemical cycling is highly sensitive to the microphysical conditions. The diffusivity of O 2 must have significantly impacted reoxidation rates. In addition, the diffusion coefficients of both Fe II and Fe III species increases with RH. Therefore, the molecular transport between both iron and oxygen reactants increases causing a highly nonlinear trend in increasing mass loss with increasing RH.

Determination of iron(III) reduction rate and iron(II) reoxidation rate by STXM/NEXAFS
In STXM/NEXAFS experiments, the freshly prepared Fe III (Cit) mixed with CA at x = 1.0 particles were irradiated to determine the Fe III (Cit) photolysis rate, as shown in Fig. 7. Each experimental data point is the average Fe III fraction from 16-36 individual particles. Fitting an exponential function, β = β 0 e j obs t , yields β 0 = 0.93 ± 0.09 and a firstorder decay rate of j obs = 0.08 ± 0.01 s −1 . The LED power at the sample was measured to be 5.9 ± 0.6 mW in total and had a Gaussian spectral profile between 361 and 374 nm at full width half maximum. When mounting the UV fiber optics and collimator lens for multiple samples, the illuminated area had a circle equivalent diameter of 5±1.5 mm. Using the absorption cross section calculated from the molar attenuation coefficient (Pozdnyakov et al., 2008), = 1.0 and propagating all uncertainties yields a photochemical reaction rate of j calc = 0.20 ± 0.12 s −1 , which is in agreement with j obs . This implies that assuming a quantum yield of 1 at these UV wavelengths is reasonable. This calculated value for the photochemical reaction rate was used in the PRAD model for  analyzing the STXM/NEXAFS experiments shown in Fig. 7. We note that our estimate of the photochemical reduction rate and agreement with data is based on the reproducibility of setting up our optical system in the X-ray vacuum chamber. In Fig. 8 (discussed below), this setup procedure was performed for RH = 40 % and RH = 50-60 % independently, and still, iron reduction was in agreement with model predictions in the first minutes after UV light was switched off. The width of the red shading in Fig. 7 is large and the UV-fiber setup was likely the largest source of error for j . A better estimate would require repeat measurements as a function of RH, e.g., to elucidate any systematic uncertainty on iron reduction reactions due to viscosity changes. However, this was not possible as usage of the X-ray beam for STXM/NEXAFS experiments was limited to a few days to complete all experiments. In a different set of experiments, we irradiated particles using a much lower power setting with j = 2.2×10 −3 s −1 in a mixed He and O 2 atmosphere and at a fixed RH for 15 min to reduce Fe III to Fe II . The UV light was then switched off to allow reoxidation in the dark while measuring β over time. Figure 8 shows β as a function of time at RH = 40 %, 50 %, and 60 %. Clearly, the Fe III fraction increased more significantly more slowly with time at drier conditions. While particles were observed to reoxidize to β = 0.7 within 2 h at 60 % RH and expected to be completely reoxidized within about 6 h according to the PRAD model simulations, no significant reoxidation occurred on this timescale for the particles exposed to only 40 % RH. Modeling the reoxidation with the PRAD model yields very satisfactory agreement indicating that the diffusivity parameterizations of the model are capturing the RH dependence of the molecular transport in the viscous matrix. Figure 7. Loss of Fe III in Fe III (Cit) / CA (molar ratio of 1) particles as a function of light exposure time determined using STXM/NEXAFS at 293.5 K. Each data point is the average Fe III fraction over about 16-36 individual particles. An exponential fit, β = β 0 e j obs t , yields for the initial condition β 0 = 0.93 ± 0.09 and a first-order decay rate of j obs = 0.08 ± 0.01 s −1 . The black shading indicates 95 % confidence on the exponential fit. The PRAD model prediction and uncertainty are given as red solid line and shading, respectively, and use a calculated decay rate of j calc = 0.20 ± 0.12 s −1 .

Determination of iron(II) reoxidation rate with single, levitated particles using EDB
As the experiments with single, levitated particles yield only bulk properties and not the oxidation state of iron citrate directly, we designed a dedicated experimental procedure to indirectly determine the reoxidation rate of Fe II . We used multiple irradiation and reoxidation repetitions as shown schematically in Fig. 9. Initially, we exposed a newly injected Fe III (Cit) / CA aqueous particle to blue laser irradiation (473 nm, 4 W cm −2 ) in pure N 2 for 500 s to ensure all Fe III was reduced through photolysis (Reactions R1 and R2 listed in Table 2). As previously described, these reactions led to CO 2 production with subsequent loss of CO 2 to the gas phase, which was observed as a shift in the Mie-resonance wavelength. This shift is shown in Fig. 2b and on an enlarged scale in Fig. 10. After irradiation, we switched the gas flow from N 2 to O 2 in the dark, and Fe II was oxidized back to Fe III over time in this period, either by ROS (Reactions R5-R8) or directly by O 2 (Reaction R9). After a defined time spent in O 2 , the gas flow was switched back to N 2 followed by irradiation to repeat the photolysis step done initially. The ratio of the Mie-resonance wavelength shift of the two photolysis steps was set to be proportional to the ratio of reoxidized Fe III /Fe tot . These two steps (i.e., photolysis in N 2 and reoxidation in O 2 ) were repeated several times, but between each irradiation the particle was exposed to O 2 for different time periods. Following this procedure we intended to map out the characteristic time for reoxidation at various RH. An example of the corresponding raw data (293.5 K and 48 % RH) and retrieved oxidation state is shown Fig. 10. Black circles indicate the first irradiation step with a fresh Fe III (Cit) / CA particle, and red circles indicate the irradiation step that followed the particle after exposure to O 2 for 45 min in the dark. Clearly, the resonance wavelength decreased more during the first irradiation than the second. Therefore, we can conclude unambiguously that the Fe III (Cit) / CA particle initially had more Fe III than what could be reoxidized in O 2 for 45 min. Quantitative scaling, however, requires knowledge of the initial Fe III fraction. Our experiments showed that long exposure (tens of hours) to O 2 Figure 9. Schematic of the experimental procedure to investigate the reoxidation rate of Fe II . Orange columns represent laser irradiation (473 nm, 4 W cm −2 ). Each irradiation takes place in pure N 2 for a period of 500 s. Grey columns mark the recovery process in pure O 2 in the dark; here the time interval is varied. For details, see text.
Figure 10. Temporal evolution of two Mie-resonance wavelength shifts and particle mass changes normalized to the wavelength and mass prior to irradiation; the experiment was done at 48 % RH and 293.5 K. The orange area marks the time with irradiation (473 nm, 4 W cm −2 ). Black: first irradiation with a fresh Fe III (Cit) / CA particle in N 2 , red: irradiation in N 2 right after the particle recovered in O 2 for 45 min in the dark. Two horizontal dashed lines are indicating the final mass remaining after CO 2 loss to the gas phase. For an explanation of the scaling of the third axis, indicating the equivalent Fe III /Fe tot ratio, see text.
yielded larger Mie-resonance shifts than those of the initial photolysis of the freshly prepared particle. This indicated that the initial Fe III fraction was less than 1.0. Hence, we accordingly normalized the Fe III fraction to the data at long (> 15 h) exposure times. For the experiment included in Fig. 10, for example, the initial Fe III fraction of the particle was 0.76, indicating that the particle was partially reduced during sample preparation. After the particle was totally photoreduced, exposure to O 2 for 45 min did not reoxidize all reduced Fe II to Fe III ; only 0.36 Fe III was recovered (as shown by two horizontal dashed lines in Fig. 10).
Another set of experiments was done by starting with freshly injected Fe II (HCit) / CA particles instead of Fe III (Cit) / CA. The only difference in experimental procedure, compared to what is described above, is that there is no first irradiation step. Instead, the particle is initially exposed to O 2 in the dark for a certain time interval, during which Fe II is oxidized only directly by O 2 (Reaction R9). Figure 11. Fe III fraction after different time intervals of recovery in O 2 from EDB observations. Fe III (Cit) / CA (molar ratio of 0.05) particles at 293.5 K with 24 % RH (red diamonds), 48 % RH (green triangles), and 65 % RH (blue squares). Black dots: experiments with Fe II (HCit) / CA particles at 48 % RH. Lines: red (24 % RH), black (48 % RH), and blue (65 % RH) are the Fe III fractions predicted using the PRAD model; the shaded areas indicate model output assuming ±2 % RH uncertainty in the EDB experimental conditions.
Afterwards, the same irradiation and recovery procedures as reported above were used. The equivalent Fe III fractions of all experiments are shown in Fig. 11. Any uncertainty in the normalization of an individual experiment will cause a corresponding uncertainty in the normalized mass loss, which made estimating the uncertainty for individual data points impossible. Nevertheless, from the complete data set, it was evident that with longer time intervals in O 2 , more Fe III was recovered. At 48 % and 65 % RH, about 10 and 3 h exposure to O 2 , respectively, was sufficient for all Fe II to be reoxidized, while at 24 % RH, the recovery even after 25 h was not yet complete. The general trend is consistent with our observation in the STXM/NEXAFS experiments (Sect. 3.2) and can be attributed to molecular diffusion limitations at lower RH: it takes more time for O 2 to diffuse into the particle and for Fe II to diffuse out to the surface of the particle to react with O 2 . It was also evident from these data that Fe II reoxidized by O 2 is as important as Fe II reoxidized by radicals and peroxides as there is no significant difference between the experiments starting from Fe II (HCit) / CA compared to those with Fe III (Cit) / CA. In addition, it indicates that both iron(II) and iron(III) can act as a photocatalyst as long as iron(II) can be oxidized to iron(III), which was also confirmed by Grgić et al. (1999).
By tuning the direct oxidation rate of Fe II (HCit) by O 2 (Reaction R9) and the diffusivity of O 2 , we are able to model the recovery rate at different RH using the PRAD model as shown in Fig. 11. There is satisfactory agreement for the larger RH, but significant underestimation of reoxidation for the experiments at 24 % RH. Our model requires the reaction rate coefficient of Reaction (R9) to be 0.05 M −1 s −1 , which is a factor of 60 smaller than the value that Gonzalez et al. (2017) estimated from their model (3 ± 0.7 M −1 s −1 ). The liquid phase diffusivity of O 2 in our model is 7.1 × 10 −19 m 2 s −1 at 24 % RH, 1.9 × 10 −15 m 2 s −1 at 48 % RH, and 2.1 × 10 −14 m 2 s −1 at 65 % RH. These diffusivities of O 2 are 2-4 orders of magnitude smaller than those of CO 2 determined in Dou et al. (2019). However, we need to stress that in the model, some of the iron-related complex equilibrium constants, their diffusion coefficients, and the Henry's law coefficient of O 2 at different water activities are highly uncertain as well, yielding to a significant uncertainty in the determination of O 2 diffusivity. For example, if the solubility of O 2 was less than what our parameters predict now, a larger O 2 diffusivity would be consistent with our data. In addition, with the total gas flow used in our experiments, it takes about 8 min for replacing the entire EDB gas volume from N 2 to O 2 , and from the response of the particle to flow condition change, we estimate an interval of about 30 min to reach full equilibration to the new gas phase conditions. Therefore, the lifetime of organic radicals needs to be reconsidered. In the PRAD model, we do not take the radical-radical self-reactions (Reaction R16) into account, which may turn out to be a significant sink for the radicals. However, it should be pointed out that the parameter set we have now is a good compromise with additional constraints from STXM/NEXAFS and CWFT experiments.
Another approximate approach to analyze the data of Fig. 11 is to use the analytical solutions for a reacto-diffusive kinetic regime. Here, O 2 taken up from the gas phase by a particle remains confined to a very thin layer below its surface compared to its size provided it reacts reasonably fast with the organic components. Under these conditions there are always pairs of reaction rate and diffusion constants representing the experiments equally well Steimer et al., 2014). In the reacto-diffusive framework, with a constrained reaction rate and Henry's law coefficient of O 2 , the diffusion coefficient of O 2 can be estimated to be 3.6 × 10 −16 and 4.4 × 10 −15 m 2 s −1 at 48 % and 65 % RH, respectively (details are given in Appendix A4). These are both 1 order of magnitude less than the values from PRAD model prediction, but still consistent with each other when considering all uncertainties.

HO . 2 production measured by CWFT experiments
The CWFT experiment allows us to investigate another aspect of the photochemistry of the Fe III (Cit) / CA system. According to Reactions (R1)-(R3) shown in Table 2, the HO .
2 radical is produced upon irradiation and will partition to the gas phase. Figure 12a shows the RH dependence of HO . 2 production, P HO .
2 , from thin films in the CWFT experiment. We observed that P HO .
2 increased with RH when the RH was in-  , respectively. (b) The maximum HO .
2 first-order loss rate determined from the result of the maximum HO .
2 concentration times HO . 2 self-reaction rate, k 5 , which was adjusted to match observations in panel (a) (black dots), or times k 5 parameterized as a function of RH for our PRAD model (Eq. A11) with confidence interval of 95 % (red line and shading), or times k 5 from Corral Arroyo et al. (2018) (dotted blue line), or k 5 = 8.3 × 10 5 M −1 s −1  (dashed blue line). creased from 13 % to 29 % by a factor of about 2. This may be expected since an increase from 13 % to 29 % RH leads to increasing molecular diffusion coefficients and faster chemical cycling (Lienhard et al., 2014;Song et al., 2016). However, at RH between 30 % and 60 %, P HO .
2 decreased with RH, with a production rate at 60 % similar to the rate under dry conditions. This is probably due to the decreasing concentrations of donors (e.g., Fe III (Cit) and Cit . ). But the decrease in concentrations with RH is too small compared to the increase in diffusion coefficients with RH. This means that there must be a strong sink of HO .
2 in the condensed phase when RH increases, which was confirmed by an increasing HO .
2 first-order loss rate as shown by black dots in Fig. 12b. When predicting P HO .
2 with the PRAD model using a constant HO .
2 with RH, as in Corral Arroyo et al. (2018), the model deviates significantly from the data and does not exhibit the observed trend of decreasing P HO .
2 for RH > 30 %. Both assumptions lead to a continuous increase of P HO .
2 with RH as indicated by dashed and dotted blue lines in Fig. 12a, with a decrease in the HO .
2 first-order loss rate as shown in Fig. 12b, which is opposite from our adjustment. Therefore, we argue that the effects of the increasing diffusivity and the stronger sink of HO 2 with RH compensate each other, making P HO .
2 almost independent of RH.

Photochemical degradation under atmospheric conditions
The PRAD model was also developed to be used in more general particle systems. After establishing a parameter set for the PRAD model framework that satisfactorily explains the experimental data obtained with three complementary experimental techniques over a wide parameter range, we used the model to predict photochemical degradation of organic aerosol particles containing carboxylate complexes. For example, an aqueous organic acid particle with a dust inclusion containing iron may exhibit a low enough pH to dissolve part of the iron of the dust inclusion and form photoreactive iron carboxylate complexes (George et al., 2012(George et al., , 2015. If we take the PRAD model as being representative for such a class of particles, we may estimate the degradation of the organic mass by volatilization of products to the gas phase. Figure 13 show examples of such PRAD model predictions: we assume that the organic mass of the particle is composed of aqueous Fe III (Cit) / CA with a molar ratio of the latter being 0.01.
Here, panel (a) illustrates size dependence and panel (b) illustrates RH dependence for particles exposed to air at 293.5 K and an actinic flux at 30 • zenith angle (Madronich et al., 1995). After 12 h exposure to sunlight, the organic mass of the particle (here, CA) was photochemically processed to CO 2 and smaller compounds of high volatility yielding a combined mass loss to the gas phase of up to 50 % depending on size and RH.
To better understand the importance of mass loss in this generalized system, we make a simple and modest comparison of mass loss in Fig. 13 on the order of 20 % over 5 h with ambient mass accumulation measured in the field. Our results are equal to a mass loss rate of about 0.4 µg m −3 (air) h −1 assuming an aerosol population with an organic mass of ∼ 10 µg m −3 (air) undergoing iron carboxylate photochemistry. This rate is much higher than the observed organic mass accumulation in ambient air masses due to photochemical aging during atmospheric transport at about 0.06 µg m −3 h −1 or 6 µg m −3 over 4 d (Zaveri et al., 2012;Moffet et al., 2012). This implies that the mass loss rates are fast enough to affect the balance between aerosol mass accumulation and loss. Hence, photochemical degradation may be very significant in Figure 13. PRAD simulated organic mass loss of Fe III (Cit) / CA (molar ratio of 0.01) for (a) particles of different radii at 40 % RH, and (b) particles with a radius of 100 nm at different RH. All simulations were performed at 293.5 K and for 12 h solar irradiation at 30 • zenith angle (Madronich et al., 1995). iron-containing organic aerosol, perhaps even more important than degradation through reactions with gas phase oxidants.
Note that our model does not fully capture the acceleration of degradation as discussed in Sect. 3.1. Therefore, we argue that the degradation presented in Fig. 13 represents a lower limit of the expected degradation under atmospheric conditions. Clearly, panel (a) indicates that smaller particles degrade significantly faster than larger ones. However, the size dependence is more complex than the expected simple square law if reactions were purely limited by condensed phase diffusion. For particles with a radius larger than 50 nm, the time to repartition 10 % of the mass to the gas phase depends almost linearly on size for these specific conditions. Panel (b) shows that photochemical processes are faster at higher RH corresponding to less viscous particles. The relative mass loss after 12 h irradiation between 20 % and 60 % RH is almost an order of magnitude larger for the highest RH compared to dry conditions because of transport limitations at low RH. However, it should be noted that the relative mass loss at higher humidities (between 40 % and 60 % RH) is quite similar. Here, the photochemical degradation is barely limited by condensed phase diffusivity, but by iron availability in the particle. These simulations show the potential of photochemical degradation under atmospheric conditions; a systematic study exploring the whole range of atmospheric conditions is beyond the scope to this work.

Conclusions
We used three complementary experimental techniques to characterize the impact of reduced mobility of aerosol constituents on photochemical degradation in highly viscous particles. As an atmospherically relevant model system, we chose aqueous Fe III (Cit) / CA particles. These three experimental techniques investigated specific aspects of this photochemical reaction system. In EDB experiments, we measured the mass loss relating to the continual production and loss of CO 2 and other volatile products. We observed very significant condensed phase degradation and strong acceleration of the degradation rate with time. Further studies are needed to quantify all atmospheric implications, but our study suggests that photochemistry in iron-containing organic aerosol will lead to a significant repartitioning of condensed phase mass to the gas phase. We used STXM/NEXAFS to directly measure iron oxidation state in situ with an environmental microreactor. These experiments yielded valuable information about where iron photochemical reduction and reoxidation reactions took place, namely only very close to the surface, and it allowed us to characterize to what degree iron compounds diffused inside single particles. We showed that O 2 uptake and diffusion into a particle is a limiting factor considering the reactions required to produce species with an oxidative potential. In addition, we found that the direct O 2 reaction with iron(II) organic complexes does occur and generates radicals inside the particle. Flow-tube experiments performed on thin Fe III (Cit) films showed continuous production of HO .
2 , revealing a radical source inside the particles driven by photochemistry.
All data were used to constrain equilibrium and kinetic parameters as well as reaction rate coefficients in a new PRAD model with sufficient complexity to allow comparison with data of all experiments simultaneously. In particular, we were able to constrain the photolysis rate of Fe III due to the use of various light sources with various spectral intensities, while capturing the photochemical reduction. In addition, we determined the HO .
2 production rate and its first-order loss rate, and the diffusivity of O 2 in aqueous Fe III (Cit) / CA system as a function of RH and Fe III (Cit) / CA molar ratio with a choice of O 2 related reaction rate coefficient and O 2 Henry's law coefficient.
Although a systematic study exploring the whole range of atmospheric conditions was beyond the scope to this work, there are some aspects of the PRAD model and certain parameters that we argue are reliable and pertinent to atmospheric aerosol photochemistry. First, coefficients in the PRAD model framework can be changed to predict mass loss rates of a different iron carboxylate complex system. We are fairly confident that diffusion coefficients of CO 2 and H 2 O can be used for atmospheric aerosol particles as these were obtained in a more targeted study . Mass loss rates in general are fairly reliable to be used in atmospheric particles as these are linked to photochemical reaction rates that have been characterized (Weller et al., 2013(Weller et al., , 2014. Finally, reoxidation rates and production of radicals are also reliable, as the system is largely reactodiffusion limited (see Appendix A4) and these rates occur on the same scales as observed mass loss rates. In our companion paper (Alpert et al., 2021), we show a detailed analysis of radical concentrations in ambient aerosol particles for a range of atmospheric conditions and iron content. However, our model still needs major improvements, such as including peroxy radical chemistry and better constraints on individual parameters such as diffusion coefficients and reaction rate constants. The overall rate may be well constrained by our experimental studies, however more targeted observations may be necessary for an accurate representation of O 2 chemistry, solubility, and molecular transport independently of each other within aerosol particles.
The chemical evolution of the organic species resulting from the continual photochemical oxidation was not the scope of the present study. A separate study focusing on individual OVOCs and condensed phase products is currently underway that will allow us to better constrain the chemical regimes and the evolution of the oxidation state of the organic fraction with time. This will also allow us to assess more details of organic peroxy radical chemistry that help explain the observed "missing" HO .
2 sink in this system. Furthermore, testing the PRAD model with different organic carboxylate ligands is desirable to broaden its applicability.
Using the PRAD model for predicting photochemical degradation for iron-containing organic aerosol under atmospheric conditions let us conclude that this pathway of repartitioning condensed phase mass to the gas phase is important and its regional and global impact should be investigated in further modeling studies. The PRAD model may serve as a basic framework for the chemistry and transport of compounds in single particles for such studies. Figure A2. Adjusted values of the HO .
2 self-reaction rate used in the PRAD model to exactly reproduce measured P HO .

A3 Parameterization of k 5
Based on the measurement of HO .
2 production as a function of RH using CWFT experiments, the HO .

A4 β estimated by a reacto-diffusive framework
It is well known that multiphase reactions can follow a reacto-diffusive kinetic regime Steimer et al., 2014). For reacto-diffusive limitations as observed in Fig. 11, the Fe III fraction as a function of time was calculated as with the assumption that β 0 = 0, i.e., there was no Fe III at t = 0: where H O 2 and D O 2 l are the Henry's law coefficient and diffusion coefficient, respectively, for O 2 in the Fe III (Cit) / CA matrix, p O 2 is the pressure of O 2 , k O 2 is the total reaction rate of O 2 , and d p is diameter of the single particle. β = 1 e when In a typical EDB experiment, d p = 20 µm, H O 2 = 3.5 × 10 −2 M atm −1 , p O 2 = 8 × 10 4 Pa = 0.789 atm, and k O 2 = 0.05 M −1 s −1 . [Fe tot ] is 0.3192, 0.2763, and 0.2345 M at 24 %, 48 %, and 65 % RH, respectively. From EDB data points in Fig. 11, we estimated that 1 e of Fe III can be fully recovered after around 1.7 h at 48 % RH and 0.45 h at 65 % RH. Thus, the diffusion coefficient of O 2 can be estimated to be 3.6 × 10 −16 and 4.4 × 10 −15 m 2 s −1 at 48 % and 65 % RH, respectively, which are both 1 order of magnitude less than the values from PRAD model prediction, but still consistent with each other when considering all uncertainties. It should be noted that in Eq. (A14), H O 2 D O 2 l k O 2 is actually the constraint, thus any uncertainty in H O 2 or k O 2 can change D O 2 l .

A5 Sensitivity of the PRAD model to various model parameters
As discussed in Sect. 2.5, we performed manual tuning of the PRAD model parameters to reach satisfactory agreement with all experimental data simultaneously. To show the sensitivity of the PRAD model results to a few of its parameters, Fig. A3 shows again the data of the photocatalytic degradation experiment at 46 % RH described in Sect. 3.1. In addition, we show the output of the model as well as model outputs obtained by varying one of the parameters by the indicated factors while keeping all other parameters constant. Clearly, the sensitivities of the model output to varying these parameters are very different: while the model output is quite sensitive to varying the rate constant for the oxidation of Fe II (HCit) (Reaction R9 in Table 2) as well as to oxygen diffusivity, the sensitivity to the equilibrium constant E8, Fe 2+ +HCit 2− Fe II (HCit), is significantly smaller and the model is basically insensitive to the equilibrium constant E5, Fe 3+ + Cit 3− Fe III (Cit). In other words, this experiment alone allows us to constrain Reaction (R9) or oxygen diffusivity as long as the other parameter is known, but does not allow us to constrain the constants for the two equilibria.
A6 The full EDB data set showing the RH effect on photocatalytic degradation efficiency Figure 6 in the main text shows the photochemical degradation up to a mass remaining ratio of 0.6 and compares its temporal evolution with the PRAD model output. Here, we  Table 2) and panel (d) shows the sensitivity to oxygen diffusivity (Table 1). Figure A4. Fe III (Cit) / CA (molar ratio 0.05) particle mass change with irradiation time (log scale) at different RH: 33 % (red), 46 % (green), and 61 % (blue). The irradiation wavelength was 375 nm, its intensity was 0.25 W cm −2 , and the experimental temperature was 293.5 K. Thick lines are EDB experimental data. Thin lines are PRAD model outputs (with ±2 % RH uncertainty shown as the shaded area).
show in Fig. A4 for completeness the full experimental data set. Clearly seen is the shift to longer time scales with decreasing RH and that the PRAD model is no longer capable of accurately capturing the mass loss once about 20 % of the initial mass is lost. However, we may conclude that even for mass remaining ratios lower than 0.5, the experimental data show a very significant shift of the degradation time scales with RH, but the total photochemical mass loss remains of the same order of magnitude, independent of RH.
Code and data availability. The data that support the findings of this study are available from the corresponding authors upon request. The PRAD model code is publicly available and accessible here: https://doi.org/10.3929/ethz-b-000451609 (Dou et al., 2020), and from the corresponding author upon request.
Author contributions. JD wrote the manuscript. MA and UKK conceptualized and planned the study. JD conducted the EDB experiments, water activity measurements, and data analysis, and the interpretation was supervised by UKK. PAA planned and conducted the STXM/NEXAFS experiments supervised by BW and MA. JD, PCA, JX, TH, CNB, and KDH also conducted STXM/NEXAFS experiments. BW, JR, and PAA conducted STXM/NEXAFS data analysis and interpretation. PCA conducted CWFT experiments, and the data analysis and interpretation was supervised by MA. FS conducted viscosity experiments and the data analysis and interpretation was supervised by PAA. BL wrote and developed the PRAD model with the assistance of JD and PAA. HH advised on the photochemical reaction mechanisms. All coauthors discussed the results and commented on the manuscript.