A multifunctional ultrathin flexible bianisotropic metasurface with miniaturized cell size

In this paper, a flexible bianisotropic metasurface possessing omega-type coupling is presented. The designed metasurface behaves differently when excited from either forward (port 1) or back (port 2) sides. It provides an absorption of 99.46% at 15.1 Gigahertz (GHz), when illuminated from port 1, whereas, on simultaneous illumination from port 2, it behaves like a partially reflective surface (PRS). Furthermore, the presented metasurface not only acts as an in-band absorptive surface (port 1) and partially reflective surface (port 2), but it also provides 97% out-of-band transmission at 7.8 GHz. The response of the presented metasurface remains the same for both transverse Electric (TE) and transverse magnetic (TM) polarized wave or any arbitrary linearly polarized wave. Additionally, the response of the metasurface is angularly stable for any oblique incidence up to 45º. The proposed ultrathin flexible metasurface with absorption, partial reflection and out-of-band transmission properties can be used in the Fabry Perrot cavity antenna for gain enhancement with radar cross-section (RCS) reduction both for passband and stop-band filtering, and conformal antenna applications.

Electromagnetic (EM) metamaterials are artificially engineered EM-structures comprising of sub-wavelength periodic i.e., dielectric or metallic structures that resonantly couple to the EM fields while exhibiting properties that are not readily found in nature 1 . As a successor of metamaterials, two-dimensional metasurfaces 2 , owing to their ease of fabrication, have attracted a lot of attention for producing nearly arbitrary wavefront and phase shift of EM waves. Due to their rare and exceptional abilities to distort phase, polarization and the amplitude of EM waves, they play a key role in antenna performance enhancement, superlenses, scattering reduction,phase manipulation, reflection control, or polarization rotation [3][4][5][6][7][8][9][10] . Another interesting area in this emerging field is the absorption of EM waves via anisotropic structures. Several such designs have been reported in the literature [11][12][13][14] , but they usually perform with a single functionality due to the incorporation of a complete ground plane. The current research impetus is not only to design compact but also multifunctional absorptive structures. The exploitation of the opportunity to design a resonant absorber, which is transparent outside of the absorption band, could open up various novel prospects in applications, for instance, in ultrathin filters for wave trapping, selective multifrequency bolometers, and sensors 15 . To attain multiple functionalities from a single metasurface, the following two techniques can be possibly used: (1) Multi-layered metasurface (2) Bianisotropic metasurface Although, a multilayer metasurface has the advantage of wider bandwidth, it exhibits some disadvantages including fabrication complexity, high cost, incident angle dependency and bulky size 16 . The only suitable, cost-effective and unexplored technique to achieve multiple functionalities in absorptive metasurfaces is bianisotropy [18][19][20][21][22][23][24][25][26][27] .
Metasurfaces that exhibit electromagnetic coupling are known as bianisotropic metasurfaces. Bianisotropic metasurfaces can be classified into four categories depending on their functionalities and internal properties. In 28 , the classification of bianisotropic materials has been discussed based on the magneto-electric and electromagnetic dyadic nature for both reciprocal media (i.e., chiral and omega) and non-reciprocal media (i.e., Tellegen and moving particle). By enabling Chiral coupling and Tellegen coupling in a metasurface, one may control the reciprocal and non-reciprocal polarization transformation. Similarly, asymmetric co-polarized reflection and transmission can be achieved through Omega and artificial moving particle coupling. Taking this into consideration, Ra' di et al. 29 discussed the theoretical modeling of omega bianisotropic metasurfaces and the design possibilities through reciprocal omega coupling. One of the key functionality that omega coupling enables in a structure is the asymmetric absorption. To attain perfect absorption when illuminated from one side and the multiple functionalities when the same metasurface is illuminated from the other side, both uniaxial reciprocal (omega) and non-reciprocal (Tellegen and AM) bianisotropic properties, are required. The research community mainly focus on reciprocal metasurfaces due to their ease of fabrication, low cost, and lack of active elements (i.e. transistors and diodes) in the structure.
Based on the theoretical analysis given by Ra' di et al 29 , the Yzadi et al 30 proposed a simple reciprocal bianisotropic metasurface consisting of metallic square patches that possess omega coupling for the very first time. The designed structure provides perfect absorption when illuminated from one side and gives asymmetric reflection when illuminated from the other side. However, its claimed performance deviates when exposed to different polarized waves (polarization between TE and TM). In addition, it provides very low out-of-band transmission of 20%. Recently in order to overcome this limitation, Fahad et al 31,32 . proposed a metasurface exhibiting polarization-insensitive behavior. Although the polarization insensitivity problem was addressed, it came with a toll of having a larger unit-cell size and low out-of-band transmission value of 10% 31,32 . The low out-of-band transmission and larger unit-cell size limit the usage of these designs for several applications.
In this paper, a miniaturized ultra-thin flexible bianisotropic metasurface is designed and discussed. The proposed metasurface is not only flexible but it also provides more than 97% out-of-band transmission power, thus, overcoming the transmission-based limitations. The maximum out-of-band transmission is achieved at 7.8 GHz, thus, enabling the band-pass filtering functionality, which is a salient feature of the proposed metasurface. Maximum absorption of 99.4% at 15.1 GHz when illuninated from port 1 and partial reflections of nearly 99% when illuminated from port 2 are achieved through the same structure. The maximum out-of-band transmission with bandpass filtering functionality, the miniaturized size, and the flexible substrate make the proposed design a suitable prospect for several applications.

