The paper presents a derivation of analytical components of S-matrices for arbitrary planar diffractive structures and metasurfaces in the Fourier domain. Attained general formulas for S-matrix components can be applied within both formulations in the Cartesian and curvilinear metric. A numerical method based on these results can benefit from all previous improvements of the Fourier domain methods. In addition, we provide expressions for S-matrix calculation in case of periodically corrugated layers of 2D materials, which are valid for arbitrary corrugation depth-to-period ratios. As an example the derived equations are used to simulate resonant grating excitation of graphene plasmons and an impact of silica interlayer on corresponding reflection curves.
We present a method of incorporating the discrete dipole approximation (DDA) method with the point matching method to formulate the T-matrix for modeling arbitrarily shaped micro-sized objects. The emph{T}-matrix elements are calculated using point matching between fields calculated using vector spherical wave functions and DDA. When applied to microrotors, their discrete rotational and mirror symmetries can be exploited to reduce memory usage and calculation time by orders of magnitude; a number of optimization methods can be employed based on the knowledge of the relationship between the azimuthal mode and phase at each discrete rotational point, and mode redundancy from Floquets theorem. A reduced-mode T-matrix can also be calculated if the illumination conditions are known.
A numerical implementation of the transition state theory (TST) is presented which can be used to calculate the attempt frequency $f_{0}$ of arbitrary shaped magnetic nanostructures. The micromagnetic equations are discretized using the finite element method. The climbing image nudged elastic band method is used to calculate the saddle point configuration, which is required for the calculation of $f_{0}$. Excellent agreement of the implemented numerical model and analytical solutions is obtained for single domain particles. The developed method is applied to compare $f_{0}$ for single phase and graded media grains of advanced recording media. $f_{0}$ is predicted to be comparable if the maximum anisotropy is the same in these two media types.
Inverse design of large-area metasurfaces can potentially exploit the full parameter space that such devices offer and achieve highly efficient multifunctional flat optical elements. However, since practically useful flat optics elements are large in the linear dimension, an accurate simulation of their scattering properties is challenging. Here, we demonstrate a method to compute accurate simulations and gradients of large-area metasurfaces. Our approach relies on two key ingredients - a simulation distribution strategy that allows a linear reduction in the simulation time with number of compute (GPU) nodes and an efficient single-node computation using the Transition-matrix (T-matrix) method. We demonstrate ability to perform a distributed simulation of large-area, while accurately accounting for scatterer-scatterer interactions significantly beyond the locally periodic approximation, and efficiently compute gradients with respect to the metasurface design parameters. This scalable and accurate metasurface simulation method opens the door to gradient-based optimization of full large-area metasurfaces.
The classical adjoint-based topology optimization (TO) method, based on the use of a random continuous dielectric function as an adjoint variable distribution, is known to be one of the most efficient optimization methods that enable the design of optical devices with outstanding performances. However, the strategy for selecting the optimal solution requires a very fine pixelation of the permittivity function of the profile under optimization. Typically, at least 28 pixels are needed while optimizing a one wavelength wide 1D metagrating. This makes it very difficult to extend TO methods to large-scale optimization problems. In this paper, we introduce a new concept of adjoint-based topology optimization that enables fast and efficient geometry based design of both periodic and aperiodic metasurfaces. The structures are built from nano-rods whose widths and positions are to be adjusted. Our new approach requires a very low number of design parameters, thus leading to a drastic reduction in the computational time: about an order of magnitude. Hence, this concept makes it possible to address the optimization of large-scale structures in record time. As a proof-of-concept we apply this method to the design of (i) a periodic metagrating, optimized to have a specific response into a particular direction, and (ii) a dielectric metalens (aperiodic metasurface), enabling a high energy focusing into a well-defined focal spot.
We investigate the use of a Genetic Algorithm (GA) to design a set of photonic crystals (PCs) in one and two dimensions. Our flexible design methodology allows us to optimize PC structures which are optimized for specific objectives. In this paper, we report the results of several such GA-based PC optimizations. We show that the GA performs well even in very complex design spaces, and therefore has great potential for use as a robust design tool in present and future applications.
Alexey A. Shcherbakov
,Yury V. Stebunov
,Denis F. Baidin andn Thomas Kampfe
.
(2017)
.
"Direct S-matrix calculation for diffractive structures and metasurfaces"
.
Alexey Shcherbakov A.
هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا