A hybridizable discontinuous Galerkin method for computing nonlocal electromagnetic effects in three-dimensional metallic nanostructures

F. Vidal-Codina, N. C. Nguyen, S. H. Oh, J. Peraire

Research output: Contribution to journalArticle

  • 1 Citations

Abstract

The interaction of light with metallic nanostructures produces a collective excitation of electrons at the metal surface, also known as surface plasmons. These collective excitations lead to resonances that enable the confinement of light in deep-subwavelength regions, thereby leading to large near-field enhancements. The simulation of plasmon resonances presents notable challenges. From the modeling perspective, the realistic behavior of conduction-band electrons in metallic nanostructures is not captured by Maxwell's equations, thus requiring additional modeling. From the simulation perspective, the disparity in length scales stemming from the extreme field localization demands efficient and accurate numerical methods. In this paper, we develop the hybridizable discontinuous Galerkin (HDG) method to solve Maxwell's equations augmented with the hydrodynamic model for the conduction-band electrons in noble metals. This method enables the efficient simulation of plasmonic nanostructures while accounting for the nonlocal interactions between electrons and the incident light. We introduce a novel postprocessing scheme to recover superconvergent solutions and demonstrate the convergence of the proposed HDG method for the simulation of a 2D gold nanowire and a 3D periodic annular nanogap structure. The results of the hydrodynamic model are compared to those of a simplified local response model, showing that differences between them can be significant at the nanoscale.

LanguageEnglish (US)
Pages548-565
Number of pages18
JournalJournal of Computational Physics
Volume355
DOIs
StatePublished - Feb 15 2018

Fingerprint

Magnetoelectric effects
Galerkin method
Galerkin methods
Nanostructures
electromagnetism
Electrons
Maxwell equations
Conduction bands
Maxwell equation
conduction bands
electrons
Hydrodynamics
simulation
hydrodynamics
Plasmons
Precious metals
noble metals
plasmons
Nanowires
excitation

Keywords

  • Hybridizable discontinuous Galerkin method
  • Hydrodynamic model for metals
  • Maxwell's equations
  • Nonlocal electrodynamics
  • Plasmonics
  • Terahertz nonlocality

Cite this

A hybridizable discontinuous Galerkin method for computing nonlocal electromagnetic effects in three-dimensional metallic nanostructures. / Vidal-Codina, F.; Nguyen, N. C.; Oh, S. H.; Peraire, J.

In: Journal of Computational Physics, Vol. 355, 15.02.2018, p. 548-565.

Research output: Contribution to journalArticle

@article{4b053452909a4dba8606b4c6224e1926,
title = "A hybridizable discontinuous Galerkin method for computing nonlocal electromagnetic effects in three-dimensional metallic nanostructures",
abstract = "The interaction of light with metallic nanostructures produces a collective excitation of electrons at the metal surface, also known as surface plasmons. These collective excitations lead to resonances that enable the confinement of light in deep-subwavelength regions, thereby leading to large near-field enhancements. The simulation of plasmon resonances presents notable challenges. From the modeling perspective, the realistic behavior of conduction-band electrons in metallic nanostructures is not captured by Maxwell's equations, thus requiring additional modeling. From the simulation perspective, the disparity in length scales stemming from the extreme field localization demands efficient and accurate numerical methods. In this paper, we develop the hybridizable discontinuous Galerkin (HDG) method to solve Maxwell's equations augmented with the hydrodynamic model for the conduction-band electrons in noble metals. This method enables the efficient simulation of plasmonic nanostructures while accounting for the nonlocal interactions between electrons and the incident light. We introduce a novel postprocessing scheme to recover superconvergent solutions and demonstrate the convergence of the proposed HDG method for the simulation of a 2D gold nanowire and a 3D periodic annular nanogap structure. The results of the hydrodynamic model are compared to those of a simplified local response model, showing that differences between them can be significant at the nanoscale.",
keywords = "Hybridizable discontinuous Galerkin method, Hydrodynamic model for metals, Maxwell's equations, Nonlocal electrodynamics, Plasmonics, Terahertz nonlocality",
author = "F. Vidal-Codina and Nguyen, {N. C.} and Oh, {S. H.} and J. Peraire",
year = "2018",
month = "2",
day = "15",
doi = "10.1016/j.jcp.2017.11.025",
language = "English (US)",
volume = "355",
pages = "548--565",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",

}

TY - JOUR

T1 - A hybridizable discontinuous Galerkin method for computing nonlocal electromagnetic effects in three-dimensional metallic nanostructures

AU - Vidal-Codina,F.

AU - Nguyen,N. C.

AU - Oh,S. H.

AU - Peraire,J.

PY - 2018/2/15

Y1 - 2018/2/15

N2 - The interaction of light with metallic nanostructures produces a collective excitation of electrons at the metal surface, also known as surface plasmons. These collective excitations lead to resonances that enable the confinement of light in deep-subwavelength regions, thereby leading to large near-field enhancements. The simulation of plasmon resonances presents notable challenges. From the modeling perspective, the realistic behavior of conduction-band electrons in metallic nanostructures is not captured by Maxwell's equations, thus requiring additional modeling. From the simulation perspective, the disparity in length scales stemming from the extreme field localization demands efficient and accurate numerical methods. In this paper, we develop the hybridizable discontinuous Galerkin (HDG) method to solve Maxwell's equations augmented with the hydrodynamic model for the conduction-band electrons in noble metals. This method enables the efficient simulation of plasmonic nanostructures while accounting for the nonlocal interactions between electrons and the incident light. We introduce a novel postprocessing scheme to recover superconvergent solutions and demonstrate the convergence of the proposed HDG method for the simulation of a 2D gold nanowire and a 3D periodic annular nanogap structure. The results of the hydrodynamic model are compared to those of a simplified local response model, showing that differences between them can be significant at the nanoscale.

AB - The interaction of light with metallic nanostructures produces a collective excitation of electrons at the metal surface, also known as surface plasmons. These collective excitations lead to resonances that enable the confinement of light in deep-subwavelength regions, thereby leading to large near-field enhancements. The simulation of plasmon resonances presents notable challenges. From the modeling perspective, the realistic behavior of conduction-band electrons in metallic nanostructures is not captured by Maxwell's equations, thus requiring additional modeling. From the simulation perspective, the disparity in length scales stemming from the extreme field localization demands efficient and accurate numerical methods. In this paper, we develop the hybridizable discontinuous Galerkin (HDG) method to solve Maxwell's equations augmented with the hydrodynamic model for the conduction-band electrons in noble metals. This method enables the efficient simulation of plasmonic nanostructures while accounting for the nonlocal interactions between electrons and the incident light. We introduce a novel postprocessing scheme to recover superconvergent solutions and demonstrate the convergence of the proposed HDG method for the simulation of a 2D gold nanowire and a 3D periodic annular nanogap structure. The results of the hydrodynamic model are compared to those of a simplified local response model, showing that differences between them can be significant at the nanoscale.

KW - Hybridizable discontinuous Galerkin method

KW - Hydrodynamic model for metals

KW - Maxwell's equations

KW - Nonlocal electrodynamics

KW - Plasmonics

KW - Terahertz nonlocality

UR - http://www.scopus.com/inward/record.url?scp=85036573731&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85036573731&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2017.11.025

DO - 10.1016/j.jcp.2017.11.025

M3 - Article

VL - 355

SP - 548

EP - 565

JO - Journal of Computational Physics

T2 - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

ER -