Photodissociation of ozone in the Hartley band: Potential energy surfaces, nonadiabatic couplings, and singlet/triplet branching ratio

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.

Wave Packet Calculation of Absolute UV Cross Section of Criegee Intermediates

Criegee intermediates, R₁R₂COO, are reactive species formed in the atmosphere through the ozonolysis of alkenes. They have an intense ultraviolet (UV) adsorption between 300 to 400 nm. However, experimentally determining the absolute cross sections is not easy. We used wave packet propagation on an one-dimensional adiabatic potential energy curve (PEC) along the OO bond to simulate the UV spectra for various Criegee intermediates. Our results showed a very fast, ∼20 fs, decay out of the Franck-Condon region. This gives justification for using the semiclassical approach which was utilized in previous studies. From the comparison of various quantum chemistry methods, we found that multireference methods can give spectra with a width and cross section reproducing the experimental results, while single reference methods tend to give narrower skewed peaks with a larger cross section. From the test using wave packet propagation on various approximated PECs and transition moment functions, we show that the Gaussian approximation within the reflection method is valid. In addition, we found that we can obtain peak positions that reproduce the experimental results by shifting those obtained by MRCI+Q, CASSCF, EOMCCSD, and TDCAMB3LYP by -0.2, -1.0, -0.3, and -0.5 eV, respectively. The Gaussian approximation using peak position, oscillator strength, and peak width from MRCI+Q is a cost-effective way to simulate the UV spectra of Crigee intermediates for which experimental determination may be hard.

Reactive, Inelastic, and Dissociation Processes in Collisions of Atomic Nitrogen with Molecular Oxygen

Collisions of atomic nitrogen with molecular oxygen have been treated with the quasiclassical trajectory method (QCT) in order to obtain a complete database of vibrationally detailed cross sections and rate coefficients for reactive, inelastic, and dissociation processes. For reaction rate coefficients, the agreement with experimental and theoretical data in the literature is excellent on the whole available interval 300-5000 K, with reliable extension to 20,000 K. For the inelastic case and for dissociation, no comparisons are available; therefore, a study of QCT reliability is proposed. In the inelastic case, it is found that "purely inelastic" and "quasireactive" collisions show not only different mechanisms but also different QCT levels of reliability at low energy. For dissociation, similar considerations bring to the conclusion that for the present collisional system, the QCT method is appropriate on the whole energy range studied. Rate coefficients for all the processes studied are provided in the electronic form.

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.

Photodissociation by Circularly Polarized Light Yields Photofragment Alignment in Ozone Arising Solely from Vibronic Interactions

We present a direct determination of photofragment alignment produced by circularly polarized light in photolysis of a planar polyatomic molecule. This alignment arises via a new mechanism involving coherent excitation of two mutually perpendicular in-plane transition dipole moment components. The alignment is described by a new anisotropy parameter γ_{2}^{'} that was isolated by a unique laser polarization geometry. The determination of the parameter γ_{2}^{'} was realized in ozone photolysis at 266 nm where dc slice images of O(^{1}D_{2}) atomic fragments were acquired. A model developed for interpretation of the photolysis mechanism shows that it can exist only in case of failure of the Born-Oppenheimer approximation when electronic and vibrational (vibronic) interactions have to be taken into account. This finding suggests that determination of the alignment parameter γ_{2}^{'} can be used as a key for direct insight into vibronic interactions in photolysis of polyatomic molecules. The results obtained for ozone photolysis via the Hartley band showed significant γ_{2}^{'} alignment but little recoil speed dependence, consistent with the notion that, as opposed to the situation for derivative coupling, under our experimental conditions, the vibronic contributions to the nonadiabatic dynamics are not dependent on recoil speed.

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.

Theory and Calculation of the Phosphorescence Phenomenon

