PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 45068.
Published online 2017 March 24. doi:  10.1038/srep45068
PMCID: PMC5364483

Simulating electric field interactions with polar molecules using spectroscopic databases

Abstract

Ro-vibrational Stark-associated phenomena of small polyatomic molecules are modelled using extensive spectroscopic data generated as part of the ExoMol project. The external field Hamiltonian is built from the computed ro-vibrational line list of the molecule in question. The Hamiltonian we propose is general and suitable for any polar molecule in the presence of an electric field. By exploiting precomputed data, the often prohibitively expensive computations associated with high accuracy simulations of molecule-field interactions are avoided. Applications to strong terahertz field-induced ro-vibrational dynamics of PH3 and NH3, and spontaneous emission data for optoelectrical Sisyphus cooling of H2CO and CH3Cl are discussed.

Molecule-field interactions are central to much current research in molecular physics. New methods are actively being developed to manipulate the internal and external degrees of freedom, ro-vibrational state populations and coherences of molecules1. Although several experimental techniques have emerged for controlling molecular dynamics, the variety of molecules to which they have been applied to is limited. Interest in a broader range of larger systems2 along with the increasing precision of instruments has created a demand for accurate theoretical approaches. At present, with a few exceptions3,4,5,6, molecule-field interactions are described by simple minimal-coupling theoretical models (e.g. rigid rotor, harmonic oscillator). More sophisticated treatments which consider all major electronic, nuclear motion, and external field effects are therefore highly desirable.

Progress in this direction is being made but very few general purpose computer programs exist for quantum mechanical modelling of molecule-field interactions. The multiconfiguration time-dependent Hartree (MCTDH) method7,8,9 is one possible option, and a robust, variational approach, RichMol, will soon be published. However, for high accuracy applications, calculations can be prohibitively expensive even for small to medium sized molecules. In this work we circumvent this issue by utilizing highly accurate data, which is freely available from the ExoMol database10,11, to model phenomena associated with polar molecules in the presence of an electric field.

A considerable amount of work goes into creating and maintaining spectroscopic databases10,11,12,13,14,15,16,17 and they contain a wealth of molecular information. From a theoretical standpoint it can take years to construct a comprehensive line list suitable for high-resolution spectroscopy. The process itself employs a combination of first principles and empirically tuned quantum mechanical methods18. Because of the enormous number of transitions which must be considered, often billions for small, polyatomic molecules, the task is computationally intensive and requires extensive use of high-performance computing facilities.

The ExoMol database10,11 contains detailed line lists for a number of important diatomic and polyatomic molecules which have been treated using nuclear motion programs specially optimized for such calculations19. Originally set up to aid the characterisation and modelling of cool stars and exoplanets10, the data has found other uses beyond atmospheric applications. These include the calculation of molecular partition functions20, radiative lifetimes of ro-vibrational states21, the study of magnetic field effects (Zeeman splittings)22, and investigating the sensitivity of spectral lines to a possible variation of the proton-to-electron mass ratio23,24.

Here we present two field-dependent applications which exploit the same high accuracy data. The two phenomena we consider are terahertz-induced two-level coherence and field-free orientation with application to PH3 and NH3, and spontaneous emission data for Sisyphus cooling of H2CO and CH3Cl. The Hamiltonians required to model these situations can be built from computed ro-vibrational transitions which are, or will soon be in the case of CH3Cl, available from the ExoMol website (www.exomol.com). Although we only discuss these two examples the external field Hamiltonian we describe is suitable for a range of molecule-field interaction applications such as double resonance spectroscopy25, frequency chirp excitation techniques26, alignment and orientation using static and laser fields27,28, molecular trapping and cooling29.

Results

Strong terahertz field induced ro-vibrational dynamics of PH3 and NH3

Recent advances in high energy laser-based terahertz (THz) technology have opened up exciting opportunities in THz coherent control, field-free orientation and alignment. Generally speaking it is considerably more challenging to create oriented, as opposed to aligned, molecular samples. Transient field-free alignment usually requires non-resonant, moderately intense laser fields to be applied several times throughout an experiment to maintain the coherent evolution of a rotational wavepacket under field-free conditions (see ref. 30 for a review). Alignment is achieved by exciting Raman transitions between ro-vibrational states with the same parity. For field-free orientation, both even- and odd-parity angular momentum states must be coherently excited. This has been demonstrated experimentally using non-resonant, intense two color laser pulses30,31,32, linearly polarized pulses with 45°-skewed mutual polarization33,34, an optical centrifuge35, and intense THz pulses36,37,38.

In contrast to non-resonant laser fields, single THz pulses contain frequency components resonant to the energy spacings between adjacent, opposite-parity rotational states with ΔJ = ±1. This gives rise to transient molecular orientation which manifests, for example, in emission of free induction decay signals36. A number of theoretical and experimental studies have reported strong THz field induced molecular alignment and orientation, rotational coherence and birefringence for diatomic (e.g. CO, HF, HBr)37,39,40 and linear triatomic (e.g. OCS)36,38,39,41 molecules. In this section we present theoretical simulations of THz-induced ro-vibrational dynamics for PH3 and NH3. Field-free alignment and orientation of phosphine is demonstrated using four time-delayed THz pulses, and by employing a combination of THz and intense non-resonant Raman pulses. For ammonia we focus on creating and manipulating the rotational coherences and populations by time-delayed THz pulse pairs.

The time-dependent quantum simulations of THz-molecule interactions were carried out using the electric field Hamiltonian in the dipole approximation as given by Eq. (2). The ro-vibrational energies εl,J and dipole matrix elements An external file that holds a picture, illustration, etc.
Object name is srep45068-m1.jpg where the polarization axis A = X, Y, Z, were extracted from the room temperature line lists for PH342 and NH343 (available from the ExoMol database; the much larger hot line lists43,44 are not necessary for the present applications). The ro-vibrational basis was truncated at the maximal rotational quantum number J = 20 and the maximal vibrational quanta vi corresponding to the ground and first excited vibrational states, i.e. vi = 0, 1. For NH3, due to strong anharmonicity of the large amplitude umbrella-type vibrational motion we extended the vibrational basis set to include vumb = 0, …, 7.

For both molecules only one quantum wavepacket simulation corresponding to the initial ground ro-vibrational state was performed. For PH3 this was the v = 0; J = 0, k = 0 state, where the quantum number k is the projection of the total angular momentum onto the molecule-fixed z axis. For NH3, due to nuclear spin statistics, the initial state was the v = 0+, J = 0, k = 0 state. In order to investigate the influence of temperature T on alignment and orientation, simulations were carried out at T = 0 K and T = 10 K for PH3. At T = 10 K a total of 81 initial states are populated. Final results were therefore averaged over all quantum trajectories from the different initial states.

