Electron beams traversing spherical nanoparticles: analytic and numerical treatment

P. Elli Stamatopoulou,Wenhua Zhao,Álvaro Rodríguez Echarri,N. Asger Mortensen,Kurt Busch,Christos Tserkezis,Christian Wolff

DOI: https://doi.org/10.1103/PhysRevResearch.6.013239

2024-04-14

Abstract:We present an analytic, Mie theory-based solution for the energy-loss and the photon-emission probabilities in the interaction of spherical nanoparticles with electrons passing nearby and through them, in both cathodoluminescence and electron energy-loss spectroscopies. In particular, we focus on the case of penetrating electron trajectories, for which the complete fully electrodynamic and relativistic formalism has not been reported as yet. We exhibit the efficiency of this method in describing collective excitations in matter through calculations for a dispersive and lossy system, namely a sphere described by a Drude permittivity. Subsequently, we use the analytic solution to corroborate the implementation of electron-beam sources in a state-of-the-art numerical method for problems in electrodynamics, the discontinuous Galerkin time-domain (DGTD) method. We show that the two approaches produce spectra in good mutual agreement, and demonstrate the versatility of DGTD via simulations of spherical nanoparticles characterized by surface roughness. The possibility of simultaneously employing both kinds of calculations (analytic and numerical) facilitates a better understanding of the rich optical response of nanophotonic architectures excited by fast electron beams.

Mesoscale and Nanoscale Physics,Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to accurately describe the energy loss and photon emission probability when the electron beam interacts with spherical nanoparticles. Specifically, the research focuses on the physical processes that occur in cathodoluminescence (CL) and electron energy - loss spectroscopy (EELS) when electrons penetrate or are close to spherical nanoparticles. The paper proposes an analytical solution based on Mie theory and compares it with the numerical method of the Discontinuous Galerkin Time - Domain (DGTD) to verify the effectiveness and consistency of these two methods in describing the optical response of nanoparticles.

### Main problems
1. **Probability of energy loss and photon emission**: The research aims to calculate and analyze in detail the probability of energy loss and photon emission when the electron beam penetrates or approaches spherical nanoparticles through analytical and numerical methods.
2. **Comparison between analytical and numerical methods**: By comparing the analytical method (based on Mie theory) and the numerical method (DGTD), verify the consistency between the two in describing the optical response of nanoparticles.
3. **Influence of surface roughness**: Explore the influence of surface roughness on the optical response of nanoparticles, especially its performance in EELS.

### Specific research contents
- **Analytical method**: Based on Mie theory, the analytical formulas for the energy loss and photon emission probability when the electron penetrates the trajectory are derived. These formulas take into account the relativistic effect and the complete electrodynamic effect.
- **Numerical method**: Use the DGTD method to simulate the electromagnetic response of the same system and compare it with the analytical results.
- **Simulation of surface roughness**: Apply the numerical method to study the optical response of nanoparticles with surface roughness, demonstrating the flexibility of the DGTD method in simulating complex geometric structures and nonlinear responses.

### Main findings
- **Consistency between analytical and numerical methods**: For far - away and penetrating trajectories, the results of the analytical method and the DGTD method show good consistency, especially in the CL spectrum. However, under the grazing trajectory, due to the finite width of the electron beam, there is a certain deviation between the two methods.
- **Accumulation of high - order multipoles**: Near the SPP energy, the accumulation of high - order multipoles is difficult to be completely resolved in the analytical and numerical methods, resulting in a large error.
- **Influence of surface roughness**: The numerical method can effectively simulate nanoparticles with surface roughness, which is more realistic under experimental conditions.

Through these studies, the paper not only provides an in - depth understanding of the optical response of nanoparticles but also provides strong theoretical support for the interpretation of experimental data.