Phosphorescence is a phenomenon of delayed luminescence that corresponds to the radiative decay of the molecular triplet state. As a general property of molecules, phosphorescence represents a cornerstone problem of chemical physics due to the spin prohibition of the underlying triplet-singlet emission and because its analysis embraces a deep knowledge of electronic molecular structure. Phosphorescence is the simplest physical process which provides an example of spin-forbidden transformation with a characteristic spin selectivity and magnetic field dependence, being the model also for more complicated chemical reactions and for spin catalysis applications. The bridging of the spin prohibition in phosphorescence is commonly analyzed by perturbation theory, which considers the intensity borrowing from spin-allowed electronic transitions. In this review, we highlight the basic theoretical principles and computational aspects for the estimation of various phosphorescence parameters, like intensity, radiative rate constant, lifetime, polarization, zero-field splitting, and spin sublevel population. Qualitative aspects of the phosphorescence phenomenon are discussed in terms of concepts like structure-activity relationships, donor-acceptor interactions, vibronic activity, and the role of spin-orbit coupling under charge-transfer perturbations. We illustrate the theory and principles of computational phosphorescence by highlighting studies of classical examples like molecular nitrogen and oxygen, benzene, naphthalene and their azaderivatives, porphyrins, as well as by reviewing current research on systems like electrophosphorescent transition metal complexes, nucleobases, and amino acids. We furthermore discuss modern studies of phosphorescence that cover topics of applied relevance, like the design of novel photofunctional materials for organic light-emitting diodes (OLEDs), photovoltaic cells, chemical sensors, and bioimaging.

Attosecond electronic and nuclear quantum photodynamics of ozone monitored with time and angle resolved photoelectron spectra

Recently we reported a series of numerical simulations proving that it is possible in principle to create an electronic wave packet and subsequent electronic motion in a neutral molecule photoexcited by a UV pump pulse within a few femtoseconds. We considered the ozone molecule: for this system the electronic wave packet leads to a dissociation process. In the present work, we investigate more specifically the time-resolved photoelectron angular distribution of the ozone molecule that provides a much more detailed description of the evolution of the electronic wave packet. We thus show that this experimental technique should be able to give access to observing in real time the creation of an electronic wave packet in a neutral molecule and its impact on a chemical process.

Photodissociation dynamics of OCS near 214 nm using ion imaging

The OCS photodissociation dynamics of the dominant S((1)D2) channel near 214 nm have been studied using velocity map ion imaging. We report a CO vibrational branching ratio of 0.79:0.21 for v = 0:v = 1, indicating substantially higher vibrational excitation than that observed at slightly longer wavelengths. The CO rotational distribution is bimodal for both v = 0 and v = 1, although the bimodality is less pronounced than at longer wavelengths. Vector correlations, including rotational alignment, indicate that absorption to both the 2(1)A' (A) and 1(1)A″ (B) states is important in the lower-j part of the rotational distribution, while only 2(1)A' state absorption contributes to the upper part; this conclusion is consistent with work at longer wavelengths. Classical trajectory calculations including surface hopping reproduce the measured CO rotational distributions and their dependence on wavelength well, though they underestimate the v = 1 population. The calculations indicate that the higher-j peak in the rotational distribution arises from molecules that begin on the 2(1)A' state but make nonadiabatic transitions to the 1(1)A' (X) state during the dissociation, while the lower-j peak arises from direct photodissociation on either the 2(1)A' or the 1(1)A″ states, as found in previous work.

The transition from the open minimum to the ring minimum on the ground state and on the lowest excited state of like symmetry in ozone: A configuration interaction study

