Nonlinear optical materials, such as transparent conductive oxides, have recently drawn a lot of attention when being integrated into metasurfaces and allowing full-optical control of the surface response. Although several methods for modeling the nonlinear materials have been proposed in the literature, most of them have the limitations on being non-dispersive and of instantaneous response. In this paper, we present a straightforward integration of an extended Drude-Lorentz model that captures the local intensity response of nonlinear materials while being dispersive and allowing for inertial response via a low-pass filtering process. This method is integrated into standard finite difference time-domain (FDTD) implementation of Maxwell’s equations and the auxiliary differential equations approach of the Drude-Lorentz model is extended via local intensity-dependent parameters. A numerical demonstration shows the response for a thin film of nonlinear material, where the parameters across the sample are time-varying with respect to the local intensity of the fields. Therefore, showing a direct feedback of the field profile to the nonlinear response of the material, which is critical when incorporating such films in resonating meta-atoms.

A simple concept to perfectly decouple the reflection magnitude and phase of a reflective metasurface is proposed. For a linearly-polarized illumination, a full 2? phase is achieved in reflection by physically displacing the unit cell normal to the surface, with respect to a fixed reference and source, while the reflection magnitude is independently controlled via an integrated resistive element inside the unit cell. This proposed concept is demonstrated using full-wave simulations, which is the first step towards an electromechanical metasurface with truly independent real-time amplitude and phase control.

We present a rigorous semi-analytical formulation for the analysis of electromagnetic (EM) scattering from a cylinder constructed from a metasurface represented using surface susceptibilities. The formulation uses the Generalized Sheet Transition Conditions (GSTCs) to represent the surface and eigenfunction expansion (EFE) of the incident and scattered fields, by exploiting their angular periodicity. Incorporating a completely general non-uniform surface formulation of 36 susceptibility components, a matrix equation is formulated that can be solved for the field harmonic coefficients. The paper illustrates the methodology with a number of examples; including two formed from a finite sized practical unit-cell exhibiting a normally oriented magnetic resonance which is compared with commercial EM full-wave solver; other examples present surfaces that have a modulated gain/loss profile (i.e. amplitude modulation) and polarization conversion. It is found that for all cases the EFE solution very accurately captures the scattered fields in both the interior and exterior regions of the surface. Detailed convergence studies are further presented including the effect of susceptibility modulation on the number of terms needed in the EFEs, to reach a correct field solution. The proposed EFE framework is further compared with a second GSTC based simulator using an Integral Equations (IE) approach and it is found that the IE-GSTC simulated fields approach that of the EFE as the surface discretization is increased. The proposed EFE approach is thus a quick, rigorous methodology which, although limited in the geometry which it can model, has many advantages for investigation into the use of GSTCs, application development and providing a baseline for simulation studies.

A ray optical methodology based on the uniform theory of diffraction is proposed to model electromagnetic field scattering from curved metasurfaces. The problem addressed is the illumination of a purely reflective uniform cylindrical metasurface by a line source, models the surface with susceptibilities and employs a methodology previously used for cylinders coated in thin dielectric layers [1]. The approach is fundamentally based on a representation of the metasurface using the General Sheet Transition Conditions (GSTCs) which characterizes the surface in terms of susceptibility dyadics. An eigenfunction description of the metasurface problem is derived considering both tangential and normal surface susceptibilities, and used to develop a ray optics (RO) description of the scattered fields; including the specular geometrical optical field, surface diffraction described by creeping waves and a transition region over the shadow boundary. The specification of the fields in the transition region is dependent on the evaluation of the Pekeris caret function integral and the method follows [1]. The proposed RO-GSTC model is then successfully demonstrated for a variety of cases and is independently verified using a rigorous eigenfunction solution (EF-GSTC) and full-wave Integral Equation method (IE-GSTC), over the entire domain from the deep lit to deep shadow.

A static metasurface reflector based on a novel coupled resonator configuration is proposed to independently control the reflection phase and magnitude of linearly polarized incident fields, and is demonstrated experimentally in the millimeter-wave Ka-band around 30 GHz. The proposed concept is illustrated using a unit cell design consisting of a rectangular ring coupled with a rectangular slot resonator backed by a grounded dielectric slab. By geometrically tuning various dimensions of the two resonators, a near-perfect amplitude-phase coverage is achieved at a fixed design frequency of 30 GHz. To demonstrate the flexible beam-forming capability of the proposed metasurface reflectors, illustrative examples of fixed beam steering with varying reflection magnitudes, and asymmetric dual-beam patterns with specified reflection magnitude, reflection angles and beam-widths, are successfully shown. Compared to the standard method based on polarization rotation and resistive loadings with discrete values, the proposed technique does not generate undesired cross-polarization field reflection, and provides a continuous magnitude tuning including full absorption, along with wide phase coverage.

