skip to main content

Explore Research Products in the NSF-PAR

Title: Fast Inverse Design of 3D Nanophotonic Devices Using Boundary Integral Methods
Recent developments in the computational automated design of electromagnetic devices, otherwise known as inverse design, have significantly enhanced the design process for nanophotonic systems. Inverse design can both reduce design time considerably and lead to high-performance, nonintuitive structures that would otherwise have been impossible to develop manually. Despite the successes enjoyed by structure optimization techniques, most approaches leverage electromagnetic solvers that require significant computational resources and suffer from slow convergence and numerical dispersion. Recently, a fast simulation and boundary-based inverse design approach based on boundary integral equations was demonstrated for two-dimensional nanophotonic problems. In this work, we introduce a new full-wave three-dimensional simulation and boundary-based optimization framework for nanophotonic devices also based on boundary integral methods, which achieves high accuracy even at coarse mesh discretizations while only requiring modest computational resources. The approach has been further accelerated by leveraging GPU computing, a sparse block-diagonal preconditioning strategy, and a matrix-free implementation of the discrete adjoint method. As a demonstration, we optimize three different devices: a 1:2 1550 nm power splitter and two nonadiabatic mode-preserving waveguide tapers. To the best of our knowledge, the tapers, which span 40 wavelengths in the silicon material, are the largest silicon photonic waveguiding devices to have been optimized using full-wave 3D solution of Maxwell’s equations.  more » « less
Award ID(s):
1849965 2047433
NSF-PAR ID:
10382209
Author(s) / Creator(s):
;
Date Published:
Journal Name:
ACS Photonics
ISSN:
2330-4022
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ( , Advanced Photonics Research)
    Optical metasurfaces consist of densely arranged unit cells that manipulate light through various light confinement and scattering processes. Due to its unique advantages, such as high performance, small form factor and easy integration with semiconductor devices, metasurfaces have been gathering increasing attention in fields such as displays, imaging, sensing and optical computation. Despite advances in fabrication and characterization, a viable design prediction for suitable optical response remains challenging for complex optical metamaterial systems. The computation cost required to obtain the optimal design exponentially grows as the design complexity increases. Furthermore, the design prediction is challenging since the inverse problem is often ill-posed. In recent years, deep learning (DL) methods have shown great promise in the area of inverse design. Inspired by this and the capability of DL to produce fast inference, we introduce a physics-informed DL framework to expedite the computation for the inverse design of metasurfaces. Addition of the physics-based constraints improve generalizability of the DL model while reducing data burden. Our approach introduces a tandem DL architecture with physics-based learning to alleviate the nonuniqueness issue by selecting designs that are scientifically consistent, with low error in design prediction and accurate reconstruction of optical responses. To prove the concept, we focus on the inverse design of a representative plasmonic device that consists of metal gratings deposited on a dielectric film on top of a metal substrate. The optical response of the device is determined by the geometrical dimensions as well as the material properties. The training and testing data are obtained through Rigorous Coupled-Wave Analysis (RCWA), while the physics-based constraint is derived from solving the electromagnetic (EM) wave equations for a simplified homogenized model. We consider the prediction of design for the optical response of a single wavelength incident or a spectrum of wavelength in the visible light range. Our model converges with an accuracy up to 97% for inverse design prediction with the optical response for the visible light spectrum as input. The model is also able to predict design with accuracy up to 96% and optical response reconstruction accuracy of 99% for optical response of a single wavelength of light as input. 
    more » « less
  2. ; ( , Annual Review of Heat Transfer)
    Chu, Wilson (Ed.)
    Two-dimensional materials (e.g., graphene and transition metal dichalcogenides) and their heterostructures have enormous applications in electrochemical energy storage systems such as batteries. A comprehensive and solid understanding of these materials’ thermal transport and mechanism is essential for practical device design. Several advanced experimental techniques have been developed to measure the intrinsic thermal conductivity of materials. However, experiments have challenges in providing improved control and characterization of complex structures, especially for low-dimensional materials. Theoretical and simulation tools, such as first-principles calculations, Boltzmann transport equations, molecular dynamics simulations, lattice dynamics simulation, and nonequilibrium Green’s function, provide reliable predictions of thermal conductivity and physical insights to understand the underlying thermal transport mechanism in materials. However, doing these calculations requires high computational resources. The development of new materials synthesis technology and fast-growing demand for rapid and accurate prediction of physical properties requires novel computational approaches. The machine learning method provides a promising solution to address such needs. This review details the recent development in atomistic/molecular studies and machine learning of thermal transport in two-dimensional materials. The paper also addresses the latest significant experimental advances. However, designing the best two-dimensional materials-based heterostructures is like a multivariate optimization problem. For example, a particular heterostructure may be suitable for thermal transport but can have lower mechanical strength/stability. For bilayer and multilayer structures, the interlayer distance may influence the thermal transport properties and interlayer strength. Therefore, the last part of this review addresses the future research direction in two-dimensional materials-based heterostructure design for thermal transport in energy storage systems. 
    more » « less
  3. ( , HOUSEHOLDER SYMPOSIUM XXI)
    Embedding properties of network realizations of dissipative reduced order models Jörn Zimmerling, Mikhail Zaslavsky,Rob Remis, Shasri Moskow, Alexander Mamonov, Murthy Guddati, Vladimir Druskin, and Liliana Borcea Mathematical Sciences Department, Worcester Polytechnic Institute https://www.wpi.edu/people/vdruskin Abstract Realizations of reduced order models of passive SISO or MIMO LTI problems can be transformed to tridiagonal and block-tridiagonal forms, respectively, via dierent modications of the Lanczos algorithm. Generally, such realizations can be interpreted as ladder resistor-capacitor-inductor (RCL) networks. They gave rise to network syntheses in the rst half of the 20th century that was at the base of modern electronics design and consecutively to MOR that tremendously impacted many areas of engineering (electrical, mechanical, aerospace, etc.) by enabling ecient compression of the underlining dynamical systems. In his seminal 1950s works Krein realized that in addition to their compressing properties, network realizations can be used to embed the data back into the state space of the underlying continuum problems. In more recent works of the authors Krein's ideas gave rise to so-called nite-dierence Gaussian quadrature rules (FDGQR), allowing to approximately map the ROM state-space representation to its full order continuum counterpart on a judicially chosen grid. Thus, the state variables can be accessed directly from the transfer function without solving the full problem and even explicit knowledge of the PDE coecients in the interior, i.e., the FDGQR directly learns" the problem from its transfer function. This embedding property found applications in PDE solvers, inverse problems and unsupervised machine learning. Here we show a generalization of this approach to dissipative PDE problems, e.g., electromagnetic and acoustic wave propagation in lossy dispersive media. Potential applications include solution of inverse scattering problems in dispersive media, such as seismic exploration, radars and sonars. To x the idea, we consider a passive irreducible SISO ROM fn(s) = Xn j=1 yi s + σj , (62) assuming that all complex terms in (62) come in conjugate pairs. We will seek ladder realization of (62) as rjuj + vj − vj−1 = −shˆjuj , uj+1 − uj + ˆrj vj = −shj vj , (63) for j = 0, . . . , n with boundary conditions un+1 = 0, v1 = −1, and 4n real parameters hi, hˆi, ri and rˆi, i = 1, . . . , n, that can be considered, respectively, as the equivalent discrete inductances, capacitors and also primary and dual conductors. Alternatively, they can be viewed as respectively masses, spring stiness, primary and dual dampers of a mechanical string. Reordering variables would bring (63) into tridiagonal form, so from the spectral measure given by (62 ) the coecients of (63) can be obtained via a non-symmetric Lanczos algorithm written in J-symmetric form and fn(s) can be equivalently computed as fn(s) = u1. The cases considered in the original FDGQR correspond to either (i) real y, θ or (ii) real y and imaginary θ. Both cases are covered by the Stieltjes theorem, that yields in case (i) real positive h, hˆ and trivial r, rˆ, and in case (ii) real positive h,r and trivial hˆ,rˆ. This result allowed us a simple interpretation of (62) as the staggered nite-dierence approximation of the underlying PDE problem [2]. For PDEs in more than one variables (including topologically rich data-manifolds), a nite-dierence interpretation is obtained via a MIMO extensions in block form, e.g., [4, 3]. The main diculty of extending this approach to general passive problems is that the Stieltjes theory is no longer applicable. Moreover, the tridiagonal realization of a passive ROM transfer function (62) via the ladder network (63) cannot always be obtained in port-Hamiltonian form, i.e., the equivalent primary and dual conductors may change sign [1]. 100 Embedding of the Stieltjes problems, e.g., the case (i) was done by mapping h and hˆ into values of acoustic (or electromagnetic) impedance at grid cells, that required a special coordinate stretching (known as travel time coordinate transform) for continuous problems. Likewise, to circumvent possible non-positivity of conductors for the non-Stieltjes case, we introduce an additional complex s-dependent coordinate stretching, vanishing as s → ∞ [1]. This stretching applied in the discrete setting induces a diagonal factorization, removes oscillating coecients, and leads to an accurate embedding for moderate variations of the coecients of the continuum problems, i.e., it maps discrete coecients onto the values of their continuum counterparts. Not only does this embedding yields an approximate linear algebraic algorithm for the solution of the inverse problems for dissipative PDEs, it also leads to new insight into the properties of their ROM realizations. We will also discuss another approach to embedding, based on Krein-Nudelman theory [5], that results in special data-driven adaptive grids. References [1] Borcea, Liliana and Druskin, Vladimir and Zimmerling, Jörn, A reduced order model approach to inverse scattering in lossy layered media, Journal of Scientic Computing, V. 89, N1, pp. 136,2021 [2] Druskin, Vladimir and Knizhnerman, Leonid, Gaussian spectral rules for the three-point second dierences: I. A two-point positive denite problem in a semi-innite domain, SIAM Journal on Numerical Analysis, V. 37, N 2, pp.403422, 1999 [3] Druskin, Vladimir and Mamonov, Alexander V and Zaslavsky, Mikhail, Distance preserving model order reduction of graph-Laplacians and cluster analysis, Druskin, Vladimir and Mamonov, Alexander V and Zaslavsky, Mikhail, Journal of Scientic Computing, V. 90, N 1, pp 130, 2022 [4] Druskin, Vladimir and Moskow, Shari and Zaslavsky, Mikhail LippmannSchwingerLanczos algorithm for inverse scattering problems, Inverse Problems, V. 37, N. 7, 2021, [5] Mark Adolfovich Nudelman The Krein String and Characteristic Functions of Maximal Dissipative Operators, Journal of Mathematical Sciences, 2004, V 124, pp 49184934 Go back to Plenary Speakers Go back to Speakers Go back 
    more » « less
  4. ; ; ; ( , ACS Photonics)
    Wave-based analog computing in the forms of inversedesigned metastructures and the meshes of Mach−Zehnder interferometers (MZI) have recently received considerable attention due to their capability in emulating linear operators, performing vector-matrix multiplication, inverting matrices, and solving integral and differential equations via electromagnetic wave interaction and manipulation in such structures. Here, we combine these two platforms to propose a wave-based metadevice that can compute scattered fields in electromagnetic forward scattering problems. The proposed device consists of two subsystems: a set of reconfigurable couplers with a proper feedback system and an inverse-designed inhomogeneous material block. The first subsystem computes the magnitude and phase of the dipole polarization induced in the scatterers when illuminated with a given incident wave (matrix inversion). The second subsystem computes the magnitude and phase of the scattered fields at given detection points (vector-matrix multiplication). We discuss the functionality of this metadevice, and through several examples, we theoretically evaluate its performance by comparing the simulation results of this device with fullwave numerical simulations and numerically evaluated matrix inversion. We also highlight that since the first section is reconfigurable, the proposed device can be used for different permittivity distributions of the scatterer and incident excitations without changing the inverse-designed section. Our proposed device may provide a versatile platform for rapid computation in various scattering scenarios. 
    more » « less
  5. ; ; ; ( , AIP Advances)

    High power sources of electromagnetic energy often require complicated structures to support electromagnetic modes and shape electromagnetic fields to maximize the coupling of the field energy to intense relativistic electron beams. Geometric fidelity is critical to the accurate simulation of these High Power Electromagnetic (HPEM) sources. Here, we present a fast and geometrically flexible approach to calculate the solution to Maxwell’s equations in vector potential form under the Lorenz gauge. The scheme is an implicit, linear-time, high-order, A-stable method that is based on the method of lines transpose (MOLT). As presented, the method is fourth order in time and second order in space, but the A-stable formulation could be extended to both high order in time and space. An O(n) fast convolution is employed for space-integration. The main focus of this work is to develop an approach to impose perfectly electrically conducting (PEC) boundary conditions in MOLT by extending our past work on embedded boundary methods. As the method is A-stable, it does not suffer from small time step limitations that are found in explicit finite difference time domain methods when using either embedded boundary or cut-cell methods to capture geometry. This is a major advance for the simulation of HPEM devices. While there is no conceptual limitation to develop this in 3D, our initial work has centered on 2D. The extension to 3D requires validation that the proposed fixed point iteration will converge and is the subject of our follow-up work. The eventual goal is to combine this method with particle methods for the simulations of plasma. In the current work, the scheme is evaluated for EM wave propagation within an object that is bounded by PEC. The consistency and performance of the scheme are confirmed using the ping test and frequency mode analysis for rotated square cavities—a standard test in the HPEM community. We then demonstrate the diffraction Q value test and the use of this method for simulating an A6 magnetron. The ability to handle both PEC and open boundaries in a standard device test problem, such as the A6, gives confidence on the robustness of this new method.

     
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email