The metastable ring structure of the ozone 1(1)A1 ground state, which theoretical calculations have shown to exist, has so far eluded experimental detection. An accurate prediction for the energy difference between this isomer and the lower open structure is therefore of interest, as is a prediction for the isomerization barrier between them, which results from interactions between the lowest two (1)A1 states. In the present work, valence correlated energies of the 1(1)A1 state and the 2(1)A1 state were calculated at the 1(1)A1 open minimum, the 1(1)A1 ring minimum, the transition state between these two minima, the minimum of the 2(1)A1 state, and the conical intersection between the two states. The geometries were determined at the full-valence multi-configuration self-consistent-field level. Configuration interaction (CI) expansions up to quadruple excitations were calculated with triple-zeta atomic basis sets. The CI expansions based on eight different reference configuration spaces were explored. To obtain some of the quadruple excitation energies, the method of Correlation Energy Extrapolation by Intrinsic Scaling was generalized to the simultaneous extrapolation for two states. This extrapolation method was shown to be very accurate. On the other hand, none of the CI expansions were found to have converged to millihartree (mh) accuracy at the quadruple excitation level. The data suggest that convergence to mh accuracy is probably attained at the sextuple excitation level. On the 1(1)A1 state, the present calculations yield the estimates of (ring minimum-open minimum) ∼45-50 mh and (transition state-open minimum) ∼85-90 mh. For the (2(1)A1-(1)A1) excitation energy, the estimate of ∼130-170 mh is found at the open minimum and 270-310 mh at the ring minimum. At the transition state, the difference (2(1)A1-(1)A1) is found to be between 1 and 10 mh. The geometry of the transition state on the 1(1)A1 surface and that of the minimum on the 2(1)A1 surface nearly coincide. More accurate predictions of the energy differences also require CI expansions to at least sextuple excitations with respect to the valence space. For every wave function considered, the omission of the correlations of the 2s oxygen orbitals, which is a widely used approximation, was found to cause errors of about ±10 mh with respect to the energy differences.

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.

Investigation of 3-fragment photodissociation of O3 at 193.4 and 157.6 nm by coincident measurements

Photodissociation of the ozone molecule at 193.4 nm (6.41 eV) and 157.6 nm (7.87 eV) is studied by fast-beam translational spectroscopy. Coincident detection of the dissociation products allows direct observation of the 3-fragment channel and determination of its kinematic parameters. The results indicate that at each wavelength, 3-fragment dissociation proceeds through synchronous concerted bond breaking, but the energy partitioning among the fragments is different. The branching fraction of the 3-fragment channel increases from 5.2(6)% at 193.4 nm to 26(4)% at 157.6 nm, in agreement with previous studies. It is shown that vibrational excitation of the symmetric stretch mode in O3 molecules created by photodetachment of O(3)(-) anion enhances the absorption efficiency, especially at 193.4 nm, but does not have a strong effect on the 3-fragment dissociation.

On the bonding nature of ozone (O3) and its sulfur-substituted analogues SO2, OS2, and S3: correlation between their biradical character and molecular properties

We investigate the bonding mechanism in ozone (O3) and its sulfur-substituted analogues, SO2, OS2, and S3. By analyzing their ground-state multireference configuration interaction wave functions, we demonstrate that the bonding in these systems can be represented as a mixture of a closed-shell structure with one and a half bonds between the central and terminal atoms and an open-shell structure with a single bond and two lone electrons on each terminal atom (biradical). The biradical character (β) further emerges as a simple measure of the relative contribution of those two classical Lewis structures emanating from the interpretation of the respective wave functions. Our analysis yields a biradical character of 3.5% for OSO, 4.4% for SSO, 11% for S3, 18% for O3, 26% for SOO, and 35% for SOS. The size/electronegativity of the end atoms relative to the central one is the prevalent factor for determining the magnitude of β: smaller and more electronegative central atoms better accommodate a pair of electrons facilitating the localization of the remaining two lone π-electrons on each of the end atoms, therefore increasing the weight of the second picture in the mixed bonding scenario (larger β). The proposed mixture of these two bonding scenarios allows for the definition of the bond order of the covalent bonds being (3-β)/2, and this accounts for the different O-O, S-S, or S-O bond lengths in the triatomic series. The biradical character was furthermore found to be a useful concept for explaining several structural and energetic trends in the series: larger values of β mark a smaller singlet-triplet splitting, closer bond lengths in the ground (1)A' and the first excited (3)A' states, and larger bond dissociation and atomization energies in the ground state. The latter explains the relative energy difference between the OSS/SOS and OOS/OSO isomers due to their different β values.

