Resonance frequency and radiative Q-factor of plasmonic and dielectric modes of small objects

The electromagnetic scattering resonances of a non-magnetic object much smaller than the incident wavelength in vacuum can be either described by the electroquasistatic approximation of the Maxwell's equations if its permittivity is negative, or by the magnetoquasistatic approximation if its permittivity is positive and sufficiently high. Nevertheless, these two approximations fail to correctly account for the frequency shift and the radiative broadening of the resonances when the size of the object becomes comparable to the wavelength of operation. In this manuscript, the radiation corrections to the electroquasistatic and magnetoquasistatic resonances of arbitrarily-shaped objects are derived, which only depend on the quasistatic current modes. Then, closed form expressions of the frequency-shift and the radiative Q-factor of both plasmonic and dielectric modes of small objects are introduced, where the dependencies on the material and the size of the object are factorized. In particular, it is shown that the radiative Q-factor explicitly depends on the multipolar components of the quasistatic modes

Magnetoquasistatic resonances of small dielectric objects

A small dielectric object with positive permittivity may resonate when the free-space wavelength is large in comparison with the object dimensions if the permittivity is sufficiently high. We show that these resonances are described by the magnetoquasistatic approximation of the Maxwell's equations in which the normal component of the displacement current density field vanishes on the surface of the particle. They are associated to values of permittivities and frequencies for which source-free quasistatic magnetic fields exist, which are connected to the eigenvalues of a magnetostatic integral operator. We present the general physical properties of magnetoquasistatic resonances in dielectrics with arbitrary shape. They arise from the interplay between the polarization energy stored in the dielectric and the energy stored in the magnetic field. Our findings improve the understanding of resonances in high-permittivity dielectric objects and provide a powerful tool that greatly simplifies the analysis and design of high-index resonators

Volume Integral Formulation for the Calculation of Material Independent Modes of Dielectric Scatterers

In the frame of volume integral equation methods, we introduce an alternative representation of the electromagnetic field scattered by a homogeneous object of arbitrary shape at a given frequency, in terms of a set of modes independent of its permittivity. This is accomplished by introducing an auxiliary eigenvalue problem, based on a volume integral operator. With this modal basis the expansion coefficients of the scattered field are simple rational functions of the permittivity of the scatterer. We show, by studying the electromagnetic scattering from a sphere and a cylinder of dimensions comparable to the incident wavelength, that only a moderate number of modes is needed to accurately describe the scattered far field. This method can be used to investigate resonant scattering phenomena, including plasmonic and photonic resonances, and to design the permittivity of the object to pursue a prescribed tailoring of the scattered field. Moreover, the presented modal expansion is computationally advantageous compared to direct solution of the volume integral equation when the scattered field has to be computed for many different values of the dielectric permittivity, given the size and shape of the dielectric body

Synthesis of resonant modes in electromagnetics

Resonant modes determine the response of electromagnetic devices, including dielectric and plasmonic resonators. Relying on the degrees of freedom that metamaterials provide, this contribution shows how to design, at will, the resonant modes of a dielectric object placed in an unbounded space. Specifically, the proposed method returns in analytical form the spatial distribution of the dielectric susceptibility tensor for which the object exhibits resonances at prescribed frequencies and spatial distribution of the polarization. Together with the synthesis of the material, two key concepts are introduced: the controlled tunability of the resonant modes and the number of essential modes, i.e. the number of modes that uniquely characterize the spatial distribution of the dielectric susceptibility. Moreover, this approach can be applied to design the resonant modes of any system where the constitutive relationship is linear and local

A Fast Matrix Compression Method for Large Scale Numerical Modelling of Rotationally Symmetric 3D Passive Structures in Fusion Devices

This paper illustrates the development of a recursive QR technique for the analysis of transient events, such as disruptions or scenario evolution, in fusion devices with three-dimensional conducting structures using an integral eddy current formulation. An integral formulation involves the solution, at each time step, of a large full linear system. For this reason, a direct solution is impractical in terms of time and memory consumption. Moreover, typical fusion devices show a symmetric/periodic structure. This can be properly exploited when the plasma and other sources possess the same symmetry/periodicity of the structure. Indeed, in this case, the computation can be reduced to only a single sector of the overall structure. In this work the periodicity and the symmetries are merged in the recursive QR technique, exhibiting a huge decrease in the computational cost. Finally, the proposed technique is applied to a realistic large-scale problem related to the International Thermonuclear Experimental Reactor (ITER)

Implementation of matrix compression in the coupling of JOREK to realistic 3D conducting wall structures

JOREK is an advanced non-linear simulation code for studying MHD instabilities in magnetically confined fusion plasmas and their control and/or mitigation. A free-boundary and resistive wall extension was introduced via coupling to the STARWALL and CARIDDI codes, both able to provide dense response matrices describing the electromagnetic interactions between plasma and conducting structures. For detailed CAD representations of the conducting structures and high resolutions for the plasma region, memory and computing time limitations restrict the possibility of simulating the ITER tokamak. In the present work, the Singular Value Decomposition provided by routines from the ScaLAPACK library has been successfully applied to compress some of the dense response matrices and thus optimize memory usage. This is demonstrated for simulations of Tearing Mode and Vertical Displacement Event instabilities. An outlook to future applications on large production cases and further extensions of the method are discussed.Comment: 19 pages, 9 figures, firstly presented at the 4th Fusion HPC Workshop in November 2023, submitted to Plasma Physics and Controlled Fusion in April 2024, First Revision after referee report in July 202