Quantum Chemical Investigation on Photochemical Reactions of

Article pubs.acs.org/JPCA

Quantum Chemical Investigation on Photochemical Reactions of Nonanoic Acids at Air−Water Interface Pin Xiao, Qian Wang, Wei-Hai Fang, and Ganglong Cui* Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China S Supporting Information *

ABSTRACT: Photoinduced chemical reactions of organic compounds at the marine boundary layer have recently attracted significant experimental attention because this kind of photoreactions has been proposed to have substantial impact on local new particle formation and their photoproducts could be a source of secondary organic aerosols. In this work, we have employed first-principles density functional theory method combined with cluster models to systematically explore photochemical reaction pathways of nonanoic acids (NAs) to form volatile saturated and unsaturated C9 and C8 aldehydes at air−water interfaces. On the basis of the results, we have found that the formation of C9 aldehydes is not initiated by intermolecular Norrish type II reaction between two NAs but by intramolecular T1 CO bond fission of NA generating acyl and hydroxyl radicals. Subsequently, saturated C9 aldehydes are formed through hydrogenation reaction of acyl radical by another intact NA. Following two dehydrogenation reactions, unsaturated C9 aldehydes are generated. In parallel, the pathway to C8 aldehydes is initiated by T1 CC bond fission of NA, which generates octyl and carboxyl radicals; then, an octanol is formed through recombination reaction of octyl with hydroxyl radical. In the following, two dehydrogenation reactions result into an enol intermediate from which saturated C8 aldehydes are produced via NA-assisted intermolecular hydrogen transfer. Finally, two dehydrogenation reactions generate unsaturated C8 aldehydes. In these reactions, water and NA molecules are found to play important roles. They significantly reduce relevant reaction barriers. Our work has also explored oxygenation reactions of NA with molecular oxygen and radical− radical dimerization reactions.

■

INTRODUCTION Air−water interface in particular near the marine boundary layer provides a unique environment where heat, momentum, water, as well as a large number of chemical species can exchange. In addition, neutral molecules, radicals, and ions meet, interact, and react at air−water interfaces to generate new molecules and eventually cause new particle formation (NPF). These exchanges and releases have a considerable and complex influence on global climate system, and regional air quality, and local secondary organic aerosol (SOA) formation. Because of these importances, physical and chemical processes at air− water interfaces have been intensively studied in the past decade.1−13 In recent years, photochemical reactions of organic compounds at the marine boundary layer, that is, air−sea interfaces, attracted a lot of experimental attention.14−20 Very recently, Rossignol et al.21 have employed ultraviolet−visible absorption and mass spectroscopy techniques to explore photochemical reactions at air−water interface coated with a monolayer of nonanoic acids. Upon direct illumination with actinic radiation, saturated and unsaturated aldehydes are seen in gas phase, while highly oxygenated products appear in © 2017 American Chemical Society

aqueous phase. In order to explain experimental observations, they have proposed two possible photoinduced reaction mechanisms to rationalize observed C9 aldehydes as shown in the left panel of Figure 1. The first pathway is initiated by an intramolecular Norrish type I reaction, that is, the CO bond fission, followed by a series of hydrogenation and dehydrogenation reactions by NA, hydroxyl, and carboxyl radicals; the second pathway starts from an intermolecular Norrish type II reaction forming a diol radical, which is followed by several dehydrogenation and hydrogen transfer reactions. Moreover, they have proposed two feasible reaction pathways for generating highly oxygenated products (see the right panel of Figure 1). In recent years, George and his co-workers have explored photosensitized reactions of nonanoic acids at air− water interfaces under different experimental conditions,20,22,23 as well as reactions with methyl vinyl ketone leading to significant amounts of isoprene.24 Received: April 2, 2017 Revised: May 16, 2017 Published: May 17, 2017 4253

DOI: 10.1021/acs.jpca.7b03123 J. Phys. Chem. A 2017, 121, 4253−4262

Article

The Journal of Physical Chemistry A

Figure 1. Experimentally proposed photochemcial reaction pathways (R1 and R2) for the formation of C9 aldehydes (left) and highly oxygenated products (right). NA· and NAH stand for CH3(CH2)7COO· and CH3(CH2)7C·(OH)OH radicals, respectively.

excited-state structures and energies of polyatomic molecules.30−35 All CASSCF and CASPT2 calculations are carried out using MOLCAS8.0.36,37 Electronic Structure Calculations in the T1 and S0 States. Geometries are optimized and characterized by frequency analysis calculations to be minima or transition states (TSs) using density functional theory (DFT) method38 with globally hybride exchange-correlation functional B3LYP with D3 dispersion correction scheme of Grimme.39−43 The 631G* basis set is used to expand total electronic wave function.44 Single-point energies of all B3LYP-D3 optimized structures are further refined by B3LYP-D3 method with a larger 6-311++G** basis set.39−43,45,46 Vertical excitation energies at the Franck−Condon point are also calculated using TD-B3LYP-D3/6-311++G** method.47 All DFT and TD-DFT calculations are performed using Gaussian 09 program package.48

