Abstract

Rotational spectroscopy can provide insights of unparalleled precision with respect to the wavefunctions of molecular systems that have relevance in fields as diverse as astronomy and biology. In this paper, we demonstrate how asymmetric molecular pure rotational spectra may be analyzed “pictorially” and with simple formulae. It is shown that the interpretation of such spectra relies heavily upon pattern recognition. The presentation of some common spectral line positions in near-prolate asymmetric rotational spectra provides a means by which spectral assignment, and approximate rotational constant determination, may be usefully explored. To aid in this endeavor we have created a supporting, free, web page and mobile web page.

1. Introduction

Microwave spectroscopy is capable of providing unique insights into the electronic structure and potential energy surfaces of molecules. It has found considerable use in the diverse fields of physics and biology, providing an invaluable method to examine the conformational landscapes of both small and large molecules. Microwave spectroscopy is a mature discipline and several recent review articles are available on the subject [13]. Yet recent technological advances [4], together with lower instrument costs mean that microwave spectroscopy is more accessible to researchers than it has ever been before [510]. The recent 2012 Microwave Newsletter [11] indicates that there are over 50 active research laboratories in locations all over the globe, with at least 400 researchers studying over 800 molecular systems. Furthermore, technology developments now allow broadband rotational spectra to be collected in a few hours, or less, meaning that undergraduate students can perform experiments in a timely fashion, that is, an undergraduate physical chemistry lab. This being so, the broadband nature of the data collected means that researchers are still faced with a complex spectrum that may not be as easily interpreted as many classroom NMR spectra or low resolution IR spectra.

The rigid rotor, as a model for a rotating linear molecule, is a subject covered in elementary quantum mechanics classes. It is unusual for classes to consider molecular asymmetry, because this is often viewed as being too complex a problem. It is the goal of this work to make this important task easier.

All types of spectroscopy require that the observed transitions, or peaks, be identified. For many types of spectroscopy spectral peak identification is associative. For example associating a specific bond stretch with a peak in an IR spectrum, or associating a methyl group with a group of peaks in an NMR spectrum. However, in many high resolution spectra, the task of transition assignment requires identifying the upper and lower quantum states involved. The spectra produced from microwave frequency studies are usually associated with the quantized rotational motions of molecules, and the task of spectral peak identification involves the latter, that is, the assignment of quantum numbers for the upper and lower states.

A rotational spectrum, like most spectra, should be viewed as a puzzle. These puzzles contain patterns, often complex in nature, that once unravelled lead to easily tractable information regarding the geometry and conformation of the subject molecule. In order for spectral patterns to be unravelled, the spectroscopist needs to know the “rules” of the game, that is, the underlying theory governing the asymmetric quantum rotor. For any given molecule these rules may be arrived at from an appropriately formed Hamiltonian operator, solving for the eigenvalues, and then either deriving or simply applying the appropriate selection rules. This could be called a bottom-up approach. But there is an alternative where, first, puzzle solving skills are invoked and spectral patterns identified. After this point the appropriate molecular properties may be introduced and discussed. This may be considered a top-down approach.

To this end several, powerful, software visualization tools have been developed to help with the quantum number assignment problem in pure rotational spectroscopy [1216]. In order to compliment these tools the provision of some common spectral patterns would be useful for the beginning student.

The aim of this paper is to present some common patterns that occur in rotational spectra, together with their approximate dependence on the magnitudes of the rotational constants, so that a top down approach is more available to interested parties. Importantly, we point out that by “patterns” we are referring to transition line positions; this work does not address transition line intensities. None of the theory introduced here is new (its over 70 years old [1719]), however, the method of presentation goes beyond that found in other sources. Through this paper it is hoped that students may be more easily introduced to rotational spectra, and that rotational spectroscopy may be more fully included in the chemistry curriculum.

2. Background

In order for spectral patterns to be identified, and in order to produce a self-contained paper, three concepts need to be introduced. These are (i) molecular rotational constants, (ii) the appropriate quantum numbers to use in transition assignment, and (iii) an overview of the selection rules with particular reference to how these selection rules relate to the geometry of the molecule. These are introduced in only a very cursory manner as they are treated thoroughly by several excellent monographs [2026].

2.1. Molecular Rotational Constants

The geometry of any molecule may be reduced to three rotational constants, , , and , which are inversely proportional to the moments of inertia about the three principal axes, , , and . By convention, the -axis system is orientated such that the rotational constant is the largest in magnitude, and is smallest. Accordingly, an important skill is to be able to locate the principal axis in such a way that it runs through, or close to, as much atomic mass as possible, followed by the principal axis such that it is perpendicular to the -axis and runs through as little atomic mass as possible. For examples of this axis system Figure 1 displays three different molecules in their principal axes. The three molecules are thioxoacetaldehyde [27] (top), hexafluoroisobutene [28] (center), and 1,1,1,3,3,3-hexafluoropropane [29] (bottom).