Photodissociation of ozone from 321 to 329 nm: the relative yields of O(3P2) with O2(X 3Σg(-)), O2(a 1Δg) and O2(b 1Σg(+))

Product imaging of O((3)P2) following dissociation of ozone has been used to determine the relative yields of the product channels O((3)P2) + O2(X (3)Σg(-)) of ozone. All three channels are prominent at all wavelengths investigated. O2 vibrational distributions for each channel and each wavelength are also estimated assuming Boltzmann rotational distributions. Averaged over wavelength in the measured range, the yields of the O((3)P2) + O2(X (3)Σg(-)), O((3)P2) + O2(a (1)Δg), and O((3)P2) + O2(b (1)Σg(+)) channels are 0.36, 0.31,and 0.34, respectively. Photofragment distributions in the spin-allowed channel O((3)P) + O2(X (3)Σg(-)) are compared with the results of quantum mechanical calculations on the vibronically coupled PESs of the singlet states B (optically bright) and R (repulsive). The experiments suggest that considerably more vibrational excitation and less rotational excitation occur than predicted by the quantum calculations. The rotational distributions, adjusted to fit the experimental images, suggest that the dissociation takes place from a more linear configuration than the Franck-Condon bending angle of 117°. The dissociation at most wavelengths results in a positive value of the anisotropy parameter, β, both in the experiment and in the calculations. Calculations indicate that both nonadiabatic transitions and intersystem crossings substantially reduce β below the nominal value of 2.

Comparison of the Huggins band for six ozone isotopologues: vibrational levels and absorption cross section

By use of the 3(1)A' ab initio potential energy surface (PES) of ozone and the multi-configuration time-dependent Hartree program for wavepacket propagation, we have determined numerous eigenstates of this state for six ozone isotopologues. These bound vibrational levels are the upper levels of the Huggins band, which covers the range from 27,000 to ~33,000 cm(-1). This study extends our previous work on the Hartley band, which was limited to the range ~32,000-50,000 cm(-1). Four isotopologues, (16)O(3), (16)O(17)O(16)O, (16)O(18)O(16)O, and (18)O(3) (noted hereafter 666, 676, 686, and 888), are symmetric, and two are asymmetric, (17)O(16)O(2) and (18)O(16)O(2) (noted hereafter 667 and 668). The PES of the 3(1)A' state has two equivalent minima of C(s) symmetry located at ~27,000 cm(-1) above the X(1)A(1) ground state. The equilibrium geometry of these two minima is r(e(1)) = 2.28 a(0), r(e(2)) = 3.2 a(0), and θ(e) = 107°. The dissociation limit of this PES, which correlates to the O((1)D) + O(2) ((1)Δ) "singlet" channel, is about 4300 cm(-1) above the two minima. For the (16)O(3) isotopologue, the 120 lowest bound eigenstates have been calculated and partially assigned up to 800 cm(-1) below the dissociation limit. The 60 lower eigenstates are easily assignable in term of three normal modes, the "long" bond (ν(1)), the bending (ν(2)), and the "short" bond (ν(3)). A new family of wave functions, aligned along the dissociation channels, appears at 3782 cm(-1) above the 3(1)A' (0,0,0) level. The 3(1)A' vibrational levels and the corresponding intensity factors from the (000), (010), (100), and (001) levels of the X(1)A(1) ground state have been calculated for the six isotopologues. The Huggins absorption cross sections of the six isotopologues have been calculated from the 3(1)A' vibrational energy levels and the corresponding intensity factors. The rotational envelope of each vibronic band has been empirically described by an ad hoc function. The ratio of the Huggins cross section of each ozone isotopologue with one of (16)O(3) provides the fractionation factor of each ozone isotopologue as a function of the photon energy. These various fractionation factors will allow predicting enrichments due to photolysis by various light sources like the actinic flux.

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.