## Abstract

Numerical optimization of photonic devices is often limited by a large design space the finite-differences gradient method requires as many electric field computations as there are design parameters. Adjoint-based optimization can deliver the same gradients with only two electric field computations. Here, we derive the relevant adjoint formalism and illustrate its application for a waveguide slab, and for the design of optical sub-wavelength gratings.

© 2014 Optical Society of America

## 1. Introduction

When designing optical structures via numerical simulations, we are often lead to compute the gradient of some target function with respect to our design variables: we may, for example, want to assess the sensitivity of our structure to manufacturing imprecisions, or optimize the structure’s performance with respect to a given metric. Computing the gradient of a function usually requires two optical simulations per design parameter. This is fine for simple systems, but intractable for many optical systems of interest (e.g. diffractive elements or photonics circuits) which can contain millions of independent parameters.

Fortunately, if the target function is specified in advance, it is possible to exploit the mathematical structure of a large class of optical simulations to compute the multivariate gradient using only two simulations, no matter how many parameters are present. Next to the primal solve, one has to solve the adjoint problem once. The adjoint solution acts as a kind of influence function regarding the respective target function and avoids to solve the primal problem for each parameter variation. This is the principle idea behind adjoint-based optimization and sensitivity analysis. This approach is well-known to several fields of application as control theory, mesh adaptation, error estimation or propagation [14] and has been used very successfully for many years in performing aerodynamic shape optimizations [5,6]. Later, it was used by electrical engineers in the microwave regime [710]. In the optical arena, the approach has been used by only a few groups [1115].

The potential of this approach to optimize complex optical structures with many design parameters can hardly be overstated. For example, optimizing optical lens gratings [16] typically requires supercomputers or very large clusters. Countless runs with slightly different parameter settings are performed to find an appropriate set of design parameters. Being able to reduce the required number of these computations significantly relaxes the computational requirements. Thus, instead of probing a huge number of parameter settings, we can probe only a few and then perform efficient local optimization in the vicinity of the best settings. Adjoint-based optimization is an elegant method to compute the gradients that are needed for local optimization. This paper shows how to reduce the computational complexity of finding the gradient tremendously: from approximately as many field computations as there are design parameters (which is often in the hundreds or thousands) to two, independently of the number of design parameters. This makes it possible to optimize rather complex photonic systems with respect to an arbitrary number of parameters on a single desktop or laptop computer.

It is worth noticing that there are several different approaches to optimizing photonic structures and devices [17]. One particularly elegant approach is inverse design [18], where the structure is optimized to fit a particular field. In any case, when an optimization approach uses local gradients, adjoints may be used to efficiently compute them.

This paper derives the relevant equations for adjoint-based optimization and sensitivity analysis for optical systems. In particular, we present a formalism that is well suited for sensitivity analysis and optimization using a Maxwell solver. Therefore, this paper provides a detailed tutorial on the adjoint method used to solve electromagnetic problems. We provide a simple example of a slab waveguide where the adjoint method can be compared to the exact analytic model. We also show that the adjoint method can be effectively applied to the optimization of a grating coupler – a task where the adjoint method has not been applied before.

## 2. Problem description

In all generality, adjoint-based optimization is used to locally optimize a target function T (x, p) which depends on the physical state of the system x, as well as a number of design parameters p. For optical gratings, for example, the physical state of the system can be described conveniently by the electric field (or sometimes just a single component thereof) at each point of the computational grid. The vector of design parameters p contains, for example, geometric parameters describing the layout of our grating. In order to determine the state of our system x(p), we solve a given vector equation R(x, p) = 0, such as Maxwell’s equations.

For our optimization, we need to know how each of the design parameters affects the value of our target function: we need to compute $dTdp$. For each design parameter pk, this derivative is

$dTdpk=∑i∂T∂xidxidpk+∂T∂pk=∂T∂xdxdpk+∂T∂pk,$
provided that the design parameters are mutually independent. What makes the evaluation of the above equation computationally costly is the the term $dxdpk$. In fact, it is important to understand that the dependency of our solution x on the design parameters p is implicit via the condition that R(x, p) = 0. Otherwise, we would be able to analytically compute the relevant derivatives with respect to our design parameters and our optimization problem would be trivial.

Computing $dxdp$ via the finite-differences approximation is numerically very inefficient:

$dxdpk≈1Δpk(x(pΔpk−x(p)),$
where Δpk is a small variation on the k-th design variable, and pΔpk denotes the vector of design parameters where Δpk was added to the k-th component only. As we can see from the above finite-difference equation, this approach demands for one reference computation of R(x, p) = 0 to compute x(p), and then one additional R(x, pΔpk) = 0 to find the necessary x(pΔpk) for each component of the gradient. To compute the full gradient with respect to n design variables, we would therefore have to solve the vector equation R = 0 a total of n + 1 times. For costly R and large numbers of design parameters, this often exceeds the available resources of computational power and time.

As we show below, it is often possible to exploit the mathematical structure of the local optimization and sensitivity analysis problem: the adjoint method allows us to compute the relevant gradient $dTdp$ performing only two non-trivial simulations. No matter how many design variables we are dealing with.

## 3. Derivation of adjoint approach

Using the adjoint method, we are able to determine $dTdp$ using only one solution of R(x, p) = 0, and one solution of the adjoint equation. As long as the latter is is of similar numerical complexity as the direct equation R, the adjoint-based optimization is likely to reduce the computational cost of $dTdp$, which we are interested in.

Our starting point is the fact that, for any set of design parameters p, and any corresponding state x, we know that the governing vector equation R vanishes,

$R(x,p)=0.$
Since this condition has to hold for any set of parameters p with corresponding x, its Taylor expansion also has to vanish. Intuitively speaking, in the solution space, where we only consider design parameters and their corresponding states, the derivative of the state function R with respect to a design parameter pk also vanishes,
$dRdpk=∂R∂xdxdpk+∂R∂pk=0.$
We can now multiply $dRdpk$ on the left by any complex vector v and subtract the product from $dTdpk$ without changing the latter gradient’s value
$dTdpk=∂T∂xdxdpk+∂T∂pk−v†(∂R∂xdxdpk+∂R∂pk).$
By regrouping the terms, we find
$dTdpk=∂T∂pk−v†∂R∂pk+(∂T∂x−v†∂R∂x) dxdpk.$
This expression is the key to the adjoint method: if we are able to find a v such that the prefactor of $dxdpk$ vanishes, we do not have to perform the expensive n + 1 solutions of our state function to find the gradient $dTdp$. In fact, we will simply have to solve R = 0 once to find our x(p). Then, we will solve our adjoint equation
$∂T∂x=v†∂R∂x,$
to find a convenient v. And these two computations then straightforwardly allow us to determine the derivative of our target function with respect to our design parameters,
$dTdpk=∂T∂pk−v†∂R∂pk.$

To illustrate the use of this powerful approach, we will now show how to apply the adjoint method in two concrete situations and derive the relevant equations for optical systems.

## 4. Optical systems

In this section, we derive the relevant equations for adjoint-based optimization and sensitivity analysis in optical systems, starting the same way as Nikolova et al. [9] did for the microwave regime. For mathematical simplicity, we consider the dispersion-free, and linear case. Nonetheless, adjoint-based optimization can also be used for non-uniform, anisotropic media. Maxwell’s equations are

$∇×H=∂D∂t+Jσ+J−∇×E=∂B∂t,$
where J denotes the density of current sources. The constitutive relations are
$B=μHD=εEJσ=σE,$
where μ denotes the permeability tensor, ε the permittivity tensor, and σ the specific conductivity tensor. These equations can be combined into
$∇×μ−1∇×E+ε∂2E∂t2+σ∂E∂t=−∂J∂t.$

In our optical gratings, the specific conductivity tensor σ vanishes. Furthermore, we write our electric field (and source) as E(t) = Re(Ēe−iωt), where Ē ∈ ℂ. Note that the permittivity ε describes the layout of our optical grating and therefore directly depends on our design parameters. We rewrite Eq (11) as

$R(E¯,p)=ℳE¯−iωJ¯=0,$
where we introduced
$ℳE¯=∇×μ−1∇×E¯−ω2εE¯.$
To compute an electric field Ē that satisfies Eq. (11), we can use any convenient Maxwell solver, such as any FDTD [19, 20] solver. The simulation returns the complex electric field amplitude at each point of the computational grid.

Our target function T is the electric field strength in a specific region Ω

$T(E¯,p)=1N∫ΩdΩ|E¯|2=E¯0†E¯,$
where Ē0 = Ē within the region Ω and zero outside, and N is some (irrelevant) normalization constant introduced for mathematical rigor, only. With this notation, we also introduce the additional hypothesis that Ē varies predominantly in norm (rather than in phase), as we vary p. This will allow us to approximate the derivative of our target function with respect to the electric field as $E¯0†$.

At this point, we are able to derive the adjoint equation for this problem. First, note that for isotropic dielectric media where ε = ε and μ = μ we have = . Second, we use Ē computed as the numerical solution of Eq. (12) in Eq. (7) and find $E¯0†=v¯†ℳ$. This equation will allow us to compute the adjoint field . Third, we rewrite this equation in a way that shows its correspondence to the direct problem Eq. (12) and conclude that

$ℳv¯−iω(1iωE¯0)=0.$

Both the direct problem Eq. (12) and the adjoint problem Eq. (15) can be simulated using the same Maxwell solver. In fact, we see that the adjoint field can be interpreted as an electric field that was created by a particular set of sources. The adjoint sources term, $1iωE¯0$, is proportional to the field of the direct solution, wherever we are optimizing the field, and zero everywhere else.

Computing the relevant gradient Eq. (8) is computationally trivial once we have the solutions of the direct problem, Ē, and of the adjoint problem, . Our target function (14) does not depend on the design parameters p. Hence, $∂T∂pk=0$. Also, the only element of the governing vector equation (12) that depends on the design parameters is the spatially varying dielectric constant ε. Therefore, $∂R∂pk=∂ℳ∂pkE¯=∂ε∂pkE¯$ and the gradient Eq. (8) is simply

$dTdpk=−v†∂ε∂pkE¯.$
Note that this expression includes simple vector-matrix multiplication and finding the derivative of our spatially varying dielectric constant ε with respect to each of the design parameters pk. However, we certainly know the current dielectric layout of the device we are trying to optimize. Hence, we can compute $∂ε∂pk$ either analytically, or by using computationally cheap numerical methods. This shows that we can ultimately compute the derivative of our target function with respect to our design parameters, $dTdpk$, with only two non-trivial computations: to find the direct solution Ē, and the adjoint solution .

## 5. Example 1: waveguide slab

Consider the text-book example of a waveguide slab of thickness h, described by the cover material index nc, the guiding film index nf, and the substrate index ns. We would like to know what modes are supported by this waveguide slab for a given wavelength λ. The vacuum wave vector is $k0=2πλ$. The characteristic equation for transverse electric (TE) modes is

$tan(hκf)=γc+γsκf (1−γcγsκf2),$
where $κf=k02nf2−β2$ is the transverse wave vector, $γc=β2−k02nc2$ and $γs=β2−k02ns2$ are the attenuation coefficients of cover and substrate, respectively, which all depend on the longitudinal wave vector β. Therefore, solving Eq. (17) means finding one or all β for which this equation holds.

If we were to ask how a specific value of β changes as we modify, for example, the height of the slab, we would usually have to solve Eq. (17) for different values of h. Worse, even, if we were to modify the indices nc, nf, and ns, too, we would have to solve that same equation over and over again to derive how β depends on them. Nowadays computers can solve this transcendental equation really easily, but for the sake of explaining the adjoint method, pretend that this is an arduous task.

To use the adjoint method, we rewrite Eq. (17) as

$R=tan(hκf)−γc+γsκf (1−γcγsκf2),$
knowing that the value of R(β, p) depends on the longitudinal wave vector, as well as a number of design parameters p = (h, nc, nf, ns)T, and vanishes for all solutions β (p). For the sake of simplicity, we shall define our target function T to be simply the longitudinal wave vector,
$T=β,$

We know from the above Eq. (6) that we can avoid these repetitive calculations if we find a v such that

$∂T∂β=1=v∂R∂β.$
In this case,
$dTdpk=dβdpk=−v∂R∂pk.$

The derivative of our characteristic equation with respect to the longitudinal wave vector is

$∂R∂β=∂∂β(tan(hκf))−∂∂β(γc+γsκf (1−γcγsκf2)),$
where
$∂∂β(tan(hκf))=−hβκfcos2(hκf),$
and, introducing $q=κf(1−γcγsκf2)$,
$∂∂β(γc+γsq)=βγc+βγsq−γc+γsq2∂q∂β,$
where
$∂q∂β=−βκf[1+γsγc+γcγs+γcγsκf2].$
Collecting all these terms, we find the solution to our adjoint Eq. (20),
$1v=−β[hκfcos2(hκf)+(γc+γs)κfκf2−γcγs(1γsγc+1+γsγc+γcγs+γcγsκf2κf2−γcγs)],$
and we know that by defining
$v=1∂R∂β$
we can compute the relevant gradient as
$dβdpk=−v∂R∂pk.$
For example, if we are interested in computing how β changes as we vary the height of the guiding film, pk = h, and we therefore need to compute
$dRdh=κfcos2(hκf)$
and finally, following Eq. (28),
$dβdh=−vκfcos2(hκf).$
As we can see in Fig. 1, the adjoint method leads to the correct gradient. However, is far more efficient for computing multi-variable optimizations than finite-differences gradients.

Fig. 1 The derivative of the longitudinal wave vector β with respect to the thickness of a waveguide slab can faithfully be reproduced using adjoints: the comparison between the direct calculation (blue line) and adjoint calculation (red dots) of $ddhβ$ shows that the two methods deliver exactly the same result. However, the adjoint method requires no further expensive solutions of the eigenvalue system in order to determine the derivatives of β with respect to other design variables. Thus, the adjoint method delivers the relevant gradients of a system depending on multiple design parameters faster than the finite-differences approach. The inset shows a schematic of the waveguide slab structure comprised of the cover material nc, the guiding film nf, and the substrate ns.

## 6. Example 2: grating coupler

Silicon grating couplers are used in photonic circuits to couple light between an on-chip waveguide and an off-chip optical fiber [21]. The exact groove pattern (width, position, depth of each groove) can be optimized to provide the best possible coupling efficiency [22]. For the sake of illustrating the adjoint method, we will restrict the design space to two variables, the uniform pitch p and duty cycle d of the grating grooves. Our simulations are assuming a waveguide thickness of 0.22μm, a groove-depth of 0.08μm, a wavelength of 1.55μm.

The vector of our design parameters is

$p=(pd)$
Figure 2 shows the effect of these design parameters on the design of our grating coupler.

Fig. 2 Our SiO2 grating coupler redirects free-space plane-waves (vertical arrows) to a guided mode inside the waveguide (horizontal arrows). We show the device at (a) a pitch of 0.58 and a duty cycle of 0.2, and (b) at a pitch of 0.68 and a duty cycle of 0.8. Both the pitch and the duty cycle are optimized to maximize the coupling between the incident light and the light inside the waveguide.

To illustrate the effect of adjoint-based optimization, we computed the full topology of our optimization problem for this example (colored background in Fig. 3). We see the areas of best coupling in dark red, and those of worst coupling in dark blue. Of course, computing the full topology is computationally impossible for most real-world applications. Here, it serves intuition and helps to verify that our algorithm works, indeed, as steepest ascent. Furthermore, it helped us select rather poor starting points for our optimization.

Fig. 3 Adjoint-based optimization convergence: from three initial parameters (stars) in parameter regions of poor optical coupling (blue, green, yellow), adjoint-based optimization evolves the design parameters to the parameter regions of optimal coupling (dark red). The computation of the relevant gradients requires only two Maxwell solutions, independently of the number of design parameters. The finite differences approach involves roughly one Maxwell solution per design variable. Thus, this figure shows that adjoints help optimize the system faster than the finite-differences approach.

To perform the actual adjoint-based optimization, we select a number of starting points that serve as initial parameters. For each point, we

1. compute the direct electromagnetic field for specific design parameters (12),
2. compute the adjoint field for the same design parameters (15),
3. combine the direct and adjoint solutions to (trivially) compute the gradient (16),
4. update the design parameters for steepest ascent,
5. repeat steps 1–4 until convergence.
Note that sensitivity analysis can be performed by executing steps 1–3 for any desired set of design parameters.

Figure 3 shows the result of the full optimization (steps 1–5), and the (usually unknown) topology of our problem. In order to compute these optimization runs, finite-differences gradient optimization would have required three non-trivial field computations per gradient. Adjoint-based optimization, in contrast, only required two non-trivial field computations to find the same gradient. Thus, adjoint based optimization is already slightly faster than finite-differences based optimization for two design parameters.

Consider a grating optimization starting at a pitch of 0.62 and a duty cycle of 0.73 (black star). Adjoint-based optimization evolves these initial parameters to a pitch of 0.60 and a duty cycle of 0.275. Figure 4 shows the result of Meep [23] simulations of real(Ē) before and after optimization. The coupling of the incident light from the top to inside the waveguide (inside the black square) has improved from 0.39 to 0.88.

Fig. 4 Electric field in and around the waveguide (outlined in light blue) before and after adjoint-based optimization: (a) the grating with the initial design parameters couples the incident light (top) poorly to the waveguide (inside black box), with a complex amplitude of 0.39 (b) After adjoint-based optimization, the grating couples the light very well to the waveguide (inside black box), with a complex amplitude of 0.88.

For n design parameters, the advantage of adjoint-based optimization would be even more striking: direct computation of the finite-differences gradient would require at least n + 1 field computations per gradient, whereas the adjoint method still only requires two.

## 7. Conclusion

We have shown how to use the adjoint method to efficiently design of optical structures. The key of this method is to compute the gradient of the target function with respect to the design parameters using only two non-trivial computations, independently of the number of design parameters. We illustrated this method in two cases: the optical waveguide slab, and a grating coupler simulated using a standard Maxwell solver. For the latter, the converged electric field Ē and the corresponding adjoint field are the only non-trivial computations required to calculate the gradient, needed to perform gradient ascent (or gradient descent) on the design parameters. We believe this method will prove itself immensely valuable in the design of various grating structures, including a host of non-periodic diffractive structures.

## Acknowledgments

A.N. thanks Hideo Mabuchi, Erling Riis, and Thomas Baer for their mentorship during his time at Stanford University. Furthermore, he thanks Charles Santori, Marco Fiorentino, Sonny Vo, and Wayne Sorin for interesting discussions on how to design optical devices. This work was supported by the SU2P Program of Stanford University and Research Councils UK.

1. J. L. Lions, Optimal Control of Systems Governed by Partial Differential Equations (Springer, 1971). [CrossRef]

2. R. Becker and R. Rannacher, “An optimal control approach to error control and mesh adaption in finite element methods,” Acta Numerica 10, 1–102 (2001). [CrossRef]

3. K. Eriksson, D. Estep, P. Hansbo, and C. Johnson, “Introduction to adaptive methods for differential equations,” Acta Numerica 4, 105–158 (1995). [CrossRef]

4. D. E. Rumelhart, G. E. Hinton, and R. J. Williams, “Learning internal representations by error propagation,” Parallel Data Processing 1, 318–362 (1986).

5. J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 1,” Journal of Aircraft 36(1), 51–60 (1999). [CrossRef]

6. J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 2,” Journal of Aircraft 36(1), 61–74 (1999). [CrossRef]

7. Y. seek Chung, Changyul-Cheon, I.-H. Park, and S.-Y. Hahn, “Optimal shape design of microwave device using fdtd and design sensitivity analysis,” Microwave Theory and Techniques, IEEE Transactions on 48, 2289–2296 (2000). [CrossRef]

8. N. Georgieva, S. Glavic, M. Bakr, and J. Bandler, “Feasible adjoint sensitivity technique for em design optimization,” Microwave Theory and Techniques, IEEE Transactions on 50, 2751–2758 (2002). [CrossRef]

9. N. K. Nikolova, H. W. Tam, and M. H. Bakr, “Sensitivity analysis with the fdtd method on structured grids,” Microwave Theory and Techniques, IEEE Transactions on 52, 1207–1216 (2004). [CrossRef]

10. N. K. Nikolova, Y. Li, Y. Li, and M. H. Bakr, “Sensitivity analysis of scattering parameters with electromagnetic time-domain simulators,” Microwave Theory and Techniques, IEEE Transactions on 54, 1598–1610 (2006). [CrossRef]

11. G. Veronis, R. W. Dutton, and S. Fan, “Method for sensitivity analysis of photonic crystal devices,” Optics Letters 29, 2288–2290 (2004). [CrossRef]   [PubMed]

12. Y. Jiao, S. Fan, and D. A. B. Miller, “Photonic crystal device sensitivity analysis with wannierbasis gradients,” Optics Letters 30, 302–304 (2005). [CrossRef]

13. P. Seliger, M. Mahvash, C. Wang, and A. F. J. Levi, “Optimization of aperiodic dielectric structures,” Journal of Applied Physics 100, 034310 (2006). [CrossRef]

14. C. M. Lalau-Keraly, S. Bhargava, O. D. Miller, and E. Yablonovitch, “Adjoint shape optimization applied to electromagnetic design,” Optics Express 21, 21693–21701 (2013). [CrossRef]   [PubMed]

15. O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014). [CrossRef]   [PubMed]

