Time-resolved Photoelectron Spectra Research Articles

Ultrafast Decay Dynamics of the 21ππ* Electronic State of N-Methyl-2-pyridone.

The ultrafast decay dynamics of N-methyl-2-pyridone upon excitation in the near-ultraviolet range of 261.5-227.9 nm is investigated using the femtosecond time-resolved photoelectron spectroscopy method. Irradiation at 261.5 nm prepares N-methyl-2-pyridone molecules with high vibrational levels in the 11ππ* state. The radiation-less decay to the ground state via internal conversion is suggested to be the dominant channel for the 11ππ* state with large vibrational excess energy, which is revealed by a lifetime of 1.6 ± 0.2 ps. As the pump wavelength decreases, we found that irradiation at 238.5 and 227.9 nm results in the population of the 21ππ* state. This is in agreement with the assignment of the vapor-phase UV absorption bands of N-methyl-2-pyridone. On the basis of the detailed analysis of our measured time-resolved photoelectron spectra at all pump wavelengths, we conclude that the 21ππ* state has an ultrashort lifetime of 50 ± 10 fs. In addition, the S1(11ππ*) state is subsequently populated via internal conversion and decays over a lifetime of 680-620 fs. The most probable whole deactivation pathway of the 21ππ* state is discussed. This experimental study provides new insights into the excitation energy-dependent decay dynamics of electronically excited N-methyl-2-pyridone.

Insights into ultrafast decay dynamics of electronically excited pyridine-N-oxide.

The ultrafast decay dynamics of pyridine-N-oxide upon excitation in the near-ultraviolet range of 340.2-217.6 nm is investigated using the femtosecond time-resolved photoelectron imaging technique. The time-resolved photoelectron spectra and photoelectron angular distributions at all pump wavelengths are carefully analyzed and the following view is derived: at the longest pump wavelengths (340.2 and 325.6 nm), pyridine-N-oxide is excited to the S1(1ππ*) state with different vibrational levels. The depopulation rate of the S1 state shows a marked dependence on vibrational energy and mode, and the lifetime is in the range of 1.4-160 ps. At 289.8 and 280.5 nm, both the second 1ππ* state and the S1 state are initially prepared. The former has an extremely short lifetime of ∼60 fs, which indicates that the ultrafast deactivation pathway such as a rapid internal conversion to one close-lying state is its dominant decay channel, while the latter is at high levels of vibrational excitation and decays within the range of 380-520 fs. At the shortest pump wavelengths (227.3 and 217.6 nm), another excited state of Rydberg character is mostly excited. We assign this state to the 3s Rydberg state which has a lifetime of 0.55-2.2 ps. This study provides a comprehensive picture of the ultrafast excited-state decay dynamics of the photoexcited pyridine-N-oxide molecule.

Tiresia: A code for molecular electronic continuum states and photoionization

The Tiresia program [1] provides access to numerically accurate solutions of the one-particle Schrödinger equation for highly excited states of complex polyatomic molecules, both bound and continuum, that cannot be described by conventional Quantum Chemistry approaches. It is based on an expansion of the required solution in a local multicentric basis set, with primitive functions built as products of a radial B-spline times a real spherical harmonic. In conjunction with Density Functional Theory (DFT), it has been extensively employed in a large variety of photoionization studies, also for rather large systems. Highly excited bound states as well as wavepacket propagation can also be accurately described. In fact, the flexibility of the basis essentially allows accurate solutions of linear operator equations, like Poisson or inhomogeneous perturbative equations, which are employed in the code. The program is parallelized with standard MPI-I instructions and makes extensive use of the Scalapack linear algebra library. Ancillary programs are available for the evaluation of photoionization cross sections and angular distributions from randomly to fully oriented molecules. Program summaryProgram Title: TiresiaCPC Library link to program files:https://doi.org/10.17632/fcrjxwgjxh.1Licensing provisions: GPLv3Programming language: Fortran77, Fortran90, MPISupplementary material: Program manual documentNature of problem: Accurate solutions for highly excited and continuum electronic states in complex polyatomic molecules. Molecular photoionization cross sections and angular distributions under high energy, high-intensity radiation pulses from Synchrotron radiation and laser sources, photoelectron imaging in pump-probe experiments, basis for electronic wavepackets under ultrafast or nonperturbative excitation.Solution method: Solution of the Schrödinger and similar linear operator equations in a finite domain is obtained via basis set expansion. Flexible basis set, that may approach practical completeness within the domain, is obtained as a multicenter set of B-spline radial functions times spherical harmonics. Accurate numerical integration is employed for the evaluation of matrix elements, and conventional diagonalization for bound states, or Galerkin approach for the full multichannel solution in the continuum. Full hamiltonian and dipole matrices in the spectral basis are available for time propagation. DFT many-body description is available, and strong correlations in the bound states may be incorporated via Dyson orbitals.Additional comments including restrictions and unusual features: The code is noted for computational efficiency, which allows fast yet reasonably accurate photoionization calculations for medium-sized molecules, allowing, e.g., calculations at many molecular geometries as required to describe time-resolved photoelectron spectra in pump-probe experiments.

