The ultra-violet photodissociation of ozone revisited

Valence Bond Alternative Yielding Compact and Accurate Wave Functions for Challenging Excited States. Application to Ozone and Sulfur Dioxide

A novel state-averaged version of ab initio nonorthogonal valence bond method is described, for the sake of accurate theoretical studies of excited states in the valence bond framework. With respect to standard calculations in the molecular orbital framework, the state-averaged breathing-orbital valence bond (BOVB) method has the advantage to be free from the penalizing constraint for the ground and excited state(s) to share the same unique set of orbitals. The ability of the BOVB method to faithfully describe excited states and to compute accurate transition energies from the ground state is tested on the five lowest-lying singlet electronic states of ozone and sulfur dioxide, among which 1¹B₂ and 2¹A₁ are the challenging ones. As the 1¹A₂, 1¹B₁, and 1¹B₂ states are of different symmetries than the ground state, they can be calculated at the state-specific BOVB level. On the other hand, the 2¹A₁ states and the 1¹A₁ ground states, which are of like symmetry, are calculated with the state-averaged BOVB technique. In all cases, the calculated vertical energies are close to the experimental values when available, and at par with the most sophisticated calculations in the molecular framework, despite the extreme compactness of the BOVB wave functions, made of no more than 5-9 valence bond structures in all cases. The features that allow the combination of compactness and accuracy in challenging cases are analyzed. For the "ionic" 1¹B₂ states, which are the site of important charge fluctuations, it is because of the built-in dynamic correlation inherent to the BOVB method. For the 2¹A₁ ones, this is the fact that these states have the degree of freedom of having different orbitals than the ground states, even though they are of like symmetry and calculated simultaneously using the newly implemented state-average BOVB algorithm. Finally, the description of the excited states in terms of Lewis structures is insightful, rationalizing the fast ring closure for the 2¹A₁ state of ozone and predicting some diradical character in the so-called "ionic" 1¹B₂ states.

Intensification of depolymerization of polyacrylic acid solution using different approaches based on ultrasound and solar irradiation with intensification studies

Depolymerization of polyacrylic acid (PAA) as sodium salt has been investigated using ultrasonic and solar irradiations with process intensification studies based on combination with hydrogen peroxide (H2O2) and ozone (O3). Effect of solar intensity, ozone flow and ultrasonic power dissipation on the extent of viscosity reduction has been investigated for individual treatment approaches. The combined approaches such as US+solar, solar+O3, solar+H2O2, US+H2O2 and US+O3 have been subsequently investigated under optimum conditions and established to be more efficient as compared to individual approaches. Approach based on US (60W)+solar+H2O2 (0.01%) resulted in the maximum extent of viscosity reduction as 98.97% in 35min whereas operation of solar+H2O2 (0.01%), US (60W), H2O2 (0.3%) and solar irradiation resulted in about 98.08%, 90.13%, 8.91% and 90.77% intrinsic viscosity reduction in 60min respectively. Approach of US (60W)+solar+ozone (400mg/h flow rate) resulted in extent of viscosity reduction as 99.47% in 35min whereas only ozone (400mg/h flow rate), ozone (400mg/h flow rate)+US (60W) and ozone (400mg/h flow rate)+solar resulted in 69.04%, 98.97% and 98.51% reduction in 60min, 55min and 55min respectively. The chemical identity of the treated polymer using combined approaches was also characterized using FTIR (Fourier transform infrared) spectra and it was established that no significant structural changes were obtained during the treatment. Overall, it can be said that the combination technique based on US and solar irradiations in the presence of hydrogen peroxide is the best approach for the depolymerization of PAA solution.

Intensified depolymerization of aqueous polyacrylamide solution using combined processes based on hydrodynamic cavitation, ozone, ultraviolet light and hydrogen peroxide

The present work deals with intensification of depolymerization of polyacrylamide (PAM) solution using hydrodynamic cavitation (HC) reactors based on a combination with hydrogen peroxide (H2O2), ozone (O3) and ultraviolet (UV) irradiation. Effect of inlet pressure in hydrodynamic cavitation reactor and power dissipation in the case of UV irradiation on the extent of viscosity reduction has been investigated. The combined approaches such as HC+UV, HC+O3, HC+H2O2, UV+H2O2 and UV+O3 have been subsequently investigated and found to be more efficient as compared to individual approaches. For the approach based on HC+UV+H2O2, the extent of viscosity reduction under the optimized conditions of HC (3 bar inlet pressure)+UV (8 W power)+H2O2 (0.2% loading) was 97.27% in 180 min whereas individual operations of HC (3 bar inlet pressure) and UV (8 W power) resulted in about 35.38% and 40.83% intrinsic viscosity reduction in 180 min respectively. In the case of HC (3 bar inlet pressure)+UV (8 W power)+ozone (400 mg/h flow rate) approach, the extent of viscosity reduction was 89.06% whereas individual processes of only ozone (400 mg/h flow rate), ozone (400 mg/h flow rate)+HC (3 bar inlet pressure) and ozone (400 mg/h flow rate)+UV (8 W power) resulted in lower extent of viscosity reduction as 50.34%, 60.65% and 75.31% respectively. The chemical structure of the treated PAM by all approaches was also characterized using FTIR (Fourier transform infrared) spectra and it was established that no significant chemical structure changes were obtained during the treatment. Overall, it can be said that the combination of HC+UV+H2O2 is an efficient approach for the depolymerization of PAM solution.

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.