Although experiments have given many useful information on photochemical reaction of nonanoic acids taking place at air−water interfaces, many atomistic mechanistic details remain elusive. For example, which reaction mechanism proposed by experiments for the formation of C9 aldehydes is more favorable energetically? Which elementary step is the ratelimiting one? How are saturated and unsaturated C8 aldehydes formed? What roles do water molecules play in these photochemical reactions? Obviously, this microscopic mechanistic information cannot be gained by experiments alone. In this regard, reliable electronic structure calculations play an irreplaceable role. Unfortunately, to date there is no theoretical work focusing on mechanistic issues of these photochemical reactions at air−water interfaces. In this work, we have employed density functional theory (DFT) method to systematically study photochemical reaction pathways of nonanoic acids that lead to saturated and unsaturated C8 and C9 aldehydes, as well as oxidation pathways to highly oxygenated products. This information certainly contributes important mechanistic insights for future experimental studies of similar photoreactions at air−water interfaces.

■

■

RESULTS AND DISCUSSION Spectroscopic Properties. In order to understand initial population of excited singlet states of NA molecule, we have explored its excited-state properties at the Franck−Condon point. Table 1 collects TD-B3LYP-D3 and B3LYP-D3 computed vertical excitation energies of NA and its two water complexes to the lowest two excited states, that is, S1 and T1. Electronic structure analysis shows that these two singlet and triplet states are of nπ* character as found in many similar

COMPUTATIONAL METHODS

Electronic Structure Calculations in the S1 State. Minimum-energy structures, minimum-energy reaction paths, and linearly interpolated internal coordinate (LIIC) paths in the S1 state are optimized using state-averaged complete active space self-consistent field (CASSCF) method in which equal weights are used. In all CASSCF calculations, an active space of 10 electrons in 8 orbitals is used (see Figure S1). Single-point energies of all CASSCF optimized structures are refined using multistate complete active space second-order perturbation theory (MS-CASPT2) approach.25,26 In MS-CASPT2 computations, the same active space, that is, 10 electrons in 8 orbitals is used; the Cholesky decomposition technique with unbiased auxiliary basis set is used for accurate two-electron integral evaluation;27 the ionization potential-electron affinity (IPEA) shift is set to zero;28 the imaginary shift technique (0.2 au) is employed to avoid intruder-state issue;29 This combined MSCASPT2//CASSCF computational protocol has been recently demonstrated to be able to give an accurate description for

Table 1. TD-B3LYP-D3/6-31G* and B3LYP-D3/6-31G* Computed Vertical Excitation Energies to the S1 and T1 States at the Franck-Condon Point and Adiabatic Excitation Energies of S1 and T1 Minima Relative to the S0 Minimum (Units in kcal/mol) NA NA-1W NA-2W

S1(1nπ*)a

S1(1nπ*)b

T1(3nπ*)a

T1(3nπ*)b

135.0 137.9 139.1

109.9 111.0 112.2

121.4 125.4 126.7

94.8 96.9 96.5

a

TD-B3LYP-D3/6-31G* and B3LYP-D3/6-31G* computed vertical excitation energies. bS1 and T1 states at the Franck-Condon point and adiabatic excitation energies. 4254


Article

The Journal of Physical Chemistry A saturated carbonyl compounds.31,34,49−57 In addition, it can be found that water hydration has small influence on vertical and adiabatic excitation energies to these two singlet and triplet states. The vertical and adiabatic excitation energies of waterhydrated NA molecules are merely blue-shifted a little compared with those of isolated NA molecule. For example, vertical excitation energies to the S1 [T1] state are estimated to be 135.0 [121.4], 137.9 [125.4], and 139.1 [126.7] kcal/mol for NA, NA-1W, and NA-2W, respectively (see complexes in Figure 2). This can be understood very well taking into account

Figure 3. CASPT2//CASSCF/6-31G* computed linearly interpolated internal coordinate path connecting both S0 and S1 minima. See text for discussion.

