Aging of α-Pinene Secondary Organic Aerosol by Hydroxyl Radicals in the Aqueous Phase: Kinetics and Products.

The reaction of hydroxyl radicals (OH) with a water-soluble fraction of the α-pinene secondary organic aerosol (SOA) was investigated using liquid chromatography coupled with negative electrospray ionization mass spectrometry. The SOA was generated by the dark ozonolysis of α-pinene, extracted into the water, and subjected to chemical aging by the OH. Bimolecular reaction rate coefficients (kOH) for the oxidation of terpenoic acids by the OH were measured using the relative rate method. The unaged SOA was dominated by the cyclobutyl-ring-retaining compounds, primarily cis-pinonic, cis-pinic, and hydroxy-pinonic acids. Aqueous oxidation by the OH resulted in the removal of early-stage products and dimers, including well-known oligomers with MW = 358 and 368 Da. Furthermore, a 2- to 5-fold increase in the concentration of cyclobutyl-ring-opening products was observed, including terpenylic and diaterpenylic acids and diaterpenylic acid acetate as well as some of the newly identified OH aging markers. At the same time, results obtained from the kinetic box model showed a high degree of SOA fragmentation following the reaction with the OH, which indicates that non-radical reactions occurring during the evaporation of water likely contribute to the high yields of terpenoic aqSOAs reported previously. The estimated atmospheric lifetimes showed that in clouds, terpenoic acids react with the OH exclusively in the aqueous phase. Aqueous OH aging of the α-pinene SOA results in a 10% increase of the average O/C ratio and a 3-fold decrease in the average kOH value, which is likely to affect the cloud condensation nuclei activity of the aqSOA formed after the evaporation of water.

The Malaprade reaction mechanism for ethylene glycol oxidation by periodic acid based on density functional theory (DFT).

A general mechanism of the Malaprade oxidative carbon-carbon bond cleavage reaction of α-glycol in the presence of periodic acid has been proposed on the basis of density functional theory (DFT) computations. Ethylene glycol and periodic acid, both in their neutral forms, have been studied as noble substrate representatives in model reactions. The proposed reaction mechanism has been constructed based on and compared with previously published experimental kinetic, spectroscopic and temperature and pH-dependent studies. In the presented theoretical mechanistic considerations, four alternative molecular transformations have been analyzed from thermodynamic and kinetic points of view. Theoretically, the predicted activation energy barriers have been compared with experimental ones published elsewhere. The presented minimum energy pathway (MEP) unveiled the shape and conformation of the intermediate and transition state structures. The three-step reaction process involves the formation of a seven-membered quasi-ring assisted by an intramolecular hydrogen-bond intermediate structure forming one I-O bond (IC1_B), a cyclic ester intermediate forming two I-O bonds (IC2_C) and the final products formed at the two very last stages (HIO3, water and two formaldehyde molecules). The computed and energetically the most favourable reaction landscape proposed in this work uniforms and refines the mechanistic proposition given by Criegee for Malaprade type of reactions and further gives a detailed molecular understanding of the reaction rate and atomic connections en route the transformation. The molecular geometries of all stationary points (intermediate and transition state structures) lying on the potential energy hypersurface have been optimized at the four alternative DFT levels under the solvation model based on the density approximation: B3PW91, CAM-B3LYP, BMK, ωB97XD. The 6-311+G(2d,p) basis set for C, O, and H atoms and both the full (DGDZVP) and Ahlrichs-Weigend1 (def2-TZVP) basis sets for iodine atoms were used during the computations.

Alkylation of Nitroarenes via Vicarious Nucleophilic Substitution - Experimental and DFT Mechanistic Studies.

Alkylation of nitroarenes via Vicarious Nucleophilic Substitution (VNS) was tested experimentally and modelled with DFT calculations. Mechanistic studies reveal intrinsic differences between reactions of archetypal carbanion precursor PhSO2 CH2 Cl, and alkyl phenyl sulfones, in which benzenesulfinate acts as a leaving group. Accordingly, for the latter precursors steric hindrance develops at the ß-elimination step, that raises energy barrier and results in the formation of byproducts.

Exploring Intra- and Intermolecular Interactions in Selected N-Oxides-The Role of Hydrogen Bonds.

