banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Detailed analysis of polyad-breaking spectroscopic Hamiltonians for multiple minima with above barrier motion: Isomerization in HO2
Rent this article for


Image of FIG. 1.
FIG. 1.

Plot of the model potential energy surface with contours representing 2000 cm−1 equipotential lines. The O–O bond length is fixed to the equilibrium value. The oxygen atoms are symmetrically placed on the x axis about zero.

Image of FIG. 2.
FIG. 2.

Comparison of various fits with the energy zero defined as the ZPE level. The isomerization barrier is located at 11 950 cm−1. Several states are labeled based on a visual inspection of the ab initio wave functions.

Image of FIG. 3.
FIG. 3.

Comparison of squared wave functions for the |3 2〉 state from several resonant fits, all of which yield good results. The physical coordinates are those of Fig. 1; however, due to the symmetry of the potential, we present a view encompassing only the X H > 0 portion. The red “X” in the exact plot indicates the location of the transition state. The wave functions for the effective Hamiltonians were constructed using a harmonic basis with frequencies obtained from the zero order fit.

Image of FIG. 4.
FIG. 4.

A comparison of squared wave functions for the |0 11〉 state. The result illustrates the negative effects on the fit by adding too many off-diagonal coupling terms. shows that addition of a single third order diagonal coupling term to a minimal set of off-diagonal resonance couplings yields physically realistic wave functions and good agreement with the energy level.

Image of FIG. 5.
FIG. 5.

A comparison of squared wave functions for the |0 13I〉 state where the “I”designates an isomerizing bend state. Although yields physically realistic results up to 13 500 cm−1, above that energy poor results are seen. Good results are obtained by including additional couplings () or a selected third order diagonal term (). The former option produces poor results for some levels (see Fig. 4), while the latter yields good results for all states below 15 000 cm−1.

Image of FIG. 6.
FIG. 6.

Comparison of squared wave functions and classical POs from the exact potential surface (a–d) and the effective Hamiltonian (e). Strong correlations are seen in (a) and (b) for bend and combination states. Panel (c) shows a strong relation to the 2:1 PO, while panel (d) illustrates the clear difference between a local PO and an isomerizing bend wave function. Although a classical PO on the exact surface cannot be found to match the isomerizing bend, the model Hamiltonians (e) shows that a 3:1 PO matches with the fit wave functions for this state.

Image of FIG. 7.
FIG. 7.

Comparison of the SOS for the exact, , and Hamiltonians, respective columns. The first column is in the physical coordinates of Fig. 1 and the second and third columns are in the action-angle variables of the semiclassical fitting Hamiltonians, hence the different appearance of the SOS. However, colors are consistent between plots and identify similar features. At low energy the phase space is simple, while at intermediate energies several resonance islands have appeared. At high energy, chaos has destroyed several of the island chains. Note that due to the sectioning plane some of the islands are missing; hence, the identity of the resonances was determined from inspection of coordinate space trajectories.

Image of FIG. 8.
FIG. 8.

A portion of the Farey sequence between 1/3 and 1/2, related to the resonances present in the system. Solid lines denote Farey neighbors (e.g., 1/3 and 1/2) and their branchings (e.g., 2/5), while vertical dotted lines divide the tree and show boundaries for which neighbors cannot exist (5/13 and 5/12). Blue members denote the 3:1 portion, while red members denote the 2:1 portion.

Image of FIG. 9.
FIG. 9.

The frequencies for the bend states (|0 n〉) and combination states (|1 m〉). The bend states show a typical anharmonic progression with no evidence of the approach to the barrier, while the combination states show a dip just above the barrier at 11 950 cm−1.

Image of FIG. 10.
FIG. 10.

Polyad phase spheres resulting from the 3:1 (top) and 2:1 (bottom) single resonance Hamiltonians. Stable fixed points associated with the stretch mode and the anharmonically perturbed bend mode are labeled as S1 and S2, respectively. The 3:1 coupling shows an additional stable fixed point, S3, associated with the isomerizing bend states. The 2:1 bifurcation results in a hyperbolic fixed point (H1) at the north pole of the 2:1 polyad phase sphere. There is also an unstable fixed point on the 3:1 sphere between S3 and the north pole.


Generic image for table
Table I.

Optimized parameters of the spectroscopic Hamiltonians.


Article metrics loading...


Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Detailed analysis of polyad-breaking spectroscopic Hamiltonians for multiple minima with above barrier motion: Isomerization in HO2