The rovibrational Aharonov-Bohm effect

Another manifestation of the Aharonov-Bohm effect is introduced to chemistry, in fact to nuclear dynamics and high-resolution molecular spectroscopy. As demonstrated, the overall rotation of a symmetric-top molecule influences the dynamics of an internal vibrational motion in a way that is analogous to the presence of a solenoid carrying magnetic flux. To a good approximation, the low-energy rovibrational energy-level structure of the quasistructural molecular ion H5+ can be understood entirely in terms of this effect.


Introduction
More than 60 years ago, following the footsteps of Ehrenberg and Siday, 1 Aharonov and Bohm 2,3 provided a detailed analysis of the significance of electromagnetic potentials in quantum theory. They made the then surprising claim that in the quantum domain these potentials have a physical and measurable effect on charged particles even in regions where the magnetic field vanishes. The findings of Aharonov and Bohm were verified experimentally. 4,5 While at the beginning the statements of these studies about fields and potentials may have seemed just counterintuitive curiosities of quantum theory, over the years the Ehrenberg-Siday-Aharonov-Bohm (usually simply referred to as the Aharonov-Bohm) effect, and its various analogues, found application in a number of different fields of molecular sciences.
Of particular importance to this study is the so-called molecular Aharonov-Bohm (MAB) effect 6-8 within Born-Oppenheimer theory, 9,10 in which the conical intersection seam, 11,12 corresponding to degenerate electronic states, has an effect on the nuclear motion analogous to the presence of a solenoid carrying magnetic flux. The implications of the MAB effect for both bound states 7,8 and scattering events 13 have been explored.
We are aware of no indication of the significance of the Aharonov-Bohm effect in nuclear dynamics unrelated to the MAB effect. In this paper we introduce the 'rovibrational Aharonov-Bohm effect', occuring in symmetric-top molecules with a separable vibrational degree of freedom (dof). As shown below, the 'rotational' dynamics influence the 'vibrational' dynamics in a way analogous to the coupling of the selected vibrational motion to the field of a magnetic solenoid. We investigate in detail whether this analogy between rovibrational nuclear dynamics and electromagnetic phenomena can help us to understand high-resolution molecular spectra and the related dynamics of molecules.
The treatment of overall rotations coupled to large-amplitude internal motions (sometimes referred to as contortions 14 ) has been the subject of intensive studies in high-resolution molecular spectroscopy. [14][15][16][17][18][19][20][21][22][23][24] Perhaps the most famous among the models is the pioneering Hougen-Bunker-Johns (HBJ) model, 17 describing the rotational-contortional dynamics of triatomic molecules with a large-amplitude dof, the bending. There exist reviews 14,16,22 on the topic of rotation-contortion-vibration Hamiltonians (the contortions include internal rotation, pseudorotation, inversion, torsion, and bending) and their utilization in the understanding of high-resolution molecular spectra. The different Hamiltonians derived, for example, the approximate principal-axis, internal-axis, and rho-axis ones, are relevant to the present discussion, as they contain different forms of kinetic energy and rotational-vibrational (Coriolis) coupling terms obtained via different rotational transformations. Nevertheless, to the best of our knowledge, the central idea of the present paper, the 'electromagnetic analogy' of such Hamiltonians leading to the rovibrational Aharonov-Bohm effect was not considered in these studies.
As an example of a molecular system exhibiting the rovibrational Aharonov-Bohm effect, the low-energy nuclear dynamics of the H + 5 molecular cation is considered. First-principles characterization of the dynamics of H + 5 , as well as of its deuterated derivatives, has been achieved and it offered a number of considerable surprises. [25][26][27] As shown below, the rovibrational Aharonov-Bohm effect dominates the low-energy nuclear dynamics of this quasistructual 28 molecular ion, which is approximately a symmetric top with a nearly-free torsional motion.

