2025-11-25T04:10:17.947558

Plasmon-assisted photoelectron emission in a model cluster using time-dependent density functional theory and the time-dependent surface-flux method

Bednov, Pervez, Barke et al.
We investigate plasmon-assisted photoelectron emission using a one-dimensional time-dependent density-functional theory (TDDFT) model. Photoelectron spectra are computed with the time-dependent surface-flux (t-SURFF) method. In addition to the expected above-threshold ionization (ATI) comb, we observe peaks that arise from long-lived plasmon oscillations and the associated electron emission occurring after the laser pulse. We further analyze the positions of these peaks and their scaling behavior with the laser intensity.
academic

Plasmon-assisted photoelectron emission in a model cluster using time-dependent density functional theory and the time-dependent surface-flux method

Basic Information

  • Paper ID: 2507.22709
  • Title: Plasmon-assisted photoelectron emission in a model cluster using time-dependent density functional theory and the time-dependent surface-flux method
  • Authors: Mikhail Bednov, Waqas Pervez, Ingo Barke, Dieter Bauer (University of Rostock, Germany)
  • Classification: quant-ph physics.atm-clus physics.comp-ph
  • Publication Date: October 14, 2025 (arXiv preprint)
  • Paper Link: https://arxiv.org/abs/2507.22709

Abstract

This study investigates plasmon-assisted photoelectron emission using a one-dimensional time-dependent density functional theory (TDDFT) model. Photoelectron spectra are calculated using the time-dependent surface-flux (t-SURFF) method. Beyond the expected above-threshold ionization (ATI) comb peaks, the research observes peaks generated by long-lived plasmon oscillations arising from correlated electron emission occurring after the laser pulse. The positions of these peaks and their scaling behavior with laser intensity are further analyzed.

Research Background and Motivation

Problem Background

Plasmonic nanostructures concentrate optical energy into deep subwavelength volumes, enabling strong-field photoemission, ultrafast electron sources, and controlled light-matter interactions at the nanoscale. Ultrafast photoemission electron microscopy (PEEM) has directly imaged the spatiotemporal evolution of plasmon dynamics, providing unique insights into how local fields drive electron emission at the nanoscale.

Theoretical Challenges

From a theoretical perspective, while real-time TDDFT has become a powerful framework for tracking collective electron motion, charge transfer, and field-driven dynamics, the exact observable for photoelectron emission in TDDFT remains unknown. The traditional interpretation of Kohn-Sham (KS) particles themselves as electrons for calculating photoelectron spectra is merely an approximation.

Research Motivation

This study aims to:

  1. Understand the physical mechanisms of plasmon-assisted photoelectron emission
  2. Verify the applicability and limitations of TDDFT in describing this phenomenon
  3. Analyze the influence of plasmon oscillations on photoelectron spectra
  4. Investigate the intensity scaling behavior of plasmon-assisted emission

Core Contributions

  1. Established a 1D TDDFT model: Constructed a one-dimensional cluster model containing 40 electrons capable of effectively simulating plasmon excitation and photoelectron emission
  2. Identified collective modes: Through linear response analysis, identified two main collective oscillation modes (ωA = 0.106, ωB = 0.156)
  3. Discovered plasmon-assisted emission peaks: Observed sharp peak structures in the photoelectron spectrum generated by post-pulse plasmon oscillations
  4. Analyzed intensity scaling behavior: Revealed that plasmon-assisted emission peaks follow specific intensity scaling laws (I^n scaling)
  5. Provided time-of-flight analysis: Distinguished temporal characteristics of direct ATI emission and plasmon-assisted emission through windowed t-SURFF method

Methodology Details

Task Definition

Investigate plasmon excitation and associated photoelectron emission in a one-dimensional electron cluster under laser fields, with particular focus on electron emission from three-photon nonlinear plasmon excitation (3ωL = ωplasmon).

Model Architecture

1. One-Dimensional Cluster Model

A one-dimensional cluster with 40 electrons with ionic potential:

V_ion(x) = -∑_{i=1}^{N_ion} 1/√((x-X_i)² + 1)

where X_i = (-(N_ion+1)/2 + i)a, a = 1.125, chosen such that the KS valence orbital energy approaches the work function of silver clusters.

2. Kohn-Sham Equations

Ground state configuration solved via KS equations:

ε_i φ_i(x) = [-1/2 ∂²/∂x² + V_KS[n](x)] φ_i(x)

