Ion Imaging Study of NO3 Radical Photodissociation Dynamics: Characterization of Multiple Reaction Pathways

Rotational Distributions and Imaging of Singlet O2 Following Spin-Forbidden Photodissociation of O3

We report REMPI spectra and velocity-mapped ion images of the O2(a1Δg) and (b1Σg+) fragments arising from the spin-forbidden photodissociation of O3 near 320 and 330 nm. The O2(a1Δg, v = 0) REMPI spectrum following a 320 nm dissociation shows enhanced peak intensity for the odd rotational states relative to the even states, which is the opposite of the trend observed by Gunthardt et al. ( J. Chem. Phys. 2019, 151, 224302) for spin-allowed dissociation at 266 nm but is consistent with the couplings between the B state and 3A' and 3A″ states calculated by Grebenshchikov and Rosenwaks ( J. Phys. Chem. A 2010, 114, 9809-9819). There are no significant differences between the ion image angular distributions of fragments in odd and even rotational states, which indicates a cold distribution of O3 and supports the explanation that the alternation in peak intensities results from a difference in the couplings. Quantitative analysis of the image angular distributions was limited due to the single laser polarization geometry accessible in one-color experiments. Radial distributions of the 320 nm images indicate a broad rotational distribution, evidenced in bimodal speed distributions with peaks corresponding to both high (j = 35-43) and low (j = 17-20) rotational states. The REMPI spectrum of O2(a1Δg) near 330 nm was collected, and while quantitative population analysis is difficult because of the perturbed resonant state, the spectrum clearly supports a broad rotational distribution as well, consistent with the images collected at 320 nm. A 2D-REMPI spectrum was collected following dissociation of O3 near 330 nm, which showed evidence of contributions from O2 fragments in both the a1Δg and b1Σg+ states. The rotational distribution for the O2(b1Σg+, v = 0) product peaks at j = 32 and is narrower than that of the O2(a1Δg) fragment, consistent with distributions reported by O'Keeffe et al. at longer dissociation wavelengths ( J. Chem. Phys. 2002, 117, 8705-8709). At smaller radii in the 2D-REMPI spectrum, there is additional signal assigned to v = 1-4 of O2(b1Σg+), with rotational distributions similar to v = 0. The vibrational distribution of the O2(b1Σg+) fragment peaks at v = 0, with populations monotonically decreasing with increasing vibrational state. Ion image angular distributions of the O2(b1Σg+) fragment and the corresponding anisotropy parameters are also reported.

Ozone Photodissociation in the Singlet Channel at 226 nm

We report the rotational state distribution and vector correlations of the O₂(a ¹Δ_g, v = 0) fragments arising from the 226 nm photodissociation of jet-cooled O₃. Consistent with previously reported trends, the rotational distribution is shifted to higher rotational states with decreasing wavelength. We observe highly suppressed odd rotational state populations due to a strong Λ-doublet propensity. The measured rotational distribution is in agreement with classical trajectory calculations for the v = 0 products, although the distribution is slightly narrower than predicted. The spatial anisotropy follows the previously observed trend of decreasing β with increasing photon energy with β = 0.72 ± 0.14 for v = 0, j = 38. As expected for a triatomic molecule, the v-j correlation is consistent with v perpendicular to j, but the measured correlation is nonlimiting due, in part, to rotational and translational depolarization. The j-dependent line width of the O₂(a ¹Δ_g) REMPI spectrum is also discussed in connection with the lifetime of the resonant O₂(d ¹Π_g) state due to predissociation via the II ¹Π_g valence state.

Photodissociation Dynamics of CO-Forming Channels on the Ground-State Surface of Methyl Formate at 248 nm: Direct Dynamics Study and Assessment of Generalized Multicenter Impulsive Models

The photodissociation dynamics of methyl formate in the electronic ground state S₀, initiated by a 248 nm-wavelength laser, is studied by direct dynamics simulations. We analyze five channels, where four of them have as products CH₃OH + CO, one leading to the formation of three fragments, H₂CO + H₂ + CO, and a channel characterized by a roaming transition state. The analysis of energy distribution among the degrees of freedom of the product and the comparison with experimental results previously published by other groups provide the ingredients to distinguish the examined dissociation pathways. The interpretation of the results proves that the characterization of dissociation mechanisms must rely on a dynamics approach involving multiple electronic states, including considerations on the features of the S₁/S₀ conical intersection. Here, we also assess the generalized multicenter impulsive model, GMCIM, that has been designed for dissociation processes with exit barriers, and the energy distribution in the products is predicted on the basis of information from the saddle points and the intrinsic reaction coordinates. Main features, advantages, limits, and future perspectives of the method are reported and discussed.