The relationship between the magnitudes of the rotational constants, that is, the shape of the molecule, governs the gross appearance of a molecule’s rotational spectra. Ray’s asymmetry parameter [18], = ()/(), provides a quantitative measure of how far a molecule is from any symmetry relations. In the symmetric prolate limit where , then = −1, whereas in the symmetric oblate limit where , = +1. Figure 2 displays the moment of inertia ellipsoids, given by ; , for the prolate, perfect asymmetric, and oblate tops [30]. These ellipsoids can be instructive because, whereas the rotation of a diatomic molecule is easily visualized, the rotation of an asymmetric molecule is less easily imagined. The moment of inertia ellipsoids provides a visual cue, helping students visualize asymmetric molecular rotation by imagining how the “ellipsoid rolls without slipping on a plane” [31]. Furthermore, one can see how the relationships between the moments of inertia affect the rotation.

Naturally, most molecules are asymmetric and, as may be observed in Figure 2, possess moment of inertia relationships that lie somewhere between the oblate and prolate limits. For many molecules is close to −1 and they are then described as near-prolate asymmetric tops. This is a very common molecular shape and this paper focuses mainly on this type of rotational spectrum.

2.2. Rotational Quantum Numbers for Asymmetric Molecules

The rotational quantum numbers for asymmetric tops are denoted as . The rotational quantum number is familiar from the linear rigid rotor; however, the and labels (not quantum numbers) require brief explanation. Consider a symmetric top, either oblate or prolate. A second quantum number, , is required in the energy expressions for these molecules. This quantum number pertains to the quantized rotation about the molecule’s symmetry axis. The quantum number can take values of through to in integer steps. However, only a term appears in the energy expression and so, excepting the levels, all levels are doubly degenerate. Now, for a prolate top, the energy levels for a specific quantum state, , increase in energy from up to . In stark contrast to this the oblate top energy levels, , increase in energy from up to .

For the asymmetric top there is no symmetry axis. This means is not a good quantum number for an asymmetric molecule. However, each energy level for the asymmetric top may be linked to (i) a value in the prolate limit, this becomes the label and also (ii) a value in the oblate limit, which becomes the label. Essentially, we may imagine that if an asymmetric top in a state were to “morph" into a prolate top, then the quantum state would be . If, on the other hand, the asymmetric state were to “morph” into an oblate top the quantum state would be .

2.3. Selection Rules

The gross selection rule for rotational spectroscopy is that the molecule must have a dipole moment. In the case of an asymmetric molecule one must consider the components of the dipole along either the and/or and/or principal axes. A key skill requires that electronegative/electropositive elements, or groups, be identified for the subject molecule within the axis system. In this way one can tell a priori along which axis or axes a dipole component will be nonzero. Figure 1 shows three examples of molecules with either a-, b-, or c-dipole components.

A molecular rotation about the principal axis places different symmetry requirements on the relation between upper and lower wavefunctions in a rotational transition compared to a rotation about either the or axes. These requirements are manifest in the relation between the upper and lower quantum states labels. It is the parity of the labels, that is, even () or odd (), that we concern ourselves with. In short, for a molecular rotation about the -axis, along which there is a dipole component, a change in parity is only allowed for the label. This is referred to as an -type transition. For a -type transition, a change in parity is only allowed in the label. For a -type transition, a change in parity must occur in the and labels simultaneously.

In regards to the changes permitted in the rotational quantum number it may be shown that . Here means (upper state, higher energy) (lower state, lower energy). Just as is the case for vibration-rotation spectroscopy, transitions of the type are referred to as -branch transitions, and by extension is the -branch, and is the -branch. In passing we note that, in our experience, Q- and R-branch transitions are by far the most commonly observed in pure rotational spectra.

Lastly, a shorthand notation for any transition, or group of transitions, is commonly used. By example, an transition is an -type, -branch transition in which the value of does not change between the lower and upper states, and the label increases by 1 between lower and upper states.

3. A Note on Centrifugal Distortion

In this work we have assumed the rotating molecule to be rigid. Real molecules do, however, undergo distortions in their average nuclear positions as rotation occurs. This is referred to as centrifugal distortion and is accounted for in a descriptive set of molecular parameters by several centrifugal distortion constants [32]. The magnitudes of these constants are, crudely, inversely proportional to the square of the molecular mass. The centrifugal distortion constants are many orders of magnitude smaller than the rotational constants. The effect of centrifugal distortion on the observed spectra is to change the frequency of a given rigid rotational transition. These effects are generally minimal, that is, a few kHz or so for low transitions, but can become significant, that is, many MHz, for high transitions. They have been ignored in this work because at low , that is, less than, say 10, centrifugal distortion does not affect the patterns of rotational transitions observed. However, with modern microwave spectrometers, the precision of measurement of rotational transitions is such that the effects of centrifugal distortion are routinely observed and must be accounted for in the Hamiltonian operator.