Photoionization cross sections and photoelectron angular distributions of molecules with XCHEM-2.0

The XCHEM code was introduced in 2017 [1] to provide an accurate description of electron correlation and exchange in the electronic continuum of molecules at the same level as complete or restricted active-space self-consistent field (CASSCF or RASSCF) methods. This has allowed for an accurate description of molecular photoionization in the region of Feshbach resonances, shake up processes in which ionization is accompanied by excitation of one or several of the remaining electrons, and interchannel couplings. The success of XCHEM for small molecules has led us to improve its performance in several aspects, which now allows for the description of resonant molecular photoionization in larger systems. In addition, we have incorporated the possibility to calculate photoelectron angular distributions in the laboratory and molecular frames, which are essential to interpret angularly resolved photoionization experiments. Here we show its performance in the N2 and pyrazine molecules. The new version of the code, XCHEM-2.0, is freely available at https://gitlab.com/xchem/xchem_public. Program summaryProgram Title: XCHEM-2.0CPC Library link to program files:https://doi.org/10.17632/t8tbk9gdt2.1Developer's repository link:https://gitlab.com/xchem/xchem_public/-/archive/2.0.0/xchem_public-2.0.0.tar.gzLicensing provisions: LGPLProgramming language: Fortran, pythonNature of problem: Understanding photoionization of molecules is a recurrent goal in Atomic, Molecular and Optical physics, even more nowadays with the advent of XUV and X-ray high-harmonic generation sources and free electron lasers. Photoionization can be employed to induce electron dynamics in molecules, e.g., by using ultrashort XUV or X-ray pulses, or as a probe of these dynamics by recording time-resolved photoelectron spectra. At low photoelectron energies, photoionization may result from various competing processes, such as direct ionization, autoionization of Feshbach and shape resonances, shake up, etc., each of these leaving its trace in the measured spectra. So, interpretation of these spectra can be a difficult task and theoretical simulations are often required to separate the different contributions. However, a method able to do so must be able to describe electron correlation in the electronic continuum beyond the Hartree-Fock (HF) and configuration interaction singles (CIS) approximations. XCHEM-2.0 has been designed to account for all these processes by providing an accurate description of electron correlation in the ionization continuum of molecules.Solution method: XCHEM-2.0 employs a close coupling formalism to describe the electronic continuum of molecules at the level of multiference configuration interaction (MRCI) methods in combination with a hybrid Gaussian-Bspline basis set (the so-called GABs basis). In particular, it makes use of a restricted active space self-consistent field (RASSCF) approximation to evaluate the bound electronic states of the remaining molecular cation using localized Gaussian functions. These localized Gaussians are then supplemented with a single centered Gaussian expansion to allow the electron to leave the nuclear environment. Finally, the basis set includes a set of B-spline functions, centered at the same position as the single-centered Gaussians, which reach the asymptotic region. Thanks to the last B-spline function, which does not vanish at the box boundary, XCHEM-2.0 can determine continuum states fulfilling any required asymptotic behavior. The present version of the code improves on the original version by incorporating a more efficient augmentation procedure to build Ne-electron configurations from (Ne−1) ones, reducing the space to store data, providing a more efficient removal of linear dependencies resulting from the over-completeness of the polycentric+GABS basis and a more efficient solution of the scattering equations, and including new routines to calculate photoelectron angular distributions in the laboratory and molecular frames.Additional comments including restrictions and unusual features: The proposed method focuses on RASSCF methodology, so part of electron correlation is still not included. Also, one- and two electron integrals involving exclusively Gaussian functions are done using a modified version of the OpenMOLCAS software, so the installation of this code is also required and limitation associated with the latter is applicable. The code is limited to include a few ionization channels (around 20) with medium angular momentum (limited to ℓ = 15 as in OpenMOLCAS).