16. D. Fattal, J. Li, Z. Peng, M. Fiorentino, and R. G. Beausoleil, “Flat dielectric grating reflectors with focusing abilities,” Nature Photonics 4, 466–470 (2010). [CrossRef]

17. V. Liu, D. Miller, and S. Fan, “Highly tailored computational electromagnetics methods for nanophotonic design and discovery,” Proceedings of the IEEE 101, 484–493 (2013). [CrossRef]

18. J. Lu and Vučković, “Inverse design of nanophotonic structures using complementary convex optimization,” Optics Express 18, 3793–3804 (2010). [CrossRef]   [PubMed]

19. R. Courant, K. Friedrichs, and H. Lewy, “Über die partiellen Differenzengleichungen der mathematischen Physik,” Mathematische Annalen 100, 32–74 (1928). [CrossRef]

20. K. Yee, “Numerical solution of initial boundary value problems involving maxwell’s equations in isotropic media,” Antennas and Propagation, IEEE Transactions on 14, 302–307 (1966). [CrossRef]

21. D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002). [CrossRef]

22. G. Roelkens, D. V. Thourhout, and R. Baets, “High efficiency silicon-on-insulator grating coupler based on a poly-silicon overlay,” Optics Express 14, 11622–11630 (2006). [CrossRef]   [PubMed]

