The triplet channel in the photodissociation of ozone in the Hartley band: Classical trajectory surface hopping analysis

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.

Origin of the "odd" behavior in the ultraviolet photochemistry of ozone

The origin of the even-odd rotational state population alternation in the ¹⁶O₂(a ¹Δ_g) fragments resulting from the ultraviolet (UV) photodissociation of ¹⁶O₃, a phenomenon first observed over 30 years ago, has been elucidated using full quantum theory. The calculated ¹⁶O₂(a ¹Δ_g) rotational state distribution following the 266-nm photolysis of 60 K ozone shows a strong even-odd propensity, in excellent agreement with the new experimental rotational state distribution measured under the same conditions. Theory indicates that the even rotational states are significantly more populated than the adjacent odd rotational states because of a preference for the formation of the A' Λ-doublet, which can only occupy even rotational states due to the exchange symmetry of the two bosonic ¹⁶O nuclei, and thus not as a result of parity-selective curve crossing as previously proposed. For nonrotating ozone, its dissociation on the excited B¹A' state dictates that only A' Λ-doublets are populated, due to symmetry conservation. This selection rule is relaxed for rotating parent molecules, but a preference still persists for A' Λ-doublets. The A''/A' ratio increases with increasing ozone rotational quantum number, and thus with increasing temperature, explaining the previously observed temperature dependence of the even-odd population alternation. In light of these results, it is concluded that the previously proposed parity-selective curve-crossing mechanism cannot be a source of heavy isotopic enrichment in the atmosphere.

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.

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.

Rovibrational energy transfer and dissociation in O2-O collisions

A set of state-specific transition rates for each rovibrational level is generated for the O2(X(3)Σ(g)(-))-O(3)P system using the quasi-classical trajectory method at temperatures observed in hypersonic flows. A system of master equations describes the relaxation of the rovibrational ensemble to thermal equilibrium under ideal heat bath conditions at a constant translational temperature. Vibrational and rotational relaxation times, obtained from the average internal energies, exhibit a pattern inherent in a chemically reactive collisional pair. An intrinsic feature of the O3 molecular system with a large attractive potential is a weak temperature dependence of the rovibrational transition rates. For this reason, the quasi-steady vibrational and rotational temperatures experience a maximum at increasing translational temperature. The energy rate coefficients, that characterize the average loss of internal energy due to dissociation, quickly diminish at high temperatures, compared to other molecular systems.

Intermediate photofragment distributions as probes of non-adiabatic dynamics at conical intersections: application to the Hartley band of ozone

Quantum dynamics at a reactive two-state conical intersection lying outside the Franck–Condon zone is studied for a prototypical reaction of ultraviolet photodissociation of ozone in the Hartley band.

Signatures of a conical intersection in photofragment distributions and absorption spectra: photodissociation in the Hartley band of ozone

Photodissociation of ozone in the near UV is studied quantum mechanically in two excited electronic states coupled at a conical intersection located outside the Franck-Condon zone. The calculations, performed using recent ab initio PESs, provide an accurate description of the photodissociation dynamics across the Hartley/Huggins absorption bands. The observed photofragment distributions are reproduced in the two electronic dissociation channels. The room temperature absorption spectrum, constructed as a Boltzmann average of many absorption spectra of rotationally excited parent ozone, agrees with experiment in terms of widths and intensities of diffuse structures. The exit channel conical intersection contributes to the coherent broadening of the absorption spectrum and directly affects the product vibrational and translational distributions. The photon energy dependences of these distributions are strikingly different for fragments created along the adiabatic and the diabatic paths through the intersection. They can be used to reverse engineer the most probable geometry of the non-adiabatic transition. The angular distributions, quantified in terms of the anisotropy parameter β, are substantially different in the two channels due to a strong anticorrelation between β and the rotational angular momentum of the fragment O2.

Ab initio quantum mechanical study of the O((1)D) formation in the photolysis of ozone between 300 and 330 nm

Spin-allowed production of O((1)D) in the near-UV photolysis of ozone is studied using ab initio potential energy surfaces and quantum mechanics. The O((1)D) quantum yield, reconstructed from the absolute cross sections for eight initial vibrational states in the ground electronic state, is shown to agree with the measurements in a broad range of photolysis wavelengths and temperatures. Relative contributions of one- and two-quantum stretching and bending initial excitations are quantified, with the contribution of the antisymmetric stretch being dominant for lambda < 330 nm. Large scale structures in the low-resolution quantum yield are shown to reflect excitations in the high-frequency short bond stretch in the upper electronic state. Spin-forbidden contribution to the O((1)D) quantum yield at wavelengths lambda > 320 nm is estimated using ab initio energies of the triplet states and their spin-orbit couplings.

New theoretical investigations of the photodissociation of ozone in the Hartley, Huggins, Chappuis, and Wulf bands

