Source author record

Steven G. Johnson

Steven G. Johnson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

86works

25topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

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.

preprint2016arXiv

Fundamental limits to optical response in absorptive systems

At visible and infrared frequencies, metals show tantalizing promise for strong subwavelength resonances, but material loss typically dampens the response. We derive fundamental limits to the optical response of absorptive systems, bounding the largest enhancements possible given intrinsic material losses. Through basic conservation-of-energy principles, we derive geometry-independent limits to per-volume absorption and scattering rates, and to local-density-of-states enhancements that represent the power radiated or expended by a dipole near a material body. We provide examples of structures that approach our absorption and scattering limits at any frequency, by contrast, we find that common "antenna" structures fall far short of our radiative LDOS bounds, suggesting the possibility for significant further improvement. Underlying the limits is a simple metric, $|χ|^2 / \operatorname{Im} χ$ for a material with susceptibility $χ$, that enables broad technological evaluation of lossy materials across optical frequencies.

preprint2016arXiv

Interaction-induced mode switching in steady-state microlasers

We demonstrate that due to strong modal interactions through cross-gain saturation, the onset of a new lasing mode can switch off an existing mode via a negative power slope. In this process of interaction-induced mode switching (IMS) the two involved modes maintain their identities, i.e. they do not change their spatial field patterns or lasing frequencies. For a fixed pump profile, a simple analytic criterion for the occurrence of IMS is given in terms of their self- and cross-interaction coefficients and non-interacting thresholds, which is verified for the example of a two-dimensional microdisk laser. When the spatial pump profile is varied as the pump power is increased, IMS can be induced even when it would not occur with a fixed pump profile, as we show for two coupled laser cavities. Our findings apply to steady-state lasing and are hence different from dynamical mode switching or hopping. IMS may have potential applications in robust and flexible all-optical switching.

preprint2016arXiv

Perfect Single-Sided Radiation and Absorption without Mirrors

Highly directional radiation from photonic structures is important for many applications, including high power photonic crystal surface emitting lasers, grating couplers, and light detection and ranging devices. However, previous dielectric, few-layer designs only achieved moderate asymmetry ratios, and a fundamental understanding of bounds on asymmetric radiation from arbitrary structures is still lacking. Here, we show that breaking the 180$^\circ$ rotational symmetry of the structure is crucial for achieving highly asymmetric radiation. We develop a general temporal coupled-mode theory formalism to derive bounds on the asymmetric decay rates to the top and bottom of a photonic crystal slab for a resonance with arbitrary in-plane wavevector. Guided by this formalism, we show that infinite asymmetry is still achievable even without the need of back-reflection mirrors, and we provide numerical examples of designs that achieve asymmetry ratios exceeding $10^4$. The emission direction can also be rapidly switched from top to bottom by tuning the wavevector or frequency. Furthermore, we show that with the addition of weak material absorption loss, such structures can be used to achieve perfect absorption with single-sided illumination, even for single-pass absorption rates less than $0.5\%$ and without back-reflection mirrors. Our work provides new design principles for achieving highly directional radiation and perfect absorption in photonics.

preprint2015arXiv

Ab-initio multimode linewidth theory for arbitrary inhomogeneous laser cavities

We present a multimode laser-linewidth theory for arbitrary cavity structures and geometries that contains nearly all previously known effects and also finds new nonlinear and multimode corrections, e.g. a bad-cavity correction to the Henry $α$ factor and a multimode Schawlow--Townes relation (each linewidth is proportional to a sum of inverse powers of all lasing modes). Our theory produces a quantitatively accurate formula for the linewidth, with no free parameters, including the full spatial degrees of freedom of the system. Starting with the Maxwell--Bloch equations, we handle quantum and thermal noise by introducing random currents whose correlations are given by the fluctuation--dissipation theorem. We derive coupled-mode equations for the lasing-mode amplitudes and obtain a formula for the linewidths in terms of simple integrals over the steady-state lasing modes.

preprint2015arXiv

All-optical three-dimensional electron pulse compression

We propose an all-optical, three-dimensional electron pulse compression scheme in which Hermite-Gaussian optical modes are used to fashion a three-dimensional optical trap in the electron pulse's rest frame. We show that the correct choices of optical incidence angles are necessary for optimal compression. We obtain analytical expressions for the net impulse imparted by Hermite-Gaussian free-space modes of arbitrary order. Although we focus on electrons, our theory applies to any charged particle and any particle with non-zero polarizability in the Rayleigh regime. We verify our theory numerically using exact solutions to Maxwell's equations for first-order Hermite-Gaussian beams, demonstrating single-electron pulse compression factors of $>10^{2}$ in both longitudinal and transverse dimensions with experimentally realizable optical pulses. The proposed scheme is useful in ultrafast electron imaging for both single- and multi-electron pulse compression, and as a means of circumventing temporal distortions in magnetic lenses when focusing ultrashort electron pulses.

preprint2015arXiv

Cavity-enhanced second harmonic generation via nonlinear-overlap optimization

We describe an approach based on topology optimization that enables automatic discovery of wavelength-scale photonic structures for achieving high-efficiency second-harmonic generation (SHG). A key distinction from previous formulation and designs that seek to maximize Purcell factors at individual frequencies is that our method not only aims to achieve frequency matching (across an entire octave) and large radiative lifetimes, but also optimizes the equally important nonlinear--coupling figure of merit $\barβ$, involving a complicated spatial overlap-integral between modes. We apply this method to the particular problem of optimizing micropost and grating-slab cavities (one-dimensional multilayered structures) and demonstrate that a variety of material platforms can support modes with the requisite frequencies, large lifetimes $Q > 10^4$, small modal volumes $\sim (λ/n)^3$, and extremely large $\barβ \gtrsim 10^{-2}$, leading to orders of magnitude enhancements in SHG efficiency compared to state of the art photonic designs. Such giant $\barβ$ alleviate the need for ultra-narrow linewidths and thus pave the way for wavelength-scale SHG devices with faster operating timescales and higher tolerance to fabrication imperfections.

preprint2015arXiv

Design of diamond microcavities for single photon frequency down-conversion

We propose monolithic diamond cavities that can be used to convert color-center Fock-state single photons from emission wavelengths to telecommunication bands. We present a detailed theoretical description of the conversion process, analyzing important practical concerns such as nonlinear phase shifts and frequency mismatch. Our analysis predicts sustainable power requirements ($ \lesssim 1~\mathrm{W}$) for a chipscale nonlinear device with high conversion efficiencies.

preprint2015arXiv

Fluctuating volume-current formulation of electromagnetic fluctuations in inhomogeneous media: incandecence and luminescence in arbitrary geometries

We describe a fluctuating volume--current formulation of electromagnetic fluctuations that extends our recent work on heat exchange and Casimir interactions between arbitrarily shaped homogeneous bodies [Phys. Rev. B. 88, 054305] to situations involving incandescence and luminescence problems, including thermal radiation, heat transfer, Casimir forces, spontaneous emission, fluorescence, and Raman scattering, in inhomogeneous media. Unlike previous scattering formulations based on field and/or surface unknowns, our work exploits powerful techniques from the volume--integral equation (VIE) method, in which electromagnetic scattering is described in terms of volumetric, current unknowns throughout the bodies. The resulting trace formulas (boxed equations) involve products of well-studied VIE matrices and describe power and momentum transfer between objects with spatially varying material properties and fluctuation characteristics. We demonstrate that thanks to the low-rank properties of the associatedmatrices, these formulas are susceptible to fast-trace computations based on iterative methods, making practical calculations tractable. We apply our techniques to study thermal radiation, heat transfer, and fluorescence in complicated geometries, checking our method against established techniques best suited for homogeneous bodies as well as applying it to obtain predictions of radiation from complex bodies with spatially varying permittivities and/or temperature profiles.

