Abstract

Recently, the volume integral equation (VIE) approach has been proposed as an efficient simulation tool for silicon photonics applications [J. Lightw. Technol. 36, 3765 (2018) [CrossRef]  ]. However, for the high-frequency and strong-contrast problems arising in photonics, the convergence of iterative solvers for the solution of the linear system can be extremely slow. The uniform discretization of the volume integral operator leads to a three-level Toeplitz matrix, which is well suited to preconditioning via its circulant approximation. In this paper, we describe an effective circulant preconditioning strategy based on the multilevel circulant preconditioner of Chan and Olkin [Numer. Algorithms 6, 89 (1994) [CrossRef]  ]. We show that this approach proves ideal in the canonical photonics problem of propagation within a uniform waveguide, in which the flow is unidirectional. For more complex photonics structures, such as Bragg gratings, directional couplers, and disk resonators, we generalize our preconditioning strategy via geometrical partitioning (leading to a block-diagonal circulant preconditioner) and homogenization (for inhomogeneous structures). Finally, we introduce a novel memory reduction technique enabling the preconditioner’s memory footprint to remain manageable, even for extremely long structures. The range of numerical results we present demonstrates that the preconditioned VIE is fast and has great utility for the numerical exploration of prototype photonics devices.

© 2019 Optical Society of America

Direct solution of the volume integral equation for electromagnetic scattering from dielectric objects with negative permittivity

A-Li Deng and Li-Ming Zhang
J. Opt. Soc. Am. A 38(7) 1015-1023 (2021)

Multigrid preconditioned conjugate-gradient method for large-scale wave-front reconstruction

Luc Gilles, Curtis R. Vogel, and Brent L. Ellerbroek
J. Opt. Soc. Am. A 19(9) 1817-1822 (2002)

Decreasing the memory of the discontinuous Galerkin volume integral equation method for scattering from inhomogeneous dielectric objects

A-Li Deng, Qi-Gang Zhu, Bing-Zhong Ren, and Li-Ming Zhang
J. Opt. Soc. Am. A 40(9) 1654-1661 (2023)

A limited-memory, quasi-Newton preconditioner for nonnegatively constrained image reconstruction

Johnathan M. Bardsley
J. Opt. Soc. Am. A 21(5) 724-731 (2004)

Multilevel Green's function interpolation method for analysis of 3-D frequency selective structures using volume/surface integral equation

Yan Shi and Chi Hou Chan
J. Opt. Soc. Am. A 27(2) 308-318 (2010)