Resolving Quantum Interference Black Box through Attosecond Photoionization Spectroscopy.

Multiphoton light-matter interactions invoke a so-called "black box" in which the experimental observations contain the quantum interference between multiple pathways. Here, we employ polarization-controlled attosecond photoelectron metrology with a partial wave manipulator to deduce the pathway interference within this quantum 'black box" for the two-photon ionization of neon atoms. The angle-dependent and attosecond time-resolved photoelectron spectra are measured across a broad energy range. Two-photon phase shifts for each partial wave are reconstructed through the comprehensive analysis of these photoelectron spectra. We resolve the quantum interference between the degenerate p→d→p and p→s→p two-photon ionization pathways, in agreement with our theoretical simulations. Our approach thus provides an attosecond time-resolved microscope to look inside the "black box" of pathway interference in ultrafast dynamics of atoms, molecules, and condensed matter.

Photodissociation of water molecule at short photon wavelengths: Dynamical studies

In our last study [J. Phys. B At. Mol. Opt. Phys. 54, 125,102 (2021).], we reported the ab initio calculation of the full-dimensional potential energy surfaces of water molecule including 9 A’ and 9 A” states in Cs symmetry. In this study, we performed additional non-adiabatic semi-classical studies based on the potential energy surfaces. Our simulation successfully repeated the near picosecond lifetime of the F∼1A′ state measured by time resolved photo-electron spectra experiment [Chinese J. Chem. Phys. 32, 53 (2019)]. We also determined the dissociation branching ratio including H + OH(X, A), H + H + O and H2+O channels. In addition, the reaction path corresponding to H2+O (1S) channel is clearly marked out, which is found in recent free-electron laser experiment [Nat. Commun. 12, 6,303 (2021)].

Ultrafast Ring Closure Reaction of Gaseous cis-Stilbene from S1(ππ*)

cis-Stilbene (cis-St) is a well-known benchmark system for cis-trans photoisomerization. cis-St also produces 4a,4b-dihydrophenanthrene (DHP) in solution with a quantum yield of less than 0.19. The ring closure reaction, however, has never been identified for gaseous cis-St, and a recent computational simulation predicted the quantum yield of DHP to be only 0.04. In the present study, we identified an ultrafast ring closure reaction of gaseous cis-St for the first time using extreme ultraviolet time-resolved photoelectron spectroscopy. Surface hopping trajectory calculations at the SA3-XMS-CASPT2(2,2) level of theory reproduce the features of the observed time-resolved photoelectron spectra and predict the cis-St:DHP:trans-St branching ratio to be 0.55:0.41:0.04, in contrast with previous estimates. The results indicate that photoexcited cis-St favors ring closure over cis-trans isomerization under the isolated condition.

Nuclear-Electron Correlation Effects and Their Photoelectron Imprint in Molecular XUV Ionisation.

The ionisation of molecules by attosecond XUV pulses is accompanied by complex correlated dynamics, such as the creation of coherent electron wave packets in the parent ion, their interplay with nuclear wave packets, and a correlated photoelectron moving in a multi-centred potential. Additionally, these processes are influenced by the dynamics prior to and during the ionisation. To fully understand and subsequently control the ionisation process on different time scales, a profound understanding of electron and nuclear correlation is needed. Here, we investigate the effect of nuclear–electron correlation in a correlated two-electron and one-nucleus quantum model system. Solving the time-dependent Schrödinger equation allows to monitor the correlation impact pre, during, and post-XUV ionisation. We show how an initial nuclear wave packet displaced from equilibrium influences the post-ionisation dynamics by means of momentum conservation between the target and parent ion, whilst the attosecond electron population remains largely unaffected. We calculate time-resolved photoelectron spectra and their asymmetries and demonstrate how the coupled electron–nuclear dynamics are imprinted on top of electron–electron correlation on the photoelectron properties. Finally, our findings give guidelines towards when correlation resulting effects have to be incorporated and in which instances the exact correlation treatment can be neglected.

Data analysis procedures for time-resolved x-ray photoelectron spectroscopy at a SASE free-electron-laser