PHOTODISSOCIATION OF OZONE IN THE HARTLEY BAND: FRAGMENT ROTATIONAL QUANTUM STATE DISTRIBUTIONS

In this paper, we have calculated the rotational state distributions following the photodissociation of ozone in the Hartley band with total angular momentum J'=1. The calculated results are obtained by using time-dependent wave packet calculations on the Sheppard–Walker potential energy surface (PES). It is found that the physically more correct treatment with J'=1 semi-quantitatively reproduces the rotational state distributions of the CARS. Compared with the previous theoretical works, which had taken J=0 on both the ground and excited potential surface, J'=1 treatment makes the rotational distributions of the fragment closer to the experimental ones.

The photodissociation of ozone: A quasi-classical approach to a quantum dynamics problem

It is extremely difficult to execute full quantum dynamics calculations on complex systems (more then three degrees of freedom) due to the exponential increase in computer resources required by these methods. Classical mechanics simulations do not suffer from this problem, but are unable to treat quantum mechanical phenomena such as non-adiabatic effects, which often play a vital role in photochemical processes. A method has been implemented for carrying out dynamical calculations using quasi-classical theory. The time dependent Schrödinger equation is solved using a swarm of trajectories treated under Newtonian laws and Tully's fewest switches trajectory surface hopping is applied to implement the surface switches. The method was applied to ozone, looking at the photodissociation that takes place after excitation into the Chappuis band of the absorption spectrum. While the goal is to treat larger systems, comparison can be made for ozone with numerically exact wavepacket calculations. The method proved successful at calculating quantities such as the rate of population transfer, but there were discrepancies in the details, especially when surface switching occurred from the lower state.

Speed dependent rotational angular momentum polarization of the O2 (aΔg1) fragment following ozone photolysis in the wavelength range 248–265nm

The translational anisotropy and rotational angular momentum polarization of a selection of rotational states of the O2 (a 1Deltag; v=0) photofragment formed from ozone photolysis at 248, 260, and 265 nm have been determined using the technique of resonance enhanced multiphoton ionization in combination with time of flight mass spectrometry. At 248 nm, the dissociation is well described as impulsive in nature with all rotational states exhibiting similarly large, near-limiting values for the bipolar moments describing their angular momentum alignment and orientation. At 265 nm, however, the angular momentum polarization parameters determined for consecutive odd and even rotational states exhibit clear differences. Studies at the intermediate wavelength of 260 nm strongly suggest that such a difference in the angular momentum polarization is speed dependent and this proposal is consistent with the angular momentum polarization parameters extracted and reported previously for longer photolysis wavelengths [G. Hancock et al., Phys. Chem. Chem. Phys. 5, 5386 (2003); S. J. Horrocks et al., J. Chem. Phys. 126, 044308 (2007)]. The alternation of angular momentum polarization for successive odd and even J states may be a consequence of the different mechanisms leading to the formation of the two O2 (a 1Deltag) Lambda doublets. Specifically, the involvement of out of plane parent rotational motion is proposed as the origin for the observed depolarization for the Delta- relative to the Delta+ state.

Probing the O2 (a 1Delta g) photofragment following ozone dissociation within the long wavelength tail of the Hartley band