4. Rotational Spectra

In the following we have selected some common transition series observed in rotational spectra for asymmetric molecules. The transitions are discussed using the notation , where the double prime indicates the lower state.

We describe them pictorially and by the provision of simple formulae. It should once again be noted that effects of centrifugal distortion have been neglected. The spectra displayed in Figures 3 through 8 were calculated exactly using the rotational constants provided in each figure using a complete Hamiltonian matrix diagonalization routine [33, 34]. The calculated spectra relate to molecules at rotational temperatures of 3 K, which is a common temperature for operating a modern microwave spectrometer, where molecules are entrained in supersonic expansions of an inert gas. The simple formulae provided have been determined through consideration of both the closed algebraic expressions for the lowest energy levels [21] and also through a posteriori examination of exact spectra. A useful compilation of empirical relations for the frequencies of rotational spectral transitions is provided in Table 7.11 in “Microwave Molecular Spectra” by Gordy and Cook [26]. However, to the best of our knowledge, many of the formulae provided here appear in print for the first time.

In regards to strategy, assignment should begin by examining the entire rotational spectrum and looking for repeating patterns. If patterns can be located they may be compared to the portions of spectra displayed in Figures 3 through 8. Once patterns have been identified the empirical relations provided below, and summarized in Table 1, allow initial approximations of the rotational constants , , and . With initial estimates of the rotational constants further patterns may be searched for and the process repeated in a bootstrap fashion.

4.1. -Branch Transitions
4.1.1. a-Types

These transitions are perhaps the easiest to locate in a broadband spectrum. The ease with which they may be located is based on their regular harmonic pattern. They occur in groups centered at approximately ()(), that is, every MHz a similar looking group is observable. Each group spans ()(). A sample of one of these groups is shown in Figure 3 for three different asymmetries. It is observed that the strongest transition in this group, the transition, moves to lower frequency, away from the center of the group, as moves further away from −1.

Within each group of transitions, and when the rotational constant is of modest magnitude compared to and , transitions involving high values of the label are observable as doublets in a characteristic pattern. The pattern contains unresolved doublets, involving the highest values, emerging from approximately ()(). The pattern progresses to higher frequencies, with doublets appearing with increased splittings as decreases in value. An example is shown in Figure 4, and the pattern is also visible in Figure 3. For low transitions, that is, , the splitting in the doublet is crudely given by . For higher transitions this estimate is approximately one order of magnitude too small.

4.1.2. -Types

(1) . The lowest frequency member of this group of transitions, that is, the transition, occurs at . Transitions in this family then progress to higher frequencies spaced by , giving a general formula of . An example of the progression is shown in Figure 5.

(2) . The lowest frequency member of this group of transitions, that is, the transition, occurs at . Transitions in this family progress to higher frequencies spaced by 2. This results in a general formula of . An example of the progression is shown in Figure 5.

(3) and . For a given these two transition types appear as characteristic doublets within the rotational spectrum. The lowest pair of transitions allowed are the and the transitions. The center frequency of the doublets in this family of transitions is given by . The splitting between the doublets is given by approximately , that is, half of , multiplied by the product of the upper and lower values. Given the dependence of the frequency on 3, these transitions often fall at higher frequencies than the previously listed transitions. For example, for the the transition falls at approximately 179 GHz. An example of the progression is shown in Figure 6.

(4) and . For a given these two transitions appear as very closely spaced doublets. For values of approaching −1, where is small, the doublets are often not resolvable. The center frequency for the doublets appearing in this family is given by . As was the case for the previous set of transitions, these doublets often fall at high frequencies, owing to their dependence on . An example of the progression is shown in Figure 6.

4.1.3. -Types

(1) . The lowest frequency member of this group of transitions, that is, the transition, occurs at . Transitions in this family then progress to higher frequencies spaced by 2, giving a general formula of . An example of the progression is shown in Figure 7. We note in passing that it is not uncommon for -type transitions and -type transitions of the type listed above to occur in similar patterns. Caution should be made to ensure that the nonzero dipole components of the molecule under study are identified prior to analysis (see Figure 1).

4.2. -Branch Transitions
4.2.1. -Types

(1) . The lowest frequency of these transitions, with =2, occurs at 3. For the transition frequencies increase according to the recursion formula . For near prolate cases is a small number and accordingly only very high transitions would fall in observable regions where they will likely be low in intensity.

