Quantum-state-dependent decay rates of electrostatically trapped Rydberg NO molecules

Nitric oxide (NO) molecules travelling in pulsed supersonic beams have been prepared in long-lived Rydberg–Stark states by resonance-enhanced two-colour two-photon excitation from the X P1/2 (v00 = 0, J00 = 3/2) ground state, through the A S (v0 = 0, N0 = 0, J0 = 1/2) intermediate state. These excited molecules were decelerated from 795 ms 1 to rest in the laboratory-fixed frame of reference, in the travelling electric traps of a transmission-line Rydberg–Stark decelerator. The decelerator was operated at 30 K to minimise effects of blackbody radiation on the molecules during deceleration and trapping. The molecules were electrostatically trapped for times of up to 1 ms, and detected in situ by pulsed electric field ionisation. Measurements of the rate of decay from the trap were performed for states with principal quantum numbers between n = 32 and 50, in Rydberg series converging to the N= 0, 1, and 2 rotational states of NO. For the range of Rydberg states studied, the measured decay times of between 200 ms and 400 ms were generally observed to reduce as the value of n was increased. For some particular values of n deviations from this trend were seen. These observations are interpreted, with the aid of numerical calculations, to arise as a result of contributions to the decay rates, on the order of 1 kHz, from rotational and vibrational channel interactions. These results shed new light on the role of weak intramolecular interactions on the slow decay of long-lived Rydberg states in NO.


Introduction
Rydberg states of atoms and small molecules, including H 2 , N 2 , O 2 , and NO, are of importance in recombination processes that occur in laboratory, atmospheric, and astrophysical plasmas. [1][2][3][4][5][6] They are also of interest in studies of chemical dynamics at low temperatures. For example, they can be employed to study ionmolecule reactions, as most recently demonstrated in the case of the H 2 + + HD -H 2 D + + H reaction, at temperatures as low as 50 mK. [7][8][9][10] In this setting the Rydberg electron shields the reaction centre from stray electric fields that would otherwise cause heating. Long-lived Rydberg states of small molecules also offer opportunities for studies of low-energy molecular scattering that are dominated at long range by resonant dipoledipole interactions, 11 and the exploration of the chemistry associated with ultra-long range Rydberg molecules 12 formed of a Rydberg atom or molecule bound to a ground-state atom or molecule through low-energy Rydberg-electron scattering. 13,14 They are also important in tests of molecular quantum mechanics and quantum chemistry, for example, through the precise determination of ionisation and dissociation energies. [15][16][17] The implementation of general methods for preparing samples of cold molecules in high Rydberg states is of interest in each of these areas, and is also foreseen as a means of producing cold samples of ground-state species following de-excitation. 18 For each of these endeavours a detailed understanding of the lifetimes and decay processes of long-lived Rydberg states is essential. The large electric dipole moments associated with Rydberg-Stark states of atoms and molecules mean that forces can be exerted on them using inhomogeneous electric fields. 19,20 This has led to the development of the methods of Rydberg-Stark deceleration [21][22][23] and the implementation of guides, [24][25][26][27] beamsplitters, 28 decelerators, [29][30][31] and traps 18,[31][32][33][34][35][36] for atoms and molecules. To implement these experimental techniques it is necessary to prepare Rydberg states with sufficiently long lifetimes, i.e., 410 ms, and small avoided crossing between individual sublevels in the presence of electric fields, that can be traversed diabatically on typical experimental timescales. Both of these requirements can be met in general for small molecules through the preparation of Rydberg states with values of |M N | Z 3 (M N is the projection of the total angular momentum vector excluding spin, -N, onto the z-axis in the laboratory frame of reference). 18 In previous experiments with H 2 , the preparation of selected |M N | = 3 Rydberg-Stark states was achieved by threephoton excitation from the X 1 S g + (v = 0, J = 0) ground electronic state using circular polarised vacuum ultraviolet, visible, and infrared laser radiation. 18,37 When these molecules were confined in a cryogenically cooled trap, n-dependent (n is the principal quantum number) decay times on the order of 100 ms, and effects of rotational channel interactions, were observed. 38 In the work described here, |M N | Z 3 Rydberg states were prepared in NO by two-photon laser photoexcitation from the ground electronic state with subsequent M N -mixing arising as a result of intramolecular interactions and effects of static and time-varying electric fields close to the time of photoexcitation. 39 Deceleration and trapping of the NO molecules in these long-lived states has permitted studies of slow decay processes over timescales of up to 1 ms that were previously inaccessible.
Studies of the lifetimes, and decay dynamics of high-n, n = 40 to 122, v + = 0 (v + is the vibrational quantum number of the NO + ion core) Rydberg states in NO have been reported previously by Vrakking and Lee. 40,41 This work was performed to aid in the interpretation of excited-state lifetimes observed in ZEKE spectroscopy experiments. 42 These measurements were performed on timescales of up to B1 ms. The results showed that np(0) and nf (2) [nc(N + ), where c is the orbital angular momentum quantum number of the Rydberg electron, and N + is the rotational quantum number of NO + ] Rydberg states in NO have short lifetimes, on the order of B1 ns and B10 ns, respectively, for this range of values of n in zero-electric-field. This indicated that fast predissociation dominated the decay of these low-c states. However, when photoexcitation was performed in the presence of sufficiently large electric fields, lifetimes on the order of 100 ns were observed. Calculations reported by Bixon and Jortner 43 showed that this lifetime enhancement could be accounted for by considering a combination of electric-field-induced c-mixing and intramolecular interactions, between each optically accessible short-lived low-c state and long-lived high-c states.
Measurements and calculations reported by Goodgame et al. 44 on the Stark effect in v + = 1 Rydberg states in NO with values of n between 10 and 20, located above the v + = 0 series limit, allowed for the study of intramolecular interactions, and of the lifetimes of the states excited to be estimated from the measured spectral linewidths. For Rydberg states in series converging to excited vibrational states of NO + , non-radiative decay by vibrational autoionisation can occur if the interaction between the Rydberg electron and the ion core couples the Rydberg state with the continuum of a lower vibrational series. Rydberg states in the v + = 1 series, that lie above the v + = 0 series limit, can therefore undergo vibrational autoionisation because of coupling to the v + = 0 continuum. From a comparison of calculated and measured spectral linewidths, lifetimes of the low-c states with v + = 1 and values of n between 10 and 20 were estimated to be between B1 ps and B50 ps. These lifetimes are dominated by fast non-radiative decay and contain contributions from both predissociation and autoionisation.
Further studies of Rydberg states in NO in electric fields were carried out by Jones et al. 45 and Patel et al. 46 to investigate effects of intramolecular interactions on the Stark structure of v + = 0 Rydberg states. These works included laser photoexcitation spectra recorded in electric fields of between 0 and 150 V cm À1 , with detection by pulsed electric field ionisation, and by monitoring predissociated fragments by photoionisation. These measurements, combined with the results of numerical calculations, allowed intramolecular couplings of optically accessible np(N + ) and nf(N + ) states with neighbouring states to be studied. They also allowed the distribution of np(N + ) and nf(N + ) character among the c-mixed Stark states in the presence of electric fields to be determined. The implementation of two complementary detection methods in these works yielded information on the importance of predissociation in the decay of the states excited, and how this changed in an electric field.
Recently, using a two-dimensional spectroscopy technique based on laser photoexcitation and delayed state-selective electric field ionisation, long-lived Rydberg states in NO were characterised in detail for values of n in the range from 40 to 100. 39 States with lifetimes in excess of 10 ms were identified and observed to ionise diabatically. These results indicated that high-|M N | c-mixed states with long lifetimes, i.e., states without any short-lived c r 3 character, were populated upon photoexcitation. That work led to the first experiments to decelerate and electrostatically trap Rydberg NO molecules. 47 The experiments and calculations reported here build upon these works, and include new high-precision decay-rate measurements for long-lived Rydberg states of different rotational series, and the identification of contributions from weak intramolecular interactions to the measured decay rates.
In the following, the experimental apparatus and techniques are described in Section 2. In Section 3, calculations of the energy-level structure and lifetimes of the Rydberg states in NO in electric fields, and considerations necessary for the implementation of Rydberg-Stark deceleration are discussed. In Section 4, the results of trap decay measurements performed for hydrogenic Rydberg-Stark states with values of n between 32 and 50 are presented and interpreted in terms of the effects of intramolecular rotational and vibrational channel interactions. Finally conclusions are drawn in Section 5.

