-
Factorized Machine Learning for Performance Modeling of Massively Parallel Heterogeneous Physical Simulations
Authors:
Ardavan Oskooi,
Christopher Hogan,
Alec M. Hammond,
M. T. Homer Reid,
Steven G. Johnson
Abstract:
We demonstrate neural-network runtime prediction for complex, many-parameter, massively parallel, heterogeneous-physics simulations running on cloud-based MPI clusters. Because individual simulations are so expensive, it is crucial to train the network on a limited dataset despite the potentially large input space of the physics at each point in the spatial domain. We achieve this using a two-part…
▽ More
We demonstrate neural-network runtime prediction for complex, many-parameter, massively parallel, heterogeneous-physics simulations running on cloud-based MPI clusters. Because individual simulations are so expensive, it is crucial to train the network on a limited dataset despite the potentially large input space of the physics at each point in the spatial domain. We achieve this using a two-part strategy. First, we perform data-driven static load balancing using regression coefficients extracted from small simulations, which both improves parallel performance and reduces the dependency of the runtime on the precise spatial layout of the heterogeneous physics. Second, we divide the execution time of these load-balanced simulations into computation and communication, factoring crude asymptotic scalings out of each term, and training neural nets for the remaining factor coefficients. This strategy is implemented for Meep, a popular and complex open-source electrodynamics simulation package, and are validated for heterogeneous simulations drawn from published engineering models.
△ Less
Submitted 6 October, 2020; v1 submitted 9 March, 2020;
originally announced March 2020.
-
Taylor-Duffy Method for Singular Tetrahedron-Product Integrals: Efficient Evaluation of Galerkin Integrals for VIE Solvers
Authors:
M. T. Homer Reid
Abstract:
I present an accurate and efficient technique for numerical evaluation of singular 6-dimensional integrals over tetrahedon-product domains, with applications to calculation of Galerkin matrix elements for discretized volume-integral-equation (VIE) solvers using Schaubert-Wilton-Glisson (SWG) and other tetrahedral basis functions. My method extends the generalized Taylor-Duffy strategy---used to ha…
▽ More
I present an accurate and efficient technique for numerical evaluation of singular 6-dimensional integrals over tetrahedon-product domains, with applications to calculation of Galerkin matrix elements for discretized volume-integral-equation (VIE) solvers using Schaubert-Wilton-Glisson (SWG) and other tetrahedral basis functions. My method extends the generalized Taylor-Duffy strategy---used to handle the singular \textit{triangle}-product integrals arising in discretized surface-integral-equation (SIE) formulations---to the tetrahedron-product case; it effects an exact transformation of a singular 6-dimensional integral to an nonsingular lower-dimensional integral that may be evaluated by simple numerical cubature The method is highly general and may---with the aid of automatic code generation facilitated by computer-algebra systems---be applied to a wide variety of singular integrals arising in various VIE formulations with various types of tetrahedral basis function, of which I present several examples. To demonstrate the accuracy and efficiency of my method, I apply it to the calculation of matrix elements for the volume electric-field integral equation (VEFIE) discretized with SWG basis functions, where the method yields 12-digit or higher accuracy with low computational cost---an improvement of many orders of magnitude compared to existing techniques.
△ Less
Submitted 7 August, 2017;
originally announced August 2017.
-
Photon Torpedoes and Rytov Pinwheels: Integral-Equation Modeling of Non-Equilibrium Fluctuation-Induced Forces and Torques on Nanoparticles
Authors:
M. T. Homer Reid,
O. D. Miller,
A. G. Polimeridis,
A. W. Rodriguez,
E. M. Tomlinson,
S. G. Johnson
Abstract:
We present new theoretical tools, based on fluctuational electrodynamics and the integral-equation approach to computational electromagnetism, for numerical modeling of forces and torques on bodies of complex shapes and materials due to emission of thermal radiation out of thermal equilibrium. This extends our recently-developed fluctuating-surface-current (FSC) and fluctuating-volume-current (FVC…
▽ More
We present new theoretical tools, based on fluctuational electrodynamics and the integral-equation approach to computational electromagnetism, for numerical modeling of forces and torques on bodies of complex shapes and materials due to emission of thermal radiation out of thermal equilibrium. This extends our recently-developed fluctuating-surface-current (FSC) and fluctuating-volume-current (FVC) techniques for radiative heat transfer to the computation of non-equilibrium fluctuation-induced forces and torques; as we show, the extension is non-trivial due to the greater computational cost of modeling radiative momentum transfer, including new singularities that must be carefully neutralized. We introduce a new analytical cancellation technique that addresses these challenges and allows, for the first time, accurate and efficient prediction of non-equilibrium forces and torques on bodies of essentially arbitrary shapes---including asymmetric and chiral particles---and complex material properties, including continuously-varying and anisotropic dielectrics. We validate our approach by showing that it reproduces known results, then present new numerical predictions of non-equilibrium self-propulsion, self-rotation, and momentum-transfer phenomena in complex geometries that would be difficult or impossible to study with existing methods. Our findings indicate that the fluctuation-induced dynamics of micron-size room-temperature bodies in cold environments involve microscopic length scales but macroscopic time scales, with typical linear and angular velocities on the order of microns/second and radians/second; For a micron-scale gear driven by thermal radiation from a nearby chiral emitter, we find a strong and non-monotonic dependence of the magnitude and even the \textit{sign} of the induced torque on the temperature of the emitter.
△ Less
Submitted 7 August, 2017;
originally announced August 2017.
-
Limits to the Optical Response of Graphene and 2D Materials
Authors:
Owen D. Miller,
Ognjen Ilic,
Thomas Christensen,
M. T. Homer Reid,
Harry A. Atwater,
John D. Joannopoulos,
Marin Soljacic,
Steven G. Johnson
Abstract:
2D materials provide a platform for strong light--matter interactions, creating wide-ranging design opportunities via new-material discoveries and new methods for geometrical structuring. We derive general upper bounds to the strength of such light--matter interactions, given only the optical conductivity of the material, including spatial nonlocality, and otherwise independent of shape and config…
▽ More
2D materials provide a platform for strong light--matter interactions, creating wide-ranging design opportunities via new-material discoveries and new methods for geometrical structuring. We derive general upper bounds to the strength of such light--matter interactions, given only the optical conductivity of the material, including spatial nonlocality, and otherwise independent of shape and configuration. Our material figure of merit shows that highly doped graphene is an optimal material at infrared frequencies, whereas single-atomic-layer silver is optimal in the visible. For quantities ranging from absorption and scattering to near-field spontaneous-emission enhancements and radiative heat transfer, we consider canonical geometrical structures and show that in certain cases the bounds can be approached, while in others there may be significant opportunity for design improvement. The bounds can encourage systematic improvements in the design of ultrathin broadband absorbers, 2D antennas, and near-field energy harvesters.
△ Less
Submitted 2 August, 2017; v1 submitted 9 May, 2017;
originally announced May 2017.
-
Non-intuitive Computational Optimization of Illumination Patterns for Maximum Optical Force and Torque
Authors:
Yoonkyung E. Lee,
Owen D. Miller,
M. T. Homer Reid,
Steven G. Johnson,
Nicholas X. Fang
Abstract:
This paper aims to maximize optical force and torque on arbitrary micro- and nano-scale objects using numerically optimized structured illumination. By developing a numerical framework for computer-automated design of 3d vector-field illumination, we demonstrate a 20-fold enhancement in optical torque per intensity over circularly polarized plane wave on a model plasmonic particle. The nonconvex o…
▽ More
This paper aims to maximize optical force and torque on arbitrary micro- and nano-scale objects using numerically optimized structured illumination. By developing a numerical framework for computer-automated design of 3d vector-field illumination, we demonstrate a 20-fold enhancement in optical torque per intensity over circularly polarized plane wave on a model plasmonic particle. The nonconvex optimization is efficiently performed by combining a compact cylindrical Bessel basis representation with a fast boundary element method and a standard derivative-free, local optimization algorithm. We analyze the optimization results for 2000 random initial configurations, discuss the tradeoff between robustness and enhancement, and compare the different effects of multipolar plasmon resonances on enhancing force and torque. All results are obtained using open-source computational software available online.
△ Less
Submitted 26 January, 2017;
originally announced January 2017.
-
Fundamental limits to optical response in absorptive systems
Authors:
Owen D. Miller,
Athanasios G. Polimeridis,
M. T. Homer Reid,
Chia Wei Hsu,
Brendan G. DeLacy,
John D. Joannopoulos,
Marin Soljačić,
Steven G. Johnson
Abstract:
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 pe…
▽ More
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.
△ Less
Submitted 6 January, 2016; v1 submitted 12 March, 2015;
originally announced March 2015.
-
On the Computation of Power in Volume Integral Equation Formulations
Authors:
A. G. Polimeridis,
M. T. H. Reid,
S. G. Johnson,
J. K. White,
A. W. Rodriguez
Abstract:
We present simple and stable formulas for computing power (including absorbed/radiated, scattered and extinction power) in current-based volume integral equation formulations. The proposed formulas are given in terms of vector-matrix-vector products of quantities found solely in the associated linear system. In addition to their efficiency, the derived expressions can guarantee the positivity of t…
▽ More
We present simple and stable formulas for computing power (including absorbed/radiated, scattered and extinction power) in current-based volume integral equation formulations. The proposed formulas are given in terms of vector-matrix-vector products of quantities found solely in the associated linear system. In addition to their efficiency, the derived expressions can guarantee the positivity of the computed power. We also discuss the application of Poynting's theorem for the case of sources immersed in dissipative materials. The formulas are validated against results obtained both with analytical and numerical methods for scattering and radiation benchmark cases.
△ Less
Submitted 12 September, 2014; v1 submitted 27 June, 2014;
originally announced June 2014.
-
Transformation Optics scheme for two-dimensional materials
Authors:
Anshuman Kumar,
Kin Hung Fung,
M. T. Homer Reid,
Nicholas Fang
Abstract:
Two dimensional optical materials, such as graphene can be characterized by a surface conductivity. So far, the transformation optics schemes have focused on three dimensional properties such as permittivity $ε$ and permeability $μ$. In this paper, we use a scheme for transforming surface currents to highlight that the surface conductivity transforms in a way different from $ε$ and $μ$. We use thi…
▽ More
Two dimensional optical materials, such as graphene can be characterized by a surface conductivity. So far, the transformation optics schemes have focused on three dimensional properties such as permittivity $ε$ and permeability $μ$. In this paper, we use a scheme for transforming surface currents to highlight that the surface conductivity transforms in a way different from $ε$ and $μ$. We use this surface conductivity transformation to demonstrate an example problem of reducing scattering of plasmon mode from sharp protrusions in graphene.
△ Less
Submitted 3 March, 2014;
originally announced March 2014.
-
Generalized Taylor-Duffy Method for Efficient Evaluation of Galerkin Integrals in Boundary-Element Method Computations
Authors:
M. T. Homer Reid,
Steven G. Johnson,
Jacob K. White
Abstract:
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 solv…
▽ More
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).
△ Less
Submitted 19 July, 2014; v1 submitted 5 December, 2013;
originally announced December 2013.
-
Fundamental limits to nanoparticle extinction
Authors:
Owen D. Miller,
Chia W. Hsu,
M. T. Homer Reid,
Wenjun Qiu,
Brendan G. DeLacy,
John D. Joannopoulos,
Marin Soljačić,
Steven G. Johnson
Abstract:
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…
▽ More
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.
△ Less
Submitted 27 November, 2013; v1 submitted 6 November, 2013;
originally announced November 2013.
-
Geometry-induced Casimir suspension of oblate bodies in fluids
Authors:
Alejandro W. Rodriguez,
M. T. Homer Reid,
Francesco Intravaia,
Alexander Woolf,
Diego A. R. Dalvit,
Federico Capasso,
Steven G. Johnson
Abstract:
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 approache…
▽ More
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.
△ Less
Submitted 16 October, 2013;
originally announced October 2013.
-
Photon Emission Rate Engineering using Graphene Nanodisc Cavities
Authors:
Anshuman Kumar,
Kin Hung Fung,
M. T. Homer Reid,
Nicholas X. Fang
Abstract:
In this work, we present a systematic study of the plasmon modes in a system of vertically stacked pair of graphene discs. Quasistatic approximation is used to model the eigenmodes of the system. Eigen-response theory is employed to explain the spatial dependence of the coupling between the plasmon modes and a quantum emitter. These results show a good match between the semi-analytical calculation…
▽ More
In this work, we present a systematic study of the plasmon modes in a system of vertically stacked pair of graphene discs. Quasistatic approximation is used to model the eigenmodes of the system. Eigen-response theory is employed to explain the spatial dependence of the coupling between the plasmon modes and a quantum emitter. These results show a good match between the semi-analytical calculation and full-wave simulations. Secondly, we have shown that it is possible to engineer the decay rates of a quantum emitter placed inside and near this cavity, using Fermi level tuning, via gate voltages and variation of emitter location and polarization. We highlighted that by coupling to the bright plasmon mode, the radiative efficiency of the emitter can be enhanced compared to the single graphene disc case, whereas the dark plasmon mode suppresses the radiative efficiency.
△ Less
Submitted 3 March, 2014; v1 submitted 9 September, 2013;
originally announced September 2013.
-
Efficient Computation of Power, Force, and Torque in BEM Scattering Calculations
Authors:
M. T. Homer Reid,
Steven G. Johnson
Abstract:
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…
▽ More
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.
△ Less
Submitted 10 July, 2013;
originally announced July 2013.
-
Nanoplasmonic near-field synthesis
Authors:
Johannes Feist,
M. T. Homer Reid,
Matthias F. Kling
Abstract:
The temporal response of resonances in nanoplasmonic structures typically converts an incoming few-cycle field into a much longer near-field at the spot where non-linear physical phenomena including electron emission, recollision and high-harmonic generation can take place. We show that for practically useful structures pulse shaping of the incoming pulse can be used to synthesize the plasmon-enha…
▽ More
The temporal response of resonances in nanoplasmonic structures typically converts an incoming few-cycle field into a much longer near-field at the spot where non-linear physical phenomena including electron emission, recollision and high-harmonic generation can take place. We show that for practically useful structures pulse shaping of the incoming pulse can be used to synthesize the plasmon-enhanced field and enable single-cycle driven nonlinear physical phenomena. Our method is demonstrated for the generation of an isolated attosecond pulse by plasmon-enhanced high harmonic generation. We furthermore show that optimal control techniques can be used even if the response of the plasmonic structure is not known a priori.
△ Less
Submitted 15 March, 2013; v1 submitted 29 August, 2012;
originally announced August 2012.
-
Fluctuation-Induced Phenomena in Nanoscale Systems: Harnessing the Power of Noise
Authors:
M. T. Homer Reid,
Alejandro W. Rodriguez,
Steven G. Johnson
Abstract:
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-i…
▽ More
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.
△ Less
Submitted 19 July, 2012; v1 submitted 17 July, 2012;
originally announced July 2012.
-
Fluctuating surface-current formulation of radiative heat transfer for arbitrary geometries
Authors:
Alejandro W. Rodriguez,
M. T. Homer Reid,
Steven G. Johnson
Abstract:
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…
▽ More
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.
△ Less
Submitted 3 July, 2012; v1 submitted 6 June, 2012;
originally announced June 2012.
-
Modeling near-field radiative heat transfer from sharp objects using a general 3d numerical scattering technique
Authors:
Alexander P. McCauley,
M. T. Homer Reid,
Matthias Krüger,
Steven G. Johnson
Abstract:
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…
▽ More
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.
△ Less
Submitted 17 July, 2011; v1 submitted 11 July, 2011;
originally announced July 2011.