The random nature of self-amplified spontaneous emission (SASE) is a well-known challenge for x-ray core level spectroscopy at SASE free-electron lasers (FELs). Especially in time-resolved experiments that require a combination of good temporal and spectral resolution the jitter and drifts in the spectral characteristics, relative arrival time as well as power fluctuations can smear out spectral-temporal features. We present a combination of methods for the analysis of time-resolved photoelectron spectra based on power and time corrections as well as self-referencing of a strong photoelectron line. Based on sulfur 2p photoelectron spectra of 2-thiouracil taken at the SASE FEL FLASH2, we show that it is possible to correct for some of the photon energy drift and jitter even when reliable shot-to-shot photon energy data is not available. The quality of pump–probe difference spectra improves as random jumps in energy between delay points reduce significantly. The data analysis allows to identify coherent oscillations of 1 eV shift on the mean photoelectron line of 4 eV width with an error of less than 0.1 eV.

Intersystem crossing and internal conversion dynamics with GAIMS-TeraChem: Excited state relaxation in 2-cyclopentenone.

Excited states relaxation in complex molecules often involves two types of nonradiative transitions, internal conversion (IC) and intersystem crossing (ISC). In the situations when the timescales of IC and ISC are comparable, an interplay between these two types of transitions can lead to complex nonadiabatic dynamics on multiple electronic states of different characters and spin multiplicities. We demonstrate that the generalized ab initio multiple spawning (GAIMS) method interfaced with the fast graphics processing unit-based TeraChem electronic structure code can be used to model such nonadiabatic dynamics involving both the IC and ISC transitions in molecules of moderate size. We carried out 1500fs GAIMS simulations leading to the creation of up to 2500 trajectory basis functions to study the excited states relaxation in 2-cyclopentenone. After a vertical excitation from the ground state to the bright S2 state, the molecule quickly relaxes to the S1 state via conical intersection. The following relaxation proceeds along two competing pathways: one involves IC to the ground state, and the other is dominated by ISC to the low-lying triplet states. The time constants describing the population transfer between the six lowest singlet and triplet states predicted by the GAIMS dynamics are in good agreement with the characteristic times of IC and ISC obtained from the analysis of the time-resolved photoelectron spectrum.

Combined Surface-Hopping, Dyson Orbital, and B-Spline Approach for the Computation of Time-Resolved Photoelectron Spectroscopy Signals: The Internal Conversion in Pyrazine.

A computational protocol for simulating time-resolved photoelectron signals of medium-sized molecules is presented. The procedure is based on a trajectory surface-hopping description of the excited-state dynamics and a combined Dyson orbital and multicenter B-spline approach for the computation of cross sections and asymmetry parameters. The accuracy of the procedure has been illustrated for the case of ultrafast internal conversion of gas-phase pyrazine excited to the 1B2u(ππ*) state. The simulated spectra and the asymmetry map are compared to the experimental data, and a very good agreement was obtained without applying any energy-dependent rescaling or broadening. An interesting side result of this work is the finding that the signature of the 1Au(nπ*) state is indistinguishable from that of the 1B3u(nπ*) state in the time-resolved photoelectron spectrum. By locating four symmetrically equivalent minima on the lowest-excited (S1) adiabatic potential energy surface of pyrazine, we revealed the strong vibronic coupling of the 1Au(nπ*) and 1B3u(nπ*) states near the S1 ← S0 band origin.

Identification of an ultrafast internal conversion pathway of pyrazine by time-resolved vacuum ultraviolet photoelectron spectrum simulations.

The internal conversion from the optically bright S2 (1B2u, ππ*) state to the dark S1 (1B3u, nπ*) state in pyrazine is a standard benchmark for experimental and theoretical studies on ultrafast radiationless decay. Since 2008, a few theoretical groups have suggested significant contributions of other dark states S3 (1Au, nπ*) and S4 (1B2g, nπ*) to the decay of S2. We have previously reported the results of nuclear wave packet simulations [Kanno et al., Phys. Chem. Chem. Phys. 17, 2012 (2015)] and photoelectron spectrum calculations [Mignolet et al., Chem. Phys. 515, 704 (2018)] that support the conventional two-state picture. In this article, the two different approaches, i.e., wave packet simulation and photoelectron spectrum calculation, are combined: We computed the time-resolved vacuum ultraviolet photoelectron spectrum and photoelectron angular distribution for the ionization of the wave packet transferred from S2 to S1. The present results reproduce almost all the characteristic features of the corresponding experimental time-resolved spectrum [Horio et al., J. Chem. Phys. 145, 044306 (2016)], such as a rapid change from a three-band to two-band structure. This further supports the existence and character of the widely accepted pathway (S2 → S1) of ultrafast internal conversion in pyrazine.

