Estimating the saturation vapor pressures of isoprene oxidation products C₅H₁₂O₆ and C₅H₁₀O₆ using COSMO-RS

• Main Authors: T. Kurtén, N. Hyttinen, E. L. D'Ambro, J. Thornton, N. L. Prisle
• Format: Article
• Language: English
• Published: Copernicus Publications 2018-12-01
• Series: Atmospheric Chemistry and Physics
• Online Access: https://www.atmos-chem-phys.net/18/17589/2018/acp-18-17589-2018.pdf

<p>We have used COSMO-RS (the conductor-like screening model for real solvents), as implemented in the COSMOtherm program, to compute the saturation vapor pressures at 298&thinsp;K of two photo-oxidation products of isoprene: the dihydroxy dihydroperoxide <span class="inline-formula">C₅H₁₂O₆</span>, and the dihydroperoxy hydroxy aldehyde, <span class="inline-formula">C₅H₁₀O₆</span>. The predicted saturation vapor pressures were significantly higher (by up to a factor of 1000) than recent experimental results, very likely due to the overestimation of the effects of intramolecular hydrogen bonds, which tend to increase saturation vapor pressures by stabilizing molecules in the gas phase relative to the liquid. Modifying the hydrogen bond enthalpy parameter used by COSMOtherm can improve the agreement with experimental results – however the optimal parameter value is likely to be system-specific. Alternatively, vapor pressure predictions can be substantially improved (to within a factor of 5 of the experimental values for the two systems studied here) by selecting only conformers with a minimum number of intramolecular hydrogen bonds. The computed saturation vapor pressures were very sensitive to the details of the conformational sampling approach, with the default scheme implemented in the COSMOconf program proving insufficient for the task, for example by predicting significant differences between enantiomers, which should have identical physical properties. Even after exhaustive conformational sampling, COSMOtherm predicts significant differences in saturation vapor pressures between both structural isomers and diastereomers. For <span class="inline-formula">C₅H₁₂O₆</span>, predicted differences in <span class="inline-formula"><i>p</i>_sat</span> between structural isomers are up to 2 orders of magnitude, and differences between stereoisomers are up to a factor of 20 – though these differences are very likely exaggerated by the overestimation of the effect of intramolecular H-bonds. For <span class="inline-formula">C₅H₁₀O₆</span>, the maximum predicted differences between the three studied structural isomers and their diastereomer pairs are around a factor of 8 and a factor of 2, respectively, when only conformers lacking intramolecular hydrogen bonds are included in the calculations. In future studies of saturation vapor pressures of polyfunctional atmospheric oxidation products using COSMOtherm, we recommend first performing thorough conformational sampling and subsequently selecting conformers with a minimal number of intramolecular H-bonds.</p>