Skip to Main content Skip to Navigation
Conference papers

3D Simulations of spatially dispersive Metals with a discontinuous Galerkin Time Domain Method

Nikolai Schmitt 1 Stéphane Lanteri 1 Claire Scheid 1
1 NACHOS - Numerical modeling and high performance computing for evolution problems in complex domains and heterogeneous media
CRISAM - Inria Sophia Antipolis - Méditerranée , JAD - Laboratoire Jean Alexandre Dieudonné : UMR6621
Abstract : We present recent advances in the development of a Finite Element Time Domain (Discontinuous Galerkin) solver for computational nanophotonics. Throughout this contribution a particular focus is put on metallic nano structures of sizes between 1 nm and 15 nm. Metal structures at these sizes are well known to show spatial dispersion which can be modeled by a nonlocal dispersion model for the electron gas. Taking such a nonlocal model into account leads to a hydrodynamic fluid equation for the free electrons in the metal. While Maxwell's equations describe the evolution of the electromagnetic fields, the additional fluid equation accounts for the material response and is strongly coupled to Maxwell's equations by means of a source current [1]. Numerically speaking, the considered 3D finite element time domain method benefits from high-order polynomial solutions on very flexible unstructured tetrahedral meshes. Additionally, working in time domain gives access to broad band frequency solution within a single simulation run due to short pulses. Exploiting distributed memory parallelism eventually allows large computational domains and hence realistic nanophotonic setups. We assess the performance of our numerical method on multiple nanophotonic scenarios like spherical dimer systems [2] and nanocubes [3]. We emphasize the importance of highly accurate numerical schemes that guarantee a powerful resolution of surface effects, which is especially indispensable for plasmonic applications since most physics happen in the vicinity of the metal surface. Moreover, we show the importance of high-order numerical methods including curvilinear tetrahedral meshes in order to properly approximate material interfaces if roundings are of concern. References: [1] Schmitt N., Scheid C., Lanteri S., Moreau A., Viquerat J., " A DGTD method for the numerical modeling of the interaction of light with nanometer scale metallic structures taking into account non-local dispersion effects "
Document type :
Conference papers
Complete list of metadatas
Contributor : Nikolai Schmitt <>
Submitted on : Thursday, January 4, 2018 - 3:59:07 PM
Last modification on : Monday, October 12, 2020 - 2:28:04 PM


  • HAL Id : hal-01675094, version 1


Nikolai Schmitt, Stéphane Lanteri, Claire Scheid. 3D Simulations of spatially dispersive Metals with a discontinuous Galerkin Time Domain Method. KWT 2017, Sep 2017, Riezlern, Austria. ⟨hal-01675094⟩



Record views