Experiment
A schematic diagram of the apparatus used in the experiments is shown in Fig. 1. A pulsed supersonic beam of NO was generated with a mean longitudinal speed of 810 ms À1 . The pulsed valve from which the beam emanated was operated at room temperature at a repetition rate of 25 Hz. After passing through a 2-mm-diameter skimmer, the beam entered the cryogenically cooled region of the apparatus containing a chip-based transmission-line Rydberg-Stark decelerator. 30,33 Laser photoexcitation of the NO molecules to Rydberg states was performed between electrodes E1 and E2, as shown in Fig. 1(a), using the ncX + 1 S + (v + = 0, N + ) ' A 2 S + (v 0 = 0, N 0 = 0, J 0 = 1/2) ' X 2 P 1/2 (v 00 = 0, J 00 = 3/2) resonance-enhanced twocolour two-photon excitation scheme 39,46,48,49 depicted in Fig. 2. Counter-propagating radiation from a pair of Nd:YAG pumped pulsed dye lasers, with full-width-at-half-maximum (FWHM) spectral widths of 0.17 cm À1 , was used for photoexcitation. The frequency-tripled output of the first laser was set to u 1 = 44193.99 cm À1 (226.275 nm), a pulse energy of B10 mJ, and linearly polarised in the z dimension to drive the A 2 S + ' X 2 P 1/2 transition. The frequency-doubled output of the second laser was tuned over the wavenumber range from u 2 = 30 420 cm À1 to 30 500 cm À1 (328.7 nm to 327.8 nm), had a pulse energy of B1.6 mJ and was linearly polarised in the y dimension to access np(0) and nf(2) Rydberg states from the intermediate A 2 S + state. The fundamental wave numbers of the dye lasers were monitored using a fibre-coupled wavelength meter with an absolute accuracy of AE0.005 cm À1 . At the time of photoexcitation, gated time-dependent electric potentials could be applied to electrode E1 to enhance population transfer to long-lived c-mixed hydrogenic Rydberg-Stark states.
After photoexcitation the excited Rydberg molecules travelled for 5.3 ms (B4 mm) before being loaded into the travelling electric traps of a chip-based transmission-line Rydberg-Stark decelerator. 47 The decelerator was operated by applying a set of 5 sinusoidally oscillating electric potentials of amplitude V 0 to the segmented centre-conductor electrodes (see Fig. 1), and a constant offset potential of À V 0 /2 to electrode E2. Molecules confined in the travelling traps were decelerated from 795 ms À1 to zero mean velocity in the laboratory-fixed frame of reference, over a distance of B100 mm, and in a time of B250 ms. In each cycle of the experiment the bunch of molecules excited was loaded into a single travelling trap. The curvature of the decelerator in the xz-plane minimised effects of collisions of decelerated and trapped molecules with the trailing component of the molecular beam. 34 To minimise blackbody radiation induced transitions during deceleration and trapping the decelerator was enclosed in a heat shield and cooled to 30 K.
The trapped molecules were detected in situ within the decelerator structure by pulsed electric field ionisation (PFI) upon the application of a pulsed potential of +500 V to the side electrodes [see Fig. 1(a)]. The resulting NO + ions were accelerated through the 2-mm-diameter aperture A1 in electrode E2, and collected on a microchannel plate (MCP) detector. Measurements were performed with rapidly switched PFI potentials (rise time B50 ns), to detect all Rydberg molecules present in the detection region above A1. To obtain information on the populations of the Rydberg states at the time of ionisation, slowly-rising pulsed potentials could be applied. These had a time-dependence of the form V Ramp (t) = ÀV max [1 À exp(Àt/t)], with the time constant t = 4.4 ms. When this slowly rising pulsed potential was applied, ionised electrons were accelerated through aperture A1 to the MCP. The shorter flight time of the electrons to the MCP detector (typically o10 ns), than the NO + ions (typically B1 ms), permitted the most direct mapping of the detection time of the electrons to the potential applied at the time of ionisation. 50 The alignment of the apparatus was optimised by performing time-of-flight (TOF) measurements with the decelerator off,  Resonance-enhanced two-colour two-photon excitation scheme employed in the preparation of Rydberg states in NO converging to the lowest vibrational state (v + = 0) of the NO + ion core. Low-n states in the v + = 1 series that are located close to the v + = 0, N + = 0 series limit are also shown. The energies of the levels in this figure were either measured in the experiments, 39 or calculated using the methods described in Section 3.1. and PFI above the on-axis aperture A2 in electrode E2. Decelerated and trapped molecules were detected by PFI above the off-axis aperture A1. PFI at times beyond those for which the trap in the decelerator in which the molecules were confined was brought to rest, allowed studies of the rate of decay of the molecules.

Stark effect in Rydberg states of NO
To aid the interpretation of the results of the experiments, calculations of the energy-level structure and lifetimes of the Rydberg states in NO in the presence of externally applied electric fields were performed. The methods used in these calculations follow those reported in ref. 44, 46, 51 and 52, and involved determining the eigenvalues and eigenvectors of the Hamiltonian matrix describing the interaction of the molecule with an external electric field.
Although the J 00 = 1/2 rotational level of the X 2 P 1/2 ground electronic state in NO can be described purely in Hund's angular momentum coupling case (a), excited rotational levels, including the X 2 P 1/2 (v 00 = 0, J 00 = 3/2) state used in the laser photoexcitation scheme implemented here, are intermediate cases between Hund's case (a) and Hund's case (b). 53,54 On the other hand, the high Rydberg states of NO, prepared here for deceleration and electrostatic trapping, are most readily described in Hund's case (d). 44 Calculations of the Stark effect in these Rydberg states were therefore performed in the |ncN + NM N i Hund's-case-(d) basis where n and c are the principal and orbital angular momentum quantum numbers of the Rydberg electron, respectively, N + is the rotational quantum number of the NO + ion core, N is the total angular momentum quantum number excluding spinÑ ¼Ñ þ þ' , and M N is the projection of -N onto the z-axis in the laboratory frame of reference. Here this z-axis is defined by the local external electric field vector at the position of the molecule. In this basis, the total Hamiltonian can be expressed as, where H (d) 0 and H (d) Stark are the field-free and Stark Hamiltonians, respectively. H (d) Stark contains all contributions arising from the presence of the external electric field.

Field-free Hamiltonian
For Rydberg states in NO, the field-free Hamiltonian, H (d) 0 in eqn (1), contains diagonal contributions that represent the energies of the unperturbed Hund's-case-(d) eigenstates. These were calculated using the Rydberg formula. Off-diagonal contributions to H (d) 0 arise from the interactions of the Rydberg electron with the electric dipole moment, electric quadrupole moment, and electric polarisability of the NO + ion core. [55][56][57] They also contain contributions from the multi-electron character of NO + which causes ss-ds configuration mixing 44,52 (s represents the l = 0 projection of the Rydberg electron orbital angular momentum vector onto the internuclear axis).
The general approach used here to construct the matrix H (d) 0 follows that presented by Goodgame et al., 44 which itself is based on the work of Fredin et al. 52 on the determination of off-diagonal elements of the field-free Hamiltonian from a first order approximation to multichannel quantum defect theory (MQDT). Off diagonal matrix elements arising from the interaction of the Rydberg electron with the electric dipole moment of the NO + ion were calculated using the method described by Bixon and Jortner. 43 To construct the Hamiltonian, H (d) 0 , it was separated into a diagonal part, H (d) diag , a part that accounts for the interaction of the Rydberg electron with the static electric dipole moment of the NO + ion core, H (d) dip , and a part that accounts for the combination of the interaction of the Rydberg electron with the electric quadrupole moment and higher order multipole moments of NO + , and the multi-electron character of the molecule, H (d) multi , i.e., The Hamiltonian H (d) diag is diagonal in the |ncN + NM N i basis with matrix elements given by, Here, W v + is the ionisation wavenumber of the Rydberg series converging to the ground electronic state of the NO + cation with vibrational quantum number v + . B þ v þ and D þ v þ are the rotational and centrifugal distortion constants for NO + . The Rydberg constant adjusted for the reduced mass of NO is R NO = 109735.31 cm À1 , and n ¼ n À m For v + = 0, the rotational constant is B + 0 = 1.987825 cm À1 , and the centrifugal distortion constant is D + 0 = 5.64 Â 10 À6 cm À1 . 46 The elements of the Hamiltonian matrix H (d) dip are given by, 43,56,57 where e is the electron charge, e 0 is the vacuum permittivity, m NO + is the electric dipole moment of NO + in a coordinate system with its origin at the centre of mass of the molecule, and hnc|r À2 |n 0 c 0 i is a radial integral. The angular integral, where the terms in the ( ) and { } parentheses represent Wigner 3J and 6J symbols, respectively, and d N,N 0 and d M N ;M 0 N are Kronecker delta functions. H (d) dip leads to interactions between Hund's-case-(d) basis states for which, The radial integral, in eqn (4), was calculated using the Gilbert-Child 58 near-threshold approximation with a finite quantum defect such that, where a NO is the Bohr radius adjusted for the reduced mass of Þ . This charge-dipole interaction is strongest between low-c states, i.e., states with non-zero quantum defects. In these calculations, a static electric dipole moment of NO + of m NO + C 0.4 D was used. This was based on the value reported from calculations, at the equilibrium internuclear separation, for the X + 1 S + state of the NO + ion. [59][60][61][62] The elements of the Hamiltonian matrix H (d) multi can be expressed as, 44 where A ðN'Þ N þ L represents the Hund's-case-(b) to Hund's-case-(d) frame transformation, and H L '' 0 are Hund's-case-(b) matrix elements.
For Rydberg states converging to the 1 S + state of NO + with L + = 0 (L + is the projection of the total electron orbital angular momentum vector of the NO + cation onto the internuclear axis), the projection of the total orbital angular momentum vector in the neutral NO molecule onto the internuclear axis is Note that the Rydberg-electron spin and spin projection are common to Hund's-case-(b) and Hund's-case-(d) so are not considered in these transformations. 63 In NO, the configuration mixing of importance in the work described here occurs predominantly between the ss (c = 0, L = 0) and ds (c = 2, L = 0) Hund's-case-(b) basis states, and therefore results in c-mixing of the Rydberg eigenstates.
In Hund's-case-(b), the corresponding matrix elements, that account for these effects of configuration mixing, are, 44,52 where y sd = À38.71 is the ss-ds mixing angle, 52 and m (b) cL are the Hund's-case-(b) quantum defects. Taking this into account, the complete sets of selection rules for the intramolecular interactions accounted for in H (d) multi are, where the Dc = AE2 interactions only occur between basis states with c = 0 and c 0 = 2, and vice versa. The Hund's-case-(b) quantum defects, m (b) cL , used in the calculations are listed in Table 1. To obtain the Hund's-case-(d) quantum defects,

Stark Hamiltonian
In an electric field -F = (0, 0, F z ), the Hund's-case-(d) Hamiltonian in eqn (2) can be extended to include contributions from the Stark effect. The resulting perturbation, H (d) Stark = eF z z, has matrix elements given by, 37,65 where the integral, This perturbation leads to non-zero off diagonal couplings between states, within the same rotational series, for which Dc = AE1, i.e., between states for which, The radial integrals hnc|r|n 0 c 0 i on the right-hand side of eqn (15) were evaluated using the Numerov method 66 with a pure Coulomb potential. Stepwise integration was performed inwards from large r to either the inner classical turning point of the Rydberg electron orbit, or the polarisability radius, r a , of the NO + ion core whichever was encountered first. The polarisability radius, r a , was taken to be r a C (a NO +) 1/3 , 66 where a NO + C 8 a.u. is the electric dipole polarisability of the NO + cation. 67,68