of these two photochemical reaction pathways and determined which one is more favorable in the view of energy. In the first reaction pathway, the formation of C9 aldehydes is proposed to start from the CO bond fission upon irradiation (see R1 in Figure 1). As discussed above, the excitation energy of 270 nm wavelength can populate the lowest excited S1 singlet state with very small probability due to small oscillator strength ( f = 0.002, CASPT2 level). In the S1 state, the direct CO bond fission that generates hydroxyl and carbonyl radicals is not efficient because the S1 CO bond fission needs to overcome a large barrier, 15.4 kcal/mol for isolated NA, and 14.5 kcal/mol for water-assisted NA-1W according to MSCASPT2//CASSCF computed S1 minimum-energy bondfission paths (see the top panels of Figure 4). Therefore, NA will first hop to the T1 state via the S1 → T1 intersystem crossing at the Franck−Condon region as already observed in many similar carbonyl compounds.34,49,52,54 In the T1 state, the CO bond fission of isolated NA from its T1 minimum 1* is still inefficient because the T1 energy of the TS(1−2) transition state, that is, 105.2 kcal/mol at B3LYPD3 level, is close to available total energy (exp. 106 kcal/mol; red lines in Figure 5). However, when a water molecule is bound with NA forming a complex 1a, the T1 energy of the TS (1a-2) transition state for the T1 CO bond fission is significantly reduced, which is 91.5 kcal/mol at B3LYP-D3 level (blue lines in Figure 5). Thus, with the assistance of one water molecule, the T1 CO bond fission becomes efficient and finally produces hydroxyl and carbonyl radicals, that is, OH and 2 radicals in Figure 5. Both water and NA molecules can react with radical species 2 to generate molecule−radical complexes 4′ and 4. However, water-assisted radical reaction, from 3′ via TS(3′-4′) to 4′, cannot compete with that assisted by NA itself, from 3 via TS(3−4) to 4. Their barriers at TS(3′-4′) and TS(3−4) are computed to be 14.9 and 12.8 kcal/mol at B3LYP-D3 level, respectively. In addition, the latter reaction is less endothermic thermodynamically. Subsequently, dehydrogenation reaction of molecule−radical complexes 4′ and 4 takes place immediately. Among them, the reaction associated with molecule−radical complex 4 is much easier than that associated with molecule− radical complex 4′ because the energy of the TS(4−5) transition state is about 10 kcal/mol lower than that of the

Figure 2. NA molecule and its two water-hydrated complexes NA-1W and NA-2W. Also shown are some selected hydrogen bond lengths (in Å) of the S0 equilibrium structures optimized by B3LYP-D3/6-31G* method. See Supporting Information for Cartesian coordinates.

that both singlet and triplet states originate from local n → π* electronic excitation of the carbonyl group. Hydrogen-bonding interaction due to water hydration to somewhat extent inhibits electronic transition of the lone-pair electron of the O atom of the carbonyl group. The S1 state cannot be directly populated at the Franck− Condon point by 270 nm wavelength, about 106 kcal/mol, because of large S0 → S1 vertical excitation energy, 135.0 kcal/ mol in Table 1. However, one can find that this excitation energy is close to the S0 → S1 adiabatic excitation energy, 109.9 kcal/mol at TD-B3LYP-D3 level in Table 1. In addition, at room temperature the NA molecule has sufficient thermal energy to fluctuate so as to deviate from its S0 minimum. Figure 3 shows linearly interpolated internal coordinate (LIIC) path connecting the Franck−Condon point of NA and its S1 minimum. It can be found that when NA molecule moves toward its S1 minimum-energy structure, which is mainly driven by the CO stretching mode and the improper torsional motion with respect to the central C atom of the COOH group, the S1−S0 energy gap will become small enough. As a result, this 270 nm wavelength can populate the S1 state of NA with a small probability. Furthermore, the S1 population is also enhanced by neat or concentrated solutions of NA, as suggested by recent experiments.21 Initiation and Formation of Unsaturated C9 Aldehydes. Rossignol et al. have experimentally proposed two photochemical reaction mechanisms for the formation of C9 aldehydes.21 Next, we have separately examined the feasibility 4255


Article


Figure 4. CASPT2//CASSCF/6-31G* computed minimum-energy reaction paths (kcal/mol) for the CO (top) and CC (bottom) bond fissions of (left) NA and (right) NA-1W in the S1 state. See text for discussion.

reaction is rather easier because of small barrier in the T1 state, which is computed to be 4.1 kcal/mol at the TS(7−8) transition state. Nevertheless, the following hydrogen migration reactions are to certain extent inhibited as demonstrated by their relatively high barriers, which are all more than 16 kcal/ mol at B3LYP-D3 level (see Figure 7). Moreover, the energies at all these transition states are more than 99 kcal/mol at B3LYP-D3 level, even 129.0 kcal/mol at TS(8−9) for 1,2 hydrogen migration. With these reasons taken together, it is not difficult to find that the R2 path is unfavorable compared with the R1 path in the formation of unsaturated C9 aldehydes. Therefore, we have refrained us from further exploration of subsequent reactions of this R2 path. Initiation and Formation of Unsaturated C8 Aldehydes. In experiments of Rossignol and co-workers, C8 aldehydes have also been observed using mass spectroscopy techniques.21 But, corresponding photoinduced reaction pathways are not proposed. In this work, we have figured out an energetically favorable reaction path producing these saturated and unsaturated C8 aldehydes, as shown in Figures 8 and 9. The formation of C8 aldehydes are initiated by the CC bond dissociation. Such bond fission is not allowed in the S1 state because it is associated with a large barrier for NA (see the bottom-left panel of Figure 4). Although water hydration can reduce this barrier to some extent, this reaction still cannot

