Relaxation dynamics in excited helium nanodroplets probed with high resolution, time-resolved photoelectron spectroscopy

A. C. LaForge *a, J. D. Asmussen b, B. Bastian b, M. Bonanomi cl, C. Callegari d, S. De e, M. Di Fraia d, L. Gorman a, S. Hartweg f, S. R. Krishnan e, M. F. Kling ghij, D. Mishra a, S. Mandal k, A. Ngai f, N. Pal d, O. Plekan d, K. C. Prince d, P. Rosenberger gh, E. Aguirre Serrata a, F. Stienkemeier f, N. Berrah a and M. Mudrich be
aDepartment of Physics, University of Connecticut, Storrs, Connecticut 06269, USA. E-mail: aaron.laforge@uconn.edu
bDepartment of Physics and Astronomy, Aarhus University, 8000 Aarhus C, Denmark
cIstituto di Fotonica e Nanotecnologie, Consiglio Nazionale delle Ricerche, 20133 Milano, Italy
dElettra-Sincrotrone Trieste S.C.p.A., 34149 Basovizza Trieste, Italy
eDepartment of Physics, Indian Institute of Technology-Madras, Chennai 600036, India
fInstitute of Physics, University of Freiburg, 79104 Freiburg im Breisgau, Germany
gDepartment of Physics, Ludwig-Maximilians-Universität Munich, 85748 Garching, Germany
hMax Planck Institute of Quantum Optics, 85748 Garching, Germany
iSLAC National Accelerator Laboratory, Menlo Park 94025, CA, USA
jDepartment of Applied Physics, Stanford University, Stanford 94305, CA, USA
kIndian Institute of Science Education and Research, Pune 411008, India
lDipartimento di Fisica Politecnico, 20133 Milano, Italy

Received 20th July 2022 , Accepted 14th November 2022

First published on 15th November 2022


Abstract

Superfluid helium nanodroplets are often considered as transparent and chemically inert nanometer-sized cryo-matrices for high-resolution or time-resolved spectroscopy of embedded molecules and clusters. On the other hand, when the helium nanodroplets are resonantly excited with XUV radiation, a multitude of ultrafast processes are initiated, such as relaxation into metastable states, formation of nanoscopic bubbles or excimers, and autoionization channels generating low-energy free electrons. Here, we discuss the full spectrum of ultrafast relaxation processes observed when helium nanodroplets are electronically excited. In particular, we perform an in-depth study of the relaxation dynamics occurring in the lowest 1s2s and 1s2p droplet bands using high resolution, time-resolved photoelectron spectroscopy. The simplified excitation scheme and improved resolution allow us to identify the relaxation into metastable triplet and excimer states even when exciting below the droplets' autoionization threshold, unobserved in previous studies.


1 Introduction

Helium (He) nanodroplets are fascinating quantum fluid clusters with distinct properties compared to other types of atomic and molecular clusters. The constituent He atoms are loosely bound to one another by extremely weak attractive London dispersion forces and their light mass implies a large zero-point energy, i.e., the emergence of collective quantum behavior. Notably, He nanodroplets evaporatively cool to the ultralow temperature of 0.37 K, where they exhibit microscopic superfluidity.1–4 Although considered inert, upon excitation or ionization, He nanodroplets can become a highly reactive environment where numerous interatomic processes can occur.5

In general, electronic excitations in nanodroplets are localized on single He atoms (He*) which tend to form void cavities (‘bubbles’) due to the repulsive interaction between the excited electron and the surrounding He.6–8 This leads to a broadening of the excited states in absorption spectra,9 although they largely retain their atomic character. We will refer to these broadened features as ‘bands’ but we label them with their corresponding atomic electronic configuration. He2* excited dimers may form, but the fraction of those directly formed by ultrafast association of an excited He* atom and a ground state He atom is very small.10,11 The main formation mechanism of He2* in large He droplets as well as in bulk superfluid He, tunneling into excited vibrational states of He2*, is much slower (∼15 μs) than any dynamics probed in the present experiment.12,13

Quasi-free electrons in an excited nanodroplet occupy states in a ‘conduction band’ ≳1 eV above the vacuum level.14,15 The corresponding electronic wavefunctions are localized in the interstitial spaces between He atoms.16 Autoionization of these highly-excited states then leads to the emission of small, singly-charged He clusters Hen+, n = 2, 3,… and low-energy electrons.17,18

The precise electronic structure of excited He nanodroplets still remains to be fully resolved.7,10 In particular, the dynamics of electron localization, relaxation, and atomic rearrangements induced by the electronic excitations are subjects of ongoing research.11 The ultrafast dynamics of excited He nanodroplets has been studied in several femtosecond pump–probe experiments employing both laser-based high-harmonic generation (HHG) sources10 and extreme-ultraviolet (XUV) free-electron lasers (FELs).19,20 Essentially, ultrafast localization of the excitation on single atomic sites was confirmed, followed by the electronic relaxation into metastable states, and the formation of nanoscale bubbles around the excited atoms on the timescale of ∼1 ps. In the first HHG-based experiments where near-infrared probe pulses were used, ejection of He* excited atoms was observed.8 Subsequent experiments by the same group showed that UV pulses are better suited for probing the full He droplet dynamics and that the ejection of a He* atom is actually a minor relaxation channel.10,21 The lowest excited atomic state of He, 1s2s3S, has only been observed at photon energies ≳ 23 eV, where He droplets autoionize and optically inaccessible states can be populated through electron–ion recombination.20,22 Static measurements have shown that even highly excited He droplet states efficiently relax into the lowest excited 1s2s1,3S atomic states.23 Eventually, He2* excimers form at the He droplet surface, detaching from the droplets over the course of vibrational relaxation.13,22