Theoretical analysis and design evolution
Theoretical background. It is evident from 28 literature that to achieve perfect absorption, the thickness of a metasurface cannot be zero, as a result of the existence of at least two induced moments e.g., the electric and the magnetic moment. These induced dipole moments are linearly dependent on the externally applied incident field and can be related using polarizability dyadic as follows: here E inc is an incident electric field and H inc is an incident magnetic field, which are related as follows: In Eq. 1, α ee ,α mm ,α em and α me are the effective electric, magnetic, electro-magneto, and magneto-electric polarizability components of the proposed metasurface. The effective polarizabilities include the effects of a single particle in an array 29 . The terms α ee ,α mm show the electric response to an incident electric field and the magnetic response to an incident magnetic field, respectively. The other two terms α me and α em represent bianisotropic coupling dyadics. Due to uniaxial symmetry, only isotropic response and rotation around z o axis is allowed, hence polarizability dyadics can be written as: (3) α ee =α co ee I t +α cr ee J t (4) α em =α co em I t +α cr em J t (5) α me =α co me I t +α cr me J t (6) α mm =α co mm I t +α cr mm J t (7) It is evident from the above relationship that by tuning one side of the structure as a perfect absorber, one may achieve special properties like asymmetric reflection and transmission from the other side. Reciprocal structures that enable such behavior are known as omega structures. They possess omega-type coupling that is responsible for enabling asymmetric reflections and perfect absorption simultaneously in the same metasurface. For the reciprocal structure, electric and magnetic polarizability dyadic are symmetric i.e., α cr ee , α cr mm = 0 , while the bianisotropic coupling dyadic meets the following conditions: where α cr em represents omega coupling dyadic and α co em represents chiral coupling dyadic. For an ultimately thin symmetric absorber, the electric and magnetic polarizabilities should be in Huygens balance 30 : Furthermore, asymmetric absorption/reflection from any structure can be realized when the following conditions are satisfied 29 In omega structure chirality α co em parameter is zero. By applying this condition in Eq. (12), the reflection coefficient can be written as shown in Eq. (19). By analyzing this relationship, one can easily perceive that by varying α cr em (omega coupling dyadic), the co-polarized reflection from the opposite side can be controlled, while maintaining total absorption from the original front side.
The co-and cross-polarizability dyadic (electric, magnetic, electro-magneto, and magneto-electro) of the proposed unit cell are calculated and shown in section III of this paper. To enable omega coupling combining two or more conductive inclusions in a one-unit cell is suggested in 30 . Using such approach to design omega bianisotropic metasurfaces not only reduces the fabrication cost and complexities but also makes them compatible with planar technologies. The presented structure is periodic in the x and y-axis, hence acts as an infinite planar array of electrically small resonant inclusion. These proposed bianisotropic unit cells are arranged in a square lattice with a period of "a" and excited by an arbitrary linearly polarized plane wave. It is worth mentioning here that the unit cell array period is kept smaller than the operating wavelength to avoid grating lobes in scattering spectra. The following key points are considered while designing the unit cell: • The proposed structure should be geometrically asymmetric, so as to behave asymmetrically on excitation from the + z and -z-directions. Such asymmetry provides the required bianisotropic response in the structure. • The proposed structure is designed using two conductive inclusions so that it enables omega coupling, and thus acting as an absorber and partially a reflective surface when excited from the + z and -z-axis respectively. • Instead of the ground plane (conventionally used in absorbers), a unique circle with outer radius "R" and inner radius "r" is imprinted on the backside, consequently, increasing the out of band transmission. Furthermore, a hole is created to enhance this out-of-band transmission.

