PCCP
View Article Online PAPER View Journal | View Issue
A semi-classical approach to the calculation of highly excited rotational energies for Cite this: Phys. Chem. Chem. Phys., 2017, 19,1847 asymmetric-top molecules
Hanno Schmiedt,a Stephan Schlemmer,a Sergey N. Yurchenko,b Andrey Yachmenevbc and Per Jensen*d
We report a new semi-classical method to compute highly excited rotational energy levels of an asymmetric-top molecule. The method forgoes the idea of a full quantum mechanical treatment of the ro-vibrational motion of the molecule. Instead, it employs a semi-classical Green’s function approach to describe the rotational motion, while retaining a quantum mechanical description of the vibrations. Similar approaches have existed for some time, but the method proposed here has two novel features. First, inspired by the path integral method, periodic orbits in the phase space and tunneling paths are naturally obtained by means of molecular symmetry analysis. Second, the rigorous variational method is Creative Commons Attribution 3.0 Unported Licence. employed for the first time to describe the molecular vibrations. In addition, we present a new robust approach to generating rotational energy surfaces for vibrationally excited states; this is done in a fully quantum-mechanical, variational manner. The semi-classical approach of the present work is applied to Received 12th August 2016, calculating the energies of very highly excited rotational states and it reduces dramatically the Accepted 8th December 2016 computing time as well as the storage and memory requirements when compared to the fullly
DOI: 10.1039/c6cp05589c quantum-mechanical variational approach. Test calculations for excited states of SO2 yield semi-classical energies in very good agreement with the available experimental data and the results of fully www.rsc.org/pccp quantum-mechanical calculations. This article is licensed under a
1 Introduction example, to limitations of the temperatures that can be reached in practice, only moderate rotational excitation can be achieved Open Access Article. Published on 09 December 2016. Downloaded 9/30/2021 3:51:21 AM. The majority of theoretical and experimental studies of molecular in this manner. With the advent of the molecular optical properties in highly excited states have focused on probing the centrifuge,1,11 generated by two powerful counter-rotating ultra- vibrational and electronic contributions to the molecular energy. fast, oppositely chirped laser pulses, it became possible to Much less is known about molecules at very high rotational prepare and control molecular ensembles at extremely high excitation. Several recent theoretical studies have demonstrated rotational excitation, even at low temperatures. The centrifuge that ultrafast rotating molecules open a new avenue for investigating has been successfully implemented to study the infrared spectrum 1,2 3–6 12,13 chemical reactivity, low-temperature collisional dynamics, of CO2 in rotational states of very high energy ( J E 220). rotational cooling and trapping,7,8 and other unusual phenomena The quantum-mechanical calculation of the rovibronic energies occurring for molecules with rotational energies comparable in of a polyatomic molecule within the Born–Oppenheimer approxi- size to their chemical bond strengths. mation (see, for example, ref. 14 and references therein) normally Until recently, excited rotational states of molecules were mostly involves the choice of a zero-order model, that is, a simplified populated in experiments by simply increasing the temperature description of the molecule leading to a Hamiltonian whose or using microwave ladder excitation techniques.9,10 Owing, for eigenvalues and eigenstates are known or easily obtainable. Historically, the zero-order model of choice has been the a I. Physikalisches Institut, Universita¨tzuKo¨ln, Zu¨lpicher Straße 77, 50937 Ko¨ln, rigidly-rotating/harmonically-vibrating molecule (see, for example, Germany ref. 15 and references therein) whose rotation is modelled as that b Department of Physics and Astronomy, University College London, Gower Street, of a rigid rotor, and whose vibration is modelled as that of a WC1E 6BT London, UK collection of uncoupled harmonic oscillators. Rotation and c Center for Free-Electron Laser Science, DESY, Notkestrasse 85, 22607 Hamburg, Germany vibration are considered as entirely independent, and so the d Theoretische Chemie, Bergische Universita¨t Wuppertal, Gaußstr. 20, eigenfunctions of the corresponding zero-order Hamiltonian 42119 Wuppertal, Germany. E-mail: [email protected] are products of well-known rigid-rotor eigenfunctions and
This journal is © the Owner Societies 2017 Phys. Chem. Chem. Phys., 2017, 19, 1847--1856 | 1847 View Article Online
Paper PCCP
harmonic-oscillator eigenfunctions.16 These zero-order eigen- capacity and, so far, they have been rarely applied to the functions are now taken as basis functions in the calculation of computation of energies and wavefunctions for highly excited the eigenstates of the actual molecule. That is, the eigen- states of molecules with more than three nuclei. However, a
functions for the Hamiltonian of the actual molecule are taken recent example of such a calculation (for H2O2) is given in as a (truncated) linear combination of the basis functions. This ref. 29. By means of the symmetry arguments mentioned ‘Ansatz’ converts the Schro¨dinger equation into a matrix above,16 it can be rigorously shown that the Hamiltonian equation: each molecular energy is given as an eigenvalue of matrix of an isolated molecule is block diagonal in the angular a matrix representing the Hamiltonian in the chosen basis momentum quantum number J (see, for example, ref. 16) and set,16 and the corresponding eigenvector contains the expansion the resulting matrix blocks can be diagonalized independently. coefficients defining the molecular wavefunction in terms of This fact can only be exploited with a basis set of eigenfunctions the basis functions. In the early days of molecular theory, of the angular-momentum operator J2, and the well-known perturbation theory represented the only possibility of diagonalizing symmetric-top-eigenfunctions16 | J,k,mi are normally used. the Hamiltonian matrix. The perturbations to the zero-order model The Hamiltonian matrix blocks are diagonal in the projection are anharmonicity and rotation-vibration coupling. The ensuing quantum number m (= J, J +1,..., J 1, J) and the matrix theoretical results, developed over many years, constitute an elements do not depend on m,16 so one can do a single extensive toolbox for the analysis of molecular spectra (again, calculation for, say, m = 0 for a given value of J. However, as see ref. 15 and 16 and references therein). The perturbation- J increases we must take into account increasingly more values theory approach runs into problem in the presence of so-called of k (= J, J +1,..., J 1, J); there are 2J + 1 values of k for resonances, when the energy separations between two or more J given. Hence, the dimension of a Hamiltonian-matrix block basis states are comparable in size to the matrix elements increases proportionally to 2J + 1 and so, for increasing J,the coupling these basis states. With the advent of high-capacity matrix dimension will eventually reach a limit at which diagonali- computers, it became possible to avoid the approximation zation is no longer feasible with the available computers.
Creative Commons Attribution 3.0 Unported Licence. inherent in the perturbation theory by simply diagonalizing With the aim of avoiding the ‘‘high-J-catastrophe’’ just described, the truncated Hamiltonian matrix numerically. This is known we explore in the present work an alternative approach to the as a variational calculation. Solution of the quantum mechanical calculation of energies for highly excited rotational states of problem by means of variational methods provides – compared molecules. It is inspired by the fact that at high rotational to the results of perturbation theory calculations – a superior excitation, as the rotational angular momentum attains ‘macro- description of molecular energies, in particular for molecules scopic’ values vastly larger than h = h/2p (where h is Planck’s in highly excited states. The variational methods have been constant), the behavior of the molecule changes from being developed during the last decade so as to provide theoretical quantum mechanical to being classical.30 This is known as the data of a quality that is, at least, close to satisfying the exacting correspondence principle. Correspondence-principle-inspired, This article is licensed under a demands of high-resolution spectroscopy.17–21 They nowadays theoretical descriptions of molecular rotation (sometimes called make use of zero-order models more sophisticated and physically semi-classical theory) were developed already in the 1970s and satisfactory than the rigidly-rotating/harmonically-vibrating 80s31–33 but, with a few exceptions,33–38 they have faded into molecule. As already mentioned, the matrices to be diagonalized oblivion. In a semi-classical approach, the rotation is described Open Access Article. Published on 09 December 2016. Downloaded 9/30/2021 3:51:21 AM. in variational calculations are truncated: the exact solution is only in terms of classical periodic trajectories on a so-called rotational obtained in the limit of an infinitely large basis set which, in energy surface (RES). The RES is a function of classical angular
practice, we must approximate by a finite number of basis momentum components ( Jx,Jy,Jz) measured along molecule-fixed functions. This finite number must be chosen large enough that axes xyz (see, for example, ref. 16). The function values are the computed eigenvalues are converged so that they do not obtained as molecular rotation-vibration energy values for the change appreciably when more basis functions are included in vibrational state under study and for given fixed values of 39–41 the calculation. As described at length in ref. 16, symmetry ( Jx,Jy,Jz). RES theory has also been developed to take into arguments can be used to facilitate a variational calculation. account the effect of tunneling between two equivalent rotational Nevertheless, the principal impediment to variational calculations, trajectories by including classically forbidden complex periodic which is most serious for highly excited states, is the substantial orbits in phase space.36 One of the purposes of the present work is to size of the basis set required to obtain converged eigenenergies review and generalize previous semi-classical RES theory by couching and eigenfunctions, and the ensuing sizes of the Hamiltonian it in the path integral formalism.42 Another – which is described for matrices whose diagonalization is often very costly in terms of thefirsttimehere–isthegenerationoftheRESbyinfullyquantum computation time and storage capacity. Several methods have mechanical calculations with the TROVE program.19,21 been developed to tackle the problem of large basis sets by A central concept of the path integral formalism is the devising appropriate ro-vibrational coordinates,21 employing quantum mechanical propagator which, in its original derivation, non-direct product contraction techniques21–23 and molecular gives the probability of finding a particle, whose position was 19 17,24 symmetry, using time-dependent basis functions, or carrying given by the coordinate value qi at the time t = 0, with the position 25 out the matrix diagonalization with low-storage iterative Lanczos coordinate value qf at a later time t. The propagator can be or filter26–28 techniques. Notwithstanding these efforts, variational expressed as an integral, which conveniently contains the methods remain extremely expensive in terms of computer Lagrangian (or the Hamiltonian) of the underlying physical
1848 | Phys. Chem. Chem. Phys., 2017, 19, 1847--1856 This journal is © the Owner Societies 2017 View Article Online
PCCP Paper
system. In very general terms, one can write the propagator makes the method of the present work a versatile tool for the
K(qi,q f,t) as: understanding of ro-vibrational molecular dynamics at high ð ð rotational excitation. This is one of the chief points of the i t KqðÞ¼; q ; t Dq exp dt 0Lðq; q_Þ ; (1) present work: our calculations directly connect the (quantized) i f qðtÞ¼q f h 0 qð0Þ¼qi rotational energies of a single molecule at large rotational speed to particular, unambiguously defined paths on the vibrationally : 2 where p is the momentum of the particle and L = mq V(q) averaged RES. This induces a clear understanding of this represents the Lagrangian with the potential energy function V(q). high-speed motion, in contrast to pure quantum studies, where N=2NQ 1 Nm 30 the eigenvalues of a Hamiltonian matrix provide energies The functional measure is defined as Dq ¼ lim dqn. N!1 2pit n¼1 more accurate than those obtained semi-classically, but The physical interpretation is rather intuitive: the probability of the where the motions associated with the quantum-mechanical
particle getting from qi to qf during the time t is given as an integral energies can only be understood (and that only to a certain
over all possible paths from qi to qf, where each individual path is extent) by complicated visualizations of the wavefunctions. The weighted by its action, equal to the time integral over the appropriate semi-classical approach allows an intuitive and definitive inter- Lagrangian function along the path in question. The action attains pretation of the molecular states in terms of quasi-classical its maximum value for the classical path but, in contrast to usual motions. classical mechanics, the path integral takes into account all other, For the proof-of-principle numerical tests of the present
less probable paths which makes it a suitable tool for describing work we have chosen the SO2 molecule since its energy level quantum mechanical processes. The strategy of weighting the paths pattern has been characterized to high angular momentum by their action is the starting point of the semi-classical description. values ( J r 80), both experimentally48 and in fully quantum- Since the path integral can also be used in a perturbative context, mechanical, variational calculations.49 For several vibrational
one can allow only for ‘‘small’’ fluctuations around the classical path states of SO2, we have calculated semi-classical rotational
Creative Commons Attribution 3.0 Unported Licence. and hence derive a semi-classical approximation of the propagator. energies which show very good accuracy when compared to It is ‘‘the ideal starting point for semiclassical calculations’’ (ref. 43, experimental data or with the results of full quantum mechanical p. 262), but is nowadays also used as the fundamental basis of calculations. By comparisons with experiment we must take into
quantum field theory, which is highly important in many fields of accountthatonlytheA1 and A2 states of SO2 – with symmetry 16 single particle and many-body physics and chemistry. In the group C2v(M) – are allowed by nuclear spin statistics, while B1 16 approach of the present work, we use the density of states – a and B2 rovibrational states are forbidden. Our calculations do physical observable – as the starting point for the discussion of the not as yet consider spin statistics. semi-classical approximation for the energy levels of the rotating The paper is organized as follows. In Section 2 we formulate molecule.ThedensityofstatesiscalculatedfromtheFourier the general semi-classical treatment of molecular rotational This article is licensed under a transform of the propagator. As the rigorous derivation of the density motion using the path integral approach. The calculation of the of states is beyond the scope of the present work, we refer to RES for a general molecule by means of TROVE is presented in textbooks such as ref. 30 and 43. The density of states exhibits poles Section 3. In Section 4 we report semi-classical calculations of
at the quantized energy levels and these poles can be calculated the rotational energies for the SO2 molecule and assess the Open Access Article. Published on 09 December 2016. Downloaded 9/30/2021 3:51:21 AM. using the semi-classical RES by identifying certain periodic orbits on accuracy of the results. Section 5 finally provides a summary it. These paths are selected by symmetry arguments and can be and a discussion of the perspectives for future studies. calculated by considering the operations of the molecular symmetry group,16 whichcanbecorrelatedwithmotionsontheRES. The resulting, new method describes the effect of tunneling 2 The semi-classical approach in a more rigorous and natural way involving a molecular At large angular momentum, we expect the rotational dynamics symmetry group16 analysis; this yields results equivalent to to be governed by the classical Euler–Lagrange equations of the original ones from ref. 36. Another novel aspect of the motion. Accordingly, a path integral approach – being a versatile proposed method is that the underlying RES is obtained quantum tool for studying classical mechanics – is well suited for solving mechanically by means of variational theory as implemented in these equations. Assuming the actual motion to take the form of computer program TROVE.19,21 This is probably the most small fluctuations around the classical path, we take the density important progress made in the present work. Indeed, it is of states r(E) to be given by the semi-classical expression long not the first time that quantum-mechanically calculated RES’s 42 known as ‘‘Gutzwillers trace formula’’ have been subjected to classical analysis,44,45 but we make a 1 X further development here in that we describe the tunneling iSaðEÞ= h rðEÞ r ðEÞþ Re TaFae : (2) effects just mentioned. In earlier works, effective Watsonian- p h a type Hamiltonians or (semi-)classical vibrational models have been used to determine the RES.33,46,47 The combination of the Here, we sum over periodic orbits a of the classical system
high accuracy first-principles quantum mechanical treatment weighted by the respective action Sa = padqa, the period Ta and a
of the (excited) vibrations by TROVE, and the semi-classical stability amplitude Fa, which here includes the Maslov index. treatment of the rotation (including the tunneling effects), The quantity r (E) is the smooth part of the density of states,
This journal is © the Owner Societies 2017 Phys. Chem. Chem. Phys., 2017, 19, 1847--1856 | 1849 View Article Online
Paper PCCP
calculated as the purely classical contribution. For a detailed possible combinations of, e.g., the number of real segments mathematical description of the semi-classical approach to between two transmissions or over the number of transmissions obtaining the density of states in the path-integral formalism, in total, is initially rather involved but, fortunately, it can be we refer the reader to ref. 30, 50 and 51. reduced to a geometric series, and a relatively simple expression In rotationally resolved molecular-spectroscopy studies, the is obtained: actual energy levels are of principal interest. In the approach of noÀÁ the present work, they can be determined at the poles of the dm iS y 1 iS roscðÞ¼Em Tr e DmðoÞ 1m tDmðtÞ e DmðoÞ oscillatory part of r(E). Furthermore, the energy levels are jGj labeled by their symmetry species (irreducible representations) (5) in the appropriate symmetry groups16 of the molecular Hamiltonian. With the symmetry labels we can obtain, for example, selection rules where o and t parametrize the real and tunneling symmetry
for state-to-state transitions. Projecting the density of states in eqn (2) elements, respectively, and the Dm are the representation matrices
onto the space of states belonging to one particular irreducible constituting Gm (1m is the unit element of Gm). The poles, i.e. the representation Gm, we obtain the symmetry reduced oscillatory part energy levels, are now found by the determinant equation: 52 of the density of states: ÀÁ det e iS Dy o tD t 0: 1 d X mð Þ 1m mð Þ ¼ (6) m iSa= h rosc ðÞEm Re TaFawmðÞga e ; (3) p h jGj a Here, o and t describe the symmetry operations parameterizing
where dm is the dimension of Gm, wm the character of the related the real and tunneling paths, respectively. Equipped with this equation and the molecular symmetry group, it becomes possible matrixP representation of the symmetry group element g and 2 jjG ¼ dm denotes the number of elements in the symmetry to find the quantized energies of any well-defined rotational energy m surface. In the following section, we will apply this strategy to group G.Furthermore,theorbitstobesummedoverarenow Creative Commons Attribution 3.0 Unported Licence. the test molecule SO2 whose energies are well known, both periodic in the symmetry reduced phase space and their initial experimentally and theoretically, up to large values of the and final points are related by a symmetry operation ga,whichisthe angular momentum quantum number J. unit element, ensuring the periodicity of the orbit. This unit element can be composed of numerous symmetry elements. To include also quantum fluctuations superimposed on the 3 The rotational energy surface classical motion, we allow for so-called beam splitting, where the paths get reflected or transmitted at their classical turning For a single vibrational state, the semi-classical approach points in phase space.36,53 The transmitted – and hence classically initially requires a RES35,36 obtained, as mentioned above, by 54 This article is licensed under a forbidden – paths are described in phase space by a Wick-rotation replacing the molecule-fixed components of the angular
to imaginary time so that the momentum becomes complex. momentum operator ˆJ =(ˆJx,ˆJy,ˆJz) by their classical, time-
Consequently, the related action is imaginary and produces an dependent analogons J =(Jx, Jy, Jz). This defines a Hamiltonian
exponential damping factor in the density of states (cf. the function H depending on ( Jx, Jy, Jz) and we make allowance Open Access Article. Published on 09 December 2016. Downloaded 9/30/2021 3:51:21 AM. discussion of instantons in double-well potentials in ref. 30). for quantum mechanics in that we fix the length of J to be 2 2 We introduce the well-known WKB reflection and transmission | J| = J( J +1) hH . This function is then used to calculate the 46,55 amplitudes r and t for the quantum fluctuations and in action SðEÞ¼ pdq of the closed orbits, where p is simply the terms of these (see Appendix for their definitions), the oscillatory solution of H(p,q)=E. The trajectories themselves are found on part of r(Em)isformallygivenby: the RES and hence it is the fundamental starting point for the d X pffiffiffiffiffiffi approach described in Section 2 when applied to the rotational m iSa roscðÞ Em Re AaTawmðÞga e ; (4) problem. jGj a The RES has been successfully used to explain qualitatively !0 1 ðaÞ ðaÞ the emergence of rotational energy clusters for large J quantum nQr nQt 2 @ 2A 33 where Aa ¼ jjri rj describes the product of numbers. In order to provide proof of principle, in the present i¼1 j¼1 work we focus on a simple molecule with no rotational-energy (a) (a) nr reflections and nt transmission probabilities, characterized clustering. The absence of energy cluster formation requires
by ri and tj respectively. An appropriate description of the paths, that the RES does not change topologically as J increases. We
combined by real (i.e., classically allowed) and tunneling (i.e., have found that the sulfur dioxide molecule SO2 satisfies this classically forbidden) segments is made by carefully choosing the requirement and, in addition, its rotation-vibration energy level symmetry group elements describing the different segments and pattern is well characterized up to high J values.48,49,56 For these combining them in a sequence corresponding to the sequence of two reasons, we have chosen to validate our approach by
factors in Aa. Thus, each path is characterized by a sequence of computing the rotational energies of SO2. the type ra 1teiaS,wherea is the number of real segments between To explain how we calculate the RES, we consider now an initially the one transmission and the end of the path (cf. eqn (11) in general ro-vibrational Hamiltonian for an N-atomic molecule ref. 46) and the action S is identical for all a. The sum over all defined in terms of the 3N 6 generalized vibrational coordinates
1850 | Phys. Chem. Chem. Phys., 2017, 19, 1847--1856 This journal is © the Owner Societies 2017 View Article Online
PCCP Paper
xk, the conjugate momentum operators pˆk = ih q/qxk, and three the RESs for several vibrational states and several values of
angular momentum operators ˆJx, ˆJy, ˆJz J may be obtained on a grid of (y, f) coordinates by solving the 1 X vibration-only eigenvalue problem for each grid point. One may H^ ¼ H^ þ G ðxÞJ^ J^ rv vib 2 ab a b think that the task of finding the RES by repeatedly solving the a;b¼x;y;z eigenvalue problem for thousands of points is quite expensive "#(7) X 3XN 6 and unlikely to be very competitive with the full variational 1 ^ þ GakðxÞp^ p^ GkaðxÞ Ja: treatment. This is partly true for low values of J but, as 2 k k a¼x;y;z k¼1 J increases the steep (2J +1)3 computational scaling of the full variational treatment is quickly outstripped by the RES calculation, Here, Hˆvib denotes the purely vibrational Hamiltonian and includes the vibrational kinetic energy operator, the Born–Oppenheimer which involves the diagonalization of matrices with dimensions potential energy function and the pseudo-potential terms (see, for independent of J. example, ref. 19 and 21). The remaining two terms (involving In the present work we employ two methods for solving the vibrational problem: the Watson-type effective-Hamiltonian the ˆJa operators) describe the rotational kinetic energy and the 57 ro-vibrational (Coriolis) coupling, respectively. In eqn (7), the method, and the rigorous variational method, as implemented in the computer program TROVE.19,21 The purpose of considering quantities Gab(x) and Gak(x) are functions of the vibrational the former method is somewhat academic; the values of the coordinates x x1,...x3N 6. Depending on the definition of x effective-Hamiltonian parameters are normally determined by and the choice of the molecule-fixed Cartesian axes, Gab(x) and 20 least-squares fitting to experimental spectroscopic data, and the Gak(x) can be evaluated analytically or numerically on a grid or as a power series expansions around a suitable reference structure.19,21 amount of such data available is typically limited to rotational As mentioned above, the rotational energy surface, required states at moderate excitation in a few vibrational states. The for determining the paths in the semi-classical approach discussed extrapolation to high rotational excitation is problematic. On the contrary, the variational approach TROVE is general and can here, is now defined by replacing in eqn (7) the operators (ˆJx, ˆJy, ˆJz)
Creative Commons Attribution 3.0 Unported Licence. be applied to molecules of an arbitrary structure. It relies on the by their classical counterparts ( Jx, Jy, Jz), here defined as qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi ab initio or spectroscopically refined potential energy surfaces 2 and makes no approximative assumptions about the kinetic Jx ¼ h JðJ þ 1Þ p cos f; (8) energy operator. The effects of anharmonicity and Coriolis qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi coupling are intrinsically included in a TROVE-generated RES. 2 Jy ¼ h JðJ þ 1Þ p sin f; (9) For more details about the TROVE approach, the interested reader is referred to ref. 19 and 21.
Jz = hp, (10)
This article is licensed under a where, for convenience, we have introduced 4 The semi-classical energies pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi p ¼ð iÞs JðJ þ 1Þ cos y: of the distortable rotor
Here, y and f are usual spherical coordinates defining the To our knowledge, Robbins et al.36,46 have provided, more than Open Access Article. Published on 09 December 2016. Downloaded 9/30/2021 3:51:21 AM. direction of J in the molecule-fixed axis system xyz, and s =0 twenty-five years ago, the very first example of the semi-classical,
or 1. A s-value of 0, with Jz being real, corresponds to a real, path-integral method discussed here being applied to a molecule. classically allowed path, whereas a value of 1, with an imaginary They treated the rotational dynamics of SF6, described by the 58 value of Jz, corresponds to a classically forbidden tunneling (effective) Hecht Hamiltonian with octahedral symmetry. Their path (see Section 2). Following ref. 33, also in the s =1 results showed good agreement with energy levels determined in 2 2 2 2 tunneling case we calculate ( Jx + Jy )/ h as J( J +1) p quantum mechanical calculations. We will first investigate how [eqn (8) and (9)], in spite of the fact that in this case, this the semi-classical method performs when applied with a RES quantity is larger than J( J + 1). obtained from a Watson-type effective Hamiltonian, and then - 19,21 We make only the replacement (ˆJx, ˆJy, ˆJz) ( Jx, Jy, Jz)inHˆ rv; apply it with a TROVE-generated RES. Our results suggest that
we do not change Hˆ vib. Thus, we obtain a Hamiltonian that it is generally useful in conjunction with any Hamiltonian retains in full the quantum mechanical description of the adequate for describing the molecular energy levels. It does vibrational motion while it depends parametrically on the two not require the Hamiltonian to depend on the angular momentum rotational variables (y, f) and on the angular momentum operators in any particular way, in contrast to the approach of
quantum number J. To our knowledge, this is the first time ref. 36 which used the ‘‘Hecht Hamiltonian’’ for SF6.Itisnot that RESs have been generated for different vibrational states in limited to any power of the perturbation expansion used to obtain a full quantum solution. the Watson-type effective Hamiltonian and it can be applied to a We use TROVE to solve the ro-vibrational Schro¨dinger vibrationally averaged RES as done in the TROVE approach. As
equation at fixed values of the parameters y A [0,p], f A [0,2p], already mentioned, our test case is the sulfur dioxide molecule SO2 and J. Each of the resulting vibrational eigenvalues E1,...En of (see Fig. 1), which is not expected to exhibit rotational-energy the Hamiltonian provides a point on the RES for the vibrational cluster formation, and for which there are extensive data available
state in question, i.e. RESi(y,f)=Ei. Following this procedure, for comparison.
This journal is © the Owner Societies 2017 Phys. Chem. Chem. Phys., 2017, 19, 1847--1856 | 1851 View Article Online
Paper PCCP
Table 1 Parameterizations of the real and tunneling paths of SO2 obtained for the quantization axes x and z and for the four possible
irreducible representations of C2v(M) (see text). The quantization conditions obtained from eqn (12) require the listed quantities to vanish for an integral value of n. Here, d is the phase shift at a classical turning point and S˜ represents the tunneling action. For exact definitions see appendix and ref. 46
Jx-Quantization Jz-Quantization p p o (12)* 8 Rx E* 8 Rz Fig. 1 The sulfur dioxide molecule SO at equilibrium with the molecule- 8 p 8 p 2 t (12) Ry (12) Ry fixed axis system xyz attached. The molecule is in the xy plane. J 1 S˜ J 1 S˜ A1 d + S ( 1) tan e 2np d + S ( 1) tan e 2np J 1 S˜ J 1 S˜ A2 d + S ( 1) tan e (2n +1)p d + S ( 1) tan e (2n +1)p J 1 S˜ J 1 S˜ In the absence of rotational energy clusters, the RES for an B1 d + S +( 1) tan e 2np d + S +( 1) tan e (2n +1)p B d + S +( 1) J tan 1 eS˜ (2n +1)p d + S +( 1) J tan 1 eS˜ 2np asymmetric rotor generally has two extrema corresponding to 2 two main axes of rotation. These axes give rise to the largest and smallest moment of inertia, respectively. A radial plot of the calculated in terms of the angular momentum along this path. In
RES shows the extrema explicitly and such plots, compared particular, for the paths closing around the Jx axis (called Jx with constant-energy spheres, are used to determine the classical quantization in Table 1), the Jz axis is chosen to define the 59 solutions of the Euler–Lagrange equations of motion. In parti- momentum p,whereasintheJz quantization case, Jx = p. cular, the intersection of the RES with a constant-energy sphere In order to find the poles of the density of states, the represents a classical path on which the angular momentum and parameterization of the paths in terms of the real and the tunneling the energy are both conserved. This corresponds to the physically segments is crucial. It depends on the chosen projection axis and is intuitive picture of a rotating body with a main axis of rotation done using two different rotation operators (cf. Table 1). The two
Creative Commons Attribution 3.0 Unported Licence. that precesses around the angular momentum vector (whose operations o and t determine the path and hence they define the direction and length are conserved, i.e., time-independent). In periodic orbits. Periodicity is ensured by the use of a combination the semi-classical approach, we project the RES onto the two of symmetry elements, so that their successive application, as main axes of rotation (cf. Fig. 2). Inspecting the classical paths described in Section 2, defines a path on the RES starting and corresponding to these two projections, we notice that, depending endinginthesamepoint.ThisisshowninFig.2,whereclassical on the energy, there is a single axis around which the path is closed. and tunneling paths are displayed on the (projected) RES. As The projection, in itself, defines a coordinate and a conjugate stated above, only for one choice of the projection these closed momentum. This conjugate momentum is paramount for carrying paths exist and hence, the two choices of (o,t) are mutually out the semi-classical analysis; the action for the path in question is exclusive. This article is licensed under a If a symmetry group is specified for the molecule under study, we follow the strategy of Robbins et al.46 and determine the poles of the density of states [eqn (4)] by solving eqn (6). For 16
Open Access Article. Published on 09 December 2016. Downloaded 9/30/2021 3:51:21 AM. SO2, the molecular symmetry group is the simple, four-element group
C2v(M)={E,(12), E*, (12)*} (11)
where16 E is the unit operation, (12) is the interchange of the two O nuclei, E* is the inversion operation and (12)* = (12)E*.
The C2v(M) group has four irreducible representations A1, A2,
B1 and B2 which are all non-degenerate (see Table A-5 of ref. 16). Thus, the quantized energy levels E of the asymmetric rotor are straightforwardly determined as the solutions of