In this article, we present an overview of the XUV-pump, UV-probe dynamics of resonantly excited He nanodroplets over a wide range of excitation energies and electron energies. Subsequently, we focus on the electronic relaxation processes occurring in the lowest 1s2p and 1s2s droplet bands.9 Our previous measurements19,20 have been limited by the spectrometer resolution and by experimental noise from XUV/UV stray light. Here, with improved resolution and a higher signal-to-noise ratio, we are able to observe new features due to relaxation in the time-resolved photoelectron spectra. In particular, when exciting the system to the lowest optically allowed droplet states, 1s2s1S and 1s2p1P, we observe efficient relaxation into triplet atomic states, indicating droplet-induced spin-relaxation, as well as relaxation into the first excimer state of He2*. Compared to previous work at higher excitation energies,10,20 the relaxation process discussed here is simpler since it only involves intra-band relaxation as opposed to a combination of inter- and intra-band relaxations. Thus, with a simpler excitation scheme and improved resolution, we are able to better distinguish weak features in the photoelectron spectrum and to probe their ultrafast dynamics on femtosecond timescales.

2 Experimental methods

The pump–probe experiment was performed at the Low Density Matter (LDM) endstation24 of the seeded FEL FERMI in Trieste, Italy.25,26 The FEL pump pulse, operated at a repetition rate of 50 Hz, was primarily tuned to a photon energy 1 near the 1s2p droplet resonance (21.6 eV)9via the seed laser and undulator gaps, yielding a pulse length tp ≈ 100 fs (full width at half maximum, FWHM).27 The FEL pulse intensity was adjusted to limit the effects of multiphoton ionization28–30 with typical values of IXUV ≈ 1 × 1010 W cm−2, as derived by the pulse energy measured upstream by gas ionization, taking into account the nominal reflectivity of the optical elements in the beam transport system.31 The UV probe pulse was obtained from a frequency-tripled Ti:Sapphire laser (2 = 4.65 eV) with a pulse intensity of IUV ≈ 1 × 1011 W cm−2.32 A tin filter of 200 nm thickness was used to suppress contributions from higher order harmonic radiation whenever the pump photon energy 1 was tuned to 21.0 eV and 21.6 eV. The cross correlation between the FEL and the probe laser was 200 fs FWHM, measured by resonant (1 + 1′) two-photon ionization of atomic He, where 1′ indicates one-photon absorption by the probe pulse.32

A supersonic jet of He nanodroplets was produced by expansion of high pressure He gas (50 bar) through a pulsed, cryogenically cooled (14 K) Even-Lavie-type nozzle (50 μm). From the expansion conditions (backing pressure and nozzle temperature), the mean droplet size was estimated to NHe ≈ 105 He atoms per droplet.33,34 After expansion, the nanodroplet beam passed through a 1 mm skimmer and traversed approximately 90 cm to the interaction region. The nanodroplet beam was perpendicularly crossed by the FEL and UV beams at the center of a high resolution (EE ≈ 50) magnetic bottle electron spectrometer.35

3 Results and discussion

When a He nanodroplet is resonantly excited by an intense XUV pulse, a wide variety of dynamic processes is initiated, including internal relaxation, interatomic energy transfer, and multiphoton ionization. To probe the dynamics, we use a time-delayed UV laser pulse which directly ionizes the excited nanodroplets and we measure the energy of emitted electrons. In this way, we track the dynamics of the various processes as a function of the pump–probe delay. Fig. 1(a) shows representative raw electron kinetic energy distributions measured for He nanodroplets that are resonantly excited into the lowest optically allowed droplet state 1s2s1S at 1 = 21.0 eV for selected pump–probe delays. Fig. 1(b) shows a map of the electron kinetic energy distributions (x-axis) as a function of the pump–probe delay time (y-axis). To better understand the spectral features, Fig. 1(c) shows an energy-level diagram of He atoms (horizontal black lines) and He droplets (gray shaded areas). The vertical arrows depict photo-excitation (red) and photo-ionization steps (black). The horizontal pink arrow indicates autoionization of two excited He atoms in the same droplet, which is discussed in detail below.
image file: d2cp03335f-f1.tif
Fig. 1 Overview of the time-resolved photoelectron spectra of He nanodroplets measured with XUV-pump and UV-probe pulses. The pump photon energy is 1 = 21.0 eV and the probe photon energy is 2 = 4.65 eV. Panel (a) shows photoelectron spectra at selected pump–probe delays on a logarithmic scale; panel (b) shows a map composed of 28 spectra at different pump–probe delays. The shaded areas in (a) indicate the regions in the electron spectrum where 1–2 UV probe photons are absorbed by the XUV-excited He droplets. The vertical dashed lines mark the energies of electrons created by ICD in He nanodroplets assuming atomic excitation energies, and 3′ + 1′-photoionization of N2. Panel (c) shows an energy level diagram illustrating transitions induced by the pump pulse (1, red vertical arrows) and the probe pulse (2, blue arrows). The ICD process, where one excited He* atom decays to the ground state and another He* is ionized, is depicted by the horizontal pink arrow. The ATI process, where an electron is promoted to high kinetic energies by absorption of one or more photons in the ionization continuum, is represented by the upper blue vertical arrow.