We review recent theoretical studies of the photodissociation of ozone in the wavelength region from 200 nm to 1100 nm comprising four major absorption bands: Hartley and Huggins (near ultraviolet), Chappuis (visible), and Wulf (near infrared). The quantum mechanical dynamics calculations use global potential energy surfaces obtained from new high-level electronic structure calculations. Altogether nine electronic states are taken into account in the theoretical descriptions: four 1A', two 1A'', one 3A' and two 3A'' states. Of particular interest is the analysis of diffuse vibrational structures, which are prominent in all absorption bands, and their dynamical origin and assignment. Another focus is the effect of non-adiabatic coupling on lifetimes in the excited states and on the population of the specific electronic product channels.

The photodissociation dynamics of ozone at 193nm: An O(D21) angular momentum polarization study

Polarized laser photolysis, coupled with resonantly enhanced multiphoton ionization detection of O(1D2) and velocity-map ion imaging, has been used to investigate the photodissociation dynamics of ozone at 193 nm. The use of multiple pump and probe laser polarization geometries and probe transitions has enabled a comprehensive characterization of the angular momentum polarization of the O(1D2) photofragments, in addition to providing high-resolution information about their speed and angular distributions. Images obtained at the probe laser wavelength of around 205 nm indicate dissociation primarily via the Hartley band, involving absorption to, and diabatic dissociation on, the B 1B2(3 1A1) potential energy surface. Rather different O(1D2) speed and electronic angular momentum spatial distributions are observed at 193 nm, suggesting that the dominant excitation at these photon energies is to a state of different symmetry from that giving rise to the Hartley band and also indicating the participation of at least one other state in the dissociation process. Evidence for a contribution from absorption into the tail of the Hartley band at 193 nm is also presented. A particularly surprising result is the observation of nonzero, albeit small values for all three rank K = 1 orientation moments of the angular momentum distribution. The polarization results obtained at 193 and 205 nm, together with those observed previously at longer wavelengths, are interpreted using an analysis of the long range quadrupole-quadrupole interaction between the O(1D2) and O2(1Deltag) species.

Semiclassical nonadiabatic dynamics based on quantum trajectories for the O(P3,D1)+H2 system

The O(3P,1D) + H2 --> OH + H reaction is studied using trajectory dynamics within the approximate quantum potential approach. Calculations of the wave-packet reaction probabilities are performed for four coupled electronic states for total angular momentum J = 0 using a mixed coordinate/polar representation of the wave function. Semiclassical dynamics is based on a single set of trajectories evolving on an effective potential-energy surface and in the presence of the approximate quantum potential. Population functions associated with each trajectory are computed for each electronic state. The effective surface is a linear combination of the electronic states with the contributions of individual components defined by their time-dependent average populations. The wave-packet reaction probabilities are in good agreement with the quantum-mechanical results. Intersystem crossing is found to have negligible effect on reaction probabilities summed over final electronic states.

Absorption spectrum and assignment of the Chappuis band of ozone

New global diabatic potential energy surfaces of the electronic states 1B1 and 1A2 of ozone and the non-adiabatic coupling surface between them are constructed from electronic structure calculations. These surfaces are used to study the visible photodissociation in the Chappuis band by means of quantum mechanical calculations. The calculated absorption spectrum and its absolute intensity are in good agreement with the experimental results. A vibrational assignment of the diffuse structures in the Chappuis band system is proposed on the basis of the nodal structures of the underlying resonance states.

Infrared spectrum of cyclic ozone: A theoretical investigation

The infrared absorption spectrum of cyclic ozone is calculated by means of a new ab initio potential energy surface, the dipole moment function, and exact quantum mechanical dynamics calculations. Five different isotopomers are considered. The absorption line for excitation of the bending fundamental near 800 cm(-1) is by far the strongest band; all other bands are more than one order of magnitude less intense. This spectral pattern as well as the isotope shifts for the various isotopomers are important for identifying cyclic ozone. Several possibilities for accessing the ring minimum of cyclic ozone are also discussed on the basis of recent electronic structure calculations.

The photodissociation of ozone in the Hartley band: A theoretical analysis

Three-dimensional diabatic potential energy surfaces for the lowest four electronic states of ozone with 1A' symmetry-termed X, A, B, and R-are constructed from electronic structure calculations. The diabatization is performed by reassigning corresponding energy points. Although approximate, these diabatic potential energy surfaces allow one to study the uv photodissociation of ozone on a level of theory not possible before. In the present work photoexcitation in the Hartley band and subsequent dissociation into the singlet channel, O3X+hnu-->O(1D)+O2(a 1Deltag), are investigated by means of quantum mechanical and classical trajectory calculations using the diabatic potential energy surface of the B state. The calculated low-resolution absorption spectrum as well as the vibrational and rotational state distributions of O2(a 1Deltag) are in good agreement with available experimental results.