Researcher profile

Steven G. Johnson

Steven G. Johnson contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
19works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2026arXiv

Differentiating through binarized topology changes: Second-order subpixel-smoothed projection

A key challenge in topology optimization (TopOpt) is that manufacturable structures, being inherently binary, are non-differentiable, creating a fundamental tension with gradient-based optimization. The subpixel-smoothed projection (SSP) method addresses this issue by smoothing sharp interfaces at the subpixel level through a first-order expansion of the filtered field. However, SSP does not guarantee differentiability under topology changes, such as the merging of two interfaces, and therefore violates the convergence guarantees of many popular gradient-based optimization algorithms. We overcome this limitation by regularizing SSP with the Hessian of the filtered field, resulting in a twice-differentiable projected density during such transitions, while still guaranteeing an almost-everywhere binary structure. We demonstrate the effectiveness of our second-order SSP (SSP2) methodology on both thermal and photonic problems, showing that SSP2 has faster convergence than SSP for connectivity-dominant cases -- where frequent topology changes occur -- while exhibiting comparable performance otherwise. Beyond improving convergence guarantees for CCSA optimizers, SSP2 enables the use of a broader class of optimization algorithms with stronger theoretical guarantees, such as interior-point methods. Since SSP2 adds minimal complexity relative to SSP or traditional projection schemes, it can be used as a drop-in replacement in existing TopOpt codes.

preprint2022arXiv

Analytical criteria for designing multiresonance filters in scattering systems, with application to microwave metasurfaces

We present general analytical criteria for the design of lossless reciprocal two-port systems, which exhibit prescribed scattering spectra $S(ω)$ satisfying $S_{22}(ω)=e^{iφ}S_{11}(ω)$, including symmetric ($S_{22}=S_{11}$) or "antimetric" ($S_{22}=-S_{11}$) responses, such as standard filters (Butterworth, Chebyshev, elliptic, etc.). We show that the non-normalized resonant (quasi-normal) modes (QNMs) of all such two-port systems couple to the input and output ports with specific unitary ratios, whose relative signs determine the position of the scattering zeros on the real frequency axis. This allows us to obtain design criteria assigning values to the poles, background response, and QNM-to-ports coupling coefficients. Filter devices can then be designed via a well-conditioned nonlinear optimization (or root-finding) problem using a numerical eigensolver. As an application, we design multiple microwave metasurfaces configured for polarization-preserving transmission, reflective polarization conversion, or diffractive "perfect anomalous reflection", to realize filters that precisely match standard bandpass or bandstop filters of various types, orders and bandwidths, with focus on the best-performing elliptic filters.

preprint2022arXiv

Free-electron-light interactions in nanophotonics

When impinging on optical structures or passing in their vicinity, free electrons can spontaneously emit electromagnetic radiation, a phenomenon generally known as cathodoluminescence. Free-electron radiation comes in many guises: Cherenkov, transition, and Smith-Purcell radiation, but also electron scintillation, commonly referred to as incoherent cathodoluminescence. While those effects have been at the heart of many fundamental discoveries and technological developments in high-energy physics in the past century, their recent demonstration in photonic and nanophotonic systems has attracted a lot of attention. Those developments arose from predictions that exploit nanophotonics for novel radiation regimes, now becoming accessible thanks to advances in nanofabrication. In general, the proper design of nanophotonic structures can enable shaping, control, and enhancement of free-electron radiation, for any of the above-mentioned effects. Free-electron radiation in nanophotonics opens the way to promising applications, such as widely-tunable integrated light sources from x-ray to THz frequencies, miniaturized particle accelerators, and highly sensitive high-energy particle detectors. Here, we review the emerging field of free-electron radiation in nanophotonics. We first present a general, unified framework to describe free-electron light-matter interaction in arbitrary nanophotonic systems. We then show how this framework sheds light on the physical underpinnings of many methods in the field used to control and enhance free-electron radiation. Namely, the framework points to the central role played by the photonic eigenmodes in controlling the output properties of free-electron radiation (e.g., frequency, directionality, and polarization). [... see full abstract in paper]

preprint2022arXiv

Inverse-Designed Meta-Optics with Spectral-Spatial Engineered Response to Mimic Color Perception

Meta-optics have rapidly become a major research field within the optics and photonics community, strongly driven by the seemingly limitless opportunities made possible by controlling optical wavefronts through interaction with arrays of sub-wavelength scatterers. As more and more modalities are explored, the design strategies to achieve desired functionalities become increasingly demanding, necessitating more advanced design techniques. Herein, the inverse-design approach is utilized to create a set of single-layer meta-optics that simultaneously focus light and shape the spectra of focused light without using any filters. Thus, both spatial and spectral properties of the meta-optics are optimized, resulting in spectra that mimic the color matching functions of the CIE 1931 XYZ color space, which links the distributions of wavelengths in light and the color perception of a human eye. Experimental demonstrations of these meta-optics show qualitative agreement with the theoretical predictions and help elucidate the focusing mechanism of these devices.

preprint2022arXiv

Multifidelity deep neural operators for efficient learning of partial differential equations with application to fast inverse design of nanoscale heat transport

Deep neural operators can learn operators mapping between infinite-dimensional function spaces via deep neural networks and have become an emerging paradigm of scientific machine learning. However, training neural operators usually requires a large amount of high-fidelity data, which is often difficult to obtain in real engineering problems. Here, we address this challenge by using multifidelity learning, i.e., learning from multifidelity datasets. We develop a multifidelity neural operator based on a deep operator network (DeepONet). A multifidelity DeepONet includes two standard DeepONets coupled by residual learning and input augmentation. Multifidelity DeepONet significantly reduces the required amount of high-fidelity data and achieves one order of magnitude smaller error when using the same amount of high-fidelity data. We apply a multifidelity DeepONet to learn the phonon Boltzmann transport equation (BTE), a framework to compute nanoscale heat transport. By combining a trained multifidelity DeepONet with genetic algorithm or topology optimization, we demonstrate a fast solver for the inverse design of BTE problems.

preprint2021arXiv

A general framework for scintillation in nanophotonics

Bombardment of materials by high-energy particles (e.g., electrons, nuclei, X- and $γ$-ray photons) often leads to light emission, known generally as scintillation. Scintillation is ubiquitous and enjoys widespread applications in many areas such as medical imaging, X-ray non-destructive inspection, night vision, electron microscopy, and high-energy particle detectors. A large body of research focuses on finding new materials optimized for brighter, faster, and more controlled scintillation. Here, we develop a fundamentally different approach based on integrating nanophotonic structures into scintillators to enhance their emission. To start, we develop a unified and ab initio theory of nanophotonic scintillators that accounts for the key aspects of scintillation: the energy loss by high-energy particles, as well as the light emission by non-equilibrium electrons in arbitrary nanostructured optical systems. This theoretical framework allows us, for the first time, to experimentally demonstrate nearly an order-of-magnitude enhancement of scintillation, in both electron-induced, and X-ray-induced scintillation. Our theory also allows the discovery of structures that could eventually achieve several orders-of-magnitude scintillation enhancement. The framework and results shown here should enable the development of a new class of brighter, faster, and higher-resolution scintillators with tailored and optimized performances - with many potential applications where scintillators are used.

preprint2021arXiv

End-to-end metasurface inverse design for single-shot multi-channel imaging

We introduce end-to-end metaoptics inverse design for multi-channel imaging: reconstruction of depth, spectral and polarization channels from a single-shot monochrome image. The proposed technique integrates a single-layer metasurface frontend with an efficient Tikhonov reconstruction backend, without any additional optics except a grayscale sensor. Our method yields multi-channel imaging by spontaneous demultiplexing: the metaoptics front-end separates different channels into distinct spatial domains whose locations on the sensor are optimally discovered by the inverse-design algorithm. We present large-area metasurface designs, compatible with standard lithography, for multi-spectral imaging, depth-spectral imaging, and ``all-in-one'' spectro-polarimetric-depth imaging with robust reconstruction performance ($\lesssim 10\%$ error with 1\% detector noise). In contrast to neural networks, our framework is physically interpretable and does not require large training sets. It can be used to reconstruct arbitrary three-dimensional scenes with full multi-wavelength spectra and polarization textures.

preprint2021arXiv

Inverse design enables large-scale high-performance meta-optics reshaping virtual reality

Meta-optics has achieved major breakthroughs in the past decade; however, conventional forward design faces challenges as functionality complexity and device size scale up. Inverse design aims at optimizing meta-optics design but has been currently limited by expensive brute-force numerical solvers to small devices, which are also difficult to realize experimentally. Here, we present a general inverse design framework for aperiodic large-scale complex meta-optics in three dimensions, which alleviates computational cost for both simulation and optimization via a fast-approximate solver and an adjoint method, respectively. Our framework naturally accounts for fabrication constraints via a surrogate model. In experiments, we demonstrate, for the first time, aberration-corrected metalenses working in the visible with high numerical aperture, poly-chromatic focusing, and large diameter up to centimeter scale. Such large-scale meta-optics opens a new paradigm for applications, and we demonstrate its potential for future virtual-reality platforms by using a meta-eyepiece and a laser back-illuminated micro-Liquid Crystal Display.

