Time-resolved measurement of surface diffusion induced by femtosecond laser pulses

Quantifying Ultrafast Energy Transfer from Plasmonic Hot Carriers for Pulsed Photocatalysis on Nanostructures

Photocatalysis with plasmonic nanostructures has lately emerged as a transformative paradigm to drive and alter chemical reactions using light. At the surface of metallic nanoparticles, photoexcitation results in strong near fields, short-lived high-energy "hot" carriers, and light-induced heating, thus creating a local environment where reactions can occur with enhanced efficiencies. In this context, it is critical to understand how to manipulate the nonequilibrium processes triggered by light, as their ultrafast (femto- to picoseconds) relaxation dynamics compete with the process of energy transfer toward the reactants. Accurate predictions of the plasmon photocatalytic activity can lead to optimized nanophotonic architectures with enhanced selectivity and rates, operating beyond the intrinsic limitations of the steady state. Here, we report on an original modeling approach to quantify, with space, time, and energy resolution, the ultrafast energy exchange from plasmonic hot carriers (HCs) to molecular systems adsorbed on the metal nanoparticle surface while consistently accounting for photothermal bond activation. Our analysis, illustrated for a few typical cases, reveals that the most energetic nonequilibrium carriers (i.e., with energies well far from the Fermi level) may introduce a wavelength-dependence of the reaction rates, and it elucidates on the role of the carriers closer to the Fermi energy and the photothermally heated lattice, suggesting ways to enhance and optimize each contribution. We show that the overall reaction rates can benefit strongly from using pulsed illumination with the optimal pulse width determined by the properties of the system. Taken together, these results contribute to the rational design of nanoreactors for pulsed catalysis, which calls for predictive modeling of the ultrafast HC-hot adsorbate energy transfer.

Ultrafast, Non-Equilibrium and Transient Heating and Sintering of Nanocrystals for Nanoscale Metal Printing

The carrier excitation, relaxation, energy transport, and conversion processes during light-nanocrystal (NC) interactions have been intensively investigated for applications in optoelectronics, photocatalysis, and photovoltaics. However, there are limited studies on the non-equilibrium heating under relatively high laser excitation that leads to NCs sintering. Here, the authors use femtosecond laser two-pulse correlation and in-situ optical transmission probing to investigate the non-equilibrium heating of NCs and transient sintering dynamics. First, a two-pulse correlation study reveals that the sintering rate strongly increases when the two heating laser pulses are temporally separated by <10 ps. Second, the sintering rate is found to increase nonlinearly with laser fluence when heating with ≈700 fs laser pulses. By three-temperature modeling, the NC sintering mechanism mediated by electron induced ligand transformation is suggested. The ultrafast and non-equilibrium process facilitates sintering in dry (spin-coated) and wet (solvent suspended) environments. The nonlinear dependence of sintering rate on laser fluence is exploited to print sub-diffraction-limited features in NC suspension. The smallest feature printed is ≈200 nm, which is ≈¼ of the laser wavelength. These findings provide a new perspective toward nanomanufacturing development based on probing and engineering ultrafast transport phenomena in functional NCs.

On-Surface Chemistry Induced by Long-Lived Excitons: (NO)2 Dissociation on C60

Solid-state excitonic excitations play an increasingly important role in optoelectronic and light harvesting processes due to their ubiquitous presence in dipolar two-dimensional materials. Here we show that long-lived solid-state excitons induce chemical reactions in adsorbed molecules and thus convert light into chemical energy. For the model system (NO)₂ dimer adsorbed on ordered c(4×4) C₆₀ films, time-of-flight measurements following UV laser excitation reveal a slow and a fast dissociative desorption channel, which are assigned to intersystem crossing and internal conversion, respectively, by time-dependent density functional theory calculations.

A practical ansatz for evaluating the electronic friction tensor accurately, efficiently, and in a nearly black-box format

