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.

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.

Topology of conical/surface intersections among five low-lying electronic states of CO2: multireference configuration interaction calculations

Multi-reference configuration interaction with single and double excitation method has been utilized to calculate the potential energy surfaces of the five low-lying electronic states (1)A1, (1)A2, (3)A2, (1)B2, and (3)B2 of carbon dioxide molecule. Topology of intersections among these five states has been fully analyzed and is associated with double-well potential energy structure for every electronic state. The analytical potential energy surfaces based on the reproducing kernel Hilbert space method have been utilized for illustrating topology of surface crossings. Double surface seam lines between (1)A1 and (3)B2 states have been found inside which the (3)B2 state is always lower in potential energy than the (1)A1 state, and thus it leads to an angle bias collision dynamics. Several conical∕surface intersections among these five low-lying states have been found to enrich dissociation pathways, and predissociation can even prefer bent-geometry channels. Especially, the dissociation of O((3)P) + CO can take place through the intersection between (3)B2 and (1)B2 states, and the intersection between (3)A2 and (1)B2 states.