Overall, numerous dynamic features can be observed in the spectrum in different ranges of the electron energy. At low electron energies ≲ 2 eV, we observe prominent features shifting to even lower energies as the delay increases. These are generated by 1 + 1′ photoionization through singly excited states of the droplets. The shifting of the peaks reflects the relaxation of the excited states into low-lying, metastable states, as previously studied for the 1s2p band19 and the 1s3p band.20 Additionally, we observe a similar feature at slightly higher kinetic energies around 5 eV, which is due to the excited state absorbing two UV photons from the probe pulse resulting in 1 + 2′ ‘above-threshold ionization’ (ATI) in the excited system. ATI of multiply excited He nanodroplets by near-infrared or visible probe pulses has recently been found to be strongly enhanced compared to isolated excited He atoms, pointing at a collective coupling effect.36 At electron energies around 3 eV, there is a time-independent peak, which we attribute to the photoionization of background nitrogen molecules by the probe pulse through a 3′ + 1′ REMPI scheme. Note that the peak cancels out when the background gas spectrum is subtracted; therefore, it is not related to the He nanodroplet beam.

The feature at electron energies around 16 eV is due to autoionization of pairs of excited He atoms in one droplet according to a process known as ICD.37,38 For this particular case, ICD occurs via the process He* + He* → He + He+ + eICD which manifests itself in the emission of electrons eICD with a characteristic energy. For the given experimental parameters, we estimate the number of resonantly excited He* atoms per He droplet to be N* = NHeσXUVtpIXUV/1 ≈ 700, where σXUV ≈ 25 Mbarn.28,30 This leads to an estimated number of photoelectrons emitted out of the lowest-lying He excited states NPIe = N*σUVtpIUV/2 ≈ 100, given that the photoionization cross sections, σUV, of these states range 4 to 7 Mbarn.39 The number of electrons created by ICD can be estimated based on the previously measured ICD efficiency, pICD ≈ 30%, for a droplet excitation rate of N*/NHe ≈ 0.7%,40 yielding NICDe = N* × pICD/2 ≈ 100. The factor 1/2 accounts for the fact that it takes a pair of He* atoms to create an eICD. Thus, we expect to observe roughly an equal number of photoelectrons and ICD electrons, in agreement with the experimental finding (Fig. 1).

At positive pump–probe delays this ICD process is partly quenched because the probe pulse depletes the population of He* excited states prior to their decay. From the rise of the ICD signal at delays ≳0.2 ps we previously inferred an ICD lifetime on the order of a few hundred fs.40 We observed only a weak dependence of the ICD lifetime on the degree of He droplet excitation which is controlled by the photon flux IXUV/1 and the droplet size NHe. This finding was rationalized by a model in which the bubbles forming around He* excitations can merge, thereby accelerating two nearby He* atoms toward each other.40 Thus, the detected ICD lifetime is mainly determined by the quantum fluid dynamics of the merging bubbles, rather than by the distance-dependent ICD rate for a bare pair of He* atoms. When the He nanodroplet is strongly excited with more intense FEL pulses, it can spontaneously evolve into a nanoplasma state by ‘collective autoionization’ (not shown in these spectra).28–30 In that case the electron spectrum is dominated by a distribution of low-energy electrons due to thermal emission.

For the remainder of this article, we will specifically focus on the relaxation dynamics of resonantly excited droplets, which are observed via 1 + 1′ photoionization at low electron energies. These measurements were performed with a high-resolution magnetic bottle spectrometer enabling us to distinguish new features, which previously remained unobserved.19Fig. 2 shows the electron spectra (y-axis) as a function of pump–probe delay (x-axis) for two photon energies, 1 = 21.0 eV and 21.6 eV in (a) and (c), respectively. The zero of the pump–probe delay, t0, is determined from the falling edge of the integrated ICD signal as it is quenched by the probe pulse. Fig. 2(e) shows the absorption spectrum of He droplets recorded by Joppien et al.9 by means of fluorescence detection along with He* energy levels and arrows indicating the relaxation pathways following excitation. The two pump photon energies, 1 = 21.0 eV and 21.6 eV, correspond to excitation into the low-energy shoulder on the main absorption band and to the maximum of this absorption band, respectively. The droplet excitation at 1 = 21.0 eV correlates with the 1s2s1S atomic state,7,9 which is inaccessible by electric dipole transition in the atom. It is more likely to occur near the surface of the droplet where the inversion symmetry is broken.19 Excitation into the main absorption peak (1 = 21.6 eV) can mostly be associated with the 1s2p1P atomic state.7,9 Additionally, Fig. 2(b) and (d) display electron spectra at longer delays when the fast relaxation dynamics have subsided. For these spectra, as well as the results in Fig. 4, we have subtracted the background to better observe the dynamic features. Specifically, we subtract, for each time delay, the electron spectrum produced by the FEL pulse alone, i.e. in the absence of the UV pulse.


image file: d2cp03335f-f2.tif
Fig. 2 Time-resolved photoelectron spectra of He droplets measured with XUV-pump and UV-probe pulses at photon energies 1 = 21.0 eV (a and b) and 1 = 21.6 eV (c and d). The maps (a and c) are composed of 28 and 31 spectra at different pump–probe delays, respectively. The white dashed lines indicate the expected electron energies for direct 1 + 1′ photoionization. Panels (b and d) show the electron spectra at fixed delays of 1.7 ps and 1.1 ps, respectively. The black horizontal lines indicate the expected electron energies for photoionization of specific excited states of quasi-free He* atoms and He2* excimers. The gray line is the multi-Gaussian fit to the experimental spectra. See the main text for additional details. Panel (e) shows the He droplet absorption spectrum from Joppien et al.9 and the energy levels populated in the course of relaxation. The red and blue arrows represent the pump and probe pulses, respectively, and the green arrows represent the relaxation channels.

We observe some similar features for the two photon energies. Near zero delay, the electron energy is at its highest and matches the combined photon energies of the XUV-pump and UV-probe pulses reduced by the ionization energy of He, 1 + 2 − IPHe = 1.1 eV for 1 = 21.0 eV (white dashed line in Fig. 2a) and 1.7 eV for 1 = 21.6 eV (white dashed line in c). For positive delays we observe a fast relaxation to lower energies. For 1 = 21.0 eV, the relaxation is so fast that it appears as an almost vertical band, meaning relaxation proceeds within the experimental time resolution (200 fs). For 1 = 21.6 eV, the relaxation is slower resulting in a clearly visible shift from 1.7 to 1.3 eV over about 300 fs. This shift is attributed to intraband electronic relaxation towards the lower edge of the droplet band correlating to the 1s2p1P state of He. Interestingly, we see a significant narrowing of this band at delays ≳200 fs converging toward a sharp peak whose position matches the 1s2p1P He atomic state, see Fig. 2(d). Simultaneously, new sharp features appear at lower electron kinetic energies. The most prominent one is an intense, nearly horizontal band at electron kinetic energy eKE ≈ 0.6 eV in Fig. 2(c). It was observed in our previous work19 and was attributed to relaxation into the 1s2s1S droplet state. This relaxation is likely mediated by a crossing of potential energy curves of the He2* excimer.41 It is accompanied by the formation of a void bubble around the localized He* excitation.42 At later times, this bubble migrates to the droplet surface and releases an excited He* atom or a He2* excimer which either remains weakly bound to the droplet surface or fully detaches into the vacuum.13 These species have previously been detected by static fluorescence measurements.22

Besides the 1s2s1S atomic state, additional features are clearly visible, which were not observed in our prior work.19 To gain a better understanding of the origin of these features we offset the vertical right axes in Fig. 2(a) and (c) to correspond to the excitation energy of the system. In this way, we can visually assign the two horizontal bands at higher kinetic energies to the 1s2p3P and 1s2p1P atomic states, respectively. Note that in a different experiment where we studied the relaxation dynamics in He nanodroplets excited to the 1s3p band,20 we did observe a shoulder on the 1s2s1S atomic peak, which we now fully resolve and attribute to the 1s2p3P atomic state. For the static measurements, shown in Fig. 2(b) and (d), horizontal lines are added to visualize the different excited states. These peaks are most prominent for 1 = 21.6 eV. For 21.0 eV, the lower intensity of the 1s2p1P atomic state could be due to its nominal excitation energy being slightly higher than the FEL photon energy.

We additionally observe another feature at lower kinetic energy, which we attribute to the 3Σ+u state of the He2* excimer photoionized into the lowest vibrational state v = 0 of the He2+ dimer ion. For photoionization into the first excited vibrational state v = 1 of He2+ we expect the electron energy to be eKE ≈ 0.17 eV, which roughly matches the position of the asymmetric feature peaked around eKE = 0.1 eV in the electron spectrum at 1 = 21.6 eV excitation. This feature may also have contributions from autoionization of superexcited He droplets;43 quasi-bound states ≲1.1 eV above the ionization energy of atomic He, Ei = 24.59 eV, or states bound by ≲1.6 eV below Ei, populated by ‘re-excitation’ out of the dark state 1s2s3S, can decay by emission of low-energy electrons.8,20 However, we refrain from a conclusive assignment of this feature.

To quantify the time evolution of the spectral features shown in Fig. 2, we fitted the delay-dependent photoelectron spectra with a multi-Gaussian function. The smooth thin lines show the best fits in Fig. 2(b) and (d). The resulting peak areas of each Gaussian component are depicted in Fig. 3(a) and (b) for the two photon energies 1 = 21.0 and 21.6 eV, respectively. In total, the fit function at 21.0 eV is the sum of 6 individual Gaussian functions to account for 3 sharp atomic lines (1s2s1S, 1s2p3P, 1s2p1P), two He2* molecular lines (v = 0, v = 1), and one broadened droplet spectral feature present in the photoelectron spectra shown in Fig. 2(b). At 1 = 21.6 eV (Fig. 2(d)), one additional peak accounting for near-zero kinetic energy electrons generated by autoionization is added to the fit. For the Gaussian functions representing the atomic components, the energy fit parameters are constrained to narrow intervals around their well-known excitation energies. Since the photon energy of the probe pulse is insufficient for photoionizing the lowest excited atomic state of He, 1s2s3S, by one-photon absorption, this state remains undetected in the present study.


image file: d2cp03335f-f3.tif
Fig. 3 Results of peak fits to the delay-dependent electron spectra recorded at pump photon energies 1 = 21.0 eV (panel a) and 2 = 21.6 eV (panel b). Shown are the integrated areas of Gaussian peaks fitted to the individual components in the electron spectra, labeled in Fig. 2. The largest ‘Droplet’ peak corresponds to the bright feature dominating the electron spectra around zero delay. The smooth lines depict fits of the data with eqn (3) and (4) from which characteristic time constants for the appearance of the individual components are inferred.

The resulting peak areas are plotted in Fig. 3 and fitted with the following simple model functions. For the droplet feature, we assume the model function

 
f(t) = Θ(tt0)[A + B[thin space (1/6-em)]exp(−(tt0)/τ)],(1)
where A and B are adjustable constants, Θ is the Heaviside step function, and τ is the decay time constant of the signal at positive delays, t > t0. For all other peaks, the model is
 
f(t) = Θ(tt0)C[1 − exp(−(tt0)/τ)].(2)

These functions are convoluted with a Gaussian cross-correlation function with a FWHM of 200 fs to account for the finite duration of the pump and probe pulses. This convolution can be carried out analytically and the resulting formulas are given explicitly in the appendix (eqn (3) and (4)). The best fits of the resulting fit functions to the experimental data are shown in Fig. 2(b) and (d) as smooth gray lines. At 1 = 21.0 eV, the 1s2p1P atomic state signal is small for all delays and is therefore omitted in the figure. At 1 = 21.6 eV (Fig. 2(b)), the 1s2p1P atomic state coincides with the lower edge of the broad droplet feature and is thus treated as a combined effect from both processes. The local maximum of the 1s2p3P peak area around 0.3 ps delay is due to a crosstalk from the dominant 1s2p1P droplet signal to the 1s2p3P signal at delays <0.6 ps. As such, these data points are excluded from the fit of the 1s2p3P signal.