preprint2015arXiv

Optimization of sharp and viewing-angle-independent structural color

Structural coloration produces some of the most brilliant colors in nature and has many applications. However, the two competing properties of narrow bandwidth and broad viewing angle have not been achieved simultaneously in previous studies. Here, we use numerical optimization to discover geometries where a sharp 7% bandwidth in scattering is achieved, yet the peak wavelength varies less than 1%, and the peak height and peak width vary less than 6% over broad viewing angles (0--90$^\circ$) under a directional illumination. Our model system consists of dipole scatterers arranged into several rings; interference among the scattered waves is optimized to yield the wavelength-selective and angle-insensitive response. Such designs can be useful for the recently proposed transparent displays that are based on wavelength-selective scattering.

preprint2015arXiv

Power--bandwidth limitations of an optical resonance

We present shape-independent upper limits to the power--bandwidth product for a single resonance in an optical scatterer, with the bound depending only on the material susceptibility. We show that quasistatic metallic scatterers can nearly reach the limits, and we apply our approach to the problem of designing $N$ independent, subwavelength scatterers to achieve flat, broadband response even if they individually exhibit narrow resonant peaks.

preprint2015arXiv

Radiative heat transfer in nonlinear Kerr media

We obtain a fluctuation--dissipation theorem describing thermal electromagnetic fluctuation effects in nonlinear media that we exploit in conjunction with a stochastic Langevin framework to study thermal radiation from Kerr ($χ^{(3)}$) photonic cavities coupled to external environments at and out of equilibrium. We show that that in addition to thermal broadening due to two-photon absorption,the emissivity of such cavities can exhibit asymmetric,non-Lorentzian lineshapes due to self-phase modulation. When the local temperature of the cavity is larger than that of the external bath, we find that the heat transfer into the bath exceeds the radiation from a corresponding linear black body at the same local temperature. We predict that these temperature-tunable thermal processes can be observed in practical, nanophotonic cavities operating at relatively small temperatures.

preprint2015arXiv

Saddle-point integration of $C_\infty$ "bump" functions

This technical note describes the application of saddle-point integration to the asymptotic Fourier analysis of the well-known $C_\infty$ "bump" function $\exp[-(1-x^2)^{-1}]$, deriving both the asymptotic decay rate $k^{-3/4} \exp(-\sqrt k)$ of the Fourier transform $F(k)$ and the exact coefficient. The result is checked against brute-force numerical integration and is extended to generalizations of this bump function.

preprint2015arXiv

Shape-Independent Limits to Near-Field Radiative Heat Transfer

We derive shape-independent limits to the spectral radiative heat-transfer rate between two closely spaced bodies, generalizing the concept of a black body to the case of near-field energy transfer. Through conservation of energy and reciprocity, we show that each body of susceptibility $χ$ can emit and absorb radiation at enhanced rates bounded by $|χ|^2 / \textrm{Im} χ$, optimally mediated by near-field photon transfer proportional to $1/d^2$ across a separation distance $d$. Dipole--dipole and dipole--plate structures approach restricted versions of the limit, but common large-area structures do not exhibit the material enhancement factor and thus fall short of the general limit. By contrast, we find that particle arrays interacting in an idealized Born approximation (i.e., neglecting multiple scattering) exhibit both enhancement factors, suggesting the possibility of orders-of-magnitude improvement beyond previous designs and the potential for radiative heat transfer to be comparable to conductive heat transfer through air at room temperature, and significantly greater at higher temperatures.

preprint2015arXiv

Three-dimensional topological photonic crystal with a single surface Dirac cone

A single Dirac cone on the surface is the hallmark of three-dimensional (3D) topological insulators, where the double degeneracy at the Dirac point is protected by time-reversal symmetry and the spin-splitting away from the point is provided by the spin-orbital coupling. Here we predict a single Dirac-cone surface state in a 3D photonic crystal, where the degeneracy at the Dirac point is protected by a nonsymmorphic glide reflection and the linear splitting away from it is enabled by breaking time-reversal symmetry. Such a gapless surface state is fully robust against random disorder of any type. This bosonic topological band structure is achieved by applying alternating magnetization to gap out the 3D "generalized Dirac points" discovered in the bulk of our crystal. The $Z_2$ bulk invariant is characterized through the evolution of Wannier centers. Our proposal--readily realizable using ferrimagnetic materials at microwave frequencies--can also be regarded as the photonic analog of topological crystalline insulators, providing the first 3D bosonic symmetry-protected topological system.

preprint2014arXiv

An improved beam waist formula for ultrashort, tightly-focused linearly, radially, and azimuthally polarized laser pulses in free space

We derive an asymptotically accurate formula for the beam waist of ultrashort, tightly-focused fundamental linearly-polarized, radially-polarized, and azimuthally-polarized modes in free space. We compute the exact beam waist via numerical cubature to ascertain the accuracy with which our formula approximates the exact beam waist over a broad range of parameters of practical interest. Based on this, we describe a method of choosing parameters in the model given the beam waist and pulse duration of a laser pulse.

preprint2014arXiv

Classical and fluctuation-induced electromagnetic interactions in micronscale systems: designer bonding, antibonding, and Casimir forces

Whether intentionally introduced to exert control over particles and macroscopic objects, such as for trapping or cooling, or whether arising from the quantum and thermal fluctuations of charges in otherwise neutral bodies, leading to unwanted stiction between nearby mechanical parts, electromagnetic interactions play a fundamental role in many naturally occurring processes and technologies. In this review, we survey recent progress in the understanding and experimental observation of optomechanical and quantum-fluctuation forces. Although both of these effects arise from exchange of electromagnetic momentum, their dramatically different origins, involving either real or virtual photons, lead to different physical manifestations and design principles. Specifically, we describe recent predictions and measurements of attractive and repulsive optomechanical forces, based on the bonding and antibonding interactions of evanescent waves, as well as predictions of modified and even repulsive Casimir forces between nanostructured bodies. Finally, we discuss the potential impact and interplay of these forces in emerging experimental regimes of micromechanical devices.

preprint2014arXiv

Generalized Taylor-Duffy Method for Efficient Evaluation of Galerkin Integrals in Boundary-Element Method Computations

We present a generic technique, automated by computer-algebra systems and available as open-source software \cite{scuff-em}, for efficient numerical evaluation of a large family of singular and nonsingular 4-dimensional integrals over triangle-product domains, such as those arising in the boundary-element method (BEM) of computational electromagnetism. To date, practical implementation of BEM solvers has often required the aggregation of multiple disparate integral-evaluation schemes to treat all of the distinct types of integrals needed for a given BEM formulation; in contrast, our technique allows many different types of integrals to be handled by the \emph{same} algorithm and the same code implementation. Our method is a significant generalization of the Taylor--Duffy approach \cite{Taylor2003,Duffy1982}, which was originally presented for just a single type of integrand; in addition to generalizing this technique to a broad class of integrands, we also achieve a significant improvement in its efficiency by showing how the \emph{dimension} of the final numerical integral may often be reduced by one. In particular, if $n$ is the number of common vertices between the two triangles, in many cases we can reduce the dimension of the integral from $4-n$ to $3-n$, obtaining a closed-form analytical result for $n=3$ (the common-triangle case).

