Home | About | Journals | Submit | Contact Us | Français |

**|**Scientific Reports**|**PMC5364483

Formats

Article sections

Authors

Related links

Sci Rep. 2017; 7: 45068.

Published online 2017 March 24. doi: 10.1038/srep45068

PMCID: PMC5364483

Alec Owens,^{a,}^{1} Emil J. Zak,^{1} Katy L. Chubb,^{1} Sergei N. Yurchenko,^{1} Jonathan Tennyson,^{1} and Andrey Yachmenev^{b,}^{2}

Received 2016 October 21; Accepted 2017 February 16.

Copyright © 2017, The Author(s)

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

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 PH_{3} and NH_{3}, and spontaneous emission data for optoelectrical Sisyphus cooling of H_{2}CO and CH_{3}Cl 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 molecules^{1}. 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 systems^{2} along with the increasing precision of instruments has created a demand for accurate theoretical approaches. At present, with a few exceptions^{3}^{,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) method^{7}^{,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 database^{10}^{,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 databases^{10}^{,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 methods^{18}. 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 database^{10}^{,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 calculations^{19}. Originally set up to aid the characterisation and modelling of cool stars and exoplanets^{10}, the data has found other uses beyond atmospheric applications. These include the calculation of molecular partition functions^{20}, radiative lifetimes of ro-vibrational states^{21}, 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 ratio^{23}^{,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 PH_{3} and NH_{3}, and spontaneous emission data for Sisyphus cooling of H_{2}CO and CH_{3}Cl. 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 CH_{3}Cl, 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 spectroscopy^{25}, frequency chirp excitation techniques^{26}, alignment and orientation using static and laser fields^{27}^{,28}, molecular trapping and cooling^{29}.

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 pulses^{30}^{,31}^{,32}, linearly polarized pulses with 45°-skewed mutual polarization^{33}^{,34}, an optical centrifuge^{35}, and intense THz pulses^{36}^{,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 signals^{36}. 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 PH_{3} and NH_{3}. 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 where the polarization axis *A*=*X, Y, Z*, were extracted from the room temperature line lists for PH_{3}^{42} and NH_{3}^{43} (available from the ExoMol database; the much larger hot line lists^{43}^{,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 *v*_{i} corresponding to the ground and first excited vibrational states, i.e. *v*_{i}=0, 1. For NH_{3}, due to strong anharmonicity of the large amplitude umbrella-type vibrational motion we extended the vibrational basis set to include *v*_{umb}=0, …, 7.

For both molecules only one quantum wavepacket simulation corresponding to the initial ground ro-vibrational state was performed. For PH_{3} 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 NH_{3}, 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*=0K and *T*=10K for PH_{3}. At *T*=10K 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 function^{39}

where *E*_{0} is the amplitude of the pulse, is the pulse duration, and is the central frequency. In all calculations we fixed to 4ps. For simulations of the mixed Raman-THz field-free orientation of PH_{3}, 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 TROVE^{6} 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.

The results of our numerical simulations of PH_{3} 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 ps^{−1} and peak field strengths at *E*_{0}=250kV/cm and *E*_{0}=1 MV/cm have been used to induce transient alignment cos^{2}*θ* (*t*) and orientation cos*θ* (*t*), where *θ* is the Euler angle. Such intense THz pulses in the few-cycle regime can be routinely generated in a convenient tabletop setup^{45}^{,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 *T*_{rev}=1/2*Bc*≈3.7 ps (*B*≈4.5cm^{−1} is the rotational constant of PH_{3} 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 Δ*τ*=*T*_{rev}.

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 cos^{2}*θ* (*t*) and cos*θ* (*t*). However, raising the temperature from *T*=0K to *T*=10K 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*=10K 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 PH_{3} 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.9ps, which corresponds to *T*_{rev}.

As a second example we explore strategies to establish control over rotational coherences and populations in NH_{3} by applying optimally timed THz pulses. Recently it was shown both theoretically and experimentally^{41} 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 NH_{3}, 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 (see Eq. (1)). The results are shown in Fig. 4 for the three rotational transitions, |0, 1→|1, 0, |1, 0→|2, 1, and |2, 1→|3, 0, where states of NH_{3} are labelled as |*J, v*_{umb}, 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 are 0.4, 1.1, and 1.6ps^{−1} for the three different rotational excitations of NH_{3}.

Figure 5 depicts the populations of the ro-vibrational wave packets of NH_{3}, initially in the |0, 1, |1, 0, and |2, 1 ro-vibrational states, following the excitation by a time-delayed THz pulse pair with *E*_{0}=500kV/cm and respective pulse central frequencies , 1.1, and 1.6ps^{−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 *T*_{rev}≈*J*/2*Bc*≈1.7, 0.8 and 0.6ps, respectively, where *B*≈6.35cm^{−1} is the rotational constant of NH_{3}.

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 . The second pulse is tuned to affect population of the |1, 0 state and transforms the two-state coherence into a new one: . 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.

To further investigate the validity of our approach we have computed the alignment and orientation of NH_{3} 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 , 0.6, and 0.9ps^{−1} each with duration 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 .

As shown in Fig. 7 it is only for the very intense 10MV/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 10MV/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 NH_{3} 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 ps^{−1} and pulse duration 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. 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, , 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.

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 50MV/cm.

Optoelectrical Sisyphus cooling is a robust and general method for producing ultracold polyatomic molecules^{48}^{,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.

Sisyphus cooling has been realised experimentally for methyl fluoride (CH_{3}F)^{49} and formaldehyde (H_{2}CO)^{50}. However, a number of symmetric top molecules with strong parallel vibrational transitions have been proposed as suitable candidates for cooling using this method^{48}. 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 H_{2}CO and CH_{3}Cl. These systems have been treated within the ExoMol framework^{19} and we refer the reader to the relevant publications for details of the ro-vibrational calculations^{51}^{,52}^{,53}^{,54}. Candidate ro-vibrational states were identified by sorting the 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 10^{1} 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.

In the work of ref. ^{50}, Sisyphus cooling of formaldehyde utilized ro-vibrational energy levels of the *v*_{1} symmetric C–H stretching mode, |*v*_{1}, 3, 3, and of the vibrational ground state, |0, 3, 3 and |0, 4, 3. Here energy levels are labelled as |Mode, *J, K* where *K*=|*K*_{A}| and we omit the quantum number *K*_{C} (see discussion in ref. ^{50}). The primary decay channel |*v*_{1}, 3, 3→|0, 3, 3 had a branching ratio of *γ*=0.75 and an estimated spontaneous decay rate (equivalent to the Einstein A coefficient) of approximately 60Hz, 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 |*v*_{1}, 3, 3 level has 41 decay channels, of which only transitions to the |0, 3, 3 and |0, 4, 3 levels were dominant. We predict a slightly smaller Einstein A coefficient for the |*v*_{1}, 3, 3→|0, 3, 3 transition and hence a slightly slower cooling rate than ref. ^{50}. However, the difference is small (≈10Hz) and our computed branching ratios are identical to those of ref. ^{50}.

For other possible candidate states in the *v*_{1} manifold, a slightly higher branching ratio and Einstein A coefficient are predicted for the |*v*_{1}, 4, 4→|0, 4, 4 transition. Alternatively, the *v*_{2} C–O stretching mode appears suitable for Sisyphus cooling but using these states would half the speed of cooling compared to the *v*_{1} 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*=100K 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.

Methyl chloride has not yet been considered for Sisyphus cooling but we expect a similar experiment to CH_{3}F^{49} could be performed. As shown in Table 2, ro-vibrational levels of the *v*_{1} symmetric CH_{3} stretching mode and vibrational ground state appear the best candidates for cooling. Similar to methyl fluoride, which had a spontaneous decay rate of about 15Hz for the |*v*_{1}, 3, 3→|0, 3, 3 transition, the respective rate in CH_{3}Cl is only marginally larger and would hence provide a similar speed of cooling. Note that we have only considered CH_{3}^{35}Cl for the present analysis but we would expect the same conclusions for CH_{3}^{37}Cl.

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*=1500K. 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 NH_{3} and H_{2}O, 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 PH_{3} and NH_{3} 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 H_{2}CO and CH_{3}Cl. In both molecules the *v*_{1} 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.

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 and files, unique to each molecule. The 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 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 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.

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

Here *ε*_{l,J} and |*l, J, m* 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 ( 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 in the ro-vibrational basis and the time-dependent classical electric field *E*_{A}(*t*), defined in the laboratory fixed axes system. To describe the radiative fields we use the form , where *A* is the polarization axis, is the pulse time profile with a maximum at *t*=*t*_{0} and duration (FWHM) *T*, and *ω*_{A} is the carrier frequency.

Calculation of the ro-vibrational matrix elements from Einstein A coefficients, extracted from the ExoMol line list, is described in detail below. The selection rules for 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*, 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 . 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

where the elements of the matrices *H*_{0} and are given by and , 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 package^{55}. Alternative and more sophisticated representations for the split-operator technique can be found in refs ^{56} and ^{57}.

We aim to derive the relationship between the Einstein A coefficients, , and the matrix elements of the electric dipole moment operator, (), defined relative to the laboratory fixed system. It is convenient to first express in terms of the dipole line strength as

where is the *f*←*i* transition frequency in units of cm^{−1}, *g*_{ns} is the nuclear spin statistical weight factor, *J*′ refers to the rotational quantum number of the final state *f*, and the units of *A*_{fi} and *S*_{fi} are *s*^{−1} and *D*^{2}, respectively.

Writing the wavefunctions for the initial and final states as products of rotational and vibrational basis functions, i.e. and , the expression for *S*_{fi} reads^{58}

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

It should be noted that the actual wavefuncitons |*ψ* are expansions in terms of |*J, k, m*|*v* 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 may be expressed as

By combining Eq. (7) with Eq. (4) and Eq. (5), it is straightforward to write down in terms of the quantities *S*_{fi} and *A*_{fi} 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 *A*_{fi}.

For the sake of simplicity the above formulas were derived using the symmetric top basis functions |*J, k, m*, 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 TROVE^{59}^{,60} and RichMol program packages, we employ the basis of symmetrized combinations of the rigid rotor functions for *k*>0 and |*J*, 0, 0, *m* 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 with *P*_{fi}=0 or 1, the final expression for the dipole matrix elements reads as

**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.

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.

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.

- 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 NH
_{3}. 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 NH
_{3}. 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 100Hz 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 LiNbO
_{3}. 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**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |