Electronic and vibrational relaxation dynamics of NH3 Rydberg states probed by vacuum-ultraviolet time-resolved photoelectron imaging

Time-resolved dynamics of high-lying Rydberg states of ammonia (NH3) prepared by using a vacuum ultraviolet (VUV) pump (∼9.3 eV) and an ultraviolet (UV) probe (∼4.7 eV) pulse are reported using photoelectron imaging detection. After photoexcitation, two main features appear in the photoelectron spectrum with vertical binding energies of ∼1.8 eV and ∼3.2 eV and with distinctly different anisotropy parameters β of ∼1.3 and ∼0.7, respectively. This information allows the unambiguous assignment of the respective Rydberg states and disentangles the induced electronic and vibrational dynamics. The combination of velocity-map imaging with femtosecond VUV and UV pulses is shown to offer an attractive approach for studying the dynamics of high-lying Rydberg states of small molecules. © 2019 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/). https://doi.org/10.1063/1.5116707., s


I. INTRODUCTION
Many spectroscopic and theoretical studies have been devoted to ammonia (NH 3 ). [1][2][3][4][5] Studies which have focused on electronically excited NH 3 have been mainly focused on the low-lying Rydberg states (∼6 eV) investigating their complex photodissociation dynamics. 4,6,7 The few studies that have involved high-lying Rydberg states have been confined to static measurements 1,8,9 because of a lack of vacuum ultraviolet (VUV) femtosecond laser sources.
In many biologically relevant molecules, such as amino acids, the sp 3 hybridized nitrogen is the central functional unit in defining their photochemistry, e.g., through valence-Rydberg mixing that causes their surprising stability to UV light. 10,11 Further understanding of the mechanisms by which these nitrogen centers are able to efficiently redistribute absorbed UV light therefore has a natural imperative that invites further study. While there are certainly differences between the photoinduced molecular dynamics of C-N bonds and N-H bonds, 12 the study of small molecules can shed light on the fundamental dynamics of moieties that could be clouded by increased complexity, in addition to having the benefit of being accessible to state-of-the-art quantum-chemical calculations without facing extreme computational costs. NH 3 makes an ideal candidate system as it is a small molecule with an sp 3 hybridized nitrogen center. Moreover, hydrogen substitution leads to the creation of nitrogen containing biomolecules, and its spectroscopy can be interpreted in terms of simple symmetry arguments. Due to this, many studies have been devoted to the spectroscopy of NH 3 . The early works studying the photodissociation dynamics were pioneered by Ashfold and co-workers, 4,13-15 who paved the way for most of the subsequent interest in the Ã state of NH 3 . Since then, many research groups contributed to this effort. [15][16][17][18][19] However, these experiments were limited to static  [20][21][22][23][24][25][26] The natural evolution of this experimental effort in understanding the photophysics of NH 3 is then to use these femtosecond UV/VUV sources and investigate the time-resolved dynamics of the high-lying Rydberg states in NH 3 .
In the present study, the dynamics of the high-lying Rydberg states of NH 3 are followed through time-and angle-resolved photoelectron spectroscopy. The states involved are fully characterized utilizing a time-resolved velocity-map-imaging (VMI) spectrometer with a VUV pump (133 nm) and a UV probe (266 nm). One of the strengths of VMI spectroscopy is that the technique allows one to capture the photoelectron spectrum while simultaneously obtaining the photoelectron angular distribution (PAD). This PAD is related to the orbital-angular-momentum character of the ejected photoelectron, which carries information about the corresponding character of the photoionized orbital. The time-and energy-dependencies of the PAD can assist the interpretation of the observed dynamics or even provide additional insights.
The presented study discusses the time-resolved Rydberg-state dynamics of NH 3 in the high-lying Rydberg states (∼9.3 eV) and their ultrafast electronic and vibrational relaxation, probed via 266 nm pulses. This wavelength was chosen because it allows the relaxation dynamics to be followed over the set of intermediate states with higher vertical binding energies than the initially prepared states. The paper is organized as follows: Section II describes the experimental setup and summarizes the typical experimental conditions, followed by a brief description of the data analysis. Section III presents the results of time-resolved photoelectron spectra as well as the time-and angle-resolved photoelectron distributions. Section IV discusses the experimental results in terms of the assignment of individual Rydberg states and their anisotropy parameters β 2 . The dynamics of these states are discussed in the framework of kinetic modeling providing the time constants for the involved processes. The conclusions are presented in Sec. V.

A. Experimental
The experimental setup has been described in detail elsewhere. 23 The experiment is driven by 1.5 mJ 800 nm pulses from a Ti:Sa regenerative amplifier (Coherent Legend Elite Duo) with a pulse duration of 35 fs operating with a repetition rate of 1 kHz. The fundamental beam is passed through a 70:30 beam splitter to create two arms of the interferometer, with the transmitted component being used to generate the 266 nm probe through third-harmonic generation in nonlinear crystals and the reflected component used to generate the 133 nm pump, as follows. The 800 nm beam is frequency doubled using a 300 μm thick β-barium borate (BBO) crystal to generate a 400 nm beam. This beam is then spectrally filtered and focused on a semi-infinite gas cell for low-order harmonic generation. The third harmonic beam centered at 133 nm (∼9.3 eV) is generated in xenon (∼10 mbar). Afterward, the linearly (and parallel) polarized beams in both arms are spectrally filtered, time delayed with respect to each other using a motorized delay stage, and noncollinearly (crossing angle <1 ○ ) focused (f = 500 mm) onto the interaction region of a VMI spectrometer.
An Amsterdam pulsed valve (MassSpecpecD BV, 1 kHz repetition rate) using a piezocantilever opening mechanism 27-29 with 150 μm orifice and 40 ○ degree conical nozzle shape for enhanced center-line intensity is used to generate the unseeded molecular beam of NH 3 (PanGas, purity 99.98%). The typical backing pressure is between 1 bar and 1.5 bar, and the pulse width of the valve is 40 μs. The molecular beam created by the nozzle is sent through a 1 mm skimmer and propagated 10 cm downstream to the interaction region where it is ionized by the combined action of the laser beams. The laser cross correlation (CC) was determined via (1 + 1 ′ )-photon nonresonant ionization of xenon and is modeled as a Gaussian function with a full-width at half-maximum (FWHM) of 92(5) fs.
For each time step, four images are recorded: two images with only the pump or probe beam, one image with both beams, and one background image with no beam. Altogether, the whole time delay axis is measured 30 times with a single-image acquisition time of 20 s. The multiple images at each time step are background corrected and averaged. The pump-probe signal is obtained by subtracting single-color images from images recorded with both beams. The pump-probe images are then inverted using the MEVELER inversion method. 30 31 All observed signals are due to the excitation by the pump pulse (∼9.3 eV), resulting in two main features decaying within 600 fs after time zero. It is evident from the TRPES that the maxima of each of the two features are reached at different time delays and shifted with respect to time zero.

III. RESULTS
The TRPES clearly shows that the entire photoelectron signal is seen below 5 eV, which can be rationalized by considering that ultrafast photoelectron spectroscopy is measuring vertical transitions, and the first vertical ionization potential of NH 3 is 10.85 eV 32 (the adiabatic ionization potential is 10.19 eV), 33 and the total photon energy is ∼14 eV which gives a maximum of PKE at ∼3.15 eV; as no signals are observed above this PKE, multiphoton processes can be excluded. To obtain further insight into the observed TRPES, a time slice at time zero is taken and modeled by a sum of Gaussians [shown in Fig. 1(b)]. From the best fit (fitting parameters shown in Table I), two Gaussians centered at 1.76 eV and 3.15 eV are necessary to reproduce the spectrum. Additional Gaussians did not result in any improvement of the fit, confirming the two-feature interpretation of the TRPES. These two features are denoted with "U" and "L," corresponding to the features originating from the higher and lower PKE signals, respectively. This is shown in Fig. 1(b) by the letters "U" and "L" in the yellow circles. Figure 2(a) shows the result of the global fit of the measured TRPES. The fit relies on a time-slice analysis, where each time slice (representing static photoelectron spectra) is modeled by a sum of Gaussians. The advantage of this method is that no specific time dependence needs to be assumed. This differs from traditional approaches which model the time evolution as a sum of exponential decays convoluted with the instrument-response function. Moreover, the employed model allows for interpretation of a complex wavepacket evolution because the Gaussians at each time step can change both their positions as well as their widths. The relative difference between the measured TRPES and the fitted TRPES is always within a relative deviation of ±10%. The fit results in a separation of the measured signal into two components, which are shown in Figs. 2(c) and 2(d) centered around PKE ∼ 1.8 eV and PKE ∼ 3.2 eV, respectively. Both components show no time-dependent change in the peak position, and the signal of both components decays to zero within less than 600 fs, albeit on two different timescales.
The time-dependent photoelectron angular distribution (PAD) is shown in Fig. 3(a) as a 2D map of the β 2 (Δt, PKE) parameter in a false-color representation. In the white area, the photoelectron signal is too weak, ≤6% of its maximum, to get a reliable β 2 parameter; therefore, this area is excluded from further analysis. Looking at Fig. 3(a), two distinct regions of the β 2 parameter can be identified, one with a high anisotropy (β 2 → 2) at PKE between 2 eV and 5 eV and a more isotropic (β 2 → 0) region at PKE between 0 eV and 2 eV. The β 2 parameters are indicative of the angular momentum components of the ionized orbital. 34 Figure 3(b) shows the time evolution of the β (i) 2 (Δt) parameter for each component i ∈ {l, u} from the TRPES, averaged over the peak width of the given component. For this, the energy-dependent distribution functions Q 2 (PKE) for each time delay Δt, which are related to the Legendre polynomials P 2 by where S(PKE, θ) is the measured photoelectron image at a given time delay as a function of energy and angle, are calculated. As the photoelectron image is, in principle, a sum of the Q 2 functions due to different ionization channels, it is possible to separate the β 2 parameters of each channel by multiplying the β 2 map with the TRPES to obtain the distribution function Q 2 . Then, the slice analysis is performed upon the Q 2 map to obtain the Q  parameter has an initial value of ∼1.5 at time zero, indicating a strong anisotropy. The parameter then decays rapidly within the first 100 fs to a value of ∼1.25. The β (l) 2 parameter is impossible to evaluate before 100 fs, due to the photoelectron signal of this band being on the level of the noise. From this point onward, it is time-independent with a value of ∼0.7.

