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

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.Comment: 23 pages; 8 figures in the png forma

Doing numerical cosmology with the Cactus code

The article presents some aspects concerning the construction of a new thorn for the Cactus code, a complete 3-dimensional machinery for numerical relativity. This thorn is completely dedicated to numerical simulations in cosmology, that means it can provide evolutions of different cosmological models, mainly based on Friedman-Robertson-Walker metric. Some numerical results are presented, testing the convergence, stability and the applicability of the code.Comment: 18 pages, 8 figures, Late