TS(4′-5) one (Figure 5). The resultant radical 5 can further evolve into final unsaturated C9 aldehyde P1 via the second dehydrogenation reaction. Different from the first dehydrogenation reaction, the reaction initiated by the OH radical is much easier than that by the 1 radical generated through intermolecular radical reaction between 2 radical and NA molecule. Figure 5 shows computed reaction path to form unsaturated C9 aldehyde from triplet NA molecule in which the most favorable one is in blue. In addition to forming α, β unsaturated aldehyde, we have explored reaction pathways to generate other kinds of unsaturated aldehydes from species 4, such as P1−X (X = 1−6). As shown in Figure 6, all these reactions are thermodynamically very exothermic and are also associated with small barriers. Among them, forming P1−6 is associated with the largest barrier, but only 7.0 kcal/mol at B3LYP-D3 level, which can be surmounted very easily concerning available total energy. Of course, among these different unsaturated aldehydes, P1 is the most stable one because of its conjugated CC and CO double bonds. Besides the above reaction path of forming unsaturated aldehydes, we have examined experimentally proposed another reaction path (see R2 in Figure 1). Different from the R1 path, the R2 path is initiated by intermolecular hydrogen extraction reaction between two NA molecules. This intermolecular 4256


Article


Figure 5. B3LYP-D3/6-311++G** computed reaction pathways to form saturated and unsaturated C9 aldehydes from triplet NA molecule initiated by the CO bond fission. See text for discussion.

Figure 6. B3LYP-D3/6-311++G** computed reaction pathways to form different unsaturated C9 aldehydes from saturated C9 aldehyde, that is, species 4. See text for discussion.

4257


Article


Figure 7. B3LYP-D3/6-311++G** computed reaction pathways to form unsaturated C9 aldehydes from initial NA molecules upon irradiation. See text for discussion.

Figure 8. B3LYP-D3/6-311++G** computed photochemical reaction pathways to form saturated C8 aldehydes from triplet NA molecules initiated by the CC bond fission. See text for discussion.

reaction process needs to overcome a barrier of 59.9 kcal/mol and thus is forbidden in the view of energy. In the second one, the radical species 10 will recombine with OH radical to generate 11. Then, from neutral molecule 11, two dehydrogenation reactions happen sequentially. These two reactions are initiated by OH radicals and are thus very facile energetically. The generated molecule 13 will convert into saturated C8 aldehyde P2 through three reaction paths, among which the

happen. Therefore, the CC bond dissociation should occur in the T1 state as a consequence of the S1 → T1 intersystem crossing at the Franck−Condon region. The CC bondfission barrier in the T1 state is computed to be 6.7 kcal/mol, which can be overcome by photon energy of 270 nm wavelength. The generated radical species 10 has two reaction channels. In the first one, it reacts with the HOC(O) radical to generate species 11 with the release of carbon oxide. This 4258


Article


Figure 9. B3LYP-D3/6-311++G** computed reaction pathways to form different unsaturated C8 aldehydes from saturated C8 aldehyde P2. See text for discussion.

Figure 10. Photochemical pathways to generate saturated and unsaturated C8 and C9 aldehydes suggested based on present studies.

NA-assisted one has a relatively smaller barrier, 13.6 kcal/mol at B3LYP-D3/6-311++G** (the other two have 56.2 and 44.8 kcal/mol; see Figure 8). This small barrier can be gotten through very easily by the “hot” system because the energy at the TS(13a-P2) transition state is significantly lower than the

T1 energy of initial triplet NA species, that is, 1* by more than 100 kcal/mol. In this reaction step, NA molecule plays an important role in reducing the barrier. Of course, the generated saturated C8 aldehyde P2 can be further converted into unsaturated C8 aldehydes through two 4259


Article