(2) . For low transitions the formula works reasonably well. Transitions are found as closely spaced groups where the level changes, the does not change, and the level increases by 1. For example, the transition will be found close by (a few tens of MHz) to the transition and so on.

4.2.2. -Types

These transitions produce very distinct patterns in rotational spectra. An example is shown in Figure 8. Transitions of the type “appear," that is, when , at frequency of . The transitions then “fan" out to lower frequencies with increasing according to the recursion formula . For molecules possessing small rotational constants, that is,  GHz, the next group of transitions, for which , occur at approximately higher in frequency.

Related transitions of the type “appear," that is, when , at at frequency of . These transitions fan out to higher frequencies with increasing according to the recursion formula . Again, for molecules possessing small rotational constants, that is,  GHz, the next group of transitions, for which , occur at approximately higher in frequency.

For the higher members of these types of transitions, the pattern observed is such that increasingly higher transitions appear at increasingly lower frequencies as doublets, whereby the two transition types mentioned immediately above, begin to overlap one another.

4.2.3. -Types

In many cases the -type Q-branch transitions appear at frequencies very similar to those outlined previously for -type Q-branch transitions. However, in our experience, molecules for which -type Q-branch transitions are important in the initial stages of quantum number assignments are very rare.

5. How Well, and under What Restrictions, Do These Simple Formulae Work?

The line positions, described by the formulae above, are summarized in Table 1 and shown in Figures 3 to 8. To test the performance of these relations the following procedure was performed. A simple FORTRAN 77 program was written (an Excel spreadsheet, or MathCad program, would achieve the same purpose) to use the relationships in Table 1 to generate the transition frequencies for a total of 75 transitions of up to 10, at different values of Ray’s asymmetry parameter, from −1 to −0.8. The frequencies generated were compared to the exact frequencies for the transitions generated by a complete Hamiltonian matrix diagonalization routine using Herb Pickett’s SPCAT program [33, 34]. For each value, and for the three cases (75 transitions), (43 transitions), and (23 transitions), a root-mean-standard deviation was calculated according to where is the number of transitions. The results are shown in Figure 9. It is not at all surprising that the simple formulae provided work poorly at high and high asymmetries. If one sets a quantitative limit of desiring the predicted transitions to fall, on average, within 50 MHz of their actual location, then it is found that it is only possible to work with high transitions, that is, up to 9, when is between −0.975 and −1. If, on the other hand, it is necessary to work with an asymmetric species where is −0.8, or even lower in magnitude, then it is found that the formulae presented produce transition frequencies within 50 MHz of the actual frequencies only for transitions involving less than or equal to 3.

This may appear very unsatisfactory. However, the situation is better than it seems when one recalls, the very thesis of this paper, that it is the pattern of transitions that is more informative in regards to quantum number assignments than the absolute frequencies. Figure 10 shows this comparison. The upper portion of Figure 10 shows a rotational spectrum calculated exactly using SPCAT. The lower portion shows a stick spectrum produced from the relationships given in Table 1 and the same rotational constants. The match is sufficient that a quantum number assignment could likely be made for over 50% of the transitions shown.

6. New Media Support

For quantum number assignment problems, students are encouraged to plot stick spectra using the formulae given in Table 1. However, in order to rapidly visualize the patterns of transitions formed by the relationships given in Table 1 both a web-based visualization tool [35] and matching smartphone mobile web page [36] are being developed. Using either medium, users can enter a set of rotational constants, press “Plot Spectrum," and then see the patterns of low transitions that result from the relationships given in Table 1. Only low transitions are shown in order to maintain a low spectral density. The transitions displayed are colorcoded so that the -type transitions are shown in red, the -type transitions are in green, and the -type transitions are in yellow. A screenshot of the smartphone application is provided in Figure 11. The display format follows closely that of the excellent AABS software package for the assignment of broadband microwave spectra developed by Professor Kisiel [12, 13]. Users are encouraged to contact the authors with comments regarding improvements.

7. Conclusions

Rotational spectra may be inspected, initially, so as to locate repeating patterns of transitions. We have presented “template” spectra indicating the types of patterns that can commonly occur for the near-prolate asymmetric molecule. Location of these patterns in experimental spectra can help leading to succesful quantum number assignments and/or initial evaluations of the subject molecule’s rotational constants.

The information provided allows students to adopt an alternative approach to considering the quantum mechanical, asymmetric rigid rotor. In this method students are actively involved early, that is, in the task of spectral assignment. Web and smartphone resources have been created to assist in this endeavor. The tools provided facilitate, and prompt, in depth discussion about the under lying theory and application of microwave spectroscopy.

Acknowledgments

We are grateful to Mr. Frank DeChirico, Mr. Bryan Harrison, and Professor S. E. Novick for proof-reading this paper and their improving comments.