Lifetime calculations
Using the eigenvectors obtained following diagonalisation of H (d) tot , the lifetimes of the Rydberg states in NO in the presence of an electric field could be determined. From the rates of decay by fluorescence, g fl(k) , and dissociation, g dis(k) , of each basis state, denoted |ki = |n k c k N + k N k M N k i, the total decay rate of an eigenstate |ii of H (d) tot is, where c ik (F z ) = hk|ii (Fz) are the coefficients of the corresponding eigenvector in the electric field F z . The fluorescence rates g fl(k) = 1/t fl(k) of the basis states were calculated by considering the zero-field lifetimes, t nc , of the hydrogenic Rydberg states averaged over all polarisations. The effect of c-mixing was then accounted for through the dependence of c ik (F z ) on the field strength. The fluorescence lifetimes can be expressed as, where A nc,n 0 c 0 are the Einstein A coefficients, and the primes denote states energetically below the |n, ci state. The values of A nc,n 0 c 0 are given by, 69 with o nc,n 0 c 0 = (E nc À E n 0 c 0 )/h , and h the reduced Planck constant. The dissociation rates, g dis(k) , required in eqn (17) are dependent upon n k and c k . This dependence can be expressed as, 43 where R dis (c k ) is the dissociation rate parameter for a basis state in the v + = 0 series with c = c k . The dissociation rate parameters for Rydberg states in the v + = 0 series that are available in the literature for values of c r 4, and were used in the calculations described here, are presented in Table 2. For higher values of c, R dis (c) was assumed to be zero.
The calculated lifetimes, t 0(i) = 1/g 0(i) , of the n = 38, N + = 2 Rydberg states in the v + = 0 series with M N = 0, 2, and 4 are shown in Fig. 3(a), (b), and (c), respectively. In each panel, the 2N + + 1 = 5 Stark manifolds can be seen. The results of all calculations shown in this figure were obtained for an electric field of 5 V cm À1 . In this field sufficient c-mixing has occurred for the calculated lifetimes to provide a reasonable representation of those of the same states in higher fields, while Stark states with different values of n remain clearly distinguishable allowing for ease of interpretation of the results. For low values of |M N |, i.e., |M N | r 2, dissociation represents the dominant contribution to the decay of the excited states and the lifetimes of the states in the different Stark manifolds overlap. For larger values of |M N |, dissociation plays a less important role, and the manifolds of states decay predominantly by fluorescence and therefore have distinguishable lifetimes.
The high dissociation decay rates of the low-c (c r 3) basis states, which are between B10 8 s À1 and B5 Â 10 9 s À1 for n = 38, dominate over the fluorescence decay rates (B10 4 s À1 for n = 38) of these states. On the other hand, for high-c basis states (c Z 4) the total decay rates are smaller, with fluorescence rates ranging from B10 2 s À1 to B10 4 s À1 and dissociation rates o3 Â 10 6 s À1 . Because of the significant reduction in the dissociation rates as the value of c increases, only the dissociation of basis states with c r 4 were considered in calculations. With this in mind, mixed states, with some high-c character and some low-c character, can have significantly reduced lifetimes compared to those of the higher-c states in zero-electric-field. This can be seen in Fig. 3(a), where all possible c values are present for M N = 0, with significant lowc character mixed into all Stark states. The lifetimes of states with M N = 0 are all o1 ms, i.e., they decay predominately by dissociation even in weak electric fields.
When |M N | 4 0 not all low-c basis states contribute to the c-mixed Stark states. However individual low-c states are not absent from all Stark manifolds with the same value of |M N |. Contributions from basis states with particular values of N + and c reduce when |M N | 4 |N + À c|, and are completely absent when |M N | 4 N + + c. Therefore, in the example case with N + = 2, to prevent c r 3 basis states mixing into at least one of the Stark manifolds populated, it is necessary for |M N | to be greater than or equal to 2. The effect this has on the lifetimes can be seen in Fig. 3(b and c). When M N = 2 [ Fig. 3(b)], only one of the N + = 2 Stark manifolds has no c r 3 character. This manifold is easily identifiable as that with lifetimes between B10 ms  Fig. 3(c)], the Stark manifolds are clearly distinguishable by the lifetimes of their component states. In the manifold with the shortest lifetimes (B10 ns) the lowest c basis states present are those for which c = 2. Whereas, in the manifold with the longest lifetimes (B1 ms) c Z 6. In between these extreme cases the lowest c basis state present lies between these values. In these data it is therefore possible to distinguish between manifolds of states that decay predominantly by dissociation and contain low-c (c r 3) character, and those with only high-c (c Z 4) character for which decay by fluorescence dominates. The asymmetry in the lifetimes of the states of the Stark manifolds with c r 3 character in Fig. 3, arises because the low-c states are not degenerate with those of higher c in zero-field.

Considerations for the implementation of Rydberg-Stark deceleration
An example of a Stark map calculated by diagonalising H (d) tot , and encompassing the wavenumber region around the fieldfree 38(2) state, is displayed in Fig. 4. The energy-level structure for states with M N = 0, 2, and 4 is displayed in Fig. 4 For Rydberg-Stark deceleration and electrostatic trapping, 18,71 c-mixed Stark states with positive Stark energy shifts, i.e., lowfield-seeking (LFS) states, are required. These states must have lifetimes greater than B10 ms, and evolve diabatically through any avoided crossings in the Stark map in the time-dependent electric fields used for deceleration and trapping. 18 These requirements demand that contributions from short-lived, non-degenerate low-c basis states are minimised. This can be achieved, as seen, e.g., from the data in Fig. 3, by populating states with high values of |M N |, i.e., |M N | Z 3. 18 The effect of increasing the value of |M N |, on the Stark energy level structure of the Rydberg states, can be seen in Rydberg states with |M N | Z 3 can be selectively populated by resonance-enhanced multi-photon excitation with circularly polarised laser radiation as achieved, e.g., in the case of H 2 . 18 They can also be populated by M N -mixing in collisions and in the presence of time-varying electric fields. 72,73 In this latter case, deceleration and electrostatic trapping acts to filter the long-lived high-|M N | components of an ensemble of excited molecules. This filtering of long-lived high-|M N | states occurs because of the finite time required for deceleration, which is only effective for molecules in long-lived Rydberg states, and the requirement to undergo diabatic traversal of avoided crossings in the deceleration fields, that only occurs with high probability for high values of |M N |.
The two-colour two-photon laser photoexcitation scheme used in the work reported here led to the population of longlived hydrogenic Rydberg-Stark states in NO through collisions and the effects of time-varying electric fields, as demonstrated previously by state-selective electric field ionisation measurements. 39 Although direct measurements of the lifetimes of the states excited were not made previously, the geometry of the apparatus meant that lifetimes in excess of 10 ms were required for the excited molecules to be detected. This observation, in conjunction with the calculated lifetimes in Fig. 3, suggest that hydrogenic Rydberg-Stark states without any low-c (c r 3) character, i.e., states with |M N | Z 3, were therefore populated.

Results
The methods of Rydberg-Stark deceleration and electrostatic trapping of NO molecules in long-lived Rydberg-Stark states that were employed in the work reported here, have allowed detailed studies of the decay of these excited states on time scales up to 1 ms. Experiments were performed to investigate:

