Smuthi stands for ‘scattering by multiple particles in thin-film systems’. It is a Python software that allows to solve light scattering problems involving one ore multiple particles near or inside a system of planar layer interfaces.
It solves the Maxwell equations (3D wave optics) in frequency domain (one wavelength per simulation).
Smuthi is based on the T-matrix method for the single particle scattering and on the scattering-matrix method for the propagation through the layered medium. This combination of methods is described in the following papers:
- Egel, Amos, and Uli Lemmer. JQSRT 148 (2014): 165-176.
- Egel, Amos, Siegfried W. Kettlitz, and Uli Lemmer. JOSA A 33.4 (2016): 698-706.
- Egel, Amos: PhD thesis (in preparation)
For spheres, the T-matrix of the individual particles is given by the Mie-coefficients. For spheroids and finite cylinders, Smuthi calls the NFM-DS to compute the single particle T-matrix. This is a Fortran software package written by A. Doicu, T. Wriedt and Y. Eremin, based on the “Null-field method with discrete sources”, see
Performance critical parts of the software are implemented in CUDA. When dealing with a large number of particles, Smuthi can benefit from a substantial acceleration if a suitable (NVIDIA) GPU is available.
For CPU-only execution, other acceleration concepts (including MPI parallelization, Numba JIT compilation) are currently tested.
Range of applications¶
Smuthi can be applied to any scattering problem in frequency domain involving
- a system of plane parallel layer interfaces separating an arbitrary number of metallic or dielectric layers.
- an arbitrary number of wavelength-scale scattering particles (currently available: spheres, spheroids, finite cylinders). The particles can be metallic or dielectric and rotated to an arbitrary orientation.
- an initial field in form of a plane wave, a beam (currently available: beam with Gaussian xy-profile) or a collection of dipole sources
Thus, the range of applications spans from scattering by a single particle on a substrate to scattering by several thousand particles inside a planarly layered medium. For a number of examplary simulations, see the examples section.
Smuthi can compute
- the 3D electric field, for example along a cut plane and save it in the form of ascii data files, png images or gif animations.
- the far field power flux of the total field, the initial field or the scattered field. For plane wave excitation, it can be processed to the form of differential scattering and extinction cross sections.
- For dipole sources, the dissipated power can be computed (Purcell effect).
The following issues need to be considered when applying Smuthi:
- Particles must not intersect with each other or with layer interfaces.
- Magnetic or anisotropic materials are currently not supported.
- The method is in principle valid for a wide range of particle sizes - however, the numerical validity has only been tested for particle diameters up to around one wavelength. For larger particles, note that the number of multipole terms in the spherical wave expansion grows with the particle size. For further details, see the hints for the selection of the multipole truncation order.
- Smuthi does not provide error checking of user input, nor does it check if numerical parameters specified by the user are sufficient for accurate simulation results. It is thus required that the user develops some understanding of the influence of various numerical parameters on the validity of the results. See the simulation guidelines.
- A consequence of using the T-matrix method is that the electric field inside the circumscribing sphere of a particle cannot be correctly computed, see for example Auguié et al. (2016). In the electric field plots, the circumscribing sphere is displayed as a dashed circle around the particle as a reminder that there, the computed near fields cannot be trusted.
- Particles with initersecting circumscribing spheres can lead to incorrect results. The use of Smuthi is therefore limited to geometries with particles that have disjoint circumscribing spheres.
- If particles are located near interfaces, such that the circumscribing shere of the particle intersects the interface, a correct simulation result can in principle be achieved. However, special care has to be taken regarding the selection of the truncation of the spherical and plane wave expansion, see the hints for the selection of the wavenumber truncation.
Smuthi was written and is maintained by Amos Egel. Email to for questions, feature requests or if you would like to contribute.
Smuthi includes contributions from the following persons:
- Adrian Doicu, Thomas Wriedt and Yuri Eremin through the NFM-DS package, a copy of which is distributed with Smuthi.
- Dominik Theobald implemented functions for the simulation of particles with arbitrary orientation. He currently works on the implementation of a plane wave based particle coupling for non-spherical particles with overlapping circumscribing spheres.
- Konstantin Ladutenko with many useful additions, including example simulations, smoother input/output, support of MPI parallel computing (currently in construction) and the option of permanent NFM-DS installation folder.
Big thanks go to Lorenzo Pattelli for designing the Smuthi logo.
Ilia Rasskazov has helped with useful comments and bug reports.