Applications of the wave packet method to resonant transmission and reflection gratings

Andrei G. Borisov,Sergei V. Shabanov

DOI: https://doi.org/10.1016/j.jcp.2004.03.013

2003-12-17

Abstract:Scattering of femtosecond laser pulses on resonant transmission and reflection gratings made of dispersive (Drude metals) and dielectric materials is studied by a time-domain numerical algorithm for Maxwell's theory of linear passive (dispersive and absorbing) media. The algorithm is based on the Hamiltonian formalism in the framework of which Maxwell's equations for passive media are shown to be equivalent to the first-order equation, $\partial \Psi/\partial t = {\cal H}\Psi$, where ${\cal H}$ is a linear differential operator (Hamiltonian) acting on a multi-dimensional vector $\Psi$ built of the electromagnetic inductions and auxiliary matter fields describing the medium response. The initial value problem is then solved by means of a modified time leapfrog method in combination with the Fourier pseudospectral method applied on a non-uniform grid that is constructed by a change of variables and designed to enhance the sampling efficiency near medium interfaces. The algorithm is shown to be highly accurate at relatively low computational costs. An excellent agreement with previous theoretical and experimental studies of the gratings is demonstrated by numerical simulations using our algorithm. In addition, our algorithm allows one to see real time dynamics of long leaving resonant excitations of electromagnetic fields in the gratings in the entire frequency range of the initial wide band wave packet as well as formation of the reflected and transmitted wave fronts.

Computational Physics,Optics

What problem does this paper attempt to address?

This paper attempts to solve two main problems:

1. **Testing a new time - domain algorithm**: First, the author aims to test a new time - domain numerical algorithm for linear passive (dispersive and absorptive) media based on Maxwell's theory. This algorithm combines the following three key elements:
- Hamiltonian formalism, which is used for evolving differential equations.
- The modified time leapfrog method, which is used for time evolution.
- The Fourier pseudospectral method combined with variable transformation, in order to improve the spatial grid resolution in specific areas (such as near the medium interface).

2. **Studying the scattering phenomenon of electromagnetic waves on gratings**: Second, the author applies this algorithm to study the scattering phenomenon of femtosecond laser pulses on resonant transmission and reflection gratings made of dispersive (Drude metal) and dielectric materials. Specifically, they focus on the scattering behavior of broadband electromagnetic pulses on these gratings, including the formation and decay of long - lived resonant excitation, as well as the formation of reflected and transmitted wavefronts.

Through these studies, the author hopes to verify that this algorithm has high precision at a relatively low computational cost and can observe the real - time dynamic changes of broadband electromagnetic pulses when they propagate in passive media. In addition, this algorithm can also help to understand the mechanisms of certain abnormal transmission phenomena (for example, the high transmittance of periodic thin - film metal gratings at specific wavelengths).

### Formula summary

- **Maxwell's equations in Hamiltonian formalism**:
\[
\frac{\partial \Psi}{\partial t}=H \Psi
\]
where \(\Psi\) is a multi - dimensional vector, containing electromagnetic induction and auxiliary fields that describe the medium response, and \(H\) is a linear differential operator (Hamiltonian).

- **Initial value problem**:
\[
\Psi(0)=\begin{pmatrix} D \\ B \end{pmatrix}
\]

- **Boundary conditions**:
\[
\nabla \cdot B(t)=\nabla \cdot D(t) = 0
\]

- **Energy conservation**:
\[
E(t)=\int dr\left(\frac{1}{2}\left[E^{2}+B^{2}+\sum_{a}(\pi_{a}^{2}+\omega_{a}^{2}\vartheta_{a}^{2})\right]\right)
\]

- **Fourier transform and non - uniform grid**:
\[
\tilde{\Psi}(k)=\sum_{n}F_{nn'}\Psi(n\Delta r)
\]

Through these formulas and methods, the author has successfully developed an efficient and accurate numerical algorithm for simulating the propagation and scattering behavior of broadband electromagnetic pulses in complex media.