Time Resolved Photoelectron Spectroscopy as a Test of Electronic Structure and Nonadiabatic Dynamics.

We compare different levels of theory for simulating excited state molecular dynamics and use time-resolved photoelectron spectroscopy measurements to benchmark the theory. We perform trajectory surface hopping simulations for uracil excited to the first bright state (ππ*) using three different levels of theory (CASSCF, MRCIS, and XMS-CASPT2) in order to understand the role of dynamical correlation in determining the excited state dynamics, with a focus on the coupling between different electronic states and internal conversion back to the ground state. These dynamics calculations are used to simulate the time-resolved photoelectron spectra. The comparison of the calculated and measured spectra allows us to draw conclusions regarding the relative insights and quantitative accuracy of the calculations at the three different levels of theory, demonstrating that detailed quantitative comparisons of time-resolved photoelectron spectra can be used to benchmark methodology.

Design and characterization of a magnetic bottle electron spectrometer for time-resolved extreme UV and X-ray photoemission spectroscopy of liquid microjets

We describe a magnetic bottle time-of-flight electron spectrometer designed for time-resolved photoemission spectroscopy of a liquid microjet using extreme UV and X-ray radiation. The spectrometer can be easily reconfigured depending on experimental requirements and the energy range of interest. To improve the energy resolution at high electron kinetic energy, a retarding potential can be applied either via a stack of electrodes or retarding mesh grids, and a flight-tube extension can be attached to increase the flight time. A gated electron detector was developed to reject intense parasitic signal from light scattered off the surface of the cylindrically shaped liquid microjet. This detector features a two-stage multiplication with a microchannel plate plus a fast-response scintillator followed by an image-intensified photon detector. The performance of the spectrometer was tested at SPring-8 and SACLA, and time-resolved photoelectron spectra were measured for an ultrafast charge transfer to solvent reaction in an aqueous NaI solution with a 200 nm UV pump pulses from a table-top ultrafast laser and the 5.5 keV hard X-ray probe pulses from SACLA.

Ultrafast Dynamics of the Isolated Adenosine-5'-triphosphate Dianion Probed by Time-Resolved Photoelectron Imaging.

The excited state dynamics of the doubly deprotonated dianion of adenosine-5'-triphosphate, [ATP-H2]2-, has been spectroscopically explored by time-resolved photoelectron spectroscopy following excitation at 4.66 eV. Time-resolved photoelectron spectra show that two competing processes occur for the initially populated 1ππ* state. The first is rapid electron emission by tunneling through a repulsive Coulomb barrier as the 1ππ* state is a resonance. The second is nuclear motion on the 1ππ* state surface leading to an intermediate that no longer tunnels and subsequently decays by internal conversion to the ground electronic state. The spectral signatures of the features are similar to those observed for other adenine-derivatives, suggesting that this nucleobase is quite insensitive to the nearby negative charges localized on the phosphates, except of course for the appearance of the additional electron tunneling channel, which is open in the dianion.

Signatures of Conical Intersection Dynamics in the Time-Resolved Photoelectron Spectrum of Furan: Theoretical Modeling with an Ensemble Density Functional Theory Method.

The non-adiabatic dynamics of furan excited in the ππ* state (S2 in the Franck–Condon geometry) was studied using non-adiabatic molecular dynamics simulations in connection with an ensemble density functional method. The time-resolved photoelectron spectra were theoretically simulated in a wide range of electron binding energies that covered the valence as well as the core electrons. The dynamics of the decay (rise) of the photoelectron signal were compared with the excited-state population dynamics. It was observed that the photoelectron signal decay parameters at certain electron binding energies displayed a good correlation with the events occurring during the excited-state dynamics. Thus, the time profile of the photoelectron intensity of the K-shell electrons of oxygen (decay constant of 34 ± 3 fs) showed a reasonable correlation with the time of passage through conical intersections with the ground state (47 ± 2 fs). The ground-state recovery constant of the photoelectron signal (121 ± 30 fs) was in good agreement with the theoretically obtained excited-state lifetime (93 ± 9 fs), as well as with the experimentally estimated recovery time constant (ca. 110 fs). Hence, it is proposed to complement the traditional TRPES observations with the trXPS (or trNEXAFS) measurements to obtain more reliable estimates of the most mechanistically important events during the excited-state dynamics.