Trap decay measurements
To measure rates of decay of electrostatically trapped Rydberg NO molecules, laser photoexcitation and Rydberg-Stark deceleration were performed as described in Section 2. When the traps in which the molecules were confined were brought to rest, at time t trap = 0 ms in the following, in situ detection by PFI was performed for times up to t trap = 1 ms. Examples of the results of such trap decay measurements, for molecules photoexcited on the 38f(2) resonance (u 2 = 30458.37 cm À1 ), can be seen in Fig. 5. There are two main components to the decay dynamics observed in these data: (i) at early trapping times a fast reduction in signal is observed, reflecting the decay of the Rydberg states predominantly populated upon photoexcitation, and (ii) at later times this reduction slows down reflecting the presence of some molecules in states with longer lifetimes that remain in the trap at late times. 47,74 The data indicated by the open squares in Fig. 5 represent measurements made upon laser photoexcitation in nominally zero electric field, V Mod = 0 V. However, because the mechanism by which long-lived Rydberg-Stark states are populated close to the time of laser photoexcitation depends on M N -mixing that occurs in time-varying electric fields, this can be enhanced by applying a rapidly-varying field in the excitation region when the laser radiation is present. The effect of this sinusoidal field modulation, generated by applying a gated sinusoidally oscillating electric potential to electrode E1 (see Fig. 1) with an amplitude of V Mod = 7.5 V and a frequency of 50 MHz for a time of 50 ns, can be seen by comparing the data recorded with (blue circles) and without (orange open squares) this modulation. The rapid modulation of the field in the excitation region resulted in a factor of B30 increase in the NO + signal detected by PFI, compared to photoexcitation performed with V Mod = 0 V, and hence an increase in the number of molecules in long-lived Rydberg-Stark states that were trapped. The corresponding increase in signal meant that upon the application of this field modulation, molecules could be observed in the trap for times up to 1 ms. Without this field modulation the NO + signal reduced to the level of the background noise after a trapping time of B500 ms. The introduction of this field modulation does not change the measured decay time constants from the trap at early times. An exponential function fitted to both sets of data in Fig. 5, for trapping times between 0 and 500 ms, resulted in decay time constants of 366 AE 6 ms and 377 AE 21 ms for the cases with and without the modulation, respectively. From this it is concluded that rapid electric field modulation close to the time of laser photoexcitation leads to the population of Rydberg-Stark states with similar characteristics to those populated as a result of effects of electric field noise, and the presence of ions. In the following, all measurements were therefore performed with the 50 MHz sinusoidal potential, of amplitude V Mod = 7.5 V, applied at the time of photoexcitation.
To measure quantum-state-dependent decay rates of the trapped molecules single exponential functions were fit to the experimental data, such as that in Fig. 5, over a range of early trapping times to obtain the decay time constant, t decay .
To validate the fitting procedure, and minimise spurious effects arising from the motion of the molecules in the trap, the range of trapping times over which the fit was performed was tested to ensure the stability of values of t decay obtained. This was done for molecules excited on the 38f(2) resonance shown in Fig. 6(a). In Fig. 6(b) the variation in the values of t decay obtained for fixed fit start times of 50 ms or 100 ms, and a range of fit end times can be seen. For end times up to B350 ms the value of t decay remains constant within the the experimental uncertainties. However at later end times the value of t decay gradually increases as the experimental data deviates from a single exponential function. This change in the decay time constant at later trapping times reflects the fact that a range of Rydberg-Stark states were populated close to the time of photoexcitation. 74 Molecules in shorter lived states decay more quickly from the trap, and therefore at later times only molecules in states with longer lifetimes remain. To minimise the effects of these decay dynamics on the values of t decay , the fit end time was chosen to be 350 ms for all measurements described below.    6(c) shows the values of t decay obtained for a range of fit start times. For trapping times r50 ms the value of t decay varied with the fit start time. This behaviour is a consequence of the motion of the molecules in the trap after it is brought to rest. 32,47 During deceleration the molecules are predominantly located ahead of the electric field minima of the travelling traps. 30 After a trap is stopped the molecules oscillate with a period of B50 ms, as they evolve to fill it. From these considerations, the start time of the fitting process was therefore set to 50 ms.
To ascertain the range of Rydberg states suitable for decelerating and electrostatically trapping, a laser photoexcitation spectrum, shown in Fig. 7, was recorded with detection after a trapping time of t trap = 350 ms. The spectral intensities of each of the features in this spectrum reflect the combined contributions from the efficiency with which long-lived LFS hydrogenic Rydberg-Stark states were populated, the deceleration and trap loading efficiencies, and the excited-state decay rates. The reduced spectral intensities at lower values of n reflect the n 2 dependence of the maximal electric dipole moments of the Rydberg states, that reduce the deceleration and trapping efficiency for lower n states. The spectral intensities in Fig. 7 also reduce for high values of n. This is a consequence of the lower threshold for electric field ionisation during deceleration and trap loading, and the lifetimes of the states. 47 The deceleration and trapping efficiency of the NO molecules depends on effects of intramolecular interactions between the optically accessible np(0) or nf(2) states, and nearby c-mixed hydrogenic Stark states. The resonances associated with transitions to the nf(2) states in Fig. 7 arise because of their close proximity to, and mixing with, the n(2) hydrogenic Stark states. 39,47 Spectral features that arise from mixing of the np(0) states and hydrogenic Stark states are also observed. Despite the np(0) states being shorter lived that the nf(2) states, 40,41 the hydrogenic states populated in both cases are long lived. Two manifolds of c-mixed hydrogenic Stark states are accessible upon excitation because of mixing with the np(0) states. These are (i) n 0 (1) (where n 0 a n) hydrogenic Stark states that mix because of accidental degeneracies with np(0) states as a result of intramolecular charge-dipole interactions. 39,43 These coincide with resonances denoted np(0) in Fig. 7. And (ii) n(0) hydrogenic Stark states that mix with the [n + 1]p(0) states in the electric field present close to the time of photoexcitation. These states are denoted n(0) in Fig. 7.
The use of state-selective PFI to detect trapped Rydberg NO molecules allowed information on the time-evolution of the excited state populations to be determined. To achieve this PFI was performed by turning off the trapping potentials, and B1 ms later applying a slowly-rising ionising potential V Ramp (t) to the side electrodes of the decelerator. The time-dependence of V Ramp (t) is shown Fig. 8(a). The electron signal obtained at a trapping time of t trap = 0 ms, for molecules excited on the 38(2) and 43 (2) resonances, are shown in Fig. 8(b)-(i). From these data a distinction between the electron detection times, and hence the ionisation fields of the molecules, can be identified. The maxima in the signals recorded following excitation on the 38(2) and 43 (2) resonances, occur at t Ramp B 5.9 ms [V Ramp (5.9 ms) B À160 V] and B3.3 ms [V Ramp (3.3 ms) B À110 V], respectively. Since the maximum in the electron signal associated with the 43(2) resonance occurs earlier than that associated with the 38(2) resonance, it is concluded that excitation on the 43(2) resonance leads to ionisation in smaller electric fields. This is expected from the n-dependence of the ionisation field and previous detailed measurements of the electric field ionisation dynamics of long-lived hydrogenic  states in NO. 39 The widths of the distributions in Fig. 8(b)-(i) reflect the combined effects of a distribution of Stark states present in the trap at the time of PFI, the inhomogeneity of the ionisation field arising as a result of the electrode geometry in the decelerator, and the motion of the molecules in this timevarying inhomogeneous electric field after the trap is switched off. The tails in the electron signals at late times in Fig. 8(b)-(i) are attributed to secondary electron generation following PFI.
The expected ratio of the ionisation fields for LFS hydrogenic Rydberg-Stark states with different values of n was determined from the ionisation field, F ion (n), of the outermost LFS Stark state, where F 0 = 2hcR NO /ea NO . For the 43(2) and 38 (2) states this ratio is F ion (43)/F ion (38) = 0.6. This value is comparable to the ratio of the values of V Ramp (t) at the maxima of the electron distribution for the molecules excited on the 43(2) and 38 (2) resonances of V Ramp (3.3 ms)/V Ramp (5.9 ms) B 0.7. It is therefore concluded that the molecules trapped in these states predominately undergoing diabatic electric field ionisation. 39 Using this state selective PFI detection scheme, information on the evolution of the excited state populations was also obtained by comparing the electron signals recorded when t trap = 0 ms and 400 ms. This was done for the 38(2) and 43 (2) resonances with the results, scaled such that the signal maximum in each case is 1, are shown in Fig. 8(b)-(ii and iii), respectively. For molecules excited on the 43(2) resonance, no significant change in the general form of the electron TOF distribution was observed between t trap = 0 ms and 400 ms. For molecules excited on the 38(2) resonance this is not the case. Under these conditions when t trap = 400 ms the maximum of the electron signal occurred at t Ramp B 4.8 ms [V Ramp (4.8 ms) B À140 V], i.e., approximately 1.1 ms before that recorded when t trap = 0 ms. This suggests that the molecules detected at late times in this case are in higher-n states. From the ratio of the corresponding values of V trap (t), this shift in the maximum of the electron signal is estimated to correspond to a change in n of B+2.
By performing numerical particle trajectory calculations assuming diabatic traversal of all avoided crossings between Stark states, 47 the typical distribution of the maximum electric fields experienced by molecules when trapped was determined. The mean of this distribution was B150 V cm À1 . As an estimate of the electric field ionisation rates of the trapped NO molecules in the experiments, it can be determined that for hydrogenic Rydberg-Stark states with n = 38, |m| = 4, k = +33 (n = 43, |m| = 4, k = +38), tunnel ionisation would occur at a rate of B10 À30 s À1 (B10 À3 s À1 ) in this field 75,76 (k = n 1 À n 2 is the difference in the parabolic quantum numbers n 1 and n 2 obtained by solving the Schrödinger equation for a pure Coulomb potential in parabolic coordinates 77 ). Hence, it is concluded that electric field ionisation arising from the motion of molecules in the trap does not affect the distribution of states with these values of n populated over time.
In a similar vein, for hydrogenic Stark states with n = 38, |m| = 4, k = +33 (n = 43, |m| = 4, k = +38), the Dn = AE1 and AE2 blackbody transitions rates at 30 K are B9 Â 10 2 s À1 and B4 Â 10 2 s À1 (B7 Â 10 2 s À1 and B3 Â 10 2 s À1 ), respectively. At this temperature these depopulation rates scale approximately with n À2 . However, because the molecules are trapped B2 mm above a 2-mm-diameter aperture in the heat shield surrounding the decelerator (A1 in Fig. 1), the blackbody field that the molecules interact with has a room temperature component. By considering the fractional solid angle subtended by the aperture at the minimum of the trap this corresponds to approximately 0.25, not accounting for modifications in the spectrum arising from the dimensions of the electrode structure. At 300 K the Dn = AE1 and AE2 blackbody depopulation rates, for n = 38, |m| = 4, k = 33 (n = 43, |m| = 4, k = 38), are B1 Â 10 4 s À1 and B5 Â 10 3 s À1 (B7 Â 10 3 s À1 and B3 Â 10 3 s À1 ), respectively. Therefore the total Dn = AE1 and AE2 blackbody depopulation rates for these states at the position of the trap are estimated to be below 0.75g BB (30 K) + 0.25g BB (300 K) and hence B3 Â 10 3 s À1 and B1 Â 10 3 s À1 (B2 Â 10 3 s À1 and B1 Â 10 3 s À1 ), respectively. These blackbody transition rates are similar for both n = 38 and n = 43. Therefore any blackbody-induced changes in the electron distribution, measured by PFI, are not expected to differ for the molecules excited to the 38(2) and 43 (2) states. Consequently, the change in the n = 38 electron distribution with trapping time [see Fig. 8(b)-(ii)] is interpreted to reflect the initial population of the Rydberg states close to the time of photoexcitation and before deceleration and trapping. Excited states with longer lifetimes remain in the trap for longer times. Hence if an ensemble of molecules is prepared in a range of states, with different lifetimes, the measured distribution of states populated in the trap will change over time. This behaviour has been discussed previously by Seiler et al. in the interpretation of experiments with electrostatically trapped Rydberg H atoms. 74 In the case of interest here, this suggests that excitation of the NO molecules on the 38f(2) resonance resulted in the population of at least two Rydberg states with different values of n(N + ), while excitation on the 43f(2) resonance lead predominantly to the population of states with a single set of values of n(N + ).

