On neglecting Coriolis and related couplings in first-principles rovibrational spectroscopy: considerations of symmetry, accuracy, and simplicity

The rotation-vibration (Coriolis) coupling contribution to variationally computed rovibrational energy levels is investigated, employing triatomic AB\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}$$\end{document}2 molecules as models. In particular, calculations are performed for H\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}_{2}$$\end{document}2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${}^{16}$$\end{document}16O, across a range of vibrational and rotational excitations, both with and without the Coriolis contribution. A variety of different embedding choices are considered, together with a hierarchy of increasingly severe approximations culminating in a generalized version of the so-called “centrifugal sudden” method. Several surprising and remarkable conclusions are found, including that the Eckart embedding is not the best embedding choice.

where the vibration-rotation (VR) coupling contribution to the kinetic energy operator (KEO)  T 1,2 -i.e.,  T VR -is often called the "Coriolis coupling" (CC) term. Since variational computation [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19] of rovibrational energy levels utilizing VR Ĥ are expensive, both in terms of memory and CPU time, it is highly useful that the rotational quantum number J is a good quantum number 20 and thus the matrix representation of Ĥ VR is block-diagonal in J. Then it is natural to consider KEO approximations that separate rotational and vibrational motions within these blocks in some form. Ideally, such approximations would apply to molecules with complicated internal motion 21-29 , as well as to semirigid ones. Hereby we explore in a joint analytical and numerical treatment-for the first time to our knowledge-the impact on accuracy and numerical efficiency that results when T VR  is neglected, employing different vibrational coordinates and different embeddings of the molecule-fixed axes. Such a study has been called the "logical next step" needed to confirm or refute theoretical hypotheses put forth in previous work 30 . Additionally, we examine what happens when parts of  T R are also discarded. Since V  depends only on vibrational coordinates, one may imagine that neglecting T VR  from Eq. (1) will "decouple" the rotational-vibrational problem, thereby leading, e.g., to a significant reduction in computational effort. As it turns out (see below), further approximations must be applied to realize such savings. Finally, we introduce a natural hierarchy of approximations, culminating in a generalized version of the centrifugal sudden (CS) method [31][32][33][34][35][36][37][38][39] , for which our analysis provides a nice theoretical framework.
In this paper special attention is given to the choice of embedding in first-principles variational rovibrational computations 40 . If the full VR Ĥ is used, converged numerical rovibrational energy levels should be independent of the choice of the vibrational coordinates and the embedding. This is, however, not the case when  T VR is neglected, since  T VR itself depends on the embedding (although not on the vibrational coordinates, as is clear from abstract operator notation). Some embeddings are therefore better than others, in terms of having smaller deviations between the exact and the approximate computed energies. Note that it is impossible to define a frame in which the Coriolis coupling vanishes over the whole configuration space 40 . We therefore explore the role of embedding on the magnitude of the energy contributions related to the Coriolis term at every rung of the approximation hierarchy. We restrict consideration to triatomic AB 2 molecules, using H 2 16 O as our canonical test case. Furthermore, a subset of just three, "linear" embeddings is considered here, as these have been found to perform best in our numerical tests 41 .
The present study builds on previous theoretical work on rotational and vibrational coordinate separation 40,[42][43][44] and coupling 2,30,45,46 . In particular, Sutcliffe and Tennyson derived general rovibrational Hamiltonians in terms of axis embeddings for triatomic molecules 47,48 . Mardis and Sibert 49 derived a Casimir-bond operator, whereby the CC term is zero at equilibrium. Wei and Carrington 50 investigated Eckart embeddings 51,52 -and their bond and bisector counterparts-for triatomics using Radau, valence, and Jacobi vibrational coordinates. The same authors derived triatomic Eckart-embedded Hamiltonian operators for valence 30 and Radau 53 coordinates. Wei and Carrington 30 investigated the properties of both operators, and compared the Coriolis coupling in the different operators focusing on the G VR tensor elements. They concluded that Eckart embedding is the best choice when T VR  is neglected 53 , a basic assumption, which has never since been questioned. Furthermore, they also claim 30 that although they made the first step at discussing the different operators' relative advantages, to really compare their efficacy it is necessary to calculate ro-vibrational energy levels 30 , a task performed during this study. In 1974, McGuire and Kouri developed the "j z conserving" centrifugal sudden approximation in a Jacobi coordinate framework, for three-atom (atom+diatom) quantum scattering calculations 32 ; in the same year, Pack published his paper on related "sudden approximation" methods 31 , specifically comparing space-fixed and body-fixed formulations. The present study generalizes the centrifugal sudden approximation.