Design evolution.
A bianisotropic metasurface is the one in which D(r,t) is anisotropically coupled to both E(r,t) and B(r,t). Similarly, H(r,t) is anisotropically coupled to both E(r,t) and B(r,t). In order to induce the omega coupling, the proposed structure is designed using asymmetric inclusions 30 , due to which it behaves asymmetrically on excitation from the + z and − z directions. The main purpose of this paper is to design a metasurface with asymmetric reflections providing maximum absorption when excited from port 1 and partial reflections when excited from port 2. The proposed design is optimized through systematic parametric analysis. The infinite unit cells along x-and y-axis are simulated in Computer Simulation Technology (CST) Studio Suite. The evolution steps of the proposed metasurface have been shown in Fig. 2a-d. These are few, but critical design steps performed to realize the final and optimized design. The resonance frequency of the proposed structure is calculated 37 using Eqs. (20) and Eq. (21)  www.nature.com/scientificreports/ In step 1, the unit cell that consists of a square patch enclosing a circle with four rectangular stubs backed with the complete ground is designed as shown in ig. 2(a). In step 2, the ground plane is replaced by an optimized square patch to enable out-of-band transmission. As it can be seen from Fig. 2b the maximum out-of-band transmission is achieved at 13.1 GHz, while the absorption at 15.1 GHz is reduced to 50%. In order to maximize both the transmission and absorption from port 1 with partial reflection from port 2, a square patch is replaced with an optimized circular patch in step 3 as shown in Fig. 3c. As a result, maximum absorption of 99.46% from port 1 and partial reflections of 97% from port 2 is achieved at 15.1 GHz, where the out-of-band transmission reduces to 65%. In the final step, a hole of 'r = 0.75 mm' is created in a unit cell design to enhance the out of band transmission 39 as shown in Fig. 2d. Hence, an absorption of 99.46% from port 1 and simultaneous partial reflections of 99% from port 2 at 15.1 GHz with maximum out of band transmission of 97% at 7.8 GHz is achieved. Unit cell equivalent model. To find the resonant frequency of the proposed structure, its LC circuit model is designed. To analyze the unit cell analogous to a transmission line, port 1 and port 2 are placed on the LC circuit model as shown in Fig. 3a. ne end of each of the two ports is connected to the upper side of the unit cell and the other end is connected to its lower side. Hence, the input signal entering port 1 and leaving port 2 experiences the same impedance as suffered by the overall metasurface. Starting from port1 to port 2, the conductive lengths, where no conductor is present on the opposite side, are modeled as an inductor, whereas the conductive lengths which have an equal metallic area on the opposite side are modeled as a capacitor. With this approach, the whole unit cell is modeled in terms of lumped elements, and its equivalent circuit is shown in Fig. 3b. The values of capacitors and inductors are calculated using parallel plate capacitance and flat wire inductance equations  Fig. 4. Figure 4a, b show the co component of the electric α co ee and the magnetic α co mm polarizabilities respectively, and Fig. 4c shows the Electromagneto and the magneto-electro response of the proposed unit cell. At the resonance frequency at 15.1 GHz, the real part of the electric polarizability does not cross zero, which is due to the non-zero static electric polarizability. In general, it can be concluded that the response is dominated by the induced electric dipole moment and therefore the mode is called electric mode 30 . It can be seen from Fig. 4c that both magnitude and phase of the cross-component of the electro-magneto ( α me ) and magneto-electro ( α em ) dyadics are equal hence the condition for omega coupling i.e. α cr em = α cr me is satisfied. It is due to the omega coupling, present at the resonance frequency, that the proposed metasurface acts as an absorber when illuminated from one side and, and like a partially reflecting surface PRS 31,32 when illuminated from the other side.
Reflection/transmission coefficient. For x polarized incident wave, the co and cross components of the reflected and transmitted wave are:  Fig. 5a,b for port 1. Similarly, for port 2, the co-and cross-components of the reflected and transmitted wave, for both x and y polarized waves are shown in Fig. 5c,d. At the resonance frequency − 15.1 GHz-the values of the co-polarized reflection coefficient (port 1) for x and y polarized incident waves are approaching -25 dB and -22 dB respectively. The existence of such strong resonances is responsible for the nearly perfect absorption from port 1. However, for port 2, the co-polarized reflection coefficient ( R yy and R xx ) remains high i.e., 99%, whereas all other reflected and transmitted components are below -10 dB. Therefore, it can be said that the proposed unit cell provides the asymmetric reflection at the resonance frequency when illuminated by EM waves from opposite sides.
Absorption. The absorption of the proposed metasurface for port 1 can be found using the following mathematical relationship 37,38 : where Here Z o is the characteristic impedance of the free space that is equal to 377 Ω and Z in is the effective impedance. To calculate the value of the effective impedance ( Z in ), constitutive parameters are extracted using susceptibilities 31,32 . Figure 6, it can be easily observed that the value of Z in is 0.955 i.e., 360 Ω. By putting these values in Eq. (15), R yy 2 comes turns out to be 0.0005, and the extracted value of T yy 2 is equal to 0.0018. By inserting the calculated values of R yy 2 and T yy 2 in Eq. (23) the absorption at 15.1 GHz becomes approximately equal to 99.7%. The absorption of port 1 for linearly polarized incident wave is calculated from the simulated results and it is shown in Fig. 7. It is evident that absorption of 99.46% is achieved at 15.1 GHz for both x and y polarized incident waves. 14 14 This exclusive property makes this metasurface a suitable candidate for those applications, where the maximum in-band absorption with non-zero out-of-band transmission is required.
Angular stability. The robust response against oblique incidences plays a crucial role in many practical applications, thus, angular stability analysis is performed. An angularly stable metasurface, is a structure for which the response of the metasurface in normal incidence remains invariant when the same EM wave is incident at oblique angles. The major factors that contribute to angular stability are given below: 1. Thickness The electrical thickness of a substrate should be chosen in such a way that it should be smaller than the wavelength at the resonant frequency so that the spatial uniform fields can be induced in the substrate, which is immutable to variation of the angle of incidence 31,32 . Hence, Rogers4350B Lo pro substrate with a thickness of 0.1016 mm (electrical thickness of 0.005 λ at the resonant frequency) is chosen for the metasurface implementation. 2. Unit cell size The overall unit cell size (p=7.5 mm, 0.37 λ) is also kept smaller than the wavelength at the frequency of interest so that the uniform spatial fields may appear across the unit cell, which again makes this structure immutable to variations of the angle of incidence. 3. Unique configuration A unique geometric configuration with optimized dimensions is used so that it becomes angularly stable. To check the angular robustness of the proposed unit cell inclusion, the angle of incidence of EM wave is varied from 0˚ to 45˚.
It is evident from Fig. 9a,b that the response of the proposed unit cell is stable up to 45˚ for both x and y polarized incident waves. For the validation of this property for the proposed metasurface, the polarization of the incident wave is rotated from 0˚ (TE) to 90˚(TM). It can be easily noted from Fig. 10b that the absorptive response of the proposed metasurface remains immutable due to the cyclic-4 (C4) symmetry i.e., 90˚ rotated symmetry and mirror symmetry presented along x, y,u, and v axes as shown in Fig. 10a.