Evidence for lambda doublet propensity in the UV photodissociation of ozone

The photodissociation of O₃ at 266 nm has been studied using velocity mapped ion imaging. We report temperature-dependent vector correlations for the O₂(a¹Δ_g, v = 0, j = 18-20) fragments at molecular beam temperatures of 70 K, 115 K, and 170 K. Both the fragment spatial anisotropy and the v-j correlations are found to be increasingly depolarized with increasing beam temperature. At all temperatures, the v-j correlations for the j = 19 state were shown to be reduced compared to those of j = 18 and 20, while no such odd/even rotational state difference was observed for the spatial anisotropy, consistent with previous measurements. We find that temperature-dependent differences in the populations and v-j correlations between the odd and even rotational states can be explained by a Λ-doublet propensity model. Although symmetry conservation should lead to formation of only the A' Λ-doublet component, and only even rotational states, out-of-plane rotation of the parent molecule breaks the planar symmetry and permits the formation of the A″ Λ-doublet component and odd rotational states. A simple classical model to treat the effect of parent rotation on the v-j correlation and the odd/even rotational population alternation reproduces both the current measurements and previously reported rotational distributions, suggesting that the "odd" behavior originates from a Λ-doublet propensity, and not from a mass independent curve crossing effect, as previously proposed.

Imaging the nonreactive collisional quenching dynamics of NO (A2Σ+) radicals with O2 (X3Σg -)

Nitric oxide (NO) radicals are ubiquitous chemical intermediates present in the atmosphere and in combustion processes, where laser-induced fluorescence is extensively used on the NO (A²Σ⁺ ← X²Π) band to report on fuel-burning properties. However, accurate fluorescence quantum yields and NO concentration measurements are impeded by electronic quenching of NO (A²Σ⁺) to NO (X²Π) with colliding atomic and molecular species. To improve predictive combustion models and develop a molecular-level understanding of NO (A²Σ⁺) quenching, we report the velocity map ion images and product state distributions of NO (X²Π, v″ = 0, J″, F_n, Λ) following nonreactive collisional quenching of NO (A²Σ⁺) with molecular oxygen, O₂ (X³Σ_g ^-). A novel dual-flow pulse valve nozzle is constructed and implemented to carry out the NO (A²Σ⁺) electronic quenching studies and to limit NO₂ formation. The isotropic ion images reveal that the NO-O₂ system evolves through a long-lived NO₃ collision complex prior to formation of products. Furthermore, the corresponding total kinetic energy release distributions support that O₂ collision coproducts are formed primarily in the c¹Σ_u ^- electronic state with NO (X²Π, v″ = 0, J″, F_n, Λ). The product state distributions also indicate that NO (X²Π) is generated with a propensity to occupy the Π(A″) Λ-doublet state, which is consistent with the NO π^* orbital aligned perpendicular to nuclear rotation. The deviations between experimental results and statistical phase space theory simulations illustrate the key role that the conical intersection plays in the quenching dynamics to funnel population to product rovibronic levels.

Invited Review Article: Photofragment imaging

Photodissociation studies in molecular beams that employ position-sensitive particle detection to map product recoil velocities emerged thirty years ago and continue to evolve with new laser and detector technologies. These powerful methods allow application of tunable laser detection of single product quantum states, simultaneous measurement of velocity and angular momentum polarization, measurement of joint product state distributions for the detected and undetected products, coincident detection of multiple product channels, and application to radicals and ions as well as closed-shell molecules. These studies have permitted deep investigation of photochemical dynamics for a broad range of systems, revealed new reaction mechanisms, and addressed problems of practical importance in atmospheric, combustion, and interstellar chemistry. This review presents an historical overview, a detailed technical account of the range of methods employed, and selected experimental highlights illustrating the capabilities of the method.