preprint2014arXiv

Robust topology optimization of three-dimensional photonic-crystal band-gap structures

We perform full 3D topology optimization (in which "every voxel" of the unit cell is a degree of freedom) of photonic-crystal structures in order to find optimal omnidirectional band gaps for various symmetry groups, including fcc (including diamond), bcc, and simple-cubic lattices. Even without imposing the constraints of any fabrication process, the resulting optimal gaps are only slightly larger than previous hand designs, suggesting that current photonic crystals are nearly optimal in this respect. However, optimization can discover new structures, e.g. a new fcc structure with the same symmetry but slightly larger gap than the well known inverse opal, which may offer new degrees of freedom to future fabrication technologies. Furthermore, our band-gap optimization is an illustration of a computational approach to 3D dispersion engineering which is applicable to many other problems in optics, based on a novel semidefinite-program formulation for nonconvex eigenvalue optimization combined with other techniques such as a simple approach to impose symmetry constraints. We also demonstrate a technique for \emph{robust} topology optimization, in which some uncertainty is included in each voxel and we optimize the worst-case gap, and we show that the resulting band gaps have increased robustness to systematic fabrication errors.

preprint2014arXiv

Speed-of-light limitations in passive linear media

We prove that well-known speed of light restrictions on electromagnetic energy velocity can be extended to a new level of generality, encompassing even nonlocal chiral media in periodic geometries, while at the same time weakening the underlying assumptions to only passivity and linearity of the medium (either with a transparency window or with dissipation). As was also shown by other authors under more limiting assumptions, passivity alone is sufficient to guarantee causality and positivity of the energy density (with no thermodynamic assumptions). Our proof is general enough to include a very broad range of material properties, including anisotropy, bianisotropy (chirality), nonlocality, dispersion, periodicity, and even delta functions or similar generalized functions. We also show that the "dynamical energy density" used by some previous authors in dissipative media reduces to the standard Brillouin formula for dispersive energy density in a transparency window. The results in this paper are proved by exploiting deep results from linear-response theory, harmonic analysis, and functional analysis that had previously not been brought together in the context of electrodynamics.

preprint2014arXiv

The effectiveness of thin films in lieu of hyperbolic metamaterials in the near field

We show that the near-field functionality of hyperbolic metamaterials (HMM), typically proposed for increasing the photonic local density of states (LDOS), can be achieved with thin metal films. Although HMMs have an infinite density of internally-propagating plane-wave states, the external coupling to nearby emitters is severely restricted. We show analytically that properly designed thin films, of thicknesses comparable to the metal size of a hyperbolic metamaterial, yield a LDOS as high as (if not higher than) that of HMMs. We illustrate these ideas by performing exact numerical computations of the LDOS of multilayer HMMs, along with their application to the problem of maximizing near-field heat transfer, to show that thin films are suitable replacements in both cases.

preprint2013arXiv

Comment on "A self-assembled three-dimensional cloak in the visible" in Scientific Reports 3, 2328

Mühlig et. al. propose and fabricate a "cloak" comprised of nano-particles on the surface of a sub-wavelength silica sphere. However, the coating only reduces the scattered fields. This is achieved by increased absorption, such that total extinction increases at all wavelengths. An object creating a large shadow is generally not considered to be cloaked; functionally, in contrast to the relatively few structures that can reduce total extinction, there are many that can reduce scattering alone.

preprint2013arXiv

Efficient Computation of Power, Force, and Torque in BEM Scattering Calculations

We present concise, computationally efficient formulas for several quantities of interest -- including absorbed and scattered power, optical force (radiation pressure), and torque -- in scattering calculations performed using the boundary-element method (BEM) [also known as the method of moments (MOM)]. Our formulas compute the quantities of interest \textit{directly} from the BEM surface currents with no need ever to compute the scattered electromagnetic fields. We derive our new formulas and demonstrate their effectiveness by computing power, force, and torque in a number of example geometries. Free, open-source software implementations of our formulas are available for download online.

preprint2013arXiv

Electromagnetic Wave Source Conditions

This chapter discusses the relationships between current sources and the resulting electromagnetic waves in FDTD simulations. First, the "total-field/scattered-field" approach to creating incident plane waves is reviewed and seen to be a special case of the well-known principle of equivalence in electromagnetism: this can be used to construct "equivalent" current sources for any desired incident field, including waveguide modes. The effects of dispersion and discretization are discussed, and a simple technique to separate incident and scattered fields is described in order to compensate for imperfect equivalent currents. The important concept of the local density of states (LDOS) is reviewed, which elucidates the relationship between current sources and the resulting fields, including enhancement of the LDOS via mode cutoffs (Van Hove singularities) and resonant cavities (Purcell enhancement). We also address various other source techniques such as covering a wide range of frequencies and incident angles in a small number of simulations for waves incident on a periodic surface, sources to excite eigenmodes in rectangular supercells of periodic systems, moving sources, and thermal sources via a Monte Carlo/Langevin approach.

preprint2013arXiv

Fluctuating surface-current formulation of radiative heat transfer: theory and applications

We describe a novel fluctuating-surface current formulation of radiative heat transfer between bodies of arbitrary shape that exploits efficient and sophisticated techniques from the surface-integral-equation formulation of classical electromagnetic scattering. Unlike previous approaches to non-equilibrium fluctuations that involve scattering matrices---relating "incoming" and "outgoing" waves from each body---our approach is formulated in terms of "unknown" surface currents, laying at the surfaces of the bodies, that need not satisfy any wave equation. We show that our formulation can be applied as a spectral method to obtain fast-converging semi-analytical formulas in high-symmetry geometries using specialized spectral bases that conform to the surfaces of the bodies (e.g. Fourier series for planar bodies or spherical harmonics for spherical bodies), and can also be employed as a numerical method by exploiting the generality of surface meshes/grids to obtain results in more complicated geometries (e.g. interleaved bodies as well as bodies with sharp corners). In particular, our formalism allows direct application of the boundary-element method, a robust and powerful numerical implementation of the surface-integral formulation of classical electromagnetism, which we use to obtain results in new geometries, including the heat transfer between finite slabs, cylinders, and cones.

preprint2013arXiv

Fundamental limits to nanoparticle extinction

We show that there are shape-independent upper bounds to the extinction cross section per unit volume of randomly oriented nanoparticles, given only material permittivity. Underlying the limits are restrictive sum rules that constrain the distribution of quasistatic eigenvalues. Surprisingly, optimally-designed spheroids, with only a single quasistatic degree of freedom, reach the upper bounds for four permittivity values. Away from these permittivities, we demonstrate computationally-optimized structures that surpass spheroids and approach the fundamental limits.

preprint2013arXiv

Geometry-induced Casimir suspension of oblate bodies in fluids