IV. DISCUSSION
NH 3 in its electronic ground state is a pyramidal molecule with a low barrier along the umbrella-inversion mode resulting in a tunneling splitting of ∼0.1 meV, which can be observed in high-energy-resolution experiments. If the tunneling splitting is resolved, the umbrella inversion is considered as "feasible" and the proper molecular-symmetry group in this case is D 3h (M). 35,36 In the present experiments, the energy resolution is not sufficient to resolve the tunneling splitting; hence, the ground state should be described in the C 3v group. However, the Rydberg states have a planar equilibrium geometry because they converge to the electronic ground state of NH 3 + which is planar. 1 Therefore, the Rydberg states are described in the D 3h group, irrespective of the energy resolution. 1,37 NH 3 is a small molecule with a central nitrogen atom in a "ligand-field" of three hydrogen atoms and is isoelectronic with Ne. Hence, it can be well described using the united-atom limit, allowing symmetry assignments of individual orbitals of NH 3 to be derived by a D 3h field splitting of the atomic orbitals of Ne (summarized in  Table II). Consequently, s orbitals are of a ′ 1 symmetry, pz orbitals are of a ′′ 2 symmetry (as they are antisymmetric with respect to the σ h plane), and px,y orbitals are e ′ . The d z 2 orbitals are totally symmetric; therefore, they have a ′ 1 symmetry. dxz,yz orbitals are degenerate and antisymmetric with respect to the σ h plane; therefore, they have e ′′ symmetry, and the d xy,x 2 −y 2 orbitals are degenerate and have e ′ symmetry because they are symmetric with respect to the σ h plane. The highest occupied molecular orbital of NH 3 is the nitrogen lone pair, having 2pa ′′ 2 symmetry. The dynamics presented here are those following single-photon excitation of one electron of the 2pa ′′ 2 orbital to a high-lying Rydberg state, converging to the cationic ground state, which is characterized by the ion core (2pa ′′ 2 ) −1 .    also Table I). In the energy range around ∼9.3 eV, there are at least three Rydberg statesẼ ′′ 1 A ′′ 2 ,Ẽ ′′′ 1 A ′′ 2 , andẼ 1 E ′′ which are accessible via one-photon excitation. 8,38 The former one is of s-character, whereas the remaining ones have d-character. Ionization of all three states would lead to a photoelectron with very similar PKE. The upper peak in the TRPES is assigned as theẼ ′′ state based on its PKE and its PAD [ Fig. 3(a)], which is consistent with a Rydberg state of s character, making contributions from other states very unlikely, because they would result in a more isotropic PAD.
Previous work by Langford et al. 1 observed four Rydberg states in the vicinity of 9.3 eV: using two-photon excitation. In contrast to single-photon excitation under which the first three states are dipole forbidden (in the atomic limit) because they are of p character, in two-photon excitation, the dipole selection rules are less restrictive and hence these states are observed. Since theẼ ′′ state is dipole forbidden in two-photon excitation, it appears as a very weak transition in the work by Langford et al. 1 The assignment of the lower band in TRPES is slightly more involved. The band is centered at 1.76 eV, which is consistent with ionization of theC ′ state. 8,38 This assignment is additionally supported by the observed PADs, which display an anisotropy consistent with a p-type Rydberg state. However, single-photon excitation of a p state is dipole forbidden in the atomic limit. This suggests that theC ′ state is populated through internal conversion, mediated by strong vibronic coupling between theẼ ′′ state and theC ′ state. This interpretation is supported by the delayed rise of the signal of theC ′ state with respect to theẼ ′′ state and also by the time-dependent change in the asymmetry parameters of the two states.
As is visible in Fig. 2, the photoelectron signals corresponding to theẼ ′′ andC ′ states evolve differently in time. This is even better visible in Fig. 4, which shows the energy-integrated signal of each peak, as obtained from the fitting procedure discussed above [see Figs. 2(c) and 2(d)]. It is clear from Fig. 4 that the signal maximum of each state is offset from time zero, located at the maximum of the cross correlation curve. A reasonable fit of theẼ ′′ state is obtained by assuming a mono-exponential decay of the population, where g(t) is the measured cross correlation function (shown in each panel of Fig. 4 as the red dotted-dashed line) and τ 1 is the decay time of the state. However, this does not provide a good fit of thẽ C ′ -state signal. Therefore, a first-order sequential kinetic model is assumed, where τR gives the rise time of the state population and τD is its decay. Figure 2(b) shows a schematic energy level diagram illustrating the proposed dynamics. Excitation of the ground state by 133 nm light prepares theẼ ′′ state, which then decays to theC ′ state due to strong vibronic coupling that is likely mediated via the umbrella mode.
The time constants obtained from the fits are shown in Fig. 4 and summarized in Table III. TheẼ ′′ signal reaches its maximum within the cross correlation time, meaning that the state is directly excited by the pump pulse. This state then decays on an ultrafast time scale with τ 1 of ∼150 fs. This time constant is shorter than the expected lifetime based on the linewidths from the high-resolution spectra. 1 The time constant must therefore represent the evolution of the photoexcited wavepacket from the initial Franck-Condon region to a configuration space where either the wavepacket cannot be energetically ionized by one photon of 266 nm or the ionization probability drops to zero for other reasons.
Contrary to theẼ ′′ state, the photoelectron signal of theC ′ state has a rise time of ∼86 fs. This observed rise is likely a consequence of ultrafast population transfer from a higher state, meaning that the initially prepared wavepacket on theẼ ′′ state populates theC ′ state via internal conversion, most likely mediated by a conical intersection. After this, the photoelectron signal decays on an ultrafast time scale of 200 fs, for similar reasons as previously discussed.
The population transfer via internal conversion is supported by the time dependence of the anisotropy parameters, as is shown in Fig. 3(b). This change is indicative of strong vibronic coupling between these two states, which can be rationalized with symmetry arguments by considering the D 3h point group. The upper state (Ẽ ′′ ) is an s-type Rydberg state, so the ejected photoelectron would be expected to display a strong anisotropy, which is seen at time zero. This anisotropy becomes weaker in the first 100 fs, before remaining   . . 150 ± 9 Monoexponential fit LowerC ′ 1 A ′ 1 86 ± 28 199 ± 40 Sequential kinetic model more or less stable at ∼1.2. The lower state (C ′ ) is a p-type Rydberg state, which would lead to s and d partial waves upon photoionization. Depending on the mixing of the two, the PAD could be expected to be nearly isotropic. The observed PAD is not isotropic and has a time-independent value of ∼0.7. This suggests that the (C ′ ) state is mixed with theẼ ′′ state, due to the population transfer. This behavior is consistent with previous observations of ultrafast internal conversion, mediated by a conical intersection. 39 The observed coupling between theẼ ′′ andC ′ states is likely due to the umbrella mode. The NH 3 umbrella mode has A ′′ 2 symmetry, meaning that this motion is allowed to mediate the relaxation processẼ ′′ →C ′ . 40 The observed change in the PAD for the (Ẽ ′′ ) state is due to this vibronic coupling between the two states, resulting in neither state being purely of s or p character, but instead having a mixed s/p character. 39,41 The assumption that the umbrella mode would be the coordinate of interest is justified by the equilibrium geometries of the ground state and Rydberg states of NH 3 , meaning that it is the umbrella mode that takes the system out of the Franck-Condon region, which would also explain the ultrafast population transfer to theC ′ state.

V. CONCLUSIONS
Our time-resolved photoelectron imaging experiment of the high-lying Rydberg states of NH 3 , accessed by a table-top VUV light source, shows an ultrafast population transfer from theẼ ′′ state to theC ′ state due to internal conversion mediated by the umbrella motion. TheẼ ′′ state is directly populated with the pump pulse and decays on the order of ∼150 fs. Afterward, the ultrafast population transfer to theC ′ state occurs within <100 fs, and the state as a whole is shown to undergo ultrafast decay within ∼200 fs. The angular information obtained by using VMI allows these states to be unambiguously assigned on the basis of the anisotropy parameters and the kinetic energy of the ejected photoelectrons, allowing one to track the evolution of the prepared wavepacket. The interpretation presented is only made possible by knowing the time evolution of the PAD as the key information is the mixed character of each state evident from the time-dependent anisotropy of the ejected photoelectrons. This demonstrates the strength of VMI in time-resolved studies, especially in cases where the complex dynamics are presented and both energy-and angular-distributions are needed for unambiguous assignments.