It is well-known that under conditions of fast electronic equilibration and weak nonadiabaticity, nonadiabatic effects induced by electron-hole pair excitations can be partly incorporated through a frictional force. However, ab initio computation of the electronic friction tensor suffers from numerical instability and usually demands a convergence check. In this study, we present an efficient and accurate interpolation method for computing the electronic friction tensor in a nearly black-box manner as appropriate for molecular dynamics. In almost all cases, our method agrees quite well with the exact friction tensor which is available for several quadratic Hamiltonians. As such, we outperform more conventional approaches that are based on the introduction of a broadening parameter. Future work will implement this interpolation approach within ab initio software packages.

Ultrafast Transient Dynamics of Adsorbates on Surfaces Deciphered: The Case of CO on Cu(100)

Time-resolved vibrational spectroscopy constitutes an invaluable experimental tool for monitoring hot-carrier-induced surface reactions. However, the absence of a full understanding of the precise microscopic mechanisms causing the transient spectral changes has limited its applicability. Here we introduce a robust first-principles theoretical framework that successfully explains both the nonthermal frequency and linewidth changes of the CO internal stretch mode on Cu(100) induced by femtosecond laser pulses. Two distinct processes engender the changes: electron-hole pair excitations underlie the nonthermal frequency shifts, while electron-mediated vibrational mode coupling gives rise to linewidth changes. Furthermore, the origin and precise sequence of coupling events are finally identified.

Understanding the Enhancement of Surface Diffusivity by Dimerization

Beyond dilute coverage, the collective diffusion of molecules might enhance material transport. We reveal an enhanced mobility of molecular dimers by separating two motions, diffusion and rotation, of CO dimers on elemental Ag(100) as well as on a dilute Cu alloy of Ag(100). From time-lapsed scanning tunneling microscopy movies recorded between 15 and 25 K, we determine the activation energy of dimer diffusion on elemental Ag(100) to be, at (40±2)  meV, considerably smaller than the one for monomer diffusion, at (72±1)  meV. The alloyed Cu atoms reduce the dimer mobility facilitating to determine their rotational barrier separately to be (39±3)  meV. Disentangling different degrees of freedom suggests that a rotational motion is at the origin of enhanced dimer diffusivity.

Adlayer structure dependent ultrafast desorption dynamics in carbon monoxide adsorbed on Pd (111)

We report our ultrafast photoinduced desorption investigation of the coverage dependence of substrate-adsorbate energy transfer in carbon monoxide adlayers on the (111) surface of palladium. As the CO coverage is increased, the adsorption site population shifts from all threefold hollows (up to 0.33 ML), to bridge and near bridge (>0.5 to 0.6 ML) and finally to mixed threefold hollow plus top site (at saturation at 0.75 ML). We show that between 0.24 and 0.75 ML this progression of binding site motifs is accompanied by two remarkable features in the ultrafast photoinduced desorption of the adsorbates: (i) the desorption probability increases roughly two orders magnitude, and (ii) the adsorbate-substrate energy transfer rate observed in two-pulse correlation experiments varies nonmonotonically, having a minimum at intermediate coverages. Simulations using a phenomenological model to describe the adsorbate-substrate energy transfer in terms of frictional coupling indicate that these features are consistent with an adsorption-site dependent electron-mediated energy coupling strength, ηel, that decreases with binding site in the order: three-fold hollow > bridge and near bridge > top site. This weakening of ηel largely counterbalances the decrease in the desorption activation energy that accompanies this progression of adsorption site motifs, moderating what would otherwise be a rise of several orders of magnitude in the desorption probability. Within this framework, the observed energy transfer rate enhancement at saturation coverage is due to interadsorbate energy transfer from the copopulation of molecules bound in three-fold hollows to their top-site neighbors.

Disentangling Multidimensional Nonequilibrium Dynamics of Adsorbates: CO Desorption from Cu(100)

