skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A DGTD method for the numerical modeling of the interaction of light with nanometer scale metallic structures taking into account non-local dispersion effects

Abstract

The interaction of light with metallic nanostructures is increasingly attracting interest because of numerous potential applications. Sub-wavelength metallic structures, when illuminated with a frequency close to the plasma frequency of the metal, present resonances that cause extreme local field enhancements. Exploiting the latter in applications of interest requires a detailed knowledge about the occurring fields which can actually not be obtained analytically. For the latter mentioned reason, numerical tools are thus an absolute necessity. The insight they provide is very often the only way to get a deep enough understanding of the very rich physics at play. For the numerical modeling of light-structure interaction on the nanoscale, the choice of an appropriate material model is a crucial point. Approaches that are adopted in a first instance are based on local (i.e. with no interaction between electrons) dispersive models, e.g. Drude or Drude–Lorentz models. From the mathematical point of view, when a time-domain modeling is considered, these models lead to an additional system of ordinary differential equations coupled to Maxwell's equations. However, recent experiments have shown that the repulsive interaction between electrons inside the metal makes the response of metals intrinsically non-local and that this effect cannot generally be overlooked. Technologicalmore » achievements have enabled the consideration of metallic structures in a regime where such non-localities have a significant influence on the structures' optical response. This leads to an additional, in general non-linear, system of partial differential equations which is, when coupled to Maxwell's equations, significantly more difficult to treat. Nevertheless, dealing with a linearized non-local dispersion model already opens the route to numerous practical applications of plasmonics. In this work, we present a Discontinuous Galerkin Time-Domain (DGTD) method able to solve the system of Maxwell's equations coupled to a linearized non-local dispersion model relevant to plasmonics. While the method is presented in the general 3D case, numerical results are given for 2D simulation settings.« less

Authors:
 [1];  [2];  [1];  [3];  [1];  [4];  [1]
  1. Inria, 2004 Route des Lucioles, BP 93, 06902 Sophia Antipolis Cedex (France)
  2. (TEMF), Schlossgartenstr. 8, 64289 Darmstadt (Germany)
  3. (France)
  4. Institut Pascal, Université Blaise Pascal, 24, avenue des Landais, 63171 Aubière Cedex (France)
Publication Date:
OSTI Identifier:
22572323
Resource Type:
Journal Article
Resource Relation:
Journal Name: Journal of Computational Physics; Journal Volume: 316; Other Information: Copyright (c) 2016 Elsevier Science B.V., Amsterdam, The Netherlands, All rights reserved.; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; INTERACTIONS; LANGMUIR FREQUENCY; NANOSTRUCTURES; NONLINEAR PROBLEMS; PARTIAL DIFFERENTIAL EQUATIONS; PLASMA; SIMULATION; WAVELENGTHS

Citation Formats