23. A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Computer Physics Communications 181, 687–702 (2010). [CrossRef]

### References

• View by:
• |
• |
• |

1. J. L. Lions, Optimal Control of Systems Governed by Partial Differential Equations (Springer, 1971).
[Crossref]
2. R. Becker and R. Rannacher, “An optimal control approach to error control and mesh adaption in finite element methods,” Acta Numerica 10, 1–102 (2001).
[Crossref]
3. K. Eriksson, D. Estep, P. Hansbo, and C. Johnson, “Introduction to adaptive methods for differential equations,” Acta Numerica 4, 105–158 (1995).
[Crossref]
4. D. E. Rumelhart, G. E. Hinton, and R. J. Williams, “Learning internal representations by error propagation,” Parallel Data Processing 1, 318–362 (1986).
5. J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 1,” Journal of Aircraft 36(1), 51–60 (1999).
[Crossref]
6. J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 2,” Journal of Aircraft 36(1), 61–74 (1999).
[Crossref]
7. Y. seek Chung, Changyul-Cheon, I.-H. Park, and S.-Y. Hahn, “Optimal shape design of microwave device using fdtd and design sensitivity analysis,” Microwave Theory and Techniques, IEEE Transactions on 48, 2289–2296 (2000).
[Crossref]
8. N. Georgieva, S. Glavic, M. Bakr, and J. Bandler, “Feasible adjoint sensitivity technique for em design optimization,” Microwave Theory and Techniques, IEEE Transactions on 50, 2751–2758 (2002).
[Crossref]
9. N. K. Nikolova, H. W. Tam, and M. H. Bakr, “Sensitivity analysis with the fdtd method on structured grids,” Microwave Theory and Techniques, IEEE Transactions on 52, 1207–1216 (2004).
[Crossref]
10. N. K. Nikolova, Y. Li, Y. Li, and M. H. Bakr, “Sensitivity analysis of scattering parameters with electromagnetic time-domain simulators,” Microwave Theory and Techniques, IEEE Transactions on 54, 1598–1610 (2006).
[Crossref]
11. G. Veronis, R. W. Dutton, and S. Fan, “Method for sensitivity analysis of photonic crystal devices,” Optics Letters 29, 2288–2290 (2004).
[Crossref] [PubMed]
12. Y. Jiao, S. Fan, and D. A. B. Miller, “Photonic crystal device sensitivity analysis with wannierbasis gradients,” Optics Letters 30, 302–304 (2005).
[Crossref]
13. P. Seliger, M. Mahvash, C. Wang, and A. F. J. Levi, “Optimization of aperiodic dielectric structures,” Journal of Applied Physics 100, 034310 (2006).
[Crossref]
14. C. M. Lalau-Keraly, S. Bhargava, O. D. Miller, and E. Yablonovitch, “Adjoint shape optimization applied to electromagnetic design,” Optics Express 21, 21693–21701 (2013).
[Crossref] [PubMed]
15. O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]
16. D. Fattal, J. Li, Z. Peng, M. Fiorentino, and R. G. Beausoleil, “Flat dielectric grating reflectors with focusing abilities,” Nature Photonics 4, 466–470 (2010).
[Crossref]
17. V. Liu, D. Miller, and S. Fan, “Highly tailored computational electromagnetics methods for nanophotonic design and discovery,” Proceedings of the IEEE 101, 484–493 (2013).
[Crossref]
18. J. Lu and Vučković, “Inverse design of nanophotonic structures using complementary convex optimization,” Optics Express 18, 3793–3804 (2010).
[Crossref] [PubMed]
19. R. Courant, K. Friedrichs, and H. Lewy, “Über die partiellen Differenzengleichungen der mathematischen Physik,” Mathematische Annalen 100, 32–74 (1928).
[Crossref]
20. K. Yee, “Numerical solution of initial boundary value problems involving maxwell’s equations in isotropic media,” Antennas and Propagation, IEEE Transactions on 14, 302–307 (1966).
[Crossref]
21. D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]
22. G. Roelkens, D. V. Thourhout, and R. Baets, “High efficiency silicon-on-insulator grating coupler based on a poly-silicon overlay,” Optics Express 14, 11622–11630 (2006).
[Crossref] [PubMed]
23. A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Computer Physics Communications 181, 687–702 (2010).
[Crossref]

#### 2014 (1)

O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]

#### 2013 (2)

V. Liu, D. Miller, and S. Fan, “Highly tailored computational electromagnetics methods for nanophotonic design and discovery,” Proceedings of the IEEE 101, 484–493 (2013).
[Crossref]

C. M. Lalau-Keraly, S. Bhargava, O. D. Miller, and E. Yablonovitch, “Adjoint shape optimization applied to electromagnetic design,” Optics Express 21, 21693–21701 (2013).
[Crossref] [PubMed]

#### 2010 (3)

J. Lu and Vučković, “Inverse design of nanophotonic structures using complementary convex optimization,” Optics Express 18, 3793–3804 (2010).
[Crossref] [PubMed]

A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Computer Physics Communications 181, 687–702 (2010).
[Crossref]

D. Fattal, J. Li, Z. Peng, M. Fiorentino, and R. G. Beausoleil, “Flat dielectric grating reflectors with focusing abilities,” Nature Photonics 4, 466–470 (2010).
[Crossref]

#### 2006 (3)

P. Seliger, M. Mahvash, C. Wang, and A. F. J. Levi, “Optimization of aperiodic dielectric structures,” Journal of Applied Physics 100, 034310 (2006).
[Crossref]

N. K. Nikolova, Y. Li, Y. Li, and M. H. Bakr, “Sensitivity analysis of scattering parameters with electromagnetic time-domain simulators,” Microwave Theory and Techniques, IEEE Transactions on 54, 1598–1610 (2006).
[Crossref]

G. Roelkens, D. V. Thourhout, and R. Baets, “High efficiency silicon-on-insulator grating coupler based on a poly-silicon overlay,” Optics Express 14, 11622–11630 (2006).
[Crossref] [PubMed]