Hot carriers at metal surfaces can drive nonthermal reactions of adsorbates. Characterizing nonequilibrium statistics among various degrees of freedom in an ultrafast time scale is crucial to understand and develop hot carrier-driven chemistry. Here we demonstrate multidimensional vibrational dynamics of carbon monoxide (CO) on Cu(100) along hot-carrier induced desorption studied by using time-resolved vibrational sum-frequency generation with phase-sensitive detection. Instantaneous frequency and amplitude of the CO internal stretching mode are tracked with a subpicosecond time resolution that is shorter than the vibrational dephasing time. These experimental results in combination with numerical analysis based on Langevin simulations enable us to extract nonequilibrium distributions of external vibrational modes of desorbing molecules. Superstatistical distributions are generated with mode-dependent frictional couplings in a few hundred femtoseconds after hot-electron excitation, and energy flow from hot electrons and intermode anharmonic coupling play crucial roles in the subsequent evolution of the non-Boltzman distributions.

How Atomic Steps Modify Diffusion and Inter-adsorbate Forces: Empirical Evidence from Hopping Dynamics in Na/Cu(115)

We followed the collective atomic-scale motion of Na atoms on a vicinal Cu(115) surface within a time scale of pico- to nanoseconds using helium spin echo spectroscopy. The well-defined stepped structure of Cu(115) allows us to study the effect that atomic steps have on the adsorption properties, the rate for motion parallel and perpendicular to the step edge, and the interaction between the Na atoms. With the support of a molecular dynamics simulation we show that the Na atoms perform strongly anisotropic 1D hopping motion parallel to the step edges. Furthermore, we observe that the spatial and temporal correlations between the Na atoms that lead to collective motion are also anisotropic, suggesting the steps efficiently screen the lateral interaction between Na atoms residing on different terraces.

Indication of non-thermal contribution to visible femtosecond laser-induced CO oxidation on Ru(0001)

We studied CO oxidation on Ru(0001) induced by 400 nm and 800 nm femtosecond laser pulses where we find a branching ratio between CO oxidation and desorption of 1:9 and 1:31, respectively, showing higher selectivity towards CO oxidation for the shorter wavelength excitation. Activation energies computed with density functional theory show discrepancies with values extracted from the experiments, indicating both a mixture between different adsorbed phases and importance of non-adiabatic effects on the effective barrier for oxidation. We simulated the reactions using kinetic modeling based on the two-temperature model of laser-induced energy transfer in the substrate combined with a friction model for the coupling to adsorbate vibrations. This model gives an overall good agreement with experiment except for the substantial difference in yield ratio between CO oxidation and desorption at 400 nm and 800 nm. However, including also the initial, non-thermal effect of electrons transiently excited into antibonding states of the O-Ru bond yielded good agreement with all experimental results.

Differences between thermal and laser-induced diffusion

A combination of femtosecond laser excitation with a low-temperature scanning tunneling microscope is used to study long-range interaction during diffusion of CO on Cu(111). Both thermal and laser-driven diffusion show an oscillatory energy dependence on the distance to neighboring molecules. Surprisingly, the phase is inverted; i.e., at distances at which thermal diffusion is most difficult, it is easiest for laser-driven diffusion and vice versa. We explain this unexpected behavior by a transient stabilization of the negative ion during diffusion as corroborated by ab initio calculations.

Cross-term selective, two-pulse correlation measurements by phase-shifted parallel modulation for analysis of a multi-photon process

We propose a new modulation technique for phase-sensitive detection to extract a cross-term signal appearing in multi-pulse experiments and demonstrate the selective detection of cross-term signal on two-pulse correlation measurements. This modulation technique can be widely applied to the selective detection in multi-beam laser experiments such as two-pulse correlation, auto-correlation, and double-resonance experiments, because the selectivity can be realized by using a conventional single-frequency mechanical chopper and 2f phase-sensitive detection. The effectiveness of the proposed technique was confirmed by two-pulse correlation measurements of two-photon-excited fluorescence from rhodamine 6G. In addition to the selective observation of a correlation peak, the technique succeeded in observing a weak correlation background. By combining asymmetric, two-pulse correlation measurements with the proposed technique, it was clarified that the background was produced by a three-photon process and assigned to fluorescence depletion caused by an up-conversion process. The results indicate that the cross-term selective, two-pulse correlation method that is based on the phase-shifted parallel modulation technique is powerful tool to find and analyze a high-order optical event buried in a low-order optical event.