Intra- and intermolecular interactions have been explored in selected N-oxide derivatives: 2-(N,N-dimethylamino-N-oxymethyl)-4,6-dimethylphenyl (1) and 5,5'-dibromo-3-diethylaminomethyl-2,2'-biphenol N-oxide (2). Both compounds possess intramolecular hydrogen bonding, which is classified as moderate in 1 and strong in 2, and resonance-assisted in both cases. Density Functional Theory (DFT) in its classical formulation as well as Time-Dependent extension (TD-DFT) were employed to study proton transfer phenomena. The simulations were performed in the gas phase and with implicit and explicit solvation models. The obtained structures of the studied N-oxides were compared with experimental data available. The proton reaction path was investigated using scan with an optimization method, and water molecule reorientation in the monohydrate of 1 was found upon the proton scan progress. It was found that spontaneous proton transfer phenomenon cannot occur in the electronic ground state of the compound 1. An opposite situation was noticed for the compound 2. The changes of nucleophilicity and electrophilicity upon the bridged proton migration were analyzed on the basis of Fukui functions in the case of 1. The interaction energy decomposition of dimers and microsolvation models was investigated using Symmetry-Adapted Perturbation Theory (SAPT). The simulations were performed in both phases to introduce polar environment influence on the interaction energies. The SAPT study showed rather minor role of induction in the formation of homodimers. However, it is worth noticing that the same induction term is responsible for the preference of water molecules' interaction with N-oxide hydrogen bond acceptor atoms in the microsolvation study. The Natural Bond Orbital (NBO) analysis was performed for the complexes with water to investigate the charge flow upon the polar environment introduction. Finally, the TD-DFT was applied for isolated molecules as well as for microsolvation models showing that the presence of solvent affects excited states, especially when the N-oxide acceptor atom is microsolvated.

Non-Covalent Forces in Naphthazarin-Cooperativity or Competition in the Light of Theoretical Approaches.

Non-covalent interactions responsible for molecular features and self-assembly in Naphthazarin C polymorph were investigated on the basis of diverse theoretical approaches: Density Functional Theory (DFT), Diffusion Quantum Monte Carlo (DQMC), Symmetry-Adapted Perturbation Theory (SAPT) and Car-Parrinello Molecular Dynamics (CPMD). The proton reaction paths in the intramolecular hydrogen bridges were studied. Two potential energy minima were found indicating that the proton transfer phenomena occur in the electronic ground state. Diffusion Quantum Monte Carlo (DQMC) and other levels of theory including Coupled Cluster (CC) employment enabled an accurate inspection of Potential Energy Surface (PES) and revealed the energy barrier for the proton transfer. The structure and reactivity evolution associated with the proton transfer were investigated using Harmonic Oscillator Model of Aromaticity - HOMA index, Fukui functions and Atoms In Molecules (AIM) theory. The energy partitioning in the studied dimers was carried out based on Symmetry-Adapted Perturbation Theory (SAPT) indicating that dispersive forces are dominant in the structure stabilization. The CPMD simulations were performed at 60 K and 300 K in vacuo and in the crystalline phase. The temperature influence on the bridged protons dynamics was studied and showed that the proton transfer phenomena were not observed at 60 K, but the frequent events were noticed at 300 K in both studied phases. The spectroscopic signatures derived from the CPMD were computed using Fourier transformation of autocorrelation function of atomic velocity for the whole molecule and bridged protons. The computed gas-phase IR spectra showed two regions with OH absorption that covers frequencies from 2500 cm-1 to 2800 cm-1 at 60 K and from 2350 cm-1 to 3250 cm-1 at 300 K for both bridged protons. In comparison, the solid state computed IR spectra revealed the environmental influence on the vibrational features. For each of them absorption regions were found between 2700-3100 cm-1 and 2400-2850 cm-1 at 60 K and 2300-3300 cm-1 and 2300-3200 cm-1 at 300 K respectively. Therefore, the CPMD study results indicated that there is a cooperation of intramolecular hydrogen bonds in Naphthazarin molecule.

Structural Characterization of Lactone-Containing MW 212 Organosulfates Originating from Isoprene Oxidation in Ambient Fine Aerosol.

Isoprene (C5H8) is the main non-methane hydrocarbon emitted into the global atmosphere. Despite intense research, atmospheric transformations of isoprene leading to secondary organic aerosol (SOA) are still not fully understood, including its multiphase chemical reactions. Herein, we report on the detailed structural characterization of atmospherically relevant isoprene-derived organosulfates (OSs) with a molecular weight (MW) of 212 (C5H8SO7), which are abundantly present in both ambient fine aerosol (PM2.5) and laboratory-generated isoprene SOA. The results obtained from smog chamber-generated isoprene SOA and aqueous-phase laboratory experiments coupled to the S(IV)-autooxidation chemistry of isoprene, 3-methyl-2(5H)-furanone, and 4-methyl-2(5H)-furanone, allowed us for the first time to fully elucidate the isomeric structures of the MW 212 OSs. By applying liquid chromatography interfaced to electrospray ionization high-resolution mass spectrometry, we firmly confirmed six positional isomers of the MW 212 OSs in PM2.5 collected from different sites in Europe and the United States. Our results also show that despite the low solubility of isoprene in water, aqueous-phase or multiphase chemistry can play an important role in the formation of OSs from isoprene. Possible formation mechanisms for the MW 212 OSs are also tentatively proposed.