triatomic AB 2 molecules: coordinates and embeddings
Consider a triatomic AB 2 molecule with a reference geometry of C v 2 point-group symmetry (Fig. 1). Upon removing the center-of-mass motion, six independent coordinates remain. Three of these are vibrational coordinates, so that in effect there is a (local) three-parameter family of possible embedding choices. Furthermore, triatomic systems are always planar (except for collinear geometries), effectively reducing the range of local embedding choices to just a one-parameter family.
All linear embeddings behave the same with respect to symmetric vibrational displacements. We thus take x ∆ and ∆y of Fig. 1 as the two symmetric vibrational coordinates, describing symmetric stretch and bend motions. Additionally, for all linear embeddings, the pure asymmetric stretch vibrational coordinate corresponds to a linear displacement by the distance ∆, in a direction for B 1 /B 2 (H 1 /H 2 in Figs. 1 and 2, reflecting the fact that our AB 2 test molecule is H 2 16 O) that is at an angle ε below/above the  x axis (see Fig. 2). Different linear embeddings therefore differ only with respect to the value of the angle ε. Only three of the many possible ε choices are considered here: Eckart embedding (EE), "Radau bisector" embedding (RBE), and "valence bisector" embedding (VBE), defined as products of a mass factor and a geometry factor:    16 O, M = 15.990 526 00 u). This is true regardless of x ∆ and ∆y, because the geometric (i.e., second) factors in the right hand side of the two equations are identical. It is the mass (i.e., the first) factor that is identical in Eqs. (2) and (3), whereas the RBE geometric factor has been modified from the EE form to include the symmetric displacements. Note that pure asymmetric stretching is the same for both EE and RBE, since this corresponds to ∆ = ∆ = x y 0. Furthermore, in the vicinity of the reference geometry, EE and RBE are locally equivalent to each other, but different from VBE. These characteristics have important repercussions, as discussed below.
Note that in order to qualitatively compare the three linear embeddings chosen, it is advantageous to choose the same set of vibrational coordinates ( x ∆ , ∆y, and ∆). This is done in this section and the next. Nevertheless, the actual first-principles numerical computations of this study were carried out using valence and Radau internal coordinates (see below). As discussed, the final computed eigenvalues-whether with or without T VR  -should not depend on the choice of vibrational coordinates, although this statement will be tested explicitly.