Decay rate measurements
Measurements of the trap decay time constants, t decay , for NO molecules in n(0), n(1), or n(2) Stark states were performed for values of n between 32 and 50. The value of t decay was determined from each of these measurements by fitting a single exponential function to the experimental data between t trap = 50 ms and 350 ms, as described in Section 4.1. The weighted mean and standard deviation of up to 6 separate measurements of t decay for each set of values of n(N + ) were calculated. These values of t decay obtained for the n(0), n(1), and n(2) states, with V 0 = 149 V, can be seen in Fig. 9. For reference, the decay time constants of the 38(2) and 43 (2) states in this figure were t decay = 346.2 AE 3.4 ms and 290.0 AE 4.7 ms, respectively.
From the measured trap decay times in Fig. 9 it is seen that there is a general trend for the measured value of t decay to decrease with increasing values of n, and therefore also increasing excitation wavenumber, u 2 . This behaviour occurs for states in each N + series studied. For the n(0) and n(1) states this trend is seen for n Z 40. For the N + = 2 series, the general reduction in the value of t decay with n begins at lower values of n. But for n = 34, 35, 40, 43, 45, and 46, deviations from this general behaviour are observed. For n = 34 and 35 the measured values of t decay are shorter than those of the surrounding states that follow the general trend. While, for n = 40, 43, 45, and 46 the measured values of t decay are longer.
In the hydrogen atom, the fluorescence lifetimes of c-mixed Rydberg states scale approximately with n 4 . 77 Therefore, the observed general reduction in t decay with n for long-lived 'hydrogenic' Rydberg-Stark states in NO in Fig. 9 indicates that intramolecular interactions with short-lived states must occur. These interactions must in general increase with n. However, for particular values of n at which deviations from this general trend occur, additional contributions from interactions with more energetically localised states must arise.