and radical−radical coupling reactions. Our present results evidently demonstrate that the pathway initiated by the CO bond fission of NA in the T1 state, and not by the intermolecular hydrogen transfer between NAs, is energetically more favorable. Alternatively, we have computationally proposed a feasible photochemical pathway for the formation of unsaturated and saturated C8 aldehydes, which is initiated by the T1 CC bond fission. Moreover, we have found that water and NA molecules play important roles in photochemical reactions of NAs at some key reaction steps. They are directly involved in water- and NA-assisted hydrogen transfer reactions in such a way as to reduce significantly related reaction barriers. We have found that all oxygenation reactions of NA radicals with molecular oxygen and radical−radical coupling reactions can easily happen. Finally, our present work opens the way for exploring photochemical reactions at air−water interfaces and could contribute useful knowledge to understanding local new particle formation near the marine boundary layer.20,22,23

sequential dehydrogenation reactions. As discussed above, both hydroxyl and carboxyl radicals can be used to extract hydrogen atoms of saturated C9 aldehydes to form unsaturated C9 ones (see Figure 5). But, it is obvious that abstraction reaction by carboxyl radical is energetically more favorable than that by hydroxyl radical. Therefore, we have merely explored dehydrogenation reactions of the C8 aldehyde P2 by carboxyl radical. As shown in Figure 9, all reactions are energetically allowed because of small barriers. In addition to the formation of C9 and C8 aldehydes, oxygenation and radical coupling reactions have been observed in experiments (see the right panel of Figure 1).21 In order to get a comprehensive understanding of photochemical reactions of NAs at air−water interfaces, we have also explored all possible pathways of oxygenation and radical coupling reactions (see Supporting Information for details). In oxygenation reactions, hydroxyl radical first abstracts a hydrogen atom from different sites of NA molecule, resulting into eight 1_X (X = 1−8) radicals and one carboxylic radical (Figure S10). All these abstraction reactions are associated with small barriers and can proceed easily. Moreover, they are thermodynamically exothermic, more than 20 kcal/mol at B3LYP-D3 level. In the following, oxygen molecules can be readily added to these eight radicals in a barrierless way. Finally, we have calculated reaction energy changes of totally more than 35 radical−radical coupling reactions and found that all radical recombination reactions are extremely exothermic thermodynamically (more than 80 kcal/ mol at B3LYP-D3 level; Figure S11). Photochemical Mechanism. On the basis of our results, we have proposed two photochemical reaction pathways generating unsaturated and saturated C8 and C9 aldehydes, respectively. Upon irradiation, NA molecules are first excited to their S1 state, which is followed by an S1 → T1 intersystem crossing process. In the T1 state, Norrish type I CO and C C bond fissions become feasible. The left panel of Figure 10 shows the most favorable reaction path to produce unsaturated and saturated C9 aldehydes. The T1 CO bond fission generates hydroxyl and acyl radicals in the same time. This step is accelerated by water molecules near air−water interfaces. The newly generated acyl radical abstracts a hydrogen atom from nearby NA molecule forming a C9 aldehyde, that is, nonanal. From this saturated aldehyde, two sequential dehydrogenation reactions lead to the formation of unsaturated C9 aldehydes. On the other hand, if the CC bond is cleaved in the T1 state, as shown in the right panel of Figure 10, the generated octyl radical will recombine with hydroxyl radical to form an octanol. Through two consecutive dehydrogenation reactions initiated by hydroxyl radical, an enol is formed. From this species, an intramolecular hydrogen transfer will result into a saturated C8 aldehyde, that is, capryl aldehyde. Of course, this reaction step is also accelerated by water molecules, in particular by NA molecules, which make the reaction barrier lower by more than 11.4 and 42.6 kcal/mol at B3LYP-D3/6-311++G** level. Similarly, final unsaturated C8 aldehydes are generated by two steps of dehydrogenation reactions.

■

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpca.7b03123. Schematic structures, oxygenation and radical−radical coupling reaction profiles and discussion, and Cartesian coordinates of optimized structures (PDF)

■

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. ORCID

Wei-Hai Fang: 0000-0002-1668-465X Ganglong Cui: 0000-0002-9752-1659 Notes

The authors declare no competing financial interest.

■

ACKNOWLEDGMENTS This work was supported by the National Key Research & Development Program of China (2016YFC0202600) and the National Natural Science Foundation of China (21522302, 21520102005, and 21421003).

■

REFERENCES

(1) George, C.; Ammann, M.; D’Anna, B.; Donaldson, D. J.; Nizkorodov, S. A. Heterogeneous Photochemistry in the Atmosphere. Chem. Rev. 2015, 115, 4218−4258. (2) Donaldson, D. J.; Valsaraj, K. T. Adsorption and Reaction of Trace Gas-Phase Organic Compounds on Atmospheric Water Film Surfaces: A Critical Review. Environ. Sci. Technol. 2010, 44, 865−873. (3) Donaldson, D. J.; George, C.; Vaida, V. Red Sky at Night: LongWavelength Photochemistry in the Atmosphere. Environ. Sci. Technol. 2010, 44, 5321−5326. (4) Donaldson, D. J.; George, C. Sea-Surface Chemistry and Its Impact on the Marine Boundary Layer. Environ. Sci. Technol. 2012, 46, 10385−10389. (5) Donaldson, D. J.; George, C. Introduction to the Focus Issue on Marine Boundary Layer: Ocean Atmosphere Interactions Processes. Environ. Sci. Technol. 2012, 46, 10383−10384. (6) Herrmann, H.; Schaefer, T.; Tilgner, A.; Styler, S. A.; Weller, C.; Teich, M.; Otto, T. Tropospheric Aqueous-Phase Chemistry: Kinetics, Mechanisms, and Its Coupling to a Changing Gas Phase. Chem. Rev. 2015, 115, 4259−4334.