preprint2021arXiv

Modified discrete Laguerre polynomials for efficient computation of exponentially bounded Matsubara sums

We develop a new type of orthogonal polynomial, the modified discrete Laguerre (MDL) polynomials, designed to accelerate the computation of bosonic Matsubara sums in statistical physics. The MDL polynomials lead to a rapidly convergent Gaussian "quadrature" scheme for Matsubara sums, and more generally for any sum $F(0)/2 + F(h) + F(2h) + \cdots$ of exponentially decaying summands $F(nh) = f(nh)e^{-nhs}$ where $hs>0$. We demonstrate this technique for computation of finite-temperature Casimir forces arising from quantum field theory, where evaluation of the summand $F$ requires expensive electromagnetic simulations. A key advantage of our scheme, compared to previous methods, is that the convergence rate is nearly independent of the spacing $h$ (proportional to the thermodynamic temperature). We also prove convergence for any polynomially decaying $F$.

preprint2021arXiv

Physics-informed neural networks with hard constraints for inverse design

Inverse design arises in a variety of areas in engineering such as acoustic, mechanics, thermal/electronic transport, electromagnetism, and optics. Topology optimization is a major form of inverse design, where we optimize a designed geometry to achieve targeted properties and the geometry is parameterized by a density function. This optimization is challenging, because it has a very high dimensionality and is usually constrained by partial differential equations (PDEs) and additional inequalities. Here, we propose a new deep learning method -- physics-informed neural networks with hard constraints (hPINNs) -- for solving topology optimization. hPINN leverages the recent development of PINNs for solving PDEs, and thus does not rely on any numerical PDE solver. However, all the constraints in PINNs are soft constraints, and hence we impose hard constraints by using the penalty method and the augmented Lagrangian method. We demonstrate the effectiveness of hPINN for a holography problem in optics and a fluid problem of Stokes flow. We achieve the same objective as conventional PDE-constrained optimization methods based on adjoint methods and numerical PDE solvers, but find that the design obtained from hPINN is often simpler and smoother for problems whose solution is not unique. Moreover, the implementation of inverse design with hPINN can be easier than that of conventional methods.

preprint2021arXiv

Towards 3D-Printed Inverse-Designed Metaoptics

Optical metasurfaces have been heralded as the platform to integrate multiple functionalities in a compact form-factor, potentially replacing bulky components. A central stepping stone towards realizing this promise is the demonstration of multifunctionality under several constraints (e.g. at multiple incident wavelengths and/or angles) in a single device -- an achievement being hampered by design limitations inherent to single-layer planar geometries. Here, we propose a general framework for the inverse design of volumetric 3D metaoptics via topology optimization, showing that even few-wavelength thick devices can achieve high-efficiency multifunctionality. We embody our framework in multiple closely-spaced patterned layers of a low-index polymer. We experimentally demonstrate our approach with an inverse-designed 3d-printed light concentrator working at five different non-paraxial angles of incidence. Our framework paves the way towards realizing multifunctional ultra-compact 3D nanophotonic devices.

preprint2020arXiv

Approaching the upper limits of the local density of states via optimized metallic cavities

By computational optimization of air-void cavities in metallic substrates, we show that the local density of states (LDOS) can reach within a factor of $\approx 10$ of recent theoretical upper limits, and within a factor $\approx 4$ for the single-polarization LDOS, demonstrating that the theoretical limits are nearly attainable. Optimizing the total LDOS results in a spontaneous symmetry breaking where it is preferable to couple to a specific polarization. Moreover, simple shapes such as optimized cylinders attain nearly the performance of complicated many-parameter optima, suggesting that only one or two key parameters matter in order to approach the theoretical LDOS bounds for metallic resonators.

preprint2020arXiv

Computational inverse design for ultra-compact single-piece metalenses free of chromatic and angular aberration

We present full-Maxwell topology-optimization design of a single-piece multlayer metalens, about 10 wavelengths~$λ$ in thickness, that simultaneously focuses over a $60^\circ$ angular range and a 23\% spectral bandwidth without suffering chromatic or angular aberration, a "plan-achromat." At all angles and frequencies it achieves diffraction-limited focusing (Strehl ratio $> 0.8$) and absolute focusing efficiency $> 50$\%. Both 2D and 3D axi-symmetric designs are presented, optimized over $\sim 10^5$ degrees of freedom. We also demonstrate shortening the lens-to-sensor distance while producing the same image as for a longer "virtual" focal length and maintaining plan-achromaticity. These proof-of-concept designs demonstrate the ultra-compact multi-functionality that can be achieved by exploiting the full wave physics of subwavelength designs, and motivate future work on design and fabrication of multi-layer meta-optics.