#### 2005 (1)

Y. Jiao, S. Fan, and D. A. B. Miller, “Photonic crystal device sensitivity analysis with wannierbasis gradients,” Optics Letters 30, 302–304 (2005).
[Crossref]

#### 2004 (2)

N. K. Nikolova, H. W. Tam, and M. H. Bakr, “Sensitivity analysis with the fdtd method on structured grids,” Microwave Theory and Techniques, IEEE Transactions on 52, 1207–1216 (2004).
[Crossref]

G. Veronis, R. W. Dutton, and S. Fan, “Method for sensitivity analysis of photonic crystal devices,” Optics Letters 29, 2288–2290 (2004).
[Crossref] [PubMed]

#### 2002 (2)

N. Georgieva, S. Glavic, M. Bakr, and J. Bandler, “Feasible adjoint sensitivity technique for em design optimization,” Microwave Theory and Techniques, IEEE Transactions on 50, 2751–2758 (2002).
[Crossref]

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### 2001 (1)

R. Becker and R. Rannacher, “An optimal control approach to error control and mesh adaption in finite element methods,” Acta Numerica 10, 1–102 (2001).
[Crossref]

#### 2000 (1)

Y. seek Chung, Changyul-Cheon, I.-H. Park, and S.-Y. Hahn, “Optimal shape design of microwave device using fdtd and design sensitivity analysis,” Microwave Theory and Techniques, IEEE Transactions on 48, 2289–2296 (2000).
[Crossref]

#### 1999 (2)

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 1,” Journal of Aircraft 36(1), 51–60 (1999).
[Crossref]

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 2,” Journal of Aircraft 36(1), 61–74 (1999).
[Crossref]

#### 1995 (1)

K. Eriksson, D. Estep, P. Hansbo, and C. Johnson, “Introduction to adaptive methods for differential equations,” Acta Numerica 4, 105–158 (1995).
[Crossref]

#### 1986 (1)

D. E. Rumelhart, G. E. Hinton, and R. J. Williams, “Learning internal representations by error propagation,” Parallel Data Processing 1, 318–362 (1986).

#### 1966 (1)

K. Yee, “Numerical solution of initial boundary value problems involving maxwell’s equations in isotropic media,” Antennas and Propagation, IEEE Transactions on 14, 302–307 (1966).
[Crossref]

#### 1928 (1)

R. Courant, K. Friedrichs, and H. Lewy, “Über die partiellen Differenzengleichungen der mathematischen Physik,” Mathematische Annalen 100, 32–74 (1928).
[Crossref]

#### Alonso, J. J.

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 1,” Journal of Aircraft 36(1), 51–60 (1999).
[Crossref]

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 2,” Journal of Aircraft 36(1), 61–74 (1999).
[Crossref]

#### Baets, R.

G. Roelkens, D. V. Thourhout, and R. Baets, “High efficiency silicon-on-insulator grating coupler based on a poly-silicon overlay,” Optics Express 14, 11622–11630 (2006).
[Crossref] [PubMed]

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### Bakr, M.

N. Georgieva, S. Glavic, M. Bakr, and J. Bandler, “Feasible adjoint sensitivity technique for em design optimization,” Microwave Theory and Techniques, IEEE Transactions on 50, 2751–2758 (2002).
[Crossref]

#### Bakr, M. H.

N. K. Nikolova, Y. Li, Y. Li, and M. H. Bakr, “Sensitivity analysis of scattering parameters with electromagnetic time-domain simulators,” Microwave Theory and Techniques, IEEE Transactions on 54, 1598–1610 (2006).
[Crossref]

N. K. Nikolova, H. W. Tam, and M. H. Bakr, “Sensitivity analysis with the fdtd method on structured grids,” Microwave Theory and Techniques, IEEE Transactions on 52, 1207–1216 (2004).
[Crossref]

#### Bandler, J.

N. Georgieva, S. Glavic, M. Bakr, and J. Bandler, “Feasible adjoint sensitivity technique for em design optimization,” Microwave Theory and Techniques, IEEE Transactions on 50, 2751–2758 (2002).
[Crossref]

#### Beausoleil, R. G.

D. Fattal, J. Li, Z. Peng, M. Fiorentino, and R. G. Beausoleil, “Flat dielectric grating reflectors with focusing abilities,” Nature Photonics 4, 466–470 (2010).
[Crossref]

#### Becker, R.

R. Becker and R. Rannacher, “An optimal control approach to error control and mesh adaption in finite element methods,” Acta Numerica 10, 1–102 (2001).
[Crossref]

#### Bermel, P.

A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Computer Physics Communications 181, 687–702 (2010).
[Crossref]

#### Bhargava, S.

C. M. Lalau-Keraly, S. Bhargava, O. D. Miller, and E. Yablonovitch, “Adjoint shape optimization applied to electromagnetic design,” Optics Express 21, 21693–21701 (2013).
[Crossref] [PubMed]

#### Bienstman, P.

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### Bogaerts, W.

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### Changyul-Cheon,

Y. seek Chung, Changyul-Cheon, I.-H. Park, and S.-Y. Hahn, “Optimal shape design of microwave device using fdtd and design sensitivity analysis,” Microwave Theory and Techniques, IEEE Transactions on 48, 2289–2296 (2000).
[Crossref]

#### Courant, R.

R. Courant, K. Friedrichs, and H. Lewy, “Über die partiellen Differenzengleichungen der mathematischen Physik,” Mathematische Annalen 100, 32–74 (1928).
[Crossref]

#### De Mesel, K.

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### DeLacy, B. G.

O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]

#### Dutton, R. W.

G. Veronis, R. W. Dutton, and S. Fan, “Method for sensitivity analysis of photonic crystal devices,” Optics Letters 29, 2288–2290 (2004).
[Crossref] [PubMed]

#### Eriksson, K.

K. Eriksson, D. Estep, P. Hansbo, and C. Johnson, “Introduction to adaptive methods for differential equations,” Acta Numerica 4, 105–158 (1995).
[Crossref]

#### Estep, D.

K. Eriksson, D. Estep, P. Hansbo, and C. Johnson, “Introduction to adaptive methods for differential equations,” Acta Numerica 4, 105–158 (1995).
[Crossref]

#### Fan, S.

V. Liu, D. Miller, and S. Fan, “Highly tailored computational electromagnetics methods for nanophotonic design and discovery,” Proceedings of the IEEE 101, 484–493 (2013).
[Crossref]

Y. Jiao, S. Fan, and D. A. B. Miller, “Photonic crystal device sensitivity analysis with wannierbasis gradients,” Optics Letters 30, 302–304 (2005).
[Crossref]

G. Veronis, R. W. Dutton, and S. Fan, “Method for sensitivity analysis of photonic crystal devices,” Optics Letters 29, 2288–2290 (2004).
[Crossref] [PubMed]

#### Fattal, D.

D. Fattal, J. Li, Z. Peng, M. Fiorentino, and R. G. Beausoleil, “Flat dielectric grating reflectors with focusing abilities,” Nature Photonics 4, 466–470 (2010).
[Crossref]

#### Fiorentino, M.