How Do Aromatic Nitro Compounds React with Nucleophiles? Theoretical Description Using Aromaticity, Nucleophilicity and Electrophilicity Indices.

In this study, we present a complete description of the addition of a model nucleophile to the nitroaromatic ring in positions occupied either by hydrogen (the first step of the SNAr-H reaction) or a leaving group (SNAr-X reaction) using theoretical parameters including aromaticity (HOMA), electrophilicity and nucleophilicity indices. It was shown both experimentally and by our calculations, including kinetic isotope effect modeling, that the addition of a nucleophile to the electron-deficient aromatic ring is the rate limiting step of both SNAr-X and SNAr-H reactions when the fast transformation of σH-adduct into the products is possible due to the specific reaction conditions, so this is the most important step of the entire reaction. The results described in this paper are helpful for better understanding of the subtle factors controlling the reaction direction and rate.

The activation of carbon dioxide by first row transition metals (Sc-Zn).

The activation of CO2 by chloride-tagged first-row transition metal anions [ClM]- (M = Sc-Zn), was examined by mass spectrometry, quantum chemical calculations, and statistical analysis. The direct formation of [ClM(CO2)]- complexes was demonstrated in the reaction between [ClM]- and neutral CO2. In addition, the reverse reaction was investigated by energy-variable collisionally induced dissociation (CID) of the corresponding [ClM(CO2)]- anions generated in-source. Five different mono- and bi-dentate binding motifs present in the ion/CO2 complexes were identified by quantum chemical calculations and the relative stability of each of these isomers was established and analyzed for all first-row transition metals based on the experimental and theoretical ion/molecule binding energies. It was found that the early first row transition metals form strong covalent bonds with the neutral CO2 molecule, while the late ones and in particular copper and zinc are weakly bonded. Using simple valence bond Lewis diagrams, the different binding motifs and their relative stabilities across the first row were described using multi-configurational self consistent field (MSCSCF) wavefunctions in a quantitative manner based on the electronic structure of the individual metals. This analysis provides an explanation for the change of the most favorite bonding motif of the transition metals with CO2 along the 1st transition metal row. The nature of the activated CO2 complex and the relationship between its stability and other structural and spectral properties was also analyzed by Principal Component Analysis (PCA) and artificial neural networks.

How Does Nucleophilic Aromatic Substitution Really Proceed in Nitroarenes? Computational Prediction and Experimental Verification.

The aim of this paper is to present a correct and complete mechanistic picture of nucleophilic substitution in nitroarenes based on the results obtained by theoretical calculations and experimental observations coming from numerous publications, reviews, and monographs. This work gives the theoretical background to the very well documented experimentally yet still ignored observations that the addition of nucleophiles to halo nitroarenes resulting in the formation of σ(H) adducts, which under proper reaction conditions can be transformed into the product of the SNArH reaction, is faster than the competing process of addition to the carbon atom bearing a nucleofugal group (usually a halogen atom) resulting in the "classic" SNAr reaction. Only when the σ(H) adduct cannot be transformed into the SNArH reaction product, SNAr reaction is observed.

Experimental versus Calculated Proton Affinities for Aromatic Carboxylic Acid Anions and Related Phenide Ions.

Herein, we present the comparison of a large set of experimentally measured proton affinity (PA) values for 65 aromatic carboxylate anions with the values calculated by using selected popular DFT (B3LYP, PBE0, and M05-2X) and composite [G3(MP2), G4(MP2)] quantum chemistry methods. The root-mean-square error (RMSE) values for the chosen methods are RMSEPBE0 =1.7, RMSEB3LYP =4.6, RMSEM05-2X =6.6, RMSEG3MP2 =6.3, RMSEG4MP2 =4.5 kJ mol(-1) . In the second part of the study, 82 PA values for substituted phenide ions and a few heteroaromatic anions were calculated. Again, very good agreement between the calculated and experimental values has been observed: RMSEPBE0 =1.9, RMSEB3LYP =4.5, RMSEM05-2X =6.3, RMSEG3MP2 =4.9, RMSEG4MP2 =5.5 kJ mol(-1) . Our results show that, for medium-sized carboxylate anions, all tested methods give reliable results and, surprisingly, much more computationally demanding composite methods do not perform significantly better than the time-efficient DFT methods.