We predict that a low-permittivity oblate body (disk-shaped object) above a thin metal substrate (plate with a hole) immersed in a fluidof intermediate permittivity will experience a meta-stable equilibrium (restoring force) near the center of the hole. Stability is the result of a geometry-induced transition in the sign of the force, from repulsive to attractive, that occurs as the disk approaches the hole---in planar or nearly-planar geometries, the same material combination yields a repulsive force at all separations in accordance with the Dzyaloshinskii--Lifshitz--Pitaevskii condition of fluid-induced repulsion between planar bodies. We explore the stability of the system with respect to rotations and lateral translations of the disks, and demonstrate interesting transitions (bifurcations) in the rotational stability of the disks as a function of their size. Finally, we consider the reciprocal situation in which the disk--plate materials are interchanged, and find that in this case the system also exhibits meta-stability. The forces in the system are sufficiently large to be observed in experiments and should enable measurements based on the diffusion dynamics of the suspended bodies.

preprint2013arXiv

High-efficiency degenerate four wave-mixing in triply resonant nanobeam cavities

We demonstrate high-efficiency, degenerate four-wave mixing in triply resonant Kerr $χ^(3)$ photonic crystal (PhC) nanobeam cavities. Using a combination of temporal coupled mode theory and nonlinear finite-difference time-domain (FDTD) simulations, we study the nonlinear dynamics of resonant four-wave mixing processes and demonstrate the possibility of observing high-efficiency limit cycles and steady-state conversion corresponding to $\approx 100$% depletion of the pump light at low powers, even including effects due to losses, self- and cross-phase modulation, and imperfect frequency matching. Assuming operation in the telecom range, we predict close to perfect quantum efficiencies at reasonably low $\sim$ 50 mW input powers in silicon micrometer-scale cavities.

preprint2013arXiv

Optical "Bernoulli" forces

By Bernoulli's law, an increase in the relative speed of a fluid around a body is accompanies by a decrease in the pressure. Therefore, a rotating body in a fluid stream experiences a force perpendicular to the motion of the fluid because of the unequal relative speed of the fluid across its surface. It is well known that light has a constant speed irrespective of the relative motion. Does a rotating body immersed in a stream of photons experience a Bernoulli-like force? We show that, indeed, a rotating dielectric cylinder experiences such a lateral force from an electromagnetic wave. In fact, the sign of the lateral force is the same as that of the fluid-mechanical analogue as long as the electric susceptibility is positive (ε>ε_{0}), but for negative-susceptibility materials (e.g. metals) we show that the lateral force is in the opposite direction. Because these results are derived from a classical electromagnetic scattering problem, Mie-resonance enhancements that occur in other scattering phenomena also enhance the lateral force.

preprint2013arXiv

Transformation Inverse Design

We present a new technique for the design of transformation-optics devices based on large-scale optimization to achieve the optimal effective isotropic dielectric materials within prescribed index bounds, which is computationally cheap because transformation optics circumvents the need to solve Maxwell's equations at each step. We apply this technique to the design of multimode waveguide bends (realized experimentally in a previous paper) and mode squeezers, in which all modes are transported equally without scattering. In addition to the optimization, a key point is the identification of the correct boundary conditions to ensure reflectionless coupling to untransformed regions while allowing maximum flexibility in the optimization. Many previous authors in transformation optics used a certain kind of quasiconformal map which overconstrained the problem by requiring that the entire boundary shape be specified a priori while at the same time underconstraining the problem by employing "slipping" boundary conditions that permit unwanted interface reflections.

preprint2012arXiv

A diameter--bandwidth product limitation of isolated-object cloaking

We show that cloaking of isolated objects is subject to a diameter-bandwidth product limitation: as the size of the object increases, the bandwidth of good (small cross-section) cloaking decreases inversely with the diameter, as a consequence of causality constraints even for perfect fabrication and materials with negligible absorption. This generalizes a previous result that perfect cloaking of isolated objects over a nonzero bandwidth violates causality. Furthermore, we demonstrate broader causality-based scaling limitations on any bandwidth-averaged cloaking cross-section, using complex analysis and the optical theorem to transform the frequency-averaged problem into a single scattering problem with transformed materials.

preprint2012arXiv

Anomalous near-field heat transfer between a cylinder and a perforated surface

We predict that the radiative heat-transfer rate between a cylinder and a perforated surface depends non-monotonically on their separation. This anomalous behavior, which arises due to near-field effects, is explained using a heuristic model based on the interaction of a dipole with a plate. We show that nonmonotonicity depends not only on geometry and temperature but also on material dispersion - for micron and submicron objects, nonmonotonicity is present in polar dielectrics but absent in metals with small skin depths.

preprint2012arXiv

Fluctuating Surface Currents: A New Algorithm for Efficient Prediction of Casimir Interactions among Arbitrary Materials in Arbitrary Geometries. I. Theory

This paper presents a new method for the efficient numerical computation of Casimir interactions between objects of arbitrary geometries, composed of materials with arbitrary frequency-dependent electrical properties. Our method formulates the Casimir effect as an interaction between effective electric and magnetic current distributions on the surfaces of material bodies, and obtains Casimir energies, forces, and torques from the spectral properties of a matrix that quantifies the interactions of these surface currents. The method can be formulated and understood in two distinct ways: \textbf{(1)} as a consequence of the familiar \textit{stress-tensor} approach to Casimir physics, or, alternatively, \textbf{(2)} as a particular case of the \textit{path-integral} approach to Casimir physics, and we present both formulations in full detail. In addition to providing an algorithm for computing Casimir interactions in geometries that could not be efficiently handled by any other method, the framework proposed here thus achieves an explicit unification of two seemingly disparate approaches to computational Casimir physics.

preprint2012arXiv

Fluctuating surface-current formulation of radiative heat transfer for arbitrary geometries

We describe a fluctuating surface-current formulation of radiative heat transfer, applicable to arbitrary geometries, that directly exploits standard, efficient, and sophisticated techniques from the boundary-element method. We validate as well as extend previous results for spheres and cylinders, and also compute the heat transfer in a more complicated geometry consisting of two interlocked rings. Finally, we demonstrate that the method can be readily adapted to compute the spatial distribution of heat flux on the surface of the interacting bodies.

preprint2012arXiv

Fluctuation-Induced Phenomena in Nanoscale Systems: Harnessing the Power of Noise

The famous Johnson-Nyquist formula relating noise current to conductance has a microscopic generalization relating noise current density to microscopic conductivity, with corollary relations governing noise in the components of the electromagnetic fields. These relations, known collectively in physics as fluctuation-dissipation relations, form the basis of the modern understanding of fluctuation-induced phenomena, a field of burgeoning importance in experimental physics and nanotechnology. In this review, we survey recent progress in computational techniques for modeling fluctuation-induced phenomena, focusing on two cases of particular interest: near-field radiative heat transfer and Casimir forces. In each case we review the basic physics of the phenomenon, discuss semi-analytical and numerical algorithms for theoretical analysis, and present recent predictions for novel phenomena in complex material and geometric configurations.

preprint2012arXiv

Optomechanical and Photothermal Interactions in Suspended Photonic Crystal Membranes