Experimental results and discussion
To validate the performance of the proposed omega-type bianisotropic metasurface, a prototype was fabricated on Rogers's 4350B Lo pro substrate having 0.1016 mm thickness. The fabricated PCB consists of 41 × 30-unit cells and its overall size is 304.8 × 228.6 mm 2 as shown in Fig. 11a,b. The holes on the PCB are created using non-PTH drilling. An experimental setup is created inside an anechoic chamber where two horn antennas (operating frequency range of 2-18 GHz) are utilized for illuminating the surface and receiving the reflected/transmitted waves, as shown in Fig. 12a,b. These antennas were placed on the same side at a 120 cm (Fresnel region) 41,42 from the surface for the reflection measurements. As far as the transmission measurements are concerned, the antennas are separated by a 240 cm distance and the metasurface is kept in the middle as shown in Fig. 12a 30 . As a reference for reflection and transmission measurements, the metallic plate and air measurements have www.nature.com/scientificreports/ been carried out before performing the actual prototype's measurements. The horn antennas are attached to the vector network analyzer (Anritsu-MS46122B) through coaxial cables. To measure co-polarized transmission ( T yy or T xx ), both transmitting and receiving antennas were placed either vertically (for T yy ) or horizontally (for T xx ). Similarly, for the measurement of the co-polarized reflection, the two antennas were placed parallel to each other as shown in Fig. 12b, for T yy both the transmitting and receiving antennas were vertical, while for T xx both antennas were placed horizontally. It is worth mentioning that only co-components of transmission/reflections are considered in the measurements because the cross-components remained below -50 dB throughout the frequency band. The experimental results are summarized in Fig. 12c-e, and they are in reasonable agreement with the simulated results. It can be verified in Fig. 12c,d, that the phenomena of asymmetric reflection and symmetric out-of-band transmission are validated from the experimental measurements. The small discrepancies can be attributed to multiple reflections from the surrounding items, manual handling errors and the substrate tolerances. Moreover, the slight deviation is due to the finite size of the fabricated metasurface as opposed to the infinite size assumption that is made for the simulation model.
A comparison of the implemented omega bianisotropic metasurface with MSFs reported in literature is summarized in Table 1. The proposed flexible metasurface slightly improved absorption when compared with previously reported omega metasurfaces and it outperforms the competition in terms of out-of-band transmission (97%). At the same time, its normalized unit cell size is the smallest (0.37 λ) and the used material is flexible in nature contributing to the multiple functionalities of the proposed metasurface.

Conclusion
In this paper, a flexible omega bianisotropic metasurface has been presented. It is shown that the designed metasurface acts as an absorber (giving maximum absorption of 99.46%) when illuminated from port 1, whereas, on simultaneous illumination deriving from port 2 it behaves like a partially reflective surface. The proposed metasurface causes very high out-of-band transmission of up to 97%. The response of the presented metasurface  www.nature.com/scientificreports/ remains the same for both TE and TM polarized waves or any arbitrary linearly polarized wave. Moreover, the response of the metasurface is stable for oblique incidences up to 45º. Such metasurfaces, with maximum absorption properties, partial reflections, and symmetric out-of-band transmission can be used in cavity antennas for gain enhancement with RCS reduction, for passband and stopband filtering applications, and due to its flexibility, it can be used in conformal antenna applications.

Data availability
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.