Atmospheric oxidation of new “green” solvents – Part 2: methyl pivalate and pinacolone

. Lab-based experimental and computational methods were used to study the atmospheric degradation of two promising “green” solvents: pinacolone, (CH 3 ) 3 CC(O)CH 3 , and methyl pivalate, (CH 3 ) 3 CC(O)OCH 3 . Pulsed laser photolysis coupled to pulsed laser-induced ﬂuorescence was used to determine absolute rate coef-ﬁcients (in 10 − 12 cm 3 molec


Introduction
Solvents are of fundamental importance to the chemical industry; however many common solvents are sourced from the petrochemical industry and are classed as hazardous pollutants.Their usage should be phased out as we move towards a net-zero carbon economy (Clarke et al., 2018;Byrne et al., 2016;Winterton, 2021).Solvents were the dominant source of volatile organic compound (VOC) emissions in England in 2019 (Jones et al., 2021) and so can have a potentially important impact on air quality.A solvent of particular concern is the monoaromatic species toluene (C 6 H 5 CH 3 ) (Calvert et al., 2002), known to pose significant human health and environmental hazards (Byrne et al., 2018) and classified as a significant air pollutant (Montero-Montoya et al., 2018;Zhang et al., 2019).This work follows on from Mapelli et al. (2022), where the atmospheric oxidation of 2,2,5,5-tetramethyloxolane (TMO) (Byrne et al., 2017) was studied in detail (Part 1).Byrne et al. (2018) similarly identified pinacolone (3,3-dimethyl-2-butanone, (CH 3 ) 3 CC(O)CH 3 ), henceforth PCO, and methyl pivalate (methyl 2,2-dimethylpropanoate, (CH 3 ) 3 CC(O)OCH 3 ), henceforth MPA, as potential replacements for hazardous volatile non-polar solvents such as toluene.Advantages include low toxicity, sustainable (biomass) sourcing and good solvation characteristics.However, whilst the atmospheric chemistry of toluene is well known and its deleterious air quality impacts have been quantified, the atmospheric breakdown of PCO or of MPA was seldom studied.No air quality considerations were included in their selection as "green" solvents (Byrne et al., 2018).
Industrially, PCO is produced for use in fungicides, herbicides and pesticides (Liu et al., 2022).In addition to direct emission from industry, the Master Chemical Mechanism (MCM; http://mcm.york.ac.uk last access: 12 December 2022) (Saunders et al., 2003;Jenkin et al., 1997), a benchmark detailed description of atmospheric chemistry, reports PCO as a product of gas-phase oxidation of both 2,2-dimethylbutane and 3,3-dimethyl-2-butanol.An important removal process from the troposphere for virtually all VOCs including PCO, is via reaction with the hydroxyl radical, OH (Reaction R1; see Fig. 1).
An alternative breakdown process for carbonyl-containing VOCs in the atmosphere is photolysis (Reaction R4), where the C = O function may absorb abundant UV light at λ > 290 nm via the n → π * electronic transition., where the H abstraction may proceed at the alpha CH 2 −H (Reaction R1a) or beta CH 2 −H (Reaction R1b, as proposed in the MCM).Similarly, two different sites for H abstraction are available for the reaction MPA + OH (Reaction R2).In the atmosphere, these radical products will react rapidly with oxygen to form peroxy radicals (RO 2 ).
PCO + hν(λ > 290 nm) → products (R4) PCO is known to have an absorption maximum at 287 nm in cyclohexane solution (Chimichi and Mealli, 1992) and at 278 nm in water (Pocker et al., 1988), though quantitative absorption cross-sections were not reported.In the absence of quantitative cross-sections or quantum yields, the MCM represents PCO photodissociation (Reaction R4) and subsequent rapid reactions with O 2 as exclusively producing (CH 3 ) 3 CO 2 + CH 3 C(O)O 2 and the photolysis rate parameters of PCO (and seven other ketones) via parameters measured for photolysis (Reaction R5) of methyl ethyl ketone (butanone, CH 3 C(O)C 2 H 6 , henceforth MEK).
The UV absorption spectrum of MPA was also measured by Wallington et al. (2001), who reported an absorption maximum at around 210 nm but no absorption at wavelengths (λ > 290 nm) where UV light is abundant in the troposphere.Reactions of OH with oxygenated volatile organic compounds (OVOCs) such as PCO and MPA are notoriously complex, with several examples displaying bimodal non-Arrhenius kinetic behaviour at atmospherically relevant temperatures (Calvert et al., 2011;Wollenhaupt et al., 2000).The lack of photochemical data for degradation of PCO and MPA in the troposphere, and OVOCs in general, is therefore concerning.Rate coefficients k 1 and k 2 were each reported just once (Wallington and Kurylo, 1987a;Wallington et al., 2001;Mapelli et al., 2022).Photolysis (Reaction R4) is potentially the dominant loss process for PCO, yet the stateof-science model MCM (Master Chemical Mechanism) represents Reaction (R4) using absorption cross-sections and quantum yields for a surrogate ketone species, MEK (Saunders et al., 2003).This paucity of data, together with the recent increase in interest around green solvents, was a prime motivation for this work, which aims to improve our knowledge of the atmospheric chemistry of the OVOC species PCO and MPA.Direct, absolute kinetic methods were used to determine temperature dependent rate coefficients for the reaction of PCO and MPA with OH.Absorption cross-sections for PCO were obtained from both laboratory and computational experiments.These fundamental photochemical parameters were then used to calculate lifetimes for the removal of PCO and MPA from the troposphere.Finally, estimated photochemical ozone creation potential (POCP E ) values were calculated according to the method proposed by Jenkin et al. (2017), allowing for air quality impact assessments of these molecules for use as green solvents.

Experimental and computational methods
Two distinct sets of experiments were used to examine solvent breakdown reactions.In Sect.2.1, a brief description of the University of York pulsed laser apparatus is provided.This setup, used to obtain absolute rate coefficients k 1 (295-485 K) and k 2 (295-485 K), was described more fully in a recent publication (Mapelli et al., 2022).Section 2.2 describes UV-Vis measurements of solution-phase absorption crosssections, σ (λ) for PCO and for MEK, a compound that has measured gas-phase absorption cross-sections and that, to date, has been used to represent PCO photolysis (as well as a range of structurally similar ketone species) in models of atmospheric chemistry (Jenkin et al., 2003;Saunders et al., 2003).The computational methods used for in silico σ (λ) calculations are presented in Sect.2.3.

Pulsed laser photolysis-laser-induced fluorescence (PLP-LIF) experiments
The absolute kinetic technique of pulsed laser photolysis (PLP), coupled to direct pulsed laser-induced fluorescence (LIF) detection of OH, was used to determine rate coefficients k 1 (295-485 K) and k 2 (295-485 K).This apparatus was described in detail in a previous publication (Mapelli et al., 2022).Briefly, output from four mass flow controllers was mixed prior to entering a 400 cm 3 Pyrex reactor equipped with thermocouples (±2 K) and capacitance manometers.Output from these calibrated analogue devices was digitised and passed to a PC to regulate and log gas temperature, pressure and flow rate (typical flow rate of 1000 sccm at 60 Torr (8.0 kPa, where 1 Torr = 0.1333 kPa).
where S 0 refers to the LIF signal at t = 0 and is proportional to the initial concentration of OH generated by Reaction (R6) and B represents the pseudo-first-order rate coefficient.
2.2 Absorption cross-sections, σ(λ), via UV-Vis experiments UV-Vis studies were carried out using a double-beam Shimadzu UV-2600 spectrometer over the wavelength range 250-400 nm at T = (298 ± 1) K and with a resolution of 1 nm.The spectra were recorded in Hellma Analytics quartz cuvettes with a path length of 10 or 2 mm.As no suitable apparatus was available at the University of York for gas-phase σ 4 (λ) determinations, all spectra recorded were in solution phase, diluted in cyclohexane (c-C 6 H 12 ).This choice was based on the work of Nakashima et al. (1982), who studied the solvent effect on the absorption spectra of ketones.Spectra in perfluorohexane (C 6 F 14 ) were reported as quasivapour spectra.In this work, use of C 6 F 14 was avoided for environmental reasons (Tsai, 2009); hence the similarly nonpolar and weakly polarisable cyclohexane was chosen to imitate C 6 F 14 .Solutions (0.03-0.3 M) were prepared using volumetric flasks and micropipettes, with the concentration uncertainty estimated at around 10 %.

Computational methods used to determine σ(λ)
Ground-state minimum-energy geometries for PCO and MEK were located with density functional theory (DFT) calculations employing the ωB97XD (Chai and Head-Gordon, 2008) exchange-correlation functional and the def2-TZVP basis set (Weigend and Ahlrichs, 2005), both in gas phase and using an implicit solvent within the IEF-PCM framework (polarizable continuum model using the integral equation formalism variant; using parameters for cyclohexane).Each stationary point located was confirmed to be a minimum by performing a frequency calculation at the very same level of theory.Two different conformers were located for MEK (in gas phase and with the implicit solvent model).
The free-energy difference between the two conformers is 0.071 eV in gas phase and 0.059 eV in solution (corresponding to a Boltzmann population ratio of 0.942 : 0.058 in gas phase and 0.909 : 0.091 in solvent at 298.15 K).Photoabsorption cross-sections were calculated following the workflow described by Prlj et al. (2022).In short, a Wigner function for uncoupled harmonic oscillator was created from the frequencies (discarding the lowest-energy mode, < 100 cm −1 ) and equilibrium geometry for each molecule and each conformer, from which 1000 different geometries were sampled.The nuclear ensemble approach (Crespo-Otero and Barbatti, 2012) was then employed to construct the photoabsorption cross-sections for each molecule (and conformer) from the transition energy and oscillator strength (S 1 ← S 0 ) of each sampled geometry obtained with linear-response time-dependent DFT (LR-TDDFT) within the Tamm-Dancoff approximation (TDA), using the same exchange-correlation functional and basis set as earlier (ωB97XD and def2-TZVP), with and without a LR-PCM (Cammi et al., 2000).This level of electronicstructure theory was benchmarked against equation-ofmotion coupled-cluster singles and doubles (EOM-CCSD) and was found to produce reliable excitation energies and oscillator strengths in earlier studies (Sarkar et al., 2021).LR-TDDFT/TDA/ωB97XD/def2-TZVP gives a vertical excitation (oscillator strength) towards S 1 of 4.61 eV (0.0000) for MEK and 4.49 eV (0.0002) for PCO, while EOM-CCSD/aug-cc-pVDZ places these transitions at 4.62 eV (0.0000) for MEK and 4.50 eV (0.0001) for PCO.The final photoabsorption cross-section of MEK (in gas phase and in solution) was obtained by averaging the photoabsorption cross-section of each conformer, weighted by the corresponding Boltzmann population.The nuclear ensemble approach (NEA) calculations were performed with Newton-X version 2.4 (Barbatti, 2013;Barbatti et al., 2014), with electronic-structure information provided by Gaussian16 Revision B.01 (Frisch et al., 2016).

Results and discussion
Results from absolute kinetic determinations of k 1 (295-485 K) and k 2 (295-485 K) are presented in Sect.3.1.Absorption cross-section data determined in this work are presented in Sect.3.2 (experimental UV-Vis results) and 3.3 (results from computational studies).In Sect.3.4, these results are discussed in full, with an emphasis on the relative importance of photolysis vs. the OH reaction route and on outstanding uncertainties.Results from all such k 1 (T ) determinations are listed in Table 1 and displayed in Arrhenius format in Fig. 3.A weighted average of the three values at close to ambient temperature yields k 1 (297 K) = (1.23 ± 0.09) × 10 −12 cm 3 molec.−1 s −1 , in good agreement with the solitary previous determination of k 1 (296 K) = (1.21± 0.05) × 10 −12 cm 3 molec.−1 s −1 (Wallington and Kurylo, 1987a).Over the range of temperatures investigated in this work, these results were adequately represented by a two-parameter Arrhenius equation, as k 1 (297-485 K) = 2.14 × 10 −12 exp(−200/T ) cm 3 molec.−1 s −1 .Also displayed in Fig. 3 are results from this work on the well-characterised reaction between ethanol (CH 3 CH 2 OH) and OH (Reaction R3), many of which were obtained in back-to-back experiments alongside those on Reaction (R1).The satisfactory agreement between the k 3 (295-447 K) determinations here and the extensive literature dataset (Dillon et al., 2005;Carr et al., 2008;Wallington and Kurylo, 1987b;Hess and Tully, 1988;Jiménez et al., 2003) lends confidence to the k 1 (T ) results from this work.A similar series of experiments was undertaken to determine k 2 (295-485 K), results from which are presented in Table 1 and depicted in Arrhenius format in Fig. 4. Control experiments were conducted at different pressures and in both N 2 and air bath gas to ensure the reactivity was not affected by these factors; again, many experiments were conducted back to back with equivalent k 6 determinations (see above).A weighted mean of the three room temperature determinations led to a value of k 2 (297 K) = (1.30± 0.01) × 10 −12 cm 3 molec.−1 s −1 .This result agrees well with the ambient temperature value reported by Wallington et al. (2001) of k 2 (296 K) = (1.20 ± 0.03) × 10 −12 cm 3 molec.−1 s −1 when considering potential systematic uncertainties, notably in [MPA] (see Sect. 3.4 below).As is evident from Fig. 4, results from this work confirmed the complex non-Arrhenius behaviour reported by Wallington et al. (2001), revealing a local minimum in k 2 (T ) close to ambient temperature.Such non-Arrhenius behaviour may be indicative of a change of rate-determining step within the mechanism, with direct hydrogen abstraction dominating at high temperatures, whilst pathways via hydrogen-bonded pre-reaction complexes play an increasingly important role at lower temperatures.A similar trend was observed for reactions of OH with other OVOCs (Wollenhaupt et al., 2000;Vasvári et al., 2001).As displayed by the dashed black line in Fig. 4, the full set of k 2 (T )  Wallington et al. (1987); the dashed black line denotes an Arrhenius fit to all k 1 (T ) results to yield k 1 (296 or 297-485 K) = 2.14 × 10 −12 exp(−200/T ) cm 3 molec.−1 s −1 ; and the solid black line represents k 1 (T ) calculated using the SAR proposed by Jenkin et al. (2018).Also displayed are k 3 (T ) values from this work (black triangles) and from literature (Dillon et al., 2005;Hess and Tully, 1989;Carr et al., 2008;Wallington and Kurylo, 1987b;Jiménez et al., 2003) together with the corresponding IUPAC evaluation (Atkinson et al., 2006), leading to the three-parameter expression k 3 (T ) = 6.70 × 10 −18 T 2 exp(511/T ) cm 3 molec.−1 s −1 , here represented with a dashed red line.

UV-Vis cross-sections σ(λ) determined in the laboratory
Figure 5 displays the solution-phase spectrum for PCO (in cyclohexane) obtained using the methods detailed in Sect.2.2 above.Previous spectra of PCO over the wavelength range of interest were not available, except for the absorption maximum recorded for PCO in water (Pocker et al., 1988) and in cyclohexane (Chimichi and Mealli, 1992) at 278 and 287 nm respectively.That such dilute, solution-phase spectra can yield useful results is demonstrated by equivalent spectra for MEK.Cross-section values similarly obtained in this work, dilute in cyclohexane, are displayed in Fig. 5 alongside the IUPAC-recommended gas-phase spectrum (Mellouki et al., 2021) reported by Martinez et al. (1992) and confirmed by Brewer et al. (2019) in more recent studies.The solution-phase determinations https://doi.org/10.5194/acp-23-7767-2023Atmos.Chem.Phys., 23, 7767-7779, 2023 consistently overestimate reported gas-phase values (e.g. at maximum value, σ (278 nm) = 6.2 × 10 −20 cm 2 molec.−1 from this work vs. σ (278 nm) = 5.8 × 10 −20 cm 2 molec.−1 in the gas phase (Martinez et al., 1992); see the Supplement).Whilst this was within the 10 % uncertainty estimated for this work, observing the spectra over at least three repeats from different solutions, the overestimation appears to be systematic (see Supplement) and was also identified in silico (see Sect. 3.3 below).The inset in Fig. 5 displays Beer-Lambert plots for data recorded in this work at 290 nm, yielding the molar absorption coefficients ε 4(290 nm) = (20.6 ± 0.2) cm −1 M −1 for PCO and ε 5(290 nm) = (16.4± 0.2) cm −1 M −1 for MEK.These coefficients can then be converted to the absorption cross-sections, leading to σ (290 nm) = (8.0 ± 0.8) × 10 −20 cm 2 molec.−1 for PCO and σ (290 nm) = (6.0 ± 0.6) × 10 −20 cm 2 molec.−1 for MEK.The UV spectra of PCO determined in this work were clearly characterised by larger cross-sections in the actinic region (see Fig. 6) and by a significantly red-shifted local maximum in absorption (286 vs. 278 nm) when compared to MEK.These observations will therefore have significant implications when revising estimates for photolysis rates of PCO, as more light is available at longer wavelengths (below).In accord with previous observations, UV-Vis spectra recorded for MPA showed no absorption over the wavelength range of interest, 290 < λ/nm < 370 (see Supplement).

UV-Vis cross-sections σ(λ) as determined in silico
Computational photochemistry can be used to estimate photoabsorption cross-sections both in gas phase and within an implicit solvent and determine what the theoretical shift in wavelength and intensity caused by cyclohexane is.The nuclear ensemble approach (NEA) (Crespo-Otero and Barbatti, 2012) was used to predict the photoabsorption crosssection for PCO and MEK (see Sect. 2.3 for the computational details), and we propose here a brief discussion of this technique as it is not commonly applied for the calculation    2).An actinic flux (the dashed green line), for conditions of ϑ = 60 • and 350 DU from the NCAR (TUV) Radiation Model, together with φ (275-380 nm) = 0.16 was used for all (Eq.2) j -value calculations.The striped blue area derives from use of σ 4 (λ) values recorded in experiments using cyclohexane solvent having applied the scaling factor of 0.93 to best reproduce gas-phase σ 4 (λ) values (see Sect. 3.5), integrating to j 4 = 2.4 × 10 −6 s −1 .Also displayed are j values for photolysis of MEK (Reaction R5) calculated using σ 5 (λ) from the gas-phase spectrum of Martinez et al. (1992) (red square area, integrated to j 5 = 0.9 × 10 −6 s −1 ).
of UV-Vis cross-sections.A common strategy to approximate a photoabsorption cross-section using quantum chemistry -the single-point approach -consists in locating the minimum-energy geometry in the ground electronic state for the molecule of interest and to calculate for this geometry the vertical transitions and oscillator strengths towards the excited electronic states.A photoabsorption cross-section can be obtained by broadening this single line with a Gaussian or Lorentzian function.The issue with the single-point approach is that it does not account properly for the quantum delocalisation of the ground-state nuclear wave function and non-Condon effects, that is, the fact that the transition dipole moment may vary with the nuclear configuration -an aspect particularly important for symmetry-forbidden transitions.
The NEA offers a rather straightforward strategy to calculate photoabsorption cross-sections that accounts for the broadening coming from quantum delocalisation and the possible enhancement in transition probability caused by non-Condon effects.We note that the NEA does not account for contributions from Franck-Condon factors, that is, overlaps of nuclear wave functions on the ground and excited electronic states, and as such does not reproduce the vibronic structure of a given electronic transition.We refer the interested reader to predominantly used to investigate the effect of a solvent (cyclohexane) on the photoabsorption cross-section of MEK and PCO.
The NEA photoabsorption cross-section calculated for MEK in gas phase shows a close agreement with the experimental cross-section in terms of intensity and centre of the S 1 ← S 0 band, which exhibits a nπ * character.Including an implicit solvent (cyclohexane) leads to a slightly more intense cross-section predicted by the NEA for MEK, again mimicking closely the shift observed experimentally when comparing the gas-phase cross-section with the MEK cross-section measured in cyclohexane.Comparing the MEK cross-section calculated in cyclohexane with that obtained for PCO shows that the latter is more intense and red-shifted, in line with the measured cross-section in cyclohexane.The difference in width between the experimental and theoretical cross-sections is caused by the approximations underlying the NEA (neglect of vibrational progressions).From a theoretical perspective, comparing the changes in photoabsorption cross-sections in gas phase versus cyclohexane for MEK and PCO shows that the gas-phase photoabsorption cross-section is weakly red-shifted (0.0275 eV on average) with respect to the cross-section in cyclohexane, and its intensity is weakened by a factor 0.86.

Estimation of photolysis rate coefficients (j values)
Having measured the absorption cross-section of PCO, the photolysis rate can then be estimated according to Eq. ( 2): where j is the first-order photolysis rate coefficient, φ stands for the quantum yield, σ is the absorption cross-section at wavelength λ and temperature T , and F is the actinic flux at wavelength λ and solar zenith angle θ (Fig. 7).As for many other organic compounds, quantum yields for PCO are unknown, and the MCM uses MEK as a surrogate with reasonably established photochemistry.Accordingly, the photolysis rate was calculated using φ = 0.16, the quantum yield determined by Pinho et al. (2005) for the wavelength range 275 to 380 nm.This value is different than that recommended by the IUPAC of φ = 0.34 (Atkinson et al., 2006), with the value from Pinho et al. (2005) optimised with environmental chamber data.Subsequently, Romero et al. (2005) reported quantum yields lower than the one estimated by Raber and Moortgat (1995), and their work stressed the importance of investigating the temperature dependence and the spectral distribution of quantum yields to avoid overestimation.The actinic flux used for the estimation of j values was obtained using the NCAR Tropospheric Ultraviolet and Visible (TUV) Radiation Model (Madronich and Flocke, 1997), with a solar zenith angle of 60 • , an O 3 column of 350 DU (Dobson units), an altitude of 0 km and the date 30 June 2015.
The integrated photolysis rate coefficient calculated via Eq.( 2) (ϕ = 0.16) for PCO gives a value of j 4 = 2.4 × 10 −6 s −1 in the considered conditions, using σ 4 (λ) values from the UV-Vis spectrum recorded in cyclohexane, scaled by the solution-to-gas scaling f 0.93 , determined as described above.Calculation of j for MEK leads to quite a different result (Fig. 6), with j 5 = 0.8 × 10 −6 s −1 estimated using the cross-sections determined in this work (see Fig. S6 in the Supplement) and j 5 = 0.9 × 10 −6 s −1 , determined using cross-sections from Martinez et al. (1992) (Fig. 6).Photolysis rates were also calculated using the parameterised φ(λ, [M], T ) formulated by Romero et al. (2005).Using the parameterised quantum yield leads to fairly similar results, with λ-dependent values larger than MCM values by a factor ∼ 1.2.Table 2 summarises the results obtained from all the available cross-sections and quantum yields.All these values are in reasonable agreement with the j 5 = 1.3 × 10 −6 s −1 calculated according to the MCM parameterisation at a similar zenith angle (Saunders et al., 2003).

Discussion
The study of Reactions (R1) and (R2) via PLP-LIF to determine the rate coefficients of these reactions led to the values reported above with their (statistical only) standard error.Considering the systemic errors, we quote more realistic values of k 1 (296 K) = (1.2 ± 0.2) × 10 −12 cm 3 molec.−1 s −1 and k 2 (296 K) = (1.3 ± 0.2) × 10 −12 cm 3 molec.−1 s −1 , which take into account the error over the estimation of [VOC].Both k 1 and k 2 are in good agreement with the literature data by Wallington and Kurylo (1987a) and by Wallington et al. (2001).Whilst this agreement was pleasing, all three studies used similar PLP methods and thus have a similar reliance on [VOC] measurements; confirmation of these results by, for example, the relative rate technique would be worthwhile in future.These results are also in agreement with the rate coefficients estimated via the structure-activity relationship (SAR) method elaborated by Jenkin et al. (2018) (Table 2).In addition, SAR product distribution for Reaction (R1) is in accordance with the one observed by Wallington and Kurylo (1987a) and recommended by the MCM (Master Chemical Mechanism; Jenkin et al., 2003;Saunders et al., 2003), which indicates Reaction (R1b) as the dominant route (Fig. 1, Table 3).
However, as is evident in Figs. 3 and 4, the temperature dependence of the SAR predictions is still unable to predict those observed experimentally for those two OVOC reactions.In the case of MPA, where results from this work together with literature data from Wallington et al. (2001) reveal a "U-shaped" non-Arrhenius k(T ) described by a fourparameter expression, the SAR would appear to underestimate the relative importance of the conventional Arrheniuslike contribution to the overall k(T ), especially at T > 300 K.The SAR does appear to account for the complex non- a j values estimated for conditions of ϑ = 60 • and 350 DU from the NCAR (TUV) Radiation Model.b Gas-phase σ (λ) by Martinez et al. (1992).c Solution-phase σ (λ) determined in this work in cyclohexane.d Scaled σ (λ) calculated scaling the solution-phase values by a factor 0.93, estimated from the ratio between the integrated gas and solution MEK spectra.e j values calculated for (275-380 nm) φ = 0.16, a constant average value estimated by Pinho et al. (2005).f j values estimated for ϕ = ϕ(λ, T ), the parameterised quantum yield formulated by Romero et al. (2005), calculated over the wavelength range 280-335 and at T = 296 K.
Arrhenius behaviour at lower temperatures, where direct H abstraction is slow, and the reaction proceeds (by analogy with similar OH + OVOC reactions) via formation of sixmembered ring complexes featuring a hydrogen bond between the carbonyl function and the OH radical.Jenkin et al. (1998) do note that additional rate coefficient data would be highly valuable for further evaluation and the constraining of the SARs for OVOC and multifunctional species.
For OH + PCO (Reaction R1), the experimental results have not provided enough evidence for any complex non-Arrhenius temperature dependence; experiments at lower temperature may yet bring further clarification.The k(T ) values from SAR so predict a non-Arrhenius trend similar to the one calculated for MPA, somewhat different from the experimental results, although in reasonable agreement in ambient temperature.Some caution should therefore be exercised when using the SAR to estimate product yields for Reactions (R1)-(R2), even at T ≈ 298 K, where the SAR appears to predict accurate k values but potentially for the wrong reasons.
The lack of absorption cross-section and quantum yield data for many oxygenated organic molecules is an obstacle to the estimation of photolysis rates.Here, UV-Vis spectra of PCO and MPA were recorded in the actinic light range (λ > 290 nm) to help clarify the role of photolysis as an alternative (to OH reaction) removal path to oxidation.Cyclohex-ane was selected as a solvent for these solution-phase spectra (Nakashima et al., 1982); tests on the well-characterised molecule indicate that such solution-phase spectra can provide reasonable approximation of the gas-phase equivalent spectra, with a minor overestimation of the gas-phase values (see Fig. 5).Our UV-Vis results indicate a ratio of 0.93 between the integrated gas-phase spectrum by Martinez et al. (1992) and the one recorded in solution in this work (See Supplement) in the spectral range 300-330 nm.Computational photochemistry calculations validated the experimental results, confirming an apparently systematic increase in photoabsorption when measuring the cross-section in a cyclohexane solution.More precisely, a similar factor of 0.86 : 1 was calculated via NEA for the gas-phase spectra to the solution-phase spectra.
Assuming then that PCO has a similar behaviour to that observed for MEK, absorption cross-sections recorded in cyclohexane solution may be used in photolysis rate determinations.In the absence of any quantum yield data for PCO, a value of φ = 0.16, the quantum yield used for MEK by the MCM, was used here.However, σ 4 (290-330 nm) values were more intense than σ 5 (290-330 nm), and the peak was significantly red-shifted (see Fig. 5) to wavelengths where light is abundant in the troposphere.Using the σ 4 photolysis rate determined in this work for PCO, we estimated j 4 = 2.4 × 10 −6 s −1 (Fig. S6), and applying a scaling fachttps://doi.org/10.5194/acp-23-7767-2023Atmos.Chem.Phys., 23, 7767-7779, 2023 tor of 0.93 to σ 4 (λ), we estimated a very similar value of j 4 = 2.4 × 10 −6 s −1 (Fig. 6).Accordingly, PCO shows a significantly larger photolysis rate than the one estimated for MEK under similar atmospheric conditions (Fig. 6), with j 5 = 0.9 × 10 −6 s −1 .As a confirmation, the spectra recorded in silico for PCO and MEK describe a similar trend, as illustrated in Fig. 5.
Clearly the use of MEK as proxy for PCO in the MCM and other chemical models is not satisfactory and is likely to lead to an underestimation of the role of photolysis as regards PCO degradation.The lack of ϕ(λ) data means that a more accurate evaluation of the photolysis rate is not possible.

Atmospheric implications and conclusions
The atmospheric chemistry of PCO and MPA has been investigated and compared to the limited number of previous studies.Both MPA and PCO react slowly with OH when compared to toluene, which reacts with OH about 5 times faster than either of these oxygenates (Mellouki et al., 2021).
From determinations of k 1 (296 K) and k 2 (296 K) in this work, we have estimated lifetimes (Eq. 3) for PCO and MPA with respect to reaction with OH in the troposphere: Using a mean tropospheric [OH] = 1.13 × 10 6 molec.cm −3 (Lelieveld et al., 2016), the calculated lifetime with respect to OH is 9 d for PCO and 9 d for MPA.These lifetimes suggest a relatively low atmospheric reactivity, especially when compared to toluene (τ = 2 d), and may allow for some dispersal or the primary emission prior to formation of ozone and other secondary pollutants.As discussed above, the UV-visible spectrum of MPA suggests that photolysis in the troposphere (λ > 290 nm) is not significant, and the main chemical loss process is therefore reaction with OH.By contrast, PCO absorbs light at 290 nm and above (Fig. 5) and hence (by analogy to the known photochemistry of other ketones) will incur photolysis losses.The photolysis rate estimated in this work (j 4 = 2.4 × 10 −6 s −1 ), using the scaled solution-phase cross-section data and the quantum yield recommended by the MCM, produces a lifetime of 5 d with respect to photolysis (τ p ) under the selected atmospheric conditions.Overall, the lifetime of PCO can be estimated using Eq. ( 4), shortening τ to 3 d (Table 4).
Although this result strongly depends on the atmospheric conditions considered, and despite the acknowledged uncertainty on the estimation of j because of lack of gas-phase absorption cross-section and quantum yield data, this shorter lifetime indicates how photolysis represents an important process for PCO and may outcompete reaction with OH as a removal pathway.Photochemical parameters obtained in this work can further be used in chemical models and in calculations of estimated photochemical ozone creation potential (POCP E ) values for environmental assessments.The method was introduced by Jenkin et al. (2017) to take size and structural features of the molecule, reactivity with OH, and the presence of suitable chromophores into account.Here, POCP E values for north-western European conditions were calculated for PCO and MPA (Table 4) using photochemical data determined in this work.The value of POCP E = 11 for MPA may be quoted with some confidence, given the good agreement in rate coefficients determined here and in previous work and the absence of an active chromophore at relevant wavelengths.By contrast, three scenarios were considered for PCO, given the remaining uncertainties in the rate of Reaction (R4).First, following the guidance in Jenkin et al. (2017), whereby one single value is recommended to account for photolysis of all aliphatic ketones, we calculated POCP E = 26 for PCO.However, this photolysis parameter appears to be based upon cross-section and quantum yield values from the MCM, determined for photolysis of MEK (Reaction R5).We cannot rule out that the (unmeasured) quantum yields for PCO photolysis are smaller than 0.16.T dependence of ϕ should also be considered and investigated, and considering only room temperature ϕ may lead to overestimation, as pointed out by Romero et al. (2005); in our second scenario we estimate that a consequently small value of POCP E = 12 may be appropriate (ϕ = 0, so no photolysis; OH loss only).A third scenario would account for the enhanced rate of photolysis of PCO due to the more intense and red-shifted σ 4 (λ) determined in this work and/or a higher quantum yield (0.16 < ϕ < 1).However, a quantitative evaluation of this consequently enhanced POCP E > 26 is beyond the scope of this work and would require detailed modelling studies.
Comparing these POCP E values to traditional hydrocarbon solvents that MPA and PCO are meant to replace (Byrne et al., 2018), these sustainable and non-toxic solvents show a much lower POCP E .Toluene, for instance, has a POCP E of 45, which greatly exceeds the POCP E scenarios estimated for PCO and MPA.Taken together, results from this work and elsewhere suggest that MPA has a moderate reactivity in the troposphere.PCO is likely more reactive due to photolysis (Reaction R4).Even if significant uncertainties remain regarding PCO photochemistry, the use of both PCO and MPA would appear to have several advantages and to not adversely impact air quality when compared to traditional harmful and non-sustainable solvents.

Figure 1 .
Figure 1.Proposed pathways for PCO + OH (Reaction R1), where the H abstraction may proceed at the alpha CH 2 −H (Reaction R1a) or beta CH 2 −H (Reaction R1b, as proposed in the MCM).Similarly, two different sites for H abstraction are available for the reaction MPA + OH (Reaction R2).In the atmosphere, these radical products will react rapidly with oxygen to form peroxy radicals (RO 2 ).

3. 1
Figure2shows a typical exponential LIF signal decay, obtained from a Reaction (R1) experiment conducted at T = 460 K under pseudo-first-order conditions of [PCO] > 1000 × [OH].The inset in Fig.2summarises results from 12 such experiments, with pseudo-first-order rate parameter, B (Eq. 1) plot vs. corresponding[PCO].The slope from a linear fit to this dataset yields the second-order rate constant k 1 (460 K).Similar experiments were conducted at different temperatures and pressures and in the presence/absence of O 2 .Results from all such k 1 (T ) determinations are listed in Table 1 and displayed in Arrhenius format in Fig.3.A weighted average of the three values at close to ambient temperature yields k 1 (297 K) = (1.23 ± 0.09) × 10 −12 cm 3 molec.−1 s −1 , in good agreement with the solitary previous determination of k 1 (296 K) = (1.21± 0.05) × 10 −12 cm 3 molec.−1 s −1(Wallington and Kurylo, 1987a).Over the range of temperatures investigated in this work, these results were adequately represented by a two-parameter Arrhenius equation, as k 1 (297-485 K) = 2.14 × 10 −12 exp(−200/T ) cm 3 molec.−1 s −1 .Also displayed in Fig.3are results from this work on the well-characterised reaction between ethanol (CH 3 CH 2 OH) and OH (Reaction R3), many of which were obtained in back-to-back experiments alongside those on Reaction (R1).

Figure 4 .
Figure 4. Arrhenius plot for the reaction of MPA with OH (Reaction R2): the filled black square data points represent k 2 (T ) determinations from this work using N 2 bath gas and open black squares using air; the filled blue circles indicate literature data from Wallington et al. (1987); and the dashed black line represents a fit to all experimental data to yield k 2 (T ) = 5.0 × 10 −11 exp(−1743/T ) + 4.3 × 10 −13 exp(276/T ).The solid red line represents k 2 (T ) calculated from the SAR proposed by Jenkin et al. (2018).

Figure 5 .
Figure 5. Solution-phase spectra of MEK (black line) and PCO (blue line) obtained dilute in cyclohexane in this work, alongside the gas-phase spectrum of MEK (dashed red line) recorded by Martinez et al. (1992) and recommended by the IUPAC.Calculated gasphase spectrum of MEK (dotted-dashed orange line), calculated spectra of PCO for the gas phase (dotted-dashed light-blue line) and the solution phase (dotted dark-blue line).The inset (top right) displays Beer-Lambert plots obtained at 290 nm in this work, yielding a molar absorption coefficient ε 4(290 nm) = (20.6 ± 0.2) cm −1 M −1 for PCO and ε 5(290 nm) = (16.4± 0.2) cm −1 M −1 for MEK.

Table 1 .
Absolute determinations of k 1 , k 2 and k 3 from this work.