The single-pulse THz electric field was represented by the analytical function39

An external file that holds a picture, illustration, etc.
Object name is srep45068-m2.jpg

where E0 is the amplitude of the pulse, An external file that holds a picture, illustration, etc.
Object name is srep45068-m3.jpg is the pulse duration, and An external file that holds a picture, illustration, etc.
Object name is srep45068-m4.jpg is the central frequency. In all calculations we fixed An external file that holds a picture, illustration, etc.
Object name is srep45068-m5.jpg to 4 ps. For simulations of the mixed Raman-THz field-free orientation of PH3, the lack of a polarizability model in the ExoMol database meant the simulations of the initial Raman excitation were carried out using a more general variational approach implemented in the TROVE6 and RichMol program packages. The resulting wavepacket was then used as the initial wavefunction in subsequent THz-driven dynamics simulations, modelled using the ExoMol line list data. Note that all fields are polarized along the Z axis in calculations.

PH3

The results of our numerical simulations of PH3 at T = 0 K and T = 10 K are shown in Figs 1 and and2,2, respectively. Four time-delayed THz pulses with central frequency An external file that holds a picture, illustration, etc.
Object name is srep45068-m6.jpg ps−1 and peak field strengths at E0 = 250 kV/cm and E0 = 1 MV/cm have been used to induce transient alignment left angle bracketcos2θright angle bracket (t) and orientation left angle bracketcosθright angle bracket (t), where θ is the Euler angle. Such intense THz pulses in the few-cycle regime can be routinely generated in a convenient tabletop setup45,46,47. As expected, increasing the number of pulses significantly enhances the degree of orientation and alignment for both temperatures and peak field strengths. The periodic behaviour of the orientation dynamics follows the quantum rotational revival pattern; strong peaks are observed at times t  2.0, 4.2, 5.7, 7.9, 9.5 ps etc., which is approximately 1/2, 1, 3/2, 2, 5/2, … of a revival time Trev = 1/2 Bc  3.7 ps (B  4.5 cm−1 is the rotational constant of PH3 and c is the speed of light). The enhancement of the alignment and orientation at delay times τ  4, 8, 12 ps etc., reveals recurrences separated by approximately Δτ = Trev.

Figure 1
Alignment and orientation of PH3 initially in the ground state at T = 0 K induced by a sequence of four time-delayed THz pulses (with central frequency An external file that holds a picture, illustration, etc.
Object name is srep45068-m55.jpg ps−1) plotted for different time delays τ between ...
Figure 2
Alignment and orientation of PH3 initially in thermal equilibrium at T = 10 K induced by a sequence of four time-delayed THz pulses (with central frequency An external file that holds a picture, illustration, etc.
Object name is srep45068-m57.jpg ps−1) plotted for different time delays τ between ...

Increasing the peak field strength of the THz pulses results in a more complex pattern of alignment and orientation with additional local minima and maxima. These secondary peaks are always present but attenuated in the weak field regime. Alignment and orientation are enhanced in stronger electric fields, hence higher absolute peak values of left angle bracketcos2θright angle bracket (t) and left angle bracketcosθright angle bracket (t). However, raising the temperature from T = 0 K to T = 10 K activates several new quantum trajectories via thermally populated states and this leads to overall decoherence of the wavepacket. This effect is responsible for the lower absolute values of the orientation and alignment parameters at T = 10 K in Fig. 2. Therefore, losses in orientation or alignment at higher temperatures can be compensated for with a stronger peak intensity of the THz pulse.

As an alternative to multiple THz pulses, in ref. 38 it was demonstrated that the degree of alignment and orientation can be improved by applying a short intense Raman pulse to coherently prepare molecules in highly excited rotational states prior to THz exposure. This allows a larger molecular population to occupy states for which transition frequencies match the peak of the THz pulse envelope. Raman excitation follows ΔJ = ±2 selection rules and hence couples states of the same parity. Therefore, when starting from an initial J = 0 state, the Raman pulse predominantly excites J = 2 states. The subsequent THz pulse induces J  J′: 0 →1, 2  3 transitions and produces a mixed-parity wavepacket, thus enabling effective net orientation. The results of our simulations for PH3 are shown in Fig. 3. The alignment pattern is insensitive to the time delay between the Raman and THz pulse and it appears the THz pulse contributes only marginally to alignment. For a given time delay the orientation peaks are separated by ≈3.9 ps, which corresponds to Trev.