The effects of electron-hole pair coupling on the infrared laser-controlled vibrational excitation of NO on Au(111)

In this work, we present theoretical simulations of laser-driven vibrational control of NO adsorbed on a gold surface. Our goal is to tailor laser pulses to selectively excite specific modes and vibrational eigenstates, as well as to favor photodesorption of the adsorbed molecule. To this end, various control schemes and algorithms are applied. For adsorbates at metallic surfaces, the creation of electron-hole pairs in the substrate is known to play a dominant role in the transfer of energy from the system to the surroundings. These nonadiabatic couplings are included perturbatively in our reduced density matrix simulations using a generalization of the state-resolved position-dependent anharmonic rate model we recently introduced. An extension of the reduced density matrix is also proposed to provide a sound model for photodesorption in dissipative systems.

Quantum corrected Langevin dynamics for adsorbates on metal surfaces interacting with hot electrons

We investigate the importance of including quantized initial conditions in Langevin dynamics for adsorbates interacting with a thermal reservoir of electrons. For quadratic potentials the time evolution is exactly described by a classical Langevin equation and it is shown how to rigorously obtain quantum mechanical probabilities from the classical phase space distributions resulting from the dynamics. At short time scales, classical and quasiclassical initial conditions lead to wrong results and only correctly quantized initial conditions give a close agreement with an inherently quantum mechanical master equation approach. With CO on Cu(100) as an example, we demonstrate the effect for a system with ab initio frictional tensor and potential energy surfaces and show that quantizing the initial conditions can have a large impact on both the desorption probability and the distribution of molecular vibrational states.

Diffusion and dimer formation of CO molecules induced by femtosecond laser pulses

We investigate two fundamental steps of a nonadiabatic surface process, the photo-induced movement and approach of CO molecules on the Cu(111) surface, at a hitherto unachieved single-molecule level through scanning tunneling microscope imaging. For the close approach of two CO molecules, we not only determine the nonadiabatic diffusion barrier (87 meV), but also discover a femto-second-laser-induced transient attraction (30 meV) of the usually repelling CO molecules.

Origin of power laws for reactions at metal surfaces mediated by hot electrons

A wide range of experiments have established that certain chemical reactions at metal surfaces can be driven by multiple hot-electron-mediated excitations of adsorbates. A high transient density of hot electrons is obtained by means of femtosecond laser pulses and a characteristic feature of such experiments is the emergence of a power law dependence of the reaction yield on the laser fluence Y approximately F(n). We propose a model of multiple inelastic scattering by hot electrons which reproduces this power law and the observed exponents of several experiments. All parameters are calculated within density functional theory and the delta self-consistent field method. With a simplifying assumption, the power law becomes exact and we obtain a simple physical interpretation of the exponent n, which represents the number of adsorbate vibrational states participating in the reaction.

Temperature-dependent femtosecond photoinduced desorption in CO/Pd(111)

The desorption of CO from a Pd(111) surface following absorption of 120 fs pulses of 780 nm light occurs on two distinct and well-separated time scales. Two-pulse correlation measurements show a fast subpicosecond decay followed by a slower, approximately 40 ps decay. Simulations based on the two-temperature model of electron and phonon heat baths within the substrate, and an empirical friction model to treat coupling to the adsorbate, support the assignment of the desorption mechanism as an electron-mediated process. The photodesorption yield and overall width of the temporal response exhibit a marked dependence on the initial surface temperature in the 100-375 K range despite the much higher transient electronic temperatures (approximately 7000 K) achieved. The observed temperature dependences can be attributed directly to variations in the initial temperature within the frictional coupling picture. Simulations of this extended data set imply that the activation barrier to photoinduced desorption is equal in magnitude to that derived from thermal desorption experiments for this system within the limits of a one-dimensional Arrhenius desorption model. The simulations also imply that the slower decay is not the result of phonon-driven desorption. Though we cannot unambiguously determine the strength of the adsorbate-phonon coupling, our results suggest that its role is to moderate the degree of the adsorbate excitation.

Adsorption-state-dependent subpicosecond photoinduced desorption dynamics