We present here an optomechanical system fabricated with novel stress management techniques that allow us to suspend an ultrathin defect-free silicon photonic-crystal membrane above a Silicon-on-Insulator (SOI) substrate with a gap that is tunable to below 200 nm. Our devices are able to generate strong attractive and repulsive optical forces over a large surface area with simple in- and outcoupling and feature the strongest repulsive optomechanical coupling in any geometry to date (gOM/2π ~ -65 GHz/nm). The interplay between the optomechanical and photo-thermal-mechanical dynamics is explored, and the latter is used to achieve cooling and amplification of the mechanical mode, demonstrating that our platform is well-suited for applications in low-power mass, force, and refractive index sensing as well as and optomechanical accelerometry.

preprint2011arXiv

Bonding, antibonding and tunable optical forces in asymmetric membranes

We demonstrate that tunable attractive (bonding) and repulsive (anti-bonding) forces can arise in highly asymmetric structures coupled to external radiation, a consequence of the bonding/anti-bonding level repulsion of guided-wave resonances that was first predicted in symmetric systems. Our focus is a geometry consisting of a photonic-crystal (holey) membrane suspended above an unpatterned layered substrate, supporting planar waveguide modes that can couple via the periodic modulation of the holey membrane. Asymmetric geometries have a clear advantage in ease of fabrication and experimental characterization compared to symmetric double-membrane structures. We show that the asymmetry can also lead to unusual behavior in the force magnitudes of a bonding/antibonding pair as the membrane separation changes, including nonmonotonic dependences on the separation. We propose a computational method that obtains the entire force spectrum via a single time-domain simulation, by Fourier-transforming the response to a short pulse and thereby obtaining the frequency-dependent stress tensor. We point out that by operating with two, instead of a single frequency, these evanescent forces can be exploited to tune the spring constant of the membrane without changing its equilibrium separation.

preprint2011arXiv

Casimir micro-sphere diclusters and three-body effects in fluids

Our previous article [Phys. Rev. Lett. 104, 060401 (2010)] predicted that Casimir forces induced by the material-dispersion properties of certain dielectrics can give rise to stable configurations of objects. This phenomenon was illustrated via a dicluster configuration of non-touching objects consisting of two spheres immersed in a fluid and suspended against gravity above a plate. Here, we examine these predictions from the perspective of a practical experiment and consider the influence of non-additive, three-body, and nonzero-temperature effects on the stability of the two spheres. We conclude that the presence of Brownian motion reduces the set of experimentally realizable silicon/teflon spherical diclusters to those consisting of layered micro-spheres, such as the hollow- core (spherical shells) considered here.

preprint2011arXiv

Casimir repulsion beyond the dipole regime

We extend a previous result [Phys. Rev. Lett. 105, 090403 (2010)] on Casimir repulsion between a plate with a hole and a cylinder centered above it to geometries in which the central object can no longer be treated as a point dipole. We show through numerical calculations that as the distance between the plate and central object decreases, there is an intermediate regime in which the repulsive force increases dramatically. Beyond this, the force rapidly switches over to attraction as the separation decreases further to zero, in line with the proximity force approximation. We demonstrate that this effect can be understood as a competition between an increased repulsion due to a larger polarizability of the central object interacting with increased fringing fields near the edge of the plate, and attractive forces due primarily to the nonzero thickness of the plate. In comparison with our previous work, we find that using the same plate geometry but replacing the single cylinder with a ring of cylinders, or more generally an extended uniaxial conductor, the repulsive force can be enhanced by a factor of approximately $10^3$. We conclude that this enhancement, although quite dramatic, is still too small to yield detectable repulsive Casimir forces.

preprint2011arXiv

Computation of Casimir Interactions between Arbitrary 3D Objects with Arbitrary Material Properties

We extend a recently introduced method for computing Casimir forces between arbitrarily--shaped metallic objects [M. T. H. Reid et al., Phys. Rev. Lett._103_ 040401 (2009)] to allow treatment of objects with arbitrary material properties, including imperfect conductors, dielectrics, and magnetic materials. Our original method considered electric currents on the surfaces of the interacting objects; the extended method considers both electric and magnetic surface current distributions, and obtains the Casimir energy of a configuration of objects in terms of the interactions of these effective surface currents. Using this new technique, we present the first predictions of Casimir interactions in several experimentally relevant geometries that would be difficult to treat with any existing method. In particular, we investigate Casimir interactions between dielectric nanodisks embedded in a dielectric fluid; we identify the threshold surface--surface separation at which finite--size effects become relevant, and we map the rotational energy landscape of bound nanoparticle diclusters.

preprint2011arXiv

Design of thin-film photonic metamaterial Lüneburg lens using analytical approach

We design an all-dielectric Lüneburg lens as an adiabatic space-variant lattice explicitly accounting for finite film thickness. We describe an all-analytical approach to compensate for the finite height of subwavelength dielectric structures in the pass-band regime. This method calculates the effective refractive index of the infinite-height lattice from effective medium theory, then embeds a medium of the same effective index into a slab waveguide of finite height and uses the waveguide dispersion diagram to calculate a new effective index. The results are compared with the conventional numerical treatment - a direct band diagram calculation, using a modified three-dimensional lattice with the superstrate and substrate included in the cell geometry. We show that the analytical results are in good agreement with the numerical ones, and the performance of the thin-film Lüneburg lens is quite different than the estimates obtained assuming infinite height.

preprint2011arXiv

Designing evanescent optical interactions to control the expression of Casimir forces in optomechanical structures

We propose an optomechanical structure consisting of a photonic-crystal (holey) membrane suspended above a layered silicon-on-insulator substrate in which resonant bonding/antibonding optical forces created by externally incident light from above enable all-optical control and actuation of stiction effects induced by the Casimir force. In this way, one can control how the Casimir force is expressed in the mechanical dynamics of the membrane, not by changing the Casimir force directly but by optically modifying the geometry and counteracting the mechanical spring constant to bring the system in or out of regimes where Casimir physics dominate. The same optical response (reflection spectrum) of the membrane to the incident light can be exploited to accurately measure the effects of the Casimir force on the equilibrium separation of the membrane.

preprint2011arXiv

Frequency-selective near-field enhancement of radiative heat transfer via photonic-crystal slabs: a general computational approach for arbitrary geometries and materials

We demonstrate the possibility of achieving enhanced frequency-selective near-field radiative heat transfer between patterned (photonic crystal) slabs at designable frequencies and separations, exploiting a general numerical approach for computing heat transfer in arbitrary geometries and materials based on the finite-difference time-domain method. Our simulations reveal a tradeoff between selectivity and near-field enhancement as the slab--slab separation decreases, with the patterned heat transfer eventually reducing to the unpatterned result multiplied by a fill factor (described by a standard proximity approximation). We also find that heat transfer can be further enhanced at selective frequencies when the slabs are brought into a glide-symmetric configuration, a consequence of the degeneracies associated with the non-symmorphic symmetry group.

preprint2011arXiv

General scaling limitations of ground-plane and isolated-object cloaks

We prove that, for arbitrary three-dimensional transformation-based invisibility cloaking of an object above a ground plane or of isolated object, there are practical constraints that increase with the object size. In particular, we show that the cloak thickness must scale proportional to the thickness of the object being cloaked, assuming bounded refractive indices, and that absorption discrepancies and other imperfections must scale inversely with the object thickness. For isolated objects, we also show that bounded refractive indices imply a lower bound on the effective cross-section.

