Limitations on bandwidth-integrated passive cloaking
Benjamin Strekha
Department of Electrical and Computer Engineering, Princeton University, Princeton, New Jersey 08544, USA
Alessio Amaolo
Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
Jewel Mohajan
Department of Electrical and Computer Engineering, Princeton University, Princeton, New Jersey 08544, USA
Pengning Chao
Department of Mathematics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
Sean Molesky
Department of Engineering Physics, Polytechnique Montréal, Montréal, Québec H3T 1J4, Canada
Alejandro W. Rodriguez
Department of Electrical and Computer Engineering, Princeton University, Princeton, New Jersey 08544, USA
Abstract
We present a general framework for the computation of structure-agnostic bounds on the performance of passive cloaks over a nonzero bandwidth.
We apply this framework in 2D to the canonical scenario of cloaking a circular object.
We find that perfect cloaking using a finite-sized isotropic cloak is impossible over any bandwidth, with the bounds scaling linearly with the bandwidth before saturating due to the finite size of the cloak and the presence of material loss.
The bounds also exhibit linear scaling with material loss in the cloak and linear scaling with the inverse of the radial thickness of the design region before saturation due to finite-size effects or the presence of material loss.
The formulation could readily find applications in the development of cloaking devices, setting expectations and benchmarks for optimal performance.
Introduction.— The ability to engineer optical cloaking of objects to make them undetectable by observers has wide-ranging implications for radar [1], obscurance [2, 3], and sensing [4] applications.
Several techniques have been proposed to achieve (near-)invisibility at a single operating optical wavelength, including transformation optics [5], scattering-cancellation via plasmonic [6] or mantle cloaking [7], waveguide cloaking [8], transmission-line networks [9], or the use of anomalous localized resonances [10], to name a few.
A fundamental and practically relevant question is whether one can achieve (near-)invisibility over a wide spectral window, e.g., the entire visible spectrum.
As pointed out by Pendry et al. using an argument based on the phase and group velocities of light in the presence of dispersion [5], and more recently by others using more formal arguments for the inevitable distortion of a pulse wave by materials with local response [11], perfect cloaking of an isolated object in vacuum over a nonzero bandwidth is impossible due to causality; the finite speed of light and the presence of dispersion limits fully effective cloaking to a single frequency.
However, the Pendry et al. and Miller proofs say little about imperfect cloaking (a nonzero but small scattering cross section).
For instance, engineering a cloak to attain near-perfect invisibility at a single frequency does not preclude “good” cloaking over some nonzero bandwidth around that frequency.
Hashemi et al. [12] presented an alternative proof that agrees with the results of Pendry et al. and Miller for the impossibility of perfect cloaking over a nonzero bandwidth for physical, causal materials, but additionally derived scaling relations with implications for design, e.g., the allowed bandwidth over which one may expect practical cloaking scales inversely with the diameter of the cloaked object.
In this article, we adapt recently developed constrained optimization approaches [13, 14, 15] to study the following question
related to bounds on cloaking performance: Given a specified cloak material along with an operating
frequency and bandwidth, what is the best possible performance of a passive cloaking device that can be formed out of that
material within a prescribed region surrounding a cloaked object?
In addition to capturing scaling behaviors, the presented formalism allows for quantitative assessments of expected performance of cloaks. Typical realizations of cloaks require anisotropic metamaterials with extreme and exotic properties, such as permittivities and permeabilities less than one and even equal to zero in some regions of the cloak [5, 16, 17].
This raises the natural question of what is the best achievable performance using the simplest materials possible.
While the numerical bounds lack the transparency and intuition of analytic expressions, they can still provide insight into scaling properties and general trends, and even lead to fairly tight limits (within an order of magnitude of structures discovered via topology optimization).
Formulation.—
In line with the majority of previous works on cloaking systems [12, 18, 19, 15], and as a proof of concept, we study the common scenario of a cloaked object illuminated by an incident plane wave emitted by a far-away source. We consider extinction, a nonnegative quantity for passive systems [20], as the fundamental figure of merit, for the following reasons. First, it is better suited for describing the practical goals of cloaks as opposed to scattered or absorbed power alone:
minimizing extinction minimizes the scattered field over both near- (absorbed power) and far-field (scattered power) domains [15], implying that an object becomes truly undetectable regardless of the position of an observer as the figure of merit approaches zero.
Second, analyticity in the upper-half of the complex frequency plane means that one can exploit the residue theorem to evaluate the spectral average over a nonzero bandwidth via a single scattering calculation [21, 12], which greatly simplifies computations.
Throughout, we use the language of scattering theory laid out in previous descriptions of electromagnetic bounds [14, 22, 23], wherein the operator describing bound polarization currents in the medium is defined by the relation
(1)
The operator represents the background Green’s function, which for
vacuum satisfies . (Note that in keeping with prior work, this definition includes an additional factor of compared to the standard convention [24]).
The operator is the scattering potential (susceptibility) relative to this background medium (whatever additional material response was not included in the definition of ), and and are defined as the incident electric fields and electric currents in the background, respectively.
When acting on incident fields, the operator produces the generated current ,
with denoting the impedance of free space.
Employing the operator relation to a scattering problem involving a prescribed incident field [14], the extinguished power may be written as,
(2)
where is left inside Im because it will be continued to the complex plane to map this quantity to a bandwidth-averaged quantity [21, 12]. Namely, since real sources emit light over a nonzero bandwidth a key figure of merit is the bandwidth average of extinct power.
Considering a Lorentzian window function centered at with a bandwidth , the average extinction power for incident plane waves,
(3)
may be evaluated by closing the contour and picking up the residue at the pole of the Lorentzian in the upper-half complex frequency plane, yielding Eq. (2) with all vectors and operators evaluated at a complex frequency and complex wave number , simplifying the calculation so that one no longer needs to consider power at each individual frequency within the Lorentzian window.
Note also that working at a complex frequency requires the evaluation of , which is mathematically equivalent to using modified materials at a real frequency [12].
To isolate the impact of the cloak on scattering from a fixed object, we let refer to the vacuum Green’s function and shift the known scattering properties of the fixed object into the background of the scattering operator describing the cloak, with denoting the Green’s function of the fixed object in isolation, so that satisfies .
In particular, let be the scattering operator of the fixed object (the object to be cloaked) in isolation, and let be the scattering operator of the designable object (the cloak) dressed by the fixed object. After some straightforward algebraic manipulations (see Appendix), we find
(4)
The first term is the extinguished power when no cloak is present while only the second term depends on the cloaked object and quantifies the interaction between the cloak and cloaked object.
Our derivation of bounds exploits the optimization procedure based on Lagrange duality laid out in Refs. [25, 13].
The loosest such bound only imposes that the optimal scattering operator satisfies the
conservation of power (optical theorem [20])
over the entire design domain, and not the full scattering equations.
Defining , the total field in the presence of only the fixed object, and so that is the induced current in the cloak, we relax the problem such that, instead of optimizing over with support in the design region, we optimize over with support in the design region, i.e., optimize over all possible polarization currents within the design region so that the bound considered here automatically takes into account all possible distributions of vacuum and the prescribed cloak material within the design region.
It should be noted, however, that this bound formalism does not yield the optimal material distribution of the cloak [13], a task that calls for more computationally demanding topology optimization methods (NP-hard problems that preclude guarantees of globally optimal solutions) [26].
Similar techniques have recently been used to derive bounds on deterministic scattering and fluctuational electrodynamic phenomena [14, 27, 23, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37].
Concretely, we wish to solve the following problem:
(5a)
(5b)
where is a projection operator into the design region , is a projection operator into a subregion contained within the design region, and .
Taking the real and imaginary parts of the complex-valued constraints enforce the conservation of resistive and reactive power, respectively, spatially integrated over a region [13].
This is a quadratically constrained quadratic program over the field and we numerically compute bounds on Eq. (5a), the primal objective, by evaluating and optimizing the corresponding concave Lagrange dual function [25].
Figure 1:
(a) Bounds on bandwidth-integrated extinction power as a function of the bandwidth.
All curves and markers refer to performance bounds and inverse designs, respectively, expressed as a ratio of the extinguished power from incident plane waves in the presence and absence of a cloak for and (bottom-right schematic).
Insets: Representative inverse designs corresponding to .
(b) Performance of structures discovered via topology optimization for , , , , and (solid) or (dashed).
Top: The cloak suppresses extinction within a frequency range around at the expense of worse performance in other parts of the spectrum.
Bottom: Real part of the total out-of-plane electric field at . Green circles outline the region of space where the object and cloak are contained, with an inner and outer radius and , respectively.
Applications.— We now exploit the framework above to numerically obtain performance bounds for passive cloaks composed of an isotropic electric susceptibility. For simplicity, we focus on 2D settings and consider the scenario of a fixed circular object of radius , a common benchmark example for cloaking [38], with the cloak restricted to an annular design region with an inner radius and outer radius [see Fig. 1(a) bottom-right schematic] for the case of incident TM (electric field out of plane) plane waves.
First, we consider the situation of illumination of broadband light from the far-field (plane waves). For a fixed we find that the bounds scale linearly with the bandwidth before saturating due to the finite size of the design domain and due to the presence of material loss, see Fig. 1(a).
As seen in the figure, enforcing the conservation of power within the design domain leads to nontrivial bounds with interesting trends and scaling behavior that are also fairly tight (within an order of magnitude of the objective values of discovered structures).
The structures found via topology optimization perform well within a bandwidth around the “center” frequency , but at the expense of more scattering in nearby regions of the spectrum, see Fig. 1(b).
This is in agreement with the findings in Ref. [39], which showed that any passive cloak made of a linear, nondiamagnetic material respecting causality always increases the scattering and extinction integrated over all frequencies compared to the original uncloaked object.
Thus, near-invisibility in a given frequency window necessarily implies significant scattered/extinct power at other frequencies, and therefore such cloaks may be detected more easily than the original uncloaked objects when illuminated by sufficiently broadband light (short pulses).
As seen in Fig. 1(b), the electric field in the lossy metallic object is expelled from and rerouted around the object with the introduction of the cloak. The dielectric object with a smaller loss, on the other hand,
may try to create zero-field regions or it may allow the electric field to penetrate the object if it is advantageous as is the case for the dielectric cloak example shown in Fig. 1(b).
Next, we investigate the performance limits of 2D cloaks of isotropic susceptibility as a function of the material loss, see Fig. 2.
The bounds (solid lines) follow trends seen in topology-optimized designs (circular markers) over a broad range of material loss values and both demonstrate, unsurprisingly, improved performance as the material loss in the cloak decreases, exhibiting linear scaling with decreasing material loss before eventually saturating. This scaling can be understood as follows.
Writing the total operator in block-form over the two regions, namely,
(6)
(7)
(8)
leads, after the use of the Woodbury formula for a matrix inverse, to a total operator expression of the form
(9)
Here, we introduced which is the scattering operator of the designable cloak dressed by the presence of the fixed object, as seen from and recognizing as the Green’s function of the fixed object in isolation in vacuum.
Likewise for but with the fixed and design regions switching roles.
Letting be the separation of the susceptibility into its real and imaginary parts, we find
(10)
and similar for but with replaced with .
Using these expansions for and in the block-form for the total scattering operator, Eq. (9), leads to the prediction that the extinction at a single frequency, , scales linearly with towards a constant for a fixed geometric structure as .
The numerical Lagrange dual bounds follow the same scaling behavior as the primal objective function for a fixed structure, demonstrating that faster scaling is not possible.
The approach to a positive constant rather than 0 as the material loss of the cloak vanishes is in agreement with the findings of Hashemi et al. [40],
which showed that if the attainable refractive index contrast is bounded from below (as is the case for a cloak of isotropic susceptibility considered here) then there is a bound on the reduction of the scattering cross section for transformation-based invisibility cloaking of an isolated object; in particular, a minimum achievable refractive index contrast greater than 0 necessarily implies a positive cross section.
A Pendry cloak can achieve perfect invisibility at a single frequency with a finite device footprint but it requires a vanishing permittivity and permeability at the inner surface of the cloak [5, 16].
As shown, however, the presented formalism can provide quantitative lower bounds supporting such qualitative observations.
Figure 3 plots limit values as a function of the radial size of the design region, for a given object radius.
Unsurprisingly, the lower bounds are monotonically nonincreasing as the allowed cloak footprint increases since any polarization current with a given performance can be contained in an enlarged design domain.
In the single frequency case, , the bounds scale inversely with the thickness of the design region, , over a notable range of thicknesses.
Hashemi et al. [40] argued that for bounded refractive indices, the cloak thickness must scale proportionally to the thickness of the object being cloaked, and the scaling of our numerical bounds demonstrates this behavior.
For positive bandwidths, , the bounds saturate to a positive value, with saturation occurring rather quickly (around ) for bandwidths .
This has interesting implications for the design of cloaks since it is typically not clear how large a device needs to be to achieve near-optimal or even reasonable performance.
For our chosen parameters, the bounds show that for nonzero bandwidths a cloak design region with a thickness equal to that of the cloaked object essentially saturates the possible performance; larger device footprints will not always lead to significant improvements in performance.
Lastly, the presented inverse designs and bounds were restricted to plane-wave sources coming from a single direction of incidence.
To consider cloaking robust to the direction of incidence, the figure of merit may be modified to encompass multiple plane wave directions (including fully angle-integrated extinction), such that
where each is a plane wave incident from a different angle.
This leads to an optimization problem of the form
(11a)
(11b)
where and .
Without the constraints between the terms, the problem reduces to a sum of decoupled optimization problems.
Although each angle is a different scattering problem,
the polarization currents induced in each scenario are generated by the same structured media and this fact is enforced by the “cross constraints”, leading to additional tightening of the bounds [41, 42].
Since the incident field is the same up to a rotation, if the object and design region are invariant under rotations then the bounds are the same for decoupled problems and any tightening should come from the cross-constraints.
To consider a more complicated scenario involving more than a single incident source, Fig. 4
shows computed bounds and inverse design performance values for a rectangular object discovered starting from random initializations for the case where , comprising horizontally and vertically propagating plane waves.
Figure 2:
Bounds (solid lines) and inverse designs (circles) relating to net extinction power (cloaking performance) for the system of Fig. 1, at a single frequency and as a function of the material loss in the cloak for and , and for both dielectric and metallic objects.
Figure 3:
Bounds on extinction power for the same system as Fig. 1 but as a function of the outer radius of the design region , for a fixed object radius , and for different values of (dashed), (dotted), and zero bandwidth (solid); for comparison, crosses, circles and stars, respectively, are performance values obtained from inverse designs. Figure 4:
Bounds on bandwidth-integrated extinction as a function of the bandwidth with either vertically or horizontally incident plane waves, or an objective function with an equally weighted average of both (bottom-right schematic).
The rectangular object has lengths and while the rectangular design region has lengths .
The susceptibility values are and .
Conclusion.—
In conclusion, we presented a formalism for computing bounds on passive, broadband cloaking systems that can provide benchmarks and set expectations for future cloaking devices, particularly those developed through the use of topology optimization.
The many variables available for tuning means an exhaustive study of their relationship to the bounds on performance is complex.
We presented a few key examples to demonstrate the capabilities of the formalism and to extract intuition for the obvious parameters one may wish to vary in experiments.
Our numerical results support the claims that broadband cloaking of electrically large objects is practically impossible using passive linear cloaks [12, 38, 18] and provide a quantitative assessment of the limits on performance.
We believe this represents an important result for the science of cloaking, suggesting that new concepts and designs, including opening to the field of nonlinear and active metamaterials [38], are necessary for increasing the cloaking bandwidth
and compensating the unwanted scattering required by linearity and causality in the presence of absorption losses [18].
In closing, we note that while our examples focused on 2D and isotropic permittivities, the formalism is valid in 3D settings, and may be easily extended to consider anisotropic permittivities and permeabilities.
Acknowledgments.—
We acknowledge the support of a Princeton SEAS Innovation Grant and by the Cornell Center for Materials Research (MRSEC) through Award DMR-1719875.
S.M. acknowledges financial support from the Canada First Research Excellence Fund via the Institut de Valorisation des Données (IVADO) collaboration.
The simulations presented in this article were performed on computational
resources managed and supported by Princeton Research
Computing, a consortium of groups including the Princeton Institute for Computational Science and Engineering
(PICSciE) and the Office of Information Technology’s
High Performance Computing Center and Visualization Laboratory at Princeton University.
Appendix A Objective function derivation
In this section, we derive an expression for the total (including the fixed and designable objects) scattering operator which isolates the contribution of the designable object.
A useful operator for this purpose is the operator defined by
(12)
From Eqs. (12) and (1), it is clear that and .
The operator takes an incident field and gives the generated current, , while the
operator takes an initial current and gives the total current, .
We will make use of the nesting property of the operator
(13)
where is the operator of the total system (including all the scatterers), is the operator of the fixed object only, and is the operator for the designable body with the fixed object contained in the background, i.e., replaced by (which is the total Green’s function of the fixed object in the background, ) in the defining relation for .
The nesting property holds since .
Using and the nesting relation Eq. (13) yields
(14)
(15)
(16)
(17)
(18)
where we made use of .
Define and .
Use to find that the total extinct power by all the scatterers, relative to a vacuum background, is given by
(19)
(20)
(21)
with subject to the fundamental relation
(22)
where is the susceptibility used for the object in the design domain (the cloak).
This constraint is similar to that in previous works [14, 28], but where the relevant Green’s function is and not the vacuum Green’s function.
Use of leads to the primal objective function presented in the main text.
Appendix B Topology optimization procedure
To calculate inverse designs, we made use of the NLopt package [43] and followed standard topology optimization algorithms [26, 44] based on the method of moving asymptotes [45] where the electric susceptibility value of each pixel in the specified design region is considered as an independent design parameter.
Each pixel in the design region is allowed to explore a continuous range of susceptibility values varying linearly between the vacuum value and the prescribed material value .
The continuous range can be interpreted, after normalization by , as a filling fraction of the pixel.
We implemented a 2D Maxwell solver which could be used for arbitrary structuring (allow the pixels within any arbitrary 2D design region to vary).
Doing a thousand function evaluations (iterations of optimization over the susceptibility profile for fixed material parameters, design region, bandwidth, etc.) in 2D takes less than 1 hour for a wavelength-scale design region.
The primary challenge in such optimization problems lies in the computational cost of evaluating the objective function and gradient, which is evaluated several hundred if not thousand times during the course of any one optimization.
Let denote the total Green’s function and let where .
After the discretization of the computational grid, then in the design region at pixel .
The gradient of the objective function with respect to the degree of freedom is given by
(23)
(24)
(25)
(26)
(27)
(28)
where is a topology optimization degree of freedom, normalized to take values in , so that the susceptibility at the -th pixel in the design region is given by , is projection onto the pixel indexed by , is the total electric field, and which for reciprocal systems becomes .
Hence, for a given incident field one only needs to perform two additional Maxwell solves to calculate and at each iteration of the optimization to obtain the required gradient Eq. (28) and the objective function
.
Details on the computational complexities encountered in 3D and potential solution methods can be found in Appendix C of Ref. [14] and the therein cited references.
References
Nicolai and Carichner [2010]
Leland M. Nicolai and Grant E. Carichner.
Fundamentals of Aircraft and Airship Design: Volume I–Aircraft Design.
American Institute of Aeronautics and Astronautics, Inc., 2010.
Li and Pendry [2008]
Jensen Li and John B. Pendry.
Hiding under the carpet: a new strategy for cloaking.
Physical Review Letters, 101(20):203901, 2008.
Ergin et al. [2010]
Tolga Ergin, Nicolas Stenger, Patrice Brenner, John B. Pendry, and Martin Wegener.
Three-dimensional invisibility cloak at optical wavelengths.
Science, 328(5976):337–339, 2010.
Alù and Engheta [2009]
Andrea Alù and Nader Engheta.
Cloaking a sensor.
Physical Review Letters, 102(23):233901, 2009.
Pendry [2006]
J. B. Pendry.
Controlling electromagnetic fields.
Science, 312(5781):1780–1782, June 2006.
ISSN 0036-8075, 1095-9203.
doi: 10.1126/science.1125907.
Alù and Engheta [2005]
Andrea Alù and Nader Engheta.
Achieving transparency with plasmonic and metamaterial coatings.
Physical Review E, 72(1):016623, 2005.
Alù [2009]
Andrea Alù.
Mantle cloak: Invisibility induced by a surface.
Physical Review B, 80(24):245115, 2009.
Tretyakov et al. [2009]
Sergei Tretyakov, Pekka Alitalo, Olli Luukkonen, and Constantin Simovski.
Broadband electromagnetic cloaking of long cylindrical objects.
Physical Review Letters, 103(10):103905, 2009.
Alitalo et al. [2008]
Pekka Alitalo, Olli Luukkonen, Liisi Jylha, Jukka Venermo, and Sergei A. Tretyakov.
Transmission-line networks cloaking objects from electromagnetic fields.
IEEE Transactions on Antennas and Propagation, 56(2):416–424, 2008.
Milton and Nicorovici [2006]
Graeme W. Milton and Nicolae-Alexandru P. Nicorovici.
On the cloaking effects associated with anomalous localized resonance.
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 462(2074):3027–3059, 2006.
Miller [2006]
David A. B. Miller.
On perfect cloaking.
Optics Express, 14(25):12457–12466, 2006.
Hashemi et al. [2012]
Hila Hashemi, Cheng-Wei Qiu, Alexander P. McCauley, J. D. Joannopoulos, and Steven G. Johnson.
Diameter-bandwidth product limitation of isolated-object cloaking.
Physical Review A, 86(1):013804, July 2012.
ISSN 1050-2947, 1094-1622.
doi: 10.1103/PhysRevA.86.013804.
Chao et al. [2022a]
Pengning Chao, Benjamin Strekha, Rodrick Kuate Defo, Sean Molesky, and Alejandro W. Rodriguez.
Physical limits in electromagnetism.
Nature Reviews Physics, 4(8):543–559, August 2022a.
ISSN 2522-5820.
doi: 10.1038/s42254-022-00468-w.
Molesky et al. [2020a]
Sean Molesky, Pengning Chao, Weiliang Jin, and Alejandro W. Rodriguez.
Global operator bounds on electromagnetic scattering: Upper bounds on far-field cross sections.
Physical Review Research, 2(3):033172, 2020a.
Jelinek et al. [2021]
Lukas Jelinek, Mats Gustafsson, Miloslav Capek, and Kurt Schab.
Fundamental bounds on the performance of monochromatic passive cloaks.
Optics Express, 29(15):24068, July 2021.
ISSN 1094-4087.
doi: 10.1364/OE.428536.
Cai et al. [2007]
Wenshan Cai, Uday K. Chettiar, Alexander V. Kildishev, and Vladimir M. Shalaev.
Optical cloaking with metamaterials.
Nature Photonics, 1(4):224–227, 2007.
Schurig et al. [2006]
David Schurig, John B. Pendry, and David R. Smith.
Calculation of material properties and ray tracing in transformation media.
Optics Express, 14(21):9794–9804, 2006.
Monticone and Alù [2016]
Francesco Monticone and Andrea Alù.
Invisibility exposed: Physical bounds on passive cloaking.
Optica, 3(7):718–724, 2016.
Cassier and Milton [2017]
Maxence Cassier and Graeme W. Milton.
Bounds on Herglotz functions and fundamental limits of broadband passive quasistatic cloaking.
Journal of Mathematical Physics, 58(7):071504, 2017.
Jackson [1999]
John David Jackson.
Classical Electrodynamics.
Wiley, New York, 3rd edition, 1999.
ISBN 978-0-471-30932-1.
Kuang and Miller [2020]
Zeyu Kuang and Owen D. Miller.
Computational bounds to light–matter interactions via local conservation laws.
Physical Review Letters, 125(26):263607, December 2020.
doi: 10.1103/PhysRevLett.125.263607.
Molesky et al. [2020b]
Sean Molesky, Prashanth S. Venkataram, Weiliang Jin, and Alejandro W. Rodriguez.
Fundamental limits to radiative heat transfer: Theory.
Physical Review B, 101(3):035408, 2020b.
Venkataram et al. [2020a]
Prashanth S. Venkataram, Sean Molesky, Weiliang Jin, and Alejandro W. Rodriguez.
Fundamental limits to radiative heat transfer: The limited role of nanostructuring in the near-field.
Physical Review Letters, 124(1):013904, 2020a.
Novotny and Hecht [2012]
Lukas Novotny and Bert Hecht.
Principles of Nano-Optics.
Cambridge University Press, Cambridge, second edition, 2012.
ISBN 978-1-139-55425-1.
Boyd and Vandenberghe [2004]
Stephen P. Boyd and Lieven Vandenberghe.
Convex Optimization.
Cambridge University Press, Cambridge, UK; New York, 2004.
ISBN 978-0-521-83378-3.
Molesky et al. [2018]
Sean Molesky, Zin Lin, Alexander Y. Piggott, Weiliang Jin, Jelena Vučković, and Alejandro W. Rodriguez.
Inverse design in nanophotonics.
Nature Photonics, 12(11):659–670, 2018.
Molesky et al. [2020c]
Sean Molesky, Pengning Chao, and Alejandro W. Rodriguez.
Hierarchical mean-field operator bounds on electromagnetic scattering: Upper bounds on near-field radiative Purcell enhancement.
Physical Review Research, 2(4):043398, December 2020c.
doi: 10.1103/PhysRevResearch.2.043398.
Chao et al. [2022b]
Pengning Chao, Rodrick Kuate Defo, Sean Molesky, and Alejandro Rodriguez.
Maximum electromagnetic local density of states via material structuring.
Nanophotonics, 12(3):549–557, 2022b.
Angeris et al. [2019]
Guillermo Angeris, Jelena Vučković, and Stephen P. Boyd.
Computational bounds for photonic design.
ACS Photonics, 6(5):1232, 2019.
doi: 10.1021/acsphotonics.9b00154.
Shim et al. [2019]
Hyungki Shim, Lingling Fan, Steven G. Johnson, and Owen D. Miller.
Fundamental Limits to Near-Field Optical Response over Any Bandwidth.
Physical Review X, 9(1):011043, March 2019.
ISSN 2160-3308.
doi: 10.1103/PhysRevX.9.011043.
Amaolo et al. [2024a]
Alessio Amaolo, Pengning Chao, Thomas J. Maldonado, Sean Molesky, and Alejandro W. Rodriguez.
Can photonic heterostructures provably outperform single-material geometries?
Nanophotonics, 13(3):283–288, 2024a.
Mohajan et al. [2023]
Jewel Mohajan, Pengning Chao, Weiliang Jin, Sean Molesky, and Alejandro W. Rodriguez.
Fundamental limits on radiative (2) second harmonic generation.
Optics Express, 31(26):44212–44223, 2023.
Strekha et al. [2024a]
Benjamin Strekha, Pengning Chao, Rodrick Kuate Defo, Sean Molesky, and Alejandro W. Rodriguez.
Suppressing electromagnetic local density of states via slow light in lossy quasi-one-dimensional gratings.
Physical Review A, 109(4):L041501, 2024a.
Amaolo et al. [2024b]
Alessio Amaolo, Pengning Chao, Thomas J. Maldonado, Sean Molesky, and Alejandro W. Rodriguez.
Physical limits on Raman scattering: The critical role of pump and signal co-design.
arXiv preprint arXiv:2403.03332, 2024b.
Venkataram et al. [2020b]
Prashanth S. Venkataram, Sean Molesky, Pengning Chao, and Alejandro W. Rodriguez.
Fundamental limits to attractive and repulsive Casimir-Polder forces.
Physical Review A, 101(5):052115, 2020b.
Strekha et al. [2022]
Benjamin Strekha, Sean Molesky, Pengning Chao, Matthias Krüger, and Alejandro W. Rodriguez.
Trace expressions and associated limits for nonequilibrium Casimir torque.
Physical Review A, 106(4):042222, 2022.
Strekha et al. [2024b]
Benjamin Strekha, Matthias Krüger, and Alejandro W. Rodriguez.
Trace expressions and associated limits for equilibrium Casimir torque.
Physical Review A, 109(1):012813, 2024b.
Fleury et al. [2015]
Romain Fleury, Francesco Monticone, and Andrea Alù.
Invisibility and cloaking: Origins, present, and future perspectives.
Physical Review Applied, 4(3):037001, 2015.
doi: 10.1103/PhysRevApplied.4.037001.
Monticone and Alù [2013]
Francesco Monticone and Andrea Alù.
Do cloaked objects really scatter less?
Physical Review X, 3(4):041005, 2013.
Hashemi et al. [2011]
Hila Hashemi, A. Oskooi, J. D. Joannopoulos, and Steven G. Johnson.
General scaling limitations of ground-plane and isolated-object cloaks.
Physical Review A, 84(2):023815, August 2011.
ISSN 1050-2947, 1094-1622.
doi: 10.1103/PhysRevA.84.023815.
Molesky et al. [2021]
Sean Molesky, Pengning Chao, Jewel Mohajan, Wesley Reinhart, Heng Chi, and Alejandro W. Rodriguez.
-Operator Limits on Optical Communication: Metaoptics, Computation, and Input-Output Transformations.
Physical Review Research, 2021.
Shim et al. [2024]
Hyungki Shim, Zeyu Kuang, Zin Lin, and Owen D. Miller.
Fundamental limits to multi-functional and tunable nanophotonic response.
Nanophotonics, 13(12):2107–2116, 2024.
Johnson et al. [2019]
S. G. Johnson et al.
The NLopt Nonlinear Optimization Package (Version 2.6.2).
August 2019.
Christiansen and Sigmund [2021]
Rasmus E. Christiansen and Ole Sigmund.
Inverse design in photonics by topology optimization: Tutorial.
JOSA B, 38(2):496–509, 2021.
Svanberg [2002]
Krister Svanberg.
A class of globally convergent optimization methods based on conservative convex separable approximations.
SIAM Journal on Optimization, 12(2):555, 2002.