Zeyou
Pan
ab,
Andras
Bodi
a,
Jeroen A.
van Bokhoven
ab and
Patrick
Hemberger
*a
aZeyou Pan, Andras Bodi, Jeroen A. van Bokhoven and Patrick Hemberger, Paul Scherrer Institute, 5232 Villigen, Switzerland. E-mail: patrick.hemberger@psi.ch
bZeyou Pan and Jeroen A. van Bokhoven, Institute for Chemical and Bioengineering, Department of Chemistry and Applied Biosciences, ETH Zurich, 8093 Zurich, Switzerland
First published on 26th January 2022
We report the absolute photoionization cross section (PICS) of fulvenone and 2-carbonyl cyclohexadienone, two crucial ketene intermediates in lignin pyrolysis, combustion and organic synthesis. Both species were generated in situ by pyrolyzing salicylamide and dectected via imaging photoelectron photoion coincidence spectroscopy. In a deamination reaction, salicylamide loses ammonia yielding 2-carbonyl cyclohexadienone, a ketoketene, which further decarbonylates at higher pyrolysis temperatures to form fulvenone. We recorded the threshold photoelectron spectrum of the ketoketene and assigned the ground state (+2A′′ ← 1A′) and excited state (Ã+2A′ ← 1A′) bands with the help of Franck–Condon simulations. Adiabatic ionization energies are 8.35 ± 0.01 and 9.19 ± 0.01 eV. In a minor reaction channel, the ketoketene isomerizes to benzpropiolactone, which decomposes subsequently to benzyne by CO2 loss. Potential energy surface and RRKM rate constant calculations agree with our experimental observations that the decarbonylation to fulvenone outcompetes the decarboxylation to benzyne by almost two orders of magnitude. The absolute PICS of fulvenone at 10.48 eV was determined to be 18.8 ± 3.8 Mb using NH3 as a calibrant. The PICS of 2-carbonyl cyclohexadienone was found to be 21.5 ± 8.6 Mb at 9 eV. Our PICS measument will enable the quantification of reactive ketenes in lignin valorization and combustion processes using photoionization techniques and provide advanced mechanistic and kinetics insights to aid the bottom-up optimization of such processes.
Due to its high reactivity, fulvenone evades detection using standard chemical analysis tools, such as GC/MS and NMR, which is the reason why fulvenone was only observed using photoionization mass spectrometry (PIMS), photoelectron spectroscopy (PES), matrix infrared spectroscopy (IR) and PEPICO detection.7,11–13 Most recently, Genossar et al. recorded the IR spectrum of fulvenone produced by salicylaldehyde pyrolysis.13 This ketene was also synthesized in situ via pyrolysis of lignin model compounds and characterized by photoelectron spectroscopy.6,12,14 We have measured the ms-TPES and photoionization spectrum of fulvenone and simulated transitions from the neutral ground 2A2 into both the ground +2A2 and the first excited Ã+2B1 cation states in the Franck–Condon approximation.15 Fulvenone has, thus, been characterized by different techniques, which helps to trace this crucial intermediate in complex reaction pathways to gain an advanced mechanistic understanding of combustion and lignin valorization processes. However, to obtain reliable kinetics and analytical data, quantification is mandatory. In photoionization measurements, the photoionization cross section (PICS) is an important measure, which relates the mass spectral signal to a concentration in a reaction mixture and can be defined using a reference signal as follows:16
(1) |
The decomposition of salicylamide 1 is initiated by deamination (R1) yielding the first intermediate, 6-carbonyl-2,4-cyclohexadien-1-one or 2-carbonyl cyclohexadienone (2, C7H4O2), which we call ketoketene 2 from here on. 2 forms fulvenone 3 by sequential decarbonylation (R1). Both CO and NH3 are produced together with fulvenone in a 1:1:1 ratio. The latter can be utilized as reference, because of its well-known PICS and ionization energy close to that of fulvenone 3.29
In this contribution, we set out to measure the absolute PICS of the fulvenone ketene utilizing salicylamide as a precursor. We will investigate the salicylamide pyrolysis mechanism to identify possible side reactions and to determine the selectivity of the ketene formation. Therefore, the ms-TPES of the ketoketene 2, the precursor of fulvenone 3, is analyzed in detail, also because of its relevance in the lignin pyrolysis chemistry.6 Potential energy surface calculations and computational rate constants will aid the analysis of the experimental findings on salicylamide pyrolysis mechanism to obtain ultimately the PICSs of fulvenone 3 and ketoketene 2.
Upon increasing the reactor temperature, the salicylamide 1 (m/z 137) and ketoketene 2 peaks decreased while the fulvenone 3 and NH3 signals increased. Between 1020 and 1140 K, fulvenone dominates the mass spectrum at 10.5 eV. Except for CO, no side products were observed. As CO is ionized above 14 eV, it is not seen in Fig. 1. As reported by Ormond et al.,6 ketoketene 2 forms fulvenone 3 by decarbonylation and subsequent ring contraction. Therefore, 2 must be fully converted into 3 for the PICS measurement, since only then are ammonia and fulvenone produced in a known, 1:1 ratio, required for ammonia to be used as reference in eqn (1). At 1020 K reactor temperature, 2 showed only a small peak at m/z 120, which disappears completely once the temperature is increased to 1140 K. The C5H44 and C6H4 isomers’ 5 peaks, however, rose simultaneously at this temperature. Note that C5H4 isomers originate from the decarbonylation of fulvenone 3 (R2)6 while C6H4 isomers were yielded by CO2 loss from 6-carbonyl-2,4-cyclohexadien-1-one 2 in (R3):
Thus, different reaction temperatures between 1020 and 1140 K were thoroughly investigated and the signals of m/z 64, 76, 92 and 120 were integrated and presented in Fig. S3 (ESI†). We also noticed that dissociative ionization of ketoketene 2 (m/z 120) yields fulvenone cations (m/z 92), too, as shown in Fig. S4 (ESI†). By increasing the reactor temperature, the concentration of salicylamide 1 in the molecular beam was lowered and the dissociative ionization contribution to fulvenone was minimized. Therefore, the photoionization spectrum (PIS) of fulvenone for the PICS measurement (see below) was recorded at 1050 K reactor temperature and 338 K sample container temperature (determined by the thermocouple outside of sample container) at a selectivity of 87%. Still, the formation of C6H4 isomers (m/z 64), such as benzyne, could not be fully supressed and needs to be further investigated, as it is in direct competition with the fulvenone 3 generation from the ketoketene 2. In fact, ketoketene was also observed during the pyrolysis of lignin model compounds, such as salicylaldehyde, and may play an important role in this process.6 Thus, in the next chapter we focus on the spectroscopic characterization of the ketoketene (2, m/z 120) and its fate at higher pyrolysis temperature to yield 3 and 5.6
Method | (+2A′′) | (Ã+2A′) |
---|---|---|
Adiabatic ionization energy (AIE)/eV | ||
Experiment: | ||
ms-TPES | 8.35 | 9.19 |
Theory: | ||
B3LYP/6-311+G(d,p) | 8.27 | 9.02 |
CBS-QB3 | 8.30 | — |
CBS-APNO | 8.29 | — |
G3 | 8.42 | — |
G4 | 8.36 | — |
CCSD/6-311+G(d,p) | 8.01 | — |
MP2/6-311+G(d,p) | 9.03 | — |
ms-TPES (band max.) | 8.35 | 9.34 |
PES11 (band max.) | 8.43 | 9.38 |
EOM-IP-CCSD/cc-pVTZ (VIE) | 8.45 | 9.53 |
Furthermore, the EOM-IP-CCSD/cc-pVTZ calculated vertical ionization energy (VIE) of 9.53 eV (see Table 1) verifies our assignment for the transitions between 9.19 and 9.64 eV computationally. Besides the ms-TPES, the PIS is shown in Fig. 2a. It exhibits a step-like onset close to the first ionization energy and linearly increases afterwards but exhibits plateaus in the 8.8–9.1 eV and 9.7–10.0 eV photon energy ranges. The constant PI signal between 8.8–9.1 eV is mirrored by the slightly decreasing ms-TPES, indicative for dropping Franck–Condon factors. To further understand the electronic structure, we carried out B3LYP/6-311+G(d,p) calculations on 2. Fig. 2b shows the three frontier orbitals of the Cs symmetric neutral ketene (1A′) and Table S2 (ESI†) summarizes the most important geometry changes upon photoionization. According to Koopmans' theorem, removing an electron from the highest occupied molecular orbital (HOMO) yields the +2A′′ cation state. The HOMO has bonding components along the C2C4 double bond as well as at the C6C5 position. At the carbonyl C3O and ketene function C7O, the HOMO possesses antibonding character. Once the ketoketene 2 is ionized, the electron density decreases, leading to an increase of bond lengths (C2C4, C6C5, and C1C7), as an electron is removed from a formally binding orbital component. In contrast, the C1–C2 and C4–C6 bond lengths are shortened, because the antibonding character of these bonds decreases. A similar pattern was also found in fulvenone 2 ionization.15 In general, the bond angles are only slightly affected upon ionization into the ground state cation.
ms-TPE spectroscopy is a unique tool to detect reactive intermediates and to selectively identify isomers in harsh environments.45,46 Its full potential can only be realized if both the electronic structure and the vibrational features in the ms-TPE spectrum are fully understood. The latter can be accurately simulated by applying the Franck–Condon principle in the double harmonic approximation at the optimized geometries of neutral and cation state, while also including the Duschinsky rotation.47 The FC predicted spectrum of the +2A′′ (red trace) and Ã+2A′ (green trace) cation states of ketoketene 2 are depicted in Fig. 2a. The 0–0 transition to the +2A′′ cation state is located at 8.35 eV, followed by excitations in several vibrational modes active upon photoionization. The lowest transition energy is 50 meV (403 cm−1) above the origin and is assigned to the v7 ring deformation mode at a theoretical value of 451 cm−1 at the B3LYP/6-311+G(d,p) level of theory. At 8.47 eV, the experimental transition energy is 1371 cm−1 which compares well to the theoretical value of 1456 cm−1 (v25) of a ring deformation mode. The third band with a shift of 0.27 eV (2177 cm−1) with respect to the origin is a CO stretching vibration with a computed transition energy of v29 = 2296 cm−1. Furthermore, we observe combination bands of v7 and v10 and v24, the latter two are assigned to further deformation ring modes.
The active vibrations mirror the geometry change upon ionization as elicited by the removal of a HOMO electron, as discussed above. Schulz and Schweig assigned the band above 8.8 eV to the excited state of 6-carbonyl-2,4-cyclohexadien-1-one cation at a IE of 9.38 eV,11 which agrees with the most-prominent band at 9.34 eV in our ms-TPES. Our computational results for the Ã+2A′ state of the keteoketene 2 cation, which finds a VIE at 9.53 eV at EOM-IP-CCSD/cc-pVTZ level of theory (see Table 1), also confirms this assignment. While the band at 9.19 eV is the origin transition into the Ã+2A′ state, the band at 9.25 eV is assigned to the v6 mode at a transition energy of 450 cm−1 at B3LYP/6-311+G(d,p) level of theory. This stretching vibration at C1–C3 and C4–C6 agrees with the change of the bond angles at C2C4–C6 and C4–C6C5 upon transition into the Ã+2A′ state and is expected from the removal of an electron from the HOMO−1. Besides the v6 vibration, the v25 mode (1422 cm−1) is active and shows a ring deformation mode dominated by the stretching of C1–C3 and C3O bonds, thereby leading to a significant change in C6C5–C3 and C5–C3–C1 angles as well as the bond length in C1–C3, C1C7, C3–C5 and C3O (see Table S2, ESI†). Above 9.4 eV, there are at least two further vibrational features assigned to combination bands of the mentioned transitions. The FC simulation of the Ã+2A′ state shows good agreement with the experimental spectrum below 9.45 eV. However, the intensity of the second band at 9.34 eV is not reproduced well. EOM-IP-CCSD/cc-pVDZ calculations for the +2A′′ and Ã+2A′ states confirmed the harmonic DFT/TDDFT FC factors. We also investigated other isomers which may contribute to the ms-TPES above 9.2 eV, which are depicted in Fig. 3 and Fig. S5 (ESI†). From the seven computationally investigated isomers, only benzpropiolactone (7-oxabicyclo[4.2.0]octa-1(6),2,4-triene-8-one) 6 is almost isoenergetic to the ketoketene 2, while the other isomers (7–12, see Fig. 3) are more than 2 eV less stable and may require much higher barriers to be formed. In addition, isomers 7–12 are likely to decompose to C2 or C3 carbon chains at higher pyrolysis temperatures, which were not observed experimentally. Furthermore, the calculated FC envelopes and ionization energies of isomers 7–12 do not match to the features observed in the ms-TPES of m/z 120 (see Fig. S5, ESI†). Yet, when investigating the ionization of benzpropiolactone 6 computationally, we found that the C–O bond of the 4-membered ring opens upon ionization to form the ketoketene 2+ cation. This large geometry change leads to unfavourable FC factors at the calculated adiabatic ionization energy of 8.31 eV (G4). The EOM-IP-CCSD/cc-pVTZ computed VIE, on the other hand, is located at 9.19 eV. Thus the differences between the experimental spectrum and the FC simulations (Fig. 2) may perhaps be explained by intensity borrowing from isomer 6, which was also observed in the matrix IR studies by Chapman et al.41,48
Fig. 3 Potential energy surface of the salicylamide pyrolysis (a) and comparison of RRKM rate constant of the CO and CO2 loss reactions (b). |
(2) |
(3) |
The photoion signal of ammonia was normalized to the PICS from Xia et al.29 and the fulvenone PIS was corrected with the same factor to yield the PICS (see Fig. 4 and Table S3, ESI†). By applying this procedure we obtain a PICS of ammonia from our PI data of 0.88 Mb at 10.35 eV, which is within the 10% error bar from the literature.29,50,51 Due to the nature of the pyrolysis source, reactive species are produced vibrationally hot, which leads to the appearance of hot and sequence bands below the adiabatic ionization energy of 8.25 eV.15 Above the IE, the PICS monotonically increases until a plateau region between 9.3 and 9.6 eV is reached, while it decreases afterwards. In the following, we discuss the error bars of our fulvenone PICS measurement. The mass discrimination in the effusive molecular beam is close to unity, and integration effects, such as the 13C isotopologue, as well as subtraction of the false coincidence signal and the photon flux correction are considered minor errors, which account for less than 5% of the absolute error. At the best measurement conditions for 3, we found about 6% unconverted ketoketene (2, at m/z 120), while the decarboxylation product benzyne accounts for up to 2% of the total ion signal (see Fig. S3, ESI†).
In addition, fulvenone decomposes to C5H4 species (see Fig. S2, ESI†), which contributes further 5% to the uncertainty, yielding a total of ca. 11% by error propagation. Considering the temperature dependence of the ion signal due to the contribution of hot and sequence bands, reproducibility effects and the uncertainty in the cross section measurements of ammonia (10%),29,50,51 leads us to a total error bar of up to 20%. The latter is indicated as grey shaded area in Fig. 4 and agrees with the typical error bar for PICS measurements in the literature.16,20
At 840 K reactor temperature, fulvenone (m/z 92) formation does not yet play a role and the ammonia signal can be used to determine the PICS of the ketoketene 2 species at m/z 120, which is depicted in Fig. 5 (red curve). Several assumptions had to be taken into account to determine this physical quantity: Dissociative ionization of the ketoketene 2 in a decarbonylation reaction above 9.5 eV produces m/z 92. The corresponding ion signal is confirmed to be due to dissociative ionization by velocity map ion images (VMIs), which are sensitive to kinetic energy release, as depicted in Fig. S4 (ESI†). Since the DPI onset of m/z 120 to m/z 92 is similar to the ionization energy of the calibrant, ammonia, both species (m/z 120 and 92) were added together to get the total PICS (red curve Fig. 5). In addition, fulvenone is also produced via pyrolysis, and was subtracted from the m/z 92 signal, utilizing the fulvenone PI curve at full conversion (Fig. 4), yielding solely the DPI signal (see Fig. 5, blue trace). It is intriguing to compare now the PICS of fulvenone 3 (Fig. 4) with the one of m/z 120 in Fig. 5 (red curve), which is about a factor of two lower at 9 eV. This may be owing to contributions of the bicyclic benzpropiolactone 6 isomer, which may account for 40% of the 120 amu population in the reaction mixture at 840 K pyrolysis temperature, due to it being only 30 meV less stable than 2 (Fig. 3 a) and the presence of a low-lying isomerization transition state . Owing to the large geometry change upon ionization, as discussed in the TPES section (AIE = 8.34 eV vs. VIE = 9.19 eV) the FC factors of isomer 6 are likely negligible below 8.9–9.0 eV photon energy. Thus, only ketoketene 2 contributes significantly to the PICS below 9 eV. However, the reference ammonia is formed in a 1:1 ratio to the analytes 2 + 6, which must be considered in eqn (1) to derive the cross section for 2 in the energy range where 6 is not expected to contribute. By scaling the PICS signal up by 1.67, according to relative abundance of 2 in 2 + 6 at 840 K, we derive an effective PICS for the ketoketene 2 at low photon energies (black curve in Fig. 5). Due to these considerations, we increased the error bars of the m/z 120 PICS measurements conservatively to about 40%, also including the discussed errors from the fulvenone measurements. We would like to point out that in all hot reactive environments (combustion & catalysis) both isomers 2 and 6 are likely to be present due to their isoenergeticity and a separation may not be possible at higher photon energies. Our experimental PICS data at m/z 120 provides a lumped PICS of 2 and 6 (red curve in Fig. 5) above 9 eV, which may still enable a semi-quantitative analysis of both isomers in reaction mixtures.
Upon increasing the temperature, ketoketene 2 decarbonylates to yield fulvenone ketene 3. We have found a second decomposition channel of the ketoketene 2, which leads to the benzyne 5 side product by decarboxylation. Potential energy surface calculations starting from salicylamide show that deamination proceeds over a barrier of less than 1.8 eV to yield the ketoketene 2. Thereafter, parallel CO2 or CO loss, yields benzyne 5 or fulvenone 3, respectively. Although the mass spectra suggest that fulvenone is the dominant product, the barrier to CO2 loss is lower than that to CO loss. Nevertheless, RRKM rate constant calculations show that the decarbonylation over a looser transition state outcompetes decarboxylation by more than one order of magnitude in the energy range where rates become commensurate with the residence time in the microreactor and in agreement with experimental observations. Equipped with this knowledge, we could determine a reliable photon-energy-dependent ionization cross section of fulvenone ketene 3 for the first time. Corresponding to different measurement uncertainties, a typical error bar of 20% was estimated. At a typical laser photon energy of 10.48 eV (3 × 355 nm) the PICS of fulvenone 3 was found to be 18.8 ± 3.8 Mb. The cross section of the ketoketene 2 could also be obtained at lower pyrolysis temperatures. However, due to non-suppressible side reactions such as dissociative ionization of the ketoketene 2, fulvenone 3 and benzpropiolactone 6 formation, we had to assume more conservative error bars for m/z 120. Nonetheless, the PICS of 2-carbonyl cyclohexadienone (ketoketene) could be estimated below the VIE of 6, by considering the thermal equilibrium between isomers 2 and 6. We have determined an effective photoionization cross section of 2 of 21.5 Mb at 9 eV, which compares with 21.9 Mb of fulvenone PICS.
Our PICS measurements enable the quantification of highly reactive ketenes in lignin valorization and combustion processes using photoionization mass spectrometric tools. This will make kinetics data accessible to determine the relative contribution of parallel reaction pathways in the lignin catalytic pyrolysis reaction mechanism, which will in turn aid the bottom up optimization of this process.46
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1cp05206c |
This journal is © the Owner Societies 2022 |