The electromagnetic analogy
Let us start our discussion with the standard spectroscopic four-dimensional (4D) rotationcontortion model Hamiltonian 15,16,24,28,29 describing the dynamics of a symmetric-top molecule with a torsional (φ, period 2π) degree of freedom (dof) coupled to the rotational dofs. In Eq. (1),p φ = −i ∂ ∂φ , a depends on the moments of inertia of the internal rotor and of the molecule about the molecular symmetry axis, the potential energyV depends solely on φ and typically has multiple minima, the three components of the overall body-fixed angular momentumĴ are denoted byĴ i , A and B are the effective rotational constants of the symmetric top, and α is the rovibrational coupling strength. The various operators appearing inĤ of Eq. (1) A common next step is to perform the so-called Nielsen transformation, 15 which leads to the apparent elimination of the rovibrational coupling. We will discuss the Nielsen transformation shortly. For now, we note that K, related to the projection of the angular momentum onto the molecule-fixed z axis, is a good quantum number. SettingĴ z = K in Eq.
(2), we get thatĤ and so the Hamiltonian is given by the sum of a 'torsional' Hamiltonian and a rotational kinetic energy contribution, given by We must point out thatĤ K tor , which we are calling the 'torsional' Hamiltonian, in fact depends on the rotational quantum number K through the appearance of the operator p φ − αK 2a . Essentially,Ĥ K tor takes the usual kinetic plus potential form, except for the fact that the torsional momentump φ has been 'modified' by the replacementp φ →p φ − αK 2a . The key insight is that this is precisely the modification one would make in order to couple the torsional motion to a magnetic field, where α 2a plays the role of an electromagnetic vector potential and K is formally identified with the electric charge. In fact, the Hamiltonian in Eq. (4) is of exactly the same form as the Hamiltonian for a particle confined to a ring encircling a solenoid which is carrying magnetic flux, 30 just as in the famous Aharonov-Bohm effect. It is as if the rotational dofs induce an effective magnetic field which couples to the torsional motion. We call this the rovibrational Aharonov-Bohm effect. It can be viewed as a special case, for symmetric tops, of the more general influence of rovibrational coupling on vibrational motion through so-called non-Abelian gauge fields (for more information on this general picture, see Refs. [31][32][33][34][35][36]. As a result, the effect of rovibrational coupling on the torsional motion can be understood by analogy with the familiar effects of a magnetic field. To the best of our knowledge, this is a new perspective on the well-studied Hamiltonian of Eq. (1). As we will see, this electromagnetic analogy allows us to get deeper insight into the rovibrational level structure, and reveals properties which are obscured by the usual separation of nuclearmotion Hamiltonians into rotational, vibrational, and rovibrational terms.

Changes of embedding
Next, let us consider the effect of a change of embedding of the molecule-fixed axes. Suppose that, for our original embedding, we have some rovibrational eigenstate χ K (φ) |JKM with the torsional wavefunction χ K satisfying Imagine now that we change our embedding choice. We only consider embeddings where the body-fixed z-axis is aligned with the symmetry axis of the molecule. At each φ, the old and new embedding must be related by a rotation about the body-fixed z-axis by some angle θ (φ). With respect to the new embedding, the torsional wavefunction χ K will becomeχ K , where Substituting this relation into Eq. (6) gives and, by comparison with Eq. (3), we see that the rotational-energy contribution looks exactly the same in the new embedding, and the only change to the torsional Hamiltonian is the of the electromagnetic vector potential by the φ-derivative of θ. We recognise this as the usual transformation law for an electromagnetic vector potential under a gauge transformation.
Thus, we see that, within the electromagnetic analogy, changes of embedding correspond to electromagnetic gauge transformations.

Eliminating rotation-vibration coupling
The electromagnetic analogy gives us insight into the important issue of the "elimination" of rotation-vibration coupling. 37 Recall that, under a change of embedding specified by θ (φ), the vector potential transforms as in Eq. (9). If we can choose our embedding so that the vector potential becomes zero, then clearly the torsional motion will be separated from the rotational motion. This seems possible, since we can just choose θ (φ) = − αφ 2a and then Eq. (9) tells us that the vector potential will become thus, with respect to this new embedding choice, the torsional Hamiltonian becomeŝ which no longer depends on the rotational quantum number K.
To see this, first recall that the torsional coordinate φ is 2π-periodic. On the other hand, the rotation angles θ (0) and θ (2π) relating the old and new embeddings at φ = 0 and φ = 2π are not the same, but differ by an angle θ (2π) − θ (0) = −απ/a. In other words, θ does not give a single-valued choice of embedding. As a consequence, the torsional wavefunction in the new embedding [see Eq. (7)],χ K , satisfies and so we see that, while the original torsional wavefunction χ K was 2π-periodic in φ, it becomes necessary to allow for modified boundary conditions, Eq. (12), for the new torsional wavefunctionχ K . These new boundary conditions 'compensate' for the multi-valuedness of the new embedding.
For example, for a molecule composed of two identical rotors (for which it turns out that α/a = 1-see the following section for an example) we have that or in other words that χ K satisfies periodic/anti-periodic boundary conditions for K even/odd. The necessity to modify the boundary conditions, as in Eq. (12), is arguably an unpleasant feature of the Nielsen transformation. It is also unnatural from the perspective of the electromagnetic analogy, in which gauge transformations ought to be single-valued. We might ask, therefore, whether it is possible to eliminate the vector potential without having to introduce modified boundary conditions for χ K . In other words, can we construct some single-valued change of embedding which transforms the vector potential away? The electromagnetic analogy makes it clear that the answer is generally 'no'. This is because the electromagnetic vector potential α/2a corresponds to the presence of a magnetic flux and, as is well known in the context of the traditional Aharonov-Bohm effect, this magnetic flux generally has a non-trivial effect on the quantum energy levels. It cannot simply be transformed away. The only exception is when α/2a = n ∈ Z is an integer, since then the single-valued embedding transformation θ (φ) = −nφ yields showing that the effect of the corresponding magnetic flux is no different to if there were no magnetic field at all. This is a flux quantization condition, as is familiar in the Aharonov-

1D(φ) torsion model
which depend on a torsional coordinate φ, which has period 2π, and R is the distance between the midpoints of the two H 2 units. Note that here, and throughout this section, we work in units where = m H = r = 1, where m H is the mass of the hydrogen atom and r is the distance between the two hydrogens forming a rotor.
The time-independent Schrödinger equation for the 1D(φ) torsion model takes the form In previous studies this equation has been split into separate rotational, vibrational and rovibrational contributions but here we have grouped the terms as suggested in the mathematical physics literature. 34 In Eq. (16),p φ = −i∂/∂φ is the momentum conjugate to the torsional coordinate φ and theĴ i are the components of the body-fixed angular momentum, while M is the moment of inertia tensor with respect to the GE molecule-fixed frame, governing the motion along the three rotational dofs, and has the explicit form Note that we have allowed for a torsional potential V (φ) in Eq. (16), which on physical grounds is assumed to satisfy for all φ (recall that φ has period 2π). In addition, the extrapotential term V 2 (φ) is a quantum contribution coming from the φ-dependence of the rovibrational G matrix, 14 which in this case takes the explicit form At this point let us pause to interpret the various terms in Eq. (16). Broadly speaking, the first term on the left-hand side should be thought of as the torsional kinetic energy; the second term is the rotational kinetic energy; the third term is the torsional potential energy, which has been modified by the extrapotential term V 2 (φ). Note that our grouping of terms in Eq. (16) is different to the usual conventions adopted in the nuclear-dynamics literature.
For example, note the appearance of the operatorp φ − 1 2Ĵ z in the first term, reminiscent of the electromagnetic coupling considered in the previous section. We will now make the connection to the previous section, and to the Aharonov-Bohm effect, more explicit.

Symmetric-top approximation
For applications to H + 5 , it is a very good approximation to take R large. 25 So let us expand in 1/R. To leading order, and Discarding terms which are of order 1 R 4 or higher, M −1 becomes diagonal and V 2 vanishes and so Eq. (16) reduces to We call this (large-R) approximation the symmetric-top approximation, since the rotational kinetic energy term is now that of a symmetric top. In fact, the Hamiltonian in Eq. (22) So, by the results of the previous section, the rovibrational energy levels are of the form where, in this case, the rotational energy contribution is and the torsional energy contribution is an eigenvalue of the torsional Hamiltonian Note that the electromagnetic couplingp φ − K 2 appearing in the torsional Hamiltonian of Eq. (26) takes a particularly simple form, due to the special value of the ratio α a = 1 in this case.
We can now take advantage of the electromagnetic analogy to better understand the torsional Hamiltonian of Eq. (26), which is simply the Hamiltonian for a charged particle, confined to a ring, encircling a solenoid carrying magnetic flux. The first insight that this analogy gives us is that the torsional energy levels only depend on whether K is even or odd, and not on the precise value of the integer K. We prove this now.
Suppose that E tor is an energy eigenvalue ofĤ K tor , with ψ a corresponding eigenstate: Then consider ψ N := exp (iN φ) ψ where N ∈ Z. We have that are all the same. Similarly, all have identical energy levels. So, in fact, the energy levels ofĤ K tor are only dependent on whether K is even or odd, and not the particular value of the integer K. In the electromagnetic anology, this can be thought of as a kind of flux quantization condition.
Another property of the Hamiltonian of Eq. (26) is that, for K odd, all torsional energy levels are doubly degenerate. This is remarkable, as we might expect the torsional potential energy V to cause a splitting of degenerate levels. For K even, this is precisely what happens -levels which are degenerate for V = 0 are split when V is turned on. However, for K odd, the special value of the magnetic flux carried by the solenoid ensures that the degeneracy of torsional energy levels persists even when V = 0. This doubling of levels is a known effect in electromagnetic theory 38 and we give a brief proof here.
Recall that V is symmetric under φ → φ+π soĤ K tor is also symmetric under φ → φ+π. As a consequence, we can choose torsional energy eigenstates to be even or odd under φ → φ+π.
Suppose that E tor is an energy eigenvalue ofĤ K tor , with ψ a corresponding eigenstate which Now define a new wavefunctionψ (φ) = exp(iKφ)ψ (−φ). Then direct calculation, remembering that V is symmetric also under φ → −φ, shows that soψ is also an energy eigenstate with the same energy eigenvalue as ψ. Note also that so, if K is an odd integer, thenψ is odd under φ → φ + π. So we have a pair of degenerate torsional energy eigenstates ψ andψ, which are clearly linearly independent because one is even and the other is odd under φ → φ + π. In other words, every even eigenstate automatically comes with an odd partner. By the same argument, every odd eigenstate comes with an even partner and so we see that all levels, for K odd, are doubly degenerate.
We will come back to these properties shortly to explain how they correspond to wellknown properties of the nuclear wavefunction in the presence of conical intersections as studied in the molecular Aharonov-Bohm effect.
In summary, the electromagnetic analogy has given us deeper insight into the energy-level structure of the rovibrational states of H + 5 . In the general case with non-zero V , states with K even have energy levels where E even tor takes one of a discrete set of non-degenerate torsional energy eigenvalues. States with K odd have energy levels where E odd tor takes one of a discrete set of doubly degenerate torsional energy eigenvalues.
H + 5 is an example of a system where the rotational and vibrational dofs are strongly coupled, so they cannot be treated independently. This is the reason why the torsional energy levels depend on the rotational quantum number K. However, we have found that the torsional energy levels only depend on K through whether it is even or odd, and not on its precise value. In particular, they only have to be computed for two values of K, say for K = 0 and K = 1, and then all the rovibrational energy levels can be deduced. This is obvious from the perspective of the rovibrational Aharonov-Bohm effect, as a kind of magnetic flux quantization condition. Even more strikingly, every torsional energy level is doubly degenerate for K odd. Again, this double degeneracy is obvious from our new perspective and corresponds to a known electromagnetic effect, 38 but it is obscured by the usual decomposition of the rotation-vibration Hamiltonian into rotational, vibrational, and rovibrational coupling terms.
Before finding the rovibrational energy levels, let us elaborate on the connection between the molecular Aharonov-Bohm effect and the rovibrational Aharonov-Bohm effect we are considering here.

Analogy with conical intersections
For chemists, the most familiar manifestation of the Aharonov-Bohm effect is probably the molecular Aharonov-Bohm effect (MAB), in which a nuclear wavefunction encircling a conical intersection of potential energy surfaces picks up a minus sign. In Born-Oppenheimer theory, 9,10 one starts by choosing an electronic ground-state wavefunction for each nuclear configuration. There is a choice of complex phase involved at each nuclear configuration, since the electronic ground-state is only defined up to a phase. Often, we demand that the electronic wavefunction is real-valued (for simplicity, let us imagine our system is timereversal invariant and that we are ignoring spin-orbit coupling). Then, there is still a choice of phase, but there are only two choices and they differ by a factor of −1 (a phase of π).
Suppose we choose a real-valued electronic ground-state wavefunction at some nuclear configuration, and we follow this choice continuously along a closed loop that encircles a conical intersection, all the time requiring our electronic ground-state wavefunction to be real valued. As is well-known, the wavefunction picks up a phase of π upon returning to the starting point. Single-valuedness of the molecular wavefunction then requires that the nuclear wavefunction also pick up a minus sign. Another way of characterizing this phase choice for the electronic wavefunctions is that this is the phase choice for which the MTB vector potential 6 vanishes.
In our nuclear-motion problem, we also have a choice: for each value of the torsional coordinate φ, we have to choose an embedding of the molecule-fixed axes. We have been working with the GE throughout this section, but we are free to transform to some other embedding by specifying a function θ (φ) which gives the rotation about the body-fixed zaxis (the symmetry axis of the symmetric top) which relates the old and the new embedding.
Specifying this rotation θ (φ) is much like specifying the phase of the electronic ground-state wavefunction in Born-Oppenheimer theory.
So let us start at φ = 0 and set θ (0) = 0. Now we imagine increasing φ and, for each φ, choosing θ (φ) so that the rotation-vibration coupling vanishes. Note that this choice is analogous to the choice of a real-valued electronic wavefunction in Born-Oppenheimer theory: in the Born-Oppenheimer case, the MTB vector potential vanishes while in this case, the vector potential associated with rovibrational coupling vanishes. In this way our original GE is transformed to a new embedding, the BE, which is depicted in Fig. 1.
It turns out that when we get to φ = 2π, the BE embedding has rotated by an angle of π (with respect to the body-fixed z-axis) compared to the BE embedding at φ = 0. This rotation by π is just like the phase of π picked up by a real-valued electronic wavefunction upon encircling a conical intersection. Now, what does this rotation by π mean for quantum states? If our rotational state |JKM has body-fixed angular momentum projection K even, then a rotation by an angle π does nothing to the state and so the state is insensitive to this rotation. If our rotational state has K odd, however, then it picks up a minus sign.
In turn, for the rovibrational wavefunction to be single-valued, we must have that for K odd the corresponding torsional states pick up a minus sign. This is just as we found in the previous section: elimination of rovibrational coupling is possible, but only if one is prepared to modify the boundary conditions on the torsional wavefunction. With respect to the BE, K even/odd torsional wavefunctions satisfy periodic/anti-periodic boundary conditions. We note that, as a consequence, nuclear motion codes will only yield correct rovibrational energy levels when working in the BE if the boundary conditions are modified compared to the GE.
Thus, the connection between our rovibrational problem and the treatment of conical intersections is clearest when we work with the BE. Let us temporarily adopt this embedding.
Then we can make two direct connections between our rovibrational problem and the familiar treatment of conical intersections.
First, for K odd, the torsional wavefunction picks up a minus sign (like when there is a conical intersection) while, for K even, the torsional wavefunction picks up no minus sign (like when there is no conical intersection) under φ → φ + 2π. This is another way of explaining why the torsional energy levels only depend on whether K is even or odd, and not on the precise value of K. It shows that K even/odd just corresponds to periodic/antiperiodic boundary conditions, i.e. conical intersection/no conical intersection.
Second, it is known that one has to introduce 'double groups' to classify nuclear states in the presence of conical intersections, 39 and these larger symmetry groups can lead to higher degeneracies than occur in the absence of conical intersections. This can be viewed as the reason for the double degeneracy of K odd torsional energy levels. Indeed, the symmetries of the torsional potential V (φ) under φ → φ+π and φ → −φ generate a symmetry group for the torsional dynamics isomorphic to C 2v , and this group only has one-dimensional irreducible representations (irreps) since it is Abelian. The double group is isomorphic to C 4v , a non-Abelian group, with two-dimensional irreps which correspond to the doubly degenerate levels.

Rovibrational energy-level structure via first-order perturbation theory
Now we return to the calculation of the rovibrational energy levels. As has been argued, we only need to compute the K odd and the K even torsional levels, E odd tor and E even tor , respectively, in order to deduce the full rovibrational energy-level structure [see Eqs. (34) and (35)]. As the torsional potential in the case of H + 5 is known to be fairly weak, we can treat the torsional potential by perturbation theory.

The K even case
The K even levels are simply the energy eigenvalues of the K = 0 torsional Hamiltonian We will treatV as a small perturbation. IgnoringV , the zeroth order eigenstates can be classified by ap φ eigenvalue k ∈ Z and have energy E even tor = 2k 2 so we get the energy eigenvalues E even tor /2 = 0, 1, 1, 4, 4, 9, 9, . . .
For small torsional potentials of the form where theṼ j are the Fourier components of the potential, we have that the matrix element ofV between states withp φ eigenvalues k and k is Thus, in the first order of perturbation theory, the perturbed energies are E even tor /2 = 0, 1 − These are our K even torsional energy levels, whereṼ j diminish quickly as j increases. Note that the E even tor levels are non-degenerate, the potential has split any degeneracies in the zeroth-order levels.

The K odd case
The K odd levels are the energy eigenvalues of the K = 1 torsional Hamiltonian Again, we will treatV as a small perturbation. IgnoringV , the zeroth order eigenstates can again be classified by ap φ eigenvalue k ∈ Z and have energy E odd tor = 2 k − 1 These are our K odd torsional energy levels. They are all doubly degenerate, the potential does not cause any splitting of degeneracies in the zeroth-order levels. In fact, we know from our earlier discussion that this degeneracy of K odd torsional energy levels is an exact result.
Even if we went to higher order in perturbation theory, there would never be any splitting of these levels due toV .
The resulting rovibrational energy-level pattern, incorporating the rotational kinetic energy E rot as well as these torsional energy levels, is depicted in Fig. 2

Beyond the symmetric-top approximation
Recall that we started our analysis of H + 5 with the 1D(φ) torsion model, before making an approximation (the large R limit) which yielded a symmetric-top Hamiltonian. This was a very useful approximation, because the symmetric-top Hamiltonian exhibits the rovibrational Aharonov-Bohm effect and can be understood by analogy with motion in an electromagnetic field. It is also a physically reasonable approximation for H + 5 . We now consider what would happen if we dropped this approximation. The answer is that the quadruple degeneracies we have just found would split into 4 = 2 + 2 and the double degeneracies would split into 2 = 1 + 1.
One way of thinking about these splittings is in terms of the symmetries of the symmetrictop rovibrational Hamiltonian in Eq. (22). Firstly, we have rotations about the symmetry axis of the symmetric top which act on the wavefunction as follows: for θ ∈ [0, 2π]. Secondly, we have feasible permutations of the protons, generated by the double transposition S = (24) (31) and the transposition T π = (21), which act on the wavefunction as follows:T Finally, we have spatial inversion which is carried out by the parity operator: One can check thatR θ andŜ generate a group R θ ,Ŝ O (2), whileR − π 2T π andPŜ generate a group R − π 2T π ,PŜ D 8 (D 8 is the dihedral group of order 8). These two groups commute with each other, and together generate the symmetry group Here, the Z 2 is generated by a common central elementR π = R − π 2T π 2 of order 2.
The irreps of G symm−top have dimensions 1, 2, and 4. In fact, this is another way of understanding why the rovibrational levels we have found have degeneracies 1, 2, and 4they must correspond to irreps of G symm−top .
If we drop the symmetric-top approximation, then the only symmetries we expect to have left over are those corresponding to spatial inversion, together with feasible permutations of the identical nuclei. So we only have the subgroup of G symm−top which is generated byŜ,T π andP . This subgroup is where the Z 2 factor is generated by parity while the D 8 factor is generated by the feasible permutations. Note that this group can be identified with G 16 in Bunker's notation. 14 Clearly G sits in G symm−top as a subgroup: The idea is that while, to a good approximation, H + 5 is a symmetric top and has symmetry group G symm−top , in reality this symmetry is actually broken down to the subgroup G. To illustrate this point further, in Table 1 we reproduce numerical results from Ref. 28 for the full 1D(φ) torsion model, compared with our perturbation theory analysis of the symmetric-top (Aharonov-Bohm) approximation. For this comparison, we need parameter values. In the units we have adopted so far in this section, for which m H = r = = 1, we take 1 2R 2 = 0.0666, while the torsional potential Fourier coefficients are given byṼ 1 = 0.376 andṼ 2 = −0.016, which corresponds to a torsional barrier of 4Ṽ 1 ≈ 1.5. For the comparisons in Table 1, where energies are given in cm −1 , our energy unit corresponds to 53.34 cm −1 .
Note that K is no longer a good quantum number in the 1D(φ) model, but because the symmetric-top symmetry is only slightly broken, it is still useful to label states by |K|, as seen in Table 1. The small splittings 4 = 2 + 2 and 2 = 1 + 1 due to the breaking of the symmetric-top symmetry are evident in Table 1. For example, looking at the J = 3 column, we see that the first row corresponds to a non-degenerate K = 0 level at 42.62 cm −1 in the symmetric top approximation. The next four rows correspond to a quadruply degenerate |K| = 1 level at 92.41 cm −1 in the symmetric top approximation, which is split as 4 = 2 + 2

Conclusions
By revisiting the problem of coupling overall rotations with an internal vibrational motion for symmetric-top molecules, we demonstrated that this coupling can be understood by analogy with the famous Aharonov-Bohm effect, with the rotational motion influencing the vibrational motion in a way anologous to the influence of a magnetic solenoid on a charged particle. The effective electric charge is identified with K, the component of the angular momentum along the symmetry axis of the molecule, while the effective magnetic flux carried by the solenoid is proportional to α, the rovibrational coupling strength. In particular, the quantum energy levels associated with the vibrational motion are affected by the rovibrational coupling in a characteristic way familiar from electromagnetic theory.
As an application, we considered the low-energy rovibrational dynamics of H + 5 , which are known from previous studies to be well-described by a model coupling overall rotations to a single torsional motion. We showed that, since H + 5 is approximately a symmetric top, the rovibrational Aharonov-Bohm effect governs the low-energy dynamics of the ion to a good approximation. Moreover, in the case of H + 5 , the value of the effective magnetic flux carried by the solenoid within this analogy is such that the torsional energy levels only depend on whether K is even or odd, and that all torsional energy levels for K odd are doubly degenerate. These latter effects have analogues in the well-known molecular Aharonov-Bohm effect, in which the presence of a conical intersection of Born-Oppenheimer potential energy surfaces influences the nuclear dynamics through the Mead-Truhlar-Berry vector potential.
The electromagnetic analogy gives us a new way of understanding the rovibrational level structure of symmetric-top molecules, as we have vividly demonstrated for H + 5 . This gain in understanding is largely due to the fact that changes of embedding of molecule-fixed axes play a transparent role in the electromagnetic analogy, corresponding to gauge transformations. This is to be contrasted with the usual decomposition of the nuclear Hamiltonian into vibrational, rotational and rovibrational terms, in which the behavior under changes of embedding is somewhat obscured.
A natural question is whether these ideas can be extended to the case of asymmetric-top molecules. This can be done, but one has to replace the effective magnetic field which couples to the vibrational dynamics with a so-called non-Abelian gauge field. This perspective on general rotation-vibration coupling is a well-developed subject, 31-34 which in a sense reduces to the electromagnetic analogy introduced here in the special case of a symmetric top. An attractive feature of the electromagnetic analogy is that the mathematical ideas involved are already familiar to chemists in the context of the molecular Aharonov-Bohm effect. The non-Abelian fields which arise in the more general case are not as familiar, although we point out that they do in fact have analogues in the chemistry literature. For example, in molecular systems with strong spin-orbit coupling one may be interested in the intersection of two potential energy surfaces, each surface corresponding to a doubly (Kramers) degenerate electronic state. This picture naturally leads to the introduction of a non-Abelian gauge field which couples to the nuclear dynamics, as explored in Refs. 8, 40 and 41. We expect that the electromagnetic analogy, and its non-Abelian generalization, will continue to yield new insights in nuclear dynamics.