classical and quantum Hamiltonians & approximations
The three contributions to the rovibrational KEO of Eq. (1) come from the blocks of the well-known covariant and contravariant G and g tensors of nuclear-motion theory 54,55 The asymmetric vibrational coordinate, ∆, and the embedding angle, ε, used to define and compare the three linear embeddings of this work. A displacement ∆ gives rise to asymmetric C s geometries. The bodyfixed z axis is normal to the molecular plane, the x and y axes vary with the embedding, depending on the value of the embedding angle, ε. It is often convenient to replace the asymmetric displacement parameter ∆ with cos µ ε = ∆ .
www.nature.com/scientificreports www.nature.com/scientificreports/ where each block in Eq. (5) is 3 3 × . For certain geometries and embeddings, it may be the case that the VR coupling, g VR , vanishes. In such cases g becomes block diagonal, so that: -- Thus, G VR also vanishes and G R is "geometric"-meaning simply that it is the inverse of the moment-of-inertia tensor, g R . More generally-i.e., when g VR and G VR are not zero-the VR coupling modifies the form of G R as follows: Thus, the three rotational constants, defined as the eigenvalues of G R , are no longer equal to the inverses of the three moments of inertia 56,57 . Note that G R as defined by Eq. (7), and alsog R 1 , are manifestly independent of the choice of vibrational coordinates. This is not true of the tensor components of G VR , which therefore serve as unreliable indicators of the true magnitude of the VR coupling, despite the fact that they have been used in this way in the past 30,53 . Instead, we propose using the tensor --G g ( ) R R 1 for this purpose. Not only is this tensor vibrational-coordinate-independent, but its Frobenius norm, F , is also independent of rotations. The structure of triatomic molecules remains always planar, leading to a block-diagonal planar (xy) and perpendicular (z) structure for g R and G R : Regardless of the embedding, g zz R is always equal to the perpendicular moment of inertia, I I I z x y = + . In contrast, the individual g xy tensor elements vary with the embedding, although the two eigenvalues of g xy are always equal to the two planar moments, I x and I y .
As to G R , there is more variability, since this tensor is in general not geometric. Nevertheless, it can be shown that for all geometries and embeddings the G xy block is always geometric-meaning that G g xy xy 1 = -, and the two planar rotational constants are always Thus, the only place where the non-geometric character of G R can manifest is in the third rotational constant, , comes about because of non-zero Coriolis coupling. It can be shown that only the z components of the Coriolis coupling-i.e., the G iz VR tensor elements, where i is the index of the vibrational coordinates-are non-zero. Note that A z need not in general be the smallest rotational constant, though certainly it is smallest in the geometric case.
For all linear embeddings, G R has the general form (9) Note , as predicted. However, even for C v 2 geometries, G zz R depends on the embedding. Specific forms may be obtained by substituting tanε to Eq. (9) via Eqs. (2), (3), or (4), as appropriate. Doing so reveals something quite special about the RBE-namely, Table 1 for numerical examples). Thus, we arrive at the first important conclusion of this study, namely that RBE is the only linear embedding for which Coriolis coupling vanishes for all C v 2 geometries. This property is in principle discernible from the form of the operator as derived in Eqs. (3)-(19) of ref. 53 ; however, it appears not to have been noticed previously. Indeed, ref. 53 even claims that the Eckart embedding is superior to RBE. For general molecules, it is well known that Coriolis coupling in the Eckart embedding always vanishes at least at one point, i.e., at the reference geometry. In the special case of AB 2 molecules, however, Coriolis coupling vanishes across a one-parameter family of geometries, defined by the pure symmetric stretch motion 53 . On the other hand, the VBE G R is not geometric-and so Coriolis coupling does not vanish-even at the reference geometry itself.
Having discussed the structure of tensor G, the next task is to construct the corresponding Hamiltonian operator and its matrix representation. The choice of vibrational basis set/vibrational coordinates is independent of the embedding and need not be considered further, except for verification purposes. For the rotational space, the usual 58 Wigner rotation function basis, JKM , can be used. This requires specification of the body-fixed axes a b ( , , ) ĉ   (with ĉ associated with K) in addition to the embedding itself. For all linear embeddings, and all exact and approximate Hamiltonians considered here, the rotational quantum numbers J and M are both rigorously good, but there is in general coupling with respect to K, giving rise to a K-block-pentadiagonal structure (i.e., non-zero matrix elements correspond to -≤ The overall parity, p 1 = ± , is also a good quantum number. Through symmetry adaptation, the corresponding = M 0 Hamiltonian matrix, ∼ H J , decouples into positive-and negative-parity symmetry blocks (replacing -≤ ≤ J K J with K J 0 ≤ ≤ ), thereby effectively reducing the basis size by a factor of two. In general, an n-fold reduction of the basis gives rise to a computational saving of n 2 . AB 2 systems are also characterized by permutation symmetry, 1 ε = ± , associated with B 1 -B 2 exchange; this leads to a further factor-of-two reduction (i.e., = n 4), with respect to even-and odd-K values.
The KEO approximations considered introduce additional Hamiltonian symmetries, which increase with the severity of the approximation 59 . The sequence of approximations is as follows. The Coriolis-free approximation (CFA) is the result of neglecting G VR . The CFA Hamiltonian then becomes a sum of pure vibration and pure rotation contributions, R -although this Hamiltonian is still not separable, because  T R depends parametrically on the vibrational coordinates. If, in addition, we set G 0 xy R = , we obtain the diagonal G R approximation (DGRA). This introduces a new permutation symmetry-together with an almost good quantum number that serves as an excellent state label. Finally, the generalized CS approximation (GCSA) is obtained by discarding all remaining K coupling in the JKM representation. As a consequence, K now also becomes a "good" quantum number. Note that the centrifugal sudden approximation has been widely used in the quantum dynamics field [31][32][33][34][35][36][37][38][39] , in the context of Jacobi and Radau coordinates and related embeddings 60 . To the best of our knowledge, we are the first to generalize the centrifugal sudden approximation for arbitrary embeddings.

Results and Discussion
Numerical determination of the rovibrational energy levels of H 2 16 O was based on the BT2 potential energy surface (PES) 61 of H 2 16 O and the GENIUSH code 3,18,19,62,63 -for the exact KEO as well as all three approximations described above, across a wide range of embeddings and rotational and vibrational excitations. In what follows, we focus only on the three linear embeddings, on J 1 = and J 10 = , and on the lowest few vibrational parent states. For all Eckart embedding calculations, results were computed using both valence and Radau internal coordinates. In all such cases, the results were found to be identical, as expected.
The Coriolis-free approximation (CFA). The Coriolis-free approximation does not lead to new symmetries; thus, it provides no significant numerical advantages over exact calculations, although it may be useful for state labeling. On the other hand, being the least severe approximation, CFA is expected to be the most accurate.  Table 1. Numerical values of the G R and G VR tensor elements for all of the embeddings studied in this paper. The only non-zero elements, the z components of the Coriolis coupling-i.e., the G iz VR tensor elements-are shown here, in cm 1 -. The numbers provided correspond to the C v 2 reference equilibrium structure of H 2 16 O, r e = 0.957 820 Å and e θ = 104.500  in valence coordinates, and several symmetrically and asymmetrically distorted geometries with either C v 2 or C s point-group symmetry. The notation "C v 2 (x, y)/C s (x, y)" refers to a symmetric/asymmetric stretch distortion of x % and a bend distortion of y % relative to the reference structure. FNGR is the Frobenius norm ( www.nature.com/scientificreports www.nature.com/scientificreports/ In the Eckart embedding, G 0 VR = for the reference geometry (usually taken as the global minimum of the PES); thus, we expect the greatest accuracy for the pure rotational states, i.e., those corresponding to the vibrational ground state, (v v v 1 2 3 ) = (0 0 0), employing the canonical ordering of the vibrations. For VBE, in contrast, ≠ G 0 V R , even at the reference geometry. Nevertheless, chemical intuition suggests the importance of the VBE picture, and the B-A-B bisector is certainly relevant given the identical B atoms. This suggests that VBE will be less accurate than EE for the lowest-lying energy levels, but may provide greater accuracy further up in the spectrum, particularly for symmetric vibrational coordinate excitations.
The numerical results on H 2 16 O bear out all of these predictions. We see from Table 1 that the EE G VR tensor elements increase from zero to 8.7 cm 1 under a symmetric vibrational coordinate displacement, whereas the corresponding VBE values vary over a narrower range, 0.7-1.7 cm -1 . Likewise, discrepancies in the computed CFA rovibrational energies (see Table 2 and Fig. 3 for J 1 = and Table 3 and Fig. 4 for = J 10) show a highly marked increase with vibrational excitation that is more pronounced for EE than for VBE. On the other hand, = J 1 EE CFA errors for (0 0 0) are as small as 0.01 cm -1 , which is remarkable. Moreover, the EE CFA description for the symmetric-stretch fundamental, (1 0 0), is nearly as accurate as for (0 0 0)-reflecting the aforementioned vanishing Coriolis coupling for this motion. As a rule, however, excitations lead to rapid growth of the errors, with the EE error for (0 1 0) increasing to 0.22 cm -1 , and that of = J 10 (0 0 0) around 1.25 cm -1 . In contrast, VBE errors for = J 1 (0 0 0) are around 0.06 cm 1 -, and for most rovibrationally excited states are significantly smaller than EE errors, especially for the bending excitations.
In Table 1, explicit G VR and G R tensor elements are provided for the Eckart embedding, for both valence and Radau vibrational coordinates, at several different geometries. The G VR /G R tensor elements are seen to be vibrational coordinate dependent/independent, as expected. Table 1 also lists values for the Frobenius norm of F ), which are independent of the choice of vibrational coordinates. The matrix elements for RBE and VBE are also shown. In order to gain a deeper understanding of the manifestation of CC in the studied embeddings, it is worth comparing G g R R

F
-values at different geometries; again, this is taken as a measure of the extent of CC.
For the reference geometry, VBE shows nonzero CC, whereas both RBE and EE show zero CC, as expected. On the other hand, RBE exhibits zero CC across all C v 2 distortions, a considerable advantage of this embedding. This is unlike EE, which has substantially larger CC, even compared to VBE, when there is bending excitation, although for pure symmetric stretching EE has zero CC, as well. For C s geometries, either RBE or EE can be better than the other embedding, depending on the actual distortion. For pure asymmetric stretch distortions, EE is better; however, the EE G g R R

F
-values increase rapidly for combinations with any bend excitation. In contrast, RBE values remain more or less constant with increasing bend excitation, reflecting the fact that bending is . The color-coded results correspond to valence bisector (VBE), Radau bisector (RBE), and Eckart (EE) embeddings. The vibrational (v v v 1 2 3 ) quantum numbers are assigned for each state and they are presented in the form of the resonance polyads, nν or ν δ + n according to the number of stretching (ν) and bending (δ) quanta, where two bending is "equivalent" to one stretching excitation. -values, we observe that these are only slightly larger than RBE for the C s geometries considered.
Of course, the "acid test" for a given embedding is the accuracy of its CFA rovibrational energy levels. Based on --G g R R 1 F values, RBE is expected to be the most accurate choice among the three embeddings considered here-particularly for symmetric stretch and bend excitations. Indeed, this is the case. As Table 2 and Fig. 3 attest, the J 1 = RBE errors can be almost as small as EE errors for those few cases where EE is the best, yet significantly smaller than EE and VBE errors in all other cases. Even for pure asymmetric stretch excitations, RBE errors are about an order of magnitude smaller than EE, as seen from Table 2. For higher J, this trend is even more evident (see Fig. 4).
Overall, all three linear embedding CFA results reproduce the exact values remarkably well. For the lowest 50 rovibrational states, the largest errors are about 1 cm -1 for = J 1, and 20 cm -1 for = J 10. Additionally, for J 1 = -for every single vibrational state and for every embedding-the Coriolis-free prediction for the 1 01 state rovibrational energy is exact. This remarkable finding relates to parity, as explained in the next subsection. the diagonal G R approximation (DGRA). The next rung in our descending hierarchy is the diagonal G R approximation. For all C v 2 geometries (including the reference geometry), G 0 xy R = . Only for large asymmetric displacements do we expect G xy R to become substantial; accordingly, only for excited asymmetric stretch states do we expect to see a large difference from CFA. In actuality, however, the DGRA energy levels are extremely close to their CFA counterparts-much more so than might be expected (see Table 2). Most discrepancies are substantially smaller-and none are significantly larger-than 0.01 cm 1 -, across the full range of vibrational excitations considered. We also note that for all embeddings the computed 1 01 DGRA levels are still exact. This is because G xy R does not contribute to the × 1 1 negative-parity block. As another important conclusion of the present study, from a practical standpoint, there is no reason not to use the diagonal G R approximation, if one is committed to throwing away Coriolis coupling anyway (at least for linear embeddings).
DGRA introduces a new permutation symmetry, which can be very useful in practice. Note that permutation affects rotation and vibration simultaneously. In the vibrational space, permutation changes the sign of the asymmetric stretch displacement, without affecting the symmetric displacements. In the rotational space, even-K rot (0 0 0) (0 1 0) (0 2 0) (1 0 0) (0 0 1) . The results are provided in cm -1 and correspond to valence bisector (VBE), Radau bisector (RBE), and Eckart (EE) embeddings, and they are given relative to the vibrational parent state. Vibrational (v v v 1 2 3 ) and rotational (rot, J K K www.nature.com/scientificreports www.nature.com/scientificreports/ corresponds to one ε value, and odd-K to the other. The effect of permutation on the two spaces is coupled solely through the G xy R contribution. Consequently, without G xy R , one obtains two independent permutation symmetries-a vibrational permutation symmetry ɛ vib , and an independent rotational permutation symmetry ɛ rot . It may be feasible to utilize this extra information during the assignment of rovibrational states.
As a practical benefit, this additional symmetry allows a further two-fold reduction in the basis size, taking us up to n 8 = . Furthermore, because DGRA is so close to CFA-which in turn does an excellent job of modeling the exact rovibrational energies-both ɛ vib and ɛ rot are nearly perfectly good quantum numbers. Indeed, the usual association for AB 2 molecules of even-v 3 vibrational quantum states with even permutation symmetry, and odd-v 3 states with odd permutation symmetry, is in fact a manifestation of vib ɛ , rather than ɛ itself.
The generalized CS approximation (GCSA). The CS approximation can drastically reduce the computational cost, especially for large J. The idea is very simple: set all off-diagonal K blocks in H J ∼ equal to zero. K then becomes a good quantum number-thus, once again, increasing the symmetry of the Hamiltonian. Since the remaining diagonal = ′ K K blocks can be diagonalized separately, the basis size of the problem reduces by a factor of n J (2 1) = + . For J 3 > , this provides a greater computational reduction than the DGRA, for which = n 8. Actually, there are some modest savings for J 2 = and J 3 = , as well, for which not all of the = n 8 irreps are realized in the diagonal G R case.
The discarded ≠ ′ K K blocks are often loosely referred to as "Coriolis coupling", although they are clearly not T VR  . In general, what these off-diagonal blocks represent depends on the embedding, as well as the particular body-fixed axis, ĉ, along which K is projected. It is convenient to choose ĉ to correspond to x , ŷ, or z , although strictly speaking any (global) orientation may be used. For planar molecules, the choice = ĉ z  may be expected to be particularly poor. All three choices,  x ĉ = , = ĉ ŷ, and ĉ z = , are considered here. Previous applications of the CS approximation have been in the context of Jacobi and Radau vector embeddings, for which ∼ H J is tridiagonal. In the present, generalized CS context, we are also discarding a block-pentadiagonal contribution, which may in principle be quite large. Note that GCSA is equivalent to the following prescription: 1. Start with the DGRA matrix representation, and identify the projection axis ĉ and the diagonal tensor elements G aa R and G bb R .  www.nature.com/scientificreports www.nature.com/scientificreports/

Replace both G aa R and G bb
Step 2 automatically results in a symmetric rotor form. It is interesting that this prescription seems to answer an "age-old" debate about how best to go from an asymmetric to a symmetric rotor form 39 : should one average two moments of inertia or two rotational constants? Evidently, the latter is the more correct approach.   16 O in Radau bisector embedding with different generalized centrifugal sudden approximations (GCSA) referring to each of the three choices of ĉ axis along which K is projected. The results are provided in cm -1 and given relative to the vibrational parent state. Vibrational [vib, v v v ( ) 1 2 3 ] and rotational (rot, J K K a c ) quantum numbers are assigned for each state. The differences of the eigenvalues compared to the diagonal G R approximation (DGRA) [∆(GCSA)] are also provided.

Scientific RepoRtS |
(2020) 10:4872 | https://doi.org/10.1038/s41598-020-60971-x www.nature.com/scientificreports www.nature.com/scientificreports/ GCSA results are presented for J 1 = in Table 4, for each of the three choices of ĉ (i.e., CS x , CS y , and CS z ), for the RBE (results for the other linear embeddings are similar). Note that = J 1 is special, in that each GCSA calculation results in one exact diagonal G R eigenvalue per vibrational parent. Although CS z follows the trend of the former approximations, by averaging G xx R and G yy R the 1 10 levels remain unchanged, and CS x yields the smallest errors. The errors are on the order of ±2.5 cm -1 for the zero-point vibration, compared to ±6.7 cm -1 and ±9.3 cm 1 for CS z and CS y . Surprisingly, CS x errors do not increase appreciably with vibrational excitation.

conclusions
This joint analytical and numerical study offers several interesting and occasionally surprising conclusions about various rotational-vibrational approximations, and associated linear embeddings, employing H 2 16 O as a canonical test system: 1. The Eckart 51 embedding is not the best choice of embedding, especially when rotational and/or vibrational excitations are significant. 2. There exist embeddings (e.g., the Radau bisector embedding) for which Coriolis coupling vanishes over an entire two-dimensional vibrational coordinate subspace (i.e., the symmetric vibrational coordinates). 3. The diagonal G R approximation should be used instead of the Coriolis-free approximation, because the computed energy levels are almost exactly as good but the benefits in terms of symmetry and numerical efficiency are substantial. 4. The diagonal G R approximation leads to two distinct and almost-perfectly-good parity quantum numbers, one for rotations and one for vibrations.
Note that conclusions (3) and (4) occur only by virtue of permutation symmetry, but should manifest for larger molecules with at least two identical nuclei. As for (2), it is well known that the Eckart embedding achieves zero Coriolis coupling at any given geometry 51 , and that this can be extended over a one-dimensional subspace using "post-Eckart" embeddings, such as the one due to Sayvetz 64 . That two-dimensional subspaces of this kind are also possible, for systems and embeddings where Coriolis coupling as a whole does not vanish, is a novel discovery, so far as we are aware. Even more surprisingly, the Radau bisector embedding that achieves this is not Eckart-based.
The full ramifications of these findings certainly merit further investigation, especially vis-à-vis larger and more general molecular systems. In an upcoming paper 41 we will address other linear and nonlinear embeddings such as those based on Jacobi coordinates, other types of triatomics such as triangular A 3 molecules, and also larger molecules. The latter class includes "floppy" systems such as ammonia with an inversion tunneling motion.