HELIOS: A surface integral equation software for light scattering in homogeneous, periodic, and stratified environments

This paper introduces HELIOS, an open-source C++ and Python software that utilizes the PMCHWT surface integral equation formulation with RWG basis functions to accurately model light scattering by particles in homogeneous, stratified, and periodic environments.

The Gist

Imagine you are trying to predict how a drop of water ripples when it hits a complex rock in a pond. Now, imagine that rock isn't just a rock—it's a tiny, intricate sculpture made of gold, sitting in a pond that might be empty, might be filled with layers of different liquids, or might be part of an infinite grid of identical ponds stretching forever.

This is the challenge of nanophotonics: figuring out how light (which acts like a wave) behaves when it hits tiny structures. Doing this by hand is impossible because the math is incredibly complex. That's where HELIOS comes in.

Think of HELIOS as a super-smart, digital crystal ball built by scientists at EPFL in Switzerland. It's a piece of free software that lets researchers simulate how light interacts with these tiny nano-worlds without needing to build them in a lab first.

Here is how it works, broken down into simple concepts:

1. The Core Idea: Painting the Surface, Not the Whole Room

Most computer programs try to simulate light by filling the entire 3D space with a grid of tiny boxes (like pixelating a whole room). This is slow and uses up a lot of computer power.

HELIOS is smarter. It uses a technique called Surface Integral Equations.

• The Analogy: Imagine you want to know how sound bounces off a statue. Instead of measuring the air inside the statue and the air outside it, HELIOS only looks at the surface of the statue. It treats the surface like a canvas and paints "currents" (imaginary flows of electricity and magnetism) onto it.
• Why it's great: By only focusing on the skin of the object, it saves massive amounts of computing power, allowing it to handle very complex shapes quickly.

2. The Three Worlds HELIOS Can Visit

HELIOS is unique because it can handle three different "worlds" or environments where these tiny objects might live:

• World A: The Empty Room (Homogeneous)

  • Scenario: A single gold nanoparticle floating in a vacuum or water.
  • How HELIOS does it: It uses a proven mathematical recipe (called PMCHWT) to solve the puzzle. It's like solving a standard jigsaw puzzle where the pieces fit together perfectly.

• World B: The Infinite Hall of Mirrors (Periodic)

  • Scenario: A photonic crystal or a metasurface, which is essentially a repeating pattern of nano-objects (like a honeycomb made of light).
  • The Problem: If you try to simulate an infinite pattern, the math goes on forever.
  • The HELIOS Trick: It uses something called Ewald's Transformation.
  • The Analogy: Imagine trying to count the echoes in a room with infinite mirrors. Instead of counting every single echo one by one, HELIOS uses a mathematical "shortcut" (a transformation) to instantly calculate the sum of all those infinite reflections. It turns a never-ending task into a quick calculation.

• World C: The Layered Cake (Stratified)

  • Scenario: A nano-object sitting on top of a silicon chip, which is sitting on a glass slide, which is sitting on a metal base.
  • The Problem: Light bounces back and forth between every layer, creating a messy web of reflections.
  • The HELIOS Trick: It uses a "Matrix-Friendly" approach and a Tabulation-Interpolation scheme.
  • The Analogy: Instead of calculating the bounces in real-time every single time, HELIOS creates a giant "cheat sheet" (a table) of how light behaves in these layers. When it needs an answer, it looks it up on the sheet and fills in the gaps. This makes it incredibly fast, even when the layers are complex.

3. How You Use It (The User Experience)

You don't need to be a math wizard to use HELIOS.

• The Engine (C++): The heavy lifting is done by a powerful C++ engine. Think of this as the high-performance race car engine under the hood.
• The Steering Wheel (Python): The user interacts with a Python interface. This is the dashboard. You tell the car where to go (set up the simulation), press the gas (run the solver), and look at the speedometer (view the results).
• The Output: It gives you beautiful pictures of where the light is bright (hotspots) or dark, and graphs showing how much light is absorbed or scattered.

4. Why Does This Matter?

Why do we need a digital crystal ball for tiny things?

• Designing Better Tech: Engineers can use HELIOS to design better solar cells, faster computer chips, or super-sensitive medical sensors before they ever spend money building them.
• Saving Time: It prevents scientists from building a prototype, realizing it doesn't work, and starting over. They can simulate it, fix the design on the computer, and then build the perfect version.

In a Nutshell

HELIOS is an open-source tool that acts as a virtual wind tunnel for light. Whether you have a single particle, an infinite grid of them, or a sandwich of layers, HELIOS uses clever math shortcuts to predict exactly how light will dance around them. It turns impossible math problems into manageable, fast simulations, helping scientists build the future of light-based technology.

Technical Summary

1. Problem Statement

Computational electromagnetics is essential for designing nanophotonic devices, ranging from isolated nanoparticles to complex periodic structures (metasurfaces, photonic crystals) and devices embedded in multilayered substrates. While Surface Integral Equation (SIE) methods are highly efficient for 3D scatterers because they only require meshing boundaries (reducing dimensionality) and inherently satisfy radiation conditions, robust open-source software capable of handling all these environments simultaneously is lacking.