From these fits we infer the characteristic time constants for each component which are summarized in Table 1, where droplet and atomic states are indicated with the superscripts (D) and (A), respectively. Overall, we observe both similarities and clear differences in the time evolution of the individual components in the electron spectra at the two photon energies 1 = 21.0 eV and 21.6 eV. First, the excited droplet state is populated within the range of overlapping pump and probe pulses and subsequently decays within a few hundred fs. At 1 = 21.0 eV, the droplet excitation, which correlates to the 1s2s1S state, reaches its maximum nearly at the zero-point of the delay (tmax ≈ 20 fs). Subsequently, it rapidly decays to zero (τ ≈ 280 fs) and evolves mostly into a well-separated narrow atomic state of the same symmetry (1s2s1S), cf.Fig. 2(a), which has a corresponding rise time τ ≈ 210 fs. These dynamics can be compared to time-dependent density functional (TDDFT) simulations of the dynamics of a 1s2s1S excited He* atom inside superfluid He.20 In these simulations, a void bubble forms around the He* whose radius was found to increase with an exponential time constant of 0.25 ps, in good agreement with the experimental values measured here. In our previous experiment where He droplets were resonantly excited into the higher-lying droplet absorption band at 1 = 23.7 eV, relaxation out of the 1s2s1S droplet excitation into the 1s2s1S atomic state was observed as a secondary process in a step-wise relaxation starting from a broad excitation correlating to the He 1s3p and 1s4p droplet states.20 There, a longer relaxation time ≈780 fs was measured. We ascribe the slower relaxation measured in those experiments to a weaker coupling of the 1s2s1S excitation to the surrounding He* when this state is populated indirectly by relaxation out of a higher level. In other words, as a high-lying state is excited, a bubble opens up around the He* over the course of electronic relaxation, thereby reducing the orbital overlap of the He* with surrounding He atoms.

Table 1 Time constants obtained by fitting peaks in the delay-dependent experimental electron spectra with a multi-Gaussian model function (third column). The first column labels the individual components in the electron spectra measured at 1 = 21.0 eV and 21.6 eV (second column). The reference values in the fourth column are from the literature
Transition 1 (eV) τ (ps) τ ref (ps) (1 = 23.8 eV)
Decay of 1S(D) 21.0 0.28 ± 0.02 0.78 ± 0.1620
1S(D)3P(A) 21.0 0.64 ± 0.17
1S(D)1S(A) 21.0 0.21 ± 0.02
1S(D)3Σ+u 21.0 0.19 ± 0.01
Decay of 1P(D)/(A) 21.6 0.27 ± 0.03 0.59 ± 0.06,20 ∼0.4521
1P(D)3P(A) 21.6 1.07 ± 0.45
1P(D)1S(A) 21.6 0.34 ± 0.03
1P(D)3Σ+u 21.6 0.73 ± 0.06


At 1 = 21.6 eV, where the He 1s2p1P droplet state is excited, the maximum intensity of the corresponding photoelectron signal is reached at a later delay time tmax = 150 fs before dropping with a similar decay time τ ≈ 270 fs as compared to the dynamics at 1 = 21.0 eV. Subsequently, the droplet state evolves continuously into the atomic 1s2p1P state which appears as a narrow band at the lower edge of the droplet feature in Fig. 2(c). The most abundant atomic state populated by relaxation is 1s2s1S with a corresponding rise time τ ≈ 340 fs. Relaxation into this state requires a decrease of the electron energy by up to 1 eV and a change of symmetry of the electronic state from 1s2p1P to 1s2s1S. This explains the wide gap between these two bright features in the time-resolved electron spectra, see Fig. 2(c), and the slower relaxation dynamics as compared to the direct excitation of the 1s2s1S droplet state at 1 = 21.0 eV Fig. 2(a). In previous experiments where the higher 1s3p/1s4p droplet band was excited, relaxation from the 1s2p1P droplet excitation to the 1s2s1S atomic state was measured to be slower, τ ≈ 590 fs.21 In another experiment using XUV pulses generated by HHG at nearly the same photon energy, τ ≈ 450 fs was measured.21 This confirms the trend that secondary steps within a relaxation cascade are slower, as the coupling of the excited state to the droplet weakens over the course of electronic relaxation and simultaneous bubble formation.

The emerging atomic-like triplet states 1s2p3P and the excimer state 3Σ+u have considerably longer rise times, see Table 1. Only the 3Σ+u state at 1 = 21.0 eV appears to be populated at the same rate as the initial relaxation of the 1s2s1S state. This may be due to its strong overlap with the 1s2s1S state in the electron spectra which introduces a large uncertainty in the peak fitting procedure. The large uncertainty of the 1s2p3P time constant at both photon energies is due to the spectral overlap of this component with the low-energy tail of the broad droplet feature at short time delays. The slower appearance time of triplet states is likely related to the fact that a spin flip is needed since optically excited states have singlet symmetry. Previous fluorescence measurements have suggested that triplet states are only formed by electron–ion recombination when exciting into the autoionizing states at > 23 eV.22 However, indications for the population of triplet states by relaxation even out of the lower-lying 1s2s and 1s2p states have been found; Penning ionization of alkali metal atoms attached to He nanodroplets occurred from both the He 1s2s1S state and, to a small extent, the 1s2s3S state, after excitation into the lower absorption band of the droplet at = 21.6 eV.23 Spin quenching was also observed for barium atoms attached to the surface of argon clusters.44 Thus, we conclude that He droplets are capable of inducing spin quenching of excited states to some degree as well, and the time scale is 600–1100 fs.