Vibrational channel interactions
The general reduction in the measured values of t decay as the value of n, and hence u 2 , increase in the data in Fig. 9, occurs because of weak intramolecular interactions with short-lived low-n electronic states in series that converge to the X + 1 S + (v + = 1, N + ) state of NO + . As can be seen in the energy-level diagram in Fig. 2, n = 7 states in the v + = 1 series lie above, but close to, the v + = 0, N + = 0 series limit. The lifetimes of these states in the v + = 1 series are dominated by fast non-radiative decay processes, i.e., a combination of predissociation and vibrational autoionisation. 44 Transitions to these states have spectral widths of B5 cm À1 , and hence they are estimated to have lifetimes of B1 ps. 44 Studies have been reported previously of the Dv + a 0 vibrational channel interactions in NO between Rydberg states in series converging to the v + = 1 series limit, and the v + = 0 continuum that give rise to autoionisation. 44,78 In the case of interest here these same vibrational channel interactions are seen to persist between states with v + = 0 and v + = 1 below both series limits.
To evaluate the effects of these vibrational channel interactions in the experiments here, calculations of t decay were performed including contributions from the v + = 0 fluorescence and dissociation decay rates, as described in Section 3.3. In these calculations the interactions between states in the v + = 0 and v + = 1 series were treated semi-empirically. The amount of low-n, v + = 1 character mixed into a v + = 0 Rydberg-Stark state |ii was described by the value of a Lorentzian function, at the off-resonance wavenumber E i /hc of the Rydberg-Stark states for which the trap decay time was measured. E v + =1 /hc is the wavenumber of the peak of the Lorentzian function, G v + =1 is the v + = 1 spectral width (G v + =1 = 1/2p c t v + =1 where t v + =1 B 1 ps), and a (N + ) is a scaling factor. In the interpretation of the experimental data in Fig. 9, E v + =1 and a (N + ) were chosen to be fit parameters. Their values depended on the dominant N + character of the state |ii. The total decay rate, g tot(i) , of the state |ii was then considered to be, where g v + =1 = 1/t v + =1 is the v + = 1 decay rate, and g 0(i) is given by eqn (17). The decay time constant of the state |ii was then t tot(i) = 1/g tot(i) . This approach allowed the v + = 1 character of the n(N + ) Stark states prepared in the experiments to be estimated. Because the data in Fig. 9 were recorded for times greater than 250 ms after laser photoexcitation, the decay rates after this time were calculated using a similar procedure to that employed by Bixon and Jortner. 43 This involved considering the population of all accessible short and long lived n(N + ) Stark states, each of which decay exponentially. The resulting exponential functions were summed, to give the total Rydberg state population, P(t), where t = 0 ms corresponds to the photoexcitation time. A single exponential function was then fit to the results of these calculations between t = 300 ms and 600 ms, i.e., between the times corresponding to t trap = 50 to 350 ms. Because the deceleration and trapping process only accepted the outer LFS Rydberg-Stark states, this procedure was implemented for the outer quarter of the Stark manifold for each value of n(N + ) to obtain the calculated time constant t calc .
Calculations of t calc were performed using fixed values of the electric field, F z , and M N . The value of F z was taken to be equal to the time-averaged field experienced by the molecules in the electrostatic trap obtained from numerical particle trajectory calculations. 47 For V 0 = 149 V this field was F z = 95 V cm À1 , and was approximately constant across the range of values of n studied. The values of M N were restricted to |M N | 4 |N + À 3| such that at least one of the 2N + + 1 Stark manifolds, for each n(N + ), contained no short-lived c r 3 character, i.e., for N + = 0, 1, and 2, |M N | Z 4, 3, and 2, respectively.
The decay time constants, t calc , where first calculated without including effects of vibrational channel interactions (i.e., a (N + ) = 0 in eqn (22)). In this case only the intrinsic v + = 0 fluorescence and dissociation decay rates were considered. For the n(2) v + = 0 states, the resulting time constants are shown for M N = 2 to 5 in Fig. 10(a). From these data it is seen that when only v + = 0 fluorescence and dissociation decay rates were included the value of t calc increases with n, as would be expected in an atom. This does not however reflect the general trend observed in the experimental data in Fig. 9.
To determine the values of E v + =1 and a (N + ) for each N + series, a fit was performed of t calc to t decay . In this fitting process the reduced chi-squared w 2 ¼ P t decay À t calc À Á 2 .
Dt  Table 3. The uncertainties on these values were calculated by considering the change required in E v + =1 or a (N + ) to increase the chi-squared to w 2 = w 2 min + 1. Comparison of the best-fit values of E v + =1 in Table 3 with the energy-level structure in Fig. 2 shows that the peaks of the Lorentzian functions, which approximately represent the positions of the v + = 1 (t v + =1 B 1 ps) states mixed into the v + = 0 Stark states, lie close to or above the v + = 0, N + = 0 series limit and are in the vicinity of the calculated energies of the (n = 7)X + 1 S + (v + = 1, N + ) states. From the values of E v + =1 and a (N + ) , in Table 3 and eqn (22) it can be estimated that the fraction of short lived (t v + =1 B 1 ps) n = 7, v + = 1 character mixed into the n(N + ) v + = 0 Stark states is between 10 À10 and 10 À9 . For example, a 10 À9 contribution from a short-lived state with a lifetime of 1 ps, mixed into a long-lived state with a lifetime of 1 ms, yields a total decay rate of 2000 s À1 and a lifetime of 500 ms. Fig. 10(b) shows the fraction of (n = 7)X + 1 S + (v + = 1, N + ) character mixed into the predominantly v + = 0, N + = 2 Rydberg states below the v + = 0, N + = 0 series limit as given by eqn (22). The shaded region corresponds to the range of mixing fractions that arise from the uncertainty in the fit parameters. As can be seen in Fig. 10(b) the fraction of this n = 7, v + = 1 character mixed into the n(2) Stark states increases by approximately a factor of 5 over the wavenumber range studied. In Fig. 10(c) the effects of these vibrational channel interactions on the values of t calc for the n(2) states (black squares) can be seen, and these results are compared to the measured values of t decay (red circles). The grey shaded region in this panel of Fig. 10 corresponds to the uncertainty in t calc arising from the uncertainty in the fit parameters. The calculated decay time constant obtained for the 38 (2), M N = 4 [43(2), M N = 4] state using this approach was t calc = 343 AE 30 ms [t calc = 279 AE 35 ms]. It can be seen in Fig. 10(c) that, in contrast to the case in which only v + = 0 fluorescence and dissociation rates were considered [ Fig. 10(a)], this type of vibrational channel interaction leads to a dependence of the excited state decay times on the value of n that follows the trend observed in the experiments. The contributions from this interaction to the excited state decay rates are B1 kHz.
To further validate this interpretation of the experimental observations as arising from vibrational channel interactions, a two state quantum mechanical model was implemented. To approximate the interaction between a v + = 0 and v + = 1 state, a 2 Â 2 Hamiltonian, H vib , was constructed, where E v + =0 and E v + =1 represent the energies of the v + = 0 and v + = 1 states, respectively, and I vib represents the strength of the interaction. Within the harmonic oscillator approximation, Dv + = 1 mixing occurs when an interaction is present that causes the Hund's-case-(b) quantum defects to have a nonzero derivative with respect to the internuclear separation R, when evaluated at the equilibrium internuclear distance R e , i.e., 44,63,[79][80][81][82] dm ðbÞ 'L dR R¼Re a0: For the Rydberg states of NO the dependence of m (b) cL on R arises as a result of the interaction of the Rydberg electron with the R-dependent multipole moments of the NO + ion core. 44,59 The off-diagonal matrix elements associated with this interaction, I vib in eqn (25), can therefore be estimated by considering the interaction between the |ncN + NM N v + = 0i and |n 0 c 0 N +0 N 0 M 0 N v +0 = 1i Hund's-case-(d) basis states. This interaction can be expressed as, 63 where, for Dv + = 1, with m red and o e the reduced mass and vibrational frequency of NO + , respectively. The units of the coupling given by eqn (27) are introduced through the density of Rydberg states as in the first order approximation to MQDT. 52,80 It is assumed that the values of the Hund's-case-(b) quantum defect at the equilibrium internuclear distance, m (b) cL | R=Re , are small because c Z 4, and hence Following this approach the interaction strength I vib was calculated to be between I vib /hc C 1 Â 10 À3 cm À1 and 5 Â 10 À2 cm À1 . The effect of this interaction on the v + = 0 and v + = 1 states was then determined by considering that (E v + =1 À E v + =0 )/hc B 100 cm À1 , and calculating the eigenvalues and eigenvectors of the Hamiltonian in eqn (25). The resulting mixed states contain 10 À10 to 10 À7 of the minority component, i.e., the absolute values of the amplitudes of the minority component in the wave functions are 10 À5 to 10 À3.5 . These mixing fractions are consistent with those expected from the measured trap decay times. Using this two-state model it can also be estimated that the vibrational channel interaction considered here would give rise to wavenumber shifts of the v + = 0 and v + = 1 states of B10 À8 cm À1 to 10 À5 cm À1 (i.e., frequency shifts of B100 Hz to 100 kHz).