While metasurfaces (MSs) are constructed from deeply-subwavelength unit cells, they are generally electrically-large and full-wave simulations of the complete structure are computationally expensive. Thus, to reduce this high computational cost, non-uniform MSs can be modeled as zero-thickness boundaries, with sheets of electric and magnetic polarizations related to the fields by surface susceptibilities and the generalized sheet transition conditions (GSTCs). While these two-sided boundary conditions have been extensively studied for single sheets of resonant particles, it has not been shown if they can correctly model structures where the two sides are electrically isolated, such as a fully-reflective surface. In particular, we consider in this work whether the fields scattered from a fully reflective metasurface can be correctly predicted for arbitrary field illuminations, with the source placed on either side of the surface. In the process, we also show the mapping of a PEC sheet with a dielectric cover layer to bi-anisotropic susceptibilities. Finally, we demonstrate the use of the susceptibilities as compact models for use in various simulation techniques, with an illustrative example of a parabolic reflector, for which the scattered fields are correctly computed using a integral equation (IE) based solver.

An Integral Equation (IE) based field solver to compute the scattered fields from spatially dispersive metasurfaces is proposed and numerically confirmed using various examples involving physical unit cells. The work is a continuation of Part- 1 [1], which proposed the basic methodology of representing spatially dispersive metasurface structure in the spatial frequency domain, k. By representing the angular dependence of the surface susceptibilities in k as a ratio of two polynomials, the standard Generalized Sheet Transition Conditions (GSTCs) have been extended to include the spatial derivatives of both the difference and average fields around the metasurface. These extended boundary conditions are successfully integrated here into a standard IE-GSTC solver, which leads to the new IEGSTC-SD simulation framework presented here. The proposed IE-GSTC-SD platform is applied to various uniform metasurfaces, including a practical short conducting wire unit cell, as a representative practical example, for various cases of finite-sized flat and curvilinear surfaces. In all cases, computed field distributions are successfully validated, either against the semi-analytical Fourier decomposition method or the brute-force full-wave simulation of volumetric metasurfaces in the commercial Ansys FEM-HFSS simulator.

A simple method to describe spatially dispersive metasurfaces is proposed where the angle-dependent surface susceptibilities are explicitly used to formulate the zero thickness sheet model of practical metasurface structures. It is shown that if the surface susceptibilities of a given metasurface are expressed as a ratio of two polynomials of tangential spatial frequencies, k||, with complex coefficients, they can be conveniently expressed as spatial derivatives of the difference and average fields around the metasurface in the space domain, leading to extended forms of the standard Generalized Sheet Transition Conditions (GSTCs) accounting for the spatial dispersion. Using two simple examples of a short electric dipole and an all-dielectric cylindrical puck unit cells, which exhibit purely tangential surface susceptibilities and reciprocal/symmetric transmission and reflection characteristics, the proposed concept is numerically confirmed in 2D. A single Lorentzian has been found to describe the spatio-temporal frequency behavior of a short dipole unit cell, while a multi-Lorentzian description is developed to capture the complex multiple angular resonances of the dielectric puck. For both cases, the appropriate spatial boundary conditions are derived.

An accelerated Integral Equations (IE) field solver for determining scattered fields from electrically large electromagnetic metasurfaces utilizing Fast Multipole Method (FMM) is proposed and demonstrated in 2D. In the proposed method, practical general metasurfaces are expressed using an equivalent zero thickness sheet model described using surface susceptibilities, and where the total fields around it satisfy the Generalized Sheet Transition Conditions (GSTCs). While the standard IE-GSTC offers fast field computation compared to other numerical methods, it is still computationally demanding when solving electrically large problems, with a large number of unknowns. Here we accelerate the IE-GSTC method using the FMM technique by dividing the current elements on the metasurface into near- and far-groups, where either the rigorous or approximated Green’s function is used, respectively, to reduce the computation time without losing solution accuracy. Using numerical examples, the speed improvement of the FMM IE-GSTC method O(N1.5) over the standard IE-GSTC O(N3) method is confirmed. Finally, the usefulness of the FMM IE-GSTC is demonstrated by applying it to solve electromagnetic propagation inside an electrically large radio environment with strategically placed metasurfaces to improve signal coverage in blind areas, where a standard IE-GSTC solver would require prohibitively large computational resources and long simulation times.