Gas-phase reactions of methyl thiocyanate with aliphatic carbanions - A mass spectrometry and computational study.

RATIONALE: Methyl thiocyanate, like other organic thiocyanates, is a molecule with many electrophilic reactive sites and it has many synthetic applications. For better understanding of the intrinsic reactivity of alkyl thiocyanates against nucleophiles it was important to study gas-phase reactions of methyl thiocyanate with carbanions differing by structure and proton affinity values. METHODS: All experiments were performed using a modified API 365 triple quadrupole mass spectrometer equipped with a TurboIonSpray electrospray ionization (ESI) source. Carbanions were generated in the ESI source by decarboxylation of the respective carboxylic acid anions. Methyl thiocyanate was delivered as a vapor with nitrogen used as a collision gas to the collision cell where the reactions take place. RESULTS: Mass spectra recorded for the gas-phase reactions of five aliphatic carbanions with methyl thiocyanate showed a variety of product ions formed via different reaction mechanisms, depending on the structure and proton affinity of the carbanion. The pathways considered are: SN 2 nucleophilic substitution, cyanophilic reaction, thiophilic reaction and proton transfer, followed in some instances by subsequent transformations. The proposed reaction pathways are supported by density functional theory (DFT) calculations. CONCLUSIONS: Our preliminary experiments showed that mass spectrometry together with quantum chemical calculations is a good tool for studying gas-phase reactions of alkyl thiocyanates with carbanions. In the gas phase all four theoretically possible products can be observed and their formation can be rationalized by the results of the modelling of the reaction energy profiles. Copyright © 2016 John Wiley & Sons, Ltd.

Competition between nucleophilic substitution of halogen (SN Ar) versus substitution of hydrogen (SN ArH)-a mass spectrometry and computational study.

The mechanism of intramolecular gas-phase reactions of N-(2-X-5-nitrophenyl)-N-methylacetamide carbanions (X=H, F, Cl) has been studied using negative ion electrospray mass spectrometry ((-)ESI-MS) technique and modelled computationally. It was proven that all three anions form cyclic σ(H) adducts, which undergo elimination of water. In the case of X=F, formation of the σ(F) adduct, leading to SN Ar reaction, was a competing process. This is the first proof that also in the gas phase formation of σ(H) adduct proceeds faster than σ(X) adduct and only when X=F, rates of these two processes are comparable. The experimental results are in full agreement with quantum chemical calculations.

Molecular reorganization of selected quinoline derivatives in the ground and excited states­Investigations via static DFT.

Quinoline derivatives are interesting objects to study internal reorganizations due to the observed excited-state-induced intramolecular proton transfer (ESIPT). Here, we report on computations for selected 12 quinoline derivatives possessing three kinds of intramolecular hydrogen bonds. Density functional theory was employed for the current investigations. The metric and electronic structure simulations were performed for the ground state and first excited singlet and triplet states. The computed potential energy profiles do not show a spontaneous proton transfer in the ground state, whereas excited states exhibit this phenomenon. Atoms in Molecules (AIM) theory was applied to study the nature of hydrogen bonding, whereas Harmonic Oscillator Model of aromaticity index (HOMA) provided data of aromaticity evolution as a derivative of the bridge proton position. The AIM-based topological analysis confirmed the presence of the intramolecular hydrogen bonding. In addition, using the theory, we were able to provide a quantitative illustration of bonding transformation: from covalent to the hydrogen. On the basis of HOMA analysis, we showed that the aromaticity of both rings is dependent on the location of the bridge proton. Further, the computed results were compared with experimental data available. Finally, ESIPT occurrence was compared for the three investigated kinds of hydrogen bridges, and competition between two bridges in one molecule was studied.

Dissociation constants of relevant secondary organic aerosol components in the atmosphere.

The presented studies focus on measuring the determination of the acidity constant (pKa) of relevant secondary organic aerosol components. For our research, we selected important oxidation products (mainly carboxylic acids) of the most abundant terpene compounds, such as α-pinene, ß-pinene, ß-caryophyllene, and δ-3-carene. The research covered the synthesis and determination of the acidity constant of selected compounds. We used three methods to measure the acidity constant, i.e., 1H NMR titration, pH-metric titration, Bates-Schwarzenbach spectrophotometric method. Moreover, the pKa values were calculated with Marvin 21.17.0 software to compare the experimentally derived values with those calculated from the chemical structure. pKa values measured with 1H NMR titration ranged from 3.51 ± 0.01 for terebic acid to 5.18 ± 0.06 for ß-norcaryophyllonic acid. Moreover, the data determined by the 1H NMR method revealed a good correlation with the data obtained with the commonly used potentiometric and UV-spectroscopic methods (R2 = 0.92). In contrast, the comparison with in silico results exhibits a relatively low correlation (R2Marvin = 0.66). We found that most of the values calculated with the Marvin Program are lower than experimental values obtained with pH-metric titration with an average difference of 0.44 pKa units. For di- and tricarboxylic acids, we obtained two and three pKa values, respectively. A good correlation with the literature values was observed, for example, Howell and Fisher (1958) used pH-metric titration and measured pKa1 and pKa2 to be 4.48 and 5.48, while our results are 4.24 ± 0.10 and 5.40 ± 0.02, respectively.