Rotational channel interactions
To interpret the trap decay time constants for the individual Rydberg states that deviate from the general trend in Fig. 9, i.e., that t decay decreases with increasing n, it is necessary to consider the role of rotational channel interactions between states with v + = 0. As discussed by Bixon and Jortner, 43 and in Section 4.1, the hydrogenic Stark states populated by photoexcitation in the experiments here are those that undergo intramolecular interactions with nearby optically accessible np(0) or nf(2) states. These interactions can be induced solely by external electricfields, as is the case for n(2) states that mix with nf(2) states, or by a combination of intramolecular charge-multipole interactions and external electric-fields, that occur, for example, in the case of n(1) states populated because of accidental degeneracies between np(0) and n 0 d(1) states (where n 0 a n), with the n 0 d(1) states further mixing with the n 0 (1) hydrogenic states. In Fig. 9(c) the measured values of t decay for the 40(2), 43(2), and 45(2) Rydberg-Stark states are larger than expected from the behaviour of the surrounding states. These n(2) hydrogenic Stark states are populated because of the effects of mixing with Table 3 Parameters associated with the Lorentzian function in eqn (22), representing the strength of the vibrational channel interactions between the short-lived (t v + =1 B 1 ps) low-n states with v + = 1 and the high Rydberg states with v + = 0, obtained by minimising the reduced chi-squared of t calc with respect to t decay for the different N + hydrogenic series separately. The fitting process was performed with the data recorded for V 0 = 149 V. Values of M N = 4 and F z = 95 V cm À1 where used in the calculations  (2) states. However, the charge-quadrupole interactions that couple states for which DN + = AE2 and Dc = 0 lead to mixing between nf(2) and n 0 f(0) states, which, in the presence of an electric field also mix with the n 0 (0) states. This type of interaction with intramolecular and external field contributions leads to the population of n 0 (0) states upon photoexcitation if they are energetically close to an nf(2) resonance. Hence, because of accidental degeneracies with n 0 (0) states, photoexcitation on the 40f(2), 43f (2), and 45f (2) resonances leads to the population of both the 40(2) and 44(0), the 43(2) and 48(0), and the 45(2) and 51(0) states, respectively. Since the n 0 (0) states have longer decay time constants than the n(2) states at a particular excitation wave number (see Fig. 9), trapped molecules in these states have larger values of t decay than those excited on resonances where only n(2) states are populated. This increase in the observed decay time constant also depends on the ratio of the n(2) and n 0 (0) state populations. Further deviations from the general trend in the values of t decay in Fig. 9 are also seen for the 34(2) and 35 (2) states. In these cases, charge-dipole interactions that couple states for which DN + = AE1 and Dc = AE1 lead to mixing of the nf(2) and n 0 c 0 (3) states. These result in the population of n 0 (3) states, in addition to n(2) states, upon photoexcitation. Consequently, excitation on the 34f(2) and 35f(2) resonances leads to the population of both the 34(2) and 32 (3), and the 35(2) and 33(3) states, respectively. The reduced decay time constants observed for n = 34 and 35 in Fig. 9(c), when compared to the surrounding states, suggests that the 32(3) and 33(3) states have shorter lifetimes than the 34(2) and 35 (2) states.
The measured decay time constant for the 46(2) state in Fig. 9, is larger than expected from the trend in the surrounding states. In this case the presence of both the optically accessible 46f(2) and 53p(0) states within the laser spectral width, and the increasing overlap of the Stark manifolds at these high values of n complicate the interpretation of this feature. However, based upon similar arguments to those outlined above, an estimate of the states that can be populated in both of these cases can still be made. Although a rigorous treatment, e.g., by combining particle trajectory calculations with MQDT, will be required to determine the relative populations of each Stark state and the efficiency with which each can be decelerated and trapped. In this instance it is considered that: The 46(2) Stark states are populated because of electricfield-induced mixing with the 46f(2) states.
The 50(1) Stark states are almost degenerate with the 46(2) states and can therefore be populated as a result of chargedipole interactions that (i) couple the 46f(2) states and 50d(1) or 50g(1) states, or (ii) couple the 53p(0) state with the 50s(1) or 50d(1) states. These are then mixed with the hydrogenic 50(1) states in the electric field.
Finally, 52(0) Stark states can be populated as a result of electric-field-induced mixing with the 53p(0) state. Since, at a particular excitation wavenumber the N + = 0 Stark states were observed to exhibit longer decay time constants than the N + = 1 and 2 Stark states (see Fig. 9), the presence of 52(0) states is expected to increase the value of t decay measured at the wave number associated with the excitation of the 46 (2) states when compared to the other surrounding states. This is consistent with the experimental observations.
The dependence of t decay on the trap depth was measured for the n(2) states in the range from n = 38 to 44, see Fig. 11. The data in this figure were recorded for three different trap depths with V 0 = 149 V (red circles), V 0 = 125 V (blue open diamonds), and V 0 = 100 V (black triangle). From these data, it is seen that in general t decay is insensitive to the value of V 0 , and hence the trap depth. This is consistent with the results of calculations for F z = 95 V cm À1 (V 0 = 149 V) and F z = 80 V cm À1 (V 0 = 125 V) in which the reduction in the electric field strength caused t calc to change by o5 ms over this range of values of n. However, two exceptions to this general trend are seen in Fig. 11. These occur for the 40(2) and 43 (2) states where the value of t decay increased as the value of V 0 was decreased.
As discussed above, photoexcitation on the 40f(2) and 43f(2) resonances results in the population of the 40(2) and 44(0), and the 43(2) and 48(0) Stark states, respectively. The increase in t decay as V 0 was decreased at these resonances indicates that in these cases the fractions of trapped molecules in the longer lived 44(0), or 48(0), states also increase as V 0 decreases. For states with n 4 40 the trapping efficiency is limited by electric field ionisation during deceleration and trapping, with the ionising field scaling with n À4 . For lower values of V 0 , the electric field experienced by the molecules during deceleration and trapping is reduced. This leads to an increase in the trapping efficiency for molecules in higher n states, because of reduced loses from electric field ionisation and the resulting increase in the effective volume of the trap. This increase in trap loading efficiency is greater for states with higher values of n. Therefore, for the measurements performed following excitation on the 40f(2) resonance, the trapping efficiency for molecules in the 40(2) states only changes a little between V 0 = 149 V and 125 V, whereas the efficiency of trapping molecules in the near degenerate 44(0) states increases as the values of V 0 is reduced. This results in an increase in the fraction of molecules in the longer-lived 44(0) states present in the trap when V 0 = 125 V, compared with V 0 = 149 V, and therefore t decay is observed to increase. The situation is similar for the 43f(2) resonance, where the trapped molecules are in the 43(2) and 48(0) states, and the fraction of molecules in the longer-lived 48(0) states increases as V 0 decreases resulting in an increase in t decay .

Conclusions
In conclusion, NO molecules have been prepared in high Rydberg states with lifetimes in excess of 100 ms, decelerated to rest in the laboratory-fixed frame of reference in the travelling traps of a transmission-line Rydberg-Stark decelerator and electrostatically trapped. Trapping was performed in a cryogenic environment to minimise the effects of blackbody radiation induced transitions. Measurements, of the decay of the molecules from the trap, were conducted in situ by pulsed electric field ionisation. The long lifetimes of the Rydberg states observed in the experiments are reminiscent of those of Rydberg states in atoms. However, over the timescales of B1 ms on which the experiments were performed, molecular effects were seen to play a significant role in the decay of the trapped molecules, and result in a breakdown of n-scaling rules typical of hydrogenic systems. These effects arise because of rotational and vibrational channel interactions which originate from the coupling of the Rydberg electron to the electric dipole, quadrupole, and higher order multipole moments of the NO + ion core. These types of channel interactions have been discussed previously in the literature, in the context of the lifetimes of Rydberg states in ZEKE spectroscopy experiments. However, ZEKE experiments are generally performed on the 1 ms timescale and many of the effects of these interactions are weak. Consequently, the quantum-statedependence of each of the interaction mechanisms identified here, which typically contribute B1 kHz to the excited state decay rates, were not observable. In the experiments reported here, with measurement times of up to 1 ms, these effects are clearly seen. Vibrational channel interactions between Rydberg states in the nX + 1 S + (v + = 0, N + ) series, and short-lived (n = 7)X + 1 S + (v + = 1, N + ) states, lead to a general reduction in the lifetimes of the nX + 1 S + (v + = 0, N + ) Rydberg states with N + = 0, 1, and 2 and values of n 4 35 as each series limit is approached. Rotational channel interactions within the v + = 0 series have been identified to specifically affect the decay time constants of some individual Rydberg states.
The effects of molecular interactions on the lifetimes of Rydberg states studied here must be carefully considered when employing general n-scaling rules in the analysis of excited state decay, e.g., following electron-ion recombination. The experimental methods developed and implemented in this work offer opportunities for similar studies with other molecules, e.g., N 2 and O 2 , provided schemes for populating long-lived Rydberg states are identified. The results reported here will provide important input for future low-energy collision experiments using Rydberg NO molecules, and the implementation of schemes to exploit the technique of Rydberg-Stark deceleration to prepare cold samples of ground state molecules.

Conflicts of interest
There are no conflicts to declare.