Specific challenges addressed by this work include:

• Penetrable Objects: Solving scattering for objects with arbitrary 3D geometries and material properties in homogeneous backgrounds.
• Periodic Structures: Handling the slow convergence of periodic dyadic Green's tensors required for infinite 2D lattices.
• Stratified Media: Accurately modeling scatterers embedded in multilayered environments, which requires evaluating complex Sommerfeld integrals and handling singularities at layer interfaces.

2. Methodology

HELIOS is an open-source software package implemented in C++ (for the computational core) with a Python interface for workflow management. It utilizes the Poggio-Miller-Chang-Harrington-Wu-Tsai (PMCHWT) formulation, a standard and stable approach for penetrable scatterers.

Key Technical Components:

• Discretization: Domain boundaries are discretized using triangular meshes. Unknown electric and magnetic surface current densities are expanded using Rao-Wilton-Glisson (RWG) basis functions. The system is solved using the Galerkin testing procedure.
• Singularity Handling: Numerical integration singularities in the Green's function are managed via a singularity subtraction technique, decomposing the function into smooth and singular components (evaluated semi-analytically).
• Periodic Environments (Mode 1):
  • Implements Floquet-periodic boundary conditions.
  • Uses Ewald's transformation to accelerate the evaluation of the infinite series associated with the periodic Green's tensor, splitting it into rapidly converging spatial and reciprocal space sums.
  • Utilizes a precomputed binary lookup table (erfc.bin) for the complementary error function to speed up calculations.
• Stratified Media (Mode 2):
  • Employs a matrix-friendly formulation for the layered media Green's tensor, separating it into primary (homogeneous) and secondary (layered interaction) terms.
  • Handles Sommerfeld integrals using optimized integration paths (semi-elliptical or parallel to the imaginary axis) based on the vertical distance between source and observation points.
  • Implements a tabulation-interpolation scheme to accelerate repeated Green's function evaluations. The integrand is decomposed into quasi-static (analytically calculated) and smooth (tabulated) parts to ensure stability and speed.
  • Includes specific line integral corrections when mesh triangles straddle layer interfaces to ensure the continuity of tangential electromagnetic fields.

3. Key Contributions

• Unified Framework: HELIOS is the first open-source SIE solver to integrate capabilities for homogeneous, periodic, and stratified environments within a single codebase.
• Modular Architecture: The software separates the heavy C++ computational core from the Python workflow, allowing for easy simulation setup, execution, and post-processing.
• Advanced Acceleration:
  • Ewald summation for periodic structures.
  • Tabulation-interpolation and quasi-static decomposition for stratified media Green's functions.
• Comprehensive Workflow: Provides tools for mesh generation (via Gmsh/COMSOL), simulation configuration, solving, and visualization of cross-sections and near-field distributions.
• Open Source: Released under the GNU GPL v3, with source code and documentation available on GitHub.

4. Results and Validation

The authors validated HELIOS through a series of benchmark examples covering all three operational modes:

• Homogeneous Background (Mode 0):
  • Gold Sphere: Results for scattering, absorption, and extinction cross-sections matched perfectly with analytical Mie theory.
  • Janus Ring: Demonstrated the ability to model inhomogeneous particles with azimuthally varying material properties (Silicon and Gold).
• Periodic Structures (Mode 1):
  • Photonic Crystal: Simulated an infinite array of dielectric pillars, successfully extracting reflectance and transmittance spectra and visualizing near-field distributions.
  • Fishnet Structure: Modeled a continuous metallic fishnet, demonstrating the solver's ability to handle complex unit cells and extract spectral features.
• Stratified Media (Mode 2):
  • Sphere on Substrate: Validated the red-shift and resonance modification of a silver sphere on a SiO2 substrate compared to free space.
  • Triangular Hole in Gold Film: Verified field confinement at aperture vertices and correct handling of bodies located entirely within intermediate layers.
  • Cylindrical Hole in Multilayer Stack: Tested the robustness of the line integral corrections and tabulation scheme for a hole penetrating multiple interfaces (SiO2, Si, Au, Air). Field continuity across interfaces was explicitly verified, confirming the accuracy of the implementation.

5. Significance

HELIOS fills a critical gap in the nanophotonics research community by providing a reliable, versatile, and open-source tool for rigorous light scattering simulations. Its significance lies in:

1. Accessibility: It lowers the barrier to entry for researchers needing high-accuracy SIE solvers for complex environments without relying on proprietary commercial software.
2. Versatility: It enables the study of realistic device configurations (e.g., metasurfaces on substrates, nanoparticles in biological layers) that were previously difficult to model with a single tool.
3. Efficiency: The implementation of advanced acceleration techniques (Ewald, tabulation-interpolation) makes computationally expensive problems (like multilayered periodic structures) feasible on standard desktop hardware.
4. Community Growth: Its open-source nature encourages collaborative development, allowing for future extensions and community-driven improvements.