skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Fast and accurate electromagnetic field calculation for substrate-supported metasurfaces using the discrete dipole approximation
Metasurface design tends to be tedious and time-consuming based on sweeping geometric parameters. Common numerical simulation techniques are slow for large areas, ultra-fine grids, and/or three-dimensional simulations. Simulation time can be reduced by combining the principle of the discrete dipole approximation (DDA) with analytical solutions for light scattered by a dipole near a flat surface. The DDA has rarely been used in metasurface design, and comprehensive benchmarking comparisons are lacking. Here, we compare the accuracy and speed of three DDA methods—substrate discretization, two-dimensional Cartesian Green’s functions, and one-dimensional (1D) cylindrical Green’s functions—against the finite difference time domain (FDTD) method. We find that the 1D cylindrical approach performs best. For example, the s-polarized field scattered from a silica-substrate-supported 600 × 180 × 60 nm gold elliptic nanocylinder discretized into 642 dipoles is computed with 0.78 % pattern error and 6.54 % net power error within 294 s, which is 6 times faster than FDTD. Our 1D cylindrical approach takes advantage of parallel processing and also gives transmitted field solutions, which, to the best of our knowledge, is not found in existing tools. We also examine the differences among four polarizability models: Clausius–Mossotti, radiation reaction, lattice dispersion relation, and digitized Green’s function, finding that the radiation reaction dipole model performs best in terms of pattern error, while the digitized Green’s function has the lowest power error.  more » « less
Award ID(s):
2045220
PAR ID:
10509799
Author(s) / Creator(s):
;
Publisher / Repository:
De Gruyter
Date Published:
Journal Name:
Nanophotonics
Volume:
12
Issue:
22
ISSN:
2192-8614
Page Range / eLocation ID:
4157 to 4173
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    A 1D Gaussian expression is derived and used as the 1D E/H incident field in the TF/SF formulation to efficiently generate plane waves in 1D, 2D, and 3D FDTD simulations. The analytic expression is simple, and it eliminates the need for computational resources to store and compute the E/H-field incident arrays and their associated absorbing boundaries. FDTD simulation results at the magic time-step in 1D, 2D, and 3D FDTD show good correlation between plane waves generated by the 1D analytic Gaussian function vs. those generated by 1D FDTD incident arrays. 
    more » « less
  2. Electromagnetic (EM) scattering may be a significant source of degradation in signal and power integrity of high-contrast silicon-on-insulator (SOI) nano-scale interconnects, such as opto-electronic or optical interconnects operating at 100 s of THz where two-dimensional (2D) analytical models of dielectric slab waveguides are often used to approximate scattering loss. In this work, a formulation is presented to relate the scattering (propagation) loss to the scattering parameters (S-parameters) for the smooth waveguide; the results are correlated with results from the finite-difference time-domain (FDTD) method in 2D space. We propose a normalization factor to the previous 2D analytical formulation for the stochastic scattering loss based on physical parameters of waveguides exhibiting random surface roughness under the exponential autocorrelation function (ACF), and validate the results by comparing against numerical experiments via the 2D FDTD method, through simulation of hundreds of rough waveguides; additionally, results are compared to other 2D analytical and previous 3D experimental results. The FDTD environment is described and validated by comparing results of the smooth waveguide against analytical solutions for wave impedance, propagation constant, and S-parameters. Results show that the FDTD model is in agreement with the analytical solution for the smooth waveguide and is a reasonable approximation of the stochastic scattering loss for the rough waveguide. 
    more » « less
  3. Abstract Optical metasurfaces have been widely investigated in recent years as a means to tailor the wavefronts of externally incident light for passive device applications. At the same time, their use in active optoelectronic devices such as light emitters is far less established. This work explores their ability to control the radiation properties of a nearby continuous ensemble of randomly oriented incoherent dipole sources via near‐field interactions. Specifically, a film of colloidal quantum dots is deposited on a plasmonic metasurface consisting of a 1D array of metallic nanoantennas on a metal film. The array is designed to introduce a linear phase profile upon reflection, and a bi‐periodic nanoparticle arrangement is introduced to ensure adequate sampling of the desired phase gradient. Highly directional radiation patterns are correspondingly obtained from the quantum dots at an enhanced emission rate. The underlying radiation mechanism involves the near‐field excitation of surface plasmon polaritons at the metal film, and their selective diffractive scattering by the metasurface into well‐collimated beams along predetermined geometrically tunable directions. These results underscore the distinctive ability of metasurfaces to control radiation properties directly at the source level, which is technologically significant for the continued miniaturization and large‐scale integration of optoelectronic devices. 
    more » « less
  4. Lorenzis, Laura (Ed.)
    Green’s function characterizes a partial differential equation (PDE) and maps its solution in the entire domain as integrals. Finding the analytical form of Green’s function is a non-trivial exercise, especially for a PDE defined on a complex domain or a PDE with variable coefficients. In this paper, we propose a novel boundary integral network to learn the domain independent Green’s function, referred to as BIN-G. We evaluate the Green’s function in the BIN-G using a radial basis function (RBF) kernel-based neural network. We train the BIN-G by minimizing the residual of the PDE and the mean squared errors of the solutions to the boundary integral equations for prescribed test functions. By leveraging the symmetry of the Green’s function and controlling refinements of the RBF kernel near the singularity of the Green function, we demonstrate that our numerical scheme enables fast training and accurate evaluation of the Green’s function for PDEs with variable coefficients. The learned Green’s function is independent of the domain geometries, forcing terms, and boundary conditions in the boundary integral formulation. Numerical experiments verify the desired properties of the method and the expected accuracy for the two-dimensional Poisson and Helmholtz equations with variable coefficients. 
    more » « less
  5. Understanding the computational power of noisy intermediate-scale quantum (NISQ) devices is of both fundamental and practical importance to quantum information science. Here, we address the question of whether error-uncorrected noisy quantum computers can provide computational advantage over classical computers. Specifically, we study noisy random circuit sampling in one dimension (or 1D noisy RCS) as a simple model for exploring the effects of noise on the computational power of a noisy quantum device. In particular, we simulate the real-time dynamics of 1D noisy random quantum circuits via matrix product operators (MPOs) and characterize the computational power of the 1D noisy quantum system by using a metric we call MPO entanglement entropy. The latter metric is chosen because it determines the cost of classical MPO simulation. We numerically demonstrate that for the two-qubit gate error rates we considered, there exists a characteristic system size above which adding more qubits does not bring about an exponential growth of the cost of classical MPO simulation of 1D noisy systems. Specifically, we show that above the characteristic system size, there is an optimal circuit depth, independent of the system size, where the MPO entanglement entropy is maximized. Most importantly, the maximum achievable MPO entanglement entropy is bounded by a constant that depends only on the gate error rate, not on the system size. We also provide a heuristic analysis to get the scaling of the maximum achievable MPO entanglement entropy as a function of the gate error rate. The obtained scaling suggests that although the cost of MPO simulation does not increase exponentially in the system size above a certain characteristic system size, it does increase exponentially as the gate error rate decreases, possibly making classical simulation practically not feasible even with state-of-the-art supercomputers. 
    more » « less