Figure 3
Alignment and orientation of PH3 initially in the ground state induced by a time-delayed Raman Gaussian pulse (E0 = 108 V/cm, ω = 800 nm,FWHM = 0.4 ps) and THz pulse (An external file that holds a picture, illustration, etc.
Object name is srep45068-m58.jpg ps ...

NH3

As a second example we explore strategies to establish control over rotational coherences and populations in NH3 by applying optimally timed THz pulses. Recently it was shown both theoretically and experimentally41 that a pair of THz pulses separated in time can induce larger population transfer (larger transient emission responses) than two interactions within a single short THz pulse.

To adjust the THz pulse parameters to the energy level structure of NH3, we first look at the Fourier intensity of a pair of THz pulses at a few characteristic frequencies as a function of the time delay τ between the pulses and the pulse central frequency An external file that holds a picture, illustration, etc.
Object name is srep45068-m7.jpg (see Eq. (1)). The results are shown in Fig. 4 for the three rotational transitions, |0, 1right angle bracket  |1, 0right angle bracket, |1, 0right angle bracket  |2, 1right angle bracket, and |2, 1right angle bracket  |3, 0right angle bracket, where states of NH3 are labelled as |J, vumbright angle bracket, and the quantum numbers k and m, which correspond to the projection of the total angular momentum onto the molecule-fixed and laboratory-fixed axes, respectively, are zero. It is evident that the optimal values of An external file that holds a picture, illustration, etc.
Object name is srep45068-m8.jpg are 0.4, 1.1, and 1.6 ps−1 for the three different rotational excitations of NH3.

Figure 4
Fourier intensity of two THz pulses separated by a delay time τ versus central frequency An external file that holds a picture, illustration, etc.
Object name is srep45068-m59.jpg of a single pulse (seeEq. (1)).

Figure 5 depicts the populations of the ro-vibrational wave packets of NH3, initially in the |0, 1right angle bracket, |1, 0right angle bracket, and |2, 1right angle bracket ro-vibrational states, following the excitation by a time-delayed THz pulse pair with E0 = 500 kV/cm and respective pulse central frequencies An external file that holds a picture, illustration, etc.
Object name is srep45068-m9.jpg, 1.1, and 1.6 ps−1. Population transfer is significantly enhanced by two properly delayed THz pulses when compared to a single, but twice as intense, THz pulse (τ = 0). Slight adjustment of the central frequency of the THz pulses increases the selectivity of excitation, suppressing (or enhancing) two-photon transitions. As expected, population modulations for the three excited rotational states with J = 1, 2 and 3, shown in the (a)–(c) panels of Fig. 5, appear with a respective revival period of Trev  J/2Bc  1.7, 0.8 and 0.6 ps, respectively, where B  6.35 cm−1 is the rotational constant of NH3.

Figure 5
Ro-vibrational wave packet populations of NH3 after excitation using a pair of THz pulses with E0 = 500 kV/cm, An external file that holds a picture, illustration, etc.
Object name is srep45068-m63.jpg ps, (a) An external file that holds a picture, illustration, etc.
Object name is srep45068-m64.jpg ps−1, (b) An external file that holds a picture, illustration, etc.
Object name is srep45068-m65.jpg ps−1, (c) An external file that holds a picture, illustration, etc.
Object name is srep45068-m66.jpg ps-1, as a function of delay time τ ...

Based on the results of Fig. 5 we can combine pairs of THz pulses to create coherences between different pairs of rotational states with even and odd parities. This is shown in Fig. 6. The first pulse at t  1.1 ps induces a J  J′:0  1 transition, creating a coherent superposition An external file that holds a picture, illustration, etc.
Object name is srep45068-m10.jpg. The second pulse is tuned to affect population of the |1, 0right angle bracket state and transforms the two-state coherence into a new one: An external file that holds a picture, illustration, etc.
Object name is srep45068-m11.jpg. Parameters of the subsequent pulses (upper inset in Fig. 6) are optimized to maintain uniform division of molecular population between the two states. This illustrates the high level of control which can be exerted over the ro-vibrational wavepacket with the use of THz pulses.

Figure 6
Population transfer between pairs (J/J′) of rotational states of NH3 driven by a sequence of five time-delayed THz pulse pairs.

To further investigate the validity of our approach we have computed the alignment and orientation of NH3 using three different Hamiltonian models. The first, given in Eq. (2), is the electric field Hamiltonian in the dipole approximation. The second Hamiltonian we use incorporates the polarizability α into the model, while the third Hamiltonian also includes the hyperpolarizability β. Latter calculations were perfomed using our recently developed variational approach RichMol. The THz field was composed of three simultaneously applied pulses with respective central frequencies An external file that holds a picture, illustration, etc.
Object name is srep45068-m12.jpg, 0.6, and 0.9 ps−1 each with duration An external file that holds a picture, illustration, etc.
Object name is srep45068-m13.jpg ps. The results are illustrated in Fig. 7 for electric field strengths up to 10 MV/cm. Calculations were performed on an initial ground vibrational state with J = 1, k = 1, m = 0 and mixed parity +/−; that is An external file that holds a picture, illustration, etc.
Object name is srep45068-m14.jpg.

Figure 7
Orientation of NH3 driven by THz pulses with electric field strengths E0 up to 10 MV/cm range.

As shown in Fig. 7 it is only for the very intense 10 MV/cm THz pulses that we see contributions from the polarizability and hyperpolarizability, otherwise their effects are negligible. In general, alignment and orientation caused by the hyperpolarizability is an order of magnitude smaller than the polarizability and contributes with opposite sign. Field strengths as strong as 10 MV/cm, where the non-linear THz effects of α and β, although noticeable, are small enough to be neglected, can therefore be modelled using the Hamiltonian in the dipole approximation in conjunction with ExoMol line list data.

Intense THz pulses offer a promising approach for controlling large-amplitude motions in floppy molecules. In Fig. 8 we investigate the effect of an intense THz-pulse-train on the dynamics of the large-amplitude vibrational coordinate of NH3 associated with the umbrella inversional motion. This was performed with 50 MV/cm intense quasi half-cycle THz pulses, obtained in Eq. (1) by setting the single-pulse central frequency An external file that holds a picture, illustration, etc.
Object name is srep45068-m15.jpg ps−1 and pulse duration An external file that holds a picture, illustration, etc.
Object name is srep45068-m16.jpg ps. The left panel of Fig. 8 shows the expectation value of the inversion coordinate ρinv (zero at planar molecular geometry) as a function of time and delay time τ between single pulses. In the region where the delay time is smaller than the single pulse duration, i.e. An external file that holds a picture, illustration, etc.
Object name is srep45068-m17.jpg ps, the inversion tunnelling is significantly extended or prohibited by the effects of the strong electric field. For many experiments on controlled molecules the presence of an intense electric field can lead to undesirable side effects. To avoid them, the delay time between THz pulses must exceed the single pulse duration, An external file that holds a picture, illustration, etc.
Object name is srep45068-m18.jpg, creating short windows in time with no strong electric field present. As seen in Fig. 8, for delay times just above the pulse duration it is possible to modulate, at least to some extent, the inversion tunnelling rate. For larger delay times the tunnelling dynamics quickly approach field-free behaviour.

Figure 8
The effect of intense quasi half-cycle THz pulses (E0 = 50 MV/cm, An external file that holds a picture, illustration, etc.
Object name is srep45068-m70.jpg ps−1, An external file that holds a picture, illustration, etc.
Object name is srep45068-m71.jpg ps) on the dynamics of the large-amplitude vibrational coordinate ρinv of NH3.

To again validate the use of ExoMol line list data with the Hamiltonian in the dipole approximation for this type of observable, we have run calculations with an extended Hamiltonian containing non-linear contributions from the polarizability and hyperpolarizability tensors. The results are shown in the right panel of Fig. 8 and confirm the validity of the ExoMol dipole-only approach to predict non-rigid dynamics with THz fields below 50 MV/cm.

Spontaneous emission data for Sisyphus cooling of H2CO and CH3Cl

Optoelectrical Sisyphus cooling is a robust and general method for producing ultracold polyatomic molecules48,49,50. This technique works by moving molecules through a closed system of trapped ro-vibrational states in the presence of an electric field. Molecules lose kinetic energy as they travel up the electric field gradient of the trapping potential. The cycle is repeated until sufficient cooling has taken place with sub-millikelvin temperatures possible. The speed and efficiency of Sisyphus cooling is dictated by the Einstein A coefficients and decay channels of the involved ro-vibrational energy levels; information which is readily available and computed to a high degree of accuracy in the ExoMol database. It is possible then to identify suitable transitions in different molecular systems which could be used for Sisyphus cooling. This type of analysis can aid future experiments if there is interest in a particular molecule, potentially leading to improved rates of cooling.

A general overview of the Sisyphus cooling level scheme is shown in Fig. 9. Starting from an initial, preferably highly populated ro-vibrational state, molecular population is transferred to an excited bridge state via an infrared (IR) transition. Spontaneous decay from the bridge state should be dominated by at most two channels to long-lived low-field seeking states. Large branching ratios are essential in this step with the Einstein A coefficient of the primary decay channel dictating the speed of the cooling cycle. As molecules penetrate the high-field region they climb up the Stark induced energy gradient and lose kinetic energy. For larger values of the quantum number M, the gradient is steeper and more energy is removed. It is important that the Stark split M sublevels are low-field seeking states (shifted towards higher energies) as this ensures the molecules remain trapped. Applied radio-frequency (RF) radiation stimulates emission to lower M quantum number states and the molecules gain a small amount of kinetic energy as they move down the electric field gradient. However, this gain is smaller than the previous loss in energy and the temperature of the molecules is reduced. A microwave (MW) field connects the two decay channels and closes the cycle. The process is repeated until sufficient cooling has occurred.

Figure 9
Overview of Sisyphus cooling scheme.

Sisyphus cooling has been realised experimentally for methyl fluoride (CH3F)49 and formaldehyde (H2CO)50. However, a number of symmetric top molecules with strong parallel vibrational transitions have been proposed as suitable candidates for cooling using this method48. The closed system of ro-vibrational energy levels should satisfy the following criteria:

  • The initial state is highly populated and long-lived. Energy levels in the vibrational ground state with rotational quantum number J in the range 2  J  4 are generally good candidates.
  • There must be a sufficiently fast (intense) IR transition from the initial state to an excited bridge state with an Einstein A coefficient large enough that population transfer can occur in the time frame of effective trapping.
  • State-selective spontaneous decay is possible from the bridge state with one or two dominant decay channels. The branching ratio of the two channels should sum to near unity.
  • Target states of spontaneous decay must be low-field seeking states with sufficiently long lifetimes to ensure the molecules remain trapped.
  • RF transitions from target states to lower M quantum number states are available.
  • MW radiation can be used to move population from the secondary decay channel back to the initial state to close the cycle.

To illustrate how the ExoMol database, or in general any reasonably “complete” spectroscopic library, can be used to identify molecular transitions for Sisyphus cooling we have performed analysis on H2CO and CH3Cl. These systems have been treated within the ExoMol framework19 and we refer the reader to the relevant publications for details of the ro-vibrational calculations51,52,53,54. Candidate ro-vibrational states were identified by sorting the An external file that holds a picture, illustration, etc.
Object name is srep45068-m19.jpg files according to upper state energy. There are a huge number of energy levels to process so to reduce the workload we only considered states involved in at least one intense transition (Einstein A coefficient of the order 101 or greater) and with J  4. Once identified, all transitions to the upper state were collected and branching ratios computed. Energy levels with two dominant decay channels or less were then analysed to check that the lower states (or target states) were low-field seeking.

H2CO

In the work of ref. 50, Sisyphus cooling of formaldehyde utilized ro-vibrational energy levels of the v1 symmetric C–H stretching mode, |v1, 3, 3right angle bracket, and of the vibrational ground state, |0, 3, 3right angle bracket and |0, 4, 3right angle bracket. Here energy levels are labelled as |Mode, J, Kright angle bracket where K = |KA| and we omit the quantum number KC (see discussion in ref. 50). The primary decay channel |v1, 3, 3right angle bracket  |0, 3, 3right angle bracket had a branching ratio of γ = 0.75 and an estimated spontaneous decay rate (equivalent to the Einstein A coefficient) of approximately 60 Hz, which determines the speed of cooling.

In Table 1 we list six collections of suitable states for cooling including those used by ref. 50 (third and forth row). Using the ExoMol database we found that the |v1, 3, 3right angle bracket level has 41 decay channels, of which only transitions to the |0, 3, 3right angle bracket and |0, 4, 3right angle bracket levels were dominant. We predict a slightly smaller Einstein A coefficient for the |v1, 3, 3right angle bracket  |0, 3, 3right angle bracket transition and hence a slightly slower cooling rate than ref. 50. However, the difference is small (≈10 Hz) and our computed branching ratios are identical to those of ref. 50.

Table 1
Transition frequencies vfi (cm−1), Einstein A coefficients Afi (s−1), and branching ratios γ for Sisyphus cooling of H2CO.

For other possible candidate states in the v1 manifold, a slightly higher branching ratio and Einstein A coefficient are predicted for the |v1, 4, 4right angle bracket  |0, 4, 4right angle bracket transition. Alternatively, the v2 C–O stretching mode appears suitable for Sisyphus cooling but using these states would half the speed of cooling compared to the v1 levels. All lower states E″ in Table 1 belong to the vibrational ground state and thus remain accessible either thermally or by methods of selective state preparation. At T = 100 K for example, lower states are thermally populated to around 30%. Note that only transitions allowed by selection rules were considered, i.e. ΔJ, ΔM = 0, ±1 and ΔK = 0.

CH3Cl

Methyl chloride has not yet been considered for Sisyphus cooling but we expect a similar experiment to CH3F49 could be performed. As shown in Table 2, ro-vibrational levels of the v1 symmetric CH3 stretching mode and vibrational ground state appear the best candidates for cooling. Similar to methyl fluoride, which had a spontaneous decay rate of about 15 Hz for the |v1, 3, 3right angle bracket  |0, 3, 3right angle bracket transition, the respective rate in CH3Cl is only marginally larger and would hence provide a similar speed of cooling. Note that we have only considered CH335Cl for the present analysis but we would expect the same conclusions for CH337Cl.

Table 2
Transition frequencies vfi (cm−1), Einstein A coefficients Afi (s−1) and branching ratios γ for Sisyphus cooling of CH3Cl.

Discussion

Comprehensive spectroscopic line lists generated as part of the ExoMol project can be used to model ro-vibrational phenomena for small polyatomic molecules in the presence of external electric fields. By doing so we avoid repeating the same, often expensive calculations. Currently the ExoMol database contains complete sets of high accuracy ro-vibrational energies, transition frequencies, Einstein A coefficients, and complex dipole moment phases for many important diatomic and small polyatomic molecules valid for temperatures up to T = 1500 K. These data can be straightforwardly utilized to model molecule-field interactions in the dipole approximation. In the future we plan to extend the ExoMol database with Raman and electric quadrupole transition moments for selected molecules, such as NH3 and H2O, thereby extending the range of possible applications.

Two illustrative examples of molecule-field interaction applications were presented. Strong THz field induced ro-vibrational dynamics of PH3 and NH3 were simulated and these represent the first high accuracy calculations on polyatomic molecules which have been reported in the literature. For optoelectrical Sisyphus cooling, suitable collections of states were identified in H2CO and CH3Cl. In both molecules the v1 symmetric stretching mode provided the fastest cooling rates and most suitable ro-vibrational energy levels for Sisyphus cooling.

Although we have chosen relatively straightforward processes to look at, our approach can simulate any case with no constraint on the availability of data (provided the molecule is in the ExoMol database). At present, we are unaware of any other line list database which stores the dipole moment phase factors. This information is always produced with the dipole line strengths or Einstein coefficients but is discarded when the line lists are compiled, thus making them unusable for modeling field-driven effects. We therefore encourage other theory-based data compilations to retain this information in future.

Methods

ExoMol line list data structure

A complete description of the ExoMol data structure along with examples was recently reported in ref. 11. Therefore we provide only a brief description that is relevant for the present work. The two main files available for download from the ExoMol website (www.exomol.com) are the An external file that holds a picture, illustration, etc.
Object name is srep45068-m20.jpg and An external file that holds a picture, illustration, etc.
Object name is srep45068-m21.jpg files, unique to each molecule. The An external file that holds a picture, illustration, etc.
Object name is srep45068-m22.jpg file contains all computed ro-vibrational energy levels (in cm−1). For polyatomics, which are the focus of the present study, each energy level has a unique state ID with symmetry and quantum number labelling. The An external file that holds a picture, illustration, etc.
Object name is srep45068-m23.jpg files, which are split into frequency windows as to be more manageable, contain all calculated molecular transitions. Upper and lower state labels, transition frequencies and Einstein A coefficients are provided.

Modelling Stark associated phenomena requires information on the complex phase of the dipole moment matrix elements (discussed below). This information, simply a + or −, is not currently available for all molecules in the ExoMol database. However, work is underway to rectify this and the complex phase for each molecular transition can be extracted from the An external file that holds a picture, illustration, etc.
Object name is srep45068-m24.jpg file. Several other molecule specific files related to atmospheric applications (e.g. pressure broadening of spectral lines) are available but since we do not use them in this work, we do not discuss them here.

Electric field interaction Hamiltonian

The complete Hamiltonian is written as a sum of the field-free ro-vibrational Hamiltonian and the molecule-field interaction Hamiltonian, which is approximated in this work as the permanent dipole interaction

An external file that holds a picture, illustration, etc.
Object name is srep45068-m25.jpg

Here εl,J and |l, J, mright angle bracket denote the field-free ro-vibrational energy and wavefunction, respectively, of a state identified by the running number l (ID number) in the line list energy file (An external file that holds a picture, illustration, etc.
Object name is srep45068-m26.jpg file), J is the rotational quantum number of the total angular momentum, and the quantum number m = J, …, J is the projection of the total angular momentum on the laboratory-fixed Z axis. The coupling strength is given by the product of the ro-vibrational matrix elements of electronic ground state electric dipole moment An external file that holds a picture, illustration, etc.
Object name is srep45068-m27.jpg in the ro-vibrational basis and the time-dependent classical electric field EA(t), defined in the laboratory fixed axes system. To describe the radiative fields we use the form An external file that holds a picture, illustration, etc.
Object name is srep45068-m28.jpg, where A is the polarization axis, An external file that holds a picture, illustration, etc.
Object name is srep45068-m29.jpg is the pulse time profile with a maximum at t = t0 and duration (FWHM) T, and ωA is the carrier frequency.

Calculation of the An external file that holds a picture, illustration, etc.
Object name is srep45068-m30.jpg ro-vibrational matrix elements from Einstein A coefficients, extracted from the ExoMol line list, is described in detail below. The selection rules for An external file that holds a picture, illustration, etc.
Object name is srep45068-m31.jpg allow coupling between ro-vibrational states with ΔJ = J  J = 0, ±1 and Δm = m  m = 0, ±1. Hence, the total Hamiltonian may be constructed in the form of a block-tridiagonal matrix, having blocks corresponding to Δm = −1, 0, or +1 in the lower, main, and upper diagonals, respectively. Using this representation, various linear algebra operations can be performed efficiently and the memory requirements reduced by use of the band matrix storage scheme. For A = Z, An external file that holds a picture, illustration, etc.
Object name is srep45068-m32.jpg can couple only states with Δm = 0, thus the total Hamiltonian becomes factorized into independent blocks for each chosen m, which are processed independently.

The time-evolution of the wavefunction is given by the time-dependent Schrödinger equation with the Hamiltonian in Eq. (2), which may be time-independent or time-dependent. For static fields the problem reduces to solving an eigenvalue equation for the Hamiltonian. For radiative fields the time-evolution of the wavefunction is described by the time-evolution operator U(t, t′) as An external file that holds a picture, illustration, etc.
Object name is srep45068-m33.jpg. Providing Δt = t  t′ is sufficiently small compared to the characteristic oscillation period of the field perturbation, U(t, t′) can be evaluated, for example, using the split-operator method as

An external file that holds a picture, illustration, etc.
Object name is srep45068-m34.jpg

where the elements of the matrices H0 and An external file that holds a picture, illustration, etc.
Object name is srep45068-m35.jpg are given by An external file that holds a picture, illustration, etc.
Object name is srep45068-m36.jpg and An external file that holds a picture, illustration, etc.
Object name is srep45068-m37.jpg, respectively. The bottleneck operation in Eq. (3) is the evaluation of the matrix exponential at each time step t = t + Δt/2, which in our case is computed using the iterative approximation based on Krylov subspace methods, as implemented in the Expokit computational package55. Alternative and more sophisticated representations for the split-operator technique can be found in refs 56 and 57.

Dipole matrix elements

We aim to derive the relationship between the Einstein A coefficients, An external file that holds a picture, illustration, etc.
Object name is srep45068-m38.jpg, and the matrix elements of the electric dipole moment operator, An external file that holds a picture, illustration, etc.
Object name is srep45068-m39.jpg (An external file that holds a picture, illustration, etc.
Object name is srep45068-m40.jpg), defined relative to the laboratory fixed system. It is convenient to first express An external file that holds a picture, illustration, etc.
Object name is srep45068-m41.jpg in terms of the dipole line strength An external file that holds a picture, illustration, etc.
Object name is srep45068-m42.jpg as

An external file that holds a picture, illustration, etc.
Object name is srep45068-m43.jpg

where An external file that holds a picture, illustration, etc.
Object name is srep45068-m44.jpg is the f  i transition frequency in units of cm−1, gns is the nuclear spin statistical weight factor, J′ refers to the rotational quantum number of the final state f, and the units of Afi and Sfi are s−1 and D2, respectively.

Writing the wavefunctions for the initial and final states as products of rotational and vibrational basis functions, i.e. An external file that holds a picture, illustration, etc.
Object name is srep45068-m45.jpg and An external file that holds a picture, illustration, etc.
Object name is srep45068-m46.jpg, the expression for Sfi reads58

An external file that holds a picture, illustration, etc.
Object name is srep45068-m47.jpg

where v and J, k, m refer to the vibrational and rotational quanta, respectively, μα (α = x, y, z) is the dipole moment in the molecule-fixed axes system, and Tσα is the Cartesian-to-spherical tensor transformation matrix

An external file that holds a picture, illustration, etc.
Object name is srep45068-m48.jpg

It should be noted that the actual wavefuncitons |ψright angle bracket are expansions in terms of |J, k, mright angle bracket|vright angle bracket as basis functions. We omit the expansion coefficients in Eq. (5) for the sake of simplicity, which does not alter the methodology presented below. Using the same variables, the matrix element of the dipole moment operator An external file that holds a picture, illustration, etc.
Object name is srep45068-m49.jpg may be expressed as

An external file that holds a picture, illustration, etc.
Object name is srep45068-m50.jpg

By combining Eq. (7) with Eq. (4) and Eq. (5), it is straightforward to write down An external file that holds a picture, illustration, etc.
Object name is srep45068-m51.jpg in terms of the quantities Sfi and Afi provided by the ExoMol line list. Because of the complex-valued expression under the modulus in Eq. (5), an additional number retaining the complex phase must be supplemented for each ro-vibrational transition f  i in a line list alongside the value of Afi.

For the sake of simplicity the above formulas were derived using the symmetric top basis functions |J, k, mright angle bracket, which yield complex phases in Eq. (5) ranging from 0 to 2π. In the variational calculations used to generate the ro-vibrational line list, which for polyatomics is performed with the TROVE59,60 and RichMol program packages, we employ the basis of symmetrized combinations of the rigid rotor functions An external file that holds a picture, illustration, etc.
Object name is srep45068-m52.jpg for k > 0 and |J, 0, 0, mright angle bracket for k = 0. Using the symmetrized basis to derive Eqs (5) and (7), the expression under the modulus will become purely imaginary, hence the corresponding complex phase factor is ±1. By defining the complex phase factor as An external file that holds a picture, illustration, etc.
Object name is srep45068-m53.jpg with Pfi = 0 or 1, the final expression for the dipole matrix elements reads as

An external file that holds a picture, illustration, etc.
Object name is srep45068-m54.jpg

Additional Information

How to cite this article: Owens, A. et al. Simulating electric field interactions with polar molecules using spectroscopic databases. Sci. Rep. 7, 45068; doi: 10.1038/srep45068 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgments

This work was supported by FP7-MC-IEF Project 629237, ERC Advanced Investigator Project 267219, and by STFC Project ST/J002925. We also thank the DiRAC@Darwin HPC cluster, which is the UK HPC facility for particle physics, astrophysics and cosmology and is supported by STFC and BIS.

Footnotes

The authors declare no competing financial interests.

Author Contributions A.Y. conceived and designed the work, and performed the terahertz induced ro-vibrational dynamics simulations. A.O., E.J.Z. and K.L.C. carried out the analysis related to optoelectrical Sisyphus cooling. A.O., E.J.Z., K.L.C., S.N.Y., J.T. and A.Y. all contributed to the writing of the manuscript.

References

  • Lemeshko M., Krems R. V., Doyle J. M. & Kais S. Manipulation of molecules with electromagnetic fields. Mol. Phys. 111, 1648–1682, doi: (2013).10.1080/00268976.2013.813595 [Cross Ref]
  • Küpper J. et al. . Instrumentation and methods: general discussion. Faraday Discuss. 171, 505–523 doi: (2014).10.1039/C4FD90017K [PubMed] [Cross Ref]
  • Luis J. M., Reis H., Papadopoulos M. & Kirtman B. Treatment of nonlinear optical properties due to large amplitude anharmonic vibrational motions: Umbrella motion in NH3. J. Chem. Phys. 131, doi: (2009).10.1063/1.3171615 [PubMed] [Cross Ref]
  • Reis H., Luis J. M., Garcia-Borrás M. & Kirtman B. Computation of nonlinear optical properties of molecules with large amplitude anharmonic motions. III. Arbitrary double-well potentials. J. Chem. Theory. Comput. 10, 236–242, doi: (2014).10.1021/ct400938a [PubMed] [Cross Ref]
  • Coudert L. H. Modeling torsional control of molecules with four-dimensional, two-dimensional, and one-dimensional approaches. Phys. Rev. A 91, 013402, doi: (2015).10.1103/PhysRevA.91.013402 [Cross Ref]
  • Yachmenev A. & Yurchenko S. N. Detecting chirality in molecules by linearly polarized laser fields. Phys. Rev. Lett. 117, 033001, doi: (2016).10.1103/PhysRevLett.117.033001 [PubMed] [Cross Ref]
  • Beck M. H., Jäckle A., Worth G. A. & Meyer H.-D. The multiconfiguration time-dependent Hartree (MCTDH) method: a highly efficient algorithm for propagating wavepackets. Phys. Rep. 324, 1–105, doi: (2000).10.1016/S0370-1573(99)00047-2 [Cross Ref]
  • Meyer H.-D. & Worth G. A. Quantum molecular dynamics: propagating wavepackets and density operators using the multiconfiguration time-dependent Hartree method. Theor. Chem. Acc. 109, 251–267, doi: (2003).10.1007/s00214-003-0439-1 [Cross Ref]
  • Meyer H.-D. Studying molecular quantum dynamics with the multiconfiguration time-dependent Hartree method. WIREs Comput. Mol. Sci. 2, 351–374 doi: (2012).10.1002/wcms.87 [Cross Ref]
  • Tennyson J. & Yurchenko S. N. ExoMol: molecular line lists for exoplanet and other atmospheres. Mon. Not. R. Astron. Soc. 425, 21–33 doi: (2012).10.1111/j.1365-2966.2012.21440.x [Cross Ref]
  • Tennyson J. et al. . The ExoMol database: molecular line lists for exoplanet and other hot atmospheres. J. Mol. Spectrosc. 327, 73–94, doi: (2016).10.1016/j.jms.2016.05.002 [Cross Ref]
  • Rothman L. S. et al. . The HITRAN2012 molecular spectroscopic database. J. Quant. Spectrosc. Radiat. Transf. 130, 4–50 doi: (2013).10.1016/j.jqsrt.2013.07.002 [Cross Ref]
  • Hill C. et al. . HITRANonline: An online interface and the flexible representation of spectroscopic data in the HITRAN database. J. Quant. Spectrosc. Radiat. Transf. 177, 4–14 doi: (2016).10.1016/j.jqsrt.2015.12.012 [Cross Ref]
  • Pickett H. M. et al. . Submillimeter, millimeter, and microwave spectral line catalog. J. Quant. Spectrosc. Radiat. Transf. 60, 883–890, doi: (1998).10.1016/S0022-4073(98)00091-0 [Cross Ref]
  • Müller H. S. P., Schlöder F., Stutzki J. & Winnewisser G. The Cologne Database for Molecular Spectroscopy, CDMS: a useful tool for astronomers and spectroscopists. J. Mol. Struct. 742, 215–227, doi: (2005).10.1016/j.molstruc.2005.01.027 [Cross Ref]
  • Jacquinet-Husson N. et al. . The 2009 edition of the GEISA spectroscopic database. J. Quant. Spectrosc. Radiat. Transf. 112, 2395–2445, doi: (2011).10.1016/j.jqsrt.2011.06.004 [Cross Ref]
  • Rey M., Nikitin A. V., Babikov Y. L. & Tyuterev V. G. TheoReTS - An information system for theoretical spectra based on variational predictions from molecular potential energy and dipole moment surfaces. J. Mol. Spectrosc. 327, 138–158, doi: (2016).10.1016/j.jms.2016.04.006 [Cross Ref]
  • Tennyson J. Accurate variational calculations for line lists to model the vibration-rotation spectra of hot astrophysical atmospheres. WIREs Comput. Mol. Sci. 2, 698–715, doi: (2012).10.1002/wcms.94 [Cross Ref]
  • Tennyson J. & Yurchenko S. N. The ExoMol project: Software for computing large molecular line lists. Int. J. Quantum Chem., 117, 92–103, doi: (2017).10.1002/qua.25190 [Cross Ref]
  • Sousa-Silva C., Hesketh N., Yurchenko S. N., Hill C. & Tennyson J. High temperature partition functions and thermodynamic data for ammonia and phosphine. J. Quant. Spectrosc. Radiat. Transf. 142, 66–74, doi: (2014).10.1016/j.jqsrt.2014.03.012 [Cross Ref]
  • Tennyson J., Hulme K., Naim O. K. & Yurchenko S. N. Radiative lifetimes and cooling functions for astrophysically important molecules. J. Phys. B: At. Mol. Opt. Phys. 49, 044002, doi: (2016).10.1088/0953-4075/49/4/044002 [Cross Ref]
  • Semenov M., Yurchenko S. N. & Tennyson J. Predicted Landé g-factors for open shell diatomic molecules. J. Mol. Spectrosc. doi: (2016).10.1016/j.jms.2016.11.004 [Cross Ref]
  • Owens A., Yurchenko S. N., Thiel W. & Špirko V. Accurate prediction of the ammonia probes of a variable proton-to-electron mass ratio. Mon. Not. R. Astron. Soc. 450, 3191–3200, doi: (2015).10.1093/mnrasl/stv869 [Cross Ref]
  • Owens A., Yurchenko S. N., Thiel W. & Špirko V. Enhanced sensitivity to a possible variation of the proton-to-electron mass ratio in ammonia. Phys. Rev. A 93, 052506, doi: (2016).10.1103/PhysRevA.93.052506 [Cross Ref]
  • Schmitz D., Shubert V. A., Patterson D., Krin A. & Schnell M. Phase dependence of double-resonance experiments in rotational spectroscopy. J. Phys. Chem. Lett. 6, 1493–1498 doi: (2015).10.1021/acs.jpclett.5b00494 [PubMed] [Cross Ref]
  • Rickes T. et al. . Efficient adiabatic population transfer by two-photon excitation assisted by a laser-induced Stark shift. J. Chem. Phys. 113, 534–546, doi: (2000).10.1063/1.481829 [Cross Ref]
  • Friedrich B. & Herschbach D. Manipulating molecules via combined static and laser fields. J. Phys. Chem. A 103, 10280–10288, doi: (1999).10.1021/jp992131w [Cross Ref]
  • Filsinger F. et al. . Quantum-state selection, alignment, and orientation of large molecules using static electric and laser fields. J. Chem. Phys. 131, 064309, doi: (2009).10.1063/1.3194287 [PubMed] [Cross Ref]
  • van de Meerakker S. Y. T., Bethlem H. L., Vanhaecke N. & Meijer G. Manipulation and control of molecular beams. Chem. Rev. 112, 4828–4878, doi: (2012).10.1021/cr200349r [PubMed] [Cross Ref]
  • Stapelfeldt H. & Seideman T. Colloquium: Aligning molecules with strong laser pulses. Rev. Mod. Phys. 75, 543, doi: (2003).10.1103/RevModPhys.75.543 [Cross Ref]
  • De S. et al. . Field-free orientation of CO molecules by femtosecond two-color laser fields. Phys. Rev. Lett. 103, 153002 doi: (2009).10.1103/PhysRevLett.103.153002 [PubMed] [Cross Ref]
  • Frumker E. et al. . Oriented rotational wave-packet dynamics studies via high harmonic generation. Phys. Rev. Lett. 109, 113901, doi: (2012).10.1103/PhysRevLett.109.113901 [PubMed] [Cross Ref]
  • Kitano K., Hasegawa H. & Ohshima Y. Ultrafast angular momentum orientation by linearly polarized laser fields. Phys. Rev. Lett. 103, 223002, doi: (2009).10.1103/PhysRevLett.103.223002 [PubMed] [Cross Ref]
  • Fleischer S., Khodorkovsky Y., Prior Y. & Averbukh I. S. Controlling the sense of molecular rotation. New J. Phys. 11, 105039, doi: (2009).10.1088/1367-2630/11/10/105039 [Cross Ref]
  • Korobenko A. & Milner V. Adiabatic field-free alignment of asymmetric top molecules with an optical centrifuge. Phys. Rev. Lett. 116, 183001, doi: (2016).10.1103/PhysRevLett.116.183001 [PubMed] [Cross Ref]
  • Fleischer S., Zhou Y., Field R. W. & Nelson K. A. Molecular orientation and alignment by intense single-cycle THz pulses. Phys. Rev. Lett. 107, 163603, doi: (2011).10.1103/PhysRevLett.107.163603 [PubMed] [Cross Ref]
  • Kitano K. et al. . Orientation of jet-cooled polar molecules with an intense single-cycle thz pulse. Phys. Rev. A 88, 061405 (2013). doi: .10.1103/PhysRevA.88.061405 [Cross Ref]
  • Egodapitiya K. N., Li S. & Jones R. R. Terahertz-induced field-free orientation of rotationally excited molecules. Phys. Rev. Lett. 112, 103002, doi: (2014).10.1103/PhysRevLett.112.103002 [PubMed] [Cross Ref]
  • Lapert M. & Sugny D. Field-free molecular orientation by terahertz laser pulses at high temperature. Phys. Rev. A 85, 063418 doi: (2012).10.1103/PhysRevA.85.063418 [Cross Ref]
  • Kitano K., Ishii N., Kanai T. & Itatani J. Selecting rotational two-level coherence in polar molecules by double terahertz pulses. Phys. Rev. A 90, 041402, doi: (2014).10.1103/PhysRevA.90.041402 [Cross Ref]
  • Fleischer S., Field R. W. & Nelson K. A. Commensurate two-quantum coherences induced by time-delayed THz fields. Phys. Rev. Lett. 109, 123603, doi: (2012).10.1103/PhysRevLett.109.123603 [PubMed] [Cross Ref]
  • Sousa-Silva C., Yurchenko S. N. & Tennyson J. A computed room temperature line list for phosphine. J. Mol. Spectrosc. 288, 28–37, doi: (2013).10.1016/j.jms.2013.04.002 [Cross Ref]
  • Yurchenko S. N., Barber R. J. & Tennyson J. A variationally computed line list for hot NH3. Mon. Not. R. Astron. Soc. 413, 1828–1834, doi: (2011).10.1111/j.1365-2966.2011.18261.x [Cross Ref]
  • Sousa-Silva C., Al-Refaie A. F., Tennyson J. & Yurchenko S. N. ExoMol line lists VIII: A hot line list for phosphine. Mon. Not. R. Astron. Soc. 446, 2337–2347 doi: (2015).10.1093/mnras/stu2246 [Cross Ref]
  • Yeh K.-L., Hoffmann M. C., Hebling J. & Nelson K. A. Generation of 10μJ ultrashort terahertz pulses by optical rectification. Appl. Phys. Lett. 90, 171121, doi: (2007).10.1063/1.2734374 [Cross Ref]
  • Stepanov A. G., Bonacina L., Chekalin S. V. & Wolf J.-P. Generation of 30 μJ single-cycle terahertz pulses at 100 Hz repetition rate by optical rectification. Opt. Lett. 33, 2497–2499, doi: (2008).10.1364/OL.33.002497 [PubMed] [Cross Ref]
  • Hirori H., doi A., Blanchard F. & Tanaka K. Single-cycle terahertz pulses with amplitudes exceeding 1 MV/cm generated by optical rectification in LiNbO3. Appl. Phys. Lett. 98, 091106 doi: (2011).10.1063/1.3560062 [Cross Ref]
  • Zeppenfeld M., Motsch M., Pinkse P. W. H. & Rempe G. Optoelectrical cooling of polar molecules. Phys. Rev. A 80, 041401, doi: (2009).10.1103/PhysRevA.80.041401 [Cross Ref]
  • Zeppenfeld M. et al. . Sisyphus cooling of electrically trapped polyatomic molecules. Nature 491, 570–573, doi: (2012).10.1038/nature11595 [PubMed] [Cross Ref]
  • Prehn A., Ibrügger M., Glöckner R., Rempe G. & Zeppenfeld M. Optoelectrical cooling of polar molecules to submillikelvin temperatures. Phys. Rev. Lett. 116, 063005, doi: (2016).10.1103/PhysRevLett.116.063005 [PubMed] [Cross Ref]
  • Yachmenev A., Yurchenko S. N., Jensen P. & Thiel W. A new “spectroscopic” potential energy surface for formaldehyde in its ground electronic state. J. Chem. Phys. 134, 244307, doi: (2011).10.1063/1.3599927 [PubMed] [Cross Ref]
  • Al-Refaie A. F., Yachmenev A., Tennyson J. & Yurchenko S. N. ExoMol line lists VIII: A variationally computed line list for hot formaldehyde. Mon. Not. R. Astron. Soc. 448, 1704–1714, doi: (2015).10.1093/mnras/stv091 [Cross Ref]
  • Owens A., Yurchenko S. N., Yachmenev A., Tennyson J. & Thiel W. Accurate ab initio vibrational energies of methyl chloride. J. Chem. Phys. 142, 244306, doi: (2015).10.1063/1.4922890 [PubMed] [Cross Ref]
  • Owens A., Yurchenko S. N., Yachmenev A., Tennyson J. & Thiel W. A global ab initio dipole moment surface for methyl chloride. J. Quant. Spectrosc. Radiat. Transf. 184, 100, doi: (2016).10.1016/j.jqsrt.2016.06.037 [Cross Ref]
  • Sidje R. B. Expokit: A software package for computing matrix exponentials. ACM Trans. Math. Softw. 24, 130–156, doi: (1998).10.1145/285861.285868 [Cross Ref]
  • Feit M. D., Fleck J. A. & Steiger A. Solution of the schrödinger equation by a spectral method. J. Comput. Phys. 47, 412–433, doi: (1982).10.1016/0021-9991(82)90091-2 [Cross Ref]
  • Bandrauk A. D. & Shen H. Exponential split operator methods for solving coupled time-dependent schrödinger equations. J. Chem. Phys. 99, 1185–1193, doi: (1993).10.1063/1.465362 [Cross Ref]
  • Bunker P. R. & Jensen P. Molecular Symmetry and Spectroscopy 2nd edn (NRC Research Press, Ottawa, 1998).
  • Yurchenko S. N., Thiel W. & Jensen P. Theoretical ROVibrational Energies (TROVE): A robust numerical approach to the calculation of rovibrational energies for polyatomic molecules. J. Mol. Spectrosc. 245, 126–140, doi: (2007).10.1016/j.jms.2007.07.009 [Cross Ref]
  • Yachmenev A. & Yurchenko S. N. Automatic differentiation method for numerical construction of the rotational-vibrational Hamiltonian as a power series in the curvilinear internal coordinates using the Eckart frame. J. Chem. Phys. 143, 014105, doi: (2015).10.1063/1.4923039 [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group