Physicochemical Properties of Terebic Acid, MBTCA, Diaterpenylic Acid Acetate, and Pinanediol as Relevant α-Pinene Oxidation Products.

The physicochemical properties and the synthesis of four α-pinene oxidation products, terebic acid, 3-methyl-1,2,3-butanetricarboxylic acid (MBTCA), diaterpenylic acid acetate (DTAA), and pinanediol, are presented in this study. The physicochemical properties encompass thermal properties, solubility in water, and dissociation constant (pK a) for the investigated compounds. It was found that terebic acid exhibits a relatively high melting temperature of 449.29 K, whereas pinanediol revealed a low melting temperature of 329.26 K. The solubility in water was determined with the dynamic method and the experimental results were correlated using three different mathematical models: Wilson, NRTL, and UNIQUAC equations. The results of the correlation indicate that the Wilson equation appears to work the best for terebic acid and pinanediol. The calculated standard deviation was for 3.79 for terebic acid and 1.25 for pinanediol. In contrast, UNIQUAC was the best mathematical model for DTAA and MBTCA. The calculated standard deviation was 0.57 for DTAA and 2.21 for MBTCA. The measured water solubility increased in the following order: pinanediol > DTAA ≥ MBTCA > terebic acid, which affects their multiphase aging chemistry in the atmosphere. Moreover, acidity constants (pK a) at 298, 303, and 308 K were determined for DTAA with the Bates-Schwarzenbach spectrophotometric method. The pK a values obtained at 298, 303, and 308 K were found to be 3.76, 3.85, and 3.88, respectively.

Correction to "Physicochemical Properties of Terebic Acid, MBTCA, Diaterpenylic Acid Acetate, and Pinanediol as Relevant α-Pinene Oxidation Products".

[This corrects the article DOI: 10.1021/acsomega.9b04231.].

Physicochemical Properties of Pinic, Pinonic, Norpinic, and Norpinonic Acids as Relevant α-Pinene Oxidation Products.

Here, the study is focused on the synthesis and determination of physicochemical properties of four α-pinene secondary organic aerosol (SOA) products: cis-pinic acid, cis-pinonic acid, cis-norpinic acid, and cis-norpinonic acid. These encompass their thermal properties, solid-liquid phase equilibria, and dissociation constant (pKa). Thermal properties, including the melting temperature, enthalpy of fusion, temperature, and enthalpy of the phase transitions, were measured with the differential scanning calorimetry technique. These SOA components exhibit relatively high melting temperatures from 364.32 K for cis-pinic acid to 440.68 K for cis-norpinic acid. The enthalpies of fusion vary from 14.75 kJ·mol-1 for cis-norpinic acid to 30.35 kJ·mol-1 for cis-pinonic acid. The solubility in water was determined with the dynamic method (solid-liquid phase equilibria method), and then experimental results were interpreted and correlated using three different mathematical models: Wilson, non-random two-liquid model, and universal quasichemical equations. The results of the correlation indicate that the Wilson equation appears to work the best for all investigated compounds, giving rise to the lowest value of a standard deviation. cis-Norpinic acid and cis-pinic acid (dicarboxylic acids) show better solubility in the aqueous solution than cis-norpinonic acid and cis-pinonic acid (monocarboxylic acids), which affect the multiphase chemistry of α-pinene SOA processes. For cis-pinonic acid and cis-norpinonic acid, also pH-profile solubility was determined. The intrinsic solubility (S0) for cis-norpinonic acid was measured to be 0.05 mmol·dm-3, while for cis-pinonic acid, it was found to be 0.043 mmol·dm-3. The acidity constants (pKa) at 298 and 310 K using the Bates-Schwarzenbach spectrophotometric method were determined. The pKa values at 298.15 K for cis-norpinonic acid and cis-pinonic acid were found to be 4.56 and 5.19, respectively, whereas at 310.15 K, pKa values were found to be -4.76 and 5.25, respectively.