The KS potential includes:

  • Ionic potential V_ion(x)
  • Hartree potential V_Hn
  • Exchange-correlation potential V_xcn ≃ V_xn = -(3n(x)/π)^(1/3)

3. Time-Dependent KS Equations

Dynamics in the laser field described by TDKS equations:

i∂_t φ_i(x,t) = [1/2(-i∂_x + A(t))² + V_KS[n](x,t)] φ_i(x,t)

4. Laser Pulse

Laser pulse with sine-squared envelope:

A(t) = a₀ cos(ωLt) sin²(ωLt/(2N_cyc)), 0 < t < Tp

Technical Innovations

1. t-SURFF Photoelectron Spectrum Calculation

In the velocity gauge, momentum space amplitudes at left and right surfaces:

b_{R,i}(k) = 1/√(2π) ∫₀^τ dt e^{ik²t/2} e^{-ik(x_R-α(t))} [(k/2+A(t))φ_i(x_R,t) - i/2 ∂_x φ_i(x_R,t)]

2. Collective Mode Identification

Collective modes identified through linear response analysis following δ-kick excitation:

P(Ω) = Ω⁴|FFT[D(t)](Ω)|², D(t) = ∫dx x n(x,t)

3. Time-of-Flight Analysis

Windowed t-SURFF method with Gaussian time windows distinguishes electron emission contributions from different times.

Experimental Setup

Numerical Parameters

  • Grid points: N = 2000
  • Spatial step: Δx = 0.5
  • Time step: Δt = 0.25
  • t-SURFF surface positions: x_L = 0.25NΔx, x_R = 0.75NΔx
  • Boundary conditions: Complex absorbing boundary

Laser Parameters

  • Amplitude: a₀ = 0.002-0.006
  • Frequency: ωL = 0.035, 0.052 (corresponding to ωA/3, ωB/3)
  • Number of cycles: N_cyc = 20

Numerical Methods

  • Real-time propagation: Crank-Nicolson scheme with predictor-corrector steps
  • Ground state solution: Imaginary time propagation
  • Hartree potential calculation: Convolution theorem

Experimental Results

Main Results

1. Collective Mode Identification

Two main collective modes identified through linear response analysis:

  • ωA = 0.106: Collective oscillation of entire electron density
  • ωB = 0.156: Higher-order oscillation mode with two nodes

2. Photoelectron Spectrum Characteristics

In the photoelectron spectrum at laser frequency ωL = ωB/3 = 0.052:

  • Conventional ATI peaks: spacing ωL
  • Plasmon-assisted emission peaks: located at E_^{(B)} = εi + nωB

3. Time-of-Flight Analysis

Windowed t-SURFF analysis reveals:

  • ATI emission occurs during laser pulse (ωLt/2π ≲ 30)
  • Plasmon-assisted emission continues after pulse, manifesting as horizontal stripes

4. Intensity Scaling Behavior

For ωL = 0.052 = ωB/3:

  • ATI peak (E = 0.037): I^3.88 scaling (theoretical expectation I^4)
  • Plasmon peak (E = 0.071): I^5.71 scaling

For ωL = 0.035 = ωA/3:

  • ATI peak: I^5.02 scaling (theoretical expectation I^5)
  • Plasmon peak: I^6.24 scaling (theoretical expectation I^6)

Physical Interpretation

The scaling behavior of plasmon-assisted emission can be understood as:

  • Plasmon excitation requires 3 laser photons: I^3 scaling
  • Electron escape by absorbing plasmon quantum: additional I^3 scaling (since 1 plasmon quantum = 3 laser photons)
  • Total scaling: I^6, consistent with numerical results

Related Work

Plasmon-Assisted Photoemission Research

  • Surface plasmon resonance-enhanced multiphoton emission
  • Light-field-controlled photoemission from plasmonic nanoparticles
  • Electron emission in strong-field nanophotonics

TDDFT Applications in Plasmon Research

  • Plasmon excitation in carbon nanotubes under linearly polarized light
  • Interaction between plasmons and single-particle excitations in metal nanoclusters
  • Real-time TDDFT simulations of plasmonic nanostructures

Photoelectron Spectrum Calculation Methods

  • Projection onto field-free continuum states method
  • Window operator method
  • Time-dependent surface-flux (t-SURFF) method

Conclusions and Discussion