■

CONCLUSIONS We have employed ab initio electronic structure calculations to explore photochemical reactions of nonanoic acids (NAs) at air−water interfaces using cluster models assisted by water and NA molecules.21 First, we have examined two experimentally proposed photochemical reaction pathways for the formation of unsaturated and saturated C9 aldehydes as well as oxygenation 4260


Article


(27) Aquilante, F.; Lindh, R.; Pedersen, T. B. Unbiased Auxiliary Basis Sets for Accurate Two-Electron Integral Approximations. J. Chem. Phys. 2007, 127, 114107−114713. (28) Ghigo, G.; Roos, B. O.; Malmqvist, P.-Å. A Modified Definition of the Zeroth-Order Hamiltonian in Multiconfigurational Perturbation Theory (CASPT2). Chem. Phys. Lett. 2004, 396, 142−149. (29) Forsberg, N.; Malmqvist, P.-Å. Multiconfiguration Perturbation Theory with Imaginary Level Shift. Chem. Phys. Lett. 1997, 274, 196− 204. (30) Cui, G. L.; Guan, P.-J.; Fang, W.-H. Photoinduced Proton Transfer and Isomerization in a Hydrogen-Bonded Aromatic Azo Compound: A CASPT2//CASSCF Study. J. Phys. Chem. A 2014, 118, 4732−4739. (31) Chang, X.-P.; Fang, Q.; Cui, G. L. Mechanistic Photodecarboxylation of Pyruvic Acid: Excited-State Proton Transfer and Three-State Intersection. J. Chem. Phys. 2014, 141, 154311−154318. (32) Chang, X.-P.; Li, C.-X.; Xie, B.-B.; Cui, G. L. Photoprotection Mechanism of p-Methoxy Methylcinnamate: A CASPT2 Study. J. Phys. Chem. A 2015, 119, 11488−11497. (33) Guan, P.-J.; Cui, G. L.; Fang, Q. Computational Photochemistry of the Azobenzene Scaffold of Sudan I and Orange II Dyes: ExcitedState Proton Transfer and Deactivation via Conical Intersections. ChemPhysChem 2015, 16, 805−811. (34) Xia, S.-H.; Liu, X.-Y.; Fang, Q.; Cui, G. L. Excited-State RingOpening Mechanism of Cyclic Ketones: A MS-CASPT2//CASSCF Study. J. Phys. Chem. A 2015, 119, 3569−3576. (35) Xie, B.-B.; Li, C.-X.; Cui, G. L.; Fang, Q. Excited-State Proton Transfer and Decay in a Hydrogen-Bonded Oxazole System: a MSCASPT2//CASSCF Study. Chin. J. Chem. Phys. 2016, 29, 38−46. (36) Karlström, G.; Lindh, R.; Malmqvist, P.-Å.; Roos, B. O.; Ryde, U.; Veryazov, V.; Widmark, P.-O.; Cossi, M.; Schimmelpfennig, B.; Neogrady, P.; et al. MOLCAS: A Program Package for Computational Chemistry. Comput. Mater. Sci. 2003, 28, 222−239. (37) Aquilante, F.; De Vico, L.; Ferré, N.; Ghigo, G.; Malmqvist, P.Å.; Neogrády, P.; Pedersen, T. B.; Pitoňaḱ , M.; Reiher, M.; Roos, B. O.; et al. MOLCAS 7: The Next Generation. J. Comput. Chem. 2010, 31, 224−247. (38) Parr, R. G.; Yang, W. T. Density-Functional Theory of Atoms and Molecules; Oxford University Press: New York, 1994. (39) Vosko, S. H.; Wilk, L.; Nusair, M. Accurate Spin-Dependent Electron Liquid Correlation Energies for Local Spin Density Calculations: A Critical Analysis. Can. J. Phys. 1980, 58, 1200−1211. (40) Lee, C.; Yang, W. T.; Parr, R. G. Development of the ColleSalvetti Correlation-Energy Formula into a Functional of the Electron Density. Phys. Rev. B: Condens. Matter Mater. Phys. 1988, 37, 785−789. (41) Becke, A. D. Density-Functional Exchange-Energy Approximation with Correct Asymptotic Behavior. Phys. Rev. A: At., Mol., Opt. Phys. 1988, 38, 3098−3100. (42) Becke, A. D. A New Mixing of Hartree-Fock and Local DensityFunctional Theories. J. Chem. Phys. 1993, 98, 1372−1377. (43) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104−154122. (44) Hariharan, P. C.; Pople, J. A. The Influence of Polarization Functions on Molecular Orbital Hydrogenation Energies. Theor. Chem. Acc. 1973, 28, 213−222. (45) Krishnan, R.; Binkley, J. S.; Seeger, R.; Pople, J. A. SelfConsistent Molecular Orbital Methods. XX. A Basis Set for Correlated Wave Functions. J. Chem. Phys. 1980, 72, 650−654. (46) Clark, T.; Chandrasekhar, J.; Spitznagel, G. W.; Schleyer, P. V. R. Efficient Diffuse Function-Augmented Basis Sets for Anion Calculations. III. The 3-21+G Basis Set for First-Row Elements, LiF. J. Comput. Chem. 1983, 4, 294−301. (47) Marques, M. A. L.; Ullrich, C. A.; Nogueira, F.; Rubio, A.; Burke, K.; Gross, E. K. U. Time-dependent Density Functional Theory; Springer: Berlin, 2006. (48) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheesem, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.;