Additionally, we have performed similar measurements at 1 = 24.5 eV which is near the He atomic ionization threshold. The results are shown in Fig. 4. The most prominent feature in the electron spectra at this photon energy for all delays is a distribution of low-energy electrons which falls off exponentially in the range eKE = 0–2 eV (not shown). This electron distribution is generated by the XUV pump pulse and is due to the autoionization of He droplets excited above their adiabatic ionization energy.43,45 The low-energy part of this distribution appears as a negative signal in Fig. 4 since the pump-only spectrum is subtracted from all shown photoelectron spectra and the probe pulse efficiently depletes this channel by photoionizing the excited-state electrons into the continuum.


image file: d2cp03335f-f4.tif
Fig. 4 Time-resolved photoelectron spectra recorded at 1 = 24.5 eV composed of 13 spectra at different pump–probe delays. The dashed line indicates the expected electron energy for direct 1 + 1′ photoionization.

On top of this large signal, we detect photoelectrons emitted by 1 + 1′ photoionization with energies eKE = 4.5 eV. As it is a weak positive signal on a large background, the signal-to-noise ratio is low and experimental artifacts are seen in the spectra as horizontal stripes. The energy resolution may have been somewhat compromised having used an additional retarding voltage near the interaction region in the attempt to reduce the contribution of low-energy electrons in the spectra. Nevertheless, a few trends are clearly discernible. Similar to the time-dependent spectra recorded at 1 = 21.0 eV and 21.6 eV discussed above, as well as the ones previously recorded at 1 = 23.7 eV,20 the initially excited droplet state (bright red distribution around t = 0, eKE = 4.5 eV) relaxes into lower-lying states, in particular the 1s2s1S band around eKE ≈ 0.7 eV, within ∼150 fs. However, in this case, a substantial amount of population remains in higher-lying states at eKE > 3 eV. As these excited states continue to undergo autoionization, depletion of these states by the probe pulse quenches the autoionization process in the entire shown range of delays (negative signal near eKE = 0). In our previous measurement at 1 = 23.7 eV,20 depletion of the autoionization signal was only visible in the range of overlapping pump and probe pulses around t = 0, in line with the observation of fast relaxation of electron energies to states with eKE < 3 eV. Such a fast reappearance time of the autoionization signal following depletion of excitation at 1 = 23.7 eV was also seen in experiments using XUV pulses generated by HHG.10,46 Thus, we conclude that droplet excitation at 1 > 23 eV (correlating to He atomic states with principal quantum numbers n ≥ 3)8 relax more slowly if populated by prior relaxation from even higher excited states. It is likely that bubble formation initiated by the electronic excitation reduces the coupling of the excited state to the droplet over the course of relaxation, thereby gradually slowing down any subsequent electronic relaxation steps.

4 Conclusions

In summary, we have performed time-resolved, high-resolution photoelectron spectroscopy of resonantly excited He nanodroplets. In the XUV-pump and UV-probe electron spectra, we observed various distinct features due to two-color photoionization, ATI, and autoionization of multiply excited droplets by ICD. In particular, thanks to the improved resolution and a higher signal-to-noise ratio compared to previous measurements of electron spectra, we could observe new, ultrafast relaxation channels of the excited droplets into triplet atomic states, indicating efficient droplet-induced spin-relaxation, as well as the formation of the first excimer state of He2*. The convergence of these features to atomic-level energies on a timescale of a few 100 fs indicates ultrafast localization of the excitation at quasi-free He* atoms and He2* excimers residing in void bubbles or ejected from the droplets. From multi-Gaussian peak fits of the delay-dependent spectra, we inferred the relaxation-time constants for the individual final states. Significantly different relaxation dynamics were seen for excitation of the two lowest optically accessible He droplet states 1s2s1S and 1s2p1P, despite them overlapping in the static absorption spectrum.9 In both cases, we find that excited He droplets partially relax into previously unobserved triplet states, indicating efficient He droplet-induced spin relaxation.

Thus, He nanodroplets were found to turn into highly reactive and dissipative systems upon resonant excitation, featuring complex ultrafast relaxation dynamics including electronic-state hopping and spin flipping, the formation of bubbles and excimers, and few-body autoionization. Future time-resolved high-resolution photoelectron spectroscopic measurements of other types of pure or doped noble-gas clusters and nanodroplets will shed more light onto the peculiar properties of superfluid He nanodroplets on the one hand, and on more general aspects of the relaxation dynamics of nanoparticles irradiated by resonant UV or XUV radiation on the other hand.

Author contributions

A. C. L., B. B., M. B., C. C., M. D. F., S. H., S. R. K., A. N., N. P., O. P., P. R., N. B., and M. M. performed the experiments. J. D. A., S. D., L. G., M. F. K., D. M., S. M., K. C. P., E. S., and F. S. aided remotely in the analysis and interpretation of the experimental results. J. D. A. analyzed the data. A. C. L., J. D. A., and M. M. wrote the manuscript with input from all the co-authors.

Conflicts of interest

The authors declare no conflicting interests.

Appendix

The fit functions used for modelling the experimental data shown in Fig. 1 are given by the convolution of the model functions eqn (1) and (2) with a Gaussian function that accounts for the cross-correlation of the pump and probe pulses,
image file: d2cp03335f-t1.tif
The standard deviation σ is related to the FWHM, which is 200 fs, by image file: d2cp03335f-t2.tif. For model function (1) the convolution can be written out as
 
image file: d2cp03335f-t3.tif(3)
For model function (2) it is
 