D. Fattal, J. Li, Z. Peng, M. Fiorentino, and R. G. Beausoleil, “Flat dielectric grating reflectors with focusing abilities,” Nature Photonics 4, 466–470 (2010).
[Crossref]

#### Friedrichs, K.

R. Courant, K. Friedrichs, and H. Lewy, “Über die partiellen Differenzengleichungen der mathematischen Physik,” Mathematische Annalen 100, 32–74 (1928).
[Crossref]

#### Georgieva, N.

N. Georgieva, S. Glavic, M. Bakr, and J. Bandler, “Feasible adjoint sensitivity technique for em design optimization,” Microwave Theory and Techniques, IEEE Transactions on 50, 2751–2758 (2002).
[Crossref]

#### Glavic, S.

N. Georgieva, S. Glavic, M. Bakr, and J. Bandler, “Feasible adjoint sensitivity technique for em design optimization,” Microwave Theory and Techniques, IEEE Transactions on 50, 2751–2758 (2002).
[Crossref]

#### Hahn, S.-Y.

Y. seek Chung, Changyul-Cheon, I.-H. Park, and S.-Y. Hahn, “Optimal shape design of microwave device using fdtd and design sensitivity analysis,” Microwave Theory and Techniques, IEEE Transactions on 48, 2289–2296 (2000).
[Crossref]

#### Hansbo, P.

K. Eriksson, D. Estep, P. Hansbo, and C. Johnson, “Introduction to adaptive methods for differential equations,” Acta Numerica 4, 105–158 (1995).
[Crossref]

#### Hinton, G. E.

D. E. Rumelhart, G. E. Hinton, and R. J. Williams, “Learning internal representations by error propagation,” Parallel Data Processing 1, 318–362 (1986).

#### Hsu, C. W.

O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]

#### Ibanescu, M.

A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Computer Physics Communications 181, 687–702 (2010).
[Crossref]

#### Jameson, A.

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 2,” Journal of Aircraft 36(1), 61–74 (1999).
[Crossref]

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 1,” Journal of Aircraft 36(1), 51–60 (1999).
[Crossref]

#### Jiao, Y.

Y. Jiao, S. Fan, and D. A. B. Miller, “Photonic crystal device sensitivity analysis with wannierbasis gradients,” Optics Letters 30, 302–304 (2005).
[Crossref]

#### Joannopoulos, J. D.

O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]

A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Computer Physics Communications 181, 687–702 (2010).
[Crossref]

#### Johnson, C.

K. Eriksson, D. Estep, P. Hansbo, and C. Johnson, “Introduction to adaptive methods for differential equations,” Acta Numerica 4, 105–158 (1995).
[Crossref]

#### Johnson, S. G.

O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]

A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Computer Physics Communications 181, 687–702 (2010).
[Crossref]

#### Krauss, T.

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### Lalau-Keraly, C. M.

C. M. Lalau-Keraly, S. Bhargava, O. D. Miller, and E. Yablonovitch, “Adjoint shape optimization applied to electromagnetic design,” Optics Express 21, 21693–21701 (2013).
[Crossref] [PubMed]

#### Levi, A. F. J.

P. Seliger, M. Mahvash, C. Wang, and A. F. J. Levi, “Optimization of aperiodic dielectric structures,” Journal of Applied Physics 100, 034310 (2006).
[Crossref]

#### Lewy, H.

R. Courant, K. Friedrichs, and H. Lewy, “Über die partiellen Differenzengleichungen der mathematischen Physik,” Mathematische Annalen 100, 32–74 (1928).
[Crossref]

#### Li, J.

D. Fattal, J. Li, Z. Peng, M. Fiorentino, and R. G. Beausoleil, “Flat dielectric grating reflectors with focusing abilities,” Nature Photonics 4, 466–470 (2010).
[Crossref]

#### Li, Y.

N. K. Nikolova, Y. Li, Y. Li, and M. H. Bakr, “Sensitivity analysis of scattering parameters with electromagnetic time-domain simulators,” Microwave Theory and Techniques, IEEE Transactions on 54, 1598–1610 (2006).
[Crossref]

N. K. Nikolova, Y. Li, Y. Li, and M. H. Bakr, “Sensitivity analysis of scattering parameters with electromagnetic time-domain simulators,” Microwave Theory and Techniques, IEEE Transactions on 54, 1598–1610 (2006).
[Crossref]

#### Lions, J. L.

J. L. Lions, Optimal Control of Systems Governed by Partial Differential Equations (Springer, 1971).
[Crossref]

#### Liu, V.

V. Liu, D. Miller, and S. Fan, “Highly tailored computational electromagnetics methods for nanophotonic design and discovery,” Proceedings of the IEEE 101, 484–493 (2013).
[Crossref]

#### Lu, J.

J. Lu and Vučković, “Inverse design of nanophotonic structures using complementary convex optimization,” Optics Express 18, 3793–3804 (2010).
[Crossref] [PubMed]

#### Mahvash, M.

P. Seliger, M. Mahvash, C. Wang, and A. F. J. Levi, “Optimization of aperiodic dielectric structures,” Journal of Applied Physics 100, 034310 (2006).
[Crossref]

#### Miller, D.

V. Liu, D. Miller, and S. Fan, “Highly tailored computational electromagnetics methods for nanophotonic design and discovery,” Proceedings of the IEEE 101, 484–493 (2013).
[Crossref]

#### Miller, D. A. B.

Y. Jiao, S. Fan, and D. A. B. Miller, “Photonic crystal device sensitivity analysis with wannierbasis gradients,” Optics Letters 30, 302–304 (2005).
[Crossref]

#### Miller, O. D.

O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]

C. M. Lalau-Keraly, S. Bhargava, O. D. Miller, and E. Yablonovitch, “Adjoint shape optimization applied to electromagnetic design,” Optics Express 21, 21693–21701 (2013).
[Crossref] [PubMed]

#### Moerman, I.

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### Nikolova, N. K.

N. K. Nikolova, Y. Li, Y. Li, and M. H. Bakr, “Sensitivity analysis of scattering parameters with electromagnetic time-domain simulators,” Microwave Theory and Techniques, IEEE Transactions on 54, 1598–1610 (2006).
[Crossref]

N. K. Nikolova, H. W. Tam, and M. H. Bakr, “Sensitivity analysis with the fdtd method on structured grids,” Microwave Theory and Techniques, IEEE Transactions on 52, 1207–1216 (2004).
[Crossref]

#### Oskooi, A. F.

A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Computer Physics Communications 181, 687–702 (2010).
[Crossref]

#### Park, I.-H.

Y. seek Chung, Changyul-Cheon, I.-H. Park, and S.-Y. Hahn, “Optimal shape design of microwave device using fdtd and design sensitivity analysis,” Microwave Theory and Techniques, IEEE Transactions on 48, 2289–2296 (2000).
[Crossref]

#### Peng, Z.

D. Fattal, J. Li, Z. Peng, M. Fiorentino, and R. G. Beausoleil, “Flat dielectric grating reflectors with focusing abilities,” Nature Photonics 4, 466–470 (2010).
[Crossref]

#### Qiu, W.

O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]

#### Rannacher, R.

R. Becker and R. Rannacher, “An optimal control approach to error control and mesh adaption in finite element methods,” Acta Numerica 10, 1–102 (2001).
[Crossref]