preprint2011arXiv

Modeling near-field radiative heat transfer from sharp objects using a general 3d numerical scattering technique

We examine the non-equilibrium radiative heat transfer between a plate and finite cylinders and cones, making the first accurate theoretical predictions for the total heat transfer and the spatial heat flux profile for three-dimensional compact objects including corners or tips. We find qualitatively different scaling laws for conical shapes at small separations, and in contrast to a flat/slightly-curved object, a sharp cone exhibits a local \emph{minimum} in the spatially resolved heat flux directly below the tip. The method we develop, in which a scattering-theory formulation of thermal transfer is combined with a boundary-element method for computing scattering matrices, can be applied to three-dimensional objects of arbitrary shape.

preprint2011arXiv

Structural anisotropy and orientation-induced Casimir repulsion in fluids

In this work we theoretically consider the Casimir force between two periodic arrays of nanowires (both in vacuum, and on a substrate separated by a fluid) at separations comparable to the period. Specifically, we compute the dependence of the exact Casimir force between the arrays under both lateral translations and rotations. Although typically the force between such structures is well-characterized by the Proximity Force Approximation (PFA), we find that in the present case the microstructure modulates the force in a way qualitatively inconsistent with PFA. We find instead that effective-medium theory, in which the slabs are treated as homogeneous, anisotropic dielectrics, gives a surprisingly accurate picture of the force, down to separations of half the period. This includes a situation for identical, fluid-separated slabs in which the exact force changes sign with the orientation of the wire arrays, whereas PFA predicts attraction. We discuss the possibility of detecting these effects in experiments, concluding that this effect is strong enough to make detection possible in the near future.

preprint2010arXiv

A novel boundary element method using surface conductive absorbers for full-wave analysis of 3-D nanophotonics

Fast surface integral equation (SIE) solvers seem to be ideal approaches for simulating 3-D nanophotonic devices, as these devices generate fields both in an interior channel and in the infinite exterior domain. However, many devices of interest, such as optical couplers, have channels that can not be terminated without generating reflections. Generating absorbers for these channels is a new problem for SIE methods, as the methods were initially developed for problems with finite surfaces. In this paper we show that the obvious approach for eliminating reflections, making the channel mildly conductive outside the domain of interest, is inaccurate. We describe a new method, in which the absorber has a gradually increasing surface conductivity; such an absorber can be easily incorporated in fast integral equation solvers. Numerical experiments from a surface-conductivity modified FFT-accelerated PMCHW-based solver are correlated with analytic results, demonstrating that this new method is orders of magnitude more effective than a volume absorber, and that the smoothness of the surface conductivity function determines the performance of the absorber. In particular, we show that the magnitude of the transition reflection is proportional to 1/L^(2d+2), where L is the absorber length and d is the order of the differentiability of the surface conductivity function.

preprint2010arXiv

Achieving a Strongly Temperature-Dependent Casimir Effect

We propose a method of achieving large temperature sensitivity in the Casimir force that involves measuring the stable separation between dielectric objects immersed in fluid. We study the Casimir force between slabs and spheres using realistic material models, and find large > 2nm/K variations in their stable separations (hundreds of nanometers) near room temperature. In addition, we analyze the effects of Brownian motion on suspended objects, and show that the average separation is also sensitive to changes in temperature . Finally, this approach also leads to rich qualitative phenomena, such as irreversible transitions, from suspension to stiction, as the temperature is varied.

preprint2010arXiv

Calculation of nonzero-temperature Casimir forces in the time domain

We show how to compute Casimir forces at nonzero temperatures with time-domain electromagnetic simulations, for example using a finite-difference time-domain (FDTD) method. Compared to our previous zero-temperature time-domain method, only a small modification is required, but we explain that some care is required to properly capture the zero-frequency contribution. We validate the method against analytical and numerical frequency-domain calculations, and show a surprising high-temperature disappearance of a non-monotonic behavior previously demonstrated in a piston-like geometry.

preprint2010arXiv

Casimir repulsion between metallic objects in vacuum

We give an example of a geometry in which two metallic objects in vacuum experience a repulsive Casimir force. The geometry consists of an elongated metal particle centered above a metal plate with a hole. We prove that this geometry has a repulsive regime using a symmetry argument and confirm it with numerical calculations for both perfect and realistic metals. The system does not support stable levitation, as the particle is unstable to displacements away from the symmetry axis.

preprint2010arXiv

Degenerate four-wave mixing in triply-resonant Kerr cavities