Time-resolved photoelectron imaging of complex resonances in molecular nitrogen.

We have used the FERMI free-electron laser to perform time-resolved photoelectron imaging experiments on a complex group of resonances near 15.38 eV in the absorption spectrum of molecular nitrogen, N2, under jet-cooled conditions. The new data complement and extend the earlier work of Fushitani et al. [Opt. Express 27, 19702-19711 (2019)], who recorded time-resolved photoelectron spectra for this same group of resonances. Time-dependent oscillations are observed in both the photoelectron yields and the photoelectron angular distributions, providing insight into the interactions among the resonant intermediate states. In addition, for most states, we observe an exponential decay of the photoelectron yield that depends on the ionic final state. This observation can be rationalized by the different lifetimes for the intermediate states contributing to a particular ionization channel. Although there are nine resonances within the group, we show that by detecting individual photoelectron final states and their angular dependence, we can identify and differentiate quantum pathways within this complex system.

Ultrafast decay dynamics of electronically excited 2-ethylpyrrole.

The excited-state decay dynamics of 2-ethylpyrrole following UV excitation in the wavelength range of 254.8-218.0 nm is investigated in detail using the femtosecond time-resolved photoelectron imaging method. The time-resolved photoelectron spectra and photoelectron angular distributions at all pump wavelengths are carefully analysed and the following picture is derived: at the longest pump wavelengths (254.8, 248.3 and 246.1 nm), 2-ethylpyrrole is excited to the S1(1πσ*) state having a lifetime of about 50 fs. At 248.3, 246.1 and 237.4 nm, another excited state of Rydberg character is excited. The lifetime of this state is ∼570 fs at 237.4 nm and becomes slightly longer at other two pump wavelengths. At the shortest pump wavelengths (230.8 and 218.0 nm), 2-ethylpyrrole is excited to a state which is tentatively assigned to the 11ππ* state, having a lifetime of 75 ± 15 and 48 ± 10 fs for the longer and shorter pump wavelengths, respectively. Internal conversion to the S1(1πσ*) state might be one of the decay mechanisms of the 11ππ* state.

Molecular fragmentation as a way to reveal early electron dynamics induced by attosecond pulses.

We present a theoretical study of the electron and nuclear dynamics that would arise in an attosecond two-color XUV-pump/XUV-probe experiment in glycine. In this scheme, the broadband pump pulse suddenly ionizes the molecule and creates an electronic wave packet that subsequently evolves under the influence of the nuclear motion until it is finally probed by the second XUV pulse. To describe the different steps of such an experiment, we have combined a multi-reference static-exchange scattering method with a trajectory surface hopping approach. We show that by changing the central frequency of the pump pulse, i.e., by engineering the initial electronic wave packet with the pump pulse, one can drive the cation dynamics into a specific fragmentation pathway. Reminiscence of this early electron dynamics can be observed in specific fragmentation channels (not all of them) as a function of the pump-probe delay and in time-resolved photoelectron spectra at specific photoelectron energies. The optimum conditions to visualize the initial electronic coherence in photoelectron and photo-ion spectra depend very much on the characteristics of the pump pulse as well as on the electronic structure of the molecule under study.

Distinction of Electron Dispersion in Time-Resolved Photoemission Spectroscopy.

While recent experiments provided compelling evidence for an intricate dependence of attosecond photoemission-time delays on the solid's electronic band structure, the extent to which electronic transport and dispersion in solids can be imaged in time-resolved photoelectron (PE) spectra remains poorly understood. Emphasizing the distinction between photoemission time delays measured with two-photon, two-color interferometric spectroscopy, and transport times, we demonstrate how the effect of energy dispersion in the solid on photoemission delays can, in principle, be observed in interferometric photoemission. We reveal analytically a scaling relation between the PE transport time in the solid and the observable photoemission delay and confirm this relation in numerical simulations for a model system. We trace photoemission delays to the phase difference the PE accumulates inside the solid and, in particular, predict negative photoemission delays. Based on these findings, we suggest a novel time-domain interferometric solid-state energy-momentum-dispersion imaging method.