#### Reid, M. T. H.

O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]

#### Remlinger, M. J.

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 1,” Journal of Aircraft 36(1), 51–60 (1999).
[Crossref]

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 2,” Journal of Aircraft 36(1), 61–74 (1999).
[Crossref]

#### Reuther, J.

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 1,” Journal of Aircraft 36(1), 51–60 (1999).
[Crossref]

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 2,” Journal of Aircraft 36(1), 61–74 (1999).
[Crossref]

#### Roelkens, G.

G. Roelkens, D. V. Thourhout, and R. Baets, “High efficiency silicon-on-insulator grating coupler based on a poly-silicon overlay,” Optics Express 14, 11622–11630 (2006).
[Crossref] [PubMed]

#### Roundy, D.

A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Computer Physics Communications 181, 687–702 (2010).
[Crossref]

#### Rumelhart, D. E.

D. E. Rumelhart, G. E. Hinton, and R. J. Williams, “Learning internal representations by error propagation,” Parallel Data Processing 1, 318–362 (1986).

#### Saunders, D.

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 1,” Journal of Aircraft 36(1), 51–60 (1999).
[Crossref]

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 2,” Journal of Aircraft 36(1), 61–74 (1999).
[Crossref]

#### seek Chung, Y.

Y. seek Chung, Changyul-Cheon, I.-H. Park, and S.-Y. Hahn, “Optimal shape design of microwave device using fdtd and design sensitivity analysis,” Microwave Theory and Techniques, IEEE Transactions on 48, 2289–2296 (2000).
[Crossref]

#### Seliger, P.

P. Seliger, M. Mahvash, C. Wang, and A. F. J. Levi, “Optimization of aperiodic dielectric structures,” Journal of Applied Physics 100, 034310 (2006).
[Crossref]

#### Soljacic, M.

O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]

#### Taillaert, D.

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### Tam, H. W.

N. K. Nikolova, H. W. Tam, and M. H. Bakr, “Sensitivity analysis with the fdtd method on structured grids,” Microwave Theory and Techniques, IEEE Transactions on 52, 1207–1216 (2004).
[Crossref]

#### Thourhout, D. V.

G. Roelkens, D. V. Thourhout, and R. Baets, “High efficiency silicon-on-insulator grating coupler based on a poly-silicon overlay,” Optics Express 14, 11622–11630 (2006).
[Crossref] [PubMed]

#### van Daele, P.

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### Veronis, G.

G. Veronis, R. W. Dutton, and S. Fan, “Method for sensitivity analysis of photonic crystal devices,” Optics Letters 29, 2288–2290 (2004).
[Crossref] [PubMed]

#### Verstuyft, S.

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### Vuckovic,

J. Lu and Vučković, “Inverse design of nanophotonic structures using complementary convex optimization,” Optics Express 18, 3793–3804 (2010).
[Crossref] [PubMed]

#### Wang, C.

P. Seliger, M. Mahvash, C. Wang, and A. F. J. Levi, “Optimization of aperiodic dielectric structures,” Journal of Applied Physics 100, 034310 (2006).
[Crossref]

#### Williams, R. J.

D. E. Rumelhart, G. E. Hinton, and R. J. Williams, “Learning internal representations by error propagation,” Parallel Data Processing 1, 318–362 (1986).

#### Yablonovitch, E.

C. M. Lalau-Keraly, S. Bhargava, O. D. Miller, and E. Yablonovitch, “Adjoint shape optimization applied to electromagnetic design,” Optics Express 21, 21693–21701 (2013).
[Crossref] [PubMed]

#### Yee, K.

K. Yee, “Numerical solution of initial boundary value problems involving maxwell’s equations in isotropic media,” Antennas and Propagation, IEEE Transactions on 14, 302–307 (1966).
[Crossref]

#### Acta Numerica (2)

R. Becker and R. Rannacher, “An optimal control approach to error control and mesh adaption in finite element methods,” Acta Numerica 10, 1–102 (2001).
[Crossref]

K. Eriksson, D. Estep, P. Hansbo, and C. Johnson, “Introduction to adaptive methods for differential equations,” Acta Numerica 4, 105–158 (1995).
[Crossref]

#### Antennas and Propagation, IEEE Transactions on (1)

K. Yee, “Numerical solution of initial boundary value problems involving maxwell’s equations in isotropic media,” Antennas and Propagation, IEEE Transactions on 14, 302–307 (1966).
[Crossref]

#### Computer Physics Communications (1)

A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: A flexible free-software package for electromagnetic simulations by the FDTD method,” Computer Physics Communications 181, 687–702 (2010).
[Crossref]

#### Journal of Aircraft (2)

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 1,” Journal of Aircraft 36(1), 51–60 (1999).
[Crossref]

J. Reuther, A. Jameson, J. J. Alonso, M. J. Remlinger, and D. Saunders, “Constrained multipoint aerodynamic shape optimization using an adjoint formulation and parallel computers, part 2,” Journal of Aircraft 36(1), 61–74 (1999).
[Crossref]

#### Journal of Applied Physics (1)

P. Seliger, M. Mahvash, C. Wang, and A. F. J. Levi, “Optimization of aperiodic dielectric structures,” Journal of Applied Physics 100, 034310 (2006).
[Crossref]

#### Mathematische Annalen (1)

R. Courant, K. Friedrichs, and H. Lewy, “Über die partiellen Differenzengleichungen der mathematischen Physik,” Mathematische Annalen 100, 32–74 (1928).
[Crossref]

#### Microwave Theory and Techniques, IEEE Transactions on (4)

Y. seek Chung, Changyul-Cheon, I.-H. Park, and S.-Y. Hahn, “Optimal shape design of microwave device using fdtd and design sensitivity analysis,” Microwave Theory and Techniques, IEEE Transactions on 48, 2289–2296 (2000).
[Crossref]

N. Georgieva, S. Glavic, M. Bakr, and J. Bandler, “Feasible adjoint sensitivity technique for em design optimization,” Microwave Theory and Techniques, IEEE Transactions on 50, 2751–2758 (2002).
[Crossref]

N. K. Nikolova, H. W. Tam, and M. H. Bakr, “Sensitivity analysis with the fdtd method on structured grids,” Microwave Theory and Techniques, IEEE Transactions on 52, 1207–1216 (2004).
[Crossref]

N. K. Nikolova, Y. Li, Y. Li, and M. H. Bakr, “Sensitivity analysis of scattering parameters with electromagnetic time-domain simulators,” Microwave Theory and Techniques, IEEE Transactions on 54, 1598–1610 (2006).
[Crossref]

#### Nature Photonics (1)

D. Fattal, J. Li, Z. Peng, M. Fiorentino, and R. G. Beausoleil, “Flat dielectric grating reflectors with focusing abilities,” Nature Photonics 4, 466–470 (2010).
[Crossref]

#### Optics Express (3)

J. Lu and Vučković, “Inverse design of nanophotonic structures using complementary convex optimization,” Optics Express 18, 3793–3804 (2010).
[Crossref] [PubMed]

C. M. Lalau-Keraly, S. Bhargava, O. D. Miller, and E. Yablonovitch, “Adjoint shape optimization applied to electromagnetic design,” Optics Express 21, 21693–21701 (2013).
[Crossref] [PubMed]