We demonstrate theoretical conditions for highly-efficient degenerate four-wave mixing in triply-resonant nonlinear (Kerr) cavities. We employ a general and accurate temporal coupled-mode analysis in which the interaction of light in arbitrary microcavities is expressed in terms a set of coupling coefficients that we rigorously derive from the full Maxwell equations. Using the coupled-mode theory, we show that light consisting of an input signal of frequency $ω_0-Δω$ can, in the presence of pump light at $ω_0$, be converted with quantum-limited efficiency into an output shifted signal of frequency $ω_0 + Δω$, and we derive expressions for the critical input powers at which this occurs. We find that critical powers in the order of 10mW assuming very conservative cavity parameters (modal volumes $\sim10$ cubic wavelengths and quality factors $\sim1000$. The standard Manley-Rowe efficiency limits are obtained from the solution of the classical coupled-mode equations, although we also derive them from simple photon-counting "quantum" arguments. Finally, using a linear stability analysis, we demonstrate that maximal conversion efficiency can be retained even in the presence of self- and cross-phase modulation effects that generally act to disrupt the resonance condition.

preprint2010arXiv

Delay-bandwidth and delay-loss limitations for cloaking of large objects

Based on a simple model of ground-plane cloaking, we argue that the diffculty of cloaking is fundamentally limited by delay-loss and delaylbandwidth/size limitations that worsen as the size of the object to be cloaked increases relative to the wavelength. These considerations must be taken into account when scaling experimental cloaking demonstrations from wavelength-scale objects towards larger sizes, and suggest quantitative material/loss challenges in cloaking human-scale objects.

preprint2010arXiv

Is the electrostatic force between a point charge and a neutral metallic object always attractive?

We give an example of a geometry in which the electrostatic force between a point charge and a neutral metallic object is repulsive. The example consists of a point charge centered above a thin metallic hemisphere, positioned concave up. We show that this geometry has a repulsive regime using both a simple analytical argument and an exact calculation for an analogous two-dimensional geometry. Analogues of this geometry-induced repulsion can appear in many other contexts, including Casimir systems.

preprint2010arXiv

Microstructure Effects for Casimir Forces in Chiral Metamaterials

We examine a recent prediction for the chirality-dependence of the Casimir force in chiral metamaterials by numerical computation of the forces between the exact microstructures, rather than homogeneous approximations. We compute the exact force for a chiral bent-cross pattern, as well as forces for an idealized "omega"-particle medium in the dilute approximation and identify the effects of structural inhomogeneity (i.e. proximity forces and anisotropy). We find that these microstructure effects dominate the force for separations where chirality was predicted to have a strong influence. To get observations of chirality free from microstructure effects, one must go to large separations where the effect of chirality is at most $\sim10^{-4}$ of the total force.

preprint2010arXiv

Numerical methods for computing Casimir interactions

We review several different approaches for computing Casimir forces and related fluctuation-induced interactions between bodies of arbitrary shapes and materials. The relationships between this problem and well known computational techniques from classical electromagnetism are emphasized. We also review the basic principles of standard computational methods, categorizing them according to three criteria---choice of problem, basis, and solution technique---that can be used to classify proposals for the Casimir problem as well. In this way, mature classical methods can be exploited to model Casimir physics, with a few important modifications.

preprint2010arXiv

Sufficient conditions for two-dimensional localization by arbitrarily weak defects in periodic potentials with band gaps

We prove, via an elementary variational method, 1d and 2d localization within the band gaps of a periodic Schrodinger operator for any mostly negative or mostly positive defect potential, V, whose depth is not too great compared to the size of the gap. In a similar way, we also prove sufficient conditions for 1d and 2d localization below the ground state of such an operator. Furthermore, we extend our results to 1d and 2d localization in d dimensions; for example, a linear or planar defect in a 3d crystal. For the case of D-fold degenerate band edges, we also give sufficient conditions for localization of up to D states.

preprint2009arXiv

Casimir forces in the time domain II: Applications

Our preceding paper introduced a method to compute Casimir forces in arbitrary geometries and for arbitrary materials that was based on a finite-difference time-domain (FDTD) scheme. In this manuscript, we focus on the efficient implementation of our method for geometries of practical interest and extend our previous proof-of-concept algorithm in one dimension to problems in two and three dimensions, introducing a number of new optimizations. We consider Casimir piston-like problems with nonmonotonic and monotonic force dependence on sidewall separation, both for previously solved geometries to validate our method and also for new geometries involving magnetic sidewalls and/or cylindrical pistons. We include realistic dielectric materials to calculate the force between suspended silicon waveguides or on a suspended membrane with periodic grooves, also demonstrating the application of PML absorbing boundaries and/or periodic boundaries. In addition we apply this method to a realizable three-dimensional system in which a silica sphere is stably suspended in a fluid above an indented metallic substrate. More generally, the method allows off-the-shelf FDTD software, already supporting a wide variety of materials (including dielectric, magnetic, and even anisotropic materials) and boundary conditions, to be exploited for the Casimir problem.

preprint2009arXiv

Difference-frequency generation with quantum-limited efficiency in triply-resonant nonlinear cavities

We present a comprehensive study of second-order nonlinear difference frequency generation in triply resonant cavities using a theoretical framework based on coupled-mode theory. We show that optimal quantum-limited conversion efficiency can be achieved at any pump power when the powers at the pump and idler frequencies satisfy a critical relationship. We demonstrate the existence of a broad parameter range in which all triply-resonant DFG processes exhibit monostable conversion. We also demonstrate the existence of a geometry-dependent bistable region.

preprint2009arXiv

Efficient Computation of Casimir Interactions between Arbitrary 3D Objects

We introduce an efficient technique for computing Casimir energies and forces between objects of arbitrarily complex 3D geometries. In contrast to other recently developed methods, our technique easily handles non-spheroidal, non-axisymmetric objects and objects with sharp corners. Using our new technique, we obtain the first predictions of Casimir interactions in a number of experimentally relevant geometries, including crossed cylinders and tetrahedral nanoparticles.

preprint2009arXiv

Efficient Terahertz Generation in Triply Resonant Nonlinear Photonic Crystal Microcavities

We propose a scheme for efficient cavity-enhanced nonlinear THz generation via difference-frequency generation (DFG) processes using a triply resonant system based on photonic crystal cavities. We show that high nonlinear overlap can be achieved by coupling a THz cavity to a doubly-resonant, dual-polarization near-infrared (e.g. telecom band) photonic-crystal nanobeam cavity, allowing the mixing of three mutually orthogonal fundamental cavity modes through a chi(2) nonlinearity. We demonstrate through coupled-mode theory that complete depletion of the pump frequency - i.e., quantum-limited conversion - is possible in an experimentally feasible geometry, with the operating output power at the point of optimal total conversion efficiency adjustable by varying the mode quality (Q) factors.

preprint2009arXiv

Ingredients of a Casimir analog computer

We present the basic ingredients of a technique to compute quantum Casimir forces at micrometer scales using antenna measurements at tabletop, e.g. centimeter, scales, forming a type of analog computer for the Casimir force. This technique relies on a correspondence that we derive between the contour integration of the Casimir force in the complex frequency plane and the electromagnetic response of a physical dissipative medium in a finite, real frequency bandwidth.

preprint2009arXiv

Non-touching nanoparticle diclusters bound by repulsive and attractive Casimir forces

We present a scheme for obtaining stable Casimir suspension of dielectric nontouching objects immersed in a fluid, validated here in various geometries consisting of ethanol-separated dielectric spheres and semi-infinite slabs. Stability is induced by the dispersion properties of real dielectric (monolithic) materials. A consequence of this effect is the possibility of stable configurations (clusters) of compact objects, which we illustrate via a "molecular" two-sphere dicluster geometry consiting of two bound spheres levitated above a gold slab. Our calculations also reveal a strong interplay between material and geometric dispersion, and this is exemplified by the qualitatively different stability behavior observed in planar versus spherical geometries.

preprint2009arXiv

Spherical cloaking using nonlinear transformations for improved segmentation into concentric isotropic coatings

Two novel classes of spherical invisibility cloaks based on nonlinear transformation have been studied. The cloaking characteristics are presented by segmenting the nonlinear transformation based spherical cloak into concentric isotropic homogeneous coatings. Detailed investigations of the optimal discretization (e.g., thickness control of each layer, nonlinear factor, etc.) are presented for both linear and nonlinear spherical cloaks and their effects on invisibility performance are also discussed. The cloaking properties and our choice of optimal segmentation are verified by the numerical simulation of not only near-field electric-field distribution but also the far-field radar cross section (RCS).

preprint2008arXiv

Nonlinear harmonic generation and devices in doubly-resonant Kerr cavities

We describea theoretical analysis of the nonlinear dynamics of third-harmonic generation ($ω\to3ω$) via Kerr ($\chithree$) nonlinearities in a resonant cavity with resonances at both $ω$ and $3ω$. Such a doubly resonant cavity greatly reduces the required power for efficient harmonic generation, by a factor of $\sim V/Q^2$ where $V$ is the modal volume and $Q$ is the lifetime, and can even exhibit 100% harmonic conversion efficiency at a critical input power. However, we show that it also exhibits a rich variety of nonlinear dynamics, such as multistable solutions and long-period limit cycles.We describe how to compensate for self/cross-phase modulation (which otherwise shifts the cavity frequencies out of resonance), and how to excite the different stable solutions (and especially the high-efficiency solutions) by specially modulated input pulses.

preprint2008arXiv

Repulsive, nonmonotonic Casimir forces in a glide-symmetric geometry

We describe a three-dimensional geometry that exhibits a repulsive Casimir force using ordinary metallic materials, as computed via an exact numerical method (no uncontrolled approximations). The geometry consists of a zippelike, glide-symmetric structure formed of interleaved metal brackets attached to parallel plates. Depending on the separation, the perpendicular force between the plates/brackets varies from attractive (large separations) to repulsive (intermediate distances) and back to attractive (close separations), with one point of stable equilibrium in the perpendicular direction. This geometry was motivated by a simple intuition of attractive interactions between surfaces, and so we also consider how a rough proximity force approximation of pairwise attractions compares to the exact calculations.

preprint2007arXiv

Chi2 and chi3 harmonic generation at a critical power in inhomogeneous doubly resonant cavities

We derive general conditions for 100 percent frequency conversion in any doubly resonant nonlinear cavity, for both second- and third-harmonic generation via chi2 and chi3 nonlinearities. We find that conversion efficiency is optimized for a certain critical power depending on the cavity parameters, and assuming reasonable parameters we predict 100 percent conversion using milliwatts of power or less. These results follow from a semi-analytical coupled-mode theory framework which is generalized from previous work to include both chi2 and chi3 media as well as inhomogeneous (fully vectorial) cavities, analyzed in the high-efficiency limit where down-conversion processes lead to a maximum efficiency at the critical power, and which is verified by direct finite-difference time-domain (FDTD) simulations of the nonlinear Maxwell equations. Explicit formulas for the nonlinear coupling coefficients are derived in terms of the linear cavity eigenmodes, which can be used to design and evaluate cavities in arbitrary geometries.

preprint2007arXiv

Computation and visualization of Casimir forces in arbitrary geometries: non-monotonic lateral forces and failure of proximity-force approximations

We present a method of computing Casimir forces for arbitrary geometries, with any desired accuracy, that can directly exploit the efficiency of standard numerical-electromagnetism techniques. Using the simplest possible finite-difference implementation of this approach, we obtain both agreement with past results for cylinder-plate geometries, and also present results for new geometries. In particular, we examine a piston-like problem involving two dielectric and metallic squares sliding between two metallic walls, in two and three dimensions, respectively, and demonstrate non-additive and non-monotonic changes in the force due to these lateral walls.

preprint2007arXiv

Optical-approximation analysis of sidewall-spacing effects on the force between two squares with parallel sidewalls

Using the ray-optics approximation, we analyze the Casimir force in a two dimensional domain formed by two metallic blocks adjacent to parallel metallic sidewalls, which are separated from the blocks by a finite distance h. For h > 0, the ray-optics approach is not exact because diffraction effects are neglected. Nevertheless, we show that ray optics is able to qualitatively reproduce a surprising effect recently identified in an exact numerical calculation: the force between the blocks varies non-monotonically with h. In this sense, the ray-optics approach captures an essential part of the physics of multi-body interactions in this system, unlike simpler pairwise-interaction approximations such as PFA. Furthermore, by comparison to the exact numerical results, we are able to quantify the impact of diffraction on Casimir forces in this geometry.

Steven G. Johnson

What is connected

Connect this record

See the researcher in context

Building this map preview

86 published item(s)

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

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

Free-electron-light interactions in nanophotonics

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

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

A general framework for scintillation in nanophotonics

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

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

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

Physics-informed neural networks with hard constraints for inverse design

Towards 3D-Printed Inverse-Designed Metaoptics

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

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

Inverse designed metalenses with extended depth of focus

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

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

Response to Comment on Widely tunable compact terahertz gas lasers

Inverse Design of Nanoparticles for Enhanced Raman Scattering

Overlapping domains for topology optimization of large-area metasurfaces

Fundamental limits to optical response in absorptive systems

Interaction-induced mode switching in steady-state microlasers

Perfect Single-Sided Radiation and Absorption without Mirrors

Ab-initio multimode linewidth theory for arbitrary inhomogeneous laser cavities

All-optical three-dimensional electron pulse compression

Cavity-enhanced second harmonic generation via nonlinear-overlap optimization

Design of diamond microcavities for single photon frequency down-conversion

Fluctuating volume-current formulation of electromagnetic fluctuations in inhomogeneous media: incandecence and luminescence in arbitrary geometries

Optimization of sharp and viewing-angle-independent structural color

Power--bandwidth limitations of an optical resonance

Radiative heat transfer in nonlinear Kerr media

Saddle-point integration of $C_\infty$ "bump" functions

Shape-Independent Limits to Near-Field Radiative Heat Transfer

Three-dimensional topological photonic crystal with a single surface Dirac cone

An improved beam waist formula for ultrashort, tightly-focused linearly, radially, and azimuthally polarized laser pulses in free space

Classical and fluctuation-induced electromagnetic interactions in micronscale systems: designer bonding, antibonding, and Casimir forces

Generalized Taylor-Duffy Method for Efficient Evaluation of Galerkin Integrals in Boundary-Element Method Computations

Robust topology optimization of three-dimensional photonic-crystal band-gap structures

Speed-of-light limitations in passive linear media

The effectiveness of thin films in lieu of hyperbolic metamaterials in the near field

Comment on "A self-assembled three-dimensional cloak in the visible" in Scientific Reports 3, 2328

Efficient Computation of Power, Force, and Torque in BEM Scattering Calculations

Electromagnetic Wave Source Conditions

Fluctuating surface-current formulation of radiative heat transfer: theory and applications

Fundamental limits to nanoparticle extinction

Geometry-induced Casimir suspension of oblate bodies in fluids

High-efficiency degenerate four wave-mixing in triply resonant nanobeam cavities

Optical "Bernoulli" forces

Transformation Inverse Design

A diameter--bandwidth product limitation of isolated-object cloaking

Anomalous near-field heat transfer between a cylinder and a perforated surface

Fluctuating Surface Currents: A New Algorithm for Efficient Prediction of Casimir Interactions among Arbitrary Materials in Arbitrary Geometries. I. Theory

Fluctuating surface-current formulation of radiative heat transfer for arbitrary geometries

Fluctuation-Induced Phenomena in Nanoscale Systems: Harnessing the Power of Noise

Optomechanical and Photothermal Interactions in Suspended Photonic Crystal Membranes

Bonding, antibonding and tunable optical forces in asymmetric membranes

Casimir micro-sphere diclusters and three-body effects in fluids

Casimir repulsion beyond the dipole regime

Computation of Casimir Interactions between Arbitrary 3D Objects with Arbitrary Material Properties

Design of thin-film photonic metamaterial Lüneburg lens using analytical approach

Designing evanescent optical interactions to control the expression of Casimir forces in optomechanical structures

Frequency-selective near-field enhancement of radiative heat transfer via photonic-crystal slabs: a general computational approach for arbitrary geometries and materials

General scaling limitations of ground-plane and isolated-object cloaks

Modeling near-field radiative heat transfer from sharp objects using a general 3d numerical scattering technique

Structural anisotropy and orientation-induced Casimir repulsion in fluids

A novel boundary element method using surface conductive absorbers for full-wave analysis of 3-D nanophotonics

Achieving a Strongly Temperature-Dependent Casimir Effect

Calculation of nonzero-temperature Casimir forces in the time domain

Casimir repulsion between metallic objects in vacuum

Degenerate four-wave mixing in triply-resonant Kerr cavities

Delay-bandwidth and delay-loss limitations for cloaking of large objects

Is the electrostatic force between a point charge and a neutral metallic object always attractive?

Microstructure Effects for Casimir Forces in Chiral Metamaterials

Numerical methods for computing Casimir interactions

Sufficient conditions for two-dimensional localization by arbitrarily weak defects in periodic potentials with band gaps