The technique of resonance enhanced multiphoton ionization (REMPI) has been used in conjunction with time-of-flight mass spectrometry (TOFMS), to investigate the dynamics of ozone photolysis in the long wavelength region of the Hartley band (301-311 nm). Specifically, both the translational anisotropy and the rotational angular momentum orientation of the O(2) (a (1)Delta(g); nu=0, J=16-20) fragments have been measured as a function of photolysis wavelength. Within this region, the thermodynamic thresholds for the formation of these products in combination with O ((1)D(2)) are approached and passed, and consequently these studies have allowed an investigation into the effects on the dynamics of slowing fragment recoil velocities and the increasing importance of vibrationally mediated photolysis. The determined beta parameters for all the J states probed follow a similar trend, decreasing from a value typical for the initial (1)B(2)<--(1)A(1) excitation responsible for the Hartley band [for example, beta=1.40+/-0.12 for the O(2) (a (1)Delta(g); J=18) fragment], to a much lower value beyond the thermodynamic threshold for the fragment's production (for example, beta=0.63+/-0.19 for the J=18 fragment following photolysis at 311 nm). This trend, similar to that observed when probing the atomic fragment in a previous set of experiments, [Horrocks et al., J. Chem. Phys. 125, 133313 (2006); Denzer et al., Phys. Chem. Chem. Phys. 16, 1954 (2006)] is consistent with the photodissociation of vibrationally excited ozone molecules beyond the threshold wavelengths and we estimate approximately 1/3 of this to be from excitation in the nu(3) asymmetric stretching mode. These observations are substantiated by the values of the beta(0) (2)(2,1) orientation moment measured, which for photolysis at 301 nm are negative, indicating that a bond opening mechanism provides the key torque for the departing O(2) fragment. The orientation moment becomes positive again for photolysis beyond threshold, however, as the increasing impulsive dissociation again begins to dominate the nature of the rotation of the departing molecular fragment. In addition, a (2+2) REMPI scheme has been utilized to probe the O(2) (a (1)Delta(g)) "low" J fragments, where the majority of the population resides following photolysis within this region. The REMPI-TOFMS technique has been used to confirm the rotational character of a spectral feature through examination of the signal line shapes obtained using different experimental geometries. The dynamical information subsequently obtained, probing the "low" J O(2) (a (1)Delta(g)) fragments on these rotational transitions, has unified previous translational anisotropy results obtained by detecting the O ((1)D(2)) atomic fragment with data for the O(2) (a (1)Delta(g); J=16-20) fragments.

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.

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.

REMPI-TOF studies of the translational anisotropy and the polarization of the O (1D2) photofragment angular momentum following ozone photolysis at 298 nm

The translational anisotropy and the polarization of the electronic angular momentum of the O ((1)D2) fragment produced from the 298 nm photodissociation of ozone have been determined using resonance enhanced multiphoton ionization (REMPI) in conjunction with time-of-flight mass spectrometry (TOFMS). The translational anisotropy parameter beta, which is necessarily averaged over the O2 co-fragment rotational distribution, is measured to be 1.08 +/- 0.04. This is consistent with that expected for the (1)B2 <-- (1)A1 transition within an impulsive model if the tangential velocity associated with the zero point motion of the bend is constricted to opening the bond angle. Molecular frame polarization parameters of rank up to k = 4 have been extracted for the O ((1)D2) fragment and the calculated m(J) populations show a strong preference for the absolute value(m(J)) = 1 states. A small coherence term is also observed, a manifestation of the nuclear geometry of the dissociating molecule and the existence of possible non-adiabatic processes in the exit channel. The orientation associated with the mapping of the photon helicity onto the O ((1)D2) electronic angular momentum distribution was observed to have been quenched. However, the parameter gamma1', which describes the contribution to the orientation from a coherent superposition of a parallel and perpendicular excitation where the photofragment angular momentum lies perpendicular to both the recoil velocity and to the transition dipole moment, was determined to be -0.06.

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.

Theory of the photodissociation of ozone in the Hartley continuum: Potential energy surfaces, conical intersections, and photodissociation dynamics

Ab initio potential energy and transition dipole moment surfaces are presented for the five lowest singlet even symmetry electronic states of ozone. The surfaces are calculated using the complete active space self consistent field method followed by contracted multireference configuration interaction (MRCI) calculations. A slightly reduced augmented correlation consistent valence triple-zeta orbital basis set is used. The ground and excited state energies of the molecule have been computed at 9282 separate nuclear geometries. Cuts through the potential energy surfaces, which pass through the geometry of the minimum of the ground electronic state, show several closely avoided crossings. Close examination, and higher level calculations, very strongly suggests that some of these seemingly avoided crossings are in fact associated with non-symmetry related conical intersections. Diabatic potential energy and transition dipole moment surfaces are created from the computed ab initio adiabatic MRCI energies and transition dipole moments. The transition dipole moment connecting the ground electronic state to the diabatic B state surface is by far the strongest. Vibrational-rotational wavefunctions and energies are computed using the ground electronic state. The energy level separations compare well with experimentally determined values. The ground vibrational state wavefunction is then used, together with the diabatic B<--X transition dipole moment surface, to form an initial wavepacket. The analysis of the time-dependent quantum dynamics of this wavepacket provides the total and partial photodissociation cross sections for the system. Both the total absorption cross section and the predicted product quantum state distributions compare well with experimental observations. A discussion is also given as to how the observed alternation in product diatom rotational state populations might be explained.

Photolysis wavelength dependence of the translational anisotropy and the angular momentum polarization of O2(aΔg1) formed from the UV photodissociation of O3