image file: d2cp03335f-t4.tif(4)
Here, image file: d2cp03335f-t5.tif denotes the error function and erfc(z) = 1 − erf(z).

Acknowledgements

We gratefully acknowledge the support of the Chemical Sciences, Geosciences and Biosciences Division, Office of Basic Energy Sciences, Office of Science, U.S. Department of Energy, Grant No. DESC0012376 and DESC0063, the German Science Foundation (DFG) through project STI 125/19-2, and the Danish Agency for Science, Technology, and Innovation through the instrument center DanScatt. A. C. L. gratefully acknowledges the support of Charles and Peggy Hardin. M. F. K. is grateful for support by the Max Planck Society via the Max Planck Fellow program. S. R. K. gratefully acknowledges the financial support from the Min. of Education, Govt. of India through the scheme for promotion of research and academic collaboration and the Institute of Eminence programmes, from the Core Research Grant, Indo-Elettra and India-DESY schemes of the Dept. Of Science and Technology, Govt. Of India and from the Max Planck Society, Germany. M. M. is grateful for financial support from the Carlsberg Foundation.

References

  1. S. Grebenev, J. P. Toennies and A. F. Vilesov, Science, 1998, 279, 2083–2085 CrossRef CAS PubMed.
  2. M. Hartmann, F. Mielke, J. P. Toennies, A. F. Vilesov and G. Benedek, Phys. Rev. Lett., 1996, 76, 4560–4563 CrossRef CAS PubMed.
  3. J. Tang, Y. Xu, A. R. W. McKeller and W. Jäger, Science, 2002, 297, 2030–2033 CrossRef CAS PubMed.
  4. N. B. Brauer, S. Smolarek, E. Loginov, D. Mateo, A. Hernando, M. Pi, M. Barranco, W. J. Buma and M. Drabbels, Phys. Rev. Lett., 2013, 111, 153002 CrossRef PubMed.
  5. M. Mudrich and F. Stienkemeier, Int. Rev. Phys. Chem., 2014, 33, 301–339 Search PubMed.
  6. K. von Haeften, T. Laarmann, H. Wabnitz and T. Möller, Phys. Rev. Lett., 2001, 87, 153403 CrossRef CAS PubMed.
  7. K. D. Closser and M. Head-Gordon, J. Phys. Chem. A, 2010, 114, 8023–8032 CrossRef CAS PubMed.
  8. O. Kornilov, O. Bünermann, D. J. Haxton, S. R. Leone, D. M. Neumark and O. Gessner, J. Phys. Chem. A, 2011, 115, 7891–7900 CrossRef CAS PubMed.
  9. M. Joppien, R. Karnbach and T. Möller, Phys. Rev. Lett., 1993, 71, 2654–2657 CrossRef CAS PubMed.
  10. M. P. Ziemkiewicz, D. M. Neumark and O. Gessner, Int. Rev. Phys. Chem., 2015, 34, 239–267 Search PubMed.
  11. K. D. Closser, O. Gessner and M. Head-Gordon, J. Chem. Phys., 2014, 140, 134306 CrossRef PubMed.
  12. J. Keto, M. Stockton and W. Fitzsimmons, Phys. Rev. Lett., 1972, 28, 792 CrossRef CAS.
  13. H. Buchenau, J. P. Toennies and J. A. Northby, J. Chem. Phys., 1991, 95, 8134–8148 CrossRef CAS.
  14. J. Jortner, N. R. Kestner, S. A. Rice and M. H. Cohen, J. Chem. Phys., 1965, 43, 2614–2625 CrossRef CAS.
  15. C. C. Wang, O. Kornilov, O. Gessner, J. H. Kim, D. S. Peterka and D. M. Neumark, J. Phys. Chem. A, 2008, 112, 9356–9365 CrossRef CAS PubMed.
  16. A. Borghesani, Ions and Electrons in Liquid Helium, Oxford University Press, Oxford, 2007 Search PubMed.
  17. O. Bünermann, O. Kornilov, D. J. Haxton, S. R. Leone, D. M. Neumark and O. Gessner, J. Chem. Phys., 2012, 137, 214302 CrossRef PubMed.
  18. D. S. Peterka, A. Lindinger, L. Poisson, M. Ahmed and D. M. Neumark, Phys. Rev. Lett., 2003, 91, 043401 CrossRef PubMed.
  19. M. Mudrich, A. C. LaForge, A. Ciavardini, P. O'Keeffe, C. Callegari, M. Coreno, A. Demidovich, M. Devetta, M. D. Fraia and M. Drabbels, et al. , Nat. Commun., 2020, 11, 1–7 CrossRef PubMed.
  20. J. D. Asmussen, R. Michiels, K. Dulitz, A. Ngai, U. Bangert, M. Barranco, M. Binz, L. Bruder, M. Danailov and M. Di Fraia, et al. , Phys. Chem. Chem. Phys., 2021, 23, 15138–15149 RSC.
  21. M. P. Ziemkiewicz, C. Bacellar, K. R. Siefermann, S. R. Leone, D. M. Neumark and O. Gessner, J. Chem. Phys., 2014, 141, 1–9 CrossRef PubMed.
  22. K. von Haeften, A. R. de Castro, M. Joppien, L. Moussavizadeh, R. von Pietrowski and T. Möller, Phys. Rev. Lett., 1997, 78, 4371–4374 CrossRef CAS.
  23. L. Ben Ltaief, M. Shcherbinin, S. Mandal, S. R. Krishnan, A. C. Laforge, R. Richter, S. Turchini, N. Zema, T. Pfeifer and E. Fasshauer, et al. , J. Phys. Chem. Lett., 2019, 10, 6904–6909 CrossRef CAS PubMed.
  24. V. Lyamayev, Y. Ovcharenko, R. Katzy, M. Devetta, L. Bruder, A. Laforge, M. Mudrich, U. Person, F. Stienkemeier and M. Krikunova, et al. , J. Phys. B: At., Mol. Opt. Phys., 2013, 46, 164007 CrossRef.
  25. E. Allaria, R. Appio, L. Badano, W. A. Barletta, S. Bassanese, S. G. Biedron, A. Borga, E. Busetto, D. Castronovo and P. Cinquegrana, et al. , Nat. Photonics, 2012, 6, 699–704 CrossRef CAS.
  26. E. Allaria, A. Battistoni, F. Bencivenga, R. Borghes, C. Callegari, F. Capotondi, D. Castronovo, P. Cinquegrana, D. Cocco and M. Coreno, et al. , New J. Phys., 2012, 14, 1–19 CrossRef.
  27. P. Finetti, H. Höppner, E. Allaria, C. Callegari, F. Capotondi, P. Cinquegrana, M. Coreno, R. Cucini, M. B. Danailov and A. Demidovich, et al. , Phys. Rev. X, 2017, 7, 021043 Search PubMed.
  28. Y. Ovcharenko, V. Lyamayev, R. Katzy, M. Devetta, A. LaForge, P. O'Keeffe, O. Plekan, P. Finetti, M. Di Fraia and M. Mudrich, et al. , Phys. Rev. Lett., 2014, 112, 073401 CrossRef CAS PubMed.
  29. A. C. LaForge, M. Drabbels, N. B. Brauer, M. Coreno, M. Devetta, M. Di Fraia, P. Finetti, C. Grazioli, R. Katzy and V. Lyamayev, et al. , Sci. Rep., 2014, 4, 3621 CrossRef CAS PubMed.
  30. Y. Ovcharenko, A. C. LaForge, B. Langbehn, O. Plekan, R. Cucini, P. Finetti, P. O'Keeffe, D. Iablonskyi, T. Nishiyama and K. Ueda, et al. , New J. Phys., 2020, 22, 083043 CrossRef CAS.
  31. C. Svetina, C. Grazioli, N. Mahne, L. Raimondi, C. Fava, M. Zangrando, S. Gerusina, M. Alagia, L. Avaldi and G. Cautero, et al. , J. Synchrotron Radiat., 2015, 22, 538–543 CrossRef PubMed.
  32. P. Finetti, A. Demidovich, O. Plekan, M. Di Fraia, R. Cucini, C. Callegari, P. Cinquegrana, P. Sigalotti, R. Ivanov and M. B. Danailov, et al. , J. Opt., 2017, 19, 114010 CrossRef.
  33. J. P. Toennies and A. F. Vilesov, Angew. Chem., Int. Ed., 2004, 43, 2622–2648 CrossRef CAS PubMed.
  34. L. F. Gomez, E. Loginov, R. Sliter and A. F. Vilesov, J. Chem. Phys., 2011, 135, 154201 CrossRef PubMed.
  35. R. Squibb, M. Sapunar, A. Ponzi, R. Richter, A. Kivimäki, O. Plekan, P. Finetti, N. Sisourat, V. Zhaunerchyk and T. Marchenko, et al. , Nat. Commun., 2018, 9, 1–7 CrossRef CAS PubMed.
  36. R. Michiels, M. Abu-samha, L. B. Madsen, M. Binz, U. Bangert, L. Bruder, R. Duim, A. Wituschek, A. C. LaForge and R. J. Squibb, Phys. Rev. Lett., 2021, 127, 093201 CrossRef CAS PubMed.
  37. L. S. Cederbaum, J. Zobeley and F. Tarantelli, Phys. Rev. Lett., 1997, 79, 4778–4781 CrossRef CAS.
  38. A. I. Kuleff, K. Gokhberg, S. Kopelke and L. S. Cederbaum, Phys. Rev. Lett., 2010, 105, 043004 CrossRef PubMed.
  39. T. Chang and T. Fang, Phys. Rev. A: At., Mol., Opt. Phys., 1995, 52, 2638 CrossRef CAS PubMed.
  40. A. C. LaForge, R. Michiels, Y. Ovcharenko, A. Ngai, J. Escartn, N. Berrah, C. Callegari, A. Clark, M. Coreno and R. Cucini, et al. , Phys. Rev. X, 2021, 11, 021011 CAS.
  41. S. L. Fiedler and J. Eloranta, J. Low Temp. Phys., 2014, 174, 269–283 CrossRef CAS.
  42. K. von Haeften, T. Laarmann, H. Wabnitz and T. Möller, Phys. Rev. Lett., 2002, 88, 4 CrossRef PubMed.
  43. D. S. Peterka, A. Lindinger, L. Poisson, M. Ahmed and D. M. Neumark, Phys. Rev. Lett., 2003, 91, 043401 CrossRef PubMed.
  44. S. Awali, M.-A. Gaveau, M. Briant, J.-M. Mestdagh, B. Soep, O. Gobert, R. Maksimenka and L. Poisson, Phys. Chem. Chem. Phys., 2016, 18, 32378–32386 RSC.
  45. D. S. Peterka, J. H. Kim, C. C. Wang, L. Poisson and D. M. Neumark, J. Phys. Chem. A, 2007, 111, 7449 CrossRef CAS PubMed.
  46. O. Kornilov, C. C. Wang, O. Bünermann, A. T. Healy, M. Leonard, C. Peng, S. R. Leone, D. M. Neumark and O. Gessner, J. Phys. Chem. A, 2010, 114, 1437–1445 CrossRef CAS PubMed.

Footnote

These authors contributed equally to this work.

This journal is © the Owner Societies 2022
Click here to see how this site uses Cookies. View our privacy policy here.