Nascent O2 (a 1Δg, v = 0, 1) rotational distributions from the photodissociation of jet-cooled O3 in the Hartley band

We report rotational distributions for the O₂ (a ¹Δ_g) fragment from the photodissociation of jet-cooled O₃ at 248, 266, and 282 nm. The rotational distributions show a population alternation that favors the even states, as previously reported for a 300 K sample by Valentini et al. [J. Chem. Phys. 86, 6745 (1987)]. The alternation from the jet-cooled precursor is much stronger than that observed by Valentini et al. and in contrast to their observations does not depend strongly on the O₂ (a ¹Δ_g) vibrational state or photolysis wavelength. The odd/even alternation diminishes substantially when the ozone beam temperature is increased from 60 to 200 K, confirming its dependence on parent internal energy. The magnitude of the even/odd alternation in product rotational states from the cold ozone sample, its temperature dependence, and other experimental and theoretical evidence reported since 1987 suggest that the alternation originates from a Λ-doublet propensity and not from a mass independent curve crossing effect, as previously proposed.

A generalized unimolecular impulsive model for curved reaction path

This work aims to introduce a generalized impulsive model for unimolecular dissociation processes. This model allows us to take into account the curvature of the reaction path explicitly. It is a generalization of the previously developed multi-center impulsive model [P.-Y. Tsai and K.-C. Lin, J. Phys. Chem. A 119, 29 (2015)]. Several limitations of conventional impulsive models are eliminated by this study: (1) Unlike conventional impulsive models, in which a single molecular geometry is responsible for the impulse determination, the gradients on the whole dissociation path are taken into account. The model can treat dissociation pathways with large curvatures and loose saddle points. (2) The method can describe the vibrational excitation of polyatomic fragments due to the bond formation by multi-center impulse. (3) The available energy in conventional impulsive models is separated into uncoupled statistical and impulsive energy reservoirs, while the interplay between these reservoirs is allowed in the new model. (4) The quantum state correlation between fragments can be preserved in analysis. Dissociations of several molecular systems including the roaming pathways of formaldehyde, nitrate radical, acetaldehyde, and glyoxal are chosen as benchmarks. The predicted photofragment energy and vector distributions are consistent with the experimental results reported previously. In these examples, the capability of the new model to treat the curved dissociation path, loose saddle points, polyatomic fragments, and multiple-body dissociation is verified. As a cheaper computational tool with respect to ab initio on-the-fly direct dynamic simulations, this model can provide detailed information on the energy disposal, quantum state correlation, and stereodynamics in unimolecular dissociation processes.

New experimental evidence to support roaming in the reaction Cl + isobutene (i-C4H8)

The reaction Cl + isobutene (i-C₄H₈) was reported by Suits et al. to proceed via, in addition to abstraction, an addition-elimination path following a roaming excursion of Cl; a near-zero translational energy release and an isotropic angular distribution observed at a small collision energy characterized this mechanism. We employed a new experimental method to further characterize this roaming mechanism through observation of the internal distribution of HCl (v, J) and their temporal behavior upon irradiation of a mixture of Cl₂C₂O₂ and i-C₄H₈ in He or Ar buffer gas. With 1–3 Torr buffer gas added to approach the condition of small collision energy, the intensities of emission of HCl (v = 1, 2) and the HCl production rates increased significantly; Ar shows a more significant effect than He because Ar quenches Cl more efficiently to reduce the collisional energy and facilitate the roaming path. According to kinetic modeling, the rate of addition-elimination (roaming) increased from k_E ≈ 2 × 10⁵ s⁻¹ when little buffer gas was present to ~1.9 × 10⁶ s⁻¹ when 2–3 Torr of Ar was added, and the branching ratio for formation of [HCl (v = 2)]/[HCl (v = 1)] increased from 0.02 ± 0.01 for abstraction to 0.06 ± 0.01 for roaming.

Competitive reaction pathways in vibrationally induced photodissociation of H2SO4