(7) Donaldson, D. J.; Vaida, V. The Influence of Organic Films at the Air-Aqueous Boundary on Atmospheric Processes. Chem. Rev. 2006, 106, 1445−1461. (8) Zhu, C.; Kumar, M.; Zhong, J.; Li, L.; Francisco, J. S.; Zeng, X. C. New Mechanistic Pathways for Criegee-Water Chemistry at the Air/ Water Interface. J. Am. Chem. Soc. 2016, 138, 11164−11169. (9) Zhu, C.; Kais, S.; Zeng, X. C.; Francisco, J. S.; Gladich, I. Interfaces Select Specific Stereochemical Conformations: The Isomerization of Glyoxal at the Liquid Water Interface. J. Am. Chem. Soc. 2017, 139, 27−30. (10) Xu, J.; Finlayson-Pitts, B. J.; Gerber, R. B. Proton Transfer in Mixed Clusters of Methanesulfonic Acid, Methylamine, and Oxalic Acid: Implications for Atmospheric Particle Formation. J. Phys. Chem. A 2017, 121, 2377−2385. (11) Miller, Y.; Thomas, J. L.; Kemp, D. D.; Finlayson-Pitts, B. J.; Gordon, M. S.; Tobias, D. J.; Gerber, R. B. Structure of Large NitrateWater Clusters at Ambient Temperatures: Simulations with Effective Fragment Potentials and Force Fields with Implications for Atmospheric Chemistry. J. Phys. Chem. A 2009, 113, 12805−12814. (12) Liu, Z.; Wu, L. Y.; Wang, T. H.; Ge, M. F.; Wang, W. G. Uptake of Methacrolein into Aqueous Solutions of Sulfuric Acid and Hydrogen Peroxide. J. Phys. Chem. A 2012, 116, 437−442. (13) Peng, C.; Jing, B.; Guo, Y. C.; Zhang, Y. H.; Ge, M. F. Hygroscopic Behavior of Multicomponent Aerosols Involving NaCl and Dicarboxylic Acids. J. Phys. Chem. A 2016, 120, 1029−1038. (14) Reeser, D. I.; George, C.; Donaldson, D. J. Photooxidation of Halides by Chlorophyll at the Air-Salt Water Interface. J. Phys. Chem. A 2009, 113, 8591−8595. (15) Reeser, D. I.; Jammoul, A.; Clifford, D.; Brigante, M.; D’Anna, B.; George, C.; Donaldson, D. J. Photoenhanced Reaction of Ozone with Chlorophyll at the Seawater Surface. J. Phys. Chem. C 2009, 113, 2071−2077. (16) Mmereki, B. T.; Donaldson, D. J.; Gilman, J. B.; Eliason, T. L.; Vaida, V. Kinetics and Products of the Reaction of Gas-Phase Ozone with Anthracene Adsorbed at the Air-Aqueous Interface. Atmos. Environ. 2004, 38, 6091−6103. (17) Mmereki, B. T.; Donaldson, D. J. Direct Observation of the Kinetics of an Atmospherically Important Reaction at the Air-Aqueous Interface. J. Phys. Chem. A 2003, 107, 11038−11042. (18) Mmereki, B. T.; Donaldson, D. J. Laser Induced Fluorescence of Pyrene at an Organic Coated Air-Water Interface. Phys. Chem. Chem. Phys. 2002, 4, 4186−4191. (19) Enami, S.; Sakamoto, Y. OH-Radical Oxidation of SurfaceActive cis-Pinonic Acid at the Air-Water Interface. J. Phys. Chem. A 2016, 120, 3578−3587. (20) Bernard, F.; Ciuraru, R.; Boréave, A.; George, C. Photosensitized Formation of Secondary Organic Aerosols above the Air/ Water Interface. Environ. Sci. Technol. 2016, 50, 8678−8686. (21) Rossignol, S.; Tinel, L.; Bianco, A.; Passananti, M.; Brigante, M.; Donaldson, D. J.; George, C. Atmospheric Photochemistry at a Fatty Acid-Coated Air-Water Interface. Science 2016, 353, 699−702. (22) Ciuraru, R.; Fine, L.; van Pinxteren, M.; D’Anna, B.; Herrmann, H.; George, C. Photosensitized Production of Functionalized and Unsaturated Organic Compounds at the Air-Sea Interface. Sci. Rep. 2015, 5, 12741−12750. (23) Tinel, L.; Rossignol, S.; Bianco, A.; Passananti, M.; Perrier, S.; Wang, X.; Brigante, M.; Donaldson, D. J.; George, C. Mechanistic Insights on the Photosensitized Chemistry of a Fatty Acid at the Air/ Water Interface. Environ. Sci. Technol. 2016, 50, 11041−11048. (24) Ciuraru, R.; Fine, L.; van Pinxteren, M.; D’Anna, B.; Herrmann, H.; George, C. Unravelling New Processes at Interfaces: Photochemical Isoprene Production at the Sea Surface. Environ. Sci. Technol. 2015, 49, 13199−13205. (25) Andersson, K.; Malmqvist, P.-Å.; Roos, B. O.; Sadlej, A. J.; Wolinski, K. Second-Order Perturbation Theory with a CASSCF Reference Function. J. Phys. Chem. 1990, 94, 5483−5488. (26) Andersson, K.; Malmqvist, P.-Å.; Roos, B. O. Second-Order Perturbation Theory with a Complete Active Space Self-Consistent Field Reference Function. J. Chem. Phys. 1992, 96, 1218−1226. 4261


