Accelerated 3D Maxwell Integral Equation Solver using the Interpolated Factored Green Function Method

Jagabandhu Paul,Constantine Sideris

2024-08-05

Abstract:This article presents an $O(N\log N)$ method for numerical solution of Maxwell's equations for dielectric scatterers using a 3D boundary integral equation (BIE) method. The underlying BIE method used is based on a hybrid Nyström-collocation method using Chebyshev polynomials. It is well known that such an approach produces a dense linear system, which requires $O(N^2)$ operations in each step of an iterative solver. In this work, we propose an approach using the recently introduced Interpolated Factored Green's Function (IFGF) acceleration strategy to reduce the cost of each iteration to $O(N\log N)$. To the best of our knowledge, this paper presents the first ever application of the IFGF method to fully-vectorial 3D Maxwell problems. The Chebyshev-based integral solver and IFGF method are first introduced, followed by the extension of the scalar IFGF to the vectorial Maxwell case. Several examples are presented verifying the $O(N\log N)$ computational complexity of the approach, including scattering from spheres, complex CAD models, and nanophotonic waveguiding devices. In one particular example with more than 6 million unknowns, the accelerated IFGF solver runs 42x faster than the unaccelerated method.

Numerical Analysis,Computational Physics

What problem does this paper attempt to address?

### What problem does this paper attempt to solve?

This paper aims to propose an efficient numerical method with a computational complexity of \(O(N\log N)\) for solving the three - dimensional Maxwell equations of dielectric scatterers. Specifically, the author uses the Boundary Integral Equation (BIE) method combined with the Interpolated Factored Green's Function (IFGF) acceleration strategy to reduce the computational cost required for each iteration.

#### Main challenges and solutions

1. **High computational complexity of dense linear systems**:
- Using the traditional hybrid Nyström - collocation method (based on Chebyshev polynomials) will result in a dense linear system, making the computational complexity of each iteration reach \(O(N^2)\).
- To reduce the computational complexity, the author introduced the IFGF acceleration strategy, reducing the computational complexity of each iteration from \(O(N^2)\) to \(O(N\log N)\).

2. **Limitations of existing acceleration methods**:
- Many existing acceleration algorithms rely on the Fast Fourier Transform (FFT), which has parallelization challenges in distributed - memory parallel computer architectures.
- The IFGF method avoids these problems by recursive interpolation (using low - order Chebyshev expansions) without relying on FFT.

3. **Application areas**:
- This method is particularly suitable for electromagnetic scattering problems, such as the design and optimization of antennas, microwaves, and photonic devices. It can significantly improve the running time and accuracy, especially in inverse design methods that require accurate field and gradient information.

#### Experimental verification

- The paper shows multiple experimental results, verifying the \(O(N\log N)\) computational complexity of this method.
- In an example with more than 6 million unknowns, the accelerated IFGF solver is 42 times faster than the unaccelerated method.

### Summary

The main contribution of this paper is that it applies the IFGF method to the full - vector 3D Maxwell problem for the first time, achieving efficient electromagnetic scattering calculations, solving the problem of excessively high computational complexity in traditional methods, and showing significant performance improvements in practical applications.