Vibrationally induced photodissociation of sulfuric acid into H2O + SO3 is investigated based on reactive molecular dynamics (MD) simulations. Multisurface adiabatic reactive MD simulations allow us to follow both, H-transfer and water elimination after excitation of the ν9 OH-stretching mode. Analysis of several thousand trajectories finds that the H2O and SO3 fragments have distinct final state distributions with respect to translational, rotational, and vibrational degrees of freedom. Rotational distributions peak at quantum numbers j ≤ 5 for water and j ≈ 60 for SO3. The final state distributions should be useful in identifying products in forthcoming experiments. Based on the MD trajectories, a kinetic scheme has been developed which is able to explain most of the trajectory data and suggests that IVR is very rapid. Typical lifetimes of the excited complex range from several 10 picoseconds to hundreds of nanoseconds, depending on the excitation level. Including temperature and pressure profiles characteristic for the stratosphere in the kinetic model shows that excitations higher than ν9 = 4 can significantly contribute to the photolysis rate. This extends and specifies earlier work in that multi-level modeling is required to understand the significance of vibrationally induced decomposition pathways of sulfuric acid in the middle atmosphere.

Two roaming pathways in the photolysis of CH3CHO between 328 and 308 nm

We attribute the two product-state distributions previously seen in CH₃CHO photodissociation to CH₃-roaming and H-roaming, unifying all previous experimental results.

Quasiclassical trajectory studies of 18O(3P) + NO2 isotope exchange and reaction to O2 + NO on D0 and D1 potentials

We report quasiclassical trajectory calculations for the bimolecular reaction (18)O((3)P) + NO2 on the recent potential energy surfaces of the ground (D0) and first excited (D1) states of NO3 [B. Fu, J. M. Bowman, H. Xiao, S. Maeda, and K. Morokuma, J. Chem. Theory. Comput. 9, 893 (2013)]. The branching ratio of isotope exchange versus O2 + NO formation, as well as the product angular distributions and energy and rovibrational state distributions are presented. The calculations are done at the collision energy of relevance to recent crossed beam experiments [K. A. Mar, A. L. Van Wyngarden, C.-W. Liang, Y. T. Lee, J. J. Lin, and K. A. Boering, J. Chem. Phys. 137, 044302 (2012)]. Very good agreement is achieved between the current calculations and these experiments for the branching ratio and final translational energy and angular distributions of isotope exchange products (16)O((3)P) + NO2 and O2 + NO formation products. The reactant (18)O atom results in (18)O(16)O but not N(18)O for the O2 + NO formation product channel, consistent with the experiment. In addition, the detailed vibrational and rotational state information of diatomic molecules calculated currently for the (34)O2 + NO formation channel on D0 and D1 states are in qualitative agreement with the previous experimental and theoretical results of the photodissociation of NO3 and are consistent with older thermal bimolecular kinetics measurements.

Exploring Multiple Potential Energy Surfaces: Photochemistry of Small Carbonyl Compounds

In theoretical studies of chemical reactions involving multiple potential energy surfaces (PESs) such as photochemical reactions, seams of intersection among the PESs often complicate the analysis. In this paper, we review our recipe for exploring multiple PESs by using an automated reaction path search method which has previously been applied to single PESs. Although any such methods for single PESs can be employed in the recipe, the global reaction route mapping (GRRM) method was employed in this study. By combining GRRM with the proposed recipe, all critical regions, that is, transition states, conical intersections, intersection seams, and local minima, associated with multiple PESs, can be explored automatically. As illustrative examples, applications to photochemistry of formaldehyde and acetone are described. In these examples as well as in recent applications to other systems, the present approach led to discovery of many unexpected nonadiabatic pathways, by which some complicated experimental data have been explained very clearly.

Automated Exploration of Photolytic Channels of HCOOH: Conformational Memory via Excited-State Roaming

To elucidate the photodissociation mechanism of HCOOH, we systematically explored reaction pathways starting from the first excited singlet state (S1) by using automated reaction path search methods. All critical points, that is, minima, transition states, minimum energy conical intersections, and minima on seam of crossing, for the S0, T1, and S1 potential energy surfaces (PESs) obtained in the present search were optimized at the CASPT2 level. The structure list obtained by the search explained all experimentally reported photolytic channels. A new mechanism for the previously suggested but unexplained conformational memory in the 193 nm photolysis is proposed, which involves two steps: partial dissociation and succeeding roaming of one of H atoms on the S1 PES, followed by intramolecular recombination on the S0 PES after radiationless transition through a conical intersection at a partially dissociated geometry. This is partially similar to the excited-state roaming recently discovered for the NO3 radical.