Article

The Journal of Physical Chemistry A Petersson, G. A. et al. Gaussian 09, Revision A.02; Gaussian, Inc.: Wallingford, CT, 2009. (49) Cui, G. L.; Li, Q.-S.; Zhang, F.; Fang, W.-H.; Yu, J.-G. Combined CASSCF and MRCI Study on Photoinduced Dissociation and Isomerization of Acryloyl Chloride. J. Phys. Chem. A 2006, 110, 11839−11846. (50) Cui, G. L.; Zhang, F.; Fang, W. Insights into the Mechanistic Photodissociation of Methyl Formate. J. Chem. Phys. 2010, 132, 034306−034310. (51) Cui, G. L.; Ai, Y. J.; Fang, W.-H. Conical Intersection Is Responsible for the Fluorescence Disappearance below 365 nm in Cyclopropanone. J. Phys. Chem. A 2010, 114, 730−734. (52) Cui, G. L.; Fang, W. Mechanistic Photodissociation of Glycolaldehyde: Insights from ab Initio and RRKM Calculations. ChemPhysChem 2011, 12, 1351−1357. (53) Cui, G. L.; Fang, W.-H. Ab Initio Based Surface-Hopping Dynamics Study on Ultrafast Internal Conversion in Cyclopropanone. J. Phys. Chem. A 2011, 115, 1547−1555. (54) Cui, G. L.; Sun, Z.; Fang, W. Adiabatic and Nonadiabatic Bond Cleavages in Norrish Type I Reaction. J. Phys. Chem. A 2011, 115, 10146−10153. (55) Shemesh, D.; Lan, Z. G.; Gerber, R. B. Dynamics of triplet-state photochemistry of pentanal: mechanisms of Norrish I, Norrish II, and H abstraction reactions. J. Phys. Chem. A 2013, 117, 11711−11724. (56) Shemesh, D.; Blair, S. L.; Nizkorodov, S. A.; Gerber, R. B. Photochemistry of aldehyde clusters: cross-molecular versus unimolecular reaction dynamics. Phys. Chem. Chem. Phys. 2014, 16, 23861− 23868. (57) Xia, S.-H.; Liu, X.-Y.; Fang, Q.; Cui, G. L. Photodissociation Dynamics of CH3C(O)SH in Argon Matrix: A QM/MM Nonadiabatic Dynamics Simulation. J. Chem. Phys. 2015, 143, 194303− 194313.

4262


Quantum Chemical Investigation on Photochemical Reactions of

Recommend Documents