Femtosecond laser excitation has been used to initiate desorption of molecular oxygen from the (111) surface of Pd and to study the adsorption-state dependence of the substrate-adsorbate coupling. The relative populations of the two chemical states, peroxo (O2(2-)) and superoxo (O2-), were varied by changing the total coverage. Two-pulse correlation measurements exhibit a dominant 400 fs response and a slower 10 ps decay that are relatively independent of the initial O2 coverage. In contrast, the photodesorption yield and the nonlinearity of the fluence dependence show a systematic coverage dependence. The coverage-independent subpicosecond response indicates that the photoinduced desorption from the two states is driven primarily by the same electron-mediated mechanism, while the coverage dependence of the yield indicates that the desorption efficiency from the superoxo state is greater than that from the peroxo state. These results are discussed in the context of the electron-phonon two-temperature model with an empirical adsorbate-electron frictional coupling that depends on both the electronic temperature and the activation energy for desorption. With a coupling strength that decreases as the activation energy decreases, the trends with varying coverage, absorbed fluence, and time delay can all be reproduced. The model is consistent with a transition from a resonantly enhanced (diabatic) regime to an adiabatic regime as the system relaxes, accounting for the biexponential correlation behavior.

An instrument to investigate femtochemistry on metal surfaces in real space

A newly established combination of a femtosecond laser with a low temperature scanning tunneling microscope is described, which facilitates one to analyze femtochemistry on metal surfaces in real space. The combined instrument enables focusing the laser to some tens of micrometers and guiding it reproducibly into the tunneling gap with the aid of in situ movable mirrors. Furthermore, a method to determine the focus size on the sample is presented. The focus size is used to calculate the electron and phonon temperatures at the surface. Despite the additional noise introduced by laser operation the vertical resolution of the microscope lies below 1 pm. The potential of the instrument is demonstrated on para-chloronitrobenzene clusters adsorbed on Au(111). Single chloronitrobenzene molecules diffuse upon femtosecond laser irradiation; some smaller clusters rotate by multiples of 30 degrees ; clusters of less compact form rearrange to close-packed clusters.

Femtosecond laser induced associative desorption of H2 from Ru(0001): Comparison of “first principles” theory with experiment

A three dimensional model based on molecular dynamics with electronic frictions is developed to describe the femtosecond laser induced associative desorption of H2 from Ru(0001)(1 x 1)H. Two molecular coordinates (internuclear separation d and center of mass distance to surface z) and a single phonon coordinate are included in the dynamics. Both the potential energy surface and the electronic friction tensor are calculated by density functional theory so that there are no adjustable parameters in the comparison of this model with the wide range of experiments available for this system. This "first principles" dynamic model gives results in semiquantitative agreement with all experimental results; nonlinear fluence dependence of the yield, isotope effect, two pulse correlation, and energy partitioning. The good agreement of theory with experiment supports a description of this surface femtochemistry in terms of thermalized hot electron induced chemistry with coupling to nuclear coordinates through electronic frictions. By comparing the dynamics with the analytical one dimensional frictional model used previously to fit the experiments for this system, we show that the success of the one dimensional model is based on the rapid intermixing of the z and d coordinates as the H-H climbs out of the adsorption well. However, projecting the three dimensional dynamics onto one dimension introduces a fluence (adsorbate temperature) dependent "entropic" barrier in addition to the potential barrier for the chemistry. This implies that some caution must be used in interpreting activation energies obtained in fitting experiments to the one dimensional model.

Real-Time Observation of Molecular Motion on a Surface

The laser-induced movement of CO molecules over a platinum surface was followed in real time by means of ultrafast vibrational spectroscopy. Because the CO molecules bound on different surface sites exhibit different C-O stretch vibrational frequencies, the site-to-site hopping, triggered by excitation with a laser pulse, can be determined from subpicosecond changes in the vibrational spectra. The unexpectedly fast motion--characterized by a 500-femtosecond time constant--reveals that a rotational motion of the CO molecules, rather than pure translation, is required for this diffusion process. This conclusion is corroborated by density functional theory calculations.