G. Roelkens, D. V. Thourhout, and R. Baets, “High efficiency silicon-on-insulator grating coupler based on a poly-silicon overlay,” Optics Express 14, 11622–11630 (2006).
[Crossref] [PubMed]

#### Optics Letters (2)

G. Veronis, R. W. Dutton, and S. Fan, “Method for sensitivity analysis of photonic crystal devices,” Optics Letters 29, 2288–2290 (2004).
[Crossref] [PubMed]

Y. Jiao, S. Fan, and D. A. B. Miller, “Photonic crystal device sensitivity analysis with wannierbasis gradients,” Optics Letters 30, 302–304 (2005).
[Crossref]

#### Parallel Data Processing (1)

D. E. Rumelhart, G. E. Hinton, and R. J. Williams, “Learning internal representations by error propagation,” Parallel Data Processing 1, 318–362 (1986).

#### Physical Review Letters (1)

O. D. Miller, C. W. Hsu, M. T. H. Reid, W. Qiu, B. G. DeLacy, J. D. Joannopoulos, M. Soljačić, and S. G. Johnson, “Fundamental limits to extinction by metallic nanoparticles,” Physical Review Letters 112, 123903 (2014).
[Crossref] [PubMed]

#### Proceedings of the IEEE (1)

V. Liu, D. Miller, and S. Fan, “Highly tailored computational electromagnetics methods for nanophotonic design and discovery,” Proceedings of the IEEE 101, 484–493 (2013).
[Crossref]

#### Quantum Electronics, IEEE Journal of (1)

D. Taillaert, W. Bogaerts, P. Bienstman, T. Krauss, P. van Daele, I. Moerman, S. Verstuyft, K. De Mesel, and R. Baets, “An out-of-plane grating coupler for efficient butt-coupling between compact planar waveguides and single-mode fibers,” Quantum Electronics, IEEE Journal of 38, 949–955 (2002).
[Crossref]

#### Other (1)

J. L. Lions, Optimal Control of Systems Governed by Partial Differential Equations (Springer, 1971).
[Crossref]

### Cited By

OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

### Figures (4)

Fig. 1 The derivative of the longitudinal wave vector β with respect to the thickness of a waveguide slab can faithfully be reproduced using adjoints: the comparison between the direct calculation (blue line) and adjoint calculation (red dots) of $d d h β$ shows that the two methods deliver exactly the same result. However, the adjoint method requires no further expensive solutions of the eigenvalue system in order to determine the derivatives of β with respect to other design variables. Thus, the adjoint method delivers the relevant gradients of a system depending on multiple design parameters faster than the finite-differences approach. The inset shows a schematic of the waveguide slab structure comprised of the cover material nc, the guiding film nf, and the substrate ns.
Fig. 2 Our SiO2 grating coupler redirects free-space plane-waves (vertical arrows) to a guided mode inside the waveguide (horizontal arrows). We show the device at (a) a pitch of 0.58 and a duty cycle of 0.2, and (b) at a pitch of 0.68 and a duty cycle of 0.8. Both the pitch and the duty cycle are optimized to maximize the coupling between the incident light and the light inside the waveguide.
Fig. 3 Adjoint-based optimization convergence: from three initial parameters (stars) in parameter regions of poor optical coupling (blue, green, yellow), adjoint-based optimization evolves the design parameters to the parameter regions of optimal coupling (dark red). The computation of the relevant gradients requires only two Maxwell solutions, independently of the number of design parameters. The finite differences approach involves roughly one Maxwell solution per design variable. Thus, this figure shows that adjoints help optimize the system faster than the finite-differences approach.
Fig. 4 Electric field in and around the waveguide (outlined in light blue) before and after adjoint-based optimization: (a) the grating with the initial design parameters couples the incident light (top) poorly to the waveguide (inside black box), with a complex amplitude of 0.39 (b) After adjoint-based optimization, the grating couples the light very well to the waveguide (inside black box), with a complex amplitude of 0.88.

### Equations (31)

$d T d p k = ∑ i ∂ T ∂ x i d x i d p k + ∂ T ∂ p k = ∂ T ∂ x d x d p k + ∂ T ∂ p k ,$
$d x d p k ≈ 1 Δ p k ( x ( p Δ p k − x ( p ) ) ,$
$R ( x , p ) = 0 .$
$d R d p k = ∂ R ∂ x d x d p k + ∂ R ∂ p k = 0 .$
$d T d p k = ∂ T ∂ x d x d p k + ∂ T ∂ p k − v † ( ∂ R ∂ x d x d p k + ∂ R ∂ p k ) .$
$d T d p k = ∂ T ∂ p k − v † ∂ R ∂ p k + ( ∂ T ∂ x − v † ∂ R ∂ x ) d x d p k .$
$∂ T ∂ x = v † ∂ R ∂ x ,$
$d T d p k = ∂ T ∂ p k − v † ∂ R ∂ p k .$
$∇ × H = ∂ D ∂ t + J σ + J − ∇ × E = ∂ B ∂ t ,$
$B = μ H D = ε E J σ = σ E ,$
$∇ × μ − 1 ∇ × E + ε ∂ 2 E ∂ t 2 + σ ∂ E ∂ t = − ∂ J ∂ t .$
$R ( E ¯ , p ) = ℳ E ¯ − i ω J ¯ = 0 ,$
$ℳ E ¯ = ∇ × μ − 1 ∇ × E ¯ − ω 2 ε E ¯ .$
$T ( E ¯ , p ) = 1 N ∫ Ω d Ω | E ¯ | 2 = E ¯ 0 † E ¯ ,$
$ℳ v ¯ − i ω ( 1 i ω E ¯ 0 ) = 0 .$
$d T d p k = − v † ∂ ε ∂ p k E ¯ .$
$tan ( h κ f ) = γ c + γ s κ f ( 1 − γ c γ s κ f 2 ) ,$
$R = tan ( h κ f ) − γ c + γ s κ f ( 1 − γ c γ s κ f 2 ) ,$
$T = β ,$
$∂ T ∂ β = 1 = v ∂ R ∂ β .$
$d T d p k = d β d p k = − v ∂ R ∂ p k .$
$∂ R ∂ β = ∂ ∂ β ( tan ( h κ f ) ) − ∂ ∂ β ( γ c + γ s κ f ( 1 − γ c γ s κ f 2 ) ) ,$
$∂ ∂ β ( tan ( h κ f ) ) = − h β κ f cos 2 ( h κ f ) ,$
$∂ ∂ β ( γ c + γ s q ) = β γ c + β γ s q − γ c + γ s q 2 ∂ q ∂ β ,$
$∂ q ∂ β = − β κ f [ 1 + γ s γ c + γ c γ s + γ c γ s κ f 2 ] .$
$1 v = − β [ h κ f cos 2 ( h κ f ) + ( γ c + γ s ) κ f κ f 2 − γ c γ s ( 1 γ s γ c + 1 + γ s γ c + γ c γ s + γ c γ s κ f 2 κ f 2 − γ c γ s ) ] ,$
$v = 1 ∂ R ∂ β$
$d β d p k = − v ∂ R ∂ p k .$
$d R d h = κ f cos 2 ( h κ f )$
$d β d h = − v κ f cos 2 ( h κ f ) .$
$p = ( p d )$