The most commonly used method for the time-domain Maxwell equations is the Finite-Difference Time-Domain method (FDTD). This is an explicit, second-order accurate method, which is used on a staggered Cartesian grid. The main drawback with the FDTD method is its inability to accurately model curved objects and small geometrical features. This is due to the Cartesian grid, which leads to a staircase approximation of the geometry and small details are not resolved at all. This thesis presents different ways to circumvent this drawback, but still take advantage of the benefits of the FDTD method. An approach to avoid staircasing errors but still retain the efficiency of the FDTD method is to use a hybrid grid. A few layers of unstructured cells are used close to curved objects and a Cartesian grid is used for the rest of the domain. For the choice of solver on the unstructured grid two different alternatives are compared: an explicit Finite-Volume Time-Domain (FVTD) solver and an implicit Finite-Element Time-Domain (FETD) solver. The hybrid solvers calculate the scattering from complex objects much more efficiently compared to using FDTD on highly resolved Cartesian grids. For the same accuracy in the solution roughly a factor of 10 in memory requirements and a factor of 20 in execution time are gained. The ability to model features that are small relative to the cell size is often important in electromagnetic simulations. In this thesis a technique to generalize a well-known subcell model for thin wires, in order to take arbitrarily oriented wires in FETD and FDTD into account, is proposed. The method gives considerable modeling flexibility compared to earlier methods and is proven stable. The results show excellent consistency and very good accuracy on different antenna configurations. The recursive convolution method is often used to model frequency dispersive materials in FDTD. This method is used to enable modeling of such materials in the unstructured FVTD and FETD solvers. The stability of both solvers is analyzed and their accuracy is demonstrated by computing the radar cross section for homogeneous as well as layered spheres with frequency dependent permittivity.

<p dir="ltr">New classes of electromagnetic modeling problems are arising in the design of quantum information processing technologies. In many cases, these problems involve looking at controlling individual quantum systems with electromagnetic fields at optical frequencies. Since the optical control fields in these applications are essentially classical fields, one promising modeling approach is to perform a semiclassical analysis where the optical fields are treated classically, but the quantum system’s dynamics are modeled quantum mechanically. Typically referred to as Maxwell-Schrödinger models, there is growing interest in these applications to solve the “Maxwell” part of the system directly in terms of the electromagnetic potentials that are used in the Schrödinger equation rather than using the more conventional electric and magnetic fields. To date, the most popular numerical method used to discretize this system of equations is the finite-difference time-domain (FDTD) method. However, these prior works are missing the necessary utilities to consider the modeling of integrated photonic systems when using FDTD methods formulated directly in terms of the magnetic vector and electric scalar potentials. In particular, these potential-based FDTD methods have not been able to model dispersive dielectric materials, which are critical in integrated photonic systems. In this work, we introduce a kind of auxiliary differential equation method for incorporating a Drude-Lorentz-Sommerfeld material model into potential-based FDTD methods. This work also shows the functionality of an absorbing boundary condition for the first time within a potential-based FDTD model to replace the particularly complex implementation of perfectly matched layers within this modeling framework. As such, the methods described in this thesis are intended to help improve the modeling capabilities of potential-based FDTD methods so that they can be used in Maxwell-Schrödinger modeling of more realistic integrated quantum photonic technologies in the future.</p>

L’objectif de cette thèse est de développer une méthode Galerkine discontinue d’ordre élevé capable de prendre en considération des simulations réalistes liées à la nanophotonique. Au cours des dernières décennies, l’évolution des techniques de lithographie a permis la création de structure géométriques de tailles nanométriques, révélant ainsi une large gamme de phénomènes nouveaux nés de l’interaction lumière-matière à ces échelles. Ces effets apparaissent généralement pour des objets de taille égale ou (très) inférieure à la longueur d’onde du champ incident. Ce travail repose sur le développement et l’implémentation de modèles de dispersion appropriés (principalement pour les métaux), ainsi que sur un large éventail de méthodes computationnelles classiques. Deux développements méthodologiques majeurs sont présentés et étudiés en détails: (i) les éléments courbes, et (ii) l’ordre d’approximation local. Ces études sont accompagnées de plusieurs cas-tests réalistes tirés de la nanophotonique. / The goal of this thesis is to develop a discontinuous Galerkin time-domain method to be able to handle realistic nanophotonics computations. During the last decades, the evolution of lithography techniques allowed the creation of geometrical structures at the nanometer scale, thus unveiling a variety of new phenomena arising from light-matter interactions at such levels. These effects usually occur when the device is of comparable size or (much) smaller than the wavelength of the incident field. This work relies on the development and implementation of appropriate models for dispersive materials (mostly metals), as well as on a large panel of classical computational techniques. Two major methodological developments are presented and studied in details: (i) curvilinear elements, and (ii) local order of approximation. This work is complemented with several physical studies of real-life nanophotonics applications.