The translational anisotropy and angular momentum polarization of the O(2)(a (1)Delta(g),v = 0;J = 15-27) molecular photofragment produced from the UV photodissociation of O(3) in the range from 270 to 300 nm have been determined using resonance-enhanced multiphoton ionization in conjunction with time-of-flight mass spectrometry. At the shortest photolysis wavelengths used, the fragments exhibit the anisotropic vector correlations expected from a prompt dissociation via the (1)B(2) <--(1)A(1) transition. Deviations from this behavior are observed at longer photolysis wavelengths with, in particular, the angular momentum orientation showing a significant reduction in magnitude. This indicates that the dissociation can no longer be described by a purely impulsive model and a change in geometry of the dissociating molecule is implied. This observation is substantiated by the variation of the translational anisotropy with photolysis wavelength. We also observe that the bipolar moments describing the angular momentum polarization of the odd J states probed are consistently lower in magnitude than those of the even J states and that this variation is observed for all photolysis wavelengths.

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

The triplet channel in the photodissociation of ozone in the Hartley band, O3 + hv-->O(3P) + O2(3sigma), is investigated by means of a classical trajectory surface hopping method using ab initio diabatic potential energy surfaces for the B and the R states. Because of the strong recoil in the R state along the breaking O-O bond, O2(3sigma) is produced with a high rotational energy. The nonadiabatic transition probability depends markedly on the coordinate along the crossing seam. As a consequence a unique correlation is found between the internuclear geometry at the crossing and the final vibrational state of O2(3sigma). The calculated distribution of the translational energy is in good accord with the measured distribution.

The Huggins band of ozone: Assignment of hot bands

The "hot bands" of the Huggins band of ozone are assigned, in both the 218 K and the 295 K spectrum. The assignment is based on intensities calculated with three-dimensional vibrational wave functions for the electronic ground state (X) and the excited state (B). The hot-band structures in the 218 K spectrum all can be assigned to transitions starting from vibrational states with one quantum of stretching excitation in the ground electronic state. The 295 K spectrum shows new structures, which are due to transitions originating from vibrational states in the X state with two quanta of excitation of the stretching modes--despite very small Boltzmann factors. All structures in the low-energy range of the 295 K spectrum, even the very weak ones, thus can be uniquely interpreted. The significance of hot bands results from the strong increase of Franck-Condon factors with excitation of the stretching modes in both the lower and/or the upper electronic states, whose equilibrium bond lengths differ significantly.

Theory of the photodissociation of ozone in the Hartley continuum; effect of vibrational excitation and O(1D) atom velocity distribution

The effect of vibrational excitation on the photodissociation cross section of ozone in the Hartley continuum is examined. The calculations make use of newly computed potential energy and transition dipole moment surfaces. The initial vibrational states of the ozone are computed using grid based techniques and the first few ab initio computed vibrational energy level spacings agree to within 10 cm(-1) with experimental values. The computed total absorption cross sections arising from different initial vibrational states of ozone are discussed in the light of the nature of the transition dipole moment surface. The computed cross section for excitation from the ground vibrational-rotational state is in good agreement with the experimentally measured cross section. Excitation of the asymmetric stretching vibration of ozone has a marked effect on both the form and magnitude of the photodissociation cross section. The velocity distributions of highly reactive O(1D) atoms arising from the photodissociation process in different wavelength ranges is also presented. The results show that the O(1D) atoms travel with a most probable translational velocity of 2.030 km s(-1) corresponding to a translational energy of 0.342 eV or 33.0 kJ mol(-1).

The Huggins band of ozone: A theoretical analysis

The Huggins band of ozone is investigated by means of dynamics calculations using a new (diabatic) potential energy surface for the 3 (1)A'(1B2) state. The good overall agreement of the calculated spectrum of vibrational energies and intensities with the experimental spectrum, especially at low to intermediate excitation energies, is considered as evidence that the Huggins band is due to the two C(s) potential wells of the 1B2 state rather than the single C2v well of the 2 (1)A'(1A1) state. The vibrational assignment of the "cold bands," based on the nodal structure of wave functions, on the whole supports the most recent experimental assignment [J. Chem. Phys. 115, 9311 (2001)]. The quantum mechanical spectrum is analyzed in terms of classical periodic orbits and the structure of the classical phase space.

The Huggins band of ozone: Unambiguous electronic and vibrational assignment

The Huggins band of ozone is investigated by means of exact dynamics calculations using a new (diabatic) potential energy surface for the (1)B(2) state. The remarkable agreement with the measured spectrum strongly suggests that the Huggins band is due to the two C(s) potential wells of the (1)B(2) state. The vibrational assignment, based on the nodal structure of wave functions, supports the most recent experimental assignment.