Main Conclusions

  1. TDDFT can describe plasmon-assisted photoemission: Despite approximations, TDDFT successfully captures the main physical features of plasmon excitation and associated electron emission
  2. Plasmon oscillations produce characteristic photoelectron peaks: Long-lived post-pulse plasmon oscillations lead to sharp characteristic peaks in the photoelectron spectrum
  3. Intensity scaling validates the physical picture: The I^6 scaling of plasmon-assisted emission supports the physical picture of "3-photon plasmon excitation + plasmon quantum-assisted electron emission"

Limitations

1. Limitations of TDDFT Approximations

  • Use of simple exchange-correlation functional (exchange term only) may overestimate plasmon-assisted emission
  • Interpretation of KS orbitals as electron wavefunctions is only approximate
  • Plasmon oscillation damping is too slow, potentially inconsistent with reality

2. Restrictions of 1D Model

  • Landau damping is suppressed in 1D; larger density of states in actual 3D systems would enhance dephasing
  • Lacks real dissipation channels such as coupling to phonons or environment
  • Model is oversimplified and may not fully reflect the complexity of actual nanostructures

3. Quantitative Accuracy Issues

  • While the scaling behavior of plasmon peaks is qualitatively correct, quantitative accuracy is uncertain
  • Similar to the second plateau in high-order harmonic generation, TDDFT may give qualitatively correct but quantitatively inaccurate results

Future Directions

  1. Develop improved exchange-correlation functionals: Plasmon-assisted electron emission serves as a stringent test case
  2. Higher-dimensional model studies: Extend to 2D or 3D models to include more realistic physical effects
  3. Many-body method benchmarking: Systematically benchmark TDDFT using high-precision many-body techniques
  4. Experimental validation: Compare with ultrafast photoemission electron microscopy experimental results

In-Depth Evaluation

Strengths

1. Methodological Innovation

  • First systematic study of plasmon-assisted photoelectron emission using TDDFT
  • Application of t-SURFF method avoids traditional method's assumption of static potential after pulse
  • Windowed t-SURFF provides time-resolved analysis of electron emission

2. Deep Physical Insights

  • Clear distinction between collective modes and single-particle transitions
  • Verification of plasmon-assisted emission physical picture through intensity scaling analysis
  • Spatiotemporal analysis reveals spatial distribution characteristics of plasmon oscillations

3. Comprehensive Analysis

  • Complete analysis chain from linear response to nonlinear photoelectron emission
  • Systematic study across multiple laser parameters and intensities
  • Honest discussion of TDDFT limitations

Shortcomings

1. High Degree of Model Simplification

  • 1D model, while computationally efficient, has limited physical realism
  • Neglects many important physical effects (spin-orbit coupling, many-body correlations, etc.)

2. Questionable Quantitative Reliability

  • Acknowledges that TDDFT may overestimate plasmon-assisted emission, but lacks quantitative assessment
  • No direct comparison with experimental data or high-precision theoretical methods

3. Limited Applicability Range

  • Primarily applicable to simple metal cluster systems
  • Applicability to complex nanostructures requires further verification

Impact

1. Theoretical Method Development

  • Provides new computational framework for theoretical research in plasmon physics
  • Promotes development of TDDFT applications in strong-field physics

2. Experimental Guidance Value

  • Provides theoretical predictions and explanation framework for ultrafast plasmon experiments
  • Aids in designing new plasmon-based electron sources

3. Methodological Contributions

  • Application of t-SURFF in TDDFT has methodological value
  • Time-of-flight analysis method generalizable to other strong-field phenomena studies

Application Scenarios

  1. Plasmonic nanostructure design: Provides theoretical guidance for optimizing plasmon-enhanced photoelectron emission
  2. Ultrafast electron source development: Provides theoretical foundation for designing plasmon-based ultrafast electron sources
  3. Strong-field physics research: Serves as theoretical tool for studying complex light-matter interactions
  4. TDDFT method development: Serves as benchmark problem for testing and improving exchange-correlation functionals

References

The paper cites 25 important references, covering:

  • Experimental studies of plasmon-assisted photoemission 1-9
  • TDDFT theory and applications 10-16
  • Photoelectron spectrum calculation methods 17-22
  • TDDFT limitations and improvements 23-25

These references provide solid theoretical foundation and experimental background, demonstrating the authors' in-depth understanding of the relevant field.