Schmitt, Nikolai, Technische Universitaet Darmstadt, Institut fuer Theorie Elektromagnetischer Felder, Scheid, Claire, University of Nice – Sophia Antipolis, Mathematics laboratory, Parc Valrose, 06108 Nice, Cedex 02, Lanteri, Stéphane, E-mail: Stephane.Lanteri@inria.fr, Moreau, Antoine, and Viquerat, Jonathan. A DGTD method for the numerical modeling of the interaction of light with nanometer scale metallic structures taking into account non-local dispersion effects. United States: N. p., 2016. Web. doi:10.1016/J.JCP.2016.04.020.
Schmitt, Nikolai, Technische Universitaet Darmstadt, Institut fuer Theorie Elektromagnetischer Felder, Scheid, Claire, University of Nice – Sophia Antipolis, Mathematics laboratory, Parc Valrose, 06108 Nice, Cedex 02, Lanteri, Stéphane, E-mail: Stephane.Lanteri@inria.fr, Moreau, Antoine, & Viquerat, Jonathan. A DGTD method for the numerical modeling of the interaction of light with nanometer scale metallic structures taking into account non-local dispersion effects. United States. doi:10.1016/J.JCP.2016.04.020.
Schmitt, Nikolai, Technische Universitaet Darmstadt, Institut fuer Theorie Elektromagnetischer Felder, Scheid, Claire, University of Nice – Sophia Antipolis, Mathematics laboratory, Parc Valrose, 06108 Nice, Cedex 02, Lanteri, Stéphane, E-mail: Stephane.Lanteri@inria.fr, Moreau, Antoine, and Viquerat, Jonathan. Fri . "A DGTD method for the numerical modeling of the interaction of light with nanometer scale metallic structures taking into account non-local dispersion effects". United States. doi:10.1016/J.JCP.2016.04.020.
@article{osti_22572323,
title = {A DGTD method for the numerical modeling of the interaction of light with nanometer scale metallic structures taking into account non-local dispersion effects},
author = {Schmitt, Nikolai and Technische Universitaet Darmstadt, Institut fuer Theorie Elektromagnetischer Felder and Scheid, Claire and University of Nice – Sophia Antipolis, Mathematics laboratory, Parc Valrose, 06108 Nice, Cedex 02 and Lanteri, Stéphane, E-mail: Stephane.Lanteri@inria.fr and Moreau, Antoine and Viquerat, Jonathan},
abstractNote = {The interaction of light with metallic nanostructures is increasingly attracting interest because of numerous potential applications. Sub-wavelength metallic structures, when illuminated with a frequency close to the plasma frequency of the metal, present resonances that cause extreme local field enhancements. Exploiting the latter in applications of interest requires a detailed knowledge about the occurring fields which can actually not be obtained analytically. For the latter mentioned reason, numerical tools are thus an absolute necessity. The insight they provide is very often the only way to get a deep enough understanding of the very rich physics at play. For the numerical modeling of light-structure interaction on the nanoscale, the choice of an appropriate material model is a crucial point. Approaches that are adopted in a first instance are based on local (i.e. with no interaction between electrons) dispersive models, e.g. Drude or Drude–Lorentz models. From the mathematical point of view, when a time-domain modeling is considered, these models lead to an additional system of ordinary differential equations coupled to Maxwell's equations. However, recent experiments have shown that the repulsive interaction between electrons inside the metal makes the response of metals intrinsically non-local and that this effect cannot generally be overlooked. Technological achievements have enabled the consideration of metallic structures in a regime where such non-localities have a significant influence on the structures' optical response. This leads to an additional, in general non-linear, system of partial differential equations which is, when coupled to Maxwell's equations, significantly more difficult to treat. Nevertheless, dealing with a linearized non-local dispersion model already opens the route to numerous practical applications of plasmonics. In this work, we present a Discontinuous Galerkin Time-Domain (DGTD) method able to solve the system of Maxwell's equations coupled to a linearized non-local dispersion model relevant to plasmonics. While the method is presented in the general 3D case, numerical results are given for 2D simulation settings.},
doi = {10.1016/J.JCP.2016.04.020},
journal = {Journal of Computational Physics},
number = ,
volume = 316,
place = {United States},
year = {Fri Jul 01 00:00:00 EDT 2016},
month = {Fri Jul 01 00:00:00 EDT 2016}
}
  • Analysis of three-dimensional problems associated with high-speed interaction of deformed solids with a solid barrier over a broad range of angles of impact constitutes a task of the highest priority. However, none of the studies of asymmetric impact, either in the two- or the three-dimensional formulation, that have been conducted up until now have focused on the subject of thermal effects. In addition, from ana analysis of studies that have been devoted to axisymmetric problems of high-speed impact, it is clear that the dependence of the characteristics of the material on temperature must be taken into account in order tomore » arrive at a correct description of the processes of strain and destruction of bodies subjected to asymmetric loading. In the present article numerical simulation fo high-speed interaction of cylindrical strikers with a rigid barrire is conducted over a broad range of angles of impact, taking into account the influence of the temperature patterns that are realized in the cylinders in the course of impact.« less
  • In the first part of this paper a discrete multi-dimensional physical-topological model of functionally integrated VLSIC elements which took heat effects into account was constructed and analyzed. In this part, the method of implementation for the model and a universal program for its implementation are described. The method is based on the solution of a truncated system of equations derived from the starting system with the help of a number of assumptions. The key assumption is that the temperature is constant over the VLSIC structure. The equations are used to represent current density, the majority charge carriers, the Fermi quasilevelsmore » of the nonequilibrium charge carriers. Methods of vector relaxation of systems are classified. The program performs two-dimensional simulation of semiconductor devices and is called UNTEMP. It is written in FORTRAN-IV for the ES series of computers.« less
  • We study the possible reasons for the disagreement between the estimates of the pion-nucleon sigma term obtained by the method of dispersion relations with extrapolation to the Cheng-Dashen point and by other methods which do not involve this extrapolation. One reason for the disagreement may be the nonanalyticity of the ..pi..N amplitude in the variable t for ..nu.. = 0. We propose a method for estimating the sigma term using the threshold data for the ..pi..N amplitude, in which the effect of this nonanalyticity is minimized. We discuss the relation between scale invariance violation and chiral symmetry breaking and givemore » the corresponding estimate of the sigma term. The two estimates are similar (42 and 34 MeV) and are in agreement when the uncertainties of the two methods are taken into consideration.« less
  • A discrete physical-topological two-dimensional model intended for calculating the elements of very large scale integrated circuits is constructed which takes into account the thermal effects which impose a fundamental limit on the increase in the degree to which such circuits can be integrated. The model is based on fundamental systems of equations for the physics of semiconductors supplemented by a heat conduction equation. The recombination-generation processes are described using the traditional Shockley-Read-Hall model. Normalizing coefficients for the temperature potentials, coordinates, charge density, diffusion, mobility, recombination-generation rate, lifetime, and thermal conductivity are calculated. A complete difference scheme is obtained which consistsmore » of the difference schemes for the Poisson equation, the equations of continuity for electrons and holes, and the analog of the heat-conduction equation.« less