preprint2020arXiv

Inverse designed metalenses with extended depth of focus

Extended depth of focus (EDOF) lenses are important for various applications in computational imaging and microscopy. In addition to enabling novel functionalities, EDOF lenses can alleviate the need for stringent alignment requirements for imaging systems. Existing EDOF lenses, however, are often inefficient or produce an asymmetric point spread function (PSF) that blurs images. Inverse design of nanophotonics, including metasurfaces, has generated strong interest in recent years owing to its potential for generating exotic and innovative optical elements, which are generally difficult to model intuitively. Using adjoint optimization-based inverse electromagnetic design, in this paper, we designed a cylindrical metasurface lens operating at ~ 625nm with a depth of focus exceeding that of an ordinary lens. We validated our design by nanofabrication and optical characterization of silicon nitride metasurface lenses (with lateral dimension of 66.66 μm) with three different focal lengths (66.66 μm,100 μm,133.33 μm). The focusing efficiencies of the fabricated extended depth of focus metasurface lenses are similar to those of traditional metalenses.

preprint2020arXiv

Is single-mode lasing possible in an infinite periodic system?

In this Letter, we present a rigorous method to study the stability of periodic lasing systems. In a linear model, the presence of a continuum of modes (with arbitrarily close lasing thresholds) gives the impression that stable single-mode lasing cannot be maintained in the limit of an infinite system. However, we show that nonlinear effects of the Maxwell-Bloch equations can lead to stable systems near threshold given a simple stability condition on the sign of the laser detuning compared to the band curvature. We examine band-edge (1d) and bound-in-continuum (2d) lasing modes and validate our stability results against time-domain simulations.

preprint2020arXiv

Modeling lasers and saturable absorbers via multilevel atomic media in the Meep FDTD software: Theory and implementation

This technical note describes the physical model, numerical implementation, and validation of multilevel atomic media for lasers and saturable absorbers in Meep: a free/open-source finite-difference time-domain (FDTD) software package for electromagnetics simulation. Simulating multilevel media in the time domain involves coupling rate equations for the populations of electronic energy levels with Maxwell's equations via a generalization of the Maxwell--Bloch equations. We describe the underlying equations and their implementation using a second-order discretization scheme, and also demonstrate their equivalence to a quantum density-matrix model. The Meep implementation is validated using a separate FDTD density-matrix model as well as a frequency-domain solver based on steady-state ab-initio laser theory (SALT).

preprint2020arXiv

Response to Comment on Widely tunable compact terahertz gas lasers

We recently demonstrated a widely tunable THz molecular laser and reported mathematical formulas and a table for comparing how various molecules would perform as such lasers (Chevalier et al., Science, 15 November 2019, p. 856-860). Here we correct the value of a single parameter used to calculate the table (see Erratum for Chevalier et al.), thereby eliminating the concerns raised by Lampin and Barbieri (Lampin et al., arXiv:2004.04422). We also show that our simplified model for the output THz power is a better approximation than the alternative one proposed in the technical comment.

preprint2019arXiv

Inverse Design of Nanoparticles for Enhanced Raman Scattering

We show that topology optimization (TO) of metallic resonators can lead to $\sim 10^2\times$ improvement in surface-enhanced Raman scattering (SERS) efficiency compared to traditional resonant structures such as bowtie antennas. TO inverse design leads to surprising structures very different from conventional designs, which simultaneously optimize focusing of the incident wave and emission from the Raman dipole. We consider isolated metallic particles as well as more complicated configurations such as periodic surfaces or resonators coupled to dielectric waveguides, and the benefits of TO are even greater in the latter case. Our results are motivated by recent rigorous upper bounds to Raman scattering enhancement, and shed light on the extent to which these bounds are achievable.

preprint2019arXiv

Overlapping domains for topology optimization of large-area metasurfaces

We introduce an overlapping-domain approach to large-area metasurface design, in which each simulated domain consists of a unit cell and overlapping regions from the neighboring cells plus PML absorbers. We show that our approach generates greatly improved metalens quality compared to designs produced using a locally periodic approximation, thanks to $\sim 10\times$ better accuracy with similar computational cost. We use the new approach with topology optimization to design large-area ($200λ$) high-NA (0.71) multichrome and broadband achromatic lenses with high focusing efficiency ($\sim 50\%$), greatly improving upon previously reported works.